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

改善單相混合勵磁磁通切換電機轉矩死區(qū)的不對稱轉子建模與測試

2015-11-15 09:18:14張宗盛王秀和楊玉波
電工技術學報 2015年6期

張宗盛 王秀和 楊玉波

(山東大學電氣工程學院 濟南 250061)

1 引言

磁通切換型磁阻電機(Flux Switching Motor,FSM)是在開關磁阻電機的基礎上發(fā)展而來的一種新型雙凸極電機[1],除了具有開關磁阻電機的一系列優(yōu)點外,還有更高的繞組利用率、更簡單的控制電路等優(yōu)點。按勵磁方式不同,分為電勵磁、永磁體勵磁及混合勵磁[2-5]。電勵磁 FSM由勵磁繞組提供勵磁磁動勢,通過調節(jié)勵磁電流可以方便地調節(jié)勵磁磁動勢,但存在勵磁電流較大、損耗較高的缺點;永磁式FSM利用永磁體提供勵磁磁動勢,不需要勵磁電流但存在氣隙磁場難以調節(jié)的問題[6,7];混合勵磁磁通切換電機(Hybrid Excitation Flux Switching Motor, HEFSM)由勵磁電流與永磁體共同提供勵磁磁動勢,所需勵磁電流較小且方便調節(jié)具有較高的研究價值。FSM 有多相及單相結構,多相FSM又以三相結構為主,通過其比較靈活的定轉子極數(shù)比的設計,可以實現(xiàn)很好的三相電流對稱性,研究多借助于交直軸數(shù)學模型等方法,分析及實驗研究表明:通過一定的控制策略,三相FSM具有良好的轉矩輸出能力和恒功率弱磁擴速能力,在交流驅動領域具有良好的應用前景[8-13]。

單相FSM的定轉子極數(shù)比恒為2:1,其結構較三相FSM更為簡單,單相也就意味著所需的控制電路最為簡單。尤其是采用混合勵磁結構時,單相HEFSM 較三相 HEFSM 電機疊片造型簡單易于裝配。但是由于受結構影響單相FSM存在定轉子極軸線重合位置,該位置下電機磁路磁阻達到極小值點,從而產(chǎn)生轉矩死區(qū)[14-16]。有不同學者就單相 FSM轉矩死區(qū)問題進行了一系列研究,歸納起來主要有兩種方法:采用輔助裝置與改善轉子疊片的形狀[17,18]。在定子或轉子上安裝輔助裝置能夠改善轉子位于轉矩死區(qū)位置時的電機起動性能,但卻會削弱電機運行時的穩(wěn)定性,同時額外的輔助裝置會增加電機制造成本增大加工難度。改善轉子疊片的形狀主要有采用不對稱氣隙、轉子疊片置孔等措施,其改善轉矩死區(qū)效果并不是很理想,且轉子疊片形狀比較復雜,加工難度提高。

本文研究了一種不對稱的轉子結構,通過該轉子可以削弱單相 HEFSM的轉矩脈動,同時還分析了帶有該轉子結構的簡化模型的建模方法;最后進行了有限元仿真分析與實驗測試。結果表明,該轉子結構能夠有效地改善單相HEFSM轉矩死區(qū)問題。

2 單相HEFSM工作原理與轉矩死區(qū)

2.1 單相HEFSM工作原理

本文所研究的單相HEFSM如圖1所示。

圖1 單相HEFSM結構示意圖Fig.1 Schematic diagram of single-phase HEFSM

線圈F1~F4串聯(lián)組成勵磁繞組F、線圈A1~A4串聯(lián)組成電樞繞組A;永磁體與勵磁繞組導體放置在同一定子槽內(nèi),永磁體沿切向充磁,勵磁電流為直流。永磁體等效成一個恒定的虛擬內(nèi)稟磁動勢Fc與一個磁阻Rm的串聯(lián)支路,勵磁電流產(chǎn)生磁動勢為Ff,在不同的轉子位置,永磁體、勵磁繞組作用于氣隙的磁動勢方向相同,兩者并聯(lián)組成勵磁磁動勢,如圖2所示。

