《Matlab在汽车振动分析(修改)课件》由会员分享，可在线阅读，更多相关《Matlab在汽车振动分析(修改)课件（61页珍藏版）》请在金锄头文库上搜索。
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时间t激励
3、函数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.5,0;001;K=5-2
6、,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二自由度系统汽车二自由系统在路面不平激励下的强迫振动。本例是在路面不平激励下的强迫振动，教材中的3.3.3节.书中有大
7、量的理论推导，得到的公式和计算结果，在Matlab里变得很简单传递函数求解函数TFESTIMATEH2,ft=TFESTIMATE(q,x2(1:n),window,noverlap,nfft,fs);Txy=tfestimate(x,y)findsa transfer functionestimateTxygiveninputsignalvectorxandoutputsignalvectory.Txy,F=tfestimate(x,y,window,noverlap,nfft,fs)returnsTxyasafunctionoffrequencyandavectorFoffrequenci
8、esatwhichtfestimateestimatesthetransferfunction.fsisthesamplingfrequencyinHz.FisthesamesizeasTxy,soplot(f,Txy)plotsthetransferfunctionestimateversusproperlyscaledfrequency.Forrealsignals,therangeofFis0,fs/2whennfftisevenand0,fs/2)whennfftisodd.Forcomplexsignals,therangeofFis0,fs).计算结果（滤波前）计算结果（滤波）频率
9、ft与传递函数H1的关系loglog Log-logscaleplotloglog(ft,abs(H1)plot(ft,abs(H1)plot(log(ft),log(abs(H1)plot(ft,H1)Warning:ImaginarypartsofcomplexXand/orYargumentsignored复数的绝对值abs(H1(80)ans=0.1192H1(80)ans=-0.0434-0.1110iplot(ft,H1)plot(ft,abs(H1)X+Yi讨论：得到了频率与传递函数的关系，有什么用？学习软件人一能之已百之，人十能之已千之。果能此道矣，虽柔比强。摘自中庸博学之，审
10、问之，慎思之，明辨之，笃行之。有弗学，学之弗能，弗措也；有弗问，问之弗知，弗措也；有弗思，思之弗得，弗措也；有弗辨，辨之弗明，弗措也。有弗行，行之弗笃，弗措也。人一能之已百之，人十能之已千之。果能此道矣，虽柔比强。摘自中庸TransferFunctionApproachtoModelingDynamicsystemTransferfunction=Thetransferfuctionofalinear,time-invariantdifferential-equationsystemisdefinedastheratiooftheLaplacetransfermoftheoutput(resp
11、onsefuntion)totheLaplaceinput(drivingfunction)undertheassumptionthatallinitialconditionsarezeroTheequationofmotionforthesystemisTakingthelaplacetransformofbothsidesofthisequationandassumingthatallinitialconditionsarezeroyieldsThelaplacetransformofg(t)givesthetransferfuncion.Therefore,thetransferfunc
12、tionandimpulse-responsefuctionofalinear,time-invariantsystemcontainsthesameinformationaboutthesystemdynamics.Impulse-ResponsefunctionSincethelaplacetransformoftheunit-impulsefunctionisunity,orX(s)=1TheinverseLaplacetransformoftheoutputequationyieldstheimpulseresponseofthesystem.Howtoobtainthesystemr
13、esponseanalytically?Partial-fractionexpansionwithMatlabThecomandr,p,k=residue(num,den)Findstheresidues,poles,anddirecttermsofapartialfractionexpansionoftheradioofthetwopolynomialsB(s)andA(s)num=25;den=1,4,25;r,p,k=residue(num,den)r=0-2.7277i0+2.7277ip=-2.0000+4.5826i-2.0000-4.5826ik=查表laplacetransfo
14、rmspairsWherey(t)ismeasuredinmetersandtinseconds.Thisequationisananalyticalsolutiontotheproblem.TransientresponseanalysiswithmatlabThissectionpresentstheMatlabapproachtoobtainingsystemresponsewhentheinputsarethetime-dormaininputssuchasthestep,impuse,andrampfunctions.Matlabusessystorepresentsuchasyst
15、emSys=tf(num,den)Forexamplenum=25;den=1,4,25;sys=tf(num,den)Transferfunction:25-s2+4s+25Producetheunit-stepresponseofthesystemnum=25;den=1,4,25;sys=tf(num,den)step(sys)ImpulseresponseTheimpulseresponseofamechanicalsystemcanbeobservedwhenthesystemissubjectedtoaverylargeforceforaveryshorttime,forinsta
16、nce,whenthemassofaspring-mass-dashpotsystemishitbyahammerorabullet.num=25;den=1,4,25;sys=tf(num,den)impulse(sys)该函数由傅立叶方程求解方法该函数由拉普拉斯函数求解 （2）应用SIMULINK，对单自由度系统在简谐激励下的振动响应进行仿真。其中(3)(3)应用应用MATLABMATLAB，分析，分析单缸活塞式发动机在惯性力作用下的振动单缸活塞式发动机在惯性力作用下的振动振动微分方程求解振动问题的关键在于建立描述系统的运动方程振动描述的一些关键量在振动响应函数的数学描述振幅固有频率初相位临界阻尼振动方程的描述自由振动强迫振动振动方程的求解方法已知条件：建立起了振动系统的微分方程求解：振动响应函数x(t)杜哈美积分法傅氏积分法拉氏变换法中文名称：中文名称：固有频率英文名称：英文名称：naturalfrequency：由系统本身的质量和刚度所决定的频率。n自由度系统一般有n个固有频率，按频率的高低排列，最低的为第一阶固有频率。有阻尼的线性系统的自由振动频率称为“阻尼固有频率”。物体作
