A method of solving optimal manoeuvre control of linear underactuated mechanical systems is presented. The nonintegrable constraints present in such systems are handled by adding dummy actuators and then by applying Lagrange multipliers to reduce their action to zero. The open- and closed-loop control schemes can be analyzed. The method, referred to as the constrained modal space optimal control (CMSOC), is illustrated in the examples of gantry crane operations.
1. Introduction
Underactuated mechanical systems have fewer independent actuators than degrees of freedom (DOFs) to be controlled [1]. Typical nonlinear examples of such systems, usually with only several DOFs, are rigid multilink robotic manipulators with passive joints or any manipulator with flexible links (described by at least one mode of vibration). Linear examples include vibrating structures with continuously distributed mass (i.e., with theoretically infinite number of DOFs to describe them) such as masts, antennas, buildings, brides, and car suspension, controlled by discrete actuators. This paper presents a method of analyzing and simulating optimal manoeuvres between two given configurations (often referred to as point-to-point manoeuvres) for linear underactuated systems. The method combines optimal control theory with computational mechanics and the finite element (FE) technique, in particular.
The number of DOFs equal to the number of actuators will be referred to as actuated (after [1]), while all remaining DOFs will be referred to as underactuated (however, all DOFs are in fact controlled). The actuated and unactuated DOFs must satisfy a number of constraints equal to the number of unactuated DOFs and resulting from the equations governing the motion of such systems. For mechanical systems we assume that these constraints may be nonintegrable (nonholonomic), meaning unactuated DOFs cannot be explicitly eliminated. Many of the techniques presented in the literature deal with underactuated problems by applying the constraints to eliminate the unactuated DOFs and then by solving the reduced fully actuated problems [2–4]. These approaches are limited to particular problems where the constraints can be simplified to a form making such mathematical manipulations possible. The method presented here is capable of dealing with any linear system, as it does not require the elimination of unactuated DOFs. Instead, the underactuated system is formulated as if it were fully actuated by adding “dummy” (zero-valued) actuators to all unactuated DOFs. The modal space is used in modelling the system motions. The method can be considered as an extension of the independent modal space control (IMSC; e.g., see [5]) into the underactuated problems, therefore it will be referred to as the constrained modal space optimal control (CMSOC) method. The system constraints resulting from underactuation are then determined by eliminating these dummy actuators. The constraints are algebraic in terms of controls but differential (nonintegrable) in terms of the DOFs. The algebraic form of the constraints is used to generate the so-called matrix of constraints, which is utilized to handle the nonintegrable constraints with the help of time-varying Lagrange multipliers. Pontryagin’s principle is used to optimize the trajectory and actuation forces.
This paper presents the CMSOC method in a general form and then explains some details of the corresponding numerical procedure in the examples of standard two or three-DOF gantry crane operations. The method is verified by recreating the closed-loop control of the two-DOF gantry crane problem obtained in [1] via applying the classical technique and the open-loop optimal control considered in [3].
2. Problem Formulation2.1. Dynamics of a General Underactuated System
The computational model for the motion of a linear mechanical system is represented by a standard form used in FE analysis:
(1)Mq¨+Cq˙+Kq=BFa,
where q and Fa are vectors of DOFs and activation forces, respectively, and M, C, and K are constant mass, damping, and stiffness matrices, respectively. In particular, (1) is suitable to model the dynamics of a range of actively controlled structural members undergoing small amplitude oscillations and finite translations. In underactuated systems na independent actuation forces are to control n>na number of DOFs. Matrix B of dimensions n×na assigns components of vector Fa to particular DOFs and obviously is not invertible if n≠na. Clearly, the actuators via (1) control all DOFs of the system. For the purpose of analysis the DOFs can be divided into actuated (qa) and unactuated (qr) ones by rearranging these equations as follows: (2a)[MaaMarMraMrr][q¨aq¨r]+[CaaCarCraCrr][q˙aq˙r]+[KaaKarKraKrr][qaqr]=[Fa0].The bottom row represents the equations constraining the actuated and unactuated DOFs in the following form:
(3)Mraq¨a+Craq˙a+Kraqa+Mrrq¨r+Crrq˙r+Krrqr=0.
The system can formally be converted to a fully actuated one by using (3) to explicitly determine vector qr in terms of qa (i.e, qr=g(qa)), and then by substituting this vector to the top row of (2a) to obtain Fa=Fa(qa,g(qa)). Unless some matrices in (3) vanish, it is not generally possible, and therefore these constraints are considered as nonholonomic.
The control task for vector Fa in (1) is to manoeuvre the system from an initial state to a final state described by the following boundary conditions (point-to-point manoeuvre):
(4)q(0)=q0,q˙(0)=q˙0,q(tf)=qf,q˙(tf)=q˙f.
It should be emphasised that no trajectory is specified in this task. A particular trajectory satisfying (1) and (4) may be determined if extra conditions are imposed on the system. We will identify such a trajectory by optimizing the performance index as discussed in the next section. Note that this problem is different from a typical trajectory tracking problem in which instead of (4) the task is specified as the system output in the form(5a)y=h(q).Several methods have been proposed to solve the inverse problems of finding the input Fa for the output y as defined by (1) and (5a), notably the servo-constraint approach [6–9] and the flatness method [10, 11]. In particular, differentiating (5a) with na outputs twice one obtains y¨=Hq¨+h-, where the size of matrix H is na×n; then square matrix HM-1B is required to be nonsingular to solve the problem. This condition does not apply in the method presented here since our output is given only in terms of (4), that is, the system’s initial and final configurations.
The set of (1) or (2a) is uncoupled when mapped into modal space, where vector of DOFs q (size n) is transformed to the equally sized vector of modal variables η=[η1⋯ηn]T. Similarly, vector F is related to an equally sized vector of modal controls U=[u1⋯un]T. These transformations are(6a)q=Φη,(6b)U=(ΦTB)Fa=B^Fa,where B^=ΦTB is the transfer matrix of size n×na between vectors Fa and U and mode shape matrix Φ=[φ1⋯φn] relates vectors q and η. The M-normalized matrix Φ, consisting of n modal shape vectors φi (each with n components), satisfies the following orthogonality conditions: (7a)ΦTMΦ=I,(7b)ΦTKΦ=Ω,where I is the unitary matrix and Ω is the diagonal matrix of ordered frequencies with the terms Ωii=ωi2. Each mode shape vector φi and frequency ωi are solutions to the eigenvalues problem (K-ωi2M)φi=0 (i=1,…,n). The above modal analysis (or operations defined by (6a)–(7b)) is carried out routinely in the FE approach, even for problems with a very large number of DOFs (large n).
The equations of motion (1) become uncoupled when applying transformations (6a) and (6b) subject to orthogonality conditions (7a) and (7b) and take the following form:
(8)Iη¨+Δη˙+Ωη=U,
where for the Rayleigh damping (i.e., C=αM+βK) the diagonal terms of Δ are Δii=2ςiωi=φiTCφi and where ςi=α/2ωi+βωi/2 are the modal damping ratios. Note that a rigid body translation, for which ωi=0, is also included in the above equation.
A continuous system, or an FE model (1) of the system described by n DOFs (where n may be a large number), can be approximated by (8) with only nm significant modes considered, where usually nm≪n. The number of significant modes that should be sufficient to represent such a system is generally problem related and depends mainly on its physical characteristics, the spatial distribution, and frequency content of the loading [12].
In the system approximated by nm modes (similarly as for the system's DOFs) one can consider na modes as actuated and nr=nm-na modes as unactuated. Then matrix Φ will be reduced to only nm columns, and transfer matrix B^ in (6a) and (6b) will be of dimensions nm×na.
In order to control all nm modes this system can be made artificially fully actuated by adding nr=nm-na dummy actuation forces (zero valued) forming subvector Fd. For the purpose of analysis vector Fa in (1) is replaced by the augmented force vector Fa'=[FaTFdT]T containing na real actuation forces forming subvector Fa and nr dummy actuators forming subvector Fd. Then, in (1) and (6a) and (6b), matrix B of dimensions n×na is replaced by matrix B′ of dimensions n×nm (this matrix assigns the component of Fa′ to particular nodes). Consequently in (6b) matrix B^ of dimensions n×na is replaced by a square matrix B^′ of dimensions nm×nm (nm modes controlled by nm actuators). The dummy actuators Fd should be placed in such a way that B^′ is nonsingular.
In the inverse dynamics based control analysis, each control Ui can be obtained from (8) by substituting the corresponding prescribed mode ηi. Then, for known vector U, the actuation forces Fa' should be determined by inverting transformation (6b) in which matrix B^′ (instead of B^) is now square and nonsingular (the dummy actuators were added to the system only to ensure that this inversion is possible). In the next step, after computing the inverse of operation (6b), the dummy actuators will be eliminated by giving them zero values. For that purpose the inverse matrix (B^′)-1, representing the mapping from modal controls U=[UaT∣UrT]T to actuation forces Fa′ for any augmented system of size nm×nm, is partitioned as follows:
(9)(B^′)-1U=Fa′⟹[B~aB~rAaAr][UaUr]=[FaFd]=[Fa0].
Square submatrix B~a is of size na×na and square submatrix Ar is of size nr×nr. To be consistent with modes classifications (actuated and unactuated), vectors Ua=[u1⋯una]T and Ur=[una+1⋯unm]T are referred to as actuated and unactuated modal controls, respectively. Given the null-valued dummy force vector Fd=[0⋯0]T (size nr×1) the bottom nr rows of operation (9) (lower partition) define constraints on the system in terms of all modal controls, in the following form:
(10)AU=AaUa+ArUr=0.
Matrix A=[Aa∣Ar] (size nr×nm) defines the system constraints written algebraically in terms of modal controls. Since (10) is homogeneous matrix A can be normalized such that the diagonal terms corresponding to controls Ua are set to unity (i.e., Aii=1). In this form A becomes independent of the choice of dummy actuators, which reflects the fact that these zero-force actuators were added somewhat arbitrarily only to facilitate the elimination process, that is, to satisfy the constraints in (9). Matrix A is discussed with more details in [13, 14].
Real actuation force(s) may be obtained from the top partition of operation (9) in terms of all modal controls in vector U. They can also be obtained in terms of only actuated modal controls in vector Ua by applying nr constraints (10) to eliminate unactuated modal controls Ur. Thus, na components of actuator forces in vector Fa can be obtained in terms of na actuated modal controls in vector Ua from the following operation:
(11)Fa=B-Ua.
Square matrix B-=B~a-B~rAr-1Aa (size na×na) is referred to as the pseudotransfer matrix, and it relates actuated modal controls to real actuator forces. Similar to the normalized constraint matrix A, the pseudotransfer matrix B- is independent of the choice of dummy actuators.
2.2. Optimal Manoeuvres of Underactuated Systems
In linear optimal control [15], the manoeuvre is optimal if, for a given task, it minimizes the performance index:
(12)J=12∫0tf((ηTQ^dη)+(η˙TQ^vη˙)+(UTR^U))dt⟶min,
where Q^d, Q^v, and R^ are matrices, with the diagonal terms Q^dii, Q^vii, and R^ii (i=1,…,nm), that are weights for the system’s potential energy, kinetic energy, and actuator work, respectively. Note that nm modal variables and modal controls are included in (12); however, these modes are not independent because of constraint (10), resulting from underactuation. Such a problem can be solved by applying Pontryagin’s principle. Here we use the procedure described in [15]. Hamiltonian H for the constrained optimization problem involving performance index (12), uncoupled equations of motion (8), and constraints (10) is defined in the following form:
(13)H=-12(ηTQ^dη+η˙TQ^vη˙+UTR^U)+Pdη˙+Pv(-Δη˙-Ωη+U)+vTAU.Pd and Pv are standard costate vectors related to modal position and velocity states (η and η˙) of a system, respectively. Vector vT=[v1⋯vnr] represents the set of time-dependent Lagrange multipliers introduced to enforce constraints (10). These multipliers play a similar role to, for example, that of the multipliers used in the servo-constraint approach [6–9] mentioned before. According to Pontryagin’s principle the costate equations take the following form:(14a)P˙d=-∂H∂η=Q^dη+ΩPv,(14b)P˙v=-∂H∂η˙=Q^vη˙-Pd+ΔPv.The Hamiltonian is stationary with respect to modal control if
(15)∂H∂U=-R^U+Pv+ATv=0.
Substituting (8) into (15) gives
(16)Pv=R^(Iη¨+Δη˙+Ωη)-ATv.
Substituting (16) into (14b) yields
(17)Pd=Q^vη˙-R^(Iη∴+Δη¨+Ωη˙)+ΔR^(Iη¨+Δη˙+Ωη)+ATv˙-ΔATv.
Finally, substituting (17) into (14a) generates the following set of optimality equations:
(18)R^η¨¨+(2ΩR^-Q^v-R^Δ2)η¨+(R^Ω2+Q^d)η-(ATv¨-ΔATv˙+ΩATv)=0.
Note that nm optimality equations (18) contain nm unknown components in η and nr unknown components in v. Therefore, additional nr constraint equations (10) are required in order to obtain all the unknown modal variable functions in vector η and Lagrange multiplier functions in vector v. However, the constraints must be written in terms of η not in terms of U (note the change in the constraints’ form from algebraic to differential). The uncoupled equations of motion (8) are substituted into algebraic constraints (10) to obtain
(19)A(Iη¨+Δη˙+Ωη)=0.
The number of nm+nr ((18) and (19)) is equal the unknown components in vectors η and v.
Boundary conditions (4) are mapped into modal space by using the inverse of transformation (6a) or through the relation η=ΦTMq (obtained by additional substitution of condition (7a)). These transformed boundary conditions are
(20)η(0)=ΦTMq0,η˙(0)=ΦTMq˙0,η(tf)=ΦTMqf,η˙(tf)=ΦTMq˙f.
For fully actuated problems, the last term (ATv¨-ΔATv˙+ΩATv) in optimality equations (18) vanishes because there are no constraints or Lagrange multipliers needed to enforce them. Therefore, a fully actuated problem involves only nm optimality equations (18) to be solved in terms of nm uncoupled modal variables in vector η.
The solution to the combined set of (18), (19), and (20) can be efficiently obtained using symbolic differential operator Dn=dn/dtn. Substituting this operator into (19) and (20) and rewriting in matrix notation give
(21)[E-E^TE~0][ην]=0orEpY=0,
where
(22)E=R^D4+(2R^Ω-Q^v-R^Δ2)D2+(R^Ω2+Q^d),E^=A(ID2+ΔD+Ω),E~=A(ID2-ΔD+Ω).
Matrix Ep contains submatrices E, E~, and -E^T. Vector Y=[ηT∣νT]T contains all unknown modal variables and Lagrange multipliers. Note that in a fully actuated case, matrix Ep in (21) consists only of submatrix E and vector Y=η.
The solution to a system described in form (21) involves the roots rl (l=1,…,4nm) of the characteristic equation for the determinant of Ep [16], where operator D is replaced by the auxiliary variable r rendering a 4nmth order polynomial. This operation is written as
(23)detEp|D→r=0.
Generally, the 4nm roots of the characteristic equation (23) take the following form:
(24)rl=±αk±iβk(k=1,…,nm,l=1,…,4nm).
The positive real numbers αk and βk characterize the response of the kth mode of motion. For nonzero, unique roots, solution vector Y consists of nm+nr components Yj that can be written in terms of 4nm independent elementary functions related to the roots (24), in the form [16]:
(25)Yj=∑k=1nm[eαkt(ckj1sin(βkt)+ckj2cos(βkt))+e-αkt(ckj3sin(βkt)+ckj4cos(βkt))eαkt(ckj1sin(βkt)+ckj2cos(βkt))],wherej=1,…,nm+nr.
Obviously, the frequency of kth mode controlled by the actuators can be interpreted as ωka=βk and its rate of active attenuation (or amplification) as ςka=αk/βk. If multiple roots and zero-valued roots are obtained from (23), then solution functions (25) must be modified to mathematically accommodate these situations. There are 4nm(nm+nr) unknown integration constants ckj1,…,ckj4 contained in the solution functions (25).
Integration constants ckj1,…,ckj4 are obtained by substituting the assumed form (25) into differential equations (18) and (19) and using the method of undetermined coefficients to generate nm+nr sets of 4nm linear algebraic equations relating these constants. By replacing one set of 4nm equations with the set of 4nm boundary conditions (20), the integration constants can be solved simultaneously. All these symbolic operations, including the determination of the roots (24) and constants in (25), can be done automatically using the MAPLE mathematical software.
For closed-loop control, asymptotically convergent solution functions are required such that the control task is met over an infinite period of time (tf→∞). The resulting number of integration constants is reduced by half, as terms involving positive exponential eαkt in the solution form (25) disappear (ckj1=ckj2=0).
To quantitatively measure the performance of closed-loop control schemes settling time tf3% is defined as the time needed for various variables to be reduced to within 3% of their initial value (i.e., e-αkmintf3%=0.03→tf3%=3.5/αkmin).
The above procedure was applied to actively suppress vibrations of a spatial antenna mast in [13] and of plane frames in [14], the cases with the number of DOFs much greater than the number of significant modes included in the analysis (i.e., with n≫nm). In both cases only oscillating modes were controlled. Here the application of the above methodology is focused on various control schemes, which are demonstrated in controlling the translational and oscillating modes of a gantry crane.
3. Dynamics and Optimal Control of the Gantry Crane System
The gantry crane problem is one of the simplest underactuated mechanical systems involving two DOFs—cart translation and suspended load rotation—and a single actuator—a cart-driving force (nm=2, na=nr=1).
The gantry crane model is shown in Figure 1. The model includes the mass of the cart M, the mass of the suspended load m, swing angle θ, gravitational acceleration g, horizontal distance a from the cart’s initial position to the origin, and length L of the massless rigid link connecting the cart and load. The task is to manoeuvre the system from an initial resting state at some nonzero horizontal distance (x=a, θ=0) to a final resting equilibrium state at the origin (x=0, θ=0) by applying time-varying force Fa. Any finite cart translations are permitted, but swings of the suspended load are assumed to be sufficiently small for a linearized model to be valid. In modal space rigid body translation for such a manoeuvre is easily separated from the oscillatory motion of the suspended load. Dummy force Fd is added to artificially make the system fully actuated and formulate the augmented gantry crane system.
Gantry crane system.
This same gantry crane model was used in several papers dealing with control or/and optimization. Notably, a Lyapunov function was used in [1] to obtain an asymptotically stable (closed-loop) control (linear and nonlinear) for attenuating disturbances (nonzero initial positions) in the system, and optimal control by applying Pontryagin's principle was considered in [3]. Results for the linearized system are of interest because they serve as a useful comparison for the controls obtained in this paper. Similar problems of controlling the plane motion of gantry cranes were presented in [10, 11] using the concept of flatness. Various aspects of controlling gantry cranes, 3D operations were considered in [9, 17–19].
The gantry crane system shown in Figure 1 and its coordinate system are chosen to mimic those used in [1]. Matrices and vectors in the general equation of motion (1) take the following forms:
(26)M=[M+m-m-mm],K=[000mgL],B=[10],B′=[110-1],q=[xLθ],Fa'=[FaFd].
To be consistent with the assumptions made in [1, 3] no dissipative effects (i.e., friction, etc.) are considered (C=0).
The initial and final conditions (consistent with [1]) take the following forms:
(27)q(0)=[a0]T,q˙(0)=q(tf)=q˙(tf)=[00]T.
The modal analysis (nm=2) gives(28a)Ω=[000(1+mM)gL],(28b)Φ=[1M+mmM(M+m)0M+mMm].The rigid body translational mode of motion is represented in (28a) by the frequency ω1=0 and the second vibrating mode (load swinging) is represented by the frequency ω2=(1+m/M)(g/L).
The uncoupled modal equations of motion (8) become:
(29)η¨1=u1,η¨2+ω22η2=u2.
The augmented system transfer matrix B^′=ΦTB′ is obtained by the appropriate substitutions from (26) and (28b) into the general partitioned form (9):
(30)[MM+mMmM+mmM+m-MmM+m][u1u2]=[FaFd]=[Fa0].
Modal controls u1 and u2 are considered actuated and unactuated, respectively. The (nr=1) constraint equation is obtained by normalizing the bottom row of matrix (B^′)-1 in (30) to obtain
(31)AU=[1-Mm][u1u2]=u1-Mmu2=0.
The constraint (31) may be applied to eliminate redundant modal control u2 from the top row operation of (30) to obtain force Fa as a function of independent modal control u1, giving
(32)Fa=B-Ua=(M+m)u1,
where B-=M+m is the pseudotransfer matrix (since na=1, this matrix has only one term). Cart-driving force Fa may be applied using open-loop control (as a known function of time) or using closed-loop control through a set of gains in full-state feedback. Both schemes will be analyzed and simulated using the CMSOC method.
The performance index (12) takes the following form, consisting of the gantry crane system’s four states (η1,η2,η˙1,η˙2) and two modal controls (u1,u2):
(33)J=12∫0tf(Q^d11η12+Q^d22η22+Q^v11η˙12+Q^v22η˙22+R^11u12+R^22u22)dt⟶min.
The nm=2 coupled optimality equations (18) take the following forms:(34a)R^11η1¨¨-Q^v11η¨1+Q^d11η1-ν¨=0,(34b)R^22η2¨¨+(2R^22ω22-Q^v22)η¨2+(R^22ω24+Q^d22)η2+Mm(ν¨+ω22ν)=0,where ν is the Lagrange multiplier used to meet the nr=1 constraint (31). The differential form (19) of constraint equation (31) is written as(35a)η¨1-Mm(η¨2+ω22η2)=0.In modal space, the boundary conditions (27) are(36a)η1(0)=aM+m,η1(tf)=0,η2(0)=η˙1(0)=η˙2(0)=0,η2(tf)=η˙1(tf)=η˙2(tf)=0.Equations (34a), (34b), and (35a) written according to form (21) (with Dn=dn/dtn) yield
(37)EpY=[E10-E^110E2-E^21E~11E~210][η1η2ν]=0,
where
(38)E1=R^11D4-Q^v11D2+Q^d11,E2=R^22D4+(2R^22ω22-Q^v22)D2+(R^22ω24+Q^d22),E^11=E~11=D2,E^21=E~21=-Mm(D2+ω22).
The characteristic equation of the system represented in (37) is obtained through operation (23), giving the 8th order polynomial equation:
(39)detEp|D→r=E1E^212+E2E^112|D→r=Mm(R^11r4-Q^v11r2+Q^d11)(r2+ω22)2+r4((R^22ω24+Q^d22)R^22r4+(2R^22ω22-Q^v22)r2+(R^22ω24+Q^d22))=0.
Eight roots may be obtained from the characteristic equation (39), which are then substituted into an appropriate assumed solution form (if the roots take the full complex form (24), then the assumed function takes form (25)) to characterize the three unknown solution functions (η1,η2,ν). This leaves twenty-four unknown integration constants to be determined by substituting the appropriate solution form into the three equations (34a), (34b), and (35a). By relating the coefficients corresponding to each of the eight independent elementary functions (i.e., in (25) each is in the form e(±αi±βi)t), eight algebraic equations are obtained for each differential equation in the set (34a),(34b) and (35a), resulting in a total of twenty-four equations in terms of twenty-four unknown integration constants ckji. However, these twenty-four equations are linearly dependent. To obtain a unique solution, any one set of eight algebraic equations (obtained from either (34a), (34b), or (35a)) must be replaced with the eight boundary conditions (36a).
The optimal actuation forces needed to drive the gantry crane from an initially disturbed position (x=a, θ=0) to the origin (x=0, θ=0) will be derived for four cases using the CMSOC method. These cases are as follows
an open-loop control that minimizes actuation forces for a fixed time interval as in [3];
a closed-loop control that mimics the control presented in [1];
a closed-loop control with response improved over that presented in [1];
a closed-loop control of the fully actuated system (two actuators).
For each case, the gantry crane’s physical parameters are chosen to match those given in [1]; namely, M=m=1kg, L=1m, g=9.8m/s2, a=-5m, and ω2=4.43s-1.
As a final case, the CMSOC method is applied to a modified three-DOF gantry crane, with an additional link-mass hinge attached to the existing model in Figure 1 and controlled by one or two actuators. This final case involves two subcases.
A closed-loop control that manoeuvres a modified gantry crane to the origin using the cart-driving force as well as a torque applied to the first rigid link (two actuators).
An open-loop control that manoeuvres the modified gantry crane to the origin using only the cart-driving force (one actuator) over a fixed time interval.
(A) Open-Loop Control of Gantry Crane Manoeuvre in a Finite Time Interval. The first control manoeuvres the gantry crane from a known initial position to the origin in a finite time interval tf in an open-loop scheme. The performance index is chosen to be consistent with that presented in [3], corresponding to the weightings R^11=R^22=1 in the general form (33) with all other weightings null valued. Thus, the optimal control minimizes
(40)J=∫0tf(u12+u22)dt=1M∫0tfFa2dt⟶min.
Performance index (40) minimizes the modal controls or the actuation force over the finite manoeuvre time tf, which is chosen as tf=4s to represent again one of the cases considered in [3]. The gantry crane’s characteristic polynomial equation (39) is simplified to
(41)(1+Mm)(r2+ω22)2r4=0.
The roots of (41) are r1,…,8=0,0,0,0,±iω2,±iω2. There are four zero roots r1,…,4=0, two imaginary roots r5,7=iω2, and two imaginary roots r6,8=-iω2. When written in form (24), these roots correspond to β2=ω2 and α1=β1=α2=0. Because of the zero roots and repeating roots, the solution functions take the following form:
(42)Yj=c1j+c2jt+c3jt2+c4jt3+(c5j+c7jt)sin(ω2t)+(c6j+c8jt)cos(ω2t).
Each solution function (42) (j=1,…,3) contains eight unknown integration constants ckj (k=1,…,8), which are determined through substitution and comparison of similar terms in any two differential equations in the set (34a), (34b), and (35a) and by substitution of the eight boundary conditions (36a). With the integration constants determined, the resulting solution functions are
(43)η1=-7.09-.104t+1.41t2-.235t3+.0235sin(4.43t)+.0151cos(4.43t),(44)η2=.144-.0719t+(.00451-.0333t)sin(4.43t)+(-.144+.0520t)cos(4.43t),(45)ν=-2.82+1.41t-.460sin(4.43t)-.295cos(4.43t).
Substituting (43) into (29) yields
(46)u1=2.82-1.41t-.460sin(4.43t)-.295cos(4.43t).
The Lagrange multiplier function v(t) (45), which represents a “modal force” enforcing the modal constraints, is not used in further analysis and is shown here only for completeness of the solution.
Mapping modal variables η1 and η2 ((43) and (44)) into DOFs x and θ via transformation (6a), the trajectories shown in Figures 2(a) and 2(b) are obtained. Optimal force Fa, shown in Figure 2(c), is obtained by substituting modal control u1 (46) into transformation (32).
Histograms of (a) cart trajectory x, (b) swing angle θ, and (c) force Fa (open loop).
x=12(η1+η2)
θ=2η2
Fa=2u1
This phase of the solution was done automatically using MAPLE. The solution procedure accepts any problem with nm modes (obtained from FE analysis for more complex structures) controlled by na≤nm actuators. The modal-to-DOF transformations for the gantry crane are indicated in Figures 2(a), 2(b), and 2(c). As shown, the open-loop control is able to perform the task in exactly four seconds, with a peak force of about 3.6 N and a maximum load swing angle of about 0.28 rad (16°). The optimal force accelerates the gantry crane over the first half of the manoeuvre (2 s) and decelerates the cart over the last half with identical, but opposite and mirrored, forces.
Similar plots for the closed-loop control presented in [1] are shown in Figure 3. This control requires an effective manoeuvre time of tf3%≈6s to reach the origin, a maximum load rotation angle of 0.73 rad (42 deg), and a maximum force of 15 N. It should be noted that this relatively large rotation angle is mentioned here (and other angles quoted in the sequel) for the purpose of comparison only.
Histograms of (a) cart trajectory, (b) swing angle, and (c) optimal force from [1].
From Figure 2 and Figure 3, one can conclude that the open-loop control performs the manoeuvre in a shorter period of time (tf=4s versus tf3%≈6s) with much smaller peak force requirements (3.6 N versus 15 N) and much smaller angles of oscillation (16° versus 42°). Also, the open-loop control brings the system to a complete stop after 4 s, while the closed-loop control produces overshoot and the system takes longer to effectively come to rest.
Calculations show that if the finite manoeuvre time for the open-loop control is extended (or shortened), the peak force requirement and maximum swing angle are reduced (or increased)—approximately proportional to tf-2. For example, if the open-loop control is modified to settle over the same effective period of time as that of the closed-loop control (tf≈6), the maximum force is reduced to approximately 1.6 N with a maximum swing of about 7 deg.
The open-loop control can always provide a faster and more efficient manoeuvre. However, it is possible only when the initial positions and manoeuvre times are known in advance. Closed-loop control is necessary if any initial configuration (unknown explicitly) is treated as disturbance, and its automatic reduction/removal is desired (the final position is at rest). Case (B) demonstrates how the CMSOC method is applied to analyze and simulate a closed-loop system that approximately produces the same dynamic responses as given in [1].
(B) Closed-Loop Control of Gantry Crane: Reproducing Control from [1]. A closed-loop control can perform the same task as that of the open-loop control (case (A)); however it does so automatically, without prior knowledge of initial conditions involved. Any disturbance from its resting configuration at the origin (x=0, θ=0) is relayed through a set of constant gains to generate the cart-driving force Fa to attenuate this disturbance.
In general, to simulate the closed-loop process analytically the manoeuvre time tf is infinite and all parameters are driven asymptotically to the origin with increasing time. For the gantry crane, this requires that all roots of the characteristic equation (39) be nonzero complex numbers in the left half of the complex plane (unlike the open-loop system of case (A), which contained zero roots and purely imaginary roots). It can be verified that the weightings Q^d11 and Q^d22 in the performance index (33) must be nonzero in order to meet these criteria.
The gantry crane control as given in [1] is closely reproduced by choosing the weightings in the performance index (33) equal to Q^v11=Q^v22=0, Q^d11=4.5, Q^d22=42, and R^11=R^22=1. The resulting characteristic polynomial equation (39) has eight roots that take form (24), with real and complex parts equal to
(47)α1=0.853,β1=0.856,α2=0.513,β2=4.46.
Note that the first actively controlled mode of frequency ω1a=0.856 s^{−1} (ω1=0 for uncontrolled system) is damped with the ratio ς1a=0.996, while the second mode of frequency ω2a=4.46 s^{−1} (ω2=4.43 s^{−1} for uncontrolled system) is damped with the ratio ς1a=0.115.
Similar to case (A), modal variables η1 and η2 are determined by substituting the parameters from (47) into the assumed solution function (25) and then solving for the unknown coefficients by comparing similar terms in two of the three optimality/constraint equations (34a), (34b) and (35a), and substituting the boundary conditions (36a). Unlike case (A), the closed-loop problem requires that only half as many integration constants must be solved because the coefficients preceding exponential growth functions (eαit) are assumed to be null valued. This gives
(48)η1=e-α1t(-7.70sin(β1t)-7.08cos(β1t))+e-α2t(.126sin(β2t)-.00564cos(β2t)),η2=e-α1t(-.565sin(β1t)+.534cos(β1t))+e-α2t(.149sin(β2t)-.534cos(β2t)),u1=e-α1t(-10.3sin(β1t)+11.3cos(β1t))+e-α2t(-2.44sin(β2t)-.686cos(β2t)).
Using the appropriate transformations (see Figure 2) the modal space variables (48) are mapped into the DOF space variables. The resulting system trajectories and the optimal force histogram are visually indistinguishable from those shown in Figure 3.
The CMSOC method can also generate the closed-loop gains from the assumed weighting coefficients to demonstrate that the gains corresponding to the solution (48) are obtained and compared with the gains used in [1].
In full-state feedback control the active force is a function of all system states in the following form:
(49)Fa=-Gdq-Gvq˙.
For the general CMSOC method, gains Gd=[g1d⋯gnmd] and Gv=[g1v⋯gnmv] correspond to the observed positions and velocities of all nm DOFs of a system. For the gantry crane, (49) takes the following form:
(50)Fa=-g1dx-g2dLθ-g1vx˙-g2vLθ˙.
By substituting known DOF trajectories (x=(1/2)(η1+η2) and θ=2η2) and the known force function (Fa=2u1) into (50) and grouping the terms related to the four independent elementary solution functions (operations are done in MAPLE automatically), one obtains: (51)e-α1t{(-14.6-5.85g1d+8.95g1v-.799g2d+.0352g2v)sin(β1t)+(16.0-4.63g1d-1.06g1v+.755g2d-1.33g2v)cos(β1t)}+e-α2t{(-3.45+.195g1d+1.56g1v+.211g2d+3.26g2v)sin(β2t)+(-.970-.374g1d+1.06g1v-.755g2d+1.33g2v)cos(β2t)}=0.Each of the bracketed terms in (51) (containing the unknown gains) must equal to zero for the equation to be true at any time. This gives four equations in terms of four unknown gains, which may be solved to obtain
(52)Gd=[g1dg2d]=[3.00.732],Gv=[g1vg2v]=[3.66-.924].
Though initial conditions were assumed in determining the trajectories x and θ and force Fa, it can be verified that gains (52) remain invariant towards any choice of these assumed conditions.
The control gains used in [1] were
(53)Gd*=[g1dg2d]=[3.0.71],Gv*=[g1dg2d]=[3.69-.87].
Comparing the gains (52) and (53) confirms that the CMSOC method is able to closely reproduce the closed-loop control in [1] by careful selection of the weighting parameters in performance index (33). However, as shown next in case (C), the performance of this closed-loop control may be improved through better selection of these weighting parameters to produce faster convergence without an increase in the required peak actuation forces.
(C) Closed-Loop Control of Gantry Crane: Improving Performance. Case (B) developed a control that closely reproduced the control given in [1] by minimizing a performance index that gave no weight (Q^v11=Q^v22=0) to states η˙1 and η˙2, representing the gantry crane’s velocity. Consequently, the control caused the gantry to gain too much speed and then overshoot its target and produce large persistent load swings. These problems are mitigated by a more careful choice of the performance index weighting parameters in (33). To demonstrate the effect these parameters have on the gantry crane’s dynamics and to illustrate how they might be meaningfully selected, several cases, labelled P1 to P5 (each with different performance indices as listed in Table 1), are considered.
Weightings for five different performance indices in form (33).
Index #
Q^d11
Q^d22
Q^v11
Q^v22
R^11
R^22
P0
0.01
0.01
0
0
1
1
P1
6
0.01
0
0
1
1
P2
6
50
0
0
1
1
P3
6
50
4
0
1
1
P4
6
50
4
50
1
1
Each case reflects a performance index which gives significant weightings to an incrementally increasing number of system states (of four possible states η1, η2, η˙1, η˙2), while holding the weighting on both modal controls (u1,u2) at unity. Case P0 gives none of the states a significant weighting, case P1 gives a significant weighting to a single state (η1), case P2 gives significant weightings to two states (η1,η2), and so on until case P4 significantly weights all four states. Table 1 summarizes how these weightings are chosen for each case. Since the gantry crane’s asymptotic convergence mathematically requires that weightings Q^d11 and Q^d22 in the index (33) are nonzero, a small value (0.01) is used instead of zero in cases P0 and P1 to demonstrate how the system behaves when these weightings are negligible. The DOF trajectories (x and θ) and force histogram (Fa) for the manoeuvres minimizing the performance indices for cases P0–P4 are presented in Table 2. The settling times tf3% of the DOFs are also listed for each case. All plots in Table 2 are shown over the first 8 s of the manoeuvre period except for P0 (30 s).
DOF responses and force histograms that optimize performance indices P0–P4.
Index #
x (m)
θ (rad)
Fa (N)
P0
P1
P2
P3
P4
Note that the first modal variable η1 primarily influences the cart’s rigid body mode of motion, while the second modal variable η2 influences the suspended load rotation. In fact there is a direct relationship between the angle of the load rotation and the second modal variable (θ=2η2) such that this angular trajectory is directly affected by varying the weights given to η2(Q^d22) and its derivative η˙2(Q^v22) in the performance index (33). Likewise, the speed at which the cart can be made to reach its target is affected through the weightings given to η1(Q^d11) and its derivative η˙1(Q^v11).
The performance index in case P0 heavily weights the modal controls u1 and u2 in comparison to modal variables η1 and η2 (100 times more) and neglects the modal velocities η˙1 and η˙2. The resulting control requires a small peak force (0.7 N), producing small maximum load swing angles (0.06 rad or 3.4 deg), but requires a very long manoeuvre time to converge to the origin (tf3%≈440s). If the weightings R^11, and R^22 were increased even further relative to the weighting Q^d11, the maximum force requirements and angular rotations would become infinitesimal while the settling times would approach infinity.
In case P1 a significant weighting value is given to the first modal variable η1(Q^d11=6), while other weightings remain unchanged from case P0. This control greatly increases the speed at which the cart reaches its target position at x=0 (~2 s), but upon reaching this position the load undergoes large swing angles (1.0 rad or 57 deg) that persist for a very long time (tf3%≈440s).
The maximum force increases significantly (17.3 N) in comparison to case P0 because the rigid body cart motion requires much larger accelerations during the initial 2 s of the manoeuvre in order to quickly attenuate η1 due to its significant weighting value.
Case P2 improves the load swing attenuation, which was poorly dampened in case P1, by including a large weighting value to the second modal variable η2(Q^d22=50) (other weightings remain the same as in the previous case). The maximum load swing angle is reduced (0.8 rad or 46 deg) and the load swinging motion is damped much more quickly (~6.5 s). The cart translation requires similar accelerations and thus approximately the same maximum force (17.3 N) is needed. By inspection, one can see that case P2 produces similar behaviour to the control given in [1] shown in Figure 3 and likewise shares the problem of target overshoot and large persistent load swings.
Case P3 reduces the tendency of the cart to overshoot the target by also giving a significant weighting to the first modal velocity η˙1(Q^v11=4). However, large persistent load swings are still present, and so convergence is not significantly improved over that produced in case P2. The maximum required force (17.3 N) remains essentially unchanged, while the load swing angles are reduced slightly (0.75 rad or 43 deg).
The performance index in case P4 includes a large weight on the second modal velocity η˙2 (Q^v22=50), while keeping all other weightings unchanged from case P3. This produces a control that reduces the magnitude of load swing angles (0.45 rad or 25.8 deg) while attenuating the swinging motion more quickly (tf3%≈4.2s). The gantry crane performs the manoeuvre in essentially a single load swing cycle, with similar initial cart accelerations and thus maximum forces (17.3 N) as in previous cases. Case P4 produces faster convergence then previous cases because, from an optimal control perspective, it incorporates all of the gantry crane’s states in the minimization by assigning all weightings in the performance index (33) with significant numerical values.
It is essential that the second modal velocity η˙2 is given a significant weighting value to yield fast convergence, because the energy of the suspended load oscillates equally between potential and kinetic energy. Since potential energy and kinetic energy are proportional to the squares of displacement and velocity, respectively, both of the corresponding states η2 and η˙2 should carry a significant and approximately equal weight in the performance index (33). Without weighting the load swing velocity state, the control focuses on eliminating swing angles but not swing velocities. However, when the load is near the bottom of its swing, its velocity is near maximum (η˙2=2θ˙→max), while its displacement is near minimum (η2=2θ→min). Therefore, the optimal force derived without consideration for the load swing velocity is unable to eliminate any significant portion of the load swing energy when the load is near the bottom of its swing (θ→0). Cases P2 and P3 failed to adequately weight η˙2, resulting in larger, more persistent load swings than in case P4.
The control produced in case P4 provides a significant improvement over the control presented in [1], as it converges more quickly to the origin, while reducing load swing magnitudes, without any increase in the required maximum forces. To complete the design of this closed-loop control, the gains are obtained from (50) in a similar fashion as in case (B), giving
(54)Gd=[g1dg2d]=[3.469.10],Gv=[g1vg2v]=[5.431.79].
Note that the gains g1d and g1v are somewhat close to the gains for the control presented in [1] (53), but gains g2d and g2v are substantially different.
(D) Closed-Loop Control of Gantry Crane: Fully Actuated Control. The CMSOC method can also be applied to fully actuated systems. To illustrate this, consider the same gantry crane system, now with both actuators Fa and Fd acting as real actuators (no dummy actuator). This situation may arise practically when a person is employed to guide the suspended load while the cart performs its translations.
Since the problem is fully actuated there, are no additional constraints on the system motion and consequently no Lagrange multipliers needed to enforce them. The optimal forces can be solved by calculating the inverse dynamics directly from (6b), which takes form (30) (except Fd≠0), written as
(55)F=B^-1U⟹[FaFd]=[MM+mMmM+mmM+m-MmM+m][u1u2].
The optimality equations in the differential operator form (21) become
(56)EpY=[E][η]=0⟹[E100E2][η1η2]=0,
where
(57)E1=R^11D4-Q^v11D2+Q^d11,E2=R^22D4+(2R^22ω22-Q^v22)D2+(R^22ω24+Q^d22).
With weightings chosen according to the performance index in case P4 (Q^d11=6, Q^d22=50, Q^v11=4, Q^v22=50, and R^11=R^22=1), the roots of the characteristic equation (39) for the system given by (56) (E1E2|D→r=0) take form (24) with the following real and imaginary parts:
(58)α1=3.62,β1=2.78,α2=1.49,β2=0.474.
For any fully actuated system, each modal variable ηi is independently controlled by a single modal control ui, resulting in uncoupled solution functions of the following form:
(59)ηi=eαit(ci1sin(βit)+ci2cos(βit)).
For the gantry crane (i=1,2) the four unknown integration constants ci1,2 are obtained by substituting the four initial conditions for ηi(0) and η˙i(0) given by (36a). As in the previous cases, the solved modal variables in form (59) are mapped into the original coordinates to obtain the DOF trajectories and optimal forces. Figure 4 shows the cart trajectory x and the optimal forces on the cart and suspended load Fa and Fd, respectively. The angular trajectory θ of the load is not shown because it remains zero (θ=0) all the time. Practically, this means that for the optimal manoeuvre the person (actuator) guiding the suspended load must simply act to prevent it from swinging. Fast convergence (tf3%=0.78s) to the origin is obtained; however the task requires relatively large maximum forces (104 N) compared to previous cases. The required actions of cart-driving force Fa and suspended load guiding force Fd are identical, as the whole gantry crane system moves as a single rigid body. If smaller forces are desired, then a larger weight may be given to the modal controls (R^11,R^22) in the performance index. Note that only a 30 N maximum force would be required to execute the manoeuvre in 1 s by applying an open-loop control scheme.
Histograms of (a) cart trajectory and (b) forces for the fully actuated gantry crane.
(E1) Closed-Loop Control of Modified Three-DOF Gantry Crane (Two Actuators). To illustrate the application of the CMSOC method to a problem of a higher dimension (nm=3), the gantry crane is modified by adding an additional link with an end load, as shown in Figure 5 (a case of a double-pendulum gantry crane in 3D was presented in [19]). In comparison to the previous cases considered the control task is unchanged except that the oscillations of the additional suspended load must also be damped. Consider the control that uses two actuators (na=2)—the standard cart-driving force Fa and torque Ta, produced by a motor fixed to the cart and applied to the first rigid link of length L1 which supports the mass m1. Dummy torque Td to be used in formulating the augmented system is applied to the second link of length L2 which carries mass m2. All other physical variables are the same as in the original gantry crane model with the exception of θ1 and θ2, which denote the angles of links of lengths L1 and L2, respectively.
Modified three-DOF gantry crane model.
The standard matrices in the augmented system’s equation of motion (1) for this new model are
(60)M=[M+m1+m2-m1-m2-m2-m1-m2m1+m2m2-m2m2m2],K=[0000g(m1+m2)L1000gm2L2],B′=[10001L1-1L1001L2].
The augmented system consists of DOF vector q=[xL1θ1L2θ2]T and force vector Fa'=[FaTaTd]T. The following numerical values are adopted: M=m1=m2=1kg, L1=L2=1m, a=-5m, and g=9.8m/s2.
The set of nm=3 equations of motion (1) are uncoupled in modal space, with matrices of ordered frequencies Ω and mode shapes φ normalized according to (7a) and (7b), taking the following forms:
(61)Ω=[000012.400046.4],Φ=[.577.577.5770.3661.3701-1].
As before, the first mode represents the rigid body mode of motion (ω12=0), while the second and third modes, with the squared frequencies ω22=12.43(rad/s)2 and ω32=46.37(rad/s)2, represent the swinging modes of the rotating link-masses.
Augmented force vector Fa' is related to modal control vector U through the inverse of transformation (6b) which is partitioned according to (9) to give
(62)B^-1U=[1.7300-1.731.37.366-.577.789-.211][u1u2u3]=[FaTaTd].
The constraint equation (nr=1) is obtained from the bottom row of (62) (Td=0) and normalized into the following form:
(63)AU=[1-1.37.366][u1u2u3]=0.
The na=2 actuation forces may be obtained directly from the top two rows of (62) in terms of all modal controls, but according to (11) these forces may be expressed in terms of two independent modal controls (chosen as u1 and u2) in the following form:
(64)B-Ua=[1.730-2.732.73][u1u2]=[FaTa].
Selecting a performance index of form (12) gives three (nm=3) optimality equations in the form (18) that, with the constraint equation (63), may be written according to (21) in the following form:
(65)EpY=[E100E^110E20E^2100E3E^31E~11E~21E~310][η1η2η3ν]=0,
where
(66)Ei=R^iiD4-(2R^iiωi2-Q^vii)D2+(R^iiωi4+Q^dii),E^ij=Aji(D2+ωi2)(i=1,…,3,j=1).
The parameters Aji in the equation above are the jth row and ith column components of the constraint matrix A given by (63). The selected weightings for the performance index are Q^d11=6, Q^d22=Q^d33=50, Q^v11=4, Q^v22=Q^v33=50, and R11=R22=R33=1.
The twelve (4nm) roots of the characteristic equation (23) are obtained in the complex form (24) and are used to generate an assumed solution of form (25) for each unknown modal variable (nm=3) and Lagrange multiplier (nr=1). Half of these roots generate exponential growth functions that are eliminated by assuming their corresponding integration constants to be zero-valued. Then through the method of undetermined coefficients, four (nm+nr) sets of six linear algebraic equations are obtained. Replacing one set by the set of six initial conditions, the unknown integration constants are obtained by solving the set of twenty-four equations (the number of equations is 2nm(nm+nr)).
The boundary conditions for this problem are the same as those chosen for the original gantry crane, written in (36a), with the additional condition that the initial and final positions and velocities of the third modal variable η3 are also zero. In other words, the manoeuvre requires a horizontal cart translation from a resting position at x=a=-5m with both links hanging vertically to the same resting position at the origin.
Figure 6 shows trajectories x, θ1, and θ2 of the three-DOF gantry crane as well as required actuation forces Fa and Ta.
(a) Cart trajectory, (b) swing angles, and (c) force/torque for the modified gantry crane.
The manoeuvre, requiring a maximum force of 29 N and a maximum torque of 18 Nm, is effectively completed after teff3%=4s. The maximum load swing angle of the first link is 0.11 rad (6.3°) and that of the second link is 0.48 rad (27.5 deg).
(E2) Open-Loop Control of Modified Three-DOF Gantry Crane (One Actuator). In order to show the case where one actuator controls three DOFs, the optimal manoeuvre for the modified gantry crane using only a single actuator—cart-driving force Fa—is investigated for an open-loop scheme. Both of the torque actuators Ta and Td (Figure 5) are treated as dummy actuators and so the inverse transformation, while identical to (62), is repartitioned in the following form:
(67)B^-1U=[1.7300-1.731.37.366-.5774.789-.211][u1u2u3]=[FaTaTd]=[Fa00].
The constraint equations (nr=2) are obtained from the two bottom rows of (67) and normalized into the following form:
(68)AU=[1-10011][u1u2u3]=0.
According to (11) the single cart-driving force may be expressed in terms of the independent modal control (chosen as u1) in the following form:
(69)B-Ua=[1.73][u1]=[Fa].
Note that matrix B- is the same for all cases involving one actuator. Choosing a performance index in form (12) gives three (nm=3) optimality equations in form (18) that, with constraint equation (68), may be written according to (21) in the following form:
(70)EpY=[E100E^11E^120E20E^21E^2200E3E^31E^32E^11E^21E^3100E^12E^22E^3200][η1η2η3ν1ν2]=0,
where
(71)Ei=R^iiD4-(2R^iiωi2-Q^vii)D2+(R^iiωi4+Q^dii),E^ij=Aji(D2+ωi2)(i=1,…,3,j=2).
The assumed weightings are R11=R22=R33=1 and Q^d11=Q^d22=Q^d33=Q^v11=Q^v22=Q^v33=0 (only the control effort is to be minimized). Consistent with the open-loop control presented in case (A) the finite manoeuvre time is chosen to be tf=4s.
The solution procedure is similar to previous examples. Figure 7 shows trajectories x, θ1, and θ2 as well as the required cart-driving force Fa. The manoeuvre requires a peak force of 4.8 N and completes the task in exactly 4 s. The maximum load swing angle of the first link is 0.19 rad (11°) and that of the second link is 0.35 rad (20°).
(a) Cart trajectory, (b) swing angles, and (c) force for modified gantry crane (open loop).
4. Conclusions
The CMSOC methodology was presented as a means of solving linear underactuated (or fully actuated) control problems. The gantry crane problem was selected to illustrate in detail various operations required for different control methodologies. As demonstrated the method can be applied to open-loop control schemes as well as closed-loop (asymptotically convergent) control schemes. In the latter case the weightings of the performance index can be translated to the gains of the full-state feedback closed-loop controllers. The operations would be identical for any similar problems with larger numbers of modes and actuators.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
FantoniI.LozanoR.MeirovitchL.BaruhH.Control of self-adjoint distributed-parameter systemsKarihalooB. L.ParberyR. D.Optimal control of a dynamical system representing a gantry craneSpongM. W.The swing up control problem for the acrobotCanfieldR. A.MeirovitchL.Integrated structural design and vibration suppression using independent modal space controlSeifriedR.Two approaches for feedforward control and optimal design of underactuated multibody systemsSeifriedR.BlajerW.Analysis of servo-constraint problems for underactuated multibody systemsBlajerW.KołodziejczykK.Control of underactuated mechanical systems with servo-constraintsBlajerW.KołodziejczykK.A geometric approach to solving problems of control constraints: theory and a DAE frameworkFliessM.LévineJ.MartinP.RouchonP.Flatness and defect of non-linear systems: introductory theory and examplesSira-RamirezH.AgrawalS. K.BatheK. -J.Modal superpositionWoodsS.SzyszkowskiW.Simulating active vibration attenuation in underactuated spatial structuresWoodsS.SzyszkowskiW.Analysis and simulation of optimal vibration attenuation for underactuated mechanical systemsKirkD. E.ZillD. G.CullenM. R.Abdel-RahmanE. M.NayfehA. H.MasoudZ. N.Dynamics and control of cranes: a reviewLeeH.-H.Modeling and control of a three-dimensional overhead craneConnorW. O.HabibiH.Gantry crane control of a double-pendulum, distributed-mass load, using mechanical wave concept