圖2 混合勵磁磁路示意圖Fig.2 Hybrid excitation magnetic circuit

勵磁電流方向不變,則等效勵磁磁動勢Fef方向不變。通以電樞電流后,勵磁、電樞磁動勢的合成磁動勢位于其中導體電流相反的勵磁繞組槽與電樞繞組槽之間定子極軸線上,電樞電流方向隨轉子位置不同而轉換,電樞電流的每次換向都會使定子合成磁動勢在空間上切換 45°,令電機轉子在直流磁場和交流磁場的共同作用下產(chǎn)生旋轉運動。

2.2 單相HEFSM轉矩死區(qū)分析

單相電勵磁FSM有勵磁、電樞繞組的串聯(lián)與并聯(lián)兩種控制方式,而對于 HEFSM不宜采用串聯(lián)控制,因為不能實現(xiàn)對勵磁電流的單獨調節(jié)功能,所以本文中對 HEFSM采用并聯(lián)控制方式。作為磁阻類電機,可采用APC(angular position control)、CCC(current chopping control)兩種控制策略,本文采用導通區(qū)間為45°的CCC控制策略,該策略有利于調節(jié)勵磁電流的大小,同時最大程度上避免繞組電流關斷帶來的轉矩死區(qū)?;诖瞬呗?,利用Ansoft Simplorer V10與 Ansoft Maxwell 14 2D搭建了HEFSM系統(tǒng)的聯(lián)合仿真模型,其中HEFSM為圖1所示的常規(guī)轉子結構,圖3為不同電流斬波之下的轉矩仿真結果曲線簇。

圖3 常規(guī)轉子轉矩仿真曲線Fig.3 Simulation torque curves of the conventional rotor

從圖 3中可以看出每條轉矩曲線均存在過零點,轉矩過零點也就是在定轉子極軸線重合的位置,即使能夠理想地進行電樞電流的換向,該位置下電機的電磁轉矩為零。并且在電機運行時存在摩擦阻力等因素,所以單相HEFSM存在轉矩死區(qū)。

3 不對稱轉子結構設計與建模

3.1 不對稱轉子結構設計

從本質上看 HEFSM的運行原理遵循“磁阻最小原理”,即磁通總要沿著磁阻最小的路徑閉合,而具有一定形狀的鐵心在移動到最小磁阻位置時,必使自己的主軸線與磁場的軸線重合。HEFSM 轉子的主軸線為轉子極軸線,磁場的軸線沿勵磁、電樞合成磁動勢的方向即定子極的軸線。當HEFSM為8/4極結構時,存在圖1所示的轉子極軸線與定子極軸線重合位置,該位置即為轉矩死區(qū)的中心點。

HEFSM的轉子結構比定子簡單,本文同樣從設計轉子結構入手來解決轉矩死區(qū)問題。本文所研究的不對稱轉子結構分為疊片數(shù)量相等的 a、b兩部分,這兩部分的極弧長度不等,如圖4a、4b所示,這兩部分沿軸向裝配時如圖 4c所示,所以該轉子結構沿徑向為非對稱結構。

圖4 不對稱轉子結構Fig.4 Asymmetric rotor structure

通過前文可知,轉子的a、b兩部分均會產(chǎn)生轉矩死區(qū),其死區(qū)的中心位置如圖4a、4b所示,當兩者不對稱裝配后,其轉矩死區(qū)中心位置不會重疊,a、b兩部分疊片軸線的夾角Δθ越大越是有利于避免兩者的轉矩死區(qū)的重合。電樞電流的換向區(qū)間依然為45°,根據(jù)上述 a、b極弧長度的選取,轉子在 45°區(qū)間內(nèi),a、b極軸線不可能與定子極的軸線同時重合,通過恰當?shù)目刂齐姌须娏鲹Q向時刻,可達到改善轉矩死區(qū)目的。

若定義τs為定子極弧長度,τslot為槽口寬度。為盡可能的增大Δθ的數(shù)值轉子 a、b兩部分的極弧長度分別取為

