土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.      

经典弹塑性力学体系下岩土材料的本构方程

剑桥模型,修正剑桥模型以及近年来出现的适用于砂土和黏土的统一硬化模型在岩土工程中有重要的科研和实用价值．在经典弹塑性力学框架下重新推导了以上这些岩土材料本构方程,对岩土工程中常用的本构模型和经典弹塑性理论框架下的岩土本构模型进行了对比,将两个体系下的相关量建立了联系．     

动力UH模型及其有限元应用

饱和砂土在循环载荷下具有复杂的应力应变关系, 通常表现出液化过程中的大变形以及往返活动性现象. 为简单有效地模拟上述特性, 在超固结UH模型的基础上, 将其扩展为可考虑砂土动力加载下的本构模型. 具体做法有3点: (1)改变屈服面椭圆长短轴之比, 将比值定义为反映应力诱导各向异性转轴斜率的函数; (2)引入旋转硬化规则, 用来反映应力诱导各向异性; (3)建立一个与旋转硬化规则以及临界状态特性相协调的统一硬化参数. 模型预测结果表明, 所提动力模型可简单、有效地用于砂土在动力载荷下应力应变关系的模拟. 最后将该动力UH模型嵌入到有限元软件中, 三维地基的动力加载模拟结果表明, 动力UH模型可方便地应用于岩土工程实践中.      

含主应力轴旋转的土体本构关系研究进展

主应力轴旋转对岩土工程的影响日益受到人们的重视．本文从纯主应力轴旋转、纯应力洛德角变化、多种因素变化、排水与不排水等方面较全面地描述了含主应力轴旋转情况下土体的基本变形特性．对当前较有影响的含主应力轴旋转的土体本构模型（基于一般应力空间的土体本构模型、运动硬化模型、边界面模型、土体弹塑性应力应变关系的完全应力增量表述等等）进行了较为系统的评述．提出了合主应力轴旋转的实验研究中存在的核心问题与建立合主应力轴旋转的土体本构关系的根本途径．      

土的应力-应变关系的一种描述模式

通常将由土的剪切试验测得的应力-应变关系曲线(q～ε1曲线)分为应变硬化型和应变软化型,文中提出了一种能同时描述应变硬化型q～ε1曲线和应变软化型q～ε1曲线的新模式,并导出了统一的切线模量表达式,进一步探讨了该应力-应变关系曲线描述模式在拟合应变硬化型曲线时的简化形式,及简化形式中参数取值方法和参数与围压的关系等方面的问题。通过与邓肯-张模型和应变软化模型的对比,结果表明该应力-应变关系曲线描述模式能更好地与试验数据吻合,且用其简化形式拟合应变硬化型曲线时,可以通过调整其中一个参数值来表达出不同的曲线形式,从而体现出土体应力-应变关系的多样性。该应力-应变关系描述模式为发展更一般的非线性弹性模型提供了一定的理论基础。     

正常固结饱和黏性土三剪有限变形弹塑性本构模型研究

基于三剪统一强度准则,利用等量代换法和坐标平移法分别推导出正常固结饱和黏性土的三剪破坏应力比,并将其与修正剑桥模型相结合,得到三剪统一屈服面方程。为反映饱和黏性土的变形非线性及大变形特性,基于有限变形理论,建立了正常固结饱和黏性土的三剪统一有限变形等量代换法和坐标平移法2种弹塑性本构模型。为验证模型的适用性,取江西红黏土制备三种不同压实度的试样,在不同围压下进行常规三轴固结不排水和固结排水压缩试验。将试验数据与本文所提2种有限变形模型及相应的2种小变形模型计算结果进行了对比分析,结果表明,随着变形的发展,有限变形模型相对于小变形模型更接近试验结果,能较好地反映黏性土因高孔隙率（低初始压实度、小围压及压缩变形前期）而产生的大变形特性,虽然小变形模型的计算偏差会随着土样初始压实度和所施加围压的增大而减小,但有限变形模型对不同压实度和围压的计算偏差均相对较小,其中,等量代换法有限变形模型在初始压实度较大或者围压较高时所得计算偏差相对最小。对所提本构模型所做真三轴分析表明,中间主应力影响系数b和初始压实度对土体的强度和变形特性具有一定的影响。主应力差、孔隙水压力和体应变与b值呈正相关性,主应力差与初始压实度呈正相关性,孔隙水压力和体应变与初始压实度呈负相关性。     

Modelling of off-road wheeled vehicles for real-time dynamic simulation

Simulation of wheel-ground and vehicle-ground interactions is very important in many applications. Achieving accuracy and efficiency is challenging for both soft and hard terrains. This is not only because of the simulation and numerical challenges, but also due to the questionable nature of the existing terrain models. For example, the most widely used terramechanics model is not a representative constitutive relation for a full range of dynamic conditions and applications, but rather a parametrization of steady state conditions. In general, the selection and development of the proper constitutive model and the parametrization of the ground properties are very challenging. Here, we present a unified framework for general wheel-ground interaction which can be used with different terramechanics models. The framework is based on a complementarity formulation and also uses the concept of kinematic constitutive relations, beside the other known concepts for modelling and parametrizing the soil properties. The framework makes it possible to consider the appropriate modelling of the terrain for a broad range of dynamic behaviours and simulation conditions. We will illustrate the material with several examples for off-road conditions.     

Modification of critical state models by Mohr–Coulomb criterion

In order to combine the Mohr–Coulomb criterion and the critical state models for soils, a transformed stress tensor based on the Mohr–Coulomb criterion is proposed. The new stress tensor is determined by a transformation that makes the Mohr–Coulomb criterion become a cone having an axis as the space diagonal in transformed principal stress space, and is applied to the Cam-clay and Sekiguchi–Ohta models, which are two typical isotropic and anisotropic hardening critical state models for soils, for improving modelling capability of describing soil behaviour in general stresses. The revised models give more accurately predicted results than the original ones.     

泡沫铝本构行为研究进展

作为泡沫金属的典型代表, 泡沫铝是一种集结构性和功能性于一体的新型材料, 随着生产工艺的提高及国民经济的发展, 泡沫铝在航空航天、交通运输、建筑工程、机械制造等领域的应用日趋深入和广泛. 工程中复杂的载荷工况, 对泡沫铝本构模型的准确性及实用性提出了更高的要求. 国内外学者对泡沫铝进行了大量的试验研究及有限元数值分析工作, 随着各种研究和数值分析本构模型的提出、验证及修正等, 人们对泡沫铝的力学性能及其应用的认识逐渐提高. 本文首先简要概述了泡沫铝力学性能的试验研究及有限元数值分析的进展, 然后重点综述泡沫铝本构行为的国内外研究进展及现状. 最后, 针对目前存在的问题, 讨论并展望了泡沫铝本构行为的研究发展趋势. 完善本构模型表征所需的特征参量, 将各向异性或横观各项同性材料假定引入模型构建体系, 明确静水压缩响应及单轴压缩响应在材料强化中的权重, 建立例如随动强化般更加真实且准确反映泡沫铝强化过程的强化准则, 将应变率效应研究成果引入本构模型中等是现有泡沫铝本构模型体系中的重要研究方向.      

Investigation of advanced strain-path dependent material models for sheet metal forming simulations

Sheet metal forming processes often involve complex loading sequences. To improve the prediction of some undesirable phenomena, such as springback, physical behavior models should be considered. This paper investigates springback behavior predicted by advanced elastoplastic hardening models which combine isotropic and kinematic hardening and take strain-path changes into account. A dislocation-based microstructural hardening model formulated from physical observations and the more classical cyclic model of Chaboche have been considered in this work. Numerical implementation was carried out in the ABAQUS software using a return mapping algorithm with a combined backward Euler and semi-analytical integration scheme of the constitutive equations. The capability of each model to reproduce transient hardening phenomena at abrupt strain-path changes has been shown via simulations of sequential rheological tests. A springback analysis of strip drawing tests was performed in order to emphasize the impact of several influential parameters, namely: process, numerical and behavior parameters. The effect of the two hardening models with respect to the process parameters has been specifically highlighted.     

Constitutive modeling of soil-structure interface through the concept of critical state soil mechanics

The behavior of soil-structure interface can be crucial to the overall response of a soil-structure system. The numerical simulation of soil-structure interaction problem requires proper modeling of the interface. The similarity between the behavior of soil and interface is first analyzed in the present paper. With this similarity, the concept of critical state soil mechanics (CSSM), which has been successfully used in the modeling of soil behavior, is used to develop a constitutive model for soil-structure interface in the framework of generalized plasticity. The model is capable of modeling strain hardening, softening, normal dilatancy and stress-path dependency of interface between sandy soil and structures during shearing. The effects of normal pressure as well as density of sand are captured in the model. The performance of the model is verified with various experimental results. The unified modeling of the behavior of interfaces with different roughness, different density of soil and different normal pressures using the concept of CSSM is also successfully attempted.     

一种黏弹塑性统一本构模型

经过对大量有关统一本构模型的文献资料分析，指出了现有统一本构模型存在的问题，并通过对材料实验数据的分析，指出了黏弹性现象在实验中的表现，并据此将黏弹性引入到弹性黏塑性统一本构模型之中，建立了黏弹塑性统一本构模型，通过模型的数值模拟证明：模型计算结果无论在变形趋势上，还是在数值精度上都与实验数据符合得很好，克服了此前统一本构模型存在的问题。黏弹塑性统一本构模型的产生将统一本构模型的产生将统一本构理论的内涵扩大到黏弹性范围，进而构造了一个黏弹塑性理论的新框架。     

含离子类溶质黏土的 C-H-M 本构模型简

在水力-力学本构模型基础上,采用化学软化公式引入孔隙水中含盐、碱类等离子类溶质对土体前期固结压力的影响,得到土体的化学-水力-力学耦合本构模型（C-H-M 模型）. 根据试验结果对化学软化公式的参数进行数值拟合,通过算例验证了C-H-M 本构模型的有效性     

On shakedown analysis in hardening plasticity

