基于M积分的脆性材料微缺陷等效损伤面积/体积表征

朱文洁; 吕俊男; 李群

doi:10.6052/0459-1879-17-378

基于M积分的脆性材料微缺陷等效损伤面积/体积表征

doi: 10.6052/0459-1879-17-378

西安交通大学航天航空学院机械结构强度与振动国家重点实验室,西安 710049

基金项目: 国家自然科学基金面上项目(11472205, 11772245), 高等学校学科创新引智计划(B18040)和中央高校基本科研业务费资助.

详细信息

作者简介:
null

作者简介：李群,教授,博士生导师,主要从事断裂与损伤力学研究. E-mail: qunli@mail.xjtu.edu.cn

中图分类号: O346.1;
计量
- 文章访问数: 1080
- HTML全文浏览量: 165
- PDF下载量: 420
- 被引次数: 0
出版历程
- 收稿日期: 2017-11-14
- 刊出日期: 2018-03-18

THE CALIBRATION OF MICRODEFECTS INDUCED EQUIVALENT DAMAGE AREA/VOLUME OF BRITTLE MATERIALS BY USING THE M-INTEGRAL

State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi’an Jiaotong University, Xi’an 710049, China

摘要

摘要: 随着脆性材料在工程中的广泛应用,对脆性材料中微缺陷进行统一的损伤水平标定,具有重要的科学研究和工程应用价值. 本研究提出一种基于M积分的材料等效损伤面积/体积标定方法, 以具有相同M积分值的圆孔面积或球孔体积来标定复杂微缺陷构型的损伤水平,从而实现不同类型微缺陷真实损伤水平的统一表征. 首先,基于Lagrangian能量密度函数推导了M积分定义式,并简述了M积分的物理意义,基于域积分方法实现二维/三维M积分的数值计算. 随后,提出了基于M积分的材料缺陷等效损伤面积/体积标定方法,以圆孔面积/球孔体积来标定复杂微缺陷材料系统的损伤水平. 最后, 针对单轴拉伸载荷作用下的二维/三维脆性体含不同缺陷构型,具体计算了椭圆孔、裂纹以及双裂纹、双孔洞、裂纹和孔洞干涉等复杂缺陷构型情况下的等效损伤面积/体积,并详细分析了缺陷之间的干涉效应及影响因素. 本研究旨在基于M积分等效方法量化脆性材料中各类微缺陷造成的损伤程度,实现损伤等级标定,有益于工程材料及结构的损伤容限设计及完整性评估.
- M积分 /
- 等效损伤面积 /
- 等效损伤体积 /
- 干涉效应
Abstract: In view of the integrity, reliability and functionality of brittle materials are substantially limited by the existence of microdefects, the calibration of materials’ damage level is of great scientific value and underlying engineering applications. An unified method of evaluating the microdefects induced equivalent damage area/volume is proposed in present study by the aid of M-integral. The damage area/volume induced by underlying multiple microdefects is assumed as equivalent to the area/volume of an individual circular/spherical void while the values of M-integral are equal for the both cases. Firstly, the analytical expression of M-integral is deduced by using the Lagrangian energy density function, the corresponding physical meaning is briefly elucidated. The domain integral method is applied to numerically calculating the M-integral for both two-dimensional (2D) and three-dimensional (3D) cases. Subsequently, the damage calibration process of arbitrary dispersed microdefects is given, the corresponding equivalent damage area for 2D defects and volume for 3D defects are defined. Finally, the elastic 2D plane and 3D body under uniaxial tensile loading condition is simulated, within which a series of different defect configurations are considered, including the singular defect (void, crack and ellipse) and the dual-defects (void-void, crack-crack, void-crack). Corresponding equivalent damage area or volume are calculated, the inherent “interactive effects” and influence factors are elucidated detailedly and quantitatively. Through the proposed damage calibration method in this study, we can estimate the damage level of any microdefects within brittle solids, the calibration process is simple and convenient, which will be beneficial to the damage tolerance design and integrity assessment of engineering structures.
- microdefects /
- M-integral /
- equivalent damage area /
- equivalent damage volume

HTML全文

参考文献(1)

[1]	Shao JF, Rudnicki JW. A microcrack-based continuous damage model for brittle geomaterials. Mechanics of Materials, 2000, 177-180(10): 607-619
[2]	嵇醒. 断裂力学判据的评述. 力学学报, 2016, 48(4): 741-753
[2]	(Ji Xing.A critical review on criteria of fracture mechanics.Chinese Journal of Theoretical and Applied Mechanics, 2016, 48(4): 741-753 (in Chinese))
[3]	Chen YH.M-integral analysis for two-dimensional solids with strongly interacting microcracks. Part 1: In an infinite brittle solid.International Journal of Solids and Structures, 2001, 38(18): 3193-3212
[4]	Chen YH.M-integral analysis for two-dimensional solids with strongly interacting microcracks. Part 2: In the brittle phase of an infinite metal/ceramic bimaterial.International Journal of Solids and Structures, 2001, 38(18): 3213-3232
[5]	Ju JW, Chen TM.Effective elastic moduli of two-dimensional brittle solids with interacting microcracks, part II: Evolutionary damage models.Journal of Applied Mechanics, 1994, 61(2): 358-366
[6]	Hu Y, Chen Y.The M -integral description for a brittle plane strip with two holes before and after coalescence.Acta Mechanica, 2009, 204(1-2):109-123
[7]	Yu MH.Advances in strength theories for materials under complex stress state in the 20th Century.Advances in Mechanics, 2004, 34(4): 529-560
[8]	GB/T1979-2001. 结构钢低倍组织缺陷评级图
[8]	(GB/T1979-2001. Standard diagrams for macrostructure and defect of structural steels (in Chinese))
[9]	Pak YE, Mishra D, Yoo SH.Energy release rates for various defects under different loading conditions. Journal of Mechanical Science and Technology, 2012, 26(11): 3549-3554
[10]	徐真军, 杨红兵.关于条形缺陷的射线照相评级探讨.无损探伤, 2010, 34(2): 28-29
[10]	(Xu Zhenjun, Yang Honbing.Analysis of the strip defect rating by using the radiography technology.Nondestructive Testing, 2010, 34(2): 28-29 (in Chinese))
[11]	汪孝欢, 牛显明, 张鹏林等. 焊缝数字图像圆形缺陷的自动评级方法. 无损检测, 2012, 34(9): 42-45
[11]	(Wang Xiaohuan, Niu Xianming, Zhang Penglin, et al.Research on automatic rating method for digital image circular weld defects.Nondestructive Testing, 2012, 34(9): 42-45 (in Chinese))
[12]	Eshelby JD, Frank FC, Nabarro FRN.XLI. The equilibrium of linear arrays of dislocations.Philosophical Magazine, 1951, 42(327): 351-364
[13]	Baxevanakis KP, Giannakopoulos AE.Finite element analysis of discrete edge dislocations: Configurational forces and conserved integrals.International Journal of Solids & Structures, 2015, 62: 52-65
[14]	Judt PO, Ricoeur A.A new application of M- and L-integrals for the numerical loading analysis of two interacting cracks.Journal of Applied Mathematics and Mechanics, 2015, 96(1): 24-36
[15]	Özenç K, Chinaryan G, Kaliske M.A configurational force approach to model the branching phenomenon in dynamic brittle fracture.Engineering Fracture Mechanics, 2016, 157: 26-42
[16]	古斌, 郭宇立, 李群. 基于构型力断裂准则的裂纹与夹杂干涉问题. 力学学报, 2017, 49(6): 1312-1321
[16]	(Gu Bin, Guo Yuli, Li Qun.Crack interacting with an individual inclusion by the fracture criterion of configurational force.Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(6): 1312-1321 (in Chinese))
[17]	Guo YL, Li Q.Material configurational forces applied to mixed mode crack propagation.Theoretical and Applied Fracture Mechanics, 2017, 89: 145-157
[18]	Tracy J, Waas A, Daly S.Experimental assessment of toughness in ceramic matrix composites using the J -integral with digital image correlation part I: Methodology and validation.Journal of Materials Science, 2015, 50(13): 4646-4658
[19]	Paxevanakis BK, Giannakopoulos AE.Finite element analysis of discrete edge dislocations: Configurational forces and conserved integrals.International Journal of Solids and Structures, 2015, 62(52-65): 52-65
[20]	李群. 材料构型力学及其在复杂缺陷系统中的应用. 力学学报, 2015, 47(2): 197-214
[20]	(Li Qun.Material configurational mechanics with application to complex defects.Acta Mechanica Sinica, 2015, 47(2): 197-214 (in Chinese))
[21]	Hu YF, Chen YH.The M-integral description for a brittle plane strip with two cracks before and after coalescence.Journal of Applied Mechanics, 2009, 76(6): 061017
[22]	于宁宇, 李群. M积分与夹杂/缺陷弹性模量的显式关系. 力学学报, 2014, 46(1): 87-93
[22]	(Yu Ningyu, Li Qun.The explicit relation between the m-integral and the elastic moduli of inclusion/damages.Acta Mechanica Sinica, 2014, 46(1): 87-93 (in Chinese))
[23]	Yu NY, Li Q.Failure theory via the concept of material configurational forces associated with the M-integral.International Journal of Solids & Structures, 2013, 50(25-26): 4320-4332
[24]	Li Q, Guo YL, Hou JL, et al.The M- integral based failure description on elasto-plastic materials with defects under biaxial loading.Mechanics of Materials, 2017, 112: 163-171
[25]	Goldstein RV, Shifrin EI, Shushpannikov PS.Application of invariant integrals to the problems of defect identification.International Journal of Fracture, 2007, 147(1-4): 45-54
[26]	Gommerstadt BY.The J and M integrals for a cylindrical cavity in a time-harmonic wave field.International Journal of Engineering Science, 2014, 83: 76-84
[27]	Meyer CR, Hutchinson JW, Rice JR. The path-independent M integral implies the creep closure of englacial and subglacial channels. Journal of Applied Mechanics, 2016, 84(1):011006(1-9)
[28]	Ma LF, Chen YH, Liu CS.On the relation between the M-integral and the change of the total potential energy in damaged brittle solids.Acta mechanica, 2001, 150(1): 79-85
[29]	Li FZ, Shih CF, Needleman A.A comparison of methods for calculating energy release rates.Engineering Fracture Mechanics, 1985, 21(2): 405-421
[30]	Moran B, Shih CF.A general treatment of crack tip contour integrals.International Journal of Fracture, 1987, 35(4): 295-310
[31]	Nikishkov GP, Atluri SN.An equivalent domain integral method for computing crack-tip integral parameters in non-elastic, thermo-mechanical fracture.Engineering Fracture Mechanics, 1987, 26(6): 851-867
[32]	Shivakumar KN, Raju IS.An equivalent domain integral method for three-dimensional mixed-mode fracture problems.Engineering Fracture Mechanics, 1992, 42(6): 935-959
[33]	Chen YZ, Lee KY.Analysis of the M-integral in plane elasticity.Journal of Applied Mechanics, 2004, 71(4): 572-574
[34]	Seo SY, Mishra D, Park CY, et al.Energy release rates for a misfitted spherical inclusion under far-field mechanical and uniform thermal loads.European Journal of Mechanics, 2015, 49: 169-182