Next Article in Journal
Effect of MetaFoundation on the Seismic Responses of Liquid Storage Tanks
Next Article in Special Issue
Explicit Solutions to Large Deformation of Cantilever Beams by Improved Homotopy Analysis Method I: Rotation Angle
Previous Article in Journal
The Influence of Sentiments of Economic Agents on Pedestrians and Vehicle Crossings along the US–Mexico Border
Previous Article in Special Issue
An Experimental Study on the Vibration Transmission Characteristics of Wrist Exposure to Hand Transmitted Vibration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Explicit Solution to Large Deformation of Cantilever Beam by Improved Homotopy Analysis Method II: Vertical and Horizontal Displacements

School of Mechanical Engineering, Hebei University of Technology, Tianjin 300401, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(5), 2513; https://doi.org/10.3390/app12052513
Submission received: 12 December 2021 / Revised: 5 January 2022 / Accepted: 6 January 2022 / Published: 28 February 2022

Abstract

:
Explicit solutions to vertical and horizontal displacements are derived for large deformation of a cantilever beam under point load at the free end by an improved homotopy analysis method (IHAM). Quadratic and cubic nonlinear differential equations are adopted to construct more proficient nonlinear equations for vertical and horizontal displacements respectively combined with their currently available nonlinear displacement equations. Higher-order nonlinear iterative homotopy equations are established to solve the vertical and horizontal displacements by combining simultaneous equations of the constructed nonlinear equations and the auxiliary linear equations. The convergence range of vertical displacement is extended by the homotopy-Páde approximation. The explicit solutions to the vertical and horizontal displacements are in favorable agreements with the respective exact solutions. The convergence ranges for a relative error of 1% by the improved homotopy analysis method for vertical and horizontal displacements increases by 60% and 7%, respectively. These explicit formulas are helpful in practical engineering design for very slender structures, such as high-rise buildings and long bridges.

1. Introduction

In view of the shortage of land worldwide, especially in metropolitan areas, high-rise building is a prioritized land-efficient architectural form developed more than 100 years ago. High-rise building is the product of the comprehensive application of modern sciences and engineering technologies. Smooth construction and sustainability of this new structural form depends on more accurate calculation theory and method than those for non-high-rise structures in light of that the extrusion of high-rise building into space could be as high as 800 m plus, and slight artefacts could be exaggerated into intractable unpredictability. Since the 21st century, beam-like structures have achieved overwhelming popularity, such as in high-rise buildings [1,2], long-span bridges [3,4], robots and aerospace; however, analysis methods for this form of structure are still those conventional Euler-Bernoulli laws for large deflection beams, which contain implicit integral to hinder general engineering application. An explicit higher-resolution analysis method is expected to tackle the high nonlinearity in high-rise buildings or long bridges.
The general theory of elastic deflection is called Elastica. To solve a general problem of elastic large deflection, the equivalent load, deflection and longitudinal displacement can be expressed by an analytic expression of elliptic integrals with free-end rotation angle as parameters [5]; indeed, the results have been reported extensively in the classical works of elastic mechanics [6,7,8]. For a cantilever beam subject to point load, the solution to this elastic problem is still an implicit solution with an elliptic integral.
In this study, the high-rise building is simplified as a large deformation of a cantilever beam subject to point load at the free end, for example, of a loading extremity, as shown in Figure 1. Given that the traditional analytical methods for solving nonlinear equations such as perturbation [9], methods introduced in [10,11], Adomian decomposition method [12], and δ -expansion method [13], in addition to the elliptical integrals mentioned above, cannot provide accurate, convergent and simple solutions, we aim to derive an explicit and accurate expression for rotation angle θ b , vertical displacement δ v and horizontal displacement δ h to minimize the potential unpredictability in the large deformation of the long cantilever-like building.
Unlike the perturbation technique, the homotopy analysis method (HAM) is independent of any small/large physical parameters. HAM provides us great freedom to choose equation-type and solution expression of the high order approximation equations. A convergence-control parameter ħ is introduced into the series solutions to guarantee the convergence. HAM distinguishes itself from other analytic approaches to be applied to complicated problems with strong nonlinearity [14,15,16,17,18,19,20,21,22,23,24,25,26,27]. The Föppl–Von Kármán’s plate equations were solved by HAM [28,29]. Especially, the steady-state resonant waves were firstly predicted by the HAM in theory [21,22,30,31], which were then confirmed experimentally [30].
Ji Wang et al. derived an explicit solution to the large deformation of a cantilever beam under point load at the free end with HAM [32]. However, the improficiency of this method is identified to be a low convergence rate due to the step-by-step iteration of solutions of linear differential equations to approximate the exact solutions of complex non-linear differential equations; moreover, it fails to serve explicit solutions of vertical and horizontal displacements for large deformation of the beam.
In our previous paper, “Explicit Solution to Large Deformation Cantilever Beam by Improved Homotopy Analysis Method I: ROTATION ANGLE”, an improved homotopy analysis method (IHAM) was proposed. In IHAM, an auxiliary nonlinear operator ε is introduced to improve the convergence rate and extend the convergence range. The improved homotopy analysis method was used to obtain the explicit solution to the rotation angle for the large deformation of a cantilever beam under point load at the free end.
For small deflection, the horizontal displacement of the cantilever beam is neglected and only the deflection and rotation angle are calculated. For large deflection, horizontal displacement is signified not to be neglected. Based on the explicit solution to the rotation angle obtained in our previous paper, explicit expressions for the deflection and horizontal displacement are derived using the improved homotopy analysis method in this paper.

2. Improved Homotopy Analysis Method

2.1. Problem Description

The bending equations of a uniform cross-section beam with large deformation are [5,6,7,8]
d θ s d s = F E I l x s
d v s d s = sin θ s
d x s d s = cos θ s
subject to boundary conditions
θ 0 = 0 ,   θ L = 0  
v 0 = 0
x 0 = 0
where s is the arc-coordinate of the neutral axis of the beam; x denotes the horizontal coordinate from the fixed end; L represents the length of the beam; F stands for the point load at the free end; E I is the bending stiffness of the beam; θ specifies the rotation angle of cross-section of the beam; and l is the horizontal distance of two ends and is unknown.
Differentiating the equation with respect to s and then using the dimensionless variables ξ = s / L , α = F L 2 / E I , X = x / L , V = v / L , Θ b = θ b / π / 2 , V b = δ v / L and U b = δ h / L , the original equations become
θ ξ + α cos θ ξ = 0
V ξ = sin θ ξ
X ξ = cos θ ξ
subject to boundary conditions
θ 0 = 0 ,   θ 1 = 0
V 0 = 0
X 0 = 0
where the prime denotes the differentiation with respect to ξ .
For infinitesimal deformation, the below linear equations are sufficiently accurate
θ ξ + α = 0
V ξ = θ ξ
X ξ = 1
subject to boundary conditions
θ 0 = 0 ,   θ 1 = 0
V 0 = 0
X 0 = 0
Accordingly, the solutions are
θ ξ = α 2 2 ξ ξ
V ξ = α 6 ξ 2 3 ξ
X ξ = ξ
The rotation angle, vertical and horizontal displacement of cross-section plane at end as denoted by θ b = θ 1 , V b = V 1 , U b = 1 X 1 , are calculated to be
θ b = α 2
V b = α 3
U b = 0
In our previous paper, “Explicit Solution to Large Deformation Cantilever Beam by Improved Homotopy Analysis Method I: ROTATION ANGLE”, the improved homotopy analysis method was proposed to solve Equations (2.3a) and (2.4a), and a series of solutions to the rotation angle are derived as [27]:
θ ξ = θ 0 ξ + n = 1 + θ n ξ
where
θ n ξ = k = 1 + a k ξ k , n = 1 , 2 ,  
In this paper, we apply the improved homotopy analysis method to derive the explicit expressions of deflection and horizontal displacement using Equations (2.3b), (2.4b), (2.3c), and (2.4c) on the basis of the derived explicit expressions of rotation angle in Equations (2.9) and (2.10).

2.2. Fundamentals of Improved Homotopy Analysis Method

