单自由度系统振动.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*t).*sin(wd*t+phi); fprintf(A=%.3g\n,A); fprintf(phi=%.3g\n,phi); elseif z==1 a1=x0; a2=v0+wn*x0; fprintf(a1=%.3g\n,a1); fprintf(a2=%.3g\n,a2); x=(a1+a2*t).*exp(-wn*t); else a1=(-v0+(-z+sqrt(z^2-1))*wn*x0*wn*x0)/2/wn/sqrt(z^2-1); a2=(v0+(z+sqrt(z^2-1))*wn*x0)/2/wn/sqrt(z^2-1); fprintf(a1=%.3g\n,a1); fprint(a2=%.3g\n,a2); x=exp(-wn*t).*(a1*exp(wn*sqrt(z^2-1)*t) +a2*exp(wn*sqrt(z^2- 1)*t)); end plot(t,x),grid xlabel(时间(s)) ylabel(位移) title(位移相对时间的关系) 叠加原理 第1章 单自由度系统振动 1.1 概述 1.2 单自由度系统振动 1.3 等效质量与等效刚度 1.4 隔振原理 1.5 等效粘性阻尼 1.6 非谐周期激励的响应 1.7 单位脉冲的响应 1.8 任意激励的响应 目的： 了解单自由度系统振动的特点、解的形式，并借助MATLAB工具求解单自由度系统振动，精确计算振动系统的响应。 掌握单自由度系统等效处理的方法，理解隔振基本原理。 了解不同的外加激励下，单自由度系统振动响应的求解方法。 1.1 概 述 实际振动系统往往是很复杂的。为了便于分析计算，在满足工程要求的前提下，可对复杂振动系统进行简化。 处理得到的最简单的系统就是单自由度振动系统。 单自由度系统振动是工程中常见的物理现象。 梁： 车辆通过桥梁时引起的振动 起重机：悬吊钢丝绳——重物的振动 用单个变量的 微分方程即可描 述的振动系统 齿轮减（变）速器：轴——齿轮的扭转振动 汽轮机、发电机：转子不平衡引起的扭转振动 讨论：请列举其他单自由度振动的实际例子（工程中、家用设备、相关课程中涉及到的）！ 1.2 单自由度系统振动 1.无阻尼自由振动 2.有阻尼自由振动 3.有阻尼受迫振动 4.MATLAB数值仿真 单自由度振动系统类型： 1.无阻尼自由振动解 单自由度无阻尼系统振动方程为： 方程的通解为： 令 若振动的初始条件： 则其解为： 其中 则固有频率为 振幅为 相位为 【例】起重机以速度 使重物 下降时，突然紧急刹车， ，钢丝绳弹性模量 求此时提升机构所受的最大的力。 已知 ，钢丝绳截面积 钢丝绳长度为 【解】紧急刹车时，钢丝绳突然停止，但此时重物具有速度 从制动的瞬间开始吊在绳上作自由振动。 初始速度为 最大位移 初始位移为 钢丝绳最大的拉伸量为 为钢丝绳刚度，由材料力学可知 钢丝绳最大拉力为 动拉力与静拉力之比为动力放大系数 结论：当紧急制动时，起重机钢丝绳中的动拉力是正常提升时的3.2134倍 2.有阻尼自由振动解 其通解为： 称为阻尼率。 ξ 1时，称为弱