比奥固结理论

太沙基固结理论只在一维情况下是精确地，对二维、三维问题并不精确。比奥（Biot） 从较严格的固结机理出发推导了准确反映孔隙压力消散与土骨架变形相互关系的三维固结 方程，一般称为真三维固结理论，而将太沙基三维方程称为拟三维固结方程。介绍饱和土体 固结的比奥理论。一比奥固结方程一）三维问题1. 平衡方程在土体中取一微分体。若体积只考虑重力， z 坐标向上为正，压力以压为正，则三维平衡微分方程为式中，丫为土的重度,上式也可以写为应力为总应力。x + Qx+k =QyQzQTQTQT亠+ +yz :QxQyQzQtQTQcxz + +=QxQyQz=0Y时.= f 其中ar = =Qx00仁X yztyzTzxTxy00QzQ QQz QyQQz0QQxQQyQQx0式中f三个方向的体积力。2. 有效应力原理根据有效应力原理，总应力为有效应力与孔隙压力u之和，且孔隙水不承受剪应力，用矩阵表示：o=l+ M u平衡方程可以写为bP br + m u = f展开即为db比 Qtdu八dxdydzo xdtdbdtduxy +y + 埒 += 0dxdydzo ydtdtdb duxz + 斗 +z += Ydxdydzo z式中学、学、du实际上是各方向的单位渗透力，此式是以土骨架为脱离体建立的 o x o y o z平衡微分方程。3. 本构方程利用本构方程中的物理方程b=D 式 1可将式中的应力用应变来表示。匕奥最初假定土骨架是线弹性体，服从广义胡克定律，则d 为线弹性矩阵，上式 1 可以写成bx=2G宀+訂 d_2o v x丿by=2Gf佥匕+匕式3Ii_2D v y 丿bz=2G 8 +8jl-2o v z丿t = g , t = g , t = gyz yz xz xz xy xy式中，G和u分别为剪切模量和泊松比。其实，物理方程并不一定要限于弹性，也可推广到线弹性体，这时d为弹塑性矩阵。4. 几何方程 再利用几何方程，将应变表示成位移。在小变形的假设下，几何方程为8 = -dco式 2式中，为位移分量。将上式展开，- Q(QQ )s =xy =-y+zxQxyzLQzQy丿Q(QQs =yy =-x-+zyQyxzL QzQx丿QqQs =zy =-x-+yzQzxyl QyQx丿应力应变符号在土力学中习惯以压为正，以拉为负，故式 2 与一般弹性力学中几何方程 的符号相反。5. 连续性方程三维固结连续性方程6. 固结微分方程 将式2代入式3，再代入式1，就得出以位移和孔隙压力表示的平衡微分方程-qt d qw + bb m u = f 式 4对于弹塑性问题，方程展开是较复杂的，这里只给出弹性问题的展开方程：-GV2 -xG Q | Q Q Q | Qu _x + + z 1 + =0QyQz 丿 Qx1 - 2u Qx ( QxG-GV2 -yQ | Q1 - 2u Qy、QxQQ | Qux +4 + z I += 0QyQz 丿 Qy-GV2 -zG Q | Q1 - 2u Qz、QxQQ | Qux +卉 +z I += -yQyQz 丿 Qz将体积应变利用式 2 以位移来表示。即s 二m = -mtqv则连续性方程Q|v = -K V2U （式5）成为w-Mt q+V2U = 0 式 6Qty其展开式为Q+ y +Q)K V7 n+ V 2u = 0y这就是以位移和孔隙压力表示的连续方程。 饱和土体中任一点的孔隙压力和位移随时间的变化，须同时满足平衡方程式 4和连续方 程式 6。将两式联立起来，并对时间取差分格式，即式7-sT d9co+ bb M u _ f -Mt 6cd-K V2u = 06ty式7 便是比奥固结方程。它是包含4个偏微分方程的微分方程组，也包含4个未知变量u、，它们都是坐标x、y、z和时间t的函数。在一定的初始条件和边界条件 xyz下，可解出这4 个变量。式 7 为联立方程，反映变形方程与渗流的耦合，也称为流固耦合 其中，平衡方程的第一项表示发生的位移所对应的力，第二项表示当前孔压所对应的力。它 们的和与外荷载平衡。连续性方程的第一项表示单位时间内位移改变所对应的体积变形，第 二项表示孔压变化所引起的渗出水量。力的平衡中又有孔压的贡献，水量平衡中又有变形的 贡献，相互耦合。需要特别指出的是：以上方程中所讲的孔隙水压力u都是指超静孔隙水压力，即荷载引 起的孔隙水压力增量。（二）轴对称问题轴对称问题的平衡微分方程为doQtq -qr +rz + ro _ fdrdzr r-rzdtdqt_ +z + rz _ fdrdzr z连续性微分方程为6s _ dqqdqv _ r + r 6t6rr6z将它们联立起来，并将应力用径向位移和竖向位移表示，流量用孔压u表示，成为以 rz位移和孔压为未知变量的联列方程组，即轴对称问题的固结微分方程。（三）二维问题对于平面变形问题，比奥固结方程仍可写为式 7，只是式中ddxddzddzdxm =11w_wxd 2d 2V2 _ + wdx dz0zz相应的应力应变为=ox ozTxz式7 在二维问题中的展开式为-GV 2 +x-GV2 -z小8x8zYxz= -m TbG 08v1 - 2u 0x+ 0u = 00xG08v1 - 2u 0z0u+=Y0Z08 K严+-0tYWV 2U = 0要解上述微分方程组，在数学上是困难的。对于轴对称和平面应变中某些简单情况，已 有人推导出了解析回答，并用以分析固结过程中的一些现象。但是对于一般的土层情况，边 界条件稍微复杂一些，便无法求得解析解。因此，从1941 年建立比奥方程以来，一直没有 在工程中广泛应用。随着计算技术的发展，特别是有限单元法的发展，真三维固结理论才重 现生命力，并开始用于工程实践。