地基液化的三维有限元研究.pdfVIP

第五届全国土动力学学术会议论文集 大连 卫·夕98石 地基液化的三维有限元研究 赵剑明 汪闻韶 张崇文 (中国水利水电科学研究院岩土所.北京 ]00044) (天津大学水利系，天淬 300072 摘 要 本文在所提出的结点等价体积和三维结点等价流tI等概念的基础上，建立了三维动力韭d 、、鲤终分近有限元法，并用此对影响地基液化的诸因素进行了研究，得出了具有工程意义的规律性结论， 可供工程场地设计参考. 基本理论 1.1基本原理 本文在前人研究基础上，利用提出的结点等价体积和三维结点等价流量的概念， 引入Biot固结理论，推导了孔隙水压力的产生、消散及地基沉降的理论公式，建立起有 限单元法的固结方程，按有效应力法进行分析计算，在动力分析中引入了履历型 Hardin-Dmevich模型来模拟土的非线性动力本构关系，在动应力全过程或每个时段求得 平均动应变之后，求出此时段末的孔隙水压力值和有效应力状态，以此时段的有效应力 确定新的模量值，作为下一时段迭代计算的初始值。若考虑孔隙水压力的消散扩散的情 况，则应由固结理论和孔隙水压力增长两个方面同时考虑进行分析，如此得出该时段消 散后的孔隙水压力值，进一步求出下一时段的初始模量，进而建立起一套较为完善实用 的研究液化的三维非线性有限元法，以此对地基液化问题做分析研究。本文采用三维八 结点等参元进行分析。 1，动力方程 根据达朗贝尔原理得到分析模型的整体运动方程为 M[]P}+PC+]}[ K[]U{}=R[, ()‘ 式中U{}为分析模型结点的位移列向量;刘和1f分别为结点位移的加速度和速度向 量:R{,}是作用于模型的荷载向量:[M],C[」和K[〕分别为分析模型总的质量矩阵、阻 尼矩阵和刚度矩阵。 2.振动孔隙水压力上升量 对于该量的分析计算主要通过建立不排水条件下孔隙水压力上升量与排水条件下 通过砂的剪切试验求得的体积变化之间的关系求解R,zl本文结合八木则男提出的孔隙 水压力上升量的实验公式来进行，以便符合不规则的应力条件。在非排水条件下，一次 循环过程中产生的孔隙水压力与排水条件下产生的体积应变的关系可写成 1盘 . 】盒 1命 Qu之一 -丁叮0口 +— — 0刀+— — d八 (2) 功二 IC 用 切 挤.A 式中ka二I‘f-竺~，AE,=,。一N ay二mdv一4 u为孔隙水压力:，‘为体积应 cv N+A 内 N + A Vk 1N+Af 变:。。为砂的体积压缩系数;d为有效应力;A伪振动反复次数;。二r1a,:二*为常数; n和A分别采用5和1,71510 1.2Biot固结方程的有限单元形式 1.单元内任一点的位移及孔隙水压力 (u=防1{u1` 七尸= (N】{P), 式中伽}任切IPzP3A二Ps1T为单元孔隙水压力向量，P为整体坐标系中单元内任一点的 孔隙水压力，P:为整体坐标系中单元第i结点的孔隙水压力口(M为形函数 、二合(1、}‘xl十。,nxl+s,c) 。二1,2.，..,8) 2.应力与孔隙水压力关系的增量形式 仁。卜A{,,)+伽AP,+伽}4P8 (4) 其中A{a}=}a}Ao,△二:AT,,AT,A:二}为总应力增量;AP为残余孔隙水压力增量:APa 为振动引起的振动孔隙水压力增量;Ad为有效应力增量:{M}=(111000]T. 3.结点等价体积和三维结点等价流量 图1是一个微元体dv=dxdydz，定义i结点的等价体积为dv二dxdydzN，见