Modeling and Investigation of the Velocity-Dependent Cutting Process with PDC Cutters Using the Discrete Element Method

,


Introduction
In the deep earth drilling process, polycrystalline diamond compact (PDC) bits with multiple fxed cutters (see also Figure 1(a)) are widely used to drill through rock formations.Te bit-rock interaction often causes severe drillstring vibrations that occur in the axial, lateral, and torsional directions.Among these unwanted vibrations, self-excited high-frequency torsional oscillations (HFTO) of the bottom hole assembly (BHA) often dominate the drill-string dynamics, leading to critical premature failures of downhole components [2].When this vibration type occurs, the BHA vibrates at a single eigenfrequency of the drill string between 50 and 400 Hz [3].Tis self-excited phenomenon is primarily caused by a velocity-dependent torque characteristic in the bit-rock interaction while drilling hard and dense formations [1].Te torque on bit (TOB) decreases as the rotary speed (i.e.RPM) increases, creating a negative slope at an operating point and thus destabilizing the drill string system due to energy infow, as illustrated in Figure 1(b).Consequently, self-excited torsional vibrations arise if the system damping is not sufcient to dissipate this energy input from the excitation.
To explain the genesis of the falling torque characteristic with increasing RPM, the frst argument is focused on whether the falling torque characteristic is primarily associated with rock properties.Richard et al. [4] claimed that the torque characteristic is not an intrinsic property of the bit-rock interface, but an indirect consequence of drill string vibrations (or stick-slip in this context).Te authors introduced a state-dependent but cutting-rate-independent bit-rock interaction model into a two-degree-of-freedom drill string model and managed to reproduce stick-slip in their simulations.Accordingly, the observed falling torque characteristic results from the reduction in depth of cut (DOC) due to the decreased intensity of axial vibrations at higher RPM.However, this theory holds true only if the modal frequency interval between the low-frequency torsional mode and the high-frequency axial mode is large enough.Tis condition is not satisfed for HFTO, where a strong coupling of torsional and axial vibrations with the same high frequency is observed [5].In this case, a more convincing explanation is that the falling characteristic is rooted in a constitutive rate dependency in the bit-rock interaction.
Since bit-rock interaction can be considered as the combination of all cutter-rock interactions on the bit, the velocity-weakening torque characteristic can be attributed to the rate-dependent rock cutting process at each cutter, i.e, the tangential cutting force decreases with increasing cutting velocity.Teoretically, this rate sensitivity could result from two potential aspects: either the reduction in tangential cutting force is a direct outcome of rock weakening at a higher velocity, or it could indirectly come from a mechanism that lowers the depth of cut (DOC) at a higher velocity because the tangential cutting force decreases with reduced DOC.In the frst case, the rate dependency is an uncoupled behavior restricted in the tangential direction, which is independent of the axial dynamics.For the latter case, the reduction of DOC over cutting velocity can be accomplished by the rate hardening of the axial cutting force, implicating a tangential-axial coupling efect.In this situation, the axial force increases with increasing cutting speed, thus creating an additional velocity-dependent lift force on the bit.Ten the drill bit is lifted if the weight on the bit (WOB) remains constant, which is an usual case in standard drilling operations.Several laboratory studies on the rock cutting process with a single PDC cutter have supported the latter hypothesis about the axial rate hardening efect, for example in [6,7].
If the velocity-dependent cutting process with PDC cutters is the root cause of the falling torque characteristic, fundamental strategies for reducing and mitigating the selfexcited HFTO can be deduced through modeling the rate hardening efect of cutting forces and investigating the infuence factors.Several studies have attempted to model the rate-dependent cutting process using diverse physical explanations.Based on the experimental observation of a shear rate strengthening efect in granular systems in [8], Pelfrene et al. [6] added a logarithmic formulation between the cutting speed and the normal stress on the wear fat into the cutting force model for a PDC cutter.With the same motivation, Matthews et al. [9] also included a velocity-dependent term in their drilling force model, achieving better agreement between simulated and laboratory test data.However, these simplifed modifcations only give a phenomenological description of the observed rate efect but are not usable to explore the infuencing factors.
Another explanation comes from the theory of poroelasticity.Te saturated rock formation at the pressurized bottom hole is a porous material.Terefore, the interaction between pore pressure dynamics and rock mechanical properties has an additional impact on the rock cutting process.With this consideration, Kolle [10] and Detournay and Atkinson [11] developed two similar analytical models to describe the dynamic confning pressure due to diferent pore pressure reductions on the shear dilatant plane at diferent cutting speeds.Tereafter, Amri et al. [12] and Chen et al. [13] independently extended the poroelastic modeling of the PDC cutting process by introducing solid- 2 Shock and Vibration fuid coupling to obtain stress and pore pressure felds in the rock.However, the rate efect caused by pore pressure difusion only applies to rocks with high permeability, which cannot explain the self-excited HFTO in hard and lowpermeable formations [5].Moreover, the application of these analytical models is limited to simple boundary conditions that are not able to describe the realistic cutterrock interface with various geometric confgurations.Instead of analytical approaches, numerical simulations are appropriate substitutions for studying the cutter-rock interaction during the rock cutting process.Compared with other numerical modeling methods, the discrete element method (DEM) is a favorable candidate due to its advantages in efectively accommodating the large deformations and fractures involved in rock cutting processes.Terefore, several studies have used DEM to simulate the rock cutting process with PDC cutters and to investigate various infuencing factors [14][15][16][17].Trough these DEM investigations, the understanding of cutter-rock interaction has been enhanced.However, to the authors' knowledge, modeling of the abovementioned rate dependency has rarely been covered in DEM cutting simulations.Helmons et al. [18] incorporated pore pressure modeling in DEM using smoothed particle techniques.However, diverse dynamic properties of rock materials can also accomplish the rate efects in the cutting process rather than being restricted to this physical mechanism that is negligible in low-permeability rocks.As reviewed by Zhang and Zhao [19], rate-dependent rock properties were commonly observed in various dynamic tests.Besides fuid diffusion, other aspects of rock properties should be explored in DEM modeling regarding rate dependencies, such as inertia efects, viscous behavior or dynamic fragmentation [19].
Tis paper aims at modeling and investigating rate dependency in the PDC cutting process based on DEM simulations developed in our previous studies [20,21].In the frst part, experimental cutting tests with a single PDC cutter are conducted and evaluated to obtain the characteristics of the rate dependency of cutting forces.Ten diferent modeling strategies of DEM for reproducing the observed rate dependency are discussed, and two potential contact models incorporating cutting-speed-dependent parameters are implemented.Interpreting the cutting simulation results from a continuum perspective explicitly reveals the relationship between rate-dependent cutting forces and the causal rock failure mechanisms.Finally, the infuences of typical operating and geometrical parameters on the rate dependency are investigated and compared with laboratory tests reported in the literature.

