冻土弹塑性本构模型研究现状

土体的本构模型是分析计算土工结构变形规律的关键。自剑桥模型提出以来,对于土体本构模型的研究,各国学者在不同环境条件下发展研究了适用于相应试验条件下的各类模型。随着寒区经济的发展,关于冻土力学性质的研究日益增多,许多学者借鉴常规融土的研究方法建立了不同条件下的冻土弹塑性本构模型。为了进一步理清冻土变形行为的特征,完善适用于冻土弹塑性行为的本构模拟理论和方法,作者总结和分析了各类冻土弹塑性本构模型的理论基础、建模方法、参数确定方法等内容,对进一步发展可准确描述冻土复杂力学行为的本构模型具有指导意义。     

A critical state two-surface plasticity model for gravelly soil-structure interfaces under monotonic and cyclic loading

An elasto-plastic model is proposed for modeling the constitutive behavior of the interface between gravelly soils and structural materials. This model is based on the two-surface plasticity formulation and it is compatible with the concept of critical state soil mechanics. The model requires the same set of eight calibration parameters for predicting the monotonic and cyclic responses of both loose and dense interfaces. The model simulates cyclic densification, shear degradation and the effects of normal pressure, soil density, and stress path. The performance of the proposed constitutive model is validated by tests data under different normal stresses and boundary conditions.     

土的量化记忆模型及其参数确定

土的真非线性分析要求追踪应力-应变历史,寻求滞回圈上应力-应变的瞬时切线模量。由于地震荷载的随机性,实际土体的应力-应变历史轨迹非常复杂。基于传统模型的骨架曲线和滞回圈表达式进行这类分析的难度很大。作者介绍的量化记忆模型是为土的动力真非线性分析提出的,它是将单调加载情况下呈非线性变化的切线模量量化,记忆为一种分段线性分布的几何结构,并通过对这种结构加以变换,以描述循环加载情况下具有滞回特性土料的动应力-应变关系。在试验的基础上,运用非线性最小二乘法,确定了量化记忆模型中的 , , , , , 6个参数,并将量化记忆模型生成的应力-应变关系曲线与试验曲线进行了比较。结果表明,量化记忆模型是有效可行的。     

饱和砂性土流变模型的试验研究

正确合理的土体本构模型确定是地面沉降模拟的关键之一。上海和江苏地区多年来野外分层标的实际观测资料表明该区含水层系统砂层变形具有流变特征,表现出变形滞后于含水层水位变化的现象。本文基于粘弹性理论,选取了常州含水层系统原状砂性土样进行单轴压缩试验,分析了砂性土流变荷载曲线,并选择现有的相关本构模型进行识别,从而确定能更好反映研究区土体变形规律的流变模型。结果表明,研究区砂性土可采用Burgers模型进行描述,反演获得的4个流变参数随荷载的增加表现出一定的规律性。     

Model identification and parameter estimation of elastoplastic constitutive model by data assimilation using the particle filter

Data assimilation, using the particle filter and incorporating the soil‐water coupled finite element method, is applied to identify the yield function of the elastoplastic constitutive model and corresponding parameters based on the sequential measurements of hypothetical soil tests and an actual construction sequence. In the proposed framework of the inverse analysis, the unknowns are both the particular parameter within the exponential contractancy model, n_E, which parameterizes various shapes for the yield function of the competing constitutive models, including the original/the modified Cam‐Clay models and in‐between models and the parameters of the corresponding constitutive model. An appropriate set, consisting of the yield function of the constitutive model and the parameters of the constitutive model, can be simultaneously identified by the particle filter to describe the most suitable soil behavior. To examine the validity of the proposed procedure, hypothetical and actual measurements for the displacements of a soil specimen were obtained for consolidated and undrained tests through a synthetic FEM computation and for consolidated and drained tests, respectively. After examining the applicability of the proposed procedure to these test results, the present paper then focuses on the actual measured data, ie, the settlement behavior including the lateral deformation of the Kobe Airport Island constructed on reclaimed land.     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

Analytical solutions for Euler‐Bernoulli Beam on Pasternak foundation subjected to arbitrary dynamic loads

