《電子技術應用》
您所在的位置:首頁 > 電源技術 > 設計應用 > IGBT模塊開關損耗計算方法綜述
IGBT模塊開關損耗計算方法綜述
電子技術應用2016年第1期
李志剛,梅 霜,王少杰,姚 芳
河北工業大學 電氣工程學院,天津300130
摘要: IGBT模塊開關損耗計算方法有基于物理方法的損耗計算法和基于數學方法的損耗計算法。 基于物理方法的損耗計算是采用軟件仿真的辦法建立相應的物理模型得到開關動態波形,計算損耗;基于數學方法的損耗計算是采用各種數學的方法建模,計算損耗。對各種開關損耗計算法進行討論,并給出其相應的優缺點。
中圖分類號: TN301
文獻標識碼: A
DOI:10.16157/j.issn.0258-7998.2016.01.001
中文引用格式: 李志剛,梅霜,王少杰,等. IGBT模塊開關損耗計算方法綜述[J].電子技術應用,2016,42(1):10-14,18.
英文引用格式: Li Zhigang,Mei Shuang,Wang Shaojie,et al. The review of IGBT module switching loss calculation method[J].Application of Electronic Technique,2016,42(1):10-14,18.
The review of IGBT module switching loss calculation method
Li Zhigang,Mei Shuang,Wang Shaojie,Yao Fang
School of Electrical Engineering,Hebei University of Technology,Tianjin 300130,China
Abstract: Switching loss calculation methods of IGBT module are divided into two types:one is based on the physics,the other is based on mathematics. In this paper,switching loss calculation method based on the physics establish the corresponding physical model by the simulation software to obtain the switch dynamic waveform and calculate the loss. Switching loss calculation method based on the mathematics use all kinds of methods of mathematical modeling for calculation of loss. Various switching loss calculation method and its corresponding advantages and disadvantages are discussed.
Key words : IGBT module;switching loss;loss calculation;loss model

0 引言

    絕緣柵型晶體管(Insulated Gate Bipolar Transistor,IGBT)是由MOSFET和功率雙極型晶體管復合而成的一種器件。IGBT既具有MOSFET的高速開關及電壓驅動特性,又有功率雙極型晶體管(BJT)的低飽和電壓特性及易實現較大電流的能力,在工業、能源、交通等場合越來越不可取代[1]。雖然在電力電子電路中,IGBT主要工作在開關狀態,但是IGBT仍然是功耗較大的電子器件,隨著開關頻率的增高,開關損耗會隨之增大而成為器件功耗的主要因素[2]

    IGBT模塊的性能與其開關特性密切相關,器件的開關特性直接決定其開關損耗,開關損耗制約著器件的工作效率的提高。而且功率器件IGBT的開關損耗可能會產生很高的熱量,引起過大的溫升,對器件的可靠性影響很大。因此,IGBT模塊的開關損耗問題一直是各國學者研究的熱點,其中如何準確估算IGBT模塊的開關損耗是研究的重點內容之一。損耗計算對系統設計、壽命預測、選擇合適的散熱系統、提高系統的可靠性很重要[3-6]

    目前,國內外有關IGBT開關損耗研究的文獻很多,這些文獻歸納總結可以看出:IGBT模塊開關損耗的計算方法主要分為基于物理方法的損耗計算法和基于數學方法的損耗計算法兩種。本文對近年來各國學者們對IGBT模塊開關損耗的計算方法進行討論,并給出其相應的應用范圍。

1 開關損耗定義

    IGBT模塊的開關瞬態電壓、電流波形及開關損耗如圖1所示。

zs-t1.gif

    IGBT模塊的開關損耗由IGBT的開關特性決定,與其集-射極間電壓Vce及集電極電流Ic有關。損耗計算公式如下所示[7]

    zs-gs1.gif

    zs-gs2.gif

其中:Pon為開通損耗,Poff為關斷損耗,ton為開通時間,toff為關斷時間,vce為集-射極間電壓,ic為集電極電流。

2 基于物理方法的開關損耗計算

    基于物理方法的IGBT模塊開關損耗計算方法是采用仿真軟件,使用電源、電容、電阻、電感及晶閘管等一些相對簡單的元件搭建器件物理模型來仿真IGBT模塊的動態特性,得到IGBT模塊的開關瞬態電流、電壓的波形,從而計算開關損耗。Hefner、Kraus、Sheng等分別建立物理模型,采用該方法詳細描述了損耗計算過程,并不斷改進。該方法的開關損耗計算流程圖如圖2所示。

zs-t2.gif

    這種計算方法的準確程度主要取決于IGBT損耗模型的精確度和模型參數的準確度,采用的物理模型越接近IGBT實際器件,模型參數越接近實際大小,仿真計算IGBT損耗值才能越接近實際損耗值。

    目前,常用的用于IGBT模塊動態仿真的軟件主要有三種:Saber、Pspice、Matlab。

    其中Saber提供的是包括Hefner模型在內的5個通用模型和各種精確的具體型號器件的專用模型[8-9],Pspice中提供詳細的器件仿真模型[10-11],這些模型基于IGBT物理結構包含了其重要的物理特征,可以描述IGBT在各種外部電路條件下的穩態和動態特性,具有很好的動態精確性。但是這些模型中多種參數值的確定對仿真影響較大,而且對一般的使用器件的用戶來說,模型參數值的確定是比較復雜和困難的。而Matlab僅提供理想的器件模型,模型通用性較好,但其可設置IGBT參數較少,在開關動態過程的精確描述方面存在很大欠缺[9,11]

    Hefner首先提出了“非準靜態近似理論”[13-17],傳統的“準靜態近似理論”忽略了電子和空穴電流相互耦合會造成集電極電流的變化,以及由于基區的快速變化造成有效輸出電容的數量級的改變,而這些對IGBT模塊物理模型相當重要,進而會影響到IGBT模塊開關損耗值,因此準靜態近似理論不再適用于對IGBT等具有電導調制效應的功率器件進行瞬態分析。Hefner模型是第一個完整的一維分析、電荷控制模型,被廣泛應用于電路仿真中,模型中包含了IGBT模塊的重要物理特征,可以描述IGBT模塊在各種外電路條件下的穩態特性,采用非線性電容,也能很好地表現器件的動態特性,與實驗結果相比,具有較準確的結果。一般可在Saber和Pspice軟件中實現。

    Kraus采用一個多項式來逼近動態過剩載流子濃度分布[18-20],更加貼近IGBT模塊內部載流子的運動,其影響IGBT模塊動態開關波形。Kraus模型是將IGBT模塊看成MOSFET和BJT組成,僅適用于NPT-IGBT的建模。MOSFET部分采用電阻和電容來表達其特征,BJT部分采用一個二極管和三個電流源來描述。Kraus模型簡單易懂,但內部理論較為復雜,很難在一般仿真軟件中實現,主要用于Saber軟件的仿真中。

    Sheng模型主要用于對D-IGBT的建模[22],采用二維載流子分布的方法描述IGBT靜態特性,同時綜合考慮了其動態特性和溫度對器件的影響,其主要用于Pspice仿真軟件中。

    基于物理方法的開關損耗計算優點是物理模型精度高(精度取決于模型的不同簡化程度和參數的準確度),針對器件的具體結構和工藝建立的仿真模型精確表示了器件的動態和靜態特性;缺點是物理模型的構建非常困難,仿真速度慢,而且模型參數的獲得也比較困難。

3 基于數學方法的開關損耗計算

3.1 基于數據手冊的開關損耗計算

    文獻[23]提出一種基于數據手冊估算IGBT模塊開關損耗的計算方法,該方法最早是典型IGBT模型生產商在其數據手冊中給出的[24],數據手冊給出典型產品的開關能量曲線如圖3所示,從中可以看出關斷能量為集電極電流的線性函數,開通能量為集電極電流的二次函數。

zs-t3.gif

    利用線性插值的方法可以得到關斷損耗:

zs-gs3-7.gif

    式中:user代表實際測量的數據,data代表數據手冊中的數據。

    利用二次插值的方法可以得到開通損耗:

zs-gs8-12.gif

    展開可以求出Aon,Bon,Con

    基于數據手冊的開關損耗計算方法的優點是直接使用數據表中數據,簡單方便,缺點是計算不精確。供應商提供的數據是基于實驗室條件下的數據,而在實際工況下運行條件與供應商的實驗工況不一致,其開關能量曲線必然有一定差異,計算損耗也會不同。

3.2 基于數學模型的開關損耗計算

    基于數學模型的IGBT模塊開關損耗計算方法就是在對實際運行條件下的大量的IGBT開關損耗數據進行一定歸納分析的基礎上,依據損耗隨各參數的變化趨勢,建立損耗與各個影響因子之間的數學關系,進而建立器件開關損耗數學模型。

    用于開關損耗建模的數學模型主要有多項式模型[25,28],冪函數模型,多項式和冪函數組合模型,多維數據庫模型[30-31],人工智能模型(如神經網絡模型,迷糊邏輯模型)等。

    冪函數模型是將開關損耗表示成電流的冪函數的形式[26],當考慮到不同的母線電壓和結溫對開關損耗的影響時,其表達式為:

    zs-gs13.gif

    式中:PSW.X為開通損耗,Vdc為實際母線電壓,Vb為母線電壓基值,Tj為結溫,Tb為基值結溫,ASW.X、BSW.X、CSW.X、DSW.X為常數,可經過大量的擬合得到。該公式同時適用于開通損耗和關斷損耗。

    冪函數模型的系數較少,擬合速度快,但擬合效果較差,在此基礎上又提出了多項式與冪函數模型相結合的建模方法[27],多項式可以提高擬合精度,但系數的無限增多又會大大降低擬合速度,因此在建立多項式與冪函數模型時要注意二者的平衡。

    神經網絡具有很好的函數逼近能力,通過對樣本的訓練,能夠很好地反映對象輸入輸出之間的映射關系,同時神經網絡又具有高度的魯棒性,因而能取得較好的建模效果[33]。建立神經網絡模型采用的是誤差反傳網絡(簡稱BP網絡),確定開關損耗模型采用的是三層網絡,結構如圖4所示,輸入層由幾個相關的變量組成如電壓、電流、結溫、門極電壓、門極電阻等,輸出層為開通損耗或者關斷損耗。最后選取的具體訓練參數越多,訓練次數越多,誤差越小。

zs-t4.gif

    神經網絡模型對大量非結構性、非精確性規律具有自適應功能,具有的信息記憶、知識推理和優化計算的特點,大大方便了開關損耗建模,其不管器件的物理機理,只對其外特性進行輸入輸出的智能學習,從而實現對IGBT輸入輸出系統的模擬,得到實況輸出的準確預測;但神經網絡的學習過程通常較慢,對突發事件的適應性差。

    采用數學模型的開關損耗計算方法的優點是規避器件的物理機理,構造模型盡可能的簡單和提高仿真速度;缺點是基于大量的實驗數據擬合,準備工作比較復雜。

3.3 基于波形擬合的開關損耗計算

    基于波形擬合的IGBT開關損耗計算方法的基本思路是:分析IGBT模塊開通、關斷時物理機理及得到的電壓、電流的波形,總結其開關暫態過程的主要特征,建立相應的簡單函數表達開通、關斷的電壓、電流波形,使該模型下獲得的開通、關斷時的電壓、電流的波形無限逼近開關暫態實測波形,然后通過電壓、電流方程的積分獲得開關損耗。

    如圖5、圖6所示分別為IGBT開通、關斷時的電壓、電流波形,采用將開關過程分段,每段分別使用簡單函數逼近實測波形。該方法可以將IGBT并聯二極管的損耗考慮在內,現代IGBT采用的快速恢復二極管的開通損耗與關斷損耗相比可以忽略不計(小于1%)[38]

zs-t5.gif

zs-t6.gif

    但是該方法的暫態波形表達式不夠完整,而且其中可能將某些變量認為恒定,因此獲得的曲線只能逼近實際波形曲線,想要真正與試驗獲得的數據一致是不可能的。

    文獻[42-43]完全忽略了開關暫態電壓和電流的拖尾過程;文獻[44]考慮了拖尾過程,但在擬合過程中,將電流上升時問、續流二極管的反向恢復電流、反向恢復時間等重要參數都視為恒定值,電壓源換流器實際工作過程中,工況是經常變化的,IGBT的開關電流也是變化的,且上述參數都是隨開關電流的變化而變化的;文獻[45]未考慮線路雜散參數影響,同時提出了大量需要由實測波形數計算得出的待定系數,將開關電流與擬合參數的關系復雜化,降低了在高壓大功率場合下的實用性;文獻[46]提出的模型擬合電壓波形與實驗電壓波形基本一致,尤其是關斷電壓上升階段和電壓過沖階段的準確擬合,電壓過沖衰減過程的衰減趨勢也能較好地擬合,但模型未計入電路雜散電容及對地電容,擬合波形沒有實驗波形中的振蕩衰減現象。

    采用曲線擬合的IGBT開關損耗計算方法優點是相對準確,且簡單實用,但是現有開關損耗模型研究存在明顯的近似處理,對開關暫態波形的描述不夠完整。雖然這些擬合波形無法與實驗波形達到絕對一致,但擬合方法也在不斷改進,近似度不斷提高,因此基于曲線擬合的IGBT開關損耗計算是絕對可行的。

4 結語與展望

    通過對IGBT模塊開關損耗的各種計算方法進行分類和分析,可大致了解其各自的優缺點如表1所示。基于物理方法計算開關損耗計算方法需要建立IGBT的等效物理模型,參數提取也比較復雜,但目前許多仿真軟件已經建立了一些器件的物理模型,大大方便了人們的應用。基于數學方法的開關損耗計算方法大都是以數據手冊和實驗數據為基礎,規避器件的復雜的物理結構分析,由于廠商提供器件資料的測試環境的確定性導致基于數據手冊的開關損耗計算方法無法獲得工況下器件的精確開關損耗值,但計算方法簡單易懂,具有通用性;基于數學模型和波形擬合的開關損耗計算方法都以實測的開關波形和開關損耗值為基礎,準備工作復雜,但模型建立后獲取開關損耗值比較簡單,計算速度快,而且精度較高。

zs-b1.gif

    現有的IGBT模塊開關損耗計算方法在其精確度和通用性上存在一定的局限性,這是由于器件的多樣性、工況的復雜性等造成的。IGBT模塊是開關器件,其開關頻率影響IGBT模塊的導通和關斷過程,因此接下來可以在開關損耗計算方法中引入開關頻率對其影響,使損耗模型更加精確。其次在基于神經網絡損耗模型成功的基礎上,可以采用其他人工智能模型對IGBT模塊開關損耗建模,獲得更精確的開關損耗值。

參考文獻

[1] 周志敏,周紀海,紀愛華.IGBT和IPM及其應用電路[M].北京:人民郵電出版社,2006.

[2] 周雒維,吳軍科,杜雄,等.功率變流器的可靠性研究現狀及展望[J].電源學報,2013(1):1-15.

[3] 熊妍,沈燕群,江劍,等.IGBT損耗計算和損耗模型研究[J].電源技術應用,2006(5):55-60.

[4] FEIX G,DIECKERHOFF S,ALLMELING J.Simple methods of calculate IGBT and diode conduction and switching losses[C].13th European Power Electronics  Conference,Barcelona:EPE,2009:1-8.

[5] MASWOOD A I.A switching loss study in SPWM IGBT inverter[C].Proceedings of 2nd IEEE International Conference on Power and Energy.Malaysia:IEEE 2008:609-612.

