馬 輝，高 昂，楊天瑞，官 宏

(東北大學(xué) 機(jī)械工程與自動(dòng)化學(xué)院，遼寧 沈陽 110819)

含止口配合的螺栓連接結(jié)構(gòu)由于其結(jié)構(gòu)簡單、安裝便捷、定心性好等特點(diǎn)，被廣泛應(yīng)用在管道、法蘭以及航空航天薄壁殼等連接結(jié)構(gòu)當(dāng)中.由于螺栓和止口結(jié)構(gòu)的存在，其結(jié)構(gòu)之間通常有多個(gè)配合接觸面，因此，結(jié)構(gòu)本身呈現(xiàn)出非連續(xù)性的特點(diǎn)，其剛度也會(huì)由于結(jié)構(gòu)的不連續(xù)性而產(chǎn)生剛度損失.

近年來，對(duì)于螺栓連接結(jié)構(gòu)的建模方法，國內(nèi)外學(xué)者展開了大量研究.Kim等[1]討論了板類連接結(jié)構(gòu)中螺栓的建模方法，分別為實(shí)體螺栓建模、蜘蛛螺栓建模、耦合螺栓建模和無螺栓建模，耦合螺栓模型在法蘭盤連接結(jié)構(gòu)之間具有更好的實(shí)用性.Luan等[2]針對(duì)圓柱殼法蘭連接結(jié)構(gòu)，通過不同剛度的線性彈簧來模擬螺栓連接結(jié)構(gòu).Yao等[3]針對(duì)航空發(fā)動(dòng)機(jī)中的螺栓結(jié)構(gòu)，采用薄層單元法進(jìn)行建模，給出關(guān)于該方法參數(shù)化建模的原理.在此基礎(chǔ)上，Ma等[4]利用試驗(yàn)方法獲得螺栓連接的切向剛度以及損耗因子，利用參數(shù)擬合的方法，建立了改進(jìn)的薄層單元模型，以此來模擬航空發(fā)動(dòng)機(jī)中的連接結(jié)構(gòu).對(duì)于螺栓連接在結(jié)構(gòu)剛度特性方面所做的研究中，Luan等[5]針對(duì)螺栓法蘭連接結(jié)構(gòu)，從靜力學(xué)的角度分別研究了軸向以及彎曲的非線性剛度特性.Wang等[6]基于剛度損失機(jī)理，提出了一種等效計(jì)算方法，通過ANSYS驗(yàn)證了方法的有效性，并且分析了螺栓連接剛度損失的影響因素.

對(duì)于含止口的螺栓連接結(jié)構(gòu)，一些學(xué)者也做了相關(guān)的研究.郭文新等[7]建立了含有止口配合的螺栓連接結(jié)構(gòu)有限元模型，通過用實(shí)驗(yàn)得到的剛度-變形量曲線來驗(yàn)證仿真模型的正確性.Liu等[8]基于帶有止口螺栓連接的法蘭結(jié)構(gòu)，采用有限元方法分析了結(jié)構(gòu)的連接剛度，對(duì)接觸面應(yīng)力分布和危險(xiǎn)點(diǎn)應(yīng)力的接觸狀態(tài)進(jìn)行討論.為了更好地研究止口螺栓結(jié)構(gòu)的連接非線性問題，Campos等[9]介紹了如何確定內(nèi)、外配合的薄壁圓筒形結(jié)構(gòu)件徑向和環(huán)向應(yīng)力的簡化方法，對(duì)模擬止口過盈提供了一定的指導(dǎo).Liu等[10]基于ANSYS完全瞬態(tài)分析法分析了含止口連接結(jié)構(gòu)的剛度非線性，通過施加簡諧載荷提取位移數(shù)據(jù)繪制出遲滯曲線.在此基礎(chǔ)上，王志等[11]通過利用詹金斯單元和彈簧單元并聯(lián)來簡化含止口的螺栓法蘭連接結(jié)構(gòu)，并分析了止口在接觸界面之間的滑移變形機(jī)理.李倫緒等[12]針對(duì)含止口的螺栓連接結(jié)構(gòu)，利用有限元數(shù)值仿真方法分析結(jié)構(gòu)非線性特征以及形成機(jī)理，并得出了結(jié)構(gòu)在拉伸載荷作用下，其剛度呈分段非線性的變化規(guī)律，以及在彎矩載荷作用下其剛度具有倒“S”形非線性的變化規(guī)律.對(duì)于含止口螺栓連接結(jié)構(gòu)的設(shè)計(jì)方面，岳偉等[13]基于止口結(jié)構(gòu)分析了結(jié)構(gòu)的受力和變形分析，利用有限元進(jìn)行非線性接觸分析，提出了可拆卸轉(zhuǎn)子穩(wěn)健性的設(shè)計(jì)方法.

基于上述文獻(xiàn)可知，目前，國內(nèi)外學(xué)者對(duì)于螺栓連接結(jié)構(gòu)的研究已較為成熟，特別是在建模方法、連接特性和接觸特性等方面的研究內(nèi)容較為豐富.與此同時(shí)，由于止口螺栓連接結(jié)構(gòu)的廣泛應(yīng)用，部分學(xué)者也針對(duì)止口結(jié)構(gòu)做了一定的研究，雖然在建模過程學(xué)者們考慮了止口結(jié)構(gòu)，但并沒有針對(duì)止口螺栓結(jié)構(gòu)在接觸特性和彎曲剛度特性的研究上進(jìn)行深入而詳盡的探討.針對(duì)這一問題，本文基于含止口配合的螺栓連接法蘭盤結(jié)構(gòu)，考慮接觸非線性，利用ANSYS仿真軟件，分別從多個(gè)角度分析了含止口結(jié)構(gòu)的接觸特性和彎曲剛度變化規(guī)律，為后續(xù)含止口結(jié)構(gòu)的結(jié)構(gòu)設(shè)計(jì)和工程應(yīng)用提供了一定的參考.

1 接觸及彎曲剛度基本理論

1.1 接觸基本理論

在含止口的螺栓連接法蘭盤結(jié)構(gòu)當(dāng)中，接觸區(qū)域非線性求解對(duì)于求解效率和精度影響較大，選擇合理的接觸算法是提高求解準(zhǔn)確性和效率的關(guān)鍵.結(jié)合ANSYS有限元軟件中常見三種接觸問題數(shù)值計(jì)算方法的優(yōu)缺點(diǎn)[14]，本文對(duì)建模過程中所涉及到的接觸問題均采用增廣拉格朗日乘子法進(jìn)行求解.

關(guān)于接觸問題的求解過程，是通過對(duì)目標(biāo)節(jié)點(diǎn)和接觸節(jié)點(diǎn)之間的侵入容差來進(jìn)行收斂性的判斷，通過多次迭代后，當(dāng)侵入容差小于設(shè)定值時(shí)，求解結(jié)束.基于將接觸問題等效成求解區(qū)域的位移場，使得系統(tǒng)的勢能在接觸邊界條件下的約束達(dá)到最小，其數(shù)學(xué)表達(dá)式為

(1)

由于接觸的高度非線性，修正系統(tǒng)的總能量與間隙函數(shù)和拉格朗日乘子矩陣緊密相關(guān)，可將式(1)中的系統(tǒng)補(bǔ)償勢能進(jìn)行展開，展開式為

(2)

(3)

1.2 彎曲剛度基本理論

航空發(fā)動(dòng)機(jī)的彎曲剛度是衡量發(fā)動(dòng)機(jī)抵抗變形能力的重要參數(shù)，結(jié)構(gòu)通過螺栓預(yù)緊力保證連接的可靠性，止口配合面在承擔(dān)彎矩載荷的同時(shí)也減小了連接界面在徑向和周向發(fā)生相對(duì)滑移的可能性.由于止口結(jié)構(gòu)的存在，當(dāng)結(jié)構(gòu)存在橫向載荷或者彎矩載荷時(shí)，對(duì)應(yīng)結(jié)構(gòu)會(huì)隨著外部載荷的作用產(chǎn)生相應(yīng)的變形，如圖1所示，即一端出現(xiàn)橫向位移 以及轉(zhuǎn)角.剛度的計(jì)算表達(dá)式為

圖1 含止口的螺栓連接法蘭盤結(jié)構(gòu)變形圖

