Next Article in Journal
Intelligent Assessment of Pavement Condition Indices Using Artificial Neural Networks
Next Article in Special Issue
Research on Hydraulic Support Attitude Monitoring Method Merging FBG Sensing Technology and AdaBoost Algorithm
Previous Article in Journal
Study on Balanced Allocation of Water Resources in the Yellow River Basin Based on Water Benefit Sharing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Path Planning and Control Method for Secondary Autonomous Cutting of Cantilever Roadheader in a Large-Section Coal Roadway

1
China Coal Huajin Group Co., Ltd., Yuncheng 043000, China
2
School of Mines, China University of Mining and Technology, Xuzhou 221116, China
3
Wangjialing Coal Mine, China Coal Huajin Group Co., Ltd., Yuncheng 043000, China
*
Authors to whom correspondence should be addressed.
Sustainability 2023, 15(1), 560; https://doi.org/10.3390/su15010560
Submission received: 9 December 2022 / Revised: 24 December 2022 / Accepted: 26 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Advanced and Sustainable Technologies for Tunnel Engineering)

Abstract

:
A cantilever roadheader is the main tunneling equipment for underground coal mine roadways. The key to the safe, efficient and intelligent development of coal enterprises is to achieve the autonomous cutting and intelligent control of the cantilever roadheader. In order to realize the automatic cutting shaping control of a large-section coal roadway, the path planning and control method of secondary automatic cutting of a cantilever roadheader were studied. The Wangjialing 12307 belt roadway was used as the engineering background, the vertical displacement law of the roadway roof under different cutting paths was simulated with the FLAC 3D software, the reasonable cutting path was determined according to the actual situation, and the underground industrial test was carried out. The simplified model and spatial position and attitude coordinate system of the roadheader were established, the kinematics of the roadheader was analyzed, and the position and attitude expression of the cutting head center in the roadway coordinate system was obtained. The simplified model of the cutting head was established, the position expression of the pick in the roadway coordinate system was derived, the position coordinate of the inflection point and the cutting step distance were determined according to the relationship between the cutting head and the roadway boundary, and the cutting path control flow was designed. Finally, the reliability of the cutting path control method was verified with a MATLAB simulation. The research works provide a theoretical foundation for path planning and control to realize “secondary autonomous cutting of cantilever roadheader”.

1. Introduction

Coal is an important energy source for the development of most countries in the world, and China’s annual demand for coal remains high. It can be said that fossil energy such as coal will remain the main energy source in China for a long time to come [1,2,3]. Coal mining mainly adopts the underground mining method, and a large number of roadways need to be excavated to ensure the efficient production of coal, so it can be said that roadway tunneling is the premise of coal mining [4,5]. The main corollary equipment for roadway excavation in coal mines in China is an anchor drilling machine and cantilever roadheader, and the cantilever roadheader integrates the functions of cutting, walking, transportation and dust reduction [6]. However, the tunneling of domestic underground roadways is mainly realized by the manual operation of the driver, and the forming accuracy of the roadway section and the efficiency of the excavation will inevitably be affected by the factors of the roadheader driver. In addition, with the increasing depth of the roadway, the working environment under the shaft is becoming worse and worse, and the personal safety of the workers during tunneling cannot be fully guaranteed. Therefore, the development direction of coal mining enterprises is to realize the automation and intelligence of tunneling equipment, and the autonomous cutting and intelligent control of the roadheader are the key technical problems that need to be solved urgently for the development of coal mining enterprises [7,8,9].
Therefore, many scholars at home and abroad have made numerous explorations on the autonomous cutting and intelligent control methods of roadheaders. Yang et al. took the cutting motor current of the roadheader during cutting as the basis for the change in cutting load, and adopted PID control and adjustment of the expansion and contraction speed of the hydraulic cylinder to propose an intelligent cutting control method for the roadheader cutting section [10]. Zhao et al. combined principal component analysis and neural network to predict the cutting efficiency of a roadheader, which improved the cutting efficiency of the roadheader [11]. Dolipski et al. studied the relationship between cutting head speed and cutting energy consumption of a roadheader [12]. Tuncer et al. used an improved genetic algorithm for the dynamic path planning of mobile robots [13]. Ebrahimabadi et al. and Jasiulek et al. used a GA-BP neural network to enable the roadheader to identify the hardness of coal or rock in front of it when cutting, so as to improve the cutting performance [14,15]. Su et al. established the cutting model of a trapezium roadway interface, determined the inflexion point position of cutting with an approximate iteration method, and simulated it with the MATLAB software, greatly improving the cutting shaping effect of the roadway [16]. Hu et al. simulated the change law of stress and displacement of the surrounding rock of a roadway under different cutting paths by using FLAC 3D, and determined the optimal cutting path in combination with the actual situation in the mine [17]. Xu et al. established the spatial pose model of a cantilever roadheader by using the D-H parameter method, deduced the spatial pose matrix of roadheader and cutting head, determined the coordinates of the inflection point of the cutting path, and finally simulated it with MATLAB [18]. Zhang et al. realized the adjustment of the pitching angle of the roadheader during tunneling by combining the particle swarm optimization algorithm and PID [19]. Yang et al. used a genetic algorithm and a mutation particle swarm optimization algorithm to optimize the path target point set, effectively improving the roadheader’s attitude planning capability [20].
Although predecessors have conducted a lot of research, it was mainly based on the single cutting of a roadheader. For large-section roadways, it is difficult for the current mainstream cantilever roadheader to cut the roadway into shape at one time, and roadways needed to be formed by secondary cutting, which required research on the control method of the roadheader’s secondary autonomous cutting roadways.
Therefore, based on the research background of the Wangjialing 12307 belt roadway, this paper studied the path planning and control method of the secondary autonomous cutting of a cantilever roadheader. Based on the geological conditions of Wangjialing, the roadway tunneling model was established using FALC 3D, the vertical variation law of the surrounding rock of the roadway under different cutting paths was simulated, the reasonable cutting path was determined in combination with the actual situation in the mine, and the underground industrial experiments were carried out. Then, the space pose model of the cantilever roadheader was established, the space pose transformation matrix of the roadheader cutting head was derived using the D-H parameter method, and the limit cutting area of the cutting head was simulated using MATLAB. Finally, the inflexion coordinates of the reasonable cutting path were determined, the cutting path control method was proposed, and the cutting path control process was simulated with MATLAB. This research provide a theoretical foundation for path planning and control to realize the “secondary autonomous cutting of cantilever roadheader”.

