Stability analysis of intermittently jointed rock slopes based on the stepped failure mode

In practical engineering,due to the noncontinuity characteristics of joints in rock slopes,in addition to plane failure,stepped sliding failure may occur for intermittently jointed rock slopes.Especially for intermittently bedding jointed rock slopes,the correlation and difference in strength parameters between joints and rock bridges,along with the various failure modes and intermittency of rock bridges,contribute to the complexity of stepped failure modes and the unpredictability of failure regions.Based on the upper-bound limit analysis method and multi-sliders step-path failure mode,considering the shear and tensile failure of rock bridges and the weakened relationship between the strength parameters of rock bridges and jointed surfaces,by introducing the modified M-C failure criterion and the formula for calculating the energy consumption of tensile failure of rock bridges,two failure mechanisms are constructed to obtain the safety factor(F_(s))of intermittently jointed rock slopes.The sequential quadratic programming method is used to obtain the optimal upper-bound solution for F_(s).The influence of multiple key parameters(slope height H,horizontal distance L,Slope angleβ,shear strength parameters of the rock bridgeφr and cr,Dimensionless parameter u,weakening coefficients of the internal friction angle and cohesion between the rock bridges and joint surfaces Kφand Kc)on the stability analysis of intermittently jointed rock slopes under the shear failure mode of rock bridges as well as under the tensile failure mode is also explored.The reliability of the failure mechanisms is verified by comparative analysis with theoretical results,numerical results,and landslide cases,and the variation rules of F_(s)with each key parameter are obtained.The results show that F_(s) varies linearly withφr and cr of the rock bridge and with K_(φ)and K_(c),whereas F_(s)changes nonlinearly with H and L.In particular,with the increase in Kφand Kc,Fs increases by approximately 52.78%and 171.02%on average,respectively.For rock bridge tensile failure,F_(s) shows a nonlinearly positive correlation withφr,cr,Kφand Kc.In particular,with the increase in Kφand Kc,Fs increases by approximately 13%and 61.69%on average,respectively.Fs decreases rapidly with increasing slope gradientβand decreasing dimensionless parameterμ.When Kφand Kc are both less than 1.0,the stepped sliding surface occurs more easily than the plane failure surface,especially in the case of tensile failure of the rock bridge.In addition,rock slopes with higher strength parameters,taller heights,and greater weakening coefficients are prone to rock bridge tension failure with lower Fs,and more attention should be given to the occurrence of such accidents in actual engineering.

A 3D microseismic data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions and its application

Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.

Strength and deformation characteristics of irregular columnar jointed rock mass: A combined experimental and theoretical study

The irregularity of jointed network poses a challenge to the determination of field mechanical param-eters of columnar jointed rock mass(CJRM),and a reasonable prediction of deformation and strength characteristics of CJRM is important for engineering construction.The Voronoi diagram and three-dimensional printing technology were used to make an irregular columnar jointed mold,and the irregular CJRM(ICJRM)specimens with different dip directions and dip angles were prepared.Uniaxial compression tests were performed,and the anisotropic strength and deformation characteristics of ICJRM were described.The failure modes and mechanisms were revealed in accordance with the final appearances of the ICJRM specimens.Based on the model test results,the empirical correlations for determining the field deformation and strength parameters of CJRM were derived using the dip angle and modified joint factor.The proposed empirical equations were used in the Baihetan Project,and the calculated mechanical parameters were compared with the field test results and those obtained from the tunneling quality index method.Results showed that the deformation parameters determined by the two proposed methods are all consistent with the field test results,and these two methods can also estimate the strength parameters effectively.

Fracture behaviors of columnar jointed rock mass using interface mechanics theorem

For a special geological structure of columnar jointed rock mass(CJRM),its mechanical properties are strongly affected by the columnar joints.To describe the fracture behaviors of CJRM using the basic theories of interface mechanics for composite materials,the interface stresses of the vertical and horizontal joints,which are the two primary joints in the CJRM under triaxial compression,are studied,and their mathematical expressions are derived based on the superposition principle.Based on the obtained interface stresses of the vertical and horizontal joints in the CJRM,the crack initiation of the joint interface in the CJRM is studied using the maximum circumferential stress theory in fracture mechanics.Moreover,based on this investigation,the fracture behaviors of CJRM are analyzed.According to the results of similar material physical model tests for the CJRM,the theoretical study is verified.Finally,the influence of the mechanical parameters of the CJRM on the joint interface stress is discussed comprehensively.

Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Rock mass quality prediction on tunnel faces with incomplete multi-source dataset via tree-augmented naive Bayesian network

