三芯光纖復合海纜拉伸的有限元建模與仿真

盧志飛，林曉波，鄭新龍，呂安強，李世強，張 杰

（1.國網(wǎng)浙江省電力公司舟山供電公司，浙江舟山 316021； 2.浙江舟山海洋輸電研究院有限公司，浙江舟山 316021；3.華北電力大學電子與通信工程系，河北保定 071003）

三芯光纖復合海纜拉伸的有限元建模與仿真

盧志飛1，2，林曉波1，2，鄭新龍1，2，呂安強3，李世強1，2，張 杰3

（1.國網(wǎng)浙江省電力公司舟山供電公司，浙江舟山 316021； 2.浙江舟山海洋輸電研究院有限公司，浙江舟山 316021；3.華北電力大學電子與通信工程系，河北保定 071003）

三芯光纖復合海底電纜（簡稱為“海纜”）在運輸、敷設和運行的過程中經(jīng)常會產(chǎn)生拉伸的機械行為，對其進行拉伸建模仿真，可獲得實體試驗難以得到的應力、應變等數(shù)據(jù)。文章通過建立海纜拉伸有限元模型，模擬了海纜的拉伸過程。首先根據(jù)海纜拉伸的力學特點對結構進行簡化，得到三芯海纜的幾何模型；然后對其施加約束及速度載荷，控制沙漏能，仿真軸向拉伸過程；最后對數(shù)據(jù)進行有效提取和處理，獲得海纜的應力與應變，為研究海纜的力學性能、判斷海纜工作狀態(tài)提供了參考。

三芯光纖復合海纜；拉伸；應力；應變；有限元

0 引 言

隨著我國島嶼間電力網(wǎng)絡的完善、通信需求的增大以及國際間信息化交流的發(fā)展，我國海底電纜（簡稱“海纜”）的用量也日益增加［1］。而海底環(huán)境的復雜性、海纜敷設不當以及人為破壞都可能對海纜的安全運行造成很大的威脅［2］。海纜在受到機械損害時通常會發(fā)生軸向拉伸，這種軸向拉伸是造成海纜機械損壞的主要原因。因此，為防止海纜破壞的發(fā)生，有必要對海纜拉伸的力學過程進行分析，了解拉伸時各層結構組件的變化趨勢。

由于海纜結構復雜，體積、重量大，進行實體試驗困難且成本高，因此，目前的研究都是利用數(shù)學方法進行建模仿真。華北電力大學張旭等人對單內(nèi)鎧光電復合海纜進行了拉伸試驗仿真［3］；陸瑩等人通過建立簡化的XLPE（交聯(lián)聚乙烯）高壓海纜模型，仿真分析了海纜受到外力造成損壞時，內(nèi)部物理量的變化［4］；武漢理工大學的李兢利用ANSYS/LSDYNA（通用顯式非線性動力分析程序）對海纜的實體抗側壓模型展開了研究，對抗鉤掛實體模型進行了分析［5］。以上研究為本文內(nèi)容奠定了基礎。

本文首先根據(jù)海纜拉伸的力學特點對海纜結構進行了簡化，然后通過對海纜施加約束以及速度載荷仿真軸向拉伸過程，最后對數(shù)據(jù)進行有效提取和處理，獲得海纜各層的應力與應變，為研究海纜的力學性能、判斷海纜工作狀態(tài)及安全預警提供參考。

1 海纜拉伸的有限元建模

1.1 海纜幾何模型建立

本文采用ZS-YJQF 4 1型3 6 k V三芯XLPE絕緣光纖復合海纜。該海纜由3個線芯及外層材料構成，結構剖面圖如圖1所示。線芯材料包括銅導體、導體屏蔽、XLPE、絕緣屏蔽、半導電阻水帶、鉛合金護套和高密度聚乙烯護套；外層材料包括填充層、扎帶、黃銅帶、鎧裝墊層、鋼絲鎧裝層和外被層。銅絲按照絞合結構生成銅導體；聚乙烯內(nèi)護套和鋼管組成光單元，兩根光單元對稱分布在填充層內(nèi)；8根通信用普通單模光纖分布在鋼管內(nèi)；外被層涂抹瀝青；鋼絲以層絞的方式構成海纜的鎧裝層。

圖1 三芯光纖復合海纜結構剖面圖

海纜結構復雜，精確地建立海纜的幾何模型十分重要。本文在盡量保證其真實結構的前提下，為降低建模難度、節(jié)約計算時間，按機械結構特點對其進行了簡化，簡化前后海纜的幾何尺寸如表1所示。

表1 海纜幾何尺寸

簡化過程描述如下：海纜拉伸時，由于導體屏蔽、絕緣屏蔽和半導電阻水層的厚度較小，機械性能較差，故將其合并到與之性能相近的XLPE層。填充層的材料性質決定了其僅起層間支撐作用，不參與拉伸受力。由于填充層被線芯分隔成了4部分，導致每部分的邊緣尺寸很小，形狀各異，很難對其進行合適的網(wǎng)格劃分。又由于本文施加的是速度載荷，即海纜拉伸位移僅與時間有關，因此本文將HDPE（高密度聚乙烯）與填充層合并，以增大填充層邊緣尺寸，降低網(wǎng)格劃分難度，提高求解速度。黃銅帶、鎧裝墊層和外被層的厚度較小，可以忽略。最終，三芯海纜簡化為銅導體、XLPE、鉛合金、光單元、填充物和鋼鎧。

三芯海纜的線芯、光單元和鎧裝鋼絲均是絞合結構，不能直接生成體。本文采用點、線、面和體的思路［6］，先建立三個線芯和光單元，再建立填充層和鎧裝，如圖2所示。海纜在拉伸過程中，每個節(jié)距變化基本相同，為了節(jié)約計算時間，同時考慮模型求解的端部效應，本文取1.33個節(jié)距，最終海纜長度為187.2 cm，建成的海纜幾何模型如圖3所示。

圖2 海纜端面和相應的螺旋線示意圖

圖3 海纜幾何模型

1.2 海纜網(wǎng)格劃分與創(chuàng)建有限元模型

銅導體、XLPE、鉛合金、光單元、填充物和鋼鎧都是三維立體單元，在拉伸過程中會產(chǎn)生位移、速度和加速度的變化。SOLID164實體單元是ANSYS/ LS-DYNA的8節(jié)點六面體單元，包含位移、速度和加速度等自由度，支持所有許可的非線性特性，并支持單點積分，適合海纜拉伸力學特性分析，因此，本文選擇了SOLID164單元。

銅導體、XLPE、鉛合金、光單元、填充物和鋼鎧用的是與應變率無關的各向同性材料，而軟件提供了3種與應變率無關的各向同性材料模型，這3種材料模型的應力-應變曲線都是用兩條直線表達，根據(jù)硬化的假設、溫度是否相關來區(qū)分。本文使用隨動硬化假設，且與溫度無關，認為二次屈服在2σy時出現(xiàn)，選取了雙線性隨動強化模型。

由于三芯海纜結構的特殊性，自動化分會得到較差的網(wǎng)格，導致計算精度降低，甚至不能完成計算。本文采用映射法實現(xiàn)銅導體、光單元和鋼鎧的網(wǎng)格劃分，通過控制體端面上的線及軸線的尺寸來完成，而XLPE、鉛合金和填充物則掃掠完成。最終，整體網(wǎng)格單元數(shù)量約為2.5萬，在保證計算精度的前提下，有效地提高了計算速度。

1.3 載荷施加、約束與接觸控制

在海纜拉伸模型中，將海纜的左端面定義成組件，定義不同的時間以及對應速度載荷值的數(shù)組參數(shù)，最后將載荷施加到該組件上。

在做拉伸實驗時，將海纜的左端面固定，即對該截面上的所有節(jié)點施加所有方向上的約束。由于三芯海纜多是絞合層，拉伸會產(chǎn)生纜體的扭轉。若不加控制扭轉角度太大，會導致海纜發(fā)生嚴重穿透，影響實驗的準確性。實際中的海纜很長，中間的一小段被拉伸后，通過對兩側海纜的約束使纜體扭轉角度較小，可忽略不計。本文對右端面進行零約束處理，也就是控制右端面在X、Y方向上的位移。

海纜在拉伸過程中，層與層的接觸會導致穿透，影響實驗的準確性。而在接觸分析中，由于海纜拉伸過程的復雜性，判斷接觸方向很困難，因此本文使用了接觸中的自動接觸。

1.4 沙漏能控制

三芯海纜結構復雜且特殊，拉伸時網(wǎng)格會產(chǎn)生大變形，嚴重時變形呈現(xiàn)鋸齒形網(wǎng)格，出現(xiàn)零能模式，這會導致從一開始沙漏能便大于內(nèi)能的10%。為確保分析的準確性，本文采用了盡可能使用均勻網(wǎng)格劃分、避免在單點載荷、全局調(diào)整關鍵字* CONTROL_HOURGLASS（設置整體的沙漏公式和沙漏系數(shù)）、控制參數(shù)為IHQ（定義沙漏粘性屬性）與QH（定義沙漏系數(shù)）等方式控制沙漏能。全局調(diào)整后XLPE、鉛合金、光單元和填充物的沙漏能基本控制在1%以內(nèi)。

2 有限元仿真結果與分析

2.1 仿真結果評價

求解海纜拉伸的有限元建模，得到海纜在拉伸后的位移圖，如圖4所示。海纜在拉伸過程中，左端面的位移最大，與本文所施加的載荷產(chǎn)生的位移相符。海纜的右端面沒有位移，符合施加的端面約束（ALL DOF）。海纜在受到軸向速度載荷后，整體直徑變小。

圖4 海纜拉伸位移圖

另外，沙漏能是衡量仿真正確性的重要指標，準確的仿真應保證沙漏能不超過內(nèi)能的10%［7］。本文中的沙漏能控制在內(nèi)能1%以內(nèi)，保證了實驗的準確性。內(nèi)能-沙漏能時間曲線如圖5所示。

