Study on an Analytical Solution to the Mechanical Response of Soil Cut by a Bionic Force Enhancement Grouser

: With the depletion of land mineral resources, people have turned their attention to the sea. As an important part of deep-sea mining systems, the technical research and development of ore collectors has always been a difﬁcult problem in many countries. According to the characteristic that buffalo hoof is suitable for walking on soft soil, a kind of bionic grouser for a deep-sea mining vehicle is designed in this paper. Through the optimization of Rankine’s passive earth pressure theory, the formula for calculating the tractive force of the grouser is obtained. The accuracy of the analytical solution is veriﬁed by ﬁnite element simulation, and the force enhancement mechanism of the bionic grouser is revealed. The results show that the design of the bionic grouser has a signiﬁcant effect on the improvement of tractive force, and the tractive force of the No. 1 bionic grouser is 17.52% higher than that of the straight grouser. On this basis, the geometric parameters of the bionic grouser proﬁle are optimized. The results show that when L is 0 mm and R is 183 mm, the force enhancement effect reaches the maximum of 27%, which provides a design basis for optimizing the grouser and improving the mining efﬁciency of the collector. Author Contributions: Conceptualization and methodology, W.M., Y.C., W.Z. and W.Q.; software, Y.C. and W.M.; validation, Y.C. and W.Z.; investigation, Y.C. and W.Q.; resources, W.M. and W.Z.; theoretical derivation, Y.C. and W.Z.; data curation, Y.C. and W.M.; writing—original draft preparation, W.M., Y.C. and W.Q.; writing—review and editing, W.M., Y.C. and W.Q.; supervision W.M. and Y.C.; administration, W.M. and W.Z.


Introduction
With the rapid economic development of countries worldwide, the demand for mineral resources is becoming more and more exuberant, which gradually makes it more and more difficult for land resources to meet the needs of development [1]. To maintain sustained economic growth, people turn their attention to the ocean [2,3]. Marine mineral resources are prosperous, especially submarine metal nodule deposits, which are extremely large, in which rare metal nodule minerals (such as manganese metal nodules) and gas hydrates are extremely abundant [4][5][6], but most of these minerals are distributed on the seafloor below 4000 m. Therefore, mining is challenging [7,8]. The complex deep-sea environment puts higher requirements on the deep-sea mining system. After the successful sea test of "Kunlong 500" developed by the team in the Pacific Ocean, the hydraulic mining system is considered the most balanced commercial mining system [9,10].
As a vital part of the deep-sea mining system, the collector shown in Figure 1 is a critical point. The characteristics required for the regular operation of the collector include the ore collection stability, energy efficiency, durability, maintainability, and controllability [11,12]. The collector's task is to efficiently walk and mine on the soft seafloor subsoil, collect, clean, and transport the ore to the ore lifting system [13,14], without causing excessive damage to the marine environment. Due to the low shear strength and obvious rheological behavior of deep-sea sediments, it is easy to slip between the grouser so that displacement. Yang et al. [19,20] derived the horizontal and vertical resistance of the tracked vehicle and used MATLAB to determine the optimal driving resistance by analyzing the vehicle parameters and soil parameters. Based on the constitutive model of the elastic model and the Mohr-Coulomb model, Abaqus was used to simulate the grouser shear soil, and the simulation results were in good agreement with the experimental results. Shin et al. [21] found that the total thrust of the tracked vehicle is slightly lower than the product of the number of single grousers and the number of grousers and carried out a series of track tests using the three-track system. The experimental results showed that the soil thrust of the rear grouser is small due to the shielding effect of the front grouser, while the influence of the rear grouser on the front is negligible. Bionics is an ancient and young discipline. Animals and plants, in order to adapt to the environment, through millions or even thousands of years of evolution, have structures that are often much more exquisite than human machinery. By studying the mechanism and anatomical structure of biology, scientists apply the developed bionic structure to engineering [22][23][24], medicine [25][26][27], and even computer fields [28], bringing new technological innovations to various industries. Sun [29] created a new type of micro air vehicle by biomimetic design of the origami wings of rhino beetles, which can restore a stable flight state from a flight collision. Wu et al. [30] designed a bionic plate by studying the microstructure of fish scales, and compared it with a smooth plate under the flow of 0.66 m/s, and found that the maximum drag reduction rate reached 2.8%.
The design of the collector grouser is crucial for improving the walking stability and mining efficiency of a collector [31,32], and therefore, improving the common phenomena In recent years, scientists have done much research on the tractive force of grousers in the process of walking on the soft subsoil. Li et al. [16] studied the tractive force of grousers under different sand content sediments and found that the higher the sand content under the same slip ratio, the greater the tractive force, and the slip rate should be controlled below 6% in order to avoid the decrease of tractive force. Using the traditional Mohr-Coulomb failure criterion and the shear speed of 15 mm/s, Janarthanan et al. [17] carried out experiments and simulations in the deep-sea soft soil environment and studied the influence of different shapes of grouser on the traction force of the collector. Li et al. [18] carried out a large number of experiments on π-shaped, T-shaped, and V-shaped grousers, and found that the traction force of the π-shaped crawler is more sensitive to width and the T-shaped one is more sensitive to height, while the structural parameters of the V-shaped grouser have little influence on it, but it has a minor vertical subsidence displacement. Yang et al. [19,20] derived the horizontal and vertical resistance of the tracked vehicle and used MATLAB to determine the optimal driving resistance by analyzing the vehicle parameters and soil parameters. Based on the constitutive model of the elastic model and the Mohr-Coulomb model, Abaqus was used to simulate the grouser shear soil, and the simulation results were in good agreement with the experimental results. Shin et al. [21] found that the total thrust of the tracked vehicle is slightly lower than the product of the number of single grousers and the number of grousers and carried out a series of track tests using the three-track system. The experimental results showed that the soil thrust of the rear grouser is small due to the shielding effect of the front grouser, while the influence of the rear grouser on the front is negligible.
Bionics is an ancient and young discipline. Animals and plants, in order to adapt to the environment, through millions or even thousands of years of evolution, have structures that are often much more exquisite than human machinery. By studying the mechanism and anatomical structure of biology, scientists apply the developed bionic structure to engineering [22][23][24], medicine [25][26][27], and even computer fields [28], bringing new technological innovations to various industries. Sun [29] created a new type of micro air vehicle by biomimetic design of the origami wings of rhino beetles, which can restore a stable flight state from a flight collision. Wu et al. [30] designed a bionic plate by studying the microstructure of fish scales, and compared it with a smooth plate under the flow of 0.66 m/s, and found that the maximum drag reduction rate reached 2.8%.
The design of the collector grouser is crucial for improving the walking stability and mining efficiency of a collector [31,32], and therefore, improving the common phenomena of skidding and subsidence in the operation of a deep-sea collector. In this paper, the bionic force enhancement design of a collector grouser is studied, which provides a necessary theoretical basis for its optimal design.

Bionic Design of the Grouser
Based on the unique shape curve of the buffalo hoof, the characteristics of the grouser that can improve the traction capacity of the collector track are analyzed. Through observation, it can be seen that the shape of the front surface of the buffalo's rear hoof can make the buffalo walk stably in the soft clay, which is an important reason for the buffalo to gain propulsion. In this study, bionic grousers were designed according to the structural characteristics of the front surface of buffalo hoof entry. Figure 2a shows the apparent morphology of the buffalo hoof, and Figure 2b shows the cross-sectional feature diagram of the bionic grouser. Where L is the distance from the curvature change point to the top of the grouser, and R is the curvature radius of the arc part. The obvious point of curvature change can be observed from the figure, which is an essential factor affecting the propulsion force and walking resistance. of skidding and subsidence in the operation of a deep-sea collector. In this paper, the bionic force enhancement design of a collector grouser is studied, which provides a necessary theoretical basis for its optimal design.

Bionic Design of the Grouser
Based on the unique shape curve of the buffalo hoof, the characteristics of the grouser that can improve the traction capacity of the collector track are analyzed. Through observation, it can be seen that the shape of the front surface of the buffalo's rear hoof can make the buffalo walk stably in the soft clay, which is an important reason for the buffalo to gain propulsion. In this study, bionic grousers were designed according to the structural characteristics of the front surface of buffalo hoof entry. Figure 2a shows the apparent morphology of the buffalo hoof, and Figure 2b shows the cross-sectional feature diagram of the bionic grouser. Where L is the distance from the curvature change point to the top of the grouser, and R is the curvature radius of the arc part. The obvious point of curvature change can be observed from the figure, which is an essential factor affecting the propulsion force and walking resistance.

Hypothetical Conditions for the Derivation of Analytical Solutions
Because the classical Rankine earth pressure theory has a wide range of applications and the calculation process is simple and easy to understand, it has been widely popularized in practical engineering. However, Rankine's earth pressure theory has very obvious limitations: the first condition of Rankine's earth pressure theory is to assume that the back of the wall is vertical and smooth and the filling behind the wall is horizontal, but there is often an inclination angle on the back of the wall in practical engineering. This requires us to adopt a new method to optimize Rankine earth pressure. Now the passive earth pressure calculation model of the bionic grouser is assumed as follows: (1) grousers are rigid and curved; (2) the surface of the deep-sea bottom soil is horizontal, and within the range of mining depth, it is homogeneous; (3) the friction in the direction of grouser thickness is not considered.

Overall Force Analysis
As shown in Figure 3, with the shearing of the grouser plate, the soil near the grouser is squeezed. According to Rankine's earth pressure theory, it is considered that the soil is in a Rankine passive earth pressure state in this area, and the soil after the grouser is shearing along the slip plane BC. The normal pressure Fp acting on the grouser is the

Hypothetical Conditions for the Derivation of Analytical Solutions
Because the classical Rankine earth pressure theory has a wide range of applications and the calculation process is simple and easy to understand, it has been widely popularized in practical engineering. However, Rankine's earth pressure theory has very obvious limitations: the first condition of Rankine's earth pressure theory is to assume that the back of the wall is vertical and smooth and the filling behind the wall is horizontal, but there is often an inclination angle on the back of the wall in practical engineering. This requires us to adopt a new method to optimize Rankine earth pressure. Now the passive earth pressure calculation model of the bionic grouser is assumed as follows: (1) grousers are rigid and curved; (2) the surface of the deep-sea bottom soil is horizontal, and within the range of mining depth, it is homogeneous; (3) the friction in the direction of grouser thickness is not considered.

Overall Force Analysis
As shown in Figure 3, with the shearing of the grouser plate, the soil near the grouser is squeezed. According to Rankine's earth pressure theory, it is considered that the soil is in a Rankine passive earth pressure state in this area, and the soil after the grouser is shearing along the slip plane BC. The normal pressure Fp acting on the grouser is the passive earth pressure, which varies with the volume of the compressed soil on the surface of the grouser. The passive earth pressure of cohesive soil is also composed of another part, that is, the earth pressure Fc caused by cohesion c is positive and does not change with depth.
The traction formula of the grouser is: where the Fp is the passive earth pressure of the grouser, and Fc is: First of all, the compression area of the surface is analyzed, and the value of passive earth pressure Fp of the grouser surface is calculated. At this time, the passive earth pressure, its sliding surface, and its stress diagram are shown in Figure  passive earth pressure, which varies with the volume of the compressed soil on the surface of the grouser. The passive earth pressure of cohesive soil is also composed of another part, that is, the earth pressure Fc caused by cohesion c is positive and does not change with depth. The traction formula of the grouser is: where the Fp is the passive earth pressure of the grouser, and Fc is: First of all, the compression area of the surface is analyzed, and the value of passive earth pressure Fp of the grouser surface is calculated. At this time, the passive earth pressure, its sliding surface, and its stress diagram are shown in Figure 3: The closed area where the ABC is located is taken as the isolation unit, and the selfweight of the isolation body under passive earth pressure is as follows: where γ is the effective unit weight of deep-sea soil, b is the depth of grouser, and SABC is the area of ABC. The calculation formula of passive earth pressure is: where pp(z) is the value of earth pressure per unit depth at depth z.

Force Analysis of the Microelement
The force acting on the microelement under passive earth pressure is shown in the Figure 4. Macroscopically, the static W, R, and f of the isolator are dW, q, and dFp, and the contact force between soil layers is expressed by the uniformly distributed load q − dq and q of the surface layer of the microelement, as shown in Figure 5. The closed area where the ABC is located is taken as the isolation unit, and the self-weight of the isolation body under passive earth pressure is as follows: where γ is the effective unit weight of deep-sea soil, b is the depth of grouser, and S ABC is the area of ABC. The calculation formula of passive earth pressure is: where p p (z) is the value of earth pressure per unit depth at depth z.

Force Analysis of the Microelement
The force acting on the microelement under passive earth pressure is shown in the Figure 4. Macroscopically, the static W, R, and f of the isolator are dW, q, and dFp, and the contact force between soil layers is expressed by the uniformly distributed load q − dq and q of the surface layer of the microelement, as shown in Figure 5. Firstly, the soil weight dW of the microelement is calculated. dW is related to the M1N1N2M2 area of trapezoidal microelements. The formula for calculating the height of the trapezoid is as follows: where θ is the dip angle between the slip surface and horizontal plane, and ψ is the angle between a tangent and the vertical direction of the grouser surface corresponding to microelement dz. ψ depends on the characteristics of the grouser surface, and its calculation formula is as follows: The length of the lower bottom edge and the upper bottom edge of the trapezoid is calculated as: In the joint Formula (5), the trapezoidal area can be expressed as:  Firstly, the soil weight dW of the microelement is calculated. dW is related to the M1N1N2M2 area of trapezoidal microelements. The formula for calculating the height of the trapezoid is as follows: where θ is the dip angle between the slip surface and horizontal plane, and ψ is the angle between a tangent and the vertical direction of the grouser surface corresponding to microelement dz. ψ depends on the characteristics of the grouser surface, and its calculation formula is as follows: The length of the lower bottom edge and the upper bottom edge of the trapezoid is calculated as: In the joint Formula (5), the trapezoidal area can be expressed as: Firstly, the soil weight dW of the microelement is calculated. dW is related to the M 1 N 1 N 2 M 2 area of trapezoidal microelements. The formula for calculating the height of the trapezoid is as follows: where θ is the dip angle between the slip surface and horizontal plane, and ψ is the angle between a tangent and the vertical direction of the grouser surface corresponding to microelement dz. ψ depends on the characteristics of the grouser surface, and its calculation formula is as follows: The length of the lower bottom edge and the upper bottom edge of the trapezoid is calculated as: In the joint Formula (5), the trapezoidal area can be expressed as: Formula (6) can be substituted for (1) to derive the soil weight of the microelement as follows: Based on the static analysis of the microelement, the equilibrium equation in the limit equilibrium state is as follows: Eliminate the (q − dq) item in the above formula and get: Insert Formula (7) into Formula (9) to obtain: Among them, the specific expression of fp (θ) is: As a result, the passive earth pressure resultant force Fp of the grouser is expressed in the form of definite integral as follows: If θ = π 4 − ϕ 2 is taken, and placed in (12), then: Formula (14) is discretized by the Simpson formula, and the analytical solution of passive earth pressure is obtained as follows:

Theoretical Data
By substituting different R and L values in Formula (16), the maximum traction force of the bionic grouser is obtained, and the results are shown in Table 1. As can be seen from the table, with the decrease of R and L, the traction force increases obviously, in which the bionic grouser, with the distance between grouser L of 30 mm and the radius of curvature R of 500 mm, has the highest tractive force. Compared with the straight plate grouser, the tractive force increases by 9.68%, and the force-increasing effect is obvious.

Simulation Model
The finite element model of a grouser cutting deep-sea bottom soil is shown in Figure 6. It is assumed that the grouser is a curved rigid body, and the velocity in the x-direction is applied on the rigid body, and a reference point is set at the upper end of the grouser, in which its traction force is easy to extract. The fixed support is applied to the bottom and left side of the soil, and the upper and right sides are free surfaces. The time length is set to 4 s and the speed to 15 mm/s [14]. Considering that the height of the bionic grouser of the second generation deep-sea collector is 130 mm, the design parameters L and R of the bionic grouser are 30 mm and 500 mm, respectively. This paper studies the traction characteristics of different bionic grousers through the numerical simulation of cutting deep-sea subsoil with different bionic grousers.

Simulation Model
The finite element model of a grouser cutting deep-sea bottom soil is shown in Figure  6. It is assumed that the grouser is a curved rigid body, and the velocity in the x-direction is applied on the rigid body, and a reference point is set at the upper end of the grouser, in which its traction force is easy to extract. The fixed support is applied to the bottom and left side of the soil, and the upper and right sides are free surfaces. The time length is set to 4 s and the speed to 15 mm/s [14]. Considering that the height of the bionic grouser of the second generation deep-sea collector is 130 mm, the design parameters L and R of the bionic grouser are 30 mm and 500 mm, respectively. This paper studies the traction characteristics of different bionic grousers through the numerical simulation of cutting deepsea subsoil with different bionic grousers.

Deep-Sea Soil Model
Through the shear strength test of deep-sea soft subsoil, many researchers, through studies on the physical and mechanical properties of deep-sea subsoil, found that its shear strength increases with the increase of depth and tends to be stable after a certain depth. It has been shown that the thin and soft subsoil has the bearing capacity to make the collector walk freely. The rate-dependent characteristics of deep-sea bottom soil mean that

