flac命令流集锦

;****************************************

group命令

group soil range z 1 2 ;定义group

prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值

model null range group soil

model elastic range group soil

plot block group range group dam ;只显示dam

range命令

range name trench x 0 1 y 0 4 z 0 2

model null range trench

model null range x=2,4 y=2,6 z=5,10

range name Big_Brick x -3 3 y -2 2 z -1 1

model elastic range Big_Brick

prop bulk 1e8 shear 1e8 range Big_Brick

range name Layer1 plane dip 0 dd 0 ori 0 0 0 above

range name Layer2 plane dip 0 dd 0 ori 0 0 0 below

range cylinder end1 x1 y1 z1 end2 x2 y2 z2 radius r

cylindrical range with one end of the cylinder axis (end1) at location (x1, y1, z1), the other end (end2) at location (x2, y2, z2), and with a cylinder radius of r '由(x1, y1, z1) 、(x2, y2, z2)两点确定旋转轴

macro命令

macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35'

macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0'

prop sand range Layer1

prop clay range layer2

macro Pt0 'p0 0 0 0'

macro Pt1 'p1 add 10 0 0'

macro Pt2 'p2 add 0 10 0'

macro Pt3 'p3 add 0 0 10'

macro Model_Size 'size 4 5 6'

macro Big_Brick 'gen zone brick Pt0 Pt1 Pt2 Pt3 Model_Size'

Big_Brick

macro 'Pt0' 'p0 15 15 15'

gen Big_Brick ; this will cause an error

group range macro 命令

Using different object types to do the same job

using a RANGE object ...

range name Big_Brick x = (-3,3) y = (-2,2) z = (-1,1)

model null range Big_Brick

using a GROUP object ...

group Big_Brick range x = (-3,3) y = (-2,2) z = (-1,1)

model null range group Big_Brick

using a macro object ...

macro Big_Brick 'x = (-3,3) y = (-2,2) z = (-1,1)'

model null range Big_Brick

gen zone brick size 10,10,10

macro SiltySand 'bulk 1.5e8 shear 0.3e8'

macro ClayeyGravel 'bulk 1.5e8 shear 0.6e8 fric 30 coh 5e6 ten 8.66e6'

model elas range z 5 10

prop SiltySand range z 5 10

model mohr range z 0 5

prop ClayeyGravel range z 0 5

;****************************************

视图操作：x y z m

ctr+r 还原

ctr+p 保存图片

ctr+z 鼠标选择

ctr+g 彩图变为灰色图

shift + （x y z m）旋转缩放

edit---copy to clipboard----粘贴到word

;****************************************

ini xdis 0 ydis 0 zdis 0

set grav 0 0 -9.81

set mech force=50 '最大不平衡力小于50n，停止计算

;****************************************

hist命令

hist reset '清除已有历史信息

set hist_rep 1 '每一个时步记录一次

hist n=5 '每5个时步记录一次

hist write 1 3 vs 2 begin 150 end 375 file xx.txt

hist gp zdisp 4,4,8

hist id=3 gp zdis 1 1 3 'id为监测变量的编号，默认的是从1开始编号。建议对监测变量进行编号，以便后处理调用

plot his 2 3 vs 4 'plots histories 2 and 3 versus history 4; history 4 plots along the abscissa.

pl his xla 'string' ylab 'string' '设置x y轴的名称

pl his xmaximum ?? xminimum ?? '设置x轴的最大最小刻度

pl his ymaximum ?? yminimum ?? '设置y轴的最大最小刻度

pl his both 'line+mark

;****************************************

log命令

set logfile xx.log

set log on

print ……

……

set log off

;****************************************

切片

pl set plane ori 0 1.5 0 norm 0 1 0 '定义剖面位置

pl con zd plane '显示变量

pl add ske

pl add dis plane

pl add axe

pl con sxx outline on ;查看水平应力云图，outline on表示显示模型的边界，默认为off ;****************************************

print gp dis range id 517 any id 533 any '输出2个节点的变形值any表示‘并集’

; 指定边界条件

