Journal Information
Vol. 8. Issue 1.
Pages 812-826 (January - March 2019)
Download PDF
More article options
Vol. 8. Issue 1.
Pages 812-826 (January - March 2019)
Original Article
DOI: 10.1016/j.jmrt.2018.06.012
Open Access
A novel failure analysis of SMA reinforced composite plate based on a strain-rate-dependent model: low-high velocity impact
Mengzhou Chang, Zhenqing Wang
Corresponding author

Corresponding author.
, Wenyan Liang, Min Sun
College of Aerospace and Civil Engineering, Harbin Engineering University, Harbin 150001, China
This item has received

Under a Creative Commons license
Article information
Full Text
Download PDF
Figures (18)
Show moreShow less
Tables (2)
Table 1. Material parameters of the three phase system.
Table 2. Comparison between the two models (10J/15J/20J).
Show moreShow less
Additional material (12)

The model of SMA reinforced composite is separated into three parts: reinforce, interphase and matrix. Taking effect of strain rate into consideration, the visco-hyperelastic model is employed in the constitutive law of interphase part. A damage model based on Hashin criterion is developed to simulate the SMA reinforced composite plate subject to low or high velocity impact. Effect of the impact velocity on the shape memory alloy reinforced composite plate under fixed boundary condition is investigated by finite element method, also as the damage state. Low velocity is applied to the model firstly to illustrate the accuracy of parameters and procedures by comparing with the experimental data. In next step, several impact velocities are applied during the simulation process to invest the impact resistance and failure mechanism. Simulation results indicate that damage model of composite is sensitive to loading speed and it is easier for the creak to grow under high velocity.

Low-high velocity impact
Numerical analysis
SMA reinforced composites
Full Text

Fiber reinforced composites have been widely used in aerospace industry due to their unique properties: high stiffness, high strength and low density [1]. However, the low impact resistance in unidirectional cross-ply fiber/matrix laminates is the main limitation in applications. Obviously, the reason for this disadvantage can be traced to the low strength, stiffness between the adjacent layers during manufacturing process which will result in delamination and fiber breakage [2]. The overall property can be improved by adding new materials in matrix, such as short glass fiber and shape memory alloy (SMA).

As for short fibers, the mechanical properties of the composites are strongly influenced by the manufacturing process, such as injection position, sample geometry, pressure and temperature during molding [3,4]. Other factors that affect the mechanical properties, such as fiber location, length, diameter and orientation have been studied by Thomason [5–7]. Also, tensile strength (or elastic modulus) of samples machined perpendicular to mold flow direction are nearly 40% lower than that of samples machined parallel to mold flow direction [8,9].

As for SMA, it has been embedded in the adjacent layer for the reinforcement purpose considering their unique properties: changing shape (elastic modulus) in accordance with temperature and stress [10–12]. SMA has the ability to change shape, to repair damage and to improve impact property according to Sun's experiment and Shariyat's analysis [13,14]. However, the interfacial debonding between SMA and matrix has been found exists extensively in composite. It is still a key problem in the development of SMA reinforced composites [15].

The new structure bounding of the conventional fiber reinforced polymer with metal layer is also found has advantage of subjecting to impact loading [16]. Other structures, such as adding particle [17,18], sandwich plate [19–22] and 3D fabric [23,24], have been developed in recent year.

SMA reinforced composites has been studied in this paper. Several works relevant to low velocity impact of SMA reinforced composites taking the assumption of uniform distribution of martensite phase in SMA has been founded [25,26]. Along with modified linearized Hertz contact law and Brison's model of SMA, Shariyat and his co-authors have proposed several theory to analysis preload effect, impact response and forced vibration of SMA reinforced composites, such as mixed order hyperbolic global-local plate theory [14,18–21]. Strain rate shows significant effect on SMA which may change phase transform process [27]. As for Glass/Epoxy composite, both strengths and modulus at tensile, compressive and in plane shear case are increased generally under higher strain rate [28–30]. However, theoretical models relating low-high velocity have rarely discussed due to the difficulties associated with the face-face contact in simulation (between SMA and matrix; fiber and matrix).

The highlight of our work is to introduce an interphase part with reasonable volume fraction rather than two surfaces. Embedding visco-hyperelastic theory in this part to charactering the effect of strain rate, the rest of the elements can be regarded as single material. Numerical simulations of tensile tests and pull-out tests of SMA reinforced composite has been conducted to verify the accuracy of the model by comparing with experimental data. Furthermore, Finite element analysis has been presented to evaluate the effect of velocity on the impact resistance of SMA reinforced composites.

