非规则砂岩节理剪切变形本构关系试验研究

岩石节理剪切变形对岩体工程的安全性和稳定性具有重要影响。为研究常法向应力下岩石节理剪切变形本构关系,采用RDS-200型岩石直剪仪对非规则砂岩节理进行了不同法向应力下的直剪试验。根据岩石节理剪切应力?位移全程曲线形状特征,将其依次划分为峰前压密阶段、线性阶段、屈服阶段和峰后软化阶段;根据剪切应力在峰后软化阶段降低幅度和速率大小,将岩石节理剪切应力?位移曲线划分为3种类型：峰后平台型、峰后缓降型和峰后跌落型。基于岩石节理剪切应力?位移曲线各阶段剪切变形特征,采用分段函数建立了岩石节理剪切变形本构模型。与其他模型相比,新提出的岩石节理剪切变形本构模型对试验数据拟合精度更高,更好地描述了岩石节理剪切应力?位移全程曲线。另外,在通过岩石节理直剪试验由经验公式确定模型参数之后,所提出本构模型可在不同法向应力下实现对不同粗糙度岩石节理剪切应力?位移曲线的预测。研究结果对岩石节理剪切变形的数值模拟和工程估计具有一定的实用价值。     

砂土与硫酸腐蚀混凝土接触面剪切试验研究

为了研究酸腐蚀作用对桩-土接触面剪切特性的影响,利用大型多功能界面剪切仪,进行了砂土与被硫酸加速腐蚀0、31、93、154 d的混凝土板之间接触面剪切试验,量测了不同腐蚀时间的混凝土与砂土接触面的剪切应力和剪切位移关系。研究结果表明:砂土与硫酸腐蚀混凝土接触面剪切应力与剪切位移符合双曲线模型,接触面剪切破坏符合摩尔-库仑准则。随着腐蚀时间的增加,接触面黏聚力减小,摩擦角、剪切强度以及达到峰值剪切应力所需的剪切位移增大,接触面的破坏趋向于土体自身剪切破坏。取未腐蚀混凝土与砂土接触面达到峰值剪应力时的剪切位移为基准,对该位移对应的不同腐蚀时间下混凝土与砂土接触面的剪切应力进行对比分析,发现腐蚀5个月后的剪切应力要小于未腐蚀时的剪切应力。按双曲线模型对试验测得的接触面剪切应力和剪切位移关系进行了拟合,获得了该模型中的参数,为进一步进行数值模拟提供参考。     

考虑围岩应变软化及采空区接触的深部煤层底板破坏分析

为研究深部煤层开采底板破坏形态,提出考虑围岩应变软化和采空区接触的FLAC3D有限差分数值方法,以河北开平煤田林西矿2023工作面底板实测导水裂隙带为工程背景,结合朗肯土压力理论定性分析,研究深部煤层底板破坏特征。结果表明:采用应变软化本构关系代替常用摩尔–库伦本构关系能够对围岩塑性破坏后的力学状态更准确表述;采用“应变软化–空–弹性”模型转变的方法,达到模拟采空区顶板垮落后应力传递的效果,弥补了以往煤层开采模拟中采空区垮落后顶底板不接触的固有缺陷;通过采空区顶底板接触与否条件下应力、位移的对比,发现采空区是否接触对数值结果影响巨大,突出考虑采空区接触的必要性;根据模拟结果中塑性剪切应变率的变化,实现了底板滑移面的三维显示,形态为斜向采空区的半包围面状结构;结合朗肯土压力理论将底板塑性区与主动区、过渡区和被动区对应,3个区破坏形式分别为剪切破坏、剪切破坏、拉张与剪切的交互破坏。提出的考虑围岩应变软化及采空区接触的FLAC3D数值方法对煤层开采模拟实现了优化,并可为其他大变形后需考虑接触的工程模拟提供参考。      

青藏冻结粉土与玻璃钢接触面本构模型研究

