Shock and Vibration
Volume 2015 (2015), Article ID 523723, 10 pages
http://dx.doi.org/10.1155/2015/523723

Research Article

Nonlinear Dynamic Characteristics and Optimal Control of SMA Composite Wings Subjected to Stochastic Excitation

Zhi-Wen Zhu,¹ Xin-Miao Li,¹ and Jia Xu²

¹Department of Mechanics, Tianjin University, Tianjin 300072, China
²Tianjin Key Laboratory of Nonlinear Dynamics and Chaos Control, Tianjin 300072, China

Received 31 May 2015; Revised 23 June 2015; Accepted 25 June 2015

Academic Editor: Nerio Tullini

Copyright © 2015 Zhi-Wen Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A kind of high-aspect-ratio shape memory alloy (SMA) composite wing is proposed to reduce the wing’s fluttering. The nonlinear dynamic characteristics and optimal control of the SMA composite wings subjected to in-plane stochastic excitation are investigated where the great bending under the flight loads is considered. The stochastic stability of the system is analyzed, and the system’s response is obtained. The conditions of stochastic Hopf bifurcation are determined, and the probability density of the first-passage time is obtained. Finally, the optimal control strategy is proposed. Numerical simulation shows that the stability of the system varies with bifurcation parameters, and stochastic Hopf bifurcation appears in the process; the reliability of the system is improved through optimal control, and the first-passage time is delayed. Finally, the effects of the control strategy are proved by experiments. The results of this paper are helpful for engineering applications of SMA.

1. Introduction

High-aspect-ratio wings are applied in uninhabited air vehicles (UAV) widely. They have a great bending and torsion deformation under the flight loads, which induce the wing’s fluttering. In this paper, a kind of shape memory alloy (SMA) composite wing is proposed to reduce the bending deformation. Shape memory alloy is a kind of smart material. It has many special properties, such as shape memory effect and damping characteristics. Moreover, SMA can be controlled by heating and cooling to achieve the structure deformation, which is helpful for reducing the wing’s fluttering and ensuring the UAV’s safety.

Several scholars have studied the application of shape memory materials in the wings. Yang et al. analyzed the characteristics of smart composite wing with SMA actuators [1]. Mirone designed adaptive airfoils and hingeless wings actuated by shape memory alloy wires [2]. Icardi and Ferrero studied an adaptive wing with shape memory alloy torsion actuators [3, 4]. Barbarino et al. controlled the wing through an SMA-based device [5]. Grigorie and Botez proposed a new adaptive controller method for SMA hysteresis modelling of a morphing wing [6]. Kim and Cho designed a simple smart wing actuator using Ni-Ti SMA [7]. Barzegari et al. discussed the aeroelastic characteristics of cantilever wing with embedded shape memory alloys [8]. Kang et al. analyzed the morphing wing mechanism using an SMA wire actuator [9]. Colorado et al. researched the biomechanics of smart wings in a bat robot using SMA actuators [10]. Although much progress has been reported, the theoretical results of the dynamic characteristics of SMA composite wings are limited. For controlling the wing’s fluttering efficiently, it is necessary to obtain the response and the dynamic characteristics of SMA composite wings.

A high-aspect-ratio SMA composite wing can be regarded as a cantilever SMA composite beam. There are also many researchers studying the dynamic characteristics of SMA beam. Lau discussed firstly the vibration characteristics of SMA beams with different boundary conditions [11]. Liew et al. studied the pseudoelastic behavior of an SMA beam by the element-free Galerkin method [12]. Zbiciak analyzed dynamic characteristics of pseudoelastic SMA beam [13]. Scarpa et al. developed spectral element formulation for SMA beams under random vibration excitation [14]. Hashemi and Khadem developed the dynamic model of SMA beam [15]. Collet et al. analyzed vibration behavior of SMA beam under dynamical loading [16]. However, the airflow force is usually random to a low-flying UAV. It means that the UAV’s high-aspect-ratio SMA composite wings are subjected to stochastic excitation. Although many achievements have been obtained, the dynamic characteristics of a cantilever SMA composite beam subjected to stochastic excitation are not abundant.

This paper aims to offer a kind of analysis method to the nonlinear dynamical characteristic of a cantilever SMA composite beam subjected to stochastic excitation. Nonlinear differential items are introduced to explain the hysteretic phenomenon of strain-stress curve of SMA, and the hysteretic nonlinear dynamic model of a cantilever SMA composite beam subjected to stochastic excitation is developed. The stochastic stability of the system is analyzed, and the system’s response is obtained. The conditions of stochastic Hopf bifurcation are determined, and the probability density of the first-passage time is obtained. Finally, the optimal control strategy to improve the system’s reliability is proposed.

2. Modelling of SMA Material

For obtaining the dynamic characteristics of a Ti-Ni SMA composite wing subjected to stochastic excitation, it is necessary to establish an accurate constitutive model of Ti-Ni SMA. The strain-stress curves of Ti-Ni SMA are presented in Figure 1, which are obtained by experiments. The SMA thin strip we used is supplied by Tianxin Company, Gansu Province, China. Its length is 12 cm, its width is 0.35 cm, and its thickness is 0.4 mm. Evidently, hysteretic phenomena are present in Ti-Ni SMA strain-stress curves. Most of the current SMA models are shown as equations with piecewise functions or double integral functions owing to the hysteretic characteristics of Ti-Ni SMA [16–22]. As a result, it is hard to analyze the dynamic characteristics of SMA structure in theory. Usually, research results based on those models can only be obtained by numerical simulation or experiment [23–27]. In [28], we proposed a Van de Pol hysteretic model about the relationship between the stress and the strain of SMA. In this paper, nonlinear differential items, which are developed from Van de Pol hysteretic model, are introduced to explain the hysteretic loop of Ti-Ni SMA strain-stress curves.

Figure 1: Strain-stress curves of Ti-Ni SMA.

Van der Pol equation is a nonlinear equation. It can be shown as follows:

The item is called Van der Pol item. becomes when is . It means that the essences of Van der Pol item are two parabolic lines symmetrical about the original point () when is constant. Initial Van der Pol model can only be used in some basic parabolic hysteretic loop, where means the loading process and means the unloading process. For those hysteresis loops whose outset is (), Van der Pol item can be deformed as , where and are the parameters which determine the width and height of the hysteretic loop [28].

However, the strain-stress curves of SMA are not pure parabolic lines. In this paper, a new differential item, which is developed from Van der Pol hysteretic model, is introduced to describe the strain-stress curves of SMA as follows:

where is the stress, is the strain, ( = 1~7) are coefficients, presents the skeleton curve of the hysteretic loop, and is the new nonlinear differential item, which describes the difference between the skeleton curve and the hysteretic loop. The item is developed from the deformed Van der Pol item . ( = 1~7) are coefficients, which can be obtained by fitting method. From (2), we can see that the SMA material has nonlinear damping.

Equation (2) is the new model developed from the model we proposed in [28], and then we should test its accuracy. In this paper, the partial least-square regression software SIMCA-P is used to test the fitting effect. The procedure is as follows:

(1)Firstly, we obtained the stress-strain data of Ti-Ni SMA materials by experiments.(2)Secondly, we choose the constructive model

 Please note that the coefficients are unknown at that time.(3)Thirdly, we input the constructive model and the stress-strain data into the partial least-square regression software SIMCA-P; then the software will calculate the variable importance (VIP) and the coefficient values of each item.(4)When the coefficient values of each item were given, the software can predict the stress-strain data according to the constructive model.(5)Finally, the software can compare the stress-strain data from prediction and that of real stress-strain curves.

The analysis results of the principal component based on the experimental data are shown in Figure 2, and the coefficient values of each item are shown in Figure 3, where VIP is the variable importance. There is no item whose VIP is evidently little, which means that all items should be reserved in the constructive model.

Figure 2: Variable importance of each item.

Figure 3: Coefficient values of each item.

The result of forecast test to (2) is shown in Figure 4, where the red lines present the experimental data and the black lines present the forecast values. Equation (2) can describe the experimental curves well.

Figure 4: Results of the forecast test for the fitting effect of (2) on strain-stress data of SMA.

3. Modelling of SMA Composite Wing

The research object is the Chinese “Wing Loong” UAV, which is shown in Figure 5. We can see that it has high-aspect-ratio wings. To control the wing’s fluttering, we paste SMA thin plate on the inner surface of the wing. The SMA thin plate is made up of several SMA thin strips. The mechanical model of a high-aspect-ratio SMA composite wing is shown in Figure 6. The composite wing can be regarded as a composite cantilever plate, where is the wing’s length; is the wing’s width, is the wing’s thickness; is the relative velocity of airflow; is Gauss white noise. is the in-plane stochastic excitation, which is related with . is standard Gauss white noise whose mean is zero and intensity is 2, . In a high-aspect-ratio wing, there are great bending and torsion deformation. Usually, the dynamic model of a high-aspect-ratio wing can be shown as a two-degree-of-freedom (TDOF) system with two variables: vertical displacement and torsion angle. For a TDOF system subjected to stochastic excitation, it is hard to analyze its dynamic characteristics in theory. In this paper, the UAV is assumed as in level flight, and the torsion angle can be regarded as constant.

Figure 5: Chinese “Wing Loong” UAV.

Figure 6: Mechanical model of a high-aspect-ratio SMA composite wing.

The boundary conditions of the composite wing can be written as follows:

where

Thus, the vibration mode can be assumed as follows:

where

Considering the complex characteristics of composite materials, we introduce Hamilton’s principle to dynamic modeling of the system. The Hamilton function can be presented as

where is kinetic energy of the wing, ; is the kinetic energy of the SMA plate, ; is the potential energy of the wing, ; is the potential energy of the SMA plate, ; is the power of the external force, ; is the density of the materials (1: wing, 2: SMA plate); is the elastic modulus; is the area of cross sections.

According to the Hamilton principle,

Thus, the nonlinear dynamic model of the SMA composite wing can be shown as follows:

According to (6), the dynamic equation of the system’s response can be solved from (10) by the Galerkin method as follows:

When the wing’s length is more than its width, which means , the coupling coefficient is close to zero, and the composite wing can be simplified as a composite beam. In this case, the dynamic equation of system response can be shown as follows:

where is the damping coefficient and

4. Dynamic Characteristics of SMA Composite Wing Subjected to In-Plane Stochastic Excitation

Equation (12) is a stochastic nonlinear differential equation. The deterministic system corresponding to (12) can be shown as follows:

Let and ; (14) can be rewritten as follows:

The nondamping autonomous system from (15) is shown as follows:

Its Hamilton function is

According to the quasi-nonintegrable Hamiltonian system theory, the Hamiltonian function converges weakly in probability to a one-dimensional Ito diffusion process [29, 30]. The averaged Ito equation about the Hamiltonian function can be shown as follows:

where is standard Wiener process and and are the drift and diffusion coefficients of Ito stochastic process, which can be obtained by the stochastic averaging method:

The averaged FPK equation of (18) is

where is the stationary probability density (SPD) function of the system’s response. Thus,

where is a normalization constant.

Equation (21) can also be rewritten as follows:

where , , and . Evidently, decreases monotonically since , , , and . Thus, the bifurcation characteristics of are only determined by the item .

Now, we can calculate the transition sets of the system response. Let and ; we obtain . For , we obtain

Since , we obtain

where , . We have the following:

(1) Bifurcation sets:

 Thus, the bifurcation set of the system is .(2) Hysteretic sets:

 Thus, the hysteretic sets of the system are and , which means .(3) Double limited sets:

 Thus, the double limited sets of the system are and , which means .

The transition sets of the system’s response are shown in Figure 7. These transfer sets divide the parameter space into three parts, and the corresponding stationary probability density functions of the system response in the three parts are also shown in Figure 6, where , , , , , and .

Figure 7: Stationary probability density (SPD) of the system response in different parameters.

From Figure 7, we can see that

(1)when , the steady-state probability density of is the maximum, and there is a crest in the SPD of the system response;(2)when , there is a crest and a loop in the SPD of the system response; although the steady-state probability density of is also the maximum, it decreases with the variation of the parameters;(3)when , there is a loop in the SPD of the system response, and the loop has the maximum SPD; it means that the system’s motion is the stochastic vibration near a periodic orbit;(4)in sum, stochastic Hopf bifurcation appears in the variation of the parameters.

5. Optimal Control to SMA Composite Wing

According to the results of the above section, the stochastic Hopf bifurcation appears in the process of varying the system parameters. Bifurcation changes the system’s motion and decreases its reliability. In this section, stochastic optimal control is introduced to improve the system’s reliability.

The dynamic model of controlled system can be shown as follows:

where is the active control force. To SMA beam, the active control force is supplied by heating and cooling. The Ito diffusion process of (28) can be shown as follows:

where , .

The reliability function of the system is determined by the dynamic programming principle shown as follows:

The optimal control force should make the reliability function the maximum, so it can be determined as follows:

The stress can be obtained when the control force is determined, where is the stress and is the area of the cross section of SMA plate. The vertical displacement can be measured by sensors, and the strain of SMA plate can be calculated. Thus, the input temperature can be determined according to Auricchio’s SMA constitutive model as follows [22]:

where is the elasticity modulus, is the inelastic strain induced by phase change, and is the thermal expansion coefficient. The input current intensity can be determined as follows:

where and are the heat transfer coefficients.

Substituting (31) into (29), we obtained

Thus, the background Kolmogorov equation (BK equation) can be shown as follows:

The initial condition is

The boundary conditions are

The probability density of the first-passage time (FPT) of the controlled system is determined as follows:

The numerical simulation results of the controlled system reliability function are shown in Figure 8, and the probability density of first-passage time is shown in Figure 9, where , , and , where is the mean intensity of control force. In Figure 8, the abscissa is time (), whose unit is second (s); the ordinate is the reliability of the optimal controlled system (), whose unit is 100%. In Figure 9, the abscissa is the first-passage time (), whose unit is second (s); the ordinate is the probability density of the system’s first-passage time (), whose unit is 100%.

Figure 8: Reliability function of the system.

Figure 9: Probability density of first-passage time.

From Figures 8 and 9, we can see the following:

(1)The reliability function decreases as time increases.(2)The system’s reliability is obviously improved when the control force increases, which means that the control strategy is effective.(3)The probability density of first-passage time increases with time. First passage means the system leaves the safe area, which causes the system’s instability.(4)First-passage time can be delayed through optimal control, which also means that the optimal control can enhance the system’s reliability.

6. Experiments

The experimental object is a 1:26 Chinese “Wing Loong” UAV model. The model is shown in Figure 10, and its parameters are shown in Table 1. All the data in Table 1 are measured by us. The model is supplied by Hongda Model Company, Beijing, China. SMA thin plate is pasted on the inner surface of the wing. The SMA thin plate is made up of several SMA thin strips. The SMA thin strips are supplied by Tianxin Company, Gansu Province, China. The SMA’s hysteresis is shown in Figure 1. Its transient temperatures are Ms = 243.4 K, Mf = 223.2 K, As = 304.7 K, and Af = 311.7 K. The system is in a low-speed wind tunnel. The stochastic airflow is produced by adding random guide vanes in the low-speed wind tunnel. The heating method for SMA plate is by heating wire, and the cooling is natural. The vibration amplitude of the wingtip is measured by vibration sensors. The control process is as follows:

(1)Firstly, we measure the vibration of the wingtip by vibration sensors and obtained the displacement and the velocity , and then we can calculate the strain .(2)Secondly, we calculate according to (30) and then obtain through integration.(3)Thirdly, we calculate the optimal control force according to (31).(4)When the optimal control force is given, we can obtain the control stress .(5)When the control stress and the strain are both determined, we can calculate the control temperature according to (32).(6)Finally, we obtained the intensity of control current according to (33).

Table 1: Model’s parameters.

Wing’s mass (g)
12.8
SMA plate’s mass (g)
2.68

Wing’s length (cm)
19.4
SMA plate’s length (cm)
12

Wing’s width (cm)
1.24
SMA plate’s width (cm)
0.7

Wing’s thickness (mm)
1.97
SMA plate’s thickness (mm)
0.4

Figure 10: 1:14 Chinese “Wing Loong” UAV model.

The experimental results of the SMA composite wing subjected to stochastic airflow are shown in Figures 11–20.

Figure 11: Wingtip’s dynamic response of an uncontrolled SMA composite wing when the mean wind velocity m/s.

Figure 12: Frequency domain of the wingtip’s dynamic response of an uncontrolled SMA composite wing when the mean wind velocity m/s.

Figure 13: Wingtip’s dynamic response of a controlled SMA composite wing when the mean wind velocity m/s.

Figure 14: Frequency domain of the wingtip’s dynamic response of a controlled SMA composite wing when the mean wind velocity m/s.

Figure 15: Voltage of the control current when the mean wind velocity m/s.

Figure 16: Wingtip’s dynamic response of an uncontrolled SMA composite wing when the mean wind velocity m/s.

Figure 17: Frequency domain of the wingtip’s dynamic response of an uncontrolled SMA composite wing when the mean wind velocity m/s.

Figure 18: Wingtip’s dynamic response of a controlled SMA composite wing when the mean wind velocity m/s.

Figure 19: Frequency domain of the wingtip’s dynamic response of a controlled SMA composite wing when the mean wind velocity m/s.

