2. 铺层设计变量与约束

铺层优化之所以比一般尺寸优化更难，根本原因在于它的设计变量不是连续实数，而是一组带有离散取值、顺序语义和工艺可行性的“铺层序列”。要把优化问题交给算法求解，第一步必须把工程语言翻译成数学语言：哪些量是设计变量，哪些量是约束，约束又如何写成可计算的不等式或等式。本节围绕铺层角度、厚度、顺序、对称性、平衡性、连续性和制造约束，建立一套可直接用于建模的数学表达框架。

一、设计变量的三种基本形式

铺层设计变量通常以三种粒度出现。第一种是单层角度变量：把每一层视为一个变量 θ_k，取值来自离散角度集合，例如工程上常用的 {0°, ±45°, 90°}，即 θ_k ∈ {−45°, 0°, 45°, 90°}。第二种是铺层组变量：把连续同角度的若干层合并为一个铺层组，变量为“角度 + 层数”，例如 [45₂/0₃/90]，这样可显著减少变量维度。第三种是厚度变量：以单层厚度 t_ply 为基准，用层数 n_k 表示厚度，即 t_k = n_k · t_ply，总厚度 T = Σ n_k · t_ply。厚度变量既可以是整数层数，也可以在一定条件下松弛为连续厚度，再在后期圆整回整数层。

需要特别强调的是顺序语义：铺层角度集合相同，但排列顺序不同，层合板的弯曲刚度矩阵 D 会完全不同。因此顺序不是可忽略的附属信息，而是与角度、厚度并列的独立设计维度。在数学上，顺序体现为变量下标 k 的排列，即 θ₁, θ₂, …, θ_N 的序列本身携带信息。

二、对称性与平衡性的数学表达

对称性约束用于消除拉-弯耦合。设总层数为 N（通常为偶数），对称铺层要求第 k 层与第 N+1−k 层角度相同：

θ_k = θ_(N+1-k),  k = 1, 2, …, N/2

满足该条件时，耦合刚度矩阵 B = 0，层合板在面内载荷下不产生弯曲变形。对称性把独立变量数从 N 降到 N/2，是最有效的降维手段之一。

平衡性约束用于消除拉-剪耦合。其数学表达为：对每一个 +θ 层，必须存在一个 −θ 层与之配对，即角度集合中正负角度成对出现：

count(+θ) = count(−θ),  对任意 θ ≠ 0°, 90°

平衡性使拉伸刚度矩阵 A 中的 A₁₆ = A₂₆ = 0，避免面内拉伸与剪切耦合。工程上常把对称与平衡合并使用，形成“对称平衡铺层”，此时 B = 0 且 A₁₆ = A₂₆ = 0，层合板行为最接近各向同性板，便于设计与分析。

三、连续性与制造约束

连续性约束限制同一角度连续出现的层数，防止过厚的同角度铺层组引起基体开裂和分层。设同一角度连续层数为 m_θ，则：

m_θ ≤ m_max,  工程常用 m_max = 4

此外还有几类常见制造约束需要写成数学形式。第一，丢层约束：相邻铺层组之间一次最多丢 1～2 层，避免刚度突变，可写为 |n_k − n_k+1| ≤ Δn_max。第二，边角约束：最外层通常要求为 ±45° 或 0°，以保护纤维端部，可写为 θ₁, θ_N ∈ {−45°, 45°}。第三，10% 规则：任一方向铺层数不少于总层数的 10%，即 n_θ / N ≥ 0.1，防止某一方向刚度过低。第四，总厚度约束：T_min ≤ Σ n_k · t_ply ≤ T_max，同时满足强度与重量指标。

四、约束体系的整体关系

把上述内容放在一起，可以看到设计变量与约束之间并非简单叠加，而是存在明确的层级与耦合关系：角度和厚度是基础变量，顺序决定刚度分布，对称性与平衡性作用于整体耦合项，连续性与制造约束作用于局部可行性。下图给出这一关系结构，帮助在建模时快速定位每条约束的作用位置。

五、建模示例：从工程描述到数学表达

假设某层合板要求：总层数 16，单层厚度 0.125 mm，角度集合 {−45°, 0°, 45°, 90°}，对称平衡，任一方向不少于 10%，同角度连续不超过 4 层，总厚度不超过 2.0 mm。将其翻译为数学模型：

变量：θ_k ∈ {−45°, 0°, 45°, 90°}, k = 1…16
       n_k ∈ Z⁺,  Σ n_k = 16
约束：θ_k = θ_(17−k)                    （对称）
       count(+45°) = count(−45°)        （平衡）
       n_θ / 16 ≥ 0.1,  θ ∈ {0°, 45°, 90°}  （10% 规则）
       m_θ ≤ 4                          （连续性）
       Σ n_k · 0.125 ≤ 2.0              （总厚）
目标：min 质量 或 max 屈曲载荷

该模型可直接交给遗传算法、模拟退火或混合整数规划求解。求解时通常先松弛连续厚度获得下界，再圆整为整数层并校验全部约束，最后用经典层合板理论计算 A、B、D 矩阵，确认对称与平衡条件是否严格满足。

总结来看，铺层设计变量与约束的数学表达，本质上是把工程经验转化为可计算的搜索空间边界。角度、厚度、顺序构成变量空间，对称、平衡、连续、丢层、10% 规则和总厚约束共同界定可行域。只有把每条约束写成明确的等式或不等式，优化算法才能在正确的空间内寻找既满足力学性能又满足制造要求的最优铺层方案。