Finite Element Model Updating of Satellite Sailboard Based on Sensitivity Analysis

'e modal analysis of a satellite sailboard finite element model is carried out to accurately investigate the response of a satellite sailboard in a complex loaded space environment through simulation.'e basic excitation vibration test of the satellite sailboard is used to perform model matching and a correlation test. Appropriate design variables are selected through sensitivity analysis. Modal analysis data and vibration table excitation test response data are used to modify the finite element model. After optimization, the orthogonality of the simulated vibration mode and experimental vibration mode is good.'e low-order frequency errors in the simulationmodel are less than 5%, the high-order errors are less than 10%, and the modal confidenceMAC values are above 0.8. 'e modal frequency and mode shape are closer to the experimental modal frequency and mode shape, respectively. 'e simulation and test acceleration response of the modified finite element model of a honeycomb panel are compared under the two conditions of sine sweep and random vibration. 'e acceleration response curves of reference points are consistent, and amplitude and frequency errors are within acceptable limits. 'e model updating effect is evident, which provides good reference for research on satellites and other aerospace products.


Introduction
At present, a large number of composite materials are used in spacecraft products [1, 2].Among them, honeycomb sandwich panels are widely used.
ey provide the advantages of high specific stiffness and specific strength and good heat insulation, vibration isolation, and impact resistance [3][4][5][6][7].Honeycomb sandwich panel structures comprise 80-90% of the shell and shell structure of spacecraft [8,9].e simulation accuracy of its finite element model has an important impact on the environmental adaptability of rocket launch conditions (acceleration overload, vibration, impact, random noise, etc.).
e dynamic parameters (modal matrix, coupling coefficient matrix, B tran /B rot , etc.,) of the sailboard assembly of satellites provide a solid foundation for the precise control of the on-orbit attitude of satellites, as it is an extended large flexible accessory.
A satellite is subjected to an extremely harsh vibration shock and noise environment during use and launch.In particular, when a spacecraft launches and reenters the atmosphere at high speeds, the resulting random high-frequency vibrations and noise vibrations within the fairing can produce a mean square response acceleration of up to 50 g for the satellite structure.erefore, it is necessary to carry out comprehensive dynamic analysis before launch to obtain accurate response prediction results for flight safety and good operation [10][11][12].
At present, the methods of satellite structural dynamics prediction are mainly the ground dynamics tests of structures, which are supplemented by finite element dynamics analysis.
e vibration test of ground shakers is indispensable in the development and design of various types of satellites.Moreover, it is an important criterion for evaluating the conformity and safety of products [13,14].However, ground tests do not fully simulate the various environments and states of a satellite during launch and on-orbit operation.In addition, the physical test period is long and expensive; this hinders the development of the satellite industry.e finite element simulation technology can conduct a comprehensive dynamic analysis of a product and obtain its response under various conditions.It is convenient to use, inexpensive, and has a short design period.Furthermore, it has become an important method of satellite product design and development and an important supplement to ground tests.
Finite element modelling consists of numerous theoretical assumptions, model simplification, distortion of connection conditions, and other factors.Additionally, satellite structures are extremely complex.Various components and auxiliary structures interact with each other, and different composite materials are widely used.For these reasons, it is quite difficult to establish an accurate and reasonable finite element model.
ere is error between finite element analysis results and experimental values.e dynamic response of models indicates that there is still a large gap between model results and actual situations. is strongly affects the reliability of finite element analysis.e design of satellite structures is bound to be highly conservative.e finite element model updating technology of satellite structures was proposed to solve this problem.An established finite element model was modified using a small amount of ground vibration test data.
en, the modified model was used for comprehensive dynamic response analysis to obtain reliable dynamic response prediction results.
In recent years, according to the method of resolving problems, a model updating method based on response surface analysis has been developed.is method transforms a finite element model updating problem from an inverse problem to a direct problem.Among numerous updating techniques, the design parametric model updating technique based on sensitivity analysis is the most mature.
e vibration mode correlation analysis between different models is an indispensable part of the model updating technology, and it can provide a basis for engineers to judge the quality of test analysis, the advantages and disadvantages of an initial finite element model, and sensor position optimization.For a certain order frequency, the degree of influence of each part of a structure is not the same. is is related to the mode distribution corresponding to the order frequency.If a certain part of a structure is in a region with large amplitude of a mode, then a slight change in this part will have a strong impact on the natural order frequency.On the contrary, if a part of a structure is close to the mode node, order frequency will not be significantly affected even if there is a large change in this part.In a project, it is expected that minor modifications to a structure will make its natural frequency meet requirements and will not affect the static strength and process requirements of the structure.In this study, the model of a solar panel is updated based on a correlation test and sensitivity analysis.A finite element model is obtained, which can accurately simulate the vibration response under complex aerospace conditions. is study provides a new idea for the precise finite element modelling of satellite solar panels and lays a strong foundation for the application research of satellite panels with realistic engineering significance.

