基于接觸面特征的螺栓聯(lián)接剛度研究

章伊華，龐奎，林丹益，麥云飛，孫曉萌，楊國玉

(1.寧波大紅鷹學(xué)院 汽車CAE應(yīng)用技術(shù)研究所, 浙江 寧波 315175; 2.力帆實(shí)業(yè)(集團(tuán))股份有限公司 摩托車研究院，重慶 400700;3.上海理工大學(xué) 機(jī)械工程學(xué)院，上海 200093; 4.EASi Engineering, Detroit 48201, MI, US)

基于接觸面特征的螺栓聯(lián)接剛度研究

章伊華1，龐奎2，林丹益1，麥云飛3，孫曉萌3，楊國玉4

(1.寧波大紅鷹學(xué)院 汽車CAE應(yīng)用技術(shù)研究所, 浙江 寧波 315175; 2.力帆實(shí)業(yè)(集團(tuán))股份有限公司 摩托車研究院，重慶 400700;3.上海理工大學(xué) 機(jī)械工程學(xué)院，上海 200093; 4.EASi Engineering, Detroit 48201, MI, US)

螺栓聯(lián)接由于其可靠和有效性在工業(yè)界得到廣泛應(yīng)用，螺栓聯(lián)接設(shè)計(jì)中的聯(lián)接剛度設(shè)計(jì)在螺栓聯(lián)接結(jié)構(gòu)中具有十分重要的作用。提出一種最佳預(yù)緊力下螺栓聯(lián)接結(jié)構(gòu)的有限元分析方法，在綜合考慮螺栓螺母與被聯(lián)接件及被聯(lián)接件之間的受力、接觸、摩擦等非線性因素的基礎(chǔ)上，采用不同結(jié)構(gòu)的系列參數(shù)模型運(yùn)算結(jié)果進(jìn)行聯(lián)接剛度分析，研究了螺栓預(yù)緊力、被聯(lián)接板厚度、螺栓公稱直徑、被聯(lián)接板外徑及螺母系數(shù)等因素對(duì)聯(lián)接剛度的影響。通過接觸面壓應(yīng)力分布實(shí)驗(yàn)驗(yàn)證和連續(xù)性假設(shè)，對(duì)聯(lián)接結(jié)構(gòu)剛度函數(shù)式進(jìn)行無量綱化推導(dǎo)，實(shí)現(xiàn)了系列參數(shù)模型無量綱化剛度數(shù)據(jù)的有效擬合，在有限元分析結(jié)果基礎(chǔ)上獲得了計(jì)算螺栓聯(lián)接結(jié)構(gòu)聯(lián)接剛度的計(jì)算公式。通過與有限元分析結(jié)果及文獻(xiàn)[4-6]結(jié)果比較，驗(yàn)證了該剛度計(jì)算公式的有效性和可靠性。

固體力學(xué)；參數(shù)模型; 接觸面積; 接觸面應(yīng)力分布; 聯(lián)接剛度; 無量綱化

0 引言

螺栓聯(lián)接由于使用方便和聯(lián)接可靠，是當(dāng)今工業(yè)中重要的和最常見的一種機(jī)械聯(lián)接方式。螺栓結(jié)構(gòu)通過施加預(yù)緊力使多個(gè)接觸面之間產(chǎn)生緊密結(jié)合的壓緊應(yīng)力，這種由壓緊應(yīng)力產(chǎn)生的聯(lián)接剛度會(huì)直接影響結(jié)構(gòu)的力學(xué)屬性和行為，這是螺栓聯(lián)接結(jié)構(gòu)區(qū)別于其它結(jié)構(gòu)的重要特點(diǎn)。如何準(zhǔn)確計(jì)算螺栓聯(lián)接結(jié)構(gòu)剛度，如螺栓螺母與被聯(lián)接件聯(lián)接剛度及被聯(lián)接件之間的聯(lián)接剛度，對(duì)于預(yù)測在工作載荷作用下機(jī)械整體結(jié)構(gòu)的可靠性、穩(wěn)定性及振動(dòng)響應(yīng)特性等都具有重要意義，這也是國內(nèi)外學(xué)者研究的重點(diǎn)。Shigley等[1]把壓應(yīng)力在被聯(lián)接件中的分布區(qū)域等效為圓錐體，認(rèn)為在垂直于軸截面方向上均勻分布，如圖1所示，由此推導(dǎo)出關(guān)于錐角的聯(lián)接剛度計(jì)算公式為

