本申請是申請?zhí)枺?01410423136.3、申請日：2014.8.26、名稱“一種基于有限元分析的輸電鐵塔桿件應力計算方法”的分案申請。

本發(fā)明涉及一種對輸電鐵塔進行結構力學分析并精確計算每根桿件應力的方法，屬于輸電線路運行狀態(tài)安全評價領域。

背景技術：

電能生產(chǎn)與傳輸是國家經(jīng)濟社會發(fā)展的命脈，輸電線路更是關系國計民生的“生命線”。輸電線路作為重要的生命線工程，為電力能源的高效供給和合理分配提供了有力保證。經(jīng)濟社會的持續(xù)健康發(fā)展，對輸電線路運行的安全和可靠性提出了更高要求。然而，由于我國幅員遼闊，輸電線路路徑長、塔架高、大量穿越深山大河，甚至是無人區(qū)，自然地理氣象條件復雜多變，輸電線路作為大電網(wǎng)的主干，其自身是一個巨大的環(huán)境災害承載體。長期運行資料表明線路安全事故大多由機械力學原因造成。雖然在線路規(guī)劃設計中，已基于鐵塔塔材理論計算強度值，按特定工況條件計算允許應力并設置安全系數(shù)，在一定程度保證線路的安全運行。但實際運行工況復雜且時變，且無法完全由設計規(guī)范所預計，部件實際運行應力可能超出安全閾值；因此實際輸電線路雖已嚴格按規(guī)范設計、施工、運行，但斷線、倒塔等安全事故仍時有發(fā)生。

針對08年特大冰雪災害，業(yè)界和學界專家普遍認為電網(wǎng)缺乏科學有效的防災減災事故預警綜合體系等是停電事故發(fā)生的重要原因之一。建立這一體系的關鍵是對輸電線路的安全狀態(tài)給予科學評價。當前相關研究主要通過應用拉力傳感器，設計在線應力監(jiān)測裝置。由于硬件成本、裝置本體易受外力破壞等影響因素，其進一步推廣應用受到了制約。開發(fā)一種零硬件成本、無環(huán)境制約、可復用在所有運行輸電線路的結構安全評價系統(tǒng)成為了電網(wǎng)安全運行的當務之急，其中對鐵塔桿件進行精確的應力分布計算是其中的關鍵。

針對輸電鐵塔的應力分布分析計算，見諸公開報導的技術路線均是對鐵塔結構進行多次簡化等效，將鐵塔所受載荷進行簡單的分解或疊加，并忽略鐵塔結構的彈性形變。此種分析方法給出的結果，無法為安全評價提供準確地應力數(shù)據(jù)。從力學分析的角度考慮，不精確。故本發(fā)明提出，對鐵塔結構進行有限元建模，考慮到鐵塔結構的彈性形變，對每一根桿件進行力學分析，并根據(jù)鐵塔結構對單根桿件的力學性質進行合理疊加，給出整體鐵塔結構的精確力學性質，以得到鐵塔結構在復雜載荷下的精確應力分布。

技術實現(xiàn)要素：

為了解決上述問題，本發(fā)明提出一種基于有限元分析的輸電鐵塔桿件應力計算方法，用于分析計算鐵塔每根桿件的應力，包括如下步驟：

步驟1：根據(jù)有限單元法建立鐵塔結構的有限元力學模型，根據(jù)得到的力學模型獲取鐵塔節(jié)點；并將鐵塔結構離散化，得到與鐵塔兩節(jié)點間桿件數(shù)目相同的單元，單元彼此之間僅靠節(jié)點(鐵塔中主材與主材、主材與斜材、斜材與斜材的交匯點視為節(jié)點)相連；

步驟2：針對每個單元和節(jié)點位移，分別生成單元剛度矩陣[k]^(e)和節(jié)點位移陣列并根據(jù)桿件之間的空間角度關系、桿件之間的連接關系，轉換疊加出鐵塔整體剛度矩陣

步驟3：將鐵塔所受均布載荷和非節(jié)點載荷等效移置到節(jié)點上，形成節(jié)點載荷陣列(其中)，并考慮到矩陣方程中的奇異性，引入鐵塔塔腿四節(jié)點的位移約束條件，求解矩陣方程，得到節(jié)點位移陣列

步驟4：計算鐵塔節(jié)點的應變、應力，獲取鐵塔結構中薄弱的桿件。

其中，步驟1中的鐵塔結構離散化方法是將鐵塔中主材與主材、主材與斜材、斜材與斜材的交匯點歸結為節(jié)點，每兩個節(jié)點之間的桿件視為一個桿單元。

步驟2中的節(jié)點位移陣列和單元剛度矩陣[k]^(e)分別為：

1)點位移陣列的表達式為：

其中，為整體坐標系中的節(jié)點位移陣列；為第1個節(jié)點的位移陣列；為第2個節(jié)點的位移陣列；以此類推為第n個節(jié)點的位移陣列；ui，vi，wi為第i節(jié)點在局部坐標系中三個方向的線位移；θxi，θyi，θzi為第i節(jié)點處截面繞三個坐標軸的轉動，θxi代表截面的扭轉，θyi，θzi分別代表截面在xz及xy坐標面內(nèi)的轉動。

2)單元剛度矩陣[k]^(e)的表達式為：

其中，[k]^(e)為桿單元在單元局部坐標系內(nèi)的剛度矩陣；a為桿單元橫截面面積；iy為在xz面內(nèi)截面慣性矩；iz為在xy面內(nèi)的截面慣性矩；ip為單元的扭轉慣性矩；l為長度；e和g分別為材料的彈性模量和剪切模量。

根據(jù)所獲得的單元剛度矩陣[k]^(e)并依據(jù)桿件之間的空間角度關系、桿件之間的連接關系，即可得到鐵塔整體剛度矩陣

步驟3中的節(jié)點載荷矩陣為：

其中，為整體坐標中所有節(jié)點載荷陣列；為整體坐標中第i個節(jié)點的載荷列陣；nxi為第i個節(jié)點的軸向力，nyi、nzi分別為第i個節(jié)點在xy及xz面內(nèi)的剪力；mxi為第i個節(jié)點的扭矩，myi、mzi為第i個節(jié)點在xz及xy面內(nèi)的彎矩。

考慮到矩陣方程中為奇異矩陣，為使得方程組有解，本發(fā)明中引入鐵塔4個塔腳與基礎連接部分固定，限制鐵塔結構的剛性位移這一約束條件，從而保證整體剛度方程具有唯一解。

步驟4中計算鐵塔節(jié)點的應變和應力公式為：

其中，σx，σy，σz為坐標軸x，y，z方向上的3個正應力分量；τxy，τyz，τzx為在xy平面，yz平面，zx平面內(nèi)的3個切應力分量；分別為鐵塔空間結構在發(fā)生形變時，節(jié)點處在坐標軸x，y，z三個方向上的產(chǎn)生的線應變分量；為鐵塔空間結構在形變時，節(jié)點處在xy平面，yz平面，zx平面內(nèi)產(chǎn)生的3個剪應變分量。；u＝u(x,y,z)，v＝v(x,y,z)，w＝w(x,y,z)分別為鐵塔空間結構在形變時，節(jié)點在坐標軸x，y，z方向的位移；e為桿件的彈性模量；μ為泊松比。

本發(fā)明的技術效果：

1)充分考慮鐵塔結構的彈性形變，對每一根桿件進行力學分析，根據(jù)鐵塔結構對單根桿件的力學性質進行疊加，使得分析結果更加精確。

2)不僅可分析單一工況下鐵塔結構受力，對分析復雜工況下的鐵塔結構受力分析也具有優(yōu)越的性能。

3)計算中引入鐵塔4個塔腳與基礎連接部分固定，限制鐵塔結構的剛性位移這一約束條件，從而保證整體剛度方程具有唯一解，巧妙地解決了矩陣方程中的奇異性矩陣的奇異性問題。

附圖說明

