the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.

A control strategy for shipboard stabilization platforms based on a fuzzy adaptive proportional–integral–derivative (PID) control architecture
Hongbo Liu
Zhaobiao Zeng
Xiaodong Yang
Yuxin Zou
Xianli Liu
To address the precision degradation of marine equipment under coupled hydrodynamic disturbances, this study develops a 6-degree-of-freedom (6-DOF) stabilization platform with a fuzzy adaptive proportional–integral–derivative (PID) control architecture. The kinematic model is established via analysis based on the virtual-work principle, complemented by Monte Carlo simulations for workspace characterization. A fuzzy inference engine dynamically adjusts PID parameters through rule-based adaptation, demonstrating superior disturbance rejection. Comparative simulations indicate a 50 % reduction in settling time (7.0 s to 3.5 s), zero overshoot, and < 0.03° steady-state tracking error under 2 Hz sinusoidal excitation. A human–machine interface (HMI) for the shipboard stabilization platform is developed using the Qt Creator framework, integrating real-time trajectory tracking and parameter tuning. The research advances marine stabilization technology through mechanical optimization via virtual-work modeling and control enhancement via fuzzy–PID synthesis. Experimental validation confirms the framework's capability to maintain sub −0.03° precision under dynamic maritime conditions.
- Article
(11932 KB) - Full-text XML
- BibTeX
- EndNote
Under combined marine wind and wave effects, ships exhibit complex multi-degree-of-freedom motions, including roll, pitch, and yaw (RPY), significantly compromising the operational accuracy of onboard equipment (Tu et al., 2022; Qiang et al., 2024). Current research on shipboard stabilization platforms focuses on structural design and control strategies (Chen et al., 2024; Liu et al., 2025; Mei et al., 2021).
In structural design, novel parallel mechanisms dominate recent advancements. Liu et al. (2023) developed a 6-degree-of-freedom (6-DOF) parallel platform with motion prediction compensation, achieving < 0.7° average attitude tracking error under Sea State IV. Tang et al. (2023) proposed a hybrid serial-parallel platform that extends attitude adjustment range by 35 % through additional serial modules. Zhang et al. (2024) established an inverse kinematics model based on a 6-UCU (where U denotes universal joint and C denotes cylinder joint) structure, demonstrating a ± 5° disturbance rejection capability in dynamic simulations. Notably, Han et al. (2023)'s 3-UPS/S (where U denotes universal joint, P denotes prismatic joint, and S denotes spherical joint) mechanism achieves a millisecond-level response time in continuous forward kinematics analysis, while Koraaa et al. (2024)'s 2-UPU (where U denotes universal joint and P denotes prismatic joint) platform verifies ± 1.2 mm actuator stroke accuracy via Simscape multibody modeling. Liu et al. (2024)'s 6-PRSS (where P denotes prismatic joint, R denotes revolute joint, and S denotes spherical joint) structure exhibits a 0.85 m3 effective workspace in multi-vibration-mode tracking, and Pashkov (2024)'s modified Gough–Stewart platform extends the linear motion range by 15 %.
Regarding control strategies, advanced algorithms continuously enhance performance (Ya et al., 2023a, b; Lei et al., 2024). Zhao et al. (2022)'s fractional-order active disturbance rejection control reduces stabilization time by 30 %, while Vu et al. (2022)'s fuzzy sliding-mode control decreases actuator saturation by 40 %. Zhou et al. (2023)'s hybrid equivalent-input-disturbance and sliding-mode method improves disturbance rejection by 25 %, while Gong et al. (2023)'s fault-tolerant control maintains < 0.5° attitude error under 20 % actuator failure. Chen et al. (2023)'s triple-loop active disturbance rejection control (ADRC) strategy reduces power consumption by 18 % while enhancing decoupling performance, and Lv et al. (2024)'s dynamic gravity-compensated PD control decreases 3-RCU (where R denotes revolute joint, C denotes cylindrical joint, and U denotes universal joint) platform sway amplitude by 62 %.
Despite these achievements, two challenges persist: (1) most control strategies rely on offline simulation without real-time interaction with multibody dynamic models; (2) parameter tuning depends on empirical trial and error and lacks online adaptive mechanisms under environmental disturbances. Our study addresses these gaps through ADAMS–MATLAB co-simulation for high-fidelity digital-twin construction, a fuzzy adaptive proportional–integral–derivative (PID) algorithm for dynamic parameter optimization, and a visualization interface based on Qt Creator. Experimental results demonstrate a 42 % faster response time and a 58 % lower steady-state error compared to conventional PID, providing a novel, real-time control solution for shipboard stabilization systems.
To mitigate wave-induced disturbances in maritime transportation, this study developed a 6-DOF parallel mechanism. The workspace was geometrically characterized through Monte Carlo simulations, supported by forward kinematics solutions. Dynamic equations were systematically derived using virtual-work principles. Experimental validation confirmed close agreement between measured platform trajectories and simulated actuator forces, demonstrating the model's efficacy in maritime motion compensation.
2.1 Kinematic modeling
The kinematic modeling begins with coordinate system establishment:
(1.) Multibody reference frames (Fig. 1) were defined as
(2.) Coordinate transformations between frames were constructed using
where PG is the absolute position in {G}, is the vessel-to-Earth rotation matrix (), is the platform-to-vessel rotation matrix , and is the translation vector from {G} to {E} origin.
(3.) Critical parameters, including geometric dimensions and mass properties, are listed in Table 1, and the coordinate relationship symbols are systematically defined in Table 2.
As shown in Fig. 2 for the distribution of the hinge position of each key connection point of the shipborne stabilized platform, the position of the origin of the {F} coordinate system under the {E} coordinate system is .
Table 2Symbol definitions and coordinate relationships.