2Mechanics of material models

The whole system is consisted by three phases: matrix, reinforce and interphase.

2.1Mechanical property of the three phases of composites2.1.1Mechanical property of matrix

For most matrix, taking epoxy resin as an example, a nearly linear elastic behavior can be founded:

where u, ɛ and σ are displacement, strain and stress (i,j = x,y and z are the reference coordinates X), respectively. λm = Emμm/1+μm1−2μm is the Lame's constant, and Gm=Em/21+μm is the shear modulus, Emandμm are the elastic modulus and Poisson's ratio, δij is the Kronecker delta (δij=1 if i=j, and δij=0 otherwise).

2.1.2Mechanical property of reinforce

Differently, mechanical property of reinforce is complicated due to the difference between the materials used. A similar stress–strain relationship as Eq. (2) can be observed if the reinforce is isotropic elastic, taking glass fiber as an example. However, the constitutive model of SMA is sensitive to temperature. Brinson's model [11,12] has been referred to most often in subsequent studies reported in the literature. Here, the constitutive law of SMA based on energy balance equations is derived considering the effects of temperature T and phase conversion between martensite and austensite can be denoted as:

where σin and ɛin are the Piola–Kirchhoff stress and Green strain of SMA; ξ is the martensite fraction charactering the phase conversion; also, Ein, Ω and Θ are the elastic modulus, transformation coefficient and thermal coefficient.

To simplify the application, some assumptions are taken as follows: the elastic modulus is a function of martensite fraction ξ; the transformation function is also connected with elastic modulus:

where Em and Ea are elastic modulus for pure martensite and austenite state, ɛL is the maximum residual strain. The temperature and stress determined equations of phase conversion between martensite and austenite can be expressed as:

(a) Phase conversion to martensite

If T<Ms and σscr<σ<σfcr

(b) Phase conversion to austenite

In Eqs. (5) and (6), parameters am=π/Ms−Mf and aa=π/Af−As. Four important parameters are introduced to character the phase conversion temperature: Ms – start temperature of martensite phase; Mf – finish temperature of martensite phase; As – start temperature of austenite phase; Af – finish temperature of austenite phase. The stress–strain relationship at any temperature T can be obtained according to Eqs. (3)–(6).

2.1.3Mechanical property of interphase

A clear visco-hyperelastic behavior has been observed in composite materials according to our previous work [31]. In a completely analogous manner, we can derive an expression relating stress σin,ijt to strain in a sample that has experienced some continuous strain history given by the function εin,ijt:

where g (or M) is relaxation modulus, t′ is new time variable, ε˙int=∂εint′/∂t′ is strain rate. Using the Boltzmann superposition principle, the stain history εint can be calculated. Also, the relaxation modulus g can be expressed using discrete relaxation spectrum as follows:
where g0, gi and ti are the related parameters which can be obtained from relaxation tests.

2.2The modified Hashin's failure criterion

Hashin's failure criterion accounts for fiber failure and matrix failure is embedded in ABAQUS using subroutine. Two failure models are considered for fiber: tensile failure and compressive failure, the damage variables can be expressed as [32]:

Tensile failure of fiber,

Compressive failure of fiber,

Similarly, tensile failure and compressive failure of matrix can also be obtained [33,34].

Tensile failure of matrix,

Compressive failure of matrix,

Using the damage variables and the related parameters, stress is decreased linearly to zeros once the failure criterion is reached. The stiffness coefficients should be re-calculated at damage state, as:

where smt and smc are the factors that control reduction in shear stiffness according to tensile and compressive failure, respectively. Parameters df=1(1dft)(1dfc) and dm=1−1−dmt1−dmc are the global damage variables charactering fiber and matrix, respectively.

In fact, the relationship between stress and strain during entire damage process can be expressed as:

where H is equivalent function of Hashin's failure criterion, obviously, effected by parameters X, Y and S. As for interphase, the values of X, Y and Z are also influenced by time, this behavior is further studied in Section 3.2.

2.3Homogenization method for the multi-phase system

