Ground Simulation Tests in Two-Dimensional Dynamic Acceleration Environment

: Acceleration environment (AE) simulation tests have been widely applied in various areas such as aviation, space, environmental medicine, medical science, biomedicine, and materials. Most existing AE simulation test methods adopt impact and ﬂight overload simulation separately, which cannot realize the synchronous controllable loading of two-dimensional (2-D) dynamic acceleration. In this paper, we propose and implement an AE simulation test method using an air-cannon-started compound centrifugal turntable. Speciﬁcally, our method emphasizes two processes, i.e., orderly transition from impact overload to steady-state overload and synchronous loading of 2-D dynamic AE. To further facilitate the method veriﬁcation in AE, we develop a specialized simulation test platform. During ﬁeld tests, axial acceleration of the unit being tested reached 20 g in 0.12 s at the launching stage and radial acceleration raised from zero g to 40 g in 0.5 s at the ﬂight stage, achieving simulation of the 2-D dynamic AE. The obtained results show the e ﬀ ectiveness of the proposed simulation method, shedding light on updated design and control of impact simulation tests on the structures of crystal oscillator circuit, acceleration sensors, and related instruments.


Introduction
Acceleration environment (AE) can generate accelerated environmental forces by means of simulation tests, which is used for performance testing and reliability experiments. AE has been widely applied in various areas such as aviation, space, environmental medicine [1], medical science [2], biomedicine [3], materials [4], and acceleration sensors [5]. The acceleration process of aircraft includes two stages: launch and maneuvering flight. The AE of aircraft is usually evaluated by the axial acceleration and the radial acceleration. The greater the axial acceleration is, the faster the speed change of the aircraft will be. The greater the radial acceleration is, the better the maneuverability of the aircraft will be. To perform ground-based tests of the acceleration sustaining ability of aircraft's inertia device and aerospace materials and to verify the strength and reliability of the system, it is necessary to investigate ground simulation test methods of two-dimensional (2-D) dynamic AE.
An aircraft ground simulation test usually uses a centrifuge to emulate a sustained AE. Kvrgic et al. [6] studied the control algorithm of a centrifuge motion simulator for pilot training. The simulator could simulate acceleration forces that act on aircraft pilots through rotations about axes. Lin et al. [7] used accelerated centrifugal force to simulate G-forces during aerial combat, and studied stress-influence factors in the centrifugal training for military aircrews. But the centrifuge cannot realize the synchronous controllable loading of 2-D dynamic acceleration, so it is unable to simulate a flight AE. A composite system of a centrifuge and a test turntable can simulate the flight AE. Test turntable is an important equipment for research on inertia components and has been widely utilized in aerospace [8]. Guan et al. [9] used a double turntable centrifuge (DTC) to test rectification error due to the composite accelerations. The DTC was utilized to produce the constant acceleration and vibration simultaneously. However, the composite system cannot be started instantaneously, thereby failing to simulate an impact AE during the launch phase.
In general, a typical dynamic impact simulation device is composed of a Hopkinson bar, a drop tower, a centrifuge, and an air cannon. The live-fire tests include artillery launches and air drops. The acceleration amplitude, pulse width, and frequency characteristics generated by each acceleration simulation test method vary, so the application scope of each method is different. Through using impact loading tests to simulate the stress environments, the strain characteristics and impact resistances of materials can be studied. Douglas et al. [10] analyzed the performance of wafer-level chips by using board-level drop impacts to produce accelerations between 10,000 g to 30,000 g with multiple impacts at the impact point, and the acceleration amplitude and mounting clearance were varied. By modifying Taylor's theory and optimizing the experimental models and parameter evaluation [11], Taylor's rod experiments have achieved a high range of strain rates. Furthermore, Taylor's rod experiments are simple, efficient, and reproducible, which become a common method of investigating material properties under high-strain-rate impact loads [12]. Brünig and Driemeier [13] proposed a specific thermo-hyperplastic constitutive model function to examine the performances of corresponding materials and numerically simulated the Taylor impact test. Shi et al. [14] designed a Hopkinson bar calibration device and studied a dynamic linear calibration method for high-g accelerometers. Using high-pressure gas as the power source and controlling the air pressure, an air cannon controlled the launch acceleration. Using the air cannon to assist the loading, it became possible to simulate a high acceleration load environment and the explosion impact environment, thereby allowing environmental testing on the structures of spacecraft and instrumentation.
In reality, the live-fire test usually requires high cost and a long cycle, which can hardly be repeated. The semi-physical simulation test can establish a similar environment in the laboratory to simulate an application environment of the aircraft. The existing AE simulation test methods generally adopt impact and flight overload simulation separately, which cannot realize the synchronous controllable loading of 2-D dynamic overload and cannot compare the real simulated impact and flight overload environment. Therefore, the objective of this paper is to propose a feasible 2-D dynamic AE simulation test method. Specifically, we focus on the key problems of loading in the impact overload environment, i.e., orderly transition from impact overload to steady-state overload and synchronous loading of 2-D dynamic AE. A specialized simulation test platform using an air cannon to start the compound centrifuge turntable is developed. Using the air cannon as the load propulsion system, the centrifuge provides a continuous AE. With vector distributions obtained by the single-axis turntable, synchronous loading of 2-D dynamic AE is implemented. This study contributes to research and development of extreme impact simulation tests on the structures of crystal oscillator circuits, acceleration sensors, and related instruments in dynamic AE.
The rest of the paper is organized as follows. In Section 2, the principle of the air-cannon-started compound centrifugal turntable simulation test method, the structure of the simulation device, and the simulation test process are introduced. In Section 3, the theoretical calculation, simulation, and testing of the simulation process in the impact overload environment are described. By controlling the pressure of the air cannon chamber, engaging the electromagnetic clutch, and adjusting the speed of the centrifuge, an orderly transition from the impact overload to the steady-state overload is achieved. In Section 4, by implementing a feed-forward proportional-integral-derivative (PID) control algorithm and the position-velocity-time (PVT) motion mode, controlling the angular displacement of the unit during testing, and applying a vector distribution of 2-D accelerations, synchronous loading of 2-D dynamic accelerations is achieved. In Section 5, a simulation test system is created to simulate 2-D AE. Finally, Section 6 gives concluding remarks.

