dynamax的FLAC3D2Tecplot程序的修改（解决部分“@”错误）
本帖最后由 cloud328 于 2010-5-20 10:28 编辑dynamax的FLAC3D2Tecplot程序在下面情况下因为改变了原来的结点单元编号会出现“@”错误:
（1）使用del删除单元，但用modelnull不会出错
（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 编辑 ] 原来的那个例子，但使用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计算的结果如下 转入比较慢，希望有新的突破！ 顶啊
:@ xiexie,正在学习中! xi e le :) 谢谢 谢谢了，不错！ 谢谢了，下来看看 mathing你好，使用了你在dynamax基础上修改的程序，我现在的模型，用了del命令，也有接触面单元，但是想得到应力等值线图，但你的程序只可以导入位移结果，能否把应力结果也导入？期待回复，谢谢！ 如果想导出应力的话，由于应力是以单元平均值为基础的，因此必须选择使用BLOCK形式的数据格式。 高手，希望可以处理fluent的结果 应力用得比较多啊~~迫切期待导出应力程序 哇！强人很多！ 应力用得比较多啊~~迫切期待导出应力程序 其实也不必一定要导入TECPLOT，FLAC3D自带的plot 和print 命令在很多情况下是可以满足后处理要求的，配合hist会更好的。 多谢，有点慢 很好，很强大，学习了~！！@~！ 下下来学习学习 谢谢，很好！
页:
[1]
2