《Aerospace》：Recent Advances in Multiaxial Shock/Vibration Environment Simulation and Damage Evaluation for Aircraft
编辑推荐：
现代航空航天飞行器在整个服役过程中处于复杂的多轴动态环境中，包括多向耦合冲击、随机振动、谐波振动以及气动弹性激励。传统的顺序单轴试验虽然便于工程实施，但无法复现运行环境的空间相干性、相位关系和非线性耦合效应，这可能导致对结构响应、疲劳损伤和机载设备功能退化的评
现代航空航天飞行器在整个服役过程中处于复杂的多轴动态环境中，包括多向耦合冲击、随机振动、谐波振动以及气动弹性激励。传统的顺序单轴试验虽然便于工程实施，但无法复现运行环境的空间相干性、相位关系和非线性耦合效应，这可能导致对结构响应、疲劳损伤和机载设备功能退化的评估产生偏差。本综述总结了航空航天飞行器多轴冲击/振动环境模拟与损伤评估的最新进展，涵盖了多轴载荷源分类、传递路径与耦合机制、动态响应分析、冲击与振动诱导损伤模型以及更新的实验测试技术。该工作进一步识别了载荷谱测量、谱矩阵复现、非线性路径识别和冲击-振动耦合损伤评估中尚未解决的挑战，并提出了未来研究方向，包括多轴服役环境数据库构建、统一响应与损伤等效准则，以及结构、设备和功能性能的综合评估框架。本综述为航空航天系统的环境设计、结构鉴定和可靠性评估提供了技术参考。
- 1.
引言
航空航天飞行器的可靠性与安全性取决于其承受复杂动态环境的能力，这些环境出现在发射、滑跑、分离、机动、着陆、舰载操作、运输和维护等过程中。多轴环境模拟对于实现真实的设计裕度、鉴定试验、基于状态的维护和任务可靠性评估至关重要。现有综述多关注通用系统的多轴振动响应与损伤，而本综述采用飞机专用的端到端处理方法，从运行源、传递路径、多轴响应到材料/结构损伤、设备功能和实验室复现进行统一分析。多轴动态载荷显著影响飞机结构强度、疲劳损伤演化、连接可靠性、航电设备功能以及乘员与任务安全，其响应并非单轴响应的简单线性叠加，而是由边界条件、载荷相位关系和模态耦合效应主导。
- 2.
多轴动态载荷特性
2.1 多轴冲击载荷
多轴冲击载荷由短时能量释放、接触约束突变、机械止挡冲击或高幅值外部激励引起，具有持续时间短、峰值高、带宽宽和方向耦合显著的特点。根据激励机理可分为四类：接触冲击（如着陆、拦阻着舰、地面碰撞）、火工冲击（如爆炸螺栓、切割器）、机构诱导冲击（如起落架收放锁定、舱门启闭）和分离冲击（如外挂投放、副油箱抛弃）。接触冲击在极短时间内改变局部接触状态，易激发结构接头、支撑结构和设备安装接口处的高频瞬态响应；火工冲击以短时、高频、高加速度的结构冲击波形式沿连接结构和外壳传播，对继电器、焊点、电子元器件等敏感部件构成威胁；机构冲击的幅值取决于作动速度、机械间隙、结构阻尼和锁定刚度；分离冲击则涉及气动干扰、执行器输入和约束突变的耦合。
2.2 多轴振动载荷
多轴振动载荷是飞机服役中持续时间长、频率成分丰富的动态载荷，主要表现为连续随机激励、周期谐波激励或叠加在宽带谱上的窄带峰值，来源包括机械激励（发动机、转子、螺旋桨）、气动激励（阵风、湍流、抖振、嗡鸣）、基础激励（跑道不平度、舰船运动、运输振动）和声学振动环境（发射声学、喷气噪声、空腔噪声）。机械激励产生与转速阶次相关的窄带峰值；气动激励具有宽带随机性和与结构变形反馈的耦合；基础激励通过支撑接口传入，其振动量级随滑跑速度变化而非平稳；声学振动环境对飞机蒙皮、隔板和设备外壳施加空间相关的分布压力波动，转化为局部振动响应。
2.3 冲击-振动耦合环境
飞行任务中的动态环境很少是纯粹的孤立冲击或振动，而是跨频段、跨方向的多个载荷源叠加，形成冲击-振动耦合环境。起飞阶段包含发动机加速、地面滑跑、抬轮、离地和起落架收起等过程，产生垂直卸载冲击、纵向加速度、俯仰姿态变化和衰减的低阶机体振动；舰载机弹射起飞则产生由弹射牵引、瞬时保持释放和前起落架回弹驱动的强耦合冲击。巡航和机动状态以宽带振动响应为主，但高g机动、阵风遭遇和气动弹性事件也会引发瞬态载荷峰值。着陆阶段，轮胎压缩和减震支柱压缩首先产生垂直冲击，静止的机轮在极短时间内被跑道摩擦加速到滚动速度，产生纵向自旋阻力，随后是轮胎和支柱结构恢复引起的弹性回弹响应；侧风条件下还会耦合垂直冲击与侧向摩擦、纵向阻力、滚转力矩和偏航力矩。
- 3.
多轴动态载荷传递路径与耦合机制
3.1 源-路径-接收者框架
源-路径-接收者（Source–Path–Receiver）框架将复杂动态环境分解为载荷源、传递路径和响应接收者三部分。传递路径分析（Transfer Path Analysis, TPA）通过界面力和频率响应函数（Frequency Response Function, FRF）建立源-路径-响应关联，用于分离物理路径的贡献值；工作TPA（Operational TPA, OTPA）利用工作响应数据诊断主导传递路径。多轴场景中，TPA不仅需检测最大单方向响应，还需考虑多点输入、多方向力和力矩、路径间相位差以及输入相干性。
3.2 从界面到结构的载荷传递
多轴动态载荷首先作用于支撑接口，经连接结构传递至主结构，最终到达设备安装接口或敏感部件。飞机结构中广泛采用螺栓、铆接和复合材料接头，这些连接区域存在预紧力、装配间隙和接触非线性。接头界面在冲击载荷下会改变载荷幅值、频率内容和传播方向，将局部冲击信号转化为由低频结构模态和高频局部波组成的复合响应。梁和壳结构作为飞机主承载介质，低频激励易激发整体弯曲、扭转、俯仰和滚转模态，高频冲击主要产生局部板壳振动和接头附近的波传播。
3.3 多轴载荷传递中的耦合效应
多轴载荷传递涉及方向、模态、频率和非线性界面行为的耦合。方向耦合表现为单方向激励引发其他方向的结构响应；模态耦合指弯曲、扭转、局部和刚体运动模态间的交互，受连接界面配置和边界条件控制；频率耦合表现为低频大位移响应与高频小振幅冲击或振动响应的叠加；非线性耦合是最具挑战性的方面，起落架减震器呈现强非线性刚度和阻尼，螺栓接头在冲击或随机振动下经历微滑移和预紧力松弛，改变结构传递特性。多轴载荷传递模型需集成非线性连接、接触和阻尼模型，而非过度简化为线性时不变系统。
3.4 多轴载荷环境重构方法
多轴环境重构旨在根据测量或模拟的结构响应反推载荷输入、主导传递路径和等效载荷谱，需要同步的多源、多点、多方向测量。载荷识别可通过逆动力学、FRF算法、TPA和OTPA技术实现。多输入多输出（Multiple-Input Multiple-Output, MIMO）建模在服役环境与实验复现间建立关键联系，传统单轴振动试验将复杂环境分解为三个正交方向分别施加，无法保留多方向激励间的相干性、相位差和力矩输入，而MIMO控制算法和谱矩阵校正能更好地逼近飞机服役环境的多轴边界条件。
- 4.
多轴动态响应分析方法
多自由度（Multi Degree of Freedom, MDOF）系统的模态模型通过有限元模型（Finite Element Model, FEM）的模态分析获得，其控制方程包含质量矩阵、阻尼矩阵和刚度矩阵。传递函数矩阵将运动随机激励与加速度响应以及应力响应关联起来，应力在频域由六个独立应力张量分量的自谱和互谱密度以6×6矩阵形式描述。
4.1 飞机多轴冲击动力学建模
多轴冲击动力学可通过三类模型描述：工程等效模型（如半正弦、梯形、锯齿波）用于试验规范和快速响应估计；有限元和多体动力学模型用于起落架落震、结构碰撞和局部应力分析；频域和响应谱模型用于火工冲击评估和高频敏感设备评定。起落架是多轴冲击输入和载荷传递的关键路径，需考虑纵向阻力、横向偏心载荷、轮胎-跑道接触、机体姿态变化和局部柔性连接。对于火工和分离冲击，建模需集成局部冲击源、连接接口阻抗、结构波传播、模态响应和设备敏感频段。
4.2 飞机多轴振动动力学建模
多轴振动动力学建模需要统一描述刚体运动、弹性振动和非定常气动特性。对于高展弦比柔性机翼，传统线性小扰动模型往往无法准确预测响应，需采用非线性模态气动伺服弹性框架，结合几何精确梁理论和ONERA非线性非定常气动模型。几何非线性是柔性飞机多轴振动建模的核心考虑，共旋梁理论可用于大变形下的动态响应和稳定性分析。对于机载设备和有效载荷，多轴振动建模必须充分考虑实验室试验边界条件与实际安装边界的差异，包括夹具柔性和多点输入关系。
4.3 多轴冲击载荷响应分析方法
冲击脉冲表征不仅需包含峰值加速度，还需考虑持续时间、速度变化、主导频率成分、能量分布和多方向同步性。冲击响应谱（Shock Response Spectrum, SRS）通过将输入冲击信号施加于一系列具有不同固有频率和固定阻尼比的单自由度系统，记录每个系统的最大响应来评估设备对频率相关冲击的敏感性。多方向冲击同步是多轴与单轴冲击响应分析的核心区别，顺序单轴冲击试验破坏了方向峰值响应的时间相关性，无法复现模态耦合和损伤累积路径。多点冲击试验中需采用随机延迟叠加和滤波叠加等方法构造满足目标SRS的响应信号。
4.4 多轴振动载荷响应分析方法
多轴振动响应分析分为时域分析、频域分析、模态分析和响应谱或疲劳谱分析。时域分析适用于瞬态冲击和非线性接触，但计算成本高；频域分析将输入功率谱密度（Power Spectral Density, PSD）矩阵转换为响应PSD矩阵，用于计算均方根响应、峰值估计、疲劳损伤和界面载荷，需考虑谱矩阵的正定性、相干性、相位和互谱特性。多轴疲劳损伤谱（Multiaxial Fatigue Damage Spectrum, MI-FDS）将传统单输入/单输出FDS理论扩展到MIMO随机环境，强调了多轴疲劳损伤评估中相位和相干性的重要性。
4.5 不确定性表征与传播
飞机多轴响应预测受服役载荷和材料属性的认知不确定性、模型形式和边界条件的不确定性以及测量/控制不确定性的影响。不确定性分析方法包括蒙特卡洛、多项式混沌和随机配置技术、区间和非概率方法以及贝叶斯更新。主要障碍是计算成本，可能需要数百或数千次重复动态分析，降阶模型和代理模型可减轻负担。不确定性也进入反演识别和试验设计，复合材料层合板损伤识别和作动器/传感器布置需报告置信区间或区间边界。
- 5.
多轴冲击/振动损伤评估
5.1 多轴冲击损伤模型
对于铝和钛结构，Johnson–Cook模型因其紧凑形式而广泛用于高速度冲击、叶片包容和鸟撞事件筛选，但经典公式对Lode角和剪切主导断裂的敏感性较弱，需要修正的Johnson–Cook、三轴度/Lode相关断裂或耦合损伤模型。Cockcroft–Latham断裂准则适用于拉伸主导载荷下的延性金属冲击失效，而嵌入应力三轴度和Lode参数的延性断裂模型更适应多轴冲击下的复杂断裂。对于复合材料结构，多轴冲击损伤模型关注各向异性损伤演化、层合板结构特征和层间界面性能。三维Hashin准则识别四种失效模式，适用于CFRP/GFRP层合板低速冲击后的面内损伤预测；Puck准则可靠描述基体开裂和纤维间失效；LaRC系列准则在纤维压缩不稳定性、基体开裂和剪切非线性评估中表现优越；内聚力模型（Cohesive Zone Model, CZM）主要用于层间分层扩展；虚拟裂纹闭合法（Virtual Crack Closure Technique, VCCT）适用于预存裂纹和分层缺陷的扩展分析。
5.2 多轴振动疲劳损伤模型
多轴振动疲劳需将具有互谱相位的六分量应力过程转换为循环和材料损伤度量。频域等效应力模型基于von Mises二次屈服准则将多轴应力PSD数据转换为等效单轴随机过程，便于与FEM随机响应分析结合。临界平面模型通过从应力张量PSD矩阵确定临界平面方向，建立等效法向应力PSD进行疲劳寿命估计，对拉伸-剪切耦合和裂纹方向预测适应性更好，但计算资源需求更高。多轴FDS将输入相关性纳入损伤潜力评估框架，适用于飞机设备、支架、舱体结构和机载电子设备的试验环境剪裁。将多轴随机激励损伤与顺序单轴激励损伤比较可建立损伤修正因子，但依赖于结构线性、平稳载荷和高斯响应等假设。
5.3 复合材料冲击后剩余强度
复合材料飞机结构的工程目标是预测冲击后的剩余承载能力，包括冲击后压缩强度、拉伸剩余强度、刚度损失及向疲劳驱动损伤的过渡。Hashin准则计算效率高，但标准公式未明确解析基体主导失效的断裂面、压缩下的纤维扭结或层间分层，其剩余强度预测受所选退化律和单元特征长度强烈控制。Puck和LaRC准则提供了对剩余强度更详细的描述，Puck的纤维间失效公式识别物理驱动的断裂面，LaRC公式对纤维压缩、扭结带形成和剪切非线性具有额外优势，但两者均需独立定义的演化律、断裂能和刚度退化方案。CZM通过牵引-分离定律显式表征分层扩展，对低速冲击和冲击后压缩分析有效；VCCT从已有裂纹前缘的能量释放率评估分层扩展，适用于实验测量或前序模拟映射的冲击损伤。
5.4 复合材料冲击疲劳退化
冲击后疲劳分析需复现初始冲击损伤及其在数千或数百万次载荷循环中的演化。冲击产生的基体裂纹、纤维断裂、分层和局部压痕改变刚度场并重新分配循环应力。Hashin模型可将冲击产生的面内损伤变量转移到后续疲劳分析中，但需附加疲劳演化律和应力比修正。Puck和LaRC公式对基体主导和压缩驱动的疲劳机制提供更高物理分辨率，疲劳内聚公式将单调牵引-分离律扩展为循环相关的界面损伤，需引入耐久阈值和模式相关扩展率。VCCT为已有冲击诱导分层的疲劳扩展提供断裂力学框架，将I型、II型和III型能量释放率与Paris型疲劳扩展律关联。
