基于混合物理论的饱和介质体应变本构模型

为解决饱和多孔介质的建模问题,采用工程混合物理论建立了饱和多孔介质体积本构理论框架。首先,假定多孔固相与流相基质体积变形功相互独立,采用Terzaghi有效球应力、孔压和流体基质压力作为本构模型的应力状态变量,获得了固相、固相基质和流相基质体应变的余能表达式。其次,根据Lade和de Boer模型试验测试数据,建立了加卸载阶段饱和多孔白塞木立方体流固两相体积本构方程,推导了固相体积切线模量、Biot切线系数和流相Biot切线模量等力学参数计算公式;分析了加载阶段固相体积切线模量、Biot切线系数、流相Biot切线模量等力学参数随Terzaghi有效球应力和孔压的变化规律。最后,根据本文体积本构模型和静力平衡方程建立了饱和多孔介质的1维固结方程,数值分析了饱和多孔白塞木立方体的固结特性,获得了固结度和沉降随时间的变化曲线。研究表明:固相体积切线模量随Terzaghi有效球应力的增大而增大,随孔压u的增大而减小。Biot切线系数介于0.42~0.95之间,随Terzaghi有效球应力和孔压的增大而减小。流体Biot切线模量随Terzaghi有效球应力的增大先减小后增大,随孔压增大而减小。孔压切线系数在大多数情况下小于1.0。考虑固相基质变形时饱和多孔介质的初始孔压不等于外荷载,因此饱和多孔介质在外荷载作用下存在瞬时沉降。本文的建模方法可用于非线性饱和多孔介质的建模和数值分析工作。