Motion Control System Design for a Flying-Type Fireﬁghting System with Water Jet Actuators

: This paper presents the design and modeling of a ﬂying-type ﬁre extinguishing system. Fire accidents present very hazardous environments, and ﬁreﬁghters are in danger of losing their lives while putting out the ﬁre. Strict safety measures should be considered to guarantee safe working conditions for ﬁreﬁghters, which is not the case every time, as fatalities and casualties are still being recorded. For this reason, a novel ﬁre extinguishing system is proposed to provide more safe ﬁreﬁghting and survivor searches. The system studied in this paper is a pilot model that consists of a water jet-based actuation system to control the ﬂying motion of the robot. The dynamic model of this ﬂying robot is derived using the actuation forces, water jet system characteristics, and related information. The mathematical system model is detailed, a sliding-mode control system and a proportional-integral-derivative controller are designed, and comparative simulation tests are carried out.


Introduction
In this research, a new fire extinguishing system that can suppress flames quickly and efficiently is studied.Fire accidents are unfortunate disasters that damage properties and endanger the lives of people, including first responders.In recent years, the risk of public safety and the difficulty of responding to such incidents is gradually increasing due to complex and large-scale buildings, such as airports, stadiums, exhibition halls, warehouses, and other vital construction areas.The risk of firefighters and victims getting secluded and losing their lives due to post-blast explosions and structures collapsing is very high, which makes it difficult to perform active firefighting tasks.
In the case of a fire onboard a ship, it is difficult to suppress the flames due to the unstable environment surrounding the ship (mainly waves), and it is best to simply direct the stream of water from a distance.And unlike land fires, the fire progresses rapidly depending on the materials constituting the ship hull.For instance, in the case of steel ships, it is difficult to approach the fire to put it out, due to high temperatures.Initially, if the fire suppression is deemed difficult, the firefighting personnel must go to the fire scene to search for and rescue survivors, a task that is carried out at risk.However, it is extremely difficult to quickly respond because of all the risk factors, which makes it difficult to reduce the degree of the damage caused by the response delay.
Various methods and technologies have been devised for rapid-fire suppression while minimizing the imminent damages bound to occur, but they are not to the level that they can respond adequately, quickly, and safely.Therefore, remote-controlled and autonomous fire extinguishing techniques became an active research area that aims to improve the efficacy of firefighting, fire suppression, and to reduce firefighter injuries and deaths.These systems may also increase the efficiency of firefighting tasks through advanced vision applications and additional sensors to support looking for survivors in low-visibility environments.
Though for small-scale fires many of the current robots developed for such outcomes are remotely controlled.For example, P. Liljeback et al. [1] suggested a water-powered actuation system enabling a fire hose to move like a snake and perform the fire intervention on its own.The shipboard autonomous firefighting robot (SAFFiR) is one of the earlier humanoid robots that helped researchers as a support tool in inspecting and suppressing fires aboard naval vessels [2].However, for large-scale fires or when the flames have already engulfed the area, deploying such technologies is impossible.Subsequently, other techniques have been developed to fight such fire and suppress it.A typical example is a ground-type firefighting robot equipped with a caterpillar track or multiple wheels to perform fire suppression tasks even in areas where firefighters cannot enter, such as obstacles and rough roads [3].Some commercial products have been launched, such as the remote-controlled Colossus robot in France [4] and Thermite RS1 and Thermite RS3 in the USA [5].However, due to the lack of mobility and clear vision in underground and complex structures, the use of such technologies is limited to specific fire sites.
Recently, research studies on indoor and outdoor fire suppression systems using unmanned aerial vehicles (UAV) [6][7][8] and drones [9,10] are being actively conducted for this purpose.However, due to the characteristics of UAVs and drones, it is difficult and even aggravating to get them close to the fire and ignition locations, due to the airflow resulting from the rotor blades that may cause the flames to spread to neighboring areas.In other words, even if a drone is used in the process of suppressing the flames, water and fire extinguishing materials must be sprayed and projected from a distance, but also closer to the fire location.For this reason, firefighter fatalities and casualties may be inevitable.Hence, the pursuit of safety is important not only for the firefighters but also for the survivors.
On the other hand, water-powered flying systems, such as personal jet packs [11][12][13], flying boards [14,15], and aerial systems [16,17], have shown their capability and flight maneuverability.Unfortunately, these systems are either manual or semi-automated.For autonomous operations, the available control approaches rely on mechanisms that can regulate both water flow and posture of the actuator nozzles.Furthermore, the motion control methods are fairly simple, such as proportional derivative (PD) control [16], proportional integral derivative (PID) control [13], and proportional control with speed feedbacks [18].As an extension of this feature, for the autonomous firefighting tasks, it is reasonable to use the available water source to actuate the motion of the flying system.
Accordingly, this study proposes an active water-powered fire suppression system that can minimize the loss of life and perform quick and efficient fire extinguishing operations.The proposed system is a device that can fly directly into the fire, remotely accessing the ignition location and extinguishing the fire with direct and precise fire extinguishing water-spraying.A study of a close concept is the Dragon Firefighter developed by the research laboratory of Tohoku University, under the supervision of Professor Satoshi Tadokoro [18][19][20].It is an aerial continuum robot with water jet propulsion and a complex actuation system.The Dragon Firefighter consists of multiple rotary nozzles distributed on two modules.The two modules and the water source are connected in series by a flexible hose with an oscillation suppression mechanism.The water jets from the nozzles propel the modules and at the same time, suppress the fire.Moreover, the objectives of controlling the water jet system are limited to stabilize [18] and rotate [20] the nozzle modules.That is, the developed nozzle units shoot constant water corresponding to the quantity supplied by the pump and rotate the nozzle angle to stabilize the head posture.Its sole objective is to extinguish a fire located on an overhand obstacle, since it focuses primarily on achieving stable flight and not on trajectory tracking, to extinguish fires located in deep locations.This is due to the actuation method used, which is not enough to control the 6-DoF motion.To overcome this limitation, an advanced control method of the actuation force or supplying additional flight force would be required.
Accordingly, this study proposes an active fire suppression system that can minimize the loss of life and perform efficient fire extinguishing operations with a method providing enough actuation force for both translational and rotational motions.The system proposed in this paper uses one head unit consisting of an actuation system with fixed water nozzles and an independent water spray unit used for fire extinguishing.The system model is formulated, and a cascade robust control scheme is designed to effectively perform 6-DoF flight maneuvers.As a result, the proposed firefighting system will be able to be in close proximity to the fire location.Therefore, the proposed device can suppress the fire with low pressurized water spray compared to the typical water extinguishing techniques that use high-pressure water that can be destructive and harmful upon human contact.In addition, it is important to establish a system that can reduce the risk of loss of life by remotely approaching the source of fire without the direct intervention of firefighters.
Therefore, the contributions of this work, in comparison with relevant studies, are as follows: • The basic model of a fire extinguishing water-spraying system that can fly into the flames and suppress the fire is proposed.The existing models of water-powered flying systems, such as in [13,16,18,21], are fairly simple.

•
The mathematical model derived in this paper evaluates and analyzes the dynamic characteristics of the mechanism in all 6-DoF motions.

•
A sliding-mode control (SMC) system is designed for motion control with the designed actuation system.Although the available automated approaches have complex actuation systems, the control methods are rather simple.Thus, the SMC is chosen since SMC techniques are well known for their robustness facing disturbances and uncertainties.

•
Simulation tests are carried out with the SMC and a PID controller.
This paper is organized as follows.Section 2 presents the concept of the novel flyingtype firefighting system, and the mathematical model is derived.A sliding mode control system is designed in Section 3. Comparative simulation tests are carried out between the designed SMC and a PID controller, and the results are listed in Section 4. Finally, conclusions are drawn in Section 5.

Concept of the Proposed System
The fire extinguishing system proposed in this study can move remotely to the fire location.This location is mainly unknown, uncertain, and most importantly, dangerous.Therefore, instead of firefighters, the system proposed is configured so that the hose and spraying nozzle can be remotely placed within the blaze to monitor the flames in real-time, while simultaneously suppressing and extinguishing the fire.
Figure 1 depicts the conceptual diagram of the flying-type fire extinguishing system under development, and the devices constituting the system are as follows:

•
Water supply hose used for motion control and fire extinguishing operation.

•
Fire extinguishing water and motion control valve.

•
Water spray nozzle and a spray direction control device.

•
Water jet-type thrust and steering device.

•
Heat-resisting and protection device.
Using the water jet actuator consist of supplying four nozzles located on both sides of the robot body (two nozzles on the left and two on the right), and the resulting repulsive force controls the flight motion of the flying fire extinguishing system.However, when the flying altitude is high and the water pressure is insufficient to keep the system afloat, or when it is necessary to increase the flying speed in a relatively open space, the additional propeller can be used to secure good and accurate performance.However, at the ignition location, where the fire extinguishing operation should take place, the resulting wind from the propeller is not desirable.Moreover, considering the high temperature of the flames on site, it is necessary to reduce as much as possible the effect of the propulsion of the propeller actuator, as it can further worsen the situation.Therefore, the method of using a water-based driving system is preferable, as it serves as an actuation system and an additional fire extinguishing option.Furthermore, using electric devices in such environments may lead to malfunctions due to high temperatures; thus, the use of a wirebased operating mechanism.In short, if the system in consideration is configured under these circumstances, the efficiency of the fire extinguishing operation can be maximized.Using the water jet actuator consist of supplying four nozzles located on both sides of the robot body (two nozzles on the left and two on the right), and the resulting repulsive force controls the flight motion of the flying fire extinguishing system.However, when the flying altitude is high and the water pressure is insufficient to keep the system afloat, or when it is necessary to increase the flying speed in a relatively open space, the additional propeller can be used to secure good and accurate performance.However, at the ignition location, where the fire extinguishing operation should take place, the resulting wind from the propeller is not desirable.Moreover, considering the high temperature of the flames on site, it is necessary to reduce as much as possible the effect of the propulsion of the propeller actuator, as it can further worsen the situation.Therefore, the method of using a water-based driving system is preferable, as it serves as an actuation system and an additional fire extinguishing option.Furthermore, using electric devices in such environments may lead to malfunctions due to high temperatures; thus, the use of a wire-based operating mechanism.In short, if the system in consideration is configured under these circumstances, the efficiency of the fire extinguishing operation can be maximized.
The proposed configuration of this flying fire extinguishing system is different from previous works with a similar concept [18][19][20].This configuration focuses mainly on a water-based driving system with a control valve that can supply enough actuation force for both translational and rotational movements.Moreover, the additional propulsion actuator with a wire-based guiding mechanism is further introduced, the detailed structure of which will be introduced in a future study.

System Modeling
The reference coordinate of the 6-DoF firefighting system is illustrated in o x y z ) in the earth fixed frame (e-frame, oxyz ) is defined by its position ( x , y , z ) and the three Euler angles (roll, pitch, yaw) ( φ ,θ ,ϕ ).The rotation matrix between the two frames is given as [22]: cos cos cos sin sin sin cos cos sin cos sin sin sin cos sin sin sin cos cos sin sin cos cos sin sin cos sin cos cos The transformation matrix between the body-fixed angular velocity The proposed configuration of this flying fire extinguishing system is different from previous works with a similar concept [18][19][20].This configuration focuses mainly on a water-based driving system with a control valve that can supply enough actuation force for both translational and rotational movements.Moreover, the additional propulsion actuator with a wire-based guiding mechanism is further introduced, the detailed structure of which will be introduced in a future study.

System Modeling
The reference coordinate of the 6-DoF firefighting system is illustrated in Figure 2. The location of the body-fixed frame (b-frame, o b x b y b z b ) in the earth fixed frame (e-frame, oxyz) is defined by its position (x,y,z) and the three Euler angles (roll, pitch, yaw) (φ,θ,ϕ).The rotation matrix between the two frames is given as [22]: The transformation matrix between the body-fixed angular velocityω b = ω xb ω yb ω zb T , and the Euler rate T is given as [23]: The schematic drawing of the 6-DoF flying type firefighting system with four water jet actuators is illustrated in Figure 3.The nozzles of the water jet actuators are connected to the sides of the head part at an angle of δ from the z b -axis and F i (i = 1, . . .4) represent the force of each actuator.Equations of translational motions and rotational motions of the head part are given as follows: where M = M h + kρAz is the mass being lifted and it consists of the mass of the head and the mass of water in the water hose connected to the head part.ρ is the density of water, A is the cross-sectional area of the water hose, and k is a correction factor.Then, Equation (3) can be rewritten as the following: Supposing that the moment of inertia matrix J is diagonal, Equation (4) becomes: In Equation ( 5), the force elements F (•)b are defined as: (7) and the torque elements in Equation ( 6) are defined as: where F d = F dx F dy F dz T and τ d = τ dx τ dy τ dz T are the disturbances induced by the wind, ground effect, and the constraints of the water-supplying hose.Moreover, since the proposed firefighting system is actuated by a water-powered propulsion system, the flow of the fluid should be considered.The derivation of the linear momentum equation of the water-powered system gives the values of the actuator forces as follows [24]: In Equation ( 9), the left-hand side is the summation of the actuator forces acting on the system.The first term on the right is the rate of change of the momentum inside the controlled volume.The last term is the net rate of the flow of the linear momentum through the controlled surface.For this system, the controlled volume is defined as the space inside the water system initialing at the water pump inlet and ending at the actuating nozzle outlet.For the simplicity of modeling, the following assumptions are made: Therefore, the actuator forces and torques are respectively derived as follows: where .m i is the mass flow rate of the water of the ith nozzle.a is the cross-sectional area of the nozzle outlet, and ρ is the density of water.Moreover, let us adopt the following inputs: U 1 corresponds to the total thrust, U 2 affects the roll orientation and the translation along the y b -axis, U 3 controls the pitch and the translation along the x b -axis, and U 4 is proportional to the yaw moment.Due to the characteristics of the water-powered system the mass flow rate range is 0 ≤ .m i ≤ . m max and one can see that the inputs are saturated: The resulting system model is based on the above-derived equations with the state vector defined as X = x y z φ θ ϕ T and .
and is given by: .. where: Remark 1.The system consists of two fully actuated channels-vertical and yaw motions-and four underactuated channels.Longitudinal and pitch motions are considered one underactuated system, while lateral and roll motions constitute the other underactuated system with input coupling.
Additionally, let the following assumption be true, which will be useful in the design process of the controller.

Control System Design
In this section, a control system is presented to address the motion control problem of the flying-type firefighting systems in the presence of parametric uncertainties and other disturbances, such as the wind, the ground effect, and constraints of the water-supplying hose.Thus, the SMC is chosen since this control theory is well known for its robustness in facing disturbances and uncertainties, as well as not requiring explicit knowledge of the controlled system.From Equations ( 14)-( 16) the system model can be rewritten as follows:

Altitude and Position Control
The objective of the altitude and position control is to bring the system to a desired position in the reference coordinate.Let us assume that the rotation matrix is slowly varying and define the position error as the difference between the actual position vector and the desired position vector denoted with a subscript d: The desired position vector is assumed to be twice differentiable.Moreover, consider the control error defined as the following: where b 2 = − sin δ/M h ρa and b 3 = cos δ/M h ρa.If the system is properly stabilized, then t 0 ω xb (ι)dι ≈ φ.Thus, the time-derivative of Equation ( 20) is given by: .e X =   .
The sliding manifold for the altitude and position control system is chosen as follows: where the n X1 and n X2 are the diagonal gain matrices of positive constants.The timederivative of the sliding manifold yields: . s X = .
Then, the control law based on the sliding-mode technique is given by: where K X > 0 is a diagonal matrix and f (s X ) = s X x s X y sgn(s X z ) T , with sgn(•) is the signum function.Based on Equation ( 24), the desired value of θ, φ and U 1 can be calculated from the following control law: Remark 3. θ d and φ d are continuous.
For stability analysis, let us consider the following Lyapunov function candidate: Then, the time-derivative of this function yields: .
where D x , D y , D z are defined as follows As a result of Assumption 4, D i , with i = x, y, z, is also bounded.The negative definiteness of .V 1 is preserved if the following conditions are jointly satisfied: This implies that the altitude and position control system is stable.Additionally, the sliding variables s Xx and s Xy are bounded by max(D x )/K Xx and max D y /K Xy , respec- tively, and the higher the controller gains are, the smaller the size of the bound.

Attitude Control
The attitude controller is responsible for keeping the 3D orientation of the system at the desired value.To design the control system, let the attitude error be as follows where Φ d is the desired attitude vector of the system and Φ is the actual attitude vector, both are represented by Euler angles.In a similar manner to altitude and position control, the desired attitude vector is also assumed to be twice differentiable.As mentioned in Section 3.1, the desired value φ d and θ d are calculated from the resulting control law of altitude and position control system, as shown in Equation (25).Then, the time-derivative of the attitude error can be calculated as follows: . e Φ = .
The sliding manifold for the attitude control system with the diagonal gain matrix n Φ as follows, The time-derivative of Equation ( 31) is given by: . s Φ = ..
Accordingly, we propose the sliding mode control law as the following: where sgn(s Φ ) = sgn(s φ ) sgn(s θ ) sgn(s ϕ ) T and K Φ is a positive definite diagonal gain matrix.Let us introduce the second Lyapunov function for stability analysis of the attitude control system as follows: Then, along with Equations ( 33) and (34), the time-derivative of the second Lyapunov function is given as: .
Therefore, by choosing K Φ such that each of its diagonal elements is greater than the upper bound of the corresponding element of ..
the stability of the attitude control system is preserved.From Assumption 4, the elements of ..
Figure 4 shows the block diagram of the complete motion control system for the flying-type firefighting system.
Therefore, by choosing Φ K such that each of its diagonal elements is greater than the upper bound of the corresponding element of τ τ τ the stability of the attitude control system is preserved.From Assumption 4, the elements Figure 4 shows the block diagram of the complete motion control system for the flying-type firefighting system.

Comparative Control Methods
For a comparative study of the motion control system, PID controllers are designed based on the simplified model of the system.By assuming that the system is properly stabilized ( 0 φ ≈ and 0 θ ≈ ) and the yaw angle is a slow-varying value ( 0 ϕ ≈


) and neglecting the uncertain terms from Equation ( 19), the simplified model is written as follows: Let us define the following state errors:

Comparative Control Methods
For a comparative study of the motion control system, PID controllers are designed based on the simplified model of the system.By assuming that the system is properly stabilized (φ ≈ 0 and θ ≈ 0) and the yaw angle is a slow-varying value ( .ϕ ≈ 0) and neglecting the uncertain terms from Equation ( 19), the simplified model is written as follows: ..
Let us define the following state errors: Then, the control laws based on the PID control technique is proposed as follows: where the computed control signals are u j = θ d , φ d , U 1 , U 2 , U 3 , U 4 and the corresponding state errors are e j = e bx , e by , e bz , e φ , e θ , e ϕ .Each of the PID controller gains is tuned by the SIMC rule for the optimal PID control of double integrating processes [25].
where the k j is the process gain which can be k j = g, −g, b, b 4 , b 5 , b 6 in respect to each u j .ι 0j is the time delays of the respective process.The tuning parameters ι cj are used as a trade-off between the performance and robustness of the control system.

