[转帖]MATLAB中曲面转入ANSYS程序(附程序)

suffer

前段时间因课题需要，写了个MATLAB程序，把MATLAB中计算的曲面导入ANSYS，以作进一步处理。

稍作修改，使其更有通用性。因水平有限，程序还有很多不足之处，请大家验证讨论。

例：

t=0:pi/10:2*pi;
[x,y,z]=cylinder(2+cos(t));

mat2ans('cylinder.lgw',x,y,z);

运行ANSYS,在命令窗口输入：/input,cylinder.lgw 得到ansys中的曲面，该过程中ansys提示有几个错误，有些样条线和关键点不能删除，那是需要保留的。

suffer

附源程序：

1. function mat2ans(filename,x,y,z)
   
2. % 把MATLAB计算的曲线曲面导入ANSYS
   
3. % author：paradiseboy
   
4. % 2004-7-13
   
5. % Henan Polytechnic University
   
6. 
   
7. % www.dytrol.com
   
8. 
   
9. if(ischar(filename)==0)
   
10. error('Invalid filename');
    
11. end
    
12. 
    
13. if(ndims(z)~=2)
    
14. error('Variable z must be a 2-dimensional array');
    
15. end
    
16. 
    
17. if any((size(x)~=size(z))|(size(y)~=size(z))) % size of x or y does not size of z
    
18. if((length(x)==1)&(length(y)==1)) % Must be specifying dx and dy, so make vectors
    
19. dx=x;
    
20. dy=y;
    
21. x=((1:size(z,2))-1)*dx;
    
22. y=((1:size(z,1))-1)*dy;
    
23. end
    
24. 
    
25. if((length(x)==size(z,2))&(length(y)==size(z,1))) % Must be specifying vectors
    
26. xvec=x;
    
27. yvec=y;
    
28. [x,y]=meshgrid(xvec,yvec);
    
29. else
    
30. error('Unable to resolve x and y variables');
    
31. end
    
32. end
    
33. 
    
34. fid=fopen(filename,'w');
    
35. 
    
36. if(fid==-1)
    
37. error(sprintf('Unable to write to %s',filename));
    
38. end
    
39. 
    
40. [m,n]=size(z);
    
41. x=reshape(x',m*n,1);
    
42. y=reshape(y',m*n,1);
    
43. z=reshape(z',m*n,1);
    
44. fprintf(fid,'/prep7\r\n'); % 进入前处理
    
45. % 创建关键点
    
46. for i=1:m*n
    
47. fprintf(fid,'k,%4.d,%.5f,%.5f,%.5f,\r\n',i,x(i),y(i),z(i));
    
48. end
    
49. % 创建样条线
    
50. for i=1:m
    
51. fprintf(fid,'FLST,3,%3.d,3\r\n',n);
    
52. for j=1:n
    
53. fprintf(fid,'FITEM,3,%4.d,\r\n',j+(i-1)*n);
    
54. end
    
55. fprintf(fid,'BSPLIN, ,P51X\r\n');
    
56. end
    
57. % 由蒙皮技术创建面
    
58. fprintf(fid,'FLST,2,%2.d,4\r\n',m);
    
59. for i=1:m
    
60. fprintf(fid,'FITEM,2,%d\r\n',i);
    
61. end
    
62. fprintf(fid,'ASKIN,P51X\r\n');
    
63. fprintf(fid,'ldele,all\r\n');
    
64. fprintf(fid,'kdele,all\r\n');
    
65. fprintf(fid,'aplot\r\n');
    
66. fprintf(fid,'NUMCMP,ALL\r\n');
    
67. fprintf(fid,'FINISH\r\n');
    
68. 
    
69. fclose(fid);

复制代码

suffer

>> z=x.*exp(-x.^2-y.^2);
>> [u,v,w]=surfnorm(x,y,z);
>> surf(x,y,z); axis equal
>> mat2ans('curve.lgw',x,y,z);

suffer

>> [x,y]=meshgrid(-3:.5:3,-3:.1:3);
>> z=peaks(x,y);
>> mat2ans('curve.lgw',x,y,z);

这个曲面特性太复杂，在ANSYS中不能由一张曲面构成。
所以导入ansys中会出错，还请高人多指点，谢谢。

心灯

我个人认为应该是对于复杂图形，matlab输出的时候无法把握住图形的全部特征，把需要的全部数据都属出来，写入到ansys命令流中，造成生成面时的错误。
不知道对不？

hunter_009

matlab中生成的曲面能不能导入autocad中去呢?

tyut_syz

太好了，非常感谢！！！

chbo76

非常感谢。:lol