Estimating shear strength of high-level pillars supported with cemented backfilling using the HoekeBrown strength criterion

Deep metal mines are often mined using the high-level pillars with subsequent cementation backfilling(HLSCB)mining method.At the design stage,it is therefore important to have a reasonable method for determining the shear strength of the high-level pillars(i.e.cohesion and internal friction angle)when they are supported by cemented backfilling.In this study,a formula was derived for the upper limit of the confining pressure σ3max on a high-level pillar supported by cemented backfilling in a deep metal mine.A new method of estimating the shear strength of such pillars was then proposed based on the Hoek eBrown failure criterion.Our analysis indicates that the horizontal stress σhh acting on the cemented backfill pillar can be simplified by expressing it as a constant value.A reasonable and effective value for σ3max can then be determined.The value of s3max predicted using the proposed method is generally less than 3 MPa.Within this range,the shear strength of the high-level pillar is accurately calculated using the equivalent MohreCoulomb theory.The proposed method can effectively avoid the calculation of inaccurate shear strength values for the high-level pillars when the original HoekeBrown criterion is used in the presence of large confining pressures,i.e.the situation in which the cohesion value is too large and the friction angle is too small can effectively be avoided.The proposed method is applied to a deep metal mine in China that is being excavated using the HLSCB method.The shear strength parameters of the high-level pillars obtained using the proposed method were input in the numerical simulations.The numerical results show that the recommended level heights and sizes of the high-level pillars and rooms in the mine are rational.

Anisotropic time-dependent behaviors of shale under direct shearing and associated empirical creep models

Understanding the anisotropic creep behaviors of shale under direct shearing is a challenging issue.In this context,we conducted shear-creep and steady-creep tests on shale with five bedding orientations (i.e.0°,30°,45°,60°,and 90°),under multiple levels of direct shearing for the first time.The results show that the anisotropic creep of shale exhibits a significant stress-dependent behavior.Under a low shear stress,the creep compliance of shale increases linearly with the logarithm of time at all bedding orientations,and the increase depends on the bedding orientation and creep time.Under high shear stress conditions,the creep compliance of shale is minimal when the bedding orientation is 0°,and the steady-creep rate of shale increases significantly with increasing bedding orientations of 30°,45°,60°,and 90°.The stress-strain values corresponding to the inception of the accelerated creep stage show an increasing and then decreasing trend with the bedding orientation.A semilogarithmic model that could reflect the stress dependence of the steady-creep rate while considering the hardening and damage process is proposed.The model minimizes the deviation of the calculated steady-state creep rate from the observed value and reveals the behavior of the bedding orientation's influence on the steady-creep rate.The applicability of the five classical empirical creep models is quantitatively evaluated.It shows that the logarithmic model can well explain the experimental creep strain and creep rate,and it can accurately predict long-term shear creep deformation.Based on an improved logarithmic model,the variations in creep parameters with shear stress and bedding orientations are discussed.With abovementioned findings,a mathematical method for constructing an anisotropic shear creep model of shale is proposed,which can characterize the nonlinear dependence of the anisotropic shear creep behavior of shale on the bedding orientation.

Failure transition of shear-to-dilation band of rock salt under triaxial stresses

Great potential of underground gas/energy storage in salt caverns seems to be a promising solution to support renewable energy.In the underground storage method,the operating cycle unfortunately may reach up to daily or even hourly,which generates complicated pressures on the salt cavern.Furthermore,the mechanical behavior of rock salt may change and present distinct failure characteristics under different stress states,which affects the performance of salt cavern during the time period of full service.To reproduce a similar loading condition on the cavern surrounding rock mass,the cyclic triaxial loading/unloading tests are performed on the rock salt to explore the mechanical transition behavior and failure characteristics under different confinement.Experimental results show that the rock salt samples pre-sent a diffused shear failure band with significant bulges at certain locations in low confining pressure conditions(e.g.5 MPa,10 MPa and 15 MPa),which is closely related to crystal misorientation and grain boundary sliding.Under the elevated confinement(e.g.20 MPa,30 MPa and 40 MPa),the dilation band dominates the failure mechanism,where the large-size halite crystals are crushed to be smaller size and new pores are developing.The failure transition mechanism revealed in the paper provides additional insight into the mechanical performance of salt caverns influenced by complicated stress states.

