FLAC3D命令流整理版.docx

1、FLAC3D命令流整理版1、怎样查看模型？答：plotgrid可以查看网格，plotgridnum可以查看节点号。2、请问在圆柱体四周如何施加约束条件？答：可以用fix.rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r2r1，其实就是选择range的时候选两个圆柱的差，即得到边界。命令流如下：fixxrangeend1100end2140rad1end1100end2140rad13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或

2、者什么作图软件绘制就行了。4、用命令建立模型后，如何显示点的坐标？答：使用ploblogrogpnumon命令5、关于gauss_dev对性质进行高斯正态分布的问题？答：根据手册上的说明：下面的命令设定一个平均摩擦角为40度，标准方差是5。则命令如下：propfriction40gauss_dev2问题：请问gauss_dev2中的2是如何计算的？如果把5改为10，则命令应如何写？40526、reflect问题问：genzoneradbrick&p0(0,0,0)p1(10,0,0)p2(0,10,0)p3(0,0,10)&size3,5,5,7&ratio1,1,1,1.5&dim142fi

3、llplotsurfgenzonereflectdip0dd90（对xy面做镜像）genzonereflectdip90dd90（对yz面做镜像）（1）dd表示y轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。首先应该按照dd的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角（及z轴负方向转过角度）（2）命令：reflectnormalxyzoriginxyz（根据法线和过一点

4、建的面）最常用，至于实际使用过程中有人问道：normal-100与normal100的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。reflectddadipboriginxyz这个就是用里建立斜面镜像的。有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字，genzonereflect7、关于两种fix是否有区别？fixzrangez-0.10.1fixxrangex-0.10.1anyx1.92.1anyfixyrangey-0.10.1anyy1.92.1any和fixzrangez-0.10.1fixx

5、rangex-0.10.1fixxrange1.92.1fixyrangey-0.10.1fixyrangey1.92.1答：是一样的，any表示并列连接关系，及and的意思。8、initialapplyset的区别答：initial初始化命令，初始化计算应力状态，定义材料参数、重力等；apply边界条件限制命令，如施加边界的力、位移等约束等；set用于设置计算的模式，控制计算过程。Ini施加的初始力在计算过程中是会变的，一般体力需要初始化，而apply施加的力、边界条件不会随着计算变化。set命令作用于整个网格(赋值全局变量),ini则初始化指定范围内网格点、域的偏移量。初始应力不仅仅是重力

6、引起的。9、Tension-pshear-ptension-nshear-n表示的含义答：在计算的循环里面,每个循环中,每个zone都依据failurecriterion处于不同的状态,shear和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态。10、history保存数据，print显示在错误的命令流中“原点位移.dat”文件中只有：“Step1Y-Displace”而没有保存下数值，怎么回事？我怎么才能获得table中的数据？输入printtabeln能得到table中的数据（正确）histgpydisp

7、0,0,0histzonesyy0,1,0histzonesyy1,1,0step3000histwrite1file原点位移.dathistwrite2file中部单元应力.dathistwrite3file边缘单元应力.dathistwrite1table1histwrite2table2histwrite3table3histwrite23vs1file单压.dat（错误）histgpydisp0,0,0histzonesyy0,1,0histzonesyy1,1,0histwrite1file原点位移.dathistwrite2file中部单元应力.dathistwrite3file边

8、缘单元应力.dathistwrite1table1histwrite2table2histwrite3table3histwrite23vs1file单压.datstep300011、plot显示x、y坐标图histn1histgpydisp0,0,0histax_strhistgpxdisp1,1,0step1500plothist2vs-1;axialstressvsaxialdisp.plothist3vs-1;circumferentialdisp.vsaxialdisp答：那2是不是代表histgpydisp0,0,0，而-1指得是histgpydisp0,0,0为什么不是1而是负值

9、？1)在plothistmvsn的形式里,m代表y轴,n代表x轴(不管m,n的正负);2)-表示对其值作mirror,如果大家使用AUTOCAD的话就明白mirror的含义了,比如上面的hist1gpydisp的值是0100,那么vs-1就变成-1000.以此类推.3)要这么做主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内。12、关于water.fis的命令注释（设置地下水位）答：（1）water.fis的功能是设置地下水位面，在二维计算中，通过watertable命令即可设置水位线，在三维计算中由于地下水位是面，因此需要作一些插值运算，以便形成水位面。这

10、个FISH程序包含三个函数parm,aux1和aux2。parm用来设定水位的边值，即设置水位的区域，aux1插值计算等值线的点，aux2调用aux1的计算结果，使用FLAC3D命令形成水位面。（2）首先设置参数值，然后定义一个table.defparmnptab1=7nprof=10endparmtable1(120)(18.2512.5)(2517)(3521)(5024)(7525)(15025)产生table1，包含7组数据点。3)在解释的两个FISH函数之前，对table,xtable,ytable作一些总的说明。water.fis程序反复使用了这三个变量，理解了这三个变量的含义，w

