flac2d使用手册.docVIP

FLAC程序使用手册 FLAC 的输入和一般的数值模拟的程序不一样, 它可以用交互的方式从键盘输入各个命令, 也可以写成命令文件, 类似于批处理, 由文件来驱动。 FLAC 命令大小写一样。所有的命令可以附带若干个关键词和有关的数值。在下面的命令解释中, 只有大写的字母起作用, 小写的字母写不写、写多少个都没有崐关系。i,j,m 和 n 开始的变量要求整型数, 否则要求实型数。实型数的小数点可崐以忽略, 但是整型数不能带小数点。数值间可以用空格隔开, 空格的数目不限,也可以用下面的分隔符隔开: ( ), / = 表示可选的参数, 输入时括号不用输入; ... 表示可以有任意个参数。 由 * 号开始到行末为注释, FLAC 在执行时不理会。 下面的 FLAC 命令按字母排列,每贴以第一个字母开头。Apply 关键词 = 数值 关键词 = 数值 ... 范围 可以有下面的关键词: Pressure 压力 XForce X-方向的力 YForce Y-方向的力 ATtach 该命令可以将一条线上的结点和另一条线上的结点互相接合在一起, 用以形成复杂的网格形状。 Call 文件名 写成的命令文件可以用 Call 命令来调用, 命令文件的最后一行必须是RETURN, 以返回到交互方式。命令文件中不能有 CALL 命令本身。 Config 关键词 FLAC 用以解平面应变问题, 但经过配置命令也可以用于解平面应力问题或轴对称问题。需要时应在形成网格之前发。关键词有: 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 命令相反, 用于放松对结点的约束。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 度 CIRCLE 园的中心为 (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 帮助命令, 可在屏幕上显示命令表。 His 关键词 . . . I=i1 J=j1 每隔 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-应力 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