Abstract

This paper presents an exact approach to investigate the flexural free vibrations of multistep nonuniform beams. Firstly, one-step beam with moment of inertia and mass per unit length varying as and was studied. By using appropriate transformations, the differential equation for flexural free vibration of one-step beam with variable cross section is reduced to a four-order differential equation with constant coefficients. According to different types of roots for the characteristic equation of four-order differential equation with constant coefficients, two kinds of modal shape functions are obtained, and the general solutions for flexural free vibration of one-step beam with variable cross section are presented. An exact approach to solve the natural frequencies and modal shapes of multistep beam with variable cross section is presented by using transfer matrix method, the exact general solutions of one-step beam, and iterative method. Numerical examples reveal that the calculated frequencies and modal shapes are in good agreement with the finite element method (FEM), which demonstrates the solutions of present method are exact ones.

1. Introduction

Beams with nonuniform cross section are widely used in various engineering fields, such as bridges, tall buildings, and helicopter rotor blades. A large number of studies can be found in literature about the free vibrations of nonuniform beams. Vibration problems of beams with nonuniform cross section are often described by partial differential equations and in most cases it is extremely difficult to find their closed form solutions. Consequently, a wide range of approximate and numerical solutions such as Rayleigh-Ritz, Galerkin, finite difference, finite element, and spectral finite element methods have been used to obtain the natural vibration characteristics of variable-section beams [14].

Huang and Li investigated the free vibration of axially functionally graded beams with variable flexural rigidity and mass density [5]. A novel and simple approach was presented to solve modal shapes and corresponding frequencies through transforming traditional fourth-order governing differential equation into Fredholm integral equation. Koplow et al. proposed an analytical solution for the dynamic response of Euler-Bernoulli beams with step changes in cross section, which was verified by experimental tests and receptance coupling methods [6]. Firouz-Abadi et al. studied the transverse free vibrations of a class of variable-cross section beams using Wentzel-Kramers-Brillouin (WKB) approximation [7]. The governing equation of motion for the Euler-Bernoulli beam including axial force distribution was utilized to obtain a singular differential equation in terms of the natural frequency of vibration and a WKB expansion series was applied to find the solution. Inaudi and Matusevich investigated longitudinal vibration problems of variable-cross section rods using an improved power series method [8]. This method introduced domain partition implementation in matrix formulation, as an alternative to other power series techniques in vibration analysis. Therefore, the method solved linear differential equations efficiently up to a desired degree of accuracy and remedies two limitations of the conventional power series method. The Adomian decomposition method (ADM) is employed to investigate the free vibrations of tapered Euler-Bernoulli beams with a continuously exponential variation of width and a constant thickness along the length under various boundary conditions [9]. Duan and Wang demonstrated the free vibration of beams with multiple step changes using the modified discrete singular convolution (DSC) [10]. The jump conditions at the steps were used to overcome the difficulty in using ordinary DSC for dealing with ill-posed problems. A transfer matrix method and the Frobenius method were adopted by J. W. Lee and J. Y. Lee to solve the free vibration characteristics of a tapered Bernoulli-Euler beam and obtain the power series solution for bending vibrations [11].

Nevertheless, besides all advantages of such numerical methods, exact solutions can provide adequate insight into the physics of the problems and convenience for parametric studies. The other advantage of exact solutions is their significance in the field of inverse problems. An exact solution can be more useful than numerical solutions to design the characteristics and damage identification of a structure.

Wang derived the closed form solutions for free vibration of a flexural bar with variably distributed stiffness but uniform mass [12]. Using a systematic approach, Abrate obtained a closed form solution of longitudinal vibration for rods whose cross section varies as [13]. Kumar and Sujith found the exact solutions for longitudinal vibration of nonuniform rods whose cross section varies as and [14]. Li presented an exact approach for free longitudinal vibrations of one-step nonuniform rod with classical and nonclassical boundary conditions [15]. The approach assumed that the distribution of mass is arbitrary, and distribution of longitudinal stiffness is expressed as a functional relation with mass distribution and vice versa. Li et al. obtained exact solutions of flexural vibration for beam-like structures whose moment of inertia and mass per unit length vary as , and , , respectively [16, 17].

The components whose moment of inertia and mass per unit length satisfy and are widely used in civil engineering, such as the bridge with cross section height varying as . This kind of beam-like structure is usually solved by dividing it into several segments whose and satisfy different distributions. The paper derived the exact general solution of one-step nonuniform beam firstly and then obtained the exact solution of multistep nonuniform beam by combining the transfer matrix method and exact solutions of one-step beam. The exact solutions of present method not only can provide convenience for parametric studies, but also are very useful for inverse problems such as damage identification.

2.1. Beams with Uniform Cross Section

The governing differential equation for undamped free flexural vibration of beam with uniform cross section can be written as [18]where is Young’s modulus, is moment of inertia, is mass per unit length, and is transverse displacement at position and time .

Assuming the beam performs a harmonic free vibration at equilibrium position, that is,where is mode shape function of uniform beam, is corresponding natural frequency.

Inserting (2) into (1) obtains

Equation (3) leads to the modal shape function of beam with uniform cross section:where , , , , are () are integration constants, .

2.2. Beams with Variable Cross Section

The governing differential equation for undamped free flexural vibration of beam with variable cross section can be written as [19]where is Young’s modulus, is bending moment of inertia at position , is mass per unit length at position , and is transverse displacement at position and time .

Assuming the beam performs a harmonic free vibration at equilibrium position, that is,where is mode shape function of beam, is natural frequency.

Substituting (6) into (5) arrives at

The moment of inertia and mass per unit length of beam are assumed to vary as where , , and are arbitrary constants but not equal to 0 and is a positive integer.

Substituting (8) into (7) yields

LetIt can be derived that

Introducing (11) into (9), one arrives at

Equation (12) can be simplified as

Let

Equation (13) can be written as

The characteristic equation of (15) is

The four roots of (16) arewhere

The modal shape function of one-step beam with variable cross section can be written as where , , , and are four undetermined coefficients.

Since , then , and , it can be found that and are two real roots in the first equation of (17). From the second equation in (17), it can be found that when , and are also two real roots and when , and are two conjugate complex roots.

Modal shapes of one-step beam are determined by the roots as follows:

(1) When , , , and are all real roots, substituting into (19), one obtains

(2) When and are real roots, while and are two conjugate complex roots, letthe roots and can be expressed aswhere .

Substituting and (22) into (19) yieldsIn (20) and (23), () have the same meanings as those in (4) and they are integration constants. The solutions of (20) and (23) are the exact general solutions only at the condition of . This is because will be equal to zero in (10) if for . And the derivation of constant term for the second equation in (10) has no mathematical meaning. Meanwhile, from the numerical point of view, small value of () is adopted and the uniform beam is approached; then, it could be solved by (20) or (23). However, compared to solutions of uniform beam by (4), the solutions of approximate uniform beam are not exact general solutions but numerical solutions.

Equations (20) and (23) can be expressed with a uniform style as when the modal shape function is (20), for , and when it is (23), for , and , .

3. Transfer Relationship for Undetermined Coefficients of Multistep Beam Modal Shapes

As shown in Figure 1, the multistep beam is divided into segments. For an arbitrary beam segment () with length , a local Cartesian coordinated system is established with the origin locating at the left end of segment, is bending moment of inertia for the th segment, and is mode shape function for the th segment, where is defined as belonging to the th segment in the local Cartesian coordinated system.

Taking the th segment, for example, the continuity of deformations, equilibrium of moments, and shear forces are satisfied at the right end of the th segment and left end of the ()th segment; that is,where prime denotes the derivative with respect to local coordinate .

Substituting modal shape functions of the th and ()th segments ((4) or (24)) into (25) yieldswhere denotes the th undetermined coefficient of modal shape function of the th segment and is the th term in modal shape function of the th segment.

Equation (26) can be rewritten into a matrix form:wherein which

