Size and shape evolution of pores in a viscoplastic matrix under compression

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.     

绿片岩三轴流变力学特性的研究(II):模型分析

首先,基于在岩石全自动流变伺服仪上得到的绿片岩三轴流变试验曲线,采用五元件线性粘弹性模型对表现为粘弹性流变特性的曲线进行了辨识,获得了绿片岩的粘弹性流变参数;然后,提出了一个新的非线性粘性元件,并将其与塑性体并联起来,得到一个新的非线性粘塑性体(NVPB),该体能充分反映岩石的加速流变特性:同时,将NVPB模型与五元件粘弹性模型串联起来,建立了一个新的岩石七元件非线性粘弹塑性流变模型。采用绿片岩加速流变全过程曲线,对提出的岩石七元件非线性粘弹塑性流变模型进行了辨识,得到了岩石七元件非线性粘弹塑性流变模型的材料参数。流变模型与试验结果的比较,显示了所建模型的正确性与合理性。     

Constitutive approach for rate‐sensitive anisotropic structured clays

This paper describes a constitutive approach to model the behavior of rate‐dependent anisotropic structured clay. Rate‐sensitivity is modeled using overstress viscoplasticity. Clay structure is treated as a viscous phenomenon whereby the viscosity of the undisturbed structured clay is initially very high and the viscosity degrades or decreases with plastic straining until the intrinsic or residual viscosity is reached. A microstructure tensor approach is used to make the structured viscosity anisotropic; whereas, the intrinsic viscosity is assumed to be isotropic. The behavior of the constitutive model is compared with the measured response of two clays (Gloucester and St. Vallier clay) from Eastern Canada during triaxial compression tests on specimens trimmed at different orientations to the vertical. The comparisons show that the constitutive framework is able to describe the anisotropic and rate‐sensitive response of both clays. The response of the model is also examined for the more general case of anisotropic consolidated triaxial compression and extension. Copyright © 2010 John Wiley & Sons, Ltd.     

A visco‐plastic constitutive model for granular soils modified according to non‐local and gradient approaches

An already available non‐associated elastic–viscoplastic constitutive model with anisotropic strain hardening is modified in order to describe both the constitutive parameter dependency on relative density and the spatio‐temporal evolution of strain localization. To achieve this latter goal, two distinct but similar approaches are introduced: one inspired by the gradient theory and one by the non‐local theory. A one‐dimensional case concerning a simple shear test for a non‐homogeneous infinitely long dense sand specimen is numerically discussed and a finite difference scheme is employed for this purpose. The results obtained by following the two different approaches are critically analysed and compared. Copyright © 2001 John Wiley & Sons, Ltd.     

无限长土坡滑移变形模拟

从建立土体应力平衡微分方程入手,在平面应变的假定下求解均质无限长土坡的应力平衡方程。通过拟定坡土的滑移速度和粘塑性应变率的关系式,以剑桥模型为依托,得出滑移变形速率的显式表达式,由此推出土坡滑动位移的计算式。因为模型计算值与实测值无论变形趋势还是大小都有不少误差,所以经修正粘滞性参数来改进模型的模拟精度。最后通过一个算例,将现场实测滑移数据和模型计算值进行比较,分析表明二者较为一致。     

常规次固结沉降计算方法的改进研究

软土具有高含水率、大孔隙比等特点,在荷载作用下超静孔隙水压力消散后会产生次固结变形。实际工程中很多工程问题是由土体的次固结引起,因此对其沉降进行计算和预测很有必要。针对常规次固结沉降计算方法的不足,进行了改进,使随时间的发展,次固结沉降趋于一个定值而不是无限增长。对Yin & Graham一维弹黏塑性模型进行了简要介绍,并进行了改进,结合改进模型提出一种次固结沉降起讫时间的确定方法。与常规计算公式相比,改进公式可行,且更加合理,可用于实际工程地基次固结沉降的计算。     

基于一致性理论的Drucker-Prager材料弹黏塑本构模型

系统地论述了一致性理论的主要内容,推导了相关的黏塑性切线模量,并在前人工作的基础上研究了类Hoffman材料的推广条件及其一般性的推广方法;将线性Drucker-Prager模型(简称D-P)推广到一致性的黏塑性模型;给出了推广后模型的向后欧拉算法的有限元离散列式;最后将数值算例和试验成果进行对比,验证了一致性的线性D-P模型的有效性.     

Mantle convection with a brittle lithosphere: thoughts on the global tectonic styles of the Earth and Venus

Plates are an integral part of the convection system in the fluid mantle, but plate boundaries are the product of brittle faulting and plate motions are strongly influenced by the existence of such faults. The conditions for plate tectonics are studied by considering brittle behaviour, using Byerlee's law to limit the maximum stress in the lithosphere, in a mantle convection model with temperature-dependent viscosity.
When the yield stress is high, convection is confined below a thick, stagnant lithosphere. At low yield stress, brittle deformation mobilizes the lithosphere which becomes a part of the overall circulation; surface deformation occurs in localized regions close to upwellings and downwellings in the system. At intermediate levels of the yield stress, there is a cycling between these two states: thick lithosphere episodically mobilizes and collapses into the interior before reforming.
The mobile-lid regime resembles convection of a fluid with temperature-dependent viscosity and the boundary-layer scalings are found to be analogous. This regime has a well defined Nusselt number–Rayleigh number relationship which is in good agreement with scaling theory. The surface velocity is nearly independent of the yield stress, indicating that the 'plate' motion is resisted by viscous stresses in the mantle.
Analysis suggests that mobilization of the Earth's lithosphere can occur if the friction coefficient in the lithosphere is less than 0.03–0.13—lower than laboratory values but consistent with seismic field studies. On Venus, the friction coefficient may be high as a result of the dry conditions, and brittle mobilization of the lithosphere would then be episodic and catastrophic.     

Damage‐viscoplastic consistency model for rock fracture in heterogeneous rocks under dynamic loading

This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd.