無粘結柔性管抗拉伸層結構分析

楊旭，孫麗萍，李鶴楠

（1哈爾濱工程大學船舶工程學院，哈爾濱150001；2挪威科技大學海洋技術系，挪威特隆赫姆7491）

無粘結柔性管抗拉伸層結構分析

楊旭1，孫麗萍1，李鶴楠2

（1哈爾濱工程大學船舶工程學院，哈爾濱150001；2挪威科技大學海洋技術系，挪威特隆赫姆7491）

柔性管抗拉伸層是復雜的空間螺旋線結構，其結構響應分析對柔性管疲勞分析、強度分析和屈曲分析有重要作用。文章基于曲梁理論，應用斜駛螺旋線假設和測地線假設兩種空間曲線公式，以空間細長桿理論及胡克定律本構方程為基礎，采用格林應變張量與第二Kirchoff應力張量度量，對深海無粘結柔性管抗拉伸層螺旋形鋼纜結構平衡方程進行了推導，編寫了分析程序。利用該程序，分析了抗拉伸層鋼纜在軸對稱載荷下和彎矩作用下的曲率變化和結構響應；同時利用三維直梁有限元模型與曲梁有限元模型建立數值模型，將程序結果與數值模擬模型結果進行了對比，證明了結果可行性。該結果可為柔性管抗拉伸層結構設計提供快速的預估計方法。

柔性管；抗拉伸層；斜駛線；測地線；結構分析

0 引言

立管和海底管道作為海洋石油工業的“生命線”，承擔著油氣水運輸的重要任務。無粘結柔性管可作為柔性立管或海底管道，具有結構布置形式靈活、順應性強、與平臺耦合較弱、安裝與回收成本低等諸多優點，在挪威北海、巴西以及美國墨西哥灣等惡劣深海環境下被廣泛運用。無粘結柔性管為多層復合結構，一般而言從內到外可分為內襯層（Internal Carcass）、內壓防護層（Internal Pressure）、內鎖壓力層（Pressure Armor）、抗磨層（Anti-wear tape）、抗拉層（Tensile Armor）和外保護層（Outer Sheath）；此外在一些層與層之間還存在用聚酯纖維制作的中介保護層。柔性管的截面形式如圖1所示，其各層結構均有特定的用途[1]。

無粘結柔性管抗拉伸層一般有兩層，由多條螺旋形高強度鋼纜繞制而成。在工作狀況下，抗拉伸鋼纜主要用于抵抗軸向拉力以及部分外壓和彎矩。由于結構形式為空間曲線，又存在多種相互接觸以及不同材料，抗拉伸層鋼纜結構分析一直是研究熱點。對于其相關研究，Hector等[2]運用ANSYS軟件模擬了柔性管受拉力與扭矩時的結構響應，以彈簧單元和面對面接觸單元模擬抗拉伸層與各層之間的接觸。將抗拉伸層作為三維直梁單元考慮，并添加大量的接觸單元。這樣做無疑會加大計算時間，同時為單元帶來了初始幾何缺陷。部分研究人員選擇了更適合于螺旋線纜的細長桿理論，將鋼纜假設為沿空間曲線的螺旋梁。Svein[3]運用細長桿理論建立了適于模擬抗拉伸層的新型曲梁單元，進行了軸對稱載荷下的結構分析；并進行了軸對稱載荷與彎矩作用下的疲勞分析。?stergaard[4]假設單根螺旋鋼纜滿足無摩擦的測地線假設，對柔性管在受壓和彎曲過程中可能產生的失穩進行了分析。Yan與Lun[5]通過加入人工系數CP對測地線假設結果進行修正，以此模擬鋼纜在彎曲過程中的摩擦力。

本文基于非線性空間細長桿理論，將螺旋抗拉伸層鋼纜考慮為空間曲線，建立了抗拉伸層解析解模型，運用MATLAB語言編寫了抗拉伸層結構響應分析程序。依據基本假設不同，建立了兩種解析解模型：一種假設變形后的抗拉伸螺旋鋼纜符合變形規律滿足測地線假設（geodesic assumption）；一種假設變形后滿足斜駛螺旋線假設（loxodromic assumption）。求解了在軸對稱載荷下以及受彎矩情況下的解析解。獲得了抗拉伸層在不同載荷下的位移與曲率。利用BFLEX軟件建立了數值模擬模型，與解析解進行了對比，在小變形階段結果較為接近。

1 理論推導

1.1 理論背景與基本假設

本文推導的抗拉伸層運動方程是以細長桿理論為基礎。其基本理論為假設空間結構是一細長體，以弧長s和時間t為基本變量，假設其他變量為其函數，建立以切向，法向，副法向為坐標軸的曲面坐標系，通過求解空間曲線的曲率變化獲得結構響應。該方法適合分析沿空間曲面變化的細長體結構。該理論最早可追溯到Kirchoff的線性細長桿理論，后由Ericksen等[6]發展為非線性理論。Washizu[7]運用格林應變張量與Kirchoff應力張量、胡克定律與虛功原理得到了線性細長體理論下的應力應變。Svein[8]此后將其延伸到大變形小應變范圍，并應用于柔性管的應力分析。本文在其基礎上進行分析。根據相關理論，采用如下基本假設：

（1）抗拉伸層螺旋鋼纜假設為細長體，根據一般柔性管的長度和抗拉伸層尺寸，這一假設是合理的。

（2）每一根鋼纜是獨立的，不考慮其與其他鋼纜間相互作用。這一假設適用于柔性管曲率小于臨界曲率時，符合工程計算的一般情況。當結構變形過大，各鋼纜直接可能產生接觸力。

（3）以一根或幾根螺旋鋼纜的結果代表整個抗拉伸層結構響應，忽略柔性管道的橢圓化。該假設可大幅減少計算時間，同時基本保證計算精度，被眾多研究學者采用。由于柔性管截面中心線對稱，所以該假設在軸對稱載荷下影響可忽略。當考慮彎矩作用時，會對計算結果有一定影響。

（4）不考慮層與層之間的摩擦作用。雖然摩擦作用在受彎矩載荷情況下被普遍認為具有重要作用，但為節約計算成本，本解析解并沒有考慮摩擦作用。可在解平衡方程時進行系數修正，獲得更符合實際的結果。

1.2 空間細長桿應力應變

假設螺旋鋼纜為一空間細長桿。其橫截面上有任意一點P，依據螺旋方向建立局部曲面坐標系，以切向量，法向量，副法線向量為基坐標向量，坐標系如圖2所示。則P點位移方程可按下式表示：

圖2 局部坐標系與總體坐標系[8]Fig.2 Local coordinate system and global coordinate system[8]

其中：ui為沿局部坐標系下位移分量，Gi為沿曲線切線，法線副法線組成的局部坐標系基向量。上角標0表示位移沿中心線變化的相關分量，在下文中將被推導。X1為弧長變量，φ（X2，X3）代表與梁橫截面相關的翹曲方程；依據圣維南假設，β為常量。由Frenet-Serret公式可得：

式中：κi代表螺旋線的曲率分量。

根據Svein[3]的研究，選用格林應變張量作為結構的應變表達式，以細長桿中同一原點建立相應的局部笛卡爾直角坐標系，若將笛卡爾坐標系與曲坐標系基向量分別設為Xi與Yi，則運用張量轉化公式可得：

其中：Eij′表示曲線坐標系下的Green應變張量，Eij表示笛卡爾坐標系下的結果。

對抗拉伸層應用如下假設：

（1）假設鋼纜的截面為中心對稱的，這在現今的柔性管設計中都可以保證。

（2）忽略格林應變中非重要的二階項。

（3）應用歐拉伯努利梁假設（平斷面假設）。

將假設代入（7）式可得：

本文中下角標“，i”表示該量對曲線坐標系下相應基向量求導。ε1表示由于軸向和橫向位移產生的軸向應變，ε2和ε3表示由于中心線旋轉產生的應變分量。△κi表示弧長變化dX1產生的i方向曲率增量。由歐拉伯努利梁基本假設，消去θ2和θ3。G為基本張量行列式，表達式如下：

