Abaqus 二次开发教程（13）：UMAT 接口契约与编译链路——DDSDDE 不能留空，Fortran 列对齐与三处编译参数必须同源

版本声明块

• 工具/软件：Abaqus 2026（2026 FD01），Abaqus/Standard 隐式求解器
• 语言/环境：Fortran 用户子程序（固定格式源码），Intel Fortran；oneAPI 2025.0 起 ifort 已移除，须改用 ifx（第三方旁证，编译器版本结论须与本机实测对齐）。执行链参数写在abaqus_v6.env（环境配置文件），详见第 10 篇
• 本文目标：拿到一份复制即可编译运行的各向同性弹性 UMAT，并把"接口契约"（参数含义、索引约定、更新责任）与"编译链路"（三段一致、编译器环境）写成可核验的判定
• 签名说明：参数名按官方文档 UMAT 参数表（argument definitions）写法列出；同一参数存在惯用短名（NDIR/NDI、NDNHR/NSHR、NDEC/NINC、SEETEN/SSE），逐条名称与顺序以官方文档为准（官方培训讲义 R6 的 Lecture 7 给出 UMAT Interface 与等温弹性示例）

一句话结论：UMAT 必须显式写入DDSDD（材料雅可比矩阵，material Jacobian，代码里写作DDSDDE），Abaqus/Standard 不会替用户做数值求导——留空或置零时，向量长度（KMAX/KMSN）为 1 的路径会把雅可比按单位矩阵处理，静默给出错误刚度；一份正确的各向同性弹性 UMAT 在NDIR=3、NDNHR=0、NSTRESS(NTENS)=6、NSTATEV(NSTATV)=0 下，DDSDDE(1,1)必须等于 129472.13 N/mm²（E=200000、ν=0.3，mm–N–s），DDSDDE(1,2)必须等于 78129.5，DDSDDE(6,6)必须等于 76923.08，且DDSDDE(1,2)与DDSDDE(2,1)必须相等。

〇、本篇要解决的认知问题

1. UMAT 的参数表里，哪些是只读输入、哪些是必须写回的输入兼输出（STRESS、STATEV、DDSDD、SEETEN、PROPS、DVAL、CEFIELD）？
2. 材料雅可比（material Jacobian，DDSDD/DDSDDE）为什么必须显式提供？留空后 Abaqus/Standard 拿什么组装单元刚度矩阵？
3. 为什么STRESS与STATEV必须更新到增量末值（end of increment）？回滚时脏数据如何污染下一个增量？
4. NDIR/NDNHR/NSTRESS的索引约定是什么？Fortran 固定格式（fixed format）列对齐为什么会直接改变结果？
5. compile_fortran、link_sl、link_exe为什么必须同源、setvars.bat与 ifx 在哪一步生效？隐式为什么对切线刚度敏感、NTENS/KMSN向量化为什么是性能陷阱？

一、机制解析

1.1 UMAT 参数表：谁是输入、谁是必须写回的输出

按方向把参数表分开，这是 UMAT 事故的第一层排查表（逐条名称与顺序以官方文档为准）：

方向写错即违约：把"输入兼输出"当只读用，等于这一增量的结果没交出去。

Abaqus/Standard（隐式）单个增量的内部循环 ① 求解器估增量步 Δt（自动增量控制，第 8 篇） ② 每个单元、每个积分点：调用一次 UMAT 输入：起点应力、应变增量 DSTR、时间、温度、PROPS 输出：增量末 STRESS、增量末 STATEV、切线刚度 DDSDD ③ 用 DDSDD 组装单元刚度矩阵 → 迭代 → 残差检查（.msg） ④ 不收敛 → 回滚到起点、缩小 Δt、重新调用（起点重置，STATEV 不自动回滚）

这条流程揭示两件事：UMAT 不是"本构库"，而是被求解器反复调用的回调函数（callback，回调），同一个增量里可能被调用十几次，每次的DSTR/DDELTA都可能不同；求解器只管平衡方程、不管材料，把任意一点的切线刚度完全外包给 UMAT，所以DDSDD是接口契约里最贵的一项。

1.2 为什么DDSDD必须显式提供——三个后果

【铁律】第 5 条：UMAT 必须显式提供DDSDD（材料雅可比），留空或置零会静默给出错误刚度与错误结果——没有报错、没有警告，只有结果不对。机制上三层：

