金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真

鄒廣平　劉　澤　劉　松　程賀章

哈爾濱工程大學(xué)，哈爾濱，150001

?

金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真

鄒廣平劉澤劉松程賀章

哈爾濱工程大學(xué)，哈爾濱，150001

摘要：提出一種金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真方法，采用正交各向異性阻尼材料模擬金屬橡膠材料。相關(guān)仿真參數(shù)由金屬橡膠減振器正弦掃頻試驗(yàn)確定。對(duì)金屬橡膠減振器進(jìn)行正弦掃頻試驗(yàn)和隨機(jī)振動(dòng)試驗(yàn)，利用有限元仿真軟件ANSYS Workbench對(duì)金屬橡膠減振器隨機(jī)振動(dòng)進(jìn)行有限元仿真模擬，并將仿真結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比分析。研究結(jié)果表明，所提出的有限元仿真方法能夠很好地對(duì)金屬橡膠減振器隨機(jī)振動(dòng)進(jìn)行仿真，且計(jì)算結(jié)果具有較高精度。

關(guān)鍵詞:金屬橡膠減振器；隨機(jī)振動(dòng)；有限元仿真；正交各向異性阻尼材料

0引言

金屬橡膠是一種干摩擦阻尼材料，在受到振動(dòng)載荷時(shí)，金屬絲之間的相對(duì)滑移產(chǎn)生干摩擦力，從而耗散振動(dòng)中的能量，具有承載能力強(qiáng)、耐高低溫、抗老化、剛度可調(diào)節(jié)等優(yōu)點(diǎn)，在工程機(jī)械減振領(lǐng)域中得到了廣泛應(yīng)用[1]。目前，國(guó)內(nèi)對(duì)金屬橡膠材料的研究主要集中在基礎(chǔ)理論、成形工藝、干摩擦減振機(jī)理以及滿足特殊要求而進(jìn)行的專項(xiàng)攻關(guān)[1-5]，對(duì)金屬橡膠減振器動(dòng)力學(xué)特性的研究也以試驗(yàn)和理論研究為主，且大多是在簡(jiǎn)諧激勵(lì)下對(duì)金屬橡膠減振器進(jìn)行動(dòng)態(tài)試驗(yàn)及理論分析[6-13]。金屬橡膠減振器的工作過程中，所受的激勵(lì)大多為脈動(dòng)推力、過載、沖擊等隨機(jī)激勵(lì)，因此，對(duì)隨機(jī)激勵(lì)下金屬橡膠減振器動(dòng)態(tài)減振力學(xué)性能的研究更為重要。目前對(duì)隨機(jī)激勵(lì)下金屬橡膠減振器減振性能的研究相對(duì)較少，且仍然以試驗(yàn)研究為主[14]。對(duì)金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真的研究目前還沒有，且振動(dòng)臺(tái)在實(shí)際應(yīng)用中可能存在試驗(yàn)?zāi)芰Σ蛔?、過試驗(yàn)或欠試驗(yàn)等問題，有限元仿真可以有效解決上述問題，因此，對(duì)金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真的研究顯得極為重要，具有重要工程應(yīng)用意義。

基于上述原因，本文在金屬橡膠減振器正弦掃頻試驗(yàn)基礎(chǔ)上，提出一種金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真方法，采用正交各向異性阻尼材料來模擬金屬橡膠材料，通過金屬橡膠減振器正弦掃頻試驗(yàn)確定金屬橡膠材料的仿真參數(shù)，利用有限元仿真軟件ANSYS Workbench對(duì)金屬橡膠減振器進(jìn)行隨機(jī)振動(dòng)仿真分析，并將仿真結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比分析，驗(yàn)證了仿真方法的可行性。

1金屬橡膠減振器振動(dòng)試驗(yàn)

1.1試驗(yàn)材料和試驗(yàn)設(shè)備

如圖1所示，本文試驗(yàn)所用金屬橡膠減振器結(jié)構(gòu)由底座、支架、緊固螺母、上下端蓋、支撐桿和2個(gè)金屬橡膠減振墊組成。振動(dòng)過程中，通過上下端蓋和支撐桿可以使金屬橡膠減振墊始終處于受壓狀態(tài)，充分發(fā)揮金屬橡膠的減振作用，且可以通過調(diào)節(jié)緊固螺母來控制預(yù)緊量。金屬橡膠減振墊為圓環(huán)形，相對(duì)密度為0.16，外徑為50 mm，內(nèi)徑為12.3 mm，試件高度為20 mm。

