Analysis and Testing of Load Characteristics for Rotary-Percussive Drilling of Lunar Rock Simulant with a Lunar Regolith Coring Bit

Based on an optimized lunar regolith coring bit (LRCB) configuration, the load characteristics of rotary-percussive drilling of lunar rock simulant in a laboratory environment are analyzed to determine the effects of the drilling parameters (the rotational velocity, the penetration rate, and the percussion frequency) on the drilling load. The process of rotary drilling into lunar rock using an LRCB is modeled as an interaction between an elemental blade and the rock.The rock’s fracture mechanism during different stages of the percussive mechanism is analyzed to create a load forecasting model for the cutting and percussive fracturing of rock using an elemental blade. Finally, a model of the load on the LRCB is obtained from the analytic equation for the bit’s cutting blade distribution; experimental verification of the rotary-impact load characteristics for lunar rock simulant with different parameters is performed. The results show that the penetrations per revolution (PPR) are the primary parameter influencing the drilling load. When the PPR are fixed, increasing the percussion frequency reduces the drilling load on the rock. Additionally, the variation pattern of the drilling load of the bit is in agreement with that predicted by the theoretical model. This provides a research basis for subsequent optimization of the drilling procedure and online recognition of the drilling process.


Introduction
During the sampling and exploration of extraterrestrial celestial bodies, humans normally employ rotary-percussive drilling to improve drilling efficiency to obtain pristine samples from a celestial body, reveal the origin and evolution of the celestial body, explore the possibility of extraterrestrial life, and effectively explore and properly use space resources [1,2].The former Soviet Union's "Luna" series of unmanned lunar sampling missions, the United States "Apollo" series of sampling projects, which included a manned moon landing, and the "Curiosity" Mars exploration mission all employed percussion-assisted rotary drilling to ensure reliable sampling during exploration [3][4][5][6].In the near future, China will complete an unmanned lunar surface sampling and return mission, that is, the "Chang' e project," which will involve collecting soil samples from the lunar surface and 2 m below it using drill sampling [7,8].When the drilling process is blocked by a large chunk of lunar rock, percussion will be used to improve the bit's ability to fracture rock.Therefore, when the drill's configuration is determined, the ability to adjust the drilling parameters to achieve the desired drilling performance is necessary to guarantee reliable lunar surface sampling.
When the former Soviet Union conducted the Luna-24 unmanned autonomous sampling mission on the lunar surface, the experience of the Luna-16/20 drilling/sampling mission was drawn upon, and a percussion function was included in the drilling machine.During the actual lunar surface sampling, the drill's alarm was triggered due to an excessive drilling load.After percussion was enabled, the drilling load decreased to a certain extent, and finally, a lunar soil sample was collected successfully and sent to Earth [9,10].Recently, many researchers have conducted extensive research on rock fracture via drilling tool percussion.Among them, Jalali and Zare Naghadehi conducted an experimental study of excavation and percussion of rock using a disc cutter tooth in a tunnel excavating machine.The result showed that the sample excavation efficiency of the disc cutter tooth improved significantly under the combined effect of cutting and percussion.Additionally, the percussion load effectively improved the penetration rate of the simulated tunnel excavating machine [11].Hashiba et al. conducted an experimental study of a single percussion rock penetration process via button bits and modified the obtained loadpenetration depth curve using numerical simulation [12].Based on the fracture criteria for viscoplastic rock, Saksala et al. conducted a simulation analysis and an experimental study of the percussive drilling of granite via a specially designed triple-button bit.The results showed that the rock fractured when the percussion velocity was ascending.Additionally, lateral compressive fracture occurred between adjacent buttons [13][14][15].Furthermore, based on the Rayleigh-Love model and a one-dimensional viscoelastic model, Tian et al. derived an equation for the axial vibration of a drill string during oil exploration and verified it experimentally.The results showed that the length of the drill string is inversely proportional to the damping at the bottom of the bore.When the string was shorter, the effect of horizontal inertia on the dynamic rigidity of the drill string was more significant [16].
Besides, acoustic emission (AE) techniques have always been adopted in identifications of microseismic source (MS) locations.Dong and Li have proposed an innovative MS/AE source location method without the need for premeasured PS wave velocities, and this method greatly weakens the influence of the big error derived from the change of the wave velocities and propagation traces [17].Then, he applies three nonlinear methodologies (such as random forests, support vector machines, and naive Bayes classifier) to discriminate between 20 earthquakes and 27 nuclear explosions.And based on the leave-one-out cross-validation, ROC curve, and test accuracy, random forests method clearly shows the best predictive power [18].Recently, Dong et al. investigate the typical characters of seismic parameters between events and blasts using the datasets of three mines in Canada and Australia, and 7 parameters are extracted as discriminant indicators to establish discriminant models using the logistic regression method.The classified accuracy of training samples, test samples, and cross-validated results have demonstrated that the discriminant models with good discriminant performance are effective and efficient methods to identify the blasts from events in real-time [19].These methodologies provide an effective way to quantify the process of rock fragmentation which are not easy to observe in rock impacting or cutting process.
In 2014, the Research Center of Aerospace Mechanism and Control (RCAMC) at Harbin University of Technology successfully designed a lunar regolith coring bit (LRCB) that exhibited excellent drilling performance in tasks such as continuous chip removal without an auxiliary medium and rock blockage breakthrough with a limited ability to drill.Compared with other geological drill bits and planetary sampling bits, it has advantages, such as a low load per unit area, for fracturing rocks [20,21].In this paper, based on an optimized LRCB configuration and the Mohr-Coulomb criterion, a model of cutting and percussive fracture for lunar rock simulant is created and expanded to the rotarypercussive drilling process to analyze the load characteristics of rotary-percussive drilling of lunar rock simulant with a coring bit.The effect of the impact load on the rock's fracture load for different drilling parameters is determined; this provides a basis for optimizing the drilling parameters for lunar surface drilling and sampling.
The structure of this paper is as follows: Section 2 mainly introduces the background of previous research, that is, the LRCB configuration and the model for the load during rotary drilling of lunar rock simulant; Section 3 describes how the interaction between the LRCB and lunar rock simulant during rotary percussion is modeled; Section 4 provides experimental verification of the rotary-percussive drilling load for lunar rock simulant, and Section 5 is the conclusion.

LRCB Configuration and
Drilling Load Analysis the coring rate of the lunar soil sample and ensuring the ability to collect original information; and (iii) being able to break lunar rock with less than grade-VI drillability without exceeding the drilling load limit [22].Therefore, the RCAMC team developed the LRCB shown in Figure 1, which employs inner and outer rows of sharp cutting blades (SCBs) (i.e., dual primary cutting blades) to increase the stability of the drilling system in the initial stage of lunar rock drilling.
Laboratorial tests showed that this bit excels at coring lunar soil and removing chips and has excellent lunar rock drilling load characteristics.The detailed configuration parameters of the LRCB are listed in Table 1.

Analysis of the Drilling
Load of the LRCB.The LRCB's configuration parameters are abstracted as an analytic geometrical model, which is shown in Figure 2. The relationships of the torque on bit (TOB) and the weight on bit (WOB) to the cutting load of the elemental blade are shown in where  _ is coordinate of the midpoint of the th elemental blade at the th row of SCBs on the I axis;  _ is coordinate of the midpoint of the th elemental blade at the th row of SCBs on the K axis; Δ _ is cutting force of the th elemental blade at the th row of SCBs along the l axis; and Δ _ is cutting pressure of the th elemental blade at the th row of SCBs along the n axis.
In (1), Δ _ and Δ _ are the cutting loads of elemental blade in the absolute coordinate system.The cutting load of where   is resultant force on the rock fragment due to the front blade surface;   is resultant force on the overcut area due to the rear blade surface;   is friction angle between the cutting blade and the dense core;   is friction angle between the cutting blade and the parent rock; g is normal vector of the front blade surface; U is chip removal velocity vector; g  is unit vector normal to the rear blade surface; and   is unit vector normal to the rear surface of the inward cutting blade.In (2), the chip removal velocity vector, U, satisfies where  is removal-cutting ratio (the ratio between the removing vector and the cutting vector);   is chip removal angle (the angle between the direction of chip removal at the front blade surface and the direction normal to the cutting blade); a is unit vector normal to the front surface of the cutting blade; and b is unit vector parallel to the cutting blade.
Based on minimum energy dissipation theory, when the system's other control parameters are constants, there must be a group of state parameters {,   } that minimize the power consumed by the front surface of the cutting blade.The obtained state parameters {,   } are substituted into (2) to obtain the cutting load on the elemental blade.For a detailed derivation of the LRCB's drilling load, please refer to [20].
Because the process of sampling the lunar surface is based on controlling the drilling velocity, the LRCB's drilling parameters are generalized as the rotational velocity,   , and the penetration rate, V Pen .The penetrations per revolution (PPR) are defined as the ratio of these two control parameters, and its actual physical significance is the progress of bit as it penetrates vertically during each rotation.When the drill bit progresses at a constant PPR, each cutting blade performs a downward spiral.After the cutting blade cuts around the bit's center of rotation for a cycle, the corresponding penetration displacement is the product of the number of cutting blades and a single-blade's cutting depth, as shown in where ℎ Pen is the real-time cutting depth of a single cutting blade during the drilling process and   is the number of cutting blades deployed around the bit.

