## Analysis of Vibration of a Pipeline Supported on Elastic Soil Using Differential Transform Method

**Mohammad Amin Rashidifar**^{1,}, **Ali Amin Rashidifar**^{2}

^{1}Department of Mechanical Engineering, Islamic Azad University, Shadegan Branch, Shadegan, Iran

^{2}Department of Computer Scince, Islamic Azad University, Shadegan Branch, Shadegan, Iran

### Abstract

In this paper, a simulation method called the differential transform Method (DTM) is employed to predict the vibration of an Euler-Bernoulli and Timoshenko beam (pipeline) resting on an elastic soil. The differential transform method is introduced briefly. DTM can easily be applied to linear or nonlinear problems and reduces the size of computational work. With this method exact solutions may be obtained without any need of cumbersome work and it is a useful tool for analytical and numerical solutions. To make clear and illustrate the features and capabilities of the presented method, different problems have been solved by using the technique and solutions have been compared with those obtained in the literature.

### At a glance: Figures

**Keywords:** differential transform method, DTM, elastic soil, vibration, beam, pipeline

*American Journal of Mechanical Engineering*, 2013 1 (4),
pp 96-102.

DOI: 10.12691/ajme-1-4-4

Received April 6, 2013; Revised May 14, 2013; Accepted June 18, 2013

**Copyright**© 2014 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Rashidifar, Mohammad Amin, and Ali Amin Rashidifar. "Analysis of Vibration of a Pipeline Supported on Elastic Soil Using Differential Transform Method."
*American Journal of Mechanical Engineering*1.4 (2013): 96-102.

- Rashidifar, M. A. , & Rashidifar, A. A. (2013). Analysis of Vibration of a Pipeline Supported on Elastic Soil Using Differential Transform Method.
*American Journal of Mechanical Engineering*,*1*(4), 96-102.

- Rashidifar, Mohammad Amin, and Ali Amin Rashidifar. "Analysis of Vibration of a Pipeline Supported on Elastic Soil Using Differential Transform Method."
*American Journal of Mechanical Engineering*1, no. 4 (2013): 96-102.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### 1. Introduction

It is obvious that the analysis, design and performance of buried structures such as pipelines require an understanding of soil-structure interaction. Soil protects and supports the buried pipelines, thus reduces the risk of possible hazards that can lead to malfunction of the structure. Without utilizing the strength of the surrounding soil, buried pipelines usually cannot resist the loads and displacements to which they are subjected. Therefore, it is important to accurately evaluate soil restraint or loading on the buried pipelines so as to increase the service life of the structure.

In geotechnical engineering, there are numerous studies in which the structural elements as pipelines, shallow foundations and piles are presented as a beam. On the other hand various types of foundation models such as Winkler, Pasternak, Vlasov, etc. have been used in the analysis of structures on elastic foundations.

The well known and widely used mechanical model is the one devised by Winkler. According to the Winkler model, the beam-supporting soil is modeled as a series of closely spaced, mutually independent, linear elastic vertical springs which provide resistance in direct proportion to the deflection of the beam. In the Winkler model, the properties of the soil are described only by the parameter k, which represents the stiffness of the vertical springs. However, due to its inability to take into account the continuity or cohesion of the soil, the Winkler model is considered as a rather crude approximation of the true mechanical behavior of the soil material. The assumption that there is no interaction between adjacent springs, also results in overlooking the influence of the soil on either side of the beam. To overcome this weakness, several two-parameter elastic foundation models have been suggested. In these models, while the first parameter represents the stiffness of the vertical spring as in the Winkler model, the second parameter is introduced to account for the coupling effect of the linear elastic springs ^{[1]}.

As stated before pipelines, shallow foundations and piles may be modeled as a beam. There are different types of beam model. The well known’s are Euler-Bernoulli and Timoshenko beam models. Euler-Bernoulli theory works well for slender beams. On the other hand this theory is not applicable for moderately short and thick beams. However, Timoshenko model evaluates the effects of transverse shear deformation and rotary inertia on the dynamic behavior of beams.

There are numerous studies on the vibration analysis of a beam resting on an elastic foundation. De Rosa ^{[2]} investigated the free vibration frequencies of Timoshenko beams on two-parameter elastic foundation, in which the differential equation of motion was deduced by following the geometrical approach. Matsunaga ^{[3]}, based on the power series expansions of displacement components, derived a fundamental set of equations of a one-dimensional higher order theory of deep elastic beam-columns resting on elastic foundations by using Hamilton’s principle. El-Mously ^{[4]} derived explicit formulae for the fundamental natural frequencies for vibration of finite Timoshenko-beams mounted on finite Pasternak-foundation by virtue of Rayleigh’s principle. Chen ^{[5]} investigated the vibration of a beam resting on an elastic foundation by using differential quadrature element method (DQEM). Chen ^{[6]} developed the DQEM free vibration analysis model of non-prismatic shear deformable beams resting on the elastic foundation. Coşkun ^{[7]} studied the response of an elastic beam on a two-dimensional tensionless Pasternak foundation that was subjected to a central concentrated harmonic load, in which the governing differential equations of the problem were solved by using the trigonometric-hyperbolic functions and suitable boundary and contiunity conditions. Chen et al. ^{[8]} studied on a mixed method which combined the state space method and the differential quadrature method for bending and free vibration of arbitrarily thick beams resting on a Pasternak elastic foundation. Maheswari et al. ^{[9]} used finite difference method for the solution of governing differential equations of the problem investigating the response of a moving load on an infinite beam resting on a reinforced granular bed on soft soil. Auciello and De Rosa ^{[10]} used the differential quadrature method and the integral variational formulation Rayleigh–Ritz method for the dynamic analysis of a foundation beam on a two parameter elastic soil in the presence of subtangential follower force. Elfelsoufi and Azrar ^{[11]} presented a model for the investigation of buckling, flutter and vibration analyses of beams using the integral equation formulation. And recently Ruta ^{[12]} applied Chebyshev series approximation to solve the problem of nonprismatic Timoshenko beam resting on a two-parameter elastic foundation.

In this study, the natural frequencies of a pipeline represented by a uniform beam resting on a Winkler and Pasternak soil are investigated by the differential transform Method (DTM). The Differential Transform Method is a semi analytical–numerical technique based on the Taylor series expansion method for solving differential equations. It is different from the traditional high order Taylor series method. The Taylor series method is computationally taken long time for large orders. However, with DTM, doing some simple mathematical operations on differential equations a closed form series solution or an approximate solution can be obtained quickly. This method was first proposed by Zhou ^{[13]} in 1986 for solving both linear and nonlinear initial-value problems of electrical circuits. Later, Chen and Ho ^{[14]} developed this method for partial differential equations and Ayaz ^{[15, 16]} studied two and three dimensional differential transform method of solution of the initial value problem for partial differential equations. Arikoğlu and Özkol ^{[17]} extended the differential transform method (DTM) to solve the integro-differential equations. Recently, the second author used the DTM method successfully to handle various kinds of rotating beam problems ^{[18, 19, 20]}.

In order to show effectiveness of the DTM method, two different problems will be considered. The first one is Euler-Bernoulli beam resting on Winkler foundation. The second example is Timoshenko beam resting on Pasternak foundation. In the next section, the governing equations and the associated boundary conditions for the above mentioned problems will be given.

### 2. The Equations of Motion

Consider Euler-Bernoulli beam resting on Winkler foundation (Figure 1). The equation of motion for this problem is given as follows.

(1) |

where k is spring constant; w is deflection (m); ρ is the mass density (kg/m^{3}); A is the cross sectional area (m^{2}); E is the Young’s Modulus (Pa) and I is the area moment of inertia about the neutral axis (m^{4}). Here x is the horizontal space coordinate measured along the length of the beam and t is any particular instant of time.

**Fig**

**ure**

