FLAC2D模拟代码解释.docx

1、FLAC2D模拟代码解释1 实例new title 新建命名grid 6,5 ;定义网格数目gen -15,-16 -15,20 60,20 60,-16 ;四角点坐标gen 0 0 0 10 100 10 100 0 i=1,201 j=1,11 ;把大的范围浓缩到小范围内prop bulk=1e9 s=5e8 d=2500 f=30 c=5e5;定义物质属性prop bulk=1e8 shear=.3e8 fric=35 ;体积模量、剪切模量、内摩擦角prop dens=1000 coh=1e10 ten=1e10 ;密度、粘聚力、抗拉强度interface 1 aside from 1,

2、10 to 20,10 bside from 1,11 to 20,11 ;定义节理，将节理放大，随后跟节理的属性m n j=11 ;interfaces,把j=11断开gen line 50.0 21.5 53.5 21.5 ;画线m ss j 33 35prop bulk=10e9 s=4.6e9 d=1400 f=30 c=5e5 i=12,20 j=33,35 ;赋予这个范围属性prop ftab 1 ctab 2 j 33 35tab 1 0 27 0.01 25 1.0 23tab 2 0 5e5 0.01 3e5 1.0 2e5 ;从m ss开始是一组固定语言，应变软化模型先不学

3、m n region i=10 20 j=10 20 ;开挖m n region 102 34 ;开挖，跟前面的gen line 组合使用，开挖一个块就会挖整个m n region 102 30 ;不定义gen line 会开挖一层ini syy=-13e8 ;初始设置，y方向为压应力fix y=1 ;固定y=1这条边界，不变形set g=9.81 ;定义重力加速度set aspect a ;调整屏幕的高宽比，使输出不变形set large ;大变形，坐标更新set small ;小变形，坐标不更新set log on off ;在当前盘上打开flac.log的文件，记录FLAC的运行进程se

4、t timeset step ;设置时间与步数的限制app syy=-12.23e6 j=129 ;使y=129边界（即上边界）受向下的力hist nstep=10 ;每十步记录一次HIS UNBAL ;设置不平衡力HIS syy i=100 j=68 ;记录此点的三向应力变化，默认编号为1HIS sxx i=100 j=68 ;HIS szz i=100 j=68 ;str ;支护，后跟属性或支护材料step ;计算步数solve ;计算solve step ;时步数目，缺省500时步solve time ;运行分钟数，缺省5分钟solve force ;不平衡力，缺省100pause ;暂

5、停continue ;继续save ;保存m m ;model mohr 摩尔模型seg ;分节点node 1 to node ;托盘的节点，分两段就有三个节点pl bou stru beam node red ;即可查看托盘节点号plot grid number ;标注坐标在格子中间plot grid gnum ;标注在点上，查看时用的，简化为pl g num 或gnuplot grid num ;可以查看节点号。plot fix ;查看边界条件plot grid mark ;查看预先留设的开挖边界pl gr if red ;查看interface及编号pl syy(sxx sdis ydi

6、s sigl)fill ;查看模型的云图，带fill 为查看云图，不带fill为查看等值线图pl gr pl ;查看塑性区图pl grid vel(disp)red ;查看模型的矢量图pl grid iface id red ;查看INTERFACE，id为其编号plo gro gnum on ;查看点的坐标plot keyword magfac value ;显示计算后模型的变形状态，value取大于0的set plot dxf 256 ;回车copy 1.dxf ;保存图片title ;定义文件名为“斜坡稳定性分析” Slope Stability Analysisplot grid 即p

7、l g ;查看SClin n x1,y1 x2,y2 ;屏幕上做扫描线，编号为1到5，起始点坐标 plot hist n ;查看位移量，hist的编号2 FLAC2D学习心得newtitleA B C D E (命名文件)grid a b (网格数ab)1，对于大规模建模时，如何利用GEN命令生成局域性的板块如果其中用到到SAME时，是如何完成建模的?1)当用GRID I J 生成网络后,程序默认原点是网格左下角.2)GEN命令,对网格坐标进行重分配时,如果新的分区坐标与GRID生成网格中坐标重合时,用SAME代替2，如何建立起模型的破坏？赋值时如何体现出来3，如何设置不同层次岩土的参数?4,

8、 GEN命令将使网格沿x,y向以1.25传于前一个网跨的比率扩展至1000个单位。注意，如果比例介于01之间，网格的尺寸将随坐标值的增加而减少。例如，发布以下命令： new grid 10，10 m e gen 100，0 100，100 0，100 0,0 rat 0.80，1.25 plot grid 你将会看见沿x,y向逐渐变小的网格。5,移动网点可以用Initial命令，把一个点或若干点从当前位置移到一个新位置。以下各命令建立了一个网格，并用Initial命令来使它变形： grid 5,5 model elastic gen 0,0 0,10 10,10 10,0（沿x,y向规定网点坐

9、标为010） Ini x12 I1 j6（把左上角水平移动2个单元） Ini x12 i6（把各右边界网点各移2个单元，注意，因为没有给j的范围，取为整个范围） 注意，可以用Initial命令，把任何网点移至任何位置。当然，各单元不能重迭。如果重迭了将给出一个为“BAD.GEOMETRY”的警告信息，直到校正了网格结构中的错误为止，flac不会继速执行的。为了得到合理的解，应把各域的高宽比的实际范围，保持在约1:10或更小。6,应力荷载 grid 10,10 m e prop s=5.7e9 b=11.1e9 d=2000 fix x y i=1 fix x y j=1 fix x y i=1

10、1 fix x y j=11 int sxx=30e6 syy=20e6(场地应力) his nstep=1 his xdis i=5 j=5(回采工作与壁处的 x位移) solve ; new excavate M null i=3.4 j=3.6(开挖非1) plot grid step 200(步入平衡) plo his 1(校核平衡) plot grid str (绘应力，网格) mod mohr i=3.4 j=3.6(用sangfill特征) mod null i=3.4 j=7(开挖非2) ini xd=0 yd=0(初始化位移，以致于看出开挖引起的变化) step 200(再次

11、步入接近平衡) plo his 1 plot grid str(绘应力及网格) plot coh(绘不同的材料类型)3 问答. FLAC3D是有限元程序吗？答：不是！是有限差分法.仿2. 最先需要掌握的命令有哪些？www.答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。3. 怎样看模型的样子?5i答：plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况？答：plo gpfix red sk5. 怎样看模型的体力分布?c答：plo fap red sk6. 怎样看模型的云图？答：位移：plo con di

12、s (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点？答：plo info仿9. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件10. 怎样调用一个文件？答：File/call或者call命令11. 如何施加面力？答：app nstress,12. 如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用

13、Ctrl+R，Ctrl+Z等快捷键13. 如何进行边界约束？答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）14. 如何知道每个单元的ID？仿真答：用鼠标双击单元的表面，可以知道单元的ID和坐标仿真15. 如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片)16. 如何保存计算结果？答：save +文件名17. 如何调用已保存的结果？答：rest +文件名；或者File / Restor仿k,Z45rD+D18. 如何暂停计算？答：Esc19. 如何在程序中进行暂停，并可恢复

14、计算？答：在命令中加入pause命令，用continue进行继续仿真仿20. 如何跳过某个计算步?答：在计算中按空格键跳过本次计算，自动进入下一步21. Fish是什么东西？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能22. Fish是否一定要学？答：可以不用，需要的时候查Mannual获得需要的变量就可以了仿23. FLAC3D允许的命令文件格式有哪些？答：无所谓，只要是文本文件，什么后缀都可以24. 如何调用一些可选模块？答：config dyn (fluid, creep, cppudm)25 .如何在圆柱体四周如何施加约束条件？答：可以用fix .

15、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轴到

16、x顺时针方向所转过的角度；dip dip angle, measured in the negative z-direction fromthe global xy-plane 这个是xy平面与z轴负方向的夹角。dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线，dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角.28.如何调用FISH文件？新手调用FI

17、SH文件出错的原因：是：文件名为 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表示pre

18、vious,即分别表示在现在和以前的循环当中处于塑性状态.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

19、 gp disp range idplot gp disp range x() y( ) z( )或plot gp disp range id32.fix x y z range x -0.1 0.1 .aS fix x range x -0.1 0.1.b请问a式与b式有何区别？a 表示 在x=0 这个平面 是固定的b 表示在x=0 这个平面 只x方向不能发生运动33.initial 与 apply 有何区别？initial初始化命令，如初始化计算体的应力状态等；apply边界条件限制命令，如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化

20、，而apply施加的边界条件不会发生变化。34.如何查看各个时段不平衡力的具体数值？采用his来记录计算，包括位移应力等 命令his unbalhis gp(zone) zdis range (0 0 0) 或者id=？导出数据命令his write n vs m begin 时步end时步 file filename.hisn表示纪录的id m表示时步要导出不平衡力的具体数值his unbalstep 100000 or solvehis write 1 vs step begin 1 end 1000 file 123.his使用上述命令就可以查看各个时步下的不平衡力的具体数值35.如何显

21、示变形轮廓线的命令？plo ske magf 10 其中10为放大系数4补充：FLAC2d命令的FAQ：1、怎样开始入门？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent,CFD,CAE,CAD,CAM1giE:R.L Up+K:U答：初步学习首先看说明书 Users gurid的Section2 Getting started和Section 3 Problem sloving with FLAC2、哪些命令需要掌握？R-EC8h ch答：grid gen prop interfa

22、ce 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建

23、模有什么区别？SimWe仿真论坛tR7Jx.V|0|答：FLAC和FLAC3D基本的建模思想大致是一样的，不同的是FLAC是首先用GRID命令把整体的网格画出来，然后用MODELL NULL命令把网格分成许多部分，再用GEN命令把这些部分一块一块搭接起来。然后再用GEN LINE，GEN ARC等命令画出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。而FLAC3D是直接用GEN 命令一块一块搭接，GEN可以用很多次。4、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。仿真分析,有限元,模拟,计算,力学,航空,航天,AN

24、SYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent &ldUun5、怎样查看边界条件？仿真分析,有限元,模拟,计算,力学,航空,航天,ANSYS,MSC,ABAQUS,ALGOR,Adina,COMSOL,FEMLAB,Matlab,Fluent3dJaK答：plot fix 2e6t6_qa:s6、怎样查看预先留设的开挖边界？答：plot grid mark。7、怎么查看模型的云图？答：pl syy (sxx xdis ydis sig1.) fill ，注意，带fill关键字就是查看云图，不带fill 关键字就是查看等直线图。8、怎样查看模型的矢量图？答：pl grid vel (disp) red9、怎样查看INTERFACE？答：pl grid iface id red print iface id ; id 是interface的编号10、FLAC怎么调整视图？tP3d,Y:N-Z答：ctrl+z 或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。