TY - JOUR AU - Gao, Zhipeng AB - Introduction A scraper conveyor constitutes the core equipment in a fully mechanized mining system, which undertakes important tasks, such as shearer haulage direction, hydraulic support movement and support and coal transportation [1]. The body of the scraper conveyor is composed of several sections of the middle trough, which is welded along the trough side and middle plate. The middle trough is one of the key components of a scraper conveyor. During operation, raw coal cut by the shearer impacts the middle trough from different heights and angles. In this process, due to the impact of falling raw coal as well as the sliding friction between raw coal, scraper and scraper chains, the surface material of the middle plate can become detached, forming impact damage and wear. Ultimately, the middle trough may be damaged. According to the literature, no less than 80% of all operation faults of scraper conveyors are caused by middle trough damage [2–4]. Therefore, prolonging the service life of the middle trough and improving the fully mechanized mining efficiency have attracted extensive attention from scholars in related fields. To improve the quality of scraper conveyors and prolong their service life, scholars have invested numerous efforts. For instance, in terms of impact, scholars have analyzed the characteristics of impact damage to conveyors under impact loads and provided a method for impact energy calculation [5, 6]. In terms of fatigue, scholars have investigated fatigue crack propagation, based on which a novel crack propagation theory has been established [7, 8]. In terms of simulation model construction, researchers have proposed a discrete element impact calculation program based on the sphere model, which provides a new impact research method [9]. In terms of discrete elements, relevant software has been employed to analyze the effects of different scraper chain speeds, different falling heights of coal particles, varying impact loads and varying scraper laying lean angles on the wear patterns of the middle trough as well as conveyor efficiency [10, 11]. Regarding design, scholars have conducted research on the effects of the scraper spacing, impact velocity of bulk materials and scraper chain speed on wear of the middle trough [12]. Considering composite factors, scholars have analyzed the characteristics of impact damage to the middle trough under the combined action of coal properties, coal impact height and conveying speed based on orthogonal tests, and the weight of each factor has been obtained [13]. The influence of the load, silicon carbide abrasive particle size, sliding speed and sliding distance on the roughness of the wear surface of cast 7075 aluminum alloy has been determined [14–16]. The influence of multiple evaluation indices, such as cutting medium mechanical parameters, roller rotatory speed and impact frequency, on the working performance of point impact cutting rollers has been analyzed with the numerical analysis method [17]. In addition, scholars have studied the wear characteristics of the middle trough of a scraper conveyor. Through Plackett–Burman experiments, it has been found that the interaction between the waste rock content, moisture content and normal load can aggravate wear of the middle trough, and a regression prediction model of the wear amount has been established [18, 19]. The operation conditions of scraper conveyors are complex, and the middle trough tends to suffer early damage under the impact of raw coal and wear of scraper chains due to raw coal. Wear of the scraper chains of the middle trough has been widely studied using finite element analysis. By changing the number of sprocket teeth, chain pitch and operation speed, this problem has been satisfactorily overcome [20, 21]. However, the middle trough is subjected to the combined action of multiple factors; i.e., this part is also affected by coal fall, scraper conveyor layout and chain spacing, in addition to friction and wear. Different combinations of these factors result in different scraper conveyor working conditions. To date, most existing studies have been single-factor studies, and there have been few studies considering the combined action of multiple factors. Although [13] is a multifactor study, it is based on the properties of raw coal. In a fully mechanized mining system, the properties of raw coal are uncontrollable factors. Therefore, the above study provides no notable significance for prolonging the service life of the middle trough. Therefore, to prolong the service life of the middle trough, the key point is to control the working conditions of the middle trough. Based on the aforementioned issues, this study aims to simulate and analyze the four influencing factors of the middle trough, i.e., roll angle, front lean angle, raw coal particle size and chain spacing, to obtain the impact damage depth and normal cumulative contact energy of the middle trough under these different factors. Based on the simulation results, orthogonal matrix analysis (OMA) [22] is employed to obtain the optimal combination working conditions. The accuracy of the results is then verified through experiments. The results of this study may provide a theoretical basis for the design and layout of the middle trough of scraper conveyors. Methods In actual practice, impact damage to the middle trough of a scraper conveyor is caused by the combined action of various factors. If a comprehensive test of multiple factors is carried out, the work burden will be tremendous, which will make the task very difficult to implement. However, the very large amount of work can be greatly reduced by using the orthogonal test method. Basic theories OMA. An orthogonal table of the above four factors at three levels is constructed (Table 1). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Orthogonal table of the four factors considering three levels. https://doi.org/10.1371/journal.pone.0266831.t001 Simulation analysis based on four factors at three levels requires 34 = 81 tests. Adoption of the orthogonal method allows some test data to replace the required, comprehensive tests. Therefore, simulation results can be efficiently obtained. The orthogonal test table is summarized in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Orthogonal test table. https://doi.org/10.1371/journal.pone.0266831.t002 In an orthogonal test with m factors at n levels, the average value of the evaluation index of the i factor at the j level is . Under the condition of , an index matrix is constructed as follows [22]: We let , and a factor-layer matrix is obtained as follows [22]: Supposing that the range of factor Ai is xi and that its weight among all factors is , a horizontal layer matrix is obtained as follows [22]: Then, a weight matrix of the evaluation indices is constructed as follows [22]: where denotes the weight of factor A at the first level among all evaluation indices at all levels. Data of the four factors at the three levels. A scraper conveyor completes sliding under the pushing force of hydraulic support. Due to unevenness of the ground, horizontal, leftward or rightward inclination of the middle trough in the transverse direction occurs, and actual working conditions, such as horizontal, uphill or downhill transportation in the vertical direction, occur. Additionally, the chains exhibit layouts of middle double chains, quasi-edge double chains, edge double chains, etc. In actual practice, impact damage to the middle slot of the scraper conveyor can also be caused by the waste rock content, apart from the transverse laying roll angle, front lean angle, raw coal particle size, and chain layout and spacing. However, due to the uncontrollable nature of the waste rock content, this factor is not considered in this study. Therefore, in this study, four factors are considered at three levels (Table 3), which include 1) the transverse laying roll angle: 0°, 5° and -5°; 2) front lean angle: 0°, 10° and -10°; 3) chain layout and spacing: middle double-chain layout with a spacing of 140 mm, quasi-edge double-chain layout with a spacing of 420 mm and edge double-chain layout with a spacing of 700 mm; and 4) coal particle size: 35 mm, 70 mm and 105 mm (0.5-, 1.0- and 1.5-fold, respectively, the basic coal particle size (70 mm)). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Four factors at the three levels for impact damage. https://doi.org/10.1371/journal.pone.0266831.t003 Simulation experiment The middle trough is impacted by falling coal, and its surface material can fall off to form impact pits. In this study, the impact damage due to raw coal to the middle plate of the middle trough is characterized by the volume of impact pits. Impact damage theory [21]. The volume of impact pits △V is calculated as follows: (1) where d denotes the contact diameter of the raw coal with the middle plate (mm), L denotes the scratch length of the impact damage (mm), h denotes the depth of the impact pit (mm) and R denotes the radius of the coal particle (mm). In the process of raw coal impacting the middle plate, the impact load on the middle plate is decomposed into the normal impact force Fn and the tangential impact force Fτ. The relationship between these forces is as follows: (2) where μ denotes the friction coefficient between the raw coal and middle plate. According to the conservation of energy, the following can be obtained: (3) where m denotes the mass of the raw coal particles (kg), vτ denotes the tangential velocity of the raw coal (m/s), and κτ denotes the tangential restitution coefficient of the impact force of the raw coal on the middle plate, among which the following applies: where a is the angle between the incident direction and the horizontal plane when raw coal impacts the middle plate (°), kn is the normal restitution coefficient of the impact of raw coal on the middle plate, and vy is the normal minimum velocity when the middle plate yields (m/s). We suppose that the yielding process of the contact surface of the middle plate satisfies the Mohr–Coulomb criteria: (4) where k is the cohesion force of the middle plate (N) and E is the equivalent elastic modulus (MPa), among which: where E1 is the elastic modulus of coal (MPa), E2 is the elastic modulus of the middle plate (MPa), δ1 is the Poisson ratio of coal, δ2 is the Poisson ratio of the middle plate and β is the internal friction angle of coal. K1, K2 and ξ are calculation parameters, as follows: Based on Eqs (1) and (4), the following can be obtained: (5) Based on Eq (5), the pit volume formed by the impact damage due to raw coal to the middle plate can be characterized by the function f. The first part of function f captures the basic characteristics of raw coal, and the second part captures the contact parameter between the raw coal and middle plate. The raw coal impacts the middle trough plate, and the impact can make the surface material of the middle trough plate fall off. Suppose that the impact wear volume dV of the raw coal acting on the middle trough plate is directly proportional to the impact energy dA. Then, the following applies: (6) where η is the dimensionless impact wear rate, which represents the wear amount per unit impact, and ψ is the impact energy conversion rate. Finite element modeling. The simulation research object in this study is an 830-series scraper conveyor. Its conveying capacity is 1000 t/h, with a conveyor length of 200 m. This scraper conveyor adopts double-chain scraper chains with a scraper spacing of 108 mm, a running speed of 1.0 m/s, and a middle trough width of 780 mm. The length of a single segment of the middle trough is 1.5 m, with five segments in total; that is, the total length of the middle trough is 7.5 m. The installed power of the conveyor is 2×(250–315) kW. Fig 1 shows the assembly scheme of the middle trough and scraper chains. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Assembly scheme of the middle trough of the scraper conveyor: 1-scraper, 2-chain, and 3-middle trough. https://doi.org/10.1371/journal.pone.0266831.g001 The raw coal is discretized to study the whole dynamic process of the impact damage due to raw coal to the middle trough of the scraper conveyor. Here, we suppose that when raw coal particles impact the middle plate of the middle trough, the contact between the raw coal particles and middle plate is subject to a point contact. Based on this supposition, spheres are used to replace raw coal particles to produce impact with the middle plate [23]. A particle impact model (Fig 2) is established in SolidWorks software (Dassault Systèmes, France), and the model is stored as an X_T file. The file is then introduced into the Explicit Dynamics module of ANSYS. In the model, single segments of the middle trough are treated as impacted objects with small balls as impacting objects. To ensure the accuracy of the simulation results and reduce the amount of computation, Hypermesh software is used to mesh the segment of the middle trough. Free meshing for tetrahedral solid elements is adopted. The unit length is set at 10 mm, with the minimal unit length at 2 mm, and the characteristic angle is set at 30°. A total of 1,283,942 units and 337,150 nodes are produced (Fig 3). Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Solid model of the middle trough (a) and raw coal sphere model (b). https://doi.org/10.1371/journal.pone.0266831.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Grid model of a single section of the middle trough. https://doi.org/10.1371/journal.pone.0266831.g003 To complete simulation in EDEM software, the material property parameters of the middle trough and raw coal impact spheres (Table 4) as well as the properties of the contact between these objects must be set [24, 25] (Table 5). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Material properties of the middle trough and raw coal impact spheres. https://doi.org/10.1371/journal.pone.0266831.t004 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Properties of the material contact. https://doi.org/10.1371/journal.pone.0266831.t005 Regarding the particle plant created during EDEM simulation, the particle generation rate should meet the production rate of the 830-series scraper conveyor studied in this study, i.e., 1000 t/h. The speed of the generated raw coal particles along the moving direction of the scraper is 0.15 m/s, the rate of raw coal particle generation is 277 kg/s with a total weight of 1500 kg, and the simulation time step is set to 20% (50.15602 s) with a total simulation time of 5 s and a storage time interval of 0.1 s [26]. The distance between the impacting raw coal spheres and the middle trough is 0.001 mm, with a speed of the impacting spheres of 4.47 m/s [27, 28]. The simulation boundary conditions are set as follows: with the fixed simulation model and constant falling height and speed of materials, the falling height of the raw coal particles is 1 m, the running speed of the scraper chain is 1.0 m/s, the scraper spacing is 108 mm, and the flow of raw coal is 277 kg/s. Validation experiment To verify the accuracy of the results obtained via OMA, we conducted an experiment. A scraper conveyor at a coal mine in Shanxi was tested by the Taiyuan Research Institute Co., Ltd., of the China Coal Science and Industry Group. The test equipment is shown in Fig 4. A power supply coil and lithium battery were set inside the scraper, and induction coils were set inside the middle plate of the trough. While the scraper conveyor was running, the scraper moved along the middle trough. When the scraper with an embedded power supply coil passed through the middle plate embedded with an induction coil, a voltage was generated at both ends of the induction coil. The voltage signal was collected by an acquisition module and transmitted to an industrial control computer. With increasing wear of the middle plate of the trough, the wear amount of the middle plate was calculated by an industrial control computer [29]. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Scraper conveyor in the validation test: 1-scraper, 2-middle trough. https://doi.org/10.1371/journal.pone.0266831.g004 Basic theories OMA. An orthogonal table of the above four factors at three levels is constructed (Table 1). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Orthogonal table of the four factors considering three levels. https://doi.org/10.1371/journal.pone.0266831.t001 Simulation analysis based on four factors at three levels requires 34 = 81 tests. Adoption of the orthogonal method allows some test data to replace the required, comprehensive tests. Therefore, simulation results can be efficiently obtained. The orthogonal test table is summarized in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Orthogonal test table. https://doi.org/10.1371/journal.pone.0266831.t002 In an orthogonal test with m factors at n levels, the average value of the evaluation index of the i factor at the j level is . Under the condition of , an index matrix is constructed as follows [22]: We let , and a factor-layer matrix is obtained as follows [22]: Supposing that the range of factor Ai is xi and that its weight among all factors is , a horizontal layer matrix is obtained as follows [22]: Then, a weight matrix of the evaluation indices is constructed as follows [22]: where denotes the weight of factor A at the first level among all evaluation indices at all levels. Data of the four factors at the three levels. A scraper conveyor completes sliding under the pushing force of hydraulic support. Due to unevenness of the ground, horizontal, leftward or rightward inclination of the middle trough in the transverse direction occurs, and actual working conditions, such as horizontal, uphill or downhill transportation in the vertical direction, occur. Additionally, the chains exhibit layouts of middle double chains, quasi-edge double chains, edge double chains, etc. In actual practice, impact damage to the middle slot of the scraper conveyor can also be caused by the waste rock content, apart from the transverse laying roll angle, front lean angle, raw coal particle size, and chain layout and spacing. However, due to the uncontrollable nature of the waste rock content, this factor is not considered in this study. Therefore, in this study, four factors are considered at three levels (Table 3), which include 1) the transverse laying roll angle: 0°, 5° and -5°; 2) front lean angle: 0°, 10° and -10°; 3) chain layout and spacing: middle double-chain layout with a spacing of 140 mm, quasi-edge double-chain layout with a spacing of 420 mm and edge double-chain layout with a spacing of 700 mm; and 4) coal particle size: 35 mm, 70 mm and 105 mm (0.5-, 1.0- and 1.5-fold, respectively, the basic coal particle size (70 mm)). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Four factors at the three levels for impact damage. https://doi.org/10.1371/journal.pone.0266831.t003 OMA. An orthogonal table of the above four factors at three levels is constructed (Table 1). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Orthogonal table of the four factors considering three levels. https://doi.org/10.1371/journal.pone.0266831.t001 Simulation analysis based on four factors at three levels requires 34 = 81 tests. Adoption of the orthogonal method allows some test data to replace the required, comprehensive tests. Therefore, simulation results can be efficiently obtained. The orthogonal test table is summarized in Table 2. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Orthogonal test table. https://doi.org/10.1371/journal.pone.0266831.t002 In an orthogonal test with m factors at n levels, the average value of the evaluation index of the i factor at the j level is . Under the condition of , an index matrix is constructed as follows [22]: We let , and a factor-layer matrix is obtained as follows [22]: Supposing that the range of factor Ai is xi and that its weight among all factors is , a horizontal layer matrix is obtained as follows [22]: Then, a weight matrix of the evaluation indices is constructed as follows [22]: where denotes the weight of factor A at the first level among all evaluation indices at all levels. Data of the four factors at the three levels. A scraper conveyor completes sliding under the pushing force of hydraulic support. Due to unevenness of the ground, horizontal, leftward or rightward inclination of the middle trough in the transverse direction occurs, and actual working conditions, such as horizontal, uphill or downhill transportation in the vertical direction, occur. Additionally, the chains exhibit layouts of middle double chains, quasi-edge double chains, edge double chains, etc. In actual practice, impact damage to the middle slot of the scraper conveyor can also be caused by the waste rock content, apart from the transverse laying roll angle, front lean angle, raw coal particle size, and chain layout and spacing. However, due to the uncontrollable nature of the waste rock content, this factor is not considered in this study. Therefore, in this study, four factors are considered at three levels (Table 3), which include 1) the transverse laying roll angle: 0°, 5° and -5°; 2) front lean angle: 0°, 10° and -10°; 3) chain layout and spacing: middle double-chain layout with a spacing of 140 mm, quasi-edge double-chain layout with a spacing of 420 mm and edge double-chain layout with a spacing of 700 mm; and 4) coal particle size: 35 mm, 70 mm and 105 mm (0.5-, 1.0- and 1.5-fold, respectively, the basic coal particle size (70 mm)). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Four factors at the three levels for impact damage. https://doi.org/10.1371/journal.pone.0266831.t003 Simulation experiment The middle trough is impacted by falling coal, and its surface material can fall off to form impact pits. In this study, the impact damage due to raw coal to the middle plate of the middle trough is characterized by the volume of impact pits. Impact damage theory [21]. The volume of impact pits △V is calculated as follows: (1) where d denotes the contact diameter of the raw coal with the middle plate (mm), L denotes the scratch length of the impact damage (mm), h denotes the depth of the impact pit (mm) and R denotes the radius of the coal particle (mm). In the process of raw coal impacting the middle plate, the impact load on the middle plate is decomposed into the normal impact force Fn and the tangential impact force Fτ. The relationship between these forces is as follows: (2) where μ denotes the friction coefficient between the raw coal and middle plate. According to the conservation of energy, the following can be obtained: (3) where m denotes the mass of the raw coal particles (kg), vτ denotes the tangential velocity of the raw coal (m/s), and κτ denotes the tangential restitution coefficient of the impact force of the raw coal on the middle plate, among which the following applies: where a is the angle between the incident direction and the horizontal plane when raw coal impacts the middle plate (°), kn is the normal restitution coefficient of the impact of raw coal on the middle plate, and vy is the normal minimum velocity when the middle plate yields (m/s). We suppose that the yielding process of the contact surface of the middle plate satisfies the Mohr–Coulomb criteria: (4) where k is the cohesion force of the middle plate (N) and E is the equivalent elastic modulus (MPa), among which: where E1 is the elastic modulus of coal (MPa), E2 is the elastic modulus of the middle plate (MPa), δ1 is the Poisson ratio of coal, δ2 is the Poisson ratio of the middle plate and β is the internal friction angle of coal. K1, K2 and ξ are calculation parameters, as follows: Based on Eqs (1) and (4), the following can be obtained: (5) Based on Eq (5), the pit volume formed by the impact damage due to raw coal to the middle plate can be characterized by the function f. The first part of function f captures the basic characteristics of raw coal, and the second part captures the contact parameter between the raw coal and middle plate. The raw coal impacts the middle trough plate, and the impact can make the surface material of the middle trough plate fall off. Suppose that the impact wear volume dV of the raw coal acting on the middle trough plate is directly proportional to the impact energy dA. Then, the following applies: (6) where η is the dimensionless impact wear rate, which represents the wear amount per unit impact, and ψ is the impact energy conversion rate. Finite element modeling. The simulation research object in this study is an 830-series scraper conveyor. Its conveying capacity is 1000 t/h, with a conveyor length of 200 m. This scraper conveyor adopts double-chain scraper chains with a scraper spacing of 108 mm, a running speed of 1.0 m/s, and a middle trough width of 780 mm. The length of a single segment of the middle trough is 1.5 m, with five segments in total; that is, the total length of the middle trough is 7.5 m. The installed power of the conveyor is 2×(250–315) kW. Fig 1 shows the assembly scheme of the middle trough and scraper chains. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Assembly scheme of the middle trough of the scraper conveyor: 1-scraper, 2-chain, and 3-middle trough. https://doi.org/10.1371/journal.pone.0266831.g001 The raw coal is discretized to study the whole dynamic process of the impact damage due to raw coal to the middle trough of the scraper conveyor. Here, we suppose that when raw coal particles impact the middle plate of the middle trough, the contact between the raw coal particles and middle plate is subject to a point contact. Based on this supposition, spheres are used to replace raw coal particles to produce impact with the middle plate [23]. A particle impact model (Fig 2) is established in SolidWorks software (Dassault Systèmes, France), and the model is stored as an X_T file. The file is then introduced into the Explicit Dynamics module of ANSYS. In the model, single segments of the middle trough are treated as impacted objects with small balls as impacting objects. To ensure the accuracy of the simulation results and reduce the amount of computation, Hypermesh software is used to mesh the segment of the middle trough. Free meshing for tetrahedral solid elements is adopted. The unit length is set at 10 mm, with the minimal unit length at 2 mm, and the characteristic angle is set at 30°. A total of 1,283,942 units and 337,150 nodes are produced (Fig 3). Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Solid model of the middle trough (a) and raw coal sphere model (b). https://doi.org/10.1371/journal.pone.0266831.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Grid model of a single section of the middle trough. https://doi.org/10.1371/journal.pone.0266831.g003 To complete simulation in EDEM software, the material property parameters of the middle trough and raw coal impact spheres (Table 4) as well as the properties of the contact between these objects must be set [24, 25] (Table 5). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Material properties of the middle trough and raw coal impact spheres. https://doi.org/10.1371/journal.pone.0266831.t004 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Properties of the material contact. https://doi.org/10.1371/journal.pone.0266831.t005 Regarding the particle plant created during EDEM simulation, the particle generation rate should meet the production rate of the 830-series scraper conveyor studied in this study, i.e., 1000 t/h. The speed of the generated raw coal particles along the moving direction of the scraper is 0.15 m/s, the rate of raw coal particle generation is 277 kg/s with a total weight of 1500 kg, and the simulation time step is set to 20% (50.15602 s) with a total simulation time of 5 s and a storage time interval of 0.1 s [26]. The distance between the impacting raw coal spheres and the middle trough is 0.001 mm, with a speed of the impacting spheres of 4.47 m/s [27, 28]. The simulation boundary conditions are set as follows: with the fixed simulation model and constant falling height and speed of materials, the falling height of the raw coal particles is 1 m, the running speed of the scraper chain is 1.0 m/s, the scraper spacing is 108 mm, and the flow of raw coal is 277 kg/s. Impact damage theory [21]. The volume of impact pits △V is calculated as follows: (1) where d denotes the contact diameter of the raw coal with the middle plate (mm), L denotes the scratch length of the impact damage (mm), h denotes the depth of the impact pit (mm) and R denotes the radius of the coal particle (mm). In the process of raw coal impacting the middle plate, the impact load on the middle plate is decomposed into the normal impact force Fn and the tangential impact force Fτ. The relationship between these forces is as follows: (2) where μ denotes the friction coefficient between the raw coal and middle plate. According to the conservation of energy, the following can be obtained: (3) where m denotes the mass of the raw coal particles (kg), vτ denotes the tangential velocity of the raw coal (m/s), and κτ denotes the tangential restitution coefficient of the impact force of the raw coal on the middle plate, among which the following applies: where a is the angle between the incident direction and the horizontal plane when raw coal impacts the middle plate (°), kn is the normal restitution coefficient of the impact of raw coal on the middle plate, and vy is the normal minimum velocity when the middle plate yields (m/s). We suppose that the yielding process of the contact surface of the middle plate satisfies the Mohr–Coulomb criteria: (4) where k is the cohesion force of the middle plate (N) and E is the equivalent elastic modulus (MPa), among which: where E1 is the elastic modulus of coal (MPa), E2 is the elastic modulus of the middle plate (MPa), δ1 is the Poisson ratio of coal, δ2 is the Poisson ratio of the middle plate and β is the internal friction angle of coal. K1, K2 and ξ are calculation parameters, as follows: Based on Eqs (1) and (4), the following can be obtained: (5) Based on Eq (5), the pit volume formed by the impact damage due to raw coal to the middle plate can be characterized by the function f. The first part of function f captures the basic characteristics of raw coal, and the second part captures the contact parameter between the raw coal and middle plate. The raw coal impacts the middle trough plate, and the impact can make the surface material of the middle trough plate fall off. Suppose that the impact wear volume dV of the raw coal acting on the middle trough plate is directly proportional to the impact energy dA. Then, the following applies: (6) where η is the dimensionless impact wear rate, which represents the wear amount per unit impact, and ψ is the impact energy conversion rate. Finite element modeling. The simulation research object in this study is an 830-series scraper conveyor. Its conveying capacity is 1000 t/h, with a conveyor length of 200 m. This scraper conveyor adopts double-chain scraper chains with a scraper spacing of 108 mm, a running speed of 1.0 m/s, and a middle trough width of 780 mm. The length of a single segment of the middle trough is 1.5 m, with five segments in total; that is, the total length of the middle trough is 7.5 m. The installed power of the conveyor is 2×(250–315) kW. Fig 1 shows the assembly scheme of the middle trough and scraper chains. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Assembly scheme of the middle trough of the scraper conveyor: 1-scraper, 2-chain, and 3-middle trough. https://doi.org/10.1371/journal.pone.0266831.g001 The raw coal is discretized to study the whole dynamic process of the impact damage due to raw coal to the middle trough of the scraper conveyor. Here, we suppose that when raw coal particles impact the middle plate of the middle trough, the contact between the raw coal particles and middle plate is subject to a point contact. Based on this supposition, spheres are used to replace raw coal particles to produce impact with the middle plate [23]. A particle impact model (Fig 2) is established in SolidWorks software (Dassault Systèmes, France), and the model is stored as an X_T file. The file is then introduced into the Explicit Dynamics module of ANSYS. In the model, single segments of the middle trough are treated as impacted objects with small balls as impacting objects. To ensure the accuracy of the simulation results and reduce the amount of computation, Hypermesh software is used to mesh the segment of the middle trough. Free meshing for tetrahedral solid elements is adopted. The unit length is set at 10 mm, with the minimal unit length at 2 mm, and the characteristic angle is set at 30°. A total of 1,283,942 units and 337,150 nodes are produced (Fig 3). Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. Solid model of the middle trough (a) and raw coal sphere model (b). https://doi.org/10.1371/journal.pone.0266831.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Grid model of a single section of the middle trough. https://doi.org/10.1371/journal.pone.0266831.g003 To complete simulation in EDEM software, the material property parameters of the middle trough and raw coal impact spheres (Table 4) as well as the properties of the contact between these objects must be set [24, 25] (Table 5). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Material properties of the middle trough and raw coal impact spheres. https://doi.org/10.1371/journal.pone.0266831.t004 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Properties of the material contact. https://doi.org/10.1371/journal.pone.0266831.t005 Regarding the particle plant created during EDEM simulation, the particle generation rate should meet the production rate of the 830-series scraper conveyor studied in this study, i.e., 1000 t/h. The speed of the generated raw coal particles along the moving direction of the scraper is 0.15 m/s, the rate of raw coal particle generation is 277 kg/s with a total weight of 1500 kg, and the simulation time step is set to 20% (50.15602 s) with a total simulation time of 5 s and a storage time interval of 0.1 s [26]. The distance between the impacting raw coal spheres and the middle trough is 0.001 mm, with a speed of the impacting spheres of 4.47 m/s [27, 28]. The simulation boundary conditions are set as follows: with the fixed simulation model and constant falling height and speed of materials, the falling height of the raw coal particles is 1 m, the running speed of the scraper chain is 1.0 m/s, the scraper spacing is 108 mm, and the flow of raw coal is 277 kg/s. Validation experiment To verify the accuracy of the results obtained via OMA, we conducted an experiment. A scraper conveyor at a coal mine in Shanxi was tested by the Taiyuan Research Institute Co., Ltd., of the China Coal Science and Industry Group. The test equipment is shown in Fig 4. A power supply coil and lithium battery were set inside the scraper, and induction coils were set inside the middle plate of the trough. While the scraper conveyor was running, the scraper moved along the middle trough. When the scraper with an embedded power supply coil passed through the middle plate embedded with an induction coil, a voltage was generated at both ends of the induction coil. The voltage signal was collected by an acquisition module and transmitted to an industrial control computer. With increasing wear of the middle plate of the trough, the wear amount of the middle plate was calculated by an industrial control computer [29]. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Scraper conveyor in the validation test: 1-scraper, 2-middle trough. https://doi.org/10.1371/journal.pone.0266831.g004 Results and discussion Within the context of explicit dynamics, research on impact can be performed only with a model under transient conditions. Therefore, to simplify the falling process of raw coal particles, the spacing between the impacting spheres and the middle trough is set to 0.001 mm, and the velocity is set to 4.470 m/s in this study. Simulation results Based on ANSYS simulation software, curves of the middle trough energy and the kinetic energy of raw coal impacting particles during the impact process are obtained (Fig 5). Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 5. Curve of the impact energy. CKE, coal kinetic energy. MTE, energy of the middle trough. https://doi.org/10.1371/journal.pone.0266831.g005 In the impact area, the velocity of the middle plate reaches a maximum value of 112.39 mm/s. From the beginning of middle plate impact to rebound, the kinetic energy of the coal particle spheres first decreases and then increases while the energy of the middle plate and trough sides gradually increases. During the impact process, the energy of the middle trough increases by h2, and the kinetic energy of the raw coal impact particles decreases by h1. At this moment, the conversion rate of the impact energy according to impact theory is calculated as follows: (7) Based on the simulation results, h1 = 2.699–2.141 = 0.558 and h2 = 0.221, ψ = 39.6% are obtained. A correction is made according to the above finite element simulation equation (Eq (6)) as follows: Based on the ANSYS simulation, from falling to landing on the middle trough, approximately 39% of the dynamic energy of raw coal is absorbed by the middle trough and transformed into impact energy. Grincova et al. [5] monitored the energy absorbed by the conveyor belt at the material receiving place of the general belt conveyor and found that the impact damage to the conveyor belt is associated with the dynamic energy loss of the falling materials. In EDEM software, the impact of raw coal on the middle trough under the conditions of the four factors at the three levels is simulated separately. An impact damage depth curve of raw coal on the middle groove under different transverse roll angles is obtained (Fig 6), and the normal cumulative contact energy curve is shown in Fig 7. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 6. Impact damage depth curve of the middle trough under the different transverse roll angles. https://doi.org/10.1371/journal.pone.0266831.g006 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 7. Normal cumulative contact energy curve of the middle trough under the different transverse roll angles. https://doi.org/10.1371/journal.pone.0266831.g007 As shown in Figs 6 and 7, the process of raw coal falling occurs before the time point of 0.5 s. Between 0.5 s and 5 s, the average impact damage depth and normal cumulative contact energy of the middle trough reach maximum values when the transverse roll angle is 0°, followed by -5° and then 5° (Table 6). The simulation results indicate that when there is a transverse roll angle of the middle trough, the load components of the normal impact of raw coal falling onto the middle trough change, and the middle trough is mainly affected by the normal impact load of the raw coal. Therefore, the occurrence of a transverse roll angle of the middle trough results in changes in the average impact damage depth and cumulative contact energy of the middle plate. In addition, as shown in Table 6, when there is a transverse roll angle of the middle trough, particularly with a positive value, the angle causes the load of raw coal to fall onto the middle trough to produce components, which reduces the contact load on the middle trough, thereby alleviating the impact damage to the middle trough. This finding is basically consistent with that reported in the literature [16], according to which an increase in contact pressure aggravates the wear of the middle trough. When the transverse roll angle is 0°, the contact load of falling coal on the middle trough is increased, as no components are produced, which leads to the most severe damage to the middle trough. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 6. Simulation results of the impact damage to the middle plate under the different transverse roll angles. https://doi.org/10.1371/journal.pone.0266831.t006 Considering the different front lean angles, the scraper conveyor is operated under different transportation conditions, such as upward, downward and horizontal transportation. Impact damage depth and normal cumulative contact energy curves are shown in Figs 8 and 9, respectively. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 8. Impact damage depth curve of the middle trough under the different front lean angles. https://doi.org/10.1371/journal.pone.0266831.g008 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 9. Normal accumulative contact energy curve of the middle trough under the different front lean angles. https://doi.org/10.1371/journal.pone.0266831.g009 When the front lean angle is positive (i.e., 10°), the scraper conveyor occurs in the upward transportation state. The falling of raw coal onto the middle trough generates a velocity component opposite to that of the scraper chain. The reduction in velocity reduces the kinetic energy of raw coal and minimizes the impact on the middle trough and the normal cumulative contact energy. When the front lean angle is negative (-10°), the scraper conveyor occurs in the downward transportation state. The falling of raw coal onto the middle trough generates a velocity component consistent with that of the scraper chain. Superposition of the velocity increases the kinetic energy of the raw coal, leading to a great impact on the middle trough. At this moment, the normal cumulative contact energy is also high. When the front lean angle is 0°, the scraper conveyor occurs in the horizontal transportation state. The falling of raw coal directly impacts the middle trough. There is no velocity component of up- or downward transportation. All kinetic energy acts on the middle trough, resulting in the maximum impact upon the middle trough. The normal cumulative contact energy also reaches a maximum value. Our findings are consistent with the result reported by Sinha and Mukhopadhyay [16]. The results are summarized in Table 7. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 7. Simulation results of the impact damage to the middle plate under the different front lean angles. https://doi.org/10.1371/journal.pone.0266831.t007 Impact damage depth and normal accumulative contact energy curves of the middle trough under the different chain spaces are shown in Figs 10 and 11, respectively. When the chain layout and spacing differ, the contact area between the raw coal and middle plate also changes. According to the analysis of Figs 10 and 11, under the edge-side double-chain layout, the average impact damage depth of the middle plate is the deepest, and the cumulative normal contact energy is the highest, followed by the quasi-edge double-chain and then middle double-chain layouts. With different chain layouts, the wear of the middle trough of the scraper conveyor also differs, and therefore, the selection of the chain layout plays a critical role in the design of a scraper conveyor [10]. Our results are summarized in Table 8. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 10. Impact damage depth curve of the middle trough under the different chain spacings. https://doi.org/10.1371/journal.pone.0266831.g010 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 11. Normal accumulative contact energy curve of the middle trough under the different chain spacings. https://doi.org/10.1371/journal.pone.0266831.g011 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 8. Simulation results of the impact damage under the different chain spacings. https://doi.org/10.1371/journal.pone.0266831.t008 Impact damage depth and normal accumulative contact energy curves of the middle trough for the different coal particle sizes are shown in Figs 12 and 13, respectively. The larger the particle size of raw coal is, the deeper the average impact damage depth of the plate and the higher the normal cumulative contact energy. Shi and Zhu [3] and Glushkova [4] investigated the influence of the particle size and moisture content of raw coal on the damage to the middle trough, and they found that a larger particle size leads to more prominent damage to the middle trough. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 12. Curve of the impact damage depth of the middle trough caused by raw coal with different particle sizes. https://doi.org/10.1371/journal.pone.0266831.g012 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 13. Curve of the normal accumulative contact energy of the middle trough caused by raw coal with different particle sizes. https://doi.org/10.1371/journal.pone.0266831.g013 Orthogonal results Based on the simulation results related to the four factors at the three levels, as summarized in Tables 1 and 6–9, the impact damage depth and normal accumulative contact energy of the middle trough under the different conditions are obtained (Tables 10 and 11, respectively; S1 File). Download: PPT PowerPoint slide PNG larger image TIFF original image Table 9. Simulation results of the impact damage caused by the different sizes of raw coal particles. https://doi.org/10.1371/journal.pone.0266831.t009 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 10. Impact damage depth (mm) under the different conditions. https://doi.org/10.1371/journal.pone.0266831.t010 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 11. Normal accumulative contact energy (kJ) under the different conditions. https://doi.org/10.1371/journal.pone.0266831.t011 In the tables, p1i denotes the sum of the impact damage depth under the four factors at the ith level, denotes the average of the impact damage depth under each factor at the ith level, p2i denotes the sum of the accumulative contact energy under the four factors at the ith level, and denotes the average of the accumulative contact energy under each factor at the ith level. Based on matrix numerical analysis [30], an impact damage depth matrix P1 and normal accumulative contact energy matrix P2 are obtained as follows: A factor layer matrix is obtained as follows: A horizontal layer matrix is obtained as follows: An evaluation index weight matrix is obtained as follows: Averages of the impact damage depth evaluation index weight matrix G1 and the normal accumulative contact energy evaluation index weight matrix G2 are calculated, and a total weight matrix is obtained as follows: Based on the weight coefficient, the weight of each factor in impact damage to the middle trough is determined: where U denotes the weight coefficient of the influencing factor and u denotes the value of the evaluation index weight matrix parameter (here, A, B, C and D). Data are introduced, and the following values are obtained: The smaller the impact damage depth and normal cumulative contact energy are, the better the effect. According to matrix numerical analysis, the influence of the four factors on the evaluation indices of the impact damage depth and normal cumulative contact energy can be ordered as C > D > A > B, namely, raw coal particle size > chain spacing > roll angle > front lean angle. The Pareto diagram of the impact damage factors of the middle trough is obtained from the weight coefficients of the influencing factors (Fig 14). Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 14. Pareto diagram of the impact damage factors of the middle trough. https://doi.org/10.1371/journal.pone.0266831.g014 Based on Tables 4 and 5 as well as the total weight matrix G, the three horizontal conditions of the four factors affecting impact damage to the middle trough can be ordered as A3 < A2 D > A > B, namely, raw coal particle size > chain spacing > roll angle > front lean angle. The Pareto diagram of the impact damage factors of the middle trough is obtained from the weight coefficients of the influencing factors (Fig 14). Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 14. Pareto diagram of the impact damage factors of the middle trough. https://doi.org/10.1371/journal.pone.0266831.g014 Based on Tables 4 and 5 as well as the total weight matrix G, the three horizontal conditions of the four factors affecting impact damage to the middle trough can be ordered as A3 < A2