《FLAC3D在岩土工程中的应用-陈育民-东南大学讲座课件》由会员分享,可在线阅读,更多相关《FLAC3D在岩土工程中的应用-陈育民-东南大学讲座课件(77页珍藏版)》请在金锄头文库上搜索。
1、FLAC3D在岩土工程中的应用,陈育民,左威龙 导师:刘汉龙 教授,河海大学岩土工程研究所 2007-4-16,东南大学交通学院岩土工程研究所,主要内容,软件介绍 动力分析 Dynamic Option 桩-土相互作用分析 Interface 隧道分析 Structure Element 流固耦合分析 Fluid-Mechanical Interaction 学习方法及经验介绍,主要内容,软件介绍 动力分析 Dynamic Option 桩-土相互作用分析 Interface 隧道分析 Structure Element 流固耦合分析 Fluid-Mechanical Interaction
2、学习方法及经验介绍,软件介绍,Fast Lagrangian Analysis of Continua 美国Itasca咨询公司开发2D程序(1986) 1990年代初引入中国 有限差分法(FDM) DOS版2.0 2.1 3.0 3.1 Itasca其他软件,软件介绍,应用 岩土力学分析,例矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定性分析、采矿巷道稳定性研究等 岩土工程、交通工程、采矿工程、水利工程、地质工程 特色 大应变模拟 完全动态运动方程使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍 动力分析功能地震工程 衬砌功能地下工程 可开发功能 FISH VC+,丰富的本构模型,软
3、件介绍,有限差分法(FDM) 古老的方法(上世纪40年代) 用差分格式转化控制方程中的微商格式 流体力学;土工渗流问题;固结 FDM 建立网格(前处理) model elas ;材料参数 prop bulk 3e6 shear 1e6 ini dens 2000 ;初始条件 fix z ran z -.1 .1 ;边界条件 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve ;求解 app nstr -10e4 ran z 3 x 1 2 y 1
4、 2 solve plo con zd ;后处理 切片功能,RUN FLAC3D,FLAC3D3.1的新特征,64-bit Version (64位计算) Parallel Processing on Multiprocessor Computers (并行算法) Two-Dimensional Grid Extrusion Tool Nodal Mixed Discretization (节点混合离散方法) Embedded Structural Liner (内置衬砌单元) 挡土墙的模拟 On-board Help File (内置的帮助系统),主要内容,软件介绍 动力分析 Dynamic
5、 Option 桩-土相互作用分析 Interface 隧道分析 Structure Element 流固耦合分析 Fluid-Mechanical Interaction 学习方法及经验介绍,完全非线性的动力分析,与等效线性方法的比较 动力荷载 动力边界条件 力学阻尼与滞回阻尼 地震波的调整 动孔压的生成,FLAC与等效线性方法,等效线性方法是岩土地震工程中的常用方法 动本构模型 等效线性模型 剪切模量的降级曲线 阻尼比随剪应变的变化 FLAC 常规模型(MC) 多种频率成分的干涉和混合 永久变形计算 弹塑性计算,动力荷载,动力输入的类型 加速度时程 速度时程 应力(压力)时程 力时程 AP
6、PLY INTERIOR (内部) TABLE FISH,Quiet边界,静态(quiet,粘性)边界 Lysmer and Kuhlemeyer(1969) 模型边界法向和切向设置独立的阻尼器 性能 对于法向p波和s波能很好的吸收 对于倾斜入射的波和Rayleigh波也有所吸收,但存在反射 人工边界仍应当足够远,Quiet边界应用,内部振动(如隧道中的列车振动问题) 动力荷载直接施加在节点上 使用Quiet边界减小人工边界上的反射 不需要FF边界 外部荷载的底部边界 软土地基上的地震荷载不适合用加速度或速度边界条件 使用应力条件t = -2Csrvs 地震底部输入的侧向边界 扭曲了入射波,q
7、uiet,quiet,quiet,Free-field边界,Cundall et al. (1980) 自由场网格与主体网格的耦合粘性阻尼器,自由场网格的不平衡力施加到主体网格边界上 设置条件 底部水平,重力方向为z向 侧面垂直,法向分别为x, y向 其他边界条件在APPLY ff之前,自由场边界示意图,Free-field边界,APPLY ff将边界上单元的属性、条件和变量全部转移ff单元上; 设置以后主体网格上的改动将不会被FF边界所响应 可存在任意的本构模型以及流体耦合(仅竖向) FF边界进行小变形计算,主体网格可大变形,FF边界上的变形要相对较小 存在attach的边界将不能设置FF边
8、界 边界上的Interface将不能连续 动力边界设置需在FF边界设置之前,力学阻尼,瑞利(rayleigh)阻尼 假设阻尼与质量、刚度的线性关系 参数确定简单(等价平均应变=60%*emax) 中心频率(共振计算,地震平均频率) 临界阻尼比 缺点:计算速度慢 局部(local)阻尼 FLAC3D的静力分析阻尼 参数简单 适合简单情况,滞回阻尼(Hysteretic Damping),模拟岩土介质的动模量衰减曲线 initial damp hysteretic name sig3 (3参数) sig4 (4参数) Hardin(1参数) (哈丁模型) Default(2参数) 优点 直接采用模
9、量降级曲线 阻尼比不会影响时间步 缺点 输出的曲线会不一致,滞回阻尼(Hysteretic Damping),低循环应变下得到的阻尼比要小于试验结果,这会导致低级的噪声,尤其在高频情况下。可以在中心频率上增加一个小量的Rayleigh阻尼(0.2%刚度比例),这样也不会降低时步; 若初始应力不为0,剪应力-剪应变曲线可能不匹配。因此在生成初始应力时就要调用Hyst阻尼; Hyst阻尼不仅会增加能量损失,还会导致在大循环应变下的平均剪切模量的降低,在输入波的基频接近共振频率的时候,可能会导致动力反应幅值的增大; Hyst阻尼之前要做一次弹性无阻尼求解,以获得发生循环应变的最大水平,若循环应变过大
10、导致剪切模量过多的降低,那么用Hyst阻尼是有问题的; 即使应变较小,使用屈服模型也会增大应变,因此若有广泛屈服的现象,则使用屈服模型,不用Hyst阻尼 参考了Itasca的中国培训资料,地震波的调整,基线校正 对于地震分析的加速度时程,其积分得到的速度和位移应归0 美国地质调查研究所 Basic Strong-Motion Accelerogram Processing Software (BAP) 对网格施加一个固定速度从而使残余的位移变为0 动力荷载的频率与单元尺寸的双向调整 高频的输入要求单元尺寸很小 一定的单元尺寸对应输入的最大频率 一般进行滤波处理 滤掉低能量的高频 FFT.FIS
11、 Origin SeismoSignal,地震波的调整,El-Centro波,FFT,修正后的时程,修正后FFT,5Hz,动孔压的生成液化,干沙剪应变循环加载试验 初始加载阶段,沙土通常先压实再膨胀。卸载时,沙土遵循与加载相似的路径,但在零应变时,有些残余体积应变存在。取决于初始孔隙率,这可能代表纯粹的压实 假定孔隙中充满水 对于常体积测试,有效应力降低,孔隙水压保持不变 对于常荷载测试,(例如,盒子上法向荷载固定),孔隙水压增加,有效应力减小 有效应力为零时发生液化,动孔压的生成液化,因此孔隙水压增加不是液化的基本原因 由于颗粒间(重组以后)的低接触力导致有效应力的减小 描述液化的模型 高级
12、模型:BSHP (边界面低塑性本构模型, Wang et al. 1990) 简单模型:MC + 体积应变增量模型 Finn模型: Byrne模型:,主要内容,软件介绍 动力分析 Dynamic Option 桩-土相互作用分析 Interface 隧道分析 Structure Element 流固耦合分析 Fluid-Mechanical Interaction 学习方法及经验介绍,Interface单元简介,用途 岩体介质中的解理、断层、岩层面 地基与土体的接触 箱、槽及其内充填物的接触 空间中无变形的固定“障碍” 原理 三角形单元(无厚度!) 8参数 三种工作模式 粘结界面 粘接滑移 库
13、伦滑动,Interface的建立(1),关键要形成同一位置的两个节点(面) “移来移去”(推荐) 建两个分开的模型 建立接触单元 通过INI * add使模型接触 注意dist的含义 NO merge, NO attach!,接触面,dist,1,2,3,4,Interface的建立(2),“导来导去” 利用expgrid, impgrid命令进行网格导出与导入 配合DELETE命令 适于内部接触面的建立,或 其他前处理工具建立的网格,Interface的建立(2),save 1.sav del ran grop 2 not Interface 1 face save 2.sav rest 1
14、.savdel ran group 2expgrid 1.fac3d rest 2.savimpgrid 1.flac3d,Interface的建立(3),GEN separate INTERFACE wrap 指定正确的group,Interface参数的确定,虚构的为了合并节点而设置的接触面 kn=ks=10* 真实的刚性接触面 如料仓下料 c,D,Tension重要,kn,ks不重要 真实的柔性接触面 断层;水力劈裂材料 试验得到参数 对于kn,ks:岩石断层10100MPa/m(粘土); 100GPa(岩石) 反分析方法:通过断层中岩石的变形与原岩的变形,例3:单桩承载力分析,软土地基
15、 bulk 1.6878E6 shear 3.6167E5 coh 15E3 fric 12 dens 1.73E3 桩体 bulk 5e9 shear 3.75e9 dens 2.5e3,0.5m,8m,10m,20m,计算过程,施加桩顶荷载,计算结果,影响因素,计算工况设计,计算结果,Ks取1,Ks取2,Ks取3,最优方案,合理步骤,单桩分析,简单网格,接触面参数,多次试算,理论、实测,加密网格,接触参数,理想结果,群桩分析,Pile结构单元,单元参数,理想结果,Pile结构单元,单元参数,理想结果,主要内容,软件介绍 动力分析 Dynamic Option 桩-土相互作用分析 Inter
16、face 隧道分析 Structure Element 流固耦合分析 Fluid-Mechanical Interaction 学习方法及经验介绍,FLAC3D中的结构单元,有限单元 梁(beam)单元 锚索(cable)单元 桩(pile)单元 锚杆: rockbolt 壳(shell)单元 格栅(geogrid)单元 土工织物;土工格栅 初衬(liner)单元,beam,cable,pile,shell,geogrid,liner,结构单元的应用,土与结构的相互作用 桩基;基坑;边坡锚固 地下硐室的支撑结构;采矿;盾构 土工织物;土工合成材料 结构不宜复杂 岩土工程软件,不宜单纯的结构分析 复杂结构的模拟很困难 结构单元仍不完善 plot显示 双向接触结构(挡土墙) 结构单元的厚度,Liner结构单元,三节点扁平有限单元 每个节点有6个自由度 3个移动,3个旋转 能够抵抗膜及弯矩荷载 能够承受主方向的拉压应力 能够模拟管片与土体之间的分离及随后