FLAC动力分析论文.docVIP

PAGE 1 非线性动力反应分析 FLAC / FLAC3D可以进行非线性动力反应分析，而且具有强大的动力分析功能。本章以FLAC3D为例，详细介绍了动力分析过程中的边界条件、阻尼形式、荷载要求等，并通过一些实例对个别问题做了详细解答。 本章要点： FLAC动力分析与等效线性方法的差别 动力分析时间步的确定方式及影响因素 动态多步的概念 动力荷载的形式及施加方法 动力边界条件的类型及适用条件 地震荷载输入的要点 三种阻尼形式的概念、参数确定及适用条件 网格尺寸的要求 输入荷载的校正 地震液化的模拟 完全非线性动力分析的步骤 MACROBUTTON MTEditEquationSection2 Equation Section 11 SEQ MTEqn \r \h SEQ MTSec \r 11 \h 11.1 概述 FLAC / FLAC3D可以进行二维或三维的完全动力分析，FLAC/FLAC3D中的动力分析功能是可选模块，需要在程序中添加动力分析模块才可以进行。 FLAC3D中在动力分析前需要采用以下的命令： CONFIG dynamic 对于FLAC，在程序开始时的Model Options对话框中选择Dynamic复选框。 FLAC / FLAC3D中的动力分析并不是只能孤立进行的，还可以与其他FLAC/FLAC3D元素进行耦合。 （1）与结构单元相耦合，可以用来进行土与结构的动力相互作用。 （2）与流体计算相耦合，可以模拟动力作用下土体孔隙水压力的上升直至土体液化。 （3）与热力学计算相耦合，可以计算热力荷载和动力荷载的共同作用。 （3）采用大变形计算模式，可以分析岩土体在动力荷载作用下发生的大变形。 FLAC和FLAC3D可以模拟岩土体在外部（如地震）或内部（如风、爆炸、地铁振动）荷载作用下的完全非线性响应，因此可以适用于土动力学、岩石动力学等学科的计算。 本章将以FLAC3D为例讨论动力计算的相关内容，FLAC的动力分析可以参照执行。 注意：FLAC和FLAC3D的动力计算十分复杂，读者在阅读本章内容之前要对FLAC3D的静力计算、流体计算十分熟悉，具体可以参阅本书的第7章和第12章的内容。 对于初次接触FLAC3D动力计算的读者，大多数都会提以下2个问题： （1）FLAC3D动力分析与一般的等效线性方法有什么区别？ （2）FLAC3D动力分析怎么会采用静力本构模型，比如Mohr-Coulomb模型？ 下面就这两个问题展开初步的讨论。 11.1.1 与等效线性方法 在岩土地震工程中，等效线性方法广泛应用于计算地基土体中波的传播及土与结构的动力相互作用。该方法已被工程师、科研人员广泛接受。而FLAC3D采用的完全非线性方法没有获得广泛使用，因此需要对这两种方法之间的差异做简要介绍。 1. 等效线性方法的特点 等效线性方法的基本原理是，假定土体是粘弹性体，参照实验室得到的切线模量及阻尼比与剪应变幅值的关系曲线，对地震中每一单元的阻尼和模量重新赋值。目前用于土动力分析的等效线性模型已有数种，根据骨干曲线的形状可以分为双直线模型、Ramberg-Osgood模型、Hardin-Drnevich模型等，其中又以Hardin模型使用最多。等效线性方法有如下的特点： 使用振动荷载的平均水平来估算每个单元的线性属性，并在振动过程中保持不变。在弱震阶段，单元会变得阻尼过大而刚度太小；在强震阶段，单元将会变得阻尼太小而刚度太大。对于不同部位不同运动水平的特性存在空间变异性。 不能计算永久变形。等效线性方法模型在加荷与卸荷时模量相同，不能计算土体在周期荷载作用下发生的剩余应变或位移。 塑形屈服模拟不合理。在塑性流动阶段，普遍认为应变增量张量是应力张量的函数，称之为“流动法则”。然而，等效线性方法使用的塑性理论认为应变张量（而不是应变增量张量）是应力张量的函数。因此，塑性屈服的模拟不合理。 大应变时误差大。等效线性方法所用割线模量在小应变时与非线性的切线模量很相近，但在大应变时二者相差很大，偏于不安全。 本构模型单一。等效线性方法本身的材料本构模型包括了应力应变的椭圆形方程，这种预设的方程形式减少了使用者的选择性，但却失去了选择其它形状的适用性。方法中使用迭代程序虽然部分考虑了不同的试验曲线形状，但是由于预先设定了模型形式，所以不能反映与频率无关的滞回圈。另外，模形是率无关的，因此不能考虑率相关性。 2. FLAC3D非线性方法的特点 FLAC3D采用完全非线性分析方法，基于显式差分方法，使用由周围区域真实密度得出的网格节点集中质量，求解全部运动方程。相对于等效线性方法而言，完全非线性分析方法主要有以下优点： 可以遵循任何指定的非线性本构模型。如果模型本身能够反映土体在动力作用下的滞回特性，则程序不需要另外提供阻尼参数