劉銳鵬,張耿斌,顏天佑
(1.廣州電力設計院,廣州 510610;2.廣州供電局有限公司輸電管理所,廣州 510620)
基于LS-DYNA的架空輸電線找形程序
劉銳鵬1,張耿斌2,顏天佑1
(1.廣州電力設計院,廣州 510610;2.廣州供電局有限公司輸電管理所,廣州 510620)
基于懸索結構初始構形分析的原理和方法,應用于架空輸電線的初始找形研究,以Visual Basic 6.0為平臺,介紹了不等高懸點、連續檔架空輸電線的找形程序,生成適用于LS-DYNA分析的APDL文件,提高了架空輸電線LS-DYNA找形工作的效率。
架空輸電線;找形;LS-DYNA;APDL
LS-DYNA有限元軟件是目前結構分析中使用最為廣泛的通用程序之一,其顯示算法特別適合于分析各種非線性結構沖擊動力學問題,如爆炸、結構碰撞、金屬加工成形、架空輸電線覆冰斷線等高度非線性的問題,同時還可以求解傳熱、流體以及流固耦合問題。因此研究架空輸電線找形在LS-DYNA上的實現具有一定的意義。由于LS-DYNA的全英文界面和相對較弱的前處理功能,在建模上需要花費較多的時間和精力。基于此,文章利用Visual Basic 6.0開發出架空輸電線找形程序,將生成的APDL文件導入LS-DYNA,便可實現對連續檔、不等高懸點的架空輸電線精確找形。下面以某大跨越線路為例,驗證了找形程序的可靠性和實用性。
1.1 位移函數[1]
采用Lagrangian方法,建立曲線坐標Oξ,令s為單元變形前任意點到中點的弧長,則ξ= 2s/l,l為單元變形前長度,即原長。架空輸電線單元模型如圖1所示。

圖1 架空輸電線單元模型
設位移

式中:

Ni(i=1,2,3)為形函數,由節點的位移連續條件求得:

同時得到坐標的變換式:

式中,X= [x y z ]T;Xe=

式 (2)便決定了單元的初始形狀。
1.2 非線性有限元法的基本原理
基本原理是在架空輸電線弦線位置創建幾何模型,采用適當的初應變和較小的彈性模量,施加自重荷載,逐步更新有限元模型,設定大變形以及應力剛化,以架空輸電線水平應力為收斂條件,當結果滿足收斂條件時退出循環迭代求解。在此線形下,恢復材料的真實彈性模量,設定較小的初應變,進行自平衡迭代求解,即求得架空輸電線在自重荷載下的初始形態[2]。
利用Visual Basic 6.0開發的架空輸電線找形程序,輸入幾何、材料、有限元、覆冰等基本參數,該程序便能自動生成適用于LS-DYNA分析的APDL命令流,實現架空輸電線找形的參數化建模。
架空輸電線找形程序流程圖如圖2所示,PRC為計算所需精度,MAX=ABS((F0-F1)/ F),其中F0為輸電線弧垂最低點處水平張力,F1為LS-DYNA迭代分析所得輸電線弧垂最低點處水平張力,ABS()為 Visual Basic 6.0中的函數。
算例:某大跨越的主跨檔距1 650 m,跨越塔采用酒杯型塔,塔高181.8 m,兩塔高差0.5 m,架空輸電線特征參數見表1。覆冰密度 0.8 g/cm3[3]。

圖2 架空輸電線找形程序流程圖

表1 架空輸電線材料特性
均勻覆冰厚度取15 mm,不均勻覆冰時將整檔架空輸電線按覆冰厚度分為5個區段,每區段水平投影長度330 m,假設同一區段內輸電線覆冰情況相同[4]。不均勻覆冰分布工況見表2。

表2 不均勻覆冰分布工況
覆冰荷載下架空輸電線的找形結果如圖3、圖4、圖5、圖6所示 (圖中虛線為架空輸電線自重荷載下初始構形,實線為覆冰后形變)。

圖3 均勻覆冰下架空輸電線形變

圖4 不均勻覆冰工況A下 架空輸電線形變

圖5 不均勻覆冰工況B下 架空輸電線形變

圖6 不均勻覆冰工況C下架空輸電線形變
1)提取均勻覆冰下架空輸電線檔距中央弧垂為116.065 m,均勻覆冰下解析法得到的弧垂理論值為115.885 m。二者找形結果差別不大,在0.2%以內,說明文中提出的找形程序是可靠的。
2)對于不均勻覆冰情況,只需要改變覆冰厚度參數,找形程序即可生成相應的APDL文件,有利于模型參數的反復修改,極大地提高了架空輸電線LS-DYNA建模、找形工作的效率。
[1] 唐建民,何署廷,段春平.懸索結構非線性有限元分析[J].河海大學學報,1998,26(6):45-46.
[2] LS-DYNA工程結構數值分析 [M].北京:人民交通出版社,2007:467-471.
[3] 李慶峰,范崢,吳穹,等.全國輸電線路覆冰情況調研及事故分析 [J].電網技術,2008,32(9):33-34.
[4] 孟遂民,單魯平.輸電線動力學分析中的找形研究 [J].電網與清潔能源,2009,25(10):43-44.
Form-finding Program Development for Overhead Transmission Line Based on LS-DYNA
LIU Ruipeng1,ZHANG Gengbin2,YAN Tianyou1
(1.Guangzhou Electric Power Design Institute,Guangzhou 510610;2.Power Transmission Management Office of Guangzhou Power Supply Bureau Corporation,Guangzhou,510620)
Based on the principles and methods of form-finding for cable structures are used in the initial form-finding of overhead transmission line.With the Visual Basic 6.0 as a platform,form-finding program for overhead transmission line was developed,which can complete form-finding for overhead transmission line of differ-height hang point and continuous spans.The APDL applied to LS-DYNA is automatically generated and improve the form-finding for overhead transmission line based on LS-DYNA efficiently.
overhead transmission line;form-finding;LS-DYNA;APDL
TM75
B
1006-7345(2014)05-0050-03
2013-10-17
劉銳鵬 (1982),男,碩士,工程師,廣州供電局有限公司設計院,主要從事輸電線路設計與運行維護方面工作 (e-mail) 30052535@qq.com。