王 雋，楊 俊，張雪冰，周 憶，肖 科

(1.武漢第二船舶設計研究所，湖北 武漢 430064;2.重慶大學 機械傳動國家重點實驗室，重慶 400030)

水潤滑橡膠軸承工作介質為水，在低速重載運行時不易建立潤滑[1]。由于橡膠材料有一定的彈性和非線性特性，在低速、重載時，往往工作在邊界潤滑狀態(tài)，造成軸－軸承之間的運動不連續(xù)，進而出現摩擦振動和噪聲。20世紀70年代開始國外就對水潤滑軸承摩擦導致的振動噪聲問題進行研究，結果表明該振動噪聲特性與軸－軸承組成的摩擦系統(tǒng)密切相關，尤其是摩擦副材料的性能和結構。本文將重點分析橡膠材料特性和橡膠層厚度對水潤滑橡膠軸承振動噪聲特性的影響。

1 水潤滑橡膠軸承振動噪聲機理和數學模型

1.1 水潤滑橡膠軸承振動噪聲機理

水潤滑橡膠軸承的振動噪聲主要是軸和軸承之間的摩擦產生的。目前摩擦噪聲機理主要有:張弛振動機理、摩擦力－速度曲線負斜率機理以及摩擦噪聲的模態(tài)耦合機理[2－3]。當運動界面靜摩擦力大于動摩擦力時，可能導致摩擦粘滑效應，產生振動噪聲;當摩擦力隨運動速度的增加而減小時，摩擦系數為負斜率，可能激發(fā)運動部件產生自激振動;當運動部件的不同固有頻率靠近時，摩擦力能使模態(tài)耦合，激發(fā)強烈振動而產生噪聲。目前被大多數人應用于實際工程問題的機理為模態(tài)耦合機理，即摩擦噪聲是由摩擦耦合誘發(fā)和軸承系統(tǒng)中各部件的模態(tài)參數匹配不當引起的系統(tǒng)不穩(wěn)定現象，從而產生自激振動[4－5]。

水潤滑橡膠軸承系統(tǒng)的穩(wěn)定性不僅取決于結構固有特性，還取決于軸和橡膠之間的耦合關系。本文應用模態(tài)耦合機理分析水潤滑橡膠軸承振動噪聲，并用復模態(tài)分析方法求解。復模態(tài)分析是通過有限元方法求解得到系統(tǒng)的復特征值，根據特征值的實部說明系統(tǒng)是否穩(wěn)定:若實部是正值，相應的虛部被認為是可能發(fā)生的摩擦噪聲的頻率;若模態(tài)具有非正實部，則該階模態(tài)將被認為是穩(wěn)定的。

1.2 水潤滑橡膠軸承振動噪聲數學模型

水潤滑橡膠軸承系統(tǒng)可簡化為如圖1所示的四自由度的軸－軸承模型。

圖1 四自由度軸－軸承系統(tǒng)模型

根據上述模型得到水潤滑橡膠軸承系統(tǒng)的四自由度運動方程(1)。

式(1)由于引入了摩擦力的作用，該系統(tǒng)運動方程的質量矩陣項與剛度矩陣項都為非對稱矩陣，因此該系統(tǒng)的特征值可能為復數，可以通過模態(tài)分析得到該系統(tǒng)的特征值。在進行模態(tài)分析前，需進行非線性靜力分析，以確定系統(tǒng)穩(wěn)態(tài)滑動的狀態(tài)，并得到那些要被添加到質量、剛度矩陣的耦合項。

在四自由度模型的基礎上，對軸承系統(tǒng)做如下簡化:慣性力對軸－軸承接觸面間法向力的影響遠小于彈性力的影響。這樣水潤滑橡膠軸承系統(tǒng)的動力學方程可表示為

其中，[M]、[C]和[K]分別為系統(tǒng)的質量矩陣、阻尼矩陣和剛度矩陣;{x}是位移向量;矩陣K[]f是接觸面間的非對稱接觸摩擦耦合剛度矩陣，它耦合了接觸面之間的法向相對位移和切向摩擦力。

其中，{}Φ為特征向量;λ=a+iω為系統(tǒng)的特征值，其虛部ω反映了振動時的固有頻率，實部a反映了系統(tǒng)運動的穩(wěn)定性。若a為正數，說明該階模態(tài)振幅隨著時間的增加會越來越大，導致系統(tǒng)運動失穩(wěn)，這樣的模態(tài)就被稱為不穩(wěn)定模態(tài)，可能導致摩擦噪聲。

1.3 水潤滑橡膠軸承系統(tǒng)有限元模型

根據水潤滑橡膠軸承結構建立圖2所示的軸承系統(tǒng)有限元模型，系統(tǒng)各部件給定的材料屬性如表1所示。

表1 軸承系統(tǒng)各部件的材料屬性表

有限元模型的邊界條件為:

(1)軸的剛度遠大于軸承剛度，軸設為剛性;

(2)約束軸x、y、z三個方向的平動自由度和x、y方向的轉動自由度;

(3)對軸承其施加x、y、z方向的約束;

(4)載荷以力的形式作用在軸的右端面，方向為y負向;

(5)進行特性條件下的分析時，接觸面之間的摩擦系數取為常數。

具體計算步驟如下[6]:

(1)在軸右端加力，進行非線性靜力計算，從而建立軸－軸承間的接觸狀態(tài);

(2)對軸施加旋轉角速度，進行非線性靜力計算;

(3)使用Lanczos法進行自然模態(tài)提取，以便獲取復模態(tài)計算所必須的投影子空間;

(4)使用子空間投影法進行復特征值提取。

圖2 軸承系統(tǒng)有限元模型及其約束關系

在使用有限元方法求解水潤滑橡膠軸承系統(tǒng)復特征值的過程中，關鍵步驟是將摩擦納入有限元模型中。采用ABAQUS/Standard提供的柔體/柔體之間的面－面接觸形式來對軸和軸承的接觸狀態(tài)進行定義。

2 計算結果和分析

根據上述模型，對不同摩擦系數及比壓時橡膠軸承的不穩(wěn)定模態(tài)進行計算分析，評估上述因素對橡膠軸承振動噪聲的影響。摩擦系數取值為0.02～0.06，軸承比壓取值為0.2 ～0.6 MPa。實際使用情況表明，橡膠軸承在低速下由于潤滑不良極易產生振動和噪聲，因此分析時取軸的轉速35 r/min。

表2 摩擦系數 f=0.02、比壓0.2 MPa時的不穩(wěn)定模態(tài)

表2和圖3為摩擦系數為0.02，軸承比壓為0.2 MPa，軸轉速為35 r/min軸承不穩(wěn)定模態(tài)，及失穩(wěn)傾向性最大時的振型圖。經分析，在3 640～3 750 Hz頻率范圍內有4個不穩(wěn)定模態(tài)，其中當頻率為3 674.6 Hz時失穩(wěn)傾向性最大，可能引起振動噪聲。

圖3 摩擦系數為0.02時的振型圖

2.1 摩擦系數對振動噪聲的影響分析

為分析摩擦系數對振動噪聲的影響，取載荷0.2 MPa，軸轉速為35 r/min，按照1.3節(jié)所述步驟進行計算，計算結果如表3和圖4～圖7所示。

表3 不同摩擦系數軸承失穩(wěn)計算結果

圖4 摩擦系數為0.03時的振型圖

圖5 摩擦系數為0.04時的振型圖

結果表明，摩擦系數越大，失穩(wěn)傾向絕對值增大，產生摩擦噪聲可能性增大。摩擦系數是導致相鄰模態(tài)重合的重要因素，并且摩擦系數的提高會導致更多不穩(wěn)定模態(tài)。

圖6 摩擦系數為0.05時的振型圖

圖7 摩擦系數為0.06時的振型圖

圖8 失穩(wěn)傾向性最大的不穩(wěn)定模態(tài)頻率值隨摩擦系數變化曲線

圖9 失穩(wěn)傾向絕對值與摩擦系數關系曲線

由圖8、圖9可以看出，隨著摩擦系數的增大，最大失穩(wěn)傾向的不穩(wěn)定模態(tài)頻率變化不大;但隨著摩擦系數的增大，失穩(wěn)傾向絕對值增大，產生摩擦噪聲的幾率增加，這與橡膠軸承的實際運行情況是吻合的。

2.2 橡膠彈性模量對振動噪聲的影響

取彈性模量分別為5×106Pa和1×107Pa的橡膠材料進行振動噪聲影響分析。分析時取摩擦系數為 0.02，軸承比壓 0.2 MPa，軸轉速為 35 r/min。分析結果如圖10所示。

圖10 兩種彈性模量橡膠材料的復模態(tài)分析圖

由圖10可以看出，當彈性模量為5×106Pa時，軸承系統(tǒng)在3 000～4 000 Hz時的不穩(wěn)定模態(tài)個數為18個;當彈性模量為1×107Pa，系統(tǒng)不穩(wěn)定模態(tài)個數為22個。兩種彈性模量橡膠材料的失穩(wěn)傾向絕對值相差不大，表明橡膠材料的彈性模量對軸承振動噪聲影響較小。

2.3 橡膠材料密度對振動噪聲的影響

選取密度分別為 1.3×103kg/m3和 1.8×103kg/m3的橡膠材料進行振動噪聲影響分析。分析時條件同2.2，分析結果如圖11所示。

圖11 兩種密度橡膠材料的復模態(tài)分析圖

由圖11可以看出，當橡膠材料密度為1.3×103kg/m3時，軸承系統(tǒng)在3 000～4 000 Hz時的不穩(wěn)定模態(tài)個數為17個;當橡膠材料密度為1.8×103kg/m3時，系統(tǒng)不穩(wěn)定模態(tài)個數為65個。遠遠多于密度較小的橡膠材料;且最大失穩(wěn)傾向的絕對值達到了1.29×10－4，因此密度較小的橡膠材料有利于減小系統(tǒng)不穩(wěn)定模態(tài)的個數。

2.4 橡膠層厚度對振動噪聲的影響

選取橡膠厚度分別為4 mm和8 mm的橡膠材料進行振動噪聲影響分析，分析工況同上，分析結果如圖12所示。

圖12 兩種橡膠層厚度軸承的復模態(tài)分析圖

分析表明，當橡膠層厚度為4 mm時，軸承系統(tǒng)在3 000～4 000 Hz時的不穩(wěn)定模態(tài)個數為15個，對應失穩(wěn)傾向的絕對值很小。當橡膠厚度為8 mm時，系統(tǒng)不穩(wěn)定模態(tài)個數為132個，對應失穩(wěn)傾向的絕對值也很大，在3 886 Hz時，達到1.31 ×10－4，即失穩(wěn)而產生噪聲的可能性很大。

上述仿真計算結果與實際船舶上進行試驗結果基本吻合[7]，船舶試驗結果表明，摩擦振動激起頻率為約3 500 Hz。因此減少橡膠厚度[7－10]有利于減小系統(tǒng)不穩(wěn)定模態(tài)的個數。

3 結論

水潤滑軸承中橡膠材料的摩擦系數、密度，以及結構設計中橡膠層厚度對水潤滑軸承的振動噪聲有較大影響。為減少橡膠軸承出現振動噪聲的幾率，易選擇摩擦系數小、密度小的軸承材料，并設計較薄的橡膠層結構。分析表明，橡膠材料彈性模量變化時軸承的失穩(wěn)傾向絕對值相差不大，因此在分析范圍內橡膠材料的彈性模量對軸承振動噪聲的影響不大。

[1]彭晉民，王家序，余江波，等.水潤滑塑料合金軸承摩擦性能實驗[J].重慶大學學報:自然科學版，2001，24(6):9 －11，35.PENG Jinmin，WANG Jiaxu，YU Jiangbo，YANG Chenyun.Experimental Study on the Frictional Performance of Water Lubricated Plastic Alloy Bearings[J].Journal of Chongqing University:Natural Science Edition，2001，24(6):9 －11，35.

[2]陳光雄，周仲榮，謝友柏.摩擦噪聲研究的現狀和進展[J].摩擦學學報，2000，20(6):478 －481.CHEN Guangxiong，ZHOU Zhongrong，XIE Youbai.Current State and Progress of the Research of Friction Induced Noise[J].Tribology，2000，Vol.(20)6:478 －481.

[3]管迪華，宿新東.制動振動噪聲研究的回顧、發(fā)展與評述[J].工程力學，2004(4):150－155.GUAN Dihua，SU Xindong.An overview on brake vibrations and noise[J].Engineering Mechanics，2004(4):150 － 155.

[4]陳光雄.金屬往復滑動摩擦噪聲的研究[D].成都:西南交通大學，2002.11.Chen Guangxiong.Investigation into Friction － Induced Noise under Metal Reciprocating Sliding Conditions[D].Southwest Jiaotong University，2002.11.

[5]H Ouyang，W Nack，Y Yuan，et al.Numerical analysis of automotive disc brake squeal:a review[J].International Journal of Vehicle Noise and Vibration，2005，1(3/4):207－231.

[6]G DLiles.Analysis of Disc Brake Squeal Using Finite Element Methods[J].1989，SAE Paper 891150.

[7]姚世衛(wèi)，王娟，王雋，等.水潤滑橡膠軸承摩擦行為及試驗研究[J].艦船科學技術，2009，31(12):32 －35.

YAO Shiwei，WANG Juan，WANG Jun，YANG Jun.Vibration and noise generation theory and experiment of water lubricated rubber bearings[J].SHIP SCIENCE AND TECHNOLOGY，2009，31(12):32－35.

[8]楊俊，王雋，周旭輝，等.水潤滑橡膠軸承結構設計[J].艦船科學技術，2011，33(8):103 －107.YangJun，WangJun，Zhouxu － hui，Yao Shi－ wei.Structure research of water－ Lubricated rubber bearings[J].SHIP SCIENCE AND TECHNOLOGY，2011，33(8):103 －107.

[9]王娟，王雋，楊俊.水潤滑橡膠軸承板條設計參數研究[J].艦船科學技術，2012，34(5).

[10]姚世衛(wèi)，楊俊，張雪冰，等.水潤滑橡膠軸承振動噪聲機理分析與試驗研究[J].振動與沖擊，2011，30(2):214－216.YAO Shiwei，YangJun，ZhangXue － bing，WangJuan.Vibration and noise mechanism analysis and tests for water lubrication rubber bearing[J].SHIP SCIENCE AND TECHNOLOGY，2009，1(12):32－35.