SOLID MECHANICS
Flexural-torsional vibration analysis of axially loaded thin-walled beam
Aleksandar Prokić^{I}; Dragan Lukić^{II}
^{I}Faculty of Civil Engineering, University of Novi Sad, Department of Civil Engineering, 24000 Suborica, Serbia, aprokic@eunet.rs
^{II}Faculty of Civil Engineering, University of Novi Sad, Department of Civil Engineering, 24000 Subotica, Serbia, drlukic.lukic@gmail.com
]]>
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_{D} the 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_{1} = z and z_{2} = z + dz (Fig. 2) subjected to external loads per unit area of midplane is considered.
]]>
At any point on the cross section z_{1} 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 asTo 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 inSetting 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^{2}
E = 68.9×10^{6} kN / m^{2}
G = 26.5×10^{6} kN / m^{2} ]]>
ρ = 2.711 kNs^{2} / m^{4}
I_{xx} = 9.26×10^{-8} m^{4}
I_{yy} = 1.77×10^{-8} m^{4}
I_{ωω} = 1.52×10^{-12} m^{6}
K = 1.64×10^{-9} m^{4}
I_{D} = 1.843×10^{-7} m^{4}
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^{2} ]]>
E = 2.164×10^{8} kN / m^{2}
G = 0.801×10^{8} kN / m^{2}
ρ = 7.8055 kNs^{2} / m^{4}
I_{xx} = 0.450×10^{-6} m^{4}
I_{yy} = 0.940×10^{-7} m^{4}
I_{ωω} = 0.1636×10^{-9} m^{6}
K = 0.140×10^{-9} m^{4}
I_{D} = 9.256×10^{-.7} m^{4}
x_{D} = 0.0
y_{D} = 0.03771 m
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).
]]>References
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. [ Links ]
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. [ Links ]
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. [ Links ]
Banerjee, J.R., Williams, F.W., 1992, "Coupled bending–torsional dynamic stiffness matrix for Timoshenko beam elements", Computers and Structures, Vol. 42, pp. 301-310. [ Links ]
]]>Banerjee, J.R., Williams, F.W., 1994, "Coupled bending–torsional dynamic stiffness matrix of an axially loaded Timoshenko beam element", International Journal of Solids and Structures, Vol. 31, pp. 749-762. [ Links ]
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. [ Links ]
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. [ Links ]
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. [ Links ]
Hashemi, S.M., Richard, M.J., 2000, "Free vibrational analysis of axially loaded bending–torsion coupled beams: a dynamic finite element", Computers and Structures, Vol. 77, No. 7, pp. 711-724. [ Links ]
]]>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. [ Links ]
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. [ Links ]
Li, J., Li, W.Y., Shen, R.Y., Hua, H.X., 2004a, "Coupled bending and torsional vibration of nonsymmetrical axially loaded thin-walled Bernoulli–Euler beams", Mechanics Research Communications, Vol. 31, No. 6, pp. 697-711. [ Links ]
Li, J., Shen, R.Y., Hua, H.X., Jin, X.D., 2004b, "Coupled bending and torsional vibration of axially loaded thin-walled Bernoulli–Euler beams including warping effects", Applied Acoustics, Vol. 65, pp. 153-170. [ Links ]
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. [ Links ]
]]>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. [ Links ]
Prokić, A., 2006, "On fivefold coupled vibrations of Thimoshenko thin-walled beams", Engineering Structures, Vol. 28, pp. 54-62. [ Links ]
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. [ Links ]
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
]]>