Next Article in Journal
Designing a Geodesic Faceted Acoustical Volumetric Array Using a Novel Analytical Method
Previous Article in Journal
Proposals for Surmounting Sensor Noises
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Differential Evolution Adaptive Fuzzy PID Control Method for Gravity Measurement Stable Platform

School of Electrical Engineering, Naval University of Engineering, Wuhan 430033, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(6), 3172; https://doi.org/10.3390/s23063172
Submission received: 21 February 2023 / Revised: 12 March 2023 / Accepted: 13 March 2023 / Published: 16 March 2023
(This article belongs to the Section Navigation and Positioning)

Abstract

:
In the platform gravimeter, the stabilization accuracy of the gravimetric stabilization platform is crucial to improve the accuracy of gravity value measurements due to its uncertainties, such as mechanical friction, inter-device coupling interference, and nonlinear disturbances. These cause fluctuations in the gravimetric stabilization platform system parameters and present nonlinear characteristics. To resolve the impact of the above problems on the control performance of the stabilization platform, an improved differential evolutionary adaptive fuzzy PID control (IDEAFC) algorithm is proposed. The proposed enhanced differential evolution algorithm is used to optimize the initial control parameters of the system adaptive fuzzy PID control algorithm to achieve accurate online adjustments of the gravimetric stabilization platform’s control parameters when it is subject to external disturbances or state changes and attain a high level of stabilization accuracy. The results of simulation tests, static stability experiments, and swaying experiments on the platform under laboratory conditions, as well as on-board experiments and shipboard experiments, all show that the improved differential evolution adaptive fuzzy PID control algorithm has a higher stability accuracy compared with the conventional control PID algorithm and traditional fuzzy control algorithm, proving the superiority, availability, and effectiveness of the algorithm.

1. Introduction

The Earth’s gravitational field is one of the essential fundamental physical fields of the Earth [1,2]. Gravitational field data are national strategic data that are irreplaceable in basic mapping, national defense and military, resource exploration, earth science, aerospace, and national integrated PNT system construction [3]. Accurate gravity measurements are obtained with the help of gravimeters, which can be divided into platform type [4,5] and strapdown type [6,7]. The platform-type gravimeters primarily consist of a gravity sensor and stabilization platform, where the gravity sensor is responsible for the accurate measurement of gravity value information and ensuring that it has a stable vertical pointing during its movement. This requires the use of a gravity measurement stabilization platform, and the installation relationship between the two is shown in Figure 1.
As shown in Figure 1, the horizontal attitude error of stable gravity measurements is θ , and the vertical deviation of the sensitive axis of the gravity sensor is also θ . In the figure, g represents the actual value of gravity, g represents the gravity value sensed by the gravity sensor, a h represents the horizontal acceleration of the carrier, and a h represents the influence of the horizontal acceleration of the page on the measured value of the gravity sensor.
Then, the gravity measurement error can be expressed as:
e g = g a h g = g cos θ a h sin θ g a h θ g θ 2 / 2
In Equation (1), the second term is a steady-state error, which exists even if the carrier does not move.
If the gravity static measurement accuracy of 1 mGal is to be achieved, the horizontal misalignment angle must be less than 4.8′ when the first item in Equation (1) is ignored. The magnitude of the first term in Equation (1) depends on the horizontal acceleration of the carrier and the extent of the error angle of the platform. If the gravity sensor’s damping mechanism is not considered, the maximum dynamic gravity measurement error should not exceed 1 mGal. When the horizontal acceleration of the carrier is a h = 0.1   m / s 2 , the error angle θ should be less than 20.6″. In actual measurement systems, the damping mechanism, such as silicone oil, is set in the gravity sensor, which gives the gravity measurement system a specific anti-interference ability against the horizontal acceleration of the carrier. In practical applications, the error angle θ does not need to be less than 20.6″. However, Equation (1) indicates that, in the case of the horizontal misalignment angle of the gravimetric measuring platform, the horizontal acceleration of the carrier has a non-negligible effect on the gravimetric measurement results.
As large-weight instrumentation, a gravimeter requires a gravimetric stabilization platform with a hefty load capacity (more than 30 kg). The related research mainly studies the control algorithm of a small load, three-axis or airborne stabilization platform. There is less research on the control of high-precision stabilization platforms with enormous load capacities. Domestic and foreign research scholars have conducted relevant studies to reduce the stabilized platform’s control system error and improve the stabilized platform’s stabilization accuracy. In Ref. [8], a robust adaptive controller in the form of expectation compensation was proposed for accurate tracking control, to reduce the influence of measurement noise and to improve the control performance of the photoelectric gyro-stabilized platform system. In Ref. [9], a new decoupling algorithm, applied to the stabilized platform for aerial remote sensing, was used to solve the problems of unknown nonlinear interference and coupling interference, which occur in tasks such as aerial photography. In Ref. [10], a combination of sliding mode control and a neural network is proposed to deal with the uncertainty disturbances in the inertially stabilized platform (ISP) system model. The article also presents an adaptive neural network to approximate the uncertainty and unknown disorders of the system and improve control performance. An anti-disturbance control scheme for an airborne radar stabilization platform based on the inverse estimation algorithm of self-rejecting control (ADRC) is used to solve the stability and clarity problems of radar imaging [11]. An improved dynamic variational differential evolution algorithm (DMDE) is proposed to optimize the PID control parameters and is validated by simulation for several standard industrial-controlled object models [12]. An adaptive fuzzy PID composite control algorithm is proposed and applied to control gyro-stabilized platforms for resisting nonlinear disturbance factors in platform control [13]. A fuzzy PID composite control is proposed to improve stability and resistance to disturbances in aero-inertial stabilized platforms [14]. A fuzzy PID intelligent control algorithm is proposed to solve the low control efficiency of small stable platforms due to frame coupling [15]. Most of the control algorithms proposed in these studies have only been tested in simulations. They are rarely applied to existing stabilized platforms, so studying an intelligent control algorithm that can be used in actual gravimetric stabilized platforms is crucial.

