研究生阶段的一些经验总结2091231.doc

研究生阶段的一些经验总结2091231目录1. 将数组中的数据存储在文本文件中12.画图的一些命令格式应用举例（有点乱，但比较有用）13. 在Matlab图形中怎样输入特殊符号?44. 将ansys计算得到的位移数据转化为加速度45. 常用函数56. 计时57. ansys的一些命令58. 自相关函数69. 从ANSYS文件中提取有用数据610. 傅里叶变换711. 梁截面的显示开关是什么？812. 提取结构反力813. /config,nres,20000814. *dim,weiyiy,array,nnn,nm915. 循环结构和选择结构916. 命令FLST917. 已知结构的前二阶固有频率f和阻尼比，求阻尼矩阵的系数918. 实验 阻尼比和频率计算101. 将数组中的数据存储在文本文件中命令：save y2.txt y1 -ASCII TABS说明：将y1这个数组中的数据以字符形式存储在文件y2.txt中2.画图的一些命令格式应用举例（有点乱，但比较有用）figure(2);plot(yyy);xlabel(结点序号);ylabel(振型改变率（%）);title(第1阶振型改变率)axis tight;grid on;加注释框plot(pl,p1,-*,pl,p2,-d,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 onplot(hh,pinlv(2,:),-ok,LineWidth,2,MarkerSize,12)hhh=legend(测点1,测点2);xlim(2 8);ylim(0 1); xlabel(工况,FontSize,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 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,-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(层数,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,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(,.% String,$pi to tau_0gg$gg,.% Position,.5 .5,.% FontSize,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/s2);在一个大的画面上，画2行2列小图，将数据充满数轴，横坐标单位是time(s)，纵坐标单位是m/(m/s2)figure(xv);产生第xv个图plot(tyyt,yy,-*);画曲线legend(j=1,j=2)；注释框ylabel(strcat(sprintf(NR%dj(0)变化率,xv),(%);上句涉及到两个函数，sprintf生成一个字符串，strcat将字符串连接起来title(sprintf(方法1的NR%dj(0)变化率,xv);xlabel(计算次序);axis tight;set(gca,XTick,tyyt);设置横坐标的刻度plot(tyyt,yy(:,1),-*k, markersize,10,LineWidth,2);在画曲线的同时，给出了线的宽度，颜色，符号的大小xlabel(损伤位置宽度改变量/mm,FontSize,18);横坐标的刻度值的标度的字符大小是18ylabel(strcat(sprintf(itNRrm_%d_j(tau_0)变化率,xv),/%),FontSize,18);在图上纵坐标上写一些希腊字母set(gca,XTick,tyyt,FontSize,18);xlim(0 12);x轴的取值范围是0到12hhh=legend(完好状态,弹性摸量减少3%,弹性摸量减少6%,弹性摸量减少9%,弹性摸量减少12%,2);将注释框的句柄给hhhaxis(-0.62 2.1 0 3)；给定了横坐标和纵坐标的取值范围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文件格式来显示字符。使用bf，it，rm表示黑体，% 斜体，正体字符，特别注意大括号 的用法。3. 在Matlab图形中怎样输入特殊符号?答: 上标用 (指数)下标用 _(下划线)调用方式为:任意字符, _任意字符, 注意的使用!希腊字母等特殊字符用 加拼音 如: alpha, beta , gamma theta Theta Gamma delta Delta xi Xi elta epsilong zeta miu nu tau lamda Lamda pi Pi sigma Sigma phi Phi psi Psi chi ommiga Ommiga geq 不等于 neq gg正负 pm 左箭头 leftarrow 右箭头 rightarrow上箭头 uparrow体会以下两例:figure, title(zeta-tsint); figure, title(xchi_alpha2(3);4. 将ansys计算得到的位移数据转化为加速度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-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;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-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)/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. 常用函数eig 求特征值zeros 0数组abs 绝对值数组n=round(t200/10); 四舍五入disp 显示h,l=size(s1s2);给出数组的行数和列数y=randn(ndzb,1);随机数yp=mean(y)；平均数hp=load(sprintf(jhp%d.txt,jjj);将文本文件中的数据给矩阵hp，否则将生成一个矩阵和文件同名。6. 计时tic toc7. ansys的一些命令/OUTPUT,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), weiyiy(1,26,1),weiyiy(1,27,1),weiyiy(1,28,1)(7E18.8E3)*CFOPEN,li, , *VWRITE,li(1,1,1), ,(E18.8E3) 8. 自相关函数function s0=zixg(s)s0=0; for j=s s0=s0+j2;end 9. 从ANSYS文件中提取有用数据clc;clear; % 注意，本文件只适用于数据列数是3，只有英文注释，且注释最后以gal为结尾的情况ch1=k17; % 要打开的文件名ch2=strcat(ch1,.asc) % 要打开的文件的扩展名是ASCch3=strcat(ch1,.mat) % 要存入的文件的扩展名是MATfp=fopen(ch2,r); % 打开待处理文件buff=fscanf(fp,%s,1) % 读取第一个字符串hang=3000000; % 数据区总行数lie=3; uii=0; % 希望得到的数据列数ss=zeros(hang,lie); % 数据要存入的矩阵 while(1) if strcmp(buff,gal)=1 % 第一次判断特征字符串 （判断是否即将离开注释区） buff=fscanf(fp,%s,1); if strcmp(buff,Write)=1 % 第二次判断特征字符串 （判断是否已经离开注释区） ss(1,1)=str2num(buff); fbuff=fscanf(fp,%f,1);ss(1,2)=fbuff; fbuff=fscanf(fp,%f,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)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在我以往动力计算中，缺省的步骤是1000，如果希望做更多的步数，就可利用这个命令，其中20000是我在某次计算中动力计算的步数。fini/cle/config,nres,20000/config,fsplit,750fini：退出当前模块/cle：进行一个新的计算/config,fsplit,750：由于结果文件过大，ansys将停止计算，因此，可以设置多个结果文件，没有结果文件是原来4G的千分之750。14. *dim,weiyiy,array,nnn,nm定义一个数组，其类型是array，其名称是weiyiy，其位数维数是nnn行，nm列*dim,bo,TABLE,nn,1,1,time,wave定义一个数组，其类型是TABLE，其名称是bo，其位数维数是nn行，1列，1页。可以插值计算，其行标名称是time，其列标名称是wave*tr