冻土与构筑物基础接触面的本构模型是分析冻土区构筑物与冻土相互作用、评价工程安全稳定性的基础和关键.为了降低切向冻胀力对基础的上拔作用,防止发生冻拔失稳破坏,青海-西藏±400 kV直流联网工程基础广泛采用玻璃钢覆盖基础表面,以消减冻胀力.为了合理描述玻璃钢与冻土接触面力学特性,采用应变直剪仪开展了青藏冻结粉土与玻璃钢基础接触面直剪试验,获取了不同冻结温度条件下接触面剪应力-位移曲线,分析总结了接触面的强度变化规律和受力变形特性,基于试验得到的接触面本构规律建立了耦合温度效应的冻结粉土与玻璃钢基础接触面剪切应力-位移关系.模型预测结果与试验结果的比较表明,预测结果与试验结果吻合良好,该模型能够较好的描述接触面的力学特性.     

基于扰动理论修正的桩-土接触面荷载传递模型 及其应用

利用硬化模型、双曲线模型、指数模型以及软化模型模拟的桩-土接触关系存在参数取值困难、误差大的问题。为深入研究土-结构接触面的强度与变形机理,基于扰动理论,假定完整状态接触单元的抗剪强度服从线弹性模型,而扰动部位则服从塑性模型,建立修正的桩-土接触面荷载传递模型。该模型参数分析表明,参数k、η对模型τ-s曲线形态影响大,而参数τ、ζ对模型τ←s曲线形态影响很小;并通过大型桩-土接触面室内直剪试验,量化接触面上剪切应力与剪切位移的关联性,进一步确定修正桩-土接触模型内部计算参数。结果表明模型τ-s曲线与试验曲线吻合较好,验证了模型的合理性。扰动桩-土接触模型既能描述桩侧应变软化也能描述硬化特性,有助于理解复杂应力条件下桩-土接触面的强度计算与变形机理。     

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

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

岩土类材料率相关性及硬化-软化特性模型研究

岩土类材料具有显著的率相关性,其应力-应变关系全过程一般表现为硬化-软化特征,研究如何描述这两类特征具有重要的意义。基于岩土类材料的物理变形机制分析,考虑到其结构性特征,采用弹黏塑性本构模型,引入反映材料结构变化的可变结构黏滞系数及其演化方程,并将其嵌入到弹黏塑性本构模型中去,建立能够反映材料结构变化的弹黏塑性本构模型。利用玄武岩常规三轴常应变率试验结果对所建立的模型进行试验验证,结果表明：所建立的弹黏塑本构模型能够完整、合理地描述岩土类材料应力-应变关系全过程以及材料本构关系的率相关性,同时还能够反映材料的峰值强度、残余强度特征。     

砂与粘土接触界面的力学特性试验

利用DSJ-2型电动四联等应变直剪仪分别进行了不同粒径砂和粘土(方案1)、不同含水量砂和粘土(方案2)的直剪试验,揭示了其接触面的力学特性。试验结果表明:法向应力(σ_n)、砂粒径大小、砂含水量对接触面的力学特性有着重要的影响。方案1和方案2接触面的抗剪强度和达到剪切应力峰值时的剪切位移随着法向应力的增大而增大,当σ_n为100 kPa时,粘土和不同粒径砂的接触面剪切应力和水平剪切位移(τ-δ)关系曲线呈现应变软化现象;方案1接触面抗剪强度随砂粒径的减小而降低,方案2在σ_n为100 kPa时,接触面抗剪强度随砂含水量增加而降低,σ_n>100 kPa时,接触面抗剪强度随砂含水量增加先降低后变大;方案2在低法向应力下,接触面抗剪强度对应的含水量敏感区间为[10%,15%],且较明显。      

岩石损伤软化统计本构模型及参数确定方法的新探讨

基于现有岩石损伤软化统计本构模型研究,通过探讨岩石损伤软化统计本构模型参数与岩石应力-应变全程曲线特征参数即峰值应力与应变的关系,建立起特定围压下模型参数与围压的解析表达式。引进岩石Mohr-Coulomb强度准则,建立不同围压下岩石应力-应变全程曲线峰值应力与围压之间的关系,再通过探讨不同岩石应力-应变全程曲线峰值应变与围压的关系,导出了具有普遍意义的不同围压下岩石峰值应变计算公式,从而建立岩石损伤软化统计本构模型参数确定的新方法,由此得到能够模拟不同围压下岩石应变软化全过程的统一损伤软化统计本构模型。该模型较同类模型具有参数少和易于确定等特点,理论计算和实测结果比较分析表明了该方法与模型的合理性。     

基于扩展Lagrange乘子的Clough接触面模型及应用

双曲线接触面本构模型能比较真实地模拟堆石材料与混凝土面板间的接触关系,扩展Lagrange乘子法能比较精确地计算接触面间接触状态。在吸取Clough双曲线接触面本构关系应用于无厚度的Goodman单元和有厚度的Desai薄层单元的成功经验的基础上,将Clough双曲线接触面模型引入扩展Lagrange乘子法进行摩擦接触问题的求解。介绍了基于扩展Lagrange乘子法的序列二次规划法（SQP）提法,详细推导了在扩展Lagrange算法的非线性接触计算中引入Clough双曲线接触面本构关系的数值列式和实施步骤。数值算例和工程实例的计算结果证明该方法是成功的,能够比较真实地模拟堆石-混凝土墙（混凝土面板）之间的接触状态。     

Numerical studies of shear banding in interface shear tests using a new strain calculation method

Strain localization is closely associated with the stress–strain behaviour of an interphase system subject to quasi‐static direct interface shear, especially after peak stress state is reached. This behaviour is important because it is closely related to deformations experienced by geotechnical composite structures. This paper presents a study using two‐dimensional discrete element method (DEM) simulations on the strain localization of an idealized interphase system composed of densely packed spherical particles in contact with rough manufactured surfaces. The manufactured surface is made up of regular or irregular triangular asperities with varying slopes. A new simple method of strain calculation is used in this study to generate strain field inside a simulated direct interface shear box. This method accounts for particle rotation and captures strain localization features at high resolution. Results show that strain localization begins with the onset of non‐linear stress–strain behaviour. A distinct but discontinuous shear band emerges above the rough surface just before the peak stress state, which becomes more expansive and coherent with post‐peak strain softening. It is found that the shear bands developed by surfaces with smaller roughness are much thinner than those developed by surfaces with greater roughness. The maximum thickness of the intense shear zone is observed to be about 8–10 median particle diameters. The shear band orientations, which are mainly dominated by the rough boundary surface, are parallel with the zero extension direction, which are horizontally oriented. Published in 2007 by John Wiley & Sons, Ltd.     

岩石应变软化本构模型建立及NR-AL法求解研究

针对岩土工程材料应变软化问题及有限元对其数值计算时切线刚度矩阵负定造成求解困难的问题进行研究。建立了基于Drucker-Prager（D-P）强度准则的岩石弹塑性应变软化本构模型,本构积分算法采用一种完全隐式返回映射算法,它具有无条件稳定和精确的特点,详细论述了如何进行本构模型的程序化求解;考虑弧长法在判断切线刚度矩阵正定性导致效率低的缺点,在弹塑性增量有限元方程的迭代计算中尝试采用Newton-Raphson法和arc-length法（NR-AL法）联合迭代求解的思路,即在结构未达到极限荷载前采用NR迭代法,而当结构接近极限荷载时转换为AL法控制迭代,从而使结构越过峰值点进入软化区直至破坏,NR-AL法汲取了2者迭代求解中具有的优势;利用C++语言对所建应变软化模型的本构求解和弹塑性增量有限元方程迭代求解过程给予程序实现,应用所编程序进行数值计算,分析了D-P理想弹塑性模型、应变软化模型、应变硬化模型计算的应力-应变曲线的区别,同时将应变软化模型计算结果与试验数据进行了对比。研究结果表明：所建应变软化本构模型可以较好地模拟岩石材料的峰后软化特性,能够揭示峰后应变软化特性和破坏机制,同时NR-AL法能够求解由于应变软化造成的负刚度问题,也克服了单独使用弧长法时判断切线刚度矩阵正定性效率低的缺点。     

Disturbed state model for sand–geosynthetic interfaces and application to pull-out tests

Successful numerical simulation of geosynthetic-reinforced earth structures depends on selecting proper constitutive models for soils, geosynthetics and soil–geosynthetic interfaces. Many constitutive models are available for modelling soils and geosynthetics. However, constitutive models for soil–geosynthetic interfaces which can capture most of the important characteristics of interface response are not readily available. In this paper, an elasto-plastic constitutive model based on the disturbed state concept (DSC) for geosynthetic–soil interfaces has been presented. The proposed model is capable of capturing most of the important characteristics of interface response, such as dilation, hardening and softening. The behaviour of interfaces under the direct shear test has been predicted by the model. The present model has been implemented in the finite element procedure in association with the thin-layer element. Five pull-out tests with two different geogrids have been simulated numerically using FEM. For the calibration of the constitutive models used in FEM, the standard laboratory tests used are: (1) triaxial tests for the sand, (2) direct shear tests for the interfaces and (3) axial tension tests for the geogrids. The results of the finite element simulations of pull-out tests agree well with the test data. The proposed model can be used for the stress-deformation study of geosynthetic-reinforced embankments through numerical simulation. Copyright © 1999 John Wiley & Sons, Ltd.     

接触角度对黄土-三趾马红土界面剪切力学特性影响研究

西北黄土高原第四系黄土广泛沉积于新近系三趾马红土之上形成粗糙接触的异质土界面,为典型的易滑层面。为探讨接触界面粗糙度对黄土-三趾马红土界面剪切力学特性影响,研制界面制样装置及剪切仪,开展简化黄土-三趾马红土界面直剪试验研究。结果表明:界面剪切破坏模式有齿间滑动、齿间滑动-齿面剪断、齿面剪断3种,界面接触角度越大,破坏模式越趋于齿面剪断,接触角度越小,破坏模式越趋于齿间滑动;界面剪切应力-剪切位移曲线演化规律表明界面脆性剪切破坏特征明显,且界面接触角度越大,峰值前剪切刚度与剪切破坏位移越大,峰值后剪切位移“跳跃”跌落现象越明显,界面脆性剪断破坏特征越显著;界面剪切过程产生明显剪胀效应,随界面接触角度增大,峰值剪胀角呈先减小而后增大趋势,反映了界面不同剪切破坏模式变化;受界面间初始黏聚强度与剪切破坏模式影响,界面抗剪强度随法向应力呈非线性变化,并受界面接触角度影响,峰值强度随界面角度增大呈先增大而后减小趋势,残余强度随界面角度增大呈增大趋势。     

填埋场衬垫系统界面本构模型研究

典型垃圾填埋场衬垫系统由土/土工膜、土工膜/土工织物等多个界面组成,这些界面抗剪强度低,成为潜在的滑裂面,容易使填埋场发生沿界面的失稳破坏。本文根据大量土/土工膜界面直剪试验的研究,针对该界面的变形机理和应变软化特性,在前人研究的基础上确定了接触面的厚度,定义了接触面的切向应变和法向应变,将界面的变形分为3个阶段:弹性阶段,塑性阶段和残余阶段。同时提出了适合于土与土工膜界面的弹塑性本构模型。该本构模型参数较少,物理意义明确,运用简单方便。最后,结合Filz接触面直剪试验和钱学德单剪试验的结果,将通过该模型得到的拟合结果与之进行了对比验证分析。研究结果表明:本文提出的土与土工膜界面弹塑性本构模型是合理可行的,且建模方便,适用于工程实践。     

The influence of the plastic potential on plane strain failure

The influence of the shape of the plastic potential in the deviatoric plane on plane strain collapse is investigated. The most commonly employed elastic‐perfect plastic models are considered, which adopt well‐known failure criteria for defining the yield and plastic potential surfaces, namely the von Mises, the Drucker–Prager, the Tresca, the Mohr–Coulomb and the Matsuoka–Nakai criteria. Finally, the conclusions are also extended to strain hardening/softening models. For simple constitutive models based on perfect plasticity, it is shown that the value of the Lode's angle at plastic collapse in plane strain conditions strongly depends on the specific failure surface adopted for reproducing the plastic potential surface. If the value of the Lode's angle at yield coincides with the failure value prescribed by the plastic potential, the stress–strain curves exhibit the typical perfect plastic behaviour with yield coinciding with failure, otherwise the stress changes after yield and the stress‐strain curves resemble those of strain hardening/softening models. The infinite strength which is in some situations exhibited by the Drucker–Prager model in plane strain condition is investigated and explained, and it is shown that this can also affect strain hardening/softening models. Copyright © 2013 John Wiley & Sons, Ltd.     

