精品久久久久中文字,日韩免费看视频三区中文字幕,国产成 人亚洲91,亚洲无玛在线观看

調(diào)節(jié)閥位于主蒸汽管道和汽輪機(jī)之間，是控制機(jī)組功率的重要通流部件。隨著機(jī)組單機(jī)容量增加和蒸汽參數(shù)的提高，汽輪機(jī)調(diào)節(jié)閥常處于極其嚴(yán)酷的工況環(huán)境中，由流體誘發(fā)調(diào)節(jié)閥的振動(dòng)所產(chǎn)生的事故屢見不鮮，已成為威脅電廠安全運(yùn)行的重要因素。隨著CFD的不斷發(fā)展和完善，數(shù)值計(jì)算因其方便性和可視化等優(yōu)點(diǎn)受到了眾多學(xué)者青睞，在分析流體誘發(fā)的調(diào)節(jié)閥振動(dòng)中扮演越來越重要的角色。但由于調(diào)節(jié)閥振動(dòng)是一種復(fù)雜的三維流固耦合現(xiàn)象，到目前為止對(duì)于該問題的機(jī)理始終沒有明確的論斷。
Araki認(rèn)為當(dāng)調(diào)節(jié)閥中閥桿系統(tǒng)的振動(dòng)在小開度時(shí)為一階縱向高頻振動(dòng)，在開度逐漸增大后主要為一階彎曲振動(dòng)。Zhang通過二維定常數(shù)值模擬指出閥內(nèi)不對(duì)稱流是引起振動(dòng)的主因。Morita的三維非定常數(shù)值模擬進(jìn)一步表明閥碟下方高壓區(qū)的周向移動(dòng)導(dǎo)致了閥桿系統(tǒng)的振動(dòng)。Morita還通過實(shí)驗(yàn)研究了閥桿在弱連接下的壓力波動(dòng)和振動(dòng)的關(guān)系，指出在中等升程時(shí)壓力波動(dòng)和閥固有頻率之間存在鎖定現(xiàn)象，在鎖定區(qū)域阻尼比接近于零且振幅較大，而在非鎖定區(qū)域阻尼比大于5%。Tecza通過有限元法確定了閥桿系統(tǒng)的固有頻率，并指出閥碟最有可能的運(yùn)動(dòng)為傾斜和擺動(dòng)，并按照閥桿擺動(dòng)頻率550Hz**擺動(dòng)角度為±0.02°進(jìn)行了瞬態(tài)力響應(yīng)分析。Yonezawa的研究表明閥碟在靈活狀態(tài)下的壓力波動(dòng)高于閥碟緊固時(shí)的波動(dòng)，將閥桿系統(tǒng)的振動(dòng)分為振動(dòng)幅值較小時(shí)的強(qiáng)迫振動(dòng)和振動(dòng)幅值較大時(shí)的自激振動(dòng)。國(guó)內(nèi)關(guān)于調(diào)節(jié)閥的數(shù)值計(jì)算研究往往不考慮閥桿系統(tǒng)振動(dòng)的影響，用數(shù)值模擬分析調(diào)節(jié)閥的閥型對(duì)振動(dòng)、流量、壓損、流型等因素的影響，而在調(diào)節(jié)閥快關(guān)和開啟特性中不得不考慮閥碟移動(dòng)，利用動(dòng)網(wǎng)格計(jì)算閥內(nèi)的瞬時(shí)流場(chǎng)的變化。但是在計(jì)算過程中由于網(wǎng)格畸變過大而導(dǎo)致計(jì)算發(fā)散一直是困擾數(shù)值計(jì)算能否順利進(jìn)行的關(guān)鍵因素。通過增加網(wǎng)格尺度減小網(wǎng)格畸變往往收效甚微，甚至有時(shí)還會(huì)因網(wǎng)格尺度過大而影響計(jì)算準(zhǔn)確性。
從這些研究可以看出，以往借助CFD分析閥內(nèi)的流型及流體誘發(fā)的調(diào)節(jié)閥閥桿系統(tǒng)振動(dòng)大都忽略振動(dòng)對(duì)閥內(nèi)流場(chǎng)的影響。但當(dāng)振動(dòng)幅值較大時(shí)，必然會(huì)對(duì)內(nèi)部流場(chǎng)的產(chǎn)生影響。本文針對(duì)以上問題，利用Fluent中的彈簧光順模型和UDF函數(shù)，對(duì)比分析了某600MW�；�調(diào)節(jié)閥的閥碟頂端振幅分別為0.6mm、1mm，及振動(dòng)頻率分別為46.4、92.8、185.6Hz時(shí)閥桿系統(tǒng)振動(dòng)對(duì)閥內(nèi)流場(chǎng)的影響。計(jì)算中提出的保護(hù)邊界層的網(wǎng)格變化方法和數(shù)值計(jì)算結(jié)果為后續(xù)數(shù)值計(jì)算中網(wǎng)格控制和是否考慮閥桿系統(tǒng)振動(dòng)對(duì)閥內(nèi)流場(chǎng)的影響提供了重要的參考。
1.數(shù)值計(jì)算的物理模型
根據(jù)相似原理對(duì)某電廠600MW調(diào)節(jié)閥進(jìn)行了�；瘜�(shí)驗(yàn)，圖1為�；�后調(diào)節(jié)閥沿進(jìn)口方向中心截面示意圖。閥桿和閥碟所組成的系統(tǒng)可看做為以提升機(jī)構(gòu)與閥桿連接處為固定點(diǎn)，頂端帶有集中質(zhì)量的長(zhǎng)細(xì)粱結(jié)構(gòu)[3]。由于調(diào)節(jié)閥在中等升程下，閥碟上壓力波動(dòng)呈現(xiàn)以低頻占優(yōu)的寬頻激勵(lì)特性，所以主要考慮閥桿閥碟所組成系統(tǒng)的一階彎曲模態(tài)對(duì)閥內(nèi)流場(chǎng)的影響。利用Workbench計(jì)算此調(diào)節(jié)閥閥桿系統(tǒng)固有頻率，結(jié)果表明此調(diào)節(jié)閥閥桿系統(tǒng)一階彎矩的固有頻率為46.4Hz。為對(duì)比閥桿系統(tǒng)按照一階彎矩振動(dòng)在不同振動(dòng)頻率下對(duì)閥內(nèi)流場(chǎng)的影響，計(jì)算了頻率分別為一階彎矩頻率及2倍、4倍一階彎矩頻率（即頻率分別為：46.4、92.8和185.6Hz）這3種方案。為對(duì)比不同幅值對(duì)閥內(nèi)流場(chǎng)的影響，計(jì)算了閥碟固定、閥碟頂端振動(dòng)幅值為0.6和1mm這3種方案。以上計(jì)算中，閥桿系統(tǒng)振動(dòng)方向?yàn)檠貓D1中進(jìn)口中心截面。為后續(xù)分析方便，定義三維坐標(biāo)所指方向如下：x軸正方向所指的沿來流方向?yàn)闁|側(cè)用E表示，x軸負(fù)方向所指為西側(cè)用W表示；y軸所指方向沿閥桿軸向；z軸所指方向垂直于沿來流中心截面，如圖1、2所示。


