? ? ? ? 大家好，今天二狗準(zhǔn)備介紹一下如何利用matlab寫一些簡單的腳本去解決一些復(fù)雜的事情，比如為大型有限元軟件Ansys服務(wù)(類比)包括數(shù)值模擬的前處理和后處理，都大有作為。

? ? ? ? 用過Ansys(沒有Ctrl+z)的童鞋知道，它有GUI和命令流倆個操作模式，優(yōu)缺點就不多說了。前處理就是建模，點線面體幾何拓撲。如果有幾千個坐標(biāo)點(XYZ)，成點很麻煩，線面體更麻煩。簡單，可回退，可傳播分享的一種方式就是利用腳本批量寫入命令流到文本文件(Txt-dat-sca-csv...)

? ? ? ?下面是導(dǎo)入點的源碼，歡迎參考:

%三維數(shù)據(jù)到ANSYS[filename,pathname,filterindex]=uigetfile({'*.csv';'*.xlsx';'*.mat'},'選擇數(shù)據(jù)文件'); str=[pathname filename]; if filename==0 % 取消return;endif filterindex==2 c_temp=importdata(str);c_close=c_temp.Sheet1;elsec_temp=importdata(str);c_close=c_temp;end[hang,lie]=size(c_close);if lie==3mingzinn=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 endfprintf(fp,'%s\r\n','!File generation using Matlab.m');fclose(fp);sound(sin(2*pi*25*(1:4000)/100));elsemsgbox('數(shù)字矩陣必須為x*3形態(tài)'); end

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

[filename,pathname,filterindex]=uigetfile({'*.sca';'*.csv';'*.xlsx';'*.mat'},'選擇數(shù)據(jù)文件'); str=[pathname filename]; if filename==0 % 取消return;endif 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==3mingzinn=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 endfprintf(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));elsemsgbox('數(shù)字矩陣必須為x*3形態(tài)(不算 *.sca)'); end

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

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

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

? ? ? ?對于Ansys，懂的看它的Log日志文件很重要，但是此文件里有很多的無效命令，就是前面加了!，感嘆號是注釋的意思，如何把這些非實質(zhì)性操作清空呢，以下是二狗寫的幾句代碼：

[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);endfprintf(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)確高效的實現(xiàn)，也可直接把其打包成Exe，若需要。嗯。為達目的，“不擇手段”。

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

? ? 歡迎大家評論，關(guān)注，轉(zhuǎn)發(fā)！

? ? ? ?您的贊賞就是我們前進的動力

總結(jié)

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

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