Fracture Analysis of Highly Flexible Adhesives: Cohesive Zone Modelling across a Wide Spectrum of Temperatures and Strain Rates.

This study focuses on the prediction of the fracture mechanics behaviour of a highly flexible adhesive (with a tensile elongation of 90%), since this type of adhesive is becoming widely used in automotive structures due to their high elongation at break and damping capacity. Despite their extensive applications, the understanding of their fracture mechanics behaviour under varying loading rates and temperatures remains limited in the literature. In addition, current prediction models are also unable to accurately predict their behaviour due to the complex failure mechanism that such bonded joints have. This study aims to determine whether a simple triangular cohesive zone model (CZM), which predefines the crack path, can reproduce the load-displacement curves of adhesives under various temperatures and strain rates. To achieve this, a calibrated CZM is used, adapting the model for reference joints and then validating it with independent test results conducted in a wide range of loading and environmental conditions. The tests were performed at speeds between 0.2 and 6000 mm/min and at three different temperatures ranging from -30 °C to 60 °C. Mode I fracture toughness was measured using the DCB (double cantilever beam) specimens. Using a simple triangular CZM may not be optimal for predicting the mechanical response of highly flexible adhesives with complex failure mechanisms and multiple crack paths. However, by correctly adjusting the cohesive zone properties for a limited set of reference conditions, it is possible to accurately predict the mechanical response of these joints across various test speeds and temperatures, significantly reducing costs and effort.

Investigation of Axial Tensile Fracture Performance of Recycled Brick Coarse Aggregate Concrete Using a Cohesion Model.

Currently, microscopic research on the tensile fracture properties of recycled brick coarse aggregate concrete has mainly adopted microscopy techniques, which can clearly observe the actual damage situations of each phase material but are unable to individually analyze the effect of a specific material factor on the tensile properties of recycled concrete. This brings much uncertainty to the practical application of recycled concrete. Therefore, this study proposes a cohesive zone model (CZM) for simulating the tensile fracture of recycled brick coarse aggregate (RBCA) concrete. To this end, the study explores the effects of various critical factors on the fracture mode and bearing capacity of recycled brick aggregate concrete, including the replacement rate of recycled brick coarse aggregate, pore structure, interfacial transition zone (ITZ) strength, mortar strength, and volume fraction of brick aggregate. The results indicate that, when the minor to major axis ratio of elliptical pores is 0.5 ≤ K < 1, the following order of influence can be observed: random convex polygonal pores, circular pores, and elliptical pores. Moreover, excessively strengthening the ITZ and mortar does not significantly enhance the tensile performance of RBCA concrete. The distribution location of aggregate has a significant impact on the crack shape of recycled concrete, as does the pore structure, due to their randomness. Therefore, this article also discusses these. These findings contribute to a comprehensive understanding of the tensile properties of recycled brick coarse aggregate and provide insights into optimizing its behavior.

A Study of the Effects of Moisture on Composite-to-Metal Double-Lap Shear Joints.

This work investigated the effects of moisture absorption treatment on composite-to-metal double-lap shear joints (DLSJs) bonded with epoxy adhesive film through experiments and simulations. The composite-to-metal DLSJ can be divided into five parts (the interface between the composite and adhesive, the interface between the adhesive and metal, the composite adherend, the metal adherend, and the adhesive layer). First, the wet-dependent properties of the adhesive and interfaces were obtained through adhesive tensile tests and GC tests, which showed that the properties of the adhesive and interfaces were significantly affected by the moist environment. Then, tensile tests of the composite-to-metal double-lap shear joints were carried out in dry and wet environments. Finally, based on the experimental investigations, a finite element (FE) model that considered cohesive damage was established for simulating damage evolution and predicting the failure loads and failure modes of the DLSJs. The results of both the experimental and numerical tests show that the DLSJ failure load decreases significantly after immersion in 95 °C water, and the major failure mode transfers from adhesive failure to interface failure. The research results provide a theoretical basis or basic data for the structural design of adhesively bonded composite-to-metal.

