网上flac3d学习笔记.docx

1、网上flac3d学习笔记网上FLAC3D学习笔记FLAC2d命令的FAQ：1、怎样开始入门？答：初步学习首先看说明书 Users gurid的Section2 Getting started和Section 3 Problem sloving with FLAC2、哪些命令需要掌握？答：grid gen prop interface ini fix set hist step solve pause cont sav 等，在学习的过程中遇到这些命令和不会的命令应该查看说明书的Command Reference，里面讲的很清楚，特别需要注意的是一定要仔细看一下Command Reference开

2、头的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等命令画

3、出要开挖的边界，在FLAC中产生网格的命令只有一个GRID，而且只能用一次。而FLAC3D是直接用GEN 命令一块一块搭接，GEN可以用很多次。4、怎样查看模型？答：plot grid 可以查看网格，plot grid num 可以查看节点号。5、怎样查看边界条件？答：plot fix 6、怎样查看预先留设的开挖边界？答：plot grid mark。7、怎么查看模型的云图？答：pl syy (sxx xdis ydis sig1.) fill ，注意，带fill关键字就是查看云图，不带fill 关键字就是查看等直线图。8、怎样查看模型的矢量图？答：pl grid vel (disp) red

4、9、怎样查看INTERFACE？答：pl grid iface id red print iface id ; id 是interface的编号10、FLAC怎么调整视图？答：ctrl+z 或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。 1.请问在圆柱体四周如何施加约束条件？可以用fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2r1，其实就是选择range的时候选两个圆柱的差，即得到边界。命令流如下：newgen zone cyl p1 1

5、0 0 p2 0 4 0 p3 0 0 1 size 12 30 12gen zone refle norm 1 0 0 gen zone refle norm 0 0 1mo ssprop den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10prop ftab 1 ctab 2 dtab 3table 1 0 45 0.05 42 0.1 40 1 40table 2 0 2e6 0.05 1e6 0.1 5e5 1 5e5table 3 0 10 0.05 3 0.1 0fix x range end1 1 0 0 end2 1 4

6、 0 rad 1 end1 0.5 0 0 end2 0.5 4 0 rad 0.5回答：加上的约束不对！fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 0.5 0 0 end2 0.5 4 0 rad 0.5回答：上面的rad 是不是相差太大了？是不是在1左右来限制了2.关于计算不收敛问题;初始网格gen zone radcyl p0 15 0 -15.75 p1 add 15 0 0 p2 add 0 6 0 p3 add 0 0 15.75 &size 3 4 8 5 dim 2.25 2.25 2.25 2.25 rat 1 1 1 1.1 f

7、illgen zone radtun size 3 4 3 5 dim 0.75 2.25 0.75 2.25 rat 1 1 1 1.1 &p0 15 0 -15.75 p1 add 0 0 -14.25 p2 add 0 6 0 p3 add 15 0 0 fillgen zone reflect dip 90 dd 270 ori 15 0 -15.75; 摩尔库仑模型model mohrmacro niyan bulk 3.2825e3 shear 1.5150e3 fric 20 coh 0.07macro meiceng bulk 2.1735e3 shear 8.8917e2 f

8、ric 19 coh 0.05prop niyan range z =-13.5 0prop meiceng range z =-15.5 -13.5prop niyan range z = -30 -15.5;为开挖定义整体部分range name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25 range name hangdao x 12.75 17.25 y 0 6 z -16.50 -15.75;进行加载fix x y z range z -30.001 -29.999app sxx=-20 range x -

9、0.001 0.001app sxx=-20 range x 29.999 30.001app syy=-20 range y -0.001 0.001app syy=-20 range y 5.999 6.001app szz=-20 range z -0.001 0.001 solve（1）model null range in model null range hangdaorange name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25range name hangdaoliangce1 x 12.749 1

10、2.751 y 0 6 z -16.50 -15.75range name hangdaoliangce2 x 17.249 17.251 y 0 6 z -16.50 -15.75range name hangdao x 12.75 17.25 y 0 6 z -16.501 -16.499apply nstress -0.05 plane dip 90 dd 0 range in ; 给巷道底部加压力app nstress -0.05 plane dip 90 dd 0 range hangdao ; 给巷道直墙部分（左侧）加压力apply nstress -0.05 plane dip

11、90 dd 0 range hangdaoliangce2 ; 给巷道直墙部分（右侧）加压力 apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce1;对点进行监测hist unbalhist gp xdisp 12.75 0 -15.75:监测左帮hist gp zdisp 15 0 -13.5 ;监测顶板hist gp zdisp 15 0 -16.5 ;监测底板solve问：第一个solve平衡了,是第二个solve不能平衡，这是为什么？如果我把支护力加大，也不能平衡？回答：我遇到的不收敛的原因，有两个：1）检查你的模型是不是

12、没有建好，相互之间是不是没有连接好；2）仔细检查你的参数 ，是不是已经发生了破坏。回答：力施加有问题，隧道底部隆起破坏，所以无法计算问：我这个模型就是该模型就是要出现巷道底臌现象，如果按照你的意思，力应该怎样施加才对？ 回答：如果无法收敛，应该用step命令，关键问题是要运行多少步。3.问：命令流如下：g 10 10m edef distortloop i (2,jzones)loop j (2,jzones)x(i,j) =x(i,j)+mag*(urand-0.5)y(i,j)=y(i,j)+mag*(urand-0.5)end_loopend_loopendprop d 1000 sh

13、1e8 bu 2e8其中：x(i,j) =x(i,j)+mag*(urand-0.5)中的urand是否是随机生成数的命令？回答：urand ：random number drawn from uniform distribution between0.0 and 1.0；grand ：random number drawn from normal distribution: mean = 0.0; standard deviation = 1.04. 问：怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？回答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用

14、surfer或者什么作图软件绘制就行了。5. 问：用命令建立模型后，如何显示点的坐标？回答：使用plo blo gro gpnum on命令6.interface问题问：我加了一个层面，初始平衡就出现了这么大的位移？是怎么原因？;No1gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 8 &size 8 8 6 rat 1.02 1.02 1.02gen zone brick p0 0 0 9 p1 10 0 9 p2 0 10 9 p3 0 0 11 &size 8 8 2 rat 1.02 1.02 1.02gen zone reflec

15、t dip 90 dd 90 origin(0,0,6) range z 0 11group layer1 range z 9 11;interface 1 face range plane norm (0,0,1) origin (0,0,8) dist 0.01ini z add -1.0 range group layer1;No2gen zone brick p0 0 0 11 p1 10 0 11 p2 0 10 11 p3 0 0 13 &size 8 8 2 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,1)

16、 range z 11 13group layer2 range z 11 13;interface 2 face range plane norm (0,0,1) origin (0,0,10) dist 0.01ini z add -1.0 range group layer2;No3gen zone brick p0 0 0 13 p1 10 0 13 p2 0 10 13 p3 0 0 14 &size 8 8 1 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,1) range z 13 14group layer

17、3 range z 13 14;interface 3 face range plane norm (0,0,1) origin (0,0,12) dist 0.01ini z add -1.0 range group layer3;No4gen zone brick p0 0 0 14 p1 10 0 14 p2 0 10 14 p3 0 0 15 &size 8 8 1 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,1) range z 14 15group layer4 range z 14 15interface

18、4 face range plane norm (0,0,1) origin (0,0,13) dist 0.01ini z add -1.0 range group layer4;No5gen zone brick p0 0 0 15 p1 10 0 15 p2 0 10 15 p3 0 0 17 &size 8 8 2 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,15) range z 15 17group layer5 range z 15 17interface 5 face range plane norm (

19、0,0,1) origin (0,0,14) dist 0.01ini z add -1.0 range group layer5;No6gen zone brick p0 0 0 17 p1 10 0 17 p2 0 10 17 p3 0 0 18 &size 8 8 1 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,1) range z 17 18group layer6 range z 17 18interface 6 face range plane norm (0,0,1) origin (0,0,16) dis

20、t 0.01ini z add -1.0 range group layer6;No7gen zone brick p0 0 0 18 p1 10 0 18 p2 0 10 18 p3 0 0 20 &size 8 8 2 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,1) range z 18 20group layer7 range z 18 20;interface 7 face range plane norm (0,0,1) origin (0,0,17) dist 0.01ini z add -1.0 rang

21、e group layer7;No8gen zone brick p0 0 0 20 p1 10 0 20 p2 0 10 20 p3 0 0 24 &size 8 8 3 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,17) range z 20 24group layer8 range z 20 24;interface 8 face range plane norm (0,0,1) origin (0,0,19) dist 0.01ini z add -1.0 range group layer8;No9gen zo

22、ne brick p0 0 0 24 p1 10 0 24 p2 0 10 24 p3 0 0 29 &size 8 8 4 rat 1.02 1.02 1.02gen zone reflect dip 90 dd 90 origin(0,0,0) range z 24 29group layer9 range z 24 29;interface 9 face range plane norm (0,0,1) origin (0,0,23) dist 0.01ini z add -1.0 range group layer9;-边界条件-fix x range x -20.1 -19.9fix

23、 x range x 19.9 20.1fix y range y -0.1 0.1fix y range y 9.9 10.1fix z range z -0.1 0.1;-初始条件-apply szz -26e6 range z 27.9 28.1initial szz -1.08e6 grad 0 0 2.6e4 range z 0 28set gravity 0 0 -10initial density 2.6e3initial sxx -45.8e6 grad 0 0 2.6e4 range z 0 28initial syy -45.8e6 grad 0 0 2.6e4 range

24、 z 0 28;-model mohr;-材料属性-;-No1-pro bulk 8e9 shear 7e9 fric 25 coh 7.5e6 ten 2.6e6 range z 0 8 ;-No2-pro bulk 7.8e9 shear 3.5e9 fric 30 coh 6.5e6 ten 2.4e6 range z 8 10;-No3-pro bulk 7e9 shear 3.0e9 fric 22 coh 5.5e6 ten 2.2e6 range z 10 12;-No4-pro bulk 6e9 shear 2.3e9 fric 26 coh 5.5e6 ten 2.0e6 r

25、ange z 12 13;-No5-pro bulk 6.2e9 shear 2.0e9 fric 28 coh 5.2e6 ten 1.2e6 range z 13 14;-No6-pro bulk 5.8e9 shear 2.5e9 fric 35 coh 6.0e6 ten 1.3e6 range z 14 16;-No7-pro bulk 9e9 shear 3.2e9 fric 30 coh 10.3e6 ten 1.2e6 range z 16 17;-No8-pro bulk 8.1e9 shear 3.1e9 fric 26 coh 8.0e6 ten 1.5e6 range

26、z 17 19;-No9-pro bulk 7.6e9 shear 2.6e9 fric 23 coh 8.5e6 ten 1.8e6 range z 19 23;-No10-pro bulk 7.5e9 shear 2.1e9 fric 26 coh 8.0e6 ten 1.9e6 range z 23 28;int 1 prop ks 8e10 kn 8e10 fric 30 dil 6 tension 2e7 bslip=on;int 2 prop ks 8e10 kn 8e10 fric 30 dil 6 tension 2e7 bslip=on;int 3 prop ks 8e10

27、kn 8e10 fric 35 dil 8 tension 2e7 bslip=onint 4 prop ks 9e11 kn 9e11 fric 32 dil 5 tension 2e7 bslip=onint 5 prop ks 9e11 kn 9e11 fric 32 dil 5 tension 2e7 bslip=onint 6 prop ks 8e11 kn 8e11 fric 33 dil 5 tension 2e7 bslip=on;int 7 prop ks 8e10 kn 8e10 fric 33 dil 5 tension 2e7 bslip=on;int 8 prop k

28、s 8e10 kn 8e10 fric 33 dil 5 tension 2e7 bslip=on;int 9 prop ks 8e10 kn 8e10 fric 33 dil 5 tension 2e7 bslip=onhist unbalstep 5000 apply szz -26e6 range z 27.9 28.1initial szz -1.08e6 grad 0 0 2.6e4 range z 0 28set gravity 0 0 -10initial density 2.6e3initial sxx -45.8e6 grad 0 0 2.6e4 range z 0 28in

29、itial syy -45.8e6 grad 0 0 2.6e4 range z 0 28initial szz -1.08e6 grad 0 0 2.6e4 range z 0 28initial density 2.6e3initial sxx -45.8e6 grad 0 0 2.6e4 range z 0 28initial syy -45.8e6 grad 0 0 2.6e4 range z 0 28这是给的模型加的初始地应力，难道这样加地应力，模型就自动平衡了吗？回答：你用ini平衡，为什么还要计算5000步？问：我的INI是用来施加初始地应力的，难道用INI那样加应力边界条件就不用step 或solve算平衡了吗？ 模型就能自动平衡吗？;-初始条件-apply szz -26e6 range z 27.9 28.1initial szz -1.08e6 grad 0 0 2.6e4 range z 0 28set gravity 0 0 -10initial density 2.6e3initial sxx -45.8e6 grad 0 0 2.6e4 range z 0 28initial syy -45