2. Cutting Path Planning

The roadway section is formed after being cut by the cutting head according to a certain cutting path, and different cutting paths will have different effects on the front coal wall and the roof of the empty roof area, thus affecting the actual efficiency of the roadheader and the personal safety of the workers. Therefore, determining a reasonable cutting path is an important means to ensure the high efficiency and safety of roadway excavation. This chapter took the Wangjialing 12307 belt roadway as the engineering background, and used the FLAC 3D software to simulate and analyze the vertical displacement of the surrounding rock of the roadway under different cutting paths.

2.1. Project Overview

The section size of the Wangjialing 12307 belt roadway is 5.6 m × 3.55 m, belonging to a large-section roadway. As the cantilever roadheader cannot shape the roadway after one cutting, the method of secondary cutting is adopted. First cut the roadway section on the left, then move the roadheader to a suitable place, and then cut the roadway section on the right. The schematic diagram of the secondary cutting roadway is shown in Figure 1.

2.2. Model Establishment

According to Saint Venant’s principle, the influence intensity of roadway excavation will decrease with an increase in the distance from the excavation point. It is generally believed that the influence range is 3–5 times the excavation radius of the roadway [21]. Accordingly, a roadway model with a size of 50 m × 10 m × 40 m was established, and the grid density near the coal seam was properly densified to improve the accuracy of the simulation. The model is shown in Figure 2.
According to the physical and mechanical test report of the Wangjialing 12307 belt roadway, the geological parameters of each coal bed were assigned. The numerical simulation geological parameters are shown in Table 1. A vertical stress of 7.35 MPa on the top boundary of the model was applied, and the Mohr Coulomb yield criterion was used for the simulation.

2.3. Simulation Analysis of Cutting Path

The current mainstream cutting paths of cantilever roadheaders include “loop” and “snake”. In order to reduce the idling time of the cutting head, four cutting paths were designed as follows, as shown in Figure 3.

2.3.1. “Snake” Cutting Path Simulation

Cutting path 1 and cutting path 2 were simulated, and the vertical displacement of roadway roof under the two cutting paths was analyzed, as shown in Figure 4 and Figure 5.
(1)
Simulation analysis of cutting path 1
As shown in Figure 4, the maximum vertical displacement of the roadway roof increases with the increase in cutting times. When the roadway is shaped, the roadway roof displacement reaches the maximum value of about 126 mm. This is because with the increase in cutting times, the two sides of the roadway are gradually formed, and the displacement of the two sides increases, leading to an increase in the roof displacement. In addition, with the gradual shaping of the two sides of the roadway, the maximum vertical displacement amplitude of the roadway roof gradually slows down, and the coal body gradually becomes stable.
(2)
Simulation analysis of cutting path 2
The cutting process under cutting path 1 was similar to that under cutting path 2, only the cutting path was opposite. Therefore, in order to reduce the space occupied by the simulation results, the roof vertical displacement nephogram after the fifth and tenth cutting was selected for analysis.
As shown in Figure 5, the maximum vertical displacement of the roadway roof increases with the increase in cutting times. When the roadway is shaped, the maximum displacement of the roadway roof reaches about 123 mm.

2.3.2. “Loop” Cutting Path Simulation

Cutting paths 3 and 4 were simulated, and the vertical displacement of the roadway roof under the two cutting paths was analyzed, as shown in Figure 6 and Figure 7.
(1)
Simulation analysis of cutting path 3
According to an analysis of Figure 6, the maximum displacement of the roadway roof after the first cut is 96 mm. After the second cut, the maximum displacement of roadway roof reaches 99 mm. After the third cut, the roadway floor is formed, the maximum displacement of the roof is 105 mm, and the maximum displacement of the floor is 40 mm. The left side is formed after the fourth cut, and the maximum displacement of the roadway roof is 107 mm, with a small change. With the increase in cutting times, the displacement of the roadway roof increases continuously, and the maximum displacement of the roadway roof after cutting is 126 mm.
(2)
Simulation analysis of cutting path 4
Similarly, as the cutting process under cutting path 4 was similar to that under cutting path 3, in order to reduce space, the roof vertical displacement nephogram with representative cutting times was selected for analysis.
According to an analysis of Figure 7, after the first cut, the roadway floor is formed and the floor displacement is small. After the third cut, the roof is formed, and the roof displacement is 95 mm. The left side is formed after the fourth cut, and the drift roof displacement is 98 mm, with a small change. With the increase in cutting times, the displacement of the roadway roof also increases correspondingly, but the displacement variation amplitude of the roadway roof is smaller than that under working conditions 1~3. The maximum displacement of the roadway roof after cutting is 120 mm.

2.4. Determination of Cutting Path

According to the statistical analysis and numerical simulation Figure 4, Figure 5, Figure 6 and Figure 7, the maximum displacement and displacement change amplitude of the roadway roof of cutting path 4 were minimal, and the roof was easy to control during roadway cutting. Combined with the actual situation of Wangjialing and the numerical simulation results, the cutting path is determined as shown in Figure 8 below.
In Figure 8, the dotted line represents the motion track of the cutting head, the black arrow represents the cutting direction of the cutting head, and the numerical sequence represents the sequence of the central position of the cutting head.
The cutting head is drilled at the left side of the roadway. After the cutting head has fully entered, the cutting arm swings from left to right to 3300 mm from the left side of the roadway, then the cutting arm swings from bottom to top and cuts along the “reverse loop” cutting path. The cantilever roadheader adjusts its position for the second cutting. The cutting head is drilled 1600 mm from the right side of the roadway. After the cutting head is fully entered, the cutting arm cuts along the “reverse loop” cutting route.

2.5. Industrial Test of Cutting Path

The industrial test of the cutting path was conducted in the belt roadway of the Wangjialing 12307 working face. According to the field measurement, after the roadway was cut and formed, the cutting section was flat and the forming effect was good. The roadheader could remain stable during the cutting process and was convenient for the roadheader driver to operate. The roadway cutting section is shown in Figure 9.

3. Kinematics Analysis of Cantilever Roadheader

