A general approach to model interfaces using existing soil constitutive models application to hypoplasticity

The modelling of interfaces is important for the holistic simulation of geotechnical structures (e.g. piles, tunnels and geogrids). For this reason, advanced constitutive interface models and numerical techniques are needed. There are few user-friendly models, and these are rarely implemented. In this paper, a new approach for advanced interface models is proposed. This is based on the assumption that the fully rough interface can be modelled considering simple-shear behaviour at the interface. A 3D soil model is used as a constitutive driver for a frictional subroutine. This minimises the effort required, and advanced interface models are available with less effort. Two different hypoplastic models are used with the new approach. The approach was verified for several aspects (e.g. mesh size dependence), and the volumetric behaviour was studied. The user-friendliness and absence of additional parameters led to more realistic simulation results. The proposed method can be extended to other modelling techniques and will improve the modelling of contacts in soil-structure interaction analysis.     

新时代岩土力学基本问题探究北大核心CSCD

“一带一路”新时代背景下大规模的工程建设,推动着岩土工程学科的繁荣和蓬勃发展。土动力学、非饱和土与特殊土力学与环境岩土工程是未来土力学的重点研究方向,有必要进行系统梳理分析。本文从土工动力测试技术、动本构关系和土-结构动力相互作用方面对土动力学的主要内容进行了阐释,然后从非饱和土有效应力原理、非饱和土本构模型及黄土力学特性方面对非饱和土与特殊土力学的基本理论进行了分析,最后对垃圾填埋处置和水土污染等环境岩土工程亟待解决的问题进行了探索,并提出了合理的应对措施。土力学学科面临着极大的发展机遇,也面临着严峻的挑战,这一系列核心问题的解决,必将为全面提升岩土工程项目建设的质量和效益作出贡献。     

基于Hardin曲线的土体边界面本构模型在ADINA软件中的实现

为了提出一种适合于岩土地震数值模拟的土体本构模型,基于土体动应力-应变关系的Hardin曲线及其在非等幅往返荷载下的Pyke修正,采用von Mises准则在偏应力平面上构造边界面,以反向加载点和当前应力点的连线在边界面上投影的比例作为硬化参数,推导了塑性硬化模量并给出该边界面本构的具体增量表述。在有限元软件ADINA中通过自定义材料的二次开发实现了该本构模型,并利用动三轴试验对该本构模型进行了验证。数值模拟与试验结果的对比表明,本构模型能如实反映土体的应力-应变关系。针对实际工程场地的地震反应,应用边界面本构模型在ADINA中进行了二维数值模拟,与SHAKE91的计算结果进行了对比,说明了该本构模型应用于岩土地震工程问题的合理性。     

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

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

Physical and Numerical Modeling of Seismic Soil-Structure Interaction in Layered Soils

The structural response of buildings subjected to seismic loads is affected by local site conditions and the interaction between the structure and the supporting soil media. Seismic centrifuge model tests were conducted on two layered clay soil profiles at 80 g field to investigate soil-structure interaction and dynamic response of foundation. Several earthquake-like shaking events were applied to the models using an electro-hydraulic shaking table to simulate linear and nonlinear soil behavior. Results showed that the foundation input motion was significantly amplified in both models, especially for weak earthquake motions. Seismic soil-structure interaction was found to have an important effect on structure response by increasing the amplification of foundation input motion. A 3D finite difference numerical model was also developed to simulate the response of centrifuge model tests and study the parameters that affect the characteristics of earthquake at the base of the structure. The results indicated that the stiffness and stratification of the soil profiles had a significant effect on modifying the foundation input motion.     

A bounding surface model for saturated and unsaturated soil-structure interfaces

In many geotechnical systems, such as reinforced slopes and embankments, soil-structure interfaces are often unsaturated. Shear behaviour of unsaturated interfaces is strongly dependent on their matric suctions, as revealed by the results of extensive laboratory tests. So far, constitutive models for unsaturated interfaces are very limited in the literature. This paper reports a new bounding surface model for saturated and unsaturated interfaces. New formulations were developed to incorporate suction effects on the flow rule and plastic modulus. To examine the capability of the proposed model, it was applied to simulate suction- and stress-controlled direct shear tests on unsaturated soil–cement, soil–steel and soil–geotextile interfaces. Measured and computed results are well matched, demonstrating that the proposed model can well capture key features of the shear behaviour of unsaturated interfaces, including suction-dependent dilatancy, stress–strain relation and peak and critical state shear strengths.     

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.     

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

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

