Next Article in Journal
Analysis of Physical and Mechanical Properties of the Mortar in the Historic Retaining Wall of the Gediminas Castle Hill (Vilnius, Lithuania)
Next Article in Special Issue
Propagating, Evanescent, and Complex Wavenumber Guided Waves in High-Performance Composites
Previous Article in Journal
Feasibility and Characterization Mortar Blended with High-Amount Basic Oxygen Furnace Slag
Previous Article in Special Issue
Damage Determination in Ceramic Composites Subject to Tensile Fatigue Using Acoustic Emission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Three-Phase Model Characterizing the Low-Velocity Impact Response of SMA-Reinforced Composites under a Vibrating Boundary Condition

College of Aerospace and Civil Engineering, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(1), 7; https://doi.org/10.3390/ma12010007
Submission received: 7 November 2018 / Revised: 14 December 2018 / Accepted: 17 December 2018 / Published: 20 December 2018
(This article belongs to the Special Issue Damage Detection and Characterization of High Performance Composites)

Abstract

:
Structural vibration induced by dynamic load or natural vibration is a non-negligible factor in failure analysis. Based on a vibrating boundary condition, the impact resistance of shape memory alloy (SMA)-reinforced composites was investigated. In this investigation, a modified Hashin’s failure criterion, Brinson’s model, and a visco-hyperelastic model were implemented into a numerical model to characterize the mechanical behavior of glass fiber/epoxy resin laminates, SMAs, and interphase, respectively. First, a fixed boundary condition was maintained in the simulation to verify the accuracy of the material parameters and procedures by a comparison with experimental data. Then, a series of vibrating boundaries with different frequencies and amplitudes was applied during the simulation process to reveal the effect on impact resistances. The results indicate that the impact resistance of the composite under a higher frequency or a larger amplitude is lower than that under a lower frequency or a smaller amplitude.

Graphical Abstract

1. Introduction

Fiber-reinforced composite materials have been widely used and investigated in recent years due to their unique properties, such as high stiffness, high strength, and low density [1]. However, the composites’ applications are limited by their weak impact resistance, especially in unidirectional cross-ply fiber/matrix laminates. Delamination between adjacent layers and debonding between the fiber and the matrix have contributed to the evolution of damage in composites, and this can be explained by the weak interfacial properties [2].
This disadvantage can be overcome by changing the structure of the material; e.g., using a short fiber instead of a long fiber. The mechanical properties of short-fiber-reinforced composites are strongly influenced by manufacturing process factors, such as the injection position, the sample’s geometry, and the pressure and temperature during the molding process [3,4]. Other factors that affect the mechanical properties, such as the fiber location, length, diameter, and orientation, have been studied by Thomason [5,6,7]. The effect of the mould flow direction (MFD) on the interfacial shear strength and the tensile strength of composites has been investigated [8,9]. The results indicate that the tensile strength (and elastic modulus) of samples machined perpendicular to the MFD are nearly 40% lower than that of samples machined parallel to the MFD.
Functional materials and new structures have been used to improve the mechanical properties of composites. Shape memory alloys (SMAs) have been embedded between adjacent layers for a reinforcement purpose considering their shape memory effect; i.e., changing shape (elastic modulus) in accordance with temperature and stress [10,11]. In our previous research, the effect of SMA positions on the damage behavior and impact response (including peak force, displacement, and energy) of laminates subjected to a low-velocity impact have been investigated [12]. In Khalili’s research, the effect of SMA type (wires, plates, strips, tubes, or layers) on dissipation of the impact energies has been studied [13,14]. Shariyat et al. [15,16] have developed a higher-order global–local hyperbolic plate theory aimed at studying the asymmetric displacement fields. Their calculations indicate that SMAs have the ability to change shape, to repair damage, and to improve the impact resistance property of composites [15,16]. Metal layers have been applied in conventional fiber-reinforced polymers to improve the impact resistance [17]. Other structures, such as particle-reinforced [18,19], sandwich plate [20,21], and three-dimensional (3D) fabric [22,23] structures, have been developed in recent years.
We note that a fixed boundary is widely used in experimental analyses. Except for low- and high-velocity impacts, the mechanical properties of a composite under an eccentric impact [24], multiple impacts [25], and compression-after-impact [26] have received much attention due to their role in engineering practices. However, a vibrating boundary condition is unavoidable when membrane structures are subjected to an impact [27] or wind flow [28]. The effect of vibration on the impact resistance of materials has rarely been discussed due to the difficulties associated with real experimental conditions. Only the vibration response during and after an impact can be observed to evaluate the impact resistance and damage state [29,30]. In Pérez’s work, damage to a carbon fiber reinforced polymer (CFRP) induced by a low-velocity impact, and its effect on the vibration response, was investigated by a micro-mechanical approach [31].
Interfacial debonding between the SMA and the matrix, which is one common failure model, is still a key problem in composites [32]. In our previous work, a three-phase model (matrix, reinforcer, and interphase) was introduced to evaluate the mechanical behavior of a fiber-reinforced composite [33,34]. In this paper, this model is further developed to match the SMA-reinforced composites. Based on this model, the effect of vibration on the impact resistance of SMA-reinforced composites is investigated through a series of frequencies and amplitudes.

2. The Three-Phase Model

In this model, SMA-reinforced glass fiber/epoxy composites are regarded as having three phases: a glass fiber/epoxy composite laminate phase, an SMA phase, and an interphase between the SMA and the laminate.

2.1. Material Property of the Glass Fiber/Epoxy Composite Laminate

Composite laminates containing glass fiber and an epoxy matrix are regarded as anisotropic materials at the macroscale. The constitutive model can be expressed as:
{ σ 11 σ 22 σ 33 σ 23 σ 31 σ 12 } = [ c 11 c 12 c 13 0 0 0 c 21 c 22 c 23 0 0 0 c 31 c 32 c 33 0 0 0 0 0 0 c 44 0 0 0 0 0 0 c 55 0 0 0 0 0 0 c 66 ] { ε 11 ε 22 ε 33 ε 23 ε 31 ε 12 }
where ε i j and σ i j (i,j = 1, 2, and 3) are the strain and stress in the ij-direction, respectively; c i j are the coefficients of the stiffness matrix, i,j = 1, 2,…6. Considering an orthotropic material, the directions are defined as: 1: the fiber’s direction; 2: the in-plane transverse direction of the fiber; and 3: the out-of-plane transverse direction of the fiber.
The stiffness matrix at the linear elastic stage can be expressed as:
c 11 = E 1 ( 1 υ 23 υ 32 ) Γ c 22 = E 2 ( 1 υ 13 υ 31 ) Γ c 33 = E 3 ( 1 υ 12 υ 21 ) Γ c 12 = E 1 ( υ 21 + υ 31 υ 23 ) Γ c 23 = E 2 ( υ 32 + υ 12 υ 31 ) Γ c 13 = E 1 ( υ 31 + υ 21 υ 32 ) Γ c 44 = E 12 c 66 = E 13 Γ = 1 / ( 1 υ 12 υ 21 υ 23 υ 32 υ 13 υ 31 2 υ 21 υ 32 υ 13 )
where E i (i = 1, 2, and 3) is the elastic modulus in the i-direction; and E i j and υ i j (i, j = 1, 2, and 3; i j ) are the shear modulus and Poisson’s ratio in the ij-direction, respectively.
The 3D Hashin’s failure criterion, which accounts for fiber failure and matrix failure, was embedded in ABAQUS using the VUMAT subroutine. In this part, two failure models (tensile failure and compressive failure) were considered for fibers with related damage variables, which can be expressed as [35,36]:
Tensile failure of the fiber:
d f t = 1 :   { ( σ 11 X T ) 2 + ( σ 12 S 12 ) 2 + ( σ 13 S 13 ) 2 1 σ 11 > 0
Compressive failure of the fiber:
d f c = 1 :   { | σ 11 X C | 1 σ 11 < 0
where X T and X C are the tensile strength and compressive strength in the longitudinal direction, respectively; S 12 and S 13 are the ultimate shear strength in the 12- and 13-direction, respectively; and d f t and d f c are the damage variables for evaluating the tensile and compressive damage to the fiber, respectively.
Similarly, the tensile failure and compressive failure of the matrix can also be obtained [35,36].
Tensile failure of the matrix:
d m t = 1 :   { ( σ 11 X T ) 2 + ( σ 12 S 12 ) 2 + ( σ 22 Y T Y C ) 2 + σ 22 Y T + σ 22 Y C 1 σ 22 + σ 33 > 0
Compressive failure of the matrix:
d m c = 1 :   { ( σ 11 X T ) 2 + ( σ 12 S 12 ) 2 + ( σ 22 Y T Y C ) 2 + σ 22 Y T + σ 22 Y C 1 σ 22 + σ 33 < 0
where Y T and Y C are the tensile strength and compressive strength in the transverse direction, respectively; and d m t and d m c are the damage variables for evaluating the tensile and compressive damage to the matrix, respectively.
Using the damage variables and the related parameters, the stress is decreased linearly to zero once the failure criterion is reached. The stiffness coefficients obtained from Equation (2) should be recalculated during the damage process, as:
c 11 = E 1 ( 1 υ 23 υ 32 ) Γ ( 1 d f )
c 22 = E 2 ( 1 υ 13 υ 31 ) Γ ( 1 d f ) ( 1 d m )
c 33 = E 3 ( 1 υ 12 υ 21 ) Γ ( 1 d f ) ( 1 d m )
c 12 = E 1 ( υ 21 + υ 31 υ 23 ) Γ ( 1 d f ) ( 1 d m )
c 23 = E 2 ( υ 32 + υ 12 υ 31 ) Γ ( 1 d f ) ( 1 d m )
c 13 = E 1 ( υ 31 + υ 21 υ 32 ) Γ ( 1 d f ) ( 1 d m )
c 44 = E 12 ( 1 d f ) ( 1 s m t d m t ) ( 1 s m c d m c )
c 55 = E 23 ( 1 d f ) ( 1 s m t d m t ) ( 1 s m c d m c )
c 66 = E 13 ( 1 d f ) ( 1 s m t d m t ) ( 1 s m c d m c )
where s m t and s m c are the factors that control the reduction in shear stiffness according to the tensile and compressive failure, respectively. The parameters d f = 1 ( 1 d f t ) ( 1 d f c ) and d m = 1 ( 1 d m t ) ( 1 d m c ) are the global damage variables characterizing the fiber and the matrix, respectively.

2.2. Material Property of the SMA

A similar stress-strain relationship to Equation (1) can be observed if the reinforced material is isotropic elastic, i.e., glass fiber or carbon fiber. However, the constitutive model of SMA is sensitive to temperature and stress. Among all of the models, Brinson’s model [37,38] has been used widely due to its high accuracy. Here, the constitutive law of SMA based on energy balance equations is derived considering the effect of the temperature T, and the phase conversion between martensite and austensite can be expressed as:
d σ s = E s ( ε s , ξ , T ) d ε s + Ω ( ε s , ξ , T ) d ξ + Θ ( ε s , ξ , T ) d T
where σ s and ε s are the Piola–Kirchhoff stress and the Green strain of SMA; ξ is the martensite fraction characterizing the phase conversion; and E s , Ω , and Θ are the elastic modulus, the transformation coefficient, and the thermal coefficient, respectively.
Some assumptions are made to simplify the application, as follows: the elastic modulus E s is a function of martensite fraction ξ ; and the transformation coefficient is also connected with the elastic modulus:
E s ( ε , ξ , T ) = E s ( ξ ) = E a + ξ ( E m E a )
Ω ( ξ ) = ε L E ( ξ )
where E m and E a are the elastic moduli for martensite and austenite; and ε L is the maximum residual strain. The equations of the phase conversion between martensite and austenite can be expressed as (determined by temperature and stress):
(a) Phase conversion to martensite:
If   T > M s   and   σ s c r + C m ( T M s ) < σ < σ f c r + C m ( T M s ) ξ S = 1 ξ S 0 2 cos { π σ s c r σ f c r ( σ σ f c r C m ( T M s ) ) } + 1 + ξ S 0 2
ξ T = ξ T 0 ξ T 0 1 ξ S 0 ( ξ S ξ S 0 )
If   T < M s   and   σ s c r < σ < σ f c r ξ S = 1 ξ S 0 2 cos { π σ s c r σ f c r ( σ σ f c r ) } + 1 + ξ S 0 2
ξ T = ξ T 0 ξ T 0 1 ξ S 0 ( ξ S ξ S 0 ) + Δ T ξ
Δ T ξ = { 1 ξ T 0 2 [ cos ( a m ( T M f ) ) + 1 ] : M f < T < M s   and   T < T 0 0 : else
(b) Phase conversion to austenite:
If   T > A s   and   C a ( T A f ) < σ < C a ( T A s ) ξ = ξ 0 2 { cos [ a a ( T A s σ C a ) ] + 1 }
ξ S = ξ S 0 ξ S 0 ξ 0 ( ξ 0 ξ )
ξ T = ξ T 0 ξ T 0 ξ 0 ( ξ 0 ξ )
where M s and M f denote the start and finish temperature of the martensite phase, respectively; A s and A f denote the start and finish temperature of the austenite phase, respectively; C m and C a are the material properties describing the relationship between stress and temperature of the martensite phase and the austenite phase, respectively; and σ s c r and σ f c r represent the critical transformation stress at the start and the finish of the transformation, respectively. In Equations (10) and (11), parameters a m = π / ( M s M f ) and a a = π / ( A f A s ) . The stress-strain relationship at an arbitrary temperature T can be obtained according to Equations (8)–(11).

2.3. Material Property of the Interphase

In our previous research, the time (or strain rate)-dependent stress-strain relationship of a glass-fiber-reinforced polymer composite has been investigated [33,34]. In this model, the stress σ i n , i j ( t ) of the interphase is a function of strain history, and can be expressed as:
σ i n ( t ) = 0 t g ( t t ) ε ˙ i n ( t ) d t
where g is the relaxation modulus, t is the new time variable, ε ˙ i n ( t ) = ε i n ( t ) / t is the strain rate. The strain history ε i n ( t ) can be obtained using Boltzmann’s superposition principle. Also, the relaxation modulus g can be expressed using the discrete relaxation spectrum, as follows:
g ( t ) = g 0 + i = 1 n g i e t t i
where g 0 , g i , and t i are the related parameters that can be obtained from relaxation tests. The parameters in the interphase part can be obtained by fitting tensile and pull–put tests with different loading speeds [33,34].

2.4. Boundary Condition

For a sample with dimensions L x × L y × L z = 100 mm × 100 mm × (n × 0.2 mm), n = 16 for our experiment, a fixed boundary is employed to investigate the accuracy of the numerical model. The frequency and the amplitude of the four sides x = 0 , x = L x , y = 0 , and y = L x are zero, and are shown in Figure 1a.
As for the vibrating boundary condition (Figure 1b), the movement of the boundaries can be expressed as:
z = A sin ( 2 π f t )
where A and f are the amplitude and the frequency of the vibration, respectively.

3. The Effect of the Fixed Boundary Condition on Impact Resistance

The fixed boundary condition is kept in the experimental process.

3.1. Composite Laminates

Simulation results of the model under a fixed boundary were investigated and compared with the experimental results. The stacking sequence of the laminate is [0°,90°]8; 0° and 90° are the glass fiber’s layer angles in the X-direction. The sample ( L x × L y × L z = 100 mm × 100 mm × 3.2 mm) was subjected to an impact from a rigid half-ball cylinder at the center of the top surface, as shown in Figure 1. The half ball’s diameter is 14 mm, and its mass is 8 kg (Steel). Two impact energies were considered in the tests (32 J and 64 J), and the corresponding impact velocities are 2.83 m/s and 4 m/s, respectively.
The material parameters were obtained by fitting with the experiments [12,34], as shown in Table 1.
SMA wires (diameter, 0.2 mm) were embedded in the middle layer laminate (between layer 8 and layer 9) at a distance of 5 mm (with a total of 21 wires for a model). In summary, four types of experiments have been conducted, as shown in Table 2.

3.2. Simulation Result: Damage State During the Impact Process

The model has been built in ABAQUS according to the details mentioned above, and the meshed result is shown in Figure 2. In ABAQUS, the ‘Rigid body’ element type of impactor was chosen; and the ‘C3D8R’ element type of laminate and interphase was chosen. The total number of elements in this model was 102,996. In the center region, the mesh net was optimized to ensure the accuracy of the model. The ‘Explicit’ step (t = 0.01 s), integrated into VUMAT (the constitutive relations of laminate, SMA, and interphase), was chosen to simulate the impact process.
Group A1: The simulation results on the composite laminate under an impact at different times (0.0015 s, 0.0030 s, 0.0045 s, 0.0060 s, 0.0075 s, and 0.0090 s) are shown in Figure 3 and Figure 4. Under the lower impact energy (32 J), a generally elastic behavior of the composite laminate can be expressed as: the deformation is increased with time t from the initial state (t = 0) to the maximum deformation (t ≈ 0.0045 s), then the deformation is decreased.
As shown in Figure 3, it is clear that the impactor has been bounced back by the composite laminate. Three layers—layer 16, layer 8, and layer 1—have been chosen to investigate the fracture morphology during the impact, as shown in Figure 4. A hole-shaped damage region can be found on several layers, especially layer 16; however, the damage region on layer 8 and layer 1 is smaller, as shown in Figure 4.
Group A2: Different from group A1, the composite is destroyed completely under the higher impact energy (64 J), as shown in Figure 5. During the simulation process, the impactor was found to move along the top layer to the bottom layer of the composite without being bounced back (t > 0.0030 s). This can also be confirmed by the experimental result, as shown in Figure 5g. The velocity of the impactor was reduced in the breakdown process and then remained as a constant.
From Figure 6c, a larger damage region can be found at the final state. This is different from group A1, and can be explained by delamination due to the friction between the impactor and layer 1 during penetration.
Groups A3 and A4: Embedding SMA alloys is an effective way to improve the impact resistance of composite laminates. As shown in Figure 7a, the SMA was stretched to a larger strain in the case of 32 J. In Figure 7b, a broken or invalid state of the SMA is obtained due to the larger strain, which is beyond the critical value. More specifically, five SMAs in the center region were chosen to demonstrate the working mechanism, as shown in Figure 8c,d.
Beyond that, the defect is obvious: the damage region of layer 8 (contact with the SMA) is larger than that of group A1 and A2, as shown in Figure 8a,b.

3.3. Simulation Result: Absorbed Energy and Contact Force

Two important items—absorbed energy and contact force—were obtained from ABAQUS and compared with the experimental data, as shown in Figure 9. The relative errors of energy comparison are 5.7%, 6.3%, 13.8%, and 7.3% for groups A1–A4, respectively. As for the force-time curves, more inflection points can be found on the simulation curves due to the breaking and removal of the element. Even so, the tendency of the simulation curves agrees well with the experimental curves. The relative errors of the force comparison are 7.3%, 8.1%, 8.9%, and 9.2% for groups A1–A4, which can also be accepted.
From the comparison in this section, the accuracy of the material parameters and the accuracy of the model used in ABAQUS can be tested. Furthermore, the constitutive relationship of SMA and the processing of the interface (between the SMA and the laminate) are also regarded as appropriate. This is the foundation of the numerical simulation in Section 4.

4. The Effect of the Vibrating Boundary Condition on Impact Resistance

The simulation of the model under the vibrating boundary condition under impact is investigated in this section. The same model ( L x × L y × L z = 100 mm × 100 mm × 3.2 mm, 16 layers) as shown Section 3.1 was subjected to an impact by a rigid half-ball cylinder with a fixed energy of 32 J at the center of the top layer. Two types of composites have been investigated: without and with SMA.

4.1. The Effect of Amplitude

In order to fully understand the influence of amplitude on the impact resistance, a low frequency f = 1000 cycles/s was maintained (10 cycles during the simulation process, ttot = 0.01 s). Several amplitudes A were chosen for the study, as shown in Table 3.
The positive value in Table 3 means that the movement direction of the boundary is contrary to the impactor’s movement direction (+z direction) at the initial state. The negative value means that the movement directions of the impactor and the boundary are the same at the initial time.
When applying an amplitude A to control the movement of the boundary, the morphologies of the composite laminate after impact are shown in Figure 10 and Figure 11. As shown in Figure 10a, the damage region is close to the size of the impactor. From Figure 10a–h, the damage regions are increased as the value of the applied amplitude increases; however, general damage with a larger area can be observed in case of A = 0.016 m, as shown in Figure 10g. It is interesting that the damage state depends on the absolute value of the amplitude rather than the value when comparing related groups, e.g., Figure 10c,d. A similar conclusion can be obtained for group C from the simulation results shown in Figure 11.
The details of the impact process for the half model in group B are shown in Figure 12. From Figure 12a, a representative impact process is shown: gradual damage with time t, similar to group A1. With increasing amplitude, more elements in the center region have been removed due to the large deformation (or severe vibration), as shown in Figure 12b–d. As for the larger amplitude, 0.016 m, a larger damage region is observed, and even separation from the main model, t = 0.001 s, as shown in Figure 12d. More importantly, the time at which the clear damage region can be observed has changed from t = 0.005 s to 0.0002 s due to the relative movement. The effect of the value of an amplitude can be further validated by comparing Figure 12b,e.
Three amplitudes are shown to demonstrate the impact process of group C: A = 0.0032 m, 0.008 m, and 0.016 m, as shown in Figure 13. From Figure 13a, layer damage along the SMA direction can be found during the impact process. This is due to the weak impact resistance of the laminate after being separated from the SMA, especially for layer 8. With increasing amplitude, a clear delamination can be observed between layer 8 and layer 9, t = 0.0005 s, as shown in Figure 13b. As for group C4, clear damage is shown at an early time, t = 0.0004 s. The damage state continues to extend even after separation due to the vibrating of the boundary, as shown in Figure 13c.
With different amplitudes, the absorbed energy and contact force are plotted against time, as shown in Figure 14. For group B1 and C1, the maximum value of absorbed energy is 32 J; however, the values at t = 0.01 s are 30.04 J and 28.6 J, respectively. As the value of the amplitude increases, the absorbed energy is decreased. For group B3, B5, and B7, the maximum value of absorbed energy is 9.24 J, 4.24 J, and 0.30 J, respectively. The maximum value of absorbed energy in the related opposite-direction groups, i.e., B3–B4, B5–B6, and B7–B8, is maintained at same level and about 20% lower. For the related group C, the absorbed energy is 3.5 times larger. For the high amplitude, the absorbed energy is close to zero, as in group B7 and B8. From t = 0 to 0.004 s, the effect of vibrations on the energy-time curve can be observed, and the energy is kept constant, as shown in Figure 14a. In Figure 14b, the force is also affected by the vibrations, and shows more dramatic changes when compared with Section 3.

4.2. The Effect of Frequency

In order to fully understand the influence of frequency on the impact resistance, a small amplitude A = 0.0032 m was maintained considering the small influence of this level. Several frequencies f were chosen for the study, as shown in Table 3.
During the simulation, f = 100 cycles/s to 500 cycles/s in groups D1 to D3, the damage states are close to each other at t = 0.01 s, and only a small hole can be found according to the fracture morphology shown in Figure 15a and Figure 16b,c. As for the higher frequency, f = 2000 cycles/s to 10,000 cycles/s, the damage states show randomicity with a larger area, as shown in Figure 15d,e.
Damage states of SMA-reinforced composite laminates are shown in Figure 15f–j. Overall, the effect of frequency on the damage state is similar to that in group D. It should be noted that the separated region in group E5 retains a more complete shape, as shown in Figure 15j.
Applying different frequencies f to the boundary, the simulation process of two groups, D5 and E5, were investigated to demonstrate the details of the damage morphologies of the composite laminate, as shown in Figure 16. For group D4, a hole-shaped damage region gradually appears at time t = 0.0025 s, which increases with the impact process, and, at time t = 0.01 s, the center region is damaged completely. For group E5, delamination is observed, except for in the hole-shaped damage region. This is mainly due to the SMA’s global enhancement effect.
The relationship between absorbed energy and time or contact force and time can be found in Figure 17. As shown in Figure 17a, the absorbed energy is decreased as the frequency increases. Considering the very low frequency f = 100/s, the maximum value of absorbed energy for group D1 and E1 is same: 32 J. Considering the high frequency f = 10,000/s, the absorbed energy for group E5 is 2.36 J, and for D5 the absorbed energy is 1.79 J, which means that the composites can barely bear the impact. As shown in Figure 17b, the maximum value of force in the case of f < 2000 cycle/s is in the range of 7–7.5 N for group D and 9–12 N for group E. For group D5 and group E5, a saltation is observed when comparing the maximum value of force between adjacent groups, which is mainly due to the transient change in velocity.
More important, the maximum value of absorbed energy and contact force for SMA-reinforced composites is generally 15–30% larger than that of pure glass-fiber-reinforced composites under the same amplitude or frequency.

4.3. Statistical Analysis of the Damage State

In Table 4, the maximum energy, Emax, the energy at time 0.01 s, Et = 0.01, the maximum force, Fmax, and the average force, Favg, of different groups are shown. An Emax = 32 J denotes a rebound behavior of the impactor. The average force is defined as:
F a v g = i = 1 F i N
where N is the total number of output data of force Fi within a time of 0.01 s. More importantly, the average value is still in accordance with A and f.
In Figure 18, the ratio between the damaged area and the whole model is plotted against time. Generally, the damaged areas of laminates for A < 0.0032 m or f < 500 cycles/s are kept at <5%. In Figure 18a, the damaged areas of the top layer (layer 16) for A > 0.0032 m or f > 500 cycles/s are kept at 30–50%. Differently, the damaged areas of the middle layer (layer 8) and the bottom layer (layer 1) under the same conditions are kept at 30–70% and 30–65%, respectively.
In Figure 19a, the average damage area was investigated by calculating 16 layers, and the results indicate a 4% average damage area for a small amplitude and frequency and nearly a 50% average damage area for a large amplitude and frequency at time t = 0.01 s. Figure 19b shows the damage state of the SMA. Damage to the SMA can be observed at an earlier time than damage to the composite laminates, which means that SMA has advantages for absorbing energy.

4.4. Mathematical Expression: Effect of Amplitude and Frequency

In this section, the relationship between velocity and amplitude (or frequency) is investigated. The results show a clearly inverse proportion, as follows:
Δ V = k a i = 0 m i A i ( for   given   f )
Δ V = k f i = 0 n i f i ( for   given   A )
where ka and mi are parameters related to amplitude; and kf and ni are parameters related to frequency. Inserting Equation (16) into the energy equation, the relationship between absorbed energy and amplitude (or frequency) is obtained as follows:
E = 1 2 m V 0 2 1 2 m ( V 0 k i = 0 m i A i i = 0 n i f i ) 2
where E and m are the absorbed energy and the mass of the impactor, respectively. Considering that F a v g = m × Δ v / Δ t , the average force is also obtained as:
F a v g = k t t o t i = 0 m i A i i = 0 n i f i
The simulation results using the abovementioned equations are shown in Figure 20 (‘+’ and ‘−’ denote the values of amplitudes; ‘Sim’ denotes the simulations in Section 4; and ‘Eq’ denotes ‘equation’). In this simulation, i = 3 was chosen to simplify the expression. The comparison indicates that both the energy and force can be predicted using Equations (17) and (18).

5. Conclusions

A three-phase model has been proposed to investigate the effect of a vibrating boundary on the impact resistance of shape memory alloy (SMA)-reinforced composite laminates. Some conclusions have been obtained based on the simulation results under different amplitudes and frequencies.
(1)
Under a smaller amplitude (A < 0.0032 m) and a lower frequency (f < 500 cycles/s), the absorbed energy and contact force of composite laminates are similar to that under a fixed boundary condition. In contrast, both a high frequency and a high amplitude can weaken the impact resistance of composite laminates, where extensive damage can be observed rather than a hole-shaped damage region.
(2)
The absolute value of amplitude has a greater influence on the impact resistance than the movement direction of the laminates at the initial time. The absorbed energy and contact force in the positive direction are about 20% larger than that in the negative direction.
(3)
Embedding an SMA can improve the impact resistance of composite laminates due to the superelasticity. In this study, embedding an SMA was found to increase the absorbed energy and contact force by about 15–30%. Also, embedding an SMA can change the damage morphology with respect to shape and proportion.
The proposed three-phase model can provide a tool for efficiently investigating the damage behavior of SMA-reinforced composite laminates.

Author Contributions

Conceptualization, J.H.; methodology, M.C.; software, M.C.; validation, M.S.; formal analysis, M.C.; investigation, M.C.; resources, M.S.; data curation, M.S.; writing (original draft preparation), F.K.; writing (review and editing), F.K.; visualization, F.K.; supervision, J.H.; project administration, J.H.; funding acquisition, J.H.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 11872157, 11472086, and 11532013.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

A amplitude
A s , A f start and finish temperature of austenite phase
a m , a a material parameters of SMA
C m , C a thermal coefficient of martensite phase and austenite phase
c i j stiffness matrix
d m t , d m c , d f t , d f c damage variables
E , E max , E t = 0.01 s absorbed energy, maximum value of energy, and absorbed energy at 0.01 s
E i , E i j elastic modulus of laminate in the i-direction, shear modulus in the ij-direction
E s , E a , E m elastic modulus of SMA, elastic modulus for martensite phase and austenite phase
F i , F max , F avg contact force, maximum value of force, average value of force
f frequency
g 0 , g i parameters of the relaxation modulus
k a , k f parameter related to amplitude, parameter related to frequency
L x , L y , L z sample dimension in the x, y, and z direction
M s , M f start and finish temperature of martensite phase
m , m i mass, parameters related to amplitude
n i parameters related to frequency
S 23 , S 13 , S 12 ultimate shear strength in the 23, 13, and 12 direction
s m t , s m c reduction factors
T temperature
t , t , t i , t t o t time, new time variable, time parameters of the relaxation modulus, total time
V 0 , Δ V initial velocity, velocity increment
X T , X C tensile strength and compressive strength in the longitudinal direction
Y T , Y C tensile strength and compressive strength in the transverse direction
ε i j strain of laminate in the ij-direction
ε s , ε L strain of SMA, maximum residual strain of SMA
ε i n , ε ˙ i n strain of interphase, strain rate of interphase
Θ thermal coefficient
σ i j stress in the ij-direction
σ s , σ i n stress of SMA, stress of interphase
υ i j Poisson’s ratio of laminate in the ij-direction
Ω transformation coefficient
ξ , ξ S , ξ T martensite fraction, martensite fraction of stress, and temperature effect
ξ 0 , ξ S 0 , ξ T 0 martensite fraction at the initial state, initial stress, and initial temperature

References

  1. Biron, M. Thermoplastics and Thermoplastic Composites, 2nd ed.; William Andrew: Waltham, MA, USA, 2013. [Google Scholar]
  2. Li, M.; Gu, Y.; Liu, Y. Interfacial improvement of carbon fiber/epoxy composites using a simple process for depositing commercially functionalized carbon nanotubes on the fibers. Carbon 2013, 52, 109–121. [Google Scholar] [CrossRef]
  3. Mortazavian, S.; Fatemi, A. Effects of fiber orientation and anisotropy on tensile strength and elastic modulus of short fiber reinforced polymer composites. Compos. Part. B Eng. 2015, 72, 116–129. [Google Scholar] [CrossRef]
  4. Mortazavian, S.; Fatemi, A. Fatigue behavior and modeling of short fiber reinforced polymer composites: A literature review. Int. J. Fatigue 2015, 70, 297–321. [Google Scholar] [CrossRef]
  5. Thomason, J.L. The influence of fibre length, diameter and concentration on the strength and strain to failure of glass fibre-reinforced polyamide 6,6. Compos. Part. A 2008, 39, 1618–1624. [Google Scholar] [CrossRef] [Green Version]
  6. Thomason, J.L. The influence of fibre length, diameter and concentration on the modulus of glass fibre reinforced polyamide 6,6. Compos. Part. A 2008, 39, 1732–1738. [Google Scholar] [CrossRef] [Green Version]
  7. Thomason, J.L. The influence of fibre length, diameter and concentration on the impact performance of long glass-fibre reinforced polyamide 6,6. Compos. Part. A 2009, 40, 114–124. [Google Scholar] [CrossRef] [Green Version]
  8. Uematsu, H.; Suzuki, Y.; Iemoto, Y. Effect of maleic anhydride-grafted polypropylene on the flow orientation of short glass fiber in molten polypropylene and on tensile properties of composites. Adv. Polym. Tech. 2017, 37, 1755–1763. [Google Scholar] [CrossRef]
  9. Zhou, Y.; Mallick, P.K. Fatigue performance of injection-molded short e-glass fiber reinforced polyamide- 6,6. II. Effects of melt temperature and hold pressure. Polym. Composite. 2011, 32, 268–276. [Google Scholar] [CrossRef]
  10. Cui, B.; Yao, J.; Wu, Y. Effect of cold rolling ratio on the microstructure and recovery properties of Ti-Ni-Nb-Co shape memory alloys. J. Alloy. Compd. 2019, 772, 728–734. [Google Scholar] [CrossRef]
  11. Jani, J.M.; Leary, M.; Subic, A. A review of shape memory alloy research, applications and opportunities. Mater. Des. 2014, 56, 1078–1113. [Google Scholar] [CrossRef]
  12. Sun, M.; Wang, Z.; Yang, B. Experimental investigation of GF/epoxy laminates with different SMAs positions subjected to low-velocity impact. Compos. Struct. 2017, 171, 170–184. [Google Scholar] [CrossRef]
  13. Khalili, S.M.R.; Botshekanan Dehkordi, M.; Carrera, E.; Shariyat, M. Non-linear dynamic analysis of a sandwich beam with pseudoelastic SMA hybrid composite faces based on higher order finite element theory. Compos. Struct. 2013, 96, 243–255. [Google Scholar] [CrossRef]
  14. Khalili, S.M.R.; Botshekanan Dehkordi, M.; Shariyat, M. Modeling and transient dynamic analysis of pseudoelastic SMA hybrid composite beam. Appl. Math. Comput. 2013, 219, 9762–9782. [Google Scholar] [CrossRef]
  15. Shariyat, M.; Hosseini, S.H. Accurate eccentric impact analysis of the preloaded SMA composite plates, based on a novel mixed-order hyperbolic global-local theory. Compos. Struct. 2015, 124, 140–151. [Google Scholar] [CrossRef]
  16. Shariyat, M.; Moradi, M.; Samaee, S. Enhanced model for nonlinear dynamic analysis of rectangular composite plates with embedded SMA wires, considering the instantaneous local phase changes. Compos. Struct. 2014, 109, 106–118. [Google Scholar] [CrossRef]
  17. Bienias, J.; Jakubczak, P.; Dadej, K. Low-velocity impact resistance of aluminium glass laminates—Experimental and numerical investigation. Compos. Struct. 2016, 152, 339–348. [Google Scholar] [CrossRef]
  18. Zhang, R.; Ni, Q.Q.; Natsuki, T. Mechanical properties of composites filled with SMA particles and short fibers. Compos. Struct. 2007, 79, 90–96. [Google Scholar] [CrossRef] [Green Version]
  19. Chawla, N.; Sidhu, R.S.; Ganesh, V.V. Three-dimensional visualization and microstructure-based modeling of deformation in particle-reinforced composites. Acta. Mater. 2006, 54, 1541–1548. [Google Scholar] [CrossRef]
  20. Susainathan, J.; Eyma, F.; De Luycker, E. Experimental investigation of impact behavior of wood-based sandwich structure. Compos. Part. A 2018, 109, 10–19. [Google Scholar] [CrossRef]
  21. Shariyat, M. A double-superposition global-local theory for vibration and dynamic buckling analyses of viscoelastic composite/sandwich plates: A complex modulus approach. Arch. Appl. Mech. 2011, 81, 1253–1268. [Google Scholar] [CrossRef]
  22. Hart, K.R.; Chia, P.X.L.; Sheridan, L.E. Mechanisms and characterization of impact damage in 2D and 3D woven fiber-reinforced composites. Compos. Part. A 2017, 101, 432–443. [Google Scholar] [CrossRef]
  23. Asaee, Z.; Taheri, F. Experimental and numerical investigation into the influence of stacking sequence on the low-velocity impact response of new 3D FMLs. Compos. Struct. 2016, 140, 136–146. [Google Scholar] [CrossRef]
  24. Shariyat, M.; Hosseini, S.H. Eccentric impact analysis of pre-stressed composite sandwich plates with viscoelastic cores: A novel global-local theory and a refined contact law. Compos. Struct. 2014, 117, 333–345. [Google Scholar] [CrossRef]
  25. Suzuki, Y.; Suzuki, T.; Todoroki, A. Smart lightning protection skin for real-time load monitoring of composite aircraft structures under multiple impacts. Compos. Part. A 2014, 67, 44–54. [Google Scholar] [CrossRef]
  26. Liu, H.; Falzon, B.G.; Tan, W. Predicting the Compression-After-Impact (CAI) strength of damage-tolerant hybrid unidirectional/woven carbon-fibre reinforced composite laminates. Compos. Part A 2018, 105, 189–202. [Google Scholar] [CrossRef]
  27. Liu, C.J.; Todd, M.D.; Zheng, Z.L. A nondestructive method for the pretension detection in membrane structures based on nonlinear vibration response to impact. Struct. Health. Monit. 2018, 17, 67–79. [Google Scholar] [CrossRef]
  28. Zhao, L.; Yang, Y. An impact-based broadband aeroelastic energy harvester for concurrent wind and base vibration energy harvesting. Appl. Energ. 2018, 212, 233–243. [Google Scholar] [CrossRef]
  29. Shahdin, A.; Morlier, J.; Mezeix, L. Evaluation of the impact resistance of various composite sandwich beams by vibration tests. Shock. Vib. 2013, 18, 789–805. [Google Scholar] [CrossRef]
  30. Li, W.; Zheng, H.; Sun, G. The moving least squares based numerical manifold method for vibration and impact analysis of cracked bodies. Eng. Fract. Mech. 2017, 190, 410–434. [Google Scholar] [CrossRef]
  31. Pérez, M.A.; Oller, S.; Felippa, C.A. Micro-mechanical approach for the vibration analysis of CFRP laminates under impact-induced damage. Compos. Part B. Eng. 2015, 83, 306–316. [Google Scholar] [CrossRef] [Green Version]
  32. Miramini, A.; Kadkhodaei, M.; Alipour, A. Analysis of interfacial debonding in shape memory alloy wire-reinforced composites. Smart. Mater. Struct. 2015, 25, 015032. [Google Scholar] [CrossRef]
  33. Chang, M.; Wang, Z.; Liang, W. A visco-hyperelastic model for short fiber reinforced polymer composites: Reinforcement and fracture mechanisms. Text. Res. J. 2017, 88, 2727–2740. [Google Scholar] [CrossRef]
  34. Chang, M.; Wang, Z.; Liang, W. A novel failure analysis of SMA reinforced composite plate based on a strain-rate-dependent model: Low-high velocity impact. J. Mater. Res. Technol. 2018. [Google Scholar] [CrossRef]
  35. Hashin, Z. Failure criteria for unidirectional composites. J. Appl. Mech. 1980, 47, 329–334. [Google Scholar] [CrossRef]
  36. Puck, A.; Schürmann, H. Failure analysis of FRP laminates by means of physically based phenomenological models. Compos. Sci. Technol. 1998, 58, 1045–1067. [Google Scholar] [CrossRef]
  37. Brinson, L.C.; Lammering, R. Finite element analysis of the behavior of shape memory alloys and their applications. Int. J. Solids. Struct. 1993, 30, 3261–3280. [Google Scholar] [CrossRef]
  38. Brinson, L.C. One-dimensional constitutive behavior of shape memory alloys: Thermomechanical derivation with non-constant material functions and redefined martensite internal variable. J. Intell. Mat. Syst. Struct. 1993, 4, 229–242. [Google Scholar] [CrossRef]
Figure 1. A schematic of the shape memory alloy (SMA)-reinforced composite samples: (a) Fixed boundary; (b) Vibrating boundary.
Figure 1. A schematic of the shape memory alloy (SMA)-reinforced composite samples: (a) Fixed boundary; (b) Vibrating boundary.
Materials 12 00007 g001
Figure 2. Modeling the impact test of the composite laminate.
Figure 2. Modeling the impact test of the composite laminate.
Materials 12 00007 g002
Figure 3. The simulation results on the composite laminates at different times (group A1), and the bottom view of the sample after impact: (a) t = 0.0015 s; (b) t = 0.0030 s; (c) t = 0.0045 s; (d) t = 0.0060 s; (e) t = 0.0075 s; (f) t = 0.0090 s; (g) bottom view of sample.
Figure 3. The simulation results on the composite laminates at different times (group A1), and the bottom view of the sample after impact: (a) t = 0.0015 s; (b) t = 0.0030 s; (c) t = 0.0045 s; (d) t = 0.0060 s; (e) t = 0.0075 s; (f) t = 0.0090 s; (g) bottom view of sample.
Materials 12 00007 g003
Figure 4. The middle section of different layers during the impact, group A1: (a) layer 16; (b) layer 8; (c) layer 1.
Figure 4. The middle section of different layers during the impact, group A1: (a) layer 16; (b) layer 8; (c) layer 1.
Materials 12 00007 g004
Figure 5. The simulation results on the composite laminates at different times, group A2: (a) t = 0.0015 s; (b) t = 0.0030 s; (c) t = 0.0045 s; (d) t = 0.0060 s; (e) t = 0.0075 s; (f) t = 0.0090 s; (g) bottom view of sample.
Figure 5. The simulation results on the composite laminates at different times, group A2: (a) t = 0.0015 s; (b) t = 0.0030 s; (c) t = 0.0045 s; (d) t = 0.0060 s; (e) t = 0.0075 s; (f) t = 0.0090 s; (g) bottom view of sample.
Materials 12 00007 g005aMaterials 12 00007 g005b
Figure 6. The middle section of different layers during the impact, group A2: (a) layer 16; (b) layer 8; (c) layer 1.
Figure 6. The middle section of different layers during the impact, group A2: (a) layer 16; (b) layer 8; (c) layer 1.
Materials 12 00007 g006
Figure 7. The simulation results on the SMA-reinforced composite laminates at different times: (a) group A3, t = 0.0050 s; (b) group A4, t = 0.0045 s.
Figure 7. The simulation results on the SMA-reinforced composite laminates at different times: (a) group A3, t = 0.0050 s; (b) group A4, t = 0.0045 s.
Materials 12 00007 g007
Figure 8. The middle section of different layers during the impact: (a) layers in group A3; (b) layers in group A4; (c) the SMA in group A3; (d) the SMA in group A4.
Figure 8. The middle section of different layers during the impact: (a) layers in group A3; (b) layers in group A4; (c) the SMA in group A3; (d) the SMA in group A4.
Materials 12 00007 g008
Figure 9. A comparison between the simulation results and the experimental results. (a) absorbed energy–time history; (b) contact force–time history.
Figure 9. A comparison between the simulation results and the experimental results. (a) absorbed energy–time history; (b) contact force–time history.
Materials 12 00007 g009
Figure 10. The fracture morphology of the top layer of the composite laminate under different amplitudes: (a) Group B1; (b) Group B2; (c) Group B3; (d) Group B4; (e) Group B5; (f) Group B6; (g) Group B7; (h) Group B8.
Figure 10. The fracture morphology of the top layer of the composite laminate under different amplitudes: (a) Group B1; (b) Group B2; (c) Group B3; (d) Group B4; (e) Group B5; (f) Group B6; (g) Group B7; (h) Group B8.
Materials 12 00007 g010
Figure 11. The fracture morphology of the top layer of the SMA-reinforced composite laminate under different amplitudes: (a) Group C1; (b) Group C2; (c) Group C3; (d) Group C4.
Figure 11. The fracture morphology of the top layer of the SMA-reinforced composite laminate under different amplitudes: (a) Group C1; (b) Group C2; (c) Group C3; (d) Group C4.
Materials 12 00007 g011
Figure 12. The middle section of the half model during the impact: (a) group B1; (b) group B3; (c) group B5; (d) group B7; (e) group B4.
Figure 12. The middle section of the half model during the impact: (a) group B1; (b) group B3; (c) group B5; (d) group B7; (e) group B4.
Materials 12 00007 g012aMaterials 12 00007 g012b
Figure 13. The middle section of the half model during the impact: (a) group C2; (b) group C3; (c) group C4.
Figure 13. The middle section of the half model during the impact: (a) group C2; (b) group C3; (c) group C4.
Materials 12 00007 g013
Figure 14. The analysis of the impact resistance of the composite laminate in group B and C: (a) absorbed energy; (b) contact force.
Figure 14. The analysis of the impact resistance of the composite laminate in group B and C: (a) absorbed energy; (b) contact force.
Materials 12 00007 g014
Figure 15. The fracture morphology of the top layer with different frequencies: (a) Group D1; (b) Group D2; (c) Group D3; (d) Group D4; (e) Group D5; (f) Group E1; (g) Group E2; (h) Group E3; (i) Group E4; (j) Group E5.
Figure 15. The fracture morphology of the top layer with different frequencies: (a) Group D1; (b) Group D2; (c) Group D3; (d) Group D4; (e) Group D5; (f) Group E1; (g) Group E2; (h) Group E3; (i) Group E4; (j) Group E5.
Materials 12 00007 g015aMaterials 12 00007 g015b
Figure 16. The cross-section of the half model during the impact: (a) group D4; (b) group E4.
Figure 16. The cross-section of the half model during the impact: (a) group D4; (b) group E4.
Materials 12 00007 g016
Figure 17. The analysis of the impact resistance of the composite laminate in group D and group E: (a) absorbed energy; (b) contact force.
Figure 17. The analysis of the impact resistance of the composite laminate in group D and group E: (a) absorbed energy; (b) contact force.
Materials 12 00007 g017
Figure 18. The statistics of the damage area in different groups: (a) layer 16; (b) layer 8; (c) layer 1.
Figure 18. The statistics of the damage area in different groups: (a) layer 16; (b) layer 8; (c) layer 1.
Materials 12 00007 g018
Figure 19. The statistics of the general damage area: (a) laminate; (b) SMA.
Figure 19. The statistics of the general damage area: (a) laminate; (b) SMA.
Materials 12 00007 g019
Figure 20. A comparative study of the statistical results and the simulation results: (a) Et = 0.01; (b) Favg.
Figure 20. A comparative study of the statistical results and the simulation results: (a) Et = 0.01; (b) Favg.
Materials 12 00007 g020
Table 1. The material parameters of unidirectional glass fiber/epoxy composite laminates.
Table 1. The material parameters of unidirectional glass fiber/epoxy composite laminates.
Mechanical ConstantsValues
Young’s modulus/GPa (E1, E2, E3)55.2, 18.4, 18.4
Poisson’s ratio ( υ 12 , υ 13 , υ 23 )0.27, 0.27, 0.43
Shear modulus/GPa (E12, E13, E23)13.8, 13.8, 13.8
Ultimate tensile stress/MPa (XT, YT, ZT)1656, 73.8, 73.8
Ultimate compressive stress/MPa (XC, YC, ZC)1656, 91.8, 91.8
Ultimate shear stress/MPa (S12, S13, S23)117.6, 117.6, 117.6
Table 2. The experimental groups.
Table 2. The experimental groups.
Stacking SequenceImpact Energy/J
Group A1[0°,90°]832
Group A2[0°,90°]864
Group A3[(0°,90°)4,SMA, (0°,90°)4]32
Group A4[(0°,90°)4,SMA, (0°,90°)4]64
Table 3. The amplitudes and frequencies used in the simulation.
Table 3. The amplitudes and frequencies used in the simulation.
Group B1Group B2Group B3Group B4Group B5Group B6Group B7Group B8
A (m)0.0016−0.00160.0032−0.00320.008−0.0080.016−0.016
f (c/s)1000
SMA
Stacking
NO
[0°,90°]8
Group C1Group C2Group C3Group C4
A (m)0.00160.00320.0080.016
f (c/s)1000
SMA
Stacking
YES
[(0°,90°)4, SMA, (0°,90°)4]
Group D1Group D2Group D3Group D4Group D5
A (m)0.0032
f (c/s)100200500200010,000
SMA
Stacking
NO
[0°,90°]8
Group E1Group E2Group E3Group E4Group E5
A (m)0.0032
f (c/s)100200500200010,000
SMA
Stacking
YES
[(0°,90°)4, SMA, (0°,90°)4]
Table 4. The statistics of absorbed energy and contact force.
Table 4. The statistics of absorbed energy and contact force.
A1B1B2B3B4B5B6A7B8
Emax (J)32323221.6718.274.303.720.300.80
Et= 0.01s (J)28.7030.0430.799.245.904.243.560.300.80
Fmax (N)7.047.747.697.927.358.396.860.861.13
Favg (N)1.931.701.810.830.650.160.130.010.02
A3C1C2C3C4
Emax (J)32323213.203.74
Et= 0.01s (J)25.1428.6031.9713.173.74
Fmax (N)6.8910.7710.5610.528.71
Favg (N)3.322.982.310.520.12
D1D2D3D4D5
Emax (J)323225.386.211.75
Et= 0.01s (J)31.9831.9223.684.691.74
Fmax (N)7.407.237.457.291.44
Favg (N)2.312.361.110.170.06
E1E2E3E4E5
Emax (J)32323211.861.79
Et= 0.01s (J)28.7729.9430.7011.861.79
Fmax (N)10.159.7911.929.022.36
Favg (N)2.952.812.690.460.05

Share and Cite

MDPI and ACS Style

Chang, M.; Kong, F.; Sun, M.; He, J. A Three-Phase Model Characterizing the Low-Velocity Impact Response of SMA-Reinforced Composites under a Vibrating Boundary Condition. Materials 2019, 12, 7. https://doi.org/10.3390/ma12010007

AMA Style

Chang M, Kong F, Sun M, He J. A Three-Phase Model Characterizing the Low-Velocity Impact Response of SMA-Reinforced Composites under a Vibrating Boundary Condition. Materials. 2019; 12(1):7. https://doi.org/10.3390/ma12010007

Chicago/Turabian Style

Chang, Mengzhou, Fangyun Kong, Min Sun, and Jian He. 2019. "A Three-Phase Model Characterizing the Low-Velocity Impact Response of SMA-Reinforced Composites under a Vibrating Boundary Condition" Materials 12, no. 1: 7. https://doi.org/10.3390/ma12010007

APA Style

Chang, M., Kong, F., Sun, M., & He, J. (2019). A Three-Phase Model Characterizing the Low-Velocity Impact Response of SMA-Reinforced Composites under a Vibrating Boundary Condition. Materials, 12(1), 7. https://doi.org/10.3390/ma12010007

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop