![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
Knowledge of the load parameters of wind turbine rotor is the key problem to eliminate the misalignment of wind turbine. In this paper, an inverse method for identifying the load parameters with imprecise information is investigated by using the probability and interval variables. The interval variables are dealt with interval analysis, only the bounds of the structure response are needed avoiding time-consuming inner loop for structure response prediction. The probability variables are dealt with a proposed iteration method combined micro genetic algorithm with the search interval advance and retreat method. Firstly, at the midpoint of the interval variables, the mean of the load parameters is identified based on the proposed iteration method; secondly, the first-order Taylor expansion is employed to gain the boundary of rotor dynamic response from predicting structure response strip derived from the existence of interval variables; thirdly, the interval of the load parameters is identified via the boundary of rotor dynamic response based on the proposed iteration method. Eventually, the accuracy and efficiency of the proposed inverse method are verified by one numerical example.
1. Introduction
The knowledge of load parameters that lead to the misalignment of the wind turbine rotor is always required to eliminate the fan misalignment [1–3]. However, its measurement with force actuators is a difficult task in many instances. Instead, vibration responses can be conveniently measured. The technique of load parameters identification, an important indirect method aimed at estimating load parameters via measured vibration responses, is often preferred in many mechanical engineering practices [4]. Most of the existing load parameters identification is executed for deterministic problems, where the structure parameters and measured vibration responses are completely known. However, due to manufacturing errors, measurement inaccuracy, the discreteness of thin silicon steel sheet material and incomplete knowledge, this information is only partly known in the real wind turbine rotor. Namely, the measured vibration responses are also variable and random mostly corrupted by the measurement errors; meanwhile, the vibration responses are a distribution strip rather than a single distribution because of the existence of the uncertain structure parameters. This kind of problem is unable to be sorted as in traditional deterministic load parameters identification. We term them as mixed uncertainty problems in this study. Due to uncertainties, the identified load parameters via vibration response strip are a distribution strip rather than a single distribution, where the inner loop is uncertain structure response prediction and the outer loop is an inverse calculation in the solving process. The accuracy of identified results can hardly be evaluated using traditional deterministic methods. Therefore, uncertainty-based load parameters identification is an emerging research topic and a more computationally challenging exercise as well.
As the most popular uncertainty modelling strategy, interval method, probability method and evidence theory method [5–7] have been widely applied in engineering for uncertain parameter identification. Hang Wei [8] developed a probability and interval estimation method to assess the probability distributions of unknown parameters through the Monte Carlo Markov Chain method (MCMC) in Bayesian theory [9–11] considering their prior information. The maximum likelihood estimation method [12,13] also has been successfully developed to identify the probability distributions of unknown parameters without considering the prior information. The maximum likelihood estimation method is to determine the unknown parameters through the relationship between the calculated and the measured response of investigated structure by the iteration method, the parameters corresponding to the maximum likelihood function are taken as the identified results, and the corresponding confidence interval is calculated [12]. Compared to MCMC, the maximum likelihood estimation method has a good prospect in the wind turbine rotor system, where the misalignment fault is the coupling influence of multiple load parameters [14], it is difficult for engineers to understand the prior information of load parameters based on experience. However, a commonly encountered difficulty in the maximum likelihood estimation method is that the phenomenon of iterative divergence rather than convergence often occurs because of the sensitivity of the initial value [13]. Therefore, the iterative method should be proposed for identifying the load parameters with probabilistic uncertainty.
Uncertain parameters are mainly described by interval or probability. Structural response strip coming from the interval uncertainty makes the nested solution of the load parameters identification problem have a large amount of computation. The maximum likelihood method to deal with the probability uncertainty is often faced with low efficiency due to the complexity of the likelihood function. Meanwhile, in practical engineering, only the bounds of load parameters are needed to eliminate the fan misalignment. In order to solve the above problems, the load parameters identification under uncertainty is improved in this paper. An inverse method is presented to identify the load parameters of the wind turbine rotor system involving probability and interval variables. The upper and lower bound of structure response coming from interval uncertainty are predicted by the Taylor expansion method and the interval analysis method; A proposed iteration method based on the search interval advance and retreat method (ARM) [15] and micro genetic algorithm (uGA) [16] is used to deal with the probabilistic uncertainty in load parameter identification of the wind turbine rotor system. In the proposed iteration method, the initial load parameter starts from zero, the iterative increment is obtained by the sensitive matrix method and maximum likelihood method, the search space is obtained by ARM, and uGA is used to minimize the deviation between the measured and calculated structure response in the corresponding search space. If the deviation is within the predefined toleration, the set of load parameters is considered to be the total solution to the load parameter identification problem, and the solution procedure ends. The proposed iteration method avoids the stringent requirements of the traditional iterative method on the search space and improves the identification accuracy of the maximum likelihood method.
2. The basic idea of load parameter identification with uncertainty
The dynamic equation of a wind turbine rotor system with misalignment fault [17] is
(1)
(1) where q,
,
are the displacement, velocity and acceleration response vectors, M,
and C respectively represent the whole structural mass, stiffness and damping matrix which can be established by using the Euler–Bernoulli beam theory in [18];
denotes the equivalent stiffness matrix of the coupling misalignment; G denotes the Gravity vector; Q denotes the unbalanced force vector and
denotes the misaligned load vector.
If the model-related matrix and force-related vector
are complete known constant;
is a n-dimensional vector of input parameters representing the external causes;
is an m-dimensional vector of output parameters representing the effects. A structural mechanics problem among the input
and the output
in Equation (1) can be conceptually expressed as follows [19]:
(2)
(2) where
represents the model-related matrix
;
represents force-related vector
;
is an m × n system characteristics matrix representing the translation process from input
to output q. It is derived from the physical laws underlying the problem under consideration. Matrix
may be dependent on the input
for the wind turbine rotor system with misalignment fault.
That is, a specific input load parameter corresponds to a corresponding output response
. According to the type of information be sought in the solution procedure from Equation (2), we traditionally treated the structural mechanics problem among the input
and the output q as either a forward problem or an identification problem. The forward problem is to find out the effect (
) from the given cause (
) combined with the known model-related matrix
and force-related vector
utilizing dynamic analysis solution techniques including FEM [20], generalized differential quadrature method [21], transfer matrix method [22], etc. The identification problem is to determine the unknown cause (
) from the given effect (
) combining with the known model-related matrix
and force-related vector
utilizing computational reverse techniques including optimization method [23], regularization method [24], etc. However, a commonly encountered difficulty in solving identification problems is that it is often unable to provide a complete known
, model-related matrix
and force-related vector
because of the existence of measurement noise and model uncertainty. That is, the effect (
) and the model-related matrix
and force-related vector
are only partly known. We term the problem is treated as a mixed uncertainty problem in this study. A deterministic identification method in a traditional sense is difficult to quantitatively evaluate the credibility of the identified results. Thus, the identified result should also present a certain extent of uncertainty.
To deal with these difficulties in identifying the load parameters of the wind turbine rotor, a mixed uncertainty analysis framework is suggested to deal with the uncertain parameters coming from the model-related matrix
, force-related vector
and the effect (
) in this paper. This brings up a new approach to solve the mix uncertainty problem, where the probability and interval theory is introduced to deal with the uncertain parameters, all the uncertain parameters in the model-related matrix
, force-related vector
and the effect (
) are as either probability or interval variables. Referring to the mathematical expression of the forward problem in [25], Equation (2) involving probability and interval variables, aimed at estimating the displacement response vectors
via a given vector
, can be expressed as follows:
(3)
(3) where
and
respectively represent a vector that collects all interval and probabilistic uncertain parameters from the model-related matrix
and force-related vector
. Probabilistic uncertain parameters obey a zero-mean Gaussian distribution with a variance
; Y is a characteristics matrix representing the uncertainty translation process from input
to output
for a wind turbine rotor system with interval and probabilistic uncertain parameters.
denotes the calculated structure response from the given input
and the established matrix Y, the
will be a response strip because of the interval uncertain parameters
;
denotes the measured response in the experiment with noise vector
drawn from a zero-mean Gaussian distribution with a variance
, the measurement noise from an experiment is described by probabilistic variables.
respectively denote calculated structure response, the measured response and measurement noise at the
specific measurement location.
is a characteristics matrix representing the uncertainty translation process from input
to output
at the
specific measurement location.
Referring to Bayesian identification [8] for uncertain structures with Gaussian distribution, the general Bayesian identification for Equation (3), aimed at estimating the via giving the m-point measured response
can be written as follows:
(4)
(4) where
is n-unknown load parameters, which would be determined;
denotes the prior probability density of
, and is often considered as constants in given prior parameter spaces;
is defined as the data likelihood, and can be regarded as a measure of how well fit the measured response data
;
is the joint posterior probability density distribution (PPDF) of
; m denotes the number of the measured point.
As the mechanism of the Bayesian identification approach, the prior knowledge of is updated with the data likelihood of
to achieve the current knowledge of
. In the present engineering practice, we are typically interested in the marginal probability density distribution (PPDF) of each unknown load parameter. Theoretically, this would require integrating Equation (4) over the other j−1 load parameters. The marginal PPDF
of the
unknown load parameters
is
(5)
(5) where
denotes the multidimensional integral over all parameters other than
; a is constant.
Equation (5) show that the uncertainty-based load parameter identification is a double-loop computational procedure, where the inner loop is executed for uncertain response prediction, the outer loop is to determine the load parameter via response
. But there are three main problems in this solution procedure. First, the marginal PPDF
will form a distribution strip rather than a single distribution due to the response strip resulting from the existing interval parameters
. Second, the multidimensional integral cannot be done directly in an implicit form due to the complexity of the likelihood function. Third, the prior probability density
is often unable to be provided as required in a traditional Bayesian approach. The traditional MCMC of the Bayesian approach to determine the load parameter
by solving Equation (4) will involve a complex multi-layer nesting solution affecting the computational efficiency. Therefore, an appropriate inverse method should be developed to deal with these difficulties.
3. A proposed inverse method with mixed uncertainty
According to the interval mathematics theory, due to the existence of interval uncertain parameters , the marginal posterior distribution
will be formed the following interval vector instead of a real number vector, the marginal PPDF in Equation (5) can be written as follows:
(6)
(6) where
denotes a closed bounded set of real numbers, the superscript s denotes an interval, the superscripts l and r denote the lower and upper bounds of interval, respectively.
and
are the lower and upper bounds of
, respectively.
and
respectively denote the lower and upper response bounds caused by all possible structural responses because of interval uncertain parameters
in terms of Equation (3);
and
respectively represent the lower and upper bounds of interval uncertain parameters
.
In this situation, the mean of the unknown load parameters is no longer the only value but has a variety of possible results, and all the results are concentrated in the interval vector. Similarly, the confidence interval of unknown load parameters is no longer the only interval but has a variety of possible results, and all the interval results are concentrated in the interval vector. The mean value and confidence interval of unknown load parameters obtained by Bayesian parameter identification are interval forms. According to the Bayesian parameter identification method, in order to obtain the maximum and minimum values of mean and confidence interval of unknown load parameters, all possible values of interval uncertain parameters
should be selected sequentially, and, the MCMC method is used to sample in the joint PPDF strip, and the maximum and minimum values of mean and confidence interval are searched in the identification results. The process of Bayesian parameter identification will involve complex multi-layer nesting solution, and the forward problem calculation in Equation (2) needs to be called repeatedly, which will inevitably lead to low efficiency of solution.
In this paper, a proposed inverse method combined with rotor dynamic response prediction and a proposed iteration method is introduced to improve the Bayesian parameter identification method, and the boundary of structural response strip coming from interval uncertain parameters is predicted by Taylor expansion and interval mathematical theory. Combining the predicted structure response bound together with the probability uncertain parameters
into the proposed iteration method to calculate the corresponding mean and interval of the unknown load parameters.
3.1. The boundary of rotor dynamic response prediction with interval uncertainty
The interval uncertain parameters is described in the form of midpoint
-radius
[26]. Which means that the possible values of each interval parameter belong to an interval. Thus, for a specific
, the possible values of the output vector
will form the following interval vector instead of a single vector:
(7)
(7) where response strip
is an interval vector.
and
are the lower bound and the upper bound of response strip
. Equation (7) created a new deterministic relation between the input load parameters
and the output interval vector
, namely a specific
corresponds to a deterministic interval vector
In other words, the unknown load parameters
can be identified via interval vector
.
When the level of uncertainty is small, the lower bound and the upper bound
can be predicted directly and explicitly by first-order Taylor expansion at the midpoint
of the interval uncertain parameters
(8)
(8) where
denotes the
interval uncertain parameters, o denotes the number of interval uncertain parameters
. The response
can be obtained directly combining the mean
together with the midpoint
of the interval uncertain parameters and the established matrix Y according to Equations (8). The partial differential equation
is transformed into an algebraic equation
using the difference method, the response increment
is obtained by a small disturbance
at the midpoint
. The output vector
,
is predicted through Equations (8) based on the interval analysis, which avoids the problem of time-consuming inner loop for structure response prediction.
According to the monotonicity analysis theory [27], the maximum and minimum values of the mean value of the unknown load parameters must correspond to the boundary of the uncertain parameters in the interval of the model. Namely, the MCMC method is used to sample and solve the boundary of
corresponding to the structural response boundary and the mean
of the unknown load parameters obtained by MCMC sampling corresponding to the midpoint
can be used as the midpoint of the mean interval. The improved framework of the Bayesian identification method based on the interval analysis can obviously reduce the computational burden. However, the problem of mass sampling also exists to obtain a high-precision joint PPDF space in the sampling process of MCMC due to probability distribution coming from the probabilistic uncertain parameters
. The maximum likelihood method can transform the Bayesian identification method into an iterative method by introducing the sensitivity matrix method that explicitly dealing with the relation between the input unknown parameter vector
and the probabilistic uncertain parameters
. The explicit method avoids the mass sampling process of MCMC, but the iterative efficiency is low coming from the unknown priori information of load parameters and the sensitivity of the initial value. Therefore, the iterative method should be proposed for the identifying the load parameters with probabilistic uncertainty.
3.2. A proposed iteration method with the probabilistic uncertainty
To overcome this iterative efficiency, a proposed iteration method combined with uGA and ARM is introduced, where the initial value starts from zero; The increment of load parameters
is obtained based on the sensitivity matrix method, which can transform the probability uncertain parameters into explicit solvable form avoiding MCMC sampling; The search interval in each iterative step is obtained by ARM; uGA is used to gradually correct the solution of load parameters
in an iterative procedure until the error between the measured response and calculated response is within the specified toleration.
Assume the variation and
, combining them together with
into Equations (3) to calculate the variation
and
of the output response of the corresponding k-measured point and u-concern point. Calculate the sensitivity matrix
centering on the present parameters
using the sensitivity matrix-based method [19], draw the sensitivity matrix
from the sensitivity matrix
according to Equation (9). The increment
of the unknown load parameter
is obtained according to the maximum likelihood method shown in Equation (10).
(9)
(9) where the subscripts k and u denote the number of measured points and concern points;
and
denote the variation of the output response
and
corresponding to the measured point and concern point, respectively.
is the variation of the unknown load parameter
;
is the variation of the known probability uncertain parameter
.
(10)
(10) where
,
,
,
,
are variances, and
are the deviations of the
parameter in the known probability uncertain parameter
and measured response
.
,
.
and
denotes the unit matrix and zero matrix, respectively.
Taking the increment of the unknown load parameter value
as the step size in the
iterative step, the corresponding search interval
is obtained by using the ARM shown in Equation (11). In the
iterative step,
is the initial step size, and
,
,
is the initial point, intermediate point and end point taken by the unknown load parameter
, respectively.
is the mean square deviation between the output response
(the measured rotor dynamic response
or predicted
,
) in the measured point and the corresponding calculated result
through Equation (11).
and
are the lower and upper bounds of the obtained search interval, respectively.
The uGA is used to search the unknown load parameter value corresponding to the minimum error between the calculated and measured response in the search interval
. If the error is within the predefined toleration, the set of load parameters
is considered to be the total solution to the problem, and the solution procedure ends. Otherwise, the load parameter value
enters the next iteration as the update of the new unknown load parameter value, and the solution procedure continues.
(11)
(11)
4. Load parameter identification process
Based on the above discussion, the proposed inverse method combining rotor dynamic response prediction together with the proposed iteration method can calculate the load parameter of the wind turbine rotor with interval and probability uncertainty according to the deterministic load parameter identification method. The solution procedure is described in Figure .
Step 1: The dynamic response of the wind turbine rotor system is measured. Determine the midpoint -radius
of the interval uncertain parameters
, the probability uncertain parameter
, the Gaussian distribution of the measurement error and the unknown load parameter
.
Step 2: Combining the midpoint of the interval uncertain parameters together with the measured rotor dynamic response
into the proposed iteration method to calculate the mean
Step 3: According to the interval analysis method, the upper bound and the low bound
of the dynamic response of the wind turbine rotor are predicted by solving Equations (8).
Step 4: Combining the lower bound of the interval uncertain parameters together with the predicted
into the proposed iteration method to calculate the lower bound
of the load parameters. similarly, obtain the upper bound
by combining the upper bound
together with the predicted
. As a result, the mean range of unknown load parameters
is obtained.
5. Identification of load parameters of the wind turbine rotor
5.1. Recognition model
The rotor system of the wind turbine is given in this section to validate the proposed method. The parameters of the conical rotor model with three equivalent disk and two sliding bearings are shown in Figure and Table .
Table 1. The parameters of the conical rotor system.
In this case, the axial load and torque load parameters are the unknown load parameters to be determined. The rotor elastic modulus
is an interval uncertain parameter. The length and density of the conical rotor (L, ) are probabilistic uncertain parameters with Gaussian distribution. The number of the measured point is 2, and the dynamic response (d1, d2) of the front and rear bearings can be measured with Gaussian measurement errors, and denote the output response
corresponding to the measured point. The dynamic response (d3, d4) of the front and rear disc is not easy to arrange the sensor and cannot be measured, denote the output response
corresponding to the concern point. Namely, there are both measurement noise and model uncertainty, the identification of load parameters
is a mixed uncertainty problem.
5.2. Recognition process
The interval uncertainty level is assumed as 3%, that is, ,
According to the identification process, the midpoint of the interval uncertain parameters
. It is assumed that the true value of load parameters
are given, respectively. Based on Equation (2), therefore, the true rotor dynamic response
can be easily obtained. In order to simulate the measured responses, noise-contaminated response
is used to identify the mean
of load parameters. It is assumed that the all-measurement results and probability uncertain parameter have Gaussian noise and obey the normal distribution with a deviation of 3%, namely, the probability uncertain parameter
are expressed as
,
, the noise contaminated response
are expressed as
,
. In this case, five individuals are selected in the uGA, the crossover factor is 0.5, the gene mutation rate is 0.003, and the maximum iterative step is 300.
5.2.1. Identification of the mean ![](//:0)
of load parameters
Table shows the calculation process of load parameter identification. The stop criterion constant is set to 2e-5, and the initial value of the given load parameter is The load parameter increment
is calculated from the sensitivity matrix, and the corresponding load parameter search interval is obtained according to the ARM. The axial load parameter interval is
and the torque parameter interval is
. In this interval, the load parameter
is obtained corresponding to the minimum error of 3.36 by the uGA, as shown in the first row of Table , and the error 3.36 is greater than the predefined toleration 2e-5, then take
as the updated initial value of the load parameter in the next iteration, and continue to run according to the proposed iteration method, as shown in the second row of Table . After 34 iterations, the error is 1.54e-05, which is within the predefined toleration and the calculation is over. The identified load parameter is
. The mean square deviation between the identified value and the real value is 0.06. As shown in Table , the proposed iteration method is feasible to identify the load parameters of the wind turbine rotor system with probabilistic uncertainty.
Table 2. Identification iterative process of load parameter with
involving probability uncertainty.
Table 3. Identified with
involving probability uncertainty.
To demonstrate the robustness of the proposed iteration method to probabilistic uncertainty, the load parameters identification results are compared with the Gaussian noise deviations of 1%, 3% and 5%. As shown in Figure and Table . It is known that the identification result is less affected by the Gaussian noise deviation, that is, the robustness of the proposed iteration method is good. Figure shows that the local optimal solution of the corresponding interval can gain in 350 iterative steps. The maximum iterative step of uGA in each interval (it is 300 in the case) is the main factor affecting the convergence speed of the proposed iteration method. However, according to the Bayesian approach, MCMC sampling with106 samples [8] must be carried out to compute all possible marginal PPDF in each interval in order to validate the accuracy of the identified results. So, the proposed iteration method improves its efficiency.
5.2.2. Identification of mixed uncertain load parameters
When the probabilistic Gaussian noise deviation is 3%, the interval parameters are considered to identify the load parameters of the wind turbine rotor system under three intervals uncertain levels. The upper and low bound of rotor dynamic response and
are predicted by the interval analysis method as shown in Table , and the load parameters identification intervals based on
and
are shown in Table . The identification iteration process is shown in Figure . The interval uncertain level affects the recognition result. The interval of load parameters identification results expands with the increase of the interval uncertain level.
Table 4. Load parameter identification result involving mixed uncertainty.
6. Conclusions
In order to reduce complex multi-layer nesting solutions coming from the interval and probability parameters, this paper proposes a new inverse analysis framework for load parameter identification in wind turbine rotor systems. A numerical example is given to demonstrate the efficiency and accuracy. Based on this study, the following conclusion can be drawn:
The uncertain parameters are characterized as interval and probability variables; the measurement noise and model uncertainty are easy to implement.
For conducting the wind turbine rotor response strip coming from the interval uncertain parameters, an interval analysis method is presented to predict the upper and lower bounds of the responses based on the first-order Taylor expansion.
For conducting the problem of mass sampling coming from the probability uncertain parameters in the Bayesian approach, the load identification problem is transformed into a parameter optimization problem to deal with probability variables. A proposed iteration method combined ARM with uGA is presented to determine the mean and interval of load parameters via the obtained response.
The level of uncertainty has a great influence on the identification results. It can be used as a reference for the uncertainty inverse problem in engineering.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Conflict of interest statement
The authors declared that they have no conflicts of interest to this work. We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.
Additional information
Funding
References
- Li HB. Fault diagnosis and treatment of rotor misalignment in wind turbines. Electr Eng. 2018;24:122–1146.
- Liu Y, Li YZ, Shi T, et al. Study on misalignment-rubbing coupling Fault of rotor system supported by oil film force. J Mech Eng. 2016;52(13):79–86.
- Long J, Wu JQ. Fault diagnosis of rotor misalignment for windmill generators. Noise Vibration Control System. 2013;33(3):222–225.
- Luo K, Shi WH, Qu JX. Multi-machine equivalent model parameter identification method for double-fed induction generator (DFIG)-based wind power plant based on measurement data. J Eng. 2017;2017:1550–1554.
- Sepahvandn K, Marburg S. Identification of composite uncertain material parameters from experimental modal data. Probab Eng Mech. 2014;37:148–153.
- Yang XQ, Xiong WL, Wang ZY, et al. Parameter identification of nonlinear multirate time-delay system with uncertain output delays. Trans Inst Meas Control. 2018;40(12):3498–3506.
- Wang C. Evidence-theory-based uncertain parameter identification method for mechanical systems with imprecise information. Comput Method Appl M. 2019;351:281–296.
- Zhang W, Liu J, Cho C, et al. A hybrid parameter identification method based on Bayesian approach and interval analysis for uncertain structures. Mech Syst Signal Process. 2015: 60–61. 853–865.
- Wang H, Zeng Y, Yu XC, et al. Surrogate-assisted Bayesian inference inverse material identification method and application to advanced high strength steel. Inverse Probl Sci Eng. 2016;24(7):1133–1161.
- Natalia CT, Tuckmantel WS, Cavalca KL, et al. Bayesian inference applied to journal bearing parameter identification. J Braz Soc Mech Sci Eng. 2017;39(8):2983–3004.
- Mao WG, Li JH, Liu GP. Unbalance parameters identification for a sliding bearing-rotor system based on genetic intelligent sampling technique. J Vib Eng. 2019;32(4):660–667.
- Aaron CZ, Langford BW, Martin FL, et al. Parameter identification of fluid line networks by frequency -domain maximum likelihood estimation. Mech Syst Signal Process. 2013;37(1-2):370–387.
- Zhang AL. The initial value research based on quasi-Newton method. Journal of Tai Yuan Normal University (Natural Science Edition). 2018;40(6):1378–1400.
- Xiao YC, Wang YJ, Ding ZT. The application of heterogeneous information fusion in misalignment fault diagnosis of wind turbines. Energies. 2018;1655(11):1–15.
- Zhang SE, Wang X, Wang K. Study of optimization of microwave cavity filter parameters based on a hybrid algorithm of trust region. Appl Mech Mater. 2015;713-715:1180–1183.
- Liu GP, Han X, Jiang C. An efficient multi-objective optimization approach based on the micro genetic algorithm and its application. Int J Mech Mater Des. 2012;8(1):37–49.
- Xia X, Ni W. A novel failure analysis and diagnosis method for hydraulic-turbine generator unit. J Vibroeng. 2016;18(6):3568–3580.
- Michael F, John ETP, Seamus DG. Dynamics of rotating machines. Cambridge (UK): Cambridge University Press; 2015.
- Liu GR, Xu YG, Wu ZP. Total solution for structure mechanics problems. Comput Methods Appl Mech Eng. 2001;19:1989–1012.
- Zhu HM, Chen WF, Zhu RP, et al. Dynamic analysis of a flexible rotor supported by ball bearings with damping rings based on FEM and lumped mass theory. J Cent South Univ. 2020;27(12):3684–3701.
- Al-Furjan MSH, Hatami A, Habibi M, et al. On the vibrations of the imperfect sandwich higher-order disk with a lactic core using generalize differential quadrature method. Compos Struct. 2021;257:1–18.
- Liu PF, Liu HJ, Wu Q, et al. Estimation of wheelset natural vibration characteristics based on transfer matrix method with various elastic beam models. Shock Vib. 2021;2021:1–15.
- Anne W. A fast subspace optimization method for nonlinear inverse problems in Banach spaces with an application in parameter identification. Inverse Prob. 2018;34(8):1–26.
- Smina D, Nabil S, Omar AA. A fractional Tikhonov regularization method for an inverse backward and source problems in the time-space fractional diffusion equations. Chaos, Solitons Fractals. 2021;150:111–127.
- Jiang C, Liu GR, Han X. A novel method for uncertainty inverse problems and application to material characterization of composites. Exp Mech. 2008;48:539–548.
- Wang WX, Xue H, Kong TZ. An efficient hybrid reliability analysis method for structures involving random and interval variables. Struct Multidiscipl Optim 61(1) 2020: 1–15.
- Jiang C H, Li WX X, et al. A hybrid reliability approach based on probability and interval for uncertain structures. J Mech Des. 2012;134(3):0310011–03100111.