论固流耦合渗透固结数学模型建立.docVIP

论固流耦合渗透固结数学模型建立摘 要：基于多相连续介质力学的理论，把土体抽象为叠合连续体，建立了有限变形下固-流耦合渗透固结问题的数学模型，模型中反映了固-流相的相互耦合作用。 ? 关键词：固流耦合； 多相介质； 有限变形； 数学模型? 中图分类号：O14文献标识码：A文章编号：1672-3198（2009）01-0284-02?? 1 基本假设? 在以固相为宿主相的三相介质固结问题中，作如下假设：? (1)平衡及临界状态以前，固相骨架为准静态的；? (2)流相在固相骨架中的渗流服从达西定律；? (3)固相骨架处理为均质各向同性体；? (4)流相为理想流体。 ? 2 有效应力原理? 以固相为宿主相的骨架的变形主要由有效应力控制，有效应力原理表示为? σ=σ′-pm（1）? 其中，σ为总应力，σ′为有效应力，都是以拉应力为正，p为流体应力，以压应力为正。许多学者提出了适用于岩土修正的有效应力原理? σ=σ′-am（2）? 其中，α为修正系数。? 3 固相骨架的应力平衡方程? 由假设1，忽略固相骨架的惯性力，利用多相介质的动量守恒定律，得到固相骨架的有效应力表示的平衡方程为：? ?gσ?se?+ρsbs-grad（p）=0（3） ? 式中，σ?se?为固相骨架的有效应力，以拉应力为正，P为平均孔隙压力，其表达式为? p=sfpf=sgp（4）? sg=sf=1（5）? 式中：sf，sf分别为水相和气相流体的饱和度，而pf，ps分别为岩土中水相压力和气相压力。方程（3）是以固相骨架为脱离体建立的平衡方程，σ?se?实际上是有效应力。由于现实状态，通常都是需要求解的状态，一般知道的是初始状态的边界条件和状态，对现实状态的初始条件和状态是未知的，是需要求解的。为了求解，还需要把方程（3）转化到初始状态下的物质描述中。? 假设体力和面力在物体的变形过程中保持不变，可以用?Lagrange?应力得到平衡方程：? T?seij，j?+ρ?s0?b?s0i?-p?,j?=0（6）? ?Lagrange?应力张量是非对称的应力张量，使用起来很不方便，把上式变换到?Kirchhoff?应力所表示的平衡方程为：? ??XjS?seki??xj?Xk+ρ?s0?b?soc?-p?，i?=0（7）? 式中，Xi为?Lagrangian?坐标，Xj为?Eulerian?坐标；ρs和ρ?s0?分别为现实构形和初始构形的固体质量密度；b0，b?s0?分别为现实构形和初始构形的外体力密度；S?se?表示固相有效的?Kirchhoff?应力。? 4 在变形多孔介质中的流相控制方程? 流体渗流运动是由流体流动的连续方程（质量守恒）、流体状态方程、流体渗流方程组成。? 假设渗流速度满足达西定律：? V=-KUgrad（p）（8）? 由多相介质的质量守恒定律，流相的连续方程为：? ρ?gm+ρmd?iv?（Vm）=?C?)m（9）? 4.1 气相控制方程? 令式（9）中m=g ，便得到气相流动的连续方程：? ρ?gg+ρgd?iv?（Vg）=?C?)g（10）? 式中：?C?)g为气体的质量增加速率，它可以反映相之间的相互作用。由V?sg?的物理意义，有：Vg=Vs+V?sg?=Vs-Kgρg? k?rg?ugBg+RsfkrfufBfgrad（pg） （11）? 代入上式（10），有：? ρ?gg+ρgd?iv?（Vg=Vs+V?sg?）= Vs-Kgρgk?rg?ugBg+RsfkrfufBfgrad（pg）=?C?)g（12）? d?iv?Kgρgk?rg?ugBg+RsfkrfufBfgrad（pg）-θ?gS- ρ?ggρg+?C?)gρg=0（13）? 其中：? θ?gS=V?s?k，k?=??t（m?Te）= ??t?u?sx?x+?u?sy?y+?u?sz?z? 为固相骨架的体积变形。? 4.2 液相控制方程? 令式（10）中m=l，便得到液相流动的连续方程：? ρ?gl+ρld?iv?（Vl）=?C?)l（14）? 式中：?C?)l为液相的质量增加速率，它可以反映相之间的相互作用。? 由V?sl?的物理意义，有：? Vl=Vs=Vsl=Vs-Klk?rl?ρlulBlgrad（pl）（15）? 代入式(14)，有：? ρ?gl+ρld?iv?Vs-Klk?rl?ρtulBlgrad（pl）=?C?)l（16）? 化简