Bilinear Cohesive Model Research Articles

Determination of fracture parameters for propellant/liner interface affect zone through measuring J-integral

The debonding of propellant/liner interface is one of the primary failure modes that impact the structural integrity of solid rocket motors. In order to study the mechanical properties of propellant/liner interface, tensile tests on rectangular specimens with various prefabricated crack lengths were conducted at different temperatures and loading rates, as well as uniaxial tensile tests on the propellant. Utilizing the digital image correlation(DIC) method, the strain field including the whole rectangular specimen and the propellant part was analyzed, and the change law of the crack length of the rectangular specimen with time was extracted by programming batch processing at different temperatures. Meanwhile, the fracture toughness of the NEPE propellant/liner interface under different loading conditions was obtained using J-integral theory. The results show that the area near the propellant/liner interface forms an interface affect zone (IAZ), where cracks grow steadily in a continuous and smooth way without significant fluctuations. The effect of temperature and loading rate on J-integral is interactive. The lower the temperature, the greater the fracture toughness of the interface and the lower the crack propagation rate. Finally, based on the bi-linear cohesive model, the numerical simulation of the interface affect zone (IAZ) in the rectangular specimen is carried out, which shows that the J-integral can effectively represent the energy required for the interface affect zone failure. A method to determine the fracture parameters for the propellant/liner interface affect zone is proposed.

Mechanical behaviour of bridge expansion joint interfaces filled with novel polyurethane-modified asphalt

The use of Asphalt Plug Joints (APJs) in bridge construction has rapidly expanded due to their distinct advantages in driving comfort, shock absorption, noise reduction, and maintenance ease. However, prolonged exposure to vehicular loads and environmental factors renders this expansion joint susceptible to interface effects and stress concentration, leading to damage, particularly at the interface layers between the expansion joint pavement and the expansion joint-cross joint plate contact areas. To address issues such as stress concentration, this study proposes a polyurethane-modified asphalt concrete mixture as the joint filling material. Pull-out and inclined shear tests are used to identify a bonding agent with superior adhesive properties. Additionally, a finite element model of pull-out specimens containing bilinear cohesive elements is established to simulate the cohesive damage process of the interface layers. Through parameter design analysis, the study examines the influence of individual parameter variations on interface damage behavior, obtaining comprehensive simulated analysis data on interface crack resistance performance. Based on interface test data, a bilinear cohesive force model and a viscoelastic material model of the asphalt concrete mixture are employed to establish a scaled-down model of the improved geometry design of the new APJ. This model simulates the interface mechanical characteristics of APJs under temperature effects and vehicle loads. The results demonstrate that the bilinear cohesive force model effectively simulates the nonlinear behavior of interface bonding slip. Furthermore, the use of SZ interface bonding agents and improvements in cross joint plate form can effectively mitigate issues such as uneven settlement and low-temperature cracking.

Thermo-mechanical coupled three-dimensional finite element simulation analysis of drilling thermoplastic braided carbon fiber composite and optimization of process parameters

Optimizing the process parameters of thermoplastic composites during drilling is paramount for mitigating tearing, delamination, and other forms of damage, while simultaneously enhancing the tensile strength and extending the long-term service life of the overall structure. This study focuses on exploring a thermal-mechanical coupling method for predicting dynamic mechanical progressive failure and determining optimal process parameters during the drilling of Braided Carbon Fiber Reinforced Polyether Ether Ketone (BCF/PEEK). Initially, a thermal conduction constitutive model of BCF/PEEK was developed based on the proposed thermal distribution ratio calculation method. Concurrently, a user-defined material subroutine VUMAT and a bilinear cohesive element model, were implemented on the ABAQUS/Explicit platform to simulate the delamination behaviors of drilling BCF/PEEK using a tapered drill-reamer. Subsequently, a comprehensive BCF/PEEK drilling experiment platform was constructed, and the simulation accuracy of the drilling finite element (FE) model was verified through temperature, thrust force, and hole-wall morphology analyses. Finally, response surface regression models were established for process parameters, and the optimal parameters were predicted and validated. The results indicate that minimum thrust force and temperature can be achieved using the tapered drill-reamer with a spindle speed set at 4878.79/3000 r/min and a feed speed of 34.04/30 mm/min, respectively, with a maximum error of only 8.78 %.

