!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/post1
file,Jobname,rst !指明从哪一个结果文件中读取数据
*get,NodeNum,NODE,0,COUNT !得到模型的所有节点数
*get,nd,NODE,0,NUM,MIN !得到模型的最小节点编号
*dim,nodes,array,NodeNum !定义一个存储节点的数组
*dim,xyz,array,NodeNum,3 !定义一个存储节点坐标的数组
*dim,NodeTemp,array,NodeNum !定义一个存储节点的数组
set,4,1 !读入第四载荷步、第一个子步的结果
*do,i,1,NodeNum,1
nodes(i)=nd !将节点编号存储在nodes数组中
xyz(i,1)=NX(nd) !在xyz(i,1)中存储节点的X坐标
xyz(i,2)=NY(nd) !在xyz(i,1)中存储节点的Y坐标
xyz(i,3)=NZ(nd) !在xyz(i,1)中存储节点的Z坐标
*Get,NodeTemp(i),NODE,nodes(i),TEMP
nd=NDNEXT(nd) !取得下一个节点编号
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*get,ElementNum,ELEM,0,COUNT !得到模型的所有单元数
*get,el,ELEM,0,NUM,MIN !得到模型的最小单元编号
*dim,elems,array,ElementNum !定义一个存储单元的数组
*dim,ndlst,array,ElementNum,8 !定义一个存储节点坐标的数组
*do,i,1,ElementNum,1
elems(i)=el !将单元编号存储在elems数组中
*do,j,1,8,1
ndlst(i,j)=NELEM(el,j) !将单元对应的节点编号存储在ndlst数组中
*enddo
el=ELNEXT(el) !取得下一个单元编号
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*cfopen,mesh,dat !输出的文件名为:mesh.dat
!写TecPlot文件的文件头
*vwrite
(’TITLE=Ansys Temperature Analysis’)
!写变量名
*vwrite
(’VARIABLES=”X”,”Y”,”Z”,”Temp”‘)
!写节点数和单元数
*vwrite,NodeNum,ElementNum
(’ZONE n=’f6.1,’ e=’f6.1,’ f=fepoint’,' et=brick’)
!写节点的坐标和温度值
*vwrite,xyz(1,1),xyz(1,2),xyz(1,3),NodeTemp(1)
(4f12.6)
!写一个空行
*vwrite
(’ ‘)
!写各单元对应的节点号
*vwrite,ndlst(1,1),ndlst(1,2),ndlst(1,3),ndlst(1,4),ndlst(1,5),ndlst(1,6),ndlst(1,7),ndlst(1,8)
(8f6.0)
!关闭文件
*CFCLOSE
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!