System Modeling and Simulation of an Unmanned Aerial Underwater Vehicle

Unmanned Aerial Underwater Vehicles (UAUVs) with multiple propellers can operate in two distinct mediums, air and underwater, and the system modeling of the autonomous vehicles is a key issue to adapt to these different external environments. In this paper, only a single set of aerial rotors with switching propulsion abilities are designed as driving components, and then a compound multi-model method is investigated to achieve good performance of the cross-medium motion. Furthermore, some additional variables, such as water resistance, buoyancy and their corresponding moments are considered for the underwater case. In particular, a critical coefficient for air-to-water switching is presented to express these gradually changing additional variables in the cross-medium motion process. Finally, the sliding mode control method is used to reduce the altitude error and attitude error of the vehicles with external environmental disturbances. The proposed scheme is tested and the model is verified on the simulation platform.


Introduction
Autonomous unmanned systems are very popular all over the world due to their broad range of applications such as surveillance, inspection, fast delivery, search and rescue, among many others [1][2][3]. In recent years, the rapid development of composite materials and the continuous upgrading of intelligent control technology have been of great help in the development of unmanned aerial vehicles (UAVs) and unmanned underwater vehicles (UUVs) [4,5]. More challenging, the unmanned aerial underwater vehicle (UAUV) has become a novel research field, which aims to design equipment that can fly in the air and navigate underwater. To enhance the characteristics of these unmanned systems, many scholars and enterprises have done lots of relevant work in fuselage shape, fuselage material, wing layout, mass-volume ratio, power system design, controller design and their intelligent applications and so on [6][7][8][9].
The system modeling and control of UAVs are investigated in some works. Many kinds of UAVs, such as fixed-wing, multi-rotor wing and flapping wing, have been developed and analyzed about their different structures and mathematical description method [10,11], and a lot of achievements have been made. Compared with traditional air vehicles, the control of multi-rotor vehicles are more stable and efficient. Thus, the multi-rotor vehicles have significant advantages in vertical take-off and landing (VTOL), precise hovering and aggressive maneuvers. In addition, transition flight modeling of a fixed-wing VTOL UAV is also investigated, and the model is specifically tailored for the design of a hover to forward flight and forward flight to hover transition control system [12]. Currently, there is

Preliminaries
In this section, the mathematical model for the UAUV is illustrated. This model is basically obtained by representing the rotor-vehicle as a non-rigid body with adjustable rotor direction evolving in six dimensional space. Our aim is to provide the hybrid dynamic models of the cross-medium vehicle working in the air, air-to-water and underwater. Thus, the system dynamics can also be decomposed into three categories, detailed in Section 3. At the same time, the Newton-Euler equations are taken into account.

Reference Frames
This section first introduces the vehicle's mechanical structure and its operation principle. Traditional hybrid vehicles have a multi-rotor structure, but the rotors are divided into two groups, which are used as power sources in the air and in the water respectively. Inspired by this model, a novel air-underwater configuration with only one set of rotors is proposed, shown in Figure 1.

Preliminaries
In this section, the mathematical model for the UAUV is illustrated. This model is basically obtained by representing the rotor-vehicle as a non-rigid body with adjustable rotor direction evolving in six dimensional space. Our aim is to provide the hybrid dynamic models of the crossmedium vehicle working in the air, air-to-water and underwater. Thus, the system dynamics can also be decomposed into three categories, detailed in Section 3. At the same time, the Newton-Euler equations are taken into account.

Reference Frames
This section first introduces the vehicle's mechanical structure and its operation principle. Traditional hybrid vehicles have a multi-rotor structure, but the rotors are divided into two groups, which are used as power sources in the air and in the water respectively. Inspired by this model, a novel air-underwater configuration with only one set of rotors is proposed, shown in Figure 1. The vehicle flies in the air like a four-rotor aircraft. However, in underwater navigation, the four rotors reverse 180 degrees through the steering gear, forming upward thrust to overcome the buoyancy of the vehicle underwater, and change the motion attitude of the vehicle by controlling the thrust of the four rotors, thus changing the trajectory of the vehicle underwater. It presents simple mechanical structure and mathematical model, good maneuverability and controllability and allows hovering.
In order to accurately describe the movement of the vehicle and establish a mathematical model of the vehicle, the body frame The vehicle flies in the air like a four-rotor aircraft. However, in underwater navigation, the four rotors reverse 180 degrees through the steering gear, forming upward thrust to overcome the buoyancy of the vehicle underwater, and change the motion attitude of the vehicle by controlling the thrust of the four rotors, thus changing the trajectory of the vehicle underwater. It presents simple mechanical structure and mathematical model, good maneuverability and controllability and allows hovering.
In order to accurately describe the movement of the vehicle and establish a mathematical model of the vehicle, the body frame O b − x b y b z b and the inertial frame O e − x e y e z e are defined as shown in Figure 2. Note that O e is an East-North-Up (ENU) coordinate system and the axes of x e , y e and z e point to the east, north and upper directions respectively. In order to express consistency for the vehicle motion, the definition of z-axis coordinates is upwards in both aerial and underwater cases in this paper. The position of the vehicle is written as ξ = x y z T and the attitude angle of the vehicle is denoted as η = φ θ ψ T in the inertial frame O e . In addition, the angle's magnitude satisfies η < η L with the upper limit η L = π/2 π/2 2π T . This just means the vehicle can tilt the fuselage in motion, but the fuselage is not allowed to reverse. In the actual cases, the upper limits for the pitch and roll angles are usually less than π/2.
It can be seen that eb R is a direction cosine matrix expressed by the Euler angle of the unmanned vehicle.

