Dynamic Model and Control of a New Underwater Three-Degree-of-Freedom Tidal Energy Converter

There is currently a growing interest in developing devices that can be used to exploit energy from oceans. In the recent past, the search for oil and gas at ever-greater depths has led to the evolution of deviceswithwhich these resources are extracted.These devices range from those that simply rest on the seabed to those that are fully floating and anchored to it.This trend can be considered as the basis needed to understand the future evolution of devices for harnessing depth renewable resources. This paper presents a simple dynamic modeling and a nonlinear multivariable control model-based system for a new three-degree-of-freedom underwater generator with which energy from depth marine currents is harnessed when reference trajectory tracking for the emersion maneuvers needed to carry out maintenance tasks is performed. The goodness of both the model and the proposed controller has been demonstrated through the development of various simulations in the MATLAB-Simulink environment. Additionally, the validation of the control algorithms was carried out by using the dynamic model offered by the simulation environment Orcina OrcaFlex (software for the dynamic analysis for offshore marine systems) through the MATLAB-OrcaFlex interface.


Introduction
The growing interest in the exploitation of marine renewable energies began several years ago, and various devices with which to harness energy from seas have therefore been conceived or developed (see [1][2][3]), their main natural energy sources being wind, waves, and marine currents.
One of the most promising sources of marine energy is the exploitation of tidal or oceanic water flows [4][5][6][7], and the industry's effort is currently focused on the so-called first generation devices [8] (fixed to the sea bottom and suitable for sites with depths below 40 m).But there are an increasing number of second-generation devices that have been conceived to be moored to the sea bottom with an expected similar trend to that which has taken place during the development of oil and gas platforms that must access resources at increasingly greater depths [9].
This evolution has led these devices to evolve from being anchored to the seabed to being located in a floating location, and the most appropriate ways and means to perform maintenance tasks have therefore also had to evolve these kinds of devices [10][11][12][13][14][15].
The successful installation of these kinds of devices for harnessing energy from depth currents can only take place once it has been proved that they are both technically and economically feasible in comparison to other traditional energy sources.One well known way in which to reduce costs is by successively automating more tasks, thus signifying less human intervention or the possibility of using the cheapest general purpose ships rather than high cost special vessels for maintenance purposes.In [16] are studied the automatic maneuvering emersion and immersion of one of these devices.
Nowadays there exist different devices with the following main alternatives for performing maintenance tasks: (i) There is the use of a servo actuated crabbing based system to move the main generation unit from the support structure.See different devices from [17,18].(iii) There is the use of a ballast management system to generate vertical forces, thus enabling the device's emersion and immersion movements to be controlled.See [24][25][26].
A new family of generators is briefly presented in this work.They use a mooring system based on buoys and wires that allow the device to be located at the desired position on the seabed and positioned at any desired depth of the layer with almost no human resources.By simply disconnecting stern wire that joins the generator to the seabed and a proper management of the ballast water that is strategically located inside the generator, it is possible to perform closed loop emersion and immersion maneuvers with three degrees of freedom (DoFs), two orientation angles and a depth control, in such a way that the process of moving from a submerged and vertically disposed state to a floating and horizontally disposed one and vice versa can be carried out fully automatically without any kind of human intervention.This paper presents a new dynamic model for this new family of second-generation tidal energy converters [8,27,28] for simulation and control design purposes [29,30].The proposed model is very simple, fully parameterized, and easily scalable with minimum computer effort.Very simple models have already been used to control complex devices with a good performance and good correspondence among simulated and experimental signals.For example, in [31], a simple model based on only one lumped mass and a special kinematic uncoupling design was used to model a threedegree-of-freedom flexible arm with complex nonlinear and time variant dynamics.This simplified dynamic model was employed as a basis for various controllers, which were used to control this robot with excellent results [32].Another different system, in this case, a stair-climbing mobility system, was also modeled with only one lumped mass [33] and this simple dynamic model has been used as the basis for various control systems (see, a.e., [34]).
In a similar manner to that shown above, the dynamic model proposed in this paper is based on the computation of the values of only four lumped masses, which are handled solely by buoying forces applied to three star-shaped distributed equilateral torpedoes on the pod (nacelle) of the generator.The model developed exhibits time-varying, nonlinear, and strongly coupled behavior.It was, meanwhile, necessary to design the control law based on the dynamic model developed in such a way that it would have a successful closed loop behavior when the underwater three-DoF tidal energy converter performs emersion and immersion maneuvers with only passive buoyancy forces.In this work, the proposed control system is characterized by its simplicity, computational efficiency, and easy implementation in a microprocessor-/microcontroller-based system.The proposed control law is composed of the following two main terms: (i) a matrix that is responsible for decoupling the open loop dynamics from one degree of freedom with regard to the others and (ii) a feedback controller based on centrifugal and Coriolis force cancellation and a simple multivariable and diagonal PD controller.
The paper is organized as follows: after providing a brief description of the proposed family of generators in Section 2, Section 3 shows how the derivation of the proposed dynamic model of the submerged device was derived when performing coupled three-degree-of-freedom motions.Section 4 proposes the solution adopted to control the system.Section 5 presents the numerical simulations obtained to validate the proposed dynamic model and the proposed control algorithm.In particular, different smooth robotic-based trajectories were used to perform uncoupled simultaneous multiple degree-of-freedom motions.Finally, Section 6 is devoted to the conclusions of the paper, along with proposals for immediate future works.

