罗山金矿节理岩体强度的PFC3D模拟及其应用

采用PFC3D数值模拟方法模拟在抗拉、抗压试验下不同倾角层理面的岩体的破坏过程,从颗粒的位移和速度分布等方面出发,阐述了节理岩体的破坏机理,并与实验室条件下获得的节理岩体强度和变形特性进行比较和分析.结果表明,单轴抗压强度在0-60°范围内随倾角α的增加而减小,在60°-90°范围内随倾角的增大而增大;单轴抗拉强度在0...     

新延安隧道层状页岩力学性能试验研究

为研究新延安隧道层状页岩的力学特性和破坏特征，分析层理面不同倾角对岩石力学参数的影响，本文分别进行了不同层理角度下的室内单轴压缩、三轴压缩和巴西劈裂试验及数值单轴压缩和三轴压缩试验。通过室内试验分析，峰值强度随倾角增加呈先减小后增大的趋势；抗压强度呈“U”型趋势，弹性模量呈“V”型变化。通过大量数值模拟试验拓展了室内试验范围并分别拟合了抗压强度和弹性模量随倾角变化关系式；围压的增大弱化了页岩各向异性，但随着围压的不断增大弱化程度减小；黏聚力对层状页岩强度的弱化程度远远大于内摩擦角；抗拉强度随倾角的增大逐渐减小；倾角不同层状页岩的破坏模式也不同，单轴压缩时，0°为张拉破坏，30°~60°时发生由弱面控制的沿层理面发生的剪切滑移破坏，90°时产生劈裂张拉破坏。三轴压缩时不同层理面倾角的页岩试件主要发生剪切破坏。     

单节理岩体抗压强度的数值模拟

为研究单节理岩体抗压强度随节理面倾角及围压变化时的变化规律,利用数值模拟软件UDEC对围压一定时不同节理面倾角的岩体进行了抗压强度实验.结果表明:围压一定,随节理面倾角的增大,岩体抗压强度先减小后增大,在倾角介于45°~60°时抗压强度达到最低值,呈典型的"U"形趋势;节理面倾角一定,随着围压的增大岩体抗压强度也随之增大,整体呈线性变化趋势,可拟合为一次函数,与理论导出式一致;随节理面倾角增大岩体抗压强度增大幅度呈先减小后增大趋势,近似于围压一定时岩体抗压强度随节理面倾角变化规律.     

基于PFC 2D的含节理岩体单轴压缩数值模拟

基于颗粒流方法模拟含节理岩体的单轴压缩破坏过程,并对不同节理倾角的岩体进行单轴压缩数值模拟,研究不同节理倾角岩体的裂纹扩展路径、破坏形式以及强度变化规律。结果表明:(1)在单轴压缩破坏过程中,岩体节理倾角不同,裂纹扩展与最终破坏形态不同,但破坏过程基本一致;(2)当节理倾角从0°到90°的变化过程中,抗压强度随节理倾角的增大先减小后增加,总体呈"U"型分布规律,并在倾角为60°时达最小值。     

节理岩体宏观力学参数尺寸效应的数值试验

基于Mont-Carlo法,建立了节理岩体的结构面网络,采用RFPA2D-DIP数值模拟软件把生成的网络图转化为有限元网格,从而建立能反映材料细观结构的数值模型.对其中6个不同尺寸的节理岩体试件进行单轴受压数值模拟试验,分析了不同尺寸节理岩体应力与应变之间的关系,以及抗压强度和变形模量的尺寸效应,并对岩体的破坏过程进行了研究,描述了1500mm×1500mm岩体的破坏过程.结果表明,随着岩体尺寸的增大其抗压强度和变形模量均减小,但变化幅度却越来越小,抗压强度随尺寸增大趋于稳定的速度更快,并得出了指数衰减的拟合关系式,不同尺寸岩体的破坏过程大致相同.该研究方法为节理岩体宏观力学参数的研究提供了一种新方法.     

板岩单轴压缩各向异性力学特征

考虑岩石材料的各向异性,选用层理比较明显的板岩进行单轴压缩试验,分析力学参数的变化规律。研究结果表明:层理是板岩力学参数、强度特性及破裂模式呈现各向异性的重要原因;板岩单轴抗压强度和弹性模量均随层理面倾角的增大呈"U"型变化趋势,单轴抗压强度的各向异性比为2.8,属中级各向异性水平,试样变形模量的各向异性参数比为1.7;岩样整体泊松比随层理面倾角的增加而增大,但不同试样、相同位置的泊松比变化规律不同;岩样不同倾角的层理面是出现不同破裂模式的重要原因,0°试样形成贯穿层理面的张拉破坏,30°和45°试样基本沿层理面形成单一剪切破坏,70°试样由平行层理面的剪切和贯穿层理面的张拉组成复合破坏,90°试样形成沿层理的张拉劈裂破坏。     

低温冻结条件下板岩破坏类型及单轴抗压强度试验研究

为揭示低温冻结作用对板岩破坏类型及抗压强度的影响,采用DX-40型低温数控试验箱、DNS100微型控制电子万能试验机进行7种不同层理倾角β和6种不同试验温度t的单轴压缩试验,对其应力-应变曲线以及单轴抗压强度、峰值应变、破坏类型的变化规律进行分析。在JAEGER层理面理论的基础上,建立以冻结温度和层理倾角为控制变量的单轴抗压强度公式,给出影响板岩破坏类型的2个极限角度β1和β2的表达式,并通过试验结果验证强度公式的正确性。研究结果表明:受低温冻结作用的影响,板岩的单轴抗压强度随温度降低呈指数增加;板岩的单轴抗压强度随倾角增加先减小后增大;在低温冻结条件下,板岩的破坏类型有3种,即当0°≤β27.0°时,板岩沿与竖直轴线呈一定角度的方向发生剪切破坏;当27.0°≤β≤82.7°时,板岩沿层理面发生剪切破坏;当82.7°β≤90.0°时,板岩沿垂直方向发生劈裂破坏。     

柱状节理岩体相似材料力学特性试验研究

通过柱状节理岩体模型试验探究白鹤滩水电站坝区柱状节理岩体各向异性特征.使用了模拟岩体内部节理的新方式——3D打印制作不同倾角六棱柱形节理试件,以树脂为节理打印材料,水泥砂浆为模型材料,通过立方体试件单轴压缩试验考察节理岩体强度随倾角0°~45°变化的各向异性特征.试验结果表明:岩体压缩强度随着节理面延伸方向与竖直方向的夹角的逐渐增大而减小,单轴压缩下柱体破坏模式主要为剪切破坏.     

层状岩体单轴压缩室内试验分析与数值模拟

为了研究层状岩体抗压强度的层面效应,通过室内试验和数值模拟研究,分析层状岩体单轴压缩情况下的应力应变响应以及强度特征。研究结果表明:随着结构面倾角的增大,层状岩体的压缩强度呈先减小后增大的趋势,含有90°倾角结构面的试样其抗压强度明显大于含有0°倾角结构面的试样对应的抗压强度;数值试验得到的分析结果与室内试验结果相同,岩样的变形特征大致分为弹性阶段、屈服阶段和破坏阶段;压缩强度均随加载速率的增大而增大,并且压缩强度与加载速率之间的关系均可通过线性方程进行拟合,相关系数较大;在实际试验过程中,应根据不同模拟工况采用不同的加载速率。     

考虑节理面法向蠕变的节理岩体蠕变模型

采用Fish语言对FLAC3D软件中的Interface单元进行修正,将Interface单元的法向刚度和切向刚度转化为时间的蠕变函数,便可得到在数值分析中考虑节理面法向蠕变的节理岩体蠕变模型.计算结果表明:节理岩体的单轴蠕变量随着节理面间距的增大而减小;当节理面的倾角小于45°时,节理岩体单轴蠕变量随着节理面倾角的增大而增大,但当节理面倾角大于45°时,轴向蠕变量反而随着节理面倾角的增大而减小;若固定节理面的切向蠕变参数,增加其法向蠕变参数,则节理岩体的轴向蠕变量也会随之减小.     

节理岩体中隧洞开挖的渐进破坏模式

节理岩体中隧洞的破坏模式与节理属性密切相关．基于一个能描述多组节理的岩体模型,采用数值模拟技术研究了节理倾角和侧压力系数对节理岩体中隧洞破坏模式的影响结果表明,节理岩体中隧洞的破坏模式受节理倾角控制．在较小节理倾角下,主要表现为肩部岩体的弯曲断裂、节理之间的分离和边墙岩体的压碎;在较大的节理倾角下,主要以边墙岩体的压碎和节理之间层状岩体的滑移为主侧压力系数亦影响隧洞破坏模式,在较小的侧压力系数下,隧洞破坏主要受节理倾角影响;随着侧压力系数的增大,其破坏主要受构造应力影响．这些结论可为隧洞支护设计和施工提供参考     

基于Hock-Brown准则的层状岩体隧道开挖响应分析

利用Hoek-Brown准则描述岩体的强度特征,并建立相应隧道开挖的数值计算模型,分析不同侧压系数时层状岩体的变形破坏特征.研究结果表明,隧道边墙水平位移随侧压系数K的增大而增大,但当K为2.0-2.5时,位移随侧压系数的变化幅度逐渐减小;顺层岩体一侧的水平位移大于逆层一侧的水平位移;当侧压系数较小时,顶板沉降量大于底板回弹量,当侧压系数较大时,二者之间的差别较小;隧道围岩的破坏形式包括剪切破坏和拉伸破坏,当侧压系数较小时,破坏区域较小,主要位于隧道的左上部;随着侧压系数的增大,围岩的破坏区域逐渐增大,但主要部位仍位于隧道的左上部,而右上部的破坏区域较小.     

多组非贯通裂隙岩石的压缩特性数值研究

裂隙岩体力学性质研究是非常复杂的课题,非贯通裂隙岩石的力学行为与裂隙的几何产状密切相关,裂纹的扩展贯通直接影响岩石的破坏。采用改进刚体弹簧方法,以随机均布Voronoi图为基础,生成不同裂隙产状的岩石数值模型,模拟含多组非贯通裂隙岩石的细观起裂和扩展过程。从细观和宏观两个方面分析了预制裂隙的贯通模式及其对岩石宏观破坏模式的影响;并得到沿平面破坏、阶梯状破坏、沿次生块体转动破坏三种破坏模式。通过分析岩石渐进破坏过程揭示其破坏机理;且定量分析了裂隙岩石在强度上表现出的各向异性特征。将数值计算结果与已有的岩石试验结果对比分析,表明数值与试验结果吻合较好。     

Particle simulation of the failure process of brittle rock under triaxial compression

In order to investigate the failure process of brittle rock under triaxial compression through both experimental and numerical approaches, the particle simulation method was used in numerical simulations and the simulated results were compared with those of the experiment. The numerical simulation results, such as fracture propagation, microcrack distribution, stress-strain response, and damage patterns, were discussed in detail. The simulated results under various confining pressures (0–60 MPa) are in good agreement with the experimental results. The simulated results reveal that rock failure is caused by axial splitting under uniaxial compression. As the confining pressure increases, rock failure occurs in a few localized shear planes and the rock mechanical behavior is changed from brittle to ductile. Consequently, the peak failure strength, microcrack numbers, and the shear plane angle increase, but the ratio of tensile to shear microcracks decreases. The damage formation during the compression simulations indicates that the particle simulation method can produce similar behaviors as those observed through laboratory compression tests.     

深埋垂直板裂结构岩体中洞室失稳破坏机制

采用RFPA数值模拟软件对深埋垂直板裂结构岩体中洞室围岩失稳破坏过程进行了模拟,研究了在不同侧压力系数条件下板裂围岩的失稳破坏特点,并与完整岩体中相同条件下的洞室的失稳破坏情况进行了对比.研究结果表明,侧压力系数对深埋垂直板裂结构洞室围岩的破坏形式具有重要影响,边墙岩柱的溃屈失稳破坏发生在侧压力系数小于1的情况;当侧压力系数大于1时,破坏集中发生在拱顶和隧道底部,边墙岩柱不发生溃屈破坏.指出梁板理论应用于板裂结构岩体洞室稳定性评价的局限性.在深埋情况下,洞室围岩的破坏以应力控制为主,结构面的影响居于次要地位.     

节理岩体爆破破坏过程数值模拟

 针对目前常用的有限元和离散元等数值方法难以客观反映岩体中存在的大量断续节理和在外力作用下岩体破碎及块体运动的不足,提出采用最近新出现的数值方法——数值流形方法以解决目前岩体爆破中存在的上述问题.数值流形方法采用数学网格与物理网格以形成求解流形单元,很容易反映岩体中存在的众多初始节理,本文采用断裂力学准则模拟岩石节理和裂纹扩展,采用DDA中的块体运动学理论模拟块体运动,并通过算例对比分析了完整岩体和节理岩体爆破破坏模式的差异.模拟结果表明,节理存在对岩体爆破破坏模式有着重要影响,且其影响程度与节理的几何分布和物理力学性质有着密切关系.     

渗流作用下柱状节理岩体强度特性试验

为研究白鹤滩水电站柱状节理玄武岩在高渗透压力作用下的强度特性,采用与水电站现场岩体物理力学性质相近的人工材料,制作具有相同节理构造的柱状节理岩体模型试样,开展了不同渗透压力的渗流应力耦合三轴力学试验,分析了渗透压力对柱状节理岩体的破坏模式、强度和渗流量的影响。试验结果表明:渗流作用下柱状节理岩体的破坏模式为沿着节理倾角的滑移破坏;渗流作用下的柱状节理岩体具有较低的强度,三轴压缩峰值强度与渗透压力呈负相关关系;应力应变曲线可分为4个阶段,渗流量随应力应变曲线呈现出阶段性特征,渗透压力对岩体强度和渗流量影响较明显。     

层状岩体强度结构面特征的数值分析

运用FLAC3D软件建立层状岩体试件模型,分析单轴压缩情况下的破坏模式和强度各向异性特征.研究结果表明:随着结构面倾角β的增大,试件抗压强度σc呈现先减小后增大的趋势;当结构面倾角为20°～30°和80°～90°时,试件抗压强度对β的灵敏度最大;当β为30°～70°时,σc变化不大;当结构面内摩擦角φj＜β＜90°时,数值计算结果和理论计算结果差别较小;当β≤φj或β=90°时,两者差别较大,数值计算结果明显小于理论计算结果;数值模拟结果能够反映出β≤φj和β=90°对应的岩体抗压强度存在一定差别,与实际情况相符.