【原创】MATLAB中曲面转入ANSYS

paradiseboy

来自动力学与控制技术论坛：http://www.dytrol.com/
  
前段时间因课题需要，写了个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提示有几个错误，有些样条线和关键点不能删除，那是需要保留的。
  
原程序请见http://www.dytrol.com/dispbbs.asp?boardID=11&ID=4550&page=1

paradiseboy

该程序还有许多不完善的地方，比如曲面特性太复杂，在ANSYS中不能成功生成曲面等问题，欢迎大家踊跃讨论啊！

kay_guo

这个程序在matlab6.1中似乎有点问题
总是说mat2ans没定义?
不知为啥?我现在也想在这方面做些尝试！

paradiseboy

kay_guo wrote:
这个程序在matlab6.1中似乎有点问题
总是说mat2ans没定义?
不知为啥?我现在也想在这方面做些尝试！

  
总于有对这方面感兴趣的了，我现在在写把matlab中封闭曲面导入ANSYS能生成实体的程序，欢迎讨论！
mat2ans是我自己写的程序，代码如下：
  
function mat2ans(filename,x,y,z)
%  把MATLAB计算的曲线曲面导入ANSYS
%  作者：xiegang
%  2004-7-13  
%  Henan Polytechnic University
  
%  www.dytrol.com
  
if(ischar(filename)==0)
     error('Invalid filename');
end
  
if(ndims(z)~=2)
     error('Variable z must be a 2-dimensional array');
end
  
if any((size(x)~=size(z))|(size(y)~=size(z)))  % size of x or y does not size of z
     if((length(x)==1)&(length(y)==1)) % Must be specifying dx and dy, so make vectors
         dx=x;
         dy=y;
         x=((1:size(z,2))-1)*dx;
         y=((1:size(z,1))-1)*dy;
     end
      
     if((length(x)==size(z,2))&(length(y)==size(z,1))) % Must be specifying vectors
         xvec=x;
         yvec=y;
         [x,y]=meshgrid(xvec,yvec);
     else
         error('Unable to resolve x and y variables');
     end
end
  
fid=fopen(filename,'w');
  
if(fid==-1)
     error(sprintf('Unable to write to %s',filename));
end
  
[m,n]=size(z);
x=reshape(x',m*n,1);
y=reshape(y',m*n,1);
z=reshape(z',m*n,1);
fprintf(fid,'/prep7\r\n');  % 进入前处理
% 创建关键点
for i=1:m*n
     fprintf(fid,'k,%4.d,%.5f,%.5f,%.5f,\r\n',i,x(i),y(i),z(i));
end
% 创建样条线
for i=1:m
     fprintf(fid,'FLST,3,%3.d,3\r\n',n);
     for j=1:n
         fprintf(fid,'FITEM,3,%4.d,\r\n',j+(i-1)*n);
     end
     fprintf(fid,'BSPLIN, ,P51X\r\n');
end
% 由蒙皮技术创建面
fprintf(fid,'FLST,2,%2.d,4\r\n',m);
for i=1:m
     fprintf(fid,'FITEM,2,%d\r\n',i);
end
fprintf(fid,'ASKIN,P51X\r\n');
fprintf(fid,'ldele,all\r\n');   
fprintf(fid,'kdele,all\r\n');   
fprintf(fid,'aplot\r\n');
fprintf(fid,'NUMCMP,ALL\r\n');   
fprintf(fid,'FINISH\r\n');
  
fclose(fid);

kay_guo

我先研究一下，我借了一本《matlab有限元分析与应用》其中有好多关于节点的计算方法。等我看完了再请教。

alexqxp

请问paradiseboy对下面这个问题怎么看。

http://www.simwe.com/forum/viewthread.php?tid=512019
http://www.simwe.com/forum/viewthread.php?tid=516785
paradiseboy很牛啊。
头衔：首席信息官CIO
  等级：校长
  文章：1779
  积分：2
  现金：2477
  魅力：0
  注册：2004-2-25

感觉动力学与控制技术论坛和这里象是兄弟帮。两边交流很多。

xiexie

非常感谢这位大侠提供原程序!!!我还没研究你的那个原文件,照着试了一把,但在ansys里出错了,说无法计算线的方向,还有线的点相距太近!以后多向大家学习!!

civilman狼人

我按照斑竹的方法进行了操作，可是遇到了问题，可以生成点和线，可是最后不能生成面．提示说：点太近或者线太近了．
请指教．我正好需要这方面的材料
谢谢斑竹的分享啊

feitianlong

不错的帖子