|
楼主 |
发表于 2004-11-13 22:23:18
|
显示全部楼层
请教用Tecplot作图问题?
不好意思,还是没搞明白,以下是我加的代码,那位高人能帮我指点一下?谢谢
1)meshA.dat文件是圆锥弹体某截面的空间网格;
2)meshB.dat文件是圆锥头的空间网格;
以下是生成网格的代码,
其中: X=R(1,J,K,I),Y=R(2,J,K,I),Z=R(3,J,K,I)为空间坐标的三个方向值;
xo(i,k),yo(i,k)为物面坐标。
当我安以下程序代码生成文件时,用tecplot作图时画的都只有一条直线,但去掉以下两行代码:
write(15,*)';TITLE="meshA.dat"';
write(15,*)';ZONE j=';,j1 ,';k=';,kk1 ,';F=';,j1*kk1
则可画出图形,但图形中第一个只作出径向的线没有周向线,第二个画的象螺旋线似的,不知咱回事好像只是把以下计算的各坐标点简单的想连???以前从没有用过这个软件,还望指教!!!谢谢!
OPEN(15,FILE=';meshA.dat';,status=';unknown';)
I=21
write(15,*)';TITLE="meshA.dat"';
write(15,*)';VARIABLES="x","y","U","V",""';
write(15,*)';ZONE j=';,j1 ,';k=';,kk1 ,';F=';,j1*kk1
DO K=1,kk1
DO J=1,j1
Z=R(3,J,K,I)
Y=R(2,J,K,I)
WRITE(15,*)Z,Y,1.0
enddo
enddo
CLOSE(15)
OPEN(9,FILE=';meshB.dat';,STATUS=';UNKNOWN';)
write(9,*) ';TITLE="meshB.dat"';
write(9,*) ';variables="x","y","z"';
write(9,*) ';zone k=';,kk1,'; j=';,j1,'; i=';,is1
DO I=1,is1
DO j=1,j1
DO k=1,Kk1
WRITE(9,*) R(1,J,K,I), R(2,J,K,I), R(3,J,K,I)
enddo
enddo
enddo
CLOSE(9) |
|