Correction of max-min approach for analyzing shock response of strongly nonlinear system

Li, Hong-wei; Wang, Jun

Journal of Vibroengineering

Browse Journal

Submit article

Published: 30 September 2015

Check for updates

Correction of max-min approach for analyzing shock response of strongly nonlinear system

Hong-wei Li¹

Jun Wang²

^{1, 2}Jiangnan University, Wuxi 214122, China

²Jiangsu Key Laboratory of Advanced Food Manufacturing Equipment and Technolog, Wuxi, 214122, China

Corresponding Author:

Jun Wang

Cite the article Download PDF

Downloads 1202

Abstract

A novel method combining max-min approach (MMA) with energy method is proposed for the derivation of the approximate analytical solution for shock response of the strong nonlinear system. MMA is used to evaluate the governing differential equation with nonlinearity terms and to obtain the approximate solution. The maximum shock response displacement of the system is obtained by energy method, and then the approximation solution by MMA is corrected. Three examples are given to illustrate the validity of the method. It is shown that the corrected MMA method always lead to a better approximate solution for different kinds of nonlinear problems. The results provides a new simple but effective method for the analysis of shock response of strong nonlinear problems.

1. Introduction

There are various kinds of nonlinear problems in the engineering field. In most cases, the constitutive relations of systems are strong nonlinear, which are usually difficult to be solved analytically. Various analytical approaches for solving nonlinear differential equations are proposed and applied in the analysis of some engineering problems, such as the homotopy analysis method [1-3], variational iteration method (VIM) [4-6], the method of harmonic balance [7-9], Adomian’s decomposition method [10-12], parameter expansion method [13-15] and max-min approach (MMA) [16-18]. Among these analytical methods, as suggested by S. S. Ganji [19] and Zeng [20], MMA is very effective for solving nonlinear oscillators.

The shock response is very different from free vibration for nonlinear systems. By combining the L-P Perturbation Method with the energy method (EM) and combining the VIM with the EM, Xiang [21] and Chen [22] get good analytical expressions of dropping shock response of nonlinear packaging systems separately. In this study, the shock response of strongly nonlinear systems is considered, and the correction method combining the approximate solution of system applying the MMA with the EM is developed.

2. MMA and EM

2.1. Basic idea of the MMA [17]

Considering the following equation:

1

\{\begin{array}{l} \ddot{u} + u f (u) = 0, \\ u (0) = 0, \dot{u} (0) = v, \end{array}

where $f (u)$ is a nonnegative function. A trial-function is chosen in the following form:

2

u = A s i n ω x,

where $ω$ is the frequency to be determined.

Observing that the square of frequency $ω^{2}$ , has:

3

f_{m i n} < ω^{2} < f_{m a x} .

Hence, it follows that:

4

\frac{f_{m i n}}{1} < ω^{2} < \frac{f_{m a x}}{1} .

According to He Chengtian’s interpolation [23, 24], one can have:

5

ω^{2} = \frac{m f_{m i n} + n f_{m a x}}{m + n} = f_{m i n} + k (f_{m a x} - f_{m i n}),

where, $m$ and $n$ are weighting factors, and $k = n / (m + n)$ .

So the frequency can be approximated as:

6

ω = \sqrt{f_{m i n} + k (f_{m a x} - f_{m i n})} .

Then, one can get the approximate solution for Eq. (1) as:

7

u = A s i n [(f_{m i n} + k (f_{m a x} - f_{m i n}))^{1 / 2} t],

where, $k$ can be easily obtained.

2.2. Basic idea of the EM

For nonlinear system without damping, the kinetic energy of the system transforms into completely elastic potential energy when the deformation of the material achieves the maximum value, $x_{m}$ . Assuming the critical impact velocity of the system is $V$ , the critical kinetic energy of the system can be written as:

8

T = \frac{1}{2} m v^{2} .

For a nonlinear system, $f (x)$ denotes the corresponding restoring force, where $f (x)$ is a nonlinear function of $x$ , depending upon only the function of x. By using the energy method (EM), one can have:

9

\frac{1}{2} m v^{2} = \int_{0}^{x_{m}} f (x) d x .

The maximum deformation $x_{m}$ can be obtained by solving Eq. (9).

3. Correction method of MMA

3.1. Example 1

Considering the duffing system, the governing equation can be written as:

10

\{\begin{array}{l} m \ddot{x} + a x + b x^{3} = 0, \\ x (0) = 0, \dot{x} (0) = V, \end{array}

where the dot denotes the differentiation with respect to the time $t$ , $m$ is the mass of product, $x$ is the displacement of the system, $a$ and $b$ denote the constant coefficient, respectively, and $V$ is the critical impact velocity. The restoring force of the system can be represented as:

11

f (x) = a x + b x^{3} .

By introducing parameters $T = \sqrt{m / a}$ , $L = \sqrt{a / b}$ , and setting $u = x / L$ and $τ = t / T$ , one can obtain the dimensionless equation as:

12

\{\begin{array}{l} u'' + u + u^{3} = 0, \\ u (0) = 0, u^{'} (0) = \sqrt{\frac{b m V^{2}}{a^{2}}}, \end{array}

where the prime denotes the differentiation with respect to the $τ$ . To apply MMA for solving Eq. (12), Eq. (12) is firstly re-written in the following form:

13

u'' + (1 + u^{2}) u = 0 .

One can choose a trial-function in the form:

14

u = A s i n Ω τ,

where $Ω$ is the frequency to be determined.

Observing that the square of frequency, $Ω^{2}$ , one can have:

15

1 < Ω^{2} < 1 + A^{2} .

Hence, it follows that:

16

\frac{1}{1} < Ω^{2} < \frac{1 + A^{2}}{1} .

According to He Chengtian’s interpolation [21, 22], one can have:

17

Ω^{2} = \frac{k + n (1 + A^{2})}{k + n} = 1 + l A^{2},

where $k$ and $n$ are weighting factors, $l = n / (k + n)$ .

So the frequency can be approximated as:

18

Ω = \sqrt{1 + l A^{2}} .

Then, we get the approximate solution of Eq. (12) as:

19

u (τ) = A s i n [(1 + l A^{2})^{1 / 2} τ] .

Combination of Eq. (12) and Eq. (19) leads to:

20

u^{''} + (1 + l A^{2}) u = l A^{2} u - u^{3} .

If, by chance, Eq. (19) is the exact solution, then the right hand side of Eq. (20) is vanishing completely. Since our approach is only an approximation to the exact solution, we set:

21

\int_{0}^{T_{0} / 4} (l A^{2} u - u^{3}) s i n Ω τ d τ = 0,

where $T_{0} = 2 π / Ω$ . Substituting Eq. (19) into Eq. (21), one can obtain:

22

l = 3 / 4 .

Finally the frequency can be obtained as:

23

Ω = \sqrt{1 + \frac{3}{4} A^{2}} .

And the approximate solution of the dimensionless equation by MMA can be written as:

24

u (τ) = A s i n [{(1 + \frac{3}{4} A^{2})}^{1 / 2} τ] .

It is easy for one to get the approximate solution of the original Eq. (10) in the form:

25

x (t) = \sqrt{\frac{a A^{2}}{b}} s i n [{(\frac{4 a + 3 a A^{2}}{4 m})}^{1 / 2} t],

26

ω = {(\frac{4 a + 3 a A^{2}}{4 m})}^{1 / 2},

where $ω$ is the frequency of the original system. The exact frequency reads:

27

ω_{e} = \frac{π}{\int_{0}^{π / 2} {\{\frac{[8 m (a^{2} + 2 m b V^{2}) - 8 m a] c o s^{2} φ}{\{2 m b V^{2} - 2 a [(a^{2} + 2 m b V^{2})^{1 / 2} - a] s i n^{2} φ - [(a^{2} + 2 m b V^{2})^{1 / 2} - a]^{2} s i n^{4} φ\}}\}}^{1 / 2} d φ} .

The velocity response of the system can be expressed as:

28

\dot{x} (t) = \sqrt{\frac{a^{2} A^{2} (4 + 3 A^{2})}{4 b m}} c o s [{(\frac{4 a + 3 a A^{2}}{4 m})}^{1 / 2} t] .

From Eq. (28), one can obtain:

29

\dot{x} (0) = \sqrt{\frac{a^{2} A^{2} (4 + 3 A^{2})}{4 b m}} .

Combining Eq. (10) and Eq. (29), the parameter $A$ is obtained. Substituting it into Eq. (26), one can get the value of frequency $ω$ .

When:

{(\frac{4 a + 3 a A^{2}}{4 m})}^{1 / 2} t = \frac{π}{2},

substituting it into Eq. (24), one can obtain:

30

x_{m} = \sqrt{\frac{a A^{2}}{b}} .

A new method, namely, the corrected MMA (CMMA) is suggested here to obtain a better solution to nonlinear equation, which integrates the MMA with the EM. According to the EM, one have:

31

\frac{1}{2} m V^{2} = \int_{0}^{x_{m}} (a x + b x^{3}) d x .

The exact maximum displacement is obtained as:

32

x_{e m} = {[\frac{(a^{2} + 2 m b V^{2})^{1 / 2} - a}{b}]}^{1 / 2} .

Combining Eq. (25), Eq. (30) and Eq. (32), the corrected solution is obtained as:

33

x_{c} = {[\frac{(a^{2} + 2 m b V^{2})^{1 / 2} - a}{b}]}^{1 / 2} s i n \{{[\frac{a + 3 (a^{2} + 2 m b V^{2})^{1 / 2}}{4 m}]}^{1 / 2} t\},

34

ω_{c} = {[\frac{a + 3 (a^{2} + 2 m b V^{2})^{1 / 2}}{4 m}]}^{1 / 2} .

To illustrate the validity of CMMA for this example, the comparison of results obtained by MMA and CMMA is given in Table 1 and Fig. 1. It is shown that the CMMA tends to give a more precise result than MMA for nonlinear equation.

Table 1Comparison of frequency corresponding to various parameters of system for Example 1

$m = a = b = 1$
$V$	0.01	0.1	1.0	10.0	100.0
$ω_{e x}$	1.0000019	1.0037230	1.2413297	3.2403053	10.0914601
$ω$	1.0000375	1.0037153	1.2247449	3.0289727	9.3329515
$ω_{c}$	1.0000375	1.0037245	1.2446036	3.2989521	10.3110945
$\|(ω_{e x} - ω) / ω_{e x}\|$	0.0000356	0.0000077	0.0133605	0.0652200	0.0751634
$\|(ω_{e x} - ω_{c}) / ω_{e x}\|$	0.0000356	0.0000015	0.0026374	0.0180992	0.0217644

3.2. Example 2

The dynamic equation of this example is shown in Eq. (35):

35

m \ddot{x} + \frac{c x}{1 + d x^{2}} = 0, x (0) = 0, \dot{x} (0) = V .

Similarly, here the dot denotes the differentiation with respect to the time $t$ , $m$ is the mass of product, $x$ is the displacement of the system, $c$ and $d$ denote the constant coefficient, respectively, and $V$ is the critical impact velocity. The restoring force of the system can be represented as:

36

f (x) = \frac{c x}{1 + d x^{2}} .

Fig. 1Comparison of MMA and CMMA with the exact solution for Example 1

By introducing parameters $T = \sqrt{m / c}$ , $L = \sqrt{1 / b}$ , and setting $u = x / L$ and $τ = t / T$ , one can obtain the dimensionless equation and initial conditions of system as following:

37

u^{''} + \frac{u}{1 + u^{2}} = 0, u (0) = 0, u^{'} (0) = \sqrt{\frac{d m V^{2}}{c}},

where the prime denotes the differentiation with respect to the $τ$ .

Eq. (37) can be rewritten in the following form:

38

u'' + (\frac{1}{1 + u^{2}}) u = 0 .

By using $u = A s i n Ω τ$ as a trial function, the maximal and minimal values of $1 / 1 + u^{2}$ are, respectively, 1 and $1 / (1 + A^{2})$ , so one can immediately obtain:

39

\frac{1}{1} < Ω^{2} < \frac{1}{1 + A^{2}} .

According to MMA, one can have:

40

Ω^{2} = \frac{k + n}{k + n (1 + A^{2})} = \frac{1}{1 + l A^{2}},

where $k$ and $n$ are weighting factors, $l = n / (k + n)$ .

So the frequency can be approximated as:

41

Ω = \frac{1}{\sqrt{1 + l A^{2}}} .

Now re-write Eq. (37) in the form:

42

(1 + l A^{2}) u^{''} + u = l A^{2} u^{''} - u^{2} u^{''} .

Setting:

43

\int_{0}^{T_{0} / 4} (l A^{2} u'' - u^{2} u'') s i n Ω τ d τ = 0,

one can obtain:

44

l = 3 / 4 .

The approximate frequency can be obtained as:

45

Ω = \frac{2}{\sqrt{4 + 3 A^{2}}} .

And the approximate solution of the dimensionless equation by MMA can be written as:

46

u (τ) = A s i n [\frac{2}{\sqrt{4 + 3 A^{2}}} τ] .

Similarly, it is easy for one to get the approximate solution of the original Eq. (35) in the form:

47

x (t) = \sqrt{\frac{A^{2}}{d}} s i n [{(\frac{4 c}{4 m + 3 m A^{2}})}^{1 / 2} t],

48

ω = {(\frac{4 c}{4 m + 3 m A^{2}})}^{1 / 2},

where $ω$ is the frequency of the original system. Its exact frequency reads:

49

ω_{e} = \frac{π}{\int_{0}^{π / 2} {\{\frac{4 m (e^{d m V^{2} / c} - 1) c o s^{2} φ}{\{d m V^{2} - c l n [1 + (e^{d m V^{2} / c} - 1) s i n^{2} φ]\}}\}}^{1 / 2} d φ} .

The velocity response of the system can be expressed as:

50

\dot{x} (t) = \sqrt{\frac{4 c A^{2}}{m d (4 + 3 A^{2})}} c o s [{(\frac{4 c}{4 m + 3 m A^{2}})}^{1 / 2} t] .

From Eq. (50), one can obtain:

51

\dot{x} (0) = \sqrt{\frac{4 c A^{2}}{m d (4 + 3 A^{2})}} .

Combining Eq. (35) and Eq. (51), the parameter $A$ is obtained. Substituting it into Eq. (48), one can get the value of frequency $ω$ .

When:

{(\frac{4 c}{4 m + 3 m A^{2}})}^{1 / 2} t = \frac{π}{2},

substituting it into Eq. (47), one can obtain:

52

x_{m} = \sqrt{\frac{A^{2}}{d}} .

By using the EM, we have:

53

\frac{1}{2} m V^{2} = \int_{0}^{x_{m}} \frac{c x}{1 + d x^{2}} d x .

Then, the exact maximum displacement is obtained as:

54

x_{e m} = {[\frac{e^{d m V^{2} / c} - 1}{d}]}^{1 / 2} .

Combining Eq. (47), Eq. (52) and Eq. (54), the corrected solution is obtained as:

55

x_{c} = {[\frac{e^{d m V^{2} / c} - 1}{d}]}^{1 / 2} s i n [{(\frac{4 c}{m + 3 m e^{d m V^{2} / c}})}^{1 / 2} t],

56

ω_{c} = {(\frac{4 c}{m + 3 m e^{d m V^{2} / c}})}^{1 / 2} .

To illustrate the validity of CMMA for this example, comparison of results is given in Table 2 and Fig. 2. It is shown that CMMA is still effective when MMA fails.

Table 2Comparison of frequency corresponding to various parameters of system for Example 2

$m = c = d = 1$
$V$	0.01	0.1	0.5	1.0	2.0
$ω_{e x}$	0.9999626	0.9962557	0.9089836	0.6702093	0.1667291
$ω$	0.9999625	0.9962429	0.9013878	0.50	–
$ω_{c}$	0.9999625	0.9962524	0.9079589	0.6610046	0.1557969
$\|(ω_{e x} - ω) / ω_{e x}\|$	0.0000001	0.0000128	0.0083564	0.2539644	–
$\|(ω_{e x} - ω_{c}) / ω_{e x}\|$	0.0000001	0.0000033	0.0011273	0.0137341	0.0655686

Fig. 2Comparison of MMA and CMMA with the exact solution for Example 2

3.3. Example 3

Considering the following example:

57

m \ddot{x} + (e + f {\dot{x}}^{2}) x = 0, x (0) = 0, \dot{x} (0) = V,

where the dot denotes the differentiation with respect to the time $t$ , $m$ is the mass of product, $x$ is the displacement of the system, $c$ and $d$ denote the constant coefficient, respectively, and $V$ is the critical impact velocity. The restoring force of the system can be represented as:

58

f (x, \dot{x}) = (e + f {\dot{x}}^{2}) x .

By introducing parameters $T = \sqrt{m / e}$ , $L = \sqrt{m / f}$ , and setting $u = x / L$ and $τ = t / T$ , one can obtain the dimensionless equation and initial conditions of system as following forms:

59

u^{''} + (1 + {u^{'}}^{2}) u = 0, u (0) = 0, u^{'} (0) = \sqrt{\frac{f V^{2}}{e}},

where the prime denotes the differentiation with respect to the $τ$ .

By using $A s i n Ω τ$ as a trial function, one can write:

60

\frac{1}{1} < Ω^{2} < \frac{A^{2} Ω^{2} + 1}{1} .

According to the MMA, one can have:

61

Ω^{2} = \frac{k + n}{k + n (1 - A^{2})} = \frac{1}{1 - l A^{2}},

where $k$ and $n$ are weighting factors, $l = n / (k + n)$ .

The frequency can be approximated as:

62

Ω = \frac{1}{\sqrt{1 - l A^{2}}} .

Now re-write Eq. (59) in the following form:

63

(1 - l A^{2}) u^{''} + u = - l A^{2} u^{''} - {u^{'}}^{2} u .

Setting:

64

\int_{0}^{T_{0} / 4} (- l A^{2} u'' - {u'}^{2} u) s i n Ω τ d τ = 0

one can obtain:

65

l = 1 / 4

The approximate frequency can be obtained as:

66

Ω = \frac{2}{\sqrt{4 - A^{2}}} .

And the approximate solution of the dimensionless equation by MMA can be written as:

67

u (τ) = A s i n [\frac{2}{\sqrt{4 - A^{2}}} τ] .

The approximate solution of the original Eq. (56) can be obtained in the form:

68

x (t) = \sqrt{\frac{m A^{2}}{f}} s i n [{(\frac{4 e}{4 m - m A^{2}})}^{1 / 2} t],

69

ω = {(\frac{4 e}{4 m - m A^{2}})}^{1 / 2},

where $ω$ is the frequency of the original system.

The velocity response of the system can be expressed as:

70

\dot{x} (t) = \sqrt{\frac{4 e A^{2}}{f (4 - A^{2})}} c o s [{(\frac{4 e}{4 m - m A^{2}})}^{1 / 2} t] .

From Eq. (70), one can obtain:

71

\dot{x} (0) = \sqrt{\frac{4 e A^{2}}{f (4 - A^{2})}} .

Combining Eq. (56) and Eq. (70), the parameter $A$ is obtained. Substituting it into Eq. (68), one can get the value of frequency $ω$ .

When:

{(\frac{4 e}{4 m - m A^{2}})}^{1 / 2} t = \frac{π}{2},

substituting it into Eq. (68), one can obtain:

72

x_{m} = \sqrt{\frac{m A^{2}}{f}} .

By using the EM, we have:

73

\frac{1}{2} m V^{2} = \int_{0}^{T / 4} (e + f V^{2} c o s^{2} ω t) x_{m} s i n ω t V c o s ω t d t .

Note that the maximum displacement, $x_{m}$ , is approximated by EM in this case.

The approximate maximum displacement by EM is obtained as:

74

x_{a m} = {[\frac{2 m V^{2}}{2 e + f V^{2}}]}^{1 / 2} .

Combining Eq. (48), Eq. (72) and Eq. (74), the corrected solution is obtained as:

75

x_{c} = {[\frac{2 m V^{2}}{2 e + f V^{2}}]}^{1 / 2} s i n [{(\frac{4 e^{2} + 2 e f V^{2}}{4 e m + f m V^{2}})}^{1 / 2} t],

76

ω_{c} = {(\frac{4 e^{2} + 2 e f V^{2}}{4 e m + f m V^{2}})}^{1 / 2} .

The exact frequency cannot be obtained for this example. We showed comparison of $x (t)$ in Figs. 3(a)-(d). These figures are related to parameters $m = e = f =$ 1, $V =$ 0.01, $m = e = f =$ 1, $V = 0.1$ , $m = e = f =$ 1, $V = 0.5$ and $m = e = f =$ 1, $V = 1.0$ , respectively. It is shown that the CMMA can give a more precise result than MMA for the nonlinear problem although the exact maximum displacement response cannot be obtained by the EM.

Fig. 3Comparison of x(t) in case of m=e=f=1 for Example 3 (Solid line: CMMA; Dash-dotted line: MMA; Symbols: Runge-Kutta)

a) $V = 0.01$

b) $V = 0.01$

c) $V = 0.5$

d) $V = 1.0$

4. Conclusion

In this paper, the max-min approach (MMA) which is widely applied for various kinds of nonlinear equations is combined with energy method (EM), a new technique called CMMA is proposed and applied to analysis the shock response of strongly nonlinear system.

Three examples were given to illustrate the validity of the method. The first one shows CMMA is better than MMA for strongly nonlinear system. And it can be seen that CMMA is still effective when the MMA fails from the second example. Although sometimes the exact maximum displacement response cannot be obtained by the EM, the approximate solution is better by CMMA and closed to the numerical solution.

Overall, it can be concluded that CMMA has a great application respect in strong nonlinear problems.

References

Liao Shijun On the homotopy analysis method for nonlinear problems. Applied Mathematics and Computation, Vol. 147, Issue 2, 2004, p. 499-513.

Search CrossRef
Zou Li, Zong Zhi, Wang L, et al. Solving the discrete KdV equation with homotopy analysis method. Physics Letters A, Vol. 370, Issue 3, 2007, p. 287-294.

Search CrossRef
Liao Shi-jun Homotopy Analysis Method in Nonlinear Differential Equations. Springer, London, 2012, p. 26-38.

Search CrossRef
He Ji-huan, Wu Xu-hong Variational iteration method: new development and applications. Computers and Mathematics with Applications, Vol. 54, 2007, p. 881-894.

Search CrossRef
Guo-cheng Wu, Dumitru Baleanu Variational iteration method for the Burgers’ flow with fractional derivatives-new Lagrange multipliers. Applied Mathematical Modelling, Vol. 37, Issue 9, 2013, p. 6183-6190.

Search CrossRef
Khader M. M. Numerical and theoretical treatment for solving linear and nonlinear delay differential equations using variational iteration method. Arab Journal of Mathematical Sciences, Vol. 19, Issue 2, 2013, p. 243-256.

Search CrossRef
Hu H., Tang J. H. Solution of a Duffing-harmonic oscillator by the method of harmonic balance. Journal of Sound and Vibration, Vol. 294, Issue 3, 2006, p. 637-639.

Search CrossRef
Al-Shyyab A., Kahraman A. Non-linear dynamic analysis of a multi-mesh gear train using multi-term harmonic balance method: sub-harmonic motions. Journal of Sound and Vibration, Vol. 279, Issue 1, 2005, p. 417-451.

Search CrossRef
Ghadimi M., Kaliji H. D. Application of the harmonic balance method on nonlinear equations. World Applied Sciences Journal, Vol. 22, Issue 4, 2013, p. 532-537.

Search CrossRef
Abbasbandy S. A numerical solution of Blasius equation by Adomian’s decomposition method and comparison with homotopy perturbation method. Chaos, Solitons and Fractals, Vol. 31, Issue 1, 2007, p. 257-260.

Search CrossRef
Rèpaci A. Nonlinear dynamical systems: on the accuracy of Adomian’s decomposition method. Applied Mathematics Letters, Vol. 3, Issue 4, 1990, p. 35-39.

Search CrossRef
Jiao Y. C., Yamamoto Y., Dang C., et al. An after treatment technique for improving the accuracy of Adomian’s decomposition method. Computers and Mathematics with Applications, Vol. 43, Issue 6, 2002, p. 783-798.

Search CrossRef
Shou Da-hua, He Ji-huan Application of parameter-expanding method to strongly nonlinear oscillators. International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 8, Issue 1, 2007, p. 121-124.

Search CrossRef
Sedighi H. M., Shirazi K. H. Vibrations of micro-beams actuated by an electric field via parameter expansion method. Acta Astronautica, Vol. 85, 2013, p. 19-24.

Search CrossRef
Molaei H., Kheybari S. A Numerical solution of classical van der Pol-Duffing oscillator by He’s parameter-expansion method. International Journal of Contemporary Mathematical Sciences, Vol. 8, Issue 15, 2013, p. 709-714.

Search CrossRef
Kaliji H. D., Ghadimi M., Barari A. Analytical solutions to non-linear mechanical oscillation problems. Journal of Vibroengineering, Vol. 13, Issue 2, 2011, p. 133-143.

Search CrossRef
He Ji-huan Max-min approach to nonlinear oscillators. International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 9, Issue 2, 2008, p. 207-210.

Search CrossRef
Ganji S. S., Ganji D. D., Davodi A. G., et al. Analytical solution to nonlinear oscillation system of the motion of a rigid rod rocking back using max-min approach. Applied Mathematical Modelling, Vol. 34, Issue 9, 2010, p. 2676-2684.

Search CrossRef
Zeng D. Q., Lee Y. Y. Analysis of strongly nonlinear oscillator using the max-min approach. International Journal of Nonlinear Sciences and Numerical Simulation, Vol. 10, Issue 10, 2009, p. 1361-1368.

Search CrossRef
Zeng De-qiang Nonlinear oscillator with discontinuity by the max-min approach. Chaos, Solitons and Fractals, Vol. 42, Issue 5, 2009, p. 2885-2889.

Search CrossRef
Xiang Hong, Tang Bo-sen Application of L-P perturbation method in dropping shock problem. Journal of Vibration and Shock, Vol. 21, Issue 1, 2002, p. 39-42.

Search CrossRef
Chen An-jun Application of variational iteration method with energy method for nonlinear equation arisen from packaging system. Mathematical Problems in Engineering, 2013.

Search CrossRef
He Ji-huan He Chengtian’s inequality and its applications. Applied Mathematics and Computation, Vol. 151, Issue 3, 2004, p. 887.

Search CrossRef
He Ji-huan, Tang Jie Rebuild of King Fang 40 BC musical scales by He’s inequality. Applied Mathematics and Computation, Vol. 165, Issue 2, 2005, p. 909-914.

Search CrossRef

About this article

Received

24 May 2015

Accepted

11 August 2015

Published

30 September 2015

SUBJECTS

Chaos, nonlinear dynamics and applications

Keywords

strong nonlinear problems

shock response

approximate analytical solution

corrected max-min approach

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.