fla3d命令集锦

;****************************************group命令group soil range z 1 2 ;定义groupprop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值model null range group soilmodel elastic range group soilplot block group range group dam ;只显示damrange命令range name trench x 0 1 y 0 4 z 0 2model null range trenchmodel null range x=2,4 y=2,6 z=5,10range name Big_Brick x -3 3 y -2 2 z -1 1model elastic range Big_Brickprop bulk 1e8 shear 1e8 range Big_Brickrange name Layer1 plane dip 0 dd 0 ori 0 0 0 aboverange name Layer2 plane dip 0 dd 0 ori 0 0 0 belowrange cylinder end1 x1 y1 z1 end2 x2 y2 z2 radius rcylindrical 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 Layer1prop clay range layer2macro 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_Brickmacro 'Pt0' 'p0 15 15 15'gen Big_Brick ; this will cause an errorgroup range macro 命令Using different object types to do the same jobusing a RANGE object ...range name Big_Brick x = (-3,3) y = (-2,2) z = (-1,1)model null range Big_Brickusing a GROUP object ...group Big_Brick range x = (-3,3) y = (-2,2) z = (-1,1)model null range group Big_Brickusing a macro object ...macro Big_Brick 'x = (-3,3) y = (-2,2) z = (-1,1)'model null range Big_Brickgen zone brick size 10,10,10macro 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 10prop SiltySand range z 5 10model mohr range z 0 5prop ClayeyGravel range z 0 5;****************************************视图操作：x y z mctr+r 还原ctr+p 保存图片ctr+z 鼠标选择ctr+g 彩图变为灰色图shift + （x y z m）旋转缩放edit---copy to clipboard----粘贴到word;****************************************ini xdis 0 ydis 0 zdis 0set grav 0 0 -9.81set 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.txthist gp zdisp 4,4,8hist 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.logset log onprint …………set log off;****************************************切片pl set plane ori 0 1.5 0 norm 0 1 0 '定义剖面位置pl con zd plane '显示变量pl add skepl add dis planepl add axepl 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 ………… th en…………else…………endifcommand…………endcommandloop n（1,21）…………endlooploop while ……………………endloopcaseof ……………………case n1…………case n2…………endcase;****************************************运算= # > < >= <=负数在运算时，要加（），以免和减号- 混淆degrad 'π/180pi 'πngp '节点总数nzone '单元体总数;****************************************单元/节点遍历p_z = zone_headloop while p_z # null…………p_z = z_next(p_z)endloopp_gp = gp_headloop 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_namecal 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 lineplot table 1 both '点、线同时显示常用编辑文本文件的方法进行表的读入与调用第1行：表的名称第2行：x1 y1第3行：x2 y2……空行' 注：在表的文本文件最后，需要有一个回车换行符，否则会出现"error reading file xxx.dat"的错误，完成的文件需要读入操作才可以供flac3d调用。