Exploring the influence of interparticle friction on critical state behaviour using DEM

Understanding the extent to which discrete element method (DEM) simulations can capture the critical state characteristics of granular materials is important to legitimize the use of DEM in geomechanics. This paper documents a DEM study that considered the sensitivity of the critical state response characteristics to the coefficient of interparticle friction (μ) using samples with gradings that are representative of a real soil. Most of the features that are typically associated with sand behaviour at the critical state were seen to emerge from the DEM simulation data. An important deviation occurs when high μ values (μ ≥ 0.5) are used, as has been the case in a number of prior DEM studies. While there is a systematic variation in the critical state behaviour with μ for μ < 0.5, when μ ≥ 0.5, the behaviour at the critical state seems to be insensitive to further increases in μ. In contrast to observations of conventional soil response, when μ ≥ 0.5, the void ratio at the critical state initially increases with increasing mean effective stress (p′). Analysis of the DEM data and use of simple models of isolated force chains enabled some key observations. When ‘floating’ particles that do not transmit stress are eliminated from the void ratio calculation, the void ratio at the critical state decreases consistently with increasing p′. There is a transition from sliding to rolling behaviour at the contact points as μ increases. Beyond a limiting value of μ, further increases in μ do not increase the buckling resistance of individual strong force chains. Copyright © 2014 John Wiley & Sons, Ltd.     

The effect of cation ordering and temperature on the high-pressure behaviour of dolomite

Synchrotron single-crystal X-ray diffraction experiments at high-pressure and high-temperature conditions were performed up to 20 GPa and 573.0(2) K on a fully ordered stoichiometric dolomite and a partially disordered stoichiometric dolomite [order parameter, s = 0.26(6)]. The ordered dolomite was found to be stable up to approximately 14 GPa at ambient temperature and up to approximately 17 GPa at T = 573.0(2) K. The P–V data from the ambient temperature experiments were analysed by a second-order Birch–Murnaghan equation-of-state giving K ₀ = 92.7(9) GPa for the ordered dolomite and K ₀ = 92.5(8) GPa for the disordered dolomite. The high-temperature data, collected for the ordered sample, were fitted by a third-order Birch–Murnaghan equation-of-state resulting in K ₀ = 95(6) GPa and K′ = 2.6(7). In order to compare the three experiments results, a third-order Birch–Murnaghan equation-of-state was also calculated for the ambient temperature experiments giving K ₀ = 93(3) GPa, K′ = 3.9(6) for the ordered dolomite and K ₀ = 92(3) GPa, K′ = 4.0(4) for the disordered dolomite. The derived axial moduli show that dolomite compresses very anisotropically, being the c-axis approximately three times more compressible than the a-axis. The axial compressibility increases as T increases, and the a-axis is the most temperature-influenced axis. On the contrary, axial compressibility is not influenced by disordering. Structural refinements at different pressures show that Ca and Mg octahedra are almost equally compressible in the ordered dolomite with K(CaO₆) = 109(4) GPa and K(MgO₆) = 103(3) GPa. On the contrary, CaO₆ compressibility is reduced and MgO₆ compressibility is increased in the disordered crystal structure where K(CaO₆) = 139(4) GPa and K(MgO₆) = 89(4) GPa. Disordering is found to increase CaO₆ and to decrease MgO₆ bond strengths, thus making stiffer the Ca octahedron and softer the Mg octahedron. Cation polyhedra are distorted in both ordered and disordered dolomites and they increase in regularity as P increases. Ordered dolomite approaches regularity at approximately 14 GPa. The increase in regularity of octahedra in the disordered dolomite is strongly affected by the very slow regularization of MgO₆ with respect to CaO₆. The phase transition to the high-pressure polymorph of dolomite (dolomite-II), which is driven by a significant increase in the regularity of both cations polyhedra and mineral crystal structure, occurs in the ordered dolomite at ambient temperature at approximately 14 GPa; whereas no clear evidences of phase transition were observed as regards the disordered crystal structure.     

高烈度地震区隧道减震模型的建立及其减震效果模型试验研究

现场调查表明：5.12汶川大地震共造成21条高速公路和16条国道、省道干线上总延长达20 417 m的24座隧道结构受损,严重影响了灾区救援和灾后重建。可见,传统的隧道抗震技术已不能保证高烈度地震区隧道结构的安全。为此,在分析隧道震害特点及原因的基础上,建立了隧道支护体系内设置减震层的减震模型,根据这一减震模型,研究了减震层刚度、地震波输入频率和减震层阻尼等参数对隧道减震效果的影响。通过三维数值模拟和室内模型试验,验证了该减震模型的减震效果。研究成果可为高烈度地震区隧道结构减震设计提供参考。     

A new SPH-based approach to simulation of granular flows using viscous damping and stress regularisation

The smoothed particle hydrodynamics (SPH) method was recently extended to simulate granular materials by the authors and demonstrated to be a powerful continuum numerical method to deal with the post-flow behaviour of granular materials. However, most existing SPH simulations of granular flows suffer from significant stress oscillation during the post-failure process, despite the use of an artificial viscosity to damp out stress fluctuation. In this paper, a new SPH approach combining viscous damping with stress/strain regularisation is proposed for simulations of granular flows. It is shown that the proposed SPH algorithm can improve the overall accuracy of the SPH performance by accurately predicting the smooth stress distribution during the post-failure process. It can also effectively remove the stress oscillation issue in the standard SPH model without having to use the standard SPH artificial viscosity that requires unphysical parameters. The predictions by the proposed SPH approach show very good agreement with experimental and numerical results reported in the literature. This suggests that the proposed method could be considered as a promising continuum alternative for simulations of granular flows.     

Modelling the effect of temperature on unsaturated soil behaviour

A simple thermohydromechanical (THM) constitutive model for unsaturated soils is described. The effective stress concept is extended to unsaturated soils with the introduction of a capillary stress. This capillary stress is based on a microstructural model and calculated from attraction forces due to water menisci. The effect of desaturation and the thermal softening phenomenon are modelled with a minimal number of material parameters and based on existing models. THM process is qualitatively and quantitatively modelled by using experimental data and previous work to show the application of the model, including a drying path under mechanical stress with transition between saturated and unsaturated states, a heating path under constant suction and a deviatoric path with imposed suction and temperature. The results show that the present model can simulate the THM behaviour in unsaturated soils in a satisfactory way.     

The combined effect of moisture and temperature on the anomalous expansion behaviour of marble

Three calcitic marble types often used as cladding material and different in lattice preferred orientation (texture), grain shape preferred orientation, grain size distribution and grain interlocking are investigated to study the combined effect of thermal cycles (day-night) and moisture on the decay of marble, particularly on the bowing phenomenon. Repetitive heating-cooling under dry conditions leads to considerable inelastic residual strain only in the first thermal cycle. The residual strain continuously increases again if water is present, whereby the moisture content after a thermal cycle has a certain impact on the decay rate. The water-enhanced thermal dilatation strongly correlates with the deterioration rate obtained from the laboratory bow test. All applied approaches reveal that the texture in combination with the grain shape preferred orientation control the intensity and anisotropy of marble deterioration, thus, the cutting direction of facade panels has to be considered in terms of durability. On-site analysis from building facades confirm the laboratory results such as the bowing tendency of different marbles, directional dependence, relevance of moisture content during cyclical heating and the loss of strength due to environmental impact.Special Issue: Stone decay hazards     

The friction and deformation behaviour of rock joints

SummaryThe Friction and Deformation Behaviour of Rock Joints The present investigation deals with the influence of joint roughness and material strength on the friction and deformation behaviour of rough indented joints. Besides the analytical determination of the joint roughness by the dilatation behaviour of the joint, friction tests on models with natural joint morphology and rock samples have been conducted.A material law is developed, describing the friction resistance and the dilatation behaviour at large deformations in their dependence on normal stress and material strength.

---

