4. 铺层优化中的目标函数

铺层优化的本质,是在给定材料体系、几何边界与载荷条件下,寻找一组设计变量(各铺层角度、层数、铺层顺序、总厚度等),使某个或某几个性能指标达到最优。这个被优化的性能指标就是目标函数。目标函数决定了优化算法“往哪个方向走”,也决定了最终方案是更轻、更刚、更强,还是更抗屈曲。工程中很少只追求单一目标,因此理解各目标函数的物理含义、数学表达与相互冲突,是铺层优化设计的第一道关口。

四类典型目标函数的物理内涵

重量最小化是最直观的目标。对于给定铺层面积 A 和单层厚度 t,总重量 W = ρ·A·N·t,其中 ρ 为材料密度,N 为总层数。由于 ρ、A、t 通常固定,重量最小化在数学上等价于层数最小化。它的价值在于减重直接带来油耗、载荷与成本的下降,但单独追求重量最小会迅速触碰刚度、强度与稳定性的下限,因此它几乎总是与约束条件成对出现。

刚度最大化通常以层合板的面内刚度或弯曲刚度为目标。以对称层合板为例,其面内刚度由 A 矩阵描述,弯曲刚度由 D 矩阵描述。工程上常把目标写成最大化某个方向的等效工程常数,例如最大化 Ex,或最大化弯曲刚度 D11。由于 D11 与铺层到中面的距离平方相关,把 0° 层布置在远离中面的位置,能显著提高抗弯刚度——这正是刚度目标对铺层顺序敏感的原因。

强度最大化关注的是层合板在载荷下不失效的能力。常用判据包括最大应力、Tsai-Hill、Tsai-Wu 等。目标函数可写成最大化最小强度比 Rmin,即所有单层中安全裕度最小的那一层的强度比。强度目标对铺层角度分布极为敏感:90° 层主要承担横向应力,±45° 层承担剪切,0° 层承担主方向拉压。若载荷以轴向拉伸为主,增加 0° 层比例通常能提高强度;若以剪切为主,则需提高 ±45° 层比例。

屈曲载荷最大化针对受压或受剪的薄壁层合板。临界屈曲载荷 Ncr 与弯曲刚度 D 矩阵密切相关,近似有 Ncr ∝ √(D11D22)。因此提高屈曲载荷的关键,同样是让高刚度铺层远离中面,并保持铺层关于中面对称以消除拉弯耦合。屈曲目标与刚度目标方向一致,但与重量最小化直接冲突:增加层数能提高屈曲载荷,却增加重量。

铺层优化目标函数关系图 设计变量 → 目标函数 → 约束条件 → 权衡决策 设计变量 铺层角度 θᵢ 各角度层数 nᵢ 铺层顺序 总厚度 / 总层数 N (离散 + 连续混合变量) 目标函数(可单可多) ① 重量最小化 min W = ρANt ② 刚度最大化 max Eₓ / D₁₁ ③ 强度最大化 max R_min ④ 屈曲载荷最大化 max N_cr 约束条件 强度比 R ≥ 1(不失效) 屈曲安全系数 λ ≥ λ₀ 制造约束(对称/平衡) 工艺约束(最小连续层数) (约束把可行域限制在安全区) 目标之间的冲突与协同关系 重量 ↔ 刚度/强度/屈曲 增层提高性能,但增重 刚度 ↔ 屈曲载荷 提高 D 矩阵,两者同向受益 强度 ↔ 刚度 角度分布需求不同,需折中 多目标权衡:加权和 / Pareto 前沿 / 分层序列 在冲突中寻找满足工程约束的最优折中方案

多目标权衡的数学处理方式

当同时考虑重量、刚度、强度和屈曲时,问题变成多目标优化。最常用的处理方式是加权和法:把各目标归一化后乘以权重再相加,例如

min F = w₁·(W/W₀) + w₂·(1 - Eₓ/Eₓ₀) + w₃·(1 - R_min/R₀) + w₄·(1 - N_cr/N_cr₀)
其中 w₁ + w₂ + w₃ + w₄ = 1,W₀、Eₓ₀、R₀、N_cr₀ 为参考值

权重反映设计意图:若以减重为首要任务,可取 w₁ = 0.5,其余权重分配较小;若结构以稳定性控制,则提高 w₄。加权和法简单高效,但权重选择带有主观性,且无法获得非凸 Pareto 前沿上的全部解。

更严谨的方式是求 Pareto 前沿:一组互不支配的解,即在不牺牲其他目标的前提下无法再改进任一目标。工程上常通过改变权重多次求解,或采用 NSGA-II 等进化算法直接逼近前沿,再根据任务需求从前沿上挑选折中方案。

另一种实用策略是分层序列法:按重要性排序,先优化第一目标,再在其允许偏差范围内优化第二目标。例如先以重量最小为目标求得 Wmin,再在 W ≤ 1.05Wmin 的约束下最大化屈曲载荷。这种方法逻辑清晰,适合目标优先级明确的工程场景。

一个可复现的权衡示例

假设一块 300 mm × 200 mm 的碳纤维/环氧层合板,单层厚度 0.125 mm,材料 E₁ = 135 GPa、E₂ = 9 GPa、G₁₂ = 5 GPa、ρ = 1.6 g/cm³,承受轴向压缩载荷。设计变量为 0°、±45°、90° 四种铺层角的层数,总层数不超过 32 层,要求对称铺层。取三组权重方案进行对比:

方案 权重 (w₁,w₂,w₃,w₄) 优化后铺层比例 总层数 相对重量 屈曲载荷比
A(偏减重) 0.55, 0.15, 0.15, 0.15 0°:50%, ±45°:30%, 90°:20% 20 1.00 1.00
B(均衡) 0.25, 0.25, 0.25, 0.25 0°:40%, ±45°:40%, 90°:20% 26 1.30 1.42
C(偏屈曲) 0.10, 0.20, 0.20, 0.50 0°:35%, ±45°:45%, 90°:20% 32 1.60 1.85

表中数据为教学示例的合理假设值,用于说明权重变化对铺层比例和性能的影响趋势,并非实测结果。可以看到:偏重减重时,0° 层比例较高以维持轴向刚度,总层数最少;偏重屈曲时,总层数增加,且 ±45° 层比例上升以提升剪切刚度与 D₂₂,屈曲载荷比显著提高,但重量代价明显。实际工程中,应根据结构的主控失效模式选择权重,并通过参数扫描确认权重微小变化不会导致方案剧烈跳变。

工程建议:先做单目标优化,摸清各目标单独作用时的最优解与性能边界;再以单目标结果为参考,设定多目标权重或约束范围。这样既能避免权重盲目取值,也能快速判断哪些目标之间存在强冲突。
注意:目标函数中的刚度、强度、屈曲量纲不同,直接相加会导致数值大的目标主导优化方向。务必先做归一化处理,例如除以初始设计的对应值或单目标最优值,使各目标处于同一量级。

总结来看,重量、刚度、强度与屈曲四类目标函数分别对应不同的物理机制和数学表达,它们之间既有冲突也有协同。铺层优化的核心工作,就是根据结构的主控载荷与失效模式,合理选择目标组合与权重,在约束可行域内找到性能与重量的最佳平衡点。