基于广义位势理论的接触面本构模型一般表述

基于广义位势理论建立的岩土体材料本构模型以及岩土体材料与结构接触面本构模型原理相通,只是前者是在三轴剪切试验条件下的三维应力空间建模,后者是在单剪试验条件下的二维应力空间建模。单剪试验条件下土与结构的接触面问题可以看作是法向与切向应力空间上的二维问题,其试验结果可以表达成由应力、应变组成的二维矢量。结合接触面力学特性,确定应力空间中的势函数以及塑性状态方程,可以推导出双重势面接触面弹塑性本构方程的一般表达式。进一步取两个势函数分别为法向应力和切向应力,建立简化双重势面接触面弹塑性模型的本构方程,该方程可直接应用于有限元等数值分析。结合试验实例对建模方法的合理性进行验证,模型拟合效果良好。研究结果表明,基于广义塑性位势理论建立接触面本构模型无需推求塑性势函数和屈服函数,可以直接得到弹塑性刚度矩阵,且建模方便。     

Review of constitutive models for describing the time dependent behaviour of soft clays

Different classes of constitutive models have been proposed to capture the time-dependent behaviour of soft soil (creep, stress relaxation, rate dependency). This paper critically reviews many of the models developed based on understanding of the time dependent stress-strain-stress rate-strain rate behaviour of soils and viscoplasticity in terms of their strengths and weaknesses. Some discussion is also made on the numerical implementation aspects of these models. Typical findings from numerical analyses of geotechnical structures constructed on soft soils are also discussed. The general elastic viscoplastic (EVP) models can roughly be divided into two categories: models based on the concept of overstress and models based on non-stationary flow surface theory. Although general in structure, both categories have their own strengths and shortcomings. This review indicates that EVP analysis is yet to be vastly used by the geotechnical engineers, apparently due to the mathematical complication involved in the formulation of the constitutive models, unconvincing benefit in terms of the accuracy of performance prediction, requirement of additional soil parameter(s), difficulties in determining them, and the necessity of excessive computing resources and time.     

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

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

动力分析中无限地基两种近似方法的应用研究

在结构-地基的动力相互作用分析所提出的众多方法中,无质量地基和黏性边界地基仍然是应用最为广泛的近似方法。针对当前结构-地基动力相互作用分析中地基模型处理中的一些问题,通过结构-地基动力相互作用分析公式的推导以及各变量物理意义的解释,详细阐明了当前考虑地基影响的动力分析中无质量地基和黏性边界地基的概念和应用范围,对当前部分研究中应用这两种地基中存在的问题进行了分析说明,为简化地基动力分析模型提供了思路。     

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.     

土-结构-流体动力相互作用的实时耦联动力试验

针对振动台试验中无限地基难以模拟和数值分析中流-固耦合作用难以计算两个难题,将最近发展的实时耦联动力试验方法引入土-结构-流体动力相互作用问题的研究。以一个渡槽结构为例,其中渡槽-水体作为物理子结构,采用振动台进行物理试验,而无限地基作为数值子结构,采用集总参数模型进行数值模拟。两个子结构之间实时交换数据,联合评估整个耦合体系的动力响应。试验结果和有限元数值模拟结果吻合良好,表明该试验方法具有较高精度。对不同特性地基土进行的试验对比分析结果表明：对于软土地基,考虑土-结构相互作用(SSI)的结构反应幅值明显减小,周期延长;随着地基土变硬,SSI效应逐渐变弱,结构反应最终收敛至刚性地基解。     

Determination of failure probabilities of flood defence systems with improved dynamic bounds method

