999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于拓展型雙曲正切函數的全橋LLC變換器大信號建模方法

2023-01-30 06:11:56孫城皓孫秋野
電工技術學報 2023年1期
關鍵詞:信號方法模型

孫城皓 孫秋野 王 睿 胡 杰

基于拓展型雙曲正切函數的全橋LLC變換器大信號建模方法

孫城皓 孫秋野 王 睿 胡 杰

(東北大學信息科學與工程學院 沈陽 110819)

連續模型的準確性對于電力電子變換器的特性分析與實時仿真十分重要。由于獨特的諧振腔結構,LLC變換器在不同工作頻率區間分別以連續電流模式和斷續電流模式運行,故而傳統方法很難以連續模型實現其多種運行模式的統一描述。為此,該文提出一種基于拓展型雙曲正切函數的全橋LLC變換器大信號建模方法。首先,分析各工作頻率區間全橋LLC變換器特性,建立一種大信號非連續模型,該模型能夠實現LLC變換器多種運行模式的統一描述。然后,構建含陡度因子/脈沖系數的拓展型雙曲正切函數,利用其將上述大信號模型連續化,從而建立了全橋LLC變換器的大信號連續模型。基于此模型,連續系統的方法可以被直接應用以實現對變換器特性分析、實時仿真和控制器設計。同時,相比于現有LLC變換器大信號連續模型,該模型具有較低的階數和更高的精確度。此外,該模型在時域中能夠提供次級側功率器件高準確性的開關信息,便于作為同步整流設計的參考。最后,仿真和實驗結果驗證了該模型的準確性和有效性。

全橋LLC變換器 連續系統 大信號模型 拓展型雙曲正切函數

0 引言

作為一種諧振型DC-DC變換器,LLC變換器因其高功率密度、高效率、低電磁干擾的特性,被廣泛應用于數據中心、航空航天、光伏、電動汽車中[1]。由于獨特的諧振腔結構,LLC變換器在不同頻率區間分別以連續電流模式和斷續電流模式運行,工作模態較為復雜[2],這使得建立準確的LLC變換器大信號連續模型較為困難。然而,連續模型的準確性對電力電子變換器的動態特性分析和實時仿真起關鍵作用[3]。因此,建立能夠準確描述LLC變換器全部運行模式的大信號連續模型具有重要意義。

傳統的LLC變換器建模通常采用基波近似法(Fundamental Harmonic Approximation, FHA)。該方法基于頻域分析,僅考慮LLC諧振變換器初、次級側方波電壓的基波分量以建立變換器穩態等效電路模型,進而繪制變換器的增益曲線并推導相關參數的設計公式[4-5]。但是,FHA方法建模的前提是將LLC變換器的輸出視為連續電流模式,這會導致變換器在工作頻率小于諧振頻率時,即變換器運行于斷續電流模式時模型準確性下降[6]。進而基于改進FHA方法的LLC變換器等效電路模型被提出,以描述LLC變換器的斷續電流模式[7-8],然而,此類基于頻域的建模方法難以描述變換器的動態特性。

為得到準確的動態模型,近年來,一些基于時域分析的建模方法被陸續應用于LLC變換器,可分為以下幾種:離散時間映射法[9-10]、廣義狀態空間平均法[11-12]和擴展描述函數法[13-17]。其中,離散時間映射法從變換器系統的各段狀態方程出發,將變換器中電力電子器件的開關周期作為采樣間隔,通過使變換器的狀態不斷地從一個采樣時間映射到下一個采樣時間來建立等效模型,由于建模過程中不存在近似,該方法的建模精度較高[9-10]。然而,這種模型的建立需要大量的數值運算,建模步驟繁瑣,計算量大,導致難以被廣泛應用。廣義狀態空間平均法則使用傅里葉級數中的開關頻率項來表示LLC變換器狀態變量的平均值,便于進行小信號擾動和線性化[11-12]。擴展描述函數法作為一種連續模型的建模方法,大致建模過程為將變換器中的狀態變量進行傅里葉級數展開,再通過諧波平衡原理建立大信號模型[13-15],由于這類模型通常具有較高的階數,一種LLC變換器等效電路的簡化方法被提出以實現模型的降階[16-17]。其中,通過添加高次諧波分量,可以提高廣義狀態空間平均法和擴展描述函數法建模的精度,但是,諧波分量次數的增加會提高狀態方程組的求解難度[18]。因此,基于S型函數的建模方法被應用于諧振變換器[19],其利用連續的S型函數將離散的開關變量近似為連續變量,從而建立諧振變換器的大信號連續模型。然而,由于LLC變換器的特性,其運行于斷續電流模式時具有兩種等效諧振拓撲[20],而基于S型函數的建模方法難以在連續模型中同時準確地描述出這兩種等效電路拓撲。因此,這種建模方法僅當LLC變換器運行于連續電流模式時較為準確。

