动力学第三章
摘要:第2章functionVTB2(m,c,k,x0,v0,tf,w,f0)%单自由度系统的谐迫振动clcwn=sqrt(k/m);z=c/2/m/wn;lan=w/wnwd=wn*sqrt(1-z^2);A=sqrt(((v0+z*wn*x0)^2+(x0*wd)^2)/wd^2);t=0:tf/1000:tf;phi1=atan2(x0*wd,v0+z*wn*x0);phi2=atan2(2*z*lan,1-lan^2);B=wn^2*f0/k/sqrt((wn^2-w^2)^2+(2*z*wn*w)^2);x=A*exp(-z*wn*t).*sin(sqrt(1-z^2)*wn*t+phi1)+B*sin(w*t-phi2);plot(t,x),gridxlabel('时间(s)')ylabel('位移')title('位移与时间的关系')functionVTB1(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-z^2);fprintf('固有频率为%.s.\n',wd);fprintf('阻尼系数为%.3g.\n',z);fprintf('有阻尼的固有频率为%.s
温馨提示：当前文档最多只能预览
5 页，若文档总页数超出了
5 页，请下载原文档以浏览全部内容。
本文档由 匿名用户 于 2020-11-11 04:17:06上传分享
第5章-发热电动力