Impact morphology characteristics and damage evolution mechanisms in CFRP laminates for hydrogen storage cylinders

In instances of impact, the preeminent load-bearing framework for on-board hydrogen storage cylinders manifests in the carbon fiber reinforced polymer (CFRP) composite layer. Investigating the morphology of impact and elucidating the mechanism governing damage evolution in CFRP is crucial for optimizing the impact resistance of the cylinder. In this work, impact tests were performed on CFRP laminates with six types of interlaminar mismatch angles, under varying impact energies. The impact damage of laminates was characterized using an extended depth-of-field 3D microscopy. A numerical model, implemented in ABAQUS/Explicit, was devised to scrutinize the mechanical properties and damage mechanisms in laminates. A subroutine (VUMAT) was crafted to effectively predict intralaminar damage, while the application of a bilinear cohesive model served to capture interlaminar damage. In addition, this study delves into the ramifications of impact energy and interlaminar mismatch angles. The results reveal that the increase of impact energy leads to more irreversible damage and energy dissipation within laminates. Compared to damage depth (D) and cross-sectional area (Ac), the assessment of damage volume (V) and surface area (As) are more reflective of energy dissipation in laminates. The laminates with mismatch angle of 0° and “helicoidal” lay-up sequence exhibit the worst impact resistance. Within the range of 24°–90°, laminates with a smaller interlaminar mismatch angle demonstrate better impact resistance. These findings lay the groundwork for optimizing the composite layer to enhance the impact resistance of hydrogen storage cylinders.

Predicting the bond stress–slip behavior of steel reinforcement in concrete under static and dynamic loadings by finite element, deep learning and analytical methods

This paper integrates finite element (FE) and deep learning (DL) methods to predict the bond stress–slip behavior of reinforcing bars in concrete under static and dynamic loading. A bilinear cohesive model is adopted to describe the bond stress–slip constitutive behavior between the steel reinforcement and the surrounding concrete. After calibrating the cohesive model with published experimental results, the influences of different working conditions on the bond strength between reinforcement and concrete are further investigated by performing a sensitivity analysis of critical cohesive parameters. Regarding their influences on the FE predictions, the cohesive parameters are classified as follows: damage initiation strength, fracture energy, and stiffness. It is found that a higher loading rate increases the bond strength without affecting the trend of the bond stress–displacement response. Furthermore, the FE results are extracted as the training samples for the established NARX DL model. After evaluating the static and dynamic bond stress–displacement responses predicted by the proposed NARX DL model, the Pearson correlation coefficient of 0.97 can be achieved with the FE predictions. Finally, an analytical model suitable for engineering applications is proposed, based on the observed trends, to reveal the dominant mechanisms. Work examples are provided for the loading rates of 0.01 mm/s and 5.00 mm/s. The proposed FE, DL and analytical models provide versatile tools to reasonably estimate the bond stress–slip behavior in reinforced concrete under static and dynamic loadings.

Failure analysis of unidirectional FRP with fiber clusters under transverse pure shear

In nearly all the failure criteria for the FRP (fiber reinforced plastic), the fundamental strength of the FRP under uniaxial tensile, compressive and pure shear stress were utilized to construct appropriate mathematical function relationship with the multiaxial stress. Among them, the strength of the FRP under transverse pure shear stress was influenced by various factors and could not be directly obtained from tests. In this paper, transverse pure shear failure behavior of FRP was studied using finite element analysis based on RVE models, the RSE and HC algorithm were utilized to generate the FRP with different fiber distribution, mainly refers the different fiber clusters which was commonly experimentally observed. The bilinear cohesive model and extended linear Drucker Prager model were used to characterize the mechanical behavior of the interface and matrix respectively. Then, constraint of periodic boundary condition was applied for the RVE models. Difference on failure fracture behavior of FRP with different fiber distribution were identified and assessed under transverse pure shear stress, while the connection between stress triaxiality and initial damage equivalent plastic strain of the matrix was also altered. On this basis, the distinction of the maximum transverse pure shear stress, the crack resistance strength and the transverse shear strength was thoroughly discussed and identified, which provided a theoretic reference for the failure analysis of the FRP.

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.