為了以連續模型準確地實現LLC變換器多種運行模式的統一描述,本文提出了一種基于拓展型雙曲正切函數的全橋LLC變換器大信號建模方法。本文的主要特點和優點如下:

1)基于符號函數、絕對值函數和本文定義的切換函數,建立了一種全橋LLC變換器大信號非連續模型。該模型能夠準確描述全橋LLC變換器斷續電流模式中兩種等效諧振拓撲的自動切換,進而實現全橋LLC變換器斷續電流模式和連續電流模式的統一描述。

2)構建了含陡度因子/脈沖系數的雙曲正切函數并給出了參數設計方法,基于此將大信號模型連續化。得益于上述拓展型雙曲正切函數,提出的連續模型可設置死區時間、器件開關延時參數,從而提高了建模準確性。此外,該模型具有較低的階數,并能提供次級側功率器件高準確性的狀態信息,便于作為同步整流的參考。

1 全工作頻率范圍統一的大信號模型

圖1 全橋LLC變換器拓撲結構

1.1 連續電流模式:

式中,sgn()為符號函數。

1.2 斷續電流模式:

2 基于拓展型雙曲正切函數的大信號連續模型

上述動態模型能夠準確地描述全橋LLC變換器在全工作頻率范圍內的不同運行模式。然而,該模型是非連續的,其中存在符號函數等多種非連續函數。對于電力電子變換器的特性分析、實時仿真和控制器設計而言,這類模型意味著連續系統方法無法使用,大信號模型的連續化十分關鍵。本文將上述非連續模型連續化,基于拓展型雙曲正切函數建立了一種全橋LLC變換器的大信號連續模型。

2.1 含陡度因子的雙曲正切函數

雙曲正切函數的值域為(-1,1),向自變量x中注入陡度因子b,本文提出一種含陡度因子b的雙曲正切函數。不同陡度因子取值的雙曲正切函數與sgn函數的曲線如圖4所示。

因此,絕對值函數同樣可以利用含陡度因子的雙曲正切函數近似逼近。

基于周期型sgn函數的和含不同陡度因子的周期型雙曲正切函數的曲線如圖5所示。顯然,與基于周期型sgn函數的不同,通過含陡度因子的周期型雙曲正切函數描述的并非標準的方波,其根據陡度因子的取值在切換過程中具有不同的斜率。

參考實際應用中對于DC-DC變換器輸出電壓紋波的要求,以其不大于1%的輸出電壓為例,將兩組開關管Q1、Q4或Q2、Q3完全導通進入穩態的近似條件定義為

LLC諧振變換器兩組開關管的開關狀態切換過程如圖6所示。

式中,為電壓紋波率。根據式(18)可求得陡度因子。實際應用中陡度因子通常取值較高,可將最后結果適當取整以簡化后續分析和計算。

2.2 含脈沖系數的雙曲正切函數

為了以連續函數準確描述LLC變換器連續模式和斷續模式的自動切換,本文構建了一種含脈沖系數的雙曲正切函數,即

式中,k為脈沖系數,其決定了函數值在0和1之間的切換速度。當脈沖系數取不同值時拓展型雙曲正切函數與函數的函數圖像如圖7所示,當脈沖系數取值較高時有。

進而,可以通過含脈沖系數的雙曲正切函數實現全橋LLC變換器的連續模式和斷續模式的統一描述。

基于上述含陡度因子和含脈沖系數的拓展型雙曲正切函數,全橋LLC變換器的大信號連續模型可表示為

本文提出的大信號連續模型的應用條件為一次側逆變拓撲為全橋、二次側整流拓撲無限制的LLC變換器,該模型能夠實現全橋LLC變換器在連續電流模式和斷續電流模式運行特性的準確描述,即在全橋LLC變換器的全工作頻率范圍內均適用。值得一提的是,為了進一步提高LLC變換器的開關頻率和功率密度,平面型變壓器開始逐漸被應用于LLC變換器中。而平面型變壓器的雜散電容對于LLC變換器的輕載特性具有一定影響。本文提出的大信號建模方法并未考慮該雜散電容。因此,本文提出的大信號連續模型在對采用平面型變壓器的LLC變換器的輕載工況進行分析時,其精確度會有所降低,這也是后續針對帶有平面型變壓器的LLC變換器建模的研究方向。

2.3 模型誤差分析

3 仿真與實驗

表1 LLC變換器關鍵參數

3.1 仿真驗證

本文利用Matlab分別建立全橋LLC變換器的基于擴展描述函數法的大信號模型和基于拓展型雙曲正切函數的大信號模型,利用Simulink搭建圖1所示的全橋LLC變換器仿真模型。

圖8~圖10分別為上述模型在小于、等于和大于諧振頻率時電流、電壓的關鍵波形,D為流過二次側同步整流器的電流。其中,基于拓展型雙曲正切函數的大信號連續模型依據式(21)~式(24)建立,基于擴展描述函數法的模型依據文獻[14]所提出的方法建立,具體模型為

圖10 全橋LLC變換器仿真模型關鍵波形()

以LLC變換器Simulink模型的仿真結果為基準對提出的大信號連續模型在不同頻率區域(大于、等于或小于諧振頻率)的方均誤差進行計算,結果見表2。可以看出,在全工作頻率范圍內,本文提出的大信號連續模型均具有較高的準確度。相比于輕載工況,重載工況下本文提出的模型對LLC變換器運行特性的描述更為準確。

仿真時間=10ms,LLC變換器大信號連續模型與本文提出的模型的計算耗時和內存占用情況見表3。雖然基于本文提出的大信號連續模型的計算耗時和資源消耗略高于文獻[14]和文獻[17]中的大信號連續模型。但從圖8~圖10中可以看出,本文提出的大信號連續模型相比于現有LLC變換器大信號連續模型則具有更高的準確度。

表2 提出模型/Simulink模型-方均根誤差

Tab.2 Proposed model/Simulink model – RMSE

表3 模型計算耗時和內存占用

圖11驗證了在輸入電壓突變條件下,提出的大信號模型的動態特性。在=3ms時,輸入電壓由600V降至400V,輸出電壓做出相應響應。由圖11可知,基于拓展型雙曲正切函數的大信號模型的瞬態響應與基于Simulink的電路仿真模型結果基本一致,即所提出LLC變換器大信號模型可以準確描述系統的動態特性。

圖11 提出的大信號模型動態特性圖

本文提出的基于拓展型雙曲正切函數的大信號模型是連續的,便于應用連續系統的方法對LLC變換器進行特性分析、實時仿真和控制器設計。雖然基于擴展描述函數法建立的連續模型與本文提出的大信號連續模型相比更易于進行線性化和小信號擾動,但其模型的階數通常較高,而本文提出的大信號模型則具有更低的階數,并具有更高的準確度。

3.2 實驗驗證

圖12 實驗環境

圖13 全橋LLC變換器實驗波形

圖14 全橋LLC變換器實驗波形

由圖13~圖15可知,基于拓展型雙曲正切函數的大信號模型和基于Simulink的模型的仿真結果與全橋LLC變換器的實際運行狀態相同,進一步驗證了本文提出大信號模型的有效性。上述仿真和實驗結果驗證了本文提出模型的準確性和有效性,所提出的基于拓展型雙曲正切函數的大信號模型實現了全橋LLC變換器全工作頻率范圍內靜/動態特性的準確描述。

圖15 全橋LLC變換器實驗波形

4 結論

為了以連續模型準確地實現對LLC變換器多種運行模式的統一描述,本文提出了一種基于拓展型雙曲正切函數的全橋LLC變換器大信號建模方法。基于此方法建立的大信號模型不包含任何非連續函數,進而可以直接應用連續系統的方法對LLC變換器進行動態特性分析、實時仿真和控制器設計。同時,該模型考慮了死區時間、器件開關延時等參數影響,能夠準確地實現LLC變換器連續電流模式與斷續電流模式的統一描述。相比于現有的LLC變換器大信號連續模型,該模型具有更低的階數,并且具有更高的準確度。此外,基于拓展型雙曲正切函數的大信號模型在時域中能夠給出LLC變換器二次電流高準確性的狀態信息,便于同步整流方案的設計。通過仿真和實驗驗證了所提出大信號建模方法的準確性和有效性。

