FLAC2D动力分析要点.pdf

FLAC2D 动力反应分析要点 FLAC 可以模拟岩土体在外部（如地震）或内部（如风、爆炸、地铁振动）荷载作用下 的完全非线性响应，可以适用于土动力学、岩石动力学等计算。 一、与等效线性方法的关系 在岩土地震工程中，等效线性方法广泛应用于计算地基土体中波的传播及土—结构的动 力相互作用；而FLAC 采用的完全非线性方法没有获得广泛使用，因此需要对这两种方法之 间的差异做简要介绍。 1.1 等效线性法的特点 等效线性方法的基本原理是，假定土体是粘弹性体，参照实验室得到的切线模量及阻尼 比与剪应变幅值的关系曲线，对地震中每一单元的阻尼和模量重新赋值。有如下的特点：  使用振动荷载的平均水平来估算每个单元的线性属性，并在振动过程中保持不变。 在弱震阶段，单元会变得阻尼过大而刚度太小；在强震阶段，单元将会变得阻尼太小而刚度 太大。对于不同部位不同运动水平的特性存在空间变异性。  不能计算永久变形。等效线性方法模型在加荷与卸荷时模量相同，不能计算土体在 周期荷载作用下发生的残余应变。  塑形屈服模拟不合理。在塑性流动阶段，普遍认为应变增量张量是应力张量的函数， 称之为“流动法则”；等效线性方法使用的塑性理论认为应变张量是应力张量的函数。因此， 塑性屈服的模拟不合理。  大应变时误差大。等效线性方法所用割线模量在小应变时与非线性的切线模量很相 近，但在大应变时二者相差很大，偏于不安全。  本构模型单一。等效线性方法本身的材料本构模型包括了应力应变的椭圆形方程， 这种预设的方程形式减少了使用者的选择性，但却失去了选择其它形状的适用性。方法中使 用迭代程序虽然部分考虑了不同的试验曲线形状，但是由于预先设定了模型形式，所以不能 反映与频率无关的滞回圈。另外，模形是频率无关的，因此不能考虑频率相关性。 1.2 FLAC 非线性方法的特点 FLAC 采用基于显式差分方法完全非线性分析方法，使用由周围区域真实密度得出的网 格节点集中质量，求解全部动力方程。相对于等效线性方法，完全非线性分析方法主要有以 下优点：  可以遵循任何指定的非线性本构模型。如果模型本身能够反映土体在动力作用下的 滞回特性，则程序不需要另外提供阻尼参数。如果采用 Rayleigh 阻尼或局部（local）阻尼， 则在动力计算中阻尼参数将保持不变。  采用非线性的材料定律，不同频率的波之间可以自然地出现干涉和叠加，而等效线 性方法做不到这一点。  采用了弹塑性模型，此程序可以自动计算永久变形。  采用合理的塑性方程，使得塑性应变增量与应力相联系。  可以方便地进行不同本构模型的比较。  可以同时模拟压缩波和剪切波的传播及两者耦合作用时对材料的影响。在强震作用 下，这种耦合作用的影响很重要，比如在摩擦型材料中，法向应力可能会动态地减小从而降 低土体的抗剪强度。 二、动力时间步 动力计算中临界计算时间步的计算如下：  Az  ∆tcrit min   （11-1） L C   d p  其中， 为p 波波速，与材料的体积模量K 和剪切模量G 有关，可以表示为： C p K +4 / 3G C （11-2 ） p ρ 为矩形区表面积， 为对角线长度， 表示遍历所有的单元，包括结构单元和 Az Ld min{} 接触面单元。 由于式（11-1 ）只是临界时间步的一个估计值，在使用中采用了一个安全系数，乘以 0.5 。因此，当采用无刚度比例的阻尼时，动力分析的时间步为： ∆td =∆tcrit 2 （11-3 ） 如果采用了