Comparison of Modal Analysis and Modal
Test of Satellite Sailboard

Structure and Finite Element Modelling of Satellite
Sailboard.A regular hexagonal aluminum honeycomb sandwich panel with complex geometry is the most common form of honeycomb panel structures.Its geometry includes upper and lower panels, a central regular hexagonal honeycomb core, and a glue layer that connects the panels and a core layer.e glue layer is small in thickness and light in weight, and it is not the main bearing part.erefore, the influence of the glue layer can be neglected in the analysis of honeycomb sandwich panels [15,16].e structure of a honeycomb sandwich panel is shown in Figure 1.
In Figure 1, a is the board length of the honeycomb sandwich panel; b is the board width; h c is the thickness of the core; the thickness of the upper panel or lower panel is h f ; t is the wall thickness of the core; and L is the length of the side.
e finite element model of a sailboard is established using PATRAN and NASTRAN (MSC Software Corporation).e model consists of upper and lower panels, a honeycomb core layer, a frame, a support base, and a compression sleeve.e upper and lower panels, frame, and intermediate cell are shell units, and the remaining are solid units.Adhesive connections are used between the parts.e expanded diagram of the model is shown in Figure 2.
e unit type and material properties of each part of the sailboard finite element model are provided in Table 1.

Results of Modal Analysis and Test.
e constrained modal analysis of the established finite element model of the satellite sailboard is carried out.e modal frequency and mode shape of each order are obtained through modal analysis [17][18][19][20][21][22]. e simulation analysis of the mode shape of the finite element model of the satellite sailboard is shown in Figure 3. e LMS test system is used to carry out a modal test on the satellite sailboard.e four corners of the sailboard are fixed on a fixture, and the fixture is fixedly connected with a vibration platform through bolts.Sixteen 4524-B acceleration sensors are uniformly distributed on the panel of the sailboard.A hammer drives the centre of the panel.A sensor port is connected to a cable for signal output, and the other end of the cable is connected with an acquisition system.e acquisition system outputs signals to a PC for data processing.e test setup is shown in Figure 4.
e arrangement of measurement sensors must be considered so that the final measurement result can easily reflect the vibration of the overall mode shape.Hence, the sensors should be arranged symmetrically and evenly on the entire satellite sailboard.e sensor arrangement in this test is shown in Figure 5. e sensors are placed at the ends and intersections of lines.
e various mode shapes of the test are obtained, and the first two mode shapes are shown in Figure 6.

Model Matching and Correlation Test.
e correlation of the model includes frequency correlation and mode shape 2 Shock and Vibration correlation.Frequency correlation represents the degree of correlation between the modal test frequency (ω EMA ) and simulated analysis frequency (ω FEA ), and it is expressed as Ideally, ε should be close to zero.e larger the difference between ε and zero, the worse the correlation between ω EMA and ω FEA and vice versa.
Modal correlation indicates the degree of correlation between the mode shape of the modal test (ϕ EMA ) and that of the modal analysis (ϕ FEA ).Modal correlation is evaluated based on modal con dence (MAC), as shown in the following equation: e MAC value is always between 0 and 1.If the value on the diagonal in the matrix is greater than or equal to 0.8, the simulation analysis mode has a good correlation with the corresponding experimental mode.In addition, if the values of the remaining positions in the matrix are less than or equal to 0.3, it implies that the modes corresponding to di erent orders have better independence and are not affected by each other.e computational logic block diagram of the correlation test is shown in Figure 7.
e rst seven modal simulation results of the satellite sailboard are compared with the frequency and mode shapes of the test results.e obtained frequency values and mode shapes are inserted into the equations for frequency correlation (1) and mode shape correlation (2).e correlation of the models is calculated, and the results are shown in Table 2.
According to the correlation calculation results, the histogram for matching the MAC values of the simulation analysis mode and experimental mode is shown in Figure 8.

Sensitivity Analysis and Model Updating
3.1.Introduction to Optimization Algorithm.It has been proved that no single optimization algorithm can be applied to all structural optimization problems and a combined optimization strategy can achieve better results.
Here, the global search capability of an intelligent optimization algorithm (particle swarm optimization: PSO) [23,24] and the local search capability based on a sensitivity analysis gradient method (NLPQL) are used to achieve a balance between optimization quality and optimization efciency [25,26].An optimized Latin hypercube test is used to design sample points to prevent running the NASTRAN nite element program every time in the iterative process.e fourth-order polynomial response surface (RMS) model of structural response is constructed for further improving the optimization e ciency.

Shock and Vibration
e PSO algorithm simulates the migration and social behaviour in the foraging process of birds.First, a system initialises a group of random particles.en, it continuously adjusts their ight direction and speed through their memory and the exchange of group information.Finally, the group reaches the optimal position.e algorithm assumes that a particle swarm is composed of N particles, each of which has a d-dimensional search space.e mathematical expressions of the particle's velocity and position at time are given by x id (t + 1) x id (t) + v id (t + 1). (3) Here, i 1, 2, . . ., N is the number of particles, d 1, 2, . . ., D is the particle search dimension, v i (t) and x i (t) are the velocity and position of the i-th particle at time t, P i is the historical best position of the i-th particle search, P g is the optimal position of the entire population, c 1 and c 2 are particle acceleration coe cients, and ξ and ς are random numbers that are uniformly distributed from 0 to 1. e above equations are for the basic PSO algorithm, whose particle ight principle is shown in Figure 9.
e ow of the PSO algorithm is shown in Figure 10.

Sensitivity Analysis.
Errors are unavoidable in engineering calculations and tests.Error sources include the following factors: (1) Sti ness Error. is is the error caused by factors such as the thickness of a sheet, the modulus of elasticity of a material, and the sti ness of a joint between parts.
(2) Quality Error. is is the error caused by factors such as structural feature simpli cation, material density, and auxiliary quality such as screws.
(3) Systematic Error. is is the error caused by theoretical assumptions, structural discrepancies, iterative algorithms, signal acquisition, and the postprocessing of results during simulation and testing.
A sensitivity analysis method can be used to nd a set of design parameters that are sensitive to a target.Let the design parameter of the initial nite element model be P and eigenvector λ be its implicit function.When parameter p changes slightly, the rst-order Taylor expansion of the r-th feature quantity λ (r) is given by [27]    Shock and Vibration e sensitivity matrix is expressed as [28] S λ x i (t)

Shock and Vibration
In the above formula, n refers to the number of design parameters to be corrected, m is the number of feature quantities considered, Δλ is the residual vector, and Δp is the amount of change in the design parameter.
erefore, the sensitivity matrix of the MAC value of the vibration correlation coe cient can be expressed as Considering the structural vibrations of n degrees of freedom, the eigenvalue equation is given by where [K] and [M] are the real-symmetric structural sti ness matrix and mass matrix of order n × n, respectively, λ r is the characteristic quantity of the r-th order of the structure, and ϕ (r) is the r-th order mode shape of the structure.e formula of eigenvalue frequency sensitivity obtained using the orthogonality condition is given by e sensitivity matrix of the mode shape correlation coe cient is expressed as where ϕ e i is the i-th order test mode of vibration, ϕ a i is the calculated j-th order mode of vibration obtained by the through nite element analysis, and superscript T indicates matrix transposition.
e sensitivity calculation formula of the mode shape correlation coe cient of the i-th modi ed design parameter, p i , is as follows: where z ϕ a j /zp i is the partial derivative of the mode shape obtained by the j-th order nite element analysis model with respect to design parameter p i .
e sensitivity analysis of the variables in Table 3 is carried out by considering the modal frequency for the rst 7 orders and the correlation MAC value of the nite element model as response targets.As the order of magnitude varies signi cantly between variables, sensitivity also varies by di erent orders of magnitude, which is normalised for convenience in comparison [29,30].
As seen in Figures 11(a) and 11(b), the modulus, density, and thickness of the panel, the modulus and density of the honeycomb, and the modulus of the frame are sensitive variables.Hence, they can be selected as model updating parameters for optimization.6

Model Updating. Model updating is achieved by adjusting the sti ness and mass distribution of the nite
Shock and Vibration element model [31].e combined PSO + NLPQL iterative optimization strategy is used to nd the nal model parameter vector, and an error function is minimised with the least amount of model parameter modi cation to perform model updating [32].e design variables are the six parameter variables determined after sensitivity analysis.e constraint condition is that the total mass of the model must not be more than 110%.e objective function has the smallest error function, that is, the error in the simulation and test frequency is close to zero, and the correlation MAC value is close to one.
e optimal model updating parameters and the comparison of the parameters before and after the updating are shown in Table 4.
e frequency and mode shape convergence of the modi ed model are shown in Figure 12.
e modi ed nite element model is obtained after optimization.e modal analysis is repeated to obtain a new natural frequency and mode shape, which are compared with the experimental values and simulation result before the updating.
e corrected results are provided in Table 5.
e corrected MAC matrix and the comparison of the frequency mode shapes for the rst 7 orders are shown in Figures 13 and 14, respectively.
After optimization, the orthogonality of the simulated mode shape and the experimental mode shape is good.e low-order frequency error in the simulation model is less than 5%, the high-order error is less than 10%, and the modal con dence MAC is above 0.8.e modal frequency and mode shape of the optimized model are closer to the frequency and mode shape of the experimental mode.e model can more accurately re ect structural dynamic characteristics.6.

