研究生階段的一些經(jīng)驗(yàn)總結(jié)2091231

1、研究生階段的一些經(jīng)驗(yàn)總結(jié)2091231目錄1. 將數(shù)組中的數(shù)據(jù)存儲(chǔ)在文本文件中12.畫圖的一些命令格式應(yīng)用舉例（有點(diǎn)亂，但比較有用）13. 在Matlab圖形中怎樣輸入特殊符號(hào)?44. 將ansys計(jì)算得到的位移數(shù)據(jù)轉(zhuǎn)化為加速度45. 常用函數(shù)56. 計(jì)時(shí)57. ansys的一些命令58. 自相關(guān)函數(shù)69. 從ANSYS文件中提取有用數(shù)據(jù)610. 傅里葉變換711. 梁截面的顯示開關(guān)是什么？812. 提取結(jié)構(gòu)反力813. /config,nres,20000814. *dim,weiyiy,array,nnn,nm915. 循環(huán)結(jié)構(gòu)和選擇結(jié)構(gòu)916. 命令FLST917. 已知結(jié)構(gòu)的前二階固有

2、頻率f和阻尼比，求阻尼矩陣的系數(shù)918. 實(shí)驗(yàn) 阻尼比和頻率計(jì)算10 1. 將數(shù)組中的數(shù)據(jù)存儲(chǔ)在文本文件中命令：save y2.txt y1 -ASCII TABS說明：將y1這個(gè)數(shù)組中的數(shù)據(jù)以字符形式存儲(chǔ)在文件y2.txt中2.畫圖的一些命令格式應(yīng)用舉例（有點(diǎn)亂，但比較有用）figure(2);plot(yyy);xlabel('結(jié)點(diǎn)序號(hào)');ylabel('振型改變率（%）');title('第1階振型改變率')axis tight;grid on;加注釋框plot(pl,p1,'-*',pl,p2,'-d&

3、#39;,pl,p3,'-p',pl,p4,'-h');xlabel('損傷單元');ylabel('頻率改變率（%）');title('前四階頻率改變率')legend('第一階','第二階','第三階','第四階');legend boxoff;axis tight;hold on;plot(hh,pinlv(1,:),'-*k','LineWidth',2,'MarkerSize',12)hold

4、onplot(hh,pinlv(2,:),'-ok','LineWidth',2,'MarkerSize',12)hhh=legend('測點(diǎn)1','測點(diǎn)2');xlim(2 8);ylim(0 1); xlabel('工況','FontSize',18); ylabel('COMAC','FontSize',18); set(gca,'XTick',2:3:12,'FontSize',18,'YTick'

5、;,0:0.1:1)set(hhh,'FontSize',16,'FontWeight','Bold'); legend('boxoff')hh=0 1 2 ;plot(pinlv(1,:),hh,'-*k','LineWidth',2,'MarkerSize',12)hold onplot(pinlv(2,:),hh,'-ok','LineWidth',2,'MarkerSize',12)plot(pinlv(3,:),hh,'

6、;-sk','LineWidth',2,'MarkerSize',12)% hhh=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,

7、9;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,yy(:,3),'-dk','markersize',

8、10,'LineWidth',2);legend('j=1','j=2','j=3')xlabel('計(jì)算次序','FontSize',18);xlabel('損傷部位彈性模量減少百分比/%','FontSize',18);ylabel(strcat(sprintf('itNRrm_%d_j(tau_0)變化率',xv),'/%'),'FontSize',18);set(gca,'XTick',tyyt,

9、'FontSize',18);xlim(0 12);clc% text(,.% 'String','$pi to tau_0gg$gg',.% 'Position',.5 .5,.% 'FontSize',16) xlabel(sprintf('itNRrm_%d_j(tau_0)',2)% Matlab圖形中允許用TEX文件格式來顯示字符。使用bf，it，rm表示黑體，% 斜體，正體字符，特別注意大括號(hào) 的用法。subplot(2,2,1);plot(zzz,s1) axis tight;xlab

10、el time(s);ylabel m/(m/s2);在一個(gè)大的畫面上，畫2行2列小圖，將數(shù)據(jù)充滿數(shù)軸，橫坐標(biāo)單位是time(s)，縱坐標(biāo)單位是m/(m/s2)figure(xv);產(chǎn)生第xv個(gè)圖plot(tyyt,yy','-*');畫曲線legend('j=1','j=2')；注釋框ylabel(strcat(sprintf('NR%dj(0)變化率',xv),'(%)');上句涉及到兩個(gè)函數(shù)，sprintf生成一個(gè)字符串，strcat將字符串連接起來title(sprintf('方法1的NR%

11、dj(0)變化率',xv);xlabel('計(jì)算次序');axis tight;set(gca,'XTick',tyyt);設(shè)置橫坐標(biāo)的刻度plot(tyyt,yy(:,1),'-*k', 'markersize',10,'LineWidth',2);在畫曲線的同時(shí)，給出了線的寬度，顏色，符號(hào)的大小xlabel('損傷位置寬度改變量/mm','FontSize',18);橫坐標(biāo)的刻度值的標(biāo)度的字符大小是18ylabel(strcat(sprintf('itNRrm_%

12、d_j(tau_0)變化率',xv),'/%'),'FontSize',18);在圖上縱坐標(biāo)上寫一些希臘字母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)；給定了橫坐標(biāo)和縱坐標(biāo)的取值范

13、圍set(hhh,'FontSize',16,'FontWeight','Bold'); 給定了注釋框圖線的字符尺寸和字符的粗度legend('boxoff')；去掉注釋框的邊框clc% text(,.% 'String','$pi to tau_0gg$gg',.% 'Position',.5 .5,.% 'FontSize',16) xlabel(sprintf('itNRrm_%d_j(tau_0)',2)% Matlab圖形中允許用TEX文件格

14、式來顯示字符。使用bf，it，rm表示黑體，% 斜體，正體字符，特別注意大括號(hào) 的用法。3. 在Matlab圖形中怎樣輸入特殊符號(hào)?答: 上標(biāo)用 (指數(shù))  下標(biāo)用 _(下劃線)調(diào)用方式為:任意字符, _任意字符, 注意的使用!希臘字母等特殊字符用 加拼音 如: alpha,    beta ,    gamma    theta    Theta Gamma    delta    Delta       xi&#

15、160;         Xi elta         epsilong      zeta       miu      nu       tau lamda       Lamda       pi    

16、60;   Pi       sigma    Sigma phi       Phi      psi      Psi       chi    ommiga      Ommiga< leq       > geq   &

17、#160; 不等于 neq      << ll     >> gg正負(fù) pm   左箭頭 leftarrow   右箭頭 rightarrow  上箭頭 uparrow體會(huì)以下兩例:figure, title('zeta-tsint'); figure, title('xchi_alpha2(3)');4. 將ansys計(jì)算得到的位移數(shù)據(jù)轉(zhuǎn)化為加速度load s1s2;a1a2=s1s2;

18、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-1 v1v2(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(2,1)-s1s2(1,2)/dt;v1

19、v2(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-1 a1a2(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)=(v1v2(2,2)-v1v2(1,2)/d

20、t;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 求特征值z(mì)eros 0數(shù)組abs 絕對(duì)值數(shù)組n=round(t200/10); 四舍五入disp 顯示h,l=size(s1s2);給出數(shù)組的行數(shù)和列數(shù)y=randn(ndzb,1);隨機(jī)數(shù)yp=mean(y)；平均數(shù)hp=load(sprintf('jhp%d.txt',jjj);將文本文件中的數(shù)據(jù)給矩陣hp，否則將生成一個(gè)矩陣和文件

21、同名。6. 計(jì)時(shí)tic toc7. ansys的一些命令/OUTPUT,apinggu0929,txt, 將輸出窗口內(nèi)容輸入到文件fini推出當(dāng)前模塊，到最高層/cle 重新開始計(jì)算/config,nres,30000/config,fsplit,750 將結(jié)果文件分成幾份，每份大小為標(biāo)準(zhǔn)的千分之750/INPUT,'0STEPG1','TXT','','', 0 *dim,li,array,nnn,1*dim,bo,TABLE,nn,1,1,time,wave*tread,bo,suiji,txt, *CFOPEN,'y

22、22z28',' ',' ' *VWRITE,weiyiy(1,22,1), weiyiy(1,23,1) ,weiyiy(1,24,1), weiyiy(1,25,1), weiyiy(1,26,1),weiyiy(1,27,1),weiyiy(1,28,1)(7E18.8E3)*CFOPEN,'li',' ',' ' *VWRITE,li(1,1,1), ,(E18.8E3) 8. 自相關(guān)函數(shù)function s0=zixg(s)s0=0; for j=s s0=s0+j2;end 9. 從ANSYS文

23、件中提取有用數(shù)據(jù)clc;clear; % 注意，本文件只適用于數(shù)據(jù)列數(shù)是3，只有英文注釋，且注釋最后以'gal'為結(jié)尾的情況ch1='k17' % 要打開的文件名ch2=strcat(ch1,'.asc') % 要打開的文件的擴(kuò)展名是ASCch3=strcat(ch1,'.mat') % 要存入的文件的擴(kuò)展名是MATfp=fopen(ch2,'r'); % 打開待處理文件buff=fscanf(fp,'%s',1) % 讀取第一個(gè)字符串hang=3000000; % 數(shù)據(jù)區(qū)總行數(shù)lie=3; uii

24、=0; % 希望得到的數(shù)據(jù)列數(shù)ss=zeros(hang,lie); % 數(shù)據(jù)要存入的矩陣 while(1) if strcmp(buff,'gal')=1 % 第一次判斷特征字符串 （判斷是否即將離開注釋區(qū)） 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

25、',1);ss(1,3)=fbuff; ui=2; while(1) buff=fscanf(fp,'%s',1); if (length(buff)=0) uii=1 break end if (abs(buff)<40) uii=1 break end ss(ui,1)=str2num(buff); fbuff=fscanf(fp,'%f',1);ss(ui,2)=fbuff; fbuff=fscanf(fp,'%f',1);ss(ui,3)=fbuff; ui=ui+1; end if uii=1 break end end

26、if uii=1 break end end if uii=1 break end buff=fscanf(fp,'%s',1) if length(buff)=0 break endendss(ui:hang,:)=; save(ch3,'ss');10. 傅里葉變換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) ax

27、is tight;xlabel time(s);ylabel m/s2;mn=1;y=y(1:8192);% plot(y)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(&#

28、39;f(Hz)')axis tight;a1,a2=max(zm);w(a2) 11. 梁截面的顯示開關(guān)是什么？/ESHAPE,1.0 Utility Menu>PlotCtrls>Style>Size and Shape12. 提取結(jié)構(gòu)反力*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在我以往動(dòng)力計(jì)算中，缺省的步驟是1000，如果希望做更多的步數(shù)，就可利用這個(gè)命令，其中20000是我在某次計(jì)算中動(dòng)力計(jì)算的步數(shù)。fini/cle/config,nres,20000/config,fsplit,750fini：退出當(dāng)前