三维Biot固结理论的一种张量形式有限元算法.pdfVIP

第43卷第9期山西建筑Vol.43No.9

2017年3月SHANXIARCHITECTUREMar.2017•83•

文章编号：1009-6825(2017)09-0083-02

三维Biot固结理论的一种张量形式有限元算法

张译心左博文

(东北林业大学土木工程学院,黑龙江哈尔滨150040)

摘要:基于Biot的假定，从连续介质力学的弹性方程开始，利用张量推导三维Biot本构方程,再根据Darcy定律推导出了控制方

程，并对三维Biot问题控制方程进行空间离散和时间离散，给出空间有限元格式以及时间差分格式，便于后续计算机求解。

关键词：三维Biot,Darcy定律,有限元，数值计算

中图分类号:TU431文献标识码:A

〇引言sij=Diju(TU⑷

Biot固结理论是岩土工程力学领域中的重要课题，是研究饱其中，为弹性系数，满足式(5)。

和土体的目前公认的流固耦合机理的理论基础。流固耦合分析i=j=k=l

是进行土工，特别是土与结构相互作用问题、开挖与填筑的施工

过程的模拟等问题深人研究的主要途径与发展方向[1]。张量理(5)

论是解决建筑学和岩土力学的一个有力的数学工具。

自1941年，Biot[2]首次提出基于严格固结机理推导的能准确

反映孔隙压力消散与土骨架变形之间耦合作用的真三维固结理^0,other

论以后，许多学者对这些方程进行了研究，并利用该理论解决了其中，£■，分别为杨氏模量，剪切模量，Poisson比，由于各向

大量的岩土工程问题。国内相关文献所介绍的Biot固结理论一

同性的基本假设，它们只有两个是独立的，并满足:C=2(二)。

般都属于Biot(1941)提出的形式，但是以张量形式推导的很少，

更重要的是近年鲜有专门研究其有限元方程形式的文献。本文1.5孔隙水压力相关项

从连续介质力学的弹性方程开始，利用张量推导三维Biot本构方因为各向同性的假设意味着对称，因此空隙水压力增加不改

程和控制方程，为之后的计算机计算求解提供理论参考。变剪切变形。再根据线弹性假设，扩展Hooke定律，得：

1Biot本构模型Sij=DijU〇-u+S^(6)

1.1Biot本构模型基本假设

其中，九，分别为空隙水压力增量，与空隙水压力相关的物

Biot本构模型基本假设为：

理常数。

1)满足各向同性;2)0—6：在固结完成时可逆;3)0—^满足

除应力外，单位体积水的增量0对于应力为线性关系。由各

线性关系;很小;5)水不可压缩;6)水中可以有气泡