橡膠阻尼式扭轉減振器固有頻率計算與測試方法的研究

上官文斌，魏玉明，趙 旭，榮俊偉，王亞杰

(1.華南理工大學機械與汽車工程學院，廣東廣州510641；2.寧波拓普集團股份有限公司，浙江寧波315800)

橡膠阻尼式扭轉減振器固有頻率計算與測試方法的研究

上官文斌1，2，魏玉明1，趙 旭1，榮俊偉2，王亞杰2

(1.華南理工大學機械與汽車工程學院，廣東廣州510641；2.寧波拓普集團股份有限公司，浙江寧波315800)

設計了橡膠試片靜動態剪切特性實驗夾具，測試得到了橡膠試片在不同壓縮比下的靜動態剪切特性。建立了表征橡膠試片靜動態剪切特性的Kelvin-Voigt本構模型、Maxwell本構模型和分數導數本構模型，提出了由實驗測得的橡膠試片剪切特性曲線來識別各模型參數的方法。給出了基于這3種本構模型計算橡膠阻尼式扭轉減振器固有頻率的方法，計算了一扭轉減振器固有頻率，并與測試值進行對比。結果表明：基于Kelvin-Voigt模型計算得到的固有頻率不隨激振振幅的變化而變化，不能表征振幅相關性，而 Maxwell模型和分數導數模型計算固有頻率能表征振幅相關性。大振幅激勵時，利用 Maxwell模型計算的固有頻率與實測值的誤差較大；在各種激勵振幅下，利用分數導數模型可較準確計算出減振器的固有頻率。

扭轉減振器；橡膠；模型參數測試與辨識；固有頻率測試

引 言

減少發動機曲軸扭轉振動的常用方法是在曲軸前端安裝減振器［1］。橡膠阻尼式扭轉減振器因為結構簡單、成本較低，而廣泛應用于乘用車發動機的曲軸減振。橡膠阻尼式扭轉減振器的固有頻率是其最重要的性能參數之一，在扭轉減振器慣量環的慣量一定的情況下，其固有頻率與橡膠圈扭轉剛度、阻尼等有關。而橡膠圈扭轉剛度、阻尼跟橡膠圈壓縮比、激振頻率、激振振幅密切相關［2］。由于橡膠阻尼式扭轉減振器固有頻率受多種因素影響，以往的橡膠阻尼式扭轉減振器的設計開發是一個不斷試驗、修正的復雜過程。目前關于橡膠阻尼式扭轉減振器固有頻率計算方面的研究工作較少。因此，研究橡膠阻尼式扭轉減振器固有頻率計算方法具有較重要意義。

在橡膠動態黏彈性力學特性的研究中，廣泛應用Kelvin-Voigt，Maxwell、分數導數等力學模型［3-5］。Berg等［6］將Maxwell力學模型應用于軌道車輛橡膠隔振器動態特性計算，并取得了良好的效果。?stberg等［7］使用分數導數模型對一橡膠隔振器進行建模，提出一種計算空心圓柱形橡膠隔振器動剛度的方法。

本文測試了用于扭轉減振器橡膠圈的橡膠試片在不同壓縮比下，靜態力與位移關系特性，測試了橡膠試片動態剪切特性與激振頻率的關系。進行了在不同角位移激勵下，橡膠阻尼式扭轉減振器固有頻率的測試。分別用Kelvin-Voigt，Maxwell，分數導數模型描述橡膠材料的黏彈性特性，通過測試得到的橡膠試片的靜動態特性，經計算分析，得到了Kelvin-Voigt，Maxwell，分數導數本構模型中的參數。利用這3種模型計算橡膠扭轉減振器的固有頻率，并和實驗值進行了對比。

1 扭轉減振器固有頻率和橡膠試片靜動態特性實驗

1.1 固有頻率測試實驗

橡膠阻尼式扭轉減振器結構和固有頻率測試實驗分別見圖1和2。

固有頻率測試時，扭轉減振器的慣量環處于自由狀態，輪轂與實驗臺的激振軸緊固聯接。試驗臺激振軸對輪轂產生恒定振幅角位移激勵。在激振軸與慣量環外緣安裝兩個加速度傳感器，測量其切向方向加速度。對慣量環上的加速度與激振軸上的加速度相除，可得到一傳遞函數，其幅值的峰值所對應的頻率為減振器固有頻率。圖3某減振器固有頻率測試曲線，由圖可見，在頻響函數達到峰值時，滯后角為-90°，測試得到的固有頻率為357 Hz。

圖1 橡膠阻尼式扭轉減振器Fig.1 Rubber damped torsional vibration absorber

圖2 橡膠阻尼式扭轉減振器固有頻率測試實驗Fig.2 Measurement of natural frequency of rubber damped torsional vibration absorber

圖3 一扭轉減振器頻率函數測試曲線Fig.3 Measured frequency response function of a torsional vibration abosrber

1.2 橡膠試片剪切試驗

扭轉減振器工作時，聯接輪轂和慣量環的橡膠圈每一小段微元均處于剪切狀態。為了表征其受力情況，得到橡膠阻尼式扭轉減振器中橡膠圈的靜、動態特性參數，本文設計了橡膠試片剪切實驗工裝，進行了與扭轉減振器相同材料的橡膠試片靜動態特性實驗。

橡膠試片剪切實驗工裝見圖4。實驗工裝包括運動塊(1)、鎖緊塊(4，7)、轉接塊(5)、夾板(2，10)、橡膠試片(8)、墊片(9)和緊固螺栓。在MTS831彈性體實驗臺上進行橡膠試片靜態剪切及動態剪切實驗。運動塊(1)與MTS彈性體實驗臺運動端聯接，鎖緊塊(4，7)通過轉接塊(5)與MTS彈性體實驗臺固定端聯接。在運動塊(1)和兩個鎖緊塊(4，7)中間放入兩個橡膠試片。兩個橡膠試片對稱布置。在兩個鎖緊塊外部套上兩個夾板(2，10)。橡膠試片壓縮量與壓縮前的厚度之比定義為壓縮比。在兩個夾板中間可安裝墊片，通過不同厚度墊片的組合可以實現不同的壓縮比。

圖4 橡膠試片剪切實驗工裝三維模型Fig.4 3D model of test fixture for measuring shear performance of a rubber specimen

試驗橡膠試片的材料和尺寸參數見表1。其中A0和d0為橡膠試片壓縮前的面積和厚度。

表1 橡膠試片材料和尺寸參數Tab.1 materials and sizes of rubber specimen

在橡膠試片準靜態試驗時，對橡膠試片的運動塊施加低頻穩態簡諧位移激勵，測試工裝固定端力隨位移激勵的變化關系曲線。為了盡量減少黏性力影響，應采用較低激振頻率和較大激振振幅。本實驗采用激振頻率為0.2 Hz，振幅為5 mm的穩態簡諧位移激勵，分別進行壓縮比從20%到50%的準靜態剪切實驗，記錄橡膠試片的力與位移關系如圖5所示。

圖5 準靜態剪切試驗力與位移關系曲線Fig.5 Measured shear force and displacement of rubber specimen

在橡膠試片動態剪切試驗時，為了盡量減小橡膠塊內部摩擦力的影響，應采用較小激振振幅。本實驗的激振振幅為0.1 mm，激振頻率范圍為5～500 Hz。動態剪切實驗時，記錄橡膠試片剪切動剛度隨激振頻率的變化關系曲線，如圖6所示。橡膠試片靜動態特性的實驗方法和數據處理方法可參考文獻［8］。

圖6 動態剪切試驗動剛度與激振頻率關系曲線Fig.6 Measured dynamic shear stiffness versus excitation frequency of rubber specimen

2 橡膠阻尼式扭轉減振器力學模型

將扭轉減振器的輪轂和慣量環簡化為轉動慣量為Jh和Ji的慣量盤，其扭振角分別為θh和θi，如圖7所示。將橡膠圈簡化為具有剛度和阻尼的黏彈性元件，采用不同方法表征黏性特性時，就得到不同力學模型。