fix x range x -0.1 0.1 any x 5.9 6.1 any '这个表示在x=0和x=6.0方向上固定x方向的位移。any的意思是表示"且"的意思，也就是说本来要写两行语句的：fix x range x -0.1 0.1，fix x range x 5.9 6.1，现在用any就可以只写一条语句。

fix y range y -0.1 0.1 any y 5.9 6.1 any

;****************************************

;****************************************

fish 主要语句

if ………… then

…………

else

…………

endif

command

…………

endcommand

loop n（1,21）

…………

endloop

loop while …………

…………

endloop

caseof …………

…………

case n1

…………

case n2

…………

endcase

;****************************************

运算

= # > < >= <=

负数在运算时，要加（），以免和减号- 混淆

degrad 'π/180

pi 'π

ngp '节点总数

nzone '单元体总数

;****************************************

单元/节点遍历

p_z = zone_head

loop while p_z # null

…………

p_z = z_next(p_z)

endloop

p_gp = gp_head

loop while p_gp # null

…………

p_gp = gp_next(p_gp)

endloop

;****************************************

string连接（保存文件）

save_file=string(n)+'_step.sav'

D3DECFile = 'D3dec_Model.dat'

FlacFile = 'Flac3D_Model.dat'

file_name='7-6_add'+string(n)+'.sav'

save file_name

cal flacfile

;****************************************

table操作详见ftd128（59/114)）

table 1 name load_settlement '创建新表格

xtable(n,s)=……对编号为n的表的第s行、x列进行赋值ytable(n,s)=……对编号为n的表的第s行、y列进行赋值plot table 1 line

plot table 1 both '点、线同时显示

常用编辑文本文件的方法进行表的读入与调用

第1行：表的名称

第2行：x1 y1

第3行：x2 y2

……

空行' 注：在表的文本文件最后，需要有一个回车换行符，否则会出现"error reading file xxx.dat"的错误，

完成的文件需要读入操作才可以供flac3d调用。采用table read 命令进行读入：

table 1 read xxx.dat

读入后，可以使用plot table 和print table命令查看生成的表文件

apply xvel 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2 '1.0表示表格中的y向数据的因子

TABLE n x1 y1 . . .

erase erases all entries in table n.

insert

name 'string' 'changes the name of table number n to 'string'. The table ID number is not changed.

read filename 'reads file, filename (in the format described below), and places it in table n.

;****************************************

将FLAC中节点以及节点对应的坐标全部输出到txt文件中？以及将每个节点的应力输出set log on

print gp pos

set log off

;****************************************

open close 用法

;close 当前文件关闭，返回值为0表示成功关闭

open(filename,wr,mode) '先打开文件，然后可以读，可以往里写内容

;===filename可以是带双引号的字符串，也可以是变量名。

;===wr必须是个整数：0-只读打开，文件要求存在；1-文件为写入打开，已经存在的文件将被覆盖（overwritten）。2-文件为写入打开，存在的文件被改写(appended to)

;===mode必须为一个整数值，0-读写fish变量，只有整型，浮点型和字符串的数据类型被传输，不传递变量名。1-读写ASCII数据。可读入一行数据，行之间用CR/LF分开。每行最大80字符。

;===2-设定为binary读模式，任何文件都以binary的读取模式打开.

;===返回值：0表示文件打开成功；1-表示文件名不是字符串；2-文件名是一个空字符串；3-wr或mode不是整型；4-mode参数错误（不是0或1）；5-wr参数错误（不是0或1）；6-不能打开所要读取的文件，比如文件不存在；7-文件已经打开；8-不是fish模式的文件

;read(ar,n) ---放入数组。读取n个记录放入数组ar，每个记录要求是一行ASCII数据，或者是一个FISH变量。数组至少有n个元素大小。

;====返回值0表示无误；-1：表示读取错误，比如到了文件的末尾；n-表示读入n行后，恰好到了文件的末尾。

;****************************************

read(ar, n) 将n个记录写到数组ar ，详细：

read(ar, n) reads n records into the array ar. The array ar must be an array of at least n elements. The returned value is:

0 requested number of lines were input without error

-1 error on read (except end-of-file)

n positive value indicates that end-of-file was encountered after

reading n lines

write(ar, n) 将数组ar的前n个记录写到文件，详细：

write(ar, n) writes n records from the first n elements of the array ar. The array ar must be an array of at least n elements. The returned value is:

0 requested number of lines were output without error

-1 error on write

n positive value (in ASCII mode) indicates that the nth element was not a string parse(s, i) This function scans the string s and decodes the ith item

pre_parse(s, i) This function scans the string s and returns an integer value according to the type of the ith item, as follows.

0 missing item

1 integer

2 float

3 string missing (unable to interpret as int or float)

type(e) 数据类型函数，1=整型，2=浮点，3=字符串，4=指针，5=数组

;****************************************

Interface 的变量命令参见ftd128.pdf（53/114）

interface 1 face range ^^^^^^^^^^^^^^^^^^^^^^^^^^

Interface 变量

Interfaces may be identified with the following functions.

i_find(id)————address of the interface with ID id; returns null if not found (pointer)

i_id(p_i) ————ID of the interface at address p_i (integer)

The interfaces may also be scanned starting at i_head and stepping through the list using i_next.

i_head————address of the first interface in the list of interfaces (pointer)

i_next(p_i)————address of the pointer to the next interface in the list (pointer)

wrap

For example, the following command would find the twinned faces between group "rock" and group "soil," and put interface elements on these "rock" faces. Only faces with centroid within the range x 50.0 75.0 would be considered.

interface 1 wrap rock soil range x 50.0 75.0

Structural Elements 的变量命令参见ftd128.pdf（64/114）

;****************************************

nd_pos( np, p, dof ) position p (p ∈{1,2} denotes current or reference position, respectively; dof-component, dof ∈{1,2,3}). The reference position is the configuration for which stiffness matrices have been formed and does not change during a small-strain analysis. The current position is updated after each timestep. During a large-strain analysis, the reference position is set equal to the current position during each large-strain update.

例子：xt0 = nd_pos(_nd, 2, 1)

yt0 = nd_pos(_nd, 2, 2)

zt0 = nd_pos(_nd, 2, 3)

nd_id(np) ID number of node np. Each node has a unique ID number.

;****************************************

详细的gp zone 命令参见ftd128.pdf（40/114）

p_gp =gp_near(x,y,z) 'address of gridpoint closest to (x, y, z)

p_z =z_near(x,y,z) 'address of zone closest to (x, y, z)

gp_head

zone_head

gp_next(p_gp)

z_next(p_z)

p_gp=find_gp(id) 'address of gridpoints with ID number id

p_z=find_zone(id) 'address of zone with ID number id

gp_group(p_gp, ind) ????

gp_region(p_gp, ind) ????

gp_zdisp(p_gp):地址为p_gp的节点的z向变形

gp_xpos(p_gp) 'x-coordinate of gridpoint px=gp_xpos(p_zp) gp_xpos(p_zp)=k*px gp_xvel(p_gp) 'x-velocity at gridpoint

gp_near(x,y,z) '得到靠近坐标(x,y,z)的节点地址

gp_id(p_gp) '节点ID号,是整数maxdisp=gp_id(p_gp)

print gp pos ran id 59 '输出59号节点的坐标位置(2.0,2.0,3.0)

注意：节点指针的循环是从ID号最小的地址开始的。

zone 函数

z_sig3(p_z) '最小主应力z_sig1(p_z) '最大主应力

z_prop(p_z,'young') 'z_prop(p_z,'young')=2e6

z_group(p_z) 'zone group name 是单元所在的组名变量(string)

z_id(p_z) 'zone ID number (integer)

z_model(p_z) ' string类型

z_prop(p_z,string) '可以赋值，可以取值

z_xcen(p_z) 'zone的质心的x坐标

z_sxx(p_z) 'zone的xx-stress

stid=2131

z_p=find_zone(stid)

xtable(3,m)=z_id(z_p)

ytable(3,m)=-z_szz(z_p)

;****************************************

def install

pnt = zone_head

loop while pnt # null

z_depth = -z_zcen(pnt)

y_mod = y_zero + cc * sqrt(z_depth)

z_prop(pnt, 'shear') = y_mod / (2.0*(1.0+P_ratio))

