Matlab给Ansys助攻

巴山学长

发布于 2019-07-15 15:14:12

8150

发布于 2019-07-15 15:14:12

大家好，今天二狗准备介绍一下如何利用matlab写一些简单的脚本去解决一些复杂的事情，比如为大型有限元软件Ansys服务（类比）包括数值模拟的前处理和后处理，都大有作为。

用过Ansys（没有Ctrl+z）的童鞋知道，它有GUI和命令流俩个操作模式，优缺点就不多说了。前处理就是建模，点线面体几何拓扑。如果有几千个坐标点（XYZ），成点很麻烦，线面体更麻烦。简单，可回退，可传播分享的一种方式就是利用脚本批量写入命令流到文本文件（Txt-dat-sca-csv...）

下面是导入点的源码，欢迎参考:

%三维数据到ANSYS
[filename,pathname,filterindex]=uigetfile({'*.csv';'*.xlsx';'*.mat'},'选择数据文件');
    str=[pathname filename]; 
if filename==0 % 取消
return;
end
if filterindex==2 
c_temp=importdata(str);
c_close=c_temp.Sheet1;
else
c_temp=importdata(str);
c_close=c_temp;
end
[hang,lie]=size(c_close);
if lie==3
mingzinn=char(filename);
ll=strfind(mingzinn,'.');
ll0=mingzinn(1:ll-1);
nametxt=[ll0,'2ANSYS','.dat'];%扩展名可改
fp=fopen(nametxt,'a');
str111='k , ';
str222=' , ';
str333=' ';
str444='LP1 = _return';
%str555='L,LP1,LP2';
%str666='LP1 = LP2';
%str777='LP2 = _return';
for i=1:1:hang
   if i==1
    fprintf(fp,'%s\r\n','/prep7');  
    fprintf(fp,'%s\r\n','CSYS,0');
    str1=[str111,num2str(i),str222,num2str(c_close(1,1)),str222,num2str(c_close(1,2)),str222,num2str(c_close(1,3))];
    fprintf(fp,'%s\r\n',str1);
   else
    str2=[str111,num2str(i),str222,num2str(c_close(i,1)),str222,num2str(c_close(i,2)),str222,num2str(c_close(i,3))];   
    fprintf(fp,'%s\r\n',str2);   
   end         
end
fprintf(fp,'%s\r\n','!File generation using Matlab.m');
fclose(fp);
sound(sin(2*pi*25*(1:4000)/100));
else
msgbox('数字矩阵必须为x*3形态');    
end

这个是导入点自动成线的源码:

[filename,pathname,filterindex]=uigetfile({'*.sca';'*.csv';'*.xlsx';'*.mat'},'选择数据文件');
    str=[pathname filename]; 
if filename==0 % 取消
return;
end
if filterindex==3 
c_temp=importdata(str);
c_close=c_temp.Sheet1;
elseif filterindex==1 
c_temp=importdata(str);
c_close=c_temp.data;  
c_close=c_close(:,2:4);
elseif filterindex==2 
c_temp=importdata(str);
c_close=c_temp;
elseif filterindex==4 
c_temp=importdata(str);
c_close=c_temp;
end
[hang,lie]=size(c_close);
if lie==3
mingzinn=char(filename);
ll=strfind(mingzinn,'.');
ll0=mingzinn(1:ll-1);
nametxt=[ll0,'-Ansys','.dat'];%扩展名可改
fp=fopen(nametxt,'a');
str111='K,, ';
str222=' , ';
str333=' ';
str444='LP1 = _return';
str555='L,LP1,LP2';
str666='LP1 = LP2';
str777='LP2 = _return';
for i=1:1:hang
   if i==1
    fprintf(fp,'%s\r\n','/prep7');  
    str1=[str111,num2str(c_close(1,1)),str222,num2str(c_close(1,2)),str222,num2str(c_close(1,3)),str333];   
    fprintf(fp,'%s\r\n',str1);
    fprintf(fp,'%s\r\n',str444);
   else
    str2=[str111,num2str(c_close(i,1)),str222,num2str(c_close(i,2)),str222,num2str(c_close(i,3)),str333];   
    fprintf(fp,'%s\r\n',str2);   
    fprintf(fp,'%s\r\n',str777);   
    fprintf(fp,'%s\r\n',str555);
    fprintf(fp,'%s\r\n',str666);
   end         
end
fprintf(fp,'%s\r\n','!File generation using Matlab.m');
fprintf(fp,'%s\r\n','!By 阿辉');
fclose(fp);
sound(sin(2*pi*25*(1:4000)/100));
else
msgbox('数字矩阵必须为x*3形态（不算 *.sca)');    
end

成面与体，看自己实际需求，也可直接在ansys里写命令流。Ansys求解后的结果该如何导出呢，举个栗子，比如各节点的x方向位移（需要把Ansys里对应的List的文本复制出来，比如到csv），源码如下:

% load a88
[hang,lie]=size(a88);
ping=[];
for i=1:1:hang
jiancha=a88(i,1);
kankan=str2num(jiancha);
if isempty(kankan)==1
    
else
    ping=[ping;kankan]; 
end
h=waitbar(i/hang);
end
delete(h);
msgbox('OK');

这样得到的就是数据文件而不是字符串，且没有空白，没有注释。

对于Ansys，懂的看它的Log日志文件很重要，但是此文件里有很多的无效命令，就是前面加了!，感叹号是注释的意思，如何把这些非实质性操作清空呢，以下是二狗写的几句代码：

[h,l]=size(log);
xin=[];
nametxt=['log1','-Ansys','.dat'];
fp=fopen(nametxt,'a');
for p=1:1:h
   jiancha= log(p,1);
   l8= strlength(jiancha);
   if l8>1
      weizhi=strfind(jiancha,'!');
      [k1,k2]=size(weizhi);
     if k1>0 && weizhi(1,1)==1
         
     else
       xin=[xin;jiancha]; 
       fprintf(fp,'%s\r\n',jiancha);  
     end
   else
       
   end 
   jdt=waitbar(p/h);
end
fprintf(fp,'%s\r\n','!File idealization using Matlab.m');%写下来源  嘿嘿
delete(jdt);
fclose(fp);
sound(sin(2*pi*25*(1:4000)/100));

之前是：