各變量按中心線位移與空間曲線曲率表示為：

采用二階Kirchoff應力張量S與上文的格林應變張量E作為應力與應變測量值，假設相互關系滿足胡克定律，則可得本構關系方程。至此，即可運用鋼纜的沿弧長的曲率變化值求出其應變值，再推導出應力值。

1.3 虛位移原理與平衡方程

本文采用虛位移原理進行求解。虛位移原理為彈性力學中的經典原理，即對結構加以微小虛擬位移獲得平衡方程。根據虛位移原理可推出結構內力虛功為：

其中：Wi表示內力虛功，A為細長螺旋桿截面積，β表示螺旋桿與管中心線產生的相對位移，d為修正摩擦系數，c為表示兩者相互摩擦關系張量。在本文中由于忽略摩擦作用，因此都可設為0。曲率變化值△κi=△κiu+△κip，△κip為螺旋鋼纜按照假定曲線（斜駛線或測地線）變化形成的曲率變化分量，△κiu為由于層與層間相互作用使鋼纜在假定曲線基礎上產生相對滑移，由此產生的曲率變化。本文中假設△κiu=0。據此公式，利用三點高斯積分公式即可獲得某一段細長桿件的內力做功。

取細長桿件的一個微小分段X1進行受力分析，如圖3所示，由此可獲得力與力矩平衡方程如下：

利用（5）式即可將其改寫為六個平衡方程。

1.4 運動限制方程

本文假設鋼纜的變化曲面為以管中心為中心線的圓柱曲面，且半徑不變。該假設在進行小變形受拉扭載荷時影響不大；在進行屈曲分析以及考慮兩端作用時會造成結果不準確，但為了便于獲得鋼纜結構響應，首先進行以上假設，在今后研究中再進行修正。由此可簡化全三維結構，為確保鋼纜移動僅

在曲面上，需加入運動限制方程：

其中：κh為圓柱表面橫向曲率分量，由歐拉方程得：

式中：κr，κc分別為中心曲面在徑向和軸向上的曲率分量。

1.5 測地線假設

測地線來源于地理學，數學上表示在空間曲面上最短的一條線，可視作曲線在彎曲空間的推廣。在結構分析中，若假設結構在曲面上沿測地線假設運動，則表示結構僅被限制于曲面內，即無摩擦作用和其他限制。以曲面中心線角度φ以及螺旋線繞中心線轉過的角度ν和中心線曲率半徑ρ表示中心線笛卡爾坐標系分量如下：

則弧長S對于笛卡爾坐標系分量的微分可表示為

若假設螺旋形鋼纜追隨測地線變化，則其長度可由微分幾何獲得。將空間曲線弧長dS積分并轉換為上述變量表示，加入Lagrange乘子λ，利用最優化理論解得

其中：C為積分常量。不妨設ν=π/2，則

由此可解得C。將φ，1以R/ρ為展開項進行泰勒展開，忽略高階項，沿弧長重新積分可得

另有沿曲線坐標關系方程

將（28）、（29）式結合，解得

其中第二項即為ν的變化值，則總相對位移可得

各方向位移分量為

上式表明了當曲面中心線存在一定的曲率半徑時，若假設鋼纜沿測地線變化，則鋼纜會產生側向滑移與軸向滑移，其滑移值與曲率半徑，螺旋角和所在方位相關。整個軸向滑移從受壓縮到受拉伸不斷變化。

事實上，若鋼纜在圓柱表面符合測地線假設移動，其會產生一定的扭轉角，產生一定的應力；由于前文中假設鋼纜法向量與曲面法向量相同，因此該應力會導致所求結果與實際結果存在一定偏差。Svein[3]討論了這種偏差的影響，指出當結構受到軸向拉力時，這種偏差會顯著減小0.1°以下，可以忽略。結構所需的軸向拉力大小取決于結構形式以及螺旋線布置角，若布置角大于30°，則所需軸向拉力大幅度減小。因本文考慮結構受到拉力與彎矩作用，故忽略這種影響，認為運動限制方程有效。

綜合上述公式解得滿足測地線假設的螺旋鋼纜曲率增量為：

若柔性管僅受拉壓載荷，則曲率半徑無窮大，管道模型為軸對稱模型。由（12）-（17）式可得螺旋線軸向應變與管中心軸軸向應變方程以及軸對稱載荷下曲率變化方程：

其中：εz表示柔性管中心軸的軸向應變。再結合空間幾何關系，最終可得曲率變化：

1.6 斜駛線假設

斜駛線的概念同樣來源于地理學，表示該空間曲線與曲面坐標軸所成角度不變，即螺旋角α不變。斜駛線與測地線在球體上的表示如圖4所示。在物理學上，可視為空間曲線結構與所在曲面無相對滑移，即摩擦力無窮大。當空間曲面彎曲變化時，由于受到軸向應力作用，原坐標系基向量G1的模也會隨之變化，設變化后的基向量為g1，由（6）式可得格林張量值，代入基向量中可得：

圖4 測地線（左）與斜駛線（右）Fig.4 Geodesic curve(left)and loxodromiccurve(right)

利用測地線假設與斜駛螺旋線假設下的曲率變化方程，可求解受到環境載荷之后的抗拉伸層鋼絲新曲率值。再由曲面坐標系下的應力應變關系，求解出各點應力，通過三點高斯積分，即可獲得抗拉伸層鋼絲彎矩和剪力。將整個過程不斷迭代，并應用新的螺旋角，即可得到抗拉伸層鋼纜結構響應的準靜態公式。

2 數值模擬實例

根據前文所述的公式，以MATLAB程序語言編寫了結構分析程序代碼TensileArmorV1.M。該程序以準靜態線性方法，分析抗拉伸層鋼纜的曲率變化，再通過上文公式導出彎矩和剪力值。程序中忽略摩擦作用以及鋼纜與管其它層間作用，假定鋼纜符合上文所述各種假設。程序流程圖如圖5所示。

圖5 結構分析程序流程圖Fig.5 Flow diagram of program TensilearmorV1.m

圖6 ABAQUS有限元模型Fig.6 Finite element model of ABAQUS

本文選擇一根柔性管作為實例分析。其抗拉伸層直徑為100 mm,抗拉伸層橫截面為長方形，厚3 mm，寬10 mm，布置角度為30°。材料為高強度碳鋼，本文假設為滿足胡克定律的彈性材料，楊氏模量E=211 GPa,屈服應力1 100 MPa，泊松比0.3。文中選擇4倍螺距長度進行分析。由于程序僅分析單根抗拉伸鋼纜作用，本文將其它層作用假設為管道位移邊界條件施加到抗拉伸層中；管道模型曲率按從0 m-1到0.1 m-1循環變化。

為對比分析不同模型差異，本文運用ABAQUS通用有限元軟件隱式計算模塊，建立了以三維直梁單元（B31）為基礎的有限元模型。模型的材料，截面參數與邊界條件與前文相同。采用SHELL單元代表整個柔性管結構。同樣采用一根等效抗拉伸層鋼纜進行分析。鋼纜按長度方向被分為100個單元，每個單元節點與對應的表面運用Tie連接，以此實現斜駛螺旋線假設。加載步驟為先勻速緩慢加入軸向拉力，再緩慢施加彎矩循環。其中軸向拉力部分為靜態分析，彎矩循環部分為動態分析。所得模型如圖6所示。

BFLEX管道分析軟件是MARINTEK水池開發的，依據空間曲梁原理與斜駛螺旋線假設以及非線性有限元框架編寫的有限元結構分析軟件。采用該軟件建立類似模型與本文結果和ABAQUS模型進行對比。BFLEX模型采用100個H353單元，該單元為曲梁單元。同時考慮管道與抗拉伸層的摩擦，摩擦系數設為0.4，采用較高摩擦系數近似模擬斜駛螺旋線的摩擦無限大假設。其他參數與ABAQUS模型相同。所得模型如圖7。

圖7 BFLEX有限元模型Fig.7 Finite element model of BFLEX

3 結果與討論

3.1 軸對稱載荷結果分析

當柔性管僅受軸向拉力時，柔性管抗拉伸鋼纜在局部坐標系下的曲率變化結果如圖8和圖9所示。軸向拉力以指定位移勻速施加在管中心軸上，從0勻速增加至1%管長。從圖中可得，切向與縱向曲率變化隨軸向拉力基本呈線性變化，三個模型趨勢相同，結果近似，證實了本文模型的可行性。采用直梁單元的ABAQUS模型曲率變化更大。BFLEX模型與本文模型由于基本方程相同，所以變化不大。差異主要是由于數值模型中采用非線性計算方法，同時布置角α會隨軸向應變不斷變化；而本文程序采用線性計算，同時假設α為一定值。徑向曲率變化值如圖10，按照BFLEX與本文假設，徑向曲率變化率為0。由于直梁模型為近似模擬螺旋曲線，故在ABAQUS模型中存在較小的曲率變化，但數量級為10-5m-1，可以忽略不計。圖11表示了抗拉伸鋼纜在拉力作用下的軸向力。圖12表示了BFLEX模型當拉力最大時抗拉伸層鋼纜關于管道中心軸的彎矩分布。ABAQUS模型與本文程序所得彎矩分布與該圖類似。本文選擇彎矩最大處進行對比分析，結果如圖13。彎矩與軸力圖表明以格林應變張量與第二Kirchoff應力張量度量并運用坐標系轉化這一方法是切實可行的。

對于計算效率的比較，本文選擇一臺Intel I7 4核處理器，32 GB內存進行計算，本文模型計算時間約10 s，BFLEX模型計算時間約100 s，ABAQUS模型計算時間約為600 s。由于采用了針對柔性管的坐標系與單元，本文模型計算時間大大減少。

圖8 軸向拉力下螺旋鋼纜縱向曲率變化Fig.8 Axial curvature increment in tensile load condition

圖9 軸向拉力下螺旋鋼纜橫向曲率變化Fig.9 Transverse curvature increment in tensile load condition

圖10 軸向拉力下螺旋鋼纜徑向曲率變化Fig.10 Radical curvature increment in tensile load condition

圖11 軸向拉力下螺旋鋼纜軸力Fig.11 Axial force of tensile armor wire

圖12 關于管中心軸彎矩分布圖Fig.12 X(m)versus Moment-Z(N*m)

圖13 關于中心軸最大彎矩對比圖Fig.13 Maximum moment about center core

3.2 彎曲載荷下結果分析

對柔性管兩端施加彎矩使其曲率從0 m-1到0.1 m-1緩慢增加，以近似于三角函數循環變化。在彎矩施加前對柔性管施加0.15%的拉應變以保證整個彎矩作用時柔性管狀態穩定，避免發生屈曲。在Bflex模型中選擇摩擦系數miu分別為0與0.4用以分別對比測地線假設與斜駛螺旋線假設的情況。在ABAQUS模型中運用tie連接將抗拉伸鋼纜與表面相連，以此來模擬斜駛螺旋線假設。同時在程序中分別算出斜駛螺旋線與測地線假設結果。

本文以曲率達到最大值時結果進行分析。數值模擬所得抗拉伸鋼纜相對于中心管軸向位移如圖14。軸向位移隨著螺旋的循環不斷循環變化，周期為一個螺距。在螺旋開始端為0，在1/4螺距處達到最大值。根據假設，斜駛線與ABAQUS模型的橫向位移為0 mm。摩擦系數為0.4的BFLEX模型較為接近斜駛線假設模型，幅值為0.4 mm。而摩擦系數為0的BFLEX模型則比測地線假設略大，幅值分別為2.5 mm與2 mm。橫向相對位移沿管長變化如圖15所示，所得結果同樣以單螺距為周期循環變化。斜駛線模型，高摩擦系數模型與ABAQUS模型結果基本相同，符合假設預期。

圖14 縱向相對位移沿管長分布圖Fig.14 Pitch times versus relative displacement(mm) in axial direction

圖15 橫向相對位移沿管長分布圖Fig.15 Pitch times versus relative displacement(mm) in transverse direction

圖16 扭轉曲率相對變化值沿管長分布圖Fig.16 Pitch times versus relative twist curvature(m-1)

圖17 橫向曲率相對變化值沿管長分布圖Fig.17 Pitch times versus relative transverse curvature(m-1)

圖16表示扭轉方向曲率相對變化值Δκ1沿管長方向分布圖。測地線曲率變化值呈余弦函數分布，周期為單位螺距，幅值約為0.05 m-1；與無摩擦BFLEX模型近似。斜駛線模型與高摩擦模型近似，相位角與測地線假設相反，幅值較小，約為0.02 m-1。ABAQUS模型變化值呈正弦函數分布，周期與其他模型相同，幅值略大，約為0.12 m-1。產生差異原因可能是ABAQUS模型采用直梁單元，其初始梁單元方向僅能代表該微段梁單元的平均方向，而不能代表起始點準確方向，因此在彎曲作用下曲率變化較為敏感。由于結果顯示為曲率相對變化值，故對實際梁運動影響不大。但在屈曲分析或疲勞分析中，需考慮管道受長時間多次彎曲循環作用，該微小差異可能導致直梁有限元結果與其他基于空間細長體結果產生不同。圖17表示橫向曲率相對變化管長分布。由于邊界作用，ABAQUS與BFLEX模型在邊界上出現小幅波動，但對結果影響較小。測地線結果與0摩擦模型相似；斜駛線模型與高摩擦模型相似，符合假設預期。所有細長桿模型結果均類似于以單螺距為周期的余弦函數。由于測地線假設摩擦為0，因此橫向曲率變化較斜駛線假設大，幅值分別為0.2 m-1與0.04 m-1。直梁模型結果曲率變化平均值約為0 m-1，與斜駛線近似，不同之處在于結果周期性不明顯，分布不穩定，這與直梁模型近似表示空間曲線結構相關。圖18為徑向曲率變化值對比圖。結果同樣較為符合，測地線曲率變化較斜駛線大。ABAQUS模型結果與斜駛線模型近似，但分布不穩定，原因與上文相同。

上述結果表明測地線假設可代表摩擦為0時抗拉伸螺旋管的位移變化；而斜駛線假設則可基本代表高摩擦系數時抗拉伸層位移變化。工程實際中摩擦系數在0與無限大之間；因此可通過該程序的兩種假設模型快速估算出抗拉伸鋼纜運動響應范圍。在同上文相同計算機上，本文程序計算時間約為100 s，BFLEX約為900 s，ABAQUS達到2 000 s；表明本文程序計算效率較高。

4 結論

本文基于測地線假設與斜駛螺旋線假設，推導了柔性管抗拉伸層螺旋鋼纜在軸對稱載荷以及循環彎矩作用下的結構曲率變化方程組。同時，應用格林應變張量，二階Kirchoff應力張量以及胡克定律，將格林應變張量在曲面坐標系下表示，并轉化為局部坐標系，得出抗拉伸層鋼纜彎矩，軸力公式。編寫了線性準靜態計算程序，選擇柔性立管實例進行了分析。運用有限元軟件建立了多個曲梁單元模型與直梁單元模型，與計算結果進行了對比，結論如下：

在柔性管受軸對稱載荷時，計算結果表明局部坐標系x和y方向曲率變化值與拉應變呈線性關系。局部坐標系徑向曲率基本為0。與通用直梁有限元模型以及曲梁單元模型對比，結果基本相同。本文假設成立，計算程序可行。

在受彎曲載荷時，斜駛線假設模型與高摩擦系數模型高度近似，證實了斜駛線假設符合高摩擦系數時螺旋鋼纜運動響應。測地線假設模型與無摩擦模型結果近似，表面測地線假設模型可用于預測低摩擦系數時抗拉伸層運動響應。結果表明摩擦力對抗拉伸鋼纜在彎曲作用下的結構響應有較大影響。直梁模型在位移上與斜駛線模型十分接近，在曲率變化上平均值與斜駛線模型接近，但分布規律較為不同；其原因主要是基本假設與單元特性造成偏差。在短時間彎曲分析中差異可忽略不計。在多次彎曲循環分析中，該差異可能引起曲梁模型與直梁模型結果不同，有待進一步驗證。

在兩種載荷分析中，本文模型都體現了較高的運行效率。工程上本程序可應用于軸對稱載荷下抗拉伸鋼纜設計;并快速分析抗拉伸鋼纜在彎曲載荷下的運動響應范圍，為抗拉伸鋼纜設計提供參考。

在今后研究中將根據實驗結果，增加修正系數建立結構響應與摩擦系數關系從而準確模擬抗拉伸層。同時，非線性動態計算方法，非線性材料模型以及螺旋線布置角的實時變化將被加入到今后的程序中。抗拉伸層與其他層之間的接觸作用在本文中沒有考慮，也將在進一步的研究中體現。

[1]楊旭,孫麗萍,艾尚茂.深水無粘結柔性管抗拉伸層屈曲問題研究進展[J].海洋工程,2013:31(1):99-106. Yang X,Sun L P,Ai S M.Recent research of deepwater unbounded flexible pipe tensile armor layer’s buckling behaviors[J].Ocean Engineering,2013:31(1):99-106.（in Chinese)

[2]Hector M,Jose S.,Carlos M,Ney R.On the coupled extensional torsional response of flexible pipes[C]//Proceedings of the ASME 2009 28th International Conference on Ocean,Offshore and Arctic Engineering.Hawaii,2009:1-13.

[3]Svein S.On stresses and fatigue in flexible pipes[D].Norway:Norwegian University of Science and Technology,1993: 2.1-2.42.

[4]?stergaard N H.On lateral buckling of armouring wires in flexible pipes[D].Denmark:Aalborg University,2012:1-10.

[5]Yangqiu Z,Lun Q.Numerical model to simulate tensile wire behavior in unbonded flexible pipe during bending[C]//Proceedings of the 26th International Conference on Ocean,Offshore and Arctic Engineering.USA,2007:1-13.

[6]Ericksen J,Truesdell C.Exact theory of stress and strain in rods and shells[J].Arch Rat Mech Anal,1957(1):295-323.

[7]Washizu K.Some considerations on a naturally curved and twisted slender beams[J].Journal of Math Phys,1964(43): 111-116.

[8]Svein S.Theoretical and experimental studies of stresses in flexible pipes[J].Computers and Structures,2011(89):2273-2291.

Analysis of tensile armor layer structural behaviors in flexible pipe

YANG Xu1,SUN Li-ping1,LI He-nan2
(1 Shipbuilding Engineering College,Harbin Engineering University,Harbin 150001,China; 2 Department of Marine Technology,Norwegian University of Science and Technology, Trondheim Norway 7491,Norway)

Tensile armor layers of flexible pipe are complex helix structures,whose stress analyses are significant to flexible pipe’s fatigue,strength and buckling behaviors.This paper is based on curved beam theory.Equilibrium equations of tensile armor wires in deepwater unbonded flexible pipe were solved with assumptions of loxodromic curve and geodesic curve,based on slender rod theory and elastic constitutive relationship.The curvatures and structural response under axisymmetric load and bend moment were analyzed.A three dimensions beam finite element model and a three dimensions curved beam element model were built by FEM software to compare with the results.The results show that the formulations and codes are reasonable for predicting structural responses of tensile armor wires.This solution provides a fast method to evaluate the stress and curvatures in design of flexible pipe’s preliminary design.

flexible pipe;tensile armor layer;loxodromic curve;geodesic curve;structural analysis

P756.2

A

10.3969/j.issn.1007-7294.2015.04.012

1007-7294（2015）04-0436-11

2014-09-22

FPSO和外輸系統數值模擬及實驗研究（2011ZX05030-006-002）；國家建設高水平大學公派研究生項目；深海工程科學與技術創新引智基地（B07019）

楊旭（1987-），男，博士研究生；孫麗萍（1963-），女，教授，博士生導師，通訊作者，E-mail:sunliping@hrbeu.edu.cn；李鶴楠（1987-），女，碩士。