A hierarchical approach for constitutive modelling of geologic materials

A hierarchical concept is proposed for the development of constitutive models to account for various factors that influence behaviour of (geologic) materials. It permits evolution of models of progressively higher grades from the basic model representing isotropic hardening with associative behaviour. Factors such as non-associativeness and induced anisotropy due to friction and cyclic loading, and softening are introduced as corrections or perturbations to the basic model. The influence of these factors is captured through non-associativeness manifested by deviation from normality of the plastic strain increments to the yield surface, F. Details of four models: isotropic hardening with associative behaviour, isotropic hardening with non-associative behavioural anisotropic hardening and strain-softening with a damage variable are presented. They are verified with respect to laboratory multiaxial test data under various paths of loading, unloading and reloading for typical soils, rock and concrete. The proposed concept is general, yet sufficiently simplified in terms of physical understanding, number of constants and their physical meanings, determination of the constants and implementation.     

An Anisotropic Model for Granular Material Based on Experiments

Soil is a heterogeneous material and most natural soil deposits show a definite stratification. The mechanical behaviour of such material is generally different in different directions, especially in the direction parallel and perpendicular to the stratification. A series of isotropic compression tests were carried out to study the behavior of granular material produced under controlled stratification in the laboratory. These tests were conducted both on cylindrical and square prismatic tri-axial specimens. It was observed that for hydrostatic loading, the strain response was different in different directions, especially in directions parallel and perpendicular to the direction of soil deposition. A definite trend of anisotropy was observed in the deformation pattern. The observed anisotropy is modeled in this paper by treating soil-dilatancy as a variable quantity. The equation of the plastic potential surface of the model which obeys a non-associated flow rule, is assumed to be dependent on three main variables confining pressure (\(\sigma_{3}\)), void ratio (e) and the angle of bedding plane orientation (δ) during deposition. The angle of bedding plane orientation (δ) was measured with respect to the direction of the major principal stress. The model has a cap yield surface in the isotropic stress direction, which is supplemented by a shear hardening Mohr–Coulomb surface in the deviator direction. This paper focuses on predicting the anisotropic strain response of stratified soil deposits subjected to isotropic compression. The proposed anisotropic model incorporates within an existing strain-hardening sand model, a modified cap yield surface and a modified plastic potential function related to the cap surface, to account for the anistropic response observed in isotropic compression tests. The two dimensional stress–strain model was extended to three dimensional Cartesian space. The strain anisotropy observed in the isotropic compression tests was predicted by the three dimensional anisotropic model proposed for granular materials.     

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

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

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

A destructuration theory and its application to SANICLAY model

Many natural clays have an undisturbed shear strength in excess of the remoulded strength. Destructuration modeling provides a means to account for such sensitivity in a constitutive model. This paper extends the SANICLAY model to include destructuration. Two distinct types of destructuration are considered: isotropic and frictional. The former is a concept already presented in relation to other models and in essence constitutes a mechanism of isotropic softening of the yield surface with destructuration. The latter refers to the reduction of the critical stress ratio reflecting the effect of destructuration on the friction angle, and is believed to be a novel proposition. Both the types depend on a measure of destructuration rate expressed in terms of combined plastic volumetric and deviatoric strain rates. The SANICLAY model itself is generalized from its previous form by additional dependence of the yield surface on the third isotropic stress invariant. Such a generalization allows to obtain as particular cases simplified model versions of lower complexity including one with a single surface and associative flow rule, by simply setting accordingly parameters of the generalized version. A detailed calibration procedure of the relatively few model constants is presented, and the performance of three versions of the model, in descending order of complexity, is validated by comparison of simulations to various data for oedometric consolidation followed by triaxial undrained compression and extension tests on two structured clays. Copyright © 2009 John Wiley & Sons, Ltd.     

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.

     

Simulation of yielding and stress–stain behavior of shanghai soft clay

In this paper, a simple bounding surface plasticity model is used to reproduce the yielding and stress–strain behavior of the structured soft clay found at Shanghai of China. A series of undrained triaxial tests and drained stress probe tests under isotropic and anisotropic consolidation modes were performed on undisturbed samples of Shanghai soft clay to study the yielding characteristics. The degradation of the clay structure is modeled with an internal variable that allows the size of the bounding surface to decay with accumulated plastic strain. An anisotropic tensor and rotational hardening law are introduced to reflect the initial anisotropy and the evolution of anisotropy. Combined with the isotropic hardening rule, the rotational hardening rule and the degradation law are incorporated into the bounding surface formulation with an associated flow rule. Validity of the model is verified by the undrained isotropic and anisotropic triaxial test and drained stress probe test results for Shanghai soft clay. The effects of stress anisotropy and loss of structure are well captured by the model.     

Mathematical modelling of monotonic and cyclic undrained clay behaviour

The proposed general analytical model describes the anisotropic, elasto-plastic, path-dependent, stress-strain-strength properties of inviscid saturated clays under undrained loading conditions. The model combines properties of isotropic and kinematic plasticity by introducing the concept of a field of plastic moduli which is defined in stress space by the relative configuration of yield surfaces. For any loading (or unloading) history, the instantaneous configuration is determined by calculating the translation and contraction (or expansion) of each yield surface. The stress-strain behaviour of clays can thus be determined for complex loading paths and in particular for cyclic loadings. The stress-strain relationships are provided for use in finite element analyses. The model parameters required to characterize the behaviour of any given clay can be derived entirely from conventional triaxial or simple shear soil test results. The model's extreme versatility is demonstrated by using it to formulate the behaviour of the Drammen clay under both monotonic and cyclic loading conditions. The parameters are determined by using solely the results from monotonic and cyclic strain-controlled simple shear experimental tests, and the model's accuracy is evaluated by applying it to predict the results of other tests such as (1) cyclic stress-controlled simple shear tests, (2) monotonic triaxial loading compression and unloading extension tests, and (3) cyclic stress- and strain-controlled triaxial tests on, this same clay. The theoretical predictions are found to agree extremely well with the experimental test results.     

A unified mechanical and retention model for saturated and unsaturated soil behaviour

The coupled mechanical and water retention elasto-plastic constitutive model of Wheeler, Sharma and Buisson (the Glasgow coupled model, GCM) predicts unique unsaturated isotropic normal compression and unsaturated critical state planar surfaces for specific volume and degree of saturation when soil states are at the intersection of mechanical (M) and wetting retention (WR) yield surfaces. Experimental results from tests performed by Sivakumar on unsaturated samples of compacted speswhite kaolin confirm the existence and form of these unique surfaces. The GCM provides consistent representation of transitions between saturated and unsaturated conditions, including the influence of retention hysteresis and the effect of plastic volumetric strains on retention behaviour, and it gives unique expressions to predict saturation and de-saturation conditions (air-exclusion and air-entry points, respectively). Mechanical behaviour is modelled consistently across these transitions, including appropriate variation of mechanical yield stress under both saturated and unsaturated conditions. The expressions defining the unsaturated isotropic normal compression planar surfaces for specific volume and degree of saturation are central to the development of a relatively straightforward methodology for determining values of all GCM parameters (soil constants and initial state) from a limited number of laboratory tests. This methodology is demonstrated by application to the experimental data of Sivakumar. Comparison of model simulations with experimental results for the full set of Sivakumar’s isotropic loading stages demonstrates that the model is able to predict accurately the variation of both specific volume and degree of saturation during isotropic stress paths under saturated and unsaturated conditions.     

A macroscopic damage‐plastic constitutive law for modeling quasi‐brittle fracture and ductile behavior of concrete

A new phenomenological macroscopic constitutive model for the numerical simulation of quasi‐brittle fracture and ductile concrete behavior, under general triaxial stress conditions, is presented. The model is particularly addressed to simulate a wide range of confinement stress states, as also, to capture the strong influence of the mean stress value in the concrete failure mechanisms. The model is based on a two‐surface damage‐plastic formulation. The mechanical behavior in different domains of the stress space is separately described by means of a quasi‐brittle or ductile material response:

相似文献   

The Shear Strength Behavior of Two Peaty Soils

This paper presents experimental research concerning the shear behavior of two types of Italian peat, one normally, the other over consolidated. Organic soils are characterized by very high compressibility and high fiber content; two features that give rise to several problems during laboratory tests. Under consideration here are the effects of fibers and over consolidation on friction angle and stress–strain behavior. These are evaluated by means of undrained compression triaxial tests, with isotropic and anisotropic (K_o) consolidation being performed on natural and remolded samples. The experimental results are also analyzed by means of a bilinear failure criterion based on soil-reinforcement interaction mechanisms.     

Three‐dimensional elasto‐plastic model for unsaturated compacted soils with different initial densities

This paper presents an elasto‐plastic model for unsaturated compacted soils and experimental results obtained from a series of suction‐controlled triaxial tests on unsaturated compacted clay with different initial densities. The initial density dependency of the compacted soil behaviour is modelled by establishing experimental relationships between the initial density and the corresponding yield stress and thereby between the initial density and the location and slope of normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure surface and the yield surface in the deviatoric plane are given by the extended SMP criterion. A considerable number of the isotropic compression, triaxial compression and extension tests on unsaturated compacted clay with different initial densities were performed using a suction‐controllable triaxial apparatus, to measure the stress–strain–volume change in different stress paths and wetting paths. The model has well‐predicting capabilities to reproduce the mechanical behaviour of specimens compacted under different conditions not only in isotropic compression but also in triaxial compression and triaxial extension. Copyright © 2003 John Wiley & Sons, Ltd.     

考虑饱和黏土初始各向异性的循环边界面塑性模型

天然土体的初始各向异性通常可对其后继循环特性产生显著影响。现有考虑循环载荷作用的土体弹塑性模型,往往采用类似修正剑桥模型的椭圆形屈服面,已有研究表明,该椭圆形屈服面因其拉伸弹性区域偏大,针对天然K0固结状态的土体,其计算精度较差。基于新近提出的广义各向同性硬化准则,在边界面方程中引入初始各向异性张量,并采用空间滑动面破坏准则（SMP）的变换应力法,建立了能考虑饱和黏土初始各向异性的循环边界面塑性模型。分别针对等压和偏压固结的饱和黏土静、动三轴试验进行模拟,结果表明,该模型能合理反映土体的初始各向异性及其后继循环动力特性。     

任意各向异性介质三维有限元航空电磁响应模拟

在航空电磁探测方法应用中,地下探测环境复杂,地下介质物性参数存在各向异性特征。如果采用常规各向同性介质模型,在资料解释过程中会产生严重的偏差。本文基于矢量有限元法开展三维任意各向异性频率域航空电磁响应模拟计算研究。首先,将总场分解成一次场和二次场,对空气介质的均匀全空间进行一次场解析计算,同时利用矢量有限元法对二次电场的双旋度方程进行求解。为提高求解速度,采用共享内存直接求解器PARDISO对大规模稀疏矩阵并行计算,大大提高了三维模型的计算速度。之后开展了4种典型目标模型的三维各向异性介质模拟：围岩各向同性-目标体各向异性（绕z轴旋转）模型;围岩各向同性-目标体各向异性（绕x轴旋转）模型;围岩各向异性-目标体各向异性（绕z轴旋转）模型;围岩各向异性-目标体各向异性（绕x轴旋转）模型。对比分析了不同模型在不同旋转角度情况下,磁场实虚分量的变化特征,进而总结了各向异性参数对航空电磁响应的影响规律和识别方法。     

A yield criterion for isotropic and cross‐anisotropic cohesive‐frictional materials

This paper proposes a yield and failure criterion for cohesive and frictional materials. The function is given by the combination of a Lode dependence for the behaviour in the deviatoric plane and a meridian function for the pressure‐dependent behaviour. A variety of shapes can be achieved with the proposed criterion including Lode dependences which are able to reproduce the behaviour of isotropic and cross‐anisotropic materials in the deviatoric plane. The criterion is validated through the comparison with experimental data based on multiaxial experimental tests on clays, sands, rocks and concrete. Finally, the convexity of the criterion is analysed and discussed. Copyright © 2009 John Wiley & Sons, Ltd.     

Implementation of DSC model and application for analysis of field pile tests under cyclic loading

The disturbed state concept (DSC) model, and a new and simplified procedure for unloading and reloading behavior are implemented in a nonlinear finite element procedure for dynamic analysis for coupled response of saturated porous materials. The DSC model is used to characterize the cyclic behavior of saturated clays and clay–steel interfaces. In the DSC, the relative intact (RI) behavior is characterized by using the hierarchical single surface (HISS) plasticity model; and the fully adjusted (FA) behavior is modeled by using the critical state concept. The DSC model is validated with respect to laboratory triaxial tests for clay and shear tests for clay‐steel interfaces. The computer procedure is used to predict field behavior of an instrumented pile subjected to cyclic loading. The predictions provide very good correlation with the field data. They also yield improved results compared to those from a HISS model with anisotropic hardening, partly because the DSC model allows for degradation or softening and interface response. Copyright © 2000 John Wiley & Sons, Ltd.     

Anisotropic UH model for soils based on a simple transformed stress method

A simple method called anisotropic transformed stress (ATS) method is proposed to develop failure criteria and constitutive models for anisotropic soils. In this method, stress components in different directions are modified differently in order to reflect the effect of anisotropy. It includes two steps of mapping of stress. First, a modified stress tensor is introduced, which is a symmetric multiplication of stress tensor and fabric tensor. In the modified stress space, anisotropic soils can be treated to be isotropic. Second, a TS tensor is derived from the modified stress tensor for the convenience of developing anisotropic constitutive models to account for the effect of intermediate principal stress. By replacing the ordinary stress tensor with the TS tensor directly, the unified hardening model is extended to model the anisotropic deformation of soils. Anisotropic Lade's criterion is adopted for shear yield and shear failure in the model. The form of the original model formulations remains unchanged, and the model parameters are independent of the loading direction. Good agreement between the experimental results and predictions of the anisotropic unified hardening model is observed. Copyright © 2016 John Wiley & Sons, Ltd.     

Intelligent model selection with updating parameters during staged excavation using optimization method

Various constitutive models have been proposed, and previous studies focused on identifying parameters of specified models. To develop the smart construction, this paper proposes a novel optimization-based intelligent model selection procedure in which parameter identification is also performed during staged excavation. To conduct the model selection, a database of seven constitutive models accounting for isotropic or anisotropic yield surface, isotropic or anisotropic elasticity, or small strain stiffness for clayey soils is established, with each model numbered and deemed as one additional parameter for optimization. A newly developed real-coded genetic algorithm is adopted to evaluate the performance of simulation against field measurement. As the process of optimization goes on, the soil model exhibiting good performance during simulation survives from the database and model parameters are also optimized. For each excavation stage, with the selected model and optimized parameters, wall deflection and ground surface settlement of the subsequent unexcavated stage are predicted. The proposed procedure is repeated until the entire excavation is finished. This proposed procedure is applied to a real staged excavation with field data, which demonstrates its effectiveness and efficiency in engineering practice with highlighting the importance of anisotropic elasticity and small strain stiffness in simulating excavation. All results demonstrate that the current study has both academic significance and practical significance in providing an efficient and effective approach of adaptive optimization-based model selection with parameters updating in engineering applications.