# A Joint Positioning and Attitude Solving Method for Shearer and Scraper Conveyor under Complex Conditions

Authors: Jiacheng Xie; Zhaojian Yang; Xuewen Wang; Shuping Wang; Qing Zhang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

In a fully mechanized coal-mining face, the positioning and attitude of the shearer and scraper conveyor are inaccurate. To overcome this problem, a joint positioning and attitude solving method that considers the effect of an uneven floor is proposed. In addition, the real-time connection and coupling relationship between the two devices is analyzed. Two types of sensors, namely, the tilt sensor and strapdown inertial navigation system (SINS), are used to measure the shearer body pitch angle and the scraper conveyor shape, respectively. To improve the accuracy, two pieces of information are fused using the adaptive information fusion algorithm. It is observed that, using a marking strategy, the shearer body pitch angle can be reversely mapped to the real-time shape of the scraper conveyor. Then, a virtual-reality (VR) software that can visually simulate this entire operation process under different conditions is developed. Finally, experiments are conducted on a prototype experimental platform. The positioning error is found to be less than 0.38 times the middle trough length; moreover, no accumulated error is detected. This method can monitor the operation of the shearer and scraper conveyor in a highly dynamic and precise manner and provide strong technical support for safe and efficient operation of a fully mechanized coal-mining face.


---

## Body

## 1. Introduction

As a key equipment of a fully mechanized coal-mining face, the shearer plays a pivotal role in the coal production process. The positioning and attitude information of the shearer determines the state of the corresponding scraper conveyor [1–3] and the corresponding hydraulic supports [4]; therefore, monitoring and determination of the shearer’s running information are of significance with regard to realizing automated operation of three machines [5, 6]. The scraper conveyor mainly completes tasks such as transporting and shipping out cutting coal, running orbit of shearer, overall pushing to the coal wall side with other corresponding equipment, and coupling with the floor in real time. Therefore, the real-time shape of the scraper conveyor is the key link between the uneven floor and the shearer position and attitude information.To date, a few attempts have been made to detect the change in the topography and fluctuation range of the floor using tilt sensors installed in the shearer body [7]. Based on real-time detection results obtained for pitch and roll angles, the cutting height was compensated in the process of shearer memory cutting [7]. However, the shearer body pitch angle was the connecting line of the two supporting sliding shoes, whose distance was four-to-six times the middle trough length. It was easy to neglect the variation in coal seam between the two supporting sliding shoes; therefore, the floor described by the shearer body pitch angle was not reliable owing to the nonsensitivity character in the uneven floor [8].Wu et al. [9, 10] obtained an accurate changing trend of the floor by installing tilt sensors in every middle trough. This improved the coupling degree of the scraper conveyor and the uneven floor. In general, in the advancing direction of the working face, the height of the roof and floor changed slowly. Mutations will only appear when encountering geological structures such as faults and folds. However, with gradual changes in the roof and floor, the accumulation of circular errors led to inaccurate positioning and attitude calculation for the shearer [11, 12]. Xu [13] put forward a three-dimensional (3D) positioning theory that adopted the strategy of integrating stability and movement. On the premise of being in the same diagonal of the supporting points of the two supporting sliding shoes, Liu and Chen [14] established a digital model of the coal roof and floor. Su et al. [15] established a mathematical model of the profile cutting of the shearer. Ge [16] proposed a 3D fine geological model that enabled the shearer to adapt to complex terrain structures such as faults and folds. Based on the real-time dynamic correction strategy, Feng [17] obtained the floor curve using a shearer kinematics model.However, the above-mentioned studies were not universal to every condition owing to the idealization of the assumed conditions and nonconsideration of the coupling relationship between the scraper conveyor and shearer.Regarding the positioning of the shearer, significant research progress has been made [18–22]. This includes the development of some accurate fusion positioning methods such as SINS and encoders [23–25], wireless sensor networks [26–28], infrared cameras [29], and geographic information systems [30], which yielded remarkable results; however, these methods are still in the theoretical research stage and are not yet used in industrial applications.To overcome this barrier, Zhang et al. [31] proposed a method to detect the layout inspection of the scraper conveyor on the basis of the running trajectory and precise positioning of the shearer. This method efficiently reproduced the theoretical results. However, the shearer prototype was driven by four small wheels, and the scraper conveyor prototype was relatively simple and could not reflect the connection between the actual shearer and actual scraper conveyor efficiently. Owing to the unique and complex characteristics of the underground environment, it was difficult to conduct a physical experiment.In some laboratories, because of the heavy equipment of the fully mechanized mining system and inability to create real uneven-floor conditions, it is difficult to verify the correctness of the method. Therefore, a prototype of the shearer and scraper conveyor must be designed to simulate the underground working conditions.In this paper, to overcome the abovementioned problems, a joint positioning and attitude solving method for shearer and scraper conveyor is investigated under complex conditions.


## 2. Theoretical Analysis

### 2.1. Related Concepts of Shearer and Scraper Conveyor

The entire coordination process of three machines in a fully mechanized coal-mining face is in accordance with the shearer location and relevant regulations. As shown in Figure1, the shearer walks on the flexible scraper conveyor and cuts the coal.Figure 1
Related concepts of shearer and conveyor scraper. (1) Left walking wheel; (2) middle pin rails; (3) connecting pin rails; (4) right walking wheel; (5) right supporting sliding shoe; (6) coal plate; (7) left supporting sliding shoe; (8) middle rough.In the process of shearer running, the coal plate comes into contact with the left and right supporting sliding shoes in the middle trough and the left and right walking wheels are meshed with the pin rails of the middle trough.Therefore, there are two meshing relationships, namely, the coupling relationship between the walking wheels and shape of pin rails and that between the supporting sliding shoes and coal plate. These two relationships directly affect the shearer pitch angle; therefore, it is necessary to analyze the coupling relationship between the two groups.The pin rails are divided into two categories: the middle pin rails and the connecting pin rails. Each pin rail connects to the corresponding middle trough with two pin shafts. The middle pin rails move with the corresponding middle trough and keep the same center position of the corresponding middle trough. Meanwhile, the connecting pin rails change correspondingly according to the horizontal inclination angle.


### 2.2. Overall Design and Layout of Sensors

Some double-axis tilt sensors are installed in each middle trough to obtain the horizontal and vertical inclination angle of each middle trough in real time (Figure2).Figure 2
Connection relationship between the shearer and scraper conveyor and sensor arrangement. (1) Hinge point of the left arm and body (characteristic pointE1); (2) tilt sensors installed in the shearer body; (3) SINS device; (4) key point of the left walking wheel (characteristic pointD1); (5) key point of the right walking wheel (characteristic pointD2); (6) hinge point of the right arm and body (characteristic pointE2); (7) key point of the right supporting sliding shoe (characteristic pointO2); (8) coal plate; (9) middle pin rails; (10) connecting pin rails; (11) key point of the left supporting sliding shoe (characteristic pointO1); (12) double-axis tilt sensor installed in every middle trough.


### 2.3. Overall Framework

In this paper, the connecting and coupling relationship between the shearer and scraper conveyor is developed under complex conditions using the sensors installed in the equipment. The general method is as follows (Figure3):(1)
The horizontal and vertical inclination angles of each middle trough are obtained using double-axis tilt sensors arranged on each middle trough; therefore, the shape of the scraper conveyor can be determined.(2)
The shapes of the coal plate and the pin rails can be obtained on the basis of the analytical results obtained for a middle trough structure.(3)
By setting the left supporting sliding shoe as the positioning point of the shearer, when the shearer is at a position corresponding to a position in the scraper conveyor, the contacting mode of the two supporting sliding shoes and the coal plate is assessed and the key point coordinate of the left supporting sliding shoe is obtained.(4)
The key point coordinate of the right supporting sliding shoe is obtained using the exhaustion method.(5)
The shearer body pitch angle is solved by connecting the key points of the two supporting sliding shoes.(6)
The key points of two walking wheels are determined by the key point coordinates of the two supporting sliding shoes, and the shearer body pitch angle is determined on the basis of the left and right walking wheels and pin rails.(7)
Taking the above-mentioned point as prior knowledge, the actual real-time shearer body pitch angle obtained using tilt sensors and SINS is reversely mapped to the shape of the scraper conveyor in the actual operation process; thus, the shearer’s walking distance and position relative to the scraper conveyor can be obtained.Figure 3
Research process.


### 2.4. Positioning and Attitude Solving Method for Shearer and Scraper Conveyor

#### 2.4.1. Positioning and Attitude Solving Method for Shearer

The shearer attitude described by the several key points shown in Figure2 can be obtained using sensors installed in the shearer body. When a vertical inclination angle exists, the coordinates of these key points can be easily calculated by converting and correcting all angles.


#### 2.4.2. Positioning and Attitude Solving Method for Scraper Conveyor