圖7 橡膠阻尼扭轉減振器的等效力學模型Fig.7 Equivalent mechanic models for rubber damped torsional vibration dampers

2.1 Kelvin-Voigt模型

將橡膠圈等效為彈簧和阻尼單元并聯，得到扭轉減振器Kelvin-Voigt力學模型，如圖7(a)所示［3］。圖中，KKV為扭轉剛度，CKV為阻尼。基于Kelvin-Voigt力學模型，扭轉減振器慣量環的振動方程為

2.2 Maxwell模型

扭轉減振器Maxwell力學模型如圖7(b)所示［3-5］，輪轂和慣量環間的相互作用力矩為彈性扭矩、黏彈性扭矩、摩擦扭矩的疊加。其中Ke為彈性扭轉剛度，KMv為黏彈性扭轉剛度，CMv為黏性阻尼，Tf為橡膠摩擦扭矩。摩擦扭矩由橡膠分子之間的摩擦產生，與頻率無關，只與橡膠變形量有關［3］。扭轉減振器Maxwell力學模型可簡化為剛度為KMeq的彈簧、阻尼為 CMeq的彈性、阻尼單元和摩擦單元并聯的等效力學模型，等效扭轉剛度和阻尼分別為：

式中 ω為作用于減振器輪轂處的激振頻率。

基于Maxwell力學模型，減振器慣量環的振動方程為

2.3 分數導數模型

橡膠阻尼式扭轉減振器分數導數力學模型如圖7(c)所示［3］。輪轂和慣量環間的相互作用力矩為彈性扭矩、黏彈性扭矩、摩擦扭矩的疊加。彈性扭矩、摩擦扭矩與Maxwell力學模型一致。黏彈性力矩在的分數導數模型中為［4］

式中 B為黏彈性扭矩系數，α為分數導數階，其取值范圍為0＜α＜1，Dα(·)表示對(·)的α階導數。

基于分數導數力學模型，扭轉減振器慣量環的振動方程為

3 橡膠本構模型參數的確定

在橡膠阻尼式扭轉減振器中，由橡膠圈提供的扭轉剛度和阻尼與橡膠材料的配方、壓縮比和橡膠圈的尺寸有關。為了消除橡膠圈的尺寸對橡膠圈材料動態性能的影響，建立橡膠圈材料的Kelvin-Voigt本構模型、Maxwell本構模型和分數導數本構模型，分別見圖8(a)，(b)，(c)所示。

3.1 Kelvin-Voigt本構模型參數的確定

在橡膠Kelvin-Voigt模型中試片剪切動剛度為

式中 kKV和cKV分別橡膠試片Kelvin-Voigt模型中的剪切剛度和剪切阻尼。

圖8 橡膠本構模型Fig.8 Rubber constitutive model

由動態剪切實驗測得的橡膠試片剪切動剛度隨激振頻率變化的函數關系曲線，根據最小二乘法擬合即可得到kKV和cKV。壓縮比為50%的動剛度測試及擬合曲線如圖9所示。

圖9 Kelvin Voigt模型動剛度測試及擬合曲線Fig.9 Measured dynamic stiffness versus frequency of rubber specimen，and the fitted curve using Kelvin-Voigt model

橡膠材料Kelvin-Voigt模型中，其剪切模量GKV可由橡膠試片的彈性剛度kKV計算得到，黏性系數ηKV可由橡膠試片剪切阻尼cKV計算得到。

圖10為單個橡膠試片剪切狀態示意圖，橡膠彈性剪切模量為

圖10 單個橡膠試片剪切狀態示意圖Fig.10 The shear deformation of a rubber specimen

式中 τKV為橡膠彈性剪切應力，γ為橡膠剪切應變，FKV為試片彈性剪切力，A為試片壓縮后與實驗工裝的接觸面積，x為試片剪切位移，d為試片壓縮后的厚度，分母中的2代表有兩塊試片同時進行試驗。

橡膠試片壓縮比為

假設橡膠試片在壓縮過程中體積保持不變，得

將式(9)和(10)代入式(8)，得

同理，由橡膠試片剪切阻尼cKV計算得到橡膠試片的Kelvin-Voigt本構力學模型的黏性系數ηKV，即

由橡膠試片動態剪切實驗計算得到的橡膠試片Kelvin-Voigt模型參數和本構模型參數如表2所示。對彈性剪切模量GKV和單元黏性系數ηKV與壓縮比的關系進行最小二乘擬合，得到如下關系式：

表2 橡膠試片Kelvin-Voigt模型和本構模型參數Tab.2 Parameters of Kelvin-Voigt model and constitutive model

3.2 Maxwell本構模型參數的確定

當橡膠試片壓縮比已知時，由準靜態實驗，測試得到作用于橡膠試片的剪切力與位移激勵關系曲線，根據參考文獻［3］，可求解得到橡膠試片的彈性剛度ke、最大摩擦力Ffmax和1/2最大摩擦力時的位移x1/2。壓縮比為20%時，準靜態實驗所測試得的一組剪切力與位移激勵函數關系曲線，如圖11所示。

圖11 剪切力與位移關系曲線Fig.11 Shear force versus displacement

由橡膠試片彈性剛度、最大摩擦力、1/2最大摩擦力時的位移，可計算得到Maxwell本構模型中的彈性剪切模量Ge、最大摩擦應力τfmax、1/2最大摩擦應力時的應變γ1/2。

在不同的試片壓縮比下，由橡膠試片準靜態剪切實驗得到的力和位移關系，計算得到的Maxwell橡膠材料本構模型參數如表3所示。對表3的計算結果，應用最小二乘方法，可擬合得到Maxwell本構模型中的參數與橡膠試片壓縮比的關系為：

在橡膠試片Maxwell模型中，當位移激勵幅值較小時，可忽略摩擦力影響，橡膠試片動態剪切剛度隨激振頻率變化的函數關系為［6］

式中 kMv，cMv分別為橡膠試片Maxwell模型中黏彈性單元的剪切剛度和剪切阻尼。利用最小二乘方法，由動態剪切實驗測得的橡膠試片剪切動剛度隨激振頻率變化關系，即可得到kMv和cMv。壓縮比為50%時，橡膠試片剪切動剛度隨激振頻率變化測試曲線與擬合曲線如圖12所示。

橡膠Maxwell本構模型中黏彈性單元的參數，GMv和ηMv，可由橡膠試片黏彈性剪切剛度 KMv和黏彈性剪切模量CMv計算得到：

表3 橡膠試片Maxwell模型摩擦單元參數、Maxwell本構模型中彈性單元與摩擦單元參數Tab.3 Parameters for friction element in Maxwell model of rubber specimen，and parameters for elastic element and friction element in Maxwell constitutive model

圖12 Maxwell模型動剛度測試及擬合曲線Fig.12 Measured dynamic stiffness versus frequency of rubber specimen，and the fitted curve using Maxwell model

在一定的壓縮比下，計算得到的橡膠試片Maxwell模型中黏彈性單元的參數如表4所示。對Maxwell本構模型黏彈性剪切模量GMv與黏性系數ηMv與壓縮比的關系進行最小二乘擬合，得到以下關系式：

表4 橡膠試片Maxwell模型中黏彈性單元參數Tab.4 Parameters for visco-elastic element in Maxwell model

3.3 分數導數本構模型參數的確定

橡膠試片分數導數本構模型中的彈性剪切模量Ge、最大摩擦應力τfmax和1/2最大摩擦應力時的應變γ1/2與Maxwell橡膠試片力學模型中的參數相同。在橡膠試片動態剪切實驗中，當位移激勵幅值較小時，分數導數橡膠試片力學模型可忽略摩擦力影響，橡膠試片動態剪切剛度隨激振頻率變化的函數關系為

