flac2d命令汉语版.doc

flac2d命令汉语版.doc

《flac2d命令汉语版.doc》由会员分享，可在线阅读，更多相关《flac2d命令汉语版.doc（15页珍藏版）》请在冰豆网上搜索。

FLAC程序使用手册

FLAC输入命令

FLAC的输入和一般的数值模拟的程序不一样,它可以用交互的方式从键盘输入各个命令,也可以写成命令文件,类似于批处理,由文件来驱动。

FLAC命令大小写一样。

所有的命令可以附带若干个关键词和有关的数值。

在下面的命令解释中,只有大写的字母起作用,小写的字母写不写、写多少个都没有崐关系。

i,j,m和n开始的变量要求整型数,否则要求实型数。

实型数的小数点可崐以忽略,但是整型数不能带小数点。

数值间可以用空格隔开,空格的数目不限,也可以用下面的分隔符隔开:

（）,/=

<>表示可选的参数,输入时括号不用输入;

...表示可以有任意个参数。

由*号开始到行末为注释,FLAC在执行时不理会。

下面的FLAC命令按字母排列。

Apply关键词=数值<关键词=数值...><范围>

可以有下面的关键词:

Pressure压力

XForceX-方向的力

YForceY-方向的力

ATtach该命令可以将一条线上的结点和另一条线上的结点互相接合在一起,用以形成复杂的网格形状。

Call文件名

写成的命令文件可以用Call命令来调用,命令文件的最后一行必须是RETURN,以返回到交互方式。

命令文件中不能有CALL命令本身。

Config关键词

FLAC用以解平面应变问题,但经过配置命令也可以用于解平面应力问题或轴对称问题。

需要时应在形成网格之前发。

关键词有:

P_STR平面应力问题

AX轴对称问题

CYCn

该命令同STEP,为执行n个时步的循环运算。

FixX<范围>

Y

XY

用此命令可以使<范围>内结点的x-或y-方向的速度保持不变。

<范围>的格式可以是I=i1,i2,J=j1,j2;i和j何者先输入没有关系。

如果要求位移固定,则必须将速度初始化为零（开始时速度的缺省值为零）。

FIX和INIXV,YV可以联合使用以提供一个刚性移动的边界条件。

如果有MARK的关键词,则只有在此范围内被标记的结点将被固定。

FReeX<范围>

Y

XY

该命令与FIX命令相反,用于放松对结点的约束。

GEnx1,y1x2,y2x3,y3x4,y4

用于在全域或局域中产生网格。

（x1,y1）的点放在左下侧,（x2,y2）、

（x3,y3）和（x4,y4）为顺时针排列。

如果点和前一个GEn命令中的相同则用SAME来代替。

如果网格不是均匀排列,则用关键词RATIO,在I和J方向的比例系数分别为ri和rj.

GEnCirclexc,ycrad

Arcxc,ycxb,ybtheta

Linex1,y1x2,y2

该命令产生园、弧或直线。

ARC弧的圆心为（xc,yc）,起始点为（xb,yb）,逆时针的角为theta度

CIRCLE园的中心为（xc,yc）,半径为rad

LINE直线的两端为（x1,y1）和（x2,y2）.

用ARC等关键词所形成的边界系被标记（参见命令MARK和UNMARK）过,标记结点所围住的区域在以后用命令INITAL,MODEL和PROPERTY时可以用关键词REGION来说明。

GenADJust

用此命令来微调网格,使之离散平滑。

该命令可以连续使用以增进平滑效果的。

Gridicoljrow

此命令用以产生数目为icol列和jrow行的网格。

对于640k内存的FLAC版本,如用摩尔-库仑的本构模型约可分两千个单元。

Help帮助命令,可在屏幕上显示命令表。

His<关键词...I=i1J=j1>

每隔NSTEP时步,记录一下关键词所示项目在i1和ji结点或单元的数值。

NSTEP的缺省值为10.用户应记住所要求His（历史）的顺序,因以后在打印或绘图时要用。

历史值在FLAC停止运行时删除,因此如果要保留的话,要用到命令HISWRITEN（见下面）。

可以用的关键词有:

Ang单元i,j内的最小主应力和X-轴所形成的夹角

PP单元i,j内的孔隙压力

SIG1单元i,j的最大主应力

SIG2单元i,j的最小主应力

SXX单元i,j的xx-应力

SYY单元i,j的yy-应力

SXY单元i,j的xy-应力

X结点i,j的x-坐标

Y结点i,j的y-坐标

XDis结点i,j的x-位移

YDis结点i,j的y-位移

Unbal最大不平衡力

XVel结点i,j的x-速度

YVel结点i,j的y-速度

XXA单元i,j三角形a的xx-应力

XYA单元i,j三角形a的xy-应力

YYA单元i,j三角形a的yy-应力

XXB单元i,j三角形b的xx-应力

XYB单元i,j三角形b的xy-应力

YYB单元i,j三角形b的yy-应力

XXC单元i,j三角形c的xx-应力

XYC单元i,j三角形c的xy-应力

YYC单元i,j三角形c的yy-应力

XXD单元i,j三角形d的xx-应力

XYD单元i,j三角形d的xy-应力

YYD单元i,j三角形d的yy-应力

His关键词

关键词可以是:

Dumpnhis将第nhis的历史写屏;

Writenhis将第nhis的历史写在文件FLAC.HIS上,该文件可在FLAC结束后打印出来。

连续执行HISWRITE命令可将结果顺序写在FLAC.HIS上,但是首次执行此命令会将以前盘上所存的同名文件冲掉;

Reset所有的历史都清除掉。

Initial关键词=值<...><范围>

某些结点值可以给初值:

PP空隙压力

Xx-坐标

Yy-坐标

SXXxx-应力

SYYyy-应力

SXYxy-应力

XDispx-位移

YDispy-位移

XVelX-速度

YVely-速度

还有几个可选用的关键词以协助实现INITIAL命令:

Mark只有标记的结点值将被初始化

Regionij

为标记结点所围住的区域将被初始化。

i,j为标记区域中的任一个单元

Varxv,yv在一定范围内的参数的变差值,xv和yv为x-方向和y-方向的变差值。

范围可以有形式I=i1,j1J=j1,j2,其中i和j的先后次序无关。

INTerfacen关键词

n关键词=值

INTERFACE为交界面命令。

网格的一部分可以和网格的另一部分通过交界面而相互作用。

界面的性质由其刚度,粘结力和摩擦力来表征。

用于本命令的关键词有:

Asidei1,j1

Bdidei1,j1

Cohesion粘结力值

Friction摩擦力值

Glued

KN法向刚度值

KS切向刚度值

TBond抗拉强度

Unglued

胶结命令GLUE将使交界面的上下两部分胶结在一起,使之不能滑动和分离,命令UNGLUE则为命令GLUE之逆。

MArk<范围>

在给定范围内的结点将被标记。

计算的进行与结点的标记与否无关,但标记的结点可以限定一个区域作为命令INITIAL,PROP和MODEL的作用域。

命令GEN将自动对结点进行标记。

<范围>的形式可以是I=i1,i2,J=j1,j2,i和j的顺序无关。

Model关键词<范围>

该命令对于给定的区域或范围赋予有关的本构模型。

关键词有: