目的 针对多相夹层材料中界面两侧性能突变导致基于均质材料假设的传统J积分路径相关性显著、界面贡献难以分离的问题,建立考虑相边界效应的多相材料裂纹断裂参数评价方法。方法 在传统域积分理论基础上,引入材料非均匀性修正项,构建考虑相边界效应的修正J积分计算框架,并基于Abaqus二次开发实现裂纹前沿自动识别、积分域构建、相界面判定及断裂参数自动化提取。通过统一求解传统域积分项与相界面修正项,实现多相夹层材料裂纹尖端能量释放率的准确计算。结果 对含材料界面的双相材料紧凑拉伸试件开展有限元仿真分析,所提出的修正J积分方法能够有效降低材料界面对积分路径选择的敏感性,修正后J积分在不同积分轮廓下的计算结果离散程度较低,最大相对偏差控制在7%以内。同时,分析了弹性模量失配条件下界面区域应力场演化特征及其对裂纹扩展驱动力的影响规律。结论 所建立的考虑相边界效应的修正J积分方法能够有效提升多相夹层材料裂纹断裂参数评价的可靠性,为涂层结构、层合复合材料及异质夹层构件的损伤评估和结构安全设计提供有效的数值方法和计算工具。
Abstract
For multiphase sandwich materials, abrupt changes in material properties across interfaces can lead to pronounced path dependence of the conventional J-integral formulated under the assumption of material homogeneity, making it difficult to isolate the contribution of material interfaces. Based on the conventional domain-integral formulation, a material-inhomogeneity correction term was introduced to establish a modified J-integral framework accounting for phase-boundary effects. An automated computational procedure was further implemented through Abaqus secondary development, enabling automatic identification of the crack front, construction of integration domains, detection of phase interfaces, and extraction of fracture parameters. By consistently evaluating the conventional domain-integral term together with the phase-interface correction term, the proposed framework enabled accurate calculation of the energy release rate at the crack tip in multiphase sandwich materials. Finite element simulations were performed on compact tension specimens of bimaterial systems containing material interfaces. The results demonstrated that the proposed modified J-integral method effectively reduced the sensitivity of the calculated fracture parameters to the selection of integration contours. The modified J-integral exhibited good contour stability, with a maximum relative deviation of less than 7% among different integration contours. In addition, the evolution of the stress field in the interfacial region under elastic-modulus mismatch and its influence on the crack-driving force were systematically investigated. The proposed modified J-integral method accounting for phase-boundary effects improves the reliability of fracture-parameter evaluation for multiphase sandwich materials and provides an effective numerical methodology and computational tool for damage assessment and structural integrity design of coating systems, laminated composites, and heterogeneous sandwich structures.
关键词
多相夹层材料 /
相边界效应 /
修正J积分 /
Abaqus二次开发 /
裂纹驱动力 /
断裂参数
Key words
multiphase sandwich materials /
phase boundary effect /
modified J-integral /
Abaqus secondary development /
crack driving force /
fracture parameter
{{custom_sec.title}}
{{custom_sec.title}}
{{custom_sec.content}}
参考文献
[1] HUTCHINSON J W, SUO Z.Mixed Mode Cracking in Layered Materials[J]. Advances in Applied Mechanics, 1992, 29: 63-191.
[2] BOGGARAPU V, GUJJALA R, OJHA S, et al.State of the Art in Functionally Graded Materials[J]. Composite Structures, 2021, 262: 113596.
[3] ERDOGAN F.Fracture Mechanics of Functionally Graded Materials[J]. Composites Engineering, 1995, 5(7): 753-770.
[4] ESHELBY J D.The Force on an Elastic Singularity[J]. Philosophical Transactions of the Royal Society of London Series A, Mathematical and Physical Sciences, 1951, 244(877): 87-112.
[5] RICE J R.A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks[J]. Journal of Applied Mechanics, 1968, 35(2): 379-386.
[6] LIEBOWITZ H.Fracture: An Advanced Treatise, Vol. 2: Mathematical Fundamentals[M]. New York: Academic Press, 1968: 191-311.
[7] ANDERSON T L.Fracture Mechanics: Fundamentals and Applications[M]. 4th ed. Boca Raton: CRC Press, 2017: 1-2.
[8] HONEIN T, HERRMANN G.Conservation Laws in Nonhomogeneous Plane Elastostatics[J]. Journal of the Mechanics and Physics of Solids, 1997, 45(5): 789-805.
[9] WEICHERT D, SCHULZ M.J-Integral Concept for Multi-Phase Materials[J]. Computational Materials Science, 1993, 1(3): 241-248.
[10] HADDI A, WEICHERT D.On the Computation of the J-Integral for Three-Dimensional Geometries in Inhomogeneous Materials[J]. Computational Materials Science, 1996, 5(1/2/3): 143-150.
[11] GU P, DAO M, ASARO R J.A Simplified Method for Calculating the Crack-Tip Field of Functionally Graded Materials Using the Domain Integral[J]. Journal of Applied Mechanics, 1999, 66(1): 101-108.
[12] CHEN J, WU L Z, DU S Y.A Modified J Integral for Functionally Graded Materials[J]. Mechanics Research Communications, 2000, 27(3): 301-306.
[13] SIMHA N K, FISCHER F D, KOLEDNIK O, et al.Inhomogeneity Effects on the Crack Driving Force in Elastic and Elastic-Plastic Materials[J]. Journal of the Mechanics and Physics of Solids, 2003, 51(1): 209-240.
[14] SIMHA N K, FISCHER F D, KOLEDNIK O, et al.Crack Tip Shielding or Anti-Shielding Due to Smooth and Discontinuous Material Inhomogeneities[J]. International Journal of Fracture, 2005, 135(1): 73-93.
[15] PARKS D M.A Stiffness Derivative Finite Element Technique for Determination of Crack Tip Stress Intensity Factors[J]. International Journal of Fracture, 1974, 10(4): 487-502.
[16] HELLEN T K.On the Method of Virtual Crack Extensions[J]. International Journal for Numerical Methods in Engineering, 1975, 9(1): 187-207.
[17] MORAN B, SHIH C F.Crack Tip and Associated Domain Integrals from Momentum and Energy Balance[J]. Engineering Fracture Mechanics, 1987, 27(6): 615-642.
[18] NIKISHKOV G P, ATLURI S N.Calculation of Fracture Mechanics Parameters for an Arbitrary Three-Dimensional Crack, by the ‘Equivalent Domain Integral’ Method[J]. International Journal for Numerical Methods in Engineering, 1987, 24(9): 1801-1821.
[19] KIM J H, PAULINO G H.Finite Element Evaluation of Mixed Mode Stress Intensity Factors in Functionally Graded Materials[J]. International Journal for Numerical Methods in Engineering, 2002, 53(8): 1903-1935.
[20] KIM J H, PAULINO G H.Mixed-Mode J-Integral Formulation and Implementation Using Graded Elements for Fracture Analysis of Nonhomogeneous Orthotropic Materials[J]. Mechanics of Materials, 2003, 35(1/2): 107-128.
[21] WALTERS M C, PAULINO G H, DODDS R H.Stress-Intensity Factors for Surface Cracks in Functionally Graded Materials under Mode-I Thermomechanical Loading[J]. International Journal of Solids and Structures, 2004, 41(3/4): 1081-1118.
[22] 白晓明, 熊夫睿, 谢海, 等. 基于交互作用积分的异种金属焊接区域裂纹扩展模拟方法[J]. 原子能科学技术, 2024, 58(5): 1109-1118.
BAI X M, XIONG F R, XIE H, et al.Numerical Simulation of Crack Propagation in Dissimilar Metal Welded Region Based on Interaction Integral[J]. Atomic Energy Science and Technology, 2024, 58(5): 1109-1118.
[23] 王必宁, 贾普荣, 矫桂琼. 用三维有限元模型计算复合材料粘接修补裂纹板的J积分[J]. 计算力学学报, 2005, 22(5): 603-607.
WANG B N, JIA P R, JIAO G Q.Computation of the J-Integral Value for Single-Edge Crack Sheet with Single-Side Bonded Composite Patch Using 3D Finite Element Model[J]. Chinese Journal of Computational Mechanics, 2005, 22(5): 603-607.
[24] BARHLI S M, MOSTAFAVI M, CINAR A F, et al.J-Integral Calculation by Finite Element Processing of Measured Full-Field Surface Displacements[J]. Experimental Mechanics, 2017, 57(6): 997-1009.
[25] 陈卫忠, 陈培帅, 王辉. 裂纹扩展过程的动态仿真技术[J]. 岩土力学, 2011, 32(S2): 573-579.
CHEN W Z, CHEN P S, WANG H.Dynamic Simulation Technology for Crack Extension Process[J]. Rock and Soil Mechanics, 2011, 32(S2): 573-579.
[26] QIU L X, YANG F, ZHANG Z N, et al.Evaluating the J-Integral of Mode I Crack in Gradient Nanocrystalline Metals with Residual Stress Induced by Surface Mechanical Treatment[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2022, 236(12): 6760-6772.
[27] LEE G B, PARK S H, JANG Y Y, et al.Development of Automatic Crack Growth Simulation Program Based on Finite Element Analysis[J]. Applied Sciences, 2022, 12(6): 3075.
[28] INTERNATIONAL A.Standard Test Method for Measurement of Fatigue Crack Growth Rates: ASTM E647- 24[S]. West Conshohocken: ASTM International, 2024.
基金
国家科技重大专项（2025ZD1600600）; 船舶CAE研发应用项目; 广东省基础与应用基础研究基金（2025A1515011471）
PDF(5196 KB)
