《FLAC代码查询》由会员分享，可在线阅读，更多相关《FLAC代码查询（15页珍藏版）》请在金锄头文库上搜索。
1、FLAC程序使用手册FLAC 输入命令 FLAC 的输入和一般的数值模拟的程序不一样, 它可以用交互的方式从键盘输入各个命令, 也可以写成命令文件, 类似于批处理, 由文件来驱动。FLAC 命令大小写一样。所有的命令可以附带若干个关键词和有关的数值。在下面的命令解释中, 只有大写的字母起作用, 小写的字母写不写、写多少个都没有崐关系。i,j,m 和 n 开始的变量要求整型数, 否则要求实型数。实型数的小数点可崐以忽略, 但是整型数不能带小数点。数值间可以用空格隔开, 空格的数目不限,也可以用下面的分隔符隔开:( ), / = 表示可选的参数, 输入时括号不用输入; . 表示可以有任意个参数。
2、由 * 号开始到行末为注释, FLAC 在执行时不理会。 下面的 FLAC 命令按字母排列。Apply 关键词 = 数值 可以有下面的关键词: Pressure 压力 XForce X-方向的力 YForce Y-方向的力ATtach 该命令可以将一条线上的结点和另一条线上的结点互相接合在一起, 用以形成复杂的网格形状。Call 文件名 写成的命令文件可以用 Call 命令来调用, 命令文件的最后一行必须是RETURN, 以返回到交互方式。命令文件中不能有 CALL 命令本身。Config 关键词 FLAC 用以解平面应变问题, 但经过配置命令也可以用于解平面应力问题或轴对称问题。 需要时应在
3、形成网格之前发。关键词有: P_STR 平面应力问题 AX 轴对称问题CYC n 该命令同 STEP, 为执行 n 个时步的循环运算。Fix X Y X Y 用此命令可以使 内结点的 x- 或 y- 方向的速度保持不变。 的格式可以是 I = i1,i2, J = j1,j2; i 和 j 何者先输入没有关系。如果要求位移固定, 则必须将速度初始化为零 (开始时速度的缺省值为零)。FIX 和 INI XV, YV 可以联合使用以提供一个刚性移动的边界条件。如果有 MARK 的关键词, 则只有在此范围内被标记的结点将被固定。FRee X Y X Y 该命令与 FIX 命令相反, 用于放松对结点的
4、约束。GEn x1,y1 x2,y2 x3,y3 x4,y4 用于在全域或局域中产生网格。(x1,y1)的点放在左下侧, (x2,y2)、 (x3,y3) 和 (x4,y4) 为顺时针排列。如果点和前一个 GEn 命令中的相同则用 SAME 来代替。如果网格不是均匀排列, 则用关键词 RATIO,在 I 和 J 方向的比例系数分别为 ri 和 rj.GEn Circle xc,yc rad Arc xc,yc xb,yb theta Line x1,y1 x2,y2 该命令产生园、弧或直线。 ARC 弧的圆心为 (xc,yc), 起始点为 (xb,yb), 逆时针的角为theta 度 CIRC
5、LE 园的中心为 (xc,yc), 半径为 rad LINE 直线的两端为 (x1,y1) 和 (x2,y2). 用 ARC 等关键词所形成的边界系被标记 (参见命令 MARK 和 UNMARK)过, 标记结点所围住的区域在以后用命令 INITAL, MODEL 和 PROPERTY时可以用关键词 REGION 来说明。Gen ADJust 用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。Grid icol jrow 此命令用以产生数目为 icol 列和 jrow 行的网格。对于 640k 内存的FLAC 版本, 如用摩尔-库仑的本构模型约可分两千个单元。Help 帮助
6、命令, 可在屏幕上显示命令表。His 每隔 NSTEP 时步, 记录一下关键词所示项目在 i1 和 ji 结点或单元的数值。 NSTEP 的缺省值为 10. 用户应记住所要求 His (历史)的顺序, 因以后在打印或绘图时要用。历史值在 FLAC 停止运行时删除, 因此如果要保留的话, 要用到命令 HIS WRITE N (见下面)。 可以用的关键词有: Ang 单元 i,j 内的最小主应力和 X-轴所形成的夹角 PP 单元 i,j 内的孔隙压力 SIG1 单元 i,j 的最大主应力 SIG2 单元 i,j 的最小主应力 SXX 单元 i,j 的 xx-应力 SYY 单元 i,j 的 yy-应
7、力 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
8、,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 关键词 关键词可以是: Dump nhis 将第 nhis 的历史写屏; Write nhis 将第 nhis 的历史写在文件 FLAC.HIS 上, 该文件可在 FLAC 结束后打印出来。连续执行 HIS WRITE 命令可将结果顺序写在 FLAC.HIS 上, 但是首次执行此命令会将以前盘上所存的同名文件
9、冲掉; Reset 所有的历史都清除掉。Initial 关键词 = 值 某些结点值可以给初值: PP 空隙压力 X x-坐标 Y y-坐标 SXX xx-应力 SYY yy-应力 SXY xy-应力 XDisp x-位移 YDisp y-位移 XVel X-速度 YVel y-速度 还有几个可选用的关键词以协助实现 INITIAL 命令: Mark 只有标记的结点值将被初始化 Region i j 为标记结点所围住的区域将被初始化。i,j 为标记区域中的任一个单元 Var xv, yv 在一定范围内的参数的变差值, xv 和 yv 为 x-方向和 y-方向的变差值。 范围可以有形式 I=i1,
10、j1 J=j1,j2, 其中 i 和 j 的先后次序无关。INTerface n 关键词 n 关键词 = 值 INTERFACE 为交界面命令。网格的一部分可以和网格的另一部分通过交界面而相互作用。界面的性质由其刚度, 粘结力和摩擦力来表征。用于本命令的关键词有: Aside i1,j1 Bdide i1,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 关键词 该命令对于给定的区域或范围赋予有关的本构模型。关键词有:
