Crack Propagation Process Research Articles

Biologically inspired adaptive crack network reconstruction based on slime mould algorithm.

The dynamic crack propagation trajectories play a crucial role in enhancing our understanding of spatial mechanisms involved in crack expansion. However, visualization of internal cracks under complex crack conditions has always been a challenge. Biological networks have been honed by many cycles of evolutionary selection pressure and are likely to yield reasonable solutions to such combinatorial optimization problems. This study applied the slime mould algorithm to improve the accuracy of internal crack localization in rocks and employed Minimum spanning tree and Gaussian mixture model to construct the crack propagation trajectories. By introducing the concept of bond length, the evolution characteristics of crack levels were effectively characterized. Research results showed that this approach effectively preserves essential crack localization information while mitigating the influence of interfering parameters, providing crack characterization results that exhibit high consistency with actual fracture patterns. The curves of cumulative bond length and relative bond length over time conform to the trend of a Growth/Sigmoidal curve. The strength of the bond was correlated with the temporal process of crack propagation. This result could be helpful for analyzing crack trajectories and predicting rock stability.

Promoted cutting and energy transition efficiency by an asymmetric cutter

Laboratory and numerical tests were performed to understand the rock breakage mechanism of a CCS (constant cross section) cutter and an asymmetric cutter. The laboratory and numerical tests consistently show that the asymmetric cutter frequently generates larger breakage areas and consumes less indentation energy than the CCS cutter, but there is a defect spacing of 70 mm. Thus, the asymmetric cutter frequently has a higher cutting efficiency than the CCS cutter. In addition, the numerical tests reveal that the CCS cutter first generates a plastic zone and subsequently forms cracks to connect with the defect tips. In this process, a large proportion of the energy is wasted in the particle friction in the plastic zone. However, the asymmetric cutter tends to form much smaller plastic zones, which makes it difficult to generate particle friction. Only a small amount of friction energy is consumed in the later crack propagation process. Thus, the asymmetric cutter has a higher energy transition efficiency than the CCS cutter.

Experimental study on the Mode l fracture toughness of frozen silty clay incorporating Digital image correlation

To investigate the effects of initial moisture content and temperature on the Mode I fracture toughness (KIc) of frozen silty clay, a series of three-point bending tests were conducted. Rectangular specimens with prefabricated cracks were tested, and digital image correlation (DIC) technology was employed to analyze the microscopic characteristics at the crack tip. The results indicate that the Mode I fracture toughness of frozen silty clay increases with decreasing temperature and increasing initial moisture content. Temperature governs the failure mode, with plastic failure predominantly occurring at high temperatures and brittle failure at low temperatures. Based on the load–displacement curves and DIC recordings, the macroscopic fracture process of the specimens can be categorized into three stages: elastic deformation, formation of the failure surface, and specimen failure. Additionally, the crack propagation process can be further divided into three stages: initiation and development of microcracks, transition from microcracks to macroscopic cracks, and rapid development of macroscopic cracks. These findings provide critical insights for slope stability research in cold regions and offer new perspectives for engineering design and construction in similar environments.

Dynamic tensile behaviors and crack propagation in fiber-reinforced cementitious composites: Experimental investigation and Peridynamic simulation

Fiber-Reinforced Cementitious Composites (FRCC) have gained significant attention in engineering applications due to their superior mechanical properties and toughness, particularly under high strain rate conditions. This study performed dynamic tensile tests on FRCC at high strain rates (35-110 s⁻¹) using a Split Hopkinson Tensile Bar (SHTB) apparatus. Additionally, a novel Peridynamic (PD) model was developed for the SHTB and FRCC system, leveraging the advanced capabilities of the emerging PD theory. The study compared and analyzed dynamic tensile strength, ultimate tensile strain, strain rate effects, failure modes, and crack development in FRCC with different fiber ratios at various high strain rates, using both experimental data and PD simulations. The results show that the PD-SHTB-FRCC dynamic model developed in this study exhibits high consistency between numerical simulations and experimental findings, effectively capturing the processes of crack initiation, propagation, and complete failure in FRCC specimens. The dynamic tensile properties of FRCC improve significantly with increased strain rates, with polyethylene (PE) fibers providing superior reinforcement compared to steel fibers. Notably, the dynamic tensile strength, peak tensile stress, and ultimate tensile strain of FRCC increase significantly with rising strain rates, with specimens containing higher PE fiber content showing a more pronounced enhancement effect. For strain rates between 42.6 s⁻¹ and 76.1 s⁻¹, considering dynamic tensile strength, ultimate tensile strain, and peak tensile stress, the optimal combination for resisting dynamic tensile loads was 1.5% PE fibers and 0.5% steel fibers. At a strain rate of 99.8 s⁻¹, a 2% PE fiber ratio alone provided the best performance.

