Numerical Simulation of the Rock SHPB Test with a Special Shape Striker Based on the Discrete Element Method

A split Hopkinson pressure bar (SHPB) system with a special shape striker has been suggested as the test method by the International Society for Rock Mechanics (ISRM) to determine the dynamic characteristics of rock materials. In order to further verify this testing technique and microscopically reveal the dynamic responses of specimens in SHPB tests, a numerical SHPB test system was established based on particle flow code (PFC). Numerical dynamic tests under different impact velocities were conducted. Investigation of the stresses at the ends of a specimen showed that the specimen could reach stress equilibrium after several wave reverberations, and this balance could be maintained well for a certain time period after the peak stress. In addition, analyses of the reflected waves showed that there was a clear relationship between the variation of the reflected wave and the stress equilibrium state in the specimen, and the turning point of the reflected wave corresponded well with the peak stress in the specimen. Furthermore, the reflected waves can be classified into three types according to their patterns. Under certain impact velocities, the specimen deforms at a constant strain rate during the whole loading process. Finally, the influence of the micro-strength ratio ( \({{\tau_{\text{c}} } \mathord{\left/ {\vphantom {{\tau_{\text{c}} } {\sigma_{\text{c}} }}} \right. \kern-0pt} {\sigma_{\text{c}} }}\) ) and distribution pattern on the dynamic increase factor (DIF) of the strength DIF were studied, and the lateral inertia confinement and heterogeneity were found to be two important factors causing the strain rate effect for rock materials.     

单轴压缩下非饱和黏土动态力学性能试验研究

黏土作为一种松散多孔材料,其波阻抗较低,采用传统分离式霍普金森压杆(SHPB)试验装置进行单轴压缩试验时,存在透射信号弱、应力不平衡等问题。针对以上问题,提出利用石英压电晶体测定试样的应力并监测试样两端的应力平衡,采用入射波整形技术,考察应变率、含水率对非饱和黏土动态力学性能的影响规律。试验结果表明,非饱和黏土试样的单轴抗压强度与应变率及含水率的关系密切,强度随着应变率的提升而上升,但随着含水率的提升而下降。通过对试验结果的分析,分别提出单轴抗压强度与含水率及应变率之间的关系式,并通过拟合确定参数,拟合结果较好。     

3D-DIC系统在岩石力学试验中的应用

介绍了全场3D-DIC系统,在试样周围均匀布置三组3D采集元,利用散斑图像的采集,尽可能在较大范围内捕捉试样表面的变形并进行分析。为了验证3D-DIC系统测量结果的可靠性,以江持安山岩为试验材料,进行了单轴压缩条件下的恒定应变速率和交替应变速率试验。通过3D-DIC系统的测量结果与应变片测量结果的对比发现,两者吻合度较好,表明3D-DIC系统能够满足岩石力学试验非接触式、可视化变形测量的需求。另外以田下凝灰岩为试验材料,通过图像分析得到了加载过程中岩石表面的全场变形信息,发现岩石变形的不均匀性和应变局部化现象,表明利用3D-DIC系统对岩石破裂过程变形测量比传统变形测量方法更具有优越性,为岩石力学试验研究提供了一种新的、便捷的测试方法。     

Numerical Investigation of the Dynamic Compressive Behaviour of Rock Materials at High Strain Rate

The dynamic compressive strength of rock materials increases with the strain rate. They are usually obtained by conducting laboratory tests such as split Hopkinson pressure bar (SHPB) test or drop-weight test. It is commonly agreed now that the dynamic increase factor (DIF) obtained from impact test is affected by lateral inertia confinement, friction confinement between the specimen and impact materials and the specimen sizes and geometries. Therefore, those derived directly from testing data do not necessarily reflect the true dynamic material properties. The influences of these parameters, however, are not straightforward to be quantified in laboratory tests. Therefore, the empirical DIF relations of rock materials obtained directly from impact tests consist of contributions from lateral inertia and end friction confinements, which need be eliminated to reflect the true dynamic material properties. Moreover, different rocks, such as granite, limestone and tuff have different material parameters, e.g., equation of state (EOS) and strength, which may also affect the DIF of materials but are not explicitly studied in the open literature. In the present study, numerical models of granite, limestone and tuff materials with different EOS and strength under impact loads are developed to simulate SHPB tests and to study the influences of EOS and strength, lateral inertia confinement and end friction confinement effects on their respective DIFs in the strain rate range between 1 and 1,000 s^?1. The commercial software AUTODYN with user-provided subroutines is used to perform the numerical simulations of SHPB tests. Numerical simulation results indicate that the lateral inertia confinement, friction confinement and specimen aspect (L/D) ratio significantly influence DIF obtained from impact tests and the inertia confinement effect is different for different rocks. Based on the numerical results, quantifications on the relative contributions from the lateral inertia confinement and the material strain rate effect on DIF of granite, limestone and tuff material compressive strength are made. The effects of friction coefficient, L/D ratio and rock type on DIF are discussed. Empirical relations of DIF with strain rate for the three rock materials representing the true material strain rate effect are also proposed.     