本文所設計裝配有不對稱轉子的 HEFSM樣機主要尺寸參數(shù)見下表。

表 HEFSM樣機的主要尺寸參數(shù)Tab. The parameters of prototype HEFSM

3.2 帶有不對稱轉子的HEFSM簡化模型

由于轉子結構的不對稱,定轉子間的氣隙形狀更為復雜,建立電機的準確模型不易實現(xiàn),因此本文建立了 HEFSM 的簡化模型。在電機結構上 HEFSM的勵磁繞組槽軛部截面積很小,從圖2所示磁路可以看出,勵磁電流逐漸增大,該部分的磁通密度會逐漸減小至零后又逐漸增大,所以該部分的磁導變化范圍受勵磁電流的影響很大。其余鐵心部分磁通密度不會過高且變化范圍較小可忽略其磁導。單相HEFSM的簡化磁路模型如圖 5所示。圖中,Gy為勵磁繞組槽軛部磁導;Ff為勵磁繞組磁動勢;Fa為電樞繞組磁動勢;每對定子極可區(qū)分為 s1、s2兩個極,分別為勵磁槽逆時針側與順時針側,s1、s2與鄰近轉子極間的氣隙磁導分別表示為Go1、Go2。

圖5 單相HEFSM的簡化磁路模型Fig.5 Simplified circuit model of single-phase HEFSM

根據(jù)磁路結構的對稱性,以結點0為磁動勢參考點,電樞繞組磁通為φa,則有結點磁路方程為

通過迭代方法求解結點磁路方程,得F1、F2和φa,進而可得式(4)、式(5);勵磁、電樞繞組均為4線圈串聯(lián)而成所以可得式(6)~式(8)。

式中,ief、Lef為等效的勵磁繞組電流、自感;La為電樞繞組自感;Lfa為互感。

利用式(4)~式(8)可得電機電磁轉矩Te為

3.3 氣隙磁導的求解

不對稱轉子結構改變了 HEFSM磁通密度分布沿軸向均勻分布的特點,借助于 3D電磁場有限元仿真軟件可以對該類轉子進行仿真耗時耗力,對分析 HEFSM帶來不便,若建立簡潔的數(shù)學模型則會極大的提高其實用價值。

圖5中的Go1、Go2可認為分別包含兩部分且是以轉子位置角θ為變量的函數(shù),故可表示為

式中,Ga1(θ)、Ga2(θ)為轉子 a 部分與定子極 s1、s2間的氣隙磁導;Gb1(θ)、Gb2(θ)為轉子b部分與定子極s1、s2間的氣隙磁導。

相同定轉子極弧長度的 8/4極單相 HEFSM與8/4極的單相開關磁阻電機(SRM)在定轉子極間氣隙結構上是一致的。而單相 SRM 不存在繞組耦合的問題,所以利用單相 SRM 電感有限元仿真結果及電感與回路磁導間的關系可快速的確定某一定子極與轉子極之間的氣隙磁導,在忽略鐵心磁導的基礎上,電感與氣隙磁導的關系為[19]

若以圖 3a、3b所示為轉子位置角參考零點,Ga1(θ)、Gb1(θ)用上述方法確定的曲線如圖6所示。

圖6 氣隙磁導有限元仿真曲線Fig.6 Air gap permeance curve of finite element analysis

若定義a部分軸線與定子軸線重合的位置為參考零點,利用傅里葉級數(shù)擬合可獲得Ga1(θ)、Gb1(θ)的表達式。由其曲線可知為偶函數(shù),所以只用傅里葉級數(shù)中的余弦分量擬合為

因此 s2與定子極間氣隙磁導Ga2(θ)、Gb2(θ)為

根據(jù)圖 3c中轉子的裝配方式,a、b兩部分軸線的夾角為

所以有

至此便可求得式(10)、式(11)氣隙磁導的數(shù)學表達式,代入式(3)后,最終將能夠給實現(xiàn)式(9)的求解,實現(xiàn)對電機轉矩的求解。

4 仿真與測試分析

4.1 有限元仿真分析

裝有不對稱轉子的 HEFSM的換向時刻應在轉子a、b兩部分軸線夾角的某條角分線與定子極軸線的重合時刻。通過表1數(shù)據(jù)可求得式(13)所示的Δθ約為12°,也就是說若以轉子a部分軸線與定子軸線重合的位置為參考點,則不對稱轉子 HEFSM的電樞電流換向時刻可較對稱轉子提前0~12°。在進行有限元仿真時設定為提前 10°。圖 7為給定勵磁電流if=3A、電樞電流ia=|4|A條件下,利用本文提出的數(shù)學模型計算與 3D有限元分析得到的不同轉子位置電磁轉矩結果。從圖中可以看出,在該激勵電流下不同轉子位置角均可產(chǎn)生電磁轉矩,存在改善轉矩死區(qū)的效果。

圖7 給定電流下電磁轉矩分析結果對比Fig.7 Torque value comparison between the model calculation and the 3D finite element analysis

4.2 聯(lián)合仿真分析

為考察控制電路對轉矩輸出的影響,利用Ansoft Simplorer V10搭建了并聯(lián)控制仿真外電路,與利用Ansoft Maxwell 14 3D搭建的電機本體仿真模型構成了聯(lián)合仿真模型。通過設置仿真電路中的開關器件等仿真元件參數(shù),能夠很好地模擬實際元件的特性,所以仿真電路會表現(xiàn)出受到開關器件延時等因素的影響,因此將以超前轉子a部分軸線與定子軸線重合位置 12°的位置為分析不對稱轉子HEFSM的參考零點。

直流電源電壓為80V;因勵磁電流的方向不變所以采用電流滯環(huán)斬波控制,滯環(huán)寬度為0.2A;而電樞電流要隨轉子位置而換向,為提高其響應速度未進行滯環(huán)設置;仿真步長為1×10-4s。圖8為控制if=3A、ia=4A,轉速為750r/min時的聯(lián)合仿真結果。

圖8 聯(lián)合仿真結果曲線Fig.8 The results of combined simulation model

由圖8b與圖7的對比可以,在外電路的控制電機連續(xù)運行過程中,轉矩的最低點約為0.25N·m較圖7中的轉矩最小值有所降低,但與圖3所示的常規(guī)轉子HEFSM的轉矩曲線相比消除了轉矩過零點,說明在電機運行過程中改善轉矩死區(qū)效果明顯。

4.3 感應電動勢測試分析

考慮到實驗的可行性,進行了電樞繞組的感應電動勢的測試,電機給定轉速為 1 000r/min,勵磁電流為零,即測試永磁體對電樞繞組的影響。圖 9為電樞繞組感應電動勢的簡化模型計算值與實驗測試波形,其中,圖9a中參考零點與聯(lián)合仿真模型中一致,從兩者的對比可以看出,本文提出的簡化模型的準確度比較高。從圖9可以推出不對稱轉子結構并沒有明顯改變電機互感導數(shù)極性不變的區(qū)間范圍,說明電機的出力與繞組的互感、自感均有關系,這改變了常規(guī)轉子 FSM 繞組自感總體上不輸出功率的情況,有利于改善電機轉矩死區(qū)。

圖9 電樞繞組感應電動勢波形Fig.9 The EMF waveforms of armature winding

4.4 轉矩測試分析

轉矩測試平臺的組成如圖 10所示。動態(tài)轉矩傳感器采用 Lonrenz Messtechnik GmbH的產(chǎn)品 DR-3000,最大偏差為 50mN·m,采樣頻率為 2 500次/s[20]。通過磁滯測功機為電機施加負載,HEFSM采用并聯(lián)控制外電路,控制電路設置與聯(lián)合模型中的仿真電路設置一致。

圖10 轉矩測試平臺示意圖Fig.10 Schematic diagram of torque test platform

測試時給定負載0.5N·m,考慮到起動時刻的摩擦阻力較大等因素,在起動時的電流控制值較大為if=ia=5A。轉矩波形的初始位置即為過零點。圖11為實驗測試結果,從圖11可以看出,電流實驗波形與仿真波形基本一致,而轉矩實驗波形與圖8b的轉矩仿真波形差別較大,實驗波形甚至出現(xiàn)了負值,其主要原因在于:轉矩的仿真波形為電機本身的理論上的轉矩輸出,不受負載的影響。轉矩傳感器測量的是兩個軸端的轉矩差值,受電機與負載之間作用力的影響。

圖11 實驗測試結果Fig.11 Experimental test waveforms

負載設置為0.5N·m恒定值,由圖8b可知即使在一個轉子步距內(nèi)電機的電磁轉矩也存在大于或小于負載轉矩的不同情況,這使得電機即使在穩(wěn)定運行時也是始終處于不斷的加、減速狀態(tài),轉速非恒定值而是在一定范圍內(nèi)波動。電機的轉動慣量要明顯小于負載的轉動慣量,電機、負載之間的相互作用力變得更為復雜。在電機轉速下降時,負載減速反應較電機有一定的滯后,此時轉矩傳感器測量值較實際值偏小,同理電機轉速上升時測量值較實際值偏大,因此轉矩的實驗波形較電磁轉矩仿真分析波形得波動更為嚴重,乃至出現(xiàn)反向轉矩值。電機在起動的瞬間電機為加速狀態(tài)負載更接近理想狀態(tài),所以初始波形較為理想。

如在圖11b中,S點為電樞電流關斷并開始換向的時刻,t1區(qū)域內(nèi)的波形出現(xiàn)負值明顯與理論曲線差異很大,出現(xiàn)較大的反向轉矩波形尖峰,而圖11a中的電流波形并不存在突變,所以出現(xiàn)很大的反向轉矩是明顯有悖于能量守恒的,其成因是該區(qū)域內(nèi)出現(xiàn)大轉動慣量的負載帶動小轉動慣量的電機正向旋轉的情況,以至傳感器兩個連接軸之間的轉矩差值即變?yōu)樨撝怠.旊姌须娏鲹Q向完畢,逐漸轉換為電機拖動負載旋轉的狀態(tài),在此過程中出現(xiàn)t2區(qū)域所示的轉矩曲線抖動情況。電機拖動負載逐漸的穩(wěn)定的過程中,電機與負載之間的相互作用力逐漸變的更為復雜:時而電機拖動負載旋轉,時而因轉動慣量的原因使得負載帶動電機旋轉,從而轉矩傳感器的輸出波形表現(xiàn)為上下波動,與理論波形的差異十分明顯。

從前文分析可知,在電樞電流換向的轉子位置點轉矩輸出能力最弱,因此對該轉子位置進行轉矩測試,看電機能否起動可驗證不同轉子結構的改善轉矩死區(qū)的效果,如圖12所示。不對稱轉子結構使得電機在電樞電流換向位置點能夠在輕微抖動后起動,說明在該位置不存在轉矩死區(qū),而對稱轉子結構則在抖動后轉矩歸零出現(xiàn)轉矩死區(qū),無法起動。

圖12 不同轉子結構的電樞電流換向點轉矩波形Fig.12 Torque test waveforms of different rotor structures at armature current reversal position

5 結論

本文研究了一種適用于改善單相 HEFSM轉矩死區(qū)的不對稱轉子結構,通過介紹該轉子的結構特點說明了其改善單相 HEFSM轉矩死區(qū)的原理,同時還介紹了一種針對該轉子的簡化模型的建模方法,最后進行了有限元仿真與聯(lián)合仿真分析與實驗測試,分析了實測結果與理論上差異的成因,并進行了特殊轉子位置的實驗對比,其結果表明本文所研究的不對稱轉子結構可有效的改善單相 HEFSM的轉矩死區(qū)問題,所建簡化模型的準確性也較高。

[1] Pollock C, Wallace M. The flux switching motor, a DC motor without magnets or brushes[C]. Proceedings of 34th Annual Meeting of the IEEE Industry Applications, Phoenix, 1999, 3: 1980-1987.

[2] Hoang E, Ahmed A H Ben, Lucidaume J. Switching flux permanent magnet polyphased machines[C].Proceeding of European Conference of Power Electronic and Applications Trondheim, Norway, 1997(3): 903-908.

[3] Zhu Z Q, Pang Y, Howe D, et al. Analysis of electromagnetic performance of flux-switching permanent magnet machines by nonlinear adaptive lumped permanent magnet circuit it model[J]. IEEE Transactions on Magnetics, 2005, 41(11): 4277-4287.

[4] Hoang E. Flux-switching dual-excitation electrical machine: EP, 20070726708[P]. 2008-12-17.

[5] 花為, 程明, 諸自強, 等. 新型兩相磁通切換型雙凸極永磁電機的靜態(tài)特性研究[J]. 電工技術學報,2006, 21(6): 70-77.

Hua Wei, Cheng Ming, Zhu Z Q, et al. Study on static characteristics of a novel two-phase flux-switching doubly-salient permanent magnet machine[J]. Transactions of China Electrotechnical Society, 2006, 21(6):70-77.

[6] Chai K S, Pollock C. Using genetic algorithms in design optimization of the flux switching motor[C].Conference of Power Electronics, Machines and Drives,2002: 540-545.

[7] Chen Y, Chen S, Zhu Z Q, et al. Starting torque of single-phase flux-switching permanent magnet motors[J]. IEEE Transactions on Magnetics, 2006, 42(10):3416-3418.

[8] 花為, 程明. 新型三相磁通切換型雙凸極永磁電機電感特性分析[J]. 電工技術學報, 2007, 22(11):21-28.

Hua Wei, Cheng Ming. Inductance characteristics of 3-phase flux-switching permanent magnet machine with doubly-salient structure[J]. Transactions of China Electrotechnical Society, 2007, 21(11): 21-28.

[9] Owen R L, Zhu Z Q, Jewell G W. Hybrid excited flux-switching permanent-magnet machines with iron flux bridges[J]. IEEE Transactions on Magnetics,2010, 46(6): 1726-1729.

[10] 許澤剛, 謝少軍, 毛鵬. 混合勵磁磁通切換電機等效磁路模型[J]. 中國電機工程學報, 2011, 31(30):124-131.

Xu Zegang, Xie Shaojun, Mao Peng. Equivalent magnetic circuit model of flux-switching hybrid excitation machine[J]. Proceedings of the CSEE, 2011,31(30): 124-131.

[11] 王宇, 鄧智泉, 王曉琳. 基于空間矢量調制的電勵磁磁通切換電機轉矩角增量控制[J]. 中國電機工程學報, 2012, 32(3): 75-82.

Wang Yu, Deng Zhiquan, Wang Xiaolin. A torque angel increment control strategy for electrical excitation flux-switching machines based on space vector modulation[J]. Proceedings of the CSEE, 2012, 32(3): 75-82.

[12] 賈紅云, 程明, 花為, 等. 基于死區(qū)補償?shù)拇磐ㄇ袚Q永磁電機定子磁場定向控制[J]. 電工技術學報,2010, 25(11): 48-55.

Jia Hongyun, Cheng Ming, Hua Wei, et al. Statorflux-oriented control for flux-switching permanent magnet motor based on dead-time compensation[J].Transactions of China Electrotechnical Society, 2010,25(11): 48-55.

[13] 劉細平, 左亮平, 刁艷美, 等. 軸向磁通切換混合勵磁電機勵磁控制系統(tǒng)[J]. 電工技術學報, 2013,28(10): 97-102.

Liu Xiping, Zuo Liangping, Diao Yanmei, et al.Excitation control system of axial flux-switching hybrid excited machine[J]. Transactions of China Electrotechnical Society, 2013, 28(10): 97-102.

[14] Bangura J F. Design of high-power density and relatively high efficiency flux-switching motor[J]. IEEE Transactions on Energy Conversion, 2006, 21(2): 416-425.

[15] 張宗盛, 王秀和, 楊玉波. 新型混合勵磁磁通切換型磁阻電機的變結構磁路模型[J]. 中國電機工程學報, 2014, 34(21): 3467-3474.

Zhang Zongsheng, Wang Xiuhe, Yang Yubo. A variable structure magnetic circuit model of novel hybrid excitation flux switching motors[J]. Proceedings of the CSEE, 2014, 34(21): 3467-3474.

[16] 張宗盛, 王秀和, 楊玉波, 等. 單相磁通切換型磁阻電機的快速建模方法[J]. 電機與控制學報, 2014,18(7): 6-11.

Zhang Zongsheng, Wang Xiuhe, Yang Yubo, et al.Speedy Modeling method for single-phase flux switching motor[J]. Electric Machines and Control, 2014,18(7): 6-11.

[17] 吳建華. 開關磁阻電機設計與應用[M]. 北京: 機械工業(yè)出版社, 2000.

[18] Chen Y, Chen S, Zhu Z Q, et al. Starting torque of single-phase flux-switching permanent magnet motors[J]. IEEE Transactions on Magnetics, 2006, 42(10):3416-3418.

[19] 王秀和. 電機學[M]. 北京: 機械工業(yè)出版社, 2009.

[20] Lorenz Messtechnik Gmbh. Operation manual for configuration and evaluation software DR-USB-VS[EB/OL]. www.lorenz-sensors.com.

主站蜘蛛池模板: 伦精品一区二区三区视频| 99视频国产精品| 亚洲永久色| 亚洲最大福利网站| h视频在线播放| 亚洲精品午夜天堂网页| 国产91丝袜| 亚洲综合色婷婷中文字幕| 日韩一区二区在线电影| 国产国语一级毛片| 女人爽到高潮免费视频大全| 好紧好深好大乳无码中文字幕| 婷婷六月激情综合一区| 国产免费久久精品99re丫丫一| 国产精品视频第一专区| 国产原创第一页在线观看| 中文字幕在线观看日本| 亚洲IV视频免费在线光看| 日本不卡视频在线| 国产欧美日韩在线一区| 久久精品无码中文字幕| 日本高清有码人妻| 美女国产在线| 麻豆国产精品视频| 久青草国产高清在线视频| 中文字幕无线码一区| 一区二区三区精品视频在线观看| 国产欧美日韩免费| 欧美午夜视频| 欧洲一区二区三区无码| 亚洲精品爱草草视频在线| 91欧美在线| 香蕉精品在线| 97在线观看视频免费| 中文字幕人妻av一区二区| 国产丰满成熟女性性满足视频 | 亚洲αv毛片| 一级毛片免费的| 色妞永久免费视频| 欧美精品成人一区二区视频一| 中文字幕在线永久在线视频2020| 99久久99视频| 国产成a人片在线播放| 日本三区视频| 激情午夜婷婷| 人妻丰满熟妇AV无码区| 国产成熟女人性满足视频| 国产亚洲精品精品精品| 亚洲综合在线网| 四虎国产精品永久一区| 超碰aⅴ人人做人人爽欧美 | 呦女精品网站| 亚洲国产欧美国产综合久久| 内射人妻无码色AV天堂| 亚洲综合经典在线一区二区| 99热这里都是国产精品| 国产精品一区二区在线播放| 国产亚洲欧美另类一区二区| 国产精品999在线| 欧美日本在线| 色精品视频| 99草精品视频| 特级做a爰片毛片免费69| 欧美成在线视频| 欧美成人精品欧美一级乱黄| 国产视频一区二区在线观看| 麻豆国产精品| 免费国产好深啊好涨好硬视频| 欧美国产综合色视频| 亚洲欧美成aⅴ人在线观看| 国产亚洲精品精品精品| 国产人碰人摸人爱免费视频 | 国产永久无码观看在线| 亚洲第一在线播放| 99久久婷婷国产综合精| 亚洲无码精品在线播放| 国产一线在线| 欧美国产日本高清不卡| 亚洲精品图区| 色网站免费在线观看| 亚洲精品图区| 亚洲精品久综合蜜|