Sound Radiation of Orthogonal Antisymmetric Composite Laminates Embedded with Pre-Strained SMA Wires in Thermal Environment

The interest of this article lies in the sound radiation of shape memory alloy (SMA) composite laminates. Different from the traditional method of avoiding resonance sound radiation of composite laminates by means of structural parameter design, this paper explores the potential of adjusting the modal peak of the resonant acoustic radiation by using material characteristics of shape memory alloys (SMA), and provides a new idea for avoiding resonance sound radiation of composite laminates. For composite laminates embedded with pre-strained SMA, an innovation of vibration-acoustic modeling of SMA composite laminates considering pre-stain of SMA and thermal expansion force of graphite-epoxy resin is proposed. Based on the classical thin plate theory and Hamilton principle, the structural dynamic governing equation and the frequency equation of the laminates subjected to thermal environment are derived. The vibration sound radiation of composite laminates is calculated with Rayleigh integral. Effects of ambient temperature, pre-strain, SMA volume fraction, substrate ratio, and geometrical parameters on the sound radiation were analyzed. New laws of SMA material and pre-strain characteristics on sound radiation of composite laminates under temperature environment are revealed, which have theoretical and engineering functional significance for vibration and sound radiation control of SMA composite laminates.


Introduction
Composite laminates as structural components have been extensively implemented in various fields of mechanical engineering such as aeronautics, astronautics, naval architecture, and automotive engineering. Composite plate offers advantages such as lightweight, high strength; but resulting in considerable vibration noise. Especially when the external excitation frequency is equivalent to or close to the natural frequency of the composite laminates, the strong structural resonance leads to more radiated noise. Shape memory alloy (SMA) reinforced composites are an extremely versatile class of materials which have the characteristics of larger internal forces, unique ability of changing its material properties, wide range of operational temperature, and high durability [1]. Though embedded with SMA wire to the composite material, the modal performance of composite laminates can be adjusted by the volume fraction and pre-strain of SMA without modifying the original shape dimensions and boundary conditions. This modification can effectively avoid the structure resonance and minimize the noise radiated into the surroundings.
Acoustic radiation of composite laminates is intimately related with the structural mode and vibration response. Over recent decades, the vibration of SMA composite plates has attracted the Materials 2020, 13, 3657 2 of 18 attention of some scholars. Rogers et al. [2] proposed utilizing SMA fibers modified structural performance of laminated plates via theoretically analyzed natural frequencies and mode shapes. Ostachowicz et al. [3] established a finite element model for predicting the first three natural frequencies subjected to several kinds of pre-strain. Malekzadeh et al. [4][5][6][7] developed the linear free vibration of hybrid laminated composite plates embedded with SMA fibers, employing the Ritz method, which lead to extracting the fundamental natural frequency, and investigated the dynamic behavior on detailed parametric analyses such as volume fraction, pre-strain, and aspect ratio. Park et al. [8] studied the vibration behavior of thermally buckled SMA composite plates using the von Karman nonlinear incremental strain-displacement relation. Furthermore, there are some investigations of the SMA composite plates with regard to vibration response, such as nonlinear vibration analysis [9][10][11], prediction of thermo-mechanical response under the combined action of thermal and mechanical loads [12,13], active vibration control of structures [14][15][16][17], and the experiment of structural acoustic characteristics [18].
Few studies about the acoustic radiation of SMA composite laminates are presented in the reported researching. However, research on vibration and acoustic radiation of conventional composite laminates has been very active. The design of material and structure parameters is an important means to improve the acoustic radiation. Comparing and analyzing sound radiation from material parameters of laminates, the conspicuousness of material composition on response, and the coupling effect of material proportion on sound radiation is often of great significance [19][20][21]. In terms of structural parameters, the research on sound radiation has developed from the design of normal geometrical parameters (aspect ratio, thickness) of laminates to more complex structural forms such as periodically reinforced structures and non-uniform laying [22][23][24][25]. In addition, as a current hotspot, the investigation of acoustic radiation of functionally graded laminates comprehensively considers the effects of materials, structures, external excitations, and the environment [26][27][28]. For this study, a new method to estimate sound power of orthogonal antisymmetric composite laminates embedded with SMA considering the pre-strain of SMA and thermal environment was presented and so the thermal expansion force and recovery force are added to the constitutive model of composite laminates. Based on the Hamilton principle and Rayleigh integration, the vibration differential equation and sound radiation calculation model of SMA composite laminates are established. The influence of SMA wires on the sound radiation power of laminates was examined by varying the pre-strain, SMA volume fraction, substrate ratio, and structural size parameters with temperature. The sound radiation results in the study provide a theoretical basis for vibration and sound radiation control of embedded SMA composite laminates.