Modeling of Multi-Rotor Vehicles
The mathematical models of four-rotor-driven vehicles are mostly described in the form of continuous differential equations. In general, it is assumed that the vehicle's center of gravity coincides with the body-fixed frame origin and each propeller rotates at the angular velocity i Ω (where the subscript i denotes the serial number of rotor, and 1, 2, 3, 4 i = ). Consequently, a force Fi is generated along or opposite to the z-direction relative to the body frame and a reaction torque Mi is produced on the vehicle body by each rotor expressed as: .
where the two positive coefficients ρ π denote the aerodynamic coefficient and the drag coefficient of the rotor respectively, and their magnitudes depend on the environmental medium density ρ , the radius of the propeller r, the thrust factor cT and torque factor Q c of the designed rotors [26]. Furthermore, in the sign function, h denotes the vehicle altitude to the air-water interface. This paper describes the attitude of the aerial underwater vehicle using the Euler angle, and the rotation matrices C φ , C θ and C ψ in the angle φ, θ and ψ are defined as: where the element s φ,θ,ψ are the sine function and c φ,θ,ψ are the cosine function of the rotation angles. Obviously, the transformation from the inertial coordinate system to the body coordinate system can be obtained as: It can be seen that R eb is a direction cosine matrix expressed by the Euler angle of the unmanned vehicle.

Modeling of Multi-Rotor Vehicles
The mathematical models of four-rotor-driven vehicles are mostly described in the form of continuous differential equations. In general, it is assumed that the vehicle's center of gravity coincides with the body-fixed frame origin and each propeller rotates at the angular velocity Ω i (where the subscript i denotes the serial number of rotor, and i = 1, 2, 3, 4). Consequently, a force F i is generated along or opposite to the z-direction relative to the body frame and a reaction torque M i is produced on the vehicle body by each rotor expressed as: where the two positive coefficients k T = c T ρr 4 π and k Q = c Q ρr 5 π denote the aerodynamic coefficient and the drag coefficient of the rotor respectively, and their magnitudes depend on the environmental medium density ρ, the radius of the propeller r, the thrust factor c T and torque factor c Q of the designed rotors [26]. Furthermore, in the sign function, h denotes the vehicle altitude to the air-water interface. The propellers represent the main source of thrust and can be designed specially in this hybrid configuration. Since the dynamic principle of the rotor structure is the same, aerial and aquatic propellers can present the same shape by optimizing the blade shape and size although the environments have different density. Therefore, the coefficients c T , c Q , the efficiency of engines and the energy delivered by the battery module have to be considered as the criteria to design the propeller shape (radius, width and curvature, et al.). Environmental density is directly proportional to the drag force and the added mass. Due to the large water density compared to the air density, it will generate high thrust forces with the same size and rotation speed to overcome the large movement resistance of the vehicles in the water. For convenience, the design process can be aided by the optimal propeller design software, such as OpenPro (OpenPro, Inc., Fountain valley, CA, USA), ShipPower (China state shipbuilding corporation limited, Shanghai, China).
The power system of the vehicle mainly includes rechargeable batteries, DC brushless motors and blades. Therefore, the magnitudes of the rotor's angular velocities decide the system's control inputs, which usually include the force and the torques created around a particular axis in body-fixed frame. In this section, the control signal can be defined as a vector T , and the components U 1~4 denote the total thrust, the roll moment, the pitch moment and the yaw moment of the vehicle respectively. According to the dynamic principle of actuators, these control inputs can further be expressed by the following equation with the square vector of the rotor's angular velocities as: where the coefficient K is a constant matrix for a symmetrical rigid multi-rotor vehicle. In addition, the arm's length is denoted as l, then the coefficient matrix K in (3) can be represented as the Kronecker product: Consequently, the dynamic model of the vehicle in frame O b can be obtained by the Newton's second law and theorem of moment of momentum as follows: where m denotes the mass, F b and M b denote the resultant force and moment acting on the multi-rotor vehicle, and v b = u v w T denotes the linear velocity and ω b = p q r T denotes the angular velocity. It is often convenient to let the origin of the frame O b coincide with the center of gravity, then the simplest form of the equations of motion is obtained when the body axes coincide with the principal axes of inertia. This implies that J = diag J x J y J z T , and J x , J y , J z are the inertia moments around x-axis, y-axis and z-axis respectively in the frame O b . If this is not the case, the body-fixed coordinates can be rotated about their axes to obtain a diagonal inertia tensor by simply performing a principal axis transformation [27]. Then we have: where (6), it also can be rewritten as: After that, by deriving the above formula, the dynamics of the vehicle can be expressed as: .. .. where denote the external disturbances for the force and the moment respectively, and the coefficient matrices A ξ , B ξ , A η can be represented as: In practical applications, the values of the position ξ and the attitude η can be measured by the sensors, such as Global Positioning Sysrtem (GPS), altimeter, gyroscope.

