Mixed time integration for the transient analysis of jointed media

A constitutive operator split method with implicit-explicit time integration is presented for the transient analysis of rigid block models of jointed media. The linear portion of the joint constitutive law is integrated by an implicit method and the non-linear, time dependent portion is integrated by an explicit method. The method features the stability of implicit procedures as well as the flexibility of explicit procedures for non-linear problems. The solutions obtained with this method are compared to the solutions obtained by the explicit central difference method; in all cases there is good to excellent agreement. For some problems, particularly for those with low frequency excitation, it is shown that the implicit-explicit method can result in a substantial savings over more conventional explicit methods.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

Identification of constitutive parameters of soil using an optimization strategy and statistical analysis

The paper at hand investigates a strategy to calibrate different constitutive models for soils via back analysis. The efficiency and reliability of the parameter identification for soil models is worked out. In order to demonstrate of how to utilise identification procedures, results from optimisation against conventional oedometer and drained triaxial compression tests on natural Pappadai clay are presented and discussed. The aim of geotechnical optimisation problems is to obtain a set of model parameter values that provide the best match between soil model simulations and appropriate measurements. For the parameter identification process, a constrained population-based algorithm is chosen, namely the Particle Swarm Optimiser. The identification is carried out in an initial step separately on each test and then simultaneously on oedometer and triaxial tests. The evaluation is performed employing three different constitutive models of varying complexity and number of constitutive parameters.A subsequent residual analysis and the computation of confidence intervals for the parameters provide valuable results to assess the quality of the identified parameters in correlation with the evaluated data. Therefore, criteria of the utility and reliability of the mathematical models for further prognosis computations can be estimated.     

敏感环境下基坑数值分析中土体本构模型的选择

数值分析已成为敏感环境下基坑工程分析的最重要手段,其关键是选择合适的土体本构模型和计算参数。在分析了岩土数值分析中常用土体本构模型特点的基础上,通过算例较系统地对比了各类模型在基坑开挖数值分析中的适用性。敏感环境下的基坑工程需重点关注墙后土体的变形,从满足工程需要和方便实用的角度出发,建议采用能考虑黏土的塑性和应变硬化特征、能区分加荷和卸荷且刚度依赖于应力水平的硬化类弹塑性模型,如MCC模型和HS模型进行分析。具体工程实例的分析,表明了硬化类弹塑性模型在敏感环境下基坑开挖数值分析中的适用性。     

Numerical study on finite element implementation of hypoplastic models

A comprehensive numerical study on finite element implementation of hypoplastic models is presented. Two crucial aspects, local integration of the constitutive equations (the local problem) and forming tangent operators for Newton–Raphson iteration (the global problem), are investigated. For solving the local problem, different integration algorithms, including explicit and implicit methods, are examined using tri-axial compression tests and incremental stress response envelopes, as well as typical boundary value problems. For solving global problems, three different ways of generating the tangent operator are compared. The numerical evidences indicate that, in terms of accuracy, efficiency and robustness, explicit methods with substepping and error control are the best choices for constitutive integration of hypoplastic models while the so-called continuum tangent operators have certain advantages over two other types of numerically-generated consistent tangent operators.     

A comparative study of different model families for the constitutive simulation of viscous clays

The simulation of the viscous behavior of some clays is of high importance in many geotechnical problems. The literature offers a vast amount of constitutive models able to simulate the rate dependence observed on these materials. Although most of these models are calibrated to very similar experimental observations and share similar definitions of material parameters, some discrepancies of their response have been detected, which are related to their mathematical formulations. In this work, the causes of these discrepancies are carefully studied. To that end, four different model families are analyzed, namely, nonstationary flow surface (NSFS) models, viscoplasticity with overstress function (OVP), viscoplasticity with Norton's power law (NVP), and visco-hypoplasticity (VHP). For the sake of a fair comparison, single constitutive models using the same set of material parameters, and following other requirements, are developed for each model family. Numerical implementations of the four resulting models are performed. Their response at different tests are carefully analyzed through simulation examples and direct examination of their constitutive equations. The set includes some basic tests at isotropic stress states and others as responses envelopes, undrained creep rupture, and an oedometer test with loading, unloading-reloading, creep, and relaxation. The article is concluded with some remarks about the observed discrepancies of these model families.     

城市固体废弃物本构模型的研究现状及展望

