Abstract

In this paper, a mathematical model for large deformation of a cantilever beam subjected to tip-concentrated load is presented. The model is governed by nonlinear differential equations. Large deformation of a cantilever beam has number of applications is structural engineering. Since finding an exact solution to such nonlinear models is difficult task, this paper focuses on developing soft computing technique based on artificial neural networks (ANNs), generalized normal distribution optimization (GNDO) algorithm, and sequential quadratic programming (SQP). The strength of ANN modeling for governing the equation of cantilever beam is exploited by the global search ability of GNDO and further explored by the local search mechanism of SQP. Design scheme is evaluated for different cases depending on variations in dimensionless end-point load . Furthermore, to validate the effectiveness and convergence of algorithm proposed technique, the results of the differential transformation method (DTM) and exact solutions are compared. The statistical analysis of performance indicators in terms of mean, median, and standard deviations further establishes the worth of ANN-GNDO-SQP algorithm.

1. Introduction

Mechanical systems, which involve nonlinearity due to large deflection of compliant mechanism, continue to be an interesting problem. Since large deformation phenomena occur often in various geotechnical practices [1]. Therefore, this topic is of practical interest and has been widely studied by many researchers. In past decades, a number of numerical methods have been developed to solve large deformation problems. The large deflection of cantilever beams was studied by Wang [2]. They developed nonlinear differential equations for postbuckling loads on the basis of Eringen’s nonlocal constitutive relation. Shooting method was used to obtain postbuckling load and the buckled shape of the beam. Framework of arbitrary Lagrangian–Eulerian (ALE) is commonly used approach to study large deformation in geo-technical engineering [3]. Based on ALE three different approaches named as the efficient ALE approach (EALE) [4, 5], interpolation technique by small strain (RITSS) [6, 7] and the Abaqus built-in coupled Eulerian–Lagrangian (CEL) method [8]. Modified Chebyshev’s polynomial is used by Schmidt and Dadeppo [9] to study large deflection of beam. Large deformation of a spring-hinged beam was investigated by Nageswara Rao [10, 11] subjected to a tip rotational concentrated and distributed load.

Ludwick’s large deformation was studied by Lee [12] to made cantilever beams with a combined loading effect of a focused load at the tip and a uniformly dispersed load over the beam length. Phungpaigram and Chucheepsakul [13, 14] used elliptic integrals to calculate exact solutions for large deflection in elastic beams with variations in arc length and inclined force. Dado and Al-Sadder [15, 16] investigated the behaviour of large deformation of prismatic and nonprismatic cantilever beams under various type of loadings. Wang [17, 18] used homotopy perturbation method (HPM) to obtain analytical solution for large deformation of beam under point load with free tip. References [19, 20] studied the longitudinal vibration analysis for microbars based on strain gradient elasticity theory. Feasibility of Adomian decomposition method for such complex nonlinear problems was studied by Tolou and Herder [21]. Mutyalarao [22] studied large deflections of a uniform cantilever beam with concentrated load at tip and having normal inclination to the deflected axis of the beam. All these recently introduced techniques have their own grains and limitations in terms of accuracy, robustness, convergence, and applicability, but they are based on well-established deterministic procedures. Complex nonlinear differential equation of large deformation of beams subjected to a concentrated load is of great interest in scientific and engineering field. The quest of finding analytical solution for the problem motivates the author to develop a soft computing technique based on feed forward artificial neural networks (ANNs).

Stochastic solvers based on computational intelligence methods using artificial neural networks (ANNs) are considered to be fundamental in pattern recognition and machine learning. In general, neural networks are widely used to solve fractional differential equations, integro-differential equations (IDEs), partial differential equations (PDEs), and ordinary differential equations (ODEs). In ANN modeling, the optimization procedure is carried by using combination of global and local search algorithms. Some recent application of stochastic algorithm are multiphase flow through porous media for imbibition phenomena [23], eye model [24], wire coating dynamics [25, 26], optimal design and temperature distribution of heat fin [27, 28], beam-column design [29], and hybrid feature analysis for diabetic retinopathy classification using fundus images [30].

