Matlab給Ansys助攻

大家好，今天二狗準備介紹一下如何利用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));

是不是很方便哦⊙∀⊙！

在學習工作中，有很多操作可以用代碼去準確高效的實現，也可直接把其打包成Exe，若需要。嗯。爲達目的，“不擇手段”。

二狗就寫這麼多，文筆不好，見諒，有錯誤的地方請告訴二狗。解釋一下，提到的Ansys是支持數據導入導出的，二狗覺得，不好用。這些操作用vb、c、Python、office的宏等等都可以實現，但是matlab，除了生孩子，無所不能，哈哈ヾﾉ≧∀≦)o