Dynamics of Cross-Medium Unmanned Aerial Underwater Vehicles (UAUVs)
By contrast with the single-medium operation environment, the cross-medium environment complicates the modeling analysis and stable continuous control of the system. The density ρ may jump aggressively from one value to another. In this case, the jump occurs when the vehicle's altitude h reaches the air-water interface, i.e., h := 0. In order to distinguish different external environments of the vehicle, a critical layer thickness of water surface is defined as ε, the value depends on the physical size of the vehicle and the model over-process requirements and ε → 0 . The mode switching for the cross-medium behavior of an aerial underwater vehicle is shown in Figure 3.
After that, by deriving the above formula, the dynamics of the vehicle can be expressed as: denote the external disturbances for the force and the moment respectively, and the coefficient In practical applications, the values of the position ξ and the attitude η can be measured by the sensors, such as Global Positioning Sysrtem (GPS), altimeter, gyroscope.

Dynamics of Cross-Medium Unmanned Aerial Underwater Vehicles (UAUVs)
By contrast with the single-medium operation environment, the cross-medium environment complicates the modeling analysis and stable continuous control of the system. The density ρ may jump aggressively from one value to another. In this case, the jump occurs when the vehicle's altitude h reaches the air-water interface, i.e., 0 := h . In order to distinguish different external environments of the vehicle, a critical layer thickness of water surface is defined as ε , the value depends on the physical size of the vehicle and the model over-process requirements and 0 → ε . The mode switching for the cross-medium behavior of an aerial underwater vehicle is shown in Figure 3. Therefore, operational modes of the cross-medium vehicles can be divided into three categories according to the current altitude value h, which denotes a relative altitude between the vehicle and the air-water interface.

Dynamics when
The aerial underwater vehicle is subjected to propeller thrust, air resistance, counter-torque generated by the propeller, and gyro moment generated by the attitude change during the movement. However, the air resistance and gyro moment can usually be ignored due to the small air density and small attitude change. Thus, only the propeller thrust, counter-torque and the gravity of the vehicle are considered during the vehicle's motion in the air.  Therefore, operational modes of the cross-medium vehicles can be divided into three categories according to the current altitude value h, which denotes a relative altitude between the vehicle and the air-water interface.

Dynamics when h > ε
The aerial underwater vehicle is subjected to propeller thrust, air resistance, counter-torque generated by the propeller, and gyro moment generated by the attitude change during the movement. However, the air resistance and gyro moment can usually be ignored due to the small air density and small attitude change. Thus, only the propeller thrust, counter-torque and the gravity of the vehicle are considered during the vehicle's motion in the air.
According to the reference frame definition in Section 2, the thrust generated by the propellers is i and the direction is upward, same to the z b axis. The motion equation of the vehicle in the inertial coordinate system can be derived as: The angular velocity of the vehicle in the inertial frame O e : Therefore, in the actual control, the motor's speed output can be calculated. Finally, combined with the Formulas (3), (9) and (10), the motion model of the vehicle when h > ε satisfies:

Dynamics when h < −ε
The unmanned vehicle navigates in the water in a similar way to the air movement. The underwater vehicle is subjected to gravity, buoyancy, propeller thrust and resistance during movement. In addition, the moments subjected include heavy moments, floating moments, propeller thrust moments, resistance moments and gyro moments generated when the attitude changes.
Underwater vehicles are normally designed to satisfy the condition that the buoyancy force is slightly larger than the weight, which is important for safety reasons, and increasing weight or upward thrust is necessary for submersion. Therefore, in our configuration, the propeller arms of the vehicle will rotate by 180 • after entering the water, so that the thrust of the propeller is reversed in an opposite direction to the z axis in the water, which is different from the convention used for some existing underwater vehicles.
In general, the motion of an underwater vehicle moving in six degrees of freedom (6DOF) at high speed will be highly nonlinear and coupled. However, in many underwater applications the vehicle will only be allowed to move at low speed, thus the Coriolis effects are smaller compared to the inertia effects of added mass, and the Coriolis effects of added mass are not considered in this paper. In addition, added mass can be seen as pressure-induced forces and moments due to a forced harmonic motion of the body which are proportional to the acceleration of the body. Since the movement of the vehicle in the water is also affected by the added mass, and the vehicle is approximately considered symmetrical, then the added mass along the three axes of the vehicle is only related to the acceleration and angular acceleration along the respective axes, and it can be expressed as: where m k represents the mass of the object in motion, m . u , m . v , m . w represent the added mass of axial motion in water, and m . p , m . q , m . r represent the added mass of rotating motion in water. In the inertial frame O e , the gravity of the vehicle in the water is in the opposite direction to the z-axis, but the buoyancy direction is the same as the z-axis. Furthermore, the matrix expression of the gravity and the buoyancy are converted to the body frame O b as: where ρ w is the density of water and V is the effective volume of the vehicle. Due to the high density of water, the resistance of the vehicle when moving in water cannot be ignored. In the frame O b , the water flow resistance force F b rw can be expressed as: where C dw is the dimensionless drag coefficient (for most underwater robots, C dw = 0.8 ∼ 1.0), S = diag S x S y S z is the area in the direction of water flow velocity. The thrust generated by the propellers is F b mw = −F b m , it's direction is contrary to the z b axis. In summary, the external force F b w of the UAUV can then be represented as: The aerial underwater vehicle changes its attitude in the water by its torque. Since the origin of the body coordinate system is located at the center of gravity of the vehicle, the gravity does not generate the moment. The buoyancy center coordinates are assumed to be located on the z b axis in the frame O b , and the buoyancy center coordinates are x f , y f , z f T in the frame O e . Then the buoyancy moment satisfies: where F b b is the buoyancy of the vehicle underwater in the body coordinate system. The thrust torque generated by the propeller in the water is similar to that in the air, but the thrust torque generated by the propeller in the water is opposite to that in the air due to the 180 • rotation of the propeller, i.e., M b mw = −M b m show in the 2~4th rows of the Formula (4). In addition, the change of the moving attitude in the underwater vehicle will produce the gyro moment as: where J rw is the moment of inertia of the underwater motor and propeller. If the aerial underwater vehicle changes its attitude in the water, it will also be affected by the resistance moment, which is proportional to the square of the angular velocity of the aerial underwater vehicle as M b rw = k x p p k y q q k z |r|r T , where k x , k y and k z are the resistance moment coefficients around the x b , y b and z b axes in the body coordinate system. In summary, the external moment of the vehicle in the water can then be represented as: If the aerial underwater vehicle fuselage is symmetrical about each coordinate system, the inertia matrix of its movement in the water can be set to J w = diag J xxw J yyw J zzw . Similarly, the translational and rotational dynamics of the UAUV in water can be obtained by the Formulas (5), (15) and (18) as:

Dynamics when −ε ≤ h ≤ ε
For water and air cross-medium unmanned vehicles, the key processes include air-water transition (surface landing) and water-air transition (surface takeoff). How to achieve a stable and reliable conversion between the water and the air is significant. The unmanned vehicle can splash into water, this method has certain advantages in both the water entering time and the falling distance. It is a suitable way for the air-water transition of the amphibious cross-medium vehicle.
The unmanned vehicle is in the semi-submersible dynamic process when −ε ≤ h ≤ ε. This kind of vehicle resolves the problem of sailing across the different mediums. The additional variables, such as water resistance, buoyancy, resistance moment, buoyancy moment and gyro moment, cannot be neglected. They are modeled for the completely underwater vehicle in Section 3.2. However, the values of these variables are related to the vehicle altitude, which are not same to the underwater case.
In order to model the dynamics of the UAUV in the environmental transition stage, a critical coefficient of the additional variables is defined as: Thus, the critical coefficient k is inverse proportional to the altitude h, and obviously its magnitude satisfies k s ∈ [0, 1]. If the vehicle decreases to below the air-water interface, then k s increases and the additional variables' function will augment gradually to the underwater case, and vice versa. Similarly, the dynamics can be deduced from the Formula (19) as:

Controller Design
In this section, a non-linear control scheme is proposed, which provides the integrated control method of the air and underwater dynamics of the UAUV. The aerial underwater vehicle is a typical under-actuated system, and it faces a nonlinear and strong coupling environment. The advantages of sliding mode control are obvious, which is used to design the altitude controller and attitude controller of this unmanned system. Figure 3 shows the integrated control block diagram for the UAUV.
The deviations of altitude and attitude are the inputs of the SMC controller, and the controller outputs are the system inputs U 1~4 . According to altitude h of the current vehicle distance from air-water interface, the dynamics of the multi-model system discussed in Section 3 is determined, and the angular velocities of the vehicle motor are calculated. Thus, a closed-loop nonlinear multi-model feedback control scheme is formed. The system outputs of the designed scheme mainly include the position and the attitude angles which can be measured directly.
In the above scheme depicted in Figure 4, the position loop is designed to attenuate the tracking position error which is expressed as ξ e = [ x e y e z e ] T . In addition, we substitute the position error ξ e into the Equation (7) and neglect the external disturbance temporarily, the following resultant error equation is deduced as: .
. where z λ is a positive constant and ( ) sat ⋅ is a continuous saturation function defined as: where δ is the thickness of the boundary layer. For our sliding mode method, when the control error enters the given boundary layer, it will slip into the vicinity of zero in the boundary layer. Thus the saturation function is used to smooth out the control discontinuity around zero to reduce undesired chattering caused by imperfection switching of the discontinuous term. The performance of robustness is dependent on the thickness of the boundary layer, that means a thicker boundary layer may not be helpful to the robustness of the system. In addition, the desired attitude can be calculated from the position errors, and the attitude error is obtained. After that, we consider the attitude loop used for attenuating the tracking attitude error where 1 k η , 2 k η are constant control parameters.

Simulations
In this section, the dynamics and the proposed control method are verified on a simulation platform of an autonomous UAUV with 4 rotors. The numerical parameters used in the simulation are a mixture of previous, similar projects and measurements of the real platform. The mass of the simulation vehicle is 1.2 kg, and the length of its arm is 0.19 m, with a small body-volume of 10 −3 m 3 , the maximum power of the motor is 202.8 W and the maximum rotational speed of the motor is 8500 rpm respectively. The detailed parameters used in the simulator are shown in Appendix A. In order to express clearly the work flow of the system, a flow diagram of the motion control process of the vehicle is presented in Figure 5.
At first we simulated the air position response and the underwater position response of the vehicle in the air and in the water single-medium environments in Sections 5.1 and 5.2, respectively. Then, to verify the cross-medium characteristics of the unmanned vehicle, a cross-media control Thus, the altitude sliding manifold of the vehicle can be defined as: where two coefficients k z 1 and k z 2 are constants. According to the SMC method, the system states are desired to remain on the manifold defined as σ z = 0, and the error dynamics can be written as k z 1 . ξ e + k z 2 .. ξ e = 0. After that, we substitute the Equation (22) into this altitude second-order equation and get the equivalent control formula as: So long as the system dynamics are accurately modeled, U eq 1 can maintain the altitude in the absence of external disturbances.
However, some unknown disturbances have to be considered in the real system, and the dynamics (22) can be represented as: ..
where d ξ denotes the external disturbance mentioned before. To keep the system states still on the defined sliding mode manifold, a feedback control law of the system is redesigned as: where λ z is a positive constant and sat(·) is a continuous saturation function defined as: where δ is the thickness of the boundary layer. For our sliding mode method, when the control error enters the given boundary layer, it will slip into the vicinity of zero in the boundary layer. Thus the saturation function is used to smooth out the control discontinuity around zero to reduce undesired chattering caused by imperfection switching of the discontinuous term. The performance of robustness is dependent on the thickness of the boundary layer, that means a thicker boundary layer may not be helpful to the robustness of the system. In addition, the desired attitude can be calculated from the position errors, and the attitude error is obtained. After that, we consider the attitude loop used for attenuating the tracking attitude error η e = φ e θ e ψ e T , Substituting the error η e into (8) and neglecting the environmental disturbance d η , we get .. η e = A η U η . Assume λ η is a positive constant, the corresponding control law obtained can be given by: where k η 1 , k η 2 are constant control parameters.

Simulations
In this section, the dynamics and the proposed control method are verified on a simulation platform of an autonomous UAUV with 4 rotors. The numerical parameters used in the simulation are a mixture of previous, similar projects and measurements of the real platform. The mass of the simulation vehicle is 1.2 kg, and the length of its arm is 0.19 m, with a small body-volume of 10 −3 m 3 , the maximum power of the motor is 202.8 W and the maximum rotational speed of the motor is 8500 rpm respectively. The detailed parameters used in the simulator are shown in Appendix A. In order to express clearly the work flow of the system, a flow diagram of the motion control process of the vehicle is presented in Figure 5.
directions. Otherwise, the water-air motion triggers the acceleration adjustment of the four rotors and reverses the rotor directions similarly. In addition, for the case h > 0.25 and h < −0.25, the vehicle operates in the air mode and in the underwater mode respectively. Furthermore, the sliding mode controllers are used so as to control the vehicle's position and attitude to infinitely close to the desired targets, and the desired attitude is calculated from the position errors of the vehicle, then the motor speeds are derived to drive the vehicle.   At first we simulated the air position response and the underwater position response of the vehicle in the air and in the water single-medium environments in Sections 5.1 and 5.2, respectively. Then, to verify the cross-medium characteristics of the unmanned vehicle, a cross-media control motion was performed in Section 5.3. Consequently, the vehicle navigation behavior is mainly divided into three cases: the air case, the underwater case and the transition case. In our platform, a critical layer thickness ε of the air-water interface is set to 0.25 m. Then for the case −0.25 < h < 0.25, the vehicle is in the transition mode. The air-water motion triggers the deceleration adjustment of the four rotors, and the rotors will reverse 180 degrees by the steering gears to change the thrust directions. Otherwise, the water-air motion triggers the acceleration adjustment of the four rotors and reverses the rotor directions similarly. In addition, for the case h > 0.25 and h < −0.25, the vehicle operates in the air mode and in the underwater mode respectively. Furthermore, the sliding mode controllers are used so as to control the vehicle's position and attitude to infinitely close to the desired targets, and the desired attitude is calculated from the position errors of the vehicle, then the motor speeds are derived to drive the vehicle.

Air Position Response
In this case, we tested a straight-line flight control of the vehicle in the air. Set the original position O(0, 0, 0) of the vehicle on the water surface and the target position at the inertia coordinates T (5,5,5) in the air, and the route is a straight way from O to T. Consequently, the position response curve and attitude response curves can be obtained, shown in Figures 6 and 7, respectively.
It can be seen that the aerial underwater vehicle can reach the specified position in a short time. Furthermore, this illustrates that although there is an overshoot of about 10% during its movement, the specified position can be reached within 5 s, and the actual path followed by the vehicle is very close to the desired straight-line path. In addition, the peak values of roll angle and pitch angle is 12.6 • and 11.5 • respectively, while the adjustment of yaw angle is relatively small.
The control inputs U i (i = 1,2,3,4) of the vehicle are shown in the following Figure 8. In the initial acceleration stage, the upward thrust force U 1 increases rapidly to about 40 N, then falls back gradually and stabilized on 12 N which exactly counteracts its own gravity and maintained hovering state after vehicle reaches the target point. Furthermore, the moments U 2-4 of the three coordinate axes are also effective within the initial 5 s, which drive the vehicle to move to the target point by attitude adjustment during its ascending process. The outputs of the actuators for the system, i.e., the rotational speeds of the vehicle motors, are also shown in Figure 9.   It can be seen that the aerial underwater vehicle can reach the specified position in a short time. Furthermore, this illustrates that although there is an overshoot of about 10% during its movement, the specified position can be reached within 5 s, and the actual path followed by the vehicle is very close to the desired straight-line path. In addition, the peak values of roll angle and pitch angle is 12.6° and 11.5° respectively, while the adjustment of yaw angle is relatively small.
The control inputs Ui (i = 1,2,3,4) of the vehicle are shown in the following Figure 8. In the initial acceleration stage, the upward thrust force U1 increases rapidly to about 40 N, then falls back gradually and stabilized on 12 N which exactly counteracts its own gravity and maintained hovering state after vehicle reaches the target point. Furthermore, the moments U2-4 of the three coordinate axes are also effective within the initial 5 s, which drive the vehicle to move to the target point by attitude adjustment during its ascending process. The outputs of the actuators for the system, i.e., the rotational speeds of the vehicle motors, are also shown in Figure 9.  It can be seen that the speeds of the motors increase rapidly to about 6800-7000 rpm, then decrease gradually to less than 4000 rpm and remain unchanged, which corresponds to the hovering state of the UAV. The speed differences of the four motors are not obvious because the upward thrust of the system is relatively large and the rotation moments are relatively small. When the UAV is in  It can be seen that the speeds of the motors increase rapidly to about 6800-7000 rpm, then decrease gradually to less than 4000 rpm and remain unchanged, which corresponds to the hovering state of the UAV. The speed differences of the four motors are not obvious because the upward thrust of the system is relatively large and the rotation moments are relatively small. When the UAV is in It can be seen that the speeds of the motors increase rapidly to about 6800-7000 rpm, then decrease gradually to less than 4000 rpm and remain unchanged, which corresponds to the hovering state of the UAV. The speed differences of the four motors are not obvious because the upward thrust of the system is relatively large and the rotation moments are relatively small. When the UAV is in the vertical take-off and landing motion and hovering state, the four motors have the same speed.

