泥质粉砂岩统一黏塑性模型的研究

统一黏塑性理论假设材料的应变仅由弹性应变和黏塑性应变组成。包含应力的黏塑性应变率表达方式,可以清晰地体现岩石材料应力增长率随应变加载速率相关的特性、在持续的偏压下的蠕变特性以及在恒定的应变下的松弛特性,但微分方程表示的本构模型在实际运用中存在计算上的困难。将无损伤Lemaitre统一黏塑性模型转换成关于应力的积分形式,并运用一致渐进展开算法得到其近似解,最后建立了基于近似解的牛顿迭代格式。开展了不同围压条件下泥质粉砂岩的三轴松弛试验,并将试验数据用于黏塑性模型的参数反分析,结果表明模型能够较好地反映泥质粉砂岩在长期作用下的力学行为特征。     

Rock dilatancy in uniaxial tests

SummaryRock Dilatancy in Uniaxial Tests A procedure to find an elastic/viscoplastic constitutive equation able to describe rock dilatancy in uniaxial unconfined creep tests is given starting from experimental data. The influence of the loading rate on the rock dilatancy can also be described with this model. It is shown that a creep of the volume takes place during tests and that this volume creep has the following property: the volume strain is oscillating before ariving at its final stable value. This property is described by the proposed mathematical model using a variable viscosity coefficient which may change sign. According to the model the volume creep takes place during loading only.With 5 Figures     

A Generalized Plasticity-Based Model for Sandstone Considering Time-Dependent Behavior and Wetting Deterioration

Based on the concept of generalized plasticity, this study proposes a constitutive model to describe the time-dependent behavior and wetting deterioration of sandstone. The proposed model (1) exhibits nonlinear elasticity under hydrostatic and shear loading, (2) follows the associated flow rule for viscoplastic deformation, (3) adopts a creep modulus that varies with the stress ratio, (4) considers the primary and secondary creep behaviors of rock, and (5) considers the effect of wetting deterioration. This model requires 13 material parameters, comprising 3 for elasticity, 7 for plasticity, and 3 for creep. All parameters can be determined easily by following the suggested procedures. The proposed model is first validated by comparison with triaxial tests of sandstone under different hydrostatic stress and cyclic loading conditions. In addition, the model is versatile in simulating time-dependent behavior through a series of multistage creep tests. Finally, to consider the effects of wetting deterioration, triaxial and creep tests under dry and water-saturated conditions are simulated. Comparison of the simulated and experimental data shows that the proposed model can predict the behavior of sandstone in dry and saturated conditions.     

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.     

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

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

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

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

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.     

一种基于修正西原模型的冻土蠕变本构关系

冻土在高应力水平下通常会表现出加速蠕变特性。经典西原（Nishihara）模型可以反映衰减蠕变和稳态蠕变,但不适用于描述加速蠕变。通过对Nishihara模型中黏性元件进行修正,将黏弹性元件中的定常的黏滞系数修正为时间的函数,并在黏塑性元件中采用非定常、非牛顿黏性元件,得到的蠕变柔量可以考虑应力水平和时间的影响,使其能够反映其加速蠕变的过程。用提出的模型分别对一维和三维应力、不同应力水平下的冻土蠕变试验数据进行了拟合,拟合结果能够反映冻土从衰减蠕变和稳定蠕变状态进入加速蠕变的过程,说明了该模型对于描述冻土在不同应力下不同蠕变状态的适用性。     

Evaluation of an anisotropic elastoplastic–viscoplastic bounding surface model for clays

An anisotropic time-dependent bounding surface model for clays is developed by generalizing a previous time-independent model that adopts a flexible bounding surface. It is based on the framework for coupled elastoplasticity–viscoplasticity for clays and Perzyna’s overstress theory. Three viscoplastic parameters were introduced and explained in detail. The model was validated against undrained creep tests for both isotropically and anisotropically consolidated clays, undrained and drained stress relaxation tests on some undisturbed clays, and undrained triaxial tests with varying strain rates on natural Hong Kong marine deposit clay. The general agreement between the model simulations and test results was satisfactory. The varying effects of lower-level parameters were discussed on the undrained multistage stress relaxation response for normally consolidated soils which had been ignored in literature. The flexibility of the model in capturing the shear strengths, which is the unique feature of the current model, was shown in the simulations of time-dependent triaxial tests on Taipei silty clay. All the simulations show that the proposed model is a relatively practical model considering both anisotropy and time dependency of clays.     

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.     

New Viscoplastic Model for Design Analysis of Tunnels in Squeezing Conditions

Summary  This paper is intended to describe the SHELVIP (Stress Hardening ELastic VIscous Plastic) model, a new viscoplastic constitutive law which has been developed to incorporate the most important features of behaviour observed in tunnels excavated in severe to very severe squeezing conditions. This model couples the elastoplastic and time-dependent behaviour by using a plastic yield surface, as frequently adopted in tunnel design analysis, and the definition of a state of overstress referred to a viscoplastic yield surface. The model is formulated in all its detailed aspects. The related analytical closed-form solution for representing triaxial creep deformations is developed. Also developed is an incremental numerical solution for describing the triaxial stress–strain behaviour under constant strain rate conditions. The model is shown to fit very satisfactorily the results of creep tests on clay shales and relaxation tests on coal specimens, as recently performed for design analysis of tunnels in squeezing conditions. Correspondence: D. Debernardi, Department of Structural and Geotechnical Engineering, Politecnico di Torino, Torino, Italy     

