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

    三維空間熱路模型及其在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中的應(yīng)用

    2015-02-24 00:46:19田井呈王飛宇卓克瓊趙朝會(huì)
    關(guān)鍵詞:數(shù)學(xué)模型溫度場(chǎng)電動(dòng)機(jī)

    田井呈, 王飛宇, 卓克瓊, 趙朝會(huì)

    (上海電機(jī)學(xué)院 電氣學(xué)院, 上海 201306)

    ?

    三維空間熱路模型及其在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中的應(yīng)用

    田井呈,王飛宇,卓克瓊,趙朝會(huì)

    (上海電機(jī)學(xué)院 電氣學(xué)院, 上海 201306)

    摘要:基于等效熱路模型的溫度場(chǎng)計(jì)算方法研究受到了廣泛的關(guān)注,針對(duì)運(yùn)用于電動(dòng)機(jī)溫度場(chǎng)計(jì)算的三維空間等效熱路模型,闡述了熱路模型的基本原理,總結(jié)了立方體T型、空心圓筒雙T型、三維立方體和扇體熱路的數(shù)學(xué)模型。結(jié)合熱路模型在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中的應(yīng)用,探討了其未來(lái)的發(fā)展趨勢(shì),對(duì)電動(dòng)機(jī)溫度場(chǎng)的計(jì)算及其散熱系統(tǒng)的設(shè)計(jì)具有一定的借鑒意義和參考價(jià)值。

    關(guān)鍵詞:電動(dòng)機(jī); 溫度場(chǎng); 三維空間等效熱路模型; 數(shù)學(xué)模型

    電動(dòng)機(jī)溫度場(chǎng)的計(jì)算方法主要有熱參數(shù)法、等效熱路法、計(jì)算機(jī)數(shù)值解法。其中,熱參數(shù)法是采用疊加原理計(jì)算各熱源單獨(dú)存在時(shí)被研究對(duì)象的溫升之和。等效熱路法假設(shè)研究對(duì)象中真實(shí)熱源和熱阻被少量的集中熱源和等值熱阻所代替,并假定兩者不取決于熱流的大小,通過(guò)繪制電動(dòng)機(jī)的等效熱路圖并結(jié)合電路的求解方法計(jì)算電動(dòng)機(jī)中的溫度分布。計(jì)算機(jī)數(shù)值解法主要有兩類,一類是以等效熱路法為基礎(chǔ)的熱網(wǎng)絡(luò)法[2-4];另一類是以計(jì)算傳熱學(xué)為基礎(chǔ)的數(shù)值計(jì)算方法[5-6],如計(jì)算流體動(dòng)力學(xué)(Computational Fluid Dynamics, CFD)和有限元分析(Finite Element Analysis, FEA)軟件。實(shí)際工程中,溫度的三維分布給其求解和分析都會(huì)帶來(lái)很大困難;而熱路法的運(yùn)用可避開三維溫度場(chǎng)直接求解的復(fù)雜性,其利用熱學(xué)中較為成熟的解析計(jì)算方法,將三維問題轉(zhuǎn)換為一維問題,結(jié)合線性電路方程的求解方法進(jìn)行求解,一直以來(lái)受到了國(guó)內(nèi)外學(xué)者的普遍關(guān)注[7-9]。為了提高熱路法的求解精度,進(jìn)一步考慮電動(dòng)機(jī)中所求區(qū)域材料導(dǎo)熱性質(zhì)的各向異性,基于三維空間熱路法的電動(dòng)機(jī)溫度場(chǎng)計(jì)算方法已成為近年來(lái)的研究熱點(diǎn)。

    文獻(xiàn)[10]中基于立方體和空心圓筒T型熱路單元,主要研究了不同轉(zhuǎn)速下旋轉(zhuǎn)電動(dòng)機(jī)氣隙的對(duì)流換熱系數(shù),建立了適用于電動(dòng)機(jī)溫度場(chǎng)計(jì)算的三維熱路模型,通過(guò)對(duì)3種不同結(jié)構(gòu)的電動(dòng)機(jī)進(jìn)行測(cè)試,證明了該模型的準(zhǔn)確性,為高功率密度電動(dòng)機(jī)的設(shè)計(jì)和應(yīng)用提供了指導(dǎo)。

    為提高T型熱路法的計(jì)算精度,文獻(xiàn)[11]中詳細(xì)推導(dǎo)了T型熱路法的基本原理,進(jìn)一步修正了T型熱路法的不足。文獻(xiàn)[12]中提出了立方體熱路單元,并將其運(yùn)用于感應(yīng)器的溫度場(chǎng)計(jì)算,利用有限元法驗(yàn)證了該模型的正確性??紤]到旋轉(zhuǎn)電動(dòng)機(jī)的結(jié)構(gòu)特點(diǎn),文獻(xiàn)[13]中提出了扇體熱路單元的概念,并檢驗(yàn)了該模型的正確性。文獻(xiàn)[14]中充分吸收了立方體和扇體熱路單元的特點(diǎn),提出了三維空間熱路法的概念,并將其運(yùn)用于一臺(tái)開關(guān)磁阻電動(dòng)機(jī)的溫度場(chǎng)計(jì)算,其結(jié)果與基于有限元法的計(jì)算結(jié)果相比較相差甚微,滿足求解精度的要求。

    由上述研究可見,三維空間熱路法主要包括立方體熱路模型和空心圓筒熱路模型。其中,由于等效條件不同,立方體熱路模型又分為僅沿直角坐標(biāo)系某軸向的單一3節(jié)點(diǎn)T型熱路模型和多軸向多節(jié)點(diǎn)的立方體熱路模型。除雙T型熱路模型外,空心圓筒熱路模型可進(jìn)一步等效為扇體熱路模型。由于三維熱路模型可針對(duì)某一指定的求解區(qū)域,根據(jù)其結(jié)構(gòu)特點(diǎn)和傳熱特性,對(duì)其利用不同的熱路模型單元進(jìn)行劃分并分析計(jì)算,故對(duì)于電動(dòng)機(jī)這種結(jié)構(gòu)復(fù)雜、多物理域相互作用的系統(tǒng),三維熱路模型的建模方式具有很大的優(yōu)勢(shì)。

    本文在詳細(xì)闡述等效熱路法基本原理的基礎(chǔ)上,總結(jié)了立方體T型熱路、空心圓筒雙T型熱路、立方體和扇體熱路的數(shù)學(xué)模型;并結(jié)合三維熱路模型在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中的應(yīng)用,探討了其未來(lái)的發(fā)展趨勢(shì),對(duì)電動(dòng)機(jī)的溫度場(chǎng)計(jì)算及其散熱系統(tǒng)的設(shè)計(jì)具有一定的參考價(jià)值。

    1熱路的基本原理及特點(diǎn)

    1.1 熱路法的基本原理

    溫度場(chǎng)的求解問題實(shí)際上是一個(gè)三維問題。等效熱路法的運(yùn)用是為了避開復(fù)雜的三維問題,轉(zhuǎn)而將其簡(jiǎn)化為一維問題求解。如圖1所示,假設(shè)熱量只沿x軸方向傳導(dǎo),立方體內(nèi)部無(wú)熱量生成且沿x軸方向上兩個(gè)平面的溫度分別為T1和T2,則可得熱量沿x軸方向上的等效熱路。該熱路具有2個(gè)端子,熱阻

    (1)

    式中,A為立方體沿x軸方向上的表面積;K為與立方體材料特性有關(guān)的等效導(dǎo)熱系數(shù);L為熱流路徑的長(zhǎng)度。由等效熱路可知,

    T1-T2=qR

    (2)

    式中,q為流入立方體的熱流量,類似于電路中的電流。

    圖1 一維傳熱模型Fig.1 One dimensional heat transfer model

    以上為熱路法的基本原理。為滿足實(shí)際工程運(yùn)用要求,等效熱路法的應(yīng)用須滿足如下條件:

    (1) 忽略所研究物體內(nèi)部的熱物理過(guò)程;

    (2) 物體本身所生產(chǎn)的熱量在物體中均勻分布;

    (3) 熱量在物體中的流動(dòng)是單向且獨(dú)立的,忽略熱量的回流;

    (4) 材料的導(dǎo)熱系數(shù)是線性的,忽略導(dǎo)熱性能的非線性影響。

    1.2 熱路法的應(yīng)用步驟

    熱路法來(lái)源于集中參數(shù)電路理論,利用傳熱學(xué)中較為成熟的解析法將所求區(qū)域集中為相應(yīng)的熱阻、熱容等參數(shù)元件,即只考慮物體外部的溫度效應(yīng)而忽略物體內(nèi)部發(fā)生的細(xì)微過(guò)程。熱路與電路的定義對(duì)比如表1所示。其中,熱容的引入是為了表征材料本身儲(chǔ)存熱量的能力,即溫度每變化1℃材料吸收或放出熱量,用于瞬態(tài)溫度計(jì)算。熱容

    Cp=c ρ V

    (3)

    式中,c為材料比熱容;ρ為材料密度;V為所分析物體的體積。

    表1 熱路與電路符號(hào)定義對(duì)比

    針對(duì)某一特定的求解區(qū)域,一般將熱路法的應(yīng)用分為以下步驟: ① 確定求解區(qū)域;② 根據(jù)求解區(qū)域選擇熱路單元;③ 計(jì)算熱路參數(shù);④ 編輯熱路圖;⑤ 求解熱路圖;⑥ 分析計(jì)算結(jié)果。

    熱路單元的選擇在整個(gè)熱路法運(yùn)用過(guò)程中尤其重要,對(duì)求解區(qū)域合理的單元?jiǎng)澐种苯記Q定了最終計(jì)算結(jié)果的精度。本文詳細(xì)敘述目前運(yùn)用的熱路單元及其數(shù)學(xué)模型。

    2三維熱路單元及其數(shù)學(xué)模型

    2.1 立方體T型熱路單元

    物體的溫度升高往往是由于物體內(nèi)部生成的熱量造成的。圖1中的立方體模型沒有考慮物體內(nèi)部產(chǎn)生的熱量,僅作為熱量傳遞的路徑。若計(jì)及物體內(nèi)部的生熱并假設(shè)其內(nèi)部生成的熱流為q,立方體的溫度變化只與其本身所產(chǎn)生的熱量相關(guān),且熱量只沿x軸方向傳導(dǎo),忽略外部環(huán)境溫度的影響,則將其作如圖2所示的等效。圖2(b)中,T為該立方體的平均溫度值。需要說(shuō)明的是,T與T′之間的熱阻為負(fù)值,為等效后的附加熱阻,若省略,將使得求解溫度值偏大[11]。通過(guò)等效熱路圖,再結(jié)合一定的邊界條件,就能估算整個(gè)立方體的平均溫度值以及沿x軸方向上兩個(gè)表面的平均溫度值。

    圖2 立方體T型熱路模型Fig.2 T-type cuboidal thermal model

    2.2 立方體三維空間熱路單元

    立方體T型熱路單元只計(jì)及單一軸向上的溫度等效計(jì)算,但是材料的導(dǎo)熱性能是各向異性的,導(dǎo)致物體中溫度的分布梯度也是各向異性的。為提高求解的精度,可進(jìn)一步考慮熱量在立方體各方向的流動(dòng),從而產(chǎn)生了如圖3所示的三維度空間熱路模型。

    圖3 立方體熱路單元Fig.3 3D cuboidal thermal model

    為區(qū)別于立方體T型熱路單元,將其稱為立方體三維空間熱路單元。圖3中:

    (4)

    (5)

    (6)

    (7)

    Cp=cρV

    (8)

    式中,Kx、Ky、Kz分別為沿x、y、z方向的傳熱系數(shù),取決于材料本身熱物理性能;Ax、Ay、Az分別為沿x、y、z方向的有效傳熱面積;Lx、Ly、Lz分別為沿x、y、z方向的傳熱路徑長(zhǎng)度;此外,Tx1、Tx2,Ty1、Ty2,Tz1、Tz2分別為立方體沿x、y、z方向上各個(gè)表面的平均溫度值;TC為整個(gè)立方體單元的平均溫度值;RC為前述引入的附加熱阻。

    2.3 空心圓筒雙T型熱路單元

    圖4所示為空心圓筒雙T型熱路模型。

    圖4 空心圓筒熱路模型Fig.4 Hollow cylinder thermal model

    圖4中,T1、T2分別為空心圓筒內(nèi)外兩個(gè)圓表面的平均溫度值,T3、T4分別為圓筒上、下兩個(gè)底面的平均溫度值。同理,TH為整個(gè)空心圓筒的平均溫度值,Ra為軸向等效傳熱路徑上的附加熱阻,Rr為徑向等效傳熱路徑上的等效熱阻,RH為整個(gè)熱路模型的附加熱阻。其中:

    RH=Ra+Rr

    式中,ka、kr分別為空心圓筒沿軸向和徑向的等效導(dǎo)熱系數(shù)。

    空心圓筒熱路單元也是一種典型的熱路單元,來(lái)源于多層圓筒壁導(dǎo)熱的解析計(jì)算,其應(yīng)用需要滿足如下假設(shè)條件:

    (1) 沿圓筒軸向和徑向的熱流路徑是獨(dú)立的;

    (2) 圓筒溫度等效為平均溫度,表征圓筒軸向和徑向的熱流流動(dòng);

    (3) 熱流是單向傳遞的,同一路徑上無(wú)回流;

    (4) 圓筒中產(chǎn)生的熱量在圓筒中均勻分布。

    2.4 扇體熱路單元

    同立方體熱路單元的發(fā)展一樣,為提高空心圓筒熱路單元的計(jì)算精度,進(jìn)一步考慮圓筒沿周向上的熱量流動(dòng),產(chǎn)生了圖5所示的扇體熱路單元。

    圖5 扇體熱路單元Fig.5 Arc segment element

    圖中:

    式中,Ka、Kc、Kr分別為沿軸向、周向、徑向的等效導(dǎo)熱系數(shù);Ta1、Ta2,Tc1、Tc2,Tr1、Tr2分別為扇體沿軸向、周向和徑向上各表面的平均溫度值;TA為扇體平均溫度值;RA為附加熱阻。同立方體單元一樣,利用等效熱路和相應(yīng)的邊界條件可估算該扇體不同表面的溫度大小和平均溫度值。

    3空間熱路法在電動(dòng)機(jī)溫度場(chǎng)分析中的應(yīng)用

    3.1 基于空間熱路法的電動(dòng)機(jī)溫度場(chǎng)計(jì)算

    1979年,Perez等[15]提出了運(yùn)用于電動(dòng)機(jī)溫度場(chǎng)計(jì)算的穩(wěn)態(tài)等效熱路圖,并對(duì)T形等效熱路單元的數(shù)學(xué)模型進(jìn)行了詳細(xì)推導(dǎo)。通過(guò)對(duì)一系列電動(dòng)機(jī)溫度測(cè)試結(jié)果的分析,證明了該建模方式在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中的可行性,但并未對(duì)電動(dòng)機(jī)軸承等效熱阻的計(jì)算提供相應(yīng)的方法。隨后,Mellor等[16]在電動(dòng)機(jī)發(fā)熱元件對(duì)應(yīng)的熱路溫度節(jié)點(diǎn)中引入熱容,計(jì)算元件吸熱和放熱瞬態(tài)過(guò)程對(duì)電動(dòng)機(jī)元件中溫度分布的影響,建立了較為完善的中小型電動(dòng)機(jī)十節(jié)點(diǎn)瞬態(tài)等效熱路模型;并通過(guò)對(duì)不同廠商設(shè)計(jì)的一系列全封閉自扇冷式感應(yīng)電動(dòng)機(jī)進(jìn)行測(cè)試,結(jié)果表明該模型能夠準(zhǔn)確預(yù)測(cè)電動(dòng)機(jī)主要部位的溫度,不足之處在于熱路中溫度節(jié)點(diǎn)數(shù)目增多,使得熱路中相關(guān)參數(shù)的計(jì)算過(guò)程變得繁瑣。

    基于立方體T型和空心圓筒雙T型熱路單元兩種基本的空間熱路單元的電動(dòng)機(jī)溫度場(chǎng)計(jì)算分析,國(guó)內(nèi)外學(xué)者做了大量研究[17-19]。但對(duì)于更高精度的三維立方體和扇體熱路單元的電動(dòng)機(jī)溫度場(chǎng)計(jì)算方法,是近年來(lái)才出現(xiàn)的。圖6所示的球形電動(dòng)機(jī)溫度熱路模型就是典型的三維空間熱路模型[20],圖中所用符號(hào)與本文所述不盡相同,相關(guān)的熱路參數(shù)符號(hào)分為4部分,即損耗(熱源)、等效傳導(dǎo)熱阻、等效對(duì)流換熱熱阻和等效輻射換熱熱阻。損耗分別為永磁體渦流損耗Pm、Pm1、Pm2,轉(zhuǎn)子鐵芯損耗PAl1、PAl2,定子鐵芯損耗PAl3,繞組銅耗PCu,線圈與鐵芯的損耗PCuAl1、PCuAl2、PCuAl3、PCuAl4,定子外殼產(chǎn)生的損耗PAl4。

    圖6 球形電動(dòng)機(jī)三維熱路模型Fig.6 Spherical motor 3D thermal model

    整個(gè)熱阻中的傳導(dǎo)熱阻分別為永磁體等效傳導(dǎo)熱阻Rm1、Rm2,轉(zhuǎn)子鐵芯和永磁體間的等效傳導(dǎo)熱阻RAl1、RAl2,定子鐵芯與線圈間的等效傳導(dǎo)熱阻R01、R02、R03和R04,線圈內(nèi)部的等效傳導(dǎo)熱阻R11、R12、R21、R22、R31、R32、R41、R42,定子鐵芯內(nèi)部的等效傳導(dǎo)熱阻R51、R52,定子鐵芯與定子外殼間的等效傳導(dǎo)熱阻RAlAl,線圈與定子外殼間的等效傳導(dǎo)熱阻RCuAl3、RCuAl4、RCuAl5、RCuAl6。

    對(duì)流換熱熱阻包括永磁體與轉(zhuǎn)子腔內(nèi)氣體間的等效對(duì)流換熱熱阻RmA,定子鐵芯與永磁體間的等效對(duì)流換熱熱阻RAlm,線圈與永磁體的等效對(duì)流換熱熱阻RCum1、RCum2、RCum3,線圈與轉(zhuǎn)子鐵芯間的等效對(duì)流換熱熱阻RCuAl1和RCuAl2,線圈間的等效對(duì)流換熱熱阻RCu1、RCu2、RCu3、RCu4,定子外殼與外部氣體間的等效對(duì)流換熱熱阻RAlA1。

    最后,計(jì)及定子外殼與外部氣體間發(fā)生的輻射換熱,引入等效輻射換熱熱阻RA1A2。文獻(xiàn)[20]中利用空間熱路模型建立了球形電動(dòng)機(jī)較為完善的等效熱路圖,并利用有限元數(shù)值解法驗(yàn)證了熱路求解結(jié)果,證明了空間熱路模型不僅能夠滿足求解精度的要求,且在很大程度上節(jié)省了計(jì)算機(jī)資源和求解時(shí)間。盡管立方體T型熱路單元和空心圓筒雙T型熱路單元已大量應(yīng)用于電動(dòng)機(jī)溫度場(chǎng)的計(jì)算,但如前所述,三維立方體和扇體熱路單元在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中的應(yīng)用較少,故文中對(duì)球形電動(dòng)機(jī)熱路模型的應(yīng)用,不僅對(duì)電動(dòng)機(jī)溫度場(chǎng)的求解具有參考價(jià)值,且對(duì)于三維立方體和扇體熱路單元的發(fā)展和應(yīng)用均具有重要的影響。

    3.2 電動(dòng)機(jī)空間熱路計(jì)算軟件Motor-CAD的應(yīng)用

    伴隨著空間熱路法在電動(dòng)機(jī)中的運(yùn)用,相應(yīng)的電動(dòng)機(jī)熱路計(jì)算軟件應(yīng)運(yùn)而生,Motor-CAD是唯一一款基于三維熱路模型制作而成的商用電動(dòng)機(jī)熱路分析軟件,主要基于T型熱路模型和空心圓筒熱路模型。該軟件給出了9種常用電動(dòng)機(jī)的熱路和參數(shù),這些參數(shù)來(lái)自于國(guó)內(nèi)外有關(guān)電動(dòng)機(jī)熱路計(jì)算方法的研究成果和大量的實(shí)際測(cè)試數(shù)據(jù)[21],計(jì)及了電樞繞組的處理、鐵芯的裝配工藝以及多種冷卻方式對(duì)電動(dòng)機(jī)溫度場(chǎng)計(jì)算的影響。由于軟件中集成的熱路僅與電動(dòng)機(jī)的制造材料和相關(guān)設(shè)計(jì)參數(shù)有關(guān),故能夠方便地對(duì)電動(dòng)機(jī)的溫度進(jìn)行參數(shù)化分析和優(yōu)化設(shè)計(jì),得到了廣泛的應(yīng)用[22-23]。

    圖7給出了Motor-CAD軟件界面。其中圖 7(a)(b)為基于Motor-CAD軟件的一臺(tái)表貼式永磁同步電動(dòng)機(jī)的徑向和軸向截面圖,在界面的左側(cè)輸入電動(dòng)機(jī)的相關(guān)結(jié)構(gòu)參數(shù),并選擇相應(yīng)的冷卻方式,可直接計(jì)算出電動(dòng)機(jī)中的溫度分布。其中,穩(wěn)態(tài)溫度分布以等效熱路圖的方式顯示(見圖7(c))。與電動(dòng)機(jī)的對(duì)稱結(jié)構(gòu)對(duì)應(yīng),等效熱路圖主要分為轉(zhuǎn)軸、轉(zhuǎn)子鐵芯、永磁磁鋼、電樞繞組、定子鐵芯、機(jī)殼和端蓋等部分,每一部分具有相應(yīng)的熱路拓?fù)浜蜔崧穮?shù);各部分熱路綜合起來(lái)構(gòu)成整個(gè)電動(dòng)機(jī)的熱路系統(tǒng);外圍部分支路為外部環(huán)境的溫度。基于上述熱路數(shù)學(xué)模型的分析可知,熱路中每一個(gè)節(jié)點(diǎn)代表電動(dòng)機(jī)相應(yīng)部位的溫度值,通過(guò)輸入電動(dòng)機(jī)的相關(guān)設(shè)計(jì)參數(shù),如定子外圓半徑、鐵芯長(zhǎng)度、定轉(zhuǎn)子槽數(shù)、各部分損耗值等,并選擇相應(yīng)的冷卻方式,熱路圖中即顯示相應(yīng)節(jié)點(diǎn)的溫度值、節(jié)點(diǎn)之間的熱阻值、損耗值等。此外,若選擇計(jì)及材料的瞬態(tài)儲(chǔ)熱特性,則可計(jì)算出熱路中各節(jié)點(diǎn)的瞬態(tài)溫度曲線(見圖7(d)),并可對(duì)電動(dòng)機(jī)的設(shè)計(jì)參數(shù)進(jìn)行溫度敏感性分析,從而可以優(yōu)化電動(dòng)機(jī)結(jié)構(gòu)。

    圖7 Motor-CAD軟件界面

    3.3 空間熱路法的發(fā)展趨勢(shì)

    從T型熱路模型、空心圓筒熱路模型到立方體和扇體熱路模型的發(fā)展來(lái)看,空間熱路法的等效程度在逐漸提高,其模型也經(jīng)過(guò)多次檢驗(yàn)和修正。但在電動(dòng)機(jī)溫度場(chǎng)計(jì)算中,大多還集中在對(duì)立方體T型熱路單元以及空心圓筒熱路單元的應(yīng)用上,而高精度的三維立方體和扇體熱路單元的運(yùn)用還較少見。盡管近年來(lái)有限元數(shù)值解法在溫度場(chǎng)計(jì)算中的運(yùn)用越來(lái)越成熟,但仍然無(wú)法對(duì)電動(dòng)機(jī)進(jìn)行全域溫度分析求解。對(duì)于電動(dòng)機(jī)這樣一個(gè)復(fù)雜的系統(tǒng)而說(shuō),合理地選用相應(yīng)的熱路單元,在一定的等效條件下,三維空間熱路法的運(yùn)用將具有更大的優(yōu)勢(shì),可作為電動(dòng)機(jī)全域溫度場(chǎng)分析與計(jì)算的有力途徑。未來(lái)對(duì)于空間熱路法的發(fā)展可歸納為以下幾個(gè)方面:

    (1) 復(fù)雜電動(dòng)機(jī)運(yùn)行系統(tǒng)的溫度場(chǎng)計(jì)算與分析。三維熱路模型運(yùn)用靈活,針對(duì)復(fù)雜的電動(dòng)機(jī)系統(tǒng)將具有很好的適應(yīng)性,如高壓大電動(dòng)機(jī)、潛水電動(dòng)機(jī)、空氣壓縮機(jī)電動(dòng)機(jī)等,這些電動(dòng)機(jī)或功率大、絕緣要求高,或運(yùn)行環(huán)境復(fù)雜、安裝空間受限,對(duì)其溫度場(chǎng)的分析計(jì)算、散熱系統(tǒng)的設(shè)計(jì)提出了較為嚴(yán)苛的要求,但目前針對(duì)這些電動(dòng)機(jī)的熱場(chǎng)分析,三維熱路法的應(yīng)用研究還未見發(fā)表。

    (2) 基于三維熱路法的大型商用軟件的開發(fā)。目前,基于空間熱路模型的電動(dòng)機(jī)溫度場(chǎng)計(jì)算軟件較為成型的只有Motor-CAD,且主要集中于普通中、小型電動(dòng)機(jī)的溫度場(chǎng)分析,雖然也可以根據(jù)熱路的數(shù)學(xué)模型利用其他編程軟件進(jìn)行編程計(jì)算,但較為耗時(shí)和不便,故開發(fā)出精度更高、使用更為靈活的大型空間熱路分析軟件,是三維熱路法發(fā)展的必然趨勢(shì)。

    (3) 在電動(dòng)機(jī)實(shí)時(shí)溫度檢測(cè)系統(tǒng)中的應(yīng)用。電動(dòng)機(jī)是一個(gè)多物理場(chǎng)相互耦合的系統(tǒng),集電磁、溫度、流體和結(jié)構(gòu)于一身,因此,電動(dòng)機(jī)的溫度檢測(cè)一直是一個(gè)有待突破的技術(shù)難點(diǎn)?,F(xiàn)有的檢測(cè)方法無(wú)法準(zhǔn)確得到某一特定區(qū)域的溫度值,且溫度是實(shí)時(shí)變化的,就更加大了電動(dòng)機(jī)溫度檢測(cè)的難度,三維熱路法具有計(jì)算速度快的優(yōu)點(diǎn),因此,集成了三維熱路模型的電動(dòng)機(jī)溫度檢測(cè)系統(tǒng)也將會(huì)是其發(fā)展的一大趨勢(shì)。

    4結(jié)語(yǔ)

    針對(duì)電動(dòng)機(jī)溫度場(chǎng)的分析和計(jì)算問題,本文運(yùn)用三維空間等效熱路模型,詳細(xì)闡述了熱路法的基本原理,總結(jié)了立方體T型熱路、空心圓筒雙T型熱路、三維立方體和扇體熱路的數(shù)學(xué)模型。最后,結(jié)合熱路模型在電機(jī)溫度場(chǎng)計(jì)算中的運(yùn)用情況,探討了熱路模型未來(lái)的發(fā)展趨勢(shì),對(duì)電動(dòng)機(jī)的溫度場(chǎng)計(jì)算及其散熱系統(tǒng)的設(shè)計(jì)具有一定的借鑒意義和參考價(jià)值。

    參考文獻(xiàn):

    [1]覃新川,文宏剛,陸凌嘯.電機(jī)散熱計(jì)算方法探討[C]∥2011中國(guó)電工技術(shù)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集.北京: 中國(guó)電工技術(shù)學(xué)會(huì),2011: 678-681.

    [2]姚望,韓雪巖,唐任遠(yuǎn).車用永磁同步電機(jī)熱計(jì)算的分析與研究[C]∥2012中國(guó)·西安微特電機(jī)及永磁電機(jī)技術(shù)創(chuàng)新與發(fā)展論壇.西安: 中國(guó)電工技術(shù)學(xué)會(huì),中國(guó)電器工業(yè)協(xié)會(huì),西安自動(dòng)化學(xué)會(huì),2012.

    [3]陳志剛.等效熱網(wǎng)絡(luò)法和有限元法在電機(jī)三維溫度場(chǎng)計(jì)算中的應(yīng)用與比較[J].中小型電機(jī),1995,22(1): 3-6,35.

    [4]王紅宇,蘇鵬聲,王祥珩.采用熱路方法和三維有限元方法計(jì)算三峽水輪發(fā)電機(jī)三維熱網(wǎng)絡(luò)熱阻參數(shù)的研究[J].大電機(jī)技術(shù),2007(1): 11-18.

    [5]鮑曉華,盛海軍,單麗,等.高壓濕式潛水電機(jī)轉(zhuǎn)子三維溫度場(chǎng)分析[J].大電機(jī)技術(shù),2013(1): 11-14.

    [6]胡田,唐任遠(yuǎn),李巖.空-空冷中型異步電機(jī)流體場(chǎng)與溫度場(chǎng)計(jì)算及分析[J].電氣工程學(xué)報(bào),2015,10(7): 11-17.

    [7]Boglietti A,Cavagnino A,Lazzari M.A simplified thermal model for variable-speed self-cooled industrial induction motor.IEEE Transactions on Industry Application,2003,39(4): 945-952.

    [8]Chin Y K, Staton D A.Transient thermal analysis using both lumped-circuit approach and finite element method of a permanent magnet traction motor[C]∥7th AFRICON Conference in Africa.[S.L.]: IEEE,2004: 1027-1035

    [9]張琪,魯茜睿,黃蘇融,等.多領(lǐng)域協(xié)同仿真的高密度永磁電機(jī)溫升計(jì)算[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,(1): 1874-1881.

    [10]Nerg J,Rilla M,Pyrhonen J.Thermal analysis of radial-flux electrical machines with a high power density[J].IEEE Transactions on Industrial Electronics,2008,55(10): 3543-3554.

    [11]Gerling D,Dajaku G.Novel lumped-parameter thermal model for electrical systems[C]∥2005 European Conference on Power Electronics and Applications.[S.L.]: IEEE,2005: 10.

    [12]Wrobel R,Mellor P H.A general cuboidal element for three-dimensional thermal modelling.IEEE Transactions on Magnetics,2010,46(8): 3197-3200.

    [13]Simpson N.Wrobel R,Mellor P H.A general arc-segment element for three-dimensional thermal modeling.IEEE Transactions on Magnetics,2014,50(2): 265-268.

    [14]Fang Q.Stippich A,Guettler M,et al.Methodical considerations for setting up space-resolved lumped-parameter thermal models for electrical machines∥2014 17th International Conference on Electrical Machines and Systems (ICEMS).Hangzhou,China: IEEE,2014: 651-657.

    [15]Perez-Arriaga I J,Kassakian J K.A stationary thermal model for smooth air-gap rotating electrical machines.Electric Machines and Power Systems,1979,3(4): 285-303.

    [16]Mellor P,Roberts D,Turner D.Lumped parameter thermal model for electrical machines of TEFC design.IEE Proceedings-B,1991,138(5): 205-219.

    [17]Kylander G.Thermal modelling of small cage induction motors[M].[S.L.]: Chalmers University of Technology,1995.

    [18]Rouhani H,Faiz J,Lucas C.Lumped thermal model for switched reluctance motor applied to mechanical design optimization[J].Mathematical and Computer Modelling,2007,45(5): 625-638.

    [19]Gerling D,Dajaku G.Thermal calculation of systems with distributed heat generation[C]∥The Tenth Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronics Systems.[S.L.]: IEEE,2006: 8,652.

    [20]Li H,Shen Y.Thermal analysis of the permanent-magnet spherical motor[J].IEEE Transactions on Energy Conversion,2015: 1-8.

    [21]Staton D,Hawkins D,Popescu M.Thermal behaviour of electrical motors-an analytical approach[EB/OL].(2015-03-12).http: ∥motor-design.com/cmsAdmin/uploads/cwieme_2009_paper_mdl.pdf.

    [22]Staton D A.Servo motor size reduction-need for thermal CAD[C]∥Drives & Controls 2001 Conference 2001.Docklands,London,UK: Motor Design Ltd, 2001: 1-10.

    [23]Deaconu G,Ghit? C,Chiril? A I,et al.Thermal study of induction machine using Motor-CAD[C]∥2010 3rd International Symposium on Electrical and Electronics Engineering(ISEEE).[S.L.]: IEEE,2010: 23-27.

    3D Space Thermal Circuit Method and Application inCalculation of Motor Temperature Field

    TIANJingcheng,WANGFeiyu,ZHUOKeqiong,ZHAOChaohui

    (School of Electrical Engineering, Shanghai Dianji University, Shanghai 201306, China)

    Abstract:This paper aims at the use of a 3D space equivalent thermal circuit model in temperature field calculation, and discusses the basic principle of thermal circuits. The mathematical models of T-type cuboidal thermal circuits, hollow cylinder double T-type thermal circuits, and 3D cuboidal and arc segment thermal circuits are introduced. In conjunction with the applications of thermal circuit models in the calculation of motor temperature fields, the development tendency is discussed. This study provides references to the calculation of motor temperature and design of cooling systems.

    Key words:motor; temperature field; 3D space equivalent thermal circuit model; mathematical model

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

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

    文章編號(hào)2095 - 0020(2015)06 -0337 - 08

    作者簡(jiǎn)介:田井呈(1990 -),男,碩士生,主要研究方向?yàn)殡姍C(jī)設(shè)計(jì)及其多物理場(chǎng)仿真,E-mail: 526837119@qq.com

    基金項(xiàng)目:上海市教育委員會(huì)科研創(chuàng)新項(xiàng)目資助(13ZZ142)

    收稿日期:2015 - 09 - 25

    猜你喜歡
    數(shù)學(xué)模型溫度場(chǎng)電動(dòng)機(jī)
    AHP法短跑數(shù)學(xué)模型分析
    活用數(shù)學(xué)模型,理解排列組合
    淺析電動(dòng)機(jī)日常維護(hù)與保養(yǎng)
    永磁同步電動(dòng)機(jī)的節(jié)能計(jì)算
    鋁合金加筋板焊接溫度場(chǎng)和殘余應(yīng)力數(shù)值模擬
    基于紋影法的溫度場(chǎng)分布測(cè)量方法
    MJS工法與凍結(jié)法結(jié)合加固區(qū)溫度場(chǎng)研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    對(duì)一個(gè)數(shù)學(xué)模型的思考
    基于KB0的電動(dòng)機(jī)軟啟動(dòng)控制系統(tǒng)
    電動(dòng)機(jī)的保護(hù)配合及CPS在民用建筑中的應(yīng)用
    大香蕉久久网| 欧美性猛交黑人性爽| 日本爱情动作片www.在线观看| 久久久色成人| 九九爱精品视频在线观看| 欧美+日韩+精品| 91久久精品电影网| www.色视频.com| 岛国在线免费视频观看| 成人二区视频| 亚洲国产精品合色在线| 久久人人精品亚洲av| 变态另类丝袜制服| 成人性生交大片免费视频hd| 99久久精品一区二区三区| 亚洲经典国产精华液单| 草草在线视频免费看| 搡老妇女老女人老熟妇| 色视频www国产| 成人欧美大片| 国产三级中文精品| 国产精品女同一区二区软件| 国产日本99.免费观看| 国产伦精品一区二区三区四那| 国产极品精品免费视频能看的| 国产av在哪里看| 国产成人福利小说| 熟女人妻精品中文字幕| 国产成人午夜福利电影在线观看| 国内少妇人妻偷人精品xxx网站| 国产一区亚洲一区在线观看| av女优亚洲男人天堂| 国产色婷婷99| ponron亚洲| 国产亚洲91精品色在线| 99热全是精品| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 亚洲av第一区精品v没综合| 啦啦啦啦在线视频资源| 天天一区二区日本电影三级| 男女做爰动态图高潮gif福利片| av在线播放精品| 国产精品久久久久久久久免| 色综合色国产| 国产伦精品一区二区三区视频9| 久久韩国三级中文字幕| 老师上课跳d突然被开到最大视频| 岛国毛片在线播放| 日韩精品有码人妻一区| 日韩高清综合在线| 啦啦啦韩国在线观看视频| 中文亚洲av片在线观看爽| 亚洲综合色惰| 一本久久中文字幕| 久久这里有精品视频免费| 99久久久亚洲精品蜜臀av| 五月玫瑰六月丁香| 久久精品夜色国产| 国产免费男女视频| 欧美日韩精品成人综合77777| 国产一区二区三区在线臀色熟女| 欧美丝袜亚洲另类| 久久久久久久久久成人| 久久99热这里只有精品18| 欧美在线一区亚洲| 国产在视频线在精品| 女同久久另类99精品国产91| 变态另类丝袜制服| 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 午夜a级毛片| 男女视频在线观看网站免费| 又黄又爽又刺激的免费视频.| 你懂的网址亚洲精品在线观看 | 一级毛片电影观看 | 中文字幕免费在线视频6| 亚洲国产精品合色在线| 97超视频在线观看视频| 一本一本综合久久| 欧美变态另类bdsm刘玥| 蜜臀久久99精品久久宅男| 国产精品精品国产色婷婷| 久久99热6这里只有精品| 欧美日本视频| 亚洲成人精品中文字幕电影| 欧美xxxx黑人xx丫x性爽| 国产 一区 欧美 日韩| 亚洲成人av在线免费| www.色视频.com| 性色avwww在线观看| 亚洲在线自拍视频| 日产精品乱码卡一卡2卡三| 欧美一区二区亚洲| 蜜桃亚洲精品一区二区三区| 亚洲国产高清在线一区二区三| 国产一区亚洲一区在线观看| 久久99热6这里只有精品| 免费看光身美女| 少妇丰满av| 少妇的逼好多水| 一本久久中文字幕| 国产成人精品一,二区 | 99在线视频只有这里精品首页| 美女被艹到高潮喷水动态| 亚洲激情五月婷婷啪啪| 欧美一区二区亚洲| av黄色大香蕉| 日韩人妻高清精品专区| 精品免费久久久久久久清纯| 亚洲天堂国产精品一区在线| 亚洲成人久久性| 国产精品爽爽va在线观看网站| 99热这里只有是精品在线观看| 不卡视频在线观看欧美| 精品一区二区免费观看| 能在线免费观看的黄片| 亚洲最大成人中文| 久久久久久久久久黄片| 麻豆精品久久久久久蜜桃| 国产精品乱码一区二三区的特点| 亚洲自偷自拍三级| 亚洲一区高清亚洲精品| h日本视频在线播放| 我的女老师完整版在线观看| 一级黄片播放器| 一个人看视频在线观看www免费| 黄片wwwwww| 日韩高清综合在线| 日韩成人伦理影院| www.色视频.com| 欧美色视频一区免费| 国产色婷婷99| 久久人妻av系列| 色哟哟·www| 最近中文字幕高清免费大全6| 91在线精品国自产拍蜜月| 国产黄色视频一区二区在线观看 | 精品人妻熟女av久视频| 国产在视频线在精品| 久久久久免费精品人妻一区二区| av.在线天堂| 亚洲成人精品中文字幕电影| 黑人高潮一二区| 天美传媒精品一区二区| 好男人在线观看高清免费视频| 成熟少妇高潮喷水视频| 亚洲无线观看免费| 精品国产三级普通话版| 国产精品一区www在线观看| 亚洲无线观看免费| 亚洲精品色激情综合| 非洲黑人性xxxx精品又粗又长| 免费看光身美女| 亚洲av.av天堂| 最近2019中文字幕mv第一页| 亚洲色图av天堂| 日韩欧美三级三区| 如何舔出高潮| 亚洲第一区二区三区不卡| 欧美性猛交黑人性爽| 人妻久久中文字幕网| av.在线天堂| 91久久精品国产一区二区三区| 青春草国产在线视频 | 青青草视频在线视频观看| 国内久久婷婷六月综合欲色啪| 一进一出抽搐gif免费好疼| 一级二级三级毛片免费看| 高清毛片免费看| 人妻夜夜爽99麻豆av| 一级毛片久久久久久久久女| 级片在线观看| 久久人人精品亚洲av| 欧美xxxx黑人xx丫x性爽| 国产欧美日韩精品一区二区| 久久这里只有精品中国| 麻豆av噜噜一区二区三区| 变态另类成人亚洲欧美熟女| 国产一区二区亚洲精品在线观看| 亚洲不卡免费看| 麻豆久久精品国产亚洲av| 欧美性感艳星| 成人综合一区亚洲| 一区福利在线观看| 久久99热6这里只有精品| 一级黄色大片毛片| 99九九线精品视频在线观看视频| 亚洲在线自拍视频| 久久亚洲精品不卡| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 18禁在线播放成人免费| 又爽又黄a免费视频| 22中文网久久字幕| 在线观看午夜福利视频| 国产午夜精品论理片| 亚洲av.av天堂| 韩国av在线不卡| 变态另类成人亚洲欧美熟女| 久久久精品欧美日韩精品| 欧美成人免费av一区二区三区| 成人特级黄色片久久久久久久| 中文字幕熟女人妻在线| 国产色婷婷99| 免费av毛片视频| 好男人在线观看高清免费视频| 欧美激情国产日韩精品一区| 亚洲内射少妇av| 午夜精品国产一区二区电影 | 性色avwww在线观看| 久久久久久久久中文| 亚洲欧洲国产日韩| 麻豆av噜噜一区二区三区| 精品熟女少妇av免费看| 性插视频无遮挡在线免费观看| 久久这里有精品视频免费| 国产高清激情床上av| 国产美女午夜福利| 免费观看在线日韩| 丝袜美腿在线中文| 欧美潮喷喷水| 男女啪啪激烈高潮av片| 两个人的视频大全免费| 麻豆一二三区av精品| 两个人的视频大全免费| 国产真实伦视频高清在线观看| 日韩欧美三级三区| 欧美xxxx黑人xx丫x性爽| 如何舔出高潮| 国产午夜精品一二区理论片| 久久热精品热| 国国产精品蜜臀av免费| 插阴视频在线观看视频| 国产久久久一区二区三区| 一个人看视频在线观看www免费| 自拍偷自拍亚洲精品老妇| 国产精品女同一区二区软件| 日本黄色片子视频| 成人永久免费在线观看视频| 91午夜精品亚洲一区二区三区| 亚洲18禁久久av| 最后的刺客免费高清国语| 亚洲久久久久久中文字幕| 久久久久国产网址| 99九九线精品视频在线观看视频| 麻豆一二三区av精品| 亚洲av.av天堂| 国产精品伦人一区二区| 亚洲成人久久性| 嫩草影院入口| 尾随美女入室| 一级毛片电影观看 | 内地一区二区视频在线| 亚洲成人中文字幕在线播放| 久久99热6这里只有精品| 不卡一级毛片| 69av精品久久久久久| 午夜视频国产福利| 在线免费观看的www视频| 嘟嘟电影网在线观看| 国产乱人视频| 简卡轻食公司| 不卡一级毛片| 国产精华一区二区三区| 国内精品一区二区在线观看| 日本免费a在线| 99久久精品一区二区三区| 国产私拍福利视频在线观看| 如何舔出高潮| 少妇人妻一区二区三区视频| 亚洲欧美中文字幕日韩二区| 久久久久久久久中文| 97人妻精品一区二区三区麻豆| 国产男人的电影天堂91| 在线观看免费视频日本深夜| 欧美bdsm另类| 亚洲婷婷狠狠爱综合网| 欧美成人精品欧美一级黄| 国产伦在线观看视频一区| 国产免费一级a男人的天堂| 国产成人91sexporn| 日韩高清综合在线| 亚洲av成人精品一区久久| 国产精品一区二区性色av| a级毛片a级免费在线| 麻豆乱淫一区二区| 日韩欧美在线乱码| 少妇人妻一区二区三区视频| 欧美日韩乱码在线| 久久精品国产亚洲av天美| 日本黄大片高清| 国产精品日韩av在线免费观看| 青春草国产在线视频 | 深夜精品福利| 亚洲国产精品久久男人天堂| 日韩中字成人| 免费看日本二区| 欧美精品一区二区大全| 欧美性猛交╳xxx乱大交人| 美女 人体艺术 gogo| 日本一二三区视频观看| 国产激情偷乱视频一区二区| 在线观看一区二区三区| 国产精品一区www在线观看| 国语自产精品视频在线第100页| 亚洲精品国产成人久久av| 久久草成人影院| 一个人看的www免费观看视频| 国产成人福利小说| 亚洲美女搞黄在线观看| 边亲边吃奶的免费视频| 青青草视频在线视频观看| 别揉我奶头 嗯啊视频| 九九在线视频观看精品| 99热只有精品国产| 美女cb高潮喷水在线观看| 免费一级毛片在线播放高清视频| 日韩人妻高清精品专区| 超碰av人人做人人爽久久| 内地一区二区视频在线| 国产午夜福利久久久久久| 一边亲一边摸免费视频| 国产中年淑女户外野战色| 欧美人与善性xxx| 久久精品国产亚洲av香蕉五月| 岛国在线免费视频观看| 黄色欧美视频在线观看| 国产精品福利在线免费观看| 两个人视频免费观看高清| 一本精品99久久精品77| 男人舔女人下体高潮全视频| 校园春色视频在线观看| 亚洲欧美成人综合另类久久久 | 啦啦啦观看免费观看视频高清| 我的女老师完整版在线观看| 国产高清三级在线| 欧美精品国产亚洲| 可以在线观看毛片的网站| 亚州av有码| 高清在线视频一区二区三区 | 男女边吃奶边做爰视频| 看片在线看免费视频| 12—13女人毛片做爰片一| 六月丁香七月| 国产老妇女一区| 一区二区三区免费毛片| 在线天堂最新版资源| videossex国产| 白带黄色成豆腐渣| 性插视频无遮挡在线免费观看| 欧美激情久久久久久爽电影| 免费搜索国产男女视频| 欧美最新免费一区二区三区| 国产成人影院久久av| 小说图片视频综合网站| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 成人性生交大片免费视频hd| 亚洲精品日韩在线中文字幕 | 色噜噜av男人的天堂激情| 99久国产av精品| 深夜a级毛片| 国产一级毛片在线| 男女下面进入的视频免费午夜| 老司机福利观看| 久久精品国产自在天天线| 老司机福利观看| www日本黄色视频网| 午夜福利成人在线免费观看| 日韩制服骚丝袜av| 男女那种视频在线观看| 亚洲内射少妇av| 久久这里有精品视频免费| 成年女人永久免费观看视频| 国产激情偷乱视频一区二区| 久久久精品94久久精品| 尾随美女入室| 亚洲av熟女| 国产白丝娇喘喷水9色精品| av卡一久久| 三级男女做爰猛烈吃奶摸视频| 成人无遮挡网站| 美女脱内裤让男人舔精品视频 | 网址你懂的国产日韩在线| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久久毛片| 禁无遮挡网站| 日韩人妻高清精品专区| 国产69精品久久久久777片| 国产一区二区三区av在线 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产 一区精品| 日本欧美国产在线视频| 国产精品电影一区二区三区| 国产乱人偷精品视频| 国产极品精品免费视频能看的| 综合色丁香网| 美女被艹到高潮喷水动态| 美女黄网站色视频| 在线免费观看不下载黄p国产| 久久精品国产亚洲av涩爱 | 一本久久中文字幕| 国产极品精品免费视频能看的| 久久久久网色| 大型黄色视频在线免费观看| 麻豆av噜噜一区二区三区| 哪个播放器可以免费观看大片| 精品日产1卡2卡| 婷婷六月久久综合丁香| 亚洲在线自拍视频| 淫秽高清视频在线观看| 日本熟妇午夜| 97热精品久久久久久| 你懂的网址亚洲精品在线观看 | 国产欧美日韩精品一区二区| 国产精品精品国产色婷婷| 搞女人的毛片| 午夜福利在线观看吧| 久久精品91蜜桃| 国产精品一区www在线观看| 日日啪夜夜撸| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香| 成人综合一区亚洲| 狂野欧美激情性xxxx在线观看| 久久久久久久久久成人| 精品久久久久久久末码| 身体一侧抽搐| 人妻系列 视频| 免费观看a级毛片全部| 亚洲图色成人| 亚洲国产欧洲综合997久久,| 在线播放无遮挡| 狠狠狠狠99中文字幕| 国产精品1区2区在线观看.| 国产亚洲精品久久久com| ponron亚洲| 日本黄色片子视频| 国产蜜桃级精品一区二区三区| 成人美女网站在线观看视频| 国产精品爽爽va在线观看网站| 内射极品少妇av片p| 亚洲高清免费不卡视频| 日韩欧美 国产精品| 国产视频内射| 搡女人真爽免费视频火全软件| 国产高清三级在线| 噜噜噜噜噜久久久久久91| 国产精品乱码一区二三区的特点| 能在线免费看毛片的网站| 国产成人精品久久久久久| 最近2019中文字幕mv第一页| www日本黄色视频网| 天堂影院成人在线观看| 有码 亚洲区| av在线亚洲专区| 免费搜索国产男女视频| 国产免费一级a男人的天堂| 人妻系列 视频| avwww免费| 日本一本二区三区精品| 亚洲精品国产成人久久av| 国产精品乱码一区二三区的特点| 亚洲精品国产av成人精品| 天天躁夜夜躁狠狠久久av| 欧美不卡视频在线免费观看| а√天堂www在线а√下载| 在线a可以看的网站| 午夜福利视频1000在线观看| 2022亚洲国产成人精品| 国产久久久一区二区三区| 日韩欧美三级三区| 国产亚洲精品av在线| 久久韩国三级中文字幕| 一区二区三区高清视频在线| 亚洲一区高清亚洲精品| 我的女老师完整版在线观看| 日韩欧美在线乱码| 六月丁香七月| 中文字幕av成人在线电影| 国产一级毛片在线| 亚洲精品色激情综合| 欧美高清成人免费视频www| а√天堂www在线а√下载| 精品不卡国产一区二区三区| 美女 人体艺术 gogo| 国产午夜精品久久久久久一区二区三区| 人妻少妇偷人精品九色| 国产三级中文精品| 国产91av在线免费观看| 久久久久久久久久黄片| 日韩视频在线欧美| 国产亚洲91精品色在线| 欧美+日韩+精品| 国产成人影院久久av| 男人舔奶头视频| 亚洲人成网站在线播| 欧美精品国产亚洲| 伊人久久精品亚洲午夜| 久久九九热精品免费| 久久99蜜桃精品久久| 级片在线观看| 最近视频中文字幕2019在线8| 国产精品1区2区在线观看.| 久99久视频精品免费| eeuss影院久久| 国产色婷婷99| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2| 桃色一区二区三区在线观看| 亚洲av.av天堂| 欧美高清成人免费视频www| 六月丁香七月| 亚洲精品自拍成人| 国产精品永久免费网站| 日韩强制内射视频| 小说图片视频综合网站| 欧美色欧美亚洲另类二区| 免费av毛片视频| 日韩在线高清观看一区二区三区| 能在线免费看毛片的网站| 2022亚洲国产成人精品| 午夜久久久久精精品| 青春草亚洲视频在线观看| 小说图片视频综合网站| 国产视频内射| 亚洲av免费在线观看| 伦精品一区二区三区| 国产精华一区二区三区| 国产真实伦视频高清在线观看| 久久久久久久亚洲中文字幕| 精品免费久久久久久久清纯| 亚洲成人久久性| 欧美成人免费av一区二区三区| АⅤ资源中文在线天堂| 欧美不卡视频在线免费观看| 亚洲人与动物交配视频| www.色视频.com| 国产精品不卡视频一区二区| 国产成人影院久久av| 2021天堂中文幕一二区在线观| 久久人妻av系列| 午夜亚洲福利在线播放| 色哟哟哟哟哟哟| 中文字幕制服av| 1024手机看黄色片| 日本色播在线视频| 中文字幕久久专区| 两个人视频免费观看高清| 国产熟女欧美一区二区| 桃色一区二区三区在线观看| 深夜精品福利| www.av在线官网国产| 国产成人aa在线观看| 久久这里只有精品中国| 在线播放无遮挡| 美女cb高潮喷水在线观看| 亚洲国产精品sss在线观看| 国产色婷婷99| 免费搜索国产男女视频| 欧美精品一区二区大全| 国产久久久一区二区三区| 亚洲最大成人av| 中文亚洲av片在线观看爽| 成人二区视频| 床上黄色一级片| 国产精品1区2区在线观看.| 免费在线观看成人毛片| 国产精品,欧美在线| 国产探花极品一区二区| 人妻少妇偷人精品九色| 免费av毛片视频| 欧美+日韩+精品| 麻豆成人午夜福利视频| 国产伦精品一区二区三区视频9| 在线播放无遮挡| 国内精品宾馆在线| 男女下面进入的视频免费午夜| 人人妻人人看人人澡| 男人的好看免费观看在线视频| 18禁黄网站禁片免费观看直播| 国产精品电影一区二区三区| 久久午夜福利片| 91aial.com中文字幕在线观看| 精品国内亚洲2022精品成人| 亚洲精品影视一区二区三区av| 国产三级在线视频| 青春草视频在线免费观看| 性插视频无遮挡在线免费观看| 免费看a级黄色片| 成人综合一区亚洲| 亚洲五月天丁香| 亚洲成人精品中文字幕电影| 欧美潮喷喷水| 一级毛片我不卡| av国产免费在线观看| 亚洲人成网站在线观看播放| 午夜福利视频1000在线观看| 26uuu在线亚洲综合色| 亚洲成av人片在线播放无| 亚洲四区av| 亚洲欧洲国产日韩| 亚洲性久久影院| 午夜爱爱视频在线播放| 国产私拍福利视频在线观看| 人妻制服诱惑在线中文字幕| 在线免费观看不下载黄p国产| 日韩三级伦理在线观看| 国产精品一区二区三区四区免费观看| 免费在线观看成人毛片| 国产亚洲91精品色在线| 三级男女做爰猛烈吃奶摸视频| 嘟嘟电影网在线观看| 亚洲精品国产成人久久av|