《动力学第三章.doc》由会员分享，可在线阅读，更多相关《动力学第三章.doc（13页珍藏版）》请在金锄头文库上搜索。
1、第2章function VTB2(m,c,k,x0,v0,tf,w,f0)%单自由度系统的谐迫振动clcwn=sqrt(k/m);z=c/2/m/wn;lan=w/wnwd=wn*sqrt(1-z2);A=sqrt(v0+z*wn*x0)2+(x0*wd)2)/wd2);t=0:tf/1000:tf;phi1=atan2(x0*wd,v0+z*wn*x0);phi2=atan2(2*z*lan,1-lan2);B=wn2*f0/k/sqrt(wn2-w2)2+(2*z*wn*w)2);x=A*exp(-z*wn*t).*sin(sqrt(1-z2)*wn*t+phi1)+B*sin(w*t-p
2、hi2);plot(t,x),gridxlabel(时间(s)ylabel(位移)title(位移与时间的关系)function VTB1(m,c,k,x0,v0,tf)%VTB1用来计算单自由度有阻尼自由振动系统的响应%VTB1绘出单自由度有阻尼自由振动系统的响应图%m为质量;c为阻尼;k为刚度;x0为初始位移;v0为初始速度;tf为仿真时间%VTB1(zeta,w,x0,v0,tf)绘出单自由度有阻尼自由振动系统的响应图%程序中z为阻尼系数;wn为固有频率n;A为振动幅度;phi为初相位clcwn=sqrt(k/m);z=c/2/m/wn;wd=wn*sqrt(1-z2);fprintf(
3、固有频率为%.3g.rad/s.n,wd);fprintf(阻尼系数为%.3g.n,z);fprintf(有阻尼的固有频率为%.3g.rad/s.n,wd);t=0:tf/1000:tf;if z1.e-6 pp0=pp; B=D*A; pp=1.0/B(3); A=B/B(3); end f=sqrt(pp)/2/pi %固有频率单位转换为Hz fprintf(fid1,%20.5f,A); %输入主振型数据 fprintf(fid2,%20.5f,f); %输入固有频率数据 D=D-A*A*M/(A*M*A*pp);endfid1=fopen(A-1,rt); %打开主振型文件A=fsca
4、nf(fid1,%f,3,3) %主振型写成矩阵fid2=fopen(B-1,rt); %打开固有频率文件f=fscanf(fid2,%f,3,1) %固有频率写成矩阵t=1:3;h1=figure(numbertitle,off,name,0,pos,50 200 420 420);bar(t,f(t,1),xlabel(频率阶级次),ylabel(Hz),title(固有频率),hold on,gridh1=figure(numbertitle,off,name,1,pos,50 200 420 420);bar(t,A(t,1),xlabel(自由度（质量块）),ylabel(振型向量)
5、,title(1阶主振型),hold on,gridpause(0.1)h1=figure(numbertitle,off,name,2,pos,50 200 420 420);bar(t,A(t,2),xlabel(自由度（质量块）),ylabel(振型向量),title(2阶主振型),hold on,gridpause(0.1)h1=figure(numbertitle,off,name,3,pos,50 200 420 420);bar(t,A(t,3),xlabel(自由度（质量块）),ylabel(振型向量),title(3阶主振型),hold on,gridpause(0.1)en
6、d%chuandijuzhen.m; %传递矩阵方法求固有频率clear all,clear closeJ1=1;J2=1;J3=2;k2=1100000;k3=1200000;k4=100000;fid=fopen(chuandi,wt); %建立（打开）速度文件M1L=0;for WN=0:0.01:2000 shita1R=1;M1R=-WN2*J1; shita2R=shita1R+1/k2*M1R;M2R=shita1R*(-WN2*J2)+(1+(-WN2*J2)/k2)*M1R; shita3R=shita2R+1/k3*M2R;M3R=shita2R*(-WN2*J3)+(1+
7、(-WN2*J3)/k3)*M2R; shita4R=shita3R+1/k4*M3R; if abs(shita4R)0.005 WN %搜索到的固有频率(rad/s),并显示 shita=shita1R;shita2R;shita3R;shita4R;%搜索到振型,并显示 bar(shita),xlabel(对应的质量块),ylabel(振型向量) pause(1.0) end fprintf(fid,%30.15f,shita4R);endfid=fopen(chuandi,rt);x=fscanf(fid,%f,1,200001);t=1:200001;plot(0.01*t,x);g
8、rid,xlabel(频率rad/s),ylabel(第四个质量块的转角(rad/s),title(用传递矩阵法求固有频率) function cdjz2%chuandijuzhen.m; %传递矩阵方法求固有频率clear all,clear closeJ1=0.5;J2=1;k2=100000;k3=100000;fid=fopen(chuandi3,wt); %建立（打开）速度文件M1L=0;for WN=0:0.01:2000 shita1R=1;M1R=-WN2*J1; shita2R=shita1R+1/k2*M1R;M2R=shita1R*(-WN2*J2)+(1+(-WN2*J2)/k2)*M1R; shita3R=shita2R+1/k3*M2R; if abs(shita3R)1.e-6 pp0=pp; B=D*A; pp=1.0/B(3); A=B/B(3); end f=sqrt(pp) fprintf(fid1,%20.5f,A); %输入主振型数据 fprintf(fid2,%20.5f,
