Numerical modeling and neural networks to identify model parameters from piezocone tests: II. Multi‐parameter identification from piezocone data

This paper completes the study presented in the accompanying paper, and demonstrates a numerical algorithm for parameter prediction from the piezocone test (CPTU) data. This part deals with a development of neural network (NN) models which are able to map multi‐variable input data onto typical geotechnical characteristics and constitutive parameters of the modified Cam clay model, which has been applied in this study. The identification procedure is designed for the coupled hydro‐mechanical boundary value problem in normally‐and lightly overconsolidated clayey soils including partially drained conditions that may appear during cone penetration. The NN models are trained with pseudo‐experimental measurements derived with the aid of the numerical model of the piezocone test, presented in the accompanying paper. Different input configurations containing CPTU measurements and some complementary data are studied with respect to the accuracy of predicted parameter values. Finally, the performance of the developed NN predictors is tested with field CPTU data which are derived from three well‐documented characterization sites, and the obtained predictions are compared with benchmark laboratory results. Copyright © 2011 John Wiley & Sons, Ltd.     

Identifying parameters controlling soil delayed behaviour from laboratory and in situ pressuremeter testing

The aim of this paper is to present a methodology for identifying the soil parameters controlling the delayed behaviour from laboratory and in situ pressuremeter tests by using an elasto‐viscoplastic model (EVP‐MCC) based on Perzyna's overstress theory and on the elasto‐plastic Modified Cam Clay model. The influence of both the model parameters and the soil permeability was studied under the loading condition of pressuremeter tests by coupling the proposed model equations with Biot's consolidation theory. On the basis of the parametric study, a methodology for identifying model parameters and soil permeability by inverse analysis from three levels of constant strain rate pressuremeter tests was then proposed and applied on tests performed on natural Saint‐Herblain clay. The methodology was validated by comparing the optimized values of soil parameters and the values of the same parameters obtained from laboratory test results, and also by using the identified parameters to simulate other tests on the same samples. The analysis of the drainage condition and the strain rate effect during a pressuremeter test demonstrated the coupled influence of consolidation and viscous effects on the test results. The numerical results also showed that the inverse analysis procedure could successfully determine the parameters controlling the time‐dependent soil behaviour. Copyright © 2008 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.     

Analysis of pressuremeter geometry effects in clay using critical state models

The conventional interpretation methods of pressuremeter testing effectively approximate pressuremeter membranes as infinitely long. As a result, the effects of the two‐dimensional geometry of pressuremeters are ignored, leading to an overestimation of soil shear strength by pressuremeter testing, as demonstrated in several previous studies. This paper presents results of a numerical study of two‐dimensional geometry effects on self‐boring pressuremeter tests in undrained clay. The results are obtained using critical state soil models with an effective stress formulation. This is in contrast to most (if not all) existing studies on pressuremeter geometry effects, which were based on perfectly plastic soil models (e.g. Yu (Cavity expansion theory and its application to the analysis of pressuremeters. DPhil Thesis, The University of Oxford, 1990), Yeung and Carter (Proc. 3rd Int. Symp. on Pressuremeters, 1990), and Houlsby and Carter (Géotechnique, 1993; 43 (4):567–576)). The present study suggests that the overestimation of soil strength due to the neglect of finite pressuremeter length is significantly affected by the soil model used in the calculations. It is found that for clays with a high overconsolidation ratio (OCR) the strength overestimation predicted using critical state soil models could be considerably smaller than that predicted using perfectly plastic soil models. The main conclusion of this numerical study is that care must be exercised before directly applying any numerically determined pressuremeter geometry correction factors in practice. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Determining soil permeability from pressuremeter tests

