2. 复合材料力学基础与Abaqus材料库

复合材料层合板的分析起点，是单层板（lamina）在材料主方向上的三维正交各向异性弹性描述。与各向同性金属只需弹性模量 E 和泊松比 ν 不同，正交各向异性材料在三个互相垂直的材料主方向 1、2、3 上具有不同的刚度，完整描述需要 9 个独立的工程常数：三个方向的弹性模量 E₁、E₂、E₃，三个泊松比 ν₁₂、ν₁₃、ν₂₃，以及三个面内剪切模量 G₁₂、G₁₃、G₂₃。对于常见的单向带单层板，通常约定 1 方向为纤维方向，2 方向为面内垂直纤维方向，3 方向为厚度方向，因此 E₁ 远大于 E₂ 和 E₃，这是复合材料"方向性强"的根本来源。

工程常数与柔度矩阵

在材料主方向坐标系下，正交各向异性材料的应力-应变关系可用柔度矩阵 S 表示，应变列向量等于柔度矩阵乘以应力列向量。柔度矩阵由工程常数直接构成，其非零项为：S₁₁ = 1/E₁，S₂₂ = 1/E₂，S₃₃ = 1/E₃，S₄₄ = 1/G₂₃，S₅₅ = 1/G₁₃，S₆₆ = 1/G₁₂，以及耦合项 S₁₂ = −ν₁₂/E₁ = −ν₂₁/E₂，S₁₃ = −ν₁₃/E₁ = −ν₃₁/E₃，S₂₃ = −ν₂₃/E₂ = −ν₃₂/E₃。这里必须注意泊松比的对称关系：νᵢⱼ/Eᵢ = νⱼᵢ/Eⱼ，因此 9 个常数中只有 3 个泊松比是独立的，另外 3 个由对称关系导出。若在 Abaqus 中输入的泊松比不满足该对称性，求解器会给出警告甚至导致刚度矩阵非正定。

工程常数描述的是"给定应力求应变"的柔度视角，而有限元求解需要的是刚度矩阵。Abaqus 内部会将工程常数转换为刚度矩阵 C = S⁻¹，再组装单元刚度。对于平面应力状态下的单层板（薄板假设 σ₃₃ = 0，且忽略面外剪切），退化为 4 个独立常数：E₁、E₂、ν₁₂、G₁₂，这正是经典层合板理论（CLT）中 Q 矩阵的输入。此时折减刚度 Q₁₁ = E₁/(1−ν₁₂ν₂₁)，Q₂₂ = E₂/(1−ν₁₂ν₂₁)，Q₁₂ = ν₁₂E₂/(1−ν₁₂ν₂₁)，Q₆₆ = G₁₂。

单层板应力-应变关系的坐标变换

上述关系仅在材料主方向 1-2-3 下成立。实际层合板中，单层板的纤维方向与结构整体坐标 x-y 通常存在铺层角 θ，因此需要将主方向刚度矩阵通过应变转换矩阵 T 变换到整体坐标：Q̄ = T⁻¹ Q T⁻ᵀ。变换后会出现耦合项 Q̄₁₆ 和 Q̄₂₆，它们把正应力与剪应变耦合起来，这正是 ±45° 铺层能产生拉剪耦合、而 0°/90° 铺层不产生耦合的原因。理解这一点，对后续判断层合板耦合行为和选择铺层顺序至关重要。

Abaqus 中的两种材料定义方式

Abaqus/CAE 的 Property 模块提供两种与本节直接相关的弹性定义类型，选择哪一种取决于分析对象和维度。

Engineering Constants（工程常数）：适用于三维实体单元或需要完整三维正交各向异性描述的场景。在 Material → Elastic → Type 中选择 Engineering Constants，Type 可选 Isotropic、Orthotropic、Anisotropic。选择 Orthotropic 后，Abaqus 要求输入 E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃，并可指定局部材料方向（通过 Assign Material Orientation 将 1-2-3 方向赋予单元）。若只输入 4 个面内常数，Abaqus 会提示需要补充面外常数，否则无法完成三维刚度组装。

Lamina（单层板）：专为壳单元和平面应力状态设计，本质是平面应力正交各向异性。在 Material → Elastic → Type 中选择 Lamina，输入 E₁、E₂、ν₁₂、G₁₂、G₁₃、G₂₃（其中 G₁₃、G₂₃ 用于横向剪切，壳单元需要）。Lamina 类型隐含 σ₃₃ = 0，因此不需要 E₃、ν₁₃、ν₂₃。它常用于复合材料壳建模，配合 Composite Layup 定义铺层角、厚度和积分点。

实操示例：定义一块 T300/环氧 单层板

以示例参数 E₁ = 135 GPa、E₂ = 9.5 GPa、ν₁₂ = 0.30、G₁₂ = 5.0 GPa、G₁₃ = 5.0 GPa、G₂₃ = 3.3 GPa 为例，在 Abaqus 中建立 Lamina 材料并赋予壳单元：

*Material, name=CFRP_Lamina
*Elastic, type=LAMINA
135000., 9500., 0.30, 5000., 5000., 3300.
*Density
1.6e-09

上述单位为 MPa（模量）、mm（长度）、t/mm³（密度），对应 Abaqus 常用的一致单位制。若改用 Engineering Constants 三维定义，则写成：

*Elastic, type=ENGINEERING CONSTANTS
135000., 9500., 9500., 0.30, 0.30, 0.45, 5000., 5000.
3300.

注意 Engineering Constants 数据行中，前 8 个数为 E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃，第二行的第一个数为 G₂₃。输入完成后，可通过 Section Assignment 指定壳厚度与积分点，再用 Assign Material Orientation 定义 1 方向。若铺层角为 45°，则材料 1 方向相对整体 x 轴旋转 45°，此时 Q̄₁₆ 非零，单层板在面内拉伸下会出现剪应变，这是后续层合板耦合分析的关键观察点。

判断材料定义是否正确，可在 Job 提交前用 Query → Material 检查刚度矩阵是否正定，或在 Datacheck 阶段查看是否有 "THE POISSON'S RATIOS DO NOT SATISFY THE SYMMETRY CONDITION" 类警告。若出现该警告，应回到 νᵢⱼ/Eᵢ = νⱼᵢ/Eⱼ 关系重新核对输入值，而不是忽略警告继续计算。