SolidWorks Simulation复合材料——层合板建模与Tsai-Wu/Hashin失效准则
一、复合材料力学基础与层合板理论
先进复合材料——尤其是碳纤维增强聚合物（CFRP）——在航空航天、赛车和高端消费电子领域的应用正以每年8%~12%的速度增长。与各向同性金属材料不同，复合材料单层的力学性能具有显著的方向性：沿纤维方向（0°）的拉伸强度可达2000~3000MPa，而横向（90°）仅有50~80MPa，两个方向刚度相差可达10~20倍。这种极端的各向异性既是复合材料可设计性的来源，也是仿真分析的难点所在。
复合材料结构的基本单元是单层（Lamina），其面内刚度由四个独立的工程常数定义：纵向弹性模量E1、横向弹性模量E2、面内剪切模量G12和主泊松比ν12。将多个不同铺层角度的单层按预定顺序堆叠并固化成形，就构成了层合板（Laminate）。铺层顺序的表达式如[45/-45/0/90/0/-45/45]s，其中s表示对称（Symmetric），总层数14层——这种对称铺层是最常用的设计形式，可避免固化降温过程中因热膨胀失配引起的翘曲变形。
经典层合板理论（Classical Laminate Theory, CLT）是复合材料力学分析的基石。CLT将层合板的力学行为统一描述为ABD矩阵：A矩阵（面内刚度）将面内合力与中面应变关联，B矩阵（耦合刚度）将面内合力与弯曲曲率关联（对称铺层使B=0），D矩阵（弯曲刚度）将弯矩与弯曲曲率关联。ABD矩阵的物理意义是：B矩阵非零意味着施加拉伸载荷会产生弯曲变形，这在非对称铺层中必须谨慎处理。最常见的[0/±45/90]准各向同性铺层可使A矩阵呈现各向同性特征（A11≈A22），但D矩阵仍保持各向异性。
单层偏轴刚度通过坐标转换得到：已知单向板在纤维坐标系下的刚度矩阵Qij，通过转换矩阵T(θ)计算偏轴刚度Q̄ij = T⁻¹QijT⁻ᵀ。这一运算贯穿层合板分析的始终——每一层的偏轴刚度矩阵贡献到整体ABD矩阵时都需基于该层的铺层角θk单独计算。
二、SolidWorks Simulation复合材料模块
SolidWorks Simulation的复合材料模块建立在壳单元架构之上，专为薄壁层合板结构设计。首先在材料库中定义单层的正交各向异性（Orthotropic）属性，包括三个方向的弹性模量E1/E2/E3、三个剪切模量G12/G23/G31、三个泊松比ν12/ν23/ν31，以及热膨胀系数α1/α2。对于T700碳纤维/环氧体系，典型值为E1=135GPa、E2=9GPa、G12=5GPa、ν12=0.3。
铺层定义通过层合板编辑器完成：逐层输入材料、铺层角度和厚度。一个典型的对称铺层[45/-45/0/90]s共8层，每层0.125mm则总厚1.0mm。求解器自动对每一层在积分点处完成平面应力计算，后处理可逐层查看应力状态——σ1（纤维方向）、σ2（横向）、τ12（面内剪切）——这是理解层合板破坏机制的核心数据。
对于需要捕捉层间应力（τ13/τ23）的厚壁或弯曲剧烈区域，实体分层建模可替代壳单元方案。实体分层模型使用一层实体单元代表一个物理单层，可在厚度方向直接输出层间剪切应力，但网格数量为壳模型的N倍（N=层数），仅在层间分层风险评估时使用。
后处理中系统自动输出三种失效度量：失效指数FI（Failure Index），FI=1表示临界状态，<1安全，>1失效；安全系数MOS（Margin of Safety），MOS=1/FI-1，正值表示还有安全裕度；强度比SR（Strength Ratio），SR=1/FI，表示当前载荷相对失效载荷的倍数。这些指标按FPF（First Ply Failure，首层失效）和LPF（Last Ply Failure，末层失效）两个层次分别报告——FPF对应最薄弱层的初始失效，LPF对应所有层依次失效后的极限承载能力，两者之间的载荷区间是复合材料的渐进损伤阶段。
三、失效准则的物理内涵与选择
最简单的失效准则——最大应力和最大应变准则——将每个应力/应变分量与对应强度值逐一比较。虽然计算简单，但完全忽略了应力分量之间的交互作用，对多轴应力状态的预测精度较差，现已较少用于正式设计评估。
Tsai-Hill准则是基于von Mises各向同性屈服准则向各向异性材料的推广，通过一个二次型表达所有应力分量的综合效应。但它不区分拉伸和压缩强度——碳纤维的纵向压缩强度（~1200MPa）通常低于拉伸强度（~2500MPa），这种拉压不对称性在Tsai-Hill中被忽略，导致压缩工况的偏不安全预测。
Tsai-Wu张量多项式准则解决了Tsai-Hill的拉压不对称问题，其完整形式为：F1σ1+F2σ2+F11σ1²+F22σ2²+F66τ12²+2F12σ1σ2 ≥ 1。强度参数F1/F2体现了拉压不对称（F1=1/Xt-1/Xc），F11/F22/F66与各方向强度平方相关，而交互系数F12是Tsai-Wu准则中最关键的参数——工程上常采用Mises类比F12≈-0.5√(F11F22)，或通过双轴试验标定。Tsai-Wu准则以其对拉压不对称性的完整处理和良好的工程精度，成为SolidWorks Simulation中默认推荐的通用失效准则。
Hashin准则将失效模式物理分离为四种：纤维拉伸（σ1≥0, (σ1/Xt)²+(τ12/S12)²≥1）、纤维压缩（σ1<0, |σ1|/Xc≥1）、基体拉伸（σ2≥0, (σ2/Yt)²+(τ12/S12)²≥1）、基体压缩（σ2<0, (σ2/2S23)²+[(Yc/2S23)²-1](σ2/Yc)+(τ12/S12)²≥1）。这种按物理机制分类的方法不仅给出了是否失效，还揭示了失效的根源机制，对改进设计有直接的指导意义——例如Hashin基体拉伸主导的失效可通过增加90°层比例来改善。
渐进损伤分析（Progressive Damage Analysis, PDA）超越了简单的"失效/不失效"二元判断。一旦Hashin检测到某层某种模式失效，系统对该层相应方向的刚度进行折减（拉伸失效残留刚度约10%~20%，压缩失效约30%~40%），然后重新分配载荷继续求解。这一迭代过程模拟了层合板的渐进式破坏路径，最终得到的LPF载荷比FPF载荷高出20%~50%不等。层间分层则需引入Cohesive Zone Model（CZM），在相邻层之间插入零厚度的粘聚区单元，其损伤起始和演化由基于能量释放率的B-K准则或Power Law控制。
四、实战案例：无人机碳纤维机翼蒙皮设计
某型固定翼无人机机翼上蒙皮采用T700碳纤维/环氧预浸料，初始铺层方案为[45/-45/0/90/0/-45/45]s，共14层，单层厚度0.125mm，总厚1.75mm。载荷工况为极限机动气动载荷（6g过载）叠加翼尖挠曲，蒙皮承受面内拉压和弯曲的组合作用。
Tsai-Wu分析显示：在6g载荷下，最大失效指数FI=0.82，出现在翼根区域上表面第5层（0°层），虽小于临界值1.0但裕度有限（安全系数MOS=0.22）。切换到Hashin准则后，基体拉伸失效指数FT=0.91是最危险模式（纤维拉伸FT=0.63，基体压缩FC=0.48），表明横向拉应力是制约设计的瓶颈。需要特别指出：FI<1不代表"完全安全"——对于需要取B基准许用值（B-basis，90%置信度/95%可靠度）的航空结构，材料的统计分散性通常要求FI≤0.7，因此FI=0.82在航空级标准下是不通过的。
优化分析从两个方向进行：调整0°层比例从初始的4层（28.6%）增加到5层（35.7%），减少±45°层相应调整以维持总厚度。优化后Tsai-Wu FI从0.82降至0.64——安全裕度提升22%，同时因减少了部分非关键方向的冗余铺层，蒙皮总重反而降低了8%。Hashin各模式失效指数也全面下降，基体拉伸FI从0.91降至0.67。
实验验证方面，制备了优化前后两种铺层方案的层合板试样，进行4点弯曲测试模拟翼面弯曲工况。实验失效载荷与仿真预测（基于渐进损伤+Hashin准则）之间的误差分别为初始方案5.3%和优化方案4.8%，验证了建模方法和失效准则选择的正确性。最终优化的蒙皮方案通过了无人机全尺寸静力试验，极限载荷下未出现任何可视损伤，满足适航审定要求。
结论：复合材料仿真成功的关键在于三个"正确"——正确的单层材料数据（特别关注拉压不对称性和统计分散性）、正确的失效准则选择（Tsai-Wu用于初步筛选，Hashin用于失效模式诊断和设计改进）、以及正确的铺层优化逻辑（通过调整0°层比例在强度和重量之间取得平衡）。
|
|
|
|
|
|
|
|
|
