网上FLAC3D学习笔记.doc

FLAC3D学习1.0版（SIMWE论坛君之枫整理） 网上FLAC3D学习笔记 建议初学者将FLAC/FLAC3D版所有帖子都浏览一遍；学有所得后，再浏览一遍，会发现又有新的收获。 第一部分（相关链接） 1. FLAC3D知识基本介绍SimWe http// http// http// id 是interface的编号 10、FLAC怎么调整视图 答ctrlz 或shiftz ,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。 第三部分（原网上学习笔记部分） 1.请问在圆柱体四周如何施加约束条件 可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2r1，其实就是选择range的时候选两个圆柱的差，即得到边界。命令流如下 new gen zone cyl p1 1 0 0 p2 0 4 0 p3 0 0 1 size 12 30 12 gen zone refle norm 1 0 0 gen zone refle norm 0 0 1 mo ss prop den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10 prop ftab 1 ctab 2 dtab 3 table 1 0 45 0.05 42 0.1 40 1 40 table 2 0 2e6 0.05 1e6 0.1 5e5 1 5e5 table 3 0 10 0.05 3 0.1 0 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 回答加上的约束不对 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 摩尔库仑模型 model mohr macro niyan bulk 3.2825e3 shear 1.5150e3 fric 20 coh 0.07 macro meiceng bulk 2.1735e3 shear 8.8917e2 fric 19 coh 0.05 prop niyan range z -13.5 0 prop meiceng range z -15.5 -13.5 prop 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.999 app sxx-20 range x -0.001 0.001 app sxx-20 range x 29.999 30.001 app syy-20 range y -0.001 0.001 app syy-20 range y 5.999 6.001 app szz-20 range z -0.001 0.001 solve （1） model null range in model null range hangdao range name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25 range name hangdaoliangce1 x 12.749 12.751 y 0 6 z -16.50 -15.75 range name hangdaoliangce2 x 17.249 17.251 y 0 6 z -16.50 -15.75 range name hangdao x 12.75 17.25 y 0 6 z -16.501 -16.499 apply 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 90 dd 0 range hangdaoliangce2 ; 给巷道直墙部分（右侧）加压力 apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce1;对点进行监测 hist unbal hist 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）检查你的模型是不是没有建好，相互之间是不是没有连接好； 2）仔细检查你的参数 ，是不是已经发生了破坏。 回答力施加有问题，隧道底部隆起破坏，所以无法计算 问我这个模型就是该模型就是要出现巷道底臌现象，如果按照你的意思，力应该怎样施加才对 回答如果无法收敛，应该用step命令，关键问题是要运行多少步。 3.问 命令流如下 g 10 10 m e def distort loop i 2,jzones loop j 2,jzones xi,j xi,jmag*urand-0.5 yi,jyi,jmag*urand-0.5 end_loop end_loop end prop d 1000 sh 1e8 bu 2e8 其中 xi,j xi,jmag*urand-0.5中的urand是否是随机生成数的命令 回答 urand random number drawn from uni distribution between 0.0 and 1.0； grand random number drawn from normal distribution mean 0.0; standard deviation 1.0 4. 问怎么能把一个PLOT的图像数据导出来以便用其他软件绘图 回答用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。 5. 问用命令建立模型后，如何显示点的坐标 回答使用plo blo gro gpnum on命令 6.interface问题 问我加了一个层面，初始平衡就出现了这么大的位移是怎么原因 ;No1 gen zone brick p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 8 interface 1 face range plane norm 0,0,1 origin 0,0,8 dist 0.01 ini z add -1.0 range group layer1 ;No2 gen zone brick p0 0 0 11 p1 10 0 11 p2 0 10 11 p3 0 0 13 interface 2 face range plane norm 0,0,1 origin 0,0,10 dist 0.01 ini z add -1.0 range group layer2 ;No3 gen zone brick p0 0 0 13 p1 10 0 13 p2 0 10 13 p3 0 0 14 interface 3 face range plane norm 0,0,1 origin 0,0,12 dist 0.01 ini z add -1.0 range group layer3 ;No4 gen zone brick p0 0 0 14 p1 10 0 14 p2 0 10 14 p3 0 0 15 No5 gen zone brick p0 0 0 15 p1 10 0 15 p2 0 10 15 p3 0 0 17 No6 gen zone brick p0 0 0 17 p1 10 0 17 p2 0 10 17 p3 0 0 18 No7 gen zone brick p0 0 0 18 p1 10 0 18 p2 0 10 18 p3 0 0 20 interface 7 face range plane norm 0,0,1 origin 0,0,17 dist 0.01 ini z add -1.0 range group layer7 ;No8 gen zone brick p0 0 0 20 p1 10 0 20 p2 0 10 20 p3 0 0 24 interface 8 face range plane norm 0,0,1 origin 0,0,19 dist 0.01 ini z add -1.0 range group layer8 ;No9 gen zone brick p0 0 0 24 p1 10 0 24 p2 0 10 24 p3 0 0 29 interface 9 face range plane norm 0,0,1 origin 0,0,23 dist 0.01 ini z add -1.0 range group layer9 ;----------边界条件------------------------ fix x range x -20.1 -19.9 fix x range x 19.9 20.1 fix y range y -0.1 0.1 fix y range y 9.9 10.1 fix z range z -0.1 0.1 ;-----------初始条件------------------------- apply szz -26e6 range z 27.9 28.1 initial szz -1.08e6 grad 0 0 2.6e4 range z 0 28 set gravity 0 0 -10 initial density 2.6e3 initial sxx -45.8e6 grad 0 0 2.6e4 range z 0 28 initial syy -45.8e6 grad 0 0 2.6e4 range 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 range 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 s