螺栓預(yù)緊過程的有限元仿真研究

桂學(xué)文 張京東 廖日東

螺栓預(yù)緊過程的有限元仿真研究

桂學(xué)文 張京東 廖日東

(北京理工大學(xué) 機(jī)械與車輛學(xué)院，北京 100081)

預(yù)緊過程直接關(guān)系到螺栓結(jié)構(gòu)的強(qiáng)度計(jì)算和松動(dòng)分析，但是螺栓的預(yù)緊過程是一個(gè)復(fù)雜的非線性問題，目前對于預(yù)緊后螺栓的受力狀態(tài)、接觸狀態(tài)以及不同加載方法的差異等方面的認(rèn)識尚不完善。本文建立了精細(xì)的螺栓六面體網(wǎng)格模型，對轉(zhuǎn)角法和預(yù)緊單元法兩種不同加載方法進(jìn)行了仿真研究。基于轉(zhuǎn)角法研究了扭矩、轉(zhuǎn)角、預(yù)緊力之間的關(guān)系，為一些特殊螺紋連接的扭拉關(guān)系研究奠定了基礎(chǔ)。并對比分析了兩種加載方法下螺栓預(yù)緊后的受力和接觸狀態(tài)。結(jié)果表明，本文采用的轉(zhuǎn)角法加載能夠準(zhǔn)確的模擬螺栓的擰緊過程，預(yù)緊單元法只能模擬軸向力的分布，并不能準(zhǔn)確考慮切向力的影響。

螺栓連接；預(yù)緊過程；扭拉關(guān)系

0 引言

螺栓是一種與螺母配合使用的螺紋緊固件，工程中螺栓預(yù)緊一般是對螺母施加扭矩，使之相對螺桿轉(zhuǎn)動(dòng)一定角度，利用螺母和螺桿上內(nèi)外螺紋的相互作用使螺母和螺桿之間發(fā)生相對軸向運(yùn)動(dòng)，夾緊被連接件從而實(shí)現(xiàn)緊固作用，這就是所謂的轉(zhuǎn)角法。在一些特殊情況下，如超大規(guī)格螺栓的應(yīng)用時(shí)，由于所需扭矩太大不方便施加，也會采用先拉伸螺桿到一定長度后，再將螺母用很小的扭矩?cái)Q緊到位的預(yù)緊方法（簡稱拉伸法）。

螺栓預(yù)緊過程研究一方面是要探明給定結(jié)構(gòu)條件下，達(dá)到給定預(yù)緊力所需螺母的轉(zhuǎn)動(dòng)角度（扭矩）或者螺栓的拉伸長度，以指導(dǎo)安裝工藝設(shè)計(jì)；另一方面是獲得預(yù)緊時(shí)螺栓和螺母的應(yīng)力狀態(tài)和接觸狀態(tài)，以支撐強(qiáng)度計(jì)算和松動(dòng)分析。

但是，由于試驗(yàn)測試的困難，人們對螺栓預(yù)緊狀態(tài)的認(rèn)識以及有關(guān)因素的影響規(guī)律并不深刻。因此，開展螺栓預(yù)緊過程的仿真研究，對加深螺栓預(yù)緊狀態(tài)的認(rèn)識，實(shí)現(xiàn)螺栓結(jié)構(gòu)的預(yù)測設(shè)計(jì)具有重要的理論意義和工程價(jià)值[1]。

已有針對螺栓預(yù)緊過程的研究主要分為對預(yù)緊力的控制和對載荷分布的研究。機(jī)械設(shè)計(jì)教材一般都給出了螺栓結(jié)構(gòu)的扭拉關(guān)系式，指出螺栓的擰緊扭矩與預(yù)緊力近似滿足線性關(guān)系[2]。這個(gè)關(guān)系式是目前螺紋緊固件結(jié)構(gòu)設(shè)計(jì)和裝配設(shè)計(jì)中最重要的依據(jù)。但是，這個(gè)關(guān)系式一般只適用于常規(guī)的三角形螺紋、矩形螺紋，對如今一些廣泛使用的新型緊固件，如自鎖螺母，其收口設(shè)計(jì)會使螺紋變形，產(chǎn)生額外的初始自鎖扭矩[3]，該關(guān)系式就不再適用。