This paper presents a methodology for identifying soil permeability from pressuremeter test. On the first part we present a numerical analysis of the permeability effects on the test results. We demonstrate that different drainage conditions arise during test, as a function of the loading rate and the soil permeability. We also studied the pore pressure dissipation during strain holding stages. Based on this analysis of these tests, we propose a general procedure to identify simultaneously mechanical parameters and permeability from pressuremeter tests with strain holding test stages. This procedure was applied on tests performed on natural Saint–Herblain clay. An apparatus called pressio‐triax was developed for this purpose. The values of the mechanical parameters as well as of the permeability value were found to agree very well with the values of the same parameters obtained from conventional laboratory tests. Copyright © 2002 John Wiley & Sons, Ltd.     

利用旁压试验的砂土初始状态反演分析

讨论了砂土旁压试验的反分析问题。建议利用旁压试验实测数据反演确定土体的初始状态而非材料参数,材料参数可以通过实验室常规试验确定。主要利用了新型的亚塑性本构模型参数不受土体状态影响的特点。主要内容包括：建立有限元数值模型,利用文献中的试验数据确定了特定砂土的材料参数;通过模拟砂土的实验室旁压试验,验证了数值模型;通过改变土体初始应力状态和相对密实度进行数值计算,并利用数值计算成果,建立了可适用于砂土初始状态反演分析的实测旁压荷载和土体中应力与相对密度的数学关系。     

Time homogenization for clays subjected to large numbers of cycles

This paper discusses the reliability and the efficiency of a time homogenization method employed to reduce the computational time during cyclic loading for two common geotechnical tests and two elastoplastic models for clays. The method of homogenization is based upon splitting time into two separate scales. The first time scale relates to the period of cyclic loading and the second to the characteristic time of the fatigue phenomenon. The time homogenization method is applied to simulate an undrained triaxial test (homogeneous stress state) and a pressuremeter test (nonhomogeneous stress state) under one‐way cyclic loading on normally consolidated clay. This method is coupled with two elastoplastic models dedicated to cyclic behavior of clay (a bounding surface plasticity model and a bubble model). Both linear and nonlinear elasticities are considered. The difficulty encountered when applying this method to models introducing nonlinear elasticity and kinematic hardening is pointed out. The performance of time homogenization related to the main parameters is numerically investigated by comparison with conventional finite element simulations. Copyright © 2012 John Wiley & Sons, Ltd.     

Optimization framework for calibration of constitutive models enhanced by neural networks

A two‐level procedure designed for the estimation of constitutive model parameters is presented in this paper. The neural network (NN) approach at the first level is applied to achieve the first approximation of parameters. This technique is used to avoid potential pitfalls related to the conventional gradient‐based optimization techniques, considered here as a corrector that improves predicted parameters. The feed‐forward NN (FFNN) and the modified Gauss–Newton algorithms are briefly presented. The proposed framework is verified for the elasto‐plastic modified Cam Clay model that can be calibrated based on standard triaxial laboratory tests, i.e. the isotropic consolidation test and the drained compression test. Two different formulations of the input data to the NN, enhanced by a dimensional reduction of experimental data using principal component analysis, are presented. The determination of model characteristics is demonstrated, first on numerical pseudo‐experiments and then on the experimental data. The efficiency of the proposed approach by means of accuracy and computational effort is also discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

SBPT测定饱和黏土不排水强度的数值分析

自钻式旁压试验（SBPT）因其扰动小、测试深度大、可以获得应力-应变、超孔隙水压力-时间等数据,在确定地基土性参数和地基承载力上有广阔的应用前景。然而由于用以解释SBPT的柱孔扩张理论（Gibson解）所采用的平面应变假设与实际旁压腔几何特征存在差异,导致试验所确定的黏土不排水剪切强度su与其他原位试验或室内试验结果存在差别。针对旁压腔几何尺寸及应变区间的选择对确定su的影响,基于修正剑桥模型,采用低渗透系数控制加载过程中不排水条件,利用有限元法模拟SBPT,建议了不同应力历史下确定su的应变区间,并给出考虑几何尺寸影响时相应应变区间上su的修正系数。     

不同固结条件下黏土数值化建模的研究