Simulation Test Method with Air-Cannon-Started Compound Centrifugal Turntable
The requirements for the simulation tests are defined as follows. The acceleration process of the unit being tested includes two sequential stages: launch and maneuvering flight. After the system receives the launch order, the axial acceleration of the unit being tested must reach the set value within 0.2 s, and the radial acceleration of the unit being tested must be zero g during the launching stage. After the system receives the flight order, the radial acceleration of the unit should reach the set values within 0.5 s while the axial acceleration of the unit being tested remains unchanged during the flight stage.
The block diagram of the air-cannon-started compound centrifugal turntable simulation test system is illustrated in Figure 1. It mainly consists of a console, a speed control system, an electromagnetic clutch, a 2-D acceleration vector distribution system, an impact loading system, and a measurement system. The console is responsible for control and test of the entire system. The speed control system provides a continuous AE. The unit being tested and the vector distributor are installed on the main arm of the centrifuge. The centrifuge is controlled by the speed control system and the electromagnetic clutch. The 2-D acceleration vector distribution system is made up of a vector distribution controller, a driver, and a distributor. With a servo motor and the turntable used as the distributer, the resultant acceleration of the main arm of the centrifuge is distributed to the axial and radial directions of the unit being tested. The impact loading system consists of an impact loading controller, an impact loading driver, and an impact loader. An air cannon is employed as the impact loader to propel the main arm of the centrifuge for instantaneous acceleration. The measurement system measures the 2-D acceleration of the unit being tested and tests the working state of the unit being tested.
Appl. Sci. 2020, 10, 910 3 of 18 synchronous loading of 2-D dynamic accelerations is achieved. In Section 5, a simulation test system is created to simulate 2-D AE. Finally, Section 6 gives concluding remarks.

Simulation Test Method with Air-Cannon-Started Compound Centrifugal Turntable
The requirements for the simulation tests are defined as follows. The acceleration process of the unit being tested includes two sequential stages: launch and maneuvering flight. After the system receives the launch order, the axial acceleration of the unit being tested must reach the set value within 0.2 s, and the radial acceleration of the unit being tested must be zero g during the launching stage. After the system receives the flight order, the radial acceleration of the unit should reach the set values within 0.5 s while the axial acceleration of the unit being tested remains unchanged during the flight stage.
The block diagram of the air-cannon-started compound centrifugal turntable simulation test system is illustrated in Figure 1. It mainly consists of a console, a speed control system, an electromagnetic clutch, a 2-D acceleration vector distribution system, an impact loading system, and a measurement system. The console is responsible for control and test of the entire system. The speed control system provides a continuous AE. The unit being tested and the vector distributor are installed on the main arm of the centrifuge. The centrifuge is controlled by the speed control system and the electromagnetic clutch. The 2-D acceleration vector distribution system is made up of a vector distribution controller, a driver, and a distributor. With a servo motor and the turntable used as the distributer, the resultant acceleration of the main arm of the centrifuge is distributed to the axial and radial directions of the unit being tested. The impact loading system consists of an impact loading controller, an impact loading driver, and an impact loader. An air cannon is employed as the impact loader to propel the main arm of the centrifuge for instantaneous acceleration. The measurement system measures the 2-D acceleration of the unit being tested and tests the working state of the unit being tested.

Vector distribution controller
Vector distribution driver Vector distributor The unit being tested Impact loading controller Impact loading driver Impact loader Main arm of the centrifuge Figure 1. Air-cannon-started compound centrifugal turntable simulation test system.

Simulation of Impact Overload Environment
The simulation test of the impact overload environment requires the axial acceleration of the unit being tested to reach the set value within 0.2 s, which is the typical characteristics of instantaneous and high impacts. When the unit being tested is driven by the centrifuge, its axial direction will point toward the axis of the arm. The normal and tangential directions of the arm are the axial and radial directions of the unit being tested. When the centrifuge accelerates, the normal acceleration of the arm increases, and the tangential acceleration fluctuates. When the centrifuge is running at a constant speed, the normal acceleration of the arm is constant and the tangential acceleration is zero g, which meet the requirements of the acceleration amplitude of the impact loading environment. However, the acceleration time lasts nearly 1 s, failing to meet the impact loading requirements. To tackle this problem, we propose a simulation method for the impact loading environment using an air cannon booster and a motor relay.

Simulation of Impact Overload Environment
The simulation test of the impact overload environment requires the axial acceleration of the unit being tested to reach the set value within 0.2 s, which is the typical characteristics of instantaneous and high impacts. When the unit being tested is driven by the centrifuge, its axial direction will point toward the axis of the arm. The normal and tangential directions of the arm are the axial and radial directions of the unit being tested. When the centrifuge accelerates, the normal acceleration of the arm increases, and the tangential acceleration fluctuates. When the centrifuge is running at a constant speed, the normal acceleration of the arm is constant and the tangential acceleration is zero g, which meet the requirements of the acceleration amplitude of the impact loading environment. However, the acceleration time lasts nearly 1 s, failing to meet the impact loading requirements. To tackle this problem, we propose a simulation method for the impact loading environment using an air cannon booster and a motor relay.

