Research on Control Method of the Power System of Stepping-Type Anchoring Equipment

To improve the roadway adaptability and control accuracy of anchoring equipment, a stepping anchoring device was designed. A permanent-magnet synchronous motor control and a harmonic suppression algorithm were integrated to optimize the dynamic control system of stepping-type anchoring equipment. The results of an experimental simulation and analysis showed that when the coefficient of coal rock hardness f = 5, 6, and 7, the pulsation coefficient of the hydraulic pump outlet pressure, hydraulic motor output speed, and pump-controlled hydraulic cylinder advance speed in the hydraulic circuit of a pump-controlled motor did not exceed 3% after the equipment based on sliding mode control (SMC) entered the steady state, while the maximum pulsation coefficient was only 32.5% of the PI control. Based on the SMC, the harmonic components of the permanent magnet synchronous motor in the power system were suppressed and compensated for. This enhanced the stiffness of the hydraulic system under motor drive. When the rock stiffness factor gradually changed from f = 5 to f = 8 and increased suddenly from f = 5 to f = 6, the pressure overshoot at the outlet of the hydraulic pump of the pump-controlled motor system was reduced from 11.19% to 7.97% and from 61.19% to 52.88%, respectively, compared with that before the optimization. It was thereby proven that SMC based on harmonic suppression can effectively reduce the system pulsation caused by the multi-factor coupling of anchoring equipment and provide technical support for the optimal control of the power system of stepping-type anchoring equipment.


Introduction
Anchoring operation is an important task in coal mining. The high rate of occurrence of geological hazards in underground coal seams and strong contingency make this operation very challenging. According to statistics [1,2], roadway support operations account for more than 60% of the total drilling time. Therefore, the improvement of the reliability and work efficiency of anchor support operations will directly enhance the speed and efficiency of roadway excavation.
Stable and efficient electro-hydraulic power system is the foundation and guarantee for improving the operation efficiency of mining equipment [3][4][5]. At present, there are many kinds of electro-hydraulic control methods [6][7][8], mainly including valve control mode based on Intelligent PI control [9,10], intelligent PID control [11,12], SMC [13][14][15][16], and pump control mode based on variable speed driver technology [17,18]. Among the anchoring equipment [19][20][21], machine-mounted rotary hydraulic roof bolters are the most widely used. The electro-hydraulic system in such equipment often employs the Ziegler-Nichols algorithm to tune the PID parameters. Although this control method is simple and easy to implement, the parameter adjustment is more complicated and the control accuracy of hydraulic systems with low damping, strong nonlinearity, and time-varying parameters is not high [22,23]. Ahn et al. [24] designed a set of adaptive robust force controllers and used the gradient descent method to adjust the controller parameters to improve the stability of the hydraulic system under time-varying parameters and large external load fluctuations. However, this system has disadvantages, such as a poor lowspeed stability, limited speed adjustment range, slow response, and high requirements for oil quality. Wang [25] proposed a hydraulic system control strategy for a roof bolter based on a variable-frequency pump-controlled hydraulic drive system and PID parameter tuning with an intelligent fusion optimization algorithm. Simulations and experiments showed that under different working conditions, the roof bolter can apply intelligent fusion optimization algorithms to adjust the working parameters of the equipment automatically according to the changes in the peak oil pressure of the buffer chamber of the variable stroke mechanism. This ensures that the roof bolter performs well under reasonable parameter-matching working conditions and that the equipment functions at maximum power. Hoang [26] analyzed the influence of system parameter uncertainty and external interference on the control of the system, and put forward a new control strategy for the control of hybrid fluid power equipment based on the fundamentals of the fast terminal sliding mode, which ensured that the tracking error converges to the origin in a limited time.
In the foundation of the analysis and research above, a variety of intelligent algorithms were applied to optimize and tune the PID control parameters in the valve control system or pump-control system of the anchoring equipment to improve the control accuracy and stability of the equipment. Since the pump-control hydraulic system is driven by a frequency-conversion permanent magnet synchronous motor (PMSM), the controller, PMSM, hydraulic pump, etc., exhibit multivariable and nonlinear strong coupling effects [27][28][29], and the system pulses tend to be more complex.
In the present work, stepping-type anchoring equipment was designed. The overall coupling mechanism of the converter-PMSM, pump-controlled motor, and pumpcontrolled hydraulic cylinder was analyzed and a mathematical model of the drive control system of the walking anchoring equipment was constructed. Furthermore, a coordinated control algorithm for PMSM harmonic suppression compensation was proposed to actively suppress the pulsation of stepping-type anchoring equipment.