圖1為220kv干字型鐵塔示意圖。

圖2為鐵塔模型中塔腳施加約束條件示意圖。

具體實施方式

一種基于有限元分析的輸電鐵塔桿件應力計算方法，可概括為四個階段：前期處理、有限元力學分析、工況荷載處理和后期處理。前期處理包括建立鐵塔結構的有限元力學模型和對整體結構進行離散化處理；有限元力學分析即對鐵塔結構的有限元力學模型進行分析，結合鐵塔結構中各桿件之間的空間角度關系和連接關系，疊加生成整體剛度矩陣；工況荷載處理即將鐵塔所受均布載荷、非節(jié)點載荷等效移置到節(jié)點上，形成節(jié)點載荷陣列；后期處理即以節(jié)點位移陣列作為未知量，與整體剛度矩陣，載荷陣列組成矩陣方程，并求解矩陣方程，給出節(jié)點應變，最終找出鐵塔結構中薄弱的桿件。該方法主要包括如下步驟：

步驟1：根據(jù)有限單元法建立鐵塔結構的有限元力學模型，根據(jù)得到的力學模型獲取鐵塔節(jié)點；并將鐵塔結構離散化，得到與鐵塔兩節(jié)點間桿件數(shù)目相同的單元，單元彼此之間僅靠節(jié)點(鐵塔中主材與主材、主材與斜材、斜材與斜材的交匯點視為節(jié)點)相連；

步驟2：針對每個單元和節(jié)點位移，分別生成單元剛度矩陣[k]^(e)和節(jié)點位移陣列并根據(jù)桿件之間的空間角度關系、桿件之間的連接關系，轉換疊加出鐵塔整體剛度矩陣

步驟3：將鐵塔所受均布載荷和非節(jié)點載荷等效移置到節(jié)點上，形成節(jié)點載荷陣列(其中)，并考慮到矩陣方程中的奇異性，引入鐵塔塔腿四節(jié)點的位移約束條件，求解矩陣方程，得到節(jié)點位移陣列

步驟4：計算鐵塔節(jié)點的應變、應力，獲取鐵塔結構中薄弱的桿件。

下面對每個步驟作進一步詳細說明：

步驟1中：根據(jù)有限單元法建立鐵塔結構的有限元力學模型，根據(jù)得到的力學模型獲取鐵塔節(jié)點；并將鐵塔結構離散化，得到與鐵塔兩節(jié)點間桿件數(shù)目相同的單元，單元彼此之間僅靠節(jié)點(鐵塔中主材與主材、主材與斜材、斜材與斜材的交匯點視為節(jié)點)相連，其實施過程為：

首先建立鐵塔結構的有限元力學模型，以鐵塔的橫擔方向作為整體坐標系的x軸，線路方向作為y軸，豎直方向作為z軸，并滿足右手定則；以桿單元所在直線作為單元局部坐標系的x軸，桿件與局部坐標系下的x軸方向重合，其正方向與整體坐標系x軸正方向一致。

輸電鐵塔作為空間桿件系統(tǒng)，用桿-梁單元(以后簡稱桿單元)進行離散，本發(fā)明中將鐵塔中主材與主材、主材與斜材、斜材與斜材的交匯點化為節(jié)點，每兩個節(jié)點之間的桿件視為一個桿單元。由于鐵塔的結構比較復雜，桿件數(shù)目比較多，若在結構離散化即劃分單元時，單元分的越小，單元數(shù)目就越多，計算時間就越長，因此考慮按照鐵塔本身的自然結構劃分，既提高了計算精度，又減少了計算工作量。

步驟2中：針對每個單元和節(jié)點位移，分別生成單元剛度矩陣[k]^(e)和節(jié)點位移陣列并根據(jù)桿件之間的空間角度關系、桿件之間的連接關系，轉換疊加出鐵塔整體剛度矩陣

對鐵塔結構離散化后，要對單元進行力學特性分析，即確定單元節(jié)點力和節(jié)點位移之間的關系。為了分析和確定這一關系，需要選擇位移模式，位移函數(shù)是單元上點的位移對點的坐標的函數(shù)，本方法用單元內(nèi)部點的坐標的多項式來表示，空間中的桿件，每個節(jié)點具有6個自由度，即桿件除了承受一維軸力、兩維剪力和兩維彎矩的作用外，還可能承受一維扭矩的作用。并且，空間桿單元承受一維軸力、兩維剪力、兩維彎矩、一維扭矩，即對應著節(jié)點的6個自由度，分別為3個方向上的線位移和在節(jié)點處截面繞3個坐標軸的轉動，因此單元內(nèi)部點的坐標的多項式可表示為δ＝k1u+k2v+k3w+k4θx+k5θy+k6θz，據(jù)此，可形成所有節(jié)點的位移陣列

其中，為整體坐標系中的節(jié)點位移陣列；為第1個節(jié)點的位移陣列；為第2個節(jié)點的位移陣列；以此類推為第n個節(jié)點的位移陣列；ui，vi，wi為第i節(jié)點在局部坐標系中三個方向的線位移；θxi，θyi，θzi為第i節(jié)點處截面繞三個坐標軸的轉動，θxi代表截面的扭轉，θyi，θzi分別代表截面在xz及xy坐標面內(nèi)的轉動。

建立單元剛度方程的基本步驟為：在假定單元位移函數(shù)的基礎上，根據(jù)彈性力學理論，來建立應變、應力與節(jié)點位移之間的關系式。然后根據(jù)虛位移原理，求得單元節(jié)點力與節(jié)點位移之間的關系，從而得出如下單元剛度矩陣[k]^(e)：

其中，[k]^(e)為桿單元在單元局部坐標系內(nèi)的剛度矩陣；a為桿單元橫截面面積；iy為在xz面內(nèi)截面慣性矩；iz為在xy面內(nèi)的截面慣性矩；ip為單元的扭轉慣性矩；l為長度；e和g分別為材料的彈性模量和剪切模量。

單元單元剛度矩陣[k]^(e)得到后，根據(jù)桿件之間的空間角度關系、桿件之間的連接關系，轉換疊加出鐵塔整體剛度矩陣其具體實施過程如下：

首先，假設局部坐標為x，y，z；整體坐標為以從局部坐標系的x軸正方向順時針轉到整體坐標系的軸正方向為正，則x軸的方向余弦為：

y軸的方向余弦為

z軸的方向余弦為

令

則整體坐標系內(nèi)的單元剛度矩陣表示為：

其中，[λ0]為節(jié)點轉換陣，[λ]為坐標轉換矩陣，[k]^(e)為桿單元在單元局部坐標系內(nèi)的剛度矩陣。

整體剛度矩陣的集成規(guī)則如下：

1)先求出每個單元的剛度矩陣[k]^(e)；

2)將的每個子塊進行換號，換成對應的整體編號；

3)將換成以后的字塊送到整體剛度矩陣中的對應位置上；

4)若在同一位置上有幾個單元的相應子塊送到時，則應進行疊加。

經(jīng)過上述步驟后，就可以得到整體剛度矩陣中的每個子塊，從而形成了整體剛度矩陣

其中，為第j個節(jié)點產(chǎn)生單位位移時，在第i個節(jié)點上引起的節(jié)點力，稱為剛度子矩陣。值得注意的是：中每個子塊都是6×6階矩陣，如果整體結構具有n個節(jié)點，那么整體剛度矩陣的階數(shù)6n×6n。

步驟3將鐵塔所受均布載荷和非節(jié)點載荷等效移置到節(jié)點上，形成節(jié)點載荷陣列(其中)，并考慮到矩陣方程中的奇異性，引入鐵塔塔腿四節(jié)點的位移約束條件，求解矩陣方程，得到節(jié)點位移陣列其具體實施過程為：

首先將鐵塔所受均布載荷、非節(jié)點載荷等效移置到節(jié)點上，形成節(jié)點載荷陣列空間中的桿件，每個節(jié)點具有6個自由度，即桿件除了承受軸力、剪力和彎矩的作用外，還可能承受扭矩的作用。并且，空間桿單元承受一維軸力、兩維剪力、兩維彎矩、一維扭矩，即對應著節(jié)點的6個自由度。輸電鐵塔的桿單元正是空間桿單元。

其中，為整體坐標中所有節(jié)點載荷陣列；為整體坐標中第i個節(jié)點的載荷列陣；nxi為第i個節(jié)點的軸向力，nyi、nzi分別為第i個節(jié)點在xy及xz面內(nèi)的剪力；mxi為第i個節(jié)點的扭矩，myi、mzi為第i個節(jié)點在xz及xy面內(nèi)的彎矩。

由于矩陣方程中為奇異矩陣，方程組無解，若要求解該方程，必須引入約束條件，限制鐵塔結構的剛性位移，保證整體剛度方程有唯一解。位移約束條件的作用是使結構上的節(jié)點的位移分量為常數(shù)值，即δi＝δ0。引入位移約束條件，就是要將δi＝δ0引入到結構總體剛度方程中。

本發(fā)明中，采用對角元素置1法，將δi＝δ0引入整體剛度矩陣

針對輸電鐵塔的4個塔腿中，與基礎連接的部分是固定端約束，因此δ0＝0；將k的第i行的主對角線元素kii置1，其余元素清零，且將第i行的載荷項ri用δ0代替，上式變?yōu)?/p>

通過置1法將位移約束條件δi＝δ0引入到整體剛度方程中，并沒有改變矩陣k和r中的各元素儲存順序，而且矩陣k仍然為對稱矩陣。

本發(fā)明計算方法中，對與基礎相連的4個塔腿節(jié)點的子矩陣采用置1法，即可代入24個位移邊界條件，消除整體剛度矩陣的奇異性，從而采用高斯消元法進行矩陣方程求解。

步驟4：計算節(jié)點應變、應力等，找出鐵塔結構中最薄弱的桿件，其實施過程如下：

由鐵塔結構的整體剛度矩陣方程求解出各節(jié)點的位移δ后，就可以得到各單元的節(jié)點位移δ^e。鐵塔結構受力后，其內(nèi)部各點將沿x，y，z三個坐標軸方向發(fā)生位移。如果各點沿x，y，z三個坐標軸方向的位移以u，v，w表示，它們是點的坐標函數(shù)，即u＝u(x,y,z)，v＝v(x,y,z)，w＝w(x,y,z)。

鐵塔結構在形變時，內(nèi)部任意一點處有3個線應變分量εx，εy，εz以及3對剪應變分量γxy＝γyx，γyz＝γzy，γzx＝γxz。由彈性力學可知，應變與位移之間的關系即幾何方程為：

式中εx，εy，εz為分別為鐵塔空間結構在發(fā)生形變時，節(jié)點處在坐標軸x，y，z三個方向上的產(chǎn)生的應變分量；γxy，γyz，γzx為鐵塔空間結構在形變時，節(jié)點處在xy平面，yz平面，zx平面內(nèi)產(chǎn)生的的3個應變分量。

鐵塔結構受到作用時，內(nèi)部任意一點處的應力狀態(tài)也是三維的，有3個正應力分量σx，σy，σz以及三對切應力分量τxy＝τyx，τyz＝τzy，τzx＝τxz。

在線彈性范圍內(nèi)，應力和應變用如下方程表示：

其中，σx，σy，σz為坐標軸x，y，z方向上的3個正應力分量；τxy，τyz，τzx為在xy平面，yz平面，zx平面內(nèi)的3個切應力分量；分別為鐵塔空間結構在發(fā)生形變時，節(jié)點處在坐標軸x，y，z三個方向上的產(chǎn)生的線應變分量；為鐵塔空間結構在形變時，節(jié)點處在xy平面，yz平面，zx平面3個平面內(nèi)產(chǎn)生的剪應變分量；e為桿件的彈性模量；μ為泊松比。在本說明書中未作詳細描述的內(nèi)容屬本領域技術人員的公知技術。