基于裂纹扩展脆性岩石动力压缩力学特性研究

动态压缩荷载作用下,脆性岩石内部动态细观裂纹扩展特性,对岩石宏观动态力学特性有着重要的影响。然而,对岩石内部动态细观裂纹扩展与宏观动态力学特性的关系研究较少。基于准静态裂纹扩展作用下的应力-应变本构模型、准静态与动态裂纹扩展断裂韧度关系、裂纹速率与应变率关系模型及应变率与动态断裂韧度关系,提出了一种基于细观力学的动态应力-应变本构模型。其中裂纹速率与应变率关系,是根据裂纹长度与应变关系的时间导数推出;应变率与动态断裂韧度关系,是根据推出的裂纹速率及应变率关系,与裂纹速率及断裂韧度关系相结合而得到。研究了应变率对应力-应变本构关系及动态压缩强度影响。并通过试验结果验证了模型的合理性。讨论了岩石初始损伤、围压、模型中参数m、ε0和R对应力-应变关系、动态压缩强度和动态弹性模量的影响。研究结果可为动态压缩荷载作用下深部地下工程脆性围岩稳定性分析提供了一定的理论支持。     

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.     

Effect of microstructure on fracture of brittle materials: Unified approach

A theoretical approach to the fracture of brittle solids based on crack opening displacement and energy rate criterion is presented. The approach allows for the prediction of elastic (Young’s modulus) and fracture (fracture strength and thermal shock) response of a brittle material containing spherical pores and polycrystalline solids containing anisotropic residual stresses.     

A cohesive plastic and damage zone model for dynamic crack growth in rate-dependent materials

Mode I steady-state dynamic crack growth in rate-dependent viscoplastic solids containing damage, under small scale yielding conditions, is analyzed based on a modified cohesive zone model. A multi-scale approach is used to describe the entire non-linear zone consisting of a plastic region and a damage region, each of which has its own constitutive law. Traction in the damage region is characterized by a softening power-law, in terms of the ultimate strength, a softening index and a rate sensitivity factor. In the plastic region, the cohesive law is assumed to be both strain hardening and rate dependent. The critical crack opening displacement at the physical crack-tip controls crack growth. The governing integral equations are derived and solved by a collocation method combined with associated boundary conditions. Numerical results are presented for the traction and opening profiles along the cohesive zone, the fracture energy and lengths of the damage and non-linear zones at different crack speeds and for different material parameters. The importance of factors, such as material softening, plastic deformation, crack speed and viscosity, is identified by parametric studies. In addition, the competition of plastic flow and material damage, and its effect on crack growth, are discussed.     

非贯通裂隙岩体损伤演化率相关性及变形特征

含非贯通裂隙岩体是自然界中岩体的主要赋存形式,其裂隙几何特征对岩体的强度及变形均产生显著影响。应变率对岩体的损伤演化及黏滞效应也具有显著的率相关性。首先,运用模型元件的方法,将非贯通裂隙岩体动态破坏过程视为具复合损伤、静态弹性特性、动态黏滞特性的非均质点组成,对黏弹性响应的Maxwell体进行改进,将细观损伤体与裂隙损伤演化的宏观损伤体根据等效应变假设并联组成宏细观复合损伤体,构建综合考虑岩体宏细观缺陷的动态损伤模型;其次,基于断裂力学及应变能理论,对岩体宏观裂隙动态扩展的能量机制进行分析,综合考虑初始裂隙应变能、裂隙动态损伤演化过程应变能、裂隙闭合应变能,得到裂隙岩体宏观动态损伤变量计算公式;最后,将模型计算结果与实验结果进行比较,模型计算结果与实验结果吻合较好,证明了模型的合理性,同时利用模型讨论了裂隙倾角、应变率、岩石性质对岩体变形特征的影响规律。     

Modelling of dynamic ductile fracture and application to the simulation of plate impact tests on tantalum

A model of dynamic damage by void nucleation and growth is proposed for elastic-viscoplastic materials sustaining intense loading. The model is dedicated to ductile materials for which fracture is caused by microvoiding. The material contains potential nucleation sites where microvoids are generated when the local pressure overcomes the nucleation pressure. A probability density function is adopted to describe the fluctuation of the nucleation pressure within the material. The void growth is described by using a hollow sphere model where micro-inertia effects are accounted for. The matrix weakening due to void growth is also included.The model has been first tested under uniaxial deformation. When the strain rate is assumed constant, the pressure inside the material has nearly a linear response up to a maximum. An analytical expression for the maximum pressure is proposed.Finite element simulations of plate impact tests have been carried out and compared to experiments on tantalum. From simulations based on the proposed model, an increase of the spall strength is observed with higher shock intensities. Therefore, the relationship between the velocity pullback and spall strength usually assumed in the literature (based on the acoustic approach) seems to be inadequate. Velocity profiles are simulated for different flyer thicknesses and different impact velocities with close agreement with experiments.     

