Next Article in Journal
Files Cooperative Caching Strategy Based on Physical Layer Security for Air-to-Ground Integrated IoV
Previous Article in Journal
FANET Routing Protocol Analysis for Multi-UAV-Based Reconnaissance Mobility Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Motion Control System Design for a Novel Water-Powered Aerial System for Firefighting with Flow-Regulating Actuators

1
Department of Mechanical System Engineering, Division of Energy Transport System Engineering, Pukyong National University, Busan 48513, Republic of Korea
2
Department of Chassis and Body, Faculty of Vehicle and Energy Engineering, Ho Chi Minh City University of Technology and Education, Ho Chi Minh City 700000, Vietnam
3
Department of Intelligent Robot Engineering, Pukyong National University, Busan 48513, Republic of Korea
*
Author to whom correspondence should be addressed.
Drones 2023, 7(3), 162; https://doi.org/10.3390/drones7030162
Submission received: 15 January 2023 / Revised: 21 February 2023 / Accepted: 24 February 2023 / Published: 26 February 2023
(This article belongs to the Section Drone Design and Development)

Abstract

:
Flying water-jet propulsion devices, such as jet boards, jet packs, and jet bikes, can execute complex flight maneuvers. However, they require the direct involvement of trained operators to control, and their applications are very limited. In this study, we design an effective controller for a novel water-powered aerial system that aims for autonomous firefighting missions, especially at or in bodies water. Unlike existing water-powered systems, an assembly of flow-regulating actuators is proposed to fully operate the system in three-dimensional space. The paper first formulates the system dynamics by coupled partial ordinary differential equations. Then, the nonlinear controller is designed to ensure the desired system motion and stability. The design takes distinct characteristics of the system, such as coupling, under actuation, and effects of the hose conveying the water, into consideration so that the system is stabilized and uniform ultimate boundedness is achieved. Computational studies in comparison with previous control methods validated the superiority and feasibility of the proposed control system.

1. Introduction

These days, water-powered aerial systems are becoming popular. Especially in hydroflight sports, it is very likely that readers have seen spectacular performances performed by these systems either online or in person. Structurally, each of these systems features three main components: a jet ski that pumps water, a water-conveying hose, and a head assembly that carries a rider. In general, their operation basics is that the water jetted from the nozzles underneath the head assembly generates thrusts and lifts off. The rider controls the system either by leaning their body (as is the case for flyboards [1,2] and hoverboards [3]) or by rotating the direction of the water nozzle (as is the case for jet packs [4] and jet bikes [3,5]) to ride. Thereby, the rider can perform complex movements, such as turning, flipping, or spinning. With this maneuverability, water-powered systems have the potential for many aerial operations. However, their inevitable drawbacks make them currently only applicable to hydroflight sports or show performances. Firstly, the systems require a substantial amount of water for propulsion. Secondly, their flight distance is limited by the length of the water-conveying hose. Especially, most of the water-powered systems are manual or semiautomated such that only trained and skilled people are able to operate them. Newly developed products can operate autonomously, thanks to the use of sensors, actuators, and automatic controllers. For example, the Zapata Aquadrone [3] is equipped with an electric pump and four electric motors rotating four thrust nozzles for self-manipulation. Nevertheless, it is designed primarily for demonstration purposes.
Water-powered aerial systems have shown potential in firefighting missions, particularly in scenarios where fires occur in water-based environments, such as crowded harbors, bridges with vehicle traffic, or boats in the middle of the sea. These scenarios present a significant challenge due to limited accessibility and potential safety hazards. For example, the current method used to put out a boat fire on the sea is demonstrated in Figure 1. It involves the use of one or more vessels to approach the scene and spray water onto the burning boat, which can be inaccurate and potentially harmful to property and people still on the scene. Moreover, terrestrial robots are not deployable, and the firefighting methods of using commercial aerial vehicles, such as helicopters and fixed-wing aircraft, also face similar challenges in these scenarios. Nevertheless, small-scale aerial systems offer several advantages, including aerial perspective, remote control capabilities, and access to difficult-to-reach areas. Hence, drones can carry fire retardants [6,7,8,9] or attach with water-conveying hoses and fire nozzles [10,11,12,13] to battle the fire. Cooperative use of multiple drones, as in [14,15], can cover a large object. However, their adoption is still limited due to multiple constraints, such as their limited load-carrying capacity, the potential for generated airflow to exacerbate the fire, and the presence of dangerous high-speed rotating blades.
Meanwhile, in these scenarios, the water supply concern for water-powered systems is no longer a limiting factor. The systems can transport to get close to the accident areas, and their head assembly only needs to take off from there, thus overcoming their limitation on flight distance. By flying above the flames and taking advantage of high ground, water-powered systems can effectively extinguish fires. One demonstration is that firefighters in Dubai in 2007 used a flyboard to bring them to the height needed to extinguish fires on a bridge over the sea [16,17]. However, firefighters need skills to operate these because of the manual systems. Additionally, firefighting at a close distance puts a threat to firefighters’ lives. In [18,19,20,21], Dragon Firefighter, an aerial robot with a water-conveying hose and water-jet propulsion, was developed for firefighting. The robot features four rotatable nozzles that are controlled to maneuver the robot itself and suppress the fire. Experimental studies have demonstrated the robot’s successful and safe suppression of flames in some simple scenarios but have also revealed its limitations in reachability and maneuverability. Therefore, the development of a different method for actuation and control is needed to enable the autonomous operation of water-powered systems.
Three different actuating methods have been proposed in the literature for controlling the motion of water-powered aerial systems. The first method involves altering the weight distribution of the head assembly in order to change its attitude angles, such as the weight-shifting mechanism proposed by Tri et al. [22]. This mimics the way a person rides his/her flyboard or hoverboard. Unfortunately, this method has some drawbacks, including an increase in weight, nonminimum phase phenomena, and uncontrollable yaw rotation. The second, and the most common actuating method, is to rotate the thrust nozzles only instead of the whole head assembly to change the direction of the thrust force generated by each nozzle. Liu et al. [23] used two actuators to drive a nozzle rotating 360 degrees, two nozzles in total, for their unmanned water-powered aerial vehicle. The flying robots with water jets in [24,25,26,27,28] used a module consisting of two fixed nozzles and two biaxial active nozzles. Additionally, the Zapata Aquadrone [3] is equipped with four rotatable thrust nozzles. These systems are fully controllable in three-dimensional space but come with several drawbacks, such as nonminimum phase characteristics [23], very limited operational space in comparison to firefighting scenarios [24,25,26,27,28], or complex structure and high requirement of waterproofing [3], etc. Finally, the last method is to regulate the water flow rate out of the thrust nozzles, for example, the flow-regulated, unmanned, water-powered, aerial vehicle in [23]. Thereby, the amplitude of the thrust forces is controllable. This mechanism is simpler and easier to control in comparison to the two other methods. The available design in [23] is also not possible to control the yaw rotation though.
In terms of control strategy, several linear ones have been proposed to achieve the desired motion of the water-powered aerial systems, for example, proportional control with speed feedback in [24,25,26,27,28], proportional–derivative control in [23], and linear quadratic integrator in [22]. These control systems are simple and only require feedback from the motion of the head assembly. However, simulations and experiments show performance is poor and greatly influenced by the water-conveying hose movement. Controllers similar to those for multi-propeller drones have also been suggested, for example, cascade proportional–integral–derivative control and cascade sliding mode control [29]. Furthermore, numerous advanced control schemes have been developed in the literature and considered different aspects of aerial systems, such as decoupled control [30] and cooperative control of multiple systems [31,32]. However, it is worth noting that the hose dynamics are not fully represented in the model of the water-powered systems in the mentioned studies, which does not help validate the control performance and preserve control system stability. For instance, [22,29] only considered that the gravitation of the hose varies with the altitude of the system. The force from the hose was formulated via a heuristic assumption in [23]. In [26,27,28], the hose was approximated by multiple links connected by elastic joints and dampers.
The hose in the abovementioned systems can be considered as a cantilevered pipe that conveys fluid under deformation. Hence, a dynamical model of it, and subsequently of the whole system, can be obtained. Theoretical models of cantilevered pipes were introduced in [33,34,35] in the form of nonlinear partial differential equations. Although the resulting models are typically complex, they can be used to determine the pipe’s deformation and vibration characteristics. Thereby, a control system that takes the partial states of the hose into consideration is able to achieve effective and robust motion performances. He et al. [36,37,38,39,40] applied this principle to stabilize flexible marine risers, which share some characteristics with a water-conveying hose. In detail, their proposed controllers need to know up to third-order partial derivatives of the riser at the control position. However, these values are challenging to obtain in practice. Alternatively, the conditions for the hose stability can be deduced from its model, allowing a control system to be designed with only the feedback from the head assembly while still preserving the stability of the whole system. Ambe et al. [24] derived the model of their complete water-jet system in the form of partial differential equations, but only for small deformation of the hose in the lateral direction. The authors proposed a proportional–derivative controller for the motion of the head and ensured energy dissipation of the system. In [41], it is shown that if the rate of water flow inside the hose is below a critical value, the water-conveying hose can be kept stable, and the total effect of the hose on the head assembly is bounded. Hence, a state-feedback controller that considers the saturation of the control inputs can preserve the system’s stability. Unfortunately, the designed control law is still insufficient in dealing with disturbances effectively.
Therefore, in this study, we first introduce a water-powered aerial system with a novel flow-regulating mechanism that makes it fully controllable in three-dimensional space. The system is aimed at firefighting missions, especially at or in bodies of water. It also can be extended for land missions by using the water source from the available fire hydrant system. The dynamics of the designed system are derived in the form of coupled partial ordinary differential equations. We then proposed a composite controller that ensures the desired altitude and attitude of the head while still preserving the overall stabilization of the whole system. Distinct characteristics of the system, such as input coupling, saturation, and motion of the hose, were considered in the control system design process. The result is that uniform ultimate boundedness is achieved. Comparative studies with two previous control strategies validated the efficiency and robustness of the proposed system. Hence, the contributions of the paper can be summarized as follows:
  • A novel water-powered aerial system with a flow-regulating mechanism is introduced and modeled.
  • A composite controller is proposed to fully drive the system in three-dimensional space and preserve the system’s stability. The uniform ultimate boundedness of the control system is proven.
  • Computational studies in comparison with two previous control strategies validated the efficiency and robustness of the proposed control system.