Rock mass quality serves as a vital index for predicting the stability and safety status of rock tunnel faces.In tunneling practice,the rock mass quality is often assessed via a combination of qualitative and quantitative parameters.However,due to the harsh on-site construction conditions,it is rather difficult to obtain some of the evaluation parameters which are essential for the rock mass quality prediction.In this study,a novel improved Swin Transformer is proposed to detect,segment,and quantify rock mass characteristic parameters such as water leakage,fractures,weak interlayers.The site experiment results demonstrate that the improved Swin Transformer achieves optimal segmentation results and achieving accuracies of 92%,81%,and 86%for water leakage,fractures,and weak interlayers,respectively.A multisource rock tunnel face characteristic(RTFC)dataset includes 11 parameters for predicting rock mass quality is established.Considering the limitations in predictive performance of incomplete evaluation parameters exist in this dataset,a novel tree-augmented naive Bayesian network(BN)is proposed to address the challenge of the incomplete dataset and achieved a prediction accuracy of 88%.In comparison with other commonly used Machine Learning models the proposed BN-based approach proved an improved performance on predicting the rock mass quality with the incomplete dataset.By utilizing the established BN,a further sensitivity analysis is conducted to quantitatively evaluate the importance of the various parameters,results indicate that the rock strength and fractures parameter exert the most significant influence on rock mass quality.

Analysis of the Elastic-Plastic Dynamic Finite Element on a Jointed Rock Mass

Aim To study the elastic plastic dynamical constitutive relations about a jointed rock mass under explosion load and its computer simulation. Methods\ Stress history is taken into account and stresses will follow changes in time during a period of explosion load. According to the principle of static force balance, the corresponding nodal concentrated force is calculated and the nodal displacement is counted. The elastic plastic dynamic finite element equations are thus obtained. Results\ A finite element method is given for a jointed rock mass under explosion load. Conclusion\ The problem of large plastic deformation for jointed rock mass on blasting was efficiently resolved through dynamic finite element analysis and the range of damages by blasting simulated, and this pushes forward the problem to engineering practice.

Effect of eccentric and inclined loading on the bearing capacity of strip footing placed on rock mass

This paper deals with the bearing capacity determination of strip footing on a rock mass in hilly area by considering the influence of inclined and eccentric loading. Applying the generalized HoekBrown failure criterion, the failure behavior of the rock mass is modeled with the help of the power cone programming in the lower bound finite element limit analysis framework. Using bearing capacity factor(Ns), the change in bearing capacity of the strip footing due to the occurrence of eccentrically inclined loading is presented. The variations of the magnitude of Ns are obtained by examining the effects of the Hoek-Brown rock mass strength parameters(uniaxial compressive strength(sci), disturbance factor(D), rock parameter(mi), and Geological Strength Index(GSI)) in the presence of different magnitudes of eccentricity(e) and inclination angle(λ) with respect to the vertical plane, and presented as design charts. Both the inclined loading modes, i.e., inclination towards the center of strip footing(+λ) and inclination away from the center of strip footing(-λ), are adopted to perform the investigation. In addition, the correlation between the input parameters and the corresponding output is developed by utilizing the artificial neural network(ANN). Additionally, from sensitivity analysis, it is observed that inclination angle(λ) is the most sensitive parameter. For practicing engineers, the obtained design equation and design charts can be beneficial to understand the bearing capacity variation in the existence of eccentrically inclined loading in mountain areas.

Evaluation of slope stability through rock mass classification and kinematic analysis of some major slopes along NH-1A from Ramban to Banihal, North Western Himalayas

The network of Himalayan roadways and highways connects some remote regions of valleys or hill slopes,which is vital for India’s socio-economic growth.Due to natural and artificial factors,frequency of slope instabilities along the networks has been increasing over last few decades.Assessment of stability of natural and artificial slopes due to construction of these connecting road networks is significant in safely executing these roads throughout the year.Several rock mass classification methods are generally used to assess the strength and deformability of rock mass.This study assesses slope stability along the NH-1A of Ramban district of North Western Himalayas.Various structurally and non-structurally controlled rock mass classification systems have been applied to assess the stability conditions of 14 slopes.For evaluating the stability of these slopes,kinematic analysis was performed along with geological strength index(GSI),rock mass rating(RMR),continuous slope mass rating(CoSMR),slope mass rating(SMR),and Q-slope in the present study.The SMR gives three slopes as completely unstable while CoSMR suggests four slopes as completely unstable.The stability of all slopes was also analyzed using a design chart under dynamic and static conditions by slope stability rating(SSR)for the factor of safety(FoS)of 1.2 and 1 respectively.Q-slope with probability of failure(PoF)1%gives two slopes as stable slopes.Stable slope angle has been determined based on the Q-slope safe angle equation and SSR design chart based on the FoS.The value ranges given by different empirical classifications were RMR(37-74),GSI(27.3-58.5),SMR(11-59),and CoSMR(3.39-74.56).Good relationship was found among RMR&SSR and RMR&GSI with correlation coefficient(R 2)value of 0.815 and 0.6866,respectively.Lastly,a comparative stability of all these slopes based on the above classification has been performed to identify the most critical slope along this road.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