Letthen,

Equation (31) represents the transfer relationship of four undetermined coefficients for modal shapes of the th and ()th segments. is the transfer matrix. In accordance with (31), transfer relationship between the first segment and the last segment can be expressed aswhere .

4. Solutions of Natural Frequencies and Modal Shapes

Based on the modal shape function of one-step beam ((4) or (24)) and transfer relationships for undetermined coefficients of modal shapes for each segment ((32)), natural frequencies and modal shapes of multistep beam can be determined by boundary conditions of the beam. The cantilever and simply supported boundaries are considered in this paper, respectively.

For cantilever multistep beam, the boundary conditions are

Inserting the modal shape functions into (33) and utilizing transfer relationship in (32), one arrives atwhere

For simply supported multistep beam, the boundary conditions are

Similar to (34), one obtainswhere

Equations (34) and (37) can be written into a uniform stylewhere the elements () in are functions of natural frequencies for multistep beam.

Existence of nontrivial solutions of (39) leads to the following frequency equation

If the frequency region [, ] for the th natural frequency of multistep beam and the expressions of modal shape functions for each segment ((20) or (23)) can be determined, then the expressions of elements in matrix can be uniquely determined. The th natural frequency of multistep beam can be calculated by (40) directly, and undetermined coefficients in modal shape function of the first segment can be solved by (39); then, the undetermined coefficients for other segments can be obtained by using transfer relationship in (32).

If the frequency region or expressions of modal shape function for each segment ((20) or (23)) cannot be determined, the natural frequencies and modal shapes can be solved by the half-interval method [20]. The detailed process for determining natural frequencies is as follows. Firstly, an initial value of the th natural frequency is assumed, and the expressions for modal shape function of each segment corresponding to can be determined ((20) or (23)); then, the coefficient matrix corresponding to is obtained and the determinant of matrix is calculated (denoted as ). Then, a new value with (e.g., ) representing the increment of is assumed, and the same calculations are repeated to determine the new determinant corresponding to (denoted as ). If and have opposite signs, there is at least one natural frequency in the interval [, ]; otherwise, let and continue the above procedures until and have opposite signs. For the next step, let ; if and have the same sign, let ; otherwise, let ; a new interval [, ] is obtained. The same calculations are repeated to determine the new interval [, ] until the rank of matrix is equal to three. The accurate values of are obtained, respectively, using the half-interval method.

5. Numerical Examples

To verify the correctness of present analytical method, several numerical examples are studied in this paper. The materials used in these examples are all the same; that is, Yong’s modulus is 3.25 × 1010 Pa and density is 2500 kg/m3. Poisson’s ratio is in a two-dimensional FEM (2D-FEM) analysis.

5.1. Reliability of the Proposed Method
5.1.1. Validation by Finite Element Method

In order to illustrate the proposed method, a two-step beam with 10 m span, as shown in Figure 2, is used as numerical model and solved by one-dimensional FEM (1D-FEM), 2D-FEM, and the present analytical method, respectively.

Cross section of this beam is rectangular with constant width 0.3 m. The heights of sections at , , and are expressed by , , and , respectively. , . And the height of each step varies nonlinearly as to : to :

For 1D-FEM. The finite element method with two-nodes beam element. Each node of the beam element has two degree-of-freedoms, that is, transverse displacement and rotation angle. The cubic Hermite polynomial is adopted as shape function for beam element. Therefore, the shape functions () can be expressed asIn (43), is the length of beam element and = 0.5 m.

Let ; stiffness matrix and mass matrix of the beam element can be evaluated aswhere is the strain matrix of beam element. Using above stiffness matrix and mass matrix, finite element analysis is carried on by the codes written in MATLAB.

For 2D-FEM. A finite element analysis is carried out using ANSYS software. A 2D-FEM model of two-step beam is built with 8-node PLANE183 elements as shown in Figure 3.

The ratios of span and height at () vary from 10 to 100. , , and are variables with , and and are constants. Relative errors of the first two-order natural frequencies obtained by 1D-FEM and 2D-FEM with respect to present analytical method are shown in Figure 4.

As can be seen from Figure 4, the relative errors between 1D-FEM and present method remain at a very small lever all the time no matter how changes, and the relative errors between 2D-FEM and the present method decrease rapidly with increasing. This is because 1D-FEM and the present method are employed without considering the effect of shear deformation and the stiffness and mass matrixes are formed by integration method for the variable moment of inertia and mass in 1D-FEM modeling. However, the shear deformation is inevitable for a two-dimensional FEM model. When reaches a certain lever, the error between 2D-FEM and the present method becomes smaller. For example, when = 60, the frequencies obtained by 2D-FEM and the present method are very close and relative error of the first two order frequencies is about −0.1%. With regard to beam-like structural free vibration, it can be treated as plane stress problem, so a 2D-FEM analysis is more accurate than a 1D-FEM analysis in theory. For easy of calculation, the beam-like structure can be treated as Euler-Bernoulli beam when reaches a certain lever and the effect of shear deformation could be ignored. And the object of this paper is not considering the effect of shear deformation, so the solutions of the proposed method are exact ones based on Euler-Bernoulli beam theory. This leads to different relative errors of 1D-FEM and 2D-FEM shown in Figure 4. It also reveals that the present method possesses favorable accuracy for the larger ratio of span and sectional height.

5.1.2. Numerical Simulation for Small Value of

From the numerical point of view, small value of () is adopted, and the uniform beam is approached. In order to verify the explanations, a one-step beam with small value of (shown in Figure 5) and a one-step uniform beam (shown in Figure 6) are used for numerical simulation by the present method ((20) and (23) for nonuniform beams, (4) for uniform beams).

One-Step Beam with Small Value of . The span is 15 m. Cross section of this beam is rectangular with constant width 0.3 m.

The heights of sections at and are = 0.25 m and = 0.249 m, respectively, and the height varies nonlinearly from to asTherefore,

One-Step Uniform Beam. The beam is also 15 m span with constant width 0.3 m, and the height is 0.249 m.

The one-step beams with different boundary conditions have been analyzed.

Case 1. A cantilever beam model.

Case 2. A simply supported beam model.

The first three natural frequencies of one-step beam with small value of () calculated by (20) or (23) are compared with those of one-step uniform beam calculated by (4); the results are listed in Tables 1 and 2.

As can be seen from the results, the relative errors between nonuniform beam with small value of and uniform beam are very small. It reveals that the uniform beam is approached when is small enough (). However, for small value of representing a uniform beam, the solutions are not exact general solutions but numerical solutions.

5.2. One-Step Beam

A one-step beam with 15 m span, as shown in Figure 5, is used for numerical simulation. Cross section of this beam is rectangular with constant width 0.3 m. The heights of sections at and are = 0.25 m and = 0.175 m, respectively, and the sectional height varies nonlinearly from to asTherefore,

The one-step beam with different boundary conditions has been analyzed.

Case 3. A cantilever beam model.

Case 4. A simply supported beam model.

The first three natural frequencies of one-step beam with different boundary conditions calculated by the proposed method and FEM are listed in Tables 3 and 4; modal shapes are shown in Figures 7 and 8.

As can be seen from the results calculated by the proposed method and FEM, relative errors for the first three natural frequencies of the one-step beam with different boundary conditions are very small. It reveals that the solutions for one-step beam with variable cross section by the proposed method are exact ones. The reasons lie in that only the differential equations of motion for one-step beam with variable cross section are used to make the solutions, and no other assumptions are introduced.

5.3. Multistep Beam

To verify the correctness of present method for free vibration analysis of multistep beams, a two-step simply supported beam with variable cross section (as shown in Figure 9) is used for numerical simulation firstly. The cross section of this beam is rectangular. The segments and have variable heights and widths. In the two-dimensional FEM analysis, a fine mesh of elements is used and the average thickness of element is adopted for element thickness.