Simulation
To illustrate the usefulness of the flying-type firefighting system model and its motion control system, simulation studies are conducted.Table 1 lists the parameters of the controlled system used in the simulation, based on the pilot model that has been under development in the Marine Cybernetics Laboratory, and is located at the Department of Mechanical System Engineering, Pukyong National University.By setting the reference trajectory as a three-dimensional spline Equation (42) [26], we can show that the control system proposed is capable of tracking complex desired trajectories from the flying-type firefighting system, even in the presence of disturbances.
x d = 10 sin(0.1t) The simulation is carried out through two scenarios: tracking the trajectory in (1) the absence of the disturbances and (2) in their presence.The disturbances used in the second simulation scenario are given as [26]: 0.1 sin(0.1t)0.1 sin(0.1t+ 0.3) 0.1 sin(0.1t− 0.1) 0.15 sin(0.3t+ 0.2) 0.15 sin(0.3t− 0.25) 0.15 sin(0.3t) With this simulation configuration, the controller gains are sufficiently tuned for the best performance possible, and their values are listed in Table 2.
The simulation results of the first scenario are presented in Figures 5-8. Figure 5 represents the three-dimensional position of the head part, Figures 6 and 7 represent the position and attitude control response.In Figures 5-7, it is clear that both controllers are able to track the reference trajectory well with a good transient response.The errors of the z position and yaw angle, in particular, are noticeably small.However, the amplitudes of other state errors using the PID control are significantly larger than with the proposed control system.Figure 8 shows the mass flow rate of each nozzle.The mass flow rate necessary to raise the altitude to 15 m and track the complex trajectory is within 2 kg/s.This is feasible, and these results can be utilized to develop the testbench and pilot model of the proposed system.12 is at a similar level to the one depicted in Figure 8 even though the effect of disturbances is considered.
As for the second simulation scenario, the disturbances in Equation ( 43) are ap plied.The results are shown in Figures 9-12.The position and attitude time response i depicted in Figures 9-11.Similar to the simulation results without disturbances, both controllers can lead the head part along the reference trajectory.However, residual vi brations appear in the state errors.The mass flow rate shown in Figure 12 is at a simila level to the one depicted in Figure 8 even though the effect of disturbances is considered

Conclusions
In this paper, the authors proposed a new flying-type firefighting system.The proposed system is actuated with a water jet system and controlled for rapid and accurate fire suppression.This paper focused mainly on the dynamic modeling of the novel flying-type firefighting system with water jet actuators and a 6-DoF motion control system design.
The detailed conceptual model of the proposed system that can fly into the flames and suppress the fire was introduced, and the 6-DoF system dynamics were derived with a water-powered propulsion system.The motion control system is designed based on the SMC technique and consists of two components.The first control component is the altitude and position control system, and the second control component is the attitude control, with a part of the desired attitude state calculated with the altitude and position controller.On another note, a PID controller was introduced to carry out a comparative study of the control performance of both control systems.
The simulations were conducted to show the usefulness of the proposed system, the accuracy of the identified parameters, and the efficiency of the motion control system.The simulation tests were conducted in the absence and the presence of disturbances.The SMC showed its superiority over the PID control system in both simulation scenarios.These simulations confirmed that the mass flow rate of the water jet actuator, identi-

Conclusions
In this paper, the authors proposed a new flying-type firefighting system.The proposed system is actuated with a water jet system and controlled for rapid and accurate fire suppression.This paper focused mainly on the dynamic modeling of the novel flying-type firefighting system with water jet actuators and a 6-DoF motion control system design.
The detailed conceptual model of the proposed system that can fly into the flames and suppress the fire was introduced, and the 6-DoF system dynamics were derived with a water-powered propulsion system.The motion control system is designed based on the SMC technique and consists of two components.The first control component is the altitude and position control system, and the second control component is the attitude control, with a part of the desired attitude state calculated with the altitude and position controller.On another note, a PID controller was introduced to carry out a comparative study of the control performance of both control systems.
The simulations were conducted to show the usefulness of the proposed system, the accuracy of the identified parameters, and the efficiency of the motion control system.The

Actuators 2021 , 25 Figure 1 .
Figure 1.Conceptual drawing of flying-type firefighting system with the water jet actuator.
Figure 2. The location of the body-fixed frame (b-frame, b b b b

Figure 1 .
Figure 1.Conceptual drawing of flying-type firefighting system with the water jet actuator.
schematic drawing of the 6-DoF flying type firefighting system with four wat jet actuators is illustrated in Figure 3.The nozzles of the water jet actuators are connec ed to the sides of the head part at an angle of δ from the b z -axis and ( 1, 4 i F i =  represent the force of each actuator.

Figure 2 .
Figure 2. Reference-coordinated description of flying type firefighting system with water jet act ator.

Figure 2 .Figure 3 .
Figure 2. Reference-coordinated description of flying type firefighting system with water jet actuator.

Figure 3 .
Figure 3. Head part schematic drawing of the proposed system.

Assumption 1 .Assumption 3 .
The water flows are steady, thus ∂ ∂t CV ρVdV = 0. Assumption 2. The flows have uniform velocity profiles, so CS The viscosity of the water is neglected in this paper.

Assumption 4 .
The elements of the vector F Dxb F Dyb F Dzb τ Dxb τ Dyb τ Dzb T are bounded.

Figure 4 .
Figure 4. Block diagram of the motion control system for the flying-type firefighting system with water jet actuators.

Figure 4 .
Figure 4. Block diagram of the motion control system for the flying-type firefighting system with water jet actuators.

Actuators 2021 , 25 Figure 5 .
Figure 5. Trajectory tracking performance of the flying-type firefighting system in three-dimensional view.

Figure 7 .Figure 7 . 25 Figure 8 .
Figure 7. Attitude control response of two motion control systems.The red dotted line is PID con trol; the solid blue line is SMC; the reference trajectory of the yaw motion is

Figure 8 .
Figure 8.The mass flow rate of each nozzle.As for the second simulation scenario, the disturbances in Equation (43) are applied.The results are shown in Figures 9-12.The position and attitude time response is depicted in Figures 9-11.Similar to the simulation results without disturbances, both controllers can lead the head part along the reference trajectory.However, residual vibrations appear in the errors.The mass flow rate shown in Figure12is at a similar level to the one depicted in Figure8even though the effect of disturbances is considered.
. The position and attitude time response is depicted in Figures 9-11.Similar to the simulation results without disturbances, both controllers can lead the head part along the reference trajectory.However, residual vibrations appear in the errors.The mass flow rate shown in Figure

Figure 9 .
Figure 9. Trajectory tracking performance of the flying type firefighting system with disturbance in three-dimensional view.

Figure 9 .
Figure 9. Trajectory tracking performance of the flying type firefighting system with disturbances in three-dimensional view.

Figure 10 .Figure 10 .
Figure 10.Position control response of two motion control systems in the presence of disturbances.The red dotted line is PID control; the solid blue line is SMC; the black dashed reference trajectory line along the z -direction is

Figure 11 .
Figure 11.Position control response of two motion control systems in the presence of disturbances.The red dotted line is PID control; the solid blue line is SMC; the reference trajectory of the yaw motion is ϕ = −0.05 .

Figure 11 . 25 Figure 12 .
Figure 11.Position control response of two motion control systems in the presence of disturbances.The red dotted line is PID control; the solid blue line is SMC; the reference trajectory of the yaw motion is ϕ d = −0.05t.

Figure 12 .
Figure 12.The mass flow rate of each nozzle.

Table 1 .
Parameters of the flying-type firefighting system used in the simulation study.

Table 2 .
Controller gains used in the simulation.