基于ANSYS的圓筒閥剛強度及屈曲分析與研究

2013-05-19 王燕 哈爾濱大電機研究所

  以某一電站圓筒閥作為計算算例,分別用有限元法和解析法對其進行剛強度及屈曲分析計算,并對如何選取單元類型,如何加邊界條件進行了探討分析,得出具有參考價值的結(jié)論。

1、前言

  圓筒閥位于水輪機固定導(dǎo)葉和活動導(dǎo)葉之間,可沿水輪機軸線方向上下移動。對于無法布置傳統(tǒng)進水閥的大型調(diào)峰調(diào)頻或多泥沙的水電站,采用圓筒閥結(jié)構(gòu)可以在停機過程中有效地保護導(dǎo)水機構(gòu)和減少電站的漏水損失。目前我國已有大潮山、漫灣、小灣、光照等水電站采用了圓筒閥結(jié)構(gòu)。性能和作用優(yōu)越的圓筒閥,在水電站設(shè)計選型中有著廣闊的發(fā)展前景。由于圓筒閥設(shè)置在頂蓋與座環(huán)上環(huán)形成的空腔內(nèi),工作時又要上下移動,所以保證圓筒閥的剛強度和穩(wěn)定性成為設(shè)計過程中的一個突出問題。

  本文以某電站圓筒閥作為算例,首先進行了筒閥剛強度計算分析,其次從對邊界條件選取、單元選擇、載荷施加方法等方面著手,提出了筒閥屈曲分析的有限元法和解析計算方法。

2、基本計算參數(shù)

  筒閥中心半徑:Rm=2575mm

  筒閥壁厚:δe=100mm

  筒閥高度:L=1310mm

  筒閥外徑:D0=5250mm

  最大靜水頭壓力:P1=0.726MPa

  升壓水頭壓力:P2=0.978MPa

  筒閥推拉桿個數(shù):6個

3、有限元法

3.1、圓筒閥剛強度分析

  進行圓筒閥剛強度分析,考核應(yīng)力值大小并找出局部最大應(yīng)力出現(xiàn)位置。

  利用圓筒閥的對稱性,切取包含一個推拉桿在內(nèi)的2π/Z(Z為推拉桿個數(shù))的扇形區(qū)域作為計算模型。本次分析采用大型結(jié)構(gòu)分析軟件ANSYS軟件完成,全部采用每個節(jié)點具有三個自由度的8節(jié)點六面體單元,即SOLID185單元,筒閥的應(yīng)力計算模型及網(wǎng)格剖分如圖1所示。

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

圖1 筒閥應(yīng)力計算模型及網(wǎng)格劃分圖

  邊界條件:為保證位移協(xié)調(diào)一致,在切開斷面處采用周期對稱邊界條件,同時在推拉桿的上端活塞處,約束下端面相應(yīng)節(jié)點的三個自由度。當圓筒閥全關(guān)閉時,承受0.726MPa的最大靜水壓力,當筒閥動水關(guān)閉時,承受0.978MPa外壓力。

  計算工況:我們主要考慮兩種工況,即最大靜水頭工況和升壓水頭工況。

  表1給出了兩種工況下的分析計算結(jié)果。圖2為升壓水頭工況下筒閥整體的應(yīng)力分布,圖3為升壓水頭工況下筒閥整體的變形分布。

表1 兩種工況下筒閥計算結(jié)果

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

圖2 升壓水頭工況下筒閥整體的應(yīng)力分布

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

圖3 升壓水頭工況下筒閥整體的變形分布

  從計算結(jié)果可看出,兩種工況下筒閥最大局部應(yīng)力分別為27.86MPa和20.68MPa,均出現(xiàn)在推拉桿與閥體接觸面上。

