李世銘， 譚 越， 王春升

(中海油研究總院， 北京 100027)

?

基于ABAQUS的管節(jié)點應(yīng)力集中系數(shù)研究

李世銘， 譚 越， 王春升

(中海油研究總院， 北京 100027)

該文基于Abaqus/Standard程序，計算得出了有限元熱點應(yīng)力修正系數(shù)。在此基礎(chǔ)之上分析了K型管節(jié)點熱點應(yīng)力沿相貫線區(qū)域變化趨勢，討論了單元類型對于熱點應(yīng)力計算結(jié)果的影響以及有限元方法與SACS程序EFT經(jīng)典理論的區(qū)別，分析了空間管節(jié)點不同平面桿件之間對于應(yīng)力集中系數(shù)計算結(jié)果的影響。

ABAQUS；應(yīng)力集中系數(shù)；熱點應(yīng)力；疲勞設(shè)計

0 引言

在固定式導(dǎo)管架平臺設(shè)計過程中，疲勞壽命一般采用S-N曲線來描述，其中S代表了節(jié)點中沿焊縫周圍最大的應(yīng)力幅值，而應(yīng)力集中系數(shù)(SCF)直接決定了實際應(yīng)力循環(huán)幅值[1,2]。因此，將管節(jié)點應(yīng)力集中現(xiàn)象的分析結(jié)果作為海洋平臺管節(jié)點疲勞壽命估算的依據(jù)，以解決管節(jié)點的安全問題，是目前海洋平臺設(shè)計需要解決的重要課題之一[2]。

目前，就導(dǎo)管架平臺疲勞評估中應(yīng)力集中系數(shù)的計算，工程設(shè)計人員普遍采用SACS軟件中的Elthymiou公式[3]。該方法根據(jù)大量有限元計算結(jié)果擬合得出了適用于T型、Y型以及K型等簡單搭接關(guān)節(jié)點的方程組。然而，對于由多個簡單管節(jié)點焊接而成的復(fù)雜管節(jié)點的計算，SACS程序的做法是將其拆分為幾個簡單管節(jié)點分別計算SCF值之后再進行疲勞評估。因此，該文基于Abaqus/Standard程序，求解管節(jié)點應(yīng)力集中系數(shù)，以探討復(fù)雜管節(jié)點中不同平面的簡單管節(jié)點之間對于SCF值的相互影響，為管節(jié)點疲勞設(shè)計提供參考。

1 有效性驗證

1.1 驗證模型

計算方法有效性驗證是后續(xù)工作的前提，該文參考DNV-RP-C203疲勞設(shè)計規(guī)范[4]中的試件1和試件2進行分析，試件幾何模型如圖1所示。

圖1 試件幾何模型

1.2 熱點應(yīng)力計算及單元選取

表1 不同單元類型試件熱點應(yīng)力計算結(jié)果

由計算結(jié)果可知，目標(biāo)值與計算值的誤差基本控制在了15%以內(nèi)，具備了較高的工程精度。而且采用Abaqus/Standard計算得到的值都大于目標(biāo)值，后續(xù)的熱點應(yīng)力計算結(jié)果將乘以f以做修正。

2 K型管節(jié)點應(yīng)力集中系數(shù)數(shù)值分析

2.1 幾何模型

該文采用的K型管節(jié)點的幾何參數(shù)如圖2所示，具體尺寸見表2。

圖2 K型管節(jié)點的幾何參數(shù)符號

表2 K型管節(jié)點尺寸

2.2 載荷施加與邊界設(shè)置

采用API RP 2A中的加載方式以方便與SACS中EFT公式計算值進行對比，如圖3所示。

圖3 邊界及載荷條件

2.3 計算結(jié)果分析

分別采用殼元(CPS8)與固體元(C3D20)進行熱點應(yīng)力計算，計算結(jié)果考慮修正系數(shù)。該文僅以面內(nèi)彎矩加載時的情況來分析有限元計算結(jié)果與SACS程序EFT公式的結(jié)果進行對比，計算結(jié)果如圖4、圖5所示。圖5中Φ為沿焊縫一周的角度，零點位于Crown Heel冠踵。

圖4 面內(nèi)彎矩加載時K型管節(jié)點計算結(jié)果

圖5 面內(nèi)彎矩加載時弦管與撐管SCF值沿焊縫變化趨勢

由計算結(jié)果可以看出，對于弦管與撐管SCF值沿焊縫變化趨勢，殼元模型與三體元模型計算結(jié)果變化趨勢基本一致。表3列出了不同單元類型受平面內(nèi)彎矩作用時撐桿和弦桿的SCF值與EFT計算值的比較結(jié)果。

表3 受平面內(nèi)彎矩作用時撐桿和弦桿的SCF值比較

結(jié)合圖4、圖5和表3可以看出，與EFT經(jīng)典理論計算結(jié)果不同，由數(shù)值計算方法得出的熱點應(yīng)力并非位于冠部跟點，而是在與跟點夾角30°左右的位置。這與API RP 2A規(guī)范[4]描述一致，也在一定程度上證實了該文數(shù)值模擬的準(zhǔn)確性。另外，殼單元模型應(yīng)力集中系數(shù)結(jié)果較實體單元模型小，這是由于殼單元并不能模擬出管節(jié)點相貫線撐桿與弦管壁厚方向的應(yīng)力變化，即忽略了壁厚的影響?？梢妼嶓w單元能夠更真實地描述管節(jié)點相貫線附近的應(yīng)力集中現(xiàn)象。

