Variable Structure Controller for Energy Savings in an Underwater Sensor Platform

This paper introduces a new variable structure controller designed for depth control of an autonomous underwater sensor platform equipped with a variable buoyancy module. To that end, the prototype linear model is presented, and a finite element-based method is used to estimate one of its parameters, the hull deformation due to pressure. To manage potential internal disturbances like hull deformation or external disturbances like weight changes, a disturbance observer is developed. An analysis of the observer steady-state estimation error in relation to input disturbances and system parameter uncertainties is developed. The locations of the observer poles according to its parameters are also identified. The variable structure controller is developed, keeping energy savings in mind. The proposed controller engages when system dynamics are unfavorable, causing the vehicle to deviate from the desired reference, and disengages when dynamics are favorable, guiding the vehicle toward the target reference. A detailed analysis determines the necessary switching control actions to ensure the system reaches the desired reference. Finally, simulations are run to compare the proposed controller’s performance with that of PID-based controllers recently developed in the literature, assessing dynamic response and energy consumption under various operating conditions. Both the VBM- and propeller-actuated vehicles were evaluated. The results demonstrate that the proposed controller achieves an average energy consumption reduction of 22% compared to the next most efficient PID-based controller for the VBM-actuated vehicle, though with some impact on control performance.


Introduction
The ocean is home to a wide array of biological entities that require in-depth study and contains abundant natural resources such as oil and gas, which need to be responsibly extracted.Addressing environmental pollution in the ocean is also a crucial global issue today.However, ocean waters are heavily influenced by currents, waves, and other dynamic factors, making the underwater environment highly complex and potentially hazardous.As a result, tasks such as underwater pipeline inspections, offshore infrastructure upkeep, seabed and deep-sea exploration, collecting samples from marine environments, and monitoring and gathering data on ecological aquatic phenomena are very complex, timeconsuming, and expensive.To tackle these challenges, it is crucial to design and develop underwater robots that meet high survivability standards while being user-friendly and capable of autonomous operation [1].
Depth control is critical in autonomous underwater vehicles (AUVs) due to its direct impact on the vehicle's ability to perform precise and effective underwater tasks [2,3].
Sensors 2024, 24, 5771 2 of 19 Accurate depth control ensures that AUVs can maintain stable navigation and positioning, which is essential for collecting reliable data in scientific research, conducting detailed underwater surveys, and executing complex industrial operations like pipeline inspections and offshore infrastructure maintenance.Furthermore, effective depth control is necessary for avoiding underwater obstacles and adapting to varying oceanographic conditions, such as changes in water density and current patterns.While precise depth control is essential, it must also consider the associated energy consumption, since energy efficiency in AUVs is vital for enhancing their operational effectiveness and extending mission duration.In this context, low-energy control strategies that can completely deactivate the actuators are of particular interest for applications requiring the sensor platform to maintain a minimal acoustic footprint.Examples of this situation occur whenever military [4,5], shipping [6,7], or seismic [8,9] activities should be traced, or, from a biology perspective, whenever whale population density should be identified [10] or sounds from fish [11], shrimp [12], or marine mammals [13] should be recorded, as they provide significant information on an ecosystem's health.Also, sea surface wind speed [14] and ice shelf calving [15] can be monitored using underwater acoustics, for which it is important that the sensor platform has the minimum self-generated sounds.
One key topic in this area is the decision of using propellers or a variable buoyancy module (VBM) for depth control.In [16], a comparison between both actuators is made using a numerical example of an underwater glider.The authors of [16] conclude that propellers can match or even outperform the efficiency of a buoyancy engine.However, in the authors' opinion, there are two shortcomings that might compromise the conclusions of the study [16]: (i) the force used to calculate the energy spent by the propeller-actuated glider does not seem comparable to the one used when considering a variable buoyancy system; (ii) the study does not include any hovering operations, in which a slight buoyancy offset may give a considerable advantage to a buoyancy engine over a propeller solution.Also, in this context, the authors of the present paper presented in [17] a detailed model of hydraulic and electromechanical variable buoyancy systems.These models were applied to numerous numerical examples based on typical mission profiles.The paper results demonstrated that using a VBS offers significant advantages over vertical propellers regarding energy consumption for hovering tasks.In fact, even at relatively shallow depths and shorter mission durations (contrary to what might be intuitively expected), the use of a VBS, particularly the electromechanical variant, proves more beneficial than relying on propellers.
Another interesting approach can be found in [18], where a hybrid passive buoyancy compensation system (HPBCS) for underwater gliders (UGs) is developed, integrating a spring accumulator, inflatable accumulator, and silicon oil bladder to address buoyancy variations.These passive elements are meant to compensate for the buoyancy change in the vehicle due to the balance between hull deformation and seawater density change with depth.This scheme increases efficiency because the active buoyancy change device does not have to spend energy compensating for the aforementioned buoyancy change.A mathematical model of the system is developed, and the design parameters are optimized using a genetic algorithm-tabu search (GA-TS) method.The performance of the HPBCS is evaluated through a case study on the Petrel-II glider, showing significant improvements in motion stability and energy efficiency compared to previous single-device methods.Although very valuable, this study is limited to energy savings due to the passive elements, not considering any gains that might come from using different control strategies in the active buoyancy changes.In another study [3], a high-accuracy buoyancy-actuated system (HA-BAS) for multimodal underwater vehicles (MUVs) is designed to optimize energy use and to provide precise depth control in deep-sea environments.The buoyancy module used in [3] is hydraulic and three different pump actuation schemes are explored.Motion simulations and sea trials confirm that two of those pump actuation schemes enhance energy savings.Although providing valuable results, the pump actuation schemes proposed in [3] appear to be based in an ad hoc procedure which might not be easily applied to other type of buoyancy control systems.
Another line of research within the topic of energy reduction is focused on the balance that must be struck between the high control effort needed for accurate depth control and the lower control effort required to conserve energy.This balance has been explored extensively in the literature.
In [19], a nonlinear backstepping depth controller for a hybrid (propeller-and buoyancydriven) AUV is developed.The controller therein developed considers that a torpedo-shaped AUV has a low inertia and drag in the roll dynamic and can thus easily experience oscillatory behavior due to propeller torque or unknown disturbances.Simulation results show good control performance in both constant and variable depth references.Unfortunately, nothing is reported about energy consumption obtained with the proposed control law.
When considering the effect that the control law has on the energy consumption of the vertical motion of the vehicle, study [20] shows an interesting control approach simulating a human intelligent control (HSIC) and S-plane control.The controller is designed for the vertical motion of a long-range AUV (LRAUV).Although the LRAUV considered in [20] has several actuators (a stern thruster, shifting mass mechanism, a buoyancy adjustment mechanism, and a set of horizontal and vertical rudders), the study developed in [20] is focused only on two of those actuators: the horizontal rudders and the shifting mass mechanism.Both actuators change the value of the pitch angle, thereby affecting the depth of the vehicle which is assumed to travel at a constant horizontal speed.The controller developed in [20] includes an allocation strategy that dynamically decides which of the actuators is put in action, and a hybrid controller that switches the control structure according to the state of the system.For large depth errors, a bang-bang controller is used.When the depth error and the pitch angle are both small, the S-plane modal control [21] is adopted and finally, when depth error and pitch angle are within a certain range, the control action is kept constant.It is shown, through simulations, that by using the proposed control strategy, the energy consumption of each actuator is smaller than the energy consumption of the corresponding actuator under a single control mode.Unfortunately, no comparison with other control strategies is made in [20], so the actual energy gains obtained are unclear.
The authors of this study have previously conducted research analyzing the impact of various control laws on the energy consumption for depth control in underwater devices.In [2], the energy consumption of two actuation mechanisms for underwater platforms: variable buoyancy and propeller-actuated devices, was examined.Using a previously developed prototype, a detailed dynamic model was created to estimate energy consumption for both systems, incorporating nonlinearities such as saturations, sensor quantization, and actuator brake models.Several PID-based controllers were designed and tested through simulations.The results showed that variable buoyancy systems can drastically reduce energy consumption-by approximately 60%-for hovering operations compared to propeller-driven systems.However, this energy efficiency comes with a trade-off of greater depth control error, though variable buoyancy systems in [2] achieved faster settling times.
This paper expands the work in [2] with the following contributions: (i) it expands the model in [2] by considering the deformation of the vehicle hull with pressure; (ii) it provides a simple procedure, based on finite element analysis, to estimate such deformation; (iii) it presents a new buoyancy disturbance observer; (iv) it presents a new variable structure controller, developed for energy savings; and (v) it demonstrates that the proposed controller can further reduce energy consumption compared to the controllers presented in [2].
This paper is organized as follows: Section 2 presents the vehicle prototype and its model, including a detailed description of how the parameter expressing the vehicle deformation with pressure was estimated.Section 3 presents a new buoyancy disturbance observer, and Section 4, a new variable structure controller developed for low energy consumption.Section 5 presents several mission simulations, using both the controller de- veloped in Section 4 and previous controllers developed in the literature.Finally, Section 6 draws the main conclusions of this work.

Prototype Description
To investigate the potential benefits of a buoyancy-driven mechanism on energy efficiency, the authors developed a prototype featuring a VBM, as previously described in the literature [22] and presented in Figure 1.It should be noted that in standard operation, the prototype works with the main axis in the vertical direction.This prototype includes a main control unit (MCU-segment 1) for autonomous operation and is tailored for integration into autonomous underwater vehicles (AUVs) using modular building blocks or, alternatively, for standalone use as a buoy for vertical profiling or hovering tasks.In [2], a propeller module (PM-segment 4) was integrated into the prototype to allow for the comparison of energy consumption in real-world applications.The prototype also includes an intermediate section for floatation foam (segment 2) and the VBM (segment 3).Its overall dimensions are 1616 mm in length, 200 mm for the outer diameter, 180 mm for the inner diameter, and its dry weight is 36 kg.Segment 1 has a length of 400 mm, segment 2 has a length of 200 mm, segment 3 has a length of 616 mm, and segment 4 has a length of 300 mm.The black end cap seen at the left side of Figure 1 has a length of 100 mm.Buoyancy adjustment is achieved by pumping seawater via a diaphragm-sealed piston mechanism [22] driven by an electric motor coupled to a mechanical transmission and spindle.This VBM enables a total volume change of approximately ±350 cm 3 and can operate down to a depth of 100 m.This prototype, which will be used in this work without any hardware changes relative to the one presented in [2], incorporates two sensors which are standard for this type of module: (i) a pressure sensor, which is used to determine the depth of the vehicle and (ii) a position sensor to measure the position of the diaphragmsealed piston mechanism.This second sensor allows the estimation of the buoyancy force as it enables the measurement of the variable volume.Only these two sensors are required for the controllers developed in subsequent sections.
observer, and Section 4, a new variable structure controller developed for low energy consumption.Section 5 presents several mission simulations, using both the controller developed in Section 4 and previous controllers developed in the literature.Finally Section 6 draws the main conclusions of this work.