Mechanical Model for Rotary-Percussive Drilling
3.1.Analysis of the Rotary-Percussive Drilling Process.In the process of lunar surface sampling, applying a certain amount of impact load to the drill effectively improves its efficiency at fracturing lunar rock simulant.According to the design of the deep drilling/sampling mechanism for use on the lunar surface, the principle underlying the operation of the percussive mechanism is as follows: when the motor rotates the cylindrical cam, it also drives the cyclic motion of a set of spring-mass elements to apply an impact load to the drill, as shown in Figure 3 [24,25].
In rotary-percussive drilling process, the cutting blades interact directly with the rock.When rotary-percussive drilling of lunar rock with a coring bit is expanded along the bit's circumferential direction, it is equivalent to rock cutting and percussive fracturing using an elemental blade along a straight line, as shown in Figure 4.
Impact-cutting and fracture of lunar rock simulant using a single blade is divided into the following stages (as shown in Figure 4): (i) the impact-cutting coupling stage: the impact load is suddenly applied to the cutting blade during forward cutting by a hammer; at the same time, the cutting blade quickly changes its original track, penetrates the rock along a direction that combines the percussion velocity and the cutting velocity, and creates a tiny impact fracture pit; (ii) the force-control cutting stage: because the effective impact time is extremely short (far less than the time required for the cam to return), the cutting blade progresses in response to a percussive spring pretightening force to cut and fracture the  rock after the impact load disappears; at the same time, the cutting depth of the rock is determined by the pretightening force of the percussive spring; and (iii) the position-control cutting stage: when the hammer and the roller start the next cycle of motion, the percussive spring pretightening force on the cutting blade disappears; at the same time, the cutting blade enters the position-control cutting stage; that is, it completely follows the predefined cutting depth to fracture the lunar rock simulant.Now, the impact load cycle is complete, and the cutting blade waits for the next impact.

Mechanical Model of the Collision between the Hammer and the Drilling Tool.
Assume the hammer's collision with the drill is an interaction between a single mass block and a semi-infinite polished rod via a spring damping element without the additional effect of reflected waves on the percussive stress, as shown in Figure 5 [26].
According to Newton's second law and the equation for the transmission of a one-dimensional stress wave transmission, the collision between the hammer and the drill should satisfy the following equations and boundary conditions: where   is hammer mass;  is elastic modulus of the drill material;  is contact area between the hammer and the drill during the collision;   is contact rigidity of the hammer and the drill during the collision; and   is equivalent damping coefficient of the hammer and the drill during the collision.The Laplace transform of ( 5) is taken, variables are substituted, and the Laplace transform is inverted to obtain where and   ,   , and  0 satisfy The impact force on the contact surface is

Model of the Impact-Cutting Coupling
Stage.In a complete impact load cycle, because the effective impact time is extremely short (several milliseconds), the horizontal motion of the cutting blade during this cycle is ignored; that is, the fracture of the lunar rock simulant due to the impact-cutting coupling is discretized as a percussive fracture process and a cutting fracture process.

