A Cart–Inverted-Pendulum System Consider the mechanical system shown in Figure 5.79, in which a uniform rod of mass M and length L is pivoted on a cart of mass m. An external force f is applied to the cart. Assume that the pendulum is constrained to move in a vertical plane, and the cart moves

Question

A Cart-Inverted-Pendulum System

Consider the mechanical system shown in Figure 5.79, in which a uniform rod of mass [latex]M[/latex] and length [latex]L[/latex] is pivoted on a cart of mass [latex]m[/latex]. An external force [latex]f[/latex] is applied to the cart. Assume that the pendulum is constrained to move in a vertical plane, and the cart moves without slipping along a horizontal line. Denote the displacement of the cart as [latex]x[/latex] and the angular displacement of the pendulum as [latex] heta[/latex].

a. Draw the necessary free-body diagram and kinematic diagram, and derive the nonlinear equations of motion.

b. Linearize the equations of motion for small angular motions, and determine the state-space form with [latex]x[/latex] and [latex] heta[/latex] as the outputs.

Accepted Answer

a. This is a mixed system with two masses. The motion of the cart is purely translational, and the inverted pendulum undergoes both translation and rotation. The free-body diagram and the kinematic diagram are given in Figure 5.80. Note that the acceleration of the mass center of the rod consists of three components. On the one hand, the pendulum moves together with the cart along the horizontal line at an acceleration of [latex]\ddot{x}[/latex]. On the other hand, it rotates about the pivot [latex]P[/latex] on the cart. The relative rotational accelerations include the tangential component [latex]\frac{L}{2} \ddot{ heta}[/latex] and the normal component [latex]\frac{L}{2} \dot{ heta}^{2}[/latex].

Applying the force equation to the whole system along the [latex]x[/latex] direction gives

[latex] \begin{gathered} + ightarrow x: \sum F_{x}=\sum_{i=1}^{2} m_{i}\left(a_{C_{i}} ight)_{x^{\prime}} \ f=m \ddot{x}+M \ddot{x}-M \frac{L}{2} \dot{ heta}^{2} \cdot \sin heta+M \frac{L}{2} \ddot{ heta} \cdot \cos heta . \end{gathered} [/latex]

The forces at the pivot are canceled out because they are internal forces between the cart and the pendulum. Applying the moment equation to the pendulum about point [latex]\mathrm{P}[/latex] results in

[latex] \begin{gathered} +\curvearrowright: \sum M_{\mathrm{P}}=I_{\mathrm{C}} \alpha+M_{\mathrm{eff}_{-} m \mathbf{a}_{C^{\prime}}} \ M g \cdot \frac{L}{2} \sin heta=\frac{1}{12} M L^{2} \ddot{ heta}+M \ddot{x} \cdot \frac{L}{2} \cos heta+M \frac{L}{2} \ddot{ heta} \cdot \frac{L}{2} . \end{gathered} [/latex]

Rearranging the two equations into the standard input-output form

[latex] \begin{gathered} (m+M) \ddot{x}+\frac{1}{2} M L \ddot{ heta} \cos heta-\frac{1}{2} M L \dot{ heta}^{2} \sin heta=f, \ \frac{1}{3} M L^{2} \ddot{ heta}+\frac{1}{2} M L \ddot{x} \cos heta-\frac{1}{2} M g L \sin heta=0 . \end{gathered} [/latex]

b. For small angular motions, [latex]\cos heta \approx 1, \sin heta \approx heta, \dot{ heta}^{2} heta \approx 0[/latex] (see Section 4.6). The linearized equations are

[latex] \begin{gathered} (m+M) \ddot{x}+\frac{1}{2} M L \ddot{ heta}=f, \ \frac{1}{3} M L^{2} \ddot{ heta}+\frac{1}{2} M L \ddot{x}-\frac{1}{2} M g L heta=0 . \end{gathered} [/latex]

The state, the input, and the output are specified as

[latex] \mathbf{x}=\left\{\begin{array}{l} x_{1} \ x_{2} \ x_{3} \ x_{4} \end{array} ight\}=\left\{\begin{array}{c} x \ heta \ \dot{x} \ \dot{ heta} \end{array} ight\}, \quad u=f, \quad y=\left\{\begin{array}{l} x \ heta \end{array} ight\} . [/latex]

We then take the time derivative of each state variable. For the first two,

[latex] \begin{aligned} & \dot{x}_{1}=\dot{x}=x_{3}, \ & \dot{x}_{2}=\dot{ heta}=x_{4} . \end{aligned} [/latex]

Note that the two linearized equations are coupled. The derivatives [latex]\ddot{x}[/latex] and [latex]\ddot{ heta}[/latex] cannot be solved using only one of the equations. To find [latex]\dot{x}_{3}[/latex] and [latex]\dot{x}_{4}[/latex], we rewrite the two linearized differential equations as