Single Cutter Test Rig.
To characterize the rate dependency in the rock cutting process, single cutter tests under hydrostatic confning pressure were conducted on a test rig [22].Te test rig mainly consists of a pressure vessel and a cutter shaft.A rock sample plate is placed in the vessel that is pressurized with fuid to simulate the confning bottom-hole pressure.Te cutter shaft is equipped with a load cell and its movement is controlled by a hydraulic servo system.
Te schematic of the cutting test inside the pressure vessel is shown in Figure 2. In these tests, the single cutter mounted on the load cell moves vertically down to penetrate the rock sample at a prescribed velocity.Meanwhile, the rock sample rotates at a constant RPM so that a stable DOC is reached after the frst revolution.Te vertical displacement of the cutter, the rotational speed, and the cutting forces in three directions (axial, tangential, and radial) are measured and recorded.Te cutting area can be geometrically calculated from the displacement signal.

Characterization of Rate Dependency.
In this study, all tests are performed with a 16 mm diameter cutter with a 0.016 inch chamfer at a 20 °back rake angle.Danby Marble is chosen as the rock specimen under 3000 psi confning pressure.To study the rate sensitivity of cutting forces, the tests are conducted at various rotational speeds and the same DOC.Te raw signals of the cutting forces are divided into a smoothed and a noise part frst using a low-pass flter.According to statistical and FFT analyses, the noise signal follows a normal distribution and demonstrates no obvious peak values in its frequency spectrum.Consequently, the noise can be neglected in further analysis.
Te smoothed cutting forces at diferent cutting speeds are shown in Figure 3. Since the cutting forces in the normal and tangential direction at the single cutter correspond to the weight on bit (WOB) and torque on bit (TOB), respectively, the test data evaluation will focus on these two directions.Te results indicate that cutting forces increase with the number of revolutions because the area of cut (AOC) keeps increasing during the cutting process with a single cutter and more energy for rock fragmentation is consumed.In addition, a distinct hardening efect with increasing RPM is found in the axial component, whereas the tangential force remains roughly independent on the cutting speed.Tis observed rate efect is consistent with the results on another rock type (Anstrude limestone) reported by Pelfrene et al. [6].
Te rise of the cutting forces in single cutter tests with unsteady AOC is not encountered while drilling with a full bit because the cutting area in front of an individual cutter overlaps with the grooves left by the previous cutter trajectories [24].To eliminate the infuence of increasing AOC, linear regression is performed based on the correlation analysis that shows that cutting forces strongly correlate with AOC among all test data.Te data from the frst revolution are ignored because the DOC does not reach a constant state yet.Following the defnitions in [25], the drilling strength (DS) and the specifc energy (SE) are obtained for each test, as shown in Figure 4. Tese two characteristic values represent the axial and tangential force components normalized to the AOC, respectively.Te cutter's aggressiveness, defned as the ratio of tangential and axial force, is given as well.
As can be seen, an obvious RPM-dependent hardening efect of DS systematically appears, while SE is nearly independent of RPM.In other words, a higher axial force is required to maintain the same DOC at a higher tangential cutting speed, while the tangential force does not change noticeably.Tis results in a falling characteristic of Shock and Vibration aggressiveness over the rotational speed.As discussed, this velocity-dependent characteristic could be the cause of selfexcited drill string oscillations.To enhance the understanding of self-excitation, this observed RPM efect will be incorporated into the DEM modeling.