(1)

式中：θ為壓應(yīng)力錐角，假設(shè)錐角為30°或45°；E為被聯(lián)接件的彈性模量；β為螺母系數(shù)；L為被聯(lián)接件聯(lián)接長度；d為螺栓公稱直徑；α為孔隙系數(shù)。該方法沒有考慮被聯(lián)接件的寬度、相互間接觸、摩擦等因素對(duì)聯(lián)接剛度的影響，且假設(shè)應(yīng)力在垂直于軸截面上均勻分布與實(shí)際不符。

圖1 螺栓聯(lián)接壓力錐角模型Fig.1 Bolt joint pressure angle model

Wileman等[2]及Musto等[3]通過理論方法與有限元分析相結(jié)合，用擬合有限元結(jié)果的方式推導(dǎo)出被聯(lián)接剛度的計(jì)算公式：

km=EdAexp (B(d/L)),

(2)

式中：A、B常量值取決于材料屬性，對(duì)于鋼材(E=206.8 GPa，泊松比μ=0.291)，常量A=0.787 15，B=0.628 73. (2)式忽略了螺母大小對(duì)剛度的影響，且當(dāng)d/L比值較大時(shí)這個(gè)曲線擬合將失效，一般適用范0.1

(3)

式中：da,gr為接觸面積直徑。

同時(shí)還推導(dǎo)出θ錐角的計(jì)算公式，da,gr和tanθ可以表示為

(4)

da,gr=βd+Ltanθ.

(5)

本文中接觸面直徑用dc表達(dá)，但文獻(xiàn)[6]中仍沿用了壓應(yīng)力值在接觸面直徑dc范圍內(nèi)均勻分布假設(shè)，因此仍產(chǎn)生較大偏差。

Sethuraman等[7]在有限元模型中把螺栓頭或螺母支承面簡化為剛性結(jié)合面和柔性結(jié)合面，并分別施加均勻軸向位移和壓力，再通過修正因子的方式模擬實(shí)際螺栓聯(lián)接結(jié)合狀況，沒有正確反映螺栓聯(lián)接結(jié)構(gòu)的真實(shí)狀態(tài)。要想獲得準(zhǔn)確的螺栓聯(lián)接剛度，首先要明確接觸面上的壓應(yīng)力分布規(guī)律和力學(xué)屬性。Marshall等[8]運(yùn)用現(xiàn)代聚焦超聲波測量技術(shù)通過對(duì)螺栓壓緊接觸面的超聲掃描，得到了比較清晰的壓應(yīng)力分布，章伊華等[9]用有限元方法通過多元混合數(shù)值仿真模型，進(jìn)一步論證了不同螺栓預(yù)緊力作用下接觸面壓應(yīng)力分布規(guī)律及對(duì)接觸面聯(lián)接結(jié)構(gòu)的影響，經(jīng)與實(shí)驗(yàn)對(duì)比驗(yàn)證獲得良好的一致性[8-9]。

本文基于螺栓聯(lián)接結(jié)構(gòu)最佳預(yù)緊力設(shè)計(jì)[10]的三維數(shù)值仿真模型，遵從螺栓預(yù)緊力作用下接觸面壓應(yīng)力分布規(guī)律的實(shí)際特征[9]，通過選用不同參數(shù)的變化計(jì)算出系列模型下的聯(lián)接結(jié)構(gòu)剛度，并以無量綱化方法給出不同系列變化的有限元分析(FEA)參數(shù)模型剛度值曲線，推導(dǎo)出其擬合曲線公式，通過對(duì)比文獻(xiàn)[4-6]的剛度曲線驗(yàn)證了模型和擬合曲線的有效性。

1 有限元模型

考察螺栓聯(lián)接結(jié)構(gòu)的接觸受力和聯(lián)接剛度關(guān)系，建立螺栓聯(lián)接結(jié)構(gòu)預(yù)緊力設(shè)計(jì)的有限元模型，如圖2所示。在結(jié)合面用接觸單元組成接觸對(duì)，并用罰函數(shù)法[11]求解其接觸關(guān)系，再以彈簧單元和其他有限單元建立混合單元模型。螺栓、螺母材料參數(shù)分別?。簭椥阅Ａ縀=2.07×105MPa，泊松比μ=0.3，屈服極限σs=924 MPa，強(qiáng)度極限σb=1 130 MPa. 被聯(lián)接件材料參數(shù)?。簭椥阅Ａ縀=2.07×105MPa，泊松比μ=0.3，屈服極限σs=355 MPa，強(qiáng)度極限σb=600 MPa. 對(duì)每組模型先后施加10 kN、20 kN預(yù)緊力，在板與板以及螺栓螺母與板之間的各個(gè)接觸面間建立摩擦接觸。考慮螺栓聯(lián)接結(jié)構(gòu)的材料屬性和表面粗糙度取摩擦系數(shù)為0.3.

根據(jù)需要分別確定不同螺栓公稱直徑d、螺母系數(shù)β、孔隙系數(shù)α、被聯(lián)接件厚度L和被連接件處部直徑da等的系列參數(shù)值，并計(jì)算和輸出結(jié)果。

2 大寬度聯(lián)接的剛度計(jì)算

2.1 參數(shù)的無量綱化

如圖1所示，按傳統(tǒng)計(jì)算方法應(yīng)力錐角完全展開所對(duì)應(yīng)的被聯(lián)接件寬度稱為大寬度被聯(lián)接件，即da>dc，以往傳統(tǒng)文獻(xiàn)一般取dc=βd+Ltanθ.

可以看到用傳統(tǒng)方法如文獻(xiàn)[4-6]計(jì)算螺栓被聯(lián)接剛度的公式不僅繁復(fù)，且產(chǎn)生誤差也很大，而在工程應(yīng)用中對(duì)螺栓聯(lián)接結(jié)構(gòu)各個(gè)不同參數(shù)取值的每種組合模型進(jìn)行有限元建模技術(shù)要求高、難度大。本文在深入研究螺栓聯(lián)接預(yù)緊力優(yōu)化設(shè)計(jì)[10]和螺栓聯(lián)接結(jié)構(gòu)接觸壓應(yīng)力分布規(guī)律[9]基礎(chǔ)上，推導(dǎo)出符合實(shí)際受力環(huán)境，更為簡捷，且能準(zhǔn)確擬合有限元結(jié)果的聯(lián)接剛度求解公式，實(shí)現(xiàn)對(duì)螺栓結(jié)構(gòu)聯(lián)接剛度的精確設(shè)計(jì)。

為此對(duì)系統(tǒng)參數(shù)進(jìn)行無量綱化處理，由于本文引用的系列模型已經(jīng)確定被壓接觸面上的壓應(yīng)力分布規(guī)律[9]，且確定所有螺栓聯(lián)接結(jié)構(gòu)遵循最佳預(yù)緊力設(shè)計(jì)原則[10]，在此前提下螺牙結(jié)構(gòu)的最大應(yīng)力為屈服前彈性狀態(tài)，這也是所有螺栓結(jié)構(gòu)的最佳工作狀態(tài)，即最大彈性力拉緊接觸狀態(tài)。在此狀態(tài)下結(jié)構(gòu)的接觸面大小和應(yīng)力強(qiáng)度分布決定了螺栓聯(lián)接剛度。因此被聯(lián)接結(jié)構(gòu)剛度只取決于本身結(jié)構(gòu)的固有特性和接觸面積區(qū)域內(nèi)應(yīng)力強(qiáng)度，外載荷不再需要進(jìn)行無量綱化處理。

設(shè)螺栓受拉壓載荷時(shí)的剛度為

(6)

(7)

由此導(dǎo)出被聯(lián)接剛度km也可以表示為比值d/L的函數(shù)形式：

(8)

式中：Em為被聯(lián)接件材料的彈性模量。

考慮到螺母系數(shù)和間隙系數(shù)的影響，并把被聯(lián)接板看作是一個(gè)外徑為βd、內(nèi)徑為αd的圓環(huán)管，先假設(shè)圓環(huán)管的應(yīng)力是均勻分布的，則其剛度可以表示為

(9)

變形為

(10)

當(dāng)然被聯(lián)接件不可能完全是一個(gè)圓環(huán)管，但是其剛度形式必然與(10)式有相似之處。假設(shè)將(β-α)d/L作為一個(gè)自變量，則聯(lián)接剛度可以表示為

(11)

圖3 (3)式和FEA結(jié)果無量綱化Fig.3 Dimensionless of Eq.(3) and FEA results

2.2 大寬度聯(lián)接剛度計(jì)算公式

應(yīng)以控制徑流系數(shù)為重點(diǎn)，強(qiáng)化地表徑流控制和人為生產(chǎn)建設(shè)活動(dòng)土石方綜合利用；強(qiáng)化雨水控制和集蓄利用，提高雨洪利用程度；適當(dāng)拓展河道空間，建設(shè)濱河（湖）綠帶；強(qiáng)化城區(qū)河湖水質(zhì)改善；加強(qiáng)小流域內(nèi)排洪水系的連通與疏浚，降低城市防洪壓力。

根據(jù)上述無量綱化方法進(jìn)行不同螺栓、不同螺母系數(shù)和孔隙系數(shù)、不同聯(lián)接厚度等多變量參數(shù)模型的有限元建模計(jì)算和數(shù)據(jù)采集，圖4為系列計(jì)算模型位移數(shù)據(jù)云圖之一，通過對(duì)螺母與被聯(lián)接件接觸面進(jìn)行位移數(shù)據(jù)采集，并對(duì)這些數(shù)據(jù)進(jìn)行積分和均勻化處理，獲得被聯(lián)接件位移剛度。

圖4 計(jì)算模型位移數(shù)據(jù)云圖Fig.4 Nephogram of displacement data of calculation model

采用2次函數(shù)式進(jìn)行曲線擬合，得到大寬度聯(lián)接剛度計(jì)算公式為

(12)

用此剛度計(jì)算公式(12)式得到的擬合曲線如圖5所示。

圖5 剛度擬合經(jīng)驗(yàn)公式曲線圖Fig.5 Fitted curve of stiffness empirical formula

需要指出的是上述剛度計(jì)算公式(12)式計(jì)算簡易方便，擯棄以往用固定錐角θ代替應(yīng)力分布關(guān)系的假設(shè)，避免由此引起的誤差。由于曲線擬合是基于FEA模型的綜合結(jié)果，可以說該公式是在da>dc時(shí)全面反映螺栓聯(lián)接結(jié)構(gòu)接觸、摩擦、接觸面分離等實(shí)際情況的。da>dc也是螺栓聯(lián)接結(jié)構(gòu)實(shí)際應(yīng)用的最普遍情況，當(dāng)聯(lián)接板的寬度直徑大于dc時(shí)，聯(lián)接剛度不再隨寬度增加而改變，此時(shí)螺栓結(jié)構(gòu)的聯(lián)接剛度趨于漸近值。

3 有限寬度聯(lián)接的剛度計(jì)算

從FEA模型曲線被聯(lián)接件寬度對(duì)聯(lián)接剛度的影響(見圖6)中可以發(fā)現(xiàn)，被聯(lián)接件寬度從da=βd到da>dc，剛度是隨da的增大而增大，最終趨近于漸近值。當(dāng)被聯(lián)接件寬度在有限區(qū)間即βd≤da≤dc時(shí)，根據(jù)螺栓聯(lián)接接觸面壓應(yīng)力分布的連續(xù)性原理[9]，本文推導(dǎo)出能夠連續(xù)擬合的剛度計(jì)算公式，并實(shí)現(xiàn)與da>dc時(shí)剛度計(jì)算公式(12)式的連續(xù)對(duì)接，這與實(shí)際螺栓聯(lián)接結(jié)構(gòu)和FEA結(jié)果的連續(xù)性一致。

圖6 被聯(lián)接件寬度對(duì)剛度的影響Fig.6 Influence of plate diameter on stiffness

當(dāng)da≤βd時(shí)螺栓聯(lián)接剛度為

(13)

式中：t是應(yīng)力分布修正系數(shù)。

當(dāng)da>dc時(shí)，剛度值則為(12)式所對(duì)應(yīng)的值。

當(dāng)βd≤da≤dc時(shí)，本文用推導(dǎo)函數(shù)來擬合上述兩部分實(shí)現(xiàn)完整的連續(xù)函數(shù)曲線。函數(shù)應(yīng)該滿足如下條件:

2) 當(dāng)da>dc時(shí)，函數(shù)值應(yīng)等于(12)式計(jì)算剛度值。

當(dāng)da>βd時(shí)，完整的聯(lián)接剛度式為

所以完整的被聯(lián)接結(jié)構(gòu)聯(lián)接剛度可以表示為

(14)

(15)

(16)

(17)

(14)式、(15)式、(16)式、(17)式即是本文推導(dǎo)出的擬合實(shí)際螺栓結(jié)構(gòu)聯(lián)接剛度的計(jì)算公式。依照(14)式～(17)式帶入相應(yīng)的不同參數(shù)螺栓結(jié)構(gòu)，計(jì)算其剛度值變化曲線，并與完全相同的數(shù)值模型曲線比對(duì)，結(jié)果如圖7所示。

圖7(a)、圖7(b)、圖7(c)顯示不同螺栓在不同聯(lián)接長度情況下剛度隨da變化曲線，可以發(fā)現(xiàn)剛度公式與FEA結(jié)果基本一致，誤差為0%～3.16%，在數(shù)值計(jì)算的工程精度范圍內(nèi)。

圖7 螺栓FEA剛度與計(jì)算剛度曲線對(duì)比圖Fig.7 Bolt joint stiffnesses calculated by FEA empirical formulas

圖8(a)是α=1.04、β=1.8、L=40 mm的M10螺栓聯(lián)接結(jié)構(gòu)FEA剛度曲線和本文剛度公式曲線與各文獻(xiàn)計(jì)算的剛度曲線比較圖，圖8(b)是α=1.067、β=1.5、L=60 mm時(shí)M12螺栓聯(lián)接結(jié)構(gòu)FEA剛度曲線和本文剛度公式曲線與各文獻(xiàn)計(jì)算剛度曲線比較圖，通過圖8(a)、圖8(b)計(jì)算結(jié)果的比較?？梢郧逦吹接杀疚膭偠扔?jì)算公式獲得的計(jì)算結(jié)果與FEA結(jié)果曲線的匹配程度比圖8中任何其他文獻(xiàn)計(jì)算結(jié)果都更為接近準(zhǔn)確，因此能有效替代以往的工程計(jì)算方法更好地滿足現(xiàn)代工程設(shè)計(jì)的高精度要求。

圖8 M10、M12的本文剛度公式、FEA與各文獻(xiàn)剛度計(jì)算結(jié)果比較Fig.8 M10, M12 stiffnesses calculated by empirical formulas, FEA and Refs.[5-6]

4 聯(lián)接長度和螺母系數(shù)對(duì)聯(lián)接剛度的影響

前述已經(jīng)推導(dǎo)和論證了不同參數(shù)下螺栓結(jié)構(gòu)聯(lián)接剛度的計(jì)算公式，并通過聯(lián)接剛度隨被聯(lián)接板最大寬度直徑da的變化曲線，證明了本文剛度計(jì)算公式的準(zhǔn)確性。由于最佳預(yù)緊力設(shè)計(jì)原則下最大受力螺牙上的承載應(yīng)力為屈服前應(yīng)力，亦是結(jié)構(gòu)中所承受的最大應(yīng)力。根據(jù)這個(gè)最大應(yīng)力原則設(shè)計(jì)的長效健康螺栓聯(lián)接結(jié)構(gòu)在材料為同一彈性模量E時(shí)，其結(jié)構(gòu)接觸面連接剛度遵循彈性剛度原則，即K=EA/L,其中主要影響因素是最大接觸壓力狀態(tài)的最大接觸面積和壓應(yīng)力強(qiáng)度分布。而文獻(xiàn)[8]和文獻(xiàn)[9]分別通過實(shí)驗(yàn)和仿真對(duì)比，有效證明了二者接觸面積和壓應(yīng)力分布的一致性，進(jìn)而證明了仿真模型接觸剛度的正確性，這也充分闡明了螺栓連接結(jié)構(gòu)的內(nèi)在力學(xué)機(jī)理。以上述理論和模型為基礎(chǔ)推導(dǎo)出的聯(lián)接剛度公式通過與FEA曲線的對(duì)比驗(yàn)證進(jìn)一步證明了其可靠性及準(zhǔn)確度。

圖9(a)和圖9(b)分別表達(dá)了M10螺栓在大寬度(da=100 mm)和有限寬度(da=22 mm)時(shí)剛度受被聯(lián)接件長度L影響下本文公式和FEA得到的結(jié)構(gòu)聯(lián)接剛度曲線，這里聯(lián)接長度L在大寬度時(shí)分別取12 mm、16 mm、20 mm、24 mm和28 mm，在有限寬度時(shí)分別取24 mm、32 mm、40 mm、48 mm和56 mm，以聯(lián)接剛度與直徑比d/L為坐標(biāo)系。圖9(a)、圖9(b)表明，隨著聯(lián)接長度L的增加，剛度將減小。圖9中FEA結(jié)果與剛度公式的誤差為1.29%～4.28%.

圖9 直徑比d/L對(duì)螺栓聯(lián)接剛度的影響Fig.9 Influence of d/L on stiffness

圖10(a)和圖10(b)分別是M10螺栓在被聯(lián)接件大寬度和有限寬度情況下，且大寬度聯(lián)接長度L=20 mm，有限寬度聯(lián)接長度L=60 mm，螺母系數(shù)β分別為1.4、1.6、1.8、2.0、2.2和2.4時(shí)本文公式和FEA所得到的聯(lián)接剛度與螺母系數(shù)β的關(guān)系。從圖10中可看出，隨著螺母系數(shù)的增加，剛度值增大，亦即在一定范圍內(nèi)隨著螺母的增大，被聯(lián)接件的壓應(yīng)力分布變大，相應(yīng)的剛度值變大。FEA結(jié)果和本文公式計(jì)算結(jié)果的誤差為0.04%～3.52%.

圖10 螺母系數(shù)對(duì)螺栓聯(lián)接剛度的影響Fig.10 Influence of nut coefficient on stiffness

5 結(jié)論

本文的螺栓聯(lián)接預(yù)緊力設(shè)計(jì)系列參數(shù)模型，充分考慮了螺栓聯(lián)接結(jié)構(gòu)實(shí)際工作狀態(tài)下材料、預(yù)緊力、摩擦等因素，通過在板與板以及螺母與板之間的各個(gè)接觸面間建立摩擦接觸對(duì)，在已論證各接觸面間壓應(yīng)力分布的有效性和可靠性基礎(chǔ)上[9]，得到系列FEA模型剛度曲線，通過無量綱化研究出新的螺栓結(jié)構(gòu)聯(lián)接剛度計(jì)算公式，并探討了螺栓結(jié)構(gòu)各幾何參數(shù)對(duì)剛度的影響，得到結(jié)論如下：

1)本文以實(shí)際接觸面壓應(yīng)力分布理論為基礎(chǔ)提出的螺栓聯(lián)接剛度計(jì)算公式，與以往文獻(xiàn)的經(jīng)驗(yàn)公式有質(zhì)的區(qū)別，對(duì)比FEA結(jié)果其計(jì)算誤差<3.5%，證明了該計(jì)算公式的有效性和實(shí)際工程應(yīng)用價(jià)值。

2)研究證明螺栓聯(lián)接剛度隨著被聯(lián)接件寬度的增加而增加，且在da>dc后逐步趨近于一個(gè)漸近值。

3)被聯(lián)接件聯(lián)接剛度隨著聯(lián)接長度的增加而減小，隨螺母系數(shù)的增大而增大，而螺栓的孔隙系數(shù)對(duì)結(jié)構(gòu)剛度影響很小。

References)

[1] Shigley J E, Mischke C R, Budynas R G. Mechanical engineering

design[M]. 7th ed. New York:McGraw-Hill, 2004.

[2] Wileman J, Choudhury M, Green I. Computation of member stiffness in bolted connection[J]. ASME Journal of Mechanical Design, 1991,113(2):432-437.

[3] Musto J C, Konkle N R. Computation of member stiffness in the design of bolted joints[J].ASME Journal of Mechanical Design,2006,128(6):1357-1360.

[4] VDI 2230 Systematische berechnung hochbeanspruchter schraubenverbindungen[S]. Berlin und K?ln, Germany:VDI,1977.

[5] VDI 2230 Systematische berechnung hochbeanspruchter schraubenverbindungen[S]. Berlin und K?ln, Germany: VDI, 1986.

[6] VDI 2230 Systematische berechnung hochbeanspruchter schraubenverbindungen zylindrische einschraubenverbindungen-Systematic calculation of high duty bolted joints-Joints with one cylindrical bolt[S]. Düsseldorf, Germany: VDI, 2003.

[7] Setheraman R，Kumar T S. Finite element based member stiffness evaluation of axisymmetric bolted joints[J].ASME Journal of Mechanical Design, 2009,131(1):1-11.

[8] Marshall M B, Lewis R, Dwyer-Joyce R S. Characterisation of contact pressure distribution in bolted joints[J]. Strain,2006,42(1):31-43.

[9] 章伊華,林丹益,楊國玉. 基于接觸理論的螺栓聯(lián)接接觸面力學(xué)特性研究[J].兵工學(xué)報(bào),2015,36(5):946-952. ZHANG Yi-hua, LIN Dan-yi, YANG Guo-yu. Research on mechanical properties of the contact surfaces in bolt joints [J]. Acta Armamentarii,2015,36(5): 946-952.(in Chinese)

[10] 章伊華,侯培海,楊國玉.汽車螺紋聯(lián)接件預(yù)緊力的優(yōu)化設(shè)計(jì)[J].東北林業(yè)大學(xué)學(xué)報(bào),2012,40(11):151-154. ZHANG Yi-hua, HOU Pei-hai, YANG Guo-yu. Optimized design of preload for tread connection of automobile[J]. Journal of Northeast Forestry University, 2012, 40(11):151-154.(in Chinese)

[11] 趙帥,王克明.螺栓聯(lián)接結(jié)構(gòu)振動(dòng)特性有限元分析方法的研究[J].沈陽航空航天大學(xué)學(xué)報(bào),2012,29(2):18-22. ZHAO Shuai,WANG Ke-ming. The study of the dynamic characteristics of a bolt-connected structure by finite element analysis method[J].Journal of Shenyang Aerospace University,2012,29(2):18-22.(in Chinese)

Study of Bolt Joint Stiffness Based on Contact Surface Characteristics

ZHANG Yi-hua1, PANG Kui2, LIN Dan-yi1, MAI Yun-fei3, SUN Xiao-meng3, YANG Guo-yu4

(1.Institute of Automotive CAE Applications, Ningbo Dahongying University, Ningbo 315175，Zhejiang，China;2.Institute of Motorcycle，Lifan Industry (Group) Co., Ltd., Chongqing 400700;3.School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093，China;4.EASi Engineering,Detroit 48201, MI, US)

Bolt joint is widely used in industry because of its reliability and effectiveness. The joint stiffness plays a very important role in bolt joint structure. A finite method of bolt joint structures under the optimal preload is proposed. The joint stiffness is analyzed for different structural parameters in consideration of the nonlinear factors,such as of force, contact property and friction. The effects of bolt pretightening force, thickness and external diameter of joint member, and nut parameters on joint stiffness are studied. The dimensionless stiffness function expression is derived through the experimental verification of contact surface compressive stress distribution and the assumption of joint stiffness continuity, and the model dimensionless stiffness data is fit. The formula for calculating bolt joint stiffness is developed based on the results of finite element analysis (FEA). The effectiveness and accuracy of joint stiffness calculation formula is verified by comparing FEA results with those from previous studies.

solid mechanics; parameter model; contact area; stress distribution of contact surface; joint stiffness; dimensionless

2015-11-02

浙江省自然科學(xué)基金項(xiàng)目(LY12E05012)

章伊華(1962—)，男，副研究員。E-mail:zhyh116@163.com

O343.3；O344.1

A

1000-1093(2017)01-0195-07

10.3969/j.issn.1000-1093.2017.01.026