Next Article in Journal
An Improved Pedotransfer Function for Soil Hydrological Properties in New Zealand
Previous Article in Journal
The DIKWP (Data, Information, Knowledge, Wisdom, Purpose) Revolution: A New Horizon in Medical Dispute Resolution
Previous Article in Special Issue
Model-Free Adaptive Nonsingular Fast Integral Terminal Sliding Mode Control for Wastewater Treatment Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trajectory Tracking Algorithm Study of Coal Mine Water Detector Drilling Bar Installation

1
School of Mechanical Engineering, Inner Mongolia University of Technology, Hohhot 010051, China
2
School of Mechanical and Electrical Engineering, Harbin Engineering University, Harbin 150001, China
3
School of Mechanical Engineering, Jiamusi University, Jiamusi 154007, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(10), 3996; https://doi.org/10.3390/app14103996
Submission received: 7 April 2024 / Revised: 4 May 2024 / Accepted: 6 May 2024 / Published: 8 May 2024
(This article belongs to the Special Issue Research and Application of Intelligent Control Algorithm)

Abstract

:
Mechanical water detection is recognized as the most reliable and safe production technology for coal mines, mainly for the detection of water hazards in pre-mining operations. Intelligent water detectors are currently the main research direction in mechanical water detection, and the automatic installation of drilling bars is the key to achieving intelligent water detection. Improving the connection accuracy in the process of installing drilling bars is an important research topic for the improvement of control links. To improve the connection accuracy of the drilling bars at the time of supplying material, we used the modified Denavit–Hartenberg method to analyze the motion gestures of the supplied material device and the Lagrange equation to establish a dynamic analysis model. We aimed at better control precision by improving the sliding mode control algorithm and at increasing the convergence rate of tracking errors with a sliding controller based on an exponential approximation law and using saturated functions instead of the symbol functions in the reaching law to weaken the vibration in the control process. We then used particle swarm optimization (PSO) to find the optimum combination parameters of the sliding mode controllers and test the performance of the sliding mode controllers before and after PSO with MATLAB/Simulink. The results showed that the optimized controller has a strong resistance to parameter fluctuations, and the system responds quickly, achieves a good performance, and improves the convergence rate of tracking errors.

1. Introduction

Detecting water is an important factor in safe production in coal mines. Mechanical water detection is recognized as the most reliable water detection technology, and its intelligent transformation is an inevitable development, in accordance with the development of science and technology, as well as the need for industrial development. At the beginning of the 21st century, Epiroc, Sandvik, Boart Longyear, and other companies entered a period of rapid development in drilling technology, developing the Smart ROC-Series intelligent drill with automated drilling holes, drilling bar loading, drilling, and other functions, the Pit Viper-Series automatic drill with remote control and automatic drilling bar switching functions, and the RDS Robot system drilling bar installation module, as shown in Figure 1. The method of using human-assisted installation and dismantling drilling bars means that the worker’s labor intensity is high, and people and machines working together can result in safety accidents [1]. An example of an artificial drilling operations scenario is shown in Figure 2. In 2019, China’s Coal Mine Safety Supervisory Authority published the “Coal Mine Robot Focus R&D Catalogue”, which requires that the water detector can be automatically installed and that drilling bar functions can be dismantled. In 2020, “the Guidelines on Accelerating the Intelligent Development of Coal Mines” clearly proposed raising the level of intelligence of coal mines to basic intelligence levels by 2025. The “14th Five-Year’ High Quality Development Guidelines for the Coal Industry”, issued in 2021, clearly proposed to significantly reduce the number of underground operators. China University of Mining and Technology and Harbin Institute of Technology jointly developed a drill robot that can automatically transport, install, and dismantle drilling bars. Chongqing Research Institute achieved fully automatic drilling of a drill Z through research and testing in the areas of automatic installation and dismantling drilling bars and automated drilling.
In the intelligent operation of mechanical water detection, how to ensure the connection accuracy of drilling bars when automatically completing the assembly tasks is the core technical issue for the development of water detection machinery [2]. The study of control algorithms for supplied material devices to optimize the trajectory tracing controller is a concrete safeguard measure to improve the connection accuracy [3,4]. Based on the theory of inverse kinematics, the application of the trajectory tracking control algorithm converts the Cartesian trajectory tracking control problem to a joint trajectory tracking control problem to study the accuracy of the control methods of the supplied material device [5].
Interference factors in the underground environment are numerous and uncertain, such as external disturbances of the controlled system, limiting the increase in tracking control accuracy [6,7]. At the same time, the supplied material device, as a nonlinear strong coupling system [8], has higher requirements for the control algorithm. The sliding mold control is essentially a special type of nonlinear control, in the design of which the sliding modules are unrelated to the parameters and disturbances of the controlled object, which can lead to effects of uncertainty and nonlinearity in most practical systems [9,10]. The good characteristics of the dynamic response make this control algorithm widely used in nonlinear control systems [11]. Therefore, this paper improves the control of the supplied material device by improving the sliding mold controller. In response to the need for the automatic injection of intelligent water detectors in coal mines, the supplied four-degrees-of-freedom material device which is the object of study here demonstrates a sliding-model control strategy based on PSO to optimize the trajectory tracking controller for the purpose of improving its joint trajectory tracking accuracy. First, we established a dynamic model of the supplied material device. Then, we designed a sliding mold controller based on the exponential approximation law and used the saturation function to replace the symbol function in the reaching law to weaken the vibration. Finally, using PSO to optimize the above controller parameters, the controller’s optimal combination of parameters is obtained. The simulation results showed that the controller optimized by PSO was better.

2. Motion Attitude Analysis Based on the MDH Method

In Figure 3 below, 3 is the base seats, 4, 5, 7, 8 are joints, 6 is the functional connectors, and 9 is the mechanical tongs, composing the supplied material device. The drilling bars are placed on the storage rack, and the supplied material device transports the new drilling bars from the storage rack to the planned installation location. A precise connection of the new drilling bars with the front drilling bars’ position can be achieved by adjusting the positional attitude of the mechanical tongs. A diagram of the relative positions is shown in Figure 3.
We preliminarily established the dynamic analysis model of the supplied material device by precisely defining its positional attitude. The supplied material device, formed by a series of functional connectors through the rotary joint, is a type of open-motion chain, and its characteristics comply with the MDH method [12]. Based on the above methods, the establishment of kinematic structural model was achieved, as shown in Figure 4. Based on the location of the rotary joint, four rectangular coordinates were established, with serial numbers of 0, 1, 2, and 3. The base coordinate corresponds to the 0 coordinate. The origin coordinate is at the bottom seat of the supplied material device, with the remaining three coordinates being the second and three rotary joints and the arm joints, respectively. The MDH coordinate model is shown in Figure 5.
Based on the coordinate model in Figure 4, the D-H parameters of the supplied material device are listed in Table 1. The meaning of the parameters is as follows: a is the length of the connecting rod, representing the distance between joint zi and zi+1 along the x axis; α is the angle of rotation, representing the changes in the rotation angle of joint zi to zi+1 around the x axis; d is the offset distance, representing the distance between joint xi and xi+1 along the z axis; θ is the angle of the connecting rod, representing the changes in the rotation angle of joint xi to xi+1 around the z axis, where x1 and x0 coincide when θ1 is equal to 0° [13].
The angular displacement of each joint was designed using the forward kinematics of the robot [14]. At the same time, the positional attitude of the end coordinates of the hand were solved in relation to the base coordinate system. We resolved the positional attitude of the hand using known linear parameters and joint variables. The homogeneous transformation matrix of the adjacent connecting rod coordinate system, based on the rigid body posture attitude transformation theory, is as follows [15]:
  i i 1 T = R P 0 1 = cos θ i sin θ i 0 a i 1 sin θ i cos α i 1 cos θ i cos α i 1 sin α i 1 d i sin α i 1 sin θ i sin α i 1 cos θ i sin α i 1 cos α i 1 d i cos α i 1 0 0 0 1
By substituting the D-H parameter values in Table 1 into Formula (1) for calculation, the coordinate transformation matrix from joint1 to joint4 of the profiling mechanism can be obtained as follows:
  1 0 T = cos θ 1 sin θ 1 0 0 sin θ 1 cos θ 1 0 0 0 0 1 0 0 0 0 1
  2 1 T = cos θ 2 sin θ 2 0 0 0 0 1 0 sin θ 2 cos θ 2 0 0 0 0 0 1
  3 2 T = cos θ 3 sin θ 3 0 L 1 sin θ 3 cos θ 3 0 0 0 0 1 0 0 0 0 1
  4 3 T = cos θ 4 sin θ 4 0 L 2 sin θ 4 cos θ 4 0 0 0 0 1 0 0 0 0 1
