一维弹黏塑性固结模型研究.doc

1 一维弹黏塑性固结模型研究摘要: 建模时考虑固结压缩的分段性, 引入参考应力状态概念、弹黏塑性屈服准则以及一些其它观点, 同时结合殷建华[3] 等提出有效应力、应变和蠕变速率的唯一性原则,建立了一个原状土的一维弹黏塑性固结模型。计算表明,该模型能适用于一维条件下任何加载方式的固结模拟, 能描述表观前期固结压力的应变率效应、次固结引起的表现前期固结压力增加等一些已被试验证实但不能被太沙基固结理论反映的现象。关键词:固结弹黏塑性屈服准则表观前期固结压力 1967 年, Bjerrum 提出了一个描述土压缩时间效应的模型[1] 。该模型首次提出次固结的表观前期固结压力效应, 并用时间线概念成功解释次固结的黏滞性行为, 但其时间线概念有若干缺陷:①瞬时时间线物理意义不清, 忽略了应力- 应变关系的时间效应, 因为土体同其它工程材料一样, 土骨架的 2 屈服应力随应变速率的增加而增加; ② Bjerrum 的瞬时压缩指加荷瞬间土颗粒间的应力增量就等于所施加荷载, 并由此而产生的瞬时沉降,不考虑水动力滞后的影响,因此按 Bjerrum 的定义, 怎样构造饱和软粘土瞬时时间线是个难点; ③ Bjerrum 的时间线是用加载持续作用时间定义的,在多级加载条件下失去参考意义, 因此根据它推导的固结模型仅能适用于单级加载条件。针对 Bjerrum 时间线概念的上述局限性, 众多的学者(Hawley Borin, 1973; Mesri Rokhsar, 1974; Magnan et al., 1979; Christie Tonks, 1985) 提出了不同的修正。实际上,各种修正方案的区别仅在于对时间线物理意义的阐述,但对时间线的描述方式基本达成共识,即不能用独立的时间变量(如加载持续时间)来定义时间线, 而应用间接时间变量, 如等效时间或蠕变速率来定义时间线, 不过各位学者对等效时间或蠕变速率所赋予的意义并不相同。 Yin Graham[2,3] 在前人工作的基础上, 给出了等效时间的物理解释和数学定义,指出“对特定加载历史的一个状态点, 其蠕变速率等于相同应力下从参考点估算的等效蠕变时间所计算的蠕变速率”,并阐明了等效时间与蠕变速率的关系,即对给定状态点,依据蠕变速率求得的时间是等效的, 同时在试验的基础上提出有效应力、应变和蠕变速率的唯一 3 性原则。可是绝大多数原状土都有一定的结构性, 其压缩曲线有明显的分段性[4] , 上述等效时间概念及有效应力、应变和蠕变速率的唯一性原则仅仅描述正常压缩阶段应力- 应变关系的特性,不能考虑再压缩阶段在固结过程中的行为,也不能反映表观前期固结压力对固结行为的影响。在本文中,考虑结构性饱和土体固结压缩过程的分段性, 同时结合殷建华等提出的有效应力、应变和蠕变速率的唯一性原则, 创建了一个与时间相关的一维固结模型。由于在构造模型中提出参考应力状态的概念和一个新的屈服准则,这些措施使本文模型适用于任何加载条件( 单级、多级和连续加载条件) ,并能描述整个固结进程中的孔隙水压及应力、应变随时间的变化规律。 1 建模基本思想根据结构性黏土的变形机理,可将压缩曲线的第一段视为弹性变形阶段, 第二段视为弹黏塑性应变硬化阶段。在弹黏塑性阶段, 有效应力、应变和黏塑性蠕变速率的相互关系具有唯一性。与文献[2] 不同, 本文用蠕变速率来描述时间线, 如图 1 所示, P0D 为加载持续为 24h 的正常压缩线, 将其定义为参考时间线, 设其上的蠕变速率为, 根据文献[2] ,,Cεa 4 为次固结系数( Δε/Δ logt) 。如果知道时间线 AC 的蠕变速率, 则在任意荷载下,从参考时间线到时间线 AC 的蠕变变形 P0C 等于。大量实验证实, 土的屈服应力( 亦即表现前期固结压力) 依赖于应变速率[3] 。表观前期固结压力与变形的速率相关的根源在于土骨架有黏滞性效应, 因此可认为表观前期固结压力主要是与黏滞性的蠕变速率相关。当实际的变形速率等于过当前应力- 应变状态点的时间线上的黏塑性蠕变速率时, 达到初始屈服状态, 这时弹性变形速率已极小, 因此可以认为达到(室内常规固结压缩试验条件下) 表观前期固结压力时的应变率与参考时间线上的蠕变速率相等,如图 1 所示。在弹性变形阶段,固结系数大,孔隙水压消散快,有效应力小于初始屈服应力, 实际的变形率大于过当前应力- 应变状态点的时间 5 图1 时间线与参考时间线线上的黏塑性蠕变速率;当实际的变形速率等于过当前应力- 应变状态点的时间线上的黏塑性蠕变速率, 达到初始屈服状态; 在弹黏塑性变形阶段, 总应变率等于弹性应变率和黏塑性蠕变速率(下文简称为蠕变速率)之和,有效应力、应变和蠕变速率的相互关系具有唯一性[2] , 实际的蠕变速率等于过当前有效应力- 应变状态点的时间