(4)

(5)

2 模型建立及靜力學(xué)特性分析

首先基于文獻(xiàn)幾何參數(shù)建立含止口的螺栓連接法蘭盤結(jié)構(gòu)實(shí)體有限元模型，并進(jìn)行了驗(yàn)證.與此同時(shí)，建立了無止口結(jié)構(gòu)的實(shí)體有限元模型，通過施加橫向載荷來對(duì)比含止口和無止口結(jié)構(gòu)對(duì)接觸特性和彎曲剛度的影響.

2.1 實(shí)體單元模型建立

基于文獻(xiàn)[8]的幾何參數(shù)建立了含止口和無止口的螺栓連接法蘭盤結(jié)構(gòu)精細(xì)化有限元模型，如圖2所示.為了驗(yàn)證本文模型的有效性，建模過程中采用與文獻(xiàn)相同的Solid186實(shí)體單元；各個(gè)接觸面之間均采用Conta174接觸單元和Targe170目標(biāo)單元，且摩擦系數(shù)均為0.15；對(duì)于螺栓預(yù)緊力的模擬則采用Prets179單元定義預(yù)緊力截面.在有限元模型中，其邊界條件設(shè)置為：左端A面為固定約束，右端B面與無質(zhì)量的中心點(diǎn)(Mass21單元)剛性耦合于P點(diǎn)，并在P點(diǎn)施加20 kN橫向載荷，分多步載荷進(jìn)行加載，并將其剛度值繪制成如圖3所示的曲線.

圖2 有限元模型

圖3 剛度-變形量曲線

如圖3所示，文獻(xiàn)[8]和本文在不同工況下的剛度-變形量曲線，對(duì)應(yīng)工況如表1所示.由圖可知，本文所對(duì)應(yīng)的4個(gè)工況曲線與文獻(xiàn)[8]曲線整體趨勢保持一致，但在數(shù)值上略有差異，其原因是，文獻(xiàn)[8]模型中雖含有止口結(jié)構(gòu)，但過盈量為0，且在止口配合界面采用默認(rèn)接觸設(shè)置.而本文中，為了后續(xù)分析過盈量對(duì)結(jié)構(gòu)靜力學(xué)特性的影響，在接觸設(shè)置中采用接觸偏移量法，即通過給定偏移量大小來設(shè)置過盈量，這是本文與文獻(xiàn)[8]存在誤差的主要來源.綜上所述，可以說明本文所建模型的有效性.

表1 工況表

2.2 接觸特性對(duì)比

對(duì)于含有10個(gè)螺栓，預(yù)緊力為2 kN的含止口和無止口模型施加20 kN橫向載荷.當(dāng)橫向載荷達(dá)到7 kN時(shí)，對(duì)于無止口模型，兩個(gè)法蘭接觸面發(fā)生了分離，即由于無止口結(jié)構(gòu)無法承擔(dān)過大橫向載荷而導(dǎo)致求解無法收斂，發(fā)生分離的單元如圖4所示.

圖4 法蘭接觸界面分離圖

為了更直觀地展示含止口結(jié)構(gòu)和無止口結(jié)構(gòu)在施加橫向載荷時(shí)接觸狀態(tài)的變化趨勢，圖5給出了施加7 kN橫向載荷時(shí)無止口和含止口結(jié)構(gòu)接觸狀態(tài)的變化云圖.由圖可知，當(dāng)橫向載荷從1 kN 增加到3 kN時(shí)，含止口和無止口結(jié)構(gòu)在法蘭接觸界面的接觸狀態(tài)云圖幾乎無變化；當(dāng)橫向載荷增加到5 kN時(shí)，可以清晰地看出無止口結(jié)構(gòu)法蘭接觸界面出現(xiàn)均勻滑移現(xiàn)象，而含止口結(jié)構(gòu)仍存在粘合區(qū)域；當(dāng)橫向載荷為7 kN時(shí)，無止口結(jié)構(gòu)的接觸界面已經(jīng)發(fā)生了不均勻滑移現(xiàn)象，而含止口結(jié)構(gòu)依然存在粘合區(qū)域.由此可以充分說明，止口結(jié)構(gòu)的存在使結(jié)構(gòu)能夠承受更大的橫向載荷.

圖5 接觸狀態(tài)變化云圖

對(duì)于施加橫向載荷后被壓區(qū)域的接觸壓力變化趨勢也是關(guān)注的重點(diǎn).由圖6可知，當(dāng)載荷從1 kN增加至7 kN時(shí)，含止口和無止口結(jié)構(gòu)在法蘭接觸面的接觸壓力云圖并無較為明顯的變化.為了更清晰地表達(dá)出接觸壓力的變化趨勢，現(xiàn)提取法蘭界面接觸壓力最大值，并繪制出隨加載時(shí)間變化的曲線.

圖6 接觸壓力變化云圖

圖7為止口對(duì)接觸壓力的影響.圖中的0～1 s為螺栓預(yù)緊力施加階段，不加以討論.在施加橫向載荷的1～1.15 s階段，即橫向載荷為3 kN，接觸壓力的最大值主要體現(xiàn)在螺栓對(duì)法蘭的壓緊作用；在1.15～1.30 s階段，即橫向載荷為3～6 kN，接觸壓力的最大值主要體現(xiàn)在橫向載荷對(duì)法蘭接觸面的影響，即隨著橫向載荷的增加，接觸壓力逐步增大.在整個(gè)加載階段過程中，含止口結(jié)構(gòu)的接觸壓力最大值均低于無止口結(jié)構(gòu).由此可以說明，止口結(jié)構(gòu)的存在，可以有效降低接觸面的最大接觸壓力.

圖7 止口對(duì)接觸壓力的影響

2.3 彎曲剛度特性對(duì)比

為了更直觀地對(duì)比含止口和無止口模型在彎曲剛度上的特性，同時(shí)對(duì)二者施加6 kN的橫向載荷，將求解后的剛度值繪制成曲線，如圖8所示.由圖可知，隨著橫向載荷的逐步增大，二者的彎曲剛度均成下降趨勢，由于接觸問題具有高度非線性，其彎曲剛度也呈現(xiàn)出明顯的非線性，且無止口結(jié)構(gòu)的非線性更為明顯.在6 kN橫向載荷時(shí)，對(duì)于含止口模型，其剛度損失值為0.736×107N/m，剛度損失比為12.98%；對(duì)于無止口模型，其剛度損失值為1.425×107N/m，剛度損失比為32.38%，由此可以說明，止口結(jié)構(gòu)對(duì)于提高結(jié)構(gòu)彎曲剛度具有顯著作用.

圖8 止口對(duì)彎曲剛度的影響

3 靜力學(xué)特性參數(shù)影響分析

基于前文所建含止口結(jié)構(gòu)有限元模型的基礎(chǔ)上，為了討論結(jié)構(gòu)在接觸特性和彎曲剛度特性的影響規(guī)律，本節(jié)分別從止口長度、止口過盈量、螺栓個(gè)數(shù)、螺栓預(yù)緊力、摩擦系數(shù)以及不均勻接觸幾個(gè)角度進(jìn)行分析.

3.1 止口參數(shù)對(duì)接觸特性的影響

為討論止口參數(shù)對(duì)接觸特性的影響，分別從止口長度和止口過盈量這兩個(gè)角度，分析法蘭盤和止口配合接觸界面的接觸壓力云圖和接觸壓力最大值的變化趨勢.

對(duì)于不同止口過盈量，法蘭接觸面接觸壓力變化云圖如圖9所示，該圖是在2 kN螺栓預(yù)緊力下，止口長度為3 mm時(shí)的接觸壓力變化云圖.由圖可知，隨著止口過盈量的不斷增加，在法蘭接觸面上的接觸壓力云圖幾乎無變化，可見增加過盈量對(duì)法蘭接觸界面接觸壓力的分布幾乎無影響.對(duì)于不同止口長度法蘭接觸面的接觸壓力云圖如圖10所示，該圖是在2 kN螺栓預(yù)緊力下，過盈量為20 μm時(shí)的接觸壓力變化云圖.由圖可知，隨著止口長度增加，法蘭界面的最大接觸壓力分布區(qū)域逐漸減小.

圖9 不同止口過盈量下法蘭界面接觸壓力云圖

