搜索
第1页共13页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个时步的循环运算。第2页共13页FixXMark范围YXY用此命令可以使范围内结点的x-或y-方向的速度保持不变。范围的格式可以是I=i1,i2,J=j1,j2;i和j何者先输入没有关系。如果要求位移固定,则必须将速度初始化为零(开始时速度的缺省值为零)。FIX和INIXV,YV可以联合使用以提供一个刚性移动的边界条件。如果有MARK的关键词,则只有在此范围内被标记的结点将被固定。FReeXMark范围YXY该命令与FIX命令相反,用于放松对结点的约束。GEnx1,y1x2,y2x3,y3x4,y4RatiorirjrangeSameSameSameSame用于在全域或局域中产生网格。(x1,y1)的点放在左下侧,(x2,y2)、(x3,y3)和(x4,y4)为顺时针排列。如果点和前一个GEn命令中的相同则用SAME来代替。如果网格不是均匀排列,则用关键词RATIO,在I和J方向的比例系数分别为ri和rj.GEnCirclexc,ycradArcxc,ycxb,ybthetaLinex1,y1x2,y2该命令产生园、弧或直线。ARC弧的圆心为(xc,yc),起始点为(xb,yb),逆时针的角为theta度CIRCLE园的中心为(xc,yc),半径为radLINE直线的两端为(x1,y1)和(x2,y2).用ARC等关键词所形成的边界系被标记(参见命令MARK和UNMARK)过,标记结点所围住的区域在以后用命令INITAL,MODEL和PROPERTY时可以用关键词REGION来说明。GenADJust用此命令来微调网格,使之离散平滑。该命令可以连续使用以增进平滑效果的。Gridicoljrow此命令用以产生数目为icol列和jrow行的网格。对于640k内存的FLAC版本,如用摩尔-库仑的本构模型约可分两千个单元。第3页共13页Help帮助命令,可在屏幕上显示命令表。HisNstep=n关键词...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-应力第4页共13页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关键词i1,j1i2,j2n关键词=值INTERFACE为交界面命令。网格的一部分可以和网格的另一部分通过交界面而相互作用。界面的性质由其刚度,粘结力和摩擦力来表征。用于本命令的关键词有:Asidei1,j1i2,j2i3,j3...第5页共13页Bdidei1,j1i2,j2i3,j3...Cohesion粘结力值Friction摩擦力值GluedKN法向刚度值KS切向刚度值TBond抗拉强度Unglued胶结命令GLUE将使交界面的上下两部分胶结在一起,使之不能滑动和分离,命令UNGLUE则为命令GLUE之逆。MArk范围在给定范围内的结点将被标记。计算的进行与结点的标记与否无关,但标记的结点可以限定一个区域作为命令INITIAL,PROP和MODEL的作用域。命令GEN将自动对结点进行标记。范围的形式可以是I=i1,i2,J=j1,j2,i和j的顺序无关。Model关键词REGION=i,j范围该命令对于给定的区域或范围赋予有关的本构模型。关键词有:ANisotropic横观各向同性弹性模型Elastic弹性,各同性本构模型Mohr-Coulomb摩尔-库伦塑性模型Mull零模型,用于开挖掉的单元SS应变软化Ubiquitous彻体节理模型NEW该命令可以不用退出FLAC重新开始一个新问题。Plot关键词开关=v...关键词...用缺省值。可以在一行上连写几个关键词以在一个图上画出几个变量。在绘图前先要用MODEL命令赋予网格以一定的本构模型。关键词及其意义为:BEam画出结构单元的几何形状Boundary画出网格的内外边界Cable画出锚束的几何形状Disp位移矢量第6页共13页E_p塑性应变的轮廓线(只限于应变软化材料)Grid画出画出纪录在第nhis个历史上的变量值PP孔隙压力的等值线RF以矢量形式画出的固定结点的反力STAte画出单元中心当前的屈服状态STress主应力矢量SXXxx-应力的等值线(全应力)SYYyy-应力的等值线(全应力)SXYxy-应力的等值线Velocity以箭头画出的速度矢量WAter水位线XDispx-位移的等值线YDispy-位移的等值线XVelx-速度YVely-速度所谓开关其本身也是关键词,用于设置图形的的某些特征,有:color如红为Red,绿为GREen,黄为Yellow等。Interval=c将等值线的间距置为c.在PLOT命令前发扫描线命令SCLIN可以画出等值线的值。Max=v在画矢量时,置箭头的最大长度为v.Noh不画标题Zero不画零值线如果要放大或缩小图形,可以在PLOT命令前用WINDOW命令来设置窗口的大小。Print关键词关键词...Region=i,j范围和绘图命令PLOT一样,结点变量的值只有在给定的材料模型和性质后才可以打印出来。关键词有:Apply打印出所施加力或压力的大小和范围Fix打印出固定x或y的结点Limits对SOLVE命令的限制值第7页共13页MArk打印出标记结点MEm占用的内存STruct打印出结构单元上有关结点的力,力矩和位移Interface打印出交界面的数据,包括结点力和单位法向矢量主要的网格关键词Bulk体积模量COhesion粘结力Den质量密度Dilation剪胀角E_p全塑性应变(应变软化模型)Friction摩擦系数SHear剪切模量XX-坐标YY-坐标PP孔隙压力JFric节理摩擦力JCoh节理粘结力JAngle节理角度NUYxy-x泊松比NUZxz-x泊松比SIG1最大主应力SIG2最小主应力SXXXX-应力SYYYY-应力SXYXY-应力Theta最小主应力与x-轴的夹角XDisX-位移YDisY-位移XModX-模量YModY-模量XVelX-速度第8页共13页YVelY-速度ASXXXX-应力(三角形A)BSXXXX-应力(三角形B)CSXXXX-应力(三角形C)DSXXXX-应力(三角形D)ASYYYY-应力(三角形A)BSYYYY-应力(三角形B)CSYYYY-应力(三角形C)DSYYYY-应力(三角形D)ASXYXY-应力(三角形A)BSXYXY-应力(三角形B)CSXYXY-应力(三角形C)DSXYXY-应力(三角形D)State塑性状态0弹性1正在屈服中2曾经屈服,现为弹性3已超过单轴抗张力4屈服并超过单轴抗张力5已经超过抗张力6彻体节理正在屈服7彻体节理过去屈服过,现为弹性状态Tables打印所存贮的表格XReactionX-反力YReactionY-反力PROp关键词=值...范围本命令为MODEL命令赋予材料的性质。下面的关键词对PROPERTY命令加以补充:Regioni,j所有为标记结点所连续包围的区域将赋予该性质。i,j单元可为标记区域中的任一单元Varxvyv第9页共13页性质可以有变差,xv和yv分别为x-方向和y-方向的变差值范围可为I=i1,i2,J=j1,j2,i和j孰者在前无妨,但范围必需在输入行的最后。各种模型所需输入的性质为:弹性(1)剪切模量(2)体积模量(3)密度摩尔-库伦(1)剪切模量(2)体积模量(3)密度(4)摩擦角(5)粘结力(6)剪胀角(任选)横观各向同性(1)剪切模量(2)x-模量(3)y-模量(4)密度(5)NUYx(6)NUYz彻体节理(1)剪切模量(2)体积模量(3)密度(4)粘结力(整体材料)(5)摩擦力(整体材料)(6)节理粘结力(7)节理摩擦力(8)节理角第10页共13页(9)整体材料的剪胀角(任选)应变-硬化/软化(1)剪切模量(2)体积模量(3)密度(4)初始摩擦角(任选)(5)初始粘结力(任选)(6)摩擦角和塑性应变的关系表的表号(7)粘结力和塑性应变的关系表的表号(8)剪胀角和塑性应变的关系表的表号如果表号给0,则取关键词COHESION,DILATION或FRICTION所给出的值。Quit该命令同STOP,FLAC终止运行。REstore文件名把以前用SAVE命令所存的文件读入内存恢复现场。RETurn该命令应是输入