flac3D蠕变基础知识 蠕 变 模 型 将flac3d的蠕变分析option进行了简单的翻译，目的是为了搞清楚蠕变过程中系统时间是如何跟真实时间对应的。
1. 简介 Flac3d可以模拟材料的蠕变特性，即时间依赖性，flac3d2.1提供6种蠕变模型： 1. 经典粘弹型模型 model viscous 2. model burger 3. model power 4. model wipp 5. model cvisc 6. 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也可以用来停止蠕变计算。）这可以用来在系统达到平衡后再开始新的蠕变计算。蠕变公式中包含时间，所以计算中时间步长对程序响应有影响。 虽然用户可以对时间步进行设置，但并不是任意的。 蠕变过程由偏应力状态控制，从数值计算的精度来讲，最大蠕变时间步长可以表示成材料粘性常数和剪切模量的比值： For the power law ----------省略。 For the WIPP law -----------省略 For the cvisc model, 上面方程应该写成： tmax = min ( ηK/GK，ηM/GM) 上标K和M分别代表Kelvin和Maxwell。 蠕变压缩的时间限制包括系统体积反应，并且估计为粘性和体积模量的比值。粘性可以表示为σ和体积蠕变压缩速率的比值。 建议利用FLAC3D作蠕变分析开始时所采用的蠕变时间步，比根据上式算得的时间tmax小两到三个数量级。通过调用SET creep dt auto on ，可以利用自动时间步自动调整。作为一项规则，时间步的最大值（SET creep maxdt ）不能超过tmax。 用来计算tmax的应力σ大小，可由蠕变开始之前的初始应力状态决定。同样，σ作为von Mises不变量，可以用FISH函数计算。 涉及体积变化响应的蠕变分析，其最大时间步长可以表示成材料粘性常数和体积模量的比值，这里粘性常数就是平均应力 和蠕变体应变率 的比值。 一般flac3d推荐使用的初始蠕变时间步长比最大时间步长（由上述公式计算得到的）约小2到3个数量级。如果使用set creep dt auto on命令，那么程序将自动调整蠕变的时间步长，同样应当记住通过命令（set creep maxdt）设置的最大蠕变时间步不能超过 。
2.3自动调整蠕变时间步长 用户可以设置蠕变时间步为一个常数值，也可以使用set creep dt auto on命令自动调节。如果时间步长自动变化，那么当最大不平衡力超过某一阀值时，它就会减小；当最大不平衡力小于某一水平时它就会增大。系统将该阀值定义为最大不平衡力和平均节点力的比值。 对于即将求解的问题，先只考虑弹性效应，通过观察接近初始平衡状态时的不平衡力就可以大概知道整个问题的不平衡力大小。 在有些算例中，又会尽量避免对时间步长的连续调整，这样的话，在某一次时间步长调整之后，可以定义一个“延迟阶段”（比如100步），在这一“延迟阶段”内不再发生进一步的时间步长调整。一般来说，为了适应开挖之类的瞬变，时间步长一开始是一个较小值，然后在模拟过程中增大，如果再发生瞬变，最好人工调小时间步长，然后让它自动增大。 要解决问题的典型的不平衡力准则，可以通过观察只有弹性作用的初始阶段接近平衡时的不平衡力来决定。很多情况下通过逐渐增大或减小时间步可以达到很好的性能。（例如默认比值lmul = 1.01和umul = 0.90).一些情况下，最好避免对时间步进行连续调整，这可能产生“干扰”。为了达到这个目的，时间步发生改变后应有一个“等待时间”（比如：100步），让系统暂停，这段时间中时间步不做进一步调整。通常时间步开始与小的数值来适应短暂的过程比如开挖，然后随模拟进程增加。如果增加一个短暂过程，则需要手工减小时间步然后在让其自动增加。命令SET creep用来设置时间步以及时间步自动调整所需要的参数。
2.4蠕变模型的指令输入 2.4.1 flac3d命令 Config creep 这个命令必须使用，以便开始蠕变分析。 History dt 时间步长 Crtime 蠕变时间 Model viscous 经典粘弹性模型 Print creep 显示蠕变模型的参数信息 Zone prop 显示分配给单元的材料属性 Property bulk（弹性体积模量，K）Shear（弹性剪切模量，G）Viscosity（动力粘性常数， ） Set creep age t 使用solve求解蠕变问题，t是蠕变时间的限值。 Set creep dt ** Set creep mindt=** maxdt=** Set creep Dt Auto on/ Auto off T就是蠕变时间步长，如果不进行赋值，则默认为mindt。Auto on auto off是自动 调整时间步长开关，默认为关；设置为开时，要通过关键字lfob，ufob，lmul，umul和latency来控制时间步长。 Lfob v 不平衡力比率低于v值，则时间步长自动增大。默认为 Ufob v 不平衡力比率超过v值，则时间步长自动减小。默认为5.0* Lmul v 不平衡力比率低于lfob，则时间步长自动乘以v值，lmul必须大于1，默认为 1.01 Umul v 不平衡力比率超过ufob，则时间步长自动乘以v值，lmul必须小于1，默认为 0.90 Latency v 延长阶段的时间步数 Maxdt v 最大蠕变时间步长，默认没有限制，但是设置的时候最好不要超过 Mindt v 最小蠕变时间步长，默认为零 On 蠕变时间开关，一旦确认蠕变分析（config creep）则默认为开 Off Time t 蠕变开始的时间，默认为t Solve age t config creep后，t为蠕变计算的时间限值 FISH中可以用的变量 crdt --------creep timestep crtime-------- creep time
2.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 麦斯韦尔粘性系数, ηM Burger-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 麦斯韦尔粘性系数, ηM
3. 蠕变的例子－FLAC3D处理二维问题（采用蠕变模型） new gen zone brick size 8 1 22 p0 0 0 -6 p1 4.2 0 -6 p2 0 1 -6 p3 0 0 5 pause ；建立模型 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.1 inter 3 face range x -.1 .1 y 0 1 z 5 -6 pause； 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 1 pause； 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 5 group exca2 range x -20 0 y 0 1 z 0 3 group mixingpiles range x 0 4.2 z -6 5 ；边界约束 fix x range x -19.9 -20.1 fix x range x 50.1 49.9 fix y fix z range z -55.1 -54.9 ============================================= config creep；启动蠕变模块 def burcal bkk=bme/(3*(1-2*b1)) ‘bkk——k(体积模量);bme——E（弹性模量）；b1——泊松比； bk1=bke/(2*(1+b1)) ‘bke—— bk2=bme/(2*(1+b1)) ‘bk2——G（剪切模量）；bme——E（弹性模量）；b1——泊松比；kshear 开尔文抗剪模量, GK bvis1=vis1 ‘vis——kviscosity 开尔文粘性系数, ηK bvis2=vis2 end set bme=1.35714 set bke=3.25714 set b1=0.35 set vis1=1.35714 set vis2=19.2857 burcal model moh range group mixingpiles；定义本构 model burger range group soil any group exca1 any group exca2 any prop b b1 ksh bke msh bme kvis vis1 mvis vis2 range group soil any group exca1 any group exca2 any prop s 1e9 b 2e9 c 2e10 range group mixingpiles；赋材料属性