Experimental and Numerical Investigation into the Mechanical Behavior of Composite Solid Propellants Subject to Uniaxial Tension.

To further explore the quasi-static mechanical characteristics of composite solid propellants at low strain rates, an investigation was conducted on the mechanical behavior and damage mechanisms of a four-component hydroxy-terminated polybutadiene (HTPB) propellant by means of experiments and numerical simulation. A uniaxial tensile test and scanning electron microscope (SEM) characterization experiment were carried out. A microstructural model, which accurately represents the mesoscopic structure, was developed via the integration of micro-CT scanning and image-processing techniques. The constructed microstructural model was utilized to conduct a numerical simulation of the mechanical behavior. The experimental results demonstrated that the maximum tensile strength increases with increasing strain rate, and the primary cause of propellant failure at low strain rates is the dewetting phenomenon occurring at the interface between the larger particles and the matrix. The maximum tensile strength is 0.48 MPa when the strain rate is 0.00119 s-1, and the maximum tensile strength is 0.37 MPa when the strain rate is 0.000119 s-1. The simulation results indicated a consistent trend in variation when comparing the simulation and experimental curves. This suggested that the established model exhibits a high level of reliability, and provides a promising approach for carrying out microstructural simulations of heterogeneous propellants in future. The mechanical behavior of the propellant can be effectively described by utilizing a mesoscopic finite element model that incorporates the superelastic constitutive model of the matrix and the bilinear cohesive model. This framework facilitates the representation of mesoscopic damage evolution, which consequently provides insights into the damage mechanism. Additionally, the utilization of such models assists in compensating for the limitations of damage evolution characterization experiments.

Theoretical analysis and numerical simulation of axial tensile behavior of pipe joints based on bi-linear cohesive bond-slip model

PurposeDue to its inexpensive production costs, low stress concentration and maintenance-friendliness, the adhesive bonded pipe joint is frequently utilized for pipe connection. However, further theoretical analysis is needed to understand the debonding failure mechanism of such bonded pipe joints under axial tension.Design/methodology/approachIn this study, based on the bi-linear cohesive zone model, the integrated closed-form solutions were derived by considering the axial stiffness ratio and failure stage to determine the relative interfacial slip, interfacial shear stress and relationship of tension–displacement in the bonded pipe joint.FindingsAdditionally, solutions for the critical bonded length and the ultimate load capacity were put forth. Besides, the numerical study was conducted to verify the theoretical solutions regarding the load–displacement relationship. The interfacial shear stress distribution at different failure stages was presented to understand the interfacial shear stress transmission and debonding process. The effect of bonded length on the ultimate load and ductility of pipe joints was also discussed.Originality/valueThe findings in this study can give a reference for the design of bonded pipe joints in their actual engineering applications.

A hybridized continuum-discontinuum modeling for investigating the fracture and energy evolution characteristics of rock mass under sequential explosive detonation using a bilinear cohesive fracture model

Modelling the dynamic mechanical response of rock mass under sequential explosive detonation, especially the fracture and energy evolution characteristics, is the key to optimizing the spatial distribution and initiation time interval of explosives. A hybridized continuum-discontinuum element method and energy statistics algorithm are implemented in this study to accurately investigate the rock dynamic response induced by sequential detonation loading. Based on the rock fracture status and the stress–strain curve of bilinear cohesive fracture model, an energy statistics algorithm considering all energy components is first proposed, and its accuracy and robustness are verified. Then, the continuum-discontinuum element method and energy statistics algorithm are adopted to simulate the dynamic response of rock mass under sequential detonation loading. The results indicate that the area of cracked interfaces generated by each explosive first increases and then decreases, which caused by explosive 2 is the largest, accounting for 38.19%. Due to the damage to rock caused by the previous detonation loading, the rock around explosive 1 mainly undergoes shear fracture, while the rock around explosives 2 and 3 mainly undergoes tensile fracture. The change trends of various energy components are different, and the explosion energy is mainly converted into friction energy and element kinetic energy.

Fracture and size effect analysis in concrete using 3-D G/XFEM and a CZM-LEFM correlation model: Validation with experiments

The existence of process zones at crack tips in quasi-brittle materials is responsible for the structural size-dependent fracture strength. Hence, accurate prediction tools should not only be capable of predicting the fracture propagation and its load–displacement attributes, but should render themselves capable of the size effect analysis as well. This paper applies a robust 3-D generalized/eXtended finite element method (G/XFEM) algorithm to fracture propagation in concrete and validates its capability to capture the size effect behavior based on the cohesive zone model (CZM). Equipped with the capability of mesh adaptivity, the framework utilized in discretizing the resulting system of equations allows the adoption of independent representation for the background solid mesh and the fracture surfaces without the need for conformity between the two. The fracture simulations on a series of geometrically similar concrete beams under three-point bending have shown very good agreement with available experimental results in terms of both load–displacement (peak load and post-peak) characteristics and size effect behavior. The study also models the experimental problems using Linear Elastic Fracture Mechanics (LEFM) to investigate the structure’s size threshold below which the applicability of LEFM diminishes. The effect of beam geometry and concrete material properties on the LEFM accuracy is assessed by calculating the deviations between LEFM and bilinear cohesive model solutions. The simulation results show that power functions best describe the effect of beam’s depth and tensile strength, and linear functions best fit the effect of elastic modulus and fracture toughness. Subsequently, the study presents and validates an empirical correlation model that modifies the LEFM solution to predict similar ultimate strength as the CZM in addition to capturing the concrete size effect law successfully.

Experimental and numerical studies on the tensile behaviors of thin-ply and thick-ply open-hole laminates

Laminated composites with thin-ply unidirectional prepregs are popular due to their higher in-situ strength as compared with thick-ply laminates due to different failure mechanisms. Thin-ply and thick-ply open-hole laminates were fabricated and the tensile failure characteristics were studied by combining experimental and numerical simulation approaches. Charge-coupled device camera and de-ply technology were used to obtain failure morphologies of open-hole specimens under tensile loading. A finite element model considering progressive damage and bilinear cohesive model was built based on ABAQUS to study the damage mechanism of open-hole laminates. In thin-ply laminates, the damage plane was relatively neat while in thick-ply laminate broom-shaped damage morphologies were observed. Based on de-ply experiments and numerical prediction, it was found that thin-ply laminates had delamination suppression effect and could restrain the propagation of delamination along the loading direction, which explains the distinct damage patterns existed in thin-ply and thick-ply laminates. Based on numerical prediction, it was found that drilling-induced delamination did not affect the tensile strength of open-hole laminates. Propagation of drilling-induced delamination was also restrained in thin-ply laminates. The results of this research suggested that thin-ply prepregs should be considered when delamination poses critical threat to the structural integrity of laminated composite structures.

Finite element analysis of damage mechanisms of composite wind turbine blade by considering fluid/solid interaction. Part I: Full-scale structure

This paper is the first one of two that studies damage mechanisms of composite wind turbine blade using numerical analysis. An important highlight is to consider fully complex failure modes of the skin and web as well as the adhesive failure between them. The main work is divided into three parts: (1) the damage model for composite laminates developed by our research group is applied to the skin and web, and the bilinear cohesive model is used to predict the adhesive failure; (2) finite element model is established by considering fluid/solid interaction to study damage mechanisms of a 59.5 m long blade; (3) the wind velocity model including the normal and extreme wind shear effects as well as the tower shadow effect is considered for fluid elements as boundary conditions. By performing fluid/solid coupling analysis using ABAQUS software, distributions of the stress, displacement and damage for the skin, web and adhesive layer as well as the fluid field information under different wind velocity are explored. In addition, numerical results are also compared with those using the blade element momentum (BEM) theory and other existing work. Some conclusions for predicting damage mechanisms and suggestions for improving damage tolerance and load-bearing ability of the blade are given.

Progressive Damage Behaviour Analysis and Comparison with 2D/3D Hashin Failure Models on Carbon Fibre-Reinforced Aluminium Laminates.

It is known that carbon fibre–reinforced aluminium laminate is the third generation of fibre metal materials. This study investigates the response of carbon fibre–reinforced aluminium laminates (CARALL) under tensile loading and three-point bending tests, which evaluate the damage initiation and propagation mechanism. The 2D Hashin and 3D Hashin VUMAT models are used to analyse and compare each composite layer for finite element modelling. A bilinear cohesive contact model is modelled for the interface failure, and the Johnson cook model describes the aluminium layer. The mechanical response and failure analysis of CARALL were evaluated using load versus deflection curves, and the scanning electron microscope was adopted. The results revealed that the failure modes of CARALL were mainly observed in the aluminium layer fracture, fibre pull-out, fracture, and matrix tensile fracture under tensile and flexural loading conditions. The 2D Hashin and 3D Hashin models were similar in predicting tensile properties, flexural properties, mechanical response before peak load points, and final failure modes. It is highlighted that the 3D Hashin model can accurately reveal the failure mechanism and failure propagation mechanism of CARALL.

Finite element analysis of competitive damage mechanisms of composite scarf adhesive joints by considering thickness effect

Cohesive zone models are demonstrated to fail to predict damage mechanisms of adhesive layer with the increase of thickness, although the load response and ultimate failure strength of adhesive joints are evaluated well to some extent. This paper attempts to explore damage mechanisms of composite scarf adhesive joints with brittle adhesive material under quasi-static loading, where the competitive relationships between the damage behaviors of adhesive layer and the interface debonding between the adhesive layer and composite laminates are investigated. For the epoxy adhesive layer, the damage initiation criteria are developed by considering the hydrostatic pressure, the tension/compression asymmetry and frictional effects, and the damage evolution law is derived using continuum damage mechanics. The bilinear cohesive model and the modified Xu and Needleman’s exponential cohesive model are used to predict the interface debonding between the epoxy adhesive layer and glassy fiber/epoxy composite laminates, respectively. FEA is performed to study two competitive damage mechanisms of composite scarf adhesive joints by discussing the effects of the scarf angle, the thickness of adhesive layer and the shape of cohesive models. Results show that 1. two approaches with/without considering damage mechanisms of adhesive layer can be chosen alternatively within the thickness of adhesive layer 0.5–1.0 mm, by weighing calculation efficiency and precision appropriately; 2. damage features of adhesive layer deserve more attention as the thickness increases beyond 1.0 mm.

Extended DDA with rotation remedies and cohesive crack model for simulation of the dynamic seismic landslide

The failure of earthquake-induced landslides generally involves four stages: tensile cracking, breaking of the substrate, shearing failure, and deposition of the sliding mass, which poses a challenge for high-accuracy simulation of such landslides. To resolve this issue, this paper presents an integrated algorithm by introducing a bilinear cohesive contact model into the proposed discontinuous deformation analysis (DDA) method. In the proposed DDA, the displacement mode of the original DDA is replaced by the multiplicative decomposition of the displacement to improve the distortion and false volume expansion of blocks with large rotations. The timely introduction of the bilinear cohesive contact model can better describe the peak-post softening behavior of rocks consisting of blocks. Finally, the extended algorithm is formulated and implemented in the source program of the DDA and then verified through typical examples. The results indicate that the proposed method can facilitate the generation of tensile cracks and effectively eliminate distortion and false volume expansion of blocks, demonstrating consistent results with the test or field investigation.

Interphase Effect on the Macro Nonlinear Mechanical Behavior of Cement-Based Solidified Sand Mixture.

This paper investigates the interphase effect on the macro nonlinear mechanical behavior of cement-based solidified sand mixture (CBSSM) using a finite element numerical simulation method. CBSSM is a multiphase composite whose main components are soil, cement, sand and water, often found in soft soil foundation reinforcement. The emergence of this composite material can reduce the cost of soft soil foundation reinforcement and weaken silt pollution. Simplifying the CBSSM into a three-phase structure can efficiently excavate the interphase effects, that is, the sand phase with higher strength, the cement-based solidified soil phase (CBSS) with moderate strength, and the interphase with weaker strength. The interphase between aggregate and CBSS in the mixture exhibits the weak properties due to high porosity but gets little attention. In order to clarify the mechanical relationship between interphase and CBSSM, a bilinear Cohesive Model (CM) was selected for the interphase, which can phenomenologically model damage behaviors such as damage nucleation, initiation and propagation. Firstly, carry out the unconfined compression experiments on the CBSSM with different artificial gradations and then gain the nonlinear stress–strain curves. Secondly, take the Monte Carlo method to establish the numerical models of CBSSM with different gradations, which can generate geometric models containing randomly distributed and non-overlapping sand aggregates in Python by code. Then, import the CBSSM geometric models into the finite element platform Abaqus and implement the same boundary conditions as the test. Fit experimental nonlinear stress–strain curves and verify the reliability of numerical models. Finally, analyze the interphase damage effect on the macroscopic mechanical properties of CBSSM by the most reliable numerical model. The results show that there is an obviously interphase effect on CBSSM mechanical behavior, and the interphase with greater strength and stiffness ensures the macro load capacity and service life of the CBSSM; a growth in the interphase number can also adversely affect the durability of CBSSM, which provides a favorable reference for the engineering practice.