For quick reference, the new improved homotopy analysis method (IHAM) is briefed. The IHAM is essentially different from the HAM. Table 1 compares the formulations of IHAM and HAM.
Substantial difference lies in the different construction methods. HAM constructs a higher order nonlinear iterative homotopy differential equation by the selected linear differential equation and the original non-linear differential equation. However, IHAM builds up a higher-order nonlinear iterative homotopy differential equation using the selected linear differential equation, the selected simple non-linear differential equation and the original non-linear differential equation (See Table 1).
Generalized Equation (2.3b) or Equation (2.3c) have the form of
N v ξ = 0 :   v ξ = f w ξ , ξ 0 , 1
where in f w is a derivable function in compliance with the boundary conditions of f w 0 = c 0 , f w 0 = c 1 , 1 2 ! f w 0 = c 2 , 1 3 ! f w 0 = c 3 ; w ξ is the selected polynomial function of an independent variable in Equations (2.9) and (2.10). The boundary condition in Equation (2.4b) or Equation (2.4c) can be expressed as
v 0 = 0
Note that it is unnecessary to assume the existence of any small or large parameters in Equation (2.11). Thus, the proposed approach is general.
Linearizing original Equation (2.11) yields
L v ξ = 0 :   v ξ = c 0 + c 1 w ξ , ξ 0 , 1  
which complies to the boundary condition of
v 0 = 0
We choose the simplest auxiliary nonlinear equations of
N 0 v ξ = 0 :   v ξ = c 0 + c 1 w ξ + c 2 ε 2 w 2 ξ + c 3 ε 1 w 3 ξ ,   ξ 0 , 1
herein, c 0 , c 1 , c 2 and c 3 are undetermined coefficients; and ε 1 0 , 1 , ε 2 0 , 1 , ε 1 2 + ε 2 2 0 . which satisfies the boundary condition of
v 0 = 0
Detailed formulation for the improved homotopy analysis method could be found in our previous paper titled “Explicit Solution to Large Deformation Cantilever Beam by Improved Homotopy Analysis Method I: ROTATION ANGLE”.

3. Vertical Displacement of Cantilever Beam by IHAM

3.1. Zero Order Deformation Equation

The nonlinear boundary value Equations (2.3b) and (2.4b) are then solved by IHAM to derive explicit form for V b . To start with, V ξ is expressed by power series of ξ , i.e.,
V ξ = k = 1 + b k ξ k
where b k is a constant coefficient. Since Equation (2.7b) from linear equation is an appropriate initial guess, we choose
V 0 ξ = α 6 ξ 2 3 ξ
The equivalent nonlinear Equation (2.3b) becomes
V θ + ε 1 3 ! θ 3 = 0
subject to boundary condition
V 0 = 0
where ε 1 is the coefficient, and ε 1 0 , 1 .
When ε 1 = 0 , Equation (3.3) becomes linear Equation (2.5b), so that
V θ = 0
When ε 1 = 1 , Equation (3.3) becomes
V θ + 1 3 ! θ 3 = 0
Equation (3.6) is an approximate expression of the second-order Taylor expansion of sine function in view of that
sin θ θ 1 3 ! θ 3
According to the original Equation (2.3b) and the equivalent Equation (3.3), we define a nonlinear operator as
N v ν ξ ; q , ψ ξ , q , q , ε 1 = d ν ξ ; q d ξ sin q ψ ξ , q
1 q ψ ξ , q + ε 1 3 ! 1 q 3 ψ ξ , q 3
where q 0 , 1 is an embedding parameter; ν ξ ; q , ψ ξ ; q are dependent on ξ and q ; ε 1 is a control parameter to represent the convergent region, and ε 1 0 ,   1 .
When q = 0 , there is
N v ν ξ ; 0 , ψ ξ , 0 , 0 , ε 1 = d ν ξ ; 0 d ξ ψ ξ , 0 + ε 1 3 ! ψ ξ , 0 3
If q = 1 , it yields
N v ν ξ ; 1 , ψ ξ , 1 , 1 , ε 1 = d ν ξ ; 1 d ξ sin ψ ξ , 1
Letting
L v ν ξ ; q = ν ξ ; q ξ
be an auxiliary linear operator, and ħ 1 , 0 denotes a nonzero auxiliary parameter (called convergence-control parameter), a homotopy is then constructed as
H v ν ξ ; q , q : = 1 q L v ν ξ ; q V 0 ξ
q ħ 1 N v ν ξ ; q , ψ ξ , q , q , ε 1
When q = 0 , we have
H v ν ξ ; 0 , 0 : = L v ν ξ ; 0 V 0 ξ
If q = 1 , there is
H v ν ξ ; 1 , 1 : = ħ 1 N v ν ξ ; 1 , 1 , ε 1 = ħ 1 d ν ξ ; 1 d ξ sin ψ ξ , 1
Thus, by enforcing
H v ν ξ ; q , q = 0
we have a family of equations
1 q L v ν ξ ; q V 0 ξ = q ħ 1 N v ν ξ ; q , ψ ξ , q , q , ε 1
subject to the boundary condition of
ν 0 ; q = 0
When q = 0 , we have
ν ξ ; 0 = V 0 ξ , ψ ξ ; 0 = θ 0 ξ
When q = 1 , Equations (3.15) and (3.16) are equivalent to the original Equations (2.3b) and (2.4b), provided that
ν ξ ; 1 = V ξ , ψ ξ ; 1 = θ ξ
Thus, as q increases from 0 to 1, ν ξ ; q evolves smoothly from the specified initial guess V 0 ξ to the solution V ξ of Equations (2.3b) and (2.4b).
With reference to Equation (3.17), by expanding ν ξ ; q into Taylor series of the embedding parameter q , we have
ν ξ ; q = V 0 ξ + n = 1 + V n ξ q n
ψ ξ ; q = θ 0 ξ + n = 1 + θ n ξ q n
where
V n = 1 n ! n ν ξ ; q q n q = 0 ,   θ n = 1 n ! n ψ ξ ; q q n q = 0  
It is worth noting that the zero order deformation Equation (3.15) contains auxiliary parameters ε 1 , ħ 1 . Auxiliary function H 1 ξ ; q , ε 1 has the form of
H 1 ξ ; q , ε 1 = 1 q ψ ξ , q + ε 1 3 ! 1 q 3 ψ ξ , q 3  
Assuming that ħ 1 and ε 1 are properly chosen, so that Equation (3.18) is convergent at q = 1 . Thus, we have
V ξ = V 0 ξ + n = 1 + V n ξ
θ ξ = θ 0 ξ + n = 1 + θ n ξ  
The governing equations of V n ξ can be deduced from the 0th-order deformation Equations (3.15) and (3.16).

3.2. High Order Deformation Equations

Substituting Equation (3.19) into Equation (3.15) and differentiating Equation (3.15) n times with respect to the embedding parameter q , then dividing by n ! and setting q = 0 , we have the n th-order deformation equation as below,
L v V n ξ χ n V n 1 ξ = ħ 1 R n v V 0 , V 1 , , V n 1 ; θ 0 , θ 1 , , θ n 1
Equation (3.21) is subject to boundary conditions of
V n 0 = 0
where
R n v = 1 n 1 ! n 1 N v ν ξ ; q , ψ ξ ; q , q , ε 1 q n 1 q = 0
and
χ n = 0 , n 1 1 ,   n > 1
Combining Equations (3.8) and (3.23), we have
R 1 ν = V 0 ξ θ 0 ξ + ε 1 6 θ 0 3 ξ
R 2 v = V 1 ξ θ 1 ξ + ε 1 2 θ 0 2 ξ θ 1 ξ
R 3 v = V 2 ξ θ 2 ξ + ε 1 2 θ 0 2 ξ θ 2 ξ + ε 1 2 θ 0 ξ θ 1 2 ξ
R 4 v = V 3 ξ θ 3 ξ + 1 6 θ 0 3 ξ ε 1 6 θ 0 3 ξ + ε 1 6 θ 1 3 ξ + ε 1 2 θ 0 2 ξ θ 3 ξ + ε 1 θ 0 ξ θ 1 ξ θ 2 ξ
R 5 v = V 4 ξ θ 4 ξ + 1 2 θ 0 2 ξ θ 1 ξ ε 1 2 θ 0 2 ξ θ 1 ξ + ε 1 2 θ 0 2 ξ θ 4 ξ + ε 1 2 θ 0 ξ θ 2 2 ξ + ε 1 2 θ 1 2 ξ θ 2 ξ + ε 1 θ 0 ξ θ 1 ξ θ 3 ξ  
These equations can be solved by symbolic calculation software such as Maple, Matlab, etc. Thereafter, the linear high-order deformation Equations (3.21) and (3.22) can be solved explicitly.
The solutions given by the IHAM contains auxiliary parameters ε 1 and ħ 1 , which control and adjust the convergence region and rate of IHAM solution series.

3.3. Vertical Displacement at Free end of Cantilever Beam

3.3.1. Accurate Solution of Vertical Displacement

From Equation (2.3), the transcendental equation used to solve V b is as follows
V b = 1 2 α E μ E ϕ , μ
where
E μ = 0 π / 2 d t 1 μ 2 sin 2 t
E ϕ , μ = 0 ϕ d t 1 μ 2 sin 2 t
herein E μ is complete elliptic integral of the second kind; E ϕ , μ is the elliptic integral of the second kind; and
ϕ = arcsin 1 2   μ
μ = 1 + sin θ b 2
From Equation (2.8b), the linear solution of the vertical displacement at the free end is
V b = α 3
When the relative error Δ 1 % , the convergent region of the linear solution is
α 0 , 0.29 ,   V b 0 , 0.097 ,   0 δ v 0.097 L

3.3.2. Control Parameter ε 1 0 , 1 for Convergent Region

Combining Equation (2.3b) with Equation (3.3) gives
sin θ = θ ε 1 3 ! θ 3
When θ = π 2 , we have
1 = π 2 ε 1 3 ! π 2 3
Accordingly, the optimal convergence parameter for the full coverage is derived as
ε 1 = 24 π 2 π 3
When n = 2 , ħ = 1 , ε = 8 / π 2 , the 2nd-order approximation of V b by IHAM is
V b = 1 3 α + 2.3159   × 10 2 ħ 1 α 3 + 1.9048 × 10 3 ħ 1 ε 1 α 3 + 9.5238 × 10 3 ħ 1 2 ε 1 α 3 2.0430 × 10 3 ħ 1 ε 1 α 5
For example, the 2nd-order approximation of V b with ħ 1 = 1 , ε 1 = 24 π 2 / π 3 is
V b = 1 3 α 3.1575 × 10 2 α 3 + 1.8053 × 10 3 α 5
When the relative error Δ 1 % , the convergent region of the 2nd-order approximation by IHAM is
α 0 , 1.5 ,   V b 0 , 0.41 ,   0 δ v 0.41 L
Table 2 shows 2nd order numerical solutions V b n = 2 by IHAM with different values ε 1 compared with exact solutions V b . The relative error in Table 2 is calculated by
Δ = V b n = 2 V b V b × 100 %
After two iterations, when ε 1 = 0 , the convergence region of V b is 0 < α 0.5 . When ε 1 = 24 π 2 / π 3 , the convergence region of V b is 0   <   α 1.5 , whence the convergence region is expanded by 200 % compared to the case with ε 1 = 0 .
When 0 < α 0.5 , ε 1 increases from 0 to 1, the series of solutions of vertical displacement converge to the exact solution with relative error to the exact solution is less than 1%, so that the effective region of ε 1 is identified to be R ε 1 = 0 , 1 .
Thus, the optimal control parameters for the convergence region of the improved homotopy analysis method is identified to be ε 1 = 24 π 2 / π 3 .

3.3.3. Control Parameters ħ 1 1 , 0 for Convergence Rate

By IHAM, when ε = 8 / π 2 , ħ = 0.5 , ε 1 = 24 π 2 / π 3 , ħ 1 = 0.5 , the 10th-order approximation of V b is
V b = 1 3 α 3.7639 × 10 2 α 3 + 7.5564 × 10 3 α 5 1.6013 × 10 3 α 7 + 2.9034 × 10 4 α 9 4.1074 × 10 5 α 11 + 4.3984 × 10 6 α 13 3.5126 × 10 7 α 15 + 2.0566 × 10 8 α 17 8.2692 × 10 10 α 19 + 1.7210 × 10 11 α 21
When the relative error Δ 1 % , the convergent regions of the 10th-order approximation by IHAM are
α 0 , 3.3 ,   V b 0 , 0.63 ,   0 δ v 0.63 L
Table 3 compares the 10th-order numerical solutions with the exact solutions. The relative error in Table 3 is calculated by
Δ = V b n = 10 V b V b × 100 %
After n = 10 iterative calculations, when ħ 1 = 1 , the convergence region of V b is 0   <   α 3 . When ħ 1 = 0.5 , the convergence region of V b is 0   <   α 3.2 . When ħ 1 = 0.3 , the convergence region of V b is 0 < α 1.4 . It is seen that when ħ 1 = 0.5 , the convergence region is 128 % larger than that with ħ 1 = 0.3 , while it is 6.7 % larger than that with ħ 1 = 1 (See Table 3).
Accordingly, the optimal control parameters ħ 1 of convergence rate for the 10th-order approximation solution to the vertical displacement is identified to be ħ 1 = 0.5 .

3.3.4. 30th-Order Approximation Solution of Vertical Displacement

When ε = 8 / π 2 , ħ = 0.1 , the 30th-order approximation of θ b by IHAM is
θ b = α 2 4.3804 × 10 2 α 3 + 7.2872 × 10 3 α 5 1.3299 × 10 3 α 7 + 2.2991 × 10 4 α 9 3.5481 × 10 5 α 11 + 4.7492 × 10 6 α 13 5.4234 × 10 7 α 15 + 5.2245 × 10 8 α 17 4.2081 × 10 9 α 19 + 2.8138 × 10 10 α 21 1.5540 × 10 11 α 23 + 7.0725 × 10 13 α 25 2.6561 × 10 14 α 27 + 8.2720 × 10 16 α 29 2.1528 × 10 17 α 31 + 4.7252 × 10 19 α 33 8.8286 × 10 21 α 35 + 1.4162 × 10 22 α 37 1.9643 × 10 24 α 39 + 2.3681 × 10 26 α 41 2.4891 × 10 28 α 43 + 2.2824 × 10 30 α 45 1.8213 × 10 32 α 47 + 1.2569 × 10 34 α 49 7.4143 × 10 37 α 51 + 3.6656 × 10 39 α 53 1.4683 × 10 41 α 55 + 4.4741 × 10 44 α 57 9.1359 × 10 47 α 59 + 9.2009 × 10 50 α 61
When the relative error Δ 1 % , the convergent region of the 30th-order approximation by IHAM is
α 0 , 5.1 ,   Θ b 0 , 0.79 ,   0 θ b 70.6 °
When ε 1 = 0 , ħ 1 = 0.3 , the 30th-order approximation of V b by HAM is
V b = 1 3 α 3.6198 × 10 2 α 3 + 6.5568 × 10 3 α 5 1.2305 × 10 3 α 7 + 2.1160 × 10 4 α 9 3.1714 × 10 5 α 11 + 4.0427 × 10 6 α 13 4.3197 × 10 7 α 15 + 3.8284 × 10 8 α 17 2.7902 × 10 9 α 19 + 1.6610 × 10 10 α 21 8.0436 × 10 12 α 23 + 3.1678 × 10 13 α 25 1.0188 × 10 14 α 27 + 2.6964 × 10 16 α 29 5.9333 × 10 18 α 31 + 1.0971 × 10 19 α 33 1.7218 × 10 21 α 35 + 2.3130 × 10 23 α 37 2.6772 × 10 25 α 39 + 2.6808 × 10 27 α 41 2.3258 × 10 29 α 43 + 1.7456 × 10 31 α 45 1.1275 × 10 33 α 47 + 6.2040 × 10 36 α 49 2.8576 × 10 38 α 51 + 1.0708 × 10 40 α 53 3.0973 × 10 43 α 55 + 6.1662 × 10 46 α 57 6.2456 × 10 49 α 59
Table 4 compares the above numerical solution with the exact solution. The relative error in Table 4 is calculated by
Δ = V b n = 30 V b V b × 100 %
When the relative error Δ 1 % , the convergent regions of the 30th-order approximation by HAM ( ε 1 = 0 ) are
α 0 , 4.6 ,   V b 0 , 0.70 ,   0 δ v 0.70 L
When ε 1 = 24 π 2 / π 3 , ħ 1 = 0.3 , the 30th-order approximation of V b by IHAM is
V b = 1 3 α 3.6198 × 10 2 α 3 + 6.6118 × 10 3 α 5 1.2700 × 10 3 α 7 + 2.2683 × 10 4 α 9 3.5802 × 10 5 α 11 + 4.8726 × 10 6 α 13 5.6396 × 10 7 α 15 + 5.5006 × 10 8 α 17 4.4899 × 10 9 α 19 + 3.0512 × 10 10 α 21 1.7202 × 10 11 α 23 + 8.0374 × 10 13 α 25 3.1168 × 10 14 α 27 + 1.0073 × 10 15 α 29 2.7301 × 10 17 α 31 + 6.2517 × 10 19 α 33 1.2190 × 10 20 α 35 + 2.0389 × 10 22 α 37 2.9432 × 10 24 α 39 + 3.6841 × 10 26 α 41 4.0095 × 10 28 α 43 + 3.7952 × 10 30 α 45 3.1164 × 10 32 α 47 + 2.2056 × 10 34 α 49 1.3297 × 10 36 α 51 + 6.6933 × 10 39 α 53 2.7186 × 10 41 α 55 + 8.3645 × 10 44 α 57 1.7196 × 10 46 α 59 + 1.7421 × 10 49 α 61
When the relative error Δ 1 % , the convergent regions of the 30th-order approximation by IHAM ( ε 1 = 24 π 2 / π 3 ) are
α 0 , 5.1 ,   V b 0 , 0.72 ,   0 δ v 0.72 L
From Table 4, the improved homotopy analysis method extends the convergence region from α 0 , 4.6 by the homotopy analysis method [32] to α 0 , 5.1 , which is an 11 % increase.
From the above comparison, it can be seen that the expansion rate of the convergence region for vertical displacement by IHAM against that by HAM increases with the increasing iterations.

