《FLAC3D在岩土工程中的应用西安.ppt》由会员分享，可在线阅读，更多相关《FLAC3D在岩土工程中的应用西安.ppt（177页珍藏版）》请在三一办公上搜索。
1、,中国科学院计算所教育中心,徐帮树 系主任/副教授山东大学岩土与结构工程研究中心山东大学城市地下空间工程系2012年1月7-9日 西安,2,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,3,软件介绍,Fast Lagrangian Analysis of Continua美国Itasca咨询公司开发2D程序(1986)1990年代初引入中国有限差分法(FDM)DOS版2.0 2.1 3.0 3.1,软件介绍使用特征,命令驱动模式：界面操作/人机交互,专一性:专为岩土工程而开发本构模型 12个（1个空/3个弹性/8个塑性）
2、计算模式 5种（静力/动力/蠕变/渗流/温度）结构单元 6种（梁/锚杆/桩/壳/格栅/衬砌）界面单元 节理/断层/虚拟物理边界机理分析 绘图/时步函数,开放性：命令驱动模式/fish/C+,软件介绍计算特征,混合离散法 离散集成法,动态松弛法：完全动态运动方程在模拟物理上的不稳定过程不存在数值上的障碍,显式差分法:非线性/大变形/刚度矩阵,6,软件介绍应用范围,岩土工程/采矿工程 土木建筑/交通/地质/石油/环境/核电/水利工程,岩土体渐进破坏和崩坍现象断层结构影响和加固系统固结过程流变现象地下存储效果变形局部化剪切带演化动力稳定分析岩土体与结构相互作用分析振动液化现象,不足之处 求解时间受网
3、格尺寸的影响很大 某些模式下的计算求解时间很长 前处理功能较弱 可以克服！!,软件介绍,8,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,快速入门图形界面,快速入门分析的基本组成部分,生成网格单元,设置边界条件,定义材料性质,设置初始条件,初始地应力平衡,加载及连续建模,求解,输出计算结果,建立分析模型,模拟求解,快速入门简单分析命令概要,快速入门文件类型,.dat 命令文件 call 记事本.fis 二次开发 call 记事本.tmp 临时文件 自动消失.sav 结果文件 restore.log 日志文件 记事本.fl
4、ac3d 网格文件 Impgrid 记事本,快速入门成果输出,图片输出,记录结果输出hist keyword x y zhist keyword id=nhist write nhist1 记事本编辑,14,快速入门一个最简单的例子,gen zon bri size 3 3 3;建立网格(前处理)model elas;材料参数prop bulk 3e6 shear 1e6ini dens 2000;初始条件fix z ran z-.1.1;边界条件fix x ran x-.1.1fix x ran x 2.9 3.1fix y ran y-.1.1fix y ran y 2.9 3.1set
5、grav 0 0-10solve;求解app nstr-10e4 ran z 3 x 1 2 y 1 2solveplo con zd;后处理切片功能,RUN FLAC3D,例2.1,快速入门收敛标准和变量解释,常用标准 不平衡力比10exp(-5)/不平衡力/典型内力自定义标准 不平衡力网格节点速度塑性区标识 shear-n/shear-p tension-n/tension-p历时曲线,例2.1,16,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,17,基本原理,有限差分法Lagrangian网格空间混合离散技术La
6、grangian格式动量平衡方程FLAC3D的求解过程FLAC3D的本构模型,18,基本原理有限差分法,古老的方法(上世纪40年代)用差分格式转化控制方程中的微商格式流体力学；土工渗流问题；固结FDM&FEM的混合求解FDM的新进展,19,基本原理Lagrangian网格,源自流体力学中的拉格朗日法跟踪流体质点的运动状态跟踪固体力学中结点，按时步用Lagrangian法研究网格节点的运动节点和单元随材料移动，边界和接触面与单元的边缘一致固体力学大变形理论,法国数学家、物理学家拉格朗日,20,基本原理空间混合离散技术,+,/2,=,21,基本原理Lagrangian格式动量平衡方程,F(t),m
7、,牛顿运动定律,对于连续体,在静力平衡条件下，加速度项为0，方程变为平衡方程,22,基本原理FLAC3D的求解过程,23,基本原理FLAC3D中的本构模型,开挖模型null 3个弹性模型各向同性弹性横观各向同性弹性正交各向同性弹性 8个弹塑性模型Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型、双屈服面塑性模型、修正剑桥模型和胡克布朗模型,24,基本原理FLAC3D中的本构模型,25,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,26,前后
8、处理,基本前后处理命令操作菜单操作外界模型的导入复杂模型的网格检查,27,菜单驱动(前处理/计算模式),命令栏,28,FLAC3D的前后处理,命令驱动(推荐)程序控制图形界面接口计算模型输出指定本构模型及参数指定初始条件及边界条件，指定结构单元指定接触面指定自定义变量及函数(FISH)求解过程的变量跟踪进行求解模型输出,前后处理功能的优点,13种zone类型后处理快捷、方便、丰富计算过程中的hist变量动态显示FISH可进行参数化模型设计单元状态的可编程计算暂停时的后处理与可保存,前后处理功能的缺点,复杂模型的建模功能不强可以编程导入其他软件形成的网格(比如：Ansys、Adina、ABQUS
9、、GeoCAD)无等值线的后处理功能(3D)可编程将.sav文件写入TecPlot/Surfer等其他后处理软件全命令操作，学习困难鼠标功能单一(双击取击点坐标),FLAC3D 网格模型,FLAC3D 网格模型,FLAC3D的建模功能不是很强大，一般用ANSYS/ABQUS建模，然后导入FLAC3D，可以通过编写接口程序实现。,FLAC3D Grid Generation with ANSYS+AutoCAD,33,外界模型的导入,复杂网格的生成难度大接口编写不同软件之间的精度差异Group的定义采用.flac3d的文件格式.flac3d文件的格式G 1 1.0e+00 1.0e+00 1.0
10、e+00Z B8 1 2 3 4 5 6 7 8ZGROUP Soil1 2 3impgrid&expgrid仅限于网格,FLAC3D 网格单元数据形式 L5.1,例5.3,35,外界模型导入实例,AutoCAD辅助建模Ansys建模Ansys模型数据导出Ansys_command-ywd.dat/APDLFlieRead input from生成node-01.dat和elem-02.dat数据格式转换myansytoflac2.exe得到flaczone.FLAC3D得到FLAC3D模型FlieImport grid flaczone.FLAC3DSav geo.sav,后处理查看分析结果
11、,云图、矢量图、曲线、数据和动画等值线图 TecPlot/SurferPlot 命令格式/快捷键分组：PLOT block group/界面方式出色背景：Ctrl+G/界面方式快捷键：缩放/旋转/平移图形输出到文档（word）图形窗口：【Edit】/【Copy to clipboard】（建议使用：矢量图）命令窗口：【File】/【Print Type】/【Jpg File】【File】/【Print Setup】/【Jpg File】【File】/【Print】/【1Base 0】,例6.1,Set plot jpgSet plo quality 100Pot hard file 1.jpg
12、,后处理查看分析结果,后处理查看分析结果,初始地应力结果后处理 Plo con/bcon应力云图PLOT con szz/界面方式 Plo con/bcon,PLOT con szz ou on;显示网格PLOT con szz ef on;有效应力PLOT con szz inter 1e4;改变云图显示增量PLOT con szz max-10e3;改变应力最大值,PLOT bcon szz PLOT con 节点/单元（插值）PLOT bcon 单元（检查）PLOT bcon prop bu,例6.1,后处理施加荷载后的计算结果,变形云图：plo con zd/szz 动态显示 plo
13、sk magf 20 变形矢量图：plo sk dis 塑性区分布：plot block state print zone state 变量监测:plot hist 2 根据坐标指定/根据ID号指定 set hist_rep 1 监测数据导出 hist write 7 vs 8 file 6-3hist.txt结果输出:print apply,fish,gp,group,history,model,zone等？set log on/set logfile 6-1.log切片操作(插值)：plo set plane ori 0 1.5 0 norm 0 1 0 plo con zd plane/
14、plo add ske,dis plane,axe动画制作,例6.2 6.3 6.4 6.5,40,Tecplot的slice切片功能.和CAD一样,可以任意切剖面出图.最大优点是可以几个剖面同时出图.,Tecplot的切片功能,41,Sufer的后处理,42,复杂网格的检查,FLAC3D生成的复杂网格attach face无接触面时检查整体模型是否存在sub-grid有接触面时给定范围进行检查其他软件导入的复杂模型网格划分的检查gen merge弹性模型model elastic求解“独立”节点“畸形”单元,例6.1,43,模型的检查,FLAC3D本身的Check功能十分有限错误提示很少十分
15、开放的工作平台 检查的基本步骤网格检查(如前所述)边界条件检查速度约束条件 plo gpfix red sk荷载条件 plo fap red sk模型检查模型赋值plo block model参数赋值plo block prop*初始应力检查(如后所述),例6.1,44,初始应力的生成,为什么要单独列出？得到所关注分析阶段之前岩土体已经存在的应力状态分析过程中出现的很多问题都与初始应力是否合理有关手册中的例子五花八门生成方法弹性求解更改强度参数的弹塑性求解设置初始应力的弹塑性求解存在水压力的初始应力生成水下建筑的初始应力生成,45,弹性求解,gen zon bri size 1 1 2m el
16、asprop bulk 3e7 shear 1e7fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0-10solve,Step=162sz=-40e3sx=-21.54e3,例7.1,浅埋和地表工程无屈服区域,46,更改强度参数的弹塑性求解,gen zon bri size 1 1 2model mohrprop bulk 3e7 shear 1e7 c 1e10 f 15 ten 1e10fix z ran z 0fix x ran x 0fix x ran x
17、1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0-10solveprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0solve,Step=163sz=-40e3sx=-21.54e3,例7.2,7,4,可能存在屈服区域计算时间长,47,设置初始应力的弹塑性求解,gen zon bri size 1 1 2model mohrprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0f
18、ix y ran y 1ini dens 2000ini szz-40e3 grad 0 0 20e3 ran z 0 2ini syy-20e3 grad 0 0 10e3 ran z 0 2ini sxx-20e3 grad 0 0 10e3 ran z 0 2set grav 0 0-10solve,Step=0sz=-40e3sx=-20e3,例7.4,加速平衡时间K0不适用深埋工程,48,存在静水压力的初始应力生成,gen zon bri size 1 1 2model mprop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10fix z ran z
19、 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 1ini dens 1500 ran z 1 2ini szz-35e3 grad 0 0 20e3 ran z 0 1ini syy-17.5e3 grad 0 0 10e3 ran z 0 1ini sxx-17.5e3 grad 0 0 10e3 ran z 0 1ini szz-15e3 grad 0 0 15e3 ran z 1 2ini syy-7.5e3 grad 0 0 7.5e3 ran z 1 2ini sxx-7.5e3
20、grad 0 0 7.5e3 ran z 1 2ini pp 10e3 grad 0 0-10e3 ran z 0 1set grav 0 0-10,Step=83sz=-33.75e3sx=-23.21e3,rd=rs nsrw,例7.5,Apply/IniP108,49,水下建筑的初始应力生成,gen zon bri size 1 1 2model mprop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 20
21、00 ran z 0 2ini szz-50e3 grad 0 0 20e3 ran z 0 2ini syy-30e3 grad 0 0 10e3 ran z 0 2ini sxx-30e3 grad 0 0 10e3 ran z 0 2ini pp 30e3 grad 0 0-10e3 ran z 0 2app nstress-10e3 ran z 2set grav 0 0-10solve,Step=83sz=-50e3sx=-32.1e3,例7.6,顶部静水压力,50,深埋工程的初始应力生成,例7.7，7.8,构造应力场和自重应力场的叠加无真正的位移边界条件和力边界条件模型内应力重分布
22、是应力与自重应力相平衡的结果，初始应力场即自重应力场李仲奎/戴荣（2002）提出快速应力边界法（S-B法）,51,初始应力检查办法,初始应力的计算时间不会“非常长”经常检查模型的响应plo con szz(syy,sxx)应力场plo con zdis(ydis,xdis)位移场plo blo sta屈服状态plo gpfix red sk速度约束条件plo fap red sk体力plo hist(unbal)不平衡力plo interface nstress(sstress)接触面单元,学习经验和建议,了解FLAC3D的使用范围、优点和局限性由简到繁，循序渐进充分利用手册了解计算中每条语句
23、的含义多做数值试验夯实知识基础相互交流，取长补短,53,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,54,FISH语言简介,软件自带的编程语言 是否一定要学？视情况而定，需要时查询FISH变量即可语法简单（变量/函数/运算符/浮点数/字符串等）xxxend_xxx注意事项语法严谨，保留字不可缩写变量可不定义，因此注意检查程序print fishtable,extra等命令使用文本文件，后缀.txt.fis,例8.1,Fish基本概念,函数和变量二者都要在fish中赋值：与常规语言类似，按运算符优先级先后顺序执行数据类型
24、规则：类型对应命名规则：不能以数字开头，不含中文，且不能包含下列符号 变量/函数名不能与保留字相冲突：不用简单的单词（P4）不能递归调用作用全局：避免在不同的函数中有相同的变量命令：PRINT/HISTORY/SET数据类型：整型/浮点型/字符型（以为分界符）/指针型,Fish主要语句,选择 Caseof 表达式.case n1.case n2.ENDCASE,条件 IF 条件表达式 then.else.ENDIF 无“并”“或”“否”,嵌套条件语句用于描述分段函数,例8.2,Fish主要语句,循环 变量数值LOOP var(exp1,exp2)ENDLOOP条件表达式LOOP while 表
25、达式ENDLOOP语法要求高关键词不能缩写,例8.3,命令COMMAND.(FLAC3D命令)ENDCOMMAMD,内置变量与函数,FLAC3D 标量：流变时间、节点总数、最大不平衡力常规标量：时钟变量、随机数、空指针、常数FLAC3D 模型变量：内存、节点、单元、接触面（P33）结构单元变量：常规、结构、节点、连接数学函数：绝对值、三角、对数表格函数和内存函数绘图函数,单元遍历和节点遍历,遍历？信息/赋值基本变量：zone_head/gp_head z_next()/gp_next(),单元遍历程序：p_z=zone_headloop while p_znull 语句 p_z=z_next(
26、p_z)endloop,节点遍历程序：p_gp=gp_headloop while p_gpnull 语句 p_gp=gp_next(p_gp)endloop,60,FISH应用,例8.7,获得最大节点位移,土体的模量随小 主应力变化,E=KPa(3/Pa),n,例8.5,分级加载与监测,例8.6,得到主应力差值云图,例8.8,编程与差错技巧,逐级编写变量函数命名规则变量无大小写区别同一类变量名用相同的字母开头字符类型及时修正语法严格 缩写 简写 标点符号PRINT fish 检查,62,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自
27、定义本构模型,63,接触面单元,应用范围原理建模方法推荐方法复杂内部接触面的设置方法参数选择单桩承载力分析挡土墙的接触面设置思路：未知问题的分析方法,64,接触面单元的用途和原则,岩体介质中的解理、断层、岩层面地基与土体的接触箱、槽及其内充填物的接触空间中无变形的固定“障碍”,大小表面相连，接触面应该建在小的表面上相邻网格密度不同，应该建在密度大的区域上接触面单元的尺寸通常应该等于或小于相连的目标面尺寸用Attach命令连接的两个表面不应再建立接触面,65,接触面的原理,三角形单元(无厚度!)库伦剪切模型，参数较多单面的，通过接触面节点与目标面建立联系三种工作模式粘结界面粘接滑移库伦滑动,66
28、,接触单元模型的建立(1),关键要形成同一位置的两个节点(面)，模拟结构物的滑移，分离“移来移去”(推荐)建两个分开的模型建立接触单元通过INI*add使模型接触注意dist的含义NO merge,NO attach!,接触面,dist,1,2,3,4,1、建立土体网格,例9.1,2、建立节理单元,3、建立柱体单元,4、桩体向下平移,70,接触单元模型的建立(2),“导来导去”利用expgrid,impgrid命令进行网格导出与导入配合DELETE命令适于内部接触面的建立，或其他前处理工具建立的网格,例9.2,72,接触单元模型的建立(3),GEN separate INTERFACE wra
29、p指定正确的group,例9.3,接触单元的参数,;切割模型法;-ngen zone brick size 3 3 3group 1 range x 1 2 y 1 2 z 2 3group 2 range group 1 notgen separate 1 int 1 wrap 1 2int 1 maxedge 0.5plo int red,74,接触面参数的确定,虚构的为了合并节点而设置的接触面Kn=ks=10*max(K+4/3G)/Dzmin真实的刚性接触面如料仓下料c,D,Tension重要，kn,ks不重要真实的柔性接触面断层；水力劈裂材料试验得到参数对于kn,ks：岩石断层101
30、00MPa/m(粘土);100GPa(岩石)反分析方法：通过断层中岩石的变形与原岩的变形,谷仓下料模拟,手册中的一个例子(ftd132.pdf)，设置接触面上粘聚力c=0，在自重作用下下滑。,例9.4,滑动前,滑动过程中,76,单桩承载力分析,软土地基bulk 1.6878E6 shear 3.6167E5 coh 15E3 fric 12dens 1.73E3桩体bulk 5e9shear 3.75e9dens 2.5e3,0.5m,8m,10m,20m,例9.5,77,计算过程,施加桩顶荷载,计算结果,78,挡土墙的接触面设置,对于未知问题的分析思路3个独立的接触面3个同ID的接触面2个独
31、立的接触面2个独立的接触面并进行底部merge,挡墙,土体,1,2,3,Wall,Soil,79,挡土墙的接触面设置(3),前两种方法的差别实质3个ID的独立接触面在相同位置产生互不影响的两个节点共同ID的接触面在相同位置自动设置为1个节点最终的结论需要您自己去判断！,3 interfaces,2 IDs,1 interface,1 ID,80,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,81,FLAC3D中的结构单元,有限单元梁(beam)单元锚索(cable)单元桩(pile)单元锚杆:rockbolt壳(shel
32、l)单元格栅(geogrid)单元土工织物；土工格栅初衬(liner)单元,beam,cable,pile,shell,geogrid,liner,82,结构单元的应用,土与结构的相互作用桩基；基坑；边坡锚固地下硐室的支撑结构；采矿；盾构土工织物；土工合成材料结构不宜复杂岩土工程软件，不宜单纯的结构分析复杂结构的模拟很困难结构单元仍不完善plot显示双向接触结构(挡土墙)FLAC3D 3.1结构单元的厚度,梁单元(beam),梁单元局部坐标系及自由度,梁单元力和力矩（包括轴力、剪力、弯矩）,梁单元(beam),本构模型：各向同性的、无屈服破坏的线弹性材料材料参数：,密度弹性模量泊松比朔性矩热膨
33、胀系统横截面积y轴的惯习矩z轴的惯性矩极惯性矩矢量Y,梁单元的建模,两种方式建立结构单元：直接通过坐标点指壑结型结构的两个端点选建节点，然后通过节点建立结构单元,例10.1,例10.2,单元id,构件cid,节点node,连接 link,梁单元(beam),案例：使用梁单元进行开挖支护问题描述：分析范围的尺寸为6m8m8m，开挖步分的尺寸为2m4m5m。开挖时采用梁单元进行支撑，分析采用支撑后模型的变形情况和梁单元的受力情况。,例10.4,计算网格,位移云图、梁轴力图,Beam单元案例,培训案例彭文斌教材ch11,荷载,位移矢量图,剪力图,锚索单元(cable),梁单元局部坐标系及自由度,锚索
34、单元本构关系,锚索单元(cable),锚索、注浆体、围岩三种材料相互作用原理,锚索单元(cable),本构模型：锚索构件为弹塑性材料，拉、压屈服，不能抵抗弯矩。材料参数：,密度弹性模量单位长度上水泥浆的粘结力水沁浆的摩擦角单位长度上水泥浆刚度水泥浆上圈周长大变形滑动标志大变形滑动容差热膨胀系统横截面积抗压强度抗拉强度,几何参数、材料参数、水泥浆特性，12个参数；10个必选参数，2个可选参数,锚索单元(cable),案例：关于预应力锚杆的模拟问题描述：对于非全长锚固、并施加预紧力的锚索，需要考虑托盘、自由段和锚固段、和预应力施加的问题。方法一：删除-建立link连接来模拟托盘方法二：通过设置极锚
35、固剂参数模拟托盘方法三：借助别的结构单元（如liner单元）来模拟托盘,例10.5,例10.6,例10.7,壳型结构单元的建模,两种方式：1、在实体单元表面建立节点、构件和单元，然后删除实体单元。2、先建立节点，然后用sel shellSel命令壳型结构的构件。,壳型结构单元的建模,壳型结构单元,三种单元形式：壳单元、格栅单元、初砌单元壳构件：假定为是均厚的3节点组成的三角形。由构件组成任意形状的壳。力学特性：壳材料结构向应、构件与网格的交互。定义为各向同性或异性的线弹性材料，无破坏极限。,壳型结构单元,密度各向同性（弹模、泊松比）；正交各向异性热膨胀系数厚度,壳构件与网格是刚性连接的,壳单元
36、(shell),彭11.10,案例：松软弹性土中开挖隧道，原始应力为1MPa，喷射厚0.2m的混凝土。混凝土保持弹性且与土壤刚性连接。,初始状态,开挖第1步,加喷射混凝土衬砌，开挖第2步,壳单元(shell),彭11.10,加衬砌，最大位移1mm,无衬砌，最大位移12mm,98,Liner结构单元,三节点扁平有限单元每个节点有6个自由度3个移动，3个旋转能够抵抗膜及弯矩荷载能够承受主方向的拉压应力能够模拟管片与土体之间的分离及随后的重新接触能够模拟管片与土体之间的摩擦相互作用,法向,切向,99,隧道与土体的相互作用,半圆隧道直径3.25m上覆土层厚度5m计算范围3r土体弹性计算(K=30MPa
37、,G=10MPa)参数化编程几何尺寸模型参数网格形状,ht,hb,r,B,100,计算步骤,模型网格,初始应力生成,施加管片,计算结果,101,管片的连接,冷连接弯矩和剪力不能直接在环与环间传递，只能通过其相邻的介质传递全连接相邻的Liner单元在连接处共用一个节点，连接处重叠单元不能发生移动或旋转 结点连接即结点间的连接在6个方向的自由度上用弹簧来模拟，每个自由度都可具有一定的特性,通缝拼接,错缝拼接,102,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,103,流-固耦合分析(单相流),基本功能理论框架计算模式（渗流
38、模式/无渗流模式）渗流边界条件，初始条件单渗流计算及渗流耦合计算,104,基本功能,不同的渗流模型和属性流体压力，涌入量，渗漏量和不渗水边界抽水井、点源、体积源饱和渗流可采用显式差分法、隐式差分法 非饱和渗流采用显式差分法渗流-固体-热的耦合流体和固体的耦合程度依赖于土体颗粒(骨架)的压缩程度，用Biot系数表示颗粒的可压缩程度。,105,计算模式,无渗流模式不设置CONFIG Fluid孔压不改变设置孔压分布INITIAL ppWATER tableWATER densityWATER table faceSET gravity手动设置干湿密度,渗流模式设置CONFIG fluid设置土体干
39、密度 渗流模型MODEL fl_isotropic MODEL fl_anisotropic MODEL fl_null 开挖单元设置WATER dens/INI fdens,12.1，12.2,无渗流模式,没有设置config fluid。可以在节点上设置孔隙水压力，改变土体有效应力状态。两种方式设置孔隙水压力ini pp,water table,例12.1,gen zon bri size 3 3 3 ini pp 20e3 grad 0 0-10e3plot con pp outline on,gen zon bri size 3 3 3 water density 1000set gr
40、av 10water table ori 0 0 3 norm 0 0 1,例12.2,无渗流模型模式，Flac3D在体积力计算时，不能自动考虑流体质量作用；手动设置水位面以上为干密度；水位面以下不饱和密度,流体参数赋值方法,参数和单位,渗透系数密度比奥系数与比奥模量（土颗粒可压缩）体积模量/孔隙率/饱和度/抗拉强度,Set fluid biot on,110,渗流边界条件,默认的边界条件是不透水边界 孔隙压力自由(不透水边界)固定孔隙水压力(透水边界)命令Fix pp(配合ini pp使用)Apply pp节点的渗流边界条件 APPLY pwell 流量interior单元的渗流边界条件AP
41、PLY vwell 流量平面的渗流边界条件 通量边界：APPLY discharge 通量（m/sec）渗漏边界：Apply leakage v1 v2(v1为孔压；v2表示渗漏系数，m3/N.sec),例12.3,12.4,111,单渗流计算及渗流耦合计算,完全耦合分析方法 两个进程：流体计算，力学计算孔压固定分析(有效应力分析)比如边坡稳定性，只考虑静水压力分布。Ini pp;water table单渗流得到孔压分布 Config fluidSet mech off流体本构模型、参数step;solve age time(只有当流体模量为真实值时，)无渗流计算孔压的力学响应 流-固耦合计算
42、,112,单渗流得到孔压分布,固定孔压分析用途：排水沟；抽水井；隧道渗流计算步骤CONFIG fluid SET mech offSET fluid implicit on/off（隐式仅完全饱和）MODEL fl_;PROP STEP;SOLVE age;SET fluid ratioSET fluid off mech onPROP biot_c 0(or set fmod 0)（避免应力变化引起孔压变化）,113,无渗流计算孔压的力学响应,不排水短期响应两种分析方法：干法和湿法干法：Ku=K+a2M两种破坏形式WATER或INI获得常孔压，不排水的c,(孔压改变较小)=0,c=cu(MK
43、+4/3G)湿法：耦合体系的短期行为（必须渗流模式）使用排水的K,c,若SET fluid off,Biot_mod(fmod)真实,例12.5,实例：载荷引起的地基土体的起孔隙水压力地基土的宽度为20m，高度为10m，采用弹性模型。地基土表面为透水边界条件，孔压固定为0。地基土表面3m的范围内缓慢施加40kPa，载荷。计算土体超孔隙水压力,孔压分布图、节点荷载plot con pp plane ou onplot add fap red plane,垂直位移、位移矢量图Plot con zdis out on planeplot add dis plane,例12.5,115,流-固耦合计算
44、,CONFIG fluid;M(Kf);K(渗透系数)真实，则FLAC3D默认耦合计算pevsevp预估流/力特征时间，力学扰动可认为瞬时发生耦合计算前先达到一个平衡状态SET fluid on mech off;SET fluid off mech on;STEPSET mech force;SET mech substep n auto;SET fluid substep m(=1)STEP:渗流步足够小,116,流固耦合的计算方法,手动调整的STEP求解SET fluid on mech offSTEP n1SET fluid off mech onSTEP n2主从进程的SOLVE求解
45、SET mech forceSET mech substep n auto(从进程)SET fluid substep m(主进程)SOLVE age自动STEP求解STEP/slove,117,渗流问题(CONFIG fluid)分析步骤,时间比例(ts,tc)稳态不排水状态相当扰动类型力学扰动孔压扰动流固刚度比Rk是否1完全耦合模式时间比例相当；力学扰动,118,心墙土坝的渗流(1),网格模型,初始孔压,例12.6,问题描述：土坝中存在一道粘土心墙，心墙渗透系数远底于土体的渗透系数，可认为是不透水材料。分三步进行：第1步，水位未上升时的应力场、孔压场,119,心墙土坝的渗流(2),第二步：
46、水位上升对模型应力影响。渗流场关闭,竖向应力,沉降,例12.7,120,心墙土坝的渗流(3),第三步：施加孔压边界，升高迎水面的孔压值，并在背水面设置透水的孔压边界条件。,例12.8,孔隙水压力、矢量图,121,真空预压的简单模拟,孔压边界条件tstc长期分析(排水)Rk1骨架很软孔压扰动进行biot_mod调整,砂层,软土层,粘土层,PVD,2m,8m,10m,例12.8,真空预压的简单模拟,粘土层表面沉降，约5cm,地表沉降，约30cm,沉降历时曲线,孔压云图,123,主要内容,软件介绍快速入门基本原理静力计算FISH语言接触面单元结构单元及应用流固耦合分析非线性动力分析自定义本构模型,概
47、述,Config dynamic与其他场的耦合结构单元耦合流体计算耦合热力学计算耦合大变形模式外部和内部荷载外部：地震内部：爆破脉冲荷载；落锤冲击荷载；交通荷载等,125,完全非线性的动力分析,动力荷载动力边界条件力学阻尼与滞回阻尼地震波的调整,动力载荷,模型边界或内部节点施加动载荷：加速度时程；速度里程；应力（压力）时程；集中力时程,加速度时程曲线,作用在炮孔壁上的爆破脉冲荷载时程曲线,动力荷载,两种加载方式：fish函数；table表。,def xxxxxx=.dytimeendapp xvel=1.0 hist xxx range,fish函数,y列均匀间隔的表,第1列：表的名称第2列：
48、数据对的个数 空格 时间间隔第3列：y列的第1个数据第4列：y列的第2个数据空行,分别给出x,y数据对的表,第1列：表的名称第2列：x1 空格 y1第3列：x2 空格 y2空行,动力时间步,临界时间步不建议放大动态多步 set dyn multi on材料密度必须设置Fish函数和变量名必须一致,案例：动态时步问题描述：土体的深度为10m，挡土墙的高度为5m。,例11.1,130,Quiet边界,静态(quiet,粘性)边界Lysmer and Kuhlemeyer(1969)模型边界法向和切向设置独立的阻尼器性能对于法向p波和s波能很好的吸收对于倾斜入射的波和Rayleigh波也有所吸收，但
49、存在反射人工边界仍应当足够远,131,Quiet边界应用,内部振动(如隧道中的列车振动问题/爆破)动力荷载直接施加在节点上使用Quiet边界减小人工边界上的反射不需要FF边界外部荷载的底部边界软土地基上的地震荷载不适合用加速度或速度边界条件使用应力条件t=-2Csrvs地震底部输入的侧向边界扭曲了入射波,quiet,quiet,quiet,案例：静态边界的例子问题描述：一根竖直的弹性杆，高50，宽1。杆的底部设置静态边界条件，杆的顶部为自由表面，杆的底部施工水平方向的应力冲击载荷。,例11.2,桩的动力学响应,134,Free-field边界,Cundall et al.(1980)自由场网格
50、与主体网格的耦合粘性阻尼器，自由场网格的不平衡力施加到主体网格边界上设置条件底部水平，重力方向为z向侧面垂直，法向分别为x,y向其他边界条件在APPLY ff之前,去掉模型底部静力约束条件施加静态边界条件施加动力荷载施加自由场边界条件,例11.3,135,力学阻尼,瑞利(rayleigh)阻尼假设阻尼与质量、刚度的线性关系参数确定简单(等价平均应变=60%*emax)中心频率(共振计算，地震平均频率)临界阻尼比（0.5%）计算速度慢局部(local)阻尼FLAC3D的静力分析阻尼参数简单适合简单情况/高频噪声,例11.4,136,滞回阻尼(Hysteretic Damping),模拟岩土介质的
