4. 层合板经典理论（CLT）与Abaqus壳单元建模
经典层合板理论（Classical Lamination Theory, CLT）是复合材料结构分析的基石。它把单层材料在厚度方向上的刚度贡献，通过铺层角度和位置加权，集成为整个层合板的等效刚度矩阵。理解CLT的核心，就是理解ABD矩阵如何描述面内刚度、耦合刚度和弯曲刚度，以及这些量如何在Abaqus壳单元中被正确表达。
ABD矩阵：层合板刚度的完整描述
CLT基于以下假设：层合板由若干正交各向异性单层按指定角度和顺序叠合而成，层间完全粘结，变形满足直法线假设。单层在材料主方向上的刚度矩阵为Q，经坐标变换到层合板整体坐标系后得到Q̄。设第k层距中面的位置为zk-1到zk，则ABD矩阵定义为：
A_ij = Σ Q̄_ij^(k) · (z_k − z_{k−1})
B_ij = ½ Σ Q̄_ij^(k) · (z_k² − z_{k−1}²)
D_ij = ⅓ Σ Q̄_ij^(k) · (z_k³ − z_{k−1}³)其中A为面内刚度矩阵，B为耦合刚度矩阵，D为弯曲刚度矩阵。B矩阵的存在意味着面内变形与弯曲变形相互耦合。对于对称层合板，B=0；对于反对称或非对称层合板，B≠0，拉伸会引起弯曲，弯曲也会引起拉伸。这正是铺层顺序影响结构响应的根本原因。
以常见的T300/环氧单层为例，假设E1=135 GPa，E2=9 GPa，G12=5 GPa，ν12=0.3，单层厚度0.125 mm。若铺层为[0/90]s，共4层，总厚度0.5 mm。由于对称，B=0。A11由两层0°和两层90°贡献，D11则与各层距中面的距离平方相关。若改为[0/90/0/90]，虽然总厚度相同，但B≠0，层合板在拉伸时会产生弯扭耦合，这是工程中需要避免的。
铺层顺序与角度定义
铺层顺序的书写规则为：从层合板底面（z=−h/2）到顶面（z=+h/2）依次列出各层角度，用斜杠分隔。下标s表示对称，即[0/90]s等价于[0/90/90/0]。下标T表示总层数，如[0/90]T表示两层。角度定义为纤维方向与层合板整体坐标系x轴的夹角，逆时针为正。在Abaqus中，这个整体坐标系通常对应壳单元的局部材料坐标系，需要与截面定义中的方向一致。
Abaqus壳单元截面定义与Composite Layup工具
在Abaqus中，复合材料层合板通常用壳单元（如S4R、S8R）模拟。壳单元的截面定义有两种方式：传统Shell Section和Composite Layup。传统方式需要手动输入各层材料、厚度、角度和积分点，操作繁琐且容易出错。Composite Layup工具则提供了图形化界面，可以直观地定义铺层区域、材料、厚度、角度和铺层顺序，并自动生成截面属性。
使用Composite Layup的典型流程如下：
- 在Property模块中创建各单层材料，定义弹性常数（Engineering Constants）和强度准则（如Hashin）。
- 创建Composite Layup，选择壳单元区域，设置铺层坐标系（通常为壳单元局部坐标系）。
- 逐层添加铺层：指定材料、厚度、角度和积分点数量。Abaqus默认每层3个积分点，对于弯曲问题可适当增加。
- 检查铺层顺序和总厚度，确保与设计一致。Composite Layup会实时显示层合板截面示意图。
- 将Composite Layup分配给相应的壳单元区域，完成截面定义。
需要注意的是，Composite Layup中的角度是相对于铺层坐标系定义的。如果壳单元的局部坐标系与整体坐标系不一致，需要先通过Assign Material Orientation或Orientation工具调整。此外，Composite Layup支持变厚度铺层和局部加强，适合复杂结构建模。
建模中的关键检查点
完成Composite Layup定义后，建议进行以下检查：第一，确认铺层顺序与设计图纸一致，特别是底面到顶面的方向；第二，检查总厚度是否等于各层厚度之和；第三，对于对称铺层，确认B矩阵是否为零（可通过Abaqus的截面查询功能查看）；第四，确认材料方向与铺层角度定义一致，避免因坐标系混淆导致角度错误；第五，对于壳单元，注意积分点数量对弯曲刚度的影响，通常每层3个积分点足够，但厚层或高精度弯曲问题可增加到5个。
CLT与Abaqus壳单元建模的关系可以概括为：CLT提供了层合板刚度的理论框架，ABD矩阵是这一框架的数学表达；Abaqus的Composite Layup工具则是将这一理论转化为有限元模型的工程手段。理解ABD矩阵的物理意义，有助于在建模时正确设置铺层顺序和角度，避免因耦合效应导致的非预期变形。对于大多数工程结构，优先采用对称均衡铺层，以消除拉弯耦合和固化翘曲，这是复合材料设计的基本原则。