The kinematics of the roadheader can be divided into two parts: cutting mechanism kinematics and body kinematics. The kinematics of the cutting mechanism mainly describes the spatial position and attitude relationship between the cutting head and the body during the cutting process, while the kinematics of the body mainly describes the spatial position and attitude relationship between the body and the roadway during the cutting process of the roadheader. The cutting mechanism and body kinematics analysis are the theoretic foundation for the roadheader to realize autonomous cutting and control, which has a significant influence on the cutting and shaping effect of the roadway section.

3.1. Establishment of Space Position and Attitude Coordinate System

The cutting head, rotary table, cantilever, hydraulic cylinder and traveling mechanism of a roadheader can be regarded as consisting of a series of translated or rotated joints and connecting rods [22]. Based on this, a simplified model of the robotization cantilever roadheader and its spatial position and attitude coordinate system were established [23], as shown in Figure 10.
In the figure: Oh-XhYhZh is the roadway coordinate system, the OhZh axis points vertically to the roadway roof, the OhXh axis points to the heading direction of the roadheader, and the OhYh axis is parallel to the roadway floor.
O0-X0Y0Z0 is the center of the gravity coordinate system of the roadheader body, which is convenient for establishing the position and attitude transformation relationship among the cutting head, body and roadway coordinate system.
0i-XiYiZi (i = 1, 2, 3, 4) is the coordinate system of the cutting mechanism of the roadheader, which is established at each link joint. 01-X1Y1Z1 is the rotary table coordinate system, 02-X2Y2Z2 is the cantilever coordinate system, 03-X3Y3Z3 is the cutting head expansion part coordinate system, and 04-X4Y4Z4 is the cutting head coordinate system.
In addition, a, d, θ and α are a joint variable. ai is the distance of translation along the OiXi axis so that Zi and Zi+1 coincide, θi is the angle of rotation around the OiZi axis to make the Xi−1 axis and Xi axis collinear, di is the distance of translation along the OiZi axis to make the Xi−1 axis and Xi axis parallel, αi is the angle of rotation around the Xi axis to make Zi and Zi+1 collinear, and d is the retraction stroke of the cutting head.

3.2. Kinematics Analysis of Cutting Mechanism

3.2.1. Forward Kinematics Solution of Cutting Mechanism

The forward kinematics problem of the cutting mechanism of the roadheader is to solve the spatial position and posture of the cutting head in different coordinate systems when the geometric parameters and joint variables of each member of the cutting mechanism is known. The D-H parameter method can be used to solve the forward kinematics problem, to realize the spatial position and orientation of the cutting head.
Based on the principle of the D-H parameter method [24], the coordinate vector {xi, yi, zi, 1}T of any point in the coordinate system Oi-XiYiZi and the coordinate vector {xi−1, yi−1, zi−1,1}T in the relative coordinate system Oi-1-Xi-1Yi-1Zi-1 have the following relationship:
{ x i 1 , y i 1 , z i 1 , 1 } T = T i i 1 { x i , y i , z i , 1 } T
where T i i 1 is the position and attitude transformation matrix of the coordinate system Oi-XiYiZi relative to the coordinate system Oi-1-Xi-1Yi-1Zi-1.
According to Formula (1), the coordinate vector {x4, y4, z4,1}T of any point in the coordinate system O4-X4Y4Z4 of the cutting head and the coordinate vector {x0, y0, z0,1}T in the coordinate system O0-X0Y0Z0 of the fuselage center of gravity have the following relationship:
{ x 0 , y 0 , z 0 , 1 } T = T 1 0 · T 2 1 · T 3 2 · T 4 3 { x 4 , y 4 , z 4 , 1 } T = T 4 0 { x 4 , y 4 , z 4 , 1 } T
According to Figure 10, the linkage parameters of the cutting mechanism were established, as shown in Table 2.
According to the parameter table of the linkage of the cutting mechanism and Formula (2), the position and attitude transformation matrix of the cutting head relative to the gravity center coordinate system T 4 0 of the fuselage was obtained:
T 4 0 = T 1 0 · T 2 1 · T 3 2 · T 4 3 = c θ 1 c θ 2 s θ 1 c θ 1 s θ 2 a 0 + a 1 c θ 1 + d + d 3 c θ 1 s θ 2 + a 2 c θ 1 c θ 2 + d 4 c θ 1 s θ 2 c θ 2 s θ 1 c θ 1 s θ 1 s θ 2 a 1 s θ 1 + d + d 3 s θ 1 s θ 2 + a 2 c θ 2 s θ 1 + d 4 s θ 1 s θ 2 s θ 2 0 c θ 2 d 1 d + d 3 c θ 2 d 4 c θ 2 + a 2 s θ 2 0 0 0 1
where c represents cos, s represents sin, T i i 1 ( i = 1 ,   2 ,   3 ,   4 ) is the posture transformation matrix of the cutting mechanism coordinate system i relative to the cutting mechanism coordinate system i−1, and
T 1 0 = [ cos θ 1 sin θ 1 0 a 0 sin θ 1 cos θ 1 0 0 0 0 1 d 1 0 0 0 1 ] , T 2 1 = [ cos θ 2 sin θ 2 0 a 1 0 0 1 0 sin θ 2 cos θ 2 0 0 0 0 0 1 ] T 3 2 = [ 1 0 0 a 2 0 0 1 ( d 3 + d ) 0 1 0 0 0 0 0 1 ] , T 4 3 = [ 1 0 0 0 0 1 0 0 0 0 1 d 4 0 0 0 1 ]
Equation (3) can be substituted into Equation (2) to obtain the coordinates {x0, y0, z0} of the cutting head center O4 (at this time, x4 = y4 = z4 = 0,) in the fuselage gravity center coordinate system O0-X0Y0Z0 as follows:
x 0 = a 0 + a 1 cos θ 1 + cos θ 1 sin θ 2 d + d 3 + a 2 cos θ 1 cos θ 2 + d 4 cos θ 1 sin θ 2 y 0 = a 1 sin θ 1 + d + d 3 sin θ 1 sin θ 2 + a 2 cos θ 2 sin θ 1 + d 4 sin θ 1 sin θ 2 z 0 = d 1 cos θ 2 d + d 3 d 4 cos θ 2 + a 2 sin θ 2

3.2.2. Inverse Kinematics Solution of Cutting Mechanism

The inverse kinematics problem of the cutting mechanism of the roadheader is to know the geometric parameters of each link of the cutting mechanism and solve the joint variables of each link when the cutting head is positioned relative to different coordinate systems. The inverse kinematics solution of the cutting mechanism is the basis for the planning and precise control of the cutting head’s moving track.
The inverse kinematics of the cutting mechanism can be solved with the inverse transformation method [23]. The position and attitude transformation matrix T 4 0 of the cutting head relative to the gravity center coordinate system of the fuselage could also be expressed as:
T 4 0 = [ n x o x a x p x n y o y a y p y n z o z a z p z 0 0 0 1 ] = [ n o a p 0 0 0 1 ]
where a is the approach vector, o is the azimuth vector, n is the normal vector, and p is the position vector. Three orthogonal unit column vectors (n, a and o) describe the position and orientation of the cutting head.
It could be seen from the analysis of Equation (5) that each variable in the pose transformation matrix of Equation (3) should be equal to the corresponding element in the matrix composed of joint variables of Equation (5). From these equations, the relationship between joint variables and n, o, a, p vectors could be obtained.
From the equality of elements (1,2) and (2,2) in Equation (3) and Equation (5), it could be obtained that θ 1 = arctan o x o y or θ 1 = θ 1 + 180 ° .
According to the equality of (3,1) and (3,3) elements in Equation (3) and Equation (5), we could obtain θ 2 = arctan n z a z or θ 2 = θ 2 + 180 ° .
According to the corresponding equality of (1,4), (2,4) and (3,4) elements of Equation (4) and Equation (9), we could obtain d = p x a 0 + ( a 1 + a 2 ) o y a 2 n x a x d 3 d 4 or d = p y a 1 o x a 2 n y a y d 3 d 4 or d = a 2 n z + p z d 1 a z d 3 d 4 . Although there were three formulas to solve d, the formula whose denominator was not 0 could be used to solve parameter d.

3.3. Kinematic Analysis of Roadheader Body

The correct expression of the position and posture of the roadheader body is the basis for the roadheader to achieve accurate and autonomous cutting. The method of describing the position and attitude of aircraft, ships and other vehicles during navigation was used for reference. The “fuselage yaw angle α”, “fuselage pitch angle β”, “fuselage rolling angle γ” and “fuselage offset L” were adopted to describe the position and attitude of the roadheader when tunneling in the roadway coordinate system [25,26]. The design centerline direction of the roadway is defined as the Xh axis; the direction vertically pointing to the roadway roof is taken as the axis Zh; the Yh axis is parallel to the roadway floor; the angle of rotation about the Zh axis is called the yaw angle α; the angle of rotation about the Yh axis is called the pitch angle β; the angle of rotation about the Xh axis is called the roll angle γ; the three-dimensional space offset distance between the fuselage center line and the roadway center line is called offset L, as shown in Figure 11.
According to the coordinate transformation theory, the fuselage center of gravity coordinate system O0-X0Y0Z0 can be obtained from the roadway coordinate system Oh XhYhZh through the following transformation: first rotate around the OhZh axis α angle, then rotate around the OhYh axis β angle, then rotate around the OhXh axis γ angle, and finally move the vector L (△X, △Y, △Z) along the OhOh axis. From this, the position and attitude transformation matrix of the fuselage gravity center coordinate system relative to the roadway coordinate system was obtained as follows:
T 0 h = Rot ( Z h , α ) · Rot ( Y h , β ) · Rot ( X h , γ ) · Trans ( X ,   Y ,   Z ) = [ c α c β s α c γ + c α s β s γ s α s γ + c α s β c γ c α c β X + ( c α s β s γ s α c γ ) Y + ( s α s γ + c α s β c γ ) Z s α c β c α c γ + s α s β s γ c α s γ + s α s β c γ s α c β X + ( c α c γ + s α s β s γ ) Y + ( s α s β c γ c α s γ ) Z s β c β s γ c β c γ s β X + c β s γ Y + c β c γ Z 0 0 0 1 ]
where c is cos and s is sin.

3.4. Kinematics Analysis of Roadheader

In order to realize the transformation from the cutting head coordinate system to the roadway coordinate system, the kinematics analysis of the roadheader was carried out. The pose transformation matrix of the cutting head coordinate system relative to the roadway coordinate system was denoted as T 4 h , then:
T 4 h = T 0 h · T 4 0
Then, the coordinate vector {xh, yh, zh, 1}T of the cutting head center O4 in the roadway coordinate system Oh-XhYhZh was:
{ x h , y h , z h , 1 } T = T 4 h · { 0 , 0 , 0 , 1 } T = T 0 h · T 4 0 · { 0 , 0 , 0 , 1 } T
In order to simplify the calculation and simulation analysis, the following assumptions were made: the roadheader kept the body stable during the secondary tunneling process, without any deflection, and only moved along the OhYh axis; the center of gravity line of the body of the roadheader at the first tunneling section was set as the center line of the roadway, that is, the gravity center line of the roadheader body coincided with the center line of the roadway during the first tunneling, as shown in Figure 12.
According to the above assumptions, α = β = γ = 0 ° and X = Y = Z = 0 for the first tunneling of the roadheader, and α = β = γ = 0 ° , X = Z = 0 , and Y = 2800 during the second tunneling of the roadheader.
Therefore, T 0 h could be simplified as:
T 0 h = [ 1 0 0 0 0 1 0 Y 0 0 1 0 0 0 0 1 ]
Equation (9) and Equation (3) can be substituted into Equation (8) to obtain the coordinate expression of the cutting head center O4 in the roadway coordinate system:
x h = a 0 + a 1 cos θ 1 + cos θ 1 sin θ 2 ( d + d 3 ) + a 2 cos θ 1 cos θ 2 + d 4 cos θ 1 sin θ 2
y h = Y + a 1 sin θ 1 + ( d + d 3 ) sin θ 1 sin θ 2 + a 2 cos θ 2 sin θ 1 + d 4 sin θ 1 sin θ 2
z h = d 1 cos θ 2 d + d 3 d 4 cos θ 2 + a 2 sin θ 2

4. Cutting Path Control

4.1. Pick Coordinate Determination

In fact, the roadway section is formed by the continuous rotation of the cutting head to drive the pick installed on the cutting head to continuously cut the coal wall. Therefore, the boundary of the roadway is formed when the rotation curve of the pick is tangent to the boundary of the roadway section. In order to obtain the curve of pick rotation, the cutting head was regarded as the envelope surface formed by a parabola rotating around the cutting head axis. The simplified model of the cutting head is shown in Figure 13.
In Figure 13, i is any point on the parabola (can be regarded as a pick), and r is the distance from point i to the cutting head axis. Let the equation of parabola be Z = aX2 + b; the trajectory equation of pick i rotation then is:
{ x = r i cos φ y = r i sin φ   z = a r i 2 + b
Let the coordinate vector of the pick i in the cutting head coordinate system O4-X4Y4Z4 be:
I 4 = { x i 4 , y i 4 , z i 4 , 1 }
Combining Formula (11), we could obtain:
{ x i 4 = r i cos φ y i 4 = r i sin φ   z i 4 = a r i 2 + b
Therefore, when the cutting head shaft length, radius and cutting depth are known, ri can be determined according to Equation (13). As the type and cutting depth of the roadheader were designed in advance during cutting, ri could be regarded as a known quantity in the subsequent calculation.
According to the D-H parameter method, the coordinate vector { x i h , y i h , z i h , 1 } of pick i in the roadway coordinate system Oh XhYhZh was:
{ x i h , y i h , z i h , 1 } = T 0 h · T 4 0 · I 4
The coordinate expression of pick i in the roadway coordinate system Oh-XhYhZh was obtained by substituting Formulas (3), (9) and (13) into Formula (14):
x i h = cos θ 1 cos θ 2 r i cos φ + sin θ 1 r i sin φ + cos θ 1 sin θ 2 ( ar i 2 + b ) + a 0 + a 1 cos θ 1 + cos θ 1 cos θ 2 ( d + d 3 ) + a 2 cos θ 1 cos θ 2 + d 4 cos θ 1 sin θ 2
y i h = Y + cos θ 2 sin θ 1 r i cos φ cos θ 1 r i sin φ + sin θ 1 sin θ 2 ( ar i 2 + b ) + a 1 sin θ 1 + ( d + d 3 ) sin θ 1 sin θ 2 + a 2 cos θ 2 sin θ 1 + d 4 sin θ 1 sin θ 2
z i h = sin θ 2 r i cos φ cos θ 2 ( ar i 2 + b ) + d 1 cos θ 2 ( d + d 3 ) d 4 cos θ 2 + a 2 sin θ 2  

4.2. Determination of Inflection Point of Cutting Path

Research shows that the inflection point position of the cutting path will have a significant impact on the forming quality of the cutting section [27]. The inflection point position of the path was determined according to the planned cutting path (Figure 8) so as to achieve accurate planning and control of the cutting path. The schematic diagram of the inflection point of the cutting path is shown in Figure 14. In the figure, L is the set cutting step, h1hupper and h1hlower are the vertical coordinates of the upper and lower boundary of the left half of the roadway, z1hupper and z1hlower are the vertical coordinates of the cutting head center when cutting the upper and lower boundary of the left half of the roadway, and b1hright is the horizontal coordinates of the bottom of the right boundary of the left half of the roadway.

4.2.1. Determination of Inflection Point Coordinates of Boundary Contour

When the cutting head was cutting the upper boundary of the left half of roadway, by substituting φ = 0 into Equation (15c) the following was obtained:
z i h = sin θ 2 r i cos θ 2 ( ar i 2 + b ) + d 1 cos θ 2 ( d + d 3 ) d 4 cos θ 2 + a 2 sin θ 2
When the cutting head was cutting the lower boundary of the left half of the roadway, by substituting φ=π into Equation (15c) the following was obtained:
z i h = sin θ 2 r i cos θ 2 ( ar i 2 + c ) + d 1 cos θ 2 ( d + d 3 ) d 4 cos θ 2 + a 2 sin θ 2
The upper ordinate h1hupper = 3550 and the lower ordinate h1hlower = 0 were substituted into Equation (16) and Equation (17) to calculate the θ2, and θ2 was then substituted into Formula (10c) to calculate the corresponding cutting head center coordinates z1hupper and z1hlower, and the vertical coordinates of inflection points 1~4 could be obtained.
In the left half of the roadway, the abscissa of the inflection point was symmetrical; therefore, only the abscissa of one inflection point needed be calculated. When the cutting head was cutting the lower right boundary of the left half of the roadway, Y = 0 and φ = π 2 were substituted into Formula (15b) to obtain:
y i h = cos θ 1 r i + sin θ 1 sin θ 2 ( ar i 2 + c ) + a 1 sin θ 1 + ( d + d 3 + d 4 ) sin θ 1 sin θ 2 + a 2 cos θ 2 sin θ 1  
The corresponding θ2 was substituted when the cutting head cut the lower boundary of the left half of the roadway and the abscissa b1hright = 1825 into Formula (18); θ1 could then be obtained, and then the abscissa of inflection points 2 and 3 could be obtained from Formula (10b). The abscissa of inflection points 1 and 4 could be obtained according to the symmetry. At this time, the coordinates of inflection points 1–4 in the left half of the roadway were calculated.
Similarly, the coordinates of the inflection point on the right half of the roadway could be obtained, but Y = 2800 needed to be substituted into the calculation.

4.2.2. Determination of Cutting Step

When the roadheader cuts, if the cutting step distance L is set too small, although the cutting accuracy of the roadway section can be improved, the cutting steps required for roadway cutting and forming will increase, resulting in an increased cutting time. If the cutting step L is set too large, although it will reduce the cutting time, the cutting accuracy of the roadway section will decline, resulting in uneven roadway sides. In order to determine the cutting step distance L, the model diagram of the error between the cutting step distance and the roadway boundary was established, as shown in Figure 15.
It could be seen from the Figure 15 that the cutting error a, the cutting step L and the cutting radius ri had the following relations:
a = r i r i 2 ( L 2 ) 2
Formula (19) showed that when the cutting radius ri remained unchanged, the cutting error a increased with an increase in the cutting step L. The cutting step L should be greater than ri but not greater than 2ri. In addition, the cutting radius ri could be taken as a known quantity, so that when the maximum cutting error a was specified, the range of cutting step L could be determined.

4.2.3. Determination of Internal Inflection Point Coordinates

In order to form a complete and smooth cutting roadway section, the four inflection points around the roadway were considered as the benchmark, a distance of step L was gradually cut inward. Then, the coordinates of inflection point 5 were (x1, y1+L), inflection point 6 was (x2−L, y2+L), inflection point 7 was (x3−L, y3−L), inflection point 8 was (x4+L, y4−L), inflection point 9 was (x8, y8−L), and inflection point 10 was (x9+L, y9).
Similarly, the coordinates of the remaining inflection points in the right half of the roadway were as follows: the coordinates of inflection point 15 were (x11, y11+L), the coordinates of inflection point 16 were ( x 11 + x 12 2 , y12+L), and the coordinates of inflection point 17 were ( x 11 + x 12 2 , y13−L).

4.3. Cutting Path Control

Based on the determined inflection point coordinates of the cutting path, the cutting path control was carried out to realize the automation of the cutting process. When the cutting head reached the specified position in the horizontal direction, the rotary oil cylinder stopped working, and the lifting oil cylinder was opened to enable the cantilever to conduct a vertical cutting movement. When the cutting head reached the specified position in the vertical direction, the lifting oil cylinder stopped working, and the rotary oil cylinder began to make the cantilever perform horizontal cutting. The control program flow is shown in Figure 16.

5. Results and Discussion

5.1. Simulation Analysis of Cutting Head’s Limit Cutting Area

The relevant dimensions of the EBZ-200 cantilever roadheader were substituted into Formula (10), and the MATLAB software was used to simulate the limit cutting area of the cutting head. The simulation results are shown in Figure 17.
The simulation results in Figure 17 show the cutting process of the cutting head. Whether it was single cutting or secondary cutting is considered, the surroundings of the roadway were flat and smooth after the cutting formation, and the shape of the limit cutting section was that of a waist drum. The single cutting limit cutting area was 4100 mm × 3900 mm, the limit cutting area of the secondary cutting was 5650 mm × 3900 mm, and the ability to achieve a 5.6 m × 3.55 m-large roadway through cutting formation proved the validity of the theory. In addition, it can be seen from Figure 17d that there were some overlapping cutting areas, which decreased with increases in the moving distance ( Y ) of the roadheader, and the corresponding limit area of the secondary cutting also increased.

5.2. Simulation of Cutting Path Control

The MATLAB software was used to conduct the cutting path control simulation, as shown in Figure 18. In the simulation, the shaft length of the cutting head was set as 1000 mm, the maximum radius was 500 mm, the cutting depth was 600 mm, and the cutting step was 700 mm.
It can be seen from Figure 18 that the simulated cutting path was consistent with the set cutting path, and the cutting error around the roadway was controlled within 1 mm, which improved the cutting accuracy and realized the planning and high-precision control of the cutting path during secondary cutting.
The cutting path planning and control method proposed in this paper for the secondary cutting of a cantilever roadheader is still applicable to the cutting path planning and control of “snake” and other large-section roadways, although it was an example of a “loop” cutting path.
In addition, when it was found that there was a fault in front of the coal wall or the dip angle of the coal seam changed, the yaw angle α, pitch angle β, roll angle γ and offset L of the roadheader body were adjusted in time. Through the adjustment of the position and attitude parameters of the roadheader body, the tunneling direction of the roadheader was changed. When the position and attitude parameters of the roadheader body changed, the position coordinates of the calculated inflection point also changed accordingly, thereby changing the cutting path.

6. Conclusions

(1)
Based on the process of secondary cutting into the roadway, the Wangjialing 12307 belt roadway was used as the engineering background, the vertical displacement law of the roadway roof under different cutting paths was simulated with the FLAC 3D software, and the reasonable cutting path of the secondary cutting roadway was determined in combination with the actual situation in the mine. The effectiveness of the cutting path was also verified through industrial tests.
(2)
According to the robot principle and the D-H parameter method, the spatial pose model and coordinate system of the roadheader were established, which provided a mathematical model and a global coordinate system for the kinematics analysis of the cutting mechanism, the body and the roadheader.
(3)
The cutting structure kinematics, body kinematics and roadheader kinematics were analyzed and solved, and the pose transformation matrix of the cutting head coordinate system relative to the roadway coordinate system and the pose expression of the cutting head center in the roadway coordinate system were obtained, which provided a kinematic theoretical basis for the roadheader’s secondary autonomous cutting control system.
(4)
The simplified model of the cutting head was established, the position and pose expression of the pick on the cutting head in the roadway coordinate system was derived, the position coordinates of the inflection point and the cutting step distance were determined according to the relationship between the cutting head and the roadway boundary. The cutting path control flow was also designed, which provided the cutting path control method for the secondary autonomous cutting control system of the roadheader.
(5)
The MATLAB software was used to simulate the limit cutting area and cutting path control process of the cutting head during the second cutting. The simulation results showed that the limit cutting area was 5650 mm × 3900 mm, and the cutting error around the roadway was controlled within 1 mm, which verified the effectiveness and reliability of the automatic cutting control method. This research provides a theoretical basis for path planning and control to realize “secondary autonomous cutting of cantilever roadheader”.

Author Contributions

J.W., Z.X., X.F. and G.S. designed the study and experiment; J.W., Z.X. and X.F. collected the data; Z.X., X.F., G.S. and H.W. conducted the data analysis; Z.X., J.W. and G.S. provided the statistical methods; Z.X., J.W., X.F. and G.S. drafted the paper; Z.X., J.W. and H.W. edited the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51874276).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data and code used or analyzed in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Milici, R.C.; Flores, R.M.; Stricker, G.D. Coal resources, reserves and peak coal production in the United States. Int. J. Coal Geol. 2013, 113, 109–115. [Google Scholar] [CrossRef]
  2. Lu, Y.; Xiao, S.; Ge, Z.; Zhou, Z.; Deng, K. Rock-breaking properties of multi-nozzle bits for tree-type drilling in underground coal mines. Energies 2016, 9, 249. [Google Scholar] [CrossRef] [Green Version]
  3. Su, H.; Zhang, H.; Qi, G.; Lu, W.; Wang, M. Preparation and CO2 adsorption properties of TEPA-functionalized multi-level porous particles based on solid waste. Colloids Surf. A Physicochem. Eng. Asp. 2022, 653, 130004. [Google Scholar] [CrossRef]
  4. Xie, Z.; Zhang, N.; Qian, D.; Han, C.; An, Y.; Wang, Y. Rapid Excavation and Stability Control of Deep Roadways for an Underground Coal Mine with High Production in Inner Mongolia. Sustainability 2018, 10, 1160. [Google Scholar] [CrossRef] [Green Version]
  5. Qi, P.; Chang, J.; Chen, X.; Wang, T.; Wu, M. Identification of Rock Properties of Rock Wall Cut by Roadheader Based on PSO-VMD-LSSVM. Front. Earth Sci. 2022, 10, 142. [Google Scholar] [CrossRef]
  6. Zhang, W.; Zhang, Y.; Zhu, F.; Liang, Z. Experimental Study on Reflection Superiority of Multisensors Measurement-While-Drilling of Roof Bolter. IEEE Trans. Instrum. Meas. 2022, 71, 1–11. [Google Scholar] [CrossRef]
  7. Ji, X.; Zhang, M.; Qu, Y.; Jiang, H.; Wu, M. Travel Dynamics Analysis and Intelligent Path Rectification Planning of a Roadheader on a Roadway. Energies 2021, 14, 7201. [Google Scholar] [CrossRef]
  8. Yang, J.; Zhang, Q.; Wang, C. Research status and development of robotization of coal mine roadheader. J. China Coal 2020, 3, 1–12. [Google Scholar]
  9. Zong, K.; Zhang, P.; Wang, P. Multifactor Analysis of roadheader’s body pose responses during the horizontal cutting process. Shock. Vib. 2018, 2018, 18. [Google Scholar] [CrossRef] [Green Version]
  10. Yang, W.; Wang, Z. Research on the adaptive control for cutting operation of roadheader. Appl. Mech. Mater. 2013, 270, 1436–1439. [Google Scholar]
  11. Zhao, L.; Wang, J. Cutting performance evaluation of a roadheader machine by PCA and RBF. Simul.-Transcation Soc. Model. Simul. Int. 2018, 94, 1129–1141. [Google Scholar] [CrossRef]
  12. Dolipski, M.; Cheluszka, P.; Sobota, P. Investigating the simulated control of the rotational speed of roadheader cutting heads, relating to the reduction of energy consumption during the cutting process. J. Min. Sci. 2015, 51, 298–308. [Google Scholar] [CrossRef]
  13. Tuncer, A.; Yildirim, M. Dynamic path planning of mobile robots with improved genetic algorithm. Comput. Electr. Eng. 2012, 38, 1564–1572. [Google Scholar] [CrossRef]
  14. Ebrahimabadi, A.; Azimipour, M.; Ali, B. Prediction of roadheaders’ performance using artificial neural network approaches (MLP and KOSFM). J. Rock Mech. Geotech. Eng. 2015, 7, 573–583. [Google Scholar] [CrossRef] [Green Version]
  15. Jasiulek, D.; Stankiewicz, K. An adaptive control system of roadheader with intelligent modelling of mechanical features of mined rock. J. KONES Powertrain Transp. 2011, 18, 197–203. [Google Scholar]
  16. Ning, S.; Su, H.; Gao, J. Research on Automatic Section Precision Forming of Boom-type Roadheader. In Proceedings of the 2015 International Conference on Intelligent Systems Research and Mechatronics Engineering, Zhengzhou, China, 11–13 April 2015. [Google Scholar]
  17. Hu, C.; Zhang, Y.; Yu, R.; Fang, X.; Xu, Z.; Wang, L.; Zhao, B. Research on the Optimization of Cutting Path of Cantilever Roadheader in Large Section Excavation. Sustainability 2022, 14, 5345. [Google Scholar] [CrossRef]
  18. Xu, Z.; Liang, M.; Fang, X.; Wu, G.; Chen, N.; Song, Y. Research on Autonomous Cutting Method of Cantilever Roadheader. Energies 2022, 15, 6190. [Google Scholar] [CrossRef]
  19. Zhang, M.; Lyu, F.; Fu, S.; Cai, X.; Zong, K.; Wu, M. Study on the pitch angle control of a robotized hydraulic roadheader using different control methods. J. Mech. Sci. Technol. 2018, 32, 4893–4901. [Google Scholar] [CrossRef]
  20. Yang, J.J.; Tang, Z.H.; Wang, X.L.; Wang, Z.R.; Wu, M. Study on the stochastic constraint environment modeling of the deviation correction path of coal roadheader. J. Min. Sci. 2019, 4, 58–65. [Google Scholar]
  21. Tsinidis, G.; Pitilakis, K.; Madabhushi, G. On the dynamic response of square tunnels in sand. Eng. Struct. 2016, 125, 419–437. [Google Scholar] [CrossRef] [Green Version]
  22. Jun, M.; Jian-gang, L.; Wei-kang, L.; Zhi-yong, H. Kinematics analysis and simulation of profiling memory cutting for tunnel robot in complex environment. Procedia Earth Planet. Sci. 2009, 1, 1411–1417. [Google Scholar] [CrossRef]
  23. Li, J.; Lian, Z.; Li, Y. Kinematics analysis and carriage positioning of roboticized roadheader. J. China Coal Soc. 2008, 5, 583–587. [Google Scholar]
  24. Craig, J.J. Introduction to Robotics: Mechanics and control, 3rd ed.; China Machine Press: Beijing, China, 2006. [Google Scholar]
  25. Tian, J. Research on Boom-Type Roadheader Auto Cutting and Profiling Control System. Ph.D. Thesis, China University of Mining and Technology (Beijing), Beijing, China, 2010. [Google Scholar]
  26. Wang, S. Research on Memory Cutting and Roadway Border Control Method of Boom-Type Roadheader. Ph.D. Thesis, China University of Mining and Technology (Beijing), Beijing, China, 2014. [Google Scholar]
  27. Guo, D. Research on Profiling Cutting Theory of Longitudinal Axis Type Roadheader. Ph.D. Thesis, Liaoning Technical University, Fuxin, China, 2007. [Google Scholar]
Figure 1. Schematic diagram of secondary cutting.
Figure 1. Schematic diagram of secondary cutting.
Sustainability 15 00560 g001
Figure 2. Roadway simulation model.
Figure 2. Roadway simulation model.
Sustainability 15 00560 g002
Figure 3. Cutting path design drawing: (a) cutting path 1; (b) cutting path 2; (c) cutting path 3; (d) cutting path 4.
Figure 3. Cutting path design drawing: (a) cutting path 1; (b) cutting path 2; (c) cutting path 3; (d) cutting path 4.
Sustainability 15 00560 g003
Figure 4. Vertical displacement of roadway roof under the cutting path in cutting path 1: (a) first cut; (b) second cut; (c) third cut; (d) fourth cut; (e) fifth cut; (f) sixth cut; (g) seventh cut; (h) eighth cut; (i) ninth cut; (j) tenth cut.
Figure 4. Vertical displacement of roadway roof under the cutting path in cutting path 1: (a) first cut; (b) second cut; (c) third cut; (d) fourth cut; (e) fifth cut; (f) sixth cut; (g) seventh cut; (h) eighth cut; (i) ninth cut; (j) tenth cut.
Sustainability 15 00560 g004aSustainability 15 00560 g004b
Figure 5. Vertical displacement of roadway roof under cutting path in cutting path 2: (a) fifth cutting; (b) tenth cutting.
Figure 5. Vertical displacement of roadway roof under cutting path in cutting path 2: (a) fifth cutting; (b) tenth cutting.
Sustainability 15 00560 g005
Figure 6. Vertical displacement of the roadway roof under the cutting path in cutting path 3: (a) first cut; (b) second cut; (c) third cut; (d) fourth cut; (e) ninth cut; (f) tenth cut.
Figure 6. Vertical displacement of the roadway roof under the cutting path in cutting path 3: (a) first cut; (b) second cut; (c) third cut; (d) fourth cut; (e) ninth cut; (f) tenth cut.
Sustainability 15 00560 g006
Figure 7. Vertical displacement of roadway roof under the cutting path in cutting path 4: (a) first cut; (b) second cut; (c) third cut; (d) fourth cut; (e) ninth cut; (f) tenth cut.
Figure 7. Vertical displacement of roadway roof under the cutting path in cutting path 4: (a) first cut; (b) second cut; (c) third cut; (d) fourth cut; (e) ninth cut; (f) tenth cut.
Sustainability 15 00560 g007
Figure 8. Cutting path planning diagram.
Figure 8. Cutting path planning diagram.
Sustainability 15 00560 g008
Figure 9. Field-measured drawing of roadway cutting section.
Figure 9. Field-measured drawing of roadway cutting section.
Sustainability 15 00560 g009
Figure 10. Spatial position and attitude coordinate system of the roadheader.
Figure 10. Spatial position and attitude coordinate system of the roadheader.
Sustainability 15 00560 g010
Figure 11. Schematic diagram of position and attitude parameters of cantilever roadheader: (a) fuselage yaw angle α; (b) fuselage offset L; (c) fuselage pitch angle β; (d) fuselage rolling angle γ.
Figure 11. Schematic diagram of position and attitude parameters of cantilever roadheader: (a) fuselage yaw angle α; (b) fuselage offset L; (c) fuselage pitch angle β; (d) fuselage rolling angle γ.
Sustainability 15 00560 g011
Figure 12. Schematic diagram of secondary cutting assumption.
Figure 12. Schematic diagram of secondary cutting assumption.
Sustainability 15 00560 g012
Figure 13. Simplified model of the cutting head.
Figure 13. Simplified model of the cutting head.
Sustainability 15 00560 g013
Figure 14. Schematic diagram of cutting path inflection point.
Figure 14. Schematic diagram of cutting path inflection point.
Sustainability 15 00560 g014
Figure 15. Cutting error model diagram.
Figure 15. Cutting error model diagram.
Sustainability 15 00560 g015
Figure 16. Control procedure flow chart.
Figure 16. Control procedure flow chart.
Sustainability 15 00560 g016
Figure 17. Simulation diagram of the cutting head limit cutting area: (a) 3-dimensional diagram of a single cutting limit cutting area; (b) YOZ plane diagram of a single cutting limit cutting area; (c) 3-dimensional diagram of the secondary cutting limit cutting area; (d) YOZ plane diagram of the secondary cutting limit cutting area.
Figure 17. Simulation diagram of the cutting head limit cutting area: (a) 3-dimensional diagram of a single cutting limit cutting area; (b) YOZ plane diagram of a single cutting limit cutting area; (c) 3-dimensional diagram of the secondary cutting limit cutting area; (d) YOZ plane diagram of the secondary cutting limit cutting area.
Sustainability 15 00560 g017aSustainability 15 00560 g017b
Figure 18. Cutting path control simulation diagram.
Figure 18. Cutting path control simulation diagram.
Sustainability 15 00560 g018
Table 1. Geological parameters of the numerical simulation.
Table 1. Geological parameters of the numerical simulation.
Number
of Layers
Thickness
/m
Poisson
Ratio
Elastic Modulus
/GPa
Internal Friction
Angle/°
Tensile Strength
/Mpa
Density
kg/m3
Cohesion
/MPa
Rock
Character
13.800.2038.74342.1024009.25Siltstone
24.170.2912.65321.0224603.50Sandy mudstone
35.990.2532.74311.6224008.16Fine-grained sandstone
43.180.2912.65321.0224603.50Sandy mudstone
56.020.353.96400.5114001.60Coal
61.050.2912.65321.0224603.50Sandy mudstone
71.090.2532.74311.6224008.16Fine-grained sandstone
814.700.2038.74342.1024009.25Siltstone
Table 2. Connecting rod parameters of cutting mechanism.
Table 2. Connecting rod parameters of cutting mechanism.
Connecting Rod iai-1αi-1diθi
1a0d1θ1
2a190°0θ2
3a290°d3+d0
40d40
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

Wu, J.; Xu, Z.; Fang, X.; Shi, G.; Wang, H. Research on Path Planning and Control Method for Secondary Autonomous Cutting of Cantilever Roadheader in a Large-Section Coal Roadway. Sustainability 2023, 15, 560. https://doi.org/10.3390/su15010560

AMA Style

Wu J, Xu Z, Fang X, Shi G, Wang H. Research on Path Planning and Control Method for Secondary Autonomous Cutting of Cantilever Roadheader in a Large-Section Coal Roadway. Sustainability. 2023; 15(1):560. https://doi.org/10.3390/su15010560

Chicago/Turabian Style

Wu, Jianjun, Ziyue Xu, Xinqiu Fang, Guangliang Shi, and Haiyan Wang. 2023. "Research on Path Planning and Control Method for Secondary Autonomous Cutting of Cantilever Roadheader in a Large-Section Coal Roadway" Sustainability 15, no. 1: 560. https://doi.org/10.3390/su15010560

APA Style

Wu, J., Xu, Z., Fang, X., Shi, G., & Wang, H. (2023). Research on Path Planning and Control Method for Secondary Autonomous Cutting of Cantilever Roadheader in a Large-Section Coal Roadway. Sustainability, 15(1), 560. https://doi.org/10.3390/su15010560

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