Finite-element modelling of interactions of needle with tympanic membrane and middle ear.

The tympanic membrane (TM) is one of the most common routes to access the middle ear and inner ear for the treatment of hearing and balance pathologies. Since the TM is a soft thin biological tissue with small dimensions, using needles seems to be among the most practical interventional approaches. In this study, we proposed a finite-element (FE) analysis of needle-TM interactions that combines a 3D model of the TM and other main middle-ear structures in gerbil, and a 2D model of needle insertion into the TM based on the cohesive zone method (CZM). The TM was modelled using a 1st-order Ogden hyperelastic material and its properties were obtained by fitting to the experimental force-displacement plots of large deformation in the TM under needle indentation. The cohesive parameters were also acquired by calibrating the puncture force against the experimental data of needle insertion into the TM. These FE models were then used to obtain the deformation behaviour of the TM and other middle-ear structures due to the insertion force applied at different locations on the TM. Moreover, we investigated the effect of the TM thickness, the geometry of the needle (i.e., diameter and tip angle), and needle material on the insertion of needles into the TM. We also studied the penetration success of deformable needles.

Modeling Interface Damage with Random Interface Strength on Asphalt Concrete Impervious Facings.

Asphalt concrete impervious facings, widely adopted as the impervious structures for rockfill dams and upper reservoirs in pumped storage power stations, typically have a multilayer structure with a thin sealing layer, a thick impervious layer, and a thick leveling bonding layer. The properties of the interfaces between these layers are crucial for the overall performance of the facings. This paper develops a model to investigate the complex interface damage behavior of the facing under static water pressure and gravity. The model considers two damage origins: one is the interface adhesion-decohesion damage, which is described by the cohesive zone model (CZM) combined with the Weibull-type random interface strength distribution, and the other is the bulk damage of each layer, described by Mazars' model. Primarily, a comparison between numerical simulation and indoor direct shear tests validates the reliability of the CZM for the asphalt concrete layer interface. Then, the damage distribution of the two interfaces is simulated, and the characteristics of the interface stress are analyzed in detail. The interface shear stresses of the ogee sections, which have different curvatures, all show an interesting oscillation between the thin sealing layer and the impervious layer, and the interface damage at this interface exhibits high heterogeneity. Furthermore, tension stress exists in the local zones of the ogee section, and the damage in this section is significantly greater than in other parts of the facings.

On the Importance of Including Cohesive Zone Models in Modelling Mixed-Mode Aneurysm Rupture.

INTRODUCTION: The precise mechanism of rupture in abdominal aortic aneurysms (AAAs) has not yet been uncovered. The phenomenological failure criterion of the coefficient of proportionality between von Mises stress and tissue strength does not account for any mechanistic foundation of tissue fracture. Experimental studies have shown that arterial failure is a stepwise process of fibrous delamination (mode II) and kinking (mode I) between layers. Such a mechanism has not previously been considered for AAA rupture. METHODS: In the current study we consider both von Mises stress in the wall, in addition to interlayer tractions and delamination using cohesive zone models. Firstly, we present a parametric investigation of the influence of a range of AAA anatomical features on the likelihood of elevated interlayer traction and delamination. RESULTS: We observe in several cases that the location of peak von Mises stress and tangential traction coincide. Our simulations also reveal however, that peak von Mises and intramural tractions are not coincident for aneurysms with Length/Radius less than 2 (short high-curvature aneurysms) and for aneurysms with symmetric intraluminal thrombus (ILT). For an aneurysm with (L/R = 2.0), the peak σ vm moves slightly towards the origin while the peak T t is near the peak bulge with a separation distance of ~ 17 mm. Additionally, we present three patient-specific AAA models derived directly from CT scans, which also illustrate that the location of von Mises stress does not correlate with the point of interlayer delamination. CONCLUSION: This study suggests that incorporating cohesive zone models into clinical based FE analyses may capture a greater proportion of ruptures in-silico.