Prototype Description
To investigate the potential benefits of a buoyancy-driven mechanism on energy efficiency, the authors developed a prototype featuring a VBM, as previously described in the literature [22] and presented in Figure 1.It should be noted that in standard operation the prototype works with the main axis in the vertical direction.This prototype includes a main control unit (MCU-segment 1) for autonomous operation and is tailored for integration into autonomous underwater vehicles (AUVs) using modular building blocks or, alternatively, for standalone use as a buoy for vertical profiling or hovering tasks.In [2], a propeller module (PM-segment 4) was integrated into the prototype to allow for the comparison of energy consumption in real-world applications.The prototype also includes an intermediate section for floatation foam (segment 2) and the VBM (segment 3).Its overall dimensions are 1616 mm in length, 200 mm for the outer diameter, 180 mm for the inner diameter, and its dry weight is 36 kg.Segment 1 has a length of 400 mm segment 2 has a length of 200 mm, segment 3 has a length of 616 mm, and segment 4 has a length of 300 mm.The black end cap seen at the left side of Figure 1 has a length of 100 mm.Buoyancy adjustment is achieved by pumping seawater via a diaphragm-sealed piston mechanism [22] driven by an electric motor coupled to a mechanical transmission and spindle.This VBM enables a total volume change of approximately ±350 cm 3 and can operate down to a depth of 100 m.This prototype, which will be used in this work without any hardware changes relative to the one presented in [2], incorporates two sensors which are standard for this type of module: (i) a pressure sensor, which is used to determine the depth of the vehicle and (ii) a position sensor to measure the position of the diaphragmsealed piston mechanism.This second sensor allows the estimation of the buoyancy force as it enables the measurement of the variable volume.Only these two sensors are required for the controllers developed in subsequent sections.The main focus of this work is to introduce a novel variable structure controller to reduce the energy consumption of the prototype, along with a disturbance observer to assist in controlling tasks by estimating the disturbances affecting the prototype.To this end, in Section 2.2. the model of the prototype required for this purpose is presented.

Prototype Model
The prototype model comprised the combination of the VBM model and the vertical motion model, as developed in [23], and is presented in Figure 2. The main focus of this work is to introduce a novel variable structure controller to reduce the energy consumption of the prototype, along with a disturbance observer to assist in controlling tasks by estimating the disturbances affecting the prototype.To this end, in Section 2.2. the model of the prototype required for this purpose is presented.

Prototype Model
The prototype model comprised the combination of the VBM model and the vertical motion model, as developed in [23], and is presented in Figure 2.
unit volume (assumed to be constant,  = 1000 kg/m ) and the acceleration of gravity ( = 9.8 ms ), respectively.Parameter  expresses the loss of volume per meter of depth due to hull deformation and will be determined in Section 2.3.This loss of volume will result in a disturbance force  _ .

Estimation of Parameter 𝜓
During an underwater mission, an AUV experiences the effects of pressure due to the water column height.This effect happens in all submerged bodies, but it is especially relevant in sealed hollow bodies with a compressible fluid inside.In fact, when a solid body is submerged, the water pressure is evenly distributed along the whole body's surface area, leading to a pressure equilibrium state in which the main load is one of compression.Though this effect may be relevant for the submerged volume of a solid body, it is reasonably safe to assume that its impact is minimal when compared to a sealed volume with a compressible fluid on the inside, such as air.In this situation, the equilibrium of pressure acting on the outside and inside of the hollow body may not be achieved without a significant amount of deformation, which is heavily dependent on the body material stiffness.Considering this, the deformation of the hull of an AUV can have In Figure 2, K 1 and T 1 are the actuator velocity dynamics steady-state gain and time constant, respectively, that were experimentally identified in [22], and A is the area of the piston in contact with sea water.When the control action U = 0, the brake is engaged, so the equivalent depth voltage U Z is zero.If U ̸ = 0, the brake is disengaged, so U Z = k z /k u Z to account for the influence of pressure due to depth Z on the actuator motion.The ratio k z /k u was determined in [2] and models the torque caused on the motor due to the force exerted by the outside pressure, as well as the motor stall torque at the applied voltage level.F net is the balance between the actual VBM buoyancy force F VBM (generated by the VBM buoyancy volume V VBM ), and the disturbances acting on the prototype F dist , whether internal (F dist_int ) or external (F dist_ext ).K 2 and T 2 are the depth dynamics parameters identified in [23].The parameters ρ and g are the water mass per unit volume (assumed to be constant, ρ = 1000 kg/m 3 ) and the acceleration of gravity (g = 9.8 ms −2 ), respectively.Parameter ψ expresses the loss of volume per meter of depth due to hull deformation and will be determined in Section 2.3.This loss of volume will result in a disturbance force F dist_int .

Estimation of Parameter ψ
During an underwater mission, an AUV experiences the effects of pressure due to the water column height.This effect happens in all submerged bodies, but it is especially relevant in sealed hollow bodies with a compressible fluid inside.In fact, when a solid body is submerged, the water pressure is evenly distributed along the whole body's surface area, leading to a pressure equilibrium state in which the main load is one of compression.Though this effect may be relevant for the submerged volume of a solid body, it is reasonably safe to assume that its impact is minimal when compared to a sealed volume with a compressible fluid on the inside, such as air.In this situation, the equilibrium of pressure acting on the outside and inside of the hollow body may not be achieved without a significant amount of deformation, which is heavily dependent on the body material stiffness.Considering this, the deformation of the hull of an AUV can have serious consequences on its buoyancy, particularly in dry segments of the vehicle.To assess the effect of pressure on the volume occupied by the dry segments of the prototype, simplified 3D models for the VBM and MCU segments were created using Solidworks (release 2020).The VBM model consists of a hollow cylinder of uniform material.Simulations using the finite elements method were then run on the VBM 3D model considering the following conditions:
Each external edge of the cylinder is not allowed to rotate about the cylinder longitudinal axis; 2.
One of the external edges is not allowed to move along the longitudinal axis; • External loads: pressure along the whole outside surface area of the cylinder.
The 3D model with the applied fixtures and loads is presented in Figure 3.
assess the effect of pressure on the volume occupied by the dry se simplified 3D models for the VBM and MCU segments were cr (release 2020).The VBM model consists of a hollow cylinde Simulations using the finite elements method were then run considering the following conditions: • Material: POM acetal copolymer; • Fixtures: 1.Each external edge of the cylinder is not allowed to r longitudinal axis; 2. One of the external edges is not allowed to move along

•
External loads: pressure along the whole outside surface are The 3D model with the applied fixtures and loads is present In Figure 3, the green arrows represent fixtures, and the red as described in the conditions above.The simulation was run levels (2.5, 5, 7.5, 10 bar) to simulate the effect of water depth presents the deformed shape of the 3D model subjected to a 10-b figure, the deformations are represented with an amplification o visible.The von Mises stress is represented with a jet colormap.T of the model was then imported back to SolidWorks as a new p volume was calculated ( ).Table 1 presents the results of the under different pressure levels and the difference to the ambient In Figure 3, the green arrows represent fixtures, and the red ones represent the loads, as described in the conditions above.The simulation was run with different pressure levels (2.5, 5, 7.5, 10 bar) to simulate the effect of water depths up to 100 m. Figure 4 presents the deformed shape of the 3D model subjected to a 10-bar pressure load.In this figure, the deformations are represented with an amplification of 40 times in order to be visible.The von Mises stress is represented with a jet colormap.The true deformed shape of the model was then imported back to SolidWorks as a new part, and the new closed volume was calculated (Vol p ). Table 1 presents the results of the VBM 3D model volume under different pressure levels and the difference to the ambient pressure volume Vol 0 .
The results presented in Table 1 show that the loss of volume due to pressure is approximately linear in the simulated pressure range.This validates the use of a parameter ψ to express the loss of volume per meter depth.An identical process was followed with the MCU.In the case of the MCU 3D model, however, an assembly of the acrylic cylinder body and two aluminum flanges was considered.The MCU 3D model is presented in Figure 5. presents the deformed shape of the 3D model subjected to a 10-bar pressure load.In this figure, the deformations are represented with an amplification of 40 times in order to be visible.The von Mises stress is represented with a jet colormap.The true deformed shape of the model was then imported back to SolidWorks as a new part, and the new closed volume was calculated ( ).Table 1 presents the results of the VBM 3D model volume under different pressure levels and the difference to the ambient pressure volume  .The results presented in Table 1 show that the loss of volume due to pre approximately linear in the simulated pressure range.This validates the u parameter  to express the loss of volume per meter depth.An identical pro followed with the MCU.In the case of the MCU 3D model, however, an assemb acrylic cylinder body and two aluminum flanges was considered.The MCU 3D presented in Figure 5.The deformed shape (amplified by a factor of 40 times) and the von Mises stress on the 3D model subjected to a 10-bar pressure load is shown in Figure 6.As performed with the VBM model, the true deformed shape was reimported to SolidWorks as a new part, and the new volume calculated.It was found that for a 10-bar pressure, the MCU loss of volume would be about 7 cm 3 .Finally, the parameter ψ, the loss of volume per meter depth, was calculated as the sum of the volume loss of the VBM ( Vol 0 − Vol p VBM ) and the volume loss of the MCU Vol 0 − Vol p MCU , at 100 m (10 bar of pressure):

Buoyancy Disturbance Observer
The buoyancy disturbance observer is based on the block diagram in Figure 7.The vehicle model in the vertical direction, shown in the upper part of Figure 7, is a second order system with parameters  ,  , and . and  have been identified in [23] and  has been estimated in Section 2 of this paper.The system model used for the observer is also second-order, but it will be assumed in the following that its parameters might not be fully known.This allows the study of the effect of this uncertainty in the observer performance.It will be considered that  is an estimate of  , while  is an estimate of  , such that: and where  is the error in the  parameter and  is the error in the  parameter.
For the purpose of the disturbance observer synthesis, the disturbance force ( ) that will be observed ( ) is considered to be composed of two parts: an external one ( _ ) caused, for instance, by dropping or collecting a weight, and an internal one ( _ ), caused by the vehicle deformation due to pressure.

Buoyancy Disturbance Observer
The buoyancy disturbance observer is based on the block diagram in Figure 7.The vehicle model in the vertical direction, shown in the upper part of Figure 7, is a second order system with parameters K 2 , T 2 , and ψ.K 2 and T 2 have been identified in [23] and ψ has been estimated in Section 2 of this paper.The system model used for the observer is also second-order, but it will be assumed in the following that its parameters might not be fully known.This allows the study of the effect of this uncertainty in the observer performance.It will be considered that K 22 is an estimate of K 2 , while T 22 is an estimate of T 2 , such that: and where ξ K is the error in the K 22 parameter and ξ T is the error in the T 22 parameter.
For the purpose of the disturbance observer synthesis, the disturbance force (F dist ) that will be observed ( Fdist ) is considered to be composed of two parts: an external one (F dist_ext ) caused, for instance, by dropping or collecting a weight, and an internal one (F dist_int ), caused by the vehicle deformation due to pressure.
Based on Figure 7, the error E OZ (s) between the vehicle depth Z(s) and the observer depth estimation Ẑ(s) can be written as The disturbance estimation error E dist (s) can be written as Replacing (4) in ( 5) leads to (6), the transfer functions between F VBM (s) and the disturbance observer error E dist (s) and between the disturbance force F dist (s) and E dist (s) can be written as The disturbance estimation error  () can be written as Replacing (4) in ( 5) leads to (6), the transfer functions between  () and the disturbance observer error  () and between the disturbance force  () and  () can be written as The dynamics between the disturbance  () and the observer error  () can be rewritten as Analyzing Equation (7), it can be concluded that the observer dynamics cannot be fully imposed as they present third-order dynamics, and the observer only has two adjustable parameters.The transfer function (7) can nevertheless be used to adequately tune the observer dynamic response since its poles are given by A particular case of Equation ( 7) occurs when there is not any error in the observer parameters  and  .In this case, the disturbance of observer dynamics does not depend on the force  () , and the transfer function between  () and  () simplifies to Equation (9): The dynamics between the disturbance F dist (s) and the observer error E dist (s) can be rewritten as Analyzing Equation ( 7), it can be concluded that the observer dynamics cannot be fully imposed as they present third-order dynamics, and the observer only has two adjustable parameters.The transfer function (7) can nevertheless be used to adequately tune the observer dynamic response since its poles are given by A particular case of Equation ( 7) occurs when there is not any error in the observer parameters K 22 and T 22 .In this case, the disturbance of observer dynamics does not depend on the force F VBM (s), and the transfer function between F dist (s) and E dist (s) simplifies to Equation ( 9): The poles in this situation are given by As seen in (10), when the observer parameters exactly match the system parameters, the observer poles are the complex poles of (8), which can be placed according to the values of K O and K OD .The upper pole of (8) appears only when the observer parameters are different from the system parameters.
It is now possible to determine the steady-state error of the observer for different F dist and F VBM inputs.Considering a step of amplitude A 1 for F dist and a step of amplitude A 2 for F VBM , the steady-state error can be calculated by the summation of the errors due to each input, leading to The result of Equation (11) shows that the steady-state error of a constant input becomes smaller as the error in the K 2 parameter becomes smaller.If this error becomes zero, then the observer can predict exactly a constant disturbance at a steady state.Note also that the result of Equation (11) does not depend on ξ T since the time constant does not affect the steady-state error.
Considering now a ramp of slope R 1 for F dist and of slope R 2 for F VBM , the steady-state error can be calculated by the summation of the errors due to each input, leading to Analyzing Equation ( 12), the error of a ramp disturbance is infinite, i.e., the observer cannot adequately track a ramp disturbance.In the particular case when there is not any error in the observer parameters K 22 and T 22 , it is clear from (9) that F VBM has no influence on the observer tracking capability.In this case, the steady-state error of the observer to an F dist ramp input with a slope R 3 can be obtained using Equation (9) as An analysis of Equation (13) shows that when the model parameters are known, the observer is able to track a ramp input, although with an error which decreases with the increase in K O .
According to Equation (8), it is not possible to fully impose a dynamic for the observer, since one of its poles is exclusively a function of the time constant T 2 and therefore, its position is set.As for the other two poles, they will depend on the values chosen for K O and K OD .A common strategy when performing pole placement is to choose all real poles to minimize undesired oscillatory behavior.Therefore, the radicand in Equation (8) must be equal to or greater than zero: Solving for K OD to determine the roots of the polynomial: Since, from Equation (8), K OD ≥ −ρg/K 22 to ensure stability of the observer, there is only one possible root: For the inequality in ( 14) to be ensured, and since K 22 2 > 0, the following condition must be met: Sensors 2024, 24, 5771 11 of 19 Solving for K O in Equation ( 17) leads to An analysis of expressions ( 16) and (18) shows that when the observer gain K O takes the minimum value of Equation ( 18), (K O = ρg 4T 22 K 22 ), then it follows by Equation ( 16) that K OD = 0. Replacing these two values in Equation (8), the critical damping situation occurs, and the poles overlap at s = −1/(2T 22 ).

Controllers
The depth control of the prototype can generically be performed using several different approaches.In [2], the authors proposed two control architectures based on PID family controllers: a single-depth controller for the PM and a cascaded control strategy consisting of a depth controller and a volume controller for the VBM.Among the controllers tested in [2], the structures that led to the smallest energy consumption were selected for datum comparison in this work: Figure 8 presents the block diagram for structures 3b and 7b.In Figure 8, Z re f is the target depth reference and E Z is the depth error, whereas V re f is the target volume reference and E V is the volume error.Figure 9 shows the block diagram for structure 11b, where I is the control action and F PM is the PM buoyancy force.

Controllers
The depth control of the prototype can generically be performed using several different approaches.In [2], the authors proposed two control architectures based on PID family controllers: a single-depth controller for the PM and a cascaded control strategy consisting of a depth controller and a volume controller for the VBM.Among the controllers tested in [2], the structures that led to the smallest energy consumption were selected for datum comparison in this work: Figure 8 presents the block diagram for structures 3b and 7b.In Figure 8,  is the target depth reference and  is the depth error, whereas  is the target volume reference and  is the volume error.Figure 9 shows the block diagram for structure 11b, where  is the control action and  is the PM buoyancy force.The full detail of each controller structure can be found in [2].This paper introduces a new variable structure controller which will be detailed in the next section.

