Abstract

When it comes to the identification of dynamic system parameters, like stiffness and damping, the systematic measurement uncertainty is mostly ignored in the frequency domain because of its complicated and elaborate uncertainty propagation. Nevertheless, the uncertainty is far from being negligible. In order to examine the importance of the systematic measurement uncertainty in the frequency domain, the propagation of systematic uncertainty caused by four different sensor error types is investigated by using Monte Carlo simulations. It is shown that an uncertainty in the bias only affects static measurements, i.e., with an excitation frequency of the system being Ω=0. Furthermore, the uncertainty of the linearity and sensitivity only impacts the magnitude of the Fourier coefficient, whereas the uncertainty in the hysteresis results only in an uncertainty in the phase. Subsequently, it is shown that the linearity and hysteresis uncertainty are completely independent in the frequency domain. The results lead to simple propagation relations applicable to various engineering, science, and industry challenges without the need for the conduction of complex Monte Carlo simulations.

1 Introduction

The determination of the dynamic characteristics of systems is of great importance in science, engineering, and industry. If the dynamic behavior of a system is known, the optimal design and control can be determined, for instance, through analytical models and their dependence on the design parameters, cf. Pelz et al. [1]. For complex mechanical systems such as hydraulic accumulators or modern turbomachinery, the analytical determination of the transfer behavior may be feasible; however, the identification of the system via experiments is usually preferred. The identification of the time-invariant system, e.g., the identification of the system inherent descriptive linear components, i.e., springs, dampers, and masses, is most commonly performed in the frequency domain [2,3].

In practice, this often but not exclusively involves measuring the response of the dynamic system to harmonic excitations, e.g., the frequency-dependent force of an air spring due to an imposed deflection, cf. Puff and Pelz [4] and Hedrich et al. [5,6]. The measurement signals are subsequently transformed into the frequency domain using fast Fourier transforms (FFT), allowing the frequency-dependent stiffness of the air spring to be determined.

While determining dynamic characteristics in the frequency domain is common practice, a proper quantification of the uncertainty of the identified quantities is often nonexistent. This is mainly due to the nontrivial propagation of the inherent random and systematic uncertainty of the measured signals from the time into the frequency domain or vice versa. Interestingly, however, for one of the two uncertainty types, i.e., the propagation of the stochastic uncertainty, often referred to as type A uncertainty according to the Guide to the expression of uncertainty in measurement (GUM) [7], into the frequency domain or within it, numerous publications are available, cf. Betta et al. [8]; Medina et al. [9]; Medina et al. [10]; Eichstädt and Wilkens [11]; and Moreland et al. [12]. In contrast, the propagation of the systematic measurement uncertainty with an unknown distribution, i.e., the type B uncertainty, has only been addressed by means of elaborate Monte Carlo simulations or propagation laws based on conservative assumptions, cf. Eichstädt et al. [13] and Kuhr [14].

Betta et al. [8] investigate the propagation of stochastic uncertainty into the frequency domain due to quantization errors, sampling frequency, jitter, and finite word length of microprocessors. The work mainly deals with the uncertainty due to the discrete Fourier transform (DFT) itself rather than the propagation of uncertain measurement signals.

Medina et al. [9] gave a general method for uncertainty propagation in parameter identification of frequency response functions. However, the authors did not address type B uncertainty when setting up their covariance matrix. Furthermore, the presented method is only suitable for the identification of linear systems. However, to analyze the uncertainty of the broad range of systems, the propagation of systematic measurement uncertainty into the frequency domain is required.

Medina et al. [10] presented a methodology to propagate experimental uncertainty applied to the identification of a simple two degree-of-freedom mechanical system. Although the suggested method has been developed independently of the types of method to evaluate uncertainty, only type A uncertainty, i.e., stochastic uncertainty, has been investigated.

Eichstädt et al. [13] use Monte Carlo simulations, according to the GUM and its supplements, to calculate dynamic measurement uncertainty. The benefit of this approach is that all types of uncertainties, i.e., type A and B, can be mapped in general. The authors identify a limitation of Monte Carlo simulations in terms of the high demands they place on available memory. In response to this challenge, they present two memory-efficient alternatives. However, the authors still use a workstation computer with 24 cores and 96 GB of RAM, which are still atypical for normal personal computers.

Eichstädt and Wilkens [11] developed a software package to implement closed formulae for the propagation of uncertainty between the time and frequency domain, according to the Joint Committee for Guides in Metrology [7]. The software provides reliable uncertainty propagation for common tasks in dynamic metrology, such as propagation between the frequency and the time domain, between amplitude/phase and real/imaginary part representation, and for estimating the dynamic measured value through deconvolution in the frequency domain. However, only uncertainty of type A is considered.

Moreland et al. [12] measured the static rotordynamic characteristics of liquid annular seals. The uncertainty quantification of the frequency-dependent stiffness and damping values is solely based on the standard deviation of the mean values of the computed complex stiffnesses within the frequency domain. The classical Gaussian uncertainty propagation is used to calculate the measurement uncertainty. It is stated that errors due to instrumentation are largely negligible, and only repetition determines the uncertainty.

Kuhr [14] takes a different approach. The author assumes that the overall measurement uncertainty calculated in the time domain can generally be transferred to the complex amplitude of the Fourier coefficients in the frequency domain. Assuming equal uncertainty of real and imaginary part, a simple relation between time and frequency domain is derived. However, the author explicitly points out that the method represents a conservative estimate of the measurement uncertainty in the frequency domain and is only approximately valid.

All of the above references focus on time-variant systems, i.e., systems without dynamic measurement errors due to the transfer functions of the equipment used. Within the presented paper, we solely focus on time-invariant systems. For further information on the propagation of dynamic measurement errors and their impact on measurement uncertainty, the interested reader is referred to the work of Vasilevskyi [15], Hessling [16,17], and Elster and Link [18,19].

It is evident that most references neglect the propagation of systematic measurement uncertainty with unknown distribution, the underlying reason, except for Monte Carlo simulations, remains unclear. One hypothesis is that the propagation of systematic measurement uncertainties is mainly neglected due to the time-consuming nature of suitable propagation methods, or the extensive memory effort and knowledge requirement of the correlation between the input quantities for covariance analysis. The extensive memory effort is particularly important when dealing with long measurement signals, i.e., where the size of the covariance matrix increases over-proportionally with the samples of a signal. To emphasize this claim, Table 1 gives the size and computation time of the covariance matrix, as well as the overall calculation time.

