# A Copula-Based and Monte Carlo Sampling Approach for Structural Dynamics Model Updating with Interval Uncertainty

Authors: Xueqian Chen; Zhanpeng Shen; Xin’en Liu
Journal: Shock and Vibration
 (2018)

---

## Abstract

As the uncertainty is widely existent in the engineering structure, it is necessary to study the finite element (FE) modeling and updating in consideration of the uncertainty. A FE model updating approach in structural dynamics with interval uncertain parameters is proposed in this work. Firstly, the mathematical relationship between the updating parameters and the output interesting qualities is created based on the copula approach and the vast samples of inputs and outputs are obtained by the Monte Carlo (MC) sampling technology according to the copula model. Secondly, the samples of updating parameters are rechosen by combining the copula model and the experiment intervals of the interesting qualities. Next, 95% confidence intervals of updating parameters are calculated by the nonparameter kernel density estimation (KDE) approach, which is regarded as the intervals of updating parameters. Lastly, the proposed approach is validated in a two degree-of-freedom mass-spring system, simple plates, and the transport mirror system. The updating results evidently demonstrate the feasibility and reliability of this approach.


---

## Body

## 1. Introduction

Finite element (FE) models that numerically solve various engineering problems can aid virtual prototyping, reduce product development cycle, and cut down the cost of performing the physical tests. However, the reliability of the simulation results by finite element modeling is not always guaranteed since FE models are the approximations of real world phenomena based on various assumptions. These assumptions may detract from the quality and accuracy of simulation results. In order to improve the accuracy of FE simulating to serve the structural design better, the FE model updating techniques are needed to develop. In the past few decades, various kinds of FE model updating approaches have been widely investigated based on the actually observed behaviors of the system. Additionally, experimental modal and vibration data are often used in FE model updating in the field of structural dynamics [1–6].In most model updating approaches, the simulations are usually deterministic where each of the updating parameters is considered to have one “true” value and the purpose of the updating procedure is to provide a deterministic estimation. In reality, there are always uncertainties in nominally identical structures, such as the structural parameter uncertainty (physical material properties, geometric parameters), the assembly joints uncertainty, and the experiment uncertainty (measurement noise, modal identification techniques, etc.). As a result, the FE model updating approaches with uncertainty have received great attentions recently. Studies have shown that the simulation results are more reliable when the uncertainties are taken into account [7], suggesting that it is necessary to consider the uncertainties during modeling and simulating [8].In FE model updating approaches with uncertainty, the updated parameters are no longer deterministic and are described as random variables. Usually, the FE model updating approaches with uncertainty can be classified into two major categories: probabilistic and nonprobabilistic approaches. In the earlier works, a probabilistic approach proposed incorporated the measurement noise into model updating [9]. Subsequently, Bayesian statistical frameworks were adopted to estimate the posterior probabilities of uncertain parameters [10–12]. However, high computational costs due to a large amount of samples required for a satisfactory estimation greatly restrain the applications of Bayesian updating approaches. As a result, surrogate models such as the Gaussian process model with the perturbation approaches and sensitivity analysis approaches have been employed in stochastic model updating to improve the efficiency [13–16]. Though, the surrogate model approaches own the superiority of computational efficiency over Monte Carlo (MC) based methods. Nevertheless, the prerequisite of small uncertainties, together with the Gaussian distribution assumption, also limits the applications to complex problems. Moreover, perturbation based predictions are sensitive to the initial estimates of parameters. Recently, an approach with the response surface models and MC simulation has been developed, which decomposed a stochastic updating process into a series of deterministic ones [17]. On the other hand, the accuracy of the probabilistic approaches depends on the estimation of the probability distribution characteristics of the structural parameters and the responses. The establishment of an accurate probability distribution function (PDF) needs lots of experiment data in the probabilistic approaches, which greatly limits its application in engineering.In nonprobabilistic approaches, the interval approach has been intensively investigated. By comparison, the experiment samples are not strictly needed in the FE model updating with interval analysis as was proposed. In the field of interval model updating (IMU), the inclusion theorem was employed to establish an interval inverse problem. And the convergence was achieved when measured responses fall into numerically predicted intervals [18–22]. Considering the easy implementation, IMU problems are usually solved within a deterministic framework where the upper and lower bounds of parameters are sought separately. For example, an IMU problem was decomposed into two deterministic constrained optimization processes where the midpoints and interval radii of parameters were separately estimated [19]. Alternatively, the vertex solution theorem is effective and cost-efficient for IMU due to its easy implementation [20], particularly in the solution of Eigen value problems [21]. But the vertex solution was valid only for particular parameterization of an FE model without the involvement of eigenvectors, which highly limits its further applications. Due to this drawback, global optimization algorithms were taken into account for more general solutions. Surrogate models such as the Kriging predictor and interval response surface were used to improve the efficiency of gradient computation and facilitate the convergence [23, 24]. So far most of IMU problems are solved within a deterministic framework since direct interval arithmetic operations are difficult to implement during inverse solutions. Therefore the upper and lower bounds of parameters should be sought separately through a deterministic inverse procedure. Additionally, global optimization of interval variables is difficult to realize due to the fact that the interval arithmetic is quite different with the traditional mathematical arithmetic.Though several probabilistic and interval model updating approaches have been developed in the past years, most of them are still complicated for implementation. Additionally, these approaches with uncertainty suffer from the challenges such as ill-condition, nonuniqueness and local optimal solution, etc. To overcome such inconvenience, an IMU approach is developed in this work based on the copula model and MC sampling. In the proposed approach, the copula model between the updating parameters and the interesting qualities is constructed firstly. Then a large amount of samples is obtained according to the copula model, and the samples are rechosen based on experiment intervals of interesting qualities. Lastly, the updating intervals of parameters are obtained by estimating on the rechosen samples with kernel density estimation (KDE). The remainder of the paper is organized as follows. In Section2, the copula-based FE model updating approach and procedure with interval uncertainty are presented. In Section 3, three examples are provided to validate the accuracy and reliability of the proposed approach. Conclusions are presented in Section 4.


## 2. Identification of Interval Parameters

The FE model updating problems are classic inverse problems in structural mechanics where the standard “forward” relationship between input and output variables of a model is inverted. The key in solving a FE model updating problem is to construct the mathematical relationship between the updating parameters and the output interesting qualities. The copula function is one of the most effective mathematical tools to determine this relationship, which expediently characterizes the correlation between the marginal functions of multivariables and the joint distribution function.


### 2.1. Brief Introduction of the Copula Function

A copula function is a general function in statistics to formulate a multivariate distribution with various statistical dependence patterns, which was presented by Sklar in 1959 [25]. Formally, a copula is a joint distribution function of standard uniform random variables. According to the Sklar’s theorem, there exists a two-dimensional copula C such that variables x1 and x2 in a real random space.(1)Fx1,x2=CF1x1,F2x2whereF(x1,x2) is a two-dimensional distribution function with marginal functions F1(x1) and F2(x2) and CF1x1,F2x2 is the copula cumulative distribution function (CDF).Equation (1) can be spread for m-dimensional variables easily, that is,(2)Fx1,x2,…,xm=CF1x1,F2x2,…,FmxmConsequently, the m-dimensional PDF is as follows:(3)fx1,x2,…,xm=∂mFx1,x2,…,xm∂x1∂x2…∂xm=∂mCF1x1,F2x2,…,Fmxm∂x1∂x2…∂xm=cF1x1,F2x2,…,Fmxm×∏i=1mfixiwherecF1x1,F2x2,…,Fmxm is the copula PDF, fx1,x2,…,xm is the united PDF for m-dimensional random variables, and fixi is the PDF of the ith random variable.At present, the general copula function types include the Gaussian copula function, t-copula function, and Archimedean copula function [26, 27]. Among them, the Gaussian copula is widely utilized because most of the parameters in the engineering satisfy the normal distribution. In the study, the Gaussian copula is adopted for the FE model. Specifically, the Gaussian copula function is constructed by multidimensional Gaussian distribution and the linear correlation parameters, and its distribution function is as follows:(4)CGau1,u2,…,um;ρ=ΦρΦ-1u1,Φ-1u2,…,Φ-1umwhereΦρ is the distribution function of the standard normal function for d-dimensional with the correlation matrix ρ, Φ-1 is the inverse function of the distribution function of the standard normal function, and ui=Fixi, i=1,2,…,m.


### 2.2. Copula-Based Approach for Model Updating with Interval Uncertainty

Firstly, the original design spaces of updating parametersxj, j=1,2,⋯,m, are assumed, and a few samples are obtained by design of experiment (DOE) approach and subsequent deterministic FE analysis on samples according to DOE. Then, the samples of the output interesting qualities yk, k=1,2,⋯,n, are obtained from the FE analysis results. Secondly, the copula model is constructed according to the samples of updating parameters and output interesting qualities, and resampling is performed to get large samples with number N for updating parameters and interesting response qualities based on the copula model. The samples falling into the experiment data space are considered to characterize the input-output relationship of the physical structure believably, and unuseful samples are needed to remove. Next, the samples of updating parameters xj are rechosen according to the experiment intervals of interesting qualities, as follows:(5)xj=xj=fy1,⋯,yn∣yk∈R,yk-≤yk≤yk+,j=1,2,⋯,m,k=1,2,⋯,nwhereR is the real space and yk+ and yk- are the upper and lower bounds of the kth output interesting quality which can be obtained from the experiment results.In practical model updating, the measured data are only a few samples in general. Reasonable interval estimation on experiment data is the precondition to obtain the reliable updated FE model. However, the KDE allows for the capture of the observed distributional structure for the random variables, without having to assume a particular parametric distribution form.Following [28], the kernel density estimator for variablex has the form(6)f^h=1nh∑i=1nKx-Xihwheren is the number of observations used to construct the estimate, K(·) is a kernel function, Xi is the ith observation, and h is the window width, or bandwidth. A typical choice for the kernel K(·) is the standard normal density and is implemented here. The choice of the window width h is usually based on the optimization of some scoring function. A least-square cross-validation score function is adopted for this work [28].The empirical CDF and 95% confidence interval (CI)[ye_95-ye_95+]of the random variable x can be obtained by KDE in Matlab that is regarded as the interval of the random variable in this work.Considering the fact that the estimation on the original intervals of the updating parameters may be inaccurate, the reliable intervals are not identified through one copula-based FE model updating procedure. In order to overcome this problem, the idea of the adaptive response surface technique is adopted for this work [29]. That is, in order to get the final updating results, multiloop on the copula-based model updating procedures may be performed.In the FE model updating procedure, the convergent criterion is that the difference of the updating parameter intervals between the(i+1)th iteration step and the ith iteration step is less than a small value, or the difference of the output interesting quality intervals between the ith iteration step results and the experiment results is less than a critical value.In order to improve the efficiency and the validity of model updating, the renewal strategy of updating parameters is as follows in each iteration step. The current intervals of updating parameters are updated according to the results of the previous iteration step, and the interval medians of the previous step are regarded as the current interval medians, and about 80% of the interval width of the previous step is regarded as the current interval width. Also, the Latin Hypercube Sample (LHS) method is suggested in the DOE, and the number of samples is not less than ten.The copula-based model updating procedure is repeated until the convergent criterion is satisfied. The flow chart for the copula-based FE model updating is outlined in Figure1.Figure 1
Flow chart of the copula-based FE model updating.


## 2.1. Brief Introduction of the Copula Function

A copula function is a general function in statistics to formulate a multivariate distribution with various statistical dependence patterns, which was presented by Sklar in 1959 [25]. Formally, a copula is a joint distribution function of standard uniform random variables. According to the Sklar’s theorem, there exists a two-dimensional copula C such that variables x1 and x2 in a real random space.(1)Fx1,x2=CF1x1,F2x2whereF(x1,x2) is a two-dimensional distribution function with marginal functions F1(x1) and F2(x2) and CF1x1,F2x2 is the copula cumulative distribution function (CDF).Equation (1) can be spread for m-dimensional variables easily, that is,(2)Fx1,x2,…,xm=CF1x1,F2x2,…,FmxmConsequently, the m-dimensional PDF is as follows:(3)fx1,x2,…,xm=∂mFx1,x2,…,xm∂x1∂x2…∂xm=∂mCF1x1,F2x2,…,Fmxm∂x1∂x2…∂xm=cF1x1,F2x2,…,Fmxm×∏i=1mfixiwherecF1x1,F2x2,…,Fmxm is the copula PDF, fx1,x2,…,xm is the united PDF for m-dimensional random variables, and fixi is the PDF of the ith random variable.At present, the general copula function types include the Gaussian copula function, t-copula function, and Archimedean copula function [26, 27]. Among them, the Gaussian copula is widely utilized because most of the parameters in the engineering satisfy the normal distribution. In the study, the Gaussian copula is adopted for the FE model. Specifically, the Gaussian copula function is constructed by multidimensional Gaussian distribution and the linear correlation parameters, and its distribution function is as follows:(4)CGau1,u2,…,um;ρ=ΦρΦ-1u1,Φ-1u2,…,Φ-1umwhereΦρ is the distribution function of the standard normal function for d-dimensional with the correlation matrix ρ, Φ-1 is the inverse function of the distribution function of the standard normal function, and ui=Fixi, i=1,2,…,m.


## 2.2. Copula-Based Approach for Model Updating with Interval Uncertainty

Firstly, the original design spaces of updating parametersxj, j=1,2,⋯,m, are assumed, and a few samples are obtained by design of experiment (DOE) approach and subsequent deterministic FE analysis on samples according to DOE. Then, the samples of the output interesting qualities yk, k=1,2,⋯,n, are obtained from the FE analysis results. Secondly, the copula model is constructed according to the samples of updating parameters and output interesting qualities, and resampling is performed to get large samples with number N for updating parameters and interesting response qualities based on the copula model. The samples falling into the experiment data space are considered to characterize the input-output relationship of the physical structure believably, and unuseful samples are needed to remove. Next, the samples of updating parameters xj are rechosen according to the experiment intervals of interesting qualities, as follows:(5)xj=xj=fy1,⋯,yn∣yk∈R,yk-≤yk≤yk+,j=1,2,⋯,m,k=1,2,⋯,nwhereR is the real space and yk+ and yk- are the upper and lower bounds of the kth output interesting quality which can be obtained from the experiment results.In practical model updating, the measured data are only a few samples in general. Reasonable interval estimation on experiment data is the precondition to obtain the reliable updated FE model. However, the KDE allows for the capture of the observed distributional structure for the random variables, without having to assume a particular parametric distribution form.Following [28], the kernel density estimator for variablex has the form(6)f^h=1nh∑i=1nKx-Xihwheren is the number of observations used to construct the estimate, K(·) is a kernel function, Xi is the ith observation, and h is the window width, or bandwidth. A typical choice for the kernel K(·) is the standard normal density and is implemented here. The choice of the window width h is usually based on the optimization of some scoring function. A least-square cross-validation score function is adopted for this work [28].The empirical CDF and 95% confidence interval (CI)[ye_95-ye_95+]of the random variable x can be obtained by KDE in Matlab that is regarded as the interval of the random variable in this work.Considering the fact that the estimation on the original intervals of the updating parameters may be inaccurate, the reliable intervals are not identified through one copula-based FE model updating procedure. In order to overcome this problem, the idea of the adaptive response surface technique is adopted for this work [29]. That is, in order to get the final updating results, multiloop on the copula-based model updating procedures may be performed.In the FE model updating procedure, the convergent criterion is that the difference of the updating parameter intervals between the(i+1)th iteration step and the ith iteration step is less than a small value, or the difference of the output interesting quality intervals between the ith iteration step results and the experiment results is less than a critical value.In order to improve the efficiency and the validity of model updating, the renewal strategy of updating parameters is as follows in each iteration step. The current intervals of updating parameters are updated according to the results of the previous iteration step, and the interval medians of the previous step are regarded as the current interval medians, and about 80% of the interval width of the previous step is regarded as the current interval width. Also, the Latin Hypercube Sample (LHS) method is suggested in the DOE, and the number of samples is not less than ten.The copula-based model updating procedure is repeated until the convergent criterion is satisfied. The flow chart for the copula-based FE model updating is outlined in Figure1.Figure 1
Flow chart of the copula-based FE model updating.


## 3. Case Studies

### 3.1. Example 1: A Two-Degree-of-Freedom Mass-Spring System

A two degree-of-freedom mass-spring system is shown in Figure2. The deterministic parameters in the system are m1=m2=1 kg and k1=1 N/m. The uncertain interval parameters are k2=[0.81.2] N/m and k3=[0.91.1]  N/m.Figure 2
A two degree-of-freedom mass-spring system.For simplicity, it is assumed that the uncertain parameters are uniformly distributed. To create such kind of uncertainty, the LHS method is used to generate twenty experiment samples. Afterwards, the experiment results of the first two natural frequencies are obtained according to the samples. Considering the effect of small sample, the 95% CI of experiment results are estimated by KDE and the estimated results are[0.96461.0338] rad.s−1 and [1.57891.8875] rad.s−1, respectively, which are regarded as the intervals of experimental results.Assume that the original intervals ofk2 and k3 are the same as [1.22.2] N/m, and the first two natural frequencies are regarded as the output interesting qualities. The intervals of k2 and k3 are identified according to the copula-based model updating flow with interval uncertainty in Figure 1. The updating results are convergent after three iteration steps.Figure3 is the scatter map between the updating parameters and the output interesting qualities when the updating results are convergent, which shows that there is strong correlation between k3 and f1and between k2 and f2 and weak correlation between k2 and f1 and between k3 and f2.Figure 3
Scatter map between the updating parameters and the output quantities of mass-spring system.
(a)
k 2 VS f1 (b)
k 3 VS f1 (c)
k 2 VS f2 (d)
k 3 VS f2Then, the updated intervals ofk2 and k3 are obtained by KDE for the resample in Figure 3, which are [0.7751.230] N/m and [0.8821.133] N/m, respectively. The comparison between the original uncertain interval and the updated interval of updating parameters is shown in Figure 4, which shows that the updated interval matches the real interval better. Because the effect of small samples on experiment data is considered in the model updating procedure, the updated interval is bigger than the real interval and the result is reasonable.Figure 4
Patch plots for the updating parameters of mass-spring system.In order to validate the updating results on the interval uncertain parameters, the copula models are reconstructed according to the original and the updated intervals ofk2 and k3. As a result, 5000 samples of k2, k3, and the first two frequencies of the system are resampled by the copula models constructed just now. The scatter map between the simulating and the experimental results is shown in Figure 5, which indicates that the frequencies of the updated model are agreement with the experimental results better.Figure 5
Scatter plots for the first two natural frequencies of mass-spring system.


### 3.2. Example 2: Interval Model Updating in Simple Plates

Impact hammer modal testing with free-free boundary conditions was conducted on thirty-three nominally identical steel plates in [30]. The nominal geometric dimensions of the plates are 564 mm(length) × 110 mm(width) × 1.45 mm(thickness). And the nominal material properties are Young’s modulus of 210 GPa, the shear modulus of 83 GPa, and the mass density of 7860 kg/m3. The statistical properties of the first five measured natural frequencies of the plates are given in Table 1.Table 1
Statistical properties of the measured frequencies of the plates.
f 1 (Hz) f 2 (Hz) f 3 (Hz) f 4 (Hz) f 5 (Hz) mean 24.12 66.92 77.65 131.97 158.80 interval [23.94 24.4] [66.5 67.53] [76.24 78.91] [131.31 133.03] [156.31 160.94] Std. 0.11 0.25 0.57 0.42 0.97The FE model of the plate is created by the SHELL181 element in ANSYS, which has 300 shell elements. The boundary is free-free in the model. The material properties of the original FE model are the nominal values. The FE model and the first five mode shapes of simulation are shown in Figure6.Figure 6
The FE model and mode shapes of steel plates: (a) FE model; (b) mode 1; (c) mode 2; (d) mode 3; (e) mode 4; (f) mode 5.
(a) (b) (c) (d) (e) (f)After investigation, the uncertainty of rectangular plates can be characterized by Young’s modulusE, the shear modulus G, and the thick T in the FE model. The model updating is to identify the interval of these three uncertain parameters E, G, and T by the experiment results.Firstly, the intervals of the first five modal frequencies are estimated by KDE according to the experiment data in Ref. [30]. Secondly, it is assumed that the original intervals ofE,G, andT are [195220] GPa, [7887] GPa, and [1.301.60] mm, respectively, and intervals ofE,G, andT are identified according to the copula-based model updating flow with interval uncertainty in Figure 1. There are 20 experimental design data in each iteration step, i.e., 20 determined FE simulations in each iteration step. The model updating of the plate is convergent after four iteration steps with the convergent indices ε1=0.002 and ε2=0.002.After updating, the estimated intervals of the three parameters wereE=[203.21205.98] GPa, G=[82.7585.20] GPa, and T=[1.441.46] mm, respectively. The comparison between the original uncertain interval and the updated interval of the three parameters is shown in Figure 7, which shows that the updated intervals of these three parameters are much less than their original uncertain interval.Figure 7
Patch plots for the updating parameters of steel plates.
(a)
E vs G (b)
E vs TIn order to validate the updating results on the interval uncertain parameters, the copula models are reconstructed on the original and the updated intervals ofE, G, and T. The 5000 samples of E, G, T, and the first five natural frequencies of plates are obtained by MC sampling. The scatter plots for the simulation and the experimental results the first five natural frequencies of plates are shown in Figure 8. The comparisons on the natural frequencies between the simulation and the experimental results are listed in Table 2. It can be seen from Table 2 and Figure 8 that the frequencies of the updated model are in better agreement with the experimental results and the mean errors of frequencies decrease from the initial [12.5112.62]% to [0.720.24]%. Thus, this experimental validation has also confirmed the feasibility of the proposed IMU approach in estimating interval parameters of steel plates.Table 2
Comparison on the natural frequencies of plates between the simulation and the experimental results.
Mode Experimental interval/Hz Original interval/Hz Error /% Updated interval/Hz Error /% 1 [23.94 24.4] [21.06 27.55] [-12.05 12.92] [23.71 24.44] [-0.95 0.15] 2 [66.5 67.53] [58.61 76.53] [-11.86 13.32] [66.84 67.83] [-0.99 0.44] 3 [76.24 78.91] [65.81 87.80] [-13.68 11.26] [75.93 78.67] [-0.41 -0.30] 4 [131.31 133.03] [116.34 151.62] [-11.40 13.98] [130.37 134.26] [-0.71 0.92] 5 [156.31 160.94] [135.12 179.67] [-11.59 8.55] [155.51 160.95] [-0.51 0.01] mean   [12.51 12.62]  [0.72 0.24]Figure 8
Scatter plots for the first five natural frequencies of steel plates.
(a)
f1 vs f2 vs f3 (b)
f3 vs f4 vs f5


### 3.3. Example 3: Interval Model Updating in Transport Mirror System

The ShenGuangIII (SGIII) facility is designed for inertial confinement fusion (ICF) high energy experiments with 48 laser beams exactly transported and oriented to target. Figure9 is the view of the beam transport system in SGIII facility target area. There are 276 transport mirror systems in the facility, and the dynamic response under ambient vibration is a key factor to affect the stability of the SGIII facility [31]. A classic transport mirror system is shown in Figure 10, which consists of mirror component and support frame.Figure 9
A view of the beam transport system in SGIII facility target area, which is reproduced from Chen X J et al. (2014).Figure 10
Sketch map of the transport mirror system.The material of the support frame is steel with nominal Young’s modulus 200 Gpa, nominal Poisson’s ratio 0.3, and the nominal density 7850 kg/m3. The material of the mirror is K9 class with nominal Young’s modulus 80 Gpa, nominal Poisson’s ratio 0.21, and the nominal density 2510 kg/m3. The transport mirror system is about 1.32 m in height.The uncertainty exists in the transport mirror system because of the installation fluctuation and the welding technology fluctuation about the support frame, which induces the uncertainty of the vibration characteristics. In order to study how these uncertain factors affect the natural frequencies of the transport mirror system, impact hammer modal testing with fixed boundary conditions was conducted on ten nominally identical transport mirror systems. The experiment mount of the transport mirror system is shown in Figure11. Three natural frequencies in the interesting frequency range were found to significantly influence the response of the transport mirror system under the work condition. Then, the first three natural frequencies of the transport mirror system should be updated before calculating the response. The corresponding natural frequencies were obtained by the modal experiments, as listed in Table 3. It is observed that the frequency variations become more obvious with the increase of the mode order.Table 3
Experiment results of the first three natural frequencies of the transport mirror system.
Case 1 2 3 4 5 6 7 8 9 10 f 1/Hz 21.573 21.642 21.675 21.472 21.583 21.683 21.657 21.797 21.644 21.462 f 2/Hz 22.213 22.223 22.31 22.087 22.168 22.293 22.283 22.407 22.232 22.121 f 3/Hz 41.328 41.432 41.58 41.203 41.416 41.56 41.557 41.807 41.497 41.237Figure 11
Experiment setup of the transport mirror system.Generally, the bolts are ignored in the structural dynamic analysis model generally, and the FE model of the transport mirror system is established by the SOLID185 and SHELL181 element in ANSYS, which is shown in Figure12. A fully fixed bottom of the support frame is used as the boundary condition. The material properties of the original FE model are the nominal values. The first three simulation mode shapes are shown in Figure 13.Figure 12
FE model of the transport mirror system.Figure 13
Mode shapes of the transport mirror system.
(a)
Mode 1 (b)
Mode 2 (c)
Mode 3Theoretical analysis finds that the installation uncertainty of the transport mirror system can be characterized by Young’s modulusE1 of the bottom part of the support frame, and the weld uncertainty of the support frame can be characterized by Young’s modulus E2 of the upper part of the support frame. The goal of FE model updating is to identify the interval of E1 and E2.Firstly, the intervals of the first three modal frequencies are estimated by KDE according to the experiment data in Table3, which are [21.4121.83] Hz, [22.0322.45] Hz, and [41.0941.88] Hz, respectively.Secondly, assuming that the original intervals ofE1 andE2 are [100150] GPa and [160200] GPa, respectively, the intervals ofE1 andE2 are identified according to the copula-based model updating flow with interval uncertainty in Figure 1. There are 12 experimental design data in each iteration step, i.e., 12 determined FE simulations in each iteration step. The model updating of the transport mirror system is convergent after four iteration steps with the convergent indices ε1=0.005 and ε2=0.005.After updating, the estimated intervals of the three parameters wereE1 = [104.46114.65] GPa andE2 = [104.46114.65] GPa, respectively. The comparison between the original uncertain interval and the updated interval of the updating parameters is shown in Figure 14, indicating that the updated intervals of the updating parameters are much less than their original uncertain interval.Figure 14
Patch plots for the updating parameters of transport mirror system.In order to validate the updating results on the interval uncertain parameters, the copula models are reconstructed on the original and the updated intervals ofE1 andE2, Next, the 5000 samples ofE1,E2, and the first three natural frequencies of the transport mirror system are obtained by MC sampling. The scatter plots for the simulation and the experimental results of the first three natural frequencies of the transport mirror system are shown in Figure 15. The comparisons on the natural frequencies between the simulation and the experimental results are listed in Table 4. It can be seen from Table 4 and Figure 15 that the frequencies of the updated model agree with the experimental results better and the mean errors of frequencies decrease from the initial [2.4910.83]% to [0.470.40]%. As a result, this experimental validation has also proved the feasibility of the proposed IMU approach in estimating the interval parameters of the transport mirror system.Table 4
Comparison on the natural frequencies of the transport mirror system between the simulation and the experimental results.
Mode Experimental interval/Hz Original interval/Hz Error /% Updated interval/Hz Error /% 1 [21.46 21.80] [22.06 24.24] [2.77 11.20] [21.38 22.91] [-0.36 0.52] 2 [22.09 22.41] [22.68 24.94] [2.69 11.31] [22.00 22.54] [-0.38 0.59] 3 [41.20 41.81] [42.03 45.98] [2.01 9.97] [40.92 41.85] [-0.68 0.09] mean   [2.49 10.83]  [0.47 0.40]Figure 15
Scatter plots for the first three natural frequencies of transport mirror system.


## 3.1. Example 1: A Two-Degree-of-Freedom Mass-Spring System

A two degree-of-freedom mass-spring system is shown in Figure2. The deterministic parameters in the system are m1=m2=1 kg and k1=1 N/m. The uncertain interval parameters are k2=[0.81.2] N/m and k3=[0.91.1]  N/m.Figure 2
A two degree-of-freedom mass-spring system.For simplicity, it is assumed that the uncertain parameters are uniformly distributed. To create such kind of uncertainty, the LHS method is used to generate twenty experiment samples. Afterwards, the experiment results of the first two natural frequencies are obtained according to the samples. Considering the effect of small sample, the 95% CI of experiment results are estimated by KDE and the estimated results are[0.96461.0338] rad.s−1 and [1.57891.8875] rad.s−1, respectively, which are regarded as the intervals of experimental results.Assume that the original intervals ofk2 and k3 are the same as [1.22.2] N/m, and the first two natural frequencies are regarded as the output interesting qualities. The intervals of k2 and k3 are identified according to the copula-based model updating flow with interval uncertainty in Figure 1. The updating results are convergent after three iteration steps.Figure3 is the scatter map between the updating parameters and the output interesting qualities when the updating results are convergent, which shows that there is strong correlation between k3 and f1and between k2 and f2 and weak correlation between k2 and f1 and between k3 and f2.Figure 3
Scatter map between the updating parameters and the output quantities of mass-spring system.
(a)
k 2 VS f1 (b)
k 3 VS f1 (c)
k 2 VS f2 (d)
k 3 VS f2Then, the updated intervals ofk2 and k3 are obtained by KDE for the resample in Figure 3, which are [0.7751.230] N/m and [0.8821.133] N/m, respectively. The comparison between the original uncertain interval and the updated interval of updating parameters is shown in Figure 4, which shows that the updated interval matches the real interval better. Because the effect of small samples on experiment data is considered in the model updating procedure, the updated interval is bigger than the real interval and the result is reasonable.Figure 4
Patch plots for the updating parameters of mass-spring system.In order to validate the updating results on the interval uncertain parameters, the copula models are reconstructed according to the original and the updated intervals ofk2 and k3. As a result, 5000 samples of k2, k3, and the first two frequencies of the system are resampled by the copula models constructed just now. The scatter map between the simulating and the experimental results is shown in Figure 5, which indicates that the frequencies of the updated model are agreement with the experimental results better.Figure 5
Scatter plots for the first two natural frequencies of mass-spring system.


## 3.2. Example 2: Interval Model Updating in Simple Plates

Impact hammer modal testing with free-free boundary conditions was conducted on thirty-three nominally identical steel plates in [30]. The nominal geometric dimensions of the plates are 564 mm(length) × 110 mm(width) × 1.45 mm(thickness). And the nominal material properties are Young’s modulus of 210 GPa, the shear modulus of 83 GPa, and the mass density of 7860 kg/m3. The statistical properties of the first five measured natural frequencies of the plates are given in Table 1.Table 1
Statistical properties of the measured frequencies of the plates.
f 1 (Hz) f 2 (Hz) f 3 (Hz) f 4 (Hz) f 5 (Hz) mean 24.12 66.92 77.65 131.97 158.80 interval [23.94 24.4] [66.5 67.53] [76.24 78.91] [131.31 133.03] [156.31 160.94] Std. 0.11 0.25 0.57 0.42 0.97The FE model of the plate is created by the SHELL181 element in ANSYS, which has 300 shell elements. The boundary is free-free in the model. The material properties of the original FE model are the nominal values. The FE model and the first five mode shapes of simulation are shown in Figure6.Figure 6
The FE model and mode shapes of steel plates: (a) FE model; (b) mode 1; (c) mode 2; (d) mode 3; (e) mode 4; (f) mode 5.
(a) (b) (c) (d) (e) (f)After investigation, the uncertainty of rectangular plates can be characterized by Young’s modulusE, the shear modulus G, and the thick T in the FE model. The model updating is to identify the interval of these three uncertain parameters E, G, and T by the experiment results.Firstly, the intervals of the first five modal frequencies are estimated by KDE according to the experiment data in Ref. [30]. Secondly, it is assumed that the original intervals ofE,G, andT are [195220] GPa, [7887] GPa, and [1.301.60] mm, respectively, and intervals ofE,G, andT are identified according to the copula-based model updating flow with interval uncertainty in Figure 1. There are 20 experimental design data in each iteration step, i.e., 20 determined FE simulations in each iteration step. The model updating of the plate is convergent after four iteration steps with the convergent indices ε1=0.002 and ε2=0.002.After updating, the estimated intervals of the three parameters wereE=[203.21205.98] GPa, G=[82.7585.20] GPa, and T=[1.441.46] mm, respectively. The comparison between the original uncertain interval and the updated interval of the three parameters is shown in Figure 7, which shows that the updated intervals of these three parameters are much less than their original uncertain interval.Figure 7
Patch plots for the updating parameters of steel plates.
(a)
E vs G (b)
E vs TIn order to validate the updating results on the interval uncertain parameters, the copula models are reconstructed on the original and the updated intervals ofE, G, and T. The 5000 samples of E, G, T, and the first five natural frequencies of plates are obtained by MC sampling. The scatter plots for the simulation and the experimental results the first five natural frequencies of plates are shown in Figure 8. The comparisons on the natural frequencies between the simulation and the experimental results are listed in Table 2. It can be seen from Table 2 and Figure 8 that the frequencies of the updated model are in better agreement with the experimental results and the mean errors of frequencies decrease from the initial [12.5112.62]% to [0.720.24]%. Thus, this experimental validation has also confirmed the feasibility of the proposed IMU approach in estimating interval parameters of steel plates.Table 2
Comparison on the natural frequencies of plates between the simulation and the experimental results.
Mode Experimental interval/Hz Original interval/Hz Error /% Updated interval/Hz Error /% 1 [23.94 24.4] [21.06 27.55] [-12.05 12.92] [23.71 24.44] [-0.95 0.15] 2 [66.5 67.53] [58.61 76.53] [-11.86 13.32] [66.84 67.83] [-0.99 0.44] 3 [76.24 78.91] [65.81 87.80] [-13.68 11.26] [75.93 78.67] [-0.41 -0.30] 4 [131.31 133.03] [116.34 151.62] [-11.40 13.98] [130.37 134.26] [-0.71 0.92] 5 [156.31 160.94] [135.12 179.67] [-11.59 8.55] [155.51 160.95] [-0.51 0.01] mean   [12.51 12.62]  [0.72 0.24]Figure 8
Scatter plots for the first five natural frequencies of steel plates.
(a)
f1 vs f2 vs f3 (b)
f3 vs f4 vs f5


## 3.3. Example 3: Interval Model Updating in Transport Mirror System

The ShenGuangIII (SGIII) facility is designed for inertial confinement fusion (ICF) high energy experiments with 48 laser beams exactly transported and oriented to target. Figure9 is the view of the beam transport system in SGIII facility target area. There are 276 transport mirror systems in the facility, and the dynamic response under ambient vibration is a key factor to affect the stability of the SGIII facility [31]. A classic transport mirror system is shown in Figure 10, which consists of mirror component and support frame.Figure 9
A view of the beam transport system in SGIII facility target area, which is reproduced from Chen X J et al. (2014).Figure 10
Sketch map of the transport mirror system.The material of the support frame is steel with nominal Young’s modulus 200 Gpa, nominal Poisson’s ratio 0.3, and the nominal density 7850 kg/m3. The material of the mirror is K9 class with nominal Young’s modulus 80 Gpa, nominal Poisson’s ratio 0.21, and the nominal density 2510 kg/m3. The transport mirror system is about 1.32 m in height.The uncertainty exists in the transport mirror system because of the installation fluctuation and the welding technology fluctuation about the support frame, which induces the uncertainty of the vibration characteristics. In order to study how these uncertain factors affect the natural frequencies of the transport mirror system, impact hammer modal testing with fixed boundary conditions was conducted on ten nominally identical transport mirror systems. The experiment mount of the transport mirror system is shown in Figure11. Three natural frequencies in the interesting frequency range were found to significantly influence the response of the transport mirror system under the work condition. Then, the first three natural frequencies of the transport mirror system should be updated before calculating the response. The corresponding natural frequencies were obtained by the modal experiments, as listed in Table 3. It is observed that the frequency variations become more obvious with the increase of the mode order.Table 3
Experiment results of the first three natural frequencies of the transport mirror system.
Case 1 2 3 4 5 6 7 8 9 10 f 1/Hz 21.573 21.642 21.675 21.472 21.583 21.683 21.657 21.797 21.644 21.462 f 2/Hz 22.213 22.223 22.31 22.087 22.168 22.293 22.283 22.407 22.232 22.121 f 3/Hz 41.328 41.432 41.58 41.203 41.416 41.56 41.557 41.807 41.497 41.237Figure 11
Experiment setup of the transport mirror system.Generally, the bolts are ignored in the structural dynamic analysis model generally, and the FE model of the transport mirror system is established by the SOLID185 and SHELL181 element in ANSYS, which is shown in Figure12. A fully fixed bottom of the support frame is used as the boundary condition. The material properties of the original FE model are the nominal values. The first three simulation mode shapes are shown in Figure 13.Figure 12
FE model of the transport mirror system.Figure 13
Mode shapes of the transport mirror system.
(a)
Mode 1 (b)
Mode 2 (c)
Mode 3Theoretical analysis finds that the installation uncertainty of the transport mirror system can be characterized by Young’s modulusE1 of the bottom part of the support frame, and the weld uncertainty of the support frame can be characterized by Young’s modulus E2 of the upper part of the support frame. The goal of FE model updating is to identify the interval of E1 and E2.Firstly, the intervals of the first three modal frequencies are estimated by KDE according to the experiment data in Table3, which are [21.4121.83] Hz, [22.0322.45] Hz, and [41.0941.88] Hz, respectively.Secondly, assuming that the original intervals ofE1 andE2 are [100150] GPa and [160200] GPa, respectively, the intervals ofE1 andE2 are identified according to the copula-based model updating flow with interval uncertainty in Figure 1. There are 12 experimental design data in each iteration step, i.e., 12 determined FE simulations in each iteration step. The model updating of the transport mirror system is convergent after four iteration steps with the convergent indices ε1=0.005 and ε2=0.005.After updating, the estimated intervals of the three parameters wereE1 = [104.46114.65] GPa andE2 = [104.46114.65] GPa, respectively. The comparison between the original uncertain interval and the updated interval of the updating parameters is shown in Figure 14, indicating that the updated intervals of the updating parameters are much less than their original uncertain interval.Figure 14
Patch plots for the updating parameters of transport mirror system.In order to validate the updating results on the interval uncertain parameters, the copula models are reconstructed on the original and the updated intervals ofE1 andE2, Next, the 5000 samples ofE1,E2, and the first three natural frequencies of the transport mirror system are obtained by MC sampling. The scatter plots for the simulation and the experimental results of the first three natural frequencies of the transport mirror system are shown in Figure 15. The comparisons on the natural frequencies between the simulation and the experimental results are listed in Table 4. It can be seen from Table 4 and Figure 15 that the frequencies of the updated model agree with the experimental results better and the mean errors of frequencies decrease from the initial [2.4910.83]% to [0.470.40]%. As a result, this experimental validation has also proved the feasibility of the proposed IMU approach in estimating the interval parameters of the transport mirror system.Table 4
Comparison on the natural frequencies of the transport mirror system between the simulation and the experimental results.
Mode Experimental interval/Hz Original interval/Hz Error /% Updated interval/Hz Error /% 1 [21.46 21.80] [22.06 24.24] [2.77 11.20] [21.38 22.91] [-0.36 0.52] 2 [22.09 22.41] [22.68 24.94] [2.69 11.31] [22.00 22.54] [-0.38 0.59] 3 [41.20 41.81] [42.03 45.98] [2.01 9.97] [40.92 41.85] [-0.68 0.09] mean   [2.49 10.83]  [0.47 0.40]Figure 15
Scatter plots for the first three natural frequencies of transport mirror system.


## 4. Conclusions

In summary, a copula-based and MC sampling FE model updating approach in the structural dynamics with interval uncertainty is presented. Because the proposed FE model updating approach is based on the copula model of input-output parameters and MC sampling, the optimization difficulty caused by the interval algorithm is avoided in the proposed FE model updating approach with interval uncertainty, and the questions of ill-condition, nonuniqueness, and local optimal solution can be avoided also. Unlike other metal-based model updating approaches, the copula-based model updating approach can deal with the nonlinear relativity between the multi-input qualities and the multioutput qualities even if the multioutput qualities are correlated. In order to get a better copula model of input-output parameters for identifying the interval parameters, the LHS approach is suggested in the DOE and the number of samples is not less than ten in the proposed FE model updating approach. Besides, the proposed approach also considers the effect of small samples and estimates the interval of experiment samples by KDE in advance, which makes the updated model more reliable for the predictions on structural dynamic responses. Three examples were used to validate the FE model updating approach proposed in the study. The updating results have proved the feasibility and reliability of this approach.The total time consumption of the FE model updating is mainly the time of the FE model simulation. For example, the FE model updating of the transport mirror system was completed after four iteration steps, and there were 12 FE model simulations in each iteration step, so only 48 deterministic FE simulations were carried out in the whole model updating with interval uncertainty. But, more than hundreds of FE simulations are needed in a determined model updating without the surrogate model, and the number of FE simulations in the model updating with uncertainty is far more than that of a deterministic model updating. Then, the model updating method proposed with interval uncertainty is efficient in the paper. Besides, after updating the FE model, the difference of the first few natural frequencies between the calculated interval and the experimental interval is less than 1%. In conclusion, the copula-based and MC sampling concept expands the application range of IMU and is useful for fast and accurate estimation of interval parameters in the FE model.



---

*Source: 3958016-2018-07-09.xml*