非饱和土力学理论的研究进展

回顾了非饱和土有效应力的发展,目前普遍认同采用两个应力变量来建立本构模型,且对基质吸力中毛细和粘吸两部分作用进行了阐述。分析了非饱和土强度问题,包括抗剪强度和抗拉强度。讨论了非饱和土的本构模型问题,包括基于净应力和基质吸力的弹塑性模型,基于Bishop有效应力和基质吸力的水力力学耦合弹塑性模型,以及双孔隙结构的模型。最后探讨了热力学方法和多孔介质理论在非饱和土中的应用,基于多孔介质理论在多场耦合条件下土体复杂的行为是当前值得研究的问题。     

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

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

Randomness and slip avalanches in gradient plasticity

While localization of deformation at macroscopic scales has been documented and carefully characterized long ago, it is only recently that systematic experimental investigations have demonstrated that plastic flow of crystalline solids on mesoscopic scales proceeds in a strongly heterogenous and intermittent manner. In fact, deformation is characterized by intermittent bursts (‘slip avalanches’) the sizes of which obey power-law statistics. In the spatial domain, these avalanches produce characteristic deformation patterns in the form of slip lines and slip bands. Unlike to the case of macroscopic localization where gradient plasticity can capture the width and spacing of shear bands in the softening regime of the stress–strain graph, this type of mesoscopically jerky like localized plastic flow is observed in spite of a globally convex stress–strain relationship and may not be captured by standard deterministic continuum modelling. We thus propose a generalized constitutive model which includes both second-order strain gradients and randomness in the local stress–strain relationship. These features are related to the internal stresses which govern dislocation motion on microscopic scales. It is shown that the model can successfully describe experimental observations on slip avalanches as well as the associated surface morphology characteristics.     

砂土的应力路径本构模型

将微元应力路径线性逼近,转变成与其充分接近且易于计算应变的等平均应力微元和等应力比微元,计算任意加荷应力路径所产生的塑性应变,建立了双屈服面的砂土应力路径本构模型．模型体现了岩土塑性理论分量屈服和非关联流动法则的要求,在p,q平面内根据双线性的屈服线确定了加卸载准则．结合广义非线性强度理论采用变换应力三维化方法简单、合理地使模型实现三维化．通过试验数据的验证表明,砂土应力路径本构模型可以合理地描述各种应力路径下砂土的变形和强度特性。     

Multi-scale analysis of subgrid stress and energy dissipation in turbulent channel flow

In present study, the subgrid scale （SGS） stress and dissipation for multiscale formulation of large eddy simulation are analyzed using the data of turbulent channel flow at Ret = 180 obtained by direct numerical simulation. It is found that the small scale SGS stress is much smaller than the large scale SGS stress for all the stress components. The dominant contributor to large scale SGS stress is the cross stress between small scale and subgrid scale motions, while the cross stress between large scale and subgrid scale motions make major contributions to small scale SGS stress. The energy transfer from resolved large scales to subgrid scales is mainly caused by SGS Reynolds stress, while that between resolved small scales and subgrid scales are mainly due to the cross stress. The multiscale formulation of SGS models are evaluated a priori, and it is found that the small- small model is superior to other variants in terms of SGS dissipation.     

Effective Stress in Unsaturated Soils: A Thermodynamic Approach Based on the Interfacial Energy and Hydromechanical Coupling

In recent years, the effective stress approach has received much attention in the constitutive modeling of unsaturated soils. In this approach, the effective stress parameter is very important. This parameter needs a correct definition and has to be determined properly. In this paper, a thermodynamic approach is used to develop a physically-based formula for the effective stress tensor in unsaturated soils. This approach accounts for the hydro-mechanical coupling, which is quite important when dealing with hydraulic hysteresis in unsaturated soils. The resulting formula takes into account the role of interfacial energy and the contribution of air?Cwater specific interfacial area to the effective stress tensor. Moreover, a bi-quadratic surface is proposed to represent the contribution of the so-called suction stress in the effective stress tensor. It is shown that the proposed relationship for suction stress is in agreement with available experimental data in the full hydraulic cycle (drying, scanning, and wetting).     

饱和土中剪切带宽度的研究

主要对饱和土在简单剪切和耦合载荷下的剪切带宽进行了理论分析. 在模型中考虑了剪胀、孔压耗散、惯性及孔隙水和土颗粒之间阻力等因素的耦合效应,将有效应力考虑为应变、应变率和孔隙压力的函数. 通过对模型的简化分析得到简单剪切和耦合应力状态下的剪切带宽度的近似计算公式. 分析表明,当土体失稳后,会出现有限厚度的剪切带, 且其形成是一个后期过程,可以忽略惯性效应. 剪切带的厚度依赖于孔隙压力特性和土体骨架的剪胀特性,但是不依赖应力状态. 最后将分析与实验结果进行了对比,在一定程度上,两者是符合的.     

饱和土中球空腔的瞬态动力响应