Material Properties of SMA in Composite Laminates
Consider a rectangular composite laminates of length L a , width L b , and thickness h in a Cartesian coordinate system as shown in Figure 1. attention of some scholars. Rogers et al. [2] proposed utilizing SMA fibers modified structural performance of laminated plates via theoretically analyzed natural frequencies and mode shapes. Ostachowicz et al. [3] established a finite element model for predicting the first three natural frequencies subjected to several kinds of pre-strain. Malekzadeh et al. [4][5][6][7] developed the linear free vibration of hybrid laminated composite plates embedded with SMA fibers, employing the Ritz method, which lead to extracting the fundamental natural frequency, and investigated the dynamic behavior on detailed parametric analyses such as volume fraction, pre-strain, and aspect ratio. Park et al. [8] studied the vibration behavior of thermally buckled SMA composite plates using the von Karman nonlinear incremental strain-displacement relation. Furthermore, there are some investigations of the SMA composite plates with regard to vibration response, such as nonlinear vibration analysis [9][10][11], prediction of thermo-mechanical response under the combined action of thermal and mechanical loads [12,13], active vibration control of structures [14][15][16][17], and the experiment of structural acoustic characteristics [18]. Few studies about the acoustic radiation of SMA composite laminates are presented in the reported researching. However, research on vibration and acoustic radiation of conventional composite laminates has been very active. The design of material and structure parameters is an important means to improve the acoustic radiation. Comparing and analyzing sound radiation from material parameters of laminates, the conspicuousness of material composition on response, and the coupling effect of material proportion on sound radiation is often of great significance [19][20][21]. In terms of structural parameters, the research on sound radiation has developed from the design of normal geometrical parameters (aspect ratio, thickness) of laminates to more complex structural forms such as periodically reinforced structures and non-uniform laying [22][23][24][25]. In addition, as a current hotspot, the investigation of acoustic radiation of functionally graded laminates comprehensively considers the effects of materials, structures, external excitations, and the environment [26][27][28]. For this study, a new method to estimate sound power of orthogonal antisymmetric composite laminates embedded with SMA considering the pre-strain of SMA and thermal environment was presented and so the thermal expansion force and recovery force are added to the constitutive model of composite laminates. Based on the Hamilton principle and Rayleigh integration, the vibration differential equation and sound radiation calculation model of SMA composite laminates are established. The influence of SMA wires on the sound radiation power of laminates was examined by varying the pre-strain, SMA volume fraction, substrate ratio, and structural size parameters with temperature. The sound radiation results in the study provide a theoretical basis for vibration and sound radiation control of embedded SMA composite laminates.

Material Properties of SMA in Composite Laminates
Consider a rectangular composite laminates of length La, width Lb, and thickness h in a Cartesian coordinate system as shown in Figure 1.  Each layer of material is composed of SMA (nickel-titanium alloy), graphite, and epoxy in which SMA fibers are distributed in laminates in the form of orthogonal anti-symmetry. The elastic modulus of the SMA appears to have a strong temperature dependence [29]. From Figure 2, the modulus is 25 to 30 GPa when the temperature is below 40 • C as SMA in the martensitic phase. Each layer of material is composed of SMA (nickel-titanium alloy), graphite, and epoxy in which SMA fibers are distributed in laminates in the form of orthogonal anti-symmetry. The elastic modulus of the SMA appears to have a strong temperature dependence [29]. From Figure 2, the modulus is 25 to 30 GPa when the temperature is below 40 °C as SMA in the martensitic phase. As the temperature increases, the modulus of SMA increases rapidly when it transforms from the martensitic phase to austenite phase. Most of SMA are in the austenite phase when the temperature reaches 60 °C, and its modulus reaches about 80 GPa. After that, the modulus increases slowly with the rising temperature.
Recovery stress produced by SMA is not only temperature dependent, but also tensile pre-strain dependent. It can be seen from Figure 3 that the recovery stress increases with the increase of prestrain after the temperature is higher than 60 °C [29]. Furthermore, the effect of pre-strain is not significant for normal or low temperature environment. To composite laminates embedded with SMA, owing to its temperature dependent material properties and pre-strain characteristics, the thermal environment should be fully considered in the process of its application. Especially through the control of material properties by As the temperature increases, the modulus of SMA increases rapidly when it transforms from the martensitic phase to austenite phase. Most of SMA are in the austenite phase when the temperature reaches 60 • C, and its modulus reaches about 80 GPa. After that, the modulus increases slowly with the rising temperature.
Recovery stress produced by SMA is not only temperature dependent, but also tensile pre-strain dependent. It can be seen from Figure 3 that the recovery stress increases with the increase of pre-strain after the temperature is higher than 60 • C [29]. Each layer of material is composed of SMA (nickel-titanium alloy), graphite, and epoxy in which SMA fibers are distributed in laminates in the form of orthogonal anti-symmetry. The elastic modulus of the SMA appears to have a strong temperature dependence [29]. From Figure 2, the modulus is 25 to 30 GPa when the temperature is below 40 °C as SMA in the martensitic phase. As the temperature increases, the modulus of SMA increases rapidly when it transforms from the martensitic phase to austenite phase. Most of SMA are in the austenite phase when the temperature reaches 60 °C, and its modulus reaches about 80 GPa. After that, the modulus increases slowly with the rising temperature.
Recovery stress produced by SMA is not only temperature dependent, but also tensile pre-strain dependent. It can be seen from Figure 3 that the recovery stress increases with the increase of prestrain after the temperature is higher than 60 °C [29]. Furthermore, the effect of pre-strain is not significant for normal or low temperature environment. To composite laminates embedded with SMA, owing to its temperature dependent material properties and pre-strain characteristics, the thermal environment should be fully considered in the process of its application. Especially through the control of material properties by Furthermore, the effect of pre-strain is not significant for normal or low temperature environment. To composite laminates embedded with SMA, owing to its temperature dependent material properties and pre-strain characteristics, the thermal environment should be fully considered in the process of its application. Especially through the control of material properties by adjusting the temperature, which has great potential for improving the vibration and sound radiation of the composite laminates.