Model Updating Verification
e geometrical centroid position of the honeycomb panel is selected as the reference point, and the simulation and experimental acceleration responses of sine sweep and random vibration are compared.
e frequency response test setup and measurement point arrangement are shown in Figure 15.
e simulation and test results of the sine sweep resonance frequency and its response value are compared in each direction.
e comparison results are shown in Table 7.
e Z direction of the satellite sailboard (perpendicular to the plane of the board) is a key consideration for the sailboard under working conditions, and it is the most typical in a complex space environment.Four measurement points are arranged on the surface of the sailboard, i.e., at the xed position of the sailboard, the middle of the sailboard close to the long side, the centre of the plane of the sailboard structure, and in the vicinity of a re connection device.e simulation and experimental random vibration response curves of the satellite sailboard at the centre measurement point in the Z direction (perpendicular to the plane of the board) are shown in Figures 17(a e gure shows that the acceleration response curves at the reference point obtained through the test and simulation are consistent and that amplitude and frequency errors are within an acceptable range.

Conclusion
Structural sensitivity analysis (SOL200) is performed based on a nite element model developed using NAS-TRAN.
e model updating method can signi cantly improve the precision of the nite element model, which is of considerable engineering and practical signi cance.During the design, the objective function and constraint conditions can be determined according to the speci c requirements of structural performance.By selecting appropriate design variables through sensitivity analysis, satisfactory structures can be obtained through optimization.
Compared with the single gradient method and intelligent optimization algorithm, the combined optimization strategy exhibits higher reliability and better universality.
is can provide balance between optimization e ciency and optimization quality and ensure the global optimization of a solution as much as possible.e RMS method is used to establish an approximate model of structural response in the implementation of the combined optimization strategy, which can replace the actual nite element calculation model and considerably reduce the number of structural heavy analysis.Based on the proposed model updating method, the nite element model of a certain satellite panel structure is modi ed using the modal test data and excitation response data of a vibration table.It is veri ed that the simulation accuracy of the model is signi cantly improved.

Data Availability
e data used in this study may be made available upon request to the authors.

Figure 3 :Figure 4 :
Figure 3: Simulation analysis of the mode shape.(a) First mode shape.(b) Second mode shape.

Figure 5 :
Figure 5: Sensor placement in the test. 4

Figure 6 :Figure 7 :
Figure 6: Mode shapes of the test.(a) First mode shape.(b) Second mode shape.

Figure 8 :
Figure 8: MAC value matching of simulation mode and test mode.

c 1
ξP i (t)x i (t)

Figure 10 :
Figure 10: Operational logic diagram of the combined optimization algorithm.
4.1.Frequency Response Test Veri cation.Vibration response analysis is performed based on the modi ed nite element model to further verify the accuracy of the modi ed

Figure 11 :
Figure 11: Sensitivity of the design parameters.(a) Sensitivity of modal frequency to design variables.(b) Sensitivity of modal shape MAC values to design variables.

Figure 12 :
Figure 12: Modal convergence of the model during the updating process.(a) Convergence of frequency error during the updating process.(b) Convergence of MAC values during the updating process.

Figure 14 :Figure 13 :
Figure 14: Comparison of the modi ed 7th order simulation mode and test mode shapes for the rst 7 orders.

Figure 16 :
Figure 16: Frequency response curves at reference points in the Z direction.(a) Simulation curves.(b) Experimental curves.

Figure 15 :
Figure 15: Frequency response test setup.(a) Measurement point arrangement in X and Y directions.(b) Measurement point arrangement in Z direction.

Table 1 :
Unit type and material properties of each part of the sailboard nite element model.

Table 2 :
Calculation results of model correlation.

Table 3 :
Initial values of the elastic parameters and dimensions of the model.

Table 5 :
Frequency error and MAC value of the model before and after updating.

Table 4 :
Comparison of initial and optimized values of model elastic parameters and dimensions.

Table 7 :
Comparison of sine sweep resonance frequency and its response value.

Table 8 :
Loading conditions for random vibration excitation.

Table 9 :
Comparison of sinusoidal sine sweep resonance frequency and its response value.