FLAC3D建模分析教学教材.ppt

FLAC3D建模分析

复杂网格的生成难度大接口编写不同软件之间的精度差异Group的定义采用.flac3d的文件格式.flac3d文件的格式G11.0e+001.0e+001.0e+00ZB812345678ZGROUPSoil123impgridexpgrid仅限于网格外部模型导入规则模型检查

数据格式*FLAC3DgridproducedbyFLAC3D*GRIDPOINTSG10.000000000e+0000.000000000e+0000.000000000e+000G21.000000000e+0000.000000000e+0000.000000000e+000G30.000000000e+0001.000000000e+0000.000000000e+000....................................*ZONESZB8112345678ZB8229571081112ZB83913101114121516....................................ZGROUPtop501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574....................................

通过内置语言FISH生成网格

通过第三方软件生成3DSHOP生成网格

通过第三方软件生成Ansys生成网格

通过第三方软件生成Surfer生成网格

后处理

通过第三方软件生成Gocad生成网格

应力初始化应力初始化的条件：具有充足的水文地质条件应力初始化的优点：由于边界及非线性因素，能够正确合理反映地质水文条件；对于非线性问题，能够明显提高计算速度。

ModelelasticModelm;inistressModelm

侧压力系数的确定

侧压力系数的确定常见土侧压力系数经验值

侧压力系数的确定

构造应力初始化(深埋情况)

接触单元的用途岩体介质中的节理、断层、岩层面基础与土体的接触箱、槽及其内充填物的接触抵触体之间的界面空间中无变形的固定“障碍”

接触单元原理一.接触单元几何几物理性质节理单元无厚度FLAC:2节点(双面)；FLAC3D:3节点(单面)“等代面积”的概念物理力学参数较多三种工作模式1.粘结界面2.粘接滑移3.库伦滑动

几何模型建立二.FLAC3D程序“移去移来”法“导去导来”法FLAC3D.V3.0及以上版本新功能

接触单元参数的确定一.刚度的选取经验估计方法：从岩体性质来估计节理的刚度(无厚度)。假定岩体的变形是由无扰动岩体的变形和节理变形两部分组成，可以写成：再假定岩体内部只有一组节理，加载方向和节理面垂直，可以用下式估计节理刚度：

接触单元参数的确定从填充材料来估计节理的刚度(有厚度且假定填充材料为弹性)：二.屈服参数的选取

接触单元参数的确定

FISH语言简介软件自带的编程语言是否一定要学？视情况而定，需要时查询FISH变量即可语法简单xxx……end_xxx注意事项与FLAC本身的关键字冲突保留字不可缩写变量可不定义，因此注意检查程序printfishtable,extra等命令使用

FISH应用介绍改变水头边界

FISH应用介绍隧洞开挖应力释放

流固耦合分析无渗流模式不设置CONFIGFluid孔压不改变INITIALppWATERtableWATERdensitySETgravityWATERtableface手动设置干湿密度设置CONFIGfluid瞬态渗流分析有效应力计算不排水计算设置土体干密度渗流模型MODELfl_isotropicMODELfl_anisotropicMODELfl_null渗流模式

安全系数计算

安全系数计算FOS=1.61FOS=2.71

耦合分析20m10mdeframpramp=min(1.0,float(st