单轴压缩下断续节理岩体动态损伤本构模型

断续节理将对工程岩体的强度及变形等力学特性产生显著影响,损伤力学中视节理为岩体的一种宏观损伤,因而采用损伤张量来刻画其对岩体的影响。目前学术界提出了用节理的几何、强度及变形等3类参数来描述节理的物理力学性质,而目前的岩体损伤张量计算方法都只涉及前2类参数,均没有涉及其变形参数即法向及切向刚度。为此,在前人研究的基础上,基于断裂及损伤理论提出了考虑节理法向及切向刚度的单轴压缩下单条断续节理引起的损伤张量计算公式,进而通过考虑节理间相互作用给出了单组单排或多排节理岩体损伤张量计算公式。其次,以岩石细观动态损伤模型为基础,结合宏细观损伤耦合观点提出了一个能够同时考虑节理几何、强度及变形参数的断续节理岩体动态损伤本构模型。最后,利用该模型讨论了节理参数及载荷应变率等对岩体动态力学特性的影响,认为节理长度减小及摩擦角增大将导致岩体动态峰值强度及弹性模量增大;岩体动态峰值强度及弹性模量则随着节理法向及切向刚度的增大分别减小或增大;而当节理法向及切向刚度按照同一比例增大时,岩体动态峰值强度及弹性模量则是增大的。岩体动态峰值强度与载荷应变率呈正相关。     

节理抗剪强度综合评价的试验研究

选取不同吻合度自然节理试样,进行干燥状态和饱和状态的JRC-JCS模型和JRC-JMS模型试验研究。结果表明,在节理粗糙度系数的定向统计测量和尺寸效应分析,以及剪切过程衰减折减的基础上,运用JRC-JCS模型综合评价节理抗剪强度可靠性较好,考虑吻合度系数的JRC-JMC模型比JRC-JCS模型具更好的预测能力。     

甘肃北山地质处置库围岩节理抗剪强度经验估算

节理抗剪强度参数是地质处置库预选地段工程地质对比和围岩稳定性分析的重要指标。本文在旧井地段英云闪长岩节理分组的基础上,运用定向统计测量方法估测节理粗糙度系数,通过评价节理粗糙度系数的尺寸效应,确定节理抗剪强度经验估算有效长度,由JRC-JCS模型求得各组节理4个方向的抗剪强度参数,并评价了地质处置库围岩节理抗剪强度的各向异性特征。     

真三轴应力状态下岩石损伤本构模型

从岩石微元破坏服从Weibull 随机分布的特点出发,建立了真三轴应力状态下岩石损伤软化统计本构关系;根据应力应变关系的几何条件,建立了本构模型参数与岩石变形过程中应力应变曲线特征参量的理论关系,从而增强了模型的适应性;最后通过试验实测数据验证了模型的合理性.     

深部节理岩体塑性损伤耦合微面模型

采用微面模型理论和损伤力学方法,建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约束模型,根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微面模型的MARC有限元子程序,对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳定问题进行了数值模拟研究. 数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体的各向异性非弹性变形和次生节理演化过程.      

Triaxial compressive behavior of rock with mesoscopic heterogenous behavior: Strain energy density factor approach

The strain energy density factor approach is used in conjunction with a micromechanics model to investigate the condition and direction of shear failure for brittle rock subjected to triaxial compression. Moderate confinement in addition to localized deformation and damage are considered. Quantified are the effects of the various geometric and load parameters that involve the interaction of microcrack, friction and the confining pressure such that the path of the wing crack is taken into account. The influence of all microcracks with different orientations are introduced into the constitutive relation. The closed-form solution for the complete stress–strain relation of rock containing microcracks is obtained. It is shown that the complete stress–strain relationship includes linear, nonlinear hardening, rapid stress drop and strain softening effects. The theoretical results show that deviation of the direction of wing cracks from the line of the pre-existing crack decreases with increasing confinement pressure and friction coefficient. Theoretical predictions and experimental results show good agreement.     

