用matlab对fluent导出的结果进行处理而且可以将数据导出的例子(转贴)

绿荷雨丝

matlab画图能力很强，而且其可编程性给处理问题带来很大方便，我尝试将fluent导出的ASCII数据导入matlab处理，原代码如下：

1. %打开文件并读入数据
   
2. fid=fopen('post','r');
   
3. name2=input('input the file name\n','s')
   
4. fgetl(fid)
   
5. data1=fscanf(fid,'%f %f %f %f %f',[5 inf]);
   
6. fclose(fid);
   
7. fid=fopen(name2,'r');
   
8. fgetl(fid);
   
9. data2=fscanf(fid,'%f %f %f %f %f',[5 inf]);
   
10. fclose(fid);
    
11. %对数据按坐标升序排列
    
12. xudata1=sortrows(data1',2);
    
13. xudata2=sortrows(data2',2);
    
14. xudata1=xudata1';xudata2=xudata2';
    
15. %找出sym和avg数据
    
16. maxr=max(xudata1(3,:));stp=0.0002;%y坐标最大值，求avg时候的x步长
    
17. avg1=zeros(2,round(1+0.02/stp));avg1(1,:)=0:stp:0.02;avg2=avg1;%对avg坐标赋初值
    
18. j=1;j2=1;jj=1;%j记录avg的当前index；j2记录avg(j)中第一个加入的xudata的index；jj记录sym的当前index
    
19. for i=1:length(xudata1)
    
20.      if xudata1(3,i)>maxr-0.00004
    
21.      sym1(1,jj)=xudata1(2,i);sym1(2,jj)=xudata1(5,i);
    
22.      sym2(1,jj)=xudata1(2,i);sym2(2,jj)=xudata2(5,i);
    
23.      jj=jj+1;
    
24.      end
    
25.       mii=fix((xudata1(2,i)+stp/2)/stp)+1;
    
26.      avg1(2,mii)=avg1(2,mii)+xudata1(5,i);
    
27.      avg2(2,mii)=avg2(2,mii)+xudata2(5,i);
    
28.      if mii>j|i==length(xudata1)
    
29.      avg1(2,j)=avg1(2,j)/(i-j2);
    
30.      avg2(2,j)=avg2(2,j)/(i-j2);
    
31.      i-j2
    
32.      j=j+1;j2=i;
    
33.      end
    
34. end
    
35. %绘图
    
36. figure(1)
    
37. plot(sym1(1,:),sym1(2,:),sym2(1,:),sym2(2,:))
    
38. figure(2)
    
39. plot(sym1(1,:),sym1(2,:),sym2(1,:),sym2(2,:),avg1(1,:),avg1(2,:),avg2(1,:),avg2(2,:));
    
40. legend('sym1','sym2','avg1','avg2');
    
41. hold on
    
42. plot(sym1(1,:),sym2(2,:)./sym1(2,:),avg1(1,:),avg2(2,:)./avg1(2,:));
    
43. legend('symrate','avgrate');
    
44. fid=fopen(['tec' name2],'w');
    
45. fprintf(fid,'TITLE = "heat transfer coefficient" \nVARIABLES = "x/d", "yta or rate" \nZONE T="sym", I=%d, F=POINT\n',length(sym2));
    
46. fprintf(fid,'%12.6f  %12.6f',sym2);
    
47. fprintf(fid,'ZONE T="avg", I=%d, F=POINT\n',length(avg2));
    
48. fprintf(fid,'%12.6f  %12.6f',avg2);
    
49. fclose(fid)
    
50. %用tecplot格式输出
    
51. fid=fopen('tecpost','w');
    
52. fprintf(fid,'TITLE = "heat transfer coefficient" \nVARIABLES = "x/d", "yta or rate" \nZONE T="sym", I=%d, F=POINT\n',length(sym1));
    
53. fprintf(fid,'%12.6f  %12.6f',sym1);
    
54. fprintf(fid,'ZONE T="avg", I=%d, F=POINT\n',length(avg1));
    
55. fprintf(fid,'%12.6f  %12.6f',avg1);
    
56. fclose(fid)

复制代码

来自：http://wlflyy.blog.hexun.com/4722405_d.html

grta

呵呵，很感兴趣，能不能把画的图贴出来啊

jiangwenxue999

请问可不可以输出数据到文本中呀？
jiangwenxue999@sina.com

hyacinth

原帖由 jiangwenxue999 于 2007-6-14 13:44 发表

                               
登录/注册后可看大图

请问可不可以输出数据到文本中呀？
jiangwenxue999@sina.com

可以，你的主题贴中已经回复

lnlyxia

为什么我用fluent —Export输出数据并以ASCII形式保存，用Matlab打不开呢？

tmjking

在matlab中要用Import导入数据才可以的，此外你可以用记事本打开看看里面的格式是不是数组，否则导入是会出错的

xiong5527

fid=fopen('post','r');
'post'    是什么格式的文件？我的fluent导出的数据文件不能用这句语言

tmjking

'post'是文件全名，无后缀名。你把你输出的数据文件全名改写'post'才能用使用这个程序。
再不明白就查matlab的帮助文件看吧

guohf

谢谢，正需要，非常不错的贴子。

chbo76

呵呵。:@)
用matlab对fluent导出的结果进行处理，要求我们也要会matlab才行啊。

taojunhust

我觉得六楼的兄弟对matlab用得很熟 啊。

tmjking

很熟说不上了，以前用过几年，也就懂些皮毛而已。
用这种东西，经验很重要，当然对于数理逻辑上的概念也很重要。

yuanweizhao

请问fluent能不能导出声学模块计算的结果呢，在EXPORT我找不到这方面的结果。另外，能讲一下贴子中，针对哪种输出呢？

seaharrier

yuanweizhao 发表于 2012-2-27 12:52

                               
登录/注册后可看大图

请问fluent能不能导出声学模块计算的结果呢，在EXPORT我找不到这方面的结果。另外，能讲一下贴子中，针对哪 ...

只要是数组就可以吧，符合一定标准格式

以前玩过一阵这个，不过是从CFX导出的数据，

现在竟忘了，

这个帖子不错，
多个软件结合应用的好例子。