A hybrid finite-discrete element method for modelling cracking processes in sandy mudstone containing a single edge-flaw under cyclic dynamic loading.

Rock mass deformation and failure are macroscopic manifestations of crack initiation, propagation, and coalescence. However, simulating the transition of rocks from continuous to discontinuous media under cyclic dynamic loading remains challenging. This study proposes a hybrid finite-discrete element method (HFDEM) to model crack propagation, incorporating a frequency-dependent cohesive-zone model. The mechanical properties of standard sandy mudstone under quasi-static and cyclic dynamic loading were simulated using HFDEM, and the method's reliability was verified through experimental comparison. The comparative analysis demonstrates that HFDEM successfully captures crack interaction mechanisms and accurately simulates the overall failure behavior of specimens. Additionally, the effects of pre-existing flaw inclination angle and dynamic loading frequency on rock failure mechanisms were investigated. The numerical results reveal that rock samples exhibit significantly higher compressive strength under dynamic loading compared to quasi-static loading, with compressive strength increasing with higher cyclic dynamic load frequencies. Furthermore, by analyzing the strength characteristics, crack propagation, and failure modes of the samples, insights into the failure mechanisms of rocks under different frequency loads were obtained. This study provides valuable insights into crack development and failure of rocks under seismic loads, offering guidance for engineering practices.

Effect of Post-Cured through Thickness Reinforcement on Disbonding Behavior in Skin-Stringer Configuration.

An experimental investigation of interlaminar toughness for post-cured through-thickness reinforcement (PTTR) skin-stringer sub-element is presented. The improvement in the crack resistance capability of skin-stringer samples was shown through experimental testing and finite element analysis (FEA) modeling. The performance of PTTR was evaluated on a pristine and initial-disbond of the skin-stringer specimen. A macro-scale pin-spring modeling approach was employed in FEA using a non-linear spring to capture the pin failure under the mixed-mode load. The experimental results showed a 15.5% and 20.9% increase in strength for the pristine-PTTR and initial-disbond PTTR specimens, respectively. The modeling approach accurately represents the overall structural response of PTTR laminate, including stiffness, adhesive strength, crack extension scenarios and progressive pin failure modes. This modeling approach can be beneficial for designing damage-tolerant structures by exploring various PTTR arrangements for achieving improved structural responses.

Influence of Failure-Load Prediction in Composite Single-Lap Joints with Brittle and Ductile Adhesives Using Different Progressive-Damage Techniques.

The usage of adhesively bonded joints, such as single-lap and double-lap joints, is increasing rapidly in aerospace composite structures as a popular alternative to bolts and rivets. Compared to the conventional joining methods such as fastening and riveting, adhesive-bonding technology better prevents damage to composite structures due to the smooth configuration and the mitigation of stress concentration around holes. In this work, the built-in progressive-damage-modeling techniques in Abaqus, including the cohesive zone model (CZM) and the virtual crack closure technique (VCCT), are used to predict the strength and progressive failure of composite single-lap joints subjected to tensile loading. Modeling of an adhesive layer by using a zero/non-zero-thickness cohesive element, cohesive surface, and VCCT is investigated, as is the effect of brittle and ductile adhesives. Two-dimensional finite-element models with different damage-modeling strategies are performed in this study. The failure-load predictions are compared with the experimental results obtained from the literature. For the ductile adhesive, the predicted failure loads using a zero/non-zero-thickness cohesive elements and a cohesive surface are all shown to be in good agreement with the experiments. However, the VCCT technique predicts higher failure loads. For a brittle adhesive, on the other hand, the predictions by zero/non-zero-thickness cohesive elements and cohesive surfaces reveal notable deviations compared to the experimental results. In contrast to the ductile adhesive, the VCCT technique is revealed to be accurate in predicting the brittle adhesive.

Computational Evaluation of the Fracture Behavior of Porous Asphalt Concrete Exposed to Moisture and Salt Erosion.

