Analytical solution for strongly nonlinear vibration of a stringer shell

Abstract

In this paper, a new approach called the global residue harmonic balance technique is applied to construct analytical solutions for nonlinear oscillators that appear in cylindrical shells. Periodic solution is analytically proved, and consequently, the relation between amplitude and frequency is obtained in an analytical form. Numerical and analytical comparisons with the Runge–Kutta method and previously existing homotopy perturbation method (HPM) are given, to show the stability, quality, and efficiency of the present approach. The GRHBM provides an outstanding agreement with the numerical results for different amplitude values, leading to sufficiently accurate solutions for nonlinear oscillators.

Introduction

The study of nonlinear vibration problems has received considerable attention because of its great importance in many areas of physical science, engineering, and various mechanical systems. Investigating the dynamic response of beams is one of the most important parts of the design process of structures. Nonlinear vibration of cylindrical shells is frequently employed in mechanical engineering; hence, it is critical to address free vibration in nonlinear form. It is extremely difficult to obtain accurate solutions to nonlinear problems, and it is sometimes even more difficult to obtain an analytical approximation. In general, the solution of nonlinear differential equations is classified into two main types, numerical and analytical methods. Traditional perturbation methods for solving nonlinear equations are restricted to weakly nonlinear situations.^1,2 To overcome this problem, many researchers have solved nonlinear problems using new analytical techniques to find approximate solutions, for example, harmonic balance method,^3–6 energy balance method,^7–9 homotopy perturbation method,^10–14 Li–He’s modified homotopy perturbation method, ¹⁵ Hamiltonian approach,^16,17 coupled homotopy–variational approach,^18–20 variational iteration algorithms^21–24 and global error minimization method, ²⁵ frequency–amplitude formulation, ²⁶ fractional complex transform, ²⁷ modification of Yao–Cheng oscillator, ²⁸ and variational approach.^29,30 Nowadays, much literature focuses on simple fractal oscillators in nonlinear systems.^31–33 Recently, some researchers have studied the analytical solutions for nonlinear vibration equations of beams and shells. For example, Bayat et al. ³⁴ used He’s variational technique, while Pakar and Bayat ³⁵ applied He’s max min approach; also, Bayat et al. ³⁶ used finite element modeling, and Bayat et al. ³⁷ used He’s variational approach to obtain the frequency of the nonlinear equation. In addition, Bayat et al. ³⁸ extended Hamiltonian approach to solve the nonlinear vibration of a stringer shell.

This paper aims to apply the higher-order global residue harmonic balance approach, which is an approximation method that combines the principles of the homotopy perturbation technique with the residue harmonic balance method. The GRHBM was first introduced by Ju and Xue,^39–42 and recently developed through References 43–46. To illustrate the effectiveness of the present method, the nonlinear governing equation of the cylindrical shell is discussed in detail. The GRHBM is very effective for solving nonlinear oscillators. The obtained results are compared with those obtained by the HPM, as well as those predicted by the numerical method. The GRHBM proposes a useful and easy approach for finding the approximate solution of nonlinear differential equations.

The originality of this paper is immediately visible in the presented technique, as well as its simplicity and speed in providing analytically accurate solutions that converge swiftly to the numerical solutions. The advantage of using the GRHBM is that it provides analytically accurate solutions to the nonlinear vibration of a stringer shell directly, avoiding complicated computations. Moreover, this approach is much faster and more successful than other known analytical approaches in the literature for determining many analytical solutions. A straightforward solution approach with high accuracy in the acquired findings from the current problem demonstrates the uniqueness, dependability, and broad application of the suggested analytical technique.

Application of the global residue harmonic balance approach

In this section, we will consider the nonlinear problem of cylindrical shell that appears in many different fields of mechanical engineering and civil.^12,38,47 A schematic view of the closed eccentrically circular cylindrical shell is illustrated in Figure 1:

$\ddot{u}+{\alpha }_1\left(u^2\ddot{u}+{\dot{u}}^{2}u\right)+{\alpha }_{2}u+{\alpha }_3{u}^3+{\alpha }_4{u}^5=0,u\left(0\right)=A,\dot{u}\left(0\right)=0,$

(1)

where the coefficients ${\alpha }_1,{\alpha }_2,{\alpha }_3$ and ${\alpha }_4$ are presented as

${\alpha }_1=0.09375n^4,{\alpha }_2={\epsilon }_1+2{\epsilon }_2{p}^{-2}+{\epsilon }_2{p}^{-4}+n^{-4},$

${\alpha }_3=0.0625+0.5n^2{\epsilon }_1-0.75,{\alpha }_4=0.25n^4,p=m_1/n,{\epsilon }_1=D_1/\left(B_1{R}^2\right),{\epsilon }_2=D_2/\left(B_1{R}^2\right),{\epsilon }_3=D_3/\left(B_1{R}^2\right).$

Figure 1.

Schematic of the closed eccentrically circular stringer shell. ¹²

For more details, see ref. 12.

Using the transformation $\tau =\omega t$ , one can obtain from equation (1).

${\omega }^2{u}^{″}+{\omega }^2{\alpha }_1\left(u^2{u}^{″}+u^{\prime 2}u\right)+{\alpha }_{2}u+{\alpha }_3{u}^3+{\alpha }_4{u}^5=0,u\left(0\right)=A,u^{\prime }\left(0\right)=0.$

(2)

In what follows, the GRHB approach, as one of the novelty approaches, is utilized to get an approximate solution for the nonlinear shaking of a stringer shell up to the second order. We will compare our results with some existing methods, including the HPM and R-K method.

First-order approximation

First, we consider

$u_{\left(0\right)}\left(\tau \right)=A\hspace{0.17em}\mathit{cos}\left(\tau \right),\hspace{1em}{\omega }^2={\omega }_0^2.$

(3)

Putting equation (3) into equation (2) yields

$\begin{array}{l}\frac{1}{16}\left(10A^5{\alpha }_4-8A^3{\alpha }_1{\omega }_0^2+12A^3{\alpha }_3+16A{\alpha }_2-16A{\omega }_0^2\right)\mathit{cos}\left(\tau \right)\\ +\frac{1}{16}\left(5A^5{\alpha }_4-8A^3{\alpha }_1{\omega }_0^2+4A^3{\alpha }_3\right)\mathit{cos}\left(3\tau \right)+\frac{1}{16}{A}^5{\alpha }_4\hspace{0.17em}\mathit{cos}\left(5\tau \right)=0.\end{array}$

(4)

Equating the coefficient of $\mathit{cos}\left(\tau \right)$ in equation (4) to zero and then solving it to obtain

${\omega }_0=\frac{\sqrt{5A^4{\alpha }_4+6A^2{\alpha }_3+8{\alpha }_2}}{2\sqrt{A^2{\alpha }_1+2}}.$

(5)

The residual part of equation (4) throws the coefficients of $\mathit{cos}\left(3\tau \right)$ and $\mathit{cos}\left(5\tau \right)$ given by

$R_0=\left(-\frac{A^3{\alpha }_1{\alpha }_2}{2+A^2{\alpha }_1}+\frac{A^3{\alpha }_3}{4}-\frac{3A^5{\alpha }_1{\alpha }_3}{4\left(2+A^2{\alpha }_1\right)}+\frac{5A^5{\alpha }_4}{16}-\frac{5A^7{\alpha }_1{\alpha }_4}{8\left(2+A^2{\alpha }_1\right)}\right)\mathit{cos}\left(3\tau \right)+\frac{1}{16}{A}^5{\alpha }_4\hspace{0.17em}\mathit{cos}\left(5\tau \right).$

(6)

Second-order approximation

We consider

$u_{\left(1\right)}\left(\tau \right)=u_0\left(\tau \right)+p{u}_1\left(\tau \right),{\omega }^2={\omega }_0^2+p{\omega }_1,$

(7)

Inserting equation (7) into equation (2), we obtain $F_1\left(\tau ,u_1,{\omega }_1\right).$

$\begin{array}{l}{F}_1\left(\tau ,u_1,{\omega }_1\right)=u_0{\alpha }_1{\omega }_0^2\hspace{0.33em}{u}_0^{\prime 2}+u_0^{″}\left({\alpha }_1{u}_0^2{\omega }_0^2+{\omega }_0^2\right)+{\alpha }_4{u}_0^5+{\alpha }_3{u}_0^3+{\alpha }_2{u}_0\\ +p\hspace{0.33em}(u_0^{\prime 2}\left({\alpha }_1{u}_1{\omega }_0^2+{\alpha }_1{u}_0{\omega }_1\right)+u_0^{″}\left({\alpha }_1{u}_0^2{\omega }_1+2{\alpha }_1{u}_1{u}_0{\omega }_0^2+{\omega }_1\right)\\ +u_1^{″}\left({\alpha }_1{u}_0^2{\omega }_0^2+{\omega }_0^2\right)+2{\alpha }_1{u}_0{\omega }_0^2{u}_0^{\prime }{u}_1^{\prime }+5{\alpha }_4{u}_1{u}_0^4+3{\alpha }_3{u}_1{u}_0^2+{\alpha }_2{u}_1)\\ +p^2(u_0^{\prime }{u}_1^{\prime }\left(2{\alpha }_1{u}_1{\omega }_0^2+2{\alpha }_1{u}_0{\omega }_1\right)+u_1^{″}\left({\alpha }_1{u}_0^2{\omega }_1+2{\alpha }_1{u}_1{u}_0{\omega }_0^2+{\omega }_1\right)\\ +u_0^{″}\left({\alpha }_1{u}_1^2{\omega }_0^2+2{\alpha }_1{u}_0{u}_1{\omega }_1\right)+{\alpha }_1{u}_0{\omega }_0^2{u}_1^{\prime 2}+{\alpha }_1{u}_1{\omega }_1{u}_0^{\prime 2}+10{\alpha }_4{u}_1^2{u}_0^3\\ +3u_1^2{u}_0{\alpha }_3)+p^3(2{\alpha }_1{u}_1{\omega }_1{u}_0^{\prime }{u}_1^{\prime }+u_1^{\prime 2}\left({\alpha }_1{u}_1{\omega }_0^2+{\alpha }_1{u}_0{\omega }_1\right)+{\alpha }_1{u}_1^2{\omega }_1{u}_0^{″}\\ +u_1^{″}\left({\alpha }_1{u}_1^2{\omega }_0^2+2{\alpha }_1{u}_0{u}_1{\omega }_1\right)+{\alpha }_3{u}_1^3+10{\alpha }_4{u}_0^2{u}_1^3)+p^4({\alpha }_1{u}_1^2{\omega }_1{u}_1^{″}\\ +5{\alpha }_4{u}_0{u}_1^4+{\alpha }_1{u}_1{\omega }_1{u}_1^{\prime 2})+p^5{\alpha }_4{u}_1^5.\end{array}$

(8)

We assume

$u_1\left(\tau \right)=a_{31}\left(\mathit{cos}\left(\tau \right)-\hspace{0.17em}\mathit{cos}\left(3\tau \right)\right),$

(9)

where $a_{31}$ is the unknown constant and will be identified later.

Consider the coefficient of $p$ and inserting equations (3), (6), and (9) into equation (8), we get the function $F_1\left(\tau ,a_{31},{\omega }_1\right)$

$\begin{array}{l}{F}_1\left(\tau ,a_{31},{\omega }_1\right)=\frac{1}{16}(25a_{31}{A}^4{\alpha }_4+24a_{31}{A}^2{\alpha }_3+16a_{31}{\alpha }_2-16a_{31}{\omega }_0^2-8A^3{\alpha }_1{\omega }_1\\ \hspace{0.33em}-16A{\omega }_1)\hspace{0.17em}\mathit{cos}\left(\tau \right)+\frac{1}{16}(-5a_{31}{A}^4{\alpha }_4+56a_{31}{A}^2{\alpha }_1{\omega }_0^2-12a_{31}{A}^2{\alpha }_3\\ -16a_{31}{\alpha }_2+144a_{31}{\omega }_0^2-8A^3{\alpha }_1{\omega }_1)\hspace{0.17em}\mathit{cos}\left(3\tau \right)+\frac{1}{16}(72a_{31}{A}^2{\alpha }_1{\omega }_0^2\\ -15a_{31}{A}^4{\alpha }_4+-12a_{31}{A}^2{\alpha }_3)\hspace{0.17em}\mathit{cos}\left(5\tau \right)-\frac{5}{16}{a}_{31}{A}^4{\alpha }_4\hspace{0.17em}\mathit{cos}\left(7\tau \right).\end{array}$

(10)

Rewrite equation (10) to be

$\begin{array}{l}{F}_1\left(\tau ,a_{31},{\omega }_1\right)=\frac{1}{16\left(A^2{\alpha }_1+2\right)}[(25a_{31}{A}^6{\alpha }_1{\alpha }_4+24a_{31}{A}^4{\alpha }_1{\alpha }_3+30a_{31}{A}^4{\alpha }_4\\ +16a_{31}{A}^2{\alpha }_1{\alpha }_2+24a_{31}{A}^2{\alpha }_3-\left(8A^5{\alpha }_1^2+32A^3{\alpha }_1+32A\right){\omega }_1)\hspace{0.17em}\mathit{cos}\left(\tau \right)\\ +(65a_{31}{A}^6{\alpha }_1{\alpha }_4+72a_{31}{A}^4{\alpha }_1{\alpha }_3+96a_{31}{A}^2{\alpha }_1{\alpha }_2+170a_{31}{A}^4{\alpha }_4\\ +192a_{31}{A}^2{\alpha }_3+256a_{31}{\alpha }_2-8A^5{\alpha }_1^2{\omega }_1-16A^3{\alpha }_1{\omega }_1)\hspace{0.17em}\mathit{cos}\left(3\tau \right)\\ +(75a_{31}{A}^6{\alpha }_1{\alpha }_4+96a_{31}{A}^4{\alpha }_1{\alpha }_3-30a_{31}{A}^4{\alpha }_4+144a_{31}{A}^2{\alpha }_1{\alpha }_2\\ -24a_{31}{A}^2{\alpha }_3)\hspace{0.17em}\mathit{cos}\left(5\tau \right)+\left(-5a_{31}{A}^6{\alpha }_1{\alpha }_4-10a_{31}{A}^4{\alpha }_4\right)\mathit{cos}\left(7\tau \right)].\end{array}$

(11)

Then, we have

$F_1\left(\tau ,a_{31},{\omega }_1\right)+R_0\left(\tau \right)=0.$

(12)

We set the coefficients of $\mathit{cos}\left(\tau \right)$ and $\mathit{cos}\left(3\tau \right)$ specified in (12) to zero. It uses nonlinear equations.

${\omega }_1=\frac{A^4{\Gamma }_1{\Gamma }_2}{16\left(A^2{\alpha }_1+2\right)\mathrm{\Delta }},a_{31}=\frac{A^3\left(A^2{\alpha }_1+2\right){\Gamma }_1}{2\mathrm{\Delta }},$

(13)

where

$\begin{array}{l}{\Gamma }_1=25A^4{\alpha }_1{\alpha }_4+24A^2{\alpha }_1{\alpha }_3+30A^2{\alpha }_4+16{\alpha }_1{\alpha }_2+24{\alpha }_3,\\ {\Gamma }_2=5A^4{\alpha }_1{\alpha }_4+8A^2{\alpha }_1{\alpha }_3-10A^2{\alpha }_4+16{\alpha }_1{\alpha }_2-8{\alpha }_3,\\ ∆=256{\alpha }_2+20A^8{\alpha }_1^2{\alpha }_4+24A^6{\alpha }_1^2{\alpha }_3+135A^6{\alpha }_1{\alpha }_4+40A^4{\alpha }_1^2{\alpha }_2\\ +156A^4{\alpha }_1{\alpha }_3+170A^4{\alpha }_4+224A^2{\alpha }_1{\alpha }_2+192A^2{\alpha }_3.\end{array}$

The second-order approximation can be expressed as

$\begin{array}{l}{u}_{\left(1\right)}\left(\tau \right)=u_0\left(\tau \right)+u_1\left(\tau \right)=(A+a_{31})\mathit{cos}\left(\tau \right)-a_{31}\hspace{0.17em}\mathit{cos}\left(3\tau \right),\\ {\omega }_{\left(1\right)}^2={\omega }_0^2+{\omega }_1,\tau ={\omega }_{\left(1\right)}t.\end{array}$

(14)

We note that the GRHBM can provide third or higher-order approximations in a similar manner.

Discussion

To confirm the validity of the GRHB method applied in this research, the results are evaluated against previous literature findings and numerical solutions. Figure 2 provides comparisons of the present results and the HPM ¹² with the numerical ones for small and large values of the amplitudes of vibration $A$ ; it is confirmed that the GRHBM presented in this work shows excellent agreement with the numerical Runge–Kutta solution than the other existing one in the literature. When comparing the estimated frequency values with analytically determined ones, it is clear that the recommended second-order approximation GRHBM provides the best results. For small amplitude values, the current approach and the homotopy perturbation method have about the same degree of accuracy. But, for the large values of the amplitude $A$ , the current approach is superior to the aforementioned analytical method. As a result, the GRHBM has a basic method that can be implemented with the necessity of simpler mathematics and is better suited for computer programming. In lower-order approximations, the GRHBM can get acceptable and accurate results.

Figure 2.

Contrasting between the global residue harmonic balance method (black line), the homotopy perturbation method ¹² (red line), and the numerical solution (blue line).

Concluding remarks

The GRHBM has been employed successfully in this study, to obtain analytically accurate solutions to the nonlinear problem of cylindrical shell. The present results are compared with the numerical ones achieved by using the HPM. It is clear from Figure 2 that the present solutions agree more closely with the numerical solutions than the homotopy perturbation method for different values of the amplitude A. Finally, we can conclude that the present method is powerful and simple, and can also be applied to other nonlinear problems without requiring complex mathematical calculations.