双向正交节理岩体的一个复合蠕变模型

将受一组正交节理切割的岩体结构视为广义宏观复合材料,用“代表元”的方法．建立了复合蠕变本构方程,这对于工程计算和评估岩体结构的力学特性,不仅具有一定的理论意义,而且具有一定的实用价值。所得结果,也适用于那些组构形式相同的复合材料．     

岩体节理非线性法向循环加载本构模型的改进

假定卸载曲线和重加载曲线在法向的起始刚度由前一次加载曲线和卸载曲线决定,从而完全确定岩体节理循环加载本构方程.对模型参数的不同取值,将已有的不同模型统一起来,通过卸载曲线渐近线内移的方式来反映滞回特性.将改进的节理弹性非线性法向变形本构模型拓展到循环加、卸载条件下,考虑节理面几何特性、充填物特性及岩石类型的复杂性,在将节理变形的非线性程度定量化描述的同时,修正以往模型对曲线形式做出硬性假定的方式.新模型的预测结果与试验结果吻合良好,验证了模型的可行性.     

无单元法应用于节理岩体

无单元法是一种新方法,其在岩体问题的应用中,至今报道甚少.根据岩体结构特点和无单元法的优点,建立了节理岩体的二维无单元模型,编制了相应的节理岩体无单元程序,并给出了算例,计算结果表明,该模型是合理可行的,同时给处理岩体节理问题带来了方便.     

互层状岩体及节理岩体的粘弹性模型

采用对互层岩体或节理岩体的不同介质进行“分离、集中”的方法得到了互层或节理岩体的二维和三维粘弹性模型。这个模型在互层或节理岩体的数值分析中是较实用的。     

节理岩体的无单元模型

根据岩体结构特点，建立了节理岩体的无单元模型，编制了二维节理岩体的无单元程序，并给出了算例。计算结果表明，该模型是合理可行的。     

损伤导致的节理岩体渗流模型

岩体内部含有大量节理裂隙,当岩体受到扰动后,这些节理裂隙将张开、扩展并贯通,形成渗流通道,这使得岩体渗透性发生急剧变化。为了描述这种由节理裂隙扩展导致的岩体渗流场发生的变化,将这些节理裂隙看作岩体内部损伤,假设渗透张量与损伤张量主轴重合。在渗透张量主空间内,认为岩体在一个主方向上的渗透系数决定于其它两个主方向上的损伤,建立了相应的损伤与渗透系数之间的关系。在外部荷载作用下,损伤将发生演化,从而导致渗透系数也发生演化。为了描述这种渗透系数演化过程,采用已有的损伤演化方程,通过渗透系数与损伤之间的定量关系,来刻划岩体扰动过程中渗透系数的变化规律。通过与已有数据对比分析,本文所发展的损伤导致的岩体渗流模型是合理的,能有效地反映节理裂隙发育对岩体渗透性的影响。     

加锚节理岩体的无单元模型

锚杆支护技术是一项综合技术，其在岩土工程中发挥着越来越重要的作用。根据岩体结构的特点和无单元法的优点，把锚杆支护应用在节理岩体的无单元模型中，编制了相应加锚节理岩体的无单元程序，并给出了算例，计算结果表明，该模型是合理可行的。     

二向应力作用下断续节理岩体的强度特性研究

本文借助断裂力学理论,着重讨论了在二向应力作用下双向交叉断续节理岩体的压剪断裂强度特性及破坏机理。结合有限元数值模拟计算,分析了节理方位,围压大小对这类岩体强度的影响及其变化规律。     

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

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

地震诱发工程岩体结构面滑移概率分析

基于非平稳随机地震动模型,编制快速合成人工地震动的程序以实现对工程地震动的随机抽样。对工程岩体地震响应进行等效线性数值计算时,定义结构面的水平等效加速度(时程),将结构面水平等效加速度(时程)的最大值超过其屈服加速度(滑移阈值)作为岩体失稳的极限条件;基于地震动随机抽样和等效线性化时程计算,提出一种随机地震动作用下工程岩体沿结构面的滑移概率和稳定可靠度的计算方法。考察某工程岩体在附近3个潜在震源影响下的结构面的滑移概率和稳定可靠性。结果表明,震级相差不大时,远场地震诱发该工程岩体的失效概率小于0.01%;近场地震作用下,工程岩体的失效概率为26.1%。     

考虑临界变形量变化的结合面法向接触刚度计算模型

考虑微凸体接触过程中临界变形量变化,建立结合面法向接触刚度计算模型.采用分形几何表征结合面的形貌特征,由表面形貌测量数据得到分形函数的分形参数D和G.采用弹塑性变形接触理论分析微凸体的接触变形特征,得到微凸体的临界变形量估计模型,其是分形参数和微凸体接触变形量的函数,进而得到考虑临界变形量变化的单个微凸体的接触刚度计算模型.由分形理论得到粗糙表面微凸体分布函数,微凸体分布函数与单个微凸体接触刚度计算模型结合得到结合面的法向接触刚度计算模型.通过对比计算结果与实验数据,验证了所提模型的合理性.     