然而，計算效率對于海工設(shè)計是一個不可忽視的因素。建立實體單元模型的最大缺點就是網(wǎng)格劃分慢，有時可能需要專門的網(wǎng)格劃分軟件,如Hypermesh、Femap等來提高網(wǎng)格質(zhì)量，此舉不可避免地降低了設(shè)計效率。所以面對設(shè)計效率與計算精度這一矛盾，采用8節(jié)點厚殼單元來提高效率也是可以的。

3 空間管節(jié)點應(yīng)力集中系數(shù)數(shù)值分析

SACS程序中EFT理論主要針對簡單管節(jié)點進行計算，即所有桿件都在一個平面或者平面間的角度不超過15°。例如，對于如圖6所示的空間管節(jié)點，SACS程序的EFT理論的具體做法是將其拆分成如圖7(a)和圖7(b)所示的K型節(jié)點和T型節(jié)點，分別計算其SCF值，再利用S-N曲線校核其壽命。為探究管節(jié)點不同平面桿件的影響，針對圖6、圖7所示模型進行述職計算，其中參照API RP 2A中K、T型簡單管節(jié)點進行加載以及邊界條件的設(shè)置。SACS程序以及有限元計算結(jié)果分別見表4、表5。為了方便比較，將圖6、圖7(a)與圖7(b)所示模型分別命名為模型1、模型2和模型3。

圖6 空間管節(jié)點模型(模型1)

圖7 空間管節(jié)點拆分結(jié)果

表4 SACS程序計算結(jié)果

節(jié)點類型桿件直徑(cm)厚度(cm)AX CRAX SDIN PLOU PLK型撐桿61．01．32．52．52．685．76弦桿76．22．23．493．492．247．74T型撐桿61．01．32．467．252．637．09弦桿76．22．24．088．482．729．51

表5 空間管節(jié)點應(yīng)力集中系數(shù)計算結(jié)果

其中：AX-CR表示軸向加載時冠點位置；AX-SD表示軸向加載時鞍點位置；IN-PL表示面內(nèi)彎矩加載；OU-PL表示面外彎矩加載。

由表4、表5計算結(jié)果可知，對于圖6所示的YZ平面的K型節(jié)點，XZ平面的T型撐桿將使其SCF值略增大；而反過來，YZ平面的K型節(jié)點的存在將大大減小XZ平面T型節(jié)點在各種加載模式下的SCF值。不同平面之間桿件的存在將影響其他平面桿件的應(yīng)力集中系數(shù)計算結(jié)果。而SACS程序中并沒有考慮平面外不同類型節(jié)點之間的相互影響，由此也將影響到疲勞壽命的評估。

4 結(jié)論

該文基于Abaqus/Standard模塊，參照最新DNV疲勞規(guī)范驗證了計算方法的有效性，并對簡單K型管節(jié)點與空間管節(jié)點進行應(yīng)力集中系數(shù)的計算，得到如下結(jié)論：

(1) 采用該文計算方法，DNV疲勞設(shè)計規(guī)范中的熱點應(yīng)力計算值與目標(biāo)值的誤差基本控制在了15%以內(nèi)，而且無論是殼單元還是實體單元，計算值都大于目標(biāo)值，即有限元熱點應(yīng)力計算皆小于1。

(2) 選用殼單元建模，計算時在減小儲存空間、縮短運算時間的基礎(chǔ)上保證了一定的精度。但其采用中面軸線建模的方式將無法真實反映關(guān)鍵內(nèi)外表面的應(yīng)力。實體單元能夠更真實的描述管節(jié)點相貫線附近的應(yīng)力集中現(xiàn)象，但是網(wǎng)格劃分的困難將導(dǎo)致設(shè)計效率降低。

(3) 相較于SACS程序中的EFT理論，采用有限元建模計算的方式能夠更加準(zhǔn)確的找到熱點應(yīng)力位置與大小。SCF經(jīng)典公式只針對簡單管節(jié)點而無法計算真實的、復(fù)雜的空間節(jié)點形式，而實際上，不同平面之間桿件的存在將影響其他平面桿件的應(yīng)力集中系數(shù)計算結(jié)果。因此，復(fù)雜空間節(jié)點的詳細(xì)數(shù)值建模很有必要。

[1] Peter W M. Design of Welded Tubular Connections-Basis and use of AWS Code Provisions[R]. Civil Engineering Consultant, Shell Oil Company, Houston, Texas, U.S.A. 1992.

[2] 楊錚,金偉良. Y型管節(jié)點應(yīng)力集中系數(shù)有限元分析[J]. 中國海洋平臺，2004，19(4)：17-21.

[3] 姜萌. 近海工程結(jié)構(gòu)物—導(dǎo)管架平臺[M]. 大連：大連理工大學(xué)出版社，2005.

[4] API RP 2A. Recommended Practice for Planning, Designing, and Constructing Fixed Offshore Platforms-Working Stress Design[S].2000.

[5] DNV-RP-C203.Fatigue strength analysis of offshore steel structures[S]. 2010.

Research on Stress Concentration FactorBased on ABAQUS

LI Shi-ming， TAN Yue， WANG Chun-sheng

(CNOOC Research Institute, Beijing 100027, China)

The paper obtained the hot spot stress correction factor by referring to DNV RP-C203. On These Foundations, variation tendency of k joints intersecting line area is analyzed; the influence of element type on calculated result of FE hot spot stress and the difference between FE method and Efthymiouequation are discussed; the influence of non-coplanar member on SCF is calculated.

ABAQUS; SCF; hot spot stress; fatigue design

2014-11-06

國家科技重大專項 “南海深水油氣開發(fā)示范工程”(2011ZX05056)。

李世銘(1987-)，男，工程師。

1001-4500(2015)03-0095-06

P75

A