Constitutive relationship of brittle rock subjected to dynamic uniaxial tensile loads with microcrack interaction effects

A micromechanical model is proposed to describe both stable and unstable damage evolution in microcrack-weakened brittle rock material subjected to dynamic uniaxial tensile loads. The basic idea of the present model is to classify the constitution relationship of rock material subjected to dynamic uniaxial tensile loads into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress-strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to consistent with the experimental results.     

Analysis of the localization of deformation and the complete stress–strain relation for mesoscopic heterogeneous brittle rock under dynamic uniaxial tensile loading

Stress redistribution induced by excavation results in the tensile zone in parts of the surrounding rock mass. It is significant to analyze the localization of deformation and damage, and to study the complete stress–strain relation for mesoscopic heterogeneous rock under dynamic uniaxial tensile loading. On the basis of micromechanics, the complete stress–strain relation including linear elasticity, nonlinear hardening, rapid stress drop and strain softening is obtained. The behaviors of rapid stress drop and strain softening are due to localization of deformation and damage. The constitutive model, which analyze localization of deformation and damage, is distinct from the conventional model. Theoretical predictions have shown to consistent with the experimental results.     

不同应变率下煤岩破坏特征及其本构模型

利用直径50 mm的分离式霍普金森压杆,对煤岩展开20～100 s?1动态应变率下的单轴冲击压缩试验,结合高速摄影分析其变形破坏特征,并建立基于Weibull统计分布和Drucker-Prager破坏准则的煤岩动态强度型统计损伤本构模型。试验结果表明:（1）煤岩动态应力-应变曲线存在明显的非线性特征,随应变率升高,动态抗压强度与弹性模量均呈线性增长且增幅显著,破坏形态由低应变率下的轴向劈裂破坏向高应变率下的压碎破坏过渡;（2）在动态应变率20～100 s?1下,煤岩破坏后碎块具有明显的分形特性,破碎块度分维值为1.9～2.2,且随着应变率的升高,煤岩破碎程度增大,碎块块度减小;（3）基于Weibull分布参数F₀、m和应变率的关系,修正煤岩的本构模型,并与试验结果进行对比,验证该模型的合理性。     

冲击载荷下盐岩的力学性能和本构关系研究

利用直径75mm的大尺寸SHPB实验装置开展盐岩的动态压缩性能实验研究,得到了其动态应力应变曲线。分析表明,在冲击压缩载荷作用下,盐岩材料有明显的损伤软化现象和应变率效应。针对实验曲线,通过引入描述材料强度随应变率强化的应变率增强因子和随不可逆变形发展而弱化的损伤因子,提出了含损伤率的相关动态本构模型,拟合得到的本构方程形式比较简单,能够较好地反映盐岩动态实验加载过程的主要特征,具有一定的工程应用价值。     

A dynamic damage constitutive model for a rock mass with non-persistent joints under uniaxial compression

Most problems faced by the practicing rock mass engineering involve the evaluation of rock mass dynamic strength and deformability. As part of a rock mass, the mesoscopic flaws such as the microcracks and the macroscopic ones such as the joints both inherently affect the rock mass dynamic strength and deformational behavior. Nearly none of the existing models can handle the co-effect of these two kinds of flaws on the rock mass dynamic mechanical behavior. This study focusses on the rock mass with multi-sets of non-persistent joints and establishes a mathematical model accounting for the anisotropy in dynamic strength and deformability induced by the joints. Accordingly, an approach incorporating the existing models or methods to enable perfect simulation of the dynamic stress-strain relationship of a rock mass is proposed, in which the joint geometrical parameters such as the joint length and dip angle, the strength ones such as the joint internal friction and the deformational ones such as the joint normal and shear stiffness can all be taken into account. In order to investigate the validity of the proposed model, a series of calculation examples have been made and the results fits very well with the theoretical ones.     

基于时域自适应算法的单向粘弹性节理岩体的等效分析

在粘弹性节理岩体的数值仿真中,必需合理地计及节理的影响.若将节理与岩体作为独立的材料组分考虑,在节理密集的情况下,时空两方面的计算开销可能是很难承受的.本文通过时域展开技术与一个简单的等效假定,得到了递推格式的粘弹性单向节理岩体的等效本构关系,及等效场的递推求解格式,并由此数值模拟了一个粘弹性节理岩体中的洞室,计算结果与考虑独立组分的ANSYS粘弹性分析进行了比较,二者的计算平均误差约为11.03%,但前者的计算效率为后者16.81倍.本文工作有可能为粘弹性节理岩体的数值模拟提供一条新的途径.     

考虑损伤的节理本构模型

本文在弹塑性损伤的理论框架内，讨论了节理等地质间断面的本构模型。这个模型能够反映节理面的损伤弱化，扩容和弹性刚度劣化等复杂特性。这个模型的另一优点是，塑性变形增量与屈服面是非正交的，但本构矩阵具有对称性。这种对称性在岩石力学的理论研究和数值分析中是至关紧要的。     

Nonlinear analysis of beam–column joints using softened truss model

The aim of this study is to expand the application of the nonlinear softened truss model for membrane elements on beam–column joints. The softened truss model employs three equations for equilibrium, three for compatibility and four equations for the constitutive laws of materials. The constitutive equations for both the concrete and steel are based on the actually observed stress–strain relationships. The model has three important attributes. The first is the nonlinear association of stress and strain. The second, and conceivably more noteworthy, is the softening of concrete in compression due to tensile strains in the perpendicular direction. The third is that the influence of the concrete tensile stresses between cracks on the average stress–strain relationship for reinforcing steel and the influence of orthogonal tensile stresses on the compression stress–strain relationship for concrete can be considered in the model. For beam–column joints, one of the most important factors influencing the behaviour is certainly the bond conditions of the beam bars. In this study, the softened truss model is expanded to take into account the influence of this important factor into account. In the revised version of the model, full strain compatibility does not exist between the steel reinforcement and the surrounding concrete and thus the factors influencing the bond-slip between concrete and reinforcement is adequately considered. The improved softened truss model is applied on 51 exterior beam–column joint tests. It is apparent from the results that the revised model gives very accurate predictions of the shear strength of joints and is an improvement on the existing version of the model proposed by Hsu.     

Effect of non-persistent joints distribution on shear behavior

The particle flow code 2D (PFC2D) is used to establish a coplanar, non-persistent joint model. Three joint distribution types, namely, both-side (type a), scattered (type b), and central (type c), are set according to their position. Numerical simulations of the direct shear test are conducted to investigate the effect of non-persistent joint distribution and connectivity on shear mechanical behavior. Simulation results are in good agreement with the analytical solutions to Jennings' criterion, and show: (1) type-c and type-b joints have high strength, whereas type-a joints have low strength. Shear strength and modulus increase with a decrease in joint persistency, and the shear displacement that correspond to shear strength increases with a decrease in persistency. (2) The brittle failure characteristics of the sample are evident when the intact rock bridge area is large. Reinforcement at both ends of the joint limits shear deformation, and shear strength can be effectively improved when joint persistency is large. The small-area dispersed reinforcement joint method cannot effectively improve shear strength. (3) The comprehensive shear strength parameters and the shear strength of the non-persistent joints can be predicted well using Jennings' criterion. Cohesion is the dominant factor that controls shear strength.     

The initiation and growth of adiabatic shear bands

A simple version of thermo/viscoplasticity theory is used to model the formation of adiabatic shear bands in high rate deformation of solids. The one dimensional shearing deformation of a finite slab is considered. For the constitutive assumptions made in this paper, homogeneous shearing produces a stress/strain response curve that always has a maximum when strain and rate hardening, plastic heating, and thermal softening are taken into account. Shear bands form if a perturbation is added to the homogeneous fields just before peak stress is obtained with these new fields being used as initial conditions. The resulting initial/boundary value problem is solved by the finite element method for one set of material parameters. The shear band grows slowly at first, then accelerates sharply, until finally the plastic strain rate in the center reaches a maximum, followed by a slow decline. Stress drops rapidly throughout the slab, and the central temperature increases rapidly as the peak in strain rate develops.     

Progressive failure constitutive model of fracture plane in geomaterial based on strain strength distribution

Progressive failure constitutive model of fracture plane in geomaterial based on strain strength distribution is proposed. The basic assumption is that strain strength of geomaterial comply with a certain distribution law in space. Failure of tensile fracture plane and shear fracture plane in representative volume element (RVE) with iso-strain are discussed, and generalized failure constitutive model of fracture plane in RVE is established considering combined effect of tension and shear. Fracture plane consists of elastic microplanes and fractured microplanes. Elastic microplanes are intact parts of the fracture plane, and fractured microplanes are the rest parts of the fracture plane whose strain have ever exceeded their strain strength. Interaction mode on elastic microplanes maintains linear elasticity, while on fractured microplanes it turns into contact and complies with Coulomb’s friction law. Intact factor and fracture factor are defined to describe damage state of the fracture plane which can be easily expressed with cumulative integration of distribution density function of strain strength. Strong nonlinear macroscopic behavior such as yielding and strain softening can be naturally obtained through statistical microstructural damage of fracture plane due to distribution of strain strength. Elastic–brittle fracture model and ideal elastoplastic model are special cases of this model when upper and lower limit of distribution interval are equal.     

超弹性材料本构关系的最新研究进展

超弹性材料是工程实际中的常用材料, 具有在外力作用下经历非常大变形、在外力撤去后完全恢复至初始状态的特征. 超弹性材料是典型的非线性弹性材料, 其性能可通过材料的应变能函数予以表征. 近几十年来, 围绕应变能函数形式的构造, 已提出许多超弹性材料本构关系研究的数学模型和物理模型, 但适用于多种变形模式和全变形范围的完全本构关系仍是该领域期待解决的重要问题. 本文从3个不同角度, 对超弹性材料本构关系研究的最新进展进行了总结和分析: (1)不同体积变化模式, 包含不可压与可压两种; (2)多变形模式, 包含单轴拉伸、剪切、等双轴以及复合拉剪等多个种类; (3)全范围变形程度, 包含小变形、中等变形到较大变形范围. 超弹性材料本构关系研究的最新进展表明, 为了全面描述具体材料的实验数据并在实际问题中应用超弹性材料, 需要建立适合于多种变形模式和全变形范围的可压超弹性材料的完全本构关系. 对实际超弹性材料完全本构关系的建立及可压超弹性材料应变能函数的构造, 笔者还提出了相应的实施步骤和研究方法.      

STATE OF THE ART OF CONSTITUTIVE RELATIONS OF HYPERELASTIC MATERIALS 1)

超弹性材料是工程实际中的常用材料, 具有在外力作用下经历非常大变形、在外力撤去后完全恢复至初始状态的特征. 超弹性材料是典型的非线性弹性材料, 其性能可通过材料的应变能函数予以表征. 近几十年来, 围绕应变能函数形式的构造, 已提出许多超弹性材料本构关系研究的数学模型和物理模型, 但适用于多种变形模式和全变形范围的完全本构关系仍是该领域期待解决的重要问题. 本文从3个不同角度, 对超弹性材料本构关系研究的最新进展进行了总结和分析: (1)不同体积变化模式, 包含不可压与可压两种; (2)多变形模式, 包含单轴拉伸、剪切、等双轴以及复合拉剪等多个种类; (3)全范围变形程度, 包含小变形、中等变形到较大变形范围. 超弹性材料本构关系研究的最新进展表明, 为了全面描述具体材料的实验数据并在实际问题中应用超弹性材料, 需要建立适合于多种变形模式和全变形范围的可压超弹性材料的完全本构关系. 对实际超弹性材料完全本构关系的建立及可压超弹性材料应变能函数的构造, 笔者还提出了相应的实施步骤和研究方法.