1、精選優質文檔-傾情為你奉上車輛懸架技術汽車平順性分析姓名：學號:班號：汽車的平順性就是保持汽車在行駛過程中乘員所處的振動和沖擊環境具有一定舒適度的性能，對于貨車還包括保持貨物完好的性能。平順性主要根據乘員主觀感受的舒適性來評價，同時也輔助以客觀評價的方法。平順性是現代高速汽車的主要性能之一。平順性分析的框圖表示：評價指標：Ø 車體質心垂直振動加速度(ACC)， 反映乘員乘坐舒適性和車體振動環境Ø 懸架動行程(SWS)，即車輪相對于車體垂直跳動的動位移，反映了車輪撞擊限位器的概率Ø 車輪相對動載(DTL)，即相對于靜平衡位置時車輪載荷的變化，衡量車輪抓地能力，反映了

2、高速車輛的行駛安全性基于以上，本文主要包括三方面內容：Ø 系統激勵分析，路面不平度，路面功率譜密度，路面等級Ø 車輛振動系統模型的搭建Ø 根據系統響應的平順性分析(一).路面激勵分析（程序見附錄一）通常把路面相對基準平面的高度q，沿道路走向長度l的變化q(l),稱為路面縱斷面曲線或路面不平度函數。路面不平度的均方根值q=10.7976mm在matlab中利用Gq,n=pwelch(q(l),0,NFFT,Ns),其中Gq為路面的功率譜密度，n=1/為路面波的空間頻率，NFFT為快速傅立葉變換的采樣點數，Ns=1/0.05=20m-1為采樣的空間頻率.畫出標準路面的

3、功率譜密度Gq=Gq(n0)*no2/n2,n0=0.1。圖中上面的直線取Gq(n0)=512mm2/m(-1),下面的直線取Gq(no)=128mm2/m(-1)。對比路面不平度系數的分級標準：路面等級 mm2/ m-1 mm下限幾何平均上限下限幾何平均上限ABCDEFGH832128512204881923276816642561024409616384655363212851220488192327682.695.3810.7721.5343.0686.13172.26344.523.817.6115.2330.4560.90121.81243.61487.225.3810.7721.5

4、343.0686.13172.26344.52689.04 該路面屬于C級路面。 路面不平度的快速傅立葉變換Q=fft(q)，求取幅值magnitude=abs(Q)在車速30Km/h時，路面波主要集中在低頻段。 （二）.simulink模型搭建利用系統動力學方程M*X+C*X+K*X=f,其中M為質量矩陣M=Mb00Mw,C為阻尼矩陣C=Cs-Cs-CsCs,K為剛度矩陣K=Ks-Ks-KsKs+Kt,f為激勵f=0Kt*Xr(三).系統響應的平順性分析（程序見附錄二）Ø 車體質心垂直振動加速度(ACC)車體質心加速度均方根值：sigma_acc_Xb =1.1007m/s2對質心

5、加速度進行功率譜密度分析：Ø 懸架動行程(SWS)懸架動行程均方根值sigma_SWS =0.0103m對懸架動行程進行功率譜密度分析：Ø 車輪相對動載(DTL)車輪動載的均方根值sigma_DTL =549.7985N對車輪動載進行功率譜密度分析附錄一clear;namein='F:Personal Documentmatlabroadroad1l.m'q=dlmread(namein);average=mean(q); q=q-average; l=0:0.05:(numel(q)-1)*0.05;plot(l,q);grid on;sigma=std(

6、q);xlabel('測量長度lm');ylabel('路面不平度度q(l)mm');title('路面不平度');Ns=1/0.05;Q=fft(q);magnitude=abs(Q);phase=angle(Q);n=(0:numel(q)-1)*Ns/numel(q);f=30/3.6*n;plot(n,phase);figure;plot(f,magnitude);xlabel('頻率/Hz');ylabel('幅值/mm');title('幅頻特性');Gq,n=pwelch(q,0,25

7、6,Ns);loglog(n, Gq) ;hold on;xlabel('空間頻率1/m');ylabel('功率譜密度mm2/(1/m)');title('路面的功率譜密度');n=0.01:0.1:10;GqBup=8*0.12./n.2;GqBlow=32*0.12./n.2loglog(n, GqBup);grid on;loglog(n, GqBlow);grid on;附錄二clear;Mb=317.5;Mw=45.4;Cs=1520;Ks=22000;Kt=;u=30; %車速30Km/hnamein='F:Personal

8、 Documentmatlabroadroad1l.m'Xr=dlmread(namein)/1000;t=0:0.05*3.6/u:0.05*3.6/u*(numel(Xr)-1);t=t'M=Mb,0;0,Mw;C=Cs,-Cs;-Cs,Cs;K=Ks,-Ks;-Ks,Ks+Kt;sim('one_fourth');acc_Xb=Acceleration.signals.values(:,1);plot(t,acc_Xb);xlabel('時間s');ylabel('車體質心加速度m/s2');title('車體質心加速度');grid on;average=mean(acc_Xb);sigma_acc_Xb=std(acc_Xb);Fs=1/(0.05*3.6/u); Gq_acc_Xb,f=pwelch(acc_Xb,0,256,Fs);fig