FLAC3D使用问答.doc
FLAC3D使用问答 看到其它板块上都有这个FAQ,也就是常见问题问答,今天抽了时间进行了整理,想到了就写下来了,因为看到很多初学者费了很多的时间,但是还是没有将常用的命令掌握,所以这个也可以作为入门的初级教材,使大家能够快速的上手,而不用为了某个小命令到处求助。 1. FLAC3D是有限元程序吗 答不是是有限差分法。 2. 最先需要掌握的命令有哪些 答需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 3. 怎样看模型的样子 答plo blo gro可以看到不同的group的颜色分布 4. 怎样看模型的边界情况 答plo gpfix red sk 5. 怎样看模型的体力分布 答plo fap red sk 6. 怎样看模型的云图 答位移plo con dis xdis, ydis, zdis 应力plo con sz sy, sx, sxy, syz, sxz out on 7. 怎样看模型的矢量图 答plo dis xdis, ydis, zdis 8. 怎样看模型有多少单元、节点 r 答plo info 9. 怎样输出模型的后处理图 答File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpe文件 10. 怎样调用一个文件 答File/call或者call命令 11. 如何施加面力 答app nstress 12. 如何调整视图的大小、角度 答综合使用x, y, z, m, Shift键,配合使用CtrlR,CtrlZ等快捷键 13. 如何进行边界约束 答fix x ran (约束的是速度,在初始情况下约束等效于位移约束) 14. 如何知道每个单元的ID 答用鼠标双击单元的表面,可以知道单元的ID和坐标 15. 如何进行切片 答plo set plane ori 点坐标 norm 法向矢量 Plo con sz frontbehind ou on plo con sz plane 显示z方向应力的切片 16. 如何保存计算结果 答save 文件名 17. 如何调用已保存的结果 答rest 文件名;或者File / Restor 18. 如何暂停计算 答ECS 19. 如何在程序中进行暂停,并可恢复计算 答在命令中加入pause命令,用continue进行继续 20. 如何跳过某个计算步 答在计算中按空格键跳过本次计算,自动进入下一步 21. Fish是什么东西 答是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能 22. Fish是否一定要学 答可以不用,需要的时候查Mannual获得需要的变量就可以了 23. FLAC3D允许的命令文件格式有哪些 答无所谓,只要是文本文件,什么后缀都可以 24. 如何调用一些可选模块 答config dyn fluid, creep, cppudm 25 .如何在圆柱体四周如何施加约束条件 可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2R1,其实就是选择range的时候选两个圆柱的差,即得到边界。 26.如何能把一个PLOT的图像数据导出来以便用其他软件绘图 回答用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。 27. 关于dd、dip的理解 dip-direction angle, measured in the global xy-plane clockwise from the positive y-axis 这个是从y轴到x顺时针方向所转过的角度;dip dip angle, measured in the negative z-direction from the global xy-plane 这个是xy平面与z轴负方向的夹角。 dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线 ,dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向,dd 指的是从y轴正方向按顺时针(clockwise)方向转向所要确定面的法线方向在xy平面上的投影的夹角,然后再确定dip,dip指的是从xy平面转向所要确定的平面的角. 28.如何调用FISH文件 新手调用FISH文件出错的原因是文件名为 fishcall.fis.txt因为在记事本里这个txt的后缀是默认的)。修改方法是将 fishcall.fis.txt里面的内容全部复制到新的记事本里,再按保存,保存时候将保存类型改为所有文件*,*,这时你便可在文件名栏内输入 fishcall.fis,类型也就正确了。 29.如何理解Tension-p shear-p tension-n shear-n表示的含义 在计算的循环里面,每个循环中,每个zone都依据failure criterion处于不同的状态,shear 和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态. 30.绘制相关曲线命令的一些解释 1在plot hist m vs n的形式里,m代表y轴,n代表x轴不管m,n的正负; 2 “-“表示对其值作“mirror“,如果大家使用AUTOCAD的话就明白“mirror“的含义了,比如上面的 hist 1 gp ydisp的值是0100,那么vs -1就变成-1000.以此类推. 3 为什么要这么做呢主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内. 31.如何查看剖面上任一点的位移值 print gp disp range x y z 或 print gp disp range id plot gp disp range x y z 或 plot gp disp range id 32. fix x y z range x -0.1 0.1 ..........................a fix x range x -0.1 0.1...........................b 请问a式与b式有何区别 a 表示 在x0 这个平面 是固定的 b 表示在x0 这个平面 只x方向不能发生运动 33.initial 与 apply 有何区别 initial初始化命令,如初始化计算体的应力状态等; apply边界条件限制命令,如施加边界的力、位移等约束等。 Initial的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而apply施加的边界条件不会发生变化。 34. 如何查看各个时段不平衡力的具体数值 采用his来记录计算,包括位移应力等 命令his unbal his gpzone zdis range 0 0 0 或者id 导出数据命令his write n vs m begin 时步end时步 file filename.his n表示纪录的id m表示时步\ Fl 要导出不平衡力的具体数值 his unbal step 100000 or solve his write 1 vs step begin 1 end 1000 file 123.hisw 使用上述命令就可以查看各个时步下的不平衡力的具体数值 35.如何显示变形轮廓线的命令 plo ske magf 10 其中10为放大系数 1、怎样开始入门n 答初步学习首先看说明书 Users gurid的Section2 Getting started和Section 3 Problem sloving with FLAC 2、哪些命令需要掌握 答grid gen prop interface ini fix set hist step solve pause cont sav 等,在学习的过程中遇到这些命令和不会的命令应该查看说明书的Command Reference,里面讲的很清楚,特别需要注意的是一定要仔细看一下Command Reference开头的Optional Range Phrases,这是讲范围的,因为很多命令的施加都需要加到模型的一部分,这里就告诉你怎样来表示这一部分。这一原则同样适合与其他5个ITASCA软件。还有一点就是六个软件(FLAC FALC3D UDEC 3DEC PFC PFC3D)的range大不相同 ,要注意区分。 3、FLAC和FLAC3D建模有什么区别 答FLAC和FLAC3D基本的建模思想大致是一样的,不同的是FLAC是首先用GRID命令把整体的网格画出来,然后用MODELL NULL命令把网格分成许多部分,再用GEN命令把这些部分一块一块搭接起来。然后再用GEN LINE,GEN ARC等命令画出要开挖的边界,在FLAC中产生网格的命令只有一个GRID,而且只能用一次。而FLAC3D是直接用GEN 命令一块一块搭接,GEN可以用很多次。 4、怎样查看模型 Gpd,{{ 答plot grid 可以查看网格,plot grid num 可以查看节点号。 5、怎样查看边界条件| 答plot fix 6、怎样查看预先留设的开挖边界 答plot grid mark。 7、怎么查看模型的云图z 答pl syy sxx xdis ydis sig1.. fill ,注意,带fill关键字就是查看云图,不带fill 关键字就是查看等直线图。 8、怎样查看模型的矢量图 答pl grid vel disp red 9、怎样查看INTERFACE 答pl grid iface id red print iface id ; id 是interface的编号 10、FLAC怎么调整视图 答ctrlz 或shiftz ,然后点击要放大部位的左上角和右下角就可以放大某个部位,如果要缩小,双击屏幕。 1.判断模型中网格点是一个还是两个的方法仿 显示模型网格点的ID,如果有两个ID号,说明是两个。节理面上的网格点应该有两个ID号。 2、改变模型应力边界时FLAC与3DEC的不同之处 在3DEC中,当改变应力边界时,模型表面的应力直接被更换,FLAC中改变应力边界时,模型表面的应力为两者应力相加。 3、对DD,DIP的理解 DD表示节理面正方向在XY面上的投影与Y轴的夹角,也等于节理面与XY面的交线与X轴的夹角; DIP表示节理面与XY面的夹角。 4 塑性状态的解释U -p表示某一区域的应力在模型运行过程中进入过屈服状态,但现在已经推出了屈服状态; -n表示某一区域正处于屈服状态Q 5 对网格点速度的理解 flac中网格点的速度表示的是单位时步网格点发生的位移,查看网格点的速度有两种方法 (1)绘制整个区域的速度图 (2)记录某一网格点的速度 6 flac中影响运行时间的因素 (1)单元材料,结构面之间的刚度对比很大; (2)划分单元尺寸对比很大,尽量避免长细比大于51的细长单元。Ww 7 阻尼的选择 (1)对于静力分析,局部阻尼比较有效,是默认状态的SET MECHANIC DANM LOCAL。 (2)组合阻尼,比局部阻尼更容易在大量网格点发生刚体运动时更好的收敛到稳定状态,使用SET MECH DAMP COMBINE来调用组合阻尼,组合阻尼对于减小动能方面不如局部阻尼有效。 8 FLAC中有多少中材料模型b内置11中本构模型 (1)零模型 用于开挖回填的区域 (2)各向同性弹性模型 (3)正交各向异性弹性模型 4)横观各向同性弹性模型 (5)德鲁克-普拉格塑性模型 (6)莫尔-库仑塑性模型,一般岩土体的力学行为 (7)节理化塑性模型,用于研究材料破坏后行为(如边破稳定,屈服煤矿开采) (8)应变硬化/软化模型-库仑塑性模型 薄板层状材料破坏后力学行为 (9)双线性应变硬化/软化节理化塑性模型临近层状土层开挖 (10)双屈服塑性模型 (11)修正的剑桥粘土模型tSw 9 初始化应力时为什么用了ini命令后还用apply命令,ini命令给模型所有区域应力赋予改定的值,但是它不能保证应力是平衡的。这导致了至少存在两个潜在的问题,第一应力可能违反该区域设定的非线性模型屈服准则,step命令运行后,模型会出现塑性流动。第二,网格边界的指定应力可能和假定的初始应力不相等,当step命令执行后网格点可能移动。为保证应力确实达到初始状态,必须使用ini和apply命令。