FLAC3D学习版-网址.doc

FLAC3D学习1.0版（SIMWE论坛君之枫整理）网上FLAC3D学习笔记建议：初学者将FLAC/FLAC3D版所有帖子都浏览一遍；学有所得后，再浏览一遍，会发现又有新的收获。第一部分（相关链接）1. FLAC3D知识基本介绍SimWe/forum/viewthread.php?tid=209662/forum/viewthread.php?tid=573644/forum/viewthread.php?tid=295354&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=164959&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=366726&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=164959&extra=page%3D6%26filter%3Dtype%26typeid%3D382.FLAC使用FAQ/forum/thread-738220-1-1.html/forum/thread-738329-1-17.html3.复杂模型建立的前处理（这么多款，总有一款适合你吧，可以说导入问题已经解决，看到此资料者，以后不要再发求接口程序帖）真推荐采用有限元软件进行复杂模型建模，网格划分，保存几何信息，导入FLAC3D。有限元软件建模可以到论坛相应版块学习。以下是有限元软件与FLAC3D的接口程序下载地址：/forum/viewthread.php?tid=724033&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=344586&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=310922&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=725347&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=629472&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=657503&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=220759&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=565490&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=548252&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=337989&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=363422&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=529506&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=491786&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=487727&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=295058&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=491786&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=355036&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=342959&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=308078&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=308038&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=302176&highlight=gocad%2Bflac4.FLAC3D后处理个人以为用tecplot进行FLAC3D结果的后处理，看下面两个帖子就可以速成了。/forum/viewthread.php?tid=283356&extra=page%3D1%26filter%3Dtype%26typeid%3D37/forum/thread-722750-1-1.html5. flac3d 命令流语法着色（值得一试）/forum/viewthread.php?tid=624862&highlight=ansys%2B%2Bflac6.初始地应力场问题（能对初始应力场问题有个初步了解）/forum/ viewthread.php?tid=208343&highlight=ansys%2B%2Bflac /forum/viewthread.php?tid=482753&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=750832&highlight=%B3%F5%CA%BC%D3%A6%C1%A6/forum/viewthread.php?tid=732925&highlight=%B3%F5%CA%BC%D3%A6%C1%A6/forum/viewthread.php?tid=356334&highlight=%B3%F5%CA%BC%D3%A6%C1%A6/forum/viewthread.php?tid=388309&highlight=%B3%F5%CA%BC%D3%A6%C1%A6/forum/viewthread.php?tid=316397&highlight=%B3%F5%CA%BC%D3%A6%C1%A6/forum/viewthread.php?tid=316397&highlight=apply%2Bini7.VB调用FLAC3D（FLAC与高级语言间的相互调用，值得参考）/forum/viewthread.php?tid=407677&highlight=ansys%2B%2Bflac/forum/viewthread.php?tid=509327&highlight=ansys%2B%2Bflac8最大主应力和最小主应力分布问题的工程实际意义/forum/viewthread.php?tid=352306&highlight=%D7%EE%B4%F3%D6%F7%D3%A6%C1%A69. 关于dd、dip及origin的问题（仔细看，琢磨琢磨，是很搞人，推荐用ori( ) nor( )的方式确定一个面，即高数书中点法式确立面的方法）/forum/viewthread.php?tid=735924&highlight=dip/forum/viewthread.php?tid=736917&highlight=dip/forum/viewthread.php?tid=728832&highlight=dip/forum/viewthread.php?tid=714259&highlight=dip/forum/viewthread.php?tid=283404&highlight=dip/forum/viewthread.php?tid=217911&highlight=dip10.自定义本构模型的开发/forum/viewthread.php?tid=630517&extra=page%3D11%26filter%3Dtype%26typeid%3D3811. 空间任意平面上固定约束的施加/forum/viewthread.php?tid=375156&highlight=dip12. 如何显示最大位移点的坐标？/forum/viewthread.php?tid=691171&page=1#pid69133913.几个快捷键命令/forum/viewthread.php?tid=564394&extra=page%3D7%26filter%3Dtype%26typeid%3D3814. 查看各个时段不平衡力的具体数值/forum/viewthread.php?tid=640573&extra=page%3D11%26filter%3Dtype%26typeid%3D3815. 在flac3D中gp_head,gp_next,zone_head 等类似的变量是怎么定义的/forum/viewthread.php?tid=543523&extra=page%3D13%26filter%3Dtype%26typeid%3D3816. 如何在剖面上显示各种应力和位移/forum/viewthread.php?tid=754532&extra=page%3D1%26filter%3Dtype%26typeid%3D3717.手册中渗流的例子（cjiao2000版主有详细说明）/forum/viewthread.php?tid=520548&highlight=apply%2Bini18. 塑性铰的问题与讨论/forum/viewthread.php?tid=727531&highlight=apply%2Bini19. 将多条曲线画在同一坐标系的方式/forum/viewthread.php?tid=410552&highlight=apply%2Bini20. 桩的参数指定问题/forum/viewthread.php?tid=406335&extra=page%3D1%26filter%3Dtype%26typeid%3D121. 手册中的PILE单元的中文翻译/forum/viewthread.php?tid=485105&extra=page%3D2%26filter%3Dtype%26typeid%3D122. geogrid单元的部分中文翻译/forum/viewthread.php?tid=488756&extra=page%3D2%26filter%3Dtype%26typeid%3D123. 结构单元中SLIDE的作用和计算方法详述/forum/viewthread.php?tid=730846&extra=page%3D3%26filter%3Dtype%26typeid%3D124. 隧道初期支护中的“钢拱架”应该如何建立/forum/viewthread.php?tid=726231&extra=page%3D3%26filter%3Dtype%26typeid%3D125. 流固耦合手册现场翻译直播/forum/viewthread.php?tid=538175&extra=page%3D3%26filter%3Dtype%26typeid%3D94&page=126.如何输出某一节点随时间的孔压变化/forum/viewthread.php?tid=460705&extra=page%3D6%26filter%3Dtype%26typeid%3D9427. 中文流固耦合资料/forum/viewthread.php?tid=483140&extra=page%3D6%26filter%3Dtype%26typeid%3D9428. 动力分析入门：一个简单的动力分析模型/forum/viewthread.php?tid=626856&extra=page%3D2%26filter%3Dtype%26typeid%3D22229. 结构单元的动力响应分析/forum/viewthread.php?tid=606447&extra=page%3D3%26filter%3Dtype%26typeid%3D22230. 动力时间步究竟如何计算？/forum/viewthread.php?tid=612563&extra=page%3D3%26filter%3Dtype%26typeid%3D22231. 共享一个隧道计算的例子/forum/viewthread.php?tid=742021&extra=page%3D2%26filter%3Dtype%26typeid%3D123第二部分（tg0215提供）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开头的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、怎样查看模型？答：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) red9、怎样查看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 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 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 fillgen 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 fric 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 -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 12.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 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）检查你的模型是不是没有建好，相互之间是不是没有连接好；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 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里处理一下，然后用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 reflect 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) 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 layer3 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 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 (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) dist 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 range 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 zone 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 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 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 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 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 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 ks 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 28initial 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.8e6 grad 0 0 2.6e4 range z 0 28回答：把上面的sxx和syy的剃度2.6e4改为9e3试一下，xx和yy方向的梯度不可能与zz方向的梯度一样的。7.关于ff边界设置的问题自由场地边界（FF边界）在动力分析中经常采用，它的原理是采用粘滞阻尼器来模拟静止边界，将自由场节点的不平衡力加到主体的网格边界上。注意：模型底部的动力边界条件应当在app ff之前，进行app ff以后底部的边界条件就自动转为自由场地边界。现举出一简单例子,比如静力模型的尺寸为10*10*10,有以下几种边界施加方式1）free x y z rang z -.1 .1apply squiet dquiet nquiet rang z -.1 .1apply nstress value hist wave rang z -.1 .1apply ff2）free x y z rang z -.1 .1apply ffapply squiet dquiet nquiet rang z -.1 .1apply nstress value hist wave rang z -.1 .13）free x y z rang z -.1 .1apply squiet dquiet nquiet rang z -.1 .1apply ffapply nstress value hist wave rang z -.1 .1可以验证,以上的施加方式得到的结果是完全一样的.所谓的自由场边界最重要的特征就是复制特征,现对以上各例子做简要分析.1）第二三两句就是在10*10的底面范围内分别设置静止边界和动力条件,当加了ff后,ff的网格的条件即从临近的maingrid的zone进行复制.所以ff的底面的位移边界和动力和maingrid完全一样.2）第二句直接加ff边界,注意这时ff底面是free的,而且动力为0,但是接下来的3.4两句仅仅给定的是z=-.1至.1的范围,导致ff的底面也包括在内了,所以作用方式和命令流1完全相同.3）同2的分析,同样可以得到结论但是对于以上三个命令流中的两句apply squiet dquiet nquiet rang z -.1 .1apply nstress value hist wave rang z -.1 .1如果在maingrid内也指定了范围的话,情形就不一样了,比如2命令流改为apply squiet dquiet nquiet rang z -.1 .1 x 0 10 y 0 10apply nstres