各向异性弹黏塑性模型在ABAQUS中的研发

通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明：(1) 对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2) 对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。     

Manuel Rocha Medal Recipient Simulating the Time-dependent Behaviour of Excavations in Hard Rock

Summary Although hard rock is not usually associated with large creep deformation, data collected from the tunnels and stopes of the deep South African gold mines illustrates significant time-dependent behaviour. Apart from application in mining, a better understanding of the time-dependent behaviour of crystalline rock is required to analyse the long term stability of nuclear waste repositories and to design better support for deep civil engineering tunnels in these rock types. To illustrate the subtle problems associated with using viscoelastic theory to simulate the time-dependent behaviour of hard rock, a viscoelastic convergence solution for the incremental enlargement of a tabular excavation is discussed. Data on the time-dependent deformation of a tunnel developed in hard rock further illustrates the limitations of the theory, as it is unable to simulate the fracture zone around these excavations. To simulate the rheology of the fracture zone, a continuum viscoplastic approach was developed and implemented in a finite difference code. This proved more successful in modelling the time-dependent closure of stopes and squeezing conditions in hard rock tunnels. A continuum approach, however, has limitations in areas where the squeezing behaviour is dominated by the time-dependent behaviour of prominent discontinuities such as bedding planes. To overcome this problem, a viscoplastic displacement discontinuity technique was developed. This, combined with a tessellation approach, leads to more realistic modelling of the time-dependent behaviour of the fracture zone around excavations. Received January 15, 2002; accepted June 3, 2002 Published online September 2, 2002     

Visco-Plastic Behaviour around Advancing Tunnels in Squeezing Rock

Summary  The visco-plastic behaviour of rocks plays a relevant role in the tunnelling works, especially for deep tunnels subjected to large initial stresses for which squeezing conditions may develop. A rheological model is discussed that accounts for visco-elastic (primary) and visco-plastic (secondary) contributions to rock creep. The effects of tertiary creep are included in the model by way of a gradual mechanical damage governed by the cumulated visco-plastic strains. The parameters of the intact rock are first identified based on laboratory test results presented in the literature. Then, after scaling them to those of the rock mass, the potential applicability of the model is tested through axisymmetric and plane strain finite element analyses of the full face excavation of a deep circular tunnel. The results are discussed with particular reference to the short term redistribution of stresses around the opening and to its convergence. The analyses show the relevant influence of tertiary creep on the tunnel closure. In addition, those based on an axisymmetric scheme turn out to be crucial for the correct long term prediction of the interaction between the rock mass and the supporting structure of the opening.     

Triaxial behaviour of transversely isotropic materials: application to sedimentary rocks

Failure and long‐term behaviour of oriented solids are studied. Transversely isotropic materials are considered and a mathematical formulation that respect the material symmetry is developed and applied to model the triaxial behaviour of sedimentary rocks. Two failure criteria and a viscoplastic constitutive model that describe, respectively, triaxial failure and triaxial creep tests are presented and discussed. The application of the developed models to describe the mechanical behaviour of Tournemire shale shows that theoretical predictions are in good agreement with the experimental data. In the present paper, the developed approach is applied to sedimentary rock materials, nevertheless, it can be generalized to any material that exhibits transverse isotropy. Copyright © 2007 John Wiley & Sons, Ltd.     

A unified thermoplastic/viscoplastic constitutive model for geomaterials

Research on the effect of temperature changes on the behaviour of geomaterials has become increasingly important in recent years. This growing interest is partially due to the recent development of high-level nuclear waste disposals. Because of the complex influence of temperature in these areas, it is necessary to understand the effects of temperature on rock-like materials and use the appropriate constitutive equations to numerically model these phenomena. In this paper, a thermoplastic/viscoplastic constitutive model is developed for this purpose. The model includes thermal softening, the evolution of the yield functions with temperature, and the effects of temperature on the time-dependent behaviour. The model performance is demonstrated by some simple test cases on Tournemire and Bure clayey rocks including triaxial compression tests and creep tests under constant temperatures. The numerical results are discussed using experimental data, which demonstrate that the model can reproduce the overall behaviour of this type of materials under deviatoric loads and non-isothermal conditions.     

Strain localization analysis by elasto-viscoplastic softening model for frozen sand

The aim of the present paper is to numerically analyse the behaviour of frozen sand by using a viscoplastic constitutive model with strain softening. A constitutive model has been developed introducing the stress history tensor which is a functional of the stress history, with respect to a generalized time measure. It is shown that Adachi and Oka's model is applicable to the results of triaxial tests on a frozen Toyoura sand at different strain rates. First, the instability of the model is discussed within the framework of bifurcation theory. The model is then implemented into a FEM code to numerically simulate the behaviour under plane strain conditions. From the numerical results, it is revealed that the formation of shear bands is possible and the characteristics of strain localization, such as shear banding, depend on the strain rates.     

A viscoplastic model including anisotropic damage for the time dependent behaviour of rock

