复杂加载条件下的砂土本构模型

试验表明,饱和砂土的应力应变关系具有显著的密度以及压力依存性,上述两点构成了描述砂土静力加载下变形特性无法忽视的因素.此外,在循环加载等复杂加载作用下,砂土还会表现出明显的应力诱导各向异性以及相变转换特性.基于在e-p空间中存在唯一的临界状态线这一基本假定,通过在e-p空间中引入当前状态点与临界状态线的距离R来作为反映密度与压力依存特性的状态参量,将变相应力比以及峰值应力比表达为状态参量的指数函数,将上述应力比参量引入到统一硬化参量中可准确地反映初始状态下围压、密度对于单调加载下应力应变关系的影响规律,能描述砂土剪缩、剪胀,应变软化、硬化等特性.采用非相关联流动法则,p-q空间中采用水滴型屈服面,塑性势面为椭圆面,松砂在单调加载下的静态液化现象也可描述.为反映循环加载下塑性体积应变的累积特性以及塑形偏应变的滞回特性,在循环加载下将状态参量R表达为应力比参量,并在硬化参数中引入描述应力诱导各向异性特性的旋转硬化部分,所提模型可有效地描述循环加载下剪切模量的衰减特性、刚度衰化性质、强度减小特性,在不排水约束作用下,则会产生往返活动性现象.通过一系列的模型模拟与试验结果对比,验证了本构模型的有效性及适用性.     

含水合物沉积物的弹塑性本构模型

土的密度对其力学特性具有明显影响.水合物以一种固相赋存于沉积物的孔隙中,使得水合物的含量和其赋存形式都会影响含水合物沉积物(GHBS)的密度,因此在研究和描述含水合物沉积物的力学性质时应考虑水合物含量和赋存形式对其密度的影响.本文基于黏土和砂土统一的本构模型(CSUH模型),首先建立水合物体积分数与压硬性参量的关系式来反映水合物对沉积物压缩规律的影响.其次,为了合理考虑水合物含量和赋存形式对沉积物密度的影响,建立了可以描述有效初始孔隙的计算式,并将其引入到状态参量中来描述水合物对沉积物剪胀性和峰值强度的影响.最后,结合CSUH模型中水滴形屈服面,建立了一个含水合物沉积物的弹塑性本构模型.通过与室内试验结果比较,验证了该模型不仅能够合理地描述不同赋存形式、不同水合物含量下含水合物沉积物的应力应变关系,而且在描述具有相同赋存形式含水合物沉积物的力学特性时,不同的水合物含量只需采用一组参数.      

恒定主应力偏转角下黏土不排水剪切特性分析

天然土的结构性、各向异性对其强度有影响.进行固结压力100kPa、主应力轴恒定偏转角为0°,15°,30°,45°,60°的恒定主应力轴偏转角的固结不排水剪切试验,分析了不同恒定主应力轴偏转角对饱和软黏土固结不排水剪切强度影响;推导了总应力控制恒定主应力轴偏转角剪切试验总应力路径,编制了总应力控制剪切三维点积分程序,对比恒定主应力轴偏转角饱和软黏土固结不排水剪切强度计算值与试验结果,表明程序能很好地反应各向异性对土剪切强度影响.     

砂土的应力空间特征面及物态模型

依据饱和砂土的剪缩剪胀性在往返荷载下的反应性状 ,揭示了饱和砂土具有初剪缩、次剪缩、剪胀、反向剪缩和弹性剪切等不同物态 ,它们在有效应力空间具有不同的路径变化。进而 ,在状态边界面和物态转换面的基础上 ,又引入了应力历史边界面和初始屈服面 ,通过有效应力状态点及应力历史边界面和初始屈服面的运动规律 ,建立了饱和砂土物态变化的数学描述 ,可区分不同的物态变化     

含水合物沉积物的弹塑性本构模型

土的密度对其力学特性具有明显影响.水合物以一种固相赋存于沉积物的孔隙中,使得水合物的含量和其赋存形式都会影响含水合物沉积物(GHBS)的密度,因此在研究和描述含水合物沉积物的力学性质时应考虑水合物含量和赋存形式对其密度的影响.本文基于黏土和砂土统一的本构模型(CSUH模型),首先建立水合物体积分数与压硬性参量的关系式来反映水合物对沉积物压缩规律的影响.其次,为了合理考虑水合物含量和赋存形式对沉积物密度的影响,建立了可以描述有效初始孔隙的计算式,并将其引入到状态参量中来描述水合物对沉积物剪胀性和峰值强度的影响.最后,结合CSUH模型中水滴形屈服面,建立了一个含水合物沉积物的弹塑性本构模型.通过与室内试验结果比较,验证了该模型不仅能够合理地描述不同赋存形式、不同水合物含量下含水合物沉积物的应力应变关系,而且在描述具有相同赋存形式含水合物沉积物的力学特性时,不同的水合物含量只需采用一组参数.     

