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 为剪切模量。

对于应力应变的体积分量，假设体积变化不受流变影响。

kk kk isoK εσσ∆+=31 （2.5）K 为体积模量，最终的应力张量就表示为偏量和球量之和：ij isodij 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。