Research Article  Open Access
Optimal Manoeuvres of Underactuated Linear Mechanical Systems: The Case of Controlling Gantry Crane Operations
Abstract
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 closedloop 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 pointtopoint 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” (zerovalued) 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 socalled matrix of constraints, which is utilized to handle the nonintegrable constraints with the help of timevarying 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 threeDOF gantry crane operations. The method is verified by recreating the closedloop control of the twoDOF gantry crane problem obtained in [1] via applying the classical technique and the openloop optimal control considered in [3].
2. Problem Formulation
2.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: where and are vectors of DOFs and activation forces, respectively, and , , and 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 independent actuation forces are to control number of DOFs. Matrix of dimensions assigns components of vector to particular DOFs and obviously is not invertible if . Clearly, the actuators via (1) control all DOFs of the system. For the purpose of analysis the DOFs can be divided into actuated () and unactuated () ones by rearranging these equations as follows: The bottom row represents the equations constraining the actuated and unactuated DOFs in the following form: The system can formally be converted to a fully actuated one by using (3) to explicitly determine vector in terms of (i.e, ), and then by substituting this vector to the top row of (2a) to obtain . Unless some matrices in (3) vanish, it is not generally possible, and therefore these constraints are considered as nonholonomic.
The control task for vector in (1) is to manoeuvre the system from an initial state to a final state described by the following boundary conditions (pointtopoint manoeuvre): 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 formSeveral methods have been proposed to solve the inverse problems of finding the input for the output as defined by (1) and (5a), notably the servoconstraint approach [6–9] and the flatness method [10, 11]. In particular, differentiating (5a) with outputs twice one obtains , where the size of matrix is ; then square matrix 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 (size ) is transformed to the equally sized vector of modal variables . Similarly, vector is related to an equally sized vector of modal controls . These transformations arewhere is the transfer matrix of size between vectors and and mode shape matrix relates vectors and . The normalized matrix , consisting of modal shape vectors (each with components), satisfies the following orthogonality conditions: where is the unitary matrix and is the diagonal matrix of ordered frequencies with the terms . Each mode shape vector and frequency are solutions to the eigenvalues problem (). 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 ).
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: where for the Rayleigh damping (i.e., ) the diagonal terms of are and where are the modal damping ratios. Note that a rigid body translation, for which , is also included in the above equation.
A continuous system, or an FE model (1) of the system described by DOFs (where may be a large number), can be approximated by (8) with only significant modes considered, where usually . 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 modes (similarly as for the system's DOFs) one can consider modes as actuated and modes as unactuated. Then matrix will be reduced to only columns, and transfer matrix in (6a) and (6b) will be of dimensions .
In order to control all modes this system can be made artificially fully actuated by adding dummy actuation forces (zero valued) forming subvector . For the purpose of analysis vector in (1) is replaced by the augmented force vector containing real actuation forces forming subvector and dummy actuators forming subvector . Then, in (1) and (6a) and (6b), matrix of dimensions is replaced by matrix of dimensions (this matrix assigns the component of to particular nodes). Consequently in (6b) matrix of dimensions is replaced by a square matrix of dimensions ( modes controlled by actuators). The dummy actuators should be placed in such a way that is nonsingular.
In the inverse dynamics based control analysis, each control can be obtained from (8) by substituting the corresponding prescribed mode . Then, for known vector , the actuation forces should be determined by inverting transformation (6b) in which matrix (instead of ) 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 , representing the mapping from modal controls to actuation forces for any augmented system of size , is partitioned as follows: Square submatrix is of size and square submatrix is of size . To be consistent with modes classifications (actuated and unactuated), vectors and are referred to as actuated and unactuated modal controls, respectively. Given the nullvalued dummy force vector (size ) the bottom rows of operation (9) (lower partition) define constraints on the system in terms of all modal controls, in the following form: Matrix (size ) defines the system constraints written algebraically in terms of modal controls. Since (10) is homogeneous matrix can be normalized such that the diagonal terms corresponding to controls are set to unity (i.e., ). In this form becomes independent of the choice of dummy actuators, which reflects the fact that these zeroforce actuators were added somewhat arbitrarily only to facilitate the elimination process, that is, to satisfy the constraints in (9). Matrix 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 . They can also be obtained in terms of only actuated modal controls in vector by applying constraints (10) to eliminate unactuated modal controls . Thus, components of actuator forces in vector can be obtained in terms of actuated modal controls in vector from the following operation: Square matrix (size ) is referred to as the pseudotransfer matrix, and it relates actuated modal controls to real actuator forces. Similar to the normalized constraint matrix , the pseudotransfer matrix 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: where , , and are matrices, with the diagonal terms , , and (), that are weights for the system’s potential energy, kinetic energy, and actuator work, respectively. Note that 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 for the constrained optimization problem involving performance index (12), uncoupled equations of motion (8), and constraints (10) is defined in the following form: and are standard costate vectors related to modal position and velocity states ( and ) of a system, respectively. Vector represents the set of timedependent Lagrange multipliers introduced to enforce constraints (10). These multipliers play a similar role to, for example, that of the multipliers used in the servoconstraint approach [6–9] mentioned before. According to Pontryagin’s principle the costate equations take the following form:The Hamiltonian is stationary with respect to modal control if Substituting (8) into (15) gives Substituting (16) into (14b) yields Finally, substituting (17) into (14a) generates the following set of optimality equations: Note that optimality equations (18) contain unknown components in and unknown components in . Therefore, additional constraint equations (10) are required in order to obtain all the unknown modal variable functions in vector and Lagrange multiplier functions in vector . However, the constraints must be written in terms of not in terms of (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 The number of ((18) and (19)) is equal the unknown components in vectors and .
Boundary conditions (4) are mapped into modal space by using the inverse of transformation (6a) or through the relation (obtained by additional substitution of condition (7a)). These transformed boundary conditions are For fully actuated problems, the last term () in optimality equations (18) vanishes because there are no constraints or Lagrange multipliers needed to enforce them. Therefore, a fully actuated problem involves only optimality equations (18) to be solved in terms of uncoupled modal variables in vector .
The solution to the combined set of (18), (19), and (20) can be efficiently obtained using symbolic differential operator . Substituting this operator into (19) and (20) and rewriting in matrix notation give where Matrix contains submatrices , , and . Vector contains all unknown modal variables and Lagrange multipliers. Note that in a fully actuated case, matrix in (21) consists only of submatrix and vector .
The solution to a system described in form (21) involves the roots () of the characteristic equation for the determinant of [16], where operator is replaced by the auxiliary variable rendering a order polynomial. This operation is written as Generally, the roots of the characteristic equation (23) take the following form: The positive real numbers and characterize the response of the mode of motion. For nonzero, unique roots, solution vector consists of components that can be written in terms of independent elementary functions related to the roots (24), in the form [16]: Obviously, the frequency of mode controlled by the actuators can be interpreted as and its rate of active attenuation (or amplification) as . If multiple roots and zerovalued roots are obtained from (23), then solution functions (25) must be modified to mathematically accommodate these situations. There are unknown integration constants contained in the solution functions (25).
Integration constants are obtained by substituting the assumed form (25) into differential equations (18) and (19) and using the method of undetermined coefficients to generate sets of linear algebraic equations relating these constants. By replacing one set of equations with the set of 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 closedloop control, asymptotically convergent solution functions are required such that the control task is met over an infinite period of time (). The resulting number of integration constants is reduced by half, as terms involving positive exponential in the solution form (25) disappear ().
To quantitatively measure the performance of closedloop control schemes settling time is defined as the time needed for various variables to be reduced to within 3% of their initial value (i.e., ).
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 ). 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 cartdriving force (, ).
The gantry crane model is shown in Figure 1. The model includes the mass of the cart , the mass of the suspended load , swing angle , gravitational acceleration , horizontal distance from the cart’s initial position to the origin, and length 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 (, ) to a final resting equilibrium state at the origin (, ) by applying timevarying force . 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 is added to artificially make the system fully actuated and formulate the augmented 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 (closedloop) 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: To be consistent with the assumptions made in [1, 3] no dissipative effects (i.e., friction, etc.) are considered ().
The initial and final conditions (consistent with [1]) take the following forms: The modal analysis () givesThe rigid body translational mode of motion is represented in (28a) by the frequency and the second vibrating mode (load swinging) is represented by the frequency .
The uncoupled modal equations of motion (8) become: The augmented system transfer matrix is obtained by the appropriate substitutions from (26) and (28b) into the general partitioned form (9): Modal controls and are considered actuated and unactuated, respectively. The constraint equation is obtained by normalizing the bottom row of matrix in (30) to obtain The constraint (31) may be applied to eliminate redundant modal control from the top row operation of (30) to obtain force as a function of independent modal control , giving where is the pseudotransfer matrix (since , this matrix has only one term). Cartdriving force may be applied using openloop control (as a known function of time) or using closedloop control through a set of gains in fullstate 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 () and two modal controls (): The coupled optimality equations (18) take the following forms:where is the Lagrange multiplier used to meet the constraint (31). The differential form (19) of constraint equation (31) is written asIn modal space, the boundary conditions (27) areEquations (34a), (34b), and (35a) written according to form (21) (with ) yield where The characteristic equation of the system represented in (37) is obtained through operation (23), giving the 8th order polynomial equation: 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 (). This leaves twentyfour 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 ), eight algebraic equations are obtained for each differential equation in the set (34a),(34b) and (35a), resulting in a total of twentyfour equations in terms of twentyfour unknown integration constants . However, these twentyfour 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 (, ) to the origin (, ) will be derived for four cases using the CMSOC method. These cases are as follows(A)an openloop control that minimizes actuation forces for a fixed time interval as in [3];(B)a closedloop control that mimics the control presented in [1];(C)a closedloop control with response improved over that presented in [1];(D)a closedloop 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, , , , , and .
As a final case, the CMSOC method is applied to a modified threeDOF gantry crane, with an additional linkmass hinge attached to the existing model in Figure 1 and controlled by one or two actuators. This final case involves two subcases.(E1)A closedloop control that manoeuvres a modified gantry crane to the origin using the cartdriving force as well as a torque applied to the first rigid link (two actuators).(E2)An openloop control that manoeuvres the modified gantry crane to the origin using only the cartdriving force (one actuator) over a fixed time interval.
(A) OpenLoop 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 in an openloop scheme. The performance index is chosen to be consistent with that presented in [3], corresponding to the weightings in the general form (33) with all other weightings null valued. Thus, the optimal control minimizes Performance index (40) minimizes the modal controls or the actuation force over the finite manoeuvre time , which is chosen as to represent again one of the cases considered in [3]. The gantry crane’s characteristic polynomial equation (39) is simplified to The roots of (41) are . There are four zero roots , two imaginary roots , and two imaginary roots . When written in form (24), these roots correspond to and . Because of the zero roots and repeating roots, the solution functions take the following form: Each solution function (42) () contains eight unknown integration constants (), 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 Substituting (43) into (29) yields The Lagrange multiplier function (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 and ((43) and (44)) into DOFs and via transformation (6a), the trajectories shown in Figures 2(a) and 2(b) are obtained. Optimal force , shown in Figure 2(c), is obtained by substituting modal control (46) into transformation (32).
(a)
(b)
(c)
This phase of the solution was done automatically using MAPLE. The solution procedure accepts any problem with modes (obtained from FE analysis for more complex structures) controlled by actuators. The modaltoDOF transformations for the gantry crane are indicated in Figures 2(a), 2(b), and 2(c). As shown, the openloop 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 closedloop control presented in [1] are shown in Figure 3. This control requires an effective manoeuvre time of 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.
(a)
(b)
(c)
From Figure 2 and Figure 3, one can conclude that the openloop control performs the manoeuvre in a shorter period of time ( versus ) with much smaller peak force requirements (3.6 N versus 15 N) and much smaller angles of oscillation (16° versus 42°). Also, the openloop control brings the system to a complete stop after 4 s, while the closedloop control produces overshoot and the system takes longer to effectively come to rest.
Calculations show that if the finite manoeuvre time for the openloop control is extended (or shortened), the peak force requirement and maximum swing angle are reduced (or increased)—approximately proportional to . For example, if the openloop control is modified to settle over the same effective period of time as that of the closedloop control (), the maximum force is reduced to approximately 1.6 N with a maximum swing of about 7 deg.
The openloop 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. Closedloop 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 closedloop system that approximately produces the same dynamic responses as given in [1].
(B) ClosedLoop Control of Gantry Crane: Reproducing Control from [1]. A closedloop control can perform the same task as that of the openloop control (case (A)); however it does so automatically, without prior knowledge of initial conditions involved. Any disturbance from its resting configuration at the origin (, ) is relayed through a set of constant gains to generate the cartdriving force to attenuate this disturbance.
In general, to simulate the closedloop process analytically the manoeuvre time 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 openloop system of case (A), which contained zero roots and purely imaginary roots). It can be verified that the weightings and 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 , , , and . The resulting characteristic polynomial equation (39) has eight roots that take form (24), with real and complex parts equal to Note that the first actively controlled mode of frequency s^{−1} ( for uncontrolled system) is damped with the ratio , while the second mode of frequency s^{−1} ( s^{−1} for uncontrolled system) is damped with the ratio .
Similar to case (A), modal variables and 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 closedloop problem requires that only half as many integration constants must be solved because the coefficients preceding exponential growth functions () are assumed to be null valued. This gives 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 closedloop 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 fullstate feedback control the active force is a function of all system states in the following form: For the general CMSOC method, gains and correspond to the observed positions and velocities of all DOFs of a system. For the gantry crane, (49) takes the following form: By substituting known DOF trajectories ( and ) and the known force function () into (50) and grouping the terms related to the four independent elementary solution functions (operations are done in MAPLE automatically), one obtains: 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 Though initial conditions were assumed in determining the trajectories and and force , it can be verified that gains (52) remain invariant towards any choice of these assumed conditions.
The control gains used in [1] were Comparing the gains (52) and (53) confirms that the CMSOC method is able to closely reproduce the closedloop 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 closedloop 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) ClosedLoop 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 () to states and , 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.

Each case reflects a performance index which gives significant weightings to an incrementally increasing number of system states (of four possible states , , , ), while holding the weighting on both modal controls () at unity. Case P0 gives none of the states a significant weighting, case P1 gives a significant weighting to a single state (), case P2 gives significant weightings to two states (), 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 and 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 ( and ) and force histogram () for the manoeuvres minimizing the performance indices for cases P0–P4 are presented in Table 2. The settling times 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).

Note that the first modal variable primarily influences the cart’s rigid body mode of motion, while the second modal variable influences the suspended load rotation. In fact there is a direct relationship between the angle of the load rotation and the second modal variable () such that this angular trajectory is directly affected by varying the weights given to and its derivative 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 and its derivative .
The performance index in case P0 heavily weights the modal controls and in comparison to modal variables and (100 times more) and neglects the modal velocities and . 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 (). If the weightings , and were increased even further relative to the weighting , 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 , while other weightings remain unchanged from case P0. This control greatly increases the speed at which the cart reaches its target position at (~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 ().
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 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 (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 . 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 (), 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 (). 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 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 and 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 (), while its displacement is near minimum (). 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 (). Cases P2 and P3 failed to adequately weight , 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 closedloop control, the gains are obtained from (50) in a similar fashion as in case (B), giving Note that the gains and are somewhat close to the gains for the control presented in [1] (53), but gains and are substantially different.
(D) ClosedLoop 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 and 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 ), written as The optimality equations in the differential operator form (21) become where With weightings chosen according to the performance index in case P4 (, , , , and ), the roots of the characteristic equation (39) for the system given by (56) () take form (24) with the following real and imaginary parts: For any fully actuated system, each modal variable is independently controlled by a single modal control , resulting in uncoupled solution functions of the following form: For the gantry crane () the four unknown integration constants are obtained by substituting the four initial conditions for and 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 and the optimal forces on the cart and suspended load and , respectively. The angular trajectory of the load is not shown because it remains zero () 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 () to the origin is obtained; however the task requires relatively large maximum forces (104 N) compared to previous cases. The required actions of cartdriving force and suspended load guiding force 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 () 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 openloop control scheme.
(a)
(b)
(E1) ClosedLoop Control of Modified ThreeDOF Gantry Crane (Two Actuators). To illustrate the application of the CMSOC method to a problem of a higher dimension (), the gantry crane is modified by adding an additional link with an end load, as shown in Figure 5 (a case of a doublependulum 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 ()—the standard cartdriving force and torque , produced by a motor fixed to the cart and applied to the first rigid link of length which supports the mass . Dummy torque to be used in formulating the augmented system is applied to the second link of length which carries mass . All other physical variables are the same as in the original gantry crane model with the exception of and , which denote the angles of links of lengths and , respectively.
The standard matrices in the augmented system’s equation of motion (1) for this new model are The augmented system consists of DOF vector and force vector . The following numerical values are adopted: , , , and .
The set of 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: As before, the first mode represents the rigid body mode of motion (), while the second and third modes, with the squared frequencies and , represent the swinging modes of the rotating linkmasses.
Augmented force vector is related to modal control vector through the inverse of transformation (6b) which is partitioned according to (9) to give The constraint equation () is obtained from the bottom row of (62) () and normalized into the following form: The 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 and ) in the following form: Selecting a performance index of form (12) gives three () optimality equations in the form (18) that, with the constraint equation (63), may be written according to (21) in the following form: where The parameters in the equation above are the row and column components of the constraint matrix given by (63). The selected weightings for the performance index are , , , , and .
The twelve () 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 () and Lagrange multiplier (). Half of these roots generate exponential growth functions that are eliminated by assuming their corresponding integration constants to be zerovalued. Then through the method of undetermined coefficients, four () 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 twentyfour equations (the number of equations is ).
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 are also zero. In other words, the manoeuvre requires a horizontal cart translation from a resting position at with both links hanging vertically to the same resting position at the origin.
Figure 6 shows trajectories , , and of the threeDOF gantry crane as well as required actuation forces and .
(a)
(b)
(c)
The manoeuvre, requiring a maximum force of 29 N and a maximum torque of 18 Nm, is effectively completed after . 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) OpenLoop Control of Modified ThreeDOF 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—cartdriving force —is investigated for an openloop scheme. Both of the torque actuators and (Figure 5) are treated as dummy actuators and so the inverse transformation, while identical to (62), is repartitioned in the following form: The constraint equations () are obtained from the two bottom rows of (67) and normalized into the following form: According to (11) the single cartdriving force may be expressed in terms of the independent modal control (chosen as ) in the following form: Note that matrix is the same for all cases involving one actuator. Choosing a performance index in form (12) gives three () optimality equations in form (18) that, with constraint equation (68), may be written according to (21) in the following form: where The assumed weightings are and (only the control effort is to be minimized). Consistent with the openloop control presented in case (A) the finite manoeuvre time is chosen to be .
The solution procedure is similar to previous examples. Figure 7 shows trajectories , , and as well as the required cartdriving force . 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)
(b)
(c)
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 openloop control schemes as well as closedloop (asymptotically convergent) control schemes. In the latter case the weightings of the performance index can be translated to the gains of the fullstate feedback closedloop 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.
References
 I. Fantoni and R. Lozano, NonLinear Control for Underactuated Mechanical Systems, Springer, 2002. View at: Publisher Site
 L. Meirovitch and H. Baruh, “Control of selfadjoint distributedparameter systems,” Journal of Guidance, Control, and Dynamics, vol. 5, no. 1, pp. 60–66, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B. L. Karihaloo and R. D. Parbery, “Optimal control of a dynamical system representing a gantry crane,” Journal of Optimization Theory and Applications, vol. 36, no. 3, pp. 409–417, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. W. Spong, “The swing up control problem for the acrobot,” IEEE Control Systems Magazine, vol. 15, no. 1, pp. 49–55, 1995. View at: Publisher Site  Google Scholar
 R. A. Canfield and L. Meirovitch, “Integrated structural design and vibration suppression using independent modal space control,” AIAA Journal, vol. 32, no. 10, pp. 2053–2060, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. Seifried, “Two approaches for feedforward control and optimal design of underactuated multibody systems,” Multibody System Dynamics, vol. 27, no. 1, pp. 75–93, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Seifried and W. Blajer, “Analysis of servoconstraint problems for underactuated multibody systems,” Mechanical Sciences, vol. 4, pp. 113–129, 2013. View at: Publisher Site  Google Scholar
 W. Blajer and K. Kołodziejczyk, “Control of underactuated mechanical systems with servoconstraints,” Nonlinear Dynamics, vol. 50, no. 4, pp. 781–791, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. Blajer and K. Kołodziejczyk, “A geometric approach to solving problems of control constraints: theory and a DAE framework,” Multibody System Dynamics, vol. 11, no. 4, pp. 343–364, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Fliess, J. Lévine, P. Martin, and P. Rouchon, “Flatness and defect of nonlinear systems: introductory theory and examples,” International Journal of Control, vol. 61, no. 6, pp. 1327–1361, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. SiraRamirez and S. K. Agrawal, Differentially Flat Systems, Control Engineering Series, Marcel Dekker, 2004.
 K. J. Bathe, “Modal superposition,” in Finite Elements Procedures, chapter 9.3, Prentice Hall, Englewood Cliffs, NJ, USA, 1996. View at: Google Scholar
 S. Woods and W. Szyszkowski, “Simulating active vibration attenuation in underactuated spatial structures,” Journal of Guidance, Control, and Dynamics, vol. 32, no. 4, pp. 1246–1253, 2009. View at: Publisher Site  Google Scholar
 S. Woods and W. Szyszkowski, “Analysis and simulation of optimal vibration attenuation for underactuated mechanical systems,” AIAA Journal, vol. 47, no. 12, pp. 2821–2835, 2009. View at: Publisher Site  Google Scholar
 D. E. Kirk, Optimal Control Theory, Dover, New York, NY, USA, 1998.
 D. G. Zill and M. R. Cullen, Advanced Engineering Mathematics, PWS Publishing Company, Boston, Mass, USA, 1992.
 E. M. AbdelRahman, A. H. Nayfeh, and Z. N. Masoud, “Dynamics and control of cranes: a review,” Journal of Vibration and Control, vol. 9, no. 7, pp. 863–908, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H.H. Lee, “Modeling and control of a threedimensional overhead crane,” Journal of Dynamic Systems, Measurement and Control, vol. 120, no. 4, pp. 471–476, 1998. View at: Publisher Site  Google Scholar
 W. O. Connor and H. Habibi, “Gantry crane control of a doublependulum, distributedmass load, using mechanical wave concept,” Mechanical Sciences, vol. 4, pp. 251–261, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 S. Woods and W. Szyszkowski. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.