In this paper, the dynamic response of an infinite beam resting on a Pasternak foundation and subjected to arbitrary dynamic loads is developed in the form of analytical solution. The beam responses investigated are deflection, velocity, acceleration, bending moment, and shear force. The mechanical resistance of the Pasternak foundation is modeled using two parameters, that is, one accounts for soil resistance due to compressive strains in the soil and the other accounts for the resistance due to shear strains. Because the Winkler model only represents the compressive resistance of soil, comparatively, the Pasternak model is more realistic to consider shear interactions between the soil springs. The governing equation of the beam is simplified into an algebraic equation by employing integration transforms, so that the analytical solution for the dynamic response of the beam can be obtained conveniently in the frequency domain. Both inverse Laplace and inverse Fourier transforms combined with convolution theorem are applied to convert the solution into the time domain. The solutions for several special cases, such as harmonic line loads, moving line loads, and travelling loads are also discussed and numerical examples are conducted to investigate the influence of the shear modulus of foundation on the beam responses. The proposed solutions can be an effective tool for practitioners. Copyright © 2017 John Wiley & Sons, Ltd.     

Advanced numerical modeling of seismic response of a propped r.c. diaphragm wall

The paper presents some results from a number of dynamic FE simulations carried out to investigate the seismic response of a propped flexible retaining wall in a dry coarse-grained soil, considering two bedrock acceleration time histories as seismic input. Two different soil plasticity models have been considered in this study: an anisotropic hardening, critical-state model for cyclic/dynamic loading of sands and the classical Mohr–Coulomb elastic-perfectly plastic model with nonassociative flow rule. The results obtained allow to highlight the main features of the seismic performance of such type of flexible retaining structures and to evaluate the effects of the constitutive assumptions made on soil behavior on the predicted wall displacements and structural loads.     

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.     

一种土的非线性弹性本构模型参数的反演方法

旨在提出一种土的非线性弹性本构模型参数反演的方法。以现今普遍实行的地基载荷试验为基础,依据遗传算法的组合优化理论,采用正演计算和遗传算法优化相结合的方式,建立了土层非线性弹性本构模型参数反演的方法;并依据某黄土场地地基载荷试验数据,实施了黄土土层非线性弹性本构模型参数反演的全过程。计算结果表明,所建立的方法可以实现土层非线性弹性本构模型中相互关联的多个参数的组合优化,并在对初始值要求较低的情况下,可以获得良好的参数反演精度。从而为土的变形特性分析和土与其中及相邻结构的共同作用分析,提供了较好的土体本构模型参数的确定方法。     

Disturbed state model for sand–geosynthetic interfaces and application to pull-out tests

Successful numerical simulation of geosynthetic-reinforced earth structures depends on selecting proper constitutive models for soils, geosynthetics and soil–geosynthetic interfaces. Many constitutive models are available for modelling soils and geosynthetics. However, constitutive models for soil–geosynthetic interfaces which can capture most of the important characteristics of interface response are not readily available. In this paper, an elasto-plastic constitutive model based on the disturbed state concept (DSC) for geosynthetic–soil interfaces has been presented. The proposed model is capable of capturing most of the important characteristics of interface response, such as dilation, hardening and softening. The behaviour of interfaces under the direct shear test has been predicted by the model. The present model has been implemented in the finite element procedure in association with the thin-layer element. Five pull-out tests with two different geogrids have been simulated numerically using FEM. For the calibration of the constitutive models used in FEM, the standard laboratory tests used are: (1) triaxial tests for the sand, (2) direct shear tests for the interfaces and (3) axial tension tests for the geogrids. The results of the finite element simulations of pull-out tests agree well with the test data. The proposed model can be used for the stress-deformation study of geosynthetic-reinforced embankments through numerical simulation. Copyright © 1999 John Wiley & Sons, Ltd.     

粒状介质三维复杂应力加载离散元数值试验方法

实际荷载条件下（如交通、地震荷载）,粒状岩土材料常受到三维复杂应力路径作用。目前,多数粒状岩土材料的本构理论和模型都基于简单应力路径加载条件下的物理试验提出,在更加复杂应力路径下的适用性则需要进一步验证。但受机械控制的限制,物理试验中无法实现很多客观存在的三维复杂应力路径加载。为了能够再现并分析三维复杂应力路径下粒状介质的力学响应,提出了一种离散元数值试验方法,该方法采用球形数值试样,通过直接控制试样边界应力达到对3个主应力大小和方向的任意控制,从而可以实现诸多物理试验中无法实现的复杂应力路径。通过与目前常见的一些物理试验进行定性对比,论证了该数值试验方法通过高精度的加载控制和测量能够再现已有物理试验现象。在此基础上,进一步开展了应力主轴的三维旋转,分析了在这种实际存在却无法通过物理试验再现的加载条件下粒状介质的变形规律,初步显示了提出的数值试验方法在深入研究三维复杂应力路径下粒状介质力学响应方面所拥有的能力和优势。     

Quantitative prediction of discrete element models on complex loading paths

The ability of discrete element models to describe quantitatively (and not only qualitatively) the constitutive behaviour of a dense sand is assessed in this paper. Two kinds of 3D discrete models are considered. Both consider spheres as elementary particles. Nevertheless, the first model implements a contact law with rolling resistance whereas the second takes into account clumps made of two spheres. The discrete models are calibrated and validated from mechanical tests performed on a dense Hostun sand with a true triaxial apparatus. The calibration is carried out from axisymmetric drained compression tests, while the validation is discussed from monotonic and cyclic stress proportional loading paths and from a circular stress path in the deviatoric stress plane. The quality of the predictions of the discrete models are evaluated by comparison with the predictions given with advanced phenomenological constitutive relations, mainly an incrementally non-linear relation. Predictions given by the discrete models are remarkable, particularly when it is put in perspective with respect to the very few number of mechanical tests required for their calibration. However, these results and conclusions were reached in enabling conditions, and some limitations of such discrete models should be kept in mind.     

Application of an advanced multi-surface kinematic constitutive soil model

The implementation of the non-linear elastic multi-surface plastic kinematic constitutive soil model ALTERNAT into a general uncoupled finite element program called ALTICA is described. The principal features of the model are discussed and its implementation into an initial stress type excess plastic stress redistribution algorithm. To show the abilities of the model and validate the algorithm, several examples are presented including the calibration of the model to a real sand with measured monotonic and cyclic properties. In the monotonic examples, the ability of ALTICA to accurately reproduce the predicted collapse load of geotechnical structures is demonstrated using the results of finite element simulations of two typical boundary value problems (with known analytical or numerical solutions). In the cyclic examples, results of cyclic liquefaction simulation using both a non-inertial ‘static’ and inertial ‘dynamic’ algorithm are presented. Copyright © 1999 John Wiley & Sons, Ltd.     

离心场中风力电机桩基复杂加载测试设备的研制

桩基础是近海风力发电中经常采用的基础形式,承受着包括低频风循环荷载在内的复杂荷载组合。为了采用离心模型试验来研究风机桩基础特性,研制了离心场中风机桩基础复杂加载测试设备。该套设备除了可以施加水平向和竖向静力荷载外,还实现了离心场中长时间低频率的水平循环加载,并能够实现多种荷载同时或者先后施加。通过一系列测试试验,验证了该套设备的有效性,满足离心模型试验中对模拟风机桩基础的复杂荷载要求。初步试验结果表明,桩顶水平位移随着水平荷载循环周数的增加而增大,但增长速率逐渐减小,至一定周数后,桩顶位移趋于稳定,土体变形主要集中在桩周一定区域内     

Development of an inverse analysis framework for extracting dynamic soil behavior and pore pressure response from downhole array measurements