**1.**Geometry of a beam on Winkler foundation

In Figure 2, a Timoshenko beam on Pasternak foundation is shown. The governing equations for this problem are given as a system of differential equations. De Rosa ^{[2]} modeled the Pasternak foundation in two different forms. These forms lead the same differential equation for Euler-Bernoulli beam; on the other hand it is not true for Timoshenko beam which is adequate for stocky beams. Following De Rosa ^{[2]}, we will call these definitions as Model I and Model II also. More information about these definitions can be found in De Rosa ^{[2]}. The first equation is the same for both Model I and Model II.

(2) |

Where G is the shear modulus.

Model I:

(3) |

Where is the constant of proportionality between bending moments and bending rotations.

Model II:

(4) |

Where is the constant of proportionality between bending moments and global rotations.

**Fig**

**ure**

**2.**Geometry of a beam on Pasternak foundation

After completing equations of motion and associated boundary conditions, we will concentrate now about free vibration analysis of the beams resting on elastic foundations.

In order to make free vibration analysis of the Euler-Bernoulli beam on the Winkler foundation, let us assume the solution is in the form of a sinusoidal variation of w(x, t) with circular frequency :

(5) |

Substituting equation (5) into equation (1), equations of motion is expressed as follows:

(6) |

Similarly, for the free vibration analysis of Timoshenko beam resting on Pasternak foundation, the solution is assumed in the form of

(7) |

(8) |

Substituting equations (7) and (8) into equations (2) and (3) for Model I; and equations (4) and for Model II, the new form of equations of motion are expressed as follows:

(9) |

Model I:

(10) |

Model II:

(11) |

### 3. Nondimensionalization

The non-dimensional parameters for the Euler-beam on the Winkler foundation are defined as

(12) |

Using these parameters, the nondimensional form of equation (6 ) can be written as:

(13) |

Following De Rosa ^{[2]}, the nondimensional parameters of Timoshenko beam on Pasternak foundation are given as follows:

(14) |

Using these parameters, nondimensional form of Equations can be written as follows:

(15) |

Model I:

(16) |

Model II:

(17) |

### 4. Differential Transformation Method

The differential transform method (DTM) is a transformation technique based on the Taylor series expansion and is a useful tool to obtain analytical solutions of the differential equations. In this method, certain transformation rules are applied and the governing differential equations and the boundary conditions of the system are transformed into a set of algebraic equations in terms of the differential transforms of the original functions and the solution of these algebraic equations gives the desired solution of the problem. It is different from high-order Taylor series method because Taylor series method requires symbolic computation is expensive for large orders. The DTM is an iterative procedure to obtain analytic Taylor Series solutions of differential equations.

Consider a function which is analytic in a domain D and let represent any point in D. The function is then represented by a power series whose center is located at. The differential transform of the function is given by

(18) |

where is the original function and is the transformed function.

The inverse transformation is defined as

(19) |

Combining Eqs.(18) and (19) gives

(20) |

Considering Eq.(20), once more it is noticed that the concept of differential transform is derived from Taylor series expansion. However, the method does not evaluate the derivatives symbolically.

In actual applications, the function is expressed by a finite series and Eq.(20) can be written as follows

(21) |

which means that

is negligibly small. Here, the value of depends on the convergence rate of the natural frequencies. Theorems that are frequently used in the transformation of the differential equations and the boundary conditions are introduced in Table 1 and Table 2, respectively.

### 5. DTM Formulation and Solution Procedure

First of all we will derive DTM form of Equation (9) which models Euler-beam on the Winkler foundation. Here we quit using the bar symbol on and instead, we use . If Table 1 is referred the following expression can be written easily.

(22) |

If Eq. (21) is arranged, a simple requrrence relation can be obtained as follows:

(23) |

DTM formulation of which deals with Timoshenko beam resting on Pasternak foundation is derived as follows.

(24) |

(25) |

(26) |

(27) |

(28) |

(29) |

Here Equation (27) denotes for Model I and Equation (28) denotes for Model II.

