Flexural-torsional vibration analysis of axially loaded thin-walled beam

TECHNICAL PAPERS

SOLID MECHANICS

Flexural-torsional vibration analysis of axially loaded thin-walled beam

Aleksandar Prokić^I; Dragan Lukić^II

^IFaculty of Civil Engineering, University of Novi Sad, Department of Civil Engineering, 24000 Suborica, Serbia, [e-mail]

^IIFaculty of Civil Engineering, University of Novi Sad, Department of Civil Engineering, 24000 Subotica, Serbia, [e-mail]

ABSTRACT

The present paper considers the flexure-torsion coupled vibrations of axially loaded thin-walled beams with arbitrary open cross section, by means of an exact solution. The effects of axial force, warping stiffness and rotary inertia are included in the present formulations. In the case of simply supported thin-walled beam, a closed-form solution for the coupled natural frequencies of free harmonic vibrations was derived by using a general solution of the governing differential equations of motion based on Vlasov theory. The method is illustrated by its application to two test examples, to demonstrate the effects of bending-torsion coupling and axial force on the dynamic behavior of thin-walled beams. Compared with those available in the relevant literature, numerical results demonstrate the accuracy and effectiveness of the proposed method.

Keywords: thin-walled beams, flexural-torsional vibration, axially loaded, free vibration

Introduction

Thin-walled beam members are widely used as basic structural elements within the fields of mechanical, civil, aeronautical engineering etc., offering a high performance in terms of minimum weight for a given strength. Because of their practical significance in engineering applications, it is essential for design engineers to evaluate the dynamic characteristics of the thin-walled beam structures accurately and so ensure that their design is reliable and safe. The effects due to axial force on the dynamic response of the thin-walled beams are especially of interest. Helicopter, turbine or propeller blades, plane and space frames and also girders of cable-stayed bridges, all could be qualified as axially loaded structures.

Due to their practical importance mentioned above, the vibration analysis of thin-walled beams have been studied by different authors and numerous approaches for calculating the free vibration frequencies have been proposed (Friberg, 1983; Dokumaci, 1987; Banerjee, 1989 and 1999; Kim et al., 1994; Tanaka and Bercin, 1999; Kollár, 2001; Arpaci and Bozdag, 2002 and 2003; Prokić, 2005 and 2006).

Relatively fewer works are available in literature toward the study of the coupled bending-torsion vibrations of axially loaded thin-walled beams. Banerjee and Williams (1992 and 1994) derived the analytical expressions for the coupled flexural-torsional dynamic stiffness matrix of an axially loaded Timoshenko beam element. But the warping stiffness was not included in their theory. The method is referred as an exact method since it is based on exact shape functions obtained from the exact solution of differential equations. As Moon-Young et al. (2003) pointed out, this analytical method, however, is sometimes inefficient because analytical operations in solving a system of simultaneous ordinary differential equations with many variables may be too complex. Hashemi and Richard (2000) presented a new dynamic finite element for the coupled bending-torsional vibration of axially loaded beams based on the closed-form solutions of the Bernoulli-Euler and St. Venant beam theories. But the warping of the cross-sections, shear deformation and rotary inertia were not included in the formulation. Taking the warping effect into account, Li et al. (2004a) carried out the free vibration analysis of an axially loaded beam with nonsymmetrical open cross section by means of dynamic transfer matrix method. Subsequently, Li et al. (2004b) extended this work to include the effects of rotary inertia and shear deflection, but limited the flexural motion to a single plane. The dynamic transfer matrix method implies mathematical procedures which are sometimes difficult to deal with determining the frequencies values of the complex transcendental characteristic equation.

The investigation in the present paper is partly motivated by the fact that natural frequencies of axially loaded thin-walled beams are often required in the design of many structures. Also, some studies have shown that the effect of axial force on natural frequencies is more pronounced than those of the shear deformation and/or rotary inertia. The proposed method is simple, rapid and accurate enough to be used in preliminary design stage and also for verifying numerical results of complex and time-consuming computer procedures. It is expected that undertaken investigation will be useful for better understanding of dynamic characteristics of thin-walled elements.

Nomenclature

a,b,c,d

= coefficients in the frequency equation

C = origin of coordinate system D = shear center e,s = curvilinear coordinates F

= area of cross section

