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

    考慮諧波電流的驅(qū)動(dòng)電機(jī)振動(dòng)噪聲優(yōu)化

    2023-06-08 09:21:52戴佳錢曾發(fā)林徐求福王佳圣

    戴佳錢,曾發(fā)林,徐求福,王佳圣

    (江蘇大學(xué) 汽車工程研究院, 江蘇 鎮(zhèn)江 212013)

    0 引言

    電動(dòng)汽車驅(qū)動(dòng)電機(jī)的振動(dòng)噪聲是整車振動(dòng)噪聲的主要組成部分[1]。驅(qū)動(dòng)電機(jī)的振動(dòng)噪聲特性與傳統(tǒng)內(nèi)燃機(jī)有較大差別,除空氣動(dòng)力噪聲與機(jī)械噪聲外,由電機(jī)永磁體磁密與電樞磁密相互作用產(chǎn)生的具有階次特性的徑向電磁力是引起電機(jī)殼體振動(dòng)并輻射產(chǎn)生噪聲的主要原因[2]。通過優(yōu)化電機(jī)徑向電磁力改善電機(jī)振動(dòng)噪聲具有十分重要的意義。

    國內(nèi)外研究對(duì)電機(jī)徑向電磁力的優(yōu)化方法主要涉及結(jié)構(gòu)優(yōu)化與電流諧波優(yōu)化。文獻(xiàn)[3]分析總結(jié)了不同極槽配合和繞組層數(shù)電機(jī)的最低徑向力波階數(shù),并指出力波階數(shù)小的極槽配合更易引起較大的振動(dòng)。文獻(xiàn)[4]推導(dǎo)了轉(zhuǎn)子分段斜極的徑向力波解析式并通過調(diào)整分段數(shù)降低了電機(jī)的振動(dòng)噪聲。然而結(jié)構(gòu)優(yōu)化一般在電機(jī)本體設(shè)計(jì)階段進(jìn)行考慮,并且無法動(dòng)態(tài)優(yōu)化徑向電磁力。優(yōu)化諧波電流可以針對(duì)性地調(diào)整徑向電磁力的時(shí)空分布,文獻(xiàn)[5]提出了一種電機(jī)電磁振動(dòng)噪聲半解析模型,并分析了諧波電流對(duì)電機(jī)振動(dòng)噪聲的影響,為針對(duì)諧波電流的優(yōu)化提供了參考。文獻(xiàn)[6]分析了注入諧波電流優(yōu)化徑向電磁力的原理,并推導(dǎo)出了永磁同步電機(jī)徑向電磁力波的諧波電流補(bǔ)償模型。文獻(xiàn)[7]基于多倍頻電流注入法的徑向振動(dòng)抑制模型,設(shè)計(jì)優(yōu)化算法實(shí)現(xiàn)了各頻次徑向振動(dòng)的綜合抑制。文獻(xiàn)[8]針對(duì)分?jǐn)?shù)槽電機(jī)中的低模數(shù)電磁力諧波,采用注入補(bǔ)償電流的方式降低了電機(jī)的電磁振動(dòng)。文獻(xiàn)[9]分析指出對(duì)整數(shù)槽永磁同步電機(jī)振動(dòng)噪聲影響最大的是空間0階次徑向電磁力波。文獻(xiàn)[10]分析了某車用永磁同步電機(jī)0階徑向電磁力波的主要來源,通過注入13次諧波電流的方法削弱了72階次噪聲。

    本文以某8極48槽永磁同步驅(qū)動(dòng)電機(jī)為研究對(duì)象,采用優(yōu)化電機(jī)電流諧波的方式降低振動(dòng)噪聲。采集電機(jī)噪聲信號(hào)并分析各階次徑向電磁力在不同轉(zhuǎn)速下對(duì)電機(jī)振動(dòng)噪聲的影響,選擇3 000 r/min轉(zhuǎn)速為研究工況,使用實(shí)驗(yàn)噪聲信號(hào)驗(yàn)證建立的磁固聲耦合仿真模型。建立考慮電流諧波的徑向電磁力解析模型,基于實(shí)驗(yàn)采樣電流分析各部分徑向電磁力時(shí)空特征。通過有限元方法仿真計(jì)算各定子齒至測(cè)點(diǎn)處的噪聲傳遞函數(shù),同時(shí)采用柯特斯公式計(jì)算各部分徑向電磁力波作用于各定子齒的集中力,結(jié)合噪聲傳遞函數(shù)與集中力基于線性疊加法得到測(cè)點(diǎn)處噪聲預(yù)測(cè)模型??紤]諧波電流對(duì)轉(zhuǎn)矩脈動(dòng)的影響,將噪聲預(yù)測(cè)模型代入遺傳算法對(duì)擬注入諧波電流的幅值與相位進(jìn)行尋優(yōu)。仿真結(jié)果表明:優(yōu)化后的諧波電流能夠在不加劇轉(zhuǎn)矩脈動(dòng)的前提下有效降低電機(jī)的振動(dòng)噪聲。

    1 驅(qū)動(dòng)電機(jī)振動(dòng)噪聲信號(hào)采集及分析

    1.1 實(shí)驗(yàn)設(shè)計(jì)與實(shí)施

    本文研究對(duì)象為某永磁同步驅(qū)動(dòng)電機(jī),其主要性能參數(shù)如表1所示。在半消聲室中將被測(cè)電機(jī)聯(lián)接至測(cè)功機(jī),并如圖1所示在距離電機(jī)外包絡(luò)面0.5 m測(cè)量面上布置傳聲器。采集電機(jī)在50%負(fù)載下從100 r/min勻加速至6 000 r/min時(shí)的瞬態(tài)振動(dòng)噪聲信號(hào),以及電機(jī)在各穩(wěn)態(tài)轉(zhuǎn)速工況下的振動(dòng)噪聲信號(hào),穩(wěn)態(tài)轉(zhuǎn)速以1 000 r/min為間隔選取1 000~6 000 r/min中的相應(yīng)轉(zhuǎn)速。

    表1 驅(qū)動(dòng)電機(jī)主要性能參數(shù)

    圖1 驅(qū)動(dòng)電機(jī)噪聲信號(hào)采集測(cè)點(diǎn)布置

    1.2 實(shí)驗(yàn)數(shù)據(jù)處理與分析

    對(duì)勻加速工況下各振動(dòng)噪聲信號(hào)進(jìn)行處理,得到全轉(zhuǎn)速振動(dòng)噪聲信號(hào)頻譜圖,各測(cè)點(diǎn)的振動(dòng)噪聲信號(hào)在頻譜圖上表現(xiàn)出相似的特征,圖2所示為Front測(cè)點(diǎn)全轉(zhuǎn)速噪聲信號(hào)頻譜圖。

    圖2 Front采集點(diǎn)噪聲信號(hào)頻譜圖

    電機(jī)噪聲存在明顯的階次特征,對(duì)于本文所研究的8極48槽驅(qū)動(dòng)電機(jī),徑向電磁力的主要階次等于極數(shù)的整數(shù)倍次,而8階次作為最低階次,主要由永磁體基波磁場(chǎng)作用產(chǎn)生[3],較難采用優(yōu)化諧波電流的方式減小,所以這種情況不予考慮。由于電機(jī)本身的模態(tài)特性,在頻譜圖上可以看到存在400~860 Hz和1 430~1 780 Hz 2個(gè)明顯的共振帶,較低頻的共振帶處于機(jī)械噪聲與較低階次的徑向電磁力活躍的頻段內(nèi),而在較高頻率的共振帶區(qū)間內(nèi),24、32、40階次的徑向電磁力所引起的噪聲有明顯的加強(qiáng)。根據(jù)各階次徑向電磁力的影響嚴(yán)重程度,本文選擇比較具有代表性的3 000 r/min轉(zhuǎn)速作為優(yōu)化研究工況。

    1.3 磁固聲仿真平臺(tái)建立與驗(yàn)證

    測(cè)點(diǎn)處的電機(jī)噪聲需要經(jīng)過磁固聲耦合仿真計(jì)算得到,本文的仿真流程如圖3所示。將 3 000 r/min轉(zhuǎn)速工況下采樣得到的電流輸入二維電機(jī)模型中計(jì)算得到徑向電磁力,施加由節(jié)點(diǎn)力法計(jì)算得到的集中力至電機(jī)定子齒,將得到的電機(jī)殼體響應(yīng)代入聲腔邊界元模型中計(jì)算得到測(cè)點(diǎn)噪聲。

    圖3 磁固聲仿真流程

    其中,二維電機(jī)有限元模型根據(jù)表1參數(shù)建立,如圖4(a)所示,電機(jī)殼體聲腔模型如圖4(b)所示。

    圖4 電機(jī)有限元及殼體聲腔邊界元模型示意圖

    將仿真計(jì)算得到的噪聲信號(hào)與實(shí)驗(yàn)采集得到的噪聲信號(hào)進(jìn)行對(duì)比,如圖5所示。仿真信號(hào)在0~800 Hz之間的低頻段與實(shí)驗(yàn)信號(hào)存在較大差異,因?yàn)樵擃l段的噪聲主要由軸承滾珠沖擊和電機(jī)轉(zhuǎn)子偏心等機(jī)械原因引起,對(duì)于研究所關(guān)注的低階次徑向電磁力所引起的噪聲,仿真信號(hào)與實(shí)驗(yàn)信號(hào)保持了較高的一致性,說明該磁固聲仿真模型具有較高的準(zhǔn)確性,可用于后續(xù)研究。

    圖5 Front測(cè)點(diǎn)處噪聲仿真信號(hào)與實(shí)驗(yàn)信號(hào)

    2 徑向電磁力解析模型建立與階次分析

    2.1 徑向電磁力解析模型

    磁勢(shì)磁導(dǎo)法是電機(jī)主要的解析建模方法之一[11]。本文研究的驅(qū)動(dòng)電機(jī)為永磁同步電機(jī),其磁動(dòng)勢(shì)由永磁體與繞組電流電樞反應(yīng)兩部分疊加得來,永磁體與單相繞組電樞反應(yīng)所引起磁動(dòng)勢(shì)的解析式可以由式(1)、式(2)分別表示。永磁體與單相繞組電樞反應(yīng)磁場(chǎng)產(chǎn)生的磁動(dòng)勢(shì)均含有空間2k-1次諧波,并且永磁體所引起磁動(dòng)勢(shì)的時(shí)間諧波來源于轉(zhuǎn)子運(yùn)動(dòng)時(shí)空間諧波的變化,而電樞反應(yīng)引起磁動(dòng)勢(shì)的時(shí)間諧波來源于繞組輸入諧波電流隨時(shí)間的變化,這使得電樞反應(yīng)所引起磁動(dòng)勢(shì)的時(shí)間屬性獨(dú)立于空間屬性。

    (1)

    (2)

    u,v=2k-1,k=0,1,2,3,…

    (3)

    永磁體與電樞反應(yīng)所引起磁動(dòng)勢(shì)的各階次幅值可由式(4)與式(5)分別表示

    (4)

    (5)

    式中:δ為氣隙長度;Br為永磁體剩磁;τm為極弧角;μ0為真空磁導(dǎo)率;kdv為繞組分布因數(shù);kpv為節(jié)距因數(shù);c為繞組層數(shù);p為極對(duì)數(shù);q為每極相槽數(shù);N為線圈匝數(shù);a為并聯(lián)支路數(shù);m為電機(jī)相數(shù)。

    為將電樞反應(yīng)所引起磁勢(shì)的時(shí)間與空間屬性解耦,聯(lián)立式(2)(5)兩式并改寫如式(6)。式(6)將電樞反應(yīng)所引起的磁動(dòng)勢(shì)理解為繞組中通入的電流與單位電流下繞組產(chǎn)生的各空間階次磁動(dòng)勢(shì)之積。

    [Fhv_unitcos(vpα+φv)]

    (6)

    式中,Fhv_unit為單位電流作用下單相繞組產(chǎn)生的各空間階次磁動(dòng)勢(shì)幅值。

    氣隙磁密由磁動(dòng)勢(shì)與氣隙磁導(dǎo)如式(7)所示相乘得到,而計(jì)算定子開槽后的氣隙磁動(dòng)勢(shì)較為復(fù)雜[12],為獲得較為精確的考慮定子開槽的磁動(dòng)勢(shì)幅值,本文施加單位電流至有限元模型中計(jì)算獲得單相繞組產(chǎn)生的各空間階次磁密。而三相繞組產(chǎn)生的磁密疊加后的表達(dá)式如式(8)所示。

    (7)

    (8)

    不同于繞組的靜態(tài)分布,永磁體磁動(dòng)勢(shì)隨轉(zhuǎn)子轉(zhuǎn)動(dòng),所以需要計(jì)算氣隙比磁導(dǎo)以得到各轉(zhuǎn)動(dòng)角度位置的永磁體磁密分布[13]。本文參考文獻(xiàn)[5]中提出的考慮徑向與切向磁密相互貢獻(xiàn),通過有限元方法計(jì)算得到開槽前后的徑向與切向磁密空間分布,根據(jù)式(9)計(jì)算得到徑向比磁導(dǎo)。忽略切向磁密的貢獻(xiàn),永磁體所引起的徑向磁密可由式(10)表示。

    (9)

    Brp(α,t)=Brps(α,t)λr(α)

    (10)

    式中:λr為徑向比磁導(dǎo);Btp為開槽后切向磁密;Brps為未開槽時(shí)的徑向磁密;Btps為未開槽時(shí)的切向磁密。

    根據(jù)麥克斯韋張量法,當(dāng)不考慮幅值較小的氣隙磁密切向分量時(shí),徑向電磁力波的表達(dá)式如式(11)所示。

    (11)

    將采樣電流分別輸入解析模型計(jì)算徑向電磁力,圖6為計(jì)算得到的氣隙徑向電磁力時(shí)空分布圖。將解析模型計(jì)算得到的結(jié)果與有限元模型進(jìn)行時(shí)間與空間2個(gè)維度上的對(duì)比,如圖7所示??梢钥吹?解析模型計(jì)算結(jié)果與有限元模型相近,證明了解析模型的準(zhǔn)確性。

    圖6 解析模型計(jì)算所得徑向電磁力時(shí)空分布圖

    圖7 空間與時(shí)間維度下解析模型與 有限元模型徑向電磁力

    2.2諧波電流對(duì)徑向電磁力貢獻(xiàn)分析

    諧波電流通過影響電樞反應(yīng)磁場(chǎng)影響電機(jī)徑向電磁力。為清楚地表示諧波電流引起的各時(shí)空階次磁場(chǎng)如何作用產(chǎn)生相應(yīng)時(shí)空階次的徑向電磁力,將式(1)與式(2)代入式(11),并將諧波電流能夠影響的徑向電磁力來源列出得到式(12)所示的電樞反應(yīng)磁場(chǎng)相互作用表達(dá)式,及式(13)所示電樞反應(yīng)磁場(chǎng)與永磁體磁場(chǎng)相互作用表達(dá)式。

    {cos((h1+h2)ω0t?(v1+v2)pα+

    (φh1+φh2)?(φv1+φv2))+

    cos((h1-h2)ω0t?(v1-v2)pθ+

    (φh1-φh2)?(φv1-φv2))+

    2cos((h1+h2)ω0t+(v1-v2)pθ+

    φh1+φh2+φv1-φv2)+

    2cos((h1-h2)ω0t+(v1+v2)pθ+

    φh1-φh2+φv1+φv2)}

    (12)

    {cos((h+u)ω0t?(v+u)pθ+

    (φh+φu)?(φv+φu))+

    cos((h-u)ω0t?(v-u)pθ+

    (φh-φu)?(φv-φu))+

    cos((h+u)ω0t?(v-u)pθ+

    (φh+φu)?(φv-φu))+

    cos((h+u)ω0t+(u-v)pθ+

    (φh+φu)+(φv-φu))+

    cos((h-u)ω0t+(u-v)pθ+

    (φh+φu)+(φv-φu))}

    (13)

    歸納得到如表2所示基波與諧波電流引起的磁動(dòng)勢(shì)對(duì)各時(shí)間階次徑向電磁力的貢獻(xiàn)情況,需要注意的是,由于幅值較小,這里沒有考慮諧波電流引起的磁動(dòng)勢(shì)之間相互作用的部分。在本文所研究的3 000 r/min轉(zhuǎn)速工況下,對(duì)測(cè)點(diǎn)噪聲影響最大的是8p及6p階次徑向電磁力,其主要由5次與7次諧波電流引起。

    表2 各磁動(dòng)勢(shì)分量對(duì)各階次徑向電磁力貢獻(xiàn)

    圖8為3 000 r/min穩(wěn)態(tài)轉(zhuǎn)速工況下電機(jī)采樣電流FFT分析頻譜圖,可以發(fā)現(xiàn)電機(jī)的主要諧波電流為5、7、11次,表3為三相電流中5、7次諧波的幅值與相位。

    圖8 3 000 r/min下PMSM A相電流頻譜

    表3 三相電流各階次諧波參數(shù)

    使用解析模型分別計(jì)算各磁動(dòng)勢(shì)產(chǎn)生的8p與6p時(shí)間階次徑向電磁力,并表示在極坐標(biāo)系上,如圖9所示。對(duì)于6p時(shí)間階次的徑向電磁力,其空間構(gòu)成較為豐富,包括0、8與24空間階次,并且各階次均有相當(dāng)?shù)暮?其中由諧波電流引起的部分包括Fa1vFa5v、Fa5vFp1、Fa1vFa7v與Fa7vFp1。而對(duì)于8p時(shí)間階次的徑向電磁力,主要包含F(xiàn)a1vFa7v與Fa7vFp1產(chǎn)生的8空間階次分量與Fp1Fp7、Fp1Fp9、Fa1vFp7、Fa1vFp9產(chǎn)生的32空間階次分量階次,并且由于是由永磁體磁勢(shì)或基波電流磁勢(shì)作用產(chǎn)生,這部分徑向電磁力無法通過調(diào)節(jié)諧波電流的方法進(jìn)行改善。

    結(jié)合徑向電磁力的極坐標(biāo)分布,可以由2個(gè)角度對(duì)其優(yōu)化,一是減小電流幅值以減小徑向電磁力幅值,二是調(diào)節(jié)諧波電流相位以優(yōu)化各徑向電磁力之間的時(shí)間相位關(guān)系。

    3 測(cè)點(diǎn)噪聲預(yù)測(cè)模型與貢獻(xiàn)量分析

    3.1 定子齒部至噪聲測(cè)點(diǎn)傳遞函數(shù)分析

    徑向電磁力波是作用在定轉(zhuǎn)子單位面積上的壓力,而作用于定子齒部的徑向電磁力對(duì)于電機(jī)的振動(dòng)有決定性的影響[12],這部分徑向電磁力可以等效為作用在各齒冠面上的離散力,如圖10所示[14]。

    圖10 徑向電磁力作用于定子齒部示意圖

    電機(jī)殼體的振動(dòng)在一定振幅范圍內(nèi)可以視作各定子齒上作用徑向電磁力所引起振動(dòng)的線性疊加,如式(14)所示,而測(cè)點(diǎn)處的噪聲響應(yīng)即疊加后所引起振動(dòng)的輻射噪聲在測(cè)點(diǎn)處的響應(yīng),如式(15)所示。于是,各齒部激勵(lì)至測(cè)點(diǎn)處噪聲傳遞函數(shù)可以表示為式(16),為獲得測(cè)點(diǎn)處響應(yīng)對(duì)于各定子齒處激勵(lì)的傳遞函數(shù),分別將單位掃頻激勵(lì)作用于各定子齒,圖11為在某定子齒處施加單位掃頻激勵(lì)時(shí)測(cè)點(diǎn)噪聲傳遞函數(shù)曲線。

    (14)

    (15)

    (16)

    式中:Fi(ω)為各定子齒上的集中力,i=1,2,…,48;HVij(ω)為殼體某點(diǎn)處的振動(dòng)響應(yīng)對(duì)于定子齒處激勵(lì)的傳遞函數(shù);HNj(ω)為噪聲測(cè)點(diǎn)處響應(yīng)對(duì)于殼體各點(diǎn)處振動(dòng)輻射的傳遞函數(shù)。

    圖11 某一定子齒至測(cè)點(diǎn)處噪聲傳遞函數(shù)曲線

    結(jié)合求得的48個(gè)齒部至噪聲測(cè)點(diǎn)傳遞函數(shù),測(cè)點(diǎn)處噪聲預(yù)測(cè)模型由式(17)表示。

    (17)

    3.2 各時(shí)空階次徑向電磁力噪聲貢獻(xiàn)量分析

    節(jié)點(diǎn)分布力法與面集中力加載法是計(jì)算徑向電磁力波對(duì)定子齒部作用力的2種方法,本文采用文獻(xiàn)[2]提及的面集中力加載法實(shí)現(xiàn)定子齒部集中力的計(jì)算。將每一個(gè)齒面的角度區(qū)間[α1,α2]四等分,并在每一個(gè)四分之一區(qū)間[βi,βi+1](i=0,1,2,3,4)應(yīng)用柯特斯公式求解徑向電磁力波在齒面上的定積分,得到的復(fù)合柯特斯公式近似解如式(18)所示。

    (18)

    式中:lr為定子軸向長度;rr為定子中心至定子齒冠半徑;y為所求定子齒面在周向上的角度;βi+0.25=0.75βi+0.25βi+1;βi+0.5=0.5βi+0.5βi+1;βi+0.25=0.25βi+0.75βi+1。

    將各磁動(dòng)勢(shì)分量產(chǎn)生的6p與8p時(shí)間階次徑向電磁力波代入上式中進(jìn)行計(jì)算,得到如圖12所示的各時(shí)空階次徑向電磁力在某一時(shí)刻作用于各齒部的集中力。根據(jù)奈奎斯特采樣定理,本文所研究的48槽定子對(duì)于24及以下的空間階次具有低通濾波效應(yīng)[14],所以對(duì)于0、8、24空間階次的徑向電磁力波分量產(chǎn)生的作用于定子齒上的集中力能夠反映出空間階次特征,而對(duì)于32空間階次的徑向電磁力波,其產(chǎn)生的集中力的空間階次表現(xiàn)為16階次。

    圖12 各磁動(dòng)勢(shì)分量產(chǎn)生的集中力在各定子齒的分布

    將各部分集中力代入噪聲預(yù)測(cè)模型中計(jì)算,得到表4所示各部分徑向電磁力對(duì)于測(cè)點(diǎn)聲壓對(duì)應(yīng)頻率的貢獻(xiàn)量。由于測(cè)點(diǎn)聲壓是由各分量引起的噪聲疊加而成,所以表中部分分量引起的聲壓幅值大于總聲壓幅值。各階次噪聲的主要分量均由較低空間階次的分量引起,對(duì)于6p階次,0空間階次分量影響最大,而對(duì)于8p階次,并不存在0空間階次含量,但此頻率接近定子共振頻率,故8空間階次所引起的噪聲聲壓也較大。

    4 基于遺傳算法的尋優(yōu)與仿真驗(yàn)證

    諧波電流注入是一種抵消因逆變器非線性等原因而產(chǎn)生的固有諧波電流的重要方法,在抑制轉(zhuǎn)矩脈動(dòng)等應(yīng)用中具有廣泛的實(shí)踐經(jīng)驗(yàn)[15]。由于諧波電流對(duì)各時(shí)空階次徑向電磁力的影響較為復(fù)雜且相互之間存在耦合,采用遺傳算法對(duì)注入諧波電流的幅值與相位進(jìn)行尋優(yōu)。優(yōu)化思路如圖13所示,將徑向電磁力解析模型與噪聲預(yù)測(cè)模型代入遺傳算法中參與尋優(yōu),同時(shí)考慮諧波電流對(duì)電機(jī)轉(zhuǎn)矩脈動(dòng)的影響。

    4.1 轉(zhuǎn)矩脈動(dòng)的解析模型

    電機(jī)的轉(zhuǎn)矩輸出可以視作齒槽轉(zhuǎn)矩tcog以及由電流與磁鏈相互作用產(chǎn)生的電磁轉(zhuǎn)矩te的線性疊加。齒槽轉(zhuǎn)矩不受諧波電流影響,且主要成分為電流頻率的12與24階次[16],可表示為式(19)所示的傅里葉級(jí)數(shù)。

    t=te+tcog

    (19)

    tcog=Tcog12cos(12ω0t+φcog12)+

    Tv24cos(24ω0t+φcog24)

    (20)

    式中:Tcog12、Tcog24分別是12階與24階齒槽轉(zhuǎn)矩的幅值;φcog12、φcog24分別是12階與24階齒槽轉(zhuǎn)矩的相位。

    圖13 注入諧波電流參數(shù)優(yōu)化思路

    對(duì)于電磁轉(zhuǎn)矩,如式(21)所示,電流基波與磁鏈基波產(chǎn)生恒定轉(zhuǎn)矩,脈動(dòng)轉(zhuǎn)矩部分主要由電流諧波與磁鏈基波,電流基波與磁鏈諧波以及電流諧波與磁鏈諧波3部分作用產(chǎn)生[17]。

    te=1.5P(LdqIdq0+Λdq0)T×Idq0+

    1.5P((LdqIdq0+Λdq0)T×idqh+

    (21)

    Λdq0= [Λ00]T

    (22)

    Idq0= [Id0Iq0]T

    (23)

    (24)

    (25)

    式中:Ldq為交直軸電感幅值;Λdq0為交直軸基波電流幅值;idqh為交直軸諧波電流;Λdq0為交直軸基波磁鏈幅值;λdqh為交直軸諧波磁鏈。

    對(duì)比轉(zhuǎn)矩解析模型與有限元計(jì)算結(jié)果,如圖14所示,兩者在轉(zhuǎn)矩脈動(dòng)特征上十分相似,最大脈動(dòng)幅值均在21.1 N·m左右,精度足夠,可以代入尋優(yōu)算法中進(jìn)行計(jì)算。

    圖14 轉(zhuǎn)矩輸出解析模型與有限元計(jì)算結(jié)果曲線

    4.2 基于遺傳算法諧波電流參數(shù)尋優(yōu)

    由于3 000 r/min工況下對(duì)噪聲聲壓影響最大的是24階次與32階次所在頻率即1 200 Hz與1 600 Hz,于是,遺傳算法的適應(yīng)度函數(shù)表示如下:

    F=P(1 200)+P(1 600)

    (26)

    式中,P(f)表示某一頻率處的聲壓。

    能夠?qū)?4階及32階噪聲產(chǎn)生影響的是5、7次諧波電流,所以優(yōu)化參數(shù)為5、7次諧波電流的幅值與相位。

    D=[Ia5,Ib5,Ic5,Ia7,Ib7,Ic7,

    φa5,φb5,φc5,φa7,φb7,φc7]

    (27)

    參數(shù)之間滿足以下關(guān)系:

    (28)

    三相電流5、7次諧波可通過Park變換轉(zhuǎn)換至d-q軸坐標(biāo)系上的6次諧波,注入的電流在d-q坐標(biāo)系下可表示如下:

    (29)

    注入諧波電流與原有電流共同產(chǎn)生的轉(zhuǎn)矩輸出為:

    (30)

    于是,為抑制轉(zhuǎn)矩脈動(dòng)增大而設(shè)計(jì)的約束條件表示如下:

    (31)

    圖15所示為以測(cè)點(diǎn)聲壓為優(yōu)化目標(biāo),基于遺傳算法對(duì)注入諧波電流幅值相位參數(shù)進(jìn)行尋優(yōu)得到的適應(yīng)度函數(shù)。最終由遺傳算法計(jì)算得到的5、7次A相諧波電流的參數(shù)如表5所示,B、C相諧波電流可由式(28)計(jì)算得到。

    圖15 遺傳算法尋優(yōu)諧波電流參數(shù)適應(yīng)度變化曲線

    表5 注入A相諧波電流最優(yōu)參數(shù)

    4.3 基于磁固聲仿真平臺(tái)的驗(yàn)證

    為驗(yàn)證注入優(yōu)化后諧波電流對(duì)各階次噪聲及轉(zhuǎn)矩脈動(dòng)的影響,將上述優(yōu)化后的諧波電流與采樣電流疊加后代入已建立的磁固聲仿真平臺(tái)中,計(jì)算得到測(cè)點(diǎn)處聲壓與此時(shí)的轉(zhuǎn)矩輸出,分別如圖16—17所示。

    表6歸納了各關(guān)注對(duì)象的優(yōu)化效果,可以看到,本文重點(diǎn)關(guān)注的6p與8p階次噪聲在注入優(yōu)化后的諧波電流后有明顯的降低,與此同時(shí),電機(jī)的轉(zhuǎn)矩脈動(dòng)也得到了有效的控制,沒有因?yàn)橹C波電流的改變而惡化。

    圖16 測(cè)點(diǎn)處響應(yīng)曲線

    圖17 諧波電流優(yōu)化前后轉(zhuǎn)矩輸出曲線

    表6 諧波電流注入前后聲壓級(jí)與轉(zhuǎn)矩脈動(dòng)

    5 結(jié)論

    1) 采集了某驅(qū)動(dòng)電機(jī)的噪聲信號(hào),分析噪聲信號(hào)的階次特征并認(rèn)為6p與8p時(shí)間階次對(duì)噪聲聲壓影響最大。選擇3 000 r/min為研究工況,建立了磁固聲仿真模型并驗(yàn)證了其在階次特征上與實(shí)驗(yàn)信號(hào)一致。

    2) 基于磁勢(shì)磁導(dǎo)法建立了徑向電磁力解析模型,通過與有限元模型計(jì)算結(jié)果的對(duì)比驗(yàn)證了其準(zhǔn)確性。根據(jù)諧波電流的階次特征分析各時(shí)空階次徑向電磁力的組成成分,并代入解析模型中計(jì)算得到各階次徑向電磁力分量的時(shí)空分布。

    3) 基于線性疊加法通過有限元方法得到定子齒至噪聲測(cè)點(diǎn)傳遞函數(shù),采用柯特斯公式計(jì)算各徑向電磁力分量作用于定子齒的集中力,實(shí)現(xiàn)各分量對(duì)噪聲測(cè)點(diǎn)聲壓的預(yù)測(cè)。基于預(yù)測(cè)模型,分析了各徑向電磁力分量對(duì)測(cè)點(diǎn)噪聲聲壓的貢獻(xiàn)量,確定了5、7次諧波電流對(duì)主要階次噪聲的影響最大。

    4) 將徑向電磁力解析模型與噪聲預(yù)測(cè)模型代入遺傳算法中,在不增大電機(jī)轉(zhuǎn)矩脈動(dòng)的前提下對(duì)擬注入的5、7次諧波電流參數(shù)進(jìn)行尋優(yōu)。仿真結(jié)果顯示,優(yōu)化后的諧波電流能夠在避免惡化轉(zhuǎn)矩脈動(dòng)的同時(shí),有效減小主要階次噪聲。

    亚洲国产精品sss在线观看 | 免费在线观看影片大全网站| 丁香六月欧美| 精品久久久精品久久久| 国产一区在线观看成人免费| 日本撒尿小便嘘嘘汇集6| 欧美激情高清一区二区三区| avwww免费| 日韩欧美在线二视频 | 韩国精品一区二区三区| 国产欧美日韩一区二区三| 国产xxxxx性猛交| 高潮久久久久久久久久久不卡| 成年版毛片免费区| 国产精品免费大片| 99热网站在线观看| a在线观看视频网站| 女人被躁到高潮嗷嗷叫费观| 国产区一区二久久| 欧美日韩视频精品一区| 午夜精品在线福利| 亚洲五月天丁香| 国产免费现黄频在线看| 99re在线观看精品视频| 超碰97精品在线观看| 亚洲国产精品合色在线| 一区二区三区国产精品乱码| 欧美中文综合在线视频| 最新在线观看一区二区三区| 三级毛片av免费| videosex国产| 亚洲成人手机| 欧美中文综合在线视频| 免费观看人在逋| 色精品久久人妻99蜜桃| 首页视频小说图片口味搜索| 少妇猛男粗大的猛烈进出视频| 亚洲精品在线美女| 91精品三级在线观看| 十八禁人妻一区二区| 夫妻午夜视频| 日本黄色视频三级网站网址 | 国产成人一区二区三区免费视频网站| 窝窝影院91人妻| 亚洲三区欧美一区| 国产一区在线观看成人免费| 国产成人精品久久二区二区91| 女人精品久久久久毛片| 国产在视频线精品| 亚洲欧美色中文字幕在线| 国产一区二区三区在线臀色熟女 | 国产高清videossex| 热99国产精品久久久久久7| 中文字幕制服av| 亚洲成人国产一区在线观看| 国产成人影院久久av| 久久久国产一区二区| 国产精品久久久久成人av| 亚洲伊人色综图| 欧美日本中文国产一区发布| 成年人黄色毛片网站| 伦理电影免费视频| 天堂俺去俺来也www色官网| 老熟女久久久| 精品国产乱子伦一区二区三区| 国产不卡一卡二| 精品久久久久久,| 久久国产乱子伦精品免费另类| 亚洲中文字幕日韩| 国产成+人综合+亚洲专区| 国产精品亚洲一级av第二区| 12—13女人毛片做爰片一| 999久久久精品免费观看国产| 热99国产精品久久久久久7| 中文字幕高清在线视频| 脱女人内裤的视频| 午夜91福利影院| 欧美精品亚洲一区二区| 超碰97精品在线观看| 无遮挡黄片免费观看| 岛国毛片在线播放| 操美女的视频在线观看| 咕卡用的链子| 国产淫语在线视频| 亚洲黑人精品在线| 亚洲avbb在线观看| 久久久精品区二区三区| 亚洲精品国产精品久久久不卡| 一进一出抽搐动态| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美一级毛片孕妇| 91九色精品人成在线观看| 精品人妻1区二区| 国产男女超爽视频在线观看| 婷婷丁香在线五月| 国产精品亚洲一级av第二区| 欧美激情极品国产一区二区三区| 性少妇av在线| 99国产精品一区二区三区| 国产一区有黄有色的免费视频| 精品免费久久久久久久清纯 | 1024香蕉在线观看| 麻豆国产av国片精品| 免费观看精品视频网站| 欧美日韩亚洲高清精品| 久久中文字幕一级| 18禁裸乳无遮挡动漫免费视频| 国产在视频线精品| 日韩欧美一区二区三区在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 亚洲aⅴ乱码一区二区在线播放 | 亚洲va日本ⅴa欧美va伊人久久| 无遮挡黄片免费观看| 亚洲av日韩在线播放| 亚洲avbb在线观看| 99久久99久久久精品蜜桃| 精品卡一卡二卡四卡免费| 午夜福利乱码中文字幕| 999久久久国产精品视频| 欧美亚洲日本最大视频资源| 在线观看一区二区三区激情| 日本精品一区二区三区蜜桃| 亚洲人成77777在线视频| 日本vs欧美在线观看视频| 国产在视频线精品| 国产精品久久久人人做人人爽| 久久久久国内视频| 免费一级毛片在线播放高清视频 | 一级毛片高清免费大全| 热99国产精品久久久久久7| 在线观看免费午夜福利视频| 国产精品国产高清国产av | 精品亚洲成a人片在线观看| 国产真人三级小视频在线观看| 精品视频人人做人人爽| 在线视频色国产色| av在线播放免费不卡| 国产有黄有色有爽视频| 伦理电影免费视频| 国产精品影院久久| 亚洲av第一区精品v没综合| 国产精品久久电影中文字幕 | а√天堂www在线а√下载 | 国产av精品麻豆| av片东京热男人的天堂| 欧美国产精品va在线观看不卡| 亚洲精品久久成人aⅴ小说| 又大又爽又粗| 亚洲国产欧美网| 久久精品国产a三级三级三级| 免费高清在线观看日韩| 免费av中文字幕在线| 一区二区三区激情视频| 黑丝袜美女国产一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美精品综合一区二区三区| 在线观看免费视频网站a站| 日韩欧美三级三区| 精品人妻在线不人妻| 可以免费在线观看a视频的电影网站| 久久影院123| 看免费av毛片| 午夜两性在线视频| 亚洲精品久久成人aⅴ小说| 叶爱在线成人免费视频播放| 免费人成视频x8x8入口观看| 男女之事视频高清在线观看| 日本a在线网址| 在线观看免费午夜福利视频| 成人影院久久| 曰老女人黄片| 国产三级黄色录像| 美女扒开内裤让男人捅视频| 久久久久久免费高清国产稀缺| 久久香蕉国产精品| 最新在线观看一区二区三区| 动漫黄色视频在线观看| 51午夜福利影视在线观看| 国产高清videossex| 亚洲全国av大片| 亚洲国产精品一区二区三区在线| 一级毛片精品| 免费少妇av软件| 欧美亚洲 丝袜 人妻 在线| 在线观看免费高清a一片| 在线观看一区二区三区激情| 在线观看免费日韩欧美大片| 欧美人与性动交α欧美软件| 99re6热这里在线精品视频| 日本wwww免费看| 国产乱人伦免费视频| www.999成人在线观看| 欧美丝袜亚洲另类 | 久久这里只有精品19| 色婷婷久久久亚洲欧美| 女性生殖器流出的白浆| 久久久久精品国产欧美久久久| 悠悠久久av| 性色av乱码一区二区三区2| 91大片在线观看| 女人高潮潮喷娇喘18禁视频| 男人操女人黄网站| 国产精品99久久99久久久不卡| 国产三级黄色录像| 成人18禁高潮啪啪吃奶动态图| 人人妻,人人澡人人爽秒播| 久久久久久久久久久久大奶| 国产激情欧美一区二区| 热99久久久久精品小说推荐| 极品人妻少妇av视频| 一级a爱视频在线免费观看| videosex国产| 久久人妻av系列| 久久久国产欧美日韩av| av福利片在线| 黑人操中国人逼视频| 一级黄色大片毛片| 欧美丝袜亚洲另类 | 嫁个100分男人电影在线观看| 免费在线观看黄色视频的| 日韩免费高清中文字幕av| 精品久久久久久,| 欧美精品一区二区免费开放| 黑丝袜美女国产一区| 三级毛片av免费| 一区二区三区国产精品乱码| 国产精品二区激情视频| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| 亚洲国产精品sss在线观看 | а√天堂www在线а√下载 | 麻豆乱淫一区二区| 欧美精品av麻豆av| 中亚洲国语对白在线视频| 一级片免费观看大全| 1024视频免费在线观看| 黄频高清免费视频| 免费观看a级毛片全部| 最新在线观看一区二区三区| 国产精品一区二区免费欧美| 91九色精品人成在线观看| 欧美 亚洲 国产 日韩一| 国产精品久久久av美女十八| 久久久久久久久久久久大奶| 高清在线国产一区| 亚洲成人免费av在线播放| а√天堂www在线а√下载 | 无遮挡黄片免费观看| 高清黄色对白视频在线免费看| 看片在线看免费视频| 日韩制服丝袜自拍偷拍| 亚洲国产看品久久| 日本精品一区二区三区蜜桃| 久久青草综合色| 久久性视频一级片| 美女 人体艺术 gogo| 国产97色在线日韩免费| 人人妻人人添人人爽欧美一区卜| av电影中文网址| 纯流量卡能插随身wifi吗| 黄色视频不卡| 亚洲人成77777在线视频| 99久久精品国产亚洲精品| 欧美最黄视频在线播放免费 | 久久国产精品影院| 19禁男女啪啪无遮挡网站| 国产精品偷伦视频观看了| 色婷婷久久久亚洲欧美| 十八禁网站免费在线| 久久久国产一区二区| 在线观看66精品国产| 色在线成人网| 99riav亚洲国产免费| 啦啦啦在线免费观看视频4| 国产av一区二区精品久久| 欧美精品一区二区免费开放| 欧美大码av| 国产精品久久久av美女十八| 成年人黄色毛片网站| 亚洲色图av天堂| 亚洲成人国产一区在线观看| 一区在线观看完整版| 久久精品亚洲精品国产色婷小说| 国产在线观看jvid| 很黄的视频免费| 99国产精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 免费不卡黄色视频| 无限看片的www在线观看| 亚洲精品美女久久久久99蜜臀| 精品一区二区三区四区五区乱码| 欧美日韩亚洲高清精品| 亚洲aⅴ乱码一区二区在线播放 | 777米奇影视久久| 激情在线观看视频在线高清 | 黄色a级毛片大全视频| 日韩中文字幕欧美一区二区| 老鸭窝网址在线观看| 精品少妇久久久久久888优播| 国产欧美日韩一区二区三| 精品熟女少妇八av免费久了| 老司机亚洲免费影院| 国产精品.久久久| 欧美成人免费av一区二区三区 | 久久香蕉国产精品| 成人三级做爰电影| av中文乱码字幕在线| 母亲3免费完整高清在线观看| 午夜成年电影在线免费观看| 天天影视国产精品| 久久香蕉激情| 久久婷婷成人综合色麻豆| 久久人妻福利社区极品人妻图片| 丁香六月欧美| av视频免费观看在线观看| 狂野欧美激情性xxxx| 中文亚洲av片在线观看爽 | 精品久久久精品久久久| 热99久久久久精品小说推荐| 国产精品秋霞免费鲁丝片| 咕卡用的链子| 久久久久精品人妻al黑| 午夜免费鲁丝| 精品人妻在线不人妻| 人妻久久中文字幕网| 欧美av亚洲av综合av国产av| 99riav亚洲国产免费| 国产精品一区二区精品视频观看| 午夜福利在线免费观看网站| 丁香六月欧美| tube8黄色片| 18在线观看网站| 人妻丰满熟妇av一区二区三区 | 精品一区二区三区视频在线观看免费 | 亚洲视频免费观看视频| av在线播放免费不卡| 国产精品自产拍在线观看55亚洲 | 日日夜夜操网爽| 国产精品98久久久久久宅男小说| 亚洲视频免费观看视频| 国产亚洲欧美98| 久久国产亚洲av麻豆专区| 深夜精品福利| 在线av久久热| 国内毛片毛片毛片毛片毛片| 人人妻人人澡人人看| 亚洲第一欧美日韩一区二区三区| 一区二区三区激情视频| 免费女性裸体啪啪无遮挡网站| 成人18禁在线播放| 日韩欧美一区二区三区在线观看 | 久久精品国产亚洲av高清一级| 色精品久久人妻99蜜桃| av片东京热男人的天堂| 中文字幕高清在线视频| av免费在线观看网站| 无遮挡黄片免费观看| 亚洲av日韩精品久久久久久密| 一级作爱视频免费观看| av一本久久久久| 国产成人精品久久二区二区免费| 岛国在线观看网站| 国产亚洲一区二区精品| 好男人电影高清在线观看| 最近最新中文字幕大全免费视频| 欧美黄色淫秽网站| 欧美日韩国产mv在线观看视频| 国产成人系列免费观看| 成人黄色视频免费在线看| 啦啦啦在线免费观看视频4| 啦啦啦 在线观看视频| 亚洲第一欧美日韩一区二区三区| 精品视频人人做人人爽| 国产精品亚洲av一区麻豆| 国产男靠女视频免费网站| 久久精品成人免费网站| 国产1区2区3区精品| 91麻豆精品激情在线观看国产 | 亚洲国产精品一区二区三区在线| 日日摸夜夜添夜夜添小说| 午夜影院日韩av| 丝袜人妻中文字幕| 婷婷成人精品国产| 久久久久久久久久久久大奶| 国产亚洲精品久久久久久毛片 | 国产免费男女视频| 人妻一区二区av| 久久人人爽av亚洲精品天堂| 在线免费观看的www视频| www.精华液| 在线永久观看黄色视频| 国产高清videossex| 亚洲第一青青草原| 50天的宝宝边吃奶边哭怎么回事| 日本黄色日本黄色录像| 亚洲av片天天在线观看| 露出奶头的视频| xxxhd国产人妻xxx| www.熟女人妻精品国产| 啦啦啦在线免费观看视频4| 三上悠亚av全集在线观看| 黄色视频,在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 一边摸一边抽搐一进一出视频| 中文字幕人妻丝袜一区二区| 国产成人一区二区三区免费视频网站| 国产精品98久久久久久宅男小说| 日韩视频一区二区在线观看| 伦理电影免费视频| 欧美亚洲日本最大视频资源| 久久久久久久久久久久大奶| 亚洲av成人一区二区三| 日韩欧美一区视频在线观看| 国产亚洲欧美精品永久| 丰满的人妻完整版| 国产成人av激情在线播放| 啦啦啦 在线观看视频| 丁香六月欧美| 亚洲人成77777在线视频| 久久性视频一级片| 变态另类成人亚洲欧美熟女 | 欧美老熟妇乱子伦牲交| 久久久国产成人免费| 亚洲第一欧美日韩一区二区三区| 男女午夜视频在线观看| 中文字幕色久视频| 女人高潮潮喷娇喘18禁视频| 成人三级做爰电影| 午夜福利视频在线观看免费| 欧美日韩亚洲国产一区二区在线观看 | 热99国产精品久久久久久7| 高清黄色对白视频在线免费看| 午夜日韩欧美国产| 国产伦人伦偷精品视频| 亚洲中文日韩欧美视频| 国产精品av久久久久免费| 中文字幕av电影在线播放| 亚洲人成77777在线视频| 久久青草综合色| 欧美精品一区二区免费开放| 国产成人免费无遮挡视频| 久久精品国产亚洲av香蕉五月 | 丝袜在线中文字幕| 亚洲av第一区精品v没综合| 国产片内射在线| 99久久人妻综合| 日本一区二区免费在线视频| 91九色精品人成在线观看| 伊人久久大香线蕉亚洲五| 黄片大片在线免费观看| av天堂久久9| 国产一区二区三区视频了| 交换朋友夫妻互换小说| 成人18禁在线播放| 999久久久精品免费观看国产| 一级片免费观看大全| 欧美精品av麻豆av| 国产精品一区二区精品视频观看| 精品国产一区二区三区久久久樱花| 国产aⅴ精品一区二区三区波| 亚洲熟妇中文字幕五十中出 | 欧美最黄视频在线播放免费 | 女同久久另类99精品国产91| 久久午夜亚洲精品久久| 777久久人妻少妇嫩草av网站| 欧美精品一区二区免费开放| 久久国产精品影院| 别揉我奶头~嗯~啊~动态视频| 国产精品一区二区精品视频观看| 午夜精品久久久久久毛片777| 九色亚洲精品在线播放| 香蕉国产在线看| 久久精品成人免费网站| 午夜视频精品福利| 亚洲一区高清亚洲精品| 午夜免费观看网址| 丝袜美足系列| 国产av一区二区精品久久| 欧美国产精品一级二级三级| 日韩三级视频一区二区三区| 亚洲 欧美一区二区三区| 久久国产精品大桥未久av| 欧美久久黑人一区二区| 亚洲人成77777在线视频| av国产精品久久久久影院| 天天添夜夜摸| 热re99久久国产66热| 亚洲aⅴ乱码一区二区在线播放 | 人妻久久中文字幕网| 亚洲精品一二三| 欧美中文综合在线视频| 成年人免费黄色播放视频| 我的亚洲天堂| 精品国产一区二区三区久久久樱花| 色播在线永久视频| 国产欧美日韩一区二区三| 欧洲精品卡2卡3卡4卡5卡区| 国产主播在线观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇一区二区三区视频日本电影| 欧美午夜高清在线| 高清在线国产一区| 国产野战对白在线观看| 国产精品亚洲av一区麻豆| 变态另类成人亚洲欧美熟女 | 久久久精品国产亚洲av高清涩受| 欧美乱码精品一区二区三区| 精品久久久久久久毛片微露脸| 黑人猛操日本美女一级片| 国产成人欧美| 制服诱惑二区| 免费在线观看影片大全网站| 久久久国产成人免费| 成人av一区二区三区在线看| 91国产中文字幕| 日本黄色视频三级网站网址 | 国产又色又爽无遮挡免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女 | 精品国产国语对白av| 一进一出好大好爽视频| 伦理电影免费视频| 亚洲av电影在线进入| 国产高清videossex| 成人国产一区最新在线观看| 精品无人区乱码1区二区| www.自偷自拍.com| 久久久久久亚洲精品国产蜜桃av| avwww免费| 最近最新中文字幕大全电影3 | 国产成人啪精品午夜网站| 亚洲国产看品久久| 欧美人与性动交α欧美精品济南到| 亚洲 欧美一区二区三区| 国产欧美日韩精品亚洲av| 中文字幕色久视频| 国产99白浆流出| 18禁美女被吸乳视频| 久久影院123| 男男h啪啪无遮挡| avwww免费| 老司机午夜十八禁免费视频| 亚洲精品在线观看二区| 黄片大片在线免费观看| 精品一区二区三区视频在线观看免费 | 丝袜美足系列| 午夜福利在线免费观看网站| 在线观看免费高清a一片| 色94色欧美一区二区| 国产精品乱码一区二三区的特点 | av欧美777| 亚洲avbb在线观看| 成人亚洲精品一区在线观看| 夜夜躁狠狠躁天天躁| 精品国产一区二区三区久久久樱花| 一本综合久久免费| 18禁美女被吸乳视频| 国产在视频线精品| 国产精品久久久久久人妻精品电影| 亚洲专区中文字幕在线| 老司机午夜十八禁免费视频| tube8黄色片| 国产成人啪精品午夜网站| 他把我摸到了高潮在线观看| 91大片在线观看| 亚洲成人免费av在线播放| 亚洲第一欧美日韩一区二区三区| 国产不卡一卡二| 国内久久婷婷六月综合欲色啪| 久久精品人人爽人人爽视色| 亚洲av片天天在线观看| 国产精华一区二区三区| 亚洲三区欧美一区| 欧美一级毛片孕妇| 国产深夜福利视频在线观看| 亚洲一区高清亚洲精品| 精品人妻熟女毛片av久久网站| 久久青草综合色| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 老司机午夜十八禁免费视频| 999久久久国产精品视频| 大片电影免费在线观看免费| 91成年电影在线观看| 色综合婷婷激情| 18禁裸乳无遮挡免费网站照片 | 日本五十路高清| 精品第一国产精品| 免费在线观看完整版高清| 757午夜福利合集在线观看| 成人亚洲精品一区在线观看| 午夜两性在线视频| 婷婷丁香在线五月| 国产精品国产av在线观看| 最近最新中文字幕大全电影3 | av中文乱码字幕在线| 亚洲情色 制服丝袜| 国产91精品成人一区二区三区| 欧美性长视频在线观看| 亚洲欧美一区二区三区黑人| 自线自在国产av| 18禁黄网站禁片午夜丰满| 久久久国产欧美日韩av| 亚洲精品av麻豆狂野| 亚洲九九香蕉| 国产欧美日韩一区二区三区在线| 国产成人av教育| 1024香蕉在线观看| 夫妻午夜视频| 欧美日本中文国产一区发布| 又紧又爽又黄一区二区| 18在线观看网站| 国产不卡一卡二| 亚洲国产中文字幕在线视频|