Salt erosion has an adverse impact on the durability of asphalt pavements. Porous asphalt concrete is particularly susceptible to the influence of salt. In this study, a finite element model was developed to investigate the fracture behavior of PAC exposed to salt erosion. The 2D heterogeneous structure of PAC was generated with an image-aided approach to computationally study the fracture behavior of PAC. Laboratory SCB tests were conducted to validate the finite element model. The simulation results of the SCB tests indicate that the peak load of the PAC decreased by 21.8% in dry-wet cycles and 26.1% in freeze-thaw cycles compared to the control group. The salt solution accelerated the degradation of the durability of PAC under both dry-wet cycles and freeze-thaw cycle conditions, which is consistent with laboratory tests. After flushing treatment before the drying phase, the peak load of the PAC in salt environments increased by 5.3% compared to that of the samples with no flushing. Salt erosion also results in a higher average value of scalar stiffness degradation (SDEG), and the damaged elements were primarily the cohesive elements in the fracture of the PAC. Additionally, the influence of crucial factors including the void content, adhesion and cohesion, and loading rate on the fracture behavior of the PAC was analyzed. As the void content increases, the average SDEG value of the cohesive elements increases and surpasses the average SDEG value of the adhesive elements at a void content of approximately 9%. The performance of the fine aggregate matrix (FAM) has a much greater impact than the FAM-aggregate interface on the durability of the PAC. And there were more damaged CZM elements with the increase in the loading rate. Salt erosion results in higher SDEG values and a larger number of cohesive damaged elements at each loading rate.

Dataset on simulated microbond tests using finite element method: Simulation cases about different geometrical influences, material behavior, damage evolution, and element meshes.

This data article provides an extensive dataset obtained from finite element (FE) simulations of microbond (MB) tests. The simulations cover a wide range of structural effects and artifacts that influence the results of the MB tests. A total of 432 simulations were performed,taking into account the various factors such as blade geometry and position, plastic behaviour of thermoset and thermoplastic droplets, material properties of the fibres, residual stresses, fracture modes at interfaces, and FE mesh sensitivity analysis. Each FE simulation consists of blade reaction force, blade displacement, fibre displacement, fibre strain and various energy metrics such as interface strain energy, total strain energy, damage energy and plastic dissipation energy. For ease of reference, the individual data files are organised in a systematic naming sequence based on the simulation matrices, detailing the specific abbreviations for each file. A user-friendly interface is also provided to read and visualisethe data from the output files in relation to the simulation matrix. For more information on the interpretation and analysis of this data, please refer to a research article entitled "Mutual dependence of experimental and data analysis features in characterization of fibre-matrix interface via microdroplets (R. Dsouza et al., 2023)".

Multiscale modeling of unfolding and bond dissociation of rubredoxin metalloprotein.

Mechanical properties of proteins that have a crucial effect on their operation. This study used a molecular dynamics simulation package to investigate rubredoxin unfolding on the atomic scale. Different simulation techniques were applied, and due to the dissociation of covalent/hydrogen bonds, this protein demonstrates several intermediate states in force-extension behavior. A conceptual model based on the cohesive finite element method was developed to consider the intermediate damages that occur during unfolding. This model is based on force-displacement curves derived from molecular dynamics results. The proposed conceptual model is designed to accurately identify bond rupture points and determine the associated forces. This is achieved by conducting a thorough comparison between molecular dynamics and cohesive finite element results. The utilization of a viscoelastic cohesive zone model allows for the consideration of loading rate effects. This rate-dependent model can be further developed and integrated into the multiscale modeling of large assemblies of metalloproteins, providing a comprehensive understanding of mechanical behavior while maintaining a reduced computational cost.

Study on the High-Temperature Interaction between Coke and Iron Ores with Different Layer Thicknesses.