圖10 不同止口長度下法蘭界面接觸壓力云圖

為了更直觀地表達(dá)法蘭接觸面接觸壓力最大值的變化情況，現(xiàn)提取圖9和圖10的最大接觸壓力值，繪制成如圖11所示的曲線.由圖可知，隨著止口長度的增加，接觸界面的壓力最大值由199.133 MPa逐漸減小為185.250 MPa；隨著止口過盈量的增加，其法蘭接觸界面壓力的最大值穩(wěn)定在188 MPa左右.

圖11 法蘭界面最大接觸壓力

對(duì)于不同止口過盈量，止口配合接觸面接觸壓力變化云圖如圖12所示，該圖是在2 kN螺栓預(yù)緊力下，止口長度為3 mm時(shí)的接觸壓力變化云圖.由圖可知，隨著止口過盈量的不斷增加，止口配合接觸面上的接觸壓力云圖幾乎無變化.對(duì)于不同止口長度，止口配合接觸面的接觸壓力云圖如圖13所示，該圖是在2 kN螺栓預(yù)緊力下，過盈量為20 μm時(shí)的接觸壓力變化云圖.由圖可知，隨著止口長度的增加，由于其配合接觸面面積增大，使得接觸壓力分布區(qū)域逐步減小.

圖12 不同止口過盈量下止口界面接觸壓力云圖

圖13 不同止口長度下止口界面接觸壓力云圖

為了更直觀地表達(dá)止口配合接觸面接觸壓力最大值的變化情況，現(xiàn)提取圖12和圖13的接觸壓力最大值，繪制如圖14所示的曲線.由圖可知，隨著止口長度的增加，接觸壓力最大值從531.586 MPa逐步減小為282.073 MPa；隨著止口過盈量的增加，其接觸壓力的最大值穩(wěn)定在392 MPa左右.

圖14 止口界面最大接觸壓力

綜上所述，無論是法蘭接觸面還是止口配合接觸面，止口過盈量的變化對(duì)接觸區(qū)域壓力云圖和接觸壓力最大值幾乎無影響；止口長度的增加有利于減小受壓區(qū)域的面積，同時(shí)也會(huì)降低接觸壓力的最大值.

3.2 關(guān)鍵參數(shù)對(duì)彎曲剛度的影響

止口結(jié)構(gòu)的存在一定程度上提高了結(jié)構(gòu)的彎曲剛度，因此對(duì)彎曲剛度影響因素的討論尤為重要.本節(jié)分別討論螺栓個(gè)數(shù)、螺栓預(yù)緊力、止口過盈量、止口長度、摩擦系數(shù)以及不均勻接觸狀態(tài)對(duì)彎曲剛度的影響.

圖15展示了含止口螺栓連接結(jié)構(gòu)在不同工況下彎曲剛度的變化趨勢.如圖15a和圖15b所示，分別展示了10個(gè)螺栓連接時(shí)止口長度為3 mm和止口過盈量為20 μm時(shí)的彎曲剛度變化三維圖.由圖可知，彎曲剛度隨著螺栓預(yù)緊力的增加而增加，在相同的螺栓預(yù)緊力下，彎曲剛度隨著止口長度和止口過盈量而增加.如圖15c和圖15d所示，分別展示了2 kN預(yù)緊力下止口長度為3 mm和過盈量為20 μm時(shí)的彎曲剛度變化三維圖.由圖可知，彎曲剛度隨著螺栓個(gè)數(shù)的增加而增加，在螺栓個(gè)數(shù)相同的情況下，彎曲剛度隨著止口長度和止口過盈量的增加而增加.如圖15e所示，展示了10個(gè)螺栓，預(yù)緊力為2 kN，止口長度為3 mm，止口過盈量為20 μm下的彎曲剛度變化圖.由圖可知，彎曲剛度隨著螺栓預(yù)緊力和摩擦系數(shù)的增加而增加.

圖15 不同工況下彎曲剛度

為了更直觀地體現(xiàn)上述5種因素的增幅情況，現(xiàn)繪制5種工況下不同因素的增幅比條狀圖，如圖16所示.增幅比函數(shù)見式(6)，其中Kmax為某一工況下的最大彎曲剛度，Kmin為某一工況下的最小彎曲剛度，η為增幅比.

圖16 5種因素增幅比

(6)

對(duì)于因素1工況止口過盈量和螺栓預(yù)緊力的增幅比分別為0.418%和5.963%，對(duì)于因素2工況止口長度和螺栓預(yù)緊力的增幅比分別為13.567%和6.058%，因此可以說明止口長度對(duì)彎曲剛度的影響大于螺栓預(yù)緊力的影響大于止口過盈量的影響；對(duì)于因素3工況止口過盈量和螺栓個(gè)數(shù)的增幅比分別為0.313%和4.207%，對(duì)于因素4工況止口長度和螺栓個(gè)數(shù)的增幅比分別為4.146%和14.010%，由此可以充分說明螺栓個(gè)數(shù)對(duì)彎曲剛度的影響大于止口長度的影響大于止口過盈量的影響；對(duì)于因素5工況螺栓預(yù)緊力和摩擦系數(shù)的增幅比分別為3.825%和10.28%，由此可以說明摩擦系數(shù)對(duì)彎曲剛度的影響要大于螺栓預(yù)緊力對(duì)彎曲剛度的影響.綜上所述，螺栓個(gè)數(shù)對(duì)彎曲剛度的影響占主導(dǎo)地位，其次是止口長度、摩擦系數(shù)和螺栓預(yù)緊力的影響，止口過盈量對(duì)彎曲剛度的影響最小.

3.3 不均勻接觸對(duì)彎曲剛度的影響

在實(shí)際法蘭接觸界面中，由于加工誤差、表面粗糙度以及安裝誤差等原因，使得法蘭接觸界面出現(xiàn)接觸不均勻的情況.為了更加貼合實(shí)際工況，本節(jié)通過對(duì)10個(gè)扇區(qū)施加不同螺栓預(yù)緊力的方式進(jìn)行模擬，具體示意圖如圖17所示，值得說明的是，4種工況的預(yù)緊力總和為20 kN.基于上述4種工況，對(duì)不同止口過盈量和止口長度進(jìn)行仿真計(jì)算，得到關(guān)于不均勻接觸的彎曲剛度曲線圖.

圖17 不均勻接觸分布圖

由圖18a可知，彎曲剛度均隨著止口過盈量的增加而增加，工況1的彎曲剛度略高于工況2，某種程度上可以認(rèn)為二者的彎曲剛度基本相同；對(duì)于工況3來說，接觸薄弱的地方為法蘭的上下兩端，相比于工況1和工況2來說，其剛度會(huì)有小幅下降；對(duì)于工況4來說，接觸薄弱的地方為法蘭的左右兩側(cè)，相比于其他工況，其彎曲剛度大幅降低.由圖18b可知，彎曲剛度隨著止口長度的增加而增加，并且隨著止口長度的增加，彎曲剛度增加越緩慢，從工況1到工況4，彎曲剛度隨著不均勻程度的加深而降低.由此，可以說明，不均勻接觸會(huì)降低結(jié)構(gòu)的彎曲剛度，同時(shí)不均勻接觸的位置也會(huì)影響結(jié)構(gòu)的彎曲剛度.

圖18 不均勻接觸對(duì)彎曲剛度的影響

4 結(jié) 論

1)相比于無止口的螺栓連接結(jié)構(gòu)，含止口配合的螺栓連接結(jié)構(gòu)能夠有效降低接觸受壓區(qū)域面積和接觸壓力最大值，同時(shí)含止口結(jié)構(gòu)能夠承受更大的橫向載荷.對(duì)于無止口的螺栓連接結(jié)構(gòu)，隨著橫向載荷的逐步增加，其彎曲剛度損失越明顯，當(dāng)橫向載荷達(dá)到一定值時(shí)，法蘭接觸界面發(fā)生分離.

2)對(duì)于彎曲剛度的影響因素，起到主導(dǎo)作用的是螺栓數(shù)目，隨后是止口長度的影響和摩擦系數(shù)的影響，然后是螺栓預(yù)緊力的影響，而止口過盈量的影響最小，良好的接觸狀態(tài)對(duì)于提高彎曲剛度有一定的作用.