[1] 孫加祥, 吳紅飛, 湯欣喜, 等. 基于整流側輔助調控的交錯并聯LLC諧振變換器[J]. 電工技術學報, 2021, 36(10): 2072-2080.

Sun Jiaxiang, Wu Hongfei, Tang Xinxi, et al. Interleaved LLC resonant converter with auxiliary regulation of rectifier[J]. Transactions of China Electrotechnical Society, 2021, 36(10): 2072-2080.

[2] 陳夢穎, 王議鋒, 陳慶, 等. 一種可變結構型高效寬增益多諧振軟開關直流變換器[J]. 電工技術學報, 2021, 36(20): 4225-4236.

Chen Mengying, Wang Yifeng, Chen Qing, et al. A variable topology multi-resonant soft-switching DC-DC converter with high efficiency and wide gain[J]. Transactions of China Electrotechnical Society, 2021, 36(20): 4225-4236.

[3] 周國華, 冷敏瑞, 李媛, 等. 開關變換器及其控制環路的建模綜述[J]. 中國電機工程學報, 2020, 40(1): 183-199, 386.

Zhou Guohua, Leng Minrui, Li Yuan, et al. A review on modeling of switching converters and their control loops[J]. Proceedings of the CSEE, 2020, 40(1): 183-199, 386.

[4] Bhat A K S. A unified approach for the steady-state analysis of resonant converters[J]. IEEE Transactions on Industrial Electronics, 1991, 38(4): 251-259.

[5] De Simone S, Adragna C, Spini C, et al. Design-oriented steady-state analysis of LLC resonant converters based on FHA[C]//2006International Symposium on Power Electronics, Electrical Drives, Automation and Motion, Taormina, Italy, 2006: 200-207.

[6] 丁超, 李勇, 姜利, 等. 電動汽車直流充電系統LLC諧振變換器軟開關電壓邊界分析[J]. 電工技術學報, 2022, 37(1): 3-11.

Ding Chao, Li Yong, Jiang Li, et al. Analysis of soft switching voltage boundary of LLC resonant converter for EV DC charging system[J]. Transactions of China Electrotechnical Society, 2022, 37(1): 3-11.

[7] Ivensky G, Bronshtein S, Abramovitz A. Approximate analysis of resonant LLC DC-DC converter[J]. IEEE Transactions on Power Electronics, 2011, 26(11): 3274-3284.

[8] Liu Jianqiang, Zhang Jiepin, Zheng T Q, et al. A modified gain model and the corresponding design method for an LLC resonant converter[J]. IEEE Transactions on Power Electronics, 2017, 32(9): 6716-6727.

[9] Agarwal V, Bhat A K S. Small signal analysis of the LCC-type parallel resonant converter using discrete time domain modeling[C]//Proceedings of 1994 Power Electronics Specialist Conference-PESC'94, Taipei, China, 1994: 805-813.

[10] Song Zhanfeng, Zhou Fengjiao. Observer-based predictive vector-resonant current control of permanent magnet synchronous machines[J]. IEEE Transactions on Power Electronics, 2019, 34(6): 5969-5980.

[11] Qin Hengsi, Kimball J W. Generalized average modeling of dual active bridge DC-DC converter[J]. IEEE Transactions on Power Electronics, 2012, 27(4): 2078-2084.

[12] Shang Fei, Wu Haowen, Niu Geng, et al. Dynamic analysis and control approach for a high-gain step-up converter for electrified transportation[J]. IEEE Transactions on Transportation Electrification, 2017, 3(3): 656-667.

[13] Menke M F, Seidel á R, Tambara R V. LLC LED driver small-signal modeling and digital control design for active ripple compensation[J]. IEEE Transactions on Industrial Electronics, 2019, 66(1): 387-396.

[14] Buccella C, Cecati C, Latafat H, et al. Observer-based control of LLC DC/DC resonant converter using extended describing functions[J]. IEEE Transactions on Power Electronics, 2015, 30(10): 5881-5891.

[15] 徐玉珍, 郭躍森, 林維明. 一種具有自動均流特性的并聯LLC諧振變換器[J]. 中國電機工程學報, 2018, 38(22): 6671-6683.

Xu Yuzhen, Guo Yuesen, Lin Weiming. A parallel LLC resonant converter with automatic current sharing characteristics[J]. Proceedings of the CSEE, 2018, 38(22): 6671-6683.

[16] Tian Shuilin, Lee F C, Li Qiang. Equivalent circuit modeling of LLC resonant converter[C]//2016 IEEE Applied Power Electronics Conference and Exposition, Long Beach, CA, USA, 2016: 1608-1615.

[17] Tian Shuilin, Lee F C, Li Qiang. Equivalent circuit modeling of LLC resonant converter[J]. IEEE Transactions on Power Electronics, 2020, 35(8): 8833-8845.

[18] Mansour A, Hajer M, Faouzi B, et al. Analysis and modeling of LLC resonant converter used in electric vehicle[C]//2019 International Conference on Advanced Systems and Emergent Technologies (IC_ASET), Hammamet, Tunisia, 2019: 357-362.

[19] Lu Yimin, Huang Xianfeng, Huang Yizheng, et al. Sigmoid function model for a PFM power electronic converter[J]. IEEE Transactions on Power Electronics, 2020, 35(4): 4233-4241.

[20] Hsu J D, Ordonez M, Eberle W, et al. LLC synchronous rectification using resonant capacitor voltage[J]. IEEE Transactions on Power Electronics, 2019, 34(11): 10970-10987.

[21] Kim J H, Kim C E, Kim J K, et al. Analysis on load-adaptive phase-shift control for high efficiency full-bridge LLC resonant converter under light-load conditions[J]. IEEE Transactions on Power Electronics, 2016, 31(7): 4942-4955.

[22] 焦健, 郭希錚, 游小杰, 等. LLC諧振變換器的改進型電流解析方法[J]. 電工技術學報, 2021, 36(23): 5002-5013.

Jiao Jian, Guo Xizheng, You Xiaojie, et al. An improved current analytical method for LLC resonant converter[J]. Transactions of China Electrotechnical Society, 2021, 36(23): 5002-5013.

[23] 袁義生, 易塵宇, 彭能. L-R復合調制T型半橋LCC諧振變換器[J]. 電工技術學報, 2022, 37(4): 892-904.

Yuan Yisheng, Yi Chenyu, Peng Neng. T-type half-bridge LCC resonant converter with L-R composite modulation[J]. Transactions of China Electrotechnical Society, 2022, 37(4): 892-904.

Large-Signal Modeling Method of Full-Bridge LLC Converter Based on Extended Hyperbolic Tangent Function

Sun Chenghao Sun Qiuye Wang Rui Hu Jie

(College of Information Science and Engineering Northeastern University Shenyang 110819 China)

The accuracy of the continuous model is very important for the characteristic analysis, real time simulation and controller design of power electronic converter. LLC resonant converter operates in continuous current mode and discontinuous current mode respectively in different operating frequency intervals due to its unique cavity structure. Therefore, it is difficult to establish a continuous model which can realize the unified description of various operating modes of LLC resonant converter. Therefore, a large signal modeling method for full bridge LLC resonant converter based on extended hyperbolic tangent function is proposed to establish an accurate large signal continuous model. Firstly, the operating characteristics of full bridge LLC resonant converters in each frequency range are analyzed. Based on this analysis, a large signal discontinuity model of full bridge LLC resonant converter is established by using the symbolic function, the absolute value function and the defined operation mode switching function. The discontinuity model can realize the unified description of continuous current mode and discontinuous current mode of full bridge LLC resonant converter. Then, two kinds of extended hyperbolic tangent functions including steepness factor and pulse coefficient are constructed, and the standardized selection methods of steepness factor and pulse coefficient are provided. By changing the steepness factor and pulse coefficient, the dead time and device switching delay time of the full bridge LLC resonant converter can be approximated with high precision, thus realizing accurate large signal modeling. Based on the extended hyperbolic tangent function, the discontinuous model is continuous, and the large signal continuous model of the full bridge LLC converter is established. Based on this model, the method of continuous system can be applied directly to realize the characteristic analysis and controller design of converter. At the same time, compared with the existing LLC resonant converter large signal continuous model, the model has lower order and higher accuracy. In addition, the large signal continuous model established in this paper can provide highly accurate switching information of the secondary power components of the full bridge LLC resonant converter in the time domain, which is convenient to be used as a reference for the design of synchronous rectifier controller. Finally, the accuracy and effectiveness of the proposed large signal continuous model are verified by building the MATLAB / Simulink simulation model and the experimental prototype with rated power of 5kW. In the full range of operating frequency, the large signal continuous model is presented with high accuracy. Compared with the light load condition, the model proposed in the heavy load condition can describe the operating characteristics of the full bridge LLC resonant converter more accurately. The computing time and resource consumption are slightly higher than the existing large signal continuous model, but the proposed large signal continuous model has higher accuracy. Although the continuous model based on the extended description function method is easier to carry out linearization and small signal perturbation than the proposed large signal continuous model, the order of the model is usually higher, while the proposed large signal continuous model proposed has a lower order. In addition, the proposed large signal continuous model can accurately describe the dynamic characteristics of the full bridge LLC resonant converter.

