Next Article in Journal
A Novel Approach for Micro-Antenna Fabrication on ZrO2 Substrate Assisted by Laser Printing for Smart Implants
Previous Article in Journal
DPS Clustering: New Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LQI Control System Design with GA Approach for Flying-Type Firefighting Robot Using Waterpower and Weight-Shifting Mechanism

1
Department of Intelligent Robot Engineering, Pukyong National University, Nam-gu, Busan 48513, Korea
2
Department of Chassis and Body, Ho Chi Minh City University of Technology and Education, Thu Duc, Ho Chi Minh City 700000, Vietnam
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(18), 9334; https://doi.org/10.3390/app12189334
Submission received: 10 August 2022 / Revised: 7 September 2022 / Accepted: 14 September 2022 / Published: 17 September 2022
(This article belongs to the Section Robotics and Automation)

Abstract

:
This study proposes a flying robot using waterpower and a novel weight-shifting mechanism, whose purpose is to be applied in firefighting tasks in water areas that are difficult to access and are suppressed by conventional firefighting methods. The sufficient amount of water in the area is used for propelling the robot, as well as fire suppression activity. A weight-shifting mechanism governs the weight distribution of the robot head in order to perform the robot motions. In this paper, the system’s dynamical characteristics are analyzed in detail through mathematical models. A linear-quadratic integrator (LQI) is designed for controlling the system motion. Additionally, the LQI is tuned with the genetic algorithm (GA) so that both system performance and robustness are optimally preserved. Simulation studies are carried out, in which the proposed system is compared with a cascade proportional-integral-derivative (PID) control system. The results validate the feasibility of the design and show the superiority of the proposed control system in motion performances. Moreover, the LQI-GA consumes 2.28% less water and uses 83.85% less kinetic energy of the actuator than the PID control system.

1. Introduction

Fire accidents are unpredictable and can spread quickly. They often cause significant damage to people and property. Therefore, the development of improved and effective fire suppression systems always attracts much interest. Automatic fire-extinguishing and sprinkler systems have been widely used and have proven their efficiency in indoor fire accidents. Nevertheless, in many cases, firetrucks and firefighters need to directly battle flames to preserve property and the lives of humans or animals. The firefighters always have to contend with multiple difficulties—for example, carrying heavy hoses and tools in and around heavy pieces of unsafe structures. The brutal conditions also put the firefighters at risk of loss of life or injury. Moreover, in some tricky conditions, such as fire accidents happening on a narrow bridge, a crowded harbor, or on the sea, it is extremely difficult for the fire trucks and firefighters to approach the location to do their mission.
Over the past decade, unmanned vehicles or remotely controlled robots have been deployed in fire accidents, not only to assist but sometimes to replace the manpower in firefighting missions. Several intelligent firefighting systems with different shapes, features, and extinguishing agents have been developed for different situations. Generally, they can be classified into two main groups according to the method of approaching the flames: terrestrial systems and aerial systems.
In chronological order, the terrestrial systems were the ones that were first studied and used in fire accidents in residential areas, industrial fields, and wildfires [1,2,3]. They can be ground vehicles or wheeled, legged, and even humanoid robots. Commercial firefighting robots [4,5,6,7] are equipped with modern technologies that greatly increase their durability, capability, and feasibility. For instance, the first commercialized firefighting robots, Thermite RS1 and RS3 (Textron, Providence, RI, USA) [4], the AirCore TAF35 (MAGIRUS, Ulm, Germany) [5], and the Shark-robotic Colossus (Shark Robotics, La Rochelle, France) [6], are wheeled robots that carry water hoses and spaying nozzles to attack the fire. They also provide reconnaissance and situational awareness in high-risk areas through camera systems, so they can operate remotely. The power engines and continuous track configurations allow them to carry a huge amount of water, overcome obstacles, and traverse difficult terrain. This enables them to battle the fire from inside or at a closer distance, which is too dangerous for a firefighter. In addition, the Firefighting Robot System (Mitsubishi Heavy Industries, Ltd., Chiyoda City, Tokyo, Japan) [7] includes a firefighting robot equipped with onboard sensors, such as RTK-GPS (SoftBank Group, Minato City, Tokyo, Japan), LiDAR (Velodyne Lidar, San Jose, CA, USA), IMU (Imugene Limited, New South Wales, Australia), and odometry, and a hose extension accessory. The completed system has the ability to self-drive directly to the scene of a fire and perform the firefighting mission effectively. However, the large dimensions and heavy weight of the wheeled robots make them hard to operate in narrow environments or high locations.
Animal-inspired robots—for instance, snake firefighting robots [8] and bug firefighting robots [9]—have been considered in these situations. The snake robot designed by Liljeback et al. [8] works as a self-propeller fire hose that can crawl into a burning narrowed space and extinguish a fire. Dumiak [9] proposed the idea of bug-inspired legged robots going around in nature to detect the location of fire ignition and suppressing it to prevent large-scale wildfires. In addition, Kim et al. [10] proposed deploying the humanoid robot SAFFiR (Robotiq, Lévis, QC, Canada) for firefighting in naval vessels. The animal-inspired and humanoid robots are flexible, so they suit the missions in narrow environments or complex structures. One drawback that would impact their use would be low visibility due to smoke and extremely harsh environments. Additionally, with high-rise buildings, traffic congestion, or water scenarios, these robots are not practical to deploy. Moreover, the operation range of terrestrial firefighting robots is generally limited.
On the other hand, aerial vehicles, especially multi-propeller drones, can easily access these locations, so they can inject fire-extinguishing agents to suppress the flames. The world’s first flying-type firefighting drones, Ehang 216-F (Ehang, Guangzhou, China) [11] and the ZHUN Walkera (Guangzhou Walkera Technology Co., Ltd., Guangzhou, China) [12], used firefighting foam and dry powder, respectively, for the task. They are specially designed for suppressing fires in high-rise buildings. However, a drone can only carry a certain amount of the firefighting agent, so it is not possible to work continuously and be ineffective in serious fire accidents. This can be overcome by continuously conveying the water from a ground hydrant or pump to the drone’s nozzle via a flexible fire hose. Alternate systems such as Aerones (Riga, Latvia) [13] and Goufei Fire-fighting UAV (Goufei Aviation, Chongquing, China) [14] were designed to overcome this issue. One drawback is that the higher the drone flies, the heavier it becomes. The higher power demanded of the propeller drone and its battery’s capacity constraint meant that the fire accident location affects both the flight range and the time of the drone. One other thing to consider is that the airflow generated by the propellers can also spread the flame to neighboring areas.
More recently, flying-type firefighting robots that use water thrust force have been considered, especially for a fire that happens in water areas. The approach was inspired by water-powered flying devices such as flying boards [15,16] and water jetpacks [17,18,19]. In [20,21,22,23], the Dragon Firefighter system generates thrust by speeding up the water flow through a nozzle assembly and uses a servos actuator system to adjust the angle of nozzles, thus adjusting the force directions. This not only lifts the hose and the system, but the water jet actuator also directs the water flow to suppress the flame. Unfortunately, rotational motions are very limited with the system due to its structure. A later flying-type firefighting system in [24] uses four nozzles for maneuvering and a separated sprinkler to suppress the fire. The movement of the proposed system is carried out by changing the flowrate of each nozzle. The development of this system is still in the initial stage. The placement of the actuation system and the location of the water-conveying hose have not been finalized yet.
Furthermore, these abovementioned water-powered firefighting systems are generally in a stage of development and completion. Thus, the conceptual designs and mechanical designs are investigated preferentially. The motion control methods, therefore, are fairly simple. In particular, a proportional controller with speed feedback is implemented in a Dragon Firefighter [20,21,22,23]. The simulations and experiments are adopted, and it is indicated that this controller is not sophisticated and cannot control the system precisely [20]. The investigation of the latter mechanical approach helps it works better with the same controller [21,22,23]. However, a more advanced control technique is also initially implemented to evaluate the feasibility of water-powered flying-type firefighting [24]. A sliding mode control technique is proposed to govern the system motions. The simulation tests indicated the robustness in tracking performance. However, the control effort itself is predicted to bring a challenge for practical applications due to the controller’s actions.
Therefore, a flying-type firefighting robot that is water-powered, automatically controlled, and has a novel weight-shifting mechanism for maneuvering is proposed in this study. The proposed robot uses the water thrust as the propulsion force, similar to the abovementioned system. The robot is aimed to be applied in firefighting tasks in water areas that have an unlimited amount of water but are difficult to access and are suppressed by conventional fire-extinguishing methods. The weight-shifting mechanism adjusts the weight distribution on the system, in a similar way to human movement on the fly board, in order to fly around. The target scenarios for the system application are visualized in Figure 1. In particular, a pump and a flexible hose convey water to generate water to the robot head. The head part generates thrust to take off and adjusts the actuator mechanism to approach the flame. The fire-extinguishing sprinkler sprays the water to suppress the fire remotely or automatically. The zoomed figure visualizes the robot head with its thrust nozzles, fire-extinguishing sprinkler, and weight-shifting mechanism.
In this paper, the system conceptual design is introduced. Mathematical models reveal the dynamical characteristics of the system. To demonstrate the maneuverability and feasibility of the proposed system, a linear-quadratic integrator (LQI) is designed such that the water flowrate and the weight-shifting mechanism are controlled to fly the system to the fire area. Moreover, an optimal approach based on the genetic algorithm (GA) tunes the controller’s gains. Thus, both motion performance and system robustness are preserved.
Accordingly, the contributions of this study in comparison to the relevant studies are as follows:
  • A novel conceptual design of the flying firefighting robot using waterpower and a weight-shifting mechanism is proposed.
  • The mathematical models of the system are formulated, and its characteristics are analyzed.
  • An LQI design with a novel GA approach is introduced.
  • Comparative simulation studies reveal the feasibility and efficiency of the proposed system.
The remainder of the paper is organized as follows. Section 2 introduces the system configuration and mathematical models. The LQI control system design with the GA approach is explained in Section 3. The simulation results are discussed in Section 4. Finally, the conclusions are drawn in Section 5.

2. System Configuration and Modeling

2.1. System Configuration

The flying firefighting robot flies thanks to the thrust forces generated from the water flow through the robot head part. In particular, a water pump placed on the ground or floating on the water, depending on the fire accident location, provides the water through a flexible hose to the inlet of the head part. The robot’s head distributes the water flow into four nozzle outlets, whose cross-sectional areas are much smaller than that of the inlet. The nozzles are fixed and pointed straight down for maximum thrust. The water flows are sped up by the nozzle and generate thrust forces. On the other side of the head, two weight blocks, WX and WY, are integrated with linear motors that allow them to move along the axes set in the longitudinal and transverse directions. When the weights are out of the central positions, the mass distribution of the head becomes unbalanced. It results in the head tilting and/or rolling, and subsequently, the robot moves transversely in the corresponding direction. In addition, the head also carries a water sprinkler for battling fire flames. The sprinkler is designed with an external nozzle ring, which atomizes water into a fine mist, which also helps to reduce heat from the flame to the devices on the robot’s head. In this conceptual design, a swivel ball joint is proposed as the connection between the water hose and the robot’s head to reduce the effect of the movement of the hose. The conceptual design of the robot head is drawn in Figure 2 below.

2.2. System Modeling

