关于航天动力学环境的最新进展与技术展望
动力学环境问题是飞行器设计中的基础与关键，国内外航天器在发射与飞行中振动、冲击造成的故障占据了相当大的比例，轻者影响精度和飞行稳定性，重者会造成失败。随着航天技术的发展，我国航天动力学取得了前所未有的成绩。同时由于信息、材料、微机械等技术的进步，飞行器也呈现出新的特点，动力学问题也出现了新的发展方向。
1 航天动力学环境预示与虚拟试验技术
1.1 动力学环境预示技术进展
动力学环境预示是航天飞行器设计的基础。近年来，航天动力学在环境预示方法上有所突破，开始摆脱对国外标准的依赖，并逐步走上工程实测数据持续修正的道路。一方面覆盖高频、中频、低频的全频域混合建模预示技术得到快速应用，为研制初期的动力学条件提供输入;同时采用型号研制累积的大量实测数据不断修正，某些产品形成了比较精确的动力学模型，制定了有别于国外与国军标的振动、冲击条件，并据此探索了多种经验外推、混合建模、非平稳数据处理等技术。
1.1.1 有限元-统计能量混合建模技术
飞行器在发射、飞行中受到振动、冲击和噪声影响，早期飞行器问题主要集中在低频部分，一般采用有限元方法就能解决。随着飞行器的发展，轻质、柔性、次级结构、设备内部结构等作用凸现，中高频问题凸现，所以近年来对FE-SEA 混合方法的研究很多。主要针对航天器等复杂系统的特点，以飞行器典型结构为主，利用板梁组合结构开展数值仿真，并细分子系统连接关系，使复杂连接随机子系统间能量传输关系具有了完整表达式。同时提出了损耗因子修正、推力修正、设备传递特性修正三级修正方法。完成了诸如卫星天线、太阳帆板、整星结构、运载器、导弹仪器舱、超声速飞行器等基础激励和混响声场动响应，并与试验数据进行了对比。中频段预示与实测结果之间吻合较好，甚至部分精度已经达到3 dB，但存在某些整体模态被遗漏、响应结果应用到复杂系统局部细分不够等问题，模型还需进一步完善。
1.1.2 振动冲击实测数据统计条件与外推技术
由于工程结构的复杂性，确定动力学环境的方法除混合建模外，还大量采用实测数据或相似产品外推等统计归纳与理论结合的方式，并取得了良好的效果。对气动激励载荷一般采用压力脉动计算或测量值施加到飞行器表面，也采用计算气动声学CAA 的方法。该方法弥补了传统CFD 宽频声波不同波长传播不能准确刻画的缺点，与CFD 紧密关联成为今后的重要方向，但☣在比例模型试验方面，需要对气动声学的相似律加以详细研究。振动冲击环境外推法是在参考型号有直接测量结果、新型号与参考型号的动力载荷自谱相似、测量位置结构近似的前提下，按照激励源缩比、结构表面质量密度缩比、飞行器直径缩比等三个基本缩比原则对新型号振动响应自谱量级进行外推的一种方法，弥补了有限元和统计能量缺点，简单实用，适用新型号研制，应用相当广泛。我国积累了大量数据，逐步形成了诸如根据设备固有频率特性确定半正弦冲击和能量谱点火冲击条件的制定方法，根据飞行数据确定振动环境条件并据此给出了马赫数从低到高的外推公式等经验公式和标准。由于数据积累系统性不够，基础数据挖掘不到位，与产品结构、性能等参数关联度不高，应用水平与国外相比还有较大差距。
1.1.3 非高斯、非平稳振动数据处理技术
飞行器在运输中因路面质量会诱发结构碰撞，飞行中因推力脉动、姿态变化等会产生非平稳振动，实测数据分析表明，其具有明显的超高斯分布特征。试验与仿真表明，超高斯分布振动会大大降低结构的疲劳寿命，采用宽带信号作为激励进行疲劳仿真分析，随峭度的增大疲劳损伤积累加快，疲劳寿命缩短。目前振动数据处理及振动试验条件均按照高斯分布进行，忽略了功率谱密度相同但时域数据分布不同对结构影响的差异，在某些特定的情况下应予以重视。
1.2 航天动力学环境预示发展方向
航天动力学环境主要方向是解决高速飞行器气动热带来的热结构动力学问题;动力学环境如何准确定位到设备甚至功能单元一级;如何建立起地面试验结果与飞行状态下的响应之间的关系等问题。
1.2.1 高速飞行器热结构动力学环境分析预示技术
高速飞行器是世界各国发展的重点，特别是长期在大气层飞行情况下与以往发动机置后不同，采用的是吸气式冲压发动机，具有飞行器结构推进一体化特点以及影响控制系统设计的气动/推进/弹性耦合效应。因此，外部气流与进气道流场同时作用于结构，有严酷的振动、噪声环境;另一方面，因气动加热、发动机高温影响，结构特性随温度变化呈现时变特征，动力学环境具有全频域以及与高温耦合的特点。因此，气动激励源作用机理、吸气式飞行器内外流激励振动环境预示与控制、流-固-温耦合物理场分析、热振声多物理场分析、热气动饲服弹性分析、复杂系统时变动力学特性分析、温度对复杂结构振动固有频率影响机理、热虚拟试验等技术成为目前研究的热点。目前在冲压发动机激励源预示技术方面取得了一定突破，在多物理场、特性分析、热饲服弹性分析等方面采用数值模拟、有限元分析与风洞试验、热模态试验等相结合的方法取得了较大进步，但对内外流场致振机理等尚不清楚，热结构试验与计算结果之间还存在很多不解之处。
1.2.2 基于设备功能的动力学环境预示技术
动力学环境预示的目标是使飞行器设计达到轻量和良好的特性。目前的预示方法采用输入结构传递响应分析方法，属于正向分析，只有结构确定后才能进行。尽管对结构优化、完善设计提供了帮助，但对结构系统设计难以发挥根本性改变。其次，有限元、统计能量方法也存在不够细化、不能预示到设备和功能单元局部的缺点，只能给出全系统级或者舱段级的环境，即动力学环境预示与目标设备/单元功能、结构设计合理性以及环境适应性评价之间的桥梁没有建立起来。因此，近年来提出了基于设备功能的反向分析方法，即逆问题。基本思路就是根据飞行器设备的功能指标敏感性要求，把功能指标分解到设备安装部位和敏感物理量上，甚至对结构某些重要位置提出控制指标要求，这样会缩短设计周期，并一次性形成设计结果优化。
目前对该方法的研究取得了初步成果。其基本过程是：首先对飞行器不同位置、相同位置不同功能设备的动力学环境敏感量进行分析，在位移、速度、加速度等18 个物理量中筛选，并分解到安装结构位置甚至到结构传递路径，作为显性指标用于约束结构设计和评估，结构完成后再进行预示，对结构目标点预示结果指标进行转化，与分解指标进行对比评估。这样，利用不同设备功能敏感的动力学环境特征量的不同，将预示结果与结构设计约束之间建立起直接联系，是未来飞行器结构设计的发展方向。该方法已初步用于惯性器件、光学导航、舵系统的动力学环境预示分析中，初步确定了功能设备敏感指标体系，并给出了指标转化方法，形成了一套较完整的思路。
1.2.3 飞行器动力学环境天地一致性问题研究
分析并模拟产品结构所在位置动力学环境的有效性是结构响应预测的关键之一。飞行器结构复杂，且具有多源输入特征，飞行中边界、气动附加刚度、载荷、噪声等影响因素很多，产品分析与地面试验的状态、模型等一致性难以得到保证，传统方法对飞行复杂环境的准确预示比较困难。同时，因振动冲击遥测通道数量有限，难以统计归纳以描述飞行器各部位状况。热试车、发动机试验、火箭撬、功能试验等仅验证了飞行状态中的一个特定状态，但地面试验数量很大。天地一致性正是综合利用地面试验、分析和飞行遥测数据，根据地面试验系统与飞行状态下结构响应之间存在的映射关系，采用基于映射关系的动力学环境预示技术，对未知部位进行预测。目前，在多个飞行器设计中得到了初步应用。主要采用地面模态试验、地面试验测量数据建立模型，消除地面边界条件影响，运用支持向量回归机构建立两者之间的映射关系模型。通过大量地面试验与遥测数据作为训练样本，得到了回归模型曲线，并在典型结构上得到了试验验证，实现了由遥测点到其他结构点的预测。应用还存在不少缺点，一是遥测点数量有限，映射模型精度较低;二是遥测飞行环境较好，遥测数据样本量不够完备，边界条件模型精度不高，需要结合仿真分析完成。
1.2.4 动力学环境虚拟试验技术
虚拟试验是随航空航天等技术发展和任务增加及大量数据积累而产生的，其意义在于：一方面产品模型包含了不同激励下的响应信息，可分析方向性结果，降低人工判断遗漏，减少实物试验数量，降低成本;另一方面，实物试验前对试验设计方案进行快速评估，确定其合理性和可行性以及安全、质量、经济性。同时对于超出试验能力范围或其他限制的，或难度大、成本高以及无法准确模拟极限载荷工况情况的，更适合采用虚拟试验技术对其进行考核验证。虚拟试验建模过程侧重于✉大量试验数据积累和不断修正的结果上，因此需要一个循序渐进不断完善的过程。国外起步较早，美国已形成基于虚拟试验验证技术的公共支撑框架TENA，以及以JMETC 为代表的大型综合试验验证系统，并在F-22、宙斯盾等武器系统中得到应用。国内取得了一些阶段性的应用成果，军工产品虚拟试验领域已形成自主知识产权的通用支撑技术框架VITA。主要困难在于数据积累及其规范性较差，比较简单的结构和单一因素的虚拟试验方法已取得良好进展，而复杂结构多物理场耦合环境方面还存在很大距离。
2 航天动力学环境试验技术
2.1 大型复杂结构模态试验技术
模态试验技术发展较成熟，在单点与多点正弦激振的基础上发展了多点随机激振和综合模态试验技术，目前已广泛用于航天等领域结构动态设计、振动控制和故障诊断中。我国航天器、运载器和发射装置等大型系统均开展结构模态分析、模态试验，研制过程中要进行结构部组件、全系统、飞行器/发射箱/发射架/平台联合模态试验，以确定结构系统特性对发射精度、飞行稳定性和设备环境适应性的影响。大型复杂结构模态试验方法主要有宽带随机激振法、纯模态法、锤击法、步进正弦法、工作模态法等。模态参数辨识方法常用的有频域直接拟合法、时域多参考点复指数法LSCE 算法、工作模态polymax 算法等。大型模态试验技术发展有高温环境下的热模态试验技术、三维/多体等非梁模型识别与应用、时变结构模态参数识别方法。
在航天产品中往往存在多体、集中质量、间隙、弱连接、柔性结构等结构非线性以及刚度、质量、阻尼等参数非线性，相应技术尚不成熟，目前多用于飞行器操纵面模态试验。方法主要有测量静刚度法、测量动刚度法、消除间隙和摩擦的橡皮绳加载法以及改变激振力的纯模态试验方法。测量静刚度法是测量结构在不同静力载荷下的刚度。测量动刚度法则是给结构施加可变的正弦激励载荷，在考虑刚度引起的非线性同时引入了阻尼引起的非线性。橡皮绳加载法则是用弹性力消除结构间隙和摩擦，然后测量结构线性化后的模态参数。改变激振力的纯模态试验方法跟线性结构纯模态法类似，通过改变激振力的频率和大小、得到不同激振力下的结构模态参数[1516]。
2.2 微机械模态试验技术
新型飞行器大量采用惯导对准、地形匹配、光学制导等技术，相应的印制电路板、辅助结构、机构、支架、内部结构的动力学特性直接影响飞行稳定性和设备精度。由于此类结构质量轻、刚度小、结构复杂，传统模态试验方法带来的附加刚度和质量与被测结构相当，精度难以保证，因此采用声激励和光测替代传统激振杆或锤击法等接触式激励和加速度传感器测量方法，可以弥补相应不足。重音喇叭发出的声能足以激发出小质量结构件的模态，ล但要解决声源与现有设备的阻抗匹配、相位匹配等问题。目前国内有采用低频高密度体积声源集成的非接触激励，加速度测量采用激光测振仪，消除了附加刚度和质量的影响。通过对比试验，两种激励方式得出的频率、振型结果一致，声激励试验结果与典型算例结果更吻合，精度更高。
2.3 多输入多输出振动试验技术
目前比较流行的振动试验做法是用多个振动台模拟飞行器运输、发射与飞行动力学环境，以模拟多激振源，并提高激振力分布均匀性或提高推力以及柔性结构的试验效果，并在军事、航空、电子等领域广泛应用，美军标和国内标准均有相关内容。同时国外SD，DP，IMV，LMS 等公司拥有多台振动控制仪商用产品。国内在20 世纪90 年代就开始了这方面的研究，在细长体双振动台试验、设备级三轴振动试验方面应用比较成熟，并形成了相关行业和国军标等标准。南航、亿恒等单位多台控制技术也逐步成熟并形成了商用产品。目前中国航空综合技术研究所等正在建设大型多振动台试验系统，但大型产品结构多样，特性各异，而多输入多输出控制方法依赖于产品结构传递特性，所以通过实测数据统计归纳而形成试验条件，包括控制矩阵对角元素设定在内的多振动台控制还需要深入研究。此外还有利用多个压电元件来模拟实际环境的分布式振动试验技术。制定合理的多台多轴的振动试验条件，针对具体试验件研究激励方案以及多轴振动与热环境耦合试验方法等问题是以后研究的重点方向。
2.4 离心复合试验技术
飞行器在机动转弯等动作会产生过载，同时有温度、振动、冲击等共同作用的复合环境。对于特定产品结构，离心复合环境会影响到产品的精度和功能。离心复合试验在国内外发展比较迅速，美国的UCDavis 大型振动离心机系统包含了水平单向、垂直/水平双向振动、试验摄像观测系统、电阻层析成像系统等，是目前振动离心机的最高水平。目前国内的离心复合试验技术通常是以离心机为主体，在机臂上安装温度箱、振动台来实现的。这种方法存在的技术难题是离心力作用下振动台动圈会发生偏转，须设计纠偏系统，同时必须考虑对振动台固有特性的影响。目前国内多为小型振动离心机，并在航天、水利水电以及核工业得到应用，但振动台推力偏小限制了试验系统的广泛应用。中物院等单位对离心复合试验系统进行了深入研究，并根据需求提出了模拟飞行过程的高动态离心试验、STS 环境模拟等试验模拟技术。
3 热结构动力学试验技术
由于高超声速飞行器速度高，动力学特性与以往有很大不同，尤其以高温耦合影响为重。在热结构动力学建模分析的基础上，根据考核目的应尽量减少耦合因素，突出主要矛盾，可以分别进行两两或多耦合因素试验验证。
3.1 热结构时变模态试验技术
由于气动加热效应等影响，高速飞行器结构受到温度梯度影响，相应的材料物性参数、结构刚度、热应力等随飞行时间变化，飞行器动力学特性随时间变化并影响到飞行控制。高温影响机理复杂，对简单结构而言，以材料弹性模量下降为主，符合结构固有频率下降规律;对复杂结构，热应力导致刚度分布变化，结构特性难以预测，各阶固有频率并不完全符合随温度升高而降低的直观规律。
结构热模态试验是获取热结构动力学特性的主要方法，但由于地面试验受到的影响因素更多，热边界影响以及结构边界受热影响的刚度变化复杂，也是国内外研究的前沿与热点，成为集控制、信号处理和结构动力学的一个交叉学科，借助线性系统和信号处理等领域的研究成果发展出一些有效的时变模态参数辨识算法。国外如NASA 早期研究成果表明，温度对结构模态参数的影响明显，但振型基本不变，均为基于傅里叶变换的模态参数辨识算法，仅是温度稳定环境结论。国内基于高超声速飞行器发展，提出了基于递推子空间方法的时变模态参数辨识算法，可以通过状态空间联系到结构的振动微分方程，找到彼此的契合点，这种方法也是参数化的模型，识别精度更高。工程中一般分两步完成，先获取不同温度稳定环境下的模态参数，以此为基础再获取结构模态参数随着温度变化的详细变化。
3.2 热振动/噪声试验
与热耦合的动力学试验用于验证不同温度、湿度分布下结构的传热特性，温度分布以及振动、噪声、疲劳特性。由于影响因素多、技术难度大，到20 世纪90年代才开始逐步完成工程实现。飞行器昂贵且材料存在不可逆因素，一般按热振动、热噪声、热疲劳等两两组合进行试验，从飞行器经历的环境剖面及飞行器局部分析，提炼需解决的主要矛盾，分析、设计多组试验，进行相对简化的环境因素组合模拟和覆盖性考核与评估。目前具有热噪声试验设备和试验能力的国家有德国、俄罗斯、美国等，根据考核的主要目的进行部件、舱段、整机试验。国内航天、航空等试验机构均建成了较大规模的热试验设施，具有系统级1300 ℃以下热振动、热强度、热疲劳、热模态试验能力，部件级650 ℃以下热噪声试验能力。热噪声试验主要针对薄蒙皮、TPS结构等对热、噪声联合作用敏感的结构。
3.3 热结构动力学参数非接触式测量技术
由于高温的影响，传统的应变、振动、压力测量手段受到限制，因此，近年来出现了适宜高温环境下测量的新技术。
非接触式振动测量技术主要是应用多普勒原理、激光干涉原理的测振方法，因其非接触、无附加质量、动态范围大和精度高，所以在热模态、热振动等试验中广泛应用。由于其存在安装对准、激光反射效果、热蒸腾效应问题，影响了振动测量精度。
非接触式热应变测量方法有基于数字图像处理的方法和数字图像相关方法等。此外还有电子散斑干涉法、云纹干涉法等方法。这些方法测量灵敏度非常高，并具有测量结果直观可视的优点。
高动态应变测量技术主要采用光纤Bragg 光栅动态应变传感器，其特点是信号频率高、信噪比高。瞬态应变测试在实用的信号解调方面还存在响应速度方面的瓶颈问题，尚需做大量工作。
压力脉动测量技术包括有以下几种：传感器外加冷却护套，但存在局部冷区大热梯度和声歪曲;光纤传声器，在727 ℃以下有良好频响特性;等离子传声器，可用于高温、高速、高焓值气流中的噪声、热流等多参数的同步测试;将高温下的噪声用波导管导出测量，已在发动机内部噪声测试等任务中得到应用。
4 气候环境作用下复杂结构动力学问题
飞行器地面阶段✯运输、执勤环境复杂，存在平台、工ฏ作/人工诱发与气候组合环境影响，其振动冲击响应不但有其特殊性，而且有时比飞行环境更为严酷。
飞行器-发射装置-平台等构成的系统动力学问题主要体现在以下方面。
1)载荷激励存在平台多样性，发射状态下风载荷激励，车辆运输中道路障碍等情况下的颠簸冲击、刹车冲击、振动，舰船平台受海浪冲击、炮击、爆炸冲击传递到产品的振动冲击响应等[24]。
2)连接结构存在诸如多联装刚弹转换、适配器、运输中的减震器阵列等特殊情况，以及车辆发射装置起竖折叠状态下的多体结构，存在结构非线性，在温度变化等情况下动力学特性会发生相应的改变。
3)发射装置、地面设备往往尺度较大，存在大尺度效应，在温度、阳光、风载荷、平台诱发等叠加或综合作用下，产生局部变形等，即存在局部微环境问题，长期作用下产生疲劳或应力腐蚀。
4.1 飞行器-发射装置动力学响应分析
地面阶段主要表现为运输中的冲击、振动，舰载爆炸冲击、风载荷等情况。近年来国内外利用NASTRAN 等软件开展了大量研究，特别是对水下爆炸冲击响应进行了深入分析，采用大质量法、动态设计分析等方法，综合考虑了分布质量惯性力、结构传播效应等因素，给出了发射装置的变形和应力，并得出大长细比产品对冲击敏感、横向响应较大、沿长度方向响应差异较大等结论。由于系统构成复杂，很难得到精确的振动应力，只能采取实测数据。通过大量实测数据统计分析，得到了直接用于工程的成果，火箭在运输车上的振动应力为正态分布或威布尔分布，对组合系统力学特性影响很大。通过比较公路、铁路、海上、飞机运输振动特点，公路的振动量级最大，且主要以低频振动为主。通过实测数据也发现，在大外载荷作用下，如在路况很差、风浪很大的情况下，由于连接结构、接触面、安装固定面的存在，可能会激发出多平台固有频率以外的高频冲击振动，并具有超高斯分布特征。
4.2 飞行器-发射装置-平台系统动力学特性研究
地面阶段存在连接状态差别大、边界不确定性大、导致的结构非线性变化等影响，诸如运输状态下的弹性边界、动态载荷下边界力变化、待机发射状态边界、车/舰发射过程中的适配器滑动状态变化、机载发射、不同温度/湿度/风载荷影响、水下发射影响等。仅部分情况一阶固有频率计算较为准确[28]，且不同测试方法模态频率差别较大。目前产品研制过程中均要完成水下模态、发射箱/架/飞行器、车辆/发射箱/飞行器、飞机/发射架/飞行器等稳定状态的特性试验，除研究飞行器在平台的性能变化外，还可以获取飞行器对平台系统的影响，如是否会与机翼、发射箱等发生共振，以及如何改进结构，降低系统各部位的动态应力。
4.3 气候-工况/平台诱发载荷组合环境下动力学特性研究
飞行器在地面阶段尽管温度比较温和，但平台诱发、气候环境的复合及累积效应突出。舰船平台上，除经历摇摆、冲击和振动外，还有阳光照射的光化学和温升效应。大型产品大尺寸效应产生的累积变形和应力不容忽视，与盐雾、潮湿等形成气候-力学复合环境，加剧了结构的破坏。地面车辆平台经历颠簸、刹车等过载和冲击以及气候环境影响。平台上大型复杂结构系统内部可能会存在着大量高于设计或整体平均设计条件的局部微环境，如大型结构中的减震器阵列在外力诱发下的不均匀分布会导致某个减震器失效并导致连锁反应整体损坏，固发中的装药界面应力等，影响飞行器寿命。系统的损坏最终可以归结为某个部位、构件、器件的故障，因此，目前的解决方法主要是把地面阶段飞行器与平台作为一个多体、柔性、非线性系统分析，确定关键薄弱环节，并根据局部结构的工况载荷应力和微气候环境等效成典型试验构件，附加少量整机级试验，评估其适应性和寿命。此外，为防止偶发事件和系统复杂性难以判断的破坏问题，普遍采用自动记录设备，可以长期监测关键部位环境数据，并与设计基线比较，预判出相应的问题。
基于上述分析，目前在该方面的发展方向主要有：气候作用下的大尺度装备动力学分析及疲劳寿命研究;气候环境下复杂结构力学特性分析技术;飞行器-发射装置-平台系统微环境分析;健康管理与诊断技术等。
5 结语
飞行器动力学环境问题会随着工程需求不断出现新的变化，但解决方法不外乎分析与试验两种方法。通过工程实践，分析预示、试验与产品功能密切相关，型号结构、功能、材料、器件的变化，会带来分析与试验方法的变化。