DEM Simulation of the Rock Cutting Process.
As a starting point for further modifcations concerning the rate dependency, the DEM base model and the rock cutting simulation are introduced in this section frst.To model a granular material such as rock, the DEM explicitly considers individual particles and their interactions.Te motion of each particle is governed by the following equation: where € u i is the translational acceleration vector and _ ω i is the angular acceleration vector for particle i, m i is the mass, Θ i is the moment of inertia, F i and M i are external forces and Te particulate interactions determine the mechanical behavior of the granular system.Te corresponding contact forces are calculated by a constitutive contact model that describes the force-displacement law at each contact.For rock materials, a bonded-particle model proposed by Potyondy and Cundall [26] is often used, which decomposes the contact forces into diferent directions and portions as denoted by the following equations: where n stands for normal direction, t for tangential direction, p for particle portion, b for bond portion, and d for damping portion.
Tis nomenclature also suggests the additional rolling resistance torque M p due to the particle irregularity, as proposed by Tergeist et al. [27].A previous study has shown that this modifcation contributes to a better agreement with triaxial compression test data.
To calibrate the parameters of the DEM model for the Danby Marble used in this study, triaxial compression simulations under several confning pressures are conducted.Te simulated stress-strain curves are compared with the experimental data on Danby Marble, as shown in Figure 5.A good agreement with the test data is achieved over a large range of strains and confning pressures by selecting proper model parameters.Several relevant calibrated parameters are listed in Table 1.More details about the model parameters and calibrations can be found in [20,21].
To investigate the rock-cutter interaction, the calibrated DEM model is applied to the implemented two-dimensional cutting simulation, as illustrated in Figure 6.Te rock specimen, consisting of bonded particles is constrained by rigid frictional walls on the right and bottom sides.Along the left and upper surfaces, a hydrostatic pressure boundary is realized through a pressure application algorithm developed by Tergeist [20].Te compressive uniform load that simulates the fuid confnement is applied to the identifed particle surface chain, except for the dry contact between the cutting chip and the cutter (red particles in Figure 6).Te profle of the PDC cutter is represented by mesh segments with a constant coefcient of friction.Te cutter movement is separately defned in two orthogonal directions.A constant horizontal velocity is applied in the tangential direction, while either a constant DOC or a constant vertical force is applied in the axial direction, depending on diferent operating conditions.In the former case, the position and velocity of the cutter are completely prescribed and are not afected by its interaction with particles.For the latter case, a PID-controlled servo wall for the cutter is adopted to achieve force equilibrium in the axial direction.Shock and Vibration

Modifcation Strategies.
To reproduce the experimentally observed rate dependency, the contact model of the DEM is further modifed.In the single cutter tests, the rate dependency is characterized by the cutting speed (i.e.RPM), whereas for the DEM contact model, the rate indicator is only refected in the relative velocities between particles.
Tese two velocity measures are on diferent geometric scales.Choosing which velocity measure to incorporate the rate dependency is, thus, a primary option for diferent modifcation strategies.
Te other attribute to be selected is the model complexity.Te least complex modifcation is the introduction of characteristic curves instead of constants for the existing parameters in the DEM contact model.In comparison, adding more rheological elements to modify the forcedisplacement law is a more complex extension, but with more degrees of freedom.Te most complex modifcation can be accomplished by considering and modeling additional physical mechanisms such as fuid-solid coupling, hydraulic fracturing, and heat transfer.Tis is usually associated with very specifc application scenarios and is beyond the scope of this study.
Trough combinations of diferent choices for these two attributes, various modeling strategies regarding the rate dependency are generated, as schematized in Figure 7. Te feasibility and specifc implementation of each strategy are further discussed.

Model Modifcation with Microscopic Rate Parameters.
Considering the complexity and intuitiveness of the implementation, the most convenient approach to extending the contact model is to introduce parameters or rheological elements related to the contact velocity between particles, namely Strategies A and B marked in Figure 7.However, the relative velocity feld during the cutting process has a nonuniform spatial distribution that depends on the particle fow driven by the cutter movement and geometry [15].Terefore, the relationship between the local contact velocities and the macroscopic cutting speed must frst be clarifed to check the feasibility of modifcations using the microscopic rate parameter.Tis request is compatible with the advantage of DEM simulations in that the positions and velocities of the discrete particles and the contact forces of the particle pairs at each numerical step are available for an in-depth feasibility analysis.
For this purpose, the normal contact forces and the relative velocities at diferent cutting speeds are compared, as shown in Figure 8. Te results are time-averaged over the whole cutting distance using the previously developed averaging approach [21] to achieve a better statistical representation.As can be seen, the higher normal forces concentrate primarily under the chamfer surface.Our previous study [21] also showed a consistent phenomenon in the volumetric stress by which a crush zone was identifed.Correspondingly, this concentration of stresses or contact  Shock and Vibration forces has a decisive infuence on the axial cutting force.Terefore, modifcations for the observed axial hardening efect should also take efect in this region.In contrast, the change of the relative velocity feld at diferent cutting speeds appears in other regions where the shear zone is formed in front of the cutter face.Tis contradiction implies that modifcations based on interparticle velocities are not suffcient to reproduce rate-dependent compression behavior.Instead, the macroscopic cutting speed should be involved as the rate indicator to modify the contact model.

Rate-Dependent Normal Bond Strength.
For modifcations with macroscopic rate parameters, the simplest way is to make the existing parameters of the bondedparticle model dependent on the cutting speed (Strategy C marked in Figure 7).For this purpose, a thorough sensitivity analysis of the contact parameters is performed to investigate their infuences on the cutting force responses.Te normal bond strength is then found to be the only

