Pre-existing Cracks Research Articles

Monte Carlo simulations of the fracture resistance of an asphalt pavement layer

The purpose of the proposed numerical model is to analyze the cracking of the wearing course in a pavement overlay, assuming a pre-existing crack that passes through the binding layer and base. The computations employed the author's simulation-based Monte Carlo material model, which describes the failure process of a Semi-Circular Bend (SCB) specimen during standard laboratory testing of asphalt concrete. A key feature of this model is the incorporation of the random nature of material parameters, allowing for the simulation of result dispersion when analyzing a sufficiently large population of samples. The proposed FEM model and obtained material data were directly applied to the numerical analysis of the pavement structure. The comprehensive computational algorithm allows for a random description of the load that induces crack propagation in the pavement wearing course, leading to the creation of a histogram that defines the range of failure load dispersion. Such supporting calculations can assist in optimizing asphalt mix design and, in the future, may allow for the estimation of pavement structure reliability.

A novel hybrid PD-FEM-FVM approach for simulating hydraulic fracture propagation in saturated porous media

To enhance the computational efficiency of fluid–solid coupling in peridynamics (PD), a hybrid modeling approach based on the classical Biot theory is proposed for simulating hydraulic crack propagation in saturated porous media. The deformation and damage of solids are described by the coupling of the finite element method (FEM) and PD. Based on Darcy’s law, the finite volume method (FVM) is used to describe fluid seepage and calculate pore water pressure. The mutual transfer of fluid pressure and solid deformation is realized through the transition layer between the solid layer and the fluid layer. Firstly, the effectiveness of the proposed method is verified by a porous media seepage simulation example. Secondly, the ability and efficiency of this method to simulate crack propagation in saturated porous media are verified by several examples of hydraulic fracturing of rock with a single pre-existing crack. Finally, the synchronous hydraulic fracturing process of rock with double cracks is simulated. The ability of this method to simulate the simultaneous propagation of multiple fractures in the rock under fluid–solid coupling is further illustrated. The aforementioned studies demonstrate that the novel hybrid PD-FEM-FVM approach not only ensures computational accuracy and effectiveness but also significantly enhances computational efficiency.

Mechanical properties and durability of lightweight cenosphere cementitious composites under coupling effect of dry-wet cycles and salt erosion

Cenospheres are suitable as supplementary cementitious materials to partially replace cement in the production of environmentally friendly lightweight cenosphere cementitious composites (LCCC). In this study, the correlation between the content and particle size of the cenosphere, the early curing system and the performance of LCCC was analyzed. The macroscopic properties and microstructural changes of LCCC under the conditions of dry-wet cycles and salt erosion were studied with different early curing conditions. The durability changes of LCCC mortar under the coupling effect of salt erosion and dry-wet cycles were studied by examining the mass loss rate, relative dynamic modulus of elasticity and salt corrosion resistance coefficient of LCCC mortar. Finally, a comparative analysis of the impact of existing cracks on the performance of LCCC under the coupling effect of dry-wet and salt erosion was conducted. The main innovation of the study includes the mechanical properties and durability of LCCC under dry-wet-salt erosion conditions, as well as the clarification of the impact of prefabricated cracks on the mechanical properties and durability of LCCC. It is found that the flexural strength and the compressive strength of specimens cured at high temperature is generally lower than that of the specimens cured under normal temperature conditions. The wet-dry salt corrosion can exacerbate the erosion the outer shell and inner wall of cenospheres and disrupt the stability of the internal interface transition zone (ITZ) in LCCC. LCCC with CsS exhibits better macroscopic performance and microstructural properties. Pre-existing cracks accelerates the erosion in the cenosphere shell and inner wall under salt solution.

Effect of long-term ageing by temperatures on mechanical behavior of Fibre Reinforced Concretes

