第三章 隧道施工三维数值模拟

《第三章 隧道施工三维数值模拟》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、兰州交通大学硕士论文第三章隧道施工三维数值模拟第三章隧道施工三维数值模拟本章主要介绍非线性有限元方程组的解法，岩体的弹塑性理论，乌鞘岭隧道F7断层施工模型和结果分析。第一节非线性有限元方程组的基本解法采用数值方法分析结构时，将结构离散化后可以得到如下的代数方程组：(3.1)式中：为结构的总刚度矩阵；为未知数（位移等）向量；为外荷载向量（）。当总刚度矩阵中的元素为常量时，式（3.1）为线性方程组，它所代表的问题为线性问题。当为变量时，例如，则式（3.1）为非线性方程组，它所描述的问题为非线性问题。材料非线性指的是当应力超过某一极限值后，应力与应变的变化不成

2、线性关系，但应变与位移的变化仍为线性关系。属于这种类型的问题称为材料非线性问题。几何非线性指的是当应变或应变速率超过某一极限值后，应变与位移的变化不成线性关系，但应力与应变的变化仍成线性关系。属于这种类型的问题称为几何非线性问题。在有些情况下，非线性问题含着材料非线性又包含着几何非线性的特征。非线性问题的最常用的求解方法是：直接求解法、牛顿法、修正的牛顿法以及增量法。3.1.1直接迭代法设在第r次迭代运算中：49兰州交通大学硕士论文第三章隧道施工三维数值模拟(3.2)式中：为残余（不平衡）力向量，且；为第次迭代中采用的总刚度矩阵，且；为第次迭代中的节点位

3、移向量。则在第次迭代中，按下式计算改进的位移向量值：(3.3)若在迭代过程收敛，则当时，精确解，。直接迭代法适用于求解很多场问题，但不能保证迭代过程的收敛。3.1.2牛顿法—切线刚度法若式（3.1）的近似解为，则可以将第次迭代所得到的改进解作为新的变量，利用泰勒级数将函数展开，取前两项之和得到：(3.4)式中：(3.5)其中为切线刚度矩阵，而：(3.6)改进的位移向量可根据位移增量向量算得，后者按下式求解：(3.7)49兰州交通大学硕士论文第三章隧道施工三维数值模拟然后进行下一步迭代，直到收敛为止。应该指出的是，只要初始刚度矩阵是对称矩阵，则切线刚度矩阵

4、将始终保持为对称矩阵。而在大变形下的割线刚度矩阵则不一定能保持这种对称性。3.1.3修正的牛顿法—初始刚度法在式（3.7）中，令：(3.8)则式（3.7）可写成：(3.9)因此，在每次迭代中不需要重新计算总切线刚度矩阵。达到收敛的迭代次数一般要多于切线刚度法，但总的计算时间并不一定增加，因为采用切线刚度法时原则上每次迭代都必须重新计算体系的切线总刚度矩阵。对于材料应变软化以及体系中塑性区域发展较大的情况，采用初始刚度法仍能取得迭代求解的收敛，而在这种情况下采用切线刚度法则难以甚至不能达到收敛。3.1.4混合法该法是切线刚度法与初始刚度法联合使用的方法。为

5、此必须采用增量加荷的方法，将总荷载分成几级，逐级加荷。在每一级荷载作用下采用一种初始刚度进行迭代运算，达到收敛后在施加下一级荷载，并采用新的切线刚度矩阵进行迭代运算。如此重复计算，直到收敛到总的荷载为止。3.1.5求解非线性问题的收敛正则判断求解非线性问题迭代运算过程是否收敛的准则如下。⑴残余力向量范数逐渐减小表示运算过程趋于收敛：49兰州交通大学硕士论文第三章隧道施工三维数值模拟(3-10)式中：为向量的二范数。⑵位移增量向量范数逐渐减小表示运算过程趋于收敛：(3-11)⑶残余力向量范数与体系外荷载节点力向量之比小于规定值，表示收敛：　　　　　　　　　

6、　(3-12)式中：为容许值，一般可以取；为体系的自由度总数。⑷节点位移增量向量的范数与节点总位移向量的范数之比小于规定值，表示收敛：　　　　　　　　　(3-13)式中：节点总位移为各级荷载作用下所产生的位移。实践经验表明，为了得到非线性问题的收敛解，采用分段加荷和分段收敛是有效的。为了加速收敛过程，还可以采用一些特殊的技巧，例如在牛顿法和修正的牛顿法中采用“超松弛”技巧。为此，可将修正的位移增量乘以某一常数，它通常等于2，这样可以提高收敛的速度[20~23]。求解非线性问题，还有很多其它的方法和技巧，目前尚难判断究竟哪一种方法最好，也没有哪一种方法能保

7、证在任何情况下都能得到收敛解。49兰州交通大学硕士论文第三章隧道施工三维数值模拟第二节岩体弹塑性有限元理论洞室围岩处于弹性状态还是塑性状态，应视岩石强度及所受应力状态来确定。通常Ⅰ、Ⅱ及Ⅲ类围岩可按线弹性处理。而Ⅳ、Ⅴ、Ⅵ类围岩，一般需按弹塑性处理，或近似模拟成非线性弹性处理。由于岩体的抗拉强度很低，因而岩体还要模拟成低抗拉材料或无拉力材料。岩体中除了有较均匀分布的节理裂隙外，还会有较大规模的软弱夹层，以及断层和破碎带等，所以，岩体中还要模拟这些不同规模的软弱结构面。由此可见，模拟岩体材料是相当复杂的，从洞室工程实际情况出发，本节只介绍各向同性的弹塑性有

8、限元理论。弹性力学和塑性力学都是连续介质力学的分支。弹性力学研究介质在弹性工作阶