Rheological properties and concentration evolution of thickened tailings under the coupling effect of compression and shear

Cemented paste backfill(CPB)is a key technology for green mining in metal mines,in which tailings thickening comprises the primary link of CPB technology.However,difficult flocculation and substandard concentrations of thickened tailings often occur.The rheological properties and concentration evolution in the thickened tailings remain unclear.Moreover,traditional indoor thickening experiments have yet to quantitatively characterize their rheological properties.An experiment of flocculation condition optimization based on the Box-Behnken design(BBD)was performed in the study,and the two response values were investigated:concentration and the mean weighted chord length(MWCL)of flocs.Thus,optimal flocculation conditions were obtained.In addition,the rheological properties and concentration evolution of different flocculant dosages and ultrafine tailing contents under shear,compression,and compression-shear coupling experimental conditions were tested and compared.The results show that the shear yield stress under compression and compression-shear coupling increases with the growth of compressive yield stress,while the shear yield stress increases slightly under shear.The order of shear yield stress from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Under compression and compression-shear coupling,the concentration first rapidly increases with the growth of compressive yield stress and then slowly increases,while concentration increases slightly under shear.The order of concentration from low to high under different thickening conditions is shear,compression,and compression-shear coupling.Finally,the evolution mechanism of the flocs and drainage channels during the thickening of the thickened tailings under different experimental conditions was revealed.

Transient response of doubly-curved bio-inspired composite shells resting on viscoelastic foundation subject to blast load using improved first-order shear theory and isogeometric approach

Investigating natural-inspired applications is a perennially appealing subject for scientists. The current increase in the speed of natural-origin structure growth may be linked to their superior mechanical properties and environmental resilience. Biological composite structures with helicoidal schemes and designs have remarkable capacities to absorb impact energy and withstand damage. However, there is a dearth of extensive study on the influence of fiber redirection and reorientation inside the matrix of a helicoid structure on its mechanical performance and reactivity. The present study aimed to explore the static and transient responses of a bio-inspired helicoid laminated composite(B-iHLC) shell under the influence of an explosive load using an isomorphic method. The structural integrity of the shell is maintained by a viscoelastic basis known as the Pasternak foundation, which encompasses two coefficients of stiffness and one coefficient of damping. The equilibrium equations governing shell dynamics are obtained by using Hamilton's principle and including the modified first-order shear theory,therefore obviating the need to employ a shear correction factor. The paper's model and approach are validated by doing numerical comparisons with respected publications. The findings of this study may be used in the construction of military and civilian infrastructure in situations when the structure is subjected to severe stresses that might potentially result in catastrophic collapse. The findings of this paper serve as the foundation for several other issues, including geometric optimization and the dynamic response of similar mechanical structures.

Anisotropic shearing mechanism of Kangding slate:Experimental investigation and numerical analysis

The shear mechanical behavior is regarded as an essential factor affecting the stability of the surrounding rocks in underground engineering.The shear strength and failure mechanisms of layered rock are significantly affected by the foliation angles.Direct shear tests were conducted on cubic slate samples with foliation angles of 0°,30°,45°,60°,and 90°.The effect of foliation angles on failure patterns,acoustic emission(AE)characteristics,and shear strength parameters was analyzed.Based on AE characteristics,the slate failure process could be divided into four stages:quiet period,step-like increasing period,dramatic increasing period,and remission period.A new empirical expression of cohesion for layered rock was proposed,which was compared with linear and sinusoidal cohesion expressions based on the results made by this paper and previous experiments.The comparative analysis demonstrated that the new expression has better prediction ability than other expressions.The proposed empirical equation was used for direct shear simulations with the combined finite-discrete element method(FDEM),and it was found to align well with the experimental results.Considering both computational efficiency and accuracy,it was recommended to use a shear rate of 0.01 m/s for FDEM to carry out direct shear simulations.To balance the relationship between the number of elements and the simulation results in the direct shear simulations,the recommended element size is 1 mm.

