Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CRACKS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electrochemical performance of ionic polymer metal composite under tensile loading

Abstract Cracks in polymer composites can lead to premature failure, which can be disastrous for polymer-based energy storage devices. Detecting these cracks is essential to guarantee the reliability and safety of such devices. However, detecting cracks in composite polymers such as ionic polymer metal composites (IPMCs) is a challenging task, which makes it difficult to ensure their performance and safety. The overall goal of this study is to investigate the effect of cracks or damage caused by tensile loading on the mechanical properties and electrochemical characteristics of IPMC based capacitors. During tensile testing, the deformation of the IPMC strips causes changes in the ion distribution and concentration in the polymer matrix, influencing the performance of the material. The measurements were conducted utilizing electrochemical impedance spectroscopy at a room temperature ( 21 ∘ C ) and frequency range of 10 KHz to 1 Hz. The method utilized in this study proved to be easy and quick with consistent results. The IPMC capacitor was found to increase its capacitance after major cracking in the Pt electrodes from high tensile mechanical loads. Furthermore, at lower frequency range (<100 Hz), the real ( ε ′ ) and imaginary ( ε ′ ′ ) part of permittivity increase with the addition of loads. This displays that the dielectric property of the material is affected due to the increasing of the loads. It is concluded that, at frequencies above 100 Hz, the permittivity is weakly load dependent.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evolution of texture and internal stresses within polycrystalline rock salt using in situ 3D synchrotron computed tomography and 3D X-ray diffraction

Over time rock salt caverns have been extensively used as reliable repositories for hazardous waste such as nuclear waste, oil or compressed gases. Undisturbed rock salt deposits in nature are usually impermeable and have very low porosity. However, rock salt formations under excavation stresses can develop crack networks, which increase their porosities; and in the case of a connected crack network within the media, rock salt may become permeable. Although the relationship between the permeability of rock salt and the applied stresses has been reported in the literature, a microscopic study that investigates the properties influencing this relationship, such as the evolution of texture and internal stresses, has yet to be conducted. This study employs in situ 3D synchrotron micro-computed tomography and 3D X-ray diffraction (3DXRD) on two small-scale polycrystalline rock salt specimens to investigate the evolution of the texture and internal stresses within the specimens. The 3DXRD technique measures the 3D crystal structure and lattice strains within rock salt grains. The specimens were prepared under 1D compression conditions and have shown an initial {111} preferred texture, a dominant {110}$\langle1\bar{1}0\rangle$ slip system and no fully connected crack network. The {111} preferred texture under the unconfined compression experiment became stronger, while the {111}$\langle1\bar{1}0\rangle$ slip system became more prominent. The specimens did not have a fully connected crack network until applied axial stresses reached about 30 MPa, at a point where the impermeability of the material becomes compromised due to the development of multiple major cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inelastic relaxation processes in amorphous sodium silicates

During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na 2 O-SiO 2 (NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na 2 O-SiO 2 (NS15) composition, the 25% Na 2 O-SiO 2 (NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na 2 O-SiO 2 compositions. Coordination analysis of the Na + ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na + ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stress Intensity Thresholds for Development of Reliable Brittle Materials

Brittle material failure in high consequence systems can appear random and unpredictable at subcritical stresses. Gaps in our understanding of how structural flaws and environmental factors (humidity, temperature) impact fracture propagation need to be addressed to circumvent this issue. A combined experimental and computational approach composed of molecular dynamics (MD) simulations, numerical modeling, and atomic force microscopy (AFM) has been undertaken to identify mechanisms of slow crack growth in silicate glasses. AFM characterization of crack growth as slow as 10 -13 m/s was observed, with some stepwise crack growth. MD simulations have identified the critical role of inelastic relaxation in crack propagation, including evolution of the structure during relaxation. A numerical model for the existence of a stress intensity threshold, a stress intensity below which a fracture will not propagate, was developed. This transferrable model for predicting slow crack growth is being incorporated into mission-based programs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multiphysics Meshfree Degradation Modeling of Energy Storage Materials with Kernel Enrichment

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. Reference: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

25 ENERGY STORAGE↗
DOE OSTI2026

Kernel Enriched Meshfree Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Laboratory of the Rockies (NLR), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. References: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Application of Convolutional Neural Networks (CNNs) to Recognize Defects in 3D-Printed Parts

Cracks and pores are two common defects in metallic additive manufacturing (AM) parts. In this paper, deep learning-based image analysis is performed for defect (cracks and pores) classification/detection based on SEM images of metallic AM parts. Three different levels of complexities, namely, defect classification, defect detection and defect image segmentation, are successfully achieved using a simple CNN model, the YOLOv4 model and the Detectron2 object detection library, respectively. The tuned CNN model can classify any single defect as either a crack or pore at almost 100% accuracy. The other two models can identify more than 90% of the cracks and pores in the testing images. In addition to the application of static image analysis, defect detection is also successfully applied on a video which mimics the AM process control images. The trained Detectron2 model can identify almost all the pores and cracks that exist in the original video. This study lays a foundation for future in situ process monitoring of the 3D printing process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Life after 30 years - a PV system in Colorado

The PV system consisted of 155 monocrystalline Si Arco M53 panels, installed in 1987. 80 were installed on a fixed tilt rack, and 75 on 5 separate dual-axis passive trackers. The system was decommissioned in 2017 due a campus re-design and building expansion at the Platte River Power Authority (PPRA) site. Prior to decommissioning a site visit established the following failure type observations: EVA browning (100%), cracked cells (60%), and cracked wires or grommets (9%, or more if replacements were made selectively on cracked parts). EVA browning is typical in modules of this vintage, whereas the occurrence of so many cracked cells is surprising and is likely weather-related. The fairly frequent occurrence of cracked wires and grommets serves as a reminder that accelerated testing of these components is important for goals of low-cost 50-year energy generation. In smaller numbers, thermal effects (7%), bent frames (1%), and cell corrosion (<1%) were observed. A sample of 20 modules was selected and brought to NREL for in-depth characterization. Because baseline measurements were not available using the nameplate rating, which can be inaccurate, revealed a mean module sample degradation of 0.7 %/year. Sporadic production data during the 30 years determined a system degradation rate closer to 1.2 %/year. Subsequently, four of these modules were re-mounted outdoors at NREL and are now approaching a total lifetime of 34 years.

durability↗
DOE OSTI2022

Examination Of IASCC In 304 Stainless Steel Core Shroud From A Commercial Boiling Water Reactor (BWR)

The 304 stainless steel core shroud and the accompanying weld material from a commercial boiling-water reactor (BWR) experienced an estimated dose of ~3.5 dpa neutron irradiation leading to irradiation-assisted stress corrosion cracking (IASCC). The resulting IASCC crack tips are examined with high-resolution analysis techniques, including scanning electron microscopy (SEM), scanning transmission electron microscopy (STEM), and atom probe tomography (APT), with select examples of directly correlated STEM/APT observations. Analytical observations reveal the crack morphology, composition, microstructure, and any associated corrosion-induced or radiation-induced segregation in the local metal grain boundaries ahead of the propagating SCC cracks. The observed oxide composition and structure were in-line with expectations, consisting of a Cr-rich inner oxide and a Fe-rich outer oxide, with the leading crack tip oxide being Cr-rich. In addition, the leading grain boundary exhibited typical radiation-induced segregation signatures that were largely unperturbed by the accompanying corrosion. Finally, the quantitative capabilities and limitations of STEM and APT in the context of corrosion fronts and grain boundary chemistry from the directly correlated measurements are discussed.

Olszta, Matthew J.↗
DOE OSTI2022

Damage Tolerance Assessment of Laser Clad Repairs of Coarse Grain Ti-6Al-4V

Background: Laser cladding offers a promising approach to repairing aviation components, with the ability to restore the geometric shape and the static and fatigue strength. Objective: The damage tolerance of a coarse grain Ti-6Al-4 V material is assessed via three sets of samples: (i) baseline, (ii) laser clad repaired, and (iii) laser cladding followed by heat treatment. Method: The specimens were subjected to fatigue loading using a marker band technique, in which cracks initiated from pre-defined laser slots and the small crack growth rates were measured by quantitative fractography. Additional characterization was performed on these samples, including hardness measurements, microstructure analysis, and multimodal residual stress analyses. Results: The region of the specimen undergoing laser cladding exhibited higher hardness values and a compressive residual stress, based on the slitting method. The samples subjected to heat treatment after cladding displayed negligible differences in grain size and minimal residual stress relaxation. Conclusions: The results demonstrated significant variability, with data demonstrating that cracks in the laser clad region (exhibiting the largest compressive residual stress) grew faster than cracks in the substrate region of the sample or in the baseline cases. Finally, these counterintuitive results are rationalized with additional analyses from synchrotron-based energy dispersive X-ray diffraction and high energy X-ray diffraction, demonstrating significant spatial variability in the residual elastic strains in the millimeter-sized prior β grains within the material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Toughening of interface networks through the introduction of weak links

The resistance of polycrystalline materials to intergranular fracture depends on the toughness of the interfaces between neighboring grains. Using finite element simulations on model, two-dimensional microstructures, we demonstrate that the steady-state crack growth resistance of a polycrystal that contains exclusively high toughness interfaces is suboptimal: its crack growth resistance may be increased by introducing a minority fraction of weak interfaces at random locations within the interface network. Furthermore, the weaker are these interfaces, the greater is their toughening effect. These weak interfaces toughen the material through crack deflection and branching, nucleation of secondary cracks, and inelastic dissipation within intact interfaces. Furthermore, our findings suggest that fracture resistance-oriented material development should not aim to eliminate all weak interfaces, but rather judiciously incorporate a low fraction of them.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

IASCC of 304 SS in BWR environments: Effects of post-irradiation annealing and surface condition

To investigate the impact of low-temperature post-irradiation annealing (PIA) on the stress corrosion cracking of neutron-irradiated 304 stainless steel, constant extension rate tests were conducted in simulated boiling water reactor normal water chemistry (BWR-NWC) and hydrogenated water chemistry (HWC) environments. Ten tensile samples, comprising five as-irradiated and five PIA specimens, were prepared by electropolishing the gauge section of electric discharge machined (EDM) samples. Here, the annealing treatment reduced the yield strength from approximately 550 MPa to around 425 MPa and significantly restored the ductility and the strain hardening capability of the alloy. Consequently, the susceptibility of this material to irradiation-assisted stress corrosion cracking (IASCC) was effectively mitigated, which is more prominent in HWC, as evident from fractography, which indicated a decreased propensity for intergranular (IG) fracture. Furthermore, it was observed that the polished surface facilitated crack initiation more readily than the EDM surface, suggesting that the EDM process suppressed crack initiation to some extent.

Crack initiation↗
DOE OSTI2024

Assessing the release, transport, and retention of radioactive aerosols from hypothetical breaches in spent fuel storage canisters

Interim dry storage of spent nuclear fuel involves storing the fuel in welded stainless-steel canisters. Under certain conditions, the canisters could be subjected to environments that may promote stress corrosion cracking leading to a risk of breach and release of aerosol-sized particulate from the interior of the canister to the external environment through the crack. Research is currently under way by several laboratories to better understand the formation and propagation of stress corrosion cracks, however little work has been done to quantitatively assess the potential aerosol release. The purpose of the present work is to introduce a reliable generic numerical model for prediction of aerosol transport, deposition, and plugging in leak paths similar to stress corrosion cracks, while accounting for potential plugging from particle deposition. The model is dynamic (changing leak path geometry due to plugging) and it relies on the numerical solution of the aerosol transport equation in one dimension using finite differences. The model’s capabilities were also incorporated into a Graphical User Interface (GUI) that was developed to enhance user accessibility. Model validation efforts presented in this paper compare the model’s predictions with recent experimental data from Sandia National Laboratories (SNL) and results available in literature. We expect this model to improve the accuracy of consequence assessments and reduce the uncertainty of radiological consequence estimations in the remote event of a through-wall breach in dry cask storage systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Study of the Printability, Microstructures, and Mechanical Performances of Laser Powder Bed Fusion Built Haynes 230

The nickel-based superalloy, Haynes 230 (H230), is widely used in high-temperature applications, e.g., heat exchangers, because of its excellent high-temperature mechanical properties and corrosion resistance. As of today, H230 is not yet in common use for 3D printing, i.e., metal additive manufacturing (AM), primarily because of its hot cracking tendency under fast solidification. The ability to additively fabricate components in H230 attracts many applications that require the additional advantages leveraged by adopting AM, e.g., higher design complexity and faster prototyping. In this study, we fabricated nearly fully dense H230 in a laser powder bed fusion (L-PBF) process through parameter optimization. The efforts revealed the optimal process space which could guide future fabrication of H230 in various metal powder bed fusion processes. The metallurgical analysis identified the cracking problem, which was resolved by increasing the pre-heat temperature from 80 °C to 200 °C. A finite element simulation suggested that the pre-heat temperature has limited impacts on the maximum stress experienced by each location during solidification. Additionally, the crack morphology and the microstructural features imply that solidification and liquation cracking are the more probable mechanisms. Both the room temperature tensile test and the creep tests under two conditions, (a) 760 °C and 100 MPa and (b) 816 °C and 121 MPa, confirmed that the AM H230 has properties comparable to its wrought counterpart. The fractography showed that the heat treatment (anneal at 1200 °C for 2 h, followed by water quench) balances the strength and the ductility, while the printing defects did not appreciably accelerate part failure.

14 SOLAR ENERGY↗
DOE OSTI2024

An experimental process parameter study on the identification of defects in additively fabricated Al6061 with laser powder bed fusion

Additively fabricated metal parts using laser powder bed fusion (L-PBF) possess sophisticated morphology due to the recurrent use of laser-induced metal powder melting and solidification. The surface and 3D morphology of these parts often include defects in the form of protrusions, depressions, pores, voids, keyholes, or cracks that are known to be influenced by laser scanning paths and layer-to-layer processing. Such inconsistent part quality hampers the extensive adoption of L-PBF. Pores and cracks are detrimental to the fatigue life of the parts and components. Quantifying and controlling part defects and optimizing processing and scanning strategy parameters adaptively in real-time through in situ monitoring systems are highly desired. This study investigates the optimization of experimental process parameters (power, scan velocity, and hatch spacing) and their effects on the cracking and porosity of Al6061 alloy using machine learning techniques. Multi-objective optimization is formulated and conducted to determine the L-PBF parameters that minimize both porosity and crack densities.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