This paper presents a new constitutive model for the time dependent mechanical behaviour of rock which takes into account both viscoplastic behaviour and evolution of damage with respect to time. This model is built by associating a viscoplastic constitutive law to the damage theory. The main characteristics of this model are the account of a viscoplastic volumetric strain (i.e. contractancy and dilatancy) as well as the anisotropy of damage. The latter is described by a second rank tensor. Using this model, it is possible to predict delayed rupture by determining time to failure, in creep tests for example. The identification of the model parameters is based on experiments such as creep tests, relaxation tests and quasi‐static tests. The physical meaning of these parameters is discussed and comparisons with lab tests are presented. The ability of the model to reproduce the delayed failure observed in tertiary creep is demonstrated as well as the sensitivity of the mechanical response to the rate of loading. The model could be used to simulate the evolution of the excavated damage zone around underground openings. Copyright © 2005 John Wiley & Sons, Ltd.     

A viscoplastic subloading soil model for rate‐dependent cyclic anisotropic structured behaviour

This paper presents a new purely viscoplastic soil model based on the subloading surface concept with a mobile centre of homothety, enabling the occurrence of viscoplastic strains inside the yield surface and avoiding the abrupt change in stiffness of the traditional overstress viscoplastic models. This is required for overconsolidated soils. The model is formulated to reproduce the soil rate‐dependent behaviour under cyclic loading (changes in loading direction) and incorporates both initial and induced anisotropy, as well as destructuring. The model shows good qualitative response to some imposed three‐dimensional stress paths under quasi‐inviscid (elastoplastic) behaviour. Some of the main time‐dependent aspects of soil behaviour that the model is capable of reproducing were also illustrated. The capability of the model to adequately reproduce the results from an undrained triaxial test performed on stiff overconsolidated clays from the Lisbon region (Formação de Benfica), with an unloading–reloading deviatoric stress cycle at constant mean stress, that incorporates a series of staggered fast loading and creep stages, was evaluated. The model was able to reproduce well the main observed aspects of the time‐dependent stress–strain response and pore pressure evolution of a stiff overconsolidated clay under complex loading. The revised and generalised viscoplastic subloading surface concept is viable and can be applied to a consistent extension to viscoplasticity, including in the interior of the yield surface, of existing elastoplastic models formulated for soils and other materials. Copyright © 2016 John Wiley & Sons, Ltd.     

A constitutive model for crushed salt

A constitutive model for crushed salt is presented in this paper. A creep constitutive model is developed first and compared with test results. The constitutive model presented here concentrates on creep deformation because saline media behave basically in a ductile and time‐dependent way. An idealized geometry is used as a common framework to obtain stress–strain macroscopic laws based on two deformation mechanisms: fluid‐assisted diffusional transfer creep and dislocation creep. The model is able to predict strain rates that compare well with results from laboratory tests under isotropic and oedometric conditions. Macroscopic laws are written using a non‐linear viscous approach, which incorporates also a viscoplastic component, based on critical state theory. The viscoplastic term is intended for non‐creep deformation mechanisms such as grain reorganization and crushing. Copyright © 2002 John Wiley & Sons, Ltd.     

Undrained behaviour of Changi sand in triaxial and plane-strain compression

The results of an experimental study of the undrained behaviour of Changi sand under axisymmetric and plane-strain conditions are presented. K₀ consolidated undrained plane-strain tests and K₀ or isotropically consolidated triaxial tests on very loose and medium dense specimens were conducted. The undrained behaviour of sand at very loose and medium dense states under plane-strain conditions was characterised and compared with that under axisymmetric conditions. It was observed that the undrained behaviour of very loose and medium dense sand under plane strain is similar to that under axisymmetric conditions. However, because of the formation of shear bands in plane-strain tests, the post-peak behaviour of medium dense sand in plane strain is different from that in triaxial tests. It was also established that an instability line for plane-strain conditions can be defined in a way similar to that for axisymmetric conditions. Using the state parameter, a unified relationship between the normalised slope of instability line and the state parameters can be established for both axisymmetric and plane-strain conditions. Using this relationship, the instability conditions established under axisymmetric conditions can be used for plane-strain conditions.     

基于应变软化指标的岩石非线性蠕变模型

研究并提出一种以应变软化指标为基础的岩石非线性蠕变模型。首先,基于岩石三轴压缩试验结果,分析提出一种峰后应变软化指标R1来描述岩石峰后力学参数软化与塑性变形间的关系。接着,通过分析加速蠕变阶段与峰后段的联系,建立加速蠕变应变软化指标R2,并以此为基础构建非线性黏塑性体。然后,将其与Hook体和Kelvin元件串联,组建非线性蠕变模型,并将模型嵌入ABAQUS有限元程序。最后,针对砂岩和泥岩三轴蠕变试验建立数值模型,模拟曲线与试验曲线吻合良好,说明所建模型适用于硬岩/软岩加速蠕变现象的模拟。其中,参数 对加速蠕变曲线形态起调节作用,显示模型模拟脆性?延性流变破坏的特点。另外,模型参数均可由常规压缩破坏试验以及蠕变试验确定,易于获取。     

Nonstationary flow surface theory for modeling the viscoplastic behaviors of soils

This paper presents a three-dimensional elastic viscoplastic model that can describe the time-dependent behaviors of soft clays. The constitutive model is formulated based on the nonstationary flow surface theory and incorporates new developments, including (i) an improved definition of the nonstationary flow surface that is capable of capturing the stress–strain behaviors under different loading paths, (ii) a unique stress–strain—viscoplastic-strain-rate equation that is able to explicitly describe the nonstationary flow surface, and (iii) a final stable state concept that identifies the final equilibrium state at the end of creep and stress relaxation, which is also used to simplify the loading criteria. The consistency condition is validated for the proposed model, and the viscoplastic multipliers are calculated by solving the consistency equations. The model performance is investigated and validated via simulation of both oedometer and triaxial tests. The numerical results demonstrate that the proposed model is able to reproduce the main viscoplastic behaviors of soils, including creep, undrained creep rupture, stress relaxation, rate effect and accumulated effect.     

Delayed plastic model for time‐dependent behaviour of materials

A delayed plastic model, based on the theory of plasticity, is proposed to represent the time‐dependent behaviour of materials. It is assumed in this model that the stress can lie outside the yield surface and the conjugate stress called static stress is defined on the yield surface. The stress–strain relation is calculated based on the plastic theory embedding the static stress. Thus, the stress–strain relation of the model practically corresponds to that of the inviscid elastoplastic model under fairly low rate deformation. The delayed plastic model is coupled with the Cam‐clay model for normally consolidated clays. The performance of the model is then examined by comparing the model predictions with reported time‐dependent behaviour of clays under undrained triaxial conditions. It is shown that the model is capable of predicting the effect of strain rate during undrained shear and the undrained creep behaviour including creep rupture. Copyright © 2001 John Wiley & Sons, Ltd.     

基于上下负荷面的弹黏塑性本构模型及应力积分算法实现

为了描述天然软土的时间相依以及结构性特征,提出了一种能考虑土体超固结和结构性的实用弹黏塑性本构模型。它以Asaoka和Hashiguchi的上下负荷面作为某一应变速率下的参考屈服面,按照相对过应力的基本思路,新引入了两个能通过不同应变速率三轴压缩试验测定的率敏性参数 和 ,建立了以当前应力、黏塑性应变以及黏塑性应变速率为状态变量的动屈服准则函数,并给出了基于Newton-Raphson迭代的应力积分算法,且成功地将其嵌入到大型有限元软件ABAQUS中。最终通过数值算例来验证模型的正确性以及应力积分算法的可靠性。结果表明：该模型能同时描述土体的率敏性、蠕变以及结构性特征,模型参数物理意义明确、易懂可测,预测结果与试验数据吻合良好,可用于复杂边值问题的有限元计算。     

Q_2黄土流变特性及其统计损伤流变模型

通过对西安地区不同湿度原状Q2黄土的单轴蠕变试验和三轴蠕变试验,得到了Q2黄土的单轴、三轴蠕变曲线和应力-应变等时曲线。试验曲线分析表明,Q2黄土具有明显的非线性流变特性。Q2黄土的粘弹性特性可采用五元件广义的Kelvin模型来描述,其线性粘塑性特性可用与一滑块并联的Maxwell模型来描述,而非线性粘塑性特性可通过损伤与塑性应变的耦合作用来反映。在考虑损伤门槛值的统计损伤及黄土粘弹塑性组合模型的基础上,应用应变等效原理,建立了考虑瞬时损伤的统计损伤流变模型,并通过拟牛顿法确定了相应的模型参数。该模型能够很好地模拟黄土流变非线性特征,该模型不仅参数较少,而且适用性较广。     

西藏贡觉粉砂质泥岩工程地质特性与蠕变强度研究

川藏铁路在穿越西藏贡觉地区时遇到三叠系粉砂质泥岩,在高地应力条件下容易发生大变形等危害。文章开展了不同围压下的岩石三轴压缩和和三轴蠕变试验,结合PFC数值模拟,研究了粉砂质泥岩在不同围压下的蠕变特性和长期强度研究,结果表明:贡觉粉砂质泥岩流变具有西原蠕变模型特征,蠕变与常规三轴试验条件下,随着围压不断增大,粉砂质泥岩试样均由拉-剪破坏向单剪破坏过渡,剪切破裂面与水平线的夹角逐渐减小,微裂纹数量减少;蠕变试验相较于常规三轴试验,由拉应力引起的压碎带影响范围更广;在高围压条件下,粉砂质泥岩更容易发生流变,随着围压的增大,轴向应变、侧向应变和体积应变均增大,微裂纹数量呈下降趋势;瞬时弹性模量及黏弹性系数与围压呈线性递增关系,黏弹性模量与围压呈对数型增长关系,黏塑性系数与围压呈指数型增长关系。在荷载长期作用下,岩石长期强度低于瞬时强度。