FLAC3D与FLAC常见命令

FLAC/FLAC3D常规问题的整理1．FLAC3D命令的FAQlakewater整理看到其它板块上都有这个FAQ，也就是常见问题问答，今天抽了时间进行了整理，想到了就写下来了，因为看到很多初学者费了很多的时间，但是还是没有将常用的命令掌握，所以这个也可以作为入门的初级教材，使大家能够快速的上手，而不用为了某个小命令到处求助。1.FLAC3D是有限元程序吗？答：不是！是有限差分法。2.最先需要掌握的命令有哪些？答：需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。3`5A5y$Z2N5R6E$i*L3.怎样看模型的样子？0D%s!R9~!E2QY6e)E;Z&A答：ploblogro可以看到不同的group的颜色分布0g)b7M.T6Hl;C4.怎样看模型的边界情况？#o-d(j8`,v/`,G答：plogpfixredsk.x!q8e7c,E3m*u8P5.怎样看模型的体力分布？答：plofapredsk6.怎样看模型的云图？2R)G#Z$K8o*S答：位移：plocondis(xdis,ydis,zdis)4c2j-l(u-B;n([应力：ploconsz(sy,sx,sxy,syz,sxz)7.怎样看模型的矢量图？.g$^b)G%ON2S8答：plodis(xdis,ydis,zdis)5P/E*_6{J%Z6z/l#u4y8.怎样看模型有多少单元、节点？答：ploinfo:m19.怎样输出模型的后处理图？12.如何进行边界约束？,z:T/N.B&\1a答：fixxran（约束的是速度，在初始情况下约束等效于位移约束）13.如何知道每个单元的ID？+o/j&T!Rd%d:h答：用鼠标双击单元的表面，可以知道单元的ID和坐标!B*Y1`6y7U2`.k'K2G4z14.如何进行切片？答：plosetplaneori(点坐标)norm(法向矢量)ploconszplane(显示z方向应力的切片)5g0b'?4E%k8}*F15.如何保存计算结果？答：save+文件名16.如何调用已保存的结果？答：rest+文件名；或者File/Restor17.如何暂停计算？^2U)F(B-D答：Esc答：无所谓，只要是文本文件，什么后缀都可以23.如何调用一些可选模块？;N5?(o5B+f8p答：configdyn(fluid,creep,cppudm)后注：这个工作很繁琐，需要的时间很多，希望广大网友能够将自己曾经遇到的常见问题在后续跟贴，也为了将这个FAQ进行很好的充实。24.如何在圆柱体四周如何施加约束条件？G(@$x(k%i可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r2r1，其实就是选择range的时候选两个圆柱的差，即得到边界。25.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？.rm6N:R+T%}#R回答：用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。26.关于dd、dip的理解：dip-directionangle,measuredintheglobalxy-planeclockwisefromthepositivey-axis这个是从y轴到x顺时针方向所转过的角度；dipdipangle,measuredinthenegativez-directionfromtheglobalxy-plane这个是xy平面与z轴负方向的夹角。dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线，dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，30.如何查看剖面上任一点的位移值？7p3C#b)y0[;[printgpdisprangex()y()z()或printgpdisprangeidplotgpdisprangex()y()z()或1T:f-F6S(hg&E5Z.~!fplotgpdisprangeid!i*}'?6P7z5}*x&C31.fixxyzrangex-0.10.1..........................afixxrangex-0.10.1................................b请问a式与b式有何区别？(fR-C#Y7W0E.|9w(`a表示在x=0这个平面是固定的b表示在x=0这个平面只x方向不能发生运动9Z!q'K/S3`&v32.initial与apply有何区别？step100000orsolve!N;[D,y8`8ehiswrite1vsstepbegin1end1000file123.his.e4I#t2S1R%gP使用上述命令就可以查看各个时步下的不平衡力的具体数值34.如何显示变形轮廓线的命令？ploskemagf10其中10为放大系数2．FLAC2d命令的FAQ：1、怎样开始入门？a#t9P(B0V,Q9^/X.|4i4s答：初步学习首先看说明书User'sgurid的Section2Gettingstarted和Section3ProblemslovingwithFLAC2、哪些命令需要掌握？答：gridgenpropinterfaceinifixsethiststepsolvepausecontsav等，在学习的过程中遇到这些命令和不会的命令应该查看说明书的CommandReference，里面讲的很清楚，特别需要注意的是一定要仔细看一下CommandReference开头的OptionalRangePhrases，这是讲范围的，因为很多命令的施加都需要加到模型的一部分，这里就告诉你怎样来表示这一部分。这一原则同样适合与其他5个ITASCA软件。还有一点就是六个软件（FLACFALC3DUDEC3DECPFCPFC3D）的range大不相同，要注意区分。3、FLAC和FLAC3D建模有什么区别？答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出来，然后用MODELLNULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。然后再用GENLINE，GENARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。而FLAC3D是直接用GEN命令一块一块搭接，GEN可以用很多次。5u9Z3M8@&I*nyO%D'u4、怎样查看模型？&n;P'`6p)A答：plotgrid可以查看网格，plotgridnum可以查看节点号。&q3@%Cw\$J/])]5、怎样查看边界条件？2D*o)a1K%~1Q答：plotfix6B)^&o5[*D6m!F0]6、怎样查看预先留设的开挖边界？%Y#q.d.o#R;E;|*S$f答：plotgridmark。7、怎么查看模型的云图？答：plsyy(sxxxdisydissig1..)fill，注意，带fill关键字就是查看云图，不带fill关键字就是查看等直线图。8、怎样查看模型的矢量图？5v3o0q8W%A2j{.j6X'b答：plgridvel(disp)red0D;e-@,N4v7]3o9、怎样查看INTERFACE？8Q'r;`.t,D(U}0g答：plgridifaceidred+z%_(t5n3Pprintifaceid;id是interface的编号10、FLAC怎么调整视图？答：ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。3FLAC使用步骤FLAC使用步骤1.FLAC程式使用前准备步骤步骤1：依比例画出所欲分析的资料在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。步骤2：换算输入资料成同一单位将现有地层资料，如Density,Bulkmodulus,Young`smodulus,tension,cohesion,frictionAngle等资料，换算成同一单位。附注：需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会造成FLAC无法计算，而产生ERROR信息。步骤3：应用公式简略计算应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。步骤4：建立x,y坐标与nodei,j之间的关系在图上距离及深度的关系，建立x,y坐标系统，再由x,y坐标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格的位置及其相对应的x,y坐标。建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行程序时，不宜使用网格太大的网格数目，因尽量使网格总数少于1000，以节省时间。2.FLAC3D程序的编写步骤：1Config________2Grid________3Model________4求起始的应力平衡(1)建立x,y坐标与网格的关系，建议使用Gen指示：Genx1,y1x2,y2,x3,y3x4,y4i=i0,i1j=j0,j1详细指令参见使用手册，FLAC程式可自动产生x,y坐标与网格的关系，但由于产生的网格坐标不易控制，将对其它的工作产生负面影响，故依使用前步骤4所建立的关系，将网格依其疏密程度需要的不同，实际控制网格的坐标。(2)设定材料性质：prop(3)设定外力：SetGrav,ApplyPressure,inisxx,Syy(4)设定边界条件：fix,free(5)求起始的应力平衡：solve(6)储存：Save5求工程的影响求出区域内的应力分布情况后，再依工程的流程及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤：(1)调出起初的应力平衡：re_____.sav(2)设定新的材料性质：model,prop(3)设定新的支撑性质：struct(4)设定新的外力apply(5)设定边界条件apply,fix(6)求工程时的应力平衡solve,step,cycle(7)储存savefile.sav4．Flac3d的命令1.在FLAC3D中用到的是体积模量K和剪切模量G，它们与杨氏模量E与泊松比υ之间的关系为：)1(2,)21(3EGEK。2.在FLAC3D中边坡的稳定性计算原理仍为折减系数法:kkcc折减前折减后折减前折减后tantan,,其中属于临界值的K即为边坡的安全系数。3.边坡稳定性的判断方法之一为最大不平衡力的收敛性，最大不平衡力的定义为：在每一个计算循环中，外力通过网格节点传递分配到体系各节点时，所有节点的外力与内力之差中的最大值。最大不平衡小于某一临界值作为边坡的收敛标准。命令流为：setmechforcevalue。4.由于FLAC3D在出图背景是灰色，图形为彩色，在论文打印时若我们需要黑白或者灰度的图形，可以使用快捷键CTRL+G或者在图形状态下执行settings/grayscale命令，将图形转换为灰度显示。5.动画的生成，flac3d可以生成动画，大致步骤为，首先命令要显示的图像，设置动画，动画开启，计算。命令流为：Inixd0yd0zd0；将各个方向的位移清零Inixv0yv0zv0；将各个方向的应力清零Appnstress-100e3ranz2.93.1x12y12；在指定体内添加应力Plotsetrot20030；将模型转动一定的可视角度Plotconszzouonmagf10;将竖向应力云图放大十倍以后输出Plotaddhist1；在动画中加入最大不平衡力，编号为1Setmovieavistep1file1.avi；定义动画，动画记录的频率为1时步一频率Moviestart；动画开始Solve；计算