The Influence of the Closure of the East Paleo-Tethys Ocean on Southern South China:Evidences from Kinematics and^(40)Ar/^(39)Ar Geochronology of the Rongxian Ductile Shear Zone in Southeastern Guangxi

The Triassic was a crucial period in the tectonic evolution of the South China Block.Research on tectonic deformation during this period provides information on intracontinental orogenic mechanisms in South China.In this study,alongside thermochronological analyses,we examine the macroscopic and microscopic structural features of the Rongxian ductile shear zone,located south of the Darongshan granite in the southeastern part of Guangxi Province,on the southern margin of South China.Sinistral shear is indicated by the characteristics of rotatedσ-type feldspar porphyroclasts,stretching lineations defined by elongated quartz grains and the orientations of quartz c-axes.LA-ICP-MS U-Pb dating of zircons from two samples of granitic mylonite and one of granite yielded ages of ca.256 Ma.Furthermore,two samples of granitic mylonite yield muscovite^(40)Ar/^(39)Ar plateau ages of 249-246 Ma.These results indicate that the Rongxian ductile shear zone resulted from Early Triassic deformation of the late Permian Darongshan granite.This deformation was likely related to the closure of the eastern Paleo-Tethys Ocean and the subsequent collision of the South China and Indochina blocks,during the early stage of the Indosinian orogeny.

Experimental investigation on shear strength deterioration at the interface between different rock types under cyclic loading

The shear strength deterioration of bedding planes between different rock types induced by cyclic loading is vital to reasonably evaluate the stability of soft and hard interbedded bedding rock slopes under earthquake;however,rare work has been devoted to this subject due to lack of attention.In this study,experimental investigations on shear strength weakening of discontinuities with different joint wall material(DDJM)under cyclic loading were conducted by taking the interface between siltstone and mudstone in the Shaba slope of Yunnan Province,China as research objects.A total of 99 pairs of similar material samples of DDJM(81 pairs)and discontinuities with identical joint wall material(DIJM)(18 pairs)were fabricated by inserting plates,engraved with typical surface morphology obtained by performing three-dimensional laser scanning on natural DDJMs sampled from field,into mold boxes.Cyclic shear tests were conducted on these samples to study their shear strength changes with the cyclic number considering the effects of normal stress,joint surface morphology,shear displacement amplitude and shear rate.The results indicate that the shear stress vs.shear displacement curves under each shear cycle and the peak shear strength vs.cyclic number curves of the studied DDJMs are between those of DIJMs with siltstone and mudstone,while closer to those of DIJMs with mudstone.The peak shear strengths of DDJMs exhibit an initial rapid decline followed by a gradual decrease with the cyclic number and the decrease rate varies from 6%to 55.9%for samples with varied surface morphology under different testing conditions.The normal stress,joint surface morphology,shear displacement amplitude and shear rate collectively influence the shear strength deterioration of DDJM under cyclic shear loading,with the degree of influence being greater for larger normal stress,rougher surface morphology,larger shear displacement amplitude and faster shear rate.

State-of-the-art on the anchorage performance of rock bolts subjected to shear load

Rock bolts are extensively utilized in underground engineering as a means of offering support and stability to rock masses in tunnels,mines,and other underground structures.In environments of high ground stress,faults or weak zones can frequently arise in rock formations,presenting a significant challenge for engineering and potentially leading to underground engineering collapse.Rock bolts serve as a crucial structural element for the transmission of tensile stress and are capable of withstanding shear loads to prevent sliding of weak zones within rock mass.Therefore,a complete understanding of the behavior of rock bolts subjected to shear loads is essential.This paper presents a state-of-the-art review of the research progress of rock bolts subjected to shear load in three categories:experiment,numerical simulation,and analytical model.The review focuses on the research studies and developments in this area since the 1970s,providing a comprehensive overview of numerous factors that influence the anchorage performance of rock bolts.These factors include the diameter and angle of the rock bolt installation,rock strength,grouting material,bolt material,borehole diameter,rock bolt preload,normal stress,joint surface roughness and joint expansion angle.The paper reviews the improvement of mechanical parameter setting in numerical simulation of rock bolt shear.Furthermore,it delves into the optimization of the analytical model concerning rock bolt shear theory,approached from the perspectives of both Elastic foundation beam theory coupled with Elastoplasticity theory and Structural mechanic methods.The significance of this review lies in its ability to provide insights into the mechanical behavior of rock bolts.The paper also highlights the limitations of current research and guidelines for further research of rock bolts.

Mechanical behaviors of backfill-rock composites: Physical shear test and back-analysis

The shear behavior of backfill-rock composites is crucial for mine safety and the management of surface subsidence.For exposing the shear failure mechanism of backfill-rock composites,we conducted shear tests on backfill-rock composites under three constant normal loads,compared with the unfilled rock.To investigate the macro-and meso-failure characteristics of the samples in the shear tests,the cracking behavior of samples was recorded by a high-speed camera and acoustic emission monitoring.In parallel with the experimental test,the numerical models of backfill-rock composites and unfilled rock were established using the discrete element method to analyze the continuous-discontinuous shearing process.Based on the damage mechanics and statistics,a novel shear constitutive model was proposed to describe mechanical behavior.The results show that backfill-rock composites had a special bimodal phenomenon of shearing load-deformation curve,i.e.the first shearing peak corresponded to rock break and the second shearing peak induced by the broken of aeolian sand-cement/fly ash paste backfill.Moreover,the shearing characteristic curves of the backfill-rock composites could be roughly divided into four stages,i.e.the shear failure of the specimens experienced:stage I:stress concentration;stage II:crack propagation;stage III:crack coalescence;stage IV:shearing friction.The numerical simulation shows that the existence of aeolian sand-cement/fly ash paste backfill inevitably altered the coalescence type and failure mode of the specimens and had a strengthening effect on the shear strength of backfillrock composites.Based on damage mechanics and statistics,a shear constitutive model was proposed to describe the shear fracture characteristics of specimens,especially the bimodal phenomenon.Finally,the micro-and meso-mechanisms of shear failure were discussed by combining the micro-test and numerical results.The research can advance the better understanding of the shear behavior of backfill-rock composites and contribute to the safety of mining engineering.

Shear band evolution and acoustic emission characteristics of sandstone containing non-persistent flaws

Direct shear tests were conducted on sandstone specimens under different constant normal stresses to study the coalescence of cracks between non-persistent flaws and the shear sliding characteristics of the shear-formed fault.Digital image correlation and acoustic emission(AE)techniques were used to monitor the evolution of shear bands at the rock bridge area and microcracking behaviors.The experimental results revealed that the shear stresses corresponding to the peak and sub-peak in the stressdisplacement curve are significantly affected by the normal stress.Strain localization bands emerged at both the tip of joints and the rock bridge,and their extension and interaction near the peak stress caused a surge in the AE hit rate and a significant decrease in the AE b value.Short and curvilinear strain bands were detected at low normal stress,while high normal stress generally led to more microcracking events and longer coplanar cracks at the rock bridge area.Furthermore,an increase in normal stress resulted in a higher AE count rate and more energetic AE events during friction sliding along the shearformed fault.It was observed that the elastic energy released during the crack coalescence at the prepeak stage was much greater than that released during friction sliding at the post-peak stage.More than 75%of AE events were located in the low-frequency band(0e100 kHz),and this proportion continued to rise with increasing normal stress.Moreover,more AE events of low AF value and high RA value were observed in specimens subjected to high normal stress,indicating that greater normal stress led to more microcracks of shear nature.

Seismic performance evaluation of hybrid coupled shear wall system with shear and flexural fuse-type steel coupling beams

Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.

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.

Photoacoustic elastography based on laser-excited shear wave

Elastography can be used as a diagnostic method for quantitative characterization of tissue hardness information and thus,differential changes in pathophysiological states of tissues.In this study,we propose a new method for shear wave elastography(SWE)based on laser-excited shear wave,called photoacoustic shear wave elastography(PASWE),which combines photoacoustic(PA)technology with ultrafast ultrasound imaging.By using a focused laser to excite shear waves and ultrafast ultrasonic imaging for detection,high-frequency excitation of shear waves and noncontact elastic imaging can be realized.The laser can stimulate the tissue with the light absorption characteristic to produce the thermal expansion,thus producing the shear wave.The frequency of shear wave induced by laser is higher and the frequency band is wider.By tracking the propagation of shear wave,Young’s modulus of tissue is reconstructed in the whole shear wave propagation region to further evaluate the elastic information of tissue.The feasibility of the method is verified by experiments.Compared with the experimental results of supersonic shear imaging(SSI),it is proved that the method can be used for quantitative elastic imaging of the phantoms.In addition,compared with the SSI method,this method can realize the noncontact excitation of the shear wave,and the frequency of the shear wave excited by the laser is higher than that of the acoustic radiation force(ARF),so the spatial resolution is higher.Compared to the traditional PA elastic imaging method,this method can obtain a larger imaging depth under the premise of ensuring the imaging resolution,and it has potential application value in the clinical diagnosis of diseases requiring noncontact quantitative elasticity.

Retraction to: Shear strength of extremely altered serpentinites based on degree of saturation (Ankara, Turkey)

The Editor in Chief has retracted this article.After publication,concerns were raised over potentially inconsistent data-for example,the author mentioned residual soil with silt fraction,but Fig.5 shows even the smallest grains are beyond silt grains.Further analysis showed that Fig.6 and Fig.9 show both loose and dense soils exhibiting very similar dilatance behaviour,which is not typical.The author shared raw data,but it was incomplete and could not be verified.The Editor in Chief,therefore,has lost confidence in the findings of this article.The author does not agree to this retraction.

Numerical investigation of friction-heating-pressurization and its control parameters in the shear band of high-speed landslides

High-speed sliding often leads to catastrophic landslides,many of which,in the initial sliding phase before disintegration,experience a friction-induced thermal pressurization effect in the bottom shear band,accelerating the movement of the overlying sliding mass.To quantitatively investigate this complex multiphysical phenomenon,we established a set of equations that describe the variations in temperature and excess pore pressure within the shear band,as well as the conservation of momentum equation for the overlying sliding mass.With a simplified landslide model,we investigated the variations of temperature and excess pore pressure within the shear band and their impacts on the velocity of the overlying sliding mass.On this basis,we studied the impact of seven key parameters on the maximum temperature and excess pore pressure in the shear band,as well as the impact on the velocity of the overlying sliding mass.The simulation results of the standard model show that the temperature and excess pore pressure in the shear band are significantly higher than those in the adjacent areas,and reach the maximum values in the center.Within a few seconds after the start,the maximum excess pore pressure in the shear zone is close to the initial stress,and the shear strength loss rate exceeds 90%.The thermal pressurization mechanism significantly increases the velocity of the overlying sliding mass.The results of parameter sensitivity analysis show that the thermal expansion coefficient has the most significant impact on the temperature and excess pore pressure in the shear band,and the sliding surface dip angle has the most significant impact on the velocity of the overlying sliding mass.The results of this study are of great significance for clarifying the mechanism of thermal pressurization-induced high-speed sliding.

Roughness characterization and shearing dislocation failure for rock-backfill interface

Shearing dislocation is a common failure type for rock–backfill interfaces because of backfill sedimentation and rock strata movement in backfill mining goaf.This paper designed a test method for rock–backfill shearing dislocation.Using digital image techno-logy and three-dimensional(3D)laser morphology scanning techniques,a set of 3D models with rough joint surfaces was established.Further,the mechanical behavior of rock–backfill shearing dislocation was investigated using a direct shear test.The effects of interface roughness on the shear–displacement curve and failure characteristics of rock–backfill specimens were considered.The 3D fractal dimen-sion,profile line joint roughness coefficient(JRC),profile line two-dimensional fractal dimension,and the surface curvature of the frac-tures were obtained.The correlation characterization of surface roughness was then analyzed,and the shear strength could be measured and calculated using JRC.The results showed the following:there were three failure threshold value points in rock–backfill shearing dis-location:30%–50%displacement before the peak,70%–90%displacement before the peak,and 100%displacement before the peak to post-peak,which could be a sign for rock–backfill shearing dislocation failure.The surface JRC could be used to judge the rock–backfill shearing dislocation failure,including post-peak sliding,uniform variations,and gradient change,corresponding to rock–backfill disloca-tion failure on the field site.The research reveals the damage mechanism for rock–backfill complexes based on the free joint surface,fills the gap of existing shearing theoretical systems for isomerism complexes,and provides a theoretical basis for the prevention and control of possible disasters in backfill mining.

Shear behavior and off-fault damage of saw-cut smooth and tension-induced rough joints in granite

The damage of rock joints or fractures upon shear includes the surface damage occurring at the contact asperities and the damage beneath the shear surface within the host rock.The latter is commonly known as off-fault damage and has been much less investigated than the surface damage.The main contribution of this study is to compare the results of direct shear tests conducted on saw-cut planar joints and tension-induced rough granite joints under normal stresses ranging from 1 MPa to 50 MPa.The shear-induced off-fault damages are quantified and compared with the optical microscope observation.Our results clearly show that the planar joints slip stably under all the normal stresses except under 50 MPa,where some local fractures and regular stick-slip occur towards the end of the test.Both post-peak stress drop and stick-slip occur for all the rough joints.The residual shear strength envelopes for the rough joints and the peak shear strength envelope for the planar joints almost overlap.The root mean square(RMS)of asperity height for the rough joints decreases while it increases for the planar joint after shear,and a larger normal stress usually leads to a more significant decrease or increase in RMS.Besides,the extent of off-fault damage(or damage zone)increases with normal stress for both planar and rough joints,and it is restricted to a very thin layer with limited micro-cracks beneath the planar joint surface.In comparison,the thickness of the damage zone for the rough joints is about an order of magnitude larger than that of the planar joints,and the coalesced micro-cracks are generally inclined to the shear direction with acute angles.The findings obtained in this study contribute to a better understanding on the frictional behavior and damage characteristics of rock joints or fractures with different roughness.

Study of hydro-mechanical behaviours of rough rock fracture with shear dilatancy and asperities using shear-flow model

The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.

Numerical simulation on the influence of plant root morphology on shear strength in the sandy soil,Northwest China

Serious riverbank erosion,caused by scouring and soil siltation on the bank slope in the lower reaches of the Tarim River,Northwest China urgently requires a solution.Plant roots play an important role in enhancing soil shear strength on the slopes to maintain slope soils,but the extent of enhancement of soil shear strength by different root distribution patterns is unclear.The study used a combination of indoor experiments and numerical simulation to investigate the effects of varying plant root morphologies on the shear strength of the sandy soil in the Tarim River.The results showed that:(1)by counting the root morphology of dominant vegetation on the bank slope,we summarized the root morphology of dominant vegetation along the coast as vertical,horizontal,and claw type;(2)the shear strength of root-soil composites(RSCs)was significantly higher than that of remolded soil,and the presence of root system made the strain-softening of soil body significantly weakened so that RSCs had better mechanical properties;and(3)compared with the lateral roots,the average particle contact degree of vertical root system was higher,and the transition zone of shear strength was more prominent.Hence,vegetation with vertical root system had the best effect on soil protection and slope fixation.The results of this study have important guiding significance for prevention and control of soil erosion in the Tarim River basin,the restoration of riparian ecosystems,and the planning of water conservancy projects.