Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 493905, 7 pages

http://dx.doi.org/10.1155/2013/493905

Research Article

Free Vibration Analysis of a Thin-Walled Beam with Shear Sensitive Material

K. B. Bozdogan¹ and D. Ozturk²

¹Kirklareli University, 39000 Kirklareli, Turkey

²Ege University, 35000 İzmir, Turkey

Received 25 June 2013; Accepted 27 August 2013

Academic Editor: Usik Lee

Copyright © 2013 K. B. Bozdogan and D. Ozturk. 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.

This paper presents a method for a free vibration analysis of a thin-walled beam of doubly asymmetric cross section filled with shear sensitive material. In the study, first of all, a dynamic transfer matrix method was obtained for planar shear flexure and torsional motion. Then, uncoupled angular frequencies were obtained by using dynamic element transfer matrices and boundary conditions. Coupled frequencies were obtained by the well-known two-dimensional approaches. At the end of the study, a sample taken from the literature was solved, and the results were evaluated in order to test the convenience of the method.

1. Introduction

In the last two decades research on the dynamics of beams has grown enormously. There are numerous studies [1–29] on the bending-torsion coupled beam. In the beams, the elastic center and the center of mass are not coincident, so the translational and torsional modes are inherently coupled as a result of this offset. Rafezy and Howson [24] proposed an exact dynamic stiffness matrix approach for the three-dimensional, bimaterial beam of doubly asymmetric cross-section. The beam comprises a thin-walled outer layer that encloses and works compositely with its shear sensitive core material.

A dynamic transfer matrix method for the free vibration analysis of a thin-walled beam of doubly asymmetric cross-section filled with shear sensitive materials is suggested in this study. The following assumptions are made in this study: the behaviour of the material is linear elastic, small displacement theory is valid, and the dynamic coupling effect of structure caused by the eccentricity between the center of shear rigidity and the flexural rigidity center is ignored in analysis.

2. Analysis

2.1. Physical Model

Figure 1 shows a uniform, three dimensional beam of length . It has a doubly asymmetric cross-section comprising a thin-walled outer layer that encloses shear sensitive material [24]. The outer layer may have the form of an open or closed section that is assumed to provide warping and Saint-Venant rigidity, while the core materials provide Saint-Venant and shear rigidity. These assumptions lead to a model in which a typical cross-section has independent centers of flexure, shear, and mass denoted by , , and , respectively [24]. For convenience, the origin of the coordinate system is located at the centre of flexure. gives the result that the axis of elastic flexure coincides with the -axis of the member. The - and -axes are subsequently aligned with the principle axes of the cross-section. The locations of the points and in the coordinate system Oxy are given by and , respectively. The resulting elastic shear and mass axes then run parallel to the -axis through and respectively. When the elastic axis of the beam does not coincide, the lateral and torsional motion of the beam will always be coupled in one or more planes [24].

(a)

(b)

Figure 1: Typical thin-walled beam [24].

2.2. Element Transfer Matrices for Planar Motion

The governing equations for th element of uncoupled thin-walled shear sensitive beam can be written as

where .

and are the flexural rigidity of the th segment in the and planes, respectively, and and are the Saint-Venant and warping torsion rigidity of the th segment about , where is the warping moment of inertia or warping constant. and are the effective shear rigidities of the core material of the th segment in and directions, respectively, and is the Saint-Venant torsional rigidity of the core material about . are the mass per unit length of the th segment, and is the polar mass radius of gyration of cross section [24].

If a sinusoidal variation of , , and with circular frequency is assumed then

where , , and are the amplitudes of the sinusoidally varying displacement.

Substituting (2) in (1) results are

When (3) is solved with respect to , , , and can be obtained as follows:

where , , , , , and can be calculated as follows:

By using (4), (5), and (6), the rotation angles in and directions , rate of twist , bending moments in and directions and bimoment , shear forces in and directions , and torque for th element can be obtained as follows:

The following equation shows the matrix form of (4), (8), (11), and (14):

For the direction, the following shows the matrix form of (5), (9), (12), and (15):

Similarly, torsional motion can be written:

At the initial point of the th element, (17), (18), and (19) can be written as follows:

When vector is solved from (20) and is substituted in (17), the following is obtained:

For , (23) can be written as

where is the element dynamic transfer matrix of the th element.

For the direction, (23) and (24) can be written as follows:

Similarly, rotation motion can be written in equations as follows:

If (24) is written successively, the displacements—internal forces relationship between the initial part and end of the beam—can be found as follows:

For and rotation motion, (29) can be written as follows:

The eigenvalue equation for a thin-walled beam filled with shear sensitive material can be established using (29), (30), and the specific boundary conditions are as follows.

(1) Clamped-Free: , , .(2) Clamped-Clamped: , , .(3) Simply-Simply:   , .(4) Free-Free: , , .(5) Clamped-Simply: , , .

In frequency equations the values of , which set the determinant to zero, are the uncoupled angular frequencies.

2.3. Coupled Frequencies

Ignoring the dynamic coupling effect of structure caused by the eccentricity between the center of shear rigidity and the geometric center the coupled frequencies of the shear torsional beam can be obtained by using uncoupled frequencies and the well-known equation as follows [28]:

3. Procedure of Computation

A program that considers the method presented in this study as a basis has been prepared in MATLAB, and the operation stages are presented below.

(1)element dynamic Transfer matrices are calculated for each element by using (24), (26), and (28).(2)System dynamic transfer matrices (see (29)–(30)) are obtained with the help of element transfer matrices. (3)The angular frequencies of uncoupled vibrations are obtained by using the boundary conditions. (4)The coupled angular frequencies are found by using (31).

4. A Numerical Example

In this part of the study two numerical examples were solved by a program written in MATLAB to validate the presented method. The results are compared with those given in the literature.

4.1. Numerical Example 1

The first example considers the beam studied by Tanaka and Bercin [11]. A typical uniform thin-walled beam has a length of 1.5 m with a doubly asymmetric cross section. The properties of the cross section are as follows:

 ,  .

The first three coupled natural frequencies of the beam are calculated by the presented method and compared with the results by Tanaka and Bercin [11] and Rafezy and Howson [24] in Table 1 for clamped-free (C-F) and simply-simply (S-S) boundary conditions.

Table 1: Coupled natural frequencies for the beam of example 1.

Natural frequencies (Hz)

BC
Proposed method
Tanaka and Bercin [11]
Rafezy and Howson [24]

C-F
17.17
27.31
59.10
17.03
27.58
59.25
17.17
27.31
59.10

S-S
44.71
75.14
164.87
41.48
74.12
164.11
44.71
75.14
164.87

4.2. Numerical Example 2

A typical continuous beam with a doubly asymmetric cross section is considered in this example (Figure 2).

Figure 2: The doubly asymmetric, continuous channel section and the cross section of beam of example 2 with warping allowed at B, C, and D but fully constrained at A.

The beam comprises a thin-walled outer layer and a shear core with the following properties between support points A and B. The typical uniform thin-walled beam has a length of 1.5 m with a doubly asymmetric cross section. The properties of the cross section are as follows:

 ,  ,  .

The shear core is omitted between points B and D, where the cross-sectional properties remain unchanged, except that , and the small change in has been ignored.

The first three coupled natural frequencies of the beam are calculated by the presented method and compared with the results of Rafezy and Howson [24] in Table 2.

Table 2: Coupled natural frequencies of the continuous beam of example 2.

Frequency number
This study
Rafezy and Howson[24]
Difference (%)

1
6.906
6.940
−0.49

2
19.763
19.796
−0.17

3
35.461
33.836
4.80

The main source of error between the proposed method and Rafezy and Howson methods is the eccentricity between the center of shear stiffness and flexural stiffness which was not taken into account in the proposed method.

5. Conclusions

This paper presents a method for a free vibration analysis of a thin-walled beam of doubly asymmetric cross section filled with shear sensitive material. In the study, first of all, a dynamic transfer matrix method was obtained for planar shear flexure and torsional motion. Then, uncoupled angular frequencies were obtained by using dynamic element transfer matrices and boundary conditions. Coupled frequencies were obtained by the well-known two-dimensional approaches. It was observed from the sample taken from the literature that the presented method gave sufficient results. The error margin of the proposed method is shown to be less than 5%. The main source of error is the eccentricity between the center of shear stiffness and flexural stiffness which was not taken into account in the proposed method.

The transfer matrix method is an efficient and computerized method which also provides a fast and practical solution since the dimension of the matrix for the elements and system never changes. Because of this the proposed method is simple and accurate enough to be used both at the concept design stage and for final analyses.

References

F. Y. Cheng, “Vibrations of timoshenko beams and frameworks,” Journal of Structural Engineering, vol. 96, no. 3, pp. 551–571, 1970.

C. Mei, “Coupled vibrations of thin-walled beams of open section using the finite element method,” International Journal of Mechanical Sciences, vol. 12, no. 10, pp. 883–891, 1970.

W. L. Hallauer and R. Y. L. Liu, “Beam bending-torsion dynamic stiffness method for calculation of exact vibration modes,” Journal of Sound and Vibration, vol. 85, no. 1, pp. 105–113, 1982.

E. Dokumaci, “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, 1987.

R. H. Gutierrez and P. A. A. Laura, “Approximate analysis of coupled flexural-torsional vibrations of a beam of non-uniform cross-section using the optimized rayleigh method,” Journal of Sound and Vibration, vol. 114, no. 2, pp. 393–397, 1987.

J. R. Banerjee, “Coupled bending-torsional dynamic stiffness matrix for beam elements,” International Journal for Numerical Methods in Engineering, vol. 28, no. 6, pp. 1283–1298, 1989.

J. R. Banerjee and F. W. Williams, “Coupled bending-torsional dynamic stiffness matrix for timoshenko beam elements,” Computers and Structures, vol. 42, no. 3, pp. 301–310, 1992.

J. R. Banerjee and F. W. Williams, “An exact dynamic stiffness matrix for coupled extensional-torsional vibration of structural members,” Computers and Structures, vol. 50, no. 2, pp. 161–166, 1994.

X. Chen and K. K. Tamma, “Dynamic response of elastic thin-walled structures influenced by coupling effects,” Computers and Structures, vol. 51, no. 1, pp. 91–105, 1994.

J. R. Banerjee, S. Guo, and W. P. Howson, “Exact dynamic stiffness matrix of a bending-torsion coupled beam including warping,” Computers and Structures, vol. 59, no. 4, pp. 613–621, 1996.

M. Tanaka and A. N. Bercin, “Free vibration solution for uniform beams of nonsymmetrical cross section using Mathematica,” Computers and Structures, vol. 71, no. 1, pp. 1–8, 1999.

S. M. Hashemi and M. J. Richard, “A dynamic finite element (DFE) method for free vibrations of bending-torsion coupled beams,” Aerospace Science and Technology, vol. 4, no. 1, pp. 41–55, 2000.

R. D. Ambrosini, J. D. Riera, and R. F. Danesi, “A modified Vlasov theory for dynamic analysis of thin-walled and variable open section beams,” Engineering Structures, vol. 22, no. 8, pp. 890–900, 2000.

L. P. Kollár, “Flexural-torsional vibration of open section composite beams with shear deformation,” International Journal of Solids and Structures, vol. 38, no. 42-43, pp. 7543–7558, 2001.

Y. Matsui and T. Hayashikawa, “Dynamic stiffness analysis for torsional vibration of continuous beams with thin-walled cross-section,” Journal of Sound and Vibration, vol. 243, no. 2, pp. 301–316, 2001.

V. H. Cortínez and M. T. Piovan, “Vibration and buckling of composite thin-walled beams with shear deformability,” Journal of Sound and Vibration, vol. 258, no. 4, pp. 701–723, 2002.

A. Arpaci, S. E. Bozdag, and E. Sunbuloglu, “Triply coupled vibrations of thin-walled open cross-section beams including rotary inertia effects,” Journal of Sound and Vibration, vol. 260, no. 5, pp. 889–900, 2003.

M. Y. Kim, H. T. Yun, and N. I. Kim, “Exact dynamic and static element stiffness matrices of nonsymmetric thin-walled beam-columns,” Computers and Structures, vol. 81, no. 14, pp. 1425–1448, 2003.

L. Jun, L. Wanyou, S. Rongying, and H. Hongxing, “Coupled bending and torsional vibration of nonsymmetrical axially loaded thin-walled Bernoulli-Euler beams,” Mechanics Research Communications, vol. 31, no. 6, pp. 697–711, 2004.

L. Jun, H. Hongxing, S. Rongying, and J. Xianding, “Dynamic response of axially loaded monosymmetrical thin-walled Bernoulli-Euler beams,” Thin-Walled Structures, vol. 42, no. 12, pp. 1689–1707, 2004.

F. Mohri, L. Azrar, and M. Potier-Ferry, “Vibration analysis of buckled thin-walled beams with open sections,” Journal of Sound and Vibration, vol. 275, no. 1-2, pp. 434–446, 2004.

A. Prokić, “On triply coupled vibrations of thin-walled beams with arbitrary cross-section,” Journal of Sound and Vibration, vol. 279, no. 3–5, pp. 723–737, 2005.

M. O. Kaya and O. O. Ozgumus, “Flexural-torsional-coupled vibration analysis of axially loaded closed-section composite Timoshenko beam by using DTM,” Journal of Sound and Vibration, vol. 306, no. 3–5, pp. 495–506, 2007.

B. Rafezy and W. P. Howson, “Exact dynamic stiffness matrix for a thin-walled beam of doubly asymmetric cross-section filled with shear sensitive material,” International Journal for Numerical Methods in Engineering, vol. 69, no. 13, pp. 2758–2779, 2007.

H. H. Chen and K. M. Hsiao, “Coupled axial-torsional vibration of thin-walled Z-section beam induced by boundary conditions,” Thin-Walled Structures, vol. 45, no. 6, pp. 573–583, 2007.

D. Ambrosini, “On free vibration of nonsymmetrical thin-walled beams,” Thin-Walled Structures, vol. 47, no. 6-7, pp. 629–636, 2009.

G. M. Vörös, “On coupled bending-torsional vibrations of beams with initial loads,” Mechanics Research Communications, vol. 36, no. 5, pp. 603–611, 2009.

B. Rafezy and W. P. Howson, “Exact natural frequencies of a three-dimensional shear-torsion beam with doubly asymmetric cross-section using a two-dimensional approach,” Journal of Sound and Vibration, vol. 295, no. 3–5, pp. 1044–1059, 2006.

F. de Borbón and D. Ambrosini, “On free vibration analysis of thin-walled beams axially loaded,” Thin-Walled Structures, vol. 48, no. 12, pp. 915–920, 2010.