圖1 調(diào)節(jié)閥結(jié)構(gòu)圖 圖2 閥碟偏向西側(cè)示意圖（h=0.6mm）
圖2為閥碟頂端**振幅h=0.6mm時(shí)的閥碟偏向東側(cè)示意圖，可以看出閥碟可能出現(xiàn)的**振幅主要由閥碟和閥碟套間隙0.2mm決定。隨著使用時(shí)間增加，閥碟表面的涂層不斷脫落此間隙可達(dá)到0.4mm，此時(shí)閥碟頂端振幅為h=1mm。由于卸載室內(nèi)流過的流量占主流流量小于1%，對(duì)主流的影響可忽略不計(jì)，因此在動(dòng)網(wǎng)格計(jì)算中不考慮由卸載室流過的流體。
2.數(shù)值計(jì)算方法和動(dòng)網(wǎng)格設(shè)置
本文湍流模型采用Fluent中基于realizablek-ε的兩方程的DES模型。DES方法結(jié)合了雷諾平均與大渦模擬技術(shù)的優(yōu)點(diǎn)，在近壁面的附面層內(nèi)采用RANS方法，用realizablek-ε兩方程模型模擬其中的小尺度脈動(dòng)運(yùn)動(dòng)；在遠(yuǎn)離物面區(qū)域?qū)⑼牧髂Ｐ秃纳㈨?xiàng)中的湍流尺度參數(shù)用網(wǎng)格尺度與一常數(shù)的乘積代替，使其起Smagorinski大渦模擬的亞格子雷諾應(yīng)力模型的作用。這樣既能在附面層內(nèi)發(fā)揮前者計(jì)算量小的優(yōu)點(diǎn)，又可以在遠(yuǎn)離物面的區(qū)域?qū)Υ蟪叨确蛛x湍流流動(dòng)進(jìn)行較好的模擬。
選取實(shí)際運(yùn)行和�；瘜�(shí)驗(yàn)中閥桿系統(tǒng)振動(dòng)較大的工況，即相對(duì)升程（調(diào)節(jié)閥的**升程與配合直徑比值）εL=11%，壓比（調(diào)節(jié)閥出口靜壓與進(jìn)口總壓的比值）ψP=0.8作為計(jì)算工況。進(jìn)口設(shè)為總壓P0=123.5kPa，溫度為298K。出口設(shè)為靜壓，其值P2=98.8kPa，溫度為295K，所選取數(shù)據(jù)均為�；瘜�(shí)驗(yàn)所測(cè)。
結(jié)合Fluent中UDF函數(shù)和彈簧光順模型可實(shí)現(xiàn)閥碟在東西截面上以不同振幅和頻率按正弦振動(dòng)。UDF函數(shù)中主要通過式（1）控制振動(dòng)。
cg_omega=Acos（BTime） （1）
式中：cg_omega為在xy平面內(nèi)轉(zhuǎn)動(dòng)角速度，參數(shù)B控制頻率，A控制振動(dòng)幅值。Fluent中彈簧光順法把網(wǎng)格的邊理想化為節(jié)點(diǎn)間相互連接的彈簧，網(wǎng)格移動(dòng)前彈簧組成的系統(tǒng)處于平衡狀態(tài)。在網(wǎng)格邊界節(jié)點(diǎn)發(fā)生位移后，會(huì)產(chǎn)生與位移成比例的力，從而由力轉(zhuǎn)化出網(wǎng)格節(jié)點(diǎn)移動(dòng)的距離。彈簧光順方法并不將網(wǎng)格質(zhì)量作為網(wǎng)格更新的判斷因素，而只是單純的將邊界運(yùn)動(dòng)作為參數(shù)擴(kuò)散至計(jì)算域中，因此無需插值，網(wǎng)格拓?fù)涫冀K不變。
網(wǎng)格在ICEM中生成且全都為六面體結(jié)構(gòu)化網(wǎng)格，總體網(wǎng)格數(shù)為800萬。網(wǎng)格未發(fā)生移動(dòng)時(shí)，在ICEM中網(wǎng)格質(zhì)量2項(xiàng)評(píng)價(jià)標(biāo)準(zhǔn)如下：Determinant2′2′2>0.6，網(wǎng)格角度Angle>18°。使用Fluent計(jì)算時(shí)，當(dāng)網(wǎng)格質(zhì)量Determinant2′2′2>0.2且網(wǎng)格角度Angle>9°就可以確保數(shù)值計(jì)算順利進(jìn)行，由此可見初始時(shí)網(wǎng)格質(zhì)量非常高。為保證閥出口為均勻流和避免回流的產(chǎn)生，將閥出口加長(zhǎng)到出口直徑的10倍，加密了喉部和閥碟下方流動(dòng)變化較為劇烈的區(qū)域網(wǎng)格，保證在此區(qū)域內(nèi)網(wǎng)格均勻變化和邊界層內(nèi)y+<=5。
圖3為沿東西方向中心截面網(wǎng)格和網(wǎng)格變形后的局部示意圖。由圖3可以看出，在閥座喉部區(qū)域網(wǎng)格很密，閥碟周圍**層網(wǎng)格厚度為0.004mm。若不加以特殊處理，閥碟移動(dòng)極易使邊界層網(wǎng)格畸變過大而使計(jì)算發(fā)散�？紤]到邊界層內(nèi)流動(dòng)對(duì)流場(chǎng)的巨大影響，同時(shí)為保護(hù)閥碟處邊界層網(wǎng)格，將包圍閥碟部分區(qū)域網(wǎng)格設(shè)置為剛體區(qū)，即在計(jì)算中這部分網(wǎng)格不會(huì)變形。而此區(qū)域外的網(wǎng)格設(shè)置為變形區(qū)，即在閥碟移動(dòng)中這部分網(wǎng)格會(huì)產(chǎn)生相應(yīng)的變形，剛體區(qū)域和變形區(qū)域的劃分如圖3所示。通過控制彈簧光順法中兩個(gè)參數(shù)：彈簧常數(shù)因子與邊界節(jié)點(diǎn)松弛因子，將網(wǎng)格中較大的變形控制在變形區(qū)域中靠近剛體區(qū)外圍，而保證靠近閥座處較密的網(wǎng)格基本不變。由圖3可以看出，當(dāng)閥碟頂端偏向東側(cè)1mm后喉部網(wǎng)格并沒有發(fā)生畸變。將閥碟頂端偏向東側(cè)1mm的網(wǎng)格重新導(dǎo)入ICEM中，其網(wǎng)格質(zhì)量為：Determinant2′2′2>0.25，網(wǎng)格角度Angle>10°，足以滿足Fluent計(jì)算中對(duì)網(wǎng)格的要求動(dòng)網(wǎng)格常用于調(diào)節(jié)閥在開啟或關(guān)閉計(jì)算中，此時(shí)由于閥碟位移較大使得網(wǎng)格變形量增大，為保證網(wǎng)格在大變形后畸變率低不得已采取較粗的網(wǎng)格，邊界層網(wǎng)格通常難以嚴(yán)格滿足湍流模型的需求。在中小升程下調(diào)節(jié)閥內(nèi)流動(dòng)往往為復(fù)雜的三維超音速流動(dòng)，粗網(wǎng)格通常難以精確的捕捉激波、射流分離和邊界層內(nèi)的諸多流動(dòng)細(xì)節(jié)，因此所得的結(jié)果往往只能定性的分析閥內(nèi)的流動(dòng)情況。本文在數(shù)值計(jì)算中所采用DES湍流模型雖對(duì)流過復(fù)雜機(jī)構(gòu)的高雷諾數(shù)分離流具有非常高的預(yù)測(cè)能力，但是DES湍流模型高度的依賴網(wǎng)格尺度、網(wǎng)格質(zhì)量和邊界層內(nèi)的網(wǎng)格布置。針對(duì)閥碟小幅振動(dòng)特性，借助Fluent中彈簧光順模型，創(chuàng)新性的將邊界層網(wǎng)格剛體化，將閥碟移動(dòng)造成的網(wǎng)格變形“轉(zhuǎn)移”到相對(duì)尺寸較大的區(qū)域，既實(shí)現(xiàn)了800萬結(jié)構(gòu)網(wǎng)格網(wǎng)格變形也滿足了湍流模型對(duì)網(wǎng)格的需求。


