单自由度系统振动

单自由度系统振动

ID:27076805

大小:555.51 KB

页数:24页

时间:2018-12-01

单自由度系统振动_第1页
单自由度系统振动_第2页
单自由度系统振动_第3页
单自由度系统振动_第4页
单自由度系统振动_第5页
资源描述:

《单自由度系统振动》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、第1章单自由度系统振动1.1概述1.2单自由度系统振动1.3等效质量与等效刚度1.4隔振原理1.5等效粘性阻尼1.6非谐周期激励的响应1.7单位脉冲的响应1.8任意激励的响应目的:了解单自由度系统振动的特点、解的形式,并借助MATLAB工具求解单自由度系统振动,精确计算振动系统的响应。掌握单自由度系统等效处理的方法,理解隔振基本原理。了解不同的外加激励下,单自由度系统振动响应的求解方法。1.1概述实际振动系统往往是很复杂的。为了便于分析计算,在满足工程要求的前提下,可对复杂振动系统进行简化。处理得到的最简单的系统就是单自由度振动系统。单自由度系统振动是工程中常见的物理现象。梁:车辆通过

2、桥梁时引起的振动起重机:悬吊钢丝绳——重物的振动用单个变量的微分方程即可描述的振动系统齿轮减(变)速器:轴——齿轮的扭转振动汽轮机、发电机:转子不平衡引起的扭转振动讨论:请列举其他单自由度振动的实际例子(工程中、家用设备、相关课程中涉及到的)!1.2单自由度系统振动1.无阻尼自由振动2.有阻尼自由振动3.有阻尼受迫振动4.MATLAB数值仿真单自由度振动系统类型:1.无阻尼自由振动解单自由度无阻尼系统振动方程为:方程的通解为:令若振动的初始条件:则其解为:其中则固有频率为振幅为相位为【例】起重机以速度使重物下降时,突然紧急刹车,,钢丝绳弹性模量求此时提升机构所受的最大的力。已知,钢丝绳

3、截面积钢丝绳长度为【解】紧急刹车时,钢丝绳突然停止,但此时重物具有速度从制动的瞬间开始吊在绳上作自由振动。初始速度为最大位移初始位移为钢丝绳最大的拉伸量为为钢丝绳刚度,由材料力学可知钢丝绳最大拉力为动拉力与静拉力之比为动力放大系数结论:当紧急制动时,起重机钢丝绳中的动拉力是正常提升时的3.2134倍2.有阻尼自由振动解其通解为:称为阻尼率。ξ<1时,称为弱阻尼状态通解为:ξ>1时,称为强阻尼状态通解为:ξ=1时,称为临界阻尼状态通解为:单自由度有阻尼系统振动方程为:注意希腊字母Ξ[ksi];ζ[zta]3.有阻尼受迫振动解瞬态响应稳态响应其解为瞬态与稳态响应两部分的叠加:瞬态响应是一个

4、有阻尼的谐振。振动频率为系统固有频率,振幅与初相位角决定于初始条件,振幅按指数规律衰减,振动持续时间决定于系统的阻尼比。稳态响应也是一个简谐振动,其频率比等于激励力的频率。振动方程为时,为谐迫振动。其解为频率比稳态响应的振幅瞬态响应的振幅相位受迫振动中阻尼比和频率比的影响分析振幅:频率比λ=1附近,受迫振动的幅值出现峰值——共振,且随着阻尼率ξ的增加,共振幅值逐渐减小;当ξ=0.7~0.8时振幅曲线平滑衰减小。相位:频率比λ=1,受迫振动的相位滞后为90°,随着阻尼率ξ的增大,相位曲线趋近线性,阻尼率ξ=0.7~0.8时,相位曲线最接近直线。工程中:尽量取阻尼率ξ=0.7~0.8注意希

5、腊字母ξ(ksi)4.MATLAB数值仿真MATLAB是MatrixLaboratory的缩写,是一种直译式的语言,易学(相比C语言)特点:强大的数值运算功能丰富的工具箱数学计算数字信号处理自动控制动态分析数据处理2D与3D绘图功能有阻尼自由振动响应计算程序functiony=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%该循环确定输入方式是

6、VTB1(m,c,k,x0,v0,tf),还是%VTB1(zeta,w,x0,v0,tf)ifnargin==5z=m;wn=c;tf=v0;v0=x0;x0=k;m=1;c=2*z*w;k=w^2;endwn=sqrt(k/m);%固有频率z=c/2/m/wn;wd=wn*sqrt(1-z^2);fprintf('固有频率为%.3g.rad/s.',wn);fprintf('阻尼系数%.3g.',z);fprintf('有阻尼的固有频率%.3g.',wd);t=0:tf/1000:tf;ifz<1A=sqrt(((v0+z*wn*x0)^2+(x0*wd)^2)/wd^2

7、);phi=atan2(x0*wd,v0+z*wn*x0);x=A*exp(-z*wn*t).*sin(wd*t+phi);fprintf('A=%.3g',A);fprintf('phi=%.3g',phi);elseifz==1a1=x0;a2=v0+wn*x0;fprintf('a1=%.3g',a1);fprintf('a2=%.3g',a2);x=(a1+a2*t).*exp(-wn*t);elsea1=(-v0+(-

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。