A macro-mesoscopic constitutive model for porous and cracked rock under true triaxial conditions

The complex mechanical and damage mechanisms of rocks are intricately tied to their diverse mineral compositions and the formation of pores and cracks under external loads.Numerous rock tests reveal a complex interplay between the closure of porous defects and the propagation of induced cracks,presenting challenges in accurately representing their mechanical properties,especially under true triaxial stress conditions.This paper proposes a conceptualization of rock at the mesoscopic level as a two-phase composite,consisting of a bonded medium matrix and frictional medium inclusions.The bonded medium is characterized as a mesoscopic elastic material,encompassing various minerals surrounding porous defects.Its mechanical properties are determined using the mixed multi-inclusion method.Transformation of the bonded medium into the frictional medium occurs through crack extension,with its elastoplastic properties defined by the DruckerePrager yield criterion,accounting for hardening,softening,and extension.MorieTanaka and Eshelby’s equivalent inclusion methods are applied to the bonded and frictional media,respectively.The macroscopic mechanical properties of the rock are derived from these mesoscopic media.Consequently,a True Triaxial Macro-Mesoscopic(TTMM)constitutive model is developed.This model effectively captures the competitive effect and accurately describes the stress-deformation characteristics of granite.Utilizing the TTMM model,the strains resulting from porous defect closure and induced crack extension are differentiated,enabling quantitative determination of the associated damage evolution.

Application of modified discontinuous deformation analysis to dynamic modelling of the Baige landslide in the Jinsha River

Accurate dynamic modeling of landslides could help understand the movement mechanisms and guide disaster mitigation and prevention.Discontinuous deformation analysis(DDA)is an effective approach for investigating landslides.However,DDA fails to accurately capture the degradation in shear strength of rock joints commonly observed in high-speed landslides.In this study,DDA is modified by incorporating simplified joint shear strength degradation.Based on the modified DDA,the kinematics of the Baige landslide that occurred along the Jinsha River in China on 10 October 2018 are reproduced.The violent starting velocity of the landslide is considered explicitly.Three cases with different violent starting velocities are investigated to show their effect on the landslide movement process.Subsequently,the landslide movement process and the final accumulation characteristics are analyzed from multiple perspectives.The results show that the violent starting velocity affects the landslide motion characteristics,which is found to be about 4 m/s in the Baige landslide.The movement process of the Baige landslide involves four stages:initiation,high-speed sliding,impact-climbing,low-speed motion and accumulation.The accumulation states of sliding masses in different zones are different,which essentially corresponds to reality.The research results suggest that the modified DDA is applicable to similar high-level rock landslides.

Global stability coefficient of large underground caverns under static loading and earthquake wave condition

Underground energy and resource development,deep underground energy storage and other projects involve the global stability of multiple interconnected cavern groups under internal and external dynamic disturbances.An evaluation method of the global stability coefficient of underground caverns based on static overload and dynamic overload was proposed.Firstly,the global failure criterion for caverns was defined based on its band connection of plastic-strain between multi-caverns.Then,overloading calculation of the boundary geostress and seismic intensity on the caverns model was carried out,and the critical unstable state of multi-caverns can be identified,if the plastic-strain band appeared between caverns during these overloading processes.Thus,the global stability coefficient for the multi-caverns under static loading and earthquake was obtained based on the corresponding overloading coefficient.Practical analysis for the Yingliangbao(YLB)hydraulic caverns indicated that this method can not only effectively obtain the global stability coefficient of caverns under static and dynamic earthquake conditions,but also identify the caverns’high-risk zone of local instability through localized plastic strain of surrounding rock.This study can provide some reference for the layout design and seismic optimization of underground cavern group.

深部原位应力环境钻进过程中岩心饼化成因机制:数值模拟与室内试验

深部岩体原位应力赋存状态一直是“黑箱”问题,基于岩心的应力测试技术已成为有效的解决手段之一。在深部原位应力环境钻进过程中时常出现岩心饼化现象,为了探索其在砂岩中的形成机制,利用PFC2D研究了不同原位应力条件下钻进过程中的裂纹分布和能量演化,获得了易发生岩心饼化的特定原位应力条件。应用自主开发的测试系统验证了实验室环境内岩心饼化所需的应力条件,并分析了岩心破裂特征、断面形貌与原位应力之间的关系。结果表明,原位应力越高,在钻进过程中岩石产生的拉裂纹越多,尤其是在孔壁和岩心根部,而岩心处裂纹从外表面向内部发展。更高的原位应力水平同样会导致更强的能量转换,进而引发岩石破裂。岩心饼化的形成需满足特定的应力条件。当径向应力(σ_(r))为最大主应力且保持45 MPa恒定、轴向应力(σa)为25 MPa或30 MPa时,岩心根部均会发生饼化现象。随着σ_(r)和σa之间差值或钻进深度的增加,岩心厚度减小,导致断面更光滑、分形维数更小,从而使岩心饼化现象更加显著。该研究可为科学阐明深部原位应力条件下岩心饼化的形成机制提供必要的技术和数据支持。

Experimental Study on Multistage Seismic Damage Process of Bedding Rock Slope:A Case Study of the Xinmo Landslide

In the early hours of June 24,2017,a major landslide event occurred in Xinmo Village,Sichuan Province,China.The landslide instantly devastated the whole village.Ten people died and 73 were missing in this major landslide event.The study area has suffered from several strong earthquakes in the past 100 y.Present studies have reported that the cumulative damage effect of the Xinmo landslide induced by earthquake is obvious.In this study,we conducted a shaking table test based on the detailed geological survey,historical seismic data,satellite optical image,unmanned aerial vehicle photography.The test result presents the characteristics of multistage seismic damage and progressive deformation process of the Xinmo landslide model,and shows that the historical earthquakes have caused serious damage to the interior of rock mass in the source area.The test also shows that the cumulative damage of the model increases with an increase in duration of earthquake loading.When the excitation intensity increases to a certain value,the damage accumulation velocity of the model suddenly increases.It reveals that frequent historical earthquake loads can be regarded as a main reason for the damage and deterioration of landslide rock mass.Damage accumulation and superposition occur in the slope.Under a long-term gravity,deformation of the slope gradually increases until catastrophic failure is triggered.The progressive deformation process of slope is summarized.Firstly,under strong earthquakes loading,a tensile fracture surface forms at the rear edge of the wavy deformation high and steep bedding slope.It reaches a certain critical depth and expands along the interlayer structural plane.Meantime,damaged fissures perpendicular to the structural plane also appear in the steep-gentle turning area of the slope.Secondly,under a coupling action of seismic loading and gravity,the interlaminar tensile crack surface at the rear edge of the slope extends to depth continuously.Meanwhile,rock fracture occurs in the steep-gentle turning area.The“two-way damage propagation”mode of the interlayer tensile crack surface occurs until the sliding surface is connected.However,due to the“locking section”effect of rock mass at the slope foot,it can still maintain a short-term stability.Thirdly,under the influences of the heavy rainfall before a landslide and the long-term gravity of the upper sliding mass,rock mass in the steep section at the slope foot breaks outward.Finally,a catastrophic landslide occurs.

Blast damage zone strength reduction method for deep cavern excavation and its application

The drill and blast(D&B)method is widely used to excavate underground spaces,but explosions generally cause damage to the rock.Still,no blast simulation method can provide computational accuracy and efficiency.In this paper,a blast equivalent simulation method called the blast damage zone strength reduction(BDZSR)method is proposed.This method first calculates the range of the blast-induced damage zone(BDZ)by formulae,then reduces the strength and deformation parameters of the rock within the BDZ ahead of excavation,and finally calculates the excavation damage zone(EDZ)for the D&B method by numerical simulation.This method combines stress wave attenuation,rock damage criteria and stress path variation to derive the BDZ depth calculation formulae.The formulae consider the initial geo-stress,and the reliability is verified by numerical simulations.The calculation of BDZ depth with these formulae allows the corresponding numerical simulation to avoid the time-consuming dynamic calculation process,thus greatly enhancing the calculation efficiency.The method was applied to the excavation in Jinping Class II hydropower station to verify its feasibility.The results show that the BDZSR method can be applied to blast simulation of underground caverns and provide a new way to study blast-induced damage.

The damage evolution behavior of polypropylene fiber reinforced concrete subjected to sulfate attack based on acoustic emission

To study the damage evolution behavior of polypropylene fiber reinforced concrete(PFRC)subjected to sulfate attack,a uniaxial compression test was carried out based on acoustic emission(AE).The effect of sulfate attack relative to time and fiber hybridization were analyzed and the compression damage factor was calculated using a mathematical model.The changes to AE ringing counts during the compression could be divided into compaction,elastic,and AE signal hyperactivity stages.In the initial stage of sulfate attack,the concrete micropores and microcracks were compacted gradually under external load and a corrosion products filling effect,and this corresponded with detection of few AE signals and with concrete compression strength enhancement.With increasing sulfate attack time,AE activity decreased.The cumulative AE ringing counts of PFRC at all corrosion ages were much higher than those for plain concrete.PFRC could still produce AE signals after peak load due to drawing effect of polypropylene fiber.After 150 d of sulfate attack,the cumulative AE ringing counts of plain concrete went down by about an order of magnitude,while that for PFRC remained at a high level.The initial damage factor of hybrid PFRC was-0.042 and-0.056 respectively after 150 d of corrosion,indicating that the advantage of hybrid polypropylene fiber was more obvious than plain concrete and single-doped PFRC.Based on a deterioration equation,the corrosion resistance coefficient of hybrid PFRC would be less than 0.75 after 42 drying-wetting sulfate attack cycles,which was 40%longer than that of plain concrete.

Optimization of Solid-Phase Extraction Conditions for Perfluorooctanoic Acid in Leachate

In order to optimize the solid phase extraction(SPE)conditions of perfluorooctanoic acid(PFOA)in the raw leachate and treated leachate,the effects of activator properties,SPE cartridge,p H value,ionic strength,and eluent properties were studied through single factor experiments.The optimal results of each single factor were obtained.Considering that the concentration of PFOA in the treated leachate is lower than that of the raw leachate,the SPE conditions of the treated leachate have been further optimized.Based on the above single-factor experiment,the main influencing factors were screened out as the volume of activator,ionic strength,and volume of eluent,and the three-factor threelevel response surface methodology(RSM)was optimized.The optimum SPE conditions of PFOA from treated landfill leachate were as follows:Activation of weak anion exchange(WAX)cartridge with 10 mL methanol,dosage of 600 mg KCl,6 m L 1%ammonia methanol eluted PFOA,the theory recovery rate is over95.67%.It has been verified that the error between the predicted value and the actual extraction recovery is small and has good repeatability.