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

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

土的三重屈服面应力应变模式

经典塑性理论把应变分成弹性应变及塑性应变两部分,不能圆满说明土的塑性变形性质。文中进一步把塑性应变分成若干部分,相应于每一部分有一个屈服面,由此形成多重屈服面的概念。根据土的变形特点,提出三重屈服面模式。三个屈服面分别称为压缩屈服面、剪切屈服面和剪胀压硬屈服面。对砂土和粘土分别建议了屈服函数的具体表达式,并引用两组试验资料加以论证。     

考虑围压效应的冻结砂土动态本构模型研究

为描述主动围压作用下冻结砂土的动态力学特性,通过在朱-王-唐模型的非线性体上串联塑性体,建立了能够考虑围压效应的冻结砂土动态损伤本构模型;分析了损伤参数对应力-应变曲线特征、屈服点、峰值应力和峰值应变的影响规律,基于冻结砂土动力学试验数据确定了模型参数;通过将模型和试验数据进行对比,并对不同试验条件下模型的预测误差进行分析,验证了模型的适用性和准确性。结果表明,损伤参数对应力-应变曲线弹性阶段和屈服点无明显影响,而对塑性阶段和破坏阶段的影响较为显著,本构模型预测的应力-应变曲线与试验结果具有较好的一致性。模型能够预测围压引起冻结砂土塑性阶段占比大和屈服点明显的特征,且能够描述围压对冻结砂土动态强度的增强效应;不同负温和主动围压条件下,模型对峰值应力和屈服强度的预测效果优于峰值应变和屈服应变。     

砂土的应力路径本构模型

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

饱和砂土的循环边界面本构模

饱和砂土在循环加卸载过程中,土体发生显著的组构变化和塑性变形的累积。试验现象分析表明,循环过程中,砂土的首次与再循环塑性模量既有区别又有联系。因此,论文引入考虑组构变化的剪胀内变量,并提出循环塑性模量的关系式,真实描述上述两种现象;进而基于已有本构模型的基础上,建立饱和砂土的循环边界面本构模型。最后,将饱和砂土的模型预测结果与三轴试验结果进行验证对比,得到较好地吻合,这表明该模型能够合理反映饱和砂土循环加卸载的变形行为。     

饱和砂土的循环边界面本构模型

饱和砂土在循环加卸载过程中,土体发生显著的组构变化和塑性变形的累积.试验现象分析表明,循环过程中,砂土的首次与再循环塑性模量既有区别又有联系.因此,论文引入考虑组构变化的剪胀内变量,并提出循环塑性模量的关系式,真实描述上述两种现象;进而基于已有本构模型的基础,建立饱和砂土的循环边界面本构模型.最后,将饱和砂土的模型预测结果与三轴试验结果进行验证对比,得到较好地吻合,这表明该模型能够合理反映饱和砂土循环加卸载的变形行为.     

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

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

剪胀性饱和砂土弹塑性模型

在修正剑桥模型基础上,本文建立的剪胀性饱和砂土弹塑性模型使用相变状态参数描述剪胀性饱和砂土剪胀特性,克服了修正剑桥模型不能直接模拟剪胀砂土力学行为这一局限性。该模型有两方面的改进:一方面,模型将剪胀应力比Md引入剪胀方程;另一方面,模型在塑性功基础上提出用与应力路径无关的硬化参数来替代修正剑桥模型中的塑性体积应变增量。通过试验验证及与修正剑桥模型计算结果对比,结果表明,该模型较适合模拟剪胀性饱和砂土的力学性能,同时也能较好地体现较密实砂土的硬化及软化现象。模型共8个参数,用常规三轴试验就可获取。     

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

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

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

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

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

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

Micromechanical modeling of the elasto-viscoplastic behavior of semi-crystalline polymers

A micromechanically based constitutive model for the elasto-viscoplastic deformation and texture evolution of semi-crystalline polymers is developed. The model idealizes the microstructure to consist of an aggregate of two-phase layered composite inclusions. A new framework for the composite inclusion model is formulated to facilitate the use of finite deformation elasto-viscoplastic constitutive models for each constituent phase. The crystalline lamellae are modeled as anisotropic elastic with plastic flow occurring via crystallographic slip. The amorphous phase is modeled as isotropic elastic with plastic flow being a rate-dependent process with strain hardening resulting from molecular orientation. The volume-averaged deformation and stress within the inclusions are related to the macroscopic fields by a hybrid interaction model. The uniaxial compression of initially isotropic high density polyethylene (HDPE) is taken as a case study. The ability of the model to capture the elasto-plastic stress-strain behavior of HDPE during monotonic and cyclic loading, the evolution of anisotropy, and the effect of crystallinity on initial modulus, yield stress, post-yield behavior and unloading-reloading cycles are presented.     

Elastic–plastic behavior of non-woven fibrous mats

Electrospinning is a novel method for creating non-woven polymer mats that have high surface area and high porosity. These attributes make them ideal candidates for multifunctional composites. Understanding the mechanical properties as a function of fiber properties and mat microstructure can aid in designing these composites. Further, a constitutive model which captures the membrane stress–strain behavior as a function of fiber properties and the geometry of the fibrous network would be a powerful design tool. Here, mats electrospun from amorphous polyamide are used as a model system. The elastic–plastic behavior of single fibers are obtained in tensile tests. Uniaxial monotonic and cyclic tensile tests are conducted on non-woven mats. The mat exhibits elastic–plastic stress–strain behavior. The transverse strain behavior provides important complementary data, showing a negligible initial Poisson's ratio followed by a transverse:axial strain ratio greater than ?1:1 after an axial strain of 0.02. A triangulated framework has been developed to emulate the fibrous network structure of the mat. The micromechanically based model incorporates the elastic–plastic behavior of single fibers into a macroscopic membrane model of the mat. This representative volume element based model is shown to capture the uniaxial elastic–plastic response of the mat under monotonic and cyclic loading. The initial modulus and yield stress of the mat are governed by the fiber properties, the network geometry, and the network density. The transverse strain behavior is linked to discrete deformation mechanisms of the fibrous mat structure including fiber alignment, fiber bending, and network consolidation. The model is further validated in comparison to experiments under different constrained axial loading conditions and found to capture the constraint effect on stiffness, yield, post-yield hardening, and post-yield transverse strain behavior. Due to the direct connection between microstructure and macroscopic behavior, this model should be extendable to other electrospun systems and other two dimensional random fibrous networks.     

剪胀性砂土边界面模型的研究

剪胀性对于砂土,尤其是中密以及密实砂土,是一个非常显著的特性。相变线是剪胀性砂土的特征曲线,能够反映砂土的围压以及初时孔隙比对变形特性的影响。本文在边界面塑性理论的框架内,把相变状态参量引入到剪胀方程以及塑性硬化模量中,建立了一个能够描述砂土剪胀性以及循环特性的本构模型。本模型采用一套参量可以模拟不同初时孔隙比、不同围压、排水(或不排水)条件下单调(或循环)加载的应力-应变特性。验证表明本模型数值计算与试验结果相吻合。     

Numerical Modeling of Hydraulic Hysteresis in Unsaturated Soils

This article presents the implementation of the constitutive model of Wheeler (Geotechnique 53(1):41–54, 2003) for coupling of hydraulic hysteresis and mechanical behavior of unsaturated soils in a fully coupled transient hydro-mechanical finite element (FE) model (computer code UNSATEX) developed by the authors. The constitutive model considers the effects of irreversible changes of degree of saturation on stress–strain behavior and the influence of plastic volumetric strains on the water retention behavior. The mathematical framework and the numerical implementation of the constitutive model are presented and discussed. The FE model is verified and validated against analytical predictions [obtained using the model of Wheeler (Geotechnique 53(1):41–54, 2003] as well as experimental results from the literature involving unsaturated soils undergoing various combinations of drying, wetting, loading, unloading, and reloading paths. Comparison of the results shows that the developed FE model can be used to predict various aspects of the behavior of unsaturated soils under drying and wetting as well as loading and unloading paths. The merits and limitations of the FE model are highlighted.     

A microstructural elastoplastic model for unsaturated granular materials

The homogenization technique is used to obtain an elastoplastic stress–strain relationship for dry, saturated and unsaturated granular materials. Deformation of a representative volume of material is generated by mobilizing particle contacts in all orientations. In this way, the stress–strain relationship can be derived as an average of the mobilization behavior of these local contact planes. The local behavior is assumed to follow a Hertz–Mindlin’s elastic law and a Mohr–Coulomb’s plastic law. For the non-saturated state, capillary forces at the grain contacts are added to the contact forces created by an external load. They are calculated as a function of the degree of saturation, depending on the grain size distribution and on the void ratio of the granular assembly. Numerical simulations show that the model is capable of reproducing the major trends of a partially saturated granular assembly under various stress and water content conditions. The model predictions are compared to experimental results on saturated and unsaturated samples of silty sands under undrained triaxial loading condition. This comparison shows that the model is able to account for the influence of capillary forces on the stress–strain response of the granular materials and therefore, to reproduce the overall mechanical behavior of unsaturated granular materials.     

Overall softening and anisotropy related with the formation and evolution of dislocation cell structures

In this work, a model, based on a representation of the dislocation cell microstructures by a non-local two-phase material with evolving microstructures, is proposed for the elastic–plastic behavior of metals under monotonic and sequential loading. The first phase represents the cell interior and the second one, the cell walls. The evolution of the microstructure is taken into account considering the cell-wall interfaces as free boundaries. Finally, the accumulation within walls of dislocations crossing the cells defines a non-local hardening process. Assuming a piecewise uniform plastic strain field and assuming ellipsoidal cells, the free energy of the system is calculated. The driving and critical forces associated with the plastic flow of the two-phases and the morphology of the cells are established. In a third part, numerical results are presented for monotonic and sequential loading. The results show an overall softening related to the destabilization of the dislocation microstructures which occurs in sequential as well as monotonic paths.     

主应力轴旋转下黏土累积变形的分数阶元件模型

在长期交通载荷作用下土体塑性累积变形本构模型对路基沉降计算至为关键.元件组合模型可以计算岩土体循环累积应变,但现有的各类元件模型未能反映饱和软黏土的主应力轴循环旋转现象.在对饱和软黏土进行等向固结条件下的主应力轴循环旋转加载试验及非等向固结下的循环扭剪试验基础上,将Abel黏壶代替Burgers模型中的Newton黏壶,得到分数阶Burgers模型;利用遗传算法优化循环塑性累积应变的Burgers模型和分数阶Burgers模型的参数,通过对比两组模型的计算值与试验值,发现分数阶模型更适合模拟计算循环载荷下饱和软黏土的累积变形.