z_prop(pnt, 'bulk') = y_mod / (3.0*(1.0-2.0*P_ratio)) pnt = z_next(pnt)

end_loop

end

set p_ratio=0.25 y_zero=1e7 cc=1e8

install

plot block prop bulk

;****************************************

用户自定义的节点、单元额外变量

config zextra n '首先进行单元额外变量的配置config gpextra m '首先进行节点额外变量的配置

gp_extra(p_gp,n)

z_extra(p_z,n)

;--- test of stress rotation ---

config gpextra 2

call fishcall.fis

gen zone brick size 1 1 1

mo el

prop she 300 bu 300

def ini_coord

pnt = gp_head

loop while pnt # null

gp_extra(pnt,1) = sqrt((gp_xpos(pnt)-xc)^ 2+(gp_zpos(pnt)-zc)^ 2)

gp_extra(pnt,2) = atan2((gp_xpos(pnt)-xc),(zc-gp_zpos(pnt)))

pnt = gp_next(pnt)

endloop

end

set xc=0 zc=0

ini_coord

;****************************************

title命令

title 'Tutorial Examples'

plot set title text 'Mesh for trench example'

;****************************************

应力梯度

ini dens 2000

ini szz -40e3 grad 0 0 20e3 ran z 0 2

ini syy -20e3 grad 0 0 10e3 ran z 0 2

ini sxx -20e3 grad 0 0 10e3 ran z 0 2

set grav 0 0 -10

ini xdis=0 ydis=0 zdis=0

hist gp xdisp 1,0,0

hist gp zdisp 0,0,2

apply syy = -20e6 grad 0,0,20e5 range y -20.1,-19.9 z 0,10

With this command, a yy-stress component is applied to boundary faces located between y = -20.1

and -19.9. The yy-stress varies linearly with z from σyy = 0 at z = 10 to σyy = -20e6 at z = 0.

S=S(0)+gx*x+gy*y+gz*z ' gx gy gz 分别对应grad 0,0,20e5 中的三个分量

;****************************************

基本形状网格

gen zone brick size 6,8,8 p0 -10, -10, -20 &

p1 10, -10, -20 &

p2 -10, 10, -20 &

p3 -10, -10, 0

plot surf

gen 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, 10

plot surf

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 fill 'The fill keyword fills the brick region with zones plot surf

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 fill

gen zone reflect dip 0 dd 90

gen zone reflect dip 90 dd 90

gen zone reflect normal (xv yv zv) origin (xv yv zv)

;

; identify the trench

gen zone radc &

dim 3 3 3 3 &

ratio 1 1 1 1.2 &

size 3 8 8 5 &

edge 10 &

p0 100 95 100 &

fill

gen zone reflect dip 90 dd 90 ori 100 100 100

gen zone reflect dip 0 dd 0 ori 100 100 100

;****************************************

cable

sel cable beg 1.0,0.4,1.5 end 5.0,0.4,1.5 nseg 4

sel cable beg 1.0,0.4,0.5 end 5.0,0.4,0.5 nseg 4

sel cable beg 1.0,1.2,1.5 end 5.0,1.2,1.5 nseg 4

sel cable beg 1.0,1.2,0.5 end 5.0,1.2,0.5 nseg 4

sel cable prop emod 2e9 ytens 1e8 xcarea 1.0 &

gr_coh 1e10 gr_k 2e9 gr_per 1.0

;****************************************

plot create GravV

plot set plane dip=90 dd=0 origin=3,4,0

plot set rot 15 0 20

plot set center 2.5 4.2 4.0

plot add bound behind

plot add bcont szz plane

plot add axes

plot show

;****************************************

plot hold hist 1 vs -2 'plot hist m vs n中m代表y轴；n代表x轴；负值是对其值的镜像pl sk magf 20 'sk 显示模型的外围网格线

pl con szz ou on magf 20

;****************************************

if x_pos = 10.0 then

new = get_mem(2)

mem(new) = head

mem(new+1) = p_gp

head = new

endif

**************************************

举个例子

def find_add ;定义fish函数find_add

head = null ;给head赋值

