Next Article in Journal
Modeling Inventory Cost Savings and Supply Chain Success Factors: A Hybrid Robust Compromise Multi-Criteria Approach
Next Article in Special Issue
Detection and Analysis of Critical Dynamic Properties of Oligodendrocyte Differentiation
Previous Article in Journal
Stochastic Approach to Investigate Protected Access to Information Resources in Combined E-Learning Environment
Previous Article in Special Issue
On the Asymptotic Behavior of Noncanonical Third-Order Emden–Fowler Delay Differential Equations with a Superlinear Neutral Term
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Approximations to the (Un)forced Pendulum–Cart System: Ansatz and KBM Methods

by
Weaam Alhejaili
1,
Alvaro H. Salas
2 and
Samir A. El-Tantawy
3,4,*
1
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
2
FIZMAKO Research Group, Department of Mathematics and Statistics, Universidad Nacional de Colombia, Manizales 500001, Colombia
3
Department of Physics, Faculty of Science, Port Said University, Port Said 42521, Egypt
4
Research Center for Physics (RCP), Department of Physics, Faculty of Science and Arts, Al-Mikhwah, Al-Baha University, Al-Baha 1988, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(16), 2908; https://doi.org/10.3390/math10162908
Submission received: 18 July 2022 / Revised: 8 August 2022 / Accepted: 10 August 2022 / Published: 12 August 2022
(This article belongs to the Special Issue Mathematical Modeling and Simulation of Oscillatory Phenomena)

Abstract

:
In the present investigation, some novel analytical approximations to both unforced and forced pendulum–cart system oscillators are obtained. In our investigation, two accurate and effective approaches, namely, the ansatz method with equilibrium point and the Krylov–Bogoliubov–Mitropolsky (KBM) method, are implemented for analyzing pendulum–cart problems.The obtained results are compared with the Runge–Kutta (RK4) numerical approximation. The obtained approximations using both ansatz and KBM methods show good coincidence with RK4 numerical approximation. In addition, the global maximum error is estimated as compared to RK4 numerical approximation.

1. Introduction

The study of the mechanism of nonlinear oscillating systems is one of the most important and recent studies at the present time, which concerns many researchers, especially in the engineering field and in the field of automated systems and robotics [1,2,3,4,5,6,7,8,9,10,11]. Nonlinear control theory (NLCT) has received a lot of attention due to its technical importance and impact on numerous sectors of applications. Robotics, for example, is a popular application of NLCT. The inverted pendulum is significant and effective for modeling and designing the robotics control system. A cart–inverted pendulum system is one of the most effective models used in robotic systems. The pendulum–cart system (PCS) is a nonlinear, underactuated system with unstable zero dynamics that must be regulated to maintain its unstable equilibrium [12,13,14,15,16]. The most popular approach for swinging up an inverted pendulum is energy control, which involves regulating the system’s energy rather than its velocity and position. In the literature and based on nonlinear analysis, intelligent control techniques have been proposed for controlling the nonlinear system [17,18,19,20,21].
There are many nonlinear evolution equations that govern pendulum oscillations, for instance, Duffing-type oscillators, Helmholtz-type oscillators, and many other oscillators related to the Duffing and Helmholtz oscillators [1,2,3,4,5,6]. On the other side, a benchmark example for underactuated mechanical systems is the PCS, which is highly nonlinear in nature and chosen as a typical experimental device in the control field. The PCS setup is suitable for testing various types of control theories and building robotics systems. The Segway is a two-wheeled balancing vehicle used for transport systems that can move forward by shifting the weight forward and can move backward by shifting the weight backward. Motivated by the potential applications of the pendulum–cart in robotic systems and many other applications, in this paper, we give an analytical solution to the PCS as shown in Figure 1 [12]. The differential equations of motion for the PCS are given by [12,13]:
R 1 = m ¯ x ¨ C + m l ( θ ¨ cos θ θ ˙ 2 sin θ ) F = 0 , R 2 = x ¨ C l cos θ + θ ¨ l 2 + g l sin θ = 0 ,
which is subjected to the following initial conditions (ICs):
θ ( 0 ) = θ 0 & θ ˙ ( 0 ) = 0 , x ( 0 ) = x 0 & x ˙ ( 0 ) = x ˙ 0 .
Here, θ θ ( t ) , x x t , and for simplicity we use m ¯ = ( M + m ) , where m and M denote the cart and pole mass, respectively, while l represents the distance from the center of rod axis rotation to the center of rod mass, g indicates the gravitational acceleration, and F expresses the excitation force. In this system, we assume that the excitation force F is constant. Numerous authors have used several approaches for solving and analyzing different types of equations of motion. For instance, the ansatz method, the Krylov–Bogoliubov–Mitropolsky (KBM) method, and He’s homotopy method were devoted to solving and analyzing the quadratically damped Duffing oscillator [22]. In addition, the nonlinear damped oscillations including the Morse oscillator and a nonlinear harmonic oscillator have been analyzed via the KBM method [23]. Moreover, the generalized KBM method was extended for investigating strong nonlinear oscillators with slowly varying parameters [24]. The (un)forced Duffing–Van der Pol oscillators have been investigated via the KBM method, He’s homotopy method, and He’s frequency–amplitude formulation [25]. The ansatz method is considered one of the most effective methods in finding many approximate analytic solutions to many evolution equations including both ordinary and partial differential equations [26,27,28]. For example, the ansatz method was used to find some approximate solutions to the families of Duffing-type equations, Helmholtz–Duffing-type equations, damped nonplanar Korteweg–de-Vries-type equations, Kawahara-type equations, and Muhammad equations, which have great credit for modeling many nonlinear phenomena that arise in many branches of science in general physics and plasma physics in particular [6,26,27,28,29,30]. Accordingly, the main objective of our investigation is to analyze the PCS (1) in order to find some analytical approximations using both the ansatz and KBM methods. For analyzing the PCS (1) via the ansatz method, we should determine the equilibrium point, while the KBM method does not need to find the equilibrium point.
This investigation is organized in the following fashion: The suggested method for analyzing the pendulum–cart system problem is discussed in detail in Section 1. The obtained approximations are discussed and verified by numerical examples in Section 2. The most important results are summarized in Section 3.

2. Mathematical Approaches for Analyzing System (1)

In this section, we try to find some approximations to system (1) using two different approaches. In the first approach, the ansatz method is employed for analyzing system (1), while in the second approach, the KBM method is applied to find an approximation to system (1).

2.1. First Approach: Ansatz Method

It is clear that system (1) is a coupled system in two independent variables x C , θ and for solving this system, we first eliminate x ¨ C from system (1), which leads to
R l m cos 2 θ l m ¯ θ ¨ cos θ F + l m θ ˙ 2 sin θ + g m ¯ sin θ = 0 .
Now, it is seen that Equation (3) is an ordinary differential equation (ode) in the variable θ . Using a Chebyshev approximation, the following approximation to the trigonometric functions is utilized for θ π / 3 :
sin θ θ θ 3 6 + θ 5 124 ,
cos θ 1 θ 2 2 + θ 4 25 .
According to the approximate formulas for sin θ and cos θ given in Equations (4) and (5), the approximate form of Equation (3) reads
θ 4 25 θ 2 2 + 1 2 l m l m ¯ θ ¨ = θ 4 25 θ 2 2 + 1 F + θ 5 124 θ 3 6 + θ θ ˙ 2 l m + g θ 5 124 θ 3 6 + θ m ¯ .
Now, for deriving an approximation to Equation (6), the following ansatz is presented:
θ = d + a cos ( w t ) , a = θ 0 d ,
where θ θ ( t ) and d indicates the equilibrium point. To determine the value of the equilibrium point d, the following MATHEMATICA command is utilized to solve system (1) at θ ˙ = 0 :
Solve R 1 = = 0 , R 2 = = 0 , x C [ t ] , θ t / / . θ t 0
which leads to
θ t = F cos θ + g m ¯ sin θ l m ¯ m cos θ 2 ,
x C [ t ] = F + g m cos θ sin θ m ¯ m cos θ 2 .
Note that x C [ t ] , θ t x ¨ C , θ ¨ . Now, solving Equation (8) or (9) using the following MATHEMATICA command,
Solve F cos θ t + g m ¯ sin θ t l m ¯ m cos θ t 2 = = 0 , θ t / / . C j _ : 0 1 / / FullSimplify
the equilibrium point for F 0 is obtained as
d θ t = tan 1 g m ¯ F 2 + g 2 m ¯ 2 , F F 2 + g 2 m ¯ 2 ,
which is equivalent to
d = tan 1 g m ¯ F .
Note that for F = 0 , we have d = 0 .
Inserting the ansatz (7) into Equation (6) yields
R = a 119040000 S cos w t + ,
with
S = 1575 a 10 l m w 2 + 75600 a 8 d 2 l m w 2 + 24094 a 8 l m w 2 + 378000 a 6 d 4 l m w 2 + 826980 a 6 d 2 l m w 2 1382000 a 6 l m w 2 + 403200 a 4 d 6 l m w 2 + 2760800 a 4 d 4 l m w 2 25088000 a 4 d 2 l m w 2 + 600000 a 4 g m + 600000 a 4 g M + 14632000 a 4 l m w 2 + 86400 a 2 d 8 l m w 2 + 1770944 a 2 d 6 l m w 2 34016000 a 2 d 4 l m w 2 + 7200000 a 2 d 2 g m + 7200000 a 2 d 2 g M + 117254400 a 2 d 2 l m w 2 + 14284800 a 2 d F x 14880000 a 2 g m 14880000 a 2 g M 59520000 a 2 l m w 2 + 190464 d 8 l m w 2 4761600 d 6 l m w 2 + 4800000 d 4 g m + 4800000 d 4 g M + 39283200 d 4 l m w 2 + 19046400 d 3 F x 59520000 d 2 g m 59520000 d 2 g M 119040000 d 2 l m w 2 119040000 d F x + 119040000 g m + 119040000 g M 119040000 l M w 2 .
Note that for vanishing S, we obtain an algebraic equation which can be utilized in determining the value of w as follows:
w 2 = S 1 S 2 ,
with
S 1 = 4800 125 a 4 g m + 125 a 4 g M + 1500 a 2 d 2 g m + 1500 a 2 d 2 g M + 2976 a 2 d F 3100 a 2 g m 3100 a 2 g M + 1000 d 4 g m + 1000 d 4 g M + 3968 d 3 F 12400 d 2 g m 12400 d 2 g M 24800 d F + 24800 g m + 24800 g M ,
and
S 2 = l 1575 a 10 m 75600 a 8 d 2 m 24094 a 8 m 378000 a 6 d 4 m 826980 a 6 d 2 m + 1382000 a 6 m 403200 a 4 d 6 m 2760800 a 4 d 4 m + 25088000 a 4 d 2 m 14632000 a 4 m 86400 a 2 d 8 m 1770944 a 2 d 6 m + 34016000 a 2 d 4 m 117254400 a 2 d 2 m + 59520000 a 2 m 190464 d 8 m + 4761600 d 6 m 39283200 d 4 m + 119040000 d 2 m + 119040000 M .
On the other side, for F = 0 , the value of w reads
w 2 = 120000 5 a 4 124 a 2 + 992 g m ¯ a 2 1575 a 8 + 24094 a 6 1382000 a 4 + 14632000 a 2 59520000 l m 119040000 l M .
Furthermore, the period of the oscillator is given by
T = 2 π w .
It remains to find the expression for x C . To this end, the first equation in system (1) is utilized as
x ¨ C = 1 m ¯ F + l m θ ˙ 2 sin θ l m θ ¨ cos θ .
Using the approximate formulas for sin θ and cos θ given in Equations (4) and (5) in Equation (18), we obtain
x ¨ C = l m m ¯ F l m + θ ˙ 2 θ 5 124 θ 3 6 + θ θ ¨ θ 4 25 θ 2 2 + θ .
The solution of Equation (19) under the ICs x ( 0 ) = x 0 and x ( 0 ) = x ˙ 0 reads
x C = x ˙ 0 t + l m 297600 m ¯ S 3 w 2 + 148800 F l m t 2 + j = 1 7 W j cos ( j w t ) ,
with
W 1 = a l m 375 a 6 + 6000 a 4 d 2 + 2480 a 4 + 6000 a 2 d 4 + 32736 a 2 d 2 74400 a 2 + 23808 d 4 297600 d 2 + 595200 595200 m ¯ , W 2 = a 2 d l m 375 a 4 + 23808 a 2 1200 d 4 + 48608 d 2 297600 1190400 ( m + M ) , W 3 = a 3 l m 75 a 4 + 3000 a 2 d 2 13640 a 2 + 6000 d 4 110112 d 2 + 223200 5356800 m ¯ , W 4 = a 4 d l m 125 a 2 + 500 d 2 4092 793600 m ¯ , W 5 = a 5 l m 225 a 2 + 3000 d 2 7688 14880000 m ¯ , W 6 = 5 a 6 d l m 142848 m ¯ , W 7 = a 7 l m 388864 m ¯ , S 3 = 375 a 6 d + 1500 a 4 d 3 372 a 4 d + 600 a 2 d 5 496 a 2 d 3 .
For F = 0 , solution (20) reduces to
x C ( t ) = 30675 a 7 l m + 370636 a 5 l m 8544375 a 3 l m + 51266250 a l m + 25633125 F x t 2 + 2 ( m + M ) t x ˙ 0 + x 0 51266250 m ¯ 75 a 6 + 496 a 4 14880 a 2 + 119040 a l m 119040 m ¯ cos ( w t ) + 15 a 4 2728 a 2 + 44640 a 3 l m 1071360 m ¯ cos ( 3 w t ) + 225 a 2 7688 a 5 l m 14880000 m ¯ cos ( 5 w t ) + a 7 l m 388864 m ¯ cos ( 7 w t ) .
Now, for examining the accuracy of the obtained approximations, the following data in the absence (existence) of the excitation force are considered: M = 2 kg, m = 1 kg, l = 1 m, g = 9.8 m/s 2 , θ 0 = 20 , θ ˙ 0 = 0 , x 0 = 1 m, and x ˙ 0 = 0 , while the force is given by
F = 0 , F = 0.01 N ,
Based on the mentioned data, the i.v.p. for the unforced case to be solved reads
θ ¨ cos θ θ ˙ 2 sin θ + 3 x ¨ = 0.01 , θ ¨ + 9.8 sin θ + x ¨ cos θ = 0 ,
and
θ ( 0 ) = 20 , θ ˙ ( 0 ) = 0 , x ( 0 ) = 1 , x ˙ ( 0 ) = 0 .
The approximate solution based on the relations (7) and (20) to this problem reads
θ = 0.349066 cos ( 3.74987 t ) ,
and
x = 1.114 0.11459 cos ( 3.74987 t ) + 0.000586335 cos ( 11.2496 t ) 8.9 × 10 7 cos ( 18.7494 t ) + 5.41 × 10 10 cos ( 26.2491 t ) .
Using the same mentioned data for the forced case ( F 0 ), we have
θ ¨ cos ( θ ) θ ˙ 2 sin ( θ ) + 3 x ¨ = 0.01 , θ ¨ + 9.8 sin ( θ ) + x ¨ cos ( θ ) = 0 ,
with
θ ( 0 ) = 20 , θ ˙ ( 0 ) = 0 , x ( 0 ) = 1 , x ˙ ( 0 ) = 0 .
The numerical values to the approximate solutions of this problem read
θ ( t ) = 0.349406 cos ( 3.74972 t ) 0.000340136 ,
and
x C ( t ) = 1.11411 + 0.00166669 t 2 0.114699 cos ( 3.74972 t ) 3.4265836921198773 × 10 6 cos ( 7.49943 t ) + 0.000588042 cos ( 11.2491 t ) + 8.68089 × 10 9 cos ( 14.9989 t ) 8.936889 × 10 7 cos ( 18.7486 t ) + 5.44996 × 10 10 cos ( 26.248 t ) .
The global maximum distance error L g in the whole domain is calculated for the unforced and forced cases as shown in Table 1.
Both solutions (26) and (27) are presented in Figure 2. In Figure 2a, the angle profile θ according to solution (26) is considered, while the profile of position x is introduced in Figure 2b. Moreover, the forced case is introduced in Figure 3, in which the angle profile θ according to solution (30) is considered, while the profile of position x given in (31) is displayed. One can see from Figure 2 and Figure 3 as well as the global maximum error given in Table 1 that there is complete congruence between the analytical and numerical approximations, which confirms the high accuracy and effectiveness of the obtained approximations.

2.2. Second Approach: KBM Method

From Equation (3), we obtain
θ ¨ θ ( t ) = F cos ( θ ) + g m sin ( θ ) + g M sin ( θ ) + l m sin ( θ ) cos ( θ ) θ ˙ 2 l m cos 2 ( θ ) m ¯ .
Note here that F can take any value for (in)dependent time, i.e., F may be a constant value or F = F t . This property is not fulfilled in the above approach, in which the excitation force must be constant because equilibrium is considered only for autonomous systems.
Using series for the right-hand side of Equation (32) yields
F cos ( θ ) + g m sin ( θ ) + g M sin ( θ ) + l m sin ( θ ) cos ( θ ) θ ˙ 2 l m cos 2 ( θ ) m ¯ = j = 0 5 Z j θ j ,
with
Z 0 = F l m m ¯ , Z 1 = g m ¯ l m m ¯ + m θ ˙ 2 m m ¯ , Z 2 = F m + m ¯ 2 l m m ¯ 2 , Z 3 = 2 l m 2 θ ˙ 2 + 5 g m m ¯ + 4 l m θ ˙ 2 m ¯ + g m ¯ 2 6 l m m ¯ 2 , Z 4 = F 5 m 2 + 18 m m ¯ + m ¯ 2 24 l m m ¯ 3 , Z 5 = 16 l m 3 θ ˙ 2 + 61 g m 2 m ¯ + 88 l m 2 θ ˙ 2 m ¯ + 58 g m m ¯ 2 + 16 l m θ ˙ 2 m ¯ 2 + g m ¯ 3 120 l m m ¯ 3 .
Then, Equation (32) can be written in the following compact form:
R = θ ¨ + g m ¯ l M θ + R θ , θ ˙ , F ,
with
R θ , θ ˙ , F = m θ ˙ 2 θ m m ¯ + Z 0 + Z 2 θ 2 + Z 3 θ 3 + Z 4 θ 4 + Z 5 θ 5 .
In order to apply the KBM method, we rewrite problem (34) in the following form of a p-problem (here, p < < 1 ) to enable us to solve this equation:
θ ¨ + g m ¯ l M θ + p R θ , θ ˙ , F = 0 , θ ( 0 ) = θ 0 and θ ( 0 ) = θ ˙ 0 .
For applying the KBM method on problem (35), we assume the solution in the following ansatz form:
θ θ p ( t ) = a cos ψ + n = 1 N 1 p n u n a , ψ + O p N ,
where each u n is a periodic function of ψ and both a a t and ψ ψ t are defined by
d a d t = n = 1 N 1 p n A n ( a ) + O ( p N ) ,
and
d ψ d t = ω 0 + n = 1 N 1 p n ψ n ( a ) + O ( p N ) ,
where ω 0 2 = g m ¯ / l M .
We define the KBM homotopy as follows:
H ( θ , p ) = θ ¨ + ω 0 2 θ + p R θ , θ ˙ , F .
The approximate analytical solution is obtained at p = 1 .
Now, to uniquely determine the values of the quantities A n and ψ n , it is required that no u n contains secularity terms. For sake of simplicity, we assume that N = 1 , and applying the KBM yields
θ p ( t ) = a cos ( ψ ) p 46080 g m ¯ M 3 Y 0 + Y 2 a 2 + Y 3 a 3 + Y 4 a 4 + Y 5 a 5 + Y 7 a 7 ,
with
Y 0 = 46080 M 3 F , Y 2 = 3840 M 2 F ( 2 m + M ) ( cos ( 2 ψ ) 3 ) , Y 3 = 240 g M 2 m ¯ ( 12 m + M ) cos ( 3 ψ ) , Y 4 = 16 M F 24 m 2 + 20 m M + M 2 20 cos ( 2 ψ ) + cos ( 4 ψ ) 45 , Y 5 = g M m ¯ 15 144 m 2 + 76 m M + M 2 cos ( 3 ψ ) + 240 m 2 + 140 m M + M 2 cos ( 5 ψ ) , Y 7 = g m m ¯ 15 m 2 + 15 m M + 2 M 2 6 cos ( 3 ψ ) + 6 cos ( 5 ψ ) + cos ( 7 ψ ) ,
while the values of a ˙ , ψ ˙ read
a ˙ = 0 , ψ ˙ = g m ¯ l M + p g m ¯ l M m a 6 15 m 2 + 15 m M + 2 M 2 + M a 4 96 m 2 + 44 m M + M 2 24 M 2 a 2 ( 4 m + M ) 384 M 3 .
We assume that a = A = constant, while the value of ψ reads
ψ = λ 384 M 3 g m ¯ l M A 6 m 15 m 2 + 15 m M + 2 M 2 + A 4 M 96 m 2 + 44 m M + M 2 24 A 2 M 2 ( 4 m + M ) + 384 M 3 t + B .
Here, the numbers λ and p are free/optimal parameters which are chosen to obtain the least residual error as possible, and of which the default values read λ = p = 1 . The value of integration constant B can be obtained from the ICs.
In order to find the value of x C ( t ) , we solve the second equation in system (1), i.e., R 2 = 0 , for x ¨ C :
x ¨ C l cos θ + θ ¨ l 2 + g l sin θ = 0 ,
We obtain
x ¨ C = θ ¨ l 2 + g l sin θ l cos θ .
Now, applying the Chebyshev approximation on the RHS of Equation (44) yields
x ¨ C = 864 g θ 3 49 π 3 18636 g θ 6125 π 34 39 l θ 2 θ ¨ 35 37 l θ ¨ .
By solving the ode (45) using RK4 or MATHEMATICA command “NDSolve” with the help of ICs x ( 0 ) = x 0 and x ( 0 ) = x ˙ 0 , we can obtain the expression for x C .
For numerical examples, the comparison between the numerical approximation using the RK4 method and the obtained approximations (40) and (45) is presented in Figure 4, Figure 5 and Figure 6 for different values of the excitation force F as follows:
F = 0 , F = 0.01 N , F = 0.1 cos 0.2 t N ,
while the other parameters have the same values given in Figure 2 and Figure 3.
The global maximum error L g is estimated for the three cases of the excitation force given in Equation (46), as illustrated in Table 2.

3. Conclusions

The equations of motion for unforced and forced pendulum–cart system oscillators have been solved and analyzed analytically using two accurate and effective approaches. In the beginning, the ansatz method was applied for deriving some effective and accurate approximations in terms of the trigonometric functions. Moreover, the KBM method was carried out to find some approximations to the problem under consideration. In the first approach, it was assumed that the excitation force F must be constant because equilibrium is considered only for autonomous systems. However, in the second approach, the excitation force F could take any value for dependent or independent time, i.e., F could be a constant value or F = F t . This feature is not fulfilled in the first approach, in which the excitation force must be constant, because equilibrium is considered only for autonomous systems. Furthermore, the analytical approximations were compared with the numerical approximations using the RK4 numerical approach. In addition, the global maximum distance error with respect to the RK4 numerical approximation was estimated. It was observed that the obtained approximations are characterized by high accuracy and convergence for long periods of time. Moreover, it was found that the accuracy of the obtained approximations using the first approach (ansatz method) is better than that of the second approach (KBM method).
The obtained results can help many authors in studying the mechanism of robotics motion and in investigating many nonlinear oscillations in plasma physics and several fields of science.

Author Contributions

Conceptualization, W.A. and A.H.S.; Data curation, W.A. and S.A.E.-T.; Formal analysis, W.A. and S.A.E.-T.; Investigation, W.A., A.H.S. and S.A.E.-T.; Methodology, A.H.S. and S.A.E.-T.; Project administration, S.A.E.-T.; Resources, S.A.E.-T.; Software, A.H.S. and S.A.E.-T.; Supervision, S.A.E.-T.; Writing – original draft, W.A. and A.H.S.; Writing – review & editing, S.A.E.-T. All authors contributed equally and approved the final manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article (more details can be requested from El-Tantawy).

Acknowledgments

The authors express their gratitude to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R229), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Albalawi, W.; Salas, A.H.; El-Tantawy, S.A.; Youssef, A.A.A. Approximate analytical and numerical solutions to the damped pendulum oscillator: Newton–Raphson and moving boundary methods. J. Taibah Univ. Sci. 2021, 15, 479–485. [Google Scholar] [CrossRef]
  2. Salas, A.H.; El-Tantawy, S.A. On the approximate solutions to a damped harmonic oscillator with higher-order nonlinearities and its application to plasma physics: Semi-analytical solution and moving boundary method. Eur. Phys. J. Plus 2020, 135, 833. [Google Scholar] [CrossRef]
  3. He, J.; El-dib, Y.; Mady, A.A. Homotopy Perturbation Method for the Fractal Toda Oscillator. Fractal Fract. 2021, 5, 93. [Google Scholar] [CrossRef]
  4. Aljahdaly, N.H.; El-Tantawy, S.A. On the multistage differential transformation method for analyzing damping Duffing oscillator and its applications to plasma physics. Mathematics 2021, 9, 432. [Google Scholar] [CrossRef]
  5. El-Tantawy, S.A.; Salas, A.H.; Alharthi, M.R. A new approach for modelling the damped Helmholtz oscillator: Applications to plasma physics and electronic circuits. Commun. Theor. Phys. 2021, 73, 035501. [Google Scholar] [CrossRef]
  6. Salas, A.H.; El-Tantawy, S.A.; Alharthi, M.R. Novel solutions to the (un)damped Helmholtz-Duffing oscillator and its application to plasma physics: Moving boundary method. Phys. Scr. 2021, 96, 104003. [Google Scholar] [CrossRef]
  7. Pedersen, H.B.; Madsen, M.L.; Andersen, J.E.V.; Nielsen, T.G. Investigation of argumental oscillations of a physical pendulum. Eur. J. Phys. 2021, 42, 025012. [Google Scholar] [CrossRef]
  8. Korytov, M.S.; Shcherbakov, V.S.; Titenko, V.V.; Belyakov, V.E. Spherical pendulum model with a moving suspension point in the problem of spatial load movement by a hoisting crane with oscillation limiting. J. Phys. Conf. Ser. 2020, 1441, 012102. [Google Scholar] [CrossRef]
  9. Johannessen, K. An approximate solution to the equation of motion for large-angle oscillations of the simple pendulum with initial velocity. Eur. J. Phys. 2010, 31, 511. [Google Scholar] [CrossRef]
  10. Butikov, E.I. Oscillations of a simple pendulum with extremely large amplitudes. Eur. J. Phys. 2012, 33, 1555. [Google Scholar] [CrossRef]
  11. Uy, R.F.; Yuan, C.; Chai, Z.; Khor, J. Wilberforce pendulum: Modelling linearly damped coupled oscillations of a spring-mass system. Eur. J. Phys. 2022, 43, 015011. [Google Scholar] [CrossRef]
  12. Owen, F. Pendulum on a Cart via Lagrangian Mechanics; Alpha Omega Engineering, Inc.: Sutton, MA, USA, 2014; Available online: http://www.aoengr.com/Dynamics/PendulumOnCart.pdf (accessed on 23 May 2014).
  13. Singla, A.; Singh, G. Real-Time Swing-up and Stabilization Control of a Cart-Pendulum System with Constrained Cart Movement. Int. J. Nonlinear Sci. Numer. Simul. 2017, 18, 525–539. [Google Scholar] [CrossRef]
  14. Neria, M.R.; Ramirez, S.H.; Moctezuma, R.G.; Juarez, A.L. Linear active disturbance rejection control of underactuatedsystems: The case of the Furuta pendulum. ISA Trans. 2014, 53, 920–928. [Google Scholar] [CrossRef] [PubMed]
  15. Kumar, V.; Jerome, J. Robust LQR controller design for sta-bilizing and trajectory tracking of inverted pendulum. Procedia Eng. 2013, 64, 169–178. [Google Scholar] [CrossRef]
  16. Ogata, K. Modern Control Engineering; Prentice-Hall: New York, NY, USA, 2000. [Google Scholar]
  17. Mori, S.; Nishihara, H.; Furuta, K. Control of unstablemechanical system- control of pendulum. Int. J. Control 1976, 23, 673–692. [Google Scholar] [CrossRef]
  18. Anderson, C.W. Learning to control an inverted pendulum usingneural network. IEEE Control Syst. Mag. 1989, 9, 31–37. [Google Scholar] [CrossRef]
  19. Boubaker, O. The inverted pendulum benchmark in nonlinearcontrol theory. Int. J. Adv. Rob. Syst. 2013, 10, 1–9. [Google Scholar]
  20. Song, Z.; Song, X.; Liu, C.; Zhao, Y. Research on real-timesimulation and control of linear 1-stage inverted pendulum. J. Comput. 2013, 8, 896. [Google Scholar] [CrossRef]
  21. Eide, R.; Egelid, M.; Stamso, A.; Karimi, H.R. LQG controldesign for balancing an inverted pendulum mobile robot. Intell. Control Autom. 2011, 2, 160–166. [Google Scholar] [CrossRef]
  22. Salas S, A.H. Analytical Approximant to a Quadratically Damped Duffing Oscillator. Sci. World J. 2022, 2022, 3131253. [Google Scholar] [CrossRef] [PubMed]
  23. Yamgoué, S.B.; Kofané, T.C. Application of the Krylov–Bogoliubov–Mitropolsky method to weakly damped strongly non-linear planar Hamiltonian systems. Int. J. Non-Linear Mech. 2007, 42, 1240–1247. [Google Scholar] [CrossRef]
  24. Cai, J. A Generalized KBM Method for Strongly Nonlinear Oscillators with Slowly Varying Parameters. Math. Comput. Appl. 2007, 12, 21–30. [Google Scholar] [CrossRef]
  25. Salas, A.H.; Albalawi, W.; El-Tantawy, S.A.; El-Sherif, L.S. Some Novel Approaches for Analyzing the Unforced and Forced Duffing–Van der Pol Oscillators. J. Math. 2022, 2022, 2174192. [Google Scholar] [CrossRef]
  26. Alyousef1, H.A.; Alharthi, M.R.; Salas, A.H.; El-Tantawy, S.A. Optimal analytical and numerical approximations to the (un)forced (un)damped parametric pendulum oscillator. Commun. Theor. Phys. 2022. Available online: https://iopscience.iop.org/article/10.1088/1572-9494/ac7bdc/pdf (accessed on 24 June 2022).
  27. Salas, A.H.; Albalawi, W.; Alharthi, M.R.; El-Tantawy, S.A. Some Novel Solutions to a Quadratically Damped Pendulum Oscillator: Analytical and Numerical Approximations. Complexity 2022, 2022, 7803798. [Google Scholar] [CrossRef]
  28. Alyousef, H.A.; Salas, A.H.; Alkhateeb, S.A.; El-Tantawy, S.A. Some Novel Analytical Approximations to the (Un)damped Duffing–Mathieu Oscillators. J. Math. 2022, 2022, 2715767. [Google Scholar] [CrossRef]
  29. El-Tantawy, S.A.; Salas, A.H.; Alharthi, M.R. On the analytical and numerical solutions of the damped nonplanar Shamel Korteweg–de Vries Burgers equation for modeling nonlinear structures in strongly coupled dusty plasmas: Multistage homotopy perturbation method. Phys. Fluids 2021, 33, 043106. [Google Scholar] [CrossRef]
  30. El-Tantawy, S.A.; Salas, A.H.; Alharthi, M.R. On the Analytical and Numerical Solutions of the Linear Damped NLSE for Modeling Dissipative Freak Waves and Breathers in Nonlinear and Dispersive Mediums: An Application to a Pair-Ion Plasma. Front. Phys. 2021, 9, 580224. [Google Scholar] [CrossRef]
Figure 1. A diagram of the pendulum–cart system.
Figure 1. A diagram of the pendulum–cart system.
Mathematics 10 02908 g001
Figure 2. The RK4 numerical approximation versus the analytical approximations ((a) Solution (26) and (b) Solution (27)) using the ansatz technique for F = 0 .
Figure 2. The RK4 numerical approximation versus the analytical approximations ((a) Solution (26) and (b) Solution (27)) using the ansatz technique for F = 0 .
Mathematics 10 02908 g002
Figure 3. The RK4 numerical approximation versus the analytical approximations ((a) Solution (30) and (b) Solution (31)) using the ansatz technique for F = 0.01 N.
Figure 3. The RK4 numerical approximation versus the analytical approximations ((a) Solution (30) and (b) Solution (31)) using the ansatz technique for F = 0.01 N.
Mathematics 10 02908 g003
Figure 4. The RK4 numerical approximation versus the analytical approximations ((a) Solution (40) and (b) Solution (45)) using the KBM method for F = 0 .
Figure 4. The RK4 numerical approximation versus the analytical approximations ((a) Solution (40) and (b) Solution (45)) using the KBM method for F = 0 .
Mathematics 10 02908 g004
Figure 5. The RK4 numerical approximation versus the analytical approximations ((a) Solution (40) and (b) Solution (45)) using the KBM method for F = 0.01 N.
Figure 5. The RK4 numerical approximation versus the analytical approximations ((a) Solution (40) and (b) Solution (45)) using the KBM method for F = 0.01 N.
Mathematics 10 02908 g005
Figure 6. The RK4 numerical approximation versus the analytical approximations ((a) Solution (40) and (b) Solution (45)) using the KBM method for F = 0.1 cos 0.2 t N.
Figure 6. The RK4 numerical approximation versus the analytical approximations ((a) Solution (40) and (b) Solution (45)) using the KBM method for F = 0.1 cos 0.2 t N.
Mathematics 10 02908 g006
Table 1. The global maximum distance error L g for the ansatz method.
Table 1. The global maximum distance error L g for the ansatz method.
F L g θ L g x C
0 0.00232 0.00076
0.01 0.002314 0.0007565
Table 2. The global maximum distance error L g for the KBM method.
Table 2. The global maximum distance error L g for the KBM method.
F L g θ L g x C
0 0.041186 0.0112186
0.01 N 0.0155981 0.0126209
0.1 cos 0.2 t N 0.0752997 0.458301
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alhejaili, W.; Salas, A.H.; El-Tantawy, S.A. Novel Approximations to the (Un)forced Pendulum–Cart System: Ansatz and KBM Methods. Mathematics 2022, 10, 2908. https://doi.org/10.3390/math10162908

AMA Style

Alhejaili W, Salas AH, El-Tantawy SA. Novel Approximations to the (Un)forced Pendulum–Cart System: Ansatz and KBM Methods. Mathematics. 2022; 10(16):2908. https://doi.org/10.3390/math10162908

Chicago/Turabian Style

Alhejaili, Weaam, Alvaro H. Salas, and Samir A. El-Tantawy. 2022. "Novel Approximations to the (Un)forced Pendulum–Cart System: Ansatz and KBM Methods" Mathematics 10, no. 16: 2908. https://doi.org/10.3390/math10162908

APA Style

Alhejaili, W., Salas, A. H., & El-Tantawy, S. A. (2022). Novel Approximations to the (Un)forced Pendulum–Cart System: Ansatz and KBM Methods. Mathematics, 10(16), 2908. https://doi.org/10.3390/math10162908

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