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

    地鐵車輛牽引逆變器熱管散熱器的溫升試驗及熱仿真

    2016-04-10 01:45:07杰,張
    中國鐵道科學 2016年3期
    關(guān)鍵詞:風罩瞬態(tài)熱管

    丁 杰,張 平

    (1.湘潭大學 土木工程與力學學院,湖南 湘潭 411105;2.南車株洲電力機車研究所有限公司 南車電氣技術(shù)與材料工程研究院,湖南 株洲 412001)

    牽引逆變器是地鐵車輛的關(guān)鍵部件,而絕緣柵雙極型晶體管(Insulated Gate Bipolar Transistor,IGBT)模塊是牽引逆變器中最為核心的器件。根據(jù)IGBT模塊的失效機理可知,其在封裝時各層材料的熱膨脹系數(shù)不一致,在長期高溫熱循環(huán)作用下可發(fā)生鋁鍵合線斷裂或脫落、硅芯片與襯板之間及襯板與基板之間的焊料層老化、柵極氧化層損壞和芯片失效等[1-3]。因此,IGBT模塊的散熱設(shè)計至關(guān)重要。

    設(shè)計地鐵車輛牽引逆變器時,大多數(shù)采用熱管散熱器加走行風冷的方式,這樣可以省去風道和風機,簡化牽引逆變器的柜體結(jié)構(gòu),減少風機故障隱患和振動噪聲的影響[4-5]。但由于地鐵線路站點多、站間距短,列車啟停、加減速頻繁,使牽引逆變器的電氣參數(shù)不斷變化,導致IGBT模塊的損耗也實時改變;加之走行風速有快有慢,使熱管散熱器的散熱能力也相應變化,因此僅計算熱管散熱器和IGBT模塊某一工況的穩(wěn)態(tài)溫升是無法準確評估熱管散熱器的性能和IGBT模塊的疲勞壽命的,有必要對其瞬態(tài)溫升進行計算。

    目前,研究熱管散熱器和IGBT模塊溫升的方法可分為試驗分析法[6]和仿真分析法[7]。試驗分析法可在實驗室或?qū)嶋H運行的車輛中進行,優(yōu)點是測試的溫升數(shù)據(jù)可信度高,缺點在于無法直接測量封裝于絕緣材料中芯片及焊料等部位的溫度。仿真分析法大多基于電氣、熱等學科建立仿真模型并進行單獨計算。電氣仿真借助Matlab/Simulink,Saber和Simplorer等軟件建立主電路仿真模型,得到IGBT模塊的電壓、電流等信息,進而計算芯片損耗。熱仿真主要有熱阻抗網(wǎng)絡(luò)、有限元、有限體積和模型降階等方法。熱阻抗網(wǎng)絡(luò)法將三維結(jié)構(gòu)的導熱物體映射成為熱阻和熱容串級連接的一維熱阻抗網(wǎng)絡(luò)模型,使用方便、計算量很小,然而結(jié)溫、殼溫和熱管散熱器溫度是虛擬的理論平均值,與實際三維結(jié)構(gòu)不同部位表現(xiàn)出明顯的溫度梯度有較大差異,且難以考慮同一熱管散熱器上不同IGBT模塊之間以及同一IGBT模塊中不同芯片之間的熱耦合效應[7]。有限元法需要將IGBT模塊和熱管散熱器的幾何結(jié)構(gòu)劃分為有限元網(wǎng)格,在熱管散熱器與冷卻介質(zhì)的接觸面設(shè)置等效對流換熱系數(shù),而且計算時用較小的時間步長才能得到較為準確的瞬態(tài)溫度場計算結(jié)果,因此對計算機資源要求較高且求解速度慢。有限體積法是計算流體動力學(Computational Fluid Dynamic,CFD)中應用最為廣泛的方法,需要將IGBT模塊、熱管散熱器和冷卻介質(zhì)劃分為有限體積網(wǎng)格,自動計算流固耦合面的對流換熱系數(shù),其計算結(jié)果比有限元法更符合實際情況,然而邊界層對流動與換熱的影響使得網(wǎng)格尺寸被限制為較小的值,導致網(wǎng)格數(shù)量大幅度增加至數(shù)千萬,與此同時,流體區(qū)域的壓力與速度耦合,需要消耗大量的計算機資源進行迭代,才能得到收斂的穩(wěn)態(tài)計算結(jié)果;并且,計算瞬態(tài)問題時,每一時間子步的迭代計算需要消耗計算機資源且計算結(jié)果文件需要數(shù)GB的硬盤存儲空間,因此,瞬態(tài)CFD計算的代價遠遠超出有限元法,極少應用于實際工程中的復雜對象。模型降階法[8-11]采用Krylov子空間、平衡截斷、正交分解、積分全等變換等算法得到階次較低的近似模型,可以在較高準確性基礎(chǔ)上極大提高計算效率。Simplorer和Caspoc等部分商業(yè)軟件包含了模型降階的功能,可以實現(xiàn)多個軟件的聯(lián)合仿真,由于商業(yè)軟件為了保證計算效率,輸出結(jié)果是多個探測點的變量變化曲線,很難輸出和云圖顯示整個物理場的數(shù)據(jù),因此不能完全滿足實際的應用需求。

    為快速、準確地分析地鐵車輛牽引逆變器熱管散熱器的性能及其IGBT模塊的溫升,本文首先搭建單個變流模塊的溫升試驗裝置,通過試驗得到可用于熱管散熱器熱仿真計算的輸入條件及驗證仿真結(jié)果準確性的試驗數(shù)據(jù);然后采用FLUENT軟件對熱管散熱器進行穩(wěn)態(tài)CFD仿真,得到其溫度場和對流換熱系數(shù)的分布;再將有限元法與模型降階法相結(jié)合,開發(fā)快速計算程序,最后通過Matlab/Simulink軟件實現(xiàn)熱與電氣的聯(lián)合仿真,以有效解決復雜工況的瞬態(tài)計算問題。

    1 溫升試驗

    1.1 試驗裝置

    某地鐵車輛牽引逆變器采用二電平電壓型直—交逆變主電路,其主要參數(shù)為:額定電壓DC750 V,額定輸出容量600 kVA,輸出電壓0~585 V,輸出頻率0~160 Hz,開關(guān)頻率500 Hz,額定輸出電流600 A,牽引工況最大輸出電流有效值884 A,制動工況最大輸出電流有效值1 156 A,額定工作點效率0.98,采用直接轉(zhuǎn)矩控制方式。當車輛處于牽引工況時,直流電經(jīng)過高壓電器箱、濾波電抗器箱分別對2個IBCM60G1型變流模塊供電,經(jīng)牽引逆變器輸出三相變壓變頻的交流電向4臺異步牽引電機供電。當車輛處于再生制動工況時,牽引逆變器將異步牽引電機輸出的三相交流電變成直流電并反饋回電網(wǎng),或由制動電阻消耗掉。三相交流電逆變及電阻制動消耗環(huán)節(jié)的開關(guān)管均為FZ1600R17KF6C-B2型IGBT模塊。

    考慮到IGBT模塊的芯片被絕緣材料封裝起來,無法對其直接進行溫度測量,因此在熱管散熱器的安裝面上布置PT100熱電阻進行溫度測量。熱管散熱器上測溫點的布置情況如圖1所示。圖中:V2—V7和V1,V8分別為實現(xiàn)三相逆變和電阻制動功能的IGBT模塊;J1—J3為溫度繼電器;1#—6#為相鄰的IGBT模塊間的測溫點;7#為J2的測溫點。

    圖1 測溫點布置情況

    在實驗室條件下進行溫升試驗的基本步驟為:首先在變流模塊組裝之前于熱管散熱器上布置測溫點,然后將其與總裝后的牽引逆變器、通風機、風罩、電源和負載(牽引電機)等相連,再將測溫點的數(shù)據(jù)線連接至溫度巡檢儀和計算機,最后進行試驗測試和數(shù)據(jù)分析。

    考慮到實驗室條件下不能真實模擬走行風且只有1臺牽引電機作為負載,因此重點研究在減小變流模塊1的輸入電流時熱管散熱器的溫升情況。在變流模塊1的翅片部分設(shè)計1個風罩,通過軸流通風機進行吹風,連接通風機的風罩部分為圓形,包裹翅片的風罩部分為矩形,風罩中間部分采用漸變的形式并設(shè)置整流隔柵,以使通風機的出風盡量均勻地吹向熱管散熱器。在風罩入口處布置1個熱電偶,用于測量牽引逆變器附近的空氣溫度。搭建的試驗裝置如圖2所示。

    圖2 試驗裝置

    1.2 試驗結(jié)果

    在風罩入口和出口處分別測得9個點的風速值,如圖3所示。由圖3可以看出:與入口處其他部位相比,中心點處的風速因通風機的阻擋而偏小,入口上部的風速高于下部,而冷卻空氣經(jīng)過通風機、風罩、熱管散熱器翅片及外罩后,出口上部的風速低于下部,出口處風速與入口處風速不對應的情況說明風罩內(nèi)的空氣流動是非常復雜的。

    圖3 風罩入口和出口處的風速測量結(jié)果(單位: m·s-1)

    根據(jù)風罩入口處的風速測量結(jié)果,可以計算出平均風速約為5.7 m·s-1。

    圖4為由各測溫點得到的溫度變化曲線。由圖4可以看出:試驗開始后0~30 min內(nèi)溫度上升速度很快,30 min后上升速度逐漸變緩;風罩入口處空氣的溫度為31.1 ℃,持續(xù)試驗220 min后升為33.8 ℃,這是由于試驗場地的空間較為狹窄,變流模塊1產(chǎn)生的熱量難以有效耗散至實驗室外部,且牽引逆變器周圍還存在其他工作中的試驗設(shè)備,導致風罩入口處的空氣溫度不斷上升;隨著試驗時間的推移,各測溫點的溫度與風罩入口處空氣溫度的差值趨于平穩(wěn),說明熱管散熱器的溫升達到了平衡狀態(tài),因此風罩入口處空氣溫度的變化在試驗誤差允許范圍內(nèi);測溫點中4#的溫度最高、3#的溫度次之,而6#的溫度最低,這主要是由于V1和V8未工作且處于冷卻空氣上游、持續(xù)工作的V4和V5處于冷卻空氣最下游所致。

    圖4 溫升試驗得到的溫度變化曲線

    2 CFD穩(wěn)態(tài)熱仿真分析

    2.1 CFD仿真模型

    熱管散熱器內(nèi)部涉及復雜的相變過程,目前還很難對其本身進行數(shù)值模擬。較為常用的方法是將熱管假設(shè)為軸向?qū)嵯禂?shù)很高、徑向?qū)嵯禂?shù)為熱管材料導熱系數(shù)的實體棒桿[4-6]。由于IGBT模塊中包含AlSiC基板、焊料、銅層、AlN基片、芯片和絕緣材料,因此結(jié)合各層材料的厚度和整個模型的復雜程度,將仿真模型中網(wǎng)格的基本尺寸設(shè)置為2 mm,遠離熱管散熱器的流體區(qū)域網(wǎng)格尺寸設(shè)置為5 mm。以風罩為參照,將被風罩包裹的變流模塊1的流體外邊界設(shè)置為壁面邊界條件;變流模塊2因為沒有被外罩包裹,故將其流體區(qū)域適度擴大,并將其流體外邊界設(shè)置為壓力出口條件。利用HyperMesh軟件建立六面體為主的高質(zhì)量網(wǎng)格,數(shù)量為3 862萬個,其中流體區(qū)域的網(wǎng)格數(shù)量為2 520萬個。

    為便于熱仿真結(jié)果與溫升試驗數(shù)據(jù)的對比,以溫升試驗條件作為熱仿真計算的輸入條件??紤]到風罩內(nèi)的空氣流速分布極為不規(guī)則,為簡化仿真條件的參數(shù)設(shè)置,將風罩入口處的風速設(shè)為5.7 m·s-1、空氣溫度設(shè)為33.8 ℃。假設(shè)冷卻空氣在風道內(nèi)的流動狀態(tài)為完全湍流,采用標準k—ε模型進行模擬。利用FLUENT軟件和DELL T7600臺式工作站對CFD仿真模型進行穩(wěn)態(tài)計算,可以得到收斂的仿真結(jié)果。

    2.2 穩(wěn)態(tài)仿真結(jié)果

    圖5為冷卻空氣在試驗裝置內(nèi)流動時不同截面處冷卻空氣的流速分布云圖。由圖5可以看出:冷卻空氣受熱管散熱器翅片、熱管和外罩的阻擋,在不同部位處的流速分布并不相等,在局部截面面積小的區(qū)域其流速可達到最大值,為16.61 m·s-1。

    圖5 不同截面的冷卻空氣流速分布云圖

    圖6為經(jīng)過87 h仿真計算得到的變流模塊1的熱管散熱器及IGBT模塊溫度場分布云圖。由圖6可以看出:因冷卻空氣主要沿z坐標軸正方向流動,處于冷卻空氣最下游的V4和V5這2個IGBT模塊的溫度最高,可達55.43 ℃;V1和V8模塊的溫度最低,這是由于它們不發(fā)熱且位于冷卻空氣的上游;1#—7#測溫點的溫度分別為39.65,45.81,49.02,48.95,45.61,39.49和45.33 ℃,其中3#點的溫度略高于4#點的,這與溫升試驗結(jié)果有區(qū)別的主要原因是試驗中通風機入口處冷卻空氣的風速不均勻,而仿真時假設(shè)入口處的冷卻空氣風速完全相等;對比仿真結(jié)果與試驗數(shù)據(jù)可知,各測溫點溫度的相對誤差不超過5%,說明仿真結(jié)果的準確性。

    圖6 熱管散熱器及IGBT模塊的溫度場分布云圖

    圖7為變流模塊1的熱管散熱器對流換熱系數(shù)分布云圖。由圖7可以看出:冷卻空氣流經(jīng)散熱片、熱管、外罩和底板時,會產(chǎn)生不同程度的對流換熱作用,從而表現(xiàn)出對流換熱系數(shù)值的差異,即對流換熱系數(shù)的分布并不一致;對流換熱系數(shù)值大的區(qū)域表示對流換熱較為充分,反之則不充分。

    圖7 熱管散熱器的對流換熱系數(shù)分布云圖

    對流換熱系數(shù)分布處于CFD模型中流體和固體相接觸的流固耦合面上,流體與固體之間的熱量通過這些面進行傳遞,因此在有限元模型只需考慮固體和流固耦合面的換熱,則對流換熱系數(shù)分布是有限元模型中非常關(guān)鍵的邊界條件。在大多數(shù)有限元法和模型降階法的工程應用中,將邊界條件即對流換熱系數(shù)設(shè)為常數(shù),這種處理方式與實際情況有差異,這也是有限元法的計算精度低于有限體積法的根本原因。

    3 快速瞬態(tài)熱仿真分析

    3.1 理論基礎(chǔ)

    圖6所示的溫度場分布是經(jīng)過較長計算時間得到的穩(wěn)態(tài)結(jié)果,而圖4所示的溫度變化曲線是通過溫升試驗得到的瞬態(tài)數(shù)據(jù),故為了更好地實現(xiàn)仿真結(jié)果與溫升試驗數(shù)據(jù)的對比,需要進行瞬態(tài)熱仿真??紤]到CFD仿真的準確性高但計算效率非常低,故在滿足計算精度的前提下為進一步實現(xiàn)熱管散熱器和IGBT模塊瞬態(tài)導熱問題的快速計算,可將圖7所示的對流換熱系數(shù)分布插值至有限元模型中,以保證有限元模型與CFD仿真模型的邊界條件一致,再結(jié)合有限元法及模型降階法中最基本和最重要的Krylov子空間法[8-11]進行程序開發(fā)。

    利用有限元法計算熱管散熱器和IGBT模塊等固體區(qū)域的溫度場時,瞬態(tài)導熱問題的微分方程為

    (1)

    式中:ρ為密度;cp為定壓比熱容;T為溫度;t為時間;λx,λy和λz為沿物體3個主方向的導熱系數(shù);q為物體內(nèi)部的熱流密度。

    有限元模型的邊界條件分為給定溫度、給定熱流密度和給定對流換熱系數(shù)3種。由微分方程等效積分形式的伽遼金提法,在空間域離散后得到包含n個有限元模型節(jié)點數(shù)目的1階常微分方程組矩陣方程如下。

    (2)

    式中:Cp為熱容矩陣;K為熱傳導矩陣;P為溫度載荷列陣;φ為節(jié)點溫度列陣。

    式(2)涉及到大規(guī)模矩陣運算,消耗的計算機資源較大,Krylov子空間法可以將物理模型的大規(guī)模狀態(tài)空間投影到以1組基矢量為特征的低維空間(降階階數(shù)r通常取5~50,主要取決于時間步長[12],該數(shù)值遠小于原有階數(shù)n),從而有效降低矩陣的維度,并在保證計算結(jié)果準確性的基礎(chǔ)上大幅度降低計算規(guī)模來提高計算效率。

    3.2 快速計算程序

    文獻[13]介紹了模型降階法在IGBT風冷散熱器上進行瞬態(tài)問題快速計算的算例,并給出了提取熱傳導矩陣及熱容矩陣、求解降階模型的部分源程序。針對算例中只能計算單一熱源、不同軟件間手動操作、計算效率存在提升空間等問題,本文進行了改進并開發(fā)了快速計算程序。具體為:針對不同的IGBT模塊芯片依次設(shè)置不同熱源的參數(shù)并提取相應的矩陣文件,降階后的矩陣使用常微分方程求解器進行求解,然后將結(jié)果投影到原模型上。通過時間步與循環(huán)功能重復熱源設(shè)置、求解和投影等步驟,自動實現(xiàn)多熱源多載荷步的瞬態(tài)計算,其中的時間步由固定步長改為變步長,大幅度減少載荷子步的計算量。借助Matlab軟件的調(diào)用功能實現(xiàn)程序的自動化,大大簡化程序使用的復雜度。使用GRUS稀疏矩陣求解器的算法代替Matlab軟件自帶的矩陣分解算法,從而進一步提高了模型降階法的計算效率。

    3.3 仿真與試驗的對比分析

    利用HyperMesh軟件,網(wǎng)格基本尺寸取8 mm,將變流模塊1的熱管散熱器和IGBT模塊劃分為如圖8所示的有限元粗糙網(wǎng)格,其中單元數(shù)目為212 125個,節(jié)點數(shù)目為165 560個。

    將圖7所示的對流換熱系數(shù)分布插值到有限元模型的流固耦合面上,且以試驗工況為輸入條件,冷卻空氣的初始溫度取33.8 ℃,利用ANSYS軟件可以在1 min內(nèi)計算出如圖9所示的穩(wěn)態(tài)溫度場分布云圖。由圖9可以看出:V4模塊的芯片溫度最高,為55.105 ℃,1#—7#測溫點的溫度分別為39.46,45.58,48.82,48.72,45.38,39.30和45.11 ℃,與圖6所示的CFD仿真結(jié)果相比,溫度云圖分布規(guī)律一致且溫度值的差異小于1%,說明將由精細網(wǎng)格CFD模型得到的對流換熱系數(shù)分布插值到粗糙網(wǎng)格有限元模型是非常有效的,這可以說明此步驟使有限元模型的計算規(guī)模急劇減小,計算速度提高了數(shù)千倍。根據(jù)牛頓冷卻定律可知,通過對流換熱系數(shù)的引入,可以將風冷、水冷等所有的對流換熱過程以簡單的規(guī)律進行描述,因此,利用較少數(shù)量的有限元粗糙網(wǎng)格和對流換熱系數(shù)插值的方法進行瞬態(tài)熱仿真,可以推廣應用到其他類型的對流換熱應用中。

    圖8 有限元網(wǎng)格

    圖9 基于對流換熱系數(shù)插值的有限元模型溫度場分布云圖

    取降階階數(shù)為5,利用快速計算程序?qū)D8所示的有限元模型進行時間范圍為0~10 000 s的瞬態(tài)計算,可用2 min左右計算出如圖10所示的熱管散熱器的溫升曲線。由圖10可以看出:在0~1 000 s時間范圍內(nèi)溫度迅速上升,之后溫度上升的速度變緩,10 000 s時刻的溫度處于完全穩(wěn)定狀態(tài);將圖10所示的仿真結(jié)果與圖9所示的穩(wěn)態(tài)溫度場結(jié)果進行比較,發(fā)現(xiàn)最高溫度和1#—7#測溫點的溫度計算結(jié)果完全一致,說明快速計算程序所用方法的可行性與準確性。

    圖10 基于快速計算程序的溫升曲線

    在時間范圍相同的條件下,利用ANSYS軟件對圖8所示的有限元模型進行瞬態(tài)計算,則需要耗時約5 h才能得到如圖11所示完全一致的溫升曲線,說明筆者開發(fā)的快速計算程序的計算準確性和計算效率很高。

    4 電氣—熱聯(lián)合仿真

    4.1 電氣仿真

    由IGBT模塊的用戶手冊可知,其損耗與電流、管壓降及結(jié)溫有關(guān),因此可將導通電壓、開通損耗等曲線擬合成二次函數(shù)[14-15]。在選用Matlab/Simulink的SimPowerSystems模型庫建立地鐵車輛主電路的電氣仿真模型時,IGBT模塊的損耗等復雜工況的計算常使用線性化的參數(shù),而不是實際的隨二次函數(shù)變化的參數(shù),這不利于反映電流周期及暫態(tài)變化對功耗波動的影響,因此下面對IGBT模塊的損耗計算方法進行改進。

    IGBT模塊的損耗來自IGBT芯片和二極管芯片。其中,對于IGBT芯片主要計算其導通狀態(tài)下的損耗以及開通和關(guān)斷過程中的開關(guān)損耗,而忽略其關(guān)斷狀態(tài)下的損耗;對于二極管芯片主要計算其導通狀態(tài)下的和反向恢復過程中的損耗,而忽略其關(guān)斷狀態(tài)下的和開通過程中的損耗。

    以單個橋臂上管IGBT芯片為例,當IGBT模塊的電流大于0時,通過擬合的二次函數(shù)計算出電流通過該模塊時的開關(guān)損耗,以下管IGBT芯片脈寬調(diào)制下降沿判斷單次開關(guān)是否結(jié)束,開關(guān)結(jié)束時則累加而得IGBT模塊的開關(guān)損耗。通過擬合的二次函數(shù)計算電流流過模塊時的導通功率,再經(jīng)過積分得到IGBT模塊導通損耗。將IGBT模塊開關(guān)損耗和導通損耗進行累加,再計算出1個周期內(nèi)的損耗平均值,即為IGBT芯片的平均損耗。以單個橋臂上管二極管芯片為例,當IGBT模塊的電流小于0時,通過查找導通功耗表得出電流通過該模塊時的導通功率值,然后再積分得到導通損耗。當橋臂輸入電流大于0時,通過擬合的二次函數(shù)計算開關(guān)損耗,以下管二極管脈寬調(diào)制上升沿判斷上管二極管是否反向恢復,二極管反向恢復時則累加而得到二極管反向恢復損耗。將二極管導通功耗和反向恢復損耗進行累加,再計算出1個周期的損耗平均值,即為二極管芯片的平均損耗。其他IGBT模塊損耗的計算與上管相同,不做贅述。

    圖11為在電機的全速范圍(電機轉(zhuǎn)速在20 s內(nèi)從0加速至300 rad·s-1)內(nèi),利用主電路電氣仿真得到的電機電磁轉(zhuǎn)矩和IGBT模塊損耗隨時間變化的曲線。

    圖11 電機全速范圍內(nèi)的主電路電氣仿真結(jié)果

    4.2 瞬態(tài)熱仿真

    以圖11(b)所示的IGBT模塊損耗數(shù)據(jù)為快速計算程序的輸入,對圖8所示的有限元模型進行瞬態(tài)計算,可在7 min內(nèi)計算出如圖12所示的溫升曲線。由圖12可以看出:V2—V7的芯片在不同時刻產(chǎn)生變化的損耗,可以直接影響其溫度的響應;熱管散熱器安裝面上的各測溫點因稍遠離產(chǎn)生損耗的芯片,相應的溫度受到的影響較小;V1和V8距離產(chǎn)生損耗的芯片較遠,相應的溫度受到的影響最小,這除與距離熱源的遠近有關(guān)外,還與熱管散熱器底板具有較大的熱容量有關(guān)。

    圖12 全速范圍的溫升曲線

    以圖11(b)所示的IGBT模塊損耗數(shù)據(jù)為ANSYS軟件的輸入,對圖8所示的有限元模型進行瞬態(tài)計算,可得到與圖12完全一致的計算結(jié)果,但計算消耗的時間長達34.7 h,說明筆者開發(fā)的快速計算程序具有很高的計算效率和準確性,與電氣聯(lián)合仿真可為解決復雜工況的瞬態(tài)熱仿真問題提供有力支持。

    5 結(jié) 語

    針對某地鐵車輛牽引逆變器的IGBT模塊和熱管散熱器的溫升問題,首先在實驗室條件下完成了溫升測試,然后以試驗結(jié)果為輸入開展了CFD穩(wěn)態(tài)熱仿真、快速瞬態(tài)熱仿真和電氣—熱聯(lián)合仿真,試驗數(shù)據(jù)及不同的仿真方法驗證了筆者開發(fā)的快速計算方法可同時滿足計算效率與計算精度的要求,為解決復雜的瞬態(tài)熱仿真問題提供了1種新的途徑。

    [1]LUTZ J, SCHLANGENOTTO H, SCHEUERMANN U, et al. Semiconductor Power Devices, Physics, Characteristics, Reliability[M]. Berlin: Springer-Verlag Berlin Heidelberg, 2011, 380-409.

    [2]CIAPPA M. Selected Failure Mechanisms of Modern Power Modules[J]. Microelectronics Reliability, 2002 (42): 653-667.

    [3]HAMIDI A, KAUFMANN S, HERR E. Increased Lifetime of Wire Bond Connections for IGBT Power Modules[C]//16th IEEE Applied Power Electronic Conference and Exposition(APEC). Anaheim: IEEE, 2001: 1040-1044.

    [4]丁杰, 唐玉兔. 地鐵牽引變流器的熱管散熱器的數(shù)值模擬[J]. 大功率變流技術(shù), 2012 (5): 31-33.

    (DING Jie, TANG Yutu. Numerical Simulation of Heat-Pipe Radiator for Metro Vehicle Converter[J]. High Power Converter Technology, 2012(5): 31-33. in Chinese)

    [5]PERPINA X, JORDA X, VELLVEHI M, et al. Long-Term Reliability of Railway Power Inverters Cooled by Heat Pipe-Based Systems[J]. IEEE Tractions on Industrial Electronics, 2011, 58(7): 2662-2672.

    [6]唐強. IGBT元件熱管冷卻傳熱性能的實驗與數(shù)值研究[D]. 蘭州: 蘭州交通大學, 2013.

    (TANG Qiang. Experimental and Numerical Analysis on Heat Transfer Performance of Heat Pipe Cooling for IGBT Component[D]. Lanzhou: Lanzhou Jiaotong University, 2013. in Chinese)

    [7]YUN Chansu, MALBERTI P, CIAPPA M, et al. Thermal Component Model for Electrothermal Analysis of IGBT Module Systems[J]. IEEE Transactions on Advanced Packaging, 2001, 24(3): 401-406.

    [8]蔣耀林. 模型降階方法[M]. 北京: 科學出版社, 2010.

    [9]WILHELMUS H A Schilders, HENK A van der Vorst, JOOST Rommes. Model Order Reduction: Theory, Research Aspects and Applications[M]. Berlin: Springer, 2008.

    [10]TAMARA Bechtold, EVGENII B Rudnyi,JAN G Korvink. Dynamic Electro-Thermal Simulation of Microsystems-a Review[J]. Journal of Micromechanics and Microengineering, 2005, 15(11): 17-31.

    [11]RUDNYI E B, KORVINK J G. Model Order Reduction for Large Scale Engineering Models Developed in ANSYS[J]. Lecture Notes in Computer Science, 2006, 3732(1): 349-356.

    [12]YANG Y J, YU C C, SHEN K Y. MEMS Heat Transfer Arnoldi-Based Macromodels and the Study of Minimum Required Orders[J]. Tamkang Journal of Science and Engineering, 2005, 8(3): 185-190.

    [13]丁杰, 唐玉兔. 模型降階方法在瞬態(tài)熱仿真中的應用[J]. 機車電傳動, 2014(5): 51-55,65.

    (DING Jie, TANG Yutu. Application Research of Model Order Reduction Method in Transient Thermal Analysis[J]. Electric Drive for Locomotives, 2014 (5): 51-55,65. in Chinese)

    [14]白保東, 陳德志, 王鑫博. 逆變器IGBT損耗計算及冷卻裝置設(shè)計[J]. 電工技術(shù)學報, 2013, 28(8): 97-106.

    (BAI Baodong, CHEN Dezhi, WANG Xinbo. Loss Calculation of Inverter IGBT and Design of Cooling Device[J]. Transactions of China Electrotechnical Society, 2013, 28(8): 97-106. in Chinese)

    [15]楊寧, 宋術(shù)全, 李紅. 高速動車組輔助變流器箱體的熱仿真設(shè)計方法[J]. 中國鐵道科學, 2013, 34(3): 87-92.

    (YANG Ning, SONG Shuquan, LI Hong. Thermal Simulation Design Method for the Auxiliary Converter Case of High-Speed EMU[J]. China Railway Science, 2013, 34(3): 87-92. in Chinese)

    猜你喜歡
    風罩瞬態(tài)熱管
    旋挖鉆機導風罩開裂故障的分析與改進
    機械制造(2021年12期)2022-01-06 06:35:54
    高壓感應電動機斷電重啟時的瞬態(tài)仿真
    防爆電機(2020年3期)2020-11-06 09:07:36
    水電站地下廠房機墩溫度荷載對風罩內(nèi)力的影響分析
    集風罩對固定污染源廢氣檢測結(jié)果的影響
    山東化工(2020年9期)2020-06-01 06:56:48
    導熱冠軍——熱管(下)
    導熱冠軍——熱管(上)
    十億像素瞬態(tài)成像系統(tǒng)實時圖像拼接
    中國光學(2015年5期)2015-12-09 09:00:39
    U型換熱管試壓胎具設(shè)計
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    DC/DC變換器中的瞬態(tài)特性分析
    日韩 亚洲 欧美在线| av在线观看视频网站免费| 极品教师在线视频| 午夜久久久久精精品| 99国产精品一区二区三区| 一个人看视频在线观看www免费| 国产精品98久久久久久宅男小说| 婷婷六月久久综合丁香| www.www免费av| 男人的好看免费观看在线视频| 欧美另类亚洲清纯唯美| 午夜免费男女啪啪视频观看 | 又爽又黄无遮挡网站| 亚洲中文日韩欧美视频| 日本五十路高清| 九色成人免费人妻av| 亚洲五月婷婷丁香| 国产大屁股一区二区在线视频| 国产精品98久久久久久宅男小说| 亚洲综合色惰| 欧美乱妇无乱码| 国产综合懂色| 悠悠久久av| 97超视频在线观看视频| 久久久国产成人免费| 欧美xxxx性猛交bbbb| 国产色婷婷99| 欧美日本视频| 美女大奶头视频| 欧美性感艳星| 又爽又黄无遮挡网站| 国产亚洲精品av在线| 91午夜精品亚洲一区二区三区 | 精品熟女少妇八av免费久了| 桃色一区二区三区在线观看| 欧美+日韩+精品| av天堂在线播放| 亚洲无线在线观看| 国产美女午夜福利| 国产精品电影一区二区三区| 人人妻人人澡欧美一区二区| а√天堂www在线а√下载| 美女高潮的动态| 在线观看舔阴道视频| 此物有八面人人有两片| 免费观看人在逋| 久久99热6这里只有精品| 日韩欧美在线乱码| 色在线成人网| 啪啪无遮挡十八禁网站| 又黄又爽又刺激的免费视频.| 国产精品久久久久久久电影| 精品久久国产蜜桃| 很黄的视频免费| 午夜免费男女啪啪视频观看 | 免费在线观看影片大全网站| 九九热线精品视视频播放| 久99久视频精品免费| 啦啦啦韩国在线观看视频| 身体一侧抽搐| 女生性感内裤真人,穿戴方法视频| 99在线视频只有这里精品首页| 淫妇啪啪啪对白视频| 国产午夜福利久久久久久| 久久久久久久久久黄片| 国产成人aa在线观看| 免费一级毛片在线播放高清视频| av天堂在线播放| 欧美性感艳星| 日本在线视频免费播放| 亚洲中文字幕一区二区三区有码在线看| 美女高潮喷水抽搐中文字幕| 国产人妻一区二区三区在| 国产伦精品一区二区三区四那| 国内少妇人妻偷人精品xxx网站| 亚洲国产欧美人成| 又爽又黄无遮挡网站| 日韩国内少妇激情av| 丁香欧美五月| 桃红色精品国产亚洲av| 国产精品一区二区免费欧美| 麻豆久久精品国产亚洲av| 亚洲av成人av| 免费在线观看亚洲国产| 看十八女毛片水多多多| 亚洲自偷自拍三级| 国产真实乱freesex| 99久久精品一区二区三区| 久久热精品热| 一个人免费在线观看电影| 男人狂女人下面高潮的视频| 一级黄色大片毛片| 99热只有精品国产| 国产中年淑女户外野战色| 欧美日本视频| 国产欧美日韩一区二区三| a在线观看视频网站| 久久国产精品人妻蜜桃| 91久久精品电影网| 天美传媒精品一区二区| 色尼玛亚洲综合影院| 中文亚洲av片在线观看爽| 日韩欧美三级三区| 精品熟女少妇八av免费久了| 国产精品亚洲美女久久久| 亚洲成人中文字幕在线播放| 一本久久中文字幕| 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 狂野欧美白嫩少妇大欣赏| 亚洲成av人片在线播放无| 午夜影院日韩av| 亚洲欧美日韩卡通动漫| 成人美女网站在线观看视频| 2021天堂中文幕一二区在线观| 制服丝袜大香蕉在线| 免费黄网站久久成人精品 | 午夜激情福利司机影院| 亚洲第一区二区三区不卡| 亚洲自偷自拍三级| 久久久国产成人精品二区| 国产精品国产高清国产av| 丰满的人妻完整版| 三级毛片av免费| 国产av麻豆久久久久久久| 最好的美女福利视频网| 亚洲精华国产精华精| 亚洲第一电影网av| 搞女人的毛片| 精品人妻偷拍中文字幕| 无遮挡黄片免费观看| 高清日韩中文字幕在线| 精品一区二区免费观看| 久久久久久国产a免费观看| 在线观看美女被高潮喷水网站 | 国产色爽女视频免费观看| 搞女人的毛片| 91久久精品电影网| 婷婷六月久久综合丁香| 九九热线精品视视频播放| www.www免费av| 国产亚洲欧美98| 高潮久久久久久久久久久不卡| 757午夜福利合集在线观看| 又粗又爽又猛毛片免费看| 婷婷六月久久综合丁香| 亚洲国产精品成人综合色| 欧美黄色淫秽网站| 天天一区二区日本电影三级| 亚洲国产欧洲综合997久久,| 精品日产1卡2卡| 日韩中文字幕欧美一区二区| 国内毛片毛片毛片毛片毛片| 国产大屁股一区二区在线视频| 在线观看av片永久免费下载| 丝袜美腿在线中文| 在线国产一区二区在线| 色综合亚洲欧美另类图片| 99精品久久久久人妻精品| 国产成+人综合+亚洲专区| 亚洲美女视频黄频| 中文字幕精品亚洲无线码一区| 少妇熟女aⅴ在线视频| 国产免费av片在线观看野外av| 麻豆成人午夜福利视频| 久久久久免费精品人妻一区二区| 欧美黄色片欧美黄色片| 久久精品国产99精品国产亚洲性色| 黄色配什么色好看| 国产三级黄色录像| 成年女人毛片免费观看观看9| 精品国内亚洲2022精品成人| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 嫁个100分男人电影在线观看| 看片在线看免费视频| 一级a爱片免费观看的视频| 亚洲五月婷婷丁香| 又黄又爽又免费观看的视频| 女生性感内裤真人,穿戴方法视频| 91麻豆av在线| 此物有八面人人有两片| 禁无遮挡网站| 欧美精品啪啪一区二区三区| 亚洲无线在线观看| 蜜桃久久精品国产亚洲av| 精品一区二区三区人妻视频| 久久6这里有精品| 久久久久久久亚洲中文字幕 | av黄色大香蕉| 美女大奶头视频| 给我免费播放毛片高清在线观看| 中文资源天堂在线| 十八禁网站免费在线| 国产亚洲欧美98| 全区人妻精品视频| 亚洲av一区综合| 乱码一卡2卡4卡精品| av天堂在线播放| 亚洲精品久久国产高清桃花| 两性午夜刺激爽爽歪歪视频在线观看| 成年女人永久免费观看视频| 亚洲天堂国产精品一区在线| 丁香欧美五月| 99热6这里只有精品| 国语自产精品视频在线第100页| 久久久久久国产a免费观看| 欧美极品一区二区三区四区| 午夜福利在线在线| 好男人在线观看高清免费视频| 99国产精品一区二区三区| 亚洲18禁久久av| 久久久国产成人精品二区| 成年女人看的毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 乱码一卡2卡4卡精品| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久久久免费视频| 91麻豆av在线| 亚洲三级黄色毛片| 久久精品91蜜桃| 此物有八面人人有两片| 床上黄色一级片| 日韩高清综合在线| 精品无人区乱码1区二区| 亚洲在线观看片| 91狼人影院| 一区二区三区高清视频在线| 无遮挡黄片免费观看| 亚洲无线在线观看| 久久九九热精品免费| 国产精品久久久久久久久免 | 又黄又爽又免费观看的视频| 搡女人真爽免费视频火全软件 | 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久久人妻蜜臀av| 亚洲精品在线美女| 久久精品国产清高在天天线| 国模一区二区三区四区视频| 亚洲第一电影网av| 成人永久免费在线观看视频| 国产极品精品免费视频能看的| 国产高清激情床上av| 噜噜噜噜噜久久久久久91| 精品熟女少妇八av免费久了| 亚洲人成网站高清观看| 久久国产精品影院| 18禁黄网站禁片免费观看直播| 琪琪午夜伦伦电影理论片6080| 狂野欧美白嫩少妇大欣赏| 色哟哟哟哟哟哟| 高清日韩中文字幕在线| 中亚洲国语对白在线视频| 久久久久精品国产欧美久久久| av天堂中文字幕网| 日本精品一区二区三区蜜桃| 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 身体一侧抽搐| 久久久久久久精品吃奶| 精品人妻一区二区三区麻豆 | 久久久久九九精品影院| 人妻久久中文字幕网| 日韩国内少妇激情av| 色在线成人网| 最近中文字幕高清免费大全6 | 亚洲最大成人手机在线| 亚洲欧美日韩高清在线视频| 精品国产亚洲在线| 美女黄网站色视频| 免费观看的影片在线观看| av在线老鸭窝| 成人国产一区最新在线观看| 综合色av麻豆| 欧美成人免费av一区二区三区| 日本免费一区二区三区高清不卡| 真人一进一出gif抽搐免费| 婷婷精品国产亚洲av| 偷拍熟女少妇极品色| 精品国内亚洲2022精品成人| 亚洲国产欧洲综合997久久,| 尤物成人国产欧美一区二区三区| 久久亚洲精品不卡| 精品久久久久久成人av| 午夜福利欧美成人| 蜜桃亚洲精品一区二区三区| 日韩人妻高清精品专区| 99精品久久久久人妻精品| 九九在线视频观看精品| 国产三级黄色录像| 亚洲黑人精品在线| 国产精品爽爽va在线观看网站| 亚洲欧美激情综合另类| 99久久99久久久精品蜜桃| 国产高潮美女av| 精品久久国产蜜桃| 欧美日韩福利视频一区二区| 欧美另类亚洲清纯唯美| 又黄又爽又免费观看的视频| 午夜影院日韩av| 成人特级av手机在线观看| 色视频www国产| 久久久久免费精品人妻一区二区| 成人三级黄色视频| 亚洲五月天丁香| 精品久久久久久久久久久久久| 看十八女毛片水多多多| 99热只有精品国产| 久久久色成人| 国产精品人妻久久久久久| 久久久久亚洲av毛片大全| 一夜夜www| 午夜福利成人在线免费观看| 婷婷亚洲欧美| 国产成年人精品一区二区| 成年免费大片在线观看| 午夜免费成人在线视频| 免费在线观看亚洲国产| 人人妻,人人澡人人爽秒播| 一级作爱视频免费观看| 麻豆久久精品国产亚洲av| 男插女下体视频免费在线播放| 国语自产精品视频在线第100页| 亚洲av第一区精品v没综合| 国产亚洲精品av在线| 亚洲av日韩精品久久久久久密| 啦啦啦观看免费观看视频高清| 我的女老师完整版在线观看| 最近最新中文字幕大全电影3| 日本 av在线| 国产精品久久久久久久久免 | 国产黄色小视频在线观看| 午夜日韩欧美国产| 听说在线观看完整版免费高清| 午夜视频国产福利| 亚洲狠狠婷婷综合久久图片| 国产精品人妻久久久久久| 色综合站精品国产| 18+在线观看网站| 亚洲国产欧美人成| 18+在线观看网站| 成年女人永久免费观看视频| 中文字幕高清在线视频| 有码 亚洲区| 欧美潮喷喷水| 成年女人永久免费观看视频| 色综合婷婷激情| 成人特级av手机在线观看| 99久久精品一区二区三区| 亚洲内射少妇av| 看免费av毛片| 国产一区二区激情短视频| 成人一区二区视频在线观看| 午夜福利欧美成人| 亚洲人成电影免费在线| 亚洲在线观看片| 夜夜躁狠狠躁天天躁| 国产爱豆传媒在线观看| 90打野战视频偷拍视频| 中出人妻视频一区二区| 日本黄色视频三级网站网址| 91久久精品国产一区二区成人| 亚洲国产精品久久男人天堂| 国产精品一区二区三区四区免费观看 | 久久热精品热| 日韩欧美免费精品| 99热这里只有是精品在线观看 | 亚洲av中文字字幕乱码综合| 午夜免费成人在线视频| 亚洲国产精品成人综合色| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人久久性| 偷拍熟女少妇极品色| 久久久国产成人精品二区| 伦理电影大哥的女人| 精品久久国产蜜桃| 亚洲欧美激情综合另类| 黄色配什么色好看| 免费电影在线观看免费观看| 一本一本综合久久| 韩国av一区二区三区四区| www.色视频.com| 欧美日韩福利视频一区二区| 国产伦在线观看视频一区| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品久久久com| 国内精品一区二区在线观看| 亚洲无线在线观看| 亚洲人成网站高清观看| 久久亚洲精品不卡| aaaaa片日本免费| 小说图片视频综合网站| 床上黄色一级片| 午夜福利成人在线免费观看| 国产午夜福利久久久久久| 色视频www国产| 我的老师免费观看完整版| 青草久久国产| 欧美区成人在线视频| 久久久久久久精品吃奶| 在线国产一区二区在线| 亚洲自偷自拍三级| 亚洲美女搞黄在线观看 | 非洲黑人性xxxx精品又粗又长| 老熟妇乱子伦视频在线观看| 精品久久久久久久久av| 国产成人欧美在线观看| 在线天堂最新版资源| 中文字幕熟女人妻在线| 蜜桃久久精品国产亚洲av| 日本黄大片高清| 日韩欧美免费精品| av在线老鸭窝| 成人三级黄色视频| 久久性视频一级片| 日本熟妇午夜| 久久人人爽人人爽人人片va | 国语自产精品视频在线第100页| 亚洲自偷自拍三级| 如何舔出高潮| 男人舔奶头视频| 99riav亚洲国产免费| .国产精品久久| 免费观看人在逋| 波多野结衣巨乳人妻| or卡值多少钱| 亚洲最大成人中文| 97热精品久久久久久| 国产极品精品免费视频能看的| 99精品久久久久人妻精品| 波多野结衣高清无吗| 成人鲁丝片一二三区免费| 色噜噜av男人的天堂激情| www.熟女人妻精品国产| 天堂√8在线中文| 国产精品电影一区二区三区| 免费看美女性在线毛片视频| 国产精品1区2区在线观看.| 十八禁人妻一区二区| 日本黄大片高清| 99热6这里只有精品| 少妇人妻一区二区三区视频| 久久久久精品国产欧美久久久| 五月玫瑰六月丁香| 国产精品不卡视频一区二区 | 国产一区二区三区视频了| 欧美xxxx性猛交bbbb| 日本与韩国留学比较| 99热这里只有精品一区| 日本精品一区二区三区蜜桃| 国产精品不卡视频一区二区 | 1000部很黄的大片| 麻豆av噜噜一区二区三区| 在线观看一区二区三区| 国产一级毛片七仙女欲春2| 乱人视频在线观看| 欧美潮喷喷水| 欧美黑人欧美精品刺激| 国产69精品久久久久777片| 国内精品美女久久久久久| 欧美在线一区亚洲| netflix在线观看网站| 男女床上黄色一级片免费看| 国产久久久一区二区三区| 国产高清视频在线播放一区| 日本a在线网址| 久久久久久久久中文| 九色成人免费人妻av| 亚洲久久久久久中文字幕| 亚洲真实伦在线观看| 给我免费播放毛片高清在线观看| 中文资源天堂在线| 九色成人免费人妻av| 久久久久久久亚洲中文字幕 | 日韩有码中文字幕| 婷婷精品国产亚洲av| 欧美在线一区亚洲| 少妇高潮的动态图| 久久6这里有精品| 成人高潮视频无遮挡免费网站| 我要搜黄色片| 欧美激情国产日韩精品一区| 欧美性感艳星| 国产精品99久久久久久久久| 亚洲性夜色夜夜综合| 成人欧美大片| 黄片小视频在线播放| 亚洲av电影不卡..在线观看| 亚洲精品日韩av片在线观看| 99久久精品国产亚洲精品| 国产一级毛片七仙女欲春2| 亚洲av一区综合| 看片在线看免费视频| 床上黄色一级片| 99精品在免费线老司机午夜| 日韩欧美三级三区| 乱码一卡2卡4卡精品| 国产91精品成人一区二区三区| 日韩欧美在线二视频| 亚洲片人在线观看| 国产精品久久电影中文字幕| 最新中文字幕久久久久| 校园春色视频在线观看| 午夜福利欧美成人| 欧美极品一区二区三区四区| 午夜福利欧美成人| 成人亚洲精品av一区二区| 色噜噜av男人的天堂激情| 欧美另类亚洲清纯唯美| 少妇丰满av| 国产精品久久电影中文字幕| 国产伦精品一区二区三区视频9| 日本三级黄在线观看| 国产精品久久久久久精品电影| 男女那种视频在线观看| 亚州av有码| 久久久精品大字幕| 一级毛片久久久久久久久女| 成人特级黄色片久久久久久久| xxxwww97欧美| 久久中文看片网| 日日摸夜夜添夜夜添av毛片 | 又爽又黄a免费视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美不卡视频在线免费观看| 能在线免费观看的黄片| 精品国产三级普通话版| 小蜜桃在线观看免费完整版高清| 老鸭窝网址在线观看| 国产爱豆传媒在线观看| 亚洲美女搞黄在线观看 | 88av欧美| 99国产精品一区二区三区| 在现免费观看毛片| 18禁黄网站禁片免费观看直播| 欧美一级a爱片免费观看看| 国产黄色小视频在线观看| 999久久久精品免费观看国产| 一级av片app| 亚洲熟妇中文字幕五十中出| 日本熟妇午夜| 国产人妻一区二区三区在| avwww免费| 天堂av国产一区二区熟女人妻| 久久亚洲真实| 欧美又色又爽又黄视频| 国产熟女xx| 久久久久九九精品影院| 国产真实伦视频高清在线观看 | 每晚都被弄得嗷嗷叫到高潮| 欧美成人免费av一区二区三区| 亚洲熟妇熟女久久| 国语自产精品视频在线第100页| 狂野欧美白嫩少妇大欣赏| 国产蜜桃级精品一区二区三区| 中文字幕av在线有码专区| 18美女黄网站色大片免费观看| 国产av在哪里看| 国产色爽女视频免费观看| 亚洲av一区综合| 色哟哟哟哟哟哟| 欧美极品一区二区三区四区| 国产精品日韩av在线免费观看| 亚洲经典国产精华液单 | 久久精品国产自在天天线| 99久久精品一区二区三区| 亚洲精华国产精华精| av在线观看视频网站免费| 国产在线精品亚洲第一网站| 男人舔奶头视频| 日韩欧美国产在线观看| 日本免费a在线| 亚洲成人久久性| 亚洲不卡免费看| 美女高潮喷水抽搐中文字幕| 国产伦精品一区二区三区四那| 偷拍熟女少妇极品色| 欧美日韩国产亚洲二区| 18禁在线播放成人免费| 国产精品99久久久久久久久| 天天躁日日操中文字幕| 在线观看免费视频日本深夜| 成人午夜高清在线视频| 免费无遮挡裸体视频| 国产精品久久视频播放| 波野结衣二区三区在线| 日本 欧美在线| 亚洲精品一卡2卡三卡4卡5卡| 啪啪无遮挡十八禁网站| 日本黄色视频三级网站网址| 国产精品嫩草影院av在线观看 | 免费看光身美女| 欧美成人一区二区免费高清观看| 国内精品美女久久久久久| 国内精品久久久久久久电影| 美女大奶头视频| 日韩欧美国产在线观看| 亚洲精品粉嫩美女一区| 亚洲一区二区三区不卡视频| 69av精品久久久久久| 全区人妻精品视频| 三级毛片av免费| 婷婷色综合大香蕉| 久久人妻av系列| 久久精品影院6| 色5月婷婷丁香| 欧美xxxx黑人xx丫x性爽| 亚洲成人久久性| 国产精品亚洲av一区麻豆| 99久久精品热视频| 免费人成在线观看视频色| 美女xxoo啪啪120秒动态图 | 欧美xxxx黑人xx丫x性爽| 精品国产亚洲在线| 久久久国产成人精品二区|