4. OptiStruct优化三要素

OptiStruct 的优化问题在数学上可以统一写成:在满足约束条件的前提下,寻找一组设计变量,使目标函数达到最优。复合材料优化之所以比金属结构优化复杂,根本原因在于它的“设计变量”不再只是板厚或截面尺寸,而是同时包含铺层角度、铺层顺序、单层厚度、区域划分等离散与连续混合的变量。理解三要素在复合材料语境下的具体表现形式,是正确建立优化模型的前提。

一、设计变量:从“尺寸”扩展到“铺层”

在复合材料优化中,设计变量按抽象层次可分为三类。第一类是尺寸变量,例如层合板总厚度、各铺层组的单层厚度、加筋条的截面尺寸,通常用 DVPREL1/DVPREL2 关联到 PSHELL 或 PCOMP 的属性字段。第二类是形状变量,用于描述铺层边界的移动或铺放区域的形状变化,在自由尺寸优化阶段常以单元密度形式出现。第三类是拓扑/铺层变量,这是复合材料特有的核心变量,包括每个区域的铺层角度(0°、±45°、90°)、铺层顺序以及某角度铺层是否存在。

OptiStruct 通过 PCOMP 卡片与铺层优化卡片(如 PCOMPG、PLY 定义)配合,把每一层的角度和厚度暴露为可优化对象。在自由尺寸优化(Free-size)阶段,设计变量是每个单元在每一铺层方向上的“材料密度”,取值 0~1,代表该处是否需要该方向的铺层;在铺层顺序优化(Shuffle)阶段,设计变量变为各角度铺层在厚度方向上的排列次序;在尺寸优化阶段,设计变量则退化为各角度铺层组的厚度。这种“先自由尺寸、再铺层顺序、最后尺寸微调”的变量分层,是复合材料优化区别于普通尺寸优化的关键。

二、目标函数:刚度、质量与多目标权衡

复合材料优化最常用的目标函数有三类。其一是柔度最小化(等价于刚度最大化),OptiStruct 中通过 MINIMIZE 卡片配合 COMPLIANCE 响应实现,适合以变形控制为主的翼面、舱门类结构。其二是质量最小化,通过 MASS 响应与 MINIMIZE 组合,是减重设计最直接的目标。其三是多目标加权,例如同时考虑柔度和质量,用 DEQATN 定义加权函数,或采用 OptiStruct 的多目标优化策略分步求解。

在复合材料中,目标函数还常与铺层工艺性挂钩。例如把“总铺层数”或“不同角度铺层数量”作为目标,以降低制造成本;或把“层间剪切柔度”作为目标,抑制分层风险。需要注意的是,目标函数必须是设计变量的连续或可微函数(对梯度类算法而言),因此铺层顺序这类离散变量通常先松弛为连续密度,再通过后续的铺层顺序优化还原为离散铺层。

三、约束条件:从应力到工艺的全链条限制

复合材料的约束条件比金属更丰富,可分为四层。第一层是响应约束,包括位移约束(如翼尖最大变形 ≤ 20 mm)、应力/应变约束(各向异性 Tsai-Wu、Tsai-Hill 或最大应变准则)、屈曲因子约束(特征值 ≥ 1.2)、频率约束(一阶频率 ≥ 15 Hz)。第二层是体积/质量约束,如总体积分数 ≤ 0.3,用于控制材料用量。第三层是制造工艺约束,这是复合材料特有的:对称铺层约束、平衡铺层约束(+45°与-45°成对出现)、连续同角度铺层数限制(如最多连续 4 层)、每层最小/最大厚度、铺层比例约束(如 0°铺层占比 30%~50%)。第四层是变量边界约束,如单层厚度 0.1~0.25 mm、角度只能取 0/±45/90。

在 OptiStruct 中,这些约束通过 DRESP1/DRESP2 定义响应,再用 DCONSTR 或 DCONADD 施加。例如用 DRESP1 提取 Tsai-Wu 失效因子,用 DCONSTR 限制其 ≤ 1.0;用 DRESP1 提取体积分数,限制 ≤ 0.3;用 MANUFACTURING 卡片施加对称、平衡、连续层数等工艺约束。约束的合理设置直接决定优化能否收敛到可制造的铺层方案。

复合材料优化三要素关系图 设计变量 → 目标函数 ← 约束条件,三者共同定义优化问题 设计变量 • 铺层角度:0° / ±45° / 90° • 单层厚度:0.1~0.25 mm • 铺层顺序:厚度方向排列 • 区域密度:0~1(自由尺寸) • 加筋截面:长宽高尺寸 卡片:PCOMP / PLY / DVPREL 目标函数 • 柔度最小化(刚度最大) • 质量最小化(减重) • 多目标加权:柔度+质量 • 铺层数最少(工艺成本) • 层间剪切柔度最小 卡片:MINIMIZE / DEQATN 约束条件 • 位移:翼尖 ≤ 20 mm • 失效:Tsai-Wu ≤ 1.0 • 屈曲因子 ≥ 1.2 • 体积分数 ≤ 0.3 • 工艺:对称/平衡/连续层 卡片:DRESP / DCONSTR / MANUFACT 复合材料优化三阶段变量演进 ① 自由尺寸优化 变量:单元密度 0~1 目标:柔度最小 约束:体积分数、应力 ② 铺层顺序优化 变量:铺层排列次序 目标:工艺性最优 约束:对称/平衡/连续层 ③ 尺寸优化 变量:各角度层厚度 目标:质量最小 约束:位移、失效、屈曲 三要素在每一阶段的具体形式不同,但始终围绕“变量—目标—约束”闭环迭代

四、三要素的协同设置方法

实际建模时,三要素不是孤立设置的,而是通过“响应—约束—目标”链条耦合在一起。推荐的操作顺序是:先定义所有需要的响应(DRESP1 提取位移、应力、体积、频率、屈曲因子;DRESP2 组合出 Tsai-Wu 因子或加权函数),再把这些响应分别挂到约束(DCONSTR)和目标(MINIMIZE)上,最后定义设计变量及其与属性的关联(DVPREL 或铺层卡片)。

以一个典型的机翼盒段蒙皮为例:设计变量为各区域 0°、±45°、90° 铺层组的厚度,初始值均取 0.125 mm;目标函数为质量最小(MASS 响应);约束包括翼尖位移 ≤ 20 mm、Tsai-Wu 失效因子 ≤ 1.0、一阶频率 ≥ 15 Hz、体积分数 ≤ 0.35,同时施加对称与平衡铺层约束。求解后若质量下降但屈曲因子低于 1.2,则说明约束设置不完整,需要补充屈曲约束后重新求解。这个“设置—求解—检查—补充约束”的循环,正是三要素协同调整的典型过程。

设置要点:复合材料优化中,约束数量往往多于金属优化。建议先用自由尺寸优化确定铺层分布,再逐步加入工艺约束和失效约束,避免一次性约束过多导致不可行或收敛困难。目标函数优先选择柔度或质量这类全局响应,局部响应(如应力)更适合作为约束而非目标。

五、常见设置错误与判断依据

初学者最容易犯的错误有三类。第一类是变量与属性脱节:定义了 DVPREL 却没有正确指向 PCOMP 的层厚度字段,导致优化迭代中厚度不变。判断方法是检查 .out 文件中设计变量的初始值与迭代值是否变化。第二类是约束过紧:例如把 Tsai-Wu 因子限制到 0.5 以下,导致可行域为空,求解器报“no feasible design”。此时应放宽约束或检查初始设计是否已违反约束。第三类是目标与约束冲突:例如同时要求质量最小和位移极小,二者相互矛盾,表现为目标函数震荡不收敛。解决方法是采用多目标加权或分步优化。

判断优化是否成功的依据包括:目标函数是否单调下降并趋于稳定;约束是否全部满足(.out 中约束违反量接近零);设计变量是否落在合理边界内(如厚度未触及上限或下限);以及最终铺层方案是否满足工艺约束。只有这四项同时成立,才能认为三要素的设置是有效的。

核心结论:设计变量决定“能改什么”,目标函数决定“往哪改”,约束条件决定“改的边界”。在复合材料优化中,三者的具体形式随优化阶段(自由尺寸→铺层顺序→尺寸)动态变化,设置时必须与铺层卡片、响应卡片和工艺约束卡片一一对应,才能得到可制造的最优铺层方案。