A Stress Analysis of a Conical Pick by Establishing a 3D ES-FEM Model and Using Experimental Measured Forces

: This paper presents a procedure for an accurate and reliable stress analysis in a conical pick used in mining operations, aiming to improve their wear resistance. This is achieved by (1) establishment of a three-dimensional (3D) edge-based smoothed ﬁnite element method (ES-FEM), and algorithms of creating the smoothing domain for accurate solution in terms of stress and strain distributions; and (2) use of experimentally measured actual forces using a full-scale rotary cutting machine. In our 3D ES-FEM model, the physical domain for the pick is ﬁrst discredited using linear triangular elements that can be generated easily for complicated geometries. The smoothing domains are then constructed based on edges of these elements in an automated fashion. In order to create the smoothing domains for the smoothed strain computation in the ES-FEM, an algorithm is presented for establishing connection between nodes, edges, faces, and elements. Each smoothing domain is bounded by a set of enclosed line-segments, besides, leading to a connectivity list for later e ﬀ ective computation. To show the e ﬀ ectiveness and accuracy of the strain energy and the displacement solution of ES-FEM, based on the actually measured forces from the laboratory rock cutting tests with a single pick, a comparison study is carried out against the standard ﬁnite element method (FEM). It can be concluded that ES-FEM has a higher convergence in energy norm and better accuracy than FEM using the same mesh from the comparison results. The 3D ES-FEM model solves the problem of the lower solution accuracy, caused by the poor quality of mesh, by using the standard FEM in solving the stress distribution of mining machinery parts, such as picks, and o ﬀ ers accurate and reliable solutions that are critical for improving the wear resistance of the pick for the mining industry.


Introduction
An imbalance in the excavation ratio has been a major obstacle affecting the safe and efficient production of coal. Especially, with the increase of rock hardness, various problems, such as pick wear and component failure, will arise, which often lead to lower efficiency, and even cause the machinery to not work. Therefore, if the mechanical properties of the cutter head for rock crushing and the wear mechanism of picks can be accurately assessed, the cutting performance of road headers will be greatly improved for higher efficiency industrial production. However, the working condition of a pick is often severe, complex, and dynamic in nature. Conical picks work usually under high impact and stress conditions, resulting in frequent failure. The failure modes include premature wear, carbide tip drop off, fracture, and normal wear [1][2][3][4][5][6][7]. In order to understand the cutting properties of a pick, researchers have carried out substantial experimental work on pick cutting force prediction. more complex, it is unrealistic to obtain higher accuracy relying on higher older elements. Nevertheless, ES-FEM can solve the problems easily, which can guarantee the nearly quadratic accuracy only using linear elements.
Thus, in this paper, for the analysis of the coal mining equipment, the ES-FEM is presented to reveal the stress and strain distributions of the conical pick, which can better master its working condition and effectively improve its wear resistance for engineers. T4 elements are employed as the background mesh and the edge-based smoothing domains are created based on it. Meanwhile, an algorithm is given for establishing connection between nodes, edges, faces, and elements. Then the smoothed strain-displacement matrices are constructed following the S-FEM theory for the smoothed stiffness matrices. Next the discrete linear algebraic system of equations is established like the procedure in FEM. What is more, the stress and strain distributions for conical picks are calculated, and the values of nodal displacement are calculated. Finally, the ES-FEM results for the conical pick with experimental measured forces were given and some comparisons with other S-FEM models and FEM were made.

Determination of the Pick Force
The cutter head operates with a composite motion, which includes rotating and swinging, and also through the conical tip of picks to break up the rocks. Due to the complex geological conditions in tunneling, it is very difficult to measure the cutting force imposed on the cutter head. Hence, the data obtained from the laboratory rock cutting tests with a single pick are usually used as the basis for calculating the forces of each cutter head pick. The prediction formulas of the pick forces obtained from the laboratory tests are more reasonable and reliable than those obtained in the theoretical and numerical simulation methods. In the paper, the cutting force acting on the pick is shown in Figure 1. The parameters present that d is the depth of cut, δ is the angle of attack, and n is the speed. F C is the cutting force, which is opposite to the direction of motion acting on the tip of the pick, and it is to separate the cuttings along the direction of the pick from the rock mass and forms a trough. F N is the normal force whose component force is perpendicular to the direction of motion of the coal rock, and it is to let the pick compress the rock and form the radiant point of the crack extension. Lastly, F S is the lateral force, from the normal force to the cutting direction according to the right-hand rule, and it is used to promote the rotation of the pick and forms a uniform wear.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 17 Nevertheless, ES-FEM can solve the problems easily, which can guarantee the nearly quadratic accuracy only using linear elements. Thus, in this paper, for the analysis of the coal mining equipment, the ES-FEM is presented to reveal the stress and strain distributions of the conical pick, which can better master its working condition and effectively improve its wear resistance for engineers. T4 elements are employed as the background mesh and the edge-based smoothing domains are created based on it. Meanwhile, an algorithm is given for establishing connection between nodes, edges, faces, and elements. Then the smoothed strain-displacement matrices are constructed following the S-FEM theory for the smoothed stiffness matrices. Next the discrete linear algebraic system of equations is established like the procedure in FEM. What is more, the stress and strain distributions for conical picks are calculated, and the values of nodal displacement are calculated. Finally, the ES-FEM results for the conical pick with experimental measured forces were given and some comparisons with other S-FEM models and FEM were made.

Determination of the Pick Force
The cutter head operates with a composite motion, which includes rotating and swinging, and also through the conical tip of picks to break up the rocks. Due to the complex geological conditions in tunneling, it is very difficult to measure the cutting force imposed on the cutter head. Hence, the data obtained from the laboratory rock cutting tests with a single pick are usually used as the basis for calculating the forces of each cutter head pick. The prediction formulas of the pick forces obtained from the laboratory tests are more reasonable and reliable than those obtained in the theoretical and numerical simulation methods. In the paper, the cutting force acting on the pick is shown in Figure  1. The parameters present that d is the depth of cut, δ is the angle of attack, and n is the speed. FC is the cutting force, which is opposite to the direction of motion acting on the tip of the pick, and it is to separate the cuttings along the direction of the pick from the rock mass and forms a trough. FN is the normal force whose component force is perpendicular to the direction of motion of the coal rock, and it is to let the pick compress the rock and form the radiant point of the crack extension. Lastly, FS is the lateral force, from the normal force to the cutting direction according to the right-hand rule, and it is used to promote the rotation of the pick and forms a uniform wear.

The Introduction of the Experimental Machine
The test is performed using a full-scale rotary cutting machine, which is currently one of the most advanced single-pick rotary cutting machines, at the National Engineering Laboratory of Coal Mining Machinery and Equipment in China. This machine consists of a test platform, a control system, a test system, and data analysis software. The structure of the test platform is shown in Figure  2. The device can sensitively monitor the cutting load, vibration, temperature, and dust amount during tests. The largest size of the sample rock that it can accommodate is 1400 × 800 × 600 mm. The rock is clamped with a special fixture to prevent it from moving under large impacts. Conical picks with wolfram carbide (WC) cemented carbide pick tip are used in the test, and the parameters are that the tip diameter is 25 mm, pick angle is 80°, external elongation of the pick is 80 mm, the handle width is 34 mm, and the handle length is 78 mm. In the whole test process, the cutting speed is 1.47

The Introduction of the Experimental Machine
The test is performed using a full-scale rotary cutting machine, which is currently one of the most advanced single-pick rotary cutting machines, at the National Engineering Laboratory of Coal Mining Machinery and Equipment in China. This machine consists of a test platform, a control system, a test system, and data analysis software. The structure of the test platform is shown in Figure 2. The device can sensitively monitor the cutting load, vibration, temperature, and dust amount during tests. The largest size of the sample rock that it can accommodate is 1400 × 800 × 600 mm. The rock is clamped with a special fixture to prevent it from moving under large impacts. Conical picks with wolfram carbide (WC) cemented carbide pick tip are used in the test, and the parameters are that the tip diameter is 25 mm, pick angle is 80 • , external elongation of the pick is 80 mm, the handle width Appl. Sci. 2019, 9, 5410 4 of 18 is 34 mm, and the handle length is 78 mm. In the whole test process, the cutting speed is 1.47 m/s, the cutting depth is 4 mm, the tool spacing is 12 mm, and the strike angle is 50 • . The pick was used to cut the rock from the left to right end with a total of seventy cutting cycles under the same line spacing.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 17 m/s, the cutting depth is 4 mm, the tool spacing is 12 mm, and the strike angle is 50°. The pick was used to cut the rock from the left to right end with a total of seventy cutting cycles under the same line spacing.

The Mean Cutting Force of Conical Picks
Through the experiment, the normal force, the cutting force, and the lateral force from the experimental data are shown in Figure 3. It is easily found that the normal force is the largest among the components of the cutting load, the cutting force is the second, and the lateral force is the smallest. In this paper, the new pick was studied, so the first nine cutting test data were extracted to obtain the average force (Table 1). In our work, we set the normal force at 5.5 kN, the cutting force at 0.6 kN, and the lateral force at 0.002 kN.

The Mean Cutting Force of Conical Picks
Through the experiment, the normal force, the cutting force, and the lateral force from the experimental data are shown in Figure 3. It is easily found that the normal force is the largest among the components of the cutting load, the cutting force is the second, and the lateral force is the smallest. In this paper, the new pick was studied, so the first nine cutting test data were extracted to obtain the average force (Table 1). In our work, we set the normal force at 5.5 kN, the cutting force at 0.6 kN, and the lateral force at 0.002 kN.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 17 m/s, the cutting depth is 4 mm, the tool spacing is 12 mm, and the strike angle is 50°. The pick was used to cut the rock from the left to right end with a total of seventy cutting cycles under the same line spacing.

The Mean Cutting Force of Conical Picks
Through the experiment, the normal force, the cutting force, and the lateral force from the experimental data are shown in Figure 3. It is easily found that the normal force is the largest among the components of the cutting load, the cutting force is the second, and the lateral force is the smallest. In this paper, the new pick was studied, so the first nine cutting test data were extracted to obtain the average force (Table 1). In our work, we set the normal force at 5.5 kN, the cutting force at 0.6 kN, and the lateral force at 0.002 kN.

Create the Edge-Based Smoothing Domains for Conical Picks
Before constructing the ES-FEM model, we uses the typical linear element T4 for 3D problems as the background mesh, which can automatically mesh and is very suitable for complex problem domains. Based on the T4 mesh, the edge-based smoothing domain can be created following the algorithm introduced below. In order to illustrate the process clearly, a simple ES smoothing domain is plotted in Figure 4 and the descriptions of elements and faces are listed in Table 2. Nodes n 1, n 2 , n 3 , n 4 made up the element 1 (e 1 ); and nodes n 1, n 2 , n 3 , n 5 made up the element 2 (e 2 ). Face1 (f 1 ) contained nodes n 1, n 3 , n 4 ; Face2 (f 2 ) contained nodes n 1, n 2 , n 3 ; Face3 (f 3 ) contained nodes n 1, n 3 , n 5 . In fact, ES smoothing domain relies on each edge so the number of smoothing domains equals the number of edges in the background. The smoothed domain will adopt the edges on the element as the indexes. Take the element in the Figure 4, for example, there are two elements e 1 and e 2 . We will discuss how to build the smoothing domain of edge n 1 n 3 . The edge of n 1 n 3 is an internal edge which is related to the different elements. Let us focus on the element e 1 , because the process is the same as that in the element e 2 . In this element, the edge locates on the faces f 1 , f 2 , and f 3 . Then connect the center of one face, the centroid of this element, and one vertex of the edge to form a surface-segment of the smoothing domain. Repeat the process and obtain six-surface segments for one element, n 3 c f1 c e1 , n 1 c f1 c e1 , n 3 c f2 c e1 , n 1 c f2 c e1 . Similarly, we can build other six-surface segments for the element e 2 . It is clearly seen that the 3D-edge-based smoothed domain is enclosed by these surfaces. Hence, the ID of the smoothing domain is the same as the ID of the edge, which is surrounded by the smoothing domain. In Figure 4, one T4 element is separated by six smoothing domains according to the edges of this element. Table 2. The details of the grid in Figure 4.