Stepping-Type Anchoring Equipment and Power System Composition
The stepping-type anchoring equipment designed to improve the roadway adaptability of the anchoring equipment is shown in Figure 1. The equipment mainly comprises a mechanical system and a power system. The mechanical system comprises subsystems such as advance support system, anchor net-laying system, walking system, and a roof bolter system. The walking system uses hydraulic cylinders for walking forward and is used for the mobile transportation of anchoring equipment; the anchor net-laying system is used to complete the net-laying work before anchoring the fully mechanized excavation face. The power system comprises a hydraulic system and an electric and control system. The hydraulic system mainly includes a quantitative hydraulic pump, an electromagnetic reversing valve, a supporting column hydraulic cylinder, a quantitative motor and pushing cylinder, and other components. The electric and control system comprises a PMSM, a controller, and its controller. The working process of the stepping-type anchoring equipment is as follows: after the equipment reaches the designated working position, the hydraulic cylinder of the walking system is supported on the ground; meanwhile, the hydraulic cylinder of the advance support system is supported on the ground, and the anchor net-laying system is extended to spread the anchor net to the roof and sidewalls of the roadway. Finally, the anchoring system adjusts the posture of the bolter drill so that it can complete anchoring operations at different positions.
To improve the ability of the stepping-type anchoring equipment to work in coal and rock roadways with different levels of hardness, the power system adopts a frequencyconversion electric motor to drive double pumps. Since the roof bolter system in this equipment is most affected by the hardness of the coal and rock, the power system of To improve the ability of the stepping-type anchoring equipment to work in coal and rock roadways with different levels of hardness, the power system adopts a frequencyconversion electric motor to drive double pumps. Since the roof bolter system in this equipment is most affected by the hardness of the coal and rock, the power system of the anchoring equipment was simplified, considering the roof bolter as the main research object. The power system and control unit are shown in Figure 2.

Mathematical Model of Stepping-Type Anchoring Equipment Power System
The optimal cutting speed and thrust speed of the roof bolter during the anchoring process of stepping-type anchoring equipment vary greatly with the hardness of the rock.   To improve the ability of the stepping-type anchoring equipment to work in coal a rock roadways with different levels of hardness, the power system adopts a frequen conversion electric motor to drive double pumps. Since the roof bolter system in t equipment is most affected by the hardness of the coal and rock, the power system of anchoring equipment was simplified, considering the roof bolter as the main research ject. The power system and control unit are shown in Figure 2.

Mathematical Model of Stepping-Type Anchoring Equipment Power System
The optimal cutting speed and thrust speed of the roof bolter during the anchor process of stepping-type anchoring equipment vary greatly with the hardness of the ro In addition, there exists a strong coupling between the transient mechanical characterist

Mathematical Model of Stepping-Type Anchoring Equipment Power System
The optimal cutting speed and thrust speed of the roof bolter during the anchoring process of stepping-type anchoring equipment vary greatly with the hardness of the rock. In addition, there exists a strong coupling between the transient mechanical characteristics of the PMSM and the hydraulic load; any change in a parameter of either will affect the working status of the entire system. Hence, it is necessary to establish a unified control equation under the coupling action of multiple factors, such as those of thecontroller-PMSM-hydraulic system. The use of such a unified control equation is essential to provide a theoretical basis for optimizing the performance of the pump-controlled hydraulic system.

