亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        環(huán)境溫度對電勵磁雙凸極發(fā)電機輸出特性的影響

        2021-04-13 03:23:46蔣思遠王開淼馮小寶
        電源學(xué)報 2021年2期
        關(guān)鍵詞:磁路磁阻鐵心

        蔣思遠,周 波,王開淼,馮小寶

        (江蘇省新能源發(fā)電與電能變換重點實驗室(南京航空航天大學(xué)),南京210016)

        多電飛機是當前航空技術(shù)發(fā)展的重要趨勢,為解決多電飛機對發(fā)電機功率的巨大需求,新型航空發(fā)電機技術(shù)受到了廣泛關(guān)注[1-4]。 電勵磁雙凸極發(fā)電機DSEG(doubly salient electro-magnetic generator)具有結(jié)構(gòu)簡單、無永磁體和高溫高速運行能力強的優(yōu)點,且勵磁可調(diào),故障時易于滅磁,尤其是發(fā)電運行中無需位置傳感器,電壓控制方便,發(fā)電可靠性高[5-8],因此DSEG 可構(gòu)成一種極具發(fā)展?jié)摿Φ男滦秃娇针娫聪到y(tǒng)。航空應(yīng)用環(huán)境對溫度適應(yīng)性要求高,大溫差的環(huán)境溫度將引起電機材料的導(dǎo)電導(dǎo)磁性能發(fā)生改變,對電機的各項運行性能產(chǎn)生較大影響,因此研究環(huán)境溫度變化對DSEG 輸出特性的影響十分必要。電機的溫度特性分析研究涉及材料科學(xué)、熱科學(xué)和電磁學(xué)等多個領(lǐng)域,本文主要對DSEG 不同溫度下的電磁特性進行分析。

        溫度對電機材料的影響主要包括對永磁體剩磁及退磁的影響、對銅導(dǎo)線電阻的影響、對鐵心磁導(dǎo)率的影響、對軸承潤滑脂及絕緣材料的影響等[9-16]。 溫度對電機電磁性能影響的研究,目前主要集中在異步電機和永磁電機的研究上。 異步電機方面,文獻[17]通過建立數(shù)學(xué)模型,分析了高溫下電阻變化對高頻異步電動機性能的影響,為高溫異步電動機的設(shè)計提供了參考;文獻[18]通過實驗和數(shù)值分析對高環(huán)境溫度下感應(yīng)電機的鐵損特性進行研究,結(jié)果表明鐵損隨著溫度的升高而降低,文中還對環(huán)境溫度變化時電機鐵心齒部磁密進行研究,發(fā)現(xiàn)齒部磁密幾乎是恒定值,與DSEG 磁路局部飽和嚴重的情況差異較大。 永磁電機方面,文獻[19]對一種可在-60~300 °C 的寬溫度范圍內(nèi)使用的永磁同步電機進行設(shè)計和性能分析,研究了溫度變化對剩余磁通密度、銅損、鐵損和風(fēng)阻損耗等的影響,并對電機在不同溫度下進行了效率計算;文獻[20-22]研究了低溫對增程式電動汽車用永磁同步電機性能的影響,并建立了測試平臺, 研究內(nèi)容包括電機起動性能、運行效率及高低溫環(huán)境下的負載能力等。

        文獻[17-22]仿真和數(shù)值分析電機溫度特性時忽略了鐵心材料性能隨溫度的改變,但與之不同的是,DSEG 由于雙凸極結(jié)構(gòu),磁場非線性嚴重,局部磁路嚴重飽和,不同溫度下鐵磁材料特性的變化將引起電機參數(shù)明顯變化,故需要考慮鐵心材料對電機性能的影響。DSEG 通過引入勵磁繞組實現(xiàn)他勵,而航空發(fā)動機轉(zhuǎn)速變化范圍大,寬轉(zhuǎn)速范圍發(fā)電運行對DSEG 勵磁機構(gòu)要求較高,環(huán)境溫度的升高也使勵磁繞組電阻增大,進一步增大勵磁損耗,帶來額外的勵磁負擔(dān),故DSEG 溫度特性分析還需考慮溫度對勵磁機構(gòu)的影響。目前針對電勵磁雙凸極電機的溫度場分布的文獻尚比較少見,文獻[23]對一臺18/12 極DSEG 進行了定子軛高優(yōu)化, 建立優(yōu)化尺寸后的電機三維熱有限元模型, 分析了所提DSEG的溫度場分布情況,以保證在安全范圍內(nèi)的最大溫升, 研制了全封閉自冷樣機DSEG 并進行了測試。而專門針對DSEG 不同環(huán)境溫度下的電機特性分析尚未檢索到相關(guān)文獻,缺乏具體的性能分析方面的深入研究,使得航空應(yīng)用環(huán)境下DSEG 的設(shè)計和應(yīng)用受到了一定的限制。

        為充分發(fā)揮DSEG 高溫運行能力強的優(yōu)勢,本文對環(huán)境溫度改變時的DSEG 發(fā)電性能進行研究。采用等效磁路法和有限元法,結(jié)合不同溫度下的鐵心材料磁導(dǎo)率曲線,考慮DSEG 不同位置時的磁路飽和程度,指出電機磁場和繞組電感隨溫度的變化規(guī)律; 對不同溫度下電機的負載特性進行分析;將不同溫度下的鐵心材料和繞組銅材進行組合,通過有限元建模,研究了溫升導(dǎo)致的電樞繞組電阻變化和鐵心磁阻的變化單一變量對電機外特性的影響規(guī)律;最后進行了電機溫度特性實驗,以驗證分析和仿真的正確性。實驗時采用直流電源靜止熱實驗的方式加熱樣機,使仿真與實驗溫度更為接近。 通過磁路計算、仿真分析和實驗的方法確保結(jié)論的有效性, 填補了DSEG 溫度特性分析方面的不足,為航空應(yīng)用環(huán)境下DSEG 的使用和設(shè)計優(yōu)化提供了參考。

        1 DSEG 結(jié)構(gòu)及工作原理

        圖1 為12/8 極DSEG 結(jié)構(gòu),定、轉(zhuǎn)子均為凸極齒槽結(jié)構(gòu),轉(zhuǎn)子上無繞組與永磁體,通過定子上的勵磁繞組進行勵磁。 DSEG 發(fā)電運行時, 三相繞組通常外接三相全橋整流電路,如圖2 所示。

        理想化DSEG 輸出電壓模型多忽略繞組電阻,但本文對不同環(huán)境溫度下DSEG 的發(fā)電特性進行研究,其電阻電感均受溫度影響較大,需考慮鐵心磁導(dǎo)率和繞組電阻率隨溫度的變化。電機的工作環(huán)境溫度相對繞組導(dǎo)體銅材料本身而言仍然屬于低溫段,一般認為在電機工作溫度下,銅材料的電阻隨溫度變化近似存在著線性關(guān)系,故考慮溫度變化的DSEG 電壓方程為

        式中:rj和rj0分別為溫度t 及標準溫度t0下的電阻,其中,j 為三相繞組和勵磁繞組j=a,b,c 和f;α 為材料在標準溫度下的電阻率溫度系數(shù),材料銅[6]取為0.003 93 ℃-1;ψa、ψb、ψc和ψf分別為三相繞組和勵磁繞組磁鏈;ia、ib、ic和if分別為三相繞組和勵磁繞組電流。

        圖1 12/8 極DSEG 結(jié)構(gòu)Fig. 1 Structure of 12/8-pole DSEG

        圖2 DSEG 功率電路Fig. 2 DSEG power circuit

        DSEG 發(fā)電運行時, 相繞組電阻隨溫度增加而增加,將使其輸出電壓降低;而勵磁繞組電阻的增大,使得相同的勵磁電流下需要更高的電壓。 鐵心材料導(dǎo)磁性能隨溫度的變化將使各繞組磁鏈發(fā)生變化,進而影響DSEG 輸出特性,。

        2 考慮溫度影響的電機磁路模型

        以A 相為例,傳統(tǒng)三相12/8 極結(jié)構(gòu)的電勵磁雙凸極電機典型位置的磁力線分布和等效磁路如圖3所示。圖中Ff和Fp分別為勵磁繞組和p(p=a,b,c)相電樞繞組產(chǎn)生的磁勢,Rsy1和Rsy2為定子軛部磁阻,Rry為轉(zhuǎn)子軛部磁阻,Rsp和Rrp分別為定子和轉(zhuǎn)子齒部磁阻。 Rg和Rsl分別為氣隙磁阻和漏磁阻,各部位的磁阻計算公式如表1 所示。

        表1 中:g 為氣隙長度;μFe為鐵心材料磁導(dǎo)率;μ0為真空磁導(dǎo)率;各處的磁路長度l 和截面積S 的下標采用和磁阻R 相同的表示方式,另外引入磁導(dǎo)率溫度系數(shù)kT來描述溫度對鐵心磁阻的影響。

        圖3 典型位置的磁力線分布和等效磁路Fig. 3 Magnetic field distribution and equivalent magnetic circuit at typical locations

        表1 磁阻計算Tab. 1 Calculation of magnetic resistance

        由于雙凸極電機定轉(zhuǎn)子鐵心和開關(guān)磁阻電機類似, 故借鑒開關(guān)磁阻電機氣隙截面積的計算公式, 考慮邊緣效應(yīng), DSEM 中一個定子極氣隙截面積[24]Sg可取為

        式中:βs和βr分別為定子極和轉(zhuǎn)子極對應(yīng)的機械角;D 為氣隙直徑;lFe為鐵心軸向長度。

        同樣,以A 相為例,對三相DSEG 典型位置的主磁路進行分析,如圖3 中粗線所示。 忽略漏磁,空載時,F(xiàn)A=FB=FC=0,根據(jù)安培環(huán)路定理,在圖3(a)顯示了A 相轉(zhuǎn)子極滑入時刻的磁路模型,勵磁磁勢Ff為

        式中:Nf和if分別為勵磁繞組的匝數(shù)和電流;ψA為A 相主磁通。

        圖3(b)顯示了A 相定、轉(zhuǎn)子極對齊時刻的磁路模型,勵磁磁勢Ff為

        文獻[9]對不同溫度下無取向電工鋼進行了系統(tǒng)實驗測試,其相對磁導(dǎo)率曲線如圖4 所示,可見,溫度對低磁感應(yīng)強度下和磁飽和情況下的硅鋼片磁導(dǎo)率幾乎沒有影響,對中等磁感應(yīng)強度下的硅鋼片磁導(dǎo)率,隨著溫度的上升,磁導(dǎo)率下降。故在分析溫度影響的電機磁路模型時,需要結(jié)合電機鐵心中的磁密分布進行分析。

        圖4 溫度對相對磁導(dǎo)率的影響Fig. 4 Effect of temperature on relative permeability

        DSEG 電樞繞組磁鏈周期性變化產(chǎn)生感應(yīng)電勢,磁鏈的最大值和最小值直接影響其輸出功率的大小。在A 相轉(zhuǎn)子極滑入和滑出時刻(如圖3(a)、(c)所示),定、轉(zhuǎn)子齒部均會局部飽和的情況。 所以轉(zhuǎn)子極開始滑入定子極階段,溫度上升時,極部的磁飽和處磁導(dǎo)率幾乎不變,可認為kT≈1;隨著定、轉(zhuǎn)子極正對面積的增加, 鐵心的磁密飽和程度降低,此時可認為kT<1,且越來越小,這意味著磁導(dǎo)率下降,磁路磁阻增大,故溫度對各相磁鏈影響最大的位置為定轉(zhuǎn)子極正對位置(如圖3(b)所示)。 滑出時的情況和滑入時類似,不再贅述。 由電感公式L=Nψ/i可知,溫度引起各相磁阻的改變也會使得各相電感受到影響,影響最大的位置會出現(xiàn)在定轉(zhuǎn)子極對齊位置。此外,軟磁材料種類繁多,目前對磁性材料性能參數(shù)隨溫度變化規(guī)律認識的不足,使得材料具體溫度的各項性能參數(shù)確定存在較大困難,雖然引入溫度磁導(dǎo)率溫度系數(shù)kT可描述溫度對鐵芯磁阻的影響, 但由于該系數(shù)的非線性和不同材料的差異性,使得難以定量分析。 下節(jié)將采用有限元法進行仿真計算,充分考慮不同溫度下軟磁材料性能的變化對電機性能的影響。

        3 電機溫度特性仿真分析

        表2 為12/8 極DSEG 的結(jié)構(gòu)參數(shù),根據(jù)其結(jié)構(gòu)參數(shù)進行有限元建模,并將不同溫度下硅鋼片的導(dǎo)磁性能和繞組銅材的導(dǎo)電參數(shù)代入有限元軟件進行分析和求解,研究DSEG 的溫度特性,有限元計算模型如圖5 所示。

        本文在進行DSEG 發(fā)電性能的仿真計算時,主要考慮溫升對電機電磁性能的影響,假設(shè)條件如下:

        (1)電機在本文涉及的溫度條件下,各部分絕緣和潤滑條件良好,能夠正常穩(wěn)定工作;

        (2)電機材料的受熱膨脹影響忽略不計,各部件在不同溫度下不發(fā)生形變;

        (3)電機機械及雜散損耗與常規(guī)溫度下額定工作時相同;

        (4)電機在高溫環(huán)境下是一個等溫體,各部分溫度相同且為環(huán)境溫度。

        表2 電機的主要結(jié)構(gòu)參數(shù)Tab. 2 Main structural parameters of DSEG

        圖5 DSEG 有限元模型Fig. 5 Finite element model of DSEG

        3.1 靜磁場與電感特性分析

        勵磁電流為5 A 時,不同溫度下的DSEG 相自感曲線及相繞組與勵磁繞組互感曲線如圖6 所示,由于二者變化趨勢幾乎一致,以相自感曲線為例進行分析,可以看出,溫度的變化對自感的最小值幾乎沒有影響,對最大值影響較為明顯,且隨著溫度的升高,最大值逐漸減小。由磁路分析可知,從轉(zhuǎn)子極開始滑入到定轉(zhuǎn)子極對齊的過程中,溫度的升高將導(dǎo)致磁路磁阻逐漸增大, 故自感隨之逐漸減小,且隨著定轉(zhuǎn)子正對面積的增加該差異趨于明顯,仿真結(jié)果與磁路分析的結(jié)論一致。

        3.2 不同溫度下DSEG 負載特性分析

        圖7 是電機勵磁電流5 A、轉(zhuǎn)速1 000 r/min 時,不同溫度下的輸出電壓和輸出功率,可以看出,隨著溫度的上升,輸出電壓曲線下移,輸出功率降低。 一方面由于溫度升高導(dǎo)致電樞繞組電阻增大, 相當于增加了發(fā)電機內(nèi)阻,影響發(fā)電機功率輸出;另一方面由于溫度升高,鐵心磁阻發(fā)生變化,進而導(dǎo)致電機自感、互感減小,引起電壓降低。

        圖7 不同溫度下的輸出電壓和輸出功率Fig. 7 Output voltage and output power at different temperatures

        為分析溫升導(dǎo)致的電樞繞組電阻變化和鐵心磁阻的變化對電機外特性的影響,以20 ℃和200 ℃的模型為基準,將20 ℃和200 ℃鐵心材料和電樞繞組電阻進行組合,分別進行仿真分析對比,有限元材料設(shè)置組合如表3 所示。

        表3 有限元材料設(shè)置組合Tab. 3 Combination of material setting based on FEM

        分兩類對有限元進行對比,分析結(jié)果如圖8 所示。其中,一類為保持電樞繞組電阻不變,改變鐵心材料,即將模型1 和3(或模型2 和4)進行對比,可以看出,輸出電壓和輸出功率曲線略微下移,原因在于鐵心材料磁阻隨溫度上升而變大,導(dǎo)致電感變化值變小,相同勵磁磁勢下的輸出電壓降低。 第二類為保持鐵心材料不變,改變電樞繞組電阻,即將模型1 和4(或模型2 和3)進行對比,可以看出,電機輸出電壓和輸出功率下降明顯, 說明溫度升高時,電樞電阻的增加是導(dǎo)致電機輸出電壓降低的主要原因。

        圖8 電阻變化和鐵心磁阻變化的影響Fig. 8 Effects of changes in winding resistance and core magnetic resistance

        4 實驗驗證

        由于理論分析和仿真條件為理想化假設(shè),僅依靠理論研究和仿真分析電機的溫度特性是不夠的。事實上,一些電機在不同環(huán)境溫度下的運行特性只有通過實驗才能得到, 故主要對12/8 三相DSEG采用全橋整流方式的溫度特性進行實驗驗證。由于實驗室條件所限,電機熱態(tài)電磁特性分析在電機溫升實驗的基礎(chǔ)上進行,分別對樣機在室溫環(huán)境下及溫升實驗后的熱狀態(tài)下進行特性分析,并作對比研究。 電機測試平臺如圖9 所示,主要包含三相12/8極DSEG、RT-LAB 系統(tǒng)控制平臺、 負載電阻箱、直流電源、采樣和驅(qū)動電路等。

        圖9 實驗平臺Fig. 9 Experimental platform

        實驗過程中, 電機外部采用熱成像儀進行測溫,電機繞組溫升采用電阻法進行測試,通過伏安法計算繞組阻值,再得出繞組溫升,表示為

        式中:θ1為測量繞組(冷態(tài))初始電阻時的溫度;θ2為熱試驗結(jié)束時繞組的溫度;θa為熱試驗結(jié)束時冷卻介質(zhì)(本文為空氣)溫度;R1為溫度為θ1(冷態(tài))時的繞組電阻;R2為熱試驗結(jié)束時的繞組電阻;k 為導(dǎo)體材料在0 ℃時電阻溫度系數(shù)的倒數(shù), 其中材料銅的kCu=235。

        4.1 基于動態(tài)溫升實驗的特性分析

        在轉(zhuǎn)速1 000 r/min、 勵磁電流6 A 時對DSEG進行了空載測溫實驗, 實驗初始環(huán)境溫度27.3 ℃,電機運行120 min 基本達到穩(wěn)態(tài)溫度。 圖10 顯示了DSEG 勵磁繞組的電阻變化和溫升曲線。 實驗過程中保持勵磁電流不變,勵磁電壓從冷態(tài)的11.89 V上升至熱態(tài)的13.47 V,上升明顯,勵磁電阻隨溫度的升高明顯變大,勵磁損耗增大。

        在轉(zhuǎn)速1 000 r/min、 勵磁電流6 A 和負載為2.875 Ω(電流10.5 A)時對DSEG 進行了負載測溫實驗,實驗初始環(huán)境溫度為21 ℃。溫升實驗后電機熱態(tài)的紅外溫度云圖如圖11 所示, 計算所得電機內(nèi)部勵磁繞組平均溫度為86.11 ℃,電樞繞組平均溫度為81.74 ℃。 冷態(tài)(實驗初始溫度)和熱態(tài)的三相DSEG 的實驗空載相電壓波形如圖12 所示。 繞組電阻隨溫度的變化并不影響電機空載相電壓波形,故空載相電壓波形主要反映鐵心材料性能隨溫度變化對電機性能的影響,對比圖12(a)和圖12(b),電機冷態(tài)和熱態(tài)空載相電壓波形變化不明顯,說明溫升對鐵心材料性能影響較小,進而對電機性能影響較小。

        圖10 空載勵磁繞組的電阻和溫度曲線Fig. 10 Resistance and temperature curves of excitation winding under no load

        圖11 熱平衡狀態(tài)紅外圖像Fig. 11 Infrared image in thermal equilibrium state

        對樣機熱狀態(tài)下的發(fā)電特性進行實驗測試,并建立了該溫度下的仿真模型。將仿真模型中鐵心材料設(shè)置為文獻[9]中100 ℃的B-H 曲線,繞組電阻設(shè)置為該實驗所得溫度下的阻值。 將冷態(tài)、熱態(tài)實驗結(jié)果和仿真結(jié)果進行了對比分析,其外特性及功率曲線如圖13 所示,仿真與實驗波形吻合度高,驗證了改變材料進行有限元建模的正確性。對比冷態(tài)實驗和20 ℃有限元模型仿真結(jié)果, 仿真外特性曲線略高于實驗曲線, 這是由于仿真模型更為理性化。熱態(tài)時的仿真外特性曲線略低于實驗曲線,這是因為仿真模型中鐵磁材料的設(shè)置為100 ℃,而溫升實驗中,鐵心溫度分布不均且低于100 ℃,仿真模型中鐵心材料溫度設(shè)置偏高導(dǎo)致磁導(dǎo)率略微下降。實驗和仿真結(jié)果均顯示熱態(tài)時的電機外特性下移,相同勵磁電流下的功率輸出減小,這是鐵心磁阻和電樞繞組電阻隨溫升變大共同作用的結(jié)果。 事實上,高溫下DSEG 若要保持輸出電壓恒定,需增大勵磁電流,進而使得勵磁損耗進一步增加,造成更大的勵磁負擔(dān),從而給高溫下DSEG 勵磁繞組的設(shè)計帶來嚴峻挑戰(zhàn)。

        圖12 實驗空載相電壓波形Fig. 12 Experimental waveforms of no-load phase voltage

        圖13 溫升實驗下外特性及功率對比Fig. 13 Comparison of external characteristics and power in temperature rise experiment

        4.2 基于直流電源靜止熱實驗的特性分析

        為了使得仿真分析與實驗驗證條件更為接近,采用直流電源靜止熱實驗的方式加熱樣機,將定子電樞繞組串聯(lián)連接并通入12 A 的直流電。 記錄直流電源所示的電壓電流, 隨時監(jiān)測繞組阻值變化,并計算繞組瞬時平均溫度,通過這種方式將電樞繞組加熱至100 ℃進行電機特性分析,并與100 ℃仿真模型進行對比分析。 圖14 為實驗和仿真的DSEG 外特性和功率曲線,可以看出仿真與實驗波形較為吻合,驗證了改變材料進行有限元建模的正確性。

        圖14 靜止熱實驗下外特性及功率曲線對比Fig. 14 Comparison of external characteristic and power curves in static thermal test

        5 結(jié)論

        本文通過磁路計算、仿真分析和實驗的方法對DSEG 不同溫度下的電機特性進行了研究,得出以下結(jié)論。

        (1)采用等效磁路法,并結(jié)合有限元法分析溫度對電機磁場和繞組電感的影響,溫度對各相磁鏈的影響最大的位置為定、轉(zhuǎn)子極正對位置。 溫度上升,電機各相電感最大值減小,最小值幾乎不變。

        (2)隨溫度上升,DSEG 外特性下移,輸出功率降低。

        (3)通過建模對比分析,繞組銅材性能隨溫度的變化對電機性能的影響大于鐵心材料的影響。

        (4)對一種三相12/8 極DSEG 進行了電機動態(tài)溫升實驗和直流電源靜止熱實驗的方式加熱樣機,并在此基礎(chǔ)上進行冷態(tài)和熱態(tài)的電機溫度特性實驗,驗證了分析和仿真的正確性。 研究環(huán)境溫度變化對DSEG 輸出特性的影響,為不同環(huán)境下工作的DSEG 設(shè)計和使用提供參考。

        猜你喜歡
        磁路磁阻鐵心
        基于分布磁路法的非晶合金高速磁浮直線電機懸浮力計算
        永磁磁阻電動機的研究
        防爆電機(2021年4期)2021-07-28 07:42:46
        異步電動機定子鐵心模態(tài)及振動響應(yīng)分析
        防爆電機(2021年3期)2021-07-21 08:12:54
        汽輪發(fā)電機階梯段鐵心損耗分析
        江西寧都:鐵心硬手守護綠水青山
        巨磁阻電渦流傳感器設(shè)計
        基于等效磁路法的永磁同步電機特性分析
        微特電機(2016年9期)2016-12-20 05:18:48
        交/直流線路并行運行對鐵心飽和不穩(wěn)定的影響
        四相開關(guān)磁阻電機的四電平DITC調(diào)速系統(tǒng)
        基于ANSYS Workbench的微型揚聲器磁路優(yōu)化分析
        拍摄av现场失控高潮数次| 日本高清乱码中文字幕| 欧美丰满熟妇bbbbbb| 色猫咪免费人成网站在线观看| 国产精品国产三级国产av创| 日本免费三片在线视频| 亚洲人成网站18禁止| 99久久精品费精品国产一区二区 | 日韩有码中文字幕在线视频| 久久精品国产字幕高潮| 亚洲乱码av中文一区二区| 亚洲精品成人av观看| 亚洲一区二区三区免费的视频 | 久久er99热精品一区二区| 无码电影在线观看一区二区三区| 亚洲一区二区女优视频| 亚洲最大成人综合网720p| 日韩人妻无码一区二区三区| 亚洲一区二区欧美色妞影院| 中文字幕人妻av四季| 中文字幕人妻熟女人妻| 亚洲欧美国产双大乳头| 少妇被粗大猛进进出出| 日本男人精品一区二区| 精品成人av一区二区三区| 国产精品一区二区资源| 一本到亚洲av日韩av在线天堂| 久久久久亚洲av成人片| 韩国19禁主播深夜福利视频| 国产一级片内射在线视频| 少妇人妻中文久久综合| 婷婷五月六月综合缴情| 999精品免费视频观看| 成人自拍偷拍视频在线观看| 亚洲av综合av成人小说| 国产精品一区二区在线观看99| 超短裙老师在线观看一区| 精品香蕉一区二区三区| 天码av无码一区二区三区四区 | 99久热在线精品视频观看| 中文字幕高清在线一区二区三区 |