流变、蠕变问题_土木相关.ppt

* ITASCA(武汉)咨询有限公司 报告人: 朱永生 2007.04.28 依泰斯卡(武汉) 内 容 概况 粘弹性体变形特征 粘弹性体流变元件 常用基本流变模型 软件内置的蠕变本构模型 用FLAC/FLAC3D解决蠕变问题 蠕变分析输入介绍 弹性理论：假定物体的形变与所受的力成正比，形变值的大小与加荷的历程无关，不会随时间变化； 塑性理论：应力-应变关系为非线性，若荷载本身不便，也是假设所得的新形变在以后维持不变，同样与时间无关； 流变理论： 荷载不变，形变值并非永远是常数； 荷载逐渐增大，形变取决于加荷速度 应力与应变关系的本构模型中包含时间。 概况 力学响应对时间的依赖型 岩石、岩体，特别是结构面(层面、节理、层理)在较高应力水平条件下的变形性态均与时间有关，可以将其看作是粘弹性体。 粘弹性体的流变特征，主要表现为蠕变、松弛、流动与长期强度等几个方面： (1) 在定常荷载作用下，应变(ε)随时间( t )而逐渐增长的现象，称为蠕变特征； (2) 当应变一定时，应力(σ)随时间 ( t )而逐渐衰减的现象，称为松弛现象； (3) 应力σ随着应变速率(d ε/d t) 逐渐增长的现象，称为流动特性； (4) 在一定的时间内，强度(τ)与时间( t )的关系，称为长期强度。 ε t I II III A B C D 加载曲线 卸载曲线 蠕变曲线 粘弹性体变形特征 流变元件 常用基本流变模型 uH uN F F 松弛模型(Maxwell体） M=H-N ε t t σ UH= UN F F 延迟模型(Kelvin体） M=H||N σ t ε t 理想弹塑性模型 EP=H-St.V F F F F N St.V 粘塑性模型 VP=N||St.V 表. 各种岩石曾采用的蠕变模型 软件内置蠕变本构模型 用FLAC/FLAC3D解决蠕变问题 简述 蠕变分析时间步长 蠕变时间步的自动调整 蠕变和动力计算 流变模型和FLAC/FLAC3D其他模型最大的不同在于模拟过程中时间概念的不同，对于蠕变，求解时间和时间步代表着真实的时间，而一般模型的静力分析中，时间步是一个人为数量，仅仅作为计算从迭代到稳态的一种手段来使用。 简述 蠕变分析时间步长 对于蠕变等时间依赖性问题，程序允许用户自定义一个时间步长，这个时间步长的默认值为零，那么材料对于粘弹性模型表现为线弹性，对于粘塑性模型表现为弹塑性。（命令set creep off也可以用来停止蠕变计算。） 建立初始平衡: 蠕变公式中包含时间，所以计算中时间步长对程序响应有影响，以上方法可以用来在系统达到平衡后再开始新的蠕变计算。 用户可以改变人为改变蠕变时间步的步长，但步长不能任意给定，必须满足差分方法收敛的步长要求。 通常，蠕变过程由偏应力状态控制，从数值计算的精度来讲，最大蠕变时间步长可以表示成粘性常数和剪切模量的比值： Power Law: Cvisc(Burger-Creep viscoplastic model): 力学过程存在压缩行为(Crushed-Salt Model，压硬和减缩)，蠕变步长由粘性常数与体积模量的比值确定： WIPP: 建议: 建议利用程序作蠕变分析开始时所采用的蠕变时间步，比根据上式算得的临界时间tmax小两到三个数量级。通过调用SET crdt auto / SET creep dt auto on ，可以利用自动时间步自动调整。作为一项规则，时间步的最大值maxdt （SET maxdt / SET creep maxdt ）不能超过tmax。 用来计算tmax的应力σ大小，可由蠕变开始之前的初始应力状态决定。同样，σ作为von Mises不变量，可以用FISH函数计算。 蠕变时间步的自动调整 蠕变时间步长可由用户给定为常值： FLAC: set crdt t FLAC3D: set creep dt t 蠕变时间步长由程序自动调节： FLAC: set crdt auto FLAC3D: set creep dt auto on/off 如果时间步长自动变化，那么当最大不平衡力超过某一阀值时，它就会减小；当最大不平衡力小于某一水平时它就会增大。 阀值的定义： FLAC：通过最大不平衡力来控制； ① Funbal≤fobl，时间步乘子lmul(default=2.0); ② Funbal≥fobu，时间步乘子umul(default=0.5)。 FLAC3D：通过最大不平衡力与平