Shock and Vibration
parameter that leads to a dominant increase in the axial cutting force when it is increased.Tis result implies an appropriate modifcation by incorporating a normal bond strength related to the cutting speed.To verify the feasibility of this modifcation, DEM cutting simulations are conducted at a constant DOC and three diferent cutting speeds (i.e.RPM), while the normal bond strength is varied in each simulation run.To keep it simple, the relationship between the normal bond strength and the cutting speed is set to be linear.Te other calibrated model parameters stay the same.Figure 9 shows the simulated time-averaged force responses and the resultant cutter aggressiveness at various cutting speeds ranging from 30 rev/min to 90 rev/min.Note that a virtual depth of 1 m is assumed in the third dimension for a physical interpretation.Te force is thus denoted as a quantity per unit length.Te shape of the cutting groove is assumed to have a negligible efect on the cutting forces.It can be seen that the axial cutting force demonstrates a more pronounced tendency to increase with increasing cutting speed, which generates a falling characteristic in aggressiveness over RPM.Tis consistency with the rate dependency observed in single cutter tests validates the good functionality and reproducibility of applying this modifcation strategy.Further parameterization reveals that the rate dependency provided by the normal bond strength with cutting speed dependency is limited because the bond strength only takes efect if the bond is not broken.However, the neighboring bonds of the particles near the cutter edge are usually all broken in the cutting simulations.Tis phenomenon was also observed in the PDC cutting experiments, where the intact rock next to the single cutter surface was pulverized into loose grains [28].

Rate-Dependent Grain Elastoplasticity.
To extend the rate dependency, more aspects should be considered by introducing additional degrees of freedom into the contact model.Preceding DEM cutting simulations [14] have shown that the energy consumed by the plastic deformation of the crushed rock under pressure is much higher than the energy dissipated on breaking the bonds between particles.In the base model of the bonded-particle method, rock plasticity is represented through bond breakage, the interparticle frictional contact and rolling resistance.Te particle size is predetermined by a uniform distribution and is thereafter assumed to be unchangeable because only elastic and viscous elements are applied in the normal contact direction.In contrast, rock grain crushing was observed by Richard et al. [29] in rock cutting tests.Tis plasticity of crushed grains is also supported by the quantitative investigation on the cuttings size [30]; the distribution of the cuttings size has a broad range from chunklike particulates to fne powder.Mendoza et al. [31] incorporated the particle-crushing effects into DEM by explicitly implementing crushable particles.However, simulating the event of a particle splitting  Shock and Vibration into smaller particles is an unstable and computationally intensive task [32].
As a substitute, an abstract formulation for elastoplasticity is deployed in the normal contact model of the grain portion.An additional module consisting of a slider element for plasticity and a second spring for hardening is connected in series with the original spring for elasticity, as illustrated in Figure 10.Te slider controls the onset of the plastic stage through a threshold value δ p for the overlap.Tis elastoplastic contact model is mathematically expressed by the following equations: where k e is the elastic stifness, k h is the hardening stifness, δ is the overlap, δ max is the maximum overlap, and δ 0 is the residual overlap at which the force decreases to zero during the unloading phase.Note that the maximum overlap δ max is a history variable that is continuously updated since the beginning of contact.Te residual overlap δ 0 is calculated by the following equation: Te corresponding piecewise force-displacement relationship is shown in Figure 11.Te energy loss due to the unrecoverable plastic deformation is dissipated by the hysteresis damping that occurs in the contact cycles throughout the loading and unloading phases.Tis energy dissipation is decoupled from the local contact velocity and can be better controlled to achieve various behaviors of grain elastoplasticity.
Tis modifed normal contact model is then integrated into the existing DEM model.Te sensitivities of the cutting force respond to both additional contact parameters (i.e.δ p and k h ) are then investigated.It is found that the axial cutting force increases faster than the tangential force as the yield threshold δ p increases, while the hardening stifness has efects in both directions in equal measure.Terefore, the observed rate efect of axial force can be reasonably reproduced by the yield threshold, depending on the cutting speed (Strategy D in Figure 7).Te threshold overlap δ p is chosen to be 2e−6 m as a reference value and increases linearly with the cutting speed.Other additional parameters for the grain elastoplasticity are selected through systematic parametric sensitivity studies to ensure their efectiveness in the considered range of cutting speed.Accordingly, a more signifcant increase of axial force with increasing cutting speed is obtained from the cutting simulations, which eventually leads to the velocity-weakening characteristic of the cutter aggressiveness, as shown in Figure 12.Tis modifcation with grain elastoplasticity takes efect mainly after bond breakages and thus complements the rock rate dependency in its full degradation during the PDC cutting process.

Causal Rock Failure Mechanisms.
In contrast to the empirical model [6], the two modifcations produce more than merely phenomenological representations of the observed rate dependency in the cutting process.Instead, they connect the macroscopic cutting speed with the microscopic contact parameters (i.e. the normal bond strength and the yield threshold). Terefore, the macroscale rate dependency is generated from interactions of microscale particles during the rock failure process in the cutting simulations.To understand the root cause, the corresponding simulation results are postprocessed from the continuum perspective using the approach in our previous study [21].Te volumetric stresses are derived from the contact forces and the locations of particles and then averaged over the whole cutting distance.
Figure 13 depicts the obtained volumetric stress distributions at two diferent cutting speeds using the two modifcations, respectively.It can be seen that the stress distribution is changed if considering the grain plasticity; the Shock and Vibration stress concentration in the crush zone below the chamfer is reduced, and more support for the cutter is created by the compressed cuttings in front of the cut face.Nevertheless, the rate dependency simulated by both modifed rate-dependent models still results from the increase in stress concentration at higher cutting speeds.Tis implies a strong correlation between the rate-hardening efect of the axial cutting force and the compression behavior of the crushed rock below the chamfer.In this way, the causal rock failure mechanisms for the simulated rate dependency are revealed.

Falling Characteristics over RPM.
Te falling characteristic over RPM is refected by the cutter aggressiveness in the above simulations with a constant DOC using the two modifed models incorporating the rate dependency.In the   drilling process with a full drill bit, the falling characteristic over RPM appears more often in the torque on the bit when WOB is held constant, as observed by some measurements both in the laboratory and in the feld [6,33,34].
To verify the falling characteristic under similar boundary conditions in the cutting process with a single cutter, DEM simulations with a constant axial force control are conducted at various cutting speeds.A PID controller is used to achieve the target force in the axial direction while the tangential cutting speed is prescribed with a constant value.Te steady-state cutting forces in both directions are averaged.With the rate-independent material model, no perceptible change occurs in the tangential direction, as shown in Figure 14.
By contrast, a falling characteristic of the averaged tangential force with respect to the RPM is obtained using both cutting-speed dependent models, as shown in Figure 15.Tis falling characteristic is caused by the incorporated rate dependency.At higher tangential speeds, the force response in the axial direction becomes larger, pushing the cutter upwards under constant force control to regain the force equilibrium.Ten the cutting depth is reduced, resulting in a smaller tangential force reaction.Terefore, the reduction in bit torque with increasing rotary speed can be facilitated by the combined axial rate-hardening efects on every single cutter.

Parameter Sensitivities.
Te full-scale test drilling [33] has exhibited that the falling torque characteristic is sensitive to WOB and the bit wear state (sharp or dull).To check the reproducibility of these parameter sensitivities in DEM simulations, the infuences of typical operating and cutter geometrical parameters on the falling characteristic are further investigated using the cutting scenario with the axial force control.For the operating parameter, the target axial force that corresponds to WOB is varied.In parallel, two diferent cutter meshes, an unchamfered sharp cutter and a dull cutter, are adopted for investigating the sensitivity to the cutter wear state.Teir infuences on the falling characteristics are evaluated by the cutter's aggressiveness.Te resultant falling characteristics using both rate-dependent models with various axial forces and cutter states are shown in Figures 16 and 17, respectively.It can be seen from Figure 16 that the simulated falling characteristic over RPM becomes more pronounced with a larger target value of the axial force because the axial rate hardening efect is intensifed at a larger axial force.Te intensifcation of the rate efect also occurs in Figure 17 with a dull cutter due to the enlarged contact area in the axial direction compared to that of a sharp cutter, resulting in a more signifcant velocityweakening characteristic on a dull cutter.Tese simulated trends are phenomenologically in good agreement with experimental measurements on full bits [33] and thereby confrm the validity of the two proposed modifcations.Te simulation results also implicate a fundamental possibility to reduce the undesired rate dependency and consequently self-excited drill string vibrations by changing operating conditions or keeping the cutter's sharpness.

