软夹层对地表地震动参数的影响.docx

1、软夹层对地表地震动参数的影响防 灾 科 技 学 院毕 业 论 文题目软夹层对地表地震动参数的影响学生姓名李凯 学号 115021132 系 别防灾工程系专 业勘查技术与工程班级115021132 开题时间2015年3月10日答辩时间2015年5月30日指导教师张宇东 职 称讲师 软夹层对地表地震动参数的影响作 者：李凯指导教师：张宇东摘要 本文以类类两种工程场地为原型，分别对软夹层的厚度和埋深对软土场地地震效应的影响进行了数值分析。利用程序对构造剖面、不同的输入地震波时程及不同的峰值加速度水平，进行了场地地震反应分析，分析软夹层对地震动峰值加速度的影响、峰值速度，反应谱最大值以及反应谱最大值所

2、对应的周期的影响、特征周期的影响等。关键词：软夹层；地表地震动；土层计算模型；场地地震反应；地震动输入；Motion Parameters soft interlayer on surface groundCandidate：Li kaiSupervisor：Zhang yudongAbstract： In this paper, two project site for the prototype, thickness and depth of the impact of soft soil site seismic effect of soft sandwich on numerical

3、analysis. Use program for different configurations profiles of different process and a different level of the input peak acceleration of seismic waves, seismic response analysis conducted to explore the soft interlayer on the impact of peak ground acceleration response spectrum maximum, and the resp

4、onse spectrum corresponding to the maximum the impact of the cycle, the impact of other features cycle.Keywords：soft interlayer；the surface ground motion; soil model; seismic response; seismic input. 目录1、引言 11.1 研究状况以及研究意义 11.2 存在的问题 11.3 本文研究内容 12、一维土层地震反应分析方法 22.1 一维场地力学模型简介 22.2一维土层地震反应线性化计算程序（RS

5、LEIBM）简介 53、数据的选取 63.1 地震动输入 63.2 场地钻孔资料 84、模型计算 104.1 软夹层厚度关系计算 模型1-3 104.2 软夹层埋深关系计算 模型4-6 125、计算结果与分析 135.1 类场地计算结果及分析 145.2 类场地计算结果及分析 186、结论 23致谢 24参考文献 251、引言1.1 研究状况以及研究意义我国作为一个地震多发国家，地震灾害很容易导致大量的人员伤亡，并且严重破坏人类的生存环境。然而现阶段关于短期地震预报仍旧没有突破性的进步，所以人们的重视更多的转向土层结构对地震地表运动的影响，并在此方面进行了诸多研究。深厚软弱的沉积层多是江河下游

6、地区或沿海地区，由于形成年代与成因类型的差异，各土层的动力特性也表现出明显的不均匀性，层状构造较为显著，且软夹层的存在现象也较为普遍。软夹层所处土层位置和厚度的不同都可能对场地地表地震动产生影响。在与场地条件有关的典型震害中可以发现，场地中与软土相关的结构对场地震动有较为直接的影响。软土是在一定的地质条件下，具有成因特性和局域性分布特征的一类特种土。它具有含水量高、孔隙比大、强度低等特点，具有流变性和结构性等特征。1.2 存在的问题工程中常见的土体一般为较长的地质历史时期期间，经由复杂的地质作用从而形成的以层状土层排列的土体。因其地质动力作用及其环境成因不同，各个土层的动力特性有明显的差异性。

7、大量的地震观测结果表明，土体中软土层的存在对场地反应有着显著的影响。软土层对地震动的影响取决于软夹层本身的厚度、软夹层的埋置深度等。但是，在我国及其他国家的结构抗震设计规范中，一般都未考虑土层结构对场地反应的影响效果。因此，在建筑抗震设计中，如遇到存在此类含有软夹层的场地时，设计人员按照现行建筑抗震设计规范进行的设计，并不能正确的反映真实场地反应对建筑的影响作用。1.3 本文研究内容本文的主要工作即是研究土层结构中软夹层厚度及埋深对地表地震动反应的影响作用。选取天然地震动El波改造成的50gal，100gal，200gal时的地震动为基岩输入地震动，分别选取、类场地作为研究软夹层厚度和埋深的土