圖3 沿東西方向中心截面初始網(wǎng)格和局部網(wǎng)格變形圖
3.數(shù)值計(jì)算結(jié)果
在閥碟下表面沿振動(dòng)的方向上設(shè)置了兩個(gè)監(jiān)測(cè)點(diǎn)A1和A2，其位置如圖1中所示。圖4為閥碟靜止?fàn)顟B(tài)下，閥碟上兩測(cè)點(diǎn)隨時(shí)間的壓力變化時(shí)域圖，可以看出兩測(cè)點(diǎn)處的壓力呈現(xiàn)無規(guī)則的波動(dòng)狀態(tài)。圖5和圖6是閥桿系統(tǒng)以頻率f=46.4Hz振動(dòng)閥碟頂端振幅分別為h=0.6mm和h=1mm時(shí)的閥碟上2個(gè)測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖，2圖中右側(cè)同時(shí)標(biāo)明了閥碟位移方向。由圖5可以看出當(dāng)閥碟頂端振幅為0.6mm時(shí)，閥碟上壓力測(cè)點(diǎn)波動(dòng)仍呈現(xiàn)隨機(jī)狀態(tài)，并未受到閥碟運(yùn)動(dòng)的影響。而由圖6可以看出，當(dāng)閥碟頂端振幅為1mm時(shí)，閥碟監(jiān)測(cè)點(diǎn)壓力明顯出現(xiàn)了隨閥碟位移變化的周期性波動(dòng)現(xiàn)象。當(dāng)閥碟運(yùn)動(dòng)到東側(cè)時(shí)，位于東側(cè)的監(jiān)測(cè)點(diǎn)A1的波動(dòng)明顯增大；當(dāng)閥碟運(yùn)動(dòng)到西側(cè)時(shí)，位于西側(cè)的監(jiān)測(cè)點(diǎn)A2的波動(dòng)明顯增大。

圖4 兩監(jiān)測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖（閥碟固定）

圖5 兩監(jiān)測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖（h=0.6mm，f=46.4Hz）

圖6 兩監(jiān)測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖（h=1mm，f=46.4Hz）
圖7為在計(jì)算中調(diào)節(jié)閥沿進(jìn)口方向中心截面上三種典型的流型�？梢钥闯鰪膬蓚�(cè)喉部流出的高速射流在閥碟下方區(qū)域某點(diǎn)撞擊后向出口流動(dòng)�？缫羲僮矒羯淞髯陨淼牟环�(wěn)定性使得喉部以后的射流處于不斷的擺動(dòng)狀態(tài)，撞擊點(diǎn)也在閥碟下方不斷移動(dòng)。當(dāng)由撞擊形成的高壓區(qū)移動(dòng)到監(jiān)測(cè)點(diǎn)附近時(shí)就會(huì)導(dǎo)致監(jiān)測(cè)點(diǎn)處壓力脈動(dòng)增加。當(dāng)閥碟處于靜止?fàn)顟B(tài)時(shí)，由兩側(cè)喉部流出的高速射流的強(qiáng)度相當(dāng)，高速射流在閥碟下方撞擊后不會(huì)出現(xiàn)射流角度發(fā)生較大偏斜的流型，撞擊位置相對(duì)遠(yuǎn)離閥碟下方，所以閥碟監(jiān)測(cè)點(diǎn)處不會(huì)出現(xiàn)明顯的壓力波動(dòng)。而當(dāng)閥碟在振動(dòng)中偏向一側(cè)時(shí)（如圖2中閥碟偏向東側(cè)為例），由西側(cè)較寬的喉部處有更多的流體流過，此時(shí)由西側(cè)喉部流出的流體具有更多的能量推動(dòng)由東側(cè)喉部流出的高速射流偏向東側(cè)，撞擊點(diǎn)也會(huì)更接近監(jiān)測(cè)點(diǎn)A1處，導(dǎo)致A1處壓力波動(dòng)增大。同樣當(dāng)閥碟偏向西側(cè)時(shí)，高速射流撞擊點(diǎn)位置也會(huì)轉(zhuǎn)移到A2附近，所以A2的壓力波動(dòng)也會(huì)相應(yīng)增加。

