24小时热门版块排行榜

返回列表

本帖产生 1 个计算强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

onesupeng

金虫 (职业作家)

计算强帖: 13
应助: 256 (大学生)
贵宾: 1.36
金币: 2336.2
散金: 9224
红花: 92
帖子: 4583
在线: 1303.8小时
虫号: 394701
注册: 2007-06-07
专业: 流体力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 计算强帖+1, 多谢指导 2013-05-13 11:39:30
zhanghuang10: 金币+5, ★★★很有帮助 2013-05-13 21:09:14

我列一个详细的tecplot的规矩给你看看吧。
打开open(unit=70,file='field.dat')
首先是结构网格的输出：
二维时：
write(70,*)'VARIABLES="X","Y","U","V","P"'
write(70,*)'ZONE F=POINT, I=',nx,', J=',ny
   do j=1,ny
   do i=1,nx
      write(70,*)xc(i),yc(j),u(i,j),v(i,j),p(i,j)
   enddo
   enddo
三维类似：
write(70,*)'VARIABLES="X","Y","Z","U","V","W","P"'
   write(70,*)'ZONE F=POINT, I=',nx,', J=',ny,' K=',nz
   do k=1,nz
   do j=1,ny
   do i=1,nx
      write(70,*)xc(i),yc(j),zc(k),u(i,j,k),v(i,j,k),w(i,j,k),p(i,j,k)
   enddo
   enddo
   enddo
非结构网格：
二维4边形单元
    write(70,*)'variables= "x" "y" "u" "v" "p" "w" '
    write(70,*)'zone f=fepoint et=quadrilateral n=',np,' e=',Nequd
    do i=1,np
   write(70,*)x(i),y(i),u(i),v(i),p(i),w(i)
    enddo
do ee=1,Nequd
   e=Equd(ee)
   write(70,*)(ep(e,j),j=1,4)
enddo

三角形单元：
write(70,*)'zone f=fepoint et=triangle n=',np,' e=',Netri
    do i=1,np
   write(70,*)x(i),y(i),u(i),v(i),p(i),w(i)
    enddo
do ee=1,Netri
   e=Etri(ee)
   write(70,*)(ep(e,j),j=1,3)
enddo

其他情况类似。多google、baidu和看手册，很容易找到答案。主要是动手！！！