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

    基于田口法和遺傳算法的取向硅鋼片直線電機(jī)優(yōu)化設(shè)計(jì)

    2023-12-28 09:14:26高子晏
    微特電機(jī) 2023年11期
    關(guān)鍵詞:硅鋼片永磁體遺傳算法

    高子晏,董 婷

    (沈陽工業(yè)大學(xué) 電氣工程學(xué)院,沈陽 110870)

    0 引 言

    在工業(yè)機(jī)器人、軍工和航空航天等領(lǐng)域,對電機(jī)的體積、推力密度和定位精度等指標(biāo)要求十分苛刻,希望電機(jī)做到輕量化的同時(shí),具有短時(shí)高過載能力[1-3]。取向硅鋼片的磁性能與軋制方向密切相關(guān),多應(yīng)用于變壓器鐵心且只應(yīng)用高磁導(dǎo)率的軋制方向[4]。文獻(xiàn)[5]采用取向硅鋼片制作永磁同步電機(jī)定子,利用其軋制方向高磁導(dǎo)率的特點(diǎn),緩解了電機(jī)定子齒部的飽和程度,提高了電機(jī)的功率密度。文獻(xiàn)[6]利用取向硅鋼片構(gòu)建了拼接齒結(jié)構(gòu)的永磁力矩電機(jī),提高了電機(jī)的轉(zhuǎn)矩密度和過載能力。文獻(xiàn)[7]利用取向硅鋼片降低了電機(jī)的鐵耗,提高了電機(jī)的效率和輸出轉(zhuǎn)矩。旋轉(zhuǎn)電機(jī)應(yīng)用取向硅鋼片多采用拼接或榫卯結(jié)構(gòu),這會提高電機(jī)的制造成本和制作難度。

    近年來,電機(jī)多目標(biāo)優(yōu)化逐漸向代理模型結(jié)合優(yōu)化算法方向發(fā)展。文獻(xiàn)[8]采用Pearson相關(guān)系數(shù)對各優(yōu)化變量進(jìn)行分層,基于Kriging模型和NSGA II對各層參數(shù)進(jìn)行多目標(biāo)優(yōu)化,降低了電機(jī)的轉(zhuǎn)矩脈動和損耗。文獻(xiàn)[9]基于田口法對組合磁極永磁同步電機(jī)的轉(zhuǎn)矩與懸浮力進(jìn)行快速尋優(yōu)。文獻(xiàn)[10]結(jié)合響應(yīng)面法和遺傳算法對圓筒型永磁直線電機(jī)進(jìn)行多目標(biāo)優(yōu)化,提高了電機(jī)的綜合性能指標(biāo)。田口法可以提高尋優(yōu)效率但全局尋優(yōu)能力較差,遺傳算法全局優(yōu)化能力強(qiáng)但存在早熟收斂導(dǎo)致無法尋得全局最優(yōu)的問題。結(jié)合田口法和遺傳算法可以充分發(fā)揮兩者的優(yōu)勢,具有更好的多目標(biāo)優(yōu)化效果。

    本文介紹了取向硅鋼片直線電機(jī)的結(jié)構(gòu)特征,構(gòu)建了該電機(jī)的磁網(wǎng)絡(luò)模型,結(jié)合該模型以電磁推力、推力波動和永磁體利用率為優(yōu)化目標(biāo),研究了一種基于靈敏度分析的分層優(yōu)化方法。針對顯著優(yōu)化變量,結(jié)合田口法與遺傳算法進(jìn)行多目標(biāo)優(yōu)化,普通優(yōu)化變量則采取單參數(shù)掃描法,針對長行程和高推力密度兩種應(yīng)用場合分別給出了兩種最優(yōu)方案。

    1 取向硅鋼片直線電機(jī)磁網(wǎng)絡(luò)模型

    1.1 取向硅鋼片直線電機(jī)

    以一臺傳統(tǒng)永磁直線同步電機(jī)為例,將其初級鐵心替換為取向硅鋼片,作為取向硅鋼片直線電機(jī)的初始方案,圖1為30P120型取向硅鋼片和DW465-50型無取向硅鋼片的磁導(dǎo)率-磁通密度曲線,表1為直線電機(jī)初始方案的結(jié)構(gòu)參數(shù)。

    表1 12槽14極直線電機(jī)初始方案結(jié)構(gòu)參數(shù)

    圖1 磁導(dǎo)率-磁通密度曲線

    將取向硅鋼片軋制方向與直線電機(jī)齒部主磁通方向保持一致,其軋制方向磁導(dǎo)率更高,剪切方向與軛部磁路方向相對應(yīng),取向硅鋼片直線電機(jī)拓?fù)浣Y(jié)構(gòu)如圖2所示。近極槽直線電機(jī)的軛部磁路較短,即使取向硅鋼片剪切方向的導(dǎo)磁能力不佳,對直線電機(jī)性能的影響仍然很小。本文直線電機(jī)額定負(fù)載電流為4.5 A,圖3為兩種材料直線電機(jī)推力-電流曲線。取向硅鋼片直線電機(jī)相比傳統(tǒng)直線電機(jī),在額定負(fù)載和3倍過載情況下推力分別提升4.03%和8.13%。

    圖2 取向硅鋼片直線電機(jī)拓?fù)浣Y(jié)構(gòu)

    圖3 推力-電流曲線

    1.2 磁網(wǎng)絡(luò)模型

    結(jié)合取向硅鋼片直線電機(jī)的拓?fù)浣Y(jié)構(gòu)和磁通路徑構(gòu)建磁網(wǎng)絡(luò)模型,其電樞和永磁體磁動勢:

    (1)

    式中:N為電樞繞組匝數(shù);i為電樞電流;hpm為永磁體高度;HC為永磁體矯頑力。

    應(yīng)用積分計(jì)算各支路磁阻,直線電機(jī)磁網(wǎng)絡(luò)模型中磁阻多為矩形,其計(jì)算公式如下:

    (2)

    式中:hcom為元件寬度;μ0和μr為真空磁導(dǎo)率和相對磁導(dǎo)率;wcom為元件寬度;lef為鐵心有效長度。

    直線電機(jī)初級齒軛連接處各向異性磁阻:

    (3)

    式中:μt和μy為軋制方向和剪切方向磁導(dǎo)率;θ為與軋制方向的夾角。

    本文構(gòu)建的取向硅鋼片直線電機(jī)磁網(wǎng)絡(luò)模型如圖4所示。有限元法計(jì)算耗時(shí)30 min,磁網(wǎng)絡(luò)法用時(shí)45 s。有限元法和磁網(wǎng)絡(luò)法求解的A相空載與負(fù)載反電動勢和推力對比如圖5和圖6所示。兩種方法的反電動勢計(jì)算結(jié)果基本吻合,兩種方法的推力計(jì)算結(jié)果為418.61 N和425.86 N,偏差為1.73%。

    圖4 取向硅鋼片直線電機(jī)磁網(wǎng)絡(luò)模型

    圖5 A相空載和負(fù)載反電動勢對比

    圖6 額定負(fù)載推力

    基于磁網(wǎng)絡(luò)模型求解直線電機(jī)的電磁特性相較于有限元法偏差在5%以內(nèi),驗(yàn)證了所構(gòu)建模型的正確性,為后續(xù)直線電機(jī)的深度優(yōu)化奠定了基礎(chǔ)。

    2 電機(jī)結(jié)構(gòu)參數(shù)的靈敏度分析

    直線電機(jī)的多目標(biāo)優(yōu)化屬于高維多物理場耦合問題,某些優(yōu)化參數(shù)與優(yōu)化目標(biāo)之間存在交叉耦合,部分參數(shù)則較為獨(dú)立且影響較小。如果對所有優(yōu)化參數(shù)均采用優(yōu)化算法,會使工作量呈指數(shù)增加,降低直線電機(jī)的優(yōu)化效率。因此,需要根據(jù)靈敏度分析對優(yōu)化參數(shù)進(jìn)行分層,針對綜合靈敏度大于0.3的顯著優(yōu)化變量,采用田口法結(jié)合遺傳算法進(jìn)行多目標(biāo)尋優(yōu),普通優(yōu)化變量則采用單參數(shù)掃描法優(yōu)化,進(jìn)而獲得取向硅鋼片直線電機(jī)的最優(yōu)設(shè)計(jì)方案。直線電機(jī)的多目標(biāo)優(yōu)化流程如圖7所示。基于磁網(wǎng)絡(luò)模型確定優(yōu)化參數(shù)的變化范圍并對其進(jìn)行靈敏度分析,后續(xù)采用有限元法進(jìn)行深度優(yōu)化。

    圖7 磁網(wǎng)絡(luò)法聯(lián)合有限元法優(yōu)化流程

    選出直線電機(jī)合適的初始優(yōu)化參數(shù),通過磁網(wǎng)絡(luò)模型確定各優(yōu)化參數(shù)合理的變化范圍,保持電流密度不變,J≈5 A/mm2,在槽面積變化時(shí)改變槽中繞組匝數(shù),保證槽滿率相同,sf≈0.7,電負(fù)荷與銅耗隨每槽導(dǎo)體數(shù)變化。表2給出了8個(gè)初始優(yōu)化參數(shù)及其變化范圍。

    表2 初始優(yōu)化參數(shù)和變化范圍

    初級鐵心長度的改變體現(xiàn)在端部槽寬度的變化,初級鐵心長度發(fā)生變化時(shí)端部齒寬度保持不變。端部齒寬度發(fā)生變化時(shí)初級鐵心長度保持不變,改變端部槽寬度。選擇取向硅鋼片直線電機(jī)的電磁推力、推力波動和永磁體利用率作為優(yōu)化目標(biāo),其目標(biāo)函數(shù):

    (4)

    式中:F′ave、F′rip和F′den表示初始模型取向硅鋼片直線電機(jī)的電磁推力、推力波動和永磁體利用率;σ1、σ2、σ3為權(quán)重系數(shù),滿足σ1+σ2+σ3=1;Fave(n)、Frip(n)和Fden(n)表示當(dāng)優(yōu)化參數(shù)為n時(shí),所對應(yīng)的優(yōu)化目標(biāo)數(shù)值。

    取向硅鋼片直線電機(jī)更適用于高推力密度場合,因此電磁推力是考慮的重點(diǎn);降低推力波動可以提高直線伺服系統(tǒng)的控制性能;隨著行程的增加,永磁體用量較大會帶來成本的激增,應(yīng)合理利用永磁材料。因此,電磁推力和推力波動的權(quán)重系數(shù)σ1和σ2別取0.5和0.3,永磁體利用率的權(quán)重系數(shù)σ3為0.2。

    永磁體利用率為電磁推力與單塊永磁體截面的比值:

    (5)

    利用磁網(wǎng)絡(luò)模型計(jì)算取向硅鋼片直線電機(jī)的電磁性能,對優(yōu)化參數(shù)進(jìn)行靈敏度分析:

    (6)

    式中:Fi(n)為當(dāng)優(yōu)化參數(shù)取值為n時(shí)優(yōu)化目標(biāo)i的數(shù)值;Fi為優(yōu)化目標(biāo)i的初始值。

    引入權(quán)重系數(shù)得到優(yōu)化參數(shù)的綜合靈敏度:

    S(n)=σ1|save(n)|+σ2|srip(n)|+σ3|sden(n)|

    (7)

    式中:save、srip和sden分別為優(yōu)化參數(shù)對應(yīng)優(yōu)化目標(biāo)的靈敏度。

    圖8和表3為8個(gè)初始優(yōu)化參數(shù)對應(yīng)優(yōu)化目標(biāo)的靈敏度和綜合靈敏度。初級齒部寬度wt、永磁體寬度wpm、永磁體高度hpm和初級鐵心長度lpri的綜合靈敏度較高,均大于0.3;齒部寬度和永磁體參數(shù)對電磁推力和永磁體利用率的影響較大,wt、hpm和lpri對推力波動有較大的影響。

    表3 優(yōu)化參數(shù)的靈敏度分析結(jié)果

    圖8 優(yōu)化參數(shù)的靈敏度分析

    綜上所述,選取wt、wpm、hpm和lpri為顯著優(yōu)化變量,采用優(yōu)化算法進(jìn)行深度優(yōu)化;軛部高度hy、槽口寬度wst、次級鐵心高度hsec和端部齒寬度wend作為普通優(yōu)化變量,僅采用單參數(shù)掃描法進(jìn)行優(yōu)化。

    各顯著優(yōu)化變量對應(yīng)目標(biāo)函數(shù)的函數(shù)值如圖9所示。函數(shù)值大于1說明該結(jié)構(gòu)參數(shù)下的電機(jī)綜合性能更優(yōu)異。為了盡可能提高取向硅鋼片直線電機(jī)的性能,并提升多目標(biāo)優(yōu)化效率,根據(jù)圖9針對性選取優(yōu)化變量的取值范圍,其中初級齒部寬度wt為5.2~5.5 mm,初級鐵心長度lpri為186~189.6 mm,永磁體寬度wpm為10~10.9 mm,永磁體高度hpm為3.1~4 mm,利于田口法尋取最優(yōu)參數(shù)組合。

    圖9 顯著優(yōu)化變量與目標(biāo)函數(shù)關(guān)系

    3 基于田口法的優(yōu)化設(shè)計(jì)

    基于上一節(jié)對各優(yōu)化參數(shù)的靈敏度分析和對優(yōu)化變量取值范圍的優(yōu)化,將顯著優(yōu)化變量作為田口正交試驗(yàn)的因數(shù),各顯著優(yōu)化變量及其水平數(shù)如表4所示。根據(jù)表5,采用田口正交矩陣僅需構(gòu)建16組試驗(yàn)方案的仿真模型?;谟邢拊ㄓ?jì)算各組試驗(yàn)的優(yōu)化目標(biāo)數(shù)值,根據(jù)優(yōu)化變量的平均值、方差和比重選取田口法的最優(yōu)方案。圖10為各優(yōu)化變量對優(yōu)化目標(biāo)的影響趨勢。

    表4 優(yōu)化變量及水平數(shù)

    表5 田口正交試驗(yàn)矩陣L16(44)

    圖10 優(yōu)化變量對性能指標(biāo)的影響

    方差可以反映一組數(shù)據(jù)各樣本點(diǎn)偏離總平均值的程度,優(yōu)化變量方差的計(jì)算公式:

    (8)

    式中:Fia為優(yōu)化目標(biāo)i的總平均值。

    取向硅鋼片直線電機(jī)各優(yōu)化變量對應(yīng)優(yōu)化目標(biāo)的比重如圖11和表6所示。

    表6 優(yōu)化變量的方差和比重

    圖11 優(yōu)化變量對應(yīng)優(yōu)化目標(biāo)的比重

    根據(jù)方差分析和比重分析,wt對電磁推力的影響最大,lpri對推力波動的影響最大,wpm對電磁推力的影響最大,hpm對永磁體利用率的影響最大。選擇優(yōu)化變量比重最大的優(yōu)化目標(biāo)對應(yīng)的最優(yōu)水平參數(shù),田口法的最優(yōu)方案為wt(1) ,lpri(4),wpm(4),hpm(1),初始方案與田口法結(jié)構(gòu)參數(shù)對比如表7所示,優(yōu)化后電機(jī)的初級鐵心長度略微增加。

    表7 田口法優(yōu)化方案結(jié)構(gòu)參數(shù)對比

    由圖10可知,wpm和hpm對電磁推力和永磁體利用率兩個(gè)優(yōu)化目標(biāo)的作用趨勢截然相反。田口法的優(yōu)化原理為將多目標(biāo)優(yōu)化簡化為多個(gè)單目標(biāo)優(yōu)化,基于該方法對取向硅鋼片直線電機(jī)進(jìn)行多目標(biāo)優(yōu)化存在弊端,無法全面考慮各優(yōu)化目標(biāo)之間的耦合矛盾關(guān)系。

    4 基于遺傳算法的優(yōu)化設(shè)計(jì)

    取向硅鋼片直線電機(jī)的優(yōu)化目標(biāo)之間存在拮抗關(guān)系,采用田口法優(yōu)化電機(jī)并沒有實(shí)現(xiàn)多優(yōu)化目標(biāo)的解耦。遺傳算法具有高效的全局搜索能力,可以解決取向硅鋼片直線電機(jī)多目標(biāo)優(yōu)化中各優(yōu)化目標(biāo)相互制約的問題。將田口法各方案的計(jì)算結(jié)果作為遺傳算法的初始種群,替代隨機(jī)生成的初始種群,可以為遺傳算法的全局優(yōu)化提供參考。

    基于上一節(jié)田口法的優(yōu)化結(jié)果和各優(yōu)化變量的比重分析,針對性選取顯著優(yōu)化變量的變化范圍,可以減少優(yōu)化變量的樣本空間,提高優(yōu)化效率,使優(yōu)化算法快速向最優(yōu)參數(shù)組合收斂。遺傳算法參數(shù)取值如表8所示,以取向硅鋼片直線電機(jī)初始方案為基準(zhǔn),對優(yōu)化目標(biāo)進(jìn)行約束:

    (9)

    表8 遺傳算法參數(shù)設(shè)置

    基于遺傳算法的取向硅鋼片直線電機(jī)優(yōu)化結(jié)果如圖12所示。由圖12可知,電磁推力與永磁體利用率近乎呈反比關(guān)系,針對不同工程實(shí)際需求,在Pareto最優(yōu)解集中選取方案,本質(zhì)為在不同優(yōu)化目標(biāo)之間折衷。由于遺傳算法優(yōu)化變量的參數(shù)值為隨機(jī)生成,考慮到實(shí)際加工的可行性,構(gòu)建尺寸數(shù)據(jù)僅保留小數(shù)點(diǎn)后一位的有限元模型,新有限元仿真結(jié)果與遺傳算法理論最優(yōu)方案可能存在差距。

    圖12 取向硅鋼片直線電機(jī)三維Pareto前沿

    以式(4)目標(biāo)函數(shù)值最大為最終優(yōu)化目標(biāo),在Pareto前沿中選取優(yōu)化方案1,為了方便對比分析遺傳算法的優(yōu)化效果,在Pareto前沿中選取永磁體結(jié)構(gòu)參數(shù)與初始方案相似的優(yōu)化方案2作為對照組,兩方案的結(jié)構(gòu)參數(shù)如表9所示。優(yōu)化方案1的結(jié)構(gòu)參數(shù)與田口法優(yōu)化結(jié)果相似,進(jìn)一步驗(yàn)證了田口法優(yōu)化的有效性。

    表9 遺傳算法優(yōu)化方案結(jié)構(gòu)參數(shù)

    在基于田口法和遺傳算法對顯著優(yōu)化變量進(jìn)行多目標(biāo)優(yōu)化的基礎(chǔ)上,對普通優(yōu)化變量采取單參數(shù)掃描法尋優(yōu),由表3的靈敏度分析結(jié)果可知,wy、wst、hsec和wend對直線電機(jī)的電磁推力和永磁體利用率的影響較低,相比之下對推力波動的靈敏度大于0.1。因此,以推力波動最小為優(yōu)化目標(biāo),對4個(gè)普通優(yōu)化變量采取單參數(shù)掃描法優(yōu)化,推力波動的優(yōu)化結(jié)果如圖13所示。

    圖13 單參數(shù)掃描法優(yōu)化結(jié)果

    由圖13可知,普通優(yōu)化變量對推力波動的影響呈現(xiàn)非線性,根據(jù)圖13中推力波動最小值,選取普通優(yōu)化變量的參數(shù)值,結(jié)合前文顯著優(yōu)化變量優(yōu)化結(jié)果,得到最終優(yōu)化方案如表10所示。田口法、優(yōu)化方案1和優(yōu)化方案2的初級鐵心長度相比初始方案都略有增加,3個(gè)優(yōu)化方案的初級齒部寬度wt相比初始方案都有所降低,直線電機(jī)的電負(fù)荷與銅耗增加。

    表10 最終優(yōu)化方案對比

    無取向、取向硅鋼片直線電機(jī)電磁性能對比如圖14和表11所示。由表11可知,優(yōu)化方案1和優(yōu)化方案2的推力波動都低于4%。優(yōu)化方案1的電磁推力為432.09 N,相比田口法降低了3.96%,但推力波動僅為3.73%,相較田口法下降了55.01%,其永磁體利用率在5種方案中最高,永磁體用量相比田口法降低了8.26%,更適合應(yīng)用于長行程直線電機(jī),可以有效降低直線電機(jī)的制造成本。優(yōu)化方案2的永磁體用量與初始方案接近,其電磁推力為456.12 N,相較于初始方案提高8.96%,相比傳統(tǒng)直線電機(jī)提升13.35%,其永磁體利用率相比初始方案提升了10.63%,該方案更適合應(yīng)用于對推力密度要求較高的場合。

    表11 優(yōu)化結(jié)果電磁性能對比

    圖14 優(yōu)化結(jié)果對比

    結(jié)果表明,基于靈敏度分析,采用田口法結(jié)合遺傳算法和單參數(shù)掃描法對各優(yōu)化變量進(jìn)行分層優(yōu)化的方法,可以有效提高取向硅鋼片直線電機(jī)的電磁性能,各優(yōu)化方案的綜合性能指標(biāo)均高于初始方案。

    5 結(jié) 語

    本文以一臺12槽14極傳統(tǒng)永磁直線同步電機(jī)為例,介紹了取向硅鋼片直線電機(jī)的拓?fù)浣Y(jié)構(gòu)。基于該電機(jī)的磁網(wǎng)絡(luò)模型,通過靈敏度分析對優(yōu)化參數(shù)進(jìn)行分層優(yōu)化,得到以下結(jié)論:

    1)取向硅鋼片直線電機(jī)可以提升直線電機(jī)的推力密度。

    2)磁網(wǎng)絡(luò)模型求解電機(jī)的電磁性能相較于有限元法偏差在5%以內(nèi),有效降低了計(jì)算的時(shí)間成本,為取向硅鋼片直線電機(jī)的設(shè)計(jì)和優(yōu)化奠定了基礎(chǔ)。

    3)采用田口法結(jié)合遺傳算法和單參數(shù)掃描法對各優(yōu)化變量進(jìn)行分層優(yōu)化,可以有效提高取向硅鋼片直線電機(jī)的綜合性能指標(biāo)。通過改變本文目標(biāo)函數(shù)的權(quán)重系數(shù),可以根據(jù)實(shí)際應(yīng)用需求選取最優(yōu)方案,為工程實(shí)踐提供參考。

    猜你喜歡
    硅鋼片永磁體遺傳算法
    冷軋硅鋼片特性研究
    銅和硅鋼片價(jià)格的變化對變壓器成本及設(shè)計(jì)方案的影響
    大型電機(jī)硅鋼片自動沖制機(jī)器人上下料系統(tǒng)改造
    考慮永磁體不可逆退磁的磁齒輪復(fù)合電機(jī)設(shè)計(jì)
    霍爾式輪速傳感器永磁體磁場均勻性測量方法研究
    基于不等厚永磁體的非均勻Halbach型PMSM氣隙磁場解析及性能研究
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測
    基于改進(jìn)的遺傳算法的模糊聚類算法
    丝瓜视频免费看黄片| 最近中文字幕高清免费大全6| 日本黄大片高清| 国产亚洲精品第一综合不卡 | 中文天堂在线官网| 国产乱人偷精品视频| av女优亚洲男人天堂| 极品人妻少妇av视频| 国产在线免费精品| 在线观看www视频免费| 亚洲精品久久午夜乱码| 久久免费观看电影| 亚洲五月色婷婷综合| 国产精品国产三级国产专区5o| 制服诱惑二区| 日本黄大片高清| 91精品一卡2卡3卡4卡| 亚洲av综合色区一区| 亚洲av福利一区| 在线天堂最新版资源| 国产精品麻豆人妻色哟哟久久| videossex国产| 男人操女人黄网站| 久久鲁丝午夜福利片| 亚洲成人一二三区av| 亚洲色图综合在线观看| 街头女战士在线观看网站| 亚洲精品一二三| 91精品国产国语对白视频| 成年女人在线观看亚洲视频| √禁漫天堂资源中文www| 久久久久久久亚洲中文字幕| 亚洲国产av影院在线观看| 国产 精品1| 亚洲人成77777在线视频| 校园人妻丝袜中文字幕| 亚洲丝袜综合中文字幕| 国产一区二区在线观看日韩| 国产亚洲一区二区精品| 亚洲精品乱码久久久久久按摩| 欧美 亚洲 国产 日韩一| 少妇高潮的动态图| 午夜91福利影院| 色94色欧美一区二区| 久久99热这里只频精品6学生| 国产av码专区亚洲av| 久久久国产欧美日韩av| 丰满少妇做爰视频| 亚洲欧美成人综合另类久久久| 成人免费观看视频高清| 又黄又爽又刺激的免费视频.| 老司机影院成人| 久久人人爽人人爽人人片va| 这个男人来自地球电影免费观看 | 久久国内精品自在自线图片| 欧美人与性动交α欧美精品济南到 | 免费av中文字幕在线| 亚洲精品aⅴ在线观看| 亚洲精华国产精华液的使用体验| 久热久热在线精品观看| 在线观看美女被高潮喷水网站| 久久久久精品久久久久真实原创| 欧美日韩一区二区视频在线观看视频在线| 国产一区有黄有色的免费视频| 男人爽女人下面视频在线观看| 一级二级三级毛片免费看| 日韩伦理黄色片| 免费大片黄手机在线观看| 三上悠亚av全集在线观看| 久久 成人 亚洲| 大片电影免费在线观看免费| 久久久午夜欧美精品| 伊人久久精品亚洲午夜| 精品亚洲成国产av| 久久99热6这里只有精品| 久久韩国三级中文字幕| 五月玫瑰六月丁香| 91国产中文字幕| 国产有黄有色有爽视频| 黄色欧美视频在线观看| 少妇 在线观看| 国产精品久久久久久av不卡| 色婷婷av一区二区三区视频| 欧美亚洲 丝袜 人妻 在线| 亚洲经典国产精华液单| 久热这里只有精品99| 亚洲欧美一区二区三区国产| 国产精品国产三级国产av玫瑰| 亚洲av国产av综合av卡| 狠狠婷婷综合久久久久久88av| 人人妻人人澡人人爽人人夜夜| 久久久欧美国产精品| 18禁观看日本| 水蜜桃什么品种好| 中文字幕精品免费在线观看视频 | 夫妻性生交免费视频一级片| 青春草视频在线免费观看| 丝袜在线中文字幕| 国语对白做爰xxxⅹ性视频网站| 精品一区在线观看国产| 精品国产乱码久久久久久小说| 寂寞人妻少妇视频99o| 欧美+日韩+精品| 黄色欧美视频在线观看| 18禁观看日本| 免费黄色在线免费观看| 久久午夜福利片| 色哟哟·www| av福利片在线| 91久久精品电影网| 国产精品.久久久| 大又大粗又爽又黄少妇毛片口| 狠狠精品人妻久久久久久综合| 一本一本综合久久| 国产av国产精品国产| 亚洲内射少妇av| 99热国产这里只有精品6| 国产熟女欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 观看av在线不卡| 91久久精品电影网| 又大又黄又爽视频免费| 久久精品国产亚洲av涩爱| videos熟女内射| 热re99久久精品国产66热6| 高清欧美精品videossex| 日日爽夜夜爽网站| 久久久国产一区二区| 欧美日韩精品成人综合77777| 天堂中文最新版在线下载| 91精品一卡2卡3卡4卡| 免费黄频网站在线观看国产| 美女视频免费永久观看网站| 日本欧美国产在线视频| 国产亚洲午夜精品一区二区久久| 伦理电影大哥的女人| 精品一区二区三区视频在线| av福利片在线| av专区在线播放| 久久国产亚洲av麻豆专区| 在线观看免费日韩欧美大片 | 少妇高潮的动态图| 国产精品国产三级国产专区5o| 一本色道久久久久久精品综合| 精品酒店卫生间| 妹子高潮喷水视频| 男女边吃奶边做爰视频| 国产永久视频网站| 高清毛片免费看| 午夜免费鲁丝| 母亲3免费完整高清在线观看 | 国产有黄有色有爽视频| 亚洲无线观看免费| av视频免费观看在线观看| 一级毛片aaaaaa免费看小| av在线app专区| 韩国av在线不卡| av不卡在线播放| 成人无遮挡网站| 久久韩国三级中文字幕| 国产精品一区二区在线不卡| √禁漫天堂资源中文www| 亚洲av福利一区| 久久久精品94久久精品| 久久午夜福利片| 国产乱来视频区| 秋霞在线观看毛片| 成人国语在线视频| 你懂的网址亚洲精品在线观看| 一个人看视频在线观看www免费| 国产精品人妻久久久久久| 最近中文字幕2019免费版| 色网站视频免费| videosex国产| 国产免费现黄频在线看| 免费高清在线观看日韩| 亚洲欧美一区二区三区黑人 | 欧美另类一区| 国产精品欧美亚洲77777| 亚洲av男天堂| 美女国产视频在线观看| 久久精品久久久久久噜噜老黄| 国产精品成人在线| 成人亚洲精品一区在线观看| 亚洲国产成人一精品久久久| av在线播放精品| 免费播放大片免费观看视频在线观看| 一个人看视频在线观看www免费| 999精品在线视频| 久久人人爽av亚洲精品天堂| 久久久午夜欧美精品| av视频免费观看在线观看| 精品一区二区三区视频在线| 大片电影免费在线观看免费| 日韩av在线免费看完整版不卡| 精品视频人人做人人爽| 狂野欧美激情性bbbbbb| 国产亚洲午夜精品一区二区久久| 男女边摸边吃奶| 亚洲人成网站在线观看播放| 秋霞在线观看毛片| 免费观看性生交大片5| 美女国产视频在线观看| 一区在线观看完整版| 制服诱惑二区| 日本午夜av视频| 日本av免费视频播放| 69精品国产乱码久久久| 欧美bdsm另类| 69精品国产乱码久久久| 午夜av观看不卡| 欧美人与性动交α欧美精品济南到 | 九九爱精品视频在线观看| 三级国产精品片| 一级黄片播放器| 51国产日韩欧美| 曰老女人黄片| 亚洲国产色片| 精品人妻熟女av久视频| 国内精品宾馆在线| 日本欧美国产在线视频| 下体分泌物呈黄色| 国产精品久久久久久精品电影小说| 在线 av 中文字幕| 在线 av 中文字幕| 国产黄色视频一区二区在线观看| 亚洲五月色婷婷综合| 性色av一级| 热re99久久精品国产66热6| 婷婷成人精品国产| 国产精品久久久久久久久免| 国产精品久久久久久久久免| 香蕉精品网在线| av线在线观看网站| 乱人伦中国视频| 一个人免费看片子| 久久久精品免费免费高清| 少妇的逼水好多| 欧美日本中文国产一区发布| 欧美激情极品国产一区二区三区 | 亚洲人成网站在线播| 高清毛片免费看| 十八禁高潮呻吟视频| 又黄又爽又刺激的免费视频.| 国产日韩欧美亚洲二区| 精品亚洲成国产av| 日韩精品免费视频一区二区三区 | 亚洲欧洲日产国产| 免费av中文字幕在线| 黑人猛操日本美女一级片| 国产亚洲最大av| 久久综合国产亚洲精品| 久久ye,这里只有精品| 亚洲国产精品一区三区| 在线观看国产h片| 中文字幕最新亚洲高清| av播播在线观看一区| 国产女主播在线喷水免费视频网站| 只有这里有精品99| 哪个播放器可以免费观看大片| √禁漫天堂资源中文www| 国语对白做爰xxxⅹ性视频网站| 夜夜骑夜夜射夜夜干| 人妻人人澡人人爽人人| 亚洲综合精品二区| 精品人妻熟女av久视频| 欧美精品亚洲一区二区| 日韩 亚洲 欧美在线| 精品一区二区三区视频在线| 日韩一区二区三区影片| 麻豆乱淫一区二区| 99热6这里只有精品| 水蜜桃什么品种好| 欧美日韩av久久| 一级黄片播放器| 建设人人有责人人尽责人人享有的| 亚洲欧美日韩另类电影网站| 日本欧美国产在线视频| 我要看黄色一级片免费的| 久久97久久精品| 黄片播放在线免费| 成人漫画全彩无遮挡| 哪个播放器可以免费观看大片| 三级国产精品欧美在线观看| www.av在线官网国产| 日韩成人伦理影院| 国产日韩欧美亚洲二区| 亚洲国产毛片av蜜桃av| 久久精品国产鲁丝片午夜精品| 亚洲精品乱码久久久v下载方式| 国产深夜福利视频在线观看| 久久久a久久爽久久v久久| 国产精品麻豆人妻色哟哟久久| 国产成人精品一,二区| 欧美一级a爱片免费观看看| 成人午夜精彩视频在线观看| 在线播放无遮挡| 亚洲在久久综合| 九草在线视频观看| 久久久国产精品麻豆| 女人久久www免费人成看片| 国产在线视频一区二区| 欧美97在线视频| 午夜免费男女啪啪视频观看| 久热这里只有精品99| 日日摸夜夜添夜夜添av毛片| 国产成人一区二区在线| 免费大片黄手机在线观看| 国产在线免费精品| 成人免费观看视频高清| av有码第一页| 肉色欧美久久久久久久蜜桃| 久久99一区二区三区| 亚洲国产精品一区二区三区在线| 交换朋友夫妻互换小说| 久久久亚洲精品成人影院| 日日摸夜夜添夜夜添av毛片| 日日啪夜夜爽| 免费看av在线观看网站| 国产 精品1| 大片免费播放器 马上看| 国产精品免费大片| 少妇被粗大猛烈的视频| 免费观看av网站的网址| 赤兔流量卡办理| 九九爱精品视频在线观看| 免费日韩欧美在线观看| 99热国产这里只有精品6| 日韩熟女老妇一区二区性免费视频| 欧美97在线视频| 欧美丝袜亚洲另类| 亚洲国产最新在线播放| 国产精品三级大全| 男人操女人黄网站| 欧美激情极品国产一区二区三区 | 插阴视频在线观看视频| 国语对白做爰xxxⅹ性视频网站| 久久韩国三级中文字幕| 精品卡一卡二卡四卡免费| 黑人欧美特级aaaaaa片| 国产视频首页在线观看| 好男人视频免费观看在线| 热99国产精品久久久久久7| 97精品久久久久久久久久精品| a级毛色黄片| 美女cb高潮喷水在线观看| 人成视频在线观看免费观看| 91aial.com中文字幕在线观看| 久久97久久精品| 国产欧美亚洲国产| 在线观看免费视频网站a站| 国产精品人妻久久久久久| 一个人看视频在线观看www免费| 国产精品久久久久久久电影| 这个男人来自地球电影免费观看 | 久久亚洲国产成人精品v| 精品人妻熟女毛片av久久网站| 国产精品久久久久久久电影| 只有这里有精品99| 成人国语在线视频| 校园人妻丝袜中文字幕| 黄色毛片三级朝国网站| 国产日韩欧美视频二区| 在线亚洲精品国产二区图片欧美 | 精品熟女少妇av免费看| 超色免费av| 伊人亚洲综合成人网| 欧美最新免费一区二区三区| a级片在线免费高清观看视频| 国产高清三级在线| 色婷婷久久久亚洲欧美| 在线观看美女被高潮喷水网站| 亚洲一级一片aⅴ在线观看| 成人手机av| 国产av码专区亚洲av| 日本黄色日本黄色录像| videos熟女内射| 国产亚洲一区二区精品| 亚洲av免费高清在线观看| 热re99久久精品国产66热6| 香蕉精品网在线| 伊人久久精品亚洲午夜| √禁漫天堂资源中文www| 久久人人爽人人爽人人片va| 亚洲av免费高清在线观看| 国产亚洲av片在线观看秒播厂| 女性生殖器流出的白浆| 观看av在线不卡| 视频中文字幕在线观看| 少妇猛男粗大的猛烈进出视频| 热re99久久精品国产66热6| 考比视频在线观看| 男女边摸边吃奶| 国产日韩一区二区三区精品不卡 | 黑人猛操日本美女一级片| 热99国产精品久久久久久7| 丝袜在线中文字幕| 18禁在线无遮挡免费观看视频| 一级黄片播放器| 成年av动漫网址| 亚洲丝袜综合中文字幕| 狠狠婷婷综合久久久久久88av| 日本色播在线视频| 青春草亚洲视频在线观看| av电影中文网址| 免费黄色在线免费观看| 51国产日韩欧美| 九九爱精品视频在线观看| 欧美变态另类bdsm刘玥| 欧美日韩综合久久久久久| 国产黄色免费在线视频| 久久久久久久久大av| 有码 亚洲区| 毛片一级片免费看久久久久| 成年人午夜在线观看视频| 在线 av 中文字幕| 麻豆精品久久久久久蜜桃| 另类精品久久| 久久国产精品男人的天堂亚洲 | 搡老乐熟女国产| 国产免费视频播放在线视频| 亚洲欧美中文字幕日韩二区| 亚洲av综合色区一区| 亚洲图色成人| 精品99又大又爽又粗少妇毛片| 亚洲精品一区蜜桃| 亚洲,欧美,日韩| 日韩av不卡免费在线播放| 日本vs欧美在线观看视频| av在线播放精品| 最近2019中文字幕mv第一页| 中文字幕制服av| 高清在线视频一区二区三区| 亚洲人成网站在线播| 简卡轻食公司| a级毛片黄视频| 女性被躁到高潮视频| 午夜福利视频精品| 一区二区三区四区激情视频| 欧美性感艳星| 日本爱情动作片www.在线观看| 男女高潮啪啪啪动态图| 亚洲情色 制服丝袜| 丰满迷人的少妇在线观看| 99九九线精品视频在线观看视频| 人人澡人人妻人| 啦啦啦视频在线资源免费观看| 日本av手机在线免费观看| 中文字幕av电影在线播放| 黄色怎么调成土黄色| 91久久精品国产一区二区成人| 中文字幕av电影在线播放| 在线观看美女被高潮喷水网站| 精品卡一卡二卡四卡免费| 丝袜在线中文字幕| 欧美97在线视频| 日本与韩国留学比较| 久久午夜综合久久蜜桃| 亚洲欧美日韩另类电影网站| 精品人妻熟女毛片av久久网站| 亚洲国产av影院在线观看| 秋霞在线观看毛片| 熟女人妻精品中文字幕| 久久国产精品男人的天堂亚洲 | 另类精品久久| 女的被弄到高潮叫床怎么办| 日韩在线高清观看一区二区三区| 国产av码专区亚洲av| 日本午夜av视频| av又黄又爽大尺度在线免费看| 日韩大片免费观看网站| 日本av手机在线免费观看| 久久精品夜色国产| 美女中出高潮动态图| 永久网站在线| 亚洲精品色激情综合| 欧美 亚洲 国产 日韩一| 免费少妇av软件| 三级国产精品欧美在线观看| 国产男人的电影天堂91| 在现免费观看毛片| 亚洲丝袜综合中文字幕| 三级国产精品欧美在线观看| 777米奇影视久久| 国产精品嫩草影院av在线观看| 国产亚洲av片在线观看秒播厂| 日韩欧美一区视频在线观看| 免费观看av网站的网址| 搡女人真爽免费视频火全软件| 日日爽夜夜爽网站| 卡戴珊不雅视频在线播放| 三级国产精品片| 精品亚洲成a人片在线观看| 不卡视频在线观看欧美| 街头女战士在线观看网站| 91在线精品国自产拍蜜月| 久久久久久久大尺度免费视频| 欧美日韩综合久久久久久| 亚洲欧美成人综合另类久久久| 黄色怎么调成土黄色| 欧美成人精品欧美一级黄| 性高湖久久久久久久久免费观看| 亚洲怡红院男人天堂| 最近手机中文字幕大全| 高清毛片免费看| videosex国产| 亚洲av.av天堂| 亚洲欧美清纯卡通| 午夜影院在线不卡| 国产免费一级a男人的天堂| 九色成人免费人妻av| 最近2019中文字幕mv第一页| 亚洲精品久久成人aⅴ小说 | 18在线观看网站| 亚洲av欧美aⅴ国产| 久久久久久久久久久丰满| 国国产精品蜜臀av免费| 国产亚洲最大av| 黄色毛片三级朝国网站| 国语对白做爰xxxⅹ性视频网站| 国产色婷婷99| 我的老师免费观看完整版| 亚洲天堂av无毛| 一二三四中文在线观看免费高清| 精品99又大又爽又粗少妇毛片| 亚洲成人av在线免费| 亚洲美女视频黄频| 啦啦啦中文免费视频观看日本| 纯流量卡能插随身wifi吗| 精品一品国产午夜福利视频| 国产女主播在线喷水免费视频网站| 激情五月婷婷亚洲| 国产日韩欧美视频二区| 亚洲精华国产精华液的使用体验| 欧美亚洲 丝袜 人妻 在线| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产av蜜桃| 久久久久久久久久人人人人人人| 爱豆传媒免费全集在线观看| 亚洲国产精品成人久久小说| 欧美老熟妇乱子伦牲交| 五月伊人婷婷丁香| 一本大道久久a久久精品| 国产乱人偷精品视频| 国产一区二区三区综合在线观看 | 美女视频免费永久观看网站| kizo精华| 男的添女的下面高潮视频| 国产精品欧美亚洲77777| 在线观看www视频免费| 久久久国产精品麻豆| 久久久国产一区二区| 最近的中文字幕免费完整| 国产精品一区二区三区四区免费观看| 观看av在线不卡| 一区二区三区乱码不卡18| 久热这里只有精品99| 亚洲av国产av综合av卡| 亚洲欧美清纯卡通| 国产高清三级在线| 欧美亚洲 丝袜 人妻 在线| 亚洲高清免费不卡视频| 一级爰片在线观看| 免费人妻精品一区二区三区视频| 精品少妇久久久久久888优播| 中文字幕亚洲精品专区| 日本免费在线观看一区| 黑人猛操日本美女一级片| 如何舔出高潮| 日本免费在线观看一区| 免费av中文字幕在线| 内地一区二区视频在线| 免费观看av网站的网址| 久久久精品94久久精品| 亚洲人成网站在线播| 你懂的网址亚洲精品在线观看| 26uuu在线亚洲综合色| 91在线精品国自产拍蜜月| 美女中出高潮动态图| 国产一区二区在线观看日韩| 亚洲四区av| 大陆偷拍与自拍| 欧美性感艳星| 一级爰片在线观看| av在线观看视频网站免费| 国产视频首页在线观看| 国产欧美日韩一区二区三区在线 | 亚洲精品国产色婷婷电影| 熟妇人妻不卡中文字幕| 免费观看av网站的网址| 精品少妇久久久久久888优播| 国产一级毛片在线| 国产免费福利视频在线观看| 国国产精品蜜臀av免费| 美女中出高潮动态图| 少妇丰满av| 亚洲国产色片| www.av在线官网国产| xxx大片免费视频| 欧美3d第一页| 肉色欧美久久久久久久蜜桃| 天堂8中文在线网| 国产成人精品一,二区| 日韩av不卡免费在线播放| 中文字幕最新亚洲高清| 亚洲人成网站在线播| 我的老师免费观看完整版| 国产av国产精品国产| 国产色婷婷99| 一区二区三区精品91| av不卡在线播放| 亚洲精品美女久久av网站| 色哟哟·www| av在线app专区| 69精品国产乱码久久久|