Suppose that the length of the middle trough isLZBC and the horizontal and vertical inclination angles of the middle trough n are αn and βn, respectively.A piecewise function of the middle troughs in theXY plane (Figure 4) can be expressed as follows:(1)f1x=xtan⁡α1,0≤x≤x1f2x=f1x1+x-x1tan⁡α2,x1<x≤x2⋮fn-1x=fn-2xn-2+x-xn-2tan⁡αn-1,xn-2<x≤xn-1fnx=fn-1xn-1+x-xn-1tan⁡αn,xn-1<x≤xn,where xi is the boundary point of the middle trough i in the X coordinate.Figure 4
Shape of a scraper conveyor.By setting the key pointO1, which is located at the p position of the middle trough k, s/LZBC=k⋯p, where s is the shearer walking length relative to the scraper conveyor, k is the serial number of the middle trough, and p is the position of the middle trough k, (xk, yk) and (xp, yp) can be calculated. Here, (xk, yk) is the coordinate of hinge joint k of the scraper conveyor and (xp, yp) is the coordinate offset of the p position of scraper conveyor k, relative to the point (xk, yk).Therefore, if the running distance iss, the coordinates can be expressed as follows:(2)xs=xk+xp=LZBC∑i=1kcos⁡αi+LZBC∗p∗cos⁡αk+1ys=yk+yp=LZBC∑i=1ksin⁡αi+LZBC∗p∗sin⁡αk+1.


### 2.5. Analysis of Coupling Positioning and Attitude Relationship between the Shearer and Scraper Conveyor

#### 2.5.1. Coupling Relationship between Supporting Sliding Shoes and Coal Plate

(
1) Contacting Modes of Supporting Sliding Shoes and Coal Plates. The shearer body pitch angle reflects the fluctuation degree between the left and the right supporting sliding shoes. Based on a theoretical analysis, we obtained three contacting modes between the supporting sliding shoes and coal plate, as shown in Figure 5:(a)
Full contact: the base line of the supporting sliding shoe is parallel to the coal plate.(b)
Semicontact: the supporting sliding shoe is at the intersection position of the two adjacent middle troughs and can only come into contact with one middle trough.(c)
Suspending: the supporting sliding shoe is at the intersection position of the two adjacent middle troughs and cannot come into full contact with any of the two middle troughs.Figure 5
Contacting model between the supporting sliding shoes and coal plate.
(a)
(b)
(c)The determination rule of the contacting mode is shown in Table1.Table 1
Determination rule of the contacting mode.
Mode
Meaning
Condition
Calculation angle
0
Full contact, the supporting sliding shoe is fully located in a middle trough
(1)Na=Nb (2) Na ≠ Nb and FloatHA[Na] = FloatHA[Nb]
Na
10
Semicontact, in the range of the middle trough Na
(1) Na ≠ Nb and Na = No1, FloatHA[Na] > FloatHA[Nb]
Na
11
Semicontact, in the range of the middle trough Nb
(1) Na ≠ Nb and Nb = No1, FloatHA[Na] > FloatHA[Nb]
Nb
2
Suspending
(1) Na ≠ Nb and FloatHA[Na] < FloatHA[Nb]
Suspending solving algorithmPointsA, B, and O are the left, right, and middle points of the base line of the supporting sliding shoes, respectively; Na, Nb, and No are the serial numbers of the middle trough that points A, B, and C belong to, respectively; and FloatHA[i] is the horizontal inclination angle of the middle trough i.(
2) Analysis of the Contacting Mode between the Support Sliding Shoes and Coal Plate. There are three contacting modes between the supporting sliding shoes and the coal plate. Taking the semicontact case, which is the most complex condition, as an example, the shearer attitude and position parameters can be obtained using the following method. This method is known as the suspending solving algorithm and its parameters are shown in Figure 6, where XA, XB, θ1, and θ2 are unknown parameters and LH and ε are structural parameters. Among them, Na=p and Nb=p+1.Figure 6
Analysis under semicontact condition.According to the relationship, we can list the following equations:(3)XB-XA=2LHcos⁡ε∗cos⁡θ1+αpXO1-XA=LHcos⁡ε+θ1+αpXB-XC/cos⁡αp+1sin⁡θ1=2LHcos⁡εsin⁡π-αp+1-αpM1=-2∗LH∗cos⁡ε∗sin⁡αp+LH∗sin⁡ε+αp-C∗cos⁡αp+1XC-XO1M2=2∗LH∗cos⁡ε∗cos⁡αp-LH∗sin⁡ε+αpXC-XO1M3=2∗LH∗cos⁡εsin⁡αp+1-αpγ=arcsin⁡M2M12+M22,where M1,M2, and M3 are the three middle variables and γ is the middle angle.Solution:(4)θ1=π2-γXA=XO1-LHcos⁡θ1+αp+βXB=XO1+2LHcos⁡β∗cos⁡θ1+αp-LHcos⁡θ1+αp+β.SoYO1 can be expressed as follows:(5)YO1=fXA+LHsin⁡θ1+αp+βNO1=pfXA+LHsin⁡θ1+αp+1+βNO1=p+1,where, for NO1, the number of middle troughs it belongs to must be determined.(
3) Shearer Body Pitch Angle. After determining the condition of the left supporting sliding shoe, the condition of the right supporting sliding shoe must be assessed.PointO2 coordinates can be solved by the following formula:(6)XO2=XO1+Ljscos⁡αjsYO2=YO1+Ljssin⁡αjs,where αjs is the shearer body pitch angle and Ljs is the shearer body length (the connection length between point D1 and point D2).There are nine possible conditions under which the contacting mode of the two supporting sliding shoes is considered simultaneously.Owing to the difficulty in calculating the condition of the right sliding shoe using a direct method, the indirect calculation method is used, as shown in Figure7.Figure 7
Flow chart of the solving method.In Figure7, S1 is the limit position of the shearer walking on the scraper conveyor.When theXO1 coordinate increases the distance to 0.8 times the length of the shearer body, the XO2 coordinate can be analyzed and the contacting mode can be assessed. Thereby, the corresponding algorithm was used to solve the problem.Based on the condition of the distance and the shearer body length, theXO2 coordinates were assessed by comparing the XO1 coordinates. If an error exists in a small range, the solution would be correct. If an error does not fall within this range, the unit operation length would be increased to the XO2 coordinates and assessment would continue until the condition was satisfied and the correct O2 point coordinates could be solved.Therefore, the shearer body pitch angle could be calculated as follows:(7)αjs=tan⁡YO2-YO1XO2-XO1.According to the shape of the scraper conveyor, the left and right supporting sliding shoes must rotate around pointsO1 and O2, respectively; thus, they affect the shearer body pitch angle.


#### 2.5.2. Coupling Relationship between Guide Sliding Shoes and the Shape of Pin Rails

(
1) Analysis of the Shape of Pin Rails. Due to a small change in the vertical inclination angle, the connecting pin rails are bent along the shape of the two adjacent middle troughs, and their pitch angle is half the sum of the horizontal inclination angles of the two adjacent middle troughs.The horizontal inclination angle of the middle pin rails is given as follows:(8)θMi=αi.The horizontal inclination angle of the connecting pin rails is given as follows:(9)θNi=αi+αi+12.The curvilinear equation of the pin rails can be expressed according to the coordinate of each axle hole; therefore, the equation of the pin rails can be expressed as follows:(10)g1x=YMXP1+x-XMXP1∗tan⁡θM1XMXP1≤x≤XNXP1g2x=YNXP1+x-XNXP1∗tan⁡θN1XNXP1<x≤XMXP2⋮g2i-1x=YMXPi+x-XMXPi∗tan⁡θMiXMXPi≤x≤XNXPig2ix=YNXPi+x-XNXPi∗tanNi⁡XNXPi<x≤XMXPi+1,where (XMXP(i),YMXP(i)) and (XNXP(i),YNXP(i)) are the coordinates of the left and right axle holes of middle trough i, respectively.(
2) Coordinate Analysis of Walking Wheels. Coupled with the curve of the pin rails, points D1 and D2 can be calculated on the basis of points O1 and O2. The shearer body pitch angle is verified and the vertical inclination angle is adjusted until the shearer body pitch angle is equal to the value calculated in Section 2.4.1. In contrast, the vertical inclination angle must be compensated.


### 2.6. Fusion Strategy of Positioning and Attitude Based on the Information Fusion Strategy

#### 2.6.1. Information Fusion Strategy