Observations from earthquakes over the past several decades have highlighted the importance of local site conditions on propagated ground motions. Downhole arrays are deployed to measure motions at the ground surface and within the soil profile, and also to record the pore pressure response within the soft soil profiles during excitation. The degradation of soil stiffness as excess pore pressures are generated during earthquake events has also been observed. An inverse analysis framework is developed and demonstrated to directly extract soil material behavior including pore water pressure (PWP) generation from downhole array measurements that can then be readily used in 1D nonlinear site response analysis. The self‐learning simulations (SelfSim) inverse analysis framework, previously developed for total stress site response analysis, is extended to extract PWP generation behavior of the soil in addition to cyclic response during ground shaking. A Neural Network based constitutive model is introduced to represent PWP generation during cyclic loading. A new analysis scheme is introduced that can use data from co‐located piezometer and accelerometer sensors. The successful performance of the proposed framework is demonstrated using four synthetic vertical array recordings. Copyright © 2013 John Wiley & Sons, Ltd.     

Single‐and multi‐objective genetic algorithm optimization for identifying soil parameters

This paper discusses the quality of the procedure employed in identifying soil parameters by inverse analysis. This procedure includes a FEM‐simulation for which two constitutive models—a linear elastic perfectly plastic Mohr–Coulomb model and a strain‐hardening elasto‐plastic model—are successively considered. Two kinds of optimization algorithms have been used: a deterministic simplex method and a stochastic genetic method. The soil data come from the results of two pressuremeter tests, complemented by triaxial and resonant column testing. First, the inverse analysis has been performed separately on each pressuremeter test. The genetic method presents the advantage of providing a collection of satisfactory solutions, among which a geotechnical engineer has to choose the optimal one based on his scientific background and/or additional analyses based on further experimental test results. This advantage is enhanced when all the constitutive parameters sensitive to the considered problem have to be identified without restrictions in the search space. Second, the experimental values of the two pressuremeter tests have been processed simultaneously, so that the inverse analysis becomes a multi‐objective optimization problem. The genetic method allows the user to choose the most suitable parameter set according to the Pareto frontier and to guarantee the coherence between the tests. The sets of optimized parameters obtained from inverse analyses are then used to calculate the response of a spread footing, which is part of a predictive benchmark. The numerical results with respect to both the constitutive models and the inverse analysis procedure are discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

软粘土地基上循环承载力的计算模型研究

分析了海洋平台基础受到静荷载和循环荷载作用时地基土体的应力状态,对循环荷载作用下计算软土地基的承载力的计算方法进行了研究.在循环三轴试验得出的计算土单元循环强度模型的基础上,提出一种计算软土地基循环承载力的拟静力弹塑性模型,以桶形基础为例,对软粘土地基的循环承载力进行了计算分析.     

A numerical hydro-geotechnical model for marine gravity structures

A new finite-volume solver named “geotechFoam” is developed within OpenFOAM®¹ for modelling soil-structure-interaction for marine gravity structures. The fully-coupled and fully-dynamic Biot’s governing equations are solved in a segregated approach. Two simplifications of the governing equations are introduced and tested. The spatial domain is composed of several zones with different material properties and/or constitutive models. A multi-yield surface plasticity model is implemented to simulate soil response under cyclic loads. Moreover, soil-structure interaction is modelled via a frictional contact model and boundary conditions accounting for skeleton-pore fluid coupling. Five benchmark cases and two sets of physical model tests are applied for the validation.     

Constitutive modeling of dense gravelly soils subjected to cyclic loading

A constitutive model for dense gravelly soils was developed to reproduce their responses under cyclic loadings. Its application aims at nonlinear dynamic analyses of earth structures involving gravelly soils, such as rockfill dams and railroad ballasts. The framework of generalized plasticity was modified to incorporate the concept of stress distance for better simulation of unloading and reloading responses. It was then combined with the theory of critical state soil mechanics to develop the constitutive model. The model has the following important features: unified simulation of particle breakage through translating critical state line, smooth transition from unloading to reloading in the stress space, and proper modeling of cyclic hysteresis, cyclic densification, and cyclic hardening of dense gravelly soils. Most of the model parameters can be obtained through simple calculation using conventional triaxial test results, and their calibration process was discussed. The model was used to simulate the cyclic responses of three gravelly soils with satisfactory accuracy. Copyright © 2014 John Wiley & Sons, Ltd.