基于互层岩体的非定常黏弹塑性蠕变模型及有限元分析

为研究互层岩体在考虑时效变形作用时的蠕变特性,提出了一种考虑应力、应变阈值的非定常黏弹塑性五元件模型。该模型能够同时描述岩体的瞬时弹性应变、黏弹性蠕变、线性黏塑性蠕变和非线性黏塑性蠕变（加速蠕变）,推导得出了其在三维应力状态下的蠕变方程;基于ABAQUS有限元软件完成了UMAT子程序的研发,并通过对比岩体的室内蠕变试验与该模型数值模拟结果来验证模型的适用性。结果表明：蠕变试验结果与数值计算结果相吻合,非定常黏弹塑性模型不仅可以精确描述减速蠕变和稳态蠕变过程,还可以较好地描述岩体的加速蠕变过程,从而验证了该模型的适用性与有效性。     

考虑蠕变效应的软岩特大断面隧道围岩抗力系数计算

根据软岩具有蠕变时效的特点,考虑岩体应力-应变全过程曲线,基于鲍格斯(Burgers)蠕变力学模型,推导出了考虑蠕变效应及初始地应力的隧道抗力系数计算公式。以沪昆高速铁路特大断面软岩隧道——明山隧道为工程依托,通过现场水压致裂试验获取隧道初始地应力以及室内岩石试验,获取围岩物理力学参数,代入该公式得到了明山隧道的围岩抗力系数。研究结果表明:该公式能够可靠的运用于具有蠕变效应的软岩隧道衬砌支护设计计算中;当围岩出现塑性流动变形时,围岩仍具有与衬砌共同承载的能力,这一研究结果对隧道衬砌设计具有十分重要的实际指导意义。     

岩石试件非线性蠕变模型及其参数确定

以改进的西原模型为基础，结合岩石的应力应变特征，对模型中的粘滞系数进行了改进，推导出的蠕变方程可描述岩石第三阶段的蠕变特性，同时对蠕变模型中的参数进行了确定；并运用MATLAB程序，结合实例作出改进蠕变方程的曲线及蠕变曲面并与通过实验所得曲线进行比较，改进的蠕变方程更能反映岩石的第三阶段蠕变变形。     

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

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

算井子地段地表岩体结构面发育程度评价

为了探讨结构面发育程度对高放废物地质处置库选址的影响,以甘肃北山算井子地段岩体为研究对象,通过调查地表岩体节理和断层几何特征,采用数理统计方法获得节理产状、密度、间距、迹长和断层产状、断裂带宽度和断层泥厚度等统计数据,提出多尺度结构面综合评价指标来评价岩体结构面发育程度。结果表明:BS22和BS23钻孔周围分别发育4个和5个节理优势组,节理产状服从正态分布;BS22和BS23钻孔周围岩体节理平均中点面密度分别为0.037 m^-2和0.035 m^-2,平均迹长分别为20.38 m和17.39 m;解译出断层F1产状为35°∠85°,延伸长度为4.49 km; BS22钻孔周围岩体结构面发育程度更低,更适合作为我国高放废物处置库建造场址。     

软岩的非线性蠕变模型及其稳定性分析

为了更明确地描述软岩蠕变过程与稳定性的对应关系;采用微分方程描述方法,建立了软岩全程应力应变过程中的非线性蠕变模型,分析了软岩的蠕变过程和对应阶段的稳定性,讨论了蠕变过程中模型参数的变化规律及其对应的稳定性;在不同的应力水平下,对应不同的蠕变过程。广义弹性摸量小于零时,对应第三阶段蠕变。在第二、三阶段蠕变,系统状态均是非稳定的。利用软岩在第三阶段蠕变的对应加速蠕变特性,可为预测、预报系统的失稳现象提供理论依据。     

一种改进的岩石黏弹塑性加速蠕变力学模型

为了全面描述岩石蠕变全过程,克服线性牛顿体不能准确描述加速蠕变的不足,在引入非线性蠕变体模型基础上,结合流变力学模型理论,定义应力与试件长期强度的比值为加速蠕变速率幂级数n,模型发生加速蠕变时的总蠕变量为蠕变特征长度εc,进而得到一种改进的能够描述岩石黏弹塑性加速蠕变的力学模型。结合东乡铜矿砂质页岩单轴压缩下分级增量循环加卸载蠕变试验,对模型参数的辨识进行解释,并将该模型的蠕变拟合曲线与实验的蠕变曲线进行对比。研究结果表明:该模型能很好地描述了岩石的加速蠕变特性。