[0045] In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention.
[0046] The application principle of the present invention will be described in detail below with reference to the accompanying drawings.
[0047] Such as figure 1 As shown, the method for system model correction for multi-point excitation vibration test provided by the embodiment of the present invention includes the following steps:
[0048] S101: By establishing a finite element model for a specific engineering object, and obtaining the modal parameters of the finite element model through finite element analysis, matching the degree of freedom with the experimental test results through model condensation, so that the finite element calculation results and the experimental test results are in the degrees of freedom Last one-to-one correspondence;
[0049] S102: Then evaluate the accuracy of the initial finite element model through modal analysis. If the error between the calculation result of the finite element model and the experimental test result is within a reasonable range, the initial finite element model is a reliable model and no model modification is required If the modal analysis is carried out and it is found that the analysis result fails to meet the engineering requirements, it means that the initially established finite element model is unreliable and the initial finite element model needs to be revised.
[0050] The system model correction method for multi-point excitation vibration test provided by the embodiment of the present invention establishes a finite element model of a specific engineering object, and obtains the modal parameters of the finite element model through finite element analysis, and then passes the model with the experimental test results The polycondensation performs the degree of freedom matching, so that the finite element calculation results and the experimental test results have a one-to-one correspondence in the degrees of freedom; then through modal analysis, the accuracy of the initial finite element model is evaluated, if the finite element model calculation results are in error with the experimental test results Within a reasonable range, it means that the initial finite element model is a reliable model and does not need to be modified. If the analysis result is found to fail to meet the engineering requirements after modal analysis, it means that the initially established finite element model is unreliable and needs To revise the initial finite element model, the detailed steps are as follows:
[0051] Step one, establishment of finite element model
[0052] The multi-point excitation vibration test system model is divided into three parts: moving coil, fixture and aircraft. Before establishing the finite element model, it is necessary to establish a geometric model. Due to the complexity of the component model, it must be simplified as follows:
[0053] 1) In order to facilitate finite element analysis, delete features such as chamfers, rounded corners, and threaded holes that have little effect on the overall response, and fill some small round holes to facilitate mesh division and calculation;
[0054] 2) In order to ensure that the quality and center of mass of each component are consistent with the original component, the model is mass-equivalent according to the given mass and center of mass distribution table and the center of gravity coordinate formula in theoretical mechanics;
[0055] 3) Material parameters are equivalently treated according to actual properties.
[0056] Step two, research on model condensation
[0057] The model condensation method realizes the degree of freedom matching by reducing the system matrix scale of the finite element model. There are mainly two types of direct method and iterative method. The direct method is represented by Guyan polycondensation method and IRS polycondensation method; the iterative method has a stepwise approximation model Polycondensation method, improved Guyan-recursive reduction method, etc. The model condensation method needs to ensure that the low-order frequencies and mode shapes of the two finite element models before and after condensation are consistent.
[0058] In the static (Guyan) reduction, it is necessary to select a set of degrees of freedom (called a set), which are the degrees of freedom reserved in the analysis. The complement of this set (called the o set) is the set of degrees of freedom to be omitted in the analysis.
[0059] Step three, based on the sensitivity analysis model correction method
[0060] In the current field of model correction, sensitivity analysis-based model correction methods and parameter selection methods to be corrected occupy a dominant position. This type of method selects the parameters to be corrected through sensitivity analysis and uses sensitivity analysis techniques to correct the parameters. Sensitivity analysis is based on Taylor expansion and is limited by computational efficiency and difficulty. At present, first-order sensitivity analysis and second-order sensitivity analysis dominate. Due to the truncation error of Taylor expansion, the sensitivity analysis method is an approximate method. In the actual sensitivity calculation process, the finite difference method can be introduced to approximate the first-order or second-order sensitivity to reduce the calculation amount of the derivative operation and improve the calculation efficiency.
[0061] In the process of model modification, three elements of optimal design need to be determined: design variables, objective function and constraint conditions.
[0062] The parameters of design variables include physical parameters and geometric parameters. In the present invention, the physical parameters of the design variables include Young's modulus of elasticity, Poisson's ratio and density; the geometric parameters of the design variables are thickness. The purpose of the present invention is that the modal calculation result of the multi-point excitation test system is consistent with the modal test result, so the objective function is the natural frequency. In order to ensure that the total mass is consistent with the actual mass of the moving parts, the constraint condition is the total mass.
[0063] The application principle of the present invention will be further described below in conjunction with specific embodiments.
[0064] (1) Establishment of finite element model of multi-point excitation vibration test system
[0065] According to the principle of quality and centroid equivalence, the original model is simplified to a certain extent in Pro/E. After the model is simplified, the quality parameters of each part of the original aircraft model are provided for quality equivalence to ensure the quality and centroid of each part before and after the simplification. Consistent with the original model. Secondly, using Pro/E to analyze the quality of a certain part of the aircraft, the parameters such as volume, mass, and center of gravity can be obtained. After the simplification is completed, import the finite element software for meshing, define reasonable boundary conditions and material properties, and finally establish a finite element model.
[0066] (2) Model condensation
[0067] In fact, it is difficult to obtain the experimental mode shape information corresponding to all finite element nodes in the experiment. Therefore, the present invention assumes that in the simulation experiment, only the lateral degrees of freedom of some nodes can be accurately measured, and the remaining degrees of freedom are correspondingly reduced degrees of freedom. Therefore, the total number of main degrees of freedom accounts for a small part of the total number of finite element degrees of freedom, which greatly reduces the calculation time and improves the calculation efficiency.
[0068] In order to verify the correctness of the model polycondensation, a modal analysis is carried out with an aircraft as an example. In the present invention, the Guyan reduction method is used for reduction, and the designated points are used as the set of A (nodes to be retained), and then the Lanczos method is used to solve the first three natural frequencies. The following table shows the three-order natural frequencies before and after the modal condensation.
[0069] Table 1 Before and after modal condensation, the third-order natural frequency of the empty shaker
[0070] Mode order Before polycondensation (HZ) After polycondensation (HZ) error Level 1 48.547 48.65 0.2% Level 2 130.1 132.74 2% Level 3 245.47 256.89 4.6%
[0071] It can be seen from Table 1 that the first three-order natural frequency error before and after condensation is within 5%, so the number of degrees of freedom of the finite element model after condensation is basically the same as the number of degrees of freedom measured by the simulation experiment, but the calculation time is greatly reduced, so for large The finite element simulation of spacecraft is very necessary for model condensation.
[0072] (3) Model correction technology based on sensitivity analysis
[0073] For the real spacecraft structure, the density and the size of the structure are relatively easy to obtain, while the Young's modulus is difficult to obtain accurately due to the complex structure, welding, and articulation. Therefore, the present invention selects the Young's modulus and thickness of all elements as the parameters to be corrected, and performs model correction in combination with the modal polycondensation method. After the model is revised, the first few natural frequencies of the shaking table, the shaking table after installing the fixture, the aircraft, and the multi-point excitation vibration test system are shown in the following table.
[0074] Table 2 Comparison of the first two natural frequencies of the empty shaker table after model correction and the experimental data
[0075] Mode order Test value After model revision error Level 1 2139.5 2079.7 2.8% Level 2 2728.5 2619.3 4%
[0076] Table 3 Comparison of the first two natural frequencies of the shaker with the experimental data after the model is revised and the fixture is installed
[0077] Mode order Test value After model revision error Level 1 1720.5 1652.4 3.95% Level 2 2388.5 2491.8 4.3%
[0078] Table 4 Comparison of the first two natural frequencies of the aircraft after model modification and the experimental data
[0079] Mode order Test value After model revision error Level 1 51.18 49.81 2.6% Level 2 124.75 129.9 4%
[0080] Table 5 Comparison of the first two natural frequencies of the aircraft and the multi-point excitation vibration test system with experimental data after model modification
[0081] Mode order Test value After model revision error Level 1 43.87 44.818 2.16% Level 2 126.04 120.2 4.6%
[0082] It can be seen from the above table that the first three-order natural frequency error after model modification is within 5%, which meets the target requirement. Therefore, the model is reliable and can be used for further finite element simulation analysis.
[0083] The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. Any modification, equivalent replacement and improvement made within the spirit and principle of the present invention shall be included in the protection of the present invention. Within range.