haiqsim 发表于 2009-5-12 01:04:53

FLAC3D计算结果—Tecplot画位移等值线详细操作步骤介绍(前辈帖子总结)

本帖最后由 haiqsim 于 2009-5-17 00:34 编辑

主要内容:
(1) 详细介绍在FLAC3D中计算完毕后,把计算结果导入Tecplot中,绘制位移等值线的操作步骤。
(2) 对dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”及相关讨论帖的整理和总结。

说明:
这个东东是小弟今日来利用Tecplot进行Flac3d后处理绘图的学习小总结,当然主要是参考各位前辈的讲解和讨论,自以为小有所得,特来分享之:)


(1) 本详细操作步骤主要参照dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”中的讲解和讨论内容,和pan801010的帖子“Tecplot 画位移等值线,比较详细的介绍!!!”
链接地址为:
http://forum.simwe.com/viewthread.php?tid=283356&highlight=tecplot

(2) 文中绘图使用的Tecplot程序,为Lakewater版主PJ的Tecplot10.3版本
链接地址为:
http://forum.simwe.com/viewthread.php?tid=306811&highlight=tecplot

(3) 文中使用的FLAC3D到Tecplot10的数据转换程序为dynamax大牛编写的Flac3d2Tecplot.dat程序的第三版,可以实现的功能为:可以导出位移和应力等信息;不再导出null zone,可以指定输出的范围(组),plot range。
dynamax编制的Flac3d2Tecplot.dat程序(第三版)的下载链接为:
http://forum.simwe.com/viewthread.php?tid=283356&extra=&highlight=tecplot&page=3 页面内的53#楼,名为325808-tecplot.rar


今天总算传完了:victory:

Tecplot画等值线的操作步骤详细介绍(第一部分),在第二楼和第四楼,前后共20个附件,大小不到2M。(有人说附件缺少第11部分,所以在13#楼重新发了part11)

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

第五、第六楼为相关附件的文字版内容;附件中除此之外有相关的操作截图。

haiqsim 发表于 2009-5-12 01:16:45

本帖最后由 haiqsim 于 2009-5-12 23:57 编辑

文件中图片比较多,所以文件较大有2M多,一个附件只能100K,所以又这么多附件:( ;假如直接帖子页面上的话有很多图片和链接,比较麻烦啦,所以搞成附件,反正不要钱,有兴趣的话就慢慢下喽。
这是Tecplot画等值线的操作步骤详细介绍(第一部分)。


晕!好像只有11个,总共有20个啊,上传附件大小有限制,郁闷:Q

haiqsim 发表于 2009-5-12 01:24:44

本帖最后由 haiqsim 于 2009-5-14 00:05 编辑

下面是第二部分:“用tecplot对flac3d计算结果进行后处理”相关讨论帖整理(第二部分)
其中包括:
       dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”和后续讨论内容;
       mathing对dynamax的FLAC3D2Tecplot程序的修改(解决部分“@”错误);
   haitang基于dynamax修改后的flac3d-tecplot命令流源文件链接地址



我晕啊,还有一日内上传文件的大小限制!:funk: :'( ,只好明天再传:(

haiqsim 发表于 2009-5-12 23:59:23

本帖最后由 haiqsim 于 2009-5-14 00:04 编辑

3# haiqsim

Tecplot画等值线的操作步骤详细介绍(第一部分),分卷12-20

又浪费了上传文件的流量,晕

haiqsim 发表于 2009-5-13 00:02:33

本帖最后由 haiqsim 于 2009-5-13 00:24 编辑

由于有附件大小限制,加上操作不熟,文件总是传不上来,只好先发文字了,以后有机会再传文件和图片。

第一部分:Tecplot画位移等值线的具体操作步骤。

1.运行FLAC3D程序,先调用FLAC3D计算结果的保存文件xxxx.sav,命令为:restorexxxx.sav
   然后,调用dynamax编制的flac3d到tecplot的转换程序Flac3d2Tecplot.dat,命令为:callFlac3d2Tecplot.dat
   则在flac3d计算结果文件的当前目录下,生成一个名为tec10.dat文件,这个文件就是将flac3d计算结果转换成符合tecplot 10的数据格式后的文件,可供tecplot 10调用。

2.在tecplot 10中,调用第一步生成的tec10.dat文件,将FLAC3D计算结果导入到tecplot 中,具体操作为:
   运行tecplot 10程序,点击File主菜单,在下拉菜单中选择Load Data File(s)选项,

   在弹出的Load Data File对话框中,选中上一步生成的tec10.dat文件,点击“打开”,

   之后一路点击默认选项,得到显示FLAC3D模型网格和坐标面的图形。

3.   使用切片工具(Iinteractive slicing tool),显示模型内部的一个指定剖面。(相当于在flac3d中设置cutting plane后,图形显示只选择on plane)

4.   点击Iinteractive slicing tool 之后,再点击下面的tool details,
      弹出如下对话框,显示切片工具的详细设置,
      先勾选Show Slices选项(显示切片);
      再选择切片的方向(X-Planes和Y-Planes等分别表示法线方向平行于X轴和Y轴等方向的切片);
      然后指定切片的位置,在Position后面输入数字,该数字为沿着所选切片的法线方向上的切片位置坐标。
      切片设置好后,点击Close,关闭设置对话框。

5.   将切片截面上的数据读出,形成除整体zone外的另一格zone,操作如下:
      选择主菜单 【Data】/ 【Extract】/ 【Current Slice】
      之后弹出如下对话框,点击Extract,
      弹出“Slice extraction successful”信息对话框,显示切片上数据提取成功。

6.   点击zone style ,弹出Zone Style设置对话框,

7.   点击Zone Style下的Contour选项,发现有两个zone,在Zone Show中关掉1号zone,GLOBAL;
      激活2号zone,Slice:Y=0.5;

8.   点击contour下的contour type,选取lines(只显示等值线)

9.   取消Mesh选项前面的勾选(不显示网格),激活Contour选项(只显示等值线),

10. 点击Contour后面的的详细设置选项卡,弹出Contour Details对话框,选择要显示的等值线变量(位移或应力,此处选位移DISP),

11. 点击等值线变量后的More选项,显示完全的Contour Details对话框;
      再点击Labels选项卡,勾选Show Labels选项(这样使得显示等值线的同时,显示等值线的数值);

12. 此时得到如下图所示的等值线图(含有坐标面网格和坐标轴),

13. 上图中背景为坐标面网格和坐标轴,若要去掉这些,需在显示屏中点击选中坐标面后,点击Object Details,弹出如下Axis Details对话框,

14. 将Axis Details对话框中,所有选项卡内的勾选都取消,即可得到干净的等值线图,目的达成。

haiqsim 发表于 2009-5-13 00:11:13

本帖最后由 haiqsim 于 2009-5-13 00:27 编辑

第二部分:对dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”及相关讨论帖的整理和总结。

1. dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”的链接地址为:
http://forum.simwe.com/viewthread.php?tid=283356&highlight=tecplot

SimWe仿真论坛 » D01:FLAC3D/FLAC2D » 【原创】用tecplot对flac3d计算结果进行后处理
说明:该flac3d-tecplot转换程序,由dynamax大人利用fish编写,总共有三个版本,不断改进。最初只能导入位移数据,画出位移等值线,并且不能剔除挖空单元(null),只能导入全部网格;第二版不再导出null zone,并且可以指定显示范围;第三版则可进一步同时导出位移和应力结果。

dynamax发表于 2004-7-11 17:45
[前后处理] 【原创】用tecplot对flac3d计算结果进行后处理

fish写的
因为flac3d中不能标注等值线图
将位移计算结果导入tecplot中进行处理 ;
可以绘出漂亮的等值线图
熟悉tecplot的可以试一下
不然得到转换后的数据文件还是不知道如何处理
给出了简单的测试实例
283356-tecplot.rar (45.85 KB) (第一版转换程序)



dynamax
切片功能,从菜单中选择,操作为:
data = >
extract = >
slice from plane

dynamax
上图显示成2D Cartesian可能更好
要显示数值的话,
plot-> Contour-> More-> Labels-> Show Labels
坐标轴可以关掉
plot-> axis-> show %-axis

dynamax
我给的例子只能导出全部网格的信息 (针对第一版转换程序而言)
如果要只导出开挖后剩下的部分
因为节点编号的关系
要另外做处理才行
否则没有办法做出像下面这样的图


我看代码有点不对,稍加修改应该没有问题,主要是把null部分的网格信息排除我们努力看看行不行。

dynamax
根据大家的意见修改了一下
不再导出null zone
可以指定plot range
; Edit the line in the file to set plot range
ran name tec_range
欢迎试用

323814-flac3d_to_tecplot.rar
(32.53 KB)
(第二版转换程序)

dynamax
在Tecplot 10之前,处理的变量只能定义在节点上

在Tecplot 10中,可以指定变量的位置在单元中心
VARLOCATION = CELLCENTERED
同时要求
DATAPACKING=BLOCK
请参考用户手册中的说明
文件头可以这么写
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=(=CELLCENTERED)

DT=(SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE )

yangmuyang
感谢dynamax:
我还有疑问就是:
"SIG1" "SIG2" "SIG3" "SXX" "SYY" "SZZ" 与 "DISP" "XDISP" "YDISP" "ZDISP"的坐标变量"X" "Y" "Z"好象不相同啊,如何将他们对应起来呢?
你上面的文件头中写到N=1134,E=318分别指的是节点数与单元数吧。
VARLOCATION=(=CELLCENTERED)
这是什么意思啊?
是固定的数字么?

dynamax
注意数据的排列方式必须是block方式
VARLOCATION=(=CELLCENTERED)
是指第8-13个变量定义在单元中心上

单击右边的 zone style
选中要显示的 zone
单击 zone show 选择 show selected only
单击 contour type 选择 lines
这样就可以绘出等值线图

dynamax
重写了一下,欢迎试用!
可以导出位移和应力等信息
位移等信息定义在节点上
tecplot和flac3d绘出的图完全一致
应力等信息定义在单元中心上
因为插值方法的不同
tecplot绘出的图和在flac3d中用average方式绘制的图比较类似
但并不完全一致,尤其在边界上`
使用的时候要注意这一点

325808-tecplot.rar (86.14 KB) (第三版转换程序)

dynamax
先做一个slice
再显示slice上的contour

dynamax
单击右边的 zone style
选中要显示的 slice单击 zone show 选择 show selected only
单击 contour type 选择 lines
这样就可以绘出等值线图
plot = > contour = > labels = > show labels
可显示等值线上面的数值

dynamax
目前最大的问题是不能做应力矢量图
幸好FLAC3D本身做得还不错
没有必要再折腾了

ruigrate
个人认为用surf做后处理,图形比较清晰好看。下面是用surf后处理得图形,可以比较一下。

dynamax
我没有用过surfer,想请教一下
用surfer能处理三维有限元数据的结果吗

ruigrate
surfer也只能处理一个截面,对截面绘制等值线图和位移三维图,上面的图就是我用surfer处理的FLAC3D结果,是复合地基的桩土界面变形三维图。

lakewater
dynamax的例子拓展了一个后处理的空间,可以自己根据需要进行修改和扩充,比如我已经进行了孔压、超孔压、孔压比的后处理代码,其它的也可以参考,但是这个原创性的工作不能否定。

jonavy
我做的图老是显示周围的图框,我怎么才能不显示图框只显示这个切片那,高人盼你指点一下谢谢!2


jonavy
解决了!!呵呵,从双击的属性里有个grid那里改改属性就可以了

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

Bob
to 楼主
如何只导入GROUP=3的网格和位移
应力啊
Bob
呵呵,用RAN NAME TEC_RANGERAN GRO 3 A GRO 4 A就行了

benjackxu
ABAQUS这个软件和tecplot不是同一类型的软件,无法比较。ABAQUS为有限元软件,tecplot是图像处理软件,至于大家为什么喜欢用tecplot画等值线,想是习惯问题,同时也更直观一些。

shufflyy
谁能把楼主的那个压缩包给解释一下怎么用,非常的急用,千里之外,万分感谢!!!

sclyyaoyao
压缩包里的是flac程序的dat文件,计算完成后在flac中call flac3d2tecplot.dat,自动会在该目录下生成tec10.dat,在tecplot中打开即可。

zkp-08
FLAC3D 计算结果导入TECPLOT进行后处理的步骤:
第一步:用FLAC3D打开已经完成的*.sav文件2 B5 T+ j' R" P; r. K+ m7 c
第二步:运行dynamax无私共享的flac3d2Tecplot.dat文件,即call flac3d2Tecplot.dat,根据你所用电脑的计算速度,可能要等待一段时间来写入数据,运行结束时会在命令提上行里显示写入成功的提示,在命令文件所在的文件夹里会发先tec10.dat数据文件(可以根据自己的需要更改文件名)。5 _) f7 k5 [+ [( S* g
第三步:打开tecplot程序,执行【File】/【Load Data Files(s)】命令打开tec10.dat数据文件,即可进行编辑了。!
再次感谢dynamax的无私奉献精神


2. mathing 对dynamax的FLAC3D2Tecplot程序的修改(解决部分“@”错误)
帖子的链接地址为:
http://forum.simwe.com/viewthread.php?tid=799634&highlight=tecplot

mathing
dynamax的FLAC3D2Tecplot程序的修改(解决部分“@”错误)

dynamax的FLAC3D2Tecplot程序在下面情况下因为改变了原来的结点单元编号会出现“@”错误:
   (1)使用del删除单元,但用model null不会出错
   (2)使用了gen separate命令/@
       解决的方法是将结点重新编号,我尝试了一下,将dynamax的程序进行了一些修改,以适用于上面两种情况,也适用于有接触面的模型。修改后的程序见附件,主要修改处:

1.将生成的tecplot文件的数据结构由原来的BLOCK形式改为POINT形式
2.去掉应力导出,只能到处位移J
3.结点重新编号
    但由于在Write Zone Connectivity这一部分中加入了一个查找函数,有大量的数据查找,因此如果单元数很多(上万),导出的速度会很慢。

使用方法:
1.如果模型中没有del 命令或gen separate,请使用dynamax的程序,速度快
2.如果有del命令或gen separate命令,但单元数比较少(几千个),可以使用附件中的程序
3.如果有del命令或gen separate命令,但单元数比较多(几万个),最好
不要全部导出,因为速度会非常非常慢,可以只导出自己关心的一小部分,例如,模型范围是z=0~100(几万个单元),如果只导处0~1这一部分(几千个单元),具体步骤如下:

rest aaa.sav
del ran z 1 100
cal Flac3d2Tecplot_A.dat
注意要使用del 命令,千万不要用model null,并且要在cal Flac3d2Tecplot_A.dat之前。
需要改进的地方:在write_zone中每个节点都要进行链表遍历查找(find_tec_num),速度太慢
Flac3d2Tecplot_A.rar
(1.33 KB)


3. haitang基于dynamax修改后的flac3d-tecplot命令流源文件,链接地址为:
http://forum.simwe.com/viewthread.php?tid=724033&page=1&authorid=54640

;; 01/12/2007
;; Modifications:

(1) 将应力输出项增加为(SIG1,SIG2,SIG3,SXX,SYY,SZZ,SXY,SYZ,SZX)
(2) 可控制输出的GROUP,若不想输出某个组,请ran name tec_range后面加入组名
Flac3D_Tecplot_average_by_zone.rar
(2.58 KB)

lmm 发表于 2009-5-13 07:42:33

楼主真是个有心人,赞一下!!!!!!!!!!

lookcity 发表于 2009-5-13 09:19:50

支持楼主   
汇总的不错

haiqsim 发表于 2009-5-14 00:03:23

7# lmm

兄台的硕士学位论文和发在岩土工程学报上的相关文章我都拜读过了,有太多值得我辈学习的地方!
没有什么高深理论上的所谓心得,于是只好干点体力活来试图做点小工作,来博大家一笑了:lol

hua320king 发表于 2009-5-14 10:03:32

感谢楼主的分享,这正是我想要的啊。非常感谢!!

hua320king 发表于 2009-5-15 16:28:32

可是怎么没有第11部分,所以解压缩不了。

a270134858 发表于 2009-5-16 17:33:48

谢谢你的无私奉献了

haiqsim 发表于 2009-5-17 00:27:05

可是怎么没有第11部分,所以解压缩不了。
hua320king 发表于 2009-5-15 16:28 http://forum.simwe.com/images/common/back.gif

2#楼不是11个附件吗?
那我再把part11重新发下好了

xinanfeng 发表于 2009-5-19 16:55:52

不错支持

绍兴人 发表于 2009-5-19 22:15:58

恩,,楼住真好!!!

yangshikou 发表于 2009-6-7 11:41:35

不错支持

yfldar 发表于 2009-6-9 14:08:11

6# haiqsim 楼主真是个有心人啊,赞一个!

baihe111 发表于 2009-6-12 10:51:27

谢谢楼主啦!

baihe111 发表于 2009-6-12 10:51:49

分享一下!!!

baihe111 发表于 2009-6-12 10:52:02

哈哈呵呵呵
页: [1] 2 3 4 5 6 7
查看完整版本: FLAC3D计算结果—Tecplot画位移等值线详细操作步骤介绍(前辈帖子总结)