p_gp = gp_head ;第一个网格结点的指针赋给p_gp

loop_while p_gp # null ;当p_gp值不为null时作循环

x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_pos

if x_pos = 10.0 then ;如果x_pos = 10.0 则(执行)

new = get_mem(2) ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址mem(new) = head ;将head类型和数值置于地址为new的fish变量

mem(new+1) =p_gp

head = new ;将new值赋给head

endif

p_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gp

endloop;结束循环

end ;结束fish函数

实际上这个fish函数为满足条件（x坐标为10的）的zone的地址开辟一定的地址空间，各地址之间存在一定的联系;找到第一个符合条件的zone地址后，用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址，并用下面的mem(new) = head ，将第一次开辟的两个变量的第一个变量存储地址head（注意第一次head＝null），第二个变量存储第一个符合条件的zone地址，并将第一个变量的地址赋予head（head＝new），第一次循环结束；下次循环，同样开辟两个变量对象，第一个变量对象记录上次循环开辟的第一个变量的地址，第二个记录第二个符合条件的zone地址，其余循环依次类推，这样子就建立了一个符合条件的zone地址链条，方便以后使用。

在调用时，用ad=head，就将最后一个循环开辟的两个变量的第一个变量的地址赋予ad，进行调用时，后找到的zone地址将被先调用，最后一直循环到最先开辟的两变量，因最先开辟的两变量的第一个变量的地址为null，因此可以控制循环结束。

;****************************************

ini state 0 ' 是把单元初始化为弹性状态，在施加初应力时可能使得某一部分单元进入了塑性状态。因此把它改回弹性。

;****************************************

id,cid是什么意思？

id是指在整个结构中的编号，而cid是指在某一类比如说cable中的编号。拿cable 中的一个单元来说，它既有自己在整个结构中的cd,又有自己在cable中的cid

;****************************************

movie命令

rest 6-1.sav

ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0

app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2

plot set rot 20 0 30

plot con szz ou on magf 10

plot add hist 1

set movie avi step 1 file 6-5.avi

movie start

solve

movie finish

;****************************************

CONFIG keyword

The following keywords apply.

cppudm C++ user-defined models (only available with the C++ user-defined model option)

creep creep material analysis (only available with creep model option)

dynamic fully dynamic analysis (only available with dynamic model option)

fluid fluid-flow analysis (see see Section 1 in Fluid-Mechanical Interaction (FLAC3D Manual))

gpextra n

n extra gridpoint variables for FISH use (see FISH REFERENCE)

thermal thermal analysis (only available with thermal model option)

zextra n

n extra zone variables for FISH use (see FISH REFERENCE)

;****************************************

结构单元

SEL keyword value

cable

cable begin . . . end . . .

id

nseg

pretension

property

cable keyword

pretension value

applies given pre-tension force to all cableSELs in the range and with an ID number of id; if id is not given, then all cableSELs in the range are considered. A positive pre-tension force places a cableSEL into tension. The given force is added to the current force being carried by each cableSEL.

property keyword value . . .

assigns the specified property to all cableSELs in the range and with an ID number of id; if id is not given, then all cableSELs in the range are considered. The following properties are available.

densitydensity (needed if dynamic mode or gravity is active)

emod Young's modulus

gr_coh grout cohesive strength (force) per unit length

gr_fric grout friction angle (degrees)

gr_k grout stiffness per unit length

gr_per grout exposed perimeter

slide large-strain sliding flag

slide_tol large-strain sliding tolerance

thexp thermal expansion coefficient

xcarea cross-sectional area

ycomp compressive yield strength (force)

ytens tensile yield strength (force)

;****************************************

dim 内部区域的尺寸edge fill

gen zone radcylinder size 5,10,10,5 &

p0 (0,0,0) p1 add (400,0,0) p2 add (0,200,0) p3 (0,0,200) &

dim (20,20,20,20) fill

del range group 2 not '模型中只保留group 2 ，其他的都删除

;****************************************

out(s) 在屏幕上显示s的字符信息，s必须是字符类型。执行成功，返回0，否则为1.

in(s) 键盘输入函数。等待用户键盘输入

;****************************************