Structural Dynamic Response
Assuming that the SMA orthogonal anti-symmetric composite laminate shown in Figure 1 contains N layers, the stress-strain relation of the kth layer can be given as [30], Equation (1): where Q (k) is the transformed reduced stiffness of the kth layer of composite laminate embedded with SMA and Q (k) is the transformed reduced stiffness of the kth layer of the substrate which is graphite-epoxy. The term σ r denotes the recovery stress in the SMA fiber direction. Other equations and terms in Equation (1) are completed in Appendix A. Utilizing the classic thin plate theory, the strains associated with the displacements are, Equation (2): where 'u 0 ', 'v 0 ' are the displacement components of the middle plane along the x, y coordinates, respectively and 'w' is the displacement components of along z coordinate. Then, the strain energy U 1 is obtained as, Equation (3): The additional strain energy U 2 which consists of the strain energy caused by recovery stress and the thermal deformation energy can be expressed as, Equation (4): where Equations (5) and (6): The damping deformation energy U 3 can be given as, Equation (7): Materials 2020, 13, 3657 w and λ is the viscous damping coefficient. The total potential energy of composite laminates can be written as, Equation (8): The kinetic energy of composite laminates can be calculated as, Equation (9): It is noteworthy in Equation (9) that the sound radiation of the composite laminates is mainly generated by bending vibration, so that the kinetic energy in the x and y directions can be neglected. For a transverse load q(x, y, t) on the surface of the laminates, the work it has done on the system can be obtained as, Equation (10): According to the Hamilton's principle, the dynamic equations of the laminates can be derived by Equation (11): Substituting Equations (8)-(10) into Equation (11), meanwhile considering the condition of orthogonal antisymmetric composite laminates, only B 11 and B 22 are not zero in B ij , and B 22 = −B 11 . Then, the vibration equation of an orthogonal antisymmetric composite laminate embedded in SMA can be derived as, Equations (12)-(14): where Equations (15)- (18): In this paper, mechanical boundary conditions of composite laminates are intended to be simply supported, Equation (19): Materials 2020, 13, 3657 6 of 18 The state variables satisfying the simply supported boundary conditions of Equation (19) are assumed as Equation (20): where U mn , V mn , and W mn are the Fourier expansion coefficients of the solution, m and n are the half wave numbers in x and y directions, and ω is the circular frequency. Similarly, the load q(x, y, t) is assumed to be expanded into a double Fourier series form, Equation (21): where q mn as the load coefficient can be given as, Equation (22): Substituting Equations (20) and (21) into Equations (12)- (14) and the vibration equation can be rewritten into a form as, Equation (23): where {X mn } = {U mn , V mn , W mn } T and {F mn } = 0, 0, q mn T are displacement amplitude vector and load amplitude vector, respectively. Moreover, the mass matrix, damping matrix, and stiffness matrix is respectively expressed as, Equation (24): in which the damping loss factor is, Equation (25): and K ij are given as, Equations (26)-(31): Materials 2020, 13, 3657 7 of 18 Multiply both sides of the Equation (23) by the inverse matrix of the coefficient matrix, Equation (23): Solving Equation (32) obtains displacement amplitude in z direction of composite laminates, Equation (33): where ω mn is the natural frequency which can be solved by the following frequency Equations (34): Substituting Equation (33) into the third equation of Equation (20), it can be obtained that the displacement components along z coordinate, Equation (35): The derivative of the displacement in Equation (35) with respect to time is utilized to obtain the amplitude of the (m, n) order modal vibration velocity of composite laminates, Equations (36):