This study examines the effect of ageing by operating temperatures on the mechanical properties of Fibre Reinforced Concretes (FRCs) when exposed during long-term periods (180 days of conditioning). There is little information in the literature about long-term mechanical performance of steel (SFRC), and especially, macro synthetic fibre-reinforced concretes (MSFRCs), which may be affected due to this exposure. One type of steel fibre and three types of macro synthetic fibres were selected to reinforce an ordinary concrete in this work. An extensive experimental campaign (with a total of 260 specimens) which includes compression and flexural tests performed at temperatures, was carried out to investigate the influence of a range of temperature from -15 to 60 °C. Additionally, the interaction between cracking conditions (up to 0.5 mm) and temperature was analysed for both type of FRC beams. The compression results obtained show higher strengths at low temperatures (25-30 %) and a decrease on the strengths at higher temperatures (15-20 %), but these variations are not crucial on the ageing of FRCs. Similar effects were detected for the strength at the peak load, and to a minor extent, for the residual flexural strengths, for all the FRCs. Specimens with a pre-existing crack (cracked condition) exposed to the studied temperature had the same results as uncracked specimens. Thus, both FRCs have the potential to obtain appropriate results for structural applications when exposed to the studied temperatures.

Investigation of the dynamic mechanical response of corroded ultra-high performance fiber reinforced concrete (UHPFRC) with initial defects

This study addresses the dynamic mechanical response of the corroded ultra-high-performance fiber-reinforced concrete (UHPFRC) with initial defects, considering the possibility of corrosion deterioration induced by various pre-existing cracks during the long-term service life. For this purpose, an integrated accelerated corrosion method and Split Hopkinson Pressure Bar (SHPB)/high-speed camera etc. techniques are employed. Results show that increasing pre-impacting damage promotes the crack density and maximum width by 32.4%–62.3 % and 1.11–1.8 times, respectively. In terms of mechanical properties, coupling damages of initial defects and corrosion have adverse effects on the dynamic mechanical response. Typical fib Model Code 2010 applies to predict the DIF evolution of the corroded UHPFRC with initial defects. Numerous shear cracks are created at an angle along the weak interface as the corroded specimens with initial defects are again subjected to axial loading, revealing the associated failure mechanism. These results shed light on the dynamic response of corroded UHPFRC containing various initial defects and the failure mechanism gives some reference to service status evaluation.

New insights into fracture and cracking simulation of quasi-brittle materials based on the NMMD model

The simulation of crack-induced failure is of vital importance for the safety and integrity assessment of solids and structures but still challenging. Being the foremost approach, fracture mechanics is plagued by the logical inconsistency in dealing with crack and non-cracked region, and to make matters worse, it has been contradicted by the ongoing experiments that the fracture energy is not a material constant irrelevant to geometry even for the unsophisticated situations. In this paper, based on the recently proposed nonlocal macro-meso-scale consistent damage (NMMD) model, it is demonstrated that the logical consistency of governing equations either on crack-tip or off crack-tip plays a crucial role in the complete solution of cracking problems. To this end, the central themes of fracture mechanics, including the definition of an idealized crack, the separated equilibrium laws and the cracking criteria attached to crack-tip, are firstly revisited. Then, after a concise description of NMMD model with emphasis on the geometrical transmission form micro- to macro-scale and the physical–mechanical translation from geometry to energy, it is theoretically proved that the geometry-based damage zone reproduces idealized crack as a limit when the nonlocal characteristic length, i.e., the influence radius in NMMD model, tends to zero. As a result, the separated equilibrium laws of fracture mechanics can be covered by NMMD model with a unified set of governing equations that does not distinguish between points at crack surface and at non-cracked region. This advantage permits NMMD model to be well-suited to analyze both crack-induced failure problems with and without pre-existing cracks in the absence of additional cracking criteria. In contrast, owing to its macro-meso-scale consistent feature, a new criterion for kinking angle of crack, the maximum stress intensity factor (SIF) of mode-I, can be derived from NMMD model. Some representative numerical examples are also presented to further confirm the above intrinsic link to fracture mechanics. In particular, the superiority of the NMMD model is demonstrated by capturing not only the non-constant fracture parameters (especially the fracture energy) but also the cracking initiation not from the pre-existing crack-tips, which are revealed by experimental results but contradict fracture mechanics. These investigations imply that the NMMD model could provide a potential avenue to understand when and how cracks nucleate and propagate based on a two-scale perspective.

Local heating at the running crack tip in bcc iron according to molecular dynamics

Abstract This study presents estimates of a possible temperature rise at the crack tip from three dimensional (3D) atomistic simulations of fracture via molecular dynamics (MD) technique. Simulations start from an initial temperature of 0 K. The pre-existing edge crack was loaded in tension mode I. Crack initiation in MD is accompanied by surface dislocation emissions and later by a cross slip of the emitted dislocations into other slip systems. This leads both to stress waves radiation and to local heating in the plastic zone created by these dislocations. The local heating at the crack tip in the surface layers reaches a level of 480–500 K at some atoms, but an average temperature in the plastic zone is lower and depends on a chosen crack tip zone size. In the bulk crystal, where no dislocation emission (i.e. no plastic zone) has been realized, no significant heating is observed at the crack tip. MD results at the free sample surface comply with experimental data for ferritic steels with a pre-existing notch, loaded (quasi-statically) in mode I, as well as with some continuum predictions.

A comprehensive study on porosity modelling and its impact on fracture behavior of edge cracked FG structures using XIGA

In the present work, the fracture analysis of functionally graded (FG) porous structure containing an edge cracked is carried out in the presence of different types of porosity distributions using extended isogeometric analysis (XIGA). Firstly, the different types of porosity distribution functions are mathematically modeled in the porous FG structure across the length of the domain. The effective properties of the porous FG structure are computed using power law. Also, an additional term of porosity is incorporated in the power law to include the effect of porosity in the FG structure. The effective properties are computed across the length of the structure in the presence of different types of porosity distributions. Further, a pre-existing crack is modeled in the domain to study its influence on the fracture behaviour of porous FG structure using XIGA. To validate the accuracy, the results for the non-porous FG structure are compared with the available results in the literature (with the analytical and numerical solution), and they are found in good agreement (percentage error in the range of 0.04–––1.78%). Moreover, the comparative study is performed to investigate the influence of different types of porosity distributions on the fracture behaviour of FG structure.

Fatigue Characteristics Analysis of Carbon Fiber Laminates with Multiple Initial Cracks

In the entire wind turbine system, the blade acts as the central load-bearing element, with its stability and reliability being essential for the safe and effective operation of the wind power unit. Carbon fiber, known for its high strength-to-weight ratio, high modulus, and lightweight characteristics, is extensively utilized in blade manufacturing due to its superior attributes. Despite these advantages, carbon fiber composites are frequently subjected to cyclic loading, which often results in fatigue issues. The presence of internal manufacturing defects further intensifies these fatigue challenges. Considering this, the current study focuses on carbon fiber composites with multiple pre-existing cracks, conducting both static and fatigue experiments by varying the crack length, the angle between cracks, and the distance among them to understand their influence on the fatigue life under various conditions. Furthermore, this study leverages the advantages of Paris theory combined with the Extended Finite Element Method (XFEM) to simulate cracks of arbitrary shapes, introducing a fatigue simulation method for carbon fiber composite laminates with multiple cracks to analyze their fatigue characteristics. Concurrently, the Particle Swarm Optimization (PSO) algorithm is employed to determine the optimal weight configuration, and the Backpropagation neural network (BP) is used to train and adjust the weights and thresholds to minimize network errors. Building on this foundation, a surrogate model for predicting the fatigue life of carbon fiber composite laminates with multiple cracks under conditions of physical parameter uncertainty has been constructed, achieving modeling and assessment of fatigue reliability. This research offers theoretical insights and methodological guidance for the utilization of carbon fiber-reinforced composites in wind turbine blade applications.

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.

Flexural behavior and numerical simulation of reinforced concrete beams with a UHPC stay-in-place formwork

To overcome the shortcomings of traditional concrete structures such as long construction cycle, high construction cost and poor durability, seven composite beams with a UHPC stay-in-place formwork (UCB) and one reinforced concrete (RC) beam were designed. An experimental study was conducted to investigate the flexural behavior of composite beams under the static load. The research parameters include formwork thickness, reinforcement rate, and formwork surface treatment. The results revealed that roughening the surface of UHPC formwork has a positive impact on enhancing the integrity of composite beams. The utilization of UHPC stay-in-place formwork can effectively improve the stiffness, cracking load, and peak load of members. Compared with that of the RC beam, the cracking load and peak load of composite beams increased by 63% ~ 103% and 6% ~ 15%, respectively, and the yielding stiffness increased by 23% ~ 41%. Different from the RC beam, new cracks occurred on one side of the initial crack's end and continued to extend upwards rather than along the original crack. The existence of a significant quantity of tiny cracks in composite beams effectively slowed down the increase in the width of pre-existing cracks in the early stage of loading. Due to the multi-crack pattern exhibited by the UHPC formwork, the strain concentration of the steel bars was effectively mitigated. When subjected to the same load level, the strain of the longitudinal steel bars in composite beams was smaller than that of the RC beam. The formula for calculating the flexural bearing capacity of composite beams is established through theoretical analysis. In addition, the numerical model of the composite beam is established by the finite element method (FEM). The influence of the contact method of the UHPC-NC interface on the flexural performance of the composite beam is explored. The supplementary analysis of the parametric study of the reinforcement rate and the UHPC formwork thickness is carried out. When using ABAQUS for the numerical analysis of the composite beam, it is appropriate to adopt the cohesion model or the Tie model for the UHPC-NC interface, and it is not appropriate to use the Coulomb friction model.

Analysis of tensile failure mechanism in intact, foliated, and cracked rocks using distinct element method: Influence of anisotropy

A thorough numerical analysis was conducted to understand the failure mechanisms and behaviour in intact rocks, foliated rocks and rocks with pre-existing cracks using the Brazilian tensile strength tests. The dominant anisotropy in rocks, especially foliated rocks like phyllites, slates and schists, as well as sedimentary rocks like sandstone, limestones and shale, leads to complex failure mechanism. The presence of cracks leads to anisotropy causing different mechanisms of failure and variation in tensile strengths of rock mass. Numerical simulations were conducted based on experimental studies of foliated phyllite and dolomitic limestone from the Tejam Group, Lesser Himalaya. The variations in failure behaviour and strength of foliated rocks were studied for samples at different foliation angle using the particle distinct element model. The simulation results show that with increase in foliation angle to the loading direction (from 0° to 90°), the peak tensile strength of rocks remains almost same for specimens with θ = 0° to 30° and then increases linearly till 90°. The influence of cracks as anisotropy was also studied by changing crack length (from 5 mm to 30 mm) and angle (from 0° to 90°). The development of failure patterns in rocks with pre-existing cracks leads to a better understanding of the failure modes and makes it possible to interpret failure behavior, pattern, and strength parameters. With increase in crack length, the tensile strength of rocks decreases. Overall, this study depicts the influence of anisotropy and microparameters on failure modes and behavior in Brazilian tensile strength tests on foliated and pre-cracked rocks.

Experimental study on the failure characteristics and mechanism between spalling failure and rockburst

To reveal the differences in the failure characteristics and mechanisms between rockburst and spalling failure, the comparison experiment between spalling failure and rockburst occurred in a circular cavity were carried out. Experiment results reveal that rockburst is accompanied by the intense ejection of fragments, whereas spalling failure features with the slow detachment of fragments. Additionally, besides the plate and flaky fragment, the rockburst fragment exhibits more pronounced blocky features. In terms of strength and deformation characteristics, the peak stress for rockburst is 1.04 times that of spalling failure, whereas the plastic deformation of spalling failure is 1.75 times that of rockburst. Furthermore, considering the micro failure and energy evolution process, three interpretations were proposed to elucidate the distinct mechanisms between rockburst and spalling failure, including: (i) The swift invasion of shear cracks triggers the buckling failure of pre-existing plate-like fractures, causing the fragment ejection observed in rockburst. When shear cracks not substantially contribute to the evolutionary process, the continuous penetrating of the pre-existing tensile cracks will culminate in the gradual delamination and spalling failure of the sidewall. (ii) The energy dissipating behavior determines the distinct mechanisms of rockburst and spalling failure. Almost all of the energy input preceding rockburst is converted into the elastic strain energy of rock units, whereas the energy input preceding spalling failure results in greater energy dissipation (about 1.85 times that of rockburst). (iii) under the framework of Mohr-Coulomb strength criterion, it has been revealed that the rock units of spalling failure adhere to the minimum failure energy criterion. However, the energy accumulated by the rock units prior to rockburst surpasses the minimum energy necessary for failure. This residual energy provides the fractured rock mass with kinetic energy, ultimately leading to a violent ejection phenomenon of rockburst.

Mechanical and Lattice Thermal Properties of Si-Ge Lateral Heterostructures.

Two-dimensional (2D) materials have drawn extensive attention due to their exceptional characteristics and potential uses in electronics and energy storage. This investigation employs simulations using molecular dynamics to examine the mechanical and thermal transport attributes of the 2D silicene-germanene (Si-Ge) lateral heterostructure. The pre-existing cracks of the Si-Ge lateral heterostructure are addressed with external strain. Then, the effect of vacancy defects and temperature on the mechanical attributes is also investigated. By manipulating temperature and incorporating vacancy defects and pre-fabricated cracks, the mechanical behaviors of the Si-Ge heterostructure can be significantly modulated. In order to investigate the heat transport performance of the Si-Ge lateral heterostructure, a non-equilibrium molecular dynamics approach is employed. The efficient phonon average free path is obtained as 136.09 nm and 194.34 nm, respectively, in the Si-Ge heterostructure with a zigzag and armchair interface. Our results present the design and application of thermal management devices based on the Si-Ge lateral heterostructure.

Study on the mixed fracture characteristics of concrete-rock Brazilian disks with different fracture angles

The fracture characteristics of cracked straight-through Brazilian discs for concrete-rock composites (CSTBD-CRC) with different fracture inclination angles (Mode I-II) were investigated to evaluate the stability of concrete reinforcement engineering. Acoustic emission (AE) and digital image correlation (DIC) techniques were employed to study crack propagation and fracture modes in CSTBD-CRC Brazil splitting tests, as well as the entropy evolution of AE and DIC. Additionally, 3D scanning technology and fractal theory were utilized to investigate the effect of the fracture mode on the fracture surface morphology. The results indicate that both the fracture energy and peak load of the CSTBD-CRC increase with increasing fracture inclination angle. Moreover, the section fractal dimensions of 0°, 10°, 15° and 25° samples are 2.01879, 2.01528, 2.01186, and 2.01072. Shear fracture mode increases as the dip angle increases and the fracture surface morphology decreases in fractal dimension. Notably, sharp increases in entropy values for both AE and DIC occurred at peak load, suggesting their use as failure criteria for CSTBD-CRC specimens. Furthermore, cracks were found to initiate at the tips of preexisting cracks within the concrete. The displacement in the x direction undergoes the most significant changes on both sides of the crack during crack growth, while the strain in the x direction experiences its greatest variations at the crack location.

Peridynamics-based model of composite lamina with progressive variations in mechanical properties

An analytical model of the fiber-reinforced composite lamina is constructed based on bond-based peridynamics (PD), which utilizes variable matrix bonds and fiber bonds to evaluate the static properties of matrix and fiber materials, respectively. The bond stiffness of the matrix is determined by trigonometric functions, and fibers are represented by constants. The proposed PD model describes the continuous variations in the static properties of composites as a function of fiber angle. Both the conventional analytical method and the PD method are used to solve the failure-free elastic deformation issue. The two methods estimate maximum errors in the displacement of material points in the x and y directions of 0.01 % and 4.86 %, respectively, indicating good consistency. The failure propagation of the preexisting crack layer under the axial tensile load is simulated and the obtained failure modes are in good agreement with experimental results. The proposed composite lamina model simulates the influence of fiber direction on overall behavior through various matrix bond stiffness. Breaking through the limitation of constant matrix bond stiffness in conventional models, is conducive to promoting the further development of bond-based composite PD models.

Effects of alkali-silica reaction on the fracture energy of mass concrete: Experimental investigations using the wedge splitting test

This work aims to study the effects of alkali-silica reaction (ASR) on the fracture energy of reactive mass concrete mixtures (maximum aggregate sizes of 38 mm and 76 mm), characteristic of an existing hydro-electric facility. More than 40 large concrete blocks were casted and stored in a controlled environmental chamber for more than 1100 days. The blocks as well as cylindrical reference cores were monitored in expansion. Wedge splitting (WS) tests with different notch to depth ratios, were performed at different ASR advancements to assess the evolution of mechanical and fracture parameters of concrete. Digital Image Correlation (DIC) technique was used to monitor the tests. Damage Rate Index (DRI) technique was also performed on cores extracted from the blocks, following the completion of the splitting tests. Important size and heterogeneity effects were observed on the assessed ASR expansion and fracture energy. By filtering these effects, it was possible to assess for the first time in literature the evolution of size-independent fracture energy with respect to intrinsic volumetric expansion of concrete. The findings suggest that fracture energy increases for moderate levels of expansion, due to crack branching within the splitting crack, interacting with pre-existing ASR-induced cracks. This phenomenon was not observed for the case of a few horizontally cast blocks where the splitting crack ran parallel to pre-existing ASR cracks, indicating an anisotropy in ASR damage effects on the fracture energy of concrete.

Integrating machine learning and the finite element method for assessing stiffness degradation in photovoltaic modules

This study introduces a novel machine learning (ML) method utilizing a stacked auto-encoder network to predict stiffness degradation in photovoltaic (PV) modules with pre-existing cracks. The input data for the training process was derived from numerical simulations, ensuring a comprehensive representation of module behavior under various conditions. The findings highlight the robust predictive capability of the model, as evidenced by its impressive R2 value of 0.961 and notably low root mean square error (RMSE) of 4.02%. These metrics significantly outperform those of other conventional methods, including the artificial neural network with R2 of 0.905 and RMSE of 9.43%, the space vector machine with R2 of 0.827 and RMSE of 17.93%, and the random forest (RF) with R2 of 0.899 and RMSE of 11.02%. Moreover, the findings suggest that the predictive dynamics of degradation are affected by the varying weight functions of different input parameters, such as climate temperature (CT), grain size (GS), material effort, and pre-crack size, as the degradation level changes. Furthermore, a geometric analysis reveals model deficiencies where significant overestimations correlate with thicker glass components, while pronounced underestimations are predominantly associated with thinner layers of polycrystalline silicon wafer and Ethylene Vinyl Acetate in the module. As a case study, it demonstrated that to maintain a constant degradation level between 1.30 and 1.32 in a PV module with components featuring consistent geometric attributes, the input parameters must be kept within specific ranges: CT ranging from 33 °C to 57 °C, GS ranging from 36 to 81 μm, material effort ranging from 0.74 to 0.81, and pre-crack size ranging from 24 to 32 μm. Therefore, this underscores that the ML model not only predicts degradation but also delineates the parameter space required to achieve a consistent output value.

Fracture Mechanism and Damage Constitutive Model of Freeze–Thaw Fissured Granite Subjected to Fatigue Loading

The failure of rock in cold regions due to repeated freeze–thaw (F-T) cycles and periodic load-induced fatigue damage presents a significant challenge. This study investigates the evolution of the multi-scale structure of fractured granite under combined freeze–thaw (F-T) cycles and periodic loading and develops a constitutive damage model. The results indicate that after F-T cycles, network cracks develop around pre-existing cracks, accompanied by block-like spalling. After applying the fatigue load, the nuclear magnetic resonance (NMR) T2 spectrum shifts to the right, significantly increasing the amplitude of the third peak. The freeze–thaw process induces a “liquid–solid” phase transition, weakening the original pore structure of the rocks and leading to meso-damage accumulation. The pores in fractured granite progressively enlarge and interconnect, reducing the rock’s load-bearing capacity and fatigue resistance. The combined effects of F-T cycles and periodic loading induce particle movement and alter fracture modes within the rock, subsequently affecting its macro-damage characteristics. The theoretical curves of the constitutive model align with the experimental data. The findings can serve as a theoretical reference for preventing and controlling engineering disasters in fractured rock masses in cold regions.

Experimental and numerical investigation on crack propagation in biomimetic nacreous composites with gradient structures

Combining interlocking structures and gradient design, this paper proposes a 2D biomimetic nacreous composite material with a gradient interlocked structure. Experimental and numerical simulation methods are employed to investigate the initiation and propagation behavior of cracks in the regular interlocking and gradient interlocked biomimetic nacreous composite materials. Samples of biomimetic nacreous composite materials with pre-existing cracks are manufactured using 3D printing technology, and uniaxial tensile tests are conducted to explore their crack propagation behavior. For biomimetic nacreous composite materials with periodic interlocking structures, the interlocking angle is found to play an important role in the toughening mechanism of composites. In biomimetic composite materials with large interlocking angles, crack inhibition effects are achieved through structural gradient design, effectively preventing catastrophic failure. Within the numerical analysis framework, a cohesive zone modeling approach is employed to represent the softer phase of the material, and the numerical simulation are validated by direct comparison with experimental results. In addition, a comprehensive numerical analysis is carried out to comprehend how the structural gradient affects the spread of cracks in these nacreous composites. This research not only illuminate the underlying deformation and toughening mechanisms intrinsic to interlocking nacreous composites but also pave the way for innovative strategies in the design of biomimetic materials through the incorporation of structural gradients.