摘要
针对复合翼无人机多飞行模态、多载荷工况带来的结构强度、刚度与轻量化矛盾，提出基于达索系统3DEXPERIENCE平台的仿真驱动设计方法。通过CATIA构建高精度数字样机，利用SIMULIA Abaqus进行强度与刚度虚拟试验，结合Tosca结构优化与Isight多学科集成，形成“虚拟试验-刚度匹配-拓扑/尺寸优化-验证”闭环。方案在满足极限载荷安全系数1.5、弯曲/扭转刚度约束的前提下，实现主要承力结构减重22%以上，同时将物理试验轮次压缩60%，显著压缩研发周期与成本。
1 引言
复合翼无人机通过固定翼与多旋翼的结合，兼具长航时巡航与垂直起降能力，在物流、巡检、应急测绘等领域得到广泛应用。然而，其结构需同时承受巡航气动弯扭、垂起升力组件集中载荷、模态转换冲击及粗暴着陆等复杂工况，既要保证足够强度与刚度防止颤振和结构屈服，又要极致轻量化以保证航程和载荷。传统“设计-制造-试验-修改”的串行模式周期长、成本高，且无法穷举多工况下的安全边界。
达索系统（Dassault Systèmes）提供的仿真驱动方案，将CATIA三维设计、SIMULIA先进有限元分析、Tosca非参数优化与Isight流程自动化统一在3DEXPERIENCE平台下，使得在数字空间完成全流程虚拟试验和自动化减重优化成为可能。本文详述该方案的技术框架、虚拟试验标准及优化流程。
2 复合翼无人机结构特征与载荷环境
典型复合翼无人机采用上单翼/中单翼布局，四旋翼臂通过管夹或快拆机构连接机身或机翼，尾部配平升力桨。主要结构组件包括：主翼、平尾/垂尾、旋翼臂杆、机身框架、电机座及起落架。材料体系以碳纤维/玻璃纤维增强树脂基复合材料为主，局部使用铝合金、工程塑料等。
关键载荷工况：
-
极限巡航对称机动：最大气动升力状态下机翼受弯、剪及显著扭转，机翼蒙皮和梁承受面外变形。
-
垂直起降最大推力：旋翼臂承受来自电机的拉力与扭矩，引起臂杆弯曲与连接区高应力集中。
-
模态转换/突风：气动与惯性载荷叠加，机翼承受额外动态弯矩。
-
粗暴着陆/跌落：以1.5~2.5 m/s下沉速度模拟起落架触地及机身惯性过载，考核结构抗冲击与剩余强度。
-
颤振与刚度约束：机翼一阶弯曲频率需远离旋翼激励频率，扭转刚度需保证舵面效率及发散安全裕度。
上述载荷工况决定了结构必须同时满足静强度、屈曲稳定性、疲劳（视设计寿命）和刚度多维度约束。
3 达索仿真驱动平台与虚拟试验体系
方案以3DEXPERIENCE平台为统一数据源，实现CAD/CAE全关联和设计变更自动同步。
-
CATIA机械与复合材料设计：实现铺层定义、几何内核无缝衔接，确保后续有限元模型精确表征层合板各向异性及芯层特性。
-
SIMULIA Abaqus：主求解器，执行线性静力、非线性接触/材料、屈曲（线性特征值与后屈曲）、模态频率响应及显式冲击分析，覆盖全部虚拟试验需求。
-
fe-safe：结合实测载荷谱完成复合材料与金属连接区的疲劳寿命评估。
-
Tosca结构优化：在给定载荷下开展拓扑优化（概念阶段）与尺寸/形状优化（详细阶段），直接输出满足制造约束的CAD模型。
-
Isight：集成载荷求解器（自研或开源）进行气动-结构耦合迭代，实现变载荷下的优化与自动化重分析。
虚拟试验框架：数字样机依次通过23个标准分析工况，包括对称拉起、俯冲、滚转、四旋翼悬停、对角电机失效、粗暴着陆等。对应每个工况施加BC（惯性释放或约束）及载荷，输出应力、应变、位移、层合板Tsai-Wu失效指数、紧固件载荷等，并自动出具虚拟试验报告与裕度云图。
4 强度与刚度虚拟试验方案
4.1 有限元模型建立
机翼蒙皮、翼肋与梁采用连续壳或常规壳单元模拟，通过复合材料铺层分区精准赋予0°/±45°/90°铺层信息；旋翼臂为碳纤维管材，采用壳或薄壁梁；连接件采用实体单元并考虑界面接触与螺栓预紧力。整机模型规模约80万~150万节点，利用Abaqus的分布式并行求解。
4.2 静强度评估标准
依据适航借鉴标准，极限载荷下结构不得出现破坏，要求安全系数≥1.5（无人飞行器建议≥1.5~2.0）。以Tsai-Wu准则判断首层失效，以渐进损伤分析确定极限承载。对于金属件，以Von Mises应力不大于材料屈服强度为判据；连接部位校核挤压、剪切裕度。给出各部件剩余强度系数云图，并识别临界工况。分析表明，旋翼臂与机身接头在“单臂失效+侧风悬停”工况应力集中，初始方案局部Tsai-Wu指数超过0.95，需局部加强或重新铺层。
4.3 刚度与模态约束
-
弯曲刚度：机翼在最大起飞重量1.5倍气动载荷下，翼尖变形应不超过半展长的5%，以维持气动外形且避免操纵反效。
-
扭转刚度：通过施加扭矩计算翼尖扭转角，要求展长扭转刚度使副翼效率保持不失速要求，并校核发散速压大于极限速压1.2倍。
-
模态振型：整机前六阶固有频率及振型，确保一阶垂向弯曲频率避开旋翼1P/2P激励带；一阶扭转频率高于弯曲以避免耦合颤振。由Abaqus频率分析可直接评估，必要时调整铺层偏角比例。
4.4 屈曲与后屈曲
对机翼上壁板、机身薄壁段施加压缩/剪切载荷，进行特征值屈曲分析获得临界载荷，并以非线性Riks分析引入初始缺陷，观察后屈曲路径，确保在极限载荷下不发生失稳穿透。
5 减重优化方案
5.1 拓扑优化——承力路径重塑
在Tosca中定义包络设计空间，给予旋翼臂接头区、机身桁架等非层合板区域较大的材料分布自由度。设置综合多工况的加权柔度最小目标，同时施加应力、位移、频率和对称/拔模制造约束。迭代约40~60轮后，获得清晰传力路径的桁架-蒙皮混合概念，相对于初始等值壁厚设计，结构体积下降35%。将该网格重新转为CATIA曲面实体，实现设计链接。
5.2 复合材料铺层尺寸与顺序优化
对机翼蒙皮、梁缘条等层合板结构，在Tosca尺寸优化中，以各区铺层厚度和角度百分比为设计变量，以所有工况Tsai-Wu指数≤1.0、翼尖挠度≤25 mm、一阶弯曲频率≥8 Hz等为约束，目标质量最小。优化后保留主要0°承载铺层，45°铺层调至面外剪切较大区域，去除低效90°铺层。经过灵敏度分析和逐步迭代，机翼结构质量减轻18.7%，翼尖挠度略有增加但仍满足约束。
5.3 夹芯与材料混杂优化
对舵面等次承力结构，采用PMI泡沫夹芯/碳纤维面板，以Tosca协同优化芯层高度和面板厚度，平衡弯曲刚度和质量。某些区域引入玻璃纤维替代碳纤维以降低脆性并减重，这些离散选择通过Isight包装的DOE+响应面进行寻优。
5.4 综合优化结果
优化后整机结构质量从初始14.2 kg降至11.0 kg，减重22.5%，同时全工况剩余强度系数均大于1.1，机翼一阶弯频提升至9 Hz，满足动力学解耦要求。旋翼臂连接区采用拓扑优化镂空并局部加厚耳片后，应力集中系数降低40%。
6 多学科集成与验证
为考虑气动载荷随结构变形重分配的影响，利用Isight搭建流固耦合迭代环：将Abaqus中变形后的几何交予自研涡格法或CFD工具重新计算压力分布，再反馈回结构分析。经3次耦合迭代收敛后，翼梢载荷有所下降，精细化后的应力再校核表明裕度进一步提升。改进后方案通过全部虚拟试验，最终进行关键部件静力物理抽检试验，应变/位移与仿真偏差小于8%，验证了虚拟试验的有效性。
7 结论
本文提出的基于达索仿真驱动的复合翼无人机结构虚拟试验与减重优化方案，将CATIA数字化样机、Abaqus多工况虚拟试验、Tosca拓扑与铺层优化、Isight工作流深度整合，形成覆盖强度、刚度、屈曲、频率及轻量化的完整闭环。方案可在无物理样机阶段预判失效风险，并通过系统级优化挖掘减重潜力超过20%，大幅降低研发成本和周期。该方案可推广至各类轻质无人机、eVTOL飞行器结构开发，是未来智能无人航空器研发的重要方法论。