研究有關(guān)因素對螺紋連接結(jié)構(gòu)扭拉關(guān)系的影響，實(shí)現(xiàn)對裝配過程更準(zhǔn)確的控制一直是研究人員努力的方向，如大連理工大學(xué)的陳德安[4]研究了幾何偏差對預(yù)緊力和擰緊扭矩的影響，侯博文[5]考慮了接觸表面形貌對擰緊工藝進(jìn)行研究，得到不同參數(shù)下不同擰緊工藝的誤差范圍，徐衛(wèi)秀等人[6]研究了摩擦因數(shù)、結(jié)構(gòu)材料對扭矩系數(shù)等的影響，為實(shí)際裝配的擰緊工藝提供了理論參考。

對于螺紋結(jié)構(gòu)載荷分布的研究，最重要的工作當(dāng)屬英國國家物理實(shí)驗(yàn)室的Sopwith[7]和東京工業(yè)大學(xué)的山本晃(Yamamoto)[8]給出的螺紋副表面軸向載荷分布的解析結(jié)果。這個(gè)結(jié)果對螺栓受拉螺母受壓的螺紋結(jié)構(gòu)強(qiáng)度設(shè)計(jì)具有重要的啟示意義。但由于該結(jié)果是基于線彈性假設(shè)的，因此難以應(yīng)用在大載荷螺紋結(jié)構(gòu)的設(shè)計(jì)計(jì)算上，另外這個(gè)結(jié)果也沒有涉及螺紋結(jié)構(gòu)所受的橫向載荷。曾攀等人[9]通過光彈實(shí)驗(yàn)，測得螺母的軸向載荷分布，與解析解吻合，但由于光彈實(shí)驗(yàn)材料的限制，只能適用于很小載荷的情況。

而實(shí)際也存在一些螺母承受拉力的情況，如使用螺紋張緊器時(shí)，大連理工大學(xué)的李志彬[10]采用離散化方法得到螺栓螺母同時(shí)受拉時(shí)的載荷分布解。華北電力大學(xué)的江文強(qiáng)等人[11]則給出受橫向載荷時(shí)螺紋副上的橫向力分布，并以此計(jì)算出螺栓連接的臨界松動(dòng)載荷，但缺乏實(shí)驗(yàn)或仿真的驗(yàn)證，且實(shí)際上軸向載荷也會造成螺紋副上的橫向力的不均勻分布，這也會影響臨界松動(dòng)載荷的計(jì)算。有限元仿真可以獲得螺紋連接接觸面的實(shí)時(shí)力學(xué)行為和相對運(yùn)動(dòng)狀態(tài)，更有效的進(jìn)行螺紋連接問題的研究。但螺栓連接結(jié)構(gòu)建模復(fù)雜，時(shí)間成本高。根據(jù)研究的側(cè)重點(diǎn)的不同，通常對模型進(jìn)行不同程度的簡化，對一些關(guān)注重點(diǎn)是整體連接結(jié)構(gòu)的計(jì)算問題，建模時(shí)可以忽略螺紋部分，或者將螺栓簡化成梁單元、剛性單元連接等模型[12]。若進(jìn)行螺栓的強(qiáng)度校核，分析螺紋牙上的承載力分布和變形情況，可以忽略螺旋升角，建立簡化軸對稱模型，將三維問題轉(zhuǎn)化為二維問題[13]。但扭拉關(guān)系、松動(dòng)等問題研究需要考慮螺旋升角產(chǎn)生的扭矩[14]，必須建立考慮升程的三維螺紋連接有限元模型。

但要建立精度和效率俱佳的三維螺紋連接有限元網(wǎng)格模型存在很大的困難。四面體網(wǎng)格劃分的網(wǎng)格質(zhì)量差，將螺紋部分獨(dú)立建模再綁定于螺桿上的方法的計(jì)算結(jié)果不準(zhǔn)確[14]。2008年日本學(xué)者Fukuko[15]提出一種一體化六面體網(wǎng)格劃分方法。該方法利用了螺紋各橫截面的輪廓曲線的一致性，基于輪廓線方程建立二維網(wǎng)格模型，再旋轉(zhuǎn)平移得到不同高度的網(wǎng)格模型，構(gòu)造出完整的六面體網(wǎng)格的螺紋有限元模型。目前學(xué)者建立的精確的螺紋連接模型多采取此種網(wǎng)格劃分方法。

有限元仿真中預(yù)緊力加載可通過多種方式實(shí)現(xiàn)，如溫度載荷法、過盈配合法[16]。這兩種方法需要試算迭代，使用較少，目前多采用的加載方法為預(yù)緊單元法。這些方法可以看作是實(shí)際中拉伸法的仿真模擬。而對實(shí)際的擰緊過程的仿真則需要使用旋轉(zhuǎn)螺母的轉(zhuǎn)角法。但由于計(jì)算效率、仿真結(jié)果的收斂性等原因，轉(zhuǎn)角法仿真的實(shí)現(xiàn)并不簡單，如文獻(xiàn)[17]基于動(dòng)力學(xué)求解器實(shí)現(xiàn)轉(zhuǎn)角法加載，得到預(yù)緊力、扭矩、轉(zhuǎn)角關(guān)系并非線性，預(yù)緊力以及轉(zhuǎn)角增加過程有“臺階”現(xiàn)象，這一仿真結(jié)果并不準(zhǔn)確。目前的研究中經(jīng)常采用預(yù)緊單元法代替擰緊過程的預(yù)緊力加載，如文獻(xiàn)[18]、[19]分別用不同的預(yù)緊力加載方法研究了橫向振動(dòng)條件下的松動(dòng)過程，沒有考慮兩種加載方法的差別，只是為了達(dá)到產(chǎn)生夾緊力的效果。

因此本文利用有限元分析軟件ABAQUS，采用兩種預(yù)緊力加載方法對螺栓的預(yù)緊受力過程進(jìn)行了仿真分析，討論了不同的預(yù)緊力施加方法對結(jié)果的影響，對比分析加載后的受力、接觸狀態(tài)的差異，為今后仿真計(jì)算應(yīng)采取何種預(yù)緊力加載方式提供參考。并基于轉(zhuǎn)角法加載研究了扭矩、轉(zhuǎn)角、預(yù)緊力之間的關(guān)系，得到了沿螺栓螺紋面的接觸應(yīng)力分布，螺紋副的軸向以及橫向載荷的分布規(guī)律。

1 螺栓連接結(jié)構(gòu)有限元建模

本文研究的是普通的三角形螺紋連接結(jié)構(gòu)，其示意圖如圖1所示，主要的尺寸參數(shù)見表1。該系統(tǒng)由螺栓、螺母和兩塊被夾緊件組成。其中螺栓規(guī)格為M8，內(nèi)外螺紋間隙0.1mm，所有材料均為合金鋼，彈性模量為200GPa，泊松比為0.3，屈服強(qiáng)度為1080MPa，切線模量為1000MPa。

借助課題組開發(fā)的螺紋結(jié)構(gòu)六面體有限元網(wǎng)格生成軟件[20-21]，建立了試驗(yàn)螺栓緊固系統(tǒng)的有限元網(wǎng)格模型，其中單元數(shù)為194400，節(jié)點(diǎn)數(shù)為208787。圖2給出了整體有限元網(wǎng)格模型和螺栓、螺母螺紋部分的局部網(wǎng)格模型。模型中存在四對接觸面：螺栓和螺母之間螺紋副的接觸面、螺母和被連接件之間的接觸面（以下簡稱支承面）、螺栓頭和被連接件之間的接觸面以及兩塊被夾緊件之間的接觸面。對這四對接觸面均分別定義法向和切向接觸模型，其中摩擦因數(shù)取為0.15。

圖1 螺栓連接結(jié)構(gòu)示意圖

Fig.1 Schematic diagram of bolted joints

圖2 螺栓連接結(jié)構(gòu)有限元模型

Fig.2 Finite element model of bolted joints

螺栓的擰緊過程可以看作是一個(gè)準(zhǔn)靜態(tài)過程，因此采用靜力學(xué)求解器。為模擬螺母擰緊過程，計(jì)算時(shí)采用了兩步加載。第一載荷步施加順時(shí)針旋轉(zhuǎn)角0.5rad以實(shí)現(xiàn)螺母的擰緊作用，達(dá)到給定轉(zhuǎn)角后保持一段時(shí)間，然后在第二載荷步釋放擰緊時(shí)的加載力矩，線性卸載至0。預(yù)緊單元法需要在光桿段設(shè)置一截面為預(yù)拉截面，定義拉伸方向，在截面處生成預(yù)緊單元，指定螺栓的預(yù)緊力即可[22]。接觸設(shè)置與轉(zhuǎn)角法相同，將預(yù)緊單元法的預(yù)緊載荷與轉(zhuǎn)角法加載得到的預(yù)緊力控制一致。下面對兩種加載方法的結(jié)果進(jìn)行討論分析。

表1 主要的尺寸參數(shù)

2 轉(zhuǎn)角法分析

2.1 扭拉關(guān)系分析

擰緊過程中加載扭矩與產(chǎn)生的預(yù)緊力之間的關(guān)系是螺栓預(yù)緊過程的一個(gè)重要特征，根據(jù)螺紋結(jié)構(gòu)的幾何特點(diǎn)和力學(xué)分析，可以將該過程等效為斜面滑塊模型，得到扭拉關(guān)系式(1)，其中為扭矩系數(shù)[23]

假設(shè)支承面上的壓力為均勻分布，則有

其中，h為螺栓孔內(nèi)徑、w為支承面的接觸外徑。

根據(jù)理論分析可知，扭矩與預(yù)緊力近似成線性關(guān)系。有限元仿真結(jié)果如圖3所示，在加載段扭矩線性加載，預(yù)緊力也線性增加，加載段結(jié)束達(dá)到最大20.84kN；保持段轉(zhuǎn)角保持不變，預(yù)緊力也保持穩(wěn)定；載荷步2的卸載段，轉(zhuǎn)角略有下降，即螺母輕微回轉(zhuǎn)，預(yù)緊力略有下降，最終為20.76kN。

圖3 隨加載歷程的預(yù)緊力和轉(zhuǎn)角變化

Fig.3 The preload and rotation angle change with the loading process

2.2 擰緊角分析

實(shí)際裝配過程中，由于各接觸面有一個(gè)初始貼合的過程，此時(shí)擰緊角與預(yù)緊力呈非線性關(guān)系。本文只考慮貼合過程后的線性段。由于螺栓頭部固定，擰緊角度即為螺母的旋轉(zhuǎn)角度。預(yù)緊過程螺栓受拉，螺母和被夾緊件受壓，由于螺母的剛度相對較大，可忽略螺母的壓縮量，則有

式中，b、c分別為螺栓和被連接件的剛度。

根據(jù)德國工程師協(xié)會VDI2230標(biāo)準(zhǔn)[24]，可以將螺栓等效為幾段圓柱串聯(lián)，如圖1所示，可分為螺栓頭、螺栓光桿段、螺紋未嚙合段、螺紋嚙合段，對每段分別計(jì)算剛度再求得螺栓總剛度為

其中，b為螺栓的彈性模量，A為各段的等效橫截面積。被連接件的剛度受其總厚度、直徑、半頂角的影響，其壓應(yīng)力在一定范圍內(nèi)呈應(yīng)力錐分布，其半頂角及剛度計(jì)算公式為

其中，為半頂角，為被夾緊件的外徑，0為被連接件的彈性模量。

剛度和螺母擰緊角度的仿真和理論解對比見表2。通過仿真計(jì)算出的b、c可以反求出一個(gè)擰緊角度，但會與加載設(shè)置的擰緊角度有所不同。這是由于螺栓和被連接件各位置的變形并不相同，仿真得到的計(jì)算剛度是在主要變形位置取平均所得，存在一定誤差。

表2 螺栓連接剛度和螺母擰緊角度

3 兩種預(yù)緊力加載方法結(jié)果對比

3.1 扭矩分量對比

預(yù)緊單元法與轉(zhuǎn)角法的受力有明顯差異，下面以斜面滑塊模型，對螺母進(jìn)行受力分析，如圖4所示。

圖4 兩種預(yù)緊方式的受力分析

Fig.4 Force analysis of two preloading methods

轉(zhuǎn)角法模擬實(shí)際轉(zhuǎn)動(dòng)過程，螺栓相對螺母向上運(yùn)動(dòng)，螺紋面上的扭矩分量t、p方向相同。且螺紋面上的摩擦為滑動(dòng)摩擦力，方向與運(yùn)動(dòng)方向相反，摩擦力和壓力的向分量tsz、tnz方向相反。而預(yù)緊單元法螺栓螺母沒有發(fā)生宏觀的相對滑動(dòng)，螺栓相對螺母有向下運(yùn)動(dòng)的趨勢，此時(shí)的t與p方向相反，大小相等，合扭矩為0。摩擦力和壓力的向分量方向相同。由于摩擦力和壓力的向分量合力即為預(yù)緊力，則預(yù)緊單元法壓力的向分量小于轉(zhuǎn)角法，其法向接觸應(yīng)力小于轉(zhuǎn)角法。

對轉(zhuǎn)角法螺母所受扭矩進(jìn)行分析，根據(jù)扭矩系數(shù)對應(yīng)的三項(xiàng)，可以將擰緊力矩分解為三個(gè)分力矩：支承面上的摩擦力矩w，螺紋配合面上的摩擦力矩t，螺紋配合面上的正壓力所形成的力矩p，也即由于螺距產(chǎn)生的扭矩[25]

