{"title":"Deformation and fracture simulation of explosive charges under low velocity impact by numerical manifold method","authors":"Rui Yue , Pengwan Chen , Youjun Ning , Ge Kang","doi":"10.1016/j.engfracmech.2025.111209","DOIUrl":"10.1016/j.engfracmech.2025.111209","url":null,"abstract":"<div><div>Numerical Manifold Method (NMM) integrates continuous and discontinuous features with mathematical and physical cover systems (MC and PC), enabling it to solve problems related to crack initiation and propagation. This study aims to analyze the fracture characteristics of explosive charges subjected to low-velocity impact, elucidating the underlying patterns in their mechanical behavior. The existing NMM framework has been innovatively extended to enable the solution of large deformations, multi-crack fields, and crack interactions. Based on the Mohr-Coulomb criterion, the novel integration of PC subdivision, element deletion, and data transmission algorithms within the NMM framework has been developed and implemented. The division technique ensures multiple tension and shear cracks can be captured, while the deletion algorithm addresses program errors caused by the distortion of the fine elements during the impact process. A simulation of the standard Steven test was conducted to evaluate stress distribution in explosive charges under low-velocity impact; consequently, the results of our NMM framework are in agreement with the LS-DYNA model. Furthermore, systematic simulations of the Steven test with varying specimen sizes and impact velocities are analyzed. The simulation results indicate that as specimen thickness increases, fragmentation intensity significantly decreases; whereas, as specimen diameter increases, the confining pressure in the impact area decreases, resulting in more severe fragmentation. The enhanced NMM developed in this study effectively simulates the dynamic deformation and fracture behavior of explosive charges.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111209"},"PeriodicalIF":4.7,"publicationDate":"2025-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144068089","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Tom Petit , Pierrick François , Bedi Aydin Baykal , Zoltán Bézi , Alexander Bakaev , Giovanni Bonny , Tomasz Brynk , Pierre-Yves Corbel , Gintautas Dundulis , Réka Erdei , Petr Gál , Vitalis Leisis , Meng Li , Timo Metzler , Evaldas Narvydas , Frédéric Péralès , Miroslav Posta , Tomoki Shinko , Katarina Siskova , Philippe Spätig , Masato Yamamoto
{"title":"Numerical investigations on size effects and mini-CT applications for Master Curve determination","authors":"Tom Petit , Pierrick François , Bedi Aydin Baykal , Zoltán Bézi , Alexander Bakaev , Giovanni Bonny , Tomasz Brynk , Pierre-Yves Corbel , Gintautas Dundulis , Réka Erdei , Petr Gál , Vitalis Leisis , Meng Li , Timo Metzler , Evaldas Narvydas , Frédéric Péralès , Miroslav Posta , Tomoki Shinko , Katarina Siskova , Philippe Spätig , Masato Yamamoto","doi":"10.1016/j.engfracmech.2025.111227","DOIUrl":"10.1016/j.engfracmech.2025.111227","url":null,"abstract":"<div><div>This paper presents a comprehensive numerical study on size effects and the application of mini-compact tension (mini-CT, or MCT) specimens within the Master Curve methodology for determining fracture toughness, particularly in irradiated nuclear materials. This study, conducted as part of the FRACTESUS project, involves collaboration among several European laboratories to estimate the efficacy of mini-CT specimens through extensive finite element modelling (FEM) and inter-laboratory simulations. The research addresses critical factors, including the consistency of FEM codes, the impact of crack length on displacement conversion factors, and the application of the Beremin model for brittle fracture analysis. The good consistency between the results obtained by the different laboratories validates the numerical approach. The comparison of the macroscopic and local mechanical fields between 1T-CT and MCT specimens highlights the in- and out-plane loss of constraint and the deterioration in the plane strain state in MCT, resulting in the apparent fracture toughness shift on the measured <span><math><msub><mi>T</mi><mn>0</mn></msub></math></span> value. The numerical analysis of this shift using the Beremin model show that: (i) numerical simulations can accurately replicate experimental results obtained with MCT specimens, (ii)<!--> <!-->a size effect is observed on the Beremin fracture parameters, (iii) the use of cross-parameter sets between geometries does not yield satisfactory results, and (iv) the <span><math><msub><mi>T</mi><mn>0</mn></msub></math></span> value, using the same parameters for both specimen geometries, is lower is the MCT compared to 1T-CT suggesting that MCT can lead to non-conservative results with respect to 1T-CT.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"325 ","pages":"Article 111227"},"PeriodicalIF":4.7,"publicationDate":"2025-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144314085","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ping Wang , Meng Li , Zhixun Wen , Hao Cheng , Yuanmin Tu , Pengfei He
{"title":"Machine learning-enabled creep life prediction model for nickel-based single crystal superalloys with consideration of structural shape effects","authors":"Ping Wang , Meng Li , Zhixun Wen , Hao Cheng , Yuanmin Tu , Pengfei He","doi":"10.1016/j.engfracmech.2025.111230","DOIUrl":"10.1016/j.engfracmech.2025.111230","url":null,"abstract":"<div><div>Significant discrepancies have been observed in the creep life of SX alloy specimens with diverse structural configurations, manifesting an evident “shape effect”. Therefore, it is of great significance to establish a unified creep life prediction model capable of accommodating a broad spectrum of structural forms. In this study, based on the data sets of creep life of specimens with different structural forms, RBF-ANN, GABP-ANN and XGBoost machine learning paradigms were used to predict the creep life of SX alloy specimens, and the prediction effects of the three models were evaluated. The analysis of mean absolute error and determination coefficient shows that the RBF-ANN model has better fitting performance and generalization ability. The contribution of various shape parameters to creep life is ranked as gauge length > cross-sectional parameters > critical cross-sectional area > SCF > perimeter. Furthermore, verified by conducting creep fracture tests under different structures, the results show that the prediction accuracy of RBF-ANN model can be controlled within the range of the two-times scatter band, which shows the effectiveness of the established model. This method provides a new idea for the life evaluation of specimens with different structural shapes and has the potential to be extended to other structural forms and mechanical properties.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111230"},"PeriodicalIF":4.7,"publicationDate":"2025-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144071293","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Tianlong Ma , Qiaoling Zhang , Yongbin Ge , Wentao Ma
{"title":"A damage-driven adaptive radial point interpolation method for phase field model of brittle fracture","authors":"Tianlong Ma , Qiaoling Zhang , Yongbin Ge , Wentao Ma","doi":"10.1016/j.engfracmech.2025.111168","DOIUrl":"10.1016/j.engfracmech.2025.111168","url":null,"abstract":"<div><div>To enhance the computational accuracy and efficiency in the analysis of phase field model (PFM) for brittle fracture, a damage-driven adaptive Radial Point Interpolation Method (DARPIM) is proposed. We employ the Polyharmonic Spline (PHS) radial basis function augmented with linear polynomial basis to approximate the displacement field as well as phase field. The coupled non-linear system of these two fields is solved by the staggered iteration scheme. Motivated by the damage evolution principle of continuum damage mechanics, we design a novel refinement criterion to improve both computational efficiency and accuracy of RPIM. In this criterion, the smeared crack zone is divided into low-, medium-, and high-damage regions based on the phase field value, and the corresponding refinement factor in each region is taken as 0.5, 1, and 2, respectively. The significant advantage of this criterion is that it can not only automatically refine nodes by tracking the crack evolution process, but also continuously adjust node density from sparse to dense within the smeared crack zone. To evaluate the accuracy, efficiency and robustness of the DARPIM, several 2D and 3D examples are simulated. The results demonstrate that the proposed method can capture the true crack paths with considerably fewer nodes, while the crack propagation paths and load–displacement curves exhibit excellent agreement with experimental data and results reported in the published literature.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111168"},"PeriodicalIF":4.7,"publicationDate":"2025-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144068087","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Shanshan Qin , Ming Zhang , Dujian Zou , Tiejun Liu , Ye Li
{"title":"Mesoscale modeling of mechanical deterioration in sulfate-attacked concrete","authors":"Shanshan Qin , Ming Zhang , Dujian Zou , Tiejun Liu , Ye Li","doi":"10.1016/j.engfracmech.2025.111229","DOIUrl":"10.1016/j.engfracmech.2025.111229","url":null,"abstract":"<div><div>This study presents a mesoscale mechanical deterioration model to investigate the chemo-mechanical degradation of concrete under sulfate attack. The model introduces sulfate-induced volumetric expansion at the microscopic level and incorporates its macroscopic equivalent expansion strain into a mechanical damage framework. A two-dimensional polygonal random aggregate structure is employed to reflect the heterogeneous microstructure of concrete and simulate damage evolution under sulfate attack. Validation against published experimental data demonstrates the model’s accuracy in capturing expansion behavior, cracking patterns, and compressive strength degradation under sulfate exposure. Simulations reveal non-uniform damage initiation at aggregate corners and propagation along aggregate–mortar interfaces, ultimately leading to macrocracking and strength loss. A continuous decline in compressive strength with increasing exposure duration confirms the model’s predictive capability. The study underscores the critical role of concrete heterogeneity in influencing ion transport, damage localization, and failure mechanisms. By distinguishing between mortar and aggregate phases, the model reflects tortuosity and dilution effects on ion diffusion and reaction product accumulation. This mesoscale framework offers mechanistic insight into the coupled transport–mechanical processes driving sulfate-induced degradation. Despite simplifications such as the exclusion of the interfacial transition zone and post-cracking transport evolution, the model provides a foundation for future refinements and supports the durability assessment of concrete structures in aggressive environments.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111229"},"PeriodicalIF":4.7,"publicationDate":"2025-05-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143949021","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Tao Liao , Zihan Zhou , Jun Lai , Ping Wang , Xicheng Feng , Zhaoguang Zheng , Kai Wang , Jingmang Xu
{"title":"A probabilistic assessment method for fatigue reliability of bonded insulated rail joint under random rolling contact","authors":"Tao Liao , Zihan Zhou , Jun Lai , Ping Wang , Xicheng Feng , Zhaoguang Zheng , Kai Wang , Jingmang Xu","doi":"10.1016/j.engfracmech.2025.111226","DOIUrl":"10.1016/j.engfracmech.2025.111226","url":null,"abstract":"<div><div>Bonded Insulated Rail Joints (BIRJs) face significant challenges due to uncertainties in load conditions and dynamic wheel-rail interactions, which affect their fatigue reliability. This study proposes a fatigue reliability assessment method for transient rolling contact in BIRJs. To reduce computational cost, a surrogate model based on the Whale Optimization Algorithm (WOA) and Artificial Neural Networks (ANN) is used, incorporating a multi-error control (MEC) strategy to enhance prediction accuracy of fatigue damage. The WOA-MEC-ANN approach is essential as it efficiently maps multiple fatigue model responses and better captures the uncertainties in wheel-rail contact, while maintaining accuracy. Shear cracks, primarily located near the endpost, are identified as the main failure mode based on the J-S fatigue model. Sensitivity analysis reveals that lateral contact position and speed are the most influential factors on fatigue life, which also significantly impact wheel-rail contact behavior. These findings provide important insights for predicting fatigue life and crack initiation in BIRJs, offering valuable references for real-world engineering practices and operational maintenance strategies.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111226"},"PeriodicalIF":4.7,"publicationDate":"2025-05-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143942462","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Ductile fracture prediction of titanium foil based on shear-modified GTN damage model considering size effect","authors":"Xiaofeng Ding , Xinyuan Wei , Xuan Zhang , Yuhua Hou , Fuqiang Zhao","doi":"10.1016/j.engfracmech.2025.111200","DOIUrl":"10.1016/j.engfracmech.2025.111200","url":null,"abstract":"<div><div>Uniaxial tensile of 0.1 mm TA1 foils with different grain size were conducted to analyze size effect on tensile behaviors. As the grain size increases, the true stress increases while the elongation and fraction of twins decrease, and the fracture mode changes from ductile to quasi-cleavage and ultimately to cleavage fracture. Thus this paper established a microscale damage model considering the size effect, and introduced the J-GTN damage model to adapt to the complex stress state, and the parameters were calibrated through tensile, shear and compression tests. VUMAT subroutine was written and embedded in Abaqus software to simulate the fracture process of the sample under various stress states. The simulation results of cracks agree well with the experimental ones. This model can accurately predict the crack initiation and the final failure, which indicates that it is feasibility and accuracy in the prediction of titanium foil microforming.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111200"},"PeriodicalIF":4.7,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143937600","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Transition from large-scale fibre bridging to short-scale via a thin plain weave interleaf","authors":"Christopher Sutcu, Ali Aravand, Zafer Kazancı","doi":"10.1016/j.engfracmech.2025.111225","DOIUrl":"10.1016/j.engfracmech.2025.111225","url":null,"abstract":"<div><div>Hybridisation of fibre architectures in composite laminates offers a means to improve damage tolerance. This has been investigated in the area of Hybrid Unidirectional Woven Composite Laminates (HUWCL) under impact loading conditions. This study examines the Mode I delamination behaviour of hybrid unidirectional and woven interfaces in Carbon Fibre Reinforced Polymer (CFRP) composites by introducing a thin plain weave (PW) interleaf at the midplane of a unidirectional (UD) laminate. Double Cantilever Beam (DCB) testing was conducted to evaluate the effect of hybridisation on Mode I fracture toughness and fibre bridging laws for two different initial crack lengths. Results reveal that hybridisation with a thin (90 g/m<sup>2</sup>) PW interleaf significantly alters the delamination process by reducing large-scale fibre bridging and concentrating energy dissipation at the crack tip. This led to a substantial increase in initiation fracture toughness, with a 147 % improvement observed in PW/UD hybrid laminates compared to non-hybridised UD controls. The hybrid interfaces also exhibited a uniform R-curve. Additionally, bridging laws demonstrated that hybridisation shifts the fibre traction response, with PW interleafed laminates producing short-scale bridging rather than large-scale fibre bridging typically observed in UD laminates. While for the non-hybridised control laminates, initial crack length changes the bridging stress distribution, indicating that the bridging law is not a material property when considering initial crack length. These findings suggest that selective hybridisation with thin PW layers offers a promising strategy to enhance fracture toughness while mitigating the stochastic nature of large-scale fibre bridging generation and its impact on propagation fracture toughness.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111225"},"PeriodicalIF":4.7,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144071296","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"On the relevance of the equivalent crack length in the determination of the ERR in modes I and II","authors":"Faustino Mujika, Ainhoa Arrese","doi":"10.1016/j.engfracmech.2025.111224","DOIUrl":"10.1016/j.engfracmech.2025.111224","url":null,"abstract":"<div><div>The equivalent crack length in a fracture test is determined by the compliance, being known the elastic properties of the specimen. In this study, asymmetric double cantilever beam in pure mode I and asymmetric end notched flexure in pure mode II are analysed numerically. In a first step, three bi-material cases of small fracture process zone are analyzed to compare the crack length obtained by an analytical approach with the input value of the numerical model. Then, an aluminium-composite joint with a significant fracture process zone has been analysed. The energy release rate has been determined by two values of the crack length: the equivalent crack length and the visual crack length. The values obtained with the equivalent crack length agree with numerical values determined by the <em>J</em>-integral. A sensitivity analysis is carried out by the Monte Carlo method to elucidate the effect of uncertainties on the values of the equivalent crack length and on the energy release rates. Finally, the sensitivity analysis is done assuming that the value of the crack length is fixed, varying only the parameters concerning the energy release rates. The coefficients of variation obtained in this second case are slightly greater than in the first case.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111224"},"PeriodicalIF":4.7,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143937597","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xuanyu Ge, Linglong Zhou, Sara Bagherifard, Mario Guagliano
{"title":"A simple and efficient implementation of explicit phase field method in ABAQUS to address complex three-dimensional fracture problems","authors":"Xuanyu Ge, Linglong Zhou, Sara Bagherifard, Mario Guagliano","doi":"10.1016/j.engfracmech.2025.111222","DOIUrl":"10.1016/j.engfracmech.2025.111222","url":null,"abstract":"<div><div>The phase field method can analyze intricate crack growth behavior based on a regularized variational framework. However, its adaptation for solving complex three-dimensional fracture problems is still a challenge mainly due to the difficulties in numerical implementation and the corresponding high computational cost. In this study, a new ABAQUS implementation of explicit phase field method is proposed. The phase field was analogized to the temperature field, and the transient thermal variables in the heat transfer equation were rederived based on the explicit phase field governing equation. The dissipated inelastic energy was leveraged to characterize the volumetric heat flux, allowing the temperature field to be updated. To ensure the numerical stability, a series of novel formulations were proposed to restrict the temperature field rate. Finally, the field variables update approaches were described and the determination of critical time increment was discussed. Using this approach, the mechanical and thermal behaviors can be defined within a single user-defined material (VUMAT) subroutine. This implementation allows for convenient utilization of most built-in functions. Two classic brittle fracture examples were simulated firstly to test the efficiency of the proposed implementation through comparison with the UMAT implicit implementation, including comparison of computational time under different processors. Then the proposed implementation was exploited to solve several complex three-dimensional quasi-brittle fracture problems, and the corresponding predicted crack pattern and load–displacement data were compared with the experimental ones to validate the accuracy of the suggested model. The results revealed that the proposed explicit implementation is significantly more efficient than the implicit implementation of phase field, while maintaining the same accuracy. This approach can be leveraged to solve complex three-dimensional fracture problems considering mode I, mixed mode I + II and mixed mode I + III failure within an acceptable computational time. The proposed framework shows promise for efficient simulation of complex brittle/quasi-brittle fracture behavior in structural components. The source code is provided (<span><span>https://github.com/xuanyge/Explicit-PFM-VUMAT.git</span><svg><path></path></svg></span>) to enable interested researchers to utilize and implement it.</div></div>","PeriodicalId":11576,"journal":{"name":"Engineering Fracture Mechanics","volume":"323 ","pages":"Article 111222"},"PeriodicalIF":4.7,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143942436","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}