Fracture behaviour of human skin in deep needle insertion can be captured using validated cohesive zone finite-element method

Medical needles have shown an appreciable contribution to the development of novel medical devices and surgical technologies. A better understanding of needle-skin interactions can advance the design of medical needles, modern surgical robots, and haptic devices. This study employed finite element (FE) modelling to explore the effect of different mechanical and geometrical parameters on the needle's force-displacement relationship, the required force for the skin puncture, and generated mechanical stress around the cutting zone. To this end, we established a cohesive FE model, and identified its parameters by a three-stage parameter identification algorithm to closely replicate the experimental data of needle insertion into the human skin available in the literature. We showed that a bilinear cohesive model with initial stiffness of 5000 MPa/mm, failure traction of 2 MPa, and separation length of 1.6 mm can lead to a model that can closely replicate experimental results. The FE results indicated that while the coefficient of friction between the needle and skin substantially changes the needle reaction force, the insertion velocity does not have a noticeable effect on the reaction force. Regarding the geometrical parameters, needle cutting angle is the prominent factor in terms of stress fields generated in the skin tissue. However, the needle diameter is more influential on the needle reaction force. We also presented an energy study on the frictional dissipation, damage dissipation, and strain energy throughout the insertion process.

Scale-span modelling of dynamic progressive failure in drilling CFRPs using a tapered drill-reamer

The focus of this study is exploring a scale-span modeling method to simulate the structural mechanical responses and dynamic progressive failure behaviors of carbon fiber reinforced plastics (CFRPs) in drilling. A dynamic progressive failure theory based on modified micromechanics of failure (MMF) criterion with new damage evolution laws was first introduced for intralaminar damage. The relationship between macroscopic stress and microscopic stress was established by the stress amplification factors (SAFs). Meanwhile, an auxiliary deleting element criterion was introduced to delete the serious distortion elements for the sake of preventing nonconvergence of the simulation model, and the bilinear cohesive model with mixed-mode failure criterion was employed to simulate the interlaminar delamination. Then, the drilling behaviors of T700S-12 K/YP-H26 CFRPs using a tapered drill-reamer (TDR) were simulated based on the established scale-span theory model which is implemented by user-defined material subroutine VUMAT built on ABAQUS/Explicit platform. Finally, a series of experiments were performed to validate correctness of the simulation results from the perspectives of hole-making quality evaluation index. Results show that the established scale-span simulation model is shown to agree reasonably well with the experiments, and different kinds of damage behavior of the prefabricated hole can be truly simulated in drilling CFRPs.

Finite Element Analysis of Adhesive Failure of Solid Composite Propellant Multi-material Structure for Underwater Intelligent Equipment

This paper is a consecutive work based on the last paper “Finite element analysis of viscoelastic/damage behaviors of composite solid propellant for underwater intelligent equipment”. This paper studies adhesive failure behaviors of multi-material structure in underwater intelligent equipment, which includes the propellant layer/lining layer/adiabatic layer/aluminum layer. First, the hyperelastic mechanical properties for ethylene–propylene–diene monomer for adiabatic layer are predicted by using the Mooney–Rivlin constitutive model, where the material parameters are fitted by combining with theoretical and experimental load curves. Second, the adhesive failure between the propellant layer and the lining layer is predicted for the designed multiple material layers under shear loads by combining with the bilinear cohesive model, the constitutive models for the propellant and adiabatic layer materials. Finally, the adhesive failure mechanisms and the load responses are studied for the shear specimens.