8、层反应模型，将不同厚度或埋深位置的软土层替代原有土层，组合出一定数量的工况模型。在上述条件下，应用程序分别计算上述各种土层分布条件下，在不同的输入地震动下的场地条件。再对所得的大量计算结果进行数据处理及分析，分别观察场地地震反应随土层的厚度及埋深的改变而变化的情况，进而得出软夹层对场地地震反应的影响规律。2、一维土层地震反应分析方法2.1 一维场地力学模型简介土体的地震反应分析即为对地震动作用下场地中各点的位移、速度、加速度及应力、应变等反应量的求解，是场地设计地震动参数的确定和抗震设计中的重要组成部分。分析过程中的两个关键点是建立场地力学模型和选择合适的动力学参数。 入射剪切波 图2.1 土

9、层地震反应分析的一维力学模型 其中，场地力学模型分为三种，分别是一维场地模型、二维场地模型和三维场地模型。如图2.1所示，一维场地模型假定土层沿两个水平方向均匀无变化，仅沿竖向有变化且层内无变化。当一些局部范围内场地条件较为均匀简单时，可以将场地介质模型简化为成层土层模型，此模型即属于一维场地模型。一维场地地震反应分析方法虽然有其局限性，但是对于局部范围内地面、土层界面及基岩界面比较平坦的场地来说，以此为模型进行土层的地震反应分析已能满足工程要求。二维场地模型适用于局部范围内，地面、土层界面、基岩界面只在一个水平方向上平坦的场地。而三维场地模型适用于在局部范围内，地面、土层界面、基岩界面在两个

10、方向均起伏较大的场地。 在一维成层场地地震反应的分析方法中，在考虑了土体的非线性效应厚，出现了等效线性化分析方法及直接时域非线性逐步积分方法。等效线性化土层地震反应分析方法是一种间接考虑土体非线性特性的方法（Idriss etal.,1968;廖振鹏等，1989），它是在频域线性波动分析方法的基础上利用非线性动力方程的等效线性化处理手段给出的。这一方法可以分成两部分，一是线性方程的频域波动求解，二是土体非线性的等效线性化处理。这一方法尤其应用的范围、条件及其不足之处。本质上讲，等效线性化方法仍属于线性计算方法，它只能从实践过程的平均意义上粗略地体现土体的非线性影响，而不能反应土体的非线性物理过

11、程。另外，等效线性化方法还存在一个适用的条件，就是非线性程度较小。为了研究土体非线性的物理过程以及土体反应将进入大非线性范围等情况，直接非线性土层地震反应分析方法被提出并得到发展，这种方法的两个关键点即为土体非线性应力-应变关系的给出及非线性动力方程的数值积分方法的选取。本文将采用一维场地模型，利用等效线性化土层地震反应分析方法对选取的两个典型剖面软夹层厚度变化时的地震反应进行分析，得到不同工况下的Amax、Vmax、Dmax、Tg及Samax，根据对以上数值的比较，得到软夹层厚度的变化对土层地震反应结果的影响。2.1.1 线弹性土层的稳态地震反应 土层模型如图2.1所示。个土层覆盖在基岩均匀

12、半无限空间之上，各覆盖层厚度、介质质量密度和剪切模量分别为，和， =1，2，下卧基岩半空间的质量密度和剪切模量为和。各层界面的编号已标示于图7.1.1中。采用局部坐标系并将z轴的坐标原点设置在各层上界面，正方向垂直向下。设第层(基岩半空间)的入射地震波位移(略去时间因子，下同)为 (2.1.1)其中，和分别为基岩的剪切波波数和波速。第层地震波位移的频域一般解可以表示成 (2.1.2)其中，分别为第层剪切波波数与波速，和分别为第层介质内上行和下行波波幅系数。求解这一波动问题的关键是传递矩阵概念。所谓传递矩阵就是把任意两层的波幅系数联系起来的矩阵，它可由相邻两层的波幅系数之间的转换关系导出。建立了

13、传递矩阵之后，依据边界条件即可求得问题的解答。2.1.2 线性滞回阻尼土层的稳态地震反应当考虑线性滞回阻尼效应时，稳态地震反应可由线弹性土层的稳态解通过简单替换求得。为此，仅需将第层介质的剪切模量换为 (2.1.3)式中为第层土介质的阻尼比。由于波速，波数，亦可通过将或分别换为或获得， (2.1.4)2.1.3 非线性土层暂态地震反应的等效线性化解法 由于在非线性条件下叠加原理不成立，土层的非线性暂态地震反应不能利用傅里叶变换通过叠加各个频率的稳态解求得。为了应用叠加原理必须引入新的假定，这就是等效线性化的假定。所谓等效线性化就是在总体动力学效应大致相当的意义上用一个等效的剪切模量及等效阻尼比

14、去替换所有不同应变振幅下的及。由于及与应变振幅无关，整个问题化为线性问题。因此，实际计算时，先假定每一土层内介质反应的初始等效动力剪切应变，利用上述方法进行反应计算，并计算出相应的土层内中点处介质的剪应变反应的最大值，而后取每一土层内层中点处介质反应的最大剪应变值乘以折减系数(常取0.65)的值作为该土层中介质的等效剪应变的计算值。比较计算所用等效剪切应变及计算所得等效剪切应变相对应的等效动力剪切模量和滞回阻尼比值，如果它们的相对误差都小于给定的允许误差(这里取0.05)，则认为土体的非线性特性的考虑满足了要求；否则，以最新计算所得等效剪切应变值取代初始等效剪切应变值，并重复上述计算过程，直到

15、相对误差都小于允许误差为止。2.2一维土层地震反应线性化计算程序（RSLEIBM）简介本程序含七个输入输出文件。四个输入数据文件：DATA.DAT，DATA1.DAT，DATA2.DAT，DATA3.DAT。输入数据文件格式均为自由格式。三个计算结果输出文件：RTEL1.DAT，RTEL2.DAT，RTEL3.DAT。其中DATA.DAT文件为计算输入加速度时程离散值文件，各条加速度时程按顺序排列。DATA1.DAT文件中的数据为与计算剖面无关的计算控制参数值。 DATA2.DAT文件为剖面参数及输入输出控制参数数据文件，根据我们的土层模型输入相关数据。DATA3.DAT文件为土的非线性曲线（

16、），G（）的数据文件。RTEL1.DAT文件为计算输入数据及计算结果保留文件。RTEL2.DAT文件为加速度时程、剪应变时程及相应的传递函数计算结果文件。RTEL3.DAT此文件为加速度反应谱、峰值加速度（速度、位移）计算结果文件，根据此结果进行对比分析。3、数据的选取3.1 地震动输入3.1.1 地震动输入特性 本文选用El波作为输入地震动，进行场地地震效应分析，用以反应改变软土层厚度对地震反应的影响。将El波改造成50gal，100gal，200gal，利用soil程序对两个分别属于、类场地的钻孔资料进行不同层剖面、不同输入地震动进行场地地震反应分析。得到的地震输入时程图 如图3.1所示：

17、（a）El波50Gal改造时程图（b）El波100Gal改造时程图（c）El波200Gal改造时程图图3.1 地震记录加速度时程曲线本文将在不同输入加速度水平下，研究地震动强度对软弱土层计算结果的影响，经改造后的地震波峰值如表3.1所示。表3.1 改造地震波峰值地震波加速度峰值（gal）El波501002003.1.2 场地类别划分标准场地条件是影响地震动特征和结构震害的重要因素。反应谱曲线的特征参数，在抗震设计规范中是通过场地分类来规定的，并以此粗略地估计不同场地条件对设计地震动的影响。目前国内外的场地分类多以土层剪切波速和覆盖层厚度这两个指标为主。建筑抗震设计规范（GB50011-2010

18、）是根据场地土类型和覆盖层厚度，将场地划分为、 、共四种类型，场地类别划分标准如表3.2。表3.2 场地类别划分与等效剪切波速（m/s）及覆盖层厚度（m）的关系等效剪切波速（ms）场 地 类 别01Vs8000800Vs5000500Vse25055250Vse150335050Vse1503315158080注：表中VS系岩石的剪切波。经计算，本文选取的两个场地分属于和类场地。3.2 场地钻孔资料本文选取了有代表性的两个软弱场地：场地一为西昌市某场地，依据各分层土剪切波速计算得该场地等效剪切波速为151 ms，且厚度为33 m，由我国建筑抗震设计规范(GB 5001l-2010)场地分类标准

19、确定该场地为类建筑场地。场地二为天津滨海某场地，依据各分层土剪切波速计算得该场地等效剪切波速为170 ms，且厚度为100 m，由我国建筑抗震设计规范(GB 5001l-2010)场地分类标准确定该场地为类建筑场地。3.2.1 类场地的钻孔资料表3.3 类场地原始土层剖面钻孔资料序号土层名称土层深度（m）土层厚度（m）波速(m/s)质量密度cm/g3土类1填土4.14.11301.982粉土6.42.31402.023粉质粘土9.93.52302.034 4粉砂10.80.92602.055粉质粘土16.25.42602.026粉土18.52.32002.027细砂22.03.52302.16

20、8砾石23.21.23802.279卵石30.16.94402.2410基岩5002.393.2.2 类场地钻孔资料表3.4 类场地原始土层剖面钻孔资料序号土层名称土层深度（m）土层厚度（m）波速(m/s)质量密度cm/g3土类1人工填土5.35.31221.982粉土7.32.01282.023粉土15.38.01592.024粉质黏土17.62.32452.035粉土19.31.72642.026粉质黏土22.12.82702.037粉质黏土26.54.43012.038粉砂31.55.03032.059粉质黏土342.53222.0310粉砂38.14.13672.0511粉砂42.24

21、.13672.0512粉砂46.24.03672.0513粉质黏土481.83292.0314粉砂53.85.83202.0515粉质黏土595.24092.0316粉质黏土64.25.24092.0317粉质黏土694.84702.0318粉砂72.43.44362.0519粉砂75.83.44362.0520粉质黏土815.24422.0321粉质黏土86.25.24422.0322粉质黏土91.35.14422.0323粉砂1008.75042.0524基岩2.393.2.3 各类土的动剪切模量比和阻尼比序号模量比阻尼比剪应变0.050.10.51510501001淤泥质黏土G/Gmax

22、0.9850.9700.8450.7300.3200.2100.0850.0580.0120.0150.0330.0550.1360.1700.2000.2052粉土G/Gmax0.9600.9300.7700.6500.3000.2000.0600.0350.0120.0170.0360.0500.0870.1050.1480.1553粉质粘土G/Gmax0.9800.9700.8400.7300.4000.2500.0700.0300.0120.0150.0370.0560.1120.1370.1700.1804卵石G/Gmax0.9900.9600.9000.8300.7100.5400

23、.3300.2200.0060.0080.0200.0330.0700.0850.1100.1205粉砂G/Gmax0.9800.9650.8850.8050.5600.4480.2200.1740.0050.0070.0200.0350.0800.1000.1200.1246细砂G/Gmax0.9200.8800.7000.5750.2600.1780.0580.0180.0150.0220.0560.0650.1040.1250.1450.1507砾石G/Gmax0.9900.9700.9000.8500.7000.5500.3200.2000.0040.0060.0190.0300.07

24、50.0900.1100.1208杂填土G/Gmax0.9600.9500.8000.7000.3000.2000.1500.1000.0250.0280.0300.0350.0800.1000.1100.1209基岩G/Gmax1.0001.0001.0001.0001.0001.0001.0001.0000.0040.0080.0100.0150.0210.0300.0360.046根据以上两组钻孔资料，在计算以上两组钻孔资料时所用到的土的动力学参数如表3.表3.5 各类土的动剪切模量比与阻尼比4、模型计算4.1 软夹层厚度关系计算 模型1-3为了研究软夹层的厚度对土层加速度反映峰值Ama

25、x、速度最大值Vmax及位移最大值Dmax、特征周期Tg及反应谱最大值Samax的影响，我们用剪切波速为120 m/s，密度为19 gm3的淤泥质粘土选用3,6,9三种厚度分别替换类场地埋深16米处的土层分别建立三个模型。具体分析模型如下所示：1、 模型一： 表4.1 模型一数据表序号土层名称土层厚度（m）波速(m/s)质量密度cm/g3土类1填土4.11301.982粉土2.31402.023粉质粘土3.52302.034 4粉砂22602.055淤泥质粘土31201.916粉土3.62002.027细砂3.52302.168砾石1.23802.279卵石6.94402.2410基岩5002

26、.392、 模型二：表4.2 模型二数据表序号土层名称土层厚度（m）波速(m/s)质量密度cm/g3土类1填土4.11301.982粉土2.31402.023粉质粘土3.52302.034 4粉砂0.92602.055淤泥质粘土61202.016粉土1.72002.027细砂3.52302.168砾石1.23802.279卵石6.94402.2410基岩5002.393、 模型三：表4.3 模型三数据表序号土层名称土层厚度（m）波速(m/s)质量密度cm/g3土类1填土4.11301.982粉土2.31402.023粉质粘土22302.034 4粉砂0.92602.055淤泥质粘土91202.016粉土1.22002.027细砂2.52302.168砾石1.23802.279卵石6.94402.2410基岩5002.394.2 软夹层埋深关系计算 模型4-6为了研究软夹层的埋深对土层加速度反应峰值Amax及其反应谱最大值Smax，速度最大值Vmax，深度Dmax的影响，以上述的类场地为研究目标将剪切波速为180 m/s，密度为19 gm3，厚度为4米的淤泥质粘土分别替换埋深为25,46.2,85的原有土层，来建立计算模型，具体分析模型如下所示：4、 模型四：表4.4 模型四数据表序