1、第七章 隧道工程设计 中的有限元法7 7.1 概述 7.2 有限元法基础 7.3 隧道围岩弹塑性 有限元分析隧 道 工 程 SUI DAO GONG CHENG 7.4 工程实例分析 7.1 概述 岩土介质作为隧道工程的对象包含着多种随机因素岩土介质作为隧道工程的对象包含着多种随机因素(例如:非均匀性和各向异性,地质构造和结构面,应力-应变的非线性本构关系,初始地应力,地下水等等),正确掌握这些因素及其变化规律非常困难。常用的数值分析方法是有限元法有限元法、有限差分法有限差分法、边界元法边界元法、变分法变分法和加权余量法加权余量法。7.2 有限元法基础 有限元法的分析过程,概括起来可分为以下几
2、个步骤:1)1)将一个受力的连续体“离散化”。2)2)按静力等效原则将外力形成等效节点力。3)3)建立作用于在每个节点上力的平衡方程式。4)4)加入位移边界条件求解方程组,得到全部未知位移,进而求得各单元的应变和应力。7.3 隧道围岩弹塑性有限元法分析 7.3.17.3.1 概述 隧道围岩的受力属于空间问题,但因计算工作量大,数据处理费事,而将其简化为二维问题进行分析也能得到令人满意的结果,本节主要讨论二维隧道围岩弹塑性有限元分析。7.3.27.3.2 非线性问题的求解方法 采用数值方法分析结构时,将结构离散化后可以得到上一节中建立的代数方程组,0PK 当总刚度矩阵K中的元素Kij为常量时,为
3、线性方程组,它所代表的问题为线性问题。当Kij为变量时,例如Kij=f(dij),则为非线性方程组,它所描述的问题为非线性问题。材料非线性材料非线性:指的是当应力超过某一限值后,应力与应变的变化不成线性关系,但应变与位移的变化仍成线性关系但应变与位移的变化仍成线性关系。属于这种类型的问题称为材料非线性问题。几何非线性几何非线性:指的是当应变或应变速率超过某一限值后,应变与位移的变化不成线性关系,但应力与应变的变化仍成线性关系。属于这种类型的问题称为几何非线性问题。7.3.37.3.3 岩土材料的弹塑性本构关系 1)1)屈服条件和破坏条件,确定材料是否塑性屈服和破坏;2)2)硬化定律,指明屈服条
4、件由于塑性应变而发生的变化;3)3)流动法则,确定塑性应变的方向;4)4)加载和卸载准则,表明材料的工作状态。岩土材料的弹塑性应力应变关系即本构关系包括以下四个组成部分:1 1几种常用的屈服准则几种常用的屈服准则 目前常用于岩土材料的屈服准则有:摩尔摩尔-库仑库仑(Mohr-(Mohr-Coulomb)Coulomb)屈服准则屈服准则,德鲁克德鲁克-普拉格普拉格(Drucke-Prager(Drucke-Prager)屈服准则屈服准则,辛克维奇辛克维奇-潘迪潘迪(Zienkiewicz-Pande(Zienkiewicz-Pande)屈服准则屈服准则等。(1)(1)摩尔摩尔-库仑库仑(Mohr
5、-Coulomb)(Mohr-Coulomb)屈服准则屈服准则 库仑式 摩尔式0ctgf0cos2sin)()(3131cf 当剪切面上的剪应力与正应力之比达到最小时材料发生屈服与破坏。最大优点是不仅能反映岩土材料的拉伸与压缩的屈服与破坏强度不同(S-D效应)与对静水压力的敏感性,而且简单实用。因此在岩土力学和塑性理论中得到广泛应用。但该准则不能反映中间主应力对屈服和破坏的影响及单纯的静水压力可以引起岩土屈服的特性,而且屈服曲面有棱角,不便于塑性应变增量的计算,这就给数值计算带来了困难。(2)(2)德鲁克德鲁克-普拉格普拉格(Drucke-Prager(Drucke-Prager)屈服准则屈服
6、准则 考虑到静水压力可以引起岩土材料的屈服,德鲁克-普拉格屈服准则为0),(1221kIJJIf2sin33sin,2sin3cos3ck 德鲁克-普拉格屈服准则可以避免摩尔-库仑屈服准则屈面在棱角处引起的数值计算上的困难,但该准则对实际破坏条件逼近较差。(3)(3)辛克维奇辛克维奇-潘迪潘迪(Zienkiewicz-Pande(Zienkiewicz-Pande)屈服准则屈服准则 为了克服摩尔-库仑屈服准则的棱边和夹角,考虑到屈服与静水压力的非线性关系和中间主应力对强度的影响,辛克维奇-潘迪提出了辛克维奇-潘迪屈服准则,0)(12ngqkppf 2 2硬化法则硬化法则 硬化法则规定材料进入塑
7、性变形后的后继屈服函数(又称加载函数或加载曲面),0),(pF(7-28)对于理想弹塑性材料,因无硬化效应,显然后继屈服函数和初始屈服函数一致,即0)(),(FFp(7-29)对于硬化材料,通常有两种硬化法则:各向同性硬化法各向同性硬化法则则和运动硬化法则运动硬化法则。3 3流动法则流动法则 流动法则规定塑性应变增量的分量和应力分量以及应力增量分量之间的关系。米赛斯(Von.Mises)流动法则假设塑性应变增量可从塑性势导出,即Qpd 当QF时,这种特殊情况为关联塑性;否则称为非关联塑性。对于关联塑性情况,Fpd(7-31)4 4加载、卸载准则加载、卸载准则 若,继续塑性加载。0d,0FF 若
8、,由塑性按弹性卸载;0d,0FF 若 ,(1)对于理想弹塑性材料,塑性加载,(2)对于硬化材料,中性变载。0d,0FF 5 5弹塑性应力应变关系弹塑性应力应变关系 若令增量弹塑性关系表达式为ddepD 可以推导出QDFADFQDDDTTep 1 1分析步骤分析步骤 1)1)确定各种计算参数。2)2)确定各个施工阶段中围岩和结构的有限元网格。3)3)分阶段计算等效节点荷载,如释放节点荷载、自重 4)4)分阶段计算围岩和结构的应力及位移。5)5)叠加各施工阶段的计算值。6)6)评估围岩和结构的强度及稳定性。2 2施工阶段施工阶段 正确的隧道围岩分析方法应该考虑到分阶段施工的特点。有限元法能够方便地
9、模拟施工过程。7.3.47.3.4 隧道围岩弹塑性有限元静力分析的 方法和步骤 一般隧道施工可划分为几个阶段,如图7-5所示:图7-5隧道施工的各阶段 7.4 工程实例分析各阶段围岩强度发挥程度等值线a)b)c)图7-8各施工阶段边界位移图a)b)c)变形比例尺01020 mmmaxmin-1.09+1.66仰拱-12.29+1.02-10.63+0.13边墙-11.40+0.42-10.15+0.34拱部二衬-19.78+0.47-18.15+0.27-8.71+1.12边墙-13.68+0.83-13.82+0.77-8.17+0.59拱部初衬第三阶段第二阶段第一阶段衬砌类型 阶段maxminmaxmin表7-2各施工阶段衬砌最大、最小主应力(单位:Mpa)2 2锚杆与围岩相互作用的弹塑性有限元分析锚杆与围岩相互作用的弹塑性有限元分析图7-9 隧道ZK46+701.5断面锚杆支护布置图7-10 x向位移等值线图图7-11 y向位移等值线图图7-12 x向应力等值线图(单位:50kPa)图7-13 y向应力等值线图(单位:50kPa)