Effects of microscale inertia on dynamic ductile crack growth

The aim of this paper is to investigate the role of microscale inertia in dynamic ductile crack growth. A constitutive model for porous solids that accounts for dynamic effects due to void growth is proposed. The model has been implemented in a finite element code and simulations of crack growth in a notched bar and in an edge cracked specimen have been performed. Results are compared to predictions obtained via the Gurson–Tvergaard–Needleman (GTN) model where micro-inertia effects are not accounted for. It is found that microscale inertia has a significant influence on the crack growth. In particular, it is shown that micro-inertia plays an important role during the strain localisation process by impeding void growth. Therefore, the resulting damage accumulation occurs in a more progressive manner. For this reason, simulations based on the proposed modelling exhibit much less mesh sensitivity than those based on the viscoplastic GTN model. Microscale inertia is also found to lead to lower crack speeds. Effects of micro-inertia on fracture toughness are evaluated.     

Photoelasto-plastic studies on brittle fracture of high-polymer solids

The mechanism of brittle fracture of high-polymer solids is experimentally investigated under one-or two-dimensional stress states by a new photoelastoplastic method suggested by the author. The application of the photoelasto-plastic method on the brittle-fracture problem is based on the principle that breaking stress can be computed in brittle fracture by the measurement of the fringe orderN _B of isochromatic lines at fracture point. Bending under three-point and four-point loads, and the plane problems, some having stress concentration and others being under contacting load, are examined by using rigid polyester cast resin containing styrol as a model specimen; and, in conclusion, the brittle fracture of high-polymer solids under one- or two-dimensional stress states is decided by the constant tensile stress, whose magnitude depends only upon the material used as a model specimen, and is larger than its ultimate tensile strength. Many kinds of factors in fracture are defined, and stress-concentration factors in fracture are compared with stress-concentration factors in elasticity. A new photoelasto-plastic simple method for the determination of stress-concentration factors in elasticity is suggested by utilization of the experimental results on this brittle fracture of high-polymer solids and is examined on the perforated plane problem having finite width under tension in comparison with theoretical analysis and the experimental results by other measuring methods.     

Zr基非晶合金JH-2模型的构建及应用

为构建Zr_62.5Nb₃Cu_14.5Ni₁₄Al₆非晶合金在高压、大应变、高应变率状态下的材料模型,采用根据实验数据理论推导和数值模拟对比反馈的方法,对材料的Johnson-Holmquist本构模型（JH-2模型）参数进行了研究:材料的静水压力-体应变关系通过平板冲击实验数据和理论推导得到;无损材料强度与应变、应变率的关系通过轴向压缩实验数据确定;材料损伤参数与破碎材料强度参数的关系通过平板冲击实验数据确定;破碎材料强度参数通过数值模拟与实验结果对比的反馈法得到。将材料模型应用于平板冲击和破片侵彻的数值模拟,通过数值模拟与实验结果对比的方式,验证材料模型的准确性。结果表明,平板冲击实验中,材料的自由面粒子速度曲线与数值模拟结果吻合度较高;破片侵彻实验中,破片对钢靶的侵彻深度、开坑孔径与数值模拟结果的一致性较好,构建的材料模型较准确反映了材料的动态力学特性。     

考虑岩石微元强度呈对数正态分布的冻融岩石损伤特性