Sound Radiation Power
For rectangular composite laminates embedded with SMA, each unit in the laminates can be regarded as a source of external vibrational acoustic radiation. The radiation sound pressure of the laminates is obtained by radiating the sound pressure of each unit, and then performing Rayleigh integration. Assuming that the amplitude of the vibration velocity of the ds in the laminates iŝ U mn (x , y ), the amplitude of the sound pressure at ds can be formulated as, Equation (37): where k, ρ 0 , and c 0 are the sound wave number, the density of the air, and the sound velocity, respectively. r is the distance from unit ds to unit ds. The sound pressure and sound intensity of all panel elements are radiated as being defined by [30], Equations (38) and (39): whereÛ * mn (x, y) is the conjugate ofÛ mn (x, y). Using Equation (39), the radiated sound power of the composite laminate can be formulated as [26], Equation (40): The radiated sound power is usually written in the form of sound power level in decibel, which is defined by, Equation (41): whereŴ 0 is the reference power andŴ 0 = 1 × 10 −12 W.

Materials and Verification
The analytical method is thus deployed to carry out several parametric studies to explore the acoustic response of composite laminates in thermal environments. The rectangular composite laminates are composed of ten layers, simply supported on all edges with dimensions of 0.4 m × 0.3 m × 0.008 m are considered as shown in Figure 1. Young's modulus of elasticity and recovery stress for SMA were taken from Figures 2 and 3, respectively. Poisson's ratio, density, and thermal expansion coefficient of SMA are 0.3, 645 kg·m 3 , and 10.26 × 10 −6 1/ • C, respectively. SMA is embedded in each layer with a volume fraction of 20%; graphite and epoxy are invoked as the substrate for a ratio of 1/9. The material properties of the substrate are given in Table 1.  In addition, the stacking sequence is [0 • /90 • /0 • /90 • /0 • ] a and the damping loss factor of the laminates is set at 0.01. To verify the code of sound radiation characteristics of composite laminates embedded with SMA, the natural frequency and the sound radiation were calculated by the present method. The first example is the natural frequency from a rectangular laminated composite plate with SMA, which is taken from Park et al. [8]. The size of the plate is 0.38 m × 0.305 m × 0.002 m with simply supported boundary conditions for all edges. Table 2 shows the comparison of the first 3 natural frequencies of the laminated composite plate subjected to temperature at 25, 45, 65, 100, and 120 • C. As can be seen from Table 2, the comparison certifies the correctness of natural frequencies that was calculated by the present method. The second example is sound radiation from antisymmetric cross-ply [0 • /90 • ] 2 laminated composite flat panel (0.5 m × 0.5 m × 0.02 m) subjected to unit harmonic point load at x = 0.125 m and y = 0.125 m which is taken from Sharma et al. [20]. Table 3 represents the material properties of the graphite-epoxy and glass-epoxy. As shown in Figure 4, the calculated result has an excellent agreement with Sharma et al. [20]. Thus, the present analytical method was a sound choice for the study.  [20]. Table 3 represents the material properties of the graphite-epoxy and glass-epoxy. As shown in Figure 4, the calculated result has an excellent agreement with Sharma et al. [20]. Thus, the present analytical method was a sound choice for the study. .

Investigating the Influence of Temperature Dependent Material Properties
Composite laminates embedded with SMA are temperature dependent. In order to depict the acoustic radiation power of the composite laminates in the thermal field, it assumes to be a uniform heating process. Figure 5 shows a three-dimensional diagram of the acoustic radiation power of composite laminates in the frequency range of 10 to 1000 Hz during the temperature rise from 15 to 155 ℃.

Investigating the Influence of Temperature Dependent Material Properties
Composite laminates embedded with SMA are temperature dependent. In order to depict the acoustic radiation power of the composite laminates in the thermal field, it assumes to be a uniform heating process. Figure 5 shows a three-dimensional diagram of the acoustic radiation power of composite laminates in the frequency range of 10 to 1000 Hz during the temperature rise from 15 to 155 • C.
As shown in Figure 5, it can be observed that there are four to five peaks of the sound power level corresponding to the natural frequencies. First and second modes shift to lower frequency and approach to each other with the heating case. Simultaneously, the frequency gap between the second mode and the third mode is widened. Moreover, the value of the sound radiation power level also varies and fewer peaks occur in the frequency range showed owing to the effect of heating. These significant effects stem from the temperature dependent material properties of the modulus and pre-strain of the SMA. To investigate further sound radiation characteristics of composite laminates embedded with SMA, representative temperatures such as 25 • C (SMA in martensitic phase) and 100 • C (SMA in austenite phase) were selected as conditions for subsequent study. As shown in Figure 5, it can be observed that there are four to five peaks of the sound power level corresponding to the natural frequencies. First and second modes shift to lower frequency and approach to each other with the heating case. Simultaneously, the frequency gap between the second mode and the third mode is widened. Moreover, the value of the sound radiation power level also varies and fewer peaks occur in the frequency range showed owing to the effect of heating. These significant effects stem from the temperature dependent material properties of the modulus and prestrain of the SMA. To investigate further sound radiation characteristics of composite laminates embedded with SMA, representative temperatures such as 25 ℃ (SMA in martensitic phase) and 100 ℃ (SMA in austenite phase) were selected as conditions for subsequent study.