Considering the two phase system generally used in simulation: reinforce and matrix, cohesive zone model is embedded in the elements between reinforce and matrix [35], as shown in Fig. 1a. A bilinear cohesive law is employed to simulate the interfacial behavior in Section 3. The new material property ‘interphase’ is introduced in the investigation, as shown in Fig. 1b.

Fig. 1.

(a) Two phase system of the composite. (b) Three phase system of the composite.


A clear distinguish are needed to charactering the real size of reinforce in two phase system, taking fiber as an example, Rr is radius of pure material. The process extending two phase system to three phase system is shown in the flow chart (Fig. 2).

Fig. 2.

An algorithm for estimation of interphase parameters.


Two important assumptions are adopted in long fiber reinforced composite: (a) fibers are distributed uniformly in the matrix; (b) mechanical property of fiber and matrix are both homogeneous isotropic. The reinforce volume fraction in two phase system []2 can be expressed as:

where S and n are the area of cross section and the number of fibers, respectively. The comparison is carried out to decide whether the two phase system is appropriate. If a large error is founded between the simulation results and experimental results, then a three phase system are need in finite element analysis (adding interphase). Volume fractions of the three phase system, []3, are as follows:
where kin=rr/Rr, kout=rin/Rr is the geometry parameters charactering the relationship between the radius of reinforce and interphase in three phase system, respectively; θ1=π−4 arccos πRr2/vr2rin2 and θ2=arccos πRr2/vr2rin2 charactering the overlap region of the interphase between two neighboring representative volume elements (RVE).

Considering that the composite is stacked by several layers with different angles along thickness direction. For a single layer, the effective material properties may be related to those of fiber and matrix through using Voigt and Reuss micromechanical rules for the longitudinal direction-1, transverse in-plane direction-2 and out of plane transverse direction-3, respectively [36].

where E¯,G¯ and υ¯ are the equivalent elastic modulus, shear modulus and Poisson's ratio of one layer, respectively.

3Numerical simulation3.1Determination of material parameters

For the sake of simply process extending the two phase system to three phase, the mechanical properties of interphase (taking elastic modulus Ein(t) as an example) are considered changing between Er and Em. The elastic modulus of glass fiber and epoxy resin are in the range of 3–4GPa and 70–80GPa, and value chosen for the study are 3.6GPa and 72.0GPa according to summary in [1,2], respectively. The equivalent elastic modulus of RVE with different configurations is investigated in this section. This can be practically implemented and is used as a basic of determination of material parameters.

The maximum and minimum value of instantaneous elastic modulus E˙rev=dErev/dt=dσrve/dεrve can be obtained according to Eqs. (7) and (8), as follows:

A practical expression can be derived using the simpler form of Eq. (8), gt=g0+g1e−t/t1, as follows:

Parameters g0,g1 and t1 can be determined according to one tensile test curve, weather other parameters are need is depends on comparison of tensile tests under different strain rates [36]. Relationship between elastic modulus and parameter vr2,kout and kin are shown in Fig. 3.

Fig. 3.

Instantaneous elastic modulus of RVE with different parameters koutandkin: (a) vr2=0.1; (b) vr2=0.3; (c) vr2=0.5.


From Fig. 3, it is obviously that the instantaneous elastic modulus of an RVE E˙rve is affected by fractions of the three phases vr3, vin3, vm3 from the original two phase system with fraction vr2. Considering a given fraction vr2, E˙rve is increased with effective elastic modulus of the interphase generally. E˙rve is mainly influenced by outside ratio kout when effective modulus of the interphase is small Ein=3.6GPa. Differently, the instantaneous elastic modulus is mainly influenced by inside ratio kin when effective modulus of the interphase is large Ein=72GPa. If 3.6GPa<Ein<72GPa, E˙rve is increased with increasing of kin or kout. Beyond that, E˙rve is also be founded increased with [vr]2 once parameters kin, kout and Ein are kept as constants. Under different loading speed, the simulation results using three phase model compared against experimental data of tensile tests of glass/epoxy composite (Fig. 4).

Fig. 4.

Tensile tests of glass/epoxy composites under various strain rates and simulation results. (a) stress vs strain; (b) tensile strength vs strain rate [28].

3.2Simulation of tensile test

As mentioned before, the tensile strength can be modified by phenomenological theory, taking tensile strength of interphase as an example:

where XTS is the tensile strength under quasi-static loading, ε˙≈0.

