圆盘内裂纹扩展断裂试验及数值模拟研究
|
摘要
鉴于含裂纹的巴西圆盘试验多集中于二维或表面裂纹,而对含纯封闭内裂纹的巴西圆盘断裂研究较少的问题,开展了含90°三维内裂纹的树脂巴西圆盘试样单轴压缩试验,分析了裂纹扩展过程特征,并基于M积分及最大拉应力准则研究了KⅠ、KⅡ、KⅢ分布规律,通过断裂力学有限元软件FRANC3D实现内裂纹扩展全过程模拟。结果表明,含90°内裂纹试样在单轴压缩下裂纹破坏形态为竖直张拉裂纹;达到峰值荷载29.13kN时试样瞬间被劈裂为两半,内裂纹扩展速度达到700mm/s,对试样扩展、破坏形态及断口特征产生显著影响;通过数值模拟得出内裂纹尖端扩展路径,短轴端扩展速率最大,且数值模拟结果与试验相符。
In view of the fact that the Brazilian disc tests with cracks are mostly concentrated on two-dimensional or surface cracks,and the less research on Brazilian disc specimen with pure closed internal crack,the uniaxial compression test of Brazilian disc specimen with 90°three-dimensional internal crack was carried out.The characteristics of crack propagation process was analyzed.Based on M integral and maximum tensile stress criterion,the distribution law of KI,KII and KIIIwas studied.The whole process of internal crack propagation was simulated by FRANC3 D.The results show that the crack failure morphology of the crack specimen with 90°under uniaxial compression is vertical tensile crack;When the peak load is 29.13 kN,the specimen is instantaneously split into two halves,and the internal crack propagation speed reaches 700 mm/s,which has a significant impact on the sample expansion,failure mode and fracture characteristics;The internal crack tip extension path is obtained by numerical simulation,and the short axis end expansion rate is the largest;Furthermore,the numerical simulation results are consistent with the test.
In view of the fact that the Brazilian disc tests with cracks are mostly concentrated on two-dimensional or surface cracks,and the less research on Brazilian disc specimen with pure closed internal crack,the uniaxial compression test of Brazilian disc specimen with 90°three-dimensional internal crack was carried out.The characteristics of crack propagation process was analyzed.Based on M integral and maximum tensile stress criterion,the distribution law of KI,KII and KIIIwas studied.The whole process of internal crack propagation was simulated by FRANC3 D.The results show that the crack failure morphology of the crack specimen with 90°under uniaxial compression is vertical tensile crack;When the peak load is 29.13 kN,the specimen is instantaneously split into two halves,and the internal crack propagation speed reaches 700 mm/s,which has a significant impact on the sample expansion,failure mode and fracture characteristics;The internal crack tip extension path is obtained by numerical simulation,and the short axis end expansion rate is the largest;Furthermore,the numerical simulation results are consistent with the test.
引文
[1] Hondros G.The Evaluation of Poisson’s Ratio and Modulus of Materials of a Low Tensile Resistence by the Brazilian(Indirect Tensile)Test with Particular Reference to Concrete[J].Australian J.Appl.Sci.,1959,10(3):243-268.
[2]黄耀光,王连国,陈家瑞,等.平台巴西劈裂试验确定岩石抗拉强度的理论分析[J].岩土力学,2015,36(3):739-748.
[3]许媛,戴峰,徐奴文,等.人字形切槽巴西圆盘岩石试样复合型断裂渐进过程数值模拟研究[J].岩土工程学报,2015,37(12):2 189-2 197.
[4]张志强,关宝树,郑道访.用直切槽巴西圆盘试样柔度法确定岩石的弹性模量[J].岩石力学与工程学报,1998,17(4):24-30.
[5]黄彦华,杨圣奇.孔槽式圆盘破坏特性与裂纹扩展机制颗粒流分析[J].岩土力学,2014,35(8):2 269-2 277.
[6] Brian Lawn.脆性固体断裂力学(第二版)[M].龚江宏译.北京:高等教育出版社,2010.
[7] BANKS-SILLS LESLIE,WAWRZYNEK PAUL A,CARTER BRUCE,et al.Methods for Calculating Stress Intensity Factors in Anisotropic Materials:PartⅡ—Arbitrary Geometry[J].Engineering Fracture Mechanics,2007,74(8):1 293-1 307.
[8] LI F Z,SHIH C F,NEEDLEMAN A.A Comparison of Methods for Calculating Energy Release Rates[J].Engfractmech,1985,21(2):405-421.
[2]黄耀光,王连国,陈家瑞,等.平台巴西劈裂试验确定岩石抗拉强度的理论分析[J].岩土力学,2015,36(3):739-748.
[3]许媛,戴峰,徐奴文,等.人字形切槽巴西圆盘岩石试样复合型断裂渐进过程数值模拟研究[J].岩土工程学报,2015,37(12):2 189-2 197.
[4]张志强,关宝树,郑道访.用直切槽巴西圆盘试样柔度法确定岩石的弹性模量[J].岩石力学与工程学报,1998,17(4):24-30.
[5]黄彦华,杨圣奇.孔槽式圆盘破坏特性与裂纹扩展机制颗粒流分析[J].岩土力学,2014,35(8):2 269-2 277.
[6] Brian Lawn.脆性固体断裂力学(第二版)[M].龚江宏译.北京:高等教育出版社,2010.
[7] BANKS-SILLS LESLIE,WAWRZYNEK PAUL A,CARTER BRUCE,et al.Methods for Calculating Stress Intensity Factors in Anisotropic Materials:PartⅡ—Arbitrary Geometry[J].Engineering Fracture Mechanics,2007,74(8):1 293-1 307.
[8] LI F Z,SHIH C F,NEEDLEMAN A.A Comparison of Methods for Calculating Energy Release Rates[J].Engfractmech,1985,21(2):405-421.