隧道全断面和三台阶对比

1、工程计算实践学号:姓名:导师:西南交通大学2017年6月目录 TOC o 1-5 h z HYPERLINK l bookmark7 o Current Document 1、选题背景3 HYPERLINK l bookmark9 o Current Document 2、概述3 HYPERLINK l bookmark11 o Current Document 有限元法概述3 HYPERLINK l bookmark13 o Current Document 有限元软件4 HYPERLINK l bookmark15 o Current Document 3、本构关系和力学参数概述5 HYP

2、ERLINK l bookmark17 o Current Document 本构模型5 HYPERLINK l bookmark25 o Current Document 基本假设和力学参数8 HYPERLINK l bookmark27 o Current Document 4、工程资料以及相关参数的确定9 HYPERLINK l bookmark29 o Current Document 工程概况9 HYPERLINK l bookmark31 o Current Document 计算说明10 HYPERLINK l bookmark33 o Current Document 5、隧道

3、开挖数值计算11 HYPERLINK l bookmark35 o Current Document 计算模型11 HYPERLINK l bookmark37 o Current Document 计算参数15 HYPERLINK l bookmark39 o Current Document 6、数值模拟计算结果17 HYPERLINK l bookmark41 o Current Document 全断面法和三台阶法的位移云图对比分析17 HYPERLINK l bookmark43 o Current Document 全断面法和三台阶法特定监测点的位移对比分析19 HYPERLINK

4、 l bookmark45 o Current Document 7小结24 HYPERLINK l bookmark47 o Current Document 附录一全断面模型ANSYS建模命令流25 HYPERLINK l bookmark49 o Current Document 附录二全断面法FLAC30计算命令流31 HYPERLINK l bookmark51 o Current Document 附录三三台阶模型ANSYS建模命令流34 HYPERLINK l bookmark53 o Current Document 附录四三台阶法FLAC30计算命令流43新奥法中的全断面法和

5、三台阶法对围岩变形影响对比1、选题背景我国是一个多山岭的国家，山丘地区面积占全国面积的四分之三。随着经济 的发展，交通基础设施建设逐步加快。特别是近年来西部大开发战略的实施和四 纵四横高铁的建设，将隧道建设推向了新的高峰，因此隧道在未来几年内将不断 增多.到08年底我国隧道总数己达五千多座，总里程三百多万米，是世界上隧 道最多的国家之一。目前，新奥法在我国的铁路隧道设计和施工中占据主要地位，尤其是长大铁 路隧道。新奥法的特点是把设计、施匚与监测二者结合为一体，通过施匚过程中 的监控量测来动态调整设计、施工参数。对于深埋长大隧道的而言，监控量测具 有重大的经济意义和实际应用价值。因为科学合理的监

6、控量测工作可以迅速、准 确的获取第一手现场数据。工程师们可以根据这些数据指导施工。然而，对于长 大隧道而言，由于其长度长、地质条件复杂，开挖和支护交错进行，围岩应力变 化和支护荷载转换变得较为复杂。在此背景下，基于一定数学模型的数值模拟技 术，以其独特的优点广泛应用于此类工程中。我国隧道建设起步较晚，虽然在设计理论和修建技术方面取得了一定的成果, 但大量工程实际表明，我国公路隧道建设中还存在很多安全隐患，一些匚程技术 问题仍有待进一步研究。本文将以葫芦山隧道为背景，通过对实际监控量测数据 的分析和数值模拟分析，研究长大山岭隧道的围岩变形及其稳定性的特点。2、概述有限元法概述伴随着隧道及地下工程

7、的高速发展，出现了越来越多的新问题，传统的研究 方法和研究理论已经无法对这些复杂的情况进行分析求解，因此有限元法的数值 模拟应运而生。有限元法以离散为基本思想，在求解诸如弹塑性问题和流变等问 题中，以及应力场位移场的藕合问题中都得到广泛的应用。有限元法在分析地下 工程的变形问题时能考虑支护结构与岩体的共同作用，对于研究围岩稳定性具有 实际意义。有限元法的基本思想是将结构物看成有限个单元组成的整体.以单元节点的 位移或者应力作为基本未知量求解。按照选取基本未知量的不同，可分为位移法、 力法和混合法。位移法的基本来知量为节点位移，力法的基本未知量为节点力， 混台法选取一部分节点位移和节点力为基本未

8、知量。在隧道与地下工程中，最常采用的有限单元法包括:线单元、面单元和体单 元三种类型。二节点线单元和在此基础上发展得到的面单兀和体单兀均属于线形 单元。三节点、多节点线单元和在此基础上发展得到的面单元和体单元均为高阶 单元。为了将结构离散化，往往采用各种类型单元的组合，如节点杆单元用来模 拟锚杆，节点梁单元用来模拟喷射混凝土层，备种线单元和体单元分别用来做二 维和三维分析。现在比较常用的有限元分析软件有ANSYS、ABAQUS、FLAC、ADINA、 MIDAS/GTSo有限元软件本文采用FLAC有限元分析软件，对葫芦山隧道进行模拟分析。FLAC3D(Three Dimensional Fas

9、t Lagrangian Analysis of Continua)是美国 Itasca Consulting Gouplnc 开发的三维快速拉格朗日分析程序，该程序能较好地模拟地质材料在达到强度极 限或屈服极限时发生的破坏或塑性流动的力学行为，特别适用于分析渐进破坏和 失稳以及模拟大变形。能够进行土质、岩石和其它材料的三维结构受力特性模拟 和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可 采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够 相应发生变形和移动(大变形模式)。FLAC3D是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土 或其他

10、材料的三维力学行为。三维快速拉格朗日分析将计算区域划分为若干四面 体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果 单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形, 这就是拉格朗日算法，这种算法非常适合于模拟大变形问题。由于无须形成刚度 矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。FLAC3D在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独 到的优点，具体有三点：（1）对模拟塑性破坏和塑性流动采用的是“混合离散法”。这种方法比有限 元法中通常采用的“离散集成法”更为准确、合理。（2）即使模拟的系统是静态的，仍采用了动态运动方程

11、，这使得FLAC3D 在模拟物理上的不稳定过程不存在数值上的障碍。（3）采用了一个“显式解”方案。因此，显式解方案对非线性的应力一应 变关系的求解所花费的时间，儿互与线性本构关系相同，而隐式求解方案将会花 费较长的时间求解非线性问题。而且，它没有必要存储刚度矩阵，这就意味着, 采用中等容量的内存可以求解多单元结构;模拟大变形问题几互并不比小变形问 题多消耗更多的计算时间，因为没有任何刚度矩阵要被修改。3、本构关系和力学参数概述3.1本构模型本文选用的计算模型为Mohr-Coulomb模型。摩尔一库伦模型是弹塑性模型的一种，它假设岩土体应力达到屈服之前是线 弹性应力应变关系，一旦发生屈服，便呈理

12、想塑性。这个模型一直以来用于分析 隧道及地下工程中的问题。摩尔一库伦模型采用的破坏准则是摩尔库伦准则和最大拉应力准则。三个主 力的大小关系为：“ 0,同时产生dU和TU;卸载时dq0,仅产生4厂0。对于A点，无论荷载q增加 或减少，都不会产生dU,仅产生U。在图1-2中，土体加载从。点逐渐到A点，则A点为屈服点。随着应变的 增加，B、C随后将成为新的屈服点。由此可见，应力状态在屈服点上，便意味 着加载时有塑性变形dU产生，卸载时只有弹性变形已当应力状态减小到屈 服点以内时，正负应力增量只引起弹性变形，总塑性应变不变。因此可以认 为屈服点与塑性应变有关。塑性应变成为屈服准则的一个内变量，在简单的

