flac3d50常用命令集锦Word格式文档下载.docx

flac3d50常用命令集锦Word格式文档下载.docx

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

setlogontruncate

setlogoff

listzoneprinrangex01y01z01;

显示指定范围内各单元的主应力，结果如下

Hist命令：

①命令编号按顺序从1开始：

eachhistoryisnumberedsequentiallyfrom1asitisenteredviatheHISTORYcommand.

②查找显示所有的his命令：

ReturntotheFlac3D>

promptandtype

listhist

foralistingofthehistoriesandtheircorrespondingnumbers.

histnstep5;

每5步记录1次。

默认是10步记录1次

histaddratio;

记录maximumunbalancedforceratio

histaddgpzdisp448;

记录节点的位移

modelnullrangex=2,4y=2,6z=5,10;

Toexcavatethetrench,enter

Thedimensionsoftheexcavationregion（i.e.,thenullmaterial）aredefinedbythex-,y-andz-coordinaterange.

setlarge;

设置大变形large-strainlogicisspecified

inixdis0ydis0zdis0;

位移归零

solve

savetrench;

保持计算结果的文件“TRENCH.F3SAV”不需要后缀

4、group命令（注意新版本的区别，groupzone）

genzonebricksize666

groupzoneTunnelrangex16y06z15;

设置group

莫尔库伦模型

modelmechnullrangegroupTunnel;

groupzoneTunnelrangecylinderend1100105100end210095100radius3

groupzoneMatrixrangegroupTunnelnot

;

namingRANGESthatarefixedinspaceandunrelated

tothemodelgeometry-representsstratigraphy...

rangenameLayer1planedip0dd0ori000above

rangenameLayer2planedip0dd0ori000below

propbulk1e8shear0.5e8coh0tens0fric35rangenrangeLayer1

propbulk1e7shear0.3e7coh1e7tens0fric0rangenrangeLayer2

plzonecolorbygroup

4、Interface命令

以下是切割模型法生成interface：

寺河模拟用的

plzone

genzonebrickp0000p120000p20400p3004size100202group1;

genzonebrickp0004p120004p20404p3008size100202group2;

genseparategroup1

interface1wrapfirstgroup1secondgroup2

interface1maxedge0.5

显示interface：

plotcreateplotPlot01

plotaddsurface

显示表面

plotaddinterfacered

显示交互面为红色

plotshow

用导来导去法，已经在5.0上运行过，很好，值得借鉴

genzonebricksize333&

p0（0,0,0）p1（3,0,0）p2（0,3,0）p3（0,0,1.5）&

p4（3,3,0）p5（0,3,1.5）p6（3,0,4.5）p7（3,3,4.5）

groupzoneBase

p0（0,0,2.5）p1（3,0,5.5）p2（0,3,2.5）p3（0,0,7）&

p4（3,3,5.5）p5（0,3,7）p6（3,0,7）p7（3,3,7）

groupzoneToprangegroupBasenot

interface1facerangeplanenorm（-1,0,1）origin（1.5,1.5,3）dist0.1

pladdinterface

pause

inizposadd-1.0rangegroupTop

4、hist命令（5.0版本需要加add！

）

histnstep=5//设置采用间隔为5等效于sethist_rep=5

histaddid=2zoneszz0,0,-1

histunbal//最大非平衡力

histaddgpzdisp448

printhist

histaddgpzdis0,0,0

histaddgpzvel0,0,0

histaddgpydisp0,0,0

监测网格点（0，0，0）点的y方向上的位移

histaddzonesyy0,1,0

监测单元（0，1，0）点的y方向上的应力

histaddzonesyy1,1,0

监测单元（1，1，0）点的y方向上的应力

histaddzoneszz0,0,-1

Listhistory

setlarge

historynstep2

historyaddfishtt

historyaddfishtheta

historyaddzonesxx（1,1,1）

historyaddzonesyy（1,1,1）

historyaddzoneszz（1,1,1）

plotaddhis345vs1

HISTORYaddfishvar;

fishvar变量的监测

5、solve命令?

?

什么意思

5、plot命令

通过图形操作界面的方式：

plzonecolorbygr;

查看各个分组

通过命令流的方式：

plotclear

plotaddhistory1linestylecolorblack

plotpostplotPlot01file127F1002.psmanual

plotcreateplotPlot02

plotaddhistory2linestylecolorblack

plotpostplotPlot02file127F1003.psmanual

plotcreateplotGravV

plotcutaddplanedip90dd0origin3,4.1,0;

截面

plotsetcenter333

plotseteye121812

plotaddzoneplanebehindplaneononplaneoff;

addlabel"

Zone"

white

plotaddbcontszzrampgrey

plotaddaxes

plotpostplotGravVfile127F1004.psmanual

reskw50

plotcreateplotszz_contour

plotcutaddplanedip90dd180origin0,18.9,0

plotsetorient1202200

plotsetcenter45450

plotaddbcontszzramprainbow

plsetmag1.0

plotcreateGravV

plotsetplanedip=90dd=0origin=3,4,0

plotsetrot15020;

thiswouldbeachievedinteractively

plotsetcenter2.54.24.0;

andisshowheresimplytoillustrate

plotaddboundbehind

plotaddbcontszzplane

plotaddaxes

这样我们就建立了一个叫“GravV”的图形，将它设为当前图层，对此图形设置一个倾角为90°

的面，（在x-y面上，约定z的负方向为下），倾向为0°

。

（在x-y面上，从y的正方向开始沿顺时针方向计量），过x-y面上一点（x=3,y=4,z=0）,在平面后加一个网格，垂直方向的应力Qzz，最后给模型加轴线，整体图形轮廓与添加图形相反，显示每个区域中心的应力计算值，每个区域的颜色直接对应基本区域的应力，当键入show,我们建立的图形被显示

PlotcreateDispcont

Pladdcontdispplanebehind

Pladdaxes

Plshow

plotsketch

绘出轮廓线

plotaddvel

添加矢量

plotaddhis12vs3

绘出历史记录n，actLoad，p_sol

cyc15000

plot

CreateTrench

Addsurfaceyellow

Addaxesblack

Show

plotcreateszz_contour

创建文件名为szz_contour（z方向的应力云图）

plotsetcent445

视图的中心坐标为（4，4，5）

plotsetrot20030

视图的旋转角度为（20，0，30）

plotsetmag1.0

视图的放大倍数为1.0

plotaddcontszzoutonshadeon

打开szz的等值线图

打开坐标轴线

显示窗口

不同的区域，材料赋值不同

plzonecolorbymodel

求解前对模型进行检查

网格检查（如前所述）

边界条件检查

速度约束条件plotgpfix

荷载条件plotfap

模型检查

模型赋值ploblockmodel

参数赋值ploblockprop***

初始应力检查（如后所述）

经常检查模型的响应

plozconszz（syy,sxx）应力场

ploconzdis（ydis,xdis）位移场

plo（add）zonecolorbystate屈服状态

plogpfix速度约束条件

plofap体力

plohist（unbal）1不平衡力（必需的）

plsetmovieactiveon

plsetmovieextensionjpg

plsetmovieinterval20

plsetmovieprefixsihekuangkaiwa

plsetmoviesizeixiy

plsetnamesihedacaigao

plbitmapview1filemyplot.png

plbitmapview1size1024768

plcreate<

viewid>

plcurrent<

plcutaddplaneddxxxdipxxx

plcutaddplanenamexxxnormalxxxoriginxxx

pldxf<

<

filefilename.dxf>

plexcel<

filefilename.xlsx>

plsvg<

filefilename.svg>

plsetbackgroundwhite;

设置背景色

5、网格的连接、对称……

不同尺寸的zone，网格用attach连接，粘起来！

6、plane命令

6、range命令

listrange

rangenametrenchx=-1,1y=0,4z=-2,2;

定义range

modelmechnullrangenrangetrench;

;

开挖excavatethetrench

7、gp命令

listgpposition;

查看节点编号、位置坐标

printgpdisprangeid105

message

（1）=string（gp_xpos（pnt））

message

（1）=message

（1）+'

'

+string（gp_zpos（pnt））

+string（abs（gp_zdisp（pnt）））

GPPNTfindgp（INT）

GPPNTgpfind（INT）

GPPNTgphead

ii=ngptotalnumberofgridpoints

更多节点变量函数，详见fishreference

8、zone命令

listzone

printzonestressrangeid2381anyid2361any2341any;

输出网格单元的应力

histzoneszz0,0,-1

ZPNTfindzone（INT）

INTzfind（INT）

ZPNTzonehead

ii=nzonetotalnumberofzones

更多zone变量函数，详见fishreference

9、config

9、边界条件BoundaryandInitialConditions

BoundaryconditionsarenormallyspecifiedwiththeAPPLYandFIX/FREEcommands,andinitialconditionsareassignedwiththeINITIALcommand

Example：

aconstant,compressivexx-stresscomponentof10MPacanbeappliedtoaboundarylocatedatx=10withthecommand

applysxx-10e6rangex10;

x=10的面上10MPa

applysyy=-20e6grad0,0,20e5rangey-20z0,10

在斜面上施加法向力

10杂项（材料力学性质）

11.M-C本构模型的选取

在进行模拟时，采用摩尔-库仑（Mohr-Coulomb）强度准则来判断矿体与顶底板的破坏机理。

该准则基本的内涵为：

岩石的剪切破坏是由岩石破裂面产生的破坏剪应力受到岩石材料的内聚力和内摩擦力的抵抗，当破裂面上的破坏剪应力大于岩石的内聚力和内摩擦力时，即Fs>

0时，材料发生剪切破坏[4]：

有限差分法（FDM）

Fish函数

1、if条件从句

2、loop条件从句

new

defderive

s_mod=y_mod/（2.0*（1.0+p_ratio））

b_mod=y_mod/（3.0*（1.0-2.0*p_ratio））

end

set@y_mod=5e8@p_ratio=0.25

@derive

list@b_mod@s_mod

genzonebricksize（2,2,2）

modelmechelastic

propbulk=@b_modshear=@s_mod

listzonepropbulk

listzonepropshear

defxxx

sum=0

prod=1

loopn（1,10）

sum=sum+n

prod=prod*n

end_loop

@xxx

list@sum@prod

defabc

ifxx>

0then

abc=33

else

abc=11

end_if

set@xx=1

list@abc

set@xx=-1

genzonebricksize1035

defplace_cables

loopn（1,5）

z_d=float（n）-0.5

command

selcablebeg（0.0,1.5,@z_d）end（7.0,1.5,@z_d）nseg7

end_command

@place_cables

plotaddzonetrans75

plotaddselgeom

modelmechmohr

propbulk1e8shear0.3e8fric35

propcoh1e3tens1e3

inidens1000

setgrav（0,0,-10）;

重力加速度

fixxyzrangez0.0

fixyrangey0.0

fixyrangey3.0

fixxrangex0.0

fixxrangex10.0

histaddunbal

savecab_str

inixdis0ydis0zdis0

histaddgpxdisp（0,1,5）

z_d=5.5-float（n）

z_t=z_d+0.5

z_b=z_d-0.5

freexrangex0.0z@z_b,@z_t

selcablebeg（0.0,0.5,@z_d）end（7.0,0.5,@z_d）nseg7

selcablebeg（0.0,2.5,@z_d）end（7.0,2.5,@z_d）nseg7

selcablepropemod2e10ytension1e8xcarea1.0...

gr_k2e10gr_coh1e10gr_per1.0

savecab_end

plotzonetrans75selcabblockforcelinewidth4

genzonebricksize121

propshear=1e8bulk=2e8...

cohes=1e5tens=1e10

fixxyzrangey0.0

applyyvel-1e-5rangey2.0

defget_ad

ad1=gp_near（0,2,0）

ad2=gp_near（1,2,0）

ad3=gp_near（0,2,1）

ad4=gp_near（1,2,1）

@get_ad

defload

load=g