flac3D蠕变基础知识03.doc

PAGE/NUMPAGES

flac３Ｄ蠕变基础知识

蠕?变?模?型

将fｌaｃ3d的蠕变分析ｏpｔion进行了简单的翻译，目的是为了搞清楚蠕变过程中系统时间是如何跟真实时间对应的。

1.简介

Ｆｌac3d可以模拟材料的蠕变特性,即时间依赖性,flａｃ３d2。1提供6种蠕变模型：

1．?????经典粘弹型模型modeｌｖｉscｏus

2.?????mｏｄｅｌburgeｒ

3。?????modelｐｏwｅr

4。?????modelwipp

5．?????modelｃviｓc

6.?????pｏwe蠕变模型结合M—C模型产生cpow蠕变模型（modeｌcｐｏｗ)

７.?????然后WIPP蠕变模型结合D—P模型产生Pｗｉpp蠕变模型（modelpwipp）；

8???????modelcwiｐp

以上模型越往下越复杂，第一个模型使用经典的ｍaxwｅll蠕变公式,第二个模型使用经典的bｕrｇeｒ蠕变公式，第三个模型主要用于采矿及地下工程,第四个模型一般用于核废料地下隔离的热力学分析，第五个模型是第二个模型的Ｍ-C扩展，第六个模型是第三个模型的M-C扩展，第七个模型是第四个模型的D—P扩展，第八个模型也是第四个模型的一种变化形式，只是包含了压硬和剪缩行为。

２.fｌac３d解流变问题

２．１简介

流变模型和flac３ｄ其他模型最大的不同在于模拟过程中时间概念的不同，对于蠕变，求解时间和时间步代表着真实的时间,而一般模型的静力分析中，时间步是一个人为数量,仅仅作为计算从迭代到稳态的一种手段来使用。

2．2fｌａc3d的蠕变时间步长

对于蠕变等时间依赖性问题,flac3d容许用户自定义一个时间步长,这个时间步长的默认值为零,那么材料对于粘弹性模型表现为线弹性，对于粘塑性模型表现为弹塑性。（命令ｓetcｒｅepoff也可以用来停止蠕变计算。)这可以用来在系统达到平衡后再开始新的蠕变计算.蠕变公式中包含时间，所以计算中时间步长对程序响应有影响。

虽然用户可以对时间步进行设置，但并不是任意的.

蠕变过程由偏应力状态控制，从数值计算的精度来讲，最大蠕变时间步长可以表示成材料粘性常数和剪切模量的比值：?????????????????

Ｆｏrｔhepowｅｒｌaｗ－——－-—－省略.??????FoｒtheWIPＰlaw—－-—-—-—-—-省略

Forｔhecviscmｏdel，上面方程应该写成:??????tmａx=min（ηK/GＫ，ηM/GM)

上标K和M分别代表Ｋelvin和Maxwelｌ。

???蠕变压缩的时间限制包括系统体积反应,并且估计为粘性和体积模量的比值.粘性可以表示为σ和体积蠕变压缩速率的比值。

???建议利用FLAC3Ｄ作蠕变分析开始时所采用的蠕变时间步，比根据上式算得的时间tmax小两到三个数量级.通过调用SETcrｅｅpdｔautooｎ,可以利用自动时间步自动调整。作为一项规则,时间步的最大值（SEＴcreeｐmaxdt)不能超过ｔｍax.????用来计算tmax的应力σ大小，可由蠕变开始之前的初始应力状态决定。同样,σ作为voｎMises不变量,可以用FIＳH函数计算。

涉及体积变化响应的蠕变分析,其最大时间步长可以表示成材料粘性常数和体积模量的比值,这里粘性常数就是平均应力和蠕变体应变率的比值.

一般flａc３d推荐使用的初始蠕变时间步长比最大时间步长（由上述公式计算得到的)约小２到3个数量级。如果使用seｔｃreepdｔautoon命令，那么程序将自动调整蠕变的时间步长，同样应当记住通过命令（sｅtcreepmaｘdt)设置的最大蠕变时间步不能超过.

2．3自动调整蠕变时间步长

用户可以设置蠕变时间步为一个常数值,也可以使用setcrｅｅpdtaｕｔoｏｎ命令自动调节。如果时间步长自动变化,那么当最大不平衡力超过某一阀值时,它就会减小;当最大不平衡力小于某一水平时它就会增大。系统将该阀值定义为最大不平衡力和平均节点力的比值。

对于即将求解的问题，先只考虑弹性效应，通过观察接近初始平衡状态时的不平衡力就可以大概知道整个问题的不平衡力大小。

在有些算例中，又会尽量避免对时间步长的连续调整，这样的话，在某一次时间步长调整之后,可以定义一个“延迟阶段”（比如100步），在这一“延迟阶段内不再发生进一步的时间步长调整.一般来说，为了适应开挖之类的瞬变，时间步长一开始是一个较小值,然后在模拟过程中增大，如果再发生瞬变，最好人工调小时间步长,然后让它自动增大。????要解决问题的典型的不平衡力准则,可以通过观察只有弹性作用的初始阶段接近平衡时的不平衡力来决定。很多情况下通过逐渐增大或减