13、应力 状态下可表示为：图3-2弹性模型应力应变曲线莫尔库伦模型的屈眼准则为：t c - c tan (p4 1 4% +.-=Sin。+ c cos(P (3.6) 乙乙塑性理论中，流动规则用来确定塑性应变增量的方向或塑性一年增量张量的 各个分量间的比例关系。塑性理论规定塑性应变增量的方向是由应力空间的塑性 势面g决定的，摩尔库伦模型中的g有两个方程描述区及国。这两个函数分别 用来定义剪切塑性流动和拉伸塑性流动，其中函数号为不相关流动法则，它不能够保证解的唯一性;国为相关联流动法则，它满足经典塑性理论要求的材料稳 定性。笃=%一咽07）_ 1 + sin 的公式中勿为膨胀角，w - j二sin

14、仅）。gS %（3-8）对于弹型性本构模型而言，硬化定律是计算一个给定的应力增量引起的塑性 应变大小的准则，硬化参数一般是塑性应变的函数，即：H = H （sp 、J ） （3-9）硬化参数实质上反映了土中颗粒间相对位置变化和颗粒破碎的量，即土的状 态和组构发生变化的情况，是一种土的状态和组构变化的内在尺度，从宏观上影 响土的应力应变关系。3.2基本假设和力学参数（1）计算基本假设岩体为理想弹塑性介质，均质、连续、各向同性；不考虑岩体变形的时间效应和地下水影响；隧道及围岩的受力和变形是平面应变问题；地层和材料一的应力一应变均在弹塑性范围内变化，地应力场由重力自动 生成；不考虑开挖爆破对围岩的损

15、伤和震动影响。（2）选取力学参数计算中围岩的物理力学指标主要包括变形参数，如弹性模量和泊松比。 还有强度指标，如粘聚力和内摩擦角。这些参数的取值对于数值分析结果影响很大。根据现场的取样实验，得到的具体参数如下:表3.1各材料物理力学参数材料类型容道 (KN/n?)弹性模量 E(GPa)泊松比口粘聚力 (MPa)内摩擦角 (p()V级围岩181.50.30.1523III级围岩2370.250.935I级围岩27210.2L640喷射混凝土24280.2锚杆78.52000.3钢拱架7X32000.34、工程资料以及相关参数的确定工程概况现选某时速80km/h的单线铁路隧道。其衬砌结构形式如图4

16、1所示，拱墙 部为C25网喷混凝土,仰拱部为C25喷射混凝土,内设工16型钢架，每为钢架 间距为0.6m,二次衬砌为C35钢筋混凝土，支护结构参数参考铁路隧道设计 规范取值。隧道处于低山区，自然坡度2025。，植被发育。地层为全风化 弱风化花岗岩，埋深15m,地下水较发育，围岩等级为IV级。图41隧道衬砌结构（单位：cm）计算说明（1）分别计算该隧道在三台阶法、全断面法两种开挖方式下的地表沉降情况。 地表沉降监测以隧道中部地表断面为目标面，比较当学子面开挖到中部断面时的 结果。（2）分别计算该隧道在三台阶法、全断面法两种开挖方式下的各位移情况。包 括A点拱顶位移，C、D点水平收敛位移及R点仰拱

17、位移，如图42。图42断面监测点示意5、隧道开挖数值计算计算模型应用有限差分软件FLAC3D对计算该单线铁路隧道采用三台阶法、全断面 法两种开挖方式地表的沉降、隧道拱顶的位移、边墙的收敛和仰拱的位移。建模 是考虑隧道开挖半径的影响范围，在宽度方向左右各取5倍隧道洞径，向下取为 隧道洞径的4倍，向上取至地表，埋深为15m,在隧道长度方向取90m。首先在AutoCAD中绘制隧道模型的二维图形，划分好区域和围岩土层等, 导入ANSYS软件中实现，建立隧道平面模型，如图51、5-2所示。1TYPE NUMLINESANSYSJUN 22 2017 10:55:58图5-1全断面开挖ANSYS模型ANS

18、YS图52三台阶法开挖ANSYS模型然后利用ANSYS该软件划分网格，划分网格的原则是靠近隧道开挖部分网 格划分较密，远离隧道开挖部分则划分较细，划分为四边形单元，建立隧道平面 模型，如图5-3、54所示。图5-3全断面开挖ANSYS模型网格划分图5-4三台阶法开挖ANSYS模型网格划分在ANSYS中利用VDRAG命令进行拉伸，利用不同的材料类型来为Hac3D5中创建不同的组。采用河海大学郑文棠博士编写的程序ANSYS-FLAC.exe实现，即可生成可读入的flac3D50模型文件，用Hac3D5读取文件后生成的三维模型如图55、5-6所示。图5-5全断面开挖flac3D50模型ZonvCol

19、orby. Group Any group 1group2I group3group4图5.6三台阶法开挖Hac3D5模型计克模型采用位移边界条件，底部边界约束竖向位移，上部为自由边界，左右两端边界约束水平位移，如图57、58所示。ZoneColorby: Group i groupl group2GP FixityColor Index: localco orI xcolor ycolor zcolorScale: 1ZoneColorby. GroupI groupl group2 groups group4 GP Fixity Color Index: localcolor xcolor

20、 ycolor | zcolor Scale 1图5-8对三台阶模型约束条件设置计算参数根据已知资料，拱墙部为C25网喷混凝土，仰拱部为C25喷射混凝土，内 设工16型钢架，每根钢架间距为0.6m,二次衬砌为C35钢筋混凝土，该地层为 全风化弱风化花岗岩，埋深15m,围岩等级为IV级，并根据铁路隧道设计规 范TB10003-2005，参数范围均取中值，可得到该隧道闱岩的物理力学计算参数 如表51。表51隧道围岩的物理力学参数围 岩级别重度 y(kN/m3)弹性模量E(GPa)泊松比V内摩擦角亚。)粘聚力c(MPa)IV21.53.650.325330.45在FLAC围岩稳定性分析中，岩体变形参

21、数采用体积模量(K)和剪切模量(G), 将弹性模量(E)和泊松比(v)通过下列公式转化成体积模量(K)和剪切模量(G)。.E . g- E 3(1-272(1+1；)该隧道的初次衬砌采用锚杆+喷射混凝土的的柔性支护，且内设1：16型钢架, 数值计算中隧道支护初次衬砌混凝土结构的物理力学参数取值如表32。表52隧道初次衬砌混凝土结构的物理力学参数类型重度 y(kN/mJ)弹性模量E(GPa)泊 松比VC25喷射混凝土2.226.590.2表52中喷射混凝土的弹性模量实际上是包括了 16工字型钢拱架的弹性模 量，本次计算中采用将16工字型钢拱架采用等效的方法折算成混凝土的弹性模 量，即C25喷射混

22、凝土弹性模量E=原C25喷射混凝土弹性模量Ei+(16工字型钢拱 架截面面积Sixl6工字型钢弹性模量E?)/混凝土截面积S其中原C25喷射混凝土弹性模量E=23GPa;16工字型钢拱架截面积Si=26131cnf;16 1字型钢弹性模量E2=206GPa;混凝土截面积 S=0.25x90=22.5m2;隧道模型纵向长度为90m,每棉钢拱架间距为0.6m,因此总共有150根钢 拱架，因此可以得到折算后的初次衬砌混凝土弹性模量E=23+(0.0026131x206x150/22.5 ) =26O9GPa0根据FLAC3D数值计算的需要，砂浆锚杆参数的选取如表53。表53砂浆锚杆参数类型电 度 (

23、kN/m0)弹性模量(GPa)抗 拉强度 (kN)刚 度(MPa)粘结力(kN/in)横截 面积(cm2)锚杆23203102010003.801在建模过程中，根据计算结果需要，只需得到地表位移和围岩位移，初次衬 砌可采用shell单元来模拟，厚度为25cm,其它物理力学参数根据表52选取, 锚杆采用梅花形布置，如图59所示。根据锚杆L2mxi.2m布置，选择每个循环 开挖进尺L2m,数值计算中在开挖步200后，应力释放一部分后喷射混凝土和 打锚杆，即支护。CableColorby: Uniform Uniform Shell Colorby: Uniform Uniform图5-9隧道开挖初

24、支结构图6、数值模拟计算结果6.1全断面法和三台阶法的位移云图对比分析采用全断面法隧道开挖时，其竖向位移和水平位移分布云图分别如图61和 图6-3所示，采用三台阶法隧道开挖时，其竖向位移和水平位移分布云图分别如 图62和图64所示。Contour Of 乙DisplacementI 9.6764E-O49.0000E.048.0000E-041 7.0000E-046.0000E.045.0000E.044.0000E.04|u 3.0000E-042.0000E-041.0000E-04 O.OOOOE*OO.-1.0000E-04I .2.0000E.04-3.0000E-0440000E

25、Q.5.0000E-04图6T全断而法胫向位移云图-5.1541E-04Contour Of Z-Displacement I 1.0084E-031.0000E-039 0000E-048.0000E-047.0000E-046.0000E-045 0000E-04I 4 0000E-043.0000E-042.0000E-041 OOOOE-04 O.OOOOE+OO1.OOOOE-04I- -2.0000E-043.0000E-0440000E3-5.0000E-045.6029E-04图62三台阶法竖向位移云图对比分析可以发现，三台阶法的拱顶胫向位移稍大与全断面法开挖隧道拱顶 的竖向位

26、移，但是三台阶法的拱底位移却元小于三台阶法开挖时拱底的隆起位移, 说明三台阶法有利于控制拱底的围岩的稳定，减小拱底隆起的程度。Contour Of X-DisplacementI3.1291E-043.0000E-04 2.5000E-04 2.0000E-04 1.5000E-04 1.0000E-04I 5.0000E-05 0.0000E*005 0000EQ51 0000E-041 5000E-04I.2 0000EXJ42 5000E-043 0000E-O4-3.1409E-04Contour Of X Displacementi 3.1268E)43.0000E-042 5000

27、E-042.0000E-04 1 5000E-04I 1.0000E-04 5 OOOOE-05 0.0000E+00 5.0000E-OS -1.0000E-04 -1.5000E-04 -2 0000E-04 2 5000E-04 -3 0000E-04 -3.1104E-O4图64三台阶法水平位移云图通过图6-3与图64的对比可以发现，全断面法与三台阶法对拱腰的水平位 移的影响相差不多，故三台阶法并不能相比于全断面法对隧道拱腰的影响小。6.2全断面法和三台阶法特定监测点的位移对比分析为具体对比全断面法和三台阶法关于特定位置位移条件的对比，现将flac3D 中关于隧道拱底、拱顶、拱肩的位移

28、监控曲线进行对比分析。|70:号：ku.u0.501.001.502.002.503.003.504.004.50StepxlOF图65全断面法监测点12 （隧道中间拱顶）的胫向位移变化曲线-0.50 1.00 -1.50-2.00-2.50-3.00-3.500.100.200.300.400.500 600.700.800.901.00Stan *1 MS图66三台阶法监测点12 （隧道中间拱顶）的竖向位移变化曲线对比图65和66分析可以知道，三台阶法拱顶的位移相比于全断面法开挖拱顶 的位移较大，但是仅仅相差0.1mm左右，在工程中基本可以忽略不计，与位移 云图中的结果也基本吻合。0.50

29、1.001.502.002.503.003.504.004.50StcikxlQ4IJ2，g.A图67全断面法监测点13 （隧道中间拱底）的修向位移变化曲线9.08.07.06.0 r;5.0?4.03.02.01.00.0 -* 7fif )ffJ0.100.200.300.400.500.600.700.800.901.00StenxlftA5图68三台阶法监测点13 （隧道中间拱底）的竖向位移变化曲线对比图67和68可知，三台阶隧道中间拱底的竖向位移先比与全断面法拱 底位移出现的时间早，这是由于三台阶法上部group2开挖的时对拱底产生影响。 而全断面法只有在隧道开挖快到隧道中间的时候才

30、会对拱底产生影响。虽说最后 隧道中间拱底的最终位移两种开挖方法都差不多，但是三台阶法拱底位移变化相 对全断面法缓慢，有助于围岩应力的缓慢释放，对于隧道施工而言更加安全。0 5 050505050 5 4 4 3 3 2 2 1 1 0 0 n SMV.A图6-9全断面法监测点14 （隧道中间左拱肩）的水平位移变化曲线图610三台阶法监测点14 （隧道中间左拱肩）的水平位移变化曲线对比分析两种方法隧道中间左拱肩的水平位以量都相差不多，说明两种方法 对其围岩稳定性的影响都不大。0.501.001.502.002.503.003.504.004.50Step x10A4图66全断面法监测点14 （隧

31、道中间表面）的水平位移变化曲线图612三台阶法监测点14 （隧道中间表面）的水平位移变化曲线对比图611和图6T2分析可以知道，三台阶法相比于全断面法对于隧道顶部表面的下沉有较大的影响，三台阶法有利于控制隧道开挖上覆土层的位移沉降,说明其对围岩的影响相比于全断面法开挖时小，三台阶法更有利于保持围岩稳定。7小结1、全断面法隧道开挖地表沉降较三台阶法隧道开挖地表沉降大，因此采用 三台阶法开挖隧道更有利于控制地表沉降值，有利于保持国岩稳定。2、采用三台阶法开挖隧道，隧道拱顶的沉降，仰拱的位移及水平位移值均 要小于全断面法隧道开挖。3、虽然三台阶法隧道开挖较全断面法对隧道围岩变形影响较小，但其工序 较

32、全断面法复杂，在围岩条件较好的情况下应灵活变更开挖方法，达到安全、经 济、高效施工。4、数值模拟分析只是一种手段，定性分析隧道开挖围岩的变形特征，本次 计算的位移变形值均较小，分析原因可能是所选取的参数保守，且在计算过程中, 开挖后应力释放时间设置过短及支护施加过快。附录一全断面模型ANSYS建模命令流/prep7!定义单元类型 et,l,mesli200,6 et,2,soli 445 mp,ex,l 365 e9!围岩mp,pixy,1,0.325mp,densJ,2150mp,ex2Q.6e9!隧道mp,prxy2O3 mp,dens,2,2100 *getJcp_mJcp,04ium4

33、nax k,l-Hcp_m,47.897957,33.07352,0 k,2+kp_m,5 8.01695 2 22,0 k,3+kp_m,40.120382,33.073289,0 k,4+kp_m,3022,0k,5+kp_m,46.114999,40.979287,0 k,6-Hcp_m,58.016952,51.77,0 k,7+kp_m,41.902952,40.979287,0 k,8+kp_m,30,51.77,0k,9-Hp_m,40.00808635.651182,0 k,10+kp_m,48.009065,35.655481,0 k,ll-Hp_ni,44.00897639

34、.603131,0 k,12+kp_m,47.897958,3307352,0 k,13+kp_m,48.659881,36.028064,0 kJ4-H lsel,a,2 lesize,all,6 all sei lsel,s,J lsel,a”2 lsel,a,8 lsel,a,ll lesize,aU”,6 all sei lesize27,10 lesize,28,10Iesize,29,10 lesize,30”,10 lesize,58,10 lesize,15,10 lesize7,10 lesize20,10 lesize,30,10 lesize,54,10 lesize,l

35、,6 lesize,50,6 1 esize,51,6 1 esize,5 2,6 lesize24,6 Iesize,49,6 !周闱土体 a,18,14,32,31 a,1422,38,32 a,2226,37,38 a,17,13,14,18 a,21252622 a,15,ll,13,17 a92325,21 a,16,12,ll,15 a,202423,19 a,33,34,12,16 a,34,35 20,12 a,35,3624,20 a,2,8,22,14 a3272,14 a,lM8,27,13a,12,3028,lla122029,30 a,20,19,109 a,10,1

36、91,5 a,52122,8 !隧洞 a,27J,82 a,28 JO,5,27 a,3029,10,28 mshape.OSmshkeyJ amesh,123,1 /view,1,1,1 /replot csys.OkJ 000,0,0,90 1,33,1000 extopt,esize,75,0 extopt5aclear,l typeNmat.l asel5ajoc,z,0 aplot vdrag,all,J3 all seivsel,s,vohi,1,20,1 vplot eslv.s eplot mpchgj .allall seivsel,s,vohi,21 vploteslv.s

37、 eplot mpchg,2,aH all sei vsel,s,vohi,22 vplot eslv.s eplot mpchg3.all all sei vsel,s,vohi,23 vplot eslv.s eplot nipchg,4,aU附录四三台阶法FLACd计算命令流；初始应力平衡model meclunohrprop bulk 3.48e9 shear 1.38e9prop fric 33 coh 450e3 dens 2150set gravity 0 0 -10fix y ran y -0.1 0.1fix y ran y-89.9 -90.1fix z ran z -0.

38、1 0.1fixx ran x -0.1 0.1fix x ran x 87.9 88.1solve elassav31.savrestore 31.sav;iiii state 0iiiixciisp 0 ydisp 0 zdisp 0iiiixvel 0 yvel 0 zvel 0liis resetliis hist rep 50;weiyanliis id 12 gpzdisp 44.01,-45,41.77liis id 13 gpzdisp 44.01,-45,32liis id 14 gpxdisp 39.45,-45,98.97his id 15 gpxdisp 47.47,-

39、45,98.97def excavationloop i (1,75)ufootbegin=l .2-1.2 *iufootend=-1.2*ishellbegin=l.2-1.2*i-0.01shellend=4.2*i+0.01commandmodel null range y ufootbegin,ufootend group group2step 200sei shell id 1 range x 39.45 48.57 y shellbegin shellend z 39 41.77sei shell id 1 prop iso =(26.59e9,0.2 ) thick = 0.2

40、5 dens 2200solveendcommandifi3ufootbeginl =4.8 -1.2 *iufootendl=3.6-1.2*iycablel=ufooteiidl4O.6ycabl=ufootendlshellbeginl=4.8-1.2*i-0.01shelleiidl=3.6-l .2 *1+0.01commandmodel null range y ufootbeginl 5ufootendl group group3step 200sei shell id 1 range x 39.45 48.57 y shellbeginl shellendl z 35 38se

41、i shell id 1 prop iso =(26.59e9,0.2 ) tliick = 0.25 dens 2200sel cable id=l begin=(40,048,ycablel,34.8689)end=(37.0511 ,ycablel ,34.7499) nseg=13sel cable id=l begin=(40.025,ycabl,35.4684)end=(37.0274,ycabl,35.3495) nseg=13sel cable id=l begin=(40,ycablel,36.073)end=(37.0024,ycablel,35.9839) nseg=13sel cable id=l begin=(40.0247,ycabl,36.671)end=(37.0093,ycabl,36.8965) nseg=13sel cable id=l begin=(40.0895ycablel,37.2668)end=(37.1208,ycablel,37.6958) nseg=13sei cable id=l begin=(40.196,ycabl,37.8564)end=(37.2484,ycabl,38.4932) n