The Korean Society of Marine Engineering
[ Original Paper ]
Journal of the Korean Society of Marine Engineering - Vol. 39, No. 5, pp.554-562
ISSN: 2234-7925 (Print) 2234-8352 (Online)
Print publication date Jun 2015
Received 26 Jan 2015 Revised 12 Feb 2015 Accepted 09 Apr 2015
DOI: https://doi.org/10.5916/jkosme.2015.39.5.554

Design and Control of a Firefight Cannon Manipulator Applying Sliding Mode Control

Mai The Vu1 ; Hyeung-Sik Choi ; Hyeon-Seung Kang2 ; Jae-Hyeon Bae3 ; Moon-G Joo4 ; Yeong-do Joo5
1Department of Mechanical Engineering, KMOU, Tel: 051-410-4969 maithevu90@gmail.com
2Department of Mechanical Engineering, KMOU, Tel: 051-410-4969 hsgangmain@naver.com
3Department of Mechanical Engineering, KMOU, Tel: 051-410-4969 doavm789@naver.com
4Department of Information and Communications Eng., Pukyong National University, Tel: 051-629-6238 gabi@pknu.ac.kr
5SMEC CO., Ltd, Tel: 053-524-5211 jyd@esmec.com

Correspondence to: Division of Mechanical and Energy System Engineering, Korea Maritime and Ocean University, E-mail: hchoi@kmou.ac.kr, Tel: 051-410-4969

Copyright © The Korean Society of Marine Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This paper describes an analysis of an architecture and control system of a firefighting cannon manipulator (FCM) composed of two joint axes and one water-shooting actuator. Because the orienting FCM motion is disturbed by the reaction force from water shooting, the water shooting force has been modeled for robust control. The dynamics model of the manipulator has been set up including the external force of water-shooting reaction on the manipulator. A PD Controller and Sliding Mode Controller have been designed and their performance been tested through simulation to track a desired trajectory under the disturbance of a watershooting reaction. The simulation shows that the performance of the Sliding Mode Controller is better than that of the PD controller.

Keywords:

Firefight cannon manipulator, Controller, Kinematics, Dynamics, Simulation

1. Introduction

In recent years, industrial robots have been evaluated for various purposes, including personal care, hazardous material treatment, fire disaster prevention, and others. Unlike industrial robots fixed at a workstation, these robots can conduct surveillance and reconnaissance by moving around the indoor and outdoor environments. This study performed a variety of tasks, including the pick-and-place task [1]. Moreover, a study of the manipulator for locomotion [2], and a study of the light weight manipulator made of carbon fiber reinforced plastic was conducted [3]. One of the best nonlinear robust controllers that can be used in uncertain nonlinear systems is the sliding mode controller (SMC), but pure SMC results in chattering in a noisy environment. This effect can be eliminated by optimizing the sliding surface slope. One study investigated a novel methodology for designing a SMC using a new heuristic search, so called "colonial competitive algorithm " to tune the sliding surface slope and the switching gain of the discontinuous part in the SMC structure [4] and another study illustrated two different control techniques for achieving finite-time convergence and continuous control [5]. These robots comprise the locomotion part and robotic arm used in different operations.

The FCM (Firefight cannon manipulator) is characterized mainly by its light weight and low degree of freedom (namely 2DOF) because it moves with a manipulator attached to the locomotion device. The features of the FMC system are shown in Figure 1.

Figure 1:

Overall view of the developed FCM

The equations of motion of the manipulator have been established and analyzed by including the external force of the water impacts on a system. The performance of the SMC compensating the water impacts have been simulated by Matlab/Simulink and compared with the one of the classical PD controller.

To control the system, the PD controller and SMC have been designed to follow the desired trajectory. Simulations have been conducted to verify the availability of generated trajectories and the performance of the designed controller.


2. Architecture of FCM

2.1 Mechanical structure

The FCM is designed to perform the tasks of fire extinguishing and is composed of an elbow, swivel and nozzle. The system includes two motors (motors 1 and 2), which are used for the pan and tilt motion and the motor control nozzle are used for fog-stream. Each motor is composed of a reduction gear, DC motor and encoder. One important point in the FCM is that the gears applied are worm gears. They are used to prevent backward rotation from the reaction of water shooting.

The structures of the manipulator are composed as follows: The manipulator joints are composed of two degrees of freedom with yawing motion of the first axis and pitching motion of the second axis, as shown in Figure 1 and 2. The first joint actuator is designed to move around in a range between 0° and 360°. The second joint actuator is in a range from -10° to 185°.

Figure 2:

Joints of the FCM

The weights of the two links are 7.55 kg and 4.06 kg each. The total weight of the manipulator system is approximately 11.61 kg. The manipulator has been developed to cope with a maximum load capacity of 123 kg.

Table 1 lists the specification of the measured and calculated FCM and Table 2 and 3 present the parameters of motors and inertial constants.

Parameters of the firefight cannon manipulator

Parameters of the motors

Inertial Constants (kgm2 )

2.2 Control system

A control system has been developed to control the FCM, and the entire control system consists of a main controller, and joint motor controller, as shown in Figure 3. The main controller based on an ARM processor plays the role of editing the control algorithm, sending the order of the operation signals to the motor controller for the robot joint actuators through RS232 communication. The main controller sends continuous motion signals to motor controllers of the joint actuators according to a trajectory planning. For the joint motion controller, a STM32F207 microprocessor is applied. The controller is built in a modular control board capable of controlling two motors and has been designed to generate the Pulse Width Modulation (PWM) signal and the direction signal for the motor and to send them from the main controller to the motor driver. The motor driver has been developed to supply the power, which is proportionally amplified with the duty ratio to the joint actuating motors according to the PWM signal and the direction signal for the joint motor. Figure 4 shows a picture of water-shooting of FCM during according to control input.

A proportional-derivative and feed-forward control algorithm is planned to apply to the joint controller of the developed system.

Figure 3:

The developed control system

Figure 4:

The water-shooting of the FCM


3. Kinematic Analysis of FCM

3.1 Forward Kinematics analysis

An analysis of the kinematics of FCM was performed on the joint actuator links based on a two-axis of FCM. This is because the one-axis of the manipulator is fixed on the ground surface, as the reference axis. The coordinates have been set up using the D-H (Denavit - Hartenberg) convention, as shown in Figure 5, and Tables 4 lists the corresponding joint link parameters of ai, αi, di, and θi Each joint link is expressed as a homogeneous transformation matrix, Ai, by multiplying the two transformation matrices. The matrix of end - effector are obtained:

A2=c12-c1s2s1c2s1-s12-c1s20c2000a1c1+a2c12-l2c12/2a1s1+a2s1c2-l2s1c2/2d1+a2s2-l2s2/21
Figure 5:

D-H coordinates of 2-axis manipulator

Denavit - Hartenberg

3.2 Inverse Kinematics analysis

Kinematic analysis was performed on joint actuator links based on a 2-axis manipulator. Since the 1-axis of manipulator was fixed on the ground surface, it was established as a reference axis. Coordinates were set up using the D-H (Denavit- Hartenberg) convention as shown in Figure 5, and then corresponding joint link parameters of ai, αi, di, and θi (where ai is called the length, αi is called the twist, di is called the offset, and θi is called the angle, respectively) were obtained. The following Equation (1) is derived after matrix transformation through consecutive multiplication. The following Equation (1) is derived after matrix transformation through consecutive multiplication.

cβsθ1+sβcθ1=sθ1+β=m2+n2+a12-a222a1m2+n2sinβ=mm2+n2,cosβ=mm2+n2(1) 

The following Equation (2) is obtained from Equation (1).

θ1=sin-1m2+n2+a12-a222a1m2+n2-β(2) 

Similarly, the following equations are obtained.

θ2=sin-1n-a1sinθ12a2-θ1(3) 

4. Dynamics Analysis of FCM

The Euler-Lagrange method is used repetitively to perform the dynamics analysis of the two-link of FCM. The dynamic equation expressing the linear and rotational motion for each link of the manipulator is induced one after another.

The dynamic equation begins with the initial condition, ω0 = a0 = ac,0 = ae,0 = [0 0 0]T, and the related equations are then obtained by consecutively increasing the i value from 1 to n.

4.1 Modeling of water-shooting pump

In this study, a disturbance of the system is expressed as an external force. For a control volume (CV), the second link of a system where the water move out is enclosed. The only mass flow across the CV is water moving to the right. Figure 6 shows the CV of second link.

Figure 6:

Control volume of the second link

The force balance is:

Fx=-F=ṁoutvout(4) 

where, ṁout is the mass flow, vout is maximum velocity, respectively. And the mass flow ṁout is given by Equation (5):

ṁout=ρvoutAout=ρvπD24(5) 

With Aout, ρ and D are section area, water density and diameter of end-effector, respectively.

The final equation of water reaction force of FCM system is:

F=-ρvπD24vout2(6) 

For the FCM system to work, it needs a water pump that can deliver a large volume of water at high pressure through a hoselike tube. This water pump will change the velocity of the water when it moves out.

In this study, maximum velocity vout = 14m / s of water has been chosen to demonstrate the robustness of SMC controller. The diameter of the shooting nozzle is D = 7cm and the water density in this case is ρ = 1000(kg/m3) . Therefore, the maximum disturbance force is 754.3(N).

4.2 Dynamics of 2 DOF FCM

In the following sections, the 2 DOF FMC dynamic equations of motion can be concisely expressed concisely as:

Mqq̈+Cq,q̇q̇+Gq=τ+JETF,(7) 

where, q: 2x1 is position vector, τ : 2x1 is vector of torques, M(q) : 2x2 is inertia matrix of the manipulator, Cq,q̇ ): 2x2 is vector of Centrifugal and Coriolis terms, G(q): 2x1 is vector of gravity terms, JET is Jacobian matrix and F is external forces, respectively.

Matrix M and C are symmetric 2x2 matrices:

Mq=m11m12m21m22Cq,q̇=c11c12c21c22

The inertia matrix is given by Equation (8),:

Mq=i=12JTiTmiJTi+JRiTIiJRi(8) 

where mi,i=1,2 is the mass of link each, JTi and , JRi represent the Jacobian of linear velocity an angular velocity of link each, respectively.

While the Centrifugal and Coriolis matrix Cq,q̇) is given by Equation (9):

Cq,q̇=MqqIn×q̇-12Mqqq̇×InT(9) 

Matrix G is given by Equation (10):

giq=Πqi(10) 

where, Π is the potential energy of the manipulator, which is just the sum of those of the two links. For each link, the potential energy is just its mass multiplied by the gravitational acceleration and the height of its center of mass.

By starting with the end conditions of fn+1 = 0 , τn+1 = 0 . The torques τi of each joint are calculated consecutively. The detailed expression of Equation (8), (9), (10), and the torques τi are described in Appendix.


5. Controller Design and Computer Simulation of FCM

In this paper, the input signals used for the simulation are θ1 = 0.1sin(πt) (rad) and θ2 = 0.1cos(πt) (rad). Both the PD controller and SMC have been designed and their performance been tested through simulations to track a desired trajectory under a disturbance due to a water reaction and without a disturbance. The simulation shows that the performance of the SMC is better than that of the PD controller.

5.1 Design PD Controllers and Computer Simulation

The SMC and PD controller have been designed to track the desired trajectories of a vehicle and joints. A PD control algorithm is applied to control the position each link of system. The PD control laws for each link are described in Equation (11).

τ=Kpe+Kdė(11) 

where ei = xdi - xi, ėi = ẋdi - ẋi, Kpi and Kdi denote a proportional and derivative gain, respectively. The gains used for the PD controller are Kp = 90000, Kd = 600, respectively. Figure 7 represents a block diagram of the PD controller with a disturbance.

Figure 7:

Block diagram of PD controller with a disturbance

