Skip to main content

Analytic and numeric analysis for deformation of non-prismatic beams resting on elastic foundations

Abstract

Background

The buckling load as well as the natural frequency under axial load for non-prismatic beam is a changeling problem. Determination of buckling load, natural frequency, and elastic deflection is very important in civil applications. The current paper used both perturbation method (PM), analytic method, and differential quadrature method (DQM), numerical method, to find buckling load and natural frequency with different end supports. The deflection of the beam resting on an elastic foundation under transverse distributed and axial loads is also obtained. Both PM and DQM are used for non-prismatic beams with rectangular and circular cross sections in the vibration analysis. The comparisons of results obtained from both PM and DQM showed perfect agreement with analytical solution for uniform beams with different end supports. The PM and DQM succeeded powerfully for investigating the buckling load as well as the natural frequency for non-prismatic beam.

Results

The percentage of relative error between DQM and PM doesn’t exceed than 5% if the gradient of rectangular section height and the gradient of circular section radius are less than 0.6. As the gradient of height and radius increase, the maximum deflection decreases and the location of maximum deflection displaced toward the smaller moment of inertia.

Conclusions

The PM has not been used for solving the problem of non-prismatic beams resting on elastic foundations subjected to transverse distributed and axial loads. The current research proved the good ability of PM as an analytical solution for a complicated problem and defined its range of accuracy as compared to DQM. Also, it introduced accurate empirical formulae to find both natural frequency and buckling load of non-prismatic beams. These empirical formulae represent a good achievement in vibration analysis of non-prismatic beams.

1 Background

The DQM is a numerical technique for initial and/or boundary value problems. It was developed by the late Richard Bellman and his associates in the early 70 s and, since then, the method has been successfully employed in a variety of problems in physical and engineering sciences. The DQM has been projected by its proponents as a potential alternative to the conventional numerical solution techniques such as the finite difference and finite element methods [1]. Non-prismatic beams resting on elastic foundations are important structural elements. The dynamic characteristics of such non-prismatic beams are of considerable importance in many designs. The vibration of non-prismatic beams is formulated as a fourth-order differential equation of variable coefficients; it is possible to determine their exact solutions [2, 3], using homotopy analysis [2] or the Green's functions method [3]. However, in most cases, the solution is obtained by approximate methods such as finite difference [4], DQM [5], boundary element method [6], the power series [7, 8], Lagrange multiplier formulation [9], coupled displacement field method [10]. Sato [11] reported the transverse vibration of linearly tapered beams using Ritz method. He studied the effect of end restraints and axial load on the natural frequencies.

The PM is a method which relies on there being a dimensionless parameter in the problem that is relatively small: ε <  < 1. The PM is used to analyze the nonlinear vibration behavior of imperfect general structures under static load. The effects on the linearized and nonlinear vibrations caused by geometric imperfections, a static fundamental state, and a nontrivial static state are included in the perturbation procedure. The theory is applied in the nonlinear vibration analysis of cylindrical shells [12]. However, the PM has not been used for solving the problem of non-uniform beams resting on elastic foundations under the action of transverse distributed load and axial load. In this paper, the PM is used to obtain the non-dimensional frequencies and buckling loads for a non-prismatic beam resting on elastic foundations. Also, the deflected shape can be obtained by PM.

In the present paper, the concepts of both DQM and PM are presented for the vibration analysis of non-prismatic Euler–Bernoulli beam. The DQM was widely used in vibration analysis of non-prismatic beams. In this paper, we used it again to solve the same problem to investigate accuracy of PM in analysis of vibration of non-prismatic beams. The generated PM results for non-prismatic beam also verified with Dua [13]. From the comparison between PM and DQM, an empirical formula is investigated to find natural frequency and buckling load of non-prismatic beams. These empirical formulae represent a good achievement in vibration analysis of non-prismatic beams. A full parametric study is then performed.

Figure 1 shows the schematic diagram for non-prismatic beam on elastic foundation subjected to axial load “n” and transverse uniform load \(p\left( x \right)\). The free vibration of a non-uniform beam is governed by the following equation

$$\frac{{\partial^{2} }}{{\partial x^{2} }}\left( {{\text{EI}}\left( x \right)\frac{{\partial^{2} Y}}{{\partial x^{2} }}} \right) - \rho A\left( x \right)\omega^{2} Y + n\frac{{\partial^{2} Y}}{{\partial x^{2} }} + kY = p\left( x \right)$$
(1)
Fig. 1
figure 1

Non-prismatic beam on an elastic foundation subjected to axial and transverse loads—rectangular and circular cross section

Using the non-dimensional parameters: \(W = \frac{Y}{L}\), \(X = \frac{x}{L}\), \(Q = \frac{{nL^{2} }}{{{\text{EI}}_{0} }}\), \(K = \frac{{kL^{4} }}{{{\text{EI}}_{0} }}\), \(P = \frac{{p(x)L^{2} }}{{{\text{EI}}_{0} }}\),\(\Omega^{4} = \frac{{\rho A(x)L^{4} \omega^{2} }}{{EI_{0} }}\), and \(S(X) = \frac{{{\text{EI}}(x)}}{{{\text{EI}}_{0} }}\).

Equation (1) by direct substitution will be reduced to

$$S\left( X \right) \cdot \frac{{{\text{d}}^{4} W}}{{{\text{d}}X^{4} }} + 2\frac{{{\text{d}}S\left( X \right)}}{{{\text{d}}X}} \cdot \frac{{{\text{d}}^{3} W}}{{{\text{d}}X^{3} }} + \frac{{{\text{d}}S^{2} \left( X \right)}}{{{\text{d}}X^{2} }} \cdot \frac{{{\text{d}}^{2} W}}{{{\text{d}}X^{2} }} + Q\frac{{{\text{d}}^{2} W\left( X \right)}}{{{\text{d}}X^{2} }} + \left( {K -\Omega ^{4} } \right) \cdot W = P$$
(2)

where \(W\) is non-dimensional transverse deflection, \(X\) is non-dimensional position along beam length, \(Q\) is non-dimensional axial load parameter, \(K\) is linear non-dimensional stiffness of elastic foundation parameter, \(P\) is non-dimensional transverse distributed load, \(\Omega\) is non-dimensional frequency parameter, \(I_{0}\) is the base mass moment of inertia, and S(X) is non-dimensional stiffness parameter of the beam.

In the current paper, the beam has a rectangular cross section of linear variable height and constant width. The height is \(h_{{\text{o}}}\) at the base and \(h_{1}\) at the end. The equations of non-dimensional height, non-dimensional cross section, and non-dimensional stiffness parameters are:

$$\begin{aligned} & H\left( X \right) = \left( {1 + \gamma X} \right) \\ & A\left( X \right) = A_{{\text{o}}} \left( {1 + \gamma X} \right) \\ & S\left( X \right) = \left( {1 + \gamma X} \right)^{3} \\ & {\text{Where}}:\gamma = \frac{{h_{1} - h_{{\text{o}}} }}{{h_{{\text{o}}} }} \\ \end{aligned}$$
(3)

Another case is considered for a beam of circular cross section of linear variable radius. The equations of non-dimensional radius, non-dimensional cross section, and non-dimensional stiffness parameters are:

$$\begin{aligned} & R\left( X \right) = \left( {1 + \gamma X} \right) \\ & A\left( X \right) = A_{{\text{o}}} \left( {1 + \gamma X} \right)^{2} \\ & S\left( X \right) = \left( {1 + \gamma X} \right)^{4} \\ & {\text{Where}}:\gamma = \frac{{R_{1} - R_{{\text{o}}} }}{{R_{{\text{o}}} }} \\ \end{aligned}$$
(4)

2 Methods

2.1 Perturbation method (PM)

The PM is employed for rectangular cross section as:

$$\begin{aligned} S\left( X \right) & = 1 + 3\gamma X + 3\gamma^{2} X^{2} + \gamma^{3} X^{3} \\ \frac{{{\text{d}}S}}{{{\text{d}}X}} & = 3\gamma + 6\gamma^{2} X + 3\gamma^{3} X^{2} \\ \frac{{{\text{d}}^{2} S}}{{{\text{d}}X^{2} }} & = 6\gamma^{2} + 6\gamma^{3} X \\ \end{aligned}$$
(5)
$$\begin{aligned} S\left( X \right) & = 1 + 4\gamma X + 6\gamma^{2} X^{2} + 2\gamma^{3} X^{3} + \gamma^{4} X^{4} \\ \frac{{{\text{d}}S}}{{{\text{d}}X}} & = 4\gamma + 12\gamma^{2} X + 6\gamma^{3} X^{2} + 4\gamma^{4} X^{3} \\ \frac{{{\text{d}}^{2} S}}{{{\text{d}}X^{2} }} & = 12\gamma^{2} + 12\gamma^{3} X + 12\gamma^{4} X^{2} \\ \end{aligned}$$
(6)

By substituting Eqs. (5, 6) in Eq. (2) and considering \(\gamma\) as perturbation parameter \(\varepsilon\)

$$\begin{aligned} & \left( {1 + 3\varepsilon X + 3\varepsilon^{2} X^{2} + \varepsilon^{3} X^{3} } \right)\frac{{{\text{d}}^{4} }}{{{\text{d}}X^{4} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots } \right) \\ & \quad \quad + 2\left( {3\varepsilon X + 6\varepsilon^{2} X + 3\varepsilon^{3} X^{2} } \right)\frac{{{\text{d}}^{3} }}{{{\text{d}}X^{3} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots } \right) \\ & \quad \quad + \left( {6\varepsilon^{2} + 6\varepsilon^{3} X} \right)\frac{{{\text{d}}^{2} }}{{{\text{d}}X^{2} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots } \right) \\ & \quad \quad + Q\frac{{{\text{d}}^{2} }}{{{\text{d}}X^{2} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots } \right) \\ & \quad \quad + \left( {K -\Omega ^{2} \left( {1 + \varepsilon X} \right)} \right)\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots } \right) \\ & \quad = \left( {P_{{\text{o}}} + \varepsilon P_{1} + \varepsilon^{2} P_{2} + \varepsilon^{3} P_{3} + \cdots } \right) \\ \end{aligned}$$
(7)
$$\begin{aligned} & \left( {1 + 4\varepsilon X + 6\varepsilon^{2} X^{2} + 2\varepsilon^{3} X^{3} + \varepsilon^{4} X^{4} } \right)\frac{{{\text{d}}^{4} }}{{{\text{d}}X^{4} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \varepsilon^{4} W_{4} + \cdots } \right) \\ & \quad \quad + 2\left( {4\varepsilon + 12\varepsilon^{2} X + 6\varepsilon^{3} X^{2} + 4\varepsilon^{4} X^{3} } \right)\frac{{{\text{d}}^{3} }}{{{\text{d}}X^{3} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \varepsilon^{4} W_{4} + \cdots } \right) \\ & \quad \quad + \left( {12\varepsilon^{2} + 12\varepsilon^{3} X + 12\varepsilon^{4} X^{2} } \right)\frac{{{\text{d}}^{2} }}{{{\text{d}}X^{2} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \varepsilon^{4} W_{4} + \cdots } \right) \\ & \quad \quad + Q\frac{{{\text{d}}^{2} }}{{{\text{d}}X^{2} }}\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \varepsilon^{4} W_{4} + \cdots } \right) \\ & \quad \quad + \left( {K - \Omega^{2} \left( {1 + \varepsilon X} \right)} \right)\left( {W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots } \right) \\ & \quad = \left( {P_{{\text{o}}} + \varepsilon P_{1} + \varepsilon^{2} P_{2} + \varepsilon^{3} P_{3} + \varepsilon^{4} P_{4} + \cdots } \right) \\ \end{aligned}$$
(8)

where \(0 \le \varepsilon \le 1\).

The solution will be presented for non-uniform rectangular cross section; for zero-order perturbation,

$$\frac{{{\text{d}}^{4} W_{{\text{o}}} }}{{{\text{d}}X^{4} }} + Q\frac{{{\text{d}}^{2} W_{{\text{o}}} }}{{{\text{d}}X^{2} }} + \left( {K -\Omega ^{4} } \right)W_{{\text{o}}} = P_{{\text{o}}}$$
(9)

For first-order perturbation,

$$\begin{aligned} & \frac{{{\text{d}}^{4} W_{1} }}{{{\text{d}}X^{4} }} + Q\frac{{{\text{d}}^{2} W_{1} }}{{{\text{d}}X^{2} }} + \left( {K - \Omega ^{4} } \right)W_{1} \\ & \quad = - 3X\frac{{{\text{d}}^{4} W_{{\text{o}}} }}{{{\text{d}}X^{4} }} - 6\frac{{{\text{d}}^{3} W_{{\text{o}}} }}{{{\text{d}}X^{3} }} + \Omega ^{4} XW_{{\text{o}}} \\ \end{aligned}$$
(10)

For second-order perturbation,

$$\begin{aligned} & \frac{{{\text{d}}^{4} W_{2} }}{{{\text{d}}X^{4} }} + Q\frac{{{\text{d}}^{2} W_{2} }}{{{\text{d}}X^{2} }} + \left( {K -\Omega ^{4} } \right)W_{2} \\ & \quad = - 6\frac{{{\text{d}}^{2} W_{{\text{o}}} }}{{{\text{d}}X^{2} }} - 12X\frac{{{\text{d}}^{3} W_{{\text{o}}} }}{{{\text{d}}X^{3} }} - 3X^{2} \frac{{{\text{d}}^{4} W_{{\text{o}}} }}{{{\text{d}}X^{4} }} \\ & \quad \quad +\Omega ^{4} XW_{1} - 6\frac{{{\text{d}}^{3} W_{1} }}{{{\text{d}}X^{3} }} - 3X\frac{{{\text{d}}^{4} W_{1} }}{{{\text{d}}X^{4} }} \\ \end{aligned}$$
(11)

For third-order perturbation,

$$\begin{aligned} & \frac{{{\text{d}}^{4} W_{3} }}{{{\text{d}}X^{4} }} + Q\frac{{{\text{d}}^{2} W_{3} }}{{{\text{d}}X^{2} }} + \left( {K -\Omega ^{4} } \right)W_{3} \\ & \quad = - 6X\frac{{{\text{d}}^{2} W_{{\text{o}}} }}{{{\text{d}}X^{2} }} - 6X^{2} \frac{{{\text{d}}^{3} W_{{\text{o}}} }}{{{\text{d}}X^{3} }} - X^{3} \frac{{{\text{d}}^{4} W_{{\text{o}}} }}{{{\text{d}}X^{4} }} \\ & \quad \quad - 6\frac{{{\text{d}}^{2} W_{1} }}{{{\text{d}}X^{2} }} - 12X\frac{{{\text{d}}^{3} W_{1} }}{{{\text{d}}X^{3} }} - 3X^{2} \frac{{{\text{d}}^{4} W_{1} }}{{{\text{d}}X^{4} }} \\ & \quad \quad +\Omega ^{4} XW_{2} - 6\frac{{{\text{d}}^{3} W_{2} }}{{{\text{d}}X^{3} }} - 3X\frac{{{\text{d}}^{4} W_{2} }}{{{\text{d}}X^{4} }} \\ \end{aligned}$$
(12)

The left-hand side of Eqs. (8–12) is the same ordinary differential equation and its complementary solution is:

$$\begin{aligned} W_{{i_{C} }} & = C_{{i1}} \cosh \left( {m_{1} X} \right) + C_{{i2}} \sinh \left( { - m_{1} X} \right) \\ & \quad + C_{{i3}} \cosh \left( {m_{2} X} \right) + C_{{i4}} \sinh \left( { - m_{2} X} \right) \\ \end{aligned}$$
(13)

where \(C_{i1} ,C_{i2} ,C_{i3} ,C_{i4}\) are constants, \(i = 0,1,2,3\) and

$$\begin{aligned} m_{1} & = \sqrt {\frac{{ - Q + \sqrt {Q^{2} - 4\left( {K -\Omega ^{4} } \right)} }}{2}} \\ m_{2} & = \sqrt {\frac{{ - Q - \sqrt {Q^{2} - 4\left( {K -\Omega ^{4} } \right)} }}{2}} \\ \end{aligned}$$
(14)

Assuming that the transverse load is a sine load across the beam length, the particular solutions for zero and first order are:

$$W_{{{\text{o}}_{{\text{P}}} }} = \frac{{P_{{\text{o}}} \sin \left( {\pi X} \right)}}{{\pi^{4} - Q\pi^{2} + \left( {K -\Omega ^{4} } \right)}}$$
(15)
$$\begin{aligned} W_{{1_{{\text{P}}} }} & = \phi _{{13}} {\text{e}}^{{m_{1} X}} + \phi _{{14}} {\text{e}}^{{ - m_{1} X}} + \phi _{{15}} {\text{e}}^{{m_{2} X}} + \phi _{{16}} {\text{e}}^{{ - m_{2} X}} \\ & \quad + \phi _{{17}} \cos (\pi x) + \phi _{{18}} x\sin (\pi x) \\ \end{aligned}$$
(16)

where