Impact Process of Air Cannon
The structure of the adopted air cannon is shown in Figure 2. The air cannon consists of a back-pressure chamber, a high-pressure chamber, a quick-opening chamber, a control valve, a piston, and a piston rod. Initially, the pressure in the quick-opening chamber is higher than the pressure in the high-pressure chamber. Thus, the high-pressure chamber is sealed by the piston, and the piston rod remains stationary. When the control valve is opened, the pressure in the quick-opening chamber drops sharply, and the piston moves to the right quickly. Under this circumstance, the high-pressure chamber is opened, the high pressure gas drives the piston rod to move to the left, thus the air cannon is fired.
Appl. Sci. 2020, 10, 910 4 of 18 The structure of the adopted air cannon is shown in Figure 2. The air cannon consists of a backpressure chamber, a high-pressure chamber, a quick-opening chamber, a control valve, a piston, and a piston rod. Initially, the pressure in the quick-opening chamber is higher than the pressure in the high-pressure chamber. Thus, the high-pressure chamber is sealed by the piston, and the piston rod remains stationary. When the control valve is opened, the pressure in the quick-opening chamber drops sharply, and the piston moves to the right quickly. Under this circumstance, the high-pressure chamber is opened, the high pressure gas drives the piston rod to move to the left, thus the air cannon is fired. The resultant force FP on the piston rod and arm and the thrust FArm on the arm are expressed as follows: where s1 = 0.02 m 2 , which is the area of the piston; s2 = 0.0188 m 2 , which is the area of the piston minus the area of the piston rod; P1 denotes the air pressure in the high-pressure chamber during the launch; P2 denotes the air pressure in the back pressure chamber during the launch; A1 is the acceleration of the piston rod in the propulsion process; θ is the rotation angle of the arm in the propulsion process; mP = 44.7 kg, which is the equivalent load of the piston rod; L is the distance between the propulsion point and the axis of the arm; x is the displacement of the piston rod; and A2 is the tangential acceleration of the arm during the propulsion process. The motion of the piston pushed by the air is assumed to follow an adiabatic expansion process of an ideal gas. The equation for the state of the gas is described below.
where V1 denotes the volume of the high-pressure chamber during the launch, M indicates the mass of the air, μg stands for the molar mass of the air, R is the gas constant, and T is the air temperature. Based on the polytropic relationship of the gas in the high-pressure chamber during the compression and expansion process [15], the polytropic index of the air during the charging and discharging processes of the pneumatic system are set to 1.2. Thus, the pressures of the high-pressure chamber and the back-pressure chamber are expressed as follows: The resultant force F P on the piston rod and arm and the thrust F Arm on the arm are expressed as follows: where s 1 = 0.02 m 2 , which is the area of the piston; s 2 = 0.0188 m 2 , which is the area of the piston minus the area of the piston rod; P 1 denotes the air pressure in the high-pressure chamber during the launch; P 2 denotes the air pressure in the back pressure chamber during the launch; A 1 is the acceleration of the piston rod in the propulsion process; θ is the rotation angle of the arm in the propulsion process; m P = 44.7 kg, which is the equivalent load of the piston rod; L is the distance between the propulsion point and the axis of the arm; x is the displacement of the piston rod; and A 2 is the tangential acceleration of the arm during the propulsion process. The motion of the piston pushed by the air is assumed to follow an adiabatic expansion process of an ideal gas. The equation for the state of the gas is described below.
where V 1 denotes the volume of the high-pressure chamber during the launch, M indicates the mass of the air, µ g stands for the molar mass of the air, R is the gas constant, and T is the air temperature. Based on the polytropic relationship of the gas in the high-pressure chamber during the compression and expansion process [15], the polytropic index of the air during the charging and discharging processes of the pneumatic system are set to 1.2. Thus, the pressures of the high-pressure chamber and the back-pressure chamber are expressed as follows: where P 10 is the initial pressure of the high-pressure chamber at the time of launch, V 10 = 0.0972 m 3 , which is the initial volume of the high-pressure chamber, P 20 is the initial pressure of the back-pressure chamber at the time of launch, and V 20 = 0.02827 m 3 , which is the initial volume of the back-pressure chamber. Substituting Equation (3) into Equation (1) yields where m Arm = 45.7 kg, which is the effective load of the arm at the impact point, including the equivalent mass of the arm, the unit being tested, and the vector distributer. As shown in Equation (4), the acceleration of the arm during the impact is determined by the initial air pressures and volumes of the back-pressure and high-pressure chambers, the areas of the piston and piston rod, the displacement of the piston rod while it is propelled, the propulsion location, and the effective loads of the arm and piston. In the tests, only the initial pressures of the back-pressure and high-pressure chamber are variables, so the initial air pressure of the high-pressure and back-pressure chambers can be adjusted to obtain different launch accelerations. As the displacement of the piston rod increases, the thrust becomes increasingly small, and the piston rod can reverse directions. To avoid reverse collision between the arm and the recoiling air cannon, the piston rod and the arm are not placed at the same level. The collision end of the piston rod is machined into a pulley. The pulley rises before and falls after the collision. When simulating an impact overload environment of 10 g, the simulated air pressures of the high-pressure and back-pressure chambers versus the displacement of the piston rod can be obtained. Figure 3 plots the relationship between the pressure and the displacement of the piston rod.
Appl. Sci. 2020, 10, 910 5 of 18 where P10 is the initial pressure of the high-pressure chamber at the time of launch, V10 = 0.0972 m 3 , which is the initial volume of the high-pressure chamber, P20 is the initial pressure of the backpressure chamber at the time of launch, and V20=0.02827 m 3 , which is the initial volume of the backpressure chamber.
Substituting Equation (3) into Equation (1) yields where mArm = 45.7 kg, which is the effective load of the arm at the impact point, including the equivalent mass of the arm, the unit being tested, and the vector distributer. As shown in Equation (4), the acceleration of the arm during the impact is determined by the initial air pressures and volumes of the back-pressure and high-pressure chambers, the areas of the piston and piston rod, the displacement of the piston rod while it is propelled, the propulsion location, and the effective loads of the arm and piston. In the tests, only the initial pressures of the back-pressure and high-pressure chamber are variables, so the initial air pressure of the highpressure and back-pressure chambers can be adjusted to obtain different launch accelerations. As the displacement of the piston rod increases, the thrust becomes increasingly small, and the piston rod can reverse directions. To avoid reverse collision between the arm and the recoiling air cannon, the piston rod and the arm are not placed at the same level. The collision end of the piston rod is machined into a pulley. The pulley rises before and falls after the collision. When simulating an impact overload environment of 10 g, the simulated air pressures of the high-pressure and backpressure chambers versus the displacement of the piston rod can be obtained. Figure 3 plots the relationship between the pressure and the displacement of the piston rod. At the displacement of 0.5 m, the air cannon is fired, which is controlled by the air pressures of the high-pressure and back-pressure chambers. The acceleration of the pulley and the tangential acceleration of the arm are plotted in Figure 4. As can be observed, with an impact overload of 10 g, 2 ms after the air cannon launching, the accelerations peak 10.74 g and then drop rapidly. The accelerations remain between 6 and 7 g for about 30 ms, after which they decrease. The entire At the displacement of 0.5 m, the air cannon is fired, which is controlled by the air pressures of the high-pressure and back-pressure chambers. The acceleration of the pulley and the tangential acceleration of the arm are plotted in Figure 4. As can be observed, with an impact overload of 10 g, 2 ms after the air cannon launching, the accelerations peak 10.74 g and then drop rapidly. The accelerations remain between 6 and 7 g for about 30 ms, after which they decrease. The entire propulsion time lasts 119.5 ms. At the end of the launching, the tangential acceleration of the arm is zero g.
Appl. Sci. 2020, 10, 910 6 of 18 propulsion time lasts 119.5 ms. At the end of the launching, the tangential acceleration of the arm is zero g.

Clutch Engagement Process
During the simulation analysis, the air cannon pushes the arm. The starting time is very short, the acceleration amplitude and loading time meet the simulation test requirements of the impact overload environment. The air cannon propulsion duration is about 120 ms, but the acceleration rise time is only 2 ms. The duration of the peak is too short. To maintain the impact state, the clutch needs to be engaged and the motor needs to be relayed to power the arm rotation. If the relay is completed at the peak of the acceleration, the acceleration of the arm can be achieved.
The driving plate of the clutch is driven by the main control motor. The load on the driven plate includes the effective loads of the reducer, the rotating shaft, the arm, the vector distributor, and the unit being tested. The power transmission system of the centrifuge is shown in Figure 5.
where Jc, ωc, and Tc are the moment of inertia, revolution speed, and output torque of the main control motor, respectively, Tf is the friction torque transmitted by the clutch, and Jd, ωd and Td are the effective moment of inertia, revolution speed, equivalent resistance torque of the driven plate, respectively. When the clutch is completely disengaged, Tf = 0 Nm, and when the clutch is fully engaged, Tf = Tc = Td. The main operation process of the clutch in the test system varies from a fully engaged sliding friction state to a fully disengaged sliding friction state. The initial speed of the driven plate is zero rad/s. When the clutch is engaged, the speeds of the driven plate are shown in Figure 6.

Driving plate
Driven plate

Control motor
Effective loads

Clutch Engagement Process
During the simulation analysis, the air cannon pushes the arm. The starting time is very short, the acceleration amplitude and loading time meet the simulation test requirements of the impact overload environment. The air cannon propulsion duration is about 120 ms, but the acceleration rise time is only 2 ms. The duration of the peak is too short. To maintain the impact state, the clutch needs to be engaged and the motor needs to be relayed to power the arm rotation. If the relay is completed at the peak of the acceleration, the acceleration of the arm can be achieved.
The driving plate of the clutch is driven by the main control motor. The load on the driven plate includes the effective loads of the reducer, the rotating shaft, the arm, the vector distributor, and the unit being tested. The power transmission system of the centrifuge is shown in Figure 5.