Note: All coordinate systems follow the right-hand rule. Rotation matrices comply with .
The position of the upper platform center of mass in {F} is {FCup} and its position under the {E} system is {ECup}, which can be calculated by
The distribution of the hinge nodes of the shipborne stabilized platform is shown in Fig. 3. The position of the ith hinge point of the upper platform under the {F} coordinate system can be expressed as FSi=(FSixFSiyFSiz)T. The position of the ith hinge point of the platform on the shipborne stabilized platform can then be expressed in the {E} coordinate system by
The position of the gimbal center point of the ith leg of the lower platform in the {E} coordinate system can be expressed as EUi=(EUixEUiyEUiz)T. The ith leg length vector is , which is obtained by solving for its upper gimbal center and lower gimbal center points in the {E} coordinate system.
From the above rotation matrix, we obtain
2.1.1 Kinematic modeling in non-inertial systems
In comparison with the previous section, the {E} coordinate system is a local coordinate system that moves with the inertial coordinate system. Thus, the center-of-mass position of the platform on the shipborne stabilized platform in the {G} coordinate system is as follows:
where GOE represents the coordinate position of the origin of the {E} coordinate system in the {G} coordinate system and ECup represents the position of the upper platform under the {E} coordinate system. The solution based on the above positions' yields the following:
Both sides of Eq. (7) are differentiated to first order in time, and the velocity of the center of mass of the upper platform in the {G} coordinate system is solved as follows:
where represents the velocity of the origin of the {E} coordinate system in the {G} coordinate system, represents the angular velocity of the {E} coordinate system in the {G} coordinate system, and represents the velocity of the upper platform under the {E} coordinate system.
The acceleration of the center of mass of the upper platform in the {G} coordinate system is solved by second-order differentiation against time on both sides:
where represents the acceleration of the origin of the {E} coordinate system in the {G} coordinate system and represents the angular acceleration of the {E} coordinate system in the {G} coordinate system. The angular velocity of the upper platform in the {G} coordinate system, i.e., the angular velocity of the {F} coordinate system in the {G} coordinate system, is as follows:
The angular acceleration of the upper platform in the {G} coordinate system is solved for both sides of Eq. (10) for first-order differentiation in time:
where represents the angular acceleration of the center of mass of the upper platform in the {G} coordinate system. The position of the center of the first leg in the {G} coordinate system is as follows:
The first-order differentiation of time on both sides of Eq. (12) is solved for the velocity of the center of the first leg in the {G} coordinate system:
The acceleration of the center of the ith leg in the {G} coordinate system is solved for the second-order differentiation of time on both sides of Eq. (13):
The position of the piston center of mass in the ith leg in the {G} coordinate system is as follows:
The first-order differentiation of time on both sides of Eq. (15) is solved for the velocity of the piston center of mass of the first leg in the {G} coordinate system:
The angular velocity of the ith leg of the shipborne stabilized platform in the {G} coordinate system is solved as follows:

Figure 7Output space and initial position comparison diagram: (a) output posture example 1; (b) output posture example 2.
The angular velocity of the ith leg in the {G} coordinate system is solved for both sides of Eq. (17) for first-order differentiation in time:
The acceleration of the piston center of mass of the ith leg in the {G} coordinate system is solved for the first-order differentiation of time on both sides of Eq. (16):
The position of the cylinder center of mass at the ith branch angle in the {G} coordinate system is as follows:
Both sides of Eq. (20) are differentiated to the first order in time, and the velocity of the cylinder center of mass of the ith leg in the {G} coordinate system is solved as follows:
The acceleration of the cylinder center of mass of the ith leg in the {G} coordinate system is solved for the second-order differentiation of time on both sides of Eq. (21):
2.1.2 Workspace parameter analysis
In this paper, a stochastic probabilistic algorithm based on a numerical method, the Monte Carlo simulation method, is used. It is particularly noteworthy that the error of this method is not associated with the dimensionality of its system and can be solved directly for problems containing statistical properties in the system, while no discretization is required for continuous problems (Hunek et al., 2022; Jing et al., 2024).
For the randomness of the Monte Carlo method, we use the rand () function for the representation, followed by the random number selection of the joint angles of the shipborne stabilized platform in its angle range. The code is as follows: . The meaning of this code is as follows: taking the angle range as the minimum value is used as the basis, and the random number in the angle range is added to get the current joint random value of this axis. Finally, the end position is obtained by the kinematic script and plotted by the plot3() function to obtain the approximate workspace. The number of steps set in this practical example is 40 000, whether it is a scatterplot or an indicator of the workspace, the accuracy of which depends on the pre-set step size. The Monte Carlo simulation method is more convenient for the spatial evaluation of the workspace, as shown in Fig. 3 for the graph of the results of this fixed-attitude workspace solution.

Figure 8Realization of the 6-DOF shipborne stabilized platform. The required degrees of freedom are obtained as follows: (a) translate in the X direction, (b) translate in the Y direction, (c) translate in the Z direction, (d) rotate around the X axis, (e) rotate around the Y axis, and (f) rotate around the Z axis.
2.1.3 Stabilized platform parametric interface simulation
Based on the above kinematic model of the shipborne stabilized platform, the simulation program is established, and the parametric motion simulation of the shipborne stabilized platform is realized through the user interface. The parameters of the stabilized platform are input, including the outer circumferential radius of the platform and the distribution angle of the platform hinge point; the motion trajectory of the platform is filled in; and the position of the inverse solution to get the rod length position curve and data of the platform is selected. The inverse solution results are shown in Fig. 4.
Similarly, by inputting the machine parameters and motion trajectory and then selecting the velocity inverse solution, the rod length velocity curve and data of the drive joint can be obtained; these output results are shown in Fig. 5. After inputting the parameters and selecting the acceleration inverse solution, the rod length acceleration curve and data of the drive mechanism can be obtained; these output results are shown in Fig. 6.
For an input yaw angle of 15°, a pitch angle of 20°, and a roll angle of 25°, the output spatial position compared with the initial position is shown in Fig. 7a. For an input yaw angle of −5°, a pitch angle of −15°, and a roll angle of −10°, the output spatial position compared with the initial position is shown in Fig. 7b.
2.2 Kinetic analysis
Based on D'Alembert's theorem, the inertial force acting on each moving part of the platform on the shipborne stabilized platform can be expressed as
where mu represents the mass of the platform on the shipborne stabilized platform.
The gravitational force applied to the upper platform can be expressed as
where g represents the acceleration of gravity.
The upper-platform inertia force can be expressed as
where FIu represents the inertia tensor of the upper platform at {G}.
where FIu represents the inertia tensor of the upper platform at {F}, . The inertia force of the ith leg piston is as follows:
where miu denotes the mass of the ith foot piston.
The gravity of the ith leg piston can be expressed as
The moment of inertia of the ith leg piston can be expressed as
where GIiu represents the inertia tensor of the ith foot piston at {G}.
where represents the inertia tensor of the first leg piston under {Qi} and represents the rotation matrix from {G} to {Qi} ().
The inertia force of the ith leg cylinder is as follows:
where mil represents the moment of inertia of the ith leg cylinder.
where GIil represents the inertia tensor of the ith foot sleeve under {G}.
where represents the inertia tensor of the ith leg under the cylinder.
2.3 Platform performance analysis
The experimental framework established 12 gimbal centers as coordinate origins, with corresponding components systematically integrated: gimbal vices positioned at each pivot node, cylindrical joints connecting motion rods to outer rails, and a central-point drive mechanism mounted on the upper platform (Fig. 8). Parameter initialization involved a two-phase simulation protocol – commencing with a 35 s dynamic analysis encompassing 4000 computational steps for platform trajectory validation, followed by a 9.36 s precision simulation executing 982 discrete steps for actuator force quantification. The implemented drive functions demonstrated 6-DOF motion compliance (Fig. 9) and enabled directional vector analysis of actuator forces (Fig. 10).
Simulation outcomes substantiated the virtual prototype's efficacy through three critical performance indicators: rapid response capability, evidenced by sub-15 s motion cycle completion across all translational axes (Fig. 9); smooth force transitions, maintaining continuous derivative profiles in all actuators (Fig. 10); operational stability, demonstrated by force curves exhibiting monotonic variation patterns. This validation establishes reliable foundational support for marine stabilization platform co-simulation studies, confirming effective coordination between dynamic responsiveness and operational steadiness.
Based on the virtual-prototype model of the shipboard stabilization platform established above, data exchange of state variables is performed to realize the joint simulation and analyze the system at the same time.
3.1 Construction of the co-simulation module
The co-simulation process of the shipborne stabilized platform is divided into three steps. The first is model building, which consists of three main aspects: geometry modeling, element division, and attribute definition. Secondly, the input and output interfaces of the built model are set up, determining the inputs and outputs. Finally, the simulation parameters are set to carry out the joint simulation, plot the simulation curves, and solve the simulation results.
3.2 Design of the control programs
The model is first initialized with parameters, and the simulation software is then interconnected, as shown in the selection window of Simulink in Fig. 11. The S-Function module encapsulates the nonlinear model of the shipborne stabilized platform developed in the previous section, while the State-Space module specifically refers to the linearized state–space representation derived from the original nonlinear multibody dynamics model through ADAMS-software-based linearization. This linearized model serves two critical purposes:
-
simplify the controller design by reducing model complexity;
-
enable joint simulation verification, particularly under conditions involving small disturbances.
The Adams_sub module contains the complete nonlinear equations and associated variables, with the S-Function module providing the interface for co-simulation. As shown in Fig. 12, this dual-model architecture integrates both representations to achieve a comprehensive dynamic analysis capability.
In the co-simulation, traditional frequency-domain analysis methods based on transfer functions (e.g., Bode plots and root locus) are challenging to apply directly, due to the nonlinear multibody dynamics model of the platform (ADAMS–Simulink coupling). Consequently, the initial PID parameter tuning adheres to a stability-first principle, employing empirical tuning methods: the baseline proportional gain is determined through the Ziegler–Nichols critical gain method based on step response characteristics, while integral and derivative times are fine-tuned through trial and error to balance dynamic response speed and overshoot suppression.
The PID parameters in Fig. 13 were optimized through Ziegler–Nichols tuning and fuzzy adaptation, governed by motor constants, control laws, hierarchical refinements, and results.
Motor constants were as follows: Ke = 2.75 N m A−1, Ra = 9.5 Ω, Ta=1.3 ms, Km = 1.52 V (rad s)−1, and Tm = 0.69 s.
Control laws were as follows: Kp⋅sign(e), , and Kd⋅d[sign(e)]dt.
Hierarchical refinements were as follows:
-
baseline from Ziegler–Nichols critical gain,
-
dynamic compensation via = 0.454 s2 rad−1,
-
49-rule fuzzy matrices ().
Results were as follows: a settling time of 0.83 s (1.2Tm) under a = 0.29 V A−1 constraint, validated through co-simulation with nonlinear dynamics.
In this system, the dynamic equations of the motor can be expressed as
where E represents the induced electromotive force during the dynamic process, u represents the voltage during the dynamic process, ia represents the current during the dynamic process, and ra represent armature resistance.
where Ce is the torque constant and n is the rotor speed (rpm).
The electromagnetic torque can be expressed as
where Cm represents the electric potential speed ratio at rated excitation.
where GD2 represents the flywheel torque of the motor (N m−1) and TL represents the load torque.
By applying the Laplace transform to the time-domain equations, we can derive the frequency-domain representation of the motor system and, ultimately, obtain its transfer function, which can be expressed as follows:
where T1 represents the armature circuit electromagnetic time constant and Tm represents the power drag electromechanical time constant. Other values in the above expression are as follows: , ra = 9.5 Ω, T1 = 1.3 × 10−3 s, Ce = 1.52 , and Tm = 0.69 s.
The Ziegler–Nichols method was applied to rectify the initial parameters of the fuzzy controller. Tm = 0.69 s, and the fact that the stabilizing platform itself is heavy indicates that the system has high inertia and responds slowly. Wave disturbances require the platform to respond within milliseconds. A higher Kp value can dramatically increase the sensitivity of the system to errors and shorten the adjustment time. T1 = 1.3 × 10−3s indicates a rapid current response. The error accumulation may cause the integral term to go out of control when the shipboard platform is activated or significantly perturbed. A lower Ki value limits the integration speed and avoids oscillations. High-frequency noises, such as wave impacts and mechanical vibrations, need to be passed through the high Kd fast suppression. The differential action is sensitive to the rate of change in the error and can correct the system behavior in advance.
Based on the Ziegler–Nichols method of calculation and then adjusted for the work requirements of the stabilized flat (the final choice of Kp0 = 20 000, Ki0 = 100, and Kd0 = 260), the corresponding fuzzy subsets are chosen as . The comparative test in Fig. 19 also verifies the validity of the initial parameter selection.
Based on the joint simulation and system design, the shipboard stabilization platform was first optimized, and the fuzzy PID controller was then designed and simulated to complete the correction of the shipboard stabilization platform control system. Finally, the research results were visualized using Qt Creator front-end software.

Figure 16The input membership functions of the fuzzy controller: (a) membership function of the system error e; (b) membership function of the system error change rate ec.
4.1 Fuzzy adaptive PID control design and simulation
The disadvantage of PID control is obvious, the controller parameters are difficult to select (Fan and Teng, 2024). The parameters have a great impact on the performance of the controller, and errors in the parameters may directly lead to system instability and control failure. When encountering nonlinear and strongly coupled systems, the model cannot be recognized, and the parameters cannot be self-regulated or adjusted in real time. Fuzzy adaptive control does not require precise design of the control object and has significant advantages over traditional PID control (Deng et al., 2022). Therefore, the fuzzy adaptive control dynamically adjusts the parameters of the PID controller (Chang et al., 2024) according to the created fuzzy rules in order to make the system more stable and, finally, compares the control effects of the two control methods according to the simulation results.
As shown in Fig. 14, in the control strategy of the control system, the platform angle deviation detected by the sensor is recorded as the system error e. The system error e and its change rate ec are used as the inputs of the fuzzy control system and are processed to obtain the quantization factors Ke and Kec, respectively. Ke and Kec are fuzzified utilizing the subordinate function library, and fuzzy reasoning is then carried out to obtain the fuzzy outputs according to the control rule library. This is followed by the clarification process, via which we finally obtain the proportional factors ΔKp, ΔKi, and ΔKd, through the proportional factors to the PID control parameters for real-time dynamic adjustment. Finally, according to the fuzzy PID controller output value, the decoupling operation is used to solve the length of each rod, to realize the stabilization and compensation effect on the platform.
4.2 Design of the fuzzy PID controller
The input variables of the systematic error and its rate of change are determined under the quantization factor. Based on the influence of PID parameters on the system control, the fuzzy rules of ΔKp, ΔKi, and ΔKd are respectively designed by combining the trial-and-error method and the design experience of PID controllers, as shown in Tables 3–5.
With respect to the effect of Kc on the system performance, the larger the Kc, the more the interval shrinks and enhances the control effect. Therefore, with a larger selection, the overshoot of the system decreases and the response speed decreases. Kc has a significant effect on the suppression of system overshoot.

Figure 17The output membership functions of the fuzzy controller: (a) membership function of the proportional constant ΔKp; (b) membership function of the integral constant ΔKi; (c) membership function of the differential constant ΔKd.
The magnitude of Ke greatly affects the dynamic characteristics of the system. The control interval of the fuzzy controller is the interval of variation in the system error e(t) when Ke=1. For Ke<1, the interval of error e(t) expands, leading to a decrease in the sensitivity of the fuzzy controller to the input and, thus, weakening the control of the offset. For Ke>1, the interval of error e(t) contracts, leading to an increase in the sensitivity of the fuzzy controller to the input and, thus, enhancing the control interval of the offset. The rising rate and overshoot of the system increase for larger Ke choices, increasing the transition time. When the output controller output is the same, if Ku is selected to be small, it will lead to a longer dynamic response time of the system. Furthermore, the larger the Ku, the better the controller control effect; the system rings quickly and is easy to overshoot. However, if Ku is selected to be too large, the system produces an oscillation.

Figure 21Step tracking and error analysis of fuzzy PID and PID: (a) step response comparison; (b) step tracking error comparison.

Figure 22Sine tracking and error analysis of fuzzy PID and PID: (a) sine response comparison; (b) sine tracking error comparison.
The output variables Kp, Ki, and Kd of the fuzzy controller are solved by fuzzy inference and inverse fuzzification of the input variables e and ec. The fuzzy output value is taken as the output variable, Kp, Ki, or Kd, corresponding to the fuzzy value in its theoretical domain, and its scale factor is determined and multiplied with the output of the fuzzy controller to obtain the increment of the PID parameters Kp, Ki, and Kd. The control variables of the system are updated to produce a suppression effect on the system, and the corresponding adjustment algorithm is as follows:
where Kp0, Ki0, and Kd0 are the initial values of Kp, Ki, and Kd , respectively.
4.2.1 Simulation of fuzzy adaptive PID control system
The model diagram of the fuzzy adaptive PID control system is shown in Fig. 15.
On the base of fuzzy set theory, the input and output data should be transformed into proper variables with details as follows:
Input variables.
-
system error e, domain of , and scaling factor Ke=100;
-
system error rate ec, domain of ° s−1, and scaling factor Kc=100.
Output Variables.
-
proportional gain increment ΔKp, domain , and scaling factor K1=0.1;
-
integral time increment Ki, domain , and scaling factor K2=0.1;
-
derivative time increment ΔKd, domain , and scaling factor K3=0.1.
This design enhances the sensitivity to dynamic disturbances through an expanded scaling factor while suppressing overshoot, thereby ensuring controller robustness under time-varying sea conditions. The input and output membership functions are shown in Figs. 16 and 17.
Once the affiliation functions of the two inputs and three outputs in the system are set, the control rules of the function will be edited as outlined in the following. The 49 rules are obtained through the fuzzy control table in the previous section. Figure 18 shows the fuzzy rule editing interface, and Fig. 19 shows the interface for defining the inputs and outputs.
From the fuzzy rules above, the surface relations of Kp, Ki, and Kd under fuzzy rules can be obtained in the setting of the fuzzy controller, as shown in Fig. 20.
The control block diagram was established using Simulink, and the membership functions and rule base were input into the fuzzy controller. Figure 21 illustrates the step response and tracking error comparison between the fuzzy control system and the conventional PID control system, further revealing performance differences under varying signals.
Under a step input, the traditional PID exhibits a settling time of 7 s. Due to fixed parameter limitations, although the steady-state error approaches zero, significant overshoot (approximately 15 %) occurs during regulation (Fig. 21a), resulting in larger tracking errors in the transient phase (Fig. 21b). In contrast, the fuzzy adaptive PID dynamically adjusts parameters, reducing the settling time by 50 % and driving the tracking error close to zero without overshoot (Fig. 21b). This demonstrates the fuzzy PID's superior error suppression capability during dynamic phases, rapid stabilization, and enhanced control precision for step signals.
As shown in Fig. 22, the dashed line represents a sine signal with an amplitude of 1 and a frequency of 2 Hz, while the solid line denotes the actual output. The traditional PID, constrained by fixed parameters and limited bandwidth, accumulates high-frequency errors, leading to phase lag and steady-state error peaks in the tracking curve. The fuzzy PID, however, employs online parameter optimization to compensate for phase and amplitude deviations in real time. Its membership function design effectively suppresses periodic disturbances, resulting in near-perfect alignment between the tracking curve and the target signal. The steady-state error peak is controlled within 0.03°, with no amplitude attenuation, demonstrating robust performance and high real-time adaptability even under time-varying signals.
The fuzzy adaptive PID outperforms traditional PID in both step and sinusoidal signal scenarios, particularly under complex operating conditions (e.g., time-varying loads and nonlinear systems). While traditional PID relies on empirical parameter tuning, the fuzzy PID reduces manual intervention costs through its adaptive mechanism. Comparative analysis of tracking and error characteristics for step and sinusoidal signals confirms that the fuzzy adaptive PID significantly surpasses conventional PID with respect to settling time, overshoot suppression, steady-state accuracy, and phase compensation. Its multi-objective optimization capability provides an advanced solution for high-precision control systems.
4.3 User interface design for fuzzy control simulation of a shipboard stabilization platform
This interface is built mainly for the fuzzy control module; therefore, one can freely edit e and ec input and Kp, Ki, and Kd output in the interface. The main interface is shown in Fig. 23.
The user can click on the “Fuzzy rule” button to freely set the rules of fuzzy control. This rule interface greatly enhances input efficiency. The user can then click on the “OK” button to save the fuzzy control rules, as shown in Fig. 24.
According to the designed fuzzy control rules, the user can click on the “Plot” button, and the program will output the response curve of the set rules. The reference of the program itself comes with a PID control response curve.
A novel 6-DOF shipboard stabilization platform integrated with a fuzzy adaptive PID control was proposed to address the challenges involved with precision equipment operation in dynamic marine environments. The main contributions of this work are summarized as follows:
-
To address the issues of limited workspace and low modeling accuracy in traditional ship stabilization platforms, a three-axis 6-DOF parallel mechanism was developed, which expands workspace and enables compensation for combined ship motions (roll, pitch, and yaw) through optimized hinge distribution with a workspace verified via Monte Carlo simulations. In addition, a dynamic modeling method for a non-inertial system based on the virtual-work principle is proposed, which solves the coupling error between the hull motion and the platform attitude.
-
The proposed fuzzy adaptive PID controller achieved a 50 % reduction in settling time (7.0 to 3.5 s) and eliminated overshoot completely, with steady-state error constrained below 0.03° under 2 Hz sinusoidal excitation. The fuzzy system provided a scalable framework for marine control applications. The fuzzy adaptive PID controller demonstrates significant superiority over the conventional PID in terms of settling time reduction, overshoot suppression, steady-state accuracy, and phase compensation. Finally, a ship stabilization platform control user interface was designed based on Qt Creator and enabled real-time motion trajectory visualization and interactive fuzzy rule configuration.
The data that support the findings of this study are available from the corresponding author upon reasonable request.
HL, ZZ, and YZ were involved in writing the first draft of the paper, HL, ZZ, and XY were engaged in the simulation analysis of the paper. HL and ZZ participated in revising and embellishing the paper. XL completed the final review of the paper.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
The authors would like to thank editor-in-chief, editor, and anonymous reviewers for their valuable reviews.
This research has been supported by the National Natural Science Foundation of China (grant no. 51805121).
This paper was edited by Zi Bin and reviewed by two anonymous referees.
Chang, S., Cao, J., Pang, J., Zhou, F., and Chen, W.: Compensation control strategy for photoelectric stabilized platform based on disturbance observation, Aerosp. Sci. Technol., 145, 108909, https://doi.org/10.1016/j.ast.2024.108909, 2024.
Chen, P., Jiang, M., Lv, X., Zhou, H., Yang, D., Zhou, Y., Jin, Z., and Peng, S.: Research on the application of inertially stabilized platform in the dynamic measurement of cold atomic gravimeter, Heliyon, 10, e23936, https://doi.org/10.1016/j.heliyon.2023.e23936, 2024.
Chen, W., Wang, S., Li, J., Lin, C., Yang, Y., Ren, A., Li, W., Zhao, X., Zhang, W., Guo, W., and Gao, F.: An ADRC-based triple-loop control strategy of ship-mounted Stewart platform for six-DOF wave compensation, Mech. Mach. Theory, 184, 105289, https://doi.org/10.1016/j.mechmachtheory.2023.105289, 2023.
Deng, J., Xue, W., Zhou, X., and Mao, Y.: On Dual Compensation to Disturbances and Uncertainties for Inertially Stabilized Platforms, Int. J. Control Autom., 20, 1521–1534, https://doi.org/10.1007/s12555-021-0022-3, 2022.
Fan, L., and Teng, L.: Design and kinematics analysis of 3-DOF ship-borne parallel stabilization platform, Ship Science and Technology, 46, 71–74, https://doi.org/10.3404/j.issn.1672-7649.2024.11.013, 2024.
Gong, Q., Liu, Z., Hu, X., Teng, Y., Li, K., and Han, G.: Robust anti-disturbance fault-tolerant control of ship-board platforms with multiplicative actuator faults and unknown disturbances, Ocean Eng., 286, 115552, https://doi.org/10.1016/j.oceaneng.2023.115552, 2023.
Han, B., Jiang, Y., Yang, W., Xu, Y., Yao, J., and Zhao, Y.: Kinematics characteristics analysis of a 3-UPS/S parallel airborne stabilized platform, Aerosp. Sci. Technol., 134, 108163, https://doi.org/10.1016/j.ast.2023.108163, 2023.
Hunek, W. P., Majewski, P., Zygarlicki, J., Nagi, Ł., Zmarzły, D., Wiench, R., Młotek, P., and Warmuzek, P.: A Measurement-Aided Control System for Stabilization of the Real-Life Stewart Platform, Sensors, 22, 7271, https://doi.org/10.3390/s22197271, 2022.
Jing, C., Ma, X., Du, H., Liu, Y., Liu, C., and Hui, Y.: Adaptive Robust Control for Six-Degree-of-Freedom Electro-Hydraulic Parallel Ship Stabilization Platform, Ijomam, I, 1–12, https://doi.org/10.17683/ijomam/issue18.12, 2024.
Koraaa, M. M., Ghamry, K. A., Abo-Elnor, M. E., and Elsherif, I. A.: Inverse Kinematics Analysis of a 2-DOF Stabilized Platform using Simulink and Simscape, J. Phys. Conf. Ser., 2811, 012003, https://doi.org/10.1088/1742-6596/2811/1/012003, 2024.
Lei, H., Tongzhen, D., Lingbin, S., Qiang, J., and Lifu, W.: Structural Design Key Points and Stability Analysis of Shipborne Self-Stabilizing Platform, in: 2024 4th Asia–Pacific Conference on Communications Technology and Computer Science (ACCTCS), 24–26 February 2024, Shenyang, China, 690–695, https://doi.org/10.1109/ACCTCS61748.2024.00128, 2024.
Liu, J., Cai, C., Song, D., Zheng, J., Guo, X., Sun, X., Hu, Z., and Li, Q.: Power output and platform motion regulation of floating offshore wind turbine using model predictive control based on coupled dynamics model, Ocean Eng., 322, 120451, https://doi.org/10.1016/j.oceaneng.2025.120451, 2025.
Liu, Y., Yuan, H., Xiao, Z., and Xiao, C.: An Offshore Self-Stabilized System Based on Motion Prediction and Compensation Control, Journal of Marine Science and Engineering, 11, 745, https://doi.org/10.3390/jmse11040745, 2023.
Liu, Y., Wang, Q., Qiu, J., Huang, S., and Tao, W.: A 6 Degrees of Freedom Platform for Adjusting Pose Based on 6-PRSS, IJST-TMech. Eng., 49, 897–916, https://doi.org/10.1007/s40997-024-00802-w, 2024.
Lv, Z., Liu, P., Ning, D., and Wang, S.: Anti-Swaying Control Strategy of Ship-Mounted 3-RCU Parallel Platform Based on Dynamic Gravity Compensation, Machines, 12, 209, https://doi.org/10.3390/machines12030209, 2024.
Mei, D. and Yu, Z.-Q.: Disturbance rejection control of airborne radar stabilized platform based on active disturbance rejection control inverse estimation algorithm, Assembly Autom., 41, 525–535, https://doi.org/10.1108/AA-10-2020-0158, 2021.
Pashkov, I. V.: Modification of the Gough–Stewart Platform with Six Degrees of Freedom, J. Mach. Manuf. Reliab., 53, 307–313, https://doi.org/10.1134/S1052618824700237, 2024.
Patra, B., Nag, A., and Bandyopadhyay, S.: Analytical determination of the optimal effective regular workspace of a 6-6 Stewart platform manipulator for a specified orientation workspace, Mech. Mach. Theory, 203, 105791, https://doi.org/10.1016/j.mechmachtheory.2024.105791, 2024.
Qiang, H., Du, L., Liu, K., Zhou, L., and Zeng, S.: Kinematics and Dynamics Modeling of 1T2R Heavy-load Parallel Stabilized Platform with Analytical Solution, Acta Armamentarii, 45, 3960–3969, https://doi.org/10.12382/bgxb.2023.1168, 2024.
Tang, G., Lei, J., Li, F., Zhu, W., Xu, X., Yao, B., Claramunt, C., and Hu, X.: A modified 6-DOF hybrid serial–parallel platform for ship wave compensation, Ocean Eng., 280, 114336, https://doi.org/10.1016/j.oceaneng.2023.114336, 2023.
Tu, Y., Du, J., and Liu, W.: Design and kinematics analysis of 3-DOF ship-borne parallel stabilization platform, Ship Science and Technology, 44, 155–157, https://doi.org/10.3404/j.issn.1672-7649.2022.05.033, 2022 (in Chinese).
Vu, M. T., Alattas, K. A., Bouteraa, Y., Rahmani, R., Fekih, A., Mobayen, S., and Assawinchaichote, W.: Optimized Fuzzy Enhanced Robust Control Design for a Stewart Parallel Robot, Mathematics, 10, 1917, https://doi.org/10.3390/math10111917, 2022.
Ya, J., Dong, R., and Yu, Z.: Control system of shipborne stable platform based on RBF neural network, Ship Science and Technology, 45, 181–185, https://doi.org/10.3404/j.issn.1672-7649.2023.09.040, 2023b.
Ya, T., Du, J., and Li, J.: PTESO based robust stabilization control and its application to 3-DOF ship-borne parallel platform, Ocean Eng., 284, 115202, https://doi.org/10.1016/j.oceaneng.2023.115202, 2023a.
Zhang, Y., Lu, D., Wang, Z., Wang, J., Pan, Q., and Su, S.: Modeling and simulation of six degrees of freedom marine stable platform, in: International Conference on Mechatronic Engineering and Artificial Intelligence (MEAI 2023), 15 December 2024, Shenyang, China, 78, https://doi.org/10.1117/12.3025491, 2024.
Zhao, J., Zhao, T., and Liu, N.: Fractional-Order Active Disturbance Rejection Control with Fuzzy Self-Tuning for Precision Stabilized Platform, Entropy, 24, 1681, https://doi.org/10.3390/e24111681, 2022.
Zhou, Y., She, J., Wang, F., and Iwasaki, M.: Disturbance Rejection for Stewart Platform Based on Integration of Equivalent-Input- Disturbance and Sliding-Mode Control Methods, IEEE/ASME Trans. Mechatron., 28, 2364–2374, https://doi.org/10.1109/TMECH.2023.3237135, 2023.