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

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.


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.

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.

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: where ε ij and σ ij (i,j = 1, 2, and 3) are the strain and stress in the ij-direction, respectively; c ij 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: 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: 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: 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 σ cr s < σ < σ cr 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 σ cr s and σ cr f 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).

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 σ in,ij (t) of the interphase is a function of strain history, and can be expressed as: ε in t dt (12) where g is the relaxation modulus, t is the new time variable, . ε in (t ) = ∂ε in (t )/∂t is the strain rate. The strain history ε in (t) can be obtained using Boltzmann's superposition principle. Also, the relaxation modulus g can be expressed using the discrete relaxation spectrum, as follows: 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].

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: where A and f are the amplitude and the frequency of the vibration, respectively. As for the vibrating boundary condition (Figure 1b), the movement of the boundaries can be expressed as: where A and f are the amplitude and the frequency of the vibration, respectively.

The Effect of the Fixed Boundary Condition on Impact Resistance
The fixed boundary condition is kept in the experimental process.

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.

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.  Figures 3 and 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.  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.  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. 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.  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. 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. 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.

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.

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. 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.

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. 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.

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 (

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 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.

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, t tot = 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 Figures 10 and 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 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.

Group E1
When applying an amplitude A to control the movement of the boundary, the morphologies of the composite laminate after impact are shown in Figures 10 and 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 and Figure10d. 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 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 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 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. 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. 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. 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.

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 Figures 15a and 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.

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 Figures 15a and 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.

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 Figures 15a and 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.  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. 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.

Statistical Analysis of the Damage State
In Table 4, the maximum energy, Emax, the energy at time 0.01s, 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: 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.  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.

Statistical Analysis of the Damage State
In Table 4, the maximum energy, E max , the energy at time 0.01 s, E t = 0.01 , the maximum force, F max , and the average force, F avg , of different groups are shown. An E max = 32 J denotes a rebound behavior of the impactor. The average force is defined as: where N is the total number of output data of force F i 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. 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.

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: where k a and m i are parameters related to amplitude; and k f and n i 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: (17) where E and m are the absorbed energy and the mass of the impactor, respectively. Considering that F avg = m × ∆v/∆t, the average force is also obtained as: 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).
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.

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.

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.

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 ij stiffness matrix d mt , d mc , d f t , d f c damage variables