研究生階段的一些經(jīng)驗總結(有點亂,希望有用)

1、研究生階段的一些經(jīng)驗總結（有點亂，希望有用）目錄Matlab 11. 將數(shù)組中的數(shù)據(jù)存儲在文本文件中 12. 畫圖的一些命令格式應用舉例（有點亂，但比較有用） 13. 在 Matlab 圖形中怎樣輸入特殊符號 ? 44. 將ansys計算得到的位移數(shù)據(jù)轉化為加速度 45. 常用函數(shù) 56. 計時 57. 自相關函數(shù) 58. 傅里葉變換 6ansys79. ansys的一些命令 710. 從 ANSYS 文件中提取有用數(shù)據(jù) 711. 梁截面的顯示開關是什么？ 812. 提取結構反力 913. /config,nres,20000 914. *dim,weiyiy,array,nnn,nm 915

2、. 循環(huán)結構和選擇結構 916. 命令 FLST 917. 已知結構的前二階固有頻率f 和阻尼比，求阻尼矩陣的系數(shù) 1018. 實驗 阻尼比和頻率計算 10Matlab1. 將數(shù)組中的數(shù)據(jù)存儲在文本文件中命令：save y2.txt y1-ASCII -TABS說明：將 y1 這個數(shù)組中的數(shù)據(jù)以字符形式存儲在文件 y2.txt 中2. 畫圖的一些命令格式應用舉例（有點亂，但比較有用）figure（2）;plot（yyy）;xlabel（ 結點序號 ）;ylabel（ 振型改變率（ %）;title（第1階振型改變率）axis tight;grid on;加注釋框plot(pl,p1,-*,pl

3、,p2,-d,pl,p3,-p,pl,p4,-h);xlabel( 損傷單元);ylabel(頻率改變率(% ) );title( 前四階頻率改變率 )lege nd(第一階,第二階,第三階,第四階);legend boxoff;axis tight;hold on;plot(hh,pinlv(1,:),-*k,LineWidth,2,MarkerSize,12)hold onplot(hh,pinlv(2,:),-ok,LineWidth,2,MarkerSize,12)hhh=legend(測點 1,測點 2);xlim(2 8);ylim(0 1);xlabel(工況,FontSize,

4、18);ylabel(COMAC,FontSize,18);set(gca,XTick,2:3:12,FontSize,18,YTick,0:0.1:1) set(hhh,FontSize,16,FontWeight,Bold);legend(boxoff)hh=0 12 ;plot(pinlv(1,:),hh,-*k,LineWidth,2,MarkerSize,12) hold on plot(pinlv(2,:),hh,-ok,LineWidth,2,MarkerSize,12) plot(pinlv(3,:),hh,-sk,LineWidth,2,MarkerSize,12)% hhh

5、=legend(工況 1,工況 4,工況 7,2);hhh=legend(工況 2,工況 5,工況 8,2);% axis tight;axis(-1 1 0 2)ylabel(層數(shù),FontSize,18);xlabel(第 2 階振型,FontSize,18); set(gca,XTick,-1:0.2:1,FontSize,18); set(gca,YTick,0:1:3,FontSize,18); set(hhh,FontSize,16,FontWeight,Bold); legend(boxoff)plot(tyyt,yy(:,1),-*k,tyyt,yy(:,2),-sk,tyyt

6、,yy(:,3),-dk,markersize,10,LineWidth,2);legend(j=1,j=2,j=3)xlabel( 計算次序 ,FontSize,18);xlabel( 損傷部位彈性模量減少百分比 /%,FontSize,18);ylabel(strcat(sprintf(itNRrm_%d_j(tau_0) 變化率 ,xv),/%),FontSize,18); set(gca,XTick,tyyt,FontSize,18);xlim(0 12);clc% text(,.Stri ng,$pi to Atau_Ogg$gg,.Position,.5 .5,.FontSize,

7、16)xlabel(sprintf(itNRrm_%d_j(tau_0),2)% Matlab 圖形中允許用 TEX 文件格式來顯示字符。使用bf，it，rm 表示黑體，% 斜體，正體字符，特別注意大括號 的用法。subplot(2,2,1);plot(zzz,s1)axis tight;xlabel time(s);ylabel m/(m/sA2);在一個大的畫面上，畫 2 行 2 列小圖，將數(shù)據(jù)充滿數(shù)軸，橫坐標單位是time(s) ，縱坐標單位是 m/(m/s2)figure(xv); 產(chǎn)生第 xv 個圖plot(tyyt,yy,-*); 畫曲線 legend(j=1,j=2) ；注釋框y

8、label(strcat(sprintf(NR%dj( t 0)變化率,xv),(%);上句涉及到兩個函數(shù)， sprintf 生成一個字符串， strcat 將字符串連接起來title(sprintf( 方法 1 的 NR%dj( t 0)變化率 ,xv);xlabel(計算次序);axis tight;set(gca,XTick,tyyt); 設置橫坐標的刻度 plot(tyyt,yy(:,1),-*k, markersize,10,LineWidth,2); 在畫曲線的同時，給出了線的寬度，顏色， 符號的大小xlabel(損傷位置寬度改變量/mm,FontSize,18);橫坐標的刻度值的

9、標度的字符大小是18ylabel(strcat(sprintf(itNRrm_%d_j(tau_0) 變化率 ,xv),/%),FontSize,18);在圖上縱坐標上寫一些希臘字母set(gca,XTick,tyyt,FontSize,18);xlim(0 12);x 軸的取值范圍是 0到 12hhh=legend(完好狀態(tài),彈性摸量減少3%，彈性摸量減少6%，彈性摸量減少9%,彈性摸量減 少 12%,2); 將注釋框的句柄給 hhhaxis(-0.62 2.1 0 3) ；給定了橫坐標和縱坐標的取值范圍 set(hhh,FontSize,16,FontWeight,Bold); 給定了注釋

10、框圖線的字符尺寸和字符的粗度legend（boxoff）;去掉注釋框的邊框 clc% text(，%Stri ng,$pi to Atau_0gg$gg,.%Positi on ,.5 .5,.%Fo ntSize,16)xlabel(spri ntf(itNRrm_%d(tau_0), 2)% Matlab圖形中允許用 TEX文件格式來顯示字符。使用 bf, it, rm表示黑體, %斜體，正體字符，特別注意大括號 的用法。3.在Matlab圖形中怎樣輸入特殊符號 ？答：上標用A（指數(shù)）下標用_（下劃線） 調用方式為：A任意字符, _任意字符,注意的使用!希臘字母等特殊字符用 加拼音如aal

11、pha,Bbeta , ygamma0 theta O ThetarGammaSdeltaADeltaExiSXineltaeepsil ongZzeta(imiuunut tau入lamdaALamdanpinPicrsigmaSSigma$phi Phipsi PsixchiwommigaQ Ommiga geq 不等于 neq gg正負 pm左箭頭 leftarrow 右箭頭 rightarrow 上箭頭uparrow體會以下兩例：figure, title(zetaA-tsi nt);figure, title(xchi_alphaA2 (3);4.將ansys計算得到的位移數(shù)據(jù)轉化為

12、加速度load s1s2;a1a2=s1s2;dt=1/208.33;for i=2:24999 a1a2(i,:)=(s1s2(i-1,:)+s1s2(i+1,:)-s1s2(i,:)*2)/dt/dt;endsave a1a2.mat a1a2plot(s1s2)figure(2);plot(a1a2);load s1s2;v1v2=s1s2;dt=1/208.33;for i=2:50000-1v1v2(i,1)=(s1s2(i+1,1)-s1s2(i-1,1)/2/dt;v1v2(i,2)=(s1s2(i+1,2)-s1s2(i-1,2)/2/dt;endv1v2(1,1)=(s1s2

13、(2,1)-s1s2(1,2)/dt;v1v2(1,2)=(s1s2(2,2)-s1s2(1,2)/dt; v1v2(50000,1)=(s1s2(50000,1)-s1s2(49999,1)/dt; v1v2(50000,2)=(s1s2(50000,2)-s1s2(49999,2)/dt;a1a2=v1v2;for i=2:50000-1a1a2(i,1)=(v1v2(i+1,1)-v1v2(i-1,1)/2/dt;a1a2(i,2)=(v1v2(i+1,2)-v1v2(i-1,2)/2/dt;enda1a2(1,1)=(v1v2(2,1)-v1v2(1,2)/dt;a1a2(1,2)=(

14、v1v2(2,2)-v1v2(1,2)/dt; a1a2(50000,1)=(v1v2(50000,1)-v1v2(49999,1)/dt; a1a2(50000,2)=(v1v2(50000,2)-v1v2(49999,2)/dt;save a1a2.mat a1a25. 常用函數(shù)eig 求特征值zeros 0 數(shù)組 abs 絕對值數(shù)組 n=round(t200/10); 四舍五入disp 顯示h,l=size(s1s2); 給出數(shù)組的行數(shù)和列數(shù)y=randn(ndzb,1); 隨機數(shù)yp=mean(y) ；平均數(shù)hp=load(sprintf(jhp%d.txt,jjj);將文本文件中的數(shù)

15、據(jù)給矩陣hp,否則將生成一個矩陣和文件同名。6. 計時tictoc7. 自相關函數(shù)function s0=zixg(s)s0=0;for j=ss0=s0+jA2;end8. 傅里葉變換chl=a1a2.mat; load(chl);qi=15000; y=a1a2(qi+1:qi+10000,2);% y=y/(100*7.47*45); t200=size(y);t200=t200(1); zz=1:t200;zzz=zz/208.33; plot(zzz,y) axis tight;xlabel time(s);ylabel m/sA2;mn=1;y=y(1:8192);% plot(y

16、)n=t200(1);x=1:n;x=x;ts=1/208.33;ws=2*pi/ts; m=n/2;mm=round(m); w=ws*(0:n-1)/n;w=w/2/pi;z=fft(y,n);% z=z.*conj(z)/n/n;zm=abs(z);zm=zm/n;% zm=zm/(2*pi);% zm=20*log10(zm);% zm=z;% zm=10*log10(zm); figure(2);plot(w(mn:mm),zm(mn:mm); xlabel(f(Hz) axis tight;a1,a2=max(zm);w(a2)ansys9. ansys 的一些命令/OUTPUT,

17、apinggu0929,txt, 將輸出窗口內容輸入到文件fini 推出當前模塊，到最高層/cle 重新開始計算/config,nres,30000/config,fsplit,750 將結果文件分成幾份，每份大小為標準的千分之 750 /INPUT,0STEPG1,TXT, 0 *dim,li,array,nnn,1*dim,bo,TABLE,nn,1,1,time,wave*tread,bo,suiji,txt, *CFOPEN,y22z28, , *VWRITE,weiyiy(1,22,1), weiyiy(1,23,1) ,weiyiy(1,24,1), weiyiy(1,25,1),

18、 weiyiy(1,26,1),weiyiy(1,27,1),weiyiy(1,28,1)(7E18.8E3)*CFOPEN,li, , *VWRITE,li(1,1,1), ,(E18.8E3)10. 從 ANSYS 文件中提取有用數(shù)據(jù)clc;clear;文注釋，且注釋最后以 gal 為結尾的情況 ch1=k17;ch2=strcat(ch1,.asc)ch3=strcat(ch1,.mat) fp=fopen(ch2,r); buff=fscanf(fp,%s,1) hang=3000000;lie=3; uii=0;ss=zeros(hang,lie);while(1)if strcmp

19、(buff,gal)=1 否即將離開注釋區(qū))% 注意， 本文件只適用于數(shù)據(jù)列數(shù)是 3，只有英% 要打開的文件名% 要打開的文件的擴展名是 ASC% 要存入的文件的擴展名是 MAT% 打開待處理文件% 讀取第一個字符串% 數(shù)據(jù)區(qū)總行數(shù)% 希望得到的數(shù)據(jù)列數(shù)% 數(shù)據(jù)要存入的矩陣% 第一次判斷特征字符串 （判斷是% 第二次判斷特征字符串 （判斷buff=fscanf(fp,%s,1);if strcmp(buff,Write)=1是否已經(jīng)離開注釋區(qū))ss(1,1)=str2num(buff);fbuff=fscanf(fp,%f,1);ss(1,2)=fbuff;fbuff=fscanf(fp,%f

20、,1);ss(1,3)=fbuff; ui=2;while(1)buff=fscanf(fp,%s,1); if (length(buff)=0) uii=1breakendif (abs(buff)PlotCtrlsStyleSize and Shape12. 提取結構反力 *GET,PAR,NODE,N,RF,FX *GET,PAR,NODE,N,RF,FY *GET,PAR,NODE,N,RF,FZ *GET,PAR,NODE,N,RF,MX *GET,PAR,NODE,N,RF,MY *GET,PAR,NODE,N,RF,MZ13. /config,nres,20000 在我以往動力計

21、算中，缺省的步驟是 1000，如果希望做更多的步數(shù)，就可利用這個命令， 其中 20000 是我在某次計算中動力計算的步數(shù)。fini/cle /config,nres,20000 /config,fsplit,750 fini ：退出當前模塊/cle :進行一個新的計算/config,fsplit,750 ：由于結果文件過大， ansys 將停止計算，因此，可以設置多個結果文件， 沒有結果文件是原來 4G 的千分之 750。14. *dim,weiyiy,array,nnn,nm定義一個數(shù)組，其類型是array，其名稱是 weiyiy，其位數(shù)維數(shù)是 nnn行，nm列*dim,bo,TABLE,n

22、n,1,1,time,wave定義一個數(shù)組，其類型是TABLE，其名稱是bo,其位數(shù)維數(shù)是nn行，1列，1頁?？梢圆逯涤嬎?，其行標名稱是time，其列標名稱是 wave*tread,bo,suiji,txt,從 suiji.txt 文件讀取數(shù)據(jù)給數(shù)組 bo *CFOPEN,s1s2, , *VWRITE,weiyiy(1,1,1), weiyiy(1,2,1) , , , (2E18.8E3) 打開一個文件,并且將變量按照指定的格式寫進文件15. 循環(huán)結構和選擇結構*do,m,1,nnn,1*enddo*if ,m,lt,20,then*else *endif16. 命令 FLSTFLST,2,56,4,ORDE,2FITEM,2,1FITEM,2,-56命令的格式是：FLST , NFIELD NARCTYPEOtype, LENGNFIELD表示NARGfe示被拾取對象的個數(shù)TYPE表示被拾取對象的類型 1 節(jié)點，2 單兀，3 關鍵點，4線，5面，6 體，Otype:數(shù)據(jù)順序 NOOR -無序；ORDER -有序；LENG數(shù)據(jù)項的行數(shù).命令格式：FITEM , NFIELD ITEM I