大家好，今天二狗準(zhǔn)備介紹一下如何利用matlab寫一些簡單的腳本去解決一些復(fù)雜的事情，比如為大型有限元軟件Ansys服務(wù)（類比）包括數(shù)值模擬的前處理和后處理，都大有作為。用過Ansys（沒有Ctrl+z）的童鞋知道，它有GUI和命令流倆個操作模式，優(yōu)缺點(diǎn)就不多說了。前處理就是建模，點(diǎn)線面體幾何拓?fù)?。如果有幾千個坐標(biāo)點(diǎn)（XYZ），成點(diǎn)很麻煩，線面體更麻煩。簡單，可回退，可傳播分享的一種方式就是利用腳本批量寫入命令流到文本文件（Txt-dat-sca-csv…）

下面是導(dǎo)入點(diǎn)的源碼，歡迎參考: %三維數(shù)據(jù)到ANSYS [filename,pathname,filterindex]=uigetfile({'*.csv';'*.xlsx';'*.mat'},'選擇數(shù)據(jù)文件');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'];%擴(kuò)展名可改 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:hangif i==1fprintf(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);elsestr2=[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('數(shù)字矩陣必須為x*3形態(tài)'); end

這個是導(dǎo)入點(diǎn)自動成線的源碼:

[filename,pathname,filterindex]=uigetfile({'*.sca';'*.csv';'*.xlsx';'*.mat'},'選擇數(shù)據(jù)文件');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'];%擴(kuò)展名可改 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:hangif i==1fprintf(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);elsestr2=[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('數(shù)字矩陣必須為x*3形態(tài)（不算 *.sca)'); end

成面與體，看自己實(shí)際需求，也可直接在ansys里寫命令流。Ansys求解后的結(jié)果該如何導(dǎo)出呢，舉個栗子，比如各節(jié)點(diǎn)的x方向位移（需要把Ansys里對應(yīng)的List的文本復(fù)制出來，比如到csv），源碼如下:

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

這樣得到的就是數(shù)據(jù)文件而不是字符串，且沒有空白，沒有注釋。

對于Ansys，懂的看它的Log日志文件很重要，但是此文件里有很多的無效命令，就是前面加了!，感嘆號是注釋的意思，如何把這些非實(shí)質(zhì)性操作清空呢，以下是二狗寫的幾句代碼： [h,l]=size(log); xin=[]; nametxt=['log1','-Ansys','.dat']; fp=fopen(nametxt,'a'); for p=1:1:hjiancha= log(p,1);l8= strlength(jiancha);if l8>1weizhi=strfind(jiancha,'!');[k1,k2]=size(weizhi);if k1>0 && weizhi(1,1)==1elsexin=[xin;jiancha]; fprintf(fp,'%s\r\n',jiancha); endelseend 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));

是不是很方便哦⊙?⊙！

在學(xué)習(xí)工作中，有很多操作可以用代碼去準(zhǔn)確高效的實(shí)現(xiàn)，也可直接把其打包成Exe，若需要。嗯。為達(dá)目的，“不擇手段”。

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

總結(jié)

以上是生活随笔為你收集整理的Matlab给Ansys助攻的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。