Investigating the Influence of SMA Pre-Strain
The influence of SMA pre-strain on the acoustic radiation of composite laminates is studied. Ensure that the substrate material ratio and SMA volume fraction are constant. The tensile pre-strains of SMA are treated as 1%, 3%, and 5%, respectively, to calculate the acoustic radiation power of composite laminates. Figures 6 and 7 show the effects of pre-strain on the sound radiation of the composite laminates with SMA in 25 and 100 ℃, respectively.  Sound radiation power of composite laminates embedded with SMA in the thermal environment.

Investigating the Influence of SMA Pre-Strain
The influence of SMA pre-strain on the acoustic radiation of composite laminates is studied. Ensure that the substrate material ratio and SMA volume fraction are constant. The tensile pre-strains of SMA are treated as 1%, 3%, and 5%, respectively, to calculate the acoustic radiation power of composite laminates. Figures 6 and 7 show the effects of pre-strain on the sound radiation of the composite laminates with SMA in 25 and 100 • C, respectively. As shown in Figure 5, it can be observed that there are four to five peaks of the sound power level corresponding to the natural frequencies. First and second modes shift to lower frequency and approach to each other with the heating case. Simultaneously, the frequency gap between the second mode and the third mode is widened. Moreover, the value of the sound radiation power level also varies and fewer peaks occur in the frequency range showed owing to the effect of heating. These significant effects stem from the temperature dependent material properties of the modulus and prestrain of the SMA. To investigate further sound radiation characteristics of composite laminates embedded with SMA, representative temperatures such as 25 ℃ (SMA in martensitic phase) and 100 ℃ (SMA in austenite phase) were selected as conditions for subsequent study.

Investigating the Influence of SMA Pre-Strain
The influence of SMA pre-strain on the acoustic radiation of composite laminates is studied. Ensure that the substrate material ratio and SMA volume fraction are constant. The tensile pre-strains of SMA are treated as 1%, 3%, and 5%, respectively, to calculate the acoustic radiation power of composite laminates. Figures 6 and 7 show the effects of pre-strain on the sound radiation of the composite laminates with SMA in 25 and 100 ℃, respectively.  No considerable changes are seen in the sound power level under different pre-strain at 25 • C; however, the peaks of sound power level shift towards high frequencies with the increase of pre-strain at 100 • C. From Figures 2 and 3, SMA in the austenite phase has higher modulus and recovery stress than that in the martensitic phase, which can explain the fact that extremely small recovery stress is generated by pre-strain at normal temperature since SMA is in the martensitic phase; but recovery stress increases gradually with rising temperature, which is due to the transformation of the martensitic phase to austenite phase. Moreover, it can be known from Equation (31) that the recovery stress has a positive correlation with the coefficient K 33 in the stiffness matrix. Therefore, the influence of recovery stress produced by tensile pre-strain of SMA at 25 • C can be ignored; nevertheless, at 100 • C, the increase of the stiffness of the composite laminates with the increase of recovery stress, results in the increase of the corresponding natural frequencies. No considerable changes are seen in the sound power level under different pre-strain at 25 ℃; however, the peaks of sound power level shift towards high frequencies with the increase of prestrain at 100 ℃. From Figures 2 and 3, SMA in the austenite phase has higher modulus and recovery stress than that in the martensitic phase, which can explain the fact that extremely small recovery stress is generated by pre-strain at normal temperature since SMA is in the martensitic phase; but recovery stress increases gradually with rising temperature, which is due to the transformation of the martensitic phase to austenite phase. Moreover, it can be known from Equation (31) that the recovery stress has a positive correlation with the coefficient in the stiffness matrix. Therefore, the influence of recovery stress produced by tensile pre-strain of SMA at 25 ℃ can be ignored; nevertheless, at 100 ℃, the increase of the stiffness of the composite laminates with the increase of recovery stress, results in the increase of the corresponding natural frequencies.

Investigating the Influence of SMA Volume Fraction
Materials volume fraction is an important component parameter of the plate materials for the acoustic properties. This section considers the effect of SMA volume fractions of 10%, 20%, 30%, and 40% on the sound radiation of composite laminates.
The variation of sound radiation power with different SMA volume fraction of the composite laminates is depicted in Figures 8 and 9.

