基于FLAC3D道路边坡地震响应研究.docx

基于FLAC3D道路边坡地震响应研究 摘要：本文通过分析地震动力学的基本理论与FLAC3D 动力分析模块中涉及的边界条件、阻尼形式、地震荷载等， 利用编写的程序对地震作用下土质边坡的地震动力响应进 行分析计算。结果表明，边坡对于地震的加速度、速度与位 移有放大效应，且边坡顶部的放大效应大于边坡底部，而边 坡的中下部位是土质边坡抗震的关键部位。 关键词：FLAC3D；地震；道路边坡 中图分类号：P315 文献标识码：A 文章编号： 1008-4428 (2012) 11-123 -04 前言 地震是地球内部介质局部发生急剧的破裂，产生的震 波，从而在一定范围内引起地面振动的现象。大地振动是地 震最直观、最普遍的表现。地震常常造成严重的人员伤亡及 建筑物破坏，地震对于道路边坡同样会带来比较严重的影 响。 FLAC3D ( Fast Lagrangian Analysis of Continua)是由美国Itasca公司采用ANSI C++语言编写开 发的快速拉格朗日差分分析软件，目前已发展到V4. 0版本。 它能够通过调整三维网格中的多面体单元来拟合实际的结 构，进行土质、岩石和其他材料的三维结构受力特性模拟和 塑性流动分析。 由于FLAC3D在模拟动荷载响应与大变形方面有一定的 优势，本文希望利用FLAC3D编写程序对地震作用下土质边 坡的地震动力响应进行分析计算。模拟边坡对于地震的加速 度、速度响应，确定道路边坡在地震荷载下的薄弱位置。 一、地震概论 （一） 地震波 地震所引起的地面振动是一种复杂的运动，它是由纵波 和横波共同作用的结果。第一类波的物理特性如声波。声波 是在空气里由交替的挤压（推）和扩张（拉）而传递。这种 类型的波叫P波，即纵波，它是首先到达的波。 地震产生的这种第二个到达的波叫S波。在S波通过时, 岩石的表现与在P波传播过程中的表现相当不同。因为S波 涉及剪切而不是挤压，使岩石颗粒的运动横过运移方向。这 些岩石运动可在一垂直向或水平面里，它们与光波的横向运 动相似。P和S波同时存在使地震波列成为具有独特的性质 的组合，使之不同于光波或声波的物理表现。 （二） 地震动力学理论 由于地面运动所产生的惯性力称为地震作用或地震力。 因为地震力不仅和地面运动有关，而且还和结构的动力特性 有关。而地面运动的特征（强度、频谱特征、持续时间等） 受震源机制、地质条件等因素的影响，结构本身又是一个十 分复杂的力学系统，加之结构和基础间的相互作用，其中许 多因素是搞不清楚的。因此地震作用的研究经历了艰难的历 程。 由于问题的复杂性，地震作用理论经历了一个由粗略到 科学的发展过程[1]:主要有：静力理论；定函数理论；反 应谱理论；时程分析理论等。随着计算机技术的发展，计算 机仿真模拟技术在地震动力研究中发挥了越来越大的作用。 二、FLAC3D地震模拟分析方法 FLAC3D可以进行非线性动力反应分析，而且具有强大的 动力分析功能，这需要在FLAC3D程序中启用动力分析模块 (CONFIG dynamic).,本节将围绕人工边界问题、地震波输 入问题和材料阻尼问题展开一系列的应用研究工作。 地震非线性动力分析方法 在地震工程中，在模拟波的传播以及结构的动态相互作 用时，等线性法已经获得广泛使用。而FLAC3D采用的是完 全非线性方法[2, 3]o 在等线性法中，假定模型中不同区域阻尼比和剪切模量 的初始值。然后使用线性化的分析方法参照试验室得出的阻 尼比及切线模量同循环剪切应变振幅间的关系曲线，对每一 单元的最大循环剪切应变加以记录。并用其确定新的阻尼和 模量值，计算过程中时常用到一些经验缩放比例系数。然后, 将阻尼比和剪切模量的新值用于新模型的数值分析。整个过 程重复数次，直至特性不发生进一步的变化。这就找到阻尼 和模量所谓的“应变相容”值，用这些值进行模拟代表实际 现场的响应。 相比之下，完全非线性法只用一次运行，这是因为求解 进程中，每个单元直接遵循应力 应变关系随时间的非线 性变化。若使用了适当的非线性法则，能自动模拟阻尼和表 观模量对应变水平的相关性。 （二）边界条件 动力学分析中，模型周围边界条件的选择是一个主要内 容，因为边界会存在波的反射，对动力分析的结果产生影响。 把分析模型的范围设置得越大，分析结果就越好，但较大的 模型会导致巨大的计算负担o FLAC3D中为求解动力问题提供 了静态（粘性）边界和自由场边界两种边界条件来减少模型 边界上的波的反射。本文在模拟过程中模型周围选用了自由 场边界条件。 对于诸如边坡之类的地面结构地震反应进行数值分析 需要对邻近基础部分区域的材料进行离散化。地震能量通常 由通过下部材料向上传播的面波来表示，在模型各边处的边 界条件必须考虑没有结构时仍存在着的自由场运动。如图 2-1所示，