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

    基于內(nèi)置溫度傳感器的碳化硅功率模塊結(jié)溫在線提取方法

    2021-07-01 05:55:16李海鵬苗軼如陳常樂黃守道
    電工技術(shù)學(xué)報 2021年12期
    關(guān)鍵詞:模型

    劉 平 李海鵬 苗軼如 陳常樂 黃守道

    基于內(nèi)置溫度傳感器的碳化硅功率模塊結(jié)溫在線提取方法

    劉 平 李海鵬 苗軼如 陳常樂 黃守道

    (湖南大學(xué)電氣與信息工程學(xué)院 長沙 410082)

    碳化硅(SiC)MOSFET具有耐壓高、開關(guān)速度快、導(dǎo)通損耗低等優(yōu)點,將越來越廣泛地應(yīng)用于高效、高功率密度場合。在這些應(yīng)用場合中,SiC MOSFET面臨著嚴(yán)峻的可靠性考驗,而結(jié)溫的在線準(zhǔn)確提取是實現(xiàn)器件壽命預(yù)測和可靠性評估的重要基礎(chǔ)。該文提出一種基于功率模塊內(nèi)置負(fù)溫度系數(shù)(NTC)溫度傳感器的器件結(jié)溫在線提取方法。首先建立考慮多芯片熱耦合效應(yīng)的內(nèi)置溫度傳感器至功率芯片的熱網(wǎng)絡(luò)模型,并建立SiC MOSFET的損耗快速計算方法;通過有限元仿真提取熱網(wǎng)絡(luò)模型的阻抗參數(shù),并驗證該熱網(wǎng)絡(luò)參數(shù)在不同邊界條件下的穩(wěn)定性。仿真和實驗結(jié)果表明,所提出的結(jié)溫在線估計方法能夠準(zhǔn)確地獲得器件的動態(tài)結(jié)溫,且熱網(wǎng)絡(luò)模型參數(shù)不受環(huán)境溫度、散熱條件等邊界條件變化的影響,適用于實際任務(wù)剖面下的結(jié)溫監(jiān)測與壽命 預(yù)測。

    負(fù)溫度系數(shù)溫度傳感器 SiC MOSFET 熱網(wǎng)絡(luò)模型 芯片溫度

    0 引言

    碳化硅(Silicon Carbide, SiC)MOSFET作為典型的寬禁帶半導(dǎo)體器件,具有高開關(guān)速度、低導(dǎo)通損耗和高熱導(dǎo)率等優(yōu)點,可以減小電感、電容和散熱器等元器件的體積,大幅提升系統(tǒng)的功率密度,非常適用于高效、高功率密度需求場合[1],近年來已逐步應(yīng)用于電氣化交通領(lǐng)域[2-3]。隨著SiC器件的推廣應(yīng)用,其健康管理和可靠性評估也越來越受到學(xué)術(shù)界和工業(yè)界的關(guān)注[4]。由電力電子系統(tǒng)可靠性調(diào)研報告可知,電力電子變流器中約有34%的失效是由功率器件引起[5],其中由熱應(yīng)力引起的器件失效占比高達(dá)55%[6]。根據(jù)試驗統(tǒng)計數(shù)據(jù)和失效機理分析,功率器件在失效前所經(jīng)歷的溫度循環(huán)周期主要由結(jié)溫擺幅Δj、最高結(jié)溫j.max和平均結(jié)溫j.ave等因素決定[7],而SiC功率模塊具有功率密度高、封裝緊湊和耐高溫等特點,其運行結(jié)溫更高且結(jié)溫波動更大。因此,器件結(jié)溫的精準(zhǔn)在線提取是其壽命預(yù)測、健康管理和可靠性評估的基礎(chǔ)[8]。

    目前,功率器件結(jié)溫提取方法主要有光學(xué)測量法、物理接觸法、熱網(wǎng)絡(luò)模型法和熱敏電參數(shù)法四類[9]。其中,熱敏電參數(shù)法和熱網(wǎng)絡(luò)模型法能夠滿足非侵入、可在線集成的要求。熱敏電參數(shù)法是通過測量對結(jié)溫敏感性強的外部電氣參數(shù)來逆向預(yù)估芯片的溫度,被認(rèn)為是極具應(yīng)用前景的結(jié)溫檢測方法[9]。然而,現(xiàn)階段的熱敏電參數(shù)法也還存在相應(yīng)的缺陷,如影響器件正常的運行[10-11]、測量參數(shù)靈敏度過低、極易受到干擾[12-14]等,應(yīng)用于復(fù)雜工況時會受到很大的限制,且目前的研究也主要聚焦于熱敏電參數(shù)的溫度特性和離線校正,較少涉及實際工況下結(jié)溫的在線測量和SiC器件的結(jié)溫提取研究[9, 15]。熱網(wǎng)絡(luò)模型法[16-18]是基于器件的實時損耗與瞬態(tài)熱阻抗網(wǎng)絡(luò)模型,通過實時計算或離線查表估算芯片結(jié)溫及其變化趨勢,優(yōu)點是無需外部測量電路,在獲得準(zhǔn)確而穩(wěn)定的熱網(wǎng)絡(luò)模型參數(shù)和損耗計算值的條件下,易于準(zhǔn)確地在線提取器件結(jié)溫。然而,傳統(tǒng)的包含功率模塊外部條件的熱網(wǎng)絡(luò)模型具有很大的局限性:以模塊基板溫度為參考點時,需要在模塊基板開槽埋設(shè)熱電偶,熱電偶的埋設(shè)位置也會對最終模塊內(nèi)部芯片熱狀態(tài)估計結(jié)果造成很大影響[19];以環(huán)境溫度為參考點時,則需要獲得導(dǎo)熱硅脂和散熱器的熱阻抗,而導(dǎo)熱硅脂的老化會導(dǎo)致硅脂層的熱阻抗發(fā)生較大的變化,冷卻風(fēng)扇故障或液冷循環(huán)系統(tǒng)故障時散熱器的熱阻抗也會發(fā)生大幅的變化,導(dǎo)致溫度估計結(jié)果產(chǎn)生較大誤差[15]。

    現(xiàn)有功率模塊基本都內(nèi)置有負(fù)溫度系數(shù)(Negative Temperature Coefficient, NTC)溫度傳感器,但NTC傳感器并不能直接準(zhǔn)確地反映器件結(jié)溫,通常只用于模塊的內(nèi)部溫度監(jiān)測和過載過溫保護(hù)[20]。目前僅有極少量文獻(xiàn)[21-23]研究了基于NTC溫度傳感器實現(xiàn)器件的動態(tài)結(jié)溫提取。文獻(xiàn)[21-22]采用實驗測試的方法獲取IGBT模塊的熱網(wǎng)絡(luò)參數(shù),但忽略了二極管的熱阻抗及芯片之間的熱耦合效應(yīng),而且缺乏不同邊界條件下的有效性驗證。文獻(xiàn)[23]建立了以NTC熱敏電阻為參考的簡化IGBT熱網(wǎng)絡(luò)模型,得到芯片估計溫度用于逆變器熱優(yōu)化控制,但其研究側(cè)重點為逆變器的熱優(yōu)化控制,并未詳細(xì)討論熱建模過程和散熱條件邊界條件對熱網(wǎng)絡(luò)模型的影響。

    鑒于此,本文提出了一種基于功率模塊內(nèi)置溫度傳感器的功率模塊結(jié)溫在線提取方法。該方法以功率模塊內(nèi)置NTC傳感器作為溫度參考點,建立考慮熱耦合效應(yīng)的內(nèi)置NTC傳感器至功率芯片的新型熱網(wǎng)絡(luò)模型。該熱網(wǎng)絡(luò)模型對模塊外部導(dǎo)熱硅脂老化和散熱器熱阻變化不敏感,結(jié)合準(zhǔn)確的損耗在線計算可以實現(xiàn)器件結(jié)溫的準(zhǔn)確在線提取。本文通過有限元仿真提取熱網(wǎng)絡(luò)模型參數(shù),結(jié)合NTC傳感器和SiC MOSFET的實時損耗計算實現(xiàn)對SiC MOSFET功率模塊內(nèi)部芯片結(jié)溫的在線提取。最后通過實驗和有限元仿真驗證了所提方法的可行性和有效性。

    1 熱網(wǎng)絡(luò)模型建模

    1.1 以內(nèi)置NTC傳感器為參考點的熱網(wǎng)絡(luò)建模

    功率模塊內(nèi)置NTC溫度傳感器并不是安裝在功率芯片表面且測量響應(yīng)速度慢,因此NTC傳感器并不能準(zhǔn)確反映芯片的溫度,本文通過研究NTC溫度與芯片結(jié)溫之間的相關(guān)性,基于傳熱原理構(gòu)建以NTC傳感器為溫度參考點的新型熱網(wǎng)絡(luò)模型,不包含參數(shù)易發(fā)生變化的底部物理層的熱網(wǎng)絡(luò),簡化了傳統(tǒng)熱網(wǎng)絡(luò)模型。

    功率模塊截面熱量流向示意圖如圖1所示。功率芯片產(chǎn)生的熱量經(jīng)過直接覆銅陶瓷基板(Direct Bonded Copper, DBC)層和導(dǎo)熱基板傳遞到散熱器,一部分熱量通過DBC和基板傳遞到NTC,在芯片和NTC之間產(chǎn)生熱耦合。圖2為功率模塊中NTC至單個功率芯片的等效熱模型,其不僅包括傳統(tǒng)熱網(wǎng)絡(luò)模型中芯片至環(huán)境的熱阻抗,同時也包含芯片至NTC傳感器的熱阻抗。

    圖1 功率模塊截面熱量流向示意圖

    圖2 包含NTC功率模塊的等效熱模型

    由圖2可知,芯片產(chǎn)生的熱量傳遞到NTC有兩條路徑,路徑1為芯片—DBC—NTC,路徑2為芯片—DBC—基板—DBC—NTC。由于芯片與NTC之間通過基板交換的熱量很小,并且芯片與NTC的熱容相對于銅基板的熱容可以忽略不計[24],因此可忽略傳遞路徑2,從而得到如圖3所示的簡化等效熱模型。考慮到實際功率模塊中DBC(陶瓷-銅鍵合基板)平面大小及其熱容,DBC層的熱阻抗th.DBC-DBC要比芯片至DBC熱阻抗th.Chip_DBC和DBC至NTC熱阻抗th.DBC_NTC大得多,DBC至環(huán)境的熱阻抗th.DBC_Ambient包含焊料層、基板、導(dǎo)熱硅脂和散熱器的所有熱阻抗,th.DBC_Ambient的值比th.DBC-DBC大得多。因此,NTC傳感器至芯片的熱阻抗th.j-NTC對外部導(dǎo)熱硅脂和散熱器熱阻抗變化不敏感,th.DBC-DBC參數(shù)也只與DBC幾何形狀以及NTC傳感器相對于芯片的距離有關(guān),不易受其他因素的影響。

    圖3 簡化后的等效熱模型

    功率模塊的多物理層熱模型可以使用Foster熱網(wǎng)絡(luò)模型來描述[24],芯片至溫度參考點之間的自熱阻抗曲線可以定義為

    式中,T()和ref()分別為芯片與參考點隨時間變化的溫度;P為芯片產(chǎn)生的功率損耗。

    在Foster熱網(wǎng)絡(luò)模型中,熱阻熱容參數(shù)值都由熱阻抗曲線擬合而來,一般可采用多階指數(shù)曲線擬合獲得兩者之間的熱阻抗參數(shù),即

    1.2 考慮熱耦合效應(yīng)的新型熱網(wǎng)絡(luò)模型

    任何產(chǎn)生功率損耗的芯片在散熱過程中都會對其周圍芯片造成影響,而當(dāng)熱量由功率損耗較小的芯片向較大的散熱器傳遞時,這種熱量橫向傳導(dǎo)、熱流路徑相互疊加的影響也更加明顯,使得周圍所有芯片的溫度升高,這種現(xiàn)象也被稱作交叉熱耦合現(xiàn)象,如圖4所示。文獻(xiàn)[25]表明,當(dāng)模塊處于逆變工作狀態(tài)時,功率芯片熱阻抗比單獨工作時要高約10%。

    圖4 芯片之間交叉熱耦合現(xiàn)象

    因此,為了保證熱網(wǎng)絡(luò)模型的準(zhǔn)確性,還需考慮熱耦合現(xiàn)象。當(dāng)對芯片施加損耗時,溫度參考點和芯片之間的耦合熱阻抗定義為

    式中,P為芯片施加的功率損耗。

    因此,對于多芯片熱源模型,以NTC為參考點時,熱網(wǎng)絡(luò)模型中的耦合熱阻抗網(wǎng)絡(luò)矩陣和自熱阻抗矩陣可分別表示為

    根據(jù)疊加原理可得計及熱耦合影響時的芯片結(jié)溫計算表達(dá)式為

    式中,j、loss、ref分別為功率模塊中各芯片結(jié)溫、各芯片功率損耗和參考點溫度,均為×1向量。將熱阻抗矩陣式(4)、式(5)中的th按實際采樣時間離散化之后,可得到用于結(jié)溫在線提取的離散化熱網(wǎng)絡(luò)。

    綜上所述,當(dāng)參考點為功率模塊內(nèi)部的NTC時,耦合熱阻抗網(wǎng)絡(luò)模型結(jié)構(gòu)如圖5所示,熱耦合的疊加點都設(shè)置在熱阻抗值較大的th.DBC-DBC端。

    圖5 基于NTC的耦合熱網(wǎng)絡(luò)模型

    1.3 SiC功率模塊應(yīng)用于電機驅(qū)動時的損耗分析

    由式(6)可知,在獲得熱阻抗矩陣參數(shù)并完成離散化之后,只要實時采樣NTC傳感器的溫度并根據(jù)電壓、電流和結(jié)溫采樣值計算器件的損耗,就能獲得模塊內(nèi)各芯片的實時結(jié)溫。SiC MOSFET工作狀態(tài)如圖6所示。圖6中,S1為SiC MOSFET,VD1為SiC肖特基勢壘二極管(Schottky-Barrier Diode, SBD)。

    圖6 SiC MOSFET工作狀態(tài)

    SiC MOSFET具有雙向?qū)ㄌ匦裕串?dāng)柵源極電壓大于開啟電壓時,無論漏源極電壓ds是正值還是負(fù)值,溝道均可以導(dǎo)通。因此,在傳統(tǒng)三相逆變系統(tǒng)應(yīng)用中,SiC MOSFET自身會參與反向續(xù)流,導(dǎo)致其損耗特性不同于IGBT。本文SiC MOSFET功率模塊的應(yīng)用場合為電機驅(qū)動逆變器,其損耗分析以A相橋臂為例。

    當(dāng)A相電流由逆變器流向電機負(fù)載時,電流流經(jīng)上橋臂S1,S1處于正向?qū)顟B(tài),反并聯(lián)的二極管VD1處于反向截止?fàn)顟B(tài)。當(dāng)A相電流由電機流向逆變器時,上橋臂將處于反向續(xù)流階段,由于S1為開通狀態(tài),而VD1的導(dǎo)通需要一定的正向電壓,續(xù)流初期的電流將全部從S1流過,該階段為S1的單獨續(xù)流階段;當(dāng)續(xù)流電流流經(jīng)S1所產(chǎn)生的導(dǎo)通壓降達(dá)到VD1的開通閾值電壓TH-VD后,S1與VD1以并聯(lián)的形式共同續(xù)流,流經(jīng)S1和VD1的續(xù)流電流由其各自導(dǎo)通電阻決定。綜上所述,在逆變器工作過程中流經(jīng)S1和VD1的電流可表示為

    圖7 A相橋臂S1和VD1的電流波形

    1.3.1 SiC MOSFET導(dǎo)通損耗計算

    當(dāng)S1導(dǎo)通時,S1的導(dǎo)通壓降DS可表示為

    因此,S1導(dǎo)通期間任意兩時間點的導(dǎo)通損耗可表示為

    當(dāng)開關(guān)頻率較高時,為了減小損耗在線計算的資源占用,可將損耗計算的計算周期設(shè)為個開關(guān)周期,采樣首尾兩個開關(guān)周期的負(fù)載電流和占空比,將采樣首尾周期負(fù)載電流的平均值作為個開關(guān)周期內(nèi)的電流值,首尾兩個開關(guān)周期的占空比平均值作為整個階段的占空比。因此,單個開關(guān)周期內(nèi)的導(dǎo)通損耗計算簡化為

    1.3.2 SiC MOSFET開關(guān)損耗計算

    SiC MOSFET的開關(guān)過程的實際波形是包含尖峰和振蕩過程的不規(guī)則波形,其開關(guān)損耗難以進(jìn)行精確的理論計算。在本文中,只采樣計算S1在個開關(guān)周期中首尾兩個開關(guān)周期的開關(guān)損耗,將個開關(guān)周期內(nèi)的S1平均開關(guān)損耗表示為

    式中,sw()sw(S)分別為首、尾兩個采樣開關(guān)周期內(nèi)器件開通損耗on、關(guān)斷損耗off之和。考慮到on和off受到器件漏極電流D、器件結(jié)溫j、直流側(cè)電壓DC以及柵極驅(qū)動電阻g的影響,其中DC、g一般為固定值,本文根據(jù)器件數(shù)據(jù)手冊中給出的不同條件下的on、off曲線線性插值構(gòu)建on、off關(guān)于D、j的二維查詢表,根據(jù)采樣得到的D、j計算相應(yīng)的開關(guān)損耗。

    1.3.3 SiC SBD 損耗計算

    SiC SBD不存在開關(guān)損耗和反向恢復(fù)損耗,只在其參與續(xù)流期間產(chǎn)生導(dǎo)通損耗,VD1的前向?qū)妷篠BD可近似認(rèn)為是流過電流的線性函數(shù),即

    因此,在兩任意時刻內(nèi)VD1的導(dǎo)通損耗計算表示為

    與S1導(dǎo)通損耗分析同理,得到VD1的導(dǎo)通損耗計算表達(dá)式為

    2 熱網(wǎng)絡(luò)模型參數(shù)提取

    為獲得以NTC傳感器為溫度參考點時各芯片的自熱阻抗和耦合熱阻抗,本文采用Ansys暫態(tài)有限元法(Finite Element Method, FEM)仿真提取熱網(wǎng)絡(luò)模型參數(shù),相應(yīng)流程如圖8所示。

    圖8 FEM仿真提取熱阻抗參數(shù)流程

    2.1 SiC模塊有限元三維建模

    本文的研究對象為CREE CCS020M12CM2三相全橋SiC MOSFET功率模塊,其內(nèi)部包含了6個SiC MOSFET芯片、6個SiC SBD和一個NTC熱敏電阻,CREE CCS020M12CM2三相全橋模塊內(nèi)部結(jié)構(gòu)如圖9所示。

    圖9 CREE CCS020M12CM2三相全橋模塊內(nèi)部結(jié)構(gòu)

    模塊結(jié)構(gòu)及材料參數(shù)如圖10所示。模塊的結(jié)構(gòu)和各層材料參數(shù)如圖10a所示,A相結(jié)構(gòu)及尺寸如圖10b所示。由文獻(xiàn)[25-26]可知,鍵合線、模塊引腳及環(huán)氧樹脂外殼對模塊內(nèi)芯片的散熱影響很小,可將模型簡化,模塊3D有限元仿真模型如圖11所示。

    圖10 模塊結(jié)構(gòu)及材料參數(shù)

    在設(shè)置仿真環(huán)境時,本文將功率模塊基板的底面與側(cè)面設(shè)置為模塊自然對流散熱形式,而芯片表面與DBC上銅層在實際中由于覆蓋了灌封硅凝膠,故認(rèn)為芯片熱量全部從模塊基板流失,其他表面視為絕熱,暫態(tài)熱仿真參數(shù)見表1。

    圖11 模塊3D有限元仿真模型

    表1 暫態(tài)熱仿真參數(shù)

    Tab.1 Parameters for transient thermal simulation

    2.2 熱阻抗參數(shù)提取

    熱阻抗參數(shù)的提取采用熱響應(yīng)擬合法,分別在每個芯片上獨立地施加階躍功率損耗作為激勵,由有限元仿真獲得目標(biāo)芯片的溫度響應(yīng)曲線,再計算暫態(tài)阻抗曲線,熱響應(yīng)曲線獲取流程如圖12所示。

    圖12 熱響應(yīng)曲線獲取流程

    2.2.1 FEM法提取熱阻抗參數(shù)的準(zhǔn)確性驗證

    為驗證所建立的有限元仿真模型的準(zhǔn)確性,本文首先提取功率模塊的結(jié)-殼熱阻抗曲線,并與功率模塊數(shù)據(jù)手冊中的結(jié)-殼熱阻抗曲線進(jìn)行對比。當(dāng)兩者誤差過大時,對有限元模型進(jìn)行修正,修正方法包括調(diào)整材料參數(shù)、微調(diào)幾何尺寸或者提高網(wǎng)格劃分的質(zhì)量。

    圖13所示為FEM仿真得到結(jié)殼熱阻抗曲線與數(shù)據(jù)手冊熱阻抗曲線的對比??梢钥闯觯瑑烧咴谡麄€時間范圍都較為接近,在暫態(tài)值部分FEM高于數(shù)據(jù)手冊,穩(wěn)態(tài)值部分FEM略小于數(shù)據(jù)手冊,但最大誤差都保持在5%以內(nèi),可以認(rèn)為所建立的模型和參數(shù)準(zhǔn)確性符合要求。

    圖13 FEM仿真與數(shù)據(jù)手冊熱阻抗對比

    2.2.2 考慮熱耦合效應(yīng)的熱阻抗參數(shù)提取

    圖14為對A相SiC MOSFET芯片M1施加20W功率損耗后的溫度分布,芯片之間的熱耦合現(xiàn)象會隨著芯片之間距離的增加而快速減弱。當(dāng)以NTC溫度為參考時,通過仿真得到的耦合熱阻抗大小與芯片距離之間的關(guān)系如圖15所示。耦合熱阻抗大小隨芯片間距的增加呈指數(shù)減小,且當(dāng)兩芯片間距大于10mm時,其熱耦合影響可忽略不計。此外,鑒于三相逆變應(yīng)用的對稱性,為簡化熱網(wǎng)絡(luò),本文只對包含NTC熱敏電阻的A相橋臂進(jìn)行研究。

    圖14 模塊溫度分布

    圖15 耦合熱阻抗與芯片間距的關(guān)系

    對A相各芯片單獨施加階躍損耗,按圖12步驟提取的以NTC為參考點的SiC MOSFET芯片M1與SiC SBD芯片D1的自熱阻抗曲線如圖16所示。

    圖16 暫態(tài)自熱阻抗及其擬合曲線

    芯片的自熱阻抗和耦合熱阻抗曲線分別使用四階、一階RC熱網(wǎng)絡(luò)來擬合,得到的熱阻抗參數(shù)見表2和表3,進(jìn)而可構(gòu)建出1.2節(jié)中的考慮熱耦合效應(yīng)的熱網(wǎng)絡(luò)模型。

    2.2.3 邊界條件對熱模型的影響

    功率模塊工作時,其芯片溫度不僅與功率損耗有關(guān),還與散熱條件、環(huán)境溫度等邊界條件有關(guān)。對于散熱條件的變化,可設(shè)置不同的表面?zhèn)鳠嵯禂?shù)(Heat Transfer Coefficient, HTC)值來表征散熱系統(tǒng)的冷卻能力[19],本文在有限元仿真中將功率模塊基板設(shè)置為不同的HTC(5~10 000W/(m2·K))以模擬不同條件的自然空氣冷卻和水冷散熱方式。此外,為討論環(huán)境溫度對熱阻抗的影響,F(xiàn)EM仿真的環(huán)境溫度設(shè)置為25~100℃,提取的芯片M1至NTC的熱阻抗曲線如圖17所示。

    表2 A相芯片熱阻抗參數(shù)

    Tab.2 Parameters of thermal impedance for phase A

    表3 A相芯片耦合熱阻抗參數(shù)

    Tab.3 Coupled thermal impedance of phase A

    由圖17可知,在不同散熱條件下,芯片M1至NTC的暫態(tài)熱阻抗幾乎不變,穩(wěn)態(tài)熱阻在傳熱系數(shù)為10 000W/(m2·K) 時只減小了0.043℃/W,變化率僅為4.7%;在不同環(huán)境溫度下,芯片M1至NTC的熱阻抗會隨著溫度上升而有所增加,但最大差值僅為0.031 9℃/W,變化率也僅為3.5%。上述結(jié)果表明,本文所提出的新型熱網(wǎng)絡(luò)模型受功率模塊的散熱條件和模塊所處的環(huán)境溫度影響很小,熱阻抗參數(shù)穩(wěn)定,更適用于長時間在線估計功率模塊內(nèi)部芯片結(jié)溫。

    3 仿真與實驗驗證

    本文采用紅外測量和有限元仿真兩種方法對所提結(jié)溫估算方法的可行性和準(zhǔn)確性進(jìn)行驗證。結(jié)溫觀測實驗平臺如圖18所示。圖18a為所搭建的實驗平臺,主要參數(shù)見表4。為了能夠使用紅外熱成像儀準(zhǔn)確地測量芯片結(jié)溫,先打開功率模塊封裝,使用化學(xué)溶劑(Adrox 2312)溶解去除功率模塊內(nèi)部的灌封硅凝膠,再對功率模塊均勻噴涂黑色耐高溫絕緣漆(MOTIP 04031)。圖18b為SiC逆變器及處理后的便于紅外測溫的SiC功率模塊。

    本文采用TI TMS320F28379 DSP作為實驗平臺控制芯片,在電機控制程序中集成了基于NTC傳感器采樣溫度的結(jié)溫在線估算功能。為了減少DSP的計算量,控制程序中的損耗計算和結(jié)溫估計模塊的采樣計算周期設(shè)置為10個開關(guān)周期,每次結(jié)溫計算結(jié)果通過CAN轉(zhuǎn)USB通信上傳至計算機,將此結(jié)溫數(shù)據(jù)與紅外熱成像觀測數(shù)據(jù)以及FEM仿真數(shù)據(jù)對比便可驗證所提結(jié)溫估算方法的可行性和準(zhǔn)確性。

    表4 逆變器和負(fù)載電機主要參數(shù)

    Tab.4 Main parameters of the inverter and motor

    3.1 不同工況下的結(jié)溫計算準(zhǔn)確性的驗證

    本文采用紅外熱成像儀對功率模塊A相芯片M1、D1和NTC傳感器溫度進(jìn)行觀測,結(jié)溫紅外定點觀測設(shè)置如圖19所示。

    圖20和圖21為不同工況下的紅外觀測結(jié)果和在線計算結(jié)果。對比可知,在150r/min、10N·m的低頻高轉(zhuǎn)矩工況下,穩(wěn)態(tài)時M1結(jié)溫的計算值和紅外測量值的平均值非常接近,約為45℃,在線計算結(jié)溫為42~49℃,紅外測量的結(jié)溫為43~47℃。A相二極管VD1和NTC傳感器溫度的計算值和測量值也很接近。在750r/min、5N·m工況下,芯片結(jié)溫的波動值較小,在線計算結(jié)溫為35~37℃,紅外測量值的波動只有1℃左右。

    圖19 結(jié)溫紅外定點觀測設(shè)置

    圖20 150r/min,10N·m工況曲線對比

    圖21 750r/min,5N·m工況結(jié)溫曲線對比

    需要說明的是,本實驗使用的紅外熱成像儀為Fotric 326Pro,其最高有效刷新頻率為16Hz,其在電機電流頻率超過16Hz(對應(yīng)電機轉(zhuǎn)速約500r/min)時觀測結(jié)溫波動范圍會有較大失真,但平均結(jié)溫是準(zhǔn)確的。鑒于此,可認(rèn)為本文所提結(jié)溫在線提取方法可行,并具有較高的準(zhǔn)確度。

    3.2 計算結(jié)果與有限元仿真對比

    本文還采用有限元仿真模擬器件的特殊運行工況,將FEM仿真得到的NTC溫度值作為參考輸入到本文構(gòu)建的以NTC為參考點的新型熱網(wǎng)絡(luò)模型中,驗證所提結(jié)溫計算方法在散熱條件突變工況和波動損耗條件下結(jié)溫提取的準(zhǔn)確性。

    3.2.1 恒定大功率損耗、散熱條件突變工況

    為了驗證在散熱條件突變工況下的器件結(jié)溫提取效果,在FEM仿真過程中將功率模塊的基板表面?zhèn)鳠?,系?shù)在120s的仿真時間內(nèi)設(shè)置為三段:0~50s,5 000W/(m2·K);55~80s,50W/(m2·K);80~120s,2 000W/(m2·K)。A相芯片M1和M2的發(fā)熱功率設(shè)置為恒定50W,D1和D2的發(fā)熱功率設(shè)置為恒定5W,模擬系統(tǒng)在重載運行條件下水冷散熱系統(tǒng)突發(fā)故障情況,恒定大功率損耗、散熱條件突變工況的結(jié)溫如圖22所示。

    圖22 恒定大功率損耗、散熱條件突變工況的結(jié)溫

    由圖22可知,在0~50s內(nèi),散熱條件良好情況下(表面?zhèn)鳠嵯禂?shù):5 000W/(m2·K),即強迫水冷),A相芯片溫度在5s內(nèi)達(dá)到穩(wěn)態(tài)值,這一階段結(jié)溫的仿真值和計算值幾乎重合;在55~80s內(nèi),表面?zhèn)鳠嵯禂?shù)突然大幅下降(可等效為水冷系統(tǒng)突然失效),芯片溫度快速上升,本文方法的計算值與有限元仿真結(jié)果非常吻合,誤差在5℃以內(nèi)。

    3.2.2 波動損耗工況

    對于功率器件而言,結(jié)溫波動Δj的快速準(zhǔn)確提取對于器件壽命預(yù)測和優(yōu)化非常重要。為了進(jìn)一步驗證所提結(jié)溫估計方法的有效性,本文針對電機運行于重載和高溫環(huán)境工況進(jìn)行分析,具體的工況條件見表5,結(jié)果如圖23~圖25所示。

    表5 運行工況

    Tab.5 Working conditions

    圖23和圖24為相同輸出轉(zhuǎn)矩、不同轉(zhuǎn)速(即不同逆變器交流輸出頻率)和環(huán)境溫度下的器件結(jié)溫,圖25為模擬惡劣運行條件下(重載且環(huán)境溫度在短時間內(nèi)大幅上升,2~2.5s內(nèi)由50℃上升到100℃)的器件結(jié)溫。由三種工況下的對比結(jié)果可知,所提出的結(jié)溫計算方法在重載及散熱環(huán)境惡劣等條件下均能夠快速準(zhǔn)確地獲取器件的動態(tài)結(jié)溫,誤差在5℃以內(nèi)。

    圖24 工況2的損耗與結(jié)溫對比結(jié)果

    圖25 工況3結(jié)溫對比結(jié)果

    4 結(jié)論

    本文提出了一種基于內(nèi)置溫度傳感器的功率模塊結(jié)溫在線提取方法。首先建立了考慮熱耦合效應(yīng)的內(nèi)置NTC傳感器至功率芯片的熱網(wǎng)絡(luò)模型,并給出了SiC功率模塊應(yīng)用于電機驅(qū)動時的損耗快速計算方法;然后通過有限元仿真提取熱網(wǎng)絡(luò)模型參數(shù),結(jié)合NTC傳感器和SiC MOSFET的實時損耗計算實現(xiàn)了對SiC MOSFET功率模塊內(nèi)部芯片結(jié)溫的在線提取。仿真與實驗對比結(jié)果驗證了所提結(jié)溫估算方法的可行性和有效性。該方法具有以下優(yōu)點:

    1)易于在線集成、結(jié)溫估算準(zhǔn)確。不需要增加額外測量電路,損耗計算和結(jié)溫估計方法易于集成到控制程序中,結(jié)溫平均值和波動值估算準(zhǔn)確。

    2)熱網(wǎng)絡(luò)參數(shù)穩(wěn)定。新型熱網(wǎng)絡(luò)模型參數(shù)對環(huán)境溫度、散熱器工作條件等變化不敏感,非常適用于變換器中功率模塊的長時間結(jié)溫監(jiān)測,并可在此基礎(chǔ)上實施熱優(yōu)化控制。

    在實際中應(yīng)用本文所提方法時,需要對不同的功率模塊進(jìn)行復(fù)雜的三維建模和有限元分析以獲得新型熱網(wǎng)絡(luò)參數(shù),該方法比較繁瑣復(fù)雜。因此,通用、準(zhǔn)確、流程化的熱網(wǎng)絡(luò)參數(shù)提取方法是本文進(jìn)一步研究的重點,更加準(zhǔn)確而實用的損耗建模也有待進(jìn)一步的深入研究。

    [1] Hamada K, Nagao M, Ajioka M, et al. SiC-emerging power device technology for next-generation electri- cally powered environmentally friendly vehicles[J]. IEEE Transactions on Electronic Devices, 2015, 62(2): 278-285.

    [2] 王學(xué)梅. 寬禁帶碳化硅功率器件在電動汽車中的研究與應(yīng)用[J]. 中國電機工程學(xué)報, 2014, 34(3): 371- 379.

    Wang Xuemei. Researches and applications of wide bandgap SiC power devices in electric vehicles[J]. Proceedings of the CSEE, 2014, 34(3): 371-379.

    [3] 趙炫, 蔣棟, 劉自程, 等. SiC功率器件在軌道交通行業(yè)中的應(yīng)用[J]. 機車電傳動, 2020(1): 38-44.

    Zhao Xuan, Jiang Dong, Liu Zicheng, et al. Appli- cation of silicon carbide power devices in rail transit[J]. Electric Drive for Locomotive, 2020(1): 38-44.

    [4] 任磊, 韋徵, 龔春英, 等. 電力電子電路功率器件故障特征參數(shù)提取技術(shù)綜述[J]. 中國電機工程學(xué)報, 2015, 35(12): 3089-3101.

    Ren Lei, Wei Zheng, Gong Chunying, et al. Fault feature extraction techniques for power devices in power electronic converters: a review[J]. Proceedings of the CSEE, 2015, 35(12): 3089-3101.

    [5] Yang Shaoyong, Bryant A, Mawby P, et al. An industry based survey of reliability in power electronic converters[J]. IEEE Transactions on Industry Applications, 2011, 47(3): 1441-1451.

    [6] Wang Huai, Liserre M, Blaabjerg F. Toward reliable power electronics: challenges, design tools, and opportunities[J]. IEEE Industrial Electronics Magazine, 2013, 7(2): 17-26.

    [7] 王學(xué)梅, 張波, 吳海平. 基于失效物理的功率器件疲勞失效機理[J]. 電工技術(shù)學(xué)報, 2019, 34(4): 717-727.

    Wang Xuemei, Zhang Bo, Wu Haiping. A review of fatigue mechanism of power devices based on physics-of-failure[J]. Transactions of China Electro- technical Society, 2019, 34(4): 717-727.

    [8] 周雒維, 王博, 張益, 等. 非平穩(wěn)工況下功率半導(dǎo)體器件結(jié)溫管理技術(shù)綜述[J]. 中國電機工程學(xué)報, 2018, 38(8): 2394-2407, 2549.

    Zhou Luowei, Wang Bo, Zhang Yi, et al. Review on junction temperature management of power semi- conductor devices under power fluctuation condition[J]. Proceedings of the CSEE, 2018, 38(8): 2394-2407, 2549.

    [9] 任磊, 沈茜, 龔春英. 電力電子電路中功率晶體管結(jié)溫在線測量技術(shù)研究現(xiàn)狀[J]. 電工技術(shù)學(xué)報, 2018, 33(8): 1750-1761.

    Ren Lei, Shen Qian, Gong Chunying. Current junction temperature online measurement techniques of power transistors in power electronic converters[J]. Transa- ctions of China Electrotechnical Society, 2018, 33(8): 1750-1761.

    [10] Khatir Z, Dupont L, Ibrahim A. Investigations on junction temperature estimation based on junction voltage measurements[J]. Microelectronics Reliability, 2010, 50(9-11): 1506-1510.

    [11] Xu Zhuxian, Xu Fan, Wang Fei. Junction temperature measurement of IGBTs using short-circuit current as temperature-sensitive electrical parameter for converter prototype evaluation[J]. IEEE Transactions on Indu- strial Electronics, 2015, 62(6): 3419-3429.

    [12] Chen H, Pickert V, Atkinson D J, et al. On-line monitoring of the mosfet device junction temperature by computation of the threshold voltage[C]// Proceedings of the 3rd IET International Conference on Power Electronics, Machines and Drives, The Contarf Castle, Dublin, 2006: 440-444.

    [13] 方化潮, 鄭利兵, 王春雷, 等. IGBT 模塊柵極電壓米勒平臺時延與結(jié)溫的關(guān)系[J]. 電工技術(shù)學(xué)報, 2016, 31(18): 134-141.

    Fang Huachao, Zheng Libing, Wang Chunlei, et al. The relationship between junction temperature and time delay of gate voltage miller plateau of IGBT module[J]. Transactions of China Electrotechnical Society, 2016, 31(18): 134-141.

    [14] 孫鵬飛, 羅皓澤, 董玉斐, 等. 基于關(guān)斷延遲時間的大功率IGBT模塊結(jié)溫提取方法研究[J]. 中國電機工程學(xué)報, 2015, 35(13): 3366-3372.

    Sun Pengfei, Luo Haoze, Dong Yufei, et al. Junction temperature extraction of high power IGBT module based on turn-off delay time[J]. Proceedings of the CSEE, 2015, 35(13): 3366-3372.

    [15] 王莉娜, 鄧潔, 楊軍一, 等. Si和SiC功率器件結(jié)溫提取技術(shù)現(xiàn)狀及展望[J]. 電工技術(shù)學(xué)報, 2019, 34(4): 703-716.

    Wang Lina, Deng Jie, Yang Junyi, et al. Junction temperature extraction methods for Si and SiC power devices-a review and possible alternatives[J]. Transa- ctions of China Electrotechnical Society, 2019, 34(4): 703-716.

    [16] 魏克新, 杜明星. 基于集總參數(shù)法的IGBT模塊溫度預(yù)測模型[J]. 電工技術(shù)學(xué)報, 2011, 26(12): 79-84.

    Wei Kexin, Du Mingxing. Temperature prediction model of IGBT modules based on lumped parameters method[J]. Transactions of China Electrotechnical Society, 2011, 26(12): 79-84.

    [17] 萬萌, 應(yīng)展烽, 張偉. 分立型功率MOSFET結(jié)溫估計的非線性熱網(wǎng)絡(luò)模型和參數(shù)辨識方法[J]. 電工技術(shù)學(xué)報, 2019, 34(12): 2477-2488.

    Wan Meng, Ying Zhanfeng, Zhang Wei. Nonlinear thermal network model and parameter identification method for junction temperature estimation of discrete power MOSFET[J]. Transactions of China Electrotechnical Society, 2019, 34(12): 2477-2488.

    [18] 李輝, 胡姚剛, 劉盛權(quán), 等. 計及焊層疲勞影響的風(fēng)電變流器IGBT模塊熱分析及改進(jìn)熱網(wǎng)絡(luò)模型[J]. 電工技術(shù)學(xué)報, 2017, 32(13): 80-87.

    Li Hui, Hu Yaogang, Liu Shengquan, et al. Thermal analysis and improved thermal network model of IGBT module for wind power converter considering solder fatigue effects[J]. Transactions of China Electrotechnical Society, 2017, 32(13): 80-87.

    [19] Wintrich A, Nicolai U, Werner T, et al. Application manual power semiconductors[M]. Nuremberg: SEMIKRON International GmbH, 2011.

    [20] Wang Ke, Liao Yongjun, Song Gaosheng. Over- temperature protection for IGBT modules[C]// International Exhibition and Conference for Power Electronics, Intelligent Motion, Renewable Energy and Energy Management, Nuremberg, Germany, 2014: 1-7.

    [21] Zhou Yu, Shi Wei, Tang Junsong, et al. Dynamic junction temperature estimation via built-in negative thermal coefficient (NTC) thermistor in high power IGBT modules[C]//IEEE Applied Power Electronics Conference and Exposition Tampa, FL, USA, 2017: 772-775.

    [22] Ma Xin, Zhao Jia, Yang Yong. A new transient thermal impedance model for estimating the dynamic junction temperature of IGBT modules[C]//Inter- national Exhibition and Conference for Power Electronics, Intelligent Motion, Renewable Energy and Energy Management, Nuremberg, Germany, 2018: 1-5.

    [23] Calzo G Lo, Lidozzi A, Solero L, et al. Thermal regulation as control reference in electric drives[C]// 15th International Power Electronics and Motion Control Conference, Novi Sad, Serbia, 2012: 1-7.

    [24] 盛永和, 羅納德P. 柯利諾. 電力電子模塊設(shè)計與制造[M]. 北京: 機械工業(yè)出版社, 2011.

    [25] Morgan A J, Xu Y, Hopkins D C, et al. Decom- position and electro-physical model creation of the CREE 1200V, 50A 3-Ph SiC module[C]//2016 IEEE Applied Power Electronics Conference and Exposition, Long Beach, CA, USA, 2016: 2141-2146.

    [26] Ayadi M, Ammous A, Ounejjar Y, et al. Thermal interaction of semiconductor devices in multi-chip modules[C]//IEEE International Conference on Systems, Man and Cybernetics, Yasmine Hammamet, Tunisia, 2002: 6-16.

    Online Junction Temperature Extraction for SiC Module Based on Built-in Temperature Sensor

    (College of Electrical and Information Engineering Hunan University Changsha 410082 China)

    Silicon carbide (SiC) MOSFET has the advantages of high switching speed and low loss, which will be more and more widely used in high efficiency and high power density applications. In these applications, SiC MOSFETs are facing severe reliability tests, and the online accurate extraction of junction temperature is essential for lifetime prediction and health assessment. This paper proposed an online junction temperature extraction method for power module based on built-in negative temperature coefficient (NTC) temperature sensor. First, the thermal network model between the built-in NTC sensor and the power chips considering the thermal coupling effect was established, and a fast loss calculation method for SiC MOSFET was proposed. Then, the thermal network parameters were extracted by finite element simulation, and the stability of thermal network parameters under different boundary conditions was verified. Simulation and experimental results show that the proposed method can accurately obtain the dynamic junction temperature, and the thermal network model parameters are not affected by changes in boundary conditions such as ambient temperature and heat dissipation conditions. It is suitable for junction temperature monitoring and lifetime prediction under real mission profiles.

    Negative temperature coefficient temperature sensor, SiC MOSFET, thermal network model, chip temperature

    TM46

    10.19595/j.cnki.1000-6753.tces.201442

    國家自然科學(xué)基金項目(51977065)和天水大型電氣傳動系統(tǒng)與裝備技術(shù)國家重點實驗室開放資金項目(SKLLDJ022019003)資助。

    2020-10-30

    2021-02-03

    劉 平 男,1983年生,博士,副教授,研究方向為電動車輛高性能電力電子變換與電機驅(qū)動控制。E-mail: pingliu@hnu.edu.cn(通信作者)

    李海鵬 男,1995年生,碩士研究生,研究方向為寬禁帶半導(dǎo)體器件驅(qū)動與可靠性。E-mail: 2099325157@qq.com

    (編輯 陳 誠)

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    久久99蜜桃精品久久| 午夜视频国产福利| 岛国毛片在线播放| av免费观看日本| 久久ye,这里只有精品| 国模一区二区三区四区视频| 极品教师在线视频| 少妇高潮的动态图| 亚洲精品,欧美精品| 亚洲内射少妇av| 国产免费一区二区三区四区乱码| 又爽又黄a免费视频| 国产一区二区三区综合在线观看 | 大片电影免费在线观看免费| 国产淫片久久久久久久久| 尤物成人国产欧美一区二区三区| 国产精品99久久久久久久久| av国产免费在线观看| 我要看日韩黄色一级片| 老司机影院成人| 精品亚洲成国产av| 日本黄大片高清| 亚洲国产av新网站| 99久久综合免费| 欧美极品一区二区三区四区| 国产欧美亚洲国产| 日韩制服骚丝袜av| 大香蕉97超碰在线| 丝袜喷水一区| 欧美精品一区二区大全| 中文精品一卡2卡3卡4更新| 国模一区二区三区四区视频| av在线app专区| 国产精品.久久久| av在线观看视频网站免费| 亚洲欧洲日产国产| 三级国产精品片| 国产精品国产av在线观看| 两个人的视频大全免费| 日韩av免费高清视频| 高清欧美精品videossex| 国产爽快片一区二区三区| 久久精品国产亚洲av涩爱| 亚洲真实伦在线观看| 日韩强制内射视频| 99热6这里只有精品| 久久婷婷青草| 日日啪夜夜撸| 深爱激情五月婷婷| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 男女边摸边吃奶| 一级a做视频免费观看| 国产老妇伦熟女老妇高清| 亚洲av成人精品一区久久| 精品国产露脸久久av麻豆| 高清日韩中文字幕在线| 久久久久网色| 久热这里只有精品99| 国产女主播在线喷水免费视频网站| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产最新在线播放| 国产免费又黄又爽又色| 久久久亚洲精品成人影院| 爱豆传媒免费全集在线观看| 在线观看国产h片| 亚洲精品久久午夜乱码| 一区二区三区免费毛片| 午夜老司机福利剧场| 在线免费十八禁| 欧美丝袜亚洲另类| 性色av一级| 又爽又黄a免费视频| 日韩欧美一区视频在线观看 | 伦精品一区二区三区| 亚洲经典国产精华液单| 高清黄色对白视频在线免费看 | 亚洲成人中文字幕在线播放| 国产亚洲91精品色在线| 偷拍熟女少妇极品色| 18禁裸乳无遮挡免费网站照片| 水蜜桃什么品种好| 青春草国产在线视频| 丝袜喷水一区| 亚洲国产色片| 久久精品国产a三级三级三级| 少妇人妻精品综合一区二区| 精品人妻熟女av久视频| 久久久久久久久大av| 日韩制服骚丝袜av| 有码 亚洲区| 午夜老司机福利剧场| 日韩亚洲欧美综合| 成年av动漫网址| 免费看不卡的av| 又爽又黄a免费视频| 亚洲精华国产精华液的使用体验| .国产精品久久| 熟女电影av网| 久久精品久久久久久噜噜老黄| 免费av中文字幕在线| 久久午夜福利片| 在线亚洲精品国产二区图片欧美 | 中文资源天堂在线| 熟女电影av网| 日本与韩国留学比较| 丰满人妻一区二区三区视频av| 亚洲天堂av无毛| 啦啦啦在线观看免费高清www| 免费观看性生交大片5| 美女国产视频在线观看| 欧美日韩视频精品一区| 精品国产露脸久久av麻豆| 日韩,欧美,国产一区二区三区| 日韩视频在线欧美| 国产无遮挡羞羞视频在线观看| 高清欧美精品videossex| 欧美成人一区二区免费高清观看| 久久久久人妻精品一区果冻| 97超视频在线观看视频| 国产精品一区二区性色av| 国产精品人妻久久久久久| 免费看日本二区| 97热精品久久久久久| 十分钟在线观看高清视频www | 99热6这里只有精品| 国产成人精品一,二区| 亚洲精品一二三| 亚洲欧美精品自产自拍| 久久99精品国语久久久| 亚洲av成人精品一二三区| 最近中文字幕高清免费大全6| 午夜福利影视在线免费观看| 欧美zozozo另类| av国产久精品久网站免费入址| 精品一区二区三区视频在线| 国产欧美日韩一区二区三区在线 | 国产男女内射视频| 久久97久久精品| 高清欧美精品videossex| 国产精品久久久久久精品电影小说 | 久久人人爽人人爽人人片va| 一区二区三区免费毛片| 狂野欧美激情性bbbbbb| 久久精品国产亚洲网站| 国产成人精品一,二区| .国产精品久久| 欧美xxⅹ黑人| 搡女人真爽免费视频火全软件| 亚洲欧美日韩卡通动漫| 久久久久久久大尺度免费视频| 狂野欧美激情性xxxx在线观看| 久久97久久精品| 免费黄色在线免费观看| 日本欧美视频一区| 伦理电影大哥的女人| 亚洲av中文av极速乱| 一本—道久久a久久精品蜜桃钙片| 黄色怎么调成土黄色| 国产精品一区www在线观看| 尤物成人国产欧美一区二区三区| 熟女av电影| 永久免费av网站大全| www.色视频.com| 18禁在线播放成人免费| 国产老妇伦熟女老妇高清| 99热这里只有精品一区| 国产又色又爽无遮挡免| 亚洲精品亚洲一区二区| 亚洲精品乱久久久久久| 亚洲经典国产精华液单| 国产欧美亚洲国产| av视频免费观看在线观看| 男女无遮挡免费网站观看| 国产免费又黄又爽又色| 日本欧美视频一区| 国产免费又黄又爽又色| 一级毛片aaaaaa免费看小| 日韩电影二区| 人妻 亚洲 视频| 国产午夜精品一二区理论片| 国产免费又黄又爽又色| 99久久精品一区二区三区| 在线观看一区二区三区激情| 免费av不卡在线播放| 国产精品偷伦视频观看了| 国产亚洲欧美精品永久| 观看美女的网站| 色婷婷久久久亚洲欧美| 少妇精品久久久久久久| 我的老师免费观看完整版| 国产精品秋霞免费鲁丝片| 国模一区二区三区四区视频| 亚洲在久久综合| 少妇人妻久久综合中文| 建设人人有责人人尽责人人享有的 | 午夜免费观看性视频| 美女主播在线视频| 狂野欧美白嫩少妇大欣赏| 久久久久国产网址| 搡老乐熟女国产| 伊人久久国产一区二区| 欧美日韩在线观看h| 国产在视频线精品| 岛国毛片在线播放| 久久影院123| 久久久久久久久久久丰满| 国产爱豆传媒在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲欧美成人精品一区二区| 各种免费的搞黄视频| 久久99热6这里只有精品| 亚洲成人av在线免费| 夫妻午夜视频| 免费黄频网站在线观看国产| 美女主播在线视频| 国产精品国产三级国产专区5o| 欧美最新免费一区二区三区| 一级毛片久久久久久久久女| 日韩一区二区视频免费看| 3wmmmm亚洲av在线观看| 亚洲国产精品成人久久小说| 精品国产一区二区三区久久久樱花 | 简卡轻食公司| 少妇高潮的动态图| 一级毛片久久久久久久久女| 亚洲av中文字字幕乱码综合| 亚洲一级一片aⅴ在线观看| 韩国高清视频一区二区三区| 91久久精品电影网| 一级毛片我不卡| 久久亚洲国产成人精品v| 国产精品一及| 美女高潮的动态| 全区人妻精品视频| 男男h啪啪无遮挡| 亚洲第一av免费看| 精品久久久久久久久亚洲| 99久久精品国产国产毛片| 99re6热这里在线精品视频| 日韩中文字幕视频在线看片 | 亚洲av男天堂| 成年av动漫网址| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 久久亚洲国产成人精品v| 97超视频在线观看视频| 国产精品久久久久久久久免| 国产黄色免费在线视频| 麻豆成人午夜福利视频| 亚洲电影在线观看av| 国产精品女同一区二区软件| 赤兔流量卡办理| 国产精品一及| 特大巨黑吊av在线直播| 亚洲真实伦在线观看| 毛片女人毛片| 伊人久久精品亚洲午夜| 我要看日韩黄色一级片| 国产av一区二区精品久久 | 嘟嘟电影网在线观看| 日本黄大片高清| 视频中文字幕在线观看| 免费看不卡的av| 日韩三级伦理在线观看| 香蕉精品网在线| 国产爱豆传媒在线观看| 日韩免费高清中文字幕av| 久久久久久久久久成人| 亚洲精品国产色婷婷电影| 久久久a久久爽久久v久久| 深夜a级毛片| 熟女人妻精品中文字幕| 午夜福利在线观看免费完整高清在| 欧美zozozo另类| 91精品一卡2卡3卡4卡| 亚洲av福利一区| 欧美高清性xxxxhd video| 在线观看免费日韩欧美大片 | 精品久久久久久久末码| 黄色配什么色好看| 一级毛片久久久久久久久女| 18+在线观看网站| 在线观看免费日韩欧美大片 | 日韩av不卡免费在线播放| 亚洲精品国产色婷婷电影| 韩国高清视频一区二区三区| 久久女婷五月综合色啪小说| 国产亚洲欧美精品永久| 亚洲人与动物交配视频| 日韩成人伦理影院| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲精品久久久com| 秋霞伦理黄片| 最黄视频免费看| 好男人视频免费观看在线| 男人狂女人下面高潮的视频| 亚洲三级黄色毛片| 岛国毛片在线播放| 纵有疾风起免费观看全集完整版| 亚洲成人中文字幕在线播放| av又黄又爽大尺度在线免费看| 亚洲av电影在线观看一区二区三区| 日韩欧美精品免费久久| 欧美成人午夜免费资源| 天堂中文最新版在线下载| 免费av中文字幕在线| 欧美成人精品欧美一级黄| 久久国产精品大桥未久av | 欧美日韩视频高清一区二区三区二| 丝瓜视频免费看黄片| 亚洲av欧美aⅴ国产| 在线观看美女被高潮喷水网站| 99久久人妻综合| 纵有疾风起免费观看全集完整版| 青青草视频在线视频观看| 成人综合一区亚洲| 我的女老师完整版在线观看| 免费黄频网站在线观看国产| 97超视频在线观看视频| 成人美女网站在线观看视频| 五月玫瑰六月丁香| 青春草视频在线免费观看| 中国三级夫妇交换| 国产成人精品一,二区| 十分钟在线观看高清视频www | 久久久精品免费免费高清| 免费不卡的大黄色大毛片视频在线观看| 亚洲图色成人| 久久久久久伊人网av| 老熟女久久久| 高清日韩中文字幕在线| 少妇精品久久久久久久| 七月丁香在线播放| 亚洲人成网站高清观看| 日日撸夜夜添| 亚洲国产精品国产精品| 色吧在线观看| 国产日韩欧美在线精品| 精品国产三级普通话版| 纯流量卡能插随身wifi吗| 精品人妻视频免费看| 国产乱来视频区| 91午夜精品亚洲一区二区三区| 简卡轻食公司| 这个男人来自地球电影免费观看 | 在线免费十八禁| 欧美三级亚洲精品| 亚洲精品日本国产第一区| 欧美日韩一区二区视频在线观看视频在线| 日本一二三区视频观看| 亚洲电影在线观看av| 亚洲综合精品二区| 亚洲电影在线观看av| 我要看黄色一级片免费的| 欧美区成人在线视频| 欧美xxxx黑人xx丫x性爽| .国产精品久久| av卡一久久| 亚洲成色77777| 美女内射精品一级片tv| 九九在线视频观看精品| 丰满少妇做爰视频| 国产免费福利视频在线观看| 乱系列少妇在线播放| 91久久精品电影网| 亚洲精品成人av观看孕妇| 91午夜精品亚洲一区二区三区| 亚洲欧美一区二区三区国产| 日日啪夜夜撸| 一级av片app| 国产真实伦视频高清在线观看| 日本黄色片子视频| 国产一区二区在线观看日韩| 国产91av在线免费观看| 少妇人妻久久综合中文| 在线观看一区二区三区激情| 国产亚洲午夜精品一区二区久久| 欧美变态另类bdsm刘玥| 亚洲国产精品一区三区| 国产欧美亚洲国产| av国产久精品久网站免费入址| 夫妻性生交免费视频一级片| videossex国产| 插逼视频在线观看| 色婷婷久久久亚洲欧美| 亚洲精品自拍成人| 免费在线观看成人毛片| 日本与韩国留学比较| 欧美3d第一页| 亚洲欧洲国产日韩| 乱码一卡2卡4卡精品| 中国三级夫妇交换| av专区在线播放| 国产亚洲最大av| 尾随美女入室| 亚洲,一卡二卡三卡| 黄色配什么色好看| 18禁裸乳无遮挡免费网站照片| 免费观看在线日韩| 日韩 亚洲 欧美在线| 午夜福利在线在线| 特大巨黑吊av在线直播| 视频区图区小说| av黄色大香蕉| 自拍欧美九色日韩亚洲蝌蚪91 | 国产爱豆传媒在线观看| 久久久久国产精品人妻一区二区| 在线观看免费日韩欧美大片 | 国产视频内射| 亚洲熟女精品中文字幕| 亚洲av在线观看美女高潮| 美女国产视频在线观看| 青春草视频在线免费观看| 国国产精品蜜臀av免费| 日日摸夜夜添夜夜添av毛片| 久久97久久精品| 国产精品一区二区在线观看99| videossex国产| 日本-黄色视频高清免费观看| 在线免费十八禁| 这个男人来自地球电影免费观看 | 建设人人有责人人尽责人人享有的 | 在线观看美女被高潮喷水网站| 亚洲国产日韩一区二区| 精品久久久久久久末码| 高清黄色对白视频在线免费看 | 中国美白少妇内射xxxbb| 嫩草影院新地址| 成人亚洲精品一区在线观看 | 亚洲色图av天堂| 一级毛片我不卡| 高清黄色对白视频在线免费看 | 中文资源天堂在线| 午夜福利视频精品| 欧美三级亚洲精品| 少妇熟女欧美另类| 久久精品国产自在天天线| av黄色大香蕉| 多毛熟女@视频| 亚洲精品一二三| av在线播放精品| av福利片在线观看| 亚洲精品日本国产第一区| 国产黄色免费在线视频| 中文字幕亚洲精品专区| 日韩大片免费观看网站| 小蜜桃在线观看免费完整版高清| 久久青草综合色| 免费av不卡在线播放| 色婷婷久久久亚洲欧美| 国产精品福利在线免费观看| 亚洲精品一区蜜桃| 少妇的逼好多水| 另类亚洲欧美激情| 国产成人免费无遮挡视频| 国产av精品麻豆| 亚洲国产精品成人久久小说| 精品人妻视频免费看| 久久精品国产鲁丝片午夜精品| av又黄又爽大尺度在线免费看| 精品一品国产午夜福利视频| 深爱激情五月婷婷| 老司机影院毛片| 五月天丁香电影| 亚洲精品,欧美精品| 亚洲av不卡在线观看| 亚洲欧洲国产日韩| 久久精品熟女亚洲av麻豆精品| 中国国产av一级| 国产探花极品一区二区| 一级毛片 在线播放| 国产美女午夜福利| 亚洲精品日韩av片在线观看| 精品人妻一区二区三区麻豆| 一级a做视频免费观看| 91狼人影院| 欧美成人a在线观看| 肉色欧美久久久久久久蜜桃| 99热6这里只有精品| 欧美日韩精品成人综合77777| 观看免费一级毛片| 欧美国产精品一级二级三级 | 汤姆久久久久久久影院中文字幕| 精品99又大又爽又粗少妇毛片| 一本—道久久a久久精品蜜桃钙片| 久热这里只有精品99| 欧美日韩视频高清一区二区三区二| 黄色日韩在线| 狂野欧美激情性bbbbbb| 少妇丰满av| 99热这里只有是精品50| 国产精品女同一区二区软件| 黑丝袜美女国产一区| 国产淫片久久久久久久久| 亚洲美女黄色视频免费看| 国产乱人偷精品视频| 精品久久国产蜜桃| 色综合色国产| 丝袜脚勾引网站| 亚洲成人中文字幕在线播放| 精品久久久噜噜| 超碰av人人做人人爽久久| 国产精品一区二区在线观看99| 下体分泌物呈黄色| 美女主播在线视频| 夜夜爽夜夜爽视频| 日韩成人av中文字幕在线观看| 亚洲内射少妇av| 国产毛片在线视频| 一区二区三区乱码不卡18| 国产黄色免费在线视频| av线在线观看网站| 成人特级av手机在线观看| 伦理电影大哥的女人| 精品一区二区三区视频在线| 男人舔奶头视频| 亚洲不卡免费看| 一级片'在线观看视频| 国产精品99久久久久久久久| 王馨瑶露胸无遮挡在线观看| 91午夜精品亚洲一区二区三区| 纯流量卡能插随身wifi吗| 啦啦啦在线观看免费高清www| 国产精品99久久久久久久久| av不卡在线播放| 欧美精品人与动牲交sv欧美| 在线亚洲精品国产二区图片欧美 | 七月丁香在线播放| 亚洲精华国产精华液的使用体验| 亚洲国产最新在线播放| 免费人妻精品一区二区三区视频| 久久精品国产鲁丝片午夜精品| 欧美另类一区| 黄色日韩在线| 99热6这里只有精品| 国产伦精品一区二区三区视频9| 成年女人在线观看亚洲视频| 国产成人免费无遮挡视频| 欧美精品一区二区大全| 晚上一个人看的免费电影| 国产男人的电影天堂91| 99国产精品免费福利视频| 日本一二三区视频观看| 国产精品无大码| 国产色爽女视频免费观看| 欧美一区二区亚洲| 高清午夜精品一区二区三区| 免费看av在线观看网站| 国产男人的电影天堂91| www.av在线官网国产| 下体分泌物呈黄色| 亚洲av成人精品一二三区| 人妻 亚洲 视频| 王馨瑶露胸无遮挡在线观看| 蜜臀久久99精品久久宅男| 97在线人人人人妻| 久久人人爽av亚洲精品天堂 | 性色av一级| 亚洲精品色激情综合| 国产一区二区三区综合在线观看 | 国产一区二区在线观看日韩| 精品一品国产午夜福利视频| 男男h啪啪无遮挡| 男人舔奶头视频| 中文字幕人妻熟人妻熟丝袜美| 成人影院久久| 人体艺术视频欧美日本| 热re99久久精品国产66热6| 大码成人一级视频| 黄色配什么色好看| 最近的中文字幕免费完整| 99热网站在线观看| 欧美+日韩+精品| 亚洲四区av| 精品久久久久久久久亚洲| a级一级毛片免费在线观看| 91精品伊人久久大香线蕉| 久久久午夜欧美精品| 国产色爽女视频免费观看| 嫩草影院新地址| 国产乱人视频| 又爽又黄a免费视频| 天堂8中文在线网| 在线看a的网站| 中文字幕亚洲精品专区| 黑人猛操日本美女一级片| 99精国产麻豆久久婷婷| 亚洲av综合色区一区| 日韩av免费高清视频| 国产精品一区二区在线观看99| 国产视频内射| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲伊人久久精品综合| 国产精品不卡视频一区二区| 日本欧美国产在线视频| 国内揄拍国产精品人妻在线| 日韩在线高清观看一区二区三区| 高清欧美精品videossex| 精品国产乱码久久久久久小说| 欧美变态另类bdsm刘玥| 一本久久精品| 美女主播在线视频| 国产精品人妻久久久久久| 国产亚洲精品久久久com| 免费播放大片免费观看视频在线观看| av视频免费观看在线观看| 国产在视频线精品| 国产一区二区在线观看日韩| 久久久欧美国产精品| 国内精品宾馆在线| 亚洲av成人精品一二三区| 国产成人一区二区在线| 精品一区二区三区视频在线| 少妇人妻久久综合中文|