圖5 內(nèi)能-沙漏能時間曲線

2.2 仿真數(shù)據(jù)分析

海纜結構復雜，各層材料參數(shù)各不相同，因此在拉伸過程中各層的應變和應力也不同，提取銅導體和光單元的應變，如圖6所示。

圖6 銅導體和光單元的應變曲線

由圖可知，拉伸過程中，導體和光單元兩端的應變變化幅度比較大，中間部分應變分布基本保持水平，這是由于本文施加了勻速載荷。在模型端部，載荷直接作用，因此會出現(xiàn)應力波動，而傳遞至中部后應力變化較平緩，因此應變分布均勻，這與海纜拉伸的實際相符。海纜在右端面固定后，施加速度載荷，海纜每一層在不同位置拉伸位移一樣，故應變也應一樣，隨著位移的增加，海纜應變增加，與圖6一致，符合實際的海纜拉伸模型。

銅導體和光單元的應力-應變曲線如圖7所示。從圖中可以看出，在拉伸過程中，銅導體的應力隨應變的增大而上升，在0～0.003 s期間，銅導體的應變隨應力以斜率6.7E5增加，而在0.003 s后，斜率比之前降低，但仍然是線性增長。光單元的應力也隨應變的增大而增大，在0～0.05 s期間，斜率較大，而在0.05 s后，斜率變小，應力隨應變變化減慢，符合彈塑性材料特性。光單元內(nèi)復合有光纖，當光單元的應變增加到足以消耗掉光纖的余長后，光纖開始受力，并產(chǎn)生應變。

圖7 銅導體和光單元的應力-應變曲線

3 結束語

本文對海纜拉伸的機械行為進行了有限元建模和仿真。通過分析實驗數(shù)據(jù)可以看出，海纜不同層的應變、應力不一樣。拉伸過程中，銅導體和光單元的應力在開始階段都隨應變的增大而線性增大，在應變達到一定數(shù)值后，應力隨應變的變化減緩，即斜率降低，符合材料特性。本文提供了三芯海纜拉伸有限元模型的建立和求解方法，為進一步研究海纜力學性能、判斷海纜工作狀態(tài)提供了參考。

［1］ 顧玉娟.海底電纜的安全性探討［C］//救撈專業(yè)委員會2003年學術交流會論文集.上海：中國航海學會救助打撈專業(yè)委員會，2003：126-127.

［2］ 張國光.海纜安全及其施工埋設技術研究［J］.海洋技術，1992，（1）：64-73.

［3］ 張旭，尹成群，呂安強，等.光電復合海纜有限元建模與仿真［J］.計算機仿真，2013，（2）：120-124.

［4］ 陸瑩，黃輝.基于分布式光纖傳感技術的高壓海纜外力損壞仿真［J］.電氣技術，2012，（12）：87-89.

［5］ 李兢.海底光纜力學特性有限元分析［D］.武漢：武漢理工大學，2009.

［6］ 呂安強，李永倩，李靜，等.用光纖應變判斷光電復合海纜錨害程度的有限元分析法 ［J］.電工技術學報，2014，（11）：261-268.

［7］ 白金澤.LS-DYNA 3D理論基礎與實例分析［M］.北京：科學出版社，2005：167-168.

Finite Element Modeling and Simulation of Three-core Optical Fiber Composite Submarine Cable Stretching

LU Zhi-fei1，2，LIN Xiao-bo1，2，ZHENG Xin-long1，2，LüAn-qiang3，LI Shi-qiang1，2，ZHANG Jie3
（1.State Grid Zhoushan Electric Power Supply Company of Zhejiang Power Corporation，Zhoushan 316021，China；2.Zhejiang Zhoushan Marine Power Research Institute Co.，Ltd.，Zhoushan 316021，China；3.Department of Electronic and Communication Engineering，North China Electric Power University，Baoding 071003，China）

Three-core optical fiber composite submarine power cable is often stretched in the process of transportation，installation and operation.The stress and strain data can be obtained by the simulation of stretching modeling，which is difficult to be obtained in the entity experiment.In this paper，the finite element model of the submarine cable is established，which is used to simulate the stretching process.We first simplify the structure of the submarine power cable according to the mechanical characteristics of stretching.The constraint，speed loads and hourglass control are then imposed to simulate the axial stretching process.Finally，the stress and strain data of the finite element model is extracted and processed to obtain the stress and strain of the submarine cable.The results provide reference for the study of mechanical properties and operating condition judgement of the submarine power cable.

three-core optical fiber composite submarine cable；stretch；stress；strain；finite element

TN818

A

1005-8788（2016）03-0029-04

10.13756/j.gtxyj.2016.03.010

2015-12-23

國家自然科學基金資助項目（51407074）；河北省自然科學基金資助項目（E2015502053）；中央高校基本科研業(yè)務費專項資金資助項目（2015ZD21）

盧志飛（1986-），男，浙江建德人。工程師，工學碩士，主要從事海底電纜運行維護技術方面的研究。