有限元固结分析中若干问题探讨.pdfVIP

东南大学学报 (哲学社会科学版·高等教育研究) 2008年第3期 有限元固结分析中若干问题探讨 介玉新 [摘 要]在用基于比奥回结理论的有限元法进行固结分析时，涉及平衡方程和连续方程的有限元格式构 造，荷载、位移和水头、流量等边界条件的处理，以及静孔压和超静孔压的计算，等问题，本文对这些问题进行 澄清，重点探讨有限元格式中以水头为未知数和以孔压为未知数两者之间的关系。对计算中所涉及的土的本 构模型，尤其邓肯一张模型中关于泊松比的限制条件也进行了初步分析。 [关键词]固结；有限元；比奥固结理论；邓肯一张模型 一、引言 比奥固结方程是土力学数值计算的基本方程。有限元法是求解比奥固结方程的强有力的数学工具。比 奥固结方程可以看作是平衡方程和连续方程的联立，以平面应变为例，可以写为： 警+努+E-0 1 警+曹迫ay+‰+E—o } (1a) 一竞未(乜差)一瓦1面a(足，器)+豪+曼罄=oJ 或 警+努+塞+R—o 1 鲁+等+考+L—o } (1b’ ．瓦1五8(也差)+去昌(五，考)+鲁一凳警=oJ 其中户为孔隙水压力，按弹性力学符号规定，压应力为负，与常规土力学规定的符号相反；元一‰IIl，h为水头。 在渗透系数不随坐标位置变化时，上式中的连续方程也可以写为 一麦(屯誊+五，孝)+象+蔓罄一。 c2a， 或 麦(屯誊托孝)+象一yS_。,aap一。 ㈨， 利用物理方程和几何方程，．-7以把式(1)进一步表示成位移分量和孔压(或水头)表示的形式‘卜51，这是位 移法有限元格式构造的基础。求解上述偏微分方程组，除了物理方程、几何方程外·还需要与平衡方程有关 的荷载和位移边界条件，以及与渗流有关的水头(孔压)和流量Eta-条件。 太沙基一伦杜立克理论也可以进行固结分析，但它假定主应力之和在固结过程中不随时间变化，显然是 不严密的。它的优点是能够把孔压计算和应力应变计算分离开来，减少计算工作量。这在计算机技术发展早 期还是右一定意义的，但在今天已没有太大意义。 介玉新清华大学水利水电工程系 [基金项目]国家重点基础研究发展计划资助(2007CB714102) ·249· 二、求解比奥固结方程的有限元格式 有限元格式可以用里兹法推导，也可以用伽辽金法推导。全量形式的有限元格式为 [K]‘{艿)‘一[K。]￡{^)。一{RF)： 1 ，。、 一[K。]打{甜t一[K，]‘{，J‘一[Kp]。{五)‘一{R。)：j 前[K]‘{艿)‘一[KdP{户)‘={RF)：1 (3b) “[K。]订{各)‘一[KI]。{户)‘一[Kp]‘{多)‘一{R。)；f 其中 EB]rED][B]dA [K]k|1 [K]峰lI[B]T{M)[N]dA [Kl卜Ⅱ陋]T[k-li-B，]dA吲一刘E-N]dA {RF)：一J[N]T{F)ds+Ⅱ[N]T{F}dA--J][B]T{M)‰ydA {R。)z—l[丙]Tu．ds{RF)；=J[N]T{F)ds+JJ[N]T{F)dA 0)T {R口)；一n丙]ru．ds--JJ-EB．]T[足]㈥dA{M)一{11 L)r为体力。对四节点四边形等参元来说， {F)={E己)