A photogrammetric approach for quantifying the evolution of rock joint void geometry under varying contact states

Accurate measurement of the evolution of rock joint void geometry is essential for comprehending the distribution characteristics of asperities responsible for shear and seepage behaviors.However,existing techniques often require specialized equipment and skilled operators,posing practical challenges.In this study,a cost-effective photogrammetric approach is proposed.Particularly,local coordinate systems are established to facilitate the alignment and precise quantification of the relative position between two halves of a rock joint.Push/pull tests are conducted on rock joints with varying roughness levels to induce different contact states.A high-precision laser scanner serves as a benchmark for evaluating the photogrammetry method.Despite certain deviations exist,the measured evolution of void geometry is generally consistent with the qualitative findings of previous studies.The photogrammetric measurements yield comparable accuracy to laser scanning,with maximum errors of 13.2%for aperture and 14.4%for void volume.Most joint matching coefficient(JMC)measurement errors are below 20%.Larger measurement errors occur primarily in highly mismatched rock joints with JMC values below 0.2,but even in cases where measurement errors exceed 80%,the maximum JMC error is only 0.0434.Thus,the proposed photogrammetric approach holds promise for widespread application in void geometry measurements in rock joints.

A statistical damage-based constitutive model for shearing of rock joints in brittle drop mode

Some rock joints exhibit significant brittleness,characterized by a sharp decrease in shear stress upon reaching the peak strength.However,existing models often fail to accurately represent this behavior and are encumbered by numerous parameters lacking clear mechanical significance.This study presents a new statistical damage constitutive model rooted in both damage mechanics and statistics,containing only three model parameters.The proposed model encompasses all stages of joint shearing,including the compaction stage,linear stage,plastic yielding stage,drop stage,strain softening stage,and residual strength stage.To derive the analytical expression of the constitutive model,three boundary conditions are introduced.Experimental data from both natural and artificial rock joints is utilized to validate the model,resulting in average absolute relative errors ranging from 3%to 8%.Moreover,a comparative analysis with established models illustrates that the proposed model captures stress drop and post-peak strain softening more effectively,with model parameters possessing clearer mechanical interpretations.Furthermore,parameter analysis is conducted to investigate the impacts of model parameters on the curves and unveil the relationship between these parameters and the mechanical properties of rock joints.Importantly,the proposed model is straightforward in form,and all model parameters can be obtained from direct shear tests,thus facilitating the utilization in numerical simulations.

On the calibration of a shear stress criterion for rock joints to represent the full stress-strain profile

Conventional numerical solutions developed to describe the geomechanical behavior of rock interfaces subjected to differential load emphasize peak and residual shear strengths.The detailed analysis of preand post-peak shear stress-displacement behavior is central to various time-dependent and dynamic rock mechanic problems such as rockbursts and structural instabilities in highly stressed conditions.The complete stress-displacement surface(CSDS)model was developed to describe analytically the pre-and post-peak behavior of rock interfaces under differential loads.Original formulations of the CSDS model required extensive curve-fitting iterations which limited its practical applicability and transparent integration into engineering tools.The present work proposes modifications to the CSDS model aimed at developing a comprehensive and modern calibration protocol to describe the complete shear stressdisplacement behavior of rock interfaces under differential loads.The proposed update to the CSDS model incorporates the concept of mobilized shear strength to enhance the post-peak formulations.Barton’s concepts of joint roughness coefficient(JRC)and joint compressive strength(JCS)are incorporated to facilitate empirical estimations for peak shear stress and normal closure relations.Triaxial/uniaxial compression test and direct shear test results are used to validate the updated model and exemplify the proposed calibration method.The results illustrate that the revised model successfully predicts the post-peak and complete axial stressestrain and shear stressedisplacement curves for rock joints.

Rock Mass Quality Rating Based on the Multi-Criteria Grey Metric Space

Assessment of rock mass quality significantly impacts the design and construction of underground and open-pit mines from the point of stability and economy.This study develops the novel Gromov-Hausdorff distance for rock quality(GHDQR)methodology for rock mass quality rating based on multi-criteria grey metric space.It usually presents the quality of surrounding rock by classes(metric spaces)with specified properties and adequate interval-grey numbers.Measuring the distance between surrounding rock sample characteristics and existing classes represents the core of this study.The Gromov-Hausdorff distance is an especially useful discriminant function,i.e.,a classifier to calculate these distances,and assess the quality of the surrounding rock.The efficiency of the developed methodology is analyzed using the Mean Absolute Percentage Error(MAPE)technique.Seven existing methods,such as the Gaussian cloud method,Discriminant method,Mutation series method,Artificial neural network(ANN),Support vector machine(SVM),Grey wolf optimizer and Support vector classification method(GWO-SVC)and Rock mass rating method(RMR)are used for comparison with the proposed GHDQR method.The share of the highly accurate category of 85.71%clearly indicates compliance with actual values obtained by the compared methods.The results of comparisons showed that the model enables objective,efficient,and reliable assessment of rock mass quality.

Optimizing profile line interval for enhanced accuracy in rock joint morphology and shear strength assessments

2D profile lines play a critical role in cost-effectively evaluating rock joint properties and shear strength.However, the interval(ΔI_(L)) between these lines significantly impacts roughness and shear strength assessments. A detailed study of 45 joint samples using four statistical measures across 500 different ΔI_(L)values identified a clear line interval effect with two stages: stable and fluctuation-discrete.Further statistical analysis showed a linear relationship between the error bounds of four parameters,shear strength evaluation, and their corresponding maximum ΔI_(L)values, where the gradient k of this linear relationship was influenced by the basic friction angle and normal stress. Accounting for these factors,lower-limit linear models were employed to determine the optimal ΔI_(L)values that met error tolerances(1%–10%) for all metrics and shear strength. The study also explored the consistent size effect on joints regardless of ΔI_(L)changes, revealing three types of size effects based on morphological heterogeneity.Notably, larger joints required generally higher ΔI_(L)to maintain the predefined error limits, suggesting an increased interval for large joint analyses. Consequently, this research provides a basis for determining the optimal ΔI_(L), improving accuracy in 2D profile line assessments of joint characteristics.

A Novel ISSA–DELM Model for Predicting Rock Mass Permeability

In pumped storage projects,the permeability of rock masses is a crucial parameter in engineering design and construction.The rock mass permeability coefficient(K)is influenced by various geological parameters,and previous studies aimed to establish an accurate relationship between K and geological parameters.This study uses the improved sparrow search algorithm(ISSA)to optimize the parameter settings of the deep extreme learning machine(DELM),constructing a prediction model with flexible parameter selection and high accuracy.First,the Spearman method is applied to analyze the correlation between geological parameters.A sample database is built by comprehensively selecting four geological indexes:burial depth,rock quality designation(RQD),fracture density characteristic index(FD),and rock mass integrity designation(RID).Hence,the defects of the sparrow search algorithm(SSA)are enhanced using the improved strategy,and the initial input weights of the DELM are optimized.Finally,the proposed ISSA–DELM model is employed to predict the permeability coefficient of rock mass in the entire study area.The results showed that the predictive performance of the model is superior to that of the DELM and SSA–DELM.Therefore,this model successfully provides insights into the distribution characteristics of rock mass permeability at engineering sites.

Inverting the rock mass P-wave velocity field ahead of deep buried tunnel face while borehole drilling

Imaging the wave velocity field surrounding a borehole while drilling is a promising and urgently needed approach for extending the exploration range of the borehole point.This paper develops a drilling process detection(DPD)system consisting of a multifunctional sensor and a pilot geophone installed at the top of the drilling rod,geophones at the tunnel face,a laser rangefinder,and an onsite computer.A weighted adjoint-state first arrival travel time tomography method is used to invert the P-wave velocity field of rock mass while borehole drilling.A field experiment in the ongoing construction of a deep buried tunnel in southwestern China demonstrated the DPD system and the tomography method.Time-frequency analysis of typical borehole drilling detection data shows that the impact drilling source is a pulse-like seismic exploration wavelet.A velocity field of the rock mass in a triangular area defined by the borehole trajectory and geophone receiving line can be obtained.Both the borehole core and optical image validate the inverted P-wave velocity field.A numerical simulation of a checkerboard benchmark model is used to test the tomography method.The rapid convergence of the misfits and consistent agreement between the inverted and observed travel times validate the P-wave velocity imaging.

Load-bearing characteristics and energy evolution of fractured rock masses after granite and sandstone grouting

Experiments on grouting-reinforced rock mass specimens with different particle sizes and features were carried out in this study to examine the effects of grouting reinforcement on the load-bearing characteristics of fractured rock mass.The strength and deformation features of grouting-reinforced rock mass were analyzed under different loading manners;the energy evolution mechanism of grouting-reinforced rock mass specimens with different particle sizes and features was investigated;the energy dissipation ratio and post-peak stress decreasing rate were employed to evaluate the bearing stability of grouting-reinforced rock mass.The results show that the strength and ductility of granite-reinforced rock mass(GRM)under biaxial loading are higher than that of sandstone-reinforced rock mass(SRM)under uniaxial loading.Besides,the energy evolution characteristics of grouting-reinforced rock mass under uniaxial and biaxial loading mainly could be divided into early,middle,and late stages.In the early stage,total,elastic,and dissipation energies were quite small with flatter curves;in the middle stage,elastic energy increased rapidly,whereas dissipation energy increased slowly;in the late stage,dissipation energy increased sharply.The energy dissipation ratio was used to represent the pre-peak plastic deformation.Under uniaxial loading,this ratio increased as the particle size increased and the pre-peak plastic deformation of grouting-reinforced rock mass became larger;under biaxial loading,it dropped as the particle size increased,and the pre-peak plastic deformation of grouting-reinforced rock mass became smaller.The post-peak stress decline rate A_(v) was used to assess the post-peak bearing performance of grouting-reinforced rock mass.Under uniaxial loading,parameter A_(v) exhibited reduction as the particle size kept increasing,and the ability of post-peak of grouting-reinforced rock mass to allow deformation development was greater,and the bearing capacity was greater;under biaxial loading,A_(v) increased with the particle size,and the ability of post-peak of grouting-reinforced rock mass to allow deformation development was low and the bearing capacity was reduced.The findings are considered instrumental in improving the stability of the roadway-surrounding rock by granite and sandstone grouting.

Failure behavior and strength model of blocky rock mass with and without rockbolts

To better understand the failure behaviours and strength of bolt-reinforced blocky rocks,large scale extensive laboratory experiments are carried out on blocky rock-like specimens with and without rockbolt reinforcement.The results show that both shear failure and tensile failure along joint surfaces are observed but the shear failure is a main controlling factor for the peak strength of the rock mass with and without rockbolts.The rockbolts are necked and shear deformation simultaneously happens in bolt reinforced rock specimens.As the joint dip angle increases,the joint shear failure becomes more dominant.The number of rockbolts has a significant impact on the peak strain and uniaxial compressive strength(UCS),but little influence on the deformation modulus of the rock mass.Using the Winkler beam model to represent the rockbolt behaviours,an analytical model for the prediction of the strength of boltreinforced blocky rocks is proposed.Good agreement between the UCS values predicted by proposed model and obtained from experiments suggest an encouraging performance of the proposed model.In addition,the performance of the proposed model is further assessed using published results in the literature,indicating the proposed model can be used effectively in the prediction of UCS of bolt-reinforced blocky rocks.

Deformation characteristics and damage ontologies of soft and hard composite rock masses under impact loading

As one of the most common occurring geological landforms in deep rock formations, the dynamic mechanical properties of layered composite rock bodies under impact loading have been widely studied by scholars. To study the dynamic properties of soft and hard composite rocks with different thickness ratios, this paper utilizes cement, quartz sand and gypsum powder to construct soft and hard composite rock specimens and utilizes a combination of indoor tests, numerical calculations, and theoretical analyses to investigate the mechanical properties of soft and hard composite rock bodies. The test results reveal that:(1) When the proportion of hard rock increases from 20% to 50%, the strength of the combined rock body increases by 69.14 MPa and 87 MPa when the hard rock face and soft rock face are loaded, respectively;however, when the proportion of hard rock is the same, the compressive strength of the hard rock face impact is 9%-17% greater than that of the soft rock face impact;(2) When a specimen of soft and hard combined rock body is subjected to impact loading, the damage mode involves mixed tension and shear damage, and the cracks generally first appear at the ends of the specimen, then develop on the laminar surface from the impact surface, and finally end in the overall damage of the soft rock part. The development rate and the total number of cracks in the same specimen when the hard rock face is impacted are significantly greater than those when the soft rock face is impacted;(3) By introducing Weibull’s statistical strength theory to establish the damage variables of soft-hard combined rock bodies, combined with the DP strength criterion, the damage model and the Kelvin body are concatenated to obtain a statistical damage constitutive model, which can better fit the full stress-strain curve of soft-hard combined rock body specimens under a single impact load.