The SINS and tilt sensors are used to measure two variables: the shearer body pitch angle and the horizontal and vertical inclination angles of every middle trough.At different temperatures and in different environments, electromagnetic interference easily affects the sensors by causing noise and error; this means that the drifting phenomenon of original data could possibly occur in a single sensor and that the true operation state of shearer and conveyor may not be accurately displayed. Thus, the information fusion algorithm was used to improve the two variables using two sensors.The theoretical values were obtained using the simulation result and the information fusion value of the middle trough obtained by two sensors, and the shearer body pitch angles were corrected and fused with the information fusion algorithm in real time.Therefore, the multisensor information fusion algorithm, which uses multiple data collected from multiple sensors at different times, marks the actual state of two devices.The premise of the adaptive algorithm is the batch algorithm, so it is necessary to explain it.The batch estimation algorithm and adaptive weighted fusion algorithm are used for calculation.(
1) Batch Estimation Algorithm. p measurement datum [γ1,γ2,…,γp], collected by one sensor at regular intervals, is divided into two groups:(1)
Whenp is odd, the two groups are [γ1,γ2,…,γ(p+1)/2] and [γ(p+1)/2,γ(p+1)/2+1,…,γp].(2)
Whenp is even, the two groups are [γ1,γ2,…,γp/2] and [γp/2+1,γp/2+2,…,γp].Taking the second case as an example, we analyzed the following.The arithmetic meanγ-1 and the mean square deviation σ1 of the first set measurements are(11)γ-1=1p/2∑i=1p/2γiσ1=1p/2-1∑i=1p/2γi-γ-1.The arithmetic meanγ-1 and the mean square deviation σ1 of the second set measurements are(12)γ-2=1p/2∑i=p/2+1pγiσ2=1p/2-1∑i=p/2+1pγi-γ-1.The batch estimationγ- and variance σ2i of the single sensor could be calculated using the following formula:(13)γ-=σ22γ-1+σ12γ-2σ12+σ22σ2=σ12•σ22σ12+σ22.The angles calculated by the above algorithms were taken as the accurate results, using which the next step was calculated and analyzed.Prior knowledge of the tilt sensor and SINS was not required, and the adaptive weighted fusion algorithm could be obtained using the value of the batch estimation angle. Working independently, every angle measured by the tilt sensor or SINS is interfered with by factors such as noise and vibration; therefore, the angle value calculated by the optimal angle is random and could be expressed as follows:(14)γm-um,σm,where um is the expected value and σm is the variance.Mutually independent of each other, the weighting factors of the tilt sensor,W1,W2,…,Wm and γ1,γ2,…,γm, are used to perform information fusion; therefore, γ, the value of integration, needs to satisfy the following relations:(15)γ=∑i=1mWiγ-i∑i=1mWii=1.The optimal weighting factor corresponding to the minimum total variance is obtained using the following formula:(16)Wi=1σi2∑i=1z1/σi2.The acquisition frequency of the sensor is determined to be 50 ms. Owing to the shearer haulage speed being generally within the range of 6–8 m/min, the walking length is small at 0.5 s. Therefore, the 10 sets of data collected on the tilt sensor and SINS (Table2) were used infusion and batch fusion, respectively; then, the accurate shearer position relative to the scraper conveyor by the adaptive weighted fusion could be obtained. In this paper, the fusion values obtained using the adaptive weighted fusion algorithm are shown in Table 3.Table 2
The shearer body pitch angle measured by tilt sensors and SINS (units: degrees).
Type
Value
Group 1
Group 2
1
2
3
4
5
6
7
8
9
10
SINS
13.52
13.61
13.63
13.67
13.53
13.49
13.52
13.67
13.69
13.63
Tilt sensor
13.69
13.70
13.90
13.84
13.84
13.69
13.71
13.82
13.86
13.87Table 3
Measured values for the SINS and tilt sensor and fusion values (units: degrees).
Tilt sensor
SINS
Group 1
Mean value
13.794
13.592
Mean square deviation
0.0088
0.0042
Group 2
Mean value
13.790
13.61
Mean square deviation
0.0071
0.0081
Batch estimation algorithm
Fusion value
13.792
13.594
Variance
5.58
e
-
5
1.39
e
-
5
Adaptive weighted fusion algorithm
Fusion value
13.664
Weighting factor
0.354
0.646In this way, a series of data is calculated as shown in Table4.Table 4
The measured value of SINS and the tilt sensor and the fusion value (units: degrees).
Number
SINS
Tilt sensor
Fusion value
(1)
13.6
13.79
13.664
(2)
13.2
13.31
13.233
(3)
14.3
14.21
14.273
(4)
12.9
12.91
12.903
(5)
12.4
13.31
12.673
(6)
5.9
6.12
5.966
(7)
8.3
8.07
8.231
(8)
4.2
4.01
4.143
(9)
5.7
5.33
5.589
(10)
4.4
4.6
4.46
(11)
0.3
0.83
0.459
(12)
1.3
1.46
1.348
(13)
8.3
8.74
8.432
(14)
1.2
1.29
1.227
(15)
−0.7
0.38
−0.376
(16)
0
0.2
0.06
(17)
−5.6
−5.63
−5.609
(18)
1.3
0.56
1.078
(19)
5.5
5.28
5.434
(20)
0.2
0.07
0.161
(21)
0.6
0.08
0.444
(22)
0.4
0.16
0.328
(23)
−8.2
−8.8
−8.38
(24)
0.7
0.15
0.535
(25)
−1
−1.42
−1.126
(26)
−0.8
−1.21
−0.923
(27)
−1.8
−1.46
−1.698
(28)
−0.2
−0.73
−0.359
(29)
0.1
−0.32
−0.026
(30)
−0.9
−1.3
−1.02
(31)
−1.3
−1.45
−1.345
(32)
−2.2
−2.6
−2.32
(33)
−6.6
−7.08
−6.744
(34)
−1.2
−1.07
−1.161
(35)
−8.4
−8.72
−8.496


#### 2.6.2. Reverse Mapping Method Based on Prior Knowledge

By considering this result obtained by the simulation as prior knowledge, the fusion value of the shearer body pitch angle, obtained using two sensors in real time, corresponds to the reverse shape of the scraper conveyor. In particular, some key inflection positions must be corrected in order to be determined according to the measured value.As shown in Figure8, the theoretical curve is first divided into some blocks corresponding to several stages, including A,B,…,M and boundary points marked as a′,b′,…,m′.Figure 8
Reverse mapping method.From prior experience, in the actual operation of the shearer, some points such asa′,b′,…,m′ are used to correct and verify the theoretical points in real time, including points a,b,…,m. Thus, every interval can be determined; then, the shearer position is reversely mapped to the shape of the scraper conveyor.


### 2.7. Planning Software Based on Unity3d

#### 2.7.1. VR Simulation Software

The models were obtained using the UG software and could access the Unity3d software through model repairing and conversion. The virtual scene was arranged according to specific rules. By integrating all the algorithms, this software could conduct various simulations under different conditions. A visual GUI was responsible for setting some simulation parameters, which included body length and structural parameters. The real-time processing datum could be output to an XML file, which was easy to analyze, as shown in Figure9.Figure 9
Interface of the Unity3d simulation software under complex conditions.The shape of the scraper conveyor could be estimated using the input parameters of every middle trough. To coordinate the virtual shape of the scraper conveyor, the walking position and the attitude of the shearer were calculated backstage in real time. The shearer position was calculated according to the virtual shearer haulage speed, which was decided by an increment, and the shearer position was reversely mapped to the shape of the scraper conveyor.


## 2.1. Related Concepts of Shearer and Scraper Conveyor

The entire coordination process of three machines in a fully mechanized coal-mining face is in accordance with the shearer location and relevant regulations. As shown in Figure1, the shearer walks on the flexible scraper conveyor and cuts the coal.Figure 1
Related concepts of shearer and conveyor scraper. (1) Left walking wheel; (2) middle pin rails; (3) connecting pin rails; (4) right walking wheel; (5) right supporting sliding shoe; (6) coal plate; (7) left supporting sliding shoe; (8) middle rough.In the process of shearer running, the coal plate comes into contact with the left and right supporting sliding shoes in the middle trough and the left and right walking wheels are meshed with the pin rails of the middle trough.Therefore, there are two meshing relationships, namely, the coupling relationship between the walking wheels and shape of pin rails and that between the supporting sliding shoes and coal plate. These two relationships directly affect the shearer pitch angle; therefore, it is necessary to analyze the coupling relationship between the two groups.The pin rails are divided into two categories: the middle pin rails and the connecting pin rails. Each pin rail connects to the corresponding middle trough with two pin shafts. The middle pin rails move with the corresponding middle trough and keep the same center position of the corresponding middle trough. Meanwhile, the connecting pin rails change correspondingly according to the horizontal inclination angle.


## 2.2. Overall Design and Layout of Sensors

Some double-axis tilt sensors are installed in each middle trough to obtain the horizontal and vertical inclination angle of each middle trough in real time (Figure2).Figure 2
Connection relationship between the shearer and scraper conveyor and sensor arrangement. (1) Hinge point of the left arm and body (characteristic pointE1); (2) tilt sensors installed in the shearer body; (3) SINS device; (4) key point of the left walking wheel (characteristic pointD1); (5) key point of the right walking wheel (characteristic pointD2); (6) hinge point of the right arm and body (characteristic pointE2); (7) key point of the right supporting sliding shoe (characteristic pointO2); (8) coal plate; (9) middle pin rails; (10) connecting pin rails; (11) key point of the left supporting sliding shoe (characteristic pointO1); (12) double-axis tilt sensor installed in every middle trough.


## 2.3. Overall Framework

In this paper, the connecting and coupling relationship between the shearer and scraper conveyor is developed under complex conditions using the sensors installed in the equipment. The general method is as follows (Figure3):(1)
The horizontal and vertical inclination angles of each middle trough are obtained using double-axis tilt sensors arranged on each middle trough; therefore, the shape of the scraper conveyor can be determined.(2)
The shapes of the coal plate and the pin rails can be obtained on the basis of the analytical results obtained for a middle trough structure.(3)
By setting the left supporting sliding shoe as the positioning point of the shearer, when the shearer is at a position corresponding to a position in the scraper conveyor, the contacting mode of the two supporting sliding shoes and the coal plate is assessed and the key point coordinate of the left supporting sliding shoe is obtained.(4)
The key point coordinate of the right supporting sliding shoe is obtained using the exhaustion method.(5)
The shearer body pitch angle is solved by connecting the key points of the two supporting sliding shoes.(6)
The key points of two walking wheels are determined by the key point coordinates of the two supporting sliding shoes, and the shearer body pitch angle is determined on the basis of the left and right walking wheels and pin rails.(7)
Taking the above-mentioned point as prior knowledge, the actual real-time shearer body pitch angle obtained using tilt sensors and SINS is reversely mapped to the shape of the scraper conveyor in the actual operation process; thus, the shearer’s walking distance and position relative to the scraper conveyor can be obtained.Figure 3
Research process.


## 2.4. Positioning and Attitude Solving Method for Shearer and Scraper Conveyor

### 2.4.1. Positioning and Attitude Solving Method for Shearer

The shearer attitude described by the several key points shown in Figure2 can be obtained using sensors installed in the shearer body. When a vertical inclination angle exists, the coordinates of these key points can be easily calculated by converting and correcting all angles.


### 2.4.2. Positioning and Attitude Solving Method for Scraper Conveyor

Suppose that the length of the middle trough isLZBC and the horizontal and vertical inclination angles of the middle trough n are αn and βn, respectively.A piecewise function of the middle troughs in theXY plane (Figure 4) can be expressed as follows:(1)f1x=xtan⁡α1,0≤x≤x1f2x=f1x1+x-x1tan⁡α2,x1<x≤x2⋮fn-1x=fn-2xn-2+x-xn-2tan⁡αn-1,xn-2<x≤xn-1fnx=fn-1xn-1+x-xn-1tan⁡αn,xn-1<x≤xn,where xi is the boundary point of the middle trough i in the X coordinate.Figure 4
Shape of a scraper conveyor.By setting the key pointO1, which is located at the p position of the middle trough k, s/LZBC=k⋯p, where s is the shearer walking length relative to the scraper conveyor, k is the serial number of the middle trough, and p is the position of the middle trough k, (xk, yk) and (xp, yp) can be calculated. Here, (xk, yk) is the coordinate of hinge joint k of the scraper conveyor and (xp, yp) is the coordinate offset of the p position of scraper conveyor k, relative to the point (xk, yk).Therefore, if the running distance iss, the coordinates can be expressed as follows:(2)xs=xk+xp=LZBC∑i=1kcos⁡αi+LZBC∗p∗cos⁡αk+1ys=yk+yp=LZBC∑i=1ksin⁡αi+LZBC∗p∗sin⁡αk+1.


## 2.4.1. Positioning and Attitude Solving Method for Shearer

The shearer attitude described by the several key points shown in Figure2 can be obtained using sensors installed in the shearer body. When a vertical inclination angle exists, the coordinates of these key points can be easily calculated by converting and correcting all angles.


## 2.4.2. Positioning and Attitude Solving Method for Scraper Conveyor

Suppose that the length of the middle trough isLZBC and the horizontal and vertical inclination angles of the middle trough n are αn and βn, respectively.A piecewise function of the middle troughs in theXY plane (Figure 4) can be expressed as follows:(1)f1x=xtan⁡α1,0≤x≤x1f2x=f1x1+x-x1tan⁡α2,x1<x≤x2⋮fn-1x=fn-2xn-2+x-xn-2tan⁡αn-1,xn-2<x≤xn-1fnx=fn-1xn-1+x-xn-1tan⁡αn,xn-1<x≤xn,where xi is the boundary point of the middle trough i in the X coordinate.Figure 4
Shape of a scraper conveyor.By setting the key pointO1, which is located at the p position of the middle trough k, s/LZBC=k⋯p, where s is the shearer walking length relative to the scraper conveyor, k is the serial number of the middle trough, and p is the position of the middle trough k, (xk, yk) and (xp, yp) can be calculated. Here, (xk, yk) is the coordinate of hinge joint k of the scraper conveyor and (xp, yp) is the coordinate offset of the p position of scraper conveyor k, relative to the point (xk, yk).Therefore, if the running distance iss, the coordinates can be expressed as follows:(2)xs=xk+xp=LZBC∑i=1kcos⁡αi+LZBC∗p∗cos⁡αk+1ys=yk+yp=LZBC∑i=1ksin⁡αi+LZBC∗p∗sin⁡αk+1.


## 2.5. Analysis of Coupling Positioning and Attitude Relationship between the Shearer and Scraper Conveyor

### 2.5.1. Coupling Relationship between Supporting Sliding Shoes and Coal Plate

(
1) Contacting Modes of Supporting Sliding Shoes and Coal Plates. The shearer body pitch angle reflects the fluctuation degree between the left and the right supporting sliding shoes. Based on a theoretical analysis, we obtained three contacting modes between the supporting sliding shoes and coal plate, as shown in Figure 5:(a)
Full contact: the base line of the supporting sliding shoe is parallel to the coal plate.(b)
Semicontact: the supporting sliding shoe is at the intersection position of the two adjacent middle troughs and can only come into contact with one middle trough.(c)
Suspending: the supporting sliding shoe is at the intersection position of the two adjacent middle troughs and cannot come into full contact with any of the two middle troughs.Figure 5
Contacting model between the supporting sliding shoes and coal plate.
(a)
(b)
(c)The determination rule of the contacting mode is shown in Table1.Table 1
Determination rule of the contacting mode.
Mode
Meaning
Condition
Calculation angle
0
Full contact, the supporting sliding shoe is fully located in a middle trough
(1)Na=Nb (2) Na ≠ Nb and FloatHA[Na] = FloatHA[Nb]
Na
10
Semicontact, in the range of the middle trough Na
(1) Na ≠ Nb and Na = No1, FloatHA[Na] > FloatHA[Nb]
Na
11
Semicontact, in the range of the middle trough Nb
(1) Na ≠ Nb and Nb = No1, FloatHA[Na] > FloatHA[Nb]
Nb
2
Suspending
(1) Na ≠ Nb and FloatHA[Na] < FloatHA[Nb]
Suspending solving algorithmPointsA, B, and O are the left, right, and middle points of the base line of the supporting sliding shoes, respectively; Na, Nb, and No are the serial numbers of the middle trough that points A, B, and C belong to, respectively; and FloatHA[i] is the horizontal inclination angle of the middle trough i.(
2) Analysis of the Contacting Mode between the Support Sliding Shoes and Coal Plate. There are three contacting modes between the supporting sliding shoes and the coal plate. Taking the semicontact case, which is the most complex condition, as an example, the shearer attitude and position parameters can be obtained using the following method. This method is known as the suspending solving algorithm and its parameters are shown in Figure 6, where XA, XB, θ1, and θ2 are unknown parameters and LH and ε are structural parameters. Among them, Na=p and Nb=p+1.Figure 6
Analysis under semicontact condition.According to the relationship, we can list the following equations:(3)XB-XA=2LHcos⁡ε∗cos⁡θ1+αpXO1-XA=LHcos⁡ε+θ1+αpXB-XC/cos⁡αp+1sin⁡θ1=2LHcos⁡εsin⁡π-αp+1-αpM1=-2∗LH∗cos⁡ε∗sin⁡αp+LH∗sin⁡ε+αp-C∗cos⁡αp+1XC-XO1M2=2∗LH∗cos⁡ε∗cos⁡αp-LH∗sin⁡ε+αpXC-XO1M3=2∗LH∗cos⁡εsin⁡αp+1-αpγ=arcsin⁡M2M12+M22,where M1,M2, and M3 are the three middle variables and γ is the middle angle.Solution:(4)θ1=π2-γXA=XO1-LHcos⁡θ1+αp+βXB=XO1+2LHcos⁡β∗cos⁡θ1+αp-LHcos⁡θ1+αp+β.SoYO1 can be expressed as follows:(5)YO1=fXA+LHsin⁡θ1+αp+βNO1=pfXA+LHsin⁡θ1+αp+1+βNO1=p+1,where, for NO1, the number of middle troughs it belongs to must be determined.(
3) Shearer Body Pitch Angle. After determining the condition of the left supporting sliding shoe, the condition of the right supporting sliding shoe must be assessed.PointO2 coordinates can be solved by the following formula:(6)XO2=XO1+Ljscos⁡αjsYO2=YO1+Ljssin⁡αjs,where αjs is the shearer body pitch angle and Ljs is the shearer body length (the connection length between point D1 and point D2).There are nine possible conditions under which the contacting mode of the two supporting sliding shoes is considered simultaneously.Owing to the difficulty in calculating the condition of the right sliding shoe using a direct method, the indirect calculation method is used, as shown in Figure7.Figure 7
Flow chart of the solving method.In Figure7, S1 is the limit position of the shearer walking on the scraper conveyor.When theXO1 coordinate increases the distance to 0.8 times the length of the shearer body, the XO2 coordinate can be analyzed and the contacting mode can be assessed. Thereby, the corresponding algorithm was used to solve the problem.Based on the condition of the distance and the shearer body length, theXO2 coordinates were assessed by comparing the XO1 coordinates. If an error exists in a small range, the solution would be correct. If an error does not fall within this range, the unit operation length would be increased to the XO2 coordinates and assessment would continue until the condition was satisfied and the correct O2 point coordinates could be solved.Therefore, the shearer body pitch angle could be calculated as follows:(7)αjs=tan⁡YO2-YO1XO2-XO1.According to the shape of the scraper conveyor, the left and right supporting sliding shoes must rotate around pointsO1 and O2, respectively; thus, they affect the shearer body pitch angle.


### 2.5.2. Coupling Relationship between Guide Sliding Shoes and the Shape of Pin Rails

(
1) Analysis of the Shape of Pin Rails. Due to a small change in the vertical inclination angle, the connecting pin rails are bent along the shape of the two adjacent middle troughs, and their pitch angle is half the sum of the horizontal inclination angles of the two adjacent middle troughs.The horizontal inclination angle of the middle pin rails is given as follows:(8)θMi=αi.The horizontal inclination angle of the connecting pin rails is given as follows:(9)θNi=αi+αi+12.The curvilinear equation of the pin rails can be expressed according to the coordinate of each axle hole; therefore, the equation of the pin rails can be expressed as follows:(10)g1x=YMXP1+x-XMXP1∗tan⁡θM1XMXP1≤x≤XNXP1g2x=YNXP1+x-XNXP1∗tan⁡θN1XNXP1<x≤XMXP2⋮g2i-1x=YMXPi+x-XMXPi∗tan⁡θMiXMXPi≤x≤XNXPig2ix=YNXPi+x-XNXPi∗tanNi⁡XNXPi<x≤XMXPi+1,where (XMXP(i),YMXP(i)) and (XNXP(i),YNXP(i)) are the coordinates of the left and right axle holes of middle trough i, respectively.(
2) Coordinate Analysis of Walking Wheels. Coupled with the curve of the pin rails, points D1 and D2 can be calculated on the basis of points O1 and O2. The shearer body pitch angle is verified and the vertical inclination angle is adjusted until the shearer body pitch angle is equal to the value calculated in Section 2.4.1. In contrast, the vertical inclination angle must be compensated.


## 2.5.1. Coupling Relationship between Supporting Sliding Shoes and Coal Plate

(
1) Contacting Modes of Supporting Sliding Shoes and Coal Plates. The shearer body pitch angle reflects the fluctuation degree between the left and the right supporting sliding shoes. Based on a theoretical analysis, we obtained three contacting modes between the supporting sliding shoes and coal plate, as shown in Figure 5:(a)
Full contact: the base line of the supporting sliding shoe is parallel to the coal plate.(b)
Semicontact: the supporting sliding shoe is at the intersection position of the two adjacent middle troughs and can only come into contact with one middle trough.(c)
Suspending: the supporting sliding shoe is at the intersection position of the two adjacent middle troughs and cannot come into full contact with any of the two middle troughs.Figure 5
Contacting model between the supporting sliding shoes and coal plate.
(a)
(b)
(c)The determination rule of the contacting mode is shown in Table1.Table 1
Determination rule of the contacting mode.
Mode
Meaning
Condition
Calculation angle
0
Full contact, the supporting sliding shoe is fully located in a middle trough
(1)Na=Nb (2) Na ≠ Nb and FloatHA[Na] = FloatHA[Nb]
Na
10
Semicontact, in the range of the middle trough Na
(1) Na ≠ Nb and Na = No1, FloatHA[Na] > FloatHA[Nb]
Na
11
Semicontact, in the range of the middle trough Nb
(1) Na ≠ Nb and Nb = No1, FloatHA[Na] > FloatHA[Nb]
Nb
2
Suspending
(1) Na ≠ Nb and FloatHA[Na] < FloatHA[Nb]
Suspending solving algorithmPointsA, B, and O are the left, right, and middle points of the base line of the supporting sliding shoes, respectively; Na, Nb, and No are the serial numbers of the middle trough that points A, B, and C belong to, respectively; and FloatHA[i] is the horizontal inclination angle of the middle trough i.(
2) Analysis of the Contacting Mode between the Support Sliding Shoes and Coal Plate. There are three contacting modes between the supporting sliding shoes and the coal plate. Taking the semicontact case, which is the most complex condition, as an example, the shearer attitude and position parameters can be obtained using the following method. This method is known as the suspending solving algorithm and its parameters are shown in Figure 6, where XA, XB, θ1, and θ2 are unknown parameters and LH and ε are structural parameters. Among them, Na=p and Nb=p+1.Figure 6
Analysis under semicontact condition.According to the relationship, we can list the following equations:(3)XB-XA=2LHcos⁡ε∗cos⁡θ1+αpXO1-XA=LHcos⁡ε+θ1+αpXB-XC/cos⁡αp+1sin⁡θ1=2LHcos⁡εsin⁡π-αp+1-αpM1=-2∗LH∗cos⁡ε∗sin⁡αp+LH∗sin⁡ε+αp-C∗cos⁡αp+1XC-XO1M2=2∗LH∗cos⁡ε∗cos⁡αp-LH∗sin⁡ε+αpXC-XO1M3=2∗LH∗cos⁡εsin⁡αp+1-αpγ=arcsin⁡M2M12+M22,where M1,M2, and M3 are the three middle variables and γ is the middle angle.Solution:(4)θ1=π2-γXA=XO1-LHcos⁡θ1+αp+βXB=XO1+2LHcos⁡β∗cos⁡θ1+αp-LHcos⁡θ1+αp+β.SoYO1 can be expressed as follows:(5)YO1=fXA+LHsin⁡θ1+αp+βNO1=pfXA+LHsin⁡θ1+αp+1+βNO1=p+1,where, for NO1, the number of middle troughs it belongs to must be determined.(
3) Shearer Body Pitch Angle. After determining the condition of the left supporting sliding shoe, the condition of the right supporting sliding shoe must be assessed.PointO2 coordinates can be solved by the following formula:(6)XO2=XO1+Ljscos⁡αjsYO2=YO1+Ljssin⁡αjs,where αjs is the shearer body pitch angle and Ljs is the shearer body length (the connection length between point D1 and point D2).There are nine possible conditions under which the contacting mode of the two supporting sliding shoes is considered simultaneously.Owing to the difficulty in calculating the condition of the right sliding shoe using a direct method, the indirect calculation method is used, as shown in Figure7.Figure 7
Flow chart of the solving method.In Figure7, S1 is the limit position of the shearer walking on the scraper conveyor.When theXO1 coordinate increases the distance to 0.8 times the length of the shearer body, the XO2 coordinate can be analyzed and the contacting mode can be assessed. Thereby, the corresponding algorithm was used to solve the problem.Based on the condition of the distance and the shearer body length, theXO2 coordinates were assessed by comparing the XO1 coordinates. If an error exists in a small range, the solution would be correct. If an error does not fall within this range, the unit operation length would be increased to the XO2 coordinates and assessment would continue until the condition was satisfied and the correct O2 point coordinates could be solved.Therefore, the shearer body pitch angle could be calculated as follows:(7)αjs=tan⁡YO2-YO1XO2-XO1.According to the shape of the scraper conveyor, the left and right supporting sliding shoes must rotate around pointsO1 and O2, respectively; thus, they affect the shearer body pitch angle.


## 2.5.2. Coupling Relationship between Guide Sliding Shoes and the Shape of Pin Rails

(
1) Analysis of the Shape of Pin Rails. Due to a small change in the vertical inclination angle, the connecting pin rails are bent along the shape of the two adjacent middle troughs, and their pitch angle is half the sum of the horizontal inclination angles of the two adjacent middle troughs.The horizontal inclination angle of the middle pin rails is given as follows:(8)θMi=αi.The horizontal inclination angle of the connecting pin rails is given as follows:(9)θNi=αi+αi+12.The curvilinear equation of the pin rails can be expressed according to the coordinate of each axle hole; therefore, the equation of the pin rails can be expressed as follows:(10)g1x=YMXP1+x-XMXP1∗tan⁡θM1XMXP1≤x≤XNXP1g2x=YNXP1+x-XNXP1∗tan⁡θN1XNXP1<x≤XMXP2⋮g2i-1x=YMXPi+x-XMXPi∗tan⁡θMiXMXPi≤x≤XNXPig2ix=YNXPi+x-XNXPi∗tanNi⁡XNXPi<x≤XMXPi+1,where (XMXP(i),YMXP(i)) and (XNXP(i),YNXP(i)) are the coordinates of the left and right axle holes of middle trough i, respectively.(
2) Coordinate Analysis of Walking Wheels. Coupled with the curve of the pin rails, points D1 and D2 can be calculated on the basis of points O1 and O2. The shearer body pitch angle is verified and the vertical inclination angle is adjusted until the shearer body pitch angle is equal to the value calculated in Section 2.4.1. In contrast, the vertical inclination angle must be compensated.


## 2.6. Fusion Strategy of Positioning and Attitude Based on the Information Fusion Strategy

### 2.6.1. Information Fusion Strategy

The SINS and tilt sensors are used to measure two variables: the shearer body pitch angle and the horizontal and vertical inclination angles of every middle trough.At different temperatures and in different environments, electromagnetic interference easily affects the sensors by causing noise and error; this means that the drifting phenomenon of original data could possibly occur in a single sensor and that the true operation state of shearer and conveyor may not be accurately displayed. Thus, the information fusion algorithm was used to improve the two variables using two sensors.The theoretical values were obtained using the simulation result and the information fusion value of the middle trough obtained by two sensors, and the shearer body pitch angles were corrected and fused with the information fusion algorithm in real time.Therefore, the multisensor information fusion algorithm, which uses multiple data collected from multiple sensors at different times, marks the actual state of two devices.The premise of the adaptive algorithm is the batch algorithm, so it is necessary to explain it.The batch estimation algorithm and adaptive weighted fusion algorithm are used for calculation.(
1) Batch Estimation Algorithm. p measurement datum [γ1,γ2,…,γp], collected by one sensor at regular intervals, is divided into two groups:(1)
Whenp is odd, the two groups are [γ1,γ2,…,γ(p+1)/2] and [γ(p+1)/2,γ(p+1)/2+1,…,γp].(2)
Whenp is even, the two groups are [γ1,γ2,…,γp/2] and [γp/2+1,γp/2+2,…,γp].Taking the second case as an example, we analyzed the following.The arithmetic meanγ-1 and the mean square deviation σ1 of the first set measurements are(11)γ-1=1p/2∑i=1p/2γiσ1=1p/2-1∑i=1p/2γi-γ-1.The arithmetic meanγ-1 and the mean square deviation σ1 of the second set measurements are(12)γ-2=1p/2∑i=p/2+1pγiσ2=1p/2-1∑i=p/2+1pγi-γ-1.The batch estimationγ- and variance σ2i of the single sensor could be calculated using the following formula:(13)γ-=σ22γ-1+σ12γ-2σ12+σ22σ2=σ12•σ22σ12+σ22.The angles calculated by the above algorithms were taken as the accurate results, using which the next step was calculated and analyzed.Prior knowledge of the tilt sensor and SINS was not required, and the adaptive weighted fusion algorithm could be obtained using the value of the batch estimation angle. Working independently, every angle measured by the tilt sensor or SINS is interfered with by factors such as noise and vibration; therefore, the angle value calculated by the optimal angle is random and could be expressed as follows:(14)γm-um,σm,where um is the expected value and σm is the variance.Mutually independent of each other, the weighting factors of the tilt sensor,W1,W2,…,Wm and γ1,γ2,…,γm, are used to perform information fusion; therefore, γ, the value of integration, needs to satisfy the following relations:(15)γ=∑i=1mWiγ-i∑i=1mWii=1.The optimal weighting factor corresponding to the minimum total variance is obtained using the following formula:(16)Wi=1σi2∑i=1z1/σi2.The acquisition frequency of the sensor is determined to be 50 ms. Owing to the shearer haulage speed being generally within the range of 6–8 m/min, the walking length is small at 0.5 s. Therefore, the 10 sets of data collected on the tilt sensor and SINS (Table2) were used infusion and batch fusion, respectively; then, the accurate shearer position relative to the scraper conveyor by the adaptive weighted fusion could be obtained. In this paper, the fusion values obtained using the adaptive weighted fusion algorithm are shown in Table 3.Table 2
The shearer body pitch angle measured by tilt sensors and SINS (units: degrees).
Type
Value
Group 1
Group 2
1
2
3
4
5
6
7
8
9
10
SINS
13.52
13.61
13.63
13.67
13.53
13.49
13.52
13.67
13.69
13.63
Tilt sensor
13.69
13.70
13.90
13.84
13.84
13.69
13.71
13.82
13.86
13.87Table 3
Measured values for the SINS and tilt sensor and fusion values (units: degrees).
Tilt sensor
SINS
Group 1
Mean value
13.794
13.592
Mean square deviation
0.0088
0.0042
Group 2
Mean value
13.790
13.61
Mean square deviation
0.0071
0.0081
Batch estimation algorithm
Fusion value
13.792
13.594
Variance
5.58
e
-
5
1.39
e
-
5
Adaptive weighted fusion algorithm
Fusion value
13.664
Weighting factor
0.354
0.646In this way, a series of data is calculated as shown in Table4.Table 4
The measured value of SINS and the tilt sensor and the fusion value (units: degrees).
Number
SINS
Tilt sensor
Fusion value
(1)
13.6
13.79
13.664
(2)
13.2
13.31
13.233
(3)
14.3
14.21
14.273
(4)
12.9
12.91
12.903
(5)
12.4
13.31
12.673
(6)
5.9
6.12
5.966
(7)
8.3
8.07
8.231
(8)
4.2
4.01
4.143
(9)
5.7
5.33
5.589
(10)
4.4
4.6
4.46
(11)
0.3
0.83
0.459
(12)
1.3
1.46
1.348
(13)
8.3
8.74
8.432
(14)
1.2
1.29
1.227
(15)
−0.7
0.38
−0.376
(16)
0
0.2
0.06
(17)
−5.6
−5.63
−5.609
(18)
1.3
0.56
1.078
(19)
5.5
5.28
5.434
(20)
0.2
0.07
0.161
(21)
0.6
0.08
0.444
(22)
0.4
0.16
0.328
(23)
−8.2
−8.8
−8.38
(24)
0.7
0.15
0.535
(25)
−1
−1.42
−1.126
(26)
−0.8
−1.21
−0.923
(27)
−1.8
−1.46
−1.698
(28)
−0.2
−0.73
−0.359
(29)
0.1
−0.32
−0.026
(30)
−0.9
−1.3
−1.02
(31)
−1.3
−1.45
−1.345
(32)
−2.2
−2.6
−2.32
(33)
−6.6
−7.08
−6.744
(34)
−1.2
−1.07
−1.161
(35)
−8.4
−8.72
−8.496


### 2.6.2. Reverse Mapping Method Based on Prior Knowledge

By considering this result obtained by the simulation as prior knowledge, the fusion value of the shearer body pitch angle, obtained using two sensors in real time, corresponds to the reverse shape of the scraper conveyor. In particular, some key inflection positions must be corrected in order to be determined according to the measured value.As shown in Figure8, the theoretical curve is first divided into some blocks corresponding to several stages, including A,B,…,M and boundary points marked as a′,b′,…,m′.Figure 8
Reverse mapping method.From prior experience, in the actual operation of the shearer, some points such asa′,b′,…,m′ are used to correct and verify the theoretical points in real time, including points a,b,…,m. Thus, every interval can be determined; then, the shearer position is reversely mapped to the shape of the scraper conveyor.


## 2.6.1. Information Fusion Strategy

The SINS and tilt sensors are used to measure two variables: the shearer body pitch angle and the horizontal and vertical inclination angles of every middle trough.At different temperatures and in different environments, electromagnetic interference easily affects the sensors by causing noise and error; this means that the drifting phenomenon of original data could possibly occur in a single sensor and that the true operation state of shearer and conveyor may not be accurately displayed. Thus, the information fusion algorithm was used to improve the two variables using two sensors.The theoretical values were obtained using the simulation result and the information fusion value of the middle trough obtained by two sensors, and the shearer body pitch angles were corrected and fused with the information fusion algorithm in real time.Therefore, the multisensor information fusion algorithm, which uses multiple data collected from multiple sensors at different times, marks the actual state of two devices.The premise of the adaptive algorithm is the batch algorithm, so it is necessary to explain it.The batch estimation algorithm and adaptive weighted fusion algorithm are used for calculation.(
1) Batch Estimation Algorithm. p measurement datum [γ1,γ2,…,γp], collected by one sensor at regular intervals, is divided into two groups:(1)
Whenp is odd, the two groups are [γ1,γ2,…,γ(p+1)/2] and [γ(p+1)/2,γ(p+1)/2+1,…,γp].(2)
Whenp is even, the two groups are [γ1,γ2,…,γp/2] and [γp/2+1,γp/2+2,…,γp].Taking the second case as an example, we analyzed the following.The arithmetic meanγ-1 and the mean square deviation σ1 of the first set measurements are(11)γ-1=1p/2∑i=1p/2γiσ1=1p/2-1∑i=1p/2γi-γ-1.The arithmetic meanγ-1 and the mean square deviation σ1 of the second set measurements are(12)γ-2=1p/2∑i=p/2+1pγiσ2=1p/2-1∑i=p/2+1pγi-γ-1.The batch estimationγ- and variance σ2i of the single sensor could be calculated using the following formula:(13)γ-=σ22γ-1+σ12γ-2σ12+σ22σ2=σ12•σ22σ12+σ22.The angles calculated by the above algorithms were taken as the accurate results, using which the next step was calculated and analyzed.Prior knowledge of the tilt sensor and SINS was not required, and the adaptive weighted fusion algorithm could be obtained using the value of the batch estimation angle. Working independently, every angle measured by the tilt sensor or SINS is interfered with by factors such as noise and vibration; therefore, the angle value calculated by the optimal angle is random and could be expressed as follows:(14)γm-um,σm,where um is the expected value and σm is the variance.Mutually independent of each other, the weighting factors of the tilt sensor,W1,W2,…,Wm and γ1,γ2,…,γm, are used to perform information fusion; therefore, γ, the value of integration, needs to satisfy the following relations:(15)γ=∑i=1mWiγ-i∑i=1mWii=1.The optimal weighting factor corresponding to the minimum total variance is obtained using the following formula:(16)Wi=1σi2∑i=1z1/σi2.The acquisition frequency of the sensor is determined to be 50 ms. Owing to the shearer haulage speed being generally within the range of 6–8 m/min, the walking length is small at 0.5 s. Therefore, the 10 sets of data collected on the tilt sensor and SINS (Table2) were used infusion and batch fusion, respectively; then, the accurate shearer position relative to the scraper conveyor by the adaptive weighted fusion could be obtained. In this paper, the fusion values obtained using the adaptive weighted fusion algorithm are shown in Table 3.Table 2
The shearer body pitch angle measured by tilt sensors and SINS (units: degrees).
Type
Value
Group 1
Group 2
1
2
3
4
5
6
7
8
9
10
SINS
13.52
13.61
13.63
13.67
13.53
13.49
13.52
13.67
13.69
13.63
Tilt sensor
13.69
13.70
13.90
13.84
13.84
13.69
13.71
13.82
13.86
13.87Table 3
Measured values for the SINS and tilt sensor and fusion values (units: degrees).
Tilt sensor
SINS
Group 1
Mean value
13.794
13.592
Mean square deviation
0.0088
0.0042
Group 2
Mean value
13.790
13.61
Mean square deviation
0.0071
0.0081
Batch estimation algorithm
Fusion value
13.792
13.594
Variance
5.58
e
-
5
1.39
e
-
5
Adaptive weighted fusion algorithm
Fusion value
13.664
Weighting factor
0.354
0.646In this way, a series of data is calculated as shown in Table4.Table 4
The measured value of SINS and the tilt sensor and the fusion value (units: degrees).
Number
SINS
Tilt sensor
Fusion value
(1)
13.6
13.79
13.664
(2)
13.2
13.31
13.233
(3)
14.3
14.21
14.273
(4)
12.9
12.91
12.903
(5)
12.4
13.31
12.673
(6)
5.9
6.12
5.966
(7)
8.3
8.07
8.231
(8)
4.2
4.01
4.143
(9)
5.7
5.33
5.589
(10)
4.4
4.6
4.46
(11)
0.3
0.83
0.459
(12)
1.3
1.46
1.348
(13)
8.3
8.74
8.432
(14)
1.2
1.29
1.227
(15)
−0.7
0.38
−0.376
(16)
0
0.2
0.06
(17)
−5.6
−5.63
−5.609
(18)
1.3
0.56
1.078
(19)
5.5
5.28
5.434
(20)
0.2
0.07
0.161
(21)
0.6
0.08
0.444
(22)
0.4
0.16
0.328
(23)
−8.2
−8.8
−8.38
(24)
0.7
0.15
0.535
(25)
−1
−1.42
−1.126
(26)
−0.8
−1.21
−0.923
(27)
−1.8
−1.46
−1.698
(28)
−0.2
−0.73
−0.359
(29)
0.1
−0.32
−0.026
(30)
−0.9
−1.3
−1.02
(31)
−1.3
−1.45
−1.345
(32)
−2.2
−2.6
−2.32
(33)
−6.6
−7.08
−6.744
(34)
−1.2
−1.07
−1.161
(35)
−8.4
−8.72
−8.496


## 2.6.2. Reverse Mapping Method Based on Prior Knowledge

By considering this result obtained by the simulation as prior knowledge, the fusion value of the shearer body pitch angle, obtained using two sensors in real time, corresponds to the reverse shape of the scraper conveyor. In particular, some key inflection positions must be corrected in order to be determined according to the measured value.As shown in Figure8, the theoretical curve is first divided into some blocks corresponding to several stages, including A,B,…,M and boundary points marked as a′,b′,…,m′.Figure 8
Reverse mapping method.From prior experience, in the actual operation of the shearer, some points such asa′,b′,…,m′ are used to correct and verify the theoretical points in real time, including points a,b,…,m. Thus, every interval can be determined; then, the shearer position is reversely mapped to the shape of the scraper conveyor.


## 2.7. Planning Software Based on Unity3d

### 2.7.1. VR Simulation Software

The models were obtained using the UG software and could access the Unity3d software through model repairing and conversion. The virtual scene was arranged according to specific rules. By integrating all the algorithms, this software could conduct various simulations under different conditions. A visual GUI was responsible for setting some simulation parameters, which included body length and structural parameters. The real-time processing datum could be output to an XML file, which was easy to analyze, as shown in Figure9.Figure 9
Interface of the Unity3d simulation software under complex conditions.The shape of the scraper conveyor could be estimated using the input parameters of every middle trough. To coordinate the virtual shape of the scraper conveyor, the walking position and the attitude of the shearer were calculated backstage in real time. The shearer position was calculated according to the virtual shearer haulage speed, which was decided by an increment, and the shearer position was reversely mapped to the shape of the scraper conveyor.


## 2.7.1. VR Simulation Software

The models were obtained using the UG software and could access the Unity3d software through model repairing and conversion. The virtual scene was arranged according to specific rules. By integrating all the algorithms, this software could conduct various simulations under different conditions. A visual GUI was responsible for setting some simulation parameters, which included body length and structural parameters. The real-time processing datum could be output to an XML file, which was easy to analyze, as shown in Figure9.Figure 9
Interface of the Unity3d simulation software under complex conditions.The shape of the scraper conveyor could be estimated using the input parameters of every middle trough. To coordinate the virtual shape of the scraper conveyor, the walking position and the attitude of the shearer were calculated backstage in real time. The shearer position was calculated according to the virtual shearer haulage speed, which was decided by an increment, and the shearer position was reversely mapped to the shape of the scraper conveyor.


## 3. Experiments and Results

### 3.1. Test Prototype

Three machines in our laboratory were selected as the research objects. The type of the scraper conveyor was SGZ764/630. The type of the shearer was MGTY250/600, and its body length was 5,327 mm.Therefore, a prototype shearer and scraper conveyor whose sizes were 13.3% of the size of the original equipment were designed and manufactured. This enabled more convenient and faster experimentation (Figure10). Using the scraper conveyor prototype, we were able to achieve the following: (1) variable shapes of the scraper conveyor could be formed; (2) in a different connection state of the middle troughs, the curve formed by the pin rails directly influenced the running trajectory of the shearer; (3) in a different connection state of the middle troughs, the contacting mode between the support sliding shoe and the coal plate could be simulated; (4) a tilt sensor was installed in the middle position of every middle trough to mark the horizontal and vertical inclination angles in real time.Figure 10
(a) Shape of the scraper conveyor; (b) the two sensors installed in the shearer body; (c) the supporting sliding shoe; (d) the walking wheel; (e) the shearer prototype.Using the shearer prototype, we could achieve the following: (1) the shearer body length could be changed; (2) coupled with the coal plate, the supporting sliding shoes could self-adapt; (3) two walking wheels were perfectly replaced by two tires, which could simulate the movement of the shearer; (4) a SINS device and a double-axis tilt sensor were installed in the position of the left supporting sliding shoe.


### 3.2. Static Experiment

The shape of the scraper conveyor prototype was placed as in Figure10(a), and tilt sensors were installed on every middle trough. Each middle trough was marked with five key points, which divided the middle trough into five parts on an average.The shearer position is successively decided at every key point belonging to the five key points of each middle trough. Series values of the shearer body pitch angle were read and recorded at every key point.The datum of every middle trough tilt angle measured by the tilt sensors and SINS was imported to the Unity3d simulation software, and two simulation curves were output. The two theoretical curves of the shearer body pitch angle measured by the Unity3d simulation software and the two actual curves of the shearer body pitch angle measured by the two sensors are shown in Figure11.Figure 11
Results of the static test.As we can see from Figure11, the variation trend of the shearer body pitch angle is basically the same as that observed in the theoretical analysis; in addition, the maximum difference is 0.53°. The positioning error of the shearer was less than 0.38 times the middle trough length.


### 3.3. Dynamic Experiment

The static experiment cannot determine the properties and measurement accuracy of the sensors in the actual process of dynamic operation. Therefore, it was necessary to conduct a dynamic experiment in order to study the dynamic operation properties of the two types of sensors under the condition in which the shearer prototype could operate along with the shape of the scraper conveyor prototype automatically.After pressing the operation button, the shearer started running, and the shearer body pitch angle in the running process was recorded using two types of sensors in real time.After selecting the shearer body length as 5,327 mm, the test was conducted five times. The comparison results of the measurement values obtained using the two types of sensors and the theoretical values obtained using the VR software are shown in Figures12 and 13.Figure 12
Comparison of the results obtained using SINS with the theoretical results.Figure 13
Comparison of the results obtained using the tilt sensors with the measurement results.The analysis showed that the tilt sensor was more fluctuant in the process of shearer dynamic operation and that it was easily disturbed by environmental noise. Moreover, owing to the interdesign of filtered characteristic, the SINS showed good seismic performance.The variation trend of the shearer body pitch angle was basically the same as that observed in the theoretical analysis. However, the deviations between the two sensors and the theoretical values were greater than those obtained in the static test. Positioning correction, caused by the numerically measured value, may lead to a location error. Therefore, it was necessary to predict and correct the result in real time using the adaptive information fusion algorithm. The curves obtained after processing are shown in Figure14.Figure 14
Theoretical values and fusion values of the two sensor measurements.According to the analysis result obtained using the two sensors, the shearer’s position relative to the shape of the scraper conveyor can be reversely inferred. After processing with the adaptive fusion algorithm, the position of the shearer could meet the high level of positioning accuracy under the static condition, which was 0.38 times the middle trough length that could be reached (Table5).Table 5
Comparison of experimental results of shearer positioning (units: a middle trough length).
Theoretical value
Shearer body pitch angle measured by the tilt sensors
Shearer body pitch angle measured by SINS
Shearer body pitch angle measured according to the fusion value
Theoretical value measured by the tilt sensors
0.73
0.59
0.42
Theoretical value measured by SINS
0.67
0.49
0.45
Theoretical value measured according to the fusion value
0.53
0.47
0.38


### 3.4. Experiments under Different Body Lengths

At different shearer body lengths, the variation trends of the shearer body pitch angle were studied. The shearer body lengths were set as 4,500, 4,900, 5,327, 5,800, and 6,300 mm, which refer to a series of specialized shearer. Under these five conditions, all the experimental results were consistent with the theoretical curves (Figure15) and two conclusions were drawn.Figure 15
Variation trend of the shearer body pitch angle under different shearer body lengths.(1) A shorter shearer body length corresponded to a more backward shearer to the shape of the scraper conveyor and was more sensitive to terrain changes; a longer shearer body length corresponded to an earlier adaptation of the shearer to terrain changes and the shearer being more insensitive to terrain changes.(2) The longer the shearer body length, the smaller the relative change in the pitching angle under the same shape of the scraper conveyor.


## 3.1. Test Prototype

Three machines in our laboratory were selected as the research objects. The type of the scraper conveyor was SGZ764/630. The type of the shearer was MGTY250/600, and its body length was 5,327 mm.Therefore, a prototype shearer and scraper conveyor whose sizes were 13.3% of the size of the original equipment were designed and manufactured. This enabled more convenient and faster experimentation (Figure10). Using the scraper conveyor prototype, we were able to achieve the following: (1) variable shapes of the scraper conveyor could be formed; (2) in a different connection state of the middle troughs, the curve formed by the pin rails directly influenced the running trajectory of the shearer; (3) in a different connection state of the middle troughs, the contacting mode between the support sliding shoe and the coal plate could be simulated; (4) a tilt sensor was installed in the middle position of every middle trough to mark the horizontal and vertical inclination angles in real time.Figure 10
(a) Shape of the scraper conveyor; (b) the two sensors installed in the shearer body; (c) the supporting sliding shoe; (d) the walking wheel; (e) the shearer prototype.Using the shearer prototype, we could achieve the following: (1) the shearer body length could be changed; (2) coupled with the coal plate, the supporting sliding shoes could self-adapt; (3) two walking wheels were perfectly replaced by two tires, which could simulate the movement of the shearer; (4) a SINS device and a double-axis tilt sensor were installed in the position of the left supporting sliding shoe.


## 3.2. Static Experiment

The shape of the scraper conveyor prototype was placed as in Figure10(a), and tilt sensors were installed on every middle trough. Each middle trough was marked with five key points, which divided the middle trough into five parts on an average.The shearer position is successively decided at every key point belonging to the five key points of each middle trough. Series values of the shearer body pitch angle were read and recorded at every key point.The datum of every middle trough tilt angle measured by the tilt sensors and SINS was imported to the Unity3d simulation software, and two simulation curves were output. The two theoretical curves of the shearer body pitch angle measured by the Unity3d simulation software and the two actual curves of the shearer body pitch angle measured by the two sensors are shown in Figure11.Figure 11
Results of the static test.As we can see from Figure11, the variation trend of the shearer body pitch angle is basically the same as that observed in the theoretical analysis; in addition, the maximum difference is 0.53°. The positioning error of the shearer was less than 0.38 times the middle trough length.


## 3.3. Dynamic Experiment

The static experiment cannot determine the properties and measurement accuracy of the sensors in the actual process of dynamic operation. Therefore, it was necessary to conduct a dynamic experiment in order to study the dynamic operation properties of the two types of sensors under the condition in which the shearer prototype could operate along with the shape of the scraper conveyor prototype automatically.After pressing the operation button, the shearer started running, and the shearer body pitch angle in the running process was recorded using two types of sensors in real time.After selecting the shearer body length as 5,327 mm, the test was conducted five times. The comparison results of the measurement values obtained using the two types of sensors and the theoretical values obtained using the VR software are shown in Figures12 and 13.Figure 12
Comparison of the results obtained using SINS with the theoretical results.Figure 13
Comparison of the results obtained using the tilt sensors with the measurement results.The analysis showed that the tilt sensor was more fluctuant in the process of shearer dynamic operation and that it was easily disturbed by environmental noise. Moreover, owing to the interdesign of filtered characteristic, the SINS showed good seismic performance.The variation trend of the shearer body pitch angle was basically the same as that observed in the theoretical analysis. However, the deviations between the two sensors and the theoretical values were greater than those obtained in the static test. Positioning correction, caused by the numerically measured value, may lead to a location error. Therefore, it was necessary to predict and correct the result in real time using the adaptive information fusion algorithm. The curves obtained after processing are shown in Figure14.Figure 14
Theoretical values and fusion values of the two sensor measurements.According to the analysis result obtained using the two sensors, the shearer’s position relative to the shape of the scraper conveyor can be reversely inferred. After processing with the adaptive fusion algorithm, the position of the shearer could meet the high level of positioning accuracy under the static condition, which was 0.38 times the middle trough length that could be reached (Table5).Table 5
Comparison of experimental results of shearer positioning (units: a middle trough length).
Theoretical value
Shearer body pitch angle measured by the tilt sensors
Shearer body pitch angle measured by SINS
Shearer body pitch angle measured according to the fusion value
Theoretical value measured by the tilt sensors
0.73
0.59
0.42
Theoretical value measured by SINS
0.67
0.49
0.45
Theoretical value measured according to the fusion value
0.53
0.47
0.38


## 3.4. Experiments under Different Body Lengths

At different shearer body lengths, the variation trends of the shearer body pitch angle were studied. The shearer body lengths were set as 4,500, 4,900, 5,327, 5,800, and 6,300 mm, which refer to a series of specialized shearer. Under these five conditions, all the experimental results were consistent with the theoretical curves (Figure15) and two conclusions were drawn.Figure 15
Variation trend of the shearer body pitch angle under different shearer body lengths.(1) A shorter shearer body length corresponded to a more backward shearer to the shape of the scraper conveyor and was more sensitive to terrain changes; a longer shearer body length corresponded to an earlier adaptation of the shearer to terrain changes and the shearer being more insensitive to terrain changes.(2) The longer the shearer body length, the smaller the relative change in the pitching angle under the same shape of the scraper conveyor.


## 4. Conclusions

In this study, a joint positioning and attitude solving method was proposed and investigated for shearer and scraper conveyor under complex conditions. The following conclusions were drawn.(1) This method can provide more precise dynamic monitoring for the operation of shearer and scraper conveyor. By obtaining the shape of the scraper conveyor in real time, the shearer can be used in advance to predict and regulate the operating attitude in the current cycle.(2) Based on this method, the cutting trajectory of the front and rear drums can be calculated in real time, providing strong support for the memory cutting method in a fully mechanized coal-mining automation face.(3) No error accumulated. Hence, this positioning method of the shearer could integrate existing positioning methods, including the SINS positioning method, infrared positioning method, walking shaft encoder positioning method, and UWB wireless sensor positioning method. These methods realize a more precise positioning for the shearer in actual composite conditions.(4) This method provides support for the 3D coordinated positioning of fully mechanized coal-mining equipment.



---

*Source: 3793412-2017-10-26.xml*