首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 187 毫秒
1.
陈卫忠  陈培帅  王辉 《岩土力学》2011,32(Z2):573-579
岩石断裂力学在岩土工程领域应用越来越广泛,其数值模拟技术是研究的主要方向之一,也是制约断裂力学发展的瓶颈。基于扩展有限元思想,依托ABAQUS用户子单元二次开发技术,引入Heaviside函数模拟不连续位移场,采用裂尖 逼近函数模拟裂尖位移场,通过J积分求解裂尖的应力强度因子(SIF),最后采用最大周向应力准则,预测裂纹扩展方向,实现了裂纹动态演化过程的扩展有限元数值模拟。通过3个经典断裂力学模型,对比了计算值与理论值,验证了所开发程序的适用性和有效性  相似文献   

2.
断裂问题的扩展有限元法研究   总被引:3,自引:0,他引:3  
扩展有限元(extended finite element method,XFEM)是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。  相似文献   

3.
模拟三维裂纹问题的扩展有限元法   总被引:4,自引:1,他引:3  
余天堂 《岩土力学》2010,31(10):3280-3285
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。  相似文献   

4.
模拟岩石压剪状态下主次裂纹萌生开裂的扩展有限元法   总被引:1,自引:0,他引:1  
师访  高峰  李玺茹  沈晓明 《岩土力学》2014,35(6):1809-1817
压剪应力状态下,岩石类材料中常见两类裂纹:翼型张拉裂纹和次生压剪裂纹。基于扩展有限元方法(XFEM),提出了模拟压剪裂纹面作用机制的扩展有限元位移增强方案,并给出了扩展有限元法分叉裂纹处理方法,分别用最大周向拉应力准则和Mohr-Coulomb准则判断张拉裂纹和压剪裂纹的萌生和扩展。基于Matlab平台编写了数值计算程序Betaxfem 2D,通过两个算例对所提方案进行了验证,所得结果与有限元法(FEM)计算结果吻合很好。模拟了单轴压缩载荷下含预制闭合裂纹试件的裂纹分叉、扩展过程,与试验结果的对比表明,所提方案可以模拟和预测岩石类材料张拉、压剪交互分叉裂纹的萌生和扩展行为。  相似文献   

5.
杨石扣  张继勋  任旭华  张道法 《岩土力学》2016,37(10):3017-3025
运用三维数值流形法(3D NMM)进行三维裂纹扩展分析,并采用C++语言编写了相应的程序。充分利用三维数值流形法模拟裂纹扩展的优势,只需要更新裂纹尖端线附近的边界环路和流形单元,不需要使用阶跃函数。根据三维数值流形法计算得到的应力结果,应用非局部求迹方法分析每个裂纹尖端的破坏状态,如果发生破坏则沿垂直于其最大主应力方向扩展。针对裂纹扩展后的不同状况,采用四边形或三角形推进法。裂纹扩展后为了使变形后的面保持平面,必须对新生成的面进行三角化分割。对诸如单边裂纹、平行钱币型裂纹和倾斜钱币型裂纹扩展问题进行数值模拟。计算结果表明,采用三维数值流形法进行裂纹扩展模拟是可行的,文中方法对裂纹尖端线非闭合和闭合的情形均适用,且文中方法对于裂纹尖端线位于单元内部的非平面裂纹扩展也是有效的。  相似文献   

6.
马文涛  李宁  师俊平 《岩土力学》2012,33(12):3795-3800
针对线弹性断裂力学问题,提出扩展径向点插值无网格法(X-RPIM)。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。  相似文献   

7.
考虑裂隙水流与岩体开裂之间的相互作用,在含裂隙单元的附加节点上引入反映裂隙局部特性的附加函数,然后基于考虑裂隙面水压力作用的虚功原理推导出了采用扩展有限元法分析水力劈裂问题的控制方程,给出了裂隙水流与岩体开裂相互作用的扩展有限元实现方法,并采用该法对岩石试件与压力隧洞在裂隙水压作用下的开裂过程进行了数值模拟。数值计算结果表明:岩石试件裂纹内水压变化规律与试验结果相吻合;裂隙长度较小时,裂隙的存在对压力隧洞周围岩体应力、位移场的影响较小,随着裂隙的扩展,其影响范围将逐步扩大;水力劈裂对隧洞的径向位移影响较小,而对环向位移的影响较大,考虑水力劈裂耦合分析得到的环向位移要大于不考虑裂隙内水压分析得到的结果。  相似文献   

8.
考虑裂隙水流与岩体开裂之间的相互作用,在含裂隙单元的附加节点上引入反映裂隙局部特性的附加函数,然后基于考虑裂隙面水压力作用的虚功原理推导出了采用扩展有限元法分析水力劈裂问题的控制方程,给出了裂隙水流与岩体开裂相互作用的扩展有限元实现方法,并采用该法对岩石试件与压力隧洞在裂隙水压作用下的开裂过程进行了数值模拟。数值计算结果表明:岩石试件裂纹内水压变化规律与试验结果相吻合;裂隙长度较小时,裂隙的存在对压力隧洞周围岩体应力、位移场的影响较小,随着裂隙的扩展,其影响范围将逐步扩大;水力劈裂对隧洞的径向位移影响较小,而对环向位移的影响较大,考虑水力劈裂耦合分析得到的环向位移要大于不考虑裂隙内水压分析得到的结果。  相似文献   

9.
余天堂 《岩土力学》2007,28(Z1):305-310
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其原理是在裂尖附近用一些奇异函数和沿裂纹面用阶跃函数加强传统有限元的基,以考虑跨过裂纹的位移场的不连续,该加强策略允许计算网格独立于不连续体几何。讨论了扩展有限元法的一些数值方面,主要包括:水平集法确定界面和加强节点与加强方式、裂尖加强范围的选择、J积分区域的确定和积分方案等。  相似文献   

10.
杨石扣  张继勋  任旭华 《岩土力学》2018,39(Z1):488-494
基于数值流形法(NMM)进行三维裂纹扩展分析研究,编写了相应的C++程序。充分发挥数值流形法在非连续变形分析领域的优势,不需要扩展有限元中的水平集和阶跃函数概念,应用数值流形法计算结果分析裂纹尖端线的破坏状况,对已有的非局部求迹法和三角形推进法进行简化和扩展,提出一种简化算法确定最终的新裂纹扩展面。应用简化算法对水平钱币型裂纹扩展问题进行数值模拟和对比分析研究。计算结果表明,基于数值流形法的三维裂纹扩展算法是可行的,采用简化处理之后,克服原有算法的不足,大大减少了新裂纹尖端线和新裂纹面的数目,降低了计算网格对新裂纹面的影响,提高了计算效率和扩大了应用范围。  相似文献   

11.
郑安兴  罗先启 《岩土力学》2018,39(9):3461-3468
危岩是三峡库区典型的地质灾害类型之一,而主控结构面受荷断裂扩展是危岩发育成灾的关键核心。将危岩主控结构面类比为宏观裂纹,利用扩展有限元法在模拟裂纹扩展方面的优势,基于考虑裂纹面水压力作用的虚功原理推导出了采用扩展有限元法分析水力劈裂问题的控制方程,给出了危岩主控结构面水力劈裂问题的扩展有限元实现方法,对重庆万州太白岩危岩主控结构面的水力劈裂进行了数值模拟分析。计算结果表明:暴雨是威胁危岩稳定性的最敏感因素,随着裂隙水压力上升,裂端拉应力会急剧升高,危岩的稳定性降低;I型裂纹扩展是危岩主要的结构面扩展形式,结构面一旦发生开裂,将处于非稳定扩展状态。  相似文献   

12.
岩体内部赋存的裂隙很多表现为折线型,为探究这类岩体的断裂机制,制备含折线型裂隙砂岩试件并对其进行单轴压缩试验。采用数字图像相关(DIC)方法计算加载过程中的变形场演化,根据新生裂纹两侧的位移差异识别裂纹类型;运用扩展有限元法(XFEM)模拟断裂过程,根据应力分布特征解释翼型裂纹起裂与扩展机制。DIC计算结果表明,新生裂纹处出现应变局部化带,裂纹两侧发生相对分离;含直线型和折线型裂隙砂岩试件的翼型裂纹分别萌生于预制裂隙端部以及折角处,这是因为裂隙几何形态会改变拉应力集中位置;含折线型裂隙砂岩试件的起裂应力小于含直线型裂隙砂岩试件,这是因为相同加载条件下前者的最大拉应力值更大;这2类试件的裂纹扩展均是由于裂纹尖端集中的拉应力引起的,裂纹依然呈张开状态;裂隙几何形态未改变试件的最终破坏模式,均表现为对角剪切破坏。  相似文献   

13.
模拟裂纹传播的新方法——无网格伽辽金法   总被引:5,自引:1,他引:5  
无网格伽辽金法采用移动最小二乘法构造位移函数,由于它脱离了单元的概念,因此特别适合岩体介质裂纹的传播分析,计算了拉伸荷载作用下裂纹尖端的应力集中,对J积分围线作了进一步的讨论;运用无网格伽辽金法模拟岩体介质中的不连续面,计算了压剪复合型裂纹的应力强度因子,运用不同的断裂准则对裂纹的传播进行了分析模型,数值结果表明该方法与实验结果符合得很好。  相似文献   

14.
王伟  付豪  邢林啸  柴波  刘波  施星宇 《地球科学》2021,46(10):3509-3519
水力压裂作为一种主要的地热能开采手段,其压裂效果除与岩体基本物理力学性质有关外,还与裂隙分布、地应力状态、压裂工程参数等密切相关.为了探究以上因素对水力压裂过程中裂缝扩展行为的影响,以冀中坳陷碳酸盐岩储层岩体为研究对象,基于扩展有限元法,建立裂缝扩展流固耦合模型,分析了水平应力差、射孔方位角、注入液排量和压裂液黏度等参数对裂缝扩展行为的影响.结果表明:单裂缝扩展时,射孔方位角越小、注入量越大、越有利于裂缝扩展;双裂缝扩展时,水平应力差增大,裂缝偏转程度变小;水力裂缝与天然裂缝相交时,较小水平应力差有利于天然裂缝开启.   相似文献   

15.
随着扩展有限元理论的深入研究,利用扩展有限元方法模拟水力压裂具有了一定的可操作性。相比于常规有限元方法,XFEM方法具有计算结果精度高和计算量小的优点。但是,如何模拟射孔孔眼、如何模拟流体与岩石相互作用以及分析水力裂缝的扩展规律仍然是难题。以研究水力压裂裂缝扩展规律为目的,建立了岩石多孔介质应力平衡方程、流体渗流连续性方程和边界条件。通过有限元离散化方法对耦合方程矩阵进行处理。通过富集函数定义初始裂缝(射孔孔眼),选择最大主应力及损伤变量D分别作为裂缝起裂和扩展判定准则,利用水平集方法模拟水力裂缝扩展过程。数值模拟结果显示:增加射孔方位角、压裂液排量和减小水平地应力差,起裂压力上升;黏度对起裂压力无明显影响。增加射孔方位角、压裂液排量、黏度和减小水平地应力差值有助于裂缝宽度的增加。增加水平地应力差值、压裂液排量和减小射孔方位角以及压裂液黏度有助于裂缝长度增加,反之亦然。基于ABAQUS的水力裂缝扩展有限元法可对不同井型和诸多储层物性参数及压裂施工参数进行分析,且裂缝形态逼真,裂缝面凹凸程度清晰,结果准确。此研究可作为一种简便有效研究水力压裂裂缝扩展规律的方法为油田水力压裂设计与施工提供参考与依据。  相似文献   

16.
Rock strengths are directly influenced by the open or closed flaws widely distributed in rock masses. Extensive studies have been conducted on the propagations of open flaws in rocks. However, few concerns are paid on the propagation of closed flaws, the influence of the surface friction on the initiation and propagation of closed flaws should be investigated systematically. In present article, the crack initiation and propagation in rock like material subjected to compressive loads have been investigated. The effects of crack surface friction on crack initiation and propagation have been quantified with the help from extended finite element method which is efficient and accurate. Based on the analysis on stress distribution and propagation patterns, following results are obtained: Firstly, minor effects are exerted by crack surface friction on the stress distribution around the flaws when the flaws inclination angle is 30° and 45°. However, as the inclination angle increases to 60°, the effects are much more significant. Secondly, as the inclination angle ranges from 30° to 60°, the most favorable angle for crack propagation is 45°. Thirdly, the initiation location and angle of the wing cracks will not be influenced by the frictions. However, the propagation length will be greatly influenced by the friction and the inclination angle.  相似文献   

17.
梁东旭  张农  荣浩宇 《岩土力学》2023,(4):1217-1229
对预制交叉裂隙岩石试样进行裂纹扩展试验,研究了裂纹萌生、扩展、聚合过程,分析了主裂隙和轴向载荷夹角及主次裂隙夹角对裂纹起裂应力和聚合应力的影响,并用混合有限元-离散元程序,即图形处理器并行化的3D Y-HFDEM代码对试验进行了仿真计算,实现了岩石破坏从连续介质向非连续介质的过渡,对裂纹类及损伤破坏模式进行了识别,捕捉到了试验中难以发现的现象。研究表明:随主裂隙与轴向载荷夹角增加,裂纹聚合区的拉伸裂纹数量增加;裂纹起裂和聚合应力与主裂隙与轴向载荷夹角成正比;主次裂隙夹角增加,岩石的破坏模式由拉伸破坏转为剪切破坏,交叉裂隙加剧岩石破碎程度;主裂隙尖端萌生扩展的拉伸-剪切混合裂缝引起的破坏在岩石破坏中占主导地位,是导致岩体失去承载能力的主控裂纹;混合有限元-离散元仿真软件GPGPU并行化的3D Y-HFDEM IDE在岩石裂纹扩展研究中具有优势,可以捕捉实验室难以发现的损伤断裂类型,可以作为岩石裂纹扩展研究的有力工具。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号