叶麓√帮约台钌.PDFVIP

维普资讯 J 叶麓√帮约台钌 第 l2誊 第 1期 岩石力学与工程学报 ro工．12 仉 l 1993年 3月 ^ e辟 J _ R。砖 M ec／sa~sicsd-I 冒 ∞ -， 84一 B8 ／ 岩体力学数值计算方法的现状与展望 周维垣 (请辙 00㈣ — 上 l 岩体结构网络的数值方法 在进行岩体力学数值计算时，一个首要的问题是建立岩体结构的数学模型。由于岩体结构 成因的复杂性和 出露条件的局限性 ，使得岩体结构的定量描述 问磨长期未能得到解决。随着岩 体结构研究的深入，计算机科学的迅猛发展 以及数据处理方法和手段的 日趋完善 ，目前 已经能 够在野外露头和地质探洞结构面实测统计的基础上，建立结构面几何参数 (分布方位、形态与 特征尺寸、密度、开合度等)的概率模型 ，运用蒙特卡洛模拟原理和方法 ，在计算机上模拟产生 符合现场岩体结构面参数统计特征的岩体结构图象和表征岩体结构面特征的各种定量指标， 从而形成了定量刻划岩体结构性质的系统的理论与方法。S．D．Priest和J．A-Hudson于 l981 年发表了结构面的扫描线统计方法 (scanlinesurvey)~13之后 ，J．A．Hudson在 1988年来华讲学 时又曾讲到 ，他们当时已采用钻孔电视照相技术来描述岩体结构特征。其具体做法是 ，在垂直 或水平钻孔中插入 电视摄像机 ，令其沿钻孔轴前移，就可以将孔壁的裂隙情况摄成图象 ，然后 输入计算机。分析这些图象 ，就可以确定岩体结构特征 。J．A．Franklin也在进行类似的研究。 M．Oda等于 1985年提出用组构张量来分析岩体结构 ，即将岩体结构面空间分布的几何参数 (产状、半径、密度等)综合成一个组构张量来表示 ，随后又提出组构张量与弹性波速及 渗透张 量的关系[2]。H．H．Enistein也对岩体结构的模拟做过系统的研究。近来还有 B．B．Mandelbrot 等用分形 (fracta1)理论，从岩体结构网络的统计 自相似性 出发 ，研究其分形维数，借以描述岩 体结构的几何特征 ，从而 由岩石的微结构网络预测岩体结构网络。国内在岩体结构特征的定量 描述方面也取得进晨 ，其成果已应用于坝基、边坡和地下采场等岩体结构分析 中[3]。 2 岩体力学的耦合分析 岩石力学中多种力学作用的耦合分析 ，是近年来为防止渗流水流引起边坡失稳而提出的 课题 。许多边坡失稳和大坝基础失稳 ，常是由于渗流与剪应力的耦合作用引起的破坏。 地下渗流场是决定岩体工程变形、破坏和失稳的基本因素 ，也是直接影响岩体工程数值分 析结果的重要原始参数。但地下渗流场的确定又十分复杂，岩体中渗流主要滑裂隙系统流动 ， 裂隙在岩体 中的产状多具方向性，分布的位置与开合度又具较大的离散性 ，致使岩体的渗透有 着明显的各向异性和因地而异的特点。另外，裂隙内单宽流量与其开度的立方成 比例，裂隙开 度又受岩体应力场控制 ，而岩体应力场又与渗流场有关 ，从而使渗流场与应力场相互耦合。 维普资讯 第 l2卷 第 1期 翊维垣 岩体力学数值计算方法的现状与展望 85 人们在继续重视野外实测的同时，对渗流模型和分析方法的研究也取得进展。在研究的渗 流模型大致有两类 (3-0：一类是离散模型一裂敷水力学等，其中叉可细分为单一裂隙水力学和 裂隙网络水力学两种。对单一裂隙水力学的研究集中在裂隙面粗糙度修正 ，如 Witherspoon (1981)建议 由典型裂隙宽度额率分布函数 出发推求水力等效宽度 tT~szlg0987)提 出用沟槽 流模型反映裂隙宽度的不均一性态。裂隙阿络水力学维二维分析后出现三维分析。具体方法 是 ，在对实测资料分析后 ，用计算机产生平面或空间裂隙同络，WRtke及Lou~(1966)提出线元 素法解平面网络 ，空 间问题中裂隙之间通过交线建立水力联系lLong(