Conclusion
In this paper, the velocity-dependent cutting process with a PDC cutter is modeled and investigated using the discrete element method.Te rock material is modeled by an aggregation of bonded particles, and the rock-cutter interaction in deep drilling environments is investigated through cutting simulations under algorithmically realized hydrostatic confnement.According to measurements in the pressurized single cutter tests, the axial cutting force is found to have a nonnegligible hardening efect with increasing cutting speed.Tis phenomenon occurring on the cutter scale is supposed to be responsible for the generally observed unfavorable falling torque characteristic in the drilling process with a PDC drill bit.
To enhance the understanding of this rate dependency, the particle model developed in our previous studies is further modifed.Potential modifcation strategies considering diferent model complexity and rate indicators on diferent scales are compared and discussed.Two specifc modifcations associated with the macroscopic cutting speed are implemented and verifed to reproduce the experimentally observed rate dependency: one modifes the normal bond strength in the existing bonded-particle model; the other extends the contact model with an abstract formulation of grain elastoplasticity and then adds the cutting speed dependency to the yield threshold.Both modifed models result in a dominant rate hardening efect of the axial cutting force in cutting simulations at a prescribed depth of cut and various cutting speeds.
In addition to qualitatively reproducing the observed rate dependency, further analysis and interpretation of the Shock and Vibration simulation results reveal the relationship between the rate-sensitive cutting forces and the causal rock failure mechanisms.Te cutting simulation is also extended by implementing a constant axial force control.Te averaged tangential force at a steady state is found to decrease with increasing cutting speed, which is consistent with the general falling torque characteristic over RPM on a full bit.Moreover, the infuences of axial force and cutter wear state on the falling characteristic are simulatively investigated in this cutting scenario.Te corresponding parameter studies indicate that the falling characteristic over RPM intensifes at a higher axial force or with a dull cutter, which is in good qualitative accordance with laboratory full-scale bit tests in the literature.Ten the infuences of studied parameters on the stability of the drill string system can be evaluated because a more pronounced falling torque characteristic leads to more instable behavior of the drill string and faster self-excitation of torsional vibrations.
Tis study provides an in-depth discussion on efectively modeling the velocity-dependent cutting process with PDC cutters using the discrete element method.Te corresponding simulation results reveal the underlying mechanisms on the single cutter scale that can give a consistent explanation for the experimental observations about the rate sensitivity of the cutting forces on a single cutter and the falling torque characteristic on a full bit.Te parameter studies enlighten the possibility of mitigating self-excited drill string vibrations by changing operating states or cutter design.Te fndings could contribute to a better cutter design that generates less torsional vibrations and a more stable drilling state, thus allowing drilling to achieve a higher rate of penetration.Future work will focus on investigating the infuences of more complex cutter designs (e.g.shaped cutters) on the rate dependency.Te rate-dependent cutting simulations will also be coupled with drill string models to study the efect of cutter-rock interaction on drill string dynamics.

Table 1 :Figure 6 :
Figure 6: Boundary conditions in the DEM cutting simulation.

Figure 7 :
Figure 7: Modifcation strategies with respect to the rate dependency.

Figure 8 :
Figure 8: Averaged normal contact force (a) and averaged relative velocity (b) with various cutting speeds (30, 60, and 90 RPM from left to right).

Figure 9 :
Figure 9: Averaged forces and aggressiveness at various cutting speeds using the model with rate-dependent normal bond strength.

Figure 10 :
Figure 10: Schematic diagram of the elastoplastic normal contact model for particle portion.

Figure 12 :Figure 11 :Figure 14 :
Figure 12: Averaged forces and aggressiveness at various cutting speeds using the model with rate-dependent grain elastoplasticity.

Figure 13 :
Figure 13: Averaged volumetric stress with 30 RPM (top) and 90 RPM (down) using models with rate-dependent normal bond strength (left) and grain elastoplasticity (right).

Figure 15 : 1 FFigure 16 :
Figure 15: Averaged forces and aggressiveness at various cutting speeds under constant axial force control using models with ratedependent (a) normal bond strength and (b) grain elastoplasticity.

Figure 17 :
Figure 17: Efect of cutter wear state using the models with cutting speed dependency on (a) normal bond strength and (b) grain elastoplasticity.