江 錦，孫華峰，秦立成，任 濤，高 巍

海洋石油工程股份有限公司，天津 300451

本文以渤海某岸電工程5 000 t組塊滑移裝船為例，采用移動荷載對裝船過程中的滑道、滑靴進行瞬態(tài)受力分析。

整個分析的難點在于載荷所加載的區(qū)域隨時間而發(fā)生變化。組塊在裝船過程中，由于牽引絞車的拖拽，重量逐步由碼頭向船上轉(zhuǎn)移。由于船上滑道承受壓力的區(qū)域也是隨著時間變化的，因此可以用施加移動荷載的方法進行瞬態(tài)受力分析?；镜姆治鲇嬎氵^程可以概括如下：其一，建立三維模型；其二，瞬態(tài)分析前處理；其三，瞬態(tài)分析設置及邊界條件；其四，瞬態(tài)分析后處理。利用ANSYS Workbench的Transient Structural模塊可以對這個過程進行仿真模擬。

1 滑道及滑靴受力分析

1.1 裝船過程中滑道受力工況分析

組塊在裝船過程中，滑道受力主要來自于上部組塊的重力。組塊通過滑靴向滑道傳遞重力，滑靴底部為木頭，裝船前一般在滑道表面涂抹黃油以減小摩擦。因此在滑道結(jié)構(gòu)受力分析中，可以忽略相切于上表面的摩擦力，而只分析沿著上表面法向方向的力。隨著組塊逐步前進，滑道上受壓的面積和力的大小也隨時間增加，直至單腿的支撐力全部由陸地轉(zhuǎn)移至滑道。提取SACS結(jié)構(gòu)計算中裝船工況下的支反力，見表1，以其中B2支撐點的最大值12 227.52 kN作為分析的輸入條件。組塊裝船過程如圖1所示，各支撐點編號如圖 2所示[1]。

圖1 組塊裝船過程示意

圖2 支點布置編號

表1 各支點反力

1.2 滑道有限元模型建立

根據(jù)圖紙建立滑靴模型和滑道模型，見圖3。為了校核滑道在滑靴移動過程中的受力狀況，建模時將滑靴擺在初始位置?；ラL10 m，單節(jié)滑道長12 m。

圖3 三維實體模型

利用Workbench進行網(wǎng)格劃分時，系統(tǒng)自動為模型指定了Solid187單元。Solid187單元為高階3維10節(jié)點固體結(jié)構(gòu)單元，具有二次位移模式，可以更好地模擬不規(guī)則模型?；ズ突乐g利用接觸單元進行模擬，采用的單元分別是CONTA174單元和TARGE170單元，見圖4。接觸類型設定為Frictional，摩擦系數(shù)按DNVGL規(guī)范（見表2） 中給定的數(shù)值取值，取0.15。

表2 DNVGL規(guī)范中摩擦系數(shù)取值

圖4 滑靴和滑道接觸

接觸面定義了滑靴和滑道之間的相互作用，在接觸面上采用相同的網(wǎng)格密度有利于分析結(jié)果更快收斂。在網(wǎng)格的局面控制中，單獨增加Contact Sizing控制模型以在接觸面上產(chǎn)生大小一致的網(wǎng)格。對滑靴頂面施加向下的力，在滑道底面施加固定約束，對滑靴的左端面施加一個2.0 m的位移約束，建立如圖5所示的有限元模型。

圖5 三維實體有限元模型

1.3 計算方法及參數(shù)設置

本文研究的計算方法使用了瞬態(tài)動力學分析，目的在于確定在受隨時間變化荷載作用下的結(jié)構(gòu)的動力學響應，分析過程涉及了非線性特征。完全法采用完整的系統(tǒng)矩陣計算瞬態(tài)響應，計算結(jié)果顯示，采用完全法完成瞬態(tài)動力分析，在三種分析方法（即完全法、縮減法、模態(tài)疊加法）中，其功能最強，可包括各類非線性特性，利用ANSYS Workbench進行分析時，需要打開結(jié)構(gòu)大變形的影響[2]。

積分時間步長的大小既影響計算效率，又影響瞬態(tài)動力分析的精度和收斂性。時間步長越小，精度越高。太大的積分時間步長會影響較高階模態(tài)的響應，從而影響計算精度。自動時間步長可按響應頻率和非線性效果自動調(diào)整求解期間的積分時間步長，對于非線性分析，自動時間步長還會適當?shù)卦黾雍奢d，并在達不到收斂時回溯到先前收斂的解。計算時利用Auto Time Stepping激活自動時間步長，具體的時間步設置如表3所示[3]。

表3 瞬態(tài)分析計算參數(shù)設置

1.4 計算結(jié)果分析

在組塊腿支反力和位移約束的作用下，滑靴從滑道的左端滑動到了右端，完成了設定的計算目標，滑靴的位移計算結(jié)果如圖6所示。

圖6 位移計算結(jié)果/mm

分別查看滑靴和滑道的應力響應結(jié)果，滑靴在分析時長期間的應力變化如圖7所示。在分析的開始階段，滑靴的最大應力達到326.48 MPa，但是與相鄰時間點的應力結(jié)果不連續(xù)，圖8中所示最大應力出現(xiàn)在滑靴上表面鋼板拼接的轉(zhuǎn)角處，此處產(chǎn)生了應力集中，但是區(qū)域很小，不應作為校核滑靴主結(jié)構(gòu)受力分析的參考應力。滑靴鋼材屈服強度為355 MPa，應力最大允許值為355 MPa×0.9=319.5 MPa。從圖中可以看出，滑靴受力在平穩(wěn)階段穩(wěn)定在94 MPa，遠小于滑靴最大許用應力[4]。

圖7 滑靴受力時間曲線

圖8 滑靴受力云圖/MPa

單獨分析滑道的受力結(jié)果如圖9所示，在整個過程中，滑道受到的最大應力為200.67 MPa，最大應力的區(qū)域是隨著滑靴的移動而變化的，大致位于滑靴中心位置的下方，這與實際的受力情況也非常吻合。按滑道最大許用應力為355 MPa×0.9=319.5 MPa進行校核，滑道的應力在整個滑動過程中沒有超過其許用應力[5]。

圖9 滑道受力云圖/MPa

2 結(jié)束語

岸電組塊結(jié)構(gòu)質(zhì)量大，滑移裝船仍是一種非常重要的選擇。針對滑移裝船中滑靴和滑道的結(jié)構(gòu)強度校核問題，采用瞬態(tài)動力分析的方法建立有限元模型，并按規(guī)范要求施加接觸摩擦系數(shù)和按實際運動施加位移約束，得到了滑靴和滑道對應不同時間的應力結(jié)果，同時得出以下結(jié)論：第一，經(jīng)過瞬態(tài)分析，測試項目中滑道和滑靴的結(jié)構(gòu)強度滿足裝船要求；第二，瞬態(tài)分析方法用于滑移裝船的結(jié)構(gòu)強度校核具有一定優(yōu)勢，可以在一次分析結(jié)果中給出結(jié)構(gòu)不同區(qū)域的動態(tài)應力響應。