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 卡片施加对称、平衡、连续层数等工艺约束。约束的合理设置直接决定优化能否收敛到可制造的铺层方案。
四、三要素的协同设置方法
实际建模时,三要素不是孤立设置的,而是通过“响应—约束—目标”链条耦合在一起。推荐的操作顺序是:先定义所有需要的响应(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 中约束违反量接近零);设计变量是否落在合理边界内(如厚度未触及上限或下限);以及最终铺层方案是否满足工艺约束。只有这四项同时成立,才能认为三要素的设置是有效的。