钟艳春，杨庆俊，包钢

(哈尔滨工业大学 机电工程学院，黑龙江 哈尔滨，150001)

气动隔振系统承载力大，有效行程大，隔振效果优异，广泛应用于仪器平台、光学平台、精密加工与检测、舰船动力设备隔振、车辆悬挂等的振动隔离[1-5]。由于非线性微分动力系统的复杂性，空气弹簧的弹性恢复力非线性对隔振系统特性的影响没有同时得到有效研究。为了简化计算过程，在大量的实际操作中，空气弹簧或被近似处理成定刚度弹簧，或忽略空气弹簧变形过程中有效面积以及其他参数改变后，由空气弹簧内气体的状态确定弹簧的状态[6]。近年来，有学者开始气动隔振系统的非线性研究，Heertjes等[7]研究了气体压缩的非线性对传递率和导纳的影响。尹万建等[8]研究了空气弹簧在单频正弦激励情况下的主共振，并认为激励幅值与非线性刚度是影响空气弹簧悬架非线性特性的主要因素，但是，没有研究系统多频激励下的非线性特性。在非线性方程方面，周一峰等[9]建立了主动隔振体的非线性动力学方程，对求解强非线性自治系统的能量迭代方法加以改进，将其用于求解强非线性非自治系统，得到了主动隔振系统周期运动响应的解析表达式和振幅-频率关系曲线。Nayfeh等[10]研究了平方和立方非线性系统在2个简谐激励共同作用下的组合共振和联合共振情况。在此，本文作者首先从基本的实验数据出发，得到空气弹簧相对载荷随高度的曲线，然后用三次多项式拟合该试验曲线，得到弹簧恢复力与高度之间的关系式。将恢复力代入载荷运动动力学方程得到系统在3个简谐激励共同作用下的非线性平衡方程，利用多尺度法[11-13]对该非线性振动方程进行求解，分析各参数变化对系统幅频特性的影响，以便为指导气动隔振系统的设计提供理论依据。

1 空气弹簧非线性模型

本文选用简单的单腔气弹簧对其进行非线性特性分析。气动隔振系统简化模型如图1所示。图1中：M为弹簧上的质量；K为空气弹簧总刚度，c为考虑空气弹簧自身阻尼在内的总阻尼；xb为基座振动，x为负载振动。

图1 气动隔振系统模型简图Fig.1 Schematic diagram of pneumatic vibration isolation system

1.1 空气弹簧载荷

选用 Firestone 1M1A-1型空气弹簧作为研究对象。空气弹簧在初始充气压力为0.7 MPa时，不同高度下的绝对载荷如图2所示。

图2 Firestone 1M1A-1型空气弹簧绝对载荷随高度变化曲线Fig.2 Absolute load curve with height of Firestone 1M1A-1 air spring

从图2可以看出：与一般弹簧的特性不同，空气弹簧的绝对弹力总是正值。但相对平衡位置，弹力可以有正有负：在平衡位置以上为负，在平衡位置以下为正。相对于平衡位置75 mm高度时的相对载荷如图3所示。

1.2 空气弹簧恢复力拟合

文献[6]指出仅用三次多项式就能非常精确地拟合试验曲线，采用太高的幂次一方面不会显著提高精度，另一方面会使拟合刚度曲线光滑性变差。因此，本文采用三次多项式对该曲线拟合。空气弹簧相对平衡位置的弹性恢复力F可以表示为：

式中：F为空气弹簧相对平衡位置的弹性恢复力；x为载荷振动；xb为基座振动；k1为空气弹簧的线性刚度；k2为二次非线性系数，k3为三次非线性系数。

以实验所得的空气弹簧相对弹力为原始数据，根据式(1)对其进行三次多项式拟合，拟合曲线如图3所示。通过曲线拟合得到的三次多项式系数分别为k1=20.5 N/m，k2=3 063 N/m2，k3=212 797 N/m3。从图3可以看出：拟合曲线和实验数据重合性好。

图3 Firestone 1M1A-1型空气弹簧相对载荷多项式拟合曲线Fig.3 Relative load polynomial fitting curve of Firestone 1M1A-1 air spring

2 多频激励下气动隔振系统非线性模型

考虑对于多频激励进行分析时，3个频率激励具有很好典型性，能很好地反映系统多频激励下的规律性，并且分析也简单明了，因此，本文选用3个激励频率进行分析。假设基座激励为：

式中：F1，F2和F3依次为3个外部激励的幅值；Ω1，Ω2和 Ω3依次为 3个外部激励的频率；θ1，θ2和 θ3依次为3个外部激励的初相角。

在其激励下，系统的动力学方程为：

式中：M为载荷质量；c为系统阻尼系数。

设y=x-xb，则方程(3)可以写为：

构造长度尺度 L=F1和时间尺度得到无量纲变量对式(4)进行无量纲化，则得到：式中

运用多尺度法对式(5)进行分析，设式(5)的近似解为：

式中：T0=τ，T2=ετ分别为不同的时间尺度变量。

因此，关于τ的导数变成关于T0和T1的偏导数的展开式，即有

式中： Dn=∂ / ( ∂ Tn)。

将式(6)～(8)代入式(5)，并展开比较ε同次幂的系数，得到：

方程(9)的通解可以写为：

将方程(11)代入方程(10)得到：

从式(12)可以看出：在 3个频率激励下，系统会出现超谐共振 2ωn≈1，3ωn≈1，亚谐共振 ωn≈2，ωn≈3，组合共振 ωm± ωn≈1，ωm± ωn≈2，ωm± ωn≈3，2ωm±ωn≈1，2ωm±ωn≈2，2(ωm±ωn)≈1，ωm± ωn± ωx≈1等。

下面讨论 ω1±ω2+ ω3≈1的情况。引进调谐参数σ，令

将其代入(ω1+ω2+ω3)T0可以写为：

由式(12)式可得，若

令式(14)和(15)α ˙ =γ˙= 0 得其稳定解。得到下列方程组：

由式(18)可得频率响应为：

联合前面各表达式，并将其代入方程(19)得到频率响应和各初始参数之间的关系为：

3 系统参数对频率响应曲线的影响

方程(20)即为系统的频率响应关系式，下面分别固定其他参数而变化其中1个参数，得到各主要参数对系统幅频曲线的影响。图4～7所示分别为各主要参数变化对系统幅频特性影响的曲线。

图4 激励幅值变化时的组合共振幅频曲线Fig.4 Frequency-response curves in differentexcitation amplitudes

图 4 中，k1=20.5 N/m，k2= -3 063 N/m2，k3=21 279 N/m3，ω1=0.18，ω2=0.34，ω3=0.49，M=280 kg，c=140 N·s/m。从图4可以看出：激励幅值对非线性响应曲线有明显影响；随着激励幅值的增大，组合共振的响应振幅也增大。并且随着激励幅值的增大，幅频曲线的骨架曲线的弯曲程度变得更大，说明激励幅值越大，系统非线性特性越明显。

图5中，其余各参数取值如下：k1=20.5 N/m，k2=-3 063 N/m2，F1=F2=F3=0.015 m，ω1=0.18，ω2=0.34，ω3=0.49，M=280 kg，c=140 N·s/m。从图5可以看出：三次非线性系数直接影响幅频响应曲线的非线性特性。当k3＞0时，幅频曲线的骨架曲线向右弯曲，当k3＜0时，幅频曲线的骨架曲线向左弯曲；k3的绝对值越大，曲线的弯曲程度越大；组合共振响应曲线的幅值也随k3的绝对值的增大而增大。

图5 三次非线性系数变化时的组合共振幅频曲线Fig.5 Frequency-response curves in different cube nonlinear parameters

图6 中，其余各参数取值如下：k1=20.5 N/m，k2=-3 063 N/m2，k3=21 279 N/m3，F1=F2=F3=0.015 m，从图6可知，当载荷变小时，这种组合共振的幅频曲线弯曲程度变大，非线性现象变得更明显。

图7中，其余各参数取值如下：k1=20.5 N/m，k2=-3 063 N/m2，k3=21 279 N/m3，F1=F2=F3=0.015 m，ω1=0.18，ω2=0.34，ω3=0.49，M=280 kg。从图7可以看出：随着阻尼系数的增大，组合共振响应曲线的幅值变小。而幅频曲线的骨架曲线不受阻尼的影响，系统的共振域也基本不变。

图6 载荷变化时的组合共振幅频曲线Fig.6 Frequency-response curves in different load masses

图7 阻尼系数变化时的组合共振幅频曲线Fig.7 Frequency-response curves in different damping coefficients

4 结论

(1) 根据实验得到的空气弹簧动载荷曲线，运用三次非线性多项式对其进行拟合，得到一定初始压力下空气弹簧相对平衡位置的非线性恢复力与高度的关系式，为对系统非线性求解奠定了基础。

(2) 在3个频率激励下系统会出现超谐共振、亚谐共振、组合共振等情况。对 ω1+ω2+ω3≈1的组合共振进行分析，得到了频率响应方程。

(3) 激励幅值越大，系统非线性现象越明显；三次非线性系数对系统非线性影响大，当k3＞0时，幅频曲线的骨架曲线向右弯曲；，当k3＜0时，幅频曲线的骨架曲线向左弯曲；载荷变小时，组合共振幅频曲线弯曲程度变大；阻尼变小时，组合共振响应曲线的幅值变大，系统的共振域不受阻尼变化的影响。

[1] Kato T, Kawashima K, Sawamoto K, et al. Active control of a pneumatic isolation table using model following control and a pressure differentiator[J]. Precision Engineering, 2007, 31(3):269-275.

[2] CHEN Ping-chang, Shih Ming-chang. Modeling and robust active control of a pneumatic vibration isolator[J]. Journal of Vibration and Control, 2007, 13(11): 1553-1571.

[3] 徐伟, 何琳, 吕志强, 等. 船舶主机气囊隔振系统动态特性分析[J]. 振动与冲击, 2007, 26(7): 122-124.XU Wei, HE Lin, LÜ Zhi-qiang, et al. Dynamic analysis of an air spring mounting system for marine main engine[J]. Journal of Vibration and Shock, 2007, 26(7): 122-124.

[4] Lee J H, Kim K J. A method of transmissibility design for dual-chamber pneumatic vibration isolator[J]. Journal of Sound and Vibration, 2009, 323(1): 67-92.

[5] Shimozawa K, TohtakE T. An air spring model with non-linear damping for vertical motion[J]. Quarterly Report of RTRI, 2008,49(4): 209-214.

[6] 陈燎, 周孔亢, 李仲兴. 空气弹簧动态特性拟合及空气悬架变刚度计算分析[J]. 机械工程学报, 2010, 46(4): 93-98.CHEN Liao, ZHOU Kong-kang, LI Zhong-xing. Dynamic characteristics fitting of air springs and numerical analysis of air suspensions with variant stiffness[J]. Chinese Journal of Mechanical Engineering, 2010, 46(4): 93-98.

[7] Heertjes M, van de Wouw B. Nonlinear dynamics and control of a pneumatic vibration isolator[J]. Journal of Vibration and Acoustics.2006, 128(4): 439-448.

[8] 尹万建, 韩鹰, 杨绍普. 空气弹簧悬架系统在强迫振动下的动力学分析[J]. 中国公路学报, 2006, 19(3): 117-121.YIN Wan-jian, HAN Ying, YANG Shao-pu. Dynamics analysis of air spring suspension system under forced vibration[J]. China Journal of Highway and Transport, 2006, 19(3): 117-121.

[9] 周一峰, 唐进元, 何旭辉. 强非线性主动隔振系统的运动响应及传递率[J]. 中南大学学报: 自然科学版, 2005, 36(3):496-500.ZHOU Yi-feng, TANG Jing-yuan, HE Xu-hui. Response and transmissibility of strong nonlinear active isolation system[J].Journal of Central South University: Science and Technology,2005, 36(3): 496-500.

[10] Nayfeh A H. Combination tones in the response of single degree of freedom systems with quadratic and cubic non-linearities[J].Journal of Sound and Vibration, 1984, 92(3): 379-386.

[11] Nayfeh A H, Mook D T. Nonlinear oscillations[M]. New York:Wiley, 1979: 161-224.

[12] 刘延柱, 陈立群. 非线性振动[M]. 北京: 高等教育出版社,2001: 83-95.LIU Yan-zhu, CHEN Li-qun. Nonlinear vibration[M]. Beijing:High Education Press, 2001: 83-95.

[13] 杨绍普, 申永军. 滞后非线性系统的分岔与奇异性[M]. 北京:科学出版社, 2003: 35-40.YANG Shao-pu, SHEN Yong-jun. Bifurcation and singularity of non-linear systems with hysteresis[M]. Beijing: Science Press,2003: 35-40.