The positional attitude transformation matrix of the mechanical tongs relative to the base coordinates, using Formulas (2)–(5), is as follows:
T 4 0 = i = 1 4   i i 1 T = R 4 0 P 4 0 0 0 0 1 = T 11 T 12 sin θ 1 L 1 cos θ 1 cos θ 2 sin θ 1 T 21 T 22 cos θ 1 L 2 sin θ 1 cos ( θ 2 + θ 3 ) T 31 T 32 0 L 2 sin ( θ 2 + θ 3 ) L 1 sin θ 2 0 0 0 1
In the formula,
T 11 = cos θ 4 cos θ 1 cos ( θ 2 + θ 3 ) sin θ 4 cos θ 1 sin ( θ 2 + θ 3 ) ;
T 12 = sin θ 4 cos θ 1 cos ( θ 2 + θ 3 ) cos θ 4 cos θ 1 sin ( θ 2 + θ 3 ) ;
T 21 = cos θ 4 sin θ 1 cos ( θ 2 + θ 3 ) sin θ 4 sin θ 1 sin ( θ 2 + θ 3 ) ;
T 22 = sin θ 4 sin θ 1 cos ( θ 2 + θ 3 ) cos θ 4 sin θ 1 sin ( θ 2 + θ 3 ) ;
T 31 = cos θ 4 sin ( θ 2 + θ 3 ) sin θ 4 cos ( θ 2 + θ 3 ) ;
T 32 = sin θ 4 sin ( θ 2 + θ 3 ) cos θ 4 cos ( θ 2 + θ 3 ) .
The MATLAB R2016b robotic toolbox can be used to construct a motion model of the supplied material device to obtain the spatial position of the terminal executor at its initial state, as shown in Figure 6.

3. Establishment of the Dynamic Model of the Supplied Material Device

Use the Lagrange equation to construct the dynamic analysis model by describing the joint control’s input force matrix and the relationship between the position, speed, and acceleration of the supplied material device. The specific process is as follows:
Based on the definition of the Lagrange function L, first establish the Lagrange function, defining the difference between the system’s kinetic energy K and potential energy P as follows [16]:
L = K P
Here, the coordinate representation method of K and P is not restricted. Defining the dynamic equation of the system, namely, the Lagrange equation, is achieved using the following formula [16]:
F i = d d t L q ˙ i L q i , i = 1 , 2 , n
In the formula, q i represents the coordinates of the system’s kinetic energy and potential energy, q ˙ i is the corresponding speed, n is the number of connecting rods, and Fi is the force matrix of the ith joint on the supplied material device. Finally, the formulation, taking into account external interference and unmodeled dynamics, is modeled as follows [17]:
M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) + d ( t ) = u ( t )
In the formula, q = q 1 , q 2 , , q n T represents the joint position vector of the supplied material device; q ˙ R n and q ¨ R n represent the angular velocity vector and angular acceleration vector, respectively; M ( q ) R n × n represents the symmetrical inertia matrix of 2 × 2; C ( q , q ˙ ) R n × n represents the Cosmic force and centrifugal matrix of 2 × 2; G ( q ) R n represents the gravitational vector of 2 × 1; d ( q ) R n represents the external interference and unmodeled dynamics of the supplied material device of 2 × 1; and u R n represents the joint-controlled force matrix vector of 2 × 1.
The expected angle qd and its derivative are bounded as a specific study premise, which is carried out by optimizing the control algorithm, achieving the positioning accuracy control of the supplied material device, while effectively shielding it from external interference.

4. Stability Design of the Controller System

4.1. Design of Controller

Given the characteristics of the upper part of the material device, as well as the sliding mold control algorithm, use the sliding mold control to optimize the trajectory tracking control of the supplied material device. The expected position of the joint angle of the supplied material device is qd, and the tracking error function is defined as follows:
e ( t ) = q d q
The design of the sliding surface is as follows:
s ( t ) = c e ( t ) + e ˙ ( t )
In the formula, c = d i a g c 1 ,   c 2 ,     c n ,   c i   >   0 . When t , the index of the tracking error of the joint corner trajectory converges; the greater the value of c is, the faster the convergence rate is. Improve the dynamic performance of the supplied material device system by optimizing the design of the diagonal matrix c. This is obtained from Formulas (10) and (11):
s ˙ ( t ) = c e ˙ ( t ) + e ¨ ( t ) = c e ˙ ( t ) + q ¨ d q ¨
Although the slide control improves the robustness of the system, at the same time, the temporal and spatial delays of the switching surface also lead to the control system vibrating. The vibration cannot be eliminated, there are only ways to weaken the vibration. In order to suppress the vibration, this paper uses the exponential approach law and replaces its symbol function with the saturation function; the expression of the reaching law is as follows:
s ˙ = ε s a t ( s ) k s , ε > 0 , k > 0
In the formula, the saturation function is s a t ( s ) = 1 , s > Δ k s , s Δ 1 , s < Δ .
Larger ε values can cause vibration, and this vibration can be weakened by adjusting the parameter ε of the reaching law to ensure the quality of the sliding module’s arrival process.

4.2. Proof of the Stability of the Supplied Material Device’s Control System

To prove the stability of the reaching law, select the Lyapunov function validation:
V = 1 2 s 2
It can be obtained from Formulas (13) and (14) that
s s ˙ = s ε s a t ( s ) k s 2 < 0
From Formula (15), it is obtained, based on saturation functions, that the exponential approach law satisfies Lyapunov’s theory of stability [18].
The system control rates obtained from Formulas (12) and (13) are as follows:
u ( t ) = C ( q , q ˙ ) q ˙ + G ( q ) + d + M [ c e ˙ ( t ) + q ¨ d + ε s a t ( s ) + k s ]