System Description
A brief description of the proposed family of generators is now presented (see [27], a.e.).These generators were conceived to extract energy from marine currents and to fulfill the following four main features: (i) being valid for operation at depths greater than 40 m, (ii) minimum installation support structure and civil works, (iii) being floating and easily transportable device, and (iv) being fully automated for emersion and immersion maneuvers.
Figure 1 shows the general view of one of the proposed devices.It is composed of a three-fixed pitch blade propeller and a central pod (gondola), and the shaft of the propeller is coupled to an electrical generator by means of a multiplier gearbox.Three main radial and symmetrically distributed columns start at the pod and end at three torpedoes of an approximately cylindrical shape.These torpedoes contain the inner ballast system used to apply hydrostatic forces to the generator.Another of the principal missions of these torpedoes, which are aligned to the direction of the current, is to minimize rotations around the main axis of the whole generator by using hydrodynamic force compensation when it is extracting energy and the propeller reaction causes heeling torques.Three structural bars starting at the torpedoes and ending in front of the propeller have been added for mooring purposes.
The main features of two of the conceived generators are briefly presented in Table 1.
Various anchoring or mooring systems can be used to place the generator in the sea at the desired position.Figure 2 shows a simple mooring system based on wires and buoys that allows the device to be positioned at the desired depth.By simply providing the generator with positive floatability, the wires become under strain and the device can be moored in both the desired position and depth and is ready to harness energy.
When employing the emersion maneuver procedure, the device has to evolve from its vertical and submerged operating position and orientation (usually denominated as posture or attitude) to its maintenance and transportation position and orientation (floating horizontally) (see Figure 3).This movement is achieved in three basic sequential steps: (1) by providing the device with zero buoyancy rather than positive buoyancy; (2) by releasing the stern wire; and (3) by managing the water in the ballast tanks.

Dynamic Model
Obtaining a dynamic model of the device shown previously under three DoFs in order to perform submerged motions is an important step for the simulation and design control algorithms needed to perform an emersion maneuver in a fully automatic closed loop mode.The proposed model is based on the following assumptions (see Figure 4 for details): (i) Only four lumped masses are considered.
(ii) Added masses are constant, well known, and those which are lumped are also considered.
(iii) Vertical translation (one DoF) dynamics and rotation (two DoFs) dynamics are fully uncoupled.
(iv) Viscous friction is modeled as being constant and fully uncoupled.
(v) The device free-surface interaction (considered semisubmergible) is not considered in this paper.
(vi) The influences of residual flow of water, waves, or wind effects are considered as external perturbations in the model.
A fixed reference frame S 0 is defined as an orthogonal reference frame whose -axis is vertical.Zero depth ( = 0) is computed at the level of the sea (free surface) and the vertical plane (, ) always contains the center of the generator.Only the depth of the center of the generator is related to this reference frame, and components  and  of the generator with regard to S 0 are not considered.The other two degrees of freedom (only rotations around axes   and   are possible) are defined with regard to the other intermediate reference frame, S  , whose origin is at the center gravity of the generator and whose plane (  ,   ) is placed on an absolutely horizontal plane (parallel to the sea surface).If the generator is positioned with null orientation (  = 0,   = 0), the  axis remains perpendicular to the plane of the generator and the   -axis is placed vertically and forms a symmetry axis that is aligned with the mass denoted as  3 .
If we consider now the actuators, which are responsible for creating the buoyancy forces, they are ideally located at the same positions at which each of these lumped masses is considered to be placed.These actuators will produce vertical component forces only that are obtained as differences between buoyancy forces (really produced by actuators) and forces due to gravity.Each of these masses (denoted and numbered with subindex  = 1, 2, 3) is formed of both real mass   and added mass  ADD and is where   is the volume fraction of the generator corresponding to the th-mass.Under static equilibrium conditions (    =   ), no forces are produced and the generator will remain in its previous state without any kind of motion.In the opposite sense, if motion of the generator is desired, each actuator will produce an incremental force that has to be computed as an increase in either mass or volume (with the opposite sign).For the sake of clarity, volume increments are considered to be responsible for creating buoyancy forces.Consider  The force vector applied to the generator with regard to the fixed reference frame S 0 (also with regard to S  ) can be expressed in a matrix form as the sum of all the vertical forces produced by actuators: ) . (4)

Forces and Torques Involved.
Let us consider a local frame S placed at the center of the central cylinder of the generator at which its center of gravity is located.This S frame is defined as follows: the -axis is perpendicular to the plane of the generator, the -axis is located in the opposite direction to the lower mass  3 (see Figure 4), and the V-axis is located at the plane of the generator and forms a right-hand frame.The central mass  CG is then located at the following local position: while the three resting lumped masses are located at the center of gravity of each of the three torpedoes with local coordinates (for the sake of simplicity, notations  ≡ cos() and  ≡ sin() are used): ) , ) , According to the definition of the second frame, S  , the relation of the orientation of both the S and S  frames can easily be obtained by computing the rotation matrix that relates the local generator orientation with regard to the S  reference system.This rotation is the result of composing the next two elementary rotations with regard to the next axes (first horizontal   -axis and then horizontal   -axis basic rotations): The positions of the three lumped masses with regard to S  can therefore be written as ) , ) , Moreover, the application of each of the forces caused by actuators (see forces  1 ,  2 , and  3 in Figure 4) to their respective centers of gravity on each torpedo leads to a torque vector (with regard to S  ), which is defined and obtained as follows: And, by substituting ( 9) in (10), one obtains As expected, from (11), it is proven that it is impossible to cause changes in the horizontal direction (rotations about  axis or -axis) with only vertical forces.

Force Conversion.
In this subsection, a matrix relation between the forces applied to the generator and a new generalized force vector is proposed.This relation makes it possible to decouple the motions of the generator in the sense that different sets of forces can be applied which then produce independent motions for each of the three degrees of freedom of the generator.By rearranging (4) and (11), the following relation is easily obtained: By substituting 60 ∘ and 60 ∘ for their respective values and rearranging (12), the following is attained: Or in a compact form, where F  = ( 1  2  3 )  denotes the buoyancy forces provided by all the actuators,  = (Γ  Γ    )  denotes the proposed vector of generalized forces, and Λ denotes the matrix that relates  to F  , which depends on   and   according to ( 13) and ( 14).In ( 14), vector  is conceived as the uncoupled generalized forces vector applied to the center of the generator, so the first component will produce a torque that is aligned to the  axis and the second component will produce a torque that is aligned to the   -axis while the third component will produce a force that is aligned to the -axis.

Dynamic Submodel for Only Vertical Movement.
The dynamics of an underwater body that is considered to be a rigid body when it performs vertical motions is well known (see [35][36][37], a.e.).As was seen in the previous section, depth  is the only variable of position that can be controlled.Using the same notation given in the previous subsection and under hypothesis shown above, the next dynamic model is proposed for only vertical motions (the effects resulting from freesurface interaction are not computed): where   =  CG +  1 +  2 +  3 denotes the total mass of the generator and is decomposed into the mass of the central cylinder  CG and the three masses of each torpedo.
For the sake of simplicity, the term  ADD denotes the total added mass resulting from the motion in a viscous fluid (see [38]) which is considered to be constant and well known.A quadratic and speed opposite friction term is computed as , where   is also considered to be constant.

Dynamic Submodel for Only Rotation Movements.
In this subsection, the Lagrange formulation allows us to obtain the equations of motion for the generator provided with only rotation movements.The two rotation angles   and   around their respective   -axis and   -axis are taken as generalized rotation coordinates q R = (    )  for rotation purposes solely.For the compotation of the Lagrange function, the kinetic and potential energies are computed beforehand.

Local Inertia Matrix. The inertia matrix in local
coordinates is obtained from the lumped mass distribution employed.Because the general rotation motion is of dimension 3, the inertia matrix results are of dimensions 3 × 3: which satisfies that it is a symmetrical and positive definite matrix.This implies that eig(J LOC ) > 0,   =   ,   =   , and   =   .In our case, the inertia matrix incorporates the effects of the added masses and takes the form (details of how the coefficients in this inertia matrix were obtained can be found in Appendix A): ) ) ) As in the case shown above, each of the masses considered is computed as the sum of its real mass fraction and its added mass, as seen in ( 1).It will be observed that if  1 =  2 , the generator becomes fully statically and dynamically equilibrated with all of its null value inertia products (  =   =   = 0).

Computation of the Rotation Kinetic Energy.
The kinetic energy of the generator when only rotations over first the   -axis and then the   -axis are considered is given by with Ω  , Ω  , and Ω  being the absolute rotational velocities of the generator with respect to S  .By substituting ( 8) and (17), one obtains

Computation of the Rotation Potential Energy.
In order to control the generator, it is necessary to obtain operation conditions that force the generator to be at an equilibrium point.Neutral buoyancy is therefore required (  =   ), and the total potential energy thus has a zero value under this condition 3.4.4.Lagrange Formulation.The kinetic and the potential energy computed previously are used, and the Lagrange function of the generator with only rotation motions with regard to the   -axis and   -axis is Lagrange's equations are then expressed by where  is the sum of external torques, which are the control torques (Γ  Γ  )  and the friction terms that have been modeled as a function F k ( q R ) with  R = (    )  .The following expression is obtained (see the coefficients of M R and C R in Appendix B): in which q R denotes the generalized rotation coordinates, q R = (    )  with regard to S G frame, M R are the inertia terms which are q R dependent, positive definite, and symmetrical, and C R includes both the centrifugal and Coriolis effects.A term  R = (    )  allows the friction terms to be modeled in a similar way to that occurring in (15).

Proposed Dynamic Model.
The complete proposed dynamic model can be obtained by joining (13), (15), and (23) and expressed in a compact form as where the complete vector of generalized coordinates q is and matrices M, C, and  are ) .( 28) ree-DoF tidal energy converter

Rotational dynamics
Translational dynamics

Torque disturbances
Force disturbance Figure 5 depicts the proposed dynamic model, along with the two submodels presented above.The set of input signals are the buoyancy forces produced by actuators F  = ( 1  2  3 )  while the output measurable signals are the twodegree-of-freedom orientation angles together with the generator depth denoted previously as a vector of the generalized coordinates q = (    )  .

Uncoupled Multivariable PD Control
The three-degree-of-freedom model proposed clearly exhibits a nonlinear, time variant, and multivariable uncoupled behavior.A new and simple control scheme based on the following three stages is proposed: (i) an uncoupled term that allows the generalized control torques to be converted into forces applied to each of the torpedoes (by means of actuators), (ii) a nonlinear uncoupling model-based term that ideally cancels out the nonlinearities of the system caused by centrifugal and Coriolis torques, and (iii) a multivariable diagonal proportional-derivative (PD) controller.

Uncoupling Matrix.
The determinant of the matrix Λ(q) defined in ( 13) and ( 14) is This determinant makes matrix Λ(q) nonsingular, with the exception of   = 0 which corresponds to the normal operation orientation.In other words, if the generator is placed in a fully vertically position, it is not possible to turn around the   -axis using vertical forces alone.In order to avoid this singularity, the position of the third actuator can be displaced a small distance   along the local -axis, thus resulting in its new local position at (see Figure 6) For relatively small values of   the angle of the new plane that conforms to the three points of application of forces with regard to the plane of the generator is (see Figure 6) and the range of the angle   is −(/2 + 2  /3) ≤   ≤ −2  /3 rather than −/2 ≤   ≤ 0, thus avoiding the singularity of matrix Λ(q) in ( 29).For the angle   is considered range −/2 ≤   ≤ /2.Then, after avoiding the singular values of the angle   , the inverse of matrix Λ(q) can be obtained symbolically, yielding ) ) ) ) ) . (32)

Proposed Control System.
The existence of the inverse of matrix Λ(q) into ( 29), (32) makes it possible to the generator in an uncoupled mode and even in open loop, by simply obtaining the forces to be applied by actuators F  from a set of desired and uncoupled generalized forces  for nonsingular values of   .The existence of external disturbances and model uncertainties between the real matrix Λ(q) and the computed matrix Λ(q) requires the addition of closed loop controllers.Equations ( 24) and ( 32) provide the following control system, which is proposed for a given desired reference q  = (      )  and which can be time dependent: where Λ(q) and Ĉ(q, q ) denote their respective computed matrices as a function of the measured output signals q.
The complete dynamics of the system controller, under the assumption of perfect cancellation, that is, Ĉ(q, q ) = C(q, q ) and Λ(q) = Λ(q), results in Since the inertia matrix M(q) is defined according to (26) and the coefficients shown in Appendix A, it is a bounded symmetric and positive definite matrix which satisfies that  min I 3×3 ≤ M(q) ≤  MAX I 3×3 for all q.Upon considering a mean inertia matrix as the closed loop controlled system with mean matrix M exhibits fully decoupled multivariable linear dynamics defined by and three independent single input/single output (SISO) PD controllers with derivative in the feedback loops can be tuned by choosing diagonal matrices of gains K D = diag(      ) and K P = diag(      ) that allow closed loop poles to be placed where desired.Equation ( 36) can then be analyzed as three SISO independent closed loop systems for any of the three degrees of freedom of the generator.Figure 7 depicts the proposed control scheme.The uncoupling matrix and the linearization term that compensates the centrifugal and Coriolis torques are clearly shown, together with an uncoupled three-dimensional diagonal PD controller with derivative component from only the output signal as proposed in (33).

Numerical Simulations Results
The proposed dynamic model (24) and the controller (33), presented above, were used to carry out numerical simulations in order to demonstrate the goodness of both the dynamic model and the proposed control in terms of generator controllability and ability to perform emersion and immersion maneuvers with only passive buoyancy . .The main reason for using a laboratory prototype is the high economic cost of a real scale device.The research group has successfully tested a 1 : 10 scale prototype [39] in a real protected sea environment.From these trials, the high acceleration values observed in open loop in some of the time intervals during the execution of the emersion maneuvers performed (see Figure 22 in [39]) could damage the rotor blades making necessary the study of the proposed closed loop control.
The first set of numerical simulations was performed on step reference signals with nominal parameters of the added mass and the friction term in order to validate both the desired transient and steady state responses and there were no coupled dynamics among the independent movements (  ,   , and ) that the generator can perform.Figures 8 to 10 The PD controller gains for each of the loops were tuned using the nominal value of matrix M in order to obtain closed loop double poles placed at  1,2 = −0.3rad/s and a unit gain that will guarantee critically damped second-order dynamics.For a given fully uncoupled open loop dynamic (second-order and type one system) of the form   () =   /( ⋅ ( +   )) and a desired closed loop dynamic of the form   () =  2 1,2 /( 2 + 2 1,2  +  2  1,2 ), for  =   ,   , , the gains of the PD controllers are easily obtained as   =  2  1,2 /  and   = (2 1,2 −   )/  .Actuator saturation is not considered here.Gaussian noises from a normal distribution (0, 3 −3 ) were added to each of the output signals, and first-order filters were used to obtain the time derivative of the output signals as L(/) =  ≈ /(0.001+ 1), where L denotes the Laplace transform.
Figures 8-10 depict that the system performs extremely well with uncoupled dynamics, unit gain, and the desired setting time response with no overshoot (the effect of the   transition at  = 20 s and the real   response at  ≥ 20 s is not perceived in the time evolution of the response of   ).
The second set of simulations were carried out assuming nonnominal terms (increments of ±20% of the nominal added mass term and null Coriolis and centripetal compensation term).Additionally, the actuators were saturated to ±5 N, thus simulating a more realistic scenario.Due to the feedback robustness provided by the PD controller (see [40], a.e.), the time responses were very similar to the nominal responses shown in Figures 8-10.Figures 11-13 depict the evolution of the errors between the nonnominal and the nominal responses.In these figures is clearly observed the coupling effect at instants  = 10 s and  = 20 s when transitions occur due to step signals.
The third set of simulations was carried out using synchronous smooth time trigonometric sixth-order Strajectories (see [41,42]) with starting and ending times of  0 = 20 to   = 120 s, respectively, for the orientation Figures 17-19 illustrate the excellent response of the system during a complete emersion sequence.It will be observed that the initial posture q 0 is equivalent to that shown in Figure 2, while the final posture q  is equivalent to that shown in Figure 3 (right), corresponding with the normal submerged operation posture for energy harnessing and the normal floating posture for maintenance tasks.Very similar results were obtained by integrating the dynamics of the laboratory prototype into the simulation environment Orcina OrcaFlex [43].OrcaFlex is one of the world's leading software packages for the design and analysis of a wide range of marine systems (riser systems, mooring systems, installation planning, towed systems, marine renewables, etc.) and is considered validated and certified software for anchor marine systems.The control signals given in (33) were computed in the runtime of Matlab which is connected to OrcaFlex through an external DLL function which is responsible for obtaining the magnitudes from the simulated responses and providing the control signals from and to OrcaFlex with a time sampling of 10 ms. Figure 20 displays a visual sequence of the response of the 3-DoF prototype when following the references given in (39).

