非饱和土固结理论的简化.docxVIP

非饱和土固结理论的简化 20世纪90年代初，国内外科学家从不同角度先后导出了非相交土壤变形、孔隙水流和气虚的联合方程，建立了非相交土壤力学的基本框架。组合方程通常包括五个未知变量，即三个位移分量、孔隙水压和孔气压。这个多变量方程的解非常复杂。因此，除了一些科学家在理论上进行了一些计算外，几乎没有使用该理论适用于实际工程的例子。为了将非饱和和土壤结构理论应用于实践，作者试图简化这一理论。在这项工作中，我们将粘土下的雨水与裂缝土的结合进行了初步开发，并进一步发展了非饱和土壤结构理论。 1 有效应力 在非饱和土力学的发展过程中有过两种研究思路, 一是以单变量的有效应力为基础, 另一是以净应力和吸力双变量为基础.对非湿陷性土类, 可以证明在一定条件下这两种研究思路是相通的. 设σm为平均应力,ua和uw分别为孔隙气压力和孔隙水压力.则球应力引起的体积变化可以写为 Δεv=m(Δσm-Δua)+mb(Δua-Δuw)(1) 式中:m和mb为压缩系数. 定义折减吸力ˉs=∫ˉχ(Δua-Δuw)(2) 式中:ˉχ=mb/m.则 (1) 式可以简化为 Δεv=mΔσ′m(3) 式中:有效应力增量Δσ′m=Δσm-Δua+Δˉs;Δˉs=ˉχΔs. 显然, 非饱和土的有效应力作为对变形有效的一个状态变量, 必须满足 (2) 式可以积分的条件, 或系数ˉχ只能是吸力 (ua-uw) 的函数.在针对强度有效的有效应力的研究中, 已得出ˉχ只与饱和度有关的结论, 而针对变形, 则尚未有过研究.但在单调加湿和单调干燥的条件下可积性大体上能保证, 这时, 单、双变量理论是等价的.所以, 在以下的讨论中, 将认为非饱和土有效应力原理是适用的. 一般采用Bishop公式表述非饱和土的有效应力 σ′=σ-ua+χ(ua-uw)(4)由(2)式积分得出ˉs=χ(ua-uw)(5) 则折减系数χ=ˉs/s;s= (ua-uw) 为实际的基质吸力.显然, 当χ是s的函数时, 有ˉχ=χ+s?χ/?s. 可由饱和土压缩-回弹曲线与非饱和土干缩-湿胀曲线的对比试验可求得折减系数χ随吸力或饱和度的变化 (见图1) .设p和s分别为同一孔隙比饱和土所受的压力和非饱和土所受的吸力, 则根据等效原则ˉs=p, 故折减系数χ=p/s, 或写成一般形式 χ=f1(s)(6) 当s与饱和度Sr的关系为已知,χ也可以表示为Sr的函数.由于吸力-饱和度曲线在失水和吸水阶段不同 (见图2) , χ的变化在失水和吸水阶段也就不同. 2 非线性和固结方程的简化 2.1 基本公式 若忽略温度影响及孔隙水中溶解气和孔隙气中蒸汽的流动, 非饱和土固结理论方程组可归纳为: 2.2 初始含气率n0 定义∶孔隙含气率na=[1-(1-ch)Sr]n(14) 则在完全不排气的情况下, 将Boyle定律ρa=ρa0(1+ua/pa) 代入 (9) 式中, 因其等号右边为0, 可得 ua=(na0na-1)pa(15) 式中:初始含气率na0=[1- (1-ch)Sr0]n0;pa为大气压力. 部分排气时, 假设单位时间内的排气量为Δqa, 并定义孔隙气的排气率ξ=Δqa/ρaΔna, 再将Boyle定律代入 (9) 式, 可得孔隙气压力的增量 Δua=-pa+uana(1-ξ)Δna(16) 当排气率ξ为常量时, 将 (16) 式积分可得 ua=[(na0na)(1-ξ)-1]pa(17) 可见, 不排气时,ξ=0, (17) 式就退化为 (15) 式;完全排气时,ξ=1,ua=0.以上假设虽然比较粗糙, 但很多情况下难以确定排气的边界条件.例如, 降雨量较大时, 排水通道被水淹没, 孔隙气只能以气泡形式冒出水面, 其边界条件就难以确定. 2.3 垂向透水系数kw 考虑到球应力引起的体积变化Δεv=-Δn和Δna=?na?Sr?Sr?s(Δua-Δuw)-?na?nΔεv, 将它们代入 (4) 式, 则得到相应的总应力增量 而 (8) 式则变为 μn?uw?t=-??xkwx?h?x-??zkwz?h?z+Sr?εv?t(21) 式中:μ=?Sr/?uw;kwx和kwz分别为水平和垂向透水系数;h=uw/ρwg+z. 对于饱和土, 相应的渗流方程为 mfn?uw?t=-??xkwx?h?x-??zkwz?h?z+?εv?t(22) 式中:mf为孔隙流体的压缩系数.当饱和度达到1时, 只要令μ=mf, (21) 式将自动退化为 (22) 式. 3 本结构模型 (1) 采用本构模型的三轴剪切试验方法,即根据重塑法求解塑性应变增量m 按照前述有效应力原理, 以下公式中的应力均指有效应力.现仍采用双硬化屈服面 F(σ,εpv,εps)=σm1-[ηα(εps)]m-p(εpv)(23) 式中:σm=13(σ1+σ