$$\begin{aligned} & \phi_{13} = \left( {c_{01} + c_{11} f_{1} (x)} \right),\phi_{14} = \left( {c_{02} + c_{12} f_{2} (x)} \right),\phi_{15} = \left( {c_{03} + c_{13} f_{3} (x)} \right) \\ & \phi_{16} = \left( {c_{04} + c_{14} f_{4} (x)} \right),\phi_{17} = f_{5} ,\phi_{18} = f_{6} \\ & f_{1} (x) = A\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)\frac{{x^{2} }}{2} - 6m_{1}^{3} x} \right] \\ & \quad + \left[ {\frac{ - A}{{2m_{1} }} + \frac{{2m_{2} C}}{{m_{1}^{2} - m_{2}^{2} }}} \right]\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)x - 6m_{1}^{3} } \right] + \left[ {\frac{A}{{4m_{1}^{2} }} - \frac{{4m_{1} m_{2} C}}{{^{{\left( {m_{1}^{2} - m_{2}^{2} } \right)^{2} }} }}} \right]\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)} \right] \\ & f_{2} (x) = - A\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)\frac{{x^{2} }}{2} - 6m_{1}^{3} x} \right] \\ & \quad + \left[ {\frac{ - A}{{2m_{1} }} + \frac{{2m_{2} C}}{{m_{1}^{2} - m_{2}^{2} }}} \right]\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)x - 6m_{1}^{3} } \right] + \left[ {\frac{ - A}{{4m_{1}^{2} }} - \frac{{4m_{1} m_{2} C}}{{^{{\left( {m_{1}^{2} - m_{2}^{2} } \right)^{2} }} }}} \right]\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)} \right] \\ & f_{3} (x) = C\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)\frac{{x^{2} }}{2} - 6m_{1}^{3} x} \right] \\ & \quad + \left[ {\frac{ - C}{{2m_{2} }} + \frac{{ - 2m_{1} A}}{{m_{1}^{2} - m_{2}^{2} }}} \right]\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)x - 6m_{1}^{3} } \right] + \left[ {\frac{C}{{4m_{2}^{2} }} + \frac{{4m_{1} m_{2} A}}{{^{{\left( {m_{1}^{2} - m_{2}^{2} } \right)^{2} }} }}} \right]\left[ {\left( {\Omega ^{4} - 3m_{1}^{4} } \right)} \right] \\ & f_{4} (x) = - C\left[ {\left( {\Omega^{4} - 3m_{1}^{4} } \right)\frac{{x^{2} }}{2} - 6m_{1}^{3} x} \right] \\ & \quad + \left[ {\frac{ - C}{{2m_{2} }} + \frac{{ - 2m_{1} A}}{{m_{1}^{2} - m_{2}^{2} }}} \right]\left[ {\left( {\Omega^{4} - 3m_{1}^{4} } \right)x - 6m_{1}^{3} } \right] + \left[ {\frac{ - C}{{4m_{2}^{2} }} + \frac{{4m_{1} m_{2} A}}{{^{{\left( {m_{1}^{2} - m_{2}^{2} } \right)^{2} }} }}} \right]\left[ {\left( {\Omega^{4} - 3m_{1}^{4} } \right)} \right] \\ & f_{5} = \frac{{4m_{1} \pi }}{{\left( {m_{1}^{2} - \pi^{2} } \right)^{2} }} + \frac{{4m_{2} \pi }}{{\left( {m_{2}^{2} - \pi^{2} } \right)^{2} }} + \frac{{6P_{{\text{o}}} \pi^{3} }}{{\pi^{4} - Q\pi^{2} + \left( {K -\Omega ^{4} } \right)}} \\ & f_{6} = \frac{{6P_{{\text{o}}} }}{{\pi^{4} - Q\pi^{2} + \left( {K -\Omega ^{4} } \right)}}\left( {\Omega ^{4} - 3\pi^{4} } \right)\left[ {\frac{{ - 2m_{1} }}{{m_{1}^{2} + \pi^{2} }} + \frac{{ - 2m_{2} }}{{m_{2}^{2} + \pi^{2} }}} \right] \\ & A = \frac{1}{{2m_{1} \left( {m_{1}^{2} - m_{2}^{2} } \right)}} \\ & C = \frac{1}{{2m_{2} \left( {m_{2}^{2} - m_{1}^{2} } \right)}} \\ \end{aligned}$$
$$W_{{2_{{\text{P}}} }} = \frac{1}{{k -\Omega ^{4} }}\left\{ \begin{gathered} \left( {u_{1} - Q\left( {u_{1}^{\prime\prime} + 2m_{1} u_{1}^{\prime} + m_{1}^{2} u_{1} } \right)} \right){\text{e}}^{{m_{1} X}} \hfill \\ + \left( {u_{2} - Q\left( {u_{2}^{\prime\prime} + 2m_{2} u_{2}^{\prime} + m_{2}^{2} u_{2} } \right)} \right){\text{e}}^{{m_{2} X}} \hfill \\ + \left( {u_{3} - Q\left( {u_{3}^{\prime\prime} + 2m_{3} u_{3}^{\prime} + m_{3}^{2} u_{3} } \right)} \right){\text{e}}^{{m_{3} X}} \hfill \\ + \left( {u_{4} - Q\left( {u_{4}^{\prime\prime} + 2m_{4} u_{4}^{\prime} + m_{4}^{2} u_{4} } \right)} \right){\text{e}}^{{m_{4} X}} \hfill \\ - Q*u_{5} \left( {\pi^{2} x\cos \left( {\pi x} \right) - 2\pi \sin \left( {\pi x} \right)} \right) \hfill \\ - Q\left( {2\pi u_{6}^{\prime} x\cos \left( {\pi x} \right) + \left( {u_{6}^{\prime\prime} - \pi^{2} u_{6} } \right)\sin \left( {\pi x} \right)} \right) \hfill \\ \end{gathered} \right\}$$
(17)

where