Table 1

Comparison of the computation time for different signal lengths using METAS unclib [20]

SamplesSignal generation time in sCovariance matrix generation time in sTime for FFT in sTotal time in sRAM space in GB
1036.60.061.17.80.0149
1048623013891.4901
1051188NaNNaNNaN149.01
SamplesSignal generation time in sCovariance matrix generation time in sTime for FFT in sTotal time in sRAM space in GB
1036.60.061.17.80.0149
1048623013891.4901
1051188NaNNaNNaN149.01

The propagation from the time into the frequency domain, i.e., the FFT, is carried out by using the advanced measurement uncertainty calculator METAS unclib [20] and the commercial calculation software matlab with 32 digits precision. It should be noted that for a typical long measurement signal, being recorded at a fairly low acquisition rate for dynamic measurements of 6000 samples per second for 5 s, cf. Sec. 4, i.e., the signal length being 105 samples, the covariance matrix could not be computed using a state-of-the-art workstation, indicated by NaN, due to its requirements of RAM space.

In light of the aforementioned considerations, two fundamental challenges emerge. First, systematic uncertainty is frequently overlooked or subjected to strong assumptions when propagated into the frequency domain, resulting in a significant under- or overestimation of the combined uncertainty of the identified system parameters. Second, time and memory limitations pose significant technical challenges to the adequate quantification of uncertainty. To address those challenges, the paper addresses the following two questions:

  • How does the systematic uncertainty affect the magnitude and phase in the frequency domain?

  • Are there simple and fast propagation laws for the systematic measurement uncertainty into the frequency domain?

To answer the questions, the main differences between error and uncertainty as well as the difference between the systematic and stochastic uncertainty are briefly recapped. Subsequently, an analytical sensor model for modeling (i) offset, (ii) sensitivity, (iii) linearity, and (iv) hysteresis errors is presented. Combined with the corresponding probability distributions, this allows the description of the systematic uncertainty. The propagation of the four different systematic sensor uncertainties into the frequency domain is investigated by using Monte Carlo simulations. The simulation results are then analyzed to derive simple propagation relations from the time into the frequency domain. It will be shown that the systematic measurement uncertainty is far from being negligible for determining the overall uncertainty of the identified system parameters. Subsequently, the derived simple propagation laws are applied to a complex parameter identification procedure, comparing the derived propagation laws to a previously presented approach.

2 Modeling Systematic Uncertainty

2.1 Errors and Uncertainty.

Prior to analyzing the propagation into the frequency domain, it is first necessary to define and distinguish between the concepts of error and uncertainty.

All measured data are uncertain due to unknown errors. An error Ei is defined as the deviation between the observed, i.e., measured, quantity yi and the real physical value xi (Joint Committee for Guides in Metrology [7])
(1)

While an error is a specific deviation, unknown measurement errors are treated statistically using probability distributions. This statistical approach allows the quantification of an uncertainty, providing a framework for analysis through statistical methods. In other words, the uncertainty δi of the measured value yi describes an interval or probability distribution within which the true value is likely to lie in.

In general, there are two distinct types of errors when it comes to the measurement of physical quantities: systematic, and therefore repeatable errors, and random, nonrepeatable errors. The former of the two types, i.e., systematic errors, can be further divided into static and dynamic errors, i.e., errors which are important for time-invariant and time-variant systems. As stated before, time-variant systems are not in the scope of this work.

Regarding the remaining category of nonrepeatable errors, namely, random errors, they can be further subdivided into numerous distinct types, including repeatability, measurement noise, time variance, and quantization errors, cf. Elhoushi et al. [21]. An overview of the different types of errors is provided in Fig. 1. The errors that are the primary focus of this paper are highlighted in gray, i.e., (i) bias errors, (ii) sensitivity errors, (iii) linearity errors, and (iv) hysteresis errors.

Fig. 1
Distinctions of possible sensor errors. The focus of this paper, i.e., systematic sensor errors according to Ref. [22], is highlighted in gray.
Fig. 1
Distinctions of possible sensor errors. The focus of this paper, i.e., systematic sensor errors according to Ref. [22], is highlighted in gray.
Close modal

For the possibility to take those errors in an experimental setup into account, i.e., addressing the distinct measurement errors as an uncertainty, the use of probability distributions necessitates. For the type A uncertainty defined by GUM [7], the distributions are known and can be treated using statistical methods. In contrast, the probability distributions of systematic measurement errors are unknown. The systematic measurement uncertainty is then defined as type B uncertainty, i.e., an uncertainty which encompasses other uncertainties.

In accordance with the static systematic errors given in Fig. 1, cf. Tränkler and Fischerauer [22], the four main systematic measurement uncertainties regarding the characteristic of a sensor are:

  • zero offset, bias, δbias,

  • sensitivity δsen,

  • linearity δlin=fn(x), and

  • hysteresis δhys=fn(x,x˙).

In this context, bias refers to a shift in the y-axis section of the sensor’s characteristic curve. The uncertainty of sensitivity pertains to the deviation of the sensor’s stated sensitivity from its actual value as specified in the data sheet. Similarly, the uncertainty of linearity represents the deviation between the real value and the measured value due to a nonlinear sensor characteristic. Finally, the uncertainty of hysteresis indicates the maximum deviation of a measured value depending on the direction from which it is approached, cf. Fig. 2.

Fig. 2
Sensor model illustrating the four types of systematic errors considered
Fig. 2
Sensor model illustrating the four types of systematic errors considered
Close modal

All of these systematic uncertainties can be reduced by calibration with much more accurate measurement systems. In general, the manufacturer’s specifications are used as the systematic uncertainty in measurements. These are usually given in a sensor’s datasheet. However, for simplicity, some manufacturers give only a general uncertainty band that covers all possible sensor errors. Table 2 gives an overview of the relative and the absolute systematic uncertainties of an exemplary Hottinger Brüel & Kjær force transducer HBK U2B\1KN with a measurement range of 1kN [23].

Table 2

Uncertainty of the force transducer U2B\1KN from HBK at 33°C ambient temperature [23] (FS = fullscale; MV = measurement value)

