FLAC3D动力分析-陈育民

1FLAC/FLAC3D基本原理与应用陈育民博士河海大学科学研究院2第三讲FLAC3D动力分析、自定义本构、结构单元3FLAC3D非线性动力分析非常复杂！SaidbyProf.PeterCundall4为什么要用FLAC做动力分析？FLAC可以模拟体系（土，岩石，结构，流体）受到的外部动力荷载（比如地震）或内部动力荷载（比如基础振动、爆炸）。可以计算塑性引起的永久变形以及孔隙水压力的消散。土动力学中常用的等效线性方法无法直接处理上述问题。5动力模拟的3个重要问题1.动力荷载与边界条件2.材料响应与阻尼3.土体液化6动力荷载动力输入的类型加速度时程速度时程应力(压力)时程力时程APPLYINTERIOR(内部)TABLEFISH7Quiet边界静态(quiet,粘性)边界LysmerandKuhlemeyer(1969)模型边界法向和切向设置独立的阻尼器性能对于法向p波和s波能很好的吸收对于倾斜入射的波和Rayleigh波也有所吸收，但存在反射人工边界仍应当足够远8Quiet边界应用内部振动(如隧道中的列车振动问题)动力荷载直接施加在节点上使用Quiet边界减小人工边界上的反射不需要FF边界外部荷载的底部边界软土地基上的地震荷载不适合用加速度或速度边界条件使用应力条件t=-2Csrvs地震底部输入的侧向边界扭曲了入射波quietquietquiet9Free-field边界Cundalletal.(1980)自由场网格与主体网格的耦合粘性阻尼器，自由场网格的不平衡力施加到主体网格边界上设置条件底部水平，重力方向为z向侧面垂直，法向分别为x,y向其他边界条件在APPLYff之前npnCvr相当于一个阻尼器10Free-field边界APPLYff将边界上单元的属性、条件和变量全部转移ff单元上；设置以后主体网格上的改动将不会被FF边界所响应可存在任意的本构模型以及流体耦合(仅竖向)FF边界进行小变形计算，主体网格可大变形，FF边界上的变形要相对较小存在attach的边界将不能设置FF边界边界上的Interface将不能连续动力边界设置需在FF边界设置之前11Free-field边界与动力荷载模型底部边界fix——施加速度或加速度荷载——刚性边界Free——施加应力时程荷载——柔性边界对于软弱的地基不适合施加速度（加速度荷载），而应当施加应力荷载2SsCvtrNotethatthereisafactorof2becausetheinputenergydividesintoadownward-&upward-propagatingwave.121.连续的非线性，表观模量随着应变的增大而降低2.对所有循环应变等级均存在滞回特性，因此导致随着循环应变的增加阻尼比增大。阻尼是率相关的。3.对于复杂波形的各个成分都产生阻尼。4.剪切应变会产生的体积应变，相应的，随着剪应变循环次数的增加体积应变逐渐积累。-3.0E+05-2.0E+05-1.0E+050.0E+001.0E+052.0E+053.0E+054.0E+05-0.1%0.0%0.1%0.2%Shearstrain%Shearstress2.材料响应与阻尼13土体在循环荷载作用下呈现出模量衰减和能量消散的特点，那么如何用非线性数值方法对其进行模拟呢？Nonlinearcharacteristicsofsoils(MartinandSeed,1979)材料响应14试验得到的阻尼比、割线模量随循环剪应变的曲线0.00010.0010.010.11ShearStrainAmplitude(%)01020304050DampingRatio(%)0.00010.0010.010.11ShearStrainAmplitude(%)0.00.10.20.30.40.50.60.70.80.91.0NormalizedShearModulus,G/GmaxMid-RangeSandCurve(Seed&Idriss,1970)SandFillInland:Friction=32,hr=0.47,Go=440SandFillunderRockDike:Friction=30,hr=0.43,Go=44015•等效线性方法是岩土地震工程中模拟波的传播的最常用的方法。•假定土体是粘弹性体，参照实验室得到的切线模量及阻尼比与剪应变幅值的关系曲线，对地震中每一单元的阻尼和模量重新赋值。等效线性方法Iterationtowardstrain-compatibleshearmodulusanddampingratio(afterKramer,1996)16等效线性方法的特点1.使用振动荷载的平均水平来估算每个单元的线性属性，并在振动过程中保持不变。在弱震阶段，单元会变得阻尼过大而刚度太小；在强震阶段，单元将会变得阻尼太小而刚度太大。对于不同部位不同运动水平的特性存在空间变异性。2.不能计算永久变形。等效线性方法模型在加荷与卸荷时模量相同，不能计算土体在周期荷载作用下发生的剩余应变或位移。3.塑形屈服模拟不合理。在塑性流动阶段，普遍认为应变增量张量是应力张量的函数，称之为“流动法则”。然而，等效线性方法使用的塑性理论认为应变张量（而不是应变增量张量）是应力张量的函数。因此，塑性屈服的模拟不合理。4.大应变时误差大。等效线性方法所用割线模量在小应变时与非线性的切线模量很相近，但在大应变时二者相差很大，偏于不安全。5.本构模型单一。等效线性方法本身的材料本构模型包括了应力应变的椭圆形方程，这种预设的方程形式减少了使用者的选择性，但却失去了选择其它形状的适用性。方法中使用迭代程序虽然部分考虑了不同的试验曲线形状，但是由于预先设定了模型形式，所以不能反映与频率无关的滞回圈。另外，模形是率无关的，因此不能考虑率相关性。17完全非线性分析方法FLAC3D采用完全非线性分析方法，基于显式差分方法，使用由周围区域真实密度得出的网格节点集中质量，求解全部运动方程。18完全非线性分析方法的特点1.可以遵循任何指定的非线性本构模型。如果模型本身能够反映土体在动力作用下的滞回特性，则程序不需要另外提供阻尼参数。如果采用Rayleigh阻尼或局部（local）阻尼，则在动力计算中阻尼参数将保持不变。2.采用非线性的材料定律，不同频率的波之间可以自然地出现干涉和混合，而等效线性方法做不到这一点。3.由于采用了弹塑性模型，因此程序可以自动计算永久变形。4.采用合理的塑性方程，使得塑性应变增量与应力相联系。5.可以方便地进行不同本构模型的比较。6.可以同时模拟压缩波和剪切波的传播及两者耦合作用时对材料的影响。在强震作用下，这种耦合作用的影响很重要，比如在摩擦型材料中，法向应力可能会动态地减小从而降低土体的抗剪强度。19使用弹塑性模型附加考虑的因素：阻尼，对于屈服面以下应力的循环体积应变积累，是循环周数与幅值的函数模量衰减，基于平均应变水平的表格最简单的弹塑性模型往往在描述累计塑性应变方面具有很好的效果，但是对于加速度放大系数的估算上效果不好。20简单的理想弹塑性本构模型仅仅在发生屈服时才会出现滞回特性strainstress注意：即使这样粗糙的模型也能够作出连续的阻尼比和模量衰减曲线。在屈服条件下会产生体积改变，但通常都是剪胀。弹塑性模型21FLAC3D中的阻尼比1.弹塑性本构模型使用瑞利（粘性）阻尼2.弹塑性本构模型使用滞后（HD）阻尼.3.复杂本构模型拥有连续的屈服应力应变关系和对应的加卸载响应（Wang，UBCSand）22瑞利阻尼瑞利阻尼最初应用于结构和弹性体的动力计算中，以减弱系统的自然振动模式的振幅。在计算时，假设动力方程中的阻尼矩阵C与刚度矩阵K和质量矩阵M有关：瑞利阻尼中的质量分量相当于连接每个节点和地面的阻尼器，而刚度分量则相当于连接单元之间的阻尼器。虽然两个阻尼器本身是与频率有关的，但是通过选取合适的系数，可以在有限的频率范围内近似获得频率无关的响应。23frequencyratioofdampingtocritical混合仅有刚度分量仅有质量分量采用叠加的方法得到的阻尼比在较大的频率范围内保持定值（3:1）Combinedcurvereachesminimumat:瑞利阻尼24假设弹性模型计算土石坝不同材料的功率谱曲线(assumingelasticmaterial)frequencyfrequencyfrequencyfrequency246810120.5001.0001.5002.0002.5003.0003.5004.0004.500(10)-09frequency246810121.0002.0003.0004.0005.000(10)-10246810120.5001.0001.5002.0002.5003.0003.500(10)-09246810120.4000.8001.2001.6002.000(10)-09246810120.5001.0001.5002.0002.500(10)-09中心频率的选择25Shake9100.20.40.60.811.20.00010.0010.010.1110strain-%ModulusreductionfactorShake91Shake910510152025300.00010.0010.010.1110strain-%Dampingratio-%Shake91根据65%最大应变来选择阻尼比和模量衰减参参数-20-15-10-505101520(10)-04-2.000-1.0000.0001.0002.0003.000(10)03阻尼比的选择26必须指定中心频率，有时须在一些充满矛盾的数据当中进行选择（场地响应或地震平均频率）刚度比例项会导致时间步减小瑞利阻尼的缺点27FLACVersion5.0andFLAC3DVersion3.0提供了滞后阻尼功能，该项阻尼独立于材料模型之外，实质是包含了一个与应变幅值相关的切线剪切模量乘子00.20.40.60.811.20.00010.0010.010.1110Cyclicstrain%Modulusreductionfactor如果割线模量是由衰减曲线来确定的，那么切线模量可以表示为：secantmodulustangentmodulusshearstressshearstrainstMMtFromSeed&Idriss(1970)ttddMMddMMssts)(oG/ttGo=small-strainshearmodulus滞后阻尼28tGMG表观应变是偏应变，根据先前反转的点进行累计得到的。使用堆栈的数据结构可以描述应变反转现象。FLAC(Version4.00)LEGEND12-Feb-0315:39step3700HISTORYPLOTY-axis:Ave.SXY(1,1)X-axis:Xdisplacement(1,2)-40-2002040(10)-05-2.000-1.0000.0001.0002.000(10)+04JOBTITLE:ItascaConsultingGroup,Inc.Minneapolis,MinnesotaUSA因此，在“迷你滞回圈”中的能量损失与主滞回圈的计算是一样的。elasticmodelwithhysteresisdamping00.20.40.60.811.20.00010.0010.010.1110Cyclicstrain%Modulusreductionfactor滞后阻尼特定的曲线可以使用滞后阻尼来描述，并在本构模型计算中根据应变的大小获得不同的切线模量。29滞后阻尼的特点1.可以直接采用动力试验中的模量衰减曲线；2.相对于瑞利阻尼而言，滞后阻尼不影响动力计算的时间步；3.可以应用于任意的材料模型，且可以与其它阻尼格式同时使用。一个缺点：模量衰减曲线一致，但阻尼比曲线存在差别00.20.40.60.81