Cracks In Materials Research Articles

T-stress extraction in arbitrarily cracked orthotropic composites with the numerical manifold method and Stroh formalism

The escalating utilization of composites over the past decades has necessitated the fracture investigation of orthotropic materials. The T-stress, namely, the first non-singular term of the normal stress parallel to the crack in William’s expansion, is of great significance for fracture analysis. In this work, the numerical manifold method (NMM) is advanced to assess the T-stress of arbitrary-shaped cracks (including non-intersecting cracks and multi-branched cracks) in two-dimensional orthotropic composites. Attributing to the bi-cover systems, the NMM can conveniently discretize the physical domain and naturally accommodate the discontinuity across crack surface. Meanwhile, the singularity at crack tip can be well captured by the wise choice of local approximation function. Through the application of interaction integral technology in the NMM postprocessing, the T-stress is extracted with the Stroh-form auxiliary fields. The accuracy of the proposed method is verified by comparing with reference solutions and then applied to arbitrarily branched and intersecting cracks. The results indicate that the present approach has convincing accuracy and also considerable convenience in T-stress evaluation. Additionally, the impacts of material orientations, crack geometries and loading conditions on the T-stress are also investigated.

Analysis, Assessment, and Mitigation of Stress Corrosion Cracking in Austenitic Stainless Steels in the Oil and Gas Sector: A Review

This comprehensive review examines the phenomena of stress corrosion cracking (SCC) and chloride-induced stress corrosion cracking (Cl-SCC) in materials commonly used in the oil and gas industry, with a focus on austenitic stainless steels. The study reveals that SCC initiation can occur at temperatures as low as 20 °C, while Cl-SCC propagation rates significantly increase above 60 °C, reaching up to 0.1 mm/day in environments with high chloride concentrations. Experimental methods such as Slow Strain Rate Tests (SSRTs), Small Punch Tests (SPTs), and Constant-Load Tests (CLTs) were employed to quantify the impacts of temperature, chloride concentration, and pH on SCC susceptibility. The results highlight the critical role of these factors in determining the susceptibility of materials to SCC. The review emphasizes the importance of implementing various mitigation strategies to prevent SCC, including the use of corrosion-resistant alloys, protective coatings, cathodic protection, and corrosion inhibitors. Additionally, regular monitoring using advanced sensor technologies capable of detecting early signs of SCC is crucial for preventing the onset of SCC. The study concludes with practical recommendations for enhancing infrastructure resilience through meticulous material selection, comprehensive environmental monitoring, and proactive maintenance strategies, aimed at safeguarding operational integrity and ensuring environmental compliance. The review underscores the significance of considering the interplay between mechanical stresses and corrosive environments in the selection and application of materials in the oil and gas industry. Low pH levels and high temperatures facilitate the rapid progression of SCC, with experimental results indicating that stainless steel forms passive films with more defects under these conditions, reducing corrosion resistance. This interplay highlights the need for a comprehensive understanding of the complex interactions between materials, environments, and mechanical stresses to ensure the long-term integrity of critical infrastructure.

Low‐velocity impact damage behavior of glass fiber reinforced composites with different crack propagation

AbstractGlass fiber‐reinforced resin matrix composites are widely used and studied for their good mechanical and impact resistance properties. In this paper, five glass fiber reinforced nylon 6 (PA6)‐based composites were prepared by mechanical blending combined with the molding process. The flexural strength, shear strength, pendulum impact strength, and low‐velocity impact properties of five composites were comparatively investigated. The crack extension modes of the composites were characterized and analyzed by light microscopy, scanning electron microscopy (SEM), and computerized X‐ray tomography (XCT), and the effects of different crack extension modes on the mechanical and impact resistance properties of the composites were researched. Experimental results show that: (1) Cracks in the composite material mainly exist in two expansion modes: transverse propagation parallel to the direction of fiber layup and longitudinal propagation perpendicular to the direction of fiber layup; (2) Different crack extension modes will cause different damage situations to the composites, however, the crack extension mode has minimal impact on the flexural, shear and pendulum impact strength of the composites; (3) The flexural and shear strengths of the composites co‐reinforced using fly ash and glass fibers were preferably 484.5 and 46.1 MPa; (4) The crack extension mode has a significant impact on low‐velocity impact performances of the composite material. The composite co‐reinforced with glass beads (GB) and glass fibers had the best impact resistance, with maximum impact force (Pmax), residual impact force (Pr), and absorbed energy (Ea) of 7015.2 N, 6273.6 N and 15.4 J, respectively.Highlights The propagation of cracks inside the composites was observed using SEM and XCT etc. A relationship between crack propagation mode and low‐velocity impact resistance of composites was established. The low‐velocity impact properties of composites are carefully analyzed.