Figure 20: Voltage of the control current when the mean wind velocity m/s.

From the experimental results, we can see that

(1)the system’s motion is the mixture of periodic motion and stochastic motion;(2)the system’s vibration increases with the mean wind velocity; when the wind velocity increases from 6 m/s to 8 m/s, the vibration amplitude of wingtip increases from 3.32 mm to 4.21 mm, which rises by 26%;(3)the vibration amplitude of an uncontrolled wingtip is 4.21 mm when the mean wind velocity m/s, and the one of a controlled wingtip is only 3.48 mm; it means that the vibration amplitude of the wingtip is reduced by 17.4% under the control; in fact, the control effect is not very evident; we think the reason is that the time-delayed characteristics of SMA are not considered in the control process, which cause the SMA actuator’s action being not very accurate; if the response time of SMA is considered and the control strategy is developed accordingly, SMA composite beam is a promising way to control the wing’s fluttering.

7. Conclusions

In this paper, a kind of high-aspect-ratio shape memory alloy (SMA) composite wing is proposed to reduce the wing’s fluttering. The nonlinear dynamic characteristics and optimal control of the SMA composite wings subjected to in-plane stochastic excitation are investigated where the great bending under the flight loads is considered. The stochastic stability of the system is analyzed, and the system’s response is obtained. The conditions of stochastic Hopf bifurcation are determined, and the probability density of the first-passage time is obtained. Finally, the optimal control strategy is proposed. Numerical simulation shows that the stability of the system varies with bifurcation parameters, and stochastic Hopf bifurcation occurs in the process; the reliability of the system is improved through optimal control, and the first-passage time is delayed. Finally, the effects of the control strategy are proved by experiments. The results of this paper are helpful for engineering applications of SMA.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the support of the Natural Science Foundation of China (NSFC) through Grants nos. 11272229, 11302144, and 11402168, the Ph.D. Programs Foundation of the Ministry of Education of China through Grant no. 20120032120006, and the Tianjin Research Program of Application Foundation and Advanced Technology through Grants nos. 13JCYBJC17900 and 14JCQNJC05300.

References

S.-M. Yang, J.-H. Han, and I. Lee, “Characteristics of smart composite wing with SMA actuators and optical fiber sensors,” International Journal of Applied Electromagnetics and Mechanics, vol. 23, no. 3, pp. 177–186, 2006.

G. Mirone, “Design and demonstrators testing of adaptive airfoils and hingeless wings actuated by shape memory alloy wires,” Smart Structures and Systems, vol. 3, no. 1, pp. 89–114, 2007.

U. Icardi and L. Ferrero, “Preliminary study of an adaptive wing with shape memory alloy torsion actuators,” Materials and Design, vol. 30, no. 10, pp. 4200–4210, 2009.

U. Icardi and L. Ferrero, “SMA actuated mechanism for an adaptive wing,” Journal of Aerospace Engineering, vol. 24, no. 1, pp. 140–143, 2011.

S. Barbarino, S. Ameduri, L. Lecce, and A. Concilio, “Wing shape control through an SMA-based device,” Journal of Intelligent Material Systems and Structures, vol. 20, no. 3, pp. 283–296, 2009.

T. L. Grigorie and R. M. Botez, “New adaptive controller method for SMA hysteresis modelling of a morphing wing,” Aeronautical Journal, vol. 114, no. 1151, pp. 1–13, 2010.

S. H. Kim and M. H. Cho, “A simple smart wing actuator using Ni-Ti SMA,” Journal of Mechanical Science and Technology, vol. 24, no. 9, pp. 1865–1873, 2010.

M. M. Barzegari, M. Dardel, A. Fathi, and M. Ghadimi, “Aeroelastic characteristics of cantilever wing with embedded shape memory alloys,” Acta Astronautica, vol. 79, no. 10, pp. 189–202, 2012.

W.-R. Kang, E.-H. Kim, M.-S. Jeong, I. Lee, and S.-M. Ahn, “Morphing wing mechanism using an SMA wire actuator,” International Journal of Aeronautical and Space Sciences, vol. 13, no. 1, pp. 58–63, 2012.

J. Colorado, A. Barrientos, and C. Rossi, “Biomechanics of smart wings in a bat robot: morphing wings using SMA actuators,” Bioinspiration and Biomimetics, vol. 7, no. 3, Article ID 036006, 10 pages, 2012.

K.-T. Lau, “Vibration characteristics of SMA composite beams with different boundary conditions,” Materials and Design, vol. 23, no. 8, pp. 741–749, 2002.

K. M. Liew, J. Ren, and S. Kitipornchai, “Analysis of the pseudoelastic behavior of a SMA beam by the element-free Galerkin method,” Engineering Analysis with Boundary Elements, vol. 28, no. 5, pp. 497–507, 2004.

A. Zbiciak, “Dynamic analysis of pseudoelastic SMA beam,” International Journal of Mechanical Sciences, vol. 52, no. 1, pp. 56–64, 2010.

F. Scarpa, M. Ruzzene, and M. R. Hassan, “Spectral element formulation for SMA beams under random vibration excitation,” Smart Materials and Structures, vol. 33, no. 8, pp. 286–293, 2004.

S. M. T. Hashemi and S. E. Khadem, “Modeling and analysis of the vibration behavior of a shape memory alloy beam,” International Journal of Mechanical Sciences, vol. 48, no. 1, pp. 44–52, 2006.

M. Collet, E. Foltête, and C. Lexcellent, “Analysis of the behavior of a shape memory alloy beam under dynamical loading,” European Journal of Mechanics A/Solids, vol. 20, no. 4, pp. 615–630, 2001.

K. Tanaka, “A thermomechanical sketch of shape memory effect: one-dimensional tensile behavior,” Res Mechanica, vol. 18, no. 2, pp. 251–263, 1986.

J. G. Boyd and D. C. Lagoudas, “A thermodynamical constitutive model for shape memory materials. Part I. The monolithic shape memory alloy,” International Journal of Plasticity, vol. 12, no. 6, pp. 805–842, 1996.

L. C. Brinson, “One-dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable,” Journal of Intelligent Material Systems and Structures, vol. 4, no. 2, pp. 229–242, 1993.

E. J. Graesser and F. A. Cozzarelli, “Proposed three-dimensional constitutive model for shape memory alloys,” Journal of Intelligent Material Systems and Structures, vol. 5, no. 1, pp. 78–89, 1994.

Y. Ivshin and T. J. Pence, “Thermomechanical model for a one variant shape memory material,” Journal of Intelligent Material Systems and Structures, vol. 5, no. 4, pp. 455–473, 1994.

F. Auricchio and J. Lubliner, “Uniaxial model for shape-memory alloys,” International Journal of Solids and Structures, vol. 34, no. 27, pp. 3601–3618, 1997.

M. A. Savi, P. M. C. L. Pacheco, and A. M. B. Braga, “Chaos in a shape memory two-bar truss,” International Journal of Non-Linear Mechanics, vol. 37, no. 8, pp. 1387–1392, 2002.

J. W. Sohn, Y. M. Han, S. B. Choi, Y. S. Lee, and M. S. Han, “Vibration and position tracking control of a flexible beam using SMA wire actuators,” Journal of Vibration and Control, vol. 15, no. 2, pp. 263–281, 2009.

H. Li, Z.-Q. Liu, and J.-P. Ou, “Experimental study of a simple reinforced concrete beam temporarily strengthened by SMA wires followed by permanent strengthening with CFRP plates,” Engineering Structures, vol. 30, no. 3, pp. 716–723, 2008.

C. G. Wang and H. F. Tan, “Experimental and numerical studies on wrinkling control of an inflated beam using SMA wires,” Smart Materials and Structures, vol. 19, no. 10, Article ID 105019, pp. 1025–1031, 2010.

M. S. Speicher, R. DesRoches, and R. T. Leon, “Experimental results of a NiTi shape memory alloy (SMA)-based recentering beam-column connection,” Engineering Structures, vol. 33, no. 9, pp. 2448–2457, 2011.

Z.-W. Zhu, W.-Y. Xie, and J. Xu, “Nonlinear dynamic characteristics of SMA simply supported beam in axial stochastic excitation,” Journal of Vibroengineering, vol. 14, no. 3, pp. 1059–1067, 2012.

W. Q. Zhu and Z. L. Huang, “Stochastic stability of quasi-non-integrable-Hamiltonian systems,” Journal of Sound and Vibration, vol. 218, no. 5, pp. 769–789, 1998.

W. Q. Zhu and Z. L. Huang, “Stochastic Hopf bifurcation of quasi-nonintegrable-Hamiltonian systems,” International Journal of Non-Linear Mechanics, vol. 34, no. 3, pp. 437–447, 1999.