In the first simulation, the PD controller is applied to the system. The trajectory of the end-effector is assigned to move from the initial position to the target position with a constant velocity for 10 seconds. Figure 8, 9, 10 and 11 present the simulation results. Figure 8 and 9 are the simulation results of manipulator without a disturbance. Figure 10 and 11 show the simulation results under a disturbance. The simulation shows that the measured trajectory is tracking the desired trajectory with good performance. Figure 9 and 11 show that the errors of two joints are quite small (i.e., without disturbance, the error of the first a second joint is 0° and approximately 7 10-3°. On the other hand, under the disturbance of a water reaction, the error of first joint varies from -0.01° to 0.01° and second joint is varied from 0.005° to 0.014°).

Figure 8:

Joint angles (PD control without a disturbance)

Figure 9:

Joint torques and position errors of the joints (PD control without a disturbance)

Under a disturbance, the input torques are larger than that without a disturbance (without disturbance, the input torques of the first and second link are approximately 11Nm and 0.3Nm, respectively, and under disturbance, the input torques of the first and second link are approximately 20Nm and 18Nm, respectively), as shown in Figure 9 and 11.

Figure 10:

Joint angles (PD control with a disturbance)

Figure 11:

Joint torques and position errors of joints (PD control with a disturbance)

5.2 Design of the SMC Controllers and Computer Simulation