Element/Face The Nodes of Element/Face The Centroid/The Center
Element 1 (e 1 ) n 1 n 2 n 3 n 4 C e1 Element 2 (e 2 ) n 1 n 2 n 3 n 5 C e2 Face1 (f 1 ) n 1 n 3 n 4 C f1 Face2 (f 2 ) n 1 n 2 n 3 C f2 Face3 (f 3 ) n 1 n 3 n 5 C f3 smoothing domain. Repeat the process and obtain six-surface segments for one element, n3cf1ce1, n1cf1ce1, n3cf2ce1, n1cf2ce1. Similarly, we can build other six-surface segments for the element e2. It is clearly seen that the 3D-edge-based smoothed domain is enclosed by these surfaces. Hence, the ID of the smoothing domain is the same as the ID of the edge, which is surrounded by the smoothing domain. In Figure 4, one T4 element is separated by six smoothing domains according to the edges of this element.  Table 2. The details of the grid in Figure 4. Through the above introduction, the constructing algorithm needs the connectivity of element, face, edge, and node. Therefore, for creating a smoothing domain quickly, we first need to get the connectivity of nodes separately with elements, faces, and edges; the connectivity of edges separately with elements, faces, and nodes; the connectivity of faces separately with elements, edges, and nodes; and the connectivity of elements separately with faces, edges, and nodes.

Construct the Smoothed Strain Field of ES-FEM
In an ES-FEM model, each of the smoothing domains Ω S k , k = 1, 2, . . . , N S in general consist N S ≥ 1 of sub-smoothing domain Ω S k,q , which is a part of the elements where the corresponding edge locates. The smoothed strain field for one smoothing domain can be obtained by boundary flux as follows: where Γ S k,q is the boundary of the qth sub-smoothing domain of smoothing domain Ω S k,q , Γ S k is the boundary of the smoothing domain Ω S k , V e j , and V e j is the volume of the jth tetrahedral element around edge k. N S is the number of the sub-smoothing domain. u(x) is the displacement matrix, which has the form of where d the displacement matrix of all nodes in ES-FEM and N(x) is the shape function matrix. L n (x) is the matrix of components of the outward normal vector on the boundary Γ S k,q , hence the form as:

Construct the ES Smoothed Strain-Displacement Matrix
Substituting Equation (2) into Equation (1), the smoothed strain is: whereB is the global "smoothed strain-displacement" matrix, which is due to the shape function matrix N(x). For ES-FEM using the four-node tetrahedron element, as shown Figure 4, the set of nodes are {n 1 , n 2 , n 3 , n 4 , n 5 } for edge k. The smoothed strain-displacement matrix B I is evaluated using: Therefore, for the internal edge of four-node tetrahedron element,B is a matrix of 6 × 15, for the border edge,B is a 6 × 12 matrix. Besides, d is a vector containing 12 elements. The Equation (6) can be further simplified to a summation form by using the Gauss quadrature technique: where n s Γ is the total number of boundary segment Γ S E,p ∈ Γ S E and x G p is the Gauss point of boundary segment of Γ S E,p ,whose area and outward unit normally are donated as A P and n h,P , respectively. Thus, B I can be given in detail as follows: Appl. Sci. 2019, 9, 5410 8 of 18 The general formulation for computing the smoothed stiffness matrix is with entries where c is the matrix of material constants.

Results
In the section, in order to illustrate the efficiency and accuracy of the present ES-FEM algorithm for conical pick problems, the numerical results of EF-FEM are listed and are compared with the solutions of other S-FEM models and FEM in stress, strain energy, and displacement.
The geometrical parameters of the practical conical pick used in our experiment are given in Figure 5. The material parameters in the experiment used are: E = 650, 000 MPa and υ = 0.3. At the same time, the normal force of F z = 5.5 kN, the cutting force of F y = 0.6 kN, and the lateral force of F x = 0.002 kN are imposed on the top surface of the conical pick, and the force area is three-quarters of the top area shown in Figure 6. First, the numbers and area of the forced elements are obtained based on the ES-FEM mesh; and then the forces of one element are imposed on three nodes through the face normal algorithm. Besides, we fully fix the bottom of the conical pick. where c is the matrix of material constants.

Results
In the section, in order to illustrate the efficiency and accuracy of the present ES-FEM algorithm for conical pick problems, the numerical results of EF-FEM are listed and are compared with the solutions of other S-FEM models and FEM in stress, strain energy, and displacement.
The geometrical parameters of the practical conical pick used in our experiment are given in   For practical problems, it is very difficult to find the analytical solution. Hence, we calculate the For practical problems, it is very difficult to find the analytical solution. Hence, we calculate the reference solutions of the problem using the standard FEM with the mesh of 65,686 T4 elements, and 197,058 degrees of freedom.
In Figure 7a, the conical pick is divided by tetrahedron elements as the background mesh. We use the connectivity of the background mesh listed above the tables and the algorithm mentioned in Section 3.1 to create an ES mesh for our numerical simulation (Figure 7b). From this figure, we can observe the smoothing domains are surrounding the edges in the background mesh. Seven background meshes are applied to different methods such as FEM, ES-FEM, FS-FEM, NS-FEM, which have the same numbers of nodes, elements, and degrees of freedom. The detailed information of these seven meshes is shown in Table 3.  For practical problems, it is very difficult to find the analytical solution. Hence, we calculate the reference solutions of the problem using the standard FEM with the mesh of 65,686 T4 elements, and 197,058 degrees of freedom.
In Figure 7a, the conical pick is divided by tetrahedron elements as the background mesh. We use the connectivity of the background mesh listed above the tables and the algorithm mentioned in  Table 3.

The Stress Solutions
Von-Mises stress, which is the fourth strength theory, can clearly describe the change of a result in the whole model through using stress contours to represent the stress distribution, and makes sure designers can quickly determine the most dangerous area in the model. Thus, Von-Mises stress is used to describe stress distribution in post-processing of mature FEM software. In order to verify the validity and accuracy of our method, we compared the Von-Mises stress value of ES-FEM with the reference value and FEM analysis results. Von-Mises stress is computed using the stress components in the form of:

The Stress Solutions
Von-Mises stress, which is the fourth strength theory, can clearly describe the change of a result in the whole model through using stress contours to represent the stress distribution, and makes sure designers can quickly determine the most dangerous area in the model. Thus, Von-Mises stress is used to describe stress distribution in post-processing of mature FEM software. In order to verify the validity and accuracy of our method, we compared the Von-Mises stress value of ES-FEM with the reference value and FEM analysis results. Von-Mises stress is computed using the stress components in the form of: Figure 8 is the reference value of Von-Mises stress σ v using the standard FEM with 197,058 degrees of freedom. Then we plot the distribution of von-Mises stress using the present ES-FEM and FEM with the mesh having 6413 T4 elements in Figure 9. It is easily seen that the Von-Mises stress σ v distributions of the ES-FEM model are similar to that of FEM. Next, the distribution of Von-Mises stress along the generator of the tip which includes the max Von-Mises stress using different methods for the conical pick are computed and shown in Figure 11. We first find that the maximum von-Mises in the Figure 11 was about 7.5 mm away from the apex of conical pick, and the von-Mises would be steadily reducing with the distance increasing from the apex. In addition, it is shown that the ES-FEM Von-Mises are closer to the reference solutions than Next, the distribution of Von-Mises stress along the generator of the tip which includes the max Von-Mises stress using different methods for the conical pick are computed and shown in Figure 11. We first find that the maximum von-Mises in the Figure 11 was about 7.5 mm away from the apex of conical pick, and the von-Mises would be steadily reducing with the distance increasing from the apex. In addition, it is shown that the ES-FEM Von-Mises are closer to the reference solutions than that of FEM. It was further proof that our ES-FEM for the practical conical pick is more accurate and   In the section, we studied the performance of the ES-FEM on the strain energy for this practical problem. First, we calculated the reference solution of strain energy using the finer mesh mentioned above, which is 2.4240 N mm ⋅ . Then, the numerical results of strain energy using different methods are presented in Table 4. The error of strain energy obtained using ES-FEM is smallest no matter which mesh is used. The according convergence curves of strain energy is plotted in Figure 12 against the degrees of freedom (DOFs). It can be found that the strain energy of ES-FEM are just a little bigger than the reference solution and quickly converge to it with the increase of DOFs. In contrast, the strain energy solution of FEM is much lower than the reference solution and it converges very slowly. It is noted that, although the ES-FEM solution fluctuates a little bit when using the coarse meshes, the value is more close to the reference solution than other methods. Besides, it is observed that the upper bounds of NS-FEM and the lower bounds of FEM still remain the same in this simulation. Next, the distribution of Von-Mises stress along the generator of the tip which includes the max Von-Mises stress using different methods for the conical pick are computed and shown in Figure 11. We first find that the maximum von-Mises in the Figure 11 was about 7.5 mm away from the apex of conical pick, and the von-Mises would be steadily reducing with the distance increasing from the apex. In addition, it is shown that the ES-FEM Von-Mises are closer to the reference solutions than that of FEM. It was further proof that our ES-FEM for the practical conical pick is more accurate and efficient among other S-FEM methods and FEM.  In the section, we studied the performance of the ES-FEM on the strain energy for this practical problem. First, we calculated the reference solution of strain energy using the finer mesh mentioned above, which is 2.4240 N mm ⋅ . Then, the numerical results of strain energy using different methods are presented in Table 4. The error of strain energy obtained using ES-FEM is smallest no matter which mesh is used. The according convergence curves of strain energy is plotted in Figure 12 against the degrees of freedom (DOFs). It can be found that the strain energy of ES-FEM are just a little bigger than the reference solution and quickly converge to it with the increase of DOFs. In contrast, the strain

The Strain Energy Solutions
In the section, we studied the performance of the ES-FEM on the strain energy for this practical problem. First, we calculated the reference solution of strain energy using the finer mesh mentioned above, which is 2.4240 N · mm. Then, the numerical results of strain energy using different methods are presented in Table 4. The error of strain energy obtained using ES-FEM is smallest no matter which mesh is used. The according convergence curves of strain energy is plotted in Figure 12 against the degrees of freedom (DOFs). It can be found that the strain energy of ES-FEM are just a little bigger than the reference solution and quickly converge to it with the increase of DOFs. In contrast, the strain energy solution of FEM is much lower than the reference solution and it converges very slowly. It is noted that, although the ES-FEM solution fluctuates a little bit when using the coarse meshes, the value is more close to the reference solution than other methods. Besides, it is observed that the upper bounds of NS-FEM and the lower bounds of FEM still remain the same in this simulation. Table 4. The error of strain energy obtained using different methods compared to the reference solution (%).   Next, in order to verify the accuracy and validity, the error norm of the strain energy of ES-FEM was calculated and compared with that of FEM, NS-FEM, and FS-FEM models. Errors of these numerical methods are analyzed using energy norms. The error norm of strain energy is calculated by:

Type of Mesh FEM ES-FEM FS-FEM NS-FEM Reference Solution
where ε is the reference solution for the strain energy, and  ε is the numerical solution for the strain energy using a numerical model. Table 5 lists seven sets of data calculated using different methods for seven meshes. It is shown that the error of ES-FEM is about 0.217 against 0.534 of FEM, which illustrates again that ES-FEM stands out clearly.  Next, in order to verify the accuracy and validity, the error norm of the strain energy of ES-FEM was calculated and compared with that of FEM, NS-FEM, and FS-FEM models. Errors of these numerical methods are analyzed using energy norms. The error norm of strain energy is calculated by: where ε is the reference solution for the strain energy, and ε is the numerical solution for the strain energy using a numerical model. Table 5 lists seven sets of data calculated using different methods for seven meshes. It is shown that the error of ES-FEM is about 0.217 against 0.534 of FEM, which illustrates again that ES-FEM stands out clearly.

The Nodes Displacement Solutions
In this part, we will examine the displacement solution for the conical pick. The reference value of nodes' displacements along y-axis u 2 and z-axis u 3 are plotted in Figure 13 and we draw the approximate displacement solutions along the y-axis u 2 and z-axis u 3 using ES-FEM in Figure 14.

The Nodes Displacement Solutions
In this part, we will examine the displacement solution for the conical pick. The reference value of nodes' displacements along y-axis 2 u and z-axis 3 u are plotted in Figure 13 and we draw the approximate displacement solutions along the y-axis 2 u and z-axis 3 u using ES-FEM in Figure 14.
From these figures it is distinctly observed that the displacement solutions of ES-FEM are reasonable and accurate.
(a) (b) Figure 13. The reference value of nodes' displacements along (a) the y-axis 2 u and (b) z-axis 3 u .
(a) (b) Figure 13. The reference value of nodes' displacements along (a) the y-axis u 2 and (b) z-axis u 3 . Table 6 shows the comparison of displacements of ES-FEM with FEM, NS-FEM, and FS-FEM. Figure 15 is plotted by the distance from the tip of the conical pick on the horizontal axis and the node's displacement on the vertical. The comparison is fair and rigorous as long as the same distribution of nodes is used. Overall, the displacements of ES-FEM are closer to the reference solutions than those of FEM. This also verifies that, regardless of strain and displacement, ES-FEM is the most accurate and stable method compared with FEM and other S-FEMs.   In order to further verify the accuracy and stability, the error of displacement of the ES-FEM model was conducted. Displacement errors and convergence rates of these numerical methods are analyzed using displacement norms. The displacement norm is defined as: where u is the reference solution of the displacements, and u is the numerical solution of the displacements using a numerical model. In order to further verify the accuracy and stability, the error of displacement of the ES-FEM model was conducted. Displacement errors and convergence rates of these numerical methods are analyzed using displacement norms. The displacement norm is defined as: where u is the reference solution of the displacements, and u  is the numerical solution of the displacements using a numerical model. Table 7

The Computational Efficiency of Solutions
In the final part, we discuss the computational efficiency of the methods for this example. The efficiency is estimated using the following formulation: where r e is the strain energy norm error of different methods compared with FEM, and r t represents CPU time compared with FEM. Tables 8 and 9 list the strain energy norm error and CPU time obtained using four methods with mesh 4 and mesh 5. It is easily found that the solution error of ES-FEM is the smallest even though the CPU time is one of the highest. Taking the two facts into consideration, the computational efficiency objectively describes the performance of different numerical methods. Hence, ES-FEM, whose computational efficiency is about 2.88 times that of FEM, is the most efficient one. Similarly, the computational efficiency of FS-FEM has been improved in comparison with FEM. Thus, the computational efficiency of ES-FEM and FS-FEM can be significantly improved when CPU time is taken into account.

Conclusions
In the paper, the efficient and ultra-accurate ES-FEM is used to simulate the cutting process of conical picks. Through analyzing the cutting characteristics of a conical pick, the normal force, the cutting force, and the lateral force are obtained from the experimental data using a full-scale rotary cutting machine as the boundary conditions in our simulation. Firstly, we present an efficient algorithm for creating edge-based smoothing domains, which generates the twelve connectivity lists, including the connectivity of face-node, face-edge, face-element, node-face, edge-face, element-face, edge-node, edge-element, element-edge, node-edge, node-element, and element-node. Then the gradient smoothing technique is employed to construct the smoothed strain gradient and displacement gradient fields. Next, the smoothed system equations are set up following the theory of S-FEM. Finally, displacement, strain energy, and stress are calculated using the ES-FEM for the conical pick. At the same time, various comparisons between ES-FEM and other methods are made to demonstrate the accuracy and effectiveness of the ES-FEM model. The concluding remark was drawn as follows: 1.
The procedure of the ES smoothing domain presented in this paper is very easy and performs very well for the complex structured parts of coal mining machines.

2.
The present method showed a better accuracy and convergence rate than FEM for this practical problem, which illustrate S-FEM is suitable for practical engineering problems.