3. 层合板理论及Abaqus铺层建模

复合材料层合板由多层单向带或织物按不同角度叠合而成，其宏观刚度不是各层刚度的简单平均，而是与每层的材料主方向、厚度位置和铺层顺序密切相关。经典层合板理论（Classical Laminate Theory, CLT）正是把三维问题降为二维中面问题的一套工程方法：它假设层间粘结完好、变形沿厚度方向线性分布、忽略层间正应力，从而把面内载荷与弯曲载荷统一到中面应变和曲率上。理解CLT和ABD矩阵，是后续在Abaqus中正确建立Composite Layup、解读层间应力和判断失效位置的前提。

一、从单层到层合板：CLT的基本假设与位移场

CLT建立在Kirchhoff薄板假设之上，核心位移场为：

u(x,y,z) = u0(x,y) - z · ∂w0/∂x
v(x,y,z) = v0(x,y) - z · ∂w0/∂y
w(x,y,z) = w0(x,y)

其中u0、v0、w0为中面位移，z为厚度方向坐标，原点取在中面。由此得到面内应变与中面应变ε0、曲率κ的关系：ε = ε0 + z·κ。这一线性关系意味着：同一层内应变沿厚度线性变化，而层与层之间由于铺层角不同，应力会出现跳跃。这正是层合板分析中“应变连续、应力不连续”的典型特征。

需要特别说明坐标系约定：本文及后续Abaqus建模中，1方向为纤维方向，2方向为面内垂直纤维方向，3方向为厚度方向；铺层角θ定义为纤维方向与结构参考轴（通常为壳单元局部1轴或实体坐标系x轴）之间的夹角，逆时针为正。全文与图示保持这一约定。

二、ABD矩阵：层合板刚度的完整表达

将各层应力沿厚度积分，可得到层合板的内力—应变关系：

[ N ]   [ A  B ] [ ε0 ]
[ M ] = [ B  D ] [ κ  ]

其中N为面内合力，M为弯矩合力。A矩阵是拉伸刚度，反映中面应变产生的面内力；D矩阵是弯曲刚度，反映曲率产生的弯矩；B矩阵是耦合刚度，反映拉伸与弯曲之间的耦合。B矩阵是否为零，直接决定层合板是否会出现“拉弯耦合”。

A、B、D各元素由每层的偏轴刚度Q̄和厚度位置计算：

A_ij = Σ Q̄_ij^(k) · (z_k - z_{k-1})
B_ij = 1/2 · Σ Q̄_ij^(k) · (z_k² - z_{k-1}²)
D_ij = 1/3 · Σ Q̄_ij^(k) · (z_k³ - z_{k-1}³)

从公式可以看出：A只与各层厚度有关，与铺层顺序无关；B和D则依赖每层到中面的距离，因此铺层顺序会显著改变B和D。这就是为什么同样材料、同样角度的层合板，仅改变叠放顺序，其弯曲刚度和耦合行为就可能完全不同。

三、铺层角度与顺序的工程影响

铺层角度决定各层偏轴刚度Q̄的数值。以T300/环氧单向带为例，典型工程参数为E1=135 GPa、E2=9 GPa、G12=5 GPa、ν12=0.3。当铺层角从0°转到90°时，Q̄11大幅下降、Q̄22上升，层合板整体拉伸刚度随之改变。若采用[0/±45/90]s这类准各向同性铺层，面内刚度接近各向同性，但层间应力分布仍不均匀。

铺层顺序的影响更隐蔽。例如[0₂/90₂]与[90₂/0₂]的A矩阵相同，但B、D不同；而[0/90]s与[0/90/0/90]虽然角度序列相似，对称性不同，B矩阵一个为零、一个非零。实际设计中，常把±45°层放在外侧以提高抗剪和抗冲击能力，把0°层放在靠近中面以承担主要轴向载荷，同时避免同一角度连续层数过多，以降低分层风险。

四、Abaqus Composite Layup 模块建模流程

在Abaqus/CAE中，Composite Layup用于为壳单元或实体单元定义多层铺层。壳铺层通过截面属性定义，实体铺层则通过分配材料方向并逐层划分网格实现。以下以壳铺层为例说明关键操作。

首先创建材料与截面。在Property模块中定义T300/环氧的弹性常数：E1=135000 MPa、E2=9000 MPa、G12=5000 MPa、G13=5000 MPa、G23=3000 MPa、ν12=0.3。随后创建Composite Section，选择Shell类型，进入Composite Layup对话框。

在Layup编辑区逐行输入铺层：Region选择壳面，Material选择已定义材料，Thickness输入单层厚度（如0.125 mm），Orientation Angle输入铺层角（0、45、-45、90），Integration Point选择3或5。Abaqus会按输入顺序从下往上堆叠，第一行对应最底层。若需要对称铺层，可勾选Symmetrize，软件自动镜像生成对称部分。

示例铺层输入（从底层到顶层）：
Layer 1: 0°   , 0.125 mm
Layer 2: 45°  , 0.125 mm
Layer 3: -45° , 0.125 mm
Layer 4: 90°  , 0.125 mm
Layer 5: 90°  , 0.125 mm
Layer 6: -45° , 0.125 mm
Layer 7: 45°  , 0.125 mm
Layer 8: 0°   , 0.125 mm
总厚度 = 1.0 mm

对于实体铺层，需先对实体划分网格，再在Mesh模块中使用Assign Material Orientation指定每层的局部坐标系，最后通过Composite Layup将各层材料赋给对应单元层。实体铺层能更真实地反映层间应力，但计算量显著增加，通常用于关注分层或层间剪切的局部细节分析。

五、结果判读与验证要点

提交计算后，应重点检查以下内容：一是层间应力分布，通过Section Points查看各层顶/底面的S11、S22、S12，判断哪一层最先达到强度极限；二是中面应变与曲率，验证是否与CLT手算结果一致；三是B矩阵是否为零，若采用对称铺层但结果出现明显拉弯耦合，应检查铺层顺序或对称设置是否正确。

以[0/±45/90]s层合板受单轴拉伸为例，0°层承担主要轴向载荷，其S11最大；90°层横向应力S22较高但刚度低；±45°层剪应力S12显著。若在Abaqus中查看各层S11，应看到0°层数值最高、90°层最低，且沿厚度方向呈阶梯状分布。若分布异常，需回到Layup定义检查角度、厚度和堆叠顺序。