G = shear modulus h_P = distance from tangent at arbitrary point on contour to pole

I_xx, I_yy, I_ωω, I_D

= geometrical properties of cross section K = Saint Venant torsion constant L = length of element M_x, M_y = bending moments M _ω = bimoment m_x, m_y, m_D, m _ω = external bending moments, torsional moment and bimoment per unit length of the beam N = axial force O = starting point (point from which s is measured) P = external force p = frequency p_x, p_y, p_z

= externally loads per unit length of the beam

= externally applied loads per unit area of midplane of the beam T_D = total torsional moment T_s = Saint Venant torsional moment t = thickness of wall U, V, Φ

= amplitudes of the transverse displacements and torsional rotation

= work of actual stresses u = vector of displacements u, v, w = displacements of shear center

u₊, v₊, w₊

= displacements of an arbitrary point of cross-section V_x, V_y

= shear forces

= work of external load and inertia forces w = axial displacement of cross section considered as rigid x, y, z = descartes coordinates Greek Symbols γ_s

= shear strain

δ = symbol of variation

ε

= strain tensor ε _z = longitudinal strain λ _n

= n π /L

ρ = density σ = stress vector σ _z = normal stress τ

= time

τ_w = shear stress uniformly distributed over thickness τ_zs

= total shear stress

τ_s = Saint Venant shear stress φ = rotation of the cross section around its shear centre ω = warping function Subscripts D = relative to shear centre

Basic Equations

A straight uniform thin-walled beam of length L with nonsymmetrical open cross-section is shown in Fig. 1. The beam consists of a linear elastic material with mass density ρ. The beam is referenced to a right-handed rectangular coordinates system x, y, z, where the z-axis is the initial elastic axis of the beam while x and y are the principal axes of the cross-section. The origin of these axes is located at the centroid C. The shear centre with coordinates x_D and y_D in Cxy is denoted by D. Furthermore, it is assumed that the beam is loaded by a given transverse forces per unit length p_x, p_y and p_z distributed along the centroidal axis, externally applied moment per unit length m_x, m_y and m_D and external distributed bimoment of intensity mω. A constant axial force P is assumed to act through the centroid of the cross-section of the thin-walled beam.

Based on the usual assumptions of Vlasov theory

the cross-section is perfectly rigid in its own plane,

the shear strains in the middle surface of the wall are negligible.

The displacements u₊, v₊ and w₊, of an arbitrary point S_* of the cross section can be described by only four components, three translations u, v and w of pole D and the cross section rotation φ about the same pole:

where x is warping function with respect to pole D.

Component deformations different from zero are given by

where e is the distance of the observed point from the middle surface measured along the normal n .

Reducing the normal stresses at the center of gravity and shear stresses at the pole D, for stress resultants the following expressions are obtained

In Eqs. (3), N represents the axial force, M_x and M_y the bending moments with respect to the x and y axis, V_x and V_y the shear forces in the x and y direction, T_Dthe torsion moment, T_s the Saint Venant torque, M_ω the bimoment and F the area of the cross section.

The equations of motion of thin-walled beam can be obtained using the principle of virtual displacements. All vector and matrix quantities are defined with respect to the right-handed rectangular coordinate system (x, y, z). The z-axis is parallel with the longitudinal centroidal axis of the beam, while x and y are arbitrarily taken.

A small element between cross sections z₁= z and z₂= z + dz (Fig. 2) subjected to external loads per unit area of midplane is considered.

At any point on the cross section z₁ acts as a stress vector

The vector of virtual displacements δu, which satisfies the necessary continuity and displacement boundary conditions, may be adopted in the same form as a vector of real displacements

Virtual displacement parameters, which for distinction from real displacements are marked with prefix δ, are arbitrary functions of coordinates and do not depend upon external loads.

The virtual work expression is

where δ = virtual work of external load and inertia forces through virtual displacements δu and δ = virtual work of actual stresses σ realized through virtual strains δε = [δε_zδγ_T].

The virtual work of the external load and inertia forces, including the second order effects of the constant axial stress , per unit length of the element is

where ρ is the density (mass per unit volume), and ü is the acceleration vector given by

A dot denotes differentiation with respect to time.

Substituting (4), (5) and (8) into (7), the following expression for δW is obtained

The virtual work of the internal load due to the corresponding variation of deformation, per unit length of the element, is

Using expressions (2) for virtual strains, where real displacement should be replaced by virtual displacement, one gets for δ :

By suitable rearrangement of (9) and (11) in accordance with virtual displacement parameters, the principle of virtual work may be expressed as

To satisfy these equations identically for any virtual displacement parameter δw_p, δu_p, δv_p, ,... it is necessary the expressions in the great brackets to vanish. Now, using the expressions for stress resultants (3), one obtains

The forces V_x , V_y and T_D can be eliminated from (13) in order to obtain four equations:

The stress resultants can be expressed directly in terms of the displacements (Prokić, 2005 and 2006). The equations are written in matrix form:

The Equations of motion can be obtained by substituting the stress resultants from (15) into (14)

To achieve the compact form the order of Eq. (14-1) is artificially raised by one. The first equation in (16), describing axial vibration, is uncoupled from the rest of the system and may be analysed independently.

The free harmonic transverse and torsional vibrations are defined by the coupled homogeneous Eqs. (16-2,3,4). The solution may be expressed in the form

where p is the radian frequency and U, V and Φ are amplitudes of the transverse displacements and torsional rotation. Substituting (17) into homogeneous Eqs. (16) yields

In the case of a beam with simply supported ends (fork supports at each end which prevent rotation and can warp freely) the end conditions are

These requirements are satisfied by taking

where C_u, C_v and C_Φ are constants and λ_n = nπ/L n = 1, 2, ......

Substituting (20) into (18) results in

Setting the determinant of the above system equal to zero:

where yields the following algebraic frequency equation of the third order

with the coefficients

Numerical Examples

The method presented in the previous section is used to calculate the natural frequencies of the thin-walled beams. For every mode the three numerical values which characterize three different types of natural frequencies: predominantly torsional, predominantly flexural in x direction and predominantly flexural in y direction are given. The particular examples are chosen to illustrate the effects of axial force on coupled bending-torsion natural vibrations of simply supported thin-walled beams and also to confirm the predictability and accuracy of the theory.

The first example considers a thin-walled beam with semi-circular cross-section, Fig. 3. This example was selected because comparative results are available in the literature, Li et al. (2004a). The geometrical and material properties of thin walled beam are given below.

L = 0.82 m

F = 3.08×10^-4 m²

E = 68.9×10⁶ kN / m²

G = 26.5×10⁶ kN / m²

ρ = 2.711 kNs² / m⁴

I_xx = 9.26×10^-8 m⁴

I_yy= 1.77×10^-8 m⁴

I_ωω = 1.52×10^-12 m⁶

K = 1.64×10^-9 m⁴

I_D = 1.843×10^-7 m⁴

x_D = 0.0

y_D = 0.0155 m

First, the axial force P is simply assumed to be zero and then the effects of a constant compressive axial force (P = 1.79 kN) on the natural frequencies are considered. The natural frequencies obtained by the present method and those found by Li et al. (2004a) are given in Table 1.

The second example examines a thin walled beam with monosymmetical channel cross section as shown in Fig. 4, with geometrical and physical properties listed below.

L = 1.28 m

F = 2.684×10^-4 m²

E = 2.164×10⁸ kN / m²

G = 0.801×10⁸ kN / m²

ρ = 7.8055 kNs² / m⁴

I_xx = 0.450×10^-6 m⁴

I_yy= 0.940×10^-7 m⁴

I_ωω = 0.1636×10^-9 m⁶

K = 0.140×10^-9 m⁴

I_D = 9.256×10^-.7 m⁴

x_D = 0.0

y_D = 0.03771 m

The computational results, both including and excluding the effect of compressive axial force (P = 2.56 kN), are given in Table 2. A very good agreement was found with published results, Li et al. (2004b), obtained by the dynamic transfer matrix method.

Conclusion

Using the principle of virtual displacements the basic governing equations of motion of an axially loaded thin walled beam which exhibits bending-torsion coupling have been derived. The effects of warping and rotatory inertia are also included. By solving the governing differential equations of motion of the beam, the analytical expressions for the coupled bending-torsional vibration of an axially loaded beam are derived in an exact sense. When the results obtained from the present theory are compared with the published results, very good agreement is observed. The method is useful particularly when better accuracy of results or higher frequencies are required.

Acknowledgement

The present work has been supported by Ministry of Science and Environmental Protection of the Republic of Serbia (Project ON174027).

Paper received 2 November 2010

Paper accepted 16 March 2012

Technical Editor: Lavinia Borges

Appendix

The values that determine geometrical properties of cross section are given by

Externally applied loads and moments per unit length of a beam are as follows

• Arpaci, A., Bozdag, E., 2002, "On free analysis of thin-walled beams with nonsymmetrical open cross-sections", Computer and Structures, Vol. 80, pp. 691-695.
• Arpaci, A., Bozdag, S.E., 2003, Sunbuloglu, E., "Triply coupled of thin-walled open cross-section beams including rotary inertia effects", Journal of Sound and Vibration, Vol. 260, pp. 889-900.
• Banerjee, J.R., 1989, "Coupled bending-torsional dynamic stiffness matrix for beam elements", International Journal for Numerical Methods in Engineering, Vol. 28, pp. 1283-1298.
• Banerjee, J.R., Williams, F.W., 1992, "Coupled bendingtorsional dynamic stiffness matrix for Timoshenko beam elements", Computers and Structures, Vol. 42, pp. 301-310.
• Banerjee, J.R., Williams, F.W., 1994, "Coupled bendingtorsional dynamic stiffness matrix of an axially loaded Timoshenko beam element", International Journal of Solids and Structures, Vol. 31, pp. 749-762.
• Banerjee, J.R., 1999, "Explicit frequency equation and mode shapes of a cantilever beam coupled in bending and torsion", Journal of Sound and Vibration, Vol. 224, pp. 267-281.
• Dokumaci, E., 1987, "An exact solution for coupled bending and torsion vibrations of uniform beams having single cross-sectional symmetry", Journal of Sound and Vibration, Vol. 119, No. 3, pp. 443-449.
• Friberg, P.O., 1983, "Coupled vibrations of beams-an exact dynamic element stiffness matrix", International Journal for Numerical Methods in Engineering, Vol. 19, pp. 479-493.
• Hashemi, S.M., Richard, M.J., 2000, "Free vibrational analysis of axially loaded bendingtorsion coupled beams: a dynamic finite element", Computers and Structures, Vol. 77, No. 7, pp. 711-724.
• Kim, M.Y., Chang, S.P., Kim, S.B., 1994, "Spatial stability and free vibration of shear flexible thin-walled elastic beams. I: analytical approach", International Journal for Numerical Methods in Engineering, Vol. 37, pp. 4097-4115.
• Kollár, L.P., 2001, "Flexural-torsional vibration of open section composite beams with shear deformation", International Journal of Solids and Structures, Vol. 38, pp. 7543-7558.
• Li, J., Li, W.Y., Shen, R.Y., Hua, H.X., 2004a, "Coupled bending and torsional vibration of nonsymmetrical axially loaded thin-walled BernoulliEuler beams", Mechanics Research Communications, Vol. 31, No. 6, pp. 697-711.
• Li, J., Shen, R.Y., Hua, H.X., Jin, X.D., 2004b, "Coupled bending and torsional vibration of axially loaded thin-walled BernoulliEuler beams including warping effects", Applied Acoustics, Vol. 65, pp. 153-170.
• Moon-Young, K., Nam II, K., Hee-Taek, Y., 2003, "Exact dynamic and static stiffness matrices of shear deformable thin-walled beam-columns", Journal of Sound and Vibration, Vol. 267, No. 1, pp. 29-55.
• Prokić, A., 2005, "On triply coupled vibrations of thin-walled beams with arbitrary cross-section", Journal of Sound and Vibration, Vol. 279, pp. 723-737.
• Prokić, A., 2006, "On fivefold coupled vibrations of Thimoshenko thin-walled beams", Engineering Structures, Vol. 28, pp. 54-62.
• Tanaka, M., Bercin, A.N., 1999, "Free vibration for uniform beams of nonsymmetrical cross-section using Mathematica", Computer and Structures, Vol. 71, pp. 1-8.

Publication Dates

• Publication in this collection
  01 Nov 2012
• Date of issue
  Sept 2012

History

• Received
  02 Nov 2010
• Accepted
  16 Mar 2012