2. System Description

2.1. System Architecture

The gravimetric stabilization platform in this paper was developed in the Marine/Aerospace Gravimeter project, which is a primary national science and technology project with a two-axis structure (pitch axis and roll axis). The mechanical structure is shown in Figure 2.
The relevant coordinate system is defined as follows [16].
  • Inertial coordinate system i : The inertial coordinate system has the center of the Earth as the origin, and the x i and y i axes are in the equatorial plane of the Earth, where the x i and z i axes are set to point from the head to the equinox and along the Earth’s rotation axis to the North Pole direction, respectively.
  • Earth coordinate system e : The Earth coordinate system has the center of the Earth as the origin, and the x e and y e axes are in the Earth’s equatorial plane, where the x e and z e axes are set to point from the head to the prime meridian and along the Earth’s rotation axis to the North Pole direction, respectively.
  • Geographic coordinate system n : The geographic coordinate system takes the center of mass of the carrier as the origin; x n , y n and z n axes are set to point east, north, and skyward from the head, respectively.
  • Carrier coordinate system b : The carrier coordinate system takes the center of mass of the carrier as the origin, and the x b , y b , and z b axes are set to point to the right side of the carrier, to the front of the carrier, and the top of the page, respectively, and form a right-handed coordinate system.
  • Outer frame (roll frame) coordinate system r : set the direction of y r axes to be the same as the direction of the y b axes, and its coordinate system can only rotate around the y b axes during the subsequent movements relative to the coordinate system b . The resulting roll angle is set to γ r .
  • The inner frame (pitch frame) coordinate system f : set the direction of the x f axes to be the same as the pointing of the x b axes, and its coordinate system can only rotate around the x r axes during the subsequent motions relative to the coordinate system r . The resulting pitch angle is set to θ f .

2.2. Principle of Isolation Carrier Movement

The structure of the gravimetric stabilization platform was designed so that the IMU is fixed directly under the gravimeter base. When the carrier moves or is otherwise disturbed, the inner and outer frames will shake slightly, resulting in the stabilization platform not coinciding with the local geographic level. Let ω n b b = [ ω n b x b   ω n b y b   ω n b z b ] T be the projection of the angular velocity of the carrier coordinate system relative to the geographic coordinate system on the carrier coordinate system; the pitch axis compensation angular velocity generated by the torque motor in the stabilization loop is θ ˙ f , and the roll axis compensation angular velocity is γ ˙ r . Since the angular velocity output from the IMU is the angular velocity of the gravimeter relative to the inertial space, the angular velocity ω i b b cannot be obtained by direct measurement from ω n b b = ω i b b ω i n b = ω i b b C n b ω i n n , where C n b denotes the attitude array of the coordinate navigation system b relative to the carrier system n , ω i b b is the gyro output, and ω i n n indicates the rotation of the n system close to the i system. This contains two components: the wheel of the navigation system caused by the course of the Earth, and the process of the inertial guidance system moving near the surface of the Earth due to the n-system process caused by the bending of the Earth’s surface. Accordingly, ω i n n = ω i e n + ω e n n .
where ω i e n = [ 0   ω i e cos L   ω i e sin L ] T ,
ω e n n = [ v N R M + h v E R N + h v E R N + h tan L ] T
where ω i e is the angular rate of Earth’s rotation; L and h are the geographic dimension and altitude, respectively; v N and v E are the northward and eastward velocities, respectively; R M and R N are the radius of principal curvature of the meridian circle and the radius of principal curvature of the Uranus circle, respectively.
The angular velocity of the transverse roll axis for the geographical coordinate system is obtained based on the relevant coordinate transformations.
ω n r r = C b r ω n b b + ω b r r = [ cos γ r 0 sin γ r 0 1 0 sin γ r 0 cos γ r ] [ ω n b x b ω n b y b ω n b z b ] + [ 0 γ ˙ r 0 ] = [ ω n b x b cos γ r + ω n b z b sin γ r ω n b y b + γ ˙ r ω n b x b sin γ r + ω n b z b cos γ r ]
where ω n r r is the projection of the angular velocity of the rotation of the r system to the n system on the r system; C b r is the coordinate system transformation matrix from the r system to the b system; ω n b b is the projection of the angular velocity of the rotation of the b system on the n system on the b system; ω b r r is the projection of the angular velocity of the rotation of the r system to the rotation of the b system on the r system; γ r is the roll angle; γ ˙ r is the rolling angular velocity; ω n b x b , ω n b y b , and ω n b z b are the components of the projection of the angular velocity of the rotation of the b system to the n system on the b system on the x , y and z axes, respectively.
Similarly, the angular velocity of the pitch axis for the geographical coordinate system can be obtained.
ω n f f = C r f ω n r r + ω r f f = [ 1 0 0 0 cos θ f sin θ f 0 sin θ f cos θ f ] [ ω n r x r ω n r y r ω n r z r ] + [ θ ˙ f 0 0 ] = [ ω n b x b cos γ r + ω n b z b sin γ r + θ ˙ f ( ω n b y b + γ ˙ r ) cos θ f + A sin θ f ( ω n b y b + γ ˙ r ) sin θ f A cos θ f ]
where A = ω n b x b sin γ r ω n b z b cos γ r ω n f f is the projection of the angular velocity of the rotation of the f system to the n system on the f system; C r f is the coordinate system transformation matrix from the f system to the r system; ω n r r is the projection of the angular velocity of the rotation of the r system on the n system on the r system; ω r f f is the projection of the angular velocity of the rotation of the f system to the rotation of the r system on the f system; θ f is the pitch angle; θ ˙ f is pitch angular velocity; ω n r x r , ω n r y r , and ω n r z r are the components of the projection of the angular velocity of the rotation of the r system to the n system on the r system on the x , y and z axes, respectively.
To maintain the gravimetric stabilization platform plane at all times within the local geographic level, the angular velocity of the inner frame relative to the geographic coordinate system must be 0 after the initial righting, and the following conditions must be met.
{ ω n f x f = 0 ω n f y f = 0
Equations (3) and (4) can be obtained by the two torque motors to produce the angular speed, meeting the following conditions.
{ θ ˙ f = ω n b z b sin γ r ω n b x b cos γ r γ ˙ r = ( ω n b x b sin γ r + ω n b z b cos γ r ) tan θ f ω n b y b
When the control result satisfies Equation (5), the gravimetric stabilization platform table surface will remain parallel to the local geographical level.

2.3. Control System Block Diagram

The servo control system is the core of the gravimetric stabilization platform. Additionally, it ensures algorithm implementation and platform stability accuracy, so a scientific and reasonable control system structure is essential. The gravity measurement stable platform’s hardware structure consists of a pitch axis and a roll axis with the same design, without considering the influence of coupling factors between the two axes. Only one of the axes can be analyzed when analyzing its control law, and the pitch axis is used as an example in the text with no further explanation. Figure 3 depicts the system composition and the principle block diagram.
Given that the stabilized platform must install heavy equipment such as an IMU and a gravimeter and has stringent real-time requirements, its control system employs a multi-stage composite control scheme, which utilizes dual position and velocity rings to achieve composite control. Figure 4 depicts the hybrid control structure diagram.
Where θ d is the pitch axis target angle, θ ˙ is the angular velocity, θ is the actual output pitch axis angle, the load includes the gravimeter and the IMU fixedly connected.

3. Improved Differential Evolutionary Adaptive Fuzzy PID Control Algorithm

The gravimetric stabilized stage is a servo-controlled system with strong nonlinearity. Due to its uncertainties, such as random disturbances, friction, mechanical coupling, and gyroscopic drift, the exact mathematical model of the stabilized platform cannot be determined. In contrast, fuzzy control is a kind of control based on fuzzy rules, which is used to realize the control of a nonlinear system without a precise mathematical model [17]. Figure 5 depicts the structure of a typical fuzzy control system.
Fuzzy control can adjust the torque motor’s control voltage according to the fuzzy rules. However, its steady-state accuracy generally does not meet the accuracy requirements of the gravimetric stabilization platform. Therefore, the control of the stabilized platform is usually realized by the PID controller, which can achieve high control and stabilization accuracy with good real-time accuracy and strong robustness when it has suitable control parameters. To meet the requirements of high precision and robust system stability, a combination of fuzzy control and PID control can realize the real-time adjustment of the control parameters using the control system according to the error and the differential of the error. Then, the adaptive fuzzy PID control of the system can be determined.
The accuracy of the initial parameter-setting of adaptive fuzzy PID control significantly impacts the control performance. In the conventional adaptive fuzzy PID control, the selection of initial control parameters is often determined based on experience or traditional parameter tuning methods, and its accuracy is inadequate. The differential evolution algorithm is an intelligent algorithm that can perform automatic optimization searches. An adaptive differential evolution algorithm based on a multivariate strategy with variable parameters (ADE-MMVP) is attempted to optimize the initial control parameters of an adaptive fuzzy PID control and realize the accurate control of a stable platform.

3.1. ADE-MMVP Algorithm

Differential Evolution (DE) was proposed by Professor Rainer Storn and Professor Kenneth Price in 1995 [18,19]. DE is mainly implemented through competition and selection among individuals in the population to optimize the population. It is a global search algorithm. The differential evolution algorithm adopts an accurate number coding and is based on the differential strategy mutation operation, as in Ref. [20], which can effectively reduce genetic complexity. At the same time, the algorithm also has a specific memory ability, so the global search ability of the algorithm is greatly improved, and the robustness is enhanced. In general, the differential evolution algorithm has the advantages of fewer undetermined parameters, a fast convergence speed, and the ability to easily obtain the global optimum. Compared with other optimization algorithms, the differential evolution algorithm has the following characteristics: [21].
(1)
The ability to deal with nonlinear, non-differentiable, and multimodal functions;
(2)
The ability to process intensive cost functions in parallel;
(3)
It is easy to use and can realize algorithm optimization while making good use of control variables;
(4)
Good convergence, which can converge to the optimal global value in continuous independent tests.
Due to the single variation strategy, fixed variation factor, and crossover probability factor of the conventional differential evolution algorithm, the algorithm is prone to slow convergence and failure to converge to the global optimum. To solve this problem, we propose an ADE-MMVP algorithm with multiple variation strategies, various factors, and crossover probability factors that can adaptively change with population evolution [22].
  • Variation strategy of ADE-MMVP algorithm
Among all DE variation strategies, DE/rand/1 is the most widely used and most beneficial for maintaining population diversity, as in Refs. [23,24] DE/best/2 has the optimal solution, which beneficial to resolving some technical problems with the algorithm and accelerate the convergence speed, but the variation strategy containing optimal information is more likely to fall into local optimum [25]. Based on this, the two variational strategies are combined using a certain ratio and adjusted by a scaling factor, and then adapted in a variety of application scenarios. The specific calculation process is as follows.
h i 1 ( t + 1 ) = x r 1 ( t ) + F ( x r 2 ( t ) x r 3 ( t ) )
h i 2 ( t + 1 ) = x best ( t ) + F ( x r 1 ( t ) x r 2 ( t ) ) + F ( x r 4 ( t ) x r 5 ( t ) )
h i ( t + 1 ) = λ × h i 1 ( t + 1 ) + ( 1 λ ) × h i 2 ( t + 1 )
where t is the current number of iterations, x r 1 ( t ) ~ x r 5 ( t ) are the randomly selected mutually dissimilar individuals in the initial population, x i ( t ) is the current individual, x best ( t ) is the current optimal individual, F is the variation factor, λ is the proportion factor of DE/rand/1 in the variation strategy.
2.
Variation factor of ADE-MMVP algorithm
The variation factor F mainly controls the search step size of the differential evolution algorithm, which affects the algorithm population’s diversity and convergence. The algorithm uses a variation factor that dynamically adjusts with the number of iterations.
x ( G ) = e 1 G m G m + 1 G
F ( G ) = F min + x ( G ) ( F max F min )
where G denotes the current number of iterations and the maximum number of iterations, G m denotes the value of the variable factor of the current number of iterations, and F max and F min represent the maximum and minimum values of the variation factor, respectively.
3.
Crossover probability factor of ADE-MMVP algorithm
The crossover probability factor C R controls the degree of participation of randomly selected individuals in the algorithm crossover, thus balancing the role of the local search and global search relationship. The following equation adjusts the algorithm crossover probability factor C R .
C R = C R min + ( C R max C R min ) G max G
where G max denotes the maximum number of iterations, C R min denotes the minimum value of the set crossover probability factor, C R max denotes the maximum value of the set crossover probability factor, and G denotes the current number of iterations.

3.2. Adaptive Fuzzy PID Control Algorithm with ADE-MMVP Optimization

To solve the problem regarding the accuracy of the initial control parameters of the adaptive fuzzy PID controller, [26] ADE-MMVP is used for automatic optimization of the control parameters, and the optimization results are transmitted to the adaptive fuzzy PID controller. The system uses a dual-loop velocity loop and position loop to achieve composite control.
In the platform mathematical model, the highest numerator and denominator of the speed loop transfer function are the same, and PI control is used in the speed loop control system to ensure the system’s stable operation. The structure of the speed loop control system is shown in Figure 6.
Structure of the position loop control system is shown in Figure 7.
Set several rules for PID controller control parameter adjustment according to the control system requirements.
  • When the system error | e | is large, a more significant proportional link coefficient should be used to speed up the response of the system. However, the proportional coefficient cannot be infinite or the system will have a massive amount of overshoot and the adjustment time of the system will increase, as the system error | e | is significant at the beginning of the system. To avoid the system control exceeding the maximum execution range of the actuator, a smaller differential coefficient should be taken at this time to speed up the system response. To avoid causing a massive overshoot to the system, the integral link should be removed when the error | e | is large and the integral coefficient K i = 0 is taken.
  • When the system error | e | is moderately large, the proportional link coefficient should be appropriately reduced to take a minor K p . This will prevent the system from having a massive amount of overshoot, resulting in system collapse. The response speed of the system and the value of the differential link are directly related, so the value of differential link coefficient K d at this time is critical. The magnitude of the integral link coefficient can be increased in this scenario, but this increase should not be too tremendous.
  • When the system error | e | is minor, to ensure the system has a good steady-state performance, this time can use a more significant differential link coefficient and integral link coefficient. To avoid system oscillation, the value of K d should be obtained appropriately.

3.3. Simulation and Analysis

Table 1 displays the ADE-MMVP algorithm parameter settings used in the simulation trials. The rectification outcomes of the ZN critical proportionality technique define the initial parameter values for the adaptive fuzzy control algorithm (AFC). In contrast, the improved differential evolution adaptive fuzzy PID control’s (IDEAFC) initial control parameters are established by following the best control parameters, as identified by the ADE-MMVP method.
The response of the speed loop under the two signals is given in Figure 8, which depicts the simulation of the speed and position loops using the sinusoidal and sawtooth wave signals, respectively. Figure 9 displays the position loop’s reaction to the two movements.
The experimental results show that the IDEAFC algorithm is superior.

4. Platform Experiments

Based on the simulation test results, this chapter presents the platform experiments that were conducted to further verify the algorithm’s superiority. The stable platform control system had a TMS320F28335 digital processor as the control core. The traditional PID control algorithm, AFC algorithm, and IDEAFC algorithm were used in laboratory static and dynamic rocking experiments (longitudinal and transverse rocking experiments, respectively). The electronic self-collimator and IMU data were collected to control each algorithm. The effect of each algorithm was compared by collecting data from the electronic self-collimator and IMU. The control effects were compared by collecting data from the electronic self-collimator and the IMU. The IDEAFC algorithm was also used for onboard experiments and shipboard experiments, and the control effects were analyzed using the data collected from the IMU.

4.1. Shaking Table Static Stability Experiment

The gravimetric stabilization platform was subjected to a shaking table static stabilization experiment in a laboratory environment. The angular offsets in the investigation were collected from IMU attitude data and TriAngle electronic self-collimator data, respectively. The physical objects of IMU and the electronic self-collimator are shown in Figure 10; the experimental environment of the platform laboratory is shown in Figure 11.
In the static experiment, the angle measured by the electronic self-collimator is the angle deviation value of the overall stabilized platform, which can truly reflect the overall accuracy of the forum. In contrast, the angle deviation value measured by the IMU can reflect the control accuracy of the control algorithm. GNSS positioning-assisted calibration was used throughout the experiment to avoid the influence of measurement errors caused by the IMU’s own drift in the test results.
When running the tests, multiple experiments were conducted, utilizing the conventional PID method, the AFC algorithm, and the IDEAFC algorithm. The set closest to the mean value was chosen as the outcome. Figure 12 depicts the IMU pitch and roll angle data gathered throughout the studies.
As seen from the comparison in Figure 12, in the experimental results obtained by IMU under static conditions, whether pitch angle or roll angle, the control effect of the IDEAFC algorithm was shown to be better than the other two. The standard deviation of the pitch angle of the IDEAFC algorithm was 0.033″, with a peak value of 0.002°; the standard deviation of the roll angle was 1.5″, with a peak value of 0.002°.
The data collected by the electronic self-collimator and IMU were compared separately to study the relationship between the overall stability accuracy and control accuracy of the platform. Since the electronic self-collimator is a high-precision optical goniometer, its reflector is small and fixed, so it can only be used under static test conditions. The relationship between the platform’s overall stability and control accuracy were evaluated from the comparison experiment under static conditions, and then used to analyze whether the stability accuracy of various algorithms meets the requirements.
In the static experiment, the pitch angle and the roll angle measured by the electronic self-collimator are shown in Figure 13.
A comparison of the control effects of the three control algorithms can be obtained from the comparison chart, as shown in Table 2.
The improvement in accuracy was calculated as follows.
Relative comparison algorithm precision improvement ratio = (comparison algorithm precision − target algorithm precision)/comparison algorithm precision.

4.2. Dynamic Rocking Experiment

The dynamic swaying experiment of the gravity measurement stable platform was conducted in the laboratory using a swaying table with adjustable longitudinal and transverse swaying amplitudes. The maximum value of longitudinal swaying was ±45°, and the total value of transverse swaying was ±30°. In the natural gravity measurement environment, whether vehicle-mounted, on the lake, at sea, or under the sea, the swaying angle generally did not exceed ±10°. To test the stability of the platform, when conducting the platform swaying experiment, the angle of both longitudinal and transverse swaying was set to ±15°. At the same time, to avoid the impact caused by the different swing speeds, the multi-grade swing experiment was conducted separately. The swing speed of each gear was different, and the period setting for each pack is shown in Table 3.
Figure 14 compares pitch angle data and roll angle data at different gears during vertical rocking.
The comparison chart shows that the control effect of the algorithm is the same under different shaking table operation gears, so the impact of other gears on the control algorithm were not considered in the subsequent comparison test.
A comparison of the data measured by IMU under the three control algorithms during the longitudinal rocking test is shown in Figure 15.
The comparison graphs show the control effects of the three control algorithms in the longitudinal rocking experiment, as shown in Table 4.
The comparison of data under the three control algorithms was measured by IMU during the swing test, as shown in Figure 16.
Using the comparison graphs, we can compare the control effects of the three control algorithms in the swing experiment, as shown in Table 5.

4.3. In-Vehicle Experiments

The gravimetric stabilization platform was tested in-vehicle to verify the stability performance of the IDEAFC algorithm in real application scenarios. The designed route in the experiment was Wuhan city road, the vehicle driving speed was 30 km/h, and an onboard experimental environment was built, as shown in Figure 17.
The vehicle-based experimental driving city route was as shown in Figure 18.
The pitch and roll angle data were collected during the experiment, as shown in Figure 19.
During the experiment, the standard deviation of the pitch angle of the gravimetric stabilized platform was 7.56″, and the standard deviation of the roll angle was 7.09″, both within 10″.

4.4. Shipboard Experiments

A gravimetric stabilization platform using an IDEAFC algorithm was used for a boat experiment conducted at Mulan Lake, Wuhan, using a small electric pleasure boat. The boat is shown in Figure 20 while docked at the experimental dock. GNSS was used to calibrate the inertial guidance during the experiment, and the GNSS antenna was installed directly above the IMU on the boat top platform, as shown in Figure 21.
The experimental waters contained light waves, with a wind speed of about 1 m/s; the speed of the boat in the lake was 15 km/h. From the docked experimental dock, sailing experiments were conducted around Mulan Lake; the whole route map is shown in Figure 22.
The pitch and roll angle data obtained throughout the cruise are shown in Figure 23.
During the whole shipboard experiment, the standard deviation of the gravity measurement stabilized platform pitch angle was 2.11″, and the standard deviation of the roll angle was 3.69″, both within 10″.

5. Conclusions

An improved differential evolutionary adaptive fuzzy PID control method was proposed to reduce the influence of multi-source perturbations on the gravimetric stabilization platform, achieve high precision control, and improve stabilization accuracy. As a comparison, conventional PID and adaptive fuzzy control were also applied to the platform control. Simulation experiments and platform tests were carried out for different algorithms. The experimental results show that the improved differential evolution algorithm can significantly improve the stabilization accuracy of the gravimetric stabilization platform. In the laboratory swing experiments, the control accuracy was enhanced by more than 25% under static conditions and more than 18% under dynamic conditions. In the onboard and shipboard experiments, the improved differential evolution adaptive fuzzy PID control algorithm pitch and roll angle standard deviation were within 10″. This lays a foundation for subsequent gravity measurements under static or dynamic conditions and creates the necessary conditions for the realization of high-precision dynamic gravity measurements.

Author Contributions

Conceptualization, H.B.; Methodology, H.H.; Validation, F.L.; Writing—review & editing, X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 41804076).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors thank Ning Mao, Ding Li and En-Fan Lin for their help in the experiment.

Conflicts of Interest

These authors declare no conflict of interest.

References

  1. Bidel, Y.; Zahzam, N.; Bresson, A.; Blanchard, C.; Cadoret, M.; Olesen, A.V.; Forsberg, R. Absolute airborne gravimetry with a cold atom sensor. J. Geod. 2020, 94, 1–9. [Google Scholar] [CrossRef] [Green Version]
  2. Kiani, M. Image Gravimetry: A new remote sensing approach for gravity analysis in geophysics. arXiv 2003, arXiv:2003.09388. [Google Scholar]
  3. Bian, H.; Xu, J.; He, H.; Wang, R.; Ma, H. The concept of resilience of national comprehensive PNT system. Geomat. Inf. Sci. Wuhan Univ. 2021, 46, 1265–1272. [Google Scholar]
  4. Yuan, Y.; Gao, J.; Wu, Z.; Shen, Z.; Wu, G. Performance estimate of some prototypes of inertial platform and strapdown marine gravimeters. Earth Planets Space 2020, 72, 89. [Google Scholar] [CrossRef]
  5. Glennie, C.L.; Schwarz, K.P.; Bruton, A.M.; Forsberg, R.; Olesen, A.V.; Keller, K. A comparison of stable platform and strapdown airborne gravity. J. Geod. 2000, 74, 383–389. [Google Scholar] [CrossRef]
  6. Wang, W.; Gao, J.; Li, D.; Zhang, T.; Luo, X.; Wang, J. Measurements and Accuracy Evaluation of a Strapdown Marine Gravimeter Based on Inertial Navigation. Sensors 2018, 18, 3902. [Google Scholar] [CrossRef] [Green Version]
  7. Jensen, T.E.; Forsberg, R. Helicopter Test of a Strapdown Airborne Gravimetry System. Sensors 2018, 18, 3121. [Google Scholar] [CrossRef] [Green Version]
  8. Wu, Y.; Yue, D. Desired compensation adaptive robust control of electrical-optical gyro-stabilized platform with continuous friction compensation using modified Lugre model. Int. J. Control. Autom. Syst. 2018, 16, 2264–2272. [Google Scholar] [CrossRef]
  9. Wang, Z. Research on Control Algorithm of Aerial Remote Sensing Stabilized Platform. Master’s Thesis, Changchun University of Technology, Jilin, China, 2021. [Google Scholar]
  10. Zou, Y.; Lei, X. A compound control method based on the adaptive neural network and sliding mode control for inertial stable platform. Neurocomputing 2015, 155, 286–294. [Google Scholar] [CrossRef]
  11. Mei, D.; Yu, Q. Disturbance rejection control of airborne radar stabilized platform based on active disturbance rejection control inverse estimation algorithm. Assem. Autom. 2021, 41, 525–535. [Google Scholar] [CrossRef]
  12. Tan, F.; Cao, L. Optimal PID control based on the improved dynamic mutation differential evolution algorithm. Control. Eng. China 2019, 26, 461–468. [Google Scholar]
  13. Zhu, Y.; Lu, Y.; Xu, J.; Cheng, X.-H. Adaptive fuzzy-PID hybrid control method of gyro stabilized platform. J. Chin. Inert. Technol. 2014, 22, 317–321. [Google Scholar]
  14. Zhou, X.; Jia, Y. Fuzzy/PID compound control for inertially stabilized platform in airborne remote sensing. Chin. J. Sci. Instrum. 2016, 37, 2545–2554. [Google Scholar]
  15. Yan, D. A method to solve the problem of low precision of micro stabilized platform caused by frame coupling-Based on Fuzzy PID. E3S Web Conf. 2021, 284, 04008. [Google Scholar] [CrossRef]
  16. Yan, G.; Wong, J. Strapdown Inertial Navigation Algorithm and Integrated Navigation Principle; Northwestern Polytechnical University Press: Xi’an, China, 2019; Chapter 4; pp. 78–79. [Google Scholar]
  17. Yin, H.; Yi, W.; Wang, K.; Guan, J.; Wu, J. Research on brushless DC motor control system based on fuzzy parameter adaptive PI algorithm. AIP Adv. 2020, 10, 105208. [Google Scholar] [CrossRef]
  18. Storn, R.; Price, K. Minimizing the real functions of the ICEC’96 contest by differential evolution. In Proceedings of the IEEE International Conference on Evolutionary Computation, Nagoya, Japan, 20–22 May 1996; pp. 842–844. [Google Scholar]
  19. Zhang, J.; Duan, J. Robot Dynamic Parameter Identification Based on Improved Differential Evolution Algorithm. J. Beijing Union Univ. 2020, 34, 48–56. [Google Scholar]
  20. Bai, X. Research on Wind Power Forecasting and Unit Dispatch of AGC; Beijing Jiaotong University: Beijing, China, 2009. [Google Scholar]
  21. Storn, R.; Price, K. Differential Evolution—A Simple and Efficient Heuristic for global Optimization over Continuous Spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  22. Chen, X.; Li, F.; He, H.; Wu, M. Optimization of PID parameter tuning for gravity stabilized platform based on improved differential evolutionary algorithm. J. Phys. Conf. Series. IOP Publ. 2021, 2029, 012107. [Google Scholar] [CrossRef]
  23. Price, K. Differential evolution. In Handbook of Optimization; Springer: Berlin, Heidelberg, 2013; pp. 187–214. [Google Scholar]
  24. Das, S.; Suganthan, P. Differential evolution: A survey of the state-of-the-art. IEEE Trans. Evol. Comput. 2010, 15, 4–31. [Google Scholar] [CrossRef]
  25. Zhang, J.; Sanderson, A. JADE: Adaptive differential evolution with optional external archive. IEEE Trans. Evol. Comput. 2009, 13, 945–958. [Google Scholar] [CrossRef]
  26. Zhou, X.; Shu, T.; Lyu, Z.; Sun, B. Sliding mode control of inertially stabilized platform based on fuzzy switching gain adjustment. Chin. J. Sci. Instrum. 2021, 42, 263–271. [Google Scholar]
Figure 1. Gravity measurement of stable platform and gravity sensor installation relationship.
Figure 1. Gravity measurement of stable platform and gravity sensor installation relationship.
Sensors 23 03172 g001
Figure 2. Stabilization platform mechanical structure.
Figure 2. Stabilization platform mechanical structure.
Sensors 23 03172 g002
Figure 3. Stabilized platform system composition and the principle block diagram.
Figure 3. Stabilized platform system composition and the principle block diagram.
Sensors 23 03172 g003
Figure 4. The hybrid control structure diagram.
Figure 4. The hybrid control structure diagram.
Sensors 23 03172 g004
Figure 5. Typical fuzzy control system structure.
Figure 5. Typical fuzzy control system structure.
Sensors 23 03172 g005
Figure 6. Structure of speed loop control system.
Figure 6. Structure of speed loop control system.
Sensors 23 03172 g006
Figure 7. Structure of the position loop control system.
Figure 7. Structure of the position loop control system.
Sensors 23 03172 g007
Figure 8. (a) Comparison curve of sine signal response of speed ring; (b) comparison curve of sawtooth wave response of speed ring gravimeter.
Figure 8. (a) Comparison curve of sine signal response of speed ring; (b) comparison curve of sawtooth wave response of speed ring gravimeter.
Sensors 23 03172 g008
Figure 9. (a) Comparison curve of sine signal response of location ring; (b) comparison curve of sawtooth wave response of location ring.
Figure 9. (a) Comparison curve of sine signal response of location ring; (b) comparison curve of sawtooth wave response of location ring.
Sensors 23 03172 g009
Figure 10. (a) IMU; (b) electronic autocollimator.
Figure 10. (a) IMU; (b) electronic autocollimator.
Sensors 23 03172 g010
Figure 11. The experimental environment of the platform laboratory.
Figure 11. The experimental environment of the platform laboratory.
Sensors 23 03172 g011
Figure 12. (a) Comparison of pitching angle data; (b) comparison of rolling angle data.
Figure 12. (a) Comparison of pitching angle data; (b) comparison of rolling angle data.
Sensors 23 03172 g012
Figure 13. (a) Comparison of static pitch angle data; (b) comparison of static roll angle data.
Figure 13. (a) Comparison of static pitch angle data; (b) comparison of static roll angle data.
Sensors 23 03172 g013
Figure 14. (a) Dynamic experimental pitch angle data; (b) dynamic experimental roll angle data.
Figure 14. (a) Dynamic experimental pitch angle data; (b) dynamic experimental roll angle data.
Sensors 23 03172 g014
Figure 15. (a) Pitch angle data comparison; (b) roll angle data comparison.
Figure 15. (a) Pitch angle data comparison; (b) roll angle data comparison.
Sensors 23 03172 g015
Figure 16. (a) Pitch angle data comparison; (b) roll angle data comparison.
Figure 16. (a) Pitch angle data comparison; (b) roll angle data comparison.
Sensors 23 03172 g016
Figure 17. Vehicle experiment environment.
Figure 17. Vehicle experiment environment.
Sensors 23 03172 g017
Figure 18. Vehicle experimentation route.
Figure 18. Vehicle experimentation route.
Sensors 23 03172 g018
Figure 19. Vehicle experimentation results.
Figure 19. Vehicle experimentation results.
Sensors 23 03172 g019
Figure 20. Experimental vessel.
Figure 20. Experimental vessel.
Sensors 23 03172 g020
Figure 21. Installation position of GNSS antenna.
Figure 21. Installation position of GNSS antenna.
Sensors 23 03172 g021
Figure 22. Navigation trajectory.
Figure 22. Navigation trajectory.
Sensors 23 03172 g022
Figure 23. Ship-borne test results.
Figure 23. Ship-borne test results.
Sensors 23 03172 g023
Table 1. Parameter setting.
Table 1. Parameter setting.
Parameter NameADE-MMVP
Number of iterations500
Population size20
Individual Dimension3
Code length20
Scale factor λ 0.35
Variable factorsMax 1.8 Min 0.4
Crossover probability factorMax 0.8 Min 0.2
Table 2. Comparison of control effects under static conditions.
Table 2. Comparison of control effects under static conditions.
Pitch Angle StdRoll Angle Std
Traditional PID control overall stability accuracy2.98″3.36″
AFC overall stability accuracy2.56″0.96″
IDEAFC overall stability accuracy1.35″0.72″
Improved relative to traditional PID control54.7%78.6%
Relative AFC improves47.3%25%
Table 3. Rolling period.
Table 3. Rolling period.
PitchRoll
Gear 140 s40 s
Gear 220 s20 s
Gear 315 s15 s
Table 4. Comparison of control effect under pitching conditions.
Table 4. Comparison of control effect under pitching conditions.
Pitch Angle StdRoll Angle Std
Traditional PID control overall stability accuracy16.00″17.27″
AFC overall stability accuracy8.92″12.23″
IDEAFC overall stability accuracy7.28″8.90″
Improved relative to traditional PID control54.5%48.5%
Relative AFC improves18.4%27.2%
Table 5. Comparison of control effects under dynamic rolling conditions.
Table 5. Comparison of control effects under dynamic rolling conditions.
Pitch Angle StdRoll Angle Std
Traditional PID control overall stability accuracy15.61″16.60″
AFC overall stability accuracy11.92″13.37″
IDEAFC overall stability accuracy6.72″9.82″
Improved relative to traditional PID control60.0%40.8%
Relative AFC improves43.6%26.6%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, X.; Bian, H.; He, H.; Li, F. An Improved Differential Evolution Adaptive Fuzzy PID Control Method for Gravity Measurement Stable Platform. Sensors 2023, 23, 3172. https://doi.org/10.3390/s23063172

AMA Style

Chen X, Bian H, He H, Li F. An Improved Differential Evolution Adaptive Fuzzy PID Control Method for Gravity Measurement Stable Platform. Sensors. 2023; 23(6):3172. https://doi.org/10.3390/s23063172

Chicago/Turabian Style

Chen, Xin, Hongwei Bian, Hongyang He, and Fangneng Li. 2023. "An Improved Differential Evolution Adaptive Fuzzy PID Control Method for Gravity Measurement Stable Platform" Sensors 23, no. 6: 3172. https://doi.org/10.3390/s23063172

APA Style

Chen, X., Bian, H., He, H., & Li, F. (2023). An Improved Differential Evolution Adaptive Fuzzy PID Control Method for Gravity Measurement Stable Platform. Sensors, 23(6), 3172. https://doi.org/10.3390/s23063172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop