4. 状态变量管理

在UMAT开发中，状态变量（Solution-Dependent State Variables，简称SDV）是连接增量步之间材料历史信息的唯一通道。ABAQUS在每次调用UMAT时，只把当前积分点的总应变、应变增量、时间增量以及上一增量步结束时的状态变量数组传入，UMAT计算完成后把更新后的状态变量写回。因此，任何需要跨增量步记忆的量——等效塑性应变、背应力张量、损伤因子、相变体积分数、粘塑性内变量——都必须通过SDV保存。状态变量管理不当，轻则结果漂移，重则收敛失败或出现“一步一个样”的不可复现结果。

一、SDV的声明与命名规范

SDV的数量由输入文件中的 *DEPVAR 关键字声明，例如 *DEPVAR, N=12 表示每个积分点保留12个状态变量。UMAT子程序中的 STATEV(NSTATV) 数组长度必须与之一致，否则会出现数组越界或数据错位。声明数量时建议留出余量，因为后期增加内变量时若修改 N，旧的重启动文件将无法直接续算。

命名规范的核心目标是“看名字就知道物理含义、张量阶次和存储顺序”。推荐采用“物理量缩写_分量”的形式，并统一大小写。例如：

• EQPLAS：等效塑性应变，标量，1个SDV；
• BACK11, BACK22, BACK33, BACK12, BACK13, BACK23：背应力张量的6个独立分量，按Voigt顺序存储；
• DAMAGE：损伤因子，标量；
• FVPREV：上一增量步的屈服函数值，用于判断是否发生卸载。

对于张量型内变量，必须固定存储顺序并在注释中写明。例如背应力采用 [11,22,33,12,13,23] 的Voigt顺序，剪切分量是否乘以2也要明确。若同一材料模型包含多个机制（如塑性+损伤+蠕变），建议按机制分块连续存放，并在代码开头用注释列出SDV索引映射表，避免后续维护时“第7个到底是背应力还是损伤”这类问题。

二、初始化：第一次调用与重启动

SDV的初始化发生在两个时机。第一是分析开始时，ABAQUS默认将 STATEV 全部置零，但并非所有内变量的物理初值都是零。例如初始损伤因子可能为0，初始背应力为0，但若模拟预应变材料，等效塑性应变初值可能不为零。此时应通过 *INITIAL CONDITIONS, TYPE=SOLUTION 在输入文件中指定初值，或在UMAT中利用 KINC（增量步号）和 KSTEP（分析步号）判断是否为第一步，在 KINC==1 时对特定SDV赋初值。

第二是重启动分析。重启动时ABAQUS会从重启动文件读入SDV，UMAT不应再次覆盖。因此初始化逻辑必须写成“仅当 KINC==1 且 KSTEP==1 时执行”，否则重启动后历史变量会被错误清零。一个稳妥的做法是在UMAT开头设置一个逻辑标志：

IF (KSTEP.EQ.1 .AND. KINC.EQ.1) THEN
    STATEV(1) = 0.0D0      ! EQPLAS
    STATEV(2:7) = 0.0D0    ! BACK stress
    STATEV(8) = 0.0D0      ! DAMAGE
END IF

注意：ABAQUS在每次调用UMAT时传入的 STATEV 已经是上一增量步结束时的值，因此上述初始化只在第一步第一次调用时生效，后续调用不会重复执行。

三、读写：在UMAT中的正确操作顺序

UMAT中SDV的读写遵循“先读旧值、再算新值、最后写回”的顺序。具体来说，进入UMAT时 STATEV 中存放的是本增量步开始时的值（即上一增量步收敛后的值），应先将需要的历史量读入局部变量，例如 EQPLAS_OLD = STATEV(1)，然后基于应变增量进行本构积分，计算出新的 EQPLAS_NEW，最后在返回前写回 STATEV(1) = EQPLAS_NEW。若在计算过程中直接修改 STATEV，一旦迭代不收敛需要重新计算，旧值已被破坏，会导致错误累积。

对于需要迭代求解的隐式本构（如径向返回算法），建议在局部变量中完成全部迭代，仅在收敛后一次性写回SDV。这样即使ABAQUS因收敛失败而重新调用UMAT，STATEV仍保持上一增量步的干净值。

四、调试输出策略

SDV的调试输出有三种常用途径。第一种是利用 WRITE(6,*) 或 WRITE(7,*) 在UMAT中直接打印，适合快速定位“某个积分点在第几步出现异常值”。但要注意，UMAT在每个积分点、每次迭代都会被调用，无节制打印会迅速撑爆 .dat 或 .msg 文件。建议加条件过滤，例如仅当 KINC 等于特定值、或当 EQPLAS_NEW 超过阈值时才输出：

IF (EQPLAS_NEW .GT. 0.5D0 .AND. KINC .EQ. 10) THEN
    WRITE(7,*) 'KINC=',KINC,' IP=',NOEL,' EQPLAS=',EQPLAS_NEW
END IF

第二种是通过 *EL PRINT 或 *NODE PRINT 将SDV输出到 .dat 文件，适合在分析结束后批量检查。例如 *EL PRINT, ELSET=ALL, FREQUENCY=1 配合 SDV 输出，可以查看每个单元每个积分点的状态变量演化。第三种是在 *OUTPUT, FIELD 中请求 SDV，写入 .odb 后在可视化模块中云图显示，适合观察损伤或等效塑性应变的分布。

调试时建议遵循“先标量、后张量；先单点、后全场”的原则。先确认 EQPLAS 和 DAMAGE 这类标量SDV的演化是否单调合理，再检查背应力张量的分量是否满足对称性和客观性。若发现SDV在卸载步出现异常跳变，优先检查是否在迭代过程中提前写入了 STATEV，或初始化逻辑是否在重启动时被错误触发。

五、一个可复用的SDV管理检查清单

状态变量管理看似琐碎，但它直接决定了本构模型能否在多个增量步之间保持物理一致性。把命名、初始化、读写和输出这四个环节固定成规范，后续增加损伤、蠕变或相变内变量时，只需扩展索引表而不必重构整个UMAT的数据流。