资源描述:
《Matlab在汽车振动分析修改》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、Matlab在振动分析中的应用刘迪辉2011-10-20大家学了游泳理论,现在我们借助MATLAB软件,来练习一下游泳!实际问题:客车的振动分析客车样车路试过程中却出现了令人意想不到的一系列振动问题,主要表现为:(1)汽车起动时发动机抖动厉害;(2)当车速在40km/h左右时,整车有共振现象;(3)当车速在85km/h左右时,整车有明显振动;(4)当车速超过118km/h时,驾驶区及方向盘有强烈振感。由于上述振动的存在,一方面大大降低了该车驾乘的舒适性和运行中的安全性;另一方面,造成一些主要总成件(如发动机、变速器、后桥等)的早期损坏;同时,也使得汽车上很多结构件出现疲劳断裂,从而进一步
2、加剧了整车或局部振动。选自王卫鸿《YBL6850C24aH》型客车振动问题及解决方案,客车技术与研究,2005.5SimulinkDemoThisdemodescribesasimplifiedhalf-carmodelthatincludesanindependentfrontandrearverticalsuspension.振动问题多自由度二自由度单自由度实际问题(1)理论方法(2)Matlab(实现理论算法)(3)有限元方法Ansys,Abaqus,Natran等(4)试验方法难易建立力学模型、微分方程求解微分方程,得到响应特性振动方程时间t响应函数x(t)质量m刚度k阻尼c时
3、间t激励函数f(t)(1)已知激励函数和响应函数,求系统固有特性(2)已知固有特性,求在一定激励条件下的响应函数汽车悬架单自由度分析例2.15质量m=2450kg的汽车,悬架总的刚度为160000N/m,减振器阻尼系数为7135.6Ns/m,求该车辆受到100kg的简谐加载时的,车身的上下运动方程.简谐激励首先得设定参数F0,w,和时间向量t,求每个时间的f(t)理论公式该函数由普通微分方程求解方法其中提问:为什么要如此参数化?方便求解和定义联系起来固有频率系统阻尼加载瞬间的瞬态响应瞬态响应(衰减振动)其中瞬态响应大家还可研究一些瞬态振动基本概念减幅系数,衰减系数,对数衰减率,衰减振动的
4、周期临界阻尼临界阻尼(汽车不振动)稳态响应稳态振动位移响应的圆频率等于激励的圆频率激振函数响应函数稳态响应例2.15中同时受到瞬态与稳态响应瞬态响应稳态响应稳态+瞬态这里主要是加载频率和固有频率相差比较大最大振幅6mm瞬态响应当固有频率和加载频率相同时,总的响应明显增大。共振!稳态响应最大振幅15mm稳态+瞬态例题11.1为什么叫幅频特性?为什么叫相频特性?图中说明了什么信息?例题11.1程序采用公式三自由度下振型求解振型modeshape;modeofvibration机械系统某一给定振动模态的振型,指在某一固有频率下,由中性面或中性轴上的点偏离其平衡位置的最大位移值所描述的图形。振型
5、是指体系的一种固有的特性。它与固有频率相对应,即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言,一个固有频率对应于一个振型。按照频率从低到高的排列,来说第一振型,第二振型等等。此处的振型就是指在该固有频率下结构的振动形态,频率越高则振动周期越小。在实验中,我们就是通过用一定的频率对结构进行激振,观测相应点的位移状况,当观测点的位移达到最大时,此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状,而是各阶振型相叠加的结果。例题11.2求固有频率和振型clc;clearall;closeallM=[200;0,1
6、.5,0;001];K=[5-2,0;-2,3,-1;0,-1,1];eigfreq(M,K)Matlab/Simulink在单自由度强迫振动分析中的应用MatlabSimulink在Matlab命令窗口,敲击simulink,可打开该模块。Simulink是图型化语言。把公式变成流程图。减少理论推导的工作量。Simulink可分析系统在一个激励函数下的响应函数等例1积分器integrater系数gain积分器后默认为x例2单自由度F0=980;m=2450;c=7135;k=160000;激励函数响应函数例3二自由度系统汽车二自由系统在路面不平激励下的强迫振动。本例是在路面不平激励下的
7、强迫振动,教材中的3.3.3节.书中有大量的理论推导,得到的公式和计算结果,在Matlab里变得很简单传递函数求解函数TFESTIMATE[H2,ft]=TFESTIMATE(q,x2(1:n),window,noverlap,nfft,fs);Txy = tfestimate(x,y)findsatransferfunctionestimateTxygiveninputsignalvectorxandoutputsignalvec