- 积分
- 0
- 注册时间
- 2010-5-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2010-5-7 22:24:46
|
显示全部楼层
来自 安徽合肥
;****************************************
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 <keyword> x1 y1 <x2 y2> <x3 y3> . . .
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 <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 <keyword> value
cable
cable begin . . . end . . .
id
nseg
pretension
property
cable <id id> keyword <range>
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 <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
not
del range group 2 not '模型中只保留 group 2 ，其他的都删除
;****************************************
out(s) 在屏幕上显示s的字符信息，s必须是字符类型。 执行成功，返回0，否则为1.
in(s) 键盘输入函数。等待用户键盘输入
;**************************************** |
|
IP卡
狗仔卡
显身卡
。。。。。。。。。。。。。。