In the present study, a novel soft computing technique is applied to find analytical series solution for large deformation of cantilever beam subjected to point load by using ANN models optimized globally with generalized normal distribution optimization (GNDO) algorithm hybrid with sequential quadratic programming (SQP) for rapid local convergence. The prominent features of the present study are given as follows.(i)Mathematical model for large deformation of cantilever beam under end point load is formulated and analyzed to study the influence of variations in dimensionless end point load .(ii)A novel soft computing paradigm is developed to model series solutions based on artificial neural networks with the generalized normal distribution optimization algorithm and the sequential quadratic programming. Our approach is named as ANN-GNDO-SQP algorithm.(iii)To validate the efficiency of the proposed technique, four cases of large deformation cantilever beam are considered. The statistical results are compared with differential transformation method (DTM) and analytical solutions.(iv)Statistical analysis of absolute errors (AE), fitness evaluation (Fit), mean absolute deviation (MAD), Theil’s inequality coefficient (TIC), root mean square error (RMSE), Nash–Sutcliffe efficiency (NSE), and error in Nash–Sutcliffe efficiency (ENSE) are presented in terms of minimum, mean, median, and standard deviation.(v)Provision of continuous solutions, computational complexity, convergence, and easy execution of the proposed methodology show the robustness and correctness of the ANN-GNDO-SQP algorithm.

2. Problem Formulation

Consider a cantilever beam with a large deformation due to end point load as shown is Figure 1. By Euler–Bernoulli beam theory, the curvature of the beam can be given aswhere represents the slope or rotation of a beam, is the natural distance from the fixed end, is Young’s modulus, is bending moment, is moment of inertia, and is the bending stiffness of cantilever beam. From Figure 1, moment for the deflected beam under end point load is given aswhere concentrated load at end point is denoted by , is the horizontal deflection, and is the distance of deflected beam from fixed point. Hence, bending equation of uniform cross sectional beam for large deformation is written as

If , then concentrated force is a dead force such as gravity, and if , then the concentrated force is perpendicular to the deflected beam at the end [31] (Figure 1).

Differentiation of (1) with respect to is given asconsidering a dimensionless parameter . Differentiating (2) with respect to and substituting in (4) will result in governing second-order differential equation for large deformation of continuer beam with end point load aswith boundary conditionswhere , and it represents the dimensionless load at end point. Angle of rotation of the beam at free end is denoted by . Furthermore, dimensionless horizontal displacement of the free end is given by [17, 32, 33] asthen, the dimensionless horizontal displacement of the free tip is

For large deformation equation (5) along with boundary conditions, (6) and (7) are given as

The analytical solution for the problem is given asat the tip .

3. Proposed Methodology

The proposed soft computing paradigm for calculating approximate solutions for mathematical model of large deformation of cantilever beam consists of two parts. In the first part, an unsupervised ANNs model is constructed in terms of input, hidden, and output layers for governing differential equation. In the second part, neurons in ANN structure are trained or tuned by hybridizing generalized normal distribution optimization (GNDO) algorithm and sequential quadratic programming (SQP).

3.1. Construction of ANN Model

Feedforward artificial neural networks (ANNs) are used to model series solutions for governing equation of large deformation of cantilever beam with end point loading. Neural network model for (5)–(7) are formulated using continuous mapping approach for the solution and its respective derivatives in terms of input, hidden, and output layer are given as follows:where , , and are real-valued vectors and are bounded, is the activation function, and represents the number of neurons in ANN structure. In the hidden layer, (13)–(15) used log sigmoid as an activation function, and then, the updated form of solution and its derivatives is given as

3.2. Formulation of Fitness Function

Fitness function or objective function for mathematical model of large deformation in cantilever beam is developed as a sum of two mean square error (MSE):where and are MSE of (5) and (6) and (7), respectively, which are defined as

For the approximate solution of (5), fitness function equation (17) is optimized by training neurons in such a way that MSEs, and , should approach to zero, and consequently, the approximate solution by proposed method will converge to exact solution.

3.3. Training of Neurons

Methodologies adopted for training of unknown neurons in ANNs structure for optimization of fitness function equation (17) are presented, which is based on hybridization of unsupervised and supervised learning of GNDO and SQP, respectively. The working mechanism of the proposed algorithm is provided in Figure 2.

3.4. Generalized Normal Distribution Optimization

Generalized normal distribution optimization (GNDO) algorithm is a novel metaheuristic technique presented by Zhang et al. [34], inspired by generalized normal distribution theory. GNDO algorithm is widely used for parameter extraction of model, unlike other metaheuristic algorithms; GNDO is easy to implement, and it only requires the essential population size and termination criteria. GNDO has a simple structure, where the position of each individual is updated by using normal distribution curve. The working strategy of GNDO algorithm is subdivided into two phases, exploitation and exploration.

3.4.1. Exploitation

Exploitation is a process of finding best solution around the search space consisting of the current positions of all individuals. Initially, model for optimization by generalized distribution model is given aswhere , , , and , are trial vector, generalized mean position, generalized standard variance, and penalty factor, respectively. Moreover, , , and are formulated aswhere is the mean position, is the current best so far, and , , , and are the random numbers between 0 and 1. Furthermore, , , and are discussed in the exploration phase.

3.4.2. Exploration

Exploration refers to the searching of population space to get best solution. Exploration of GNDO is based on three randomly selected individuals as follows:where and are the trail vectors, is the adjustment parameter, and and are the random numbers between 0 and 1, which are subjected to the standard normal distribution. Trail vectors are computed as follows:where , and are integers. It is worth mentioning that GNDO algorithm is inspired by the relationship between normal distribution law and traditional teaching phenomena, search process of metaheuristics, and group teaching phenomena, respectively. GNDO has been applied to study the parameter extraction of photovoltaic models [34].

3.5. Sequential Quadratic Programming

The best performance (weights) obtained by GNDO algorithm is refined by the process of hybridization with efficient local search technique known as sequential quadratic programming using MATLAB toolbox setting. SQP is one of the powerful methods for numerical solution of constrained nonlinear optimization problems. It was developed in 1963 and further refined in 1970 [35]. SQP has been applied to a number of problems, which proves its power, accuracy, and efficiency. Nocedal and Wright [35] discuss SQP in detail and also give a mathematical formulation for various large-scale numerical optimization problems. Some recent applications of SQP are numerical solution for transient heat conduction problem [36], profile error evaluation of free-form surface [37], nonlinear model predictive control [38], OPF problem in DC grids [39], Bagley–Torvik systems arising in fluid mechanics [40], and optimal design of heating system in rapid thermal cycling blow mold [41].

3.6. Hybrid ANN-GNDO-SQP Algorithm

Necessary details of the procedural steps for proposed algorithm are given as follows.Step 1. Initialization of GNDO: Initial weights are created randomly from population space with number of entries equal to number of neurons in ANN structure. Mathematical formulation is given aswhere , , and are the real values of unknown neurons in ANN model. Parameter setting for GNDO algorithm is given in Table 1.Step 2. Fitness evaluation: Objective function equation (17) is evaluated to calculate the fitness value for the problem using the weights generated in the previous step.Step 3. Criteria for termination: Execution of GNDO is stopped when any of the following criteria is satisfied.Objective value ie Function tolerance ie Predefined number of iterations is achievedIf stopping criteria are fulfilled, then go to step V otherwise continue.Step 4. Storage: Store the global best weight corresponding to minimum fitness value and time taken for the execution.Step 5. Hybridization: Global best weights obtained by GNDO for minimization of (21) are considered as an initial guess for SQP to start the procedure.Step 6. Fitness evaluation: SQP starts the supervised learning, update the weights, and evaluate the fitness function until the following terminations conditions are satisfied.Objective value, i.e., .Predefined number of iterations is achieved.Step 7. Storage: Store the best weight, minimum fitness value, and time taken for the execution by SQP and the total time by GNDO-SQP in seconds.

Repeat the procedure from steps I–VII for a sufficient large number of independent runs to generate a large dataset for reliable statistical analysis.

4. Performance Indices

In this section, the performance of design scheme for solving mathematical model of large deformation of cantilever beam is examined by incorporating performance indicators in terms of mean absolute deviation (MAD), Theil’s inequality coefficient (TIC), root mean square error (RMSE), and Nash–Sutcliffe efficiency (NSE). Mathematical formulations of these indicators are given as follows [23].where is the analytical solution and represents the approximate solution by proposed algorithm. denotes the grid points.

5. Numerical Simulation and Discussion

In this section, different cases of (5) are considered to study the effect of variations in dimensionless end point load on large deflection of cantilever beam. The following cases are considered. Case I: , Case II: , Case III: , and Case IV: . The formulation of fitness functions for each case is given as follows:

In this paper, the mathematical model of large deflected cantilever beam under end point load is investigated by developing soft computing technique. ANNs-based fitness function is constructed to model approximate solutions, which are further optimized by using hybridization of GNDO and SQP algorithms. To briefly study the efficiency and behaviour of proposed technique, it is executed for 100 times. Results obtained by ANN-GDDO-SQP algorithm are compared with exact solution and differential transformation method [33]. Approximate solutions along with absolute errors obtained by the proposed algorithm for different cases of deflected cantilever beam are shown in Figure 3. Figures 4(a) and 4(b) illustrate the influence of variations in dimensionless end point load on horizontal distance and rotational angle of the beam at the tip or free end, respectively. Convergence of fitness values, MAD, TIC, RMSE, and ENSE for each case of large deflected cantilever beam during 100 independent runs is shown in Figures 5 and 6. Figure 7 represents the boxplots analysis for each case of deflected cantilever beam. The graphical illustration of mean absolute values, global values for fitness function, and performance indicators for each case study are presented in Figure 8.

It can be seen that the approximate solutions overlap the exact solution with minimum errors that show the accuracy of proposed technique. Table 2 dictates the comparison of approximate solutions and analytical solutions for different cases depending on variations in dimensionless point load at free end. Tables 3 and 4 show that absolute errors in best solutions for case I-IV lie around to , to , to and to with standard deviations to , to , to and to , respectively. Table 5 dictates that mean or global values of fitness function for each case study are , , , and . Also, minimum values of MAD, TIC, RMSE, and ENSE lie around , , , and with standard deviations around to , to , to and to , respectively. Analyses based on the computational complexity of the design scheme for obtaining solution to equations (25)–(28) are dictated in Table 6. Weights in ANN structure for best solution of each case obtained by proposed algorithm are presented in Table 7 and graphically shown in Figure 9.

6. Conclusion

In this paper, we have analyzed the mathematical model of large deflected beam under variational point load and the free end. Furthermore, we present a new soft computing technique to calculate approximate solution for different cases depending on variations in end point load. We conclude our finding as follows.(i)A new soft computing evolutionary algorithm is developed in which the strength of artificial neural networks are utilized to model an approximate series solution and combined with hybridization of generalized normal distribution optimization algorithm and sequential quadratic programming. The proposed algorithm is named as ANN-GNDO-SQP algorithm.(ii)ANNs-based fitness functions are constructed and optimize with the ANN-GNDO-SQP algorithm to obtain an overlapping solutions with minimum absolute errors as shown in Figure 3.(iii)It can be seen that increasing end-point load increases the angle of deflection, horizontal distance, and rotational angle of the cantilever beam.(iv)Convergence graphs and boxplots of 100 independent executions for fitness evaluation, MAD, TIC, RMSE, and ENSE show the stability of proposed algorithm.(v)Extensive graphical and statistical analysis along with complexity analysis of the proposed algorithm for solving large deflected cantilever beam shows the correctness and robustness of ANN-GNDO-SQP algorithm.

Approximate solutions for cases I–IV are as follows:

Abbreviations

ANNs:Artificial neural networks
GNDO:Generalized normal distribution optimization
MAD:Mean absolute deviation
TIC:Theil’s inequality coefficient
NSE:Nash–Sutcliffe efficiency
ENSE:Error in Nash–Sutcliffe efficiency
SQP:Sequential quadratic programming
RMSE:Root mean square error
DTM:Differential transform method
:Distance of deflected beam
:Concentrated load
:Unknown real valued neurons in ANNs
:End-point load
:Curvature
:Rotation of beam
:Distance from fix point
:Inertia
:Bending moment
:Bending stiffness
:Horizontal deflection of beam
:Dimensionless parameter
:Trail vector
:Mean position
:Standard variance
:Trail vector
:Mean position
:Standard variance.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the teaching research and reform project of Shanghai DianJi University (No. 6).