• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    某無刷勵(lì)磁機(jī)通風(fēng)冷卻數(shù)值模擬研究

    2016-07-16 03:00:37路義萍張東學(xué)孫博王佐民李夢(mèng)啟
    關(guān)鍵詞:數(shù)值模擬

    路義萍, 張東學(xué), 孫博, 王佐民, 李夢(mèng)啟

    (1.哈爾濱理工大學(xué) 機(jī)械動(dòng)力工程學(xué)院,黑龍江 哈爾濱 150080;2.哈爾濱電氣動(dòng)力裝備有限公司,黑龍江 哈爾濱 150040)

    ?

    某無刷勵(lì)磁機(jī)通風(fēng)冷卻數(shù)值模擬研究

    路義萍1,張東學(xué)1,孫博1,王佐民1,李夢(mèng)啟2

    (1.哈爾濱理工大學(xué) 機(jī)械動(dòng)力工程學(xué)院,黑龍江 哈爾濱 150080;2.哈爾濱電氣動(dòng)力裝備有限公司,黑龍江 哈爾濱 150040)

    摘要:為了研究隱極同步電動(dòng)機(jī)配套用無刷勵(lì)磁機(jī)的通風(fēng)方案中冷卻空氣的流動(dòng)路徑及冷卻特點(diǎn),依據(jù)流體動(dòng)力學(xué)原理理論,采用有限體積法,建立了某新型無刷勵(lì)磁機(jī)相對(duì)地面結(jié)構(gòu)的1/2為物理模型,研究了無刷勵(lì)磁機(jī)通風(fēng)冷卻方案的流動(dòng)路徑、流場(chǎng)和溫度場(chǎng),并給出了通風(fēng)方案的冷卻特點(diǎn)及流場(chǎng)和溫度場(chǎng)的分布特點(diǎn)。計(jì)算結(jié)果表明,該勵(lì)磁機(jī)方案的最高溫度位置出現(xiàn)在靠近整流盤側(cè)的轉(zhuǎn)子線棒端部,絕緣均不超溫。與某電機(jī)研究所采用等效風(fēng)路法所得該方案相同工況下計(jì)算值相比,該通風(fēng)方案數(shù)值模擬進(jìn)風(fēng)量誤差為5.6%,為進(jìn)一步優(yōu)化該勵(lì)磁機(jī)冷卻風(fēng)路系統(tǒng)提供了理論依據(jù)。

    關(guān)鍵詞:無刷勵(lì)磁機(jī);通風(fēng)方案;數(shù)值模擬;有限體積法;熱流場(chǎng)

    0引言

    近年來我國(guó)西氣東輸一線和二線管道工程已建成并投入輸氣,95%以上天然氣通過長(zhǎng)輸管線輸送到內(nèi)地,沿管道需增壓,壓縮機(jī)由高速隱極同步電動(dòng)機(jī)驅(qū)動(dòng),為長(zhǎng)輸管道持續(xù)加壓。與隱極同步電動(dòng)機(jī)配套的無刷勵(lì)磁機(jī)為機(jī)組的重要部分。無刷勵(lì)磁機(jī)不僅僅是實(shí)現(xiàn)能量轉(zhuǎn)換必不可少的部分,對(duì)隱極同步電動(dòng)機(jī)機(jī)組的安全、高效率持續(xù)運(yùn)行也尤為重要。無刷勵(lì)磁機(jī)的絕緣或者是整流盤上二極管工作環(huán)境溫度超過允許極限,都將影響勵(lì)磁機(jī)的壽命和機(jī)組的安全運(yùn)行,機(jī)組研發(fā)過程中,勵(lì)磁機(jī)與主機(jī)部分的通風(fēng)冷卻問題同等重要。因此,了解勵(lì)磁機(jī)內(nèi)部的通風(fēng)冷卻方式,對(duì)指導(dǎo)勵(lì)磁機(jī)的設(shè)計(jì)及運(yùn)行必不可少。

    近十幾年來,國(guó)內(nèi)外學(xué)者一直重視對(duì)勵(lì)磁機(jī)的研究。例如,文獻(xiàn)[1]介紹了空冷系列汽輪發(fā)電機(jī)無刷勵(lì)磁機(jī)的設(shè)計(jì)理念、結(jié)構(gòu)特色、關(guān)鍵工藝及投產(chǎn)應(yīng)用情況,文獻(xiàn)[2]對(duì)大容量勵(lì)磁機(jī)離心式風(fēng)扇的結(jié)構(gòu)做了改進(jìn),文獻(xiàn)[3]介紹了100MW汽輪發(fā)電機(jī)及勵(lì)磁機(jī)通風(fēng)系統(tǒng)的改進(jìn),文獻(xiàn)[4]對(duì)勵(lì)磁機(jī)的流場(chǎng)及溫度場(chǎng)進(jìn)行了數(shù)值模擬和分析,文獻(xiàn)[5]對(duì)11相無刷勵(lì)磁機(jī)進(jìn)行了額定工況及故障工況的仿真研究,文獻(xiàn)[6]對(duì)無刷勵(lì)磁機(jī)同步電動(dòng)機(jī)異常振動(dòng)進(jìn)行了故障分析,文獻(xiàn)[7]對(duì)三相或多相無刷勵(lì)磁機(jī)及旋轉(zhuǎn)整流系統(tǒng)故障診斷進(jìn)行了研究,文獻(xiàn)[8]給出了用非破壞性方法來檢測(cè)無刷勵(lì)磁機(jī)的整流盤故障,文獻(xiàn)[9]說明了旋轉(zhuǎn)無刷勵(lì)磁機(jī)到靜態(tài)無刷勵(lì)磁機(jī)的主要?jiǎng)?lì)磁系統(tǒng)的改造過程。由于整流二極管處于高速旋轉(zhuǎn)的狀態(tài),在實(shí)際運(yùn)行中,二極管發(fā)生損壞的情況時(shí)有發(fā)生[10]。當(dāng)勵(lì)磁機(jī)運(yùn)行時(shí),如果經(jīng)常發(fā)生磁極線圈過熱燒損現(xiàn)象[11],將嚴(yán)重威脅機(jī)組的安全穩(wěn)定運(yùn)行。目前,針對(duì)隱極電動(dòng)機(jī)用無刷勵(lì)磁機(jī)的熱流場(chǎng)研究方法,與大型發(fā)電機(jī)、風(fēng)力發(fā)電機(jī)、汽輪發(fā)電機(jī)的一樣,大多數(shù)采用有限元法[12-14]、熱網(wǎng)絡(luò)法[15-16]、有限體積法[17-18]等,這幾種方法在當(dāng)今應(yīng)用比較廣泛,文獻(xiàn)[19]中BOGLIETTI A等綜合分析了熱網(wǎng)絡(luò)法、有限元法、有限體積法在電機(jī)溫度場(chǎng)計(jì)算時(shí)的優(yōu)缺點(diǎn),但在國(guó)內(nèi)外用實(shí)驗(yàn)方法來解決這類問題的文獻(xiàn)較為少見。

    本文研究的隱極同步電動(dòng)機(jī)用無刷勵(lì)磁機(jī)屬于國(guó)外引進(jìn)技術(shù)消化吸收預(yù)研產(chǎn)品,其通風(fēng)方案中冷卻空氣的流動(dòng)路徑及冷卻特點(diǎn)還處于探索階段。本文以計(jì)算流體動(dòng)力學(xué)理論作為基礎(chǔ),主要研究了出口位置在整流盤正對(duì)上方的方案,對(duì)其流場(chǎng)及溫度場(chǎng)進(jìn)行了CFD數(shù)值計(jì)算,得到了流場(chǎng)和溫度場(chǎng)的模擬結(jié)果,并驗(yàn)證了其準(zhǔn)確性。

    1物理模型

    研究的無刷勵(lì)磁機(jī)轉(zhuǎn)軸由高強(qiáng)度合金鋼鍛件加工而成,轉(zhuǎn)軸上銑出鍵槽及通風(fēng)槽。電樞鐵心是由高強(qiáng)度、抗疲勞的硅鋼片疊壓而成,兩端用壓圈固緊。電樞繞組由成型的線圈構(gòu)成,絕緣股線繞制成線圈之后,電樞繞組端部用無緯帶綁扎環(huán)支撐和固定。該無刷勵(lì)磁機(jī)主要由定子、轉(zhuǎn)子、轉(zhuǎn)軸、雙側(cè)整流盤及機(jī)座構(gòu)成,見圖1。定轉(zhuǎn)子沿軸向共分布6個(gè)風(fēng)溝,風(fēng)溝的寬度為10 mm,靠近旋轉(zhuǎn)整流盤側(cè)風(fēng)溝標(biāo)號(hào)為1,其他2~6見圖1。

    1~6—定子風(fēng)溝;7—機(jī)殼;8—機(jī)座;9—轉(zhuǎn)軸;10—雙側(cè)整流盤;11—機(jī)座環(huán)板;12—定子;13—轉(zhuǎn)子;14—二極管;15—軸套;16—絕緣墊板。圖1 計(jì)算域物理模型圖Fig.1 Computational domain of the physical model

    在建模及計(jì)算的過程中考慮了如下幾點(diǎn):1) 由于勵(lì)磁機(jī)通風(fēng)結(jié)構(gòu)的對(duì)稱性,故計(jì)算域選取整機(jī)沿軸的周向180°,即相對(duì)地面結(jié)構(gòu)的1/2;2)用Solid works軟件建立材料物性不同的固體模型;3)將無刷勵(lì)磁機(jī)的流體區(qū)域置于旋轉(zhuǎn)坐標(biāo)系和靜止坐標(biāo)系中。這使得三維流場(chǎng)、溫度場(chǎng)計(jì)算結(jié)果更加準(zhǔn)確可靠。

    2數(shù)學(xué)模型及求解條件

    2.1數(shù)學(xué)模型

    通過計(jì)算得出入口處空氣的雷諾數(shù)Re遠(yuǎn)大于2 300,處于湍流狀態(tài)。計(jì)算域內(nèi)空氣視為不可壓縮流體,利用Fluent軟件,選擇多重參考系計(jì)算方法。其轉(zhuǎn)子、整流盤固體區(qū)域中的空氣處于旋轉(zhuǎn)坐標(biāo)系中,其他流體部分位于三維固定直角坐標(biāo)系中。

    在溫度場(chǎng)計(jì)算求解過程中采用質(zhì)量、動(dòng)量與能量守恒方程及標(biāo)準(zhǔn)k-ε兩方程湍流模型[4],其通用控制方程如下:

    div(ρUφ)=div(Γφgradφ)+Sφ。

    (1)

    式中:ρ為密度;φ為通用變量;Γφ、Sφ分別為廣義擴(kuò)散系數(shù)和廣義源項(xiàng)。

    2.2計(jì)算條件

    在額定工況下,轉(zhuǎn)子轉(zhuǎn)速為3 120 r/min,某電機(jī)研究所通過電磁場(chǎng)計(jì)算得到了各種損耗數(shù)值,見表1。由表1可算得計(jì)算域內(nèi)各部分熱源數(shù)值。空氣所受重力與其他力相比較小,可以忽略不計(jì)。入出口風(fēng)筒均接至主機(jī)風(fēng)扇前端,因此,其表壓力均為0 Pa,空氣入口溫度按照電機(jī)規(guī)范取40℃。物理模型外周界機(jī)殼為壁面邊界,中心剖面轉(zhuǎn)子部分為周期性邊界。在計(jì)算過程中,匝間絕緣材料、銅繞組、軸均為各向同性材料,物性參數(shù)為常數(shù),熱導(dǎo)率λ分別為0.22、387.6、31.8 W/(m·K),而鐵心疊片為各向異性材料,熱導(dǎo)率λ沿疊片徑、切向的數(shù)值都為42.5 W/(m·K),軸向(厚度方向)為0.57 W/(m·K)。

    表1 勵(lì)磁機(jī)內(nèi)部件損耗值

    物理模型建完以后,需對(duì)其網(wǎng)格進(jìn)行劃分。在近壁面區(qū)采用標(biāo)準(zhǔn)的壁面函數(shù)法進(jìn)行處理,近壁面網(wǎng)格y+滿足標(biāo)準(zhǔn)壁面函數(shù)要求;方程組采用分離、隱式求解,壓力速度耦合方程組采用SIMPLE算法,方程均采用二階迎風(fēng)格式離散,經(jīng)多次網(wǎng)格調(diào)整,最終獲得計(jì)算域的穩(wěn)態(tài)湍流流動(dòng)及傳熱耦合計(jì)算網(wǎng)格獨(dú)立收斂解。

    3計(jì)算結(jié)果及分析

    3.1速度場(chǎng)計(jì)算結(jié)果及分析

    設(shè)計(jì)階段,廠商認(rèn)為定子上方風(fēng)筒為勵(lì)磁機(jī)出風(fēng)筒,現(xiàn)場(chǎng)實(shí)驗(yàn)及數(shù)值模擬實(shí)驗(yàn)發(fā)現(xiàn),現(xiàn)實(shí)與想象相反,整流盤上方的出口為熱風(fēng),屬于空氣出口。經(jīng)過流場(chǎng)反演,得到這種同軸布置3個(gè)旋轉(zhuǎn)部件(雙側(cè)整流盤+轉(zhuǎn)子)的電機(jī)內(nèi)部空氣流動(dòng)路徑,見圖2。圖2給出了兩個(gè)特征截面極角10°、120°速度矢量圖。

    圖2 極角10°與120°截面位置速度矢量圖Fig.2    Velocity vector diagram of polar angle    10° and 120° section

    由圖2可知冷卻空氣的具體走向,具體如下:取自于主機(jī)風(fēng)扇前端的流道中的冷空氣經(jīng)定子上方的進(jìn)風(fēng)口進(jìn)入定子鐵心背部,大部分經(jīng)過定子鐵心邊段風(fēng)溝通過機(jī)座環(huán)板,進(jìn)入兩側(cè)的空氣區(qū)內(nèi);另一小部分進(jìn)入定子風(fēng)溝中向下流動(dòng),與經(jīng)轉(zhuǎn)軸通風(fēng)槽到達(dá)轉(zhuǎn)子風(fēng)溝甩出的空氣在氣隙旋轉(zhuǎn)射流混合,之后到左右兩側(cè)空氣區(qū)中。從定子邊段鐵心風(fēng)溝甩出靠近風(fēng)溝1的空氣和從轉(zhuǎn)子風(fēng)溝1甩出的在整流盤引力作用下經(jīng)氣隙沿軸向流出的空氣,被右側(cè)整流盤吸入,冷卻旋轉(zhuǎn)整流盤二極管;另一部分自軸與機(jī)座間的空隙進(jìn)入外罩與機(jī)座之間,冷卻機(jī)座與外罩后被左側(cè)整流盤吸入,冷卻其旋轉(zhuǎn)整流盤二極管。冷卻空氣從左右雙側(cè)整流盤上部出風(fēng)口流出,進(jìn)入主機(jī)風(fēng)扇前端,在主機(jī)內(nèi)循環(huán)。分析速度矢量圖可知:計(jì)算域內(nèi)的轉(zhuǎn)子風(fēng)溝空氣區(qū)、整流盤空氣區(qū)等的速度非常高,最高達(dá)到81.9 m/s,位置在旋轉(zhuǎn)半徑較大的整流盤處,氣隙區(qū)靠近轉(zhuǎn)子外表面處空氣受到壁面旋轉(zhuǎn)剪切粘性力的帶動(dòng),速度也較大。在機(jī)座、機(jī)殼和定子風(fēng)溝2~4附近的冷卻空氣速度一般較低,低于5.6 m/s。在旋轉(zhuǎn)部件周圍會(huì)形成大小不一的渦旋,產(chǎn)生摩擦生熱,不利于空氣的流動(dòng)和各部件的冷卻。

    計(jì)算結(jié)果表明,轉(zhuǎn)子內(nèi)部旋轉(zhuǎn)半徑越大,空氣運(yùn)動(dòng)速度越大,定子風(fēng)溝內(nèi)部流動(dòng)復(fù)雜,能夠體現(xiàn)該種風(fēng)路布置的勵(lì)磁機(jī)內(nèi)部流動(dòng)特征,故圖3給出了該方案定子風(fēng)溝3(Z=25 mm,r=200~295 mm)沿周向0~90°內(nèi)的速度矢量圖補(bǔ)充說明勵(lì)磁機(jī)內(nèi)速度分布特征。

    圖3 定子風(fēng)溝3速度矢量圖Fig.3    Velocity vector diagram of wind duct 3    of stator part

    由圖3可知,氣隙部分上部靠近定子側(cè)與下部靠近轉(zhuǎn)子部分的速度大小和方向顯著不同,界面明顯,靠近轉(zhuǎn)子外徑處的空氣在粘性剪切應(yīng)力作用下,一同旋轉(zhuǎn),形成速度較大的流體邊界層,最大速度為47.9 m/s,隨著氣隙沿半徑增大方向,速度數(shù)值逐漸減小,速度梯度很大,且在與定子齒部相接處,空氣部分流進(jìn)、部分流出。一小部分向Y軸正方向進(jìn)入定子線棒周圍較小截面的鐵心齒部風(fēng)溝中,空氣速度數(shù)值大多在10.5 m/s以下,相對(duì)較?。豢拷B接筋(工字鋼)兩邊來自齒部的空氣速度較大,空氣流動(dòng)較復(fù)雜,定子線棒周邊的齒部與軛部中心位置空氣有明顯的渦流形成,冷卻空氣圍繞固體線棒旋轉(zhuǎn),與背部相連處的定子風(fēng)溝上方冷卻空氣有進(jìn)有出,定子鐵心軛部上邊緣外為定子鐵心背部空氣區(qū),風(fēng)溝軛部中每條定子線棒上方的空氣流通截面積沿徑向逐漸增大,定子疊片連接筋向心布置,因而形成中心局部速度較小。

    3.2壓力場(chǎng)計(jì)算結(jié)果及分析

    圖4給出了該方案極角30°、極角150°特征截面的壓力分布云圖進(jìn)行對(duì)比分析。

    由圖4可知,總體而言,轉(zhuǎn)軸上方轉(zhuǎn)子槽空氣區(qū)處的壓力最低,壓力值為-2 721 Pa,緊挨機(jī)座左側(cè)入口的整流盤靠近軸套處空氣區(qū)負(fù)壓較低,壓力值-1 886~-1 052 Pa之間,吸風(fēng)作用較強(qiáng),極角150°截面右側(cè)整流盤二極管下方靠近軸套和絕緣墊板之間局部空氣區(qū)壓力最高,壓力值為1 450 Pa,由于其間隙較小,故會(huì)出現(xiàn)憋風(fēng)現(xiàn)象,壓力較高,最低值與最高值壓差4 171 Pa,機(jī)殼與機(jī)座之間空腔中負(fù)壓較低,數(shù)值在-1 052~-635 Pa之間,轉(zhuǎn)子風(fēng)溝空氣區(qū)壓力在-1 886~-218 Pa之間,沿徑向壓力逐漸升高,主要由于轉(zhuǎn)軸、轉(zhuǎn)子鐵心和整流盤高速旋轉(zhuǎn),形成負(fù)壓區(qū),使定子上方進(jìn)風(fēng)口處的空氣順利進(jìn)入定轉(zhuǎn)子風(fēng)溝,勵(lì)磁機(jī)內(nèi)空氣入出口無靜壓差,除定子風(fēng)溝2~4入出口壓力差很小,由于定轉(zhuǎn)子風(fēng)溝中間通過厚度僅3 mm的環(huán)形氣隙空間連接,氣隙中空氣與軸向定轉(zhuǎn)子兩端空氣區(qū)存在壓差,導(dǎo)致氣隙內(nèi)空氣主要沿軸向流動(dòng)。

    圖4 極角30°與150°截面位置壓力分布云圖Fig.4    Static pressure contours of polar angle    30° and 150° section

    3.3三維溫度場(chǎng)分析

    以勵(lì)磁機(jī)該方案極角10°、120°特征截面的溫度分布云圖5為例,說明勵(lì)磁機(jī)定轉(zhuǎn)子、絕緣、空氣溫度分布特點(diǎn)。

    圖5 極角10°與120°截面位置溫度分布云圖Fig.5    Temperature contours of polar angle   10° and 120° section

    由圖5可看出,最高溫度位置出現(xiàn)在靠近整流盤的轉(zhuǎn)子線棒端部處,數(shù)值為107.1℃,顯然,與之相鄰的轉(zhuǎn)子端部線棒絕緣溫度遠(yuǎn)低于許用溫度130℃,而且整流盤二極管附近空氣溫度范圍為70~73℃。定子峰值溫度出現(xiàn)在線棒端部靠近整流盤部分,數(shù)值93℃,定子線棒溫度分布左右不對(duì)稱,另一端部溫度較低,數(shù)值范圍為76~80℃,左右溫度差為13℃左右,原因是遠(yuǎn)離整流盤端的邊段鐵心風(fēng)溝中冷卻空氣流量大。定子鐵心背部靠近空氣入口的部分溫度較低,靠近定子線棒附近齒部溫度較高。原因是冷空氣剛進(jìn)入定子鐵心背部時(shí),溫度較低,齒部除空氣溫度高外,位置方面靠近溫度較高的定子線棒。機(jī)殼與機(jī)座之間,部分靠近出口處的空氣溫度較高,溫度范圍為66~70℃,機(jī)座下方溫度較低,溫度范圍為63~66℃。

    圖6為該方案補(bǔ)充和更詳細(xì)展示轉(zhuǎn)子線棒的溫度分布,由圖可知,轉(zhuǎn)子線棒部分的峰值溫度為107.1℃,位置同上,遠(yuǎn)離整流盤側(cè)轉(zhuǎn)子線棒端部溫度為91~99℃,靠近整流盤側(cè)轉(zhuǎn)子線棒端部的溫度為97~107℃,位于遠(yuǎn)離整流盤側(cè)的轉(zhuǎn)子線棒端部溫度低于靠近整流盤側(cè)的,主要原因是靠近整流盤側(cè)有轉(zhuǎn)子線圈架緊鄰,線圈架材料導(dǎo)熱性能差,另外下部軸空氣槽入口部分被轉(zhuǎn)子引線阻擋,該側(cè)槽入口進(jìn)風(fēng)量較少,導(dǎo)致轉(zhuǎn)子線棒直段的溫度從左到右逐漸降低。此外,被空氣直接沖刷冷卻的風(fēng)溝位置處的繞組溫度低于周邊與鐵心疊片相鄰部分的繞組溫度,見圖6。

    圖6 轉(zhuǎn)子線棒溫度分布云圖Fig.6 Temperature contour of rotor winding bar

    3.4結(jié)果準(zhǔn)確性分析

    勵(lì)磁機(jī)內(nèi)部旋轉(zhuǎn)流場(chǎng)與空氣動(dòng)力場(chǎng)很復(fù)雜,目前不能獲得解析解。為了驗(yàn)證結(jié)果的準(zhǔn)確性,一方面,將該方案的體積流量為0.283 m3/s,與某電機(jī)研究所采用等效風(fēng)路法所得該方案相同工況下的總通風(fēng)量0.3 m3/s做了對(duì)比,誤差為5.6%比較準(zhǔn)確。另一方面,數(shù)值模擬出的空氣入出口流動(dòng)方向與實(shí)驗(yàn)結(jié)果測(cè)量結(jié)果一致。

    4結(jié)論

    本文通過對(duì)同軸上布置雙側(cè)整流盤和轉(zhuǎn)子三個(gè)旋轉(zhuǎn)部件、且整流盤和定子上方均布置風(fēng)筒的新型無刷勵(lì)磁機(jī)三維湍流流場(chǎng)和溫度場(chǎng)的數(shù)值計(jì)算反演,得出以下結(jié)論:

    1)冷卻空氣從定子上方風(fēng)筒入口進(jìn)入,大部分從定子鐵心邊段風(fēng)溝流出,從中部定子風(fēng)溝流出較少,在鐵心齒部線棒周圍及軛部空氣均形成渦流。從轉(zhuǎn)子風(fēng)溝甩出的空氣,在整流盤負(fù)壓引力作用下經(jīng)氣隙沿軸向流出。

    2)轉(zhuǎn)軸上方空氣槽中空氣負(fù)壓最低,緊挨機(jī)座入口的整流盤靠近軸套處空氣區(qū)負(fù)壓較低,吸風(fēng)作用較強(qiáng)。

    3)該方案勵(lì)磁機(jī)最高溫度位置出現(xiàn)在靠近整流盤側(cè)的轉(zhuǎn)子線棒端部,左右兩側(cè)溫度不對(duì)稱,轉(zhuǎn)子線棒直段的溫度從左到右逐漸降低,絕緣均不超溫,并且整流盤的工作環(huán)境溫度較低,符合要求。

    參 考 文 獻(xiàn):

    [1]孫維炎,黃浩.空冷系列汽輪發(fā)電機(jī)上的無刷勵(lì)磁機(jī)研制[J]. 上海大中型電機(jī), 2009,(2):14-16.

    SUN Weiyan, HUANG Hao. Study on the brushless exciter of air cooling steam-turbine generator[J]. Shanghai Medium and Large Electrical Machines, 2009,(2):14-16.

    [2]李治, 張燕, 王庭山, 等. 大容量勵(lì)磁機(jī)離心式風(fēng)扇的結(jié)構(gòu)改進(jìn)[J]. 電機(jī)技術(shù), 2014,(5):1-5.

    LI Zhi, ZHANG Yan, WANG Tingshan, et al. Structure improvement of centrifugal blower for large capacity exciter[J]. Electrical Machinery Technology,2014,(5):1-5 .

    [3]秦海平, 高滿林. 100 MW汽輪發(fā)電機(jī)及勵(lì)磁機(jī)通風(fēng)系統(tǒng)改進(jìn)[J]. 山西電力, 2002,(4):43-44.

    QIN Haiping, GAO Manlin. Improvement of ventilation system of 100 MW turbine generator and exciter[J]. Shan Xi Electric Power, 2002,(4): 43-44.

    [4]路義萍, 付余, 張萬全, 等. 一種新型無刷勵(lì)磁機(jī)的三維熱流場(chǎng)分析[J]. 中國(guó)電機(jī)工程學(xué)報(bào),2014, 34(27):4699-4705.

    LU Yiping, FU Yu, ZHANG Wanquan, et al. Analysis of three dimensional thermal and flow fields for a new brushless exciter[J]. Proceedings of the CSEE,2014,34(27):4699-4705.

    [5]賈小川, 勵(lì)慶孚, 陳世坤, 等. 11相無刷勵(lì)磁機(jī)故障工況數(shù)字仿真研究(一)[J]. 中國(guó)電機(jī)工程學(xué)報(bào),1998,18(2):149-152.

    JIA Xiaochuan, LI Qingfu, CHEN Shikun, et al. Analysis of fault cases of 11-phase brushless exciter(Part Ⅰ)[J]. Proceedings of the CSEE,1998,18(2):149-152.

    [6]蘭曉東, 席海濤. 無刷勵(lì)磁同步電動(dòng)機(jī)異常振動(dòng)故障分析[J]. 防爆電機(jī),2014,49(3):48-50.

    LAN Xiaodong, XI Haitao. Analysis on abnormal vibration fault of brushless excitation synchronous motor [J]. Explosion-proof electric machine,2014,49(3):48-50.

    [7]張超, 夏立. 發(fā)電機(jī)旋轉(zhuǎn)整流器故障的分形和動(dòng)態(tài)測(cè)度診斷[J]. 電機(jī)與控制學(xué)報(bào),2009,13(1):6-10.

    ZHANG Chao, XIA Li. Fault diagnosis of generator rotating rectifier based on fractal and dynamics[J]. Electric Machines and Control,2009,13(1):6-10.

    [8]MCARDLE M G, MORROW D J. Noninvasive detection of brushless exciter rotating diode failure[J]. IEEE Transactions on Energy Conversion,2004,19(2):378- 382.

    [9]RICHARD C S. Brushless rotating exciter conversion to main field static exciter system[C]∥Pulp and Paper Industry Technical Conference,Conference Record of 2013 Annual IEEE,2013,(172):1-5.

    [10]郝亮亮, 王善銘, 邱阿瑞, 等. 多相無刷勵(lì)磁系統(tǒng)旋轉(zhuǎn)整流器故障的仿真與識(shí)別[J].電工技術(shù)學(xué)報(bào),2012,27(4):58-62.

    HAO Liangliang, WANG Shanming, QIU Arui, et al. Simulation and recognition for rotary rectifier fault of multi-phase brushless excitation system[J]. Transactions of China Electro-Technical Society ,2012,27(4):58-62.

    [11]朱榮, 孫暉. 阿爾斯通TKJ 型無刷勵(lì)磁機(jī)磁極線圈燒損原因及處理對(duì)策[J]. 大電機(jī)技術(shù),2008(5):56-58.

    ZHU Rong, SUN Hui. The causes of and countermeasures against the burnt accidents of the magnetic pole coils of the TKJ-type blushless exciter of Alstom[J]. Large Electric Machine and Hydraulic Turbine,2008(5):56-58.

    [12]鮑曉華, 程曉巍, 方勇, 等. 大型異步電機(jī)定子端部繞組電磁力的研究[J]. 電機(jī)與控制學(xué)報(bào),2013,17(10):27-32.

    BAO Xiaohua, CHENG XiaoWei, FANG Yong, et al. Research on electromagnetic force in the stator end winding of large induction motor[J].Electric Machines and Control,2013,17(10):27-32.

    [13]林鵬, 戈寶軍, 陶大軍, 等. 多相電機(jī)定子繞組改進(jìn)繞組函數(shù)的有限元分析[J]. 電機(jī)與控制學(xué)報(bào),2015, 19(2):53-65.

    LIN Peng, GE Baojun, TAO Dajun, et al. Finite element analysis of improved winding function of stator winding in multiphase alternator[J]. Electric Machines and Control,2015,19(2):53-65.

    [14]路義萍, 潘慶輝, 韓家德, 等. 多風(fēng)路空冷汽輪發(fā)電機(jī)PIV流場(chǎng)測(cè)試實(shí)驗(yàn)臺(tái)研發(fā)[J]. 哈爾濱理工大學(xué)學(xué)報(bào),2014,19(4):32-37.

    LU Yiping, PAN Qinghui, HAN Jiade, et al. Development of experimental setup of PIV flow field test for multiple ventilation ducts of turbo generator with air coolant[J]. Journal of Harbin University of Science and Technology, 2014,19(4):32-37.

    [15]張建忠, 姜永將. 基于等效熱網(wǎng)絡(luò)法的定頻雙轉(zhuǎn)子永磁風(fēng)力發(fā)電機(jī)的熱分析[J]. 電工技術(shù)學(xué)報(bào),2015,30(2):88-96.

    ZHANG Jianzhong, JIANG Yongjiang. Thermal analysis of constant frequency double rotor permanent magnet generator based on equivalent thermal network method[J]. Transactions of China Electro Technical Society,2015,30(2):88-96.

    [16]李翠萍, 柴鳳, 程樹康. 冷卻水流速對(duì)汽車水冷電機(jī)溫升影響研究[J]. 電機(jī)與控制學(xué)報(bào),2012,16(9):1-8.

    LI Cuiping, CHAI Feng, CHENG Shukang. Research on the effects of cooling water velocity on temperature rise of the water-cooled motor in electric vehicles[J]. Electric Machines and Control,2012,16(9):1-8.

    [17]路義萍, 孫雪梅, 杜鵬, 等. CFD設(shè)置對(duì)凸極同步電動(dòng)機(jī)轉(zhuǎn)子溫度場(chǎng)影響[J]. 哈爾濱理工大學(xué)學(xué)報(bào),2015,20(3):92-95.

    LU Yiping, SUN Xuemei,DU Peng,et al. Effect of CFD settings on the temperature field of the rotor with salient pole of synchronous motor[J]. Journal of Harbin University of Science and Technology, 2015,20(3):92-95.

    [18]王芳,高斯博,湯文俠,等. 大型水氫氫冷卻汽輪發(fā)電機(jī)定子溫度分布[J]. 哈爾濱理工大學(xué)學(xué)報(bào), 2015, 20(1):26-30.

    WANG Fang, GAO Sibo, TANG Wenxia, et al.Temperature distribution of large water-hydrogen-hydrogen cooled turbo-generator stator[J]. Journal of Harbin University of Science and Technology, 2015, 20(1):26-30.

    [19]BOGLIETTI A,CAVAGNINO A,STATON D,et al.Evolution and modern approaches for thermal analysis of electrical machines[J]. IEEE Transactions on Industrial Electronics,2009,56(3):871-882.

    (編輯:劉素菊)

    Numerical simulation of ventilation cooling system for brushless exciter

    LU Yi-ping1,ZHANG Dong-xue1,SUN Bo1,WANG Zuo-min1,LI Meng-qi2

    (1.School of Mechanical & Power Engineering, Harbin University of Science and Technology, Harbin 150080;2.Harbin Electric and Power Equipment Company Limited, Harbin 150040, China)

    Abstract:In order to study the flow path and cooling characteristic of cooling air of ventilation scheme for the brushless exciter of non-salient pole synchronous motor, the finite volume method was adopted, according to the computational fluid dynamics (CFD) principles. A half physical model of a new brushless exciter structure that is relative to the ground was established. The flow path, flow and temperature fields of ventilation scheme of the brushless exciter were researched. The cooling characteristic of ventilation scheme and distribution feature of flow and temperature fields were given. The results show that the position of the highest temperature of the brushless exciter appear in the end of rotor bar. And the insulation temperature of brushless exciter is not overheating. Compared with the value that is calculated by the institute of electronic machinery technology through Equivalent Network Method, the errors of total volume flow rate is only 5.6%, under the same condition. The results provide theoretical basis for further optimizing the air ducts structure of the exciter.

    Keywords:brushless exciter; ventilation scheme; numerical simulation; finite volume method; flow and temperature fields

    收稿日期:2015-07-09

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(50976027);哈爾濱市科技創(chuàng)新人才基金(2013RFXXJ032)

    作者簡(jiǎn)介:路義萍(1965—),女,博士,教授,博士生導(dǎo)師,研究方向?yàn)楦黝愲姍C(jī)中流動(dòng)與傳熱數(shù)值模擬;

    通訊作者:路義萍

    DOI:10.15938/j.emc.2016.06.004

    中圖分類號(hào):TM 311

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1007-449X(2016)06-0026-06

    張東學(xué)(1990—),男,碩士,研究方向?yàn)闊o刷勵(lì)磁機(jī)中的流動(dòng)與傳熱數(shù)值模擬;

    孫博(1980—),男,博士,講師,研究方向?yàn)閭鳠崤c流動(dòng)數(shù)值模擬;

    王佐民(1971—),男,碩士,副教授,研究方向?yàn)闊犭娹D(zhuǎn)換;

    李夢(mèng)啟(1966—),男,碩士,高級(jí)工程師,研究方向?yàn)殡姍C(jī)設(shè)計(jì)。

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場(chǎng)信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測(cè)
    鉆孔灌注樁樁底沉渣對(duì)樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測(cè)試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場(chǎng)數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    欧美乱妇无乱码| 一夜夜www| 亚洲熟妇中文字幕五十中出| 欧美黄色片欧美黄色片| 床上黄色一级片| 精品乱码久久久久久99久播| 国产午夜精品久久久久久| 久久久久久久午夜电影| 国产伦人伦偷精品视频| 天天躁日日操中文字幕| 免费搜索国产男女视频| 一本一本综合久久| 88av欧美| 哪里可以看免费的av片| 免费看美女性在线毛片视频| 这个男人来自地球电影免费观看| 三级毛片av免费| 99久久国产精品久久久| 久久久久性生活片| 国产精品久久电影中文字幕| 国产成人av激情在线播放| 99久久精品一区二区三区| 熟妇人妻久久中文字幕3abv| 精品不卡国产一区二区三区| 波多野结衣高清无吗| 2021天堂中文幕一二区在线观| 日本色播在线视频| 亚州av有码| 国产国拍精品亚洲av在线观看| 舔av片在线| 女人十人毛片免费观看3o分钟| 18禁在线无遮挡免费观看视频| 欧美性感艳星| videossex国产| 在现免费观看毛片| 国产v大片淫在线免费观看| 国产av一区在线观看免费| 看片在线看免费视频| 麻豆成人av视频| 看十八女毛片水多多多| 欧美日韩一区二区视频在线观看视频在线 | 搡女人真爽免费视频火全软件| 丰满人妻一区二区三区视频av| 欧美精品国产亚洲| 国产精品精品国产色婷婷| 久久6这里有精品| 插阴视频在线观看视频| 秋霞在线观看毛片| 内地一区二区视频在线| 国产探花在线观看一区二区| 天天一区二区日本电影三级| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成年av动漫网址| a级毛片免费高清观看在线播放| 男女下面进入的视频免费午夜| 黄片wwwwww| 久久久久精品久久久久真实原创| 免费播放大片免费观看视频在线观看 | 国产乱人视频| 麻豆久久精品国产亚洲av| 欧美一区二区精品小视频在线| 一区二区三区高清视频在线| 1000部很黄的大片| 亚洲电影在线观看av| 一级毛片aaaaaa免费看小| 成年版毛片免费区| 免费播放大片免费观看视频在线观看 | 热99re8久久精品国产| 国产 一区 欧美 日韩| 久久鲁丝午夜福利片| 亚洲精品日韩在线中文字幕| 国产亚洲午夜精品一区二区久久 | 搡老妇女老女人老熟妇| 青春草国产在线视频| 精品久久久久久久人妻蜜臀av| 在线a可以看的网站| 国产精品av视频在线免费观看| 观看美女的网站| 三级毛片av免费| 看片在线看免费视频| 我的老师免费观看完整版| 亚洲av男天堂| 18禁裸乳无遮挡免费网站照片| 亚洲久久久久久中文字幕| 伊人久久精品亚洲午夜| 精品久久久久久久久亚洲| 日韩制服骚丝袜av| 欧美成人a在线观看| 蜜桃亚洲精品一区二区三区| 亚洲在久久综合| 亚洲欧美精品专区久久| 日本av手机在线免费观看| 欧美3d第一页| 99久久精品一区二区三区| 国产免费一级a男人的天堂| 热99在线观看视频| 好男人视频免费观看在线| 色尼玛亚洲综合影院| 国产白丝娇喘喷水9色精品| 久久99蜜桃精品久久| 亚洲av电影在线观看一区二区三区 | 亚洲av成人精品一二三区| 男女视频在线观看网站免费| 国产av一区在线观看免费| 一个人看的www免费观看视频| 91av网一区二区| 18禁在线播放成人免费| 超碰97精品在线观看| 日日摸夜夜添夜夜添av毛片| 村上凉子中文字幕在线| 亚洲熟妇中文字幕五十中出| 精品一区二区三区视频在线| 亚洲丝袜综合中文字幕| 国产黄片美女视频| 亚洲成人精品中文字幕电影| 最近视频中文字幕2019在线8| 嫩草影院精品99| 日本欧美国产在线视频| 久久久a久久爽久久v久久| 99久久九九国产精品国产免费| 国产三级在线视频| 色综合站精品国产| 女人被狂操c到高潮| 小蜜桃在线观看免费完整版高清| 两个人的视频大全免费| 欧美性猛交黑人性爽| 国产麻豆成人av免费视频| 久久精品熟女亚洲av麻豆精品 | 在线观看av片永久免费下载| 国产亚洲最大av| 国产人妻一区二区三区在| 美女被艹到高潮喷水动态| 大话2 男鬼变身卡| av又黄又爽大尺度在线免费看 | 国产精品人妻久久久久久| 精品一区二区三区人妻视频| 亚州av有码| 亚洲欧美日韩高清专用| 人人妻人人看人人澡| 在线观看av片永久免费下载| 久久精品影院6| 欧美xxxx黑人xx丫x性爽| 成年av动漫网址| 精品国产露脸久久av麻豆 | 国产精品日韩av在线免费观看| 精品无人区乱码1区二区| 国产成人一区二区在线| 亚洲欧美成人综合另类久久久 | 日韩欧美 国产精品| 国产一区二区亚洲精品在线观看| 国产亚洲一区二区精品| 色综合色国产| 亚洲av免费在线观看| av女优亚洲男人天堂| 国产精品电影一区二区三区| 2021天堂中文幕一二区在线观| 国产免费视频播放在线视频 | 欧美三级亚洲精品| 狂野欧美激情性xxxx在线观看| 成人二区视频| 天堂中文最新版在线下载 | 亚洲精品乱码久久久久久按摩| 美女xxoo啪啪120秒动态图| 亚洲av不卡在线观看| 欧美日韩综合久久久久久| av免费观看日本| 久久精品国产自在天天线| 国产黄a三级三级三级人| 最近中文字幕高清免费大全6| 精品一区二区三区视频在线| 亚洲图色成人| 久久鲁丝午夜福利片| 国产精品综合久久久久久久免费| 日本wwww免费看| 狂野欧美激情性xxxx在线观看| 成人无遮挡网站| 久久久久久久国产电影| www.av在线官网国产| 国产成年人精品一区二区| 国产黄色视频一区二区在线观看 | 夜夜爽夜夜爽视频| 只有这里有精品99| 国语对白做爰xxxⅹ性视频网站| 欧美潮喷喷水| 久久婷婷人人爽人人干人人爱| 亚洲欧美中文字幕日韩二区| 成人高潮视频无遮挡免费网站| 我的女老师完整版在线观看| 麻豆精品久久久久久蜜桃| 国产黄片美女视频| 久久久久久国产a免费观看| 国产一级毛片在线| 成人毛片a级毛片在线播放| 久久精品久久精品一区二区三区| 国产探花在线观看一区二区| 黄片无遮挡物在线观看| 国产免费视频播放在线视频 | 日本黄色视频三级网站网址| or卡值多少钱| 干丝袜人妻中文字幕| a级一级毛片免费在线观看| 亚洲av日韩在线播放| 有码 亚洲区| 中文字幕免费在线视频6| 国内精品一区二区在线观看| 搡女人真爽免费视频火全软件| 91久久精品国产一区二区三区| 亚洲av电影不卡..在线观看| 国产精品一区二区在线观看99 | 亚洲成人中文字幕在线播放| 天堂√8在线中文| 老司机影院成人| 天堂网av新在线| 中文欧美无线码| 欧美一区二区亚洲| 亚洲人成网站在线播| 欧美一区二区精品小视频在线| 亚洲欧洲日产国产| 水蜜桃什么品种好| 国产黄片美女视频| 国产私拍福利视频在线观看| 日本欧美国产在线视频| 欧美日本视频| 午夜a级毛片| 欧美xxxx性猛交bbbb| 国产成人午夜福利电影在线观看| 欧美一区二区精品小视频在线| 亚洲欧美日韩卡通动漫| 尤物成人国产欧美一区二区三区| 天天躁夜夜躁狠狠久久av| 亚洲最大成人av| 国产视频内射| 在线天堂最新版资源| 又粗又硬又长又爽又黄的视频| 国产老妇伦熟女老妇高清| 国产成人freesex在线| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美中文字幕日韩二区| 欧美激情国产日韩精品一区| 久久亚洲精品不卡| 热99re8久久精品国产| 国产午夜精品论理片| 美女黄网站色视频| 非洲黑人性xxxx精品又粗又长| 成人亚洲欧美一区二区av| 色吧在线观看| 国内精品宾馆在线| 国产高清不卡午夜福利| 国产黄a三级三级三级人| 青春草国产在线视频| 欧美日本视频| 国产精品日韩av在线免费观看| 三级国产精品片| 寂寞人妻少妇视频99o| 97在线视频观看| 免费在线观看成人毛片| 精品国产露脸久久av麻豆 | 亚洲美女视频黄频| 国产精品一区www在线观看| 国产精品国产三级专区第一集| 天美传媒精品一区二区| 99久久人妻综合| 熟女电影av网| av视频在线观看入口| 男女那种视频在线观看| 日本色播在线视频| 国产精品久久视频播放| 亚洲人成网站高清观看| 搡女人真爽免费视频火全软件| 午夜福利视频1000在线观看| 亚洲欧美日韩东京热| 亚洲国产精品专区欧美| 日韩亚洲欧美综合| 精品酒店卫生间| 免费不卡的大黄色大毛片视频在线观看 | 狠狠狠狠99中文字幕| 国产精品永久免费网站| 亚洲成人av在线免费| 欧美一区二区亚洲| 国产综合懂色| 色噜噜av男人的天堂激情| 一区二区三区免费毛片| 国产精品野战在线观看| 国产精品,欧美在线| 国产淫片久久久久久久久| 亚洲av二区三区四区| 国产成年人精品一区二区| 国产精品蜜桃在线观看| 成人毛片60女人毛片免费| 九九热线精品视视频播放| 亚洲在线观看片| 日本免费一区二区三区高清不卡| 精品久久久久久久久av| 国产精品久久久久久久久免| 免费看a级黄色片| 日韩av在线大香蕉| 麻豆成人av视频| 亚洲乱码一区二区免费版| 少妇的逼好多水| 免费观看的影片在线观看| 99热全是精品| 国产成人午夜福利电影在线观看| 免费不卡的大黄色大毛片视频在线观看 | 一区二区三区免费毛片| 九草在线视频观看| 国产黄a三级三级三级人| 丝袜喷水一区| 欧美日本亚洲视频在线播放| 简卡轻食公司| 亚洲国产精品合色在线| 黑人高潮一二区| 免费黄色在线免费观看| 91狼人影院| 久久久精品94久久精品| 欧美极品一区二区三区四区| 深夜a级毛片| 亚洲高清免费不卡视频| 99热精品在线国产| 精品国产露脸久久av麻豆 | 午夜爱爱视频在线播放| 永久免费av网站大全| 午夜老司机福利剧场| 日韩欧美在线乱码| 亚洲国产欧美人成| 国产91av在线免费观看| 日本午夜av视频| 大香蕉97超碰在线| 国产白丝娇喘喷水9色精品| 久久久久久国产a免费观看| 成人综合一区亚洲| eeuss影院久久| 亚洲人成网站在线观看播放| 2021天堂中文幕一二区在线观| 少妇人妻精品综合一区二区| 神马国产精品三级电影在线观看| 激情 狠狠 欧美| 国产欧美日韩精品一区二区| 日韩三级伦理在线观看| 亚洲av电影不卡..在线观看| 精品久久久久久久久亚洲| 午夜福利视频1000在线观看| 欧美区成人在线视频| 久久久精品94久久精品| 日韩欧美国产在线观看| 亚洲美女视频黄频| 少妇人妻精品综合一区二区| 边亲边吃奶的免费视频| 91精品伊人久久大香线蕉| 免费看光身美女| 亚洲av成人av| 丝袜喷水一区| 哪个播放器可以免费观看大片| 国产亚洲一区二区精品| 精品人妻视频免费看| 日本一二三区视频观看| 亚洲国产高清在线一区二区三| 日韩强制内射视频| 日韩在线高清观看一区二区三区| 国产精品国产三级国产av玫瑰| 久久久久久久国产电影| eeuss影院久久| 乱系列少妇在线播放| 国产成人精品久久久久久| 一区二区三区四区激情视频| 麻豆成人av视频| 99在线人妻在线中文字幕| 亚洲在线观看片| 看非洲黑人一级黄片| 性色avwww在线观看| 欧美激情久久久久久爽电影| 久久精品综合一区二区三区| 啦啦啦韩国在线观看视频| 国产午夜精品久久久久久一区二区三区| 丝袜美腿在线中文| 亚洲欧美成人精品一区二区| 欧美性猛交黑人性爽| 别揉我奶头 嗯啊视频| 久久久久精品久久久久真实原创| 午夜爱爱视频在线播放| 中文字幕av成人在线电影| 国产成人a∨麻豆精品| 国产午夜精品论理片| 国产精品人妻久久久影院| 内地一区二区视频在线| 亚洲乱码一区二区免费版| 国产黄色小视频在线观看| 啦啦啦韩国在线观看视频| av在线观看视频网站免费| 能在线免费观看的黄片| 黄色日韩在线| 噜噜噜噜噜久久久久久91| 99热这里只有精品一区| 自拍偷自拍亚洲精品老妇| 性色avwww在线观看| 国产精品不卡视频一区二区| 国产精品人妻久久久影院| 亚洲精品456在线播放app| 少妇人妻一区二区三区视频| 人妻系列 视频| 久久精品91蜜桃| 丰满少妇做爰视频| 久久久久久久久久成人| 偷拍熟女少妇极品色| 乱码一卡2卡4卡精品| 久久99蜜桃精品久久| 亚洲欧美日韩高清专用| 色网站视频免费| 视频中文字幕在线观看| 日本一本二区三区精品| 欧美精品国产亚洲| 免费黄网站久久成人精品| 我的老师免费观看完整版| 久久久国产成人免费| 人妻制服诱惑在线中文字幕| 色综合站精品国产| 特级一级黄色大片| 欧美精品国产亚洲| 亚洲av成人av| 最近最新中文字幕大全电影3| 免费电影在线观看免费观看| 久久久精品大字幕| 色哟哟·www| 精品久久久久久久久久久久久| 久久这里只有精品中国| 国产大屁股一区二区在线视频| 黄片无遮挡物在线观看| 午夜福利成人在线免费观看| 国产极品天堂在线| 欧美3d第一页| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 日本免费一区二区三区高清不卡| 日本一二三区视频观看| 五月伊人婷婷丁香| 毛片女人毛片| 九草在线视频观看| 最后的刺客免费高清国语| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品一及| 国产乱人偷精品视频| 精品人妻偷拍中文字幕| 国产精品久久久久久av不卡| 国产精品电影一区二区三区| 日本猛色少妇xxxxx猛交久久| 老女人水多毛片| 久久这里只有精品中国| 国产免费又黄又爽又色| 国产极品天堂在线| 亚洲av熟女| 久久精品国产亚洲av天美| 三级经典国产精品| 亚洲精品,欧美精品| 最近中文字幕高清免费大全6| 51国产日韩欧美| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品合色在线| 高清日韩中文字幕在线| 欧美色视频一区免费| 日韩欧美在线乱码| 久久国产乱子免费精品| 一区二区三区乱码不卡18| 精品久久久久久久人妻蜜臀av| 久久婷婷人人爽人人干人人爱| 日本一二三区视频观看| 亚洲欧美日韩东京热| 亚洲精品影视一区二区三区av| 欧美变态另类bdsm刘玥| 综合色丁香网| 又黄又爽又刺激的免费视频.| 一级毛片电影观看 | 纵有疾风起免费观看全集完整版 | 国产精品无大码| 中文字幕亚洲精品专区| 国产三级在线视频| 国产片特级美女逼逼视频| 亚洲av二区三区四区| 99在线视频只有这里精品首页| 国产av码专区亚洲av| 少妇猛男粗大的猛烈进出视频 | 国产黄色视频一区二区在线观看 | 国产精品久久久久久精品电影小说 | 神马国产精品三级电影在线观看| 国产精品一区www在线观看| 国产一区二区三区av在线| 精品无人区乱码1区二区| ponron亚洲| 如何舔出高潮| 欧美日本亚洲视频在线播放| 91精品伊人久久大香线蕉| 国产免费一级a男人的天堂| 午夜免费激情av| 亚洲性久久影院| 欧美zozozo另类| 国产免费男女视频| 禁无遮挡网站| 九九在线视频观看精品| 少妇人妻精品综合一区二区| 少妇被粗大猛烈的视频| 中文字幕免费在线视频6| 国产男人的电影天堂91| 欧美变态另类bdsm刘玥| 国产成人福利小说| 人妻制服诱惑在线中文字幕| 一个人免费在线观看电影| 中文亚洲av片在线观看爽| 久久韩国三级中文字幕| 久久热精品热| 2021少妇久久久久久久久久久| 亚洲精品,欧美精品| 亚洲高清免费不卡视频| 又粗又爽又猛毛片免费看| 日本黄色片子视频| 久久99蜜桃精品久久| 亚洲精品一区蜜桃| 日本午夜av视频| 一卡2卡三卡四卡精品乱码亚洲| 蜜臀久久99精品久久宅男| 久久精品国产亚洲网站| 日韩欧美精品v在线| 中文字幕人妻熟人妻熟丝袜美| 美女国产视频在线观看| 国产精品不卡视频一区二区| 听说在线观看完整版免费高清| 亚洲成av人片在线播放无| a级毛片免费高清观看在线播放| 日韩在线高清观看一区二区三区| 一二三四中文在线观看免费高清| 三级男女做爰猛烈吃奶摸视频| av视频在线观看入口| 午夜福利网站1000一区二区三区| 中文乱码字字幕精品一区二区三区 | 青春草视频在线免费观看| 99在线人妻在线中文字幕| 亚洲国产精品合色在线| 欧美色视频一区免费| 亚洲精品aⅴ在线观看| 国产精品久久久久久久久免| 搡女人真爽免费视频火全软件| 亚洲av中文字字幕乱码综合| 成年女人永久免费观看视频| 九草在线视频观看| 六月丁香七月| 久久欧美精品欧美久久欧美| 啦啦啦啦在线视频资源| 国产高清视频在线观看网站| 国产精品麻豆人妻色哟哟久久 | 国产免费视频播放在线视频 | 三级国产精品欧美在线观看| 亚洲精品456在线播放app| 五月伊人婷婷丁香| 国产精品野战在线观看| 久久久久久久午夜电影| 一级黄色大片毛片| 久久精品国产亚洲网站| 婷婷色综合大香蕉| 国产伦精品一区二区三区四那| 国产视频首页在线观看| 国产av码专区亚洲av| 一本一本综合久久| 在线观看一区二区三区| 在线播放国产精品三级| 免费大片18禁| 内射极品少妇av片p| av女优亚洲男人天堂| 26uuu在线亚洲综合色| 国产av一区在线观看免费| 欧美性感艳星| www.av在线官网国产| 国产成人aa在线观看| 国产av在哪里看| 久久久久久久亚洲中文字幕| 国产精品三级大全| 亚洲成人av在线免费| 18禁裸乳无遮挡免费网站照片| 国产午夜精品久久久久久一区二区三区| 午夜福利视频1000在线观看| 人妻夜夜爽99麻豆av| 国产精品一区二区性色av| videossex国产| 内射极品少妇av片p| 成人漫画全彩无遮挡| 亚洲精品一区蜜桃| 亚洲欧美成人综合另类久久久 | 色尼玛亚洲综合影院| 中文欧美无线码| 国语对白做爰xxxⅹ性视频网站| 国产免费一级a男人的天堂| 国产精品久久久久久av不卡| 啦啦啦韩国在线观看视频| a级毛片免费高清观看在线播放| 国产又黄又爽又无遮挡在线| 亚洲成人中文字幕在线播放| 国产成人aa在线观看| av在线蜜桃| 成人漫画全彩无遮挡| 在线观看美女被高潮喷水网站| 国产综合懂色| 国产不卡一卡二| 熟女人妻精品中文字幕| 国产乱人视频| 99热这里只有是精品在线观看| 看十八女毛片水多多多| 日日干狠狠操夜夜爽| 日韩成人伦理影院| 丰满乱子伦码专区| 日韩在线高清观看一区二区三区| 在线免费观看的www视频| 国产v大片淫在线免费观看| 国产三级在线视频| 午夜精品一区二区三区免费看| 国产精品三级大全| 在线播放无遮挡| 午夜精品在线福利| 久久久亚洲精品成人影院| 尤物成人国产欧美一区二区三区|