Uncertainty typeRelative uncertaintyAbsolute uncertainty
Offset due to temperatureδbias0.1%FS1N
Sensitivity due to temperatureδsen0.05%MV0.5Nat33°Catx=1kN
Linearityδlin0.1%FS1N
Hysteresisδhys0.75%FS0.075N
Uncertainty typeRelative uncertaintyAbsolute uncertainty
Offset due to temperatureδbias0.1%FS1N
Sensitivity due to temperatureδsen0.05%MV0.5Nat33°Catx=1kN
Linearityδlin0.1%FS1N
Hysteresisδhys0.75%FS0.075N

2.2 Systematic Measurement Errors.

Prior to analyzing the propagation of the systematic measurement uncertainty into the frequency domain through Monte Carlo simulations, it is first necessary to derive a suitable error model. The distribution of this error model represents the uncertainty of the sensor. In general, a sensor has a characteristic which outputs a measured value y=fn(x) depending on its input x. This characteristic is often assumed to be strictly linear, i.e., the transfer function simplifies to
(2)
with the sensitivity of the sensor s and the offset o. Regardless, this linear characteristic is, naturally, subject to a certain degree of uncertainty. To model the real, i.e., error-prone sensor characteristic, it can be assumed that the systematic error can be specified as an additional term in the transfer function of the sensor. This additional term depends on the input value x and its first time derivative x˙
(3)
The additional term encompasses the four distinct categories of the static, systematic measurement errors Ei, i.e., the error of the bias, sensitivity, linearity, and hysteresis, cf. Tränkler and Fischerauer [22]
(4)
The modeling of the first error, i.e., the bias, as well as the second one, i.e., the sensitivity error, is straightforward. Here, the bias is modeled by a simple constant value, Cbias, whereas the error of the sensors’ sensitivity is given by a constant, Csens times the input value x
(5)
The values of the constants without the corresponding probability function are typically given in the sensor’s data sheet, cf. Table 2. To be able to propagate the uncertainty and not only perform an interval estimation, the GUM [7] proposes the assumption of a uniform distribution U(±Ci) as a conservative estimation. The uncertainty of the bias and the sensitivity is given by
(6)
and
(7)

respectively. To determine the 95% confidence interval from this assumption δ95, the specified interval of the uniform distribution is reduced by 5%, i.e., δ95=0.95/2[max(C)min(C)]. In the following, all uncertainties are given within the 95% confidence interval. Therefore, the index 95 is neglected.

In contrast to the first two errors, the modeling of the remaining ones, i.e., the systematic error of the linearity as well as the systematic error of the hysteresis, is far from being straightforward. Starting with the linearity error, it is obvious that it is solely a function of the input value Elin=fn(x) rather than the first time derivative x˙. However, its particular relationship is unknown, and the literature does not provide sufficient methods or hints for its modeling purpose. Generally, there are an infinite number of possible functions for description purposes. In this paper, we assume that the function fn(x) can be decomposed into a sum of sine functions, similar to the Fourier decomposition of periodic signals. This assumption is reasonable because harmonic excitations are typically used as inputs for identifying systems’ frequency response functions. Although an alternative approach using a sum of arbitrary powers, similar to a Taylor series, is possible, it is not pursued further. Instead, the linearity error is modeled by
(8)
Here, the parameters ai,wi, and φi, respectively, the amplitude, wavelength, and phase shift, are to be considered as statistically distributed independent quantities, whereas M is the number of spectral components. The sensor characteristic combined with this model generates an arbitrary characteristic curve of the sensor, which always produces the same systematic error when repeating the measurements. To ensure that the error remains within the specified limits, either the distribution of the amplitudes can be varied or the error itself can be scaled. Here, the error is scaled to the limits of the specified interval to obtain a conservative estimate within the limits δlin
(9)
Finally, the error of hysteresis is modeled. Similar to the linearity error, the error in hysteresis is represented by an unknown function. However, it additionally depends on the direction of measurement and the measure x, specifically the first time derivative x˙. The directional dependence is described by the signum function of the time derivative of the measure, sign(x˙). To achieve continuous signals, the signum function is approximated using the hyperbolic tangent, tanh(x˙). Consequently, the hysteresis error is modeled and scaled similarly to the linearity error
(10)

Figure 2 gives an overview of the modeled sensor errors. The thin lines represent the ideal sensor, whereas the thick lines give the influence of the error. The dashed lines show the uncertainty band given by the data sheet.

2.3 Parametrization of the Linearity and Hysteresis Error Model.

The initial question that arises is how the input parameters of the modeled linearity and hysteresis errors should be distributed to achieve a uniform distribution to fulfill the assumption of the Joint Committee for Guides in Metrology [7] and therefore derive the uncertainty.

To conduct this analysis, an input signal x is used where each point on the sensor characteristic is sampled at the same frequency, i.e., number of occurrences. Sampling a harmonic signal at equidistant intervals results in an arcsine distribution, which is unsuitable here, as it disproportionately sampled the edges of the sensor characteristic. Instead, a triangular signal is selected as a uniformly distributed input.

The triangular signal f(t) is defined as follows:
(11)

A full period T0 is used so that the characteristic is traversed in both ascending and descending directions. To simplify the analysis, only linearity errors are initially considered.

Subsequently, a Monte Carlo simulation with N=105 deviations from the linear characteristic, i.e., Eq. (2), is conducted, and the frequency distribution of the error is analyzed. The distributions of the parameters are chosen to obtain a uniform distribution of the linearity error.

Table 3 shows the initial parameters of the error model with the selected distributions.

Table 3

Initial distributions used for the error model

ParameterDistribution
aiUniform:
U[0, 1]
φiUniform:
U[0,2π]
wiHalf-normal:
|N(0,1.7)|
ParameterDistribution
aiUniform:
U[0, 1]
φiUniform:
U[0,2π]
wiHalf-normal:
|N(0,1.7)|

In this model, the amplitude ai and phase φi span the full range of possible values, while the frequency wi theoretically has no upper limit. To reflect the decreasing probability of higher frequencies, a half-normal distribution is used, ensuring only positive frequencies are included. The effects of distribution parameters and the number M of spectral components wi on the error are subsequently evaluated.

The uncertainty limits for linearity and hysteresis are set to δlin=δhys=0.01, and the errors are added to the ideal sensor signal using Eq. (4) with sensor sensitivity s=1.

Figure 3 presents three exemplary resulting error distributions. As the variance of the half-normal distribution for wi increases, the resulting distribution transitions from an arcsine to a normal distribution.

Fig. 3
Frequency distributions from the Monte Carlo simulation of the error model, based on the initial parameters in Table 4
Fig. 3
Frequency distributions from the Monte Carlo simulation of the error model, based on the initial parameters in Table 4
Close modal

Varying the distributions of wi and ai alone does not yield a uniform error distribution. To address this, an approach based on the Fourier decomposition of a function that is uniformly distributed under equidistant sampling is explored. This approach couples wavelength and amplitude as ai=fn(wi), reducing the model’s variables to two: wavelength wi and phase φi.

Using the periodic triangular function f(t) from Eq. (11) meets this uniform distribution requirement. Its Fourier transform yields a nonlinear relationship between amplitude ai and wavelength wi
(12)

When this relationship is used, distributions of the error are obtained, cf. Fig. 4. This ensures that the error is distributed uniformly.

Fig. 4
Frequency distributions from the Monte Carlo simulation of the error model, based on the correlation between amplitude and wavelength ai=1/wi
Fig. 4
Frequency distributions from the Monte Carlo simulation of the error model, based on the correlation between amplitude and wavelength ai=1/wi
Close modal

For further investigation, the number of spectral components M and the distribution of wavelengths wi have a significant impact on the output distribution. To choose the appropriate parameters, a full factorial simulation design with M:115 and wi:|N(0,1500)| is conducted. The maximum of the error distribution with 50 bins is used as the evaluation metric. For a uniform distribution, the maximum is 2%. For all other distributions, the maximum is higher. The closer the determined maximum is to 2%, the closer the error is to a uniform distribution.

Figure 5 shows the results depending on the distribution of wavelengths for a constant number of spectral components.

Fig. 5
Dependence of the error distribution on the distribution of the wavelength wi and the number of spectral components M
Fig. 5
Dependence of the error distribution on the distribution of the wavelength wi and the number of spectral components M
Close modal

Generally, higher values of σ(wi) reduce the maximum frequency error, approaching 2% for large σ(wi) values over 500, indicating a uniform distribution.

Furthermore, two curves are noteworthy. If only one spectral component M=1 is used, it is not possible to achieve a uniform distribution. Moreover, with spectral components of M=2, this is not feasible, despite being closer to a uniform distribution than the previously considered curve. Therefore, M3 is chosen.

With an increasing number of spectral components, it is observed that the slope of the curves becomes flatter. This means that for M1>M2, σ1(wi)>σ2(wi) must also apply to achieve the same results. After this analysis, M=3 is chosen. This is also advantageous considering the increasing simulation times and increasing number of spectral components.

Finally, the combined hysteresis and linearity errors are examined. Figure 6 shows the resulting distribution, which is triangular within an error range of −0.02 to 0.02.

Fig. 6
Frequency distribution from the Monte Carlo simulation of the combined linearity and hysteresis error model
Fig. 6
Frequency distribution from the Monte Carlo simulation of the combined linearity and hysteresis error model
Close modal

This result is expected, as the modeled addition of errors results in a convolution of the two uniform error distributions ElinEhys, confirming their statistical independence.

In summary, this analysis determines the input parameter distributions of the modeled linearity and hysteresis errors. For following analysis and the propagation into the frequency domain, an amplitude-to-wavelength dependency ai=1/wi2 with parameters M=3 and σ(wi)=500 is employed.

3 Propagation Into the Frequency Domain

The availability of a description or model for each of the four types of errors, respectively, the corresponding uncertainty, now permits the analysis of the propagation of uncertainties into the frequency domain. Prior to commencing the Monte Carlo simulation, it is prudent to undertake a preliminary examination to ascertain whether its application is necessary for all four types of uncertainty. Starting with the uncertainty of the bias and sensitivity, it becomes evident that both uncertainties can be propagated directly into the frequency domain, i.e., without the need for a proper Monte Carlo simulation
(13a)
(13b)

Here, the ˜ indicated frequency domain variables. It should be noted that the uncertainty of the bias as well as the uncertainty of the sensitivity only affect the magnitude of the Fourier coefficients. This becomes evident when taking the linearity of the Fourier transform into account, cf. Brigham [24]. The underlying probability distributions remain unchanged.

3.1 Monte Carlo Simulation.

It is evident that even when the uncertainties in both bias and sensitivity are determined using the uniform distribution and propagated directly into the frequency domain, cf. Eq. (13), the uncertainties related to the sensor’s linearity and hysteresis, as well as their propagation into the frequency domain, however, are by no means straightforward and sophisticated Monte Carlo simulations are needed.

Monte Carlo simulations have been widely used since the 19th century for various applications, cf. Metropolis et al. [25], Barbu and Zhu [26], etc. According to the GUM [7], Monte Carlo simulations are a suitable method for analyzing the propagation of an uncertainty with a previously unknown distribution. The simulations use numerous calculations to estimate the model output and its distribution with stochastic distributed input parameters [26]. Within this paper, the Monte Carlo simulations are carried out using the commercial software matlab. Each simulation is designed to emulate a real measurement using the aforementioned sensor model. The number of iterations within the Monte Carlo simulation has to be large enough to ensure that the resulting probability distributions no longer change with the number of iterations. This can either be guaranteed by an iterative procedure or by a very high number of iterations. For simple models, the GUM [7] suggests a number of N=106 iterations to get accurate results, which will be put to test at the end of the section.

The simulations are carried out in four steps, shown in Fig. 7.

Fig. 7
Flow chart of the Monte Carlo simulations carried out
Fig. 7
Flow chart of the Monte Carlo simulations carried out
Close modal
First, a generic, monofrequent, harmonic test signal with the offset o and frequency f=Ω/(2π) is generated
(14)

It is assumed that the signal covers the full measurement range of the modeled sensor, i.e., the amplitude â is chosen accordingly. Subsequently, the actual Monte Carlo simulation then starts with N iterations. Second, an error is added to the sensor characteristic curve, creating an erroneous signal. The signal is transferred into the frequency domain in the following step by using matlab’s FFT algorithm. Since the signal length is a multiple of the period, no window function is used. For the statistical analysis, only the coefficient of the frequency of the test signal Y(Ω) is considered. In the final step, the results are evaluated using the standard deviations of the magnitude σ(|Y(Ω)|) and the phase σ([Y(Ω)]).

Table 4 gives an overview of the used parameters of the discrete time vector t and the harmonic excitation.

Table 4

Parameter of the time and signal vector

ParameterValue
Sample frequencyfsample=1kHz
DurationL=1s
Signal frequencyf=10Hz
Offseto=10
Amplitudeâ=1
PhaseΦ=45deg
ParameterValue
Sample frequencyfsample=1kHz
DurationL=1s
Signal frequencyf=10Hz
Offseto=10
Amplitudeâ=1
PhaseΦ=45deg

Finally, the suggested number of Monte Carlo iterations N by the Joint Committee for Guides in Metrology [7] has to be confirmed. Figure 8 shows the standard deviation of the magnitude σ(|Y(Ω)|) of the Fourier coefficient depending on the number of Monte Carlo iterations N. Here, the number of iterations varies between N=103 and N=107.

Fig. 8
Convergence of the magnitude of the Fourier coefficient σ(|Y(Ω)|) over the number of Monte Carlo iterations N
Fig. 8
Convergence of the magnitude of the Fourier coefficient σ(|Y(Ω)|) over the number of Monte Carlo iterations N
Close modal

It exhibits that the evaluation value, i.e., the magnitude of the Fourier coefficient σ(|Y(Ω)|), converges and hardly changes from N106 iterations onward. Therefore, the number of iterations in every following investigation is set to N=106, which is in line with the GUM [7] recommendation. Subsequently, the propagation of the linearity uncertainty and a hysteresis uncertainty is investigated.

3.2 Propagating the Linearity and Hysteresis Uncertainty Into the Frequency Domain.

With the error model and the Monte Carlo simulations available, the propagation of the linearity and hysteresis uncertainty in the frequency domain can be examined, to derive an equally simple propagation law as for the uncertainty of the bias and sensitivity. To do so, all influencing factors of the Monte Carlo simulations must be examined. These factors include:

  • the input parameters and their distributions within the error model,

  • the parameters of the test signal used, and

  • the boundaries of the uncertainty intervals δlin and δhys.

The input parameters of the linearity and hysteresis error model (i) are already set through the investigation of the model itself, cf. Sec. 2.3. The found parameters for the error model meet the assumptions made by the GUM [7].

Initially, the other influencing factors (ii) are varied individually to determine their impact on the results. In the subsequent parameter determination, conservative estimates are aimed for. Finally, the uncertainty intervals (iii) δlin and δhys are examined through a sensitivity analysis. From this analysis, simple and fast propagation rules into the frequency domain can be derived.

Using the obtained input parameters for the distribution of the error model, the influence of the test signal and its parameters are investigated. An influence is expected for the amplitude â of the test signal for two reasons. On the one hand, the amplitude of the input signal corresponds to the measurement range of the sensor. If this changes, it is expected that the wavelengths of the error must also be adjusted. On the other hand, this can also be mathematically justified. To generate the error, the input signal x is used, as shown in Eq. (8). There, the amplitude of the signal is used, according to Eq. (14). To obtain consistent results, the relationship
(15)

must apply. This relationship can also be seen in Fig. 9.

Fig. 9
Influence of the amplitude â of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Fig. 9
Influence of the amplitude â of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Close modal

All other parameters of the signal have no influence on the result in the frequency domain. This indicates that it is an adequate representation of systematic errors. Especially, the parameters of the test signal should have no influence on the results. The variations of the remaining parameters can be found in the  Appendix.

It is noteworthy that the sampling rate fab or the length of the signal L do not affect the results. Both quantities increase the number of measurement points when increased. If these were stochastic errors, the uncertainty would decrease as a result, cf. Betta et al. [8] and Eichstädt and Wilkens [11]. However, systematic errors cannot decrease by increasing the sampling rate or the length of the signal, as they are inherent to the measurement system. From this, it can be concluded that the chosen form of modeling concerns solely systematic errors.

After the influences (ii) have been examined, the distributions of the results in frequency domain are shown in Fig. 10. The distributions are oriented toward the magnitude and phase. A cross is formed in the complex number plane. This suggests that magnitude and phase can be considered independently. If they were dependent on each other, an elliptical shape would be expected. The distributions obtained with the chosen parameters do not correspond to a uniform distribution, although this is present in the time domain, cf. Sec. 2.3.

Fig. 10
Results of the Monte Carlo simulation with the final determined parameters; 0.5% of the simulations carried out are shown in the complex number plane
Fig. 10
Results of the Monte Carlo simulation with the final determined parameters; 0.5% of the simulations carried out are shown in the complex number plane
Close modal

Finally, a sensitivity analysis determines the influence of the uncertainty intervals δlin and δhys on the distribution of magnitude and phase in the frequency domain. From this, the propagation of systematic uncertainty due to linearity and hysteresis errors into the frequency domain is determined. Both boundaries are varied in steps of 0.001 within a range from 0 to 0.01. A full factorial design is used.

The 95% confidence interval cannot be determined from the standard deviation, as the distribution in the frequency domain could not be described by known distributions in the time domain. Instead, the two quantiles Q2.5 and Q97.5 are determined, and thus the width of the 95% confidence interval is extracted. The uncertainty is given by
(16)
(17)

Figure 11 shows the results for the uncertainty. Both metrics are viewed against the input parameters, i.e., the uncertainty due to linearity and hysteresis errors. The independence of magnitude and phase is confirmed. The linearity error only leads to a change in the magnitude of the Fourier coefficient, while the hysteresis error only affects the phase. The combination of both errors has no impact on this result. The two errors are independent in the frequency domain, which is not the case in the time domain, see Fig. 6.

Fig. 11
Influence of the uncertainty intervals of the linearity and hysteresis errors on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Fig. 11
Influence of the uncertainty intervals of the linearity and hysteresis errors on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Close modal
From the two linear correlations in Fig. 11, simple propagation rules for the uncertainty due to linearity error δlin and hysteresis error δhys into the frequency domain can be derived
(18)
(19)
The overall systematic uncertainty resulting from all four considered error sources of the magnitude in the frequency domain is
(20)
whereas the systematic uncertainty of the phase yields
(21)
Applying the results to the given exemplary force sensor, cf. Table 2, we obtain the systematic uncertainty in the frequency domain
(22)
(23)

4 Application to Experimental Parameter Identification

