如何能把结果与节点坐标显示在一个文件中

a190189563

  我想把结果与节点的坐标显示在一个文件中，例如我要得到综合位移的结果，但是在的到的结果列表中，没有节点的坐标，我想把结果和坐标显示在一个文件中，因为我要做表格。谢谢

Rainyboy

查询数据的话用*GET命令，写文件的话用*VWRITE，实际上都是“RTM”能解决的问题。

bbmingming

回复 2 # Rainyboy 的帖子

能不能具体点，讲结果和坐标同时显示出来。

Rainyboy

回复 3 # bbmingming 的帖子

恩，找不到以前用过的代码了，本来想写个例子，临时有点事……所以得晚一些时候再弄了，不过翻翻手册实现这个是不难的。

Chelsea

参考下面的，下面的命令流是导出坐标温度的，自己修改以下

1. ********************************输出结果命令流********************************
   
2. /post1
   
3. file,Jobname,rst                       !指明从哪一个结果文件中读取数据
   
4. *get,NodeNum,NODE,0,COUNT                    !得到模型的所有节点数
   
5. *get,nd,NODE,0,NUM,MIN                       !得到模型的最小节点编号
   
6. *dim,nodes,array,NodeNum                     !定义一个存储节点的数组
   
7. *dim,xyz,array,NodeNum,3                     !定义一个存储节点坐标的数组
   
8. *dim,NodeTemp,array,NodeNum                  !定义一个存储节点的数组
   
9. set,4,1                                      !读入第四载荷步、第一个子步的结果
   
10. *do,i,1,NodeNum,1
    
11.     nodes(i)=nd                              !将节点编号存储在nodes数组中
    
12.     xyz(i,1)=NX(nd)                          !在xyz(i,1)中存储节点的X坐标
    
13.     xyz(i,2)=NY(nd)                          !在xyz(i,1)中存储节点的Y坐标
    
14.     xyz(i,3)=NZ(nd)                          !在xyz(i,1)中存储节点的Z坐标
    
15.     *Get,NodeTemp(i),NODE,nodes(i),TEMP
    
16.     nd=NDNEXT(nd)                            !取得下一个节点编号
    
17. *enddo
    
18. *get,ElementNum,ELEM,0,COUNT                  !得到模型的所有单元数
    
19. *get,el,ELEM,0,NUM,MIN                        !得到模型的最小单元编号
    
20. *dim,elems,array,ElementNum                   !定义一个存储单元的数组   
    
21. *dim,ndlst,array,ElementNum,8                 !定义一个存储节点坐标的数组
    
22. *do,i,1,ElementNum,1
    
23.     elems(i)=el                               !将单元编号存储在elems数组中
    
24.     *do,j,1,8,1
    
25.         ndlst(i,j)=NELEM(el,j)              !将单元对应的节点编号存储在ndlst数组中
    
26.     *enddo
    
27.     el=ELNEXT(el)                             !取得下一个单元编号
    
28. *enddo
    
29. *cfopen,mesh,dat                              !输出的文件名为：mesh.dat
    
30. !写TecPlot文件的文件头
    
31. *vwrite
    
32. ("TITLE=Ansys Temperature Analysis")
    
33. !写变量名
    
34. *vwrite
    
35. ("VARIABLES="X","Y","Z","Temp"")
    
36. !写节点数和单元数
    
37. *vwrite,NodeNum,ElementNum
    
38. ("ZONE n="f6.1,"   e="f6.1,"    f=fepoint","    et=brick")
    
39. !写节点的坐标和温度值
    
40. *vwrite,xyz(1,1),xyz(1,2),xyz(1,3),NodeTemp(1)
    
41. (4f12.6)
    
42. 
    
43. !写一个空行
    
44. *vwrite
    
45. (" ")
    
46. !写各单元对应的节点号
    
47. *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)
    
48. (8f6.0)
    
49. !关闭文件
    
50. *CFCLOSE
    
51. finish
    
52. ********************************输出结果命令流 ********************************

复制代码