Free vibration and stability analysis of axially functionally graded tapered Timoshenko beams with classical and non-classical boundary conditions

https://doi.org/10.1016/j.compositesb.2011.01.017Get rights and content

Abstract

Free vibration and stability analysis of axially functionally graded tapered Timoshenko beams are studied through a finite element approach. The exact shape functions for uniform homogeneous Timoshenko beam elements are used to formulate the proposed element. The accuracy of the present element is considerably improved by considering the exact variations of cross-sectional profile and mechanical properties in the evaluation of the structural matrices. Carrying out several numerical examples, the convergence of the method is verified and the effects of taper ratio, elastic constraint, attached mass and the material non-homogeneity on the natural frequencies and critical buckling load are investigated.

Introduction

Functionally graded (FG) materials have received great interest from engineers and researchers due to their specific characteristics such as high stiffness and thermal resistance. FG materials are formed by varying percentage content of materials in any desired spatial direction; as a result, the mechanical properties of the new material would show gradation in that direction which could be formulated using different theories among which the most renowned ones are exponential [1] and power laws [2]. Reviewing the literature, it is understood that most of the previous works on FG beams have considered the gradation of material properties in thickwise direction [3], [4], [5] and relatively few researchers [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16] have studied the structural behavior of those types of FG beams whose material properties vary in lengthwise direction. Recent works on axially FG beams have considered Euler–Bernoulli beam theory (EBT). It is well understood that EBT [17], [18], [19], [20], [21] neglects the effects of shear deformation and rotary inertia while Timoshenko beam theory (TBT) takes these effects into account. Therefore, long slender beams could be efficiently modeled by EBT since the flexural behavior is dominant while application of EBT leads to great inaccuracy in modeling of short thick beams, where the shear deformations are of more significance. Moreover, it is well known that TBT presents a more realistic model of the beam in determination of the higher modes of vibration.

Non-prismatic beams, i.e. those beams whose cross-sectional profile changes gradually or abruptly along their length are of great importance in different fields of engineering due to their ability in meeting the architectural and aesthetical needs and optimizing the weight and strength of the structure. This particular advantage of non-prismatic members, i.e. optimization of weight and strength, plays a very important role in construction and performance of aerospace structures. The key problem in analysis of tapered Timoshenko beams is the presence of variable coefficients in the governing differential equations introduced by variable cross-sectional area and moment of inertia. Due to this problem, there are closed-form solutions for neither free vibration nor stability of Timoshenko beams with an arbitrarily variable cross-section. Thus, numerical methods have been used among which finite element method has gained a more prominent position [22], [23], [24], [25]. Analysis of tapered Timoshenko beams with arbitrary distributions of material properties along the beam axis are more complicated than the tapered Timoshenko beams since the variable material properties, in addition to the previously mentioned variable coefficients, show up in the governing differential equations. There are relatively few works [26], [27], [28] in the literature on the axially FG tapered Timoshenko beams. Tong et al. [26] modeled a tapered non-homogeneous Timoshenko beam as an assemblage of several uniform homogeneous segments, but they used their method to study the vibration characteristics of tapered Timoshenko beams with constant material properties. Sorrentino et al. [27] studied the effects of generalized damping on the non-homogeneous Timoshenko beams. Sklyar and Szkibiel [28] studied the controllability of a slowly rotating non-homogeneous Timoshenko beam.

To the authors’ best knowledge, the stability analysis of axially FG tapered Timoshenko beams has not been performed previously and there is still a gap in the literature on the structural analysis of these beams and the effects of material non-homogeneity and taper ratio on their structural behavior. In this paper, finite element method is applied to study the free vibration and stability of these structural members. The most important part of a finite element formulation is the selection of shape functions, i.e. functions which interpolate the displacement field within an element. In this paper, the exact shape functions of a uniform homogeneous Timoshenko beam element are used to formulate the axially FG tapered Timoshenko beam element. This set of shape functions will propose an inaccuracy in the results; however this inaccuracy is lessened by considering the exact variation of cross-sectional profile and the mechanical properties in evaluation of the structural matrices. The present element could be used for modeling Timoshenko beams with any type of cross-sectional variations and gradations of material properties along the beam element; hence it could be used for most of engineering applications dealing with such beams.

The objective of the paper is to introduce a beam element which could be used to fill the present gap in the literature on axially FG tapered Timoshenko beams and present an insight into the structural behavior of these structural members. In the following, the structural model is firstly explained and afterwards the finite element formulations are presented. Finally numerical examples are provided to check the competency of the proposed element.

Section snippets

Structural model

Following the Timoshenko beam theory, the axial and transverse displacement fields are respectively given asU(x,y,z,t)=-zθ(x,t)W(x,y,z,t)=w(x,t)where x,y and z are the spatial coordinates as shown in Fig. 1, t is time, θ and w are respectively the bending rotation and the transverse displacement. Assuming that the beam element is subject to a constant compressive load P, the axial and shear strains are respectively obtained using Eq. (1a), (1b) as,εxx=-zθx+12wx2γxz=wx-θThe stresses could

Finite element formulation

As it was previously mentioned, Eqs. (6a), (6b) and (7) do not have closed-form solutions due to the variable coefficients, thus the exact shape functions could not be determined. In this paper, the exact shape functions of homogeneous and isotropic uniform Timoshenko beams are used. These shape functions are given as [29].

Transverse Shape Functions:ψ1=1μ[μ-12Λη-(3-2η)η2],ψ2=lμ[η(1-η)2+6Λη(1-η)],ψ3=1μ[(3-2η)η2+12Λη],ψ4=-lμ[η2(1-η)+6Λη(1-η)],

Bending Rotation Shape Functions:φ1=6lμ(η-1)η,φ2=1μ(μ-4

Numerical examples

In this section, several numerical examples are provided to demonstrate the competency of the present element. In this paper, the effect of axial compressive load on natural frequencies is not studied; hence in free vibration analysis P is set to zero.

In order to facilitate the presentation of results, the following dimensionless parameters are introduced asr=I0A0L2,ϖ=ωρ0A0L4E0I0,P¯=PcrL2E0I0,Kt=ktL3E0I0,Kr=krLE0I0,α=Imρ0A0L3,β=Mmρ0A0Lwhere L is length of the whole beam, kt and kr are

Conclusions

The free vibration and stability analysis of axially FG tapered Timoshenko beams are carried out through a finite element approach. The structural matrices of the proposed element are evaluated using the shape functions of homogeneous uniform Timoshenko beam elements and the exact variations of the cross-sectional profile and distributions of the mechanical properties. Numerical examples are carried out for a tapered Timoshenko beam which is pure zirconia at the left end and pure aluminum at

Acknowledgment

The financial support provided by the vice chancellor for research of University of Tehran (Grant No. 8102033/1/03) is gratefully acknowledged.

References (32)

Cited by (262)

View all citing articles on Scopus
View full text