4. 复合材料失效准则与约束

复合材料铺层优化的目标，不只是让结构“不坏”，而是让每一层的应力状态在给定载荷下都落在可接受的安全边界内。OptiStruct 在复合材料优化中提供了多种失效准则，用于把层内应力转换成可比较的失效指数（Failure Index, FI），再通过约束把 FI 控制在 1 以下。理解这些准则的数学形式、适用场景以及 OptiStruct 中的输出方式，是正确设置约束和解读结果的前提。

一、失效准则的物理含义与数学形式

复合材料单层在面内通常处于平面应力状态，应力分量包括沿纤维方向的 σ₁、垂直纤维方向的 σ₂ 以及面内剪切 τ₁₂。不同失效准则对这些分量的组合方式不同，因而对同一应力状态会给出不同的失效指数。

Tsai-Hill 准则是最早被广泛使用的二次准则之一，其表达式为：

FI = (σ₁/X)² - (σ₁σ₂/X²) + (σ₂/Y)² + (τ₁₂/S)²

其中 X、Y、S 分别为纵向强度、横向强度和面内剪切强度。该准则假设拉压强度相同，因此对拉压性能差异明显的材料会带来偏差。它的优点是形式简单、参数少，适合初步筛选。

Tsai-Wu 准则在 Tsai-Hill 基础上引入了线性项，能够区分拉伸和压缩强度的不同：

FI = F₁σ₁ + F₂σ₂ + F₁₁σ₁² + F₂₂σ₂² + F₆₆τ₁₂² + 2F₁₂σ₁σ₂

其中 F₁、F₂ 由拉压强度差决定，F₁₁、F₂₂、F₆₆ 由各方向强度决定，F₁₂ 为交互项，通常通过双轴试验或经验公式确定。Tsai-Wu 对拉压不对称材料更准确，但 F₁₂ 的取值会显著影响结果，工程中常取 F₁₂ = -0.5√(F₁₁F₂₂) 作为近似。

Hashin 准则则采用分模式判断，把失效区分为纤维拉伸、纤维压缩、基体拉伸、基体压缩四种模式，每种模式有独立的表达式。例如纤维拉伸失效可写为：

FI_ft = (σ₁/X_T)² + (τ₁₂/S)²  (σ₁ ≥ 0)

基体拉伸失效为：

FI_mt = (σ₂/Y_T)² + (τ₁₂/S)²  (σ₂ ≥ 0)

Hashin 的优势在于能指出“哪一层、以什么模式失效”，对铺层优化中的铺层顺序调整和局部加强更有指导意义。OptiStruct 支持上述准则，并允许用户按分析需求选择。

二、OptiStruct 中的实现方式

在 OptiStruct 中，失效准则的指定主要通过 PCOMP 或 PCOMPG 卡片的 FT 字段完成。FT 字段是一个整数选项，不同数值对应不同准则。例如 FT=1 对应 Tsai-Hill，FT=2 对应 Tsai-Wu，FT=3 对应 Hashin，FT=4 对应最大应力，FT=5 对应最大应变。材料强度参数则在 MAT8 卡片中定义，包括 XT、XC、YT、YC、S 等字段。

一个典型的 PCOMP 定义片段如下：

PCOMP   1       0.0     0.0     1       0.0     FT=3
+       0.125   0.0     1       0.125   45.0   1
+       0.125   -45.0   1       0.125   90.0   1

其中 FT=3 表示该铺层使用 Hashin 准则。MAT8 中需给出对应的强度值：

MAT8    1       150000.0        9000.0  0.3     0.0
+       1500.0  1200.0  50.0    200.0   80.0

这里 XT=1500 MPa，XC=1200 MPa，YT=50 MPa，YC=200 MPa，S=80 MPa。注意 MAT8 的字段顺序在不同版本中可能略有差异，应以求解器文档为准。

如果需要在优化中把失效指数作为约束，通常先用 DRESP1 定义响应，选择 FAILURE 类型，并指定层号或全部层；再用 DCONSTR 设置上限为 1.0。这样优化器会在迭代中自动调整铺层厚度或角度，使所有层的 FI 不超过 1。

三、失效指数的输出与解读

OptiStruct 求解完成后，失效指数会出现在 .out 文件的 FAILURE INDEX 表中，按单元和层号列出。同时，HyperView 中可以查看 CFAILURE 结果云图，直观显示高 FI 区域。对于 Hashin 准则，还可以分别查看纤维失效和基体失效的指数，便于判断失效模式。

解读 FI 时需要注意几点。第一，FI 是无量纲比值，FI=1 表示刚好达到准则定义的失效边界，FI<1 表示安全，FI>1 表示失效。第二，不同准则的 FI 数值不可直接横向比较，例如同一应力状态下 Tsai-Hill 的 FI 可能为 0.8，而 Tsai-Wu 的 FI 可能为 0.9，这并不矛盾，因为两者对拉压差异的处理不同。第三，FI 接近 1 的区域往往是优化中最需要关注的区域，可能需要增加铺层厚度或调整铺层角度。

四、准则选择与工程判断

三种准则各有适用场景。Tsai-Hill 适合快速筛选和初步设计，参数少、计算快；Tsai-Wu 适合拉压性能差异明显的材料，如碳纤维复合材料；Hashin 适合需要区分失效模式的详细分析，尤其是铺层优化中需要针对性调整铺层顺序和角度时。

在实际工程中，还可以把不同准则的结果进行对比。如果 Tsai-Hill 和 Tsai-Wu 给出的高 FI 区域一致，说明该区域确实危险；如果两者差异较大，则需要检查材料强度参数是否准确，尤其是拉压强度差异是否被正确输入。

最后，失效指数的输出应与铺层优化目标结合。OptiStruct 的铺层优化通常以质量最小为目标，以 FI ≤ 1 为约束。优化完成后，应检查最终方案的 FI 分布是否均匀，避免局部 FI 过高导致安全裕度不足。如果 FI 最大值远小于 1，说明约束过松，可以进一步减重；如果 FI 最大值接近 1，则说明约束起控制作用，需要确认该区域是否允许接近极限状态。