基于UDEC平台的裂隙岩体边坡渗流应力耦合分析

根据裂隙岩体边坡渗流的基本特征和基本理论,阐述了离散元软件UDEC渗流—应力耦合的基本计算理论。利用UDEC作为分析平台建立数学模型,在考虑到渗流场与应力场耦合的作用下对裂隙岩体边坡的稳定性进行了全面分析。对边坡的应力、位移、渗流、裂隙分布等方面进行了研究,重点分析了在渗流与应力耦合的作用下边坡的变形破坏规律及自身稳定性。结果表明利用离散元软件UDEC分析流固耦合是一种有效、可行的方法,同时也为裂隙岩体边坡流固耦合方面的研究提供较高的工程应用价值。     

裂隙岩质边坡渗流与非连续变形耦合过程分析

 裂隙岩体中的渗流–应力耦合作用是岩质边坡失稳的重要因素之一。离散裂隙网络(DFN)模型用于研究裂隙岩体渗流,具有概念简单、效率高、适用性强的优点,是研究裂隙岩体渗流问题最为有效的手段之一。非连续变形分析(DDA)方法是专门针对裂隙岩体的非连续特性提出的一种变形场求解方法,能够更加真实地刻画工程岩体。将DFN模拟和DDA方法结合起来,提出基于DDA-DFN的渗流–应力耦合模型,给出考虑裂隙渗流情况下岩体块体系统的瞬时平衡方程,用于研究裂隙岩体变形对渗流的影响和渗流–应力耦合作用下裂隙岩体的变形破坏特征。利用该耦合模型,对一大型水利水电工程边坡稳定性进行分析。结果表明,水库蓄水后,地下水大幅度抬升,渗流–应力耦合作用加剧,导致边坡裂隙岩体中的关键部位发生大变形甚至破坏,进而触发边坡的进一步失稳。实例分析验证了这种方法用于边坡稳定性分析的有效性。     

岩石破裂过程渗透性质及其与应力耦合作用研究

为了分析岩石细观结构特性变化引起渗透性演化对宏观力学行为的影响,并进行渗流应力耦合作用下岩石破裂机制的研究。基于细观损伤力学和Biot经典渗流力学,建立了岩体损伤非线性本构方程和渗透率关系模型,开发出岩石破裂过程渗流-应力耦合分析系统(coupling system of flow and solid in rock failure process analysis,简称F-RFPA2D),拓宽了原有程序RFPA2D的研究领域。这个系统能够对裂纹在萌生、扩展过程中渗透率演化规律及其渗流-应力耦合机制进行模拟分析,把流固耦合问题的研究从应力状态深入到破坏过程。围绕岩石破裂过程中渗透性的演化规律及其渗流-应力耦合作用机理这一课题,开展以下方面的研究工作: (1) 对经典Biot渗流力学做了进一步的考察,验证了建立耦合渗流方程的主要假设,讨论了各种渗流与应力耦合方程及数学模型的适用条件,通过不同深度岩体渗透率工程试验研究,分析了连续介质模型耦合渗流方程参数的物理意义、适用性和测试方法。 (2) 通过岩石应力应变-渗透率全过程实验研究,从细观结构特征揭示出岩石应力峰值前后的渗透性演化规律。基于逾渗理论,通过引入突跳系数这一概念,建立了描述岩石破裂过程的渗流-应力-损伤关系方程。 (3) 秉承RFPA2D关于岩石材料的细观非均匀性的基本思想,在统计     

裂隙岩体非饱和水力应力耦合的不连续介质模型研究

目前对于裂隙岩体饱和水力应力耦合的研究取得了一些进展,但在很多工程领域不能简单地采用饱和渗流分析,而是要考虑岩体饱和-非饱和渗流、应力耦合作用对工程岩体的强度和稳定性的重要影响.因此在总结众多学者对裂隙岩体水力耦合研究成果的基础上,根据DDA力学计算和非饱和渗流计算原理,提出了基于非连续介质方法的--DDA方法的非饱和水力应力耦合模型;并给出了降雨入渗工况下的边坡水力耦合算例.计算结果表明,边坡稳定性随着降雨入渗时间的增加而减小,降雨强度越大,边坡稳定系数的降幅越大;考虑水力耦合时的边坡稳定性要小于不考虑水力耦合时的边坡稳定性,且在同一时刻,若降雨强度越大,考虑水力耦合与不考虑水力耦合的稳定系数差值越大.仿真试验和工程应用表明其计算成果是符合实践规律的,由此说明了所提出的水力耦合模型能正确反映裂隙岩体的水力学特性,验证了该模型是可行有效的,可付诸于实践.     

裂隙岩体渗流与应力耦合分析

岩体变形与水力特性主要决定于岩体中的裂隙分布、密度和尺寸。本文建议的岩体渗流与应力耦合分析方法以裂隙渗流理论和变形本构关系为基础。文末给出重力坝及坝基渗流与应力耦合数值分析算例。     

基于岩体地下流固耦合理论的研究

本文通过对渗流场、应力场与温度场三场耦合的研究得出三场耦合作用的"完全耦合"模式;通过对孔隙介质流固耦合研究得出在此种情况下应力与应变的关系;通过对裂隙岩体流固耦合研究得出了单一裂隙渗流与应力的关系;分析了裂隙岩体渗流与应力耦合的不同数学模型,得出各种数学模型的优缺点以及适用范围。     

岩石破裂过程渗流与应力耦合分析

通过对经典Biot渗流力学做进一步的考察 ,讨论了耦合渗流方程参数的物理意义 ,在Biot基本方程的基础上 ,增加一个反映渗透系数和孔隙变化率关系的耦合方程 ,并结合原有的RFPA弹性损伤本构方程 ,引入渗透率突跳系数 ξ这一概念 ,提出了岩石损伤演化过程渗流 -应力耦合方程 ,开发出岩石破裂过程渗流 -应力耦合分析系统F -RFPA2D,实例分析表明 ,这个系统能够对裂纹的萌生、扩展过程中渗透率演化规律及其渗流 -应力耦合机制进行模拟分析     

含裂隙岩石渗流力学特性研究

岩体中裂隙的存在严重影响着岩体的渗流特性。为了解不同载荷作用对含裂隙岩体渗流性能的影响规律,利用高精度渗流应力耦合三轴试验系统,对含裂隙砂岩和粉砂岩加载及卸载作用下的渗流特性进行试验研究。试验结果表明:(1)加载试验过程中,随着载荷的增大,试样裂隙隙宽逐渐减小,渗透率随之逐渐减小,渗透率与有效围压呈负指数关系;(2)卸载过程中,随着载荷的减小,岩石渗透率逐渐回升,但回升路径明显低于原始路径,路径不重合表明试样中裂隙的变形具有塑性变形的特征。根据试验结果,建立渗透率与有效围压的关系式,并确定关系式中的待定参数。在试验及理论研究的基础上,通过数值模拟分析试样裂隙面渗透率及渗流速度的变化规律。     

考虑非饱和渗流过程的岩体变形规律分析

蓄水初期谷幅变形异常是我国特高拱坝面临的新问题,考虑岩体从非饱渗流和到饱和渗流转变过程中的变形机制亟待研究。基于岩体从非饱和渗流到饱和渗流的过渡过程,推导岩体在非饱和渗流–应力耦合作用下的控制方程,建立非饱和岩体水力耦合模型,并应用于圆柱试件和边坡入渗过程的变形分析。结果表明:岩体从非饱渗流和到饱和渗流转变过程中,饱和度的改变是引起岩体变形的重要因素之一。研究成果对揭示特高拱坝蓄水初期坝体及山体异常变形机制具有重要的指导意义。     

岩石水压致裂过程的耦合分析

岩石的水压致裂过程是典型的渗流与应力耦合作用的破坏过程。在经典Biot渗流力学基本方程的基础上，基于弹性损伤理论，建立了岩石损伤演化过程的渗流-应力耦合模型，运用岩石损伤破裂过程渗流-应力耦合分析系统F-RFPA^2D，对水压致裂过程中裂纹的萌生、扩展、渗透率演化规律及渗流-应力耦合机制进行模拟分析，初步揭示了岩石水压致裂过程的失稳力学行为。     

裂隙岩体渗流场与应力场耦合研究进展与展望

在分析岩体多介质渗透特性的基础上，提出了岩体多介质渗透类型，综述了裂隙岩体渗流场与应力场耦合模型的研究进展，并提出耦合模型进一步研究的方向。     

考虑渗流特性的岩体结构面分形特性研究

 裂隙的连通性和密度是影响岩体渗流特性的重要因素。从岩体渗流研究的需要出发,对计算机模拟的岩体裂隙网络,应用分形几何理论,提出考虑裂隙连通性和密度影响的岩体结构面信息维数的计算方法,建立信息维数与岩体渗透系数的关系,进而可以用信息维数和岩体结构面几何参数来直接推求各向异性裂隙岩体的渗透系数张量。工程算例表明：(1) 考虑渗流应力耦合作用时,用容量维数计算的渗透系数比用信息维数计算的值高出2倍多,说明用容量维数计算岩体渗透系数会夸大裂隙岩体的渗透能力;(2) 信息维数能较好地反映裂隙密度对渗流的影响。信息维数越大,表明岩体内连通裂隙数量越多,因而岩体的渗透性就更大一些。     

高压引水隧洞陡倾角断层岩体高压压水试验研究

结合一抽水蓄能电站高压岔管区高压压水试验,详细介绍了具有陡倾角断层岩体的高压压水试验方法及试验成果。探讨了岩体高压下水力劈裂区扩展的波动理论和渗透破坏特性,明确了岩体与土体渗透破坏形式差异。试验首次对岩体内的渗透压力和变形进行了同步测量,以了解高压条件下岩体的渗透变形特性。在高孔隙压力作用下,岩体的变形发展试验成果表明,岩体应力和稳定性分析过程中孔隙压力的作用不可忽视,因此采用耦合理论对高孔隙压力环境下的工程岩体进行分析有助于降低围岩失效风险。     

开采对底板岩体渗透性的影响

受奥灰水威胁，煤层“带压开采”的关键是必须搞活底板岩体的破坏机理，确定底板岩体的破坏深度。而确定底板岩体的破坏深度通常是依据其渗透性的改变进行的。底板岩体渗透性的改变和底板岩体的应力场密切相关。利用现场实测数据；岩石力学参数和位移反分析的方法，确定了赵各庄矿1237工作面的底板应力场，进而确定了底板岩体的渗透系数，最后，具体分析了底板应力场和渗透系数的关系。     

裂隙岩体渗流场与损伤场的耦合分析

通过引入渗透压力附加柔度张量的概念以及考虑渗透张量的演化行为，详细研究了裂隙岩体渗流与损伤变形的相互作用机理。在此基础上建立了基于两场耦合的裂隙岩体渗流损伤模型。工程算例结果表叫，当存在地下流体流动时，渗流对岩体变形破坏以及损伤变形对岩体渗流场的影响是不可忽略的。     

开挖动态卸荷对节理岩体渗透特性的影响研究

 利用应力波理论分析高地应力条件下节理岩体边坡开挖过程中的动态卸荷效应，探讨卸荷松动后节理岩体渗透特性的变化规律。研究结果表明，高地应力作用下节理岩体边坡开挖过程中的动态卸荷松动效应显著，其影响必须考虑。开挖荷载幅值、岩体弹性模量和节理分布间距等因素通过影响节理开度而影响岩体的渗透特性；岩体中完整母岩与节理岩体的结合部位往往是开挖松动较大的部位，其渗透系数变化明显。对于等间距平行节理组垂直切割的直立坡岩体，若岩体的弹性模量由坡外向坡内呈线性增大，计算得到的节理岩体渗透系数分布符合幂函数衰减规律(由坡外向坡内)；若节理岩体的弹性模量由坡外向坡内无变化，则开挖动态卸荷松动现象不会发生。     

地下硐室围岩应力状态及相关参数测量结果的分析与评价

在深埋及压力硐室工程中,采用水压致裂法测定围岩三维应力状态,宜首先判定围岩应力的分布状况。在具有应力松弛圈、应力集中区与不受开挖影响的正常应力分布域情况下,只有分别进行计算分析,才能更为真实地揭示出硐壁围岩的三维应力分布特征。对于确定压力硐室围岩自身承载能力的水力劈裂测试成果,只有结合裂隙性状、岩层结构状况等进行全面分析,才能准确评价围岩的抗载强度。由于在高压力作用下,岩体中存在的软弱结构面有可能被张裂或扩展,从而改变岩层的透水性,因此只有按照工程运行状态下围岩实际承受的压力进行高压压水测试,才能得到岩层透水性的可靠资料。     

规则岩体结构面的蠕变特性研究

 岩体结构面的蠕变特性研究对解决岩石力学实际问题具有很重要的意义。通过规则齿型结构面在剪切应力条件下的双轴蠕变试验,研究不同爬坡角的结构面在不同法向应力时的蠕变特性,并根据蠕变试验结果,得到结构面剪切蠕变、法向蠕变随法向应力、剪应力、爬坡角的变化规律和时效性,并选取Burgers模型来描述蠕变的剪切变形特性。在模型试验的基础上采用三维有限元程序对具有锯齿状结构面的水泥砂浆块体在荷载作用下的受力性状进行模拟分析。通过有限元分析,探讨规则齿状结构面的应力分布规律。利用ANSYS计算结果,通过最大正应力理论对结构面的破坏进行判别,从而验证在剪切试验中,齿面以拉断形式破坏。     

地下硐室围岩应力状态及相关参数 测量结果的分析与评价

 在深埋及压力硐室工程中，采用水压致裂法测定围岩三维应力状态，宜首先判定围岩应力的分布状况。在具有应力松弛圈、应力集中区与不受开挖影响的正常应力分布域情况下，只有分别进行计算分析，才能更为真实地揭示出硐壁围岩的三维应力分布特征。对于确定压力硐室围岩自身承载能力的水力劈裂测试成果，只有结合裂隙性状、岩层结构状况等进行全面分析，才能准确评价围岩的抗载强度。由于在高压力作用下，岩体中存在的软弱结构面有可能被张裂或扩展，从而改变岩层的透水性，因此只有按照工程运行状态下围岩实际承受的压力进行高压压水测试，才能得到岩层透水性的可靠资料。     

锦屏二级水电站深埋引水隧洞稳定性研究

裂隙岩体的渗流应力耦合机理是岩体水力学的研究重点。基于多孔介质有效应力原理,以体积应变为切入点,建立了考虑岩体变形特性的渗透系数动态演化方程和弹塑性耦合损伤本构模型,并将所建立的本构模型编制程序导入到大型商用有限元程序ABAQUS,研究锦屏二级水电站深埋引水隧洞运营期围岩和衬砌的受力和变形特征。研究成果表明:考虑渗透系数动态演化的力学模型可以更好地反映隧洞运营期衬砌的外水压力特征,比不考虑渗透系数动态演化的计算结果更加合理,由于锦屏4个引水隧洞的间距较大,隧洞开挖和运行放水所产生的裂隙损伤演化区没有形成贯通区域。