卸荷条件下岩石变形特征及本构模型研究

基于岩石试件的卸荷试验,研究了卸荷条件下岩石的应力—应变全过程曲线和破裂特征,研究表明：卸荷过程中岩石向卸荷方向回弹变形强烈、扩容显著、脆性破坏特征明显;卸荷条件下岩石破坏具有较强的张性破裂特征,各种级别的张裂隙发育,其剪性破裂面追随张拉裂隙发展。基于试验认识,将卸荷条件下岩石的本构模型分为4段：弹性段、卸荷屈服段、峰后脆性段及残余理想塑性段;通过体积应变εv将Griffith和Mohr Coulomb屈服准则结合起来,建立卸荷条件下岩石的屈服准则,并在应变空间建立相应的卸荷非线性屈服段本构方程;运用应变软化理论,建立卸荷岩石峰后脆性段本构方程。     

Method of Numerical Modeling for Constitutive Relations of Clay

INTRODUCTIONThe mechanical responses of soils are more com-plicated compared with metals.By comparing thephysical and mechanical properties of the metals withsoils,Lade(1988)found that there are17differentpoints between the metals and soils,which differfrom metals inthree basic mechanical characteristics:pressure sensitivity,dilatancy,and dependence ofstress path.Wang(2004)proposed a principle of theinteraction between plastic volume and shear strains,that is,there are two relatively inde…     

两体接触面剪切力学行为的三维数值分析

研究岩石-岩石节理或混凝土-岩石接触面的力学行为,对于评价岩体的稳定性及工程体与岩基的稳定性有重要的意 义。采用ANSYS有限元软件构建了在宽度方向一致的三维粗糙表面,模拟计算了Bandis等在玄武岩试件上的直剪试验结 果,并将数值计算结果与试验结果对比,验证了ANSYS在分析两体剪切力学行为方面的可行性。此外对于有规则起伏的锯齿界面,同样进行了三维直剪试验的数值模拟,获到了不同界面和应力边界条件下的剪切强度,并与Patton剪胀公式理论值进行比较,从理论上再次验证了ANSYS模拟直剪试验的适用性。同时重点分析了直剪试验下的剪切位移曲线、接触面的剪应力分布、剪胀作用;并初步探讨了ANSYS软件自身设定凝聚力对接触面力学性能的影响。     

An elastoplastic constitutive model for frozen saline coarse sandy soil undergoing particle breakage

The mechanical property of frozen saline sandy soil is complicated due to its complex components and sensitivity to salt content and confining pressure. Thus, a series of triaxial compression tests were carried out on sandy samples with different Na₂SO₄ contents under different confining pressures to explore the effects of particle breakage, pressure melting, shear dilation and strain softening or hardening. The test results indicate that the stress–strain curves exhibit strain softening/hardening phenomena when the confining pressures are below or above 6 MPa, respectively. A shear dilation phenomenon was observed in the loading process. With increasing confining pressure, the strength firstly increases and then decreases. By taking into consideration the changes between the grain size distributions before and after triaxial compression tests, a failure strength line incorporating the influences of both particle breakage and pressure melting is proposed. In order to describe the deformation characteristics of frozen saline sandy soil, an elastoplastic incremental constitutive model is established based on the test results. The proposed model considers the plastic compressive, plastic shear and breakage mechanisms by adopting the non-associated flow rule. The breakage mechanism can be reflected by an index related to the initial, current and ultimate grain size distributions. The hardening parameters corresponding to compressive and shear mechanisms consider the influence of particle breakage. Then the effect of particle breakage on both the stress–strain and volumetric strain curves is analyzed. The calculated results fit well with the test results, indicating that the developed constitutive model can well describe the mechanical and deformation features of frozen saline sandy soil under various stress levels and stress paths. In addition, the strain softening/hardening, contraction, high dilation and particle breakage can be well captured.