2. 复合材料力学基础与Ansys实现映射
复合材料层合板之所以难算，根本原因在于它的刚度不是"一个数"，而是随方向变化的张量。理解这一点，才能理解为什么 Ansys Mechanical 里要填那么多工程常数、为什么要指定铺层角度、为什么结果里会出现耦合变形。本节把各向异性弹性、正交各向异性工程常数、经典层合板理论三条线索串起来，并落到 Mechanical 的材料数据组织方式上。
一、从各向同性到各向异性：刚度矩阵的自由度
各向同性材料（如钢）的弹性行为只需两个独立常数：弹性模量 E 和泊松比 ν，刚度矩阵 C 中大量元素相等或为零。正交各向异性材料（单向带、织物单层）有三个互相垂直的弹性对称面，独立常数增加到 9 个：三个方向模量 E₁、E₂、E₃，三个泊松比 ν₁₂、ν₁₃、ν₂₃，三个剪切模量 G₁₂、G₁₃、G₂₃。若进一步简化为横观各向同性（纤维在 2-3 平面内随机分布），E₂=E₃、ν₁₂=ν₁₃、G₁₂=G₁₃，独立常数降为 5 个。
工程上最常用的是平面应力假设下的单层板：厚度方向应力 σ₃、τ₁₃、τ₂₃ 近似为零，于是正交各向异性退化为 4 个独立常数 E₁、E₂、ν₁₂、G₁₂。这正是 Mechanical 中"正交各向异性弹性"材料模型在二维壳单元下的实际输入项。注意 ν₂₁ 不是独立输入，它由互等关系确定：ν₂₁ = ν₁₂·E₂/E₁。填错这一条，刚度矩阵就不对称，求解器会报错或给出非物理结果。
二、坐标约定：材料主方向与铺层角
单层板的 1 方向沿纤维，2 方向垂直纤维，3 方向为厚度方向。层合板整体坐标用 x、y、z 表示，铺层角 θ 定义为材料 1 轴与整体 x 轴的夹角，逆时针为正。这个符号约定必须在建模、铺层定义、结果查看三处保持一致，否则 ±45° 铺层的耦合效应会完全反向。
单层刚度从材料主方向转到整体坐标，靠的是应力/应变转轴矩阵 T：
Q11 = Q11m·c^4 + 2(Q12m+2Q66m)·s^2c^2 + Q22m·s^4
Q22 = Q11m·s^4 + 2(Q12m+2Q66m)·s^2c^2 + Q22m·c^4
Q12 = (Q11m+Q22m-4Q66m)·s^2c^2 + Q12m·(c^4+s^4)
Q16 = (Q11m-Q12m-2Q66m)·c^3s - (Q22m-Q12m-2Q66m)·cs^3
Q26 = (Q11m-Q12m-2Q66m)·cs^3 - (Q22m-Q12m-2Q66m)·c^3s
Q66 = (Q11m+Q22m-2Q12m-2Q66m)·s^2c^2 + Q66m·(c^4+s^4)
其中 c=cosθ, s=sinθQ16、Q26 就是拉伸-剪切耦合项。θ=0° 或 90° 时它们为零；θ=±45° 时最大。这就是为什么 ±45° 铺层能显著提高扭转刚度，也是层合板"拉剪耦合"的数学来源。
三、经典层合板理论：A、B、D 三块矩阵
经典层合板理论（CLT）假设：层间完美粘接、直法线保持、厚度方向应力可忽略。把每层的 Q̄ 沿厚度积分，得到层合板本构：
[N] [A B] [ε⁰]
[ ] = [ ] [ ]
[M] [B D] [κ ]
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 是面内刚度，D 是弯曲刚度，B 是拉弯耦合刚度。对称铺层（如 [0/90]s）的 B 矩阵为零，层合板不产生固化后的翘曲；反对称铺层（如 [0/90]₂）B 不为零，拉伸会伴随弯曲。这是判断铺层方案是否合理的第一道检查。
以 T300/环氧单层为例（示例参数，非实测）：E₁=135 GPa，E₂=9.5 GPa，ν₁₂=0.30，G₁₂=5.0 GPa，单层厚 0.125 mm。四层 [0/90]s 板总厚 0.5 mm，中性面在几何中面。用上述公式可算出 A₁₁≈(135+9.5)/2×0.5≈36.1 kN/mm 量级，而 D₁₁ 因 z³ 项对 0° 层更敏感，明显高于同厚度各向同性板。这些量级关系决定了：面内载荷看 A，弯曲/屈曲看 D，耦合变形看 B。
四、映射到 Ansys Mechanical 的材料数据组织
在 Engineering Data 中，复合材料单层用"Orthotropic Elasticity"定义：输入 E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃。若用 Shell 单元配合 ACP 或 Mechanical 的 Layered Section，则只需平面应力下的 E₁、E₂、ν₁₂、G₁₂，其余项由求解器按平面应力条件处理。密度、热膨胀系数、强度极限（若做失效判据）分别在不同材料属性页填写。
铺层信息不在材料里，而在几何的 Section 或 ACP 的 Modeling Ply 中定义：每层的材料、厚度、角度、参考方向。Mechanical 的 Layered Section 用表格逐层录入，角度以参考坐标系为基准；ACP 则用 Rosette 定义参考方向，再在 Ply 中给角度。两者最终都生成 A、B、D 矩阵并传给求解器。
五、本节知识关系图
下图把"材料常数 → 单层刚度 → 层合板 A/B/D → Mechanical 数据组织"的映射关系集中呈现，便于对照检查输入项与物理量之间的对应。