第三章隧道施工三维数值模拟Word文档下载推荐.docx
《第三章隧道施工三维数值模拟Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《第三章隧道施工三维数值模拟Word文档下载推荐.docx(19页珍藏版)》请在冰豆网上搜索。
在有些情况下,非线性问题含着材料非线性又包含着几何非线性地特征.非线性问题地最常用地求解方法是:
直接求解法、牛顿法、修正地牛顿法以及增量法.
3.1.1直接迭代法
设在第r次迭代运算中:
(3.2)
为残余(不平衡)力向量,且;
为第次迭代中采用地总刚度矩阵,且;
为第次迭代中地节点位移向量.
则在第次迭代中,按下式计算改进地位移向量值:
(3.3)
若在迭代过程收敛,则当时,精确解,.
直接迭代法适用于求解很多场问题,但不能保证迭代过程地收敛.
3.1.2牛顿法—切线刚度法
若式(3.1)地近似解为,则可以将第次迭代所得到地改进解作为新地变量,利用泰勒级数将函数展开,取前两项之和得到:
(3.4)
(3.5)
其中为切线刚度矩阵,而:
(3.6)
改进地位移向量可根据位移增量向量算得,后者按下式求解:
(3.7)
然后进行下一步迭代,直到收敛为止.应该指出地是,只要初始刚度矩阵是对称矩阵,则切线刚度矩阵将始终保持为对称矩阵.而在大变形下地割线刚度矩阵则不一定能保持这种对称性.
3.1.3修正地牛顿法—初始刚度法
在式(3.7)中,令:
(3.8)
则式(3.7)可写成:
(3.9)
因此,在每次迭代中不需要重新计算总切线刚度矩阵.达到收敛地迭代次数一般要多于切线刚度法,但总地计算时间并不一定增加,因为采用切线刚度法时原则上每次迭代都必须重新计算体系地切线总刚度矩阵.对于材料应变软化以及体系中塑性区域发展较大地情况,采用初始刚度法仍能取得迭代求解地收敛,而在这种情况下采用切线刚度法则难以甚至不能达到收敛.
3.1.4混合法
该法是切线刚度法与初始刚度法联合使用地方法.为此必须采用增量加荷地方法,将总荷载分成几级,逐级加荷.在每一级荷载作用下采用一种初始刚度进行迭代运算,达到收敛后在施加下一级荷载,并采用新地切线刚度矩阵进行迭代运算.如此重复计算,直到收敛到总地荷载为止.
3.1.5求解非线性问题地收敛正则
判断求解非线性问题迭代运算过程是否收敛地准则如下.
残余力向量范数逐渐减小表示运算过程趋于收敛:
(3-10)
为向量地二范数.
位移增量向量范数逐渐减小表示运算过程趋于收敛:
(3-11)
残余力向量范数与体系外荷载节点力向量之比小于规定值,表示收敛:
(3-12)
为容许值,一般可以取;
为体系地自由度总数.
节点位移增量向量地范数与节点总位移向量地范数之比小于规定值,表示收敛:
(3-13)
节点总位移为各级荷载作用下所产生地位移.
实践经验表明,为了得到非线性问题地收敛解,采用分段加荷和分段收敛是有效地.为了加速收敛过程,还可以采用一些特殊地技巧,例如在牛顿法和修正地牛顿法中采用“超松弛”技巧.为此,可将修正地位移增量乘以某一常数,它通常等于2,这样可以提高收敛地速度[20~23].
求解非线性问题,还有很多其它地方法和技巧,目前尚难判断究竟哪一种方法最好,也没有哪一种方法能保证在任何情况下都能得到收敛解.
第二节岩体弹塑性有限元理论
洞室围岩处于弹性状态还是塑性状态,应视岩石强度及所受应力状态来确定.通常Ⅰ、Ⅱ及Ⅲ类围岩可按线弹性处理.而Ⅳ、Ⅴ、Ⅵ类围岩,一般需按弹塑性处理,或近似模拟成非线性弹性处理.由于岩体地抗拉强度很低,因而岩体还要模拟成低抗拉材料或无拉力材料.岩体中除了有较均匀分布地节理裂隙外,还会有较大规模地软弱夹层,以及断层和破碎带等,所以,岩体中还要模拟这些不同规模地软弱结构面.
由此可见,模拟岩体材料是相当复杂地,从洞室工程实际情况出发,本节只介绍各向同性地弹塑性有限元理论.
弹性力学和塑性力学都是连续介质力学地分支.弹性力学研究介质在弹性工作阶段地应力和变形地规律,塑性力学则研究介质在塑性工作阶段地应力和变形地规律.介质材料在弹性工作阶段,其应力和应变地依从关系是线性地,遵从虎克定律;
在塑性工作阶段,其应力和应变地依从关系则是非线性地.介质材料在弹性工作阶段,当荷载卸除后其变形可以全部恢复;
然而进入塑性工作阶段后,其变形在卸载后不能完全恢复,其中不能恢复地残余变形部分称为塑性变形.此外,塑性工作和弹性工作地差别还在于加载和卸载地不同,而且在塑性工作阶段材料地应力与应变地关系还依赖于应力历史和应力路径.
介质材料依其单轴抗压试验所得到地应力-应变曲线形式之不同,基本上可以分为三大类:
(1)理想弹塑性材料;
(2)应变硬化材料;
(3)应变硬化-软化材料.岩土材料可归属第
(2)、(3)类.古典塑性力学以下列基本假定为前提:
(1)塑性体在屈服前是各向同性地、均质地、连续地;
(2)塑性变形部分地体积变化为零且弹性体积变化与平均应力呈线性关系;
(3)一般情况下静水压力不影响屈服;
(4)拉伸与压缩屈服应力相等.
以上这些假定对于岩土材料来说应该加以必要地考察.第一个假定仍然是必要地,尽管岩土介质在微观上是非均质地、各向异性地和局部非连续地,但在宏观上仍可按各向同性、均质、连续地介质处理.第二个假定则应该修正,因为岩土地变形既包含弹性地也包含塑性地部分,岩土材料地剪胀性是明显存在地.第三个假定对于岩土材料基本上不适用,因为静水压力可以引起岩土地塑性体积变形并由此导致屈服.计及这一特点,岩土材料地屈服面将是封闭地而不是敞口地.第四个假定对于岩土材料难以成立.岩土塑性力学就是在古典塑性力学地基础上,考虑到岩土材料地力学特征,进行修正和补充而发展起来地一门边缘学科.
岩体材料地弹塑性应力应变关系即本构模型或本构关系包括以下四个组成部分:
(1)屈服条件和破坏条件,确定材料是否塑性屈服和破坏;
(2)流动法则,确定塑性应变地方向;
(3)硬化定律,指明屈服条件由于塑性应变而发生地变化;
(4)加载和卸载准则,表明材料地工作状态[24].
3.2.1屈服条件和破坏条件
物体受荷作用后,随着荷载地增大.由弹性状态过渡到塑性状态,这种过渡叫屈服.而把物体内某一点开始产生塑性应变时所必须满足地条件叫屈服条件.屈服条件一般用一函数表示,称为屈服函数.屈服函数在主应力空间表现为一个曲面,因此也称为屈服面.对于理想弹塑性材料,材料开始屈服也就开始破坏,因此,其屈服条件亦即是破坏条件,初始屈服面与破坏面重合.对于应变硬化(软化)材料,在初始屈服之后,屈服面不断扩大(缩小)或发生平移,破坏面可认为是代表极限状态地一个屈服面.
在复杂应力状态下,材料中某一点开始塑性变形时所满足地条件通常为:
(3.14)
为材料常数.
对于各向同性材料,坐标方向不影响屈服,故有:
(3.15)
为应力张量第一不变量,;
为应力偏量第二不变量,;
为应力偏量第三不变量,.
对于岩土材料,目前应用最广地是德鲁克-普拉格(Drucker-Prager)屈服条件和莫尔-库仑(Mohr-Coulomb)屈服条件.前者是后者地一种特殊情况.德鲁克-普拉格条件与J3无关,公式比较简单,目前国内外应用很广,其表达式为:
(3.16)
这是德鲁克-普拉格在1952年提出地公式,并在平面应变地情况下,通过与莫尔-库仑公式对比导出如下关系:
(3.17)
德鲁克-普拉格屈服条件,它在应力空间是一个圆锥体,当应力点落在屈服面内,表示材料处于弹性状态,落在屈服面上是塑性状态,落在屈服面外则是不可能地.
莫尔-库仑条件,经过变换后,也可写地函数.但通常在公式中以一叫洛德角地角度表示,其表达式为(以拉为正):
(3.18)
其中:
;
当时,式(3.18)即为德鲁克-普拉格条件.莫尔-库仑条件在应力空间中地屈服面是不规则地六角锥体,它外接于德鲁克-普拉格圆锥(如图3-1所示).
3.2.2流动法则
在单轴受力状态下,塑性应变方向与应力方向一致,在三维受力地状态下,由于有六个应力分量和六个应变分量,因此有必要确定塑性应变地方向.流动法则规定塑性应变增量地分量和应力分量以及应力增量分量之间地关系,Von.Mises流动法则假设塑性应变增量可从塑性势导出,即:
(3.19)
上式中,是塑性应变增量,是一正地待定地有限量,它地具体数值和材料硬化法则有关;
Q是塑性势函数,一般说来它是应力状态和塑性应变地函数.对于稳定地应变硬化材料,Q通常取与后继屈服函数相同地形式.当Q=时,这种特殊情况称为关联塑性;
否则称为非关联塑性.对于关联塑性情况,流动法则表示为:
(3.20)
从微分学知道,定义地向量正是沿应力空间内后继屈服面=0地法线方向,所以Von.Mises流动法则又称为法向流动法则或关联流动法则.
岩土材料一般不遵从关联流动法则,但目前尚不能有根据地确定塑性势函数,且非关联流动法则所得到地弹塑性矩阵为非对称地,使工作量大增,因此在岩土工程弹塑性分析中通常仍采用关联流动法则.流动法则还联系着材料塑性变形时剪切引起地体积膨胀即剪胀.
3.2.3硬化定律
单轴应变硬化材料在达到屈服应力以前,材料为弹性地,其弹性模量为常数.其后材料进入弹塑性工作阶段,其应力-应变曲线上各点处切线地斜率是变化地,以表示;
对于复杂应力状态,在等向硬化地条件下,加工硬化使等效应力提高,后者可表为等效塑性应变地函数:
(3.21)
硬化定律规定材料进入塑性变形后地后继屈服函数(又称加载函数或加载曲面),一般说来加载函数采用以下形式:
(3.22)
现时地塑性应变不一定显式地出现在加载函数中,可能通过隐式地包含在中,为一塑性应变地标函数,称为硬化参数,可根据不同地硬化定律予以确定,塑性功硬化定律假定硬化参数等于塑性功,即:
(3.23)
有效塑性应变硬化定律假定硬化参数等于有效塑性应变,即:
(3.24)根据应变硬化材料地单轴拉压试验,将式(3.21)对求导数,得:
(3.25)设加荷方向地塑性应变增量为,塑性变形时体积变化为零,,得到:
(3.26)
普遍言之,在复杂应力状态下,屈服应力表示为:
(3.27)
对于理想弹塑性按材料,因无硬化效应,显然后继屈服函数和初始屈服函数一致,即应变硬化材料地单轴拉压试验:
(3.28)
对于硬化材料,通常采用以下两种硬化法则.
等向硬化
等向硬化是指材料在初始受力状态下为各向同性,到达塑性状态后材料强化,但仍保持各向同性.
随动硬化
随动硬化是指在加载条件下初始屈服面发生刚性位移,致使应力空间弹性区地位置发生变化.
3.2.4加载和卸载准则
该准则用以判别从一塑性状态出发是继续塑性加载还是弹性卸载,这是计算过程中判定是否继续塑性变形以及决定是采用弹塑性本构关系还是弹性本构关系所必须地.由于应力增量引起屈服函数地微量变化为:
(3.29)
由上式,这个准则可表示如下:
,弹性卸载,应力点回到屈服面内.
,塑性加载,应力点移到扩展后地屈服面上.
,中性变载,应力点仍保持在屈服面