Experimental study on the effect of water absorption level on rockburst occurrence of sandstone

To investigate the mechanism of rockburst prevention by spraying water onto the surrounding rocks,15 experiments are performed considering different water absorption levels on a single face.High-speed photography and acoustic emission(AE)system are used to monitor the rockburst process.The effect of water on sandstone rockburst and the prevention mechanism of water on sandstone rockburst are analyzed from the perspective of energy and failure mode.The results show that the higher the ab-sorption degree,the lower the intensity of the rockburst after absorbing water on single side of sand-stone.This is reflected in the fact that with the increase in the water absorption level,the ejection velocity of rockburst fragments is smaller,the depth of the rockburst pit is shallower,and the AE energy is smaller.Under the water absorption level of 100%,the magnitude of rockburst intensity changes from medium to slight.The prevention mechanism of water on sandstone rockburst is that water reduces the capacity of sandstone to store strain energy and accelerates the expansion of shear cracks,which is not conducive to the occurrence of plate cracking before rockburst,and destroys the conditions for rockburst incubation.

基于氮气吸附和压汞定量表征砂岩吸水孔隙结构演化规律

为了探究水对岩石孔隙结构的作用机制,本文依托山东万福煤矿砂岩地层,针对不同浸泡时间作用下的砂岩试样,开展扫描电镜、氮气吸脱附和高压压汞联合孔隙结构测试实验,分析微孔(<2 nm)、介孔(2~50 nm)和大孔(>50 nm)的孔隙结构变化规律。结果表明:砂岩孔隙结构类型为狭窄型和楔形结构,孔隙之间具有连通成缝的能力。砂岩吸水过程中孔隙度先减少后增加。微孔占比先减少后增加,介孔和大孔占比先增加后减少。根据广义分形维数和谱宽定量评估孔隙结构,得出其与孔隙质量正相关,即分维差值和谱宽越大,孔隙结构越复杂;而与岩石孔隙度呈负相关,即较大孔隙度对应较小的分维数差值和谱宽。

Intermittent disturbance mechanical behavior and fractional deterioration mechanical model of rock under complex true triaxial stress paths

Mechanical excavation,blasting,adjacent rockburst and fracture slip that occur during mining excavation impose dynamic loads on the rock mass,leading to further fracture of damaged surrounding rock in three-dimensional high-stress and even causing disasters.Therefore,a novel complex true triaxial static-dynamic combined loading method reflecting underground excavation damage and then frequent intermittent disturbance failure is proposed.True triaxial static compression and intermittent disturbance tests are carried out on monzogabbro.The effects of intermediate principal stress and amplitude on the strength characteristics,deformation characteristics,failure characteristics,and precursors of monzogabbro are analyzed,intermediate principal stress and amplitude increase monzogabbro strength and tensile fracture mechanism.Rapid increases in microseismic parameters during rock loading can be precursors for intermittent rock disturbance.Based on the experimental result,the new damage fractional elements and method with considering crack initiation stress and crack unstable stress as initiation and acceleration condition of intermittent disturbance irreversible deformation are proposed.A novel three-dimensional disturbance fractional deterioration model considering the intermediate principal stress effect and intermittent disturbance damage effect is established,and the model predicted results align well with the experimental results.The sensitivity of stress states and model parameters is further explored,and the intermittent disturbance behaviors at different f are predicted.This study provides valuable theoretical bases for the stability analysis of deep mining engineering under dynamic loads.

Analysis of the interaction between bolt-reinforced rock and surface support in tunnels based on convergence-confinement method

To investigate the interaction of the bolt-reinforced rock and the surface support,an analytical model of the convergence-confinement type is proposed,considering the sequential installation of the fully grouted rockbolts and the surface support.The rock mass is assumed to be elastic-brittle-plastic material,obeying the linear Mohr-Coulomb criterion or the non-linear Hoek-Brown criterion.According to the strain states of the tunnel wall at bolt and surface support installation and the relative magnitude between the bolt length and the plastic depth during the whole process,six cases are categorized upon solving the problem.Each case is divided into three stages due to the different effects of the active rockbolts and the passive surface support.The fictitious pressure is introduced to quantify the threedimensional(3D)effect of the tunnel face,and thus,the actual physical location along the tunnel axis of the analytical section can be considered.By using the bolt-rock strain compatibility and the rocksurface support displacement compatibility conditions,the solutions of longitudinal tunnel displacement and the reaction pressure of surface support along the tunnel axis are obtained.The proposed analytical solutions are validated by a series of 3D numerical simulations.Extensive parametric studies are conducted to examine the effect of the typical parameters of rockbolts and surface support on the tunnel displacement and the reaction pressure of the surface support under different rock conditions.The results show that the rockbolts are more effective in controlling the tunnel displacement than the surface support,which should be installed as soon as possible with a suitable length.For tunnels excavated in weak rocks or with restricted displacement control requirements,the surface support should also be installed or closed timely with a certain stiffness.The proposed method provides a convenient alternative approach for the optimization of rockbolts and surface support at the preliminary stage of tunnel design.

基于大变形升阻锚索深埋巷道围岩稳定性控制研究

研发了一种新型大变形升阻锚索支护结构,通过串联6个壁厚5 mm含环形V槽的直管型钢结构实现锚索吸能让压,吸能率为44.55 J/g。拉拔试验结果表明新型锚索极限变形量(221.58 mm)与普通传统锚索相比增加了4.87倍,具有大变形、增阻及高吸能等优良性能。通过离散元UDEC程序的二次开发,模拟了新型大变形增阻锚索支护阻力−变形曲线的双弹性阶段。采用新型锚索对深埋巷道围岩支护效果数值模拟结果表明,与无支护及传统锚索支护相比,顶板下沉量分别减小了88.33%与74.64%,稳定状态下最大主应力峰值分别减小了23.15%与18.34%,塑性区厚度分别减小了37.50%与62.50%。最后将新型大变形增阻锚索成功应用于真实巷道,两帮最大收敛变形量约为300 mm,实现了较好的支护效果。

An improved strain-softening constitutive model of granite considering the effect of crack deformation

This paper presents an improved strain-softening constitutive model considering the effect of crack deformation based on the triaxial cyclic loading and unloading test results.The improved model assumes that total strain is a combination of plastic,elastic,and crack strains.The constitutive relationship between the crack strain and the stress was further derived.The evolutions of mechanical parameters,i.e.strength parameters,dilation angle,unloading elastic modulus,and deformation parameters of crack,with the plastic strain and confining pressure were studied.With the increase in plastic strain,the cohesion,friction angle,dilation angle,and crack Poisson's ratio initially increase and subsequently decrease,and the unloading elastic modulus and the crack elastic modulus nonlinearly decrease.The increasing confining pressure enhances the strength and unloading elastic modulus,and decreases the dilation angle and Poisson's ratio of the crack.The theoretical triaxial compressive stress-strain curves were compared with the experimental results,and they present a good agreement with each other.The improved constitutive model can well reflect the nonlinear mechanical behavior of granite.

基于切削能密度的岩体强度-裂隙随钻识别方法

岩体强度与裂隙参数是反映工程岩体质量的基本参数,其准确获取是进行地下工程围岩分级与支护合理设计优化的前提。传统的岩体强度测试方法需要对现场围岩取芯并运至实验室测试,测试结果难以反映工程现场环境下的岩体力学性质,而对于岩体强度与裂隙参数的原位测试方法研究较少。本文基于能量守恒定律,建立了岩石切削能密度与随钻参数的关系式,构建了岩体等效强度随钻反演模型(ES-DP模型),系统开展了岩体数字钻进试验。结果表明,相对于传统测试方法,岩体等效强度随钻测试结果平均差异率为2.4%,验证了ES-DP模型对岩体等效强度测试的有效性。在此基础上,建立了岩体裂隙参数随钻识别模型,与实际测量结果相比,该模型测得的岩体裂隙位置与宽度平均精度分别为1.8和1.6 mm,测试精度较高。基于上述研究,本文提出了一种岩体强度-裂隙随钻识别方法,为实现地下工程围岩强度参数与结构特征的原位实时测试提供了新方法。

Velocity analysis of supersonic jet flow in double-cone ignition scheme

In the double-cone ignition schemes(DCIS), the deuterium–tritium target shell is ablated and compressed by a highpower nanosecond laser in Au-cones to generate plasmas. Under the actions of spherically symmetric compression and acceleration along the Au cone, they will be ejected out of the cone mouth and collide with each other. The plasmas experience conversion from kinetic energy to internal energy at the vertex of the geometric center of two Au cones that are symmetric to each other, because of which high-density fusion plasmas are preheated. This key physical process has undergone experimental verification on the Shenguang-II upgraded facility in China. Apparently, the improvement and optimization of the velocity of plasmas in hypersonic jet flow at the cone mouth are crucial for the success of the DCIS. In the DCIR7 experiment of the Shenguang-II upgraded facility, a velocity yield of approximately 130–260 km/s was achieved for the plasmas at the cone mouth, with a result of nearly 300 km/s based on numerical simulation. In this paper, theoretical analysis is performed as regards the process, in which target shells are ablated and compressed by laser to generate highvelocity plasmas ejected through jet flow. Based on this analysis, the formula for the velocity of plasmas in supersonic jet flow at the cone mouth is proposed. This study also provides measures that are more effective for improving the kinetic energy of plasmas and optimizing energy conversion efficiency, which can serve as theoretical references for the adjustment and optimization of processes in subsequent experiments.

Effect of surfactant frequently used in soil flushing on oxygen mass transfer in micro-nano-bubble aeration system

In-site soil flushing and aeration are the typical synergetic remediation technology for contaminated sites.The surfactant present in flushing solutions is bound to affect the aeration efficiency.The purpose of this study is to evaluate the effect of surfactant frequently used in soil flushing on the oxygen mass transfer in micro-nano-bubble(MNB)aeration system.Firstly,bio-surfactants and chemical surfactants were used to investigate their effects on Sauter mean diameter of bubble(dBS),gas holdup(ε),volumetric mass-transfer coefficient(kLa)and liquid-side mass-transfer coefficient(kL)in the MNB aeration system.Then,based upon the experimental results,the Sardeing's and Frossling's models were modified to describe the effect of surfactant on kL in the MNB aeration.The results showed that,for the twenty aqueous surfactant solutions,with the increase in surfactant concentration,the value of dBS,kLa and kL decreased,while the value ofεand gas-liquid interfacial area(a)increased.These phenomena were mainly attributed to the synergistic effects of immobile bubble surface and the suppression of coalescence in the surfactant solutions.In addition,with the presence of electric charge,MNBs in anionic surfactant solutions were smaller and higher in number than in non-ionic surfactant solutions.Furthermore,the accumulation of surfactant on the gas-liquid interface was more conspicuous for small MNB,so the reduction of kL in anionic surfactant solutions was larger than that in non-ionic surfactant solutions.Besides,the modified Frossling's model predicted the effect of surfactant on kL in MNB aeration system with reasonable accuracy.

Mechanical behaviour of fiber-reinforced grout in rock bolt reinforcement

Grouted rock bolts subject to axial loading in the field exhibit various failure modes,among which the most predominant one is the bolt-grout interface failure.Thus,mechanical characterization of the grout is essential for understanding its performance in ground support.To date,few studies have been conducted to characterize the mechanical behaviour of fiber-reinforced grout(FRG)in rock bolt reinforcement.Here we experimentally studied the mechanical behaviour of FRG under uniaxial compression,indirect tension,and direct shear loading conditions.We also conducted a series of pullout tests of rebar bolt encapsulated with different grouts including conventional cementitious grout and FRG.FRG was developed using 15%silica fume(SF)replacement of cement(by weight)and steel fiber to achieve highstrength and crack-resistance to overcome drawbacks of the conventional grout.Two types of steel fibers including straight and wavy steel fibers were further added to enhance the grout quality.The effect of fiber shape and fiber volume proportion on the grout mechanical properties were examined.Our experimental results showed that the addition of SF and steel fiber by 1.5%fiber volume proportion could lead to the highest compressive,tensile,and shear strengths of the grout.The minimum volume of fiber that could improve the mechanical properties of grout was found at 0.5%.The scanning electron microscopy(SEM)analysis demonstrated that steel fibers act as an excellent bridge to prevent the cracks from propagating at the interfacial region and hence to aid in maintaining the integrity of the cementitious grout.Our laboratory pullout tests further confirmed that FRG could prevent the cylindrical grout annulus from radial crack and hence improve the rebar’s load carrying capacity.Therefore,FRG has a potential to be utilized in civil and mining applications where high-strength and crack-resistance support is required.

二氧化碳静态爆破的边坡振动场地效应和能量分析

进行真三轴CO_(2)静态爆破压裂实验和现场测试,获得压力曲线、试样裂隙分布、现场振动曲线和岩体破坏特征,并通过分析压力曲线和裂缝扩展特征来研究致裂原理,同时还利用窗函数傅里叶变换提取振动特征,以评估破岩现场边坡的结构安全。研究结果表明:(1)CO_(2)静态爆破的峰值振动速度在8 m处衰减至满足建筑物保护要求,比传统的钻爆法更安全;(2)CO_(2)静态爆破引起的振动波形的主频高于钻爆法的,且随着距离的加长,振动能量的分布逐渐分散;(3)钻爆法的振动速度预测公式适用于CO_(2)静态爆破的,且带高程的预测公式更为准确。根据压力曲线、试样裂隙分布和现场岩体损伤,提出CO_(2)静态爆破的HIG破岩模型,为新技术的研究提供了基础框架。在实验中没有采集到温度数据,今后迫切需要对聚能剂的瞬时放热特性和温度对CO_(2)状态的影响进行研究。

A modified smoothed particle hydrodynamics method considering residual stress for simulating failure and its application in layered rock mass

Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.

Instability mechanism of mining roadway passing through fault at different angles in kilometre-deep mine and control measures of roof cutting and NPR cables

The angle α between the fault strike and the axial direction of the roadway produces different damage characteristics. In this paper, the research methodology includes theoretical analyses, numerical simulations and field experiments in the context of the Daqiang coal mine located in Shenyang, China. The stability control countermeasure of "pre-splitting cutting roof + NPR anchor cable"(PSCR-NPR) is simultaneously proposed. According to the different deformation characteristics of the roadway, the faults are innovatively classified into three types, with α of type I being 0°-30°, α of type II being 30°-60°, and α of type III being 60°-90°. The full-cycle stress evolution paths during mining roadway traverses across different types of faults are investigated by numerical simulation. Different pinch angles α lead to high stress concentration areas at different locations in the surrounding rock. The non-uniform stress field formed in the shallow surrounding rock is an important reason for the instability of the roadway. The pre-cracked cut top shifted the high stress region to the deep rock mass and formed a low stress region in the shallow rock mass. The high prestressing NPR anchor cable transforms the non-uniform stress field of the shallow surrounding rock into a uniform stress field. PSCR-NPR is applied in the fault-through roadway of Daqiang mine. The low stress area of the surrounding rock was enlarged by 3-7 times, and the cumulative convergence was reduced by 45%-50%. It provides a reference for the stability control of the deep fault-through mining roadway.

Dynamic Response of Foundations during Startup of High-Frequency Tunnel Equipment

The specialized equipment utilized in long-line tunnel engineering is evolving towards large-scale,multifunctional,and complex orientations.The vibration caused by the high-frequency units during regular operation is supported by the foundation of the units,and the magnitude of vibration and the operating frequency fluctuate in different engineering contexts,leading to variations in the dynamic response of the foundation.The high-frequency units yield significantly diverse outcomes under different startup conditions and times,resulting in failure to meet operational requirements,influencing the normal function of the tunnel,and causing harm to the foundation structure,personnel,and property in severe cases.This article formulates a finite element numerical computation model for solid elements using three-dimensional elastic body theory and integrates field measurements to substantiate and ascertain the crucial parameter configurations of the finite element model.By proposing a comprehensive startup timing function for high-frequency dynamic machines under different startup conditions,simulating the frequency andmagnitude variations during the startup process,and suggesting functions for changes in frequency and magnitude,a simulated startup schedule function for high-frequency machines is created through coupling.Taking into account the selection of the transient dynamic analysis step length,the dynamic response results for the lower dynamic foundation during its fundamental frequency crossing process are obtained.The validation checks if the structural magnitude surpasses the safety threshold during the critical phase of unit startup traversing the structural resonance region.The design recommendations for high-frequency units’dynamic foundations are provided,taking into account the startup process of the machine and ensuring the safe operation of the tunnel.

Parameters Optimization and Performance Evaluation of the Tuned Inerter Damper for the Seismic Protection of Adjacent Building Structures

In order to improve the seismic performance of adjacent buildings,two types of tuned inerter damper(TID)damping systems for adjacent buildings are proposed,which are composed of springs,inerter devices and dampers in serial or in parallel.The dynamic equations of TID adjacent building damping systems were derived,and the H2 norm criterion was used to optimize and adjust them,so that the system had the optimum damping performance under white noise random excitation.Taking TID frequency ratio and damping ratio as optimization parameters,the optimum analytical solutions of the displacement frequency response of the undamped structure under white noise excitation were obtained.The results showed that compared with the classic TMD,TID could obtain a better damping effect in the adjacent buildings.Comparing the TIDs composed of serial or parallel,it was found that the parallel TIDs had more significant advantages in controlling the peak displacement frequency response,while the H2 norm of the displacement frequency response of the damping system under the coupling of serial TID was smaller.Taking the adjacent building composed of two ten-story frame structures as an example,the displacement and energy collection time history analysis of the adjacent building coupled with the optimum design parameter TIDs were carried out.It was found that TID had a better damping effect in the full-time range compared with the classic TMD.This paper also studied the potential power of TID in adjacent buildings,which can be converted into available power resources during earthquakes.

Digital monitoring of rotary-percussive drilling with down-the-hole hammer for profiling weathered granitic ground

Rock and geotechnical engineering investigations involve drilling holes in ground with or without retrieving soil and rock samples to construct the subsurface ground profile.On the basis of an actual soil nailing drilling for a slope stability project in Hong Kong,this paper further develops the drilling process monitoring(DPM)method for digitally profiling the subsurface geomaterials of weathered granitic rocks using a compressed airflow driven percussive-rotary drilling machine with down-the-hole(DTH)hammer.Seven transducers are installed on the drilling machine and record the chuck displacement,DTH rotational speed,and five pressures from five compressed airflows in real-time series.The mechanism and operations of the drilling machine are elaborated in detail,which is essential for understanding and evaluating the drilling data.A MATLAB program is developed to automatically filter the recorded drilling data in time series and classify them into different drilling processes in sub-time series.These processes include penetration,push-in with or without rod,pull-back with or without rod,rod-tightening and rod-untightening.The drilling data are further reconstructed to plot the curve of drill-bit depth versus the net drilling time along each of the six drillholes.Each curve is found to contain multiple linear segments with a constant penetration rate,which implies a zone of homogenous geomaterial with different weathering grades.The effect from fluctuation of the applied pressures is evaluated quantitatively.Detailed analyses are presented for accurately assess and verify the underground profiling and strength in weathered granitic rock,which provided the basis of using DPM method to confidently assess drilling measurements to interpret the subsurface profile in real time.

Damage constitutive model of lunar soil simulant geopolymer under impact loading

Lunar base construction is a crucial component of the lunar exploration program,and considering the dynamic characteristics of lunar soil is important for moon construction.Therefore,investigating the dynamic properties of lunar soil by establishing a constitutive relationship is critical for providing a theoretical basis for its damage evolution.In this paper,a split Hopkinson pressure bar(SHPB)device was used to perform three sets of impact tests under different pressures on a lunar soil simulant geopolymer(LSSG)with sodium silicate(Na_(2)SiO_(3))contents of 1%,3%,5%and 7%.The dynamic stressestrain curves,failure modes,and energy variation rules of LSSG under different pressures were obtained.The equation was modified based on the ZWT viscoelastic constitutive model and was combined with the damage variable.The damage element obeys the Weibull distribution and the constitutive equation that can describe the mechanical properties of LSSG under dynamic loading was obtained.The results demonstrate that the dynamic compressive strength of LSSG has a marked strain-rate strengthening effect.Na_(2)SiO_(3) has both strengthening and deterioration effects on the dynamic compressive strength of LSSG.As Na_(2)SiO_(3) grows,the dynamic compressive strength of LSSG first increases and then decreases.At a fixed air pressure,5%Na_(2)SiO_(3) had the largest dynamic compressive strength,the largest incident energy,the smallest absorbed energy,and the lightest damage.The ZWT equation was modified according to the stress response properties of LSSG and the range of the SHPB strain rate to obtain the constitutive equation of the LSSG,and the model’s correctness was confirmed.

双线盾构隧道开挖引起既有铁路路基沉降及盾构掘进参数研究:工程案例

为揭示盾构掘进参数对既有铁路路基沉降的影响规律,本文依托北京地铁16号线丰益桥南站至丰台站双线盾构隧道下穿既有铁路线路工程进行案例分析。将盾构掘进过程分为盾构穿越前、盾构穿越中和盾构穿越后三个阶段。实测数据表明,路基沉降在盾构穿越前逐渐增大,在盾构穿越过程中急剧增大,在盾构穿越后趋于稳定。将盾构管片及时封闭成环,施作补偿注浆、二次注浆,有利于控制路基沉降。此外,盾构掘进参数对控制路基沉降具有至关重要的作用。研究结果表明:盾构掘进速度越快引起的路基沉降越大,继而引起较大的土仓压力,并伴随刀盘扭矩和盾构总推力的骤减;由于盾构过程中出土量易于控制,因此其对路基沉降影响不大;然而,如果同步注浆量控制不当,将导致既有路基向上隆起。本研究为盾构开挖下穿多条既有铁路线路引起的路基变形控制以及盾构掘进参数设置提供参考和指导。

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.

Method of fabricating artificial rock specimens based on extrusion free forming(EFF)3D printing

Three-dimensional(3D)printing technology has been widely used to create artificial rock samples in rock mechanics.While 3D printing can create complex fractures,the material still lacks sufficient similarity to natural rock.Extrusion free forming(EFF)is a 3D printing technique that uses clay as the printing material and cures the specimens through high-temperature sintering.In this study,we attempted to use the EFF technology to fabricate artificial rock specimens.The results show the physico-mechanical properties of the specimens are significantly affected by the sintering temperature,while the nozzle diameter and layer thickness also have a certain impact.The specimens are primarily composed of SiO_(2),with mineral compositions similar to that of natural rocks.The density,uniaxial compressive strength(UCS),elastic modulus,and tensile strength of the printed specimens fall in the range of 1.65–2.54 g/cm3,16.46–50.49 MPa,2.17–13.35 GPa,and 0.82–17.18 MPa,respectively.It is capable of simulating different types of rocks,especially mudstone,sandstone,limestone,and gneiss.However,the simulation of hard rocks with UCS exceeding 50 MPa still requires validation.