“用tecplot对flac3d计算结果进行后处理”相关讨论帖整理(第二部分)

1、文档来源为：从网络收集整理.word版本可编辑.欢迎下载支持 第二部分：对dynamax大牛的原创主题帖 【原创】用tecplot对flac3d计算 结果进行后处理”及相关讨论帖的整理和总结。 1. dynamax大牛的原创主题帖 【原创】用tecplot对flac3d计算结果进行后 处理”的链接地址为： SimWe 仿真论坛 ？ D01: FLAC3D/FLAC2D ?【原创】用 tecplot 对 flac3d 计算 结果进行后处理 说明：该flac3d-tecplot转换程序，由dynamax 大人利用fish编写，总共有 三个版本，不断改进。最初只能导入位移数据，画出位移等值线，并且不

2、能剔 除挖空单元（null），只能导入全部网格；第二版不再导出null zone，并且可 以指定显示范围；第三版则可进一步同时导出位移和应力结果。 dynamax 发表于 2004-7-11 17:45 前后处理【原创】用tecplot对flac3d计算结果进行后处理 fish写的 因为flac3d中不能标注等值线图 将位移计算结果导入tecplot中进行处理； 可以绘出漂亮的等值线图 熟悉tecplot的可以试一下 不然得到转换后的数据文件还是不知道如何处理 给出了简单的测试实例 283356-tecplot.rar （45.85 KB）（第一版转换程序） dyn amax 切片功能，从菜单

3、中选择，操作为： data = extract = slice from pla ne dyn amax 上图显示成2D Cartesian可能更好 要显示数值的话， plot- Con tour- More- Labels- Show Labels 坐标轴可以关掉 plot- axis- show %-axis dyn amax 我给的例子只能导出全部网格的信息（针对第一版转换程序而言） 如果要只导出开挖后剩下的部分 因为节点编号的关系 要另外做处理才行 否则没有办法做出像下面这样的图 我看代码有点不对，稍加修改应该没有问题，主要是把null部分的网格信息排除 我们努力看看行不行。 dyn

4、amax 根据大家的意见修改了一下 不再导出null zone 可以指定plot range ;Edit the line in the file to set plot range ran n ame tec_ra nge 欢迎试用一 笔323814-flac3d to tecplot.rar (32.53 KB)(第二版转换程序) dyn amax 在Tecplot 10之前，处理的变量只能定义在节点上 在Tecplot 10中，可以指定变量的位置在单元中心 VARLOCATION = CELLCENTERED 同时要求 DATAPACKING=BLOCK 请参考用户手册中的说明 文件头可

5、以这么写 TITLE = FLAC3D to Tecplot 10 VARIABLES = X Y Z DISP XDISP YDISP ZDISP SIG1 SIG2 SIG3 SXX SYY SZZ ZONE T=GLOBAL N=1134, E=318, ZONETYPE=FEBrick DATAPACKING=BLOCK VARLOCATION=(8-13=CELLCENTERED) DT=(SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE ) yan g

6、muya ng 感谢 dynamax : 我还有疑问就是： SIG1 SIG2 SIG3 SXX SYY SZZ 与DISP XDISP YDISP ZDISP 的坐标变量X Y Z好象不相同啊，如何将他们对应起来呢？ 你上面的文件头中写到N=1134， E=318分别指的是节点数与单元数吧。 VARLOCATION=(8-13=CELLCENTERED) 这是什么意思啊？ 8-13是固定的数字么？ dyn amax 注意数据的排列方式必须是block方式 VARLOCATION=(8-13=CELLCENTERED) 是指第8-13个变量定义在单元中心上 单击右边的 zone style 选

7、中要显示的 zone 单击 zone show 选择 show selected only 单击 con tour type 选择 lines 这样就可以绘出等值线图 dyn amax 重写了一下，欢迎试用！ 可以导出位移和应力等信息 位移等信息定义在节点上 tecplot 和flac3d绘出的图完全一致 应力等信息定义在单元中心上 因为插值方法的不同 tecplot绘出的图和在flac3d中用average方式绘制的图比较类似 但并不完全一致，尤其在边界上 使用的时候要注意这一点 325808-tecplot.rar (86.14 KB)(第三版转换程序) dyn amax 先做一个slic

8、e 再显示slice上的con tour dyn amax 单击右边的 zone style 选中要显示的 slice 单击 zone show 选择 show selected only 单击 con tour type 选择 lines 这样就可以绘出等值线图 plot = con tour = labels = show labels 可显示等值线上面的数值 dyn amax 目前最大的问题是不能做应力矢量图 幸好FLAC3D本身做得还不错 没有必要再折腾了 ruigrate 个人认为用surf做后处理，图形比较清晰好看。下面是用surf后处理得图形, 可以比较一下。 dyn amax

9、我没有用过surfer，想请教一下 用surfer能处理三维有限元数据的结果吗 ruigrate surfer也只能处理一个截面，对截面绘制等值线图和位移三维图，上面的图就是 我用surfer处理的FLAC3D结果，是复合地基的桩土界面变形三维图。 dyn amax 的例子拓展了一个后处理的空间，可以自己根据需要进行修改和扩 充，比如我已经进行了孔压、超孔压、孔压比的后处理代码，其它的也可以参考， 但是这个原创性的工作不能否定。 jonavy 我做的图老是显示周围的图框，我怎么才能不显示图框只显示这个切片那，高人盼 你指点一下谢谢! jonavy 解决了 !呵呵,从双击的属性里有个grid那里

10、改改属性就可以了 haita ng 今天由于出图需要，再把dynamax大哥的帖子瞻仰了一遍，又受到很大启发。 以前剖分的单元主要是六面体和四面体，所以问题不大，现在我的模型很多是三 棱柱。 wedge，而且有的是pyrami。好像导入后就出现了一些问题，都变成了四面体。 早上研究了一下，按照datapacking = point的格式已经实现了各种单元的导入， 位移可以输出了，可是应力是体现在单元中心的。 我就想让对单元进行平均，分 配到节点，但还没实现。看了帖子后，这才发现datapacking = block是个好办法。 由于没有tecplot的用户手册，很多代码不了解意思，dynama

11、x大哥的ftp 有一直上不去，所以请dynamax大哥发一份tecplot的用户手册到我邮箱.谢谢 啦。鉉 Bob to楼主如何只导入GR0UP=3的网格和位移应力啊 Bob 呵呵，用 RAN NAME TEC_RANGE RAN GRO 3 A GRO 4 A 就行了 benjackxu ABAQUS这个软件和tecplot不是同一类型的软件，无法比较。ABAQUS为有限 元软件，tecplot是图像处理软件，至于大家为什么喜欢用tecplot画等值线， 想是习惯问题，同时也更直观一些。 shufflyy 谁能把楼主的那个压缩包给解释一下怎么用，非常的急用，千里之外，万分感 谢！ ！ scl

12、yyaoyao 压缩包里的是flac程序的dat文件，计算完成后在flac中call flac3d2 tecplot .dat， 自动会在该目录下生成tec10.dat，在tecplot中打开即可。 zkp-08 FLAC3D计算结果导入TECPLOT进行后处理的步骤： 第一步：用FLAC3D打开已经完成的*.sav文件 第二步：运行 dynamax 无私共享的 flac3d2Tecplot.dat 文件，即 call flac3d2Tecplot.dat ， 根据你所用电脑的计算速度， 可能要等待一段时间来写入数据， 运行结束时会在命令提上行 里显示写入成功的提示，在命令文件所在的文件夹里会

13、发先tec10.dat数据文件(可以根据 自己的需要更改文件名)。 第三步：打开tecplot 程序，执行【File】/【Load Data Files(s)】命令打开tec10.dat数据 文件，即可进行编辑了。 再次感谢dynamax的无私奉献精神 2. mathing 对dynamax的FLAC3D2Tecplot程序的修改(解决部分 “ ”误) 帖子的链接地址为： mathi ng dynamax的FLAC3D2Tecplot程序的修改(解决部分 “ ”误) dynamax的FLAC3D2Tecplot程序在下面情况下因为改变了原来的结点单 元编号会出现“ 错昔误： (1) 使用del

14、删除单元，但用model null不会出错 (2) 使用了 gen separate 命令 解决的方法是将结点重新编号，我尝试了一下，将dynamax的程序进行了 一些修改，以适用于上面两种情况，也适用于有接触面的模型。修改后的程 序见附件，主要修改处： 1. 将生成的tecplot文件的数据结构由原来的BLOCK形式改为POINT形式 2. 去掉应力导出，只能到处位移 3. 结点重新编号 但由于在 Write Zone Connectivity这一部分中加入了一个查找函数，有大量的 数据查找，因此如果单元数很多(上万)，导出的速度会很慢。 使用方法： 1. 如果模型中没有del命令或gen

15、separate，请使用dynamax的程序， 速度快 2. 如果有del命令或gen separate命令，但单元数比较少(几千个)，可 以使用附件中的程序 3. 如果有del命令或gen separate命令，但单元数比较多(几万个)，最 好 不要全部导出，因为速度会非常非常慢，可以只导出自己关心的一小部分， 例如，模型范围是z= 0100 (几万个单元)，如果只导处01这一部分(几 千个单元)，具体步骤如下： rest aaa.sav del ran z 1 100 cal Flac3d2Tecplot_A.dat 注意要使用del命令，千万不要用 model null，并且要在cal

16、Flac3d2Tecplot_A.dat 之前。 需要改进的地方：在write_zone中每个节点都要进行链表遍历查找 (find_tec_num )，速度太慢 蹇Flac3d2Tecplot A.rar (1.33 KB) 3. haitang 基于dynamax修改后的flac3d tecplot命令流源文件，链接地址为： ;01/12/2007 ;Modificati ons: (1) 将应力输出项增加为(SIG1,SIG2,SIG3,SXX,SYY,SZZ,SXYSYZ,SZX) (2) 可控制输出的GROUP,若不想输出某个组，请ran name tec_range后面 加入组名 F

17、Iac3D Tecplot average by zon e.rar (2.58 KB) Haita ng 的修订 ;orig ini by dyn amaxsimwe, modified by Haita ngHohai ;01/12/2007 ;Modificati ons: ；(1)将应力输出项增加为 ;(SIG1,SIG2,SIG3,SXX,SYY,SZZ,SX YSYZ,SZX) ；(2)可控制输出的GROUP,若不想输出某个组， ;请ran name tec_range 后面加入组名 J J ;18/12/2006 ;Modificati ons: ;(1) adapted to all types of zones ;(brick,pyramid,wedge,tetra ) ;(2) If all the zones are tetra, use the format ;ZONETYPE=FETETRAHEDRON dyn amax 的修订 ；9/14/2004(第三版) ;Modificati ons: ; cha nge data pack ing mode from point to block ;(2) add zon e-related info to be plotted J J ;9/11/2004(第二版) ;Modifica