To analyze the motion characteristics of the system, the Earth-fixed reference frame OXYZ and body-fixed frame Obxbybzb are defined, as shown in Figure 2. The former has an X-axis pointing to the magnetic north and a Z-axis pointing up. The latter is placed at the center of the mass of the head, in which its xb-axis is along the longitudinal direction of the rigid body and points forward. The yb-axis follows the lateral direction of the robot’s head, and finally, the zb-axis points upward. Let ω x b ω y b ω z b T denote the angular velocity of the body-fixed frame orientation. The body axes rotate from the reference frame axes via a series of Euler angle rotations. They are roll, pitch, and yaw, respectively, and they are denoted by the vector ϕ θ ψ T . The kinematics of the designed robot in the reference frame are written in terms of the Euler angles, as follows:
X ˙ Y ˙ Z ˙ ϕ ˙ θ ˙ ψ ˙ = R x ˙ b y ˙ b z ˙ b ω x b ω y b ω z b T , R = R b E O 3 O 3 T b E
where O m × p denotes an m-by-p matrix of zeros, R b E is the rotational matrix from the body-fixed frame to the reference frame [25], and T b E is the transformation matrix between the angular velocity and the Euler angle rate [26], respectively.
For the robot dynamics, the free body diagram of the head part is illustrated in Figure 3. The hose conveying water and the actuating weights are isolated from the head part. The water hose is replaced by its mass being lifted M p at the center of the head and a reaction force F0 at the water inlet port. The actuating weights are represented by the reaction forces due to their motions and position. Additionally, the thrust force generated by the ith nozzle outlet is denoted by Fi (i = 1, 2, 3, 4). The masses of the actuating weights WX and WY are mwx and mwy, respectively. Meanwhile, the mass of the remaining head part with its containing water is M h . Additionally, when the firefighting robot gets close to the flame and the extinguishing nozzle activates, a reaction force Fs appears and acts on the system.
Considering the head part without the two actuating weights as a rigid body, thanks to the isolation of the free-body diagram, the dynamic equations for translational and rotational motions in the body-fixed frame are well-known as:
M h + M p v ˙ b + ω b × v b = F b J ω ˙ b + ω b × J ω b = τ b
where v b and ω b are the velocities of the translational and rotational motions in the body-fixed frame. The mass in the translational motions consists of the mass of the robot’s head M h and that of the water pipe being lifted M p . For the sake of simplicity, M p is assumed to linearly vary with the robot’s flight altitude via a correction factor k such that M p = k ρ A Z . ρ is the density of the water, and A is the internal cross-sectional area of the hose. ρ A Z is the mass of a water cylinder of area A and height Z. This assumption was also experimentally verified in the literature [19]. Moreover, the correction factor can be validated based on later experiments. For the rotational motions, J is the inertia matrix of the head, which is assumed to be diagonal to reflect the objective of the technical design of the fabricated model. That is, J = d i a g J x x , J y y , J z z . In Equation (2), M h + M p v ˙ b and J ω ˙ b are the differentiation of the system linear momentum and the moment-of-momentum, respectively. The cross-product terms ω b × v b and ω b × J ω b represent the Coriolis phenomenon.
On the right-hand side, the total forces and torques acting on the robot’s head are given by F b and τ b , respectively. They include the gravitation force, the water thrust (Ft, τt), actuating forces ( F w x , F w y ) and torques ( τ w x , τ w y ) from the weight-shifting mechanism, reaction forces from the hose F h and the water sprinkler F s , and other external disturbances ( F d , τ d ), such as the wind gust, ceiling, and ground effect. That is:
F b = F t F w x F w y + M h + M p R b E T 0 0 g T F h F s + F d τ b = τ t τ w x τ w y + τ d
The water thrust can be obtained by applying Newton’s second law on the water flowing through the head part. That is, the time rate of change of the linear momentum and moment-of-momentum of the water equals the summation of external forces and torques, respectively, acting on the water [24]. To obtain simple but sufficiently reliable dynamics of the fluid, assume that the water flow through the head is frictionless and steady [27]. The water thrust force and torque are written as follows:
F t = 0 0 1 ρ A m ˙ 0 2 + 1 ρ a i = 1 4 m ˙ i 2 + p p A + i = 1 4 m ˙ i ω b × r i τ t = i = 1 4 m ˙ i r i × ω b × r i
In which m ˙ 0 and m ˙ i are the mass flowrates at the water inlet port and the ith nozzle outlet, respectively. a is the cross-sectional of each nozzle outlet. The term p p A represents the water force acting at the water inlet due to the gauge pressure p p at this area. Furthermore, the vector ri determines the position of the ith nozzle outlet in the robot’s body frame.
On the other hand, the weight-shifting mechanism affects the robot’s head by the reaction forces and torques resulting from the motion of its two weight blocks. Equations of their motion are written as follows:
m w j v ˙ w j + ω b × v w j = F w j + m w j R b E T 0 0 g T , J w j ω ˙ b + ω b × J w j ω b = τ w j + m w j r w j × R b E T 0 0 g T , j = x , y
The translational motion of each weight consists of its displacement ( r w x and r w y ) observed from the robot’s head and the motion of the head itself. Thereby, the translational velocities of the weight WX and WY are given by:
v w x = r ˙ w x 0 0 T + v b + r w x × ω b , r w x = r w x 0 h w x T , v w y = 0 r ˙ w y 0 T + v b + r w y × ω b , r w y = 0 r w y h w y T
where r ˙ w x and r ˙ w y are their respective relative velocities, and h w x and h w y are their corresponding vertical distances from the coordinate origin.
On the order hand, the weights rotate with the rotation of the head part carrying them. Their moments of inertia taken at the body coordinate vary with their relative displacements. In other words:
J w x = m w x h w x 2 0 h w x r w x 0 h w x 2 + r w x 2 0 h w x r w x 0 r w x 2 , J w y = m w y h w y 2 + r w y 2 0 0 0 h w y 2 h w y r w y 0 h w y r w y r w y 2
Moreover, consider that four nozzles are perfectly symmetric, and because the head is kept stabilized in operation, neglect the nozzles’ elevation difference. The continuity equation indicates that the water flows through the nozzles are equal:
m ˙ 1 = m ˙ 2 = m ˙ 3 = m ˙ 4 = 1 4 m ˙ 0
The thrust force of the system is generated by these flows, so the different behaviors of the nozzle—for example, when the head tilts or rolls—are always able to be considered as input disturbances.
In summary, by substituting Ft and τt from Equation (4), Fwx, Fwy, τwx, and τwy from Equations (5)–(7), and m ˙ i from Equation (8), the equations of motion in all directions are obtained as follows:
x ¨ b = g sin θ + 1 a 1 h m ˙ 0 ω y b + d 1 , y ¨ b = g cos θ sin ϕ 1 a 1 h m ˙ 0 ω x b + d 2 , z ¨ b = 1 ρ a 1 1 A + 1 4 a m ˙ 0 2 g cos θ cos ϕ + d 3 , ω ˙ x b = 1 σ 1 m w x h x + m w y h y g cos θ sin ϕ m w y r y g cos θ cos ϕ ω x b h 2 + l 2 m ˙ 0 + w 1 , ω ˙ y b = 1 σ 2 m w x r x g cos θ cos ϕ + m w x h x + m w y h y g sin θ ω y b h 2 + l 2 m ˙ 0 + w 2 , ω ˙ z b = 1 σ 3 2 l 2 ω z b m ˙ 0 + w 3
where:
a 1 = M + m w x + m w y , σ 1 = J x x + m w x h x 2 + m w y r y 2 + h y 2 , σ 2 = J y y + m w y h y 2 + m w x h x 2 + r x 2 , σ 3 = J z z + m w x r x 2 + m w y r y 2
and:
d 1 = 1 a 1 m w x r ¨ w x + y ˙ b ω z b z ˙ b ω y b + 1 a 1 ω ˙ y b m w y h y + m w x h x + m w y ω z b r ˙ w x + ω ˙ z b m w y r w y F h x F s x + F d x , d 2 = 1 a 1 m w y r ¨ w y + ω x b z ˙ b ω z b x ˙ b + 1 a 1 + ω ˙ x b m w x h x + m w y h y m w x ω z b r ˙ w y m w x ω ˙ z b r x F h y F s y + F d y , d 3 = ω y b x ˙ b ω x b y ˙ b + 1 a 1 m w x ω y b r ˙ w x + ω ˙ y b r w x m w y ω x b r ˙ w y + ω ˙ x b r w y + p p A F h z F s z + F d z , w 1 = 1 σ 1 m w x h x + m w y h y y ¨ b m w y r y z ¨ b + ω ˙ z b m w x r x h x + ω x b ω y b m w x r x h x h i = 1 4 m ˙ i y ˙ b + m w y h y r y ω y b 2 ω z b 2 + ω y b ω z b m w x h x 2 + r x 2 m w x r x 2 + m w y h y 2 m w y r y 2 J z z + J y y + τ d x , w 2 = 1 σ 2 m w x h + m w y h y x ¨ b + ω z b ω x b m w x h x 2 + m w x r x 2 + m w y r y 2 m w y h y 2 + r y 2 J x x + J z z + m w x r x h x ω z b 2 ω x b 2 + m w x r x z ¨ b + a 6 + ω ˙ z b m w y h y r y + x ˙ b h i = 1 4 m ˙ i ω x b ω y b m w y h y r y + τ d y , w 3 = 1 σ 3 ω x b ω y b m w x h x 2 + r x 2 + m w x h x 2 + m w y h y 2 + r y 2 m w y h y 2 J y y + J x x + ω x b ω z b m w y h y r y m w x r x h x m w x r x y ¨ b + m w y r y x ¨ b + m w x ω ˙ x b r x h x + m w y ω ˙ y b r y h y m w y r y g sin θ m w x r x g cos θ sin ϕ + τ d z
Equation (11) lists the remaining terms of the system model in addition to Equation (9), including the accelerations of the actuation system, the higher-order terms, the input-output coupling terms, the effects of the water sprinkler and the water hose, and the external disturbance F d T τ d T T . It is worth noting that if the control inputs are the water flowrate and the position of the two weights, the yaw rotation is uncontrollable. Fortunately, it is stable, and the system can still operate with uncontrolled yaw rotation—for example, as in [28]. Furthermore, the effect of the acceleration of two weight blocks leads the system to a non-minimum phase phenomenon. That is, when a weight block moves in one direction, the robot initially moves backward before tilting and moving in that same direction.

3. Control System Design

3.1. LQI Control Law Design

The objective of the control system is to drive the robot from an initial position to follow the desired trajectory with a small tracking error, even in the presence of external disturbances and non-minimum phase phenomena. A servo control scheme with state feedback and an integral action is able to fulfill both requirements. The LQI is the servomechanism with optimal gain matrices that minimize the quadratic cost of the closed-loop control system. The cost function consists of the system states, integrator outputs, and control inputs with the corresponding coefficient matrices that can be easily tuned to reflect the weight of each element. Therefore, the LQI control system is designed based on the linearization of the system model in (9)–(11). Select a local equilibrium state at the robot that is unrotated and stabilized about the desired altitude z ¯ b . It is easily seen that the robot has to inject a constant flowrate m ¯ ˙ 0 such that the thrust force is equal to the total gravitational force at this position. Meanwhile, the two weights in the weight-shifting mechanism remain in their original position. Additionally, the yaw motion is uncontrollable but stable. Therefore, by denoting the state vector as X b = x b y b z b z ¯ b φ x φ y x ˙ b y ˙ b z ˙ b ω x b ω y b T , with φ x and φ y angular positions about xb- and yb-axes, and the control input vector as U = m ˙ 0 m ¯ ˙ 0 r w x r w y T , the system can be linearized as follows:
X ˙ b = A X b + B U + D , Y b = C X b = x b y b z b z ¯ b T
The matrices A, B, and C are given by:
A = O 5 I 5 0 0 0 0 g 0 0 0 0 h m ¯ ˙ 0 a 1 0 0 0 g 0 0 0 0 h m ¯ ˙ 0 a 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 m w x h x g + m w y h y g σ 1 0 0 h m ¯ ˙ 0 σ 1 0 m ¯ ˙ 0 h 2 + l 2 σ 1 0 0 0 0 0 m w x h x g + m w y h y g σ 2 h m ¯ ˙ 0 σ 2 0 0 0 m ¯ ˙ 0 h 2 + l 2 σ 2 , B = O 7 × 3 2 1 ρ A + 1 4 ρ a m ¯ ˙ 0 M ¯ + m w x + m w y 0 0 0 0 m w y g J x x + m w x h x 2 + m w y r ¯ y 2 + h y 2 0 m w x g J y y + m w y h y 2 + m w x h x 2 + r ¯ x 2 0 , C = I 3 O 3 × 7
I n is the identity matrix size n and D is the vector that contains the remaining nonlinear elements and unmodelled disturbances.
Moreover, given a desired trajectory r = X d Y d Z d T in the Earth-fixed frame, from the kinematic relationship in (1), the trajectory can be represented in the robot body frame as follows:
r b = R b E T r 0 t R ˙ b E ( τ ) T r ( τ ) d τ
An augmented variable of tracking error is defined as follows:
X ˙ i = r b Y b
Then, the system model, in combination with the augmented variable, is written as:
X ˙ b X ˙ i X ˜ ˙ = A O 10 × 3 C O 3 A ¯ X b X i X ˜ + B O 3 B ¯ U + D r b D ˜
A state-feedback control law is designed in the form of:
U = K f K i X ˜
where the optimal matrix gains are computed by the LQI control method:
K f K i K = R 1 B ¯ T P
such that the control law (17) minimizes the following quadratic cost function:
min U J U = 0 X ˜ T Q X ˜ + U T R U d t
Q and R are positive-definite weighting matrices of the state and input, respectively. P is the positive definite and symmetric matrix, which satisfies the algebraic Riccati Equation (20):
A ¯ T P + P A ¯ P B ¯ R - 1 B ¯ T P + Q = 0
The schematic drawing of the designed control system is depicted in Figure 4.
Assume that the vector of disturbances and unmodelled dynamics D ˜ is bounded by a non-negative value γ such that D ˜ 2 γ . To analyze the robustness of the closed-loop system, a Lyapunov function candidate is chosen as follows:
V = X ˜ T P X ˜
Taking the time-derivative of V and applying the control law in (17) and the algebraic Riccati equation in (20) results in the following:
V ˙ = X ˜ T A ¯ B ¯ K T P + P A ¯ B ¯ K X ˜ + D ˜ T P X ˜ + X ˜ T P D ˜ = X ˜ T P B ¯ R 1 B ¯ T P + Q X ˜ + D ˜ T P X ˜ + X ˜ T P D ˜ λ min P B ¯ R 1 B ¯ T P + Q X ˜ 2 + 2 γ λ max P X ˜ 2
where λ min ( . ) and λ max ( . ) indicate the minimum and maximum eigenvalues of the corresponding matrix. In order to preserve the system stability, V ˙ must be negative-definite. From Equation (22), one can see that V ˙ < 0 whenever the state satisfies the following condition:
X ˜ 2 > 2 λ max P λ min P B ¯ R 1 B ¯ T P + Q γ
which implies the boundedness of X ˜ .

3.2. GA Approach for Optimally Tuning the Controller Parameters

One can see that minimizing the right-hand side of Equation (23) would result in the smallest boundedness of the tracking errors, thus enhancing the robot’s performance and stability. However, better performance usually comes along with larger control energy. In the case of the proposed firefighting robot, the control energy includes the water pump energy and, especially, the kinetic energies of the two weight blocks in the weight shifting mechanism. From Equation (11), it is worth noting that the fast movement of the actuating weights makes the non-minimum phase phenomenon become worse and even causes the system to be unstable. In fact, if the acceleration of the actuating weights is separated from the vector D ˜ of the unmodelled dynamics, the boundedness γ of D ˜ can be considered as follows:
γ γ w + γ d , γ w = 1 a 1 m w x r ¨ w x m w y r ¨ w y 2 ,   γ d = D ˜ 1 a 1 m w x r ¨ w x m w y r ¨ w y O 1 × 1 T 2
where γ w contains the weights’ acceleration and γ d bounds the remainder. Minimizing the effect of γ w while minimizing the tracking error strikes a balance between enhancing the system performance and reducing the non-minimum phase phenomenon. In other words, let a cost function be proposed as follows:
min Q , R J Q , R = 2 λ max P λ min P B ¯ R 1 B ¯ T P + Q 1 + β m w x 0 t r ¨ w x 2 ( τ ) d τ + m w y 0 t r ¨ w y 2 ( τ ) d τ         r w x l , r w y l ;   0 m ˙ 0 m ˙ ¯ max otherwise
in which β = γ w / γ d is a trade-off between the performance and the non-minimum phase phenomenon. The conditions r w x l , r w y l , and 0 m ˙ 0 m ˙ ¯ max consider the saturation of the control input of the actual system. In particular, the left term of the cost is globally valid. Meanwhile, the other depends on the operating scenario. One can see that they are largest when the robot has to follow step-type references. Therefore, this scenario is going to be taken in the optimization process.
The objective of the optimization process is to find the controller matrices Q and R such that the cost function in (25) is minimized in the mentioned scenarios. The GA [29] is implemented for solving the optimization problem. The process is depicted in the following steps, as well as in Figure 5.
Step 1: Initialization
The first step of optimization is to create an initial population of the Q and R matrices. A pair ( Q i , R i ) of matrices, where i 1 , N , N is the population size, is an individual.
Step 2: Fitness evaluation
Each individual from the previous step is substituted in Equations (18) and (20) to calculate the P and K matrices, respectively. A simulation with the step-type trajectory is conducted with the above matrices, and the cost function (25) is obtained after finishing. A fitness function, which is an evaluation tool to determine the existence of the individual after each optimization generation, can also be defined as the inverse of the cost.
Step 3: Decimal Encoding
Each individual ( Q i , R i ) is encoded into a series of chromosomes si by using a decimal encoder.
Step 4: Natural selection
These fitness values are arranged in ascending order based on the linear ranking selection, in which the best individual is placed last and the worst one is in the first position. The jth individuals are selected for the next step via its accumulative probability pj, as described in Algorithm 1a. The higher a fitness value, the greater the probability of the individual being selected
Step 5: Crossover
Several individuals are randomly chosen to carry out the crossover. The number of conducted individuals is determined by a predefined crossover probability κ. Two crossover points in the chromosome series are selected randomly. A two-point crossover method, as in Algorithm 1b, swaps the chromosomes between two points from one individual to another.
Step 6: Mutation
As in Algorithm 1c, the mutation process arbitrarily changes one or more chromosomes, with a small probability η, in an individual to avoid premature convergence.
Step 7: Decode
The series of chromosomes si is decoded and returns the corresponding matrices Q i and R i .
Step 8: Fitness evaluation
The dependent matrices are computed in the same manner as in Step 2. The simulation is carried out to evaluate the best fitness value.
Step 9: Stop condition
The optimization process will be stopped if there is no change to the best fitness for a specified number of generations. Otherwise, the new loop will be continued.
Algorithm 1. The evolutionary process within the GA optimization
a.Linear Ranking Selection
   for k = 1:N
   r = rand(0,N);
     if (pj–1r < pj) then Ik = Ij
     end
   end
b.Two-point Crossover
   for m = 1:N
     k1 = rand(0,l); k2 = rand(0,l)
     if (κ > rand (0,1)) then
     Im(1:k1) := p1(1:k1); Im(k1 + 1:k2) := p2(k1 + 1:k2); Im(k2 + 1:k) := p1(k2 + 1:k)
     end
   end
c.Uniform Mutation
   for m = 1:N
     for k = 1:l
       if (η > rand(0,1)) then Im(k) = rand(0,9)
       end
     end
   end

4. Simulation Results and Discussion

4.1. Simulation

To evaluate the dynamical characteristics of the designed firefighting robot and validate the feasibility of the proposed LQI control system, two scenarios reflecting two common motion trajectories are considered. One is a step-type trajectory, which is similar to the case where the firefighting system moves from one point to another to battle individual flames. The other is a continuous flying path replicating the way in which the system overcomes obstacles and approaches the fire site. In simulation studies, each scenario is carried out with and without the influence of external forces and torques. These external terms represent the unformulated effects in Section 2, including the effects of the water sprinkler F s and the water hose F h and external disturbances F d T τ d T T . The latter include the effects of the wind gust, ceiling, ground effect, etc. in the practical operation of the system, but they are neglected in the modeling process. This is to validate the robustness of the proposed control system when operating the fire-extinguishing sprinkler or facing disturbances. The formulations of the external forces and torques are adopted from the previous publication of Nguyen and Hong [30], where they were also used to test the robustness of an aerial system. As mentioned in the previous section, the step-type trajectory is also used in the optimization process of finding the controller’s parameters.
Moreover, a cascade proportional—integral—derivative (PID) control structure from [24] is also compared to the proposed LQI control with the GA approach, in which the outer PID controllers use the horizontal tracking errors in X- and Y-directions to compute the robot’s required roll and tilt rotations, respectively. They are the references for the inner controllers, where the control inputs r w x and r w y are correspondingly obtained. In addition, another PID controller computes the required flowrate m ˙ 0 from the altitude control error of the system. Those control laws are mathematically represented in the following formulation:
U P I D = K p e + K I 0 t e d τ + K D d d t e , e = x d θ d y d ϕ d z d T x θ y ϕ z T , U P I D = θ d r w x ϕ d r w y m ˙ 0 T
The system model given by (9)–(11) is carried into the simulation. The parameters of the designed concept are shown in Table 1. The equilibrium point for linearization is at z ¯ b = 10 [ m ] height with respect to an injected flowrate m ¯ ˙ 0 = 2.7 [ kg / s ] to maintain this altitude for the system. The weighting matrices Q and R are found by the optimization process with the GA after 331 generations. The convergence of the cost function during the process is depicted in Figure 6. The feedback control gains Ki and Kf are computed as in Equations (18) and (20) with the previously obtained Q and R. They are shown in Table 2, where the PID controllers’ gains are also presented.

4.2. Results and Discussion

In the first scenario, the firefighting robot takes off, up to a height of 10 [m], and then moves horizontally to a length of 2 [m] in each of the X- and Y-directions, respectively. The geometry of the desired path is shown in Figure 7, along with the tracking performances provided by the two considered controllers: the proposed LQI-GA and the cascade PID. Subsequently, Figure 8 and Figure 9 depict the system response in every direction. The corresponding control inputs are given in Figure 10. For all of them, the left panels show the results when there is no disturbance and vice versa for the ones on the right.
In particular, Figure 7 shows that the designed robot can perform flying tasks as required. Moreover, the proposed LQI-GA manages the robot to closely follow the trajectory. The proposed controller is also robust against the disturbances, as the tracking performance is insignificantly affected in Figure 7b. The time responses in the positioning of the LQI-GA control system are also smooth, without any overshoot. The zoomed figures in Figure 8a prove that the non-minimum phase phenomena are smaller than those in the PID control system; even the PID controllers provide much slower performances. The external disturbances also heavily influence the PID-controlled system, as seen in the figures on the right. Figure 9 indicates that the PID controllers give the robot’s head less tilting and rolling. Nevertheless, the orientations are smoother with the LQI-GA controller, either with or without the presence of the disturbances. The yaw rotation of the robot is generally stable, without any control effort needed.
The smooth control inputs of the proposed controller explain the smooth response and the reduction of the non-minimum phase phenomena in the system. The smooth displacements rwx and rwy from the LQI-GA controller indicate the small acceleration of the weight blocks, thus resulting in smaller non-minimum phase phenomena. In contrast, the rough movements of the blocks in the PID control system lead to the fluctuation in orientations and the significant non-minimum phase phenomena.
In the second scenario, the firefighting robot follows a continuous trajectory while also maintaining system stabilization. The complex trajectory, in this case, is given by:
r = 10 e t 100 sin 0.1 u cos 0.01 u 0.1 t sin 0.05 t 0.1 t
The simulation results are presented in Figure 11, Figure 12, Figure 13 and Figure 14 in the same manner as the previous scenario. The tracking paths in Figure 11 show that the proposed control system remains the better one. Its resulting motion is closed to the reference, even with the presence of disturbances. Meanwhile, the PID controllers give a lot of overshoots and are easily influenced by external forces and torques. Interestingly, the PID control system follows the trajectory in real-time, while the LQI-GA is lagging behind, as seen in their time responses in Figure 12. However, Figure 13 and Figure 14 clarify that the PID control system requires a lot of rotational corrections and control efforts to achieve this result. They are important factors in evaluating the feasibility of the system in a real application. The integral of square values (ISVs) presented in Table 3 shows the consumed kinematic energies of the water and the weight-shifting actuators, in which the proposed LQI-GA is significantly more efficient than the other.
Mathematical analyses of the stability and performance of the proposed system have been given in the previous Section. In particular, the robustness of the closed-loop system has been proved via the Lyapunov function candidate (21), and the boundedness of the control error has been given in (23). The GA-based tuning controller gains minimize the cost function (25), which is equivalent to minimizing the control error, even in the presence of disturbances, and minimizing the variation of the weight-shifting actuation system. The proper performances in all simulation scenarios verify the theoretical analyses. The water flows smoothly, especially in comparison with the turbulent flow, due to the control of the nozzles’ flow in the study [24]. This is a great feature since it helps prevent the flutter instability of the water-conveying hose and enhances the total stability of the system.
Meanwhile, PID is a popularly standard control technique, and the cascade control structure has been adopted for aerial systems in many studies—for example, in the references [19,24,30] of the manuscript. The controlled quality depends on the PID controller parameters, which are not easily tuned in cascade control. The result is that the PID-controlled system seems to be slow in tracking the step-type reference but fast and overshot in tracking the sinusoid trajectory. Factors such as the saturation of the inputs and the objective of the reduction of the non-minimum phase phenomenon are difficult to incorporate into the design of the PID controllers. Thus, the PID control system requires a lot of rotational corrections to maintain the tracking performance, and noticeable chattering appears in the PID control inputs. That leads to significant control efforts and consumed water in the PID control system. On average, the LQI-GA consumes 2.28% less water and uses 83.85% less kinetic energy of the actuator than the PID control system.

