3. 弹性本构模型选型

弹性本构模型是材料力学行为描述中最基础、也是工程中使用频率最高的一类模型。它的核心假设是：应力与应变之间存在一一对应的可逆关系，卸载后变形完全恢复。但“弹性”本身并不是一个单一模型，而是一族模型。选型的关键不在于哪个模型更“高级”，而在于当前工程对象的材料对称性、载荷路径、分析目标是否落在该模型的适用边界之内。

一、三类弹性假设的数学结构与适用边界

最一般形式的线弹性本构关系由广义胡克定律给出：σ_ij = C_ijkl ε_kl，其中四阶刚度张量 C_ijkl 最多有 81 个分量。利用应力张量和应变张量的对称性，独立分量降至 36 个；再利用应变能密度与加载路径无关这一热力学条件，C_ijkl = C_klij，独立常数进一步降为 21 个。这 21 个常数描述的是最一般的各向异性线弹性材料，例如某些单晶、木材顺纹与横纹差异极大的区域、以及增材制造中具有明显方向性的沉积态材料。

当材料内部存在一个弹性对称面时，独立常数从 21 个降到 13 个（单斜）；存在三个互相正交的弹性对称面时，即得到正交各向异性模型，独立常数降为 9 个：E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃。正交各向异性是工程中最常见的各向异性形式，层合复合材料单层、拉拔金属棒材、轧制钢板、木材、蜂窝芯材都近似属于此类。其柔度矩阵形式为：

[ε1]   [1/E1    -ν21/E2   -ν31/E3      0        0        0    ][σ1]
[ε2]   [-ν12/E1   1/E2    -ν32/E3      0        0        0    ][σ2]
[ε3] = [-ν13/E1  -ν23/E2    1/E3       0        0        0    ][σ3]
[γ23]  [  0        0         0      1/G23      0        0    ][τ23]
[γ13]  [  0        0         0        0      1/G13      0    ][τ13]
[γ12]  [  0        0         0        0        0      1/G12  ][τ12]

注意柔度矩阵必须对称，因此存在互等关系：ν₁₂/E₁ = ν₂₁/E₂，ν₁₃/E₁ = ν₃₁/E₃，ν₂₃/E₂ = ν₃₂/E₃。这意味着 9 个常数中只有 3 个主泊松比是独立的，另外 3 个由互等关系导出。工程中若直接输入 6 个互不相关的泊松比，会导致刚度矩阵非正定，求解器报错或结果失真。

当材料在三个方向性能完全一致时，正交各向异性退化为各向同性模型，独立常数只剩 2 个：弹性模量 E 和泊松比 ν，剪切模量由 G = E/[2(1+ν)] 导出。各向同性假设适用于大多数金属在无显著织构、无冷作方向性、晶粒随机取向的情况，也适用于均质各向同性聚合物和陶瓷。

二、选型判断：从材料状态到模型层级

选型的第一步不是查手册，而是判断材料在当前分析尺度下是否真的表现出方向性。一个实用的判断流程是：先看材料是否存在可辨识的制造方向或微观结构方向。轧制钢板在轧制方向与横向的屈服强度和弹性模量差异通常小于 5%，若分析目标只是整体刚度或低应力弹性变形，各向同性模型往往足够；但若进行冲压回弹预测，塑性各向异性才是主导，弹性阶段仍可用各向同性。相反，单向碳纤维增强环氧树脂层合板单层，沿纤维方向模量约 135 GPa，垂直纤维方向约 9 GPa，相差一个数量级以上，此时必须使用正交各向异性模型，否则层间应力与铺层顺序效应完全无法体现。

第二步是判断载荷路径是否激活了方向差异。即使材料本身正交各向异性，若结构只承受单一方向拉伸且横向应力可忽略，一维线弹性模型 σ = Eε 也可能满足精度。例如拉索、杆系结构中只关心轴向刚度时，只需该方向的等效弹性模量。但一旦出现多轴应力状态、弯扭组合或热膨胀约束，方向耦合项就会起作用，必须回到完整正交各向异性描述。

第三步是检查参数是否可获取。正交各向异性需要 9 个常数，若试验条件只允许做单轴拉伸，则无法完整标定。此时工程上常采用两种折中：一是用横观各向同性模型（如单向复合材料，独立常数降为 5 个：E_L、E_T、ν_LT、ν_TT、G_LT），二是用各向同性等效模量做初步估算，但必须明确这会在多轴载荷下引入误差。

三、工程参数获取的实际操作

以单向玻璃纤维/环氧层合板单层为例，说明正交各向异性参数的获取过程。假设纤维方向为 1 轴，面内垂直纤维方向为 2 轴，厚度方向为 3 轴。典型教学示例参数如下（非实测值，仅用于说明流程）：

输入有限元软件时，需注意泊松比互等关系校验。例如 ν₂₁ = ν₁₂ × E₂/E₁ = 0.28 × 12/45 ≈ 0.0747，ν₃₁ = ν₁₃ × E₃/E₁ ≈ 0.0747，ν₃₂ = ν₂₃ × E₃/E₂ = 0.40。若软件要求输入全部 6 个泊松比，必须按此关系填写，否则刚度矩阵可能非正定。一个快速检查方法是：计算柔度矩阵的行列式，若为负或接近零，说明参数组合不物理。

四、常见误用与边界条件

最常见的误用是把正交各向异性参数直接用于各向同性求解器，或反过来把各向同性参数填入正交各向异性模型。前者会丢失方向耦合，后者会人为引入不存在的方向差异。另一个高频问题是忽略材料主方向与全局坐标系的夹角。正交各向异性模型的 1、2、3 轴必须与材料主方向对齐，若纤维方向与全局 X 轴成 30°，则需要在材料定义中指定局部坐标系，或通过坐标变换将刚度矩阵旋转到全局系。旋转后的刚度矩阵会出现拉剪耦合项，这是物理上真实存在的效应，不能忽略。

适用边界方面，线弹性模型本身要求应力水平低于屈服强度或基体开裂起始点。对于金属，弹性阶段通常只到 0.2% 残余应变对应的应力；对于复合材料，基体微裂纹可能在远低于极限强度的水平就开始累积，此时弹性模型仍可用于刚度预测，但强度预测必须切换到损伤或塑性模型。此外，线弹性假设要求小变形，当应变超过约 2%～5% 时，几何非线性效应显著，需改用大变形框架下的弹性模型或超弹性模型。

最后需要强调：弹性本构模型选型不是一次性决策，而是随分析目标变化的。同一块层合板，在模态分析中可用正交各向异性弹性；在低速冲击中，弹性模型只能用于初始响应阶段，后续必须引入损伤；在热应力分析中，还需补充热膨胀系数 α₁、α₂、α₃，且这些系数同样具有方向性。选型的本质，是让模型的对称性、参数数量与当前工程问题的物理复杂度和可用数据相匹配。