In this paper, we describe the computational framework of a novel method for solving the challenging problem of probabilistic finite elements. The method is called Improved Dynamic Bounds (IDB) and was developed recently to improve the efficiency of the dynamic bounds. The IDB is used in finite element numerical models to calculate time-dependent failure analyses of structures. In applications, the IDB can speed up the overall simulation process by several orders of magnitude. In applications controlled by two influential variables (e.g, two-dimensional problem), the computational efficiency is improved by a factor of 769 according to Rajabalinejad (2009). Applications of IDB indicate the method is most efficient for problems where the number of influential variables are limited. This is often the case for geotechnical and coastal flood defence systems. The IDB method is applied in this paper to the 17th Street Flood Wall, a component of the flood defence system (levee infra-structure) that failed during the Hurricane Katrina, to calculate the failure probability of an I-wall.     

基于耦合分析法的地铁隧道抗震研究

在土体与结构动力相互作用理论的基础上,提出用有限元和无限元耦合分析法对地铁隧道结构进行抗震分析。为考虑相互作用的非线性特点,选用Ramberg-Osgood骨架曲线来模拟土体本构关系的非线性特性,同时选用GOODMAN单元来模拟土与结构接触面上的非线性接触状态。以双洞口矩形隧道结构为例进行了动力反应的分析与比较,得到了一些有益的规律和结论。     

Experimental Assessment of the Stress–Strain Behaviour of Leighton Buzzard Sand for the Calibration of a Constitutive Model

A number of constitutive models are nowadays implemented in numerical codes which simulate the stress–strain behaviour of soil from very small to large strain. In this paper, the mechanical behaviour of Leighton Buzzard sand (grade E), used worldwide for physical modelling, has been thoroughly characterized by laboratory testing along several stress paths. Tests were aimed at calibrating a constitutive model, that allows considering stiffness nonlinearities in a wide range of strains, in the framework of isotropically hardening plasticity. As a validation, the results of dynamic centrifuge tests on a layer of the same sand were compared with finite element predictions.     

海岸岩土工程的物理与数值模拟方法

采用物理模型试验来再现实际现象,验证数值模拟方法的有效性;同时采用数值模拟来弥补离心模型试验的缺陷,揭示模型试验规律,评价实际海岸岩土-结构系统在不同工况下的变形与强度稳定性及其变化.并简要地介绍了针对防波堤土工织物加筋垫层作用机理研究、黄骅港外航道防沙堤工程安全性评价、上海洋山港浅置式大圆筒结构试验段工程安全性评价三个具体项目进行评价的成果,初步展现了物理模型试验与数值模拟方法在海岸岩土工程的机理研究、方案论证、结构与施工优化、结构安全性评价等方面所能发挥的重要作用和广阔的应用前景.     

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.     

On the time-dependent behaviour of glacial sediments: a geotechnical approach

Soft sediments can deform under both changes in stress and changes in strain rate, as well as during sustained loading (e.g. creep). In glacial conditions, soft subglacial sediments can contribute to glacier movements, and when the ice/sediment system is coupled, the rate-sensitivity of the sub-layer can affect the velocity of the glacier. These issues have been the object of studies by several geologists. In parallel, the effects of time and strain rate on the behaviour of soils have been studied for many years by geotechnical engineers, to estimate the long-term performance of ground structures. This paper applies the knowledge acquired in soil mechanics to a geological problem: results on two glacial sediments of different origins, obtained from advanced geotechnical laboratory tests, are presented and analysed in the light of recent advances in soil mechanics. Independently of time effects, the test data show that the pre-failure behaviour is not elastic but that the sediments develop plastic strains before reaching their peak strength, with the stiffness gradually degrading. Results from tests performed with variable strain rates suggest that the time-dependent (viscous) behaviour of the sediments may be governed by their mineralogy, in particular the presence of plastic fines. The viscous behaviour observed in the plastic tills seems to affect the pre-failure behaviour only, for example the yield surface, stiffness, strength, but the criteria of failure, e.g. critical or residual angle of failure are found to be constant and independent of strain rate. The test data also suggest that the viscous behaviour may be related to stress level, being more prevalent at low stresses. These results, obtained under laboratory-controlled conditions, over a large range of strains, pre- and post-peak strength, are a necessary step towards developing constitutive models for subglacial sediments rigorously. It is only after this has been achieved that such models should be implemented in finite element code to analyse ice/sediment systems, with the aim to understand glacier movements better.