大型水輪機球閥有限元分析邊界條件探討

1、大型水輪機球閥有限元分析邊界條件探討*呂桂萍1趙冰燃2(1. 水力發(fā)電設備國家重點實驗室哈爾濱 150040;2. 廈門大學 廈門 361005摘要:本文通過對球閥不同邊界條件的有限元對比分析計算,給出了不同邊界條件下的分析計算結果,同時,根據球閥的具體結構形式和受力特點,最終給出了球閥有限元分析計算邊界條件,對今后球閥有限元分析具有重要的指導意義。關鍵詞:球閥邊界條件力學模型中圖分類號:TG156Research on Boundary Condition to Use FEA for Sphere ValveLU Guiping1 ZHAO Bingran2(1.State Key Lab

2、oratory of Hydropower Equipment, Harbin 150040;2. Xiamen University , Xiamen 361005 Abstract:The results are given by different boundary condition and different mechanics model in this page. According to the structure characteristic of sphere valve, the last mechanics model and boundary condition to

3、 analysis sphere valve are given in this page. It is very important to analysis stress level for sphere valve.Key words:Sphere Valve Boundary Condition Mechanics Mode0 前言球閥是高水頭電站以及抽水蓄能機組必備的成套產品,也是水輪發(fā)電機組的關鍵部件之一。近年來,隨著抽水蓄能機組在國內的迅速建立,對球閥的有限元法分析,掌握球閥在各種工況下的應力和位移水平顯得十分必要。眾所周知,有限元分析結果的可信度是建立在合理的力學模型和邊界條件基

4、礎上的,因此,合理地建立球閥有限元分析的力學模型和邊界條件,是掌握球閥各個工況下的應力和位移水平的關鍵。本文以某電站的球閥結構為依據,采用目前通常采用的幾種球閥有限元分析的力學模型和邊界條件進行對比分析,并給出分析結果。從結果中找出各種文件力學模型以及邊界條件的優(yōu)劣,最終根據球閥的工作狀態(tài),給出比較合理的球閥有限元分析力學模型和邊界條件。國家重點基礎研究計劃(973計劃(2010C B736200-81 有限元基本方程首先應引入形函數的概念。形函數是定義在單元內部的坐標的連續(xù)函數,它應滿足下述條件: 對應于節(jié)點i的形函數有 Ni=1;Nj=0 (ij;能保證用它定義的位移函數在單元間邊界上連續(xù)

5、和單值;形函數應包含足夠的線性項;對于一個單元,形函數應滿足等式:1=i N這些條件都是保證由此形函數導出的單元,能夠收斂于正確解的必要條件。三維母單元是(,坐標系中的2×2×2的正六面體(如圖1所示。單元邊界是六個平面 =±1,=±1,=±1那么在局部坐標系,下,其單元的形函數表達式為:1(1(1(81iiiiN+=459其中i i i 為第i 個節(jié)點的局部坐標值,i=1,2,8(1 圖1 三維母單元 3圖2 三維子單元作為空間結構中的三維子單元,可以直接以結構分析所用的坐標系作為全局坐標系,如圖2中的xyz 坐標。如果已知子單元各節(jié)點的全局

6、坐標xi、yi、z i,那么子單元內某點(,的全局坐標可由下式確定:=ii i i i i z N z y N y x N x ,(,(,(,(,(,(2位移函數可表示為:=mi ii mi ii mi ii w N w v N v u N u 111,(,(,(,(,(,(3或寫成矩陣形式:=m m m m m m w v u w v u N N N N N N N N N w v u 11121212100000000000000 (4記作N f =(5為便于分析,可將(2.21式寫成分塊形式:=mi ii N f 1(6式中m 為單元節(jié)點數;Ti iii w v u = (7表示節(jié)點位移

7、向量;=i i ii N N N N 00000 (8必須注意,坐標變換式和位移函數的描述都采用了形函數,但它們的含義是不同的。在坐標變換時,是利用形函數以節(jié)點坐標來內插單元點坐標,屬于幾何圖形的變換;而在描述位移時,則是通過形函數以節(jié)點位移來內插單元的位移場,它描述的是單元的變形形態(tài)。在小位移彈性理論中,對于三維問題,應變-位移關系可表示成:=w v u z xz y y x z y x zx yz xy z y x 00000000 (9記作f L = (10將(6 式代入得:460imi i m i i B N L =11(11其中=x N zN y N z N x N y Nz N y

8、 N x N B i i i i i i i i i i 000000000 (12上式出現了形函數對全局坐標系的微商,但形函數是由計算坐標系給出的。運用復合函數的微商法則+=+=+=z z N y yN x xN N z z N y y N x x N N zz N y y N x x N N i iii i iii i iii寫成矩陣形式為:=z N y N x N J z N y N xN z y x z y x z yxN N N i i i i i i i i i (13 J為Jacobin 矩陣。=m mm m m m z y xz y x z y x N N N N N N N

9、N N J (222111212121 (14記作:xyz DNK J =(13 式可得:=i i i i i i N N N J z N y N x N 1 (15應力可表示成:=mi i i B D D 1 (16+=G G G G G G D 000000000000222 (17zx yz xy z y x = (181(2;21(1(+=+=EG E (19以(11和(16式代入應變能表達式:dv u vT=21令TT m T T 21= (20m B B B B 21= (21則有k u T21=(22dvB D B k vT = (23k即為三維等參元的剛度矩陣,可表示為:461

10、=mm m m m m k k k k k k k k k k212222111211 (24其中子矩陣:dvB D B k vs T r rs = (25r,s=1,2,3,m由于Bi 是計算坐標的函數,因此上式應化成對計算坐標系進行積分,即:=d d d J B D B K s Tvr rs (26r,s=1,2,3,m應用卡氏第一定理,由(2.38式可得等參單元的節(jié)點平衡方程:F K = (27式中F 是與相對應的節(jié)點力向量。2 分析算例2.1 球閥基本參數公稱直徑: 3500.0 m 工作壓力: 2.619 MPa 試驗壓力: 4.84 MPa 水頭: 262.0 m 流量: 175.

11、0 m3s-1 2.2 力學模型及邊界條件球閥位于上下游連接管之間,通過游離法蘭將上下游連接管和球閥連接。通常情況下,下游連接管又稱為伸縮節(jié)。對球閥進行剛強度分析時,主要考慮一下3種工況:正常運行工況;工作密封工況,此時球閥在下游側關閉; 檢修密封工況,此時球閥在下游側關閉; 選取上下游連接管游離法蘭以及球閥作為分析計算模型,如圖3所示。在上述分析模型中,分別考慮一下3中邊界條件:上、下游連接管與球閥閥體之間采用一體的結構形式; 圖3. 球閥有限元分析計算模型上、下游連接管與球閥閥體之間不是一體,即法蘭接觸面試兩個面,但是在螺栓連接處,螺栓直徑大小的圓上是一個面,兩個法蘭之間呈現自由狀態(tài);上、

12、下游連接管與球閥閥體之間不是一體,即法蘭接觸面試兩個面,但是在螺栓連接處,螺栓直徑大小的圓上是一個面,兩個法蘭之間采用接觸單元。在分析時,為防止產生剛體位移,約束上游連接管與混凝土連接斷面上的所有節(jié)點的X;Y;Z 自由度,同時約束地腳把合平面相應節(jié)點垂直于地面方向的自由度。載荷:正常工況下,上下游連接管以及球閥內腔承受水壓力;工作密封工況,由于活門關閉,作用在活門上的水壓力通過樞軸傳遞給球閥的樞軸孔處,此時上游連接管和球閥承受水壓力,而下游連接管沒有水壓力;檢修密封工況,由于活門關閉,作用在活門上的水壓力通過樞軸傳遞給球閥的樞軸孔處,此時上游連接管承受水壓力,而球閥和下游連接管沒有水壓力。分別

13、選取以下兩種分析模型對球閥進行分析。 2.3 有限元分析計算結果本次分析采用ANSYS 軟件完成。選取每個節(jié)點具有3個自由度的10節(jié)點四面體單元,即SOLID92單元。在考慮接觸問題的力學模型中,分別采用131293個節(jié)點和81785個單元進行模擬。表1給出了分析計算結果。從表1不難看出,在正常工況下,3種邊界條件得出的應力以及位移基本相當,但在其他兩種工況下,邊界條件1和邊界條件2及3相差較大;而邊界條件2和3相差較小。但是,在邊界條件2中,由于沒有考慮接觸單元,可能存在兩個法蘭面之間出現咬合的情況,這與實際結構是不吻合的。4623 結論本文通過不同的邊界條件進行對比計算發(fā)現:在相同的力學模

14、型下,不同的邊界條件對球閥的位移和應力結果有一定的影響。考慮到工程實際情況,建議今后對球閥進行有限元分析時采用本文提到的力學模型和第3種邊界條件。參考文獻1哈爾濱大電機研究所編制.水輪機設計手冊.北京:機械工業(yè)出版社.1976年.2R.Pavillet 16 Spherical valves for three pumpstorage power stations in China. Alstom Power Hydro.3 Ph. BAUD SPHERICAL V ALVE PLUG AND MOBILERING CALCULATION. Alstom Power Hydro.4 Ph.BA

15、UD SPHERICAL V ALVE BODY CALCULATION.Alstom Power Hydro.5 D BARDINET. SPHERICAL V ALVE FATIGUEANALYSIS TECHNICAL NOTICE. Alstom PowerHydro.作者簡介:呂桂萍(通信作者,女,1963年出生,博士,長期從事水輪發(fā)電機組大部件剛強度研究工作。E-mail:pedestrian2001附表1. 分析計算結果正常工況工作密封檢修密封邊界條件位移方向位移/mm應力/MPa最大應力位置位移/mm應力/MPa最大應力位置位移/mm應力/MPa最大應力位置R -0.571.78 -4.033.95-4.904.701Z -1.430 145.21上游連接管處-0.