Constitutive modelling of granular materials using a contact normal-based fabric tensor

This paper presents a fabric tensor-based bounding surface model accounting for anisotropic behaviour (e.g. the dependency of peak strength on loading direction and non-coaxial deformation) of granular materials. This model is developed based on a well-calibrated isotropic bounding surface model. The yield surface is modified by incorporating the back stress which is proportional to a contact normal-based fabric tensor for characterising fabric anisotropy. The evolution law of the fabric tensor, which is dependent on both rates of the stress ratio and the plastic strain, rules that the material fabric tends to align with the loading direction and evolves towards a unique critical state fabric tensor under monotonic shearing. The incorporation of the evolution law leads to a rotational hardening of the yield surface. The anisotropic critical state is assumed to be independent of the initial values of void ratio and fabric tensor. The critical state fabric tensor has the same intermediate stress ratio (i.e. b value) and principal directions as the critical state stress tensor. A non-associated flow rule in the deviatoric plane is adopted, which is able to predict the non-coaxial flow naturally. The stress–strain relation and fabric evolution of model predictions show a satisfactory agreement with DEM simulation results under monotonic shearing with different loading directions. The model is also validated by comparing with laboratory test results of Leighton Buzzard sand and Toyoura sand under various loading paths. The comparison results demonstrate encouraging applicability of the model for predicting the anisotropic behaviour of granular materials.

     

Damage–viscoplastic consistency model with a parabolic cap for rocks with brittle and ductile behavior under low‐velocity impact loading

This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd.     

Clay minerals in fluvial shaley sandstone of upper Triassic reservoir (TAGS)—Toual field SE Algeria: identification from wireline logs and core data

The creep property of rock under cyclic loading is very important in civil engineering. In order to establish a novel constitutive equation for rock under cyclic loading, a fractional-order viscoplastic body under cyclic loading was constructed based on fractional-order viscous element. A fractional-order visco-elastoplastic model (FVEPM) for rock was established by connecting constructed fractional-order viscoplastic body with Burgers model. The model was a Burgers model when the maximum value of cyclic loading was less than the critical strength of rock; otherwise, it was a FVEPM which can be used to reflect the transient, steady-state, and tertiary creep phases of rock. The cyclic loading was decomposed into a static load and a cyclic loading with a zero average stress. According to rheological mechanics theory, the rheology constitutive equation of rock under the static load can be derived. According to viscoelastic mechanics theory, the constitutive equation under cyclic loading with a zero average stress was established by introducing the variation parameters of energy storage and energy dissipation compliance caused by rock damage and fracture. Finally, a new dynamic constitutive equation of rock cyclic loading can be obtained by superimposing the constitutive equation under static load and cyclic loading with a zero average stress. Compared with existing test results of rock under cyclic loading, the proposed constitutive model can be used to describe the creep characteristics of rock under cyclic loading and reflect the presented fluctuation of strain curve of rock under cyclic loading.     

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

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

三轴压缩下粉砂质泥岩蠕变本构模型研究

依据粉砂质泥岩三轴压缩蠕变试验结果,分析了岩石的蠕变力学特性。选取线性黏弹性Burgers模型来描述岩石的衰减蠕变和稳定蠕变特性。引入非线性黏塑性元件,将其与Burgers模型串联起来,建立一个新的六元件非线性黏弹塑性Burgers蠕变本构模型。采用非线性Burgers模型来描述岩石的加速蠕变特性。应用Levenberg-Marquardt算法对蠕变试验曲线进行辨识,得出各级应力水平下岩石轴向与径向的蠕变模型参数。对比试验曲线和模型拟合曲线,二者吻合较好,表明线性Burgers模型与非线性Burgers模型可以准确描述岩石3个蠕变阶段的力学特性。模型参数辨识结果表明,在同一级应力水平下岩石轴向与径向的蠕变模型参数各自独立并不统一,模型参数取值时应考虑岩石蠕变的各向异性以及轴向与径向蠕变量、蠕变速率大小的差异。在前7级应力水平下,模型参数应取以轴向应变辨识得到的Burgers模型参数;在第8级应力水平下,模型参数应取以径向应变辨识得到的Burgers模型参数;在破裂应力水平下,模型参数应取以径向应变辨识得到的非线性Burgers模型参数。     

Anisotropic viscoplastic modeling of rate‐dependent behavior of clay

The objective of this study is to derive an effective stress‐based constitutive law capable of predicting rate‐dependent stress–strain, stress path and undrained shear strength and creep behavior. The flow rule used in the MIT‐E3 model and viscoplasticity theory is employed in the derivation. The model adopts the yield surface capable of representing the yield behavior of the Taipei silty clay and assumes that it is initially symmetric about the K₀‐line. A method is then developed to compute the gyration and expansion of the loading surface to simulate the anisotropic behavior due to the principal stress rotation after shear. There are 11 parameters required for the model to describe the soil behavior and six of them are exactly the same as those used in the Modified Cam‐clay model. The five additional parameters can be obtained by parametric studies or conventional soil tests, such as consolidation tests, triaxial compression and extension tests. Finally, verification of the model for the anisotropic behavior, creep behavior and the rate‐dependent undrained stress–strain and shear strength of the Taipei silty clay is conducted. Copyright © 2010 John Wiley & Sons, Ltd.