CN104833466A - Spacecraft ground test and on-orbit micro-vibration mechanical environment mapping method - Google Patents
Spacecraft ground test and on-orbit micro-vibration mechanical environment mapping method Download PDFInfo
- Publication number
- CN104833466A CN104833466A CN201510219792.6A CN201510219792A CN104833466A CN 104833466 A CN104833466 A CN 104833466A CN 201510219792 A CN201510219792 A CN 201510219792A CN 104833466 A CN104833466 A CN 104833466A
- Authority
- CN
- China
- Prior art keywords
- orbit
- spacecraft
- finite element
- ground test
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种航天器地面测试与在轨微振动力学环境映射方法,该方法首先建立地面测试航天器有限元模型以及在轨航天器有限元模型;并进行模态分析后提取频率、振型数据,确定两个模型的频率、振型的一一对应关系;建立和确定减缩模型的正确性;最后通过BP网络实现从地面微振动测试状态的频率、振型到在轨状态的频率和振型之间的映射;根据映射得到的在轨模型的频率、振型,进行在轨模型的动力学响应分析。本发明消除了地面微振动测试状态下空气、重力、悬挂约束对地面微振动测试的影响,实现地面测试对实际在轨状态微振动特性的预示。同时该方法可以实现地面微振动测试数据和在轨测试数据的相互比较,验证地面微振动测试的有效性。
The invention discloses a method for mapping a spacecraft ground test and on-orbit microvibration dynamics environment. The method firstly establishes a finite element model of the ground test spacecraft and an on-orbit spacecraft finite element model; type data to determine the one-to-one correspondence between the frequency and mode shape of the two models; establish and confirm the correctness of the reduction model; Mapping between mode shapes; analyze the dynamic response of the on-orbit model according to the frequency and mode shape of the on-orbit model obtained from the mapping. The invention eliminates the influence of air, gravity and suspension constraints on the ground micro-vibration test in the ground micro-vibration test state, and realizes the prediction of the micro-vibration characteristics of the actual on-orbit state by the ground test. At the same time, this method can realize the mutual comparison between the ground micro-vibration test data and the on-orbit test data, and verify the validity of the ground micro-vibration test.
Description
技术领域technical field
本发明是一种航天器地面测试与在轨微振动力学环境方法,通过此方法,实现地面测试状态对实际在轨状态微振动特性的预示。The invention relates to a space vehicle ground test and on-orbit micro-vibration dynamics environment method, through which the prediction of the micro-vibration characteristics of the ground test state to the actual on-orbit state is realized.
背景技术Background technique
随着社会经济的发展,高分辨率航天器无疑是航天器发展的方向,如美国的KH系列军事观察卫星,从KH-1到KH-13其分辨率从12m提高到0.05m。深空探测遥感航天器与对地观测卫星相比,其分辨率要高出1~2个数量级,如哈勃空间望远镜(0.1角秒,1990年)。下一代空间望远镜詹姆斯·韦伯太空望远镜分辨率达0.004角秒。With the development of society and economy, high-resolution spacecraft is undoubtedly the direction of spacecraft development, such as the KH series of military observation satellites in the United States, from KH-1 to KH-13, its resolution has increased from 12m to 0.05m. Compared with earth observation satellites, the resolution of remote sensing spacecraft for deep space exploration is 1 to 2 orders of magnitude higher, such as the Hubble Space Telescope (0.1 arcsecond, 1990). The next-generation space telescope, the James Webb Space Telescope, has a resolution of 0.004 arcseconds.
微振动指航天器在轨运行期间,星上转动部件高速转动、大型可控构件驱动机构步进运动、变轨调姿期间推力器点火工作、大型柔性结构进出阴影冷热交变诱发扰动等都会使星体产生一种幅值较小、频率较高的抖动响应。大多数航天器都存在微振动扰动源。由于微振动力学环境效应幅值小、频率高,对大部分航天器不会产生明显影响,通常予以忽略。但对高精度航天器将严重影响有效载荷指向精度、稳定度及分辨率等重要性能指标,所以在高分辨率航天器设计中必须考虑微振动的影响。Micro-vibration refers to the high-speed rotation of the rotating parts on the star, the stepping motion of the large-scale controllable component drive mechanism, the ignition of the thruster during the orbit change and attitude adjustment, and the disturbance induced by the alternating cold and heat of the large-scale flexible structure when it enters and exits the shadow. Make the star produce a jitter response with a smaller amplitude and higher frequency. Micro-vibration disturbance sources exist in most spacecraft. Due to the small amplitude and high frequency of the micro-vibration dynamics environmental effect, it will not have a significant impact on most spacecraft and is usually ignored. However, for high-precision spacecraft, it will seriously affect important performance indicators such as payload pointing accuracy, stability and resolution, so the influence of micro-vibration must be considered in the design of high-resolution spacecraft.
由于空间飞行器在轨工作时所处的动力学环境极其复杂,加之在轨测试的成本高,且姿态控制系统对微振动响应无法测控,因此目前对航天器微振动的研究主要采用数值模拟和地面微振动测试两种方法。根据国外公开的文献开展调研工作,目前各国规模较大的且技术较成熟的地面微振动测试平台主要有Honeywell公司的SCT地面微振动测试台,JPL实验室的MPI地面微振动测试台以及SSL实验室的OT地面微振动测试台。然而地面测试和在轨航天器力学环境仍然存在很大的差异,地面微振动测试环境中的重力场、空气、约束(悬挂装置)等因素可能会使地面测试结果与在轨航天器微振动特性出现较大差别。因此,地面微振动测试结果仅仅能用于评估,并不能准确分析航天器在轨微振动特性。Due to the extremely complex dynamic environment in which the spacecraft works in orbit, the high cost of on-orbit testing, and the inability to measure and control the response of the attitude control system to micro-vibrations, the current research on spacecraft micro-vibrations mainly uses numerical simulation and ground simulation. There are two methods of micro-vibration testing. According to the research work carried out in foreign literature, the currently large-scale and mature ground micro-vibration test platforms in various countries mainly include Honeywell's SCT ground micro-vibration test bench, JPL laboratory's MPI ground micro-vibration test bench and SSL experiment. The OT ground micro-vibration test bench in the laboratory. However, there are still great differences between the ground test and the on-orbit spacecraft mechanical environment. Factors such as gravity field, air, constraints (suspension devices) and other factors in the ground micro-vibration test environment may make the ground test results different from the micro-vibration characteristics of the on-orbit spacecraft. There is a big difference. Therefore, the ground micro-vibration test results can only be used for evaluation, and cannot accurately analyze the on-orbit micro-vibration characteristics of spacecraft.
为了得到在轨航天器的微振动特性,而航天器结构复杂,难以得到航天器微振动的解析解,因此目前主要采用数值模拟的方法,美国等科研机构对此开展了大量的研究。MIT空间系统实验室开发了微振动集成建模与综合评价分析软件DOCS;NASA开发了能够进行抖振和结构/热/光学分析系统IME。虽然目前数值模拟能够在一定程度上得到航天器的微振动特性,但存在计算效率差和应用范围窄等问题。In order to obtain the micro-vibration characteristics of the spacecraft in orbit, and the structure of the spacecraft is complex, it is difficult to obtain the analytical solution of the micro-vibration of the spacecraft. Therefore, the method of numerical simulation is mainly used at present, and scientific research institutions such as the United States have carried out a lot of research on this. The MIT Space Systems Laboratory has developed a micro-vibration integrated modeling and comprehensive evaluation analysis software DOCS; NASA has developed a chattering and structural/thermal/optical analysis system IME. Although the current numerical simulation can obtain the micro-vibration characteristics of the spacecraft to a certain extent, there are problems such as poor calculation efficiency and narrow application range.
发明内容Contents of the invention
本发明提供一种航天器地面测试与在轨微振动力学环境映射方法,所述方法消除地面测试状态下空气、重力、约束对微振动特性的影响,实现地面测试对实际在轨状态微振动特性的预示。The invention provides a space vehicle ground test and on-orbit micro-vibration dynamics environment mapping method. The method eliminates the influence of air, gravity and constraints on the micro-vibration characteristics under the ground test state, and realizes the ground test on the micro-vibration characteristics of the actual on-orbit state. foreshadowing.
本发明提供的映射方法包括以下步骤:The mapping method provided by the present invention comprises the following steps:
(1)建立模拟地面微振动测试力学环境的地面测试航天器有限元模型以及在轨航天器有限元模型;(1) Establish the finite element model of the ground test spacecraft and the finite element model of the in-orbit spacecraft to simulate the mechanical environment of the ground micro-vibration test;
(2)通过对地面测试航天器有限元模型和在轨航天器有限元模型分别进行模态分析,提取频率、振型数据,分别建立地面测试模态坐标减缩模型和在轨模态坐标减缩模型,并确定地面测试航天器有限元模型和在轨航天器有限元模型的频率、振型的一一对应关系。(2) Through the modal analysis of the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft, the frequency and mode shape data are extracted, and the ground test modal coordinate reduction model and the on-orbit modal coordinate reduction model are respectively established , and determine the one-to-one correspondence between the frequency and mode shape of the finite element model of the ground test spacecraft and the finite element model of the in-orbit spacecraft.
(3)根据地面测试模态坐标减缩模型和地面测试航天器有限元模型的相应对比,确定地面测试模态坐标减缩模型的正确性;根据在轨模态坐标减缩模型和在轨航天器有限元模型的响应对比,确定在轨模态坐标减缩模型的正确性;(3) According to the corresponding comparison between the ground test modal coordinate reduction model and the ground test spacecraft finite element model, determine the correctness of the ground test modal coordinate reduction model; according to the on-orbit modal coordinate reduction model and the on-orbit spacecraft finite element model Comparing the response of the model to determine the correctness of the on-orbit modal coordinate reduction model;
(4)通过BP(Back Propagation)网络实现考虑重力、约束、空气力学环境因素的从地面测试状态的频率、振型到在轨状态的频率和振型之间的映射,得到在轨航天器有限元模型的频率、振型。(4) Through the BP (Back Propagation) network, the mapping between the frequency and mode shape of the ground test state to the frequency and mode shape of the on-orbit state is realized considering gravity, constraints, and aerodynamic environmental factors, and the limited The frequency and mode shape of the meta-model.
(5)根据映射得到的在轨航天器有限元模型的频率、振型,进行在轨航天器有限元模型的动力学响应分析。(5) According to the frequency and mode shape of the finite element model of the on-orbit spacecraft obtained by mapping, the dynamic response analysis of the finite element model of the on-orbit spacecraft is performed.
本发明提供的映射方法优点在于:The mapping method provided by the present invention has the advantages of:
实现了航天器地面测试与在轨微振动力学环境映射,消除了地面微振动测试状态下空气、重力、约束对地面微振动测试的影响,实现地面测试对实际在轨状态微振动特性的预示。同时该映射方法可以实现地面微振动测试数据和在轨测试数据的相互比较,验证地面微振动测试的有效性。Realized the mapping of spacecraft ground test and on-orbit micro-vibration dynamics environment, eliminated the influence of air, gravity, and constraints on ground micro-vibration test under the ground micro-vibration test state, and realized the prediction of the micro-vibration characteristics of the actual on-orbit state by ground test. At the same time, the mapping method can realize the mutual comparison between the ground micro-vibration test data and the on-orbit test data, and verify the validity of the ground micro-vibration test.
附图说明Description of drawings
图1是本发明映射方法的流程图;Fig. 1 is the flowchart of mapping method of the present invention;
图2是本发明航天器在轨和地面测试力学环境比较;Fig. 2 is the comparison of the mechanical environment of the spacecraft in orbit and the ground test of the present invention;
图3是本发明三级映射的流程图;Fig. 3 is the flowchart of three-level mapping of the present invention;
图4是本发明实例地面模态坐标减缩模型与有限元模型的时域响应对比;Fig. 4 is the time domain response comparison of the ground modal coordinate reduction model of the example of the present invention and the finite element model;
图5是本发明实例地面微振动时域响应相对误差随时间变化图;Fig. 5 is the time-domain response relative error of the ground micro-vibration of the example of the present invention, which varies with time;
图6是本发明实例在轨模态坐标减缩模型与有限元模型的微振动时域响应对比;Fig. 6 is a comparison of micro-vibration time-domain responses between the on-orbit modal coordinate reduction model and the finite element model of the example of the present invention;
图7是本发明实例在轨微振动时域响应相对误差随时间变化图;Fig. 7 is a time domain response relative error graph of micro-vibration in orbit in the example of the present invention;
图8是本发明实例预测在轨微振动时域响应;Fig. 8 is the time domain response of on-orbit micro-vibration predicted by the example of the present invention;
图9是本发明实例预测在轨响应与减缩模型响应相比相对误差随时间变化图。Fig. 9 is a time-dependent diagram of the relative error of the predicted on-orbit response compared with the reduced model response according to the example of the present invention.
具体实施方式Detailed ways
下面结合附图和实施例对本发明进行详细说明。The present invention will be described in detail below in conjunction with the accompanying drawings and embodiments.
本发明提供一种航天器地面测试与在轨微振动力学环境映射方法,如图1所示流程,所述映射方法包括如下步骤:The present invention provides a method for mapping spacecraft ground testing and on-orbit microvibration dynamics environment, as shown in the flow chart in Figure 1. The mapping method includes the following steps:
(1)建立航天器有限元模型,包括地面测试航天器有限元模型和在轨航天器有限元模型;(1) Establish the finite element model of the spacecraft, including the finite element model of the ground test spacecraft and the finite element model of the spacecraft in orbit;
建立在轨航天器有限元模型:根据给定典型航天器结构参数,进行在轨航天器有限元模型建立。建模完成后通过进行参数型模型修正如材料刚度、弹簧刚度以及模态阻尼比等,使调整后的数学模型尽可能的全面反映航天器结构的动力学特性(如频响函数、固有频率等)。Establish the finite element model of the on-orbit spacecraft: According to the given typical spacecraft structure parameters, the finite element model of the on-orbit spacecraft is established. After the modeling is completed, parametric model corrections such as material stiffness, spring stiffness, and modal damping ratio are performed to make the adjusted mathematical model fully reflect the dynamic characteristics of the spacecraft structure (such as frequency response function, natural frequency, etc. ).
建立地面测试航天器有限元模型:由于在进行地面测试时,需要对卫星增加边界条件模拟,如进行吊绳悬挂,同时会受到重力、空气的影响,因此在参数型模型修正完成后的在轨航天器有限元模型基础上,采用梁单元模拟建立吊绳,重力的影响可通过增加预应力来模拟,空气的影响可通过附加质量来模拟。Establish the finite element model of the ground test spacecraft: Since it is necessary to simulate the boundary conditions of the satellite during the ground test, such as hanging by a sling, it will be affected by gravity and air at the same time, so after the parametric model correction is completed, the in-orbit Based on the finite element model of the spacecraft, the beam element is used to simulate the establishment of the suspension rope. The influence of gravity can be simulated by adding prestress, and the influence of air can be simulated by adding mass.
(2)建立地面测试模态坐标减缩模型和在轨模态坐标减缩模型,确定地面测试航天器有限元模型和在轨航天器有限元模型的频率、振型的一一对应关系。(2) Establish the ground test modal coordinate reduction model and the on-orbit modal coordinate reduction model, and determine the one-to-one correspondence between the frequency and mode shape of the ground test spacecraft finite element model and the on-orbit spacecraft finite element model.
多自由度阻尼系统满足下面方程:The multi-degree-of-freedom damping system satisfies the following equation:
其中,M、K、C分别为质量矩阵、刚度矩阵、阻尼矩阵,u(t)分别表示t时刻的加速度、速度、位移,u0、分别表示初始位移、初始速度,f(t)为t时刻的外力。因此,地面测试航天器有限元模型与在轨航天器有限元模型之间的映射关系,可以看成两模型质量矩阵、刚度矩阵、阻尼矩阵之间的映射关系。以固有振型矩阵Φ引入坐标变换u(t)=Φq(t),q(t)代表模态坐标,则方程(1)转换为:Among them, M, K, and C are mass matrix, stiffness matrix, and damping matrix, respectively, u(t) represent the acceleration, velocity and displacement at time t respectively, u 0 , represent the initial displacement and initial velocity respectively, and f(t) is the external force at time t. Therefore, the mapping relationship between the finite element model of the ground test spacecraft and the finite element model of the in-orbit spacecraft can be regarded as the mapping relationship between the mass matrix, stiffness matrix and damping matrix of the two models. The coordinate transformation u(t)=Φq(t) is introduced by the natural mode matrix Φ, and q(t) represents the modal coordinates, then the equation (1) is transformed into:
将选取的固有振型矩阵关于模态质量归一化,则有,Normalize the selected natural mode shape matrix with respect to the modal mass, then,
Mq=ΦTMΦ=I,Kq=ΦTKΦ=diag[ω2],Cq=ΦTCΦ (3)M q =Φ T MΦ=I, K q =Φ T KΦ=diag[ω 2 ], C q =Φ T CΦ (3)
Mq、Cq、Kq分别为模态坐标下的质量矩阵、阻尼矩阵、刚度矩阵,I表示单位阵,diag[ω2]表示以各阶频率ω平方为对角线的对角阵,将阻尼近似处理为比例阻尼,此时方程(2)转换为n个单自由度阻尼系统:M q , C q , and K q are the mass matrix, damping matrix, and stiffness matrix under the modal coordinates, respectively, I represents the unit matrix, and diag[ω 2 ] represents the diagonal matrix whose diagonal is the square of frequency ω of each order, The damping is approximately treated as proportional damping, and equation (2) is transformed into n single-degree-of-freedom damping systems:
qj、ζj、ωj分别表示第j个模态坐标的广义位移、广义速度、广义加速度、阻尼比、频率,将方程组(4)按照频率从小到大排列,取前m(1<m<n)阶对应的方程组,引入状态向量X,q j , ζ j and ω j represent the generalized displacement, generalized velocity, generalized acceleration, damping ratio, and frequency of the j-th modal coordinate respectively. Arrange the equations (4) in ascending order of frequency, and take the first m (1<m<n ) order corresponding to the system of equations, introducing the state vector X,
X=[x11,x22,…,x1j,x2j,…,x1m,x2m] (5)X=[x 11 ,x 22 ,…,x 1j ,x 2j ,…,x 1m ,x 2m ] (5)
其中则式(4)可以改写成以矢量形式表达的含有2m个一阶微分方程的微分系统:in Then equation (4) can be rewritten as a differential system expressed in vector form containing 2m first-order differential equations:
X'=f(X) (6)X'=f(X) (6)
通过Runge-Kutta数值积分求解,得到系统在物理坐标下的位移为:Solved by Runge-Kutta numerical integration, the displacement of the system in the physical coordinates is obtained as:
由以上推导过程可知,建立地面测试航天器有限元模型与在轨航天器有限元模型之间的映射关系,等效于建立地面测试航天器有限元模型与在轨航天器有限元模型之间特征值、特征向量的映射关系。针对地面测试航天器有限元模型和在轨航天器有限元模型,利用Nastran二次开发语言DMAP,分别提取地面测试航天器有限元模型和在轨航天器有限元模型的频率、振型数据,分别建立各自的地面测试模态坐标减缩模型和在轨模态坐标减缩模型。对地面测试航天器有限元模型和在轨航天器有限元模型分别进行模态分析,其中地面测试航天器有限元模型采用预应力模态分析,对比确定地面测试航天器有限元模型和在轨航天器有限元模型的频率、振型存在一一对应关系。From the above derivation process, it can be seen that establishing the mapping relationship between the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft is equivalent to establishing the characteristic relationship between the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft The mapping relationship between value and feature vector. For the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft, the Nastran secondary development language DMAP is used to extract the frequency and mode shape data of the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft, respectively. Establish their respective ground test modal coordinate reduction models and on-orbit modal coordinate reduction models. The finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft are respectively subjected to modal analysis. There is a one-to-one correspondence between the frequency and mode shape of the finite element model of the device.
(3)响应计算,确定地面测试模态坐标减缩模型和在轨模态坐标减缩模型的正确性。(3) Response calculation to determine the correctness of the ground test modal coordinate reduction model and the on-orbit modal coordinate reduction model.
对地面测试模态坐标减缩模型和地面测试航天器有限元模型进行时域响应分析,通过分析比较地面测试模态坐标减缩模型和地面测试航天器有限元模型的时域响应;对在轨模态坐标减缩模型和在轨航天器有限元模型进行时域响应分析,通过分析比较在轨模态坐标减缩模型和在轨航天器有限元模型的时域响应;如果时域响应结果显示响应分析误差低于20%,则表明所述的地面测试模态坐标减缩模型或在轨模态坐标减缩模型是正确的。The time domain response analysis of the ground test modal coordinate reduction model and the ground test spacecraft finite element model is carried out, and the time domain response of the ground test modal coordinate reduction model and the ground test spacecraft finite element model is analyzed and compared; the on-orbit modal Coordinate reduction model and on-orbit spacecraft finite element model for time-domain response analysis, and compare the time-domain responses of on-orbit modal coordinate reduction model and on-orbit spacecraft finite element model through analysis; if the time-domain response results show that the response analysis error is low If it is less than 20%, it indicates that the ground test modal coordinate reduction model or the on-orbit modal coordinate reduction model is correct.
(4)航天器微振动力学映射关系的建立。(4) Establishment of spacecraft fretting dynamics mapping relationship.
在轨航天器与地面测试力学环境具有较大区别,如图2所示,在轨航天器的工作环境为真空、无重力的自由飞行,地面测试中的航天器受到重力、约束(是指地面测试航天器有限元模型中采用的吊绳悬挂约束)、空气环境噪声因素的影响,需考虑这些因素对频率和振型的影响量级,通过采用BP神经网络建立映射关系。通过BP神经网络实现从地面测试航天器的频率、振型到在轨航天器的频率和振型之间的映射。The mechanical environment of the on-orbit spacecraft is quite different from that of the ground test. As shown in Figure 2, the working environment of the on-orbit spacecraft is vacuum and free flight without gravity. The spacecraft in the ground test is subject to gravity, constraints (referring to the ground To test the influence of the sling suspension constraint used in the spacecraft finite element model) and the air environment noise factors, it is necessary to consider the magnitude of the influence of these factors on the frequency and mode shape, and establish a mapping relationship by using the BP neural network. The mapping from the frequency and mode shape of the ground test spacecraft to the frequency and mode shape of the orbiting spacecraft is realized through the BP neural network.
将复杂力学环境映射分解成为三级映射,如图3所示。首先根据地面测试得到的频率、振型数据(通过仅考虑空气影响的映射关系)得到无空气状态(即包含重力、约束影响)的频率、振型数据。同理,通过仅考虑重力的映射关系,得到无空气无重力(包含约束状态)的频率、振型数据。最后根据仅考虑约束的映射关系得到无约束的在轨的频率、振型数据。三级映射具体实现步骤如下:The complex mechanical environment mapping is decomposed into three-level mapping, as shown in Figure 3. Firstly, the frequency and mode shape data of the airless state (that is, including gravity and constraint effects) are obtained according to the frequency and mode shape data obtained from the ground test (by considering only the mapping relationship of air influence). In the same way, by only considering the mapping relationship of gravity, the frequency and mode shape data of no air and no gravity (including the constraint state) are obtained. Finally, the unconstrained on-orbit frequency and vibration mode data are obtained according to the mapping relationship that only considers constraints. The specific implementation steps of the three-level mapping are as follows:
a)考虑空气影响的映射关系研究;a) Research on the mapping relationship considering the influence of air;
运用BP神经网络方法建立考虑空气影响的从地面测试到无空气状态固有频率和振型的映射关系。不考虑重力、约束的变化,仅考虑空气的影响,计算得到不同空气密度情况下的频率、振型数据,根据这些数据建立映射关系,通过此映射关系得到在轨(无空气,即空气密度为0的状态)状态下数据。The BP neural network method is used to establish the mapping relationship from the ground test to the natural frequency and mode shape of the airless state considering the influence of air. Regardless of the changes in gravity and constraints, only the influence of air is considered, and the frequency and mode shape data under different air densities are calculated, and the mapping relationship is established based on these data. Through this mapping relationship, the on-orbit (no air, that is, the air density is 0 state) state data.
b)考虑重力影响的映射关系研究;b) Research on the mapping relationship considering the influence of gravity;
运用神经网络模型建立考虑重力影响的从地面测试到在轨状态固有频率和振型的映射关系。在考虑了空气的基础上,考虑重力的影响,计算得到不同重力加速度情况下的频率、振型,根据这些数据建立映射关系,通过此映射关系得到在轨(0g)状态下数据。The neural network model is used to establish the mapping relationship from the ground test to the natural frequency and mode shape of the on-orbit state considering the influence of gravity. On the basis of considering the air and the influence of gravity, the frequency and mode shape under different gravitational accelerations are calculated, and the mapping relationship is established based on these data, and the data in the orbit (0g) state is obtained through this mapping relationship.
c)考虑约束影响的映射关系研究;c) Research on the mapping relationship considering the influence of constraints;
运用神经网络模型建立考虑约束影响的从地面测试到在轨状态固有频率和振型的映射关系。在考虑了空气、重力的基础上,考虑约束的影响,计算得到不同吊绳截面积情况下的频率、振型数据。根据这些数据建立映射关系,由此映射关系得到在轨(无约束,等效于吊绳截面积为0mm2的状态)状态下数据。The neural network model is used to establish the mapping relationship from the ground test to the natural frequency and mode shape of the on-orbit state considering the influence of constraints. On the basis of air and gravity, and considering the influence of constraints, the frequency and mode data of different suspension rope cross-sectional areas are calculated. Based on these data, a mapping relationship is established, from which the data in the state of on-orbit (unconstrained, equivalent to the state where the cross-sectional area of the suspension rope is 0 mm 2 ) is obtained.
(5)在轨响应计算。根据映射得到的在轨航天器有限元模型的频率、振型,进行在轨航天器有限元模型的动力学响应计算。可以对比实际在轨数据,验证此映射方法和地面测试的有效性。(5) On-orbit response calculation. According to the frequency and mode shape of the finite element model of the on-orbit spacecraft obtained through mapping, the dynamic response calculation of the finite element model of the on-orbit spacecraft is performed. The effectiveness of this mapping method and ground testing can be verified by comparing with actual in-orbit data.
对本发明实例地面测试航天器有限元模型进行预应力模态分析,第七到十阶频率为140.904Hz、146.906Hz、149.918Hz、160.582Hz。在轨状态前四阶频率为140.852Hz、146.830Hz、149.841Hz、160.504Hz。同时振型也都一一对应,因此可以对地面测试航天器有限元模型和在轨航天器有限元模型的频率和振型建立映射关系。The prestressed modal analysis is performed on the finite element model of the ground test spacecraft of the present invention, and the seventh to tenth order frequencies are 140.904Hz, 146.906Hz, 149.918Hz, and 160.582Hz. The first four frequencies in orbit state are 140.852Hz, 146.830Hz, 149.841Hz, 160.504Hz. At the same time, the mode shapes are also in one-to-one correspondence, so a mapping relationship can be established between the frequencies and mode shapes of the ground test spacecraft finite element model and the in-orbit spacecraft finite element model.
对地面测试航天器有限元模型和在轨航天器有限元模型进行响应分析,分别取前20阶弹性模态对地面测试航天器有限元模型和在轨航天器有限元模型进行减缩,地面测试模态坐标减缩模型与地面测试航天器有限元模型的时域响应如图4所示,地面测试时域响应相对误差随时间变化如图5所示,图6是在轨模态坐标减缩模型与在轨航天器有限元模型的微振动时域响应对比,图7是在轨微振动时域响应相对误差随时间变化图。由图5、图7可知相对误差小于20%,因此对地面测试和在轨航天器有限元模型的减缩是正确的。Response analysis was performed on the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft, and the first 20 order elastic modes were respectively taken to reduce the finite element model of the ground test spacecraft and the finite element model of the on-orbit spacecraft. The time-domain response of the modal coordinate reduction model and the ground test spacecraft finite element model is shown in Figure 4, and the relative error of the ground test time-domain response changes with time is shown in Figure 5. The comparison of the micro-vibration time-domain response of the finite element model of the orbiting spacecraft. Figure 7 is a graph of the relative error of the on-orbit micro-vibration time-domain response versus time. It can be seen from Fig. 5 and Fig. 7 that the relative error is less than 20%, so the reduction of the finite element model of the ground test and the spacecraft in orbit is correct.
表1为频率映射计算结果,空气这两列表示的是仅考虑空气影响情况下的在轨与地面振动测试相比前四阶频率的相对变化以及通过映射得到频率与实际在轨频率的相对误差,同理是仅考虑重力加速度和约束的情况,最后两列为同时考虑空气、重力加速度、约束三种力学环境因素影响的相对变化以及通过三级映射的得到频率值与在轨模型频率值相比的相对误差。由表可知,与复杂力学环境下的相对变化相比,相对误差数量级小3个以上数量级,因此频率映射是有效的。Table 1 shows the calculation results of frequency mapping. The two columns of air represent the relative change of the first four orders of frequency between the on-orbit and ground vibration tests and the relative error between the frequency obtained through mapping and the actual on-orbit frequency when only the air influence is considered. , the same is the case where only gravity acceleration and constraints are considered. The last two columns are the relative changes of the three mechanical environment factors including air, gravity acceleration, and constraints, and the frequency value obtained through the three-level mapping is compared with the frequency value of the on-orbit model. The relative error of the ratio. It can be seen from the table that compared with the relative change in the complex mechanical environment, the magnitude of the relative error is more than 3 orders of magnitude smaller, so the frequency mapping is effective.
表1频率映射计算结果Table 1 Calculation results of frequency mapping
对于振型映射能得到同样的结论。表2为前四阶振型最大位移相对变化和误差。与复杂力学环境下的相对变化相比,相对误差小2以上数量级,因此振型的映射也是有效的。The same conclusion can be drawn for mode shape mapping. Table 2 shows the relative change and error of the maximum displacement of the first four vibration modes. Compared with the relative change in the complex mechanical environment, the relative error is more than 2 orders of magnitude smaller, so the mapping of the mode shape is also effective.
表2前四阶振型最大位移映射计算结果Table 2 Calculation results of the maximum displacement mapping of the first four modes
根据得到的在轨航天器有限元模型的频率、振型,同步骤3,即可进行在轨航天器有限元模型的动力学响应计算。如图8所示为预测在轨时域响应,图9为预测在轨响应与模态坐标减缩模型响应相比相对误差随时间变化图,由图9可以看出响应预测相对误差在7%以内,因此预测是正确的。According to the obtained frequency and mode shape of the finite element model of the on-orbit spacecraft, as in step 3, the dynamic response calculation of the finite element model of the on-orbit spacecraft can be performed. Figure 8 shows the predicted on-orbit time-domain response, and Figure 9 shows the relative error versus time of the predicted on-orbit response compared with the modal coordinate reduction model response. From Figure 9, it can be seen that the relative error of the response prediction is within 7% , so the prediction is correct.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219792.6A CN104833466B (en) | 2015-04-30 | 2015-04-30 | A Spacecraft Ground Test and On-orbit Microvibration Dynamics Environment Mapping Method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219792.6A CN104833466B (en) | 2015-04-30 | 2015-04-30 | A Spacecraft Ground Test and On-orbit Microvibration Dynamics Environment Mapping Method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104833466A true CN104833466A (en) | 2015-08-12 |
CN104833466B CN104833466B (en) | 2017-05-10 |
Family
ID=53811490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510219792.6A Expired - Fee Related CN104833466B (en) | 2015-04-30 | 2015-04-30 | A Spacecraft Ground Test and On-orbit Microvibration Dynamics Environment Mapping Method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104833466B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106469256A (en) * | 2016-09-05 | 2017-03-01 | 重庆大学 | The determination method of the prestress annular film maximum stress of uniform load lower band hard core |
CN106599480A (en) * | 2016-12-16 | 2017-04-26 | 中国科学院长春光学精密机械与物理研究所 | Modification method of space camera on-orbit micro-vibration simulation model |
CN107941441A (en) * | 2017-11-14 | 2018-04-20 | 北京卫星环境工程研究所 | Determine the method that the in-orbit border of simulation influences the in-orbit dynamics of spacecraft |
CN108061660A (en) * | 2017-10-23 | 2018-05-22 | 上海卫星工程研究所 | The in-orbit thrust real-time calibration method of satellite engine based on linearly coupled measurement |
CN111444648A (en) * | 2020-01-15 | 2020-07-24 | 北京理工大学 | A Fast Estimation Method for Structural Dynamic Characteristics of Modular Spacecraft |
CN113340404A (en) * | 2021-06-21 | 2021-09-03 | 上海航天计算机技术研究所 | Low-vacuum modal measuring device for flexible antenna |
CN114813004A (en) * | 2022-06-29 | 2022-07-29 | 中国飞机强度研究所 | Loading attitude control method in airplane vertical fin vibration fatigue test |
CN117194959A (en) * | 2023-08-17 | 2023-12-08 | 西安交通大学 | Method and related device for identifying structural operation modal parameters of slender-body aircraft |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1006049A2 (en) * | 1998-12-03 | 2000-06-07 | Trw Inc. | Method and apparatus for integration and testing of satellites |
CN103488812A (en) * | 2013-08-12 | 2014-01-01 | 上海卫星工程研究所 | Satellite in-orbit micro-vibration simulation modeling correction method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH064252B2 (en) * | 1985-06-17 | 1994-01-19 | 株式会社明治ゴム化成 | Hose continuous vulcanization method |
-
2015
- 2015-04-30 CN CN201510219792.6A patent/CN104833466B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1006049A2 (en) * | 1998-12-03 | 2000-06-07 | Trw Inc. | Method and apparatus for integration and testing of satellites |
CN103488812A (en) * | 2013-08-12 | 2014-01-01 | 上海卫星工程研究所 | Satellite in-orbit micro-vibration simulation modeling correction method |
Non-Patent Citations (3)
Title |
---|
李华 等: "一个用于载荷地面测试的航天器仿真模型", 《电子设计工程》 * |
蒋国伟 等: "卫星微振动虚拟仿真技术研究及应用", 《计算机测量与控制》 * |
蒋国伟 等: "某卫星微振动建模与仿真", 《航天器环境工程》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106469256A (en) * | 2016-09-05 | 2017-03-01 | 重庆大学 | The determination method of the prestress annular film maximum stress of uniform load lower band hard core |
CN106469256B (en) * | 2016-09-05 | 2019-02-05 | 重庆大学 | Method for Determining Maximum Stress of Prestressed Ring Membrane with Hard Core under Uniform Load |
CN106599480A (en) * | 2016-12-16 | 2017-04-26 | 中国科学院长春光学精密机械与物理研究所 | Modification method of space camera on-orbit micro-vibration simulation model |
CN108061660A (en) * | 2017-10-23 | 2018-05-22 | 上海卫星工程研究所 | The in-orbit thrust real-time calibration method of satellite engine based on linearly coupled measurement |
CN107941441A (en) * | 2017-11-14 | 2018-04-20 | 北京卫星环境工程研究所 | Determine the method that the in-orbit border of simulation influences the in-orbit dynamics of spacecraft |
CN107941441B (en) * | 2017-11-14 | 2019-12-03 | 北京卫星环境工程研究所 | Determine the method that in-orbit boundary influences the in-orbit kinetic characteristics of spacecraft of simulating |
CN111444648A (en) * | 2020-01-15 | 2020-07-24 | 北京理工大学 | A Fast Estimation Method for Structural Dynamic Characteristics of Modular Spacecraft |
CN111444648B (en) * | 2020-01-15 | 2022-02-08 | 北京理工大学 | Method for quickly estimating structural dynamic characteristics of modular spacecraft |
CN113340404A (en) * | 2021-06-21 | 2021-09-03 | 上海航天计算机技术研究所 | Low-vacuum modal measuring device for flexible antenna |
CN114813004A (en) * | 2022-06-29 | 2022-07-29 | 中国飞机强度研究所 | Loading attitude control method in airplane vertical fin vibration fatigue test |
CN114813004B (en) * | 2022-06-29 | 2022-09-13 | 中国飞机强度研究所 | Loading attitude control method in airplane vertical fin vibration fatigue test |
CN117194959A (en) * | 2023-08-17 | 2023-12-08 | 西安交通大学 | Method and related device for identifying structural operation modal parameters of slender-body aircraft |
Also Published As
Publication number | Publication date |
---|---|
CN104833466B (en) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104833466B (en) | A Spacecraft Ground Test and On-orbit Microvibration Dynamics Environment Mapping Method | |
Yang et al. | Fuel-optimal control for soft landing on an irregular asteroid | |
CN103970964B (en) | Flexible satellite modal parameter in-orbit identification method | |
CN103471613A (en) | Parameter simulation method for inertial navigation system of aircraft | |
CN106844880A (en) | A kind of satellite momentum wheel disturbs shake experiment and data interpretation method | |
Gasbarri et al. | Control-oriented modelization of a satellite with large flexible appendages and use of worst-case analysis to verify robustness to model uncertainties of attitude control | |
Muralidharan | Stretching directions in cislunar space: Stationkeeping and an application to transfer trajectory design | |
Capristan et al. | Range Safety Assessment Tool (RSAT): An analysis environment for safety assessment of launch and reentry vehicles | |
Gao et al. | An optimal control approach to robust control of nonlinear spacecraft rendezvous system with θ-D technique | |
Norman et al. | Integrals of motion for planar multibody formations with internal forces | |
Schroeder | A Comprehensive Entry, Descent, Landing, and Locomotion (EDLL) Vehicle for Planetary Exploration | |
Nicholas | Attitude and formation control design and system simulation for a three-satellite CubeSat mission | |
Matsumoto et al. | FEM-based High-fidelity Solar Radiation Pressure Analysis | |
Kazemba et al. | Dynamic stability analysis of blunt body entry vehicles through the use of a time-lagged aftbody pitching moment | |
Johansson | Experimental and computational evaluation of capabilities of predicting aerodynamic performance for a mars helicopter rotor | |
Zanardi et al. | Spin stabilized satellite’s attitude analytical prediction | |
Cohan | Integrated modeling and design of lightweight, active mirrors for launch survival and on-orbit performance | |
Zuiker | Linear covariance analysis of atmospheric entry for sample return mission | |
Tagliani | Design and Validation of a Flexible Spacecraft Model for Attitude Control Applications. | |
Mosquera Alonso | State-and parameter estimation for spacecraft with flexible appendages using unscented kalman filters | |
Elshaal | DESIGN ANALYSIS OF THE STRUCTURAL SYSTEM OF ALAINSAT-1 CUBESAT BY SIMULATION AND EXPERIMENTAL VALIDATION | |
Ripoll | Assessment of Spacecraft Aerodynamics under Uncertainty toward the Optimization of Atmospheric Entry Trajectories | |
Abtahi et al. | Chaos and Bifurcation in Roto-Translatory Motion of Gyrostat Satellite | |
Albert | Aerocapture, Entry, and Co-Delivery in Uncertain Planetary Atmospheres | |
Solari | Reinforcement Learning Guidance for Pinpoint Landing on Mars |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170510 |
|
CF01 | Termination of patent right due to non-payment of annual fee |