flac命令流

1、FLAC3D常见命令：1.FLAC3D是有限元程序吗？答：不是！是有限差分法。2.最先需要掌握的命令有哪些？答：需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。3.怎样看模型的样子？答：ploblogro可以看到不同的group的颜色分布4.怎样看模型的边界情况？答：plogpfixred5.怎样看模型的体力分布？答：plofapred6.怎样看模型的云图？答：位移：plocondis(xdis,ydis,zdis)应力：ploconsz(sy,sx,sxy,syz,sxz)7.怎样看模型的矢量图？答：plodis(xdis,ydis,zdis)8.怎样看模型有多少单元、节点？答：priinfo9.怎样输出模型的后处理图？答：File/Printtype/Jpgfile，然后选择File/Print，将保存格式选择为jpe文件10.怎样调用一个文件？答：File/call或者call命令10.如何施加面力？答：appnstress11.如何调整视图的大小、角度？答：综合使用x,y,z,m,Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键12.如何进行边界约束？答：fixxran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID和坐标14.如何进行切片？答：plosetplaneori(点坐标)norm(法向矢量)ploconszplane(显示z方向应力的切片)15.如何保存计算结果？答：save+文件名.16.如何调用已保存的结果？答：rest+文件名；或者File/Restore17.如何暂停计算？答：Esc18.如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，用continue进行继续19.如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步20.Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21.Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了22.FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块？答：configdyn(fluid,creep,cppudm)24.如何在圆柱体四周如何施加约束条件？答：可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？答：用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。26.关于dd、dip的理解：答：首先你要确定你定义面的法线方向，dip是倾角的意思（xoy面和定义面的夹角），法线向上dip是锐角，法线向下dip时钝角。dd是y轴顺时针转到法线在xoy面上的投影线时，所转过的角度。27.如何调用FISH文件？答：新手调用FISH文件出错的原因：是：文件名为fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的）。修改方法是：将fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入fishcall.fis，类型也就正确了。28.如何理解Tension-pshear-ptension-nshear-n表示的含义？答：在计算的循环里面,每个循环中,每个zone都依据failurecriterion处于不同的状态,shear和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态.29.绘制相关曲线命令的一些解释在plothistmvsn的形式里,m代表y轴,n代表x轴(不管m,n的正负);-表示对其值作mirror,如果大家使用AUTOCAD的话就明白mirror的含义了30.如何查看剖面上任一点的位移值？答：printgpdisprangex()y()z()或printgpdisprangeidplotgpdisprangex()y()z()或plotgpdisprangeid31.a:fixxyzrangex-0.10.1b:fixxrangex-0.10.1请问a式与b式有何区别？答：a表示在x=0这个平面是固定的；b表示在x=0这个平面只x方向不能发生运动32.initial与apply有何区别？答：初始化命令，如初始化计算体的应力状态等；apply边界条件限制命令，如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。33.如何查看各个时段不平衡力的具体数值？答：采用his来记录计算，包括位移应力等命令hisunbalhisgp(zone)zdisrange(000)或者id=？导出数据命令hiswritenvsmbegin时步end时步+filefilename.his要导出不平衡力的具体数值hisunbalstep100000orsolvehiswrite1vsstepbegin1end1000file123.his使用上述命令就可以查看各个时步下的不平衡力的具体数值34.如何显示变形轮廓线的命令？ploskemagf10其中10为放大系数35、怎样开始入门？答：初步学习首先看说明书User'sgurid的Section2Gettingstarted和Section3ProblemslovingwithFLAC36、哪些命令需要掌握？答：gridgenpropinterfaceinifixsethiststepsolvepausecontsav等，在学习的过程中遇到这些命令和不会的命令应该查看说明书的CommandReference，里面讲的很清楚，特别需要注意的是一定要仔细看一下CommandReference开头的OptionalRangePhrases，这是讲范围的，因为很多命令的施加都需要加到模型的一部分，这里就告诉你怎样来表示这一部分。这一原则同样适合与其他5个ITASCA软件。还有一点就是六个软件（FLACFALC3DUDEC3DECPFCPFC3D）的range大不相37、FLAC和FLAC3D建模有什么区别？答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出来，然后用MODELLNULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。然后再用GENLINE，GENARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。而FLAC3D是直接用GEN命令一块一块搭接，GEN可以用很多次。38、怎样查看模型？答：plotgrid可以查看网格，plotgridnum可以查看节点号。39、怎样查看边界条件？答：plotfix40、怎样查看预先留设的开挖边界？答：plotgridmark。41、怎么查看模型的云图？答：plsyy(sxxxdisydissig1..)fill，注意，带fill关键字就是查看云图，不带fill关键字就是查看等直线图。42、怎样查看模型的矢量图？答：plgridvel(disp)red43、怎样查看INTERFACE？答：plgridifaceidredprintifaceid;id是interface的编号44、FLAC怎么调整视图？答：ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。45、判断模型中网格点是一个还是两个的方法显示模型网格点的ID，如果有两个ID号，说明是两个。节理面上的网格点应该有两个ID号。46、改变模型应力边界时FLAC与3DEC的不同之处在3DEC中，当改变应力边界时，模型表面的应力直接被更换，FLAC中改变应力边界时，模型表面的应力为两者应力相加。47、flac中影响运行时间的因素（1）单元材料，结构面之间的刚度对比很大；（2）划分单元尺寸对比很大，尽量避免长细比大于5：1的细长单元。48FLAC中有多少中材料模型内置11中本构模型%_4[W/N,nR5q（1）零模型用于开挖回填的区域（2）各向同性弹性模型（3）正交各向异性弹性模型（4）横观各向同性弹性模型（5）德鲁克－普拉格塑性模型（6）莫尔－库仑塑性模型一般岩土体的力学行为（7）节理化塑性模型用于研究材料破坏后行为（如边破稳定，屈服煤矿开采）（8）应变硬化/软化模型－库仑塑性模型薄板层状材料破坏后力学行为（9）双线性应变硬化/软化节理化塑性模型临近层状土层开挖（10）双屈服塑性模型（11）修正的剑桥粘土模型49、对不平衡力的理解flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小，该值为百分数，且从不会减小到零，因此，模型不会达到绝对稳定的平横状态，可以根据要求的精度设定R值，模型默认的R值是1e-5。注意R值小仅表示所有网格点的力都平衡，但模型还可能发生塑性流动，通过plotvel命令查看网格点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。50、简单命令解释：plotaddcontdispouton时时显示云图plotblockgroupgpnumon显示节点编号plotinterredske看接触面的效果.delrangegroup1not不删除group信息delrangegroup2删除group信息inizadd数值大小rangegroupname顺z方向移动建立接触面interface1facerangex(范围）y(范围）z(范围）或者interface1facerangeplaneorig(xyz坐标）norm（法线方向）watertableclear清除水位线plsetbackwhite;设置背景setplotbitmapsize(800,510);设置图片大小plotsetcaptionsize37;设置图例数字大小(default=35范围10-50)plotsetcenter100040fixxyzrangex-0.10.1表示在x=0这个平面是固定的fixedboundaryfixxrangex-0.10.1表示在x=0这个平面只x方向不能发生运动rollerboundary51、fish命令的解释：defaux1定义函数loopn(1,nptab1)开始一个循环，从1循环到变量nptab1rr=xtable(1,n)从table1的x列中读第n个数付给变量rrzz=ytable(1,n)从table1的y列中读第n个数付给变量zzloopk(1,nprof)开始第二个循环，从1循环到变量nprofcase_ofk下面是一个选择结构，利用变量k判断iFalfa=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当k＝1的时候执行下面语句xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rrytable(n+10,k)=-40.0tablen+10的y列第k个值等于变量-40case2当k＝2的时候执行下面语句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执行什么定义的函数52.在固定面上施加位移时，fix与apply不能同时使用，应该用ini和f