Modeling and Dynamical Behavior of Rotating Composite Shafts with SMA Wires

A dynamical model is developed for the rotating composite shaft with shape-memory alloy (SMA) wires embedded in.The rotating shaft is represented as a thin-walled composite of circular cross-section with SMA wires embedded parallel to shaft’s longitudinal axis. A thermomechanical constitutive equation of SMAproposed by Brinson is employed and the recovery stress of the constrained SMA wires is derived. The equations of motion are derived based on the variational-asymptotical method (VAM) and Hamilton’s principle. The partial differential equations of motion are reduced to the ordinary differential equations of motion by using the Galerkin method. The model incorporates the transverse shear, rotary inertia, and anisotropy of composite material. Numerical results of natural frequencies and critical speeds are obtained. It is shown that the natural frequencies of the nonrotating shaft and the critical rotating speed increase as SMA wire fraction and initial strain increase and the increase in natural frequencies becomes more significant as SMA wire fraction increases. The initial strain of SMA wires appears to have marginal effect on dynamical behaviors of the shaft. The actuation performance of SMA wires is found to be closely related to the ply-angle.


Introduction
Composite materials have found the increased applications for replacement of the conventional metallic materials in the rotating flexible shaft employed for drive shafts of helicopters, steam, and gas turbines.This is likely attributed to high stiffness and strength/weight ratios of composite shaft compared with its metallic counterparts.The development trend in design of light-weight composite shafts is towards higher operating speeds, which gives rise to the problems of high vibration amplitude and stability.Seeking the solution of these problems has caused great research effort [1] in the dynamic of composite rotor.
A review on the literature in this area has shown that composite shafts have high whirling resistance capability and are less susceptible to dynamic instability associated with metallic shafts [2].Several attempts to develop mathematical models of spinning composite shafts are reported in the literature.These models include the shaft models based on shell theories [3], or beam theories combined with the strain-displacement relations of the shell theories [4], or a thin-walled beam theory [5].Song et al. [5] developed the composite thin-walled shaft model based on a thinwalled beam theory of Rehfield [6].This model was used to investigate the natural frequencies and stability of the system subject to the variation of the axial edge load and the lamination angle of the composite layer.However, since Rehfield's formulation is known to be nonasymptotically correct [7], there is no guarantee for consistent accuracy on the results.
SMA composites are a new class of materials that have the ability to change both their stiffness and their elastic properties [8].SMA composites consist of SMA actuators embedded into a matrix material or in a fiber-reinforced composite.This stiffness modification occurs as a result of thermally induced martensite phase transformation of SMA actuators embedded in composite structures.Several studies have been done on [9][10][11][12][13][14][15] combining the advantages of both the composite material and the SMA to build smart composite shafts.Baz and Chen [9] have proposed a smart shaft, which can actively stiffen in response to increased rotational speed or increased amplitudes of vibration.Their results showed activating the Nitinol wires results in reducing the amplitude of its vibration amplitudes by about 50%.Gupta [10] investigated the combined effect of embedding the SMA wires in a rotor shaft and change of support stiffness using the SMA on rotor critical speeds.Sawhney and Jain [11] carried out fabrication and experimental investigations on the fiber-reinforced composite shaft embedded with the SMA wires.Baz and Chen [12] studied the static and thermal characteristics of SMA-reinforced composite drive shafts.They found that the Nitinol wires can play a role in enhancing the torsional stiffness.Tylikowski [13] studied the dynamic stability of globally activated simply supported hybrid shells consisting of symmetrical balanced angle-ply laminated classical plies and symmetrically laminated active plies with axially oriented SMA fibers.Their results indicate that the activation significantly increases the critical value of angular velocity.Tylikowski and Hetnarski [14] analyzed the stability criterion of the shaft equilibrium.The shaft is treated as a thin symmetrically laminated shell containing both the conventional fibers and the activated SMA fibers.The results indicated that the SMA activation significantly increases the stability domains of the shaft.Gupta et al. [15] designed an experimental setup to embed prestrained SMA wires in the fiber-reinforced composite shaft.Experimental results showed a noticeable increase in the natural frequency of the composite shaft due to activation of SMA wires.
In our dynamics study, the rotating shaft is represented as a thin-walled composite of circular cross-section embedded with SMA wires.It is the purpose of the present work to study the effect of increase in stiffness and tension in wires due to phase recovery stresses when wires are activated on the rotor dynamic characteristics such as natural frequencies and the critical speeds.The modeling of the rotating shaft is based on the VAM by Berdichevsky et al. [7].The thermomechanical constitutive equation of SMA proposed by Brinson [16] is employed to establish the constitutive equation of the SMA reinforced composite shaft.To determine the rotating shaft's dynamical characteristics, the Galerkin approach is carried out here to approximate the motion equations by a system of ordinary differential equations.Based on these approximate equations the dynamical characteristics of the rotating shaft systems are then calculated.Finally, the influence of SMA activation on dynamical behaviors of the rotating composite shaft has been examined.

Coordinate Systems and Basic
Assumptions.The composite shaft of length , thickness ℎ, and the radius of curvature  rotating with constant rate Ω.As seen in Figure 1, (, , ) is the rotating coordinate system, (, , ) is the inertial coordinate system, and (, , ) is the local coordinate system.The coordinate systems (, , ) and (, , ) have the common origin  located in the geometric center.It is assumed that at  = 0, the axes of the two systems coincide.
In our study, the following assumptions are made.(1) The composite shaft is characterized as a slender thin-walled beam [7]; the geometric dimensions are such that  ≪ , ℎ ≪ , ℎ ≪ ; (2) transverse shear effects are considered, where , ℎ, , and  denote the length, the thickness, the radius of curvature, and the maximum cross-sectional dimension of the cylinder, respectively; (3) the SMA wires are embedded at an interlayer of the shaft and arranged parallel to its longitudinal axis.

Equations of Motion.
Berdichevsky et al. [7] have shown that VAM is an asymptotically correct theory which can be used effectively for the analysis of tubular composite thinwalled.However, VAM does not account for the effects of transverse shear because of which it may generate inaccurate predictions of the rotating composite shaft.In the present work, VAM has been refined to include effects of transverse shear.
The warping functions (, , ) can be assumed as follows: In the above equation, the functions  1 (),  2 (),  3 (), () are associated with physical behavior for the axial strain, the bending curvatures, and the torsion twist rate, respectively.The primes in (2) denote differentiation with respect to .
In ( 1) and ( 2), the expressions of   (, ),   (, ) are in the following form: Based on the displacement representations (1), (2), and (3), and using the linear strain-displacement relations [7], while referring to [17], the strain of the composite shaft can be expressed as The position vector of an arbitrary point on the crosssection of the deformed shaft is r = ( +  2 )i + ( +  3 )j + ( +  1 )k.Here, i, j, k are unit vectors of the coordinate systems (, , ).By taking the time derivatives of unit vectors and adopting the assumption of constant rate , the velocity vector of an arbitrary point is The equations of motion for the rotating composite shaft are derived based on Hamilton's principle, which can be expressed as [18] ∫ where are the strain energy and the kinetic energy of the composite shaft, respectively.
Here  0 and  1 are two arbitrary instants of time,   and   are the stress components and the strain components, respectively,  is the virtual work of the external forces,  is the mass density, and  is the variation operator.
In order to display the elastic coupling between vertical bending and horizontal bending, a special ply-angle distribution referred to as circumferentially uniform stiffness-CUS configuration [19], achieved by skewing angle plies with respect to the longitudinal beam axis meeting the condition () = (−), () = (−), is considered.
By employing Hamilton's principle and taking into account the results with a longitudinal compressive force in [20], which is generated by SMA wires activation, the equations of motion involving CUS configuration in terms of displacements can be written as in which

