1. 引言
两相饱和土介质的动力特性远比单项固体介质复杂,因此对土体动力特性的研究往往都视其为简单的单相介质。显然,采用两相饱和土介质模型,即有效应力方法比采用单相固体介质来研究饱和土体的动力响应更为合理,更具有实际意义 [1] [2] [3] 。与总应力方法相比,有效应力方法不但提高了计算精度,更加合理地考虑了动力作用过程中土动力性质的变化,而且还可以预测动力作用过程中孔隙水压力的变化过程、土体液化的可能性及土层软化对场地自振周期及地面振动反应的影响。
本文利用流固耦合单元,分别采用普通地震波和长周期地震波,对某典型深覆盖饱和软土场地进行地震响应分析,研究了基于饱和介质理论的深软场地的地震瞬态响应问题,得出了一些有意义的结论。
2. 分析方法
Biot固结理论于1941年首次提出,是可以准确反映孔隙压力消散与土骨架变形之间耦合作用的真三维固结理论 [4] 。Biot固结方程可表示为:
(1)
式中,
为向量或二阶张量散度算子,
为有效应力张量,
为Biot系数,
为孔隙压力,
为二阶恒等张量,
为体力,
为弹性体积应变,
为Biot模量,
为流量矢量,
为流源。流量与孔压的关系符合达西定律:
,其中
为二阶渗透系数张量,
为梯度算子。
3. 场地土层模型与地震波选取
3.1. 场地土层
选取某典型深覆盖土层场地 [5] (表1),利用流固耦合单元建立有限元模型 [6] ,研究饱和两相介质在地震作用下的响应结果。饱和土体的抗液化强度和振动孔隙水压力与振前应力状态有关,应首先计算初始应力场,得到各个单元的初始应力水平。地震动时程作用持续时间很短,可认为地表位置的孔隙水来不及排出,在分析中采用不排水边界。
3.2. 地震波选取
地震动的选取对计算结果的影响不容忽视,深覆盖软土场地通常对长周期地震动亦有放大作用。本
![](Images/Table_Tmp.jpg)
Table 1. Material parameters of deep soil layer site
表1. 典型深覆盖场地土层的力学参数
文从太平洋地震工程研究中心和日本地球科学与防灾研究中心的强震记录库中选取5条不同频率成分的地震波(见表2),并将峰值统一调整为100 gal,其中普通地震动时程记录包括El Centro波和Taft波,长周期分量丰富的地震动时程纪录包括Loma Prieta波、Chi-Chi波和Tokachi-oki (十胜冲)波。图1中分别列出了5条地震波的加速度反应谱,从图中可以看出,普通地震波反应谱谱值主要集中在0~2 s之间,2 s后谱值迅速下降;长周期地震波反应谱谱值分布则更为宽泛,2 s后长周期部分依然具有较大的谱值,Tokachi-oki波的谱值长周期特征更为明显,周期到达4 s后谱值仍然保持在0.5左右。
4. 深覆盖饱和土层场地地震响应分析
图2中列出了单相介质自由场和饱和两相介质自由场在不同地震动作用下场地沿深度方向的加速度幅值。可以看出,两种介质得到的曲线规律基本一致,但加速度幅值结果不尽相同。这可能是由于振动过程中孔隙水压力变化对计算体系刚度的影响造成的。为明显起见,将两种分析方法得到的地表水平向加速度反应幅值和位移反应幅值列入表3中。可以看出,饱和两相介质自由场的地表响应幅值均大于单相介质自由场。相差最大的是Taft波作用下的结果,加速度幅值和位移幅值分别达到了单相介质场地的1.21倍和1.09倍。但无论是基于单相介质还是饱和两相介质,长周期地震动作用下的加速度和位移反应均大于普通地震波,而且长周期对于地表位移反应的放大效应也强于加速度反应,深覆盖土层的滤波效应明显。
流固耦合单元随着地震作用产生的孔压比是影响场地地震响应结果的重要指标,图3中列出了场地不同土层的孔压比结果曲线。限于篇幅,图中仅列出了第3层淤泥质黏土、第4层黏性土和第6层粉细砂中点位置处的孔压时程曲线。由图3可知,埋深越浅的位置孔压比越大,由于场地土层有效应力随着观测点位置埋深增大而增大,而土层的地震响应又随着埋深的增大而减小,这与实际的地震震害调查经验相一致。
![](Images/Table_Tmp.jpg)
Table 2. Seismic wave records selected
表2. 文中所选的典型地震动记录
![](//html.hanspub.org/file/11-2750396x27_hanspub.png)
Figure 1. Acceleration response spectrum of seismic wave records selected
图1. 选取地震动的加速度反应谱
地震动持时也对孔压比有明显的影响,Tokachi-oki波由于持时是其它地震波的两倍(80 s),流固耦合单元的孔压比也相比其它地震动时程有明显的积累。孔压比增长的速率,与地震波能量在整个时间历程内的分布也有直接关系,幅值出现得晚,孔压比增加得也较缓。
将输入的地震动幅值统一调整为220 gal,得到的场地土层的孔压比时程曲线列于图4中。可以看出随着地震输入能量的增大,孔压比也明显增加。Tokachi-oki地震时程结束后,淤泥质黏土层的孔压比达到了0.56。
5. 结论
1) 饱和土体场地地震反应分析时,基于流固耦合单元的两相饱和土介质模型可以考虑孔隙水的影响,比采用单相固体介质来研究饱和土体的动力反应更为合理,更具有实际意义;
2) 输入地震动时程的频率特性对结果有明显影响,长周期地震动作用下的场地地震加速度和位移反应均大于普通地震波;
3) 孔压比时程曲线受土层埋深与地震动强度、持时和能量分布的影响。埋深较浅的土层孔压比较大,在强震作用下的饱和土场地地震响应分析时,应充分考虑地震动特性综合分析;
![](Images/Table_Tmp.jpg)
Table 3. Amplitude of seismic responses of site surface
表3. 二维场地地震反应地表结果幅值
4) 由于采用了不可压缩的流固耦合单元,饱和土体场地地震反应过程中无法考虑土体震陷和地震动峰值过后孔隙水的消散过程,需要在后续工作中进一步研究。