汽车减振器设计与特性仿真
上QQ阅读APP看书,第一时间看更新

2.2 单质量车身振动及特性

2.2.1 单质量车身振动微分方程

图2-6所示为分析车身振动的单质量系统模型,该模型由车身质量m2和弹簧刚度为k、减振器阻尼系数为c的悬架组成,q为路面不平度函数,它是以沿路前进方向的坐标x为参数的随机过程。

取车身垂直位移坐标z的原点在静力平衡位置,可得到系统运动的微分方程为

978-7-111-37673-6-Chapter02-7.jpg

978-7-111-37673-6-Chapter02-8.jpg

这是单自由度系统随机基础位移激励问题,系统的响应是由两部分组成的,即由系统自由振动响应和强迫振动响应叠加组合而成,也就是说,系统振动响应是振动微分方程的齐次方程的解和非齐次方程特解之和。

978-7-111-37673-6-Chapter02-9.jpg

图2-6 单质量车身振动模型

2.2.2 单质量系统的自由振动响应

978-7-111-37673-6-Chapter02-10.jpg978-7-111-37673-6-Chapter02-11.jpg,当激励q=0时,则由式(2-8)可得单质量系统的自由振动微分方程为

978-7-111-37673-6-Chapter02-12.jpg

式中,p为系统的固有圆频率;而阻尼对系统的影响,取决于np的比值ξ,称ξ为阻尼比,即

978-7-111-37673-6-Chapter02-13.jpg

汽车悬架系统阻尼比ξ通常在0.25左右,属于小阻尼。因此,振动系统的齐次微分方程的解,也就是车身自由衰减的振动响应为

978-7-111-37673-6-Chapter02-14.jpg

式中,A是由初始条件所决定的常数,978-7-111-37673-6-Chapter02-15.jpgφ是由初始条件决定的初相角,978-7-111-37673-6-Chapter02-16.jpg

由式(2-11)可知,有阻尼的自由衰减振动,车身质量m2以有阻尼的固有频率978-7-111-37673-6-Chapter02-17.jpg振动,而振幅却按Ae-ξpt规律衰减,如图2-7所示。

由上述可知,阻尼对自由振动有如下两方面的影响。

978-7-111-37673-6-Chapter02-18.jpg

图2-7 自由衰减振动曲线

1.阻尼使固有频率降低

如果无阻尼自由振动系统原固有频率为978-7-111-37673-6-Chapter02-19.jpg,则弱阻尼悬架系统的固有频率p′

978-7-111-37673-6-Chapter02-20.jpg

因此,可知小阻尼悬架系统的固有频率p′随阻尼比ξ的增大而降低。

由于汽车悬架系统的阻尼比ξ约为0.25,因此,阻尼是悬架系统的固有频率仅下降了3%左右,可以忽略不计,所以,工程上小阻尼振动系统的固有圆频率p′,可以近似地认为等于无阻尼振动系统的固有圆频率p,即p′p。因此,车身振动的固有圆频率和固有频率,分别为

978-7-111-37673-6-Chapter02-21.jpg

2.阻尼决定振幅的衰减程度

设相邻两振幅分别为AiAi+1(图2-7),它们的比值η称为减幅系数

978-7-111-37673-6-Chapter02-22.jpg

式中,n为衰减系数。n越大表示阻尼越大,振幅衰减也就越大。令978-7-111-37673-6-Chapter02-23.jpg为对数衰减率,因此可得

978-7-111-37673-6-Chapter02-24.jpg

所以,由式(2-16)可得振动系统的阻尼比为

978-7-111-37673-6-Chapter02-25.jpg

ξ≪1,则由式(2-16)得,lnη≈2πξ,因此由式(2-17)可知,小阻尼车辆悬架系统的阻尼比可近似表示为

978-7-111-37673-6-Chapter02-26.jpg

2.2.3 单质量系统在简谐激振力下的响应

由于阻尼会使自由振动逐渐衰减,最后达到完全停止。工程上一些能持续下去的振动必定有外加能源,以弥补阻尼所消耗的能量,使系统的振动不会衰减。因此,工程上多采用强迫振动的响应。

若简谐激振力ft)=F0sinωt,则根据牛顿第二定律,可得单质量振动系统在简谐激振力ft)=F0sinωt作用下的振动微分方程,可表示为

978-7-111-37673-6-Chapter02-27.jpg

式(2-19)可写为

978-7-111-37673-6-Chapter02-28.jpg

式中,978-7-111-37673-6-Chapter02-29.jpg978-7-111-37673-6-Chapter02-30.jpg978-7-111-37673-6-Chapter02-31.jpg

振动微分方程式(2-20)的解包括两部分:齐次方程的通解z1和方程的特解z2,即

z=z1+z2

由2.2.2节可知,在弱阻尼(ξ<1)的情况下,有阻尼自由振动齐次方程的解z1

978-7-111-37673-6-Chapter02-32.jpg

式(2-21)代表的是一种衰减振动,只在振动开始的一段时间内才有意义,故为瞬态振动。在一般情况下实际工程意义不大,可以不予考虑。

振动微分方程式(2-20)的特解z2,代表系统在简谐激振下所产生的强迫振动,它是一种持续的等幅振动,故为稳态振动。设特解z2

z2=Zsin(ωt-ψ) (2-22)

式中,Z为振动响应的幅值;ω为激振力圆频率,也是振动响应的圆频率;ψ为响应滞后于激励的相位差。

又因为

978-7-111-37673-6-Chapter02-33.jpg

978-7-111-37673-6-Chapter02-34.jpg

将式(2-22)~式(2-24)代入微分方程(2-20),可得

-ω2Zsin(ωt-ψ)+2ξpωZcos(ωt-ψ)+p2Zsin(ωt-ψ)=Fsin ωt(2-25)

利用三角函数关系得

Fsin ωt=Fsin[(ωt-ψ)+ψ]=Fcos ψsin(ωt-ψ)+Fsin ψcos(ωt-ψ)(2-26)

比较式(2-26)和式(2-25),由于对任何瞬时t都成立,故sin(ωt-ψ)和cos(ωt-ψ)前的系数必须分别相等,即

p2-ω2Z=Fcos ψ

2ξpωZ=Fsin ψ

因此,可得

978-7-111-37673-6-Chapter02-35.jpg

式中,978-7-111-37673-6-Chapter02-36.jpg为频率比;978-7-111-37673-6-Chapter02-37.jpg,为系统的最大静位移。

因此,强迫振动的稳态解为

978-7-111-37673-6-Chapter02-38.jpg

由上述强迫振动解可见:在简谐激振力作用下,强迫振动响应为也简谐振动,其频率与激振频率ω相同,但相位角滞后ψ,这是由于阻尼存在的关系。振幅Z与相位差ψ都只与系统固有特性及激振力的性质有关,而与初始条件无关。

由式(2-29)可得ZZ0之比ββ称为放大因子,为

978-7-111-37673-6-Chapter02-39.jpg

放大因子β代表稳态振幅X与激振力幅F0作用于弹簧上的静位移Z0之比。β值不仅随λ而变,而且还随ξ值而变。

在不同的阻尼比ξ的情况下,放大因子β与频率比λ的关系以及相位角ψλ的关系,如图2-8和图2-9所示。其中,图2-8所示为幅频响应曲线,而图2-9所示为相频响应曲线。

978-7-111-37673-6-Chapter02-40.jpg

图2-8 幅频响应曲线

978-7-111-37673-6-Chapter02-41.jpg

图2-9 相频响应曲线

1)当λ≪1,即激振频率ω远小于系统的固有频率p时,无论阻尼大小如何,β接近于1,即振幅近似等于激振力幅值F0作用下的静变形X0。故在低频区内,振幅X主要由弹簧刚度控制。此时,相位差ψ≈0,即位移与激振力接近于同相位。

2)当λ≫1,即激振频率ω远大于系统的固有频率p时,β趋近于0。因为激振力方向改变太快,振动物体由于惯性来不及跟随,几乎停止不动。故在高频区内,振幅X主要决定于系统的惯性。这一特性正是隔振和惯性传感器的理论依据。相位差ψ≈π,即在高频范围内位移与激振力接近于反相位。

3)当λ≈1时,即ω接近p,振幅Z急剧增加,β趋向βmax,这种现象称为共振。严格地讲,βmax发生在978-7-111-37673-6-Chapter02-42.jpg处,但通常ξ2≪1,故ω=p时系统发生共振。由式(2-18)可以看出,振幅Z达到最大值978-7-111-37673-6-Chapter02-43.jpg时,由式(2-28)得978-7-111-37673-6-Chapter02-44.jpg

可见在共振时,振幅最大值Zmax与阻尼比ξ的值有关,ξ越小,则Zmax将越大;在ξ→0时,Zmax可达到无穷大。但共振时的ψ值与阻尼比ξ的值无关,不论ξ为何值,共振时的ψ总是978-7-111-37673-6-Chapter02-45.jpg,这是共振的一个重要特征。从分析幅频响应与相频响应所引出的共振现象,是传统的共振试验法测定系统固有频率的理论基础。

2.2.4 单质量系统在单位谐波函数激励下的响应

单位谐波函数激励为复数形式的单位幅值简谐激振力,即fct)=eiωt=cosωt+isinωt,则单质量系统的振动微分方程为

978-7-111-37673-6-Chapter02-46.jpg

单质量系统在单位谐波函数激励下的复数形式的响应为zct)。由于复数激振力和复数响应既是t的函数,又是ω的函数,故可令复数响应与复数激振力之比为Hω),即

978-7-111-37673-6-Chapter02-47.jpg

Hω)被称为频率响应函数。它是一个由系统特性参数所确定的,表示系统在单位幅值的简谐激振力fct)=eiωt作用下所产生的振幅。对于简谐激励,若知道系统的频率响应函数,便可由式(2-32)可求得振动系统的输出响应,即

zct)=Hωfct)(2-33)

根据式(2-31),可得单位简谐激振力作用下的响应为

978-7-111-37673-6-Chapter02-48.jpg

将上述三式代入式(2-32),两边消去eiωt,即得频率响应函数为

978-7-111-37673-6-Chapter02-49.jpg

式中,λ为频率比,ξ为阻尼比。

频率响应函数的模为978-7-111-37673-6-Chapter02-50.jpg,称为幅频特性。

频率响应函数的相位差角为978-7-111-37673-6-Chapter02-51.jpg,称为相频特性。

将复数形式的简谐激振力F0eiωt代入式(2-33),则复数形式的响应为

zc=Hω)e-iψF0eiωt=F0Hω)eiωt-ψ) (2-35)

若实际激振力为正弦函数F0sin ωt,则实际响应取复数形式响应的虚部,得实际解为

978-7-111-37673-6-Chapter02-52.jpg

若实际激励为余弦函数F0cos ωt,则取复数形式响应的实部,得实际解为

978-7-111-37673-6-Chapter02-53.jpg

2.2.5 单质量系统振动响应的傅里叶积分法

激励函数ft)的傅里叶积分形式为

978-7-111-37673-6-Chapter02-54.jpg

式(2-37)右端的积分运算称为激励ft)的傅里叶变换,式(2-36)相应地称为傅里叶逆变换。由式(2-36)和式(2-37)所联系的两个量ft)和Fω)称为一个傅里叶变换对。

通常响应函数zt)可以用傅里叶积分式(2-36)表示为

978-7-111-37673-6-Chapter02-55.jpg

式中,

978-7-111-37673-6-Chapter02-56.jpg,是响应zt)的傅里叶变换。

可以把非周期函数看成是由无数个复振幅为978-7-111-37673-6-Chapter02-57.jpg的谐波分量所组成,于是,根据式(2-35)求出对应于每个谐波分量的响应后,再根据线性系统的叠加原理,就可求得系统的响应

978-7-111-37673-6-Chapter02-58.jpg

比较式(2-38)和式(2-39),得

Xω)=HωFω) (2-40)

它表示输出和输入傅里叶变换之比,等于频率响应函数Hω),简称频响函数。这与在简谐激振力作用下的输出与输入关系式相同。这说明频率响应函数能表示系统的动态特性。在简谐激振力的作用下,线性单质量系统的频率响应函数为

978-7-111-37673-6-Chapter02-59.jpg

它的模978-7-111-37673-6-Chapter02-60.jpg,它的虚部与实部之比为相位角978-7-111-37673-6-Chapter02-61.jpg,分别确定系统的幅频特性和相频特性,能全面反映系统的传递特性。

2.2.6 单质量车身在路面激励下的振动响应

对式(2-8),通常关心其稳态随机响应,它取决于路面不平度函数随机激励qx)和系统的频率响应特性函数Hω)。由上可知,系统频率响应函数Hωz-q为系统的振动响应z的傅里叶变换与激励q的傅里叶变换之比,即

978-7-111-37673-6-Chapter02-62.jpg

式中,Zω)为响应zt)的傅里叶变换;Qω)为激励qt)的傅里叶变换。

对式(2-8)进行傅里叶变换,可得单质量车身在路面激励下响应的频响函数为

978-7-111-37673-6-Chapter02-63.jpg

式中,978-7-111-37673-6-Chapter02-64.jpg为阻尼比;λ为频率比,978-7-111-37673-6-Chapter02-65.jpgω为路面激励的圆频率;978-7-111-37673-6-Chapter02-66.jpg为系统固有圆频率。

由式(2-43)得为单质量车身在路面激励下的幅频特性和相频特性。

幅频特性为

978-7-111-37673-6-Chapter02-67.jpg

相频特性为

978-7-111-37673-6-Chapter02-68.jpg

汽车在具有一定幅值的正弦波路面上行驶,即路面激励为

qt)=asin ωt

则单质量车身在路面激励下的响应为

978-7-111-37673-6-Chapter02-69.jpg

路面激励qt)=asin ωt为正弦,所以系统的实际响应为

978-7-111-37673-6-Chapter02-70.jpg

式中,978-7-111-37673-6-Chapter02-71.jpg为幅值Z,即路面激励响应的幅值为

978-7-111-37673-6-Chapter02-72.jpg

如果路面激励以速度978-7-111-37673-6-Chapter02-73.jpg=bsin ωt来表达,用上面同样的推导方法可得

978-7-111-37673-6-Chapter02-74.jpg

若以加速度978-7-111-37673-6-Chapter02-75.jpg来表达,则有

978-7-111-37673-6-Chapter02-76.jpg

Hωz-q可以得到单质量系统的幅频特性曲线,如图2-10所示。

由频响函数式(2-44)和图2-10可知

1)当频率比λ=1时,系统出现共振,幅频特性达到最大,即共振时的幅值

978-7-111-37673-6-Chapter02-77.jpg

2)在低频段(0≤λ≤0.75),Hωz-q略大于1,不呈现明显的动态特性,阻尼比对低频段的影响不大。

978-7-111-37673-6-Chapter02-78.jpg

图2-10 单质量系统的幅频特性曲线

3)在共振段(0.75<λ978-7-111-37673-6-Chapter02-79.jpg),Hωz-q出现峰值,将输入激励放大,增大阻尼比ξ,可使共振峰值明显降低。

4)在高频段(978-7-111-37673-6-Chapter02-80.jpg),当978-7-111-37673-6-Chapter02-81.jpg时,

Hωz-q=1,系统响应与阻尼比ξ无关;当978-7-111-37673-6-Chapter02-82.jpg时,Hωz-q<1,对输入位移有衰减作用,且阻尼比减小对减振有利。