3.3.5. Homotopy-Páde Approximation

By application of the homotopy-Páde acceleration technique [27], the convergence region can be greatly enlarged. For example, when n = 30 , ε 1 = 0 , ħ 1 = 0.3 , the P a d e 11 , 11 homotopy-Páde approximation of V b in Equation (3.44) by HAM can be written as
V b = α 3 f α g α
f α = 1 + 0.31708 α 2 + 6.1434 × 10 2 α 4 + 6.0404 × 10 3 α 6 + 4.1539 × 10 4 α 8 + 5.5063 × 10 6 α 10
g α = 1 + 0.42567 α 2 + 8.7988 × 10 2 α 4 + 1.0914 × 10 2 α 6 + 8.0640 × 10 4 α 8 + 2.8143 × 10 6 α 10
When the relative error Δ 1 % , the convergent regions of the P a d e 11 , 11 homotopy-Páde approximation by HAM are
α 0 , 5 ,   V b 0 , 0.71 ,   0 δ v 0.71 L
When n = 30 , ε 1 = 24 π 2 / π 3 , ħ 1 = 0.3 , the P a d e 15 , 15 homotopy-Páde approximation of V b in Equation (3.47) by IHAM are
V b = α 3 f α g α
f α = 1 + 0.33026 α 2 + 7.1207 × 10 2 α 4 + 8.9035 × 10 3 α 6 + 9.2792 × 10 4 α 8 + 5.3226 × 10 5 α 10 + 2.5321 × 10 6 α 12 + 1.9082 × 10 8 α 14
g α = 1 + 0.43885 α 2 + 9.9028 × 10 2 α 4 + 1.4763 × 10 2 α 6 + 1.5584 × 10 3 α 8 + 1.1572 × 10 4 α 10 + 5.5647 × 10 6 α 12 + 1.3297 × 10 7 α 14
When the relative error Δ 1 % , the convergent regions of the P a d e 15 , 15 homotopy-Páde approximation by IHAM are
α 0 , 8.1 ,   V b 0 , 0.79 ,   0 δ v 0.79 L
Table 5 shows that numerical solutions V b P of homotopy-Páde approximation by IHAM compares with exact solutions V b . The relative error in Table 5 is calculated by
Δ = V b P V b V b × 100 %
From Equations (3.50) and (3.52), the improved homotopy analysis method extends the convergence region from α 0 , 5 to α 0 , 8.1 for the homotopy-Páde approximate solution compared with the homotopy analysis method [32], which is a 62 % increase.
Figure 2 compares the 30th-order iteration solution and the exact solution of the vertical displacement. From Figure 2, it can be observed that
  • the difference between the linear solution and the exact solution is remarkable;
  • the 30th-order improved homotopy analysis solution for ε 1 = 24 ( π 2 ) / π 3 is much closer to the exact solution than does the 30th-order homotopy analysis solution for ε 1 = 0 ;
  • the convergence range and rate have been greatly improved by the expression by homotopy-Páde approximation;
  • the Páde approximation solution P a d e 15 , 15 by IHAM is more approximating to the exact solution than does the Páde approximation solution P a d e 11 , 11 by HAM.

4. Horizontal Displacement of Cantilever Beam by IHAM

4.1. Zero Order Deformation Equation

The nonlinear boundary value Equations (2.3c) and (2.4c) are then solved by IHAM to give explicit expression for U b . To start with, X ξ is expressed by the power series of ξ , so that
X ξ = k = 1 + c k ξ k
where c k is a constant coefficient. Linear Equation (2.7c) is specified as initial guess. Accordingly, we choose
X 0 ξ = ξ
The equivalent nonlinear equation of Equation (2.3c) is
X 1 + ε 2 2 ! θ 2 = 0
subject to boundary condition of
X 0 = 0
where ε 2 is the correction coefficient of the convergent region, and ε 2 0 , 1 .
When ε 2 = 0 , Equation (4.3) becomes the linear Equation (2.5c), so that
X 1 = 0
When ε 2 = 1 , Equation (4.3) becomes
X 1 + 1 2 ! θ 2 = 0
Approximation expression of the second-order Taylor expansion of cosine function can be written as
cos θ 1 1 2 ! θ 2
According to the original Equation (2.3c) and the equivalent Equation (4.3), we define a nonlinear operator as
N u μ ξ ; q , ψ ξ , q , q , ε 2 = d μ ξ ; q d ξ cos q ψ ξ , q + ε 2 2 ! 1 q 2 ψ ξ , q 2
where q 0 , 1 is an embedding parameter; μ ξ ; q , ψ ξ ; q are the functions dependent on ξ and q ; ε 2 is a control parameter for the convergent region, and ε 2 0 , 1 .
When q = 0 , there is
N u μ ξ ; 0 , ψ ξ , 0 , 0 , ε 2 = d μ ξ ; 0 d ξ 1 + ε 2 2 ! ψ ξ , 0 2
If q = 1 , it yields
N u μ ξ ; 1 , ψ ξ , 1 , 1 , ε 2 = d μ ξ ; 1 d ξ sin ψ ξ , 1
Specifying
L u μ ξ ; q = μ ξ ; q ξ
as an auxiliary linear operator, and denoting ħ 2 1 , 0 as a nonzero auxiliary parameter (convergence-control parameter), a homotopy is constructed as
H u μ ξ ; q , q : = 1 q L u μ ξ ; q X 0 ξ
q ħ 2 N u μ ξ ; q , ψ ξ , q , q , ε 2
Accordingly, when q = 0 , we have
H u μ ξ ; 0 , 0 : = L u μ ξ ; 0 X 0 ξ
If q = 1 , it becomes
H u μ ξ ; 1 , 1 : = ħ 2 N u μ ξ ; 1 , 1 , ε 2 = ħ 2 d μ ξ ; 1 d ξ cos ψ ξ , 1
Thus, by enforcing
H u μ ξ ; q , q = 0
we have a family of equations
1 q L u μ ξ ; q X 0 ξ = q ħ 2 N u μ ξ ; q , ψ ξ , q , q , ε 2
subject to boundary conditions
μ 0 ; q = 0
When q = 0 , we have
μ ξ ; 0 = X 0 ξ ,   ψ ξ ; 0 = θ 0 ξ
When q = 1 , Equations (4.16) and (4.17) are equivalent to the original Equations (2.3c) and (2.4c), provided that
μ ξ ; 1 = X ξ ,   ψ ξ ; 1 = θ ξ
Thus, as q increases from 0 to 1, μ ξ ; q evolves from the known initial guess solution X 0 ξ to the exact solution X ξ of Equations (2.3c) and (2.4c).
With reference to Equation (4.10), by expanding μ ξ ; q into Taylor series of the embedding parameter q , we have
μ ξ ; q = X 0 ξ + n = 1 + X n ξ q n
ψ ξ ; q = θ 0 ξ + n = 1 + θ n ξ q n
where
X n = 1 n ! n μ ξ ; q q n q = 0 ,   θ n = 1 n ! n ψ ξ ; q q n q = 0
It is worth noting that the zero order deformation Equation (4.16) contains auxiliary parameters ε 2 , ħ 2 . And the auxiliary function H 2 ξ ; q , ε 2 has the form of
H 2 ξ ; q , ε 2 = ε 2 2 ! 1 q 2 ψ ξ , q 2
Assuming that ħ 2 and ε 2 are properly chosen, so that Equation (4.19) is convergent at q = 1 , thus, we have
X ξ = X 0 ξ + n = 1 + X n ξ
The governing equations X n ξ can be deduced from 0th-order deformation Equations (4.16) and (4.17).

4.2. High Order Deformation Equations

Substituting Equation (4.20) into Equation (4.16) and differentiating (4.16) n times with respect to the embedding parameter q , then dividing by n ! and setting q = 0 , we have the n th-order deformation equation as below
L u X n ξ χ n X n 1 ξ = ħ 2 R n u X 0 , X 1 , , X n 1 ; θ 0 , θ 1 , , θ n 1
Equation (4.24) is subject to boundary condition of
X n 0 = 0
where
R n u = 1 n 1 ! n 1 N u μ ξ ; q , ψ ξ ; q , q , ε 2 q n 1 q = 0
and
χ n = 0 , n 1 1 ,   n > 1
Combining Equations (4.8) and (4.26), we have
R 1 u = X 0 ξ 1 + ε 2 2 θ 0 2 ξ
R 2 u = X 1 ξ + ε 2 θ 0 ξ θ 1 ξ
R 3 u = X 2 ξ + 1 2 θ 0 2 ξ ε 2 2 θ 0 2 ξ + ε 2 2 θ 1 2 ξ + ε 2 θ 0 ξ θ 2 ξ
R 4 u = X 3 ξ + θ 0 ξ θ 1 ξ ε 2 θ 0 ξ θ 1 ξ + ε 2 θ 1 ξ θ 2 ξ + ε 2 θ 0 ξ θ 3 ξ
R 5 u = X 4 ξ + θ 0 ξ θ 2 ξ + 1 2 θ 1 2 ξ 1 24 θ 0 4 ξ ε 2 θ 0 ξ θ 2 ξ + ε 2 θ 1 ξ θ 3 ξ + ε 2 θ 0 ξ θ 4 ξ ε 2 2 θ 1 2 ξ + ε 2 2 θ 2 2 ξ
The right-hand side term R n u of the above equations can be derived by symbolic calculation software such as Maple, Matlab, etc. Thereafter, explicit form of solution to the linear high-order deformation Equations (4.24) and (4.25) have been derived.