Study on Electrochemical Characteristics of Crystal Structure Changes Effects of Silicon Anode Materials for Lithium Ion Batteries

Graphite is widely used as a representative anode material in commercialized lithium-ion secondary batteries. However, with a theoretical capacity of 370 mAh/g, it is difficult to manufacture high-capacity/high-density lithium-ion secondary batteries. Therefore, research is actively underway to find high-capacity substances to replace graphite. Silicon anode materials are currently attracting attention as next-generation high-capacity cathode materials. Silicon forms Li4.4Si when fully charged through an alloying/dealloying process with lithium during the reaction. Graphite is intercalation/deintercalation in reaction with lithium and LiC6 is formed when fully charged. That's why silicon has a theoretical capacity of 3800 mAh/g, about 10 times that of graphite.However, silicon is a high-capacity electrode material that overcomes the theoretical capacity limit of graphite, and has the disadvantage of experiencing volume expansion of about 400% during repeated charge/discharge cycles. Problems such as material cracks due to volume expansion can affect the safety of the battery. Various studies are being conducted to solve silicon-related challenges to effectively utilize silicon materials. Previous studies have explored the use of binders or additives to mitigate volume expansion. In addition, various types of silicon materials such as Si alloy and SiOx are being developed to suppress volume expansion through matrix formation. Nevertheless, despite the development of silicon materials, structural changes in the charging and discharging process of silicon remain unclear.This study aims to analyze the structural changes that occur during the charging and discharging process of silicon. The experiment evaluates the electrochemical characteristics of crystalline silicon and amorphous silicon and compares the behavior and structural changes during charging/discharging. The volumetric expansion measurements are performed with electrochemical characterization, and ex-situ X-ray diffraction (XRD) and transmission electron microscope (TEM) analyses are performed for comparison with crystal structural changes.

Operando chemo-mechanical evolution in LiNi0.8Co0.1Mn0.1O2 cathodes.

Ni-rich LiNi x Co y Mn z O2 (NCMxyz, x+y+z=1, x≥0.8) layered oxide materials are considered the main cathode materials for high-energy-density Li-ion batteries. However, the endless cracking of polycrystalline NCM materials caused by stress accelerates the loss of active materials and electrolyte decomposition, limiting the cycle life. Hence, understanding the chemo-mechanical evolution during (de)lithiation of NCM materials is crucial to performance improvement. In this work, an optical fiber with με resolution is designed to in operando detect the stress evolution of a polycrystalline LiNi0.8Co0.1Mn0.1O2 (P-NCM811) cathode during cycling. By integrating the sensor inside the cathode, the stress variation of P-NCM811 is completely transferred to the optical fiber. We find that the anisotropy of primary particles leads to the appearance of structural stress, inducing the formation of microcracks in polycrystalline particles, which is the main reason for capacity decay. The isotropy of primary particles reduces the structural stress of polycrystalline particles, eliminating the generation of microcracks. Accordingly, the P-NCM811 with an ordered arrangement structure delivered high electrochemical performance with capacity retention of 82% over 500 cycles. This work provides a brand-new perspective with regard to understanding the operando chemo-mechanical evolution of NCM materials during battery operation, and guides the design of electrode materials for rechargeable batteries.

Study on flexural and self-repairing properties of shape memory alloy concrete beams

The generation and expansion of cracks in concrete materials have always been a problem that needs to be solved, and an increasing number of scholars are beginning to hope that structures will have bionic functions of self-diagnosis and self-repairing of damage. In this study, shape memory alloy (SMA) concrete beams are prepared. The flexural and self-repairing properties of SMA concrete beams were investigated by a four-point bending test. The effects of various factors on the bearing capacity, ductility, and crack width before and after repair were analyzed. The experimental results indicate that the increase of SMA fiber content and the complexity of the SMA bar end shape can significantly enhance the flexural bearing capacity and repair performance of the beam. The increase in the bond length of SMA bars can enhance the bending resistance of the beam, but it will reduce its self-repairing performance. The repair rate of SMA on concrete cracks is 45.7 %-76.5 %, and the load recovery rate is between 49.4 % and 80.6 %. The flexural load capacity calculation formula is proposed, and the calculation results closely match the test values, providing a reference for the application of SMA.

X-IGA Used for Orthotropic Material Crack Growth.

In this paper, we propose a new approach for numerically simulating the growth of cracks in unidirectional composite materials, termed extended isogeometric analysis, evaluating the maximum stress intensity factor and T-stress. To validate our approach, we used a small anisotropic plate with two edge cracks, beginning with formulating the governing equations based on the energy integral method, Stroh's Formula, and the Elastic Law describing the behaviour of anisotropic materials, while considering boundary conditions and initial states. A MATLAB code was developed to solve these equations numerically and to post-process the tensile stress and the stress intensity factor (SIF) in the first mode. The results for the SIF closely match those obtained using the extended finite element method (X-FEM), with a discrepancy of only 0.0021 Pa·m0.5. This finding underscores the credibility of our approach. The extended finite element method has demonstrated robustness in predicting crack propagation in composite materials in recent years, leading to its adoption by several widely used software packages in various industries.

Efficiently assessing the early-age cracking risk of cementitious materials with a mini temperature stress testing machine

Temperature Stress Testing Machine (TSTM) is a universal testing tool for many properties relevant to early-age cracking of cementitious materials. However, the complexity of TSTMs require heavy lab work and thus hinders a more thorough parametric study on a range of cementitious materials. This study presents the development and validation of a Mini-TSTM for efficiently testing the autogenous deformation (AD), viscoelastic properties, and their combined results, the early-age stress (EAS). The setup was validated through systematic tests of EAS, AD, elastic modulus, and creep. Besides, the heating/cooling capability of the setup was examined by tests of coefficient of thermal expansion by temperature cycles. The results of EAS correspond well to that of AD, which qualitatively validates the developed setup. To quantitatively validate the setup, a classical viscoelastic model was built, based on the scenario of a 1-D uniaxial restraint test, to predict the EAS results with the tested AD, elastic modulus, and creep of the same cementitious material as the input. The predicted EAS matched the testing results of Mini-TSTM with good accuracy in 6 different cases. The viscoelastic model also provided quantitative explanations for why variations in early AD do not influence the EAS results. The testing and modelling results together validate the developed Mini-TSTM setup as an efficient tool for studying early-age cracking of cementitious materials. At the end, the potential limitations of the Mini-TSTM are discussed and its applicability for concrete with aggregate size up to 22 mm is demonstrated.

Enhancing machining efficiency in hybrid metal matrix composites through EDM parameter optimization via grey relational analysis

In summary, this study utilized stir casting to develop three distinct aluminium hybrid composites incorporating boron carbide and alumina. These composites were identified as Al7075 + 12%B4C + 6%Al2O3, Al7075 + 6%B4C + 6%Al2O3, and Al7075 + 6%B4C + 12%Al2O3. The novelty of this research lies in the unique combination of materials and their proportions, offering enhanced properties yet unexplored in this specific configuration. The research focused on machining these composite materials using EDM and employed a Taguchi L27 orthogonal array design to plan the machining tests. Key process parameters considered were gap voltage (V), duty cycle (%), current (A), and different percentages of reinforcement content. This aspect adds a significant novel dimension to the study, as the impact of EDM machining parameters on hybrid composites with this particular composition has not been extensively explored previously. The study applied Grey Relational Analysis to evaluate the impact of EDM machining process parameters on critical output responses, specifically material removal rate (MRR) and tool wear rate (TWR). SEM images revealed redeposited molten material, craters, and surface cracks, emphasizing the importance of controlling discharge energy and current levels for desired surface quality. An optimal configuration was identified peak current 2 A, pulse duration 100 μs, duty factor 20%, Composition 2 for improved surface attributes. These qualitative and quantitative insights inform process optimization strategies for EDM machining of aluminium hybrid composites.

Extension of Tsai-Hill failure concept for mixed-mode I/II fracture investigation of orthotropic materials considering T-stress effects

An essential study in the discourse surrounding composite materials is the investigation of their fracture mechanics due to the intricate nature of composite structures, it is necessary to take into account the presence of cracks and faults throughout the design process. In this research, within the context of linear elastic fracture mechanics, the fracture criterion for orthotropic materials under mixed mode II/I loading is presented by expanding the Tsai-Hill failure criterion. According to the crack, it may be aligned or perpendicular to the fibers; the criterion is stated for both cases. Also, the effects of the non-singular stress component on crack growth have been investigated. The mechanical and fracture data of Russian pine wood was extracted through tensile testing and the finite element (FE) method. The validity of the proposed criterion, in comparison with the experimental findings extracted in this research, has been proved. The precision of the previous proposed fracture criteria has been reviewed in comparison with the proposed criterion. The appropriate behavior of the fracture limit curves compared to the experimental data shows the efficacy of the proposed criterion to predict the initiation and propagation of cracks in orthotropic materials.

3D internal crack propagation in brittle solids under non-uniform temperature field: Experimental and numerical simulation

In the field of engineering, elevated temperatures often induce thermal stress. Especially in brittle materials, the propagation of internal cracks under such thermal stress is frequently the primary cause of failure. We utilize the 3D internal laser-engraved crack technology to create internal cracks at various locations inside the half-Brazilian disk specimens. By precisely controlling the power of the heating device, we achieved temperature non-uniformity of the external heat source, resulting in uneven heating of the material’s surface. Simultaneously, by adjusting the geometry of the material and using semi-Brazilian disk specimens, we generated uneven temperature gradients and thermal capacities along different heat conduction paths. The internal cracks at different locations also affect the distribution of thermal stress in the material, further influencing the formation of the non-uniform temperature fields (NUTF). Through the above methods, we successfully conducted physical experiments on the 3D propagation of internal cracks in brittle materials under NUTF. Additionally, we utilized J-integral and the finite element method for thermal analysis to conduct thermo-mechanical coupled 3D numerical simulations. Furthermore, we employed a polarising microscope to analyze the microscopic morphology of 3D internal cracks. Finally, we revealed the mechanism of 3D internal crack propagation under NUTF. The results show that prefabricated internal cracks in group A produce “n-shaped” crack at a height of 15 mm and “u-shaped” cracks at heights of 25 mm and 35 mm. Meanwhile, in group B, the eccentric prefabricated internal cracks ultimately produce “s-shaped” cracks. Through analysis, it was determined that the differential formation occurred due to relative sliding between the upper and lower surfaces of the crack under longitudinal shear stress induced by NUTF. Group A cracks have smooth surfaces, which are mode I-II cracks. group B generated “lance-like” fractures with obvious “anti-binary tree”, which are mixed mode I-II-III cracks. These findings provide an experimental and theoretical basis for the study of 3D internal crack propagation modes under NUTF and deepen the understanding of fracture mechanics in brittle materials.

Automated Surface Crack Detection in Historical Constructions with Various Materials Using Deep Learning-Based YOLO Network

ABSTRACT Cultural heritage (CH) constructions involve the use of diverse masonry materials. Under natural and human influences, masonry materials can undergo various types of damages, with crack damages being most prevalent. Developing a robust model using deep learning (DL) capable of detecting cracks in various CH materials is crucial. In this study, we compared the performance of the DL method You Only Look Once (YOLO) object detection network based on images in different masonry materials (stone, brick, cob, and tile) with that in a modern material (concrete). The dataset used in the study comprised 1213 bricks, 1116 concretes, 955 cobs, 882 stones, and 208 tile images. YOLOv5 architecture, transfer learning, and object detection models were utilized for detecting cracks and compare their performance in different materials. This study represents the first comparison of this kind using an original dataset. The DL model achieved mean average precision values of 94.4%, 93.9%, 92.7%, 87.2%, 83.4%, 81.6%, and 70.3% for concrete; concrete and cob, cob; stone; stone and brick; brick; and tile, respectively. The findings of this study indicate considerable potential for the widespread use of DL techniques in identifying cracks from images across various CH materials and assisting inspection professionals in damage surveys.

Molecular dynamics simulation and experimental study of the material machinability characteristics and crack propagation mechanisms for fused silica double nanoscratches