Mechanical Model of the Percussive Penetration of the Elemental Blade into the Lunar Rock Simulant.
As with a single impact, the interaction between the cutting blade and the rock is extremely short.The load,  Per , on the rock due to the cutting blade and the penetration depth, ℎ Per , satisfy the following relation (the variation of the rock's contact rigidity during the loading/offloading process is ignored): where  represents the contact rigidity of the rock.An analysis of the stress at the instant of percussive fracture of the lunar rock simulant via the cutting blade is shown in Figure 6.Based on the fluctuation dynamics, the percussion is transmitted as a fluctuation via the cutting blade to the rock.On the contact surface between the cutting blade and the rock, the wave is transmitted and reflected.Based on force equilibrium and fluctuation dynamics, the relationship between the incident/reflected stress and the velocity is where V Per is velocity of the cutting blade's percussive penetration into the rock;   () is cutting blade's incident impact force on the rock;  Pen is cutting blade's pressure on the rock;  is wave damping of the drill.
After substituting variables and rearranging (11), a differential equation or the percussive fracture of lunar rock simulant by a cutting blade is as follows: When the initial condition  Per | =0 = 0 is substituted into the incident percussion equation ( 9), the impact force that the cutting blade exerts to fracture the lunar rock simulant is where  * , the impact coefficient of the cutting blade on the rock, is equal to /.Because the percussive fracture region of the rock is much smaller than the width of the cutting blade, the impact force of the cutting blade on the rock calculated in the previous section is treated as a linear load.For a linear load  Per /  , the internal stress state of the rock (a semi-infinite plane, as shown in Figure 6) is where   is polar radius of pole  in a polar coordinate system;   is polar angle of pole  in a polar coordinate system;   is normal stress on a microelement of the rock along the polar radius;   is normal stress on a microelement of the rock along the polar angle; and   is shear stress on a microelement of the rock along the polar radius/angle.As shown in Figure 6, the shear stress   = 0; therefore,   and   are main stresses on the microelement of the lunar rock simulant.Under a linear load, shear cracks in rock normally spread along tracks that maintain constant angles with respect to the maximum principle stress.For a linear load, these tracks constitute two logarithmic spirals, as shown in where  0 is the distance between the intersection point of the two logarithmic spirals and pole  and   is the angle between the tangent at any point along the shear crack's track and the maximum principle stress at that point.
Next, based on equation of state for the stress on any microelement in an object under external load and ( 14) and ( 15), the shear stress,   , and the normal stress,   , at any point on the logarithmic spiral are Additionally, rock fracture satisfies the Mohr-Coulomb yield theorem; that is,   =   ⋅ tan  + .Equation ( 16) is substituted, and ( 17) is obtained after rearrangement.
Equation (17) suggests that when the rock surface is under a certain load, both sides of the load develop two shear cracks, as shown in Figure 6.These cracks satisfy relationship given by (15).The distance between the intersection of the two cracks and pole  is  0 .This area is defined as the compressive dense area formed after the impact load is applied to rock.After  0 is obtained, a complete profile of the impact fracture pit on the lunar rock simulant is obtained via (15).

Analysis of the Cut Fracture Characteristics of Lunar
Rock Simulant during Impact.After an impact, two impact fracture pits are formed before and after the load is applied to the lunar rock simulant.While the cutting blade continues cutting forward, the lunar rock simulant continues to fracture along the shear fracture previously formed by the impact load.To simplify the stress analysis, the shear fracture surface of the rock during cutting is treated as a plane whose starting and end points are where the profile of the impact fracture pit (a logarithmic spiral) intersects the cut and uncut surfaces of the lunar rock simulant.When the impact and cutting are coupled, three regions are formed in front of the cutting blade (as shown in Figure 7); that is, OAGF is the compact fracture region formed by the impact; ABCD is the percussive fracture fragment formed by the impact; and BEC is the shear fracture region formed by cutting.The equilibrium equations are as follows: where   is resultant force on the cutting blade;   is horizontal component of the resultant reaction force exerted by the parent rock on a rock fragment in the fracture region; and   is vertical component of the resultant reaction force exerted by the parent rock on a rock fragment in the fracture region.
In (19),   satisfies following relationship: where  is stress distribution coefficient; ℎ  is actual cut depth;   is cutting blade width;  sh is shear fracture angle;  s cohesion of the rock;  is inner friction angle of the rock;   is angle between the percussive fracture fragment and the cutting direction; and   is friction angle between the percussive fracture fragment and the parent rock.
For the other parameters and variables in Figure 7 and (19), please refer to the stress analysis for the cutting and fracturing of lunar rock simulant using an elemental blade [17].Before the resultant force on the cutting blade,   , is calculated, the resultant force exerted by the parent rock on a rock fragment in fracture region,   , should be calculated.Assume uniformly distributed force, , in logarithmic spiral AB satisfies the following condition: where  0 is the stress constant.The polar function in (15) is converted to rectangular coordinates as follows:

Shock and Vibration
The arc length,   , between any two polar angles   and  +1 in the logarithmic spiral in the impact fracture pit is The resultant force along any arc is where   is the polar angle corresponding to the midpoint of the arc between adjacent polar angles   and  +1 ;   ≈ (  +  +1 )/2.The horizontal and vertical components of the resultant reaction force exerted by the parent rock on a rock fragment in fracture region are where   is the angle between the normal vector/reaction force at the midpoint of arc microelement   and the horizontal direction; Equations ( 23) and ( 25) are substituted into (24) to obtain the horizontal and vertical components of the resultant reaction force exerted by the parent rock on a rock fragment in the fracture region, Equation ( 26) is substituted into (18) to obtain where When impact and cutting are coupled, the cutting force and cutting pressure required by a single cutting blade to fracture lunar rock simulant are as follows: (29)

Fracture Load Analysis for Lunar Rock Simulant in the
Force-and Position-Control Cutting Stages.The depth of a cut in lunar rock simulant during the force-control cutting stage is determined by the percussive spring pretightening force.If the percussive spring pretightening force exceeds the cutting pressure required to cut and fracture rock, the actual cut depth exceeds the initial predefined cut depth.Otherwise, the cutting blade fractures the lunar rock simulant according to predefined cut depth.Moreover, the cutting load in the force-control cutting stage is affected by both the percussive spring pretightening force and the impact load.In the impactcutting coupling stage, the lunar rock simulant develops an impact fracture pit during the impact.While the cutting blade continues cutting forward, the rock fractures along the trail of the impact fracture pit.Therefore, part of the load in the forcecontrol cutting stage is due to the impact pit in the fractured lunar rock.Additionally, when the cutting velocity is fixed, horizontal cut displacement in the force-control cutting stage is primarily affected by the percussion frequency.A higher percussion frequency leads to a smaller cut displacement; that is, a higher percussion frequency leads to a smaller displacement in the force-control cutting stage.It takes more time for the cutting blade to cross an impact fracture pit, and the load required to fracture lunar rock simulant increases when there is an impact pit.The average load required to fracture lunar rock simulant in the force-control cutting stage is where  Cut_cra is cutting force required to complete the first cutting cycle after the impact fracture pit is formed;  Pen_cra is cutting pressure required to complete first cutting cycle after the impact fracture pit is formed;  cra is weight coefficient of the required cutting load to complete the first cutting cycle after the impact fracture pit is formed;  Cut_pre is cutting force required to cut the rock with a percussive spring pretightening force;  Pen_pre is cutting pressure required to cut the rock with a percussive spring pretightening force;  pre is weight coefficient of the cutting load for the pretightening force.In (30), the weight coefficient satisfies the following condition: where  cra is time required for the cutting blade to cross the impact fracture pit;  for is time required for the cutting blade to complete the force-control cutting stage;  cra is horizontal displacement of the impact fracture pit;   is percussion frequency; V Cut is horizontal velocity of the cutting blade;  1 is idle movement angle of the cam in percussive mode; and  cb is number of bulges in the cam in percussive mode.The cutting load characteristics of lunar rock simulant under force-control cutting conditions are essentially the same as those with a predefined cut depth.Therefore, the calculation of the cutting load of lunar rock simulant under force-control cutting conditions given in (30) is based on (2).
When the percussive cam changes from being idle to traveling, the hammer moves upward along the convex profile and compresses the percussive spring to accumulate energy for next impact load.At the same time, the hammer is no longer pressed to the top of cutting blade's post, and the cutting blade is temporarily free.As the cut continues, the cutting blade is quickly lifted by the uncut surface of the lunar rock simulant and clings to the upper boundary of the tool post.At the same time, the lunar rock simulant is cut and fractured according to the predefined cut depth.

Experimental Verification
The experimental verification of the characteristics of the load for rotary-percussive drilling with an LRCB is divided into two stages.First, the impact-cutting load characteristics are tested using a single-blade linear cutting load test platform to verify the characteristics of the load on the cutting blade in different stages.Next, a rotary-percussive drilling test of the lunar rock simulant is performed via a comprehensive test platform for the characteristics of the drill to compare the variation trends of the drilling load before and after the impact load is applied and to analyze and verify the characteristics of the LRCB drilling load for different drilling parameters.Based on the specifications in the Chang' e plan for experimental verification of the LRCB's ground functionality, marble of grade-VI drillability is selected to simulate lunar rock.

Test Environment
Overview.The single-blade linear cutting load test platform is shown in Figure 8(a).The cutting blade installed in the mobile plate is driven by a servo actuator to make a horizontal cut.The front panel of the mobile plate is connected to the penetrating mechanism to make the cutting blade penetrate vertically.Additionally, the percussive mechanism (which is driven according to the theory shown in Figure 3) is attached to the top of the cutting blade's post to apply an impact to the cutting blade.A magnescale sensor is installed on the left side of the front panel to monitor the cutting blade's penetration displacement and fluctuations in the cut in real-time.A rock holder is installed on the platform rack by means of a six-axis force sensor to measure the cutting and impact loads on the rock in real-time during the cutting/percussion process [27].
The comprehensive drilling/sampling characteristics test platform is shown in Figure 8(b).The bit is installed in a rotary-percussive drill via a 1 m hollow auger.Guide rails are installed on both sides of the drill box and connected to the penetration mechanism via the drive chain.The drill and the penetration mechanism drive the LRCB during rotarypercussive drilling.Tension and torque sensors are installed in the drive chain and along the axis of the auger, respectively, to measure the drilling pressure and the torque on the coring bit in real-time.A magnescale sensor is installed between the drill box and the guide rails to measure the penetration displacement [28].

Experimental Verification of the Cutting Impact Load
Characteristics for Lunar Rock Simulant.According to the theoretical model described in the previous sections, the impact load on the lunar rock simulant is related to the collision contact rigidities of the hammer and cutting blade's post and between the cutting blade and the lunar rock simulant.The lunar rock simulant's impact load transmission characteristics are tested to obtain these two contact rigidity coefficients.The pattern in the lunar rock simulant's variations when it is impacted with different percussion energies is tested to calibrate the two contact rigidities during collision in the model.The variation pattern of the maximum impact load as a function of the percussion energy exerted by the cutting blade on the lunar rock simulant is shown in Figure 9.In the test, the percussive spring pretightening force is adjusted to vary the percussion energy.As the percussion energy increases from 0.0138 J to 0.8027 J, the maximum impact load on the lunar rock simulant increases from 623.2 N to 3,598 N. Using the method of undetermined coefficients, the collision contact rigidities between the hammer and the cutting blade's post and between the cutting blade and the lunar rock simulant are obtained; they are listed in Table 2.
The other parameters of the model are obtained from a test of the rock's mechanical characteristics and a discrete element simulation of the single-blade linear cutting process.
Because the effect of the cut configuration on the transmission loss of the impact load is insignificant, a 4 mm wide cutting blade with a −13 ∘ front rake angle and a 0 ∘ blade rake angle is used for the percussive cutting and fracture of lunar rock simulant via a single blade in the load characteristics test.When the cut depth is 0.033 mm, the cutting velocity is 26.1 mm/s, the percussion energy is 0.28 J, the percussion frequency is 4 Hz, and the cutting pressure load as a function of time for the lunar rock simulant in a single impact cycle is as shown in Figure 10(a).The solid blue line represents the test load.After one cycle of impact loading, the cutting pressure on the lunar rock simulant peaks instantaneously, .When the percussive cam begins to push it, the cutting blade enters the position-control cutting stage (the light grey area).When the pretightening force on the tool post is instantaneously released, the cutting blade is lifted rapidly by the cutting area of the rock and progresses according to the predefined cut depth.The cutting pressure stabilizes at approximately 23.1 N. Because there are assembly gaps and elastic deformation between the rack and the tool post, the rigidity of the system and the overcut depth, Δℎ  , vary during the cutting process.
For different cut depths and percussion frequencies, the average cutting impact load on the lunar rock simulant is shown in Figure 11.Figures 11(a) and 11(b) show curved surfaces representing the variation of the load on the lunar rock simulant during the force-control cutting stage.In this stage, the cutting load is inversely proportional to the percussion frequency.This is mainly because as the percussion frequency increases, the proportion of the fracture load to which the impact-cutting coupling condition applies increases, and a smaller shear fracture angle becomes sufficient to complete the fracture in the lunar rock simulant under the impactcutting coupling condition.Therefore, the average load gradually decreases.Moreover, with the increase in the cut depth, the cutting load becomes clearly divided into two stages.
When the cut depth is less than 0.1 mm, no variation of fracture load on the lunar rock simulant is apparent.This is because when there is a percussive spring pretightening force, a cut with a depth of 0.1 mm is made on the lunar rock simulant; for comparison, when cut depth exceeds 0.1 mm, the fracture load increases significantly with the cut depth, and the percussive spring pretightening force never plays a dominant role.Figures 11(c) and 11(d) show curved surfaces that represent the variation in the load on the lunar rock simulant in the position-control cutting stage.When the cut is relatively shallow (less than 0.1 mm deep), the fracture load on the rock increases with the percussion frequency; as the cut depth increases, the effect of the percussion frequency on the fracture load on the rock decreases.Figures 11(e) and 11(f) show the variation trend of the average impact-cutting load on the lunar rock simulant as a function of the percussion frequency and the cut depth.In general, the fracture load on the lunar rock simulant increases significantly with the cut depth and decreases slightly as the percussion frequency increases.The variation trend of the test results matches that of the theoretical calculation.

Analysis and Experimental Verification of the Load Characteristics of Rotary-Percussive Drilling with an LRCB.
The load characteristics of rotary-percussive drilling with an LRCB for lunar rock simulant are investigated using a comprehensive drill characteristics test platform.The variation trends of the load on the drill bit before and after percussion is applied are compared to analyze the auxiliary effect of the impact load on the fracture of lunar rock simulant.In the test, the LRCB's percussion energy input is set to 2.6 J.
Based on the LRCB's configuration, the percussion energies gained by the different elemental blades in a percussion cycle are 0.196 J, 0.196 J, 0.182 J, and 0.294 J.When the coring bit rotates at 30 rpm, the penetration rate is 3 mm/min, and the percussion frequency is 20 Hz; the variations in the LRCB drilling load before and after percussion is applied are as shown in Figure 12.
In the initial drilling stage, the LRCB penetrates at PPR of 0.1 mm/r, and the percussive mechanism is not activated (light yellow area).During this time, the drilling load gradually increases with the drilling depth, and the rotary torque and the drilling pressure stabilize at 3.71 Nm and 482.3 N, respectively.When the drilling depth reaches 30 mm, the percussive mechanism is activated (light grey area).This shows that, under an impact load, the rotary torque and drilling pressure on the LRCB both decrease and eventually stabilize at 3.08 Nm and 348.6 N, respectively, which shows that the percussive mechanism effectively reduces the drilling load required to penetrate lunar rock simulant.
Figure 13 shows variation trend of the average drilling load of the LRCB for different PPR as the percussion frequency increases from 0 Hz to 20 Hz.The bit's rotational velocity is set to 50 rpm; the PPR during drilling (0.04-0.
Figure 13 shows that, once the impact load has been introduced, the drilling load of the LRCB decreases by different amounts and decreases further as the percussion frequency increases.When the percussion frequency increases to 20 Hz and the PPR is 0.04 mm/r and 0.2 mm/r, the rotary torque required to fracture the lunar rock simulant is reduced by 0.86 Nm and 0.56 Nm, respectively; the drilling pressure required to fracture the lunar rock simulant is decreased by 85.8 N and 116 N, respectively.Additionally, the variation trend of the drilling load as a function of the PPR and the percussion frequency matches the theoretical prediction.
For percussion frequency 20 Hz and PPR 0.1 mm/r, the variation trend of the drilling load of the LRCB is shown in Figure 14.As the rotational velocity increases, both the TOB and the WOB increase.This is because the increase in the bit's rotational velocity decreases the number of impact loads applied to the lunar rock simulant by the bit in each rotation cycle of the drilling process, that is, the decrease in the blow per revolution decreases the proportion of fracture due to impact-cutting coupling for the lunar rock simulant.Therefore, the drill should provide a larger TOB and a higher WOB.

Conclusions
In this paper, the process of penetrating lunar rock simulant via LRCB rotary-percussive drilling is analyzed.The following conclusions are reached: (1) A load prediction model for lunar rock simulant penetration via LRCB rotary-percussive drilling is created.It effectively predicts the characteristics of the load on lunar rock simulant during rotarypercussive drilling within known parameter ranges (the rotational velocity is between 30 and 100 rpm; the penetration rate is between 3 and 20 mm/min; the percussion frequency is between 0 and 20 Hz; and the single percussion energy is 2.6 J) for a fixed drill configuration.These parameter ranges cover the adjustable ranges of the lunar surface sampling mechanism.
(2) The effects of the PPR, the percussion frequency, and the rotational velocity on the lunar rock simulant's drilling load are determined.The PPR is the primary parameter influencing the drilling load.When the PPR is fixed, introducing an impact load effectively decreases the drilling load on the lunar rock simulant.Additionally, the drilling load is inversely proportional to the percussion frequency.When the percussion frequency reaches its upper limit, reducing the rotational velocity increases the blow per revolution, improves the bit's efficiency at penetrating and fracturing rock in a cycle, and further decreases the drilling load on the rock.
To simplify the model, impact load loss induced by the bit's configuration is not considered in this paper.A relevant investigation will be conducted in subsequent research.

Figure 1 :Figure 2 :
Figure 1: Three-dimensional configuration illustration of the LRCB and the actual LRCB.

Figure 3 :
Figure 3: Three-dimensional configuration of the percussive mechanism.

Figure 4 :Figure 5 :
Figure 4: Fracture process for lunar rock simulant with impact-cutting coupling.

Figure 6 :
Figure 6: Stress analysis of the percussive fracture of lunar rock simulant via the cutting blade.

Figure 7 :
Figure 7: Stress analysis for lunar rock simulant fracture via an elemental blade by means of coupled impact and cutting.

Figure 8 :
Figure 8: Single-blade linear cutting load test platform (a) and comprehensive drill characteristics test platform (b).

Figure 9 :
Figure 9: Maximum impact load on the lunar rock simulant with the percussion energy increase.

Figure 10 ( 3 Pen − 41.71ℎ 2
b) shows the variation trend of the impact load on the lunar rock simulant versus time.It shows that the theoretical variation curve for the impact load essentially matches the experimental one.The average penetration distance is defined as the calibration value of the fluctuation in the cut depth, Δℎ  .The fluctuations in the cut depth as a function of the cut depth and the percussion frequency during the force-and positioncontrol cutting stages are obtained via fitting as follows:Δℎ _for = (39450 ⋅ ℎ 2 Pen + 22290 ⋅ ℎ Pen + 6883) × 10 −6 , Δℎ _pos = (−21.11+ 764.8ℎPen + 1.315  − 5924ℎ 2 Pen + 1.119ℎ Pen ⋅   + 20210ℎ Pen ⋅   ) × 10 −4 .

Figure 10 :
Figure 10: Variation trend of the resistant force on the lunar rock simulant in a single impact cycle.

Figure 11 :Figure 12 :Figure 13 :
Figure 11: Trends of the average impact-cutting load on the lunar rock simulant versus the percussion frequency and the cut depth.

Figure 14 :
Figure 14: Trend of the LRCB load in rotary-percussive drilling with the increased rotational velocity and fixed PPR.

Table 2 :
Characteristic parameters in the rotary-percussive drilling load prediction model for lunar rock simulant.