5. Parameter Optimization of the Controller Based on PSO

Since the sliding surface coefficient c and the parameters of the reaching law ε and k can affect the tracking effect of the supplied material device when designing the slider controller, it is necessary to optimize the design of the above three parameters. The efficient search capability of PSO helps find better solutions, and its global search feature can find the best solution in a larger range and avoid falling into the local best. Therefore, this paper uses particle swarm optimization, seeking the optimal combination of the above three parameters to improve the robustness of the control system. PSO is a method to search for the best result through a large number of parameters to optimize an optimal result, with structural simplicity and rapid convergence [19]. The optimal combination of parameters can be found through crowd search, and track tracking control of the supplied material device can be achieved based on the assumption that the target function J obtains the minimum value. The specific design steps are as follows [20]:

5.1. Determination of the Initial State of the Parameter

Suppose that an initial controller particle group PSO is randomly composed of n particles in a D-dimensional space, i.e., P S O = p 1 ,   p 2 ,     p n , and p i = ( c o n 1 p i , c o n 2 p i , c o n 3 p i ) represents the three parameters of the controller to be optimized. The particle iteration has the corresponding flight speed v p i and position x p i and needs to meet v max < v p i < v max , x max < x p i < x max .

5.2. Determination of Controller Parameters Based on PSO

During the iteration, the particle dynamically adjusts the speed and position through individual and global extremums to constantly approach the optimal target at a certain rate [21]. In the t + 1 iteration, the speed and position of the particle update are as follows:
v i d ( t + 1 ) = ω p s o v i d ( t ) + c 1 r 1 ( P i d ( t ) x i d ( t ) ) + c 2 r 2 ( G i d ( t ) x i d ( t ) ) x i d ( t + 1 ) = x i d ( t ) + v i d ( t + 1 )
In the formula, P i d t is the individual extremum of particle i based on its own flight experience, G i d t is the global extremum obtained from the group of particles’ flight experience, ωpso is the inertia weight factor, and c1, c2 is the acceleration constant.

5.3. Determine Fitness Function

Use the fitness function to evaluate the properties of the particles [22], using the control system as the object of study to find a set of optimal combinations of parameters for the sliding surface coefficient c, as well as the approximation law parameters ε and k, so that the tracking error of the system is minimized. Based on the above requirement, take the sum of the errors of the two joints as the function to find the optimal parameter for the controller; the objective functions are as follows:
J = e 1 ( t ) + e 2 ( t )

5.4. Algorithm Flow

In this algorithm, the input parameter is the particle group’s parameter and the acceleration constant of the sliding mold controller of the supplied material device, and the output result is the optimal parameter of the sliding mode controller. The specific processes are shown in Figure 7.
First, initialize the algorithm parameter and set the number of iterations, t, and number of sliding mode controller group particles, n. While randomly defining the initial position, flight speed, and acceleration constant of the particles, assume that the individual extremum of each particle is its current optimal position, and that the global extremum is the best particle position in the individual extremum. Calculate the fitness of each particle in the group using Formula (18); then, compare the particle fitness values to the individual extremum and the global extremum. If the current position of the particle is superior to the individual extremum, update the individual extremum, and if the particle swarm position is superior to the global extremum, update the global extremum. Obtain the particle in question’s up-to-date location and flight speed based on the updated individual extremum, the global extremum, and Formula (17). Finally, if it is found that the number of iterations exceeds the set maximum number, the iteration is stopped and the optimal solution is obtained, i.e., the optimal parameter of the sliding mode controller for the trajectory of the supplied material device; otherwise, Formula (18) is used to calculate the fitness of each particle in the group.
After the MATLAB/Simulink R2016b simulation iteration, the optimized results of PSO are obtained, as shown in Figure 8 and Figure 9.
As shown in Figure 8, after 15 iterations, the PSO has found the global extremum, at which point the target function achieves the minimum value. As shown in Figure 9, the red line represents the c value; the green line represents the k value; and the blue line represents the ε value, so the controller parameters have reached the optimum: c = 4.27, k = 1.5, ε = 2.1.

6. Analysis of the Trajectory Tracking Control of the Joint

Examples of two joints of a four-degrees-of-freedom serially supplied material device are given in this paper. Build a dynamic model of the supplied material device, then design and optimize the sliding mode controller, and ultimately optimize the controller parameters using PSO. Simulate and analyze the effect of the sliding mode controller before and after PSO on the trajectory tracking control accuracy of the supplied material device. The control structure of the supplied material system and the construction of the simulation module are shown in Figure 10 and Figure 11.
The matrix parameters in the dynamic model are as follows:
M ( q ) = P + m 2 l 2 2 + 2 m 2 l 1 l 2 cos q 2 m 2 l 2 2 + m 2 l 1 l 2 cos q 2 m 2 l 2 2 + m 2 l 1 l 2 cos q 2 m 2 l 2 2 ;
C ( q , q ˙ ) = m 2 l 1 l 2 q ˙ 2 sin q 2 m 2 l 1 l 2 ( q ˙ 1 + q ˙ 2 sin q 2 ) m 2 l 1 l 2 q ˙ 1 sin q 2 0 ;
G ( q ) = ( m 1 + m 2 ) g l 1 cos q 2 + m 2 g l 2 cos ( q 1 + q 2 ) m 2 g l 2 cos ( q 1 + q 2 )
In the formula, p = (m1 + m2)l12. The actual trajectories of the two connecting rods are q1 and q2, respectively. Create a dynamic model of the supplied material device in MATLAB/Simulink and simulate the process, with the simulation parameter set to m1 = m2 = 0.5 kg, and L1 = L2 = 0.8 m. They are the mass and length of the two connecting rods, respectively. The control effect of the sliding mode controller on the track tracking of the supplied material device under two different parameter determination methods can be determined.
As shown in Figure 12, the sliding mold controller with the manually determined parameter traces the expected trajectory of joint3 in 2.24 s. As shown in Figure 13, after optimization of the controller parameter by PSO, joint3 traces the expected trajectory in 1.60 s under the control of that controller. It can be seen by comparing Figure 12 and Figure 13 that the new controller model improves the control of joint3 after optimizing its parameters by PSO, resulting in a 36.6% reduction in the time taken to trace the expected trajectory.
As shown in Figure 14, the sliding mold controller with the manually determined parameter traces the expected trajectory of joint4 in 2.52 s. As shown in Figure 15, after optimization of the controller parameter by PSO, joint4 traces the expected trajectory in 2.03 s under the control of that controller. It can be seen by comparing Figure 14 and Figure 15 that the new model controller improves the control of joint4 after optimizing its parameters by PSO, resulting in a 19.4% reduction in the time taken to trace the expected trajectory.
As shown in Figure 16, the sliding mold controller with the manually determined parameter traces the expected speed of joint3 in 2.62 s. As shown in Figure 17, after optimization of the controller parameter by PSO, joint3 traces the expected speed in 2.19 s under the control of that controller. It can be seen by comparing Figure 16 and Figure 17 that the new model controller improves the control of joint3 after optimizing its parameters by PSO, resulting in a 16.4% reduction in the time taken to trace the expected speed.
As shown in Figure 18, the sliding mold controller with the manually determined parameter traces the expected speed of joint4 in 2.91 s. As shown in Figure 19, after optimization of the controller parameter by PSO, joint4 traces the expected speed in 2.03 s under the control of that controller. It can be seen by comparing Figure 18 and Figure 19 that the new model controller improves the control of joint4 after optimizing its parameters by PSO, resulting in a 30.2% reduction in the time taken to trace the expected speed.
As shown in Figure 20, the slider mold controller uses 2.8126 s to converge the track tracking error of the supplied material device to 0.0012 when the parameter is manually determined. As shown in Figure 21, after optimization of the controller parameter by PSO, the slider mold controller uses 2.0321 s to converge the track tracking error of the supplied material device to 0.0010 under the control of that controller. It can be seen by comparing Figure 20 and Figure 21 that the new model controller has better controls after optimizing its parameters by PSO, resulting in smaller eventual convergence errors and a shorter usage time.
In summary, the slide mold controller after the optimization of parameters with PSO has several advantages in terms of the position tracking control and speed tracking control of the joints, and it reduces the convergence time of tracking errors. Details of this are shown in Table 2.

7. Conclusions

In order to improve the connection accuracy of drilling bars in the automatic drilling process of intelligent water detectors, the trajectory tracking control of the drilling device was studied. First, the MDH method described the positional attitude of the supplied material device and was then used to construct the dynamic model based on the Lagrange equation. The sliding mold controller was designed on the basis of the dynamics model. The controller used the exponential approximation law and replaced the symbol function with the saturated function to weaken the inevitable vibration in the sliding mold control. Finally, using PSO to optimize the parameters of the controller, compared to the manually determined controller parameters, the combination of parameters obtained by PSO leads to sliding mold controller having a better effect on joints3 and -4 of the supplied material device, which can be summarized by the following points:
(1)
The time of the expected trajectory of joint3’s tracking decreased by approximately 36.6%. The time of the expected trajectory of joint4’s tracking decreased by approximately 19.4%. The above-detailed description of an optimized sliding mold controller increases the speed of the expected position on the tracking of the joints.
(2)
The expected speed of joint3’s tracking decreased by approximately 16.4%. The expected speed of joint4’s tracking decreased by approximately 30.2%. The above-detailed description of an optimized sliding mold controller decreases the expected speed of the tracking of the joints.
(3)
The sliding mold controller, before optimizing, used 2.8126 s to converge the overall tracking error of the supplied material device to 0.0012. Compared to the above effects, it took 2.0321 s to converge the overall tracking error of the supplied material device to 0.0010. The above data show that the sliding mold controller improves the convergence rate of the control error after optimizing the parameters using PSO.