Accordingly, the remainder of the paper is organized as follows. Section 2 describes the design and operation of the system. Mathematical models of the water-conveying hose, the head assembly, and the whole system are derived in Section 3. The proposed controller is presented in Section 4. In Section 5, the comparative results are shown and discussed. Finally, conclusions are drawn in Section 6.

2. System Description

As indicated by its name, the proposed system uses water as a source of propulsion and firefighting. At or in bodies of water, the enormous amount of water is used. In other fire sites, water can be taken from fire hydrants. The proposed system utilizes a water pump that pulls water through a flexible hose and delivers it to a head assembly, where it flows through a manifold before being distributed to nozzles for generating thrust and performing maneuvering and firefighting actions. The head assembly, which contains major components of the system, such as nozzles, landing gear, instruments, actuators, etc., connects to the hose via a swivel ball joint to isolate the torques from the hose. The head is designed such that the water inlet is located at its mass center, while the water-conveying hose has a uniform cross-section to ensure a consistent flow rate at any point. The design of the water-powered firefighting system is illustrated in Figure 2a.
On the head assembly, there are four nozzles for generating thrust. Each nozzle has an outlet with a much smaller cross-sectional area than the head inlet. This structure accelerates the water flow, hence increasing its momentum and generating thrust. The four nozzles are numbered 1 to 4, with nozzle 1 being located at the front left and the remaining nozzles following in a counterclockwise direction. Each nozzle injects water downward and slightly out to the lateral sides, as shown in the figure. This arrangement gives several advantages over the other configuration in [1,3,4,5,23,42], including increased flight stability, ability to achieve yaw rotation, and decoupling of control input in the longitudinal and pitch directions. The system also features another nozzle located at the front of the head that sprays the water forward in fine droplets to cover and suppress a flame. The posture of this fire nozzle is slightly downward, and its axial line goes through the head’s center to prevent unwanted torque generation.
Three motorized valves, labeled as a, b, and c, regulate the water flow distribution among the four thrust nozzles. An additional on/off valve, labeled as valve d, is used to fully open or close the fire hose nozzle. Apparently, the fire hose nozzle only needs to open after approaching the fire. The corresponding flow diagram is shown in Figure 2b. At the neutral position of the motorized valves, the water flow is equally distributed among the four thrust nozzles. Manipulating one or more of these valves results in a difference in the flow rates and thus creates torques. For instance, opening valve a toward valve b increases the water flow rate and thrust of nozzles 1 and 3 and decreases those of nozzles 2 and 4. This results in a torque around the vertical axis passing through the center of the head, thereby producing a yaw rotation.

3. System Modeling

In order to realize the maneuverability of the proposed system, its mathematical model is derived in this section. Initially, a system of coordinates is defined, and the kinematic relationships among them are derived. Next, the nonlinear dynamics of the water-conveying hose and the head part are considered. Finally, a governor equation with initial and boundary conditions represents the completed model of the system.

3.1. Coordinates and System Kinematics

Two Cartesian coordinates and a curvilinear one are defined regarding the motion of the water hose and the system head. One Cartesian coordinate is the Earth-fixed reference frame XYZ, and the other is the body-fixed frame of the head assembly XbYbZb. The latter is placed at the mass center of the head with the water inside it. The Xb-, Yb-, and Zb-axes follow the longitudinal, lateral, and vertical of the head, with the Xb-axis pointing forward and the Zb-axis pointing upward. The curvilinear coordinate starts at the pump outlet and is along the center line of the water-conveying hose. This is depicted in Figure 3. Forces acting on the head part are also shown in this figure.
In the reference frame, the hose element at a coordinate s at a time t can be described by the partial vector X ( s , t ) = X ( s , t ) Y ( s , t ) Z ( s , t ) T . Let the dots be denoted for time derivatives, and the primes stand for derivatives with respect to coordinate s. Then, the velocity of a hose element is X ˙ ( s , t ) , and the velocity of the water flowing through that is computed by:
V f ( s , t ) = X ˙ ( s , t ) + V h ( t ) X ( s , t )
where V h ( t ) is the axial velocity of the water flow in the hose.
It is noted that X ( 0 , t ) defines the position of the hose inlet, and X ( L , t ) , with L being the length of the hose, is for the hose outlet. In other words, X ( L , t ) is also the position of the head assembly. Additionally, the completed motion of the head can be represented by:
X ˙ ( L , t ) φ ˙ h ( t ) = R E b ( t ) O 3 O 3 R E u l e r b ( t ) v b ( t ) ω b ( t )
In which:
  • v b ( t ) = v x b ( t ) v y b ( t ) v z b ( t ) T : vector of translational velocity with respect to the body-fixed frame.
  • φ ( t ) = ϕ ( t ) θ ( t ) ψ ( t ) T : vector of three Euler angles in XYZ order.
  • ω b ( t ) = ω x ( t ) ω y ( t ) ω z ( t ) : vector of angular velocities.
  • R E b ( t ) = c ψ c θ c ψ s θ s ϕ s ψ c ϕ c ψ s θ c ϕ + s ψ s ϕ s ψ c θ s ψ s θ s ϕ + c ψ c ϕ s ψ s θ c ϕ c ψ s ϕ s θ c θ s ϕ c θ c ϕ , with c ( ) = cos ( ) and s ( ) = sin ( ) , are the transformation matrix from the head-fixed frame to the reference frame. R E u l e r b ( t ) = 1 c θ c θ s ϕ s θ c ϕ s θ 0 c ϕ c θ s ϕ c θ 0 s ϕ c ϕ is that from the angular velocities to the Euler angles’ rate.
  • O m × n : m-by-n zero matrix.
For the water flow at an inlet/outlet of the head assembly, its velocity can be obtained as follows:
V j b ( t ) = V j ( t ) n ^ j b ( t ) + ω b ( t ) × r j b + v b ( t )
with j = h , l , i i = 1 ~ 4 indicating the head inlet, the outlet of the fire hose nozzle, and the outlet of the ith thrust nozzle, respectively. The second subscript b means that the value is in the body-fixed frame of the head assembly. r j b is the position vector of the j inlet/outlet, and n ^ j b ( t ) is the unit-length normal vector of the corresponding cross-sectional area.
Remark 1:
X ( s , t ) is the vector tangent with the hose center line at coordinate s, so it is also the normal vector of the cross-section at this location. Subsequently, X ( L , t ) = n ^ h ( t ) .

3.2. Water-Conveying Hose Model

In order to realize the dynamics of the water-conveying hose, let the hose be isolated from the head assembly and replaced by corresponding reaction forces. Then, it can be considered that the hose has a free end. A modified version of Hamilton’s principle for the assembly of the hose and the water inside is written as follows:
t 1 t 2 δ W h ( t ) d t + δ t 1 t 2 L h ( t ) d t = t 1 t 2 m f V h X ˙ ( L , t ) + V h ( t ) X ( L , t ) δ X ( L , t ) d t
where δ is the Hamilton’s variational operator. δ W s ( t ) is the virtual work performed by external forces on the hose. They are listed as follows:
δ W s = 0 L F ( s , t ) T δ X Δ ( s , t ) d s , X Δ ( s , t ) = X ( s , t ) X ( s , 0 ) , F ( s , t ) = μ s X ˙ ( s , t ) + F r ( s , t ) δ Z ( s , t ) A h / π F p ( t ) δ s L + F 0 ( t ) δ s p h A h n ^ h ( t ) δ s L
In which μ s X ˙ ( s , t ) , with μ s being a matrix of positive definite coefficients, is the damping on the hose. F p and F 0 are the reaction forces due to isolating the head and the total force holding the hose’s first end. F r is the total interactive force between the hose and its surroundings. For the sake of simplicity, complex aerodynamic forces, such as drag, wind, ground effect, etc., were neglected. Hence, F r contains the contact forces (normal and frictional forces) with the ground in onshore applications and includes buoyance force, hydrodynamic force, and seabed interaction force in offshore applications. p h A h is the magnitude of the force due to the exit pressure at the hose outlet. The Dirac delta functions δ indicate the locations that the forces act on, i.e., ( ) = 0 .
The second term in Equation (4) is the Lagrangian including kinetic energy T and potential energy thanks to the hose’s bending strain B and gravitation G. That is:
L s ( t ) = T ( t ) B ( t ) G ( t )
In detail, they are obtained as follows:
T ( t ) = 1 2 m p 0 L X ˙ ( s , t ) T X ˙ ( s , t ) d s + 1 2 m f 0 L V f ( s , t ) T V f ( s , t ) d s , B ( t ) = 1 2 E I 0 L κ 2 ( s , t ) d s , κ 2 ( s , t ) = X ( s , t ) s T X ( s , t ) , G ( s , t ) = m p + m f g 0 L Z ( s , t ) d s
where m p and m f are the masses per unit length of the hose and the water inside, respectively. E I is the flexural rigidity of the hose, and κ is the curvature of the deformed centerline. g is the gravitational acceleration.
The right-hand term in Equation (4) is the virtual momentum transport across the open surface at the end of the hose. Additionally, we have the following relations with the mass flow rate of the water m ˙ h ( t ) inside the hose:
V h ( t ) = m ˙ h ( t ) / ρ A h , m f = ρ A h
with ρ being the water density and A h being the hose’s cross-sectional area.
Substitute Equations (5)–(8) into Equation (4) and neglect the acceleration/deceleration of the water flow inside the hose. After a complicated but straightforward calculation, the governor equation of the water-conveying hose is obtained:
m p + ρ A h X ¨ ( s , t ) = μ s X ˙ ( s , t ) m ˙ h 2 ( t ) X ( s , t ) / ρ A h 2 m ˙ h ( t ) X ˙ ( s , t ) E I . X ( s , t ) + F r ( s , t ) δ Z s A h / π + m p + ρ A h g
with the boundary and initial conditions given by:
F p ( t ) + p h A h n ^ h ( t ) E I . X ( L , t ) = 0 , F 0 ( t ) + m f V h ( t ) X ˙ ( 0 , t ) + E I . X ( 0 , t ) = 0 , X ( L , t ) = 0 , X ( 0 , t ) = 0
Remark 2:
The three-dimensional dynamics (9) of the water-conveying hose can be reduced to two partial motions thanks to the inextensible condition of the hose as follows [33,34,35]:
X ( s , t ) T X ( s , t ) = 1 , or d d t X ( s , t ) T X ( s , t ) = 0  
Remark 3:
Consider the system operating in onshore missions, and the interactive force F r consists of the normal contact force with the ground and frictional forces generated by the motion of the hose on the ground surface. They can be modeled by:
F r ( s , t ) = g m p + ρ A h μ s i g n X ˙ ( s , t ) μ s i g n Y ˙ ( s , t ) 1
with μ being the frictional coefficient and s i g n ( ) being the signum functional.

3.3. Head Assembly Model

For the sake of modeling, the head assembly with its contained water was treated as a rigid body. This was performed thanks to the isolation of the hose and the nondeforming control volume from the water inlet to the nozzle outlet. The Newton–Euler formulation gives the dynamics of the head with respect to the body-fixed frame as follows:
M h v ˙ b ( t ) + ω b ( t ) × v b ( t ) = F h b ( t ) + M h g b + p h A h n ^ h b ( t ) + F p b ( t ) , J h ω ˙ b ( t ) + ω b ( t ) × J h ω b ( t ) = τ h b ( t )
M h is the mass of the head assembly, including its contained water. J h is its matrix of inertia. The cross-product terms appear thanks to time derivatives with respect to the rotating frame. The total force acting on the head is due to the momentum change of the water flowing through it F h b ( t ) ; the gravitation M h g b ; the water pressure on the inlet cross-section p h A h n ^ h b ( t ) ; and the reaction due to the hose isolation F p b ( t ) . Meanwhile, the torque acting on the head is only due to the thrust at the nozzles’ outlets τ h b , because the water inlet is at the mass center of the head. The force F p b ( t ) and the torque τ h b were computed by:
F h b ( t ) = i = 1 4 m ˙ i ( t ) V i b ( t ) m ˙ l ( t ) V l b ( t ) + m ˙ h ( t ) V h b ( t ) , τ h b ( t ) = i = 1 4 r i b × m ˙ i ( t ) V i b ( t ) r l b × m ˙ l ( t ) V l b ( t )
Moreover, let us define a mapping variable u j 1 , 1 , j = a , b , c , between the opening of the three-way valve j and the resulting flow rates after it. u j = 0 means that valve j is at a neutral position and the water is distributed equally to its two outlets. u j = 1 indicates that the valve is fully opened to the left side in the diagram in Figure 1b. For example, u a = 1 means all the water entering valve a will pass to valve b, and if u b = 1 , the water entering valve b will flow to nozzle 1 completely. In contrast, u j = 1 represents that the valve is fully opened to the right side in the diagram. Hence, the complete mapping can be written as follows:
m ˙ 1 ( t ) = 1 4 1 + u b ( t ) 1 + u a ( t ) m ˙ h ( t ) m ˙ l ( t ) , m ˙ 3 ( t ) = 1 4 1 u b ( t ) 1 + u a ( t ) m ˙ h ( t ) m ˙ l ( t ) , m ˙ 2 ( t ) = 1 4 1 + u c ( t ) 1 u a ( t ) m ˙ h ( t ) m ˙ l ( t ) , m ˙ 4 ( t ) = 1 4 1 u c ( t ) 1 u a ( t ) m ˙ h ( t ) m ˙ l ( t ) ,
In order to obtain the representation of the head dynamics in the reference frame, substitute Equations (14) and (15), the boundary conditions (10), and the kinematics (1) into Equation (13). By neglecting higher-order terms and applying the approximation of small pitch and roll angles, the mathematical model of the head assembly is derived as in Equations (16) and (17).
X ¨ ( L , t ) = 1 M h R E b ( t ) 0 s δ w c δ h s δ U Y ( t ) U Z ( t ) + 1 ρ A l c δ l 0 s δ l m ˙ l 2 ( t ) + h m ˙ h m ˙ l 0 0 0 h m ˙ h m ˙ l 0 0 0 0 R E u l e r 1 ( t ) φ ˙ ( t ) + g + 1 M h m ˙ h 2 ( t ) ρ A n n ^ h ( t ) + 1 M h E I . X ( L , t ) + d X ( t ) φ ¨ ( t ) = R E u l e r ( t ) J h 1 U Y ( t ) U X ( t ) U ψ ( t ) + 0 h l c δ l l l s δ l ρ A l m ˙ l 2 ( t ) 0 + R ˙ E u l e r ( t ) ω b ( t ) + d φ ( t )
with
U X ( t ) = l c δ 4 A n ρ m ˙ h ( t ) m ˙ l ( t ) 2 u b ( t ) u c ( t ) , U Y ( t ) = 1 4 A n ρ w c δ h s δ m ˙ h ( t ) m ˙ l ( t ) 2 u b ( t ) + u c ( t ) , U Z ( t ) = c δ 4 A n ρ m ˙ h ( t ) m ˙ l ( t ) 2 , U ψ ( t ) = l s δ 2 A n ρ m ˙ h ( t ) m ˙ l ( t ) 2 u a ( t )
In these equations, l, w, and h are the distances in the longitudinal, lateral, and vertical directions, from the coordinate origin to the outlets of the thrust nozzles, respectively. Similarly, ll, wl, and hl are those regarding the fire extinguishing nozzle. δ is the magnitude of the acute angle between the thrust nozzles and the vertical axis. δ l is the acute angle between the fire extinguishing nozzle and the longitudinal axis Xb. Finally, d X ( t ) and d φ ( t ) are the remaining disturbances and higher-order terms acting on the head’s translational and rotational movements, respectively.
Remark 4:
The system is underactuated because four control inputs control the 6-DoF motion of the head, as well as maintain the stability of the flexible hose. Moreover, there is a coupling of input in the lateral motion and roll rotation, with respect to the body-fixed coordinate, of the head.

4. Control System Design

In a firefighting mission, the head assembly follows a given trajectory to approach its working position at a proper distance to a flame and then sprays the water to suppress the flame. During the process, the motion of the head should be kept stabilized, and the fluctuation of the water-conveying hose should be minimized. Furthermore, because the system is underactuated, it is not possible to control every DoF independently. In fact, the head’s longitudinal and lateral translations are coupled with its pitch and roll rotations, respectively. Hence, a composite control law is designed in this section and takes the above characteristics into account. Because it is difficult to obtain the partial states of the hose in practical operation, the proposed control law uses the feedback from the head assembly and the water supply pump. Additionally, a command filter is introduced to ensure the proposed control law is applicable.

4.1. Control Law Design

The control law to be designed consists of the pump’s required flow rate m ˙ h ( t ) , the opening of the valves u j 1 , 1 , j = a , b , c , and two virtual functions ϕ d ( t ) and θ d ( t ) of the head’s desired roll and pitch rotations. For the system head, the references of the trajectory and attitude are given by X d = X d Y d Z d T and ψ d ( t ) . Firstly, filtered references were introduced to guarantee they are twice differentiable.
X ^ ¨ d ( t ) φ ^ ¨ d ( t ) + 2 ζ ω n X ^ ˙ d ( t ) φ ^ ˙ d ( t ) + ω n 2 e f ( t ) = 0 , e f ( t ) = X d ( t ) φ d ( t ) X ^ d ( t ) φ ^ d ( t )
with φ d ( t ) = ϕ d ( t ) θ d ( t ) ψ d ( t ) T . In Equation (18), the hats indicate the filtered values, and e f ( t ) is the vector of difference between them and the unfiltered values. ζ is the diagonal matrix of damping ratios, and ω n is that of the natural frequencies of the filters. Then, a vector of tracking errors of the head was denoted with respect to the newly filtered references as follows:
e 1 ( t ) = X ( L , t ) φ ( t ) X ^ d ( t ) φ ^ d ( t )
To compensate for this error, a virtual control law for the head’s velocity was designed as in (20):
η 1 ( t ) = X ^ ˙ d ( t ) φ ^ ˙ d ( t ) κ 1 e 1 ( t ) κ 2 f 1 ( t )
with f 1 ( t ) = 0 t e 1 ι d ι , κ 1 and κ 2 are the positive definite and diagonal matrices. In addition, a composite error was considered as follows:
e 2 ( t ) = I 3 β ( t ) R E u l e r 1 ( t ) O 3 I 3 X ˙ ( L , t ) φ ˙ ( t ) η 1 ( t ) e u ( t )
with In and On being the identity matrix and zero matrix, respectively, of size n. In addition to the difference between the actual and desired velocities of the head assembly (i.e., X ˙ T ( L , t ) φ ˙ T ( t ) T and η 1 ( t ) ), this error representation also takes into account coupling and the saturation of the control inputs. In detail, β ( t ) in Equation (21) is aimed to decouple the control input U Y ( t ) that appears in the lateral motion. e u ( t ) is the first-order filtered error due to saturated control inputs. They are given by:
β ( t ) = J xx s δ M h w c δ h s δ c ϕ s ψ c ψ s ϕ s θ 0 0 c ϕ c ψ s ψ s ϕ s θ 0 0 s ϕ c θ 0 0 , e ˙ u ( t ) = κ u e u ( t ) e 1 ( t ) κ s f 1 ( t ) + Δ u X ( t ) Δ u φ ( t ) , Δ u X ( t ) = 1 M h R E b ( t ) 0 0 U Z ( t ) R ¯ E b ( t ) 0 0 U ¯ Z ( t ) , Δ u φ ( t ) = R E u l e r ( t ) J h 1 U Y ( t ) U ¯ Y ( t ) U X ( t ) + U ¯ X ( t ) U ψ ( t ) + U ¯ ψ ( t )
κ u and κ s in the dynamics of e u ( t ) are positive definite diagonal matrices. κ s is chosen such that κ 2 κ 1 κ s κ s = 0 . The bar above a variable indicates the computed values from the controller, i.e., without saturations, to distinguish from the true value (without the bars) performed on the system.
Then, the dynamics of the error e 2 ( t ) can be obtained by substituting the head dynamics in Equation (16) and the terms in Equation (22) into its time derivative:
e ˙ 2 ( t ) = U ¯ ( t ) + F ω ( t ) τ ω ( t ) ω b ( t ) + F l ( t ) τ l ( t ) m ˙ l 2 ( t ) + g + 1 M h m ˙ h 2 ( t ) ρ A n n ^ h ( t ) + 1 M h E I . X ( L , t ) O 3 × 1 + d X ( t ) β ( t ) R E u l e r 1 ( t ) d φ ( t ) d φ ( t ) X ^ ¨ d ( t ) φ ^ ¨ d ( t ) + κ 1 e 2 ( t ) κ 1 κ 1 κ 2 e 1 ( t ) κ 1 κ 2 f 1 ( t ) + κ 1 + κ u e u ( t )
with
U ¯ ( t ) = 1 M h R ¯ E b ( t ) O 3 O 3 R E u l e r ( t ) J h 1 O 1 × 2 U ¯ Z ( t ) U ¯ Y ( t ) U ¯ X ( t ) U ¯ ψ ( t ) T , F ω ( t ) = 1 M h R E b ( t ) h m ˙ h m ˙ l 0 0 0 h m ˙ h m ˙ l 0 0 0 0 + β ( t ) β ˙ ( t ) , F l ( t ) = 1 M h ρ A l R E b ( t ) c δ l 0 s δ l T h l c δ l l l s δ l ρ A l β ( t ) J h 1 0 1 1 T , τ ω ( t ) = R ˙ E u l e r ( t ) , τ l ( t ) = h l c δ l l l s δ l ρ A l R E u l e r ( t ) J h 1 0 1 1 T
Then, a saturation-free control law was designed as follows:
U ¯ ( t ) = X ^ ¨ d ( t ) φ ^ ¨ d ( t ) F ω ( t ) τ ω ( t ) R E u l e r 1 φ ˙ ( t ) F l ( t ) τ l ( t ) m ˙ l 2 ( t ) g O 3 × 1 + κ 1 κ 1 κ 2 2 I e 1 ( t ) + κ 1 κ 2 2 κ s f 1 ( t ) κ 1 + κ 3 e 2 ( t ) κ 1 + κ u e u ( t ) + O 1 × 2 m p + ρ A h g M h ς ^ ( t ) Z L , t A h / π O 1 × 3 T , ς ^ ˙ ( t ) = 1 κ ς Z L , t A h / π e Z 2 ( t )
with κ 3 as positive definite and diagonal matrices and κ ς as a positive constant value. From the above control vector, we derived the control inputs defined in (17), as well as the required roll and pitch rotations of the head, as follows:
U ¯ Z ( t ) = M h I 3 O 3 U ¯ ( t ) 2 , ϕ d ( t ) = a r c s i n M h U ¯ Z ( t ) s ψ c ψ O 1 × 4 U ¯ ( t ) , θ d ( t ) = a r c t a n M h c ψ s ψ O 1 × 4 U ¯ ( t ) O 1 × 2 1 O 1 × 3 U ¯ ( t ) , U ¯ Y ( t ) U ¯ X ( t ) U ¯ ψ ( t ) = J h R E u l e r 1 ( t ) O 3 I 3 U ¯ ( t )
As a result, the required flow rate from the pump and the opening of the flow-regulating valves were obtained by mapping (27). The overall structure of the proposed control system is depicted in Figure 4.
m ˙ ¯ h ( t ) = m ˙ l ( t ) + 4 A n ρ c δ U ¯ Z ( t ) , u ¯ a ( t ) = 2 A n ρ U ¯ ψ ( t ) l s δ m ˙ ¯ h ( t ) m ˙ l ( t ) 2 , u ¯ b ( t ) = 2 A n ρ w c δ h s δ U ¯ X ( t ) + l c δ U ¯ Y ( t ) l c δ w c δ h s δ m ˙ ¯ h ( t ) m ˙ l ( t ) 2 , u ¯ c ( t ) = 2 A n ρ w c δ h s δ U ¯ X ( t ) + l c δ U ¯ Y ( t ) l c δ w c δ h s δ m ˙ ¯ h ( t ) m ˙ l ( t ) 2

4.2. Stability Analysis

Theorem. 
For the water-powered aerial system whose dynamics are represented by Equations (9), (10), (16), and (17), the control law (25)–(27), with a proper choice of gains κ ς and κ i , ( i = 1 , 2 , 3 , u ) , guarantees the uniform ultimate boundedness of the entire control system.
Proof of the Theorem.
To begin with, we considered the following Lyapunov function candidate:
V ( t ) = V 1 ( t ) + V 2 ( t ) + V u ( t ) + V s 1 ( t ) + V s 2 ( t ) + V ς ( t ) + Λ ( t ) V 1 ( t ) = 1 2 s ( t ) T s ( t ) , s ( t ) = e 1 ( t ) + κ s f 1 ( t ) , V 2 ( t ) = 1 2 e 2 ( t ) T e 2 ( t ) , V u ( t ) = 1 2 e u ( t ) T e u ( t ) , V s 1 ( t ) = m p + ρ A h 2 L 2 M h 2 0 L X ˙ ( s , t ) T X ˙ ( s , t ) d s , V s 2 ( t ) = E I m p + ρ A h L 2 M h 2 0 L X ( s , t ) s T X ( s , t ) d s , V ς ( t ) = 1 2 κ ς g m p + ρ A h M h ς ^ ( t ) ς 1 2 , Λ ( t ) = m p + ρ A h M h e 2 ( t ) T 0 L X ˙ ( s , t ) d s O 3 × 1
For proof of the Theorem, the following Lemma 1 and Lemma 2 were considered.
Lemma 1:
The Lyapunov function candidate V ( t ) in (28) has the upper and lower boundaries established as follows.
0 λ min Q 1 z ( t ) 2 2 + V s 2 ( t ) + V ς ( t ) + Λ ( t ) V ( t ) λ max Q 2 z ( t ) 2 2 + V s 2 ( t ) + V ς ( t ) + Λ ( t )  
with:
z ( t ) = s ( t ) T e 2 ( t ) T e u ( t ) T 0 L X ˙ ( s , t ) 2 d s T , Q 1 = 1 2 d i a g 1 , 1 , 1 , m p + ρ A h 2 M h 2 , Q 2 = 1 2 d i a g 1 , 1 , 1 , m p + ρ A h 2 L M h 2
In Equation (30), the notations λ and 2 stand for the eigenvalue and Euclidean norm of , respectively.
Proof of Lemma 1:
From the Cauchy–Schwarz inequality and Holder’s inequality for integral, one can derive:
m p + ρ A h 2 2 M h 2 0 L X ˙ ( s , t ) 2 d s T 0 L X ˙ ( s , t ) 2 d s V s 1 ( t ) m p + ρ A h 2 L 2 M h 2 0 L X ˙ ( s , t ) 2 d s T 0 L X ˙ ( s , t ) 2 d s
Substituting inequalities (31) into Equation (28) gives the result as follows:
z ( t ) T Q 1 z ( t ) + V s 2 ( t ) + V ς ( t ) + Λ V ( t ) z ( t ) T Q 2 z ( t ) + V s 2 ( t ) + V ς ( t ) + Λ
Moreover:
V 2 ( t ) + V s 1 ( t ) + Λ 1 2 e 2 ( t ) T e 2 ( t ) + m p + ρ A h 2 2 M h 2 0 L X ˙ ( s , t ) d s T 0 L X ˙ ( s , t ) d s + m p + ρ A h M h e 2 ( t ) T 0 L X ˙ ( s , t ) d s O 3 × 1 = 1 2 e 2 ( t ) + 0 L X ˙ ( s , t ) d s O 3 × 1 T e 2 ( t ) + 0 L X ˙ ( s , t ) d s O 3 × 1 0
Hence, Lemma 1 is obtained. □
Lemma 2:
With the proposed control law, the derivative of V ( t ) in (28) along the trajectory of the system satisfies:
V ˙ ( t ) z ( t ) T P 1 z ( t ) , z ( t ) 2 2 σ
with P 1 as a positive definite matrix, σ as a positive definite constant, and:
z ( t ) = s ( t ) T e 2 ( t ) T e u ( t ) T 0 L X ˙ ( s , t ) 2 d s T
Proof of Lemma 2:
By taking the integration of both sides of the governor Equation (9) and adopting the initial conditions of the hose, the third-order differential term E I . X ( L , t ) can be obtained as follows:
E I . X ( L , t ) = m ˙ h 2 ( t ) X ( L , t ) / ρ A h 2 m ˙ h ( t ) X ˙ ( L , t ) o L m p + ρ A h X ¨ ( s , t ) μ s X ˙ ( s , t ) F r ( s , t ) δ Z s A h / π m p + ρ A h g d s
The total weight of the hose that the system must lift increases with its height and becomes a noticeable influence. From the properties of a line integral in a vector field and by adopting the mean value theorem for definite integrals [43], this weight can be represented in the following way:
0 L m p + ρ A h g + F r Z ( s , t ) δ Z s A h / π d s = g m p + ρ A h Z L , t A h / π O 2 × 1 d s d z ( Z a , t )
with Z a ( t ) A h / π , Z L , t . Let us write d s ( Z a , t ) / d z = ς 1 + ς 2 ( t ) , where ς 1 is the nonzero constant of the mean value of d s ( Z a , t ) / d z , and ς 2 ( t ) is the remainder to make the equality (39) true.
Taking the time derivative of the Lyapunov candidate function, substituting the system’s dynamics and above control law, and taking integrals by parts result in (38).
V ˙ ( t ) = s ( t ) T κ 1 κ s s ( t ) e 2 ( t ) T κ 3 e 2 ( t ) e u ( t ) T κ u e u ( t ) m p + ρ A h L M h 2 0 L X ˙ ( s , t ) T μ s X ˙ ( s , t ) d s m ˙ h m p + ρ A h L M h X ˙ ( L , t ) T X ˙ ( L , t ) m p + ρ A h 2 L M h 2 g 0 L 1 μ X ˙ ( s , t ) + Y ˙ ( s , t ) d s + ε ( t ) ,
with
ε ( t ) = s ( t ) T β ( t ) R E u l e r 1 ( t ) φ ˙ ( t ) + e 2 ( t ) T d 2 ( t ) + e 2 ( t ) T P 2 0 L X ˙ ( s , t ) d s O 3 × 1 + 2 m ˙ h ( t ) M h e 2 ( t ) T X ˙ ( L , t ) O 3 × 1 + e u ( t ) T Δ u X ( t ) Δ u φ ( t ) + 0 L X ˙ ( s , t ) d s O 3 × 1 T d s ( t ) + m p + ρ A h M h 2 Δ X ˙ T m ˙ h 2 ( t ) ρ A h X ( L , t ) + E I . X ( L , t ) ,
and
d 2 ( t ) = d L ( t ) F ˜ r ( t ) O 3 × 1 , d L ( t ) = d X ( t ) β ( t ) R E u l e r 1 ( t ) d φ ( t ) d φ ( t ) , d s ( t ) = m p + ρ A h L 2 M h d L ( t ) O 1 × 2 m p + ρ A h 2 L g M h 2 O 1 × 3 T , Δ X ˙ = 0 L X ˙ ( s , t ) d s L X ˙ ( L , t ) , P 2 = m p + ρ A h L 2 M h κ 3 L M h μ s O 3 O 3 O 3
Here, one can see that except ε ( t ) , all the elements in Equation (38) are negative definite. In Equations (39) and (40), F ˜ r ( t ) is the remaining interactive force between the hose and its surroundings after being compensated by the adaptive element ς ^ ( t ) in the control law, and Δ X ˙ is the difference between the head’s speed and the average speed of the hose, due to its flexibility. Then, by adopting Young’s inequality for ε ( t ) in (39), the differential V ˙ ( t ) was obtained as follows:
V ˙ ( t ) z ( t ) T P 1 z ( t ) m ˙ h ( t ) M h m p + ρ A h L 1 X ˙ ( L , t ) T X ˙ ( L , t ) m p + ρ A h 2 L M h 2 g 0 L 1 μ X ˙ ( s , t ) + Y ˙ ( s , t ) d s + σ , P 1 = d i a g κ 1 κ s 1 2 I , κ 3 1 2 P 2 T P 2 1 2 + m ˙ h ( t ) L 2 M h I , κ u 1 2 I , λ μ s 2 + L 2 L , σ = max 1 2 β ( t ) ω b ( t ) 2 2 + 1 2 d 2 ( t ) 2 2 + 1 2 Δ u X ( t ) Δ u φ ( t ) 2 2 + 1 2 d s ( t ) 2 2 + 1 2 d 2 ( t ) 2 2 + σ Δ X ˙ ( t )
By choosing the control matrices such that P 1 is positive definite, we thereby completed the proof of Lemma 2. □
Lemma 1 and 2 imply the uniform ultimate boundedness of the control system. From the right inequality of (31), we have:
z ( t ) 2 2 ε V 1 ( t ) 4 λ max Q 2 ε + V s 2 ( t ) + V ς ( t ) + Λ ( t )
From the left inequality of (31):
λ min Q 1 z ( t ) 2 2 + V s 2 ( t ) + V ς ( t ) + Λ ( t ) 4 λ max Q 2 ε + V s 2 ( t ) + V ς ( t ) + Λ ( t )
Then, the ultimate bound b with respect to z ( t ) can be taken as (44) regardless of either the hose’s curvature or the head’s fly height.
b = λ max Q 2 λ min Q 1 2 ε
This completes the proof for the Theorem. □

5. Results and Discussion

5.1. Comparative Simulation Studies

