Matlab在汽车振动分析(修改)

Matlab在振动分析中的应用刘迪辉2011-10-20大家学了游泳理论，现在我们借助MATLAB软件，来练习一下游泳！实际问题：客车的振动分析•客车样车路试过程中却出现了令人意想不到的一系列振动问题,主要表现为:(1)汽车起动时发动机抖动厉害;(2)当车速在40km/h左右时,整车有共振现象;(3)当车速在85km/h左右时,整车有明显振动;(4)当车速超过118km/h时,驾驶区及方向盘有强烈振感。•由于上述振动的存在,一方面大大降低了该车驾乘的舒适性和运行中的安全性;另一方面,造成一些主要总成件(如发动机、变速器、后桥等)的早期损坏;同时,也使得汽车上很多结构件出现疲劳断裂,从而进一步加剧了整车或局部振动。•选自王卫鸿《YBL6850C24aH》型客车振动问题及解决方案，客车技术与研究，2005.5SimulinkDemo•Thisdemodescribesasimplifiedhalf-carmodelthatincludesanindependentfrontandrearverticalsuspension.振动问题•多自由度•二自由度•单自由度•实际问题•（1）理论方法•（2）Matlab（实现理论算法）•(3)有限元方法Ansys,Abaqus,Natran等•（4）试验方法难易•建立力学模型、微分方程•求解微分方程，得到响应特性振动方程)(tfkxxcxm时间t响应函数x(t)质量m刚度k阻尼c时间t激励函数f(t)（1）已知激励函数和响应函数，求系统固有特性（2）已知固有特性，求在一定激励条件下的响应函数汽车悬架单自由度分析•例2.15质量m=2450kg的汽车，悬架总的刚度为160000N/m,减振器阻尼系数为7135.6Ns/m,求该车辆受到100kg的简谐加载时的，车身的上下运动方程.0()()()sinmxtcxtkxtFt•简谐激励wtFtfsin)(005101520253035404550-1-0.8-0.6-0.4-0.200.20.40.60.81tf(t)f(t)=F0sin(wt)首先得设定参数F0,w,和时间向量t,求每个时间的f(t)理论公式0()()()sinmxtcxtkxtFtwtqxpxpxsin22mcp2该函数由普通微分方程求解方法其中mkp2mFq0提问：为什么要如此参数化？方便求解和定义联系起来mkpmpc2固有频率系统阻尼加载瞬间的瞬态响应瞬态响应（衰减振动）)sin(1tpAexpt其中21'pp00.511.522.533.544.55-0.6-0.4-0.200.20.40.60.8tx(t)x(t)瞬态响应00x200x100x00.511.522.533.544.55-0.6-0.4-0.200.20.40.60.811.2tx(t)x(t)00.511.522.533.544.55-0.8-0.6-0.4-0.200.20.40.60.811.2tx(t)x(t)00.511.522.533.544.55-0.8-0.6-0.4-0.200.20.40.6tx(t)x(t)200x)sin(1tpAexpt大家还可研究一些瞬态振动基本概念减幅系数，衰减系数，对数衰减率，衰减振动的周期临界阻尼00.511.522.533.544.5500.0020.0040.0060.0080.010.0120.0140.0160.0180.02tx(t)x(t)临界阻尼（汽车不振动）msNcc/39600稳态响应)sin(2wtXx稳态振动位移响应的圆频率等于激励的圆频率05101520253035404550-8-6-4-202468x10-3tx(t)x(t)05101520253035404550-1000-800-600-400-20002004006008001000tf(t)f(t)=F0sin(wt)激振函数响应函数稳态响应05101520253035404550-1000-800-600-400-20002004006008001000tf(t)f(t)=F0sin(wt)05101520253035404550-8-6-4-202468x10-3tx(t)x(t))sin(2wtXx例2.15中同时受到瞬态与稳态响应•瞬态响应•稳态响应02468101214161820-8-6-4-202468x10-3tx(t)x(t)02468101214161820-4-202468x10-3tx(t)x(t)02468101214161820-8-6-4-202468x10-3tx(t)x(t)稳态+瞬态这里主要是加载频率和固有频率相差比较大018.0pw•最大振幅6mm•瞬态响应00.511.522.533.544.55-4-202468x10-3tx(t)x(t)00.511.522.533.544.55-0.02-0.015-0.01-0.00500.0050.010.0150.02tx(t)x(t)00.511.522.533.544.55-0.02-0.015-0.01-0.00500.0050.010.0150.02tx(t)x(t)当固有频率和加载频率相同时，总的响应明显增大。共振！•稳态响应1pw•最大振幅15mm稳态+瞬态例题11.10()()()sinmxtcxtkxtFt•为什么叫幅频特性？•为什么叫相频特性？图中说明了什么信息？例题11.12220)2()1(1XX程序采用公式212tanacr三自由度下振型求解•振型•modeshape;modeofvibration机械系统某一给定振动模态的振型，指在某一固有频率下，由中性面或中性轴上的点偏离其平衡位置的最大位移值所描述的图形。•振型是指体系的一种固有的特性。它与固有频率相对应，即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。•振型与体系实际的振动形态不一定相同。•振型对应于频率而言，一个固有频率对应于一个振型。按照频率从低到高的排列，来说第一振型，第二振型等等。此处的振型就是指在该固有频率下结构的振动形态，频率越高则振动周期越小。在实验中，我们就是通过用一定的频率对结构进行激振，观测相应点的位移状况，当观测点的位移达到最大时，此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状，而是各阶振型相叠加的结果。例题11.2求固有频率和振型•clc;clearall;closeall••M=[200;0,1.5,0;001];•K=[5-2,0;-2,3,-1;0,-1,1];••eigfreq(M,K)•Matlab/Simulink在单自由度强迫振动分析中的应用MatlabSimulink•在Matlab命令窗口，敲击simulink，可打开该模块。•Simulink是图型化语言。把公式变成流程图。减少理论推导的工作量。•Simulink可分析系统在一个激励函数下的响应函数等例1uxx2dxxdxdxx•积分器integrater•系数gain•积分器后默认为x例2单自由度0()()()sinmxtcxtkxtFtF0=980;m=2450;c=7135;k=160000;0()()()sinmxtcxtkxtFt•激励函数•响应函数例3二自由度系统01221222xxkxxcxm0122122111qxkxxkxxcxm•汽车二自由系统在路面不平激励下的强迫振动。本例是在路面不平激励下的强迫振动，教材中的3.3.3节.书中有大量的理论推导，得到的公式和计算结果，在Matlab里变得很简单01221222xxkxxcxm0122122111qxkxxkxxcxm传递函数求解函数TFESTIMATE•[H2,ft]=TFESTIMATE(q,x2(1:n),window,noverlap,nfft,fs);•Txy=tfestimate(x,y)findsatransferfunctionestimateTxygiveninputsignalvectorxandoutputsignalvectory.•[Txy,F]=tfestimate(x,y,window,noverlap,nfft,fs)returnsTxyasafunctionoffrequencyandavectorFoffrequenciesatwhichtfestimateestimatesthetransferfunction.fsisthesamplingfrequencyinHz.FisthesamesizeasTxy,soplot(f,Txy)plotsthetransferfunctionestimateversusproperlyscaledfrequency.Forrealsignals,therangeofFis[0,fs/2]whennfftisevenand[0,fs/2)whennfftisodd.Forcomplexsignals,therangeofFis[0,fs).•计算结果（滤波前）•计算结果（滤波）频率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.1192•H1(80)•ans=•-0.0434-0.1110i•plot(ft,H1)•plot(ft,abs(H1))•X+Yi•讨论：得到了频率与传递函数的关系，有什么用？学习软件人一能之已百之，人十能之已千之。果能此道矣，虽柔比强。——摘自《中庸》•博学之，审问之，慎思之，明辨之，笃行之。•有弗学，学之弗能，弗措也；•有弗问，问之弗知，弗措也；•有弗思，思之弗得，弗措也；•有弗辨，辨之弗明，弗措也。•有弗行，行之弗笃，弗措也。•人一能之已百之，人十能之已千之。果能此道矣，虽柔比强。•摘自《中庸》TransferFunctionApproachtoModelingDynamicsystem•Transferfunction=)()()(sXsYinputoutputsGThetransferfuctionofalinear,time-invariantdifferential-equationsystemisdefinedastheratiooftheLaplacetransfermoftheoutput(responsefuntion)totheLaplaceinput(drivingfunction)undertheassumptionthatallinitialconditionsarezero•Theequationofmotionforthesystemis)(tfkxxcxmTakingthelaplacetransformofbothsidesofthisequationandassumingthatallinitialconditionsarezeroyieldskcsmssFsXsG21)()()(•Thelaplacetransformofg(t)givesthetransferfuncion.Therefore,thetransferfunctionandimpulse-responsefuctionofalinear,time-invariantsystemcontainsthesa