有限元軟件可以直接提取得到接觸面的法向和切向扭矩，仿真得到加載扭矩和各分力矩隨加載歷程的變化如圖5所示，可以看到在載荷步1的加載段，各分力矩與加載扭矩成線性關(guān)系。當(dāng)加載扭矩逐漸卸載時(shí)，Tw將明顯下降。當(dāng)Tf

Fig.5 Torque variation under rotation angle method

載荷步1結(jié)束時(shí)刻仿真與理論值對比見表3，此時(shí)加載扭矩為34.89N·m。仿真結(jié)果與理論解相對誤差在3%以內(nèi)。

表3 擰緊過程仿真與理論值對比

仿真得到預(yù)緊單元法的各扭矩如圖6所示，加載完成時(shí)p=3.81N·m，略小于轉(zhuǎn)角法。由p可以求出摩擦力和壓力的向分量

得

而對轉(zhuǎn)角法，這兩個(gè)分量分別為21kN，-0.24kN。

Fig.6 Torque components of preload element method

3.2 應(yīng)力分布對比

轉(zhuǎn)角法卸載階段結(jié)束后的應(yīng)力如圖7所示，圖7 a)為螺栓的VonMises應(yīng)力分布，縱向應(yīng)力分布近似為軸對稱，螺栓螺桿中心區(qū)域應(yīng)力較小，在螺紋牙根部有明顯的應(yīng)力集中現(xiàn)象。圖8、圖9為螺紋牙和橫截面上的應(yīng)力分布，可以看出螺紋小徑應(yīng)力較大，即對應(yīng)螺栓的螺紋牙根部，而螺栓的螺紋牙頂、螺母的螺紋牙根部應(yīng)力較小。結(jié)合螺栓螺母的裝配位置，以支承面為起始位置，沿軸向每移動(dòng)一個(gè)螺距距離記為一圈嚙合螺紋。螺栓最下面幾圈的螺紋牙根部等效應(yīng)力最大，這些螺紋包括未嚙合螺紋和嚙合螺紋的前幾圈，未嚙合螺紋牙不發(fā)生接觸，嚙合螺紋的前三圈承力較大，這些位置的軸向拉力都較大，導(dǎo)致VonMises應(yīng)力較大，是容易發(fā)生強(qiáng)度破壞的位置。

圖7 轉(zhuǎn)角法的應(yīng)力分布

Fig.7 Stress distribution of rotation angle method

圖8 螺紋牙上的應(yīng)力

Fig.8 Stress distribution on thread

圖9 橫截面上的應(yīng)力分布

Fig.9 Stress distribution on cross section

轉(zhuǎn)角法的法向接觸應(yīng)力分布如圖10所示，可以看到只有螺栓的螺紋下表面有接觸應(yīng)力，只有該區(qū)域和螺母的螺紋上表面發(fā)生了接觸，其他區(qū)域并沒有接觸。如圖11定義一條沿螺紋面螺旋上升的路徑，得到該路徑上節(jié)點(diǎn)的法向接觸應(yīng)力分布曲線，如圖12所示。第一圈螺紋開始接觸應(yīng)力較小，并迅速增大達(dá)到最大值，之后隨著軸向位置增加，法向接觸應(yīng)力不斷減小。

圖10 法向接觸應(yīng)力分布

Fig.10 Normal contact stress distribution

圖11 路徑定義

Fig.11 Path definition

圖12 定義路徑上的接觸應(yīng)力分布

Fig.12 Contact stress distribution on the path

預(yù)緊單元法的VonMises等效應(yīng)力、扭轉(zhuǎn)切應(yīng)力云圖如圖13所示。與圖7對比可知，預(yù)緊單元法的扭轉(zhuǎn)切應(yīng)力遠(yuǎn)遠(yuǎn)小于轉(zhuǎn)角法，最大VonMises應(yīng)力略大于轉(zhuǎn)角法。但由于軸向拉力是主要影響因素，VonMises應(yīng)力相差不大。對于一般的以合金鋼為材料的螺栓，可以采用第四強(qiáng)度準(zhǔn)則進(jìn)行強(qiáng)度校核分析，因此可以使用預(yù)緊單元法進(jìn)行強(qiáng)度計(jì)算，提高計(jì)算效率。但對于松動(dòng)問題，必須了解兩種方法的加載方式帶來的差異，預(yù)緊單元法不能準(zhǔn)確模擬螺栓的拉扭復(fù)合作用，受力狀態(tài)、接觸狀態(tài)的不同會影響松動(dòng)過程的計(jì)算。

