找回密码
 注册
查看: 5177|回复: 4

关于AIAA高升力Workshop的网格问题

[复制链接]
发表于 2010-9-12 12:08:40 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
我下载了这个网站( http://hiliftpw.larc.nasa.gov/ )提供的网格数据,想拿来参考一下,结果发现按照PLOT3D格式根本打不开(结构OVERSET),不知道那位高人知道怎么打开该网站提供的网格数据,请不吝赐教!先谢谢了!
发表于 2010-9-13 11:45:57 | 显示全部楼层

关于Plot3D网格文件格式

我看了一下。
如果您下载的是Str-Overset-A-v1等中的数据文件的话,文件格式是Plot3d, double-precision, blanked, big endian,没错。这也是NASA的Overflow软件采用的格式。
您用Tecplot中Plot3D格式可以读取,只是需要将其中的File Structure选项手工设置为Overflow风格。
如果您要自己编程序读取,我还可以告诉您具体编码方法。
下图是Str-Overset-A-v1中Config1之coarse网格,供参考。
Str-Overset-A-v1-Config1-coarse.jpg
 楼主| 发表于 2010-9-13 13:20:30 | 显示全部楼层

回复 2# 白文 的帖子

非常感谢白文老师的回复,真的是太感谢了!如果方便的话麻烦您帮我看看我写的这个程序那里有问题呢?死活就是读不出来,急死了!
   dimension imax(35),jmax(35),kmax(35)
       dimension x(300,110,200),y(300,110,110),z(300,110,110),
     1                    iblank(300,110,110)
       open(1,file='grid.in',form='unformatted')
       open(2,file='mesh.dat')
       write(2,*)'variables=x,y,z,iblank'
       iunit = 1         
       read(iunit)nblock
       read(iunit)(imax(n),jmax(n),kmax(n),n=1,nblock)         
       do 10 n=1,nblock
       read(iunit) (((x(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
     1            ,(((y(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
     2            ,(((z(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))
     3       ,(((iblank(i,j,k),i=1,imax(n)),j=1,jmax(n)),k=1,kmax(n))         
       write(2,*)'zone i=',imax(n),' j=',jmax(n),' k=',kmax(n)
       do k = 1,kmax(n)
       do j = 1,jmax(n)
       do i = 1,imax(n)
         write(2,*)x(i,j,k),y(i,j,k),z(i,j,k),iblank(i,j,k)
       end do
       end do
       end do
10     continue
       close(1)
       close(2)
       end
发表于 2010-9-13 13:59:15 | 显示全部楼层

可能还是文件格式的问题

将那个
open(1,file='grid.in',form='unformatted')
改作
open(1,file='grid.in',form='unformatted',convert='big_endian')
试一下。
老师就不敢当了!
 楼主| 发表于 2010-9-13 15:34:31 | 显示全部楼层
白文老师:非常感谢,已经按照您的方法修改成功了,我确实是不懂这个big_endian!
再次感谢!
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表