FLACFLAC3D基础与应用全部陈育民.pptx

GeoHohaiFLAC/FLAC3D基础与应用基础与应用陈育民陈育民 副教授副教授河海大学土木与交通学院河海大学土木与交通学院Email:1河海大学，河海大学，20112011年年1111月月2424日日GeoHohai几个问题什么是什么是FLAC？为什么要用为什么要用FLAC？FLAC能做什么？能做什么？FLAC为何这么流行？为何这么流行？怎么学怎么学FLAC？GeoHohai什么是FLAC？Fast Lagrangian Analysis of Continua3连续介质拉格朗日算法快速GeoHohai为什么要用FLAC？4Abaqus：745条Plaxis：80条ADINA：310条GeoHohaiFLAC能做什么？岩土工程中的绝大多数问题岩土工程中的绝大多数问题土力学、岩石力学、防灾减灾、隧道、地下空间等土力学、岩石力学、防灾减灾、隧道、地下空间等采矿工程中的大部分问题采矿工程中的大部分问题水工结构中的部分问题水工结构中的部分问题结构工程结构工程国国际通用的岩土工程通用的岩土工程专业分析程序分析程序5GeoHohaiFLAC为何这么流行？Charles Fairhurst美国工程院、瑞典皇家工程院美国工程院、瑞典皇家工程院院士，国际岩石力学学科院士，国际岩石力学学科和岩石力学学会创始人之一，历任国际岩石力学学会主和岩石力学学会创始人之一，历任国际岩石力学学会主席和副主席，国际岩石力学学会席和副主席，国际岩石力学学会MULLER奖、美国岩石力奖、美国岩石力学学会终生成就奖获得者。学学会终生成就奖获得者。Peter Cundall美国工程院、英国皇家工程院院士，国际资深计算岩石美国工程院、英国皇家工程院院士，国际资深计算岩石力学学家。力学学家。6GeoHohai怎么学习FLAC？7问：好学不？答：好学！问：难不？答：难！GeoHohai报告安排第一讲：基本介绍、静力分析、前后处理第一讲：基本介绍、静力分析、前后处理第二讲：接触面、第二讲：接触面、FISH语言、流固耦合分析语言、流固耦合分析第三讲：动力分析、自定义本构、结构单元第三讲：动力分析、自定义本构、结构单元第四讲：第四讲：FLAC（2D）基本介绍与应用实例）基本介绍与应用实例讨论讨论8GeoHohai第一讲FLAC3D基本介绍、静力分析、前基本介绍、静力分析、前后处理后处理9GeoHohai软件介绍Fast Lagrangian Analysis of Continua美国美国Itasca咨询公司开发咨询公司开发2D程序程序(1986)1990年代初引入中国年代初引入中国有限差分法有限差分法(FDM)3D版本：版本：DOS版版2.0 2.1 3.0 3.14.02D版本：版本：DOS版版4.05.06.07.010GeoHohaiFLAC1.大应变大应变、小应变计算模式、小应变计算模式.2.丰富的本构模型、提供自定义的本构模型功能丰富的本构模型、提供自定义的本构模型功能3.接触面可以模拟不同材料的接触接触面可以模拟不同材料的接触4.流固耦合实现土体的固结与渗流流固耦合实现土体的固结与渗流5.拥有各种功能的结构单元类型，拥有各种功能的结构单元类型，模拟土与结构的相互作用模拟土与结构的相互作用6.强大的强大的动力分析功能动力分析功能.7.流变分析，拥有粘弹性模型和粘塑性模型流变分析，拥有粘弹性模型和粘塑性模型8.热力学分析热力学分析.11ShearstrainratecontoursGeoHohaiFLAC3D12-与FLAC类似，是FLAC的三维版本-与FLAC拥有相同的优点upstreamdownstreamGeoHohai基本特点内置材料模型内置材料模型连续介质非线性，大应变模拟连续介质非线性，大应变模拟显式解题方案，显式解题方案，为不稳定物理过程提供稳定解为不稳定物理过程提供稳定解界面或滑动面用来模拟可产生滑动或分离的离散面，界面或滑动面用来模拟可产生滑动或分离的离散面，从而模拟断层，节理或摩擦边界从而模拟断层，节理或摩擦边界内置材料模型丰富：内置材料模型丰富：n零模型零模型,n三个弹性模型三个弹性模型(各向同性，横观各向同性和正交各向同性，横观各向同性和正交各向异性各向异性),n八个八个朔性朔性模型模型(德鲁克德鲁克-布拉格布拉格,摩尔摩尔-库伦库伦,应变硬应变硬化化/软化，单一节理，双线性应变硬化软化，单一节理，双线性应变硬化/软化单一软化单一节理节理,双屈服，修正剑桥粘土，霍克双屈服，修正剑桥粘土，霍克-布朗布朗)13隧道工程GeoHohai可选模块可选模块包括可选模块包括:n 热力学热力学,热热-力学耦合力学耦合,热热-流体流体-力力学耦合包括热传导和对流学耦合包括热传导和对流;n 粘弹粘弹,粘朔性粘朔性(蠕变蠕变)材料模型材料模型;n动力学分析动力学分析,并可以模拟静边界和并可以模拟静边界和自由域自由域 n使用使用C+定义自己的模型定义自己的模型14核废料储存中的热力学研究问题核废料储存中的热力学研究问题GeoHohaiFLAC/FLAC3D基本原理基本原理FLAC/FLAC3D利用有限差分，显示方案，动态利用有限差分，显示方案，动态松弛方法模拟连续体的非线性力学行为：松弛方法模拟连续体的非线性力学行为：n即使对准静态问题，程序仍然求解完整的动力学方程。即使对准静态问题，程序仍然求解完整的动力学方程。这种方法的好处在于可以为物理非稳定过程例如塌方提这种方法的好处在于可以为物理非稳定过程例如塌方提供稳定解；供稳定解；n在在 “松弛松弛”方法中，使用阻尼来吸收动能以模拟系统方法中，使用阻尼来吸收动能以模拟系统的的“静态静态”反应。反应。这种方法可以用比其它方案如解矩这种方法可以用比其它方案如解矩阵法更为真实有效地模拟塌方问题。阵法更为真实有效地模拟塌方问题。15GeoHohaiLagrangian法源自流体力学中的拉格朗日法源自流体力学中的拉格朗日法跟踪流体质点的运动状态跟踪流体质点的运动状态跟踪固体力学中结点，按时步用跟踪固体力学中结点，按时步用Lagrangian法研究网格法研究网格节点的运动节点的运动节点和单元随材料移动，边界和接触面与单元的节点和单元随材料移动，边界和接触面与单元的边缘一致边缘一致固体力学大变形理论固体力学大变形理论16法国数学家、物理学家拉格朗日 GeoHohai混和离散技术混和离散技术FLAC混和离散混和离散17+/2=每个每个 为常应力为常应力/应变应变:体积应变由整个四边形算出体积应变由整个四边形算出.应变偏量则有两个三角形应变偏量则有两个三角形 和和 分别算出分别算出(混合离散混合离散 过程过程)解题过程中网格坐标按照解题过程中网格坐标按照“拉格朗日方式更新拉格朗日方式更新”(网格随材料移动网格随材料移动),且为且为显式显式(一个时步内局部变化不会影响邻域一个时步内局部变化不会影响邻域)GeoHohaiFLAC3D混和离散混和离散18+/2=GeoHohaiFLAC3D混和离散19n结构域离散为可由四面体单元组合形成的五面体或六面体等单元;n以 为基本单元(常应力、常应变);n体应变的计算：;n偏应变的计算：.GeoHohai动态松弛动态松弛法动态松弛法 在动态松弛法中，网格点根据牛顿运动定律运动在动态松弛法中，网格点根据牛顿运动定律运动.网格点的速度与该网格点的速度与该点的不平衡力呈正比点的不平衡力呈正比.这种求解方法所决定的一系列位移将把系统这种求解方法所决定的一系列位移将把系统带入平衡状态，或表明破坏模式带入平衡状态，或表明破坏模式.在动态松弛法中有两个因素很重要在动态松弛法中有两个因素很重要:1.1.时步的选择时步的选择2.2.阻尼效应阻尼效应 20GeoHohai显式算法显式解与隐式解的比较显式解与隐式解的比较21显式显式,逐时推进逐时推进隐式隐式,静态静态1.无需进行反复迭代来实现非线无需进行反复迭代来实现非线性本构关系性本构关系.2.类似问题求解时间呈类似问题求解时间呈 N3/2 规规律增长律增长3.物理非稳定性不会引起数值物理非稳定性不会引起数值不稳定性不稳定性.4.因为无需储存矩阵，用较小因为无需储存矩阵，用较小内存即可模拟大尺度问题内存即可模拟大尺度问题.5.对大位移、大应变问题同样适对大位移、大应变问题同样适合，无需额外的计算合，无需额外的计算.1.需进行反复迭代来实现非线性需进行反复迭代来实现非线性本构关系本构关系 2.类似问题求解时间呈类似问题求解时间呈 N2 甚至甚至 N3规律增长规律增长.3.难以模拟物理非稳定性问题难以模拟物理非稳定性问题.4.需存储刚度矩阵，需克服相关需存储刚度矩阵，需克服相关的带宽问题，需要的内存较大的带宽问题，需要的内存较大.5.对大位移、大应变问题需进行对大位移、大应变问题需进行大量的计算大量的计算.GeoHohaiNew Features in FLAC Version 6.0使用使用Intel Fortran compiler拥有更快的计算速度拥有更快的计算速度自动网格重画功能，解决自动网格重画功能，解决 bad-geometry 问题问题.新的模拟颗粒土材料的硬化模型新的模拟颗粒土材料的硬化模型更新的通用网格生成工具更新的通用网格生成工具22GeoHohaiNew Features in FLAC3D Version 3.11.多处理器的并行计算功能多处理器的并行计算功能 2.新新结构单元类型结构单元类型“Embedded Liner”提供两个方提供两个方向的接触作用，可以很好地模拟挡土墙向的接触作用，可以很好地模拟挡土墙3.对对四面体单元采用新的混合离散方法四面体单元采用新的混合离散方法“Nodal Mixed Discretization”提供塑性问题更精确的解提供塑性问题更精确的解答答4.64位程序位程序5.包含包含命令手册、命令手册、FISH手册和应用实例的帮助手册和应用实例的帮助23GeoHohaiNew Features in FLAC3D Version 4.0模拟颗粒状材料的硬化模型模拟颗粒状材料的硬化模型自动网格重画功能，解决自动网格重画功能，解决 bad-geometry 问题问题.改进的改进的interface更快的渗流计算更快的渗流计算更新的动力计算功能更新的动力计算功能24GeoHohaiLagrangian格式动量平衡方程25F(t)m牛顿运动定律对于连续体在静力平衡条件下，加速度项为0，方程变为平衡方程GeoHohai自由落体的模拟26G=mgS=1/2gt2=20m命令流：命令流：config dyngen zon bri size 1 1 1ini x mul 0.1 y m 0.1 z m 0.1model elasprop bulk 3e8 shear 1e8ini dens 1000set grav 0 0-10solve age 2GeoHohai自由落体的模拟(movie)27GeoHohaiFLAC3D中模型术语中模型术语28节点gridpoint：节点zone：单元boundary：边界GeoHohaiFLAC3D的求解过程29平衡方程(动量方程)应力应变关系(本构模型)Gauss定律单元积分应变率速度节点力新的应力对所有的网格节点对所有单元GeoHohaiFLAC3D中的本构模型开挖模型开挖模型null 3个弹性模型个弹性模型各向同性弹性各向同性弹性横观各向同性弹性横观各向同性弹性正交各向同性弹性正交各向同性弹性 8个弹塑性模型个弹塑性模型Drucker-Prager模型、模型、Morh-Coulomb模型、应变硬化模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化软化模型、遍布节理模型、双线性应变硬化/软化遍布软化遍布节理模型、修正剑桥模型和胡克布朗模型节理模型、修正剑桥模型和胡克布朗模型30GeoHohaiFLAC3D中的本构模型31GeoHohai一个最简单的例子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 grav 0 0-10solve ;求解求解app nstr-10e4 ran z 3 x 1 2 y 1 2solveplo con zd ;后处理后处理切片功能切片功能32RUN FLAC3DGeoHohai分析问题的过程33建立网格建立网格初始条件初始条件边界条件边界条件初始应力平衡初始应力平衡外荷载外荷载求解求解前处理前处理后处理后处理GeoHohaiFLAC3D的文件格式的文件格式n保存文件保存文件(*.sav)含有所有状态变量和用户定义条件的二进制文件含有所有状态变量和用户定义条件的二进制文件 n数据文件数据文件(*.dat)数据文件由用户创建的一种数据文件由用户创建的一种ASC格式的文件，格式的文件，它包括一系列的用于描述所分析问题的它包括一系列的用于描述所分析问题的FLAC3D命令命令 nFISH文件文件(*.fis)FISH程序文件程序文件nFLAC3D文件文件(*.flac3d)FLAC3D的网格信息文件的网格信息文件n历史记录文件历史记录文件(*.his)记录输入输出历史值的文件记录输入输出历史值的文件 n图形文件图形文件 图形文件图形文件(各种标准格式各种标准格式)n电影文件电影文件(*.dcx)AVI或或PCX图像文件，这些图像文件可以当作电图像文件，这些图像文件可以当作电影放映影放映34GeoHohai初始应力的生成为什么要单独列出？为什么要单独列出？分析过程中出现的很多问题都与初始应力是否合理有分析过程中出现的很多问题都与初始应力是否合理有关关手册中的例子五花八门手册中的例子五花八门是所有后续分析的基础是所有后续分析的基础！生成方法生成方法弹性求解弹性求解更改强度参数的弹塑性求解更改强度参数的弹塑性求解设置初始应力的弹塑性求解设置初始应力的弹塑性求解存在水压力的初始应力生成存在水压力的初始应力生成水下建筑的初始应力生成水下建筑的初始应力生成35GeoHohai弹性求解模型尺寸单元数量密度KGu112(m3)112200030MPa10MPa0.3536gen zon bri size 1 1 2m elasprop 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-10solveStep=162s sz=-40e3s sx=-21.54e3GeoHohai更改强度参数的弹塑性求解模型尺寸单元数量密度KGcfu112(m3)112200030MPa10MPa10kPa150.3537gen 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 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0-10solveprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0solveStep=163s sz=-40e3s sx=-21.54e3Or:solve elasticGeoHohai设置初始应力的弹塑性求解模型尺寸单元数量密度KGcfu112(m3)112200030MPa10MPa10kPa150.3538gen 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 0fix 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-10solveStep=0s sz=-40e3s sx=-20e3GeoHohai存在水压力的初始应力生成(1)模型尺寸单元数量饱和密度KGcfu水位线孔隙率112(m3)112200030MPa10MPa10kPa150.351m0.539genzonebricksize112modelmohrpropbulk3e7shear1e7coh10e3fri15ten0fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000ranz01inidens1500ranz12iniszz-35e3grad0020e3ranz01inisyy-22.5e3grad0015e3ranz01inisxx-22.5e3grad0015e3ranz01iniszz-30e3grad0015e3ranz12inisyy-15e3grad007.5e3ranz12inisxx-15e3grad007.5e3ranz12inipp10e3grad00-10e3ranz01setgrav00-10solveStep=0s sz=-35e3s sx=-22.5e3r rd=r rs nsr rfGeoHohai存在水压力的初始应力生成(2)模型尺寸单元数量饱和密度KGcfu水位线孔隙率112(m3)112200030MPa10MPa10kPa150.351m0.540config fluidgen zon bri size 1 1 2model elasprop bu 3e7 sh 1e7ini dens 1500model fl_isoini fdens=1000 fmod 0prop por 0.5set grav 0 0-10water dens 1500water table face 0 0 1,0 1 1,1 1 1,1 0 1 ini pp 10e3 grad 0 0-10e3 ran z 0 1ini szz-30e3 grad 0 0 15e3 ran z 1 2.set fluid offsolveStep=142s sz=-40e3s sx=-24e3单元数较少产生的误差单元数较少产生的误差GeoHohai水下建筑的初始应力生成模型尺寸单元数量饱和密度KGcfu水位线112(m3)112200030MPa10MPa10kPa150.353m41gen 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 2000 ran z 0 2ini szz-50e3 grad 0 0 20e3 ran z 0 2ini syy-40e3 grad 0 0 15e3 ran z 0 2ini sxx-40e3 grad 0 0 15e3 ran z 0 2ini pp 30e3 grad 0 0-10e3 ran z 0 2app nstress-10e3 ran z 2set grav 0 0-10solveStep=0s sz=-50e3s sx=-40e3GeoHohai前后处理42GeoHohai前后处理基本前后处理基本前后处理命令操作命令操作菜单操作菜单操作dd&dipattach&merge外界模型的导入外界模型的导入复杂模型的网格检查复杂模型的网格检查43GeoHohaiFLAC3D的前后处理命令驱动命令驱动(推荐推荐)程序控制程序控制图形界面接口图形界面接口计算模型输出计算模型输出指定本构模型及参数指定本构模型及参数指定初始条件及边界条件，指定结构单元指定初始条件及边界条件，指定结构单元指定接触面指定接触面指定自定义变量及函数指定自定义变量及函数(FISH)求解过程的变量跟踪求解过程的变量跟踪进行求解进行求解模型输出模型输出44GeoHohai菜单驱动(计算模式)45命令栏GeoHohai菜单驱动(Plot)46GeoHohaiFLAC3D的前处理的前处理47FLAC3D网格生成的关键特征：网格生成的关键特征：FLAC3D 是命令驱动.使用FLAC3D内置基元进行形状组合可形成复杂网格.用户自定义FISH 函数可以用来修改基元网格以创建更为复杂的网格.第三方软件导入.GeoHohaiFLAC3D网格基元48块体块体退化块体退化块体楔体楔体金字塔金字塔四面体四面体柱体柱体径向块体径向块体径向隧道径向隧道径向柱体径向柱体柱状壳体柱状壳体柱状交叉柱状交叉隧道交叉隧道交叉GeoHohai49GeoHohai50块体块体径向柱体径向柱体柱状交叉柱状交叉GeoHohaiTecplot的后处理51GeoHohai52GeoHohaiTecplot后处理53GeoHohai54GeoHohai55GeoHohai56GeoHohaiTecplot的云图和等值线57TecplotTecplot中云图的效果与中云图的效果与flac3dflac3d的的shade onshade on的效果差的效果差不多不多,但是其出三维等值线的功能是卓越的但是其出三维等值线的功能是卓越的.GeoHohai58GeoHohaiTecplot的切片功能59TecplotTecplot的的sliceslice切片切片功能功能.和和CADCAD一样一样,可以任意可以任意切剖面出切剖面出图图.最大优最大优点是可以点是可以几个剖面几个剖面同时出图同时出图.GeoHohai60整体和截面的整体和截面的数据以单元形数据以单元形式存在列表中式存在列表中,可自由选择和可自由选择和组合多个单元组合多个单元出图出图.GeoHohaidd 和 dipSimwe.Com上有近百的讨论贴上有近百的讨论贴地质上的概念，倾向和倾角地质上的概念，倾向和倾角建议用建议用ori和和norm代替代替后处理的切片功能后处理的切片功能plo set plane ori(*,*,*)norm (*,*,*)plo con zdis plane概念清晰，使用快捷概念清晰，使用快捷其他应用其他应用reflect网格网格建立水面建立水面61xyzorinormGeoHohaiattach 和 mergeattach用于连接次节点用于连接次节点 (sub-grid)允许网格存在一定随意性允许网格存在一定随意性要求成比例要求成比例(1:2,1:3等等)可用于模型的检查可用于模型的检查慎用慎用merge用于节点的合并用于节点的合并外来模型导入的精度差异外来模型导入的精度差异tolerance的设置的设置62sub-grid 4:25:2tol121mergGeoHohai外界模型的导入复杂网格的生成难度大复杂网格的生成难度大接口编写接口编写不同软件之间的精度差异不同软件之间的精度差异Group的定义的定义采用采用.flac3d的文件格式的文件格式.flac3d文件的格式文件的格式G 1 1.0e+00 1.0e+00 1.0e+00Z B8 1 2 3 4 5 6 7 8ZGROUP Soil1 2 3impgrid&expgrid仅限于仅限于网格网格633w zones 3secGeoHohai复杂网格的检查FLAC3D生成的复杂网格生成的复杂网格attach face无接触面时检查整体模型是否存在无接触面时检查整体模型是否存在sub-grid有接触面时给定范围进行检查有接触面时给定范围进行检查其他软件导入的复杂模型其他软件导入的复杂模型网格划分的检查网格划分的检查gen merge弹性模型弹性模型model elastic求解求解“独立独立”节点节点“畸形畸形”单元单元64GeoHohai模型的检查FLAC本身的本身的Check功能十分有限功能十分有限错误提示很少错误提示很少十分开放的工作平台十分开放的工作平台 (自由落体自由落体)检查的基本步骤检查的基本步骤网格检查网格检查(如前所述如前所述)边界条件检查边界条件检查速度约束条件速度约束条件 plo gpfix red sk荷载条件荷载条件 plo fap red sk模型检查模型检查模型赋值模型赋值plo block model参数赋值参数赋值plo block prop*初始应力检查初始应力检查(如后所述如后所述)65GeoHohai初始应力检查办法初始应力的计算时间不会初始应力的计算时间不会“非常长非常长”经常检查模型的响应经常检查模型的响应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)接触面单元接触面单元66GeoHohai第二讲FLAC3D接触面、接触面、FISH语言、流固语言、流固耦合分析耦合分析GeoHohai接触面（interface）GeoHohai接触面单元应用范围应用范围原理原理建模方法建模方法推荐方法推荐方法复杂内部接触面的设置方法复杂内部接触面的设置方法参数选择参数选择单桩承载力分析单桩承载力分析挡土墙的接触面设置挡土墙的接触面设置思路：未知问题的分析方法思路：未知问题的分析方法GeoHohai接触面单元的用途网格不连续网格不连续岩体介质中的解理、断层、岩层面岩体介质中的解理、断层、岩层面地基与土体的接触地基与土体的接触箱、槽及其内充填物的接触箱、槽及其内充填物的接触空间中无变形的固定空间中无变形的固定“障碍障碍”GeoHohai接触面的原理三角形单元三角形单元(无厚度无厚度!)参数较多参数较多三种工作模式三种工作模式粘结界面粘结界面粘接滑移粘接滑移库伦滑动库伦滑动GeoHohai接触单元模型的建立(1)关键关键要形成同一位置的两个节点要形成同一位置的两个节点(面面)“移来移去移来移去”(推荐推荐)建两个分开的模型建两个分开的模型建立接触单元建立接触单元通过通过INI*add使模型接触使模型接触注意注意dist的含义的含义NO merge,NO attach!接触面dist1234GeoHohai接触单元模型的建立(2)“导来导去导来导去”利用利用expgrid,impgrid命令进行网格导出与导入命令进行网格导出与导入配合配合DELETE命令命令适于适于内部接触面内部接触面的建立，或的建立，或其他前处理其他前处理工具建立的网格工具建立的网格GeoHohai接触单元模型的建立(3)GEN separate INTERFACE wrap指定正确的指定正确的groupGeoHohai接触面参数的确定虚构的为了合并节点而设置的接触面虚构的为了合并节点而设置的接触面Kn=ks=10*max(K+4/3G)/D Dzmin真实的刚性接触面真实的刚性接触面如料仓下料如料仓下料c,D,Tension重要，重要，kn,ks不重要不重要真实的柔性接触面真实的柔性接触面断层；断层；水力劈裂水力劈裂材料材料试验得到参数试验得到参数对于对于kn,ks：岩石断层：岩石断层10100MPa/m(粘土粘土);100GPa(岩石岩石)反分析方法反分析方法：通过断层中岩石的变形与原岩的变形：通过断层中岩石的变形与原岩的变形GeoHohai单桩承载力分析软土地基软土地基bulk 1.6878E6 shear 3.6167E5 coh 15E3 fric 12dens 1.73E3桩体桩体bulk 5e9shear 3.75e9dens 2.5e30.5m 8m 10m 20m GeoHohai计算过程施加桩顶荷载施加桩顶荷载计算结果计算结果GeoHohai影响因素7.50E+03fricgk3coh/0.70.7fric10g10k23.00E+0420100e100e1cohfrickskn水平因素GeoHohai计算工况设计1233(9)3123(8)2313(7)2132(6)1322(5)3212(4)3331(3)2221(2)1111(1)cohfrickskn水平因素工况GeoHohai计算结果Ks取取1 Ks取取2 Ks取取3 96.9%1233(9)22.4%3123(8)54.1%2313(7)96.9%2132(6)7.1%1322(5)37.8%3212(4)96.9%3331(3)42.9%2221(2)37.8%1111(1)Ra差异差异度度cohfrickskn水平因素水平因素GeoHohai最优方案2222111122101122最优方最优方案案17.3%6.8%72.8%11.9%极差极差52.4%52.7%96.9%57.8%k364.6%59.2%24.1%47.3%k247.3%52.4%43.2%59.2%k1cohfrickskn水平因水平因素素GeoHohai合理步骤单桩分析单桩分析简单网格简单网格接触面参数接触面参数多次试算多次试算理论、实测理论、实测加密网格加密网格接触参数接触参数理想结果理想结果群桩分析群桩分析Pile结构单元结构单元单元参数单元参数理想结果理想结果Pile结构单元结构单元单元参数单元参数理想结果理想结果GeoHohai挡土墙的接触面设置对于未知问题的分析思路对于未知问题的分析思路3个独立的接触面个独立的接触面3个同个同ID的接触面的接触面2个独立的接触面个独立的接触面2个独立的接触面并进行底部个独立的接触面并进行底部merge挡墙挡墙 土体土体 123WallSoilGeoHohai挡土墙的接触面设置(2)Z=5.16cmZ=1.14cmZ=不能初始平衡不能初始平衡 Z=0.35cm？GeoHohai挡土墙的接触面设置(3)X=1.95cmX=1.91cmX=1.47cm？GeoHohai挡土墙的接触面设置(3)前两种方法的差别实前两种方法的差别实质质3个个ID的独立接触面在的独立接触面在相同位置产生互不影相同位置产生互不影响的两个节点响的两个节点共同共同ID的接触面在相的接触面在相同位置自动设置为同位置自动设置为1个个节点节点最终的结论需要您自最终的结论需要您自己去判断！己去判断！3 interfaces,2 IDs 1 interface,1 ID GeoHohaiFISH语言GeoHohaiFISH语言简介软件自带的编程语言软件自带的编程语言是否一定要学？是否一定要学？视情况而定，需要时查询视情况而定，需要时查询FISH变量即可变量即可语法简单语法简单xxxend_xxx注意事项注意事项与与FLAC本身的关键字冲突本身的关键字冲突保留字不可缩写保留字不可缩写变量可不定义，因此注意检查程序变量可不定义，因此注意检查程序print fishtable,extra等命令使用等命令使用GeoHohai一个最简单的FISH程序def abc abc=1+2*3 abcd=1.0/2.0endabcprint fish数据格式GeoHohai函数与变量1.都可以在都可以在FISH函数中进行赋值，赋值操作与常规的编程函数中进行赋值，赋值操作与常规的编程语言类似，按照运算符的优先级先后顺序来执行。语言类似，按照运算符的优先级先后顺序来执行。2.函数和变量的赋值遵守数据类型的规则，即整型的计算结函数和变量的赋值遵守数据类型的规则，即整型的计算结果为整型，浮点型的计算结果为浮点型，因此读者在进行果为整型，浮点型的计算结果为浮点型，因此读者在进行除法运算、开方运算时都需要将数据类型设置为浮点型，除法运算、开方运算时都需要将数据类型设置为浮点型，数字尽量使用小数点以保证运算正确。数字尽量使用小数点以保证运算正确。3.变量和函数名的命名规则是不能以数字开头，不能含有中变量和函数名的命名规则是不能以数字开头，不能含有中文，并且不能包含如下的字符。文，并且不能包含如下的字符。.,*/+-=#();4.变量和函数名不能与变量和函数名不能与FLAC3D、FISH的保留字相冲突，不的保留字相冲突，不要采用过于简单的单词，比如要采用过于简单的单词，比如a，hist等，这些都与保留字等，这些都与保留字相冲突。相冲突。5.即使程序中存在与保留字相冲突的变量，即使程序中存在与保留字相冲突的变量，FLAC3D也不会也不会提供任何提示，所以提醒读者在编制提供任何提示，所以提醒读者在编制FISH程序时程序时尽量使尽量使用较长的、复杂的变量和函数名用较长的、复杂的变量和函数名。GeoHohai函数与变量5.对变量进行赋值时，不能使用当前函数的函数名放在对变量进行赋值时，不能使用当前函数的函数名放在“=”的右边，比如采用下面的定义的右边，比如采用下面的定义abcd=abc+1.0在在FISH程序执行时会提出错误，因为这样会形成递归调程序执行时会提出错误，因为这样会形成递归调用，这种调用方式在用，这种调用方式在FISH程序中是不允许的。程序中是不允许的。6.变量和函数的作用是变量和函数的作用是全局全局的，在命令中的任何地方修改的，在命令中的任何地方修改变量的值都会立即生效，因此在实际应用中尽量避免不变量的值都会立即生效，因此在实际应用中尽量避免不同的函数中含有相同的变量，因为这样可能会造成赋值同的函数中含有相同的变量，因为这样可能会造成赋值错误，并难以检查。错误，并难以检查。7.在在FLAC3D中可以用如下的命令来引用中可以用如下的命令来引用FISH函数和变量函数和变量PRINT用于查看函数和变量的数值；用于查看函数和变量的数值；HISTORY命令可以对函数和变量的数值进行记录；命令可以对函数和变量的数值进行记录；SET命令用于变量的赋值。命令用于变量的赋值。GeoHohai主要语句选择语句选择语句CASEOF 表达式表达式默认语句默认语句CASE n1表达式的值为表达式的值为n1时的语句时的语句CASE n2表达式的值为表达式的值为n2时的语句时的语句ENDCASEGeoHohai主要语句条件语句条件语句IF 条件表达式条件表达式 THENELSEENDIFFISH中条件运算符没有“并”、“或”、“否”这样的符号ifaa1.0ifaa2.0执行语句endifendif表达“1aamaxdisp_valuemaxdisp_value=disp_gpmaxdisp_gpid=gp_id(p_gp)endifp_gp=gp_next(p_gp)endloopendfind_max_dispprintmaxdisp_valuemaxdisp_gpidGeoHohaiFISH的编写习惯第一步第二步def abcendabcdef abc p_gp=gp_head loop while p_gp#null p_gp=gp_next(p_gp)endloopendabc第三步第四步def abc p_gp=gp_head loop while p_gp#null command endcommand p_gp=gp_next(p_gp)endloopendabcdef abc p_gp=gp_head loop while p_gp#null command app nstress endcommand p_gp=gp_next(p_gp)endloopendabcGeoHohaiFISH与建模两个圆形隧道的连接部分两个圆形隧道的连接部分 变直径的隧道部分变直径的隧道部分 GeoHohaiFISH检查主要采用主要采用PRINT