FLAC3D5.0详细命令
[实用参考]Flac3d-5.0常用命令集锦.doc
![[实用参考]Flac3d-5.0常用命令集锦.doc](https://img.taocdn.com/s3/m/dee971461711cc7930b7163b.png)
建模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlyt hefilenameneedbeenteredwiththe CALL command.②不在同一个文件夹，全路径：Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfol der\file.dat).Undo;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negative z-direction.（垂直向下！常用的）下下面面这这代代码码，，是是沿沿着着--y y方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续p110,-10,-20...p2-10,10,-20...p3-10,-10,0plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体p110,-10,-20p2-10,10,-20...p3-10,-10,0p410,10,-20...p5-10,10,10p610,-10,0...p710,10,10plotcurrentplotPlot01plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01y01z01;显示指定范围内各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：eachhistoryisnumberedsequentiallyfrom1asitisenteredviathe HISTORY co mmand.②查找显示所有的his命令：ReturntotheFlac3D>promptandtypelisthist foralistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。
FLAC3D5.0模型及输入参数说明
1.1模型参数代码可参考manual中各个章节的command命令及说明，注意单位。
用prop 赋值。
1.1.1各向同性弹性模型1.1.2横向同行弹性模型1.1.3正交各向异性弹性模型1.1.4德鲁克-普拉格模型1.1.5摩尔-库伦模型1.1.6多节理模型1.1.7应变硬化/软化模型1.1.8双线性应变硬化/软化多节理模型1.1.9D-Y模型1.1.10修正剑桥模型1.1.11纯动力学模型1.1.12经典粘弹性模型1.1.13粘弹性模型v1.0 可编辑可修改4mkshear Maxwell切边模量，G M5mviscosity Maxwell动力粘度，ηM1.1.14二分幂律模型二分幂律模型的材料参数（Power Law –MODEL mechanical power）1a_1常数，A12a_2常数，A23bulk弹性体积模量，K4n_1指数，n15n_2指数，n26rs_1参考应力，σ1ref7rs_2参考应力，σ2ref8shear弹性剪切模量，G1.1.15蠕变模型蠕变模型材料参数（WIPP Model –MODEL mechanical wipp）1act_energy活化能，Q2a_wipp常数，A3b_wipp常数，B4bulk弹性体积模量，K5d_wipp常数，D6e_dot_star临界稳定状态蠕变率，7gas_c气体常数，R8n_wipp指数，n9shear弹性剪切模量，G10temp温度，T下列参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变1.1.16Burger、蠕变组合材料模型1.1.17幂律模型8n_2指数，n29rs_1参考应力，σ1ref10rs_2参考应力，σ2ref11shear弹性剪切模量，G12tension抗拉强度，σt1.1.18粘塑形模型粘塑形模型的材料参数（WIPP-Creep Viscoplastic Model –MODEL mechanical pwipp）1act_energy活化能，Q2a_wipp常数，A3b_wipp常数，B4bulk弹性体积模量，K5d_wipp常数，D6e_dot_star临界稳定状态蠕变率，7gas_c气体常数，R8kshear材料参数，KΦ9n_wipp指数，n10kdil材料参数，q k11kvol材料参数，qΦ12shear弹性切变模量，G13temp温度，T14tension抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率3es_plastic累积塑性切应变4et_plastic累积塑性拉应变1.1.19碎盐变形模型碎盐变形模型的材料参数（Crushed-Salt Model –MODEL mechanical cwipp）1act_energy活化能，Q2a_wipp常数，A3b_f最终体积模量，K f4b_wipp常数，B5b0蠕变压实系数，B06b1蠕变压实系数，B17b2蠕变压实系数，B28bulk弹性体积模量，K9d_f最终密度，ρf10d_wipp常数，D11e_dot_star临界稳定状态蠕变率，12gas_c气体常数，R13n_wipp指数，n14rho密度，ρ15s_f最终切变模量，G f16shear弹性切变模量，G17temp温度，T以下计算参数可以显示、绘图和通过fish访问1frac_d当前碎片密度，ρd2s_g1蠕变压实参数，G3s_k1蠕变压实参数，K1.1.20均质流体模型均质流体模型的材料参数1permeability等方向渗透性，k1.1.21各向异性流体模型1.1.22均质热导模型1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。
（完整版）FLAC3D命令流（整理版）
（完整版）FLAC3D命令流（整理版）1、怎样查看模型？答：plot grid 可以查看⽹格，plot grid num 可以查看节点号。
2、请问在圆柱体四周如何施加约束条件？答：可以⽤fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把⼀个PLOT的图像数据导出来以便⽤其他软件绘图？答：⽤set log on 命令，把数据导出来，转到excel⾥处理⼀下，然后⽤surfer或者什么作图软件绘制就⾏了。
4、⽤命令建⽴模型后，如何显⽰点的坐标？答：使⽤plo blo gro gpnum on命令5、关于gauss_dev对性质进⾏⾼斯正态分布的问题？答：根据⼿册上的说明：下⾯的命令设定⼀个平均摩擦⾓为40度，标准⽅差是±5％。
则命令如下：prop friction 40 gauss_dev 2问题：请问gauss_dev 2中的2是如何计算的？如果把±5％改为±10％，则命令应如何写？40×5％＝26、reflect问题问：gen zone radbrick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fillplot surfgen zone reflect dip 0 dd 90 （对xy⾯做镜像）gen zone reflect dip 90 dd 90 （对yz⾯做镜像）（1）dd表⽰y轴正向顺时针到那条射线的夹⾓，dip表⽰对称参照⾯与xy平⾯的夹⾓，对称参照⾯与xy平⾯的夹⾓在xy平⾯的投影是⼀条射线。
[实用参考]Flac3d-5.0常用命令集锦
![[实用参考]Flac3d-5.0常用命令集锦](https://img.taocdn.com/s3/m/0be3be44f46527d3240ce044.png)
建模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlPthefilenamen eedbeenteredwiththe CALL command. ②不在同一个文件夹，全路径：Otherwise,thefilemaPbecalledbPspecifPingits completepath(e.g.,c:\mPfolder\file.dat ). Undo ;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr ;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negative z -direction.（垂直向下！常用的）下下面面这这代代码码，，是是沿沿着着--P P 方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续 p110,-10,-20... p2-10,10,-20... p3-10,-10,0 plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体 p110,-10,-20p2-10,10,-20... p3-10,-10,0p410,10,-20... p5-10,10,10p610,-10,0... p710,10,10plotcurrentplotPlot01 plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01P01z01;显示指定范围内各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：eachhistorPisnumberedsequentiallP from1asitisenteredviathe HISTORP command.②查找显示所有的his命令：ReturntotheFlac3D>promptandtPpelisthistforalistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。
(新)Flac3D5.0笔记(注释)
FLAC3D5.0导入dat、txt文件的方法
FLAC3D5.0导入.sav文件的方法
输入命令流：Restore 文件名.sav
FLAC3D5.0导入.f3prj文件的方法
视图窗口关闭后，如何再次呈现？输入命令plot zone
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名
做切片
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面，双击I nterface显示接触面
接触面属性项被激活后，可以更改属性隐藏网格显示
输出网格
如果点击了透明，将不显示Information
只要不点击透明，则显示Information
等同于
点击显示节点坐标
fish简介
之前要有（空格）之前加是Flac3D3.0和5.0区别而在之间没有（空格）
按下列步骤取消fish语言前要加
点击fish
软件的刷新如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果
Origin原点normal法向量（该点与原点的坐标差）
对于结构单元如桩，若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。
flac3d5.0软件隧道支护与开挖命令流
flac3d5.0软件隧道支护与开挖命令流;-----------------------------------------------------; ---- Excavation and Support for a Shallow Tunnel ---;-----------------------------------------------------new ;新建项目set fish autocreate offtitle 'Excavation and Support for a Shallow Tunnel' ;定义题目; generate primitive components of grid; concrete liner - upper tunnelgen zon cshell p0 0 0 0 p1 7 0 0 p2 0 51 0 p3 0 0 5.5 &dim 5 5 5 5 size 2 51 10group zone 'concrete liner';; upper tunnelgen zon cylinder p0 0 0 0 p1 5 0 0 p2 0 51 0 p3 0 0 5 &size 5 51 10group zone tunnel range group 'concrete liner' not;; lower tunnel & linergen zone brick p0 0 0 -4.5 p1 add 7 0 0 p2 add 0 51 0 p3 add 0 0 4.5 &size 7 51 3;; surrounding rock (8 primitives)gen zon radcyl p0 0 0 0 p1 27 0 0 p2 0 51 0 p3 0 0 25 &dim 7 5.5 7 5.5 size 5 51 10 8 rat 1 1 1 1.3;gen zone brick p0 7 0 -4.5 p1 27 0 -15 p2 add 0 51 0 p3 7 0 0 &p4 27 51 -15 p5 7 51 0 p6 27 0 0 p7 27 51 0 &size 8 51 3 ratio 1.3 1 1;gen zone brick p0 0 0 -15 p1 add 27 0 0 p2 add 0 51 0 p3 0 0 -4.5 &p4 27 51 -15 p5 0 51 -4.5 p6 7 0 -4.5 p7 7 51 -4.5 &size 7 51 8 rat 1 1 0.7692307692307692;gen zon brick p0 0 0 25 p1 add 27 0 0 p2 add 0 51 0 p3 add 0 0 10 &size 5 51 2;gen zon bric p0 27 0 25 p1 add 17 0 0 p2 add 0 51 0 p3 add 0 0 10 &size 2 51 2 rat 2 1 1;gen zon bric p0 27 0 -15 p1 add 17 0 0 p2 add 0 51 0 p3 add 0 0 40 &size 2 51 8 rat 2 1 1;gen zon bric p0 27 0 -40 p1 add 17 0 0 p2 add 0 51 0 p3 add 0 0 25 &size 2 51 2 rat 2 1 0.5;gen zon bric p0 0 0 -40 p1 add 27 0 0 p2 add 0 51 0 p3 add 0 0 25 &size 7 51 2 rat 1 1 0.5;; assign names to groups of zonesgroup zone rock range group 'concrete liner' not group tunnel not;; assign Mohr-Coulomb material modelmodel mech mohrpro bulk 50e6 she 18e6 fric 20 coh 25e3 ten 0 dil 0 range z 25 35pro bulk 4e8 she 1.5e8 fric 20 coh 50e3 ten 5e3 dil 3 range z -50 25; assign boundary conditions ;施加边界条件，后面可以直接修改为具体的数字fix x range x -.1 .1fix x range x 43.9 44.1fix z range z -40.1 -39.9fix y range y -.1 .1fix y range y 50.9 51.1; assign initial stress state ;初始应力状态set grav 0 0 -10ini density 2200ini szz -770e3 grad 0 0 22000ini sxx -770e3 grad 0 0 22000ini syy -385e3 grad 0 0 11000 ;施加初始应力; monitor variables in model ;模型中变量的监控hist add unbal ;监控不平衡力hist add gp zdisp 0 0 5.5hist add gp xdisp 7 0 0hist add gp zdisp 0 0 0hist add gp zdisp 0 0 35hist add gp zdisp 0 30 5.5hist add gp xdisp 7 30 0hist add gp zdisp 0 30 0hist add gp zdisp 0 30 35hist add gp zdisp 0 12 35hist add gp zdisp 0 18 35hist add gp zdisp 0 24 35hist add gp zdisp 0 36 35hist add gp zdisp 5 30 35hist add gp zdisp 10 30 35 ;监控以上这些点的x及z方向位移变化;sav geom1;def conc_parm ;定义支护参数，parm即parameter，参数的意思global bmc = 20.7e9 ;定义体积模量为全局变量，b代表bulk，m代表modulus，c代表concreteglobal smc = 12.6e9 ;定义剪切模量为全局变量，s代表shear，m代表modulus，c代表concreteend@conc_parm;; define the locations of cable patterns 1, 2 and 3;def cab_parm ;定义锚杆参数global x_b = get_array(4,3) ;定义数组（4,3），即锚杆的位置global z_b = get_array(4,3)global y0 = -3 ;将锚杆的位置以数组的形式表示出来（x,y,z）x_b(1,1) = 0.8x_b(2,1) = 2.1x_b(3,1) = 3.5x_b(4,1) = 5.5z_b(1,1) = 5.5z_b(2,1) = 2.4z_b(3,1) = 4.7z_b(4,1) = 1.5x_b(1,2) = 0.8x_b(2,2) = 0.8x_b(3,2) = 3.5x_b(4,2) = 5.5z_b(1,2) = 0.6z_b(2,2) = 4.0z_b(3,2) = 2.4z_b(4,2) = 0.6x_b(1,3) = 0.8x_b(2,3) = 2.6x_b(3,3) = 5.0x_b(4,3) = 3.5z_b(1,3) = 2.4z_b(2,3) = 4.0z_b(3,3) = 3.0z_b(4,3) = 0.6 ;对每根锚杆的位置进行赋值enddef inip(iidx) ;定义初始锚杆位置global x1 = x_b(1,iidx)global x2 = x_b(2,iidx)global x3 = x_b(3,iidx)global x4 = x_b(4,iidx)global z1 = z_b(1,iidx)global z2 = z_b(2,iidx)global z3 = z_b(3,iidx)global z4 = z_b(4,iidx)end@cab_parm;; install initial cables ;安装初始锚杆;def ins_cab ;定义要初始安装的锚杆，其中ins_cab表示install initial cablesglobal iidx ;定义全局变量iidxglobal cab_seg ;定义锚杆划分单元数global cab_seg_m ;锚杆长度（有待进一步确认）loop iidx (1,3)inip(iidx)cab_seg = cab_seg_m-3*(3-iidx)global y1 = 0.global y2 = float(cab_seg)commandsel cable id @iidx begin @x1 @y1 @z1 end @x1 @y2 @z1nseg @cab_segsel cable id @iidx begin @x2 @y1 @z2 end @x2 @y2 @z2 nseg @cab_segsel cable id @iidx begin @x3 @y1 @z3 end @x3 @y2 @z3 nseg @cab_segsel cable id @iidx begin @x4 @y1 @z4 end @x4 @y2 @z4 nseg @cab_segsel cable pro emod 45e9 xcarea 1.57e-3 gr_per 1.0 &yten 25e4 gr_k 17.5e6 gr_c 20e4 range id @iidx ;施作初始锚杆end_commandend_loopendset @cab_seg_m 15@ins_cab; install pre-support concrete ;预支护;sel shell id 10 group rock range cyl end1 0 0 -1.5 end2 0 1 -1.5 rad 7.4 &cyl end1 0 0 -1.5 end2 0 1 -1.5 rad 6.7 not &z -0.1 6sel shell prop isotropic 10.5e9,0.25 thickness 0.3 density 2500 def monitglobal ipt_surf = gp_near(0,30,35) ;地表global ipt_crown = gp_near(0,30,5.5) ;拱顶global ipt_spring = gp_near(7,30,0)end@monitsave m_initable 1 name 'ground surface at tunnel center line'table 2 name 'tunnel crown' ;隧道拱顶table 3 name 'tunnel sidewall' ;定义表格的名字，隧道边墙;; FISH function to control excavation and support sequence def excavy0 = y0+3local cut_i = y0/3+1global cutloop cut (cut_i,16)local cut_cur = cutlocal ii = out(' EXCAVATION STEP ' + string(cut))y0 = 3*(cut-1)y1 = y0+3global yp0 = y0+1global yp1 = y1+1global ys0 = yp0-3global ys1 = yp1-3global yc0 = y0-3global yc1 = y1-3global id_ = 10; id_ = 10*(cut+1) ; use if shells unconnectedcommand; install pre support concretesel shell id @id_ group rock &**************************@yp1-1.5rad7.4&*********************@yp1-1.5rad6.7not&z -0.1 6sel shell prop isotropic 10.5e9,0.25 thickness 0.3 density 2500 &ran y @yp0 @yp1; excavate next cutmodel mech null range group tunnel y @y0 @y1model mech null range group 'concrete liner' y @y0 @y1; delete-cables in the excavated areasel delete cable range id 1 y @y0 @y1sel delete cable range id 2 y @y0 @y1sel delete cable range id 3 y @y0 @y1end_commandlocal cut_1 = cut-1iidx=int(cut_1-3*(cut_1/3))+1y2=min(y1+15,51)inip(iidx)ii = out(' CABLE BOLT PATTERN '+string(iidx))commandsel delete cable range id @iidx; install new cablessel cable id @iidx begin @x1 @y1 @z1 end @x1 @y2 @z1 nseg @cab_seg_msel cable id @iidx begin @x2 @y1 @z2 end @x2 @y2 @z2 nseg @cab_seg_msel cable id @iidx begin @x3 @y1 @z3 end @x3 @y2 @z3 nseg @cab_seg_msel cable id @iidx begin @x4 @y1 @z4 end @x4 @y2 @z4 nseg @cab_seg_msel cable pro emod 45e9 xcarea 1.57e-3 gr_per 1.0 &yten 25e4 gr_k 17.5e6 gr_c 20e4 ran id @iidx; shotcretesel shell prop isotropic 10.5e9,0.25 thickness 0.5 density 2500 &ran y @ys0 @ys1end_commandif cut > 1 thencommand; concrete linermodel mech el range group 'concrete liner' y @yc0 @yc1prop bulk @bmc sh @smc range group 'concrete liner' y @yc0 @yc1end_commandend_ifcommandstep 3000end_command; store displacements in tables ;将位移储存在表格中xtable(1,cut) = 3.0 * cutytable(1,cut) = gp_zdisp(ipt_surf)xtable(2,cut) = 3.0 * cutytable(2,cut) = gp_zdisp(ipt_crown)xtable(3,cut) = 3.0 * cutytable(3,cut) = gp_zdisp(ipt_spring)commandsave m1end_commandif cut=5 thencommandsave m1_15end_command end_ifif cut=9 then command save m1_27 end_command end_ifif cut=10 then command save m1_30 end_command end_ifend_loopend@excav return。
FLAC3D V5.0技术参数
Itasca(武汉)咨询有限公司高级、三维连续介质力学方法新版主要功能概述计算效率获取更多信息，请访问官方网站： (+86) 027 –8786 2349itasca@建模技术引擎随机裂隙网络(Discrete Fracture Network ，DFN)岩土工程分析、支护结构设计自V3.1版本（2008）至今，FLAC 3D 并行计算技术已历经近5年的发展并趋于成熟，同步于计算机硬件技术更新的日新月异，V5.0改进、优化多核并行计算内核，使得程序分析效率得到较大幅度提升，与以往版本不同地，新版并行分析支持流体和结构单元计算，进一步满足大规模工程分析对计算效率的苛刻要求。
与客户诉求一致地，新版建模技术引擎倾注了FLAC 3D 技术团队更多的辛劳，并预示着Itasca 将进一步加大对Itasca 全线产品前处理技术的开发力度。
新版FLAC 3D 建模引擎植入两项技术改进，分别为CAD 辅助建模技术及其模型接口平台KUBRIX ，前者指嵌入CAD 绘图环境并与DXF 、STL 文件有效兼容以辅助建模，后者指专业前处理工具KUBRIX ，包括KUBRIX GEO 和KUBRIX HEX 。
Itasca 三维程序产品全线嵌入随机裂隙网络技术以强化地质结构面模拟能力，即FLAC 3D 、3DEC 和PFC 系列。
作为采用连续介质力学方法等效处理非连续问题的最新手段，该项技术将进一步拓展FLAC 3D 程序的应用广度和深度，使其成为岩土工程领域某些非常规、高难度问题的研究平台，如压裂设计、岩体力学特性的不均匀性、各向异性及其尺寸效应等。
后处理其他多种后处理方式改进、扩展功能：(1)强化鼠标控制操作；(2)新增PostSript 、VRML 、SVG 、Excel 格式输出接口，且图片输出精度得到提升；(3)多项张量后处理，包括第二不变量、Von Mises 应力等；(4)等值面创建和模型结果与CAD 对象间数据插值、喷涂渲染功能，特别适用于工程方案评估、比较分析。
FLAC3D5.0模型及输入参数说明
1.1模型参数代码可参考manual中各个章节的command命令及说明，注意单位。
用prop 赋值。
1.1.1各向同性弹性模型16 ttable 塑性拉应变-抗拉强度的表号下列参数可以显示、绘图和通过fish访问1 es_plastic 塑性切应变2 et_plastic 塑性拉应变3 ff_count 检测切应变反向的数4 ff_cvd 体应变，εvd经典粘弹性模型的材料参数（Classical Viscoelastic （Maxwell Substance）–MODEL mechanical viscous）1 bulk 弹性体积模量，K2 shear 弹性剪切模量，G3 viscosity 动力粘度，η粘弹性模型的材料参数（Burgers Model –MODEL mechanical burgers）1 bulk 弹性体积模量，K2 kshear Kelvin弹性剪切模量，G K3 kviscosity Kelvin动力粘度，ηK4 mkshear Maxwell切边模量，G M5 mviscosity Maxwell动力粘度，ηM二分幂律模型的材料参数（Power Law –MODEL mechanical power）1 a_1 常数，A12 a_2 常数，A23 bulk 弹性体积模量，K4 n_1 指数，n15 n_2 指数，n26 rs_1 参考应力，σ1ref7 rs_2 参考应力，σ2ref8 shear 弹性剪切模量，G蠕变模型材料参数（WIPP Model –MODEL mechanical wipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star临界稳定状态蠕变率，7 gas_c 气体常数，R8 n_wipp 指数，n9 shear 弹性剪切模量，G10 temp 温度，T下列参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率Burger、蠕变组合材料模型的材料参数（Burgers-Creep Viscoplastic Model –MODEL mechanical cvisc）1 bulk 弹性体积模量，K2 cohesion 聚力，c3 density 密度，ρ4 dilation 剪胀角，Ψ5 friction 摩擦角，Φ6 kshear Kelvin弹性剪切模量，G K7 kviscosity Kelvin粘度，ηK8 shear 弹性剪切模量，G9 tension 抗拉强度，σt10 mviscosity Maxwell动力粘度，ηM下列计算参数可以显示、绘图和通过fish访问1 es_plastic 累积塑性切应变2 et_plastic 累积塑性拉应变幂律模型的材料参数（Power-Law Viscoplastic Model –MODEL mechanical cpower）1 a_1 常数，A12 a_2 常数，A23 bulk 弹性体积模量，K4 cohesion 聚力，c5 dilation 剪胀角，Ψ6 friction 摩擦角，Φ7 n_1 指数，n18 n_2 指数，n29 rs_1 参考应力，σ1ref10 rs_2 参考应力，σ2ref11 shear 弹性剪切模量，G12 tension 抗拉强度，σt粘塑形模型的材料参数（WIPP-Creep Viscoplastic Model –MODEL mechanical pwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star临界稳定状态蠕变率，7 gas_c 气体常数，R8 kshear 材料参数，KΦ9 n_wipp 指数，n10 kdil 材料参数，q k11 kvol 材料参数，qΦ12 shear 弹性切变模量，G13 temp 温度，T14 tension 抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率3 es_plastic 累积塑性切应变4 et_plastic 累积塑性拉应变碎盐变形模型的材料参数（Crushed-Salt Model –MODEL mechanical cwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_f 最终体积模量，K f4 b_wipp 常数，B5 b0 蠕变压实系数，B06 b1 蠕变压实系数，B17 b2 蠕变压实系数，B28 bulk 弹性体积模量，K9 d_f 最终密度，ρf10 d_wipp 常数，D11 e_dot_star临界稳定状态蠕变率，12 gas_c 气体常数，R13 n_wipp 指数，n14 rho 密度，ρ15 s_f 最终切变模量，G f16 shear 弹性切变模量，G17 temp 温度，T以下计算参数可以显示、绘图和通过fish访问1 frac_d 当前碎片密度，ρd2 s_g1 蠕变压实参数，G3 s_k1 蠕变压实参数，K均质流体模型的材料参数1 permeability 等方向渗透性，k2 porosity 孔隙率，n（默认时，n=0.5）各向异性流体模型的材料参数1 fdd k1-k2的平面倾向2 fdip k1-k2的平面倾角1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。
FLAC3D5.0模型及输入参数说明
1.1模型参数代码可参考manual中各个章节的command命令及说明，注意单位。
用prop 赋值。
15 tension 抗拉强度，σt16 ttable 塑性拉应变-抗拉强度的表号下列参数可以显示、绘图和通过fish访问1 es_plastic 塑性切应变2 et_plastic 塑性拉应变3 ff_count 检测切应变反向的数4 ff_cvd 体应变，εvd1.1.12经典粘弹性模型经典粘弹性模型的材料参数（Classical Viscoelastic （Maxwell Substance）–MODEL mechanical viscous）1 bulk 弹性体积模量，K2 shear 弹性剪切模量，G3 viscosity 动力粘度，η1.1.13粘弹性模型粘弹性模型的材料参数（Burgers Model –MODEL mechanical burgers）1 bulk 弹性体积模量，K2 kshear Kelvin弹性剪切模量，G K3 kviscosity Kelvin动力粘度，ηK4 mkshear Maxwell切边模量，G M5 mviscosity Maxwell动力粘度，ηM1.1.14二分幂律模型二分幂律模型的材料参数（Power Law –MODEL mechanical power）1 a_1 常数，A12 a_2 常数，A23 bulk 弹性体积模量，K4 n_1 指数，n15 n_2 指数，n26 rs_1 参考应力，σ1ref7 rs_2 参考应力，σ2ref8 shear 弹性剪切模量，G1.1.15蠕变模型蠕变模型材料参数（WIPP Model –MODEL mechanical wipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star临界稳定状态蠕变率，7 gas_c 气体常数，R8 n_wipp 指数，n9 shear 弹性剪切模量，G10 temp 温度，T6 e_dot_star临界稳定状态蠕变率，7 gas_c 气体常数，R8 kshear 材料参数，KΦ9 n_wipp 指数，n10 kdil 材料参数，q k11 kvol 材料参数，qΦ12 shear 弹性切变模量，G13 temp 温度，T14 tension 抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率3 es_plastic 累积塑性切应变4 et_plastic 累积塑性拉应变1.1.19碎盐变形模型碎盐变形模型的材料参数（Crushed-Salt Model –MODEL mechanical cwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_f 最终体积模量，K f4 b_wipp 常数，B5 b0 蠕变压实系数，B06 b1 蠕变压实系数，B17 b2 蠕变压实系数，B28 bulk 弹性体积模量，K9 d_f 最终密度，ρf10 d_wipp 常数，D11 e_dot_star临界稳定状态蠕变率，12 gas_c 气体常数，R13 n_wipp 指数，n14 rho 密度，ρ15 s_f 最终切变模量，G f16 shear 弹性切变模量，G17 temp 温度，T以下计算参数可以显示、绘图和通过fish访问1 frac_d 当前碎片密度，ρd2 s_g1 蠕变压实参数，G3 s_k1 蠕变压实参数，K均质流体模型的材料参数1 permeability 等方向渗透性，k2 porosity 孔隙率，n（默认时，n=0.5）1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。
FLAC3D5.0详细命令
FLAC3D5.0详细命令1、plot zones and gridpoint velocities directly, simply type plot zone velocity2、plot linecol or colorspecifies the color of the line.on/off turns line drawing for this option on or off.stylespecifies the linestyle used.widthispecifies the width of the line. This value must be between 0 and 10. A value of 0 indicates the smallest possible line (1 pixel).3、plot textNot all of the following keywords are available in all contexts. col orcolorspecifies which color the text will be drawn in.fa mily stringspecifies the name of the font family (for instance, “Ar- ial”or “Times New Roman”). (Note that underscores must be used in place of spaces in the name.)on/off turns text drawing on or off.4、plot create plot Myplotplot text on textplot text "this is my plot";（怎样显示文字？不会）5、保存图片plot create plot GravV;plot cut add plane dip=90 dd=0 origin=3,4,0plot cut add plane origin=119,90,110 nor 0 1 0plot set orient 120, 220, 0plot set center 45, 45, 0plot add zone plane behindplaneonplane offplot add bcontszzplot add axesplot post view GravV filename "103F1068"plot post filename "721F01F3-fft-input.ps" manualplot post view GravV filename "GravV.ps" manualplot post plot Plot01 file 127F1017.ps manual上面保存的是.ps的图片，不容易打开，我知道用photoshop可以打开，麻烦。
flac3d--5.0常用命令集锦
建模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：If the called file is located in the same folder as the FLAC3D project file, then only the file name needbe entered with the CALL command.②不在同一个文件夹，全路径：Otherwise, the file may be called by specifying its complete path (e.g., c:\myfolder\file.dat).Undo ;撤销上一条命令2、创建旋转缩放视图3、建模命令model mech mohr ;莫尔库伦模型model mech elastic ;弹性模型set grav 0, 0, -9.81 ; 重力加速度 negative z-direction. （垂直向下！常用的）gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续p1 10, -10, -20 ...p2 -10, 10, -20 ...p3 -10, -10, 0plot zonegen zone brick size 6,8,8 p0 -10, -10, -20 ... ;不规则六面体p1 10, -10, -20 p2 -10, 10, -20 ...p3 -10, -10, 0 p4 10, 10, -20 ...p5 -10, 10, 10 p6 10, -10, 0 ...p7 10, 10, 10plot current plot Plot01plot clearplot zoneUndo ;撤销命令set logfile 127X1001.txtset log on truncateset log offlist zone prin range x 0 1 y 0 1 z 0 1 ;显示指定范围内各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：each history is numbered sequentially from 1 as it is entered via the HISTORY command.②查找显示所有的his命令：Return to the Flac3D> prompt and typelist histfor a listing of the histories and their corresponding numbers.hist nstep 5 ;每5步记录1次。
FLAC3D5.0模型及输入全参数说明书
1・1模型参数代码可参考manual 中各个章节的comma nd 命令及说明，注意单位。
用 prop 赋值。
1.1.1各向同性弹性模型1.1・4德鲁克•普拉格模型1.1.5摩尔•库伦模型1.1.6多节理模型1.1.7应变硬化/软化模型1.1.18粘塑形模型1.2模型适用说明遍布节理模型适用于Mohr・Coulomb材料来明确显示力在各个方向上的差异性。
双线性软化应变遍布节理模型综合了软化应变Mohr・Coulomb模型和遍布节理模型，这种模型包含面向矩阵和遍布节理的一个双线性断裂点集。
改进的Cam-clay模型反映了形变度和抗破坏能力对体积变化的影响。
Mohr-Coulomb模型最适用于一般工程研究，同时，Mohr-Coulomb的聚力和摩擦角参数相对于地质工程材料的其它属性，更容易获得。
软化应变和遍布节理塑性模型实际上是Mohr・Coulomb 模型的变形，这些模型如果在附加材料参数的值较高时将得出与Mohr-Coulomb模型同样的结果。
Druck-Prager模型是一个相对于Mohr・Coulomb模型的破坏标准的简化体，但是它一般不适于用来描述地质工程材料的破坏情况。
它主要是用来把FLAC3D与其它一些有Druck-Prager模型但却没有Mohr-Coulomb模型的数学软件作比较。
在摩擦力为零的时候请注意，此时Mohr・Coulomb模型退化为Tresca模型，而Druck-Prager模型退化为Von Mises模型。
Druck-Prager模型和Mohr・Coulomb模型是计算起来效率最高的塑性模型，而其它的塑性模型在计算时却需要更多的存和额外的时间。
例如，塑性应变不能在Mohr・Coulomb模型中直接计算出来（参见附录G ）。
如果需要计算塑性应变，则必需要用应变软化模型。
这种模型主要是用于破坏后的情况对工程影响重大的工程活动中，如弯曲柱、开采塌落或回填研究。
Druck-Prager模型、Mohr・Coulomb模型、遍布节理模型、应变软化Mohr-Coulomb模型和双线性应变软化遍布节理模型中的拉伸破坏标准是相同的。
FLAC3D5.0模型及输入参数说明
材料模型及关键字说明（model ）MODEL 关键字说明力学模型anisotropic 横向同性弹性模型cam-clay 修正剑桥模型doubleyield D-Y模型drucker 德鲁克普拉格模型elactic 各向同性弹性模型mohr 摩尔-库伦模型null 空模型orthotropic 正交各向异性弹性模型ss 应变硬化/软化模型subliquitous 双线性应变硬化/软化多节理模型ubiquitous 多节理模型finn 纯动力学模型蠕变模型burger Burger材料粘弹性模型cpower 幂律材料模型cvisc Burger、蠕变组合材料模型cwipp 岩盐变形模型power 二分幂律模型pwipp 粘塑性模型viscous 经典粘弹性模型wipp 蠕变模型流体模型fl_anistropic 各向异性流体模型fl_isotropic 各向同性流体模型fl_null 空流体模型热力模型th_anistropic 均质热导模型th_null 空热导模型Example:model mohr range x=1,10 y=1,10 z=1,10 （group name）prop bulk 25.0e9 shear 5.2e9 coh 0.23e6 dil 20.0 frict 38.0 tension 0.0e61.1模型参数代码可参考manual中各个章节的command命令及说明，注意单位。
用prop 赋值。
各向同性弹性模型各向同性弹性材料参数（Isotropic Elastic –MODEL mechanical elastic）1 bulk 体积模量，或松散系数K2 shear 切边模量，G横向同行弹性模型横向同行弹性模型材料参数（Transversely Isotropic Elastic –MODEL mechanical anisotropic）1 dd 同性平面的倾向2 dip 同性平面的倾角3 e1 同性平面的弹性模量4 e3 垂直同性平面的弹性模量5 g12 切变模量6 nu12 同性平面内施力时的泊松比7 nu13 垂直同性平面内施力时的泊松比正交各向异性弹性模型正交各向异性弹性模型材料参数（Orthotropic Elastic –MODEL mechanical orthotropic）1 dd 轴1-‘2’所定义平面的倾向2 dip 轴1-‘2’所定义平面的倾角3 e1 1’方向弹性模量4 e2 2’方向弹性模量5 e3 3’方向弹性模量6 g12 平行于轴1-‘2’平面的切变模量7 g13 平行于轴1-‘3’平面的切变模量8 g23 平行于轴2-‘3’平面的切变模量9 nu12 沿2’方向施力，1’方向的泊松比10 nu13 沿3’方向施力，1’方向的泊松比11 nu23 沿3’方向施力，2’方向的泊松比12 nx 轴1-‘2’所定义平面单位法线x分量13 ny 轴1-‘2’所定义平面单位法线y分量14 nz 轴1-‘2’所定义平面单位法线z分量15 rot 旋转角德鲁克-普拉格模型德鲁克-普拉格模型的材料参数（Drucker-Prager –MODEL mechanical drucker）1 bulk 弹性体积模量，K2 ksnear 材料参数，kφ3 qdil 材料参数，qφ4 qvol 材料参数，qΨ5 shear 弹性切变模量，K6 tension 抗拉强度，σt摩尔-库伦模型摩尔-库伦模型的材料参数（Mohr-Coulomb –MODEL mechanical mohr）1 bulk 弹性体积模量，K2 cohesion 内聚力，c3 dilation 剪胀角，Ψ4 friction 内摩擦角，Φ5 shear 弹性切边模量，G6 tension 抗拉强度，σt多节理模型多节理模型的材料参数（Ubiquitous-Joint –MODEL mechanical ubiquitous）1 bulk 弹性体积模量，K2 cohesion 内聚力，c3 dilation 剪胀角，Ψ4 friction 内摩擦角，Φ5 jcohesion 节理内聚力，c j6 jddirection 弱面dip方向（倾向）7 jdilation 节理剪胀角，Ψj8 jdip 弱面dip角度（倾角）9 jfriction 节理摩擦角，Φj10 jnx 弱面单位法线x分量11 jny 弱面单位法线y分量12 jnz 弱面单位法线z分量13 jtension 抗拉强度，σsj14 shear 弹性切边模量，G15 tension 抗拉强度，σt应变硬化/软化模型应变硬化/软化模型的材料参数（Strain-Hardening/Softening –MODEL mechanical ssoftening）1 bulk 弹性体积模量，K2 cohesion 内聚力，c3 ctable 塑性剪切应变-内聚力的表号4 diation 剪胀角，Ψ5 dtable 塑性剪切应变-剪胀角的表号6 friction 内摩擦角，Φ7 ftable 塑性剪切应变-摩擦角的表号8 shear 弹性切边模量，G9 tension 抗拉强度，σs10 ttable 塑性剪切应变-抗拉强度的表号双线性应变硬化/软化多节理模型双线性应变硬化/软化多节理模型的材料参数（Bilinear Strain-Hardening/Softening Ubiquitous-Joint –MODEL mechanical subiquitous）1 bijoint =0，为线性节理，默认=1，为双线性节理2 bimatrix =0，线性矩阵=1，双线性矩阵3 bulk 弹性体积模量，K4 c2table 塑性剪切应变-内聚力c2的表号5 cjtable 节理塑性剪切应变-节理内聚力c j1的表号6 cj2table 节理塑性剪切应变-节理内聚力c j2的表号7 cohesion 内聚力，c18 co2 内聚力，c29 ctable 塑性剪切应变-内聚力c1的表号10 d2table 塑性剪切应变-剪胀角Ψ2的表号11 di2 剪胀角，Ψ212 dilation 剪胀角，Ψ113 djtable 节理塑性剪切应变-节理剪胀角Ψj1的表号14 dj2table 节理塑性剪切应变-节理剪胀角Ψj2的表号15 dtable 塑性剪切应变-剪胀角Ψ1的表号16 f2table 塑性剪切应变-摩擦角Φ2的表号17 fjtable 节理塑性剪切应变-节理摩擦角Φj1的表号18 fj2table 节理塑性剪切应变-节理摩擦角Φj2的表号19 fr2 摩擦角Φ220 friction 摩擦角Φ121 ftable 塑性剪切应变-摩擦角Φ1表号22 jc2 节理内聚力c j123 jcohesion 节理内聚力c j224 jddirection 弱面dip方向（倾向）25 jdilation 节理剪胀角Ψj126 jdip 弱面dip角度（倾角）27 jd2 节理剪胀角，Ψj128 jfriction 节理摩擦角，Φj129 jf2 节理摩擦角，Φj230 jnx 弱面单位法线x分量31 jny 弱面单位法线y分量32 jnz 弱面单位法线z分量33 jtension 节理抗拉强度，σsj34 shear 弹性切边模量，G35 tension 抗拉强度，σs36 tjtable 节理塑性剪切应变-节理抗拉强度σtj的表号37 ttable 节理塑性剪切应变-节理抗拉强度σt的表号下列参数可以显示、绘图或者fish访问1 es_plastic 塑性切应变2 et_plastic 塑性拉应变3 etj_plastic 节理塑性拉应变4 esj_plastic 节理塑性切应变D-Y模型D-Y模型的材料参数（Double-Yield –MODEL mechanical doubleyield）1 bulk 弹性体积模量，K2 cap_pressure cap压力，p c3 cohesion 内聚力，c4 cptable 塑性体应变-cap压力的表号5 ctable 塑性切应变-内聚力的表号6 dilation 剪胀角，Ψ7 dtable 塑性切应变-剪胀角的表号8 ev_plastic 塑性体应变总量9 friction 内摩擦角，Φ10 ftable 塑性切应变-摩擦角的表号11 multiplier 当前塑性-cap模量与弹性体积和切变模量的倍数，R12 shear 最大弹性切变模量，G13 tension 抗拉强度，σs14 ttable 塑性拉应变-抗拉强度的表号下列计算参数可以显示、绘图和通过fish访问1 es_plastic 累积塑性切应变2 et_plastic 累积塑性拉应变3 ev_plastic 累积塑性体应变修正剑桥模型修正剑桥模型的材料参数（Modified Cam-Clay –MODEL mechanical cam-clay）1 bulk_bcund 最大的弹性体积模量，K max2 cv 初始容积，v c3 kappa 弹性膨胀线斜率，k4 lamda 常态固结线斜率，λ5 mm 摩擦常数，M6 mpc 预固结压力，p c07 mp1 预固结压力，p18 mv_1 指定在参考压力常态固结线的容积v A9 poisson 泊松比，ν10 shear 弹性剪切模量，G下列参数可以显示、绘图以及fish访问1 bulk 体积模量，K2 cam_cp 当前平均有效应力3 cam_ev 累积总容积应变4 camev_cp 累积塑性容积应变5 cq 当前平均差分应力纯动力学模型纯动力学模型的材料参数（config dynamic, model mech finn）1 bulk 弹性体积模量，K2 cohesion 内聚力，c3 ctable 弹性切应变-内聚力的表号4 dilation 剪胀角，Ψ5 dtable 塑性切应变-剪胀角Ψ的表号6 ff_c1 常量，c17 ff_c2 常量，c28 ff_c3 常量，c39 ff_c4 常量，c410 ff_latensy 反向之间的最小时间步数11 ff_switch =0：Martin（1995）公式=1：Byme（1991）公式12 friction 内摩擦角，Φ13 ftable 塑性切应变-摩擦角的表号14 shear 最大弹性切变模量，G15 tension 抗拉强度，σt16 ttable 塑性拉应变-抗拉强度的表号下列参数可以显示、绘图和通过fish访问1 es_plastic 塑性切应变2 et_plastic 塑性拉应变3 ff_count 检测切应变反向的数4 ff_cvd 体应变，εvd经典粘弹性模型经典粘弹性模型的材料参数（Classical Viscoelastic （Maxwell Substance）–MODEL mechanical viscous）1 bulk 弹性体积模量，K2 shear 弹性剪切模量，G3 viscosity 动力粘度，η粘弹性模型粘弹性模型的材料参数（Burgers Model –MODEL mechanical burgers）1 bulk 弹性体积模量，K2 kshear Kelvin弹性剪切模量，G K3 kviscosity Kelvin动力粘度，ηK4 mkshear Maxwell切边模量，G M5 mviscosity Maxwell动力粘度，ηM二分幂律模型二分幂律模型的材料参数（Power Law –MODEL mechanical power）1 a_1 常数，A12 a_2 常数，A23 bulk 弹性体积模量，K4 n_1 指数，n15 n_2 指数，n26 rs_1 参考应力，σ1ref7 rs_2 参考应力，σ2ref8 shear 弹性剪切模量，G蠕变模型蠕变模型材料参数（WIPP Model –MODEL mechanical wipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star 临界稳定状态蠕变率，ϵ̇ss∗7 gas_c 气体常数，R8 n_wipp 指数，n9 shear 弹性剪切模量，G10 temp 温度，T下列参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率Burger、蠕变组合材料模型Burger、蠕变组合材料模型的材料参数（Burgers-Creep Viscoplastic Model –MODEL mechanical cvisc）1 bulk 弹性体积模量，K2 cohesion 内聚力，c3 density 密度，ρ4 dilation 剪胀角，Ψ5 friction 内摩擦角，Φ6 kshear Kelvin弹性剪切模量，G K7 kviscosity Kelvin粘度，ηK8 shear 弹性剪切模量，G9 tension 抗拉强度，σt10 mviscosity Maxwell动力粘度，ηM下列计算参数可以显示、绘图和通过fish访问1 es_plastic 累积塑性切应变2 et_plastic 累积塑性拉应变幂律模型幂律模型的材料参数（Power-Law Viscoplastic Model –MODEL mechanical cpower）1 a_1 常数，A12 a_2 常数，A23 bulk 弹性体积模量，K4 cohesion 内聚力，c5 dilation 剪胀角，Ψ6 friction 内摩擦角，Φ7 n_1 指数，n18 n_2 指数，n29 rs_1 参考应力，σ1ref10 rs_2 参考应力，σ2ref11 shear 弹性剪切模量，G12 tension 抗拉强度，σt粘塑形模型粘塑形模型的材料参数（WIPP-Creep Viscoplastic Model –MODEL mechanical pwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_wipp 常数，B4 bulk 弹性体积模量，K5 d_wipp 常数，D6 e_dot_star 临界稳定状态蠕变率，ϵ̇ss∗7 gas_c 气体常数，R8 kshear 材料参数，KΦ9 n_wipp 指数，n10 kdil 材料参数，q k11 kvol 材料参数，qΦ12 shear 弹性切变模量，G13 temp 温度，T14 tension 抗拉强度，σt以下计算参数可以显示、绘图和通过fish访问1 e_prime 累积主蠕变应变2 e_rate 累积主蠕变应变率3 es_plastic 累积塑性切应变4 et_plastic 累积塑性拉应变碎盐变形模型碎盐变形模型的材料参数（Crushed-Salt Model –MODEL mechanical cwipp）1 act_energy 活化能，Q2 a_wipp 常数，A3 b_f 最终体积模量，K f4 b_wipp 常数，B5 b0 蠕变压实系数，B06 b1 蠕变压实系数，B17 b2 蠕变压实系数，B28 bulk 弹性体积模量，K9 d_f 最终密度，ρf10 d_wipp 常数，D11 e_dot_star 临界稳定状态蠕变率，ϵ̇ss∗12 gas_c 气体常数，R13 n_wipp 指数，n14 rho 密度，ρ15 s_f 最终切变模量，G f16 shear 弹性切变模量，G17 temp 温度，T以下计算参数可以显示、绘图和通过fish访问1 frac_d 当前碎片密度，ρd2 s_g1 蠕变压实参数，G3 s_k1 蠕变压实参数，K均质流体模型均质流体模型的材料参数1 permeability 等方向渗透性，k2 porosity 孔隙率，n（默认时，n=0.5）各向异性流体模型各向异性流体模型的材料参数1 fdd k1-k2的平面倾向2 fdip k1-k2的平面倾角3 frot k1轴和倾角矢量的转角4 h1 k1方向的渗透性5 h2 k2方向的渗透性6 h3 k3方向的渗透性均质热导模型均质热导模型的材料参数1 conductivity 等方向传热系数，K2 spec_heat 比热容，C v1.2模型适用说明遍布节理模型适用于Mohr-Coulomb材料来明确显示力在各个方向上的差异性。
flac3d5.0常用命令集锦
建 模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：If the called file is located in the same folder as the FLAC3D project file, then only the file name need be entered with the CALL command.②不在同一个文件夹，全路径：Otherwise, the file may be called by specifying its complete path (e.g., c:\myfolder\file.dat ).Undo ;撤销上一条命令2、创建 旋转 缩放视图3、建模命令model mech mohr ;莫尔库伦模型model mech elastic ;弹性模型set grav 0, 0, -9.81 ; 重力加速度 negative z -direction. （垂直向下！常用的）下下面面这这代代码码，，是是沿沿着着--y y 方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续p1 10, -10, -20 ...p2 -10, 10, -20 ...p3 -10, -10, 0plot zonegen zone brick size 6,8,8 p0 -10, -10, -20 ... ;不规则六面体p1 10, -10, -20 p2 -10, 10, -20 ...p3 -10, -10, 0 p4 10, 10, -20 ...p5 -10, 10, 10 p6 10, -10, 0 ...p7 10, 10, 10plot current plot Plot01plot clearplot zoneUndo ;撤销命令set logfile 127X1001.txtset log on truncateset log offlist zone prin range x 0 1 y 0 1 z 0 1 ;显示指定围各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：each history is numbered sequentially from 1 as it is entered via the HISTORY command.②查找显示所有的his命令：Return to the Flac3D> prompt and typelist histfor a listing of the histories and their corresponding numbers.hist nstep 5 ;每5步记录1次。
Flac3D5.0操作方法笔记
FLAC3D5.0导入dat、txt文件的方法
FLAC3D5.0导入.sav文件的方法
输入命令流：Restore 文件名.sav
FLAC3D5.0导入.f3prj文件的方法
视图窗口关闭后，如何再次呈现？输入命令plot zone
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名
做切片
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面，双击I nterface显示接触面
接触面属性项被激活后，可以更改属性
隐藏网格显示
输出网格
如果点击了透明，将不显示Information
只要不点击透明，则显示Information
等同于
点击显示节点坐标
fish简介
@之前要有（空格）之前加@是Flac3D3.0和5.0区别
而在之间没有（空格）
按下列步骤取消fish语言前要加@
点击fish
软件的刷新
如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果
Origin原点normal法向量（该点与原点的坐标差）
对于结构单元如桩，若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。
02-0_FLAC3D5.0基本操作方法vs应用流程
1.1、建模—Generate
网格形状
1.1、建模—Generate
brick
gen zone brick size 6 8 8 plot zone
1.1、建模—Generate
dbrick
gen zone dbrick size 6 6 6 p0 0 0 0 p1 6 0 0 p2 0 6 0 … p3 0 0 6 p4 6 6 0 p5 0 6 6 p6 6 0 6 plot zone
1.1、建模—Extruder
新建图形名称
导入DXF文件 在绘图区域描画出平面模型 生成平面网格并调整 材料分组
切换侧面并划分
生成三维模型
1.1、建模—Extruder
删除图形
选择 节点 中间点 网格生成
网格调整
查看分组 正面切换 侧面切换 模型生成
1.1、建模—Extruder
elastic, isotropic model 各向同性线弹性模型
Mohr-Coulomb model 摩尔库伦模型 Hoek-Brown model 霍克布朗模型 Druck-Prager model 德鲁克普拉格模型 null model常被用于开挖过程的模拟，elastic model赋予材料的线弹 性力学行为，mohr model赋予材料塑性力学行为
38
3.3、根据求解结果进行后处理
• • • • 云图（or内力图） 变形图 监测的关键点数据 等等。
4、实例操作 ----地表基坑的开挖
中间点 线型： 1.line 2.curve 3.arc ① ② 网格 ③ 节点 ②
1.1、建模—Extruder
平面图和纵向 图上的每个对 象（点、线、 方块）都有自 己的属性可供 调整。！！
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1、plot zones and gridpoint velocities directly, simply typeplot zone velocity2、plot linecol or colorspecifies the color of the line.on/off turns line drawing for this option on or off.stylespecifies the linestyle used.width ispecifies the width of the line. This value must bebetween 0 and 10. A value of 0 indicates the smallestpossible line (1 pixel).3、plot textNot all of the following keywords are available in all contexts.col or colorspecifies which color the text will be drawn in.fa mily stringspecifies the name of the font family (for instance, “Ar-ial”or “Times New Roman”). (Note that underscoresmust be used in place of spaces in the name.)on/off turns text drawing on or off.4、plot create plot Myplotplot text on textplot text"this is my plot";（怎样显示文字？不会）5、保存图片plot create plot GravV;plot cut add plane dip=90 dd=0 origin=3,4,0plot cut add plane origin=119,90,110nor 0 1 0plot set orient 120, 220, 0plot set center 45, 45, 0plot add zone plane behindplane onplane offplot add bcont szzplot add axesplot post view GravV filename "103F1068"plot post filename "721F01F3-fft-input.ps" manualplot post view GravV filename "GravV.ps" manualplot post plot Plot01 file 127F1017.ps manual上面保存的是.ps的图片，不容易打开，我知道用photoshop可以打开，麻烦。
下面保存为bitmap图片：plot create plot GravVplot add zoneplot add geom location addlabel "Polygon" lightcyanplot set movieextension bmp;(bmp or jpeg/jpg/png/ppm/xbm/xpm);设置输出图片类型实际用时我用：plot set movieextension bmp或plot set movieextension pngplot bitmap filename "GravV.png" size 1024 768;输出图片，名为GravV.png，并控制大小6、显示cad保存的.dxf多面图元在flac3d500中并在左面bar上加Polygon字样plot add geom location addlabel "Polygon" lightcyan7、切片plot cut add plane name string origin x y z nor x y zplot set movieex tension bmp(or jpeg/jpg/png/ppm/xbm/xpm)plot cut add plane origin x y z nor x y z ;plot cut add plane origin 5100 -9460 -40 nor 1 0 08、显示不同破坏准则的模型单元体plot bcontour property propname young9、显示单元plot add sel geom10、按组显示单元体plot add zone colorby group addlabel "Default" rgb (255,255,255) addlabel "Count 1" rgb (160,160,160)lighting off range z 0 3.9或plot add zgroup11、显示cad图元：plot add geometry contour12、显示坐标轴plot add axes13、显示模型表面szz应力云图,力都用下面的命令plot add isozone szz---显示内部竖直压应力：pl zcon szz14、显示模型表面的z向位移，只要是位移或速度或加速度都用下面的plot add isosurface zdis15、显示块体的表面szzplot bcontour szz16、显示锚杆力plot add sel cabblock force line width 4plot cont disp ramp grey reverse on trans 50 backface onplot add sel cabblock force line width 5 ramp grey reverse on17、清除视图plot clear18、改变hist线属性plot add hist 1 linestyle color black19、用刚才视图窗口plot current plot Plot01plot clearplot zone colorby model20、range 的用法; --- range definitions ---range name drift x -39.63 0 y -19.51 -15.24 z -1.57 2.09range name room x -14.63 -9.14 y -15.24 15.24 z -1.57 3.92 range name bhole cylinder end1 0 0 0 end2 -9.14 0 0 rad 0.46 range name empty union nrange drift nrange room nrange bhole model mech null range nrange emptysolvesave wipe1_a21、切片plot create plot DispCont1plot set job onplot set center (0,0,0) eye (-2,0,0) orientation (90,90,0)plot cut add plane origin (-0.1,0,0) normal (1,0,0) caption on plot add cont disp ramp grey reverse on wireframe off fill &plane onplane on &maximum 0.1 minimum 0 interval 0.01plot create plot Model1plot set job onplot set center (-18,-6,-23) eye (60,-72,62) orientation (130,-50,0) plot cut add plane origin (0,0,0) normal (0,0,1) caption onplot add zone colorby uniform addlabel "Zone" white &plane onplane on behindplane on &backface offplot create plot Model2plot set job onplot set center (-30,1.5,-7) eye (94,-102,86) orientation (130,-50,0) plot add zone colorby uniform addlabel "Zone" whiteplot create plot Model3plot set job onplot set center (-94,53,-56) eye (6,-4,10) orientation (120,-60,0) plot cut add plane origin (0,0,0) normal (0,0,1) caption onplot add zone colorby uniform addlabel "Zone" white &plane onplane on behindplane on &backface off22、清零initial xdisp 0 ydisp 0 zdisp 0initial xvel 0 yvel 0 zvel 023、历史history add id 1 gp pp 2 18 1 ;i=3, j=19history add id 2 gp pp 10 8 1 ;i=12, j=9history add id 3 gp pp 15 21 1 ;i=17, j=22history add id 4 gp pp 15 39 1 ;i=17, j=40history add id 5 fltime24、求解; --- reach new steady-state fluid flow ---set fluid time=0.0set fluid on mech onset mech force 0 ratio 1e-4set fluid substep 100set mech substep 50 auto ; mechanical cycles are slaved to fluid cycles solve age 5e925、删除单元体delete zone range x 0 10 z 38 4026、打锚杆; install horizontal strutsel beam id 2 begin 0 1 40 end 10 1 40 nseg 1sel beam id 2 property density 3000.0 emod 4.0E9 nu 0.30 & xcarea 1.0 xciy=0.083 xciz=0.083 xcj=0.0;list sel node pos range x -0.1 0.1 y 0.9 1.1 z 39.9 40.0;list sel node pos range x 9.9 10.1 y 0.9 1.1 z 39.9 40.0; beam node on plane of symmetry free to move in zsel node fix x y xr yr zr range id 154;list sel node fix range id 154sel join beam id 2 liner id 1; get the link id of the beam node that is connected to the liner;list sel node link range id 155;list sel link attach range id 308; set up a pin connectionsel link attach xrdir free range id 308sel link attach yrdir free range id 308sel link attach zrdir free range id 30827、组的定义; (wall is impermeable)group zone wall range x @_x1 @_x2 y 0 @_y2 z @_z1 @_z3group zone wall range x 0 @_x1 y @_y1 @_y2 z @_z1 @_z3prop poro 0 range group wall28、历史显示;disp contoursplot create plot disp_contplot set job onplot set viewtitle text 'Displacement contours [m]' onplot add contour displacement ramp greyscale ...minimum 0.0 ...maximum 1.0e-2 ...interval 1.0e-3 reversed onplot set orientation (58.295,231.78,-2.0495) ¢er (37.193,22.662,35.162) &eye (-66.413,-58.915,116.62) &radius 64.203;plot 20; historiesplot create plot disp_histplot set job onplot set viewtitle text 'Displacements histories near mid-height of the wall' on plot add his 100 linestyle color hsv (0,0,225) ...101 linestyle color hsv (0,0,150) ...102 linestyle color hsv (0,0,75) ...103 linestyle color hsv (0,0,0) ...xaxis min 0.0 ...label "Steps [-]" ...yaxis label "Displacement [m]"29、保存模型、组、结构体、位移云图等到图片plot post plot model filename "167F0005-model.ps" manualplot post plot groups filename "167F0006-layers.ps" manualplot post plot sel_geom filename "167F0019-sels.ps" manualplot post plot disp_cont filename "167F0014-disp5.ps" manualplot post plot disp_hist filename "167F0020-hist.ps" manualplot post plot tiebacks_forces filename "167F0010-forces5.ps" manual30、显示单元体在面的法线方向那侧plot zgroup plane frontplane31、显示塑性区plot zone colorby state32、显示单元体在plane的一侧按类型显示单元：plot zone colorby type;he geometric type of the zone (brick, tetra,wedge, etc.) plot zone plane b;Objects in back (opposite the normal(s)) ofthe current cutting plane are shown.。