【经验法则】判断雅可比是否与本构一致的客观办法：在 UMAT 内对DSTR每个分量加 ε=1e-8 再跑同一套更新逻辑，把差分结果与写入的DDSDDE逐元素对比，相对误差 > 1e-6 判不一致——官方讲义 Workshop 2（Tangent Stiffness）就是这一对照（R6），第 15 篇用DQN做同类检查。

【最佳实践】显式雅可比要写成当前状态的切线，不是初始弹性矩阵：塑性段的DDSDDE应当是与应力更新算法（return mapping，返回映射）同源的一致切线（consistent tangent）；用连续体弹塑性切线替代一致切线，在强硬化或大应变增量下会拖慢收敛。

过渡句：雅可比解决"刚度对不对"，另一半契约是"数据更新到哪儿"——这一半更容易出错，因为它不报错。

1.3 增量末值与回滚：STATEV是最常见的事故点

【铁律】第 6 条：UMAT 必须把应力与状态变量更新到增量末值（end of increment），否则下一个增量读到脏数据。

增量 N 起点：STRESS(起点) ← Abaqus 传入；STATEV(起点) ← 上一增量写入的值 UMAT：STRESS、STATEV 是"输入兼输出"，必须就地写回增量末值 末点：求解器检查残差 收敛 → 末值成为增量 N+1 的起点 不收敛 → 回滚到增量 N 起点重算，但 STATEV 不会自动回滚

后果分两类。NSTATEV=0（本例）风险较小，只需保证STRESS被更新。NSTATEV>0（累积塑性应变、损伤变量、背应力）时，某次调用只更新一部分（例如只在加载分支写、卸载分支在RETURN前没写），STATEV会把中间值永久带进后续增量：损伤从 0.2 跳到 0.35、累积塑性应变出现负值，都是不可复现的漂移。【最佳实践】每次进入 UMAT 先把输入值按NSTATEV复制到本地数组，再计算、最后统一写回；卸载分支也必须显式写回（哪怕写回原值）。

另一处易漏：SEETEN（增量末弹性应变能密度）与SPD、SCD（塑性耗散能、蠕变耗散能）是可选输出，不写不会报错，但.dat与 ODB 的能量类输出会缺项（第 9 篇）。写入规则以官方文档为准。

1.4 向量与矩阵索引约定：6 个分量而不是 3×3

NDIR（正应力方向数，短名NDI）、NDNHR（剪切分量数，短名NSHR）、NSTRESS（总分量数，短名NTENS）三个数决定所有数组的形状，约定是"一维展开"而不是张量：

（前 3 行的顺序是 Abaqus 长期沿用的工程约定；逐单元类型的确切顺序以官方文档为准。）两个推论：DDSDDE是NSTRESS×NSTRESS方阵而非NDIR×NDIR；剪应力索引从NDIR+1起，写成DO I=1,NDIR+NDNHR才完整覆盖。平面应力NSTRESS=3时照搬 6 分量写法会越界读写——Fortran 里不报错，只给出错误数字。

【经验法则】按物理张量计算的本构必须手动映射回一维索引：4→12、5→23、6→31。顺序写错时主应力符号仍可能"看着对"，但剪切耦合错——这是排查成本最高的一类错误。

1.5 Fortran 固定格式：列对齐直接改变结果

固定格式（fixed format）是编译期语义而非风格：compile_fortran的选项决定固定或自由格式，用户子程序按惯例写成固定格式，列规则如下：

列号 1 2 3 4 5 6 7 72 73+ 注释 C │ 续行标志（第 6 列） │ 语句起始 └── 语句从第 7 列起，到第 72 列为止 ──┘

三种写法直接改变结果：

【铁律】第 5 条与第 6 条之所以是铁律而不是建议，是因为违反后没有运行期报错：Fortran 隐式类型与固定格式截断都会让"看起来正确的代码"安静地给出错误数字。

【最佳实践】用IMPLICIT NONE显式声明所有变量；长签名拆多行时续行标志写在第 6 列、语句从第 7 列起、每行不超过 72 列；若编译器按自由格式编译（如 ifx 的-free选项，选项名以本机与官方文档为准），必须与compile_fortran的格式选项一致，不能一半固定一半自由。

1.6 编译链路：三段必须同源

第 10 篇的执行链三段，在用户子程序这里更严格：

