Pub Date : 2025-06-10DOI: 10.1007/s10704-025-00862-2
T. Corre, F. Hild, V. Lazarus
Sharp kinks may be observed under shear loading or in materials containing weak directions, such as those produced by additive manufacturing. A better understanding of the fracture of these materials, both theoretically and experimentally, is required to deploy them in structural applications. This study focuses on the measurement of stress intensity factors (SIFs) around a sharp kink using digital image correlation (DIC). The performances of two DIC-based techniques, namely, integrated-DIC and post-processing of DIC-measured displacement fields, are assessed on a benchmark test using fused deposit modeling capabilities, and are compared to a reference finite element solution. It is shown that Williams’ expansion remains valid on a large enough region around the crack to extract reliable SIFs even very close to the crack kink. Both techniques are very trustworthy, provided the SIF identification zone is carefully defined to exclude the kink zone of influence.
{"title":"Stress intensity factor determination along a kinked crack path by DIC analyses","authors":"T. Corre, F. Hild, V. Lazarus","doi":"10.1007/s10704-025-00862-2","DOIUrl":"10.1007/s10704-025-00862-2","url":null,"abstract":"<div><p>Sharp kinks may be observed under shear loading or in materials containing weak directions, such as those produced by additive manufacturing. A better understanding of the fracture of these materials, both theoretically and experimentally, is required to deploy them in structural applications. This study focuses on the measurement of stress intensity factors (SIFs) around a sharp kink using digital image correlation (DIC). The performances of two DIC-based techniques, namely, integrated-DIC and post-processing of DIC-measured displacement fields, are assessed on a benchmark test using fused deposit modeling capabilities, and are compared to a reference finite element solution. It is shown that Williams’ expansion remains valid on a large enough region around the crack to extract reliable SIFs even very close to the crack kink. Both techniques are very trustworthy, provided the SIF identification zone is carefully defined to exclude the kink zone of influence.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145143108","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-06-05DOI: 10.1007/s10704-025-00858-y
GaoFei Ji, LingHui Hu
This study proposes a fatigue life prediction method combining small-sample data expansion with the Weibull distribution function, incorporating the first order reliability factor (FOSM) to improve accuracy. Using Generalized Polynomial Chaos Expansion (GPC) and Latin Hypercube Sampling (LHS), small-sample fatigue data is expanded, followed by enhancing the two-parameter Weibull model with FOSM. Results show the generalized polynomial chaotic expansion method and Latin hypercube sampling are used to obtain the probability density curve when the stress level is 350 MPa, and the original data are all on this probability density curve, indicating that the expansion method is more credible. High prediction precision within a 1.5 × error range, with logarithmic safety life linearly related to stress level and decreasing with higher failure probability.
{"title":"Fatigue life prediction method based on polynomial chaotic expansion and Weibull distribution","authors":"GaoFei Ji, LingHui Hu","doi":"10.1007/s10704-025-00858-y","DOIUrl":"10.1007/s10704-025-00858-y","url":null,"abstract":"<div><p>This study proposes a fatigue life prediction method combining small-sample data expansion with the Weibull distribution function, incorporating the first order reliability factor (FOSM) to improve accuracy. Using Generalized Polynomial Chaos Expansion (GPC) and Latin Hypercube Sampling (LHS), small-sample fatigue data is expanded, followed by enhancing the two-parameter Weibull model with FOSM. Results show the generalized polynomial chaotic expansion method and Latin hypercube sampling are used to obtain the probability density curve when the stress level is 350 MPa, and the original data are all on this probability density curve, indicating that the expansion method is more credible. High prediction precision within a 1.5 × error range, with logarithmic safety life linearly related to stress level and decreasing with higher failure probability.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145142685","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-28DOI: 10.1007/s10704-025-00851-5
Jianing Xie, Krishnaswamy Ravi-Chandar
The rapid development of 3D printing of 316L stainless steel thin-walled structures obtained by direct energy deposition has generated an increased interest in the mechanical properties of such materials for use in applications; in particular, failure models are needed to ensure structural reliability. We consider the response of uniaxial tension specimens, with and without notches, to characterize the constitutive and failure behavior of the material. Specifically, we use numerical simulations of the notched tension experiment, achieved with a simple power-law strain hardening model and a failure criterion based on attaining a triaxiality-dependent critical strain-to-failure, to demonstrate that this model is capable of reproducing the material behavior accurately.
{"title":"Simulation of ductile failure of a notched-tension specimen of 3D printed 316L stainless steel","authors":"Jianing Xie, Krishnaswamy Ravi-Chandar","doi":"10.1007/s10704-025-00851-5","DOIUrl":"10.1007/s10704-025-00851-5","url":null,"abstract":"<div><p>The rapid development of 3D printing of 316L stainless steel thin-walled structures obtained by direct energy deposition has generated an increased interest in the mechanical properties of such materials for use in applications; in particular, failure models are needed to ensure structural reliability. We consider the response of uniaxial tension specimens, with and without notches, to characterize the constitutive and failure behavior of the material. Specifically, we use numerical simulations of the notched tension experiment, achieved with a simple power-law strain hardening model and a failure criterion based on attaining a triaxiality-dependent critical strain-to-failure, to demonstrate that this model is capable of reproducing the material behavior accurately.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145145525","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-28DOI: 10.1007/s10704-025-00860-4
Shucheta Shegufta, Michael Zaiser
A particular failure mode of highly porous brittle materials consists in the propagation of cracks under uniaxial compressive loads. Such ’anticracks’ have been observed in a range of materials, from snow and porous sandstone to brittle foams. Here we present a computational model for the formation and propagation of anticrack-type failure in porous materials within the general computational framework of bond-based peridynamics. Random porosity is represented, on a scale well above the characteristic pore size, by random bond deletion (dilution disorder). We apply our framework to experimental data on anticrack propagation in silicate foams.
{"title":"Peridynamics based model of anticrack-type fracture in brittle foams","authors":"Shucheta Shegufta, Michael Zaiser","doi":"10.1007/s10704-025-00860-4","DOIUrl":"10.1007/s10704-025-00860-4","url":null,"abstract":"<div><p>A particular failure mode of highly porous brittle materials consists in the propagation of cracks under uniaxial compressive loads. Such ’anticracks’ have been observed in a range of materials, from snow and porous sandstone to brittle foams. Here we present a computational model for the formation and propagation of anticrack-type failure in porous materials within the general computational framework of bond-based peridynamics. Random porosity is represented, on a scale well above the characteristic pore size, by random bond deletion (dilution disorder). We apply our framework to experimental data on anticrack propagation in silicate foams.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10704-025-00860-4.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145145524","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We introduce the fast convolution-based method (FCBM) for a peridynamic correspondence (cPD) model to simulate finite plastic deformations and ductile fracture in 3D. The cPD model allows the direct use of classical finite plasticity constitutive ductile failure models, like the Johnson–Cook (J-C) model used here. We validate the FCBM for the cPD model against experimental results from the literature on ductile failure in Al2021-351 alloy samples of various geometries. Notably, calibration of elastic and hardening material parameters is made only using the experimental data from the simplest geometry, a smooth round bar, and only up to the necking point. We then use that calibrated model beyond necking, through full failure, and for all the different sample geometries. The performance (speedup and memory allocation) of the new method is compared versus the meshfree method normally used to discretize PD models for fracture and damage. The proposed method leads to efficient large-scale peridynamic simulations of finite plastic deformations and ductile failure that are closer to experimental measurements in terms of displacement and plastic strain at failure than previous FEM-based solutions from the literature.
{"title":"A fast convolution-based method for peridynamic models in plasticity and ductile fracture","authors":"Farzaneh Mousavi, Siavash Jafarzadeh, Florin Bobaru","doi":"10.1007/s10704-025-00849-z","DOIUrl":"10.1007/s10704-025-00849-z","url":null,"abstract":"<div><p>We introduce the fast convolution-based method (FCBM) for a peridynamic correspondence (cPD) model to simulate finite plastic deformations and ductile fracture in 3D. The cPD model allows the direct use of classical finite plasticity constitutive ductile failure models, like the Johnson–Cook (J-C) model used here. We validate the FCBM for the cPD model against experimental results from the literature on ductile failure in Al2021-351 alloy samples of various geometries. Notably, calibration of elastic and hardening material parameters is made only using the experimental data from the simplest geometry, a smooth round bar, and only up to the necking point. We then use that calibrated model beyond necking, through full failure, and for all the different sample geometries. The performance (speedup and memory allocation) of the new method is compared versus the meshfree method normally used to discretize PD models for fracture and damage. The proposed method leads to efficient large-scale peridynamic simulations of finite plastic deformations and ductile failure that are closer to experimental measurements in terms of displacement and plastic strain at failure than previous FEM-based solutions from the literature.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 3","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-05-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144100376","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-10DOI: 10.1007/s10704-025-00852-4
Daniela V. Klein, Magnus Boåsen, Pål Efsing, Jonas Faleskog
The fracture surfaces of 49 SE(B) toughness tests performed on five different geometries, were carefully investigated by SEM imaging and cross-section analysis. The specimens were extracted from a large multi-pass weld in T-S orientation. The failure characteristics were associated with three distinctly different zones of the weld. Transgranular fracture occurred primarily in the reheated zone and in the as-welded zone with a dendritic microstructure inclined relative to the crack plane. With a dendritic microstructure aligned with the crack plane intergranular fracture occurred. The toughness of the as-welded zone with a low inclination angle, was significantly lower than that obtained in the other two weld zones. Due to the relatively large size of the zones compared to the fracture process zones of the tests, it is appropriate to characterize the failure behavior as large-scale heterogeneity. Weakest-link modeling may be applied locally in each weld zone, giving rise to three different sets of model parameters. A new calibration technique is introduced and used to fit a local weakest-link model to the toughness distribution curves of the individual zones.
{"title":"Assessment of large-scale heterogeneity due to toughness variations in a multipass weld: brittle failure mechanisms and modeling","authors":"Daniela V. Klein, Magnus Boåsen, Pål Efsing, Jonas Faleskog","doi":"10.1007/s10704-025-00852-4","DOIUrl":"10.1007/s10704-025-00852-4","url":null,"abstract":"<div><p>The fracture surfaces of 49 SE(B) toughness tests performed on five different geometries, were carefully investigated by SEM imaging and cross-section analysis. The specimens were extracted from a large multi-pass weld in T-S orientation. The failure characteristics were associated with three distinctly different zones of the weld. Transgranular fracture occurred primarily in the reheated zone and in the as-welded zone with a dendritic microstructure inclined relative to the crack plane. With a dendritic microstructure aligned with the crack plane intergranular fracture occurred. The toughness of the as-welded zone with a low inclination angle, was significantly lower than that obtained in the other two weld zones. Due to the relatively large size of the zones compared to the fracture process zones of the tests, it is appropriate to characterize the failure behavior as large-scale heterogeneity. Weakest-link modeling may be applied locally in each weld zone, giving rise to three different sets of model parameters. A new calibration technique is introduced and used to fit a local weakest-link model to the toughness distribution curves of the individual zones.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 2","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-05-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10704-025-00852-4.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143932325","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-30DOI: 10.1007/s10704-025-00850-6
Saiedeh Marashi, Hamidreza Abdolvand
This paper aims to numerically investigate the nucleation and propagation of microcracks in dual phase Zirconium (Zr) containing both Hexagonal Close-Packed (HCP) α-Zr and Body Centered Cubic (BCC) β-Zr crystals. For this purpose, a subroutine that incorporates different damage criteria is coupled with a crystal plasticity finite element model to investigate the effects of crystals elastic and plastic anisotropy. Attention is given to the role of the BCC β-phase in the crack nucleation of notched zirconium polycrystals. First, the maximum shear strain accumulated on the predominant slip system is used as the crack initiation criterion. The modeling results reveal that for single phase HCP α-grains cracks lie on the prismatic planes, but for dual phase α/β cases, cracks may lie on either basal or prismatic planes depending on the α/β crystal orientations, and the adjacent β-phase features such as its thickness or distance from the notch. Moreover, numerical results indicate that the presence of thin layered β-phase hinders crack propagation, regardless of its geometrical or crystallographic features. The performance of other damage criteria is also discussed. Lastly, it is shown that in comparison to α-grains undergoing cyclic loads, the crack propagation rate is reduced in β-crystals.
{"title":"The role of β-phase on crack nucleation and propagation in dual phase zirconium polycrystals: a crystal plasticity finite element modeling","authors":"Saiedeh Marashi, Hamidreza Abdolvand","doi":"10.1007/s10704-025-00850-6","DOIUrl":"10.1007/s10704-025-00850-6","url":null,"abstract":"<div><p>This paper aims to numerically investigate the nucleation and propagation of microcracks in dual phase Zirconium (Zr) containing both Hexagonal Close-Packed (HCP) α-Zr and Body Centered Cubic (BCC) β-Zr crystals. For this purpose, a subroutine that incorporates different damage criteria is coupled with a crystal plasticity finite element model to investigate the effects of crystals elastic and plastic anisotropy. Attention is given to the role of the BCC β-phase in the crack nucleation of notched zirconium polycrystals. First, the maximum shear strain accumulated on the predominant slip system is used as the crack initiation criterion. The modeling results reveal that for single phase HCP α-grains cracks lie on the prismatic planes, but for dual phase α/β cases, cracks may lie on either basal or prismatic planes depending on the α/β crystal orientations, and the adjacent β-phase features such as its thickness or distance from the notch. Moreover, numerical results indicate that the presence of thin layered β-phase hinders crack propagation, regardless of its geometrical or crystallographic features. The performance of other damage criteria is also discussed. Lastly, it is shown that in comparison to α-grains undergoing cyclic loads, the crack propagation rate is reduced in β-crystals.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 2","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143892676","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-04-28DOI: 10.1007/s10704-025-00846-2
Abdallah Salama, Ahmed Elsayed, Atef Eraky, Rania Samir
This paper investigates the application of the Virtual Element Method (VEM) for simulating crack propagation in 2D marble rock under linear elastic fracture mechanics (LEFM) conditions. The inherent mesh flexibility of VEM is leveraged by employing an adaptive mesh refinement (AMR) strategy based on recovery by compatibility in patches (RCP) for triangular, quadrilateral, and even polygonal meshes. The accuracy and efficiency of crack path prediction are enhanced by calculating stress intensity factors (SIFs) and T-stress through the interaction domain integral method coupled with the Generalized Maximum Tangential Stress (GMTS) criterion. The effectiveness of this approach is validated using three distinct marble rock specimens with varying material properties and initial crack configurations: semi-circular bend (SCB) Harsian Marble, center-cracked circular disk (CCCD) limestone, and edge-cracked triangular (ECT) Neyriz Marble. The GMTS criterion, incorporating three parameters (KI, KII, and T), precisely predicts crack initiation and propagation directions, demonstrating its superiority for mixed-mode fractures.
{"title":"Adaptive virtual element method with RCP for mixed-mode fracture analysis of marble rocks using GMTS criterion","authors":"Abdallah Salama, Ahmed Elsayed, Atef Eraky, Rania Samir","doi":"10.1007/s10704-025-00846-2","DOIUrl":"10.1007/s10704-025-00846-2","url":null,"abstract":"<div><p>This paper investigates the application of the Virtual Element Method (VEM) for simulating crack propagation in 2D marble rock under linear elastic fracture mechanics (LEFM) conditions. The inherent mesh flexibility of VEM is leveraged by employing an adaptive mesh refinement (AMR) strategy based on recovery by compatibility in patches (RCP) for triangular, quadrilateral, and even polygonal meshes. The accuracy and efficiency of crack path prediction are enhanced by calculating stress intensity factors (SIFs) and T-stress through the interaction domain integral method coupled with the Generalized Maximum Tangential Stress (GMTS) criterion. The effectiveness of this approach is validated using three distinct marble rock specimens with varying material properties and initial crack configurations: semi-circular bend (SCB) Harsian Marble, center-cracked circular disk (CCCD) limestone, and edge-cracked triangular (ECT) Neyriz Marble. The GMTS criterion, incorporating three parameters (KI, KII, and T), precisely predicts crack initiation and propagation directions, demonstrating its superiority for mixed-mode fractures.</p></div>","PeriodicalId":590,"journal":{"name":"International Journal of Fracture","volume":"249 2","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10704-025-00846-2.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143883666","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}