大型岩体开挖施工过程及其支护效果的三维模拟分析.pdfVIP

大型岩体开挖施工过程及支护效果的三维模拟分析 蒋宇静张晚兵江崎哲郎 (九州：七学环境系统科学研窀中心福冈 日本 摘要运用数值技术进行岩体F挖模拟已成为施工设计中彳i可缺少的手段，时确立台 ；里的施工方案起着重要作用。但以往的研究中存在着模拟过程不尽切台实际的一些问题． ：#文结台欺岩隧道开挖及三峡船俐开挖实例，进行了与真实开挖最大程度接近的模拟计 算。在几种不同模拟过程的结果对比的基础上，阐明了施工过程的正确模拟方法． ：是键词 三维数值模拟施工过程支护犬型岩体开挖 1 分步开挖分步支护的基本想法 我们知道，地下工程的失稳主要是由于开挖工作引起的应力重新分布超过围岩强度或 固岩过分变形而造成的，而采用适当的开挖步骤和台理的支护手段，将有效地改变应力平 衡过程，进而达到阻止破坏发q、碱少圈岩变形的目的．而对于不同岩性而占，其应力平 衡过程叉具有不同的特点；硬岩在被开挖后，在较短时间内达到应力平衡。故对于硬岩而 言，主要是分析开挖过程中塑性区南否出现，分布范围及规模、开挖面的位移变化较大 位移出现的部位，以便有钳对性地在各步逐次开挖时即紧接着分步支护加固已开挖面，这 样对鲢r束刚开挖岩面的变形与，屯力的进一步释放、防范塑性变形和拉应力的扩展将起到有 效的作用。而对于软岩而言，7芝吱应力平衡所需时间较长，同时还存在一个最佳支|户时间 川题，因而甩?力平衡和位移变一艺历程作为施工过程所关心的问题应予以重点考虑，纾步开 挖、分步支护的模拟分析及其欢粟评价对确立合理的施工方案将起决定性作用； 2兰维施工模拟的基本特点 2．1基本原理 本：研究三维施工模拟分析时采用的FLAC，。摄：基于FLAC2。的基本原理进行三维扩展而 设计融：的软似模块，故其设计雹想、应用及解析特点与二维FLAC计算有着紧密的联系和 相似之处。F二～c程序采用的Lagrangian法是～种分析连续体非线性充变形问题的数值方 法。这种力法遵循连续介质的假设．轨嗣差分格式，按时步积分求解，采用拖带座标系，随 着结构形状的j￡化不断更新座标，允订链续介质有大的变形，Lagrangian法用差分求解，需 划分刚格，与有限元法的不同之处是网格分物理网格和数学网格且互为影射。所划分的网 格只娶有序可以有小规则的形状， 差分方程就其求解的数值方法而言可以分淘两大类：一类是隐式方法；另一类是显式 方法．隐式方法要解一个线性联立力程组，所有的未知数一次求解，因而需要存储一个大 的矩阵，对于非线性问题迁要多次分段线性化后用增量法求解，显式法则不用解方程组． 580 2．2计算过程的实现 大变形差分解析法采用显式时间函数来求解．在每个时步内按下列顺序来进行计算： ①由节点速度求出新的应变速度；②运用结构方程式由应变速度和应力计算出各时步内 的新的应力；③运用运动方程式由应力和失衡力算出新的质点速度和位移。 equation ／／—_|Equilibrium／／] of l、 noitauqmotion)、E(|f velocities flewstresses new and Orforces dispheement 图1 大变形差分显式求解法的计算循环 如表1所示，在各个时步内进行上述循环计算。当计算达到平衡状态时，失衡力趋 近于0。如果有塑性流动产生，则失衡力趋近某一数值。 模型内各个单元的应变速度张量由下式算出： 争一击善(vln：j_十V”∥” (1) 为了计算在微小时间△t时各个单元的应力增量结构方程式使用了增量形H： △6。=△g，．+A6：j (2) 式中：Alg。l=Hi(6。地。)，△g。为相互旋转应力差分，H：为设定函数。△￡