橫撐結(jié)構(gòu)壓扭極限強度的簡化計算方法研究

李鈞暉，趙 南，顧學(xué)康，吳劍國

(1.浙江工業(yè)大學(xué)，杭州310023；2.中國船舶科學(xué)研究中心，江蘇 無錫 214082）

0 引 言

撐桿是連接超大型浮體結(jié)構(gòu)下浮體的細(xì)長結(jié)構(gòu)，是平臺結(jié)構(gòu)的關(guān)鍵構(gòu)件，如圖1所示。在風(fēng)、浪、流的作用下，撐桿存在多種失效模式，其在組合載荷作用下的極限承載能力是業(yè)界關(guān)心的熱點問題。目前對于組合荷載的研究主要集中在船體梁上，Tatsumi[1]研究了船體梁在彎矩和局部荷載作用下的極限強度；吳嘉蒙[2]研究了集裝箱船在縱向彎曲以及扭轉(zhuǎn)下的極限強度；傅宇[3]采用有限元方法研究了船體梁在彎剪扭及其組合作用下的極限強度。目前國內(nèi)外關(guān)于超大型浮體結(jié)構(gòu)在復(fù)雜載荷作用下的結(jié)構(gòu)極限強度研究依然較少，顧學(xué)康[4]通過數(shù)值模擬和試驗相結(jié)合的方法，研究了超大浮體中連接器結(jié)構(gòu)在復(fù)雜荷載下的極限強度；顧海英[5]利用非線性有限元方法，計算了單一荷載下橫撐的極限強度，沒有考慮荷載的組合情況；張劍波[6]、姜峰[7]、楊鵬[8]等先后開展過半潛平臺結(jié)構(gòu)極限強度研究，考慮了整體結(jié)構(gòu)在組合荷載下的極限強度，但沒有考慮壓桿這一構(gòu)件本身可能發(fā)生的局部失效；趙南[9-10]通過橫撐結(jié)構(gòu)試驗和非線性有限元方法，給出了撐桿結(jié)構(gòu)在壓扭組合作用下的失效模式和極限強度的近似擬合曲線，并通過求解橫撐截面的翹曲函數(shù)的方法，計算壓扭下的極限強度，由于公式較為復(fù)雜，對設(shè)計計算人員要求較高。鑒此，有必要給出一套高效、簡便的組合荷載作用下橫撐結(jié)構(gòu)極限強度的計算方法。

本文以超大浮體中橫撐結(jié)構(gòu)為研究對象，針對軸向壓力和扭矩組合作用下的荷載工況，基于有限元計算結(jié)果和一定的假設(shè)，提出一套高效計算橫撐極限強度計算方法，并與有限元和試驗結(jié)果相對比，驗證計算方法的可行性，給出壓扭組合下橫撐結(jié)構(gòu)極限承載力曲線。

1 撐桿結(jié)構(gòu)的有限元應(yīng)力分析

針對一截面半徑為1.5 m，一跨長度為2.5 m，外圍曲板厚度為20 mm，中心十字板厚度為24.6 mm，骨材尺寸為500 mm×200 mm×15 mm×20 mm 的撐桿[9]進行有限元應(yīng)力分析。撐桿結(jié)構(gòu)主要包括縱骨、弱橫隔板、強橫隔板、相應(yīng)的縱隔板以及圓柱外殼和兩段的喇叭口等結(jié)構(gòu)，相應(yīng)的有限元模型以及去除外殼和縱隔板的內(nèi)部結(jié)構(gòu)見圖1（b），模型板格間劃分8個網(wǎng)格，縱骨劃分6個網(wǎng)格，面板劃分2個網(wǎng)格?？傮w的網(wǎng)格數(shù)量為10萬左右，其中四邊形單元采用S4R5，三角形采用STRI3單元模擬。

圖1 撐桿結(jié)構(gòu)示意圖Fig.1 Schematic diagram of brace strut

對于組合荷載下的極限強度計算，目前主流的計算思路為：先給定一個荷載，如在壓扭組合的工況下，給定一個軸向壓力F，求此壓力作用下，剩余扭矩的極限強度T，通過改變初始壓力的大小，得到不同壓力作用下扭矩的極限強度，獲得壓扭組合荷載下極限強度的包絡(luò)線。

有限元中邊界條件為一端固定，一端設(shè)置參考點施加組合荷載。Paik[11]研究指出組合載荷極限強度的計算結(jié)果與荷載施加順序和所取的路徑無關(guān)。因此，計算中分析步為兩步：第一步為靜力-通用計算方法，計算給定壓力F作用下的橫撐截面應(yīng)力分布；第二步為靜力-Risk法，基于第一步中的正應(yīng)力分布結(jié)果，迭代求得扭矩的極限強度T。圖2給出了壓縮、扭轉(zhuǎn)、壓扭組合等3種工況下極限狀態(tài)的應(yīng)力云圖。

圖2 三種工況下橫撐結(jié)構(gòu)極限狀態(tài)的應(yīng)力云圖Fig.2 Stress nephogram of the limit state of brace structure under three working conditions

從圖2（c）可以看出，在純扭下，外圍曲板的剪應(yīng)力較大，十字板格和加強筋的剪應(yīng)力非常小，可以忽略不計；從圖2（e）可以看出，壓扭組合作用下十字板格存在一定的剪應(yīng)力，但其剪流通過截面形心，不產(chǎn)生扭矩，而加強筋的剪應(yīng)力依然非常小，產(chǎn)生的扭矩可以忽略不計。因此，可以認(rèn)為不論哪種工況，外圍曲板承擔(dān)了所有的扭矩。從圖2（b）和（f）中可以看出，純壓時，整個截面均勻受壓，然而在壓扭作用下，中心板格和加強筋的正應(yīng)力比外圍曲板的正應(yīng)力大很多，外圍曲板承擔(dān)的壓力相對較小。

2 簡化計算的假設(shè)與流程

與船體梁彎曲極限強度相似，橫撐壓扭組合下的極限強度計算也是以橫向框架具有足夠的強度和剛度為前提，認(rèn)為橫向主要支撐構(gòu)件之間的板格和扶強材的屈曲或屈服的破壞導(dǎo)致橫撐整體垮塌。如Smith 法計算彎矩極限強度時，計算受壓區(qū)單元屈曲時，以一跨長度為單元計算長度。而且Paik[11]的有限元計算結(jié)果表明一跨模型和一艙模型的扭轉(zhuǎn)極限強度值基本相同，因此本文程序所計算的橫撐結(jié)構(gòu)均為一跨橫撐結(jié)構(gòu)模型?；趽螚U組合載荷作用下極限狀態(tài)的有限元結(jié)果，為了簡化極限狀態(tài)的計算過程，本文作出如下假設(shè)：

（1）抗扭能力的假設(shè)

橫撐的扭矩全部由外圍曲板板格承擔(dān)，忽略主要支撐構(gòu)件、扶強材的抗扭能力，僅將其作為曲板板格的支撐。

（2）軸向壓力分布的假設(shè)

軸力首先由中間十字板格與加強筋承受，當(dāng)軸力大于十字板格和加強筋的極限軸向受壓承載力時，外圍曲板板格才參與承擔(dān)軸向壓力。

（3）轉(zhuǎn)角相同的假設(shè)

假設(shè)撐桿橫截面上各點的轉(zhuǎn)角相同，橫截面上各點的剪應(yīng)變與該點到截面形心的距離成正比。

根據(jù)上述假設(shè)，給出撐桿結(jié)構(gòu)在組合載荷作用下的極限承載力計算流程如下：

（1）根據(jù)荷載類型的不同，將橫撐截面離散成一系列無相互作用的獨立單元。對軸向壓力作用，根據(jù)加強筋的位置，將外圍曲板離散成扶強材單元，中心十字為板單元，如圖2（a）所示。對扭矩作用，只考慮外圍曲板受扭，將兩個加強筋中間的曲板作為曲板板格單元，如圖2（b）所示。

（2）計算十字架的極限軸向受壓承載力：

式中，ti為第i個板格凈厚度，單位為mm；bi為第i個板格寬度，單位為mm；n為十字架所劃分的板格數(shù)目；σcxi為第i板格的極限受壓強度?？筛鶕?jù)《鋼質(zhì)海船入級規(guī)范》[12]中四邊簡支板格受壓屈曲的計算公式（為簡便起見，略去下標(biāo)i）：

式中，ReH-P為板格的屈服強度，Cx為屈曲折減因子，

式中，ψ為應(yīng)力比，取1；λ為板格的參考長細(xì)比，

λc為板格的臨界長細(xì)比，

K為屈曲因子，

σE-p為板格的參考應(yīng)力，單位為N/mm2，

式中，E為材料彈性模量，單位為N/mm2；其他符號含義見公式（1）。

（3）加強筋極限軸向受壓承載力為

式中，A1i為第i個加強筋的橫截面面積；σc1i為第i個加強筋的受壓極限強度，可以根據(jù)《鋼質(zhì)海船入級規(guī)范》中梁柱屈曲失效模式計算，具體計算公式如下（為簡便起見，略去下標(biāo)i）：

式中，σE為扶強材歐拉屈曲臨界應(yīng)力，計算公式如下：

式中，I為扶強材慣性矩，As為扶強材橫截面面積，L為扶強材跨長。

（4）根據(jù)施加力F的大小，判斷確定曲板的正應(yīng)力：

若F＜F1+F2，根據(jù)假設(shè)（2），認(rèn)為軸向壓力完全由十字板格和加強筋承擔(dān)，不考慮軸力對曲板扭轉(zhuǎn)極限強度的影響。

若F＞F1+F2，認(rèn)為軸向壓力由十字板格和加強筋、以及曲板共同承擔(dān)，根據(jù)假設(shè)（2），此時曲板所受壓應(yīng)力σax為

式中，A為曲板的總面積。此時應(yīng)考慮正應(yīng)力對曲板受剪極限強度折減，組合荷載下曲板極限強度計算公式為

式中，Cτ為剪切屈曲折減因子，

式中，Cax為曲板板格受壓屈曲折減因子，計算公式如下：

式中，λ為曲板的長細(xì)比，同式（4），而其中σE-P改為曲板板格的臨界應(yīng)力σE-bP，K改為曲板屈曲因子Kb，計算公式為

式中，R為曲板的半徑，tp為曲板的凈厚度，d為曲板的長度，單位均為mm。

（5）根據(jù)假設(shè)（3），通過步驟（4）得到外圍曲板的極限受剪承載力，依據(jù)文獻[14]中扭矩極限強度計算方法，通過構(gòu)造板格扭轉(zhuǎn)過程中剪應(yīng)力與剪應(yīng)變函數(shù)，迭代計算得到橫撐的受扭極限強度T。亦可通過文獻[15]的直接計算方法，計算得到橫撐的受扭極限強度T。文獻[13]中有關(guān)扭轉(zhuǎn)過程中剪應(yīng)力與剪應(yīng)變函數(shù)為

式中，γE為單元應(yīng)變，γy為單元達(dá)到屈服時的應(yīng)變，由下式確定：

式中，G為材料剪切彈性模量，β為板的柔度，計算公式為

式中：tp為板格厚度，單位為mm；b為板格短邊長度，單位為mm。

（6）計算給定軸壓力F作用下的極限扭矩T；改變施加壓力F的大小，即可獲得壓扭組合荷載下極限強度的包絡(luò)線。

計算流程如圖4所示。

圖4 壓扭組合作用下橫撐結(jié)構(gòu)極限強度的失效方程計算流程Fig.4 Flow chart of failure equation calculation for ultimate strength of transverse bracing structure under combined action of compression and torsion

3 結(jié)果對比

為驗證所提出計算方法的精度，應(yīng)用上述計算方法，編制了計算程序，計算純壓下的極限強度和計算文獻[9]的壓扭組合作用下橫撐極限強度失效方程，并與趙南論文[9-10]中試驗結(jié)果和實際橫撐結(jié)構(gòu)的有限元結(jié)果進行了對比。試驗結(jié)果對比如表1所示。實際橫撐結(jié)構(gòu)純扭以及壓扭組合工況下的計算結(jié)果如表2所示，壓扭組合作用下極限強度包絡(luò)線見圖5。不同軸壓力作用下?lián)螚U結(jié)構(gòu)載荷—位移曲線及承載能力曲線見圖6。

圖5 壓扭組合作用下實際橫撐結(jié)構(gòu)極限強度包絡(luò)線的對比Fig.5 Comparison of ultimate strength envelopes of actual brace structures under combined action of compression and torsion

圖6 不同軸壓力作用下?lián)螚U結(jié)構(gòu)載荷—位移曲線及承載能力曲線Fig.6 Load-displacement curve and interaction curve of brace strut under different compressions

表1 壓扭組合下試驗橫撐模型極限強度結(jié)果對比Tab.1 Comparison of ultimate strength of transverse brace model under compression and torsion combination

表2 不同壓扭組合下實際橫撐結(jié)構(gòu)極限強度的對比Tab.2 Comparison of ultimate strength of actual brace structures under different compression and torsion combinations

續(xù)表2

可見本文所提的計算方法精度較高，最大誤差在軸力等于1.076×108N 處，最大誤差約為8%，其余各點均小于此誤差。相比有限元方法和求翹曲函數(shù)的方法，本文方法簡便、高效，且精度較高。

為更好地應(yīng)用于實際工程中，對撐桿的外圍壁厚進行敏感性分析。采用有限元和本文所提方法分別改變外圍壁厚和十字板壁厚，進行了橫撐在軸向壓力8×107N 作用下極限扭矩敏感性分析，結(jié)果見表3。

表3 不同外圍壁厚及十字板格厚橫撐極限強度對比（單位：107N·m）Tab.3 Comparison of ultimate strength of transverse bracing with different thicknesses of outer wall and cross plate

從表3 可以看出，本文計算方法與有限元方法結(jié)果一致，精度較高；扭轉(zhuǎn)極限強度對外圍壁厚度更加敏感，驗證了假設(shè)（1）。

4 結(jié) 語

本文在超大浮體中橫撐結(jié)構(gòu)壓扭組合作用下有限元分析的基礎(chǔ)上，提出了撐桿結(jié)構(gòu)在組合載荷作用下極限狀態(tài)的假設(shè)，依據(jù)組合荷載下曲板極限狀態(tài)下失效方程提出了一套簡化的計算方法。從與有限元和試驗結(jié)果的對比可以看出，本文所提方法具有較高的準(zhǔn)確性，可以減少大量有限元建模和計算時間，可用于超大浮體橫撐的設(shè)計。