Mathematical Model of Permanent Magnet Synchronous Motor
To improve the control accuracy of the power system of stepping-type anchoring equipment, considering the distortion magnetic field of the permanent magnet-synchronous electric motor, the dead time of the inverter, and the tube voltage drop, a large number of harmonic components are introduced into the three-phase current of the main circuit [30,31]. Therefore, a mathematical model of PMSM with motor distortion magnetic field and a harmonic current was built [32,33] as follows: Here, ω e = p n ω p where u d and u q are the two-phase rotating coordinate system d-q axis voltages (V); R is stator phase resistance (Ω); i d and i q are two-phase rotating coordinate system d-q axis currents (A); ψ d and ψ q are Flux linkage in the d-q rotating coordinate system (Wb); ω e is the electrical angular velocity (rad/s); L d and L q are two-phase rotating coordinate system d-q axis inductances (H); ψ f is the rotor flux linkage (Wb); θ e is the rotor electrical angle (rad); ω e is the electrical angular velocity (rad/s); ω p is the mechanical angular velocity (rad/s); p n is the number of pole pairs; J T is the moment of inertia of the motor (kg·m 2 ); B T is the damping coefficient of the motor shaft (N·m·s/rad); T e is the electromagnetic torque (N·m), and T L is the PMSM load torque (N·m).

Mathematical Model of Hydraulic System and Its Load
To facilitate the control of the rotary torque and feed rate of the roof bolter and improve the stability of the system, a "one machine with two pumps" system was adopted as the power source. A theoretical model of the hydraulic motor, hydraulic cylinder and load was established, and the influence of multiple factors on the stability of the system was explored.
where T pm/pg is the input torque of the pump-controlled motor system and pump-controlled hydraulic cylinder system (N·m); D pm/pg is the pump displacement of the pump-controlled motor system and pump-controlled hydraulic cylinder system (m 3 /rad); p pm/pg is the outlet pressure of the pump-controlled motor system and pump-controlled hydraulic cylinder system (MPa); η pm/pg is the mechanical efficiency of the pump-controlled motor system and pump-controlled hydraulic cylinder system; Q pm/pg is the pump-controlled motor system and pump-controlled hydraulic output flow of the cylinder system pump (m 3 /s); ω pm/pg is the pump speed of the pump-controlled motor system and pump-controlled hydraulic cylinder system (rad/s); λ p is the leakage coefficient of the pump (m 3 ); Q m is the input flow of the hydraulic motor (m 3 /s) V pm is the volume of the enclosed cavity between the pump and the hydraulic motor (m 3 ); V pm is the volume of the enclosed cavity between the pump and the hydraulic cylinder (m 3 ); β e is the elastic modulus of the oil volume in the enclosed cavity (N/m 2 ); D m is the displacement of the hydraulic motor (m 3 /rad); ω m is the speed of the hydraulic motor (rad/s); C pm/pg is the leakage coefficient of the pump-controlled motor system and pump-controlled hydraulic cylinder system; C m is the leakage coefficient of the hydraulic motor; C g is the leakage coefficient of the hydraulic cylinder; J m is the equivalent moment of inertia of the hydraulic motor output shaft (kg·m 2 ); η mm is the mechanical efficiency of the hydraulic motor; B m is the viscous damping coefficient of the hydraulic motor output shaft (N·m·s/rad); T Lm is the hydraulic motor load Torque (N·m); ν g is the feed speed of the hydraulic cylinder (m/s); A is the effective area of the hydraulic cylinder (m 2 ); m d is the equivalent mass of the hydraulic cylinder (kg); B g is the damping coefficient of the hydraulic cylinder (N·m·s/rad), and F g the is hydraulic cylinder load (N).
According to the drilling characteristics and rock breaking mechanism of hydraulic roof bolter, the drilling energy per unit volume of rock (i.e., the rock crushing work ratio), E/V, and the required drilling torque and propulsion for different rock hardness conditions can be expressed as follows [34,35]: where E 1 is the work done by the roof bolter to push the hydraulic cylinder within time ∆t (J); E 2 is the work done by the roof bolter rotating the motor within time ∆t (J); V is the rock-breaking volume of the roof bolter within time ∆t (m 3 ); D z is the diameter of the borehole (m); A z is the cross-sectional area of the borehole (m 2 ), and A is the active area of the piston rod of the hydraulic propulsion cylinder (m 2 ).
where c is the cutting speed constant; f is the rock's hardness coefficient; D is the outside diameter of the drill (mm); D 0 is the inner diameter of the drill (mm); "2.5~4" is the range of correction coefficients of the formula; σ is the compressive strength of the rock (MPa), and ∆ is the limit of drilling (mm).
Equations (3) and (4) show that the values of the rock crushing work ratio E/V required when the hydraulic roof bolter drills into different rocks differ from the best hydraulic motor speed for drilling different rocks. This leads to corresponding differences in the pressure and flow of the hydraulic cylinder and hydraulic motor circuit, especially under working conditions wherein the hardness coefficients of soft and hard rocks differ greatly; these differences will have a greater impact on the stable operation of the motor. Therefore, it is necessary to study the control method used for the frequency-conversion hydraulic roof bolter.

Research on the Control Method of Power System of Stepping-Type Anchoring Equipment
To improve the working ability of the stepping-type anchoring equipment in complex coal and rock roadways where there are large differences in the rock layer hardness, its power system adopts frequency-conversion vector control to ensure that there is a good matching relationship between the torque, speed, and rock hardness of the roof bolt.

PI Control of Pump-Controlled Hydraulic System
The variable frequency vector control method can independently control the motor current to realize the decoupling of the magnetic field and the torque and can effectively improve the mechanical characteristics of PMSM. Hence, the variable frequency control system of the stepping-type anchoring equipment adopts PI as the speed control loop for characteristic research. The dynamic system based on the control system of PI (CSPI) is shown in Figure 3.

Research on the Control Method of Power System of Stepping-Type Anchoring Equipment
To improve the working ability of the stepping-type anchoring equipment in com plex coal and rock roadways where there are large differences in the rock layer hardnes its power system adopts frequency-conversion vector control to ensure that there is a goo matching relationship between the torque, speed, and rock hardness of the roof bolt.

PI Control of Pump-Controlled Hydraulic System
The variable frequency vector control method can independently control the moto current to realize the decoupling of the magnetic field and the torque and can effective improve the mechanical characteristics of PMSM. Hence, the variable frequency contr system of the stepping-type anchoring equipment adopts PI as the speed control loop fo characteristic research. The dynamic system based on the control system of PI (CSPI) shown in Figure 3.

Sliding Mode Control of Pump-Controlled Hydraulic System
Linear controls such as PI control are easily affected by internal parameters and load fluctuations. In contrast, sliding mode control (SMC) is a discontinuous control strategy with time-varying switching characteristics. This process is unaffected by system parameters and external disturbances. The SMC of PMSM is described by [32]: where ω ref is reference speed of motor (rad/s), n ref = 30ω ref /π; i q * is the reference current of the q-axis (rad/s), and in the formula of the paper, ()' is the derivative of ().
In the paper, SMC is adopted as the system control method and the stability of vector control based on SMC is the basis of the stability of this control system. Therefore, the Lyapunov function is supplemented for stability analysis of the control system to demonstrate the validity and reliability of the content studied in this paper. The proof process is as follows [32]: The sliding mode surface function is defined as: Lyapunov function is defined as: Pursuant to Lyapunov's stability decision theorem, the requisite conditions for system stability are t = 0, V' < 0, i.e., ss' < 0: The paper employs the exponential convergence law method, which yields: Hence, this control system is asymptotically stable.
The dynamic system based on the control system of SMC (CSSMC) is shown in Figure 4.

Experimental Simulation and Analysis
To study the change law of the dynamic system during the anchoring operation of the stepping-type anchoring equipment and the control effect of the dynamic system based on the CSPI and CSSMC, MATLAB and AMESim were used for joint simulation. The simulation parameters are listed in Table 1.

Experimental Simulation and Analysis
To study the change law of the dynamic system during the anchoring operation of the stepping-type anchoring equipment and the control effect of the dynamic system based on the CSPI and CSSMC, MATLAB and AMESim were used for joint simulation. The simulation parameters are listed in Table 1.
The simulation conditions were set as: operating time t = 0-6 s, roadway rock hardness f = 5; t = 6-9 s, rock hardness f gradually increased from 5 to 6 and 8, respectively. According to Formula (4), the input curve of the simulation system, based on the expected torque and the expected speed of the hydraulic motor, and the expected feed force of the hydraulic cylinder is obtained, as shown in Figure 5. The experimental simulation results are shown in Figures 6-8. Table 1. Basic parameters of PMSM and the hydraulic system [27,30].

Symbol
Physical Significance Value The simulation conditions were set as: operating time t = 0-6 s, roadway rock hard ness f = 5; t = 6-9 s, rock hardness f gradually increased from 5 to 6 and 8, respectively According to Formula (4), the input curve of the simulation system, based on the expected torque and the expected speed of the hydraulic motor, and the expected feed force of th hydraulic cylinder is obtained, as shown in Figure 5. The experimental simulation result are shown in Figures 6-8.
where xmax, xmin, x are the maximum, minimum, and average values of the variables evaluated for pulsation.     8 show that, after applying a constant drive load (rock hardness coefficie f = 5), the hydraulic pump outlet pressure overshoot by the CSSMC is lower than by t CSPI. After the pressure, motor speed, and hydraulic cylinder propulsion speed enter steady state, the pulsation coefficients are reduced from 6.46%, 1.12%, and 18.18% (by t CSPI) to 2.10%, 0.69%, and 0.72% (by the CSSMC). This finding suggests that the propose CSSMC can improve the dynamic quality of the system and that the physical quantity the system output is smoother compared with the traditional CSPI.
To further study the performance of the hydraulic system based on the CSSMC, t hydraulic cylinder and hydraulic motor load of the anchoring equipment were used provide a step signal input to simulate the operating conditions of the anchoring equi ment under conditions of larger changes in the hardness of the rock layers: for t = 0-6 the roadway rock hardness is f = 5; for t = 6-11 s, the roadway rock hardness is f = 6-7. T change curves for outlet pressure of the hydraulic pump, the output speed of the hydra lic motor, and the forward speed of the hydraulic cylinder in the hydraulic circuit of th pump control motor are shown in Figures 9-11  The speed pulsation coefficient of the hydraulic motor is defined as δ nm and the velocity pulsation coefficient of hydraulic cylinder is defined as δ ng : The pulsation coefficient of the outlet pressure pulsation is defined as follows δ p a : where x max , x min , x are the maximum, minimum, and average values of the variables evaluated for pulsation. Figures 6-8 show that, after applying a constant drive load (rock hardness coefficient f = 5), the hydraulic pump outlet pressure overshoot by the CSSMC is lower than by the CSPI. After the pressure, motor speed, and hydraulic cylinder propulsion speed enter a steady state, the pulsation coefficients are reduced from 6.46%, 1.12%, and 18.18% (by the CSPI) to 2.10%, 0.69%, and 0.72% (by the CSSMC). This finding suggests that the proposed CSSMC can improve the dynamic quality of the system and that the physical quantity of the system output is smoother compared with the traditional CSPI.
To further study the performance of the hydraulic system based on the CSSMC, the hydraulic cylinder and hydraulic motor load of the anchoring equipment were used to provide a step signal input to simulate the operating conditions of the anchoring equipment under conditions of larger changes in the hardness of the rock layers: for t = 0-6 s, the roadway rock hardness is f = 5; for t = 6-11 s, the roadway rock hardness is f = 6-7. The change curves for outlet pressure of the hydraulic pump, the output speed of the hydraulic motor, and the forward speed of the hydraulic cylinder in the hydraulic circuit of the pump control motor are shown in Figures 9-11, respectively.      When the coefficient of rock hardness changes abruptly from f = 5 to f = 6, the hyd lic pump outlet pressure of the pump-controlled motor system based on the CSPI flu ates the most, with an overshoot of 67.87%, while the pulsation coefficients of press motor speed, and hydraulic cylinder propulsion speed after entering the steady state come 5.87%, 1.32%, and 17.02%, respectively. By the CSSMC, the hydraulic pump o pressure, the speed of hydraulic motor, and the speed of hydraulic cylinder are 0. 0.18%, and 0.31%, respectively (see . When the rock hardness coefficie abruptly changed from f = 5 to f = 7, the parameters of the hydraulic control system ch in a similar pattern, which will not be repeated here. To further compare and analyz pulsation of each physical quantity after the system enters the steady state in Figure  11, the statistical results are shown in Table 2.

Performance Improvement Analysis of the Power Control System of Stepping-T Anchoring Equipment
From Figures 6-11, it can be seen that the dynamic quality of the system can be proved based on the CSSMC compared with the CSPI. However, it was also found When the coefficient of rock hardness changes abruptly from f = 5 to f = 6, the hydraulic pump outlet pressure of the pump-controlled motor system based on the CSPI fluctuates the most, with an overshoot of 67.87%, while the pulsation coefficients of pressure, motor speed, and hydraulic cylinder propulsion speed after entering the steady state become 5.87%, 1.32%, and 17.02%, respectively. By the CSSMC, the hydraulic pump outlet pressure, the speed of hydraulic motor, and the speed of hydraulic cylinder are 0.39%, 0.18%, and 0.31%, respectively (see . When the rock hardness coefficient is abruptly changed from f = 5 to f = 7, the parameters of the hydraulic control system change in a similar pattern, which will not be repeated here. To further compare and analyze the pulsation of each physical quantity after the system enters the steady state in Figures 6-11, the statistical results are shown in Table 2.

Performance Improvement Analysis of the Power Control System of Stepping-Type Anchoring Equipment
From Figures 6-11, it can be seen that the dynamic quality of the system can be improved based on the CSSMC compared with the CSPI. However, it was also found that under the condition of a sudden load change, the pulsation overshoot of the system is large, meaning that it is necessary to further optimize its control strategy to improve the kinetic energy stiffness of the hydraulic control power system.

Harmonic Suppression Compensation Coordinated Control of Permanent Magnet Synchronous Motor
According to the working characteristics of PMSM, the mechanical characteristics of PMSM periodically pulsated change due to the harmonic current of the converter, the tube resistance voltage drop, etc., while the pulsating component mainly exists in the form of the 5th and 7th harmonics [36,37]. Therefore, to realize the current control of the d-axis and q-axis of the 5th and 7th harmonics in the system, a synchronous rotating coordinate axis system corresponding to the 5th and 7th harmonics was established. The coordinate transformation matrix from the d-q coordinate system to 5d-q and 7d-q is thus as follows [36][37][38]: By substituting Equation (9) into Equation (1), we obtain the following relation: From Equation (10), we can see that in the 5d-q coordinate system, only the 5th-order current is the DC component, while the rest are AC components. A low-pass filter can be used to extract the 5th harmonic component. The 7th harmonic component can also be extracted similarly [37]. Meanwhile, according to the tracking characteristics of the PI controller without static error and, given i d5 = 0, i q5 = 0, i d7 = 0, and i q7 = 0, the current harmonics pass through the PI regulator, and then the coordinates of 5d-q and 7d-q are transformed into a d-q coordinate system through Equation (15). Then, the corresponding voltage compensation can be obtained [36,37]: To enhance the robustness of the control system, we propose a cooperative control algorithm that combines the SMC and harmonic suppression compensation. The principle of collaborative control system based on SMC and harmonic suppression (CCSSMC-HS) is shown in Figure 12.

Control Method Improvement Results and Analysis
In order to compare and analyze the performance of the control strategy based o CCSSMC-HS, the experimental simulation parameters and process are kept the sa those used and completed in the previous section. The results are shown in Figure  18.

Control Method Improvement Results and Analysis
In order to compare and analyze the performance of the control strategy based on the CCSSMC-HS, the experimental simulation parameters and process are kept the same as those used and completed in the previous section. The results are shown in Figures 13-18.

Control Method Improvement Results and Analysis
In order to compare and analyze the performance of the control strategy based on CCSSMC-HS, the experimental simulation parameters and process are kept the sam those used and completed in the previous section. The results are shown in Figures 18.      15 show that the change curve of each physical quantity of the hydra system based on the CSSMC-HS coincides with the curve based on the CSSMC. The r tive deviation between the time when they enter the steady state and the pulsation co cient under the steady state are both less than 2%. While the pump controlling the draulic system was based on the CSSMC (the rock hardness factor gradually chan from f = 5 to f = 8), the hydraulic pump outlet pressure overshoot was reduced from 11.1 to 7.97% compared with that before the optimization, which reduced the impact of load change on the system stability.   Figures 13-15 show that the change curve of each physical quantity of the hydra system based on the CSSMC-HS coincides with the curve based on the CSSMC. The r tive deviation between the time when they enter the steady state and the pulsation co cient under the steady state are both less than 2%. While the pump controlling the draulic system was based on the CSSMC (the rock hardness factor gradually chan from f = 5 to f = 8), the hydraulic pump outlet pressure overshoot was reduced from 11. to 7.97% compared with that before the optimization, which reduced the impact of load change on the system stability.  Figures 13-15 show that the change curve of each physical quantity of the hydraulic system based on the CSSMC-HS coincides with the curve based on the CSSMC. The relative deviation between the time when they enter the steady state and the pulsation coefficient under the steady state are both less than 2%. While the pump controlling the hydraulic system was based on the CSSMC (the rock hardness factor gradually changed from f = 5 to f = 8), the hydraulic pump outlet pressure overshoot was reduced from 11.19% to 7.97% compared with that before the optimization, which reduced the impact of the load change on the system stability.       The change law of each physical quantity of the hydraulic system based on the CSSMC-HS in Figures 16-18 is the same as that shown in Figures 13-15. When the rock stiffness factor increased suddenly from f = 5 to f = 6, the pressure overshoot at the outlet of the hydraulic pump of the pump-controlled motor system was reduced from 61.19% to 52.88%, compared with that before the optimization. The reason for this is that as the PMSM operates at a high speed, there is a sudden change in the load torque and the influence of the harmonic current in the system is significantly enhanced, which in turn decreases the rigidity of the hydraulic system driven by the motor. For the CSSMC-HS, the harmonic components are suppressed and compensated for; hence, the rigidity of the hydraulic system driven by the motor is enhanced. That is, the proposed vector control strategy based on the CSSMC-HS can effectively reduce the impact of load disturbance on the system's performance. The change law of each physical quantity of the hydraulic system based on CSSMC-HS in Figures 16-18 is the same as that shown in Figures 13-15. When the r stiffness factor increased suddenly from f = 5 to f = 6, the pressure overshoot at the ou of the hydraulic pump of the pump-controlled motor system was reduced from 61.19% 52.88%, compared with that before the optimization. The reason for this is that as PMSM operates at a high speed, there is a sudden change in the load torque and the fluence of the harmonic current in the system is significantly enhanced, which in t decreases the rigidity of the hydraulic system driven by the motor. For the CSSMCthe harmonic components are suppressed and compensated for; hence, the rigidity of hydraulic system driven by the motor is enhanced. That is, the proposed vector con strategy based on the CSSMC-HS can effectively reduce the impact of load disturbance the system's performance.

Conclusions
(1) Following the pump control mechanism of PMSM and the system theory of interaction between anchoring equipment and coal rocks, a mathematical model of the power control system of step-type anchoring equipment was established, and the influence of coal rock hardness factor f on the stability of the power system of the anchoring equipment was analyzed. (2) A control strategy applicable to the drive system of walking anchoring equipment based on the CSSMC was also proposed. Compared with the CSPI, the dynamic characteristics of the hydraulic pump outlet pressure, hydraulic motor output speed, and pump-controlled hydraulic cylinder advance speed in the pump-controlled motor hydraulic circuit of this control method were greatly improved, while the maximum pulsation coefficient of each variable after the system entered the steady state was less than 3%. (3) Furthermore, the CSSMC was optimized by suppressing and compensating the harmonic components of PMSM in the power system. For the CSSMC-HS, when the coefficient of coal and rock hardness in action gradually changed from f = 5 to f = 8 and abruptly changed from f = 5 to f = 6, the hydraulic pump outlet pressure overshoot fell from 11.19% to 7.97% and from 61.19% to 52.88%, respectively, compared to the levels prior to optimization. This optimized control strategy was therefore able to mitigate the impact of load disturbance on the system's performance.
Author Contributions: Conceptualization and methodology, G.S. and Y.G.; theoretical analysis, simulation and writing-original draft preparation, G.S., P.W. and D.Z.; funding acquisition, G.S., Y.G., P.W. and G.C.; writing-review and editing, D.Z. All authors have read and agreed to the published version of the manuscript.