Model Updating Using Frequency Response Functions Based on Sherman–Morrison Formula

Model updating plays an important role in dynamics modeling with high accuracy, which is widely used in mechanical engineering. In this paper, a model updating method using frequency response function (FRF) is proposed based on Sherman–Morrison formula, in which only the initial FRFs and parameter perturbations are employed to calculate the sensitivity avoiding repeated finite element (FE) analyses and improving the computational efficiency. Firstly, the sensitivity of FRFs to the design parameters is calculated by Sherman–Morrison formula based on the QR decomposition of the system dynamic stiffness matrix variation after parameter perturbations, then the influence of damping on the amplitude of FRFs is considered to select an appropriate frequency range, and finally conduct the model updating according to the sensitivity of the FRFs. By employing simulation examples of a truss and a solar wing and the experiment of an aluminum frame, the updating error is still within ±1.00% in the condition of 5% random white noise, which shows the proposed method has high accuracy and a certain anti-noise capability. When only a few numbers of frequency points are selected near the resonance peak of the FRFs, the result shows that selecting the appropriate frequency range and points can reduce the computational cost. The results of the experiment study show that the proposed method can effectively identify the structural parameters. The above results verify the feasibility and effectiveness of proposed model updating method using FRFs.


Introduction
Finite element model with high accuracy plays an important role in structural dynamics design, which is the premise of dynamic response prediction, vibration control, and structural optimization [1][2][3]. Due to the idealization hypothesis of structural simplification, the errors of structural parameters, and the errors of finite element (FE) discretization, the FE models may not accurately reflect the dynamics performance of the structures [4][5][6]. With the rapid development of numerical simulation in engineering applications, model updating methods gradually revealed its importance, which makes the dynamic characteristics of the theoretical model as consistent as possible with the actual structure.
As a bond between FE models and actual structures, model updating has become a powerful tool and which has been paid extensive attention [7,8]. Many types of approaches have been developed for a variety of problems. Shabbir et al. [9] proposed a method based on particle swarm optimization to update the FE model by using sequence niche technology, which greatly improved the reliability of finding the global minimum; Behmanesh et al. [10] proposed a new probabilistic FE model updating method based on hierarchical Bayesian modeling for identification of structural parameters 2 of 23 under varying environmental conditions. Taking the strain mode as the target response, Zhan and Guo et al. [11] constructed the target function to be updated by using the strain mode frequency and modal confidence errors, the parameters were updated by using the Genetic Algorithm. Aiming at the model updating of the nonlinear structure, Wang and Ren et al. [12] took the instantaneous characteristics, including instantaneous frequency and instantaneous amplitude, as the objective function of the nonlinear model updating, and used the annealing global optimization method to obtain the optimal value of the nonlinear parameters. Zang et al. [13] proposed a model updating method that separated the non-linear characteristics of subsystem from the parameter uncertainties to establishing a representative linear mathematical model and applied it to a structural dynamics verification problem which developed by Sandia National Laboratory. The above methods focus on the recent hot spots to carry out research, frequency response function (FRF)'s based model updating is an alternative method with some advantages, FRFs can adequately reflect the dynamic characteristics of a structure, using which as the feature will get rid of modal identification.
Model updating based on FRFs avoids the error caused by modal fitting, and match of the predicted modal shapes with the measured data are not necessary, so it is suitable for structural updating with relatively dense modal distribution. Shadan et al. [14,15] proposed a new structural damage detection method based on FRFs data and natural frequencies, adopted pseudo-linear sensitivity equation to accurately identify the location and severity of damage; Esfandiari et al. [16,17] proposed an estimation method for structure mass and stiffness considering damping effect and using the FRF data for FE model updating; Pradhan et al. [18] proposed a FE model updating method that can identify and update the damping matrix of the system, and studied the updating capability of different damping levels through three different experimental beam structures; Li et al. [19] proposed an iterative method combining model updating method with model reduction technology, effectively improved the efficiency of the updating; Hong et al. [20] successfully updated the stiffness, mass and damping parameters between the analytical and the experimental FRFs through the laboratory test of the four-story shear frame structure; Wang et al. [21,22] conducted the model updating of local nonlinear structure based on the initial FE model of the structure and the primary harmonic response data obtained from low amplitude and high amplitude excitation. The above scholars have done a lot of research on the model updating based on FRFs. However, the selection of parameters plays a significant role in the progress of model updating. The sensitivity of the FRFs has great potential for a vibration analysis, which is not involved in the above methods. Reasonable selection of parameters can effectively improve the quality and efficiency of model updating, which is inseparable from the sensitivity analysis of the parameters, so model updating method based on the sensitivity of FRFs has more important significance.
The model updating based on the sensitivity of FRFs attracted wide attention in the early 1990 s [23]. Sipple et al. [24] proposed a new FE model updating method that uses numerical sensitivity instead of analytical sensitivity to solve inverse problems, which does not need model reduction and data extension, and has good robustness, stability and fault tolerance; Weng et al. [25] proposed a new iterative substructuring method to accurately obtain the characteristic solution and characteristic sensitivity of the structure, which is suitable for model updating of large and complex structures; Jiang et al. [26] expressed the thermal dependent elastic constants and thermal expansion coefficients as intermediate function by temperature independent variables, and used the perturbation method to calculate the sensitivity for parameter identification; On the basis of the basic concept of unconstrained optimization problem and three regularization solutions, Mohammad et al. [27] proposed a novel sensitivity based FE model updating method, and simultaneously updated the element mass matrix and stiffness matrix of the FE model; Esfandiari et al. [28] proposed a FE model updating method using incomplete strain data in frequency domain to detect variations in stiffness and mass parameters. The common sensitivity analysis methods using FRFs are the finite difference method [29] or the combined approximations (CAs) method [30]. These methods are based on obtaining the stiffness matrix, damping matrix, and mass matrix of the structures. The disadvantage of these methods is that repeated analysis is required in the process of structural parameter identification, and they are not applicable under the premise of unknown dynamic stiffness matrix. For example, the modal parameters of the structure can be obtained through experiments, but the complete information of the dynamic stiffness matrix is difficult to obtain. The Sherman-Morrison formula was first proposed by Sherman and Morrison and it can be used to represent the varying of matrix; Turan et al. [31] used the Sherman-Morrison formula of singular value decomposition to calculate the optimal perturbation of the sensitivity value of the extreme value and the design variable, thus achieving the response change that is difficult to obtain with the response derivatives; Zhu et al. [32][33][34][35] used the Sherman-Morrison formula for mass elimination of multi-transducers, and the revised FRFs results were in good agreement with the target values.
In this study, a model updating method using FRFs based on Sherman-Morrison formula is proposed, only the initial FRFs and parameter perturbations are employed to calculate the sensitivity avoiding repeated FE analyses and calculations, which improves the computational efficiency of model updating. The detail formulation of sensitivity analysis and model updating based on Sherman-Morrison formula is introduced in Section 2, the implementation procedure of the model updating is described in Section 3, and the numerical examples of a planar truss and a solar wing and the experiment of aluminum solar panel frame are used to verify the proposed method in Section 4. The conclusions are summarized in Section 5.

Theoretical Analysis
This section verifies that the model updating method using FRFs based on Sherman-Morrison formula can obtain the sensitivity of the FRFs to the parameters only by initial FRFs and parameter perturbations, which can avoid repeated finite element analyses and calculations.

Sensitivity Analysis Based on SMW
The dynamic equation of the multi-degrees-of-freedom system in the frequency domain can be expressed as where M is the mass matrix, K is the stiffness matrix, C is the damping matrix, x(ω) and F(ω) are the Fourier transformation of displacement response and excitation force, respectively. The dynamic stiffness matrix Z(ω) of the system is The relation between FRFs H(ω) and dynamic stiffness matrix Z(ω) is If the FRFs data H(ω) of the system is a function of a design parameter P m , then the sensitivity S of the FRFs data H(ω) of the system is the first partial derivative of the FRFs with respect to the design parameter is ∂H ∂P m , which be expressed as According to the formula, the variation of FRFs is the variation of dynamic stiffness. For a specific frequency ω r , it can be obtained from Equation (4) that the first derivative of dynamic stiffness with respect to a variable P m is: Appl. Sci. 2020, 10, 4985 4 of 23 Equation (5) can be expressed in differential form as where δZ, δM, δC, and δK represent the variations of dynamic stiffness matrix, mass matrix, damping matrix, and stiffness matrix, respectively. However, the damping of the structure is related to the material performance of individual members and the difference of FE modeling. The influence of damping on the FRF amplitude is significant in resonance and decreases rapidly as the damping moves away from resonance. Therefore, the selection of appropriate frequency range can mitigate the impact of damping on the model updating results.
We assume that a system of n-degrees of freedom is subject to mass perturbation δm i at the ith-degree of freedom, and stiffness perturbation δk i at the jth-degree of freedom [32], the perturbations of the mass matrix δM and the stiffness matrix δK can be expressed as We can rewrite the perturbations of the mass matrix δM and the stiffness matrix δK as Where u and v represent the nodes on each side of the unit. G uv and X uv are expressed as Appl. Sci. 2020, 10, 4985 5 of 23 Therefore, the variation of the system dynamic stiffness matrix δZ can be expressed as The inverse for the variation of dynamic stiffness matrix δZ which is the variation of FRFs δH is obtained by applying the Sherman-Morrison formula where L is a 2N × 2N matrix given by By substituting Equation (13) into Equation (4), the sensitivity of the FRFs with respect to perturbation parameters can be obtained. Therefore, the sensitivity of FRFs to the perturbation parameters can be directly solved by only the initial FRFs and parameter perturbations to avoid repeated FE analyses and calculations.
The FRFs H*(ω) of the system changed by parameter perturbation can be obtained according to Sherman-Morrison formula: Therefore, the sensitivity of FRFs based on Sherman-Morrison formula can be expressed as where Z*(ω) is the dynamic stiffness matrix changed by parameter perturbation, ∆Z(ω) is the variation of the dynamic stiffness matrix, which is caused by the change quantity ∆P, and ∆Z(ω) can be decomposed into two matrices U(ω) and V T (ω), I is the n-order identity matrix.

Frequency Response Function (FRF)'s Model Updating Based on Sensitivity Analysis
The purpose of the model updating is to minimize the difference of FRFs between the experimental and the FE model. The difference between the FRFs of experimental model H*(ω) and the FRFs of FE model H(ω) are used to construct the objective function as follows: where δ is the difference between the experimental FRFs and the corresponding analysis values, which is called the residual vector. Equation (4) is rewritten as: Appl. Sci. 2020, 10, 4985 6 of 23 and δ = ∆H(ω) = −H * (ω) ∆K + jω r ∆C − ω r 2 ∆M H(ω) (20) where the changes of the total mass matrix ∆M and the total stiffness matrix ∆K of the structure can be expressed by the sum of partial derivatives of each updating element: where nm and ns are the total number of mass and stiffness updating units, and ∆P i M represents the change of the ith mass updating parameter, ∆P j S represents the change of the jth stiffness updating parameter. By substituting Equations (21) and (22) into Equation (20), and setting ∆C = 0, which means ignoring the updating of damping parameter Equation (23) is rewritten as follows where S M and S S are the sensitivity matrix of mass and stiffness, respectively, S is the sensitivity matrix. ∆P M and ∆P S are the variation of mass and stiffness parameters, and ∆P is the vector of all the variation of stiffness and mass parameters. The model updating of parameters is an iterative process of convergence, and the regularization process is required to solve Equation (24). Generally, the following estimation criteria are obtained by adding constraint information for the variation of updating parameters: where λ 2 is the regularization parameter whose value can be determined by the L-curve method. It can be obtained that the estimated variation of the updating parameters according to Equation (25): Therefore, the iterative formula of the updating parameters is: where n and n + 1, respectively, represents the values corresponding to the analysis model in the iteration of step n and n + 1. Therefore, the core of the model updating is how to establish the sensitivity matrix S. It is not difficult to find from Equations (4) and (5) that for multi-degree-of-freedom systems, the solution of the sensitivity matrix problem finally boils down to the sensitivity of the system mass and stiffness to the design parameters.
Conducting the model updating progress by the obtained sensitivity S, and the FRFs was iteratively updated by using Equation (27) to obtain the finite element model consistent with the actual structure. A simulation example of a truss and a solar wing is used to verify the effectiveness of this method in this paper.

Frequency Range Selection
The FRFs of the real experimental model and the FE model around the nth natural frequency is shown in Figure 1. The modeling errors will cause the FRFs of the same mode order to move along the coordinate axis to a certain extent, but the basic shape does not change significantly. The approximate estimation of the FRFs Hd(ω) near the resonance region can be expressed as where , , , where, ωn and ωdn are the nth natural frequency of the real experimental model and the FE model, respectively, δωn is the frequency variation of the real experimental model and the FE model, and m is the order of the natural frequency taken by simulation.
The frequency range of model updating is expressed as ' ω '. According to Equations (28) and (29), Hd(ω) on ' ω ' can be approximately expressed as After repeated iterative calculations, the structural parameters of the FE model will be constantly updated and gradually approach to the parameters of the real experimental model. Therefore, in Equations (28) and (29), the initial frequency ωn and FRFs H(ω) will be constantly updated iteratively and converge to the natural frequency and FRFs of the real experimental model.
At present, most of the damping models cannot accurately simulate the actual damping effect but will increase the model errors and affect the result of model updating. Therefore, the damping parameter is not updated in this paper, so the influence of structural damping on the result of model updating should be reduced. From the analysis in Section 2.2, the damping has little effect on the frequency response amplitude when far away from the resonance peak. Therefore, the selection of appropriate frequency range can alleviate the influence of damping on the results of model updating, and the influence of damping can be ignored when the model is updated based on FRFs. The frequency close to resonance would be affected by damping, so the small region of about 4Hz around resonance was not considered in the model updating process [14].
If the two adjacent resonant peaks are too close to each other, the variation of FRFs in the frequency band between them cannot be guaranteed to decrease monotonically in the updating process, so this region is excluded to avoid the amplitude nonlinear problem. In each iteration, the frequency range also needs to be adjusted according to the change of the resonance point.  The approximate estimation of the FRFs H d (ω) near the resonance region can be expressed as

Implementation Procedure
where δω n = ω n − ω dn , n = 1, · · · , m where, ω n and ω dn are the nth natural frequency of the real experimental model and the FE model, respectively, δω n is the frequency variation of the real experimental model and the FE model, and m is the order of the natural frequency taken by simulation.
The frequency range of model updating is expressed as 'ω'. According to Equations (28) and (29), H d (ω) on 'ω' can be approximately expressed as After repeated iterative calculations, the structural parameters of the FE model will be constantly updated and gradually approach to the parameters of the real experimental model. Therefore, in Equations (28) and (29), the initial frequency ω n and FRFs H(ω) will be constantly updated iteratively and converge to the natural frequency and FRFs of the real experimental model.
At present, most of the damping models cannot accurately simulate the actual damping effect but will increase the model errors and affect the result of model updating. Therefore, the damping parameter is not updated in this paper, so the influence of structural damping on the result of model updating should be reduced. From the analysis in Section 2.2, the damping has little effect on the frequency response amplitude when far away from the resonance peak. Therefore, the selection of appropriate frequency range can alleviate the influence of damping on the results of model updating, and the influence of damping can be ignored when the model is updated based on FRFs. The frequency close to resonance would be affected by damping, so the small region of about 4Hz around resonance was not considered in the model updating process [14].
If the two adjacent resonant peaks are too close to each other, the variation of FRFs in the frequency band between them cannot be guaranteed to decrease monotonically in the updating process, so this region is excluded to avoid the amplitude nonlinear problem. In each iteration, the frequency range also needs to be adjusted according to the change of the resonance point. Figure 2 is a flowchart of the model updating method using FRFs based on Sherman-Morrison formula, with the steps described below:
Establish the FE analysis model, and extract the mass matrix M 1 and the stiffness matrix K 1 of the original model to obtain the dynamic stiffness matrix Z(ω) of the system, and calculate the FRFs, that is H(ω); 2.
Select several parameters that have an impact on the dynamic characteristics of the structure as the initial variables P 1 , and P 1 become P 2 by perturbations, and the variation is ∆P; extract the mass matrix M 2 and the stiffness matrix K 2 of the changed model to obtain the variation of dynamic stiffness matrix of the system, that is ∆Z(ω), and recalculate the FRFs H*(ω) after parameter perturbations by using Sherman-Morrison formula; 3.
Select the appropriate frequency range which is away from the resonance for model updating and obtain the sensitivity S of the FRFs to the designed parameters. 4.
Conduct the model updating by the obtained sensitivity S, and the FRFs was iteratively updated by using Equation (27), obtain the variation P n+1 − P n , if the variables P n are converged, stop; otherwise set n = n + 1, go to step (2); the stop criterion is P n+1 − P nn ≤ 10 −6 ; parameter perturbations by using Sherman-Morrison formula; 3. Select the appropriate frequency range which is away from the resonance for model updating and obtain the sensitivity S of the FRFs to the designed parameters. 4. Conduct the model updating by the obtained sensitivity S, and the FRFs was iteratively updated by using Equation (27), obtain the variation Pn+1 − Pn, if the variables Pn are converged, stop; otherwise set n = n + 1, go to step (2); the stop criterion is

Case Study
The simulation examples of the truss and the solar wing and the experiment of aluminum solar panel frame are used to verify the feasibility of the model updating method using FRFs based on Sherman-Morrison formula in this section.

Truss Structure
In this section, the FE model of truss structure is used for model updating, and considering the results under the interference of noise at the same time, so as to verify the feasibility of the model updating method using FRFs based on Sherman-Morrison formula.

Finite Element (FE) Modeling
The truss structure model is shown in Figure 3, numerically modeled by the 'Matlab' programming software. To simply to verify the feasibility of the method, each beam is discretized into a beam element, there are 21 discrete beam elements and 12 nodes, each node contains translational degrees of freedom (DOFs) in X direction and Y direction. Due to the Node 1 is fixed and Node 12 is fixed on the hinge support, so the translational DOFs x 1 and y 1 in X direction and Y direction of Node 1 and translational DOFs y 12 in Y direction of Node 12 are locked, so the total DOFs of the truss structure is 21. It is assumed that the DOFs of the target model that need to be updated and the FE model are matched, so there is no problem of model reduction and modal expansion. Therefore, the degrees of freedom of both the FE model and the target model is 21. The elastic modulus of all beam elements in the initial FE analysis is E = 70 Gpa, the density of the truss-element material is ρ = 2700 kg/m 3 , the length of each beam is L = 1 m, and the cross-sectional area of each truss element is A = 1 × 10 −4 m 2 .
into a beam element, there are 21 discrete beam elements and 12 nodes, each node contains translational degrees of freedom (DOFs) in X direction and Y direction. Due to the Node 1 is fixed and Node 12 is fixed on the hinge support, so the translational DOFs x1 and y1 in X direction and Y direction of Node 1 and translational DOFs y12 in Y direction of Node 12 are locked, so the total DOFs of the truss structure is 21. It is assumed that the DOFs of the target model that need to be updated and the FE model are matched, so there is no problem of model reduction and modal expansion. Therefore, the degrees of freedom of both the FE model and the target model is 21. The elastic modulus of all beam elements in the initial FE analysis is E = 70 Gpa, the density of the truss-element material is ρ = 2700 kg/m 3 , the length of each beam is L = 1 m, and the cross-sectional area of each truss element is A = 1 × 10 −4 m 2 . It is assumed that the FE model of the truss structure has modeling errors, and the elastic modulus of beam element 1, 6, 11, 14, and 21 of the target model is 1.1, 1.2, 1.2, 1.2, and 1.1 times of the FE model, which is bold in red in Figure 3. It is assumed that the simulation experiment is a single point excitation, and the excitation point and response point in the target model is the same one. Using the above data for FE analysis, the frequencies of the first 7 order and FRFs of Node 4 and Node 10 are selected as the target value for model updating. Because the planar truss system is modeled in undamped mode, the influence of damping on FRFs at resonance peak need not be considered in the updating progress.

Comparison of Sensitivity Analysis
To further more intuitively illustrate the accuracy and efficiency of the proposed method in sensitivity calculating, we make 10 −1 , 10 −2 , and 10 −3 of the perturbations in the elastic modulus respectively and compare the proposed method with the analytic derivatives method of Equation (4) to obtain the results shown in Figure 4. It is assumed that the FE model of the truss structure has modeling errors, and the elastic modulus of beam element 1, 6, 11, 14, and 21 of the target model is 1.1, 1.2, 1.2, 1.2, and 1.1 times of the FE model, which is bold in red in Figure 3. It is assumed that the simulation experiment is a single point excitation, and the excitation point and response point in the target model is the same one. Using the above data for FE analysis, the frequencies of the first 7 order and FRFs of Node 4 and Node 10 are selected as the target value for model updating. Because the planar truss system is modeled in undamped mode, the influence of damping on FRFs at resonance peak need not be considered in the updating progress.

Comparison of Sensitivity Analysis
To further more intuitively illustrate the accuracy and efficiency of the proposed method in sensitivity calculating, we make 10 −1 , 10 −2 , and 10 −3 of the perturbations in the elastic modulus respectively and compare the proposed method with the analytic derivatives method of Equation (4) to obtain the results shown in Figure 4.  The following error index O is used to quantify the accuracy of the sensitivity of the FRFs: where N is the number of points on the curves. Si Exact and Si Proposed are the sensitivities of the FRFs obtained by analytical derivatives method of Equation (4)    The following error index O is used to quantify the accuracy of the sensitivity of the FRFs: where N is the number of points on the curves. S i Exact and S i Proposed are the sensitivities of the FRFs obtained by analytical derivatives method of Equation (4) and the proposed method, respectively. A smaller O indicates a more accurate result.
The results in Figure 4 show that under different perturbations, the sensitivity calculated by the proposed method compares reasonably well with that calculated by the analytical derivatives. The small figure in the blue dashed box in Figure 4 is a local magnification comparison of the sensitivity curves calculated by the proposed method and the analytical derivatives method under different perturbations. With the decrease of the perturbation value, the sensitivity calculated by the proposed method is closer to that of the analytic derivatives method.
As shown in Table 1, the error is only 0.01% for the perturbation of 10 −3 , which shows that the proposed method can carry out sensitivity analysis sufficient accurately for small variations in the design parameters. It can also be found that although the sensitivity calculated by the two different methods is basically the same, the calculation efficiency of the proposed method in this paper is higher than that of the analytic derivatives method as is shown in Table 2. This advantage will be of great use for model updating of large-scale engineering structures, which can save much computation time. Due to the inevitable interference of various noises in the real experiment, in order to simulate the real experiment state and test the anti-noise interference ability of the method proposed in this paper, 5% random white noise was added to the target value which was obtained in Section 4.1.1 in the form of H**(ω) = H*(ω) (1 + 0.05 µ), where the µ is a random number of −1~1, H*(ω) is the target value obtained in Section 4.1.1, and H**(ω) is the FRFs after adding random white noise. Then, the parameter updating formula shown in Equation (27) was used for 15 iterative computations. The frequency comparison between the FE model and the target model before and after modal updating is shown in Table 3, and the comparison results of the updated parameters are shown in Table 4. The convergence process of updated parameters is shown in Figure 5. The comparison of FRFs before and after updating is shown in Figures 6 and 7.

Parameter Updating
Due to the inevitable interference of various noises in the real experiment, in order to simulate the real experiment state and test the anti-noise interference ability of the method proposed in this paper, 5% random white noise was added to the target value which was obtained in Section 4.1.1 in the form of H**(ω) = H*(ω) (1 + 0.05 μ), where the μ is a random number of −1~1, H*(ω) is the target value obtained in Section 4.1.1, and H**(ω) is the FRFs after adding random white noise. Then, the parameter updating formula shown in Equation (27) was used for 15 iterative computations. The frequency comparison between the FE model and the target model before and after modal updating is shown in Table 3, and the comparison results of the updated parameters are shown in Table 4. The convergence process of updated parameters is shown in Figure 5. The comparison of FRFs before and after updating is shown in Figures 6 and 7.

Parameter Updating
Due to the inevitable interference of various noises in the real experiment, in order to simulate the real experiment state and test the anti-noise interference ability of the method proposed in this paper, 5% random white noise was added to the target value which was obtained in Section 4.  1, and H**(ω) is the FRFs after adding random white noise. Then, the parameter updating formula shown in Equation (27) was used for 15 iterative computations. The frequency comparison between the FE model and the target model before and after modal updating is shown in Table 3, and the comparison results of the updated parameters are shown in Table 4. The convergence process of updated parameters is shown in Figure 5. The comparison of FRFs before and after updating is shown in Figures 6 and 7.   The model updating process uses the frequencies of the first 5 order of the target value and their corresponding FRFs, while the fifth and sixth frequencies and their corresponding FRFs do not directly participate in the updating, but are used to evaluate the accuracy of the model updating method which is proposed in this paper.
Where Errorb and Errora is the relative error between the frequency values of the FE model before and after the updating and the target value, respectively.
As can be seen from Table 3, under the excitation of 5% random white noise, there is a significant difference between the frequency of the FE model and the target model before the modal updating, and the maximum frequency error has reached −3.26%. This indicates that the reliability of the FE model is questionable, so it is necessary to update the initial FE model to improve the accuracy and reliability. After using the Sherman-Morrison formula proposed in this paper, the maximum error of the frequency of the first 7 order after modal updating is only −0.50%, which verifies that the model updating method using FRFs based on Sherman-Morrison formula has a good anti-noise performance. The model updating process uses the frequencies of the first 5 order of the target value and their corresponding FRFs, while the fifth and sixth frequencies and their corresponding FRFs do not directly participate in the updating, but are used to evaluate the accuracy of the model updating method which is proposed in this paper.
Where Error b and Error a is the relative error between the frequency values of the FE model before and after the updating and the target value, respectively.
As can be seen from Table 3, under the excitation of 5% random white noise, there is a significant difference between the frequency of the FE model and the target model before the modal updating, and the maximum frequency error has reached −3.26%. This indicates that the reliability of the FE model is questionable, so it is necessary to update the initial FE model to improve the accuracy and reliability. After using the Sherman-Morrison formula proposed in this paper, the maximum error of the frequency of the first 7 order after modal updating is only −0.50%, which verifies that the model updating method using FRFs based on Sherman-Morrison formula has a good anti-noise performance.
Where Error b and Error a is the relative error between the elastic modulus of the FE model before and after the updating and the target value, respectively.
The model updating results of elastic modulus are shown in Table 4, under the excitation of 5% random white noise, the maximum relative error of the corrected elastic modulus is changed from 20.00% to −0.85%. Compared with the result in noiseless environment, the maximum error of the elastic modulus is within 1%, which greatly improves the reliability of the FE model and verifies the feasibility of the model updating method in this paper.
It can also be seen from the iterative convergence diagram of elastic modulus E of the truss in Figure 5, under the excitation of 5% random white noise, after 15 iterative computations, each parameter tends to converge, and the maximum error of each parameter is within ±1.00% after convergence. Therefore, the results verify that the model updating method using FRFs based on Sherman-Morrison formula is still applicable in noisy environment.
As shown in Figures 6 and 7, the FRFs before and after updating are compared with the target value under 5% random white noise. There is obviously a big gap between the FRFs of the FE model before updating and the target model under 5% random white noise, and the FRFs data of the FE model almost completely coincides with the target value after updating, which shows that the model updating method using FRFs based on Sherman-Morrison formula can effectively reduce noise interference.
Each iteration is based on the sensitivity of the response at all the frequency points involved in the model updating. However, too many frequency points involved in the updating will greatly increase the computation and affect the efficiency of model updating. The change of FRFs is often discontinuous at the resonance peak, which especially happens in the low-damping structure. Therefore, the iterative process will not converge when the response of the frequency points close to the resonance peak are used for updating. So, it is generally necessary to select the frequency point carefully before model updating.
Taking the FRFs of Node 10 under 5% random white noise as an example, three frequency points selection cases are presented are given as follows, and the influence of them on the results of model updating is compared by observing the convergence of the elastic modulus of each element.

•
Case 1: Select one frequency point every 10 Hz in the frequency band from 0 to 400 Hz and select 40 points in total. • Case 2: Select 40 points randomly in the frequency band from 0 to 400 Hz. • Case 3: Select 8 frequency points at the interval of 0.5 Hz from 10 Hz to the left of each resonance peak in the frequency band from 0 to 400 Hz and select 40 points in total. Figure 8 shows the comparison of elastic modulus updating results of three cases. Since the frequency points selected by Case 1 and Case 2 do not avoid the resonance peak, so the FRFs are discontinuous which causes the elastic modulus of the two cases is not convergent. While only after 7 iterations, the elastic modulus of Case 3 is convergent. However, due to the noise pollution of 5% in the FRFs data, the result does not converge to 0%, but the maximum error is within ±1.5%.
updating method using FRFs based on Sherman-Morrison formula can effectively reduce noise interference.
Each iteration is based on the sensitivity of the response at all the frequency points involved in the model updating. However, too many frequency points involved in the updating will greatly increase the computation and affect the efficiency of model updating. The change of FRFs is often discontinuous at the resonance peak, which especially happens in the low-damping structure. Therefore, the iterative process will not converge when the response of the frequency points close to the resonance peak are used for updating. So, it is generally necessary to select the frequency point carefully before model updating.
Taking the FRFs of Node 10 under 5% random white noise as an example, three frequency points selection cases are presented are given as follows, and the influence of them on the results of model updating is compared by observing the convergence of the elastic modulus of each element.

•
Case 1: Select one frequency point every 10 Hz in the frequency band from 0 to 400 Hz and select 40 points in total. • Case 2: Select 40 points randomly in the frequency band from 0 to 400 Hz. • Case 3: Select 8 frequency points at the interval of 0.5 Hz from 10 Hz to the left of each resonance peak in the frequency band from 0 to 400 Hz and select 40 points in total. Figure 8 shows the comparison of elastic modulus updating results of three cases. Since the frequency points selected by Case 1 and Case 2 do not avoid the resonance peak, so the FRFs are discontinuous which causes the elastic modulus of the two cases is not convergent. While only after 7 iterations, the elastic modulus of Case 3 is convergent. However, due to the noise pollution of 5% in the FRFs data, the result does not converge to 0%, but the maximum error is within ±1.5%.  By comparing the results of model updating using all the FRFs data, it can be seen that only a small amount of FRFs data is needed to complete the model updating, which indicates that the correct selection of frequency points can effectively improve the computational efficiency of model updating, and further proves that the proposed method has high computational efficiency and accuracy. Furthermore, only a small amount of FRFs data can be selected to complete the model updating quickly, the result is not very perfect compared with the result of using all frequency points, which indicates that increasing the frequency point can correspondingly increase the information contained in the FRFs data, which is an effective means to improve the quality of model updating.

Solar Wing
In this section, the FE model of solar wing is used for model updating to verify the feasibility of the model updating method using FRFs based on Sherman-Morrison formula.

FE Modeling
A complete FE model of the solar wing was modeled by the FE analysis software 'Nastran' based on the size diagram of the solar wing model as shown in Figure 9. The solar wing is an up-down axisymmetric structure, consisting of two plates with identical materials and geometric shapes. Each plate is divided into 12 shell elements numbered 1~24, and the two plates are connected by beam elements numbered 25~34. The solar wing model consists of 34 elements and 47 nodes, each node contains translational and rotational DOFs in the direction of X, Y, and Z. Since Node 44, 45, and 46 are fixed, and the rest are hinged, the translational and rotational DOFs in the X, Y, and Z directions of these three nodes are all locked. Therefore, the total DOFs of the solar wing structure is 264. The material parameters of the finite element model are as follows: the elastic modulus of the beam is E 1 = 200 Gpa, the density is ρ 1 = 7850 kg/m 3 , the cross-sectional area of the beam is A 1 = 4 × 10 −4 m 2 , and the poisson ratio is µ 1 = 0.32. The elastic modulus of the solar wing is E 2 = 70 Gpa, the density is ρ 2 = 2700 kg/m 3 , the thickness is t = 4 × 10 −3 m, and the poisson ratio is µ 2 = 0.34. To simulate the actual damping of the model, the Rayleigh proportional damping model is adopted, as is shown below where C represents the damping matrix of the system, M represents the mass matrix of the system, K represents the stiffness matrix of the system, a and b represent the damping coefficients. In this case, the actual damping coefficients are assumed to be a = 0.2 and b = 0.0005.

FE Modeling
A complete FE model of the solar wing was modeled by the FE analysis software 'Nastran' based on the size diagram of the solar wing model as shown in Figure 9. The solar wing is an up-down axisymmetric structure, consisting of two plates with identical materials and geometric shapes. Each plate is divided into 12 shell elements numbered 1~24, and the two plates are connected by beam elements numbered 25~34. The solar wing model consists of 34 elements and 47 nodes, each node contains translational and rotational DOFs in the direction of X, Y, and Z. Since Node 44, 45, and 46 are fixed, and the rest are hinged, the translational and rotational DOFs in the X, Y, and Z directions of these three nodes are all locked. Therefore, the total DOFs of the solar wing structure is 264. The material parameters of the finite element model are as follows: the elastic modulus of the beam is E1 = 200 Gpa, the density is ρ1 = 7850 kg/m 3 , the cross-sectional area of the beam is A1 = 4 × 10 −4 m 2 , and the poisson ratio is μ1 = 0.32. The elastic modulus of the solar wing is E2 = 70 Gpa, the density is ρ2 = 2700 kg/m 3 , the thickness is t = 4 × 10 −3 m, and the poisson ratio is μ2 = 0.34. To simulate the actual damping of the model, the Rayleigh proportional damping model is adopted, as is shown below (32) where C represents the damping matrix of the system, M represents the mass matrix of the system, K represents the stiffness matrix of the system, a and b represent the damping coefficients. In this case, the actual damping coefficients are assumed to be a = 0.2 and b = 0.0005.  Table 5. The modeling errors of the target model that needs to be updated and the FE model are shown in Table 6. There are initial errors of 30.00%, 20.00%, 15.00%, and 20.00% in E1, E2, ρ1, and A1, respectively. It is assumed that the simulation experiment is a single point excitation, and the excitation point and response point in the target model is the same one. The frequency of the  Table 5. The modeling errors of the target model that needs to be updated and the FE model are shown in Table 6. There are initial errors of 30.00%, 20.00%, 15.00%, and 20.00% in E 1 , E 2 , ρ 1 , and A 1 , respectively. It is assumed that the simulation experiment is a single point excitation, and the excitation point and response point in the target model is the same one. The frequency of the first 7 order and FRFs of Node 13 and Node 26 are obtained as the target value to conduct the model updating progress.

Parameter Updating
The elastic modulus of the connecting beams and the solar wing, and the density and the cross-sectional area of the connecting beams are updated in turn with the strategy of step by step. Meanwhile, 5% random white noise was added to the target value of the FRF data obtained in The model updating process uses the frequencies of the first 5 order of the target model and their corresponding FRFs, while the fifth and sixth frequencies and their corresponding FRFs do not directly participate in the updating, but are used to evaluate the accuracy of the model updating method which is proposed in this paper. Combined with the rules in Section 2.3, the frequency of the model updating is selected as the frequency point near the resonance region that is shown in Figure 10. Remove the low frequency band containing the first two resonance regions, and to avoid the problem of amplitude nonlinearity, the frequency points between the fourth and fifth resonances are ignored.  Where Errorb and Errora is the relative error between the frequency values of the FE model before and after the updating and the target value, respectively.
As can be seen from Table 5, under the excitation of 10% random white noise, the maximum relative error of the frequency is changed from −5.49% to −1.10%, which verifies the effectiveness and anti-noise performance of the proposed method in this paper.
Where Errorb and Errora is the relative error between the parameters of the FE model before and after the updating and the target value, respectively.
The model updating results of parameters are shown in Table 6, under the excitation of 10% random white noise, the maximum relative error of the updated parameters is changed from 30.00% to −0.24%, which greatly improves the reliability of the FE model and verifies the feasibility of the model updating method in this paper.
As shown in Figures 11 and 12, due to modeling errors between the target model and the FE model, there is obviously a big gap between the FRFs of the FE model before updating and the target model under 5% random white noise. While after 20 iterative computations, the error between the obtained FRFs and the target value is quite close, which shows that the model updating method using FRFs based on Sherman-Morrison formula can effectively reduce noise interference. At the same time, the influence of damping on the amplitude of FRFs is effectively reduced by choosing the appropriate frequency range. Therefore, the results verify that the model updating method using Where Error b and Error a is the relative error between the frequency values of the FE model before and after the updating and the target value, respectively.
As can be seen from Table 5, under the excitation of 10% random white noise, the maximum relative error of the frequency is changed from −5.49% to −1.10%, which verifies the effectiveness and anti-noise performance of the proposed method in this paper.
Where Error b and Error a is the relative error between the parameters of the FE model before and after the updating and the target value, respectively.
The model updating results of parameters are shown in Table 6, under the excitation of 10% random white noise, the maximum relative error of the updated parameters is changed from 30.00% to −0.24%, which greatly improves the reliability of the FE model and verifies the feasibility of the model updating method in this paper. Figures 11 and 12, due to modeling errors between the target model and the FE model, there is obviously a big gap between the FRFs of the FE model before updating and the target model under 5% random white noise. While after 20 iterative computations, the error between the obtained FRFs and the target value is quite close, which shows that the model updating method using FRFs based on Sherman-Morrison formula can effectively reduce noise interference. At the same time, the influence of damping on the amplitude of FRFs is effectively reduced by choosing the appropriate frequency range. Therefore, the results verify that the model updating method using FRFs based on Sherman-Morrison formula has a good anti-noise performance.

As shown in
anti-noise performance of the proposed method in this paper.
Where Errorb and Errora is the relative error between the parameters of the FE model before and after the updating and the target value, respectively.
The model updating results of parameters are shown in Table 6, under the excitation of 10% random white noise, the maximum relative error of the updated parameters is changed from 30.00% to −0.24%, which greatly improves the reliability of the FE model and verifies the feasibility of the model updating method in this paper.
As shown in Figures 11 and 12, due to modeling errors between the target model and the FE model, there is obviously a big gap between the FRFs of the FE model before updating and the target model under 5% random white noise. While after 20 iterative computations, the error between the obtained FRFs and the target value is quite close, which shows that the model updating method using FRFs based on Sherman-Morrison formula can effectively reduce noise interference. At the same time, the influence of damping on the amplitude of FRFs is effectively reduced by choosing the appropriate frequency range. Therefore, the results verify that the model updating method using FRFs based on Sherman-Morrison formula has a good anti-noise performance.  Figure 11. Comparison of FRFs data measured at Node 13 under 5% random white noise. Figure 11. Comparison of FRFs data measured at Node 13 under 5% random white noise.

Modal Test
To investigate the applicability of the proposed method in practical engineering, a simplified aluminum solar panel frame is designed at the dimensions shown in Figure 13. This only considers the frame components and ignores the influence of supporting structure and hinge. The connection between plate and plate is directly integrated.

Modal Test
To investigate the applicability of the proposed method in practical engineering, a simplified aluminum solar panel frame is designed at the dimensions shown in Figure 13. This only considers the frame components and ignores the influence of supporting structure and hinge. The connection between plate and plate is directly integrated.

Modal Test
To investigate the applicability of the proposed method in practical engineering, a simplified aluminum solar panel frame is designed at the dimensions shown in Figure 13. This only considers the frame components and ignores the influence of supporting structure and hinge. The connection between plate and plate is directly integrated. The piezoelectric accelerometers were adopted in this modal test, which with a mass of 20 g, charge sensitivity of 6.05 pC/ms 2 , and frequency range from 0.5 Hz to 5 kHz. The CL-YD-303 hammer was also adopted, the reference sensitivity is 3.99 pC/N. CRAS V7.0, a vibration and dynamic signal acquisition and analysis system developed by Nanjing AnZheng Software Engineering Company, was used as the analytical instrument, and MaCras was used as the modal analysis software.
The layout of test site is shown in Figure 14a. The end of the solar panel frame is fixed with a support, and the two accelerometers are installed on Position 1 and Position 2, where the mass of Accelerometer 1 is m1 = 20 g and the mass of Accelerometer 2 is m2 = 20 g. Then make a single point excitation at Position 3 with the hammer to measure the FRFs within 0~30 Hz of the solar panel frame. The FRFs obtained after 5 excitations are fitted by MaCras. The piezoelectric accelerometers were adopted in this modal test, which with a mass of 20 g, charge sensitivity of 6.05 pC/ms 2 , and frequency range from 0.5 Hz to 5 kHz. The CL-YD-303 hammer was also adopted, the reference sensitivity is 3.99 pC/N. CRAS V7.0, a vibration and dynamic signal acquisition and analysis system developed by Nanjing AnZheng Software Engineering Company, was used as the analytical instrument, and MaCras was used as the modal analysis software.
The layout of test site is shown in Figure 14a. The end of the solar panel frame is fixed with a support, and the two accelerometers are installed on Position 1 and Position 2, where the mass of Accelerometer 1 is m 1 = 20 g and the mass of Accelerometer 2 is m 2 = 20 g. Then make a single point excitation at Position 3 with the hammer to measure the FRFs within 0~30 Hz of the solar panel frame. The FRFs obtained after 5 excitations are fitted by MaCras. The FE model of the solar panel frame was modeled by the FE analysis software 'Nastran' as shown in Figure 14b, which consists of 86 shell elements and 148 nodes. The mass points were established at Node 1 and Node 2 to conduct the model updating and both the mass identification of the two accelerometers through mass perturbation of the mass points.

Parameter Updating
Take the FRFs h13 measured by Accelerometer 1 as an example. Combined with the rules in Section 2.3, the frequency of the model updating is selected as the frequency point near the resonance region, and to avoid the problem of amplitude nonlinearity, the frequency points of the third resonances are ignored as is shown in Figure 15.  The FE model of the solar panel frame was modeled by the FE analysis software 'Nastran' as shown in Figure 14b, which consists of 86 shell elements and 148 nodes. The mass points were established at Node 1 and Node 2 to conduct the model updating and both the mass identification of the two accelerometers through mass perturbation of the mass points.

Parameter Updating
Take the FRFs h 13 measured by Accelerometer 1 as an example. Combined with the rules in Section 2.3, the frequency of the model updating is selected as the frequency point near the resonance region, and to avoid the problem of amplitude nonlinearity, the frequency points of the third resonances are ignored as is shown in Figure 15.
shown in Figure 14b, which consists of 86 shell elements and 148 nodes. The mass points were established at Node 1 and Node 2 to conduct the model updating and both the mass identification of the two accelerometers through mass perturbation of the mass points.

Parameter Updating
Take the FRFs h13 measured by Accelerometer 1 as an example. Combined with the rules in Section 2.3, the frequency of the model updating is selected as the frequency point near the resonance region, and to avoid the problem of amplitude nonlinearity, the frequency points of the third resonances are ignored as is shown in Figure 15. Where Errorb and Errora is the relative error between the frequency values of the FE model before and after the updating and the experimental value, respectively.
As can be seen from Table 7, after several iterative computations, the maximum relative error of the frequency is changed from 9.13% to 2.99%. Although the error between the second order Where Error b and Error a is the relative error between the frequency values of the FE model before and after the updating and the experimental value, respectively.
As can be seen from Table 7, after several iterative computations, the maximum relative error of the frequency is changed from 9.13% to 2.99%. Although the error between the second order frequency and the experimental value increased slightly after updating, the average error of the first four orders decreased from 7.03% to 1.72%, which verifies the effectiveness of the proposed method in this paper. The identification results of accelerometer mass are shown in Table 8. Since the accelerometer is attached with wires, the identification results will be affected to some extent, so it can be considered that the identification results are within the ideal range, which greatly improves the reliability of the FE model and verifies the feasibility of the model updating method proposed in this paper. As shown in Figure 16, due to the additional mass of the accelerometers between the real experimental model and the FE model, there is obviously a big gap between the FRFs of the FE model before updating and the experimental structure. While after several iterative computations, the error between the obtained FRFs and the real experimental value h 13 is quite close, which shows that the proposed method has a good model updating capability and engineering application.
proposed method has a good model updating capability and engineering application.   Using modal assurance criteria (MAC) to verify the mutual independence and consistency of the modal shapes between the first 4 order of the FE model before and after updating and the corresponding experimental model: The ϕu, ϕv are the modal vector of the FE model and the experimental model, respectively, which are all the column vector, u, v is the modal order.
Theoretically, if the two models are the same physical modal, the MAC value is 1, if the two models are not the same physical modal or if the modal vectors are independent of each other, the MAC value is 0. The MAC value will be affected due to the nonlinearity of the structure, external noise interference of the measured data and improper data processing. It is generally believed that if MAC value is greater than 0.9, the two models are related, and if MAC value is less than 0.05, the two models are unrelated. Table 9 shows the comparison of modal shapes between the first 4 order of the FE model before and after updating and the corresponding experimental model. The first modal of the solar panel frame is first-order bending around X axis; the second modal is first-order torsion around Y axis; the third modal is second-order bending around X axis; the fourth modal is second-order torsion around Y axis. The blue line in the modal shape diagram represents the fixed side. Using modal assurance criteria (MAC) to verify the mutual independence and consistency of the modal shapes between the first 4 order of the FE model before and after updating and the corresponding experimental model: The φ u , φ v are the modal vector of the FE model and the experimental model, respectively, which are all the column vector, u, v is the modal order.
Theoretically, if the two models are the same physical modal, the MAC value is 1, if the two models are not the same physical modal or if the modal vectors are independent of each other, the MAC value is 0. The MAC value will be affected due to the nonlinearity of the structure, external noise interference of the measured data and improper data processing. It is generally believed that if MAC value is greater than 0.9, the two models are related, and if MAC value is less than 0.05, the two models are unrelated. Table 9 shows the comparison of modal shapes between the first 4 order of the FE model before and after updating and the corresponding experimental model. The first modal of the solar panel frame is first-order bending around X axis; the second modal is first-order torsion around Y axis; the third modal is second-order bending around X axis; the fourth modal is second-order torsion around Y axis. The blue line in the modal shape diagram represents the fixed side.  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids Appl. Sci. 2020, 10, x FOR PEER REVIEW 21 of 23 It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids Appl. Sci. 2020, 10, x FOR PEER REVIEW 21 of 23 It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids  It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids Appl. Sci. 2020, 10, x FOR PEER REVIEW 21 of 23 It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids Appl. Sci. 2020, 10, x FOR PEER REVIEW 21 of 23 It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method. It can be seen from Figure 18 that, after several iterative computations, the MAC value of the first 4 orders after updating are all improved, indicating that the FE modal shapes after updating are well matching with the experimental ones, which further verifies the accuracy requirement of the proposed method.

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids

Conclusions
A model updating method using FRFs based on Sherman-Morrison sensitivity is proposed in this paper. The sensitivity of FRFs is calculated using Sherman-Morrison formula which avoids recalculation of FE dynamic analysis. The simulation examples of truss and solar wing and the experiment of aluminum solar panel frame are implemented to verify the accuracy of the proposed method. The following conclusions can be drawn:

1.
Only the initial FRFs and parameter perturbations are employed for the acquisition of sensitivity, which avoids recalculation of FE dynamic analysis, and improves the efficiency of analysis and calculation at the same time. Only a small amount of FRFs data is needed to complete the model updating process.

2.
With the effect of 5% random white noise, the parameter error after updating is within ±1.00%, which shows that the proposed method has a certain anti-noise capability. The experiment study shows that the proposed method can effectively identify the structural parameters.

3.
Selecting appropriate frequency range and points can effectively improve the stability and efficiency of model updating. In addition, the method can consider multiple parameters simultaneously and has great potential for engineering application.