Investigating the Influence of SMA Volume Fraction
Materials volume fraction is an important component parameter of the plate materials for the acoustic properties. This section considers the effect of SMA volume fractions of 10%, 20%, 30%, and 40% on the sound radiation of composite laminates.
The variation of sound radiation power with different SMA volume fraction of the composite laminates is depicted in Figures 8 and 9. No considerable changes are seen in the sound power level under different pre-strain at 25 ℃; however, the peaks of sound power level shift towards high frequencies with the increase of prestrain at 100 ℃. From Figures 2 and 3, SMA in the austenite phase has higher modulus and recovery stress than that in the martensitic phase, which can explain the fact that extremely small recovery stress is generated by pre-strain at normal temperature since SMA is in the martensitic phase; but recovery stress increases gradually with rising temperature, which is due to the transformation of the martensitic phase to austenite phase. Moreover, it can be known from Equation (31) that the recovery stress has a positive correlation with the coefficient in the stiffness matrix. Therefore, the influence of recovery stress produced by tensile pre-strain of SMA at 25 ℃ can be ignored; nevertheless, at 100 ℃, the increase of the stiffness of the composite laminates with the increase of recovery stress, results in the increase of the corresponding natural frequencies.

Investigating the Influence of SMA Volume Fraction
Materials volume fraction is an important component parameter of the plate materials for the acoustic properties. This section considers the effect of SMA volume fractions of 10%, 20%, 30%, and 40% on the sound radiation of composite laminates.
The variation of sound radiation power with different SMA volume fraction of the composite laminates is depicted in Figures 8 and 9.  The plot Figure 8 reveals that at 25 • C, the modes shift towards low frequencies with an increase of SMA volume fraction. Modulus of SMA in the martensitic phase at ordinary temperature is less than that of substrate modulus as shown in Figure 2. The increase in SMA volume fraction leads to a decrease in the overall stiffness of the laminates, which bring about the decrease of natural frequencies.
When the temperature is 100 • C, the modulus of SMA in austenite phase is greater than the substrate. As the SMA volume fraction increases, the stiffness of the laminates increases, so the first two modes shift towards high frequencies. Besides, third and fourth modes shift slightly towards low frequencies. Materials 2020, 13, x FOR PEER REVIEW 11 of 17 Figure 9. Sound radiation power of composite laminates embedded with SMA subjected to different SMA volume fraction at 100 ℃.
The plot Figure 8 reveals that at 25 ℃, the modes shift towards low frequencies with an increase of SMA volume fraction. Modulus of SMA in the martensitic phase at ordinary temperature is less than that of substrate modulus as shown in Figure 2. The increase in SMA volume fraction leads to a decrease in the overall stiffness of the laminates, which bring about the decrease of natural frequencies. When the temperature is 100 ℃, the modulus of SMA in austenite phase is greater than the substrate. As the SMA volume fraction increases, the stiffness of the laminates increases, so the first two modes shift towards high frequencies. Besides, third and fourth modes shift slightly towards low frequencies.

Investigating the Influence of Substrate Ratio
Considering the larger modulus of graphite in the substrate, the volume fraction ratio of graphite to epoxy studied in this section determines the modulus of the substrate. The volume fraction of SMA is 20%, so the volume fraction of substrate is 80%. A comparative study was carried out with the volume fraction ratio of graphite to epoxy (Vf/Vm) as 1/9, 2/8, 3/7, and 4/6.
The result in Figures 10 and 11 illustrate influence of the substrate ratio on the sound radiation of composite laminates.

Investigating the Influence of Substrate Ratio
Considering the larger modulus of graphite in the substrate, the volume fraction ratio of graphite to epoxy studied in this section determines the modulus of the substrate. The volume fraction of SMA is 20%, so the volume fraction of substrate is 80%. A comparative study was carried out with the volume fraction ratio of graphite to epoxy (V f /V m ) as 1/9, 2/8, 3/7, and 4/6.
The result in Figures 10 and 11 illustrate influence of the substrate ratio on the sound radiation of composite laminates. The plot Figure 8 reveals that at 25 ℃, the modes shift towards low frequencies with an increase of SMA volume fraction. Modulus of SMA in the martensitic phase at ordinary temperature is less than that of substrate modulus as shown in Figure 2. The increase in SMA volume fraction leads to a decrease in the overall stiffness of the laminates, which bring about the decrease of natural frequencies. When the temperature is 100 ℃, the modulus of SMA in austenite phase is greater than the substrate. As the SMA volume fraction increases, the stiffness of the laminates increases, so the first two modes shift towards high frequencies. Besides, third and fourth modes shift slightly towards low frequencies.

Investigating the Influence of Substrate Ratio
Considering the larger modulus of graphite in the substrate, the volume fraction ratio of graphite to epoxy studied in this section determines the modulus of the substrate. The volume fraction of SMA is 20%, so the volume fraction of substrate is 80%. A comparative study was carried out with the volume fraction ratio of graphite to epoxy (Vf/Vm) as 1/9, 2/8, 3/7, and 4/6.
The result in Figures 10 and 11 illustrate influence of the substrate ratio on the sound radiation of composite laminates.  It is observed that the corresponding peaks of the sound power level shift towards the higher frequency domain when the volume fraction of graphite increases in Figure 10. The reason is that graphite has a larger modulus than that SMA in the martensitic phase which significantly affects the stiffness of the laminates at 25 • C. As can be seen from Figure 11, at 100 • C, increasing the volume fraction of graphite in the substrate has little effect on the first mode of the composite laminates, which is due to the influence of the modulus of SMA in the austenite phase on the stiffness increase. However, after the second mode, as the proportion of graphite increases, the peaks of sound power gradually shift towards the high frequency, i.e., the higher the volume fraction of graphite in the substrate, the higher the high-order modal frequency can be increased. Materials 2020, 13, x FOR PEER REVIEW 12 of 17 Figure 11. Sound radiation power of composite laminates embedded with SMA subjected to different substrate ratio at 100 ℃.
It is observed that the corresponding peaks of the sound power level shift towards the higher frequency domain when the volume fraction of graphite increases in Figure 10. The reason is that graphite has a larger modulus than that SMA in the martensitic phase which significantly affects the stiffness of the laminates at 25 ℃. As can be seen from Figure 11, at 100 ℃, increasing the volume fraction of graphite in the substrate has little effect on the first mode of the composite laminates, which is due to the influence of the modulus of SMA in the austenite phase on the stiffness increase. However, after the second mode, as the proportion of graphite increases, the peaks of sound power gradually shift towards the high frequency, i.e., the higher the volume fraction of graphite in the substrate, the higher the high-order modal frequency can be increased.

Investigating the Influence of Geometrical Properties
The geometric properties of composite laminates are studied from the thickness and aspect ratio.

Investigating the Influence of Geometrical Properties
The geometric properties of composite laminates are studied from the thickness and aspect ratio. Figures 12 and 13   It is observed that the corresponding peaks of the sound power level shift towards the higher frequency domain when the volume fraction of graphite increases in Figure 10. The reason is that graphite has a larger modulus than that SMA in the martensitic phase which significantly affects the stiffness of the laminates at 25 ℃. As can be seen from Figure 11, at 100 ℃, increasing the volume fraction of graphite in the substrate has little effect on the first mode of the composite laminates, which is due to the influence of the modulus of SMA in the austenite phase on the stiffness increase. However, after the second mode, as the proportion of graphite increases, the peaks of sound power gradually shift towards the high frequency, i.e., the higher the volume fraction of graphite in the substrate, the higher the high-order modal frequency can be increased.

Investigating the Influence of Geometrical Properties
The geometric properties of composite laminates are studied from the thickness and aspect ratio.   From Figures 12 and 13, the trend of mode shift caused by thickness variation is consistent for 25 and 100 • C. With the increase of thickness, the first and second mode shifts towards high frequencies; fewer peaks occur in the frequency range showed for thickness h > 0.012 m in both normal and high temperature, which is due to the fact that some higher-order modes have shifted above 1000 Hz. In addition, the increase of the thickness of composite laminates leads to a reduction of sound radiation, which is significantly below 400 Hz; however, the thickness has no significant effect on sound radiation in the high frequency region. From Figures 12 and 13, the trend of mode shift caused by thickness variation is consistent for 25 and 100 ℃ . With the increase of thickness, the first and second mode shifts towards high frequencies; fewer peaks occur in the frequency range showed for thickness h > 0.012 m in both normal and high temperature, which is due to the fact that some higher-order modes have shifted above 1000 Hz. In addition, the increase of the thickness of composite laminates leads to a reduction of sound radiation, which is significantly below 400 Hz; however, the thickness has no significant effect on sound radiation in the high frequency region.
Similar to the plate thickness, aspect ratio is another structural geometric parameter of composite laminates. In order to study the effect of increasing aspect ratio on sound radiation, the length of the laminates is increased and the width is shortened while maintaining a constant area, i.e., 0.  Figure 14 and Figure 15 indicate that the aspect ratio has a significant effect on the acoustic radiation, in which the first mode shifts towards high frequencies with the increase of the aspect ratio.  Similar to the plate thickness, aspect ratio is another structural geometric parameter of composite laminates. In order to study the effect of increasing aspect ratio on sound radiation, the length of the laminates is increased and the width is shortened while maintaining a constant area, i.e., 0.  Figures 14 and 15 indicate that the aspect ratio has a significant effect on the acoustic radiation, in which the first mode shifts towards high frequencies with the increase of the aspect ratio. From Figures 12 and 13, the trend of mode shift caused by thickness variation is consistent for 25 and 100 ℃ . With the increase of thickness, the first and second mode shifts towards high frequencies; fewer peaks occur in the frequency range showed for thickness h > 0.012 m in both normal and high temperature, which is due to the fact that some higher-order modes have shifted above 1000 Hz. In addition, the increase of the thickness of composite laminates leads to a reduction of sound radiation, which is significantly below 400 Hz; however, the thickness has no significant effect on sound radiation in the high frequency region.
Similar to the plate thickness, aspect ratio is another structural geometric parameter of composite laminates. In order to study the effect of increasing aspect ratio on sound radiation, the length of the laminates is increased and the width is shortened while maintaining a constant area, i.e., 0.  Figure 14 and Figure 15 indicate that the aspect ratio has a significant effect on the acoustic radiation, in which the first mode shifts towards high frequencies with the increase of the aspect ratio.  The foremost reason is that the bending stiffness of composite laminates increases with the ratio of length to width under the condition of a certain area, which increases the modal frequency.
As mentioned above, the observations that the variation of thickness and aspect ratio affects sound radiation is determined by the variation of structural stiffness. However, temperature changes stiffness of laminates by affecting material properties, which is not a factor influencing the stiffness from structural geometric dimension. Therefore, the variation of structural geometric parameters has no remarkable law on acoustic radiation between normal and high temperature, but the trends are comparable. The foremost reason is that the bending stiffness of composite laminates increases with the ratio of length to width under the condition of a certain area, which increases the modal frequency.
As mentioned above, the observations that the variation of thickness and aspect ratio affects sound radiation is determined by the variation of structural stiffness. However, temperature changes stiffness of laminates by affecting material properties, which is not a factor influencing the stiffness from structural geometric dimension. Therefore, the variation of structural geometric parameters has no remarkable law on acoustic radiation between normal and high temperature, but the trends are comparable.

