European Journal of Scientific Research ISSN 1450-216X Vol.33 No.4 (2009), pp.630-641 © EuroJournals Publishing, Inc 2009 http://www.eurojournals.com/ejsr.htm Optimal Tracking with Sway Suppression Control for a Gantry Crane System M.A Ahmad Control and Instrumentation Research Group (COINS), Faculty of Electrical and Electronics Engineering Universiti Malaysia Pahang, Lebuhraya Tun Razak 26300, Kuantan, Pahang, Malaysia E-mail: mashraf@ump.edu.my Tel: +609-5492366; Fax: +609-5492377 R.M.T Raja Ismail Control and Instrumentation Research Group (COINS) Faculty of Electrical and Electronics Engineering Universiti Malaysia Pahang, Lebuhraya Tun Razak 26300, Kuantan, Pahang, Malaysia E-mail: rajamohd@ump.edu.my Tel: +609-5492366; Fax: +609-5492377 M.S Ramli Control and Instrumentation Research Group (COINS) Faculty of Electrical and Electronics Engineering Universiti Malaysia Pahang, Lebuhraya Tun Razak 26300, Kuantan, Pahang, Malaysia E-mail: syakirin@ump.edu.my Tel: +609-5492366; Fax: +609-5492377 N.M Abdul Ghani Control and Instrumentation Research Group (COINS) Faculty of Electrical and Electronics Engineering Universiti Malaysia Pahang, Lebuhraya Tun Razak 26300, Kuantan, Pahang, Malaysia E-mail: normaniha@ump.edu.my Tel: +609-5492366; Fax: +609-5492377 M.A Zawawi Control and Instrumentation Research Group (COINS) Faculty of Electrical and Electronics Engineering Universiti Malaysia Pahang, Lebuhraya Tun Razak 26300, Kuantan, Pahang, Malaysia E-mail: mohdanwar@ump.edu.my Tel: +609-5492366; Fax: +609-5492377 Optimal Tracking with Sway Suppression Control for a Gantry Crane System 631 Abstract This paper presents investigations into the development of hybrid control schemes for input tracking and anti-swaying control of a gantry crane system A nonlinear overhead gantry crane system is considered and the dynamic model of the system is derived using the Euler-Lagrange formulation To study the effectiveness of the controllers, initially a Linear Quadratic Regulator (LQR) control is developed for cart position control of gantry crane This is then extended to incorporate input shaper control schemes for anti-swaying control of the system The positive input shapers with the derivative effects are designed based on the properties of the system Simulation results of the response of the manipulator with the controllers are presented in time and frequency domains The performances of the hybrid control schemes are examined in terms of level of input tracking capability, swing angle reduction, time response specifications and robustness to parameters uncertainty in comparison to the LQR control The effects of derivative order of the input shaper on the performance of the system are investigated Finally, a comparative assessment of the control techniques is presented and discussed Keywords: Gantry crane, sway control, input shaping, LQR controller Introduction The main purpose of controlling a gantry crane is transporting the load as fast as possible without causing any excessive swing at the final position However, most of the common gantry crane results in a swing motion when payload is suddenly stopped after a fast motion [1] The swing motion can be reduced but will be time consuming Moreover, the gantry crane needs a skilful operator to control manually based on his or her experiences to stop the swing immediately at the right position The failure of controlling crane also might cause accident and may harm people and the surrounding Various attempts in controlling gantry cranes system based on open loop system were proposed For example, open loop time optimal strategies were applied to the crane by many researchers such as discussed in [2,3] They came out with poor results because open loop strategy is sensitive to the system parameters (e.g rope length) and could not compensate for wind disturbances Another open loop control strategies is input shaping [4,5,6] Input shaping is implemented in real time by convolving the command signal with an impulse sequence The process has the effect of placing zeros at the locations of the flexible poles of the original system An IIR filtering technique related to input shaping has been proposed for controlling suspended payloads [7] Input shaping has been shown to be effective for controlling oscillation of gantry cranes when the load does not undergo hoisting [8, 9] Experimental results also indicate that shaped commands can be of benefit when the load is hoisted during the motion [10] On the other hand, feedback control which is well known to be less sensitive to disturbances and parameter variations [11] is also adopted for controlling the gantry crane system Recent work on gantry crane control system was presented by Omar [1] The author had proposed proportionalderivative PD controllers for both position and anti-swing controls Furthermore, a fuzzy-based intelligent gantry crane system has been proposed [12] The proposed fuzzy logic controllers consist of position as well as anti-sway controllers However, most of the feedback control system proposed needs sensors for measuring the cart position as well as the load swing angle In addition, designing the swing angle measurement of the real gantry crane system, in particular, is not an easy task since there is a hoisting mechanism This paper presents investigations into the development of hybrid control schemes for input tracking and anti-swaying control of a gantry crane system Hybrid control schemes based on feedforward with LQR controllers are investigated To demonstrate the effectiveness of the proposed 632 M.A Ahmad, R.M.T Raja Ismail, M.S Ramli, N.M Abdul Ghani and M.A Zawawi control schemes, initially a LQR controller is developed for cart position control of gantry crane system This is then extended to incorporate the positive input shapers for swing control of the gantry crane This paper provides a comparative assessment of the performance of hybrid control schemes with different derivative order of input shapers Gantry Crane System The two-dimensional gantry crane system with its payload considered in this work is shown in Figure 1, where x is the horizontal position of the cart, l is the length of the rope, θ is the sway angle of the rope, M and m is the mass of the cart and payload respectively In this simulation, the cart and payload can be considered as point masses and are assumed to move in two-dimensional, x-y plane The tension force that may cause the hoisting rope elongate is also ignored In this study the length of the cart, l = 1.00 m, M = 2.49 kg, m = 1.00 kg and g = 9.81 m/s2 is considered Figure 1: Description of the gantry crane system Modelling of the Gantry Crane This section provides a brief description on the modelling of the gantry crane system, as a basis of a simulation environment for development and assessment of the active sway control techniques The Euler-Lagrange formulation is considered in characterizing the dynamic behaviour of the crane system incorporating payload Considering the motion of the gantry crane system on a two-dimensional plane, the kinetic energy of the system can thus be formulated as T= 1 Mx& + m( x& + l&2 + l θ& + x&l& sin θ + x&lθ& cos θ ) 2 (1) The potential energy of the beam can be formulated as U = − mgl cos θ (2) To obtain a closed-form dynamic model of the gantry crane, the energy expressions in (1) and (2) are used to formulate the Lagrangian L = T − U Let the generalized forces corresponding to the generalized displacements q = { x , θ } be F = { F x ,0} Using Lagrangian’s equation Optimal Tracking with Sway Suppression Control for a Gantry Crane System 633 ⎛ ∂L ⎞ ∂ L ⎜ ⎟− j = 1,2 = Fj ⎜ ∂q& j ⎟ ∂q j ⎝ ⎠ the equation of motion is obtained as below, d dt (3) Fx = ( M + m ) &x& + ml (θ&&cos θ − θ& sin θ ) + ml&θ& cos θ + m&l&sin θ lθ&& + 2l&θ& + &x& cos θ + g sin θ = (4) (5) In order to eliminate the nonlinearity equation in the system, a linear model of gantry crane system is obtained The linear model of the uncontrolled system can be represented in a state-space form as shown in equation (6) by assuming the change of rope and sway angle are very small x& = Ax + Bu y = Cx [ (6) ] T with the vector x = x θ x& θ& ⎡0 ⎢0 ⎢ mg A = ⎢0 M ⎢ ( M + m) g ⎢ − ⎢⎣ Ml and the matrices A and B are given by 0⎤ ⎡ ⎤ ⎢ ⎥ 1⎥⎥ ⎢ ⎥ ⎥, C = [1 0 0] , 0⎥ , B = ⎢ ⎢ M ⎥ ⎥ ⎢ ⎥ ⎥ 0⎥ ⎢⎣− Ml ⎥⎦ ⎦ D = [0] (7) Linear Quadratic Regulator (LQR) Control Scheme A more common approach in the control of manipulator systems involves the utilization linear quadratic regulator (LQR) design [13] Such an approach is adopted at this stage of the investigation here In order to design the LQR controller a linear state-space model of the gantry crane system was obtained by linearising the equations of motion of the system For a LTI system x& = Ax + Bu (8) the technique involves choosing a control law u = ψ (x) which stabilizes the origin (i.e., regulates x to zero) while minimizing the quadratic cost function ∞ J = ∫ x (t ) T Qx (t ) + u (t ) T Ru (t ) dt (9) where Q = Q T ≥ and R = R T > The term “linear-quadratic” refers to the linear system dynamics and the quadratic cost function The matrices Q and R are called the state and control penalty matrices, respectively If the components of Q are chosen large relative to those of R , then deviations of R from zero will be penalized heavily relative to deviations of u from zero On the other hand, if the components of R are large relative to those of Q then control effort will be more costly and the state will not converge to zero as quickly A famous and somewhat surprising result due to Kalman is that the control law which minimizes J always takes the form u = ψ ( x) = − Kx The optimal regulator for a LTI system with respect to the quadratic cost function above is always a linear control law With this observation in mind, the closed-loop system takes the form x& = ( A − BK ) x (10) and the cost function J takes the form ∞ J = ∫ x(t )T Qx(t ) + (− Kx(t ))T R(− Kx(t ))dt (11) 634 M.A Ahmad, R.M.T Raja Ismail, M.S Ramli, N.M Abdul Ghani and M.A Zawawi ∞ J = ∫ x(t )T (Q + K T RK ) x(t )dt (12) Assuming that the closed-loop system is internally stable, which is a fundamental requirement for any feedback controller, the following theorem allows the computation value of the cost function for a given control gain matrix K Input Shaping Control Schemes The design objectives of input shaping are to determine the amplitude and time locations of the impulses in order to reduce the detrimental effects of system flexibility These parameters are obtained from the natural frequencies and damping ratios of the system The input shaping process is illustrated in Figure The corresponding design relations for achieving a zero residual single-mode swaying of a system and to ensure that the shaped command input produces the same rigid body motion as the unshaped command yields a two-impulse sequence with parameter as (13) t1 = 0, t2 = π , A1 = , A2 = K ωd 1+ K 1+ K where − ζπ 1−ζ K =e , ωd = ωn − ζ ( ω n and ζ representing the natural frequency and damping ratio respectively) and tj and Aj are the time location and amplitude of impulse j respectively The robustness of the input shaper to errors in natural frequencies of the system can be increased by solving the derivatives of the system swaying equation This yields a four-impulse sequence with parameter as 3π π 2π , t4 = , t3 = t1 = 0, t2 = ωd A1 = A3 = ωd ωd 1 + 3K + 3K + K 3K 2 + 3K + 3K + K where K as in (13) , A2 = 3K + K + 3K + K , A4 = K3 (14) + 3K + 3K + K A1 * A2 Time Unshaped Input Input Shaper Amplitude Amplitude Figure 2: Illustration of input shaping technique Time Shaped input Implementation and Result In this investigation, hybrid control schemes for tracking capability and sway angle suppression of the gantry crane system are examined Initially, a Linear Quadratic Regulator (LQR) is designed This is then extended to incorporate input shaping scheme for control of sway angle of the hoisting rope The tracking performance of the Linear Quadratic Regulator applied to the gantry crane system was investigated by firstly setting the value of vector K and N which determines the feedback control Optimal Tracking with Sway Suppression Control for a Gantry Crane System 635 law and for elimination of steady state error capability respectively Using the lqr function in the Matlab, both vector K and N were set as K = [700 428.3885 3.8267 1.4066 ] and N = [700] The natural frequencies were obtained by exciting the gantry crane system with an unshaped reference input under LQR controller The input shapers were designed for pre-processing the trajectory reference input and applied to the system in a closed-loop configuration, as shown in Figure Figure 3: Block diagram of the hybrid control schemes configuration LQR Controller Desired input Input shaper Shaped input N + Gantry Crane System - Output responses K 6.1 LQR Controller In this work, the input is applied at the cart of the gantry crane The cart position of the gantry crane is required to follow a trajectory within the range of ± m The first three modes of swing angle frequencies of the system are considered, as these dominate the dynamic of the system The responses of the gantry crane system to the unshaped trajectory reference input were analyzed in time-domain and frequency domain (spectral density) These results were considered as the system response to the unshaped input under tracking capability and will be used to evaluate the performance of the input shaping techniques The steady-state cart position trajectory of +4 m for the gantry crane was achieved within the rise and settling times and overshoot of 1.372 s, 2.403 s and 0.20 % respectively It is noted that the cart reaches the required position from +4 m to -4 m within s, with little overshoot However, a noticeable amount of swing angle occurs during movement of the cart It is noted from the swing angle response with a maximum residual of ±1.4 rad Moreover, from the PSD of the swing angle response the swaying frequencies are dominated by the first three modes, which are obtained as 0.3925 Hz, 1.177 Hz and 2.06 Hz with magnitude of 33.02 dB, -9.929 dB and -22.86 dB respectively 6.2 Hybrid Controller In the case of hybrid control schemes, a ZV (two-impulse sequence) and ZVDD (four-impulse sequence) shapers were designed for three modes utilising the properties of the system With the exact natural frequencies of 0.3925 Hz, 1.177 Hz and 2.06 Hz, the time locations and amplitudes of the impulses were obtained by solving equations (13) and (14) For evaluation of robustness, input shapers with error in natural frequencies were also evaluated With the 30% error in natural frequency, the system swaying frequencies were considered at 0.5103 Hz, 1.5301 Hz and 2.678 Hz for the three modes of swaying frequencies Similarly, the amplitudes and time locations of the input shapers with 30% erroneous natural frequencies for both the ZV and ZVDD shapers were calculated For digital implementation of the input shapers, locations of the impulses were selected at the nearest sampling time 636 M.A Ahmad, R.M.T Raja Ismail, M.S Ramli, N.M Abdul Ghani and M.A Zawawi The system responses of the gantry crane system to the shaped trajectory input with exact natural frequencies using LQR control with ZV and ZVDD shapers are shown in Figure Table summarises the levels of swaying reduction of the system responses at the first three modes in comparison to the LQR control Higher levels of swaying reduction were obtained using LQR control with ZVDD shaper as compared to the case with ZV shaper However, with ZVDD shaper, the system response is slower Hence, it is evidenced that the speed of the system response reduces with the increase in number of impulse sequence The corresponding rise time, setting time and overshoot of the cart trajectory response using LQR control with ZV and ZVDD shapers with exact natural frequencies is depicted in Table It is noted that a slower cart trajectory response with less overshoot, as compared to the LQR control, was achieved To examine the robustness of the shapers, the shapers with 30% error in swaying frequencies were designed and implemented to the gantry crane system The analysis from Figure shows that the swaying angles of the system were considerable reduced as compared to the system with LQR control However, the level of swaying reduction is slightly less than the case with exact natural frequencies Table summarises the levels of swaying reduction with erroneous natural frequencies in comparison to the LQR control The time response specifications of the cart trajectory with error in natural frequencies are summarised in Table It is noted that the response is slightly faster for the shaped input with error in natural frequencies than the case with exact frequencies However, the overshoot of the response is slightly higher than the case with exact frequencies Significant swaying reduction was achieved for the overall response of the system to the shaped input with 30% error in natural frequencies, and hence proved the robustness of the input shapers Optimal Tracking with Sway Suppression Control for a Gantry Crane System Figure 4: Response of the gantry crane with exact natural frequencies (a): Cart trajectory LQR-ZV Shaper LQR-ZVDD Shaper Horizontal position of the cart (m) -1 -2 -3 -4 -5 10 15 20 25 30 Time (s) (b): Sway angle LQR-ZV Shaper LQR-ZVDD Shaper Swing angle of the rope (rad) 1.5 0.5 -0.5 -1 -1.5 -2 10 15 20 25 30 Time (s) (c): Power spectral density 40 LQR LQR-ZV Shaper LQR-ZVDD Shaper Power Spectral Density (dB) 20 -20 -40 -60 -80 -100 0.5 1.5 Frequency (Hz) 2.5 637 M.A Ahmad, R.M.T Raja Ismail, M.S Ramli, N.M Abdul Ghani and M.A Zawawi Figure 5: Response of the gantry crane with erroneous natural frequencies (a): Cart trajectory LQR-ZV Shaper LQR-ZVDD Shaper Horizontal position of the cart (m) -1 -2 -3 -4 -5 10 15 20 25 30 Time (s) (b): Sway angle LQR-ZV Shaper LQR-ZVDD Shaper Swing angle of the rope (rad) 1.5 0.5 -0.5 -1 -1.5 -2 10 15 20 25 30 Time (s) (c): Power spectral density 40 LQR LQR-ZV Shaper LQR-ZVDD Shaper 20 Power Spectral Density dB 638 -20 -40 -60 -80 -100 0.5 1.5 Frequency Hz 2.5 Optimal Tracking with Sway Suppression Control for a Gantry Crane System 639 6.3 Comparative Performance Assessment By comparing the results presented in Table 1, it is noted that the higher performance in the reduction of swaying of the system is achieved using LQR control with ZVDD shaper This is observed and compared to the LQR control with ZV shaper at the first three modes of swaying frequency For comparative assessment, the levels of swaying reduction of the hoisting rope using LQR control with both ZV and ZVDD shapers are shown with the bar graphs in Figure The result shows that, highest level of swaying reduction is achieved in hybrid control schemes using the ZVDD shaper, followed by the ZV shaper for all modes of swaying frequency Therefore, it can be concluded that the LQR control with ZVDD shapers provide better performance in swaying reduction effect as compared to the LQR control with ZV shapers in overall Comparisons of the specifications of the cart trajectory responses of hybrid control schemes using both ZV and ZVDD shapers are summarised in Figure for the rise times and settling times It is noted that the differences in rise times of the cart trajectory response for the LQR control with ZV and ZVDD shapers are negligibly small However, the settling time of the cart trajectory response using the LQR control with ZV shaper is faster than the case using the ZVDD shaper It shows that, by incorporating more number of impulses in hybrid control schemes resulted in a slower response Comparison of the results shown in Table for the shaping techniques with error in natural frequencies reveals that the higher robustness to parameter uncertainty is achieved with the LQR control with ZVDD shaper For both case of the ZV and ZVDD shapers, errors in natural frequencies can successfully be handled This is revealed by comparing the magnitude level of swaying frequency of the system in Figure Comparisons of the cart trajectory response using LQR control with ZV and ZVDD shapers with erroneous natural frequencies are summarised in Figure The results show a similar pattern as the case with exact natural frequencies The system response with ZVDD shaper provides slightly slower responses than the ZV shaper Figure 6: Level of swaying reduction with exact and erroneous natural frequencies using ZV and ZVDD shapers 70 Level of reduction (dB) 60 50 LQR-ZV Shaper 40 LQR-ZVDD Shaper 30 LQR-ZV Shaper (error) 20 LQR-ZVDD Shaper (error) 10 Mode Mode Mode of Vibration Mode 640 M.A Ahmad, R.M.T Raja Ismail, M.S Ramli, N.M Abdul Ghani and M.A Zawawi Figure 7: Rise and settling times of the cart trajectory with exact and erroneous natural frequencies using ZV and ZVDD shapers Time (s) LQR LQR-ZV Shaper LQR-ZVDD Shaper LQR-ZV Shaper (error) LQR-ZVDD Shaper (error) Rise time Table 1: Frequenc y Exact Error Settling time Level of swaying reduction of the hoisting rope and specifications of cart trajectory response for hybrid control schemes Attenuation (dB) of swaying angle Types of shaper Mode Mode Mode ZV ZVDD ZV ZVDD 12.99 38.45 3.94 18.97 20.29 51.20 0.09 25.91 26.26 64.18 12.26 40.17 Specifications of cart trajectory response Rise time Settling Overshoot (s) time (s) (%) 2.287 3.707 0.08 3.449 6.285 0.00 2.004 3.333 0.15 2.786 5.241 0.03 Conclusion The development of hybrid control schemes for input tracking and swaying suppression of a gantry crane system has been presented The hybrid control schemes have been developed based on LQR with feedforward control using the ZV and ZVDD input shapers technique The proposed control schemes have been implemented and investigated within the simulation environment of the overhead gantry crane system The performances of the control schemes have been evaluated in terms of level of input tracking capability, swaying reduction, time response specifications and robustness Acceptable input tracking capability and swaying suppression have been achieved with both control strategies A comparison of the results has demonstrated that the LQR control with input shaping using ZVDD shapers provide higher level of swaying reduction as compared to the cases using ZV shapers In term of speed of the responses, ZVDD shapers results in a slower tracking response with less overshoot It has also demonstrated that input shaping technique is very robust to error in natural frequencies especially with higher number of impulses It is noted that the proposed hybrid controllers are capable of reducing the system swaying while maintaining the input tracking performance of the crane Optimal Tracking with Sway Suppression Control for a Gantry Crane System 641 References [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] Lueg, P 1936 Process of silencing sound oscillations US Patent 043 416 Omar, H.M 2003 Control of gantry and tower cranes Ph.D Thesis, M.S Virginia Tech Manson, G.A 1992 Time-optimal control of and overhead crane model Optimal Control Applications & Methods, Vol 3, No 2, pp 115-120 Auernig, J and Troger, H 1987 Time optimal control of overhead cranes with hoisting of the load Automatica, Vol 23, No 4, pp 437-447 Karnopp, B.H., Fisher, F.E and Yoon, B.O 1992 A strategy for moving mass from one point to another Journal of the Franklin Institute, Vol 329, pp 881-892 Teo, C.L., Ong, C.J and Xu, M 1998 Pulse input sequences for residual vibration reduction Journal of Sound and Vibration, Vol 211, No 2, pp 157-177 Singhose, W.E., Porter, L.J and Seering, W 1997 Input shaped of a planar gantry crane with hoisting Proc of the American Control Conference, pp 97-100 Feddema, J.T 1993 Digital Filter Control of Remotely Operated Flexible Robotic Structures American Control Conference, San Francisco, CA, Vol 3, pp 2710-2715 Noakes, M.W and Jansen, J.F 1992 Generalized Inputs for Damped-Vibration Control of Suspended Payloads Robotics and Autonomous Systems, 10(2): p 199-205 Singer, N., Singhose, W and Kriikku, E 1997 An Input Shaping Controller Enabling Cranes to Move without Sway ANS 7th Topical Meeting on Robotics and Remote Systems, Augusta, GA Kress, R.L., Jansen, J.F and Noakes, M.W 1994 Experimental Implementation of a Robust Damped-Oscillation Control Algorithm on a Full Sized, Two-DOF, AC Induction MotorDriven Crane 5th ISRAM, Maui, HA, pp 585-592 Wahyudi and Jalani, J 2005 Design and implementation of fuzzy logic controller for an intelligent gantry crane system Proceedings of the 2nd International Conference on Mechatronics, pp 345- 351 Ogata, K 1997 Modern Control Engineering Prentice-Hall International, Upper Saddle River, NJ ... input shapers Optimal Tracking with Sway Suppression Control for a Gantry Crane System Figure 4: Response of the gantry crane with exact natural frequencies (a) : Cart trajectory LQR-ZV Shaper LQR-ZVDD... the value of vector K and N which determines the feedback control Optimal Tracking with Sway Suppression Control for a Gantry Crane System 635 law and for elimination of steady state error capability... Implementation and Result In this investigation, hybrid control schemes for tracking capability and sway angle suppression of the gantry crane system are examined Initially, a Linear Quadratic Regulator