For segment , the heights of sections and are = 0.36 m and = 0.24 m, respectively; the widths of sections and are = 0.3 m and = 0.2449 m, respectively. The height and width from to vary asTherefore,

For segment , the heights of sections and are = 0.24 m and = 0.38 m, respectively; the widths of sections and are = 0.2449 m and = 0.3082 m, respectively. The height and width from to vary asTherefore,

The first three natural frequencies of two-step simply supported beam are calculated by the proposed method and FEM, which are listed in Table 5. The modal shapes are shown in Figure 10.

Secondly, a three-step simply supported beam with variable cross section, as shown in Figure 11, is used for numerical simulation. Cross section of this beam is rectangular with constant width 0.3 m. The segments - and - are variable sections, and - is uniform cross section.

For segment , the heights of sections and are = 0.3 m and = 0.2 m, respectively. The height of section from to varies asTherefore,

For segment , the heights of sections and are = 0.2 m and = 0.3 m. The height of section from to varies asTherefore,

The first three natural frequencies of three-step simply supported beam calculated by the proposed method and FEM are listed in Table 6, and the modal shapes are shown in Figure 12.

From Tables 5 and 6 and Figures 10 and 12, for two-step and three-step simply supported beams, it can be found that the first three natural frequencies and modal shapes calculated by the proposed method are in good agreement with those calculated by FEM. It reveals that the solutions for multistep beam with variable cross section by the proposed method are exact ones. This is because only transfer matrix method and exact solutions of one-step beam with variable cross section are used to solve the equations of vibration for multistep beam with variable cross section, and no other assumptions are introduced.

6. Conclusions

In this paper, an exact approach to investigate the flexural free vibrations of multistep nonuniform beams is presented. The differential equation for flexural free vibration of one-step beam with and is reduced to a four-order differential equation with constant coefficients by using appropriate transformations. And the general solutions of one-step beam are obtained, whose modal shape functions are found to have two kinds of expressions. Then, combining transfer method, iterative method, and the general solutions of one-step beam, the solving method for natural frequencies and modal shapes of multistep beam is formulated. Numerical simulations on multistep nonuniform beam are used to verify the feasibility. The following conclusions can be obtained:(1)The comparison of relative errors among the present analytical method and finite element methods reveals that the solutions of proposed method are exact ones based on Euler-Bernoulli beam theory. The present method is suitable for the larger ratio of span and sectional height and possesses favorable accuracy.(2)The relative errors between nonuniform beam with small value of and uniform beam are very small. It reveals that for small value of representing a uniform beam, the solutions are not exact general solutions but numerical solutions.(3)Numerical studies of one-step cantilever beam and simply supported beam indicate that natural frequencies and modal shapes calculated by the proposed method are very close to the FEM results, which demonstrates the solutions of one-step beam are exact ones. And numerical examples of two-step and three-step simply supported beams show that the calculated frequencies and modal shapes are also in good agreements with FEM results, which also demonstrate the solutions of presented method for multistep beam with variable cross section are exact ones. The reasons lie in that only transfer matrix method and general solution of one-step beam are used to obtain the solution equation of natural vibration characteristics, and no other assumptions are introduced.

The exact solutions of proposed method can provide adequate insight into the physics of problems and convenience for parametric studies. Furthermore, they are very useful for inverse problems such as structural damage identification. Therefore, it is always desirable to obtain exact solutions for such problems. However, when the ratio of span and sectional height is small, relative errors between FEM and present method may be rather larger. The deficiency of this paper is not considering shear deformation in vibration analysis of nonuniform beams, leading to larger error for the small ratio of span and sectional height. This also has become the topic which the author further deliberated from now on.

Competing Interests

The authors declare that they have no competing interests regarding the publication of this paper.

Acknowledgments

The authors express their appreciation for financial support of the National Natural Science Foundation of China under Grants nos. 51478203 and 51408258; China Postdoctoral Science Foundation funded project (nos. 2014M560237 and 2015T80305); Fundamental Research Funds for the Central Universities and Science (JCKYQKJC06), and Technology Development Program of Jilin Province.