Controllers
The depth control of the prototype can generically be performed using several different approaches.In [2], the authors proposed two control architectures based on PID family controllers: a single-depth controller for the PM and a cascaded control strategy consisting of a depth controller and a volume controller for the VBM.Among the controllers tested in [2], the structures that led to the smallest energy consumption were selected for datum comparison in this work: Figure 8 presents the block diagram for structures 3b and 7b.In Figure 8,  is the target depth reference and  is the depth error, whereas  is the target volume reference and  is the volume error.Figure 9 shows the block diagram for structure 11b, where  is the control action and  is the PM buoyancy force.The full detail of each controller structure can be found in [2].This paper introduces a new variable structure controller which will be detailed in the next section.The full detail of each controller structure can be found in [2].This paper introduces a new variable structure controller which will be detailed in the next section.

Proposed Variable Structure Controller
To minimize the energy consumption of the VBM, a novel variable structure controller (VSC) was devised, as presented in Figure 10.The underlying idea behind the VSC is to make use of the systems' inherent dynamics, whenever it is favorable to the prototype depth control.In other words, switching on the VBM is only required whenever the prototype is moving away from the depth target reference.When the prototype is moving in the direction that minimizes the depth error, it should be switched off.In this way, the VBM is expected to be switched off for a considerable amount of time throughout the mission, therefore saving energy.
To minimize the energy consumption of the VBM, a novel variable structure controller (VSC) was devised, as presented in Figure 10.The underlying idea behind the VSC is to make use of the systems' inherent dynamics, whenever it is favorable to the prototype depth control.In other words, switching on the VBM is only required whenever the prototype is moving away from the depth target reference.When the prototype is moving in the direction that minimizes the depth error, it should be switched off.In this way, the VBM is expected to be switched off for a considerable amount of time throughout the mission, therefore saving energy.The switch block, detailed below, is responsible for deciding whether it is necessary to switch on the VBM or not. Figure 11 presents a schematic to explain the switch block decision process.
In the shaded regions of the schematic in Figure 11, the VBM is switched off, and in the remaining regions, it is switched on.If there are no disturbances present,  is sufficient to characterize the prototype depth since it is the only force that can change the prototype vertical motion equilibrium.In real-world conditions, there are countless possible disturbances that may occur during an underwater mission, such as load lifting, buoyancy trimming errors, and hull deformation, among others.This means that it is possible for a situation to occur where  should lead to an expected output in terms of prototype acceleration, but the existence of  does not lead to the expected prototype movement.To take into account disturbances, the vertical axis of Figure 11's schematic represents the net force  ( =  +  ).Parameter  is a depth deadband, meaning there is an admissible depth error in which the controller is switched off (| |  ⟹  = 0).Considering the system starts in region α, the depth error  and net force  are positive ( 0,  0).In this situation, the target depth reference is deeper than the current depth (  ) but the positive buoyancy would lead to a decrease in depth, which would increase  .For this reason, the VBM should be switched on in this case, so that  = −( ) •  .Since  will directly act on the derivative of  and so, on the derivative of  (please check Figure 2), the value of  will eventually become negative and the depth error will eventually decrease.A similar line of thought can be followed for region δ of Figure 11, while keeping in mind that in this region, the depth error and the net force are both negative.In regions β and γ, the depth error and the net volume have opposite signs; therefore, despite the VBM being turned off, the natural dynamics of the system will inevitably lower the error, leading the system to region Ω.The system cannot remain at a steady state inside region Ω, so it will eventually move to the other regions.From the description above, it is possible to infer that parameter  has a direct influence on the steady-state error.The switch block, detailed below, is responsible for deciding whether it is necessary to switch on the VBM or not. Figure 11 presents a schematic to explain the switch block decision process.Table 2 presents the different switch block decisions and control actions  according to each region in the schematic of Figure 11.

Simulation Results
In this section, the VSC performance presented in Section 4 will be assessed.To that end, the performance of the VSC will be tested against (i) the PID-type controllers presented in [2] leading to the lowest energy consumption of the VBM (controllers 3b and 7b of [2]) and (ii) the PID-type controller leading to the lowest energy consumption of the prototype actuated by propellers (controller 11b of [2]).To simulate the VSC, it is also necessary to determine adequate gains for the disturbance observer.In the simulations, it In the shaded regions of the schematic in Figure 11, the VBM is switched off, and in the remaining regions, it is switched on.If there are no disturbances present, F VBM is sufficient to characterize the prototype depth acceleration, since it is the only force that can change the prototype vertical motion equilibrium.In real-world conditions, there are countless possible disturbances that may occur during an underwater mission, such as load lifting, buoyancy trimming errors, and hull deformation, among others.This means that it is possible for a situation to occur where F VBM should lead to an expected output in terms of prototype acceleration, but the existence of F dist does not lead to the expected prototype movement.To take into account disturbances, the vertical axis of Figure 11's schematic represents the net force F net (F net = F VBM + F dist ).Parameter a is a depth deadband, meaning there is an admissible depth error in which the controller is switched off ( |E Z | < a =⇒ U = 0 ).Considering the system starts in region α, the depth error E Z and net force F net are positive (E Z > 0, F net > 0).In this situation, the target depth reference is deeper than the current depth (Z re f > Z) but the positive buoyancy would lead to a decrease in depth, which would increase e Z .For this reason, the VBM should be switched on in this case, so that U = −sign(E Z )•U switch .Since U will directly act on the derivative of V VBM and so, on the derivative of F VBM (please check Figure 2), the value of F net will eventually become negative and the depth error will eventually decrease.A similar line of thought can be followed for region δ of Figure 11, while keeping in mind that in this region, the depth error and the net force are both negative.In regions β and γ, the depth error and the net volume have opposite signs; therefore, despite the VBM being turned off, the natural dynamics of the system will inevitably lower the error, leading the system to region Ω.The system cannot remain at a steady state inside region Ω, so it will eventually move to the other regions.From the description above, it is possible to infer that parameter a has a direct influence on the steady-state error.
Table 2 presents the different switch block decisions and control actions U according to each region in the schematic of Figure 11.

Simulation Results
In this section, the VSC performance presented in Section 4 will be assessed.To that end, the performance of the VSC will be tested against (i) the PID-type controllers presented in [2] leading to the lowest energy consumption of the VBM (controllers 3b and 7b of [2]) and (ii) the PID-type controller leading to the lowest energy consumption of the prototype actuated by propellers (controller 11b of [2]).To simulate the VSC, it is also necessary to determine adequate gains for the disturbance observer.In the simulations, it was assumed there was no error in the estimation of parameters K 2 and T 2 .The observer gains were chosen so that all poles are real.Table 3 presents the values chosen for the simulation.Table 4 lists the observer pole locations.

Gain
Value Unit The controller gains for the controllers presented in [2] are the same as the ones used in [2].Since those controller gains were chosen to minimize the energy consumption of the prototype, the same premise was considered for the VSC.As such, the parameters for the VSC were experimentally tuned to keep a low energy consumption while maintaining an acceptable time response.The chosen parameters were U switch = 8 V and a = 0.1 m.
The simulation trials were run in Matlab Simulink (R2019a) with the same depth target reference signal and the same external buoyancy force disturbances considered in [2].Three disturbance scenarios were applied: in type 1 trials, no disturbance was considered.In type 2 trials, a constant disturbance of 1.75 N was added.In type 3 tests, the disturbance increased in discrete 0.35 N steps: starting from 1.75 N at time t = 0 s, it increased by 0.35 N at t = 1800 s, then increased by 0.35 N every 1800 s until it reached 2.8 N. At t = 12,600 s, it then decreased by 0.35 N every 1800 s until it reached 2.1 N at t = 16,200 s.This disturbance value was then held until t = 18,000 s.In this work, the internal disturbance of the hull deformation was added with the parameter ψ = 3.4 × 10 −7 m 3 /m, as estimated in Section 2.3.The energy consumption results for each controller are presented in Table 5.  show the simulated prototype depth during the three tests with each controller.the prototype, the same premise was considered for the VSC.As such, the parameters for the VSC were experimentally tuned to keep a low energy consumption while maintaining an acceptable time response.The chosen parameters were  = 8 V and  = 0.1 m.The simulation trials were run in Matlab Simulink (R2019a) with the same depth target reference signal and the same external buoyancy force disturbances considered in [2].Three disturbance scenarios were applied: in type 1 trials, no disturbance was considered.In type 2 trials, a constant disturbance of 1.75 N was added.In type 3 tests, the disturbance increased in discrete 0.35 N steps: starting from 1.75 N at time t = 0 s, it increased by 0.35 N at t = 1800 s, then increased by 0.35 N every 1800 s until it reached 2.8 N. At t = 12,600 s, it then decreased by 0.35 N every 1800 s until it reached 2.1 N at t = 16,200 s.This disturbance value was then held until t = 18,000 s.In this work, the internal disturbance of the hull deformation was added with the parameter  = 3.4 × 10 −7 m 3 /m, as estimated in Section 2.3.The energy consumption results for each controller are presented in Table 5. Figures 12-14 show the simulated prototype depth during the three tests with each controller.Figure 12.Simulated prototype depth for each controller in type 1 trials: zref is the reference position and zj is the depth obtained with controller j in Table 5. Figure 12.Simulated prototype depth for each controller in type 1 trials: z ref is the reference position and z j is the depth obtained with controller j in Table 5.
Sensors 2024, 24, x FOR PEER REVIEW 15 of 19 . Simulated prototype depth for each controller in type 2 trials: zref is the reference position and zj is the depth obtained with controller j in Table 5. Figure 13.Simulated prototype depth for each controller in type 2 trials: z ref is the reference position and z j is the depth obtained with controller j in Table 5.
Figure 13.Simulated prototype depth for each controller in type 2 trials: zref is the reference position and zj is the depth obtained with controller j in Table 5.
Figure 14.Simulated prototype depth for each controller in type 3 trials: zref is the reference position and zj is the depth obtained with controller j in Table 5.
To evaluate the performance in step i (i = 1:5) of the controller j (j = 1:4), the performance metrics outlined in [2] and listed below were used.
1. Overshoot:  , .2. Maximum error magnitude  , : this metric is defined as the maximum error after the third inflection point of the prototype for each step i and controller j. 3. Absolute settling time  , : the time required for the prototype to reach and stay within its  , band. 4. Relative settling time  , : the time required for the prototype to reach the widest band across all control structures (  , ). 5. Maximum error after new external disturbance:  , .
Table 6 presents the average values of these metrics for all steps i and each controller j.  5.
To evaluate the performance in step i (i = 1:5) of the controller j (j = 1:4), the performance metrics outlined in [2] and listed below were used.
Maximum error magnitude Ezm i,j : this metric is defined as the maximum error after the third inflection point of the prototype for each step i and controller j.

3.
Absolute settling time tssa i,j : the time required for the prototype to reach and stay within its Ezm i,j band.4.
Relative settling time tssr i,j : the time required for the prototype to reach the widest band across all control structures (max j Ezm i,j ). 5.
Maximum error after new external disturbance: Ezd i,j .
Table 6 presents the average values of these metrics for all steps i and each controller j.Considering the results presented in Tables 5 and 6, it is clear that when considering the average energy consumption for all trials, the prototype actuated by controller 3 leads to the lowest values, requiring 22% less energy when compared to the next best-performing controller (controller 2).However, it yields the worst control performance, since on average, it leads to the highest overshoot, error magnitude, and relative settling time.In performance terms, it could be argued that controller 4 is the best, given the minimal average overshoot and error magnitude.Nevertheless, its energy consumption is the highest in the type 2 and 3 trials and it is very slow, with the highest average absolute settling time.This high average absolute settling time is also relevant to understand this combination's high error after a new external disturbance, since it takes more time for the controller to react to the disturbance, during which the prototype depth error is increasing.The most balanced combination may be controller 1, since it leads to a low average overshoot, error magnitude, error after new external disturbance, and acceptable average absolute and relative settling times.In comparison, with a similar energy consumption, controller 2 leads to a bigger average overshoot and error magnitude.However, it is also the solution with the average absolute and relative settling times.
Finally, the effectiveness of each controller in ensuring the sensor platform's viability for underwater acoustic measurements was evaluated.This evaluation involved counting the duration each controller remained switched off for at least 60 s.Table 7 presents the accumulated total time for each trial and Table 8 shows the percentages of the total running time for each trial.It was observed that the controller proposed in this work remained switched off for at least 60 s for more than 75% of each trial duration.This indicates that the proposed controller is highly suitable for applications requiring a minimal acoustic signature.Notably, these results are comparable to those obtained for gliders (approximately 87% off time in [24]).However, while the trials in [24] were based on sawtooth-profiled trajectories, this study employed hovering operations, which pose greater challenges in keeping the actuators switched off, especially in the presence of disturbances.Finally, it should be remarked that (i) controller 1 presents much less switched-off time than the proposed controller and (ii) controllers 2 and 4 could not be applied whenever the sensor platform was used for acoustic measurements, as they never switch off.

Conclusions
This paper focused on the development of a novel variable structure controller for depth control of an autonomous underwater device prototype using a variable buoyancy module.The parameters of the prototype linear model were presented and the prototype hull deformation due to pressure was estimated using a finite element method.The presence of possible disturbances acting on the prototype, such as hull deformation, led to the need to estimate said disturbances.For this reason, a disturbance observer was developed.The estimation error was analyzed according to the input disturbance and system parameter uncertainty.Additionally, the location of the observer poles was determined.
Considering the prototype linear model, a plane of state-dependent variables was defined and split into different regions.The proposed variable structure controller is switched on when the system's natural dynamics are unfavorable (when the vehicle will tend to move away from the desired reference) and is switched off when the system's natural dynamics are favorable (when the vehicle will tend to move toward the desired reference).A detailed analysis of the controller was performed to determine the necessary switching control action that guarantees that the system is driven to the desired reference.
The proposed controller was tested through a simulation along with PID-based controllers to evaluate their dynamic response and energy consumption across different operation conditions.Both VBM-and propeller-actuated vehicles were considered.The findings indicate that the proposed controller requires 22% less average energy consumption compared to the next best-performing controller (a PID-based controller used on the VBMactuated vehicle).However, this energy reduction comes with a significant cost in terms of control performance in most of the proposed metrics.It can nevertheless be of important use whenever the mission only requires coarse-performance hovering tasks.An interesting feature of the proposed controller is that it is switched off for long periods of time, thereby being very silent, which is useful for underwater acoustic monitoring.
Future work will focus on sea trials of the proposed controller.Also, the development of controllers which can further reduce the time during which the system is switched on to reduce the acoustic trace of the sensor platform, will be pursued.

Figure 1 .
Figure 1.A picture of the prototype.

Figure 1 .
Figure 1.A picture of the prototype.

Figure 3 .
Figure 3. VBM 3D model with fixtures and loads for simulation.

Figure 3 .
Figure 3. VBM 3D model with fixtures and loads for simulation.

Figure 4 .
Figure 4. Deformed shape of the VBM module and von Mises stress.

Figure 4 .
Figure 4. Deformed shape of the VBM module and von Mises stress.

Figure 5 .
Figure 5. MCU 3D model with fixtures and loads for simulation.

Figure 5 .
Figure 5. MCU 3D model with fixtures and loads for simulation.

Figure 6 .
Figure 6.Deformed shape of the MCU model and von Mises stress.

Figure 6 .
Figure 6.Deformed shape of the MCU model and von Mises stress.

Figure 7 .
Figure 7. Block diagram of the vehicle model and observer.

•
Structure 3b: a PID depth controller and a PI volume controller (VBM); • Structure 7b: an I-PID depth controller and a PI volume controller (VBM); • Structure 11b: an I-PD depth controller (PM).

•
Structure 3b: a PID depth controller and a PI volume controller (VBM); • Structure 7b: an I-PID depth controller and a PI volume controller (VBM); • Structure 11b: an I-PD depth controller (PM).

•
Structure 3b: a PID depth controller and a PI volume controller (VBM); • Structure 7b: an I-PID depth controller and a PI volume controller (VBM); • Structure 11b: an I-PD depth controller (PM).

Figure 11 .
Figure 11.Schematic of the switch block decision process.

Figure 14 .
Figure 14.Simulated prototype depth for each controller in type 3 trials: z ref is the reference position and z j is the depth obtained with controller j in Table5.

Table 1 .
Volume of the VBM 3D model under different pressure levels.

Table 1 .
Volume of the VBM 3D model under different pressure levels.

Table 2 .
Switch block decision table.

Table 2 .
Switch block decision table.

Table 3 .
Disturbance observer minimum and chosen gains.

Table 5 .
Energy consumption results for the simulated trials.

Table 5 .
Energy consumption results for the simulated trials.

Table 6 .
Performance metrics for each control structure.

Table 7 .
Switched-off time for each trial.

Table 8 .
Switched-off time for each trial (percentages of the total running time for each trial).