Acoustic emission and 4D X-ray micro-tomography for monitoring crack propagation in rocks

Acoustic emission (AE) and 4D X-ray computed tomography (4D XCT) were used simultaneously to study crack initiation and propagation in two different types of quartz-rich sandstones during the four-point bending experiments. Statistical analysis of the AE response indicated the failure mechanisms and their dynamics. The characteristic changes observed in the AE response defined the timing of the bending interruptions for XCT scanning to reveal the development of the crack. It was possible to quantitatively describe the developing cracks in their dimensions and volume and relate this information to the rate of decrease in the post-peak region of the material response. It could be concluded that the combination and concurrent use of AE and XCT techniques represents a highly effective and reliable instrument for observation, description, analysis of the crack propagation process, and rock disintegration in detail at a microscale level. With regard to the specific sandstones studied, Mšené sandstone is softer, respectively, less brittle, while Kocbeře sandstone is characterised by a more brittle behaviour accompanied by an AE signal with higher amplitudes compared to those of Mšené.

Study on the impact of grouting reinforcement on the mechanical behavior of non-penetrating fracture sandstone

The grouting reinforcement technology is an essential method to enhance the mechanical performance of fractured rock masses and the effectiveness of reinforcement varies with different grouting materials. To further understand the mechanical improvement capabilities of each grout and the reinforcement mechanisms at the grout-rock interface, this study prepared samples with different grouting materials (sulphoaluminate cement (SAC), ultra-fine cement (UFC), and epoxy resin (EPR)) and the uniaxial compression tests were conducted. Based on these tests, the macro and micro mechanical characteristics of different grouting samples were revealed using particle image velocimetry (PIV), acoustic emission (AE), scanning electron microscopy (SEM), and nuclear magnetic resonance (NMR). The results indicate that grouting helps improve the mechanical performance and deformation resistance of fractured rock masses. It effectively limited lateral displacement of the samples, reduced stress concentration at fracture tips, enhanced shear effects during sample fracture, and altered the crack propagation process and failure modes. Compared to the fractured samples, the peak strength of SAC, UFC, and EPR samples increased by 17.8 %, 23.4 %, and 28.3 %, and the elastic modulus increased by 14.3 %, 7.9 %, and 24.8 %, respectively. Among these, the EPR samples exhibited a similarity in parameter indicators to intact samples of over 85 %, making EPR the optimal grouting material. The degree of grout-rock fusion is the primary factor influencing grouting reinforcement effectiveness. SAC is covering-type cement, UFC is embedded cement, EPR is a fusion material, and the fusion-type materials are more beneficial for improving the mechanical performance of fractured rocks.

Structural characteristics and evolution in the crack propagation process of NiCrAlY/YSZ gradient coatings through thermal shock at 900 °C

In this study, four different thicknesses of NiCrAlY/8YSZ gradient coatings were fabricated using atmospheric plasma spraying (APS) technology. The failure behavior and lifetime of coatings were subsequently compared and analyzed through thermal shock experiments. The results revealed a pronounced correlation between the initiation and propagation of internal cracks within the gradient coatings and the thickness of the gradient region at a thermal shock temperature of 900°C. Both horizontal and vertical cracks were observed to propagate and interconnect, with horizontal crack propagation being the primary cause of coating failure. Notably, samples featuring a gradient region thickness of 150μm exhibited superior thermal shock resistance, suggesting promising prospects for enhanced coating durability.

A coupled scaled boundary finite element and phase-field algorithm for seismic loading

Seismic-induced damage, integral to the safety evaluations of major engineering projects, persists as a key focus of research worldwide. Based on the Scaled Boundary Finite Element Method (SBFEM) and the Phase-Field Method (PFM), a coupled algorithm tailored for reciprocal loading was introduced in this paper, integrating strategies including "closure constraints," "numerical threshold strategy," and "subdomain block integration." Adopting object-oriented principles, a universal coupling solution framework has been built and seamlessly embedded within GEODYNA, a self-developed finite element software system. The accuracy was validated through rigorous benchmark tests. The entire process of crack initiation, propagation, and dynamic opening-closing cycles in the Koyna concrete gravity dam was reproduced. Furthermore, the effect of mesh size and computational timestep on the structural seismic response, crack localization, and the extent of damage in the dam were explored. The outcomes demonstrated that the SBFEM-PFM coupling algorithm performs effectively and meets the engineering precision criteria for seismic evaluations and reinforcement analyses of crucial structures.

Study on Multi-Crack Damage Evolution and Fatigue Life of Corroded Steel Wires Inside In-Service Bridge Suspenders

The parallel steel wires used in arch bridge suspenders experience random corrosion damage on their surfaces during service. Corrosion damage, including micro-cracks, pitting, and a combination of both, leads to significant stress concentration under axial loading, which affects the performance of the steel wires. The change in the stress field caused by surface damage alters the stress intensity factor at the crack tip, and the presence of adjacent crack tips significantly amplifies the stress intensity factor, thereby accelerating crack propagation. The development of small surface damages in the steel wires is difficult to control and observe through experiments. By utilizing finite element methods for simulation, it is possible to intuitively analyze the crack propagation process, the trend of stress changes at the crack tip, and the interaction between damages. Numerical simulation results based on Paris’ law indicate that corrosion pits have a certain impact on the stress intensity factor at the crack tip. The propagation process of coplanar double cracks is highly sensitive to the initial crack size and the distance between adjacent crack tips. When the crack spacing is less than the crack depth, the stress intensity factor at the adjacent crack tips exhibits significant amplification. Based on this phenomenon, the coplanar double-crack system can be simplified to a complete single crack for analysis. By comparing the fatigue life of the double-crack system with that of the equivalent single crack, the effectiveness of the simplification rule has been validated.

Experimental and Numerical Research on Fracture Properties of Mass Concrete Under Quasi-Static and Dynamic Loading

The dynamic fracture behavior of mass concrete is crucial to the dynamic analysis and safety evaluation of concrete dams subjected to strong earthquake shocks in the framework of fracture mechanics. In the presented research, cylindrical specimens with a ring of preset cracks were cast by three-graded mass concrete, and direct tension tests were performed with two loading rates considered, i.e., 10−6/s for quasi-static loading and 10−3/s for dynamic loading. The load–crack mouth opening displacement (P-CMOD) curves were obtained, from which the fracture toughness, fracture energy, and characteristic length of the mass concrete were obtained. In this process, the influence of the eccentricity in the tests was compensated by the numerical modeling of the tests. Next, the crack propagation process of the mass concrete was modeled using the extended finite element method. From the test results, it is found that, under quasi-static loading, the crack generally propagates along the interface between the aggregates and the matrix, while, under dynamic loading, more aggregates are fractured. As compared to the case of quasi-static loading, the energy absorption capacity, fracture energy, and fracture toughness increase for dynamic loading, while the characteristic length decreases. Moreover, the numerically predicted P-CMOD curves agree reasonably well with the experimental measurements.

An efficient moving-mesh strategy for predicting crack propagation in unidirectional composites: Application to materials reinforced with aligned CNTs

This paper presents an efficient numerical approach for reproducing the process of crack propagation inside unidirectional composites subjected to general loading conditions, with special reference to epoxy materials enhanced with embedded aligned CNTs. This approach involves a traditional FE framework improved by the Moving Mesh (MM) technique based on the Arbitrary Lagrangian-Eulerian (ALE) formulation and the Interaction Integral Method (M−integral). The MM serves as a powerful numerical tool to simulate the discrete crack advance with minimal remeshing, thus reducing the computational complexities. Instead, the M−Integral method, formulated for generally anisotropic materials, has been employed to extract the mixed-mode Stress Intensity Factors (SIFs), which are necessary to define the crack onset condition and propagation direction on the basis of the modified Maximum Hoop Stress Criterion. The proposed strategy includes the extended rule of mixtures to evaluate the homogenized elastic properties of nano-reinforced composites. The validity of the proposed methodology has been assessed through comparisons with experimental data and numerical results available in the literature.

A data-driven prediction for concrete crack propagation path based on deep learning method

Accurate and efficient prediction of concrete crack propagation path facilitates timely maintenance. To avoid errors arising from assumptions in physical models and parameter measurements, a data-driven prediction method for concrete crack propagation path is proposed. The spatio-temporal prediction approach is presented by unifying the conditional Generative Adversarial Network and Long Short-Term Memory Network (UcGLS-Net). The UcGLS-Net has the advantage of enabling temporal (by Long Short-Term Memory Network, LSTM) and spatial information (by the conditional Generative Adversarial Network, cGAN) in the concrete crack propagation process to be treated simultaneously and distinctly thus improving prediction accuracy. Notably, results demonstrate that UcGLS-Net still works in crack propagation path prediction even with unknown aggregate distribution. Moreover, the effects of historical input quantity and prediction intervals on prediction accuracy of the UcGLS-Net are studied and the network framework is optimized. Due to high efficiency, this study provides a new adaptable and real-time prediction method for concrete crack propagation path, offering support for the design, damage prediction, assessment, and maintenance of concrete structures.

Experimental study of type-I crack propagation in rock monitored by fiber Bragg grating

The occurrence of rock mass engineering disasters can be reduced significantly by obtaining the key information in the process of rock internal fracture accumulation and evolution through advanced monitoring means and taking certain preventive measures. Through fiber Bragg grating monitoring tests of the type-I crack propagation in rock, the spectral characteristics of the fiber Bragg grating in the process of rock failure were analyzed, the bandwidth expansion was extracted to characterize the local tensile strain during the crack propagation process, and the damage variable based on the bandwidth expansion was calculated to describe the type-I crack propagation process. The research results showed the following: (1) The fiber Bragg grating generated a non-uniform local tensile strain during type-I crack propagation. The bandwidth expansion of the wave spectrum was positively correlated with the local tensile strain. The change in the characteristics of the wave spectrum could provide the local real strain in the process of crack propagation. (2) During the steady propagation of the type-I crack, the bandwidth increased obviously, while the number of acoustic emission events was relatively small, and most of them were in the process zone. The gentle increase in the bandwidth corresponded to the unstable growth stage of the type-I crack, which could be used as a precursor of rock sample instability. During loading, the maximum tensile strain of the type-I crack was 6 %, while the maximum shear strain reached 1.8 %, primarily attributed to tensile failure. The average value of fracture toughness KIC, based on the equivalent linear elastic fracture mechanics model, was determined to be 1.21 MPa·m1/2, and the average value of fracture energy Gf was calculated as 57.7 N/m. The variation of peak reflectivity in fiber Bragg grating was associated with the closure of pre-existing defects and the propagation of microcracks in rocks, resulting in unstable shear strain at crack extension sites. (3) The damage variable D defined by the bandwidth broadening began to become concentrated and accumulated in the crack steady growth stage, and it increased slowly in the crack intensification and evolution stage. The grating was sensitive to changes when it was pasted to the surface at 90°, but it could easily fracture or become unable to capture spectral signals in time. Based on a comprehensive comparison, the monitoring effect was the best when the grating was pasted at 60° relative to the crack.

Fatigue life prediction for Q420qFNH weathering steel welded joints considering the effect of multiple cracks

Multiple semi-elliptical cracks frequently develop at the weld toes of welded joints. These interfering multiple cracks significantly accelerate crack growth rates and reduce the fatigue life of a welded joint compared to a single crack. This study proposes a method for predicting the fatigue life of a welded joint with multiple cracks, considering weld toe amplification effect, multiple crack interference amplification effect, and crack closure effect. The method is validated through fatigue tests on Q420qFNH weathering steel welded joints. The main factors influencing the interference amplification effect are investigated, and theoretical formulas for the interference amplification factor of the stress intensity factor are derived using the finite element method and the superposition principle. The initiation, aggregation, and morphological evolution during the propagation process of multiple cracks is accurately simulated, and the fatigue life of both cruciform-welded and butt-welded joints is predicted. Additionally, a quantitative analysis is conducted to assess the impact of the number of fatigue cracks on fatigue strength. The results indicate that the interference amplification effect of cracks depends on their distance and mainly affects the surface point where the crack front is closest. The proposed fatigue life prediction method can accurately estimate the fatigue life of Q420qFNH steel welded joints considering the effect of multi-cracks. The interference and coalescence of multiple cracks significantly reduce the fatigue life of welded joints. An increase in the number of fatigue cracks notably decreases the fatigue strength.

Hydro-mechanical modeling of cohesive crack propagation of concrete lining in high internal pressure tunnels

High pressure tunnels with concrete lining have been extensively utilized in project practice. However, due to the characteristic of concrete being susceptible to cracking under tension, the lining inevitably develops cracks under high internal water pressure, posing a serious threat to the operation of tunnels. This study aims at developing a hydro-mechanical numerical model of cohesive crack propagation of concrete lining in high internal pressure tunnels. In this regard, the determination of cohesive element parameters is elucidated, the contact simulation within the software ABAQUS is improved to accurately characterize the interface between lining and surrounding rock, and the numerical calculation process in ABAQUS is realized using indirect coupled method. The simulation results align well with the physical model test and engineering monitoring data, demonstrating that the proposed method can accurately simulate the hydraulic interactions of high pressure tunnel. Additionally, a comparison with calculation models employing tie constraints to simulate the lining-surrounding rock interface is conducted. Finally, comparison with traditional continuum method reveals that while both methods exhibit consistent overall trends. It is recommended to choose the proposed method when describing the discontinuous propagation process of cracks, which cannot be simulated by the continuum analysis method.

Study on the influence of slit pipe wall thickness on the rock breaking effect of slit charge explosion

AbstractIn this study, the influence of slit pipe wall thickness on shaped charge blasting efficiency was investigated. The influence of different slit pipe wall thicknesses on the dynamic behavior of slotted charge blasting was analyzed by a combination of experimental observations and numerical simulations. Polymethyl methacrylate (PMMA) was used as a rock simulation material to prepare slit pipes with different slit pipe wall thicknesses (0.5, 1.0, and 1.5 mm). A dynamic caustics system is used to capture the crack propagation process. ANSYS/LS‐DYNA was used for the numerical simulation analysis. The experimental results show that there is a positive correlation between the slit pipe wall thickness and the crack propagation length, velocity, and dynamic stress intensity factor. With increasing slit pipe wall thickness, the directional propagation and energy concentration ability of cracks are significantly improved, but an increase in the slit pipe wall thickness will also promote the propagation of secondary cracks. In addition, the numerical simulation results support the experimental observations and reveal the existence of an optimal slit pipe wall thickness that can maximize the blasting efficiency while maintaining structural integrity. This study not only deepens the understanding of the physical laws of the blasting process but also provides a practical reference for the optimal design of slotted charges in engineering applications.

Effect of grain boundary ferrite morphologies on impact toughness of X80 girth weld

The X80 girth weld microstructure formed by high-strength acicular ferrite and low-strength grain boundary ferrite was studied to investigate the effect of grain boundary ferrite morphologies on impact toughness. The results show that due to the difference in direction of heat loss caused by the narrow groove, the grain boundary ferrite morphology of the weld metal center and the weld metal edge are continuous network structure and parallel layered structure, respectively. The continuous network structure grain boundary ferrite causes strain localization, resulting in the formation of a large number of secondary cracks in the crack propagation process, which seriously deteriorates the toughness.

Toughening Mechanism in Nanotwinned Boron Carbide: A Molecular Dynamics Study.

Boron carbide ceramics are potentially ideal candidates for lightweight bulletproof armor, but their use is currently limited by their low fracture toughness. Recent experimental results have shown that sintered samples with high twin densities exhibit high fracture toughness, but the toughening mechanism and associated crack propagation process of nanotwinned boron carbide at the atomic scale remain a mystery. Reported here are molecular dynamics simulations with a reactive force field potential to investigate how nanoscale twins affect the fracture toughness of boron carbide ceramics. The results show that the strength disparity on either side of a twin boundary is the fundamental reason for the toughening effect; the twin boundary impedes crack propagation only when the crack moves to a region of higher fracture strength. The fracture toughness of nanotwinned boron carbide is greatly affected by the angle between the twin boundary and the prefabricated crack. At an angle of 120°, the twin boundary provides the maximum toughening effect, enhancing the toughness by 32.72%. Moreover, phase boundaries-another common structure in boron carbide ceramics-have no toughening effect. This study provides new insights into the design of boron carbide ceramics with high fracture toughness.

Crack propagation process in double-flawed granite under compression using digital image correlation method and numerical simulation

The existence of flaws seriously weakens the rock strength, directly affects the crack expansion morphology, and indirectly affects the slope stability. In this study, uniaxial compression tests were carried out on double-flawed granite to investigate the effects of flaw angle on its compressive strength and crack coalescence based on DIC technology and PFC2D numerical simulation. The result indicates that the UCS values of flawed specimens at angles of 15°, 30°, 45°, and 60° are approximately 25.5–51.7% lower than that of the intact specimen. The failure strength of the double-flawed sample increases with the increase of flaw inclination angle, and the fracture morphology shifts from no expansion to split expansion. When the flaw tip strain of each flaw sample exceeded 0.6%, the stress concentration was generated at the flaw tip, and the sample began to appear macroscopic large cracks. DIC technique can well observe the crack initiation and propagation process, recording inclined shear strain bands at flaw tips, tensile strain bands in the middle of flaws, and tensile shear O-shaped strain ring in rock bridge areas. Numerical simulations using PFC2D software were carried out and showed a good agreement with the physical results. In addition, the effect of structural inclination on specimen failure strength is clearly explained by the theory of compressive damage based on the projected size of flaws. The deformation of rock mass is not a simple material deformation but is composed of material deformation and structural deformation. These experimental and numerical results enhance our understanding of crack initiation and coalescence characteristics, aiding in the analysis of rock structure stability in scenarios such as excavated underground openings, slopes, and tunneling construction, where pre-existing cracks or step-path fractures are pivotal to structural integrity.

Investigation on mixed mode I/II crack propagation in nitrate ester plasticized polyether propellant: Experimental and numerical study

The fracture of solid propellant is predominantly attributed to the existence of mixed mode cracks, so it is essential to investigate the the fracture behavior of solid propellant with mixed mode I/II crack. This paper presents fracture characteristics of nitrate ester plasticized polyether (NEPE) propellant under different crack inclination angles (β = 30°–90°). Based on the combination of a drawing machine and a high-speed camera, the mechanical response, crack propagation velocity and crack-path morphology were investigated. The critical equivalent stress intensity factor Keqc was calculated to assess the fracture toughness of the NEPE propellant, and a potential simplified criterion related to the stress intensity factor was proposed. The experimental results demonstrated that the NEPE propellant with mixed mode I/II crack exhibited blunting fracture phenomena during crack propagation, resulting in fluctuating crack propagation velocity. As the crack inclination angle decreases, the fracture toughness of the NEPE propellant increases and then decreases, and the value of Keqc reaches its maximum at β = 45°. Furthermore, numerical studies based on bond-based peridynamic (BBPD) were performed by modeling the crack propagation process of the NEPE propellant, including the crack phase field diagram and the load–displacement curve of the NEPE propellant. The simulation results were then compared with the experiments.