Abstract

This paper analyzed the model of swinging oscillation of a solid circular sector arising in hydrodynamical machines, electrical engineering, heat transfer applications, and civil engineering. Nonlinear differential equations govern the mathematical model for frequency oscillation of the system. Furthermore, a computational strength of Cascade neural networks (CNNs) is utilized with backpropagated Levenberg–Marquardt (BLM) algorithm to study the oscillations in angular displacement , velocity , and acceleration . A data set for the supervised learning of the CNN-BLM algorithm for different angles and radius are generated by Runge–Kutta (RK-4) method. The BLM algorithm further interprets the dataset with log-sigmoid as an activation function for the solutions’ validation, testing, and training. The results obtained by the design scheme are compared with Akbari–Ganji’s (AG) method. The rapid convergence and quality of the solutions are validated through performance indicators such as mean absolute deviations (MAD), root means square error, and error in Nash–Sutcliffe efficiency (ENSE). The statistics demonstrate the design scheme’s applicability and efficiency to highly singular nonlinear problems.

1. Introduction

Nonlinear oscillation and its behavior is an important topic in applied physics, mathematics, and mechanical engineering that has piqued the interest of many scientists from the dawn of human recognition of the equations of motion [1]. In vibrations theory, an extremely unique and important place is held by oscillatory systems that contain fundamental nonlinearities. On the other hand, it is of the utmost significance to comprehend their behavior, not only from an academic standpoint but also in view of the countless potential applications [24]. The systems of complex nonlinear and chaotic responses are characterized by a unique equilibrium position, a strong nonlinearity described by a monotone increasing mooring restoring force, and a fluid structure [57]. The dynamical systems, including spherical pendulum, harmonic oscillator model of aromaticity, spring systems, symmetric, and biased hardening duffing oscillators, produce the nonlinear and modulated responses of multidegree exhibits the periodic vibration, saturation, thermal resonance, and waves [8, 9]. Fox and Goulbourne [10, 11] carried out experiments to study primary and superharmonic resonances on the nonlinear oscillation of dielectric elastomers subjected to applied static pressure and dynamic voltage. A. Alibakhshi [12] investigated the static and dynamic response of the free and forced vibrations of a micro/nanobeam made of a hyperelastic material incorporating strain-stiffening, size effect, and moderate rotation.

The oscillation of a rigid rod over a circular surface is considered a classical oscillator whose mathematical model was first investigated by Gaylord [13] in 1979. Nonlinearity is an inherent property of every differential equation that describes a physical or biological phenomenon. Solving linear differential equations is quite straightforward and has been successfully implemented. On the other hand, the methods for solving nonlinear differential equations (NDEs) are not as readily available. Generally, NDEs do not have a precise solution. In general, finding exact and semianalytical solutions for such problems is challenging because of the nonlinearity in the elastic and damping components of the governing equations (14). NDEs have been the subject of all-embracing studies in various branches of nonlinear science and engineering. Poincare and Lyaponuv [15] were credited for being the first to solve the differential equations governing nonlinear processes. Their approaches were based on the existence of a values (large or small) parameter in the nonlinear equations, which allows the solution to be expanded as a power series whose components are generated by the integer powers of that small parameter. Using the series in the nonlinear model and equating the coefficients of the small parameters with like powers results in a system of linear ordinary differential equations whose solution caused the final result of the governing equation in closed form. This basic concept allowed the research community to solve numerous nonlinear scientific problems with a small or large parameter. The procedure was named the perturbation method. The method’s only flaw is that it requires the existence of a small parameter in the equations, even though we know that the proposed parameter does not exist in many types of differential equations. A series of analytical and semiexact techniques that do not require a small parameter was proposed in the middle of the twentieth century to overcome the problem. Some well-known techniques used for the numerical and analytical solution of nonlinear problems are Adomian decomposition method (ADM) [16, 17], energy balance method (EBM) [18], variational iteration method (VIM) [1921], extended direct algebraic method (EDAM) [22], Adams–Bashforth method [23], differential transform method [24], Hamiltonian approach [25, 26], Backlund method [27] rational harmonic balance method [28], min-max approach (MMA) [29, 30], power series technique [31], one-step hybrid block method [32], amplitude-frequency formulation (AFF) [33, 34], and modified homotopy perturbation method (MHPM) [3537]. Ebrahimi Khah [38] studied the motion of a rigid rod rocking back over a circular surface using He’s energy balance method. Galerkin method and modal analysis technique is utilized by G. Sheng [39, 40] to study the dynamic stability and nonlinear vibrations of the stiffened functionally graded (FG) cylindrical shell in a thermal environment. Reference [41] investigates the vibration analysis of a rigid rod using the Hamiltonian approach. Abul-Ez [42] employed the hybridization of the iteration perturbation method and variational iteration method to study the duffing oscillation of nonlinear oscillators. The above-cited literature shows that all of these techniques have been effectively applied to study the solutions and behavior of the nonlinear oscillatory models, but besides their advantages, it is observed with keen interest that such methods are gradient-based and require prior information about the problem. Prior information includes the choice of a small parameter, initial guess, continuity, and differentiability of the function. To overcome these drawbacks, gradient-free stochastic computing techniques based on the approximation ability of artificial neural networks are designed.

The researcher’s community has been widely employing the concepts of artificial neural network-based numerical computing techniques for exploring and exploiting linear/nonlinear mathematical systems. Some recent applications of the black-box stochastic methodologies include the approximate solutions for a mathematical model of chaotic base secure communication systems [43], nonlinear problems arising in heat transfer [44], thermal radiations of nanofluid [45, 46], nonlinear restoring moment and damping effects of ships [47], and wire coating dynamics [48]. These reported articles motivate authors to utilize the computational strength of artificial neural networks (ANNs) for the numerical treatment and analysis of nonlinear oscillation over a circular sector. The innovative insights of the presented study are summarized as follows:(i)A mathematical model for the oscillation of a homogenous solid circular sector object is analyzed to study the influence of variations in angles and radius on angular frequency, velocity, and acceleration.(ii)A novel application of artificial intelligence-based cascade neural networks via backpropagated Levenberg–Marquardt algorithm is presented effectively for the numerical solution of the nonlinear oscillator.(iii)Approximate solutions for different scenarios of the nonlinear problem obtained by the design CNN-BLM algorithm are compared with analytical solutions by Ranga–Kutta method and Akbari–Ganji’s method.(iv)The proposed algorithm is executed multiple times to test the results’ convergence, accuracy, and effectiveness. The performance of CNN-BLM is further validated by mean absolute deviations (MAD), root means square errors, and errors in Nash–Sutcliffe efficiency (ENSE).

2. Swinging Oscillation of a Solid Circular Sector Object

In this section, an oscillation of homogenous circular sector of angle and radius over a solid surface is considered as shown in Figure 1. The circular object is allowed to roll on a fixed solid surface producing an oscillatory motion (back and forth) with no sliding effect, it is evident that will change, and the system will satisfy the following nonlinear differential equation [49, 50]:

Introducing the dimensionless geometrical parameterwhere is the height of the mass center. Using (2) in (1) will give

Subjected to the initial conditionswhere is vibrational amplitude, denotes the angular displacement, and is the dimensionless time variable.

3. Design Methodology

3.1. Cascade Neural Networks

An artificial neural network (ANN) is a collection of interconnected, basic components known as neurons with multiple inputs and a single output, and each neuron represents a mapping. The output of a neuron is a function of the sum of its inputs which is generated with the help of the activation function. A neural network architecture that only contains the input and output layer is called perceptron [51]. The sum of weighted signals in perceptron is transmitted directly to the output layer from the input layer. In neural network modeling, multilayer perceptron (MLP), also called feed-forward neural network (FNN) comes with an additional layer between the input and output layer known as the hidden layer. In MLP, the signals that enter the hidden layer are processed by the transfer function to obtain the desired results. In FNN, the connection between the input and output layer is indirect while in perceptron, they possess a direct relation. The direct network between input and output layers formulated by the combination of MLP and perceptron is named as cascade neural networks (CNNs) [52]. The basic structure of CNN is based on two ideas. First, is to build an architecture of cascade by adding new neurons with their connection to all inputs and previously hidden neurons without changing configuration at each layer. The second idea is to minimize the residual error by training only the newly created neurons by fitting their weights. The network is updated with new neurons and the performance increases.

Initially, the cascade neural networks start the process with inputs and a single targeted date without hidden neurons. The output neuron is connected to each input using weights such as with a standard sigmoid function. The Log-sigmoid is one of the best nonlinear normalized functions, which is differentiable and continuous; therefore, it gives a smooth gradient preventing jumps in the output values. The output of the network is presented as

The new neurons formed are added to the network, and each neuron is connected to hidden neurons. A suitable algorithm trains the neurons for output. The basic advantage of CNN is that no structure for the network is predefined, it automatically built its architecture from the training data. Secondly, the learning process of CNN is fast compared to other networks because each neuron is trained independently.

3.2. Learning Procedure and Performance Measures

This section discusses the working and training procedure of neurons in CNN architecture. The learning strategy of the neurons in CNN is based on two phases. Initially, a network is presented with inputs of the reference data set of 1001 points generated by the Ranka–Kutta method. Secondly, the inputs and the weights in CNN are trained by an appropriate algorithm Levenberg–Marquardt with a log-sigmoid activation function for the output. The architecture of CNN along with the workflow of the design algorithm is shown in Figure 2.To examine the accuracy and effectiveness of the results of the CNN-BLM algorithm for oscillations of the nonlinear circular sector, performance indices are defined in terms of mean square error (MSE) of fitness function of the model, regression , absolute errors (AE), mean absolute deviations (MAD), root mean square error (RMSE), and error in Nash–Sutcliffe efficiency (ENSE). Formulation of these indices is given as [53]here, shows the number of grid points, , and are the reference, approximate, and mean of solution at input. The desired value of MSE, AE, MAD, RMSE, and ENSE for perfect fitting is equal to zero while the value of and NSE should be one.

4. Numerical Experimentation and Discussion

In this section, the artificial intelligence-based stochastic numerical technique is implemented on a nonlinear mathematical model of oscillations in a circular section given by equations (3)−(4) to study the angular displacement, velocity, and acceleration. Different cases of the problem based on variations in and are given in Table 1.

Approximate solutions for angular displacement of solid circular sector obtained by the designs scheme for different values of radius and angle are shown in Figure 3. It can be seen that by reducing the angle with a constant semicircular radius, the frequency of oscillation increases. Also, reducing the semicircular radius causes an increase in the frequency of the oscillation. It is worthy to note that the frequency of oscillation determined by the design scheme is higher than the frequency obtained by Akbari–Ganji’s (AG) method for different cases of scenario 1 as shown in Table 2. Table 3 compares the proposed algorithm’s approximate solutions with numerical solutions by the RK-4 method. Results of CNN-BLM for different scenarios of (3) overlap the analytical solution with minimum absolute errors, as shown in Table 4.

To study the smoothness of solutions, absolute errors for 1001 points are plotted in Figure 4. For scenario 1, the mean values of absolute errors lie around to while the maximum values of absolute errors lie around to at and 100. Similarly, the mean absolute errors for scenario 2 and 3 lies around to and to , respectively. Furthermore, to study the angular displacement and velocity, graphs of against are shown though Figure 5. Apart from the graphs of the solutions, Figure 5 depicts that the oscillation trends are somehow harmonic and stable. It can be seen that velocity reaches its peak when the angle or semicircular decreases. Surface plots are shown through Figure 6 for angular velocity and acceleration in the homogenous circular sector. The oscillations are all stable because neither successive increases in the amplitudes are noticed in the displacement time diagram, nor open curves have been detected in the phase plane.To study the performance of the design scheme, the CNN-BLM algorithm is executed for 20 independent runs. Results of performance function in terms of mean square errors are dictated in Tables 5 and 6 and plotted through Figure 7. Mean values of MSE lie around to . Figure 8 shows the convergence of mean absolute deviations of the solutions at each step size. It can be seen that the mean value of MAD for each case lies around to . Tables 7 and 8 dictates the statistics of root mean square error (RMSE) and error in Nash–Sutcliffe efficiency (ENSE) in terms of minimum, mean and standard deviations. Additionally, regression plots are shown in Figure 9. It can be seen that the values of regression for each case are equal to unity which shows the efficiency, correctness, and accuracy of the results obtained by the CNN-BLM algorithm for oscillations in the homogenous circular sector.

5. Conclusions

In this paper, an artificial intelligent strength of cascade neural networks with backpropagated Levenberg–Marquardt algorithm is utilized to study the swinging oscillation of a solid circular sector over a solid surface. The governing relation of the system was presented in the form of the differential equation. The proposed algorithm is implemented to study the effect of variations in semicircular radius and semicircular angle on angular displacement, velocity, and acceleration of the system. The results dictate that by reducing the angle with a constant semicircular radius, the frequency of oscillation increases, in these cases, higher velocities are achievable and the phase-plane ellipse height is greater. In addition, reducing the semicircular radius also causes an increase in the frequency of the oscillation. The approximate solution for angular displacement obtained by the CNN-BLM algorithm is compared with analytical solutions and approximate solutions by Akbari–Ganji’s (AG) method. Absolute errors demonstrate the effectiveness of the results as the solutions by the CNN-BLM algorithm overlap the numerical solution up to 8 and 9 decimal places. Statistics of MAD, RMSE, MSE, and ENSE lie between to , to , to , and to which reflects the accuracy of solutions and validates the worth, efficiency, and stability of the proposed algorithm.

In the future, the authors aim to implement cascade neural networks to study the solutions of nonlinear fractional differential equations modeling physical problems such as biochemical reactors, skin cancer, tumor growth, physiological temperature regulations, drying processes, and chaotic behaviors in warless communications.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors equally contributed to this manuscript and approved the final version.