The tensile properties of glass fiber reinforced composites have been studied extensively [1,36]. On the basic of it, effect of loading speed on tensile properties of composites has been experimentally invested by researchers [28]. The results indicate that the elastic modulus, tensile strength and failure strain are increased with strain rate, and an exponential law can be used to model the effect, as shown in Eq. (20).

Simulation results of the three phase model extending fiber fraction [vr]2=0.5 with parameters kin and kout (parameters are shown in Table 1). For different strain rates, from 0.017/s to 85/s, there is a good agreement between simulation results and experimental data. Similarly, simulation results of tensile strength are also acceptable comparing with experimental results.

Table 1.

Material parameters of the three phase system.

  kin[−]  kout[−]  g0[GPa]  g1[GPa]  t1 [s]  a  b 
Glass-epoxy  0.90  1.15  3.6  68.4  0.995  0.24  0.78 
SMA-epoxy  0.96  1.20  3.6  47.4  10.376  0.30  0.92 
3.3Simulation of pull-out test

Pull-out tests of SMA (Nitinol) reinforced matrix(epoxy resin) have been implemented in our lab. The surface of SMA wire was washed with acetone to remove dust particles on surface before embedded them in epoxy. Epoxy resin used was epoxy vinyl resin with an accelerator and hardener in a mass ratio of 1:2%:0.2%, the matrix was uniformly stirred, and then the matrix was removed into a vacuum oven for drying at room temperature for 20min to ensure that the resin bubbles are excluded. The matrix was poured into the mold and cured for 24h in the vacuum oven at room temperature.

Two series of pull-out test were performed to evaluate the interphase property, and the radius of SMA are dSMA=0.5mm and 1.0mm, respectively. As shown in Fig. 5, the length and radius of the three phase model in simulation with cylindrical shape before pull out is 20mm and 10mm, respectively.

Fig. 5.

Illustration of the three phase model used in simulation.


From Fig. 6a, one can conduce that simulation results match well with the experimental results both for dSMA=0.5mm and dSMA=1.0mm. With embedding SMA, dSMA=0.5mm, tensile force is increased quickly with pull displacement (when it is small than 1mm) then kept in steady state with small vibration. The force–displacement relationship indicates typical phase transformation from austensite to martensite with weak interphase effect. Differently, as for embedded SMA with dSMA=1.0mm a typical force–displacement relationship can be observed: force increased quickly to peak value then decreased to a low level suddenly due to the general debonding. Mises stresses of the three phase model shown in Fig. 6b and c, and the comparison reveal that mechanical property of interphase is affected by geometry.

Fig. 6.

Experimental results and simulation results of pull-out tests: (a) force vs displacement; (b) stress state at displacement 1mm, dSMA=0.5mm; (c) stress state at displacement 1mm, dSMA=1.0mm.

4Impact resistance of SMA reinforced composite4.1Numerical model for the simulation

A layer including glass fiber and epoxy resin is regarded as anisotropic material. The stacking sequence of the glass fiber is [0°, 90°]6, and the total thickness of the sample is about 2.4mm (0.2mm for each layer). For each layer without SMA the volume fraction of fiber glass is about 34%, and this is the basic for determining material properties of each layer according to Eq. (17). Between layer 3–4 and layer 9–10, SMA wires are embedded after washing on surface. The distance of two adjacent SMAs is 5mm, the number of SMA at each layer is 15 (total 15×2=30), as shown in Fig. 7a.

Fig. 7.

(a) Illustration of the impact experiment of SMA reinforced composite. (b) Three phase model in the simulation.


In the impact experiment, a sample with dimension Lx×Ly×Lz=100mm×100mm×2.4mm is impacted by a rigid half ball-cylinder on the center of the top surface. Marginal areas on the up and down surfaces of sample are clamped to ensure the fixed boundary condition, as shown in Fig. 7b. The diameter of the half ball of the impactor is 10mm, and the mass of the impactor is 4.2kg (Steel). In the simulation, only the impacted region are considered to simply the finite element model, also, a ball (diameter 10mm) is used to simulate the impactor based on two principles: same mass and same geometry of impact surface, as shown in Fig. 7b.

Three impact energies are considered in the low velocity impact tests 10J, 15J and 20J, and the corresponding impact velocities are 2.18m/s, 2.67m/s and 3.09m/s. As for high velocity, three energies have been considered:100J, 1000J and 10,000J, and the corresponding impact velocities are 6.89m/s, 21.78m/s and 68.9m/s.

The two phase model and three phase model are used in the simulation to fully understand the mechanism of the impact process. However, it is quite unrealistic modeling each glass fiber in ABAQUS considering the complexity both for the two models. Homogenization method is used in each layer both for the two models.

Two phase model: the interface between layer-layer and matrix-SMA are characterized using cohesive zone model (bilinear cohesive law) in ABAQUS.

Three phase model: the interface between layer-layer of three phase model is characterized using subroutine of ABAQUS (VUMAT). As for the interphase between matrix and SMA, a novel part is added to simulate the interphase.

4.2Numerical investigation of SMA reinforced composites under impact.4.2.1Low velocity impact

The simulation results of three aspects (velocity, displacement and energy) of impactor using two models (three phase and two phase) are plotted and compared with the experimental data, as shown in Fig. 8. In Table 2, the relative errors of the two models are shown, from which we can also concludes that the three phase models shows a better result.

Fig. 8.

Comparison between the simulation results and experimental results: (a) noSMA-velocity; (b) noSMA-displacement; (c) noSMA-energy; (d) SMA-velocity; (e) SMA-displacement; (f) SMA-energy.

Table 2.

Comparison between the two models (10J/15J/20J).

  V−noSMA[m/s]  V−SMA[m/s]  D−noSMA[mm]  D−SMA[mm]  E−noSMA[J]  E−SMA[J] 
Maximum error (two)  0.67/0.57/0.96  0.97/0.86/1.49  2.76/1.90/4.21  4.16/2.97/6.69  3.32/4.46/4.14  4.65/6.23/6.19 
Maximum error (three)  0.35/0.15/0.23  0.40/0.16/0.54  1.47/1.51/2.19  1.46/1.37/1.13  2.59/1.16/1.88  2.19/0.86/2.13 

Under lower impact energy, 10J, a generally elastic behavior of the composite laminate can be draw as: deformation is increased with time t at initial state (t<5.6ms) until the maximum deformation is reached, then deformation is decreased (t>5.6ms). This behavior can be founded in the composite with SMA or without SMA, as shown in Figs. 9a, d and 10a, b, e and f. The simulation results of the two models for this case are both acceptable. Whole dynamic impact process of composite without SMA (10J) can be found in videos: Two-phase-noSMA-10J and Three-phase-noSMA-10J; as for SMA reinforced composite, the simulation results using the two models can be found in videos: Two-phase-SMA-10J and Three-phase-SMA-10J.

Fig. 9.

Stress state of simulation results of the composites laminates without SMA using different models.

Fig. 10.

Stress state of simulation results of SMA reinforced composites laminates under low velocity impact using different models.


Differently, the composite is destroyed completely under higher energy, 20J, as shown in Figs. 9c, f and 10d, h. During this simulation, the impactor is founded moving along top layer to the bottom layer of the composite without bounced back. However, the velocity of impactor is reduced due to the friction of the hole of the composite during the breakdown process. This behavior can also be founded in the composite with SMA or without SMA using the same parameters from tensile tests, the simulation process of the two phase model cannot demonstrate this behavior.

As for impact energy 15J and composites without SMA, no clear bounded back behavior has been clearly founded in the test. However, the impactor shows bounded back during this impact process due to the embedding of SMA. This indicates that the resistance property of composite is increased by using SMA, increasing the superelasticity of the composites. Furthermore, the absorbed energy has increased from 14.6J to15J.

4.2.2High velocity impact

In this section, three impact energies have been considered: 100J, 1000J and 10,000J, the related velocities are 6.89m/s, 21.79m/s and 68.9m/s. The simulation times are chosen properly to balance the computational cost and computational accuracy: Ttot=5ms, 1.6ms and 0.5ms for 100J, 1000J and 10,000J, respectively. As shown in Fig. 11a–c, the ratio of velocity to the original velocity is decreased with the increasing of the velocity, both for two phase model and three phase model and both for composite with SMA or without SMA. When comparing the last three figures in Fig. 11, it can be concluded that the absorbed energy is increased with velocity under the range 2.18m/s to 68.9m/s, gradually. As shown in Fig. 12, the absorbed energy of SMA reinforced composite is increased from 10J to 54.56J using the three phase model (19.38J for two phase model); the absorbed energy is smaller for composites without SMA, 22.28J and 30.92J, respectively. Generally, the energy absorbed by SMA reinforced composites is larger than that without SMA, increasing 35–80%.

Fig. 11.

Simulation results of composites under high velocity impact: (a) velocity-100J; (b) velocity-1000J; (c) velocity-10,000J; (d) energy-100J; (e) energy-1000J; (f) energy-10,000J.

Fig. 12.

The absorbed energies of composites under different impact energies.


More details about the simulation process can be founded in Fig. 13. The comparison between the two models indicate that with increasing impact velocity the damage state of composite is more obviously using the three phase model. As shown in Fig. 13h, the damage state is clearly demonstrated in the regions lack of SMA and boundary. Without SMA, the dynamic damage process of composites (100J and 10,000J) using the two models can be found in videos: Two-phase-noSMA-100J, Three-phase-noSMA-100J, Two-phase-noSMA-10,000J and Three-phase-noSMA-10,000J; also, the dynamic damage process of SMA reinforced composites using the two models can be found in videos: Two-phase-SMA-100J, Three-phase-SMA-100J, Two-phase-SMA-10,000J and Three-phase-SMA-10,000J.

Fig. 13.

Stress state of simulation results of SMA reinforced composites laminates under high velocity impact using different models.

4.3Analysis of damage state4.3.1Damage state of different layers after impact

Comparing with the two models, a less integrity can be kept using the three phase model due to the effect of strain rate on the modified Hashin's failure criterion. Without SMA, a hole shape damage can be founded on several layers for both models, also, an arbitrarily damage on the layer 1 can also be founded. With SMA, the layers between the SMA layers (3–4 and 9–10) are damaged more extensive considering the debonding of SMA wires, as shown in Fig. 14.

Fig. 14.

Fracture morphology of the different layers of composite laminate with impact energy 100J: (a) layer 1 (bottom); (b) layer 3; (c) layer 7; (d) layer 10; (e) layer 12 (upper).


The damage area can be separated to two parts: hole area and arbitrarily damage area. The two areas of different layers are shown in Fig. 15. As for composites without SMA, the hole area is generally larger than the arbitrarily damage area for energy level 10J, 15J, 20J, 100J and 1000J; as for highest impact energy 10,000J, the arbitrarily damage area is larger than the hole area, especially for surface layer 1 and 12. When embedding SMAs, the arbitrarily damage area is much larger than the previous case.

Fig. 15.

A statistic of hole area and arbitrarily damage area of the composites.

4.3.2Damage state during impact process

In Fig. 16, the impact processes of composites using different models are demonstrated. A distinct bending behavior along the impact direction during the impact process, when t<2ms, can be observed. After deviating from composites, the velocity of impactor is kept constant. As the same time, the composite is vibrating along the opposite impact direction.

Fig. 16.

Cross-section of the layer 12 during impact, 100J: (a) two phase-noSMA; (b) three phase-noSMA; (c) two phase-SMA; (d) three phase-SMA.


As mentioned before, the displacements of points along Z direction show a parabola shape when the impact energy is 10J. As for higher impact energy 100J and 10,000J, the simulation results indicate a vibrating from negative direction to positive direction with amplitude −9 to 9mm, as shown in Fig. 17. Comparing with composites without SMA, SMA reinforced composites shows larger amplitude after impact, larger 11–17%.

Fig. 17.

The displacements of points on layer 12 during impact process using three phase model.

4.4Working mechanism and failure mechanism of SMA

As mentioned before, SMA sustain load extensively at early time during impact, as shown in Fig. 10a. The working mechanism is mainly focused on the larger elastic modulus and larger strain of SMA. However, the working mechanism is strongly restricted by the interface behavior between SMA and matrix both for two phase model and three phase model.

The characters ‘B’ and ‘U’ denote the SMA in bottom layer (3–4) and upper layer (9–10), respectively. Effective SMA ratio means subtracting the deleted parts of SMA in two phase model. Differently, in three phase model, this denotes the percentage of interphase without deleting. The existence of interphase is the only reason that the SMA working effectively. The effective ratio is decreased with the increasing of impact energy, as shown in Fig. 18. The effective of upper layer is generally lower than that of bottom layer both for the two models. Effective ratio in three phase model is larger comparing with two phase model, that is, the reason why this model demonstrates more absorbed energy in Section 4.2.

Fig. 18.

A statistic of the effective SMA ratio during the impact process: (a) two phase model; (b) three phase model.


The effect of velocity (from low to high) on impact resistance of SMA reinforced composite laminate has been investigated based on a novel strain-rate-dependent model. Three main parts are used to charactering the whole system: matrix, reinforce and interphase. A modified 3D finite element model based on Hashin's failure criterion was employed in ABAQUS to study the destruction process. Some of the simulation results can be concluded as follows:

  • Comparison between the simulation results and our previous work with fixed boundary shows that the parameters and the simulation process of the three phase model are acceptable, meanwhile, a larger error existed in the two phase model. The relative error between the simulation results (energy or force) and the experimental data is controlled smaller than 10% for the three phase model.

  • The absorbed energy is increased with increasing the impact energy (from 10J to 10,000J), and kept at 25–60J for different models.

  • Embedding SMA wires can enhance the resistance property of composites, especially for low velocity.

  • As for SMA reinforced composites under high impact velocity, extensive damage morphology of the composite can be observed both in the center hole areas and the arbitrarily damage areas. This can be explained by the effect of strain rate on failure mechanism of interphase.

  • The damage state of composite after high velocity impact is increasing due to the vibrating behavior.

Conflicts of interest

The authors declare no conflicts of interest.


The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China within the project with Grant Nos. 11472086 and 11532013.

Appendix A
Supplementary data

The following are the supplementary data to this article:


M. Biron.
Thermoplastics and thermoplastic composites.
2nd ed., William Andrew, (2013),
P.K. Mallick.
Fiber-reinforced composites: materials, manufacturing and design.
3rd ed., CRC Press, (1993),
S. Mortazavian, A. Fatemi.
Effects of fiber orientation and anisotropy on tensile strength and elastic modulus of short fiber reinforced polymer composites.
Compos Part B Eng, 72 (2015), pp. 116-129
S. Mortazavian, A. Fatemi.
Fatigue behavior and modeling of short fiber reinforced polymer composites: a literature review.
Int J Fatigue, 70 (2015), pp. 297-321
J.L. Thomason.
The influence of fiber properties of the performance of glass-fibre-reinforced polyamide-6,6.
Compos Sci Technol, 59 (1999), pp. 2315-2328
J.L. Thomason.
The influence of fibre length, diameter and concentration on the modulus of glass fibre reinforced polyamide 6,6.
Compos Part A, 39 (2008), pp. 1732-1738
J.L. Thomason.
The influence of fibre length, diameter and concentration on the impact performance of long glass-fibre reinforced polyamide 6,6.
Compos Part A, 40 (2009), pp. 114-124
Z. Wang, Y. Zhou, P.K. Mallick.
Effects of temperature and strain rate on the tensile behavior of short fiber reinforced polyamide-6.
Polym Compos, 23 (2002), pp. 858-871
Y. Zhou, P.K. Mallick.
A non-linear damage model for the tensile behavior of an injection molded short E-glass fiber reinforced polyamide-6,6.
Mater Sci Eng A, 393 (2005), pp. 303-309
H.K. Cho, J. Rhee.
Nonlinear finite element analysis of shape memory alloy (SMA) wire reinforced hybrid laminate composite shells.
Int J Nonlinear Mech, 47 (2012), pp. 672-678
L.C. Brinson.
One-dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable.
J Intell Mater Syst Struct, 4 (1993), pp. 229-242
L.C. Brinson, R. Lammering.
Finite element analysis of the behavior of shape memory alloys and their applications.
Int J Solids Struct, 30 (1993), pp. 3261-3280
M. Sun, Z. Wang, B. Yang, X. Sun.
Experimental investigation of GF/epoxy laminates with different SMAs positions subjected to low-velocity impact.
Compos Struct, 171 (2017), pp. 170-184
M. Shariyat, S.H. Hosseini.
Accurate eccentric impact analysis of the preloaded SMA composite plates, based on a novel mixed-order hyperbolic global-local theory.
Compos Struct, 124 (2015), pp. 140-151
A. Miramini, M. Kadkhodaei, A. Alipour, M. Mashayekhi.
Analysis of interfacial debonding in shape memory alloy wire-reinforced composites.
Smart Mater Struct, 25 (2015), pp. 015032
J. Bienias, P. Jakubczak, K. Dadej.
Low-velocity impact resistance of aluminium glass laminates – experimental and numerical investigation.
Compos Struct, 152 (2016), pp. 339-348
R. Zhang, Q.Q. Ni, T. Natsuki, M. Iwamoto.
Mechanical properties of composites filled with SMA particles and short fibers.
Compos Struct, 79 (2007), pp. 90-96
N. Chawla, R.S. Sidhu, V.V. Ganesh.
Three-dimensional visualization and microstructure-based modeling of deformation in particle-reinforced composites.
Acta Mater, 54 (2006), pp. 1541-1548
M. Shariyat.
A double-superposition global-local theory for vibration and dynamic buckling analyses of viscoelastic composite/sandwich plates: a complex modulus approach.
Arch Appl Mech, 81 (2011), pp. 1253-1268
M. Shariyat, S.H. Hosseini.
Eccentric impact analysis of pre-stressed composite sandwich plates with viscoelastic cores: a novel global-local theory and a refined contact law.
Compos Struct, 117 (2014), pp. 333-345
M. Shariyat, M. Moradi, S. Samaee.
Enhanced model for nonlinear dynamic analysis of rectangular composite plates with embedded SMA wires, considering the instantaneous local phase changes.
Compos Struct, 109 (2014), pp. 106-118
M. Shariyat, R. Jafari.
A micromechanical approach for semi-analytical low-velocity impact analysis of a bidirectional functionally graded circular plate resting on an elastic foundation.
Meccanica, 48 (2013), pp. 2127-2148
Z. Asaee, F. Taheri.
Experimental and numerical investigation into the influence of stacking sequence on the low-velocity impact response of new 3D FMLs.
Compos Struct, 140 (2016), pp. 136-146
Z. Asaee, S. Shadlou, F. Taheri.
Low-velocity impact response of fiberglass/magnesium FMLs with a new 3D fiberglass fabric.
Compos Struct, 122 (2015), pp. 155-165
M. Meo, F. Marulo, M. Guida, S. Russo.
Shape memory alloy hybrid composites for improved impact properties for aeronautical applications.
Compos Struct, 95 (2013), pp. 756-766
P.L. Vachon, V. Brailovski, P. Terriault.
Impact-induced damage and damage propagation under flexural load in TiNi and Kevlar-stitched carbon/epoxy laminates.
Compos Struct, 100 (2013), pp. 424-435
M. Kadkhodaei, R.K.N.D. Rajapakse, M. Mahzoon, M. Salimi.
Modeling of the cyclic thermo-mechanical response of SMA wires at different strain rates.
Smart Mater Struct, 16 (2007), pp. 2091-2101
M.M. Shokrieh, M.J. Omidi.
Tension behavior of unidirectional glass/epoxy composites under different strain rates.
Compos Struct, 88 (2009), pp. 595-601
M.M. Shokrieh, M.J. Omidi.
Compressive response of glass–fiber reinforced polymeric composites to increasing compressive strain rates.
Compos Struct, 89 (2009), pp. 517-523
M.M. Shokrieh, M.J. Omidi.
Investigation of strain rate effects on in-plane shear properties of glass/epoxy composites.
Compos Struct, 91 (2009), pp. 95-102
M. Chang, Z. Wang, W. Liang.
A visco-hyperelastic model for short fiber reinforced polymer composites: reinforcement and fracture mechanisms.
Z. Hashin.
Failure criteria for unidirectional composites.
J Appl Mech, 47 (1980), pp. 329-334
A. Puck, H. Schürmann.
Failure analysis of FRP laminates by means of physically based phenomenological models.
Compos Sci Technol, 58 (1998), pp. 1045-1067
M. Meo, E. Antonucci, P. Duclaux, M. Giordano.
Finite element simulation of low velocity impact on shape memory alloy composite plates.
Compos Struct, 71 (2005), pp. 337-342
C.S. Jarali, D.R. Mahapatra.
Modeling the interface effect of shape memory alloy composite materials.
Multidiscip Model Mater Struct, 6 (2010), pp. 257-283
K.K. Chawla.
Composite materials: science and engineering.
3rd ed., Springer, (2013),
Journal of Materials Research and Technology

Subscribe to our newsletter

Article options
Supplemental materials
Cookies policy
To improve our services and products, we use cookies (own or third parties authorized) to show advertising related to client preferences through the analyses of navigation customer behavior. Continuing navigation will be considered as acceptance of this use. You can change the settings or obtain more information by clicking here.