干涉配合連接解析分析

姜杰鳳，黃小東

(1.杭州師范大學錢江學院, 浙江 杭州 310036； 2.浙江大學機械工程學系，浙江 杭州 310027)

?

干涉配合連接解析分析

姜杰鳳1，黃小東2

(1.杭州師范大學錢江學院, 浙江 杭州 310036； 2.浙江大學機械工程學系，浙江 杭州 310027)

摘要：干涉配合可有效提高結構疲勞壽命，將軸銷視為彈性體，帶孔板件視為彈塑性體，進行干涉配合連接的理論推導.在彈性范圍內，推導了軸孔接觸應力和最大彈性干涉量公式；在彈塑性范圍內得出塑性區(qū)和干涉量的關系式.同時建立有限元模型，進行數(shù)值模擬.通過實例分析，比較理論推導和有限元結果，相對偏差在可接受范圍，理論公式得到驗證.

關鍵詞：干涉配合；銷插入；彈塑性；干涉量；有限元

過盈配合是一種常用連接技術，在機械領域孔、軸連接中得到了非常廣泛的應用.干涉螺栓連接即螺栓桿和板孔存在過盈，作為一種提高結構疲勞壽命的技術，在飛機裝配中得到應用與發(fā)展[1-4].

國內外學者在干涉配合方面進行了大量的理論推導計算.雅克維茨等[5](前蘇聯(lián))推導了干涉量小于0.4%時板上孔周的應力方程，因材料仍處于彈性范圍內，只適合于低應力水平的結構；在彈塑性范圍內的干涉配合計算相當復雜，劉莉等[6]將干涉配合螺接簡化為平面應力問題進行了解析計算，得出塑性區(qū)尺寸可作為判斷最優(yōu)干涉量的一種度量；吳森等[7]運用平面應變理論進行了應力分析，但是將螺栓作為剛體處理，沒有考慮其彈性變形情況.隨著有限元技術的日漸成熟，Pedersen[8]、Paredes[9]、Zhang[10]、Ozturk[11]等都建立了干涉配合連接的有限元模型，進行數(shù)值模擬計算，并將其和理論推導結果進行比較，結果接近.Ozturk指出對于較為復雜的干涉問題，有限元解比理論解更貼合實際.

之前研究大都將軸銷作為剛體處理，本文視軸銷(螺栓)為彈性變形體，利用厚壁筒理論，對孔板彈性和彈塑性區(qū)別分析；建立干涉銷靜態(tài)安裝有限元模型并進行數(shù)值模擬；通過實例，比較和分析理論解和有限元結果.

1干涉配合彈塑性力學分析

彈塑性力學中厚壁筒理論[12]，當直徑較大的軸銷或螺栓桿(b1)插入孔(a2)中時，存在如圖1干涉配合彈塑性模型.干涉量是干涉配合連接中最重要的參數(shù)，常以相對干涉量的形式被使用，在這里可定義為：

(1)

式中δ為實際干涉量尺寸，b1為軸銷半徑，a2為孔半徑(軸銷中的量用下角標1表示，板孔上的量用下角標2表示).

圖1　干涉配合彈塑性模型Fig. 1　Elastic-plastic model of the interference fit

1.1　彈性變形

當軸銷和孔壁都發(fā)生彈性變形時，如圖1(b)所示，參考文獻[12]中有關彈性理論進行推導.

當軸銷與孔的接觸應力為p時，軸銷的徑向應力、切向應力和位移為：

(2)

(3)

上式為彈性變形的極限干涉量.

1.2　彈塑性變形

當孔板發(fā)生彈塑性變形時，如圖1(c)所示，軸銷與孔壁接觸應力為pp，rp為塑性區(qū)半徑，q為孔板彈性區(qū)和塑性區(qū)交界處的應力(和彈性變形中相同或相似的物理量在塑性量中加’表示).

以上兩q表達式相等，得

(4)

根據(jù)r1=r2=a2=b1假設，則軸銷和孔板的總變形為

(5)

將式(4)代入式(5)，則相對干涉量為

(6)

上式為干涉量和塑性區(qū)半徑之間的解析方程式.

2單層板干涉配合的有限元模型

大型商用有限元軟件ABAQUS，對解決固體力學中的非線性問題具有明顯的優(yōu)勢[13].本文中干涉配合問題涉及到彈塑性和復雜接觸等非線性問題，故借助ABAQUS軟件進行有限元建模和分析.

2.1　幾何模型

為提高計算效率，計算模型采用二維軸對稱結構；為了便于進行有限元計算又不失真，對模型進行必要的簡化.軸銷前端導入部為10°角的倒角，有利于壓入導向并減小壓入力.帶孔板厚度為5mm，外半徑30mm(如銷半徑3mm，外半徑相當于10倍銷半徑)，用較大的板表示無限大板，板件孔入口適當?shù)菇?軸銷(或螺栓桿)和板件初始的裝配模型如圖2所示.選擇增強沙漏控制的CAX4R單元劃分網格.

2.2　接觸對定義

在軸銷干涉配合模擬安裝中存在著明顯的接觸問題，為防止由于材料變形而發(fā)生一種材料嵌入另一種材料中去，須定義主、被動面的接觸對.本例中，軸銷光桿及導入表面作為主動面，孔壁為被動面.其中，接觸算法采用了庫侖摩擦模型，τ=μσ(其中，μ為摩擦系數(shù)).理想條件下，接觸表面摩擦系數(shù)設為0，即忽略摩擦，只考慮導入部對孔壁的擠壓作用.

圖2　有限元模型Fig. 2　Finite element model

圖3　孔開始塑性變形時應力分布Fig. 3　Stress distribution when hole occurs plastic deformation

2.3　邊界條件施加

施加邊界條件情況如圖2所示，在板件底面上，UZ=0；疊層板右端面，UR=0；在軸銷對稱軸上，UR=0，URZ=0；對軸銷上邊施加距離約束以實現(xiàn)軸銷勻速壓下，UZ=-7mm.

2.4　有限元模擬

建立靜態(tài)分析步，使用默認分析步時間1 s.在邊界條件及載荷相同的情況下，通過修改模型中孔徑完成不同干涉量下軸銷靜態(tài)壓入模擬.

3算例與比較

本例選擇Ti-6Al-4V鈦合金螺栓和7050-T7451鋁合金帶孔板(孔徑6 mm)作為研究對象.兩種材料的性能參數(shù)[14]如表1所示.有限元計算中，采用了Mises本構方程.此外，在數(shù)值模擬中，設孔徑不變，通過修改螺栓直徑，來獲得一定大小的干涉量.

表1　兩種材料物理力學性能

3.1　彈性極限干涉量

從式(3)可以看出，最大彈性干涉量和孔徑(a2)沒有關系.將表1中的彈性模量和波松比代入式(3)，得到最大彈性干涉量I=0.33%.

有限元模型中孔徑為6mm (即a2=3mm)時，更改軸銷半徑b1，獲得相應干涉量，進行多次模擬.當干涉量為0.5%(即b1=3.015mm)時，有限元模擬結果如圖3所示，當銷軸進入穩(wěn)定安裝狀態(tài)時，孔壁處應力達到470.9(剛剛超過屈服極限470MPa)，孔開始屈服，所以有限元模擬的最大彈性干涉量為I=0.5%，大于理論推導解.

3.2　塑性區(qū)尺寸

式(6)表明孔發(fā)生彈塑性變形時干涉量和塑性區(qū)半徑的關系，和孔徑(a2)相關.當a2=3mm，b2=30mm時，代入兩種材料的彈性模量、波松比、屈服極限等已知量，公式變?yōu)椋?/p>

(7)

式(7)為非線性方程，借助Matlab軟件進行求解，給定干涉量可獲得相應塑性區(qū)半徑.利用有限元模型在大于0.5%某干涉量下進行模擬計算，獲得應力分布，并測量板件中塑性區(qū)尺寸.例如，干涉量為1.5%時，應力分布如圖4所示，塑性區(qū)半徑為5.82 mm.

圖4　I=2%時應力分布Fig. 4　Stress distribution (I=2%)

圖5　a2=3mm時塑性區(qū)尺寸和干涉量關系Fig. 5　Relation between plastic zone radius andinterference fit size (a2=3mm)

a2=3mm時，將理論計算和有限元結果中的干涉量和塑性區(qū)尺寸關系表示于圖5中.由圖5可見，相同干涉量時，有限元結果中塑性區(qū)尺寸小于理論解.干涉量1%~1.5%時，塑性區(qū)相差較大，干涉量2%~2.5%時，最為接近，達到干涉量3%時，差距又變大.

理論解和有限元結果存在一定差距，對于彈性極限干涉量，相差值為0.17%，對于塑性區(qū)尺寸，最大相對偏差為12%.原因可能為理論推導和有限元模型的簡化條件不完全相同，有限元結果測量位置等.理論解和有限元解的偏差值在可接受范圍內，因此可利用理論公式對彈性極限干涉量、塑性區(qū)尺寸進行估算和預測.

4結論

利用厚壁筒理論，對干涉配合問題進行了理論推導，給出了彈性變形范圍內孔壁接觸應力、彈性極限干涉量，另外得到了彈塑性變形范圍內干涉量和塑性區(qū)半徑的關系方程.

建立了有限元模型并進行模擬，通過實例比較了理論解和有限元結果.理論解和有限元解的偏差在可接受范圍內，故可利用理論推導公式對干涉配合問題中的彈性極限干涉量和塑性區(qū)尺寸進行初步估算.

參考文獻:

[1] 中國航空科學技術研究院.飛機結構抗疲勞-斷裂強化工藝手冊[M].北京：航空工業(yè)出版社，1993.

[2] CHAKHERLOU T N, MIRZAJANZADEH M, ABAZADEH B, et al. An investigation about interference fit effect on improving fatigue life of a holed single plate in joints [J]. European Journal of Mechanics A/Solids,2010,29(4):675-682.

[3] KIM S Y, HENNIGAN D J, KIM D, et al. Fatigue enhancement by interference-fit in a pin-loaded glass fiber-reinforced plastics laminate [J]. Proc IMech E Part C: Journal of Mechanical Engineering Science,2012:226(6):1437-1446.

[4] 姜杰鳳.高鎖螺栓干涉安裝及其對螺接結構力學性能的影響[D].杭州：浙江大學，2014.

[5] 雅柯維茨.飛機長壽命螺栓連接和鉚接技術[M].張國梁，譯.北京：航空工業(yè)出版社，1991.

[6] 劉莉，何慶芝，瞿履和.螺接接頭最優(yōu)干涉問題的研究[J].航空學報，1989，10(12)：640-644.

[7] 吳森.干涉配合緊固件孔的彈塑性工程分析[J].南京航空學院學報，1990，22(4)：17-24.

[8] PEDERSEN P. On Shrink Fit Analysis and Design [J]. Computational Mechanics,2006,37(2):121-130.

[9] PAREDES M, NEFUSSI N, SARTOR M. Study of an interference fit fastener assembly by finite element modeling, analysis and experiment [J]. International Journal of Interactive Design and Manufacture,2012,6(3):171-177.

[10] ZHANG Y, MCCLAIN B, FANG X D. Design of interference fits via finite element method [J]. International Journal of Mechanical Science,2000,42(9):1835-1850.

[11] OZTURK F. Finite-element modelling of two-disc shrink fit assembly and an evaluation of material pairs of discs [J]. Proc IMech E Part C: Journal of Mechanical Engineering Science,2011,225(C2):263-273.

[12] 徐秉業(yè).應用彈塑性力學[M].北京:清華大學出版社,1995.

[13] 趙騰倫.ABAQUS 6.6在機械工程中的應用[M].北京:中國水利水電出版社，2007.

[14] 宋小龍，安繼儒.中外金屬材料手冊[M].北京:化學工業(yè)出版社,2008.

第15卷第1期2016年1月杭州師范大學學報(自然科學版)JournalofHangzhouNormalUniversity(NaturalScienceEdition)Vol.15No.1Jan.2016

Analytical Analysis of Interference Fit Connection

JIANG Jiefeng1, HUANG Xiaodong2

(1.Qianjiang College, Hangzhou Normal University, Hangzhou 310036, China； 2.Department of Mechanical Engineering,

Zhejiang University, Hangzhou 310027, China)

Abstract:Interference fit can effectively improve the fatigue life of aircraft structure. In this paper, the pin is treated as an elastic body, the plate with a hole as the elastic-plastic body, and the interference fit connection is analyzed. In the elastic range, the contact stress between the hole and pin, and the maximum elastic interference value figure are deduced. In the elastic-plastic deformation range, the relation between the plastic zone radius and interference value is obtained. At the same time, the finite element model is established and numerical simulation is conducted. Through an example, theoretical and finite element results are compared. The relative difference is acceptable, and the theoretical figures are verified.

Key words:interference fit; pin insertion; elastic-plastic deformation; interference value; finite element (FE)

文章編號:1674-232X(2016)01-0062-05

中圖分類號:TH131

文獻標志碼:A

doi:10.3969/j.issn.1674-232X.2016.01.012

通信作者：姜杰鳳(1981—)，男，講師，博士，主要從事CAD/CAE/CAM技術.E-mail: jiang19jie81@zju.edu.cn

收稿日期：2015-05-14