Clutch Engagement Process
During the simulation analysis, the air cannon pushes the arm. The starting time is very short, the acceleration amplitude and loading time meet the simulation test requirements of the impact overload environment. The air cannon propulsion duration is about 120 ms, but the acceleration rise time is only 2 ms. The duration of the peak is too short. To maintain the impact state, the clutch needs to be engaged and the motor needs to be relayed to power the arm rotation. If the relay is completed at the peak of the acceleration, the acceleration of the arm can be achieved.
The driving plate of the clutch is driven by the main control motor. The load on the driven plate includes the effective loads of the reducer, the rotating shaft, the arm, the vector distributor, and the unit being tested. The power transmission system of the centrifuge is shown in Figure 5.
where Jc, ωc, and Tc are the moment of inertia, revolution speed, and output torque of the main control motor, respectively, Tf is the friction torque transmitted by the clutch, and Jd, ωd and Td are the effective moment of inertia, revolution speed, equivalent resistance torque of the driven plate, respectively. When the clutch is completely disengaged, Tf = 0 Nm, and when the clutch is fully engaged, Tf = Tc = Td. The main operation process of the clutch in the test system varies from a fully engaged sliding friction state to a fully disengaged sliding friction state. The initial speed of the driven plate is zero rad/s. When the clutch is engaged, the speeds of the driven plate are shown in Figure 6.  .
where J c , ω c , and T c are the moment of inertia, revolution speed, and output torque of the main control motor, respectively, T f is the friction torque transmitted by the clutch, and J d , ω d and T d are the effective moment of inertia, revolution speed, equivalent resistance torque of the driven plate, respectively. When the clutch is completely disengaged, T f = 0 Nm, and when the clutch is fully engaged, T f = T c = T d . The main operation process of the clutch in the test system varies from a fully engaged sliding Appl. Sci. 2020, 10, 910 7 of 18 friction state to a fully disengaged sliding friction state. The initial speed of the driven plate is zero rad/s. When the clutch is engaged, the speeds of the driven plate are shown in Figure 6. At the beginning, the driven plate keeps stationary, the clutch is disengaged, and the motor is idling. Therefore, Td, Tf, and Tc are all zero Nm, and the motor's speed ωc is a constant. When the clutch is engaged, the difference in the speeds of revolution becomes large. Tf is gradually increased, the revolution speed (ωd) of the driven plate increases, and the output torque (Tc) of the motor also increases. Because Tf is less than Tc, the motor accelerates. The difference in the speeds of revolution increases. Tf increases faster than Tc, and the driven plate accelerates. When Tf > Tc, the motor decelerates, the difference in the revolution speeds becomes increasingly small, and Tf gradually decreases. When Tf = Tc, the clutch is fully engaged, the motor drives the load to accelerate and stabilize after attaining the set revolution speed.
In the simulation tests, the arm is started by the air cannon. Ideally, the engagement process will be completed when the propelled acceleration reaches the maximum, and the tangential acceleration of the arm after launching will be equal to the set normal acceleration of the arm. Since the propulsion time is very short, ωd should be less than ωc, and the arm always accelerates during the engagement process. If the timing is missed, the engagement occurs too late, the arm decelerates, the impact acceleration is less than zero g, and the impact overload cannot be simulated. Therefore, it is necessary to select a suitable clutch and proper engagement timing. By controlling the pressures of the air cannon chambers and the revolution speed of the motor, the rapid transition between the impact overload and the continuous acceleration can be achieved.

Speed Control System Based on Active Disturbance Rejection Control
When the clutch is engaged, the output torque of the motor is increased rapidly, and the serious load disturbance is produced. In order to improve the precision of the speed control system, it becomes necessary to compensate for the disturbance. Active disturbance rejection control (ADRC) considers perturbations to the model and parameters and unmeasurable external disturbances as the total disturbance of the system. Based on the feedback control of the system error, the total disturbance is observed and corrected, which alleviates the conflict between the response speed and overshoot in the PID controls. ADRC does not rely on any precise mathematical model of a system, and it significantly improves the dynamic response characteristics of the controller through a proper nonlinear combination. The load disturbance can be estimated and compensated by a speed controller based on active disturbance rejection control (ADRC). Many researchers have studied the application of the speed controller based on ADRC. In order to improve pointing accuracy and rotation speed, Li et al. [16] designed the inner and outer loops of an ADRC for the antenna pointing control of a large flexible satellite system. Kang et al. [17] designed an idle speed controller based on ADRC to compensate for varying engine load and friction torque in passenger car diesel engines. At the beginning, the driven plate keeps stationary, the clutch is disengaged, and the motor is idling. Therefore, T d , T f , and T c are all zero Nm, and the motor's speed ω c is a constant. When the clutch is engaged, the difference in the speeds of revolution becomes large. T f is gradually increased, the revolution speed (ω d ) of the driven plate increases, and the output torque (T c ) of the motor also increases. Because T f is less than T c , the motor accelerates. The difference in the speeds of revolution increases. T f increases faster than T c , and the driven plate accelerates. When T f > T c , the motor decelerates, the difference in the revolution speeds becomes increasingly small, and T f gradually decreases. When T f = T c , the clutch is fully engaged, the motor drives the load to accelerate and stabilize after attaining the set revolution speed.
In the simulation tests, the arm is started by the air cannon. Ideally, the engagement process will be completed when the propelled acceleration reaches the maximum, and the tangential acceleration of the arm after launching will be equal to the set normal acceleration of the arm. Since the propulsion time is very short, ω d should be less than ω c , and the arm always accelerates during the engagement process. If the timing is missed, the engagement occurs too late, the arm decelerates, the impact acceleration is less than zero g, and the impact overload cannot be simulated. Therefore, it is necessary to select a suitable clutch and proper engagement timing. By controlling the pressures of the air cannon chambers and the revolution speed of the motor, the rapid transition between the impact overload and the continuous acceleration can be achieved.

