汽车悬架双质量系统的传递特性仿真研究.docxVIP

Frequency (rad/s) Frequency (rad/s) 实例13汽车悬架双质量系统的传递特性仿真研究 根据汽车理论可知悬架双质量系统微分方程为m2z2 c(z2 z1) k(z2 z 根据汽车理论可知悬架双质量系统微分方程为 m2z2 c(z2 z1) k(z2 z1) 0 gzi c(乙 Z2) k(zi Z2) kt(zi q) (13-1) 0 对式(13-3 )和式(13-4 )进行拉氏变换并整理，可得 2 z2(m2s cs k) z1(cs k) 2 z1(m1s cs k kt) z2(cs k) qkt (13-2) 由上式可得Z2和?之间的传递函数为 G(S) cs k m2 s2 cs k (13-3) 2另 A cs k, A 2 另 A cs k, A2 m2s cs k kt, 2 gs cs k kt,将(13-3)代入(13-2),可得 到Z1与路面激励q的传递函数为 G2(s) G2(s) A kt 2 A3A2 A1 A?kt N (13-4) 13.1车身位移Z2与路面激励位移q的传递函数 现在可分析车轮与车身双质量系统的传递函数。由式（ 13-3）（ 13-4 ）相乘可以得到车 身位移Z2与路面激励位移q的传递函数为 Z2 Ajkt G（s） - G（s）G2（s） - - （ 13-5） q N MATLAB多项式计算函数由于传递函数分母为高阶多项式相乘，计算量比较大，因此可利用 求出分母N的系数。具体程序如下： MATLAB多项式计算函数 m2=317.5; m仁 45.4; k=22000; kt=192e3; c=1.5e3; a仁[c k]; a2=[m2 c k]; a3=[m1 c k+kt]; n仁con v(a3,a2); N仁poly2sym( n1); n 2=c on v(a1,a1); N2=poly2sym( n2); nn=N1-N2; pretty( nn); a1=[c*kt k*kt]; den=[28829/2,544350288000000,4224000000]; sys=tf(a1,de n); w=0.1:.1:100; figure⑴ [h,w1]=freqs(a1,de n, w); freqs(a1,de n, w); 运行可得到传递函数表达式以及传递函数的频率响应特性图: Tran sfer function: 2.88e008 s + 4.224e009 1.441e004 sA4 + 544350 sA3 + 6.894e007 sA2 + 2.88e008 s + 4.224e009 由图（ 由图（ 13-1 ）中的频率响应幅值（ Magnitude ）可以清楚的看到，在 0.1 到 100rad/s 的频率围， 有两个明显的共振峰值， 由汽车理论可知， 这两个共振峰值是由汽车双质量系统 两个主频率 1和 2决定的。下面进一步具体计算汽车双质量系统两个主频率 1 和 2。具体 程序如下： 200Frequency (rad/s)-20010—?、1000-100010i10102 200 Frequency (rad/s) -200 10 — ?、 100 0 -100 0 10 i 10 10 2 m2=317.5; m1=45.4; k=22000; kt=192e3; c=1.5e3; a1=[c*kt k*kt]; den=[28829/2,544350288000000,4224000000]; sys=tf(a1,den); w=0.1:.1:100; [h,w1]=freqs(a1,den,w); freqs(a1,den,w); title(z2/z0) mag=abs(h); margin(sys); j=0; for i=1:length(w)-1; if (mag(i+1)-mag(i))0(mag(i)-mag(i-1))0; j=j+1; magmax(1)=mag(i); wm(1)=w(i); elseif(mag(i+1)-mag(i))0abs(mag(i)-mag(i-1))0.0001; j=j+1; magmax(2)=mag(i); wm(2)=w(i); end end wm1=wm(1)/2/pi; wm2=wm(2)/2/pi; magm1=magmax(1); magm2=magmax(2); i1=wm1/f0; i2=wm2/f0; fprintf（ 汽车双质量悬架系统主频率 1（低频） w1=%3.3f HZ \n,wm1） 汽车双质量悬架系统主频率 1（低频） w1=1.210 HZ fprintf（