Conclusions
Vibro-acoustic radiation of composite laminates which embeds SMA orthogonal antisymmetric has been developed. The constitutive equation, energy method, and Hamilton principle are utilized to derive the governing equations, and sound radiation is calculated by employing the Rayleigh integral approach. From the current investigation, the following conclusions are made: (1) The pre-strain and volume fraction of SMA in composite laminates have a significant impact on vibro-acoustic radiation. The higher the pre-strain of the SMA at high temperatures, the modes shift towards high frequencies. However, the volume fraction of SMA has opposite effects on the modalities of acoustic radiation between low and high temperatures. That is, as the volume fraction of SMA increases, modes shift towards high frequencies at elevated temperatures, and the opposite occurs at low temperatures. (2) Graphite in the base material has a larger modulus which affects the stiffness of the laminates.
The effect of increasing volume fraction of graphite at low temperature is opposite to that of SMA on sound radiation. At this point, the modulus of graphite is larger than that of SMA, at which point the graphite content dominates the overall stiffness of the laminates. (3) Several geometric parameter studies on the effects of geometrical properties on sound radiation of composite laminates are carried out. It is noticed that composite laminates embedded SMA have similar characteristics of sound radiation with homogeneous laminates in the temperature range. Influence of geometric parameters on sound radiation is mainly dependent on structural stiffness rather than material properties, resulting in insensitivity to temperature.
The above conclusion fully demonstrates that the composite laminates embedded with SMA wires have excellent ability to adjust the mode frequency and change sound radiation characteristics. This study provides a basis for avoiding modal resonance and reducing vibration and sound radiation of composite laminates from the point of view of material design.

Conclusions
Vibro-acoustic radiation of composite laminates which embeds SMA orthogonal antisymmetric has been developed. The constitutive equation, energy method, and Hamilton principle are utilized to derive the governing equations, and sound radiation is calculated by employing the Rayleigh integral approach. From the current investigation, the following conclusions are made: (1) The pre-strain and volume fraction of SMA in composite laminates have a significant impact on vibro-acoustic radiation. The higher the pre-strain of the SMA at high temperatures, the modes shift towards high frequencies. However, the volume fraction of SMA has opposite effects on the modalities of acoustic radiation between low and high temperatures. That is, as the volume fraction of SMA increases, modes shift towards high frequencies at elevated temperatures, and the opposite occurs at low temperatures. (2) Graphite in the base material has a larger modulus which affects the stiffness of the laminates.
The effect of increasing volume fraction of graphite at low temperature is opposite to that of SMA on sound radiation. At this point, the modulus of graphite is larger than that of SMA, at which point the graphite content dominates the overall stiffness of the laminates. (3) Several geometric parameter studies on the effects of geometrical properties on sound radiation of composite laminates are carried out. It is noticed that composite laminates embedded SMA have similar characteristics of sound radiation with homogeneous laminates in the temperature range. Influence of geometric parameters on sound radiation is mainly dependent on structural stiffness rather than material properties, resulting in insensitivity to temperature.
The above conclusion fully demonstrates that the composite laminates embedded with SMA wires have excellent ability to adjust the mode frequency and change sound radiation characteristics. This study provides a basis for avoiding modal resonance and reducing vibration and sound radiation of composite laminates from the point of view of material design.