下载文档原格式
一种蠕行控制方法叫什么引言蠕行（creep）是材料学中常见的现象，指的是材料在长时间受到持续作用力时发生的无反弹性变形。
蠕行现象在许多领域都存在，如工程结构设计、材料加工和地质工程等。
为了控制蠕行现象，人们提出了许多方法。
本文将介绍一种“蠕行控制方法”。
蠕行控制方法的原理目前通常使用的蠕行控制方法主要分为两大类：应力降低法和材料改性法。
- 应力降低法：通过减小材料的受力强度，使其承受更小的应力从而降低蠕变速率。
例如，在工程结构设计中，可以采用增大材料截面尺寸、增加结构支撑等方式来降低结构的应力水平，从而减缓蠕变速率。
- 材料改性法：通过改变材料的组织结构，调控材料的形变行为，从而控制蠕行现象。
例如，通过合金化、陶瓷增强、纤维增强等手段来改变材料的微结构，增强材料的抗蠕变性能。
然而，以上方法都有其局限性。
应力降低法需要在设计初期耗费大量的时间和精力，并且难以在运营期进行调整；材料改性法则需要以提高材料强度为代价，无法得到完美的控制效果。
因此，人们需要一种更加高效和实用的蠕行控制方法。
一种新的蠕行控制方法：多物理场耦合法近年来，随着科学技术的发展和理论模拟能力的提高，人们提出了一种新的蠕行控制方法，即多物理场耦合法。
这种方法基于多物理场的相互作用，通过控制不同物理场的耦合行为来实现对蠕行现象的控制。
多物理场指的是同时存在于同一介质或系统中的多个物理场，例如电场、热场、磁场、应力场等。
而耦合行为则是指这些物理场之间的相互影响和相互作用。
通过调控不同物理场的大小、分布和时间演化等参数，可以改变材料内部的微观状态，从而影响材料的蠕变行为。
多物理场耦合法的具体实施方式可以通过数值模拟、实验验证和实际应用等不同途径来进行。
其中，数值模拟是其中的重要手段之一。
通过建立适当的物理模型和数值算法，可以模拟不同物理场的耦合行为，并通过改变参数来控制蠕行现象。
实验验证可以通过实验室试验和实际工程应用来进行，通过对不同物理场的加载和监测，验证多物理场耦合法对蠕行的控制效果。
材料的蠕变机制了解材料在高温下的变形机制材料的蠕变机制：了解材料在高温下的变形机制材料的蠕变（Creep）是指在高温和应力作用下，材料会逐渐发生形变和软化的现象。
蠕变机制的了解对于材料在高温环境下的设计和应用非常重要。
本文将探讨材料在高温下的变形机制，以及蠕变现象对材料性能的影响。
1. 引言随着工程领域对材料性能要求的提高，对于材料在极端条件下的变形行为的研究变得愈发重要。
高温下材料的蠕变现象是一项关键的考察对象，本文旨在深入探讨材料在高温环境下的变形机制，并分析其对工程应用的意义。
2. 材料在高温下的变形机制2.1 晶体滑移晶体滑移是一种经典的材料变形机制，在高温条件下尤其突出。
晶体原子通过在晶格面上滑移来改变位置，从而导致材料的塑性变形。
高温会增加晶体内的原子动能，使得晶体滑移更容易发生。
2.2 相互扩散在高温下，材料中的原子会因为动能提高而显示出相互扩散的特性。
相互扩散会导致晶界和晶体内部的松弛，进而引起材料的蠕变变形。
2.3 晶体再结晶晶体在高温和应力作用下可能发生再结晶，即原本晶体中的晶粒重新排列和重组。
这种再结晶过程也是材料蠕变的一个重要机制。
3. 蠕变对材料性能的影响3.1 强度和韧性的降低材料蠕变使得材料在高温下的强度和韧性均降低。
原子的扩散会导致材料晶界的松弛，晶体结构的破坏以及晶粒的再结晶，这些因素会使材料的强度降低。
3.2 变形速率蠕变行为与应力的大小、温度的高低以及时间的长短密切相关。
蠕变变形速率随温度的升高而增加，随应力的增加而加快。
这对一些高温工程应用中材料的耐久性和设计有着重要影响。
4. 变形机制的研究与应用4.1 实验方法通过高温实验设备对材料进行蠕变试验，可以模拟材料在高温下的变形行为。
通过测量变形曲线、变形速率和变形温度等参数，可以获取材料的蠕变特性。
4.2 材料改进与设计通过深入研究材料的蠕变机制，可以针对不同的应用需求改进和设计材料。
例如，在航空航天领域，对发动机叶片等高温结构部件材料的蠕变机制进行研究，以提高其抗蠕变性能。
土层蠕变的计算公式土层蠕变（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模型的弹簧常数。
这些土层蠕变的计算公式是经过实验验证和统计分析得出的经验公式，可以较为准确地预测土壤蠕变变形的发展趋势。
各个论坛均有关于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相当。
蠕变模型将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）：内耗是指聚合物材料在受力或形变时，由于分子内部摩擦和相互作用而产生的能量损耗。
内耗可以导致材料的能量耗散和温升。
聚合物材料的内耗通常与材料的分子结构、聚合度和温度等因素有关。
内耗在聚合物的动态力学性能和阻尼特性中起着重要作用。
这些现象在聚合物工程和材料科学中具有重要的应用。
研究和了解聚合物的蠕变、应力松弛、滞后和内耗行为对于设计和开发具有特定力学性能和可靠性的聚合物制品非常重要。
蠕变应变（Creep Strain）是指物质在长时间（很长的时间段）下持续受力时所发生的变形现象。
蠕变应变主要发生在高温、高应力条件下，例如金属、陶瓷、塑料等材料在高温环境下承受恒定或逐渐增加的应力。
在蠕变过程中，材料会逐渐发生塑性变形，而且这种变形具有时间依赖性。
总应变（Total Strain）是指材料在外部作用下发生的总体变形量。
总应变可以分解为弹性应变和塑性应变两个部分。
●弹性应变（Elastic Strain）是指材料在外部加载后发生的可恢复的变形。
当外部加
载移除时，材料会完全回复到原始状态，没有永久性变形。
弹性应变与斯托克斯定
律相符。
●塑性应变（Plastic Strain）是指材料在超过弹性极限的应力条件下发生的不可恢复
性变形。
塑性应变会导致材料的内部结构发生改变，无法通过去除应力来完全恢复
原始形状。
在材料的应力-应变曲线上，弹性应变位于材料的线性阶段，而塑性应变则位于非线性阶段。
总应变等于弹性应变加上塑性应变。
总结起来，蠕变应变是指材料在长时间和高应力条件下的变形现象，而总应变则是材料在外部作用下的总体变形量，包括弹性应变和塑性应变。
flac3D蠕变基础知识分类：岩土蠕变 | 标签：FLAC3D creep2009-06-09 18:37 阅读(1422)评论(0)收集了一些FLAC3D的蠕变基础知识，希望对有需要的人起到帮助作用，欢迎下载！蠕变模型将flac3d的蠕变分析option进行了简单的翻译，目的是为了搞清楚蠕变过程中系统时间是如何跟真实时间对应的。
1. 简介Flac3d可以模拟材料的蠕变特性，即时间依赖性，flac3d2.1提供6种蠕变模型：1. 经典粘弹型模型model viscous2. model burger3. model power4. model wipp5. 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. flac3d解流变问题2.1简介流变模型和flac3d其他模型最大的不同在于模拟过程中时间概念的不同，对于蠕变，求解时间和时间步代表着真实的时间，而一般模型的静力分析中，时间步是一个人为数量，仅仅作为计算从迭代到稳态的一种手段来使用。
2.2 flac3d的蠕变时间步长对于蠕变等时间依赖性问题，flac3d容许用户自定义一个时间步长，这个时间步长的默认值为零，那么材料对于粘弹性模型表现为线弹性，对于粘塑性模型表现为弹塑性。
（命令set creep off也可以用来停止蠕变计算。
蠕变模型将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。
