水泥搅拦桩加固地基地震放大效应的有限元分析.pdfVIP

水泥搅拌桩加固地基地震放大效应的有限元分析 郭海燕 孟凡顺(青岛海洋大学) 吴世明(浙江大学) I概述 油罐由于运输关系，大部分建在沿海地区的软土地基上，软土地基上构筑物在最近几 次大地震中曾发生过大规模的破坏．美国Loma Prieta大地震便是典型的实例。这次地震引 起大量建在该海湾软土地基上的结构物(包括油罐)和交通设施的破坏。现场测试和理论 分析表明，场地地震放大效应和地基液化是结构物破坏的主要原因。 南京炼油厂石阜桥油罐区位于西霞山北麓长江岸滩地带，场地上部为淤泥质粘土，下 部为细砂，属于软弱场地，其工程地基状况和地基环境同美国Lomal：b-ieta地震情况比较接 近。该罐区5万立方米油罐扩建工程的地基加固采用了浙江大学岩土所提出的单轴深层水 泥搅拌桩和加筋垫层联合处理油罐软基方案。对于该软基处理方案．抗震加固效果如何必 须作出合理地评估。由于实验条件所限，对场地地震放大效应及抗震加固效果的评估，目 前大多采用理论计算方法，文献【t1采用剪切粱法，将士壤材料动力非线性性质等效线性化． 化成一个常微分方程组，然后用Wilson．0法求解，评价该地基的抗震加固效果。等效线性 化法计算简便，工程应用较广，但该法本身有较大程度地近似，计算结果有一定的误差． 本文将采用有限元法，按照土体的非线性应力一应交模式，即将Hadin的双曲线模式 作为骨架曲线，加卸荷载采用Pyke等修正的masin窖二倍法则．选用粘性边界作为侧面人 工边界，编制非线性地基动力有限元计算程序．对南京炼油厂5万立方米油罐地基进行地 震反应计算分析．研究其抗震加固效果。 2非线性地基动力有限元计算方法 2．I有限元计算公式 根据Hamilton变分原理，可导出动力有限元基本方程四 (1) 【M】{“}-卜【cl{u}H足】{“)={，(f)) 式中：fM】为总体质量矩阵；【q=fc】‘+【q功总体阻尼矩阵，【c】‘为瑞霄阻尼矩阵，【q2 为粘性边界阻尼矩阵：【K】为总体刚度矩阵：{u)为位移矢量：{f(O)为等效结点荷载向量· 这里假定土体在地震荷载作用下，正应力与正应变之间关系是线弹性的。动态剪应力剪应 变关系采用非线性模式。 一565 2．2土体在循环剪切荷载作用下的本构关系 在等幅循环剪切荷载作用下．土体应力应变的曲线轨迹是一系列滞回圈，滞回圈的顶点 轨迹骨架曲线可用Hardin和Drnevichl3]提出的双曲线表示： T：Goy、／(1+y／Y。) (2) 式中Go为初始切线模量；Y，=T．I／Go为参考应变：T。。为土的抗剪强度。 考虑土体的非线性，关键是确定卸载与再加载时应力一应变滞回曲线的形状从而确定各 时刻的切线剪切模量。目前广泛采用的是iasing二倍法，即把骨架曲线的尺度放大一倍作 为卸载和再加载线段。但Masing二倍法得出较大的滞回圈，其阻尼比与实测阻尼比不符，为 此，本文采用PyhI．I修正的iasing曲线如下： (3) f_々掣IL ‘1咐一儿)／nr， 一 ^=|±l—fc，f。I (4) 式中(T。、Y。)为转折点坐标。正、负号按加、卸荷分别取用。 在非等幅循环荷载(地震荷载)下尚应满足加、卸荷准则。即；如果}r9I。f。，卸荷 或再加荷曲线与当前骨架曲线相遇时，应力点将沿当前骨架曲线前进；如果If9l心，。则 把骨架曲线放大2倍作为加载、卸载曲线；如果卸载再加载曲线超出与最大应力比点对应 的修正后iamg曲线的加载分枝上，则应力点就转到该曲线的分枝上． 2．3动力方程组的求解及程序验证 本文采用预校正形式的Ncwmark方法啪，对动力方程组进行逐步积分，求解离散系统 的动力响应，在微机上编写了二维地基动 X 力育限元计算程序．