Deep-Sea Soil Model
Through the shear strength test of deep-sea soft subsoil, many researchers, through studies on the physical and mechanical properties of deep-sea subsoil, found that its shear strength increases with the increase of depth and tends to be stable after a certain depth. It has been shown that the thin and soft subsoil has the bearing capacity to make the collector walk freely. The rate-dependent characteristics of deep-sea bottom soil mean that the deepsea bottom soil shows different mechanical properties when the grouser of the collector travel at different speeds. According to the physical and mechanical properties of deep-sea bottom soil, scholars at home and abroad continue to study and improve it, in which many researchers use the extended Drucker-Prager model to simulate the soil cutting process [33][34][35]. However, few researchers have considered the rate-dependent characteristics of deep-sea subsoil. In this paper, deep-sea subsoil's rate-dependent characteristics are studied, the rate-dependent parameters in the finite element software are calibrated with triaxial compression experiments, and the rate-dependent extended Drucker-Prager constitutive model is obtained. The soil parameters and rate-related parameters are shown in Tables 2 and 3.

Numerical Simulation Data
The force-increasing performance of a bionic grouser was studied by simulating a straight plate grouser and bionic grouser in the soil under the same conditions. The simulation results are shown in Table 4. Compared with straight plate grousers, bionic grousers have a greater traction force. The maximum traction force of a bionic grouser with a grouser tip distance L of 30 mm and curvature radius R of 500 mm is 17.52% higher than that of a straight plate grouser, which has a significant force-increasing effect.

Effect of L on the Tractive Force of a Grouser
As can be seen from Figure 7, the grouser traction force increases with the decrease in grouser tip distance L. This is because the grouser is composed of a straight plate part with a grouser tip distance L and a curved surface part with a curvature radius of R, and the grouser height is fixed, so when the grouser tip distance L decreases, the soil volume in the compression area increases. The length of the sliding surface of the soil also increases, and, when it is sheared, the sliding resistance of the soil and the traction provided to the grouser is greater.
in the compression area increases. The length of the sliding surface of the soil also increases, and, when it is sheared, the sliding resistance of the soil and the traction provided to the grouser is greater.  Figure 8 shows the relationship between the radius of curvature R of the curved part of the grouser and the maximum traction force Fmax under different grouser tip distance L. As can be seen from the figure, Fmax increases with the decrease of R. This is because when the grouser walks on the simulated soil, the bionic curved surface of the grouser shears the soil in the compression area. On the one hand, with the decrease of the radius of curvature and the distance from the top of the grouser, the volume of the soil in the slip zone increases, and the greater the tractive force provided by the soil to the grouser; on the other hand, because of the curvature part of the grouser in the process of walking, the soil in the slip zone produces vertical component force upward and forward horizontal component force, while the vertical component makes the soil in the slip zone denser, thus increasing the soil shear strength and increasing the grouser traction force.

Effect of R on the Tractive Force of the Grouser
In addition, it is obvious that the simulation results are slightly higher than the experimental results, but both are within the allowable fault tolerance range. The possible reason for this phenomenon is that the soil accumulation effect is considered in the numerical simulation. With the advance of the grouser, more and more soil accumulates, the resistance of the grouser becomes greater and greater, and the corresponding traction force is also greater.

Accuracy Verification of the Analytical Solution
By comparing the maximum traction force of the grouser with the same profile parameters in Figures 7 and 8, it is found that the grouser parameters of the maximum traction force are the same, in which the value of the analytical solution is 182.4 N, and that of the numerical simulation is 198.694 N. Furthermore, the maximum and minimum errors of the maximum traction force measured by the analytical solution and the finite element method are 8.2% and 0.1%, respectively, and the average error is 2.132%, which shows In addition, it is obvious that the simulation results are slightly higher than the experimental results, but both are within the allowable fault tolerance range. The possible reason for this phenomenon is that the soil accumulation effect is considered in the numerical simulation. With the advance of the grouser, more and more soil accumulates, the resistance of the grouser becomes greater and greater, and the corresponding traction force is also greater.

