基于ANSYS的螺栓聯(lián)接有限元計算

劉雪霞， 李杰， 孫萌

（1.大連交通大學(xué)，遼寧大連 116028；2.大連華銳重工集團股份有限公司，遼寧 大連 116013）

0 引言

螺栓聯(lián)接是鋼結(jié)構(gòu)聯(lián)接的主要方式之一，具有施工簡單、拆裝方便等優(yōu)點。根據(jù)受力情況，可將螺栓分為剪力螺栓和拉力螺栓，其中剪力螺栓依靠螺栓桿承擔(dān)剪力和桿件孔壁承擔(dān)壓力來傳力，拉力螺栓則依靠螺桿受拉來傳遞平行于螺桿的外力。目前，設(shè)計人員主要采用設(shè)計規(guī)范中的傳統(tǒng)算法對螺栓聯(lián)接結(jié)構(gòu)進(jìn)行校核，以指導(dǎo)設(shè)計，但傳統(tǒng)算法無法詳細(xì)反映聯(lián)接結(jié)構(gòu)的應(yīng)力分布和變形情況，因此本文借助有限元分析軟件ANSYS，對螺栓聯(lián)接結(jié)構(gòu)進(jìn)行整體模擬分析。

1 有限元計算方法

螺栓聯(lián)接屬于裝配架構(gòu)的接觸問題，在數(shù)值模擬計算中，能夠準(zhǔn)確定義接觸對和施加預(yù)緊力是螺栓聯(lián)接計算的關(guān)鍵。在以往解決類似裝配結(jié)構(gòu)的有限元計算中，通常將結(jié)構(gòu)簡化成獨立結(jié)構(gòu)分別進(jìn)行計算以避免接觸非線性問題，但這種方式的計算精度往往難以保證。

ANSYS軟件作為大型通用有限元分析軟件，在解決高度非線性的接觸問題上具有很強的功能，它提供了多種接觸分析的單元類型，能夠解決所有接觸問題，其求解精度取決于模型接觸區(qū)域網(wǎng)格尺寸及單元接觸剛度值，同時ANSYS軟件還針對螺栓等有預(yù)緊力結(jié)構(gòu)提供了特殊的單元（PRETS179）和求解方法，因此該軟件能夠精確模擬螺栓聯(lián)接結(jié)構(gòu)。

1.1 建立計算模型

選擇螺栓聯(lián)接結(jié)構(gòu)由蓋板、主鋼板和螺栓組（含螺栓、螺母）組成，見圖1中所示。螺栓受到預(yù)緊力和橫向工作載荷。螺栓規(guī)格為M20，5個螺栓橫向均布，間距為80 mm；螺栓孔直徑為φ20 mm；蓋板厚度為10 mm；主鋼板厚度為20 mm；橫向載荷為N=100 kN；螺栓預(yù)緊載荷為P=5 000 N。

圖1 螺栓聯(lián)接結(jié)構(gòu)簡圖

在ANSYS前處理模塊中完成有限元建模，蓋板、主鋼板和螺栓組模型均選用的是八節(jié)點六面體實體單元SOLID45，彈性模量 E=206 MPa，密度為 7 850 kg/m3，泊松比為0.3。網(wǎng)格劃分采用體掃略SWEEP方式生成六面體網(wǎng)格，如圖2中所示，該模型共有單元53 640個，節(jié)點41 950個。

圖2 有限元模型

1.2 接觸定義

ANSYS軟件中支持3種接觸方式：點-點接觸，點-面接觸，面-面接觸。在接觸定義之前首先要確定模型中的哪些部分可能會接觸，若是點接觸，則模型的對應(yīng)組元是節(jié)點，若是面接觸，則模型的對應(yīng)組元是單元，有限元模型通過指定的接觸單元來識別可能的接觸對，接觸單元是覆蓋在分析模型接觸面上的一層單元，有限元模型會通過所設(shè)定的接觸單元來識別可能的接觸對。

本文采用的是三維實體單元建模并確定可能接觸的部分都是面，故選用CONTA174單元和TARGE170單元分別模擬“接觸”面和“目標(biāo)”面。在蓋板和主鋼板之間、螺栓頭與蓋板之間、螺母與蓋板之間、螺栓桿與螺栓孔之間創(chuàng)建面-面接觸對，接觸面摩擦因數(shù)取0.3。

1.3 預(yù)緊力施加

整個結(jié)構(gòu)的單元劃分完成后，定義螺栓的預(yù)緊截面，利用PSMESH命令：定義預(yù)緊截面位置、編號和名稱，并指定預(yù)緊力方向。

施加預(yù)緊力：利用SLOAD命令選擇已定義的預(yù)緊截面，給定預(yù)緊力數(shù)值。預(yù)緊力在第一載荷步（Loadstep1）中施加并計算，然后在第二載荷步（Loadstep2）中鎖定，圖 3所示為其中一螺栓上的預(yù)緊單元。

圖3 預(yù)緊單元（PRETS179）

2 計算結(jié)果分析

2.1 傳統(tǒng)法計算結(jié)果

設(shè)計規(guī)范中指出，剪力螺栓在受力以后，當(dāng)外力不大時，由構(gòu)件間的摩擦力來傳遞外力。當(dāng)外力繼續(xù)增大而超過極限摩擦后，構(gòu)件之間出現(xiàn)相對滑移，螺桿開始接觸構(gòu)件的孔壁而受剪，孔壁則受壓。當(dāng)聯(lián)接處于彈性階段，由于設(shè)計時螺桿與螺孔之間存在著間隙，導(dǎo)致螺栓群中各螺栓承受的剪力也不同。

為考慮到這種不均勻分布，設(shè)計規(guī)范中規(guī)定：當(dāng)聯(lián)接長度l≤15d0時，可認(rèn)為橫向載荷N由每個螺栓平均承受；當(dāng)l1>15d0時，應(yīng)將螺栓的承載力設(shè)計值乘以折減系數(shù)η=1.1-l1/（150d0）≥0.7。式中：l1為兩端螺栓間的中心距離；d0為螺栓孔直徑。

對于本文中選擇的螺栓聯(lián)接結(jié)構(gòu)l1≈15d0，根據(jù)設(shè)計規(guī)范，可認(rèn)為軸心力N由每個螺栓平均承受，每個螺栓上的剪應(yīng)力為 τ=4×100 000÷（5×π×202）=64 MPa。

2.2 有限元計算結(jié)果

有限元數(shù)值模擬計算各螺栓剪應(yīng)力云圖及位移云圖分別見圖4、圖5中所示。蓋板螺栓孔及主鋼板螺栓孔處擠壓應(yīng)力見圖6、圖7所示。

2.3 兩種計算結(jié)果比較

根據(jù)兩種算法結(jié)果，可得出以下結(jié)論：1）在橫向載荷作用下，兩種計算方法得到的剪應(yīng)力分布趨勢一致，即兩端大中間小，近似于按螺栓群中心對稱分布，所以兩種算法所確定的最危險螺栓一致；2）利用設(shè)計規(guī)范公式計算所得的剪應(yīng)力比有限元計算結(jié)果數(shù)值偏大，可以認(rèn)為傳統(tǒng)算法偏于保守。

圖4 螺栓剪應(yīng)力云圖

圖5 螺栓橫向位移云圖

圖6 蓋板螺栓孔擠壓應(yīng)力云圖

圖7 主鋼板螺栓孔擠壓應(yīng)力云圖

3結(jié)語

本文利用有限元軟件ANSYS對整個螺栓聯(lián)接結(jié)構(gòu)進(jìn)行數(shù)值模擬分析，充分考慮結(jié)構(gòu)的接觸特性和應(yīng)力集中等特點，其計算結(jié)果能夠更好地反映螺栓的承載特性，并與設(shè)計規(guī)范中傳統(tǒng)公式的計算結(jié)果互相驗證，進(jìn)一步肯定了兩種方法的正確性。

［1］ GB50017-2003 鋼結(jié)構(gòu)設(shè)計規(guī)范［S］.

［2］ 金晶，吳新躍，鄭建華.螺栓聯(lián)結(jié)預(yù)應(yīng)力施加方法改進(jìn)研究［J］.海軍工程大學(xué)學(xué)報，2010，22（2）：20-24.

［3］ 李會勛，胡迎春，張建中.利用ANSYS模擬螺栓預(yù)緊力的研究［J］.山東科技大學(xué)學(xué)報，2006，25（1）：57-59.

［4］ 葉先磊，史亞杰.ANSYS工程分析軟件應(yīng)用實例［M］.北京：清華大學(xué)出版社，2003.