The sliding mode method is based on the idea of keeping the scalar quantity, s, which is a weighted sum of the position error (q - qd), the velocity error(q̇ - q̇d, and (not required) the acceleration error(q̇ - q̇d), at zero as reported elsewhere [6]. Here, the expression of s is chosen as expressed in Equation (12):

s=q̇+q̇d+λq-qd(12) 

where λ > 0 is the weight parameter.

Therefore, the task of the controller is to take s to zero. When s approaches zero, the position error (and velocity error, also) also approaches zero too. Therefore, and thus, trajectory tracking is performed. Once s is zero, to keep it at this value, the derivative of s is expected to be zero. From Equation (12), the expression of s Ṡ can be deduced easily as follows:

ṡ=q̈+q̈d+λq̇-q̇d(13) 

Substituting the expression of q̈ deduced from Equation (7) into Equation (13):

ṡ=1Mqτ+JETF-Cq,q̇q̇-Gq-qd̈+λq̇-q̇d(14) 

So, if derivative of the scalar quantity Ṡ = 0, we have:

τ=Mqqd̈-Mqλq̇-q̇d+Gq+Cq,q̇q̇-JETF(15) 

In addition, the actual control law, which can be robust to uncertainties, is the function chosen as follows:

τ=τ-Ksatsϕ(16) 

where, sat(.) is the saturating function, ㅣsat(.)ㅣ≤ 1, φ > 0 is the boundary layer thickness, K is the design parameter chosen so that ṡ ≤ -η ㅣsㅣ < 0, with η is a strictly positive constant.

The parameters, φ, λ and K used for the SMC controller are φ=1, λ = [800; 800]T and K = [500; 400]T, respectively. To verify the good performance of the designed SMC, the trajectories were have been created for links and the simulation of trajectory tracking was has been carried out. Figure 12 presents the block diagram of SMC with a disturbance.

Figure 12:

A block diagram of the SMC with a disturbance

In the second simulation, the SMC was applied to the system. The trajectory of the end-effector was assigned to move from the initial position to the target position with a constant velocity for 10 seconds. Figure 13 - 16 present the simulation results. Figure 13 and 14 show simulation results of the manipulator without a disturbance. Figure 15 and 16 present the simulation results under disturbance. The simulation result shows that the measured trajectory is tracking the desired trajectory with good performance. Each joint angle is still controlled well with a disturbance. The simulation results show the validity of the SMC algorithms. Figure 14 and 16 show that the errors of the two joints are very small (i.e., without a disturbance, the error of the first and second joint is 0° and from -0.25×10-5 degrees to 0.25×10-5 degree. On the other hand, under a disturbance of the water reaction, the error of first joint varies from -2.2×10-3 degrees to -2.2×10-3 degrees and that of the second joint varies from -1.87×10-3 degree to -1.87×10-3 degree).

Figure 13:

Joint angles (SMC without a disturbance)

Figure 14:

Joint torques and position errors of joints (SMC without a disturbance)

Figure 15:

Joint angles (SMC with a disturbance)

Under a disturbance, the input torques are larger than that without a disturbance case (without a disturbance, the input torques of the first and second link are approximately 11Nm and 0.3Nm, respectively, and under a disturbance, the input torques of the first and second link are approximately 20Nm and 18Nm, respectively), as shown in Figure 14 and 16.

Figure 16:

Joint torques and position errors of joints (SMC with a disturbance)

5.3 Analyzed quality of the controller based on combining the data between the PD Controllers and SMC

The performance of the PD Controller and SMC have been compared. Figure 17 and 18 show the simulation results of the PD Controller and SMC Controller for both cases without a disturbance and under the presence of a disturbance. According to the results, the SMC shows good trajectory tracking along the target joint angle.

Figure 17:

Position errors of the joints (PD control and SMC without a disturbance)

Without a disturbance, the errors of the joints are very small, as shown in Figure 17. Under a disturbance, the simulation results show that the position errors are larger than the case without, as shown in Figure 18. Based on the results of comparing the performance of the SMC and PD controllers, the superior performance of the SMC in reducing the position errors is confirmed.

Figure 18:

Position errors of the joints (PD control and SMC with a disturbance)


6. Conclusions

In this paper, a firefighting cannon manipulator (FCM) has been developed for extinguishing fires. The architecture of the developed mechanical system composed of two joint axes and one water-shooting actuator and control system for exact water shooting for fire has been described.

For the control joint of a FCM, the dynamics model of FCM composed of two joint axes and one water-shooting actuator has been set up including the external force of the water-shooting reaction on the FCM. Because the FCM motion is affected by the reaction force from water shooting, the water shooting force has been modeled.

To control the FCM accurately along the desired trajectory under an external force of a water-shooting reaction, a robust controller has been designed using the SMC theory. Through a number of computer simulations applying the SMC and PD controller, the results of position errors and applied input torques have been presented. When both controllers are compared, the performance of the sliding mode controller is better in reducing the position errors under disturbance of water-shooting reaction without the extra input torques of the joints than that of the PD controller.

For further research, experiments applying the proposed controller to actual FCM will be carried out.

Acknowledgments

This study was supported by National Emergency Management Agency (2012-NEMA10-004-01010002-2012). In addition, it was supported by the development of a multi-legged underwater walking robot project under the Ministry of land, infrastructure, and transport.

References

  • R. R. Murphy, “Activities of the rescue robots at the World Trade Center from 11-21 September 2001”, IEEE Robotics and Automation Magazine, 11(2), p50-61, (2004).
  • T. Frost, C. Norman, S. Pratt, and B. Yamauchi, “Derived performance metrics and measurements compared to field experience for the PackBot”, Proceedings of the 2002 PerMIS Workshop, p201-208, (2002).
  • J. R. White, T. Sunagawa, and T. Akajima, “Hazardous duty robots experiences and needs”, Proceedings of the IEEE/RSJ International Workshop on Intelligent Robots and Systems, p262-267, (1989). [https://doi.org/10.1109/iros.1989.637917]
  • A. Jalali, F. Piltan, M. Keshtgar, and M. Jalali, “Colonial competitive optimization sliding mode controller with application to robot manipulator”, International Journal of Modern Education and Computer Science, 5(7), p50-56, (2013). [https://doi.org/10.5815/ijisa.2013.07.07]
  • M. Bhavea, S. Janardhananb, and L. Dewana, “A finitetime convergent sliding mode control for rigid underactuated robotic manipulator”, Proceedings of the Systems Science & Control Engineering: An Open Access Journal, 2(1), (2014).
  • J. J. Slotine, and W. Li, Applied Nonlinear Control, Prentice Hall, (1991).

Appendix

Inertia matrix M

m11= Ix22 + Iy1 + m2a1+a2cq2-l22cq22+ a12m1-Ix22c2q2;m12=m21= Ix22s2q2;m22=m2a22- m2a2L2+ m24L22+ Ix2 cq22+ Iz2;

Matrix of Centrifugal and Coriolis terms C

c11= q̇2Ix22s2q2- m2a2sq2-L22sq2a1+a2cq2-L22cq2;c12= q̇1Ix22s2q2- m2a2sq2-L22sq2a1+a2cq2-L22cq2+ Ix2q̇2c2q2;c21= -q̇1Ix22s2q2- m2a2sq2-L22sq2a1+a2cq2-L22cq2;c22= -Ix2q̇2cq2sq2;

Matrix of gravity terms G

g1=0;g2=L2gm2cq2;

Matrix of torques τ

τ1= -Ix2c2q2q̇22- 2q̇112Ix2s2q2+m212L2sq2-a2sq2a1-12L2cq2+a2cq2q̇2++q̈1Ix22+Iy1+m2a1-12L2cq2+a2cq22+a12m1-12Ix2c2q2+12q̈2Ix2s2q2+Fsq1a1-12L2cq2+a2cq2;τ2= 12Ix2s2q2+m212L2sq2-a2sq2a1-12L2cq2+a2cq2q̇12+12Ix2s2q2q̇22+q̈214m2L22-m2L2a2+m2a22+Ix2cq22+Iz2+12q̈1Ix2s2q2-12F cq1sq2L2-2a2+a2gm2cq2;

where c denotes cos while s denotes sin.

Figure 1:

Figure 1:
Overall view of the developed FCM

Figure 2:

Figure 2:
Joints of the FCM

Figure 3:

Figure 3:
The developed control system

Figure 4:

Figure 4:
The water-shooting of the FCM

Figure 5:

Figure 5:
D-H coordinates of 2-axis manipulator

Figure 6:

Figure 6:
Control volume of the second link

Figure 7:

Figure 7:
Block diagram of PD controller with a disturbance

Figure 8:

Figure 8:
Joint angles (PD control without a disturbance)

Figure 9:

Figure 9:
Joint torques and position errors of the joints (PD control without a disturbance)

Figure 10:

Figure 10:
Joint angles (PD control with a disturbance)

Figure 11:

Figure 11:
Joint torques and position errors of joints (PD control with a disturbance)

Figure 12:

Figure 12:
A block diagram of the SMC with a disturbance

Figure 13:

Figure 13:
Joint angles (SMC without a disturbance)

Figure 14:

Figure 14:
Joint torques and position errors of joints (SMC without a disturbance)

Figure 15:

Figure 15:
Joint angles (SMC with a disturbance)

Figure 16:

Figure 16:
Joint torques and position errors of joints (SMC with a disturbance)

Figure 17:

Figure 17:
Position errors of the joints (PD control and SMC without a disturbance)

Figure 18:

Figure 18:
Position errors of the joints (PD control and SMC with a disturbance)

Table 1

Parameters of the firefight cannon manipulator

Link 1 2
Length of Link, Li (mm) 282.81 275.56
Length of mass center (mm) 141.4 137.78
Mass, Mi (kg) 7.55 4.06
Motion R (Yaw) R (Pitch)
Limit of angles (degrees) 0°~360° -10°~185°

Table 2

Parameters of the motors

No. Axis Motor Gear ratio Worm wheel
1 1st axis 272763 EC-max30, 60 (W), 24 (V) 23:1 50:1
1 2nd axis 272768 EC-max30, 40 (W), 24 (V) 23:1 50:1

Table 3

Inertial Constants (kgm2 )

Link Mass (Mi) (kg) Inertial moment (kgm2)
Ixx Iyy Izz Ixy = Ixy = Ixy
1 7.55 0.050 0.045 0.045 0
2 4.06 0.026 0.015 0.015 0

Table 4

Denavit - Hartenberg

Link ai (mm) αi di (mm) θi
1 100 90 282.81 θ1
2 275.56 0 0 θ2