2. 材料本构模型的基本数学框架

材料本构模型的任务，是在给定变形历史下，把应变（或变形率）与应力（或应力率）联系起来。要写出一个可用的本构关系，必须先回答四个基础问题：用什么量描述变形与受力？这些量在刚体转动下是否保持不变？关系式写成全量还是增量？不可逆变形如何通过内变量和耗散来描述？本节围绕这四个问题建立统一框架。

一、应力与应变的度量

描述变形最直接的是变形梯度 F，它把参考构形中的微元线段映射到当前构形。由 F 可构造右 Cauchy-Green 张量 C = F^TF 和 Green-Lagrange 应变 E = (C − I)/2。这两个量都定义在参考构形上，且对刚体转动不敏感，因此适合作为超弹性等路径无关模型的自变量。小变形情形下，E 退化为工程应变张量 ε = (∇u + ∇u^T)/2。

应力度量必须与应变度量配对，才能保证功率共轭。与 E 配对的是第二类 Piola-Kirchhoff 应力 S；与变形率张量 D 配对的是 Cauchy 应力 σ。二者通过 S = J F⁻¹σF^−T 联系，J = det F。工程中常用的名义应力（第一类 Piola-Kirchhoff 应力 P = JσF^−T）与 F 共轭。选型时若混用不共轭的应力应变对，能量计算和切线刚度都会出错。

二、客观性原理

客观性（材料标架无差异性）要求：在叠加刚体运动后，本构关系形式不变。设当前构形叠加转动 Q，则 Cauchy 应力按 σ* = QσQ^T 变换。若本构写成 σ = f(F)，则必须满足 f(QF) = Qf(F)Q^T。直接以 σ 的时间导数 σ̇ 构造率型关系会破坏客观性，因为 σ̇ 中混入了刚体转动贡献。解决办法是采用客观应力率，如 Jaumann 率 σ^▽ = σ̇ − Wσ + σW，其中 W 是自旋张量。这样 σ^▽ 在刚体转动下与 σ 同步旋转，可用于金属塑性等率型模型。

三、全量格式与增量格式

全量格式直接给出 σ = g(ε)，应力只取决于当前应变，与路径无关。线弹性 σ = C:ε 和超弹性 S = ∂W/∂E 都属于此类。它的优点是简单、能量守恒、易于求导得到切线模量；缺点是无法描述加载历史相关的塑性、损伤和蠕变。

增量格式把应力增量与应变增量联系起来：dσ = C^ep:dε，其中 C^ep 是弹塑性切线模量。它沿加载路径积分，能反映历史效应。有限元实现中，增量格式通常写成返回映射算法：先按弹性预测 σ^tr = σ_n + C:Δε，再判断是否屈服，若屈服则把应力拉回屈服面并更新内变量。全量与增量并非互斥：超弹性可写成全量，但其切线刚度在数值上仍以增量形式使用。

四、内变量与耗散表述

不可逆过程需要内变量记录历史。常见内变量包括等效塑性应变 p、背应力 α（随动硬化）、损伤变量 d、黏塑性过应力等。以等向硬化 von Mises 塑性为例，屈服函数 f = q − (σ_y0 + H p) ≤ 0，其中 q 为等效应力，H 为硬化模量。塑性流动时 p 按 ṗ = λ̇ 累积，λ̇ 为塑性乘子。

热力学第二定律要求耗散非负：D = σ:ε̇^p − A·α̇ ≥ 0，其中 A 是与内变量共轭的热力学力。对关联塑性，流动方向由屈服函数梯度给出，ε̇^p = λ̇ ∂f/∂σ，内变量演化也由 λ̇ 驱动。这样构造的模型自动满足最大耗散原理，数值上更稳定。

五、选型时的框架检查清单

面对一个具体材料，可按以下顺序核对框架是否自洽：第一，变形是否大转动？若是，应选参考构形量或客观率，避免用 σ̇ 直接积分。第二，路径是否相关？若相关，必须用增量格式并定义内变量演化方程。第三，应力应变对是否功率共轭？混用会导致能量不守恒。第四，耗散是否非负？关联塑性可自动满足，非关联或损伤模型需单独验证。第五，切线模量是否与返回映射一致？数值实现中常用一致切线模量以保证 Newton 迭代收敛。

以某低碳钢单轴拉伸为例（示例参数，非实测）：E = 210 GPa，ν = 0.3，σ_y0 = 235 MPa，H = 1.2 GPa。采用增量格式，给定应变增量 Δε = 0.002。弹性预测应力增量为 420 MPa，预测后等效应力超过屈服面，进入塑性修正。返回映射后实际应力增量约为 235 + 1.2×10³×p 对应的值，塑性应变增量由一致性条件 λ̇ = (∂f/∂σ : C : Δε) / (∂f/∂σ : C : ∂f/∂σ + H) 确定。通过比较预测与修正后的应力，可判断该步是否屈服以及硬化量是否合理。