Speed Control System Based on Active Disturbance Rejection Control
When the clutch is engaged, the output torque of the motor is increased rapidly, and the serious load disturbance is produced. In order to improve the precision of the speed control system, it becomes necessary to compensate for the disturbance. Active disturbance rejection control (ADRC) considers perturbations to the model and parameters and unmeasurable external disturbances as the total disturbance of the system. Based on the feedback control of the system error, the total disturbance is observed and corrected, which alleviates the conflict between the response speed and overshoot in the PID controls. ADRC does not rely on any precise mathematical model of a system, and it significantly improves the dynamic response characteristics of the controller through a proper nonlinear combination. The load disturbance can be estimated and compensated by a speed controller based on active disturbance rejection control (ADRC). Many researchers have studied the application of the speed controller based on ADRC. In order to improve pointing accuracy and rotation speed, Li et al. [16] designed the inner and outer loops of an ADRC for the antenna pointing control of a large flexible satellite system. Kang et al. [17] designed an idle speed controller based on ADRC to compensate for varying engine load and friction torque in passenger car diesel engines.
ADRC consists of three main components: a tracking differentiator (TD), an extended state observer (ESO), and a nonlinear state error feedback law (NSEFL). The TD arranges the transition process for the system and extracts the differential signal. ESO is the core of ADRC and can estimate the total disturbance of the system in real time and compensate for it. The outputs of the TD and ESO are compared to obtain the systematic error. Based on the error signal, the NSEFL produces a control strategy for the controlled object.
Based on the torque and motion equations of the three-phase asynchronous motor, the differential equation of the speed loop can be obtained as follows: .
where ω is the rotor speed, ρ is the number of the pole pairs of the motor, L m is the stator-rotor mutual inductance, J is the system moment of inertia converted to the motor shaft, L r is the rotor inductance, ψ r is the rotor magnetic flux, i Ts is the torque component of the stator current, and T L is the load torque.
The nonlinear state error feedback control law is The control quantity i Ts obtained by the nonlinear state error feedback control can be written as Substituting Equation (11) into Equation (7) yields the following formula After the ADRC is adjusted, the coupling term is no longer included in the speed state equation and the total disturbance of the system is canceled. The structure of the speed outer loop ADRC system is illustrated in Figure 7. Similarly, the first-order ADRC regulations for the excitation current loop and the torque current loop can also be obtained. The differential equation for the rotor flux is shown as follows: where Rr is the rotor resistance and iMs is the excitation component of the stator current. Equation (13) does not contain a coupling term, and the rotor flux outer loop can be regulated by a PI controller. The first-order ADRC regulators are applied to the speed outer loop, the excitation current loop, and the torque current loop. Thus, the asynchronous motor speed control system based on the ADRC is designed. Simulation and analysis of the speed control system based on ADRC are conducted. The parameters of the asynchronous motor are shown in Table 1. The appropriate parameters of the speed control system should be set before the simulation. The discretization time values are all set to 100 ms in this system. The regulator coefficients of the speed outer loop are: β21 = 1000, β22 = 16,000, a3 = 0.5, δ3 = 0.01, β3 = 10. The regulator coefficients of the excitation current loop are: β01 = 10,000, β02 = 120,000, a1 = 0.5, δ1 = 0.01, β1 = 100. The regulator coefficients of the torque current loop are: β11 = 10,000, β12 = 120,000, a2 = 0.5, δ2 = 0.01, β2 = 100. The regulator coefficients of the rotor flux outer loop are: Kp = 2, TI = 10.  Similarly, the first-order ADRC regulations for the excitation current loop and the torque current loop can also be obtained. The differential equation for the rotor flux is shown as follows: where R r is the rotor resistance and i Ms is the excitation component of the stator current. Equation (13) does not contain a coupling term, and the rotor flux outer loop can be regulated by a PI controller. The first-order ADRC regulators are applied to the speed outer loop, the excitation current loop, and the torque current loop. Thus, the asynchronous motor speed control system based on the ADRC is designed. Simulation and analysis of the speed control system based on ADRC are conducted. The parameters of the asynchronous motor are shown in Table 1. The appropriate parameters of the speed control system should be set before the simulation. The discretization time values are all set to 100 ms in this system. The regulator coefficients of the speed outer loop are: β 21 = 1000, β 22 = 16,000, a 3 = 0.5, δ 3 = 0.01, β 3 = 10. The regulator coefficients of the excitation current loop are: β 01 = 10,000, β 02 = 120,000, a 1 = 0.5, δ 1 = 0.01, β 1 = 100. The regulator coefficients of the torque current loop are: β 11 = 10,000, β 12 = 120,000, a 2 = 0.5, δ 2 = 0.01, β 2 = 100. The regulator coefficients of the rotor flux outer loop are: K p = 2, T I = 10. In the simulation test of the system operation, the motor is started without a load. When the clutch is engaged, the load disturbance is produced, and the torque transmitted depends on the speed difference between the driving plate and the driven plate.
The load torque is shown as follows: where J d = 10 Kgm 2 , which is the load moment of inertia, α = 3 rad/s 2 , which is the load angular acceleration. The set motor speed is 200 rpm, and the motor relay is set as 0.4 s with a load of 30 Nm added abruptly. The simulation results of the motor speed controlled by ADRC are shown in Figures 8  and 9. As can be observed, the system response does not exhibit overshoot and the rise time is 0.037 s. When the load is suddenly applied at 0.4 s, the motor speed drops to 196.9 rpm. After 0.01 s, the speed quickly returns to a stable state. The load is driven by the air cannon before the clutch is engaged, and it is driven by the motor after the clutch is engaged. The load is rotating in the same direction as the motor and the speed difference between the load and the motor is small, so the change of the speed will become small and the recovery time will get short. Accompanied with the improved disturbance-rejection performance of the system, an orderly transition between the impact overload and steady-state overload can be achieved.
In the simulation test of the system operation, the motor is started without a load. When the clutch is engaged, the load disturbance is produced, and the torque transmitted depends on the speed difference between the driving plate and the driven plate.
The load torque is shown as follows: where Jd = 10 Kgm 2 , which is the load moment of inertia, α = 3 rad/s 2 , which is the load angular acceleration. The set motor speed is 200 rpm, and the motor relay is set as 0.4 s with a load of 30 Nm added abruptly. The simulation results of the motor speed controlled by ADRC are shown in Figures 8 and  9. As can be observed, the system response does not exhibit overshoot and the rise time is 0.037 s. When the load is suddenly applied at 0.4 s, the motor speed drops to 196.9 rpm. After 0.01 s, the speed quickly returns to a stable state. The load is driven by the air cannon before the clutch is engaged, and it is driven by the motor after the clutch is engaged. The load is rotating in the same direction as the motor and the speed difference between the load and the motor is small, so the change of the speed will become small and the recovery time will get short. Accompanied with the improved disturbance-rejection performance of the system, an orderly transition between the impact overload and steady-state overload can be achieved.

Loading Test of Impact Overload Environment
Both the air cannon firing and the electromagnetic clutch engagement commands are sent by the digital output module to the driving system controlled by the relay. The electromagnetic clutch is controlled by the carbon brush, and the air cannon is controlled by the electric solenoid valve. The response speeds of the two mechanisms are different.
The revolution speed of the arm is set as 105 rpm. After the air cannon is fired, the electromagnetic clutch is immediately relayed. The measured impact acceleration of the arm is shown in Figure 10. The steady-state normal acceleration of the arm is 11 g. The acceleration rise

Loading Test of Impact Overload Environment
Both the air cannon firing and the electromagnetic clutch engagement commands are sent by the digital output module to the driving system controlled by the relay. The electromagnetic clutch is controlled by the carbon brush, and the air cannon is controlled by the electric solenoid valve. The response speeds of the two mechanisms are different.
The revolution speed of the arm is set as 105 rpm. After the air cannon is fired, the electromagnetic clutch is immediately relayed. The measured impact acceleration of the arm is shown in Figure 10. The steady-state normal acceleration of the arm is 11 g. The acceleration rise time is 74.56 ms. The loading time of the launch acceleration is 819 ms, which does not meet the simulation test requirements for the impact overload environment. The electromagnetic clutch has a longer response time than the solenoid valve, causing too long loading time. As can be measured, the electromagnetic clutch has a response time that is 250 ms slower than that of the electric solenoid valve. Therefore, the loading time of the impact acceleration can be regulated by adjusting the operating time sequence of the air cannon launch and the clutch engagement. After the clutch has been engaged for 250 ms, the air cannon is fired.

Loading Test of Impact Overload Environment
Both the air cannon firing and the electromagnetic clutch engagement commands are sent by the digital output module to the driving system controlled by the relay. The electromagnetic clutch is controlled by the carbon brush, and the air cannon is controlled by the electric solenoid valve. The response speeds of the two mechanisms are different.
The revolution speed of the arm is set as 105 rpm. After the air cannon is fired, the electromagnetic clutch is immediately relayed. The measured impact acceleration of the arm is shown in Figure 10. The steady-state normal acceleration of the arm is 11 g. The acceleration rise time is 74.56 ms. The loading time of the launch acceleration is 819 ms, which does not meet the simulation test requirements for the impact overload environment. The electromagnetic clutch has a longer response time than the solenoid valve, causing too long loading time. As can be measured, the electromagnetic clutch has a response time that is 250 ms slower than that of the electric solenoid valve. Therefore, the loading time of the impact acceleration can be regulated by adjusting the operating time sequence of the air cannon launch and the clutch engagement. After the clutch has been engaged for 250 ms, the air cannon is fired. The unit being tested is fixed on the arm. The axial direction of the unit being tested points to the normal direction of the arm, and the radial direction of the unit being tested points to the tangential direction of the arm. So, the tangential acceleration of the arm is the radial acceleration of the unit, and the normal acceleration of the arm is the axial acceleration of the unit. The impact overload acceleration is set as 20 g. After 2.5 s of the air cannon being fired, the motor accelerates, and the axial acceleration of the unit being tested is set as 40 g. The 2-D acceleration test results of the unit being tested are shown in Figure 11. The air cannon propels the arm and produces a powerful impact. The radial acceleration of the unit being tested rises instantaneously. After the clutch is engaged, the motor drives the arm to rotate at a constant revolution speed, while the radial acceleration of the unit being tested is zero g. After 2.5 s, the motor is accelerated, and the axial acceleration of the unit being tested reaches 40 g. The unit being tested is fixed on the arm. The axial direction of the unit being tested points to the normal direction of the arm, and the radial direction of the unit being tested points to the tangential direction of the arm. So, the tangential acceleration of the arm is the radial acceleration of the unit, and the normal acceleration of the arm is the axial acceleration of the unit. The impact overload acceleration is set as 20 g. After 2.5 s of the air cannon being fired, the motor accelerates, and the axial acceleration of the unit being tested is set as 40 g. The 2-D acceleration test results of the unit being tested are shown in Figure 11. The air cannon propels the arm and produces a powerful impact. The radial acceleration of the unit being tested rises instantaneously. After the clutch is engaged, the motor drives the arm to rotate at a constant revolution speed, while the radial acceleration of the unit being tested is zero g. After 2.5 s, the motor is accelerated, and the axial acceleration of the unit being tested reaches 40 g. As shown in Figure 11, the loading time of the impact overload acceleration is 120 ms, and the loading time of the flight overload acceleration is 335 ms, which satisfy the speed requirements of the axial acceleration of the unit being tested. However, during the acceleration of the arm, the tangential acceleration occurs, which does not meet the requirements of the simulation test for the radial acceleration of the unit being tested. Therefore, the vector shaft angle needs to be adjusted to track or decompose the resultant acceleration of the arm to achieve dynamically controllable loading of the As shown in Figure 11, the loading time of the impact overload acceleration is 120 ms, and the loading time of the flight overload acceleration is 335 ms, which satisfy the speed requirements of the axial acceleration of the unit being tested. However, during the acceleration of the arm, the tangential acceleration occurs, which does not meet the requirements of the simulation test for the radial acceleration of the unit being tested. Therefore, the vector shaft angle needs to be adjusted to track or decompose the resultant acceleration of the arm to achieve dynamically controllable loading of the radial acceleration of the unit being tested.

Synchronous Loading of 2-D Dynamic AE
As shown in Figure 11, the radial acceleration of the unit being tested is generated during the acceleration process. The shorter the acceleration duration is, the greater the radial acceleration amplitude will become. The radial acceleration is zero g at the stable and constant speed. The simulation test requires the axial and radial accelerations to be simultaneously applied to the unit being tested with controllable amplitudes. To this end, we propose a loading method of the centrifugal turntable AE based on the 2-D vector distribution. To implement the synchronous loading of the 2-D dynamic AE, when the centrifuge is accelerating, the vector axis rotation angle of the unit being tested is adjusted by the 2-D vector distributor. The tangential and normal accelerations generated by the centrifuge are then combined or decomposed to the axial direction and radial direction of the unit being tested. Therefore, the 2-D dynamic vector distribution for the accelerations can be obtained.

Centrifugal Turntable AE Loading System Based on 2-D Vector Distribution
The centrifugal turntable simulation test system includes an arm control subsystem, a 2-D vector distribution subsystem, and an acceleration test subsystem. The 2-D vector distribution subsystem and the acceleration test subsystem are mounted on the arm. The arm is driven by the variable-frequency adjustable-speed motor. Both the unit being tested and the accelerometer are fixed on the turntable. The turntable is driven by the servo motor. The acceleration of the unit being tested is fed back to the computer through the accelerometer. The control and test signals of the arm are transmitted through the conductive slip ring.
Based on the axial and radial acceleration requirements of the unit being tested, the combined acceleration of the arm can be calculated, and the rotation of the arm can be controlled. When the arm is accelerated, the servo motor controls the angular displacement of the turntable and the unit being tested, and the resultant acceleration of the arm can be decomposed along the axial and radial directions of the unit being tested to achieve the synchronous loading of the 2-D dynamic AE. If the axial direction of the unit being tested is along the direction of the resultant acceleration of the arm, the radial acceleration of the unit being tested is zero g, which simulates the impact overload environment. By decomposing the resultant acceleration of the arm in two dimensions and keeping the axial acceleration of the unit being tested constant, the flight overload environment can be simulated.

