凝固计算后，提取厚度，总是写入文件不成功
我在凝固计算后，提取厚度，想写入文件，总是不成功，但不知道错误在哪，谁能帮我看一下，谢谢了这样的一个1/4圆的一个二维非稳态凝固。需要得到如图所示的那条红线随时间的温度变化，那条红线上的温度小于某一值即表示该处为凝固坯壳开始位置，获得该点的位置信息，写入文件。我写的如下：/post1
/uis,msgpop,3
*set,x1,0
*set,x2,0.19
*set,y1,0
*set,y2,0
*set,xdivd,190+1
*dim,x,array,86 !定义一个临时存储节点坐标2的数组
*dim,cout,array,xdivd
*dim,tem,array,86 !定义一个临时存储节点温度2的数组
*dim,xy,array,xdivd,2 !定义一个存储节点坐标1的数组
*dim,NodeTemp,array,xdivd !定义一个存储节点温度1的数组
*dim,xx,array,86 !转换为mm单位存储节点X坐标
*dim,ttem,array,86
*cfopen,huaigang,dat !输出的文件名为：huaigang.dat
!写TecPlot文件的文件头
*vwrite
('TITLE=counter')
*vwrite
('VARIABLES= "Time","X"')
*do,n,1,86,1!整个计算时间
set,,,,,1
path,1,2,,xdivd-1 !定义路径
ppath,1,,x1,y2
ppath,2,,x2,y2
pdef,temp,temp
paget,pathdat,table
*do,k,1,xdivd,1
xy(k,1)=pathdat(k,4,1) !在xy(k,1)中存储点的X坐标
xy(k,2)=y2 !在xy(k,2)中存储点的Y坐标
NodeTemp(k)=pathdat(k,5,1) !在NodeTemp(k)中存储点的温度
*enddo
padele,all
pathdat(1,1,1)=
*do,k,1,xdivd,1
*if,NodeTemp(k),lt,1420,and,NodeTemp(k),gt,0,then !获取有效节点
tem(n)=NodeTemp(k)
x(n)=xy(k,1)
*endif
*enddo
*vwrite,xdivd
('ZONE I='f6.1', f=point')
xx(n)=x(n)*1000
ttem(n)=tem(n)
*vwrite,time,xx(n)
(3f12.6)
xy(1,1)=
NodeTemp(1)=
tem(1)=
x(1)=
*CFCLOSE
*enddo
这之中有做过调整，但总是不能成功写入，谁能帮我看下，谢谢了
请问楼主你解决了没有，最近我也再用数据提取，也是总不成功，希望能从你这学习 lansejixian 发表于 2013-5-13 16:55 static/image/common/back.gif
请问楼主你解决了没有，最近我也再用数据提取，也是总不成功，希望能从你这学习 ...
还没有，不知道什么原因
页:
[1]