FLAC3D与FLAC常见命令

FLAC-命令流及其解释⼤全1. FLAC3D是有限元程序吗？答：不是！是有限差分法。

2. 最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。

3. 怎样看模型的样⼦？答：plo blo gro可以看到不同的group的颜⾊分布4. 怎样看模型的边界情况？答：plo gpfix red sk5. 怎样看模型的体⼒分布？答：plo fap red sk6. 怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis)应⼒：plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的⽮量图？答：plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点？答：plo info9. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe⽂件10. 怎样调⽤⼀个⽂件？答：File/call或者call命令10. 如何施加⾯⼒？答：app nstress11. 如何调整视图的⼤⼩、⾓度？答：综合使⽤x, y, z, m, Shift键，配合使⽤Ctrl+R，Ctrl+Z等快捷键12. 如何进⾏边界约束？答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）13. 如何知道每个单元的ID？答：⽤⿏标双击单元的表⾯，可以知道单元的ID和坐标14. 如何进⾏切⽚？答：plo set plane ori (点坐标) norm (法向⽮量)plo con sz plane (显⽰z⽅向应⼒的切⽚)15. 如何保存计算结果？答：save +⽂件名16. 如何调⽤已保存的结果？答：rest +⽂件名；或者File / Restor17. 如何暂停计算？答：Esc18. 如何在程序中进⾏暂停，并可恢复计算？答：在命令中加⼊pause命令，⽤continue进⾏继续19. 如何跳过某个计算步？答：在计算中按空格键跳过本次计算，⾃动进⼊下⼀步20. Fish是什么东西？答：是FLAC3D的内置语⾔，可以⽤来进⾏参数化模型、完成命令本⾝不能进⾏的功能21. Fish是否⼀定要学？答：可以不⽤，需要的时候查Mannual获得需要的变量就可以了22. FLAC3D允许的命令⽂件格式有哪些？答：⽆所谓，只要是⽂本⽂件，什么后缀都可以23. 如何调⽤⼀些可选模块？答：config dyn (fluid, creep, cppudm)后注：这个⼯作很繁琐，需要的时间很多，希望⼴⼤⽹友能够将⾃⼰曾经遇到的常见问题在后续跟贴，也为了将这个FAQ进⾏很好的充实。

建模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlyt hefilenameneedbeenteredwiththe CALL command.②不在同一个文件夹，全路径：Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfol der\file.dat).Undo;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negative z-direction.（垂直向下！常用的）下下面面这这代代码码，，是是沿沿着着--y y方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续p110,-10,-20...p2-10,10,-20...p3-10,-10,0plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体p110,-10,-20p2-10,10,-20...p3-10,-10,0p410,10,-20...p5-10,10,10p610,-10,0...p710,10,10plotcurrentplotPlot01plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01y01z01;显示指定范围内各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：eachhistoryisnumberedsequentiallyfrom1asitisenteredviathe HISTORY co mmand.②查找显示所有的his命令：ReturntotheFlac3D>promptandtypelisthist foralistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。

FLAC3D 常用命令流;模型镜像gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面gen zone reflect normal 0 0 -1 ;沿z轴镜像;绘图控制pl <add> contour szz outline on;在模型中显示位移-应变曲线hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror";对模型进行压缩实验的方法;即在模型两侧施加相反方向的速度ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1;修改模型的坐标值ini x add -100 y add -100 z add -100;显示云图的同时也显示模型网格轮廓plot add cont disp outline on <gradient/average>;gradient更精确;输入角度、弧度方法pi=π，90°为90.0*degraddef set_valsptA = 25.0 * sin(pi/2);ptA=25.0ptB = 25.0 * cos( 60.0*degrad );ptB=12.5ptC=pi;ptC=3.1415926endset_valsprint ptA ptB ptC;施加结构单元方法sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &end2=(0.0,25.0,0.0) radius=24.5 notplot add sel geom black black cid on scale=0.03sel node init zpos add -25.0;如何显示某一平面plot create name_planeplot set plane origin 3 4 0 normal 1 0 0plot add cont disp plane behind shade onplot add sel geom blackplot add axes redplot current name_planeplot show;如何施加法向应力并显示所施加的力apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力plot grid green fap black;显示施加的力;获取某网格点的ID，坐标值def aaget_id=gp_id(gp_near(3,0,0))endaaprint get_idprint gp position range id 3;保存命令流文件set log onset logfile hua.log;文件存放的路径与调用的文本路径相同......set log off;fish函数的说明与解释def find_add ;定义fish函数find_addhead = null ;给head赋值p_gp = gp_head ;第一个网格结点的指针赋给p_gploop_while p_gp # null ;当p_gp值不为null时作循环x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_posif x_pos = 10.0 then ;如果x_pos = 10.0 则(执行)new = get_mem(2) ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址mem(new) = head ;将head类型和数值置于地址为new的fish变量mem(new+1)=p_gp ;第二个变量存储符合条件的zone地址head = new ;将new值赋给headendifp_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gp endloop;结束循环end ;结束fish函数;实际上这个fish函数为满足条件（x坐标为10的）的zone的地址开辟一定的地址空间，各地址之间存在一定的联系;找到第一个符合条件的zone地址后，用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址，并用;下面的mem(new) = head ，将第一次开辟的两个变量的第一个变量存储地址head（注意第一次head＝null），第二个变量存储第一;个符合条件的zone地址，并将第一个变量的地址赋予head（head＝new），第一次循环结束；下次循环，同样开辟两个变量对象，第;一个变量对象记录上次循环开辟的第一个变量的地址，第二个记录第二个符合条件的zone地址，其余循环依次类推，这样子就建立了;一个符合条件的zone地址链条，方便以后使用。

flac常用命令（FLAC common commands）FLAC3D common commands1. What are the first commands you need to master? Answer: commands such as gen, ini, app, PLO, solve, such as modeling, initial conditions, boundary conditions, postprocessing, and solving are required.2. How do I output the post-processing diagram of the model? Answer: File/Print, type/Jpg, file, and then select File/Print, save the format as JPE file.3. How do I call a file? Answer: File/call or call command 4, how to exert surface force? Answer: app nstress 5, how to adjust the size and angle of view? Answer: use x, y, Z, m, Shift keys, and use Ctrl+R, Ctrl+Z and other shortcut keys together. 6, how to carry out boundary constraints? Answer: fix x ran (the constraint is speed, in the initial case, the constraint is equivalent to the displacement constraint). 7, how do you know the ID of each unit? Answer: double-click the surface of the cell with the mouse to know the ID and coordinates of the unit. 8, how to slice? Answer: PLO, set, plane, ori (point coordinates), norm (normal vector), PLO, con, SZ, plane (showing the slice of the stress in the direction of the z) 9, how to save the result of calculation? Answer: save + file name 10, how to call saved results? Answer: rest + file name; or File / Restor 11, how to pause the calculation? Answer: Esc 12, how to pause in the program, and can restore the calculation? Answer: add the pause command to the command and proceed with continue. In our step by step solution, we want to get the result of a certain process, not wait until all the requests are finished, and then we can correct the error when the distribution is wrong, not wait for the result to come out. 13, how do you skip a calculation step? Answer: in the calculation, press the space bar to skip thiscalculation, automatically go to the next 14, Fish what is it? Does Fish have to learn? Answer: it's the built-in language of FLAC3D that can be used to parameterize models and perform commands that cannot be performed by themselves. Fish doesn't have to learn. When you need it, check Mannual and get the variables you need. What is the command file format allowed by 15 and FLAC3D? Answer: does not matter, as long as it is a text file, what suffix can be. 16. How do I call some optional modules? Answer: config dyn (fluid, creep, cppudm) 17, how to impose constraints around the cylinder? You can use fix, ran, cylinder, end1, end2, radius, R1, cylinder, end1, end2, radius, R2, not, R219, about DD, dip dip-direction angle, measured in: the global xy-plane clockwise from the positive y-axis from y to X is the shaft clockwise turn angle; dip dip angle, measured in the negative Z-direction from the global xy-plane this is the angle between the XY plane and Z axis; dip the angle between the reference plane and the XY plane symmetry, symmetry plane and XY angle reference plane is a ray projection in the XY plane; DD refers to the positive direction of axis Y clockwise direction steering angle (clockwise) projection to the normal direction is determined in the XY plane, and then determine the dip. Refers to dip XY from the plane to the plane defined by the angle. 20, range command: 1 spherical range, is a center, part of the different radius between the force applied to the sphere, with the constraint, command, command: range annulus center XYZ radius R1 R2 (Center) (two radius), want to express is spherical, r1r2 value difference can be small, such as the radius 2.9 3.1 when the range is a radius of 3 plane. The order of force is applied to the torus of a cylinder, exerting forceor restraint.范围缸end1 XYZ end2 XYZ半径R（端是表示柱体断面的中心点坐标）范围平面以上表示的是Z大于0的上部结构范围面吹表示的是Z小于0的下部结构范围DD DIP B XYZ表示的是走向角度起源，倾向角度B 的面DD表示，当面与xoy面相交时的垂线与Y轴正向的夹角，顺为正，说白话，就是面的指向方向面是朝那个方向的意思。

FLAC3D常用命令总结第一篇：FLAC3D常用命令总结1、模型建好如何检验所建模型是否有问题注意：第一步建模一定不要建太大的模型，可以建非常小的模型，确保模型的准确性，这样可以尽快开始后边的步骤，等后边计算分析步骤快速做完验证没有大问题，再将准确的模型建好换上，重新计算得到更加准确的结果。

2、如何检验应力边界、位移边界、加载速率等是否准确Plot add gpfix：看位移边界Plo add fap：看应力边界，看加载方向是否准确 Plot add vel：看速率，加载速率是否准确3、如何查看运算结果Plo con szz outline on：看z向应力云图，并显示边界网格 Plo add disp：看位移Plo con zdisp outline on magf 5：看z向位移云图，变形扩大5倍 Plo add bo ：添加边界 Plo blo state：看是否屈服Plo add stensor：查看大主应力方向Plo con smin：小主应力Plo con smax：大主应力Plo con ssi（shear strain increment）：查看剪切应变增量云图Plo bcont smin：查看单元小主应力4、如何做切面并查看切面结果在编辑“Edit“菜单中设置切面位置，在上述3命令后加上plane （pl）即可。

如：Plo con smin pl：查看切面小主应力5、后处理将模型导入到tecplot中做切面，加等值线等，对结果进一步处理，得到更漂亮的一些结果。

其他重要命令总结：1、如何在圆柱体四周如何施加约束条件？可以用fix...ran cylinder end1 end2 radius r1 cylinder end1end2 radius r2 not，其中r22、如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

1. apply（缩写：app）可用来定义边界条件及初始条件：1）添加应力格式1：apply szz -0.3395e6 range z -0.1 0.1格式1：apply szz -0.3395e6 range group pile格式3：apply szz -0.3395e6 range z -0.1 0.1 group pile格式4：apply nstress 数值range z 2.9 3.1（或3）x 1 2 y 1 22）以一定速度施加位移边界格式1：apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step3）添加边界条件格式：apply szz <数值> grad <梯度> range <范围>示例1：apply szz -1e9 grad 0 0 8.3e5 range z 0 120示例2：apply szz -0.6e6 range z 0.05 0.15 group pile注：<数值>是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。

2. range(缩写：ran)通过range功能，可以使命令作用在一定指定范围的目标上；如果一个命令没有使用range来确定范围，则命令对整个模型有效。

1)利用坐标指定一定的范围格式1：range z 0 1格式2：range z 2.9 3.1 x 1 2 y 1 22)利用分组来指定范围格式：range group 13)以上两种的复合格式：range z -0.1 0.1 group pile4）利用id号来指定一定的范围格式：range id 0 10该命令后跟起始id和结束id，这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。

例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

（完整版）FLAC3D命令流（整理版）1、怎样查看模型？答：plot grid 可以查看⽹格，plot grid num 可以查看节点号。

2、请问在圆柱体四周如何施加约束条件？答：可以⽤fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把⼀个PLOT的图像数据导出来以便⽤其他软件绘图？答：⽤set log on 命令，把数据导出来，转到excel⾥处理⼀下，然后⽤surfer或者什么作图软件绘制就⾏了。

4、⽤命令建⽴模型后，如何显⽰点的坐标？答：使⽤plo blo gro gpnum on命令5、关于gauss_dev对性质进⾏⾼斯正态分布的问题？答：根据⼿册上的说明：下⾯的命令设定⼀个平均摩擦⾓为40度，标准⽅差是±5％。

则命令如下：prop friction 40 gauss_dev 2问题：请问gauss_dev 2中的2是如何计算的？如果把±5％改为±10％，则命令应如何写？40×5％＝26、reflect问题问：gen zone radbrick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fillplot surfgen zone reflect dip 0 dd 90 （对xy⾯做镜像）gen zone reflect dip 90 dd 90 （对yz⾯做镜像）（1）dd表⽰y轴正向顺时针到那条射线的夹⾓，dip表⽰对称参照⾯与xy平⾯的夹⾓，对称参照⾯与xy平⾯的夹⾓在xy平⾯的投影是⼀条射线。

建 模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：If the called file is located in the same folder as the FLAC3D project file, then only the file name need be entered with the CALL command.②不在同一个文件夹，全路径：Otherwise, the file may be called by specifying its complete path (e.g., c:\myfolder\file.dat ).Undo ;撤销上一条命令2、创建 旋转 缩放视图3、建模命令model mech mohr ;莫尔库伦模型model mech elastic ;弹性模型set grav 0, 0, -9.81 ; 重力加速度 negative z -direction. （垂直向下！常用的）下下面面这这代代码码，，是是沿沿着着--y y 方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续p1 10, -10, -20 ...p2 -10, 10, -20 ...p3 -10, -10, 0plot zonegen zone brick size 6,8,8 p0 -10, -10, -20 ... ;不规则六面体p1 10, -10, -20 p2 -10, 10, -20 ...p3 -10, -10, 0 p4 10, 10, -20 ...p5 -10, 10, 10 p6 10, -10, 0 ...p7 10, 10, 10plot current plot Plot01plot clearplot zoneUndo ;撤销命令set logfile 127X1001.txtset log on truncateset log offlist zone prin range x 0 1 y 0 1 z 0 1 ;显示指定围各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：each history is numbered sequentially from 1 as it is entered via the HISTORY command.②查找显示所有的his命令：Return to the Flac3D> prompt and typelist histfor a listing of the histories and their corresponding numbers.hist nstep 5 ;每5步记录1次。

建模1、调用文件：①文件与工程在同一个文件夹，只写文件名即可：Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlPthefilenamen eedbeenteredwiththe CALL command. ②不在同一个文件夹，全路径：Otherwise,thefilemaPbecalledbPspecifPingits completepath(e.g.,c:\mPfolder\file.dat ). Undo ;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr ;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negative z -direction.（垂直向下！常用的）下下面面这这代代码码，，是是沿沿着着--P P 方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续 p110,-10,-20... p2-10,10,-20... p3-10,-10,0 plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体 p110,-10,-20p2-10,10,-20... p3-10,-10,0p410,10,-20... p5-10,10,10p610,-10,0... p710,10,10plotcurrentplotPlot01 plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01P01z01;显示指定范围内各单元的主应力，结果如下Hist命令：①命令编号按顺序从1开始：eachhistorPisnumberedsequentiallP from1asitisenteredviathe HISTORP command.②查找显示所有的his命令：ReturntotheFlac3D>promptandtPpelisthistforalistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。

flac3d常用命令师姐整理。

pl sel rec sres my surfx 0 1 0 ra y 0 24 scale=0.05pl add fap ra cyl end1 0 0 -17 end2 0 42 -17 r 3.1pl add sel geom id on;查局部坐标pr sel node local ra id=1pr sel rec surfpl sel geom node off nodesys on scale=0.06pl sel geom node off surfsys on scale=0.06pl sel pile fo fx; 看桩轴力pl bl gr id on ra y 20 21;查看zoneIDpl sur gp on ra y 20 24 ;监测点位移，找点gpidpl sel geom cid on ra y 22 24 ;监测衬砌内力，找衬砌cidpl sel geom cid on node off ra y 21 23 ;监测桩号，找桩cidpl bl state shear-n ;看屈服pl add gpfix ;查固定边界条件pl add su magf 10 ;变形后位移pr gp disp ra id=12067;查看某点位移pl bcon prop bu;无刚度错误错误pr sel node link ra id ;找node连接的link号pr sel link tar ra id ;找link的目标节点pr sel link sou ra id ;找link的源节点pr sel link attachpr sel link c lind k;查link的刚度pr sel rec sr my;pl sel disp magf 100;pl sel pile fo fx;pl sel pile mo my;pl con zdis。