Accuracy Verification of the Analytical Solution
By comparing the maximum traction force of the grouser with the same profile parameters in Figures 7 and 8, it is found that the grouser parameters of the maximum traction force are the same, in which the value of the analytical solution is 182.4 N, and that of the numerical simulation is 198.694 N. Furthermore, the maximum and minimum errors of the maximum traction force measured by the analytical solution and the finite element method are 8.2% and 0.1%, respectively, and the average error is 2.132%, which shows that the derivation of the analytical solution has high calculation accuracy.

Force Enhancement Mechanism of the Bionic Grouser Cutting Deep-Sea Soil
As shown in Table 5, to study the force enhancement mechanism of the bionic grouser in cutting deep-sea bottom soil, the straight plate grouser and No. 1 bionic grouser (L = 30 mm, R = 500 mm) were selected for finite element cutting simulation, and the soil stress cloud diagrams at 1.2 s, 2 s and 4 s were obtained.

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27.0% compared with 166.3 N of straight grouser, and the force-increasing effect is significant.

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27.0% compared with 166.3 N of straight grouser, and the force-increasing effect is significant.

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27.0% compared with 166.3 N of straight grouser, and the force-increasing effect is significant.

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27.0% compared with 166.3 N of straight grouser, and the force-increasing ef-In the process of cutting deep-sea bottom soil with a grouser, the slip zone in front of the grouser is formed as shown by ABC in Figure 9, and the plastic deformation of this part of the soil is larger, which is shown in dark color in the cloud image. From the comparison of the cloud images of the straight plate grouser and bionic grouser in Table 5, it can be clearly seen that the plastic deformation region of the bionic grouser is larger. That is, the soil volume in the slip zone increases, and the length of the soil slip line also increases. When the soil is sheared, the sliding resistance is more excellent, and the traction force provided to the grouser is greater.

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27.0% compared with 166.3 N of straight grouser, and the force-increasing effect is significant. On the other hand, when the distance from the top of the grouser L decreases, the B point extends in the opposite direction of movement, and the soil directly below the AB curve will increase the upward vertical component when the grouser is shearing, resulting in a compaction effect between the grouser and the soil, resulting in greater friction between the soil and the grouser, thus increasing the tractive force of the grouser.

Force Enhancement Design of the Bionic Grouser
The previous data shows that the maximum traction force Fmax of the bionic grouser has apparent regularity with grouser tip distance L and curvature radius R. To explore this law, we used the Mesh function of MATLAB to make a three-dimensional image of the analytical solution. It can be seen from Figure 10 that the traction force increases with the decrease in grouser tip distance L and curvature radius R, and when L is 0 mm and R is 183 mm, the maximum traction force reaches the optimal value of 211.2 N, which increases by 27.0% compared with 166.3 N of straight grouser, and the force-increasing effect is significant.

Conclusions
According to the bionic structure of buffalo hoof, taking the distance L from the top of the bionic grouser to the point of curvature change and the radius of bionic curvature R as parameters, the influence equation of traction force on the bionic grouser profile parameters is established by optimizing Rankine earth pressure theory. The analytical solution of track parameters under the action of maximum traction force is obtained. The process of cutting deep-sea bottom soil with a grouser with different grouser profile parameters is studied by numerical simulation. The simulation results are close to the analytical solution, and the force-increasing mechanism of the bionic grouser is analyzed combined with the displacement cloud map. The correctness of the analytical solution is further verified. According to the three-dimensional variation curve of the grouser profile parameters and the maximum traction force, we found that the maximum traction force of the

Conclusions
According to the bionic structure of buffalo hoof, taking the distance L from the top of the bionic grouser to the point of curvature change and the radius of bionic curvature R as parameters, the influence equation of traction force on the bionic grouser profile parameters is established by optimizing Rankine earth pressure theory. The analytical solution of track parameters under the action of maximum traction force is obtained. The process of cutting deep-sea bottom soil with a grouser with different grouser profile parameters is studied by numerical simulation. The simulation results are close to the analytical solution, and the force-increasing mechanism of the bionic grouser is analyzed combined with the displacement cloud map. The correctness of the analytical solution is further verified. According to the three-dimensional variation curve of the grouser profile parameters and the maximum traction force, we found that the maximum traction force of the bionic grouser increases with the decrease of the distance between the grouser and the radius of curvature. The reasonable selection of grouser profile parameters (L = 0 mm, R = 183 mm) can effectively improve the tractive force of the grouser, provide theoretical support for the deep-sea collector's grouser design, and significantly improve the mining efficiency.
The future work will focus on the derivation of the new analytical solution, which will consider the asymptotic slip line and the soil arch effect and be verified by experiments.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.