圖1　金屬橡膠減振器結(jié)構(gòu)示意圖

試驗(yàn)使用的振動(dòng)試驗(yàn)裝置是江蘇東菱振動(dòng)ES-050-120/LT0202型號(hào)振動(dòng)試驗(yàn)機(jī)，由振動(dòng)控制儀、功率放大器和振動(dòng)發(fā)生機(jī)組成。如圖2所示，試驗(yàn)時(shí)，先在計(jì)算機(jī)上的振動(dòng)測(cè)試軟件中設(shè)置試驗(yàn)的相關(guān)參數(shù)。振動(dòng)控制儀將信號(hào)傳送給功率放大器，信號(hào)經(jīng)過功率放大器放大后傳輸給振動(dòng)發(fā)生機(jī)，使臺(tái)面振動(dòng)并給予減振器穩(wěn)定的振動(dòng)載荷。在減振器和振動(dòng)臺(tái)面上各放一個(gè)壓電式加速度傳感器，來測(cè)量減振前后的加速度響應(yīng)。測(cè)試時(shí)，把傳感器連接到振動(dòng)控制儀上，振動(dòng)控制儀與計(jì)算機(jī)相連，通過計(jì)算機(jī)測(cè)量并儲(chǔ)存試驗(yàn)數(shù)據(jù)。

圖2　振動(dòng)試驗(yàn)原理圖

1.2金屬橡膠減振器正弦掃頻試驗(yàn)

試驗(yàn)掃頻范圍為5～900 Hz，振動(dòng)臺(tái)前期激振控制方式為位移控制，當(dāng)輸入加速度激勵(lì)a達(dá)到設(shè)定值后改為加速度控制。對(duì)金屬橡膠減振器進(jìn)行加速度激勵(lì)為1g～7g的正弦掃頻試驗(yàn)，得到圖3所示的金屬橡膠減振器傳遞率曲線。

圖3　金屬橡膠減振器傳遞率曲線

1.3金屬橡膠減振器隨機(jī)振動(dòng)試驗(yàn)

隨機(jī)振動(dòng)試驗(yàn)采用平均控制方式，試驗(yàn)控制譜如圖4所示，輸入加速度功率譜密度分為3個(gè)階段，上升斜率為3dB/oct，下降斜率為-3dB/oct，隨機(jī)振動(dòng)頻率范圍20～2000 Hz。對(duì)金屬橡膠減振器進(jìn)行振動(dòng)量級(jí)(功率譜密度曲線平譜階段的功率譜密度值)0.01g2/Hz～0.03g2/Hz隨機(jī)振動(dòng)試驗(yàn)，得到不同振動(dòng)量級(jí)下金屬橡膠減振器的均方根加速度，見表1。

圖4　隨機(jī)振動(dòng)試驗(yàn)控制譜

振動(dòng)量級(jí)(g2/Hz)0.0100.0150.0200.0250.030均方根加速度3.22g3.69g4.06g4.38g4.65g

2金屬橡膠材料仿真參數(shù)識(shí)別

采用正交各向異性阻尼材料來模擬金屬橡膠材料，需要輸入的材料參數(shù)包括密度、各方向彈性模量、剪切模量、泊松比和阻尼比。試驗(yàn)中，金屬橡膠減振器為單方向振動(dòng)，為提高計(jì)算精度，除振動(dòng)方向(Z方向)的彈性模量以外，其他方向的彈性模量和所有剪切模量均設(shè)置為10 GPa。在極限載荷范圍內(nèi)，金屬橡膠材料在成形方向受壓變形時(shí)，側(cè)變形為零，所有泊松比均設(shè)置為0，材料的阻尼比和Z方向彈性模量由正弦掃頻試驗(yàn)確定。

隨機(jī)振動(dòng)不同振動(dòng)量級(jí)所對(duì)應(yīng)的輸入均方根加速度、固有頻率和傳遞率峰值如表2所示，其中，固有頻率fn和傳遞率峰值β0由與隨機(jī)振動(dòng)輸入均方根加速度大小相等的加速度激勵(lì)的正弦掃頻試驗(yàn)確定，由于金屬橡膠減振器正弦掃頻試驗(yàn)的激勵(lì)加速度為1g～7g，因此，采用線性插值方法計(jì)算各振動(dòng)量級(jí)所對(duì)應(yīng)的固有頻率fn和傳遞率峰值β0。

表2　不同振動(dòng)量級(jí)下的固有頻率和傳遞率峰值

2.1金屬橡膠材料Z方向的彈性模量

仿真中的金屬橡膠材料既可以承受壓力也可以承受拉力，圖1所示的金屬橡膠減振器，相當(dāng)于2個(gè)金屬橡膠減振墊并聯(lián)。

根據(jù)胡克定律，單個(gè)金屬橡膠減振墊剛度為

k=EzS/h

(1)

式中，k為金屬橡膠減振墊剛度；Ez為Z方向彈性模量；S為金屬橡膠試件承載面積；h為金屬橡膠試件高度。

振動(dòng)系統(tǒng)中，2個(gè)金屬橡膠減振墊并聯(lián)，剛度為原來2倍，根據(jù)振動(dòng)理論，金屬橡膠減振器的固有頻率為

(2)

其中，m為振動(dòng)系統(tǒng)負(fù)載總質(zhì)量，包括重物負(fù)載質(zhì)量(1 kg)、上下端蓋質(zhì)量、緊固螺絲質(zhì)量以及支撐桿質(zhì)量，約為1.2 kg。則Z方向彈性模量為

(3)

2.2金屬橡膠材料阻尼比

將減振器的振動(dòng)過程看成簡(jiǎn)諧激勵(lì)下的單自由度受迫振動(dòng)系統(tǒng)，由振動(dòng)理論可知，在穩(wěn)態(tài)強(qiáng)迫振動(dòng)過程中，振動(dòng)系統(tǒng)的傳遞率為

(4)

式中，λ為頻率比，λ=ω/ωn；ωn為系統(tǒng)圓頻率；ξ為阻尼比。

對(duì)式(4)求極值，令?β/?λ=0，得β0處對(duì)應(yīng)的頻率比：

(5)

將式(5)代入式(4)得到傳遞率β的最大值：

(6)

ξ較小時(shí)，上式可以簡(jiǎn)化為

(7)

則阻尼比

(8)

3金屬橡膠減振器隨機(jī)振動(dòng)仿真

3.1金屬橡膠減振器仿真模型及材料參數(shù)

圖5為金屬橡膠減振器有限元仿真模型的剖面圖，該模型與試驗(yàn)所用減振器完全一致，減振器主體外殼為鋁合金材料，質(zhì)量塊、支撐桿和下端蓋材料為鋼，具體材料參數(shù)如表3所示。

圖5　金屬橡膠減振器仿真模型

材料密度ρ(kg/m3)彈性模量E(GPa)泊松比μ鋼78502100.30鋁合金2770710.33

根據(jù)式(3)、式(8)對(duì)表2中數(shù)據(jù)進(jìn)行處理，得到隨機(jī)振動(dòng)各振動(dòng)量級(jí)下金屬橡膠材料的彈性模量Ez和阻尼比ξ，計(jì)算結(jié)果如表4所示。

表4　金屬橡膠材料仿真參數(shù)

3.2部件接觸定義和約束設(shè)置

金屬橡膠減振器相互連接部件間接觸關(guān)系均為綁定關(guān)系，非連接部件間接觸關(guān)系均為自由關(guān)系。對(duì)減振器底面施加固定約束，由于金屬橡膠減振器沿單軸方向做往復(fù)上下振動(dòng)，為保證模擬的準(zhǔn)確性，對(duì)其他零部件的側(cè)面施加X方向和Y方向的位移約束，使減振器沿Z軸方向振動(dòng)。

3.3模態(tài)分析

采用自由網(wǎng)格劃分方法對(duì)金屬橡膠減振器模型整體進(jìn)行網(wǎng)格劃分，利用ANSYS Workbench的模態(tài)分析模塊對(duì)金屬橡膠減振器模型進(jìn)行模態(tài)分析，各振動(dòng)量級(jí)所對(duì)應(yīng)的金屬橡膠減振器仿真固有頻率(系統(tǒng)一階固有頻率)見表5。

表5　模態(tài)分析仿真結(jié)果

3.4隨機(jī)振動(dòng)仿真及結(jié)果分析

采用與隨機(jī)振動(dòng)試驗(yàn)相同的加速度譜，在減振器底部施加Z向隨機(jī)激勵(lì)，對(duì)金屬橡膠減振器模型進(jìn)行隨機(jī)振動(dòng)分析。表6所示為振動(dòng)量級(jí)0.01g2/Hz～0.03g2/Hz的金屬橡膠減振器Z向1-σ(概率68.269%)隨機(jī)振動(dòng)仿真結(jié)果，其中，減振器底部為輸入均方根加速度仿真結(jié)果，質(zhì)量塊頂部為系統(tǒng)的均方根加速度響應(yīng)仿真結(jié)果。從表6可以看出，輸入均方根加速度的誤差均在5%以內(nèi)，均方根加速度響應(yīng)的誤差均在10%以內(nèi)，仿真誤差較小，由此可見，所提出的仿真方法對(duì)金屬橡膠減振器隨機(jī)振動(dòng)具有很好的仿真效果。

表6　仿真結(jié)果與試驗(yàn)結(jié)果誤差分析

4結(jié)論

(1)提出了一種金屬橡膠減振器隨機(jī)振動(dòng)有限元仿真方法，采用正交各向異性阻尼材料來模擬金屬橡膠，材料的仿真參數(shù)由金屬橡膠減振器的正弦掃頻試驗(yàn)確定，推導(dǎo)了仿真材料Z向彈性模量和阻尼比的計(jì)算公式。

(2)利用ANSYSWorkbench對(duì)金屬橡膠減振器進(jìn)行了隨機(jī)振動(dòng)的仿真，并將仿真結(jié)果與試驗(yàn)結(jié)果進(jìn)行了對(duì)比分析。結(jié)果表明：輸入均方根加速度的誤差不超過5%，均方根加速度響應(yīng)的誤差不超過10%，仿真結(jié)果具有較高的精度，所提出的有限元仿真方法對(duì)金屬橡膠減振器隨機(jī)振動(dòng)具有很好的仿真效果。

參考文獻(xiàn)：

[1]趙程，賀躍進(jìn)，張恒. 金屬橡膠的應(yīng)用研究[J]. 噪聲與振動(dòng)控制，2006(5)：45-47.

ZhaoCheng，HeYuejin，ZhangHen.ApplicationResearchofMeta-rubber[J].NoiseandVibrationControl，2006(5)：45-47.

[2]敖宏瑞，姜洪源，夏宇宏，等. 金屬橡膠干摩擦阻尼系統(tǒng)動(dòng)態(tài)性能分析方法的研究[J]. 中國(guó)機(jī)械工程，2003，14(23)：2053-2056.

AoHongrui，JiangHongyuan，XiaHongyu，etal.ResearchontheDescriptionofDifferentialEquationofMetalRubberIsolationSystem[J].ChinaMechanicalEngineering，2003，14(23)：2053-2056.

[3]王新，朱梓根. 環(huán)形金屬橡膠減振器[J]. 航空動(dòng)力學(xué)報(bào)，1997，12(2)：143-145.

WangXin,ZhuZigen.ARing-likeMetalRubberDamper[J].JournalofAerospacePower，1997，12(2)：143-145.

[4]王平，張國(guó)玉，高玉軍，等. 金屬橡膠減振器在機(jī)載光電吊艙復(fù)合減振系統(tǒng)中的應(yīng)用[J]. 振動(dòng)與沖擊，2014，33(5)：193-198.

WangPing，ZhangGuoyu，GaoYujun，etal.ApplicationofMetal-rubberDampersinVibrationReductionSystemofanAirborneElectro-opticalPod[J].JournalofVibrationandShock，2014，33(5)：193-198.

[5]李宇燕，黃協(xié)清. 承載面積對(duì)金屬橡膠材料本構(gòu)關(guān)系的影響[J]. 振動(dòng)、測(cè)試與診斷，2010，30(5)：544-546.

LiYuyan,HuangXieqing.EffectofLoadedAreaonConstitutiveRelationofMetal-rubber[J].JournalofVibrationMeasurement&Diagnosis，2010，30(5)：544-546.

[6]MaYanhong，LiangZhichao，WangHong，etal.TheoreticalandExperimentalSteady-stateRotorDynamicsofanAdaptiveAirFilmDamperwithMetalRubber[J].JournalofSoundandVibration，2013，332：5710-5726.

[7]ZhangDayi，MaYanhong，HongJie，etal.DynamicMechanicalBehaviorofNickel-basedSuperalloyMetalRubber[J].MaterialsandDesign，2014，56：69-77.

[8]梅小龍，趙俊生，張保成，等. 壓縮長(zhǎng)度對(duì)金屬橡膠減振器剛度的影響[J]. 噪聲與振動(dòng)控制，2013(1)：204-207.

MeiXiaolong，ZhaoJunsheng，ZhangBaocheng，etal.EffectofCompressedLengthontheStiffnessofMetalRubberVibrationAbsorber[J].NoiseandVibrationControl，2013(1)：204-207.

[9]楊春香，周易，張虎. 金屬橡膠動(dòng)態(tài)隔振性能的實(shí)驗(yàn)研究[J]. 航空學(xué)報(bào)，2006，27(3)：536-539.

YangChunxiang，ZhouYi，ZhangHu.ResearchonDynamicPerformanceofMetalRubberDamper[J].ActaAeronauticaetAstronauticaSinica，2006，27(3)：536-539.

[10]ZhangB，LangZQ，BillingsSA，etal.SystemIdentificationMethodsforMetalRubberDevices[J].MechanicalSystemsandSignalProcessing，2013，39(1/2)：207-226.

[11]謝振宇，牟偉興，竇忠才，等. 金屬橡膠環(huán)和磁懸浮阻尼器對(duì)磁軸承轉(zhuǎn)子系統(tǒng)不平衡振動(dòng)的影響[J]. 中國(guó)機(jī)械工程，2010，21(6)：635-638.

XieZhenyu，MouWeixing,DouZhongcai，etal.InfluenceofMetalRubberAnnulusesandMagneticDamperonUnbalanceVibrationofActiveMagneticBearingSystem[J].ChinaMechanicalEngineering，2010，21(6)：635-638.

[12]ZhangRuihua，JiangHongyuan，ZhaoKeding.DampingCharacteristicsAnalysisofSqueezeFilmDamperwithMetalRubber[J].JournalofHarbinInstituteofTechnology，2006，13(2)：146-150.

[13]YanHui，ZhangWenjing，JiangHongyuan，etal.EnergyDissipationofaRing-likeMetalRubberIsolator[J].Chin.Phys.B，2014，23(4)：202-205.

[14]閆輝，姜洪源，劉文劍，等. 金屬橡膠隔振器隨機(jī)振動(dòng)加速度響應(yīng)分析[J]. 物理學(xué)報(bào)，2010，59(6)：4065-4070.

Yan Hui，Jiang Hongyuan，Liu Wenjian，et al. Analysis of Acceleration Response of Metal Rubber Isolator under Random Vibration[J]. Acta Physica Sinica，2010，59(6)：4065-407.

(編輯張洋)

收稿日期：2015-06-17

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目(11372081)

中圖分類號(hào):V252.1

DOI：10.3969/j.issn.1004-132X.2016.14.020

作者簡(jiǎn)介：鄒廣平，男，1963年生。哈爾濱工程大學(xué)航天與建筑工程學(xué)院教授、博士研究生導(dǎo)師。主要研究方向?yàn)閷?shí)驗(yàn)力學(xué)及輕量化結(jié)構(gòu)振動(dòng)與沖擊。發(fā)表論文100余篇。劉澤，男，1989年生。哈爾濱工程大學(xué)航天與建筑工程學(xué)院碩士研究生。劉松，男，1993年生。哈爾濱工程大學(xué)航天與建筑工程學(xué)院碩士研究生。程賀章，男，1988年生。哈爾濱工程大學(xué)航天與建筑工程學(xué)院碩士研究生。

Finite Element Simulation of Metal Rubber Damper Random Vibration

Zou GuangpingLiu ZeLiu SongCheng Hezhang

Harbin Engineering University，Harbin，150001

Abstract:A finite element simulation method of random vibration for metal rubber damper was put forward. The orthotropic anisotropy damping material was used to simulate metal rubber, and the related simulation parameters were determined through the sine sweep experiments of the metal rubber damper. The finite element simulation of metal rubber damper random vibration was done via ANSYS Workbench. Simulation results was compared with the experimental ones. Results show that the proposed finite element simulation method can well simulate metal rubber damper random vibration and the calculation results with high accuracy.

Key words：metal rubber damper; random vibration; finite element simulation; orthotropic anisotropy damping material