5. Conclusions

This paper proposed a flying-type firefighting robot that is operated by waterpower and the weight shifting mechanism. The aim of this robot is to be efficiently used for accidents in water areas and other places where conventional firefighting methods are ineffective and the water sources are inexhaustibly available. The mathematical model of the flying firefighting robot was retrieved in detail.
The LQI law with the GA approach was proposed for the autonomous operation of the firefighting system. The LQI preserved the system’s stability. The GA optimization helps to find the gain matrices that minimize both the tracking error and the non-minimum-phase phenomenon. The simulation tests carried two scenarios, where the LQI-GA control system was compared to the cascade PID one from the previous study. The simulation results show the superiority of the proposed control system in terms of performance, robustness, and control effort.
At present, this conceptual design itself found the limits of a lack of generated yawing torque and an increase in the total mass due to two weights of the weight-shifting mechanism. For future works, the limits will be overcome by adding an actuator that is responsible for yaw motion control and using a unique weight for both lateral and horizontal motions. A prototype model is then going to be fabricated and experimented with for further validation. Other approaches for maneuvering the water-powered system are also being considered, such as regulating the injected flowrates or the direction of the nozzle outlets. Moreover, the effect of the water hose will be considered in further research to observe a system’s behavior fully. The system configuration can also be applied for aerial vehicles and other applications.

Author Contributions

Conceptualization, C.-T.D., T.H. and Y.-B.K.; methodology, C.-T.D. and T.H.; validation, T.H. and Y.-B.K.; formal analysis, T.H.; investigation, C.-T.D.; writing—original draft preparation, C.-T.D. and T.H.; writing—review and editing, T.H. and Y.-B.K.; supervision, 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).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

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. Miyazawa, K. Fire robots developed by the Tokyo fire department. Adv. Robot. 2002, 16, 553–556. [Google Scholar] [CrossRef]
  2. Kitagawa, A.; Tsukagoshi, H. Utilization of fluid power for rescue operations and robots. Adv. Robot. 2002, 16, 513–516. [Google Scholar] [CrossRef]
  3. Zhang, Q.; Ke, G. Kinematic analysis of fire-fighting robot under the impact of waterflow recoil force. In Proceedings of the 2015 IEEE International Conference on Mechatronics and Automation (ICMA), Beijing, China, 2–5 August 2015; pp. 264–268. [Google Scholar] [CrossRef]
  4. Howe & Howe Thermite. Available online: https://www.howeandhowe.com/civil/thermite (accessed on 7 July 2022).
  5. EmiControls. Available online: https://www.emicontrols.com/media/emicontrols_datasheet_fire_aircoretaf35_us_digital.pdf (accessed on 7 July 2022).
  6. Shark Robotics Colossus. Available online: https://www.shark-robotics.com (accessed on 7 July 2022).
  7. Mitsubishi Heavy Industries. Available online: http://www.mhi.com/news/190325.html (accessed on 7 July 2022).
  8. Liljeback, P.; Stavdahl, O.; Beitnes, A. SnakeFighter—Development of a Water Hydraulic Fire Fighting Snake Robot. In Proceedings of the 2006 9th International Conference on Control,, Automation, Robotics and Vision, Singapore, 5–8 December 2006; pp. 1001–1006. [Google Scholar]
  9. PopSci+. Available online: https://www.popsci.com/scitech/article/2008-03/firefighting-robot/ (accessed on 7 July 2022).
  10. Kim, J.-H.; Jo, S.; Lattimer, B.Y. Feature Selection for Intelligent Firefighting Robot Classification of Fire, Smoke, and Thermal Reflections Using Thermal Infrared Images. J. Sens. 2016, 2016, 8410731. [Google Scholar] [CrossRef]
  11. Ehang. Available online: https://www.ehang.com/ehang216f/ (accessed on 7 July 2022).
  12. Zhun. Available online: https://www.walkera.com/index.php/Goods/info/id/71 (accessed on 7 July 2022).
  13. Aerones. Available online: https://theindexproject.org/award/nominees/4690 (accessed on 7 July 2022).
  14. Goufei Aviation. Available online: https://mikeshouts.com/first-live-fire-fighting-drill-with-fire-fighting-drones/ (accessed on 7 July 2022).
  15. Li, R. Personal Propulsion Devices with Improved Balance. Spain Patent 2 659 042, 8 November 2014. [Google Scholar]
  16. Robinson, B. Water Propelled Flying Board. U.S. Patent 9,145,206, 30 October 2015. [Google Scholar]
  17. Zapata, F. Maneuvering and Stability Control System for Jet-Pack. U.S. Patent 2014/0103165, 17 April 2014. [Google Scholar]
  18. Li, R. Personal Propulsion Device 2007. U.S. Patent 7,900,867, 8 March 2011. [Google Scholar]
  19. Naidoo, J.E. Modelling, Design and Analysis of a Water Jetpack Powered by an Autonomous Underwater Vehicle (AUV) System. Master’s Thesis, University of KwaZulu-Natal, Pietermaritzburg, South Africa, 2017. [Google Scholar]
  20. 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]
  21. 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]
  22. 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]
  23. 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]
  24. 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]
  25. Niku, S.B. Introduction to Robotics: Analysis, Control, Applications; John Wiley & Sons: Hoboken, NJ, USA, 2020; ISBN 1119527627. [Google Scholar]
  26. Khalil, W.; Dombre, E. Modeling, Identification and Control of Robots; Kogan Page Science: London, UK, 2004; Chapters 5 and 12; pp. 85–115 and 291–311. [Google Scholar]
  27. Frank, M. White. In Fluid Mechanics, 8th ed.; SI UNITS; McGraw-Hill Education: New York, NY, USA, 2016. [Google Scholar]
  28. Liu, X.; Zhou, H. Unmanned Water-Powered Aerial Vehicles: Theory and Experiments. IEEE Access 2019, 7, 15349–15356. [Google Scholar] [CrossRef]
  29. Sivanandam, S.N.; Deepa, S.N. Introduction to Genetic Algorithms; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  30. Nguyen, X.-M.; Hong, S.K. Robust Backstepping Trajectory Tracking Control of a Quadrotor with Input Saturation via Extended State Observer. Appl. Sci. 2019, 9, 5184. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The applications of the proposed firefighting robot.
Figure 1. The applications of the proposed firefighting robot.
Applsci 12 09334 g001
Figure 2. The conceptual design of the head part of the flying firefighting robot.
Figure 2. The conceptual design of the head part of the flying firefighting robot.
Applsci 12 09334 g002
Figure 3. The dimension notation of the head part of the flying firefighting robot.
Figure 3. The dimension notation of the head part of the flying firefighting robot.
Applsci 12 09334 g003
Figure 4. The proposed LQI control system.
Figure 4. The proposed LQI control system.
Applsci 12 09334 g004
Figure 5. Flow chart of the GA optimization process.
Figure 5. Flow chart of the GA optimization process.
Applsci 12 09334 g005
Figure 6. Convergence of the cost function.
Figure 6. Convergence of the cost function.
Applsci 12 09334 g006
Figure 7. Tracking performances with the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Figure 7. Tracking performances with the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g007
Figure 8. Time responses of the robot position with the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Figure 8. Time responses of the robot position with the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g008
Figure 9. Time responses of the robot orientation with the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Figure 9. Time responses of the robot orientation with the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g009
Figure 10. Control inputs in tracking the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Figure 10. Control inputs in tracking the step-type trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g010
Figure 11. Tracking performances with the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Figure 11. Tracking performances with the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g011
Figure 12. Time responses of the robot position with the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Figure 12. Time responses of the robot position with the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g012
Figure 13. Time responses of the robot orientation with the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Figure 13. Time responses of the robot orientation with the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g013
Figure 14. Control inputs in tracking the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Figure 14. Control inputs in tracking the continuous trajectory. (a) Without disturbances; (b) With disturbances.
Applsci 12 09334 g014
Table 1. Parameters of the proposed firefighting robot.
Table 1. Parameters of the proposed firefighting robot.
ParameterDescriptionValueUnit
M h Mass of the head part3kg
JMoment of inertia of the head part d i a g 0.0071 ,   0.0085 ,   0.0133 Kg·m2
AInner diameter of the water hose0.05m
aInner diameter of the nozzle outlets0.01m
kCorrection factor of the water hose1.5-
l, hDimensions of the head part0.25, 0.1m
mwx, hxWeight and vertical location of WX0.5, 0.13 kg, m
mwy, hyWeight and vertical location of WY0.5, 0.08kg, m
Table 2. Controllers’ gains.
Table 2. Controllers’ gains.
ControllerController’s Gains
LQI-GA K f = 0 0 0 0 2.1669 2.7294 0 0 0 0 1.7038 3.5594 0 0 0 0 0 0 40.6451 28.6725 0 0 1.6960 3.8174 0 0 46.4869 39.9021 0 0 , K i = 0 0 0.6355 0.3487 0 0 0 0.3276 0
PID K P = d i a g 0.01 , 60 , 0.01 , 60 , 5 , K I = d i a g 0.0001 , 0.05 , 0.0001 , 0.01 , 0.27 , K D = d i a g 0.075 , 150 , 0.075 , 150 , 14
Table 3. Control effort in the second scenario.
Table 3. Control effort in the second scenario.
ISVWithout DisturbancesWith Disturbances
PIDLQI-GAPIDLQI-GA
0 t m ˙ 0 2 d t [ kg 2 / s 2 ] 8.0152 × 1057.8164 × 1058.0024 × 1057.7944 × 105
0 t r ˙ w x 2 d t [ m 2 / s 2 ] 0.47380.08510.47600.00850
0 t r ˙ w y 2 d t [ m 2 / s 2 ] 0.47380.08510.48590.0938
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/app12189334

AMA Style

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. Applied Sciences. 2022; 12(18):9334. https://doi.org/10.3390/app12189334

Chicago/Turabian Style

Dinh, Cao-Tri, Thinh Huynh, and Young-Bok Kim. 2022. "LQI Control System Design with GA Approach for Flying-Type Firefighting Robot Using Waterpower and Weight-Shifting Mechanism" Applied Sciences 12, no. 18: 9334. https://doi.org/10.3390/app12189334

APA Style

Dinh, C. -T., Huynh, T., & Kim, Y. -B. (2022). LQI Control System Design with GA Approach for Flying-Type Firefighting Robot Using Waterpower and Weight-Shifting Mechanism. Applied Sciences, 12(18), 9334. https://doi.org/10.3390/app12189334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop