Investigation of an Underwater Vectored Thruster Based on 3RPS Parallel Manipulator

Autonomous underwater vehicles (AUVs) are important and useful tool platforms in exploring and utilizing ocean resource. However, the effect of control surfaces would decrease even invalid complete in this condition, and it is very hard for conventional AUVs to perform detailed missions at a low forward speed. ,erefore, solving this problem of AUVs becomes particularly important to increase the application scope of AUVs. In this paper, we present a design scheme for the vectored thruster AUV based on 3RPS parallel manipulator, which is a kind of parallel manipulator and has advantages of compact structure and reliable performance. To study the performance and characteristics of the proposed thrust-vectoring mechanism, a series of works about corresponding kinematic and dynamic analysis have been performed through the theoretical analysis and numerical simulation. In the part of kinematics, the inverse, forward kinematics, and workspace analysis of the thrust-vectoring mechanism is presented, and the numerical simulations are accomplished to prove the feasibility and effectiveness of this design in AUVs. In order to further verify feasibility of the thrust-vectoring mechanism, based on the considerations of various affecting factors, a dynamic model of the designed thrust-vectoring mechanism is established according to theoretical analysis, and the driving forces of the linear actuator are presented through a series of numerical simulations. In addition, a control scheme based on PID algorithm is proposed for the designed vectored thruster with considering various affecting factors and the application environment. Meanwhile, the control scheme is also established and verified inMATLAB SimscapeMutibody. A series of numerical simulations of the thrust-vectoring mechanism prove the feasibility of the vectored thruster. According to equipping the designed vectored thruster, the AUVs can overcome the limit of weakening the control ability at zero or low forward speeds, and this improvement also expands the application of it, which has been scaled greatly.


Introduction
Over the last few decades, due to exhausting of resources and energy, human beings are bearing with a series of survival predicaments and development challenges [1,2]. Because of the lack of land-based resource and the continuing need for all kinds of resources, an increasing number of countries and scientists have paid more and more attention on the exploitation and utilization of resources [3][4][5]. In the present, most of the water available on Earth exists in the oceans, yet only a small part of this vast resource has been explored [2]. e ocean has vast areas and is rich in all kinds of natural resources, such as marine life, oil, natural gas, and minerals. Additionally, the ocean not only contains a lot of marine recourse but also brings a lot of traffic convenience around the world. With the progress of society and economy, the development of mineral resource has become an inevitable trend. Exploring and exploiting the oceans has become the principal development strategy of every country in the world. However, the nature environment of the ocean is too harsh to explore, and the advanced technology has been rapidly developed in recent years, such as autonomous underwater vehicles (AUVs), remotely operated vehicles (ROVs) [6][7][8], and unmanned marine vehicles (USVs) [9]. In addition, other different techniques have also been used for controlling all kinds of robots, such as proportional integral derivate [10], fuzzy control [11,12], and sliding model control [13].
AUVs have become a main tool for surveying below the sea due to the great improvement in their performance and advancement in underwater research. rough equipping a large quantity of advanced instruments and equipment, AUVs are capable of accomplishing applications including scientific, commercial, and military tasks such as exploration of oceans [14,15], oceanography mapping surveys [16][17][18], the collecting ocean environment information [19][20][21], and searching and rescuing for shipwrecks [22,23] and debris from the missing airplanes [24,25]. With the expanding area of applications, the design of AUVs needs to meet the higher demand continuously. Although it has made great steps in AUVs performance, the new emerging technologies and demands for exploiting oceans have attracted critical mass of scientists and engineers to undertake the research of AUVs.
Conventional AUVs are designed equipped with a main propeller and control surfaces at the tail cone for propulsion and control [26][27][28]. ese conventional AUVs can fulfill the work well under normal conditions. When the conventional AUVs need to complete exploration tasks with a lower speed in a complex and unknown underwater environment, the control capability of AUV depends heavily on the control faces made up of fin and rudder. e velocity of AUV is relatively low or zero because of the demand of practical problems, such as scanning topographic map, taking photographs, and monitoring marine observation data. However, they are unable to perform detailed inspection missions at zero or low forward speeds due to the control faces which become ineffective in this condition [28,29]. e cause of this problem is that the generation of control forces from control surfaces depends on forward speeds of AUVs [28,30]. erefore, this disadvantage has further development and application of conventional AUVs greatly.
ere are some approaches to solve this problem, such as installing additional thrusters to provide additional control forces for controlling AUVs [28,[31][32][33][34], but this method would result in the problems of complicated structure and increasing energy consumption. Its complex structure, adverse working environment, and so on causes the decrease of reliability of the whole AUV control system. Installation and maintenance of additional thrusters would significantly increase energy expenditure or energy carrier for sailing.
e AUVs equipped with vectored thruster do not require the use of fin and rudder for controlling at all. Since this kind of AUVs driven by a vectored thruster, the control forces are generalized force components produced by vectored thruster, and these forces only depend on the rotational speed of the propeller. erefore, the AUVs equipped with vectored thruster are independent from any control forces generated by control faces, and the controllability of vehicle is markedly improved and obtains a better good application effect. So, the vectored thruster AUVs are capable of accomplishing detailed missions at a low forward speed.
In the research areas of this field, some companies and research institutions have made progress in theory study and application of vectored thrusters [29,36,38]. Among above research AUVs, Bluefin and MBARI have achieved great successes and provided considerable experience in the use and study of vectored thrusters. More importantly, the engineering practice of Bluefin and MBARI shows that this method can raise control efficiency greatly and also reduce the possibility of losing control at low speed. However, the existing design of the vectored thruster is almost designed based on serial mechanism; this kind of mechanism has the disadvantages of complex structure, low bearing capacity, and high moment inertia. So, based heavily on practicalities of serial mechanism, the existing design of the vectored thruster is too bulky and complex to use for AUVs. Considering the restrictions of application environment and structure size, it is crucial to choose a suitable mechanism for designing new vectored thruster AUVs.
On the contrary, compared with other commonly used mechanical structures, parallel manipulators have numerous advantages, such as small size, compact and reasonable structure, reliable performance, fast response, high positioning precision, high stability, high sensitivity, high stiffness, and better dynamic performance [39][40][41][42]. ose merits of parallel manipulators make the device have high popularization value and use value, such as medical and industrial robots, flight simulator, and mechanical device. Inspired by various applications of parallel manipulators, the idea of vectored thruster AUVs based on parallel manipulator is generated. In the field of thrust-vectoring mechanism research, many scientists and engineers have made great contributions to the development of vectored thruster AUVs based on parallel manipulators. e full deflection vectored thruster is based on the spatial linkage and universal joint proposed by Cavallo and Michelini [43]; the authors designed a 3-SPS-S parallel manipulator with passive constraining spherical joints to drive the underwater vehicle [30]. e above thrusters currently have some problems, such as the structure is relatively complex, and the motion real-time resolving method and the dynamics model for the vectored thruster are difficult.
With comparing structure characteristics of different kinds of parallel manipulators and considering actors of application environment, 3RPS is chosen from various parallel manipulators as the thrust-vectoring mechanism mainly various advantages, including its compact structure, high position tracking precision, and fast response speed.
is parallel mechanism is a strong coupled nonlinear structure, so its motion control is too complex to use more widely [44][45][46][47]. Despite its advantages, the 3RPS parallel manipulator also needs to overcome some problems that would restrict the development and application of the thrust-vectoring mechanism. rough reading and analyzing the domestic and foreign related literature, various methods of kinematics and dynamics for 3RPS parallel manipulators have been presented [44,45,[48][49][50].
On the basis of the above considerations, the design concept of vectored thruster which is made up of 3RPS parallel manipulators is introduced. e vectored thruster based on 3RPS parallel manipulators has terse structure, convenient operation, convenient installation, steady working system, and wide adjustable range. Using this method, the AUVs are able to provide the vectored thrust effectively and efficiently. More than anything, the AUVs equipped with vectored thruster are able to complete a variety of the complex tasks at a comparably low forward speed.
In this paper, the structural design of the vectored thruster based on 3RPS is introduced briefly. In order to satisfy the design requirements and study the motion characteristics of vectored thruster, the kinematics and dynamics model of the thrust-vectoring mechanism are established, and the related simulation is presented to verify feasibility of the scheme. Finally, a control scheme for the vectored thruster is designed and simulated in Matlab. e theoretical analysis and numerical simulations prove that the proposed vectored thruster based on a 3RPS parallel manipulator can effectively realize the function of providing the required vectored thrust for thrust-vectoring propulsion.

The General Design of Vectored Thruster
e configuration of the whole AUV equipped with the designed vectored thruster based on 3RPS parallel manipulator is presented, as shown in Figure 1. Due to the existence of the vectored thruster, the AUVs do not need any more rudders to provide control forces. e force generated by the vectored thruster is used as control force for AUV's yaw and pitch motion. Consequently, the tilt angle is important one of the criteria to assessing the performance of the proposed vectored thruster. However, the space in the stern of AUVs is limited, and the tilt angle range of the designed vectored thruster is also limited. Referring to a literature review [29,51], the duct propeller's tilt angle is limited to plus or minus 15°in our design. In addition, the vectored thruster contains a duct that can be used for protecting the propeller form damaging and enhancing flow capability.
In terms of structure design, considering the specific requirements of application environment and the stability of system, we adopted the modular design for vectored thruster AUV. e designed vectored thruster is mounted on the stern of an AUV as an integral and independent, which is adopted for convenient installation and maintenance. e designed vectored thrust duct propeller system mainly contains the duct propeller and the thrust-vectoring mechanism. A whole structure model of the vectored thruster AUV based on 3RPS parallel manipulator is built up, as shown in Figure 2.
At present, the duct propeller is the most widely used form of propulsion device for underwater robots. A duct propeller is mainly composed of an annular wing and a propeller.
ere are many underwater vehicles equipped with duct propellers, for the extraordinary performance of improving the propulsive efficiency and avoiding cavitation conditions [52]. is kind of propeller is able to provide the thruster from zero to cruising speed more effectively. Just because of an effectively increased thrust in the condition of a low forward speed, the duct propeller is widely used in various marine vessels, such as AUVs and ROVs.
Furthermore, underwater environments are very complex and harsh; propellers are very easily destroyed by underwater animals and plants, waves, even currents, and other uncertainties. Hence, the existence of a duct can protect the propeller against damage from the underwater environment during all kinds of missions. Moreover, since the duct also can generate the thrust during the voyage, the duct is an important source of control force for AUV's yaw and pitch motion.
In our design, the duct propeller is driven by a main electric motor installed in the rotating body, which is aligned with the holes of the duct's inner shaft with fastening screws. In order to simplify the unnecessary transmission structure and reduce the redundant weight, the main motor has been bedded on the rotating body with a duct propeller. It is very clean and efficient to take direct connects with the propeller and the main motor.
is installation mode of the main motor and rotating body can improve space utilization significantly and reduce the weight of vectored thruster effectively. And this modularity makes the duct propeller system easy to maintain and debug on the whole vectored thruster control system.
As the implementing actuator of the vectored thruster, the thrust-vectoring mechanism is fundamental to the overall system for its basic functions.
ere are many methods on how to realize thrust vectoring, and each method has its own advantages and disadvantages. Considering the limited space of AUVs' tailcone and the harsh operation condition, it is central to choose an appropriate mechanism structure that can complete the design function of achieving vector control effectively for AUVs. Comparing to the serial mechanism, parallel manipulators have many inherent superiorities, such as small size, compact and reasonable structure, reliable performance, fast response, high positioning precision, high stability, high sensitivity, high stiffness, and better dynamic performance.
Integrating practical application environment of AUVs and based on the application background of various parallel manipulators, 3RPS parallel manipulator is chosen as the thrust-vectoring mechanism after analyzing various mechanical structures. In accordance with this notion, a novel thrust-vectoring mechanism based on the 3RPS parallel manipulator for AUVs is designed, as shown in Figure 3. e thrust-vectoring mechanism is designed based on 3-RPS manipulator, which has a top rotting platform, a fixed base, and three identical sets of driving limbs and joints. Each driving limb has an actuating prismatic joint (P) attached to the fixed base by a rotational joint (R) and connected to the platform by a spherical joint (S) [53][54][55]. Each limb is driven by a linear actuator. us, the length of limb could be changed within the operation range, and the top platform would rotate when the length of limbs changes with a certain law of motion. e vectored thruster is installed on the tail horizontally, which connected with the shell of the AUV via fastening screws. e duct propeller can rotate around the center of the top platform horizontally and vertically, and the thrust generated by the vectored thruster will drive the vehicle forward or changes the direction of movement. e 3RPS parallel manipulator has two rotational and one translational degree of freedom (DOF). It is superfluous to have the translational DOF for the thrust-vectoring mechanism, the redundant DOF needs to be constrained by motion control, and the other two rotational DOFs are used to realize functions of the thrust-vectoring mechanism. In addition, the translational DOF of 3RPS parallel manipulator will make the top rotating body bump into the shell of the AUV. So, the importance of the study on redundant DOF of the 3RPS parallel manipulator is obvious for the actual application of the thrustvectoring mechanism.
Since the vectored thruster could generate required control forces for controlling AUVs motion, there is no need to have extra rudders as conventional AUVs. e component of the thrust as control forces is dependent on the deflection angle and the thrust of the vectored thruster. erefore, the research on deflection angle of the vectored thruster is essential for controlling the motions of AUVs. However, it is very difficult to measure the tilt angle of the vectored thruster directly because the limited space and underwater environment is not suitable for installing sensors to measure. Another common and efficient approach to get the rotation angles is using the kinematic analysis method, which can obtain the tilt angle by measuring the lengths of the three limbs. Based on this kinematic method, tilt angle information about the vectored thruster can be obtained via relative calculation with the lengths of three limbs, which can be measured directly by length sensors installed in actuators.
In order to realize precision and stable positioning control of the proposed vectored thruster, the design of the automation control system is fundamental to achieve objectives. Hence, establishment of kinematic and dynamic models for the thrust-vectoring mechanism based on the 3RPS parallel manipulator is significant to achieving perfect control of the vectored thruster based on the above analysis.

Kinematic Analysis of the Thrust-
Vectoring Mechanism e thrust-vectoring mechanism is designed based on the 3-RPS parallel manipulator, which is composed of a base plate, a rotating platform, and three uniformly distributed driving   limbs, as shown in Figure 3. According to the needs of analyzing the motion of the top rotating platform, two Cartesian coordinate systems with associated symbols have been established in the 3RPS parallel manipulator and shown in Figure 4. e reference frame O-xyz, which is the global coordinate system, is fixed to the center of the immovable base and the z-axis normal to the fixed base. Similarly, the reference frame P-ijk denotes the local coordinate system located on center point P of the rotating platform, whose j-axis is perpendicular to the bottom surface of the platform. e moving sides of driving limbs (linear actuators) are connected to the upper rotating platform through three sphere joints that is fixed directly to the center of the top platform, while the other sides of the limbs are connected to the base with three revolute joints that are symmetrical about the center of base. A 1 , A 2 , and A 3 are the connected points between the fixed base and the driving legs (linear actuators), B 1 , B 2 , and B 3 denote the points of the revolute joints. It should be mentioned that A 1 B 1 , A 2 B 2 , and A 3 B 3 are perpendicular to the fixed base because A 1 B 1 , A 2 B 2 , and A 3 B 3 represent revolute joints with a certain height. A reference frame O′-xyz is established with respect to the plane formed by points B 1 , B 2 , and B 3 , and this plane parallels with the fixed base A 1 , A 2 , and A 3 . e connected points between the moving parts and the rotating platform are expressed as C 1 , C 2 , and C 3 . e radius of the fixed base and the top platform are defined as r 1 and r 2 , O and P denote the center points of the base and top platform, respectively. L 1 , L 2 , and L 3 denote the lengths of three linear actuators between the top platform and the fixed base.
As we can see in Figure 4, in the global reference frame O-xyz, the center point of the equilateral triangle made up of three points A 1 , A 2 , and A 3 is expressed as O, and the radius of the fixed base is defined as OA 1 � OA 2 � OA 3 � r 1 . Hence, the location of A i in global reference frame O-xyz can be described as follows: Similarly, B 1 , B 2 , and B 3 denote the axes of revolution of the revolute joints with a certain height h r , and this plane parallel is with the fixed base. Hence, the locations of point B i can be denoted as follows: A local coordinate system P-ijk is established on the rotating platform bottom surface of the thrust-vectoring mechanism, in which the origin point P is the circumcenter of triangle C 1 , C 2 , and C 3 . So, the locations of connection point between the linear actuators and the top platform can be described as follows: From Figure 4, p as a position vector denotes the translation vector from the center point O to point P of top rotating platform in global reference frame O-xyz. To describe the relative motion between top rotating platform and fixed base, a rotation matrix about frame P-ijk with respect to the fixed base reference frame O-xyz needs to be established. e position vector p and the rotation matrix R are defined as follows: where s (·) � sin (·), c (·) � cos (·), and α, β, and c denote the rotation angles about the k-axis, j-axis, and i-axis, respectively. e thrust-vectoring mechanism only needs two rotational DOFs to realize its functionality; the 3RPS parallel manipulator has one more translational DOF that is redundant. In order to present the condition of the top rotating platform of the thrust-vectoring mechanism, the rotation angles are also important parameters that need to be defined. According to the need of the thrust-vectoring mechanism, the displacement between centers of the base and the top platform are set as h. A generalized vector q is established to describe the position and orientation of the top rotating platform in the global reference frame as follows: where P x , P y , and c are associated with the rotation angles α and β. Based on P x � r 2 (sαsβcc − cαsc), P x � r 2 (sαsβcc − cαsc), P x � r 2 (sαsβcc − cαsc), P x � r 2 (sαsβcc− cαsc), P y � (r 2 /2)[(cβcc − sαsβsc) − cαcc], L P � ������ � P 2 x + P 2 y , and c � tan −1 [sβsα/(sβ + cα)], so P x and P y also can be written as follows: where δ Px , δ Py , and δ Ly denote the influence factor of P x , P y , and L P . According to the application needs to be designed in Mathematical Problems in Engineering this design, the numerical simulations about δ Px , δ Py , and δ LP and the angle of rotation of c can be calculated, and the simulation results are obtained and plotted in Figure 5. From Sections 3.1 and 3.3, the inverse, forward kinematic analyses, and workspace analysis of the thrust-vectoring mechanism are performed, and numerical simulations are conducted to validate the accuracy and reliability accordingly.

Inverse Position Analysis of rust-Vectoring Mechanism.
In this design, inverse position analysis of the thrust-vectoring mechanism is carried out to establish the mapping relations between the position and orientation of top moving platform and the lengths of three driving linear actuators. Referring to Figure 4, the length of linear actuator L i with respect to the fixed base reference frame O-xyz can be written as where a i (c i ) denotes the vector from point O (P) to point A i (C i ) in frame O-xyz (P-ijk). e length change of the ith limb can be calculated as where L ave is the initial length of the linear actuators at the tilt angle α � β � 0°. According to (7) and (8), the length changes of the three linear actuators can be obtained with related parameters presented in Table 1, and the dimension parameters of vectored thruster are calculated through CAD software. e results of the length change of ith linear actuator are plotted in Figure 6.
To further study the relationship between the tilt angles of the top rotating platform and the lengths of three linear actuators, the top moving platform moves according to α s � π/9 · sin(t) rad and β s � π/9 · cos(t) rad. When the top rotating platform moves according to α � α s and β � β s , based on the kinematic analysis of the thrustvectoring mechanism above, the length of linear actuators is plotted in Figure 7.

Forward Position Analysis of
rust-Vectoring Mechanism. Similarly, the forward position analysis of the thrust-vectoring mechanism is established to study the mapping relations between the lengths, three linear actuators, and the position and orientation of top moving platform. e position and orientation of the top moving platform is obtained according to the given length of the ith linear actuator.
Referring to Figure 4, the position vector of point C i on the top rotating platform in global frame O-xyz can be expressed as where B i C i denotes the vector of the ith linear actuator, which can be expressed as where θ 1i is the angle between the actuator and the fixed base. Since points C 1 , C 2 , and C 3 form an equilateral triangle in the top rotating platform, based on the theory in geometry, the relationship of C 1 C 2 C 3 can be determined by The moving platform The fixed base  Mathematical Problems in Engineering Hence, the position of points C i in reference frame O-xyz can be obtained through (9)- (11).
Because P is the center of the circumcircle generated by three points C 1 C 2 C 3 in the top moving platform, the position of center point P in global reference frame O-xyz can be expressed as where p of the center point P of the top rotating platform, which can be used to describe the position of point P with three components (P x , P y , P z ). e position vector OC i of points C i is regarded as a known parameter when the length of linear actuator is given; thus, the three equations about the point of P can be established and calculated.
To further investigate the relationship between the lengths of linear actuators with the position and orientation of the top platform, different lengths of the linear actuator are used for the forward kinematic analysis. In this simulation, the linear actuator L 1 , L 2 ∈ [310 340] mm and L 3 � 310, 325, 340 mm. According to equations (9)-(12), the orientation and position vector of the top moving platform can be calculated directly, as shown in Figure 8.

Workspace Analysis of rust-Vectoring Mechanism.
Due to the available space of AUV is limited, it is important to analyze the workspace of the thrust-vectoring   Mathematical Problems in Engineering mechanism for optimizing structure design and improving performance. According to the kinematic analysis mentioned above, all the positions and orientations of the top rotating platform can be obtained by changing the lengths of actuators. Considering the motion characteristics of the thrust-vectoring mechanism and constraint on available space, the workspace analysis is mainly referring to study the tilt angle and angle change of the revolute joint and spherical joint of the thrust-vectoring mechanism in this paper. In this section, θ r and θ t denote, respectively, the rotational angles of the revolute joint and spherical joint. e schematic diagram of revolute joint and spherical joint is presented, as shown in Figure 9. e tilt angle θ ri and angle change ∆θ ri of the revolute joint at point B i can be defined by θ ri � arccos s ri · e 1 s ri e 1 , where e 1 is the direction vector of the y-axis in global reference frame O-xyz and θ r−avg is the initial angle generated by the linear actuator and the fixed base at the tilt angle α � β � 0 rad. e tilt angle θ ti and angle change Δθ ti of the spherical joint at point C i can be expressed as θ ti � arccos s ti · t 1 s ti t 1 , where t 1 is the direction vector of the k-axis in local reference frame P-ijk and θ t-avg is the initial angle between the linear actuators with the rotating platform at the tilt angle α � β � 0 rad.
Similarly, to study the relations between the tilt angles of revolute and spherical joints and the lengths of linear actuators, the tilt angles of revolute joint θ ri and spherical joint θ ti by forward kinematic analysis are performed with linear actuator length L 1 , L 2 ∈ [310 340] mm and L 3 � 310, 325, 340 mm. e simulation results are plotted in Figure 10.
e tilt angle of the three linear actuators and the top rotating platform has close relation, the platform moves according to α s and β s mentioned above. When the top rotating platform moves according to α � α s and β � β s (α s � π/9 · sin(π · (t/36) rad and β s � π/9 · cos(π · t/36) rad), the tilt angle and angle change of linear actuators are plotted in Figure 11.

Dynamic Analysis of Thrust-Vectoring Mechanism
For improving the dynamic performance and control accuracy of the designed vectored thruster, it is greatly important to analyze the dynamics model. Since the thrust-vectoring mechanism is designed based on the 3RPS parallel manipulator, which includes three closedloops kinematic chains, it is very complicated to perform the dynamic analysis of the thrust-vectoring mechanism. According to the theoretical analysis and practical needs, the dynamic model of the thrust-vectoring mechanism based on the 3RPS parallel manipulator is established. e schematic diagram of the dynamic analysis model of the 3RPS parallel manipulator is represented graphically, as shown in Figure 12.
Referring to Figure 12, dynamic formulation of the force and moment balances on the linear actuator of the 3RPS parallel manipulator can be expressed as follows:  Mathematical Problems in Engineering thrust-vectoring mechanism. According to the forces and moments distributions analysis of the proposed thrustvectoring mechanism, the stress conditions of the top platform can be represented as Figure 13. As shown in Figure 13, forces and moments are applied on the connection points C i generated from many respects, such as the linear actuator, the top rotating platform, and duct propeller. Based on the definition of forces and moments at point C i above, the dynamical equations of the top rotating platform can be expressed as where v p and _ v P are the velocity and acceleration and ω P and _ ω P denote the angular velocity and angular acceleration at the center of the top platform, respectively. F Ci denotes actuating force from the linear actuator along the direction of actuator. Because the top platform is an axisymmetric structure, F g−P and F b−P are the gravity and buoyancy from the top platform, and l g−P and l b−P denote the distance from the center point P generated by connection points C 1 C 2 C 3 to the center of mass and buoyancy of the top rotating platform, respectively. F e and M e are the external force and moment mainly generated from the propeller and the duct in our paper. Referring to Figure 13, the external force F e and external moment M e can be defined as where F prop and F duct denote the force generated by the propeller and the duct of the thrust-vectoring mechanism and M e denotes the moment generated by the duct, respectively. e thrust vector F prop is produced by the propeller and can be expressed as where T p denotes the thrust produced by propeller and based on standard propeller theory [28,31], T p � K T ρn 2 p D 4 . K T , ρ, n p , and D denote the thrust coefficient, the water density, the rotation speed of the propeller, and the propeller diameter, respectively. In this AUV, the duct propeller has been widely adopted to protect from damage and improve the propulsive efficiency by being enclosed by a duct. To further investigate the dynamic model of the vectored thruster, it is clearly necessary to considerate the effect on lift and drag generated by the duct. e force generated by the duct applied to the platform can be expressed as where L and D denote the lift and drag of the duct, which can be calculated by CFD software. Because the duct rotates around the center of the duct in use, a transformation matrix R d is established to convert the duct frame into the body frame, and the matrix R d can be described as where α and β represent the tilt angles of the ducted propeller.
Referring to equations (17)- (20), the force F duct and moment M duct generated by the duct that are applied on the platform of the thrust-vectoring mechanism can be calculated as Finally, based on the above analysis and according to Figures 12 and (15)- (21), the force balance along the leg direction can be expressed as where F i denotes the force produced by the linear actuator to complete the key components of drive function and s i is the unit vector of the ith actuator. In order to do a better research on the effect of motion on the vectored thruster, the numerical dynamic simulation on the thrust-vectoring mechanism has been developing. Some parameters used in the simulation, such as the dimension parameters of vectored thruster, are calculated through CAD software, and other parameters can be obtained by in [56]. All parameters of the thrust-vectoring mechanism are given in Table 2.
Based on the abovementioned theory analysis and parameters, the analysis formulations in Section 4 have been implemented in MATLAB.
When only considering the thrust produced by the duct propeller and the top platform moving according to α � α s and β � 0, the length changes and the driving forces of linear actuator can be calculated, and the results are plotted in Figures 14 and 15.
In the dynamic analysis above, the gravity, buoyancy, drag, and torque of the vectored thruster is ignored. When the weights of the rotating platform and the three actuators are taken into account only, the driving forces of the linear actuator can be calculated and plotted in Figure 16.
In addition, the buoyancy of actuator depends on the length changes of actuators, and the buoyancy of linear actuator can be calculated by the length change. When the top platform moves according to the designed trajectory, the buoyancy is created with the movement of the vectored thruster and the result of buoyancy is shown in Figure 17.
As we can see in Figure 17, the buoyant forces of actuators change with the length changes of actuators, but the magnitude of the change is relatively small. Hence, the buoyant forces of actuators can be approximately equal to 4.7 N in this next calculation. Generally, the emphases of dynamic analysis of the parallel manipulator for research are mainly focused on the gravity, the external force, and torques. However, the buoyancy is also an important factor that affects the driving force of the actuators because the designed vectored thruster is used in underwater vehicles. In addition, the buoyancy of the other parts can be directly calculated by CAD software.
With considering the buoyant forces of actuators as shown in Figure 17 and buoyancy of other parts, the driving force of the ith actuator can be recalculated and plotted in Figure 18.  Comparing the driving forces in Figures 15, 17, and 18, it is obvious that the gravity and buoyancy of the whole vectored thruster have an effect on the value of the driving force greatly.
In order to investigate the performance and characteristic of the vectored thruster more fully, some influential factors of the duct have been considered in the following simulations. Referring to (19)- (21), the parameters of the duct are important to calculate in the numerical dynamic simulation. In this paper, we use CFD simulation to get the dynamics parameters of the duct. Figure 19 shows the drag, lift, and moment of the duct with the angle of attack range between 0°and 20°and maintains the flow velocity at 2.5 m/s, respectively.
Based on the dynamics parameters of the duct, as shown in Figure 19, the force and moment generated by the duct are calculated and presented in Figure 20.
e calculation results show that the influential factors of the duct play a very important role in calculating the driving forces of actuators; hence, the influential factors of the duct should be considered in a calculation schedule.
To improve the performance of the vehicle motion control, the tilt angles α and β play a crucial role in the vectored thruster AUVS. With considering the usage of the vectored thruster, a control scheme using the PID method is developed for the designed thrust-vectoring mechanism, as shown in Figure 21. PID algorithm is the most widely used control methods in all kinds of application fields for its effectiveness and practicability. Based on the control scheme shown in Figure 21, a control model of the proposed thrustvectoring mechanism is developed by Matlab and Simscape Mutibody, as shown in Figures 22 and 23.
To further investigate the performance of proposed control model, related numerical simulations are carried out with top platform moves according to � α s and β � 0. e related parameters are given in Table 2, the length responses of linear actuators can be obtained, as shown in Figure 24. e simulation results from Figures 14 and 24 show that the designed PID controller for the thrust-vectoring mechanism is fast, effective, and able to achieve the expected goal commendably. Based on the designed controller and considering the influence of factors as analyzed above, which includes extra forces and moments from duct and buoyant forces of the vectored thruster, the driving force of linear actuators is presented in Figure 25.
Comparing the driving force in Figures 18 and 25, the driving forces of actuators have some similarities in change trend at the same time, but the maximum and minimum of driving forces are distinctly different. Hence, it is concluded that the lift, drag, and torque of the duct propeller are important influential factors to the dynamic model of the vectored thruster. Accordingly, for the purposes of optimizing the structure and decreasing the dimension, it is of significance to choose the appropriate linear actuator for driving the proposed thrust-vectoring mechanism by studying the dynamic analysis with considering all kinds of influence factor.

Conclusion
In this paper, a design scheme for the vectored thruster based on a 3RPS parallel manipulator is proposed to solve the effect of the control surface weakening problems. Parallel manipulators have several advantages over the mechanical structure and are suitable for various application fields, such as compact and reasonable structure, fast response, and high positioning precision. Because of the merits of itself, this type mechanical structure is used to design the thrustvectoring mechanism considering the limited space and hash environment. Additionally, a duct propeller is adopted   as power source of the proposed vectored thruster, which is installed on the top rotating platform with a main motor as a whole structure, thereby this design ensures compact structure, reliable motion, and high propulsive efficiency. Because the control force is provided by the component force of thrust rather than common rudders, the proposed vectored thruster AUVs have the abilities to complete all kinds of certain tasks and operations at a low forward speed.
In order to make sure the designed vectored thruster can run efficiently and stably, studying and developing the control system is fundamental to implement the design function of the vectored thruster. Owing to the importance of control system, related theoretical research about kinematics and dynamics of the thrust-vectoring mechanism is carried out to establish the motion model. In the kinematic analysis, the inverse and forward kinematics of the thrustvectoring mechanism is presented, and the numerical simulations are accomplished to prove the feasibility and effectiveness of this design. In the section of workspace analysis, the study of the tilt angles of revolute and spherical joints is also carried out to make sure the motion platform can implement its designed function in limited motion space. In order to further verify feasibility of the thrustvectoring mechanism, based on the considerations of various affecting factors, a dynamics model of the designed thrust-vectoring mechanism is established according to theoretical analysis, and the driving forces of the linear actuator are presented though a series of numerical simulations. In addition, a control scheme based on PID algorithm is proposed for the thrust vector control system on the existing work basis, and a control model is established using Simsacpe Mutibody, and the simulation results proved the feasibility of the proposed control scheme, which can effectively realize the goal of controlling the thrust-vectoring mechanism.
According to the above, the designed vectored thruster is able to provide the vectored thrust effectively and efficiently, and the AUVs equipped with vectored thruster are able to complete a variety of the complex tasks at a comparably low forward speed.
In the future research, a series of numerical simulations and theoretical study are carried out to investigate hydrodynamic performance of this vectored thruster AUV. On this basis, a prototype of this designed vectored thruster will be developed and experimental test will be carried out to verify the principles design. Moreover, the corresponding control system of the vectored thruster as a part of the AUV will be developed and tested in pools or open water to check its performance.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.