式中 b為分數導數橡膠模型中橡膠試片黏彈性系數，α為分數導數橡膠模型中橡膠試片分數導數階。由動態剪切實驗測得的橡膠試片剪切動剛度隨激振頻率變化的函數關系曲線，最小二乘法擬合得到b和α。壓縮比為50%時，橡膠試片剪切動剛度隨激振頻率變化測試曲線與擬合曲線如圖13所示。

圖13 分數導數模型動剛度測試及擬合曲線Fig.13 Measured dynamic stiffness versus frequency of rubber specimen，and the fitted curve using fractional derivative model

橡膠材料本構模型中的單元黏彈性系數m可由橡膠試片中的黏彈性系數b計算得到

計算得到的分數導數模型和分數導數本構模型中的常數見表5。對單元黏彈性系數m和分數導數階α與壓縮比的關系，進行最小二乘擬合，擬合關系式為：

表5 分數導數模型中黏彈性單元模型參數Tab.5 Parameters for visco-elastic element in fractional derivate model

4 橡膠阻尼式扭轉減振器固有頻率計算

4.1 扭轉減振器扭轉剛度與阻尼的計算

通過橡膠試片準靜態剪切實驗和動態剪切實驗得到了3種不同橡膠模型的參數。各種模型中的模型參數與橡膠試片尺寸無關。由3種橡膠本構模型的參數，和扭轉減振器尺寸參數，可計算得到扭轉減振器的剛度和阻尼，進而用于扭轉減振器固有頻率的計算。

任意輪轂形狀的橡膠阻尼式扭轉減振器如圖14所示，扭轉剛度K與橡膠彈性剪切模量G的關系 為［9］

式中 r1(x)為沿橡膠圈軸向變化的輪轂外半徑，r2(x)為沿橡膠圈軸向變化的慣量環內半徑，L為橡膠圈寬度。分別將橡膠試片Kelvin-Voigt模型中彈性剪切模量GKV、Maxwell力學模型彈性剪切模量Ge和黏彈性剪切模量GMv、分數導數力學模型彈性剪切模量Ge帶入式(29)中，即可得到扭轉減振器Kelvin-Voigt模型扭轉剛度KKV、Maxwell模型彈性扭轉剛度Ke和黏彈性扭轉剛度KMv、分數導數模型彈性扭轉剛度Ke。

圖14 任意輪轂形狀橡膠阻尼式扭轉減振器Fig.14 Cross-section of a torsional vibration absorber

扭轉減振器阻尼C與橡膠材料的黏性系數η的關系為

將橡膠試片Kelvin-Voigt本構模型單元黏性系數ηKV、Maxwell本構學模型單元黏性系數ηMv帶入式(30)中，即可得到扭轉減振器Kelvin-Voigt力學模型阻尼CKV、Maxwell力學模型黏性阻尼CMv。

由公式(2)可得橡膠阻尼式扭轉減振器Maxwell力學模型等效扭轉剛度KMeq，由公式(3)可得橡膠阻尼式扭轉減振器Maxwell力學模型等效阻尼CMeq。

任意輪轂形狀的減振器最大摩擦扭矩Tfmax與最大剪切應力τfmax的關系為

式中 r(x)為橡膠圈內外半徑均值。

4.2 基于Kelvin-Voigt模型計算固有頻率

當扭轉減振器中的橡膠圈的動態特性用Kelvin-Voigt模型來表征時，扭轉減振器慣量環的振動方程為公式(1)，其頻率的計算公式為

式中 fKVn為橡膠阻尼式扭轉減振器Kelvin-Voigt模型中無阻尼固有頻率，ζ為阻尼比。由式(32)可得，橡膠阻尼式扭轉減振器Kelvin-Voigt力學模型固有頻率未考慮幅值相關性，即不同角位移激勵計算得到的減振器固有頻率均相等。

4.3基于Maxwell模型計算固有頻率

當扭轉減振器中的橡膠圈的動態特性用Maxwell模型來表征時，扭轉減振器慣量環的振動方程為公式(4)。由于Maxwell模型中的等效剛度KMeq和等效阻尼CMeq均與頻率相關，因此不能利用公式(4)來直接求解固有頻率。可采用在輪轂處加一激勵，當振幅一定時，頻率在一定的范圍內變化，求解扭轉減振器慣量環的振動響應，其峰值即為扭轉減振器的固有頻率。

本文采用中心差分法，求解慣量環的響應。應用有限差分替代位移對時間的導數，采用等時間步長Δt，有：

式中 n為迭代步數。

由式(33)和(34)對慣量環的振動方程(4)在時域進行離散化，并化簡，得到慣量環轉角迭代計算公式：

摩擦扭矩的迭代計算公式計算方法見參考文獻［3］。計算中，慣量環轉角迭代公式和摩擦扭矩迭代公式的初始條件為

作用于扭轉減振器輪轂處位移激勵為單一頻率的簡諧分量，即

式中 Ae為輪轂處角位移激勵幅值，fe為激勵頻率。

綜上所述，基于Maxwell模型，計算扭轉減振器固有頻率步驟為：

(1)當輪轂激勵角振幅為Ae，激勵頻率為fe時，由式(33)～(37)計算得到扭轉減振器的慣量環轉角時域響應；

(2)當慣量環轉角位移的時域響應穩定后，其幅值為該位移激勵下的慣量環穩態響應幅值；

(3)激勵振幅Ae不變，改變激勵頻率fe，可計算不同激勵頻率下的慣量環穩態響應幅值，其中慣量環最大的穩態響應幅值對應的激勵頻率，即為在位移激勵振幅Ae下，橡膠阻尼式扭轉減振器的固有頻率。

(4)改變激勵振幅Ae，重復計算步驟(1)～(3)，即可計算不同位移激勵振幅下橡膠阻尼式扭轉減振器的固有頻率。

4.4 基于分數導數模型計算固有頻率

由于分數導數模型中，表征橡膠黏彈性特性的常數B和α與頻率相關，因此不能直接求解固有頻率。可利用與4.3節類似的方法求解扭轉減振器的固有頻率。

對表征橡膠圈的黏彈性力的公式(5)進行時域離 散 化［2］，得

式中 r0為橡膠圈內外半徑均值，L0為橡膠圈的寬度，h0為橡膠圈內外半徑差。Aj+1為Grunwald系數，定義為

式中 Γ(·)為Gamma函數。

由式(33)，(34)和(38)，對式(6)進行時域離散化，化簡得到慣量環轉角迭代公式為

4.5 固有頻率的計算值與測試值對比

某款橡膠阻尼式扭轉減振器的尺寸和慣量參數如表6所示。由3種不同的模型計算得到的該減振器固有頻率與測試得到的固有頻率對比如圖15所示。

表6 橡膠阻尼式扭轉減振器尺寸和慣量環慣量參數Tab.6 Parameters for a torsional vibration damper

圖15 不同激振角度下的減振器固有頻率計算值與測試值對比Fig.15 Measured and calculated natural frequency of a torsional vibration absorber under different oscillation angles

5 結 論

利用建立的扭轉減振器的模型，和橡膠材料的本構模型，進行了一款橡膠阻尼式扭轉減振器固有頻率計算，并與測試值進行了對比。分析結果表明，Kelvin-Voigt模型計算固有頻率不隨激振振幅的變化而變化，不能很好地反應橡膠的動態特性。Maxwell模型可表征減振器固有頻率隨激振振幅變化而變化，可較為準確地計算小振幅激勵下橡膠阻尼式扭轉減振器固有頻率，但大振幅激勵下固有頻率的計算值與實測值的誤差較大。分數導數力學模型能夠比較準確地計算橡膠阻尼式扭轉減振器在各種激勵振幅下的固有頻率，其相對誤差均在10%以內。

［1］Shangguan Wen-Bin，Pan Xiao-Yong.Multi-mode and rubber-damped torsional vibration absorbers for engine crankshaft systems［J］.International Journal of Vehicle Design，2008，47(1/2/3/4)：176—188.

［2］Tarrago Garcia MJ，Kari L.Frequency and amplitude dependence of the axial and radial stiffness of carbonblack filled rubber bushings［J］.Polymer Testing，2007，26：629—638.

［3］Sjoberg M，Kari L.Non-linear behavior of a rubber isolator system using fractional derivatives［J］.Vehicle System Dynamics，2002，37(3)：217—236.

［4］Pritz T.Analysis of four-parameter fractional derivative model of real solid materials［J］.Journal of Sound and Vibration，1996，195(1)：103—115.

［5］吳杰，上官文斌.采用黏彈性分數導數模型的橡膠隔振器動態特性的建模及應用［J］.工程力學，2008，25(1)：161—166.WU Jie，SHANGGUAN Wenbin.Modeling and applications of dynamic characteristics for rubber isolators using viscolastic fractional derivative model［J］.Engineering Mechanics，2008，25(1)：161—166.

［6］Berg M.Anon-linear rubber spring model for rail vehicle dynamics analysis［J］.Vehicle System Dynamics，1998，30(4)：197—212.

［7］?stberg M，Coja M，Kari L.Dynamic stiffness of hollowed cylindrical rubber vibration isolators-The waveguide solution［J］.International Journal of Solids and Structures，2013，50(10)：1 791—1 811.

［8］Shangguan Wen-Bin，Lu Zhen-Hua.Experimental study and simulation of a hydraulic engine mount with fully coupled fluid structure interaction finite element analysis model［J］.Computers&Structures，2004，82(22)：1 751—1 771.

［9］Tarago G MJ，Kari L，Vinolas J，et al Torsion stiffness of a rubber bushing：a simple engineering design formula including the amplitude，dependence［J］.The Journal of Strain Analysis for Engineering Design，2007，42：13—21.

Astudy on method of calculation and measurement for natural frequency of torsional vibration rubber dampers

SHANGGUAN Wen-bin1，2，WEI Yu-ming1，ZHAO Xu1，RONG Jun-wei2，WANG Ya-jie2
(1.School of Mechanical&Automotive Engineering，South China University of Technology，Guangzhou 510641，China；(2.Ningbo Tuopu Group Co.Ltd.，Ningbo 315800，China)

The natural frequency is one of the most performance parameters of a rubber damped torsional vibration dampers(TVD).The natural frequency is determined by the static and dynamic shear properties of the rubber ring of the TVD，and the moment of inertia of the inertia ring.In this paper a special fixture is designed and used for measuring static and dynamic performance of a rubber shear specimen under different compression ratios.Three constructive models，Kelvin-Voigt model，the Maxwell model and the fractional derivative model，for characterizing the rubber viscous properties are described，and the method for obtaining the model parameters in the three constructive models are developed using the measured dynamic performance of a rubber shear specimen.The natural frequency of a TVD is calculated using the three models to characterize the rubber ring of the TVD，and the calculated frequencies are compared with the measurement.It is demonstrated that the natural frequency of the TVD calculated by the Kelvin-Voigt model is constant when the excitation amplitude is increased，which is inconsistent with the measured data.The calculated natural frequency using the Maxwell and the fractional derivative model to describe the rubber ting of a TVD is the excitation amplitude dependent，which is consist with the trend of measured data.But the relative errors between the calculation using Maxwell model and the measurement under the large amplitude excitation are large，and is not be accepted in engineering application.The natural frequency estimated from the fractional derivative model agrees well with the measured data for different excitation amplitudes，and the relative error between the calculation and measurement is less than 10%.

torsional vibration damper；rubbers；measurement and identification of rubber parameters；natural frequency measurement

U464.133+3

A

1004-4523(2015)04-550-10

10.16385/j.cnki.issn.1004-4523.2015.04.007

上官文斌(1963—)，男，博士，教授，博士生導師。E-mail：sgwb@163.com

2013-08-10；

2015-03-30

國家自然科學基金資助項目(51275175，11472107)