Shock and Vibration
Herein, where   denotes the local stretching stiffness.
denotes the element of the transformed stiffness matrix of the th layer.
In the case of CUS configuration, it can be seen that ( 7) reveals different elastic couplings which consist of the coupling between flapwise bending and chordwise bending, the coupling between flapwise transverse shear and chordwise bending, and the coupling between chordwise transverse shear and flapwise bending. And in which and Q is the stiffness matrix,  = cos ,  = sin , and  is the fiber ply-angle.And In addition, the mass terms  4 ,  5 , and  6 in ( 7) are expressed as The modulus and Poisson's ratio of the shaft can be determined by considering the mixture rule of composition where the subscript  and  denote the composite matrix and SMA wire, respectively, and  is volume fraction.Herein, where  sma ,  Total express the total cross-sectional areas of the embedded SMA wires and the composite shaft, respectively, and where  and  represent the number and diameter of the SMA wires, respectively.The terms  SMA and  Δ appearing in (7) denote the axial tension due to thermal expansion and SMA wires activation, respectively.
The axial tensions  SMA and  Δ are given by where   ,   , and   denote the reduced thermal expansion coefficients in the th layer of composite medium, and where   represents the recovery stress which can be determined analytically.By eliminating from (7) the quantities  55 (  2 −    ) and  66 (  3 −    ), using the relations   (, ) =   2 (, ) and   (, ) =   3 (, ), the equations of motion for unshearable shaft can be obtained.The results are not presented in this paper for the sake of simplicity.

The Recovery Stress of the Constrained SMA Wires
Based on the one-dimensional model of SMA proposed by Brinson [16] and assuming all SMA wires are fully constrained, the expressions for the recovery stress of SMA wires during heating and cooling can be written as follows, respectively: (i) during heating, (ii) during cooling, where  denotes the martensite fraction, () denotes the elastic modulus of SMA, Θ denotes the thermal elastic modulus,  denotes temperature,  0 is the reference temperature, and Ω() denotes the phase transformation coefficient.The subscript 0 denotes initial state.   and    denote the start and finish temperatures of Austenite in stress, and    and    denote the start and finish temperatures of martensite in stress.  is the martensite fraction induced by stress.The first expression in (21a) is used for SMA in the initial martensite state, while the third expression among them is used for SMA in 100% austenite state, and the second one among them is used for SMA in the phase transformation state from martensite to Austenite.
The transformation processes from martensite to austenite ( → ) and from austenite to martensite are given as follows ( → ), respectively: where   and   denote the start and finish temperature points during the phase transformation from austenite to martensite and   and   denote the start and finish temperature points during the phase transformation from martensite to austenite, respectively.  and   are material constant of SMA wires which determine the influence of the stress on the transformation temperature. cr  and  cr  denote the start and finish stress points of phase transformation, respectively, and the subscript "s0" denotes martensite fraction induced by stress corresponding to initial state.
In (21a) and (21b),  1 and  2 denote the start and finish stress for the transformation ( → );  3 denotes the start stress for the transformation ( → ).
The phase transformation coefficient Ω() and the elastic modulus   () can be expressed, respectively, as where   and   denote Young's modulus in the martensite and austenite phase, respectively, and   is recovery strain limit.

Solution Method of Motion Equations
In order to find the approximate solution of the rotating composite shaft, the quantities  2 (, ),  3 (, ),   (, ), and   (, ) are assumed in the form where Substituting (24) into the governing (7) and applying Galerkin procedure, the following governing equations in matrix form can be found: where (,  = 1, 2, . . ., ) , (,  = 1, 2, . . ., ) . (28) Herein, the generalized displacements can be defined as From the matrix Equation ( 26), the eigenvalue problem of the rotating composite shaft can be expressed as where ) .(31) Lower frequencies are primarily concerned in order to investigate the dynamic behavior of rotating composite shaft system.If there is a value of  for which any increase in  will not significantly change the eigenvalues associated with these lower frequencies, then the eigenvalues are considered to be converged at this value of .In (26) the stiffness matrix is made up of three components: the elastic stiffness of the shaft with stiffness coefficients   (see K 1 in (27)), the geometric stiffness that accounts for the axial tension due to thermal expansion and SMA wires activation, and the stiffness K Ω due to the shaft rotation.It should be remarked that the stiffness matrix of the shaft is reduced by the rotational stiffness matrix K Ω .However, on the other hand, if the total axial tension of the shaft, which is the sum of  SMA and  Δ , is high enough such that ( SMA −  Δ ) is positive then the stiffness of the shaft can been enhanced.Thus, it can be seen that the SMA wires activation plays an important role in controlling stiffness of the shaft and maintaining its stability at higher rotating speed.

Numerical Results and Discussion
The numerical calculations are performed by considering the shaft made of graphite-epoxy whose elastic characteristics are listed in Table 1.The shaft has geometrical characteristics as  = 0.127 m,  = 2.023 m, ℎ = 0.381 mm.The SMA wires used in the numerical simulations are listed in Table 2.
The numerical results will be given in terms of the normalized natural frequencies and rotating rate that are defined by  * = / 0 , Ω * = Ω/ 0 , where the normalizing factor  0 = 138.85rad/s is the fundamental frequency of the nonrotating shaft with  = 0 ∘ and corresponds to the case of the absence of the SMA activation ( SMA = 0) and of thermal effect ( Δ = 0).
The natural frequencies of a cantilever composite shaft obtained for without shear deformation using the present model together with those obtained in [22] are shown in Table 3 for different rotating speeds.A perfect agreement of numerical results with those in [22] can be seen.In order to examine the influence of the number of mode shape functions used in the solution of the equation on the accuracy of the results, the numerical results of natural frequency are shown in Table 4 for an increasing number of mode shape functions.From Table 4, it can be seen that to obtain an accurate result of the first three natural frequencies, no more than six mode shape functions are required.This indicates clearly that the convergence of the present model is quite good.
Based on mode convergence examination, it is found that  = 6 gives suitably converged eigenvalues.So for all results given in this paper,  = 6 unless otherwise noted.
Figure 2 shows the variations of the first three natural frequencies versus rotating speed for without SMA activation and with SMA activation ( = 30 ∘ ,  = 50 ∘ C,   = 0.24,  0 = 0.067).From this figure it clearly appears that when Ω * = 0, a single zero-speed mode natural frequency is obtained.The reason is that for the circular cross-sectional Ply angle () LWF, V s = 0.12 UWF, V s = 0.12 LWF, V s = 0.24  shaft, the frequencies in transversal and lateral bending of each mode coincide.When Ω * ̸ = 0, the natural frequency "splits" into two distinct branches of bending vibration due to the gyroscopic effect, namely, the upper and lower frequency branches.The branch corresponding to the higher frequencies is associated with the up-whirling frequency (UWF) motion.Similarly, the branch corresponding to the lower frequencies is associated with the low-whirling frequency (LWF) motion.The minimum rotating rate at which the LWF becomes zero is referred to as the critical rotating speed that corresponds to the dynamical instability of the rotating shaft.In addition, this figure clearly indicates that SMA actuation can play a significant role in increasing the frequencies of the nonrotating shaft and postponing the occurrence of the whirling instability.
Figure 3 shows the variations of the first three natural frequencies versus rotating speed for different fractions of SMA wires at a temperature above   ( = 30 ∘ ;  = 50 ∘ C;  0 = 0.067).It is obviously seen that the natural frequencies of the nonrotating shaft and the critical rotating speed increase with increasing SMA wire fraction, because increasing of the fraction (or number) of the SMA wires results in generating a higher axial tension  SMA inside the rotating shaft.It thus is apparent that the increase of the SMA wires which induces the increase of bending stiffness brings the increase of overall natural frequencies of the shaft.
Figure 4 shows that the variation of the first three natural frequencies as function of the temperature for different SMA wire fractions (Ω * = 5 rad/s;  = 30 ∘ ;  0 = 0.067).Figure 4 is plotted for both UWF and LWF.The temperature is increased from 0 ∘ C to 100 ∘ C and then decreased from 100 ∘ C to 0 ∘ C. The phase transformation from the martensite to the austenite ( → ) is induced in the SMA wires during heating, and the phase transformation from the austenite to the martensite ( → ) is induced during cooling.From Figure 4, the following phenomena can be observed: (1) the curves of natural frequencies and temperature are typical hysteresis loops in a thermal cycle, and the development of the natural frequencies can be divided into two stage, that is, the increase stage and the decline stage; (2) the natural frequencies increase with the increase in temperature during the phase transformation from the martensite to the austenite but decline with the decrease in temperature during the phase transformation from the austenite to the martensite; (3) the increase of the SMA wire fraction is accompanied by the shift of these hysteresis loops towards higher natural frequencies.
Figure 5 shows the variation of the first three natural frequencies with the ply-angle for different SMA wire fractions (Ω * = 5 rad/s;  = 50 ∘ C;  0 = 0.067).The results show that the effect of SMA wire fraction on the natural frequencies is substantial when the ply-angle is located near 90 ∘ .In addition, from Figure 5 it can be found that the trend of variation of natural frequencies with ply-angle is the same as that of the shaft without SMA actuation.Figure 5 also shows that as the SMA wire fraction increases, the natural frequencies increase.
Figure 6 shows the variations of the first three natural frequencies versus rotating speed for the different initial strains of SMA wire ( = 30 ∘ ,  = 50 ∘ C,   = 0.36).It can be seen that the natural frequencies of the rotating shaft increase with increasing of the initial strain of SMA wire.
Figure 7 shows the variation of the first three natural frequencies and temperature for the different initial strains of SMA wire ( = 30 ∘ ;  = 50 ∘ C;   = 0.36).The results show that there is an obvious change in the shapes of natural frequencies and temperature curves as the initial strain of SMA wire is increased.
Figure 8 shows the variation of the first three natural frequencies with the ply-angle for the different initial strains of SMA wire (Ω * = 5rad/s;  = 50 ∘ C;   = 0.36).It can be noted that the effect of the initial strain of SMA wire on the natural frequencies is similar to that previously presented for the SMA wire fraction.However, the SMA wire fraction has a clear influence on the modes, while the initial strain of SMA wire appears to have a marginal effect on the modes.
Figure 9 shows the variation of the first three natural frequencies with the ply-angle for different rotating speeds (  = 0.24;  = 30 ∘ ;  = 50 ∘ C;  0 = 0.067).The results show that the increase of rotating speed is accompanied by both an upward shift of the UWF branches and by a downward shift of the LWF branches, which are similar to that in connection with the results of Figures 2, 3, and 6.
Figure 10 shows the variation of the first three natural frequencies versus ply-angle during heating and cooling at 50 ∘ C (  = 0.12; Ω * = 5rad/s;  0 = 0.067).For the same temperature, there are two natural frequencies which correspond to the heating and cooling processes, respectively (see Figure 4).Moreover, when rotating speed Ω * ̸ = 0 two different frequencies are produced.Thus, we can obtain totally four different curves at temperature  = 50 ∘ C.  The effect of the fraction and the initial strain of SMA wires on the critical rotating speed during heating at 50 ∘ C is presented in Tables 5 and 6, respectively.
The significant effects of the SMA wire fraction on the critical rotating speed can be observed in Table 5.However, as it appears from Table 6, the initial strain of SMA wires does not influence the critical rotating speed significantly.

Conclusion
A model was presented for the study of the vibration and stability of rotating thin-walled composite shaft embedded with SMA wires.The presented model was used to predict the natural frequencies and dynamical stability.Emphasis is placed on the study of the effects of SMA activation on the above dynamical behavior.From the present analysis and the numerical results, the following main conclusions can be obtained.
(1) The developed model provides means of predicting the natural frequency and critical rotating speed of rotating composite thin-walled shaft with SMA wires actuation.
(2) SMA wires activation can significantly postpone the occurrence of the whirling instability and increase the critical rotating speed through SMA phase transformation.
(3) The SMA wire fraction and ply-angle affect the actuation performance of SMA significantly.There is an obvious increase in the natural frequencies and the critical rotating speed as the SMA wire fraction is increased.
(4) The initial strain of SMA wire appears to have marginal effect on the critical rotating speed of the rotating shaft during heating.

Figure 1 :
Figure 1: Schematic drawing of the composite thin-walled beam of a circular cross-section.

Table 3 :
Comparison of the natural frequencies of a cantilever composite shaft without shear deformation.

Table 6 :
The critical rotating speed corresponding to the different initial strains of SMA wire ( = 50 ∘ C,   = 0.36).