4.3. Horizontal Displacement at Free End of Cantilever Beam

4.3.1. Accurate Solution to Horizontal Displacement

Accordingly, in Equation (2.3c), the transcendental equation to solve U b is as follows
U b = 1 2 sin θ b α
With reference to Equation (2.8c), the linear solution of the horizontal displacement at the free end is
U b = 0

4.3.2. Explicit Solution to Horizontal Displacement by IHAM

The solutions by IHAM contains two auxiliary parameters ε 2 and ħ 2 , which control and adjust the convergence region and rate of IHAM solution series.
For example, the 10th-order approximation of U b with n = 10 , ε = 8 / π 2 , ħ = 0.5 , ε 2 = 8 / π 2 , ħ 2 = 0.4 is derived as
U b = 6.6128 × 10 2 α 2 1.2132 × 10 2 α 4 + 2.6029 × 10 3 α 6 5.0108 × 10 4 α 8 + 7.4259 × 10 5 α 10 8.1166 × 10 6 α 12 + 6.5090 × 10 7 α 14 3.7900 × 10 8 α 16 + 1.5743 × 10 9 α 18 3.4811 × 10 11 α 20
When the relative error Δ 1 % , the convergent regions of the 10th-order approximation by IHAM are
α 0 , 2.8 ,   U b 0 , 0.24 ,   0 δ h 0.24 L
The 20th-order approximation of U b with n = 20 , ε = 8 / π 2 , ħ = 0.3 , ε 2 = 8 / π 2 , ħ 2 = 0.4 is expressed as
U b = 6.6663 × 10 2 α 2 1.2647 × 10 2 α 4 + 3.0553 × 10 3 α 6 7.8913 × 10 4 α 8 + 1.9565 × 10 4 α 10 4.2939 × 10 5 α 12 + 7.8795 × 10 6 α 14 1.1689 × 10 6 α 16 + 1.3828 × 10 7 α 18 1.3058 × 10 8 α 20 + 9.9286 × 10 10 α 22 6.1404 × 10 11 α 24 + 3.1142 × 10 12 α 26 1.3002 × 10 13 α 28 + 4.4587 × 10 15 α 30 1.2419 × 10 16 α 32 + 2.7429 × 10 18 α 34 4.5778 × 10 20 α 36 + 5.1643 × 10 22 α 38 2.9055 × 10 24 α 40
When the relative error Δ 1 % , the convergent region of the 20th-order approximation by IHAM is
α 0 , 3.4 ,   U b 0 , 0.29 ,   0 δ h 0.29 L
The 30th-order approximation of U b with n = 30 , ε = 8 / π 2 , ħ = 0.1 , ε 2 = 8 / π 2 , ħ 2 = 0.3 is
U b = 6.6665 × 10 2 α 2 1.1891 × 10 2 α 4 + 2.4037 × 10 3 α 6 4.6144 × 10 4 α 8 + 7.8466 × 10 5 α 10 1.1452 × 10 5 α 12 + 1.4110 × 10 6 α 14 1.4525 × 10 7 α 16 + 1.2395 × 10 8 α 18 8.7140 × 10 10 α 20 + 5.0250 × 10 11 α 22 2.3731 × 10 12 α 24 + 9.1975 × 10 14 α 26 2.9421 × 10 15 α 28 + 7.8331 × 10 17 α 30 1.7531 × 10 18 α 32 + 3.3309 × 10 20 α 34 5.4217 × 10 22 α 36 + 7.6171 × 10 24 α 38 9.2885 × 10 26 α 40 + 9.8639 × 10 28 α 42 9.1295 × 10 30 α 44 + 7.3481 × 10 32 α 46 5.1117 × 10 34 α 48 + 3.0383 × 10 36 α 50 1.5131 × 10 38 α 52 + 6.1056 × 10 41 α 54 1.8743 × 10 43 α 56 + 3.8564 × 10 46 α 58 3.9126 × 10 49 α 60
When the relative error Δ 1 % , the convergent regions of the 30th-order approximation by IHAM are
α 0 , 4.6 ,   U b 0 , 0.37 ,   0 δ h 0.37 L
Table 6 compares the above numerical solutions with the exact solutions. The relative error in Table 6 is calculated by
Δ = U b n = 30 U b U b × 100 %
From Table 6, IHAM extends the convergence regions from α 0 , 4.2 by HAM to α 0 , 4.6 , which is 7.1 % increase.
Figure 3 shows the relationship between horizontal displacement and concentrated force at the free end. From Figure 3, it is observed that
  • the difference between the linear solution and the exact solution is remarkable;
  • from n = 10 result of Equation (4.32), n = 20 result of Equation (4.34) and n = 30 result of Equation (4.37), the convergence regions of horizontal displacement increases with increasing iterations.

5. Discussions and Conclusions

5.1. Error Analysis of Vertical Displacement

The influence of the control parameters ε 1 and ħ 1 on the error of vertical displacement after n = 5 iterations with α = 2.0 is then discussed.
After five iterations, the improved homotopy analysis method is adopted to solve the Equations (3.21) and (3.22) to yield the solution to the original Equations (2.3b) and (2.4b) as below
V b = 0.33333 α + 0.096738 ħ 1 α 3 + 0.10197 ħ 1 2 α 3 + 0.052108 ħ 1 3 α 3 + 0.011580 ħ 1 4 α 3 + 0.028571 ε 1 ħ 1 α 3 + 0.085714 ε 1 ħ 1 2 α 3 + 0.095238 ε 1 ħ 1 3 α 3 + 0.047619 ε 1 ħ 1 4 α 3 + 9.5238 × 10 3 ε 1 ħ 1 5 α 3 6.1695 × 10 3 ħ 1 α 5 3.3120 × 10 3 ħ 1 2 α 5 8.2799 × 10 4 ħ 1 3 α 5 5.8320 × 10 3 ε 1 ħ 1 α 5 8.1553 × 10 3 ε 1 ħ 1 2 α 5 4.5967 × 10 3 ε 1 ħ 1 3 α 5 1.0215 × 10 3 ε 1 ħ 1 4 α 5 + 2.5955 × 10 4 ħ 1 α 7 + 7.4156 × 10 5 ħ 1 2 α 7 + 6.8257 × 10 4 ε 1 ħ 1 α 7 + 4.3913 × 10 4 ε 1 ħ 1 2 α 7 + 1.0978 × 10 4 ε 1 ħ 1 3 α 7 7.4415 × 10 6 ħ 1 α 9 4.2843 × 10 5 ε 1 ħ 1 α 9 1.2241 × 10 5 ε 1 ħ 1 2 α 9 + 1.4103 × 10 6 ε 1 ħ 1 α 11
When α = 2.0 , Equation (5.1) becomes
V b = 0.66667 + 0.60589 ħ 1 + 0.71928 ħ 1 2 + 0.39037 ħ 1 3 + 0.092637 ħ 1 4 + 0.11027 ε 1 ħ 1 + 0.47469 ε 1 ħ 1 2 + 0.62886 ε 1 ħ 1 3 + 0.34826 ε 1 ħ 1 4 + 0.076190 ε 1 ħ 1 5
When α = 2.0 , the exact solution is V b = 0.49346 . When the iteration is five times, the relative error between the result by IHAM and the exact solution can be expressed as
Δ V = V b n = 5 V b V b × 100 %
When ħ 1 = 0.68 , Equation (5.2) becomes
V b = 0.48432 + 0.010163 ε 1
According to Equations (5.3) and (5.4), the relative error curve of control parameter Δ V ε 1 is a broken line. When the relative error is less than Δ V 1 % , the convergence interval is R ε 1 = 0.4 , 1 , which is the effective region of ε 1 .
Figure 4 shows the relative errors curve Δ V ε 1 of Equation (2.3b) when α = 2.0 , ħ 1 = 0.68 . From Figure 4, when ε 1 = 24 π 2 / π 3 , the relative error is the minimum.
When ε 1 = 24 π 2 / π 3 , Equation (5.2) becomes
V b = 0.66667 + 0.70333 ħ 1 + 1.1387 ħ 1 2 + 0.94605 ħ 1 3 + 0.40038 ħ 1 4 + 0.067325 ħ 1 5
According to Equations (5.3) and (5.5), the relative error curve of control parameter Δ V ħ 1 is a quintic power function. When the relative error is less than Δ V 1 % , the convergence interval is R ħ 1 = 1 , 0.58 , which is the effective region of ħ 1 .
Figure 5 gives the relative errors curve Δ V ħ 1 of Equation (2.3b) when α = 2.0 , ε 1 = 24 π 2 / π 3 . From Figure 5, when ħ 1 = 0.68 , the relative error is the minimum.

5.2. Error Analysis of Horizontal Displacement

The influence of the control parameter ε 2 and ħ 2 on the error of horizontal displacement after 5 iterations when α = 2.0 is then discussed.
After five iterations, the improved homotopy analysis method is adopted to solve the Equations (4.24) and (4.25) to derive the solution to the original Equations (2.3c) and (2.4c). Thus,
U b = 0.20000 ħ 2 α 2 0.20000 ħ 2 2 α 2 0.066667 ħ 2 3 α 2 0.13333 ε 2 ħ 2 α 2 0.46667 ε 2 ħ 2 2 α 2 0.60000 ε 2 ħ 2 3 α 2 0.33333 ε 2 ħ 2 4 α 2 0.066667 ε 2 ħ 2 5 α 2 + 0.010871 ħ 2 α 4 + 3.7741 × 10 3 ħ 2 2 α 4 + 0.017903 ε 2 ħ 2 α 4 + 0.027904 ε 2 ħ 2 2 α 4 + 0.017361 ε 2 ħ 2 3 α 4 + 3.7741 × 10 3 ε 2 ħ 2 4 α 4 2.6973 × 10 4 ħ 2 α 6 1.6042 × 10 3 ε 2 ħ 2 α 6 1.1329 × 10 3 ε 2 ħ 2 2 α 6 2.6973 × 10 4 ε 2 ħ 2 3 α 6 + 8.2194 × 10 5 ε 2 ħ 2 α 8 + 2.1630 × 10 5 ε 2 ħ 2 2 α 8 1.8596 × 10 6 ε 2 ħ 2 α 10
When α = 2.0 , Equation (5.6) becomes
U b = 0.64333 ħ 2 0.73961 ħ 2 2 0.26667 ħ 2 3 0.33042 ε 2 h 2 1.4872 ε 2 ħ 2 2 2.1395 ε 2 ħ 2 3 1.2729 ε 2 ħ 2 4 0.26667 ε 2 ħ 2 5
When α = 2.0 , the exact solution is U b = 0.16064 . After 5 time iterations, the relative error between the results by IHAM method and the exact solution can be written as
Δ U = U b n = 5 U b U b × 100 %
When ħ 2 = 0.465 , Equation (5.7) becomes
U b = 0.16603 0.0065196 ε 2
According to Equations (5.8) and (5.9), the relative error curve of control parameter Δ U ε 2 is a broken line. When the relative error is less than Δ U 1 % , the convergence region is R ε 2 = 0.58 , 1 , which is the effective region of ε 2 .
Figure 6 shows the relative errors curve Δ U ε 2 of Equation (2.3c) when α = 2.0 , ħ 2 = 0.465 . From Figure 6, when ε 2 = 8 / π 2 , the relative error is the minimum.
When ε 2 = 8 / π 2 , Equation (5.7) becomes
U b = 0.91116 ħ 2 1.9451 ħ 2 2 2.0009 ħ 2 3 1.0318 ħ 2 4 0.21615 ħ 2 5
According to Equations (5.8) and (5.10), the relative error curve of control parameter Δ U ħ 2 is a quintic polynomial curve. When the relative error is less than Δ U 1 % , the convergence interval is R ħ 2 = 0.75 , 0.42 , which is the effective region of ħ 2 .
Figure 7 depicts the relative errors curve Δ U ħ 2 of Equation (2.3c) when α = 2.0 , ε 2 = 8 / π 2 . From Figure 7, when ħ 2 = 0.465 or ħ 2 = 0.66 , the relative error is the minimum.

5.3. Conclusions

In this paper, an improved homotopy analysis method (IHAM) is applied to derive complete expressions of vertical and horizontal displacements of a cantilever beam under point load at the free end. In fact, the proposed method is a semi-analytical method since many numerical algorithms are used.
The linear solution is only valid for small α . For large α , the nonlinear effect of the equation cannot be neglected. A family of the explicit solutions of the rotation angle, the vertical and horizontal displacements of a cantilever beam under a point load at the free end are obtained by IHAM, which are summarized in Table A1 for 30th-order approximation. The current explicit solutions with different parameters and solution schemes match well with the exact solutions from elliptical integrals; however, our solutions for the current problem by the proposed IHAM, especially Equations (3.36) and (4.31), are explicit and simple compared with the exact solution from the elliptical integrals which requires solving a transcendental equation. The current solution is easy to calculate with the explicit polynomial expressions, which should be, thus, prioritized for practical engineering applications with minimal requirements on the calculation and computation.
By solving the exact bending equation of a cantilever beam, we demonstrate the superior proficiency of IHAM over a class of existing equations in solid mechanics. The improved homotopy analysis method has great superiority compared over the traditional homotopy analysis method in view of the very much increased rate and precision of convergence as well as the substantially enlarged convergence range.
The convergence ranges of the rotation angle and vertical displacement are extended by the homotopy-Páde approximation, whereas it does not avail horizontal displacement.
In summary, the 30th-order improved homotopy analysis solutions and exact solutions of the angular, vertical and horizontal displacements of large deflection cantilever beams by Equations (3.42), (3.47) and (4.35) are compared in the following Appendix A table. It can be concluded from this study that the proposed IHAM can be applied to solve a general strong nonlinear problem in structural analysis.

Author Contributions

Conceptualization and methodogoly, Y.L. and X.L.; validation and data curation, C.X. and S.H.; formal analysis, Y.L.; investigation, X.L.; writing—original draft preparation, C.X.; writing—review and editing, X.L.; visualization, S.H.; supervision, Y.L.; project administration, Y.L.; funding acquisition, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 11972145) and Tianjin Excellent Science and Technology Commissioner Support Project (Grant No. 20YDTPJC0080).

Institutional Review Board Statement

No application.

Informed Consent Statement

No application.

Data Availability Statement

No application.

Acknowledgments

This work was funded by the National Natural Science Foundation of China (Grant No. 11972145) and Tianjin Excellent Science and Technology Commissioner Support Project (Grant No. 20YDTPJC0080).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. 30th-order numerical solutions Θ b n = 30 , V b n = 30 , U b n = 30 by IHAM compared with exact solutions Θ b , V b and U b .
Table A1. 30th-order numerical solutions Θ b n = 30 , V b n = 30 , U b n = 30 by IHAM compared with exact solutions Θ b , V b and U b .
α Rotation Angle
Θ b
Vertical Displacement
V b
Horizontal Displacement
U b
Exact Solutions n = 30
ε = 8 / π 2
h = 0.1
Exact Solutions n = 30
ε 1 = 24 π 2 / π 3
h 1 = 0.3
Exact Solutions n = 30
ε 2 = 8 / π 2
h 2 = 0.3
0000 00
0.30.0947190.0947510.0989910.0990380.00589940.0059052
0.60.185090.185300.192350.192660.0224880.022563
0.90.268010.268530.276210.276990.0470340.047300
1.20.342060.342910.349010.350300.0764000.076933
1.50.407140.408250.410980.412710.107940.10872
1.80.463940.465200.463260.465340.139810.14076
2.10.513420.514770.507320.509660.170850.17190
2.40.556610.557980.544550.547090.200460.20153
2.70.594440.595800.576180.578880.228350.22939
3.00.627720.629050.603250.606050.254420.25540
3.30.657140.658430.626580.629450.278700.27961
3.60.683270.684520.646840.649750.301280.30210
3.90.706590.707740.664550.667460.322280.32310
4.20.727490.728210.680140.682820.341800.34302
4.50.746300.746640.693970.695710.359990.36042
4.80.763290.706290.707240.37695
5.10.778700.717350.723490.39279
5.40.792720.727310.40762
5.70.805520.736350.42153
6.00.817230.744570.43459
111111

