1. 各向异性材料本构建模工程导论
在工程材料中，各向同性只是一种理想化假设。绝大多数实际材料——纤维增强复合材料、轧制金属板材、增材制造构件、木材、骨骼、单晶高温合金——其刚度、强度、热膨胀乃至导电导热性能都随方向变化。这种方向依赖性就是各向异性。它并非缺陷，而是材料内部微结构（纤维取向、晶粒织构、层片排列、孔隙形态）在宏观力学行为上的直接映射。理解并量化这种映射，正是本构模型存在的理由。
一、各向异性现象：从微结构到方向相关响应
以最常见的单向碳纤维/环氧层合板为例。沿纤维方向（记为1轴）拉伸时，纤维承担主要载荷，弹性模量可达 130～180 GPa；而垂直纤维方向（2轴）拉伸时，载荷主要由基体承担，模量往往只有 8～12 GPa。两者相差一个数量级以上。更关键的是，当沿与纤维成 45° 方向加载时，试件不仅产生轴向伸长，还会出现明显的拉剪耦合——正应力诱发出剪应变。这是各向同性材料绝不会出现的现象，也是各向异性本构必须用完整张量描述的根本原因。
工程中常见的各向异性层次可归纳为：
- 横观各向同性：存在一个对称轴，垂直于该轴的平面内性能相同，如单向带、单向纤维复合材料。
- 正交各向异性：存在三个相互垂直的弹性对称面，如正交铺层层合板、轧制板材。
- 一般各向异性：无任何对称面，如某些单晶、编织复合材料，独立弹性常数最多可达21个。
方向相关性的工程后果非常直接：用各向同性模型去预测一个正交各向异性层合板的刚度，误差可能超过 50%；在屈曲、振动频率、疲劳寿命预测中，这种误差会被进一步放大。因此，各向异性本构建模不是学术上的精细修饰，而是结构安全与轻量化设计的前提。
二、本构模型的定义与作用
本构模型（constitutive model）是描述材料应力—应变—温度—时间之间关系的数学表达。它回答一个核心问题：给定当前的变形状态（应变、应变率、温度、损伤变量等），材料会输出怎样的应力？在有限元求解中，本构模型是每个积分点上的“材料心脏”，它把几何变形转化为节点内力，驱动整个平衡迭代。
一个完整的各向异性本构模型通常包含三个层次：
- 弹性（或粘弹性）关系：用刚度矩阵
C或柔度矩阵S建立应力与应变的线性映射，如σ = C : ε。正交各向异性材料有 9 个独立工程常数：E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃。 - 强度准则与损伤演化：判断材料何时失效、以何种模式失效（纤维拉伸、纤维压缩、基体开裂、分层），并据此退化刚度。常用准则包括 Tsai-Wu、Hashin、Puck 等。
- 非线性与路径相关行为：塑性、蠕变、粘塑性、损伤累积等，需要引入内变量和演化方程。
本构模型的作用可以概括为“三座桥”：连接微观机制与宏观响应、连接材料试验与结构仿真、连接设计需求与制造工艺。没有可靠的本构模型，复合材料结构的虚拟验证就失去了物理基础。
三、课程知识地图与建模工作流
本课程围绕“从现象到模型、从模型到仿真、从仿真到验证”的主线展开，知识地图可划分为四个模块：基础理论（张量、坐标变换、弹性对称性）、本构构建（刚度矩阵、强度准则、损伤演化、非线性）、参数获取（试验设计、参数识别、反演标定）、工程应用（有限元实现、层合板分析、多尺度耦合、验证与确认）。
建模工作流则是一条可复现的工程链路：明确材料体系与载荷工况 → 确定各向异性类型与对称性 → 选择本构形式 → 设计试验并获取参数 → 编写或配置材料子程序 → 在有限元中验证单元级响应 → 与试验数据对比并修正 → 应用于结构级分析。每一步都有明确的输入、输出和判断依据。
四、从本构到仿真：一个可执行的起点
为了让上述工作流不停留在概念层面，这里给出一个最小可执行示例。假设我们要为一层单向碳纤维/环氧建立三维正交各向异性弹性模型，已知工程常数为：E₁ = 135 GPa，E₂ = E₃ = 10 GPa，ν₁₂ = ν₁₃ = 0.30，ν₂₃ = 0.45，G₁₂ = G₁₃ = 5.0 GPa，G₂₃ = 3.5 GPa。在有限元材料卡片中，这些常数需要转换为刚度矩阵分量。以 1-2 平面为例，柔度矩阵 S 的独立分量为：
S11 = 1/E1 = 7.407e-12 Pa⁻¹
S22 = 1/E2 = 1.000e-10 Pa⁻¹
S12 = -ν12/E1 = -2.222e-12 Pa⁻¹
S66 = 1/G12 = 2.000e-10 Pa⁻¹通过 C = S⁻¹ 求逆，可得 C11 ≈ 136.2 GPa，C22 ≈ 10.1 GPa，C12 ≈ 3.03 GPa，C66 = 5.0 GPa。将这些值写入材料子程序后，对单单元施加 1% 的 1 方向应变，应得到 σ₁ ≈ 1.362 GPa；施加 1% 的 2 方向应变，应得到 σ₂ ≈ 0.101 GPa。如果仿真输出与手算一致，说明弹性部分实现正确；若不一致，应优先检查坐标变换矩阵和刚度矩阵的对称性。