采用工程上通用的饱和土力学模型,考虑流体与固全之间的耦合作用,利用Ｌａｐｌａｃｅ变换求解了饱和土中球空腔的瞬态动力响应问题,得到了变换域内的解析解,借助数值Ｌａｐｌａｃｅ换求解了饱和土中球空腔瞬态动力响应的位移、应力及孔压的变化规律,为分析地下结构瞬态动力响应提供了一种有效的方法,模型符合工程实际,有一定的工程应用价值。     

Saturated and unsaturated behaviour modelling of Meuse–Haute/Marne argillite

Poromechanical behaviour modelling of the Callovo–Oxfordian argillite under saturated and partially saturated conditions is proposed using the equivalent stress concept. In comparison with the previous works on this rock, the particular form of the yield criterion and the plastic flow potential proposed here help to better describe the rock behaviour in tension–stress paths. The evolution of the poroelastic parameters due to the induced cracks is also considered in a simple way. Due to its physical nature, different from classical soils, the partially saturated behaviour of this rock could not be correctly described by any of the pre-existing theories usually used for the partially saturated porous media (soils). Based upon experimental results on this rock and developing some ideas proposed by other authors, an extension of the saturated elastic–plastic model in unsaturated field is proposed. The key hypothesis of this extension is the evolution of the Biot’s coefficient as a function of the suction, justified by laboratory results and micromechanical analyses. The predictions of the model in saturated and partially saturated conditions are compared with laboratory results and a good general agreement is found.     

等p路径十砂土本构关系的归一化特性及数值建模方法

用神经网络对土体进行建模能反映应力路径相关性、反映土的剪胀剪缩以及反映体应力、剪应力对体应变、剪应变的交互影响,因而成为一种比较理想的建模方式.能否在样本有限的情况下获得精度比较高的本构模型正是主要的研究目的.通过研究中密砂在等p路径下的三轴试验曲线,发现其应力-应变关系曲线在常规应力范围内具有归一化特性.选择合适的归一化指标对砂土三轴试验数据进行归一化,以归一化的试验数据为训练样本进行神经网络训练,得到了比较理想的砂土的神经网络本构模型.本构模型仿真值与试验值符合较好,表明所给出的建模方法是合理的.提出的建模方法可以在所有试验数据的基础上自动实现概率寻优,能有效降低噪声信号的干扰、减小试验数据的分散造成的影响.     

A Gas-Kinetic Scheme for Turbulent Flow

RANS simulations may not provide accurate results for all flow conditions. The interaction between a shock wave and a turbulent boundary layer is an example which may still be difficult to simulate accurately. Beside the inability to reproduce physical phenomena such as shock unsteadiness, the argument is put forward that the conventional numerical schemes, based on the Navier-Stokes equations, may be unable to generate a physically consistent turbulent stress tensor in the presence of large unresolved scales of motion. A large ratio between unresolved and resolved scales of motion, a sort of Knudsen number based on turbulent fluctuations, might introduce inaccuracies for which the turbulence model is not accountable. In order to improve the accuracy of RANS simulations, researchers have suggested various ad-hoc modifications to standard turbulence models which limit eddy viscosity or the turbulent stress tensor in the presence of strong gradients. Gas-kinetic schemes might be able to improve RANS predictions in shocklayers by removing or limiting the errors caused by the large scales ratio. These schemes are a class of their own; in the framework of a finite-volume or finite-elements discretizations, they model the numerical fluxes on the basis of the Boltzmann equation instead of the Navier-Stokes equations as is conventionally done. In practical terms, these schemes provide a higher accuracy and, more importantly, an in-built “multiscalar” mechanism, i.e. the ability to adjust to the size of unresolved scales of motion. This property makes them suitable for shock-capturing and rarefied flow. Gas-kinetic scheme may be coupled to a conventional RANS turbulence model; it is shown that the turbulent stress tensor is naturally adjusted as a function of the unresolved-to-resolved scales ratios and achieves a higher physical consistency than conventional schemes. The simulations shown - well-known benchmark cases with strong shock-boundary layer interactions - have been obtained with a standard two-equation turbulence model (k- ω). It is shown that the gas-kinetic scheme provides good quality predictions, where conventional schemes with the same turbulence model are known to fail.     

速度估计模型在大涡模拟中的应用

介绍了速度估计模型的基本思想及其在物理空间的实现。速度估计模型依靠湍流大尺度的非线性作用估计出小尺度,从而可直接求解亚格子应力项而不需要额外的模型。本文采用该模型对不同雷诺数下的各向同性衰减湍流进行了模拟,并与直接数值模拟、理论分析和其他亚格子模型的大涡模拟结果进行了比较。也初步考察了网格分辨率、不同精度的紧致格式对大涡模拟结果的影响。     

Mode I and mixed mode crack-tip fields in strain gradient plasticity

Strain gradients develop near the crack-tip of Mode I or mixed mode cracks. A finite strain version of the phenomenological strain gradient plasticity theory of Fleck-Hutchinson (2001) is used here to quantify the effect of the material length scales on the crack-tip stress field for a sharp stationary crack under Mode I and mixed mode loading. It is found that for material length scales much smaller than the scale of the deformation gradients, the predictions converge to conventional elastic-plastic solutions. For length scales sufficiently large, the predictions converge to elastic solutions. Thus, the range of length scales over which a strain gradient plasticity model is necessary is identified. The role of each of the three material length scales, incorporated in the multiple length scale theory, in altering the near-tip stress field is systematically studied in order to quantify their effect.     

单向压缩状态下上海地区软土的蠕变变形与次固结特 性研究

基于上海地区原状软土的一维固结试验,讨论了在不同加荷比、加荷速率及超载预压情况下软土的蠕变变形特性和次固结特性。试验结果表明:蠕变变形随应力水平、加荷速率的增大而增大。相同应力水平下,加荷比较大的土样蠕变变形较小。蠕变变形为非线性,随着应力水平的增加,蠕变变形从粘弹性到粘塑性转变,加荷比较大时,这种转变更明显。对试验结果进行曲线拟和,确定Metchell一维经验模型参数,拟和曲线与试验结果吻合较好。通过对次固结阶段试验数据的线性拟合可知,对于正常固结土,次固结系数Ca值受应力水平和加荷比的影响不大,而与加荷速率密切相关,加荷速率较大时,Ca值随荷载的增大急剧增大后略有降低。超载预压可有效减小次固结系数,从而有利于减小工后次固结沉降。     

MICROMECHANICS ANALYSIS FOR UNSATURATED GRANULAR SOILS

This paper aims at establishing an anisotropic stress expression for unsaturated pendular-state granular soils. Using the second-order fabric tensor, we formulate a micromechanics scheme of soils with statistically averaging method, and reveal that the macroscopic average stress of unsaturated granular soils in pendular-state is not isotropic. Not only is the stress from contact forces anisotropic due to the fabric, but also the capillary stress is directional dependent, which is different from the common point that the capillary stress is isotropic. The capillary stress of unsaturated pendular-state granular soils is determined by the orientation distribution of contact normals, so it is closely related to the initial and induced anisotropy of soils. Finally, DEM numerical simulations of triaxial compression tests of pendular-state soils at different degrees of saturation are used to verify the existence of above anisotropy of stresses.     

CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$

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

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

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

土-地铁隧道动力相互作用模型试验分析

本文对土的动力特性本构模型以及土-结构的接触面模型进行了较详细的探讨,最后分别选用非线性等价粘弹性模型和Ramberg-Osgood模型。分析表明,在进行模型试验设计时,不仅材料要满足相似性原理,而且在设计边界条件和地震波输入时也应以上述模型为参照。通过对一圆形隧道的地震反应分析,表明拱顶处的最大动应力是圆形隧道出现裂痕的主要原因,进而分析了孔隙率对结构横向变形的影响,分析结果与实际工程情况较为吻合。本文的研究成果可为地铁隧道的抗震设计提供一定的参考意见。     

One-Dimensional Consolidation of Layered and Visco-Elastic Soils under Arbitrary Loading

ＩｎｔｒｏｄｕｃｔｉｏｎＡｆｔｅｒＴｅｒｚａｇｈｉ,ａｌｏｔｏｆｓｃｈｏｌａｒｓｈａｖｅｄｏｎｅｏｕｔｓｔａｎｄｉｎｇｒｅｓｅａｒｃｈｏｎｔｈｅｐｒｏｂｌｅｍｓｏｆｏｎｅ＿ｄｉｍｅｎｓｉｏｎａｌｃｏｎｓｏｌｉｄａｔｉｏｎ .Ｗｉｌｓｏｎ ( 1974) [1] ａｎｄＢａｌｉｇｈ( 1978) [2 ] ｅｔａｌ.ｈａｖｅｓｔｕｄｉｅｄｔｈｅｃｏｎｓｏｌｉｄａｔｉｏｎｕｎｄｅｒｔｈｅｒｅｃｔａｎｇｕｌａｒｌｏａｄｉｎｇ .Ｉｎ 1988,ＷＵｅｔａｌ.[3] ｗｏｒｋｅｄｏｕｔｔｈｅｇｅｎｅｒａｌｓｏｌｕｔｉｏｎｓｏｆｏｎｅ＿ｄｉｍｅｎ…     

ADVANCES IN TRUE TRIAXIAL TEST OF UNSATURATED SOILS

研究非饱和土的强度是非饱和土理论及其工程应用的首要问题,但因复杂应力状态下基质吸力的控制技术、量测技术和吸力平衡时间等原因,非饱和土真三轴试验研究的进展迟缓. 在介绍非饱和土刚性真三轴仪、柔性真三轴仪和刚柔复合型真三轴仪的基础上,重点分析已有非饱和土刚性和柔性真三轴试验结果,总结现有非饱和土真三轴仪及其试验研究的不足,同时指出Mohr-Coulomb 强度准则和外接圆Drucker-Prager准则对非饱和土真三轴试验结果的不适用性. 应开展更多不同种类非饱和土的完整真三轴试验研究,特别是非饱和黏性土;应结合非饱和土的应力状态变量和强度特性,建立符合工程实际受力状况的非饱和土真三轴强度准则,完善非饱和土的理论基础.