基于最小耗能原理的黏弹塑性损伤模型

在基于“应变等效假定”和“余应变能等效假定”两种不同损伤力学基本假定的黏塑性损伤本构方程基础上,将最小耗能原理应用于固体介质的黏弹塑性问题,给出了损伤材料的黏塑性应变流变速率、黏塑性损伤发展速率和黏弹-塑性材料积累应变强弱化速率的非线性演化发展方程,建立了两个不同的黏弹塑性损伤模型（A和B）,并进行了试验验证。作为算例,最后运用黏塑性损伤模型A,对地震荷载作用下建于岩石基础上的某混凝土重力坝的动力损伤问题进行了分析,结果表明该模型应用于混凝土重力坝的动力损伤分析是可行的。     

基于内变量和张量函数表示定理的本构方程

针对各向同性材料,基于张量函数表示定理,建立了本构关系的张量不变性表示,其中,3个不可约基张量取决于应力的0～2次幂,且相互正交,3个系数由塑性应变增量和应力的不变量表示。基于塑性应变增量的不变量定义内变量,本构关系归结为确定内变量的演化。使用张量函数表示定理,给出了内变量演化方程的一般表达式,它取决于应力不变量的增量,因而与主轴旋转无关。讨论了如何根据试验资料和引入适当的假定,确定具体的演化方程。通过与塑性势理论和多重屈服面理论进行比较,表明所建模型是这些理论的最一般表示,且简捷直观、使用方便。     

广义Drucker-Prager强度准则

基于摩擦性材料的试验研究,在已有研究成果的基础上,提出了具有Drucker-Prager准则表达形式的广义非线性强度理论,即广义Drucker-Prager强度理论。该强度理论在偏平面上的强度函数是介于Drucker-Prager准则和Mohr-Coulomb准则外包线的光滑外凸曲线,在子午面上的强度函数为线性函数。该强度理论用一个表达式来统一描述材料的强度特性,包含或逼近了现有主要的非线性单一强度理论,并且在应力空间中偏应力和静水压力分离,容易与具体的应力应变模型相结合。     

Thermodynamic‐based model for coupling temperature‐dependent viscoelastic,viscoplastic, and viscodamage constitutive behavior of asphalt mixtures

Based on the continuum damage mechanics, a general and comprehensive thermodynamic‐based framework for coupling the temperature‐dependent viscoelastic, viscoplastic, and viscodamage behaviors of bituminous materials is presented. This general framework derives systematically Schapery‐type nonlinear viscoelasticity, Perzyna‐type viscoplasticity, and a viscodamage model analogous to the Perzyna‐type viscoplasticity. The resulting constitutive equations are implemented in the well‐known finite element code Abaqus via the user material subroutine UMAT. A systematic procedure for identifying the model parameters is discussed. Finally, the model is validated by comparing the model predictions with a comprehensive set of experimental data on hot mix asphalt that include creep‐recovery, creep, uniaxial constant strain rate, and repeated creep‐recovery tests in both tension and compression over a range of temperatures, stress levels, and strain rates. Comparisons between model predictions and experimental measurements show that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

岩石弹脆性分维损伤本构模型

本文定义岩石构元中破裂面的分维值为各向同性损伤变量,而各个方向上裂纹面的累加量定义为各向异性损伤变量,并根据裂纹发育特征提出了损伤变量演化方程,从而建立起岩石脆性变形破坏过程的分维损伤本构模型。最后,利用该模型对大理岩单轴压缩应力应变曲线进行了模拟,结果说明本文提出的模型是较为合理的。     

A method of time-dependent analysis using elastic solutions for nonlinear materials

The formulation of viscoelastic solutions from elastic equations using the ‘correspondence principle’ and an inverse Laplace transform has been discussed extensively in the literature. Because this method has been developed, many time-dependent solutions can be obtained from closed form elastic solutions and conditions have been delineated in which the ‘quasi-elastic’ approximation of the viscoelastic solution is within acceptable tolerance. This communication shows the feasibility of the application of these methods to formulate approximate nonlinear viscoelastic solutions with nonlinear stress-strain materials, and for want of a specific nonlinear model to demonstrate this, the hyperbolic model was selected. The ‘power law’ is used to model the relaxation modulus of the viscoelastic materials. There are five related development that are discussed here using a simple numerical example to illustrate each of them and they are: (1) a linear elastic solution, (2) a linear viscoelastic solution, (3) a nonlinear elastic solution, (4) a nonlinear viscoelastic solution and finally, (5) a ‘regression’ approximation of the nonlinear viscoelastic solution which is suggested by the series form of the elastic solution. All of these are related to one another and each provides an acceptably accurate solution of the problem it addresses. The latter is of particular practical interest since it can be used to provide answers to problems involving nonlinear viscoelastic materials while requiring only very small calculation times. The problem used as an example is the calculation of the displacement of a circular hole in an infinite plate made of a material with a nonlinear time-dependent stress-strain relationship. The nonlinear elastic form of the solution was developed by matching results from nonlinear finite element analysis.     

基于内时理论的含瓦斯煤岩损伤本构模型研究

从内时理论出发,通过内蕴时间的重新构造,在Helmholtz自由能中引入损伤变量,利用连续介质不可逆热力学的基本原理推导出了含瓦斯煤岩的内时损伤本构方程,并给出了确定本构模型中各参数的方法。同时,对在三轴应力条件下所建立的本构方程进行了求解分析与试验验证。结果表明,该本构方程能很好地描述含瓦斯煤岩在损伤发生前后的变形特点,并能有效地反映含瓦斯煤岩的非线性和剪胀扩容等物理力学现象。     

基于不同弹性本构方程的塑性区宽度及封孔长度研究

煤层注水是煤矿开采过程中重要的增透卸压手段,而封孔长度是影响注水压裂效果的一个关键因素。为了确定出合理的封孔长度,根据Lippmann煤层失稳基本理论,选取了6种不同的弹性区本构方程,引入非线性动态瓦斯压力分布的表达式,构建了含瓦斯煤层力学理论模型,得到煤层垂直、水平应力和塑性区宽度的解析式,之后结合高家堡的现场情况,对比分析6种弹性本构方程下的煤层应力分布曲线及塑性区宽度,并进行了现场封孔试验。结果表明：封孔长度在9.5 m时孔洞漏水量小、未出现异常现象,封孔效果良好,最终确定出了合理的封孔长度为9.5 m;相对于其他3种弹性本构方程,选取第i、v、vi种弹性本构方程对封孔长度进行理论计算分析,更加符合现场实际情况。     

混凝土混合硬化黏塑性损伤统本构模型

传统的黏塑性统一本构模型只能用于金属类材料的本构分析,而不适用于混凝土的本构分析.基于此,以不可逆热力学理论为基础,选定了涉及混合硬化参量的Helmholtz自由能函数,推导出混凝土的黏塑性损伤耦合本构模型.通过对混凝土硬化模型的分析,构造了含运动硬化和等向硬化内变量的非相关流动势函数,又推导出流动方程以及内变量演化方程.所建模型对混凝土试验曲线的数值模拟显示,其能够正确描述混凝土的率相关性、非弹性体积膨胀特性、加载过程的软硬化特性以及由断裂和损伤引起的应力软化现象.     

基于堆石坝竣工期沉降观测数据的材料非线性本构模型参数反演

为了高效和准确地确定堆石料的非线性本构模型参数,提出了基于响应面方法的参数反演方法。采用有限元方法数值模拟了堆石坝的分层填筑过程。建立了堆石坝变形观测点沉降的响应面函数,确定了多项式响应面函数的系数。根据堆石坝竣工期变形观测数据和确定的响应面函数,采用优化方法反演确定了堆石料本构模型参数。工程实际应用结果表明,该方法具有较高的计算效率,预测的堆石坝沉降变形与现场观测值吻合较好。     

基于改进Bingham模型的软岩参数非定常三维非线性黏弹塑性蠕变本构研究

传统线性元件组合模型难以描述岩石非线性加速蠕变阶段特性,给工程应用带来了巨大困难。提出一种新的应变触发非线性黏壶元件,并将其与Bingham模型串联,建立了可以描述岩石等速和加速蠕变特性的一维本构模型,进而推导出相应的三维蠕变本构模型,再将三轴压缩过程中岩石弹性模量的衰减方程引入该三维蠕变本构方程,得到一个能反映岩石蠕变全过程特性的三维非线性黏弹塑蠕变本构模型,根据泥岩蠕变试验数据对该模型参数进行了辨识。模型理论计算值与泥岩蠕变试验数据的对比表明改进的Bingham模型不仅可以充分反映软岩的等速蠕变特性,还可以很好地描述软岩的衰减和加速蠕变规律,且模型元件少,组合形式简单,为软岩非线性蠕变本构模型研究提供了新的借鉴。     

考虑温度效应的岩石损伤内时本构关系

本文根据岩石类材料受载时的损伤破坏机理,首次引入两个损伤变量 D_1和 D_0用以反映外载荷体量和偏量的作用,在内时理论的基础上建立起考虑温度效应的岩石损伤本构方程,并导出具有明确物理意义的本构方程增量型关系式。在不同温度下用大理岩试样进行了压缩实验和声发射实验。根据实验所得结果,确定了本构方程中材料的所有参数,得到了大理岩材料受温度影响其性质发生变化的一些有意义的结论。     

Incrementalization of a single hardening constitutive model for frictional materials

The governing equations for an elasto‐plastic constitutive model for frictional materials such as soil, rock, and concrete are presented, and the incremental form is indicated in preparation for implementation of the model in a user‐defined module for finite element calculations. This isotropic, work‐hardening and ‐softening model employs a single yield surface, it incorporates non‐associated plastic flow, and its capability of capturing the behaviour of different types of frictional materials under various three‐dimensional conditions has been demonstrated by comparison with measured behaviour, as presented in the literature. The incrementalization procedure is indicated and the resulting equations for the single hardening model are presented together with parameters for a dense sand. Following the implementation of the model, these parameters are used for evaluation of different integration schemes as presented in a companion paper by Jakobsen and Lade (Int. J. Numer. Anal. Meth. Geomech. 2002; 26 :661). Copyright © 2002 John Wiley & Sons, Ltd.     

