Damage Modes Research Articles

Low-velocity impact response of a novel bionic turtle shell back armor sandwich structure

Turtle shell dorsal armor comprises a typical lightweight laminated composite structural material, which is constructed by cuticle and bone layer, with high specific strength and high fracture toughness. The red-eared slider turtle, a widely studied species of turtle, has a rich hierarchical structure. In this paper, according to the structure of red-eared slider turtle shell dorsal armor and the characteristics of each layer as well as by taking into full consideration of the fine structure of cuticle and ventral cortex of the shell dorsal armor, two bionic foam silicone rubber sandwich structures of different levels were prepared. Low-velocity impact test and finite element simulation analysis were carried out to study the mechanical response, specific absorbed energy, and damage mode of the two bionic sandwich structures under different impact energies, hammerhead shapes, and hammerhead radii. By comparing the mechanical response and specific energy absorption with the existing ordinary sandwich structure, we find that the two sandwich structures modeled in this paper exhibit better impact resistance than the ordinary sandwich structure. The shape and radius of the hammer head have a significant effect on the energy absorption characteristics of the sandwich structure, and the damage mode under the hemispherical hammer head is mainly dominated by tensile tearing damage, while the damage mode of extrusion fracture is more obvious under the impact of the conical hammer head. With the increase of the radius of the hammer head, the absorbed energy and specific energy absorption of the sandwich structure increases.

The transplant rejection response involves neutrophil and macrophage adhesion-mediated trogocytosis and is regulated by NFATc3

The anti-foreign tissue (transplant rejection) response, mediated by the immune system, has been the biggest obstacle to successful organ transplantation. There are still many enigmas regarding this process and some aspects of the underlying mechanisms driving the immune response against foreign tissues remain poorly understood. Here, we found that a large number of neutrophils and macrophages were attached to the graft during skin transplantation. Furthermore, both types of cells could autonomously adhere to and damage neonatal rat cardiomyocyte mass (NRCM) in vitro. We have demonstrated that Complement C3 and the receptor CR3 participated in neutrophils/macrophages-mediated adhesion and damage this foreign tissue (NRCM or skin grafts). We have provided direct evidence that the damage to these tissues occurs by a process referred to as trogocytosis, a damage mode that has never previously been reported to directly destroy grafts. We further demonstrated that this process can be regulated by NFAT, in particular, NFATc3. This study not only enriches an understanding of host-donor interaction in transplant rejection, but also provides new avenues for exploring the development of novel immunosuppressive drugs which prevent rejection during transplant therapy.

Finite element analysis of the effect of crystallinity on low-velocity impact performance of poly (aryl ether ketone) composites

A new copolymer poly (aryl ether ketone) (PAEK) resin was used as a matrix to prepare carbon fiber (CF) reinforced composites with different crystallinity by slow and rapid cooling. Finite element (FE) models were constructed according to micromechanics, and the progressive damage method based on the three-dimensional (3D) PUCK damage criterion was used to analyze the low-velocity impact damage mechanism of the CF/PAEK composites. The experimental results show that the composites with lower crystallinity have higher impact resistance and smaller impact damage area. In the FE simulated impact damage evolution process, the damage of the impact side is mainly compressive damage of the matrix and the fiber, while that of the impact back side is mainly tensile damage of the matrix. Moreover, due to their high interfacial bonding strength, the main damage modes of rapid-cooled composites with lower crystallinity are severe matrix tensile damage and slight fiber tensile damage in the impact back side. However, the slow-cooled composites show more delamination. In addition, for the damage evolution of compression after impact, the FE simulation shows that fiber compression, matrix tension, matrix compression and delamination expand along the transverse to the compression load. Furthermore, the in-plane damage expansion occurs earlier in the rapid-cooled composites, which also have graver interlaminar damage.

Rotational Stiffness Investigation and Parametric Analysis of a Novel Assembled Joint in Lattice Shells

Although there are currently many types of lattice shell joints with different characteristics, assessing the flexural capacity of lattice shell joints is always a great challenge. In this paper, a fan-shaped assembled joint and a welded joint for comparison were subjected to bending tests to investigate the flexural behavior and rotational stiffness of the assembled joint. The strain distribution, load–displacement curve, moment–rotation curve, and damage modes of key parts were analyzed to determine the vulnerable parts of the joints. Our test results show that, with an initial rotational stiffness of about one third of that of the welded joint, the assembled joint specimen exhibits the obvious characteristics of a semi-rigid joint. The finite element analysis results were in good agreement with the experimental results. The results of our parametric analysis show that the rotational stiffness and ultimate moment of the assembled joint increase with increases in the spacing of the bolts and the number of bolts. The performance of the high-strength bolts had a significant influence on the flexural stiffness of the assembled joints. The spacing of the bolts and the number of bolts for the assembled joint are suggested to be greater than the height of the member section and more than three, respectively. The proposed theoretical formula can approximately simulate the initial rotational stiffness of the joint. More in-depth investigations are required in the future for assessing the mechanical behavior of FSA joints subjected to combined bending–compression loads.