Until now, the presented uncertainty quantification is limited to the propagation of the systematic sensor uncertainty of an arbitrary sensor. This section presents the results of the aforementioned methodology applied to a complex system comprising various sensors, including eddy current and hall sensors. Subsequently, the uncertain sensor signals are employed to identify a system of complex system parameters, specifically the stiffness, damping, and inertia coefficients of a liquid-filled annular gap, cf. upper part of Fig. 12, within the frequency domain. Accordingly, the annular gap flow test rig at the Chair of Fluid Systems at Technische Universität Darmstadt, the employed identification procedure, and the previously utilized uncertainty quantification are briefly described. The uncertainty quantification presented in Sec. 3.2 is then employed as a substitute for the former one. Ultimately, the calculated uncertainty of the identified parameter is compared to a previously presented method by Kuhr [14].

Fig. 12
The annular gap test rig at the Chair of Fluid Systems at Technische Universität Darmstadt, cf. Kuhr [14]
Fig. 12
The annular gap test rig at the Chair of Fluid Systems at Technische Universität Darmstadt, cf. Kuhr [14]
Close modal

4.1 Experimental Setup, Identification Procedure, and Approach to Uncertainty Quantification.

Figure 12 shows the annular gap flow test rig for the experimental investigation of the static and dynamic force characteristics of generic fluid-filled annuli, cf. Kuhr et al. [27] and Kuhr [14]. In this setup, either journal bearings or annular seals are investigated not only for their static characteristics, e.g., their induced load carrying capacity or generated pressure drop across the gap, but also for their dynamic behavior, i.e., their stiffness, damping, and inertia properties. The test rig is primarily composed of two active magnetic bearings, which serve the rotor support as well as an intrinsic displacement, excitation, and force measurement system. The displacements of the rotor are measured using two eddy current sensors in an 90deg arrangement at both the entrance and exit of the lubrication gap. Each sensor has a user-defined measuring range of 1mm and an uncertainty in linearity of δx,lin<±2.4μm. As the manufacturer does not specify any further uncertainties regarding bias, sensitivity, and hysteresis, these are assumed to be negligible or included in the linearity uncertainty. The static and dynamic forces generated by the magnetic bearing are obtained by measuring the magnetic flux density B˜ in the air gap between the rotor and the active magnetic bearing by using hall sensors on each pole of each electromagnet. The hall sensors have a measuring range of ±2T and an uncertainty in linearity of δB,lin=±0.1%B. Similar to the eddy current sensor, the remaining sensor uncertainties in bias, sensitivity, and hysteresis are assumed to be negligible. The force applied per pole is proportional to the magnetic flux density FHB2. While the identification procedure and uncertainty analysis used in the initial reference are outlined briefly below, a more comprehensive description of the test rig, its calibration, and the actual measurement task can be found in the work of Kuhr et al. [27] and Kuhr [14].

Since the stiffness, damping, and inertia characteristics of annular gaps cannot be directly measured, they must be identified using parameter estimation methods, such as linear and quadratic regressions. The approach detailed by Kuhr [14] employs four linearly independent whirling motions to excite the rotor monofrequently at predefined precessional frequencies ωi and amplitudes Δe and Δγ: (i) translational excitation in the direction of rotor rotation; (ii) translational excitation against the direction of rotor rotation; (iii) angular excitation in the direction of rotor rotation; and (iv) angular excitation against the direction of rotor rotation. The induced dynamic forces and moments acting on the rotor FX, FY and MX, MY along with the rotor motion in the four degrees-of-freedom, X,Y and αX,βY, are measured using the aforementioned hall and eddy current sensors. All sensor data are collected simultaneously via a multifunction I/O module operating at a sampling rate of 6000Hz per channel for a duration of 5 s.

To identify the stiffness, damping, and inertia characteristic of the annular gap, the equation of motion
(24)
is first transferred into the frequency domain, yielding the complex equation of motion
(25)

Here, Kij are the complex stiffness coefficients depending on the precessional frequency ω. F=DFT(F) and M=DFT(M) are the Fourier transformations of the induced forces and moments on the rotor, whereas Di=DFT(i˜) are the Fourier transforms of the translational and angular excitation amplitudes.

The real part of the complex stiffness R(Kij) contains the stiffness and inertia coefficient, whereas the imaginary part I(Kij) contains the damping coefficients
(26a)
(26b)

The stiffness and inertia coefficients are obtained by applying a quadratic regression if the rotor is excited at different precessional frequencies ωi. In contrast, a linear regression gives the damping of the annulus.

To quantify the measurement uncertainty of the identified parameters, Kuhr [14] follows the guidelines set forth in the GUM [7]. However, for the propagation from the time into the frequency domain, i.e., from Eq. (25) to Eq. (26), it is assumed that the total measurement uncertainty of the signal δ(x), i.e., the Euclidean norm of the systematic δsys(x) and stochastic measurement uncertainty δstat(x), can be assigned to each discrete value of the recorded time signal, i.e., each of the 6000 samples. Furthermore, it is assumed that the overall measurement uncertainty can be regarded as an uncertainty in the sensitivity, i.e., the slope of the characteristic, and the real and imaginary part of the Fourier coefficients have the same magnitude of uncertainty after applying the discrete Fourier transform. This leads to a straightforward yet highly conservative propagation law for the propagation of uncertainty from the time domain to the frequency domain. The uncertainty of the stiffness, damping, and inertia coefficients is then calculated following Gaussian uncertainty propagation based on the equations of the identification procedure, i.e., the quadratic and linear regression.

To apply the simple propagation laws presented in this paper, to the complex parameter identification for the stiffness, damping, and inertia coefficients, the assumptions employed by Kuhr are discarded, and the genuine sensor signals are regarded as analogous to the generic signal utilized in the Monte Carlo analysis. Strictly speaking, the simple propagation laws, i.e., Eqs. (18) and (19), are applied directly at the transformation from the time into the frequency domain, i.e., from Eq. (25) to Eq. (26). The subsequent procedure of uncertainty propagation of the transformed signals according to GUM remains unaffected.

4.2 Comparison of the Uncertainty of the Identified Parameters.

Figure 13 depicts the identified coefficients for an annular gap with length L=92.95mm at a water temperature of 35°C, spanning a relative eccentricity range 0.1ε:=e˜/h¯˜0.8. Here, e˜ represents the eccentricity, and h¯˜ denotes the mean gap height of the annulus. The operating conditions of the annulus, i.e., the volume flow through the gap, the angular velocity of the rotor, as well as the velocity of the fluid flow in front of the annulus, are kept constant throughout the measurements. The corresponding dimensionless operating conditions, i.e., the modified Reynolds number, the flow number, and the preswirl ratio, are Reφ*=0.031, ϕ=0.7, and Cφ|z=0=0.5. Here, the makers represent the experimental results including uncertainty bound for a confidence interval of 95%. The solid lines are model predictions of the rotordynamic coefficients obtained by the clearance-averaged pressure model [28]. The left side of the figure shows the mean values of the results, while the figures on the right side compare the calculated measurement uncertainty obtained by the method used by Kuhr [14] and the one presented here.

Fig. 13
Influence of the improved uncertainty quantification for the identification of the dynamic characteristics, i.e., selected stiffness, damping, and inertia coefficients, of an annular gap. Left: Mean values of the identified parameters by Kuhr [14]; right: uncertainty interval given by Kuhr [14] compared to the presented results.
Fig. 13
Influence of the improved uncertainty quantification for the identification of the dynamic characteristics, i.e., selected stiffness, damping, and inertia coefficients, of an annular gap. Left: Mean values of the identified parameters by Kuhr [14]; right: uncertainty interval given by Kuhr [14] compared to the presented results.
Close modal

Considering that 16 stiffness, damping, and inertia coefficients are determined in the original publication, only the change in the uncertainty interval of the direct and cross-coupling stiffness KYYandKYX, the direct and cross-coupling damping CXX and CXY, as well as the direct and cross-coupling additional mass MXX and MXY is given here.

It is evident that the presented uncertainty quantification has a severe impact on the size of the error bars on the dynamic coefficients. In contrast to the previous method, where the results are clearly visible and comparable to the actual measured value, e.g., for the identification of the direct and cross-coupled added masses MXX and MXY, the results presented here are almost indiscernible, obscured by the markers. The new uncertainty quantification is a more realistic one compared to the worst-case scenario used previously. This is mainly due to the reduction of the impact of the systematic uncertainty. Due to the long measurement period of 5 s, the influence of the statistic measurement uncertainty is negligible.

5 Conclusion

This paper investigates the propagation of four different categories of systematic errors into the frequency domain. Overall, the two research questions:

  • What types of systematic error affect magnitude and phase in the frequency domain?

  • Are there simple rules for propagating systematic measurement uncertainty into the frequency domain?

are investigated. Regarding the first of the two questions, offset errors affect only static measurements, i.e., Ω=0, and can therefore easily be transferred into the frequency domain. Since Fourier analysis are typically used for dynamic transfer functions, offset errors can be neglected. Linearity errors and sensitivity errors affect only the magnitude of a Fourier coefficient, whereas hysteresis errors result only in the phase of the Fourier coefficient. The linearity error and the hysteresis error are completely independent. This also applies if the assumptions about the distributions and the error model are changed. Another important note is that the uncertainty does not disappear at high sampling rates, proving that the model is a true systematic uncertainty model. Regarding the second research question, simple conservative propagation relations are provided for the uncertainty of the linearity error and the hysteresis error into the frequency domain
(27)
(28)

Author Contribution Statement

Manuel Rexer: Conceptualization, Methodology, Software, Formal analysis, Visualization, Writing—original draft, Writing—review and editing, Validation. Peter F. Pelz: Funding acquisition, Resources. Maximilian M. G. Kuhr: Conceptualization, Investigation, Supervision, Visualization, Writing—original draft, Writing—review and editing.

Funding Data

  • Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) (Project No. 57157498—SFB805; Funder ID: 10.13039/501100001659).

  • Federal Ministry for Economic Affairs and Energy (BMWi) due to an enactment of the German Bundestag (Grant No. 03EE5036B; Funder ID: 10.13039/501100006360).

  • Industrial collective research programme supported by the Federal Ministry for Economic Affairs and Energy (BMWi) through the AiF (German Federation of Industrial Research Associations e.V.) due to an enactment of the German Bundestag (IGF No. 21029 N/1; Funder ID: 10.13039/501100002723).

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Appendix: Further Influences of the Input Signal on the Results

Fig. 14
Influence of the length L of the test signal on the standard deviation of the magnitude and phase of the Fourier coefficient Y(Ω)
Fig. 14
Influence of the length L of the test signal on the standard deviation of the magnitude and phase of the Fourier coefficient Y(Ω)
Close modal
Fig. 15
Influence of the offset o of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Fig. 15
Influence of the offset o of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Close modal
Fig. 16
Influence of the phase Φ of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Fig. 16
Influence of the phase Φ of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Close modal
Fig. 17
Influence of the excitation frequency Ω of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Fig. 17
Influence of the excitation frequency Ω of the test signal on the standard deviation of magnitude and phase of the Fourier coefficient Y(Ω)
Close modal

References

1.
Pelz
,
P. F.
,
Groche
,
P.
,
Pfetsch
,
M. E.
, and
Schaeffner
,
M.
, eds.,
2021
,
Mastering Uncertainty in Mechanical Engineering
(Springer Tracts in Mechanical Engineering), 1st ed.,
Springer
,
Cham, Switzerland
.
2.
Pintelon
,
R.
,
2012
,
System Identification: A Frequency Domain Approach
, 2nd ed.,
Wiley and IEEE Press
,
Hoboken, NJ and Piscataway, NJ
.
3.
Isermann
,
R.
,
2005
,
Mechatronic Systems: Fundamentals
,
Springer
,
London
.
4.
Puff
,
M.
, and
Pelz
,
P. F.
,
2009
,
Entwicklung einer Prüfspezifikation zur Charakterisierung von Luftfedern
(FAT-Schriftenreihe),
Verband der Automobilindustrie (VDA)
,
Berlin, Germany
.
5.
Hedrich
,
P.
,
Lenz
,
E.
,
Brötz
,
N.
, and
Pelz
,
P. F.
,
2018
, “
Active Pneumatic Suspension for Future Autonomous Vehicles: Design, Prove of Concept and Hardware-in-the-Loop Simulations
,”
Proceedings of the 11th International Fluid Power Conference on Fluid Power Networks
, Aachen, Germany, Mar. 19–21, pp.
352
365
.10.18154/RWTH-2018-224797
6.
Hedrich
,
P.
,
Lenz
,
E.
, and
Pelz
,
P. F.
,
2018
, “
Minimizing of Kinetosis During Autonomous Driving
,”
ATZ Woldwide
,
120
(
7–8
), pp.
68
75
.10.1007/s38311-018-0081-0
7.
Joint Committee for Guides in Metrology
,
2008
,
Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement
, Bureau International des Poids et Mesures (BIPM), Sèvres, France.
8.
Betta
,
G.
,
Liguori
,
C.
, and
Pietrosanto
,
A.
,
2000
, “
Propagation of Uncertainty in a Discrete Fourier Transform Algorithm
,”
Measurement
,
27
(
4
), pp.
231
239
.10.1016/S0263-2241(99)00068-8
9.
Medina
,
L.
,
Ruiz
,
R.
, and
Díaz
,
S.
,
2008
, “
A Simple Approach to Determine Uncertainty Bounds on Bearing Rotordynamic Coefficients Identification
,”
ASME
Paper No. GT2008-51200.10.1115/GT2008-51200
10.
Medina
,
L. U.
,
Diaz
,
S. E.
, and
LLatas
,
I.
,
2010
, “
Measurement Uncertainty Propagation in Frequency Domain Parameter Identification of Mechanical Systems
,”
ISMA 2010
, Leuven, Belgium, Sept. 20–22, pp.
5129
5142
.https://past.isma-isaac.be/downloads/isma2010/papers/isma2010_0083.pdf
11.
Eichstädt
,
S.
, and
Wilkens
,
V.
,
2016
, “
GUM2DFT—A Software Tool for Uncertainty Evaluation of Transient Signals in the Frequency Domain
,”
Meas. Sci. Technol.
,
27
(
5
), p.
055001
.10.1088/0957-0233/27/5/055001
12.
Moreland
,
J. A.
,
Childs
,
D. W.
, and
Bullock
,
J. T.
,
2018
, “
Measured Static and Rotordynamic Characteristics of a Smooth-Stator/Grooved-Rotor Liquid Annular Seal
,”
ASME J. Fluids Eng.
,
140
(
10
), p.
101109
.10.1115/1.4040762
13.
Eichstädt
,
S.
,
Link
,
A.
,
Harris
,
P.
, and
Elster
,
C.
,
2012
, “
Efficient Implementation of a Monte Carlo Method for Uncertainty Evaluation in Dynamic Measurements
,”
Metrologia
,
49
(
3
), pp.
401
410
.10.1088/0026-1394/49/3/401
14.
Kuhr
,
M.
,
2023
, “
Identification of the Dynamic Force and Moment Characteristics of Annular Gaps Using Linear Independent Rotor Whirling Motions
,”
Mech. Syst. Signal Process.
,
187
, p.
109936
.10.1016/j.ymssp.2022.109936
15.
Vasilevskyi
,
O. M.
,
2015
, “
A Frequency Method for Dynamic Uncertainty Evaluation of Measurement During Modes of Dynamic Operation
,”
Int. J. Metrol. Qual. Eng.
,
6
(
2
), p.
202
.10.1051/ijmqe/2015008
16.
Hessling
,
J. P.
,
2008
, “
Dynamic Metrology—An Approach to Dynamic Evaluation of Linear Time-Invariant Measurement Systems
,”
Meas. Sci. Technol.
,
19
(
8
), p.
084008
.10.1088/0957-0233/19/8/084008
17.
Hessling
,
J. P.
,
2011
, “
Propagation of Dynamic Measurement Uncertainty
,”
Meas. Sci. Technol.
,
22
(
10
), p.
105105
.10.1088/0957-0233/22/10/105105
18.
Elster
,
C.
, and
Link
,
A.
,
2008
, “
Uncertainty Evaluation for Dynamic Measurements Modelled by a Linear Time-Invariant System
,”
Metrologia
,
45
(
4
), pp.
464
473
.10.1088/0026-1394/45/4/013
19.
Elster
,
C.
, and
Link
,
A.
,
2009
, “
Analysis of Dynamic Measurements: Compensation of Dynamic Error and Evaluation of Uncertainty
,”
Advanced Mathematical and Computational Tools in Metrology and Testing
, World Scientific Publishing, Singapore, pp.
80
89
.
20.
Zeier
,
M.
,
Hoffmann
,
J.
, and
Wollensack
,
M.
,
2012
, “
Metas.UncLib—A Measurement Uncertainty Calculator for Advanced Problems
,”
Metrologia
,
49
(
6
), pp.
809
815
.10.1088/0026-1394/49/6/809
21.
Elhoushi
,
M.
,
Georgy
,
J.
,
Noureldin
,
A.
, and
Korenberg
,
M. J.
,
2017
, “
A Survey on Approaches of Motion Mode Recognition Using Sensors
,”
IEEE Trans. Intell. Transp. Syst.
,
18
(
7
), pp.
1662
1686
.10.1109/TITS.2016.2617200
22.
Tränkler
,
H.-R.
, and
Fischerauer
,
G.
,
2014
,
Das Ingenieurwissen: Messtechnik
,
Springer Berlin Heidelberg
,
Berlin, Heidelberg, Germany
.
23.
HBK,
U2B Kraftaufnehmer Datenblatt
,” Hottinger Brüel and Kjær, Darmstadt, Germany, accessed Feb. 21, 2025, https://www.hbm.com/fileadmin/mediapool/hbmdoc/technical/B00389.pdf
24.
Brigham
,
E.
,
1974
,
The Fast Fourier Transform
(Prentice-Hall Signal Processing Series),
Prentice Hall
, Englewood Cliffs, NJ.
25.
Metropolis
,
N.
,
Rosenbluth
,
A. W.
,
Rosenbluth
,
M. N.
,
Teller
,
A. H.
, and
Teller
,
E.
,
1953
, “
Equation of State Calculations by Fast Computing Machines
,”
J. Chem. Phys.
,
21
(
6
), pp.
1087
1092
.10.1063/1.1699114
26.
Barbu
,
A.
, and
Zhu
,
S.-C.
,
2020
,
Monte Carlo Methods
,
Springer Singapore
,
Singapore
.
27.
Kuhr
,
M. M. G.
,
Lang
,
S. R.
, and
Pelz
,
P. F.
,
2022
, “
Static Force Characteristic of Annular Gaps—Experimental and Simulation Results
,”
ASME J. Tribol.
,
144
(
11
), p.
111804
.10.1115/1.4054792
28.
Kuhr
,
M. M. G.
,
Nordmann
,
R.
, and
Pelz
,
P. F.
,
2023
, “
Dynamic Force and Moment Characteristics of Annular Gaps—Simulation Results and Evaluation of the Relevance of the Tilt and Moment Coefficients
,”
ASME J. Tribol.
,
145
(
1
), p.
011801
.10.1115/1.4055180