3.2、圓筒閥屈曲分析

  屈曲分析是一種用于確定結(jié)構(gòu)開始變得不穩(wěn)定時的臨界載荷和屈曲結(jié)構(gòu)發(fā)生屈曲響應(yīng)時的模態(tài)形狀的技術(shù)。圓柱形薄殼在均勻外壓作用下,壁內(nèi)應(yīng)力變?yōu)閴簯?yīng)力,則當外壓達到臨界值時,薄殼的圓形平衡就變?yōu)椴环(wěn)定,會突然變成長圓形。由屈曲導(dǎo)致結(jié)構(gòu)的不穩(wěn)定,對于結(jié)構(gòu)的設(shè)計來說是致命的,因為通常在結(jié)構(gòu)強度還遠沒達到極限時就發(fā)生了屈曲。所以對圓筒閥進行屈曲分析是非常必要的。

  屈曲問題的有限元分析方法大致有兩類:一類是通過特征值分析計算屈曲載荷,根據(jù)是否考慮非線性因素對屈曲載荷的影響。這類方法又細分為線性屈曲和非線性屈曲分析。另一類是利用結(jié)合Newton-Raphson迭代的弧長法來追蹤確定加載方向,追蹤失穩(wěn)路徑的增量非線性分析方法能有效的分析非線性屈曲和失穩(wěn)問題。本文只針對第一種方法中的線性屈曲分析進行討論,這種方法通過提取使線性系統(tǒng)剛度矩陣奇異的特征值來獲得結(jié)構(gòu)的臨界失穩(wěn)載荷及失穩(wěn)模態(tài)。

  進行筒閥的屈曲分析計算時,取整個閥體作為計算模型。對單元和邊界條件的選取分別有如下幾種方案。

  單元選。

  (1)分析時采用每個節(jié)點具有三個自由度的8節(jié)點六面體單元,即SOLID185單元。

  (2)采用每個節(jié)點具有三個自由度的20節(jié)點六面體單元,即SOLID186單元。

  通過多個算例對比分析得知,采用不同的單元類型,其計算結(jié)果也不相同,而且單元長度的大小對計算結(jié)果有一定的影響。一般情況下,單元網(wǎng)格剖分的越小,其計算結(jié)果越接近于GBl50)1998的結(jié)果,所以建議采用SOLID186單元。圖4為有限元分析模型。

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

圖4 筒閥屈曲分析計算模型

  邊界條件選取方法:

  (1)在接力器活塞桿處約束軸向位移,并約束閥體中心半徑處0°,90°,180°,270°四個節(jié)點的切向位移。

  (2)在圓筒閥的底部,約束中心半徑處0°,90°,180°,270°位置節(jié)點的切向位移和軸向位移。

  (3)在圓筒的外側(cè)施加1MPa(一個單位載荷)的壓力,不用施加實際載荷。因為由屈曲分析計算出的特征值,表示屈曲載荷系數(shù)。所以,若施加的是單位載荷,則該特征值就表示實際的屈曲載荷。

  經(jīng)過比較分析,建議今后對圓筒閥進行屈曲分析時采用第二種邊界條件約束方法,并按方法(3)施加載荷。

  通過有限元屈曲分析可知,圓筒閥屈曲載荷為3.495MPa,而實際圓筒閥承受的最大外壓為0.978MPa,臨界壓力安全系數(shù)為3.57。因此,該圓筒閥在外壓力作用下,不會發(fā)生失穩(wěn)情況,圖5為筒閥的有限元屈曲分析結(jié)果。

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

圖5 筒閥的屈曲分析結(jié)果

4、解析法計算圓筒閥平均應(yīng)力、變形及許用外壓

  周向平均應(yīng)力:

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

  徑向變形:

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

  許用外壓:

基于ANSYS的圓筒閥剛強度及屈曲分析與研究

  式中系數(shù)B值根據(jù)GB150-1998中的圖6-2、6-4中查取。

  采用解析法計算得到的是閥體的平均應(yīng)力,其應(yīng)力及變形水平與有限單元法計算結(jié)果接近。

5、結(jié)語

  (1)通過對某電站圓筒閥剛強度分析,發(fā)現(xiàn)筒閥最大應(yīng)力位置通常出現(xiàn)在推拉桿與閥體接觸面上。

  (2)邊界條件選取、單元選擇、載荷施加方法等對屈曲分析結(jié)果是有影響的。建議在今后進行筒閥屈曲分析時,參照本文提出的方法進行計算。

  (3)從有限元和解析計算結(jié)果可以看出本文提出的計算方法真實可靠,從而為今后對圓筒閥進行剛強度和屈曲分析提供了強有力的技術(shù)支持。