umat.f（固定格式源码） ① compile_fortran → 目标文件（.obj） ← 须已加载 setvars.bat ② link_sl → Abaqus/Standard 执行链；link_exe → Abaqus/Explicit 执行链 两者必须指向同一工具链，混用 ifort 目标文件 + oneAPI 库 → 链接期缺符号 ③ 求解器加载编译产物 → 每次增量调用 UMAT

【铁律】第 7 条：编译用户子程序前必须先加载编译器环境（oneAPIsetvars.bat），且compile_fortran、link_sl、link_exe三处一致。三处一致的含意是"同一个安装、同一套运行时库"：目标文件由 ifx 生成、链接却走 ifort 的库路径时，缺符号出现在链接期而非编译期——编译阶段只看单个源文件，链接阶段才需要运行时库与 Abaqus 的接口对象。

版本结论（标注生效版本）：Intel oneAPI 2025.0 起 ifort 已从工具链移除，须改用 ifx（第三方旁证，Intel 社区与发布说明口径）。此结论只影响compile_fortran里的可执行文件名，具体文件名与参数以官方文档和本机abaqus_v6.env为准；若机器上仍残留 ifort 路径，link_sl、link_exe、compile_fortran必须同时改到 ifx。

1.7 隐式为什么对切线刚度如此敏感

隐式每增量的成本由迭代次数决定，迭代次数由雅可比质量决定。纯弹性材料一次即收敛，因为DDSDDE就是本构本身；非线性材料的差别全在切线矩阵上。

调用次数在.sta与.msg里可读：.sta给出每增量的迭代次数与增量大小变化，.msg给出残差明细与 Force/Residual 对比（第 10 篇的判定口径）。【经验法则】跑 100 增量时平均迭代次数超过 5，先怀疑雅可比而不是网格或容差。

向量化（vectorization）是同一件事的另一面：KMAX/KMSN是每次调用可处理的积分点向量长度（名称以官方文档为准），写成逐点循环、向量长度为 1 就拿不到批量调用的收益；数组按(NTENS, KMAX)形态传入时，模块级或SAVE变量必须按调用实例隔离，否则线程并行会互相踩写。【最佳实践】先对齐单核结果再开并行。

二、完整代码与逐行剖析

2.1umat_iso_elastic.for（Fortran 固定格式，复制即可编译）

C ============================================================ C UMAT：各向同性线弹性（isotropic linear elasticity） C 单位制：mm - N - s【铁律】第 1 条：单位整套显式一致 C 固定格式：注释第 1 列；续行标志第 6 列；语句第 7 至 72 列 C ============================================================ SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT, 2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, 3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, 4 ORT,CEFIELD,AVDTIME,DVAL,NOEL,NPT,LAYER,KS,KN,KK,KL, 5 AMPS,DIAG,NSTEP,NDEC,DSPECIF,IORT) IMPLICIT NONE C --- 输入兼输出（必须就地写回增量末值）【铁律】第 6 条 --- DOUBLE PRECISION STRESS(*), STATEV(*), DDSDDE(NTENS,NTENS) DOUBLE PRECISION STRAN(*), DSTRAN(*), SSE, SPD, SCD, RPL DOUBLE PRECISION DDSDDT, DRPLDE, DRPLDT C --- 输入（只读）--- DOUBLE PRECISION TIME(2), DTIME, TEMP, DTEMP DOUBLE PRECISION PROPS(*), DVAL(*), COORDS(*), PREDEF(*), DPRED(*) DOUBLE PRECISION CEFIELD(*), DIAG(*), DDVAL(3,3), ANISOM(6), ORT(3,3) DOUBLE PRECISION DROT(3,3), PNEWDT, AVDTIME, DSPECIF DOUBLE PRECISION ENERGY INTEGER NDI, NSHR, NTENS, NSTATV, NPROPS, NSTEP, NDEC, IORT INTEGER NOEL, NPT, LAYER, KS, KN, KK, KL DOUBLE PRECISION E, NU, G, Kbulk, TR, DEVI(6) INTEGER I, J CHARACTER*80 CMNAME C (1) 材料参数：PROPS(1)=E、PROPS(2)=NU（由 *USER MATERIAL 传入） E = PROPS(1) NU = PROPS(2) IF (E .LE. 0.0D0 .OR. NU .LT. 0.0D0 .OR. NU .GT. 0.5D0) THEN WRITE(*,'(A,2D20.6)') 'UMAT: bad E or NU in PROPS', E, NU PNEWDT = 0.5D0 RETURN END IF G = E / (2.0D0 * (1.0D0 + NU)) ! 剪切模量：DDSDDE 的对角后三项用 G Kbulk = E / (3.0D0 * (1.0D0 - 2.0D0 * NU)) ! 体积模量：仅用于体积项，非对角不共用 G C (2) 雅可比必须显式写入：Abaqus/Standard 不会为 UMAT 数值求导【铁律】第 5 条 DO I = 1, NTENS DO J = 1, NTENS DDSDDE(I,J) = 0.0D0 END DO END DO DO I = 1, NDI DDSDDE(I,I) = Kbulk + 2.0D0 * G ! 正应力方向：λ + 2G END DO DO I = NDI + 1, NDI + NSHR DDSDDE(I,I) = G ! 剪切分量：从 NDI+1 起算，索引不可错位 END DO DO I = 1, NDI DO J = 1, NDI IF (I .NE. J) DDSDDE(I,J) = Kbulk - (2.0D0/3.0D0) * G ! λ = K - 2G/3 END DO END DO C DDSDDE 是 NSTRESS×NSTRESS（即 NTENS×NTENS）方阵，不是 NDIR×NDIR；对称性检查见练习 1 C (3) 应变增量与应力：更新到增量末值（end of increment）【铁律】第 6 条 TR = 0.0D0 DO I = 1, NTENS TR = TR + STRAN(I) END DO DO I = 1, NTENS DEVI(I) = DSTRAN(I) IF (I .LE. NDI) DEVI(I) = DEVI(I) - TR / 3.0D0 ! 偏应变：减掉体积部分 STRESS(I) = STRESS(I) + 2.0D0 * G * DEVI(I) + Kbulk * TR END DO C (4) 状态变量：本模型 NSTATV=0，无状态变量；若有，必须显式写回增量末值 IF (NSTATV .GT. 0) THEN DO I = 1, NSTATV STATEV(I) = STATEV(I) END DO END IF C (5) 可选能量输出：不写不会报错，但 .dat 与 ODB 的能量类输出会缺项（以官方文档为准） C 弹性应变能密度按增量累加：G*(偏应变增量点积) + 0.5*K*(体应变)^2 ENERGY = 0.0D0 DO I = 1, NTENS ENERGY = ENERGY + DEVI(I) * DSTRAN(I) END DO SSE = SSE + G * ENERGY + 0.5D0 * Kbulk * TR * TR C SPD、SCD、RPL、DDSDDT、DRPLDE、DRPLDT 保持 0（本模型无耗散、无温度依赖） RETURN END

注意：签名尾部的DROT、PNEWDT、DSPECIF、IORT属于不同版本的接口扩展，逐条名称与顺序以官方文档为准；本例只使用PROPS、DSTRAN、NTENS、NDI、NSHR、NSTATV、NSTEP、NDEC、NOEL、NPT、PNEWDT与DDSDDE。【最佳实践】正式使用前按官方 UMAT Interface 页逐条核对签名，再按本机compile_fortran的格式选项调整列宽。

逐行剖析关键段：

• 第 (1) 段：PROPS由输入文件的*USER MATERIAL, CONSTANTS=2提供，长度是NPROPS；DVAL是求解器算好的派生量，两者只读。把参数写在PROPS里而不是硬编码，第 12 篇的参数表才能直接驱动 UMAT。这里的防御性检查给出PNEWDT=0.5D0（请求求解器把增量步长建议减半），是用户唯一能影响增量控制的字段。
• 第 (2) 段：先清零、再写对角、再写非对角，顺序不能颠倒——非对角分支只覆盖I, J ≤ NDI，先写会覆盖对角线。λ（Lamé 常数）用Kbulk - 2G/3表达比E*NU/((1+NU)(1-2NU))更稳：NU → 0.5时Kbulk发散，这个写法把不可压缩极限的风险摆出来，而不是藏在一串代数里。
• 第 (3) 段：STRESS是"输入兼输出"，用STRESS(I) = STRESS(I) + 2G·偏应变 + K·体应变就地累加，等于把起点值带到增量末。若写成STRESS(I) = 新值，必须确保STRESS仍是起点值——这正是"没更新到增量末值"事故的常见写法差异。
• 第 (4) 段：NSTATV=0时这段看似冗余，但保留显式写回是【最佳实践】：后续若增加状态变量（累积塑性应变、损伤），这段就是必须存在的骨架。
• 第 (5) 段：能量项是可选输出，把SEETEN/SSE写成累加而不是赋值是官方参数表对可选输出的常见约定；是否必须累加、是否要求用更新前的应变，以官方文档为准。

【铁律】第 7 条：编译前必须加载编译器环境。setvars.bat设置PATH与运行时库路径，compile_fortran只在已加载的环境里找编译器；未加载时编译阶段报"找不到编译器"、链接阶段报缺符号。

2.2abaqus_v6.env与提交命令

工作目录 C:\work\umat_case：umat_case.inp（含 *USER MATERIAL, CONSTANTS=2、model=USER）、 umat_iso_elastic.for（上一节源码）、abaqus_v6.env（与 .inp 同目录，作用于本目录作业） # abaqus_v6.env（三处必须同源：同一个 oneAPI 安装、同一个编译器） compile_fortran = "ifx /assume bytsint -f77 -MD -c $source_file -obj $object_file" link_sl = "ifx /assume bytsint -f77 -MD -dll $object_file $dll_file" link_exe = "ifx /assume bytsint -f77 -MD -exe $object_file $exe_file" # 环境文件只写工具链；编译器运行时路径由 setvars.bat 提供

提交前须在同一会话先加载编译器环境：

call "C:\Program Files (x86)\Intel\oneAPI\setvars.bat" abaqus job=umat_case user=umat_iso_elastic.for cpus=1

判定链路是否真走通（三处证据，按第 10 篇口径）：

2.3DDSDDE的预期值（对照表，用于自检）

E = 200000.0、NU = 0.3，单位 mm–N–s：G = 76923.077、Kbulk = 115384.615、λ =K − 2G/3 = 25641.026。

对称性（DDSDDE(I,J) = DDSDDE(J,I)）必查：不对称时不报错，但残差方向错，隐式收敛变慢甚至不收敛。

三、常见报错与排查

1. 现象：作业"跑完"但位移比理论值大 3~5 倍，且没有任何报错。
根因：DDSDD留空或置零，向量长度为 1 的路径把雅可比当单位矩阵处理（1.2 节），刚度被静默替换。解法：按第 (2) 段显式写入，并用 2.3 节对照表逐元素核对；与数值差分结果对比，相对误差 > 1e-6 判失败。

2. 现象：编译通过，链接期报缺符号（undefined symbol），或运行期报找不到UMAT。
根因：compile_fortran用 ifx 生成目标文件，link_sl/link_exe仍指向 ifort 路径（oneAPI 2025.0 起 ifort 已移除），三处不同源；或提交时未加载setvars.bat，运行时库路径缺失。解法：三处指向同一安装与编译器，并在加载setvars.bat后的同一会话提交；用 2.2 节三处证据表定位断链段。报错原文以官方文档与本机日志为准。

3. 现象：结果随网格变细单调偏离，主应力方向看似正确、剪切耦合不对。
根因：NSTRESS的一维索引与物理 3×3 张量混用（1.4 节），或平面应力下按 6 分量读写越界（Fortran 不报错）。解法：在子程序开头打印NDI、NSHR、NTENS、NSTATV的实际值，按打印值确定循环上界，并把映射关系写成注释表。

4. 现象：状态变量出现负值或跳变（如累积塑性应变 0.2 → 0.35 后回落）。
根因：卸载分支在RETURN前没把STATEV更新到增量末值，或收敛失败回滚后STATEV未被回滚（1.3 节）。解法：每次调用先复制到本地数组、末尾统一写回；打印NSTEP、NDEC、NOEL、NPT，用"同一增量被重复调用"识别回滚。【经验法则】每增量调用次数超过 5，先查状态变量写回，再查网格。

5. 现象：单核跑完 60 秒，改成cpus=8后反而变慢，或报子程序不支持并行。
根因：按线程并行调用时，模块级或SAVE共享变量被不同实例踩写；向量化要求数组按(NTENS, KMAX)形态传入，逐点循环写法拿不到收益。解法：先对齐单核结果（第 10 篇cpus判定），再把共享状态改成调用内局部数组；官方讲义的 threads 与 MPI 章节（R6 的 Appendix 5、Workshop 5/6）列出并行调用的写法约束，具体限制以官方文档为准。

四、动手练习

练习 1：用解析解验证DDSDD与应力更新（容差）。
在E=200000、NU=0.3下取单轴应变增量（DSTRAN(1)=1e-3、其余为 0）。判定标准：DDSDDE(1,1)=129472.13、DDSDDE(1,2)=78129.5、DDSDDE(4,4)=76923.08（相对误差 < 1e-6）；由DDSDDE算出的应力增量与直接按本构算的增量末应力逐元素相对误差 < 1e-9；DDSDDE(I,J) − DDSDDE(J,I)全为 0。任一不满足则第 (2) 段有索引或公式错误。

练习 2：把状态变量加进契约（可观察输出判定）。
把NSTATEV设为 1，在 UMAT 中把DSTRAN(1)累加到STATEV(1)，并在每次调用时打印NSTEP、NDEC、NOEL、NPT。判定标准：单轴 10 增量后STATEV(1)与解析累加值 0.01 的相对误差 < 1e-6；打印的调用次数 ≥ 增量数，出现同一增量多次调用即为回滚，此时STATEV必须仍能从起点重建。

练习 3：把编译链路写成可数判定（三处证据）。
在一个空目录放置abaqus_v6.env（2.2 节版本）与.for源码，提交含*USER MATERIAL, CONSTANTS=2的作业。判定标准：目录出现.obj与链接产物；.log或.msg出现一次 UMAT 打印的行；.sta中每增量迭代次数为 1。未加载setvars.bat重跑时至少一段失败——这条对照就是【铁律】第 7 条的可复现验证。

五、小结与下一篇预告

本篇把 UMAT 的接口契约拆成四条硬约束：DDSDD必须显式提供、STRESS与STATEV必须更新到增量末值、NDIR/NDNHR/NSTRESS索引约定不可混用、固定格式列对齐与三处编译参数必须同源且须先加载setvars.bat。第 14 篇（VUMAT 与显式用户材料）对比显式与隐式的增量格式：VUMAT 按块（block）处理多个单元，不需要切线刚度，但需要密度与步长的显式契约与可调用单元类型限制。雅可比在显式里不再是收敛工具，而是性能与内存布局问题——四条契约中只有"单位整套显式一致"和"更新到增量末值"仍然成立。第 15 篇用DQN做精度与DDSDD一致性检查。

本篇认知问题回显（FAQ）

Q1：UMAT 的参数表里，哪些是只读输入、哪些是必须写回的输入兼输出（STRESS、STATEV、DDSDD、SEETEN、PROPS、DVAL、CEFIELD）？

A：NSTEP、NDEC、DDELTA、TIME、CEFIELD、PROPS、DVAL是只读输入，STRESS、STATEV、DDSDD必须写回增量末值，SEETEN是可选输出。

Q2：材料雅可比（material Jacobian，DDSDD/DDSDDE）为什么必须显式提供？留空后 Abaqus/Standard 拿什么组装单元刚度矩阵？

A：Abaqus/Standard 不为 UMAT 数值求导，留空时向量长度为 1 的路径把雅可比按单位矩阵处理；E=200000、ν=0.3下DDSDDE(1,1)必须是 129472.13 N/mm²。

Q3：为什么STRESS与STATEV必须更新到增量末值（end of increment）？回滚时脏数据如何污染下一个增量？

A：本次写入的值就是下一增量的起点，回滚只恢复起点数据而STATEV不自动回滚，未写回的分支（典型是卸载分支）会把中间值永久带进后续增量。

Q4：NDIR/NDNHR/NSTRESS的索引约定是什么？Fortran 固定格式（fixed format）列对齐为什么会直接改变结果？

A：3D 实体单元按NDIR=3、NDNHR=0、NSTRESS=6一维展开（4=12、5=23、6=31），DDSDDE是 6×6 方阵；固定格式续行标志必须在第 6 列、语句从第 7 列到第 72 列。

Q5：compile_fortran、link_sl、link_exe为什么必须同源、setvars.bat与 ifx 在哪一步生效？隐式为什么对切线刚度敏感、NTENS/KMSN向量化为什么是性能陷阱？

A：compile_fortran只产目标文件、运行时库在链接期才用到，三处须同源且先加载setvars.bat（oneAPI 2025.0 起 ifort 移除须用 ifx）；隐式迭代方向由DDSDD决定，KMSN/KMAX为 1 时失去向量化收益。