圖7 調(diào)節(jié)閥沿進(jìn)口中心面Ma圖
由圖4—6可以看出，當(dāng)閥碟頂端的振幅小于等于0.6mm時(shí)，數(shù)值模擬此調(diào)節(jié)閥內(nèi)流動(dòng)可以不考慮閥桿系統(tǒng)振動(dòng)對(duì)流場(chǎng)的影響。但當(dāng)閥碟頂端振幅達(dá)到1mm時(shí)，閥內(nèi)的流動(dòng)明顯受到閥桿系統(tǒng)振動(dòng)的影響，此時(shí)須對(duì)其加以考慮。閥碟頂端振幅主要由閥碟和閥碟套間隙決定，當(dāng)這一間隙足夠大時(shí)，閥桿系統(tǒng)的振動(dòng)就會(huì)引起調(diào)節(jié)閥內(nèi)流場(chǎng)周期性波動(dòng)，形成流固耦合振動(dòng)，這種振動(dòng)往往嚴(yán)重威脅調(diào)節(jié)閥的安全。
圖8和圖9為閥碟頂端振幅為h=1mm，振動(dòng)頻率分別為f=92.8Hz和f=185.6Hz時(shí)的閥碟兩測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖。當(dāng)閥桿系統(tǒng)按照92.8Hz振動(dòng)時(shí)，兩測(cè)點(diǎn)壓力脈動(dòng)仍隨閥碟位移有明顯的周期變化，但壓力波動(dòng)已出現(xiàn)了滯后現(xiàn)象。當(dāng)按照185.6Hz振動(dòng)時(shí)，閥碟測(cè)點(diǎn)處雖出現(xiàn)了較大幅度的壓力脈動(dòng)，但和位移已沒有明顯的周期對(duì)應(yīng)關(guān)系。這是由于兩側(cè)高速射流的擺動(dòng)存在慣性，隨著振動(dòng)頻率的增加兩側(cè)射流難以及時(shí)跟隨閥碟一起運(yùn)動(dòng)，所以在監(jiān)測(cè)點(diǎn)處也不會(huì)呈現(xiàn)明顯的隨閥碟移動(dòng)的壓力波動(dòng)。

圖8 兩監(jiān)測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖（h=1mm，f=92.8Hz）

圖9 兩監(jiān)測(cè)點(diǎn)壓力波動(dòng)時(shí)域圖（h=1mm，f=185.6Hz）
4.結(jié)論
本文利用FLUNT中動(dòng)網(wǎng)格模型和UDF函數(shù)，考慮了閥碟在不同振動(dòng)幅值和振動(dòng)頻率下對(duì)閥內(nèi)流場(chǎng)的影響，得到如下結(jié)論：
1）在動(dòng)網(wǎng)格變形中，采取固定閥碟邊界層處附近的網(wǎng)格使網(wǎng)格在尺度相對(duì)較大的位置移動(dòng)變形。避免了以往動(dòng)網(wǎng)格變形中由于邊界層網(wǎng)格畸變過大而導(dǎo)致的計(jì)算發(fā)散。這種高質(zhì)量的網(wǎng)格變化方法對(duì)調(diào)節(jié)閥內(nèi)利用動(dòng)網(wǎng)格模擬閥碟移動(dòng)對(duì)閥內(nèi)流場(chǎng)的影響有重要意義。
2）動(dòng)網(wǎng)格計(jì)算表明當(dāng)閥碟頂端振動(dòng)幅值小于等于0.6mm時(shí)，閥碟上壓力波動(dòng)為隨機(jī)波動(dòng)；當(dāng)閥碟頂端振動(dòng)幅值達(dá)到1mm時(shí)，閥碟上壓力波動(dòng)呈現(xiàn)隨閥碟位移周期性變化的特點(diǎn)，此時(shí)數(shù)值計(jì)算中應(yīng)考慮到閥桿系統(tǒng)振動(dòng)對(duì)流場(chǎng)的影響。
3）隨著閥桿系統(tǒng)振動(dòng)頻率的增加，由于主流流型變化存在慣性，閥碟上的壓力波動(dòng)隨閥碟位移出現(xiàn)滯后現(xiàn)象直至周期性波動(dòng)消失。數(shù)值模擬了閥桿系統(tǒng)按照不同振幅和振動(dòng)頻率對(duì)閥內(nèi)流場(chǎng)的影響，對(duì)后續(xù)數(shù)值模擬中是否考慮振動(dòng)對(duì)閥內(nèi)流場(chǎng)的影響提供了參考。
文章鏈接：http://www.klevalve.com/jsxw/1055.htm


資訊來源：無錫科萊恩流體控制設(shè)備有限公司