11、ater.fis这段程序就解决了大半问题了。Itasca的所有软件都使用了相同的数据结构储存数据，称之为链表结构，类似于关系型数据库中的数据链接方法。如果有数据结构的知识，或掌握关系型数据库的SQL语言，这三个变量的含义就不难理解了，只不过这个操作与真正的数据库操作相比，只有两个数据变量x,y。我用数据库的术语作解释：table就是产生一个“表”，xtable,ytable的作用就是产生一个新表(Create)以及对表中的数据进行取出(Select)、插入(Insert)、替换(Update)操作.4）FISH函数aux1的作用是在table1的基础上产生出另外个table，table的ID从

12、11到17，每个table有10个数据集。alfa的计算方法是为适合(fit)曲面而采用的坐标点算法。defaux1loopn(1,nptab1)rr=xtable(1,n)zz=ytable(1,n)loopk(1,nprof)case_ofkalfa=0.5*(pi+.1)*float(k-2)/float(nprof-2)xtable(n+10,k)=rr*cos(alfa)ytable(n+10,k)=rr*sin(alfa)case1xtable(n+10,k)=rrytable(n+10,k)=-40.0case2xtable(n+10,k)=rrytable(n+10,k)=0.

13、0end_caseend_loopend_loopendaux1当执行这个FISH函数后，生成(nptab1+10,nprof)组数据，这些数据为Aux2调用。注意：本段程序产生的数据是为那个特定的曲面边坡制作的，不是一个generalpurpose的LIB5）FISH函数aux2的作用是使用WATERTABLEface命令产生水面，把table1的ytable值作为z坐标，aux1产生的数据作为x,y坐标。defaux2commandsetgrav00-10waterdens1000end_commandloopn(1,nptab1-1)loopk(1,nprof-1)xx1=xtable(

14、n+10,k)yy1=ytable(n+10,k)zz1=ytable(1,n)xx2=xtable(n+10,k+1)yy2=ytable(n+10,k+1)zz2=zz1xx3=xtable(n+11,k+1)yy3=ytable(n+11,k+1)zz3=ytable(1,n+1)xx4=xtable(n+11,k)yy4=ytable(n+11,k)zz4=zz3commandWATERTABLEfacexx1,yy1,zz1xx2,yy2,zz2xx3,yy3,zz3WATERTABLEfacexx1,yy1,zz1xx3,yy3,zz3xx4,yy4,zz4end_commande

15、nd_loopend_loopendaux213、FLAC结果数据的导出plsetbackwhite;设置背景setplotbitmapsize(800,510);设置图片大小plotsetcaptionsize37;设置图例数字大小(default=35范围10-50)mainwinsize11position00.8;commandwindowplotsetwindowsize0.80.8position0.10.0;graphicswindowplotsetcenter100040根据情况更改14、一段fish命令的解释defaux1定义函数loopn(1,nptab1)开始一个循环，从

16、1循环到变量nptab1rr=xtable(1,n)从table1的x列中读第n个数付给变量rrzz=ytable(1,n)从table1的y列中读第n个数付给变量zzloopk(1,nprof)开始第二个循环，从1循环到变量nprofcase_ofk下面是一个选择结构，利用变量k判断alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2)xtable(n+10,k)=rr*cos(alfa)ytable(n+10,k)=rr*sin(alfa)case1当k1的时候执行下面语句xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rrytabl

17、e(n+10,k)=-40.0tablen+10的y列第k个值等于变量-40case2当k2的时候执行下面语句xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rrytable(n+10,k)=0.0tablen+10的y列第k个值等于变量0end_case结束选择结构end_loop结束第二个循环end_loop结束第一个循环end函数定义结束aux1执行什么定义的函数15、有关命令的解释（1）applynvel0.1planedip60dd270rangeleft_bounapply施加边界条件，plane法向速度为0.1，plane面的位置：倾角60（在xoy平面

18、从-Z轴方向逆时针旋转），dd走向270度，从+y方向绕z轴顺时针旋转；（2）iniszz-5e4grad00-1e4ini初始边界条件，szz=-5e4+x*0+y*0+-z*-1e4,注：-5e4是range范围起点坐标的SZZ值，x,y,z所施加应力网格点的坐标。16、渗流示例;-参数部分bulkmodulus,K390MPa体积模量（土质比较硬）shearmodulus,G280MPa剪切模量soildrydensity,d1200kg/m3土的密度waterdensity,w1000kg/m3水的密度walldensity,wal1500kg/m3档土墙的密度permeability

19、,k10?12m2/Pa-s渗透系数porosity,n0.3孔隙率fluidbulkmodulus,Kf2.0GPa流体的体积模量;-以下是命令流;-;Excavationinasaturatedsoil;-configfluid;设置流体算法;-geometricalmodel-建立模型genzonebrickp11200p20120p30012size121212rat111groupsoil;设置土组groupexcavaterangex04y04z05;设置开挖范围组groupwal1rangex45y05z07;设置档土墙1组groupwal2rangex04y45z07;设置档土

20、墙2组groupwallrangegroupwal1anygroupwal2any;档土墙1和2合并为一个档土墙组;-fluidflowmodel-流体模型modelfl_iso;各项同性流体propperm1e-12poro0.3;设置渗透系数和？比or？率inifdensity1e3;初始化水密度inifmod2.0e9ftens-1e-3;初始化水的体积模量及？modelfl_nullrangegroupexcavate;开挖部分没有水流modelfl_nullrangegroupwall;档土墙不渗水inipp0grad001e4;初始化孔隙水压力fixpprangez-0.10.1;

21、约束应该是整个模型的底部吧？fixpprangex-0.14.1y-0.14.1z4.95.1;开挖部分也给约束住了，而其他的地方不管;-mechanicalmodel-力学模型modelelas;弹塑性propbul3.9e6shea2.8e6;设置体积模量和剪切模量，这两个必须，有一次我只设置了体积模量，就错的一塌糊涂modelnullrangegroupexcavate;力学模型的开挖inidensity1.2e3;初始化密度，应该是土的密度inidensity1.5e3rangegroupwall;初始化墙（档土墙）的密度，fixxrangex-.1.1;施加约束，据说差分法约束0的时

22、候，必须取一个-0.100.1,范围fixxrangex11.912.1;施加约束，或者叫做边界条件fixyrangey-.1.1;施加约束fixyrangey11.912.1fixzrangez11.912.1;initialtotalstresses出世总应力iniszz0grad00-1.5e4;z方向的应力大一些inisxx0grad00-1.2e4;水平方向的，小一些inisyy0grad00-1.2e4;applynstress0grad00-1.2e4rangex0.04.0y3.94.1z0.05.0applynstress0grad00-1.2e4rangex3.94.1y0

23、.04.0z0.05.0applynstress-7.5e4rangex0.04.0y0.04.0z4.95.1;-setting-设置setgravity0010;设置重力加速度;-initialstate-;初始化状态solveforce1;checkinitialequilibrium求解平衡;-histories-记录sethist_rep40;记录步为40histfltime;记录时间histgppp006;记录某点的孔隙水压力histgpxdis400;水平位移histgpxdis402;记录水平位移histgpxdis420histgpzdis005;记录竖向位移histgpzd

24、is205histgpzdis405histgpzdis225histgpzdis425histgpzdis445histgpzdis1001histgpzdis1002;-excavation-开挖setfluidoff;关闭水;applyporepressureatwalls在墙上施加孔隙水压力applynstress0grad00-1.e4rangex0.04.0y3.94.1z0.05.0applynstress0grad00-1.e4rangex3.94.1y0.04.0z0.05.0applynstress-5.e4rangex0.04.0y0.04.0z4.95.1solve;f

25、orce1;求解saveexc1.sav;保存;-drainage-排水applyremovenstress;撤掉刚才的压力defrelaxsetup;定义一个释放函数的参数step0=stependrelaxsetup;调用defrelax;定义一个释放函数rstep=step-step0ifrstepgp0_dispgp0_disp=gp_dispendifp_gp=gp_next(p_gp)endloop;找最大值的坐标p_gp=gp_headss=0loopwhilep_gp#nullgp_disp=gp_xdisp(p_gp)*gp_xdisp(p_gp)gp_disp=gp_dis

26、p+gp_ydisp(p_gp)*gp_ydisp(p_gp)gp_disp=gp_disp+gp_zdisp(p_gp)*gp_zdisp(p_gp)gp_disp=sqrt(gp_disp)ifgp_disp=gp0_dispss=ss+1xtable(1,ss)=gp_xpos(p_gp)ytable(1,ss)=gp_ypos(p_gp)endifp_gp=gp_next(p_gp)endloopendget_gp_maxdisp18、如何显示计算后模型的变形状态？答：plskmagfac*(一个放大的数字)19、如何查看剖面上任一点的位移值答：printgpdisprangex()y()z()或printgpdisprangeidplotgpdisprangex()y()z()或plotgpdisprangeid20、group命令划分斜面答：groupgroupnamecolourrangeplaneabovedd*dip*origin*planebelowdd*dip*origin*dd和dip就是两个斜面的倾向和倾角,origin是这个斜面通过的某个点,这个命令