冻融环境下准脆性岩石损伤力学特性对寒区工程建设有重要影响。由于准脆性岩石存在特征长度,其尺度效应不再符合Weibull尺度效应统计理论。本文以准脆性岩石为研究对象,考虑准脆性岩石微单元强度服从对数正态分布的条件下,依据D-P破坏准则,运用损伤力学理论建立了围压作用下冻融准脆性岩石的损伤本构模型,对冻融准脆性岩石损伤力学特性进行了探讨。研究表明,（1）采用对数正态分布建立的冻融准脆性岩石损伤力学本构模型与试验结果一致,该模型反映出冻融环境下准脆性岩石总损伤沿应变、围压和冻融循环路径相互作用,能够揭示出准脆性岩石损伤演化的本质特征。（2）冻融与围压共同作用下,准脆性岩石应力-应变变化关系可划分为线弹性阶段、非线性强化阶段、应力减弱阶段以及应变衰退阶段。（3）当准脆性岩石所受围压一定时,冻融次数是使准脆性岩石冻融受荷总损伤增大的关键因素。冻融循环次数对准脆性岩石的损伤劣化起促进作用;当准脆性岩石所受冻融循环次数不变时,冻融受荷总损伤随围压增大而减小,即围压对准脆性岩石力学性能劣化起抑制作用。相关研究成果对寒区岩体工程安全稳定性评价有一定的理论参考价值。     

Modeling dynamic brittle behavior of materials with circular flaws or pores

Compressive failure of brittle materials is driven primarily by crack growth from pre-existing flaws in the material. These flaws, such as grain boundaries, pores, preexisting cracks, inclusions and missing grains, are randomly spaced and have a range of possible shapes and sizes. The current work proposes a micromechanics-based model for compressive dynamic failure of brittle materials with circular pore flaws, which incorporates both the number density and the size distribution of flaws. Results show that the distribution of flaw sizes is very important, particularly at moderate strain rate, since analyses based solely on the mean flaw size overpredict strength. Therefore, in order to increase dynamic strength at low to moderate strain rates, it is most effective to control the presence of large flaws. At very high strain rates, however, crack growth is activated even in small flaws and therefore controlling the total number density rather than the size of the flaws is effective for increasing dynamic strength. Finally, the model shows that neglecting very small flaws in the pore population may not have significant effects on the results in many cases, suggesting that the model is a useful tool for identifying a minimum resolution required for experimental characterization of microstructure.     

PMMA膨胀环动态拉伸碎裂实验研究

在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150\mathrm{~}500s^{-1}范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.      

考虑微裂纹相互作用的的岩石微观力学弹塑性损伤模型研究

本文建立基于微裂纹扩展的岩石弹塑性损伤微观力学模型。用自洽方法考虑裂隙间相互影响,压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的微观损伤,基于应变能密度准则用Newton迭代法求复合型断裂的翼裂纹扩展长度,并采用微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映塑性内变量对硬化函数的影响;建立岩石的弹塑性损伤本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性损伤模型的程序。从围压和微裂隙长度等因素分析弹塑性损伤模型的岩石的损伤和宏观塑性特性。     

A rate-dependent damage model for brittle materials based on the dominant crack

A rate-dependent, continuum damage model is developed for brittle materials under dynamic loading. This model improves on the approach (ISOSCM) of [Addessio, F.L., Johnson, J.N., 1990. A constitutive model for the dynamic response of brittle materials. Journal of Applied Physics 67, 3275–3286] in several respects. (1) A new damage surface is found by applying the generalized Griffith instability criterion to the dominant crack (having the most unstable orientation), rather than by averaging the instability condition over all crack orientations as done previously. The new surface removes a discontinuity in the damage surface in ISOSCM when the pressure changes sign. (2) The strain due to crack opening is more consistent with crack mechanics, with only the tensile principal stresses contributing to the crack opening strain. This is achieved by incorporating a projection operator in the equation for the crack opening strain. One consequence of incorporating the projection operator is a prediction of shear dilatancy, which is not accounted for in ISOSCM. (3) The evolution of damage, which is based on the energy-release rate for the dominant crack, has a physical basis, whereas in the previous approach the damage growth rate was assumed to be an exponential function of the distance from the stress state to the damage surface without specific physical justification.An implicit algorithm has been developed so that a larger time step can be used than with the explicit algorithm used in ISOSCM. The numerical results of a silicon carbide (SiC) ceramic under several loading paths (hydrostatic tension/compression, uniaxial strain, uniaxial stress, and shear) and strain rates are presented to illustrate the main features of the model.     

Localization of deformation and stress–strain relation for mesoscopic heterogeneous brittle rock materials under unloading

Stress redistribution induced by excavation of underground engineering and slope engineering results in the unloading zone in parts of surrounding rock masses. The mechanical behaviors of crack-weakened rock masses under unloading are different from those of crack-weakened rock masses under loading. A micromechanics-based model has been proposed for brittle rock material undergoing irreversible changes of their microscopic structures due to microcrack growth when axial stress is held constant while lateral confinement is reduced. The basic idea of the present model is to classify the constitution relation of rock material 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 closed-form explicit expression for the complete stress–strain relation of rock materials containing cracks under unloading is obtained. The results show that the complete stress–strain relation and the strength of rock materials under unloading depend on the crack spacing, the fracture toughness of rock materials, orientation of the cracks, the crack half-length and the crack density parameter.     

Modelling of compressive behaviour of concrete-like materials at high strain rate

Uniaxial compression tests are the most common tests for characterizing the strength of concrete-like materials. The dynamic compression strength of concrete-like material is typically obtained by Split Hopkinson Pressure Bar (SHPB) tests. The increase in material strength under dynamic loading is usually attributed to the strain rate effect and modelled with a dynamic increase factor (DIF). However, it was observed by some researchers that the radial inertial confinement caused apparent increase of dynamic strength of concrete-like specimen in SHPB tests. They attributed the material strength increase to this inertial effect, instead of the strain rate effect. In the present study, numerical analyses are performed to investigate the compressive behaviour of concrete-like material at high strain rates. A homogeneous macroscale model and a heterogeneous mesoscale model are developed in the study. In the macroscale model, the material is assumed to be homogeneous and isotropic. In the mesoscale model, the test sample is modelled as a three-phase composite consisting of aggregate, mortar matrix and interfacial transaction zone (ITZ) between the aggregate and the mortar matrix. The aggregate is assumed to be circular and the ITZ is modelled as a thin boundary around the aggregate. In the both models, the materials are assumed to be insensitive to the strain rate first. Therefore, the obtained strength enhancement is only due to the inertial confinement. Strain rate sensitive material properties are then used in the two models in the calculations. Numerical simulations of the concrete samples under compression at different strain rates are carried out. The relative contribution of the inertial effect and the strain rate effect on the compressive strength DIF is examined based on the numerical results. The failure process of concrete specimen is also studied.     

动力损伤后的脆性岩石静力蠕变断裂模型研究

应力波动力扰动下脆性岩石的静力蠕变特性,对深部地下工程围岩变形的评价有重要的实践意义.动力载荷作用导致的局部细观裂纹损伤严重影响脆性岩石蠕变力学行为.基于细观裂纹扩展与应力关系模型、动力扰动损伤演化函数、静动力载荷演化路径函数与黏弹性本构模型,提出一种应力波动力扰动下脆性岩石蠕变断裂特性的宏细观力学模型.其中动力损伤通过控制岩石内部细观裂纹数量变化实现.模型描述了应力波动力扰动下岩石的应变时间演化曲线,解释了岩石动力扰动下蠕变失效特性.研究了不同应力波幅值及周期影响下的脆性岩石应变-时间关系曲线,并通过试验结果验证了模型的合理性.讨论了动力损伤变化形式,突变发生时刻,突变量的大小对岩石蠕变失效特性的影响.分析了应力波幅值、周期对岩石动态动力损伤效应以及蠕变失效特性的影响.主要研究结果:动力损伤的变化值越大,岩石蠕变失效发生时间越短.冲击载荷扰动期间,动力损伤发生的时刻及增加的形式,对动力扰动后的岩石应变及蠕变破坏时间影响很小.动力损伤变化量随应力波幅值增加、周期减小而加速增大.应力波幅值越大、周期越小,岩石发生蠕变失效时间越短.     

混凝土类材料动态压缩强度在多维应力状态下的应变率效应

对混凝土类材料动态压缩应变率效应研究的发展及问题进行了概述,对比不同应力状态下混凝土类材料动态压缩应变率效应的表现特征,揭示了不同加载路径下实测动态强度提高系数的显著差异。研究表明,在高应变率下,基于初始一维应力加载路径的试件将因横向惯性效应导致的侧向围压而演化至多维应力状态,传统霍普金森杆技术无法获得高应变率下基于真实一维应力路径的动态强度提高系数,在强度模型中直接应用实测数据将过高估计材料的动态强度。鉴于应变率效应的加载路径依赖性,将仅包含应变率的强度提高系数模型扩展至同时计及应变率和应力状态的多维应力状态模型,并结合Drucker-Prager准则在强度模型中给予了实现。针对具有自由和约束边界试件开展的数值霍普金森杆实验表明,多维应力状态下的应变率效应模型可以考虑应变率效应随应力状态改变的特点,从而准确预测该类材料的动态压缩强度。研究结果可为正确应用霍普金森杆技术确定脆性材料的动态压缩强度提供参考。     

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

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