浅谈ANSYS系统在隧道结构计算中应用条件.docVIP

浅谈ANSYS系统在隧道结构计算中应用条件 　　摘要：在公路隧道设计与施工中，为了提前判断在开挖和支护工程中隧道的结构安全性，隧道结构计算的数值研究方法就成为了一种重要的设计依据和施工控制措施。本文提供了一种方法，即利用ANSYS软件模拟隧道结构在开挖个步骤中的计算模式与应用条件。 　　关键词： 隧道结构ANSYS模拟 　　隧道的结构分析是利用工程力学原理，选取合理的介质，通过相似模型体系对其结构进行计算，具体过程一般通过两个途径来进行，其一是利用相似性理论，采取合理的相似系数，在室内通过模型试验来模拟实际的工程问题。其二是数值计算，这种方法伴随着计算机的发展有了长足的进步。目前，伴随着岩土力学的发展，再加上计算机的普遍使用及其性能的不断提高，有限元法成为发展最迅速的用于隧道结构分析的数值计算方法。 　　有限元法先将结构分解为有限的小单元，在每一个单元上，利用弹性力学、弹塑性力学等力学理论建立力学性能参数之间的关系，然后根据位移或者应力协调条件把这些小单元组合起来，求出整体结构的力学特征。因为有限元法是利用矩阵代数方法求解方程组，而矩阵代数建立的方程组非常方便与计算机的存储与求解，所以，有限元法非常适用于分析复杂的地下结构。 　　1模型的建立 　　利用ANSYS来模拟隧道开挖过程，有两种建模方法，一个是建立真三维的模型，三维模型不仅可考虑围岩的流变特性，还能考虑开挖和支护的空间效应，能保证较好的计算精度。但是建模复杂，计算时间长，且费用较高。另一种建模方法是建立二维模型，即按平面应变问题来处理，隧道在长度方向的尺寸比横截面的尺寸大得多 ,在忽略掘进的空间效应及岩石流变效应的影响时 ,计算模型取为平面应变是可行的。另外，可以通过各阶段相应的初始应力释放系数来考虑开挖过程和支护时间早晚对围岩及支护受力的影响。本文采用后者建立有限元模型。 　　相对于整个岩体而言，开挖所引起的应力重分布的区域是有限的，因而要确定计算模型的范围。实践和理论分析表明，对于地下洞室开挖后的应力应变，仅在洞室周围距洞室中心点3～5倍洞室开挖宽??（或高度）的范围内存在实际影响。在3倍宽处的应力变化一般在10%以下，在5倍宽处的应力变化一般在3%以内。所以，计算宽度可确定在3～5倍洞室开挖宽度（或高度）。本文所采用的模型水平方向上隧道两边的长度均取洞跨的5倍为限 ,即计算模型的水平宽度为开挖隧道跨度的7倍；垂直方向上，隧道下方的距离为洞高的3倍 ,而隧道上方按实际地形尺寸。 　　根据隧道围岩的物理力学性质,在本文有限元计算中,围岩材料采用符合Drucker-Prager屈服条件的(D-P)材料来模拟，以计算隧道结构与地层在开挖过程中发生的非线性变形特性。D-P材料是一种理想弹塑性材料，屈服面不随着塑性应变的增加而改变。Drucker-Prager屈服条件是在Mises屈服条件的基础上，再考虑静水应力影响的屈服条件。 　　即： （式～1） 　　式中： 是应力第一不变量,； 　　 是应力偏量第二不变量， 　　 和K是与材料的内摩擦角 和内聚力C有关的材料系数，其数值为: 　　 （式～2） 　　我们认为混凝土和钢材的变形在弹性范围内，采用线弹性模型。在计算过程中,围岩、喷射混凝土和模注混凝土均采用ANSYS程序中的PLANE42加以模拟。对锚杆和钢拱架的模拟下文单独叙述。 　　2边界条件和初始应力场 　　岩体边界上的位移边界条件通常两侧边界约束其水平方向的自由度，竖直自由度不约束；上边界条件一般取为实际的自由边界；下边界的两个位移自由度采取全约束。在内侧的开挖边界上，施加每一步开挖相应的释放荷载。 　　初始应力场由自重应力和构造应力两部分叠加而成。西南交大的王明年等人针对三车道扁坦隧道，通过模型试验和数值分析，得出了“在静水应力场作用下 ,三车道扁坦隧道的极限承载力最高；在自重应力场作用下 ,三车道扁坦隧道处于不利的受力状态”的结论。所以有限元计算时按不利情况考虑，即只考虑自重应力而忽略构造应力。ANSYS具体有两种操作方法：一是把重力作为单独的一个荷载步。在分析的第一步，首先计算出模型在重力作用下的应力场，以此应力场作为后续分析的初始应力场。应当注意的是后面每一步分析得到的位移场都叠加了重力引起的根本不存在的位移场，所以如果需要位移场的结果，要扣除第一步重力引起的位移场。第二种方法是先在弹性范围内单独对模型施加重力荷载，计算得到初始应力场并将各节点的应力写成一个初始应力文件。然后重新开始分析，将之前得到的初始应力文件读入，作为荷载施加到模型上，最后进行接下来的计算。这种方法计算得到的应力场和位移场直接就是所要的结果，不需要另行处理，但要求前后两次计算的模型必须完全一致。 　　3 施工过程的模拟 　　利用ANSYS程序的荷载步