Interpretation of torsional shear results for nonlinear stress–strain relationship

Torsional tests are frequently used in material testing to determine shear stress–strain relationships. In the tests, torque is applied to a solid or hollow cylindrical specimen, whose twist angle is recorded as the response. While the raw test results are readily available in terms of torque–twist sets, converting them into intrinsic stress–strain relationships is not always straightforward, owing to the nonlinear stress–strain behavior and the non‐uniform radial distribution of the stress and strain fields. It follows that the approaches currently used to interpret torsional data are essentially empirical based. The present paper addresses this data interpretation issue from a more rigorous perspective. The paper includes two parts: the first part concerns hollow cylindrical specimens under torsion, where an analytical‐based recursive procedure is developed; and the second part deals with solid torsional specimens, for which an analytical solution given by Nadai is directly employed. Examples pertaining to these two cases are examined and the resulting shear stress–strain curves are compared with those based on a number of popular empirical methods. The proposed analytical approaches can serve as a benchmark for checking the accuracy and validity of those empirical propositions. It is worth noting that, for brittle materials exhibiting softening behavior, only the analytical method gives the correct interpretation. Copyright © 2007 John Wiley & Sons, Ltd.     

不同应力路径下土石混填体变形力学特性大型三轴试验研究

为更加全面地掌握土石混填体的变形、力学特性,采用大型三轴剪切仪分别在常规路径、等应力比路径和等p路径下对不同含石量的土石混填体进行了一系列试验研究。试验结果表明,在低围压下土石混填体的应变软化特征并不明显,其黏聚力普遍较低,而内摩擦角则较高,随试样的含石量从25%增加到70%,其内摩擦角和变形模量亦会随之增加。土石混填体在常规路径和等p路径下的应力-应变关系曲线比较相似,其所求得的抗剪强度指标亦比较接近;而等应力比路径下的应力-应变关系曲线的形态和发展过程均与之存在明显差别。土石混填体在常规路径下具有低压剪胀和高压剪缩的特性,含石量越高,低压剪胀性越明显,反之,则高压剪缩性更为显著;等应力比路径下试样的体积应变与轴向应变之间近似呈线性关系,其斜率随应力比增加而减小;等p路径下试样则是以剪胀性为主,尤其在平均主应力较低时更为显著。     

Experimental creep tests and prediction of long-term creep behavior of grouting material

Rock reinforcement systems, such as resin-grouted rock bolts, display complex creep behavior because both grouting materials and bolts show time-dependent behavior. In this paper, only the time-dependent behavior of grouting material was investigated, in which creep tests of grouting material was conducted in triaxial compression apparatus at room temperature. The test specimens were provided from the Araldite epoxy resin used in rock reinforcement. We attempt to predict long-term creep parameter using triaxial creep tests and to define time-dependent characteristics of the bounding material. In short-term creep tests, three different axial and confining stress levels were applied in steps to each specimen. The transient creep for all the stress levels were described by power function which fit properly to time–strain curves. The maximum difference between the proposed model and experimental long-term creep strain was less than 7.1 %. It was observed that the creep rate of a grouting material specimen directly depends on the deviator of stress (i.e., σ₁–σ₃).     

动载作用下端锚锚固体力学响应特征研究

为研究动载作用下巷道围岩失稳破坏机制,基于SHPB试验系统,研究了不同冲击气压（0.2、0.3、0.4、0.5、0.6 MPa）下预应力端锚锚固体的应力波传播规律及层裂破坏特征,采用图解法计算了锚固体试件的层裂强度及应变率,建立了锚固体层裂损伤模型并采用高速摄像机进行了验证。结果表明：不同冲击气压下试件应力波峰值应变均呈指数形式衰减,且空间衰减幅值及衰减指数与冲击气压呈正相关关系;锚固体试件层裂首先发生在自由端附近,沿反射波传播方向层裂厚度逐渐增加,0.3 MPa冲击气压以上,试件每两处层裂面中间产生一道新层裂;在18～33 s⁻¹应变率范围内,应变率率效应占主导地位,层裂强度最高达69 MPa,而应变率为41 s⁻¹时,层裂强度降低至17 MPa;锚固体自由端在动载冲击后出现层裂闭合现象,但在锚固端由于锚固界面残余黏聚力作用,层裂位置未发生闭合现象。     

Uniaxial compression and tension tests of anthracite and loading rate dependence of peak strength

The mechanical characteristics of coal under uniaxial compressive stress and uniaxial tensile stress are key factors in investigations of the stability of galleries and coal faces, and are vital considerations for the efficient design of coalmines, disaster prevention, and environmental preservation. These mechanical characteristics have long been a topic of research, but remain insufficiently understood. In this study we performed uniaxial compression tests and uniaxial tension tests on coal, with particular attention to two concerns. The first was to measure the loading rate dependence of the peak strength of coal. Coal shows high sample-to-sample scatter in strength; test methods that are successful in comparatively homogeneous rock types cannot be used for coal. We therefore used an alternative test method that was recently developed by the authors. Anthracite samples were subjected to alternating slow and fast strain rates. Measured variations in stress during this process were used to estimate the loading rate dependence of peak strength. The second objective was to obtain complete stress–strain curves for coal under uniaxial tensile stress. It is difficult to hold samples secure during a tensile test; consequently, such curves have yet to be obtained. This study presents a successful application of the authors' method to the analysis of coal samples, yielding a complete stress–strain curve under tensile stress. The two methods presented here hold promise for application not only to anthracite but also to a wide variety of coals. It is possible to derive the values of the constants used in constitutive equations from the obtained experimental results. Once these equations are determined, they can be incorporated into finite-element software to investigate various time-dependent behaviors of coal and aid in the efficient design of coalmines and the prevention of subsidence and mine disasters.     

基于应变空间硅藻质软岩的软化本构模型

三轴试验结果表明,软岩具有显著的应变软化特征。正常固结软岩表现为峰值后较缓的软化特征,而超固结软岩峰值强度后的应力-应变曲线呈陡降软化特征。应变空间表述的弹塑性理论在解决大应变和软化问题时比应力空间表述的弹塑性理论更具有优越性。基于应变空间的基本的弹塑性本构方程式,采用Mises剪切屈服准则及相关联流动法则,导出固结不排水三轴条件下的应力-应变本构关系式,并采用不同的硬化函数表达式对软岩在不同围压下的应力-应变曲线进行数值模拟。结果表明,应变空间的弹塑性理论能较好地模拟软岩应变软化特征,其中硬化函数的确定是关键问题之一。通过研究,提出了用于固结不排水状态下正常固结软岩的硬化函数形式。     

煤岩两体模型变形破坏数值模拟

采用拉格朗日元法,在弹性岩石与弹性-应变软化煤体所构成的平面应变两体模型的上、下端面上不存在水平方向摩擦力条件下,模拟了模型的破坏过程、岩石高度对模型及煤体全程应力-应变曲线、煤体变形速率、煤体破坏模式及剪切应变增量分布的影响。结果表明,当模型的全程应力-应变曲线达到峰值时煤体内部的剪切带图案已经十分明显,在模型的应变硬化阶段,煤体中的应变局部化可视为模型失稳破坏的前兆,随岩石高度的增加,模型应力-应变曲线的软化段变得陡峭,这与单轴压缩条件下的解析解在定性上是一致的;煤体应力-应变曲线的软化段变得平缓,煤体消耗能量的能力增强;弹性阶段煤体的变形速率降低;煤体内部的剪切应变增量增加。煤体应力-应变曲线的软化段的斜率、弹性阶段煤体的变形速率、煤体内部的剪切应变增量及塑性耗散能都受岩石高度的影响,说明了岩石几何尺寸对煤体的影响(煤岩相互作用)是不容忽视的。     

A fractional order creep constitutive model of warm frozen silt

A series of triaxial creep tests were conducted on warm frozen silts extracted from Qinghai–Tibet Plateau at temperature of ?1.5 °C under confining pressures of 0.5, 1.0, and 2.0 MPa, respectively. The applied test stress levels were 30, 50, 60, and 70% of triaxial shear strength, respectively. The test results indicate that the creep strain increases with the increase in applied stress level and there is a stress threshold, based on which the test results can be classified into two types of creep strain curves. The creep strain curve only includes primary and secondary creep stages when the stress level is less than the threshold value. When the stress level exceeds the threshold value, the creep strain velocity gradually increases and the specimen quickly fails in tertiary creep stage. Based on the creep test results, a fractional order rheological element model is established for warm frozen silt, which is also generalized from uniaxial stress state to the three-dimensional stress state. From the analysis on the features of the stress threshold, a creep strength criterion is also proposed simultaneously. Comparing the calculated results of the warm frozen silt with the tested ones, it is found that the predicted results of the proposed model are in good agreement with the test results. In the proposed fractional order model, the relationship between the damage factor and time is established to describe the damage degree of the specimen. Compared with the existing creep constitutive model of frozen soil, the proposed fractional order model has advantages of fewer model parameters, higher simulation precision and wider applicability in analyzing the mechanical properties of warm frozen silt.     

基于复合损伤的节理岩体动态本构模型研究

基于运用霍普金森压杆（SHPB）装置对节理岩体动载试验得出的数据,从节理面倾角、贯通度、厚度、组数、填充物及应变率等不同方面分析各因素对节理岩体力学特性的影响,通过对节理岩体在高应变率下的损伤机制和破坏形式进行分析的基础上,基于复合损伤理论,对广义Bingham模型进行改进,构造了节理岩体材料在不同应变率下动态响应的本构模型,对模型计算和试验结果的比较表明,该模型能够很好地描述节理岩体动荷载下初始弹性变形阶段、稳态塑性变形阶段和加速变形破坏阶段的应力-应变关系,并且理论与试验结果吻合较好,从而证明了该模型的正确性和合理性。     

盐岩蠕变特性和本构关系研究

对金坛盐岩蠕变规律进行了试验研究,分析了不同应力状态和温度对盐岩蠕变特性的影响,通过试验获取了蠕变曲线和蠕变特征参数。根据其蠕变曲线特征,用黏弹塑性理论建立了盐岩蠕变本构模型。试验结果表明,稳态蠕变率是偏应力和温度的非线性函数关系,并提出盐岩的稳态蠕变本构方程。通过试验数据拟合得到本构方程中的参数,将理论和试验结果进行比较,结果表明两者的吻合性较好。     

Dynamic Characteristics of Granite Subjected to Intermediate Loading Rate

Summary A large diameter split Hopkinson pressure bar (SHPB) has been developed. This equipment is briefly described, together with a shaped striker that initiates a half-sine incident waveform to obtain the complete stress-strain relationship of the Bukit Timah granite at medium strain rate. Good constant strain rate was derived, and the dynamic complete stress-strain curves and energy absorption of the granite were measured at a strain rate between 20 and 60 per second. Repeated impact between 60–90% of the static strength of the granite was also conducted. Results from the tests show that the cumulative damage of the granite depends on the peak stress of the dynamic loads with a fixed duration. The dynamic fracture strength of the granite loaded at medium strain rate is directly proportional to the cube root of the strain rate. For the granite loaded at this strain rate, Youngs modulus is unchanged. Energy absorption of the samples loading to fragmentation determined its fragmented size distribution. At high strain rate, the rock possesses large energy absorption and the particle size of the fragments is much smaller.     

应力路径和干湿状态对堆石料颗粒破碎的影响研究

堆石料颗粒破碎是引起高土石坝变形的重要因素。在大坝填筑、蓄水期,堆石料的应力路径和干湿状态均是变化的。通过大型三轴试验,系统地研究了应力路径和干湿状态对堆石料颗粒破碎规律的影响。试验结果表明：（1）相同初始条件下,按不同应力路径达到同一轴向应变停机时测定的颗粒破碎率是不同的,等围压?3试验产生的颗粒破碎最大,等平均主应力p试验的次之,等最大主应力?1试验的最小,但不同应力路径下的颗粒级配演化规律是一致的。（2）相同初始条件下,湿样的颗粒破碎率明显高于干样,且二者的差距随着围压的增大而增大,不同干湿条件下的颗粒级配演化规律同样是一致的。（3）建立的考虑母岩强度的颗粒破碎率与塑性功的关系可以较好地统一描述不同应力路径及干湿状态下的颗粒破碎。该研究成果可为建立复杂应力路径及干湿变化条件下考虑颗粒破碎的堆石料弹塑性本构模型提供依据。     

基于考虑动力本构的连续−非连续方法的单轴压缩岩样开裂过程模拟

动载作用下岩石的破坏规律研究对于众多地质灾害的机制分析和预防具有重要的理论及实际意义。鉴于数值模拟研究的优势,应大力发展适于岩石动力断裂过程模拟的数值方法。在自主开发的拉格朗日元与离散元耦合连续−非连续方法的基础上,采用朱−王−唐本构模型取代了广义胡克定律,发展了考虑动力本构的连续−非连续方法,其正确性通过模拟不同加载速度时砂岩试样的单轴压缩试验进行了验证。通过统计裂缝区段数目随着岩样的纵向应变的演化规律,并监测岩样左、右对称线上多个测点的最小主应力的演化规律,开展了不同加载速度时单轴压缩花岗岩试样的变形−开裂过程研究,阐明了岩样的开裂机制。研究发现,剪裂缝以雁列式展布,整体上形成剪切带。随着时步数目的增加,各测点的最小主应力均呈波动下降−震荡上升的变化趋势。震荡上升阶段对应岩样的应变软化阶段。测点分离后最小主应力的震荡幅度较大,这是由于节点分离和单元接触激发了较大的应力波。剪切带尖端的最小主应力集中会使测点发生剪切分离。当岩样的三角块向下楔入时,下方测点的应力状态类似于紧凑拉伸试验进而发生拉伸分离。     

动载确定方法对岩石动态断裂韧度测试的影响

为了考察不同方法确定的动态载荷对测试岩石动态断裂韧度的影响,在SHPB压杆系统上动态冲击直径80 mm的大理岩圆孔裂缝平台巴西圆盘,获得了弹性压杆上的应力波形,间接计算得到3种不同的作用在圆盘端部的动态载荷。将载荷输入ANSYS动态有限元模型中,求得了相应的动态应力强度因子,并根据试验-数值分析方法确定了岩石的动态断裂韧度测试值。结果表明,在加载速率约为4.0×104 MPa•m1/2/s的条件下,采用三波法确定的大理岩的平均动态断裂韧度为 3.92 MPa•m1/2,采用一波法比三波法计算的结果偏低11.22%,采用二波法比三波法计算的结果偏高20.15%,3种方法得到的结果差异较大。应力波在传播过程中,通过圆盘表面和预制裂缝面发生散射,部分能量不断发生释放是造成圆盘试件两端加载载荷不相等的主要原因。三波法是3种方法中比较理想的动态载荷确定方法,但需要考察试件的动态应力平衡性。     

坝基碎屑岩三轴蠕变特性及长期强度试验研究

基于物理和常规变形特性分析,认为某水电站坝基碎屑岩力学特性异常复杂,故采用岩石三轴伺服仪开展了系列的蠕变力学特性试验研究。首先,探讨了轴向、侧向及体积蠕变特性和速率变化规律;其次,分析了蠕变对偏应力-应变特性影响并开展了破坏岩样微细观电镜扫描试验;最后,在采用等时曲线簇法确定长期强度的基础上,以实际延性扩容为依据,认为岩石侧向体积扩容速率大于轴向压缩速率的临界点为快速蠕变破坏的标志,提出确定长期强度的新方法。结果表明,碎屑岩表现出显著的蠕变特性,蠕变曲线主要分为初始衰减和稳态蠕变两阶段,且最后一级应力施加后呈现加速蠕变现象;应力增加导致蠕变变形加剧,最终破坏表现出轴向压缩变形过大、体积延性扩容明显、稳态蠕变速率较大等特点,且蠕变速率与应力符合指数函数关系。岩样长期强度与三向稳态蠕变速率交点和侧向体积扩容应力阀值基本一致,为常规强度的54%～80%。试验结果旨在为相关岩石工程长期稳定分析及蠕变模型构建提供可靠的依据。     

基于连续—非连续方法的地质体材料变形—拉裂过程模拟——以岩样紧凑拉伸试验为例

地质体材料易发生拉裂,为了有效模拟地质体材料的变形-拉裂过程,自主研发了一种连续-非连续方法。该方法实质是拉格朗日元法与虚拟裂缝模型的耦合,既能较准确模拟应力应变场,又能较准确模拟连续介质向非连续介质转化的复杂过程。以岩样紧凑拉伸试验为例开展变形-拉裂过程研究,得到以下结果。紧凑拉伸岩样的变形-拉裂过程:在岩样的Ⅴ形缺口尖端附近出现最大主应力集中现象;节点发生分离,虚拟或真实裂缝扩展,最大主应力始终集中于虚拟裂缝的尖端位置;岩样被拉裂成两部分。最大不平衡力发生1次突增对应着1个节点的分离。在峰值之前,岩样的载荷-位移曲线表现出了硬化现象;随着岩样尺寸的增加,应力-应变曲线的峰值有所下降,这与Bazǎnt的尺度律相一致,且峰后应力-应变曲线的陡峭程度增大。目前针对紧凑拉伸试验的模拟结果是合理的,由此在一定程度上说明了提出的连续-非连续方法在连续介质向非连续介质转化模拟方面的突出能力。