Fused silica is a high-purity, hard and brittle, difficult-to-machine amorphous material that is widely used in the aerospace, optical instruments, semiconductors, microelectronics and other important fields. However, the influence of adjacent scratches on material deformation and subsurface damage during ultraprecision machining by abrasive particles is not well understood. This study aims to investigate the influence of micron-spaced double scratches on the deformation and damage for fused silica through molecular dynamics simulations and double nanoscratch experiments. An amorphous model of fused silica is created using a two-step simulated annealing method with three-directional periodic boundary conditions. Variable loaded double scratches with different nanometer distances are simulated on the surface of the model, and variable loaded double scratch experiments with varying micron-sized distances between scratches are conducted for fused silica wafers. The simulation and experimental results reveal that the deformation of fused silica during nanoscale double scratching involves both elastic and plastic deformation, and the surface and subsurface damage primarily consists of radial, lateral, Hertzian, and median cracks. The distance between the double scratches affects the load variation, surface and subsurface crack propagation, extent of material damage and removal rate, and surface roughness. Examination of the material subsurfaces using the FIB-TEM sample preparation and testing technique shows that the distance between the double scratches also influences the subsurface damage to the material.

Numerical design of composite material crack arrestors for long-distance gas pipelines

Long-distance gas pipelines are the most common means of transporting natural gas resources. However, the propagation of longitudinal ductile fractures in gas pipelines can result in catastrophic accidents, which must be avoided as much as possible. Crack arrestors can prevent or at least limit the rapid spread of longitudinal cracks by constraining the outwall of pipelines. Carbon fiber-reinforced composite material is one of the best candidates for crack arrestors because of its excellent strength and toughness. However, due to the significant variations in composite material performance, clear design criteria for determining the size of crack arrestors have not yet been developed. This paper presents a comprehensive design procedure for determining the geometrical parameters of composite material crack arrestors based on a finite element model of dynamic pipeline fractures under crack arrestor constraints. Pipeline crack propagation is accomplished using the cohesive zone model, and composite material failure is defined by Hashin criteria. A series of numerical simulations are conducted to assess the ability of crack arrestors to stop a growing crack in a gas pipeline. Ultimately, the minimum wall thickness and axial length of the effective crack arrestor are determined. The proposed design procedure enables the numerical design of composite material crack arrestors and supports safe pipeline operation.

ICDW-YOLO: An Efficient Timber Construction Crack Detection Algorithm.

A robust wood material crack detection algorithm, sensitive to small targets, is indispensable for production and building protection. However, the precise identification and localization of cracks in wooden materials present challenges owing to significant scale variations among cracks and the irregular quality of existing data. In response, we propose a crack detection algorithm tailored to wooden materials, leveraging advancements in the YOLOv8 model, named ICDW-YOLO (improved crack detection for wooden material-YOLO). The ICDW-YOLO model introduces novel designs for the neck network and layer structure, along with an anchor algorithm, which features a dual-layer attention mechanism and dynamic gradient gain characteristics to optimize and enhance the original model. Initially, a new layer structure was crafted using GSConv and GS bottleneck, improving the model's recognition accuracy by maximizing the preservation of hidden channel connections. Subsequently, enhancements to the network are achieved through the gather-distribute mechanism, aimed at augmenting the fusion capability of multi-scale features and introducing a higher-resolution input layer to enhance small target recognition. Empirical results obtained from a customized wooden material crack detection dataset demonstrate the efficacy of the proposed ICDW-YOLO algorithm in effectively detecting targets. Without significant augmentation in model complexity, the mAP50-95 metric attains 79.018%, marking a 1.869% improvement over YOLOv8. Further validation of our algorithm's effectiveness is conducted through experiments on fire and smoke detection datasets, aerial remote sensing image datasets, and the coco128 dataset. The results showcase that ICDW-YOLO achieves a mAP50 of 69.226% and a mAP50-95 of 44.210%, indicating robust generalization and competitiveness vis-à-vis state-of-the-art detectors.

Rock crack initiation triggered by energy digestion

The critical value of rock failure is determined by irreversible deformation (inelastic deformation, damage, and other internal dissipation) processes and external conditions before rock failure. Nevertheless, a thorough explanation of the mechanism causing cracks in rock material has not yet been provided. The strain energy theory is applied in this work to assess the initiation of rock cracks and investigate the relationship between energy digestion and rock strength. Firstly, the uniaxial compression test was conducted on sandstone samples under quasi-static loading conditions and the results of energy evolution, non-linear cumulative digestion, and stored ultimate energy were obtained. Then, a novel algorithm for assessing the initiation of rock cracks has been put forth. The concept of energy digestion index (EDI), which is the ratio of digested energy over the external loading energy, has been developed to characterize the energy absorption capacity of rock material. The result shows a relationship between the maximum growth rate of energy digestion and the increasing rate of variable elasticity modulus and crack initiation. The mechanical characteristics and peak strength of the rock material are negatively correlated with the EDI. By monitoring the digested energy status, an evaluation of the residual strength is introduced based on the relationships, which will initiate further research into in-situ monitoring and failure prediction.

A lift-off suppression algorithm to improve the accuracy of ACFM in crack detection of ferromagnetic materials

Ferromagnetic pressurized components are widely used in petroleum, aerospace and nuclear industries, and their safety and stability are the key to the reliable operation of pressure equipment. In the process of crack detection of ferromagnetic materials by ACFM technology, the changes in lift-off distance caused by coating corrosion and irregularity of the surface will greatly affect the extraction of useful detection signals and reduce the detection accuracy. To overcome the drawback, a lift-off suppression algorithm for crack detection of ACFM is proposed. First, a differential strategy is proposed to discover the time-domain lift-off point of intersection (LOI) from the raw detection signals of ACFM for ferromagnetic materials. Then, the corresponding Bx and Bz signals are demodulated based on the amplitude of the LOI. The simulation and experimental results of the surface cracking of ferromagnetic materials by ACFM show that the novel Bx signals obtained according to the proposed method vary in a small range for different lift-off distances. The lift-off invariant property of the LOI enables the method to improve the interference of the lift-off effect on the conventional Bx signals and to enhance the detection accuracy of crack depth of ferromagnetic materials. The novel Bz signals perform consistently with the conventional Bz signals, and are virtually unaffected by the changes in lift-off distance.

Micro-surface crack inspection with arc-shaped bi-grating laser ultrasound spectroscopy

This paper introduces a novel method for detecting and characterizing initial micro-surface crack in metal structures and materials, with a particular focus on initial fatigue cracks. The proposed technique utilizes an arc-shaped double-spaced grating mask to spatially modulate laser beams, forming an arc-shaped bi-grating laser source on the surface of the specimen. This configuration excites two focused surface waves that have the same focus point but differ in frequency and propagate in opposite directions. By monitoring surface wave signals as they pass through defects and analyzing the relative changes in the amplitude of the frequency components in the spectrum of the reflected and transmitted waves, the initial defects can be identified. Both simulation and experimental results validate the effectiveness of this method in detecting tiny artificial surface cracks, including those shorter than 0.3 mm. In addition, this method has been successfully applied to the detection of actual initial fatigue cracks. Experimental findings demonstrate that focused surface waves enable highly sensitive detection of sub-millimeter fatigue cracks.

A generalized Rivlin–Thomas theory for nucleation of cracks in viscoelastic materials

Recent fracture experiments using the Rivlin–Thomas pure-shear geometry in rubbery viscoelastic materials have suggested nucleation of cracks at a critical stretch, independent on loading rate. For linear material using the rate-independent cohesive model theory, the critical elongation between grips (in a quite general testing geometry) at nucleation is instead monotonically decreasing with rate of up to the square root of the ratio between the instantaneous and relaxed elastic modulus of the material. Shrimali & Lopez-Pamies have made the assumption of a constant critical stretch to develop a theory which contrasts therefore with classical models. However, we further generalize the Rivlin–Thomas theory by assuming an arbitrary relation between nucleation stretch and loading rate, which therefore includes both models as limit cases. We present a simple case of a Double Cantilever Beam (DCB) geometry with linear standard materials, to show an analytical example.

Numerical modeling of shrinkage induced cracking in cementitious materials with three-dimensional simulation and phase-field method

This study investigates the damage and cracking of concrete induced by drying shrinkage. While previous studies have primarily focused on two-dimensional configurations, this research employs three-dimensional phase-field simulations of concrete samples undergoing the drying process. The distribution of inclusions is derived from tomographic images of real samples. The analysis includes the evolution of the damage field and the propagation of cracked zones. The study introduces the formulation and implementation of the numerical method and applies it to elucidate the cracking process resulting from drying shrinkage in concrete composites with various types of inclusions. Numerous numerical simulations were conducted to accurately depict the cracking process induced by drying shrinkage in the samples, evaluating the effects of different inclusions with varying stiffness on the concrete’s cracking phenomenon and comparing the impact of various energy degradation methods on crack simulation. Moreover, the study analyzes the influence of the spatial distribution of the inclusions in the three-dimensional simulation and compares the numerical results with experimental observations.