FLAC动力分析.docVIP

非线性动力反应分析 FLAC / FLAC3D可以进行非线性动力反应分析，而且具有强大的动力分析功能。本章以FLAC3D为例，详细介绍了动力分析过程中的边界条件、阻尼形式、荷载要求等，并通过一些实例对个别问题做了详细解答。 本章要点： FLAC动力分析与等效线性方法的差别 动力分析时间步的确定方式及影响因素 动态多步的概念 动力荷载的形式及施加方法 动力边界条件的类型及适用条件 地震荷载输入的要点 三种阻尼形式的概念、参数确定及适用条件 网格尺寸的要求 输入荷载的校正 地震液化的模拟 完全非线性动力分析的步骤 11.1 概述 FLAC / FLAC3D可以进行三维完全动力分析FLAC/FLAC3D中的动力分析功能是可选模块，需要在程序中添加动力分析模块才可以进行。 FLAC3D中在动力分析前需要采用以下的命令： CONFIG dynamic 对于FLAC，在程序开始时的Model Options对话框中选择Dynamic复选框。 FLAC / FLAC3D中的动力分析并不是只能孤立进行的，还可以与其他FLAC/FLAC3D元素进行耦合。 （1）与结构，可以用来进行土与结构的动力相互作用。FLAC和FLAC3D可以模拟岩土体在外部（如地震）或内部（如风、爆炸、地铁振动）荷载作用下的完全非线性响应，因此可以适用于土动力学、岩石动力学等学科的计算。 本章FLAC的动力分析可以参照执行。 注意：FLAC和FLAC3D的动力计算十分复杂，读者在阅读本章内容之前要对FLAC3D的静力计算、流体计算十分熟悉，具体可以参阅本书的第7章和第12章的内容。 对于初次接触FLAC3D动力计算的读者，大多数都会提以下2个问题： （1）FLAC3D动力分析与一般的等效线性方法有什么区别？ （2）FLAC3D动力分析怎么会采用静力本构模型，比如Mohr-Coulomb模型？ 下面就这两个问题展开初步的讨论。 11.1.1 与等效线性方法的关系 在岩土地震工程中，等效线性方法广泛应用于计算地基土体中波的传播及土与结构的动力相互作用。该方法已被工程师、科研人员广泛接受。而FLAC3D采用的完全非线性方法没有获得广泛使用，因此需要对这两种方法之间的差异做简要介绍。 1. 等效线性方法的特点 大应变时误差大。等效线性方法所用割线模量在小应变时与非线性的切线模量很相近，但在大应变时二者相差很大，偏于不安全。 本构模型单一。等效线性方法本身的材料本构模型包括了应力应变的椭圆形方程，这种预设的方程形式减少了使用者的选择性，但却失去了选择其它形状的适用性。方法中使用迭代程序虽然部分考虑了不同的试验曲线形状，但是由于预先设定了模型形式，所以不能反映与频率无关的滞回圈。另外，模形是率无关的，因此不能考虑率相关性。 2. FLAC3D非线性方法的特点FLAC3D动力计算采用的本构模型 Coulomb是静力本构模型，不适合用于动力分析，而应当采用更合适的Hardin模型。 其实这是对FLAC3D动力计算的误解。FLAC3D的原理是求解动力方程，所以从其算法上来说，不管是进行静力分析还是动力分析，其实质都是求解运动方程。只是对于静力分析而言，采用了特定的阻尼方式以达到快速收敛的目的。所以，有的场合将FLAC3D的静力分析方法称为“拟动力方法”。相应的，FLAC3D在进行动力分析时，通过求解动力方程理所当然地可以得到合适的动力问题解答。对于本构模型的选择，主要是描述单元的应力-应变关系，如果是弹塑性的，则考虑的是单元的屈服准则、流动法则等。 等效线性方法考虑土体的滞后性常常是通过将骨干曲线进行变换，比如Masing二倍法，而在FLAC3D的动力分析中，滞后性是通过阻尼来考虑，通过设置合适的阻尼形式和阻尼参数，同样可以描述土体在动力作用下的滞回曲线和滞回圈。 因此，FLAC3D动力分析中采用的本构模型可以选取任意模型，其参数也是对应静力本构模型的参数，关键是要设置合适的阻尼形式、阻尼参数、边界条件等，这些内容将在本章的后续内容中进行讲解。 11.2 动力时间步 动力计算中临界计算时间步的计算如下： （11-1） 其中，为p波波速，与材料的体积模量K和剪切模量G有关，可以表示为： （11-2） 为四面体子单元（sub-zone）的体积，为与四面体子单元相关的最大表面积，表示遍历所有的单元，包括结构单元和接触面单元。 由于式（11-1）只是临界时间步的一个估计值，因此在使用中采用了一个安全系数，乘以0.5。因此，当采用无刚度比例的阻尼时，动力分析的时间步为： （11-3） 如果采用了刚度比例的阻尼，那么为了保持数值稳定性，时间步必须减小。Belytschko（1983）提出了一个临界时间步的公式，其中考虑了刚度比例阻尼的影响。 （11-4） 其中，为系统