FLAC3D常用命令.doc

1. apply（缩写：app）可用来定义边界条件及初始条件：1）添加应力格式1：apply szz -0.3395e6 range z -0.1 0.1格式1：apply szz -0.3395e6 range group pile格式3：apply szz -0.3395e6 range z -0.1 0.1 group pile格式4：apply nstress 数值 range z 2.9 3.1（或3） x 1 2 y 1 22）以一定速度施加位移边界格式1：apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step3）添加边界条件格式：apply szz grad range 示例1：apply szz -1e9 grad 0 0 8.3e5 range z 0 120示例2：apply szz -0.6e6 range z 0.05 0.15 group pile注：是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。2. range(缩写：ran)通过range功能，可以使命令作用在一定指定范围的目标上；如果一个命令没有使用range来确定范围，则命令对整个模型有效。1)利用坐标指定一定的范围格式1：range z 0 1格式2：range z 2.9 3.1 x 1 2 y 1 22)利用分组来指定范围格式：range group 13)以上两种的复合格式：range z -0.1 0.1 group pile4）利用id号来指定一定的范围格式：range id 0 10该命令后跟起始id和结束id，这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。5）命名一个范围（需要先命名这个范围）格式：range name 示例：range name intersected_zones x 5 8 y 3 7引用这个范围的格式：内容+range nrange 示例：range nrange intersected_zone6）命名一个柱形范围格式：同4）示例1：range name bhole cylinder end1 0 0 0 end2 -9.14 0 0 rad 0.46示例2：range name interface1 cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.36 & cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.34 not7）any和notAny跟在range命令后表示“和”的意思，即使若干个定义的范围同时起效；not跟在range后表示“除了”的意思，即除此之外的范围示例1：group wall range group 1 any group 2 any ;把组1和组2定义为新组，组名为wall示例2：group wall range group 1 not ;把除了组1之外的所有部分定义为新组，组名为wall.示例3：range name interface1 cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.36 & cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.34 not3. set用来设置重力、收敛标准、给某些变量赋值1）设置重力格式：set gravity 0 0 -102）设置收敛标准格式1：set mech ratio :用R无量纲值确定收敛标准（默认10e-5）格式2：set mech force ;绝对收敛标准，体系最大不平衡力小于某一临界值，此局限性较大。以上两种格式也可以进行耦合：示例：set mech force 0 ratio 1e-4格式3：自定义收敛标准，比如某一关键点的位移或速度达到某一临界值（遇到这种情况进行补充）3）在FLAC3D中给FISH语言中的变量赋值格式：set 变量=value 如：set hh=-224）设置大变形格式：set large 4）设置监测历史变量时的监测间隔格式1：set hist_rep ；默认的应该是10步一个间隔。示例：set hist_rep 1格式2 hist nstep 示例：hist nstep 54. history（缩写：hist）设置监测点要进行监测的项目，不设置编号的话会从1自动进行编号。监测点必须用id或坐标这两种方法之一指定，新版本后需加add，不加add其实也可以。1）监测整个模型的不平衡力格式：hist add id unbal2）设置节点监测项目格式1：hist add gp zdis 0 0 6（监测节点z方向位移） hist add gp zvel 0 0 6（监测节点z方向速度）格式2：hist add id gp zdis 或者监测点id示例1：hist add id 2 gp zdis 1 1 3示例2：hist add id 1 gp zdis id 1 ；监测1号节点z方向的位移，这个监测点编号是13）设置单元监测项目格式1：hist add zone szz 0 0 -1(监测节点zz方向应力)格式2：hist add id zone szz 或者监测单元id示例1：hist add id 7 zone szz 1.5 1.5 2.54）设置监测历史变量时的监测间隔格式1：set hist_rep ；默认的应该是10步一个间隔示例：set hist_rep 1格式2 hist nstep 示例：hist nstep 55）将计算结果写入到txt文件夹中格式：hist write skip begin end file 路径/文件名.txt示例1：hist write 7 vs 8 file C:UserszhangzhaoDesktop7-8.txt示例2：hist write 7 file C:UserszhangzhaoDesktop7.txt示例3：hist write 1 2 3 4 5 skip 20 begin 8300 end 8400 file C:UserszhangzhaoDesktop7.txt ;把1-5号监测点在8300到8400步之间的数据，按照每20步一个的记录间隔输出到7.txt文件中5.model用来指定模型的本构模型，以下格式中mech可以省略格式1：model mech mohr ;定义全模型的本构模型格式2：model mech mohr range group ；格式3：model mech null range group ；格式4：model mech null range z 0 16.property(缩写：pro)给模型材料赋予参数格式：pro bulk shear coh fric ten range 示例1：prop bulk 5.6420e6 shear 1.3789e6 coh 9.77e3 fric 15.36 range group 1示例2：prop bulk 6e6 she 2.2e6 fric 10 coh 5000 range z 17.3 20.17.group用来进行模型的分组1）建模时直接将这部分命名为一个分组格式1：gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 group 1格式2：gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 fill group 1 ；(对填充区域进行有别于外围材料的分组命名)格式3：gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 group 1 fill group 2 ；(外围区域命名为1，填充区域进行有别于外围材料的分组命名2)2）对一定范围内的数据命名分组格式：group 1 range ;用range定义范围3）分开命令为一个组示例：group dam range x 5 7 z -5 0group dam range id 201 any id 211 any id 231 any id 241 anygroup dam range id 202 any id 212 any id 222 any id 232 any 以上区域命名了一个分组dam。可以不同的区域命名为同一个组名，最后就会是一个组了。8.attach进行网格的连接，使用该命令时对各网格连接面上的单元尺寸有限制，要求它们之间的比率成整数倍格式：attach face range ;range后跟定义范围的关键词9.reflcet(缩写：ref)以某对称轴生成镜像网格格式1：gen zone ref ;以z=0平面为对称轴镜像生成网格格式2：gen zone ref dip 90 ; 以y=0平面为对称轴镜像生成网格格式3：gen zone ref dip 90 dd 90 ; 以x=0平面为对称轴镜像生成网格格式4：gen zone ref dip 90 dd 90 range z 20 35 ; 可以加上范围，在z为2035范围内以x=0平面为对称轴镜像生成网格。其中：dip为平面在空间坐标系下沿着z轴负向与xy平面所成的夹角；dd为平面的dip方向角，即坐标系的xy平面从y轴正向顺时针旋转到指定方向的夹角。如果不注明代表值为0。格式5：gen zone reflect normal 1 0 0 ori 0 0 0格式：normal x y z origin x y z:其中normal为指定平面的法向方向，origin为法线起点位置。这种方法比较直观，相比于用dd和dip定义平面，更容易记住，也不易出错。10.fix固定/free放松进行模型边界条件的设置1）约束x/y/z方向的位移，其实约束的是速度格式；fix range 示例1：fix z range z -0.1 0.1示例2：fix y range z -0.1 0.1示例3：fix x z range z -0.1 0.1 group 1示例4：fix x ；也可以不加范围，把所有的地方都固定2)fix pp 格式：fix pp 给定孔隙水压力边界条件3）给定温度边界条件格式：fix temperature 11.initial(简写ini)用来设置初始值1）用来设置密度格式：ini dens range 示例1：ini dens 1830 range group 1ceng示例2：ini dens 2000 range z 0 12）初始量归零格式：ini range 示例1：ini xdis 0 ydis 0 zdis 0 xvel 0 yvel 0 zvel 0示例2：ini szz 0 sxx 0 syy 03）初始状态为0格式：ini state 示例：ini state 0一般的程序地应力平衡后，初始化一般要进行塑性区、位移、速度清零ini state 0ini xdis 0 ydis 0 zdis 0ini xvel 0 yvel 0 zvel 04）设置初始孔隙水压力格式：ini pp grad range 示例：ini pp 30e3 grad 0 0 -10e3 range z 0 1注：是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。5）设置初始应力格式：ini szz grad range 示例：ini szz -35e3 grad 0 0 20e3 range z 0 1注：是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。6）将单元体移动格式：ini add range 示例：ini zposition add -6.0 range group pile12.save将计算状态进行保存，会保存在当前文件下1.常规格式：save 文件名.sav示例1：save 6-1.sav2.自动循环rest 8-4.savtable 1 name load_settlement ;可要可不要，不要会自动从1编号def add_loadp_gp=gp_near(2,1,3)loop n (1,5)app_load=n*(-1000e3)filename=7-6_add_step+string(n)+.savcommandapp nstress app_load range z 2.9 3.1 x 1 2 y 1 2solvesave filenameendcommandxtable(1,n)=app_loadytable(1,n)=gp_zdisp(p_gp)endloopendadd_loadsave 8-6.sav13.solve对模型进行计算1、格式1：solve ;直接计算，默认最大平衡应力比小于或等于1e-5时，计算结果获得稳定值2、格式2：solve ratio ； 设置最大不平衡应力比，计算达到这个条件即终止示例：solve ratio 1e-53、格式3：step ；设置计算终止步数4、渗流求解格式1：solve ；得到渗流平衡状态下的孔压场，渗流平衡的概念类似于力学平衡，表示节点上的不平衡流量比达到了程序设置的最小值。格式2：step ；在计算中完成相应的步数以后程序自动结束。格式3：solve age ；时间是流动时间，这个时间并不一定是真实的时间，只有当所有的流体参数都为真实时才是真实的时间。14.plot详见FLAC3D快速入门及简单实例第四章FLAC3D图形显示进行绘图，进行绘图前最好指定一个绘图文件格式：plot create view(或者plot) 1）绘制模型单元、建立的界面等格式1：plot zone示例1：plot zone trans 75 ;绘制透明度为75的单元体示例2：plot zone fill off ;单元体不进行填充示例3：plot zone trans 75 fill off ；以上两种的组合示例4：plot zone colorlist red示例5：plot zone colorlist red range nrange 示例6：plot interface colorlist red示例程序：newgen zone brick size 10 10 10range name intersected_zones x 5 8 y 3 7range name union_zones x 5 8 y 3 7plot create view intersected_zones ;（创建一个intersected_zones的绘图纸）plot add zone colorlist yellow trans 80plot add zone colorlist red range nrange intersected_zoneplot add axesplot create view union_zones ;(也可通过plot create plot union_zones)plot add zone colorlist yellow trans 80plot add zone colorlist red range nrange union_zonesplot add axes格式2：plot interface2）绘制云图（应力、变形，节点与单元的量都可绘制）格式1：plot zcon 示例1：plot zcon szz Plot zcon sxy格式2：plot con 格式3：plot zcon ramp grey ;灰度显示示例1：plot con zdis3）绘制块云图（一般只表示应力的量）格式：plot bcon 示例1：plot bcon szz示例2：plot bcon pro bulk4）绘制矢量图（最好结合1）内容带上网格）格式：plot dis；绘制变形矢量图（）示例：plot dis scale 0.1 ;scale后跟的数值越大，矢量图的箭头越大。5）绘制监测变量曲线图格式1：plot hist 示例：plot hist 2格式2：plot hist vs ；两个监测变量之间的关系示例：plot hist 7 vs 8 (历史变量7位横坐标、8位纵坐标)6）plot add 的应用格式：plot add +内容 ；把要绘制的图形绘制在前面的图形上示例：plot zone trans 75 fill off Plot add dis相当于在单元图形上绘制了变形矢量图7）按组显示格:1：plot zgroup range 示例：plot zgroup range z 0 5格式2：plot zone range group 1其他常见的plot命令Plot gpfixPlot gpPlot boundaryPlot zone trans 75Plot clearPlot interface colorlist red8）绘制表格数据格式1：plot table id style both格式2：plot table id style line格式3：plot table id style mark示例1：plot table 1 style both9）绘制孔压云图节点：plot con gpp单元：plot zcon zpp15.restore(简写为rest)调用保存好的sav文件，并在此状态的基础上进行下面的操作格式：rest 文件名.sav示例：rest 6-1.sav16.state单元状态1）显示单元状态格式：plot zone colorby state2）输出单元状态格式：print zone state17.print输出计算结果1）输出节点、单元数据格式1：print zone/gp ;输出单元数据格式2：print gp position格式3：print zone stress格式4：print gp dis格式5：print gp dis range id 517 any 533 any格式6：print gp pos range id 59 ;id号为59的节点坐标信息2）输出接触面的数据格式1：print interface ;接触面数据输出格式2：print interface stress ；输出接触面上的nstress/stress等格式3：print interface dis格式4：print interface position格式5：print interface prop3）输出监测变量格式1：print hist ;输出监测点信息格式2：2）将计算结果输入到指定的文件(也可参考4hist)格式：（set log on Set logfile C:UserszhangzhaoDesktop6-3.log Print zone stress Set log off）以上的程序将单元应力数据输出到6-3文件中3）输出表格数据格式：print table id示例：print table 118.plane创建切片，显示剖面图Oringin/normal和dip/dd在定义平面上的作用是等效的1）建立剖面最后先建立一个绘图纸格式：plot create view 2）创建剖面格式：plot cut add plane ori norm 示例：plot cut add plane ori 1.5 1.5 1.5 norm 0 1 03）绘制剖面上的内容格式：plot add (plane on) ;括号内可省略示例1：plot add zone fill off ;绘制剖面上的网格，建议先绘制这个。示例2：plot add zcon szz示例3：plot add dis3）绘制剖面外的内容格式：plot add plane b(f) b代表与norm反向一侧，f代表正向侧。示例1：plot add zone plane b示例2：plot add zone fill off plane f示例程序：plot create view plane4plot cut add plane ori 1.5 1.5 1.5 norm 0 1 0plot add zone fill off plane oplot add zcon szz plane onplot add zone plane bplot add zone fill off plane f19.输出图片将界面上显示的内容以图片形式显示出来1）输出ps照片格式：plot post view file 文件路径/名称.格式示例：plot post view plot04 file C:UserszhangzhaoDesktop3.bmp2）输出bitmap(位图)图片 注：常用格式bmp png格式：plot bitmap view file 文件路径/名称.格式示例： plot bitmap view plot04 file C:UserszhangzhaoDesktop7.bmp以上示例中的view plot04是可以省略的，但是输出的就为当前界面的内容。20.delete删除指定范围内的网格格式：delete zone range Delete和model null是不一样的，delete是将网格信息删除，model null是只是将网格赋予了空属性，即乘以一个很小的数，后面还可以通过model将这些属性恢复。21.interface(int)建立接触关系，详见FLAC3D快速入门及简单实例第9章9.8节、FLAC/FLAC3D基础与工程实例第6章1）移来移去法格式：Interface face range 示例1：interface 5 face range cylinder end1 0 0 0 end2 0 0 -15.1 radius 0.16 & cylinder end1 0 0 0 end2 0 0 -15.1 radius 0.14 not多加上0.1,本来是0-15示例2：interface 2 face range cylin end1 (2.,2.,-19.9) end2 (2.,2.,-20.1) radius .26多加上0.1,本来半径是0.252）切割模量法关键命令：gen separation group 1 Int 1 wrap group1 group2示例程序：newgen zone brick size 3 3 3group 1 range x 1 2 y 1 2 z 1 2group 2 range group 1 not gen seperate group 1 ；将group 1 与主体模型分开int 1 wrap 1 2 ；建立接触面plot interface3)导来导去法示例程序：newgen zone brick size 3 3 3group 2 range x 1 2 y 1 2 z 1 2group 1 range group 2 notsave 1.savdelete zone range group 1 interface 1 face range x 0.9 1.1 y 1 2 z 1 2interface 1 face range x 1.9 2.1 y 1 2 z 1 2interface 1 face range x 1 2 y 0.9 1.1 z 1 2interface 1 face range x 1 2 y 1.9 2.1 z 1 2interface 1 face range x 1 2 y 1 1 z 0.9 1.1interface 1 face range x 1 2 y 1 1 z 1.9 2.1save 2.savrest 1.savdelete zone range group 2expgrid C:UserszhangzhaoDesktopgroup2.flac3drest 2.savimpgrid C:UserszhangzhaoDesktopgroup2.flac3dplot interfaceplot add zone fill off4）设置界面参数格式：interface prop kn ks coh fric ten 示例：interface 1 prop kn 2.097e13 ks 2.097e13 coh 30000 fric 2022.expgrid与impgrid网格的导出与导入格式：expgrid 文件名.flac3dimpgrid 文件名.flac3d示例：expgrid C:UserszhangzhaoDesktopgroup2.flac3dimpgrid C:UserszhangzhaoDesktopgroup2.flac3d23.list1)打印节点、单元状态格式1：list zone或者print zone格式2：list gp position 或者 print gp position2）打印单元数据等格式1：list zone stress格式2：list zone prop bulk/shear3）打印有关节点数据格式1：list gp格式2：list gp position格式3：list gp dis4）打印有关fish函数格式1：list var ;如果是函数变量，计算函数，并输出数值格式2：list fish ;按其名称的字母顺序输出所有的当前值，可用来检查5）输出表格数据格式：list table id示例：list table 16）孔压的输出节点：print/list gp pp 单元：print/list zone pp后面可以加上范围24.generate详细参考FLAC3D快速入门及简单实例第40页、FLAC/FLAC3D基础与工程实例第5章25.pause/continue这是两个控制命令流（.txt文件）读入进程的命令。Pause的作用是让程序在此处暂停，可以在命令行中输入continue让程序继续运行。当命令流有问题时但却不知道在哪里的时候，可以通过在不同的位置插入pause来确定问题的所在。26.structure elements结构单元，建立方法详见FLAC/FLAC3D基础与工程实例第157页。1）创建线性结构单元（beam cable pile）线性结构单元是通过sel来创建的格式1：sel pile id 1 begin 0 0 0 end 0 0 10 nseg 4格式2：示例程序如下newsel node id 1 0 0 0sel node id 2 2 0 0sel node id 3 4 0 -1sel node id 4 5 0 -2sel beamsel id 1 cid 1 node 1 2sel beamsel id 1 cid 2 node 2 3sel beamsel id 1 cid 3 node 3 4plot sel geometryplot add ax2）创建面结构单元（shell geogrid liner）面结构单元也是sel来创建的格式：3）显示格式1：plot sel geometry ;将所有单元都显示出来格式2：plot sel pile ;显示某种特定的单元格式3：plot sel geometry shrink 0.9 marktype node colorby id ;这种格式只能应用在格式1geometry中，因为marktype只有在这种条件下有。格式4：plot sel pile shrink 0.9 line width 4 colorby id ;这种格式只能应用在格式2pile中，因为line width 4只有在这种条件下有。4）结构单元之间的连接通过sel link 命令进行设置在实体单元内部或表面创建结构单元时，程序会自动对所有结构节点创建了node-zone连接。5）结构单元的边界与初始条件结构单元所有的边界条件和初始条件都是通过结构节点来设置的，采用sel node 命令来实现，详见FLAC/FLAC3D基础与工程实例第162页下面6）7）8）内容详见FLAC3D快速入门及简单实例7.8节。6）结构单元的赋值问题格式：sel type peop ；type可以是beam/cable/pile/shell/geogrid/liner示例：sel pile prop xcarea 6e-3 xciz 200e-6 xciy 200e-6 xcj 07）监测变量详见FLAC3D快速入门及简单实例7.8节8）输出各个节点的变量格式1：list sel beam type ;type为具体的监测项目格式2：list sel cable type ;type为具体的监测项目格式1：list sel pile type ;type为具体的监测项目格式1：list sel recover type ;shell 单元，type为具体的监测项目格式1：list sel geogrid type ;type为具体的监测项目格式1：list sel liner type ;type为具体的监测项目27.config格式1：config dyn ;打开动力响应格式2：config fluid ;打开渗流模式28.显示某些赋值的量格式1：plot zcon/bcon prop bulk格式2：list zone prop bulk/shear29.call格式1：call find-add.f3fis格式2：注：1. FLAC3D及FISH语言中是不分大小写的，全部会被转换为小写；2. FLAC3D中，速度以位移除以时步表示；3. 在FLAC3D总不进行设置的量代表为0，不设置范围的代表全模型。常用的计算步骤（也可进行局部的调整，如开挖与回填问题）：1、画网格2、材料本构模型与材料参数（先本构，再参数）3、边界条件与初始条件4、计算（初始平衡与后续计算）5、后处理（处理数据、作图）可以用来检查编程正确性的命令1、attach2、gen merge3、plot bcon szz