Abaqus 二次开发教程13UMAT 接口契约与编译链路——DDSDDE 不能留空Fortran 列对齐与三处编译参数必须同源版本声明块工具/软件Abaqus 20262026 FD01Abaqus/Standard 隐式求解器语言/环境Fortran 用户子程序固定格式源码Intel FortranoneAPI 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代码里写作DDSDDEAbaqus/Standard 不会替用户做数值求导——留空或置零时向量长度KMAX/KMSN为 1 的路径会把雅可比按单位矩阵处理静默给出错误刚度一份正确的各向同性弹性 UMAT 在NDIR3、NDNHR0、NSTRESS(NTENS)6、NSTATEV(NSTATV)0 下DDSDDE(1,1)必须等于 129472.13 N/mm²E200000、ν0.3mm–N–sDDSDDE(1,2)必须等于 78129.5DDSDDE(6,6)必须等于 76923.08且DDSDDE(1,2)与DDSDDE(2,1)必须相等。〇、本篇要解决的认知问题UMAT 的参数表里哪些是只读输入、哪些是必须写回的输入兼输出STRESS、STATEV、DDSDD、SEETEN、PROPS、DVAL、CEFIELD材料雅可比material JacobianDDSDD/DDSDDE为什么必须显式提供留空后 Abaqus/Standard 拿什么组装单元刚度矩阵为什么STRESS与STATEV必须更新到增量末值end of increment回滚时脏数据如何污染下一个增量NDIR/NDNHR/NSTRESS的索引约定是什么Fortran 固定格式fixed format列对齐为什么会直接改变结果compile_fortran、link_sl、link_exe为什么必须同源、setvars.bat与 ifx 在哪一步生效隐式为什么对切线刚度敏感、NTENS/KMSN向量化为什么是性能陷阱一、机制解析1.1 UMAT 参数表谁是输入、谁是必须写回的输出按方向把参数表分开这是 UMAT 事故的第一层排查表逐条名称与顺序以官方文档为准参数方向含义STRESS输入兼输出进来是增量起点应力必须写回增量末值STATEV输入兼输出状态变量必须写回增量末值DDSDD代码里DDSDDE输出材料雅可比切线刚度NSTRESS×NSTRESSSEETEN短名SSE可选输出增量末弹性应变能密度SPD/SCD可选输出塑性、蠕变耗散能增量STRAN/DSTRDSTRAN输入起点应变与应变增量NSTEP/NDEC输入增量编号与增量计数定义以官方文档为准DDELTA输入与DSTR同长度的分量增量TIME输入TIME(1)总时间、TIME(2)增量起点时间步长DTIMENDIR/NDNHR/NSTRESS输入正应力方向数、剪切分量数、总分量数PROPS/NPROP输入材料参数与个数来自*USER MATERIAL, CONSTANTSnDVAL输入求解器算好的派生量各向同性时对应 E、G、KCEFIELD输入场变量数组USDFLD写入的值在此读回方向写错即违约把输入兼输出当只读用等于这一增量的结果没交出去。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材料雅可比留空或置零会静默给出错误刚度与错误结果——没有报错、没有警告只有结果不对。机制上三层层面发生了什么可观察后果组装刚度矩阵Abaqus/Standard 用DDSDDE生成单元刚度矩阵它不会为 UMAT 做数值差分求导置零后矩阵奇异或退化迭代方向完全错误向量长度为 1 的路径KMAX/KMSN每次调用处理的积分点数为 1 时部分版本把雅可比按单位矩阵处理而非报错位移场被放大数倍网格越粗误差越大收敛判据失效隐式迭代靠雅可比预测残差下降方向预测矩阵与实际本构不一致残差降不到容差内.msg里出现残差停滞措辞以官方文档为准表现为反复切步、最大迭代次数用尽【经验法则】判断雅可比是否与本构一致的客观办法在 UMAT 内对DSTR每个分量加 ε1e-8 再跑同一套更新逻辑把差分结果与写入的DDSDDE逐元素对比相对误差 1e-6 判不一致——官方讲义 Workshop 2Tangent Stiffness就是这一对照R6第 15 篇用DQN做同类检查。【最佳实践】显式雅可比要写成当前状态的切线不是初始弹性矩阵塑性段的DDSDDE应当是与应力更新算法return mapping返回映射同源的一致切线consistent tangent用连续体弹塑性切线替代一致切线在强硬化或大应变增量下会拖慢收敛。过渡句雅可比解决刚度对不对另一半契约是数据更新到哪儿——这一半更容易出错因为它不报错。1.3 增量末值与回滚STATEV是最常见的事故点【铁律】第 6 条UMAT 必须把应力与状态变量更新到增量末值end of increment否则下一个增量读到脏数据。增量 N 起点STRESS(起点) ← Abaqus 传入STATEV(起点) ← 上一增量写入的值 UMATSTRESS、STATEV 是输入兼输出必须就地写回增量末值 末点求解器检查残差 收敛 → 末值成为增量 N1 的起点 不收敛 → 回滚到增量 N 起点重算但 STATEV 不会自动回滚后果分两类。NSTATEV0本例风险较小只需保证STRESS被更新。NSTATEV0累积塑性应变、损伤变量、背应力时某次调用只更新一部分例如只在加载分支写、卸载分支在RETURN前没写STATEV会把中间值永久带进后续增量损伤从 0.2 跳到 0.35、累积塑性应变出现负值都是不可复现的漂移。【最佳实践】每次进入 UMAT 先把输入值按NSTATEV复制到本地数组再计算、最后统一写回卸载分支也必须显式写回哪怕写回原值。另一处易漏SEETEN增量末弹性应变能密度与SPD、SCD塑性耗散能、蠕变耗散能是可选输出不写不会报错但.dat与 ODB 的能量类输出会缺项第 9 篇。写入规则以官方文档为准。1.4 向量与矩阵索引约定6 个分量而不是 3×3NDIR正应力方向数短名NDI、NDNHR剪切分量数短名NSHR、NSTRESS总分量数短名NTENS三个数决定所有数组的形状约定是一维展开而不是张量单元类型NDIRNDNHRNSTRESS索引顺序本系列全篇统一3D 实体C3D8、C3D8R、C3D4306111, 222, 333, 412, 523, 631平面应力/平面应变单元213111, 222, 312壳单元S4、S4R积分点按层336111, 222, 333, 412, 523, 631轴对称单元314111, 222, 333, 412前 3 行的顺序是 Abaqus 长期沿用的工程约定逐单元类型的确切顺序以官方文档为准。两个推论DDSDDE是NSTRESS×NSTRESS方阵而非NDIR×NDIR剪应力索引从NDIR1起写成DO I1,NDIRNDNHR才完整覆盖。平面应力NSTRESS3时照搬 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 列为止 ──┘三种写法直接改变结果写法后果INTEGER KMAX, DV(2), DRV(2)从第 1 列开始整行被当注释声明消失后续引用按隐式类型处理I–N 开头视为 INTEGER续行标志写在第 5 列或第 7 列不是续行上一行在 72 列被截断数组声明不完整编译期直接报错语句超出第 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就是本构本身非线性材料的差别全在切线矩阵上。场景每增量调用次数经验法则说明线弹性、单积分点1雅可比精确一次迭代收敛一致切线与更新算法同源2~4常见残差快速下降连续体弹塑性切线近似4~10强硬化或大增量时更容易停滞DDSDDE置零或错索引无法收敛残差不降触发最大迭代次数或增量缩减到底调用次数在.sta与.msg里可读.sta给出每增量的迭代次数与增量大小变化.msg给出残差明细与 Force/Residual 对比第 10 篇的判定口径。【经验法则】跑 100 增量时平均迭代次数超过 5先怀疑雅可比而不是网格或容差。向量化vectorization是同一件事的另一面KMAX/KMSN是每次调用可处理的积分点向量长度名称以官方文档为准写成逐点循环、向量长度为 1 就拿不到批量调用的收益数组按(NTENS, KMAX)形态传入时模块级或SAVE变量必须按调用实例隔离否则线程并行会互相踩写。【最佳实践】先对齐单核结果再开并行。二、完整代码与逐行剖析2.1umat_iso_elastic.forFortran 固定格式复制即可编译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 ! 剪切分量从 NDI1 起算索引不可错位 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) 状态变量本模型 NSTATV0无状态变量若有必须显式写回增量末值 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, CONSTANTS2提供长度是NPROPSDVAL是求解器算好的派生量两者只读。把参数写在PROPS里而不是硬编码第 12 篇的参数表才能直接驱动 UMAT。这里的防御性检查给出PNEWDT0.5D0请求求解器把增量步长建议减半是用户唯一能影响增量控制的字段。第 (2) 段先清零、再写对角、再写非对角顺序不能颠倒——非对角分支只覆盖I, J ≤ NDI先写会覆盖对角线。λLamé 常数用Kbulk - 2G/3表达比E*NU/((1NU)(1-2NU))更稳NU → 0.5时Kbulk发散这个写法把不可压缩极限的风险摆出来而不是藏在一串代数里。第 (3) 段STRESS是输入兼输出用STRESS(I) STRESS(I) 2G·偏应变 K·体应变就地累加等于把起点值带到增量末。若写成STRESS(I) 新值必须确保STRESS仍是起点值——这正是没更新到增量末值事故的常见写法差异。第 (4) 段NSTATV0时这段看似冗余但保留显式写回是【最佳实践】后续若增加状态变量累积塑性应变、损伤这段就是必须存在的骨架。第 (5) 段能量项是可选输出把SEETEN/SSE写成累加而不是赋值是官方参数表对可选输出的常见约定是否必须累加、是否要求用更新前的应变以官方文档为准。【铁律】第 7 条编译前必须加载编译器环境。setvars.bat设置PATH与运行时库路径compile_fortran只在已加载的环境里找编译器未加载时编译阶段报找不到编译器、链接阶段报缺符号。2.2abaqus_v6.env与提交命令工作目录 C:\work\umat_caseumat_case.inp含 *USER MATERIAL, CONSTANTS2、modelUSER、 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 jobumat_case userumat_iso_elastic.for cpus1判定链路是否真走通三处证据按第 10 篇口径证据位置读法编译发生了目录出现umat_iso_elastic.obj与链接产物只有.for没有.obj 没走compile_fortran子程序被调用.log/.msg里的调用记录本例用WRITE打印一次都没打印 没挂上子程序user拼写或路径雅可比被使用.sta的迭代次数线弹性应为每增量 1 次每增量 1 次且残差不降回到第 (2) 段查DDSDDE2.3DDSDDE的预期值对照表用于自检E 200000.0、NU 0.3单位 mm–N–sG 76923.077、Kbulk 115384.615、λ K − 2G/3 25641.026。元素预期值N/mm²公式DDSDDE(1,1)、DDSDDE(2,2)、DDSDDE(3,3)129472.13λ 2G K 2GDDSDDE(1,2)、DDSDDE(2,3)、DDSDDE(1,3)及其转置78129.5K − 2G/3 λDDSDDE(4,4)、DDSDDE(5,5)、DDSDDE(6,6)76923.08G其余非对角正应力与剪切之间0各向同性无耦合对称性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 秒改成cpus8后反而变慢或报子程序不支持并行。根因按线程并行调用时模块级或SAVE共享变量被不同实例踩写向量化要求数组按(NTENS, KMAX)形态传入逐点循环写法拿不到收益。解法先对齐单核结果第 10 篇cpus判定再把共享状态改成调用内局部数组官方讲义的 threads 与 MPI 章节R6 的 Appendix 5、Workshop 5/6列出并行调用的写法约束具体限制以官方文档为准。四、动手练习练习 1用解析解验证DDSDD与应力更新容差。在E200000、NU0.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-9DDSDDE(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.env2.2 节版本与.for源码提交含*USER MATERIAL, CONSTANTS2的作业。判定标准目录出现.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一致性检查。本篇认知问题回显FAQQ1UMAT 的参数表里哪些是只读输入、哪些是必须写回的输入兼输出STRESS、STATEV、DDSDD、SEETEN、PROPS、DVAL、CEFIELDANSTEP、NDEC、DDELTA、TIME、CEFIELD、PROPS、DVAL是只读输入STRESS、STATEV、DDSDD必须写回增量末值SEETEN是可选输出。Q2材料雅可比material JacobianDDSDD/DDSDDE为什么必须显式提供留空后 Abaqus/Standard 拿什么组装单元刚度矩阵AAbaqus/Standard 不为 UMAT 数值求导留空时向量长度为 1 的路径把雅可比按单位矩阵处理E200000、ν0.3下DDSDDE(1,1)必须是 129472.13 N/mm²。Q3为什么STRESS与STATEV必须更新到增量末值end of increment回滚时脏数据如何污染下一个增量A本次写入的值就是下一增量的起点回滚只恢复起点数据而STATEV不自动回滚未写回的分支典型是卸载分支会把中间值永久带进后续增量。Q4NDIR/NDNHR/NSTRESS的索引约定是什么Fortran 固定格式fixed format列对齐为什么会直接改变结果A3D 实体单元按NDIR3、NDNHR0、NSTRESS6一维展开412、523、631DDSDDE是 6×6 方阵固定格式续行标志必须在第 6 列、语句从第 7 列到第 72 列。Q5compile_fortran、link_sl、link_exe为什么必须同源、setvars.bat与 ifx 在哪一步生效隐式为什么对切线刚度敏感、NTENS/KMSN向量化为什么是性能陷阱Acompile_fortran只产目标文件、运行时库在链接期才用到三处须同源且先加载setvars.batoneAPI 2025.0 起 ifort 移除须用 ifx隐式迭代方向由DDSDD决定KMSN/KMAX为 1 时失去向量化收益。
