一维弹黏塑性固结模型研究论文

《一维弹黏塑性固结模型研究论文》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、一维弹黏塑性固结模型研究论文摘要：建模时考虑固结压缩的分段性，引入参考应力状态概念、弹黏塑性屈服准则以及一些其它观点，同时结合殷建华3等提出有效应力、应变和蠕变速率的唯一性原则，建立了一个原状土的一维弹黏塑性固结模型。计算表明，该模型能适用于一维条件下任何加载方式的固结模拟，能描述表观前期固结压力的应变率效应、次固结引起的表现前期固结压力增加等一些已被试验证实但不能被太沙基固结理论反映的现象。关键词：固结弹黏塑性屈服准则表观前期固结压力1967年，Bjerrum提出了一个描述土压缩时间效应的模型1。

2、该模型首次提出次固结的表观前期固结压力效应，并用时间线概念成功解释次固结的黏滞性行为.freel的瞬时压缩指加荷瞬间土颗粒间的应力增量就等于所施加荷载，并由此而产生的瞬时沉降，不考虑水动力滞后的影响，因此按Bjerrum的定义，怎样构造饱和软粘土瞬时时间线是个难点；③Bjerrum的时间线是用加载持续作用时间定义的，在多级加载条件下失去参考意义，因此根据它推导的固结模型仅能适用于单级加载条件。针对Bjerrum时间线概念的上述局限性，众多的学者(Ha2,3在前人工作的基础上，给出了等效时间的物理解释

3、和数学定义，指出“对特定加载历史的一个状态点，其蠕变速率等于相同应力下从参考点估算的等效蠕变时间所计算的蠕变速率”，并阐明了等效时间与蠕变速率的关系，即对给定状态点，依据蠕变速率求得的时间是等效的，同时在试验的基础上提出有效应力、应变和蠕变速率的唯一性原则。可是绝大多数原状土都有一定的结构性，其压缩曲线有明显的分段性4，上述等效时间概念及有效应力、应变和蠕变速率的唯一性原则仅仅描述正常压缩阶段应力-应变关系的特性，不能考虑再压缩阶段在固结过程中的行为，也不能反映表观前期固结压力对固结行为的影响。在本

4、文中，考虑结构性饱和土体固结压缩过程的分段性，同时结合殷建华等提出的有效应力、应变和蠕变速率的唯一性原则，创建了一个与时间相关的一维固结模型。由于在构造模型中提出参考应力状态的概念和一个新的屈服准则，这些措施使本文模型适用于任何加载条件（单级、多级和连续加载条件），并能描述整个固结进程中的孔隙水压及应力、应变随时间的变化规律。1建模基本思想根据结构性黏土的变形机理，可将压缩曲线的第一段视为弹性变形阶段，第二段视为弹黏塑性应变硬化阶段。在弹黏塑性阶段，有效应力、应变和黏塑性蠕变速率的相互关系具有唯一性

5、。与文献2不同，本文用蠕变速率来描述时间线，如图1所示，P0D为加载持续为24h的正常压缩线，将其定义为参考时间线，设其上的蠕变速率为，根据文献2，，Cεa为次固结系数(Δε/Δlogt)。如果知道时间线AC的蠕变速率，则在任意荷载下，从参考时间线到时间线AC的蠕变变形P0C等于。大量实验证实，土的屈服应力(亦即表现前期固结压力)依赖于应变速率3。表观前期固结压力与变形的速率相关的根源在于土骨架有黏滞性效应，因此可认为表观前期固结压力主要是与黏滞性的蠕变速率相关。当实际的变形速率等于过当前应力-应变

6、状态点的时间线上的黏塑性蠕变速率时，达到初始屈服状态，这时弹性变形速率已极小，因此可以认为达到（室内常规固结压缩试验条件下）表观前期固结压力时的应变率与参考时间线上的蠕变速率相等，如图1所示。在弹性变形阶段，固结系数大，孔隙水压消散快，有效应力小于初始屈服应力，实际的变形率大于过当前应力-应变状态点的时间图1时间线与参考时间线线上的黏塑性蠕变速率；当实际的变形速率等于过当前应力-应变状态点的时间线上的黏塑性蠕变速率，达到初始屈服状态；在弹黏塑性变形阶段，总应变率等于弹性应变率和黏塑性蠕变速率（下文简

7、称为蠕变速率）之和，有效应力、应变和蠕变速率的相互关系具有唯一性2，实际的蠕变速率等于过当前有效应力-应变状态点的时间线上的蠕变速率。为了记录变形过程中应力-应变状态的变化，本文将当前应力-应变状态点在参考时间线上的“投影”，称为参考状态点，它表示当前应力状态的参考应力状态，如图1中初始状态点O的参考状态点为P0。对单级加载，初始状态和弹性变形状态在参考时间线上的参考状态点相同，将此参考状态点称为初始参考状态点，该点的有效应力等于室内常规压缩试验条件下测得的表观前期固结压力（亦即该试验条件下土骨架的

8、初始屈服应力）。在弹黏塑性变形阶段，变形处于黏塑性应变硬化状态，随着新的屈服产生，参考应力状态（点）也不断发生变化。2一维弹黏塑性团结模型推导初始应力状态对土的应力-应变性状有很大影响，而且加载方式不同，初始应力状态会发生变化。下文分两种情况来推导一维弹黏塑性本构方程。2.1单级加载条件下的一维弹黏塑性固结模型推导2.1.1弹性阶段变形计算在有效应力达到相应于表观前期固结压力的屈服应力之前，变形是弹性的，应力-应变关系为ε=Cεrlog(σ/σ0)，对此式求时间导数：