下载文档原格式
收稿日期:2006Ο10Ο26基金项目:国家重点基础研究发展规划(973)资助项目(2002C B412707);国家自然科学基金重点资助项目(5053110)作者简介:孟国涛(1978—),男,四川大邑人,博士研究生,主要从事工程地质、岩石力学与工程方面的研究.岩质高边坡开挖变形的三维离散单元法分析孟国涛1,徐卫亚1,郑文棠1,石安池2(1.河海大学岩土工程研究所,江苏南京 210098;2.中国水电顾问集团华东勘测设计研究院,浙江杭州 310014)摘要:对一典型的反倾向层状结构岩质高边坡进行了非连续非线性三维数值模拟.数值模型中考虑了复杂的地质条件,同时模拟了进水口高边坡的控制性结构面和陡倾节理组,实现了边坡、引水隧洞的分步开挖及锚固.运用三维离散单元法(3DEC )进行了混合不连续和蜕化连续2种方式的求解,在比较不同求解方式所得计算位移的基础上,探讨了结构面对计算结果的影响,分析了进水口边坡在开挖、加固作用下的变形状况和稳定性.关键词:岩石力学;边坡工程;开挖;变形;三维离散单元法(3DEC )中图分类号:T U457 文献标识码:A 文章编号:1000Ο1980(2007)04Ο0393Ο05岩体是经历过变形,遭受过破坏,并赋存于一定地质环境中的地质体,其属性与人工材料有着根本的不同.因此,岩体的固有特性可以概况为“DI ANE ”[1],即非连续性、各向异性、非均质性和非线性.通常条件下,高边坡工程中的卸荷松弛低应力区在总体范围上占主导地位,岩体赋存于应力水平不高的地质环境中,结构面对岩体的潜在破坏方式起到决定性作用,此类边坡问题常表现为结构控制型问题.此外,边坡岩体结构面为成岩建造、构造改造、表生演化的产物,其特征与地质历史有关,多条结构面的几何特征通常存在统计意义上的“优势性”.当此优势性相对显著时,岩体的各向异性特征是不可忽略的.因此,从充分尊重研究对象的基本属性出发,非连续力学数值分析方法就成了解决此类问题的正确途径[2].三维离散元程序3DEC 是处理结构控制型岩体工程问题最为成熟的技术之一[3].该程序不但允许有限位移和离散体的转动及脱离,而且在计算过程中可以自动判别块体之间可能出现的新的接触关系,因此它可以方便地实现对复杂结构体变形破坏的模拟.3DEC 针对问题的性质提供了3种求解方式:蜕化成连续力学解法(类似于F LAC 3D )、完全不连续解法和混合不连续解法.本文借助3DEC 内嵌语言FISH ,采用混合不连续解法对白鹤滩水电站进水口边坡开挖变形进行了离散单元法计算,并将计算结果与蜕化连续力学解进行了比较.1 工程地质条件及工程作用1.1 边坡工程地质条件拟建金沙江白鹤滩水电站位于金沙江下游攀枝花至宜宾河段.河谷的深切侵蚀在坝区右岸形成了高差达600m 级的连续边坡(图1).地层产状总体平缓,层面反坡内倾,地层为二迭系上统峨眉山玄武岩组,以块状玄武岩、杏仁状玄武岩和斑状玄武岩为主.上覆三叠系下统飞仙关组(T 1f ),为一套紫红色河湖相沉积的砂岩、泥岩及砾岩,与下伏二迭系呈微角度不整合接触.岸坡高陡,组成工程边坡的岩石主要为玄武岩和层间错动带之凝灰岩.进水口边坡无区域性控制性结构面通过,主要构造类型为层间错动带、断层、层内错动带和大量的构造裂隙组.其中,区内控制性构造结构面主要为层间错动带及断层,此类结构面贯通性好,延伸贯穿整个坝区,构成坝区构造及岩体结构的基本构架,如图1,2所示.层间错动带C 2～C 11为凝灰质(岩)建造,产状N30°～50°E ,SE ∠15°～25°,具有一定厚度,发育于玄武岩各岩流层之间,凝灰岩易于破碎风化崩解.控制性断层有F 16和F 4,结构面泥化、潮湿.第35卷第4期2007年7月河海大学学报(自然科学版)Journal of H ohai University (Natural Sciences )V ol.35N o.4Jul.2007图1 进水口边坡整体模型Fig.1 I ntegral model of inletslope图2 主要地质结构面Fig.2 Main geological structural plane 进水口边坡岩体统计优势结构面为1组层内错动带和2组陡倾裂隙组,层内错动带的产状为N30°～54°E ,SE ∠6°～18°,陡倾裂隙组的产状分别为N15°～40°E ,NW ∠70°～86°和N30°～50°W ,SW ∠83°～90°.其产状“优势性”明显,岩体各向异性特征显著.图3 开挖模型Fig.3 Excavation model1.2 进水口边坡开挖进水口布置于大寨沟口,自然边坡地形高陡,坡度多在65°～75°之间.进水口位于该陡壁中下部,底板高程约735m ,分11级开挖,开挖高度294m ,如图3所示.由于层间、层内错动带发育,NE 和NW 向断层和节理与层间、层内错动带的组合对边坡稳定不利.此外,大寨沟距坝址较近,对右岸进水口有不利影响,存在高边坡稳定问题.开挖后反倾向岩体的变形特征、陡倾裂隙组的影响以及加固支护方案的合理性是工程建设的重大技术问题.该水电站进水口边坡是一典型的结构控制性问题,应运用非连续力学数值分析方法进行分析.因此,本文应用3DEC 程序建立了进水口开挖高边坡以及引水洞组的联合三维数值模型,并通过此三维数值模型来研究坡、洞在开挖过程中的变形和稳定性问题[4Ο5].2 计算模型和计算参数重视地质结构面对边坡稳定性的影响,考虑结构面的作用是岩质高边坡分析的一个基本思想.由于工程问题的大型混合不连续求解对硬件配置要求太高,出于计算效率的考虑,本文仅取进水口段进行非连续数值分析.2.1 计算模型的建立进水口高边坡的三维计算区域包括了进水口开挖边坡的主要部分.引水隧洞考虑1号～9号机的9个引水洞的进水段.计算范围长840m ,宽240m ,最高点高程为1158m ,底面高程为550m ,模型如图4所示.由于离散元中的结构面有明确的地质意义,即每一条节理都对应于节理地质编录过程中的一个采样参数,因此,节理编录资料可以直接移植到模型中.3DEC 程序对节理的处理是按条进行的,需要指定每一条节理的倾向、倾角、节理面任一点的坐标和节理的连通率.本次研究开展了现场平硐节理详细地质编录,根据测线法原理求得了节理采样的倾向、倾角和位置等关键参数[3].利用3DEC 的隐藏(HI DE )功能,按区域生成长大节理;进而根据已经形成的模型块体尺寸和一般节理的长度,选择合适的连通率参数值,实现对迹长的模拟.生成的1组层内错动带和2组陡倾裂隙组如图4(c )所示.2.2 本构与参数计算中岩体采用M ohr 2C oulomb 准则弹塑性本构方程,层间错动带及控制性断层等结构面采用C oulomb 2slip 模型.493河海大学学报(自然科学版)第35卷图4 3DEC 计算模型Fig.4 Three 2dimensional discrete element calculation model按照蜕化连续力学求解方式进行边坡变形、稳定计算时,岩体中的非控制次要结构面(节理)通过岩体的宏观力学参数来反映,因此在进行边坡岩体分层概化时,主要考虑了平硐勘测所得结构面线密度特征[6].如图4(b )所示:0～40m 为弱风化上段,岩体多为碎裂～镶嵌结构;40～100m 为弱风化下段,岩体多为次块状～块状结构;≥100m 为微新岩体,岩体多为块状～整体结构.计算参数见表1,2.表2 结构面参数取值T able 2 P arameters of structural plane 类型法向刚度/(MPa ・m -1)切向刚度/(MPa ・m -1)黏聚力/MPa 内摩擦角/(°)层间带500600112412F 1650010001082616F 460015001152818层内错动带120040001183510陡倾节理组1120045001153810陡倾节理组2150060001183810表1 岩体物理力学参数T able 1 Physio 2mech anical p arameters of rock m ass 岩体类型变形模量/G Pa 泊松比密度/(g ・cm -3)黏聚力/MPa 内摩擦角/(°)抗拉强度/MPa 玄武岩弱上601272140153817015玄武岩弱下1001252150184717110微新1501222171155315115粉砂岩401302140143510018凝灰岩101352120112616010 考虑与不考虑节理条件的计算模型,其岩体完整程度不同,力学参数也不同.若模型中考虑了陡倾裂隙组,则应适当提高节理切割块体的岩体力学参数值,即岩体的参数值应大于岩体质量分级及参数估计得到的参数值,小于完整岩块的参数值.譬如:将Ⅳ级结构面作为单独地质单元模拟时,岩体的参数即为岩块与Ⅴ级结构面组合的岩体等效参数;而仅将Ⅲ级结构面进行单独模拟时,岩体的参数应该是岩块及Ⅳ级、Ⅴ级结构面组合的等效参数.对于复杂结构岩体,可以采用岩石力学试验与数值方法相结合的直接分析法来获取不同尺度下的宏观岩体力学等效参数[7Ο8];而对于简单的结构岩体,可采用式(1),(2)估算岩体力学参数,且这种估算可通过编写的FISH 命令流来实现.E r =E m K n sK n s -E m (1)G r =G m K s s K s s -G m(2)式中:E m ,E r ———岩体、岩块的弹性模量;G m ,G r ———岩体、岩块的剪切模量;K n ,K s ———结构面的法向、切向刚度;s ———节理的平均间距.3 计算结果3.1 位移特征由图5可知:无节理条件下,计算位移的一致性、连续性较好,最大位移为3616mm.节理条件下,计算位移不连续特征明显,最大位移为4916mm ,位移量级较无节理情况普遍增大.二者的位移特征亦有所差异.无节理条件下,开挖引起的变形深度较大;节理条件下,影响深度较小.无节理条件下,最大位移发生在下部;节理条件下,最大位移发生在中上部位.其原因为:节理条件下,位移特征593第4期孟国涛,等 岩质高边坡开挖变形的三维离散单元法分析693河海大学学报(自然科学版)第35卷图5 计算位移矢量图Fig.5 V ector diagram of calculated displacement更多地受控于结构面的分布,而陡倾角裂隙在边坡中上部的密度更大,其各向异性特征致使中上部岩体下滑趋势更加明显.3.2 结构面影响区内控制性结构面F16断层和层间错动带共同切割形成的断块,具备倾倒—滑塌破坏的边界条件.但由图5可知,该块体未表现出破坏的变形趋势,可视为潜在的不稳定块体.作为次一级结构面,陡倾角裂隙组在边坡开挖的中上部所起的作用相对突出,主要原因是该部位开挖深度不大,浅层破裂岩体未被挖除,结构面的切割有利于个别块体的失稳.总体上,结构面的作用还没有改变边坡的总体变形状态,各断块岩体尚未出现破坏、滑离的位移特点,边坡的总体变形趋势尚保持在正常状态.3.3 锚固效果锚索的主要作用是限制由非连续结构面切割而成的不稳定块体的局部变形,在蜕化近似连续模型中,锚索的作用不明显.在含节理模型中,锚固的作用相对显著.加锚使得水平向最大位移有所降低,同时使塑性区减小.在下部锚索为主的强力加固区,由于锚索水平施加,与近似水平的层间(内)错动带交角很小,进而造成部分锚索置于软弱的层间带中,因此施加的效果不甚明显.中上部浅部加固兼柔性支护区,由于模型中多采用锚杆单元,其长度小,且未设置预应力,所以锚杆单元对岩体变形的约束作用不明显,中上部断块岩体依然表现出较大位移特征.总之,依据节理条件下的计算结果,针对中上部岩体的实际地质条件,在开挖过程中有针对性地调整加固方案,及时实施锚索支护以限制结构面张开,对维持边坡稳定有重要作用.4 结 束 语本文应用三维离散单元法对白鹤滩水电站右岸进水口高边坡进行了非连续非线性数值模拟.模拟结果表明:考虑节理的不连续和各向异性分布特征的计算位移比不考虑节理条件时的计算位移大,中上部岩体的下滑趋势也更明显;结构面尚未改变边坡总体变形状态,边坡整体稳定;锚固作用明显,但尚需针对实际地质条件进行中上部岩体的加固,以保证局部岩体的稳定.本文还将蜕化连续与混合不连续2种求解方法的计算结果进行了对比分析.结果表明:在低应力条件下,结构控制性边坡问题中,节理网络几何特征的大量简化模拟会使计算结果出现显著偏差[9];考虑岩体的非连续性和各向异性,用三维离散单元法进行非连续力学数值分析,能够准确地评价节理条件下工程作用造成的岩体变形与破坏,进而正确评价高边坡的整体稳定与局部稳定.参考文献:[1]H ARRIS ON J P,H UDS ON J A.Engineering rock mechanics[M].Ox ford:Pergam on Press,2000.[2]J I NGLan 2ru.A review of techniques ,advances and outstanding issues in numerical m odelling for rock mechanics and rock engineering[J ].International Journal of R ock Mechanics &M ining Sciences ,2003,40(3):283Ο353.[3]IT ASC A.Three dimensional distinct element code :user ’s guide[K].M inneapolis :I tasca C onsulting G roup ,Inc ,1999.[4]张贵科,徐卫亚,石安池.三峡船闸高边坡及中隔墩变形影响因素分析[J ].河海大学学报:自然科学版,2006,34(3):315Ο320.[5]徐卫亚,宋晓晨,周维垣.水电站进水口岩石高边坡及坝坡与洞室相互作用的三维数值分析[J ].岩石力学与工程学报,2004,23(16):2712Ο2717.[6]巫德斌,徐卫亚.基于H oek 2Brown 准则的边坡开挖岩体力学参数研究[J ].河海大学学报:自然科学版,2005,33(1):89Ο93.[7]MI N K B ,J I NG Lan 2ru.Numerical determination of the equivalent elastic com pliance tens or for fractured rock masses using the distinctelement method[J ].International Journal of R ock Mechanics &M ining Sciences ,2003,40(6):795Ο816.[8]张贵科.节理岩体正交各向异性等效力学参数与屈服准则研究及其工程应用[D].南京:河海大学,2006.[9]朱焕春,BRUM MER R ,ANDRIE UX P.节理岩体数值计算方法及其应用(一):方法与讨论[J ].岩石力学与工程学报,2004,23(20):3444Ο3449.Analysis of excavation 2induced high rock slope deformationwith 3D discrete element codeMENG G uo 2tao 1,XU Wei 2ya 1,ZHENG Wen 2tang 1,SHI An 2chi 2(1.G eotechnical Research Institute o f Hohai Univer sity ,Nanjing 210098,China ;2.East China Hydropower Investigation and Design Institute ,CHECC ,Zhejiang 310014,China )Abstract :Three 2dimensional discontinuous and nonlinear numerical simulation was performed of a typical high rock slope with anti 2dip strata structure.With com plicated geological condition taken into account ,the control structural plane and anti 2dip joint set of the high slope at inlets were simulated ,and the excavation and anchoring of the slope and diversion tunnel in stages were realized numerically.M oreover ,3D discrete element code (3DEC )was used to s olve the deformation of the high slope in the m odes of mixed discontinuousness and degenerate continuousness.The in fluence of structural plane on the calculated result was discussed ,and the stability of the high slope at inlets after excavation and rein forcement was analyzed by com paris on of the calculated results in different calculation m odes.K ey w ords :rock mechanics ;slope engineering ;excavation ;deformation ;3D discrete element code (3DEC )793第4期孟国涛,等 岩质高边坡开挖变形的三维离散单元法分析。
高精度的离散元计算算法及其在岩土力学中的应用研究离散元计算算法是一种基于单元的数值分析方法，适用于解决固体和结构体系的动力学问题。
近年来，随着计算机技术的不断进步，离散元计算算法得到了广泛应用，并在不同领域取得了显著的成果。
在岩土力学中，高精度的离散元计算算法可以为土体的工程问题提供较为精确的预测和分析。
一、离散元计算算法的基本原理离散元计算算法是以单元为基本计算单元的计算方法。
在工程计算中，固体被划分为若干个小的单元，每个单元之间存在着相互作用力，通过对这些单元的受力分析和解算，可以得到整个体系的行为。
离散元计算算法的主要特点是计算对象的离散化，并对离散化单元之间的相互作用力进行分析，从而确定体系的位移、速度、加速度和应力等参数。
离散元计算算法的核心在于求解离散化单元之间的相互作用力。
这些相互作用力来源于接触面接触的分布力、视为完整的体积力、作用于断开箍筋的非连续的力以及强制施加的边界条件等等。
通过对这些力的分析，可以得到整个体系的受力情况，从而解算出位移、速度、加速度和应力等参数。
二、高精度的离散元计算算法在岩土工程中的应用离散元计算算法在岩土力学中的应用主要涉及岩土体的力学性能分析、岩土体的动力响应分析和岩土体的工程可靠性评估等方面。
1. 岩土体的力学性能分析通过离散元计算算法对岩土体的受力行为进行分析，可以得到岩土体的点位移和刚度等参数。
这对于了解岩土体的强度和变形特性非常有帮助，并可以提供参考数据，以制定施工方案和设计方案。
2. 岩土体的动力响应分析岩土体在地震、爆炸、振动等外界环境的作用下会发生动态响应。
通过离散元计算算法对岩土体动态响应的分析，可以预测岩土体在不同工况下的位移、速度、加速度和应力等参数，为岩土工程的稳定性和安全性评估提供依据。
3. 岩土体的工程可靠性评估离散元计算算法可以为岩土体的工程可靠性评估提供数据基础。
通过使用离散元计算算法，可以对不同工况下岩土体的力学性能进行评估，从而为岩土工程的可行性和稳定性评估提供全面而准确的数据依据。
《基于离散单元法的反铲液压挖掘机挖掘阻力研究》篇一一、引言反铲液压挖掘机作为现代工程建设的核心设备之一，其挖掘阻力的研究对于提高挖掘效率、降低能耗以及优化设备设计具有重要意义。
传统的挖掘阻力研究方法往往依赖于经验公式和模型简化，难以准确反映实际工作过程中的复杂力学行为。
近年来，随着计算机技术的发展，离散单元法在岩土工程领域得到了广泛应用。
本文旨在基于离散单元法，对反铲液压挖掘机的挖掘阻力进行深入研究，以期为相关领域的研究和应用提供理论支持。
二、离散单元法概述离散单元法是一种数值模拟方法，主要用于研究颗粒介质、块状介质等散体材料的力学行为。
该方法通过将研究对象划分为多个离散的单元，分析各单元之间的相互作用力、运动状态等，从而模拟整体的运动过程。
在反铲液压挖掘机挖掘阻力的研究中，离散单元法可以有效地模拟土壤的颗粒特性和力学行为，为挖掘阻力的研究提供更为准确的模型。
三、反铲液压挖掘机挖掘阻力研究1. 模型建立本研究采用离散单元法建立反铲液压挖掘机的三维模型。
首先，根据实际设备的结构参数和尺寸，将挖掘机划分为多个离散的单元。
然后，根据土壤的物理特性，如颗粒大小、密度、内摩擦角等，建立土壤的离散单元模型。
最后，将两个模型进行耦合，构建完整的挖掘过程仿真模型。
2. 仿真分析在仿真分析中，本研究主要关注挖掘过程中的挖掘阻力和土壤的变形行为。
首先，通过设定不同的挖掘速度、挖掘深度和土壤条件，模拟不同的挖掘工况。
然后，分析各工况下挖掘阻力的大小及其变化规律。
同时，观察土壤的变形过程，分析土壤的应力分布和位移变化。
3. 结果与讨论通过仿真分析，我们得到了反铲液压挖掘机在不同工况下的挖掘阻力数据。
结果表明，挖掘阻力的主要来源于土壤的剪切力和支撑力。
随着挖掘深度的增加和速度的提高，挖掘阻力呈现上升趋势。
此外，土壤的物理特性如内摩擦角和粘聚力也对挖掘阻力产生重要影响。
在实际工作中，应根据不同的土壤条件和工况，合理调整挖掘机的操作参数，以降低挖掘阻力、提高工作效率。
离散元法在岩土工程中的应用作者：贾书耀张红忠郭超赵栋来源：《西部论丛》2018年第10期摘要：离散单元法是基于不连续性假设的方法，适合求解岩体的不连续性问题，本文就岩体的连续性，离散单元的发展历史和基本原理，以及在岩土工程中的应用给出综述。
关键词：离散单元基本原理不连续数值模拟0 引言岩体是经历过变形，遭受过破坏，具有一定岩石的成分和结构，并赋存于一地质环境中的地质。
在传统研究岩体力学的方法中，通常是将岩体看做连续介质，进而用弹塑性力学的方法进行分析和计算。
但是，在处理存在节理的岩体变形问题时，发现将岩体看做连续介质始终不能符合实际情况。
因此，有学者提出了能够模拟不连续介质的离散单元法。
离散单元法是指将所研究的区域块体划分成一个个互相分离的凸多边形块体单元，各单元之间的连接可以是角角接触，角边接触，角面接触，边边接触，边面接触或者是面面接触，而且随着单元的平动，转动或者变形，允许调整各个单元之间的接触关系。
最终，块体单元可以达到静态平衡或者动态平衡。
1 离散单元法的历史及发展离散单元法的概念最先是由Cundall于1971年提出，主要用于求解节理和块体问题的准静力和动态解。
基于此理论，1978年Cundall和Strack研发出了用于模拟平面颗粒介质的球形单元计算程序 BALL，模拟结果比较符合实测结果，证明了离散单元法对不连续介质本构关系的适用性。
1983年Dowding又进一步提出了边-边（面-面）接触模型，该模型认为单元发生边-角模型接触后也可以分离，且不需要在每个时步内检索单元的接触关系。
1988年Gilbert 研发出了面面接触模型，并给出了基于空间立体单元的基本方程和算法。
同年，Cundall发表了三维离散单元法的基本理论，对块体间的接触判断算法做了详细的论述。
随后日本学者 Kawai 等人基于固体力学提出了新的离散模型，块体为刚性单元，相互之间用弹簧单元链接，可以模拟二维和三维条件下不连续面附近的弹塑性应力场分布。
离散单元法及其在岩土力学中的应用离散单元法（Discrete Element Method, DEM）是一种基于颗粒力学理论的数值模拟技术，广泛应用于工程、地质、材料科学等领域。
该方法通过将介质划分为离散的单元体，并对这些单元体之间的相互作用进行模拟，从而研究材料在各种外力作用下的力学行为。
离散单元法在岩土力学中的应用，能够有效地解决传统有限元方法难以处理的复杂工程问题，特别是在模拟颗粒介质的变形、破坏和流动过程方面表现出显著优势。
离散单元法的基本原理是将介质视为由大量相互作用的离散颗粒组成的系统。
这些颗粒不仅可以自由移动，还可以通过接触力相互作用。
方法通过对每一个颗粒进行力学分析，并考虑颗粒间的接触力、摩擦力、黏附力等，建立数学模型来描述颗粒系统的整体行为。
相较于传统的连续介质力学模型，离散单元法能够更真实地模拟颗粒材料的力学特性，适用于复杂的颗粒介质问题。
在岩土力学中，离散单元法被用于解决各种涉及颗粒介质的工程问题。
岩土工程中的很多问题，如土体稳定性分析、矿山开采、地下工程设计等，都涉及到岩土体在外力作用下的变形和破坏。
离散单元法通过模拟岩土体的颗粒间相互作用，能够为这些问题提供精确的数值解答。
在土体稳定性分析中，离散单元法能够有效地模拟土体在不同荷载条件下的变形和破坏过程。
传统的有限元方法在处理土体大变形问题时可能存在计算效率低和结果不稳定的问题。
而离散单元法通过颗粒间的接触模型，能够模拟土体在荷载作用下的复杂变形过程，分析土体在不同工况下的稳定性，并为工程设计提供依据。
在矿山开采中，离散单元法被广泛应用于模拟矿体的开采过程、评估开采对周围岩体的影响。
通过构建矿体的离散颗粒模型，能够模拟矿体的破裂、变形以及开采过程中可能出现的地表沉降等现象。
这种模拟能够帮助工程师优化开采方案，减少对环境的影响，保证矿山开采的安全与可持续发展。
离散单元法在地下工程设计中的应用也非常重要。
在隧道开挖、地下储层开发等工程中，岩土体的力学行为对工程的安全性和经济性具有重要影响。
一、概述岩石动力学是地球科学中重要的一个分支，研究的是岩石在外部力作用下的力学性质和变形规律。
岩石动力学的研究对工程建设、地质灾害防治、矿山开采等领域具有重要的理论和实践价值。
随着计算机技术和数值方法的不断发展，多尺度离散数值计算方法成为了岩石动力学研究的热点之一。
二、岩石动力学多尺度离散数值计算方法1. 多尺度概念及应用多尺度概念是指在不同尺度下对同一物理系统进行研究，通过层次化的方法，将不同尺度下的物理过程相互通联起来，达到全面、深入理解系统行为的目的。
在岩石动力学研究中，多尺度方法可以将微观尺度的颗粒、微裂纹行为与宏观尺度的整体岩石变形通联起来，为力学行为的研究提供了新的途径。
2. 离散数值计算方法离散数值计算方法是一种在计算机上通过数值模拟的方法对物理问题进行研究的技术。
在岩石动力学研究中，离散数值计算方法可以将连续性的岩石介质离散为有限的物质点或元素，通过求解动力学方程来模拟岩石的变形和破裂行为。
3. 多尺度离散数值计算方法多尺度离散数值计算方法是将多尺度概念和离散数值计算方法相结合的一种研究手段。
它通过对不同尺度下的物理过程进行离散数值模拟，使得微观尺度的颗粒或微裂纹行为能够影响宏观尺度的整体岩石变形行为，从而提高了模拟的准确性和真实性。
三、工程验证1. 工程案例选择在岩石动力学多尺度离散数值计算方法的工程验证中，通常会选择一些具有代表性的工程案例来进行模拟和验证。
例如岩石坡体的稳定性分析、隧道周围岩体的变形预测、岩石爆破工程的影响评估等。
2. 模拟结果与实测数据对比进行工程验证时，研究人员通常会将多尺度离散数值计算方法得到的模拟结果与实测数据进行对比分析。
通过比较模拟结果与实际工程情况的差异，可以评估多尺度离散数值计算方法在工程应用中的可靠性和准确性。
3. 参数灵敏性分析在工程验证过程中，还需要对多尺度离散数值计算方法的参数进行灵敏性分析。
通过改变模型中的一些关键参数，观察其对模拟结果的影响，从而进一步验证多尺度离散数值计算方法的有效性和稳定性。
离散单元法在岩体力学中的应用摘要：岩体是一种地质材料，岩体的力学性质具有各向异性，高度非均质，不连续性等特点。
为了解决工程中遇到的岩体力学问题，数值模拟是岩体力学中常用的手段。
由于岩体中存在大量节理面，基于非连续介质理论的离散单元法更加适合于岩体力学。
本文主要介绍了离散元法块体元的基本原理，以及其在岩石力学中应用范围和应用过程中的典型问题。
最后，提出一些个人见解。
关键词：离散元，非连续介质，岩体力学，数值模拟一．引言当前，我国正处在一个基础建设的繁盛时期，在水利水电，核电，矿山，隧道，地下工程等各领域都会遇到地质环境复杂的岩石力学问题。
为了解决工程中遇到的问题，对于岩体的力学性质有一个较为准确的把握，数值模拟是一个广泛应用的方法。
岩土力学中常用的数值计算方法可以分为两大类。
一类基于连续介质的理论。
如有限元方法，有限差分法，边界元法等。
特别是有限元和有限差分法，应用极为广泛。
连续介质方法对于处理断层、节理、裂隙这样的不连续结构面具有一定的局限性，只能处理为数不多的不连续结构面，例如，在有限元中，岩体中的节理被看作是特殊的节理单元[2]；在有限差分中，岩体中的节理被看作滑移面；在有限元与边界元的耦合中，岩体中的节理被看作是边界面单元。
在这些方法中，对于节理的处理都是小数量、小位移的，因此，对多结构面的不连续介质不适合用连续介质方法模拟，而应采用非连续介质方法进行模拟。
于是离散单元法应运而生。
离散单元法是Cundall 于1971年提出的[3]。
该法将结构面切割的岩体视为复杂的块体的集合体，允许各个块体平移或者转动，甚至相互分离。
离散元法以受裂缝切割或分立的块体为出发点，块和块之间的相互作用在角和面上有接触，角点可以有较大的位移。
在某些情况下如滑坡或冒顶时，岩块可以滑动甚至脱离母体而自由下落。
二．离散单元法原理介绍离散元法的单元从几何形状上分类可分为块体元和颗粒元两大类,本文主要介绍块体元在岩石力学中的应用。
文章编号:1007-2993(2005)04-0177-05离散元法及其在岩土工程中的应用综述王卫华 李夕兵(中南大学资源与安全工程学院,湖南长沙 410083) 【摘 要】 离散元法是基于不连续性假设的数值方法,它特别适合于求解节理岩体中的非连续性问题。
在介绍离散元法基本原理的基础上,着重对离散元法在岩土工程领域的应用现状作了叙述和分析,并对其发展趋势进行了探讨。
【关键词】 离散元法;岩土工程;数值方法【中图分类号】 T B 115A Review on Fundamentals of Distinct Element Methodand Its Applications in Geotechnical EngineeringWang Weihua Li Xibing(School of Resources and Safe ty Engineering ,Central South U niversity ,Changsha Hunan 410083China )【Abstract 】 The Distinct Element M ethod (DEM )is a discontinuum -based numerical method especially applicable to solve the discontinuity problems in jointed rock mass .Firstly the fundamentals of DEM are introduced ,and then its applications in geo -technical engineering are summarzied emphatically ,and finally the development trends of DEM are discussed .【Key Words 】 Distinct Element M ethod ;geo -technical eng ineering ;numerical method0 引 言岩体是一种具有不连续性、非均质性、各向异性和非线性的天然地质体[1～3]。
基于离散元方法的岩石力学性质研究引言离散元方法是一种用于模拟和预测固体材料和结构中颗粒行为的数值方法。
在岩石力学领域，离散元方法可以模拟岩体的损伤、变形、断裂和破坏等过程，为岩石勘探、开采和工程建设提供重要的理论支持。
本文将从离散元方法的基本原理、岩石物理力学性质的离散元模拟、岩石断裂和动力学性能的离散元模拟三个方面，讨论基于离散元方法的岩石力学性质研究。
离散元方法的基本原理离散元方法是一种基于颗粒行为的数值方法，其基本原理是用离散颗粒来表示材料或结构，通过颗粒之间的相互作用来模拟材料或结构的力学行为。
离散元方法的模拟区域通常被划分为小的离散块，每个块都包含若干个颗粒，颗粒之间通过刚性或弹性的联系相连。
离散元方法的模拟过程可以分为两个阶段：初始化和时间步进。
在初始化阶段，首先确定岩石模型的几何形状、颗粒形状和颗粒间联系的类型以及岩石样本的初始状态，然后确定时间步长和重力加速度等物理参数。
在时间步进阶段，根据颗粒间的相互作用力计算颗粒的运动和变形情况，更新颗粒的位置和速度，进而计算岩石模型的宏观力学行为。
岩石物理力学性质的离散元模拟岩石物理力学性质是岩石力学研究的基础，其离散元模拟主要包括密度、弹性模量和泊松比的模拟。
在密度的模拟中，通过将岩石模型划分为若干小块，计算每个块的体积和质量，并根据质量和体积的比值计算密度；在弹性模量和泊松比的模拟中，通过施加一定大小和方向的应力，计算颗粒的变形，进而计算岩石的应力应变关系。
离散元模拟可以较好地模拟岩石物理力学性质的变化规律，特别是在不均质岩体和含裂隙岩体中的应用效果更加显著。
此外，结合数值最优化等方法可以对岩石物理性质进行反演分析，以获得更加准确的模型参数和物理特征。
岩石断裂的离散元模拟岩石断裂是岩石力学中的重要问题，其离散元模拟已成为一种常用的手段。
岩石断裂的离散元模拟主要包括单轴压缩、双轴压缩、剪切等岩石试验和不同尺度下的岩石断裂模拟。
在离散元模拟中，通过施加一定大小和方向的应力，模拟岩石的压缩和剪切过程，进而分析岩石的变形和破坏机制，比如岩体的剪切带形成、脆性破裂和粘性破裂等。
岩土力学与工程离散单元法岩土力学与工程离散单元法引言岩土力学与工程离散单元法是研究土壤和岩石力学性质以及工程结构行为的重要方法之一。
这种方法通过将土壤或岩石划分为离散的单元，以模拟其力学行为，可以更准确地分析和评估土壤和岩石的变形和破坏。
本文将从简单到复杂、由浅入深地介绍岩土力学与工程离散单元法的基本原理、应用领域和发展趋势。
1. 基本原理岩土力学与工程离散单元法的基本原理是将土壤或岩石划分为离散的单元，该方法的关键是通过数学模型和计算方法来模拟单元之间的相互作用和力学行为。
每个单元具有自己的力学特性和应力变形关系，通过计算每个单元的力学行为并综合所有单元的结果，可以得到整个土壤或岩石体系的力学响应。
离散单元法的主要特点是考虑了土壤或岩石的非均质性和非线性性，可以模拟土壤或岩石的各种力学行为，如弹性、塑性、破坏等。
离散单元法还可以考虑复杂的边界条件和加载方式，使得模拟结果更加真实和准确。
2. 应用领域岩土力学与工程离散单元法在土木工程领域具有广泛的应用。
以下是几个典型的应用领域：2.1 地基工程地基工程是岩土力学与工程离散单元法的重要应用领域之一。
通过离散单元法，可以模拟土壤的本构关系、变形特性和破坏机理，对地基工程进行稳定性分析、变形预测和破坏评估。
这对于选择合适的地基处理方法、设计土方工程和预测工程变形具有重要意义。
2.2 岩土工程岩土工程是另一个重要的应用领域。
离散单元法可以模拟岩石的力学行为和破坏模式，对岩石边坡稳定性、岩体固结、地下岩石开挖等问题进行分析和评估。
这对于岩土工程设计和施工的安全性和效益具有重要意义。
2.3 地震工程地震工程是离散单元法的又一个重要应用领域。
通过模拟岩土体在地震作用下的动力响应和地震破坏过程，可以评估地震对工程结构的影响和风险，并优化地震设计和抗震措施。
离散单元法在地震工程中的应用对于提高工程结构的抗震性能具有重要意义。
3. 发展趋势随着计算机技术的不断发展和计算能力的提高，岩土力学与工程离散单元法正朝着更加精确、高效和多尺度的方向发展。