Author Contributions

Conceptualization and Supervision, Z.C.; Methodology and Writing—original draft, S.L.; Project administration, H.G.; Writing—Review and Editing, J.Q.; Software and Validation, Y.Z. and S.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by New Drilling Bars Used in Coal Mine Install Automatic and Fast, Grant Number: 2045150001000048, a Joint Fund Project of the Natural Science Fund of the Autonomous Region of Inner Mongolia, Grant Number: 2023LHMS05029, and Shandong Province’s key research and development project. Grant Number: 2021JMRH0302.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy concerns.

Acknowledgments

We wish to acknowledge Jianguo Qin and Sijia Guo for acquiring funding, calculating data, and contributing during the article revision process. We wish to acknowledge Haixia Gong for providing funding for this project. We wish to acknowledge Zhaoxia Cui and Yunhe Zou for providing data support for this project. We wish to acknowledge Shufang Li for analyzing the data.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yao, Y.F.; Li, X.P.; Zhang, G. Development of An Automatic Loading and Unloading Drilling Bars Device for Coal Mine Pipe Drilling Machines. Coal Min. Mach. 2017, 38, 91–93. [Google Scholar]
  2. Xu, J.; Bone, G.M. Actuators for Improving Robotic Arm Safety While Maintaining Performance: Comparison Study. Actuators 2024, 13, 69. [Google Scholar] [CrossRef]
  3. Sun, Y.; Li, C.; Qin, H.; Deng, Z.; Chen, Z. Robust Neural Network-based Tracking Control for Unmanned Surface Vessels Under Deferred Asymmetric Constraints. Int. J. Robust Nonlinear Control 2021, 32, 2741–2759. Available online: https://onlinelibrary.wiley.com/doi/10.1002/rnc.5520 (accessed on 25 August 2023). [CrossRef]
  4. Subedi, D.; Tyapin, I.; Hovland, G. Dynamic Modeling of Planar Multi-link Flexible Manipulators. Robotics 2021, 10, 70. Available online: https://www.mdpi.com/2218-6581/10/2/70 (accessed on 25 August 2023). [CrossRef]
  5. Meng, Q.; Lai, X.; Yan, Z. Motion Planning and Adaptive Neural Tracking Control of An Uncertain Two-link Rigid–flexible Manipulator with Vibration Amplitude Constraint. IEEE Trans. Neural Netw. Learn. Syst. 2021, 33, 3814–3828. Available online: https://ieeexplore.ieee.org/document/9352495 (accessed on 25 August 2023). [CrossRef] [PubMed]
  6. Ding, R.Q.; Wang, Z.; Chen, M.; Xu, B.; Liu, Z. High-precision Trajectory Tracking of Model-based Hydraulic Mechanical Arm. Mech. Eng. J. 2023, 59, 298–309. Available online: https://kns.cnki.net/kcms2/article/abstract?v=Skeo7MzZydaoEFtRJB1uHkxXeLHhY6_I0CoQNhO3tlQp2cOi9NrAWbjET1qU-USvaTIpiKYK5f3LIRILIOXHe3h_s2mgACTO_0OrH7IGhlN0HsNkAfGDkI0H0t8c1q88R--95mBxmJljH8XeEzAm4g==&uniplatform=NZKPT&language=CHS (accessed on 25 August 2023).
  7. Pujol-Vázquez, G.; Acho, L.; Gibergans-Báguena, J. Commuted PD Controller for Nonlinear Systems: Glucose–Insulin Regulatory Case. Appl. Sci. 2023, 13, 8129. Available online: https://www.mdpi.com/2076-3417/13/14/8129 (accessed on 25 August 2023). [CrossRef]
  8. The Cost Simulation Benchmark-Description and Simulator Manual; Office for Publications of the European Community: Luxembourg, 2001.
  9. Chi, K.H.; Hsiao, Y.F.; Chen, C.C. Robust Feedback Linearization Control Design for Five-Link Human Biped Robot with Multi-Performances. Appl. Sci. 2022, 13, 76. [Google Scholar] [CrossRef]
  10. Feng, J.; Wang, W.; Zeng, H.B. Integral Sliding Mode Control for a Class of Nonlinear Multiagent Systems with Multiple Time-varying Delays. IEEE Access 2024, 12, 10512–10520. [Google Scholar] [CrossRef]
  11. Yu, L.; Huang, J.; Luo, W.; Chang, S.; Sun, H.; Tian, H. Sliding-mode Control for PMLSM Position Control—A review. Actuators 2023, 12, 31. [Google Scholar] [CrossRef]
  12. Li, F.; Zhu, X.; Cao, J.; Yao, B. Kinematics Modeling and Workspace Analysis of a 5-DOF Hydraulic Manipulator. In Proceedings of the 5th International Conference on Control, Automation and Robotics (ICCAR), Beijing, China, 19–22 April 2019; pp. 713–718. [Google Scholar] [CrossRef]
  13. Deng, J.M.; Li, G.D.; Meng, Q.M.; Shen, H.P.; Li, W.H. Modeling and Compensation of Pose Error of 2-RPaRSS Parallel Mechanism Based on D-H Matrix. Mach. Des. Res. 2020, 36, 91–95, 101. Available online: https://webofscience.clarivate.cn/wos/alldb/full-record/CSCD:6686260 (accessed on 25 August 2023).
  14. Zhuo, Z.; Cheng, Z. SCARA Modeling and Simulation based on SimMechanics and Solidworks. In Proceedings of the 2021 IEEE 4th Advanced Information Management, Communicates, Electronic and Automation Control Conference (IMCEC), Chongqing, China, 18–20 June 2021; IEEE: Piscataway, NJ, USA, 2021; Volume 4, pp. 1559–1563. [Google Scholar] [CrossRef]
  15. Wang, C.; Liu, D.; Sun, Q.; Wang, T. Analysis of Open Architecture 6R Robot Forward and Inverse Kinematics Adaptive to Structural Variations. Math. Probl. Eng. 2021, 2021, 4516109. [Google Scholar] [CrossRef]
  16. Cao, G.; Liu, Y.; Jiang, Y.; Zhang, F.; Bian, G.; Owens, D.H. Observer-based Continuous Adaptive Sliding Mode Control for Soft Actuators. Nonlinear Dyn. 2021, 105, 371–386. Available online: https://link.springer.com/article/10.1007/s11071-021-06606-w (accessed on 25 August 2023). [CrossRef]
  17. Li, H.; Liu, H.B. Mechanical Arm Adaptive Fuzzy Sliding Mold Control Based on Interference Observer. Manuf. Technol. Mach. 2023, 115–120. Available online: https://1951.mtmt.com.cn/cn/article/doi/10.19287/j.mtmt.1005-2402.2023.03.015 (accessed on 25 August 2023).
  18. Fu, W.; Wen, H.; Huang, J.H.; Sun, B.; Chen, J.; Chen, W.; Feng, Y.; Duan, X. Underwater Mechanical Arm Adaptive Slide Control Based on Nonlinear Dynamics Model Compensation. J. Tsinghua Univ. Sci. Technol. 2023, 63, 1068–1077. Available online: http://jst.tsinghuajournals.com/article/2023/4336/1687862613628-1587825620.htm (accessed on 25 August 2023).
  19. Chen, G.Q.; Lv, S.B.; Li, G.S.; Dai, J.; Yang, Z.F. Study of Semi-automotive Semi-active Blurry PI Control Suspension Optimized by Particle Swarm Optimization. J. Henan Polytech. Univ. Nat. Sci. 2020; 69–79. [Google Scholar]
  20. Liu, J.K. Intelligent Control. Electron. Ind. Publ. 2017, 116–125. [Google Scholar]
  21. Li, Q.Y.; Cao, Q.S.; Wang, T.T. Parameter Optimization of The Fleet Sliding Mold Controller Based on The Particle Swarm Optimization. Comput. Eng. Des. 2022, 43, 808–813. [Google Scholar]
  22. Jin, T.; Dong, X.C.; Li, Y.N.; Ren, L.; Fan, P.P. Enhanced Particle Swarm Optimization Algorithm to Optimize Fractional Gradient PID Controller Parameters. Comput. Appl. 2019, 39, 796–801. [Google Scholar]
Figure 1. Drilling bar installation and dismantling module.
Figure 1. Drilling bar installation and dismantling module.
Applsci 14 03996 g001
Figure 2. Artificial drilling operations scenario.
Figure 2. Artificial drilling operations scenario.
Applsci 14 03996 g002
Figure 3. Illustration of drilling bars’ placement. 1—Drilling bars. 2—Storage rack. 3—Base seats. 4—Joint1. 5—Joint2. 6—Functional connector. 7—Joint3. 8—Joint4. 9—Mechanical tongs. 10—The front drilling bars. 11—Chuck. 12—Support frame.
Figure 3. Illustration of drilling bars’ placement. 1—Drilling bars. 2—Storage rack. 3—Base seats. 4—Joint1. 5—Joint2. 6—Functional connector. 7—Joint3. 8—Joint4. 9—Mechanical tongs. 10—The front drilling bars. 11—Chuck. 12—Support frame.
Applsci 14 03996 g003
Figure 4. Kinematic structural model of the supplied material device.
Figure 4. Kinematic structural model of the supplied material device.
Applsci 14 03996 g004
Figure 5. MDH coordinate system model.
Figure 5. MDH coordinate system model.
Applsci 14 03996 g005
Figure 6. Initial spatial position of the mechanical tongs.
Figure 6. Initial spatial position of the mechanical tongs.
Applsci 14 03996 g006
Figure 7. Optimization process for sliding mold controller parameters based on PSO.
Figure 7. Optimization process for sliding mold controller parameters based on PSO.
Applsci 14 03996 g007
Figure 8. Iterative change curve of fitness function.
Figure 8. Iterative change curve of fitness function.
Applsci 14 03996 g008
Figure 9. Curves of parameters’ variation.
Figure 9. Curves of parameters’ variation.
Applsci 14 03996 g009
Figure 10. The control structure of the supplied material system.
Figure 10. The control structure of the supplied material system.
Applsci 14 03996 g010
Figure 11. Construction of simulation module.
Figure 11. Construction of simulation module.
Applsci 14 03996 g011
Figure 12. Position tracking curve of joint3.
Figure 12. Position tracking curve of joint3.
Applsci 14 03996 g012
Figure 13. Position tracking curve of joint3 after optimization.
Figure 13. Position tracking curve of joint3 after optimization.
Applsci 14 03996 g013
Figure 14. Position tracking curve of joint4.
Figure 14. Position tracking curve of joint4.
Applsci 14 03996 g014
Figure 15. Position tracking curve of joint4 after optimization.
Figure 15. Position tracking curve of joint4 after optimization.
Applsci 14 03996 g015
Figure 16. Speed tracking curve of joint3.
Figure 16. Speed tracking curve of joint3.
Applsci 14 03996 g016
Figure 17. Speed tracking curve of joint3 after optimization.
Figure 17. Speed tracking curve of joint3 after optimization.
Applsci 14 03996 g017
Figure 18. Speed tracking curve of joint4.
Figure 18. Speed tracking curve of joint4.
Applsci 14 03996 g018
Figure 19. Speed tracking curve of joint4 after optimization.
Figure 19. Speed tracking curve of joint4 after optimization.
Applsci 14 03996 g019
Figure 20. Tracking error curve.
Figure 20. Tracking error curve.
Applsci 14 03996 g020
Figure 21. Tracking error curve after optimization.
Figure 21. Tracking error curve after optimization.
Applsci 14 03996 g021
Table 1. D-H parameter table of the supplied material device.
Table 1. D-H parameter table of the supplied material device.
iai−1 (mm)αi−1 (°)di (mm)θi (°)
10000
20−900−90
3L1000
h(4)L2000
Table 2. Optimized pre- and post-situation comparison.
Table 2. Optimized pre- and post-situation comparison.
Time(s)
Time Spent on
Position Tracking
before Optimizing
Time Spent on
Position Tracking after Optimizing
Saved TimeTime Spent on Speed Tracking
before Optimizing
Time Spent on Speed Tracking after
Optimizing
Saved Time
joint32.242.160.082.622.190.43
joint42.522.030.492.912.030.88
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Qin, J.; Li, S.; Gong, H.; Cui, Z.; Zou, Y.; Guo, S. Trajectory Tracking Algorithm Study of Coal Mine Water Detector Drilling Bar Installation. Appl. Sci. 2024, 14, 3996. https://doi.org/10.3390/app14103996

AMA Style

Qin J, Li S, Gong H, Cui Z, Zou Y, Guo S. Trajectory Tracking Algorithm Study of Coal Mine Water Detector Drilling Bar Installation. Applied Sciences. 2024; 14(10):3996. https://doi.org/10.3390/app14103996

Chicago/Turabian Style

Qin, Jianguo, Shufang Li, Haixia Gong, Zhaoxia Cui, Yunhe Zou, and Sijia Guo. 2024. "Trajectory Tracking Algorithm Study of Coal Mine Water Detector Drilling Bar Installation" Applied Sciences 14, no. 10: 3996. https://doi.org/10.3390/app14103996

APA Style

Qin, J., Li, S., Gong, H., Cui, Z., Zou, Y., & Guo, S. (2024). Trajectory Tracking Algorithm Study of Coal Mine Water Detector Drilling Bar Installation. Applied Sciences, 14(10), 3996. https://doi.org/10.3390/app14103996

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

Article Metrics

Back to TopTop