Coke plays a key role as the skeleton of the charge column in BF. The gas path formed by the coke layer in the BF has a decisive influence on gas permeability. At high temperatures, the interface between coke and ore undergoes a melting reaction of coke and a reduction reaction of ore. The better the reducibility of the ore, the more conducive it is to the coupling reaction of ore and coke. The melting loss reaction of coke becomes more intense, and the corresponding strength of coke will decrease, which will affect the permeability of the blast furnace and is not conducive to the smooth operation of the blast furnace. Especially with a deterioration in iron ore quality, BF operation faces severe challenges, which makes it necessary to find an effective way to strengthen BF operation. In this study, a melting-dropping furnace was used to develop and clarify the high-temperature interaction between coke and iron ores with different layer thicknesses. The influencing factors were studied by establishing a gas permeability mathematical model and observing the metallographic microscope images of samples after the coke solution loss reaction. The relationships between coke layer thickness, distribution of gas flow, and pressure drop were obtained. The results showed that, under certain conditions, the gas permeability property of a furnace burden has been improved after the coke layer thickness increased. Upon observing the size of coke particles at the interface reaction site, the degree of melting loss reaction can be determined. A smaller particle size indicates more melting loss reaction. A dripping eigenvalue for molten metal was introduced to evaluate the dynamic changes in the comprehensive dripping properties of molten metal of furnace burden, which showed that the dripping eigenvalue for the molten metal could deteriorate because of the unruly thickness and the coke layer thickness should be limited through considering the operational indicators of the blast furnace.

Using Machine Learning and Finite Element Analysis to Extract Traction-Separation Relations at Bonding Wire Interfaces of Insulated Gate Bipolar Transistor Modules.

For insulated gate bipolar transistor (IGBT) modules using wire bonding as the interconnection method, the main failure mechanism is cracking of the bonded interface. Studying the mechanical properties of the bonded interface is crucial for assessing the reliability of IGBT modules. In this paper, first, shear tests are conducted on the bonded interface to test the bonded interface's strength. Then, finite element-cohesive zone modeling (FE-CZM) is established to describe the mechanical behavior of the bonded interface. A novel machine learning (ML) architecture integrating a convolutional neural network (CNN) and a long short-term memory (LSTM) network is used to identify the shape and parameters of the traction separation law (TSL) of the FE-CZM model accurately and efficiently. The CNN-LSTM architecture not only has excellent feature extraction and sequence-data-processing abilities but can also effectively address the long-term dependency problem. A total of 1800 sets of datasets are obtained based on numerical computations, and the CNN-LSTM architecture is trained with load-displacement (F-δ) curves as input parameters and TSL shapes and parameters as output parameters. The results show that the error rate of the model for TSL shape prediction is only 0.186%. The performance metric's mean absolute percentage error (MAPE) is less than 3.5044% for all the predictions of the TSL parameters. Compared with separate CNN and LSTM architectures, the proposed CNN-LSTM-architecture approach exhibits obvious advantages in recognizing TSL shapes and parameters. A combination of the FE-CZM and ML methods in this paper provides a promising and effective solution for identifying the mechanical parameters of the bonded interfaces of IGBT modules.

Modelling of Fatigue Delamination Growth and Prediction of Residual Tensile Strength of Thermoplastic Coupons.

Thermoplastic composites are continuously replacing thermosetting composites in lightweight structures. However, the accomplished work on the fatigue behavior of thermoplastics is quite limited. In the present work, we propose a numerical modeling approach for simulating fatigue delamination growth and predicting the residual tensile strength of quasi-isotropic TC 1225 LM PAEK thermoplastic coupons. The approach was supported and validated by tension and fatigue (non-interrupted and interrupted) tests. Fatigue delamination growth was simulated using a mixed-mode fatigue crack growth model, which was based on the cohesive zone modeling method. Quasi-static tension analyses on pristine and fatigued coupons were performed using a progressive damage model. These analyses were implemented using a set of Hashin-type strain-based failure criteria and a damage mechanics-based material property degradation module. Utilizing the fatigue model, we accurately foretold the expansion of delamination concerning the cycle count across all interfaces. The results agree well with C-scan images taken on fatigued coupons during interruptions of fatigue tests. An unequal and unsymmetric delamination growth was predicted due to the quasi-isotropic layup. Moreover, the combined models capture the decrease in the residual tensile strength of the coupons. During the quasi-static tension analysis of the fatigued coupons, we observed that the primary driving failure mechanisms were the rapid spread of existing delamination and the consequential severe matrix cracking.