研究了不同初始固结压力条件对黏土本构关系的影响。进行了2组不同固结条件下的排水压缩试验,即等压固结条件和K0固结条件。采用数值建模方法,以 和 为硬化参数的双屈服面模型为框架,通过对试验数据空间的反演获得两种不同初始条件下的黏土弹塑性本构关系,并在整个应力场 可视化。将该数值模型嵌入有限元程序对三轴压缩试验进行数值模拟,并与试验曲线及相同路径下的剑桥模型和Duncan-Chang模型关系曲线进行对比。结果表明,数值建模方法有较好的容错性,能很好地反映应力历史的影响,能模拟真实应力路径。     

Optimization techniques for identifying soil parameters in geotechnical engineering: Comparative study and enhancement

A comparative study of optimization techniques for identifying soil parameters in geotechnical engineering was first presented. The identification methodology with its 3 main parts, error function, search strategy, and identification procedure, was introduced and summarized. Then, current optimization methods were reviewed and classified into 3 categories with an introduction to their basic principles and applications in geotechnical engineering. A comparative study on the identification of model parameters from a synthetic pressuremeter and an excavation tests was then performed by using 5 among the mostly common optimization methods, including genetic algorithms, particle swarm optimization, simulated annealing, the differential evolution algorithm and the artificial bee colony algorithm. The results demonstrated that the differential evolution had the strongest search ability but the slowest convergence speed. All the selected methods could reach approximate solutions with very small objective errors, but these solutions were different from the preset parameters. To improve the identification performance, an enhanced algorithm was developed by implementing the Nelder‐Mead simplex method in a differential algorithm to accelerate the convergence speed with strong reliable search ability. The performance of the enhanced optimization algorithm was finally highlighted by identifying the Mohr‐Coulomb parameters from the 2 same synthetic cases and from 2 real pressuremeter tests in sand, and ANICREEP parameters from 2 real pressuremeter tests in soft clay.     

A coupled analysis of cavity and pore volume changes for pulse tests conducted in soft clay deposits

This paper presents a new interpretation method for pulse tests, a field permeability test that allows for rapid measurements of hydraulic conductivity k in aquitards. This new method applies to soft clay deposits. To initiate a pulse test, a known volume of water is injected into the sand filter of a monitoring well isolated by a packer. The resulting pressure increase yields an outward movement of the sand filter cavity wall. After presenting the usual interpretation methods and their limits, this paper proposes a new interpretation method based on a coupled analysis of the pressure and displacement fields in the soil using the Biot–Darcy formulation. A series of analytical and numerical non‐dimensional velocity graphs, normalized plots of the mean hydraulic head difference versus its rate of change, are given. For a linear elastic material, these type curves show relatively small variations with the sand filter aspect ratio and the Poisson ratio of the tested clay. The type curves are also found to be independent of the clay compressibility (m_v) and k, an important result. A series of pulse tests conducted in a soft marine clay deposit near Montreal, Canada, are interpreted with the proposed method. The hydraulic conductivity values calculated from these tests are closely correlated with independent estimates obtained using long‐duration variable‐head tests. Compared with previous interpretation methods, the proposed method allows soil volume changes to be reconciled with cavity expansion phenomena and the range of type curves available for the interpretation of test data to be constrained. Copyright © 2013 John Wiley & Sons, Ltd.     

A rigorous semi‐analytical solution for undrained cylindrical cavity expansion in critical state soils

This paper presents a generalized, rigorous and simple large strain solution for the undrained expansion of a vertical cylindrical cavity in critical state soils using a rate‐based plasticity formulation: the initial stress field is taken as anisotropic, that is with horizontal stresses that differ from the vertical stress, and the soil is assumed to satisfy any two‐invariant constitutive model from the critical state (Cam‐clay) family; no simplifying assumption is made during the mathematical derivation; calculating the effective stresses around the cavity requires the solution of a nonlinear equation by means of the Newton–Raphson method in combination with quadrature. Cavity expansion curves and stress distributions in the soil are then presented for different critical state models (including the modified Cam‐clay model). The solution derived can be useful for estimating the instantaneous response of saturated low‐permeability soils around piles and self‐boring pressuremeters and can serve as trustworthy benchmark for numerical analysis codes. Copyright © 2016 John Wiley & Sons, Ltd.     

In Situ Pore Water Pressure Dissipation Testing of Marine Clay under Reclamation Fills

In situ dissipation tests provide a means of evaluating the in situ coefficient of horizontal consolidation and horizontal hydraulic conductivity of soft clays. Dissipation tests by means of piezocone (CPTU), dilatometer (DMT), self-boring pressuremeter (SBPT) and BAT permeameter (BAT) were utilized in the characterization of the coefficient of horizontal consolidation and horizontal hydraulic conductivity of Singapore marine clay at Changi in a land reclamation project. Dissipation tests were carried out prior to reclamation as well as after ground improvement with vertical drains to compare the changes in the coefficient of horizontal consolidation and horizontal hydraulic conductivity prior to and after ground improvement. Tests were carried out in a vertical drain area as well as in an adjacent untreated control area after 23 months of surcharge loading, for comparison purposes. The purpose of this research is to determine the horizontal consolidation parameters of Singapore marine clay prior to reclamation as well as after 23 months of surcharge loading with and without vertical drains by means of in situ dissipation tests.     

Computer-aided calibration of a soil plasticity model

A non-linear optimization technique based on the quasi-Newton approach is employed to back-calculate certain model parameters of a simple, bounding surface, soil plasticity model from in situ pressuremeter data. The theoretical response corresponding to a given set of parameters is generated by finite element analysis. A semi-analytical procedure is developed for the accurate and efficient evaluation of the gradient of objective function with respect to the model parameters of interest. The BFGS update is used to update the Hessian. Results of a series of numerical experimentation using artificial pressuremeter responses is first reported and discussed. A set of laboratory cavity expansion data is then used to calibrate the constitutive model.     

Effects of transport of pore water and material heterogeneity on strain localization of fluid‐saturated gradient‐dependent viscoplastic geomaterial

The purpose of the present paper is to clarify the effects of permeability and initial heterogeneity on the strain localization of fluid‐saturated cohesive soil modelled by a strain gradient‐dependent poro‐viscoplastic constitutive model. The effects of permeability and gradient parameters on the growth rate of the fluctuation were obtained by a linear instability analysis. Deformation behaviour of clay specimens modelled as a viscoplastic model with a second order strain gradient during shear was numerically analysed by a soil–water coupled FEM under both globally undrained and partially drained conditions. It was found that the deformation pattern and the stress–strain curve greatly depend on the permeability, the drainage conditions and the initial non‐homogeneous properties. Copyright © 2005 John Wiley & Sons, Ltd.     

一种适用于饱和黏土循环动力分析边界面塑性模型的隐式积分算法

修正了传统隐式回映算法,建立了适用于饱和黏土循环动力分析的边界面塑性模型的完全隐式积分格式。该模型基于无弹性域概念和临界状态理论,采用各向同性、运动硬化准则、旋转的边界面,并引入表征土体结构损伤和重塑程度的损伤变量以反映循环载荷作用下饱和黏土的各向异性、刚度、强度软化及塑性变形累积等特征。针对等压固结 和偏压固结 的饱和高岭黏土的不排水三轴试验进行模拟,采用不同的应变增量步长进行计算,并与试验数据对比,结果表明,修正隐式回映算法应用于该类边界面模型的合理性、积分格式的精确性和稳定性;另外,结合有限元软件自动时间步长的增量迭代解法,对饱和黏土应力控制的不排水动三轴试验进行预测,结果表明,修正的适用于该边界面的塑性模型隐式回映算法可以得到比较合理的数值分析结果,能够反映饱和黏土的循环刚度的退化和强度的弱化等动力特性。     

Analysis of undrained cavity expansion in elasto‐plastic soils with non‐linear elasticity

A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.