1.1常见问题及其解答
Gen separate 不能被识别
答:原因是FLAC3D版本不行,我用3.0的版本不能。
1. FLAC3D是有限元软件吗?
答:不是,是有限差法软件。
2. FLAC3D最先需要掌握的命令有哪些?
答:需要掌握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)
7. 怎样看模型的矢量图?
答:plo dis (xdis, ydis, zdis)
8. 怎样看模型有多少单元、节点?
答:print info
9. 怎样输出模型的后处理图?
答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpg文件。
10. 怎样调用一个文件?
答:使用菜单File/call或者call命令。
11. 如何施加面力?
答:app nstress ran <……>
12. 如何调整视图的大小、角度?
答:综合使用x, y, z, m, Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键。
13. 如何进行边界约束?
答:fix x ran <……>(约束的是速度,在初始情况下约束等效于位移约束)
14. 如何知道每个单元的ID?
答:使用鼠标双击单元的表面,可以知道单元的ID和坐标。
15. 如何进行切片?
答:plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 16. 如何保存计算结果?
答:save filename(文件名可自定义)
17. 如何调用已保存的结果?
答:使用菜单File/call或者命令rest filename(文件名可自定义)。
18. 如何暂停计算?
答:运行中使用Esc命令。
19. 如何在程序中进行暂停,并可恢复计算?
答:在命令中加入pause命令,键入continue命令后可恢复计算。
20. 如何跳过某个计算步?
答:在计算中按空格键可跳过本次计算,自动进入下一步。
21. FISH是什么?
答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能。22. FISH是否一定要学?
答:可以不用,需要的时候查Manual获得需要的变量就可以了。
23. FLAC3D允许的命令文件格式有哪些?
答:只要是符合FLAC3D格式要求的文本文件,无论是什么后缀名,都可以为FLAC3D调用。24. 如何调用一些可选模块?
答:使用命令config dyn (fluid, creep, cppudm)。
25. 如何使用gauss_dev对符合高斯正态分布的材料参数进行赋值?
答:假定某材料的摩擦角均值为40度,标准差是2,则命令如下:prop friction 40 gauss_dev 2
26. FISH函数中是否能调用“.sav”文件?
答:不能。FLAC3D中规定,new和restore命令不允许出现在FISH函数中,因为new和restore
命令会将原有存储信息清除掉。
27. initial 与apply 有何区别?
答:initial初始化命令,如初始化计算体的应力状态等;apply边界条件限制命令,如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而apply施加的边界条件不会发生变化。
28. FLAC3D动力分析中是如何计算永久变形的?
答:FLAC3D采用动态运动方程求解动力方程,因此采用弹塑性本构模型可以计算永久变形。而土动力学常用的粘弹性模型由于没有考虑土体的塑性,因此不能计算永久变形。
29. 对于初学者而言,是学习FLAC还是FLAC3D?
答:FLAC有较好的图形化操作界面,而FLAC3D目前只能通过命令流来操作,从学习难度上来说,FLAC要简单一些,不过复杂的三维问题还是需要使用FLAC3D才能解决。FLAC和FLAC3D的某些命令和分析方法类似,读者在学习过程中可以相互借鉴。
30. interface建模命令中的dist关键词是否表示接触面的厚度?
答:FLAC3D中的interface是没有厚度的,dist关键词表示的是接触面建模时选择范围时的容差,表示该范围内的“面”上将被赋予interface单元。
31. 初始应力场计算中位移场和速度场是否都要清零?
答:是的。一般,FLAC和FLAC3D中位移场和速度场的清零命令都是同时使用的。
32. 加了fix边界,再使用apply施加应力边界有效吗?
答:无效。fix和apply都是边界条件,两者不能混用,fix的作用是固定节点的速度,只要用户不更改这个速度,在计算中都会保持不变。
33. solve age后面跟随的时间是真实的时间吗?
答:FLAC和FLAC3D在动力、渗流、流变模式下才有真实的时间,时间的单位默认为秒,也可以根据读者使用的量纲进行调整。
34. FLAC3D中主应力大小是怎么规定的?
答:FLAC和FLAC3D中的大小主应力是根据应力的数值大小来规定的,并且规定压为负,而土力学中一般规定压为正,所以FLAC3D中的大小主应力z_sig1(p_z)、z_sig2(p_z)和z_sig3(p_z)分别对应于土力学中的小主应力、中主应力和大主应力,在使用时要注意区别。
35. FISH函数中dof的含义是什么?
答:一些关于结构单元的FISH函数中常常出现dof变量,该变量表示的是自由度,如nd_pos( np, p, dof )函数中dof?{1,2,3}表示结构节点的三个方向的自由度。
36. 怎么在不规则的面上施加水压力?
答:设置合理的水压力梯度和作用范围,使用apply nstress命令即可。因为apply施加的应力边界条件是作用在“边界”上的,所以程序会在用户设置的range范围中自动寻找“边界”,而不管这个“边界”有多么复杂,而且nstress表示力作用的方向是垂直于“边界”,该关键词可以保证水压力的作用方向始终垂直于作用面。
37. hist记录的数据如何转到Excel?
答:使用类似如下的命令:hist write 1 vs 2 file 1.xls 可以将历史记录ID为1和2的对应关系输入到文件1.xls中,然后用Excel打开进行编辑、处理。
1.2常见错误提示及其解决办法
本节汇集了FLAC3D在使用过程中常见的错误提示,并根据不同的提示总结了出错原因和解决办法。FLAC3D程序自身的检查功能不多,但也有一部分错误提示,读者也可以根据软件提示的内容迅速找到错误的原因,并予以修正。
1. Bad type (pointer) conversion
出错原因:在编写FISH函数时,某些变量的赋值错误所致。解决办法:仔细检查FISH函数中的变量赋值情况,尤其注意涉及到指针、FISH自有变量的赋值等语句。
2. Gridpoints 19801 and 19803 have identical coordinates in zone 9703
出错原因:在同一个单元内的两个节点有相同的坐标,这可能是由于将其他软件建立的模型导入到FLAC3D时两个软件的节点坐标精度差异所导致的。解决办法:使用attach face来合并相关节点,或者重新检查模型。
3. Memory allocation error
出错原因:可能是网格划分的过多,超过了计算的内存所致。解决办法;减小网格数量或者加大计算机的内存。
4. Mesh primitives does not conform to node numbering convention
出错原因:在建模时各节点坐标设置的顺序与FLAC3D中基本网格形状不一致。解决办法:检查建模时p0~p12等节点坐标,使其符合FLAC3D的要求。
5. Source node 2 already has a link!
出错原因:在结构单元计算中,对已存在连接的节点进行设置时会出现此类错误。解决办法:检查需要设置连接的结构节点,确保已有连接已被删除才能设置新的连接。
6. Timestep rejected by module
出错原因:一般是由于结构单元的密度没有赋值造成的。解决办法:用命令print shell prop dens来显示结构单元的密度,查看是否所有单元都已经赋值。如有遗漏,应重新赋值。
7. The model name does not exist
出错原因:可能是由于模型名称输入错误,或者在调用某些可选模块(如渗流、动力)的模型时没有设置相应的Config。解决办法:检查模型名称是否输入有误,在可选模块下检查是否已设置相应的Config。
8. Unrecognized parameter 3 (***)
出错原因:命令输入时存在错误的参数,且出错的是命令中的第3个参数。解决办法:检查出错命令的具体位置,找到第3个参数,进行修改。
9. Viscous damping too high
出错原因:在进行UDM编写动力方面的本构时可能会遇到这种错误。解决办法:可能由于粘性函数偏大造成的,时间步的增大会导致粘性函数值的增大。在FLAC3D程序中,如果这个值大于1,那么就会出现这个错误。
10. Zero stiffness in grid-point xxx
出错原因:“0刚度”一般是由于材料参数未正确赋值所致。解决办法:仔细检查计算中的材料赋值命令,看是否有遗漏,如使用以下命令来显示体积模量(bulk)参数的赋值情况。plot block prop bulk
11. Zero volume tet in zone xxx
出错原因:一般是在计算分析中使用了大变形模式(set large)。由于在大变形模式计算过程中,节点坐标会随时步自动,这样有时会导致网格畸形,而无法进行下去。解决方法:慎用大变形模式,大变形模式适用于粘结力较小材料的开挖过程模拟。因此,一般问题的模拟过程,宜采用小变形模式(默认变形模式)进行。即使是大变形问题的初始应力场,也应采用小变形模式生成,再视后续工况的具体情况确定是否改为大变形模式。
1.3学习经验和建议
以作者与大多数FLAC/FLAC3D使用者的交流,总结了几点软件学习方面的经验和建议,希望可以对读者提供帮助。
1. 了解FLAC和FLAC3D的适用范围、优点和局限性
任何一种方法都是有一定的适用范围,并不能不能解决所有问题,这学要读者对所使用工具的优点和局限性有清醒的认识。数值模拟的最终目的是为工程问题的诊断和解决提供服务的,需根据问题的本质选择合适的方法和工具;而非“膜拜”和迷信某种方法,机械地用它去套工程,本末倒置。
2. 由简到繁,循序渐进
遵循“由简到繁,循序渐进”的学习方法,切忌盲目求大求全,期望一口气吃成胖子。学习时,可进行少量单元的简单数值试验来理解软件的特点和功能,积累一定的经验后再进行复杂的数值模拟试验。
3. 充分利用手册
手册是最权威的软件说明书,一定要充分利用。尽管FLAC和FLAC3D的手册编制顺序不一定适合中国读者的思维习惯,但应尽量养成查阅手册的习惯,做到常翻常新。手册中的例子大多都是为了说明某个特定的问题而设定的,因此在讲述该问题时往往会忽略与该问题无关的一些细节,比如参数选择等,因此读者在学习手册时不要“迷信”某个特定的例子,也不要“纠缠”于某些无关的细节,而是要从这些例子中掌握分析问题的基本方法。
4. 了解计算中每一条语句的含义
初学者由于对FLAC和FLAC3D软件了解的不多,在计算时往往会直接套用软件手册或教科书中的例子,而对例子中某些语句的含义并不是真正的了解,这些“不明其意”的语句往往是造成计算结果不合理的原因。这里建议读者在使用FLAC和FLAC3D程序时,要对自己编写的命令文件中的每一条语句都有清晰的认识和了解,这就要求读者要勤查手册、注重平时的积累。
5. 多做“数值试验”
FLAC和FLAC3D程序功能强大,内容众多,在分析具体问题时,读者往往会遇到如法解决的新问题,这些问题在软件手册或教科书中都很难找到答案,这时读者应该多做一些小的算例,开展数值试验,从而了解程序的功能,达到解决问题的目的。
6. 使用“?”
FLAC3D的命令很多,在初学者看来,记住数量可观的各种命令及语句格式是一件很困难的事情,事实也的确如此。幸运的是,FLAC3D在命令窗口中提供了“?”功能,无论在命令的什么位置都可以插入“?”字符,让系统告诉你接下来可以应该输入的是哪些关键字或变量。
7. 夯实知识基础
FLAC和FLAC3D的计算结果和中间时步表现出一些不合实际的结果,需要读者具有足够的专业和数学知识进行判断与解释。因此,决定FLAC和FLAC3D使用水平高低的决定性因素取决于使用者的专业素养、工程经验和数理知识。因此加强专业知识、数学和力学的学习,夯实知识基础十分重要。
8. 相互交流,取长补短
FLAC和FLAC3D命令、关键词和变量繁多,个人学习难免顾此失彼,因此强交流,与他人共__享学习经验是提高FLAC和FLAC3D应用水平的一个捷径。互联网的出现,为大家提供了一个讨论和共享的平台,读者可以在相互间的交流、争论中取长补短,共同提高。
遇到失意伤心事,多想有一个懂你的人来指点迷津,因他懂你,会以我心,换你心,站在你的位置上思虑,为你排优解难。
一个人,来这世间,必须懂得一些人情事理,才能不断成长。就像躬耕于陇亩的农人,必须懂得土地与种子的情怀,才能有所收获。
一个女子,一生所求,莫过于找到一个懂她的人,执手白头,相伴终老。
即使芦花暖鞋,菊花枕头,也觉温暖;即使粗食布衣,陋室简静,也觉舒适,一句“懂你”,叫人无怨无悔,愿以自己的一生来交付。
懂得是彼此的欣赏,是灵魂的轻唤,是惺惺相惜,是爱,是暖,是彼此的融化;是走一段很远的路,蓦然回首却发现,我依然在你的视线里;是回眸相视一笑的无言;是一条偏僻幽静的小路,不显山,不露水,路边长满你喜爱的花草,静默无语却馨香盈怀,而路的尽头,便是通达你心灵的小屋……
瑟瑟严冬,窗外雪飘,絮絮自语说了这多,你可懂我了吗?若你知晓,无需说话,只报一声心灵的轻叹,那,便是我的花开春暖。
你相不相信,人生有一种念想,不求奢华不求结果,不求你在我身边,只愿有一种陪伴暖在心灵,那,便是懂得。
有人懂得是一种幸福,懂得别人是一种襟怀,互为懂得是一种境界。
懂得,真好!
实例分析命令: 1. X ,Y ,Z 旋转 Shift+ X ,Y ,Z 反向旋转 Gen zone ……;model ……;prop ……(材料参数);set grav 0,0,-9.81(重力加速度) plot add block group red yellow 把在group 中的部分染成红色和黄色 plot add axes black 坐标轴线为黑色;print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b (设置初始密度,有时不同层密度不同);ini ……(设置初始条件);fix ……(固定界面) set plot jpg ;set plot quality 100 ;plot hard file 1.jpg 图像输出(格式、像素、名称) plot set magf 1.0视图的放大倍数为1.0;plo con szz z 方向应力云图 2. ini z add -1 range group one 群one 的所有单元,在z 方向上向下移动1m ;然后合并 命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体,1e-5是容差 3. (基坑开挖步骤):Step 1: create initial model state (建立初始模型)Step 2: excavate trench (开挖隧道) 4. group Top range group Base not 定义(群组Base 以外的为)群组Top 5. plot blo gro 使得各个群组不同颜色显示 6. (两个部分间设置界面;切割法):gen separate Top 使两部分的接触网格分离 为两部分;interface 1 wrap Base Top 在(Base 和Top )这两部分之间添加接触单元;plot create view_int 显示,并创建标题view_int ;plot add surface 显示表面;plot add interface red 界面颜色红色 7. (简单的定义函数及运行函数)new ;def setup 定义函数setup ;numy = 8定义常 量numy 为8;depth = 10.0 定义depth 为10;end 结束对函数的定义;setup 运行函数setup 8. (隧道生成)上部圆形放射性圆柱及下部块体单元体的建立,然后镜像。 9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度?— 10. 渐变应力施加:apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60 dd 270 origin .1 0 0;施加法向应力:apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 0 11. d ip dd 确定平面位置使用:(纠结) 12. p rint gp position range id=14647 输出节点坐标 13. a pply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中,z 为变量,所以xx σ为:65=-1010+10xx z σ?? ;原点(0,0,0) 14. f ree x range x -.1 .1 z 6.9 10.1放松x=0 平面上,z=7,10 这一部分在x 方向的约 束(可以在此处产生破坏) 15. 体积模量K 和剪切模量G 与杨氏模量及泊松比v 之间的转换关系如下: =3(1-2v)E K G=2(1+v) E 16. 一般而言,大多数问题可以采用FLAC 3D 默认的收敛标准(或称相对收敛标准),即当体 系最大不平衡力与典型内力的比率R 小于定值10-5;(也可由用户自定义该值,命令:
建模 1、调用文件: ①文件与工程在同一个文件夹,只写文件名即可:Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlyt hefilenameneed beenteredwiththe CALL command. ②不在同一个文件夹,全路径:Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfol der\file.dat). Undo;撤销上一条命令 2、创建旋转缩放视图 3、建模命令 modelmechmohr;莫尔库伦模型 modelmechelastic;弹性模型 setgrav0,0,-9.81;重力加速度negative z-direction.(垂直向下!常用的) 下下面面这这代代码码,,是是沿沿着着--y y方方向向的的重重力力加加速速度度,,注注意意区区别别!!!!!!!! genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续 p110,-10,-20... p2-10,10,-20... p3-10,-10,0 plotzone
genzonebricksize6,8,8p0-10,-10,-20...;不规则六面体 p110,-10,-20p2-10,10,-20... p3-10,-10,0p410,10,-20... p5-10,10,10p610,-10,0... p710,10,10 plotcurrentplotPlot01 plotclear plotzone Undo;撤销命令 setlogfile127G1001.tGt setlogontruncate setlogoff listzoneprinrangeG01y01z01;显示指定范围内各单元的主应力,结果如下 Hist命令: ①命令编号按顺序从1开始:eachhistoryisnumberedsequentiallyfrom1asitisenteredviathe HISTORY co mmand. ②查找显示所有的his命令:ReturntotheFlac3D>promptandtype listhist foralistingofthehistoriesandtheircorrespondingnumbers. histnstep5;每5步记录1次。默认是10步记录1次
1、最先需要掌握的命令有哪些? 答:需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 2、怎样输出模型的后处理图? 答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpe文件。 3、怎样调用一个文件? 答:File/call或者call命令 4、如何施加面力? 答:app nstress 5、如何调整视图的大小、角度? 答:综合使用x, y, z, m, Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键。 6、如何进行边界约束? 答:fix x ran (约束的是速度,在初始情况下约束等效于位移约束)。 7、如何知道每个单元的ID? 答:用鼠标双击单元的表面,可以知道单元的ID和坐标。 8、如何进行切片? 答:plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 9、如何保存计算结果? 答:save +文件名 10、如何调用已保存的结果? 答:rest +文件名;或者File / Restor 11、如何暂停计算? 答:Esc 12、如何在程序中进行暂停,并可恢复计算? 答:在命令中加入pause命令,用continue进行继续。 在我们分步求解中想得到某一个过程中的结果,不用等到全求完,还可以在分布求解错误的时候就进行改正,而不是等到结果出来。 13、如何跳过某个计算步? 答:在计算中按空格键跳过本次计算,自动进入下一步 14、Fish是什么东西?Fish是否一定要学?
答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能。Fish可以不用学,需要的时候查Mannual获得需要的变量就可以了。 15、FLAC3D允许的命令文件格式有哪些? 答:无所谓,只要是文本文件,什么后缀都可以。 16、如何调用一些可选模块? 答:config dyn (fluid, creep, cppudm) 17、如何在圆柱体四周如何施加约束条件? 可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2 Flac3D 中文手册 FLAC3D的计算模式中是否需要做孔压分析取决于是否采用config fluid命令。 1 无渗流模式(不使用config fluid) 即使不使用命令config fluid,仍然可以在节点上施加孔压。这种模式下,孔压将保持为常量。如果采用塑性本构模型的话,材料的破坏将由有效应力状态来控制。 节点上的孔压分布可由initial pp命令或water table命令来设定。如果采用water table命令,由程序自动计算水位线以下的静水孔压分布。此时,必须施加流体密度(water density)和重力(set gravity)。流体密度值和水位位置可以用命令print water显示。如果水位线是由face关键字来定义的,则可用命令plot water命令显示水位。 这两种情况,单元的孔压都由节点孔压值平均求出,并在本构模型计算中用作有效应力。这种计算模式下,体积力中不反映流体的出现:用户必须根据水位线以上或以下相应地指定干密度和湿密度。使用命令print gp pp和priint zone pp可分别得到节点或单元孔压。plot contour pp命令可绘出节点孔压云图。 2 渗流模式(使用config fluid) 如果使用命令config fluid,则可进行瞬时渗流分析,孔压改变和潜水面的改变都可能出现。在config fluid模式下,有效应力计算(静态孔压分布)和非排水计算均被执行。除此之外,还可进行全耦合 分析,这种情况下,孔压改变将使固体产生变形,同时体积应变反过来影响孔压的变化。 如果采用渗流模式,单元孔压仍由节点孔压平均求出。但这种模式,用户只能指定干密度(不论是水位以上还是以下),因为FLAC3D将流体的影响考虑到了体积力的计算中。 采用渗流模式时,渗流模型必须施加到单元上,使用命令model fl_isotropic模拟各向同性渗流,model fl_anisotropic模拟各向异性渗流,model fl_null模拟非渗透物质。注意,力学模型为空的单元并不代表渗流模型为空。 流体性质(参数)可施加到单元或节点上。各向同性渗透率、孔隙率、比奥系数和非排水热系数等单元流体性质由命令property施加。 对于各向同性渗流,渗透率通过perm关键字赋予。对各向异性渗流,渗透率的3个主值采用关键字k1,k2,k3赋予,主方向由关键字fdip,fdd,frot确定。渗透率的主方向服从右手系统。fdip和fdd分别为k1和k2确定的平面的倾向和倾角。frot为k1轴和倾角矢量的旋转角。如果不特别指定,比奥系数默认为1,孔隙率默认为0.5。节点的渗流性质由命令initial指定。这些性质包括流体重度、流体体积模量、比奥模量、流体抗拉强度和饱和度。每种性质在空间上都可以变化。流体重度也可以用water命令给出。 在渗流模式里,有必要知道可压缩性被定义在以下两种参数中:(1)比奥系数和比奥模量;(2)流体体积模量和孔隙率。第一种 FLAC3D常见命令与使用技巧 1、FLAC3D常见命令: 是有限元程序吗答:不是!是有限差分法。 2.最先需要掌握的命令有哪些 答:需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 3.怎样看模型的样子答:plo blo gro可以看到不同的group的颜色分布 4.怎样看模型的边界情况答:plo gpfix red 5.怎样看模型的体力分布答:plo fap red 6.怎样看模型的云图答:位移:plo con dis (xdis, ydis, zdis)应力:plo con sz (sy, sx,sxy, syz, sxz) 7.怎样看模型的矢量图答:plo dis (xdis, ydis, zdis) 8.怎样看模型有多少单元、节点答:pri info 9.怎样输出模型的后处理图 答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpe文件 10.怎样调用一个文件答:File/call或者call命令 10.如何施加面力答:app nstress 11.如何调整视图的大小、角度答:综合使用x, y, z, m, Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键 12.如何进行边界约束答:fix x ran(约束的是速度,在初始情况下约束等效于位移约束) 13.如何知道每个单元的ID答:用鼠标双击单元的表面,可以知道单元的ID和坐标 14.如何进行切片 答:plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 15.如何保存计算结果答:save +文件名. 16.如何调用已保存的结果答:rest +文件名;或者File / Restore 17.如何暂停计算答:Esc 18.如何在程序中进行暂停,并可恢复计算答:在命令中加入pause命令,用continue进行继续 19.如何跳过某个计算步答:在计算中按空格键跳过本次计算,自动进入下一步 20. Fish是什么东西 答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能 1、首先rest123.sav,然后再callFLAC3d2Tecplop生成Tec10.dat 2、打开Tecplot程序,执行【File】/【Load Data File】命令,打开数据文件对话 框,选择刚刚生成的tec10.dat文件,并执行【打开】命令,在出现的Select Initial Plot 对话框中,单击OK按钮。并在Tecplot左上角显示的3维笛卡尔坐标系下选择【3D Cartasion 】 3、在Tecplot中选中【Data】→【Extract】→【Slice from Plane】,打开切 面对话框。在对话框中我们可以选中x=6/y=5/z=7等切片,并单击Extract. 4、在Tecplot中选中【File】→【White Data File】,打开【White Data File Options】对话框,在对话框中选中ASCII、Point。并在下面的Zone(s)选项中选中你所要的2 SLc=3,在Variable(s)中同时选中1 x(m),2 y(m)及你所要的13 szz(pa)等,并单击OK按钮。保存的文件自己最好独自命名(如tec-90.dat格式) 5、以txt格式打开刚才命名的tec-90.dat,把里面的三列数复制到txt里面后 再导入Excel里面,同时在Excel里面把szz的值变成相反的值保存在D列,(txt里面注意不是三列的都不要,只复制三列数据,把前面的和后面的都去掉)。 6、打开surfer,【网格】→【数据】,选中刚才生成的123.xls(即刚才生成的Excel 数据,选中相应的A、B和D列),并单击打开。对于【网格化数据】选项单则单击【确定】,则会生成与123.xls相对应的123.grd(即与刚才生成的Excel 文件名相同)。 7、在【地图】→【表面图】单击,在【打开网格对话框】中选中生成的123.grd 文件,并单击打开。 8、右击图片选中【属性】并打击打开,可以对网格、颜色进行设置(颜色谱选 安装文件夹中的rainbow.CLR即可,具体是在打开的【属性对话框】→【常规】→【材质颜色】→【上部】单击上部中的颜色,在颜色谱中选择【载入】,在打开里面选中安装文件夹里面的rainbow.CLR即可,在多次确认) 9、此外也可以在surfer里面对其他参数进行设置。 快 速 入 门(GETTING STARTED) 版本:flac3d 3.0版(FTD127) 翻译:一米 2009.06 声 明 现在市面上关于FLAC3D软件的教材寥寥无几,在学习的过程中,主要还是参考软件本身的使用手册,虽然读英文版手册有些吃力,但是它论述非常详细,我觉得是用户最好的教材。我在边看手册的时候边做了翻译,目前为止翻译完成了本部分的内容(略去了部分内容和例子),还翻译了命令手册的前半部分内容,等翻译完成了,也会和网友共享,但是像本人这类英语水平一般的人做这样的翻译工作是比较辛苦的,我也不确定是否有毅力完成命令手册下半部分的内容。虽然这样的工作比较艰难,但我觉得还是学到了不少东西,手册是最原始,最翔实的基础教材,看明白了手册,运用软件才会游刃有余。 由于本人专业水平和英语能力的限制,存在问题是在所难免的,有的地方甚至可能曲解了原意。考虑到时间因素,译文的措辞没有细细斟酌,还请网友谅解。如果发现译文中的错误,还请广大读者斧正。 一米 2 快速入门 这一部分将向初次使用flac3d的用户介绍软件的基本使用方法。主要有以下内容:软件的安装与启动;用软件分析解决问题的步骤,在每一步的操作中,都有简单例题来说明该步骤具体是如何操作的。 如果你对软件比较熟悉,但是现在很少用它来处理问题,那么这部分的内容(尤其2.7节)能很好的帮你回顾软件操作的要点。本部分3.3节全面详细的介绍了如何进行问题的求解。 Flac3d支持命令驱动和图形菜单驱动两种模式*。在本手册中大部分的算例都采用了命令驱动模式。我们认为这种模式能给用户提供操作软件最清晰的思路。在1.1节中我们就已经提到了命令驱动模式使得flac3d在分析求解工程问题时成为了一个功能强大的“多面手”。然而这种模式让新用户,或者长时间未接触软件的老用户用起来有点不那么容易。命令行必须用键盘输入,可以直接输入到软件的命令窗口,或者先保存为数据文件,再通过软件的相关命令进行读取。Flac3d能识别超过40个主命令和400多个附属的关键词。 本部分主要包括以下内容: 1 在2.1节,手把手的教你们如何在自己的电脑上安装和启动flac3d软件。 2 在2.2节,用一些简单的教学案例帮组用户熟悉一些常用的命令。 3 在用户建立自己的模型并进行分析计算之前,有必要先了解flac3d的一 些基本知识。在2.3节讲述了flac3d的基本术语;在2.4节主要说明了有 限差分网格的定义规则;而在2.5节阐述了输入命令的基本句法。 4 在2.6节,阐述了flac3d的特点,比如创建、命名和使用对象,以方便 用户进行问题的求解 5在2.7节,一步步的指导用户如何建模和分析问题,每一个步骤都分开论述,并提供简单的例子帮助用户理解。 6 2.8节-2.10节分别论述了系统的符号约定、单位体系和精度限制 7 2.11节说明了软件中各种类型文件的创建和使用。 8 2.12节对图形菜单操作模式进行了简介。 *:对于初级用户来说一般图形菜单驱动模式只进行图形输出或者文件操作。本章节的最后一部分将向用户展示如何使用图形菜单驱动模式来操作软件。 2、2 三维数值模拟方法及其原理 2、2、1 FLAC3D工程分析软件特点 FLAC3D就是由美国Itasca Consulting Group, Inc、为地质工程应用而开发得连续介质显式有限差分计算机软件。FLAC即Fast Lagrangian Analysis of Continua 得缩写。该软件主要适用于模拟计算岩土体材料得力学行为及岩土材料达到屈服极限后产生得塑性流动,对大变形情况应用效果更好。 FLAC3D程序在数学上采用得就是快速拉格朗日方法,基于显式差分来获得模型全部运动方程与本构方程得步长解,其本构方程由基本应力应变定义及虎克定律导出,运动平衡方程则直接应用了柯西运动方程,该方程由牛顿运动定律导出。 计算模型一般就是由若干不同形状得三维单元体组成,也即剖分得空间单元网络区,计算中又将每个单元体进一步划分成由四个节点构成得四面体,四面体得应力应变只通过四个节点向其它四面体传递,进而传递到其它单元体。当对某一节点施加荷载后,在某一个微小得时间段内,作用于该点得荷载只对周围得若干节点(相邻节点)有影响。利用运动方程,根据单元节点得速度变化与时间,可计算出单元之间得相对位移,进而求出单元应变,再利用单元模型得本构方程,可求出单元应力。在计算应变过程中,利用高斯积分理论,将三维问题转化为二维问题而使其简单化。在运动方程中,还充分考虑了岩土体所具有得粘滞性,将其视作阻尼附加于方程中。 FLAC3D具有一个功能强大得网格生成器,有12种基本形状得单元体可供选择,利用这12种基本单元体,几乎可以构成任何形状得空间立体模型。 FLAC3D主要就是为地质工程应用而开发得岩土体力学数值评价计算程序,自身设计有九种材料本构模型: (1)空模型(Null Model) (2)弹性各向同性材料模型(Elastic, Isotropic Model) (3)弹性各向异性材料模型(Elastic, anisotropic Model) (4)德拉克-普拉格弹塑性材料模型(Drucker-Prager Model) (5)莫尔-库伦弹塑性材料模型(Mohr-Coulomb Model) ;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点);沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl FLAC 3D基础知识介绍 一、概述 FLAC(Fast Lagrangian Analysis of Continua)由美国Itasca公司开发的。目前,FLAC有二维和三维计算程序两个版本,二维计算程序V3.0以前的为DOS版本,V2.5版本仅仅能够使用计算机的基本内存64K),所以,程序求解的最大结点数仅限于2000个以内。1995年,FLAC2D已升级为V3.3的版本,其程序能够使用护展内存。因此,大大发护展了计算规模。FLAC3D是一个三维有限差分程序,目前已发展到V3.0版本。 FLAC3D的输入和一般的数值分析程序不同,它可以用交互的方式,从键盘输入各种命令,也可以写成命令(集)文件,类似于批处理,由文件来驱动。因此,采用FLAC程序进行计算,必须了解各种命令关键词的功能,然后,按照计算顺序,将命令按先后,依次排列,形成可以完成一定计算任务的命令文件。 FLAC3D是二维的有限差分程序FLAC2D的护展,能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型,在外力作用下,当材料发生屈服流动后,网格能够相应发生变形和移动(大变形模式)。FLAC3D采用的显式拉格朗日算法和混合-离散分区技术,能够非常准确的模拟材料的塑性破坏和流动。由于无须形成刚度矩阵,因此,基于较小内存空间就能够求解大范围 的三维问题。 三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法,它可以模拟岩土或其他材料的三维力学行为。三维快速拉格朗日分析将计算区域划分为若干四面体单元,每个单元在给定的边界条件下遵循指定的线性或非线性本构关系,如果单元应力使得材料屈服或产生塑性流动,则单元网格可以随着材料的变形而变形,这就是所谓的拉格朗日算法,这种算法非常适合于模拟大变形问题。三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程,并应用了混合单元离散模型,可以准确地模拟材料的屈服、塑性流动、软化直至大变形,尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。 FLAC-3D(Three Dimensional Fast Lagrangian Analysis of Continua)是美国Itasca Consulting Goup lnc开发的三维快速拉格朗日分析程序,该程序能较好地模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学行为,特别适用于分析渐进破坏和失稳以及模拟大变形。它包含10种弹塑性材料本构模型,有静力、动力、蠕变、渗流、温度五种计算模式,各种模式间可以互相藕合,可以模拟多种结构形式,如岩体、土体或其他材料实体,梁、锚元、桩、壳以及人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩、界面单元等,可以模拟复杂的岩土工程或力学问题。 FLAC3D采用ANSI C++语言编写的。 二、FLAC3D的优点与不足 1、怎样查看模型? 答:plot grid 可以查看网格,plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件? 答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2 1、模型建好如何检验所建模型是否有问题 注意:第一步建模一定不要建太大的模型,可以建非常小的模型,确保模型的准确性,这样可以尽快开始后边的步骤,等后边计算分析步骤快速做完验证没有大问题,再将准确的模型建好换上,重新计算得到更加准确的结果。 2、如何检验应力边界、位移边界、加载速率等是否准确 Plot add gpfix:看位移边界 Plo add fap:看应力边界,看加载方向是否准确 Plot add vel:看速率,加载速率是否准确 3、如何查看运算结果 Plo con szz outline on:看z向应力云图,并显示边界网格 Plo add disp:看位移 Plo con zdisp outline on magf 5:看z向位移云图,变形扩大5倍 Plo add bo :添加边界 Plo blo state:看是否屈服 Plo add stensor:查看大主应力方向 Plo con smin:小主应力 Plo con smax:大主应力 Plo con ssi(shear strain increment):查看剪切应变增量云图 Plo bcont smin:查看单元小主应力 4、如何做切面并查看切面结果 在编辑“Edit“菜单中设置切面位置,在上述3命令后加上plane(pl)即可。如:Plo con smin pl:查看切面小主应力 5、后处理 将模型导入到tecplot中做切面,加等值线等,对结果进一步处理,得到更漂亮的一些结果。 其他重要命令总结: 1、如何在圆柱体四周如何施加约束条件? 可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2 快速入门 (GETTING STARTED) 制作:xxxx 2010年12月 2.1 安装启动程序 2.1.1 系统要求 安装运行flac3d 需要的系统最低配臵如下: 处理器:时钟频率至少为1GHZ,处理器的主频越高,那么flac3d 的计算速度将越快。 硬盘:安装软件至少需要12MB 的硬盘空间。如果装载了在线的用户手册,那么还需16MB 的空间。(注意默认情况下,安装软件时会自动装载用户手册)。除此之外,还需要至100MB 的硬盘空间来存储分析计算时生成的各种文件。内存-启动软件至少需要3MB 的内存。在建模过程中,软件所占用的内存,会不断的发生变化(见表2.1)WINDOW 操作系统还限定了软件建模时占用的内存不能超过2GB。 显示器:推荐1024×768 分辨率,16 位彩色显示器。 操作系统:FLAC3D 是32 位操作系统的应用程序,所以基于intel 技术的WINDOWS 98 及以 上操作系统均支持软件的安装和使用。 输出设备:默认情况下,系统图形会输出到系统打印机上。也可以复制到剪贴板上,或者保 存为格式化的文件,这里所说的格式包括:加强型图元文件格式和位图文件(PCX/BMP/JPEG)。用户可以使用set plot 命令来指定输出的形式及格式。 2.1.2 软件的安装 (略) 2.1.3 组件 软件的可执行文件为“F3300.EXE”。FLAC3D 是使用VC++ 7.0 编写的。除了可执行程序外, 还需要两套动态链接库(DLL 文件),一套用来接入和存取各种各样的图形;另一套提供内臵 的各种本构模型。 2.1.4 应用程序和图形处理设备 在使用FLAC3D 时,各种应用软件和图形处理设备会起到很大的辅助作用。 编辑器:任何以ASCII 码为标准格式的文本编辑器都可以用来创建FLAC3D 的数据文件。但是必须要注意一些“先进”的文档编辑器(如WordPerfect, Word等软件),这些编辑器会把格式说明信息编译成标准输出格式,这些说明信息并不能被FLAC3D 识别,所以导入这类文档时会出现 错误。FLAC3D 输入的数据文件必须是标准ASCII 码形式的文件。 图形输出设备:FLAC3D 支持很多种类型的图形处理设备,默认情况下,生成的图形可以“Plot hardcopy”命令来连接到系统默认的打印机以便输出。(或者通过FLAC3D 主窗口中FILE 菜单 栏下的print-view 来设定)“Plot clipboard”命令可以将显示的图形,存放到WINDOWS 剪贴板 上(没有任何文件生成)。该图形接着就可以以加强型图元文件格式被粘贴到其它兼容该格式 的WINDOWS 应用程序中去。“Set plot metafile”命令可以将图形以加强型图元格式存盘,以便 作为计算的参考或日后插入到文档中去。 通过命令:Set plot +关键词(pcx, bitmap, bmp 或者jpg)可以存储为许多图像格式(pcx,bmp,jpeg 等)。输出的这些位图的分辨率由命名行:Set plot FLAC3D快速入门 及简单实例 李佳宇编 LJY指南针教程 前言 FLAC及FLAC3D是由国际著名学者、英国皇家工程院院士、离散元的 发明人Peter Cundall博士在70年代中期开始研究的,主要面对岩土工程的通 用软件系统,目前已经在全球70多个国家得到广泛应用,在岩土工程学术界 和工业界赢得了广泛的赞誉。前国际岩石力学会主席 C.Fairhurst(1994)对 FLAC程序的评价是:“现在它是国际上广泛应用的可靠程序。” 我从研二(2010年)开始接触FLAC3D,最初的原因是导师要求每一个人至 少学会一个数值计算软件,而他嘴里每天念叨最多的就是FLAC,自己当时对数 值计算一无所知,便答应老师要学会FLAC3D。第一次打开软件界面,我心里 就凉了大半截,面对着一个操作界面跟记事本无异的所谓“功能强大”的岩土工 程专业软件,半点兴趣也提不起来。年底,从项目工地回到学校准备论文开题, 老师对我的开题报告非常不满意,当着全教研室师生的面,劈头盖脸大批一顿, 第二天又找谈话。在巨大的压力和强烈的自尊心驱使下,我硬着头皮开始啃 FLAC3D,一个半月之后,终于有了初步的计算结果,对老师有个交代,我也能 回家过年了。 前面这一段过程可能是大多数FLAC3D初学者的必经阶段,或者是即将 开始软件学习的人惧怕的事情。毫无疑问,FLAC3D极其不友好的界面是阻碍 初学者前进的很大障碍,当然还包括它是一个全英文的软件。但是当你费尽周折 的走进FLAC3D的世界,你就会发现它独特的魅力,比如简洁的界面,快捷的 命令流操作,高效的计算方法,不易报错等等。另外一个拿不上台面的优点就是 它非常小巧,包括Manual在内一共才几十兆大小,而且已经被破解成绿色版, 只要把它和命令流装进U盘,你就可以随便找一个身边功能最强大的电脑开始 计算了,如果你有过ANSYS、ABAQUS等大型软件痛苦的安装经历,你便能 毕业之后,本以为不用再接触数值计算,但工作需要使得我又一次开始与理解“绿色版”的含义,当然还请大家尊重知识产权,支持正版。 FLAC3D进行亲密接触。我的领导给了我很多新思路和很大的支持,如今我的 水平比研究生时有了不少提高。于是,我想把我的经验总结成文,希望对初学 者起 1.1常见问题及其解答 Gen separate 不能被识别 答:原因是FLAC3D版本不行,我用3.0的版本不能。 1. FLAC3D是有限元软件吗? 答:不是,是有限差法软件。 2. FLAC3D最先需要掌握的命令有哪些? 答:需要掌握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) 7. 怎样看模型的矢量图? 答:plo dis (xdis, ydis, zdis) 8. 怎样看模型有多少单元、节点? 答:print info 9. 怎样输出模型的后处理图? 答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpg文件。 10. 怎样调用一个文件? 答:使用菜单File/call或者call命令。 11. 如何施加面力? 答:app nstress ran <……> 12. 如何调整视图的大小、角度? 答:综合使用x, y, z, m, Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键。 13. 如何进行边界约束? 答:fix x ran <……>(约束的是速度,在初始情况下约束等效于位移约束) 14. 如何知道每个单元的ID? 答:使用鼠标双击单元的表面,可以知道单元的ID和坐标。 15. 如何进行切片? 答:plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 16. 如何保存计算结果? 答:save filename(文件名可自定义) 17. 如何调用已保存的结果? 答:使用菜单File/call或者命令rest filename(文件名可自定义)。 18. 如何暂停计算? 答:运行中使用Esc命令。 19. 如何在程序中进行暂停,并可恢复计算? 答:在命令中加入pause命令,键入continue命令后可恢复计算。Flac3D中文流体计算
FLAC3D常见命令与使用技巧
FLAC3d后处理步骤
flac3D中文使用手册
FLAC3D原理..
FLAC3D 高手总结
FLAC3D基础知识介绍解析word版本
FLAC3D命令流(挺实用)
FLAC3D常用命令总结
flac3d中文使用指导
FLAC3D快速入门及简单实例
FLAC3D错误提示解决办法