Structural Damping Effects on Dynamic Instability of Subtangentially Loaded and Shear Deformable Beck’s Columns

A frequency equation of externally and internally damped and shear-flexible cantilever columns subjected to a subtangentially follower force is analytically derived in a dimensionless form with relation to the linear instability theory of Beck’s columns. Some parametric studies are then performed with variation of two damping coefficients under the assumption of Rayleigh damping. Based on the analysis results, it is demonstrated that three damping cases in association with flutter loads of Beck’s columns can be selected including one case representative of structural damping. Finally, stabilitymaps of shear-flexible and dampedBeck’s columns are constructed for the three damping cases and discussed in the practical range of damping coefficients and shear parameters. In addition, flutter loads and time history analysis results are presented using dimensionless FE analysis and compared with exact solutions.


Introduction
The dynamic instability problem of a cantilever column subjected to a follower force has been well known and various interesting related topics have been intensively studied by many researchers since Beck [1] solved this problem analytically.It is worth referring to the monographs by Ziegler [2], Bolotin [3], and Leipholz [4] who addressed the static and dynamic stability of a nonconservative system.
According to the linear stability theory of nonconservative systems, external damping tends to increase flutter loads but columns under small internal damping lose their stability at the flutter load drastically deceased.In this case, Beck's columns subjected to follower forces slightly larger than internally damped flutter loads become unstable in form of oscillations with a slow growth of amplitude, which is sometimes called a quiet flutter.This destabilizing effect of small internal damping on the stability of nonconservative systems, Ziegler's paradox, has been one of attractive research topics [5][6][7][8][9].Particularly, Langthjem and Sugiyama [10] and Elishakoff [11] published survey papers regarding the dynamic stability of columns subjected to follower loads.
The effects of shear deformation and rotary inertia on the stability of Beck's columns have been investigated by many scholars [12][13][14][15][16][17][18][19].Most have investigated shear effects on the divergence and flutter loads of nonconservative systems using FE formulations, but those studies neglected the effects of damping.Some researchers [20][21][22][23][24][25][26][27] noticed that external and internal (viscoelastic) damping of nonconservative systems can be treated as Rayleigh damping and explored the destabilizing effects of internal damping using an analytical approach or FE analysis.
To the authors' knowledge, it is judged that a closedform solution of Beck's problem considering not only internal and external damping effects but also the effects of shear deformation and rotary inertia has not been proposed to date, and its stability behavior in the practical range of internal and external damping has not been reported.The important points presented in this paper are summarized as follows: (1) A linear stability theory of subtangentially loaded, damped, and shear-flexible Beck's columns is first formulated in a dimensionless form.
(2) A frequency equation and a buckling equation of Timoshenko cantilever beams subjected to both a subtangentially follow force and internal and external damping forces are then derived in the closed form to determine the damped flutter load and divergence load, respectively.
(3) Based on the well-known fact that internal and external damping can be transformed to Rayleigh damping, it is shown that the damping coefficients can be effectively determined using proportional damping and two damping coefficients can be reasonably estimated for real systems through a parametric study of damping coefficients in association with flutter loads of Beck's columns.
(4) Finally, stability maps of shear deformable and damped Beck columns are newly constructed using the analytical solution, which is compared with the FE solution in the practical range of damping coefficients and shear parameters.

Analytical Formulation for Instability Theory of Damped Beck's Columns
In this section, a linear stability theory of damped and shear deformable Beck's columns is formulated using an analytical approach and its frequency equation is derived in the closed form.

Stability Theory of Damped and Shear Deformable Beck's
Columns. Figure 1 shows a prismatic cantilever column subjected to a subtangential follower force  at the tip end  in which the direction of the force changes subtangentially according to the deformed column axis.The coordinate  is measured along the centroidal axis of the column and the cross-section properties are constant throughout the column length.
The extended Hamilton's principle for shear deformable Beck's columns including external and internal damping forces and the nonconservative parameter  can be expressed as where  is the total length of the column; V,  is the lateral displacement of the column and the bending rotation of its cross section; ,   are the sectional area and the effective shear area;  is the mass per unit length; ,   are the flexural and the shear rigidity; (≡ /√/) is slenderness ratio; (≡ /(   2 )) is dimensionless shear parameter;  1 and  2 are the external and internal (viscoelastic) damping coefficients, respectively;  is the nonconservativeness parameter denoting subtangentiality; the shear angle is related to the bending angle by  = V  − .
The first three terms in the square brackets in (1) are the kinetic energy, the elastic strain energy of the system, and the potential energy due to the axial force , respectively; and the fourth, fifth, and last terms denote the works done by nonconservative damping and the follower end force.Also, the symbol  denotes the first variation,  represents time, and  1 and  2 are the integration limits.
For convenience, the following dimensionless variables including  and  are introduced: Equation ( 1) can be then rewritten using (2a), (2b), (2c), (2d), (2e), (2f), and (2g) as Now invoking Hamilton's principle leads to the following dimensionless equation of motion and boundary conditions: in which the geometric and natural boundary conditions are * (0) = 0, (4d) Hamilton's principle of Beck's columns considering shear deformation and rotary inertia based on Engesser's buckling theory has been presented by Attard et al. [19] in which damping effects were neglected.Recently it is worth noting that Elishakoff et al. [28] have discussed Timoshenko beam theory from a historical perspective.

Frequency Equation of Damped and Shear Deformable
Beck's Columns.In order to analytically derive the frequency equation of Beck's columns, the vibrating mode shape may be assumed to be harmonic.Then substituting V * =  Ω * ( * ) and  * =  Ω * ( * ) into (4a) and (4b) and eliminating  lead to where the dimensionless frequency Ω can be a complex number defined as in which  and  are the real and imaginary parts of the frequency.
Finally the general solution to (5) is obtained with consideration of (4a) and (4b) as where In contrast, the boundary conditions of (4c)-(4f) can be rewritten with respect to ,   as follows: Substituting (7a), (7b), (7c), (7d), (7e), and (7f) into (8) and invoking the four boundary conditions, a set of homogeneous equations with respect to , , , and  is obtained.Now for a nontrivial solution to exist, the determinant of the homogeneous equation must be zero.Finally we obtain a frequency equation of shear-flexible and damped Beck's columns subjected to a subtangential follower force in the analytical form as where In particular, if both shear deformation and rotary inertia effects are neglected ( = 0.0,  = ∞), (9) is reduced to In contrast, the frequency is zero in the case of a static divergence system.So the solution to (5) is in which the divergence loads can be evaluated from the following buckling equation: Detinko [29] presented a generalized version of (11) considering additionally the effects of a tip mass but neglected shear effects.However, to the authors' knowledge, the frequency equation ( 9) and the buckling equation ( 13) of shearflexible, subtangentially loaded, and internally and externally damped Beck's columns have first been derived in the closed form.
The critical flutter and divergence loads can be determined by constructing double eigencurves of  * − Ω 2 based on (9).Generally, the stability of Beck's columns depends on the location of Ω on the complex plane.As the follower force  * is increased, the system is stable if Ω stays in the lefthand half-plane ( ≤ 0).In an undamped nonconservative system, Ω remains on the pure imaginary axis ( = 0) at first, but the two frequencies coalesce and the real part  changes negatively around the undamped flutter load  * fl .The flutter loads  * fl  of the damped system are calculated when  transits from negative to positive values with the oscillatory part  being nonzero.On the other hand, static divergence occurs when  becomes positive and  is equal to zero.Therefore, the divergence loads  *  can be directly calculated from the root of (13).
Later it is required to analytically evaluate the natural frequencies of a shear-flexible cantilever beam without both damping forces and follower forces.To this end, (7c) and (7d) are simplified as 2 . ( The resulting frequency equation for free vibration of the shear-flexible cantilever is

FE Formulation of Damped and Shear-Flexible Beck Column
A FE formulation of Beck's column is briefly presented based on the dimensionless energy expression of (3) in this section, and Rayleigh damping is discussed in relation to internal and external damping terms in the next section.
Figure 2 shows Beck's column subjected to a follower force  and a lateral impulse force  * () at the tip end which is modeled using 20 two-node shear-flexible Hermitian beam elements.The element model considered here consists of two nodes, and each node has two degrees of freedom with element length .The nondimensional vertical displacement V * and rotation  * of a typical point within the element can be related to the nodal displacements using the third-order Hermitian interpolation polynomial as follows: where it is noted that the element length  in the lowercase "" in ( 17) is calculated by dividing the dimensionless column length of 1.0 by the element number.Now, substituting the above interpolation functions into (3) and integrating it, the equations of motion are expressed as where the 4 × 4 element mass matrix m e and elastic and geometric stiffness matrices k e , k g are where The load correction stiffness matrix k nc is not null only for the element including the tip end, as shown in (19c), and the damping matrix is discussed in the next section.The matrix formulation presented is similar to that by Attard et al. [19] but was reproduced for the sake of completeness.First, the critical divergence load  *  is evaluated from the following equilibrium equations of Beck's column: In the case of the undamped system, by putting U =  Ω * H, the equations of motion are written as As  * increases, stability is lost when two consecutive eigenvalues of Ω become equal at a finite critical value of The values in parentheses denote the shear parameter and frequencies evaluated for the column given in Table 1.
* .Beyond this value, known as the undamped flutter load  * fl , the perturbed motion of the system displays explosively diverging oscillations with increasing amplitudes, while the critical eigenvalues become complex conjugates.
Finally, the flutter load  * fl  of the damped flutter system is determined by adding the damping forces to (22) and solving the quadratic eigenproblem:

Estimation of Damping Coefficients for Stability Analysis of Beck's Columns
In deriving (18), it should be noted that integrating the damping-related terms in (3) leads to = u e ( * 1 m e +  * 2 k e ) u e (24) which means that the damping matrix in (9) due to external and internal damping terms can be expressed in the form of Rayleigh damping as follows: In other words, the external and internal damping coefficients  * 1 ,  * 2 are directly connected to Rayleigh damping coefficients.Hence, the damping coefficients can be easily determined, under the assumption of being proportional damping, as follows: where  * 1 ,  * 2 denote the first and second dimensionless natural frequencies of a cantilever beam, respectively, and  * 1 ,  * 2 can be obtained by solving simultaneous equation (26).It is worth noting from (26) that  * 1 ,  * 2 under the assumption of Rayleigh damping increase in direct proportion to the damping ratio .Additionally, in the case of isotropic materials with a rectangular section, the following relationship between the slenderness ratio  and the shear parameter  can be obtained: where   = 10(1 + ])/ (12 + 11]), ] = 0.3.Table 1 presents the material and geometric properties of Beck's column model and the corresponding dimensionless data, and Table 2 lists dimensionless natural frequencies and damping coefficients for  chosen to be 0.0, 0.005, and 0.02 in this study.For a column length of 1.0 m, the shear parameter  and the slenderness ratio  of the cantilever column model are evaluated as 2.55 × 10 −5 and 346.4,respectively, in which  is so small that it is expected that the shear effect can be neglected.Table 2 shows that the two natural frequencies calculated for  = 0.0, 2.55 × 10 −5 are almost identical, which means that  and  can be regarded as infinity and zero, respectively.Values in round brackets denote damped flutter loads in the case of  = 0.5.
To evaluate shear and rotary inertia effects on divergence and flutter loads, the three slenderness ratios corresponding to three shear parameters are determined using (27) and displayed in Table 1.In addition, the dimensionless natural frequencies and damping coefficients corresponding to  = 0.001, 0.02 are calculated from ( 15) and ( 26), and the results are presented in Table 2. Conclusively, after determining dimensionless damping coefficients, the dimensionless flutter and divergence loads of Beck's column model can be determined from the closed-form solutions of ( 9) and ( 13) and from the FE solutions of ( 23) and (21).
To explore the external and internal damping effects on the flutter behaviors of Beck's columns effectively, three representative cases of external and internal damping coefficients, namely, the extreme case of external damping only ( * 1 in the case of 2% proportional damping,  * 2 = 0), the extreme case of small internal damping only ( * 1 = 0,  * 2 in the case of 0.1% proportional damping), and the practical case of external and internal damping ( * 1 ,  * 2 in the case of 2% proportional damping), have been selected based on Table 2, as follows: In this, Cases 1, 2, and 3 have been chosen to investigate the effects of external damping only, very small internal damping only, and realistic damping containing both external and internal damping forces, respectively.Here, the damped flutter loads of Beck's columns corresponding to Cases 1, 2, and 3 are denoted by  * fl 1 for  = 0.02,  * fl 2 for  = 0.001, and  * fl 3 for  = 0.02, respectively.In the remaining part of this section, a parametric study is performed to demonstrate that the three damping cases can indeed be regarded as representative values of internal and external damping.
Table 3 shows how the damped flutter loads  * fl  of shearrigid Beck's columns vary for  = 1.0 as the two damping coefficients increase independently, where values in round brackets denote the damped flutter loads in the case of  = 0.5.From Table 2, it is observed that the damped flutter load tends to decease with introduction of internal damping, but external damping alleviates this tendency, which has been well known.In particular, it should be noticed that flutter loads hardly alter in the practical range when the ratio of two damping coefficients  * 1 / * 2 is maintained constant as 10.0 and 100.0 (marked by the superscripts ‡ and † , resp.).
To investigate this tendency further, Table 4 shows fluctuations of  * fl  for  = 0.0, 0.02 with increasing damping ratio , where values in round brackets denote the flutter loads for  = 0.5.It should be noted that the damped flutter loads in the second, third, and fourth columns of Table 4 are evaluated based on the three damping conditions of taking  * 1 with  * 2 being set as zero,  * 2 with  * 1 being set as zero, and both  * 1 and  * 2 from two coefficients  * 1 ,  * 2 calculated for  ranging from 0.01 to 100%.Also note that the three conditions correspond to the  * 1 / * 2 ratio being constant at ∞, 0.0, and 77.5.Interestingly, Table 4 demonstrates that overall the damped flutter loads are hardly sensitive to the proportional damping ratio  but strongly depend on the damping ratio,  * 1 / * 2 , except for the last low corresponding to 100% damping ratio.That is, it is observed from the second column that the flutter loads  * fl 1 with consideration of external damping only remain around the undamped flutter load  * fl , 20.05 (16.05) and 12.12 (10.95) for  = 0.0, 0.02 despite the increased external damping.Also, the flutter loads  * fl 2 of the third column with increasing internal damping only are nearly equal to the damped flutter load, 10.94 (9.870) and 8.083 (7.529) under the extreme condition of very small internal damping only.In particular, the fourth column shows that the damped flutter loads  * fl 3 under both damping coefficients are around 17.03 (14.20) and 10.90 (9.969) at a damping ratio  of up to 10%.Here it should be mentioned that Bolotin and Zhinzher [20] and Herrmann and Jong [21] noticed the importance of the ratio between internal and  Values in round brackets denote the damped flutter loads in the case of  = 0.5.
external damping, and, recently, Kirillov and Seyranian [22] demonstrated this dependence of damped flutter loads on the damping ratio in an analytical form.From these observations, the following summary can be made: (1) Two damping coefficients  * 1 ,  * 2 increase in direct proportion to the damping ratio  under the assumption of Rayleigh damping.

Stability Maps of Damped and Shear-Flexible Beck's Columns
Stability maps of damped and shear-flexible Beck's columns subjected to a subtangential follower force are presented and compared in this section.Stability maps of subtangentially loaded and damped Beck's column for  = 0.0, 0.005, and 0.02 have been constructed from analytical solutions using Mathematica [30] as the subtangentiality  is varied in the range of 0.0 ≤  ≤ 1.0. is divided into five ranges according to the instability characteristics (Table 5), and the resulting stability maps are presented in Figures 3-6.Additionally, Table 6 lists the divergence and flutter loads of Beck's column at increments and important transition points of  related to Figures 3-6.The figures and Table 6 indicate that not only the critical divergence loads but also the undamped and damped flutter loads are greatly reduced as the shear parameter  is increased.Firstly, Figure 3 shows that an undamped Beck's column with  = 0.0, 0.005, and 0.02 may be stable or unstable depending on the range of  as follows: (1) Range 1: it loses its stability by divergence at load levels higher than  * 1 .(2) Range 2: it is first stable with increasing follower force but becomes unstable due to divergence at  * 1 , then due to dynamic flutter at  * fl1 , and again due to divergence at load levels higher than  * fl2 .instability behaviors can be addressed depending on  as follows: (1) Ranges 1 and 2: it is the same as those of the undamped system.
(2) Ranges 3 and 4: it is stable initially, then loses its stability due to divergence at  * 1 , and becomes unstable due to quiet flutter at  * 2 and then due to violent flutter at  * fl1 .Finally it is again governed by divergence at much higher values of  * fl2 .(3) Range 5: it loses its stability due to quiet flutter at a damped flutter load  * fl 2 and due to violent flutter at  * fl1 .The instability behaviors for Case 3 (external and internal damping) are as follows: (1) Ranges 1 and 2: it is the same as those of the undamped system.
(2) Range 3: overall it is the same as those of Range 3 of Case 2.
(  7 shows the analytical solutions with FE solutions using 10, 20, and 40 elements for critical loads of shearflexible Beck's columns for  = 0.4 and 1.0.Clearly, FE solutions for shear deformable Beck's columns converge to the analytical solutions as the total number of elements is increased. Finally, it is meaningful to investigate how not only the dimensionless frequencies Ω (≡  * +  * ) but also dynamic   flutter behaviors vary using a time history analysis by the FE procedure with increasing internal and external damping coefficients.
Table 8 shows the dimensionless frequencies of shearrigid Beck's columns subjected to follower forces 1% more or less than the damped flutter loads, which are  * fl 1 (= 20.05),  * fl 2 (= 10.94), and  * fl 3 (= 17.03) corresponding to three cases of external damping only, internal damping only, and external and internal damping for  = 0.001, 0.02, and 0.08.Also, for time history analysis of the system, Beck's column model is discretized using 20 elements and subjected to a follower force and a small lateral load at the tip end is taken into account (see Figure 2), and a time history analysis is performed.Data of the beam model are shown in Table 1 and the time increment Δ (Δ * ) is fixed to be 0.0001 s (0.000913).Three follower forces, 1.01 * fl 2 under internal damping only, 1.01 * fl 3 under external and internal damping, and 1.01 * fl 1 under external damping only for  = 0.001, 0.02, and 0.08, act with a constant magnitude and the small disturbance  ( * ) of 0.01 N (0.001) vanishes after 0.0001 s elapses, as shown in It is apparent from the sign of the real part of frequencies in Table 8 that the system is stable under a follower force less than the damped flutter load but unstable under a follower force more than the damped load.In addition, it is noted from Table 8 and Figure 8 that internal damping makes the nonconservative system destabilized while external damping has a stabilizing effect on the system.Particularly, it is found that the positive real parts of the frequencies marked as * * in Table 8 become very small in the case of a 0.01% damping ratio.Related to them, the results of dynamic analyses in Figures 8(a) and 8(b) suggest flutter behaviors with very slowly increasing amplitudes, which is why these dynamic phenomena are called quiet flutter.
Lastly, three damped Beck columns subjected to a lateral impulse force are analyzed under the same conditions as in Figure 8, except that the shear parameter  is increased and the damping ratio is fixed at 2%.The results of time

Figure 2 :
Figure 2: FE modeling of Beck's column subjected to a follower force and a lateral impulse load.

( 2 )
The damped flutter loads are closely dependent on the  * 1 / * 2 ratio.(3) The damped flutter loads hardly vary when two coefficients  * 1 ,  * 2 are increased in the practical range of  * 1 / * 2 ratio maintained at a specific value.(4) Cases 1 and 2 represent two extreme cases of external and internal damping only and the flutter loads of shear-rigid Beck's columns for the two cases will be similar to 20.05 (16.05) and 10.94 (9.870), respectively, despite the increase in the corresponding damping coefficients.(5) Particularly, Case 3, which corresponds to a 2% damping ratio under Rayleigh damping, represents the most realistic damping case and the damped flutter loads for Case 3 are unchanged 17.03, (14.20) with increasing damping ratio .

Figure 3 :
Figure 3: Stability regions of shear deformable Beck's column with variation of .

Figure 6 :
Figure 6: Stability regions of Beck's column for  = 0.02 with variation of .

) Range 4 :
it is the same as Range 3 until the follower force reaches  * 2 but shows restabilization between  * 2 and the damped flutter load  * fl 3 and quiet flutter instability between  * fl 3 and  * fl1 .(4) Range 5: it loses its stability due to quiet flutter at  * fl 3 and to violent flutter at  * fl1 .Table

Figure 7 .
Figure 7. Figures 8(a)-8(c) show the results of time history analyses of damped Beck's columns under a lateral impulse force and three follower forces.It is apparent from the sign of the real part of frequencies in Table8that the system is stable under a follower force less than the damped flutter load but unstable under a follower force more than the damped load.In addition, it is noted from Table8and Figure8that internal damping makes the nonconservative system destabilized while external damping has a stabilizing effect on the system.Particularly, it is found that the positive real parts of the frequencies marked as * * in Table8become very small in the case of a 0.01% damping ratio.Related to them, the results of dynamic analyses in Figures8(a) and 8(b) suggest flutter behaviors with very slowly increasing amplitudes, which is why these dynamic phenomena are called quiet flutter.Lastly, three damped Beck columns subjected to a lateral impulse force are analyzed under the same conditions as in Figure8, except that the shear parameter  is increased and the damping ratio is fixed at 2%.The results of time

Table 1 :
Material and geometric properties of the Beck's column model.

Table 2 :
Dimensionless natural frequencies and damping coefficients of a cantilever beam.

Table 3 :
Damped flutter loads of shear-rigid Beck's columns with variation of two damping coefficients for  = 1.0.

Table 5 :
Ranges of  according to instability characteristics.

Table 6 :
Divergence and damped flutter loads of Beck's columns for  = 0.0, 0.005, and 0.02 with variation of .

Table 7 :
Flutter and divergence loads of shear-flexible Beck's columns by analytical and FE solutions for  = 0.4, 1.0.