圖13 預(yù)緊單元法的應(yīng)力分布

Fig.13 Stress distribution of preloaded element method

3.3 載荷分布對比

由Yamamoto法[7]得到螺栓的螺紋嚙合部分的軸向載荷分布公式為

式中

式中，b、n為螺栓和螺母的橫截面面積；b、n為螺栓和螺母的彈性模量，b、n分別為外螺紋和內(nèi)螺紋的剛度系數(shù)。提取不同橫截面上的軸向力，相距一個(gè)螺距的橫截面的軸向力之差即為一圈螺紋牙上的承力。對比各圈螺紋牙的載荷占比如圖14所示。同理可以得到不同加載方式下軸向載荷的分布，如圖15所示。

圖14 軸向載荷分布占比

Fig.14 Proportion of axial load distribution

圖15 軸向載荷分布

Fig.15 Axial load distribution

可以看出，兩種加載方式嚙合螺紋承受的軸向載荷分布相似，沿軸向先迅速增大，達(dá)到峰值后逐漸減小。在第二圈載荷占比達(dá)到最高，但第一圈載荷占比明顯低于理論值。文獻(xiàn)[26]中也得到了相似結(jié)論。這是因?yàn)槲覀円灾С忻鏋槠鹗键c(diǎn)，由圖10、圖12所示的接觸應(yīng)力分布可知，第一圈螺紋起始段的接觸并不完全，只有部分面積有接觸壓力，且接觸應(yīng)力沿向迅速增大。起始段嚙合面積小，接觸應(yīng)力小，承受的軸向載荷偏小，而理論認(rèn)為各圈螺紋的嚙合狀態(tài)相同，導(dǎo)致載荷占比第一圈最大。后幾牙的分布規(guī)律相似，占比都不斷減小。前三牙承載仍然是軸向載荷的主要占比。

同理，得到橫截面上的，方向的受力分布，如圖16所示，圖中F、F分別為單位長度橫截面上所受的沿和方向的分力，F為其合力。預(yù)緊單元法的面力遠(yuǎn)小于轉(zhuǎn)角法，兩種方法的曲線都呈衰減的簡諧分布，一個(gè)周期對應(yīng)一個(gè)螺距距離。合力F變化趨勢與軸向載荷類似，隨著軸向位置增加先增大后減小，但方向隨位置不斷變化。分力分布的周期性與螺紋形狀的有關(guān)。每個(gè)高度只有螺紋嚙合部分承力。而從橫向看，由于螺旋升角的影響，螺紋的橫截面并不是嚴(yán)格繞軸的軸對稱，橫截面是由小徑、中徑、大徑三段圓弧組成，小徑的位置隨著軸向位置的增大繞著軸轉(zhuǎn)動(dòng)?？梢园l(fā)現(xiàn)其合力方向指向該高度上的螺牙小徑所在方向，因此合力方向也就隨之轉(zhuǎn)動(dòng)。由于力的簡諧分布，使得整個(gè)螺紋牙的橫向力不斷相互平衡抵消，而不是像軸向載荷累加，因此對于整個(gè)螺栓，其所受橫向力很小，只有軸向拉力是主要因素。

圖16 橫截面內(nèi)的橫向力分布

Fig.16 Transverse force distribution in cross section

4 結(jié)束語

螺栓的預(yù)緊過程是研究螺栓強(qiáng)度和松動(dòng)過程的關(guān)鍵。本文考慮了螺紋細(xì)節(jié)，建立了精細(xì)的六面體結(jié)構(gòu)化網(wǎng)格模型，對比了兩種不同的預(yù)緊力加載方式，結(jié)論如下：1）基于轉(zhuǎn)角法仿真模型，可以準(zhǔn)確模擬標(biāo)準(zhǔn)螺栓的擰緊過程，得到螺栓受力情況、應(yīng)力分布以及變形情況，并得到扭矩-預(yù)緊力-轉(zhuǎn)角三個(gè)工程中的重要參數(shù)之間的關(guān)系，與經(jīng)典扭拉關(guān)系式吻合，為今后一些特殊螺紋連接（如高鎖螺母、楔形螺母）的扭拉關(guān)系、轉(zhuǎn)角關(guān)系研究奠定了基礎(chǔ)。2）對比了轉(zhuǎn)角法和預(yù)緊單元法，結(jié)果表明對于強(qiáng)度校核，為了提高計(jì)算效率可以采用預(yù)緊單元法，但是預(yù)緊單元法沒有準(zhǔn)確考慮螺栓預(yù)緊過程中的切向載荷作用，兩種方法的受力狀態(tài)有所差別，對于松動(dòng)問題的研究，仍有必要基于轉(zhuǎn)角法得到螺栓螺母的初始受力狀態(tài)和接觸狀態(tài)。

[1] 王楨, 張京東, 丁曉宇, 等.螺紋連接結(jié)構(gòu)松動(dòng)分析方法研究[J].強(qiáng)度與環(huán)境, 2020, 47(6): 8-14.[Wang Zhen, Zhang Jingdong, Ding Xiaoyu, et al.Technique research on looseness analysis of threaded fasteners [J].Structure & Environment Engineering, 2020, 47(6): 8-14.]

[2] 濮良貴, 陳國定, 吳立言.機(jī)械設(shè)計(jì)(第九版)[M].北京: 高等教育出版社, 2013.

[3] 吳南星, 成飛, 廖達(dá)海, 等.高鎖螺栓連接力學(xué)參數(shù)關(guān)系的實(shí)驗(yàn)研究[J].機(jī)床與液壓, 2017, 45(1):77-80.[Wu Nanxing, Cheng Fei, Liao Dahai, et al.Experimental study on the relationship between the mechanical parameters of hi-lock bolt connection [J].Machine Tool & Hydraulics, 2017,45(1): 77-80.]

[4] Chen D, Ma Y, Hou B.Tightening behavior of bolted joint with non-parallel bearing surface[J].International Journal of Mechanical Sciences, 2019: 153-154.

[5] 侯博文.真實(shí)接觸表面形貌下螺栓連接結(jié)構(gòu)的力學(xué)性能[D].大連理工大學(xué), 2020.

[6] 徐衛(wèi)秀, 王淑范, 楊帆, 等.考慮螺紋細(xì)節(jié)的螺栓預(yù)緊過程仿真分析研究[J].宇航總體技術(shù), 2018, 2(5): 50-56.[Xu Weixiu, Wang Shufan, Yang Fan, et al.Study of bolt pre-tightening process based on simulation with thread details [J].Astronautical Systems Engineering Technology, 2018, 2(5): 50-56.]

[7] Sopwith D G.The distribution of load in screw threads[J].Iust.Mech.Engrs.Appl.Mech.Proc,1948, 159: 373-383.

[8] 山本晃.螺紋聯(lián)接的理論與計(jì)算[M].上海: 上?？茖W(xué)技術(shù)文獻(xiàn)出版社, 1984.

[9] 陳海平, 曾攀, 方剛, 等.螺紋副承載的分布規(guī)律[J].機(jī)械工程學(xué)報(bào), 2010,46(9): 171-178.[Chen Haiping, Zeng Pan, Fang Gang, et al.Load distribution of bolted joint [J].Journal of Mechanical Engineering, 2010, 46(9): 171-178.]

[10] 李志彬.螺栓連接結(jié)構(gòu)松動(dòng)失效機(jī)理研究[D].大連理工大學(xué), 2019.

[11] 江文強(qiáng), 墨澤, 安利強(qiáng), 等.考慮螺紋柔性的螺栓連接臨界松動(dòng)載荷計(jì)算方法[J].機(jī)械工程學(xué)報(bào), 2020, 56(15): 238-248.[Jiang Wenqiang, Mo Ze, An Liqiang, et al.Computing method of bolted joint critical loosening load with flexible thread [J].Journal of Mechanical Engineering, 2020, 56(15): 238-248.]

[12] 賀李平, 任雪梅.螺栓聯(lián)接的有限元建模方法研究[J].北京理工大學(xué)學(xué)報(bào),2020,40(12):1275-1281.[He Liping, Ren Xuemei.Study on finite element modeling method of bolted joints [J].Transactions of Beijing Institute of Technology, 2020, 40(12):1275-1281.]

[13] Liao Ridong, Sun Yujuan, Zhang Weizheng.Nonlinear analysis of axial-load and stress distribution for threaded connection [J].Chinese Journal of Mechanical Engineering, 2009, 22(6): 869-875.

[14] Gong Hao, Liu Jianhua, Ding Xiaoyu.Study on the mechanism of preload decrease of bolted joints subjected to transversal vibration loading [J].Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture, 2019, 233(12): 2320-2329.

[15] Fukuoka Toshimichi, Nomura Masataka.Proposition of helical thread modeling with accurate geometry and finite element analysis[J].Journal of Pressure Vessel Technology, 2008, 130(1).

[16] 徐超, 余紹蓉, 鄭曉亞, 等.機(jī)械螺栓法蘭連接的有限元力學(xué)模型比較研究[J].機(jī)械設(shè)計(jì)與制造, 2009(6): 37-39.[Xu Chao, Yu Shaorong, Zheng Xiaoya, et al.Comparative evaluation of finite element models for mechanical structures with bolted joints [J].Machinery Design & Manufacture, 2009(6): 37-39.]

[17] 王燦, 丁振軍, 韓帥帥, 等.基于精確螺紋建模的螺栓緊固工藝研究[J].工程機(jī)械, 2020, 51(2): 31-39.[Wang Can, Ding Zhenjun, Han Shuaishuai, et al.Research of bolt fastening process based on accurate thread modeling [J].Construction Machinery and Equipment, 2020, 51(2): 31-39.]

[18] 王傳華, 王文莉, 林清源, 等.基于精確建模的橫向振動(dòng)工況下螺栓松動(dòng)機(jī)理研究[J].宇航總體技術(shù), 2018, 2(4): 43-51,64.[Wang Chuanhua, Wang Wenli, Lin Qingyuan, et al.Study on bolt loosening mechanism in vibration condition based on precise model [J].Astronautical Systems Engineering Technology, 2018, 2(4): 43-51,64.]

[19] 胡陽, 姜東, 王旻睿, 等.橫向載荷作用下螺栓連接松動(dòng)過程研究[J].振動(dòng)、測試與診斷,2020,40(6):1091-1098.[Hu Yang, Jiang Dong, Wang Minrui, et al.Study on loosening process of bolted joints under transverse load [J].Journal of Vibration, Measurement & Diagnosis, 2020, 40(6): 1091-1098.]

[20] 侯世遠(yuǎn).螺紋聯(lián)接松動(dòng)機(jī)理研究[D].北京理工大學(xué),2015.

[21] 鄔杰.螺紋連接結(jié)構(gòu)橫向松動(dòng)壽命預(yù)測模型研究[D].北京理工大學(xué), 2019.

[22] 陳巖.螺栓松動(dòng)的失效機(jī)理以及對整體結(jié)構(gòu)力學(xué)行為的影響研究[D].大連理工大學(xué),2019.

[23] Motosh Nabil.Load Distribution on threads of titanium tension nuts and steel bolts[J].Journal of Engineering for Industry, 1975, 97(1): 162-166.

[24] Ingenieure V D .Systematic calculation of high duty bolted joints joints with one cylindrical bolt [J].Berlin : VDI,2003:171.

[25] Nassar Sayed A, Yang Xianjie.Novel formulation of the tightening and breakaway torque components in threaded fasteners [J].Journal of Pressure Vessel Technology, 2007, 129(4): 653-663.

[26] 潘瑞鵬.摩擦系數(shù)對螺栓連接緊固性能的影響研究[D].北京理工大學(xué), 2020.

Finite Element Simulation of Bolt Pre-Tightening Process

GUI Xue-wen ZHANG Jing-dong LIAO Ri-dong

(School of Mechanical and Vehicular Engineering，Beijing Institute of Technology, Beijing 100081, China)

The pre-tightening process is directly related to the strength calculation and loosening analysis of bolted joints.However, the pre-tightening process of bolts is a complex nonlinear problem.At present, the understanding of the stress state, contact state and the differences of different loading methods of bolts after pre-tightening is still unclear.In this paper, a fine bolt hexahedron mesh model is established, and two different loading methods, rotation angle method and preload element method, are simulated.Based on the rotation angle method, the relationship among torque, angle and preload is studied, which lays a foundation for the study of some special threaded connections.The stress and contact state of bolts after preloading under the two loading methods are compared and analyzed.The results show that the rotation angle loading method used in this paper can accurately simulate the tightening process of bolts, and the preloading element method can only simulate the distribution of axial force, but cannot accurately consider the influence of tangential force.

Bolted joints; Pre-tightening process; Relationship between the preload and tightening torque

TH122

A

1006-3919(2021)06-0009-08

10.19447/j.cnki.11-1773/v.2021.06.002

2021-06-23；

2021-09-17

科學(xué)挑戰(zhàn)專題項(xiàng)目（TZ018007）

桂學(xué)文（1998—），男，碩士研究生，研究方向：動(dòng)力機(jī)械結(jié)構(gòu)強(qiáng)度；（100081）北京市海淀區(qū)中關(guān)村南大街5號北京理工大學(xué)發(fā)動(dòng)機(jī)西樓301.