Underwater Position Response
In this case, we tested a flight motion of the vehicle in the air. Set the original position O(0, 0, 0) of the vehicle on the water surface and the underwater target position at the inertia coordinates T(−1, −1, −1), and the route is a straight way from O to T. Consequently, the position response curve and attitude response curves can be obtained, as shown in Figures 10 and 11, respectively. It can be seen that the aerial underwater vehicle can reach the set point steadily in the water. Because the water density is 1000 kg/m 3 , which is about 1000 times larger than the air density, it can produce high thrust with reduced motor speed. Furthermore, the peak values of roll angle and pitch angle are 16 • and 24.2 • respectively, while the adjustment of yaw angle is relatively small during the vehicle's movement in water as shown in Figure 11.
The control inputs U i (i = 1,2,3,4) of the vehicle are shown in the following Figure 12. In the initial acceleration stage, the upward thrust force U 1 increases rapidly to about 30 N, then falls back gradually and stabilizes on 3.5 N which is used to counteract the effects of its own gravity, buoyancy and water resistance. Furthermore, the moments U 2-4 of the three coordinate axes are also effective within the initial 5 s, which drive the vehicle to move to the target point by attitude adjustment during its descending process. The outputs of the actuators for the system, i.e., the rotational speeds of the vehicle motors are also shown in the Figure 13. It can be seen that the speeds of the motors increase rapidly to about 1410 rpm, then decrease gradually to less than 500 rpm and remained unchanged. However, the motor speeds are slower (one eighth) than those in the air.  24.2° respectively, while the adjustment of yaw angle is relatively small during the vehicle's movement in water as shown in Figure 11.    The control inputs Ui (i = 1,2,3,4) of the vehicle are shown in the following Figure 12. In the initial acceleration stage, the upward thrust force U1 increases rapidly to about 30 N, then falls back gradually and stabilizes on 3.5 N which is used to counteract the effects of its own gravity, buoyancy and water resistance. Furthermore, the moments U2-4 of the three coordinate axes are also effective within the initial 5 s, which drive the vehicle to move to the target point by attitude adjustment during its descending process. The outputs of the actuators for the system, i.e., the rotational speeds of the vehicle motors are also shown in the Figure 13. It can be seen that the speeds of the motors increase rapidly to about 1410 rpm, then decrease gradually to less than 500 rpm and remained unchanged. However, the motor speeds are slower (one eighth) than those in the air.

Aerial Underwater Vehicle Cross-Media Response
In this case, we test a cross-medium motion near the air-water interface, and we set the desired motion trajectory: when the altitude of the vehicle satisfies | | 0.25 z ≥ . Otherwise, the aerial underwater vehicle is in vertical lift transition. As

Aerial Underwater Vehicle Cross-Media Response
In this case, we test a cross-medium motion near the air-water interface, and we set the desired motion trajectory: when the altitude of the vehicle satisfies | | 0.25 z ≥ . Otherwise, the aerial underwater vehicle is in vertical lift transition. As

Aerial Underwater Vehicle Cross-Media Response
In this case, we test a cross-medium motion near the air-water interface, and we set the desired motion trajectory: x d = sin(0.5t), y d = sin(0.25t), z d = sin(0.25t)+0.25 · sign(z) when the altitude of the vehicle satisfies |z|≥ 0.25 . Otherwise, the aerial underwater vehicle is in vertical lift transition. As can be seen from the simulation result of the three-dimensional cross-medium trajectory shown in Figure 14, the vehicle first flies an arc in the air, then lands and enters the water from the point O, and dives through an arc trajectory, returning to the vicinity of the starting point. Moreover, point O is the incidence point on the water surface, and point P a and point P b are the vehicle state adjustment points respectively, and the vehicle moves vertically when its altitude is within the range [−0.25, 0.25]. Therefore, the vehicle can travel across the medium and track the predetermined trajectory along the arrow directions. In addition, Figure 15 shows that the attitude angles are adjusted to be small when 13 s t ≤ at initial stage. However, the roll angle and pitch angle fluctuate slightly at t = 13 s and t = 35 s, which corresponded to the vehicle's critical state adjustments of water-entering and water-exiting respectively. However, at t = 22.5 s, these two attitude angles suddenly increase to 0.4-0.5 rad and last for about 12.5 s, which corresponds to the vehicle's motion in water. In addition, Figure 15 shows that the attitude angles are adjusted to be small when t ≤ 13 s at initial stage. However, the roll angle and pitch angle fluctuate slightly at t = 13 s and t = 35 s, which corresponded to the vehicle's critical state adjustments of water-entering and water-exiting respectively. However, at t = 22.5 s, these two attitude angles suddenly increase to 0.4-0.5 rad and last for about 12.5 s, which corresponds to the vehicle's motion in water.
The control inputs U i (i = 1,2,3,4) of the vehicle are shown in the following Figure 16. In the initial stage, the upward thrust force U 1 is adjusted and stabilized on 12 N until t = 13 s. Then, the thrust force decreases to near zero and the vehicle's arms are in the state of steering adjustment, and then the force gradually increases to about 16N at t = 20 s and lasts for a period of time. Thereafter, the reverse thrust force adjustment is made at t = 35 s, and the final effluent operation is completed. Furthermore, the moments U 2-4 of the three coordinate axes increase at about t = 20 s and the maximum is not more than 1 Nm, which drives the vehicle to adjust its attitude during the underwater motion. The rotational speeds of the vehicle motors are also shown in Figure 17. The motor speeds are approximately 4020 rpm and 560 rpm for the air condition and underwater condition respectively. In addition, all the rotational speeds of the motor have to undergo a process of first deceleration and then acceleration in the process of critical state adjustment. In addition, Figure 15 shows that the attitude angles are adjusted to be small when 13 s t ≤ at initial stage. However, the roll angle and pitch angle fluctuate slightly at t = 13 s and t = 35 s, which corresponded to the vehicle's critical state adjustments of water-entering and water-exiting respectively. However, at t = 22.5 s, these two attitude angles suddenly increase to 0.4-0.5 rad and last for about 12.5 s, which corresponds to the vehicle's motion in water. The control inputs Ui (i = 1,2,3,4) of the vehicle are shown in the following Figure 16. In the initial stage, the upward thrust force U1 is adjusted and stabilized on 12 N until t = 13 s. Then, the thrust force decreases to near zero and the vehicle's arms are in the state of steering adjustment, and then the force gradually increases to about 16N at t = 20 s and lasts for a period of time. Thereafter, the reverse thrust force adjustment is made at t = 35 s, and the final effluent operation is completed. Furthermore, the moments U2-4 of the three coordinate axes increase at about t = 20 s and the maximum is not more than 1 Nm, which drives the vehicle to adjust its attitude during the  Figure 17. The motor speeds are approximately 4020 rpm and 560 rpm for the air condition and underwater condition respectively. In addition, all the rotational speeds of the motor have to undergo a process of first deceleration and then acceleration in the process of critical state adjustment.

Conclusions
This study evaluates the cross-medium motion of a quadrotor-like UAUV. The goal is to model the vehicle and simulate the autonomous control process both in the air and in the underwater environments. A hybrid vehicle configuration with only 4 rotors is proposed, and the thrust direction

Conclusions
This study evaluates the cross-medium motion of a quadrotor-like UAUV. The goal is to model the vehicle and simulate the autonomous control process both in the air and in the underwater environments. A hybrid vehicle configuration with only 4 rotors is proposed, and the thrust direction changes by the steering gears installed at the arm junctions. The simple mechanism makes the vehicle more flexible and reliable. The use of a single set of rotors in both mediums greatly reduces cost, weight, and complexity compared to other aerial underwater vehicle concepts. Without losing generality, the body frame and the inertial frame are defined to describe the position and attitude of the vehicle, and the transformation between these two frames is expressed as a direction cosine matrix of the vehicle's Euler angles.
The thrust force and the reaction torque generated by each rotor are proportional to the square of the propeller's angular velocities. Furthermore, the aerodynamic coefficient and the drag coefficient depend on the medium density, the radius of the propeller, the thrust and the torque coefficient. In our configuration, the aerial and aquatic propellers present the same shape, and the large water density generates high thrust forces with the same size and rotation speed to overcome the large movement resistance of the vehicles in the water. After that, the Quaternary control force and torque of the vehicle are defined. The system dynamics of the vehicle in the body frame are summarized and reconstructed as a second-order system based on Newton's second law and the theorem of the moment of momentum.
For the cross-medium motion the of aerial underwater vehicle, a critical layer thickness of the water surface is first defined and the dynamics model switches according to the vehicle's altitude h to the air-water interface. In our configuration, the thrust of the propeller is downward in the air to counteract the forces of gravity and acceleration, and is upward to submerge underwater. For the air environments, the air resistance and gyro moment are usually ignored to simplify the vehicle's dynamics expression due to the small air density and small attitude change. However, the buoyancy, the resistance and the corresponding resistance moments and gyro moments generated when the attitude changes have to be considered in the dynamics due to the large water density. Besides, the added mass is considered and its magnitude is related to the acceleration and angular acceleration along the respective axes. Then the external force and moment are derived and represented. In particular, for the air-water and water-air transitions, a critical coefficient is presented to model the changed mass, force and moment in the cross-medium motion process. In addition, the sliding mode control method is used to design the altitude controller and attitude controller of this unmanned system, and a continuous saturation function is performed to reduce undesired chattering caused by imperfection switching of the method. Finally, the deduced dynamics and proposed control scheme are tested on the simulation platform of an autonomous UAUV in different scenarios.
In future work, it is important to apply our modeling and control method to the real cross-medium experimental platform, and study how to improve the vehicle's maneuverability and stability, as well as investigating the influence of wind, wave, current and other disturbances on the system.