Conclusions
This paper proposes a new dynamics model that can be used to control a family of submarine electrical generators that was conceived to harness energy from marine currents.The submarine generator is provided with three degrees of freedom by which it is possible to perform closed loop emersion and immersion maneuvers.The dynamic model proposed is based solely on the definition of four lumped masses placed on a plane and a matrix that permits the forces from the three buoyancy actuators to be converted into two torques and a vertical force that are responsible for the rotations and the vertical displacement of the whole system.The model developed exhibits a nonlinear, strongly coupled, and timevarying behavior between buoyancy forces produced by the three actuators and the posture magnitudes measured using a depth sensor and a three-dimensional inclinometer on which the angle rotation around the -axis is ignored.As a first conclusion, we should note that the dynamic model developed in this research is sufficiently precise to describe the underwater 3-DoF tidal energy converter motions and sufficiently simple to be used in the control law design.
The proposed control scheme with which to control the generator using only depth and orientation measurements is, meanwhile, based on three stages: (a) the definition of an uncoupling matrix is used to decouple the generator motions for nonsingular values of orientation; this matrix can be easily computed in real-time and allows the generator to be handled in an uncoupled manner, even in open loop; (b) it is a nonlinear uncoupling-based model for the compensation of the centrifugal and Coriolis torques; and (c) it is a proportionalderivative (PD) control action.The closed loop dynamics was chosen by simply designing proportional and derivative matrix gains.The proposed control allows the generator to perform fully automatic closed loop emersion/immersion maneuvers from a submerged position with vertical orientation to a floating posture with a horizontal orientation by using time-varying smooth reference trajectories.A last obtained conclusion is that the control law implemented in this work demonstrates that it is a simple control strategy, which is computationally efficient and easily implementable in a microprocessor-/microcontroller-based system.
Finally, based on more than satisfactory numerical simulations achieved in this research, various experimental branches are now the focus of our attention and are detailed as follows.Our intention is (a) to finalize the construction of a real small-sized laboratory prototype; (b) to use realtime experiments to validate the proposed research; and (c) to study different control strategies and new trajectory generations in order to improve the quality of the closed loop emersion/immersion maneuvers.These will be the main topics of our future research.

A. Coefficients of the Inertia Matrix J LOC
The local positions of each of the three masses located at their respective center of gravity on each torpedo and which are taken from ( 6) are used to obtain the equivalent inertia matrix J LOC with regard to the local reference system.It will be noted that the central mass  CG does not produce any inertia value in accordance with its local position given by ( 5) B. Obtaining Matrices M(q R ) and C(q, q R ) Equation ( 23) was obtained from (19)-( 22).This appendix shows how the coefficients of matrices M(q R ) and C(q R , q R ) were obtained.Firstly, partial derivatives with regard to

2 Mathematical
Problems in Engineering(ii) There is the use of elevation and placement by means of floating cranes.See [19-23].

Figure 1 :
Figure 1: Main appearance of a generator for unidirectional currents.

Figure 3 :
Figure 3: Final stages of the emersion maneuver.
illustrate the time responses of each of the independent closed loop systems of each of the degrees of freedom of the generator when the following step reference signals were used (units in rad for angles, m for depth, resp.): from 0 to  3 at  =10 s,   from − 0.05 to −  2 at  = 20 s,   from − 20 to − 1 at  = 20 s.

Figures 13 -Figure 17 :Figure 18 :Figure 19 :
Figures[13][14][15][16] illustrate that the generator exhibits a good correspondence with both the desired responses and the simulated responses obtained.Finally, the fourth set of simulations corresponds to a complete sequence of sixth-order S-trajectories to achieve an emersion maneuver.Figures 17 to 19 depict both the references and the time responses obtained for each of the independent movements of the generator.The desired sequence is defined below (units in rad for angles, m for depth):

Table 2 :
Main parameters for simulation (per each of the cylindrical torpedoes).