基于临界状态的砂土本构模型研究

砂土孔隙比及所受压力是其力学特性的重要影响因素. 本文基于砂土临界状态线特性分析,采用以e-(p/p_a)ξ平面内的线性关系描述其等向压缩线. 通过对比分析两种不同压缩线函数 与临界状态线函数之间的关系提出更适合描述砂土在等向压缩下的参考压缩线,并给出了基于参考压缩线的等向硬化规律. 建议了适用于 描述砂土剪切特性的屈服面函数,并给出利用等向压缩和等p路径确定屈服面形状参数μ的方法. 将不同应力比对应的压缩线作为砂土状态参量参考线,以获取潜在强度M_f与特征状态应力比M_c,进而描述砂土压缩与剪切特性;基于等向压缩与等p路径建立了当前应力比与状态参量参考线之间的相关关系,从而实现了砂土状 态参量参考线由参考压缩线向临界状态线平稳过渡. 建立的砂土本构模型共11个参数,均能够通过常规土工试验或经验获取. 基于模型预测与Toyoura砂的等向压缩、三轴不排水剪切试验及排水剪切试验的对比结果,本文建立的砂土本构模型很好地描述了Toyoura 砂在不同孔隙比和不同压力下的压缩与剪切特性.      

超固结饱和黏性土的弹塑性本构模型及三轴试验模拟

针对饱和超固结黏性土现有的下加载面修正剑桥模型破坏应力比为定值,不能反映土体在不同应力状态下强度特性存在差异这一问题,为了更好地反映应力状态变化对破坏应力比的影响,基于三剪统一强度准则确定超固结土的破坏应力比。在此基础上得到了超固结土的弹塑性本构模型,该模型的特点是能够描述土体受力时的中间主应力效应、应力区间效应、拉压差影响。基于该本构模型作了常规三轴压缩条件下超固结饱和黏性土的模拟结果与试验结果对比,结果表明:该模型能够反映不同超固结比下土的变形、剪胀、孔隙水压力变化特性和规律。利用该本构模型模拟了排水和不排水条件下的真三轴压缩试验,得到了相应的应力-应变特性曲线。     

复杂加载条件下的砂土本构模型

试验表明,饱和砂土的应力应变关系具有显著的密度以及压力依存性,上述两点构成了描述砂土静力加载下变形特性无法忽视的因素. 此外,在循环加载等复杂加载作用下,砂土还会表现出明显的应力诱导各向异性以及相变转换特性. 基于在e--p空间中存在唯一的临界状态线这一基本假定,通过在e--p空间中引入当前状态点与临界状态线的距离R来作为反映密度与压力依存特性的状态参量, 将变相应力比以及峰值应力比表达为状态参量的指数函数,将上述应力比参量引入到统一硬化参量中可准确地反映初始状态下围压、密度 对于单调加载下应力应变关系的影响规律,能描述砂土剪缩、剪胀,应变软化、硬化等特性. 采用非相关联流动法则,p--q空间中采用水滴型屈服面,塑性势面为椭圆面,松砂在单调加载下的静态液化现象也可描述. 为反映循环加载下塑性体积应变的累积特性以及塑形偏应变的滞回特性,在循环加载下将状态参量R表达为应力比参量,并在硬化参数中引入描述应力诱导各向异性特性的旋转硬化部分,所提模型可有效地描述循环加载下剪切模量的衰减特性、刚度衰化性质、强度减小特性,在不排水约束作用下,则会产生往返活动性现象. 通过一系列的模型模拟与试验结果对比,验证了本构模型的有效性及适用性.      

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

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

饱和超固结黏性土的三剪弹塑性本构模型研究

针对饱和超固结黏性土现有下加载面修正剑桥模型中破坏应力比为定值、土体黏聚力为零,以及不能准确反映不同应力状态下土的强度差异这些问题,基于三剪统一强度准则以及应力坐标平移法得到了扩展破坏应力比,其特点是能更好地反映应力状态变化以及土体黏聚力的影响。在此基础上提出了饱和超固结黏性土的三剪弹塑性本构模型,该模型的特点是能描述土体受力时的中间主应力效应,应力区间效应和拉压差影响,同时也能更好地考虑土体黏聚力的影响。基于该模型对ABAQUS软件进行了二次开发,并利用其模拟了饱和超固结黏性土在排水和不排水条件下的真三轴和常规三轴压缩试验特性。对常规三轴压缩条件下土体力学特性作了模拟和试验结果对比。结果表明所提模型能很好地反映不同超固结比下土体的变形、剪胀、孔隙水压力变化特性。     