[6] Wu Rui,Wen Jialiang.A power loss calculation method of IGBT three-phase SPWM converter[C].Proceeding of 2012 InternationalConference on Intelligent Systems Design and Engineering Applications.2012:1180-1183.

[7] 王博宇.混沌SPWM控制AC-DC變換器IGBT的損耗計算方法及溫升研究[D].北京:北京交通大學,2014.

[8] RAJAPAKSE A D,GOLE A M,WILSON P L.Electromagnetic transient simulation models for accurate representation of switching  losses  and  thermal  performance  in  power electronic systems[J].IEEE Trans on Power Delivery 2005,20(1):319-327.

[9] 夏加寬,韓文澤,顧立天.Saber和Matlab在IGBT動態特性仿真中的應用研究[J].電子世界,2013(23):68-69.

[10] 穆懷宇.基于PSpice的IGBT建模與損耗仿真分析[D].哈爾濱:哈爾濱理工大學,2013.

[11] 蘇麟,梅軍,鄭建勇,等.PSPICE和Matlab在IGBT動態仿真中應用[J].電力自動化設備,2004(6):58-61.

[12] MIYAKE M,UENO M,FELDMANN V,et al.Modeling of SiC IGBT turn-off behavior valid for over 5-kV circuit simulation[J].IEEE Trans.on Power Delivery 2012,60(2):622-629.

[13] HEFNER A H.Analytical modeling  of devicr-circuit interactions for the power insulated atr bipolar transistor(IGBT)[J].IEEE Transactions on lrulractryApplications,1990,26:995-1005.

[14] 張日文,李肇基.IGBT的傳輸特性分析[J].微電子學,1996(5):8-12.

[15] HEFNER A R.An investigation of the drive circuit requirements for the power insulated gate bipolar transistor(IGBT)[J].IEEE Transactions  on  Power Electronics,1991(6):208-219.

[16] HEFNER A R.A dynamic electro-thermal model for the IGBTS[J].IEEE Transactions on Industry Applicationa,1999,30(2):394-405.

[17] 張薇琳,張波,丘東元.電力電子開關器件仿真模型分析和比較[J].電氣應用,2007(9):64-67.

[18] KRAUS R,HOFFINANN K.Analytical model of IGBTs with int symposium power semiconductor devices ICs.IEEE IGBT.IEEE Transactions on Power Low Emitter Efficiency.Proceedings,30-34.

[19] KRAUS R,TURKES P,SIGG J.Physics based models of power semiconductor devices for the circuit simulator spice[C].Procedings lEEE PESC'98,1998:1726-1731.

[20] KRAUS R,Mattaush H J.Status and trends of power semiconductor device models for circuit simulation[J].IEEE Transactions on Power Electronics,1998,13:452-465.

[21] SHENG K,WILLIAMS B W,FINNEY S J.A review of IGBT molels[J].Transactions on Power Electronics,2000,15:1250-1266.

[22] SHEN Y,XIONG Y,JIANG J,et al.Switching loss analysis and modeling of power semiconductor devices base on an automatic measurement system[C].In Proceedings of the 2006 IEEE International Symposium on Industrial Electronics,2006,853858.

[23] BIERHOFF  M  H,FUCHS F W.Semiconductor losses in voltage source and current source IGBT converters based on  analytical  derivation[C].Proceedings  of  35th  Annual IEEE Power Electronics Specialists Conference.Aachen,Germany:IEEE,2004:2836-2842.

[24] Calculation of Major IGBT Operating Parameters[R].Inflneon Application Notes,http://www.infineon.com.

[25] 王瑞.IGBT功耗計算及其優化的研究[J].現代電子技術,2014(17):122-124.

[26] BLAABJERG F,PEDERSEN J,JAEGER U.Evaluation of modern IGBT-modules for Hard-switched AC/DC/AC Converters[J].IAS′95.1995,12:997-1005.

[27] Shen Yanqun,Xiong Yan,Jiang Jian,et al.Switching loss analysis and modeling of power semiconductor devices base on an automatic measurement system[J].Industrial Electronics,2006(2):853-858.

[28] 胡建輝,李錦庚,鄒繼斌,等.變頻器中的IGBT模塊損耗計算及散熱系統設計[J].電工技術學報,2009(3):159-163.

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

[30] 何湘寧,楊兵建,吳巖松,等.大功率三電平逆變器開關動態特性在線測試方法研究[J].中國電機工程學報,2012(24):23-29,6.

[31] 何湘寧,吳巖松,羅皓澤,等.基于IGBT離線測試平臺的功率逆變器損耗準在線建模方法[J].電工技術學報,2014(6):1-6.

[32] 陳娜.中高壓功率IGBT模塊開關特性測試及建模[D].杭州:浙江大學,2012.

[33] 陳娜,何湘寧,鄧焰,等.IGBT開關特性離線測試系統[J].中國電機工程學報,2010,30(12):50-55.

[34] 劉璐,朱國榮,陳皓,等.三相四線制VIENNA整流器損耗分析與熱設計[J].電工技術學報,2014,S1:282-290.

35-47略.

此內容為AET網站原創,未經授權禁止轉載。
亚洲一区二区欧美_亚洲丝袜一区_99re亚洲国产精品_日韩亚洲一区二区
国产欧美精品日韩| 国产一区二区三区四区在线观看| 欧美在线观看网址综合| 一区二区三区四区精品| 亚洲精品美女| 亚洲国产高清视频| 亚洲国产高清在线| 亚洲高清视频在线观看| 国产精品精品视频| 国产欧美一区二区精品性| 久久午夜视频| 久久视频在线视频| 久久综合九色九九| 美女露胸一区二区三区| 开心色5月久久精品| 久久久久综合一区二区三区| 久久久久久夜| 久久五月激情| 另类激情亚洲| 欧美国产亚洲精品久久久8v| 欧美精品www在线观看| 欧美激情在线狂野欧美精品| 欧美激情网站在线观看| 欧美精品性视频| 欧美日韩国产成人在线观看| 欧美日韩美女在线观看| 欧美性猛交xxxx乱大交退制版| 国产精品久久久久aaaa樱花| 国产精品久久久久久五月尺| 国产精品青草久久| 国产人成一区二区三区影院| 国产综合香蕉五月婷在线| 伊人影院久久| 亚洲三级影院| 亚洲性xxxx| 欧美一级专区| 亚洲欧洲在线观看| 一区二区三区成人精品| 亚洲一区影院| 久久精品二区三区| 久久久中精品2020中文| 欧美电影资源| 国产精品成人播放| 国产一区亚洲一区| 亚洲国产精品久久久久秋霞影院| 99国产麻豆精品| 性做久久久久久免费观看欧美| 久久精品久久99精品久久| 亚洲精品美女91| 亚洲欧美日韩国产精品| 久久精品一区二区三区四区 | 在线观看免费视频综合| 亚洲精品一区在线| 性欧美在线看片a免费观看| 亚洲激情电影在线| 亚洲自拍16p| 裸体歌舞表演一区二区| 欧美日韩在线精品一区二区三区| 国产女同一区二区| 亚洲国产日韩综合一区| 中文精品视频| 91久久久久久国产精品| 亚洲综合色激情五月| 久久一区二区三区四区| 国产精品a久久久久| 一区二区三区亚洲| 亚洲一区二区精品在线| 亚洲精品免费看| 欧美在线视频一区二区三区| 欧美精品免费看| 国模私拍视频一区| 一区二区三区波多野结衣在线观看| 欧美在线三区| 亚洲一区二区三区乱码aⅴ| 久久综合久久久| 国产精品免费一区二区三区观看| 在线成人性视频| 午夜久久久久久| 一本色道88久久加勒比精品 | 狂野欧美激情性xxxx| 国产精品久久久久久久app| 亚洲第一精品夜夜躁人人躁| 亚洲欧美日韩精品综合在线观看| 亚洲美女诱惑| 久久亚洲一区二区三区四区| 国产精品国产三级国产普通话蜜臀 | 亚洲视频在线看| 欧美福利视频一区| 国内外成人在线视频| 亚洲一区精品视频| 99精品欧美一区二区蜜桃免费| 久久久久久久91| 国产精品美女一区二区在线观看 | 亚洲久久视频| 亚洲日产国产精品| 久久最新视频| 国产色视频一区| 亚洲无玛一区| 亚洲视频在线观看三级| 欧美久久99| 亚洲国产精品久久久久秋霞不卡 | 91久久久精品| 亚洲国产精品ⅴa在线观看| 欧美资源在线观看| 国产精品视频免费在线观看| 日韩一区二区精品| 一本一本久久| 欧美精选午夜久久久乱码6080| 亚洲福利视频网站| 亚洲第一综合天堂另类专| 久久成人18免费观看| 国产嫩草一区二区三区在线观看 | 久久精品亚洲乱码伦伦中文 | 欧美精品久久天天躁| 亚洲高清免费在线| 亚洲三级网站| 欧美国产日韩在线观看| 亚洲第一视频| 日韩亚洲在线| 欧美日本二区| 一本色道久久综合亚洲精品不卡 | 欧美高清视频一区二区| 在线视频国内自拍亚洲视频| 亚洲高清电影| 美女精品网站| 亚洲国产日韩在线一区模特| 亚洲看片一区| 欧美日本一道本| 99热精品在线| 小黄鸭精品密入口导航| 国产精品一区久久久| 午夜精品一区二区三区在线视| 欧美一区三区三区高中清蜜桃| 国产免费成人| 亚洲盗摄视频| 欧美成人精品影院| 亚洲精品自在久久| 亚洲欧美综合精品久久成人| 国产精品综合av一区二区国产馆| 欧美一区1区三区3区公司| 久热精品视频在线观看| 亚洲国产成人tv| 日韩视频免费观看高清在线视频| 欧美日韩精品一区| 亚洲图片欧美日产| 久久成人免费| 伊人精品成人久久综合软件| 亚洲美女啪啪| 欧美视频精品在线观看| 亚洲一二三四久久| 久久蜜桃资源一区二区老牛| 亚洲丰满在线| 亚洲一区二区三区欧美| 国产精品欧美激情| 亚洲国产成人av好男人在线观看| 欧美成人一区在线| 日韩一级在线观看| 久久成人国产精品| 亚洲国产欧美在线人成| 亚洲一级高清| 国产一区日韩二区欧美三区| 亚洲欧洲精品成人久久奇米网| 欧美日韩国产专区| 亚洲欧美不卡| 欧美激情国产高清| 亚洲一区在线播放| 蜜臀99久久精品久久久久久软件 | 久久精品一级爱片| 亚洲人妖在线| 欧美中文字幕在线播放| 亚洲国产成人精品女人久久久 | 亚洲免费人成在线视频观看| 久久这里有精品15一区二区三区 | 一本色道久久综合狠狠躁的推荐| 欧美在线在线| 亚洲国产另类 国产精品国产免费| 亚洲在线一区| 亚洲第一精品在线| 午夜激情综合网| 国外成人在线| 亚洲综合色噜噜狠狠| 亚洲一区影院| 极品尤物av久久免费看| 日韩午夜免费| 国产欧美一区二区精品性| 亚洲精品一区二| 国产日韩精品在线观看| 亚洲第一天堂av| 国产精品国产三级国产a| 亚洲国产综合在线| 国产美女精品视频免费观看| 亚洲九九精品| 国产亚洲欧美中文| 亚洲一级黄色av| 亚洲国产精品小视频| 久久精品中文字幕一区| 一本色道久久综合亚洲精品不 | 国产综合18久久久久久| 亚洲婷婷综合色高清在线| 依依成人综合视频|