Control of Servo Tracking System for Angular Position of Unit Being Tested
Because of the elastic deformation of the transmission mechanism, the servo system of the axial angular position of the unit being tested is a very complicated high-order system. The servo motor drives the turntable through the conveyor belt, which negatively affects the dynamic response characteristics of the system. To suppress the interference caused by load changes, measures are required to form fully closed-loop control. In addition, control strategies and system rigidity both needs to be improved.
Feed-forward control is an open-loop control strategy for a certain disturbance parameter. It can make use of the characteristics of the object to improve the response speed of the system. However, it does not compensate for other disturbances. This problem can be solved by a closed-loop system using a compound control strategy with a combination of feed-forward and feedback control. The compound control strategy can not only reduce the dynamic influence of the disturbance on the controlled parameters by feed-forward control, but also reduce the steady error of the system by feedback control. Many researchers have studied the application of the compound control strategy based on feed-forward compensation in the position servo control system. To compensate for system nonlinearity, Sarkar et al. [18] applied feed-forward PID tracking controls to the design of an electro-hydraulic system. Fathi et al. [19] used a feed-forward PID compound controller to design a vision-based closed-loop laser solid freeform fabrication system. As shown in Figure 12, the position compound controller uses a control strategy with speed/acceleration feed-forward compensation and integral separation PID. The speed/acceleration feed-forward compensation makes the controlled signals of the system insensitive to external disturbances and changes in the internal parameters of the system. The speed feed-forward controller can reduce the tracking error caused by damping. The acceleration feed-forward controller can reduce or eliminate the tracking error caused by the system inertia, thus improving the disturbance-rejection performance of the system. The parameters of the feed-forward PID controller includes a proportional coefficient (K P ), derivative coefficient (K d ), integral coefficient (K i ), velocity feed-forward coefficient (K vff ), acceleration feed-forward coefficient (K aff ), and integral separation coefficient K l . When the error is less than the set value, K l = 1 and the integral loop nodes take effect; otherwise, K l = 0 and the integral loop nodes do not take effect, the position error is no longer integrated, and the integral saturation is suppressed.
needs to be improved.
Feed-forward control is an open-loop control strategy for a certain disturbance parameter. It can make use of the characteristics of the object to improve the response speed of the system. However, it does not compensate for other disturbances. This problem can be solved by a closed-loop system using a compound control strategy with a combination of feed-forward and feedback control. The compound control strategy can not only reduce the dynamic influence of the disturbance on the controlled parameters by feed-forward control, but also reduce the steady error of the system by feedback control. Many researchers have studied the application of the compound control strategy based on feed-forward compensation in the position servo control system. To compensate for system nonlinearity, Sarkar et al. [18] applied feed-forward PID tracking controls to the design of an electrohydraulic system. Fathi et al. [19] used a feed-forward PID compound controller to design a visionbased closed-loop laser solid freeform fabrication system.
As shown in Figure 12, the position compound controller uses a control strategy with speed/acceleration feed-forward compensation and integral separation PID. The speed/acceleration feed-forward compensation makes the controlled signals of the system insensitive to external disturbances and changes in the internal parameters of the system. The speed feed-forward controller can reduce the tracking error caused by damping. The acceleration feed-forward controller can reduce or eliminate the tracking error caused by the system inertia, thus improving the disturbancerejection performance of the system. The parameters of the feed-forward PID controller includes a proportional coefficient (KP), derivative coefficient (Kd), integral coefficient (Ki), velocity feedforward coefficient (Kvff), acceleration feed-forward coefficient (Kaff), and integral separation coefficient Kl. When the error is less than the set value, Kl = 1 and the integral loop nodes take effect; otherwise, Kl = 0 and the integral loop nodes do not take effect, the position error is no longer integrated, and the integral saturation is suppressed.  Consistent with Figure 12, the feed-forward compound PID control rules can be expressed as follows:  Figure 12, the feed-forward compound PID control rules can be expressed as follows:

Consistent with
where u 1(t) is the feed-forward output of the velocity and acceleration, u 2(t) is the output of the position loop, u 3(t) is the output of the velocity loop, u (t) is the output of the PID controller, cp (t) is the set position, v (t) is the actual angular speed of the servo motor, and e (t) is the error of the set position from the actual position.

Vector Axial Angular Displacement Motion Trajectory Planning
Hermite interpolation can be used to eliminate the issues that arise when the interpolation conditions contain derivatives. Therefore, the motions of the axial positions of the unit being tested are fitted using Hermite interpolation polynomials, which can guarantee that the displacement and velocity at the interpolation points are equal to those of the ideal motion. Each interpolation segment of the PVT motion mode is fitted with a two-point cubic Hermite interpolation polynomial. Only by specifying the position, velocity and running time of the segment and calculating the motion trajectory that meets the constraint conditions, the more compact control of the trajectory can be realized.
The angular displacement motions of the unit being tested shown in Figure 13 are divided into six segments. Three cubic Hermite interpolations in the different intervals are used to ensure the following point: the function values (displacement) and the first derivative values (velocity) at the interpolation points are exactly the same as the values of the motion trajectory. The velocity curve can be obtained by taking the derivative of the displacement curve. The cubic Hermite interpolation function corresponding to the first segment of the curve is calculated as follows: where h 0 (x), h 1 (x), H 0 (x), and H 1 (x) are the interpolation basis functions.
where u1(t) is the feed-forward output of the velocity and acceleration, u2(t) is the output of the position loop, u3(t) is the output of the velocity loop, u(t) is the output of the PID controller, cp(t) is the set position, v(t) is the actual angular speed of the servo motor, and e(t) is the error of the set position from the actual position.

Vector Axial Angular Displacement Motion Trajectory Planning
Hermite interpolation can be used to eliminate the issues that arise when the interpolation conditions contain derivatives. Therefore, the motions of the axial positions of the unit being tested are fitted using Hermite interpolation polynomials, which can guarantee that the displacement and velocity at the interpolation points are equal to those of the ideal motion. Each interpolation segment of the PVT motion mode is fitted with a two-point cubic Hermite interpolation polynomial. Only by specifying the position, velocity and running time of the segment and calculating the motion trajectory that meets the constraint conditions, the more compact control of the trajectory can be realized. The angular displacement motions of the unit being tested shown in Figure 13 are divided into six segments. Three cubic Hermite interpolations in the different intervals are used to ensure the following point: the function values (displacement) and the first derivative values (velocity) at the interpolation points are exactly the same as the values of the motion trajectory. The velocity curve can be obtained by taking the derivative of the displacement curve. The cubic Hermite interpolation function corresponding to the first segment of the curve is calculated as follows: The cubic Hermitian interpolation polynomial is as follows: Similarly, the cubic Hermitian interpolation polynomials corresponding to the other five segments can be obtained. The displacement curve after interpolation is depicted in Figure 14.
The cubic Hermitian interpolation polynomial is as follows: Similarly, the cubic Hermitian interpolation polynomials corresponding to the other five segments can be obtained. The displacement curve after interpolation is depicted in Figure 14. The cubic Hermitian interpolation must extract the first-order derivative of the motion curve at the interpolation point, which is smoother than the motion trajectory of the cubic spline interpolation and is suitable for S-shaped acceleration/deceleration controls of high-speed machining. As shown in Figure 11, the impact load is achieved in 0.12 s. The launch AE requires that the axial acceleration of the unit to be increased from zero g to 20 g in 0.12 s, and the radial acceleration is zero g. To simulate this AE, the initial position of the unit is vertical with the arm. During the acceleration, the servo motor controls the vector shaft rotation angle so that the axial direction of the unit is along the resultant acceleration direction of the arm. The shaft angular displacement-tracking curves based on the PVT mode are plotted in Figure 15. As can be observed, the unit is rotated 90° from the initial position in a pre-planned trajectory within 0.12 s. When the revolution speed of the arm reaches a steady state, the axial direction of the unit is along the normal direction of the arm, the axial acceleration is rω 2 , and the radial acceleration of the unit is maintained at rdω/dt = 0 g. The cubic Hermitian interpolation must extract the first-order derivative of the motion curve at the interpolation point, which is smoother than the motion trajectory of the cubic spline interpolation and is suitable for S-shaped acceleration/deceleration controls of high-speed machining. As shown in Figure 11, the impact load is achieved in 0.12 s. The launch AE requires that the axial acceleration of the unit to be increased from zero g to 20 g in 0.12 s, and the radial acceleration is zero g. To simulate this AE, the initial position of the unit is vertical with the arm. During the acceleration, the servo motor controls the vector shaft rotation angle so that the axial direction of the unit is along the resultant acceleration direction of the arm. The shaft angular displacement-tracking curves based on the PVT mode are plotted in Figure 15. As can be observed, the unit is rotated 90 • from the initial position in a pre-planned trajectory within 0.12 s. When the revolution speed of the arm reaches a steady state, the axial direction of the unit is along the normal direction of the arm, the axial acceleration is rω 2 , and the radial acceleration of the unit is maintained at rdω/dt = 0 g.  Figure 16 shows a photograph of the built simulation test system. Mechanically, the impact and flight 2-D dynamic AE was simulated using an air cannon propulsion, a motor relay, and a vector axial angle tracking. The system has been applied to the performance test of aerospace components. The load time and the amplitude of the 2-D acceleration were set according to the performance index of the unit being tested. The axial acceleration of the impact overload was set as 20 g and the radial acceleration was zero g. The axial acceleration of the flight overload was maintained at 20 g, while the radial acceleration was increased from zero g to 40 g within 0.5 s.  Figure 16 shows a photograph of the built simulation test system. Mechanically, the impact and flight 2-D dynamic AE was simulated using an air cannon propulsion, a motor relay, and a vector axial angle tracking. The system has been applied to the performance test of aerospace components. The load time and the amplitude of the 2-D acceleration were set according to the performance index of the unit being tested. The axial acceleration of the impact overload was set as 20 g and the radial acceleration was zero g. The axial acceleration of the flight overload was maintained at 20 g, while the radial acceleration was increased from zero g to 40 g within 0.5 s. flight 2-D dynamic AE was simulated using an air cannon propulsion, a motor relay, and a vector axial angle tracking. The system has been applied to the performance test of aerospace components. The load time and the amplitude of the 2-D acceleration were set according to the performance index of the unit being tested. The axial acceleration of the impact overload was set as 20 g and the radial acceleration was zero g. The axial acceleration of the flight overload was maintained at 20 g, while the radial acceleration was increased from zero g to 40 g within 0.5 s. The measured accelerations of the unit being tested with the synchronous loading are plotted in Figure 17. The simulation test process involved the following six steps. The measured accelerations of the unit being tested with the synchronous loading are plotted in Figure 17. The simulation test process involved the following six steps. Figure 16 shows a photograph of the built simulation test system. Mechanically, the impact and flight 2-D dynamic AE was simulated using an air cannon propulsion, a motor relay, and a vector axial angle tracking. The system has been applied to the performance test of aerospace components. The load time and the amplitude of the 2-D acceleration were set according to the performance index of the unit being tested. The axial acceleration of the impact overload was set as 20 g and the radial acceleration was zero g. The axial acceleration of the flight overload was maintained at 20 g, while the radial acceleration was increased from zero g to 40 g within 0.5 s.   The measured accelerations of the unit being tested with the synchronous loading are plotted in Figure 17. The simulation test process involved the following six steps. (1) The unit being tested was fixed on the vector turntable and the initial position of the unit being tested was vertical with the arm. The axial direction of the unit being tested pointed to the tangential direction of the arm, and the radial direction of the unit being tested pointed to the normal direction of the arm.

Verification Tests
(2) The electromagnetic clutch was disengaged, the asynchronous motor controlled based on ADRC accelerated to the set speed without the load, and the arm of the centrifuge remained stationary.
(3) The air cannon received the launch order and the simulation of the launch AE of the unit being tested began at 1.5 s. The air cannon pushed the centrifuge arm, starting the movement of the arm instantaneously, and then the impact overload was achieved. The radial acceleration of the unit being tested coincided with the tangential acceleration of the arm. Meanwhile the electromagnetic clutch was engaged, and the vector turntable was controlled by the servo motor based on the PVT mode. The orderly transition between the impact overload and the steady-state overload was achieved. The axial direction of the unit being tested was along the direction of the resultant acceleration of the arm. When the clutch was fully engaged, the arm was driven by the motor. The axial acceleration reached 20 g within 120 ms and the radial acceleration briefly reached a maximum of 10 g during the axial acceleration. The radial acceleration dropped sharply. As could be seen, the maximum of the radial acceleration of the unit being tested was 30 g without the synchronous loading in Figure 11, while it was 10 g with the synchronous loading in Figure 17. The maximum of the radial acceleration was greatly reduced, and this amount was safe and acceptable. It was verified that the acceleration could be effectively regulated by synchronous loading of 2-D dynamic acceleration.
(4) The unit being tested was rotated 90 • and the axial direction of the unit being tested pointed to the pivot of the arm. The axial direction of the unit being tested pointed to the normal direction of the arm, and the radial direction of the unit being tested pointed to the tangential direction of the arm. The system ran stably and maintained the launch state. The radial acceleration of the unit was reduced to zero g, and the axial acceleration was unchanged.
(5) The unit being tested received the flight order and the simulation of the flight AE began at 3 s. The centrifuge accelerated again, and the servo motor controlled the rotation angle of the unit being tested. By decomposing the resultant acceleration of the arm in two dimensions and maintaining the axial acceleration of the unit being tested constant, synchronous loading of 2-D dynamic acceleration was realized. The radial acceleration of the unit was increased from zero g to 41.7 g in 0.5 s, and the axial acceleration was 21.67 g.
(6) The system ran stably and retained the flight state. After the system became stabilized, the axial acceleration of the unit was 20 g, and the radial acceleration was 40 g, satisfying the launch and flight simulation test requirements for the 2-D dynamic AEs.

Conclusions
Performing ground assessment of key components is impossible due to lack of 2-D dynamic AEs on missile on-board equipment, aviation materials, acceleration sensors, weapon safety systems, or other mechanical systems. To solve this problem, we have investigated simulation methods of the 2-D dynamic AEs. By combining theoretical analysis, computer simulations, and simulation tests of real objects, we propose and develop an air-cannon-started hybrid centrifugal turntable simulation test method and the simulation test system. Simulation results reveal that the proposed test method achieves the synchronous loading of the 2-D dynamic AE and is able to effectively simulate impact and flight in extreme AEs. The test principle and testbed can be potentially extended to a variety of extreme impact simulation tests in dynamic AEs.
In the future, we plan to improve and optimize the motion control and position tracking control strategies using the developed simulation testbed in dynamic AEs. To this end, it is necessary to examine new impact excitation methods to reduce energy losses, optimize the system structure, and increase the test accuracy to meet the requirements of higher and more stable simulation test performance.