Study of Simplified Calculation Method for Seismic Settlement of Soft Soil Based on Triaxial Test
-
摘要: 本文通过对天津滨海新区典型软土的静、动力学试验,分析了软土残余应变的变化规律,提出了结合地区震陷经验系数的软土震陷简化计算方法,并利用1976年唐山7.8级地震中天津塘沽地区软土震陷观测资料对计算方法进行了验证。结果表明,本方法用于估算自由地表震陷量是可行的,对中小工程中软土震陷量的估算具有一定的参考价值。Abstract: Based on the static and dynamic mechanical test of soft soil in Tianjin, we analyzed change law of residual strain of soft soil. The estimation method based on static algorithm and empirical coefficient is developed to compute the seismic settlement of soft soil. By comparing the calculated results with the measured data in the 1976 Tangshan earthquake, we proved that this method is feasible. It can be used to estimate the seismic settlement of soft soil in small and middle scale engineering project.
-
Key words:
- Soft soil /
- Dynamic triaxial test /
- Residual stain /
- Calculation method
-
引言
利用夯锤坠落产生的冲击力对软弱场地进行加固,能够有效降低场地压缩性、提高其承载力(张庆国,2003)。强夯施工过程中,夯锤落地瞬间产生脉冲振动,并以波的形式从锤击点呈放射状向外传播。夯击振动影响人体健康,当振动强度达到一定量级时,会对结构安全产生不利影响。因此,厘清强夯振动衰减规律和影响范围,有助于评价施工作业区建筑物结构安全性和强夯法的应用。
国内外针对强夯法的应用研究主要侧重于场地加固机理,针对振动衰减规律的研究相对较少。同时,现有的研究成果中,相当一部分是针对特定工程进行的现场振动监测,其研究的通用性受到限制,很难形成成熟完善的衰减理论和计算方法(程祖峰等,2007;安惠泽等,2010)。
本文依托国家电网安徽换流站场地测试项目,对项目所在区域的粉土场地进行了强夯振动规律现场测试,重点关注不同能级夯锤作业时地面振动的峰值与频谱信息,研究其振动特性与振动衰减规律,建立粉土场地强夯振动衰减计算公式,并应用振动监测数据开展结构安全性研究,确定了强夯施工区建筑最小安全距离。研究成果为黄泛区夯击振动作业区的结构安全性评价和强夯设计计算提供参考。
1. 概况
1.1 试验区工程地质条件
夯击振动测试区位于安徽省宣城市,工程场地为构造剥蚀丘陵,测试区地表土层主要为冲洪积相砂卵石层和坡残积相粘性土,基岩为风化程度不等的泥质粉砂岩,主要土层的物理、力学指标见表 1。宣城地区地震基本烈度为Ⅵ度,根据《建筑抗震设计规范》(GB 50011—2010)(中华人民共和国住房和城乡建设部等,2010)规定,场地所在地层形成年代为第四纪晚更新世(Q3)及以前,因此可不考虑场地液化影响。
表 1 岩土层的物理、力学指标Table 1. The physical and mechanical properties of soil layers岩土层名称 埋深/m 天然含水量W/% 天然重度γ/kN·m-3 液限WL/% 塑性指数IP/% 天然孔隙比e 直剪 标贯试验数N 承载力特征值fak/kPa 凝聚力C/kPa 摩擦角φ/° 粉质粘土 0—4.1 23.8 19.0 36.9 16.9 0.75 45 17.5 14 230 粉质粘土夹卵石 4.1—10.2 26.0 20.5 34.7 16.3 0.61 45 14.1 15 235 黏土 10.2—27.5 27.8 18.8 43.3 19.6 0.85 41 13.4 9 200 泥质粉砂岩 27.5—未见底 22.3 1.2 强夯施工参数
换流站强夯工作区采用能级为6000kN·m和2000kN·m的2种夯锤进行地基处理(图 1)。按照试夯方案,首先采用能级6000kN·m的夯锤(锤重37.5×103kg,落距17.2m)进行点夯,点夯完成后进行场地整平;然后采用能级2000kN·m的夯锤(锤重19.0×103kg,落距11.0m)进行补充点夯,以对地基进行加固。
2. 夯击振动影响测试试验
2.1 振动测试目的
众所周知,受场地阻尼及能量耗散等综合作用的影响,强夯振动幅值随震源距的增大而减小。为确保强夯作业区周围建筑物的安全,选取典型场地进行强夯振动监测,依据监测结果研究强夯振动衰减规律,探讨强夯振动对建筑物安全性的影响及最小安全距离的确定方法。
2.2 测试方案
(1)测点的选取
为研究夯锤作业时地面振动的衰减规律及对周围建筑的影响,参照《地震动力特性测试规范》(国家技术监督局等,1998)进行现场地基振动测试研究。对2000kN·m和6000kN·m能级夯锤分别进行振动监测,探讨强夯振动衰减规律,建立适用于粉土场地的强夯振动衰减计算公式。2个能级强夯振动测试的测点布设方案分别如图 2(a)、(b)所示。2种能级振动监测均设置5个测点,其中2000kN·m能级振动测试的测点距分别为10m、27m、47.6m、70.6m和91.4m,6000kN·m能级振动测试的测点距分别为29.5m、50.7m、68.3m、70.6m和91.4m,每个测点均布设水平及竖向加速度传感器。
(2)传感器及数据采集系统
本次振动测试选用中国地震局工程力学研究所研制的941B型加速度数据采集器及配套信号放大器,配接SigLab数据采集器对夯击作业时产生的振动进行同步监测。941B型传感器直接输出与加速度成比例的电信号,经过放大器调制、滤波后,传至SigLab转化为数字信号并存储在计算机上。经过现场比选,采用1档(加速度档)进行现场测试工作,941B型号1档对应灵敏度为0.3V·s2/m;与941B振动传感器配接的信号放大器通频带为0.025—35Hz,进行数据分析时,将放大器放大倍数与传感器灵敏度相乘,得到总灵敏度。实际记录器采集得到的信号电压最大值除以总灵敏度,得到以m/s2为单位的加速度峰值。
3. 测试结果与分析
3.1 强夯振动衰减规律
对各测点的振动信号进行数据处理,可得到各点加速度时程曲线的峰值Ap,对测试数据进行傅立叶谱分析,可获取强夯振动的频谱信息。以6000kN·m能级测点1、3的测试结果分析强地面运动的一般规律,2测点的加速度时程曲线及傅立叶谱如图 3所示。
振源距29.5m(测点1)处得到的峰值加速度均值为0.388m/s2,振源距68.3m(测点3)处得到的加速度峰值为0.077m/s2,仅为测点1处加速度峰值的19.85%。由图 3(a)可见,2测点振动持续时间均约0.5s;由图 3(b)可见,夯锤作业时夯击产生的冲击波振动频率集中在1—30Hz内,测点1的2个卓越峰值点对应的频率分别为7.3Hz(0.13s)和16.4Hz(0.06s),测点3的2个卓越峰值点对应的频率分别为7.3Hz和13.31Hz,其中f=7.3Hz对应的振动影响最大。后续计算振动速度、位移时,采用该卓越频率进行计算。图 4、5进一步给出了不同能级夯击作业时,地震振动随距离的衰减规律。
6000kN·m能级夯锤作业时,引起的地面振动最大加速度随振源距的衰减规律为:y=1.33e-0.042x,其中x为振源距(m),相关系数R=0.9942。而2000kN·m夯锤引起的地面振动的衰减规律为y=0.04+4.80e-0.131x,相关系数R=0.9977。
由图 4(a)、(b)可知,无论哪个能级的强夯振动,振动衰减均呈现近处快、远处慢的特点,衰减规律曲线近似可用指数函数表达;对比2条衰减曲线(图 5)可以看出,振源距相同的情况下,6000kN·m能级的振动衰减更快,表明振动能量越大,加速度峰值衰减越快。
表 2、3给出了2种能级夯锤作业时,各测点的水平及竖向加速度峰值信息。
表 2 6000kN·m夯锤引起的振动Table 2. The vibration caused by 6000kN·m pounder距夯点距离/m 竖向加速度/m·s-2 水平向加速度/m·s-2 双向平均值/m·s-2 速度/mm·s-1 实测值 均值 实测值 均值 29.5 0.420
0.320
0.470
0.3000.378 0.410
0.370
0.330
0.4800.398 0.388 6.950 50.7 0.180
0.120
0.120
0.120
0.140
0.1000.130 0.180
0.160
0.210
0.1400.173 0.152 2.720 68.3 0.077
0.077
0.0750.076 0.110
0.087
0.0870.095 0.086 1.540 70.6 0.073
0.0680.071 0.100
0.0740.087 0.079 1.420 91.4 0.012
0.014
0.0220.016 0.028
0.037
0.0380.034 0.025 0.450 表 3 2000kN·m夯锤引起的振动Table 3. The vibration caused by 2000kN·m pounder距夯点距离/m 竖向加速度/m·s-2 水平向加速度/m·s-2 双向平均值/m·s-2 速度/mm·s-1 实测值 均值 实测值 均值 10.0 0.980
0.730
0.8080.839 1.580
1.810
2.0901.827 1.333 23.880 27.0 0.121
0.095
0.0710.096 0.209
0.250
0.2900.250 0.173 3.100 47.6 0.040
0.025
0.061
0.045
0.052
0.0590.047 0.048
0.070
0.077
0.084
0.100
0.1150.082 0.065 1.160 70.6 0.035
0.033
0.0370.034 0.034 0.610 91.4 0.004
0.006
0.0060.005 0.007
0.008
0.0070.007 0.006 0.110 李福民等(2002)和罗辉才等(2011)的研究成果表明,夯击产生的振动在其有效影响范围内,振动主频率变化不大,故本次测试并未对振动卓越频率随振源距的变化开展专门分析。由测试数据可知,在振源距较小的测点上,水平向振动加速度峰值明显大于竖向加速度峰值,随着振源距的增大,水平向与竖向加速度峰值差异越来越小。
测试点i的振动速度峰值可由公式Vi=Ai/2πf求得,其中Ai为强夯引起的振动波在测点i的加速度振幅。6000kN·m能级夯锤作业在距作业点30m处,产生的振动速度峰值为6.95mm/s,50m处的速度峰值为2.72mm/s;2000kN·m能级夯锤作业在距作业点10m处,产生的振动速度峰值为23.88mm/s,50m附近处的速度峰值为1.16mm/s。
3.2 最小安全施工距离
(1)确保正常工作和生活的最小安全施工距离
原冶金部《机器动荷载作用下建筑物承重结构的振动计算和隔震设计规程》(冶金工业建筑研究总院,1990)中,规定了工作和生活环境的允许振动速度限值,见表 4。
表 4 工作和生活环境的允许振动速度Table 4. The suggested vibration velocity of working and living environment环境类别 允许振动速度V/mm·s-1 白天(6时至20时) 夜晚(20时至次日6时) Ⅰ 0.25 0.13 Ⅱ 0.50 0.25 Ⅲ 1.00 0.50 Ⅳ 2.80 2.80 注:环境类别Ⅰ类,对环境振动要求特别严格,如医院、学校等;Ⅱ类,对环境振动要求比较严格,如宿舍等生活区;Ⅲ类,允许有轻微的振动感觉,但不影响精神集中,如一般的办公室等公共场所;Ⅳ类,车间范围内的非操作区。 参照表 4规定的允许振动速度对夯锤作业振动影响进行评价。在不考虑夜间施工的前提下,当6000kN·m夯击作业时,距离作业点70m范围内的振动速度峰值超过Ⅰ、Ⅱ、Ⅲ类环境的允许振动速度;当2000kN·m夯锤作业时,距离作业点50m范围内的振动速度峰值超过Ⅰ、Ⅱ、Ⅲ类环境的允许振动速度,不适宜生活或办公。
(2)确保结构安全的最小安全施工距离
我国《建筑抗震设计规范》(以下简称《规范》)确定的基本原则是“小震不坏,中震可修,大震不倒”(中华人民共和国住房和城乡建设部等,2010),具体操作时,作为基本设防烈度的中震在抗震设计、验算中却没有直接体现出来,而是采用设防烈度的常遇地震进行设计或验算,而中震对应的地震动加速度峰值是小震对应加速度峰值的3倍。在确定强夯作业最小安全距离时,若采用设防烈度的中震作为判定标准,强夯作业的最小安全距离将被缩短,使得最小安全距离之外存在峰值大于小震对应加速度峰值的振动。因此,采用中震(设防烈度)作为评判依据不能完全满足结构在线弹性阶段“小震不坏”的设计要求。
《规范》中规定,时程分析时所有地震动加速度时程曲线最大值的取值依据为设防烈度下多遇地震(小震)和罕遇地震(大震)对应的地震动加速度峰值,其中多遇地震对应的地震动加速度峰值用于验算结构弹性极限状态,罕遇地震对应的地震动加速度峰值用于验算结构的塑性极限状态。故采用设防烈度下多遇地震(小震)对应的地震动加速度峰值作为结构最小安全距离评判标准更为准确。
对建筑物进行时程分析时,输入加速度时程的最大值可按表 5选取。通过对比可以发现,当6000kN·m能级夯击作业时,距作业点30m处的振动加速度幅值为38.8cm/s2,相当于对建筑物输入Ⅶ度多遇地震的强震作用,夯击振动持续时间虽短,但重复次数多,对建筑物的破坏作用同样不可轻视,对30m范围内Ⅵ度设防建筑及未设防建筑具有较高的地震危险性,建筑物虽不至倒塌,但局部破坏的风险较高;距作业点50m处的振动加速度幅值降为15.2cm/s2,略低于Ⅵ度多遇地震的强度,据此可判定此振动对建筑物基本无害。因此,6000kN·m能级夯击作业时,确保建筑物结构安全的最小距离应为50m。当2000kN·m能级夯击作业时,距作业点27.2m处振动加速度幅值为17.3cm/s2,略低于Ⅵ度多遇地震的强度,此时的振动可判定为对建筑物基本无害,因此,2000kN·m能级夯击作业时确保建筑物结构安全的最小距离应为30m。
表 5 时程分析所用地震加速度时程曲线的最大值(cm/s2)Table 5. The Maximum values for the seismic acceleration of ground motion used in time-history analysis (cm/s2)地震类型 Ⅵ度 Ⅶ度 Ⅷ度 Ⅸ度 多遇地震 18 35(55) 70(110) 150 罕遇地震 220(310) 400(510) 620 注:括号内的数值分别用于设计基本地震加速度为0.15g和0.30g的地区。 3.3 降低振动影响的措施
按照波在介质内的传播形式,可将强夯产生的冲击波分为体波(P波、S波)和面波(瑞利波、勒夫波)2种。体波率先到达,但携带的能量较少,同时体波的衰减速度快;面波传播速度慢,衰减速度相对较慢,且携带的能量大,对建筑物的影响大。面波中,瑞利波携带的能量占其总波能的2/3,对激振源附近地表附着物的破坏力较大(杨成林,1993;雷学文等,2002;陆伟东等,2002),因此,减小瑞利波能量是降低夯击作业对建筑物影响的有效途径。
根据瑞利波的传播特性,在实际工程中多采用隔震沟来减小振动影响。影响隔震沟减震效果的因素很多,大致可归为隔震沟深度、宽度、截面形状等。孔令伟等(1996)通过理论分析与试验研究,得出瑞利波在垂直向的影响深度约为1.0—1.5倍的波长λR,而夯击作业时λR通常为8—14m;前人通过大量现场试验、数值模型分析得出:隔震沟深度对减震效果影响显著,沟宽度对减震效果影响不大;对于相同宽度、深度的减震沟,当其位于作业点与保护对象的中间位置时,减震效果最佳(Wood,1968;邓亚虹等,2007;余德运等,2011)。在实际工程中,综合考虑具体施工条件与经济因素,可将减震沟深度设置为4m左右。
4. 结论
通过对粉质土场地夯击作业开展现场振动测试,研究了2种能级夯击作业的振动衰减规律及对建筑物的影响,得出如下结论:
(1)对于粉土场地,强夯振动衰减规律可由形如y=Aebx的指数函数近似表达。本文测试采用的2种能级强夯作业产生的振动衰减均呈现进处快、远处慢的规律。
(2)夯击振动为典型冲击荷载,振动持续时间多在0.5s内,其振动频率范围通常在20Hz内,同时,振源距对振动卓越频率的影响不显著。对于结构而言,10Hz以下的低频振动对建筑物影响最大。
(3)夯击能越大,强夯振动衰减速度越快。与2000KN·m夯锤比较,振源距大致相当时,6000KN·m夯锤产生的振动衰减速度更快。
(4)距离夯击作业点超过50m时,2个能级夯锤作业产生的振动加速度峰值均在18cm/s2以下,低于Ⅵ度设防条件结构时程分析对应的地震动加速度限值,不会对建筑物造成破坏,可将建筑物最小安全距离定为50m。
-
表 1 天津软土的邓肯参数
Table 1. Duncan parameters of soft soil in Tianjin
土类 Ks/kPa Ns 摩擦角/° 粘聚力/kPa 破坏比 淤泥质土 1237 0.465 12 20 0.478 淤泥质粉质粘土 2930 0.39 12 44 0.687 表 2 50kPa动应力下各组拟合成果
Table 2. Fitting result of residual strain under 50kPa
试样深度/m 拟合参数 相关系数 A b 6.5 0.0014 1.72 0.9941 7.0 0.00018 1.95 0.987 7.5 0.000035 2.47 0.997 9.5 0.00010 2.53 0.9987 11.5 0.00010 2.31 0.995 表 3 60kPa动应力下各组拟合成果
Table 3. Fitting result of residual strain under 60kPa
试样深度/m 拟合参数 相关系数 A b 6.0 0.0040 2.28 0.9975 7.0 0.00040 1.96 0.9949 7.5 0.0040 1.77 0.9985 11.5 0.00010 2.52 0.995 12.5 0.00060 2.142 0.9988 13.5 0.0020 2.14 0.997 表 4 不同震级下等效振动次数
Table 4. Equivalent vibration frequency under different magnitude
震级 6.5 7 7.5 8 振次 8 12 20 30 表 5 天津地区软土残余应变
Table 5. Residual strain of soft soil in Tianjin
序号 深度 层厚 土类 7.5级地震影响下残余应变 8级地震影响下残余应变 1 4.7 4.7 填土 - - 2 5.5 0.8 填土 - - 3 8.3 2.8 淤粉 0.01 0.2 4 9 0.7 粉粘 0.1 0.3 5 12 3 淤粘 0.1 0.3 6 15.8 3.8 淤粘 0.2 0.35 7 17.5 1.7 粉粘 - - 8 20 2.5 粉粘 - - 表 6 土层震前模量
Table 6. Soil modulus before earthquake
序号 深度/m 层厚/m 土类 大主应力/kPa 小主应力/kPa 震前模量/kPa 1 4.7 4.7 杂填土 42 28 9571 2 5.5 0.8 素填土 50 34 18103 3 8.3 2.8 淤粉粘 76 50 14632 4 9 0.7 粉粘 81 54 13723 5 12 3 淤粘 103 69 14428 6 15.8 3.8 淤粘 131 88 15186 7 17.5 1.7 粉粘 147 98 16207 8 20 2.5 粉粘 169 113 14232 表 7 土层震后模量
Table 7. Soil modulus after earthquake
序号 深度/m 层厚/m 土类 震前模量/kPa 7.5级地震影响下割线模量/kPa 8级地震影响下割线模量/kPa 7.5级地震影响下软化模量/kPa 8级地震影响下软化模量/kPa 1 4.7 4.7 杂填土 9571 - 9571 9571 9571 2 5.5 0.8 素填土 18103 - 18103 18103 18103 3 8.3 2.8 淤粉粘 14632 251667 14632 13828 6765 4 9 0.7 粉粘 13723 27033 13723 9103 5439 5 12 3 淤粘 14428 34433 14428 10168 6392 6 15.8 3.8 淤粘 15186 21903 15186 8968 6861 7 17.5 1.7 粉粘 16207 - 16207 16207 16207 8 20 2.5 粉粘 14232 - 14232 14232 14232 表 8 地震影响下自由地表震陷
Table 8. Seismic settlement of free ground under earthquake motion
序号 深度/m 层厚/m 土类 7.5级地震影响下各层震陷量/m 7.5级地震影响下总沉震陷量/cm 8级地震影响下各层震陷量/m 8级地震影响下总震陷量/m 1 4.7 4.7 填土 - 4.9 - 12.9 2 5.5 0.8 填土 - - 3 8.3 2.8 淤粉 0.001 0.025 4 9 0.7 粉粘 0.003 0.009 5 12 3 淤粘 0.013 0.039 6 15.8 3.8 淤粘 0.032 0.056 7 17.5 1.7 粉粘 - - 8 20 2.5 粉粘 - - 表 9 不考虑模量软化的自由场地沉降量
Table 9. Free ground settlement without consideration of soften model
序号 深度/m 层厚/m 土类 竖向应力/kPa 动应力/kPa 震前模量/kPa 应变值/% 沉降值/cm 1 4.7 4.7 杂填土 42.3 20.6 9571 0.007 20 2 5.5 0.8 素填土 50.3 24.2 18103 0.004 3 8.3 2.8 淤粉粘 75.5 34.9 14632 0.008 4 9 0.7 粉粘 81.1 37.1 13723 0.009 5 12 3 淤粘 103.3 45.1 14428 0.010 6 15.8 3.8 淤粘 131.42 54.0 15186 0.012 7 17.5 1.7 粉粘 146.72 58.5 16207 0.013 8 20 2.5 粉粘 169.22 64.6 14232 0.016 表 10 天津滨海地区自由场地软土震陷量修正系数
Table 10. The settlement correction coefficient of free soft soil site in Tianjin
地面峰值加速度/g 震级/M 修正系数C 0.25 7.5 0.2 8 0.6 0.4 7.5 0.25 8 0.5 -
陈国兴, 李方明, 从卫民, 2004.多层建筑物地基震陷的简化计算方法及其影响因素分析.防灾减灾工程学报, 24(1):47-52. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXK200401008.htm 陈建峰, 俞松波, 叶铁锋等, 2008.软土地基加筋石灰土路堤离心模型试验研究.岩石力学与工程学报, 27(2):287-293. http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200802012.htm 刘恢先, 1986.唐山大地震震害-第二册.北京:地震出版社. 江席苗, 2009. 汶川地震地基基础震害调查研究. 上海: 同济大学. 李冬, 陈培雄, 吕小飞等, 2011.软土地震震陷研究现状综述.工程抗震与加固改造, 33(2):130-135. http://youxian.cnki.com.cn/yxdetail.aspx?filename=SCSA201704303&dbname=CJFDPREP 孟上九, 袁晓铭, 孙锐, 2002.建筑物不均匀震陷机理的振动台实验研究.岩土工程学报, 24(6):747-751. http://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200206016.htm 吴跃东, 2008. 软土震陷计算方法及影响研究. 南京: 河海大学. 王建华, 要明伦, 1993.饱和软粘土振动弱化特性的研究.水利学报, 12:37-43. doi: 10.3321/j.issn:0559-9350.1993.06.006 谢君斐, 石兆吉, 郁寿松等, 1988.液化危害性分析.地震工程与工程振动, 8(1):61-77. http://www.cnki.com.cn/Article/CJFDTOTAL-DGGC198404006.htm 谢君斐, 石兆吉, 1981.神头电厂地基震陷的初步分析.中国科学院工程力学研究所. 杨石红, 刘静蓉, 刘金珠等, 1997.软弱地基土层震陷简化计算方法研究.世界地震工程, 13(2):53-61. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC702.007.htm 郁寿松, 石兆吉, 1989.土壤震陷试验研究.岩土工程学报, 11(4):35-44. http://www.cnki.com.cn/Article/CJFDTOTAL-YTGC198904004.htm 中国科学院工程力学研究所, 1979.海城地震震害.北京:地震出版社. 张建毅, 张宇东, 徐国栋, 王强, 2012.软土震陷判别方法分析及工程应用.世界地震工程, 28(3):53-59. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201203012.htm 周健, 蔡宏英, 许朝阳, 2000.软粘土地基震陷分析.工程抗震, (1):40-42, 39. http://www.cnki.com.cn/Article/CJFDTOTAL-GCKZ200001011.htm Lee K. L., 1974. Seismic Permanent Deformation in Earth Dams. School of Engineering and Applied Science, University of California, Los Angeles, California, December. Monismith C. L., Ogawa N., Freeme C. R., 1975. Permanent deformation characteristics of subgrade soils due to repeated loading. Transp. Res. Rec. No.537. Transportation Research Board. Washington D.C.:1-17. https://trid.trb.org/view.aspx?id=511499 Narikawa M., Ohshima Y., Nishimura I., et al., 1999. Settlement and seismic behavior of florting foundation on thick cohesive strata. In:Seco e Pointo, ed. Earthquake Geotechnical Engineering. Rotterdam:A A Balkema.421-426. Seed H. B. and Clarence K. C., 1966. Clay Strongth under Earthquake loading Condition, Journal of the Soil Mechanics and Foundation Division, ASCE, 1992, SMZ, March. -