ZusammenfassungDas Reibungs- und Verformungsverhalten von Klüften Die vorliegende Untersuchung war dem Einfluß der Kluftflächenunebenheit und der Materialfestigkeit auf das Reibungs- und Verformungsverhalten rauher, verzahnter Kluftflächen gewidmet.Neben der analytischen Bestimmung der Kluftflächenunebenheit anhand des Dilatationsverhaltens der Kluftfuge wurden Reibungsversuche an Modellkörpern mit natürlicher Kluftmorphologie und an Gesteinen durchgeführt. Es wird ein Stoffgesetz entwickelt, mit welchem der Reibungswiderstand und das Dilatationsverhalten in Abhängigkeit von der Normalspannung und der Materialfestigkeit bei längeren Gleitwegen beschrieben werden kann.

---

RésuméLe comportement des fissures au frottement et aux déformations La présente étude traite de l'influence de la rugosité d'une diaclase et de la résistance de la roche sur le comportement au frottement et aux déformations d'une fissure âpre et endentée. La rugosité de la fissure a été déterminée d'après la dilatation maximale. Des éssais de frottement ont été exécutés à l'aide de modèles ayant une surface de fissure naturelle et à l'aide d'échantillons de roche. Une loi a été développée qui décrit la resistancé au frottement et le comportement à la dilatation de la fissure pour de grandes déformations en fonction de la contrainte normale et de la résistance des matériaux.

---

List of Symbols h dilatation - i deviation angle, dilatation angle - s shear deformation - residual friction angle - normal stress - _F tensile strength - shear stress With 15 Figures     

饱水和天然状态下页岩滞后效应及阻尼特性研究

运用RLW-2000型微机伺服岩石三轴试验机,对饱水状态和天然状态页岩在不同围压下进行三轴循环加卸载试验,分析了两种含水状态页岩的力学特性和滞后效应,并基于能量原理讨论了阻尼比演化规律。试验结果表明:随循环次数增加,累积残余应变逐渐增加,相对残余应变先降低,后趋于稳定区域,直至破坏前急剧增加;饱水页岩的加卸载变形模量均比天然页岩小,加载变形模量整体比卸载变形模量小;在加载和卸载阶段,均出现应变始终滞后于应力的现象;饱水页岩的滞后效应比天然页岩更明显。提出了考虑滞后效应的能量计算方法,较以往能量计算的误差更小。最后基于能量原理对阻尼比的计算公式进行修正,发现饱水页岩的阻尼比比天然页岩大。阻尼比的变化较好地反映了页岩损伤机制,可作为预判页岩失稳破坏的重要指标。     

The combined effect of pH and temperature on smectite dissolution rate under acidic conditions

The main goal of this paper is to propose a new rate law describing the combined effect of pH (1 to 4.5) and temperature (25 to 70 °C) on smectite dissolution rate, under far from equilibrium conditions, as a step towards establishing the full rate law of smectite dissolution under acidic conditions. Dissolution experiments were carried out using non-stirred flow-through reactors fully immersed in a thermostatic water bath held at a constant temperature of 25.0°C, 50.0°C or 70.0°C ± 0.1°C. Smectite dissolution rates were obtained based on the release of silicon and aluminum at steady state. The results show good agreement between these two estimates of smectite dissolution rate. Low Al/Si ratios were obtained in experiments that were conducted at pH ≥4. These low Al/Si ratios are explained by precipitation of gibbsite and/or diaspore.Dissolution rate increases with temperature and decreases with increasing pH. Dissolution rates of experiments in which ΔG_r ≤ −21 kcal mol ⁻¹, are not affected by deviation from equilibrium. Dissolution rates in most experiments are not affected by the addition of up to 0.3 M NaNO₃ to the input solution.A simple model is used to describe the combined effect of pH and temperature on smectite dissolution rate. According to this model, dissolution rate is linearly proportional to the concentration of adsorbed protons on the mineral surface, and proton adsorption is described using a Langmuir adsorption isotherm. All experimental results at pH <4 were fitted to the model using a multiple non-linear regression. The resulting rate law is:

(A1)     

Experimental and modelled mechanical behaviour of a rock fracture under normal stress

Summary  The mechanical and hydromechanical behaviour of isolated rock joints is of prime importance for a correct understanding of the behaviour of jointed rock masses. This paper focuses on the mechanical behaviour of a fracture under normal stress (fracture closure), using approaches based on both experimentation and modelled analysis. Experimental closure tests were carried out by positioning four displacement transducers around a fracture, leading to results which tended to vary as a function of transducer location. Such variations can be explained by the non-constant void space distribution between both walls of the fracture. The present study focuses on the importance of transducer location in such a test, and on the significant role played, in terms of mechanical response, by the morphology of the fracture surfaces. An analytical mechanical model is then developed, which takes into account the deformation of surface asperities and of the bulk material surrounding the fracture; it also includes the effects of mechanical interaction between contact points. The model is validated by simulating the behaviour which is very similar to experimental observations. Various parametric studies (scale effect, spatial distribution of contact points) are then carried out. The study of scale effects reveals a decrease in the normal stiffness with increasing fracture size. Finally, analysis of the role of various mechanical parameters has shown that the most influential of these is Young’s modulus corresponding to the bulk material surrounding the joint. Many applications, such as geothermal fluid recovery from fractures, could benefit from these results. Correspondence: Antoine Marache, Université Bordeaux 1, GHYMAC, Av. des Facultés, 33405 Talence Cedex, France     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

基于原位地应力测试及流变模型的深部泥页岩储层地应力状态研究

深部泥页岩储层地应力状态的准确确定是页岩气等非常规能源高效开发的关键。综合基于原位地应力测试获得水平最小主应力,建立基于流变模型的地应力剖面,应用成像测井技术确定水平最大主应力方向等,是准确确定泥页岩储层地应力的有效方法。将该研究思路应用于陕西汉中SZ1井,利用水压致裂原地应力测试方法获得储层水平最小主应力值范围为32~41 MPa;利用偶极声波测井数据获得岩石力学参数,结合地壳应变率和储层埋藏史,建立了SZ1井地应力剖面,结果表明牛蹄塘组1950~2025 m深度范围内水平主应力差介于10~15 MPa,水平最小主应力值范围为28~41 MPa,水平最大主应力值范围为47~49 MPa,预测得到的水平最小主应力值与实测结果具有较好的一致性。原地应力实测及流变模型预测结果揭示SZ1井地应力为正断型(Sv＞SH＞Sh)或正断型与走滑型相结合的应力状态(Sv≈SH＞Sh)。水平主应力差随伽玛值的升高而变小,表明地应力剖面与地层岩性具有较好的对应关系。基于成像测井揭示的钻孔诱导张裂隙分布特征,SZ1井水平最大主应力方向约为N74°W,与区域构造应力场方向基本一致。相关结论为准确认识SZ1井目标层地应力状态,以及后期水平井布设及压裂控制等提供了重要依据。     

挤土作用下Y型桩侧摩阻力分布模型研究

Y型沉管灌注桩是在传统圆形桩基础上通过改变截面形状以增大桩土接触面积、增加侧摩阻力进而增加承载力的一种新桩型。由于桩身截面的特殊性,在振动沉管的过程中Y型沉管灌注桩截面不同处对土有不同程度的挤压作用,其侧摩阻力分布性状比圆形桩复杂,若将Y型桩侧摩阻力看作沿桩身矩形分布、沿桩周均匀分布与实际不符,工程经验亦表明其不适。通过孔扩张理论与Y型截面的解析方程相结合,推导Y型桩侧摩阻力沿桩周及桩身的不均匀分布模型,同一截面最大侧摩阻力为最小侧摩阻力的1.24~1.75倍,同时分析外包圆半径R、模板弧度θ、桩长L对侧摩阻力的影响。将侧摩阻力沿Y型截面周边的9个区间积分,然后沿桩长积分,得到Y型桩总侧摩阻力,结合端阻力模型计算了Y型桩单桩极限承载力,与静载试验数据吻合。     

The nature and patterns of damping of stress waves in soils

This paper considers the main patterns of damping of stress waves in disperse and semi-rock soils. We determined two strain thresholds, beyond which the damping ratio qualitatively changes: in the region of relative deformations up to 10⁻⁴ the damping ratio can be considered to be constant and is caused by non-ideal elasticity and viscous losses; beyond this strain threshold the damping monotonically increases due to initial plastic soil deformations; in the region of moderate-level shear deformations the increase of the damping ratio slows for relative shear deformations higher than 0.002, primarily because of the rapid decrease of the fraction of the damping due to non-ideal soil elasticity and heat losses in the course of degradation of structural bonds.     

温度和有效应力对砂岩渗透性影响的试验研究

为探讨温度和有效应力对渗透性的影响,在不同温度水平和不同有效应力水平下进行了砂岩的渗透率试验。根据试验结果,初步分析了温度和有效应力对砂岩渗透率及其水力传导系数的影响。试验结果表明,砂岩的渗透率及其水力传导系数都是温度和有效应力的函数:在温度一定的条件下,砂岩的渗透率和水力传导系数均随有效应力的增加而呈负指数规律减小;在有效应力一定的条件下,砂岩的渗透率随温度的升高而减小,但水力传导系数与温度之间的关系函数并非单调函数。因此,在工程实践中不考虑温度和有效应力对岩石 (体)渗透性的影响是不合适的。     

初始剪应力与加荷速率共同作用下饱和软黏土孔压模型试验研究

通过GDS循环三轴试验系统,对循环荷载作用下饱和软黏土的孔压变化规律进行了研究,分析了循环应力比,初始剪应力与振动频率对孔隙水压力的影响。研究结果表明：随着循环应力比的增加,孔压发展速度增快。循环荷载作用下饱和软黏土存在临界循环应力比,通过孔压也可以确定其值大小。在循环初期,孔压率较大,随着循环时间的增加,孔压率逐渐减小。随着循环应力比的增加,孔压率增加。振动频率对孔压比-循环次数关系影响明显,随着频率的增加,孔压比减小;然而,当振动频率大于1 Hz时,振动频率对孔压比-时间与孔压率-时间关系影响不明显。随着初始剪应力的增大,孔压增加。初始剪应力对应变率具有显著影响;随着初始剪应力的增加,应变率增加。在对数坐标下,孔压率与时间呈线性关系。在上述试验基础上建立了孔压率与时间关系表达式,通过积分得到了循环荷载作用下饱和软黏土的动孔压模型     

温度和含水状态对岩石劈裂强度影响的试验研究

岩土类材料的强度强烈依赖于含水率和温度条件.以兰州兰山砂岩为研究对象,分别进行了不同温度（+20℃、-5℃、-10℃）和不同含水率（干燥、天然含水率和饱和状态）条件下岩石的巴西圆盘劈裂试验.试验结果表明,常温下,岩石抗拉强度随含水率增加而急剧减小,试样饱和时软化、崩解,丧失承载能力;在不同负温条件下,天然含水状态的试样抗拉强度最大,干燥状态下最小;在-10~+20℃范围内,干燥岩石强度随温度升高,抗拉强度增大,含水岩石均是随温度升高,抗拉强度减小,但饱和岩石在-5℃时抗拉强度最大.含水率和温度对岩样强度的影响存在临界值,超过临界值,岩样强度随上述因素反向变化.试验结论为岩土类材料劈裂强度的标准化测试及其工程应用提供了重要的基础数据参考.     

Effects of hydraulic shear stress and rate of erosion on the magnitude,degree, and rate of collapse

Arid regions worldwide are plagued by collapsible soils. Collapsible soil is characterised by the sudden decrease in volume that occurs when it is subjected to inundation under constant stress. This volume change manifest itself as drastic and unpredicted foundation settlement, which may lead to further catastrophic failure of the supported structures. Collapse settlement is the term applied to the additional settlement of a foundation due to wetting of the underlying soils. The results of an experimental investigation of the effects of the saturation of soil with water, kerosene, and crude oil, and of the effects of the fluid head on the magnitude, degree, and rate of collapse of the underlying soil are presented in this paper. Soil erodibility is presented in terms of the applied hydraulic shear stress and the rate of erosion. The relationship between soil erosion and the magnitude and rate of collapse is examined. Empirical methods for the prediction of the magnitude and rate of collapse of a soil saturated with the test fluids and subjected to a hydraulic constant head are proposed.