$$\begin{aligned} u_{1} & = - 3x\phi_{1} - 6\phi_{7} + x\Omega ^{4} \phi_{13} - c_{01} \left( {3x^{2} m_{1}^{4} + 12xm_{1}^{3} + 6m_{1}^{2} } \right) \\ u_{2} & = - 3x\phi_{2} - 6\phi_{8} + x\Omega ^{4} \phi_{14} - c_{02} \left( {3x^{2} m_{2}^{4} + 12xm_{2}^{3} + 6m_{2}^{2} } \right) \\ u_{3} & = - 3x\phi_{3} - 6\phi_{9} + x\Omega ^{4} \phi_{15} - c_{03} \left( {3x^{2} m_{3}^{4} + 12xm_{3}^{3} + 6m_{3}^{2} } \right) \\ u_{4} & = - 3x\phi_{4} - 6\phi_{10} + x\Omega ^{4} \phi_{16} - c_{04} \left( {3x^{2} m_{4}^{4} + 12xm_{4}^{3} + 6m_{4}^{2} } \right) \\ u_{5} & = - 3\phi_{5} - 6\phi_{11} +\Omega ^{4} \phi_{17} + \frac{{12\pi^{3} P_{o} }}{{}} \\ u_{6} & = - 3\phi_{6} - 6\phi_{12} +\Omega ^{4} \phi_{18} + \frac{{12\pi^{3} P_{o} }}{{}} \\ \phi_{1} & = m_{1}^{4} \phi_{13} + c_{11} \left( {4m_{1}^{3} f_{1}^{\prime} + 6f_{1}^{\prime\prime} } \right) \\ \phi_{2} & = m_{2}^{4} \phi_{14} + c_{12} \left( {4m_{2}^{3} f_{2}^{\prime} + 6f_{2}^{\prime\prime} } \right) \\ \phi_{3} & = m_{3}^{4} \phi_{15} + c_{13} \left( {4m_{3}^{3} f_{3}^{\prime} + 6f_{3}^{\prime\prime} } \right) \\ \phi_{4} & = m_{4}^{4} \phi_{16} + c_{14} \left( {4m_{4}^{3} f_{4}^{\prime} + 6f_{4}^{\prime\prime} } \right) \\ \phi_{5} & = \pi^{4} f_{5} - 4\pi^{3} f_{6} \\ \phi_{6} & = \pi^{4} f_{6} \\ \phi_{7} & = m_{1}^{3} \phi_{13} + 3c_{11} \left( {m_{1}^{2} f_{1}^{\prime} + m_{1} f_{1}^{\prime\prime} } \right) \\ \phi_{8} & = m_{2}^{3} \phi_{14} + 3c_{12} \left( {m_{2}^{2} f_{2}^{\prime} + m_{2} f_{2}^{\prime\prime} } \right) \\ \phi_{9} & = m_{3}^{3} \phi_{15} + 3c_{13} \left( {m_{3}^{2} f_{3}^{\prime} + m_{3} f_{3}^{\prime\prime} } \right) \\ \phi_{10} & = m_{2}^{3} \phi_{16} + 3c_{14} \left( {m_{4}^{2} f_{4}^{\prime} + m_{4} f_{4}^{\prime} } \right) \\ \phi_{11} & = - \pi^{3} f_{6} \\ \phi_{12} & = \pi^{3} f_{5} - 3\pi^{2} f_{6} \\ \end{aligned}$$
$$\begin{aligned} W_{{3_{{\text{P}}} }} & = T_{1} \left[ {{\text{e}}^{{m_{1} X}} G_{1} + {\text{e}}^{{ - m_{1} X}} G_{2} + {\text{e}}^{{m_{2} X}} G_{3} + {\text{e}}^{{ - m_{2} X}} G_{4} } \right] \\ & \quad+ {\text{Re}} \left( {{\text{e}}^{i\pi X} T_{1} G_{5} } \right) + {\text{Im}} \left( {{\text{e}}^{i\pi X} T_{1} G_{6} } \right) + {\text{Re}} \left( {{\text{e}}^{i\pi X} T_{1} G_{7} } \right) + {\text{Im}} \left( {{\text{e}}^{i\pi X} T_{1} G_{8} } \right) \\ & \quad+ {\text{Im}} \left( {{\text{e}}^{i\pi X} T_{1} G_{9} } \right) \\ \end{aligned}$$
(18)
$$\begin{aligned} T_{1} & = \alpha_{1} \left( {1 - \alpha_{2} \left( {D + m_{1} } \right)^{2} + \alpha_{3} \left( {D + m_{1} } \right)^{4} } \right),\alpha_{1} = \frac{1}{{K -\Omega ^{4} }} \\ G_{1} & = F_{1} \psi_{1} + F_{2} \phi_{13} + F_{2} c_{01} ,G_{2} = F_{4} \psi_{2} + F_{5} \phi_{14} + F_{6} c_{02} \\ G_{3} & = F_{7} \psi_{3} + F_{8} \phi_{15} + F_{9} c_{03} ,G_{4} = F_{10} \psi_{4} + F_{11} \phi_{16} + F_{12} c_{04} \\ G_{5} & = F_{13} \psi_{5} ,G_{6} = F_{13} \psi_{6} ,G_{7} = F_{14} \phi_{17} ,G_{8} = F_{14} \phi_{18} \\ G_{5} & = \frac{{P_{o} }}{{\pi^{4} - Q\pi^{2} + \left( {K -\Omega ^{4} } \right)}}F_{15} \\ F_{1} & = - 3X\left( {D + m_{1} } \right)^{4} - 6\left( {D + m_{1} } \right)^{3} + \Omega^{4} X \\ F_{2} & = - 3X^{2} \left( {D + m_{1} } \right)^{4} - 12X\left( {D + m_{1} } \right)^{3} - 6\left( {D + m_{1} } \right)^{3} \\ F_{3} & = - X^{3} \left( {D + m_{1} } \right)^{4} - 6X^{2} \left( {D + m_{1} } \right)^{3} - 6X\left( {D + m_{1} } \right)^{3} \\ F_{4} & = - 3X\left( {D - m_{1} } \right)^{4} - 6\left( {D - m_{1} } \right)^{3} + \Omega^{4} X \\ F_{5} & = - 3X^{2} \left( {D - m_{1} } \right)^{4} - 12X\left( {D - m_{1} } \right)^{3} - 6\left( {D - m_{1} } \right)^{3} \\ F_{6} & = - X^{3} \left( {D - m_{1} } \right)^{4} - 6X^{2} \left( {D - m_{1} } \right)^{3} - 6X\left( {D - m_{1} } \right)^{3} \\ F_{7} & = - 3X\left( {D + m_{2} } \right)^{4} - 6\left( {D + m_{2} } \right)^{3} + \Omega^{4} X \\ F_{8} & = - 3X^{2} \left( {D + m_{2} } \right)^{4} - 12X\left( {D + m_{2} } \right)^{3} - 6\left( {D + m_{2} } \right)^{3} \\ F_{9} & = - X^{3} \left( {D + m_{2} } \right)^{4} - 6X^{2} \left( {D + m_{2} } \right)^{3} - 6X\left( {D + m_{2} } \right)^{3} \\ F_{10} & = - 3X\left( {D - m_{2} } \right)^{4} - 6\left( {D - m_{2} } \right)^{3} + \Omega^{4} X \\ F_{11} & = - 3X^{2} \left( {D - m_{2} } \right)^{4} - 12X\left( {D - m_{2} } \right)^{3} - 6\left( {D - m_{2} } \right)^{3} \\ F_{12} & = - X^{3} \left( {D - m_{2} } \right)^{4} - 6X^{2} \left( {D - m_{2} } \right)^{3} - 6X\left( {D - m_{2} } \right)^{3} \\ F_{13} & = - 3X\left( {D + i\pi } \right)^{4} - 6\left( {D + i\pi } \right)^{3} +\Omega ^{4} X \\ F_{14} & = - 3X^{2} \left( {D + i\pi } \right)^{4} - 12X\left( {D + i\pi } \right)^{3} - 6\left( {D + i\pi } \right)^{3} \\ F_{15} & = - X^{3} \left( {D + i\pi } \right)^{4} - 6X^{2} \left( {D + i\pi } \right)^{3} - 6X\left( {D + i\pi } \right)^{3} \\ \end{aligned}$$

The total non-dimensional lateral deflection is:

$$W = W_{{\text{o}}} + \varepsilon W_{1} + \varepsilon^{2} W_{2} + \varepsilon^{3} W_{3} + \cdots$$
(18)

The natural frequency and buckling load can be calculated as:

For a free non-prismatic beam on elastic foundation:

$$m_{1} = m_{2} = \pm i\beta_{\Omega } ,m_{3} = m_{4} = \pm \beta_{\Omega } ,\beta_{\Omega }^{4} =\Omega ^{4} - K$$
(19)

For a non-prismatic beam subjected to axial load:

$$m_{1} = m_{2} = 0,m_{3} = m_{4} = \pm i\beta_{Q} ,\beta^{2} = Q$$

Applying the boundary conditions to obtain four linear equations in the constants \(C_{i1} ,C_{i2} ,C_{i3} ,C_{i4}\) in the form of:

$$\left[ G \right]_{4*4} \left[ C \right]_{4*1} = 0$$
(20)

To have a non-trivial solution of Eq. (20), determinant of \(\left[ G \right]_{4*4}\) must be equal to zero. This will lead the critical value of \(\beta,\) and then, we can find the non-dimensional frequency and buckling load parameters for non-prismatic beam as

$$\Omega =\Omega _{{\text{o}}} \left( {1 + \varepsilon + \epsilon^{2} + \epsilon^{3} } \right)^{{{\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 4}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{$4$}}}} \quad {\text{Rectangular}}$$
(21)
$$\Omega =\Omega _{{\text{o}}} \left( {1 + \varepsilon + \epsilon^{2} + \epsilon^{3} + \epsilon^{4} } \right)^{{{\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 4}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{$4$}}}} \quad {\text{Circular}}$$
(22)
$$Q = Q_{{\text{o}}} \left( {1 + \varepsilon + \epsilon^{2} + \epsilon^{3} } \right)\quad {\text{Rectangular}}$$
(23)
$$Q = Q_{{\text{o}}} \left( {1 + \varepsilon + \epsilon^{2} + \epsilon^{3} + \epsilon^{4} } \right)\quad {\text{Circular}}$$
(24)

where \(\Omega _{{\text{o}}} \;{\text{and}}\;Q_{{\text{o}}}\) are the non-dimensional frequency and buckling load parameters for uniform beam.

2.2 Differential quadrature method (DQM)

The DQM requires to divide the domain of the problem into N pointes. Then, the derivatives at any point are approximated by a weighted linear summation of all the functional values along the domain, as follows [14,15,16]:

$$f_{x} \left( {x_{i} } \right) = \left. {\frac{{{\text{d}}f}}{{{\text{d}}x}}} \right|_{{x_{i} }} = \sum\limits_{j = 1}^{N} {A_{ij} \cdot f\left( {x_{j} } \right),\quad i = 1,2, \ldots ,N}$$
(25)

where Aij represented the weighting coefficient and N is the number of grid points in the whole domain. The weighting coefficient can be determined by making use of Lagrange formula as follows:

$$g_{k} \left( x \right) = \frac{M(x)}{{\left( {x - x_{k} } \right) \cdot M^{(1)} \left( {x_{k} } \right)}},\quad k = 1,2, \ldots N$$
(26)

where

$$M(x) = \left( {x - x_{1} } \right)\left( {x - x_{2} } \right) \ldots \left( {x - x_{N} } \right)$$
(27)
$$M^{(1)} \left( {x_{i} } \right) = \prod\limits_{k = 1,k \ne i}^{N} {\left( {x_{i} - x_{N} } \right)}$$
(28)

By applying Eq. (28) at N grid points, the following algebraic formulations to compute the weighting coefficients Aij can be obtained.

$$A_{ij} = \frac{1}{{x_{i} - x_{j} }}\prod\limits_{k = 1,k \ne i,j}^{N} {\frac{{\left( {x_{i} - x_{k} } \right)}}{{\left( {x_{j} - x_{k} } \right)}}} ,\quad j \ne i$$
(29)
$$A_{ii} = \sum\limits_{k = 1,k \ne i}^{N} {\frac{1}{{x_{i} - x_{k} }}}$$
(30)

For calculating the weighting coefficients of mth order

$$\left\lfloor {A^{(m)} } \right\rfloor = \left\lfloor {A^{(1)} } \right\rfloor \cdot \left\lfloor {A^{(m - 1)} } \right\rfloor = \left\lfloor {A^{(m - 1)} } \right\rfloor \cdot \left\lfloor {A^{(1)} } \right\rfloor \quad m = 2,3, \ldots ,N - 1$$
(31)

The accuracy of DQM is affected by choosing of the number, N, and type of the sampling points, xi. The optimal selection of the sampling points in the vibration problems, is the normalized Gauss–Chebyshev–Lobatto points [12, 13],

$$X\left( i \right) = \frac{1}{2}\left[ {1 - \cos \left( {\frac{i - 1}{{N - 1}}\pi } \right)} \right],\quad i = 1,2,3, \ldots ,N$$
(32)

Applying the DQM scheme to the non-dimensional governing Eq. (8) yields:

$$\begin{aligned} & S\left( {X_{i} } \right)\left( {\mathop \sum \limits_{j = 1}^{N} C_{ij}^{\left( 4 \right)} W_{j} } \right) + 2*S^{\left( 1 \right)} \left( {X_{i} } \right)\left( {\mathop \sum \limits_{j = 1}^{N} C_{ij}^{\left( 3 \right)} W_{j} } \right) + S^{\left( 2 \right)} \left( {X_{i} } \right)\left( {\mathop \sum \limits_{j = 1}^{N} C_{ij}^{\left( 2 \right)} W_{j} } \right) \\ & \quad + Q\left( {\mathop \sum \limits_{j = 1}^{N} C_{ij}^{\left( 2 \right)} W_{j} } \right) + \left( {K - \Omega^{4} } \right)W_{i} = P\left( {X_{i} } \right)\quad i = 1,2,3, \ldots ,N \\ \end{aligned}$$
(33)

where \(W_{i}\) is the functional value at the grid points \(X_{i}\), \(C_{ij}^{\left( 2 \right)}\),\(C_{ij}^{\left( 3 \right)}\) and \(C_{ij}^{\left( 4 \right)}\) are the weighting coefficient matrices of the second-, third- and fourth-order derivatives. The result is a linear system of equations of N unknowns that can be written in a matrix form as

$$\left[ {G_{o} } \right]_{N*N} \left[ w \right]_{N*1} = \lambda \left[ {R_{o} } \right]_{N*N} \left[ w \right]_{N*1}$$
(34)

Four boundary conditions are needed to solve the problem. For clamped and simply supported end conditions, the discrete boundary conditions using the DQM can be written as:

$$W_{1} = 0$$
(35)
$$\sum\limits_{k = 1}^{N} {C_{1,k}^{{(n_{0} )}} \cdot W_{k} = 0}$$
(36)
$$W_{N} = 0$$
(37)
$$\sum\limits_{k = 1}^{N} {C_{N,k}^{{(n_{1} )}} \cdot W_{k} = 0}$$
(38)

where n0, n1 may be taken either 1 or 2; value of 1 represents simply supported and value of 2 represents clamped end. the boundary condition can be written in a matrix form as

$$\left[ {G_{1} } \right]_{4*N} \left[ w \right]_{N*1} = \left[ {R_{1} } \right]_{4*N} \left[ w \right]_{N*1}$$
(39)

Combining governing Eqs. (33) and boundary conditions (34), we get

$$\left[ {\begin{array}{*{20}c} {\left[ {G_{{\text{o}}} } \right]_{N*N} } \\ {\left[ {G_{1} } \right]_{4*N} } \\ \end{array} } \right]\left[ w \right]_{N*1} = \upsilon \left[ {\begin{array}{*{20}c} {\left[ {R_{{\text{o}}} } \right]_{N*N} } \\ {\left[ {R_{1} } \right]_{4*N} } \\ \end{array} } \right]\left[ w \right]_{N*1}$$
(40)

; introducing the Lagrange multiplier approach [17], Eq. (40) can be modified to square matrix as

$$\left[ {\begin{array}{*{20}c} {\left[ {G_{{\text{o}}} } \right] } & {\left[ {G_{1} } \right]^{T} } \\ {\left[ {G_{1} } \right] } & 0 \\ \end{array} } \right]_{{\left( {N + 4} \right)*\left( {N + 4} \right)}} \left[ {\begin{array}{*{20}c} w \\ l \\ \end{array} } \right]_{{\left( {N + 4} \right)*1}} = \upsilon \left[ {\begin{array}{*{20}c} { \left[ {R_{{\text{o}}} } \right] } & {\left[ {R_{1} } \right]^{T} } \\ {\left[ {R_{1} } \right]} & 0 \\ \end{array} } \right]_{{\left( {N + 4} \right)*\left( {N + 4} \right)}} \left[ {\begin{array}{*{20}c} w \\ l \\ \end{array} } \right]_{{\left( {N + 4} \right)*1}}$$
(41)

where \(l\) is the Lagrange multiplier and \(\upsilon\) is the Eigen-value to be obtained. Rewrite Eq. (41) in the form

$$\left[ {\begin{array}{*{20}c} {\left[ {G_{{\text{o}}} } \right] } & {\left[ {G_{1} } \right]^{T} } \\ {\left[ {G_{1} } \right] } & 0 \\ \end{array} } \right]\left[ {\begin{array}{*{20}c} { \left[ {R_{{\text{o}}} } \right] } & {\left[ {R_{1} } \right]^{T} } \\ {\left[ {R_{1} } \right]} & 0 \\ \end{array} } \right]^{ - 1} \left[ {\begin{array}{*{20}c} w \\ l \\ \end{array} } \right] = \gamma \left[ {\begin{array}{*{20}c} w \\ l \\ \end{array} } \right]$$
(42)

Solve the Eigen-value problem in Eq. (42) to find non-dimensional frequency parameter and non-dimensional buckling load parameter. This approach can be used to easily implement any boundary condition without modifying weighting coefficients of the DQM.

3 Results

3.1 Verification for prismatic beams

The present PM and DQM results for uniform beam with different end support, listed in Table 1, are identical as compared with the analytical solution in [15]. It is clear that there is a very good agreement between the present results and the previous analytic solution for uniform beam. On the other hand, the results of non-dimensional frequency obtained from the present DQM for non-uniform beam, listed in Table 2, are compared to numerical values obtained in [11] giving a good agreement with error less than 0.5%.

Table 1 Comparison between present DQM, PM, and analytical values for the first three non-dimensional frequencies of uniform beams
Table 2 Comparison between DQM and numerical values for the first three non-dimensional frequencies of non-uniform beams

3.2 Results of vibration parameters for non-prismatic beams

The DQM is applied with different numbers of grid points. The optimum number of grid points is found to be 15 points. After 15 grid points, the results of DQM are fixed [16]. The DQM showed good ability for the analysis of non-prismatic beam with different end supports. The results of non-dimensional frequency and buckling parameters are obtained using the present PM and DQM for clamped–clamped ends, and different values of \(\gamma\) are presented in Tables 3, 4, 5 and 6. The results showed agreement between results for \(< 0.6\). If \(\gamma = 1\), the relative error reached about 20%. The accuracy of the PM depends on the value of \(\gamma\), so as \(\gamma\) increases, the effect of neglected terms affects the results of PM.

Table 3 Results of PM and DQM for the non-dimensional frequencies of clamped–clamped non-prismatic rectangular beam
Table 4 Results of PM and DQM for the non-dimensional frequencies of clamped–clamped non-prismatic circular beam
Table 5 Comparison between PM and DQM for the non-dimensional buckling parameter of clamped–clamped non-prismatic rectangular beam
Table 6 Results of DQM and PM for the non-dimensional buckling parameter of clamped–clamped non-prismatic circular beam

Figure 2 presents the percentage relative error between the PM and DQM for the first non-dimensional frequency of non-prismatic beam calculated by the following equation:

$$\% {\text{Absolute}}\;{\text{relative}}\;{\text{error}} = \left| {\frac{{\Omega _{{{\text{DQM}}}} -\Omega _{{{\text{PM}}}} }}{{\Omega _{{{\text{DQM}}}} }}} \right| \times 100$$
Fig. 2
figure 2

Relative error between PM and DQM for the first non-dimensional frequency of non-prismatic beam

3.3 Deflection for non-prismatic beams

The deflection of non-prismatic beam with circular cross section under the effect of constant vertical load is plotted for two different cases of end supports and different values of \(\gamma\) as shown in Figs. 3 and 4. As \(\gamma\) increases, the maximum deflection decreases and the location of maximum deflection displaced toward the smaller moment of inertia as predicted.

Fig. 3
figure 3

Deflection of clamped–clamped non-prismatic beam

Fig. 4
figure 4

Deflection of clamped–simply supported non-prismatic beam

4 Discussion

Results of DQM and PM for uniform beams, as shown in Table 1, showed perfect agreement with analytical solution for non-dimensional frequency parameter. DQM results for non-uniform beams, as shown in Table 2, proved its perfect ability in prediction of non-dimensional frequency parameter as compared with numerical results given in [11]. The error for listed cases is less than 0.1%. Good results of DQM required fifteen discretization points.

Both methods are used to predict non-dimensional frequency and buckling parameters for non-uniform beams with rectangular or circular cross sections, as shown in Tables 3, 4, 5 and 6. All cases have clamped–clamped ends. The results showed good ability of PM in prediction of non-dimensional frequency and non-dimensional buckling load when the height or the radius increasing ratios are less than 60%. The results of PM are not reliable for higher increasing ratio of height or radius.

Deflection results for clamped–clamped (Fig. 3), clamped–simply supported (Fig. 4), non-uniform beams subjected to a constant uniform load showed good ability of PM in prediction of lateral deflection and the effect of radius increasing ratio on the location of maximum lateral deflection along the beam length. As the increasing ratio of radius increases, the location of maximum lateral deflection moves toward the position of maximum height. Maximum deflection decreases as the increasing ratio of radius increases.

5 Conclusions

The present paper presents two mathematical techniques to investigate the buckling load as well as the natural frequency under axial and transverse load for non-prismatic beams. The deflection of the beam resting on an elastic foundation under transverse distributed and axial loads is also obtained for different end supports. From the present study, the following points could be concluded:

  • The PM and DQM succeeded powerfully for investigating the buckling load as well as the natural frequency for non-prismatic beam.

  • The percentage of relative error between DQM and PM doesn’t exceed than 5% for \(\gamma < 0.6\).

  • As \(\gamma\) increases, the maximum deflection decreases and the location of maximum deflection displaced toward the smaller moment of inertia.

Availability of data and material

Not applicable.

Abbreviations

PM:

Perturbation method

DQM:

Differential quadrature method

References

  1. Bert CW, Malik M (1996) Differential quadrature method in computational mechanics: a review. ASME Appl Mech Rev 49(1):1–28

    Article  Google Scholar 

  2. Hassan HN, El-Tawil MA (2012) A new technique of using homotopy analysis method for second order nonlinear differential equations. Appl Math Comput 219:708–728

    Google Scholar 

  3. Kukla S, Zamojska I (2005) Application of the Green’s function method in free vibration analysis of non-uniform beams. Sci Res Inst Math Comput Sci 1(4):87–94

    Google Scholar 

  4. Yeh Y-L, Jang M-J, Wang C-C (2006) Analyzing the free vibrations of a plate using finite difference and differential transformation method. Appl Math Comput 178:493–501

    Google Scholar 

  5. Taha MH, Abohadima S (2008) Mathematical model for vibrations of non-uniform flexural beams. EngMech 15(1):3–11

    Google Scholar 

  6. Xiao JR (2001) Boundary element analysis of unilateral supported reissner plates on elastic foundations. Comput Mech 27:1–10

    Article  Google Scholar 

  7. Qaisi MI (1997) A power series solution solution for the non-linear vibration of beams. J Sound Vib 199(4):587–594

    Article  Google Scholar 

  8. Zamojska I (2014) Solution of differential equation for the Euler Bernoulli beam. J Appl Math Comput Mech 13(4):157–162

    Article  Google Scholar 

  9. Cekus D (2012) Free vibration of a cantilever tapered Timoshenko beam, Scientific Research of the Institute of Mathematics and Computer. Science 4(11):11–17

    Google Scholar 

  10. Korabathina R, Koppanati MS (2016) Linear free vibration analysis of tapered Timoshenko beams using coupled displacement field method. Math Models Eng 2(1):27–33

    Google Scholar 

  11. Sato K (1999) Transverse vibrations of linearly tapered beams with ends restrained elastically against rotation subjected to axial force. Int J Mech Sci 22:109–115

    Article  Google Scholar 

  12. Jansen EL (2008) A perturbation method for nonlinear vibrations of imperfect structures: application to cylindrical shell vibrations. Int J Solids Struct 45(3–4):1124–1145

    Article  Google Scholar 

  13. Du H, Lim MK, Lin RM (1994) Application of generalized differential quadrature method to structural problems. Int J Numer Eng 37:1881–1896

    Article  Google Scholar 

  14. Chen S, Du H (1997) A generalized approach for implementing general boundary conditions in the GDQ free vibration analysis of plates. Int J Solids Struct 34:837

    Article  Google Scholar 

  15. Chen S, Du H (1997) Implementation of clamped and simply supported boundary conditions in the GDQ free vibration analysis of beams and plates. Int J Solids Struct 34:819

    Article  Google Scholar 

  16. Chang S (2000) Differential Quadrature and Its Application in Engineering. Springer

    Google Scholar 

  17. Wilson EL (2002) Three Dimensional Static and Dynamic Analysis of Structures, Computers and Structures. Inc

    Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

MA team leader who initiates the concept of the paper derived the governing equations. SA performed the analytical solution by perturbation method. ME applied the numerical method of differential quadrature and wrote the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Mohamed Elshabrawy.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Elshabrawy, M., Abdeen, M.A. & Beshir, S. Analytic and numeric analysis for deformation of non-prismatic beams resting on elastic foundations. Beni-Suef Univ J Basic Appl Sci 10, 57 (2021). https://doi.org/10.1186/s43088-021-00144-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s43088-021-00144-5

Keywords