Research on damage evolution mechanisms under compressive and tensile tests of plain weave SiCf/SiC composites using in situ X-ray CT

Abstract To investigate the damage evolution and failure mechanisms of fiber-reinforced composite materials under complex conditions, this study conducted in situ X-ray computed tomography (CT) compression and tensile tests on plain weave two-dimensional woven SiC/SiC composite materials. The obtained CT in situ image data captured the behavior of materials during loading and after failure. Using the image reconstruction of CT data, the actual microstructure and damage evolution of the material under six consecutive loading levels were accurately revealed. Three-dimensional visualization models of the composite material were established using image processing software to analyze the damage evolution under compression and tension, and the failure mechanisms were compared. The results showed that the compression and tension failure mechanisms of SiC/SiC composite materials were similar, with the transverse cracking of the matrix being the first mode of damage, followed by delamination between layers and longitudinal matrix cracking of fiber bundles. Specifically, in terms of compression failure, the strength of the fiber bundle itself has a greater influence, and fiber fracture is the main cause of ultimate material failure. On the other hand, the primary cause of tensile failure is the presence of porosity defects generated during material fabrication. Consequently, the tensile material fails earlier and can withstand lower loads.

Characterizing damage patterns and evolution in Multi-Hole GLARE laminates under tensile load via integrated AE and DIC techniques

In this study, a combined numerical simulation and experimental approach was employed to investigate the mechanical behavior and damage mechanisms of multi-hole GLARE laminates with varying off-axis angles under static tensile loading. Digital image correlation techniques were used to monitor changes in the strain field near the open hole region. The acoustic emission parameters were used to recognize damage patterns and evolution. The progressive damage evolution processes of various GLARE components were investigated using a 3D finite element model (FEM) based on the Hashin strain criterion and cohesive zone modeling. The results indicate that GLARE experiences damage modes such as aluminum alloy breakage, matrix cracking, fiber/matrix debonding, delamination, and fiber breakage. The frequency ranges associated with the five damage modes are as follows: [0–50 kHz], [100–175 kHz], [175–220 kHz], [220–300 kHz], and [300–400 kHz]. The dominant damage mode throughout the tensile phase is matrix cracking, but the ultimate failure of the specimen is determined by aluminum alloy breakage and fiber breakage. In addition, increasing the off-axis angle of the multi-hole leads to a shift in the final failure mode from tensile-dominated to tensile shear-dominated, while also slightly reducing the overall initial stiffness of the specimen.

Study on the low‐velocity impact response of the repaired carbon fiber‐reinforced plastic laminates: Bonded versus scarf repair

AbstractExperimental and simulation studies of the response of bonded‐repaired and scarf‐repaired carbon fiber‐reinforced polymer laminates under the low‐velocity impact (LVI) were carried out. The center‐ and eccentric‐repaired laminates were prepared for experiments. The experimental results verify the validity of the finite element model in terms of the force response and different modes of damage. The effect of impact energy and the variation of repair position on the impact response of the laminate was investigated based on the benchmark model. The bonded repaired laminates have higher peak impact forces and lower energy absorption than the scarf‐repaired laminates; the damage of the bonded repaired laminate is mainly intralaminar damage and delamination damage of the patch; the adhesive of the scarf repair is prone to failure, leading to poor impact resistance. The relative positions of the impact point and the patch determine the impact response of a laminate repaired at different locations. The impact point at the edge of the patch has a detrimental effect on the repaired laminate. The impact point vertical position near the bottom hole edge of the parent plate in the scarf‐repaired case harms the impact resistance of the laminate.Highlights The low‐velocity impact response of bonded versus scarf‐repaired laminates was compared. Bonding repairs provide superior impact resistance to damaged laminates than scarf repairs. The risk of adhesive failure in scarf repair is significantly higher than in bonded repair. The relative position of the damage (repair) to the impact point critically influences the impact response. Impacts near the edge of the initial damage hole can cause severe damage to scarf‐repaired laminates.

Tensile damage self-monitoring of carbon fiber/epoxy 3D braided composites with electrical resistance method

Electrical resistance method (ERM) has great potential in self-monitoring structural health for carbon fiber composites due to carbon fiber conductivity. Here we report the correlation between tensile damages and electrical resistance changes of 3D braided composites subjected to monotonic and cyclic tensile loading. Experimental and finite element methods were utilized to reveal the electrical resistance variation and damage evolution under monotonic tensile loading. We found that the electrical resistance changes under monotonic tensile loading could be segmented into three regions, namely, linear growth region, transition region and exponential growth region, corresponding to different damage modes. Synchronous and reversible electrical resistance changes at a low strain level (12% of fracture strain) cyclic loading proved the stability and repeatability of self-monitoring capability. While, the abrupt electrical resistance changes at a high strain level (75% of fracture strain), indicated catastrophic failure of the composites. Such kind of correspondence between tensile damages and electrical resistance change under tensile loading enables ERM to achieve the in-situ structural health monitoring of 3D braided composites without relying on external sensors.

Thermal-mechanical coupling in drilling high-performance CFRP: Scale-span modeling and experimental validation

The focus of this study is to explore a scale-span thermal–mechanical coupling method for predicting the dynamic mechanical progressive failure behaviors and predicting thermal damage during the drilling process of Carbon Fiber Reinforced Plastic (CFRP). Initially, a thermal conduction constitutive model of drilling CFRP was developed based on the proposed thermal distribution ratio calculation method. Meanwhile, a dynamic span-scale progressive damage constitutive model for CFRP, incorporating modified micromechanics failure criterion with bilinear damage evolution laws, was proposed, and a bilinear cohesive model, including three damage modes, is employed to simulate interlaminar delamination. Subsequently, a user-defined material subroutine VUMAT was implemented on the ABAQUS/Explicit platform to simulate the thermal–mechanical coupling behaviors of drilling T700S-12K/YPH-23 CFRP using a twist drill bit. Finally, a comprehensive information monitoring platform for CFRP drilling experiments was established to validate the accuracy of the simulation results by considering drilling temperature, thrust force, and hole-wall morphology. The results demonstrate excellent agreement between the established thermal–mechanical coupling span-scale model and the experimental data. Furthermore, the simulation effectively captures the various damage behaviors and thermal conduction phenomenon, that occur during the intact drilling process.

Dynamic response and energy absorption of aluminum foam sandwich under low-velocity impact

Aluminum foam sandwich (AFS) are widely used in energy absorption because of their light weight and excellent energy dissipation capabilities. In this study, energy absorption, deformation modes, and dynamic response of AFS under low-velocity impact are investigated using experimental and numerical approaches. Dynamic impact tests are performed utilizing a drop hammer impact test on AFS with different core densities, face sheet thicknesses, and impact energy. Based on the 3D Voronoi foam model, a full-scale finite element model (FEM) is developed to simulate the mechanical response of AFS under low-velocity impact, and verified by experimental results. The contributions of different AFS components to energy absorption at various impact velocities are further analyzed, along with the effects of face sheet thickness distribution. The results indicate that the performance of the face sheet and core and the impact energy have a remarkable influence on the low-velocity impact response and damage modes of the AFS. The type of structure configuration, “thin top and thick bottom,” enables the AFS better play its energy absorption while ensuring the overall light weight. This study provides a reference for the design and optimization of the face sheet thickness of AFS when it is used as an energy-absorbing member and improves their design efficiency.

Damaging effects of earth fissures on high-speed railway bridges: Insights from physical experiments and numerical simulations

Earth fissure is a typical geological phenomenon related to surface rupture. With the rapid development of high-speed railways (HSRs) in recent years, earth fissure disasters pose a serious challenge to the construction and subsequent safe operation of HSRs. Here, we aimed to investigate the damaging effect of earth fissure activities on HSR bridges spanning the earth fissure zone. Taking the Shuixiu earth fissure of Da-Xi HSR as the prototype, a large-scale laboratory physical experiment with geometric model scale 1:20 was carried out to reveal the deformation and failure characteristics of HSR bridges under earth fissure activity. In addition, the damage influence range and damage mode of HSR bridges crossing earth fissures are further determined by numerical simulation. The results show that the earth fissure activity mainly has damaging effects on the structure of HSR bridges near the earth fissure zone, which is mainly manifested as uneven settlement of reinforced concrete pile, bending deformation of box girder and cracking of track slab. The width of the influence zone of Shuixiu earth fissure on the HSR bridge is determined to be 32.4 m. Moreover, the pile-end pressure on the hanging wall side of the earth fissure gradually decreases with the earth fissure activity, while that on the footwall side increases significantly. Most importantly, it is emphasized that a practical long-term monitoring and early warning system of HSR bridge structures should be established to prevent the hazards caused by earth fissure activities. The above understanding will benefit to the design, construction and operation of Da-Xi HSR line, and also provide significant guidance and reference to the construction of HSR bridges in other earth fissure prone areas.

Enhancing Structural Health Monitoring with Acoustic Emission Sensors: A Case Study on Composites under Cyclic Loading.

This study conducts an in-depth analysis of the failure behavior of woven GFRP under cyclic loading, leveraging AE sensors for monitoring damage progression. Utilizing destructive testing and AE methods, we observed the GFRP's response to varied stress conditions. Key findings include identifying distinct failure modes of GFRP and the effectiveness of AE sensors in detecting broadband frequency signals indicative of crack initiation and growth. Notably, the Felicity effect was observed in AE signal patterns, marking a significant characteristic of composite materials. This study introduces the Ibe-value, based on statistical parameters, to effectively track crack development from inception to growth. The Ibe-values potential for assessing structural integrity in composite materials is highlighted, with a particular focus on its variation with propagation distance and frequency-dependent attenuation. Our research reveals challenges in measuring different damage modes across frequency ranges and distances. The effectiveness of Ibe-values, combined with the challenges of propagation distance, underscores the need for further investigation. Future research aims to refine assessment metrics and improve crack evaluation methods in composite materials, contributing to the field's advancement.

Study of Low-Velocity Impact Behavior of Hybrid Fiber-Reinforced Metal Laminates.

In this paper, the low-velocity impact behavior and damage modes of carbon/glass-hybrid fiber-reinforced magnesium alloy laminates (FMLs-H) and pure carbon-fiber-reinforced magnesium alloy laminates (FMLs-C) are investigated using experimental, theoretical modeling, and numerical simulation methods. Low-velocity impact tests were conducted at incident energies of 20 J, 40 J, and 60 J using a drop-weight impact tester, and the load-displacement curves and energy-time curves of the FMLs were recorded and plotted. The results showed that compared with FMLs-C, the stiffness of FMLs-H was slightly reduced, but the peak load and energy absorption were both greatly improved. Finally, a finite element model based on the Abaqus-VUMAT subroutine was developed to simulate the experimental results, and the damage modes of the metal layer, fiber layer, and interlayer were observed and analyzed. The experimental results are in good agreement with the finite element analysis results. The damage mechanisms of two kinds of FMLs under low-velocity impacts are discussed, providing a reference for the design and application of laminates.

Assessment of tensile damage mechanism of open‐hole GLARE laminates based on acoustic emission and digital image correlation techniques

AbstractGLARE laminates have emerged as a favored option for the construction of fuselage and wing skins in large airliners, owing to their exceptional mechanical characteristics. Nevertheless, the incorporation of open‐hole designs poses a challenge as it disrupts the continuity of the laminates, resulting in stress concentration and subsequent damage. To investigate the impact of various layup orientations and hole sizes on the tensile properties of open‐hole GLARE laminates, this study conducted axial tensile tests. Additionally, the tensile damage process was monitored using DIC and AE techniques, enabling the identification of damage patterns and the analysis of their evolution. The results demonstrate a noteworthy decline in the ultimate strength and failure strain as the size of the opening increases. Moreover, the retention rate of failure strain displays a marked sensitivity to the layup orientation. In conjunction with observations made through DIC and SEM, the k‐means++ algorithm successfully clustered peak frequencies, thereby revealing distinct damage patterns and their corresponding frequency ranges as aluminum alloy damage [0–90 kHz], matrix cracking [104–174 kHz], fiber/matrix debonding [175–224 kHz], interlaminar delamination [234–300 kHz], and fiber fracture [304–469 kHz]. AE cumulative counts were utilized to evaluate the progression of individual damage modes. The results emphasize that matrix cracking demonstrates the most substantial cumulative counts, whereas damage to the fibers and aluminum alloy noticeably affects the load‐carrying capability of the laminate. Furthermore, the fibers/matrix debonding and interlaminar delamination, exhibit heightened susceptibility to layup orientation and hole size.Highlights The tensile performance evaluation of GLARE laminates with open‐hole varying layup orientations and hole sizes was investigated. A combination of DIC and AE techniques was used to monitor the tensile damage process. Damage pattern was identified based on the observations of DIC and SEM. An approach based on Pearson's correlation coefficient and k‐means++ clustering algorithm was used for damage pattern recognition in GLARE laminates. AE cumulative counts were employed to assess the evolution of each damage mode.

Seismic performance test of bolted U-shaped steel assembled shear wall

A new type of assembled shear wall with energy dissipation and seismic damping function was proposed for the damage mechanism and weak points of typical assembled shear wall structures under strong earthquake effects. Using a combination of model tests and numerical simulations, four new assembled shear wall specimens with a scaling ratio of 1:1.54 and a shear-to-span ratio of 1.52 were designed and fabricated. The corresponding seismic performance tests were carried out. The effects of the number of bolts, axial compression ratio and longitudinal reinforcement ratio of edge members on the damage mode, hysteresis performance, bearing capacity, deformation performance, stiffness degradation and energy dissipation capacity of the specimens were systematically analyzed. The test results show that the damage mode of each specimen is basically the same as that of the common shear wall with the same shear-to-span ratio, which is bending-shear type damage However, the new assembled shear wall has more excellent hysteresis performance and energy dissipation capacity. Friction between the steel plates of the horizontal joints consumes a lot of energy and further inhibits wall damage. The energy dissipation value at the point of destruction of the wall is significantly higher than that of ordinary cast-in-place wall. The energy consumption value of the specimen PFSW is about 6 times that of a cast-in-place wall. When the number of bolts decreases, hysteresis performance decreases and wall deformation increases. The decrease of axial pressure ratio or longitudinal reinforcement ratio of edge members leads to the decrease of bearing capacity and increase of ultimate displacement of shear wall. Finally, the finite element model of the corresponding specimens was established by ABAQUS software. The simulation results were in good agreement with the test results, which indicated the correctness of the proposed model. And it can be applied to the analysis of new assembled shear walls.

Investigation of low‐velocity impact behaviors of polymer composites reinforced with different natural fiber fabrics

AbstractIn this study, the low‐velocity impact behavior of polymer composites produced separately by reinforcing with natural fiber fabrics such as linen, hemp, jute, and bamboo was investigated for three different thicknesses (2.5 ± 0.1, 5 ± 0.1, and 7.5 ± 0.1 mm) by vacuum infusion method. Low‐speed impact tests were performed at three different energy levels of 20, 35, and 50 J by striking the geometric center of the composite samples with a hemispherical impactor. The impact strengths of the composites were evaluated by comparing the contact force‐displacement, contact force‐time graphs, and impacted sample images. The types of damage occurring in each experimental sample were evaluated. Increasing the thickness had a positive effect on the impact strength of the composites. When the thickness of the composite plates increases, the peak contact force increases, the displacement decreases and capacity to carry impact loads increases due to the increased stiffness of the composites. However, except for the damaged 5 ± 0.1 mm thick bamboo and jute samples, the energy absorption ability of the samples generally increases with the increase of impact energy. Bamboo, jute, and hemp reinforced samples with a thickness of 2.5 ± 0.1 mm were completely perforated under 20 J impact energy. Also, it was determined that linen reinforced composites have the best impact resistance among natural fabric reinforced composites.Highlights Low‐velocity impact behavior of natural fiber fabric reinforced PMCs was studied. The impact strength of PMCs increased as the thickness increased. Linen reinforced composites had the best impact resistance. The main damage modes are matrix cracking and delamination at low‐impact energies. As the impact energy increases, perforation and fiber breakage or separations are observed.

Physical and Mechanical Properties and Damage Mechanism of Sandstone at High Temperatures

The physical and mechanical properties of rocks change significantly after being subjected to high temperatures, which poses safety hazards to underground projects such as coal underground gasification. In order to investigate the effect of temperature on the macroscopic and microscopic properties of rocks, this paper has taken sandstone as the research object and conducted uniaxial compression tests on sandstone specimens at different temperatures (20–1000 °C) and different heating rates (5–30 °C/min). At the same time, the acoustic emission (AE) test system was used to observe the acoustic emission characteristics of the rock damage process, and the microstructural changes after high temperature were analyzed with the help of a scanning electron microscope (SEM). The test results show that the effect of temperature on sandstone is mainly divided into three stages: Stage I (20–500 °C) is the strengthening zone, the evaporation of water and the contraction of primary fissures, and sandstone densification is enhanced. In particular, the compressive strength and elastic modulus increase, the macroscopic damage mode is dominated by shear damage, and the fracture micromorphology is mainly brittle fracture. Stage II (500–600 °C) is the transition zone, 500 °C is the threshold temperature for the compressive strength and modulus of elasticity, and the damage mode changes from shear to cleavage damage, and the sandstone undergoes brittle–ductile transition in this temperature interval. Stage III is the physicochemical deterioration stage. The changes in the physical and chemical properties make the sandstone compressive strength and modulus of elasticity continue to decline, the macroscopic damage mode is mainly dominated by cleavage damage, and the fracture microscopic morphology is of a more toughness fracture. The effect of different heating rates on the mechanical properties of sandstone was further studied, and it was found that the mechanical properties of the rock further deteriorated under higher heating rates.

A novel constitutive model for the concrete face slab of a CFRD and the numerical simulation of its seismic behavior

AbstractWith the rising demand for hydropower energy, China has built numerous concrete‐faced rockfill dams (CFRDs) with considerable heights. The potential failure pattern of the concrete face slab is the key to the safety of the anti‐seepage system. However, engineers commonly employ an elastic relationship to model the concrete face slab for simplicity, which induces inaccurate estimation of the face slab's stresses, especially under seismic conditions. On the other hand, the existing elastoplastic models for concrete often require complex derivation of the necessary formula and thus are difficult to apply to finite element (FE) analysis. Under the comprehensive action of the upstream water, supporting rockfill, and surrounding mountains, the damage mode of the quasibrittle concrete slab is not a tensile mode but a compressive shearing mode. Therefore, it is urgent to establish an elastoplastic model that is adapted to the actual working environment of the slab and easy to apply. To this end, a generalized plasticity model is developed within the phenomenology modeling concept, which is validated at both the experimental level and application level. The results indicate that three‐dimensional stress‒strain behaviors, such as the strength nonlinearity, prepeak dilatancy, and postpeak softening, can be effectively captured. In addition, the seismic characteristics of a 200‐m high CFRD are numerically investigated. The findings show that the slab's stress accumulation can be reflected during the earthquake and that the final stress distribution conforms to the field observations. The conclusions of this paper provide a reliable reference for the seismic resistance design of CFRDs.

The thickness effect of rubbery nanofibrous mat on modes I–II fracture mechanism of composite laminates

This study investigates the effect of the interleaving nanofibers, made of NBR/PCL blend, on the interlaminar fracture toughness of carbon/epoxy laminates. Different nanomat thicknesses, ranging from 20 μm to 120 μm, were tested at Mode-I and mode-II and results were compared to the non-modified laminates. Acoustic Emission (AE) technique was used to assess the influence of interleaving nanofibers on dominant damage modes of the specimens, i.e., matrix cracking, fiber/matrix debonding, and fiber breakage. Moreover, the damage mechanism and the nanofiber toughening contribution were investigated by means of crack path and surface analysis. Results indicated that the optimum nanomat thickness for mode-I is 40 μm (GI,R = + 333 %), while for mode-II is 20 μm (GII,R = + 43 %). The study also confirmed by AE the significant impact of nanofibers on various damage modes, especially during mode-I loading.

Multi-scale numerical analysis of damage modes in 3D stitched composites

3D stitching composites have already sparked great interest among researchers due to their ability to effectively improve the weak inter-laminar properties in composite laminate structures. In this paper, a systematic analysis on mode I fracture and low-velocity impact damage of 3D stitched composites has been demonstrated via a novel multi-scale model which is capable of efficiently characterizing the distribution of stitching threads. The damage modes including intra-laminar cracks, inter-laminar delamination and stitching thread breakage are identified via a combination of continuum damage mechanics (CDM), cohesive zone model (CZM) element and truss element. This novel modeling approach enables an association analysis of stitching thread breakage status and other damage morphologies, promoting an in-depth study for the influence of stitching thread on damage evolution, which has not been reported yet. The introduction of stitching threads effectively suppresses the sharp propagation of delamination after peak load, avoiding the catastrophic destruction during mode I fracture. The mechanism of damage inhibition from stitching threads is well revealed by the numerical simulation. Furthermore, the capability of 3D stitched composites to resist low-velocity impact damage is effectively predicted in terms of peak load, damage area and damage morphology, finding that the increase of stitch density inhibits the propagation of matrix crack, as well as delamination. The universality and the expanding application of the model are also evaluated, demonstrating that it is suitable for various 3D stitched composites.