The extension of classical shakedown theorems for hardening plasticity is interesting from both theoretical and practical aspects of the theory of plasticity. This problem has been much discussed in the literature. In particular, the model of generalized standard materials gives a convenient framework to derive appropriate results for common models of plasticity with strain-hardening. This paper gives a comprehensive presentation of the subject, in particular, on general results which can be obtained in this framework. The extension of the static shakedown theorem to hardening plasticity is presented at first. It leads by min-max duality to the definition of dual static and kinematic safety coefficients in hardening plasticity. Dual static and kinematic approaches are discussed for common models of isotropic hardening of limited or unlimited kinematic hardening. The kinematic approach also suggests for these models the introduction of a relaxed kinematic coefficient following a method due to Koiter. Some models for soils such as the Cam-clay model are discussed in the same spirit for applications in geomechanics. In particular, new appropriate results concerning the variational expressions of the dual kinematic coefficients are obtained.     

土水特征曲线的通用数学模型研究

土水特征曲线对于研究非饱和土的物理力学特性至关重要。根据土水特征曲线可以确定非饱和土的强度、体应变和渗透系数 ,甚至可以确定地下水面以上水份分布。由于土体物理力学特性的差异 ,导致描述其土水特征曲线的数学模型也各不相同。因此 ,建立土水特征曲线的通用数学表达式 ,显得尤为必要。本文对土水特征曲线数学模型进行了研究 ,依据这些模型的数学表达式形式 ,将其划分为 4种类型。分别由这 4种类型的数学模型推导出具有统一表达式形式的土水特征曲线通用数学模型 ,并运用陕北高原黄土土水特征曲线试验数据对通用数学模型进行了研究     

Time-dependent ratchetting experiments of SS304 stainless steel

The time-dependent strain cyclic characteristics and ratchetting behaviours of SS304 stainless steel were investigated by uniaxial/multiaxial cyclic loading tests at room and elevated temperatures (350 and 700 °C). The effects of loading rate, peak/valley strain or stress holds, ambient temperature and non-proportional loading path on the cyclic softening/hardening and ratchetting behaviours of the material were discussed. It is shown that: the cyclic deformation of the material presents remarkable time-dependence at room temperature and 700 °C; the cyclic hardening feature and ratchetting strain depend significantly on straining or stressing rate, hold-time, ambient temperature and the non-proportionality of loading path; the time-dependent ratchetting is resulted from the slight opening of hysteresis loop and visco-plasticity together, and the viscosity is a dominating factor at 700 °C; at 350 °C, abnormal rate-dependence and quick shakedown of ratchetting are observed due to the dynamic strain aging of the material at this temperature. Some significant conclusions are obtained, which are useful to construct a constitutive model to describe the time-dependent ratchetting behaviour of the material. It is also stated that the unified visco-plastic constitutive model discussed here cannot provide reasonable simulation to the time-dependent ratchetting at 700 °C, especially to that with certain peak/valley stress hold, since the effect of the high viscosity on time-dependent ratchetting cannot be properly described by using a unified visco-plastic flow rule.     

考虑塑性应变记忆恢复的循环棘轮本构模型研究

在统一粘塑性循环本构理论框架下,以Ohno-Abdel-Karim非线性随动硬化模型为基础,建立了一个循环本构模型。模型通过引入塑性应变幅值记忆效应,并在塑性应变记忆项中加入恢复系数,提高了对循环硬化材料单轴棘轮行为的预言能力。将模型应用于316L不锈钢单轴棘轮行为的描述中,模拟不同平均应力、应力幅值下的棘轮应变,均与实验数据吻合较好,证明本文改进的本构模型能合理地描述循环硬化材料的单轴棘轮行为。     

Plastic collapse in non-associated hardening materials with application to Cam-clay

This paper presents a variational formulation for the analysis of plastic collapse conditions for a class of hardening materials that accounts for some non-associated flow laws such as the modified Cam-clay model of soils. In this framework, classical statical and kinematical principles of limit analysis do not hold. The variational principle is formulated for the general class of materials whose flow equations are derived from a kind of generalized potentials named bipotentials by de Saxcé.The plastic collapse phenomenon for hardening materials is considered first and formulated as a system of equations. In particular, the case of the usual modified Cam-clay model is analyzed. The paper follows with the proposal of a minimization principle whose solution is then related to the solution of the plastic collapse problem. We demonstrate the use of this minimum principle in a simple example of triaxial compression of a modified Cam-clay material. Finally, we discuss the particular form of the proposed variational formulation for the case of associated plasticity.     

冲击载荷下盐岩的力学性能和本构关系研究

利用直径75mm的大尺寸SHPB实验装置开展盐岩的动态压缩性能实验研究,得到了其动态应力应变曲线。分析表明,在冲击压缩载荷作用下,盐岩材料有明显的损伤软化现象和应变率效应。针对实验曲线,通过引入描述材料强度随应变率强化的应变率增强因子和随不可逆变形发展而弱化的损伤因子,提出了含损伤率的相关动态本构模型,拟合得到的本构方程形式比较简单,能够较好地反映盐岩动态实验加载过程的主要特征,具有一定的工程应用价值。