• <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)特性分析
    日本一二三区视频观看| 久久久a久久爽久久v久久| av在线亚洲专区| 成人毛片a级毛片在线播放| 国产一区二区在线观看日韩| 狠狠狠狠99中文字幕| 韩国av在线不卡| 欧美潮喷喷水| 成年女人看的毛片在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产高清国产精品国产三级 | 国产午夜精品久久久久久一区二区三区| 长腿黑丝高跟| 日本欧美国产在线视频| 亚洲丝袜综合中文字幕| 久久久亚洲精品成人影院| 国产一区亚洲一区在线观看| 亚洲av免费高清在线观看| 国产综合懂色| 亚洲电影在线观看av| 九九久久精品国产亚洲av麻豆| 69av精品久久久久久| 婷婷色av中文字幕| 国产毛片a区久久久久| 国产 一区精品| 色哟哟·www| 国产淫片久久久久久久久| 又黄又爽又刺激的免费视频.| .国产精品久久| 久久精品久久久久久噜噜老黄 | 成年女人看的毛片在线观看| 亚洲国产欧美人成| 久久久久性生活片| 搞女人的毛片| 午夜a级毛片| 岛国在线免费视频观看| 欧美另类亚洲清纯唯美| 久久久久国产网址| 亚洲av男天堂| 亚洲内射少妇av| 亚洲欧美清纯卡通| 99热精品在线国产| 日韩精品青青久久久久久| 久久久久久国产a免费观看| 亚洲欧美中文字幕日韩二区| 国产精品野战在线观看| 欧美又色又爽又黄视频| 18禁在线播放成人免费| 欧美极品一区二区三区四区| 精品不卡国产一区二区三区| 国产男人的电影天堂91| 国产麻豆成人av免费视频| 看片在线看免费视频| 两个人的视频大全免费| 国产国拍精品亚洲av在线观看| 一本一本综合久久| 看黄色毛片网站| 日韩av在线免费看完整版不卡| 床上黄色一级片| 2021少妇久久久久久久久久久| videos熟女内射| 九九在线视频观看精品| 精品久久久久久久久av| 久久精品久久久久久久性| 欧美日韩国产亚洲二区| 少妇人妻一区二区三区视频| 男女视频在线观看网站免费| 亚洲人成网站高清观看| 亚洲精品亚洲一区二区| 日日干狠狠操夜夜爽| 极品教师在线视频| 国产熟女欧美一区二区| 国产精品三级大全| 国产中年淑女户外野战色| 波多野结衣巨乳人妻| 欧美又色又爽又黄视频| 少妇猛男粗大的猛烈进出视频 | 久久草成人影院| 日韩国内少妇激情av| 深夜a级毛片| 一夜夜www| 久久久国产成人精品二区| 日韩国内少妇激情av| eeuss影院久久| 亚洲aⅴ乱码一区二区在线播放| 国产91av在线免费观看| 久久久久久久久久久免费av| 黑人高潮一二区| 久久久国产成人免费| 国产欧美日韩精品一区二区| 91aial.com中文字幕在线观看| 日本爱情动作片www.在线观看| 极品教师在线视频| 国产老妇女一区| 国产一区亚洲一区在线观看| 桃色一区二区三区在线观看| 岛国毛片在线播放| 精品一区二区免费观看| 国产真实乱freesex| 少妇的逼水好多| 汤姆久久久久久久影院中文字幕 | 午夜免费激情av| 久久久久久久亚洲中文字幕| av在线观看视频网站免费| 一级毛片久久久久久久久女| 女的被弄到高潮叫床怎么办| 欧美高清成人免费视频www| 男的添女的下面高潮视频| 国内精品美女久久久久久| 国产精品不卡视频一区二区| 91久久精品电影网| 高清日韩中文字幕在线| 美女内射精品一级片tv| 亚洲国产精品久久男人天堂| 日本熟妇午夜| 国产亚洲5aaaaa淫片| 国产视频首页在线观看| 中文字幕精品亚洲无线码一区| 中文字幕人妻熟人妻熟丝袜美| 亚洲成av人片在线播放无| ponron亚洲| 精品久久久久久久久av| 能在线免费观看的黄片| 小蜜桃在线观看免费完整版高清| 亚洲精品aⅴ在线观看| 国产男人的电影天堂91| 99国产精品一区二区蜜桃av| 99热这里只有是精品在线观看| 欧美变态另类bdsm刘玥| 一区二区三区高清视频在线| 国产欧美另类精品又又久久亚洲欧美| 一区二区三区四区激情视频| 男人狂女人下面高潮的视频| 天堂av国产一区二区熟女人妻| 精品不卡国产一区二区三区| 日本色播在线视频| 成年女人永久免费观看视频| 99久久无色码亚洲精品果冻| 少妇熟女欧美另类| 黄色欧美视频在线观看| 亚洲av二区三区四区| 神马国产精品三级电影在线观看| 26uuu在线亚洲综合色| av在线天堂中文字幕| 国产亚洲av嫩草精品影院| 亚洲av不卡在线观看| 久久亚洲国产成人精品v| 亚洲欧美精品自产自拍| 搡女人真爽免费视频火全软件| 一级二级三级毛片免费看| 久久久久网色| 亚洲av成人精品一区久久| 黄色配什么色好看| 亚洲国产欧洲综合997久久,| 一卡2卡三卡四卡精品乱码亚洲| 99视频精品全部免费 在线| 亚洲欧洲日产国产| 别揉我奶头 嗯啊视频| 大话2 男鬼变身卡| 久久韩国三级中文字幕| 99热精品在线国产| 久久婷婷人人爽人人干人人爱| 大香蕉97超碰在线| 少妇的逼水好多| 在线观看一区二区三区| 简卡轻食公司| 欧美日韩综合久久久久久| 久久精品影院6| 色视频www国产| av天堂中文字幕网| or卡值多少钱| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 成人综合一区亚洲| 97人妻精品一区二区三区麻豆| 男女视频在线观看网站免费| 国产亚洲精品av在线| 久久这里只有精品中国| 午夜激情福利司机影院| 国产高清不卡午夜福利| 欧美性猛交黑人性爽| 中文字幕制服av| 国产免费福利视频在线观看| 一二三四中文在线观看免费高清| 赤兔流量卡办理| 久久久久精品久久久久真实原创| 中文字幕精品亚洲无线码一区| 国产白丝娇喘喷水9色精品| 又粗又硬又长又爽又黄的视频| 色视频www国产| 九九在线视频观看精品| 嫩草影院新地址| 国产高清三级在线| 少妇被粗大猛烈的视频| 51国产日韩欧美| 欧美激情久久久久久爽电影| 综合色av麻豆| 国产日韩欧美在线精品| 日本午夜av视频| 99久久精品国产国产毛片| 99久久精品一区二区三区| 男人舔女人下体高潮全视频| 国产精品国产三级国产专区5o | 综合色丁香网| 黄色欧美视频在线观看| 日本wwww免费看| 国产人妻一区二区三区在| 不卡视频在线观看欧美| 亚洲精品国产成人久久av| av线在线观看网站| 国产 一区精品| 欧美色视频一区免费| 国内精品一区二区在线观看| 狂野欧美激情性xxxx在线观看| 欧美成人午夜免费资源| 婷婷色麻豆天堂久久 | 久久久久久九九精品二区国产| 亚洲aⅴ乱码一区二区在线播放| 精品人妻一区二区三区麻豆| 国产探花极品一区二区| 国产一区二区在线观看日韩| 狂野欧美白嫩少妇大欣赏| 白带黄色成豆腐渣| 国产精品女同一区二区软件| 国产淫语在线视频| 国产精品,欧美在线| 国产69精品久久久久777片| 免费观看在线日韩| 一区二区三区免费毛片| 欧美激情国产日韩精品一区| 国产视频首页在线观看| 中文字幕久久专区| 美女黄网站色视频| 国产精品久久电影中文字幕| 中文欧美无线码| 精品欧美国产一区二区三| 看免费成人av毛片| 中国国产av一级| 最近最新中文字幕大全电影3| 精品人妻视频免费看| 一本一本综合久久| 全区人妻精品视频| 91在线精品国自产拍蜜月| 久久这里有精品视频免费| 国产不卡一卡二| 狂野欧美白嫩少妇大欣赏| 精品国产露脸久久av麻豆 | 日韩高清综合在线| 网址你懂的国产日韩在线| 狂野欧美激情性xxxx在线观看| 真实男女啪啪啪动态图| 精品一区二区免费观看| 国产精品久久久久久久久免| 亚洲精品亚洲一区二区| 亚洲国产精品成人久久小说| 成年女人看的毛片在线观看| 日韩欧美精品免费久久| 欧美不卡视频在线免费观看| 久久99热这里只频精品6学生 | 日韩大片免费观看网站 | 婷婷色麻豆天堂久久 | 精品人妻视频免费看| 成人性生交大片免费视频hd| 日本黄大片高清| 国产成人a∨麻豆精品| 中文资源天堂在线| 久久久久久久午夜电影| 亚洲在线自拍视频| 我的女老师完整版在线观看| 精品久久久久久成人av| www日本黄色视频网| 非洲黑人性xxxx精品又粗又长| 国产午夜福利久久久久久| 丝袜美腿在线中文| 不卡视频在线观看欧美| 视频中文字幕在线观看| 国产伦理片在线播放av一区| 九九热线精品视视频播放| 男女那种视频在线观看| 精品一区二区免费观看| 色尼玛亚洲综合影院| 美女大奶头视频| 午夜激情欧美在线| 久久精品人妻少妇| 乱系列少妇在线播放| 亚洲中文字幕日韩| 日韩欧美精品免费久久| 人人妻人人澡人人爽人人夜夜 | 国产精华一区二区三区| 在线a可以看的网站| 一个人免费在线观看电影| 精品久久久久久久人妻蜜臀av| 国产私拍福利视频在线观看| 尤物成人国产欧美一区二区三区| 亚洲美女视频黄频| 亚洲精品乱码久久久v下载方式| 亚洲乱码一区二区免费版| 精品一区二区三区视频在线| 国产精品一区二区性色av| 能在线免费看毛片的网站| 大香蕉久久网| 免费搜索国产男女视频| 三级经典国产精品| 色播亚洲综合网| 色5月婷婷丁香| 久久久久久久午夜电影| 国产精品野战在线观看| 亚洲国产欧美人成| 精品久久久久久久人妻蜜臀av| 亚洲av福利一区| 直男gayav资源| 色噜噜av男人的天堂激情| 欧美3d第一页| 亚洲国产精品sss在线观看| 麻豆乱淫一区二区| 亚洲真实伦在线观看| 午夜爱爱视频在线播放| 久久热精品热| 少妇人妻精品综合一区二区| 欧美又色又爽又黄视频| 国产视频内射| av播播在线观看一区| 国产中年淑女户外野战色| 国产爱豆传媒在线观看| 欧美xxxx性猛交bbbb| 亚洲av日韩在线播放| 性色avwww在线观看| h日本视频在线播放| 久久精品久久精品一区二区三区| 黄色配什么色好看| 99久久无色码亚洲精品果冻| 波野结衣二区三区在线| 亚洲高清免费不卡视频| 成人亚洲欧美一区二区av| 天天躁日日操中文字幕| 亚洲av二区三区四区| 国产免费视频播放在线视频 | 国内揄拍国产精品人妻在线| 日本三级黄在线观看| 日韩精品有码人妻一区| 内地一区二区视频在线| 国内揄拍国产精品人妻在线| 国产精品永久免费网站| 国产高清三级在线| 亚洲精品成人久久久久久| 成人鲁丝片一二三区免费| 波多野结衣高清无吗| 日韩 亚洲 欧美在线| 欧美一区二区精品小视频在线| 麻豆成人午夜福利视频| 日韩av不卡免费在线播放| 精品人妻熟女av久视频| 亚洲精品久久久久久婷婷小说 | 国产成人精品久久久久久| videossex国产| 人人妻人人澡人人爽人人夜夜 | 亚洲成人久久爱视频| 国产精品一区二区性色av| 国产精品熟女久久久久浪| 免费观看在线日韩| av女优亚洲男人天堂| 麻豆久久精品国产亚洲av| 日韩亚洲欧美综合| 22中文网久久字幕| 在线免费观看不下载黄p国产| 身体一侧抽搐| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 麻豆乱淫一区二区| 亚洲国产精品sss在线观看| 亚洲精品456在线播放app| 亚洲成人中文字幕在线播放| 高清毛片免费看| 亚洲怡红院男人天堂| 亚洲久久久久久中文字幕| av播播在线观看一区| 毛片女人毛片| 人人妻人人澡人人爽人人夜夜 | 舔av片在线| 99久久精品一区二区三区| 国产精品日韩av在线免费观看| av专区在线播放| 久久久久久大精品| 高清视频免费观看一区二区 | 99热这里只有精品一区| 色综合亚洲欧美另类图片| 国产精品永久免费网站| 亚洲av男天堂| 精品国内亚洲2022精品成人| 日韩强制内射视频| 色5月婷婷丁香| 天美传媒精品一区二区| 网址你懂的国产日韩在线| videossex国产| 全区人妻精品视频| 久久久色成人| 久久精品人妻少妇| 久久久国产成人精品二区| www日本黄色视频网| 成年av动漫网址| 日韩制服骚丝袜av| 国产日韩欧美在线精品| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 免费观看性生交大片5| 伦理电影大哥的女人| 18禁在线无遮挡免费观看视频| 精品国产露脸久久av麻豆 | 一级av片app| 色噜噜av男人的天堂激情| 亚洲国产精品久久男人天堂| 国产成年人精品一区二区| 啦啦啦韩国在线观看视频| 男女国产视频网站| 伊人久久精品亚洲午夜| 18禁在线无遮挡免费观看视频| 在线天堂最新版资源| 午夜精品在线福利| 亚洲精品日韩在线中文字幕| 久久久国产成人精品二区| 26uuu在线亚洲综合色| 精品久久久久久久久久久久久| 黑人高潮一二区| 99久久精品一区二区三区| 舔av片在线| 免费看光身美女| 亚洲精品影视一区二区三区av| 插逼视频在线观看| 水蜜桃什么品种好| 国产真实伦视频高清在线观看| 特大巨黑吊av在线直播| 亚洲伊人久久精品综合 | 精品一区二区免费观看| 青春草视频在线免费观看| 国产高清国产精品国产三级 | av女优亚洲男人天堂| 男人舔奶头视频| 91精品国产九色| 日韩精品有码人妻一区| 成人美女网站在线观看视频| 91久久精品国产一区二区三区| 内射极品少妇av片p| 六月丁香七月| 午夜福利网站1000一区二区三区| 日本黄色片子视频| 草草在线视频免费看| 国产 一区 欧美 日韩| 国产男人的电影天堂91| 尾随美女入室| av卡一久久| av.在线天堂| 精品一区二区三区视频在线| 亚洲欧美精品专区久久| 中文字幕亚洲精品专区| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 亚洲欧美日韩卡通动漫| 夜夜看夜夜爽夜夜摸| 夫妻性生交免费视频一级片| 亚洲精品国产av成人精品| 22中文网久久字幕| 亚洲欧美精品综合久久99| 99热网站在线观看| 国产午夜精品一二区理论片| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩一区二区视频在线观看视频在线 | 成年女人看的毛片在线观看| 亚洲人与动物交配视频| 欧美bdsm另类| 欧美97在线视频| 听说在线观看完整版免费高清| 欧美不卡视频在线免费观看| 99在线视频只有这里精品首页| 蜜臀久久99精品久久宅男| 精品99又大又爽又粗少妇毛片| 国产精品国产三级国产av玫瑰| 成人午夜精彩视频在线观看| 久久6这里有精品| 九九在线视频观看精品| 一级黄片播放器| 最新中文字幕久久久久| 成人综合一区亚洲| 日韩成人av中文字幕在线观看| 精品不卡国产一区二区三区| 免费黄色在线免费观看| 久久久久精品久久久久真实原创| 99热精品在线国产| av在线播放精品| 欧美三级亚洲精品| 日本黄色片子视频| 最近最新中文字幕免费大全7| 成人欧美大片| 狠狠狠狠99中文字幕| 天堂中文最新版在线下载 | 最近中文字幕高清免费大全6| 插逼视频在线观看| 国产视频内射| 我的女老师完整版在线观看| 91精品一卡2卡3卡4卡| 美女内射精品一级片tv| 日韩av在线免费看完整版不卡| 久久久久久久国产电影| 亚洲av成人av| 亚洲在久久综合| 毛片女人毛片| 精品人妻视频免费看| 久久人妻av系列| 成人漫画全彩无遮挡| 最近中文字幕2019免费版| 狂野欧美白嫩少妇大欣赏| 久久精品久久久久久噜噜老黄 | av在线播放精品| 欧美激情在线99| 国产精品久久久久久久电影| 国产视频内射| 最近最新中文字幕免费大全7| kizo精华| 久久久久九九精品影院| 久久99蜜桃精品久久| 国产又黄又爽又无遮挡在线| 免费大片18禁| 永久网站在线| 一级av片app| 亚洲国产欧美在线一区| 国产成人freesex在线| 在现免费观看毛片| 亚洲国产精品sss在线观看| 又黄又爽又刺激的免费视频.| 欧美一区二区亚洲| 日本-黄色视频高清免费观看| 桃色一区二区三区在线观看| 日韩欧美三级三区| 久久精品国产亚洲av涩爱| 日本欧美国产在线视频| 亚洲最大成人av| 国产黄片美女视频| 一级黄片播放器| 日日撸夜夜添| 欧美+日韩+精品| 男插女下体视频免费在线播放| 蜜臀久久99精品久久宅男| 中文字幕人妻熟人妻熟丝袜美| 国产精品精品国产色婷婷| 美女脱内裤让男人舔精品视频| 国产精品一二三区在线看| 国产精品国产三级国产av玫瑰| 伦精品一区二区三区| 午夜精品在线福利| 少妇被粗大猛烈的视频| a级一级毛片免费在线观看| 嫩草影院入口| av在线观看视频网站免费| 国产色爽女视频免费观看| 天美传媒精品一区二区| 成人毛片a级毛片在线播放| 精品久久久久久成人av| 91久久精品国产一区二区三区| 最近中文字幕2019免费版| 桃色一区二区三区在线观看| 国产一区亚洲一区在线观看| 欧美激情久久久久久爽电影| 国产麻豆成人av免费视频| 色综合亚洲欧美另类图片| 久久久久免费精品人妻一区二区| 2021少妇久久久久久久久久久| 亚洲精品成人久久久久久| av专区在线播放| av国产久精品久网站免费入址| 欧美日本亚洲视频在线播放| 夜夜爽夜夜爽视频| 卡戴珊不雅视频在线播放| 欧美又色又爽又黄视频| 亚洲人成网站高清观看| 99国产精品一区二区蜜桃av| 寂寞人妻少妇视频99o| 国产91av在线免费观看| 日韩强制内射视频| 国产午夜福利久久久久久| 午夜精品国产一区二区电影 | 中文欧美无线码| 免费观看精品视频网站| 六月丁香七月| 亚洲自拍偷在线| 青春草视频在线免费观看| 亚洲精品影视一区二区三区av| 人人妻人人澡欧美一区二区| 午夜久久久久精精品| 嫩草影院新地址| 久久精品人妻少妇| 真实男女啪啪啪动态图| 日日摸夜夜添夜夜爱| 午夜福利在线在线| 国产日韩欧美在线精品| 国产精品久久久久久av不卡| 日本与韩国留学比较| 天堂网av新在线| 91精品伊人久久大香线蕉| 国产三级在线视频| 一级毛片aaaaaa免费看小| 亚洲国产最新在线播放| 亚洲成人中文字幕在线播放| 波多野结衣巨乳人妻| 禁无遮挡网站| 欧美bdsm另类| 国产伦理片在线播放av一区| 精品一区二区三区视频在线| 美女国产视频在线观看| 久久鲁丝午夜福利片| 51国产日韩欧美| 国产精品人妻久久久影院| 亚洲经典国产精华液单| 免费看av在线观看网站| 97在线视频观看| 国产极品精品免费视频能看的| 亚洲不卡免费看|