References

  1. He, G.K. High-rise building structure facing the 21st century. Build. Sci. 2002, 11, 1–4. (In Chinese) [Google Scholar]
  2. Hu, S.D.; Fang, E.H. Analysis of the world’s tallest 100 buildings. Constr. Technol. 2004, 35, 540–542. (In Chinese) [Google Scholar]
  3. Xiang, H.F. Prospect of world bridge engineering in the 21st century. J. Civ. Eng. 2000, 33, 1–6. (In Chinese) [Google Scholar]
  4. Qin, R. Long-Span Bridge Structure; Science Publishing House: Beijing, China, 2008. [Google Scholar]
  5. Gere, J.M.; Timoshenko, S.P. Mechanics of Materials; PWS Publishing Company: Boston, MA, USA, 1997. [Google Scholar]
  6. Love, A.E.H. The Mathematics Theory of Elasticicy, 4th ed.; New York Dover Publications: Mineola, NY, USA, 1944. [Google Scholar]
  7. Frisch-Fay, R. Flexible Bars; Butter Worths: London, UK, 1962. [Google Scholar]
  8. Yu, T.X.; Zhang, L.Z. Theory of Plastic Bending and Its Applications; Science Publishing House: Beijing, China, 1992. [Google Scholar]
  9. Nayfeh, A.H. Perturbation Methods; Wiley: New York, NY, USA, 2000. [Google Scholar]
  10. Novozhilov, V.V. Foundations of the Nonlinear Theory of Elasticity; Graylock Press: New York, NY, USA, 1953. [Google Scholar]
  11. Luo, A.C.J. Nonlinear Deformable-Body Dynamics; Higher Education Press: Beijing, China; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  12. Adomian, G. Nonlinear stochastic differential equations. J. Math. Anal. Appl. 1976, 55, 441–452. [Google Scholar] [CrossRef] [Green Version]
  13. Karmishin, A.V.; Zhukov, A.T.; Kolosov, V.G. Methods of Dynamics Calculation and Testing for Thin-Walled Structures; Mashinostroyenie: Moscow, Russia, 1990. (In Russian) [Google Scholar]
  14. Abbasbandy, S. The application of homotopy analysis method to nonlinear equations arising in heat transfer. Phys. Lett. A 2006, 360, 109–113. [Google Scholar] [CrossRef]
  15. Liang, S.; Jeffrey, D.J. Approximate solutions to a parameterized sixth order boundary value problem. Comput. Math. Appl. 2010, 59, 247–253. [Google Scholar] [CrossRef] [Green Version]
  16. Ghotbi, A.R.; Omidvar, M.; Barari, A. Infiltration in unsaturated soils—An analytical approach. Comput. Geotech. 2011, 38, 777–782. [Google Scholar] [CrossRef]
  17. Nassar, C.J.; Revelli, J.F.; Bowman, R.J. Application of the homotopy analysis method to the Poisson Boltzmann equation for semiconductor devices. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 2501–2512. [Google Scholar] [CrossRef]
  18. Aureli, M. A framework for iterative analysis of non-classically damped dynamical systems. J. Sound Vib. 2014, 333, 6688–6705. [Google Scholar] [CrossRef]
  19. Sardanyés, J.; Rodrigues, C.; Januário, C.; Martins, N.; Gil-Gómez, G.; Duarte, J. Activation of effector immune cells promotes tumor stochastic extinction: A homotopy analysis approach. Appl. Math. Comput. 2015, 252, 484–495. [Google Scholar] [CrossRef] [Green Version]
  20. Zou, K.; Nagarajaiah, S. An analytical method for analyzing symmetry-breaking bifurcation and period-doubling bifurcation. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 780–792. [Google Scholar] [CrossRef] [Green Version]
  21. Xu, D.; Lin, Z.; Liao, S.; Stiassnie, M. On the steady-state fully resonant progressive waves in water of finite depth. J. Fluid Mech. 2012, 710, 379–418. [Google Scholar] [CrossRef] [Green Version]
  22. Liu, Z.; Liao, S.J. Steady-state resonance of mulendle wave interactions in deep water. J. Fluid Mech. 2014, 742, 664–700. [Google Scholar] [CrossRef]
  23. Cheng, J.; Zhu, S.P.; Liao, S.J. An explicit series approximation to the optimal exercise boundary of American put options. Commun. Nonlin. Sci. Numer. Simul. 2010, 15, 1148–1158. [Google Scholar] [CrossRef]
  24. Liao, S.; Tan, Y. A general approach to obtain series solutions of nonlinear differential equations. Stud. Appl. Math. 2007, 119, 297–354. [Google Scholar] [CrossRef]
  25. Liao, S.; Zhao, Y. On the method of directly defining inverse mapping for nonlinear differential equations. Numer. Algor. 2016, 72, 989–1020. [Google Scholar] [CrossRef] [Green Version]
  26. Liao, S.J. Beyond Perturbation: Introduction to the Homotopy Analysis Method; Chapman & Hall/CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  27. Sajid, M.; Hayat, T. Comparison of HAM and HPM methods for nonlinear heat conduction and convection equations. Nonlinear Anal. Real World Appl. 2008, 9, 2296–2301. [Google Scholar] [CrossRef]
  28. Van Gorder, R.A. Analytical method for the construction of solutions to the Föppl–Von Kármán equations governing deflections of a thin flat plate. Int. J. Non-Linear Mech. 2012, 47, 1–6. [Google Scholar] [CrossRef]
  29. Zhong, X.X.; Liao, S.J. Analytic Solutions of Von Kármán Plate under Arbitrary Uniform Pressure—Equations in Differential Form. Stud. Appl. Math. 2017, 138, 371–400. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, Z.; Xu, D.L.; Li, J.; Peng, T.; Alsaedi, A.; Liao, S.J. On the existence of steady-state resonant waves in experiments. J. Fluid Mech. 2015, 763, 1–23. [Google Scholar] [CrossRef] [Green Version]
  31. Liao, S.; Xu, D.; Stiassnie, M. On the steady-state nearly resonant waves. J. Fluid Mech. 2016, 794, 175–199. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, J.; Chen, J.-K.; Liao, S. An explicit solution of the large deformation of a cantilever beam under point load at the free end. J. Comput. Appl. Math. 2008, 212, 320–330. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Large deformation of beam under point load at the free end.
Figure 1. Large deformation of beam under point load at the free end.
Applsci 12 02513 g001
Figure 2. Comparison of vertical displacement at the free end V b with n = 30 .
Figure 2. Comparison of vertical displacement at the free end V b with n = 30 .
Applsci 12 02513 g002
Figure 3. Horizontal displacement at free end.
Figure 3. Horizontal displacement at free end.
Applsci 12 02513 g003
Figure 4. Relative errors curve Δ V ε 1 of Equation (2.3b) when α = 2.0 , ħ 1 = 0.68 .
Figure 4. Relative errors curve Δ V ε 1 of Equation (2.3b) when α = 2.0 , ħ 1 = 0.68 .
Applsci 12 02513 g004
Figure 5. Relative errors curve Δ V ħ 1 of Equation (2.3b) when α = 2.0 , ε 1 = 24 π 2 / π 3 .
Figure 5. Relative errors curve Δ V ħ 1 of Equation (2.3b) when α = 2.0 , ε 1 = 24 π 2 / π 3 .
Applsci 12 02513 g005
Figure 6. Relative errors curve Δ U ε 2 of Equation (2.3c) when α = 2.0 ,   ħ 2 = 0.465 .
Figure 6. Relative errors curve Δ U ε 2 of Equation (2.3c) when α = 2.0 ,   ħ 2 = 0.465 .
Applsci 12 02513 g006
Figure 7. Relative errors curve Δ U ħ 2 of Equation (2.3c) when α = 2.0 , ε 2 = 8 / π 2 .
Figure 7. Relative errors curve Δ U ħ 2 of Equation (2.3c) when α = 2.0 , ε 2 = 8 / π 2 .
Applsci 12 02513 g007
Table 1. Comparison of improved homotopy analysis method against homotopy analysis method.
Table 1. Comparison of improved homotopy analysis method against homotopy analysis method.
Homotopy Analysis Method
(HAM)
Improved Homotopy Analysis Method (IHAM)
Original N ξ = 0 and
Auxiliary L ξ = 0
Higher order H ξ ; q = 0
① Original N ξ = 0 and
Auxiliary N 0 ξ = 0
New   N ¯ ξ ; q = 0
② New   N ¯ ξ ; q = 0 and
Auxiliary L ξ = 0
Higher order H ¯ ξ ; q = 0
Notes: N nonlinear operator, L linear operator, N 0 auxiliary nonlinear operator, N ¯ new nonlinear operator, H homotopy operator, ξ independent variable, q embedding parameter.
Table 2. 2st-order numerical solutions V b n = 2 by IHAM ( ħ 1 = 1 ) with different values ε 1 compared with exact solutions V b .
Table 2. 2st-order numerical solutions V b n = 2 by IHAM ( ħ 1 = 1 ) with different values ε 1 compared with exact solutions V b .
a Exact Solution ε 1 = 0
(HAM)
ε 1 = 24 π 2 / π 3
(IHAM)
ε 1 = 1
(IHAM)
V b V b Δ V b Δ V b Δ
00000000
0.10.0332950.0333100.0446%0.0333020.0194%0.0333010.0161%
0.20.0663640.0664810.176%0.0664140.0756%0.0664050.0623%
0.30.0989910.0993750.388%0.0991520.163%0.0991230.133%
0.40.130980.131850.669%0.131330.272%0.131260.219%
0.50.162140.163771.00%0.162780.390%0.162650.309%
0.60.192350.195001.38%0.193320.504%0.193100.389%
0.70.221480.225391.76%0.222810.599%0.222470.445%
0.80.249450.254812.15%0.251090.658%0.250600.461%
0.90.276210.283122.50%0.278050.666%0.277380.425%
1.00.301720.310172.80%0.303560.611%0.302690.322%
1.10.325980.335843.02%0.327550.480%0.326460.145%
1.20.349010.359983.14%0.349930.264%0.348610.115%
1.30.370820.382453.14%0.370670.0427%0.369110.461%
1.40.391470.403122.97%0.389730.442%0.387970.892%
1.50.410980.421842.64%0.407140.933%0.405211.40%
1.60.429410.438472.11%0.422931.51%0.420891.99%
1.70.446820.452891.36%0.437172.16%0.435102.62%
1.80.463260.464940.361%0.449972.87%0.448003.30%
1.90.478790.474480.899%0.461463.62%0.459753.98%
20.493460.481392.44%0.471844.38%0.470584.64%
1101010
Table 3. 10th-order numerical solutions V b n = 10 by IHAM with different values of ħ 1 ( ε 1 = 24 π 2 / π 3 ) compared with exact solutions V b .
Table 3. 10th-order numerical solutions V b n = 10 by IHAM with different values of ħ 1 ( ε 1 = 24 π 2 / π 3 ) compared with exact solutions V b .
α Exact Solution h 1 = 1 h 1 = 0.5 h 1 = 0.3
V b V b Δ V b Δ V b Δ
00000000
0.20.0663640.0663650.000951%0.0663680.00522%0.0663850.0316%
0.40.130980.130980.00286%0.131000.0183%0.131130.120%
0.60.192350.192360.00418%0.192410.0339%0.192830.252%
0.80.249450.249460.00464%0.249570.0480%0.250470.410%
1.00.301720.301740.00486%0.301900.0602%0.303480.583%
1.20.349010.349020.00369%0.349260.0713%0.351680.765%
1.40.391470.391450.00312%0.391780.0803%0.395190.950%
1.60.429410.429330.0201%0.429770.0829%0.434291.14%
1.80.463260.463050.0457%0.463610.0750%0.469361.32%
2.00.493460.493130.0665%0.4937400572%0.500821.49%
2.20.520420.520140.0539%0.520610.0371%0.529081.66%
2.40.544550.544720.0319%0.544710.0296%0.554571.84%
2.60.566190.567450.222%0.566460.0475%0.577672.03%
2.80.585670.588460.476%0.586100.0741%0.598642.21%
30.603250.607090.636%0.603400.0243%0.617452.35%
3.20.619180.626011.10%0.618700.0777%0.634282.44%
3.40.633650.683937.94%0.644051.64%0.655423.44%
3.60.646840.7618117.8%0.7250412.1%
3.80.65890
4.00.66996
1101010
Table 4. 30th-order numerical solutions V b n = 30 by IHAM ( ε = 8 / π 2 , ħ = 0.1 , ħ 1 = 0.3 ) compared with exact solutions V b .
Table 4. 30th-order numerical solutions V b n = 30 by IHAM ( ε = 8 / π 2 , ħ = 0.1 , ħ 1 = 0.3 ) compared with exact solutions V b .
α Exact Solution ε 1 = 0
(HAM)
ε 1 = 24 π 2 / π 3
(IHAM)
V b V b Δ V b Δ
000000
0.30.0989910.0990380.0481%0.0990380.0482%
0.60.192350.192660.160%0.192660.162%
0.90.276210.276970.275%0.276990.282%
1.20.349010.350250.355%0.350300.370%
1.50.410980.412620.400%0.412710.422%
1.80.463260.465220.423%0.465340.449%
2.10.507320.509530.437%0.509660.462%
2.40.544550.546990.449%0.547090.467%
2.70.576180.578830.459%0.578880.467%
3.00.603250.606090.470%0.606050.464%
3.30.626580.629590.479%0.629450.457%
3.60.646840.650010.490%0.649750.450%
3.90.664550.667940.510%0.667460.438%
4.20.680140.683380.475%0.682820.393%
4.50.693970.692740.177%0.695710.251%
4.80.706290.684063.15%0.707240.135%
5.10.717350.657588.33%0.723490.856%
5.40.727310.761204.66%0.739881.73%
5.70.736350.685836.86%
6.00.744570.4495939.6%
11010
Table 5. Numerical solutions V b P of homotopy-Páde approximation by IHAM. ( n = 30 , ε = 8 / π 2 , ħ = 0.1 , ħ 1 = 0.3 ) compared with exact solutions V b .
Table 5. Numerical solutions V b P of homotopy-Páde approximation by IHAM. ( n = 30 , ε = 8 / π 2 , ħ = 0.1 , ħ 1 = 0.3 ) compared with exact solutions V b .
α Exact Solution ε 1 = 0
(HAM)
P a d e 11 , 11
ε 1 = 24 π 2 / π 3
(IHAM)
P a d e 15 , 15
V b V b Δ V b Δ
000000
0.50.162140.162340.120%0.162340.121%
10.301720.302650.307%0.302670.315%
1.50.410980.412620.400%0.412710.422%
20.493460.495590.433%0.495720.459%
2.50.555660.558170.452%0.558260.468%
30.603250.606100.471%0.606050.463%
3.50.640390.643600.501%0.643250.446%
40.669960.673790.571%0.672720.412%
4.50.693970.698990.723%0.696450.358%
50.713790.720961.00%0.715870.292%
5.50.730420.741021.45%0.732140.235%
60.744570.760072.08%0.746160.213%
6.50.756760.778772.91%0.758670.252%
70.767370.797543.93%0.770240.374%
7.50.776700.816645.14%0.781320594%
80.784980.836246.53%0.792220.922%
8.50.792390.856428.08%0.803191.36%
90.799060.877259.79%0.814391.92%
9.50.805100.8987111.6%0.825932.59%
100.810610.9208013.6%0.837883.36%
11010
Table 6. 30th-order numerical solutions U b n = 30 by IHAM ( ε = 8 / π 2 , ħ = 0.1 , ħ 2 = 0.3 ) compared with exact solutions U b .
Table 6. 30th-order numerical solutions U b n = 30 by IHAM ( ε = 8 / π 2 , ħ = 0.1 , ħ 2 = 0.3 ) compared with exact solutions U b .
α Exact Solution ε 2 = 0
(HAM)
ε 2 = 8 / π 2
(IHAM)
U b U b Δ U b Δ
000000
0.30.00589940.00590630.116%0.00590520.0993%
0.60.0224880.0225780.402%0.0225630.335%
0.90.0470340.0473590.690%0.0473000.565%
1.20.0764000.0770650.870%0.0769330.696%
1.50.107940.108940.923%0.108720.723%
1.80.139810.141040.886%0.140760.681%
2.10.170850.172230.804%0.171900.610%
2.40.200460.201880.705%0.201530.531%
2.70.228350.229730.603%0.229390.456%
3.00.254420.255700.504%0.255400.387%
3.30.278700.279850.411%0.279610.324%
3.60.301280.302260.324%0.302100.270%
3.90.322280.323160.274%0.323100.257%
4.20.341800.343340.448%0.343020.356%
4.50.359990.363711.03%0.360420.121%
4.80.376950.375250.449%0.361864.00%
5.10.392790.3192018.7%0.3044722.5%
5.40.407620.02701193.4%0.1492463.4%
5.70.421530.1090374.1%
6.00.43459
11010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Li, X.; Xie, C.; Huo, S. Explicit Solution to Large Deformation of Cantilever Beam by Improved Homotopy Analysis Method II: Vertical and Horizontal Displacements. Appl. Sci. 2022, 12, 2513. https://doi.org/10.3390/app12052513

AMA Style

Li Y, Li X, Xie C, Huo S. Explicit Solution to Large Deformation of Cantilever Beam by Improved Homotopy Analysis Method II: Vertical and Horizontal Displacements. Applied Sciences. 2022; 12(5):2513. https://doi.org/10.3390/app12052513

Chicago/Turabian Style

Li, Yinshan, Xinye Li, Chen Xie, and Shuhao Huo. 2022. "Explicit Solution to Large Deformation of Cantilever Beam by Improved Homotopy Analysis Method II: Vertical and Horizontal Displacements" Applied Sciences 12, no. 5: 2513. https://doi.org/10.3390/app12052513

APA Style

Li, Y., Li, X., Xie, C., & Huo, S. (2022). Explicit Solution to Large Deformation of Cantilever Beam by Improved Homotopy Analysis Method II: Vertical and Horizontal Displacements. Applied Sciences, 12(5), 2513. https://doi.org/10.3390/app12052513

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop