Model Updating Method Based on Kriging Model for Structural Dynamics

Model updating in structural dynamics has attracted much attention in recent decades. And high computational cost is frequently encountered during model updating. Surrogate model has attracted considerable attention for saving computational cost in finite element model updating (FEMU). In this study, a model updating method using frequency response function (FRF) based on Kriging model is proposed. +e optimal excitation point is selected by using modal participation criterion. Initial sample points are chosen via design of experiment (DOE), and Kriging model is built using the corresponding acceleration frequency response functions.+en, Kriging model is improved via new sample points using mean square error (MSE) criterion and is used to replace the finite element model to participate in optimization. Cuckoo algorithm is used to obtain the updating parameters, where the objective function with the minimum frequency response deviation is constructed. And the proposed method is applied to a plane truss model FEMU, and the results are compared with those by the second-order response surface model (RSM) and the radial basis function model (RBF). +e analysis results showed that the proposed method has good accuracy and high computational efficiency; errors of updating parameters are less than 0.2%; damage identification is with high precision. After updating, the curves of real and imaginary parts of acceleration FRF are in good agreement with the real ones.


Introduction
e accurate finite element model (FEM) is the basis of reflecting structural dynamics characteristics and guiding the optimization design of the structure.However, due to the modelling error and other uncertain factors, there is always a certain discrepancy between the FEM and the accurate model.Difference always exists between the mathematical model and experimentally measured one.erefore, it is a general practice to update the theoretical model by using experimental measurements in order to obtain a more accurate model.Finite element model updating (FEMU) can be defined as the process of tuning a FEM to minimize the discrepancy between the measured and computed responses of the structure [1].FEMU has attracted significant attention from the structural engineering community because of its applications in structural dynamics [2].e current mainstream FEMU methods can be divided into two categories: FEMU based on modal parameters and FEMU based on frequency response function (FRF).e former needs modal parameter identification, which inevitably introduces the identification error, and requires that the measured degree of freedom (DOF) be consistent with the DOF of the analytical model [3].FEMU based on FRF does not need to carry out modal analysis and mode matching, so it is suitable for structures with relatively dense modal distribution [4].In addition, the FRF contains much information, which increases the available data and has a wider range of applications [5].
e results of structural model updating are important for engineers to accurately model similar types of structures in the future.Updated numerical models can be used in numerous engineering fields (e.g., structural vibration control, force identification, vehicle load identification of bridge systems, and structural damage detection) [6].One of the main limitations of the FEMU in the engineering application is the computational efficiency.In theory, the higher the density of the finite element mesh, the higher the accuracy of the structure response prediction and the more reliable the model updating result.However, a fine finite element mesh will result in a significant increase in the computational time of the individual theoretical model response.Methods for constrained nonlinear optimization are typically used to solve the inverse problem of FEMU considering an objective function describing the discrepancies between the predicted FEM response and the measured one.Complex relationship exists between the model parameters to be identified and the objective function.e objective function may include many local minima.Gradient-based methods are commonly used for FEMU [7].However, they might be trapped in local minima and their solution heavily depends on the starting point (i.e., initial guess of the model parameters).In order to avoid this issue, different global optimization algorithms have been used for FEMU [2].e computational costs of the above methods are high.Surrogate model can provide a cheaper but lowerfidelity solution.e surrogate model is gaining popularity as a way of fast developed approximation for timeconsuming simulations [8].
e surrogate model does not need the complicated calculation as the FEM does.It is an approximate model constructed directly using the relationship between inputs (design parameters) and outputs (structural responses).e commonly used surrogate models include response surface method (RSM), neural network (NN), radial basis function (RBF), support vector regression (SVR), Kriging model, and so on [9][10][11][12].Kriging model is an equivalent model based on Kriging interpolation.Unlike other surrogate models, Kriging model can not only give the prevaluation of unknown function but also get the error estimate of the prevaluation.Only a small number of samples can accurately describe the relationship between the inputs and outputs of the structure, which is widely used in the field of evaluation of machine tools [13], structural reliability analysis [14], and model updating [15].In the field of FEMU, Zhang and Guo [16] applied Kriging theory to FEM confirmation to predict the response of the structure.Jensen et al. [17] integrated an adaptive metamodel into a FEMU using dynamic response data.e updating technique is combined with an adaptive surrogate model based on Kriging interpolation of the measure-of-fit function.Yang et al. [18] presented a new surrogate model-based FEMU taking advantage of the measured FRFs.To improve the accuracy of the surrogate model, Han et al. [19] proposed a novel formulation of gradient-enhanced surrogate model, called weighted gradient-enhanced Kriging, in combination with the adjoint method to ameliorate the curse of dimensionality for highdimensional surrogate model and design optimizations.Leifsson et al. [20] presented a method for multiobjective design optimization with variable-fidelity models and response surface surrogates, and the underlying low-fidelity model is corrected using parameterized output space mapping.
However, most methods for constructing Kriging model need to take some samples and to build a Kriging model with certain precision based on the samples and their corresponding responses.en, the Kriging model takes place of the FEM to carry out iterative optimization.If the precision of the constructed Kriging model is too low, it may result in poor or even failed optimization results.In this case, samples must be reextracted and the Kriging model must be reconstructed.In addition, most of the methods have chosen the frequency as the response of the surrogate model.However, as an important part of model updating, FRF can provide more structural dynamic information.Unfortunately, FRF is rarely used as the response of surrogate model.erefore, in this paper, the Kriging model with good simulation effect is applied to the model updating based on FRF data.e initial Kriging model is constructed by the parameters to be modified and the acceleration frequency response functions.In order to improve its accuracy, new sample points are added with the mean square error (MSE) criterion, and the Kriging model is updated iteratively.en, the Kriging model which meets the accuracy requirement is applied to model updating of a truss structure.e Cuckoo algorithm with better global optimization ability is used to optimize the parameters.e numerical example verifies the effectiveness of the proposed method.e remainder of this paper is organized as follows: Section 2 gives an overview of the Kriging model.Section 3 presents the model updating based on the Kriging model, including the MSE criterion, the objective function, and the model updating process.Section 4 provides case studies.Finally, Section 5 provides the conclusions.

Kriging Model
Kriging model is considered as the best linear unbiased estimation to the real computer model.It is a semiparametric interpolation technique which estimates the unknown information at one point according to the known information [21].Nowadays, it has become a popular method for approximating deterministic computer model [22].
A Kriging model is a surrogate model based on a stochastic process.e model includes linear regression part and nonparametric part.And the latter can be considered as the realization of a stochastic process.For a given set of sample data (input) X � [x 1 , x 2 , . . ., x n ] T and the observed response (output) Y � [y 1 , y 2 , . . ., y n ] T , the expression of Kriging model is where β � [β 1 , β 2 , . . ., β p ] T is the coefficient vector to be estimated; f(x) � [f 1 (x), f 2 (x), . . ., f p (x)] T is the polynomial vector of the sample x ; p is the number of f l (x); n is the number of sample points; and z(x) represents error and is assumed to be a stochastic process that follows a normal distribution N(0, σ 2 ) with a zero mean and standard deviation σ generally and has nonzero covariance estimated by where x i and x j are two sample points; σ2 is the variance of z(x); and R(x i , x j ) is the correlation function characterizing 2 Shock and Vibration the correlation between any two sample points.Different R(x i , x j ) may induce different approximation accuracy of the constructed Kriging model.To estimate the stochastic process z(x), assuming that the true response surface of Kriging model is continuous, any two points will tend to have the same value as the distance approaches zero and it is the same for z(x) of two points.us, the correlation between z(x) of any two sample points can be expressed as a function of their spatial distance.e most widely used Gaussian correlation model is adapted [18]: where x k i and x k j are the kth components of the two sample points x i and x j , respectively; m denotes the number of design variables; and θ k is the correlation coefficient which controls the decay rate of correlation on different dimensions.And then the matrix of correlation functions can be expressed as When training the parameters in the Kriging model, the maximum likelihood estimation is usually used.e likelihood function of the sample point can then be expressed as where |R| is the determinant of R, which is a function of θ k ; Y is a column vector which is composed of the outputs (responses) of training samples; and F is a matrix of vector f(x) from each sample point.e least squares estimations of β and σ 2 can be expressed as Substituting equation ( 6) and equation (7) in equation ( 5) and ignoring the constant term, the logarithmic form of the maximum likelihood function can be expressed as Both  σ and R are the functions of θ k .And they can be obtained from the maximum likelihood estimation of the response.So, by solving the maximum value of the above expression, the value of θ k can be determined.
Stochastic distribution z(x) (in equation ( 1)) is sufficient to indicate the changing trend of the function response, so the regression part f T (x) is regarded as a constant term.en, the unique unknown parameter in the Kriging model is θ k , which can be obtained through optimization algorithm.
After the Kriging model is constructed, the next step is to predict the response value of the point to be measured.For any point x 0 to be measured, its response value  y(x 0 ) can be expressed as e prediction accuracy can be evaluated by estimating MSE of the predicted value, which can be estimated as where r T (x 0 ) is the row vector of correlation function between each sample point and the point x 0 to be measured:

Improvement of Kriging Model.
e accuracy of Kriging model directly affects the results of model updating.Constructing a high-precision Kriging model requires a large number of sample points.However, the number of sample points is the key to restrict the computational efficiency of Kriging model.e prediction of each new point needs to compute the correlation between the new sample point and all sample points.erefore, in order to ensure the computational efficiency, the number of sample points should not be too large.But small number of sample points may lead to low accuracy of the model.In order to make Kriging model have a good compromise between computational efficiency and accuracy, it is necessary to add new sample points in effective area according to certain rules (infill criteria).
At present, the infill criteria mainly include improved expectation criterion (EI criterion), improved probability criterion (PI criterion), confidence lower bound criterion (LCB criterion), and mean square error criterion (MSE criterion) [12].In this paper, the MSE criterion is selected to add new sample points for updating the Kriging model.e estimation of root mean square error (RMSE) provided by the Kriging model is used to guide the addition of new sample points, thereby improving the global accuracy of the model.Relative to the optimization of model updating, the optimization here is referred as suboptimization.In short, it is to use a reasonable optimization algorithm to solve the following optimization problem so as to get a new sample point: of R; r is a correlation vector, consisting of correlation functions between unknown points and all known sample points; F is a matrix of polynomial vector for each sample point; and F T is the transpose of F.
It is necessary to choose a suitable convergence criterion to terminate the surrogate infill process.Here, the maximum number of new samples is taken as the termination criterion.
e solution for the location of the maximum RMSE is finished by Cuckoo algorithm.
e main steps of updating the Kriging model can be described as follows: e RMSE distribution of one-dimensional Schwefel function is shown in Figure 1.As shown in Figure 1(a), the sample point corresponding to the maximum RMSE is around 0.8. Figure 1(b) is obtained after adding this sample point; it can be seen that the precision of the surrogate model is greatly improved.erefore, it is an effective method to update the Kriging model by using the optimization algorithm step by step to select the point with the largest RMSE as the new sample point.

Objective Function.
Essentially, model updating is an optimization problem.Minimizing the objective function, the modified values of design parameters are obtained by solving the optimization problem.In this paper, the objective function is established based on the differences of the acceleration frequency response function between the experimental model and Kriging model.
In the frequency domain, the general mathematical representation of n DOFs system is expressed as where M, K, and C are n × n matrices of mass, stiffness, and damping, respectively.X(ω) is n × 1 displacement response vector subjected to the applied force vector of F(ω).ω is the excitation frequency.e displacement vector X(ω) can be achieved by where H(ω) is a n × n FRF matrix defined as e objective function here is defined as where  y(ω i ) represents the prediction response of the updated Kriging model at ω i ; max F denotes the maximum observation frequency; and H(ω i ) represents the response of the structure at ω i .e peak value of the FRF is usually large at some frequency points, which will affect the prediction accuracy of the Kriging model, so the logarithmic (log (FRF)) is taken as the response.
Many optimization algorithms can be used to solve equation (17).e Cuckoo algorithm [23] has strong global search ability due to its unique Lévy performance.Moreover, it has the advantages of less parameters, simple operation, easy implementation, and good optimization performance.
erefore, the Cuckoo algorithm is used to solve the parameters to be modified.

Process of Model Updating.
As mentioned above, when the Kriging model is introduced into the structural FEMU, the Kriging model replaces the FEM to participate in optimization, and thus FRF model updating is fulfilled.
e main steps of model updating method based on Kriging model can be summarized as follows: Step 1: select the excitation point and measurement point and determinate the parameters to be modified.
Step 2: generate initial sample points of the modified parameters by design of experiment (DOE).
Step 3: run the finite element analysis program to compute the response of the sample points and construct the initial Kriging model.
Step 4: find the point with the maximum RMSE of the current Kriging model by the Cuckoo algorithm and add it into the set of sample points.
Step 5: compute the responses of the new sample points and update the initial Kriging model.
Step 6: check whether the maximum number of new samples is satisfied.If it is, go to Step 7. Otherwise, go back to Step 4 and continue adding new points.

Numerical Example
A plane truss structure (Figure 3) is taken as an example to verify the model updating method.e truss has 14 nodes and 25 DOFs.e elasticity modulus of each bar element is 210 GPa and the density is 7850 kg•m −3 .Select the first fourth modes as the modes of interest.

4
Shock and Vibration

Selection of Excitation Points and Measurement Points.
Modal participation criterion is used for selecting excitation points.Modal participation [24] is applied to evaluate the contribution of each DOF to the excitation modes of the structure.Modal participation can be expressed as where subscripts p and q denote the output DOFs and input DOFs, respectively; r represents number of modes; N 0 represents the number of DOFs; and A pqi represents the residue.
Assuming that the structure is proportional damping, the FRF can be expressed as where ϕ pi and ϕ qi represent the (i, p)th and (i, q)th elements the modal matrix, respectively; ω i denotes the undamped natural frequency at the ith order mode; and ξ i represents the damped coefficient at the ith order mode.A pqi in equation ( 7) can be expressed as e contribution of the qth DOF to the excitation of all modes of interest can be expressed as In single-input and single-output modal response test, structural DOF with maximum W q can be selected as the optimal excitation point.W qi and W q of each DOF of the structure are shown in Figure 4. e maximum W q is at the 20th DOF in the Y-direction of node 11.So, the 20th DOF is chosen as the best excitation point in the test.e DOFs in the Y-direction make the most important contribution to the first four modes of excitation. is is consistent with the characteristics of the plane truss structure.
ere are also many methods for selecting measurement points, but most of them are multipoint selection.According to the characteristics of the model, a DOF in the Y-direction should be selected as the sensor measurement point.e

Shock and Vibration
fourth DOF in the Y-direction at the 3rd node is chosen as the measurement point.

Test Design and Initial Sample Selection.
In this paper, the FEM shown in Figure 3 is used as the test model, and the corresponding FEM is obtained by deviating from the values of the parameters to be modified in the test model.Considering that the test data will be affected by noise in the actual data measurement, 5% random noise is added to the response of the test model.
e elasticity modulus and material density are chosen as parameters to be modified and their values are deviated.e elasticity modulus is increased by 10%.e material density is reduced by 10%.e FEM parameters are shown in Table 1.
Before constructing the surrogate model, selecting the sample points is the first step.For the global optimization problem, a better method is to select a set of sample points through DOE. e method proposed in this paper has no strict requirements for the number of initial sample points, which is different from the traditional Kriging model.Considering the efficiency of adding new sample points, the number of initial samples should not be too small.In other words, the optimization efficiency of the improved Kriging model based on our method is not obviously dependent on the number of initial sample points.Here, the Latin hypercube sampling (LHS) is used to sample the two parameters (elasticity modulus and material density).Finally, 40 samples are extracted.

Construction and Verification of Kriging Model.
On the basis of the 40 extracted samples and their corresponding acceleration FRFs, the initial Kriging model is constructed.
en, the new samples are added according to the MSE criterion introduced in Section 3, and the maximum number of additional sample points is 40.
en, the improved Kriging model is built.e FRFs at the fourth DOF of the FEM are predicted by using Kriging model and the improved Kriging model, respectively.e results are shown in Figure 5 (RMSE1 denotes RMSE of Kriging model; RMSE2 denotes RMSE of the improved Kriging model).Both the Kriging model and the improved Kriging model have good prediction accuracy.However, the improved Kriging model performs better in predicting the peak value of the curve, and   6.
e overlap between the sample response and the predicted response from the improved Kriging model is good, and the maximum RMSE is less than 0.05.e improved Kriging model has good approximation accuracy to the FEM.

Model Updating.
e improved Kriging model mentioned above is used to substitute the FEM to optimize iteratively.Assume that the test parameters are within the interval of finite element parameter values.
e Cuckoo algorithm is used to find the optimum iteratively.e number of nests is 40, and the maximum number of iterations is 100.In order to prove the stability of the algorithm, the iterative convergence curve is shown in Figure 7 when the algorithm runs 100 times.e Cuckoo algorithm is stable and converges before the number of iterations reaches 60.
e difference between the optimal value and the worst value is also very small.
In order to compare with the optimization effect of the updated Kriging model, the second-order RSM and the RBF are constructed based on the same sample points.e average values and average errors of the modified parameters are shown in Table 2. Except for the second-order RSM, all the other three methods obtain good results, but the parameter values by the improved Kriging model are more accurate, and the average error is the smallest.
All the algorithms are coded in Matlab 2014b.e operating system is Windows 10.Simulation hardware is a PC with 3.50 GHz Intel Pentium (R) G4560 and 12.00 GB memory, and all the algorithms are repeated 5 times.e total computation time of model updating using these surrogate models is counted, and the average time is shown in Table 3.
e RSM takes the least time.e RBF needs the longest time.And the Kriging model needs longer time than RSM, but its computational efficiency is much higher than RBF.
In summary, the Kriging model has the best accuracy and needs less time.e accuracy of RSM is too poor.e computation time of RBF is too long.RSM and RBF cannot be compromised between accuracy and time.
is shows that the proposed method can improve the computational efficiency while satisfying the accuracy and has little dependence on the number of sample points.It is noted that there are some errors between the Kriging model and the improved Kriging model, which are not caused by the insufficient optimization ability of the optimization algorithm, but by the prediction error of the Kriging model itself.So, only if the Kriging model is accurate enough, it can be used to modify the structure model to reduce the computational cost and get the accurate updating results.
e FRFs predicted by the Kriging model, the updated Kriging model, RBF, and RSM are shown in Figure 8.
Comparison of updated Kriging model with Kriging model, RSM, and RBF is shown in Figures 8(a e FEM is updated with the average modified values obtained from the updated Kriging model.By comparing the real and imaginary FRF curves before and after modification at the measurement points, the updating effect of the proposed method is further verified. e comparison curve between real and imaginary parts of FRF before and after updating is shown in Figure 9. e FRFs of the updated FEM (whether real or imaginary part), can coincide well with the FRFs of the test model.

Damage Identification.
In order to further verify the proposed method, the above mentioned truss structure elements 8, 12, 15, 17, 21, and 25 are assumed to be damage elements.eir stiffness of the elements is reduced by 10%, 10%, 20%, 20%, 30%, and 30%, respectively.e elasticity modulus of each damaged unit is identified.e number of samples is set to 500.e maximum number of new additional sample points is 100.
e maximum number of iterations of the Cuckoo algorithm is 300.
e damage identification errors by the Kriging model, the updated Kriging model, the second-order RSM, and RBF under  4. e effects of different measurement points on damage identification accuracy are different.When sensors are evenly arranged at all Y-direction DOFs, 4 surrogate models have better identification accuracy.e identification error of the updated Kriging model is less than 0.5%.For a single measurement point, except for 14-Y measurement point, the identification accuracy of the four methods for unit 12 at other measurement points is poor; at 5-Y measurement point, the surrogate models have larger identification error.
e identification accuracy of the updated Kriging model is better than that of the other surrogate models on the whole.For the identification of elasticity modulus of the truss damage elements, the computation time of the Kriging model is about twice that of RBF and 0.04 times that of RSM.
FEMU based on Kriging model can make a good compromise between computational accuracy and time.e accuracy of damage identification is not only affected by the surrogate models but also affected by the selection of measurement points.

Conclusions
In this paper, the updated Kriging model is used to modify the parameters and to identify the damage of the structure.
e conclusions are as follows:    Shock and Vibration 9 (1) For modal response test, the selection of excitation point and measurement point is very important.For single-input and single-output test, modal participation criterion can be used to select the excitation point, and the response point can be selected according to the structural characteristics.(2) e updated Kriging model obtained by using the MSE criterion has better accuracy than the Kriging model, and the accuracy of the model is more stable.(3) Compared with the second-order RSM and RBF, the updated Kriging model can improve the computational efficiency while satisfying the computational accuracy and has less dependence on the number of sample points.(4) e parameter errors of the updated model are less than 0.2%.When there are enough measurement points, the parameter errors of damage identification are less than 0.5%.(5) When the optimization algorithm has enough optimization performance, the parameter modification errors often depend on the accuracy of the surrogate models, which are independent of the optimization algorithm, and when the surrogate models are accurate enough, the parameter modification errors will tend to zero.
Further research is needed on the selection of multiple measurement points for model updating.In the following research, modal response test will be carried out.

( 1 )
Construct the initial Kriging model according to the theory in Section 2. (2) Compute the maximum RMSE value of the current Kriging model by the Cuckoo algorithm and add it into the set of sample points.(3) Run the finite element analysis program to compute the responses of the new sample points and update the initial Kriging model until the set termination condition is reached.
) as the target; the updated Kriging model is used to replace the FEM to participate in optimization, and the modified values of parameters are obtained by the Cuckoo algorithm.e flow chart of model updating is shown in Figure 2.
) and 8(b), respectively (RMSE1 denotes RMSE of Kriging model; RMSE2 denotes RMSE of the updated Kriging model; RMSE3 denotes RMSE of RSM; RMSE4 denotes RMSE of RBF).e prediction value of the updated Kriging model is more accurate than others.e form of FRFs of the test model and the FEM does not change, and the peak values are very close, only causing movement along the frequency axis.

Figure 8 :
Figure 8: Log(FRF) of experimental model.(a) Comparison of the updated Kriging model with the Kriging model.(b) Comparison of the updated Kriging model with RBF and RSM.

Table 1 :
Parameters of test model and FEM.

Table 2 :
Parameters and errors of model updating.

Table 4 :
Identification error of elasticity modulus of each damage unit.Note.K represents the Kriging model; U-K represents the updated Kriging model; 3-Y represents the measurement points arranged in the Y direction at the 3rd node.