The solution procedure of DTM will be shown for Model I with cantilever condition and the values are set to χ=40, r = 12.5, α = 10, S = 1.

and values for can now be evaluated in terms of , and . These values were achieved by using the Mathematica computer package.

where and are constants. Substituting all and terms into boundary condition expressions, i.e. Eqs. (27) and (28), the following equation is obtained.

(30) |

where are polynomials of corresponding to term.

When Eq.(30) is written in matrix form, we get

(31) |

The eigenvalue equation is obtained from equation (31) as follows

(32) |

Solving equation (31), we get where . Here, is the estimated eigenvalue corresponding to . The value of is obtained by the following equation:

(33) |

where is the tolerance parameter

If equation (37) is satisfied, then we have eigenvalue . In general, are conjugated complex values, and can be written as . Neglecting the small imaginary part , we have the natural frequency. In this study the value of n=50 was enough.

### 6. Numerical Examples

The computer package Mathematical is used to solve recurrence relations with associated boundary conditions. In the first example Euler beam resting on Winkler soil with different boundary conditions are investigated. Results are compared with Chen ^{[5]} study in which he used differential quadrature element method (DQEM). Chen ^{[5]} set all values to unity, such as I=E=A=ρ=1, hence λ=1. In Table 3, first three natural frequencies for simply supported beam are compared with Chen ^{[5]} results and exact solution. Very good agreement is observed with exact solution. The first five natural frequencies for fixed-fixed beam and cantilever beam are presented in Table 4 and Table 5 respectively. Again good a greement with Chen ^{[5]} results are demonstrated.

In the second example, Timoshenko beam resting on Pasternak foundation are investigated and results are compared with the ones De Rosa ^{[2]}. In Figure 3, the first nondimensional frequency μ vs. χ, of a simply supported beam is shown for =100, r=12.5 and S_{w}′= = 1 and 2.5. The difference between Model I and II is more evident for small χ values. For large χ values μ converges to Euler-Bernoulli limiting case. The same curves are given for a cantilever beam in Figure 4, but for different α value, which is now α=10. Again the discrepancies are high for small χ values. In Table 6-9 we presented nondimensional frequencies only limiting values for K_{TR}, such as K_{TR}=0 is cantilever beam case, on the other hand K_{TR}=∞ is clamped-simply supported beam case. Results are given for first five natural frequencies, both for Model I and Model II. We could make comparisons only for first three modes of De Rosa ^{[2]}. Good agreement are shown between the results.

**Fig**

**ure**

**3.**Simply supported beam with α = 100 and r=12.5. First nondimensional frequency as a function of

**Fig**

**ure**

**4.**Cantilever beam with α = 10 and r=12.5. First nondimensional frequency as a function of

### 7. Conclusions

In this study, we mainly concerned application of Differential Transformation Method (DTM) to predict the vibration of an Euler-Bernoulli and Timoshenko beam (pipeline) resting on an elastic soil. The study showed that by using DTM one can easily obtain solution for the problems described above. DTM eliminates cumbersome computational work needed by Taylor series. When the comparisons are made with the studies in literature, a very good agreement was observed.

### References

[1] | Avramidis IE, Morfidis K. 2006. ‘Bending of beams on three-parameter elastic foundation. International Journal of Solids and Structures’. volume 43. Page 357-375. | ||

In article | CrossRef | ||

[2] | De Rosa MA. 1995. ‘Free vibration of Timoshenko beams on two-parameter elastic foundation’. Journal of Computers and Structures.volume 57(1). Page 151-156. | ||

In article | CrossRef | ||

[3] | Matsunaga H. 1999. ‘vibration and bucklig of deep beam-columns on two-parameter elasti foundatins, Journal of Sound and Vibration’.volume 228(2). page 359-376. | ||

In article | CrossRef | ||

[4] | El-Mously M. 1999. ‘Fundamental frequencies of Timoshenko beams mounted on pasternak foundatıon’. Journal of Sound and Vibration. Volume 228(2). Page 452-457. | ||

In article | CrossRef | ||

[5] | Chen CN. 2000. ‘Vibration of prismatic beam on an elastic foundation by the differential quadrature element method’. Journal of Computers and Structures. volume 77. Page 1-9. | ||

In article | CrossRef | ||

[6] | Chen CN. 2002. ‘dqem vibration analyses of non-prismatic shear deformable beams resting on elastic foundations’. Journal of Sound and Vibration. Volume 255(5). Page 989-999. | ||

In article | CrossRef | ||

[7] | Coşkun İ. 2003. ‘The response of a finite beam on a tensionless Pasternak foundation subjected to a harmonic load’. European Journal of Mechanics A/Solids. Volume 22. Page 151-161. | ||

In article | CrossRef | ||

[8] | Chen WQ, Lü CF, Bian ZG. 2004. ‘A mixed method for bending and free vibration of beams resting on a Pasternak elastic foundation’. Applied Mathematical Modelling. volume28. Page 877-890. | ||

In article | CrossRef | ||

[9] | Maheshwari P, Chandra S, Basudhar PK. 2004. ‘Response of beams on a tensionless extensible geosynthetic-reinforced earth bed subjected to moving loads’. Journal of Computers and Geotechnics. volume 31. Page 537-548. | ||

In article | CrossRef | ||

[10] | Auciello NM, De Rosa MA. 2004. ‘Two approaches to the dynamic analysis of foundation beams subjected to subtangential forces’. Journal of Computers and Structures. Volume 82. Page 519-524. | ||

In article | CrossRef | ||

[11] | Elfelsoufi Z, Azrar L. 2005. ‘Buckling, flutter and vibration analyses of beams by integral equation formulations. Computers and Structures’. Volume 83. Page 2632-2649. | ||

In article | CrossRef | ||

[12] | Ruta P. 2006. ‘The application of Chebyshev polynomials to the solution of the nonprismatic Timoshenko beam vibration problem’. Journal of Sound and Vibration. Volume 296. Page 243-263. | ||

In article | CrossRef | ||

[13] | Zhou JK. 1986. ‘Differential Transformation and its Application for Electrical Circuits’. China: Wuhan, Huazhong University Press. | ||

In article | PubMed | ||

[14] | Chen CK, Ho SH. 1999. ‘Solving partial differential equations by two-dimensional differential transform method. Applied Mathematics and Computation’. Volume 106. Page 171-179. | ||

In article | CrossRef | ||

[15] | Ayaz F. 2003. ‘On the two-dimensional differential transform method’. Journal of Applied Mathematics and Computation 2003; 143: 361-374. | ||

In article | CrossRef | ||

[16] | Ayaz F. 2004. ‘Solutions of the system of differential equations by differential transform method’. Journal of Applied Mathematics and Computation. Volume 147. Page 547-567. | ||

In article | CrossRef | ||

[17] | Arıkoğlu A, Özkol İ. 2005. ‘Solution of boundary value problems for integro-differential equations by using differential transforms method’. Journal of Applied Mathematics and Computation. Volume 168. Page 1145-1158. | ||

In article | CrossRef | ||

[18] | Özdemir Ö, Kaya MO. 2006. ‘Flapwise bending vibration analysis of a rotating tapered cantilever Bernoulli–Euler beam by differential transform method’. Journal of Sound and Vibration. Volume 289. Page 413-420. | ||

In article | CrossRef | ||

[19] | Kaya MO. 2006. ‘Free Vibration Analysis of Rotating Timoshenko Beam by Differential Transform Method’. Journal of Aircraft Eng Aerosp Tec. Volume 78(3). Page 194-203. | ||

In article | CrossRef | ||

[20] | Özdemir Ö, Kaya MO. 2006. ‘Flexural Vibration Analysis of Double Tapered Rotating Euler-Bernoulli Beam by Using the Differential Transform Method’. Journal of Meccanica. | ||

In article | |||