复杂应力路径下堆石体本构模型比较验证

目前对堆石体应力-应变关系特性的研究主要建立在常规三轴试验的基础上,对复杂应力路径上堆石体本构模型的验证工作尚不充分。利用糯扎渡高心墙堆石坝主堆石料复杂应力路径大型三轴试验成果,对国内常用的堆石体本构模型--邓肯-张EB模型、清华非线性解耦KG模型和沈珠江双屈服面模型进行了比较和验证。结果表明,邓肯-张EB模型通常会夸大堆石体的体积压缩变形,其加卸载准则无法正确判别一些复杂应力路径的加卸载状况,清华KG模型和沈珠江双屈服面模型对堆石体复杂应力路径的适应性相对较好。     

基于应变能密度理论的岩石破裂数值模拟方法研究

采用双线性应变软化本构模型结合能量耗散原理建立了损伤本构方程,并通过应变能密度理论建立了细观单元岩石破坏的能量判别准则。当某一单元所存储的应变能超过某固定值时,单元进入损伤状态,同时单元的损伤程度随着能量耗散的增加而增加,损伤单元的材料属性也随之改变,直到变为具有一定残余强度的单元。随着荷载增加,单元损伤的程度变大,当单元储存的应变能超过所建立的能量判别准则时,定义单元破坏,随着破坏单元的数目不断增多,破坏单元相互连通形成宏观裂纹,实现了利用线性计算完成非线性计算的过程,避免了数值计算在单元断裂时的奇异性,模拟了岩石的峰后破裂行为。对上述算法利用FLAC中的FISH语言开发了岩石破裂化计算程序,并将该程序成功应用于巴西劈裂和中间裂隙拉伸试验的破裂模拟过程中,其模拟结果与相应的理论和试验结果吻合较好,说明该方法对于模拟岩石破裂过程的正确性和可行性。     

A stability criterion for elasto‐viscoplastic constitutive relationships

In this paper, the onset of mechanical instability in time‐sensitive elasto‐viscoplastic solids is theoretically analyzed at the constitutive level and associated with the occurrence of ‘spontaneous accelerations’ under stationary external perturbations. For this purpose, a second‐order form of Perzyna's constitutive equations is first derived by time differentiation, and a sufficient stability condition is identified for general mixed loading programs. These loading conditions are in fact the most general in both laboratory tests and real boundary value problems, where a combination of certain stress and strain components is known/prescribed. The theoretical analysis leads to find precise stability limits in terms of material hardening modulus. In the case of constitutive relationships with isotropic strain‐hardening, no instabilities are possible while the hardening modulus is larger than the so‐called ‘controllability modulus’ defined for (inviscid) elasto‐plastic materials. It is also shown that the current stress/strain rate may also directly influence the occurrence of elasto‐viscoplastic instability, which is at variance with elasto‐plastic inviscid media. Copyright © 2015 John Wiley & Sons, Ltd.     

Instability and ill-posedness in the deformation of granular materials

A linear stability analysis of the partial differential equations of granular flow is performed. The constitutive relations include (i) elastic effects, (ii) non-associated flow rules and (iii) shear-strain hardening. Conditions for the equations to be well-posed and to be stable are derived. It turns out that there are two qualitatively different kinds of instability; which appears first depends on the parameters of the constitutive relations. If one kind appears first, then in a constitutive test a period of inhomogeneous deformation is predicted to occur before the formation of shear bands. If the other kind appears first, then shear bands are predicted to form in an approximately homogeneous sample. Some constitutive experiments are analysed from this viewpoint, and the analysis suggests some new experimental work.     

Isogeometric analysis of THM coupled processes in ground freezing

An isogeometric analysis (IGA) based numerical model is presented for simulation of thermo-hydro-mechanically (THM) coupled processes in ground freezing. The momentum, mass and energy conservation equations are derived based on porous media theory. The governing equations are supplemented by a saturation curve, a hydraulic conductivity model and constitutive equations. Variational and Galerkin formulation results in a highly nonlinear system of equations, which are solved using Newton-Raphson iteration. Numerical examples on isothermal consolidation in plane strain, one-dimensional freezing and heave due to a chilled pipeline are presented. Reasonably good agreements were observed between the IGA based heave simulations and experimental results.     

岩石类材料的能量基率相关弹塑性损伤模型

岩石类材料的动态本构模型研究是岩石动力学理论研究的基础问题。相比于其他的非线性理论,损伤力学理论已证实可成功地模拟岩石类材料的应变软化和渐进破坏等特征,可用于解释其静、动态破坏机制。现有的通过理论推导得出的动态本构中大部分并未考虑损伤因素,而通过维象学试验方法建立起的动态本构则缺乏损伤力学理论基础。为此基于弹性余能等效原理和损伤力学的基本概念,并结合动态试验,推导建立了岩石类材料的率相关弹塑性损伤模型。并通过与文献中的试验资料对比,证明了模型的有效性,为岩石类材料组成的结构的动力响应研究奠定了基础。