基于CLoE亚塑性模型的密砂循环荷载下力学行为模拟

亚塑性模型为模拟颗粒材料的非线性力学行为提供了一条新途径,特别是CLoE亚塑性模型在模拟应变局部化时具有一定的优势。然而该模型在模拟小幅应力-应变循环时表现出一定的锯齿效应。为了克服该效应,基于颗粒间应变张量的概念发展了修正的CLoE亚塑性模型以正确模拟循环荷载下密砂的力学行为。此外,为保证单调荷载作用下修正模型与原模型预测结果的一致性,改进了颗粒间应变率及颗粒间最大应变的定义。数值算例表明:（1）修正模型保留了克服锯齿效应的优点。（2）修正模型能够反映不同振幅条件下的卸载刚度。（3）在大振幅循环条件下,滞回圈的面积随着循环次数增加而增大。（4）修正模型能够保证单调加载条件下所得结果与原模型的一致性。（5）修正模型可以反映材料的疲劳破坏机制。     

Gudehus-Bauer亚塑性本构模型研究

无黏性土的应力-应变关系可以用Gudehus-Bauer亚塑性本构模型来模拟,该模型强调应力增量的大小和方向不仅与当前应力状态有关,而且还取决于当前应变增量的大小和方向。为分析其与传统弹塑性理论的不同之处,对Gudehus-Bauer理论的线性项和非线性项进行了研究,并对不同初始孔隙比下Gudehus-Bauer亚塑性模型的应力-应变关系进行了探讨。结果表明Gudehus-Bauer亚塑性模型不用把应变分为弹性和塑性部分就能考虑不可逆变形,并能体现密砂的剪胀特性和应变软化特性以及松砂的剪缩特性和应变硬化特性。     

堆石料Gudehus-Bauer亚塑性本构模型改进及参数确定方法

研究了Gudehus-Bauer亚塑性本构模型和模型参数的求取方法。采用侧限压缩试验曲线求取模型参数颗粒硬度hs和指数n。根据模型方程的推导,建立了拟合指数 和 与围压之间的关系,并提出了新的拟合参数。考虑到堆石料具有明显的剪胀、剪缩性,在Gudehus-Bauer模型线性项中增加了主要控制体积应变项 ,以改善模型对堆石料体积应变曲线的描述。采用堆石料大型侧限压缩试验、常规三轴试验分别验证了新的拟合参数和改进后的Gudehus-Bauer亚塑性本构模型。与堆石料试验成果比较,提出的新拟合参数与改进后的Gudehus-Bauer亚塑性本构模型可以较好地模拟堆石料的应力-应变特性,并较好地改善了堆石料体变曲线的模拟结果。对改进后的模型作了常规三轴加、卸载模拟,模拟结果反映了改进的Gudehus-Bauer亚塑性本构模型具有一定的卸载适应性。     

一类新型的散粒型土体本构理论

详细地介绍了国际上近30才发展起来的具有独特建模思想的一类新型的散粒型土体本构建模理论--亚塑性本构理论。将该理论与经典的弹塑性理论进行比较后发现,亚塑性理论的本构表达式自动隐含了弹塑性理论中相对应的一些基本概念和假定,无需额外引入,因此增加了模型的客观性。以Gudehus-Bauer亚塑性模型为例,分析了模型在临界状态时的本构特征。此外,通过三轴试验和循环剪切试验的亚塑性数值模拟表明,亚塑性本构模型能很好地反映无黏性散粒型土体的非弹性、非线性及剪胀（剪缩）性等主要应力-应变特性。     

可液化倾斜场地的侧向扩展机制分析

针对已完成的液化水平自由场大型振动台试验,采用OpenSees有限元平台,建立了振动台试验的数值分析模型,并验证了该模型的可靠性。基于此,建立了地基整体倾斜的自由场数值模型,重点讨论了场地土体的非循环地震响应和液化侧向扩展机制。结果表明,建立的数值模型可以有效地模拟可液化场地的地震反应。倾斜场地中,可液化松砂与上覆非液化层界面处具有显著的相对位移,饱和砂土的应变累积从松砂层浅层开始,逐步向深层发展,超孔隙水压力增长和土体非循环应变累积未表现出耦合的现象,场地的中部土体控制非循环横向位移的发展。另一方面,土体液化过程中,当沿滑动面的剪切应力小于初始静剪应力时,液化侧向扩展启动。此时饱和松砂层的剪应力比在0.04～0.06范围内,略小于初始静剪应力比。此外,还发现液化诱导侧向扩展需要场地具有一定的倾斜度（大于0.5°）;土的侧向位移符合余弦分布模式;随着场地倾斜度的增大,可液化深层土对整体侧向位移的贡献更显著。     

强度参数对初始随机缺陷岩样全部变形的影响

在单轴平面应变压缩条件下, 采用FLAC模拟了初始内聚力及内摩擦角对具有随机材料缺陷岩样轴向、侧向、体积变形及由侧向应变及轴向应变计算得到的泊松比的演变的影响。采用编写的FISH函数于试样内部规定随机缺陷并计算其全部变形特征。密实的岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则, 破坏之后呈现应变软化—理想塑性行为。随着强度参数(初始内聚力及内摩擦角)的增加, 应力峰值及对应的侧向应变的值提高, 体积应变的峰值增加, 应力峰值所对应的计算得到的泊松比稍有增加。由于被剪切带所分割的毗邻块体之间的相对滑动, 纵然扩容角为零, 试样在峰后的变形阶段, 仍然可以观察到体积膨胀现象, 这与作者的理论分析结果一致。当强度参数降低时, 侧向应变—轴向应变曲线、体积应变—轴向应变曲线及计算得到的泊松比—轴向应变曲线的峰前非线性部分变得短暂。剪切带倾角的数值解的上限低于Coulomb理论, 下限在Roscoe理论附近波动, Arthur理论的预测结果与本文剪切带倾角的数值解更接近.      

考虑颗粒破碎对特征孔隙比影响的堆石体亚塑性本构模型

颗粒破碎是影响堆石体强度和变形特性的主要问题之一。相比于砂土,堆石料在较低的应力水平下就会发生严重的颗粒破碎,因此,在进行堆石体力学特性及本构模型研究时必须考虑颗粒破碎的影响。同时,堆石体在受力过程中孔隙比是变化的,而传统本构模型不能使用一组参数模拟不同孔隙比的同种材料。因此,以能够考虑应力水平和土体孔隙比影响的Gudehus-Bauer亚塑性本构模型为基础,考虑堆石体有别于砂土的孔隙比变化特征,提出了考虑堆石破碎的亚塑性本构模型。亚塑性理论是目前可最大限度地减少人为假定的一种本构理论,颗粒材料在不同特征应力路径下,破碎造成的过度变形量不同;但相同应力水平、不同特征应力路径下孔隙比已不满足Gudehus-Bauer亚塑性本构模型中提出的等比例变化规律。据此,结合考虑颗粒破碎的临界状态理论和堆石体常规三轴试验和循环加载试验结果,提出了考虑颗粒破碎堆石体特征孔隙比的表达式,并将其引入到Gudehus-Bauer亚塑性本构模型中,建立了考虑颗粒破碎的堆石体亚塑性本构模型,提出了模型参数的确定方法。与堆石体试验结果对比表明,该本构模型可以较好地模拟其力学与变形特性。     

基于Cosserat理论的重力坝深层抗滑稳定分析

重力坝深层滑动主要表现为沿缓倾角的软弱结构面形成滑移通道,滑移通道内应变积聚且应变梯度急剧不连续,是典型的应变局部化现象,采用经典连续介质理论进行数值模拟时存在病态的有限元网格依赖性。引入Cosserat连续体理论作为正则化机制,提出了基于Cosserat理论的Mohr-Coulomb弹塑性模型,考虑非关联的流动法则,在经典塑性理论框架下采用向后Euler隐式积分算法进行应力更新。采用ABAQUS的自定义单元接口（UEL）进行二次开发,进行了平面应变条件下单轴受压的数值验证。数值模拟结果表明,该模型能保证应变局部化问题的正定性。基于Cosserat理论的重力坝深层抗滑稳定分析结果表明,采用超载法进行重力坝渐进破坏过程模拟时,基于经典连续体理论的模拟结果有较大的网格依赖性,而且结果偏于安全,而采用Cosserat连续体理论的结果对网格密度不敏感。     

岩石试件端面摩擦效应数值模拟研究

试件端面摩擦效应直接影响试件内的塑性等效应变、侧向位移的分布和单元应力应变曲线。本文运用ANSYS中的接触单元模拟了平面应变状态下端面摩擦效应对塑性等效应变、侧向位移和单元应力应变曲线的影响,得到了不同摩擦系数时塑性等效应变及侧向位移的渐进变化形式。当接触面摩擦较小时,塑性等效应变图案为上下两个X形网络,侧向位移上下分布均匀;当接触面摩擦增大时,塑性等效应变网络向中部靠拢并且明显增大,侧向位移上下分布不均匀,中部较上下端面位移大;当试件端面侧向位移被限制,即摩擦力很大时,塑性等效应变网络变为一个X形局部化带,侧向位移分布更加不均匀,中部明显隆起。     

考虑湿化效应的堆石料Gudehus-Bauer 亚塑性模型应用研究

在扼要分析堆石料主要力学特性的基础上,在亚塑性理论框架内,采用扩展的Gudehus-Bauer亚塑性模型对堆石料因含水率发生变化引起的湿化行为进行建模预测。扩展模型除能考虑含水率对颗粒硬度的影响外,还能反映土体密实度、应力状态和含水率等对堆石料非线性和非弹性变形的影响。将该模型应用于心墙坝初次蓄水时应力变形的亚塑性计算分析中。结果表明,扩展模型能较合理地反映心墙堆石坝的湿化变形规律。     

Mohr–Coulomb MiniCLoE model Uniqueness and localization studies,links with normality rule

This paper is devoted to a parametric study of a plane Mohr–Coulomb CLoE model. As CLoE models are designed with a consistency condition, it is possible to define a normality condition and to study its consequences. The positiveness of the second order work which implies the uniqueness of the solution of a small strain boundary value problem is studied firstly. Then the localization criterion is also studied. It is proved that normality has consequences similar to those for classical elasto plastic models. However if induced anisotropy is introduced in the hypoplastic CLoE model, some conclusions are no longer true. Finally plane strain experimental data are used to identify the parameters of the model. Copyright © 2002 John Wiley & Sons, Ltd.     

A hypoplastic model for granular soils incorporating anisotropic critical state theory

It is well known that soil is inherently anisotropic and its mechanical behavior is significantly influenced by its fabric anisotropy. Hypoplasticity is increasingly being accepted in the constitutive modeling for soils, in which many salient features, such as nonlinear stress-strain relations, dilatancy, and critical state failure, can be described by a single tensorial equation. However, within the framework of hypoplasticity, modeling fabric anisotropy remains challenging, as the fabric and its evolution are often vaguely assumed without a sound basis. This paper presents a hypoplastic constitutive model for granular soils based on the newly developed anisotropic critical state theory, in which the conditions of fabric anisotropy are concurrently satisfied along with the traditional conditions at the critical state. A deviatoric fabric tensor is introduced into the Gudehus-Bauer hypoplastic model, and a scalar-valued anisotropic state variable signifying the interplay between the fabric and the stress state is used to characterize its impact on the dilatancy and strength of the soils. In addition, fabric evolution during shearing can explicitly be addressed. Modifications have also been undertaken to improve the performance of the undrained response of the model. The anisotropic hypoplastic model can simulate experimental tests for sand under various combinations of principle stress direction, intermediate principal stress (or mode of shearing), soil densities, and confining pressures, and the associated drastic effect of different principal stress orientations in reference to the material axes of anisotropy can be well captured.     

适用于砂土循环加载分析的边界面塑性模型

基于临界状态土力学框架,建立了一个适用于砂土排水循环加载的边界面塑性模型。采用了考虑虚拟峰值应力比的偏应变硬化准则,初始加载阶段应力点位于边界面上,反向加载阶段以历史最大屈服面作为边界面,同时实现了对密砂软化现象的模拟和对历史所受最大应力的记忆。边界面采用修正的椭圆形,引入考虑密度与应力水平的状态相关剪胀函数,采用非相关联流动法则和以应力反向点作为映射中心的径向映射准则。模型仅有10个参数,通过常规三轴试验即可确定,并且使用一套参数可以模拟不同围压、密度的单调和循环加载情况。分别对饱和砂土的单调、循环排水三轴试验进行模拟,结果表明,该模型能够合理地反映饱和砂土排水条件下的应力-应变特性。     

Directional response of a reconstituted fine‐grained soil—Part II: performance of different constitutive models

In this paper, the performance of different advanced constitutive models for soils is evaluated with respect to the experimentally observed behaviour of a soft reconstituted clay subject to a wide range of loading directions, see (presented in the companion paper). The models considered include a three‐surface kinematic hardening elastoplastic model; the CLoE hypoplastic model; a recently proposed K‐hypoplastic model for clays, and an enhanced version of the same model incorporating the concept of intergranular strain. A clear qualitative picture of the relative performance of the different models as a function of the loading direction is obtained by means of the incremental strain response envelopes. The definition of suitable error measures allows to obtain further quantitative information in this respect. For the particular initial conditions and loading programme considered in this study, the kinematic hardening and the enhanced K‐hypoplastic models appear to provide the best performance overall. Copyright © 2006 John Wiley & Sons, Ltd.     

Performance of mixed and enhanced finite elements for strain localization in hypoplasticity

Displacement and mixed finite element formulations of shear localization in materials are presented. The formulations are based on hypoplastic constitutive laws for soils and the mixed enhanced treatment involving displacement, strain and stress rates as independently varied fields. Included in these formulations are the standard displacement method, the three‐field mixed formulation, the enhanced assumed strain method and the mixed enhanced strain method. Several numerical examples demonstrating the capability and performance of the different finite element formulations are presented. The numerical results are compared with available experimental data for Hostun RF sand and numerical results for Karlsruhe sand on biaxial tests. Copyright © 2011 John Wiley & Sons, Ltd.     

Elastoplastic model with loading memory surfaces (LMS) for monotonic and cyclic behaviour of geomaterials

A new elastoplastic model called loading memory surface based on the critical state concept and the multi‐surface framework is proposed for geomaterials. The model uses a hypoelastic formulation and two plastic mechanisms. The formulations of the model are made in three‐dimensional stress–strain space and work under both monotonic and cyclic loadings. A newly introduced formalism makes it possible to obtain the cyclic response directly from the monotonic loading one. This formalism gives a three‐dimensional generalization of the well‐known Masing rule. The model has been validated against test results of Hostun sand under several conditions: monotonic and cyclic, drained and undrained, tests in compression and in extension, and at different confining pressures and different densities. Copyright © 2014 John Wiley & Sons, Ltd.     

Hypoplastic model for sands with loading surface

Although the hypoplastic models for sands have exhibited good predictive capability in monotonic loading, they are not able to reproduce memory effects and predict excessive plastic accumulation under cyclic loading. To overcome these issues, a loading surface has been incorporated into a hypoplastic model. This surface is capped and has two hardening variables. Notions from the bounding surface plasticity were borrowed in order to formulate the hardening functions. With this novel model, some salient features can be described: the model can account for the accumulation of plastic deformation, a memory effect is provided by the new surface, and stress-induced anisotropy effects observed in sands are successfully simulated. A short calibration guide of the parameters is given, and some simulations for Hostun RF loose sand and Toyoura sand are presented.