dynamax的FLAC3D2Tecplot程序的修改（解决部分“@”错误）

mathing

dynamax的FLAC3D2Tecplot程序在下面情况下因为改变了原来的结点单元编号会出现“@”错误:
（1）使用del删除单元，但用model  null不会出错
（2）使用了gen separate命令
解决的方法是将结点重新编号，我尝试了一下，将dynamax的程序进行了一些修改，以适用于上面两种情况，也适用于有接触面的模型。修改后的程序见附件，主要修改处：
  1.将生成的tecplot文件的数据结构由原来的BLOCK形式改为POINT形式
  2.去掉应力导出，只能到处位移
  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），速度太慢

[ 本帖最后由 mathing 于 2007-9-24 19:30 编辑 ]

mathing

原来的那个例子，但使用del命令

new
gen zone brick size 6 1 20 p1=3.0,0.0,0.0 ratio 0.9 1.0 0.97
gen zone brick size 20 1 20 p0=3.0,0.0,0.0 p1=20.0,0.0,0.0 &
    ratio 1.08 1.0 0.97
ini z mul 0.5
model mohr
prop bul 2.e8 shea 1.e8 cohesion 1.e5
prop friction 0. dilation 0. tension 1.e10
fix x range x -.1 .1
fix x y z range z -.1 .1
fix x y z range x 19.9 20.1
fix y
fix x y z range x -.1 3.1 z 9.9 10.1
free x range x 2.9 3.1 z 9.9 10.1
ini zvel -0.5e-5 range x -.1 3.1 z 9.9 10.1
del ran z 6 9 x 2 5
cyc 300

使用原来的程序，会出现“@”错误，但使用修改后的程序，导出的结果以及flac3d计算的结果如下

男人

转入比较慢，希望有新的突破！

wangyi2002598

顶啊

lovemm

xiexie,正在学习中!

zql2001

xi e le

arnoldray

谢谢

humble

谢谢了，不错！

xiarilawage

谢谢了，下来看看

wangminsuning

mathing你好，使用了你在dynamax基础上修改的程序，我现在的模型，用了del命令，也有接触面单元，但是想得到应力等值线图，但你的程序只可以导入位移结果，能否把应力结果也导入？期待回复，谢谢！

lookcity

如果想导出应力的话，由于应力是以单元平均值为基础的，因此必须选择使用BLOCK形式的数据格式。

tao787

高手，希望可以处理fluent的结果

276774950

应力用得比较多啊~~迫切期待导出应力程序

cwww_com

哇！强人很多！

rainman1v1

应力用得比较多啊~~迫切期待导出应力程序

dutt

其实也不必一定要导入TECPLOT，FLAC3D自带的plot 和print 命令在很多情况下是可以满足后处理要求的，配合hist会更好的。

okpzy

多谢，有点慢

disciple_

很好，很强大，学习了~！！@~！

岩土小弟

下下来学习学习

xinshe_cgs

谢谢，很好！