FLAC3D动力分析-供QQ群.pptVIP

FLAC / FLAC3D基本原理与应用 陈育民 博士 第三讲 FLAC3D动力分析、自定义本构、结构单元 FLAC3D非线性动力分析 为什么要用FLAC做动力分析？ FLAC 可以模拟体系（土，岩石，结构，流体）受到的外部动力荷载（比如地震）或内部动力荷载（比如基础振动、爆炸）。 可以计算塑性引起的永久变形以及孔隙水压力的消散。 土动力学中常用的等效线性方法无法直接处理上述问题。 动力模拟的3个重要问题 动力荷载与边界条件 材料响应与阻尼 土体液化 动力荷载 动力输入的类型 加速度时程 速度时程 应力(压力)时程 力时程 APPLY INTERIOR (内部) TABLE FISH Quiet边界 静态(quiet,粘性)边界 Lysmer and Kuhlemeyer(1969) 模型边界法向和切向设置独立的阻尼器 性能 对于法向p波和s波能很好的吸收 对于倾斜入射的波和Rayleigh波也有所吸收，但存在反射 人工边界仍应当足够远 Quiet边界应用 内部振动(如隧道中的列车振动问题)? 动力荷载直接施加在节点上 使用Quiet边界减小人工边界上的反射 不需要FF边界 外部荷载的底部边界? 软土地基上的地震荷载不适合用加速度或速度边界条件 使用应力条件t = -2Csrvs 地震底部输入的侧向边界? 扭曲了入射波 Free-field边界 Cundall et al. (1980) 自由场网格与主体网格的耦合粘性阻尼器，自由场网格的不平衡力施加到主体网格边界上 设置条件 底部水平，重力方向为z向 侧面垂直，法向分别为x, y向 其他边界条件在APPLY ff之前 Free-field边界 APPLY ff将边界上单元的属性、条件和变量全部转移ff单元上； 设置以后主体网格上的改动将不会被FF边界所响应 可存在任意的本构模型以及流体耦合(仅竖向) FF边界进行小变形计算，主体网格可大变形，FF边界上的变形要相对较小 存在attach的边界将不能设置FF边界 边界上的Interface将不能连续 动力边界设置需在FF边界设置之前 Free-field边界与动力荷载 模型底部边界 fix——施加速度或加速度荷载——刚性边界 Free——施加应力时程荷载——柔性边界 对于软弱的地基不适合施加速度（加速度荷载），而应当施加应力荷载 2. 材料响应与阻尼 材料响应 试验得到的阻尼比、割线模量随循环剪应变的曲线 等效线性方法 等效线性方法的特点 使用振动荷载的平均水平来估算每个单元的线性属性，并在振动过程中保持不变。在弱震阶段，单元会变得阻尼过大而刚度太小；在强震阶段，单元将会变得阻尼太小而刚度太大。对于不同部位不同运动水平的特性存在空间变异性。 不能计算永久变形。等效线性方法模型在加荷与卸荷时模量相同，不能计算土体在周期荷载作用下发生的剩余应变或位移。 塑形屈服模拟不合理。在塑性流动阶段，普遍认为应变增量张量是应力张量的函数，称之为“流动法则”。然而，等效线性方法使用的塑性理论认为应变张量（而不是应变增量张量）是应力张量的函数。因此，塑性屈服的模拟不合理。 大应变时误差大。等效线性方法所用割线模量在小应变时与非线性的切线模量很相近，但在大应变时二者相差很大，偏于不安全。 本构模型单一。等效线性方法本身的材料本构模型包括了应力应变的椭圆形方程，这种预设的方程形式减少了使用者的选择性，但却失去了选择其它形状的适用性。方法中使用迭代程序虽然部分考虑了不同的试验曲线形状，但是由于预先设定了模型形式，所以不能反映与频率无关的滞回圈。另外，模形是率无关的，因此不能考虑率相关性。 完全非线性分析方法 FLAC3D采用完全非线性分析方法，基于显式差分方法，使用由周围区域真实密度得出的网格节点集中质量，求解全部运动方程。 完全非线性分析方法的特点 可以遵循任何指定的非线性本构模型。如果模型本身能够反映土体在动力作用下的滞回特性，则程序不需要另外提供阻尼参数。如果采用Rayleigh阻尼或局部（local）阻尼，则在动力计算中阻尼参数将保持不变。 采用非线性的材料定律，不同频率的波之间可以自然地出现干涉和混合，而等效线性方法做不到这一点。 由于采用了弹塑性模型，因此程序可以自动计算永久变形。 采用合理的塑性方程，使得塑性应变增量与应力相联系。 可以方便地进行不同本构模型的比较。 可以同时模拟压缩波和剪切波的传播及两者耦合作用时对材料的影响。在强震作用下，这种耦合作用的影响很重要，比如在摩擦型材料中，法向应力可能会动态地减小从而降低土体的抗剪强度。 使用弹塑性模型 附加考虑的因素： 阻尼，对于屈服面以下应力的循环 体积应变积累，是循环周数与幅值的函数 模量衰减，基于平均应变水平的表格 弹塑性模型 FLAC3D中的阻尼比 弹塑