淺海井組計量平臺設計含4張CAD圖

淺海井組計量平臺設計含4張CAD圖,淺海,計量,平臺,設計,cad 第一章 前言 隨著社會的發(fā)展和科學技術的進步，人類社會對能源的需求越來越大。陸地上的油氣資源經(jīng)過長時期大規(guī)模的開發(fā)之后已日趨枯竭，油氣勘探與開發(fā)漸漸轉向了資源豐富的海洋，并形成了投資高、風險大并且高新技術密集的海洋工程產(chǎn)業(yè)。海洋石油開發(fā)是海洋資源開發(fā)利用的一部分。目前世界上己有39個國家(或地區(qū))從事近海石油開發(fā)，22個國家(或地區(qū))從事近海天然氣開發(fā)。 我國海域遼闊，其中大陸架面積約有110萬平方公里，渤海，黃海，東海和南海都有大面積的沉積盆地，其中具有油氣勘探價值的面積在60萬平方公里以上，即一半以上的海域有寶貴的石油。預測的石油儲量達250億噸，這是我國海上石油天然氣開發(fā)的豐富資源基礎。我國從1957年便已開始進行海洋石油勘探開發(fā)，從1979年實行改革開放以來，我國的海上石油開發(fā)更是進入高速發(fā)展期，到2000年生產(chǎn)能力已達到2000萬噸：預計到2005年，我國的海上原油生產(chǎn)能力會達到4000萬噸。目前海洋石油己成為我國重要的原油生產(chǎn)基地。 在海洋開發(fā)中，尤其是對淺海的開發(fā)中，各國廣泛采用導管架平臺。導管架平臺具有如下特點：(1)平臺的支承結構是以圓鋼管為主要構件的鋼結構，因此結構受力狀態(tài)較好。(2)由于樁是通過導管架腿柱打入海床，因此在惡劣的海洋環(huán)境中，樁打的準、打的直。(3)打樁作業(yè)大大簡化，保證平臺結構的整體穩(wěn)定性。(4)平臺可以在陸上分塊預制，海上組裝，保證施工質量，節(jié)約投資。(5)導管架平臺設計、制造及安裝技術成熟，實踐經(jīng)驗多，適用性強。 隨著計算機技術的發(fā)展進步，SAP、MSC/NASTRAN、ANSYS等各種結構分析軟件被廣泛的應用到平臺及導管架的設計和強度計算中，并且發(fā)揮了重要的作用。其中，ANSYA作為世界著名的美國ANSYS軟件公司的最具盛名的CAE軟件，其在海洋工程領域的應用更為通用廣泛，它提供了支持圓管形構件的流體靜力、動力效應的Pipe59單元，能夠有效模擬海洋環(huán)境中的導管架結構，計算構件在海水中受浮力、波浪力、流力等載荷的影響。在淺海平臺結構的研究和設計中，ANSYS已經(jīng)得到成功的運用。 1.1 國外海洋工程發(fā)展狀況 世界石油開發(fā)的歷史已有200多年， 1990－1995年間，除美國以外全世界共安裝了703座平臺，其中83座為半潛式、張力腿式和可移動生產(chǎn)平臺，41各國家安裝了370多座水深不超過60m的淺水平臺。 長期以來，樁基導管架平臺是世界海洋石油生產(chǎn)中采用最廣泛的一種結構。以墨西哥灣海域的平臺發(fā)展為例，1978年建造的“Cognac”鉆井平臺以極端的颶風載荷控制設計，導管架分三段建造；1981年水深285m的“Cerveza”平臺則使用更為先進的整體制造，簡化了結構且成本大大降低。而北海的環(huán)境條件比墨西哥灣還要惡劣很多，海水腐蝕也更為嚴重，1974年建造的一座平臺，其導管架所有的樁基在四角上，除主樁腿外周圍又打入數(shù)根環(huán)繞樁，以增大抗傾覆能力。 1.2 我國海洋工程發(fā)展現(xiàn)狀 我國海域遼闊，大陸架面積約有110萬平方公里，管轄海域近300萬平方公里。為開發(fā)利用海洋，我們建設了大量的海洋工程，其中用于油氣開發(fā)的海洋平臺100多座。海洋環(huán)境十分復雜和惡劣，風、海浪、海流、海冰和潮汐時時作用于結構，同時還受到地震和海嘯的威脅。從1966年渤海建造第一座鋼質鉆井平臺到1980年對外合作勘探開發(fā)的15年里，我國自主設計建造了11座固定式鉆井平臺、7座固定式平臺、1座自升式平臺、1座單點系泊系統(tǒng)和3艘工程船舶。 雖然我國海洋石油開發(fā)較晚，但近年來通過對外合作，引進國外先進技術，加上自己研究開發(fā)，已在該領域取得了顯著的進步。2004年，國內(nèi)最大石油平臺——渤海南堡35－2油田開發(fā)項目平臺組塊建造工程開工，南堡35－2CEP/WHPB平臺總重達到了12000多噸。其中CEP平臺組塊長64米、寬59米、高20.6米，重達8000余噸，WHPB平臺也達到了4000噸。 本文介紹了淺海9米水深海上石油計量導管架平臺的設計，主要包括以下幾方面內(nèi)容： 1． 學習ANSYS在平臺結構中的應用； 2． 平臺選型、主尺度和構件尺寸的確定； 3． 環(huán)境載荷計算及各工況載荷的組合； 4． 建立ANSYS模型； 5． 對平臺進行靜動力強度校核； 6． 地基承載力計算 第二章 環(huán)境條件和設計依據(jù) 2.1 環(huán)境條件 設計水深： 9米 潮位以黃海平均海平面為基準 校核高水位（50年重現(xiàn)期）： 3.08米 設計高水位： 1.48米 設計低水位： -0.69米 校核低水位： -2.32米 波浪： 風浪為主，涌浪次之，強浪向NNE-ENE 校核高水位最大可能波高： 7.2米 對應波浪周期： 8.6秒 設計高水位最大可能波高： 6米 對應波浪周期： 8.6秒 海流： 方向：ESE，WNW 最大可能流速： 1.6m/s 最大實測流速： 0.98 m/s 風： 強風向：NW-NNW，NNE-ENE 風速： 冰作用下，設計風速28m/s 波浪作用下，設計風速51.5m/s 海冰： 設計冰厚（50年一遇）： 0.45m 抗壓強度： 2244kPa 腐蝕和磨損： 飛濺區(qū)構件腐蝕裕量： 3mm 飛濺區(qū)定義標高： -2.69m～+5.84m 冰接觸區(qū)構件磨損量： 1mm 冰接觸區(qū)標高： -1.14m～+1.93m 地質資料： 表2.1地基土壤的物理學性質 土壤 層號 土壤 名稱 深度 m 水下容重 KN/m3 剪切強度 C(Kpa) 摩擦角 ?（度） 1 非常軟的粘土 0－2 5.0 3.0 2 粉砂質粘土 2－4 7.0 8.0 4－5.15 8.5 40.0 5.15－7.55 8.0 18.0 3 細砂 7.55－11.35 8.5 0 30 4 硬的細砂質粘土 11.35－15.25 10.0 80.0 5 粉砂 15.25－18.3 10.5 25 6 粉砂質細砂 18.3－23.4 10.0 30 7 粉砂質粘土 23.4－27.8 10.0 50.0 8 粉砂 27.8－31.1 9.0 30 2.2 設計依據(jù) 平臺用途和主要功能： 本文所設計的導管架平臺為淺海井組計量平臺，包括甲板、導管架、樁基等部分。平臺上部設備總重量為200噸。 平臺用鋼材： 平臺所用鋼材為：樁基礎和導管架均用，甲板使用Q-235-A。按照規(guī)范要求，鋼材的屈服應力為315Mpa，許用應力取189 Mpa；Q-235-A鋼材的屈服應力為235 Mpa，許用應力取141 Mpa。 使用年限： 平臺的使用年限為15年。 依據(jù)規(guī)范： 中國船級社《淺海固定平臺規(guī)范》（2003） 第三章 平臺選型和主尺度 設計要求：平臺甲板面積17m×12m，承載能力200t。 考慮環(huán)境條件，決定采用摩擦樁基礎四腿導管架平臺： 平臺甲板高程：9.00 m 工作點高程：5.50 m 導管架尺寸：上部工作點處（EL+5.50m）為9m×8m 底部（EL-9.00m）為11.9m×10.9m 樁腿導管直徑為?1340×25，斜度為1/10。 導管架設3層水平橫撐，潮差帶不設斜撐，水平外圍橫撐與斜撐尺寸均為?610×20，水平內(nèi)圍橫撐尺寸為?400×20。 導管架底部設置防沉板，防沉板厚度為8 mm。 選取摩擦樁樁徑為1.2m，則根據(jù)CCS規(guī)范鋼管樁壁的最小厚度t按下式計算： t=6.35+D/100 mm 式中 D一樁徑，mm 則取樁壁厚為30 mm > t=6.35+1200/100=18.35 mm。 樁入土深度為30m，總長為46.3m。 平臺甲板采用板、梁結構，面積為17m×12m，甲板板厚為 8mm ，Y向設5根主梁，X向設4根，主梁為600×300×25工字梁。 表3.1 導管架結構重量計算 名稱 截面面積（㎡） 長度（m） 密度 數(shù)量 質量（kg） 導管?1340×25 0.103 15.58 7800 4 50067.888 EL+4.00水平外圍橫撐 0.037 7.23 7800 4 8346.312 EL-2.00水平外圍橫撐 0.037 8.43 7800 4 9731.592 EL-9.00水平外圍橫撐 0.037 9.83 7800 4 11347.752 EL+4.00水平內(nèi)圍橫撐 0.024 5.43 7800 4 4065.984 EL-2.00水平內(nèi)圍橫撐 0.024 6.27 7800 4 4694.976 EL-9.00水平內(nèi)圍橫撐 0.024 7.27 7800 4 5443.776 斜撐?610×20 0.037 7.34 7800 8 16946.592 合 110644.872 第四章 環(huán)境條件計算 本章分別計算風浪、風冰兩種工況下X和Y向的載荷。 4.1 風載荷 4.1.1 計算公式 作用于平臺上的風載荷按下式計算： 式4.1 式中 p—風壓； A—結構垂直于風向的投影面積。 載荷作用在上述投影面積的形心位置。 其中，結構所承受的風壓為： 式4.2 為基本風壓，標準高度為海面上10m 式4.3 式中 —風壓的高度系數(shù)； —構件的形狀系數(shù)。 、取值見表4.1和表4.2 表4.1 高度系數(shù) 表4.2 形狀系數(shù) h(m) CH 構件形狀 Cs 0~15.3 1.00 球 0.4 15.3~30.5 1.10 圓柱 0.5 30.5~46.0 1.20 大平板 1.0 46.0~61.0 1.30 鉆井架 1.25 61.0~76.0 1.37 甲板以下暴露的梁和桁材 1.3 76.0~91.5 1.43 孤立結構 1.5 其中h（m）為構件距離海平面的高度 4.1.2 計算結果 表4.3 風、浪、流工況下X向風載計算結果：(風速51.5m/s) 名稱 底部位置 b或r h Ch Cs P0 P（Pa） F(kN) 作用點 M(kN*M) 工作間 9.0 13 5 1 1 1625.83 1625.83 105.679 11.50 1215.308 主梁 8.2 17 0.8 1 1 1625.83 1625.83 22.111 8.60 190.157 連接構件 5.5 1.2 3.5 1 1 1625.83 812.915 3.414 7.25 24.753 樁腿 5.0 1.2 0.5 1 1 1625.83 812.915 0.487 5.25 2.560 導管架 0 1.34 5 1 1 1625.83 812.915 5.473 2.51 13.752 導管架橫撐 4.0 0.61 8.2 1 1 1625.83 812.915 4.082 4 16.328 合 141.248 合 1462.860 合力作用點高度=1462.860/141.248=10.35m 等效載荷： F=141.248kN M=141.248X(10.13-9)=160.175kN*M 作用點距靜水面9米 表4.4 風、浪、流工況下Y向風載計算結果：（風速51.5m/s） 名稱 底部位置 b或r h Ch Cs P0 P（Pa） F(kN) 作用點 M(kN*M) 工作間 9.0 9 5 1 1 1625.83 1625.83 73.162 11.5 841.367 主梁 8.2 12 0.8 1 1 1625.83 1625.83 15.607 8.6 134.228 連接構件 5.5 1.2 3.5 1 1 1625.83 812.915 3.414 7.25 24.753 樁腿 5.0 1.2 0.5 1 1 1625.83 812.915 0.487 5.25 2.560 導管架 0 1.34 5.025 1 1 1625.83 812.915 5.473 2.51 13.752 導管架橫撐 4.0 0.61 8.232 1 1 1625.83 812.915 4.082 4 16.328 合 102.228 合 1032.991 合力作用點高度=1032.991/102.228=10.10m 等效載荷 F=102.228kN M=102.228X(10.10-9)=112.450kN*m 作用點距靜水面9米 表4.5 風、冰、流工況下X向風載計算結果：（風速28m/s） 名稱 底部位置 b或r h Ch Cs P0 P（Pa） F(kN) 作用點 M(kN*M) 工作間 9 13 5 1 1 480.6 480.59 31.238 11.5 359.241 主梁 8.2 17 0.8 1 1 480.6 480.59 6.536 8.6 56.209 .連接構件 5.5 1.2 3.5 1 1 480.6 240.295 1.009 7.25 7.316 樁腿 5 1.2 0.5 1 1 480.6 240.295 0.144 5.25 0.756 導管架 0 1.34 5.025 1 1 480.6 240.295 1.618 2.5125 4.065 導管架橫撐 4 0.61 7.228 1 1 480.6 240.295 1.059 4 4.237 合 41.605 合 431.828 合力作用點高度=431.828/41.605=10.38m 等效載荷 F=41.605kN M=41605.3X(10.38-9)=57.415kN*m 作用點距靜水面9米 表4.6 風、冰、流工況下Y向風載計算結果：（風速28m/s） 名稱 底部位置 b或r h Ch Cs P0 P（Pa） F(kN) 作用點 M(kN*M) 工作間 9 9 5 1 1 480.59 480.59 21.626 11.5 248.705 主梁 8.2 12 0.8 1 1 480.59 480.59 4.613 8.6 39.677 連接構件 5.5 1.2 3.5 1 1 480.59 240.295 1.009 7.25 7.316 樁腿 5 1.2 0.5 1 1 480.59 240.295 0.144 5.25 0.756 導管架 0 1.34 5.025 1 1 480.59 240.295 1.618 2.5125 4.065 導管架橫撐 4 0.61 7.228 1 1 480.59 240.295 1.059 4 4.237 合 30.071 合 304.760 合力作用點高度=304.760/30.071=10.134m 等效載荷 F=30.071kN M=30071.14X(10.134-9)=34.100kN*M 作用點距靜水面9米 注：為了在進行強度校核時簡化計算并且減少應力集中，根據(jù)力的平移法則，將風載等效為一個集中力加一個彎矩，作用于導管架頂部，作用點距靜水面9米。 表格中“底部位置”均指構件底部距海平面的距離，“合力作用點高度”亦是指合力作用點距海平面的高度。 4.2 冰載荷 4.2.1 計算公式 作用于平臺上的冰載荷按下式計算： 式4.4 式中：—樁柱形狀系數(shù)，對圓截面柱采用0.9； —局部擠壓系數(shù)； —樁柱與冰層的接觸系數(shù)； —樁柱寬度（或直徑）； —冰層計算厚度。 4.2.2 計算結果 根據(jù)CCS規(guī)范擠壓系數(shù)K1=2.5；接觸系數(shù)K2=0.45； 由環(huán)境條件：＝2244kPa； 由圖紙資料：b=1.34m h=0.45m。 則單個樁腿所受的冰載荷為： 0.9×2.5×0.45×1.34×0.45×2244=1370kN 考慮群樁產(chǎn)生的遮蔽效應和堵塞作用，受遮蔽的樁腿受到的冰載荷為： 0.3＝1370×0.3＝411 kN 根據(jù)CCS規(guī)范，當樁腿之間距離L>8D時不考慮“群樁效應”，D為樁的直徑。 有 8D=8×1.2＝9.6m 依據(jù)圖紙，側向樁泥線處最小間距為11.2m>9.6m，所以不考慮“群樁效應”。 4.3 波浪、流載荷 4.3.1 計算公式 對小尺度圓形構件，垂直于其軸線方向單位長度上的波浪力， 當D/L≦0.2（D為圓形構件直徑，m；L為設計波長，m）時，可按Morison公式計算： N/m 式4.5 式中： ——海水密度，kg/m3 ； ——垂直于構件軸線的阻力系數(shù)。必要時，應盡量由試驗確定。在實驗資料不足時，對圓形構件可取=0.6~1.0； ——慣性力系數(shù)，應盡量由試驗確定，在實驗資料不足時，對圓形構件可取2.0； ——水質點相對于構件的垂直于構件軸線的速度分量，m/s， 為其絕對值，當海流和波浪聯(lián)合對平臺作用時，為水質點的波浪速度矢量與海流速度矢量之和在垂直于構件方向上的分矢量； ——水質點相對于構件的垂直于構件軸線的加速度分量，m/s2。 當只考慮海流作用時，圓形構件單位長度上的海流載荷按下式計算： N/m 式4.6 式中： ——阻力系數(shù)； ——海水密度，kg/m3； ——設計海流速度，m/s； ——單位長度構件垂直于海流方向的投影面積，m2/m。 設計海流速度采用平臺使用期間可能出現(xiàn)的最大流速。 4.3.2 計算結果 在ANSYS程序中，提供了支持圓管形構件的流體靜力、動力效應的Pipe59單元，能夠有效模擬海洋環(huán)境中的導管架結構，所以需要在ANSYS模型建立之后，將有關波浪和海流參數(shù)填入water table表格中，程序將根據(jù)所選用的波浪理論對使用了Pipe59單元的結構進行波浪力及流力的計算。 考慮到所給出的環(huán)境資料（H/d>0.2），選用斯托克斯五階波進行計算，斯托克斯五階波相關公式如下： 波面方程為： 式4.7 迭代求L、方程為： 式4.8 其中為深水波長 式4.9 速度勢方程為： 式4.10 相位角為： 式4.11 波形系數(shù)： 速度勢函數(shù)： 其中 、、、為系數(shù)。 4.4 甲板設備載荷 該平臺為石油計量平臺，其設備總重約為200噸，計算時甲板載荷等效為4個各50噸的質量單元，作用于樁腿頂端的節(jié)點之上；平臺結構的自重通過輸入z軸方向的重力加速度（9.8m/）由ANSYS程序自動生成。 4.5 工況和載荷組合 4.5.1 波流工況 1）X方向上（0°方向），波浪力+海流力+風力+結構自重+甲板設備重 2）Y方向上（90°方向），波浪力+海流力+風力+結構自重+甲板設備重 4.5.2 海冰工況 1）X方向上（0°方向），冰力+海流力+風力+結構自重+甲板設備重 2）Y方向上（90°方向），冰力+海流力+風力+結構自重+甲板設備重 第五章 建立數(shù)學模型及ANSYS程序和單元的相關說明 5.1 有限元法基本思路 有限元法是一種利用電子計算機求解結構靜、動態(tài)力學特征問題的數(shù)值解法，其基本思路是： 1.把很復雜的結構拆分為若干個形狀簡單的單元，這些單元一般要小到可以用簡單的數(shù)學模型來描述特征參數(shù)在其中的分布，這一過程即離散。 2.通過對單元的分析來建立各特征參數(shù)之間的關系方程，即單元分析。彈性力學中，單元分析的任務是建立聯(lián)系應變和節(jié)點位移分量的方程，同時研究單元的節(jié)點力和節(jié)點位移之間的關系，以及把作用在單元中間的外載荷轉化為節(jié)點載荷。 3.在單元分析的基礎上，利用平衡條件和連續(xù)條件將各個單元拼裝成整體結構。對整體在確定邊界條件下進行分析，從而得到整體的參數(shù)關系方程，即整體矩陣方程。這一過程為整體分析。 4.求解整體矩陣方程，即可得到各種參數(shù)在整體結構中的分布。 5.2 ANSYS有限元分析軟件概況： 1970年Dr .John A.Swanson成立了Swanson Analysis System，Inc.(SASI)，后經(jīng)重組改稱ANSYS眾司。ANSYS是世界CAE行業(yè)最著名的公司之一，其總部位于美國賓西法尼亞州的匹茲堡。三十年來，ANSYS公司一直致力于設計分析軟件的開發(fā)、維護以及售后服務等，并不斷吸取最新的計算方法和計算機技術，始終領導著國際分析仿真和優(yōu)化技術的新潮流，為全球工業(yè)界所認同，擁有十分廣泛的用戶群。 作為新一代有限元分析軟件的代表，ANSYS抓住了現(xiàn)代設計方法對CAE技術提出的新要求，即以廣泛的多物理場仿真分析的功能，承擔起虛擬樣機這一設計核心的角色。ANSYS軟件可以實現(xiàn)多場及多場藕合分析，具有強大的非線性分析功能。它是一種可以實現(xiàn)前后處理、求解及多場分析統(tǒng)一數(shù)據(jù)庫的一體化大型FEA軟件，并可支持從微機、工作站到巨型機的所有平臺以及所有平臺的并行計算。 5.3 單元特性 本文采用了以下ANSYS單元模型進行計算： 5.3.1 Beam188 單元特性 Beam188 單元適合于分析從細長到中等粗短的梁結構，該單元基于鐵木辛哥梁結構理論，并考慮了剪切變形的影響。 Beam188 是三維線性（2 節(jié)點）或者二次梁單元。每個節(jié)點有六個或者七個自由度，自由度的個數(shù)取決于KEYOPT(1)的值。當KEYOPT(1)＝0（缺?。r，每個節(jié)點有六個自由度；節(jié)點坐標系的x、y、z 方向的平動和繞x、y、z 軸的轉動。當KEYOPT(1)=1 時，每個節(jié)點有七個自由度，這時引入了第七個自由度（橫截面的翹曲）。這個單元非常適合線性、大角度轉動和/并非線性大應變問題。 Beam188/beam189 可以采用sectype、secdata、secoffset、secwrite 及secread 定義橫截面。本單元支持彈性、蠕變及素性模型（不考慮橫截面子模型）。這種單元類型的截面可以是不同材料組成的組和截面。 5.3.2 Pipe16單元特性 圖5.1 Pipe16單元模型 Pipe16是一種單軸單元，具有拉壓、扭轉、和彎曲性能。 該單元在兩個結點有6個自由度：沿節(jié)點X，Y，Z方向的平移和繞結點X，Y，Z軸的旋轉。該單元基于三維梁單元(BEAM4)，包含了根據(jù)對稱性和標準管幾何尺寸進行的簡化。 圖5.1描述顯示了該單元的幾何形狀，節(jié)點位置和坐標系。單元的數(shù)據(jù)輸入包括2個或3個節(jié)點，管的外部直徑(OD) 和管壁厚度(TKWALL)，應力增量系數(shù)(SIF)與撓曲系數(shù)(FLEX)，內(nèi)部流體密度(DENSFL)，外部絕緣層密度(DENSIN)與厚度(TKIN)，允許侵蝕厚度(TKCORR)，絕緣表面積(AREAIN)， 管壁質量(MWALL)， 管的軸向剛度(STIFF)，基于轉子動力學的自旋頻率(SPIN)，和各向同性材料性質。 該單元的X軸為從I結點到J節(jié)點的方向。當單元由2個結點組成時，單元的Y軸被自動設成平行于整體坐標系的X-Y面。參見圖16.1 PIPE16幾何描述。在單元平行于Z軸的情況下（或在0.01%坡度范圍內(nèi)），該單元的Y軸是與整體坐標系的Y軸（如圖示）平行。用戶若想人為控制單元X軸的方向，需定義第3個節(jié)點。如果使用了第3節(jié)點(K)，則K和I、J一起定義了一個包括單元X軸和Y軸的平面 (如圖示)。繞管圓周的輸入與輸出位置定義為：若沿單元Y軸，為0度；類似的，沿Z軸，為90度。 應力增量系數(shù)(SIF)影響彎曲應力。若 KEYOPT(2) = 0，則應力增量系數(shù)在I(SIFI)節(jié)點和J(SIFJ)節(jié)點末端輸入；若KEYOPT(2) = 1，2，或3，則軟件按照T形接頭自行計算并確定。當SIF的值小于1.0時取1.0。撓曲系數(shù)(FLEX)被分為截面的轉動慣量中從而生成一個計算彎曲剛度所需的修正轉動慣量。FLEX的缺省直為1.0，但也可以輸入任何一個正數(shù)。 單元的質量可由管壁材料，外部絕緣體和內(nèi)部流體計算得到。絕緣體和流體決定了單元的質量矩陣。允許侵蝕厚度用來計算應力。一個確定的管壁其質量是一個常數(shù)，不用考慮其計算值。一個非零的絕緣體面積是一常數(shù)，不用考慮其計算值(有管壁外直徑和長度得來)。一個非零的剛度也是一常數(shù)，不考慮管的軸向剛度。 5.3.3 Pipe59單元特性 pipe59單元是一種可承受拉、壓、彎作用，并且能夠模擬海洋波浪和水流的單軸單元。單元的每個節(jié)點有六個自由度，即沿x，y，z方向的線位移及繞X，Y，Z軸的角位移。除了本單元的單元力包括水動力和浮力效應，單元質量包括附連水質量和內(nèi)部水質量，其余與單元pipe16相似。pipe59還可以模擬纜索單元，和link8相似。這個單元還適合剛度硬化和非線性大應變問題。 圖5.2 Pipe59 單元模型 pipe59輸入數(shù)據(jù)： 圖5.2給出了單元的幾何圖形、節(jié)點位置及坐標系統(tǒng)。本單元輸入數(shù)據(jù)包括：兩個節(jié)點，管外徑，壁厚以及一些荷載和慣性信息；各向同性材料屬性；外部附著物（包括冰荷載和生物附著物）；材料粘滯系數(shù)用來計算外部流體的雷諾系數(shù)。 單元的x 軸方向為i 節(jié)點指向j 節(jié)點，y 軸方向按平行x－y 平面自動計算，其他方向如上圖所示。對于單元平行與z 軸的情況（或者斜度在0.01％以內(nèi)）， 圖5.3 波、流作用方向 單元的y 軸的方向平行與整體坐標的y 軸（如圖5.3）。被認為為0度的單元的外部環(huán)境輸入或輸出沿y 軸分布，就和90度的單元的外部環(huán)境輸入或輸出沿z 軸分布類似。 KEYOPT(1)用來消除抗彎剛度將管單元轉換為纜索單元，如果構件存在扭轉，就可以用KEYOPT(1)=2來解決。KEYOPT(2)可以定義質量矩陣是團聚質量矩陣還是一致質量矩陣，可用于長柔結構，常用來分析帶扭轉的纜索結構。 對于海浪，海流和水密度通過water motion table輸入，如果不輸入，就不會考慮水的作用。雖然文章中用“水”的不同性質，事實上還可用于其他液體性質描述，不同的曳力系數(shù)和溫度數(shù)據(jù)也可以通過此表輸入。 此外，單元長度不能為0，此外，外徑必須大于0，內(nèi)徑不能小于0。水面附近的單元長度相對波浪長度應小。單元的兩個端點不能同時處于泥面以下，并且如果積分點位于泥面以下，就會忽略水動力。如果單元位于水平面上，可以不考慮PIPE59 Water Motion Table作用。在用縮減法進行瞬態(tài)動力分析時，應當注意。因為這種分析會忽略單元上的荷載向量。 5.3.4 MASS21單元特性 MASS21 是一個具有六個自由度的點元素： 即 x， y， 和 z 方向的移動和繞 x， y， 和 z 軸的轉動。每個方向可以具有不同的質量和轉動慣量。 圖5.4 Mass21 單元模型 此質點元素由一個單一的節(jié)點來定義， 此單元的坐標系統(tǒng)可以平行于全局的笛卡爾坐標系統(tǒng)也可以使用節(jié)點坐標系統(tǒng) (KEYOPT(2)). 在節(jié)點坐標系的操作中可以看到對節(jié)點坐標系元素的操作討論。做大偏轉（a large deflection analysis）分析時元素坐標系相對于節(jié)點坐標系會產(chǎn)生旋轉?？梢酝ㄟ^KEYOPT(3)選項來禁止轉動慣量的產(chǎn)生或使元素退化為二維形式。 如果一個元素只需一個輸入，那末這個輸入將運用到其它所有方向。此元素的坐標系顯示如上圖所示。 使用該單元的假設與限制：質量單元在靜態(tài)解中無任何效應，除非具有加速度或旋轉載荷或慣性解除，如果質量輸入具有方向性，則質量僅輸出用x方向表示。 5.4 建立模型 5.4.1 建模準備工作 模型樁腿與導管架等效厚度計算： 結構模型化的過程中，出于計算簡潔和方便建模的目的，利用抗彎剛度等效原理，將樁腿的壁厚等效到外層的導管架上。 根據(jù)材料力學有： 式5.1 式中，W為構件的抗彎截面系數(shù) 式5.2 對于管單元，慣性矩為 式5.3 —外徑 d—內(nèi)徑 根據(jù)靜力等效原理，要使得構件的最大應力相同，需有： 式5.4 即： 式5.5 在此，設為導管架樁腿導管的外徑，為樁腿的外徑，由于是將樁腿的壁厚等效到外層的導管架上，即 代入上式則有： 將　　　　　　　　　　代入可得： 由＝ 式5.6 代入構件尺寸有 則可得到等效壁厚為 t＝(1340—1236.78)/2＝51.6mm 5.4.2 建模步驟 １）創(chuàng)建關鍵點 根據(jù)原始圖紙，計算初各關鍵點坐標，并將其輸入ANSYS程序。本平臺幾何模型共有關鍵點75個，一些關鍵點是由ANSYS程序“Preprocessor →Modeling →Move/Modify →Keypoints”操作鏡像而來。部分關鍵點坐標如下表所列，由于其數(shù)量較多，其余不一一列舉。 表5.1 部分關鍵點坐標 NO. X，Y， Z LOCATION 1 5.950000 5.450000 -9.000000 2 -5.950000 5.450000 -9.000000 3 -5.950000 -5.450000 -9.000000 4 5.950000 -5.450000 -9.000000 5 5.100000 4.750000 -2.000000 6 -5.100000 4.750000 -2.000000 7 -5.100000 -4.750000 -2.000000 8 5.100000 -4.750000 -2.000000 9 4.650000 4.150000 4.000000 10 -4.650000 4.150000 4.000000 ２）通過關鍵點創(chuàng)建直線，完成幾何模型的建立 通過ANSYS程序“Preprocessor →Modeling →Create →Lines →lines →Straight lines”的操作，在相關的關鍵點之間創(chuàng)建直線，則生成下圖。 圖5.5 平臺幾何模型 ３）定義材料及單元屬性 新建單元類型： 該模型共包括mass21、pipe16、pipe59、beam188四種單元類型，此操作的菜單路徑為“Preprocessor →Element Type→Add/Edit/Delete… →Add… →select” 定義單元實常數(shù)： 選擇菜單路徑“ Preprocessor →Real Constants… →Add… →select Type * → OK→input ”，選擇需要添加實常數(shù)的單元類型。 本平臺模型共創(chuàng)建了6種單元常量（Set 1~Set 6），各單元常量及其單元屬性如下。 表5.2 Pipe 59及pipe 16實常數(shù)輸入 編號 Set 2 Set 3 Set 4 Set 5 Set 6 單元類型 Pipe 59 Pipe 59 Pipe 59 Pipe 59 Pipe 16 外徑DO(m) 1.34 1.2 0.61 0.4 1.2 壁厚TWALL(m) 0.0516 0.03 0.02 0.02 0.03 表5.3 Pipe 59的其他常數(shù)值 拖曳力系數(shù)CD 慣性系數(shù)CM 附加質量比CI 浮力比CB 內(nèi)部流體密度DENSO 1.4 2 1 1 1028 表5.4 Mass 21實常數(shù) X方向質量MASS-X X方向質量MASS-X X方向質量MASS-X 50000 50000 50000 (建模過程中，在四根樁腿頂部各加一個50噸的質量單元，用于模擬總重為200噸的平臺上部模塊) 表5.5 Beam 188實常數(shù)輸入（甲板主梁采用工字梁，建模時選取工字截面） 工字梁下翼板寬度W1 工字梁上翼板寬度W2 工字梁高度W3 0.3 0.3 0.6 工字梁下翼板厚度t1 工字梁上翼板厚度t2 工字梁腹板厚度t3 0.025 0.025 0.025 ４）定義材料屬性 定義材料的彈性模量和泊松比： 平臺所選用鋼材的彈性模量和泊松比分別為2.06×（單位為N/㎡）和0.3，菜單路徑 “Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX：2.06e11， PRXY：0.3 → OK” 定義材料的密度： 鋼材的密度為7800kg/，菜單路徑為““Preprocessor →Material Props →Material Models →Structural →Density→input DENS：7800→ OK” 定義water table： 根據(jù)所給出的環(huán)境條件，創(chuàng)建water table表格。 相關參數(shù)如下： KWAVE(波浪理論) 2 DEPTH（水深） DENSW（海水密度） θw（波向角） 9 1028 － 定義海流： Z(j)（海流深度） W(j)（海流速度） θd(j)（流向角） 0 1.6 － 定義波浪： A(i)（波高） τ(i)（周期） φ(i)（相位角） 6 8.2 － 其中，為了保守起見，θd(j)（流向角）和θw（波向角）取同樣大小，兩者的基線均為X軸正向。 Water table 相關說明 Pipe 59單元波浪載荷計算有四種波浪理論（KWAVE）可供選擇，見表5.6： 表5-6 波浪理論選項 波浪理論選項（KWAVE） 相應的波浪理論 0 深度衰減經(jīng)驗修正的微幅波理論 1 Airy波理論 2 Stokes五階波理論 3 流函數(shù)波浪理論 ５）定義重力加速度 　　菜單路徑為“Solution →Define Loads →Apply →Structural →Inertia →Gravity →Global”，在ACELZ（Z向加速度）項中填入9.8。 ６）劃分網(wǎng)格 １．把單元類型賦給幾何模型 對于導管架、樁腿、以及甲板主梁，選擇“Preprocessor →Meshing →Mesh Tool→(Size Controls) lines： Set →”，拾取相應的線定義屬性。 對于質量點單元，選擇“Preprocessor →Meshing →Mesh Tool→(Size Controls) key pionts： Set →”定義屬性。 ２．選擇自由網(wǎng)格劃分 在“Mesh Tool…”對話框中，選擇復選框Smart Size，激活其下的滾動條，調整劃分網(wǎng)格的精度為“6”。選擇“mush”分別對點和線進行網(wǎng)格劃分。 ７）施加邊界條件 菜單路徑“Solution →Define Loads →Apply →Structural →Displacement →On key pionts ”，拾取模型最下面的四點，定義為固支，即UX=UY=UZ=ROTX=ROTY=ROTZ=0. 注：模型將樁腿簡化為泥面以下6倍樁徑固支。 最終生成模型如圖5.6： 圖5.6 平臺模型 第六章 平臺整體結構靜力分析 6.1 平臺位移 使用ANSYS結構分析軟件4.5節(jié)中各工況進行結構靜力分析，由后處理器可輸出結構位移和應力。各工況平臺最大節(jié)點位移見下表，各工況位移圖見下。 表6.1 各工況下結構最大位移（m） 工況 主要載荷和作用方向 X方向位移 Y方向位移 合位移 波流工況1） 波流力+風力0° 0.50874E-02 -0.18288E-02 0.53106E-02 波流工況2） 波流力+風力90° 0.35284E-02 0.22232E-02 0.38146E-02 海冰工況1） 冰流力+風力0° 0.24684E-01 0.42651E-03 0.24863E-01 海冰工況2） 冰流力+風力90° 0.35935E-02 0.20161E-01 0.20443E-01 各工況位移圖如下： 圖6.1 波流工況1）位移圖 圖6.2 波流工況2）位移圖 圖6.3 海冰工況1）位移圖 圖6.4 海冰工況2）位移圖 由以上結論可得： 平臺頂部的最大水平X位移發(fā)生在海冰1工況下，為2.47 cm； 平臺頂部的最大水平Y位移發(fā)生在海冰2工況下，為2.02 cm； 平臺頂部的最大水平合位移發(fā)生在海冰1工況下，為2.49 cm。 6.2 結構強度校核 使用ANSYS結構分析軟件可對平臺的各工況進行結構靜力分析，后處理器可輸出結構位移和應力。各工況平臺最大單元應力見下表，各工況應力云圖見下。 表6.2 各工況平臺最大單元應力 工況 主要載荷和作用方向 最大單元應力(Pa) 波流工況1） 波流力+風力0° 0.305E+08 波流工況2） 波流力+風力90° 0.269E+08 海冰工況1） 冰流力+風力0° 0.121E+09 海冰工況2） 冰流力+風力90° 0.114E+09 圖6.5 波流工況1）應力云圖 圖6.6 波流工況2）應力云圖 圖6.7 海冰工況1）應力云圖 圖6.8 海冰工況2）應力云圖 部分節(jié)點單元編號如圖6.9： 圖6.9 部分節(jié)點單元編號 各工況下部分構件單元應力比如下： 表6.3 波流工況2）各構件單元應力比 構件名稱 型常數(shù) 單元編號 屈服強度(MPa) 計算應力(MPa) 應力比 導管架導管 ?1340×25 445 315 7.33 0.023 428 315 2.94 0.009 438 315 5.83 0.018 407 315 2.62 0.008 樁腿 ?1200×30 1013 315 8.49 0.027 1012 315 21.27 0.068 1004 315 9.82 0.031 1005 315 19.53 0.062 EL+4.00外橫撐 ?610×20 319 315 1.75 0.006 347 315 3.29 0.010 EL-2.00外橫撐 ?610×20 193 315 5.84 0.019 207 315 22.34 0.071 209 315 22.66 0.072 223 315 4.99 0.016 EL-9.00外橫撐 ?610×20 65 315 0.27 0.001 95 315 1.91 0.006 斜撐 ?610×20 549 315 12.17 0.039 564 315 13.84 0.044 477 315 9.59 0.030 463 315 7.74 0.025 EL+4.00內(nèi)橫撐 ?400×20 609 315 1.17 0.004 623 315 1.79 0.006 EL-2.00內(nèi)橫撐 ?400×20 654 315 4.37 0.014 668 315 1.34 0.004 EL-9.00內(nèi)橫撐 ?400×20 729 315 3.12 0.010 743 315 1.27 0.004 表6.4 波流工況2）各構件單元應力比 構件名稱 型常數(shù) 單元編號 屈服強度(MPa) 計算應力(MPa) 應力比 導管架導管 ?1340×25 445 315 9.3 0.030 428 315 2.62 0.008 438 315 2.31 0.007 407 315 5.46 0.017 樁腿 ?1200×30 1013 315 4.49 0.014 1012 315 18.64 0.059 1004 315 3.21 0.010 1005 315 20.51 0.065 EL+4.00外橫撐 ?610×20 348 315 3.64 0.012 379 315 4.36 0.014 EL-2.00外橫撐 ?610×20 224 315 3.7 0.012 238 315 3.47 0.011 241 315 5.99 0.019 255 315 2.53 0.008 EL-9.00外橫撐 ?610×20 96 315 3.68 0.012 128 315 1.15 0.004 斜撐 ?610×20 549 315 7.92 0.025 564 315 9.95 0.032 477 315 5.97 0.019 463 315 4.52 0.014 EL+4.00內(nèi)橫撐 ?400×20 609 315 1.66 0.005 623 315 0.89 0.003 EL-2.00內(nèi)橫撐 ?400×20 654 315 4.47 0.014 668 315 0.96 0.003 EL-9.00內(nèi)橫撐 ?400×20 729 315 3.05 0.010 743 315 1.09 0.003 表6.5 海冰工況1）各構件單元應力比 構件名稱 型常數(shù) 單元編號 屈服強度(MPa) 計算應力(MPa) 應力比 導管架導管 ?1340×25 445 315 12.33 0.039 428 315 2.98 0.009 438 315 10.47 0.033 407 315 4.59 0.015 樁腿 ?1200×30 1013 315 73.05 0.232 1012 315 71.59 0.227 1004 315 30.84 0.098 1005 315 63.06 0.200 EL+4.00外橫撐 ?610×20 319 315 3.31 0.011 347 315 5.91 0.019 EL-2.00外橫撐 ?610×20 193 315 32.09 0.102 207 315 118.55 0.376 209 315 118.58 0.376 223 315 23.61 0.075 EL-9.00外橫撐 ?610×20 65 315 2.97 0.009 95 315 4.34 0.014 斜撐 ?610×20 549 315 70.74 0.225 564 315 68.7 0.218 477 315 10.01 0.032 463 315 8.34 0.026 EL+4.00內(nèi)橫撐 ?400×20 609 315 2.55 0.008 623 315 3.82 0.012 EL-2.00內(nèi)橫撐 ?400×20 654 315 5.26 0.017 668 315 4.62 0.015 EL-9.00內(nèi)橫撐 ?400×20 729 315 4.33 0.014 743 315 3.68 0.012 表6.6 海冰工況2）各構件單元應力比 構件名稱 型常數(shù) 單元編號 屈服強度(MPa) 計算應力(MPa) 應力比 導管架導管 ?1340×25 445 315 3.52 0.011 428 315 6.33 0.020 438 315 10.84 0.034 407 315 10.17 0.032 樁腿 ?1200×30 1013 315 40.03 0.127 1012 315 46.32 0.147 1004 315 55.58 0.176 1005 315 86.11 0.273 EL+4.00外橫撐 ?610×20 348 315 5.76 0.018 379 315 8.42 0.027 EL-2.00外橫撐 ?610×20 224 315 23.14 0.073 238 315 33.13 0.105 241 315 112.83 0.358 255 315 19.43 0.062 EL-9.00外橫撐 ?610×20 96 315 3.65 0.012 128 315 4.52 0.014 斜撐 ?610×20 549 315 8.17 0.026 564 315 10.48 0.033 477 315 57.16 0.181 463 315 53