ABAQUS蠕变问题计算流程：本构选型、参数设置与避坑指南

简介：《ABAQUS蠕变问题计算流程》是一份面向科研人员与工程师的有限元分析参考资料，系统梳理了材料长期受力下蠕变行为的仿真思路。内容以文档形式呈现，先介绍静态加载阶段的操作要点，再重点讲解蠕变计算流程，包括Power-law model等蠕变模型的选择、材料参数的拟合方法、蠕变分析步的设置以及蠕变应变容差（建议1E-4～1E-6）对增量步的影响，兼顾计算精度与效率。资源包为1个doc文档，共1个文件，大小约198KB，结构紧凑、重点突出，适合正在学习ABAQUS蠕变分析或希望快速上手建模的读者。该资源已有1157人学习。通过学习这份资料，读者可以理清从初始应力场到蠕变模拟的完整技术路线，掌握材料蠕变特性参数在CAE中的设置方法，理解不同积分算法与容差设置对结果稳定性的影响，从而减少反复试错，提升仿真分析效率。

1. ABAQUS蠕变问题计算流程：不是换一套材料参数那么简单

某公司有一个高温法兰模拟项目X，A同学图省事，从手册里抄了一组Norton蠕变参数直接填进材料卡，提交任务后算了整晚都不收敛，最后强行放开增量步，出来的等效蠕变应变超过200%。问题不在参数本身，而在“计算流程”——ABAQUS里的蠕变分析不是选个模型、填几个数就能跑完，它牵涉本构选型、单位制时间基准、增量步策略、应力如何传入子程序这几个环节，任何一环脱节，结果都废。这篇笔记把整套流程拆开讲，从本构原理到最小可复现算例，再到我踩过的几个坑，帮做结构强度分析的工程师少走夜路。

2. 蠕变本构与ABAQUS内置模型的选型逻辑

2.1 三阶段蠕变与工程构件的稳态假设

蠕变是材料在持续载荷和温度作用下随时间发生的非弹性变形。单轴蠕变曲线上通常分三段：减速段、稳态段、加速段，加速段很快进入断裂。工程上绝大多数高温承压构件的设计寿命都落在稳态段，减速段机制复杂、持续时间短，加速段不允许进入。所以ABAQUS蠕变问题计算流程的起点，就是选一个能描述稳态段的本构模型。

在ABAQUS里，材料模块的Creep子选项下有几个内置蠕变模型：幂律模型（power law）包含时间硬化形式和应变硬化形式，以及一些特殊模型。它们数学上长得像，但行为差异很大。很多人第一次做蠕变分析时直接选Time Hardening，因为界面默认就是它，但这恰恰是后续一系列问题的根源。

2.2 三个模型的数学形式与选型差异

Norton律是工程上最广泛使用的模型，参数少，单轴蠕变数据就能拟合出A和n。但Norton律有一个隐含假设：等效应力恒定。在ABAQUS里跑一个真实结构，应力会随时间松弛和重分配，如果还按Norton律算，蠕变变形势必偏大。尤其是螺栓、法兰这类预紧连接结构，应力松弛是主要失效模式，Norton律会给出偏危险的结果。

Time Hardening模型引入了时间t的幂函数，能描述蠕变率随时间衰减的过程，适合温度、载荷在时间轴上变化不大、整体应力水平维持恒定的构件。它的问题在于：同一计算里如果存在卸载和再加载，时间t已经很大，蠕变率会被严重低估。真实材料不记得钟表时间，它记得的是内部蠕变应变状态——这是Strain Hardening模型存在的根本原因。

Strain Hardening把时间依赖转换成蠕变应变的函数，对历史更敏感，数值上需要处理蠕变应变为零时的奇点。ABAQUS内置了对这类模型的下限保护，但我一般仍会建议在材料参数里配合蠕变截断值使用，避免程序在纯弹性阶段硬算。从我的项目经验看，结构件蠕变分析首选Strain Hardening，算例验证性分析用Norton律就够，Time Hardening反而容易在卸载工况上给出迷惑性结果。

2.3 材料参数输入表与三个必调参数

落到ABAQUS材料卡里，核心参数就三类：蠕变系数C1、应力指数C2、时间指数C3。注意ABAQUS没有内置单位制，全靠使用者保证一致性。用常见的N-mm-s单位制时，应力单位是MPa，长度是mm，时间是s；手册数据如果以小时h为单位，C1必须做时间基准换算，否则整个分析结果会放大3600倍——这个坑我在第5章专门展开。

C1和C2是从实验曲线上拟合出来的，C3才是需要工程师根据工况判断的调参项。很多资料直接写“时间硬化指数取-0.7”，但我的经验是：如果工况里有明显的应力松弛，C3取接近0的负值能保证收敛；更精细的做法是直接切换到Strain Hardening形式。三个必调参数里，除了上述本构参数，还有分析步里的时间长度——这是很多新手最容易忽略的，蠕变分析的时间不是计算伪时间，而是构件真实的服役时长。

3. ABAQUS蠕变计算的两条实现路线：内置参数与CREEP子程序

3.1 路线A：内置Creep参数的路子，适合材料数据完善的老牌号

内置模型的最大优势是简单：在材料模块里直接填参数，不需要写一行代码，也不涉及子程序编译和调试。适合三类场景：材料牌号有成熟蠕变数据、分析目的是结构校核与报告交付、等效应力在分析过程中变化幅度不超过20%。

我处理模拟项目X最初版本时，就是按这条路线走的。当时的温度恒定在550°C，载荷稳定，等效应力变化不大，用Time Hardening加一组拟合参数，几分钟就算完了。inp里对应的材料卡长这样：

*Material, name=CREEP-STEEL *Elastic 200000.0, 0.3 *Creep, law=TIME HARDENING 1.0e-14, 5.0, -0.5, 0.0

这里C1=1.0e-14对应单位时间蠕变率基数，C2=5.0是应力指数，C3=-0.5是时间指数，最后一个0.0是温度项指数，常温下不起作用。要注意：这个参数组合是针对模拟项目X的示例数据，真实工程必须用材料实验数据拟合，不能直接抄。

内置模型的参数敏感性很直观。C1每增大一个量级，同一时间点的蠕变应变就增大一个量级；C2则决定应力翻倍时蠕变率翻多少倍。调试时最有效的方法是做一个单轴单元验证算例，把CEEQ提取出来和解析解对拍，确认参数方向正确后再上真实结构。

3.2 路线B：CREEP用户子程序，应力路径复杂时的正解

当工况包含温度循环、多轴应力重分布、或者材料数据根本不是标准幂律曲线时，内置模型就不够用了。CREEP用户子程序真正的优势不是“复杂”，而是你可以自己访问状态变量、自己写蠕变增量公式，理论上可以模拟任意本构。

但这里有一个很多新手不知道的坑：CREEP子程序的接口里没有直接给出应力张量。很多人写好子程序，编译通过，任务提交后却发现蠕变应变恒为零，就是因为子程序里根本拿不到等效应力。我现在的做法是配一个USDFLD子程序，在每一增量步开头用GETVRM取出应力分量，算好米塞斯等效应力存入状态变量，CREEP子程序再读取这个状态变量来计算蠕变增量。

SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT, 1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT, 2 LAYER,KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP, 3 MATLAYO,LACCFLA) INCLUDE 'ABA_PARAM.INC' CHARACTER*80 CMNAME,ORNAME DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(2,3), 1 T(3,3),TIME(2),COORD(*),JMAC(*),JMATYP(*) CHARACTER*3 FLGRAY(15) DIMENSION ARRAY(15),JARRAY(15) REAL S11,S22,S33,S12,S13,S23 C CALL GETVRM('S',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP, 1 MATLAYO,LACCFLA) S11=ARRAY(1) S22=ARRAY(2) S33=ARRAY(3) S12=ARRAY(4) S13=ARRAY(5) S23=ARRAY(6) C 计算米塞斯等效应力并存入第一个状态变量 STATEV(1)=SQRT(((S11-S22)**2+(S22-S33)**2+(S33-S11)**2 1 +6.0*(S12**2+S13**2+S23**2))/2.0) RETURN END

对应的CREEP子程序读这个状态变量：

SUBROUTINE CREEP(DECRA,DESWA,STATEV,TIME,DTIME,DTEMP,DTEMPME, 1 PREDEF,DPRED,CMNAME,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC, 2 NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA) INCLUDE 'ABA_PARAM.INC' CHARACTER*80 CMNAME DIMENSION DECRA(5),DESWA(5),STATEV(*),TIME(2),PREDEF(*), 1 DPRED(*),COORD(*),JMAC(*),JMATYP(*) C A=1.0D-14 XN=5.0D0 SIG=STATEV(1) C DECRA(1)=A*SIG**XN*DTIME DESWA(1)=A*SIG**XN RETURN END

DECRA(1)是当前增量步的等效蠕变应变增量，DESWA(1)是等效蠕变应变率。子程序里A和XN的取值只是示意，真实使用必须换成材料数据拟合结果。提交时用命令：

abaqus job=creep_job user=creep.for cpus=4 -interactive

使用子程序时我一般会先做一个单单元模型跑通调试，确认状态变量传递无误后，再换成完整结构网格。否则完整模型动辄几小时起步，子程序里一个数组越界就要浪费一整天。

3.3 我的选型判断：什么情况打死不写子程序

经常有工程师问：我要学子程序吗？我的回答是：先问工况需不需要。如果结构简单、温度恒定、只需要最大蠕变变形量来校核，内置模型完全够用，写子程序是自找麻烦。子程序适合三类人：材料本构需要自定义的、温度场随时间剧烈变化的、以及做科研需要参数反演的。

另一个判断维度是团队维护成本。子程序不是写完就完了，换版本、换机器都要重编译，Fortran和C的接口细节各家还不一样。如果项目交付周期只有两周，内置模型能覆盖需求，我强烈建议用内置。子程序是手段不是目的，它不是蠕变计算流程里必需的一环。

4. 实操：从材料卡到结果输出的完整ABAQUS设置

4.1 开始前的三项检查：单位制、对称性、加载顺序

蠕变分析最怕的是模型建好了、材料卡也写了，结果方向全错。我在每次建模前固定做三项检查，这三项检查大约十分钟，能规避后续三小时的排查。

第一项是单位制一致性。ABAQUS没有单位制，一套模型里只能有一套自洽单位。我常用N-mm-s单位制，应力单位是MPa，密度用t/mm³，时间用秒。如果数据手册里蠕变率单位是1/h，C1换算成1/s时必须除以3600；如果模型时间单位是h，那弹性模量、密度、蠕变参数全都要按小时体系重新折算。混用单位是新手翻车第一来源。

第二项是几何对称性。蠕变分析计算成本比弹塑性分析高很多，因为增量步受时间尺度限制。能取1/4模型就不要取整圈，能用轴对称就不要建三维。模拟项目X里的法兰，我当时用了二维轴对称模型，计算时长从整圈三维模型的4小时降到20分钟。

第三项是加载顺序。蠕变率对应力路径极其敏感，先加载后升温还是先升温后加载，结果可以差出百分之几十。常见做法是先用Static, General建立载荷和温度场，再创建一个Visco分析步让蠕变在这个步内发生。如果温度场本身随时间变化，必须先跑传热分析，再把温度场作为预定义场导入结构分析。

4.2 用内置模型跑通最小算例的完整配置

以一个单轴拉伸蠕变算例为例，我习惯手写inp而不是点界面，因为能完整控制每一个关键词。最小算例包括一个单元、一组弹性参数、一组蠕变参数和一个Visco分析步：

*Node 1, 0.0, 0.0 2, 1.0, 0.0 3, 1.0, 1.0 4, 0.0, 1.0 *Element, type=CPE4 1, 1, 2, 3, 4 *Material, name=CREEP-STEEL *Elastic 200000.0, 0.3 *Creep, law=STRAIN HARDENING 1.0e-14, 5.0, 0.0, 0.0 *Step, name=CreepStep, inc=1000 *Visco, cetol=0.005 1000000.0 *Boundary 1, 2, 2 4, 2, 2 1, 1, 1 2, 1, 1 *Cload 2, 1, 1000.0 3, 1, 1000.0 *Output, field, freq=10 *Node Output U *Element Output CEEQ, S *End Step

这个例子里，材料卡在原有弹性参数基础上加了蠕变参数，分析步长度设为1000000秒，约11.5天，cetol=0.005是蠕变容差，控制每个增量步内蠕变应变增量的收敛精度。边界条件约束了刚体位移，上表面施加总载荷1000N，等效于一个单轴拉伸应力状态。

参数说明：C3=0表示当前用的是退化后的Norton律，C1=1.0e-14在应力水平100MPa时蠕变率约1e-9每秒，积累一天后应变约0.01%，处于合理范围。算完后提取单元输出CEEQ和S，检查CEEQ是否随时间线性增长，如果是，说明参数设置方向正确。

4.3 CREEP子程序的编译与任务提交

使用子程序时，建模和材料定义与内置模型一致，区别在于材料卡里不再写*Creep关键词，而是留空，本构行为完全由子程序通过材料点调用。提交命令在4.2里给过，这里补充两个环境层面的细节。

第一个是环境变量设置。ABAQUS调用子程序需要Fortran编译器，不同版本对应的编译器版本不同，装错会导致链接阶段报错。我一般在提交前先跑一个最简单的creep子程序算例，验证编译环境通不通，不直接在正式模型上试。

第二个是子程序文件的编码问题。Windows下用记事本写的Fortran文件如果带中文注释，编译时可能出现乱码报错。我的习惯是子程序文件全部用英文路径、英文文件名，注释写英文或拼音，彻底避开编码坑。

abaqus job=creep_job user=creep.for cpus=4 -interactive

提交后观察log文件，如果出现“CREEP subroutine is active”之类的提示，说明子程序已成功挂载。如果log里没有任何相关字样，检查是否忘了在job命令里加user参数。

4.4 增量步与收敛参数怎么给

蠕变分析对增量步的敏感度远超弹塑性分析，因为蠕变增量本身和DTIME成正比。一个常见的错误是沿用静力分析的默认增量步设置，初始增量步给到总时间的10%，结果前三步就不断报收敛错误。

这些参数可以在inp里直接写，也可以界面操作。我的经验是：初始增量步宁可给得小一点，让程序自己往大调，也比一开始就给大了然后反复回退效率高。cetol从默认值调到0.005，增量步总数通常增加20%~30%，但结果的振荡明显减少。

5. 蠕变计算避坑指南：五个让我翻车的现场

5.1 时间单位翻车：一年变一天

现象：手册数据给的是每小时蠕变率，模型单位制是秒，直接把C1填进去，出来的蠕变应变比预期大了3600倍。

原因：C1的量纲里包含时间单位1/h，填进以秒为单位的模型后，程序把它解释成1/s，相当于放大了3600倍。

解决：把C1除以3600再做换算。更稳妥的方式是先确认目标分析步的时间单位，再把所有材料和载荷参数统一到同一套单位制下。我现在做蠕变分析前，固定先在纸上列一套单位清单，把所有参数的物理单位写清楚再动手。

5.2 子程序里取不到应力，蠕变应变为零

现象：子程序编译通过、任务正常算完，CEEQ云图全为0，查log文件没有任何错误。

原因：CREEP子程序接口参数里没有应力张量，很多人直接写SIG=STATEV(1)，但STATEV(1)从未被赋值，默认就是0，蠕变增量自然为0。

解决：按照3.2节的方案，用USDFLD配合GETVRM取应力，存入状态变量后再由CREEP读取。调试时在这个环节加两个辅助状态变量输出，看看拿到的米塞斯应力有没有量级异常。

5.3 应力松弛算不出来，螺栓预紧力不降

现象：用Time Hardening模型算螺栓蠕变，分析结束时预紧力纹丝不动，和实验趋势完全相反。

原因：Time Hardening模型把蠕变率绑定在时间t上，螺栓的应力虽然松弛了，但时间一直在走，蠕变率被模型“冻结”在初始水平，无法反映应力松弛后的减速效应。

解决：换成Strain Hardening模型或自定义子程序。这个案例让我彻底放弃了在应力松弛场景里使用Time Hardening模型的做法，现在结构件分析一律先用Strain Hardening跑一版。

5.4 初始增量步给太大，任务一直不收敛

现象：任务提交后log文件里全是“Time increment is reduced”，增量步反复回退，最终分析步走不动。

原因：初始增量步设为总时间的10%，第一步蠕变增量超过程序容差，程序不断缩小步长但始终无法稳定。

解决：初始增量步改到总时间的1/1000，cetol设到0.005，任务立刻顺畅。记住一个规律：蠕变分析里增量步的尺度由蠕变率决定，不是由网格尺寸决定。

5.5 温度场加载顺序反了，结果差30%

现象：同一个模型，先加载后升温得到的结果和先升温后加载差出30%以上，两边都没报错，但谁都说不清哪个对。

原因：蠕变率对温度高度敏感，先加载后升温意味着蠕变在低温低应力阶段就开始积累，而先升温后加载则是先膨胀再受力，应力路径完全不同。

解决：先准确判断工况是温度恒定还是温度变化，恒定温度就直接在所有增量步里保持同一温度场，变化温度必须把传热分析的结果按时间节点导入。不要让结构分析去拟合一个未经传热计算的温度场。

6. 单轴蠕变解析验证法与我的收尾习惯

6.1 用Norton律解析解与有限元结果对拍

验证蠕变计算流程是否正确，最直接的方式是做一个单轴单元算例，把CEEQ提取出来，与Norton律解析解对比。恒定应力σ₀下，Norton律的蠕变应变表达式是ε_cr = A·σ₀^n·t。

误差随累积应变略微增大，是因为数值积分对时间增量的离散误差在积累，正常情况下几种误差都应在2%以内。超过3%就先查单位制，再查增量步设置。

6.2 验证时怎么控制变量

一个算例里有太多变量时，出了问题根本找不到根源。我现在对拍的固定做法是：先用纯弹性算例验证单位制，再用带蠕变参数的算例验证本构模型，最后才让子程序参与测试。每一步只改一个变量，测试时间三分钟，换来的是后续排查效率大幅提升。

另一个控制变量的维度是结果输出频率。不要每增量步都输出，默认每10个增量步输出一次足够看趋势，需要精细分析局部时再单独调输出频率。

6.3 我现在的收尾习惯

每个蠕变模型跑完后，我固定做两件事：第一件是看一眼log文件里每增量步的蠕变增量值，如果出现连续回退迹象就回查初始步长；第二件是把CEEQ云图和解析解或实验数据对拍，确认数量级没问题再进下一轮。这个习惯帮我挡掉了至少三次单位制翻车事故。如果你只记住一句话：蠕变计算流程里，单位制一致性和增量步策略比本构模型的精确度更容易毁掉结果，先跑通最小算例，再上完整模型，这个顺序不能乱。希望帮到你。

本文还有配套的精品资源，点击获取