飞行器结构动力学仿真验证

问题的来源

客户承担的某型复合材料无人机机翼结构验证任务，在试验中出现了局部屈曲：加载到设计载荷的某一阶段后，翼面靠近连接区的位置出现可见变形，与预期不符。要回答“为什么在这里失稳”，必须先能复现试验过程，而原先使用的国外软件采用的是壳单元加绑定约束的简化建模方式，连接区的接触状态与铺层的损伤起始都没有被显式表达，计算出的应力分布无法解释试验现象。

复核任务的约束很明确：不能重新做试验，只能靠计算回溯；结论要能对上试验测点数据，否则无法支撑设计更改；而型号进度不允许无限期地反复试算。

难点不在算得动，而在算得对

复合材料机翼的建模有三个绕不开的技术点。

第一是连接区的接触非线性。螺栓连接处存在装配间隙与预紧力，加载过程中接触面会分离、滑移，刚度随载荷变化。简化模型用绑定约束把两层“粘死”，等于假设刚度恒定，在接近失稳的工况下偏差极大。要准确复现，必须显式建模接触对，而这会带来收敛困难——接触状态在迭代中反复跳变，常规求解器很容易在某个载荷步停滞。

第二是损伤起始与演化。复合材料的失效不是一次性断裂，而是基体开裂、分层、纤维断裂逐步累积。这里采用 Hashin 准则判断损伤起始，并引入渐进损伤演化描述承载力随损伤扩展的退化过程。难点在于材料参数需要标定，直接采用手册值往往与实测偏差较大。

第三是模型规模与求解时间的矛盾。全尺度模型要覆盖机翼全部细节，单元数量大；而回溯分析需要批量扫描多个工况，单次求解太慢就失去了工程意义。

建模与求解的组织方式

做法是把问题拆成两层：先用连接区的局部细节模型与试验件数据对齐，标定接触参数与损伤参数；参数可信之后再放进全尺度模型，避免在整机上反复试错。

• 全尺度模型采用壳与实体混合建模，蒙皮、翼梁用壳单元，连接区与厚截面区域用实体单元，兼顾精度与规模
• 螺栓连接显式建立接触对，按装配工艺设定预紧力，允许接触面分离与滑移
• 铺层按实际铺层顺序定义，损伤模型覆盖基体拉伸／压缩、纤维拉伸／压缩四种模式
• 静力与模态两类工况做成参数化模板，批次提交、自动汇总，替代逐工况手工操作

实施中的两次回头

第一次回头出现在损伤参数标定阶段。初始采用手册值计算，失稳载荷明显高于试验值，说明模型“太结实”。回到局部模型逐项对比后确认，基体压缩的断裂能取值偏大，导致损伤演化过于缓慢。按试验件数据重新标定后，失稳载荷与试验吻合。

第二次回头是网格密度。连接区接触面附近的网格偏粗时，接触压力的分布出现明显锯齿，收敛也随之变差。加密之后收敛性反而改善——接触问题的非线性迭代对网格质量比线性问题敏感得多，这一点在方案初期被低估了。

结果与沉淀

最终单工况求解时间由 6.5 小时降至 1.4 小时，提速 4.6 倍，使批量工况扫描成为可行做法；模态频率与敲击试验对比误差 2.7%，损伤起始位置与试验观察到的位置一致，回溯结论直接支撑了连接区的设计更改。

更重要的是形成了可复用的东西：一套连接区建模准则（接触定义方式、预紧力设定、网格密度要求），一套材料损伤参数的标定流程，以及静力与模态的参数化模板。这套东西随后在该院三个型号系列上复用，新项目的建模时间不再从零开始。