FLAC3D常用命令Word文档下载推荐.docx

FLAC3D常用命令Word文档下载推荐.docx

《FLAC3D常用命令Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《FLAC3D常用命令Word文档下载推荐.docx（18页珍藏版）》请在冰豆网上搜索。

指定id为1到10的单元为各向同性弹性本构。

5）命名一个围（需要先命名这个围）

rangename<

自己起一个名字>

<

示例：

rangenameintersected_zonesx58y37

引用这个围的格式：

容+rangenrange<

自己起的那个名字>

rangenrangeintersected_zone

6）命名一个柱形围

同4）

rangenamebholecylinderend1000end2-9.1400rad0.46

rangenameinterface1cylinderend1000end200-0.9radius0.36&

cylinderend1000end200-0.9radius0.34not

7）any和not

Any跟在range命令后表示“和”的意思，即使若干个定义的围同时起效；

not跟在range后表示“除了”的意思，即除此之外的围

groupwallrangegroup1anygroup2any;

把组1和组2定义为新组，组名为wall

groupwallrangegroup1not;

把除了组1之外的所有部分定义为新组，组名为wall.

示例3：

3.set

用来设置重力、收敛标准、给某些变量赋值

1）设置重力

setgravity00-10

2）设置收敛标准

setmechratio<

value>

:

用R无量纲值确定收敛标准（默认10e-5）

setmechforce<

;

绝对收敛标准，体系最大不平衡力小于某一临界值，此局限性较大。

以上两种格式也可以进行耦合：

setmechforce0ratio1e-4

自定义收敛标准，比如某一关键点的位移或速度达到某一临界值（遇到这种情况进行补充）

3）在FLAC3D中给FISH语言中的变量赋值

set变量=value如：

sethh=-22

4）设置大变形

setlarge

4）设置监测历史变量时的监测间隔

sethist_rep<

数目>

；

默认的应该是10步一个间隔。

sethist_rep1

格式2histnstep<

histnstep5

4.history（缩写：

hist）

设置监测点要进行监测的项目，不设置编号的话会从1自动进行编号。

监测点必须用id或坐标这两种方法之一指定，新版本后需加add，不加add其实也可以。

1）监测整个模型的不平衡力

histaddid<

编号>

unbal

2）设置节点监测项目

histaddgpzdis006（监测节点z方向位移）

histaddgpzvel006（监测节点z方向速度）

给定该监测变量的编号>

gpzdis<

坐标>

或者监测点id

histaddid2gpzdis113

histaddid1gpzdisid1；

监测1号节点z方向的位移，这个监测点编号是1

3）设置单元监测项目

histaddzoneszz00-1（监测节点zz方向应力）

zoneszz<

或者监测单元id

histaddid7zoneszz1.51.52.5

默认的应该是10步一个间隔

5）将计算结果写入到txt文件夹中

histwrite<

监测变量的编号>

skip<

记录间隔>

begin<

监测开始的步数>

end<

监测结束步数>

file路径/文件名.txt

histwrite7vs8fileC:

\Users\zhangzhao\Desktop\7-8.txt

histwrite7fileC:

\Users\zhangzhao\Desktop\7.txt

histwrite12345skip20begin8300end8400fileC:

\Users\zhangzhao\Desktop\7.txt;

把1-5号监测点在8300到8400步之间的数据，按照每20步一个的记录间隔输出到7.txt文件中

5.model

用来指定模型的本构模型，以下格式中mech可以省略

modelmechmohr;

定义全模型的本构模型

modelmechmohrrangegroup<

组名>

；

modelmechnullrangegroup<

modelmechnullrangez01

6.property（缩写：

pro）

给模型材料赋予参数

probulk<

值>

shear<

coh<

fric<

ten<

propbulk5.6420e6shear1.3789e6coh9.77e3fric15.36rangegroup1

propbulk6e6she2.2e6fric10coh5000rangez17.320.1

7.group

用来进行模型的分组

1）建模时直接将这部分命名为一个分组

genzoneradcylinderp000-18p180-18p200-19.1p308-18size84815ratio1111.15group1

genzoneradcylinderp000-18p180-18p200-19.1p308-18size84815ratio1111.15fillgroup1；

（对填充区域进行有别于外围材料的分组命名）

genzoneradcylinderp000-18p180-18p200-19.1p308-18size84815ratio1111.15group1fillgroup2；

（外围区域命名为1，填充区域进行有别于外围材料的分组命名2）

2）对一定围的数据命名分组

group1range<

>

用range定义围

3）分开命令为一个组

groupdamrangex57z-50

groupdamrangeid201anyid211anyid231anyid241any

groupdamrangeid202anyid212anyid222anyid232any

以上区域命名了一个分组dam。

可以不同的区域命名为同一个组名，最后就会是一个组了。

8.attach

进行网格的连接，使用该命令时对各网格连接面上的单元尺寸有限制，要求它们之间的比率成整数倍

attachfacerange<

range后跟定义围的关键词

9.reflcet（缩写：

ref）

以某对称轴生成镜像网格

genzoneref;

以z=0平面为对称轴镜像生成网格

genzonerefdip90;

以y=0平面为对称轴镜像生成网格

genzonerefdip90dd90;

以x=0平面为对称轴镜像生成网格

genzonerefdip90dd90rangez2035;

可以加上围，在z为20～35围以x=0平面为对称轴镜像生成网格。

其中：

dip为平面在空间坐标系下沿着z轴负向与xy平面所成的夹角；

dd为平面的dip方向角，即坐标系的xy平面从y轴正向顺时针旋转到指定方向的夹角。

如果不注明代表值为0。

格式5：

genzonereflectnormal100ori000

normalxyzoriginxyz:

其中normal为指定平面的法向方向，origin为法线起点位置。

这种方法比较直观，相比于用dd和dip定义平面，更容易记住，也不易出错。

10.fix固定/free放松

进行模型边界条件的设置

1）约束x/y/z方向的位移，其实约束的是速度

格式；

fix<

限制变形的方向>

fixzrangez-0.10.1

fixyrangez-0.10.1

fixxzrangez-0.10.1group1

示例4：

fixx；

也可以不加围，把所有的地方都固定

2）fixpp

fixpp<

给定孔隙水压力边界条件

3）给定温度边界条件

fixtemperature<

11.initial（简写ini）

用来设置初始值

1）用来设置密度

inidens<

inidens1830rangegroup1ceng

inidens2000rangez01

2）初始量归零

ini<

量>

inixdis0ydis0zdis0xvel0yvel0zvel0

iniszz0sxx0syy0

3）初始状态为0

inistate<

状态>

inistate0

一般的程序地应力平衡后，初始化一般要进行塑性区、位移、速度清零

inixdis0ydis0zdis0

inixvel0yvel0zvel0

4）设置初始孔隙水压力

inipp<

初始值>

inipp30e3grad00-10e3rangez01

5）设置初始应力

iniszz<

iniszz-35e3grad0020e3rangez01

6）将单元体移动

方向>

add<

距离>

inizpositionadd-6.0rangegrouppile

12.save

将计算状态进行保存，会保存在当前文件下

1.常规

save文件名.sav

save6-1.sav

2.自动循环

rest8-4.sav

table1nameload_settlement;

可要可不要，不要会自动从1编号

defadd_load

p_gp=gp_near（2,1,3）

loopn（1,5）

app_load=n*（-1000e3）

filename='

7-6_add_step'

+string（n）+'

.sav'

command

appnstressapp_loadrangez2.93.1x12y12

solve

savefilename

endcommand

xtable（1,n）=app_load

ytable（1,n）=gp_zdisp（p_gp）

endloop

end

add_load

save8-6.sav

13.solve

对模型进行计算

1、格式1：

solve;

直接计算，默认最大平衡应力比小于或等于1e-5时，计算结果获得稳定值

2、格式2：

solveratio<

计算精度值>

设置最大不平衡应力比，计算达到这个条件即终止

solveratio1e-5

3、格式3：

step<

步数>

设置计算终止步数

4、渗流求解

solve；

得到渗流平衡状态下的孔压场，渗流平衡的概念类似于力学平衡，表示节点上的不平衡流量比达到了程序设置的最小值。

在计算中完成相应的步数以后程序自动结束。

solveage<

时间>

时间是流动时间，这个时间并不一定是真实的时间，只有当所有的流体参数都为真实时才是真实的时间。

14.plot

详见《FLAC3D快速入门及简单实例》第四章FLAC3D图形显示

进行绘图，进行绘图前最好指定一个绘图文件

plotcreateview（或者plot）<

名字>

1）绘制模型单元、建立的界面等

plotzone

plotzonetrans75;

绘制透明度为75的单元体

plotzonefilloff;

单元体不进行填充

plotzonetrans75filloff；

以上两种的组合

plotzonecolorlistred

示例5：

plotzonecolorlistredrangenrange<

围名称>

示例6：

plotinterfacecolorlistred

示例程序：

new

genzonebricksize101010

rangenameunion_zonesx58y37

plotcreateviewintersected_zones;

（创建一个intersected_zones的绘图纸）

plotaddzonecolorlistyellowtrans80

plotaddzonecolorlistredrangenrangeintersected_zone

plotaddaxes

plotcreateviewunion_zones;

（也可通过plotcreateplotunion_zones）

plotaddzonecolorlistredrangenrangeunion_zones

plotinterface

2）绘制云图（应力、变形，节点与单元的量都可绘制）

plotzcon<

表示应力的量>

plotzconszz

Plotzconsxy

plotcon<

表示变形的量>

表示的量>

rampgrey;

灰度显示

plotconzdis

3）绘制块云图（一般只表示应力的量）

plotbcon<

plotbconszz

plotbconprobulk

4）绘制矢量图（最好结合1）容带上网格）

plotdis；

绘制变形矢量图（）

plotdisscale0.1;

scale后跟的数值越大，矢量图的箭头越大。

5）绘制监测变量曲线图

plothist<

给监测变量设置的id编号>

plothist2

编号1>

vs<

编号2>

两个监测变量之间的关系

plothist7vs8（历史变量7位横坐标、8位纵坐标）

6）plotadd的应用

plotadd+容；

把要绘制的图形绘制在前面的图形上

plotzonetrans75filloff

Plotadddis

相当于在单元图形上绘制了变形矢量图

7）按组显示

格:

1：

plotzgrouprange<

plotzgrouprangez05

plotzonerangegroup1

其他常见的plot命令

Plotgpfix

Plotgp

Plotboundary

Plotzonetrans75

Plotclear

Plotinterfacecolorlistred

8）绘制表格数据

plottableidstyleboth

plottableidstyleline

plottableidstylemark

plottable1styleboth

9）绘制孔压云图

节点：

plotcongpp

单元：

plotzconzpp

15.restore（简写为rest）

调用保存好的sav文件，并在此状态的基础上进行下面的操作

rest文件名.sav

rest6-1.sav

16.state

单元状态

1）显示单元状态

plotzonecolorbystate

2）输出单元状态

printzonestate

17.print

输出计算结果

1）输出节点、单元数据

printzone/gp;

输出单元数据

printgpposition

printzonestress

printgpdis

printgpdisrangeid517any533any

格式6：

printgpposrangeid59;

id号为59的节点坐标信息

2）输出接触面的数据

printinterface;

接触面数据输出

printinterface<

stress；

输出接触面上的nstress/stress等

dis

position

prop

3）输出监测变量

printhist;

输出监测点信息

2）将计算结果输入到指定的文件（也可参考4hist）

（setlogon

SetlogfileC:

\Users\zhangzhao\Desktop\6-3.log

Printzonestress

Setlogoff）

以上的程序将单元应力数据输出到6-3文件中

3）输出表格数据

printtableid

printtable1

18.plane

创建切片，显示剖面图

Oringin/normal和dip/dd在定义平面上的作用是等效的

1）建立剖面

最后先建立一个绘图纸

plotcreateview<

绘图纸名字>

2）创建剖面

plotcutaddplaneori<

norm<

方向矢量>

plotcutaddplaneori1.51.51.5norm010

3）绘制剖面上的容

plotadd<

容>

（planeon）;

括号可省略

plotaddzonefilloff;

绘制剖面上的网格，建议先绘制这个。

plotaddzconszz

plotadddis

3）绘制剖面外的容

planeb（f）b代表与norm反向一侧，f代表正向侧。

plotaddzoneplaneb

plotaddzonefilloffplanef

plotcreateviewplane4

plotaddzonefilloffplaneo

plotaddzconszzplaneon

19.输出图片

将界面上显示的容以图片形式显示出来

1）输出ps照片

plotpostview<

要输出的绘图纸名称>

file文件路径/名称.格式

plotpostviewplot04fileC:

\Users\zhangzhao\Desktop\3.bmp

2）输出bitmap（位