Based on the CAD model of the designed system, its specifications were determined and are shown in Table 1. In order to validate the feasibility of the proposed system, the mathematical models presented in Section 2 and the control law introduced in Section 4.1 were utilized in the simulations. Moreover, a linear state-feedback controller in [41] and a cascade sliding mode controller (SMC) from [29] were employed as comparative controllers. They represent two different approaches to controlling a water-powered aerial system with a flexible hose. For the state-feedback controller, the stability conditions of the hose were considered in terms of the saturated water flow rate. This saturation along with the limitations of the valve position was accounted for by optimally computing the control matrices. In contrast, no knowledge of the hose’s dynamics is required for the design of the cascade SMC. The water-conveying hose was only a lumped disturbance and a mass increasing proportionally with the system attitude. The controller was designed to be robust when facing these perturbations. Formulations of the state-feedback and cascade SMC laws are given in [29] and [41], respectively. Additionally, the parameters of the proposed control law and these two comparative controls are given in Table 2.
The simulations were conducted using Mathworks Simulink and Simscape Multibody. However, it was difficult to solve the partial differential equation (PDE) of the hose. Instead, the PDE (9) with the boundary conditions (10) was approximated using the finite difference method [44], resulting in a series of ordinary differential equations (ODEs) as in [41]. Two types of references were considered: step and sinusoidal. The former mimicked a point-to-point tracking scenario, where the system was required to approach several fire sites and complete firefighting tasks. The latter simulated continuous motion, such as following a desired trajectory, avoiding obstacles, or extinguishing fires in a wide area. For each reference, the system’s maneuverability was validated for each controller in two cases, without and with firefighting actions, resulting in a total of four scenarios. The simulation results are shown in the following sections and are also visually demonstrated in Supplementary Video S1.

5.2. Step-Type Reference Tracking

In this scenario, the system first takes off to a desired attitude of 3 [m] and maintains that height. It then tracks the desired setpoints in the corresponding horizontal plane and lands back on the ground at the end. Hence, the references of translational motions are in the form of step signals, and the setpoints are easily seen in Figure 5a. Moreover, a third-order filter was added for each signal to reflect the desired movements in practical operation. Thus, the resulting reference with respect to time is as in Figure 5b. During the process, the yaw angle of the head assembly also followed a sinusoidal trajectory, simulating the movement of widely spraying water to cover the flame. The resulting motions of the three control systems, including the head trajectory and postures of the head assembly and the hose, are shown in three-dimensional geometry in Figure 5a. The motion of the head assembly in every direction over time is given in Figure 5b. Moreover, their control inputs are subsequently illustrated in Figure 5c. In Figure 5d, every second, the posture of the hose is drawn once and then lined up side-by-side along the time axis. In this figure, the Y*-axis represents the Y-coordinate being spread out with time. As shown, all the considered control systems can reach the required positions, but the two comparative ones exhibit several visible deviations from the reference. In detail, one can see that the cascade SMC system can track the references in real time. It gives a good tracking performance in the Z-direction, slight overshoot in the ψ-direction, and smooth transient responses in the horizontal translations. However, the controller must make allotted corrections in its control inputs, especially those for flow-regulating valves b and c as seen in Figure 5c. Even so, the resulting motions in the horizontal plane are poor. This is because the opening and closing of the valve lead to the fluctuation of the water flow and generate additional unwanted disturbance acting on the head assembly. As shown in Figure 5d, the more even distribution of the hose’s postures, the better its stability. In this sense, one can see the fluctuation of the hose due to this control action. On the other hand, the state-feedback control requires less control effort and provides smoother motions over time, but the motions tend to fall behind the desired trajectory. The horizontal motions are still not robust as their steady-state error can easily increase under the influence of other factors, although it is better than the results from the cascade SMC. For example, one can see in Figure 5b, between 20 [s] and 30 [s], while tracking the reference in the X-direction, that the cascade SMC and state-feedback controls cannot maintain their position in the Y-direction. Meanwhile, the proposed control system has good and real-time tracking performance in all directions during the whole simulation. The proposed system compensates effectively for the influence of couplings and disturbances so that the resulting translation and heading angles closely follow their corresponding references. It is also worth noting that due to the underactuated nature, the translations cannot be achieved without controlling the roll and pitch of the system. Hence, the proposed control system adjusts more in these two rotations compared to the comparative systems, but by doing so, it retains superior tracking performance in the horizontal plane. In terms of the control inputs, the proposed controller generates a large correction in the initial state, especially for the control signal of the valves. This is because of the low water flow rate in this state resulting in a small thrust force that requires the head assembly to lean more to correct its initial position and heading direction. After that, the adjustment of these inputs is still slightly larger than that of the state-feedback control (to achieve better motion performance) but significantly smaller than that of the cascade SMC. This results in better stability of the hose with the proposed control and the state-feedback control, in comparison to the cascade SMC.
In the second scenario of this reference type, the system’s performance was investigated in its application, i.e., the fire hose nozzle sprays water to perform firefighting actions. The obtained results are shown side-by-side in Figure 6a,d. The fire extinguishing flow has a rate of 1 [kg/s], being initiated at 15 [s], then lasting 10 [s] every 10 [s] interval. Not only does this flow generate an additional thrust acting on the head assembly, but it also changes the water supply leading to a fluctuation of the hose. The latter adds disturbances to the head and may even lead to flutter instability of the hose. The performance of the cascade SMC system is found to be unstable due to these disturbances. Meanwhile, the design of the two others takes the hose stability into account and hence preserves the stability of the whole system even though the fluctuation of the hose is still larger than in the first scenario, due to the change in the water supply for the firefighting action as mentioned. Moreover, the state-feedback controller can preserve the system stability but cannot reject disturbances effectively. Therefore, the resulting motions, especially the translations, are significantly affected. In contrast, the proposed control system retains good performance, with slight effects of disturbances visible in the X- and Y-directions.

5.3. Sinusoidal Trajectory Following

In this case, the head assembly is required to follow a continuous trajectory described by:
X d = 2 sin ( 0.15 t ) , Y d = 3 sin ( 0.3 t ) , Z d = 2 + 0.3 sin ( 0.2 t ) + 0.03 t
The required heading angle is in the form of a rectangular pulse signal, initialized at 11 [s], with an amplitude of π / 6 [rad], lasting 10 [s] every 10 [s] interval. No fire extinguishing flow was considered in the first scenario with this reference. Figure 7a,d depict the simulation results in the same manner as the previous section. Similar results to those observed in the previous trajectory type can also be seen in this case. That is, the motions of the state-feedback control system fall behind the desired trajectories; the cascade SMC tracks the reference in real time but contains allotted deviation, and the proposed system remains the best solution. In the initial stage, all the control systems struggle to track their reference, due to the lack of water supply. Furthermore, all control systems experience challenges in tracking the translational motion smoothly, as the thrust force required to track the varying altitude changes. As seen in Equation (17), all acting torques are powered by this thrust. Hence, the flow-regulating valves must constantly adjust to retain the desired torques. In addition, the complex trajectory also leads to irregular shapes of the hose, such as twisting and releasing at around 45 [s] of the simulation. This obviously affects the motion of the head assembly and the stability of the hose, as visible in Figure 7a,b,d, and Supplementary Video S1. On the other note, the performance of the cascade SMC in the yaw motion is inferior to other controllers, despite being the second best controller in the previous reference. The cascade SMC quickly reaches the desired heading angle but cannot maintain a small steady-state error. In contrast, the state-feedback control system has a slow transient response but a small steady-state error. Meanwhile, the proposed system achieves excellent performance in the yaw direction, both transient and steady-state responses, as shown at the bottom right of Figure 7b.
In a similar intention to the second scenario with the step-type reference, the impact of firefighting actions was evaluated for this trajectory. The fire extinguishing flow was identical to the previous case, and its impact is visualized in the system performance in Figure 8. Accordingly, the cascade SMC is most affected followed by the state-feedback control and the proposed control that most effectively compensates for this firefighting action. Additionally, in Table 3, two metrics are presented to evaluate the performance of the three control systems in the four simulations: the average root mean square errors (RMSEs) of the head motions and the integral of the square of the control input values (ISVs). The proposed control system shows superior positioning performance, as evidenced by a significant reduction in RMSE of 87% and 90% when compared to the cascade SMC and state-feedback control systems, respectively. However, for orientation, the average RMSE is not as good because the pitch and roll of the head must be alternated in order to position the system. When considering the yaw rotation only, the proposed system has a significantly better tracking performance, with an RMSE of 0.007 compared to 0.0162 and 0.124 in the cascade SMC and state-feedback systems, respectively. The ISV of the flow rate indicates that the proposed system consumes more water than the state-feedback control system but less than the SMC one. Fortunately, it is not a concern in water-based environments, which are the target application scenario of the proposed firefighting system. In the proposed system, the flow-regulating valves require more control actions than those with the state-feedback control but are significantly less intense than the cascade SMC system. From the results of the four simulation studies, the designed water-powered aerial system shows its capability, and the proposed control law demonstrates its superiority compared to the two control systems.

6. Conclusions

In this study, a novel water-powered aerial system with flow-regulating actuators and an effective control law was developed and validated. The system was inspired by hydroflight equipment and is capable of autonomous firefighting missions, particularly in or at bodies of water areas. The system comprises a water pump, a flexible hose, and a head assembly that can perform motions in three-dimensional space. The mathematical model was established to describe the dynamics of the whole system, and several dynamical characteristics were deduced from it. The proposed control law incorporates these characteristics into its design to result in the uniform ultimate boundedness of the system’s motion. Its effectiveness and robustness were evaluated via simulation studies, in which its motion performance and stability were compared to those controlled by different control approaches. Accordingly, the proposed system demonstrated its maneuverability and superiority. The system prototype will be fabricated, and experimental studies will be conducted for further evaluation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/drones7030162/s1, Video S1: Supplementary Video S1.

Author Contributions

Conceptualization, Y.-B.K. and T.H.; methodology, Y.-B.K. and T.H.; software, Y.-B.K.; validation, Y.-B.K. and T.H.; formal analysis, Y.-B.K.; investigation, T.H.; resources, Y.-B.K.; writing—original draft preparation, T.H.; writing—review and editing, Y.-B.K. and T.H.; visualization, T.H.; supervision, Y.-B.K.; funding acquisition, Y.-B.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. 2022R1A2C1003486).

Data Availability Statement

Not applicable.

Acknowledgments

This work was also supported by the National Research Foundation (NRF), South Korea under Project BK21 FOUR (Smart Convergence and Application Education Research Center).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Robinson, B. Water Propelled Flying Board. U.S. Patent US9145206B1, 29 September 2015. [Google Scholar]
  2. Li, R. Personal Propulsion Device. U.S. Patent US7735772B2, 21 August 2007. [Google Scholar]
  3. Water—Zapata. Available online: https://www.zapata.com/en/water/ (accessed on 21 September 2022).
  4. Zapata, F. Maneuvering and Stability Control System for Jet-Pack. U.S. Patent US20140103165A1, 17 April 2014. [Google Scholar]
  5. Jetovator—Rise the Hose. Available online: https://www.jetovator.com/ (accessed on 21 September 2022).
  6. Qin, H.; Cui, J.Q.; Li, J.; Bi, Y.; Lan, M.; Shan, M.; Liu, W.; Wang, K.; Lin, F.; Zhang, Y.F.; et al. Design and implementation of an unmanned aerial vehicle for autonomous firefighting missions. In Proceedings of the 2016 12th IEEE International Conference on Control and Automation (ICCA), Kathmandu, Nepal, 1–3 June 2016; IEEE: Piscataway, NJ, USA, 2016; pp. 62–67. [Google Scholar]
  7. Spurny, V.; Pritzl, V.; Walter, V.; Petrlik, M.; Baca, T.; Stepan, P.; Zaitlik, D.; Saska, M. Autonomous Firefighting Inside Buildings by an Unmanned Aerial Vehicle. IEEE Access 2021, 9, 15872–15890. [Google Scholar] [CrossRef]
  8. EHang|Smart City Management—EHang 216F. Available online: https://www.ehang.com/ehang216f/ (accessed on 21 September 2022).
  9. High-Rise Firefighting with Dry Powder Fire Extinguishing—Guofei UAV. Available online: http://www.guofei-uav.com/index.php?c=article&a=type&tid=43 (accessed on 21 September 2022).
  10. Chaikalis, D.; Tzes, A.; Khorrami, F. Aerial Worker for Skyscraper Fire Fighting using a Water-Jetpack Inspired Approach. In Proceedings of the 2020 International Conference on Unmanned Aircraft Systems (ICUAS), Athens, Greece, 1–4 September 2020; pp. 1392–1397. [Google Scholar]
  11. Chaikalis, D.; Evangeliou, N.; Tzes, A.; Khorrami, F. Design, Modelling, Localization, and Control for Fire-Fighting Aerial Vehicles. In Proceedings of the 2022 30th Mediterranean Conference on Control and Automation (MED), Vouliagmeni, Greece, 28 June–1 July 2022; pp. 432–437. [Google Scholar]
  12. Lee, S.; Ng, W.; Liu, J.; Wong, S.; Srigrarom, S.; Foong, S. Flow-Induced Force Modeling and Active Compensation for a Fluid-Tethered Multirotor Aerial Craft during Pressurised Jetting. Drones 2022, 6, 88. [Google Scholar] [CrossRef]
  13. Viegas, C.; Chehreh, B.; Andrade, J.; Lourenço, J. Tethered UAV with Combined Multi-rotor and Water Jet Propulsion for Forest Fire Fighting. J. Intell. Robot. Syst. 2022, 104, 21. [Google Scholar] [CrossRef]
  14. Ausonio, E.; Bagnerini, P.; Ghio, M. Drone Swarms in Fire Suppression Activities: A Conceptual Framework. Drones 2021, 5, 17. [Google Scholar] [CrossRef]
  15. Harikumar, K.; Senthilnath, J.; Sundaram, S. Multi-UAV Oxyrrhis Marina-Inspired Search and Dynamic Formation Control for Forest Firefighting. IEEE Trans. Autom. Sci. Eng. 2019, 16, 863–873. [Google Scholar] [CrossRef]
  16. Dubai Firefighters Use Jetpacks to Aid High-Speed Response—Euronews. Available online: https://www.youtube.com/watch?v=H1VIxUTzb1I (accessed on 21 September 2022).
  17. Dubai Media Office. Available online: https://twitter.com/DXBMediaOffice/status/822797935176024065 (accessed on 20 February 2023).
  18. Ando, H.; Ambe, Y.; Ishii, A.; Konyo, M.; Tadakuma, K.; Maruyama, S.; Tadokoro, S. Aerial Hose Type Robot by Water Jet for Fire Fighting. IEEE Robot. Autom. Lett. 2018, 3, 1128–1135. [Google Scholar] [CrossRef]
  19. Ando, H.; Ambe, Y.; Yamaguchi, T.; Yamauchi, Y.; Konyo, M.; Tadakuma, K.; Maruyama, S.; Tadokoro, S. Fire extinguishment using a 4 m long flying-hose-type robot with multiple water-jet nozzles. Adv. Robot. 2020, 34, 700–714. [Google Scholar] [CrossRef]
  20. Ando, H.; Ambe, Y.; Yamaguchi, T.; Konyo, M.; Tadakuma, K.; Maruyama, S.; Tadokoro, S. Fire Fighting Tactics with Aerial Hose-type Robot “Dragon Firefighter”. In Proceedings of the 2019 IEEE International Conference on Advanced Robotics and its Social Impacts (ARSO), Beijing, China, 31 October–2 November 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 291–297. [Google Scholar]
  21. Ando, H.; Ambe, Y.; Yamaguchi, T.; Konyo, M.; Tadakuma, K.; Maruyama, S.; Tadokoro, S. Activation of Hose by Multiple Water-Jet Injection Nozzles for Fire Fighting. Proc. JSME Annu. Conf. Robot. Mechatron. (Robomec) 2018, 2018, 2A1-M02. [Google Scholar] [CrossRef]
  22. Dinh, C.-T.; Huynh, T.; Kim, Y.-B. LQI Control System Design with GA Approach for Flying-Type Firefighting Robot Using Waterpower and Weight-Shifting Mechanism. Appl. Sci. 2022, 12, 9334. [Google Scholar] [CrossRef]
  23. Liu, X.; Zhou, H. Unmanned Water-Powered Aerial Vehicles: Theory and Experiments. IEEE Access 2019, 7, 15349–15356. [Google Scholar] [CrossRef]
  24. Ambe, Y.; Yamauchi, Y.; Konyo, M.; Tadakuma, K.; Tadokoro, S. Stabilized Controller for Jet Actuated Cantilevered Pipe Using Damping Effect of an Internal Flowing Fluid. IEEE Access 2022, 10, 5238–5249. [Google Scholar] [CrossRef]
  25. Yamaguchi, T.; Ambe, Y.; Ando, H.; Konyo, M.; Tadakuma, K.; Maruyama, S.; Tadokoro, S. A Mechanical Approach to Suppress the Oscillation of a Long Continuum Robot Flying With Water Jets. IEEE Robot. Autom. Lett. 2019, 4, 4346–4353. [Google Scholar] [CrossRef]
  26. Ambe, Y.; Ando, H.; Konyo, M.; Tadokoro, S. Design of a simple controller for the flying hose type robot with water jet. In The Proceedings of JSME annual Conference on Robotics and Mechatronics (Robomec), Fukushima, Japan, 10 – 13 May 2017; The Japan Society of Mechanical Engineers: Shinjuku-ku, Tokyo, Japan, 2017; p. 1P1-P02. [Google Scholar]
  27. Yamauchi, Y.; Ambe, Y.; Konyo, M.; Tadakuma, K.; Tadokoro, S. Realizing Large Shape Deformations of a Flying Continuum Robot with a Passive Rotating Nozzle Unit That Enlarges Jet Directions in Three-Dimensional Space. IEEE Access 2022, 10, 37646–37657. [Google Scholar] [CrossRef]
  28. Yamauchi, Y.; Ambe, Y.; Konyo, M.; Tadakuma, K.; Tadokoro, S. Passive Orientation Control of Nozzle Unit with Multiple Water Jets to Expand the Net Force Direction Range for Aerial Hose-Type Robots. IEEE Robot. Autom. Lett. 2021, 6, 5634–5641. [Google Scholar] [CrossRef]
  29. Lee, D.-H.; Huynh, T.; Kim, Y.-B.; Soumayya, C. Motion Control System Design for a Flying-Type Firefighting System with Water Jet Actuators. Actuators 2021, 10, 275. [Google Scholar] [CrossRef]
  30. Samadikhoshkho, Z.; Lipsett, M. Decoupled Control Design of Aerial Manipulation Systems for Vegetation Sampling Application. Drones 2023, 7, 110. [Google Scholar] [CrossRef]
  31. Shi, C.; Wang, K.; Yu, Y. Expandable Fully Actuated Aerial Vehicle Assembly: Geometric Control Adapted from an Existing Flight Controller and Real-World Prototype Implementation. Drones 2022, 6, 272. [Google Scholar] [CrossRef]
  32. Yu, Y.; Shi, C.; Shan, D.; Lippiello, V.; Yang, Y. A hierarchical control scheme for multiple aerial vehicle transportation systems with uncertainties and state/input constraints. Appl. Math. Model. 2022, 109, 651–678. [Google Scholar] [CrossRef]
  33. Chen, W.; Dai, H.; Wang, L. Three-dimensional dynamical model for cantilevered pipes conveying fluid under large deformation. J. Fluids Struct. 2021, 105, 103329. [Google Scholar] [CrossRef]
  34. Wadham-Gagnon, M.; Paı¨doussis, M.P.; Semler, C. Dynamics of cantilevered pipes conveying fluid. Part 1: Nonlinear equations of three-dimensional motion. J. Fluids Struct. 2007, 23, 545–567. [Google Scholar] [CrossRef]
  35. Paidoussis, M.P. Fluid-Structure Interactions: Slender Structures and Axial Flow; Academic Press: Cambridge, MA, USA, 1998; Volume 1, ISBN 008053175X. [Google Scholar]
  36. He, W.; He, X.; Ge, S.S. Vibration Control of Flexible Marine Riser Systems with Input Saturation. IEEE/ASME Trans. Mechatron. 2016, 21, 254–265. [Google Scholar] [CrossRef]
  37. He, W.; Zhang, S.; Ge, S.S. Boundary Control of a Flexible Riser With the Application to Marine Installation. IEEE Trans. Ind. Electron. 2013, 60, 5802–5810. [Google Scholar] [CrossRef]
  38. He, W.; He, X.; Sam Ge, S. Modeling and Vibration Control of a Coupled Vessel-Mooring-Riser System. IEEE/ASME Trans. Mechatron. 2015, 20, 2832–2840. [Google Scholar] [CrossRef]
  39. He, W.; Sun, C.; Ge, S.S. Top tension control of a flexible marine riser by using integral-barrier lyapunov function. IEEE/ASME Trans. Mechatron. 2015, 20, 497–505. [Google Scholar] [CrossRef]
  40. Ge, S.S.; He, W.; How, B.V.E.; Choo, Y.S. Boundary Control of a Coupled Nonlinear Flexible Marine Riser. IEEE Trans. Control. Syst. Technol. 2010, 18, 1080–1091. [Google Scholar] [CrossRef]
  41. Huynh, T.; Lee, D.-H.; Kim, Y.-B. Study on Actuator Performance Evaluation of Aerial Water-powered System for Firefighting Applications. Appl. Sci. 2023, 13, 1965. [Google Scholar] [CrossRef]
  42. Li, R. Personal propulsion devices with improved balance. U.S. Patent US9849980B2, 18 September 2014. [Google Scholar]
  43. Sahoo, P.K. Some results related to the integral mean value theorem. Int. J. Math. Educ. Sci. Technol. 2007, 38, 818–822. [Google Scholar] [CrossRef]
  44. Leader, J.J. Numerical Analysis and Scientific Computation; CRC Press: Boca Raton, FL, USA, 2022; ISBN 1000540391. [Google Scholar]
Figure 1. Putting out a boat fire.
Figure 1. Putting out a boat fire.
Drones 07 00162 g001
Figure 2. The proposed aerial water-powered firefighting system with flow rate control: (a) Proposed system concept drawing and (b) Diagram of water flow in the head part.
Figure 2. The proposed aerial water-powered firefighting system with flow rate control: (a) Proposed system concept drawing and (b) Diagram of water flow in the head part.
Drones 07 00162 g002
Figure 3. Coordinate definition and forces analysis.
Figure 3. Coordinate definition and forces analysis.
Drones 07 00162 g003
Figure 4. Structure of the proposed control system.
Figure 4. Structure of the proposed control system.
Drones 07 00162 g004
Figure 5. Motion performance in the case of step-type reference tracking. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Figure 5. Motion performance in the case of step-type reference tracking. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Drones 07 00162 g005
Figure 6. Motion performance in the case of step-type reference tracking with firefighting action. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Figure 6. Motion performance in the case of step-type reference tracking with firefighting action. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Drones 07 00162 g006aDrones 07 00162 g006b
Figure 7. Motion performance in the case of sinusoidal trajectory following. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Figure 7. Motion performance in the case of sinusoidal trajectory following. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Drones 07 00162 g007aDrones 07 00162 g007b
Figure 8. Motion performance in the case of sinusoidal trajectory following. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Figure 8. Motion performance in the case of sinusoidal trajectory following. (a) Tracking performance in 3-dimensional geometry; (b) Time responses of translations and rotations; (c) Three-dimensional partial motion of the water-conveying hose; and (d) Control inputs.
Drones 07 00162 g008aDrones 07 00162 g008b
Table 1. Pilot model specifications and finite difference method parameters.
Table 1. Pilot model specifications and finite difference method parameters.
ParameterNotationValueUnit
Mass M h 3.10 kg
Inertia matrix J h 0.035 0.001 0.003 0.001 0.034 0.001 0.003 0.001 0.034 kg . m 2
Dimension l w h T 0.12 0.23 0.03 T m
Cross area of water inlet A h 0.025 2 π m 2
Cross area of thrust nozzle outlet A n 0.002 2 π m 2
Inclined angle of thrust nozzle δ π / 6 rad
Saturation of water flow rate m ˙ h , min , m ˙ h , max 0 , 7.5 kg / s
Saturation of valve position u i , min , u i , max 1 , 1
Cross area of flexible hose A h 0.025 2 π m 2
Stiffness of flexible hose 2.866 N . m / rad
Damping coefficient of hose 0.143 N . m . s / rad
Water density ρ 1000 kg / m 3
Cross area of fire hose nozzle A l 0.01 2 π m 2
Dimension of fire hose nozzle l w h T , δ 0.18 0 0.07 T , π / 18 m, rad
Table 2. Control system parameters.
Table 2. Control system parameters.
Cascade sliding mode control n X 1 = d i a g { 1.5   1.5   2.5 } , n X 2 = d i a g { 0.5   0.5   0.5 } , K X = d i a g { 15   20   120 } , n Φ = d i a g { 1.5   1.5   1.5 } , K Φ = d i a g { 25   25   75 }
Linear state-
feedback control
K f = 0 0 2.544 0 0 0 0 0 0 0 0 0.628 1.659 0.64 0 1.492 3.923 0 1.659 0.64 0 1.492 3.923 0 0 0 0.755 0 0 0 0 0 0 0 0 0.064 0.711 0.271 0 0.056 0.156 0 0.711 0.271 0 0.056 0.156 0 , K g = 0 0 2.113 0 0 0 0 0.48 1.264 0.488 0 0 1.264 0.488 0 0
Proposed control K 1 = d i a g { 21   32   75   50   50   110 } ; K 3 = d i a g { 3.6   4   5   10   20   36 } ; K s = d i a g { 1   4   4   1   1   20 } ; K u = d i a g { 4   4   10   4   4   10 } ; K ς = 50 ; ω n = 50 I 6 ; I 6
Table 3. Comparison of average performance indexes.
Table 3. Comparison of average performance indexes.
Control SystemRMSEISV
PositionOrientationFlow RateValve Operation
Cascade SMC control0.3280.141425.7250.505
State-feedback control0.4280.146414.9790.091
Proposed control0.0440.107422.8790.141
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huynh, T.; Kim, Y.-B. Motion Control System Design for a Novel Water-Powered Aerial System for Firefighting with Flow-Regulating Actuators. Drones 2023, 7, 162. https://doi.org/10.3390/drones7030162

AMA Style

Huynh T, Kim Y-B. Motion Control System Design for a Novel Water-Powered Aerial System for Firefighting with Flow-Regulating Actuators. Drones. 2023; 7(3):162. https://doi.org/10.3390/drones7030162

Chicago/Turabian Style

Huynh, Thinh, and Young-Bok Kim. 2023. "Motion Control System Design for a Novel Water-Powered Aerial System for Firefighting with Flow-Regulating Actuators" Drones 7, no. 3: 162. https://doi.org/10.3390/drones7030162

APA Style

Huynh, T., & Kim, Y. -B. (2023). Motion Control System Design for a Novel Water-Powered Aerial System for Firefighting with Flow-Regulating Actuators. Drones, 7(3), 162. https://doi.org/10.3390/drones7030162

Article Metrics

Back to TopTop