第1章自由度系统振动.pptVIP

有阻尼自由振动响应计算程序 function y=VTB1 (m,c,k,x0,v0,tf) %m为质量;c为阻尼;k为刚度;x0为初始位移;v0为初始速度;tf为仿真时间 %VTB1(zeta,w,x0,v0,tf)绘出单自由度有阻尼自由振动系统的响应图 %zeta为阻尼系数;ωn为固有频率 %程序中z为阻尼系数;A为振动幅度;phi为初相位 clc %该循环确定输入方式是VTB1(m,c,k,x0,v0,tf),还是%VTB1(zeta,w,x0,v0,tf) if nargin==5 z=m;wn=c;tf=v0;v0=x0;x0=k;m=1;c=2*z*w;k=w^2; end wn=sqrt(k/m); %固有频率 z=c/2/m/wn; wd=wn*sqrt(1-z^2); fprintf(固有频率为%.3g.rad/s.\n,wn); fprintf(阻尼系数%.3g.\n,z); fprintf(有阻尼的固有频率%.3g.\n,wd); t=0:tf/1000:tf; if z1 A=sqrt(((v0+z*wn*x0)^2+(x0*wd)^2)/wd^2); phi=atan2(x0*wd,v0+z*wn*x0); x=A*exp(-z*wn