creep蠕变基础知识
一种蠕行控制方法叫什么
一种蠕行控制方法叫什么引言蠕行（creep）是材料学中常见的现象，指的是材料在长时间受到持续作用力时发生的无反弹性变形。
蠕行现象在许多领域都存在，如工程结构设计、材料加工和地质工程等。
为了控制蠕行现象，人们提出了许多方法。
本文将介绍一种“蠕行控制方法”。
蠕行控制方法的原理目前通常使用的蠕行控制方法主要分为两大类：应力降低法和材料改性法。
- 应力降低法：通过减小材料的受力强度，使其承受更小的应力从而降低蠕变速率。
例如，在工程结构设计中，可以采用增大材料截面尺寸、增加结构支撑等方式来降低结构的应力水平，从而减缓蠕变速率。
- 材料改性法：通过改变材料的组织结构，调控材料的形变行为，从而控制蠕行现象。
例如，通过合金化、陶瓷增强、纤维增强等手段来改变材料的微结构，增强材料的抗蠕变性能。
然而，以上方法都有其局限性。
应力降低法需要在设计初期耗费大量的时间和精力，并且难以在运营期进行调整；材料改性法则需要以提高材料强度为代价，无法得到完美的控制效果。
因此，人们需要一种更加高效和实用的蠕行控制方法。
一种新的蠕行控制方法：多物理场耦合法近年来，随着科学技术的发展和理论模拟能力的提高，人们提出了一种新的蠕行控制方法，即多物理场耦合法。
这种方法基于多物理场的相互作用，通过控制不同物理场的耦合行为来实现对蠕行现象的控制。
多物理场指的是同时存在于同一介质或系统中的多个物理场，例如电场、热场、磁场、应力场等。
而耦合行为则是指这些物理场之间的相互影响和相互作用。
通过调控不同物理场的大小、分布和时间演化等参数，可以改变材料内部的微观状态，从而影响材料的蠕变行为。
多物理场耦合法的具体实施方式可以通过数值模拟、实验验证和实际应用等不同途径来进行。
其中，数值模拟是其中的重要手段之一。
通过建立适当的物理模型和数值算法，可以模拟不同物理场的耦合行为，并通过改变参数来控制蠕行现象。
实验验证可以通过实验室试验和实际工程应用来进行，通过对不同物理场的加载和监测，验证多物理场耦合法对蠕行的控制效果。
材料的蠕变机制了解材料在高温下的变形机制
材料的蠕变机制了解材料在高温下的变形机制材料的蠕变机制：了解材料在高温下的变形机制材料的蠕变（Creep）是指在高温和应力作用下，材料会逐渐发生形变和软化的现象。
蠕变机制的了解对于材料在高温环境下的设计和应用非常重要。
本文将探讨材料在高温下的变形机制，以及蠕变现象对材料性能的影响。
1. 引言随着工程领域对材料性能要求的提高，对于材料在极端条件下的变形行为的研究变得愈发重要。
高温下材料的蠕变现象是一项关键的考察对象，本文旨在深入探讨材料在高温环境下的变形机制，并分析其对工程应用的意义。
2. 材料在高温下的变形机制2.1 晶体滑移晶体滑移是一种经典的材料变形机制，在高温条件下尤其突出。
晶体原子通过在晶格面上滑移来改变位置，从而导致材料的塑性变形。
高温会增加晶体内的原子动能，使得晶体滑移更容易发生。
2.2 相互扩散在高温下，材料中的原子会因为动能提高而显示出相互扩散的特性。
相互扩散会导致晶界和晶体内部的松弛，进而引起材料的蠕变变形。
2.3 晶体再结晶晶体在高温和应力作用下可能发生再结晶，即原本晶体中的晶粒重新排列和重组。
这种再结晶过程也是材料蠕变的一个重要机制。
3. 蠕变对材料性能的影响3.1 强度和韧性的降低材料蠕变使得材料在高温下的强度和韧性均降低。
原子的扩散会导致材料晶界的松弛，晶体结构的破坏以及晶粒的再结晶，这些因素会使材料的强度降低。
3.2 变形速率蠕变行为与应力的大小、温度的高低以及时间的长短密切相关。
蠕变变形速率随温度的升高而增加，随应力的增加而加快。
这对一些高温工程应用中材料的耐久性和设计有着重要影响。
4. 变形机制的研究与应用4.1 实验方法通过高温实验设备对材料进行蠕变试验，可以模拟材料在高温下的变形行为。
通过测量变形曲线、变形速率和变形温度等参数，可以获取材料的蠕变特性。
4.2 材料改进与设计通过深入研究材料的蠕变机制，可以针对不同的应用需求改进和设计材料。
例如，在航空航天领域，对发动机叶片等高温结构部件材料的蠕变机制进行研究，以提高其抗蠕变性能。
蠕变理论及寿命评估
蠕变断裂寿命(持久寿命)成为寿命预测中最关键的 因素之一. 因此相关设计规范规定高温机械设备必须采 用材料的长期持久强度来进行设计. 应用最广泛的蠕变 持久寿命预测方法, 是以拉森- 米勒法(简称L - M 法)为 代表的时间- 温度参数。
Your company slogan
发生蠕变所需的应力都可以很低，甚至远低于高 温屈服强度。而发生蠕变的温度则是相对的，蠕变在 低温下也会产生，但只有在约比温度（T/Tm）高于0.3 时才较显著，所以通常称为高温蠕变。
碳钢温度超过300°C，合金温度超过400°C时， 就必须考虑蠕变影响。
Your company slogan
蠕变曲线
Your company slogan
（二）冶炼工艺的影响 各种耐热钢及高温合金对冶炼工艺的要求较高，
因为钢中的夹杂物和某些冶金缺陷会使材料的持久 强度极限降低。高温合金对杂质元素和气体含量要 求更加严格，常存杂质除S、P外，还有铅、锡、砷、 锑、铋等，即使其含量只有十万分之几，当其在晶 界偏聚后，会导致晶界严重弱化，而使热强性急剧 降低，并增大蠕变脆性。
Your company slogan
5）将每一个样品的试验温度T、应力σ和老化程度即 老化因子Ca代入 T(C+logtr)＝f(σ)+Ca中即可得到该试样 的蠕变使用寿命。
Your company slogan
Your company slogan
Your company slogan
Your company slogan
(二)扩散蠕变
在较高温度（(T／Tm)远超过0.5)下的一种蠕变变形机理。
在高温下大量原子和空位定向移动造成的。
土层蠕变的计算公式
土层蠕变的计算公式土层蠕变（Creep）是指土壤在长时间荷载作用下，由于颗粒重排和应力重新分布，导致土层体积发生持续变形的现象。
土层蠕变可以导致土体的应力松弛和稳定性变化，对工程结构的设计和施工具有重要影响。
1.永久变形速率法永久变形速率法是通过测定土壤的蠕变速率来计算土壤蠕变变形的方法。
根据实验数据统计分析，蠕变速率与应力水平、时间的关系可以用以下经验公式表示：ε=ασ^nt^m式中，ε为土壤蠕变变形，σ为荷载应力，t为荷载时间，α为经验系数，n和m为经验指数。
2. Kelvin弹簧模型法Kelvin弹簧模型是一种经典的线性黏弹性模型，将土体视为由弹簧和阻尼器组成的系统，并假设弹簧和阻尼器之间有一定的松弛时间。
Kelvin弹簧模型法通过描述弹簧和阻尼器之间相互作用的力学关系，来计算土壤的蠕变变形。
Kelvin弹簧模型中，弹簧代表土体的弹性变形，阻尼器代表土体的黏性变形。
Kelvin弹簧模型的计算公式可以表示为：ε=(1-e^(-λt))/λσ式中，ε为土壤蠕变变形，σ为荷载应力，t为荷载时间，λ为Kelvin模型的弹簧常数。
3. Burgers弹簧模型法Burgers弹簧模型是在Kelvin模型基础上发展出的一种复杂的弹簧模型，将土壤的蠕变变形分解为由弹性分量和黏性分量组成的两个部分。
Burgers弹簧模型法通过分析弹簧和阻尼器之间的相互作用，来计算土壤的蠕变变形。
Burgers弹簧模型中，弹性分量相当于Kelvin模型中的弹簧，而黏性分量相当于Kelvin模型中的阻尼器。
Burgers弹簧模型的计算公式可以表示为：ε=ε_e+ε_vε_e=C_1(1-e^(-λ_1t))/λ_1σε_v=C_2(1-e^(-λ_2t))/λ_2σ式中，ε为土壤蠕变变形，ε_e为土壤的弹性变形，ε_v为土壤的黏性变形，C_1和C_2为常数，λ_1和λ_2为Burgers模型的弹簧常数。
这些土层蠕变的计算公式是经过实验验证和统计分析得出的经验公式，可以较为准确地预测土壤蠕变变形的发展趋势。
聚合物的粘弹性之蠕变分析
（2）恒速蠕变（又称稳态蠕变）阶段－bc。 这一阶段的特点是蠕变速率几乎保持丌变。 一般所指的金属蠕发速率，就是以这一阶段 蠕发速率表示。
（3）加速蠕变阶段－cd 在此阶段随着时间的延长，蠕发速率逐渐 增大， 至d 点产生蠕发断裂。
温度不应力对蠕变曲线的影响： 在应力较小或温度较低时，蠕变第二阶 段持续时间较长，甚至丌出现第三阶段。 反之，蠕发第二阶段很短，甚至消失， 很短时间内就断裂。
蠕变断裂断口的宏观特征为： （1）在断口附近产生塑性发形， 在发形区域附近有很多裂纹，使 断裂机件表面出现龟裂现象。
（2）由于高温氧化，断口表面往往被 一层氧化膜所覆盖。
（3）蠕变断裂微观特征：为冰糖状花 样的沿晶断裂形貌。
蠕发现象 蠕发发形机理 蠕发机制图谱 蠕发断裂 抗蠕发材料设计
• 蠕变第一阶段：由于蠕发发形逐渐产生应变 硬化，使位错源开动的阻力及位错滑秱阻力 增大，使蠕发速率丌断降低。 • 蠕变第二阶段：因应发硬化収展，促迚动态 回复，使金属丌断软化。当应变硬化不回复 软化达到平衡时，蠕发速率为一常数。
(二)扩散蠕变
(二)扩散蠕变
扩散蠕变：是在较高温度(约比温度(T／Tm) 进超过0.5)下的一种蠕发发形机理。 它是在高温下大量原子和空位定向移动造成 的。 在丌叐外力情冴下，原子和空位的秱动无方 向性，因而宏观上丌显示塑性发形。 但当叐拉应力σ作用时，在多晶体内产生丌 均匀的应力场。
1、刃型位错的攀移（图）
2、交滑移
b
b
b
b b
刃位错攀移克服障碍的几种模型： 可见，塞积在某种障碍前的位错通过热 激活可以在新的滑秱面上运动（a）， 戒不异号位错相遇而对消（b），戒形 成亚晶界（c），戒被晶界所吸收（d）。
ANSYS蠕变
各个论坛均有关于creep的问题，但似乎都没有弄明白，为此这里专门再研讨研讨。
初步理解：V1.Creep在金属材料中被称为蠕变，而在混凝土材料中被称为徐变，二者规律各异，但本质是一样的，所以讨论中统一使用“徐变”。
V2.Ansys中给出了Creep两种方式，即隐式和显式的;隐式通过tb命令的tbopt指定方程号，显式则通过tbopt=0来调用，这里仅讨论隐式的。
V3.通过ansys给出的方程可以看出，εcr=f(σ,ε,t)。
依照混凝土的相关知识，一般的εcr=σ/E×υ(t,τ)=F(σ,t)，显然ansys的徐变能够包罗混凝土的徐变问题。
或者讲用ansys的creep是能够解决混凝土徐变问题的。
尽管其给出的徐变方程与我们常用的不同，但似乎通过数值曲线的模拟，然后再确定哪些系数应该是可行的；实在不行不是有用户自定义的吗，所以应该不是问题。
初步问题:q1.在隐式徐变中，1～12对应的哪个是"徐变"，而哪个是"徐变率"呢(对时间)?好像help中的与中文手册中的不一样！q2.通过徐变方程可以看出，可以是变化的应力，是否用ansys就可解决变应力时的徐变呢？果真这样的话就太好了，因为常规计算这东西太麻烦了。
q3.对混凝土材料究竟选择哪个更合适呢？或者更接近呢？(可以规范的υ(t,τ)而论)。
q4.因为混凝土徐变受环境、材料、应力、时间等因素的影响，所以对于不同的混凝土其徐变系数是不同的，这里不必具体讨论。
事实上由于你没有足够的时间来进行徐变试验，尽管你知道与规范的不同，恐怕也得采用规范给出的徐变系数公式。
q5.给出的方程中与时间有关的话，程序中也要用time，此时是否方程中time的单位与apdl 中time的单位相同即可？只要你统一了就行？q6.混凝土徐变是采用主徐变方程？还是主＋二期徐变方程？初步例题：结合ansys手册中的例子，给出如下例题，以抛砖引玉！!两端固定的受压板finish/clearb=100h=30n1=20n2=5/prep7et,1,42mp,ex,1,3e4mp,prxy,1,0.2tb,creep,1,,,6 !定义6号方程tbdata,1,3.728e-6,1,-0.84,0rect,0,b,0,hlsel,s,loc,y,0lesize,all,,,n1lsel,s,loc,x,0lesize,all,,,n2mshkey,1mshape,0,2damesh,allfinish/solunsel,s,loc,x,0d,all,allnsel,s,loc,x,bd,all,ux,-0.1allsel,allbfunif,temp,900time,7300 !20yearkbc,1nsubst,500outres,all,allallsel,allsolvefinish/post1pldips,1/post26nsol,2,100,u,yplvar,2esol,3,55,,s,xprod,4,3,,,,,,-1plvar,4esol,5,55,,epel,xprod,6,5,,,,,,-1plvar,6esol,7,55,,epcr,xprod,8,7,,,,,,-1plvar,8!======================================!该例子中设：!Φ(t)=0.7*t^0.16(t单位为天)!εcr=σ/E*0.7*t^0.16=2.33e-5*σ*t^0.16!===================================================== ==============================e1.从图中可以看出，εcr7300≈0.92e-3=920με,而初始弹性应变为1000με左右,故剩余80με则应力在2.4MPa左右，与elem55-σx相当。
creep蠕变基础知识
蠕变模型将flac3d 的蠕变分析option 进行了简单的翻译，目的是为了搞清楚蠕变过程中系统时间是如何跟真实时间对应的。
2.1 简介Flac3d 可以模拟材料的蠕变特性，即时间依赖性，flac3d2.1提供6种蠕变模型： 1. 经典粘弹型模型 model viscous 2. model burger 3. model power 4. model wipp 5. model cvisc6. powe 蠕变模型结合M-C 模型产生cpow 蠕变模型（model cpow ）7. 然后WIPP 蠕变模型结合D-P 模型产生Pwipp 蠕变模型（model pwipp ）； 8 model cwipp以上模型越往下越复杂，第一个模型使用经典的maxwell 蠕变公式，第二个模型使用经典的burger 蠕变公式，第三个模型主要用于采矿及地下工程，第四个模型一般用于核废料地下隔离的热力学分析，第五个模型是第二个模型的M-C 扩展，第六个模型是第三个模型的M-C 扩展，第七个模型是第四个模型的D-P 扩展，第八个模型也是第四个模型的一种变化形式，只是包含了压硬和剪缩行为。
2.2蠕变模型描述2.2.1只介绍经典粘弹型模型即maxwell 蠕变公式牛顿粘性的经典概念是应变率正比于应力，对于粘性流变应力应变关系以近似于弹性变形的方式发展。
粘弹型材料既有粘性又有弹性，maxwell 材料就是如此，在一维空间它可以表示为一根弹簧（弹性常数κ）连接一个粘壶（粘性常数η），它的力-位移增量关系可以写成：ηκμFF+=∙∙（2.1）式中∙μ是速度，F 是力，设力的初始值为F ，增量值为F '经过一个t ∆时间步，式（2.1）可以写成ηκμ2F F tFF t+'+∆-'=∆∆ （2.2）这就是中心差分公式。
解F '得21)(C C F F μκ∆+='（2.3）ηκ211t C ∆-=ηκ2112tC ∆+=式（2.3）写成偏应力与应变增量的关系()212C G C dij dij dij εσσ∆+=（2.4）上式中：ijij ij dij δεεε∆-∆=∆31ij ij ij dij δσσσ31-=η211t G C ∆-=η2112t G C ∆+=这里，ij ε∆为应变增量张量分量，ij σ为初始应力张量分量，G 为剪切模量。
聚合物的蠕变,应力松弛,滞后和内耗
在聚合物科学中，蠕变、应力松弛、滞后和内耗是与聚合物材料的力学行为相关的术语。
蠕变（Creep）：蠕变是指在持续受到应力的情况下，聚合物材料会随着时间的推移发生形变。
蠕变是一个时间依赖的现象，即应力施加时间越长，形变越明显。
蠕变通常由于聚合物链的重新排列和滑移引起，导致聚合物结构的变化。
蠕变是一种可逆现象，当去除应力时，材料会回弹至原始形态。
应力松弛（Stress relaxation）：应力松弛是指在一定的应变条件下，聚合物材料所受的应力会随着时间的推移逐渐减小。
这是因为聚合物链在应力作用下发生重排，使得材料内部的应力逐渐减小。
与蠕变不同，应力松弛通常是在给定应变条件下观察到的。
滞后（Hysteresis）：滞后是指聚合物材料在循环加载和卸载的过程中，其应力和应变之间存在的差异。
在加载期间，聚合物会表现出较高的应力响应，但在卸载期间，应力并不完全消失。
这种差异是由于聚合物链的结构重排和能量耗散引起的。
滞后现象常见于高分子弹性材料，如弹簧和橡胶。
内耗（Internal friction）：内耗是指聚合物材料在受力或形变时，由于分子内部摩擦和相互作用而产生的能量损耗。
内耗可以导致材料的能量耗散和温升。
聚合物材料的内耗通常与材料的分子结构、聚合度和温度等因素有关。
内耗在聚合物的动态力学性能和阻尼特性中起着重要作用。
这些现象在聚合物工程和材料科学中具有重要的应用。
研究和了解聚合物的蠕变、应力松弛、滞后和内耗行为对于设计和开发具有特定力学性能和可靠性的聚合物制品非常重要。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
蠕变模型将flac3d 的蠕变分析option 进行了简单的翻译，目的是为了搞清楚蠕变过程中系统时间是如何跟真实时间对应的。
2.1 简介Flac3d 可以模拟材料的蠕变特性，即时间依赖性，flac3d2.1提供6种蠕变模型： 1. 经典粘弹型模型 model viscous 2. model burger 3. model power 4. model wipp 5. model cvisc6. powe 蠕变模型结合M-C 模型产生cpow 蠕变模型（model cpow ）7. 然后WIPP 蠕变模型结合D-P 模型产生Pwipp 蠕变模型（model pwipp ）； 8 model cwipp以上模型越往下越复杂，第一个模型使用经典的maxwell 蠕变公式，第二个模型使用经典的burger 蠕变公式，第三个模型主要用于采矿及地下工程，第四个模型一般用于核废料地下隔离的热力学分析，第五个模型是第二个模型的M-C 扩展，第六个模型是第三个模型的M-C 扩展，第七个模型是第四个模型的D-P 扩展，第八个模型也是第四个模型的一种变化形式，只是包含了压硬和剪缩行为。
2.2蠕变模型描述2.2.1只介绍经典粘弹型模型即maxwell 蠕变公式牛顿粘性的经典概念是应变率正比于应力，对于粘性流变应力应变关系以近似于弹性变形的方式发展。
粘弹型材料既有粘性又有弹性，maxwell 材料就是如此，在一维空间它可以表示为一根弹簧（弹性常数κ）连接一个粘壶（粘性常数η），它的力-位移增量关系可以写成：ηκμFF+=∙∙（2.1）式中∙μ是速度，F 是力，设力的初始值为F ，增量值为F '经过一个t ∆时间步，式（2.1）可以写成ηκμ2F F t F F t +'+∆-'=∆∆ （2.2） 这就是中心差分公式。
解F '得21)(C C F F μκ∆+=' （2.3）ηκ211tC ∆-= ηκ2112t C ∆+=式（2.3）写成偏应力与应变增量的关系()212C G C d ij d ij d ij εσσ∆+=（2.4）上式中：ij ij ij d ij δεεε∆-∆=∆31ij ij ij d ij δσσσ31-=η211tG C ∆-= η2112t G C ∆+=这里，ij ε∆为应变增量张量分量，ij σ为初始应力张量分量，G 为剪切模量。
对于应力应变的体积分量，假设体积变化不受流变影响。
kk kk iso K εσσ∆+=31（2.5）K 为体积模量，最终的应力张量就表示为偏量和球量之和：ij iso dij ij δσσσ+=该模型要求输入材料属性K 和G 和粘性常数，在剪应力作用下材料连续流变，在球应力作用下它表现为弹性。
2.3 flac3d 解流变问题 2.3.1简介流变模型和flac3d 其他模型最大的不同在于模拟过程中时间概念的不同，对于蠕变，求解时间和时间步代表着真实的时间，而一般模型的静力分析中，时间步是一个人为数量，仅仅作为计算从迭代到稳态的一种手段来使用。
2.3.2 flac3d 的蠕变时间步长对于蠕变等时间依赖性问题，flac3d 容许用户自定义一个时间步长，这个时间步长的默认值为零，那么材料对于粘弹性模型表现为线弹性，对于粘塑性模型表现为弹塑性。
（命令set creep off 也可以用来停止蠕变计算。
）这可以用来在系统达到平衡后再开始新的蠕变计算。
蠕变公式中包含时间，所以计算中时间步长对程序响应有影响。
虽然用户可以对时间步进行设置，但并不是任意的。
蠕变过程由偏应力状态控制，从数值计算的精度来讲，最大蠕变时间步长可以表示成材料粘性常数和剪切模量的比值：Gt crη=∆maxFor the power law . . . . .省略 For the WIPP law . . . . . 省略For the Burger-creep viscoplastic model, 上面方程应该写成： tmax = min ( ηK/GK ，ηM/GM)上标K 和M 分别代表Kelvin 和Maxwell 。
蠕变压缩的时间限制包括系统体积反应，并且估计为粘性和体积模量的比值。
粘性可以表示为σ和体积蠕变压缩速率的比值。
蠕变压缩的最大蠕变时间步如下：(2-127) 建议利用FLAC3D 作蠕变分析开始时所采用的蠕变时间步，比根据上式算得的时间tmax 小两到三个数量级。
通过调用SET creep dt auto on ，可以利用自动时间步自动调整（见2.3.3）。
作为一项规则，时间步的最大值（SET creep maxdt ）不能超过tmax 。
参见2.5节例题。
用来计算tmax 的应力σ大小，可由蠕变开始之前的初始应力状态决定。
同样，σ作为von Mises 不变量，可以用FISH 函数计算。
涉及体积变化响应的蠕变分析，其最大时间步长可以表示成材料粘性常数和体积模量的比值，这里粘性常数就是平均应力σ和蠕变体应变率cv ∙ε的比值。
一般flac3d 推荐使用的初始蠕变时间步长比最大时间步长（由上述公式计算得到的）约小2到3个数量级。
如果使用set creep dt auto on 命令，那么程序将自动调整蠕变的时间步长，这在下一小节2.3.3中介绍，同样应当记住通过命令（set creep maxdt）设置的最大蠕变时间步不能超过cr t。
max2.3.3自动调整蠕变时间步长用户可以设置蠕变时间步为一个常数值，也可以使用set creep dt auto on命令自动调节。
如果时间步长自动变化，那么当最大不平衡力超过某一阀值时，它就会减小；当最大不平衡力小于某一水平时它就会增大。
系统将该阀值定义为最大不平衡力和平均节点力的比值。
对于即将求解的问题，先只考虑弹性效应，通过观察接近初始平衡状态时的不平衡力就可以大概知道整个问题的不平衡力大小。
在有些算例中，又会尽量避免对时间步长的连续调整，这样的话，在某一次时间步长调整之后，可以定义一个“延迟阶段”（比如100步），在这一“延迟阶段”内不再发生进一步的时间步长调整。
一般来说，为了适应开挖之类的瞬变，时间步长一开始是一个较小值，然后在模拟过程中增大，如果再发生瞬变，最好人工调小时间步长，然后让它自动增大。
要解决问题的典型的不平衡力准则，可以通过观察只有弹性作用的初始阶段接近平衡时的不平衡力来决定。
很多情况下通过逐渐增大或减小时间步可以达到很好的性能。
（例如默认比值lmul = 1.01和umul = 0.90).一些情况下，最好避免对时间步进行连续调整，这可能产生“干扰”。
为了达到这个目的，时间步发生改变后应有一个“等待时间”（比如：100步），让系统暂停，这段时间中时间步不做进一步调整。
通常时间步开始与小的数值来适应短暂的过程比如开挖，然后随模拟进程增加。
如果增加一个短暂过程，则需要手工减小时间步然后在让其自动增加。
命令SET creep用来设置时间步以及时间步自动调整所需要的参数。
关键词见2.4.12.4蠕变模型的指令输入2.4.1 flac3d命令Config creep这个命令必须使用，以便开始蠕变分析。
History dt 时间步长Crtime 蠕变时间Model viscous 经典粘弹性模型Print creep 显示蠕变模型的参数信息Zone prop 显示分配给单元的材料属性Property bulk 弹性体积模量，KShear 弹性剪切模量，GViscosity 动力粘性常数，ηSet creep age t使用solve求解蠕变问题，t是蠕变时间的限值。
Dt tAuto onAuto offT就是蠕变时间步长，如果不进行赋值，则默认为mindt。
Auto on auto off是自动调整时间步长开关，默认为关；设置为开时，要通过关键字lfob，ufob，lmul，umul和latency来控制时间步长。
Lfob v不平衡力比率低于v值，则时间步长自动增大。
默认为310-Ufob v不平衡力比率超过v值，则时间步长自动减小。
默认为5.0*310-Lmul v不平衡力比率低于lfob，则时间步长自动乘以v值，lmul必须大于1，默认为1.01Umul v不平衡力比率超过ufob，则时间步长自动乘以v值，lmul必须小于1，默认为0.90Latency v延长阶段的时间步数Maxdt v最大蠕变时间步长，默认没有限制，但是设置的时候最好不要超过cr tmaxMindt v最小蠕变时间步长，默认为零On 蠕变时间开关，一旦确认蠕变分析（config creep）则默认为开OffTime t蠕变开始的时间，默认为tSolve age tconfig creep后，t为蠕变计算的时间限值FISH中可以用的变量crdt creep timestepcrtime creep time2.4.2 property经典粘弹性（Maxwell介质）— MODEL viscous(1) bulk 弹性体积模量, K(2) shear 弹性抗剪模量, G(3) viscosity 动态粘滞度, ηBurger’s Model — MODEL burger(1) bulk 弹性体积模量, K(2) kshear 开尔文抗剪模量, GK(3) kviscosity 开尔文粘性系数, ηK(4) mshear 麦斯韦尔剪切模量, GM(5) mviscosity 麦斯韦尔粘性系数, ηMBurger-Creep 粘塑性模型— MODEL cvisc(1) bulk 弹性体积模量, K(2) cohesion 内聚力, c(3) dilation 膨胀角, ψ(4) friction 内摩擦角, φ(5) kshear 开尔文抗剪模量, GK（6）kviscosity 开尔文粘性系数, ηK(7) mshear 麦斯韦尔剪切模量, GM（8）tension 极限抗拉强度, σt(9) mviscosity 麦斯韦尔粘性系数, ηM2.5蠕变的例子－FLAC3D处理二维问题（采用蠕变模型）newgen zone brick size 8 1 22 p0 0 0 -6 p1 4.2 0 -6 p2 0 1 -6 p3 0 0 5pause ；建立模型inter 1 face range x 4.1 4.3 y 0 1 z -6 5；加3个接触面inter 2 face range x 0 4.2 y 0 1 z -5.9 -6.1inter 3 face range x -.1 .1 y 0 1 z 5 -6pause；gen zone brick size 30 1 22 p0 4.2 0 -6 p1 50 0 -6 p2 4.2 1 -6 p3 4.2 0 5 ratio 1.05 1 1 pause；gen zone brick size 1 20 10 p0 0,0 0 p1 0 1 0 p2 -20 0 0 p3 0 0 5 ratio 1 1.1 1 pause；gen zone brick size 12 20 1 p0 0 0 0 p1 0 0 -6 p2 -20 0 0 p3 0 1 0 ratio 1 1.1 1 pause；gen zone brick size 30 20 1 p0 4.2 0 -6 p1 50 0 -6 p2 4.2 0 -55 p3 4.2 1 -6 ratio 1.05 1.1 1pause；gen zone brick size 8 20 1 p0 0 0 -6 p1 4.2 0 -6 p2 0 0 -55 p3 0 1 -6 ratio 1 1.1 1 pause；gen zone brick size 20 20 1 p0 0 0 -6 p1 0 0 -55 p2 -20 0 -6 p3 0 1 -6 ratio 1.1 1.1 1 pause；group soil ；定义组group exca1 range x -20 0 y 0 1 z 3 5group exca2 range x -20 0 y 0 1 z 0 3group mixingpiles range x 0 4.2 z -6 5；边界约束fix x range x -19.9 -20.1fix x range x 50.1 49.9fix yfix z range z -55.1 -54.9config creep；启动蠕变模块def burcalbkk=bme/(3*(1-2*b1))bk1=bke/(2*(1+b1))bk2=bme/(2*(1+b1))bvis1=vis1bvis2=vis2endset bme=1.35714set bke=3.25714set b1=0.35set vis1=1.35714set vis2=19.2857burcalmodel moh range group mixingpiles；定义本构model burger range group soil any group exca1 any group exca2 anyprop b b1 ksh bke msh bme kvis vis1 mvis vis2 range group soil any group exca1 any group exca2 anyprop s 1e9 b 2e9 c 2e10 range group mixingpiles；赋材料属性ini dens=1900 range group mixingpiles；定义密度ini dens=1800 range group soil any group exca1 any group exca2 anyinter 1 prop fric 40 kn 1e8 ks 1e8 ten 0；给接触面赋材料属性inter 2 prop fric 40 kn 1e8 ks 1e8 ten 0inter 3 prop fric 40 kn 1e8 ks 1e8 ten 0set grav 0 0 -10；重力加速度ini syy -4.5e4 grad 0 0 0.9e4；初始应力平衡ini sxx -4.5e4 grad 0 0 0.9e4ini szz -9e4 grad 0 0 1.8e4pause;初试平衡求解前solvepause;初试平衡求解结果ini xdis 0 ydis 0 zdis 0 ；初始位移为零set creep mindt=1e-2 maxdt=1set creep on；开启蠕变计算his unbalhis gp xdis 0. 0.5 5his crtimemodel null range group exca1；开挖第一步solve age 6；蠕变计算pause;开挖第一步求解结果model null range group exca2；开挖第二步solve age 12 ;开挖第2步求解结果2.6开挖后围岩随时间变化的例子model cvisc(或其它的蠕变本构模型) proper 各种蠕变参数hist id 1 crtime 记录蠕变时间hist id 2 gp xdisp 一点的三维坐标hist id 3 gp ydisp 一点的三维坐标hist id 4 gp zdisp 一点的三维坐标plot hist 2 vs 1plot hist 3 vs 1plot hist 4 vs 1。