Effect of dilatancy on the strain localization of water-saturated elasto-viscoplastic soil

It is well known that geomaterials such as soils exhibit an increase in volume during shearing deformation, referred to as dilatancy. Dilatancy is a typical property of such granular materials as soils and is closely related to changes in the microstructure. Normally consolidated clay exhibits negative dilatancy or contractancy, namely, a decrease in volume during shearing. On the other hand, overconsolidated clay shows positive dilatancy, namely, an increase in volume during shearing. The aim of the present paper is to study the effects of the microstructure, such as dilatancy and permeability, on the strain localization of water-saturated clay using an elasto-viscoplastic constitutive model. Based on the non-linear kinematic hardening theory and a Chaboche type of viscoplasticity model, an elasto-viscoplastic model for both normally consolidated and overconsolidated clays is proposed; the model can address both negative and positive dilatancies. Firstly, the instability of the model under undrained creep conditions is analyzed in terms of the accelerating creep failure. The analysis shows that clay with positive dilatancy is more unstable than clay with negative dilatancy. Secondly, a finite element analysis of the deformation of water-saturated clay is presented with focus on the numerical results under plane strain conditions. From the present numerical analysis, it is found that both dilatancy and permeability prominently affect shear strain localization behavior.     

Extended yield condition of soils with tensile yield strength and rotational hardening

An extended yield condition of soils applicable to the range of negative pressure is proposed. A closed surface is formulated in a simple mathematical form not limited to the range of positive pressure. A hardening function of plastic volumetric strain is adopted, based on the linear relation of both logarithms of volume and pressure extended to the negative range of pressure for the isotropic consolidation. Isotropic softening, due to the deviatoric plastic strain rate, and anisotropic hardening, due to rotation of yield surface (rotational hardening) are also incorporated into the extended yield condition.     

Changes in some mechanical properties of a loamy soil under the influence of mechanized forest exploitation in a beech forest of central Belgium

Modification of some soil mechanical properties (penetration resistance and consolidation pressure) induced by vehicle compaction during mechanized forest exploitation was studied in an acid and loamy leached forest soil of the loessic belt of central Belgium. In situ penetration tests and laboratory Bishop–Wesley cell tests were undertaken for the two main soil horizons of a beech high-forest, i.e. the eluvial E horizon (5–30 cm depth) and the underlying clay-enriched B_t horizon (30–60 cm depth). Both undisturbed and wheel-rutted soil areas were studied (E and B_t horizons vs. E_g and B_tg horizons).

Results show that: The experimental overconsolidation pressure of the eluvial reference horizon (E) is about 50 kPa higher than the value calculated from soil overburden pressure; this probably results from suction action during dry periods. The clay-enriched reference horizon (B_t) shows the same trends. In wheel-rutted areas, seven years after logging operations, the E_g horizon memorizes only 14.5% of the wheel induced stress due to forest machinery.

In the compacted B_tg horizon, the experimental overconsolidation pressure represents 96% of the exerted theoretical stresses due to harvesting actions. The good recording of the exerted stresses, after seven years, can be explained by: (1) The B_tg depth which keeps it from seasonal variations i.e. from desiccation–moistening or freeze–thaw cycling; (2) amorphous and free iron accumulation inducing a “glue” effect of the B_tg soil matrix, which could stabilize the soil structure and prevent recovery to initial conditions. These results provide clear evidence that on loessic materials, soil compaction due to logging operations leads to modifications in both physical (bulk density, total porosity) and mechanical (penetration resistance and consolidation pressure) soil properties.     

电解质吸附对水基带电表面黏着的影响

为研究界面黏着,利用表面力仪在高离子强度电解质水溶液中测量负电云母表面之间的法向力.在纯水中,因范德华吸引,两表面跳跃至直接接触(间距0?),分离两表面时测得界面黏着力为-46.7 m N/m.在0.1 mol/L K₂SO₄中,因K₊离子牢固吸附于云母,表面间距稳定于5?,黏着力仅-2.9 m N/m;在0.1 mol/L Ca(NO₃)₂中,云母表面吸附的Ca₂⁺离子产生显著短程水合排斥,但在较低载荷下解吸附,导致两表面直接接触,黏着力高达-40.7 m N/m.加入Ca(OH)₂于0.1 mol/L K₂SO₄仅产生微弱短程排斥,黏着状态几乎与纯K₂SO₄水溶液中相同.聚电解质PCE从0.1 mol/L K₂SO₄中吸附于云母表面且诱导远程空间位阻排斥,但在中等载荷下解吸附,导...     

复杂加载下混凝土的弹塑性本构模型

混凝土材料在不同应力路径下或复杂加载条件下会表现出差异性显著的应力应变关系,在小幅循环加载条件下,其应力应变关系会表现出类似于弹性变形的滞回曲线.在不同应力水平下,混凝土的应力应变关系以及破坏特性都具有静水压力相关特点,即随着静水压力增大,各向异性强度特性弱化.此外,混凝土受压及受拉破坏机理不同,因而对应于混凝土硬化损伤亦有不同,即可分为受压硬化损伤,受拉硬化损伤及两者的混合硬化损伤类型.基于Hsieh模型,对该模型进行了三点改进.（1）针对小幅循环加载下混凝土无塑性变形的试验规律,而模型中在应力水平较低的循环加载条件下始终存在塑性变形的预测问题,采用在边界面模型框架下,设置了应力空间的弹性域,初始屈服面与后续临界状态屈服面几何相似的假定.（2）基于广义非线性强度准则将原模型采用变换应力方法将其推广为三维弹塑性本构模型,采用变换后模型可合理的考虑不同应力路径对于子午面以及偏平面上静水压力效应形成的影响,并避免了边界面应力点奇异问题.（3）分别对拉压两种加载损伤模式建议了相应的硬化参数表达式,可分别用于描述上述加载中产生的应变软化及强度退化行为.基于多种加载路径模拟表明:所建立的三维弹塑性本构模型可合理地用于描述混凝土的一般应力应变关系特性.      

非比例循环塑性内时本构关系研究

在Valanis的内时本构理论的框架中,引入内结构张量以反映由于非比例加载而引起金属材料的附加等向强化及异向强化效应,同时提出材料强化程度的度量采用沿路径法线方向的塑性应变分量来描述.这些考虑的有效性已经通过用所建模型对304不锈钢材料在一些典型非比例循环加载路径下的响应进行的理论预测得到了验证;将该模型应用于U71Mn材料室温单轴棘轮行为描述中,结果显示内结构张量的引入不仅能较好地反映应变控制下的非比例附加效应,而且也能较好地反映应力控制下塑性应变的累积及变化率.     

Anisotropic hardening in single crystals and the plasticity of polycrystals

Based on the dislocation structures developed during plastic deformation, an anisotropic hardening law is developed to describe the latent hardening behavior of slip systems under multislip. This theory incorporates the concept of isotropic hardening, kinematic hardening, and the two-parameter representation; it automatically includes the strength differential between the forward and reversed slips and between the acute and obtuse cross slips. The self-hardening modulus of a slip system is found to be “associated” with the latent hardening law involved, and, based on some experimental evidence, two specific sets of self-hardening modulus are suggested. An important feature of this associated modulus is that the slip system with a soft latent hardening (e.g., the reversed system with a Bauschinge effect) will have an enhanced self-hardening modulus. This newly developed hardening law, together with its associated latent hardening moduli, is then applied to examine the strain-hardening behavior of a polycrystal. Although crystals with a stronger latent hardening will, in general, also lead to a stranger strain-hardening for the polycrystal, the stress-strain behavior of the polycrystal using the kinematic hardening law of single crystals is found to be not necessarily softer than that using the isotropic hardening law. Within the range of experimentally measured latent hardening ratio of slip systems, the anisotropic theory is also used to calculate the motion of yield surface of a polycrystal. The general results, employing four selected types of anisotropic hardening, all show the essential features of experimental observations by Phillips and his co-workers. The application is highlighted with a reasonably successful quantitative modeling of initial and subsequent yield surfaces of an aluminum.     

An anisotropic elastic–viscoplastic model for soft clays

Experimental evidences have shown deficiencies of the existing overstress and creep models for viscous behaviour of natural soft clay. The purpose of this paper is to develop a modelling method for viscous behaviour of soft clays without these deficiencies. A new anisotropic elastic–viscoplastic model is extended from overstress theory of Perzyna. A scaling function based on the experimental results of constant strain-rate oedometer tests is adopted, which allows viscoplastic strain-rate occurring whether the stress state is inside or outside of the yielding surface. The inherent and induced anisotropy is modelled using the formulations of yield surface with kinematic hardening and rotation (S-CLAY1). The parameter determination is straightforward and no additional experimental test is needed, compared to the Modified Cam Clay model. Parameters determined from two types of tests (i.e., the constant strain-rate oedometer test and the 24 h standard oedometer test) are examined. Experimental verifications are carried out using the constant strain-rate and creep tests on St. Herblain clay. All comparisons between predicted and measured results demonstrate that the proposed model can successfully reproduce the anisotropic and viscous behaviours of natural soft clays under different loading conditions.