bio固结方程的推导及其应用.docxVIP

bio固结方程的推导及其应用 1 以应力及应变张量形式的存 固固界面理论是岩石力学的研究对象，是固固界面机制的理论基础。自1941年,Biot提出固结方程以后,许多学者对这些方程进行了研究,并利用该理论解决了大量的工程问题。 通过对现有文献资料的调研,发现对于Biot固结方程的形式,至今尚未形成统一的认识。Biot固结方程包括4部分,即岩土体单元的平衡方程、固体骨架的本构方程、岩土体单元的连续性方程和Darcy定律。根据这些方程,给定初始条件和边界条件后就可以对问题进行求解。岩土体单元的平衡方程用张量形式可以表示为(对于应力和应变的符号,本文以拉为正): 式中:?为梯度算子;σe为有效应力张量;p为流体压力;α为有效应力系数(标量或二阶张量);F为体力矢量。对上述平衡方程,学者间不存在异议。在弹性条件下,固体骨架的本构方程即为虎克定律,用张量可以表示为 式中:ε为应变张量;σe为弹性应力增量;C为4阶弹性系数张量。对此本构方程,学者间不存在异议。Darcy定律描述水头梯度与渗流速度的关系,一般表示为 式中:q为渗流速度;k为渗透系数;h为水头,对此方程,学者间也不存在争议。对于岩土体单元的连续性方程(下称连续性方程),在固体颗粒和流体都不可压缩的条件下,Biot给出的形式为 式中:εv为固体骨架的体积应变,文献采用了该形式的连续性方程,但部分学者提出了不同的看法,争议的焦点是式(4)的右端是否应该包含有效应力系数α。 Biot是在小变形假设条件下通过引入应变能密度函数,利用偏导数的求导法则推导出连续性方程,虽然此过程中没有用到质量守恒方程,但其推导思路是严密的,按照该思路应该能得到正确的连续性方程。既然不同的学者在式(4)所示的连续性方程上存在分歧,说明Biot在推导过程中可能存在失误。本文将首先回顾Biot推导土体连续性方程的过程,对此过程中可能出现的问题进行探讨,并尝试利用Biot的思路推导出正确的连续性方程。 2 固结方程的方法 为了对流体的连续性方程进行讨论,本节将首先回顾Biot初期推导固结方程的方法,找出其在推导过程中存在的错误。在此基础上,根据Biot的思路重新推出连续性方程,并将该方程与通过质量守恒定律得到的方程进行对比,从而确定固结方程中的连续性方程的正确形式。 2.1 固结基本方程 式(4)所示的连续性方程的是Biot在文献中首次给出的,并在文献中重新给出了推导过程,尽管这两篇文献最后推导出了相同的固结基本方程,但两者的推导过程是有差别的。下面将概述上述两篇文章的推导思路。 文献并没有引入有效应力的概念。采用的基本变量为土体的7个动力分量和7个运动分量。7个动力分量为土体的6个应力(总应力)分量σxx、σyy、σzz、τxy、τyz、τxz以及水压力增量p,7个运动分量为土体骨架的6个应变分量εxx、εyy、εzz、γxy、γyz、γxz以及单位体积土体中水的体积改变量θ。求解固结问题时,所用的基本场变量为土体骨架的位移ux、uy、uz和水压力增量p。 文献中,在小变形的条件下,假设运动分量与动力分量间是线性关系,并且在各向同性条件下剪应力对水的体积变化没有影响,则有(本文符号规定以拉为正,这可能造成本文给出的公式和相关文献中的符号有差异,但这不会影响公式的意义): 式中:E和ν为固体骨架的弹性模量和泊松比;G为固体骨架的剪切模量;H、H′和R为参数。假设土体单元的弹性应变能密度函数为 为了推导公式的方便,假设土体单元仅受到静水压力σ1的作用,此时有σxx=σyy=σzz=σ1,τxy=τyz=τxz=0,则有: 式中:εv为土体骨架的体积变形,即εv=εxx+εyy+εzz。 此时,由式(5)可得 式中:,为土体单元的体积模量。 由式(7)可得 则有: 由式(8)和式(10)可以得到:H=H′。如果令H1=1H′=Kα,K为骨架的体积模量,α为有效应力系数,则式(5)可以表示为 式中:;则表示流体压力引起的流体和固体颗粒的体积变形量。 假设固体颗粒和流体不可压缩,则Q=∞,此时单位土体内流体的体积改变量必然等于经过单位土体边界上流入或流出的流体的体积,设土体单元的渗流速度为q,由式(11)的第7式可得 由此,得到了式(4)所示的连续性方程。 文献得到的固结基本方程和文献是一样的,但其推导过程与文献稍有不同,该文引入了有效应力的概念,并以孔隙率作为有效应力系数,采用的基本变量仍然是7个动力分量和7个运动分量,但此处的7个应力分量的应力部分已经从前文的土体的应力分量换成了土体骨架的应力分量,即用有效应力分量代替了总应力分量,而水压则变成了流体承受的压力,即σ=fp,其中f为孔隙率,p为流体压力;同时,7个运动分量则变为固体骨架的变形和流体的体积应变。求解固结问题时,所用的基本场变量与前文不同,为土体骨架的位