Cohesive Properties of Bimaterial Interfaces in Semiconductors: Experimental Study and Numerical Simulation Using an Inverse Cohesive Contact Approach.

Examining crack propagation at the interface of bimaterial components under various conditions is essential for improving the reliability of semiconductor designs. However, the fracture behavior of bimaterial interfaces has been relatively underexplored in the literature, particularly in terms of numerical predictions. Numerical simulations offer vital insights into the evolution of interfacial damage and stress distribution in wafers, showcasing their dependence on material properties. The lack of knowledge about specific interfaces poses a significant obstacle to the development of new products and necessitates active remediation for further progress. The objective of this paper is twofold: firstly, to experimentally investigate the behavior of bimaterial interfaces commonly found in semiconductors under quasi-static loading conditions, and secondly, to determine their respective interfacial cohesive properties using an inverse cohesive zone modeling approach. For this purpose, double cantilever beam specimens were manufactured that allow Mode I static fracture analysis of the interfaces. A compliance-based method was used to obtain the crack size during the tests and the Mode I energy release rate (GIc). Experimental results were utilized to simulate the behavior of different interfaces under specific test conditions in Abaqus. The simulation aimed to extract the interfacial cohesive contact properties of the studied bimaterial interfaces. These properties enable designers to predict the strength of the interfaces, particularly under Mode I loading conditions. To this extent, the cohesive zone modeling (CZM) assisted in defining the behavior of the damage propagation through the bimaterial interfaces. As a result, for the silicon-epoxy molding compound (EMC) interface, the results for maximum strength and GIc are, respectively, 26 MPa and 0.05 N/mm. The second interface tested consisted of polyimide and silicon oxide between the silicon and EMC layers, and the results obtained are 21.5 MPa for the maximum tensile strength and 0.02 N/mm for GIc. This study's findings aid in predicting and mitigating failure modes in the studied chip packaging. The insights offer directions for future research, focusing on enhancing material properties and exploring the impact of manufacturing parameters and temperature conditions on delamination in multilayer semiconductors.

Prediction of Aging Performance of Particle-Filled Polymer Composites Based on RVE Model.

In this study, the aging performance of particle-filled polymer composites (PFPCs) under thermo-oxidative conditions was investigated on multiple scales. High-temperature-accelerated tests were conducted to analyze the effects of aging time and temperature. A representative volume element (RVE) model was established for the PFPCs using a random particle-filling algorithm. A predictive model for the crosslink density was conducted based on the closed-loop chain reaction of polymer oxidation. According to the theory of polymer physics, the relation between the crosslink density and matrix modulus was determined. The particle/matrix interface in the RVE model was represented by the cohesive zone model (CZM). The parameters of the CZM were determined by the inversion techniques. Then, a comprehensive multiscale RVE model was constructed, which was applied to predict the modulus and dewetting strain of the aged PFPCs. The predicted results show good agreement with the test results, which verifies the reliability of our model.

Cohesive Zone Modeling of the Interface Fracture in Full-Thermoplastic Hybrid Composites for Lightweight Application.

With the increasing demand for lightweight and high-performance materials in the automotive and aerospace industries, full-thermoplastic hybrid composites have emerged as a pivotal solution, offering enhanced mechanical properties and design flexibility. This work aims to numerically model the fracture strength in full-thermoplastic hybrid composites made by forming and overmolding organosheets. The mode I fracture was investigated by modeling the behavior of T-joint specimens under a tensile test following the cohesive zone modeling (CZM) approach. The sample was designed to replicate the connection between the laminate and the overmolded part. Double cantilever beam (DCB) specimens were manufactured with organosheets and tested to mode I opening to determine the interlaminar fracture toughness. The fracture toughness out of the mode I test with DCB specimens was used to define the CZM parameters that describe the traction-separation law. Later, due to the particular geometry of the T-join specimens that under tensile load work close to pure mode I, the cohesive parameters were determined by inverse analysis, i.e., calibrating the theoretical models to match experimental results. The fracture resistance T-joint specimens appeared dependent on the fiber-bridging phenomenon during the delamination. In particular, the presence of fiber-bridging visible from the experimental results has been replicated by virtual analyses, and it is observed that it leads to a higher energy value before the interface's complete breakage. Moreover, a correspondence between the mode I fracture toughness of the DCB specimen and T-joint specimens was observed.

Impact Fracture Simulation of Laminated Glass Based on Thick Shell Elements and a Cohesive Zone Model.

Laminated glass is extensively used in automotive windshields, making it crucial to have a comprehensive understanding of its fracture mechanism to ensure driver and pedestrian safety in various windshield impact scenarios. Current research on the cohesive zone model of glass impact failure has encountered challenges related to accuracy and computational efficiency. This paper addresses these issues by utilizing the finite element software LS-DYNA, which integrates a cohesive zone model and thick shell (Tshell) elements to simulate and analyze the impact failure process of laminated glass. The combination of Tshell and cohesive elements was validated using a DCB example. Subsequently, the proposed method was applied to simulate the impact damage on an automobile's front windshield, providing valuable insights from the obtained results. Finally, the influence of curvature, the number of layers, and the thickness ratio of each layer were investigated, leading to some valuable conclusions. Firstly, an increase in the thickness of the upper glass layer correlates with a decrease in the peak acceleration of the dummy-head model due to the ductility of PVB material. Secondly, when a curvature exists, the arched configuration of the windshield promotes higher resistance against impact, consequently leading to increased peak acceleration.

Experimental and Numerical Study of Healing Effect on Delamination Defect in Infusible Thermoplastic Composite Laminates.

The integrity of delaminated composite structures can be restored by introducing a thermally-based healing effect on continuous fiber-reinforced thermoplastic composites (CFRTPC). The phenomenon of thermoplastics retaining their properties after melting and consolidation has been applied by heating the delaminated composite plates above their glass transition temperature under pressure. In the current investigation, the composite is comprised of Methyl methacrylate (MMA)-based infusible lamination resin combined with benzoyl peroxide initiator, which polymerizes into a Polymethyl methacrylate (PMMA) matrix. For the reinforcement, unidirectional 220 gr/m2 glass filament fabric was used. Delamination damage is artificially induced during the fabrication of laminate plates. The distributed delamination region before and after thermally activated healing was determined by using non-destructive testing with active thermography. An experimental approach is employed to characterize the thermal healing effect on mechanical properties. Experimentally determined technological parameters for thermal healing have been successfully applied to repair delamination defects on composite plates. Based on the compression-after-impact (CAI) test methodology, the intact, damaged, and healed composite laminates were loaded cyclically to evaluate the healing effect on stiffness and strength. During the CAI test, the 3D digital image correlation (DIC) technique was used to measure the displacement and deformation fields. Experimental results reveal the difference between the behavior of healed and damaged specimens. Additionally, the numerical models of intact, damaged, and healed composite laminates were developed using the finite element code LS-Dyna. Numerical models with calibrated material properties and tie-break contact constants provide good correlation with experimental results and allow for the prediction of the mechanical behavior of intact, damaged, and healed laminated plates. The comparison analysis based on CAI test results and modal characteristics obtained by the 3D Laser Doppler Vibrometer (Polytec GmbH, Karlsbad, Germany) proved that thermal healing partially restores the mechanical properties of damaged laminate plates. In contrast, active thermography does not necessarily indicate a healing effect.