Full-bridge LLC converter, continuous system, large-signal model, extended hyperbolic tangent function

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

TM46

孫城皓 男,1996年生,博士研究生,研究方向為諧振變換器的建模與控制。E-mail:1910243@stu.neu.edu.cn

孫秋野 男,1977年生,教授,博士生導師,研究方向為能源互聯網的建模與優化運行、電力電子化電力系統等。E-mail:sunqiuye@ise.neu.edu.cn(通信作者)

國家重點研發計劃變革性課題(2018YFA0702200)和國家自然科學基金(62073065)資助項目。

2021-06-30

2022-04-21

(編輯 赫蕾)

猜你喜歡
信號方法模型
一半模型
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
重要模型『一線三等角』
完形填空二則
重尾非線性自回歸模型自加權M-估計的漸近分布
基于FPGA的多功能信號發生器的設計
電子制作(2018年11期)2018-08-04 03:25:42
3D打印中的模型分割與打包
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
基于LabVIEW的力加載信號采集與PID控制
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
主站蜘蛛池模板: 国产成人精品三级| 天堂在线亚洲| jizz亚洲高清在线观看| 国产噜噜在线视频观看| 亚洲天堂精品在线| 亚洲swag精品自拍一区| 亚洲五月激情网| 欧美日韩国产综合视频在线观看| 欧美视频在线不卡| 欧美影院久久| 97超碰精品成人国产| 久久久久久午夜精品| 欧美精品亚洲精品日韩专区va| 美女内射视频WWW网站午夜 | 亚洲成人网在线播放| 亚洲天堂免费| 国产一区二区三区在线观看免费| 亚洲欧美人成人让影院| 久久频这里精品99香蕉久网址| 人妻无码一区二区视频| 91在线精品免费免费播放| 久久香蕉国产线看观看亚洲片| 日韩精品一区二区三区swag| 99精品视频在线观看免费播放| 亚洲中文制服丝袜欧美精品| 久久国产精品夜色| 久久久久无码精品国产免费| 强乱中文字幕在线播放不卡| 欧美精品三级在线| 伊人福利视频| 中国特黄美女一级视频| 亚洲欧洲日本在线| 日韩一区二区在线电影| 欧美三级自拍| 久久久久久久久亚洲精品| 亚洲欧美极品| 丁香六月综合网| 国产在线欧美| 国产va免费精品观看| 日韩欧美中文字幕在线精品| 青青草欧美| 国产精品永久在线| 国产精品视频观看裸模| www.av男人.com| 国产成人a在线观看视频| 无码aaa视频| 国产成人在线小视频| 国产精品网址在线观看你懂的| 亚洲男人的天堂久久精品| 日韩一区精品视频一区二区| 欧美97欧美综合色伦图| 免费毛片全部不收费的| 91福利免费视频| 香蕉伊思人视频| 欧美伦理一区| 人妻精品久久无码区| 免费无遮挡AV| 国产理论一区| 一级毛片免费不卡在线| 国产欧美在线观看精品一区污| 干中文字幕| 综合久久五月天| 狠狠ⅴ日韩v欧美v天堂| 日本国产在线| 亚洲中文字幕手机在线第一页| 免费一级毛片在线播放傲雪网| 亚洲aⅴ天堂| 在线网站18禁| 欧美亚洲国产视频| 91免费观看视频| 成人欧美日韩| 国产在线八区| 久久人人97超碰人人澡爱香蕉| 亚洲天堂网在线播放| 久草性视频| 国产Av无码精品色午夜| 99视频精品全国免费品| 91久久国产综合精品女同我| 国产成人久久777777| 91九色国产在线| 午夜无码一区二区三区| 亚洲欧美精品在线|