城市固体废弃物(MSW)是由多种无机和有机成分组成的复合材料,具有高压缩性和可降解性。随着时间的变化,基本的物理力学特性发生变化,与传统土体相比,性质差异较大。城市固体废弃物的本构模型对于垃圾填埋场中土工问题的研究极其重要,考虑不同因素的影响,系统总结了国内外一些学者对MSW本构模型的研究现状,指出已有的研究成果主要考虑了荷载引起的瞬时变形、机械蠕变、生物降解、纤维物质的加筋作用等,但这些都是考虑单个因素或者将几个因素简单叠加的本构模型,并未考虑不同因素的耦合作用。结合土体本构模型的发展特点,指出了今后城市固体废弃物本构模型的研究方向和发展趋势：应当综合考虑各种因素对MSW应力-应变关系的影响,从工程实践出发,通过岩土工程、环境工程、化学工程和生物工程等学科的交叉寻求突破,形成考虑生化反应-骨架变形-水气运移-溶质迁移耦合作用的理论框架。     

对岩土工程数值分析的几点思考

首先,介绍了笔者对我国岩土工程数值分析现状的调查结果及分析,然后,分析了采用连续介质力学分析岩土工程问题的关键,并讨论分析了岩土本构理论发展现状,提出对岩土本构理论发展方向的思考,最后对数值分析在岩土工程分析中的地位作了分析。分析表明,岩土工程数值分析结果是岩土工程师在岩土工程分析过程中进行综合判断的重要依据之一;采用连续介质力学模型求解岩土工程问题的关键是如何建立岩土的工程实用本构方程;建立多个工程实用本构方程结合积累大量工程经验才能促使数值方法在岩土工程中由用于定性分析转变到定量分析。     

NGI‐ADP: Anisotropic shear strength model for clay

Many geotechnical problems involve undrained behavior of clay and the capacity in undrained loading. Most constitutive models used today are effective stress based and only indirectly obtain values for the undrained shear strength. To match the design profiles of undrained shear strengths, in active (A), direct simple shear (D) and passive (P) modes of loading are complicated. This paper presents the elastoplastic constitutive model NGI‐ADP which is based on the undrained shear strength approach with direct input of shear strengths. Consequently, exact match with design undrained shear strengths profiles is obtained and the well‐known anisotropy of undrained shear strength and stiffness is accounted for in the constitutive model. A non‐linear stress path‐dependent hardening relationship is used, defined from direct input of failure strains in the three directions of shearing represented by triaxial compression, direct simple shear and triaxial extension. With its clear input parameters the model has significant advantages for design analysis of undrained problems. The constitutive model is implemented, into finite element codes, with an implicit integration scheme. Its performance is demonstrated by a finite element analysis of a bearing capacity problem. Copyright © 2011 John Wiley & Sons, Ltd.     

Procedure,selection and application of plasticity models for a soil

A five-step procedure involving mathematical formulation, identification and determination of parameters and verification is presented for development and selection of appropriate and reliable constitutive law(s) for geologic media. Comprehensive analyses are performed toward determination of an appropriate law for a (artificial) soil. The most suitable model is obtained by critical evaluation of four different plasticity models; here verification and comparisons of predictions with observations from laboratory tests, and with those from two boundary value problems are used as the basis of the selection. The model thus selected is found to be appropriate for applications to relevant practical problems.     

土结构性本构模型研究现状综述

土本构模型的建立是一个重要而又复杂的问题,到目前为止,国内外学者们已提出数以百计的土本构模型,诸多文献也对这些模型进行了评述和归纳。然而这些土本构模型多是在扰动土或砂土的基础上发展和建立起来的,它们难以描述由于土结构性引起的各种非线性行为,其计算结果与实际情况相差甚远。天然土体一般都具有一定的结构性,所以有必要建立考虑土结构性影响的土本构模型。针对这个现实,目前有些学者已基于各种理论和方法,提出了一些可以考虑土结构性影响的土本构模型,并得了较好的应用。但在目前的文献中还很少有对土的结构性本构模型研究进行归纳,基于此,本文简要介绍了一下目前土的结构性本构模型研究现状,并提出了这些本构模型在应用中所存在的问题。     

在ABAQUS中开发实现Duncan-Chang本构模型

美国HKS公司开发的软件产品ABAQUS是目前世界上最强大的非线性有限元分析工具之一。该软件包括了众多材料本构模型,但尚缺少国内外土工数值分析中广泛采用的Duncan-Chang本构模型。这一缺憾影响了ABAQUS软件在土工分析中的应用范围。本文旨在介绍ABAQUS中开发Duncan-Chang材料本构模型实现方法,给出了开发过程的概要,完成了两个典型常规三轴压缩模型问题数值测试。结果表明：在ABAQUS中增加Duncan-Chang材料本构模型后,不仅可以充分地利用该软件强大的非线性求解平台,而且还有可能完成复杂土工应力应变的有限元数值分析问题,同时,具有计算速度快、计算精度高和前后处理快捷方便的优点,极大地降低了土工分析程序开发的难度,并减少了维护工作量。     

ParCYCLIC: finite element modelling of earthquake liquefaction response on parallel computers

This paper presents the computational procedures and solution strategy employed in ParCYCLIC, a parallel non‐linear finite element program developed based on an existing serial code CYCLIC for the analysis of cyclic seismically‐induced liquefaction problems. In ParCYCLIC, finite elements are employed within an incremental plasticity, coupled solid–fluid formulation. A constitutive model developed for simulating liquefaction‐induced deformations is a main component of this analysis framework. The elements of the computational strategy, designed for distributed‐memory message‐passing parallel computer systems, include: (a) an automatic domain decomposer to partition the finite element mesh; (b) nodal ordering strategies to minimize storage space for the matrix coefficients; (c) an efficient scheme for the allocation of sparse matrix coefficients among the processors; and (d) a parallel sparse direct solver. Application of ParCYCLIC to simulate 3‐D geotechnical experimental models is demonstrated. The computational results show excellent parallel performance and scalability of ParCYCLIC on parallel computers with a large number of processors. Copyright © 2004 John Wiley & Sons, Ltd.     

岩石统计渗流模型和统计损伤本构模型研究

采用Eshelby等效夹杂方法建立考虑渗流的岩石损伤本构模型是一种有效方法,但相关文献目前还极少。利用细观力学的Eshelby等效夹杂方法, 探索了考虑渗流和损伤的Helmholtz自由比能函数的确定,用连续介质损伤力学方法建立了相应的考虑渗流的岩石损伤统计本构模型,提出了考虑渗流和损伤过程的岩石破坏准则,建议了考虑损伤和应变的岩石渗透系数演化方程。与试验结果比较表明,所建立模型是合理的。     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.     

Numerical methods for the analysis of tunnel performance in squeezing rocks

Summary Some comments are presented on the finite element analysis of the time dependent effects that develop when a tunnel is driven in a rock mass characterized by a viscous behaviour. First, the so called swelling and squeezing phenomena are described considering in particular rocks containing clay minerals. Subsequently the discussion is focused on the squeezing behaviour, i.e. on the time dependent increase of the shear deformation which develops with minor volume changes. Some linear and non-linear viscous constitutive laws are discussed, as well as their use in the solution of boundary value problems. Finally, the effects taking place around tunnels driven into squeezing rocks are illustrated through some numerical analyses based on these constitutive models.     

An evaluation on the degradation evolutions in three constitutive models for bonded geomaterials by DEM analyses

This paper presents a numerical evaluation on the degradation evolutions in three well-known constitutive models for bonded geomaterials using the Distinct Element Method (DEM). A series of isotropic, constant stress ratio and biaxial compression tests on the bonded geomaterials were carried out by a two-dimensional (2D) DEM code, NS2D, for this purpose. The constitutive models examined are the disturbed state constitutive model (DSC model), the superloading yield surface model and Nova’s model. First, a microscopic interpretation of the degradation variable, which is used to describe the degradation evolution in the models, was presented based on the micromechanics theory to compare the degradation evolutions used by the models with the DEM results. Then, the macromechanical responses of the numerical specimens in the DEM tests were examined in comparison with the experimental data on the artificially bonded sands. Finally, the degradation evolutions observed in the DEM tests were compared with the results predicted by the constitutive models. The study shows that the DEM results are able to capture the main features of the artificially bonded sands. The degradation evolutions used in the three constitutive models are all overally in agreement with the DEM results. In addition, two simple evolution equations for the degradation variable, which can well describe the degradation evolution observed in the DEM tests, were presented and discussed.     

节理卸荷岩体非线性本构关系再研究

对地质材料模型、石膏模型、砂浆模型、进行了模型试验,通过对模型试验得出的数据进行分析、拟合,得出初步的节理卸荷岩体的增量本构。针对其不足,从参数选择、物理意义等方面进一步进行修正得出一改进的增量本构关系,并在有限元分析程序ADINA中接入该本构关系,对清江隔河岩电站厂房高边坡三维模型进行有限元分析。     

Calibration of cyclic constitutive models by oscillating functions

Most structures are subjected to more cyclic loads during their life time than static loads. These cyclic action could be a result of either natural or man-made activities and may lead to soil failure. In order to understand the response of the foundation and its interaction with these complex cyclic loadings, various researchers have over the years developed different constitutive models. Although a lot of research is being carried out on these relatively new models, little or no details exist in literature about the model-based identification of the cyclic constitutive parameters which to a large extent govern the quality of the model output. This could be attributed to the difficulties and complexities of the inverse modeling of such complex phenomena. A variety of optimisation strategies are available for the solution of the sum of least-squares problems as usually done in the field of model calibration. However, for the back analysis (calibration) of the soil response to oscillatory load functions, this article gives insight into the model calibration challenges and also puts forward a method for the inverse modeling of cyclic loaded foundation response such that high-quality solutions are obtained with minimum computational effort.