土钉支护系统稳定性数值研三

PAGEPAGE435.土钉作用机理的有限元模拟有限元法自50年代发展至今，已成为解决复杂岩土力学问题的有力工具。有限元法能比较容易地处理各种复杂的几何形状和统一处理各种类型的边界条件，解决难以用解析法求解的力学问题。特别是近些年来，随着计算机应用技术的发展与普及化，人们开发了大量实用的有限元计算机软件，其中既有大型通用软件包，如ADINA、SAP94、MARC等，又有适用于岩土问题的专用软件（如FINAL、FLAC、VIPLEF等），使得有限法已成为解决实际工程问题的一种有效方法。5.1有限元法原理及基本方程【61,62,63】有限元法是通过变分原理（或加权余量法）和分区插值的离散化处理把基本支配方程转化为线性代数方程，把求解待解域内的连续场函数转化为求解有限个离散点（节点）处的场函数。由于有限元法把求解整个求解域内满足支配方程的连续函数问题，转化成在各个单元上寻找合适的插值函数，使它在单个单元上满足支配方程的问题，且单元的形状及材料性能计算公式远比整个求解域简单得多，因此该法较易实现。有限元法根据未知变量的设置类型可分为三种方法。即一种是以单元节点位移作为基本变量的分析方法，称为位移有限元法（位移型有限元），该法主要采用最小位能原理或虚位移原理进行分析；另一种是以单元节点力作为基本未知量的分析方法，称为力法有限元（平衡型有限元），该法常采用最小余能原理进行分析；第三种是以部分节点位移和部分节点力作为未知量的分析方法，称为混合型的有限单元法（混合型有限元法）。由于位移型有限元法在计算机上更易实现复杂问题的系统化，且便于求解和更易推广到非线性和动力效应等其它方面，因此位移型有限元更为广泛地被采用。在此利用虚位移原理来推导有限元法的基本方程。5.1.1虚位移原理如果弹性体在外力作用下出于平衡状态，则任何从平衡位置开始的虚位移，外力在其上所作的虚功等于该虚位移在弹性体内所引起的虚应变能。反之，如果任何从协调变形状态开始的虚位移，都有虚功等于虚应变能，则该变形状态就是弹性体在外力作用下的平衡状态。虚位移是一种约束允许的无限微小的可能位移，它在弹性体内是连续的，在边界面上满足位移边界条件。设所研究的弹性体内受给定体力的作用，它的整个边界面可分为已知面力边界和已知位移边界两部分。在上受给定面力的作用，在上有给定位移。虚位移原理的数学表达式可用下面的虚功方程表示(5.1)式中是外力（包括体力,,和面力,,)在虚位移上所作的虚功，其计算公式为涝(5茎.2)府是弹性体由于讯虚位移引起的奸虚应变能，其豪计算式为却(5睁.3)蠢虚位移在弹性依体内引起的虚多应变满足下面虹的几何关系仍丛士词臭宋薯厕贪(5.4)助创柴教威将奶(协5.4)蜂代入宗(壁5.3)璃式，经过积分枪得垄利用震Green糖公式上式化为贩式中社分别是边界法员线与和坐标轴间夹角映的方向余弦。戴考虑到在给定追位移的边界最上，玩,够则得尤疗群旱渐(5.5)踩根据续,斗将烟(诵6.1.5)厨式代入，可得截如下矩阵表达龟式得钟(搁5.6)节式中责情是虚应变列阵勿，拔是虚位移列阵恢，皂是应力列阵，爪是体力列阵，鸭是面力列阵。稍其中殊由于在上述推置导过程中，没跳有涉及材料的轰本构关系，因武此虚位移原理吃即可以用于线鸡弹性材料，又背可以用于非线屋性材料。昂对于平面应力潜问题，设在力宏的边界陆上给定面力列驶阵通，在定位移边妈界蒙上给定位移列编阵达，板内有与平业面平行的体力冤，虚功方程为低弓(夫5.7)艘式中传许，降，就为虚应变列阵屯，纪为虚位移列阵揪，读为板的厚度。思将材料的应力括应变关系得代入到上式中靠得杯闹(5.8)深鼠对于平面应变概问题，将平面玉应变状态的弹关性系数矩阵代惨人到陪(羞5.8)奥式中，其虚功齐方程仍和旧(毁5.8)略式相同。逆5.1.2银转单元刚度方程栋的建立谦根据虚功方程师(扁5.8)哪式可知，单元今的虚应变能的荡计算公式为嗽界会凤(5.9)牙式中异为单元隐的面积，夹为单元厚度，腥为弹性系数矩旬阵。阵根据单元几何谣矩阵冷、单元应变号、节点位移列桥阵懂的关系有额王叉骤(5.10倾)那将式吨(债5.10)雪代入到垫(举5.9)德中得辈遭驻(5.1骤1)上式可简写成枣夕五化(5.12素)式中从禁泉(斜5.13)谁称为单元刚度妙矩阵教由式驻(役5.8)暖可知单元敬的外力虚功计摇算公式为种始秃(5.14)畜式中第一项积勇分是体力在虚坑位移上所作的比虚功，第二项翁积分是面力在丘虚位移上所作跳的虚功。扣为形函数矩阵送。湖如果计算单元钉不是在边界单什元或在边界上飞没有面力的作亦用，在第二项边积分为零。琴(6.1.超14)撕式可以简写成镜伍看意(5搞.15)式中识(朗5.16)指从上式可以看近出，列阵弦的力学意义是丈单元体力和面鹿力的等效节点独荷载列阵，因筒为它与节点虚肯位移列阵痛的乘积是单元课体力和面力的掏虚功。群式恢(立5.13)抓和式撕(侵5.16)驳是单元刚度矩贪阵和单元荷载漫列阵算式的一舟般形式。坏如果只研究单曲元垂的平衡状态，寸把周围单元对派它作用的边界郑面力和单元内填部的体力用单齐元节点荷载列却阵摇等效表示，则晃由虚位移原理剪有道或寿享贩骡掌聚由于单元节点萝虚位移列阵婶为非零的任意切微小的列阵，绍上式可以写成某抹歇损(5.17)饰上式即为单元祖刚度方程。场5.2歇深基坑开挖土谱钉支护的有限祖元法分析卧教土钉支护技术略一种柔性主动出制约机制的支胁护技术，它通答过插筋（土钉驴）改善了临近岩开挖面附近土岂体的强度，从晒而达到支护的题目的。由于该顽支护技术涉及驾到钉土相互作欠用、群钉效应吵、土拱效应、客应力路径问题劲，因此其机理靠较刚性的被动投支护技术复杂音得多。对于大届多数十分重要纯的问题很难利凑用解析方法予酸以解决。例如写面层上土压力竹的大小问题、残潜在滑裂面的扒问题、支护结愚构的变形问题酒，土钉有效长勤度的确定、土我钉最佳间距的剂确定等。而这萝些问题是工程驻师们最关心的口问题，这些问弱题对于土钉支植护结构的优化鹊设计和安全稳锈定性起着决定廉性的作用。还论有一些土钉支烫护问题也是无控法用数学力学烤上的解析方法赴能够解决的，渠例如支护结构输面层的位移为支何有时为三角盼形的变形模式躺，有时为抛物其线的变形模式思；支护结构为告什么有时呈现搂滑动破坏，有握时表现为水平菌移动，显然，径这些问题与基痰坑开挖支护时吨的应力场、应倍变场的演化有纪十分密切的关姨系，要研究这锋些问题就必须灿研究与基坑开晌挖支护时的应山力场、应变场危的演化规律，磁这就需借助于迎数值模拟的方弓法来研究。有岔限元法是岩土没工程问题中最聋常用、最成熟联的一种数值分导析方法。在深叮基坑支护中，隶有限元分析常译将土体和支护霸体都作为变形券体，分析开挖叨卸载条件下土角体和支护体各限点的应力与位柱移，因此和常贞规设计计算方搞法相比，有限炉元法是一种比槐较精确、适用虎性广、更为合弃理的设计计算掘方法。它除了匀可以从整体上虫分析支护结构叙及周围土体的冤应力与位移形股状，还可适用率于动态模拟计训算，不仅为事趣前设计与方案挠比较而且亦为吐信息反馈施工踢管理提供适时即处理手段。从漠原理上体来讲建，常规设计方削法中存在的问林题在有限元方悦法中都可不同衰程度地得到解地决醒【扑71慧】狗。肃土钉支护技术商可以看成是一览种原位加筋技米术思【泽5,10,6破4唱】驾，人们对加筋舍土已经成功地缠应用了有限元疯法和离散元法或数值模拟。在禽有限元分析中寇，文献【纤66丘】中引入了复税合单元模型扯(赚compos太iteun番itcel嘱lmode算l),莲该模型将土体病和加筋当作一矩种复合土体，芹作为一种材料察来处理。在离烈散单元法中，帽加筋和土体被阅作为相互作用把的离散单元。杜现在有限元窝【兼4,50,6怀4,65,6撕7,68,7壤0钟】壤和离散单元法请【卧52潮】乏数值模拟方法抚也已成功地应瓣用到土钉支护趁研究中。在有消限元模拟中，逮大多都将土钉狡支护作为平面皇应变问题来考体虑，严格地来局讲，土钉支护曾的数值模拟是佣一个三维问题受，只有当土钉婶的水平间距为机零时，土钉支映护采可采用平源面应变分析。索但由于将土钉使支护近似地作揭为平面问题，轮基本上能满足影工程实际的需猎要，且计算的枝误差在可接受栗的范围之内，炕因此用将土钉稻支护作为平面晓问题来考虑是纹切实可行的眨【声4蓝】田。爽5.2.1笔北本构关系的选眠取惑在高层建筑深像基坑开挖与支懒护的数值分析祝方面，目前国招内外已有不少周的用于分析高肌层建筑深基坑奴开挖与支护的眠有限元计算机絮程序，其数值始分析从早先的闹沿用线弹性理坑论过渡到采用柄弹塑性理论计享算，近年来，伏由于计算技术归的发展，考虑准岩土材料和支杠护结构的流变自问题的数值分隶析已经逐步地商应用到实际工晚程计算中，尤酱其在软土地区飞显得更为重要数。但对于土钉叫支护结构来说搞，流变问题不灯会对计算结果条产生多大的影厉响，因为土钉准施工一般适用疯于较好的硬土久层中，且土钉托的施工速度极狼为迅速，因此昂采用非线性弹言性模型或弹塑孟性分析足以满傲足工程实际的酿要求。报非线性弹性模貌型和弹塑性模险型都反映了土遇体的非线性应裹力角-职应变关系特性瓶。非线性弹性抖模型比较简单面，便于考虑土灾体的初始应力碎状态，计算参敞数容易确定，恶有现成的试验动规程可供使用涝，能为一般工椅作技术人员掌轮握，便于推广津，常用的非线销性弹性模型为衬Duncan愤-chang偏模型胃【衔72-74籍】所。线弹性模型左虽然简单，使击用方便，但只感能在一定荷载察条件下适用，纠当应力水平较休高时，土体中钉还会产生塑性乌变形，这时就橡要考虑弹塑性奥应力应变模型摧。目前应用广眨泛的弹塑性应振力桶-掏应变模型主要左有娘Drucke劫r-Prag麦er丝模型、雄Cambri被dge医模型、奸Lade收模型、雄Lade-D桥uncan谋模型、帽子模冻型等置【委72-74距】炊。弹塑性模型遍的基本概念是念：应力在屈服取前只有可恢复证的弹性变形，夹达到屈服后，疯应变由弹性应石变和塑性应变嘉两部分组成，委即削由于深基坑土港钉支护的过程杀是一个分步开朽挖分步支护的扩过程，土体的屯开挖与支护导触致沿开挖线上坐各点的应力释豆放以及相应土题体的变形及应惧力场的变化，贱要正确地模拟寺这种施工工况泻，在计算时应较该采用增量形朵式的本构关系林。用增量关系兵计算时，弹塑陵性模型应变间降的关系可表示俘为爬弹性应变增量酿按照弹性模型秤计算，弹性模掌量与泊松比根舞据卸载与再加管载曲线确定，功塑性应变增量礼根据塑性模型贫理论计算。塑建性模型理论包旅括屈服条件与胃破坏条件，流猾动法则与加工牢硬化定律三部砖分组成。也Drucke园r-Prag锡er勤弹塑性本构模救型是最早提出觉的适用于岩土估材料的弹塑性佣本构模型之一碗，是建立在考栗虑了静水压力仍的广义涌Mises颜屈服准则或序Drucke截r-Prag悲er即屈服准则的基努础上的。它的著最大优点在于域采用简单的方瓦法考虑了静水行压力对屈服强怨度的影响，而到且模型参数少舱，计算简单。伤同时该模型考馋虑了岩土材料绕的剪胀性与扩六容性。因此它张是岩土工程中运常采用的模型惯之一。试验表萄明，土体的剪棋胀性对于土钉缎的作用影响很而大牛【谜6句】念，因此本文选调用了能反映土锅体剪胀性与扩旱容性的遇Drucke欣r-Prag螺er独弹塑性本构模客型。金Drucke屑r-Prag朵er送破坏准则可表蚂示为挪【仪72,74碍】现猜随(5.1驳8)灭式中，闷为历D-P岩准则材料常数非，按照平面应偿变条件下的应番力和塑性变形倘条件，可表示督为枯雅梯辫(5.握19)屿式妥(临5.19)提中面分别为岩土体弟的粘聚力和内张摩擦角。式敞(春5.18)臭中累分别表示应力中状态第一不变傅量与偏应力张匪量第二不变量扭，可分别表示邮为稿根据弹塑性本锻构关系的增量盒形式有代危参(鞭5.20)辉其中附为弹塑性矩阵之，可表示为补【债62,72统】叠屡肚钩(5.21)驻疤塑性矩阵就可由增量理论经导出。假定采梦用相关联流动乡法则，即塑性衔势函数重等于屈服函数巷,井便可以导出普恋遍弹塑性本构券关系矩阵管己绕较(5.22舅)馅上式不仅适用过于各向同性硬盐化或软化材料垄，同时适用于瓣理想塑性材料诊(察)羽。本文采用芦来建立弹塑性袭有限元计算格足式，以改变参溉数制值的方式来近柔似考虑其硬化政（软化）特性屯。即将土体的阻内摩擦角畜修改为剪胀角拢(梦)粒，适当调整如值，就可以使丈计算的剪胀量百与实测的剪胀平量比较接近或支一致。文献【蛙75章】给出了习惯示上采用的风值的关系。沉5.2.2标辅破坏判据固冻破坏判据有两蹦种，一种是根荐据直Drucke馋r-Prag碰er葡破坏准则，即区公式歼(长5.18).览另一种破坏判跳据为拉破坏，耕认为土体具有猜极低的抗拉强牙度，且在受拉追开裂后不再能喇够承受拉应力想。对于土体的舒受拉破坏，从花数值分析的角喜度来说，可以盒认为是塑性的芳一种特殊形式择，土体受拉破蹦坏的应力钩-阀应变关系具有辈与塑性软化类迫似的特性，即幸初始受拉时具壁有线弹性特征评，受拉破坏后垫其抗拉强度为躲零，因此可以结按照塑性分析侨的类似格式，训进行受拉破坏议的非线性分析寿【挂62碎】旗考虑拉破坏的钉条件可写成裳透(投5.23)鲁式中约为主拉应力贿奉为了判断是否奇发生拉破坏，倚通常须按主应离力方向判别。疤主应力与主应减变的关系仍可眨表示为如下的单形成式号:四熊馋起垃对于弹性有焰厕;置宽请涨广对于塑性有软捞根据拉破坏条扒件，仿照弹塑望性本构关系的哗推导，可导出群如下的拉破就坏框“充弹塑性矩竞阵催”嘴【沸62按】医:听当夹为受拉主应力脉且有胜，则有烦属叉(5.骂24)论当埋均受拉且有首况由两条件联立讯推导出弹塑性教矩阵少为挣:差蚁询(会5.25)泪当照均受拉时，则窝该单元完全破荷坏，不能承受侨任何拉力。则婚取残余强度湿值。对于各向谦同性材料，其鼻强度指标可由亮三轴压缩试验专得到，其初始压强度常及残余强度值馅分别对应于峰祥值得强度及残讯余强度，由相灶应的莫尔包络豆线给出。迫亏配对于平面应变中问题，弹塑性唐矩阵可以表示闭为久流败董(怀5.26)但对于平面应力劝问题，弹塑性族矩阵可以表示高为疯盼育志范(5.2冠7)摄5.2.3安疯接触面性能模协拟艇【躲62,76葱】建由于土体与土猾钉的变形模量邀相差很大，因拣此在两者的界疗面上常伴随着距较大的剪应力歪。剪应力的发限挥程度直接影挂响着土钉的轴江向拉力及支护哈结构的内部稳柏定性，因此研址究接触面的性宿能就显得极为炕重要。然而结微构物与土的接炒触面模拟问题挠是难以表达的畅，过去对于土仍体与结构物接走触面特征的研都究，大多采用森下列两种极端字假设之一锯【亡80虹】杰：瞒(欣1)麻接触面十分粗闲糙，土与结构兼物之间无滑动缸可能；割(励2)骑接触面十分光拢滑，不可能产凑生剪应力以阻增止结构物与土誉体间的相对位齿移。泰总体来说，接长触面单元包括竹两种材料模型扣：层状材料单迈元俩(辜薄单元林)翼和接触摩擦型汇单元炕(糟无厚度单元对)蚊。其中快Goodma延n葡单元是一种应朱用较为广泛的戒无厚度接触摩泪擦型单元雨【躺76拆】纯，该单元能模握拟接触面的滑宜动和张裂。薄易层单元是由延Desai年等嘱【扯78艳】康提出的。该接搜触面单元的的跟一个最大特点前是弹性范围内赛引入了独立的祝剪切模量，可喝以考虑滑动破宅坏，拉裂破坏饭。单元厚度的叨选取对于计算湾结果影响较大暖，单元厚度太蔑大，会从物理宋上带来误差；础单元厚度太小蒸，又会从数学达上带来误差况【泡79取】员，如何较为合锐理地确定单元泥厚度，使得两政方面的误差均蚕不大成为该方令法的关键。目谁前还无较好的察方法，只能根培据已有的经验忆来选取。淡对于接触面性润能模拟，本文饮采用改进了的丽Goodma葵n堤无厚度四节点烤平面节理单元澡【雾62,76渔】睬。如图陡6.2-1钳所示。该单元迹承受切线方向多的剪应力盾与法线方向的膨法向应力地，接触面单元草的应力对-掠应变关系对于石局部坐标离可写成司(伤5.28)常式中盲闲为接触面单元朱的切向及法向消刚度；言，楚为单元两侧对径应点的相对切街向及法向位移潮。处假定位移沿节奉理单元长度呈中线性变化，则披可导出单元对牛于局部坐标差的单元刚度矩茎阵为景尝(5.嫩29)丝刚度矩阵向整娃体坐标转化为尚种将境(5.30唇)其中霜沉赔(侨5.31)努咳万美(5拨.32)剖子块监为正交矩阵。方狂上述接触单元辽的缺点是由于罩单元无厚度，胳计算中可能发摔生单元上下面稠相引互迫“诵嵌拳入汪”厕的现象，必须贿对这种嵌入量仍作人为的限制相，以免导致过符大误差。例如猪对于受压情况屡，为避免两种比材料的重叠，霸计算时常取一万个很大的法向袋劲度系数，这记样不可避免地品会给计算结果葱带来一定的误即差棍【垄79背】芦。此外，当接骗触单元的上下怠面发生相对位凯移时，也将产迎生误差。因此梨Goodma庄n锹于液1976政年对该单元作牺了修正，考虑查到相对转角的染影响。假定以仇接触面单元的哑中心作为计算谷点，则其应力荣-碑应变关系表示源为桨甚(5件.33)垃式中，败，确为接触面单元娘中点力矩，盏为相对转角，孕以单元的社3-4过节点的平面逆业时针转为正。子假定位移沿长桌度线性变化，耻可导出接触面柳单元刚度矩阵烂为偏(5.形34)屿式中，胀为单元长度，溪为沿纵向的厚没度，对整体坐们标的转换同式捷(级5.30)-免式乡(谈5.32)恳。震5.2.4挤嚷土钉杆单元模晕拟谊对土钉结构中气刚筋可采用一膜维杆单元进行价模拟。图根5.2涉表示一维杆单抚元桐，节点为寒。登假定杆件的节志点力与杆端位序移之间的关系促式为秩表妖(晃5.35)其中，结设杆单元臭的长度为购，弹性模量为肿，横截面面积颤为俱，与挨轴夹角为罪，如图知5.2池所示，杆件发叛生的位移产，此时抄。根据几何关山系，杆单元缩气短了府，根据材料力艰学理论，该杆庆的轴向压力为迈。由该轴力在念轴和蚊轴的投影，得狗出相应的四个屑节点力为，，栽这是当屡时的四个节点贿力，因而就是鲜单元刚度矩阵浑中的第一列的本4较个元素。同样基可以得出县，柳，锹时相应的节点握力，从而得出齿该矩阵中其余胶三列的元素。淋这样就得出甜由此可见，杆先单元刚度矩阵诵为一个秃4棋×泄4盛阶矩阵。暮5.3纸分步开挖荷载早及位移增量计幕算翅由于土体自重需及地面超载的亿作用，土体处日于一定的初始骑地应力状态，估土体内的应力糟为超载产生的永应力与初始地倍应力之和。由押于初始地应力晶的存在，开挖滴卸载必然引起逃沿开挖边界上炉各点的应力释雄放，从而引起谎周围土体的变蓬形和应力场的渣变化。开挖后风的应力场可看困作初始应力场疑一系列递减过纷程的结果。在激第一步开挖后肠应力和位移分奇别为餐、皮分别为初始应个力和位移，砖、寿分别为第一步援开挖后的应力友和位移增量。泄在第龙步开挖后的应记力和位移分别贼为腊以上是根据开猾挖边界上各点贷的初始应力场熊和位移计算各档边界节点避上尊“妹等效释放荷邀载轨”垃及节点位移。衔对于深基坑工劝程，开挖前佣，将开挖面作纵为应力自由面朱，开挖前自由方面上各点处于伟平衡状态，开悼挖后，被挖去他单元的作用力渡消失，土体失宿去平衡，导致拣周围土体的变饭形和应力场的亏变化。对于开摧挖荷载的计算答可采用霞Mana揉【罩111沉】此提出的下面方闷法计算各阶段跃开挖荷载。氏式中繁为被挖去的单爬元数；雪为单元应变矩慰阵；枪为单元应力矢积量。勾正求得开挖等效佣荷载后，将其行反转方向施加通于系统，相当柄于在原有应力葵场上迭加一组敢反向节点力，护将挖去单元对疑开挖面的作用重相抵消。这一葱模拟开挖边界队效果，确定释拨放荷载的方法湖，最先是由乔轰ncan誉等人提出的，回称堡为绝“粉反转应力释放圈法黑”蒜，是目前开挖注效果模拟中常丹用的方法。根天据漂“捏反转应力释放顶法紫”若计算出来的位逗移就是基坑开亿挖后的位移值尽。类有限元法模拟防中，对于被挖嫁掉单元的处理翁有两种方法。澡较常用的方法乱是将每一步开奥挖掉的单元作哑为轨“迷空单喇元睬”位，令其刚度接选近于零。把开吵挖部分单元赔以腊“姓空单殊元拳”突代替，可能导菌致方程康的较“吴病充态软”饰，为此，常把词与被挖去的节除点相对应的方炸程从总刚度方排程中消去，即岭令这些节点的新位移为零，并资修改其方程。趴另一对开挖单昼元处理的方法网是，将开挖掉欺的单元和节点遵除去，对单元民网络重新排列丰；重新排列单资元和节点的编洗码，则单元和洞编码会形成不稿规则的排列，廉从而给计算带塘来不便，因此流该方法应用较顷少，本文采用蹲的是前一种方存法。晨5.4物支护过程模拟瓣为了模拟支护宇过程，在确定脱离散化网格时蓝，必须考虑各丢步施工的工况眉特征。土钉支藏护过程可以简蹈单地描述为分密步开挖所腊—护熊加筋支护各血—济基分步开挖的过笔程。首先是开她挖土体，土体丛开挖完后，再都进行土钉鲁的怕“土锚舅固适”脑以及喷射面层吗。因此面层单芬元和土钉单元婆是在某一次支若护时被生成而轰参与运算的，侧这些单元在离殃散化时就该考蚁虑进去。对于槽土钉单元采用蛛一维杆单元，滴面层单元采用悠板单元来模拟铜，支护结构和盆土体间的接触住面单元采用修欠正了的挺Goodma贼n佳单元。艳对于每一步支炼护，把与该步拢分支护所对应匪的单元拒(坊开挖后已普是把“系空单区元构”旦)血重新赋给支护呈材料的几何参攀数及力学参数樱。例如，对于姜面层施工模拟则而言，将铺设美面层处对应的附开挖单元赋给惠钢筋混凝土的皂力学参数；对吩于土钉，在设半置土钉处增加狱土钉单元，并性赋给其相应的著几何参数芽(尤土钉的倾角、唤长度、直径等罚)技和力学参数望(迫抗压强度、弹颠性模量等延)译。把每一施工桨步骤的计算结锦果迭加，即得虏最终的结果。嚷严格来说，土雪钉支护问题是得一个三维问题逼，当