[latex] \left[\begin{array}{cc} m+M & \frac{1}{2} M L \ \frac{1}{2} M L & \frac{1}{3} M L^{2} \end{array} ight]\left\{\begin{array}{c}\ddot{x}\\ddot{ heta }\end{array} ight\}=\left\{\begin{array}{c} f \ \frac{1}{2} M g L heta \end{array} ight\}, [/latex]

from which [latex]\ddot{x}[/latex] and [latex]\ddot{ heta}[/latex] can be solved using Cramer's rule:

[latex] \begin{aligned} & \ddot{x}=\frac{(1 / 3) M L^{2} f-(1 / 4) M^{2} L^{2} g heta}{(1 / 12) M L^{2}(M+4 m)}, \ & \ddot{ heta}=\frac{-(1 / 2) M L f+(M+m)(1 / 2) M g L heta}{(1 / 12) M L^{2}(M+4 m)} . \end{aligned} [/latex]

Simplifying the equations gives

[latex] \begin{aligned} & \dot{x}_{3}=\frac{4}{M+4 m} u-\frac{3 M g}{M+4 m} x_{2}, \ & \dot{x}_{4}=\frac{-6}{L(M+4 m)} u+\frac{6(M+m) g}{L(M+4 m)} x_{2} . \end{aligned} [/latex]

Thus, the state-space form is

[latex] \begin{aligned} \left\{\begin{array}{l} \dot{x}_{1} \ \dot{x}_{2} \ \dot{x}_{3} \ \dot{x}_{4} \end{array} ight\} & =\left[\begin{array}{cccc} 0 & 0 & 1 & 0 \ 0 & 0 & 0 & 1 \ 0 & -\frac{3 M g}{M+4 m} & 0 & 0 \ 0 & \frac{6(M+m) g}{L(M+4 m)} & 0 & 0 \end{array} ight]\left\{\begin{array}{l} x_{1} \ x_{2} \ x_{3} \ x_{4} \end{array} ight\}+\left[\begin{array}{c} 0 \ 0 \ \frac{4}{M+4 m} \ -\frac{6}{L(M+4 m)} \end{array} ight] u \ \mathbf{y} & =\left[\begin{array}{llll} 1 & 0 & 0 & 0 \ 0 & 1 & 0 & 0 \end{array} ight]\left\{\begin{array}{c} x_{1} \ x_{2} \ x_{3} \ x_{4} \end{array} ight\}+\left[\begin{array}{l} 0 \ 0 \end{array} ight] u . \end{aligned} [/latex]

Note that in Part (a), Equation 5.31

[latex]\sum{\bf M_{p}}=I_{\mathrm{p}}\alpha+m{\bf r}_{\mathrm{C/P}} imes{\bf a_{P}}.[/latex] (5.31)

instead of Equation 5.41

[latex]\sum M_{\mathrm{P}}=I_{\mathrm{C}} \alpha+M_{ ext {eff } \_\mathbf{ma}_{\mathrm{C}}}[/latex] (5.41)

can also be applied to the pendulum to derive the second equation of motion. The related free-body diagram and the kinematic diagram are shown in Figure 5.81. Summing all externally applied moments about the pivot [latex]\mathrm{P}[/latex], which moves together with the cart with an acceleration of [latex]\ddot{x}[/latex], gives

[latex] \begin{gathered} +\curvearrowright: \sum \boldsymbol{M}_{\mathrm{P}}=l_{\mathrm{P}} \boldsymbol{\alpha}+m \mathbf{r}_{\mathrm{C} / \mathrm{P}} imes \mathbf{a}_{\mathrm{P}} \ M g \cdot \frac{L}{2} \sin heta=\frac{1}{3} M L^{2} \ddot{ heta}+M \frac{L}{2} \ddot{x} \sin \left(\frac{\pi}{2}- heta ight) \ \frac{1}{3} M L^{2} \ddot{ heta}+\frac{1}{2} M L \ddot{x} \cos heta-\frac{1}{2} M g L \sin heta=0 . \end{gathered} [/latex]

Question 5.13: A Cart–Inverted-Pendulum System Consider the mechanical syst...

Related Answered Questions

Verified Answer:

Verified Answer:

A Two-Degree-of-Freedom Mass–Spring System Consider the two-degree-of-freedom mass–spring system shown in Figure 5.106. The mathematical model of the system is given by a set of ordinary differential equations m1x1 + (k1 + k2)x1 – k2x2 = 0, m2x2 – k2x1 + k2x2 = 0, where m1 = m2 = 5 kg, ...

Verified Answer:

Verified Answer:

Verified Answer:

A Single-Degree-of-Freedom Gear–Train System For the gear–train system shown in Figure 5.95a, derive the differential equation of motion. The mass moments of inertia of the two gears about their respective fixed centers are IC1 and IC2. ...

Verified Answer:

Verified Answer:

Verified Answer:

Verified Answer:

A Pure Rolling Disk Consider the system shown in Figure 5.60a, in which a uniform disk of mass m and radius r rolls on a horizontal surface. A translational spring of stiffness k is attached to the disk. Assuming that there is no slipping between the disk and the surface, derive the differential ...

Verified Answer: