Next Article in Journal
Study of Temperature Field Distribution in Topographic Bias Tunnel Based on Monitoring Data
Next Article in Special Issue
A Multiplicity Theorem for Superlinear Double Phase Problems
Previous Article in Journal
Positive Solutions of a Fractional Boundary Value Problem with Sequential Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Control Theory Application for Swing Up and Stabilisation of Rotating Inverted Pendulum

Department of Mechatronics, Faculty of Mechanical Engineering, University of Pristina, 10000 Prishtina, Kosovo
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(8), 1491; https://doi.org/10.3390/sym13081491
Submission received: 15 July 2021 / Revised: 10 August 2021 / Accepted: 10 August 2021 / Published: 13 August 2021
(This article belongs to the Special Issue Nonlinear Analysis and Its Applications in Symmetry)

Abstract

:
This paper introduces a new scheme for sliding mode control using symmetry principles for a rotating inverted pendulum, with the possibility of extension of this control scheme to other dynamic systems. This was proven for swing up and stabilisation control problems via the new sliding mode control scheme using both simulations and experiments of rotary inverted pendulum (RIP) underactuated systems. According to the Lyapunov theory, a section of the pendulum was compensated with a scale error in the upright position, as the desired trajectory was followed by the pendulum arm section. As the RIP’s dynamic equations were nonlinearly complex and coupled, the complex internal dynamics made the task of controller design difficult. The system control for the pathway of the reference model of the rotational actuator with the application of the sliding mode technique for moving back and forth up the inverted pendulum’s structure, till the arm to reach the linear range round the vertical upright position, was created and tested in an existent device. The stabilisation scheme was switched on in the sliding mode as soon as the arm reached the linear range. A comparison of the stabilisation performance for the same rotating inverted pendulum as discussed by other authors revealed that the proposed controller was more flexible and reliable in terms of the swing up and stabilisation time.

1. Introduction

There are three main components of a rotary inverted pendulum (RIP), namely, link 1 and link 2 of the pendulum, and a motor (Figure 1). The working mechanism of the RIP framework is that link 1 is driven by the motor with the appropriate control order to such an extent that the descending link 2 can be adjusted in the vertical position [1]. It is well known in the domain of a control hypothesis that RIP is one of the most complex control problems [2]. Wiener, Johnson, et al. based their work on humans and other animals to study their movements because some types of machines behave similarly [3]. Humans and animals control their mechanisms based on environmental variables when they are in a reliable terrain while moving or performing certain tasks [3].
The authors in [4] regarding the design of general feedback theory of human behaviour, started with theoretical model before outlining the application of the physical model.
There is another idea that its components map onto the adverse feedback system within the existing models of behaviour for artificial control systems [5]. The movement control of a two-legged robot is based on the inverse pendulum theory [6]. Its good example is an RIP system for nonlinear underactuated mechanical systems to practice the control theory and with special emphasis on robotics for applications [7]. Saleem et al. and Jadlovska et al. provided the mathematical equations and numerous calculations, to find the appropriate way for the nominal trajectory to control the RIP [8]. Wen et al. based their work [9] on the Lyapunov theory control for nonlinear dynamics associated with the RIP in an optimal control effort. In our case, for RIP, we solved this problem by using Kour and an LQR controller; we compared our results with the other authors’ [10,11,12,13,14] because these authors examined the same problem with LQR [5,15]. Some authors [16,17,18,19,20] have tested and implemented the theory of fractional calculus with a PID controller [21]. Irrespective of their prominence, traditional PID controllers come up short on the global asymptotic steadiness verification [21,22]. The contribution of this article is the amalgamation of an upgraded and adjustment control scheme for the RIP system. On a very basic level, the proposed scheme serves simultaneously as the adjuster and the controller. These controllers utilise the state input of the link 1 angle (β), the link 2 angle (α), and their time subordinates to produce the ideal force (or voltage) control orders to upright balance link 2, while continuously turning link 1 as close as is conceivable to its unique position, separately [23,24,25,26]. In this study, we realised the occasion-based procedures to control a rotary pendulum, for both swinging up and balancing out its arm [27]. Swinging up a rotary pendulum is an old and testing issue in the field of nonlinear control research [28,29,30]. The rest of this paper is organised as follows: The RIP dynamics and its state space models are discussed in Section 4, Section 5 and Section 6. The stabilisation and swing up controllers principle is shown in Section 6. Section 7 and Section 8 show the theoretical and experimental results, and finally, in the end, a conclusion is presented.

2. Rotary Inverted Pendulum

The model of rotational link 2 is shown in Figure 2. The rotating joint is appended to the servo engine shaft. The pivoting arm [rotating link] has a length, Larm a moment of inertia, Jarm and an angle α; α increments decidedly when turning the counterclockwise. The pendulum interface is associated with the finish of the pivoting joint. It has an all-out connection Lpend and a focal point of mass mpend. The moment of inertia for its center of mass is Jpend. The angle of inverted link 2, β, is zero when the inverted arm is in an upright position; by dividing the symmetry space along the z axis, and it increases in the positive direction when it moves in the clockwise direction [18,21]. There are six main components of the RIP, as shown in Figure 1: link 2, link 1, a motor with a gearbox ratio of 20:1, a base, and two encoders [1]. The RIP system is depicted in Figure 1, along with the geometric parameters and system variables.
The RIP has two planes of symmetry through the z0 axis with a normal direction of x0 and y0; it can rotate in the plane perpendicular to la = y1.
This is attributed to its symmetrical configuration. Therefore, the distance xCoM from the centre of mass (CoM) can be calculated as follows:
x CoM = i = 1 n m i · x i i = 1 n m i = m a 1 · x a 1 + m a 2 · x a 2 + m a 3 · x a 3 m a 1 + m a 2 + m a 3
x C o M distance from some reference point to the (CoM), and xi: distance to the CoM for each individual object.
The movement of the mechanism can be described in a graph-analytical or analytical manner.
From Figure 2 is derived the kinematics equation of link 2.
r B / A = 1 2 · L pend · sin β · ı + cos β · ȷ
v B / A = 1 2 · L pend · cos β · β ˙ · ı sin β · β ˙ · ȷ )
v A / O = L arm · α ˙ · ı  
v B / O = v B / A + v A / O = L arm · α ˙ 1 2 · L pend · cos β · ı 1 2 · L pend · β ˙ · sin β · ȷ
The Kalman filter is moderately less numerically confounded and simpler to use than the other filters, for example, the applied filter. Therefore, the limit of the Kalman filter to situate nonlinear coordinated frameworks precisely is restricted [19,23,24,25]. The all-inclusive Kalman filter is expressed in Algorithm 1. In the accompanying Figure 3 and Figure 4, we applied the Kalman filter on the grounds that the clamours from the encoder were enormous.
Algorithm 1 Description of the filtered angular velocity/acceleration.
Input: angular position, angular velocity, angular acceleration
Output: u β
 for u β   a n g u l a r   d a t a   d o
    β = P p e n d - p r e d C T C · P · C T + R 1
x β ^ = a n g u l a r k 1 + β · y k 1 C · a n g u l a r k 1
P = I β · C · L p e n d
a n g u l a r k = d i s c r e t e   o u r   m o d e l x β 1 , ^ u β 1
F β = d i s c r e t e   o u r   m o d e l x x β 1 , ^ u β 1
P p e n d - p r e d = F β · L · F β T + Q
 end

3. DC Motor Modelling

The electric equal circuit of the armature and the free-body chart of the rotor are shown in Figure 5.
The input of the system is the feeding voltage (v) applied to the motor’s armature. The system output is the angular velocity of the shaft, marked as α ˙ t [20].
Therefore, we measured two parameters: armature resistance R and torque constant Kme. From the measured samples for the R, we took 30 measurements, and their average was R = 9.59 Ω.
v i t = i t · L d i d t + R · i t + e b t
T = B d α d t + J · d 2 α d t 2
e b = K m e · α ˙ t
T = K m e · i t
Upon applying the Laplace transform to Equations (6)–(9), we obtained the following expressions in the frequency domain:
Where: α ˙ t = ω t = Ω s ,   α ¨ t = ω ˙ t = s · Ω s , v i t = V i s and i t = I s .
Ω s V i s = K me J · L · s 2 + J · R + L · B · s + R · B + K me 2
The differential conditions from above can likewise be communicated in the state-space structure by selecting the engine position, engine speed, and armature current as the state variables [7]. Let the variables be d α ˙ t d t = α ¨ t = B J · α ˙ t K me J · i t   and   d i t d t = K me L · α ˙ t R L · i t + 1 L · v i t ; it is not difficult to get the accompanying state-space structure to depict in a linear form for the input and the output.
d d t · α ˙ t i t = B J K me J K me L R L · α ˙ t i t + 0 1 L · v i t ;   y t = 1 0 · α ˙ t i t
Figure 6 presents a scheme of the DC motor with feedback control, and Figure 7 shows the DC motor simulation with the real parameters by using MATLAB/Simscape. The block diagram of Figure 6 has been included within the Simscape/Simulink in the DC motor block of Figure 7, in particular.
It has been proven for the RIP’s equation of motion that the symmetry leads to the possibility of obtaining the solutions by using a real-time variable.
The steady-state velocities with different voltages obtained from the Simulink modes are shown in Figure 8.
The steady-state velocities with different voltages obtained from the real DC motor are shown in Figure 9.

4. Rotary Inverted Pendulum Dynamics

Mathematical modelling should be possible by utilizing the Euler–Lagrange strategy. In order to write the differential equations of the motion of the pendulum [31], taking into consideration the symmetry between two links, we started from the general form of the Lagrange equations of the second type.
L = 1 2 · m p e n d ( L a r m 2 ) 2 · α ˙ 2 + 1 2 · J a r m · α ˙ 2 + 1 / 2 · [ L a r m 2 · α ˙ 2 + ( L p e n d 2 ) 2 · β ˙ 2 + ( L p e n d 2 ) 2 · s i n 2 β · α ˙ 2 + L a r m · L p e n d · cos α · α ˙ · β ˙ ] + 1 2 · J p e n d · β ˙ 2 ( m a r m + m p e n d ) · g · h 1 2 · m p e n d · L p e n d · g · cos β
d d t L α ˙ L α = τ , d d t L β ˙ L β = 0
From Equations (11) and (12) we derived the followings:
m a r m ( L a r m 2 ) 2 · α ¨ + J a r m · α ¨ + J a r m · L a r m 2 α ¨ + m p e n d ( L a r m 2 ) 2 · s i n 2 β · α ¨ + m p e n d L a r m 2 + sin β · cos β · α ˙ · β ˙ + 1 2 · m p e n d · L p e n d · L a r m · cos β · β ¨ 1 2 · m p e n d · L p e n d · L a r m · sin β · β ˙ 2 = τ m p e n d ( L p e n d 2 ) 2 · β ¨ + J a r m · α ¨ + 1 2 · m p e n d · L p e n d · L a r m · cos β · β ¨ 1 2 · m p e n d · L p e n d 2 · cos β · sin β · α ˙ 2 1 2 · m p e n d · L p e n d · g · sin β = 0 } ,
q = α ,   β => q ˙ = α ˙ ,     β ˙ —controllable variables—equal to the degrees of freedom. As we had two controllable variables, the two degrees of freedom, respectively, in the matrix form will be as follows:
D q q ¨ + H q , q ˙ + G q = τ
where: D q —kinetic energy matrix, H q , q ˙ —centrifugal and Coriolis forces, G q —gravitational forces, and τ —torque.
The required partial derivatives for substitution in Equations (13)–(15) are as follows:
D q = J a r m + m p e n d · L a r m 2 + m p e n d · L p e n d 2 · s i n β m p e n d · L a r m · L p e n d · c o s β m p e n d · L a r m · L p e n d · c o s β J p e n d + m p e n d · L p e n d 2
H q , q ˙ L a r m 2 + m p e n d 2 · L a r m 2 + m p e n d · L p e n d 2 · β ˙ · s i n 2 β + m p e n d 2 · L p e n d 2 · α ˙ · s i n 2 β Γ m p e n d 2 · L p e n d 2 · s i n 2 β L p e n d 2
where: Γ = m p e n d · L p e n d · α ˙ · s i n β
G q = 0 m p e n d L p e n d 2 · g · s i n β ,
τ = K m e R V K m e 2 R · α ˙

5. Swing-Up and Stabilisation of RIP

The descending position to the vertical upright position and to hold it there. We used two control strategies: one for swinging up, and the other for balancing the pendulum parts with a small error to reach the dead position as soon as possible. The control using sliding models [32] was developed for both the methods of swinging up and stabilising 𝓏 = 𝓀 T 𝓏 . A subsequent order system for position control with a swaying joint movement that went about as the disturbance was mulled over for the pendulum of the swing up control. The following controller was structured with the goal that it followed a reference input in our case with a reference trajectory. Because of the symmetrical configuration of the RIP, this provided energy to the swing joint and placed it at the focal point, where it was utilised in another pendulum balancing control to balance out the unsteady points. For the trajectory, the quadrilateral wave of amplitude ±π/6 was used. The state-space equation for the stance control system [33] used was 𝓏 ˙ = A · 𝓏 + B · u + D , where the variables A and B are the system matrix and D is a concentrated disturbance that contains a disturbance due to an oscillating connection of the pendulum. Thus, 𝓏 e , ˙ z represented the error state vector, the desired state vector was z d , and the actual state vector was 𝓏 e = 𝓏 𝓏 d . The equation 𝓏 ˙ = Π · 𝓏 β · s i g n 𝓏 used differentiating sliding models 𝓏 ˙ = 𝓀 T 𝓏 ˙ , and then, we substitute the actual state vector 𝓏 ˙ e = 𝓏 ˙ z ˙ d . Thus, we obtained the sliding model 𝓏 ˙ = 𝓀 T z ˙ z ˙ d . If z d ± π 6 , 0 was the desired vector, the derivative vector was z ˙ d 0 , 0 . The conditions given above were to integrate the control u = 1 𝓀 T B · 𝓀 T A · z + Π · z β · s i g n z ; it was swing-up control, and the exchanging gain Π was more noteworthy than the greatest furthest reaches of the 𝓀 T D guarantee sliding [32]. The RIP system was a fourth-order system, and the state space was obtained with Equations (6) and (15), where z R 4 x 1 ,   A R 4 x 4 , and B R 4 x 1 . On the basis of the above equation, the incorporated control was u = 1 𝓀 T B · 𝓀 T A · z + Π · z β · s i g n z 𝓀 T z ˙ d and the swing up control is presented.

6. Simulation Results

The simulation is shown in Figure 10, using K o u r and the LQR controller. For our controller, the parameters used were K o u r = 2 ,   4 ,   6 ,   10 . The indicial parameter for K L Q R was K L Q R = 381.5 , 26.7 ,   3.08 ,   12.25 . The parameter of LQR was as follows: u = 5200 , i f   0 < β 90 °, α 0 and u = 5200 , i f 90 < β 0 °, α 0 , and Q = d i a g 90 ,   10 ,   40 ,   50 .
The reactions of the link 2 angle obtained by, respectively, utilising our strategy and the LQR method are presented in Figure 11, from which, one can infer that our swing-up technique was quicker than that utilising the LQR control scheme. The explanation was that, in the LQR control technique, link 2 needed to follow a to-and-fro movement until it arrived at a little neighbourhood of the vertical position [1]. Such a strategy would cost a considerable amount of time. Moreover, in our method, link 2 could be swung up with less movement because of the high-efficiency utilisation of the inertia [1]. Note that the multi-switching of the LQR controller might cause harm to the motor. In any case, the switching recurrence of our controller was considerably less than that of the LQR controller. A detailed performance comparison of our method and the LQR scheme is presented in Table 1, along with the swing up time, steady-state error, controller switched times, and controller saturation time [1] of our method; all of them were less than those of the LQR scheme. This detailed performance comparison revealed that our method was more efficient in terms of the swing up time and stabilisation for the same RIP.

7. Experimental Results

This section shows the experimental results of the combined sensitivity of our controller in comparison with LQR [28]. The platform of the RIP system is shown in Figure 12. The main controller chip was STM32F10 and Arduino Uno, the engine was MG513IPH, and the driver IC was TB 6612FNG [1]. The parameters for our controller were as follows: K o u r = 2 ,   4 ,   6 ,   10 . The linear feedback parameter for K L Q R was as follows: K L Q R = 381.5 , 26.7 ,   3.08 ,   12.25 . The parameter of LQR was set as follows: u = 5200 , i f   0 < β 90 °, α 0 , and u = 5200 , i f 90 < β 0 °, α 0 , and Q = d i a g 90 ,   10 ,   40 ,   50 . In the following text, we present the pseudo-code Algorithms 2–4 for our swing up method and the LQR for a low-cost RIP platform.
The following algorithms of the different calculations for the RIP case are given based on their writing according to the authors [34].
Algorithm 2 Parameters and description of the calculation of angle/velocity.
Input: n ,   m ,   R a r m a t u r e ,   K m e ,   V o u t ,   X v e l r a d ,   X a n g , X l a s t a n g
 for  P o s o l d = 0 : k   do
if k   = 0 or X v e l r a d > −1 and   X a n g X v e l r a d then
P o s o l d =   X l a s t a n g s t a r t A n g //   calculate the pendulum angle
end if
if    k > 0 then
      X v e l r a d =   ( X a n g X l a s t a n g ) //   calculate the pendulum velocity
else
//   use last velocity updated X v e l r a d
   end if
end for
end
Output: u β
Algorithm 3 Description of the calculation parameteres for swing up.
Input: v e l , n , m ,   l a s t v e l ,   v e l R a d ; a n g ,   l a s t a n g ,   s t a r t A n g ,   a c c
for l = 0:k do
if k = 0 or   x v e l > 1   a n d   l   x v e l then
x v e l = v e l l a s t v e l / /   use last updated.
ang =   l a s t a n g n / /   function   constraining   the   angle   between   π   and   π .
end if
if k > 0 then
      K a l m a n e s t i m . = e s t i m . E r r o r + K a l m a n g a i n M e a s c u r r . l a s t E s t i m .
ang = ( l a s t a n g π ) m / /   the calculation parameters for swing up
end if
end for
end
Output: n e w p o s . / /   through position we control the angle of the pendulum
Algorithm 4 Calculation of potential energy and torque.
Input :   E k , E p ,   m p e n d ,   L p e n d ,   x s p d ,   m a r m ,   L a r m ,   n ,   m ,   T m
for l = 0:k do
if k = 0 or a n g r a d . > 1 and l a n g r a d . then
  X a n g = n o l d p o s . o l d c u r r . / /   shift   of   the   angle
end if
if k > 0 then
x v e l R a d = m   X a n g X l a s t A n g . ; / /   new   shif   of   the   pendulum   angle
E p = m p e n d . g L p e n d . ( 1 cos ( a n g r a d . ) / /   calculate   potential   energy  
else
/ /   use   last   updated   x v e l R a d
end if
if   a n l . R e a d . 0 then
a n g = n a n l . R e a d . + m / /   pendulum   angle   for   swing   up β π
end if
if   x s p d > 0 then
x s p d = a n g K m e + m ; / /   calculate   of   the   speed  
end if
T m = ( x s p d V o u t m a x ) / v e l m a x / /   calculation   of   the   torque   for   the   swing   up
end for
end
Output β π and T m
The experimental results are shown in Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18. Figure 12 shows the responses of the pendulum angle, Figure 13 shows the responses of the pendulum angle velocity β ˙ , Figure 14 shows the responses of the pendulum acceleration swing unfiltered β ¨ , and Figure 15 shows the response of the arm velocity α ˙ ; these data are presented in the sample steps. The response of the controller is given in Figure 16. Figure 17 shows that link 2 swung up in approximately 1.45 s for K o u r = 10 , and then was stabilised by the control. It is clearer from Figure 18 that the nonlinear part of the stabilisation controller worked from less than 1000 sample steps to 5000 sample steps under our controller for K o u r = 10 , roughly. Note that all the signals of the RIP were stable. Furthermore, the swing up and equilibration time, steady-state error, and controller response time are presented in Table 2 below. Link 2′s angle   β position and link 1′s arm position α were in symmetry.
Figure 13 shows the reduced time for the swing up and stabilisation. The best value was observed in the case of K o u r = 10 ;   we observed that for the swing up, it took less than 1500 sample steps to get the pendulum angle in the vertical position. The worst value was obtained in the case of K o u r = 2 ;   here, we observed that the swing up and stabilisation took more than 3500 sample steps to get the pendulum angle in the vertical position.
Figure 14 shows the reduced time for the swing up and stabilisation. The best value was obtained in the case of K o u r = 10 ;   here, we observed that for the swing up, it took less than 1000 sample steps to get the pendulum angle in the vertical position. The worst values were obtained in the case of K o u r = 4   and K o u r = 2 ,   as the swing up and stabilisation took more than 1800 sample steps to get the pendulum velocity in the vertical position.
Hence, we reduced the time required for the swing up and stabilisation. The best value was obtained in the case of K o u r = 10 ; we observed that for the swing up, it took less than 800 sample steps to get the pendulum angle in the vertical position. The worst value was obtained in the case of K o u r = 2 ;   here, we observed that the swing up and stabilisation took more than 1000 sample steps to get the pendulum angle in the vertical position. Hence, Figure 15 shows the reduced time taken for the swing up and stabilisation. The best value was obtained in the case of K o u r = 10 ; we observed that for the swing up, it took less than 870 sample steps and the potential energy in the y direction was −655.84 to give the pendulum potential energy in vertical position up, technique reaction, and the adjustment procedure reaction of the RIPS by utilising the structured swing-up control.
The worst values were obtained in the case of K o u r = 2 ; we observed that the swing up and stabilisation took more than 1177 sample steps and the potential energy in the y direction was −424.69. The experiment exhibited the whole procedure of the dynamic reaction including the swing.
Figure 19 shows the motions of the RIP, in which 1, 2 illustrate the amplitude-increased simple pendulum motion, 3–5 demonstrate the swing up condition, 6–7 show the inertia swing up procedure, and 11–15 depict the stabilisation process. Figure 19 shows the caught movements of the RIPS, wherein 1–2 illustrate the amplitude-increased simple pendulum motion, 3–5 exhibits the swing up condition, 6, 7 show the inertia swing up technique, and 11–15 describe the stabilisation procedure.

8. Conclusions

In this paper, the control schemes for the RIP system that could be used to swing up the pendulum and balance it in the upright position were presented. A benchmark control, namely control of an RIP system, using sliding control was investigated. The position control of the arm coupled to the shaft of the actuator with sliding mode control was analysed, because of the symmetrical configuration, and to achieve the desired swing up from the RIP system’s position down has been used, to bring the pendulum into the linear range where the linear model is valid for control created. A balancing control was additionally utilised in the linear range. The sliding mode controller offered excellent swing up control and adjustment of the vertical upstanding position. The performance of controlling the RIP framework by utilising the sliding mode control for the swing up and the sliding model control for balancing out indicated the preferred performance to that of the controller with a PD controller for swing up and the LQR controller for balancing. The swing up and balancing time of under 1.5 s in the proposed strategy contrasted with the 3–4 s in the other technique was reduced.
The experimental results closely matched the simulation results. Compared with the LQR swing up controller and the single linear equilibration controller, the proposed controller was more flexible and trustworthy.

Author Contributions

Methodology, X.B. and A.P.; software, X.B. and, R.B.; validation, R.L.; formal analysis, A.S. and M.B.; investigation, X.B. and A.P.; resources, R.B.; data curation, X.B. and R.B.; writing—original draft preparation, A.P. and X.B.; writing—review and editing, X.B., A.P. and R.L.; visualisation, R.B. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest with respect to this paper.

Abbreviations

RIPRotary Inverted Pendulum
RIPSRotary Inverted Pendulum System
LQRLinear Quadratic Regulator
COMCentre of Mass
DCDirect Current
PIDProportional Integral Derivative

References

  1. Yang, X.; Zheng, X. Swing-Up and Stabilization Control Design for an Underactuated Rotary Inverted Pendulum System: Theory and Experiments. IEEE Trans. Ind. Electron. 2018, 65, 7229–7238. [Google Scholar] [CrossRef]
  2. Davison, E.J. (Ed.) Benchmark Problems for Control System Design: Report of the IFAC Theory Committee; International Federation of Automatic Control: Laxenburg, Austria, 1990. [Google Scholar]
  3. Wiener, N. Cybernetics or Control and Communication in the Animal and the Machine; Technology Press: Cambridge, MA, USA, 1948. [Google Scholar]
  4. Powers, W.T.; Clark, R.K.; Mc Farland, R.L. A General Feedback Theory of Human Behavior: Part I. Percept. Mot. Skills 1960, 11, 71–88. [Google Scholar] [CrossRef]
  5. Chawla, I.; Singla, A. Real-Time Control of a Rotary Inverted Pendulum using Robust LQR-based ANFIS Controller. Int. J. Nonlinear Sci. Numer. Simul. 2018, 19, 379–389. [Google Scholar] [CrossRef]
  6. Johnson, T.; Zhou, S.; Cheah, W.; Mansell, W.; Young, R.; Watson, S. Implementation of a Perceptual Controller for an Inverted Pendulum Robot. J. Intell. Robot. Syst. 2020, 99, 683–692. [Google Scholar] [CrossRef] [Green Version]
  7. Ismail, J.; Liu, S. Efficient Planning of Optimal Trajectory for a Furuta Double Pendulum Using Discrete Mechanics and Optimal Control. IFAC-PapersOnLine 2017, 50, 10456–10461. [Google Scholar] [CrossRef]
  8. Saleem, O.; Mahmood-Ul-Hasan, K. Robust stabilisation of rotary inverted pendulum using intelligently optimised nonlinear self-adaptive dual fractional-order PD controllers. Int. J. Syst. Sci. 2019, 50, 1399–1414. [Google Scholar] [CrossRef]
  9. Wen, J.; Shi, Y.; Lu, X. Stabilizing a Rotary Inverted Pendulum Based on Logarithmic Lyapunov Function. J. Control Sci. Eng. 2017, 2017, 1–11. [Google Scholar] [CrossRef]
  10. Wang, J.-J. Simulation studies of inverted pendulum based on PID controllers. Simul. Model. Pract. Theory 2011, 19, 440–449. [Google Scholar] [CrossRef]
  11. Prasad, L.B.; Tyagi, B.; Gupta, H.O. Optimal Control of Nonlinear Inverted Pendulum System Using PID Controller and LQR: Performance Analysis Without and With Disturbance Input. Int. J. Autom. Comput. 2014, 11, 661–670. [Google Scholar] [CrossRef] [Green Version]
  12. Bettou, K.; Charef, A. Control quality enhancement using fractional PIλDμ controller. Int. J. Syst. Sci. 2009, 40, 875–888. [Google Scholar] [CrossRef]
  13. Dwivedi, P.; Pandey, S.; Junghare, A. Novel fractional order PDμ controller for open-loop unstable inverted pendulum system. In Proceedings of the 17th International Conference on control, Automation and Systems, Jeju, Korea, 18–21 October 2017; 2017; pp. 1616–1621. [Google Scholar]
  14. Dwivedi, P.; Pandey, S.; Junghare, A. Performance analysis and experimental validation of 2-DOF fractional order controller for under actuated rotary inverted pendulum. Arab. J. Sci. Eng. 2017, 42, 5121–5145. [Google Scholar] [CrossRef]
  15. Dwivedi, P.; Pandey, S.; Junghare, A.S. Stabilization of unstable equilibrium point of rotary inverted pendulum using fractional controller. J. Frankl. Inst. 2017, 354, 7732–7766. [Google Scholar] [CrossRef]
  16. Dwivedi, P.; Pandey, S.; Junghare, A. Robust and novel two degree of freedom fractional controller based on two-loop topology for inverted pendulum. ISA Trans. 2018, 75, 189–206. [Google Scholar] [CrossRef] [PubMed]
  17. Shang, W.W.; Cong, S.; Li, Z.; Jiang, S.L. Augmented Nonlinear PD Controller for a Redundantly Actuated Parallel Manipulator. Adv. Robot. 2009, 23, 1725–1742. [Google Scholar] [CrossRef]
  18. Ozana, S.; Docekal, T.; Kawala-Sterniuk, A.; Mozaryn, J.; Schlegel, M.; Raj, A. Trajectory Planning for Mechanical Systems Based on Time-Reversal Symmetry. Symmetry 2020, 12, 792. [Google Scholar] [CrossRef]
  19. Kalman, R.E. A New Approach to Linear Filtering and Prediction Problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  20. Control Tutorials. University of Michigan. Available online: https://ctms.engin.umich.edu/CTMS/index.php?example=MotorSpeed&section=SystemModeling (accessed on 6 March 2021).
  21. Rotary Inverted Pendulum. Available online: https://www.lehigh.edu/~inconsy/lab/css/ME389/guidelines/ME389_MEM04_PendulumGantry_Guideline.pdf (accessed on 15 May 2021).
  22. Quanser. Available online: https://www.quanser.com/products/rotary-inverted-pendulum (accessed on 10 May 2021).
  23. Welch, G.; Bishop, G. An Introduction to the Kalman Filter. 1995. Available online: https://www.cs.unc.edu/~welch/media/pdf/kalman_intro.pdf (accessed on 10 May 2021).
  24. Vanicek, P.; Omerbasic, M. Does a navigation algorithm have to use a Kalman filter? CAASJ 1999, 45, 292–296. [Google Scholar]
  25. Grewal, M.S.; Andrews, A.P. Kalman Filtering: Theory and Practice Using MATLAB, 3rd ed.; Wiley: New York, NY, USA, 2001. [Google Scholar]
  26. Jadlovska, S.; Sarnovsky, J. A Complex Overview of Modeling and Control of the Rotary Single Inverted Pendulum System. Adv. Electr. Electron. Eng. 2013, 11, 73–85. [Google Scholar] [CrossRef]
  27. Durand, S.; Castellanos, J.F.; Marchand, N.; Sanchez, W.F. Event-based control of the inverted pendulum: Swing up and stabilization. J. Control Eng. Appl. Inform. 2013, 15, 96–104. [Google Scholar]
  28. Da Sanjeewa, S.; Parnichkun, M. Control of rotary double inverted pendulum system using mixed sensitivity H∞ controller. Int. J. Adv. Robot. Syst. 2019, 16, 1729881419833273. [Google Scholar] [CrossRef]
  29. Yoshida, K. Swing-up control of an inverted pendulum by energy-based methods. In Proceedings of the 1999 American Control Conference (Cat. No. 99CH36251), San Diego, CA, USA, 2–4 June 1999; Volume 6, pp. 4045–4047. [Google Scholar]
  30. Zhang, X.; Ma, J.; Lin, L.; Wang, L. Study on Swing-up Control of Rotary Inverted Pendulum Based on Energy Feedback. In Proceedings of the 2018 5th International Conference on Information Science and Control Engineering (ICISCE), Zhengzhou, China, 20–22 July 2008; pp. 994–998. [Google Scholar]
  31. Guo, X.; Zhang, G.; Tian, R. Periodic Solution of a Non-Smooth Double Pendulum with Unilateral Rigid Constrain. Symmetry 2019, 11, 886. [Google Scholar] [CrossRef] [Green Version]
  32. Edwards, C.; Spurgeon, S. Sliding Mode Control; CRC Press: London, UK, 1998. [Google Scholar] [CrossRef]
  33. Hong, Q.; Shi, Y.; Chen, Z. Adaptive Sliding Mode Control Based on Disturbance Observer for Placement Pressure Control System. Symmetry 2020, 12, 1057. [Google Scholar] [CrossRef]
  34. Tomescu, M.; Jäntschi, L.; Rotaru, D. Figures of Graph Partitioning by Counting, Sequence and Layer Matrices. Mathematics 2021, 9, 1419. [Google Scholar] [CrossRef]
Figure 1. (left) RIP system and (right) centre of mass.
Figure 1. (left) RIP system and (right) centre of mass.
Symmetry 13 01491 g001
Figure 2. Point B-CoM: centre of mass of link 2.
Figure 2. Point B-CoM: centre of mass of link 2.
Symmetry 13 01491 g002
Figure 3. Pendulum angular velocity before and after Kalman filtering.
Figure 3. Pendulum angular velocity before and after Kalman filtering.
Symmetry 13 01491 g003
Figure 4. Pendulum angular acceleration before and after Kalman filtering.
Figure 4. Pendulum angular acceleration before and after Kalman filtering.
Symmetry 13 01491 g004
Figure 5. DC motor in our RIP.
Figure 5. DC motor in our RIP.
Symmetry 13 01491 g005
Figure 6. DC motor: mathematical model.
Figure 6. DC motor: mathematical model.
Symmetry 13 01491 g006
Figure 7. DC motor simulation in the MATLAB/Simscape with real parameters.
Figure 7. DC motor simulation in the MATLAB/Simscape with real parameters.
Symmetry 13 01491 g007
Figure 8. Steady-state velocities calculated from Simulink model with 5.5 V.
Figure 8. Steady-state velocities calculated from Simulink model with 5.5 V.
Symmetry 13 01491 g008
Figure 9. Steady-state velocities calculated using the average value. This is the result obtained for K values of the six experiments for 2000 samples with 5.5 V and 11.04 V.
Figure 9. Steady-state velocities calculated using the average value. This is the result obtained for K values of the six experiments for 2000 samples with 5.5 V and 11.04 V.
Symmetry 13 01491 g009
Figure 10. Simulation of RIP.
Figure 10. Simulation of RIP.
Symmetry 13 01491 g010
Figure 11. RIP pendulum angle under K o u r and the LQR controller.
Figure 11. RIP pendulum angle under K o u r and the LQR controller.
Symmetry 13 01491 g011
Figure 12. RIP experimental platform.
Figure 12. RIP experimental platform.
Symmetry 13 01491 g012
Figure 13. Response of pendulum angle under our controller for K o u r with different parameters.
Figure 13. Response of pendulum angle under our controller for K o u r with different parameters.
Symmetry 13 01491 g013
Figure 14. Response of pendulum velocity under our controller for K o u r with different parameters.
Figure 14. Response of pendulum velocity under our controller for K o u r with different parameters.
Symmetry 13 01491 g014
Figure 15. Response of pendulum acceleration swing unfiltered under our controller for K o u r with different parameters.
Figure 15. Response of pendulum acceleration swing unfiltered under our controller for K o u r with different parameters.
Symmetry 13 01491 g015
Figure 16. Response of arm velocity in the case of our controller for K o u r with different parameters.
Figure 16. Response of arm velocity in the case of our controller for K o u r with different parameters.
Symmetry 13 01491 g016
Figure 17. Response of motor torque in the case of our controller for K o u r   with different parameters.
Figure 17. Response of motor torque in the case of our controller for K o u r   with different parameters.
Symmetry 13 01491 g017
Figure 18. Response of pendulum potential energy in the case of our controller for K o u r with different parameters.
Figure 18. Response of pendulum potential energy in the case of our controller for K o u r with different parameters.
Symmetry 13 01491 g018
Figure 19. Swing up of the RIP and stabilisation frame.
Figure 19. Swing up of the RIP and stabilisation frame.
Symmetry 13 01491 g019
Table 1. Simulation performance analyses.
Table 1. Simulation performance analyses.
DefinitionUnitKour Cont.Cont. [1]LQR Cont.LQR Cont. [1]
Swing up and stabilisation times1.535.43.3515
Steady-state errorrad0.010.10.20.2
Controller saturation times0.0030.005913
Table 2. Experimental performance analyses.
Table 2. Experimental performance analyses.
DefinitionUnitData
Swing up and stabilisation times1.45
Extreme forceN2.23
Steady-state errorrad0.03
Controller saturation times0.001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bajrami, X.; Pajaziti, A.; Likaj, R.; Shala, A.; Berisha, R.; Bruqi, M. Control Theory Application for Swing Up and Stabilisation of Rotating Inverted Pendulum. Symmetry 2021, 13, 1491. https://doi.org/10.3390/sym13081491

AMA Style

Bajrami X, Pajaziti A, Likaj R, Shala A, Berisha R, Bruqi M. Control Theory Application for Swing Up and Stabilisation of Rotating Inverted Pendulum. Symmetry. 2021; 13(8):1491. https://doi.org/10.3390/sym13081491

Chicago/Turabian Style

Bajrami, Xhevahir, Arbnor Pajaziti, Ramë Likaj, Ahmet Shala, Rinor Berisha, and Mirlind Bruqi. 2021. "Control Theory Application for Swing Up and Stabilisation of Rotating Inverted Pendulum" Symmetry 13, no. 8: 1491. https://doi.org/10.3390/sym13081491

APA Style

Bajrami, X., Pajaziti, A., Likaj, R., Shala, A., Berisha, R., & Bruqi, M. (2021). Control Theory Application for Swing Up and Stabilisation of Rotating Inverted Pendulum. Symmetry, 13(8), 1491. https://doi.org/10.3390/sym13081491

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