肖全山，趙應(yīng)龍，金 著

(1. 海軍工程大學(xué) 振動(dòng)與噪聲研究所，湖北 武漢 430033；2. 船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430033)

0 引 言

橡膠材料具有高彈性，減振效果好，被廣泛應(yīng)用于各類船舶裝置中，起到衰減和吸收裝置振動(dòng)噪聲的作用。橡膠減振器是船舶降噪減振系統(tǒng)中的關(guān)鍵部件，其性能的好壞對(duì)船舶艙室內(nèi)部生活環(huán)境、機(jī)械裝置的使用壽命有重要的影響。靜剛度為評(píng)價(jià)橡膠減振器性能的重要指標(biāo)，其對(duì)減振器傳遞率有著重要的影響，過高的靜剛度會(huì)使得減振器的傳遞率變得很高，而過低的靜剛度又使減振器承載性能太差。目前廣泛以靜剛度為橡膠減振器設(shè)計(jì)的參考指標(biāo)。通常設(shè)計(jì)人員先大體上確定減振器結(jié)構(gòu)后進(jìn)行模具設(shè)計(jì)，制作出減振器的樣品，然后對(duì)樣品的靜剛度進(jìn)行檢測(cè)，如減振器未達(dá)到要求，再對(duì)模具和橡膠膠料硬度進(jìn)行調(diào)節(jié)。通過周而復(fù)始的試驗(yàn)，最終到達(dá)減振器設(shè)計(jì)要求[1]。

運(yùn)用試驗(yàn)反復(fù)驗(yàn)證的方法縱然能實(shí)現(xiàn)減振器的優(yōu)化生產(chǎn)，但此過程需要耗費(fèi)大量的人力物力，成本較高。因此，近年來廣泛應(yīng)用有限元仿真法進(jìn)行減振器的設(shè)計(jì)生產(chǎn)[2-3]，很大程度節(jié)約了成本。

本文選取某型船用減振器為研究對(duì)象，通過理論推導(dǎo)得到靜剛度計(jì)算模型，利用Pro/E三維繪圖軟件建立橡膠減振器幾何模型，再結(jié)合有限元仿真軟件Abaqus對(duì)橡膠減振器的靜剛度計(jì)算結(jié)果進(jìn)行驗(yàn)證，并利用靜剛度計(jì)算模型進(jìn)行結(jié)構(gòu)尺寸的設(shè)計(jì)，很大程度上縮短了減振器設(shè)計(jì)周期。

1 靜剛度模型

橡膠減振器在各類船舶裝置上使用廣泛，具有性能穩(wěn)定、減振效果優(yōu)良等特點(diǎn)，此外，橡膠減振器的尺寸范圍大，在狹窄的艙內(nèi)空間內(nèi)使用仍然不受限制[4]。靜剛度作為確保減振系統(tǒng)具有良好力學(xué)性能的關(guān)鍵參數(shù)，應(yīng)保持在一個(gè)相對(duì)恒定的水平。由于金屬和橡膠的粘結(jié)強(qiáng)度高于橡膠強(qiáng)度本身，因此，橡膠減振器的靜剛度主要由橡膠元件的性能決定。

圖1為某型橡膠減振器的結(jié)構(gòu)簡圖，虛線填充部分表示橡膠材料，承擔(dān)著整個(gè)減振器的減振功能，無填充部分為金屬件，起著限位和支撐作用。

如圖1所示，橡膠受到切向應(yīng)力和正向應(yīng)力，產(chǎn)生正應(yīng)變和切應(yīng)變，在正應(yīng)力下，橡膠材料剛度為[5-6]：

式中：E為橡膠材料的彈性模量；mz為正向形狀系數(shù)；A為受壓面積；H為橡膠材料原始厚度。在切應(yīng)力下，橡膠材料的剪切剛度為：

式中：G為剪切模量；mx為切向形狀系數(shù)；S為受剪面積；L為原始高度。

對(duì)橡膠元件進(jìn)行分析，得到垂直于橡膠—金屬接觸面的正向剛度為

表1 橡膠減振器初始尺寸Tab. 1 The initial size of rubber absorber

2 靜剛度仿真

2.1 材料本構(gòu)模型及參數(shù)

橡膠減振器的彈性體為丁腈橡膠，其具有耐油性好、阻尼大、耐熱性高等優(yōu)點(diǎn)，廣泛作為一般動(dòng)力裝置的減振材料。丁腈橡膠的應(yīng)力應(yīng)變曲線程典型的非線性，其本構(gòu)模型較為復(fù)雜，不能簡單地像描述線性材料那樣描述橡膠材料。關(guān)于橡膠材料本構(gòu)模型橡膠已有大量研究[7-8]，在小變形情況下，可采用多項(xiàng)式本構(gòu)模型來模擬超彈屬性，其應(yīng)變能可表示為：

式中：W（I1，I2）為應(yīng)變能函數(shù)；Cij為Rivlin系數(shù)；I1和I2分別為第一和第二Green應(yīng)變不變量；J為橡膠變形前后體積比；Di決定橡膠材料是否可壓縮。如果采用兩參數(shù)的Mooney-Rivlin本構(gòu)模型，則式（6）變?yōu)?/p>

彈性常數(shù)為：

式中：G為橡膠剪切模量；E為橡膠彈性模量；ν為泊松比，當(dāng)丁腈橡膠看作不可壓縮材料時(shí)，J=1；，此時(shí)，E=3G。當(dāng)橡膠材料參數(shù)C10和C01已知時(shí)，可通過式（8）得到橡膠材料的彈性模量E和剪切模量G，并作為靜剛度理論計(jì)算的已知量輸入。

本文采用Mooney-Rivlin本構(gòu)模型來描述橡膠材料，采用線性本構(gòu)來描述金屬件，材料參數(shù)如表2所示。

表2 材料本構(gòu)模型及參數(shù)Tab. 2 The parameter of material constitutive model

2.2 不同載荷下位移仿真計(jì)算

由于橡膠減振器在生產(chǎn)過程中，橡膠材料與金屬件表面經(jīng)過特殊的硫化處理，其接觸面粘接強(qiáng)度甚至比橡膠材料本身強(qiáng)度更高，因此，橡膠與金屬件接觸面可定義為三向固定的邊界條件，得到如圖2所示的有限元分析模型，橡膠部分單元用雜交型單元C3D8H模擬，金屬部分用減縮型C3D8R單元模擬。

在橡膠減振器上表面建立耦合的約束，使整個(gè)上表面受到合力與耦合點(diǎn)相同，在耦合約束點(diǎn)施加向下的載荷，對(duì)橡膠減振器進(jìn)行靜力分析，得到不同載荷下橡膠減振器的載荷—位移關(guān)系（見圖3），從圖中可以看出，位移和載荷約成線性關(guān)系，說明橡膠在小變形下可以看成線彈性材料。圖4表示橡膠元件在11 kN下的等效應(yīng)力，從圖中可以看出，橡膠元件在與金屬接觸面上受到的應(yīng)力最大，因此橡膠與金屬面的粘結(jié)力決定著整個(gè)橡膠減振器的剛度上限。

2.3 靜剛度仿真計(jì)算

在橡膠減振器上表面耦合約束點(diǎn)施加2 kN和10 kN的載荷，計(jì)算得到圖5和圖6橡膠減振器的位移云圖，從圖中可以看出，橡膠減振器的最大變形發(fā)生在橡膠件下端面，上表面垂向位移為0.56 mm和2.556 mm，可計(jì)算得到靜剛度為3 636.4 N/mm和3 912.4 N/mm。理論計(jì)算結(jié)果為3 726.8 N/mm。靜剛度理論結(jié)果和仿真結(jié)果誤差范圍均在5%以內(nèi)，橡膠減振器的仿真結(jié)果與理論結(jié)果吻合較好，可以判斷所建立的靜剛度理論計(jì)算模型基本無誤。

3 結(jié)構(gòu)尺寸設(shè)計(jì)

當(dāng)橡膠減振器幾何尺寸H，B，L恒定時(shí)，調(diào)整橡膠元件與水平基座的角度Φ，可得到如圖減振器的靜剛度變化曲線，從圖7可以看出，當(dāng)角度Φ增大時(shí)，減振器的垂向靜剛度逐漸減小。當(dāng)橡膠減振器的幾何尺寸B，L，Φ恒定時(shí)，調(diào)整橡膠元件的厚度H，得到靜剛度與橡膠元件厚度H的關(guān)系（見圖8）。

當(dāng)橡膠元件厚度增大時(shí)，橡膠減振器剛度急劇下降，繼續(xù)增加橡膠元件厚度，減振器剛度下降趨勢(shì)漸緩，靜剛度約恒定在2 000 N/mm。此外，橡膠元件為減振器減振部件，增加橡膠元件厚度能顯著提高減振效率，但由于船艙內(nèi)空間狹窄，橡膠減振器高度受限，本文在設(shè)計(jì)橡膠減振器外形結(jié)構(gòu)尺寸時(shí)，需要滿足3點(diǎn)：1）減振器靜剛度必須達(dá)到使用要求；2）限位金屬基座高度不能太低，否則失去限位功能；3）橡膠垂向高度不能超過金屬限位基座高度的1/2，否則當(dāng)減振器失效破壞時(shí)，垂向擺動(dòng)過大，使機(jī)器損壞，但也不能低于限位器高度的1/4，否則減振性能將大大降低。根據(jù)上述3點(diǎn)要求，建立橡膠減振器設(shè)計(jì)約束條件，對(duì)橡膠減振器進(jìn)行設(shè)計(jì)。

圖9為橡膠減振器橫截面輪廓，根據(jù)橡膠減振器設(shè)計(jì)要求及幾何尺寸關(guān)系，建立橡膠如下減振器的約束方程

式中：s為限位器高度，且。式（9）為橡膠垂向高度與金屬限位器高度關(guān)系，式（10）假定靜剛度值為4 000 N/mm。由于船艙內(nèi)實(shí)際空間狹小，機(jī)器設(shè)備重量較大，以初始尺寸為設(shè)計(jì)基礎(chǔ)，假定橡膠減振器的靜剛度和限位器高度，利用約束方程求解在規(guī)定條件下的減振器尺寸最優(yōu)值。

如圖10所示，圖中三維網(wǎng)格曲面為靜剛度變化趨勢(shì)曲面，兩深色曲線為約束條件，在兩曲線之間與三維曲面相交部分所形成的有限長度曲線為滿足假定剛度和限位器高度的橡膠減振器尺寸，從圖10中可以看出，可將將此曲線簡化為線性線段，得到滿足約束條件下的橡膠元件厚度H與Φ的關(guān)系：

通過計(jì)算式（11），可得到在滿足假定靜剛度和限位器高度條件下，不同Φ對(duì)應(yīng)的橡膠元件厚度H。從而為橡膠減振器設(shè)計(jì)提供參考。

4 結(jié) 語

本文以某船用橡膠減振器為研究對(duì)象，通過理論分析得到減振器的靜剛度計(jì)算公式，利用有限元軟件建立橡膠減振器仿真分析模型，計(jì)算初始結(jié)構(gòu)尺寸下橡膠減振器的靜剛度，得到理論計(jì)算結(jié)果與仿真計(jì)算結(jié)果誤差范圍在5%以內(nèi)，證明建立的橡膠減振器靜剛度理論模型滿足要求。根據(jù)理論模型，對(duì)橡膠減振器進(jìn)行結(jié)構(gòu)尺寸設(shè)計(jì)，得到在額定靜剛度和限位器高度下，橡膠減振器中橡膠元件厚度和與水平角度關(guān)系，可為此類型橡膠減振器的設(shè)計(jì)生產(chǎn)提供參考。

參考文獻(xiàn)：

[1]第海龍. 某型重卡橡膠減振器優(yōu)化設(shè)計(jì)[D]. 長沙: 湖南大學(xué),2011.DI Hai-long. The optimization of rubber spring of heavy-duty truck[D]. Changsha: Hunan University, 2011.

[2]IAVORNIC C, PRAISACH Z I, VASILE O, et al. Study of stress and deformation in elastomeric isolation systems using the Finite Element Method[C]//International Conference on Siqnal Processing, Computational Geometry & Artificial Vision, 2011:239–244.

[3]WENBIN S, ZHENGHUA L. Finite element analysis of static elastic characterization of the rubber isolators in automotive dynamic system[J]. Sae Transaction, 2003(112): 185–193.

[4]李中郢, 盧正人. 金屬橡膠減振器組合剛度特性研究[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2005(10): 23–25.LI Zhong-ying, LU Zheng-ren. Research on combined stiffness characteristic of metal rubber damper[J]. Journal of Harbin Institute of Technology, 2005(10): 23–25.

[5]嚴(yán)濟(jì)寬. 機(jī)械振動(dòng)隔離技術(shù)[M]. 上海: 上海科學(xué)技術(shù)文獻(xiàn)出版社, 1985.

[6]程昌鈞, 朱媛媛. 彈性力學(xué)[M]. 上海: 上海大學(xué)出版社, 2005.

[7]孫延奎. 橡膠減振器的設(shè)計(jì)與優(yōu)化[D]. 青島: 青島科技大學(xué),2016.

[8]余天超, 劉夫云, 孫永厚. 基于ABAQUS的橡膠減振器參數(shù)化優(yōu)化設(shè)計(jì)[J]. 計(jì)算機(jī)仿真, 2013, 30(10): 119–202.YU Tian-chao, LIU Fu-yun, SUN Yong-hou. Rubber shock absorber based on ABAQUS parametric optimization design [J].Computer Simulation , 2013, 30(10): 119–202.