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

    基于熵權(quán)與混合代理模型的永磁驅(qū)動(dòng)器的優(yōu)化設(shè)計(jì)

    2016-07-16 03:22:07李召王大志時(shí)統(tǒng)宇鄭迪
    關(guān)鍵詞:熵權(quán)法

    李召, 王大志, 時(shí)統(tǒng)宇, 鄭迪

    (東北大學(xué) 信息科學(xué)與工程學(xué)院,遼寧 沈陽(yáng),110004)

    ?

    基于熵權(quán)與混合代理模型的永磁驅(qū)動(dòng)器的優(yōu)化設(shè)計(jì)

    李召,王大志,時(shí)統(tǒng)宇,鄭迪

    (東北大學(xué) 信息科學(xué)與工程學(xué)院,遼寧 沈陽(yáng),110004)

    摘要:針對(duì)永磁驅(qū)動(dòng)器(PMD)的結(jié)構(gòu)設(shè)計(jì)問(wèn)題,提出一種基于改進(jìn)熵權(quán)法結(jié)合混合代理模型的優(yōu)化設(shè)計(jì)方法。首先利用基于交叉驗(yàn)證誤差的最優(yōu)加權(quán)法,將響應(yīng)曲面法、克里金法以及支持向量機(jī)回歸結(jié)合起來(lái),構(gòu)建PMD的參數(shù)變量與響應(yīng)變量之間的混合代理模型;然后引入改進(jìn)的熵權(quán)法,將PMD的多指標(biāo)轉(zhuǎn)化為單一綜合指標(biāo),并建立其優(yōu)化的數(shù)學(xué)模型,通過(guò)自適應(yīng)權(quán)重粒子群優(yōu)化算法求解;最后對(duì)結(jié)果進(jìn)行有限元仿真分析和實(shí)驗(yàn)室仿真平臺(tái)驗(yàn)證。研究結(jié)果表明,所提出的優(yōu)化設(shè)計(jì)方法優(yōu)于其它方法,得到的PMD結(jié)構(gòu)參數(shù)合理有效,較好的實(shí)現(xiàn)了PMD的多目標(biāo)優(yōu)化設(shè)計(jì)。

    關(guān)鍵詞:永磁驅(qū)動(dòng)器;混合代理模型;熵權(quán)法;最優(yōu)加權(quán)法;自適應(yīng)權(quán)重粒子群算法

    0引言

    永磁驅(qū)動(dòng)器(permanent magnet drive,PMD)作為一種新興的節(jié)能調(diào)速設(shè)備,利用永磁材料所產(chǎn)生的磁力作用,實(shí)現(xiàn)了運(yùn)動(dòng)和轉(zhuǎn)矩?zé)o接觸傳遞,且具有維護(hù)方便、減小振動(dòng)傳遞、無(wú)諧波污染和容忍對(duì)中誤差等優(yōu)點(diǎn),已經(jīng)廣泛應(yīng)用于風(fēng)機(jī)、壓縮機(jī)以及泵類等相關(guān)傳動(dòng)系統(tǒng)中[1]。當(dāng)前,國(guó)內(nèi)外學(xué)者對(duì)PMD的研究主要集中在運(yùn)用解析法和數(shù)值法對(duì)其電磁場(chǎng)進(jìn)行建模、計(jì)算以及性能分析等方面,并且取得了大量的研究成果[2-5]。對(duì)PMD結(jié)構(gòu)參數(shù)的優(yōu)化方面研究卻不多,如文獻(xiàn)[6]、文獻(xiàn)[7]和文獻(xiàn)[8]側(cè)重于討論不同優(yōu)化方法針對(duì)PMD的單一性能的優(yōu)化效果,在這些文獻(xiàn)中性能指標(biāo)都是通過(guò)解析法獲得,計(jì)算過(guò)程復(fù)雜且不準(zhǔn)確。代理模型理論的出現(xiàn)解決了上述問(wèn)題,且被認(rèn)為是解決復(fù)雜工程設(shè)計(jì)優(yōu)化問(wèn)題的有效途徑[9]。目前,已經(jīng)有研究者嘗試將單一代理模型引入到PMD的設(shè)計(jì)中,如文獻(xiàn)[10]和文獻(xiàn)[11]分別引入響應(yīng)曲面(response surface method,RSM)和支持向量機(jī)(support vector machine regression,SVR)對(duì)PMD的偏心磁極結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)。

    文中以盤式PMD為研究對(duì)象,著重關(guān)注其主體結(jié)構(gòu)的參數(shù),提出采用基于交叉驗(yàn)證誤差的最優(yōu)加權(quán)法構(gòu)建其混合代理模型,將RSM、Kriging和SVR方法結(jié)合起來(lái),按照交叉驗(yàn)證誤差平方和最小的原則進(jìn)行權(quán)值計(jì)算,求得最優(yōu)加權(quán)系數(shù),建立了PMD最優(yōu)加權(quán)組合模型,即優(yōu)化設(shè)計(jì)的指標(biāo)函數(shù);同時(shí),引入了一種改進(jìn)的熵權(quán)法,將重量和渦流損耗指標(biāo)轉(zhuǎn)化為單一的綜合指標(biāo),并利用粒子群優(yōu)化方法求取最佳結(jié)構(gòu)參數(shù)。通過(guò)ANSYS仿真分析和測(cè)試平臺(tái)驗(yàn)證,表明文中提出的基于熵權(quán)與混合代理模型的方法較好的實(shí)現(xiàn)了PMD的優(yōu)化設(shè)計(jì)目標(biāo)。

    1混合代理模型理論

    單個(gè)代理模型預(yù)測(cè)能力的差異性給工程人員的選取帶來(lái)了一定的風(fēng)險(xiǎn)與挑戰(zhàn),而基于多個(gè)代理模型構(gòu)造混合代理模型的方法可以降低這些風(fēng)險(xiǎn),同時(shí)提高模型的精確度,將其表示為

    (1)

    式中:yh為混合代理模型的響應(yīng)預(yù)測(cè)值,m是單個(gè)代理模型的個(gè)數(shù),μi和yi分別為第i個(gè)模型的權(quán)系數(shù)和響應(yīng)預(yù)測(cè)值,設(shè)計(jì)變量x為矢量,式中權(quán)重系數(shù)之和等于1。因此,單個(gè)模型的選取以及權(quán)重系數(shù)的計(jì)算是影響混合代理模型性能的重要因素。

    1.1單個(gè)代理模型的選取

    混合代理模型中單個(gè)模型的選取應(yīng)當(dāng)滿足互補(bǔ)原則,由于多項(xiàng)式響應(yīng)面模型適用于低維數(shù)據(jù),Kriging模型適用于高維數(shù)據(jù),而SVR適用于小樣本數(shù)據(jù),因此,將這3種單代理模型結(jié)合起來(lái)充分的利用它們各自的優(yōu)勢(shì)來(lái)構(gòu)造混合代理模型。

    1)多項(xiàng)式響應(yīng)面模型。多項(xiàng)式響應(yīng)面是目前廣泛應(yīng)用的一種代理模型,在工程中已經(jīng)使用30多年。選擇常用的具有較高準(zhǔn)確性的二階多項(xiàng)式響應(yīng)面模型,則n個(gè)設(shè)計(jì)變量的二階多項(xiàng)式響應(yīng)面模型可表示為[12]

    (2)

    式中:Y為響應(yīng)變量;β為回歸系數(shù),其元素個(gè)數(shù)為(n+1)(n+2)/2;ξ為一隨機(jī)誤差向量。

    2)Kriging模型。Kriging模型采用了插值方法,它擬合的曲面通過(guò)所有的樣本點(diǎn),是一種估計(jì)方差最小的無(wú)偏估計(jì)模型,其響應(yīng)預(yù)測(cè)值可以表示為[13]

    (3)

    (4)

    β=(ITM-1I)-1(ITM-1I)。

    (5)

    式中:M為相關(guān)系數(shù)矩陣,θk(θk≥0)為相關(guān)系數(shù)待估參數(shù)向量θ的第k個(gè)元素,r(X)為n×1相關(guān)系數(shù)向量,且ri(X)=R(X,Xi),i=1,…,m。

    3)SVR模型。SVR是支持向量機(jī)在回歸領(lǐng)域的應(yīng)用,是基于小樣本統(tǒng)計(jì)學(xué)習(xí)理論的一種通用機(jī)器學(xué)習(xí)方法。采用分類近似支持向量回歸機(jī)模型,如下所示[14]:

    (6)

    (7)

    式中:αi為拉格朗日系數(shù),b為值域,m為樣本個(gè)數(shù),zi為類別標(biāo)簽,取值為+1或者-1,ε為容許誤差。核函數(shù)選擇徑向核函數(shù)(RBF),如式(7)所示,核函數(shù)相關(guān)參數(shù)的確定和容量誤差的選擇采用交叉驗(yàn)證法。

    1.2權(quán)重系數(shù)計(jì)算

    權(quán)重系數(shù)μi是構(gòu)造混合代理模型的關(guān)鍵,提出了一種基于交叉驗(yàn)證誤差的最優(yōu)加權(quán)法,該方法不需要額外的樣本信息,整個(gè)計(jì)算過(guò)程可以通過(guò)求解二次規(guī)劃的最優(yōu)問(wèn)題得到[15]:

    采用LOO(Leave one out)法計(jì)算第i個(gè)代理模型在樣本j的預(yù)測(cè)誤差為

    (8)

    對(duì)應(yīng)的誤差矩陣為

    E=[(eij)m×l][(eij)m×l]T。

    (9)

    得到混合代理模型的預(yù)測(cè)誤差及其平方和分別為

    (10)

    (11)

    以交叉驗(yàn)證誤差平方和最小為準(zhǔn)則,通過(guò)求解下面的模型來(lái)選擇混合權(quán)值,可以表示為

    (12)

    引入拉格朗日乘子,求解二次規(guī)劃問(wèn)題,得到

    (13)

    2熵值權(quán)重法理論

    對(duì)PMD的結(jié)構(gòu)參數(shù)進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì)前,需要確定產(chǎn)品的綜合評(píng)價(jià)指標(biāo),即將多指標(biāo)轉(zhuǎn)化為單一的綜合指標(biāo),在這個(gè)過(guò)程中,權(quán)重系數(shù)的確定是關(guān)鍵環(huán)節(jié),直接影響最終結(jié)構(gòu)參數(shù)的合理性和準(zhǔn)確性。引入了一種改進(jìn)的熵值權(quán)重法確定客觀權(quán)重,該方法修正了傳統(tǒng)熵權(quán)法在所有熵值都趨近于1時(shí)熵權(quán)與熵值傳遞的信息不一致的問(wèn)題[16],步驟如下:

    1)構(gòu)造指標(biāo)初始矩陣A。設(shè)有k個(gè)評(píng)價(jià)指標(biāo)和l個(gè)樣本,則指標(biāo)初始矩陣為

    (14)

    2)初始矩陣的標(biāo)準(zhǔn)化。為消除各指標(biāo)量綱、數(shù)量級(jí)不同對(duì)評(píng)價(jià)結(jié)果的影響,對(duì)初始指標(biāo)矩陣標(biāo)準(zhǔn)化處理,得到矩陣B。

    3)計(jì)算指標(biāo)信息熵值和信息效用值。第j項(xiàng)指標(biāo)的信息熵和信息效用值分別如式(15)和式(16)所示。

    (15)

    dj=1-ej,

    (16)

    (17)

    (18)

    (19)

    3永磁驅(qū)動(dòng)器熵權(quán)與混合代理模型的建立

    單組盤式結(jié)構(gòu)的PMD如圖1所示,主要包括導(dǎo)體轉(zhuǎn)子、永磁轉(zhuǎn)子以及絲杠調(diào)節(jié)裝置,各部分具體功能以及工作原理,可以參考相關(guān)文獻(xiàn)[17],在此不再贅述。結(jié)合眾多文獻(xiàn)關(guān)于PMD結(jié)構(gòu)參數(shù)對(duì)其性能影響的闡述,選擇PMD的主體結(jié)構(gòu)參數(shù)(永磁體長(zhǎng)x1、永磁體寬x2、永磁體高x3,銅盤厚度x4)為設(shè)計(jì)參數(shù),輸出轉(zhuǎn)矩、渦流損耗和重量為優(yōu)化設(shè)計(jì)的指標(biāo)。利用熵權(quán)與混合代理模型法對(duì)PMD結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì)可以歸納為如下過(guò)程:

    圖1 單組盤式結(jié)構(gòu)的示意圖Fig.1 Schematic of single group disk structure

    1)建立PMD的主體有限元模型。采用ANSYS軟件建立PMD的三維有限元模型[17],如圖2所示。為求接近設(shè)備運(yùn)行的實(shí)際情況,設(shè)定原動(dòng)機(jī)為三相異步電動(dòng)機(jī),其額定功率為15 Kw,額定轉(zhuǎn)速為1 455 r/min,表1為模型固定不變的模型參數(shù)。

    圖2 PMD的三維有限元仿真模型Fig.2 3D finite element simulation model of PMD

    參數(shù)數(shù)值參數(shù)數(shù)值導(dǎo)體材料(黃銅)H62永磁材料(釹鐵硼)N38H永磁體數(shù)量18氣隙長(zhǎng)度/mm5永磁體矯頑力(A/n)87500鋼盤內(nèi)徑/mm90鋼盤外徑/mm140銅電阻率/(ohm/m)7.1E-08永磁體相對(duì)磁導(dǎo)率1.05鋼盤厚度/mm10

    2)獲取樣本數(shù)據(jù)。正交試驗(yàn)設(shè)計(jì)法是一種在工程中被廣泛應(yīng)用的獲取樣本點(diǎn)的實(shí)驗(yàn)方法,該方法獲得的數(shù)據(jù)分布均勻,結(jié)論可靠,因此引入該方法選取50組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行后續(xù)代理模型的建立,表2是實(shí)驗(yàn)中采用的因素-水平表。

    表2 結(jié)構(gòu)因素水平表

    3)建立混合代理模型。根據(jù)前述代理模型理論,分別建立PMD輸出轉(zhuǎn)矩y1、渦流損耗y2和重量y3關(guān)于優(yōu)化設(shè)計(jì)參數(shù)的三種代理模型及其混合代理模型。通過(guò)最優(yōu)加權(quán)法得到混合權(quán)重系數(shù)為

    W=(w1,w2,w3)T=

    (0.2822,0.4475,0.2703)T。

    (20)

    利用均方根誤差(RMSE)作為代理模型的一個(gè)誤差分析指標(biāo),其計(jì)算公式如式(21)所示。

    (21)

    如表3所示為RSM模型、Kriging模型、SVR模型以及混合代理模型關(guān)于PMD的輸出轉(zhuǎn)矩、渦流損耗以及重量的統(tǒng)計(jì)誤差RMSE。由表3可以看出,混合代理模型對(duì)指標(biāo)函數(shù)的擬合精度要高于單一代理模型,因此,用來(lái)構(gòu)造PMD優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型更有優(yōu)勢(shì)。

    表3不同方法的統(tǒng)計(jì)參數(shù)比較

    Table 3Statistical parameters(RMSE)of different

    models

    響應(yīng)量RSMKrigingSVRHybrid輸出轉(zhuǎn)矩1.245381.055191.187210.65346渦流損耗1.078650.975861.144170.57421重量0.653160.745380.639910.45355

    4)計(jì)算PMD設(shè)計(jì)的綜合指標(biāo)。根據(jù)前述改進(jìn)的熵值權(quán)重理論方法得到綜合指標(biāo)的權(quán)重系數(shù)值,最終得到如下式所示的綜合指標(biāo)函數(shù)

    f(x)=δ1f2(x)+δ2f3(x)=0.643f2(x)+0.357f3(x)。

    (22)

    5)建立優(yōu)化設(shè)計(jì)的最終數(shù)學(xué)模型。構(gòu)建精確完備的約束條件是建立優(yōu)化對(duì)象數(shù)學(xué)模型的關(guān)鍵因素,PMD的約束條件主要有結(jié)構(gòu)約束,如式(24)所示以及性能約束[8],如式(25)所示。則整個(gè)PMD數(shù)學(xué)模型可以描述為:

    miny=f(x),

    (23)

    (24)

    (25)

    式(24)是PMD的結(jié)構(gòu)約束條件,根據(jù)實(shí)際的工程需得到。式(25)的第一個(gè)式子表示優(yōu)化之后PMD的啟動(dòng)轉(zhuǎn)矩要大于標(biāo)準(zhǔn)模型的啟動(dòng)轉(zhuǎn)矩,Tsmin表示待設(shè)計(jì)PMD的啟動(dòng)轉(zhuǎn)矩,假定該值為25 N·m,第二個(gè)式子表示優(yōu)化設(shè)計(jì)的PMD的輸出最大轉(zhuǎn)矩要大于給定的最大轉(zhuǎn)矩Tmax,假定該值為90.15 N·m,誤差范圍在1%左右認(rèn)為合理,否則認(rèn)為設(shè)計(jì)失敗。

    6)選用合適的優(yōu)化方法進(jìn)行優(yōu)化計(jì)算。利用自適應(yīng)權(quán)重PSO算法[10]對(duì)步驟5)中得到的數(shù)學(xué)模型進(jìn)行求解。該算法的優(yōu)點(diǎn)在于:無(wú)論在搜索初期還是后期都具有較好的局部和全局搜索能力,同時(shí)由于加速因子隨迭代次數(shù)變化,會(huì)加快粒子在搜索后期的收斂速度。算法當(dāng)中加速因子的基數(shù)α0=0.5,慣性權(quán)重的常數(shù)ω0=0.5,最大迭代次數(shù)kmax=100,種群規(guī)模N=50。

    4仿真分析與驗(yàn)證

    根據(jù)文獻(xiàn)[10]和文獻(xiàn)[11]中的PMD多目標(biāo)優(yōu)化設(shè)計(jì)方法,得到不同的優(yōu)化設(shè)計(jì)結(jié)果,與文中方法得出的結(jié)果作對(duì)比,如表4所示。

    表4 不同方法得到的結(jié)構(gòu)參數(shù)比較

    由表4可以看出,此方法得到的PMD性能指標(biāo)較其它兩種方法都得到了提高,而且經(jīng)過(guò)計(jì)算可知,與其它方法的最優(yōu)指標(biāo)相比,此方法使得PMD的渦流損耗降低約5%,同時(shí)輸出轉(zhuǎn)矩提高了約6%且PMD的重量降低了約10%,因此,很顯然優(yōu)化設(shè)計(jì)效果要好于其它兩種方法。

    使用三維有限元仿真軟件ANSYS仿真,得到3種結(jié)構(gòu)參數(shù)下PMD的性能擬合曲線,如圖3所示。從圖中可以看出,3種方法均滿足前述性能約束條件。相比較而言,此方法對(duì)PMD的性能最好,使得PMD在減少重量的情況下(意味著減少了材料使用),提高了最大輸出轉(zhuǎn)矩,同時(shí)啟動(dòng)轉(zhuǎn)矩幾乎沒(méi)有受到影響。

    銅盤的渦流特性在PMD中特別重要,一是因?yàn)榕ぞ氐膫鬟f需要銅盤的渦流產(chǎn)生磁場(chǎng),二是PMD的損耗主要以渦流產(chǎn)生熱量的形式從銅盤上釋放出來(lái)。圖5所示為此設(shè)計(jì)的PMD中銅盤表面渦流密度的矢量圖,可以看出,渦流分布是較為均勻合理,主要集中在銅盤表面(面對(duì)永磁體側(cè))的薄層中,集膚效應(yīng)明顯且渦流在銅盤表面形成了18個(gè)回路,這樣有利于減少損耗。

    圖3 優(yōu)化之后PMD的性能Fig.3 Performance of optimized PMD

    圖4 銅盤表面渦流密度的矢量圖Fig.4    Vector diagram of eddy current density    on copper disk surface

    如果結(jié)構(gòu)設(shè)計(jì)不合理,運(yùn)行時(shí)PMD的永磁體會(huì)逐漸發(fā)熱升溫,使其性能降低,因此有必要分析PMD的熱量分布情況,如圖6所示為此設(shè)計(jì)的PMD銅盤表面熱量分布的等值圖??梢钥闯觯瑹崃糠植寂c渦流分布的規(guī)律類似,熱量主要分布在銅盤表面薄層中并且沿軸向由內(nèi)至外逐漸減??;熱量與渦流回路相對(duì)應(yīng),主要集中在相鄰回路的徑向路徑上,這樣有利于散熱。

    動(dòng)態(tài)時(shí)磁路中除了永磁場(chǎng)外,還存在渦流感應(yīng)磁場(chǎng),如圖6所示為此設(shè)計(jì)的PMD在動(dòng)態(tài)時(shí)的磁通密度矢量圖,可以看出,磁通從永磁體的N極出發(fā),經(jīng)氣隙、銅盤表面的感應(yīng)磁場(chǎng),到達(dá)相鄰永磁體的S極,最后通過(guò)永磁體的背襯鋼盤構(gòu)成閉合的主磁路;與靜態(tài)時(shí)相比,使得PMD的磁通密度得到了大幅度的提高,有利于產(chǎn)生更高的傳動(dòng)轉(zhuǎn)矩。

    圖5 銅盤表面熱量分布的等值圖Fig.5    Contour diagram of heat distribution    on copper disk surface

    圖6 動(dòng)態(tài)時(shí)PMD的磁通密度矢量圖Fig.6    Vector diagram of dynamic magnetic flux density    on adjustable permanent magnetic coupler

    圖7 PMD樣機(jī)示意圖Fig.7 Schematic of PMD model machine

    圖8 PMD實(shí)驗(yàn)平臺(tái)的結(jié)構(gòu)圖Fig.8 Structure of PMD test bench

    為了進(jìn)一步驗(yàn)證文中所設(shè)計(jì)PMD的真實(shí)性能狀況,按照上面得到的結(jié)構(gòu)參數(shù)制造了如圖7所示的簡(jiǎn)易樣機(jī),同時(shí),按照?qǐng)D8所示的結(jié)構(gòu)搭建PMD實(shí)驗(yàn)室仿真平臺(tái)。通過(guò)實(shí)驗(yàn),測(cè)得該P(yáng)MD的啟動(dòng)轉(zhuǎn)矩為26.8 N·m,最大輸出轉(zhuǎn)矩近似為107.5 N·m,這滿足了式(25)的性能約束條件,與ANSYS仿真的結(jié)果基本吻合,表明此方法是可行的。

    5結(jié)論

    針對(duì)永磁驅(qū)動(dòng)器的結(jié)構(gòu)設(shè)計(jì)問(wèn)題,文中提出了一種基于熵權(quán)與混合代理模型的優(yōu)化設(shè)計(jì)方法。該方法的優(yōu)點(diǎn)在于:基于混合代理模型理論建立了有關(guān)PMD渦流損耗、重量以及輸出轉(zhuǎn)矩的回歸模型,且精確度要高于單一的RSM、Kriging以及SVR模型,避免了用復(fù)雜的解析法求取優(yōu)化設(shè)計(jì)的指標(biāo)函數(shù);通過(guò)引入改進(jìn)的熵權(quán)法,將多指標(biāo)轉(zhuǎn)化為單指標(biāo),建立了基于PMD綜合指標(biāo)函數(shù)的數(shù)學(xué)模型,使其優(yōu)化求解更加方便準(zhǔn)確。仿真分析結(jié)果表明:該方法得到的PMD結(jié)構(gòu)參數(shù)可靠有效,與其它方法相比,優(yōu)化設(shè)計(jì)效果有了明顯改善,為PMD的優(yōu)化設(shè)計(jì)提供了一種可行的方法。

    參 考 文 獻(xiàn):

    [1]ANNA WANG,WANG J,SHI G,et all.Study on nonlinear regression modeling methods of the permanent magnet drive[C]//International Conference on Electronics and Optoelectronics,July 26-29,2011,Dalian,China.2011: 87-90.

    [2]WANG Jian,LIN Heyun,FANG S,et all.A general analytical model of permanent magnet eddy current couplings[J].IEEE Transcations on magnetics,2014,50 (1):1-9.

    [3]JACOB KRIZAN A,SCOTT SUDHOFF D.A design model for salient permanent-magnet machines with investigation of saliency and wide-speed-range performance[J].IEEE Transcations on energy conversion,2013,28(1):95-105.

    [4]S M,MIRSALIM M,VAEZZADEH S,et all.Analytical modeling and analysis of axial-flux interior permanent-magnet couplers[J].IEEE Transcations on industialI electronics,2014,51 (11):5940-5947.

    [5]楊超君,鄭武,李志寶.可調(diào)速異步盤式磁力聯(lián)軸器的轉(zhuǎn)矩計(jì)算及其影響因素分析[J].電機(jī)與控制學(xué)報(bào),2012,16(1):85-91.

    YANG Chaojun,ZHENG Wu,LI Zhibao.Torque calculation of speed-adjustable asynchronous disk type of magnetic coupling and its impact factors analysis[J].Electric Machinesand Control,2012,16(1):85-91.

    [6]SMITH A C,EIWAKEEL A S,WALLACE A.Formal design optimization of PM drive couplings[C]// Proceedings of 37th IAS Annual Meeting IEEE Conference of the Industry Applications,Oct 13-18,2002,New York,USA,2002: 205-211.

    [7]WANG A,WANG J,WU B.Structural optimization of the permanent magnet drive based on artificial neural network and particle swarm optimization[C]//IEEE Conference of the Third Intelligent Human-Machine Systems and Cybernetics,August 26-27,2011,Hangzhou,China 2011:70-75.

    [8]AMGED S,EIWAKEEL.Design optimization of PM couplings using hybrid Particle Swarm Optimization-Simplex Method (PSO-SM) Algorithm[J].Electric Power Systems Research,2014,116:29-35.

    [9]YAN L,MATTHEW C.Improving surrogate-assisted variable fidelity multi-objective optimization using a clustering algorithm[J].Applied Soft Computing 2014,24:482-493.

    [10]石松寧,王大志,時(shí)統(tǒng)宇.基于RSM的永磁驅(qū)動(dòng)器偏心磁極的多目標(biāo)優(yōu)化[J].儀器儀表學(xué)報(bào),2014,35(9):1963-1971.

    SHISongning,WANG Dazhi,SHI Tongyu.Multi-objective optimization of eccentric magnet pole for permanent magnet drive based on response surface methodology[J].Journal of Scientific Instrumentation,2014,35(9):1963-1971.

    [11]石松寧,王大志,時(shí)統(tǒng)宇.永磁驅(qū)動(dòng)器偏心磁極的優(yōu)化設(shè)計(jì)[J].東北大學(xué)學(xué)報(bào):自然科學(xué)版,2014,35(8):1078-1082.

    SHI Songning,WANG Dazhi,SHI Tongyu.Optimal design of eccentric magnetic pole for permanent magnet drive[J].Journal of Northeastern University: Natural Science Edition,2014,35(8):1078-1082.

    [12]VAFAEESEFAT A.Optimization of composite pressure vessels with metal liner by adaptive response surface method[J].Journal of Mechanical Science and Technology,2011,25 (11):2811-1816.

    [13]黃章俊,王成恩.基于Kriging模型的渦輪盤優(yōu)化設(shè)計(jì)方法[J].計(jì)算機(jī)集成制造系統(tǒng),2010,16(5):905-911.

    HUANG Zhangjun,WANG Chengen.Turbine discs optimization design based on Kriging model[J].Computer Integrated Manufacturing Systems,2010,16(5):905-911.

    [14]司利云,林輝.基于支持向量機(jī)的開(kāi)關(guān)磁阻電機(jī)轉(zhuǎn)子位置在線建模[J].電機(jī)與控制學(xué)報(bào),2014,18(4):72-78.

    SI Liyun,LIN Hui.Online modeling based on support vector machine for rotor position signal of switched reluctance motors[J].Electric Machinesand Control,2014,18(4):72-78.

    [15]謝開(kāi)貴,何斌,楊萬(wàn)年.組合預(yù)測(cè)權(quán)系數(shù)的確定[J].預(yù)測(cè),1998,17(7):151-154.

    XIE Kaigui,HE Bin,YANG Wannian.Solving the weights of combination forecasting[J].Forecasting,1998,17(7):151-154.

    [16]歐陽(yáng)森,石怡理.改進(jìn)熵權(quán)法及其在電能質(zhì)量評(píng)估中的應(yīng)用[J].電力系統(tǒng)自動(dòng)化,2013,37(21):156-164.

    OUYANG Sen,SHI Yili.A new improved entropy method and its application in power quality evalution[J].Automation of Electric Power System,2013,37(21):156-164.

    [17]王旭,王大志,劉震,等.永磁調(diào)速器的渦流場(chǎng)分析與性能計(jì)算[J].儀器儀表學(xué)報(bào),2012,33(1):156-160.

    WANG Xu,WANG Dazhi,LIU Zhen,et al.Eddy current field analysis and performance calculations for adjustable permanent magnetic coupler[J].Journal of Scientific Instrumentation,2012,33(1):156-159.

    (編輯:賈志超)

    Optimization design for permanent magnet drive based on entropy-weight and hybrid surrogate model

    LI Zhao,WANG Da-zhi,SHI Tong-yu,ZHENG Di

    (School of Information Science & Engineering, Northeastern University, Shenyang 110004, China)

    Abstract:Permanent magnet drive (PMD) design problems were studied and an optimization design approach based on hybrid surrogate model and improved entropy-weight was put forward.Firstly by combining the response surface method (RSM),Kriging method and support vector machine regression (SVR),a hybrid surrogate model of optimizing parameters and response variables was built,in which the optimal weighting method based on cross validation error is used; Secondly the improved entropy-weight was introduced to transform the multi-indexes of PMD into a single composite one,and the mathematical model was built,which was solved by adaptive weight particle swarm optimization particle swarm optimization(AWPSO); Finally the parameters results were analyzed with finite element method and experimental validation.The results show that the proposed method is superior to other methods,and the structural parameters are reasonable and effective.And the expected results are derived.

    Keywords:permanent magnet drive; hybrid surrogate model; entropy-weight method ; optimal weighting method ; adaptive weight particle swarm optimization particle swarm optimization

    收稿日期:2015-01-31

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(61433004);中國(guó)高校基本科研業(yè)務(wù)費(fèi)(N150403005,L1504011)

    作者簡(jiǎn)介:李召(1987—),男,博士研究生,研究方向?yàn)橛来膨?qū)動(dòng)器的建模和優(yōu)化設(shè)計(jì); 王大志(1963—),男,教授,博士生導(dǎo)師,研究方向?yàn)殡娏ο到y(tǒng)與電力傳動(dòng); 時(shí)統(tǒng)宇(1988—),男,博士研究生,研究方向?yàn)橛来膨?qū)動(dòng)器的建模及設(shè)計(jì);

    通訊作者:李召

    DOI:10.15938/j.emc.2016.06.013

    中圖分類號(hào):TH 133

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1007-449X(2016)06-0102-07

    鄭迪(1988—),女,碩士研究生,研究方向?yàn)橛来膨?qū)動(dòng)器的建模及設(shè)計(jì)。

    猜你喜歡
    熵權(quán)法
    高職機(jī)電專業(yè)學(xué)生數(shù)學(xué)能力的調(diào)查及對(duì)策
    考試周刊(2016年103期)2017-01-23 15:58:59
    大學(xué)周邊健身房滿意度調(diào)查報(bào)告
    基于熵權(quán)法的京津冀區(qū)域信息化協(xié)調(diào)發(fā)展規(guī)律模型及其應(yīng)用
    基于熵權(quán)法的“互聯(lián)網(wǎng)+”農(nóng)業(yè)發(fā)展影響因素權(quán)重確定
    商(2016年34期)2016-11-24 20:05:32
    基于熵權(quán)法的西安市外向型經(jīng)濟(jì)發(fā)展綜合評(píng)價(jià)研究
    城市與區(qū)域經(jīng)濟(jì)可持續(xù)發(fā)展
    商(2016年27期)2016-10-17 04:35:55
    蘭州市城關(guān)區(qū)銀行顧客滿意度調(diào)查研究
    商(2016年17期)2016-06-06 08:38:10
    福建省高新技術(shù)產(chǎn)業(yè)競(jìng)爭(zhēng)力評(píng)價(jià)研究
    “互聯(lián)網(wǎng)+”時(shí)代的出租車資源配置研究
    商(2016年12期)2016-05-09 10:19:39
    基于熵權(quán)TOPSIS模型的出版企業(yè)核心競(jìng)爭(zhēng)力評(píng)價(jià)
    x7x7x7水蜜桃| 欧美丝袜亚洲另类 | xxxwww97欧美| 两人在一起打扑克的视频| x7x7x7水蜜桃| 成人鲁丝片一二三区免费| 午夜福利成人在线免费观看| 亚洲不卡免费看| 亚洲第一电影网av| av黄色大香蕉| 国产av麻豆久久久久久久| 免费av观看视频| 国产一级毛片七仙女欲春2| 成年女人看的毛片在线观看| 色av中文字幕| 校园春色视频在线观看| 国产中年淑女户外野战色| 精品欧美国产一区二区三| 免费观看的影片在线观看| 久久久久久久精品吃奶| 免费看美女性在线毛片视频| 少妇裸体淫交视频免费看高清| 免费搜索国产男女视频| 一本久久中文字幕| 欧美乱妇无乱码| 国产美女午夜福利| 日本黄色片子视频| 亚洲av电影不卡..在线观看| 制服人妻中文乱码| 国内精品久久久久精免费| 欧美最黄视频在线播放免费| 特级一级黄色大片| 麻豆成人av在线观看| 国产高清视频在线播放一区| 90打野战视频偷拍视频| 两人在一起打扑克的视频| 亚洲国产欧美网| 97碰自拍视频| 制服丝袜大香蕉在线| e午夜精品久久久久久久| 99久久99久久久精品蜜桃| 精品电影一区二区在线| 免费大片18禁| 99热这里只有精品一区| 久久国产精品影院| 在线免费观看的www视频| 中国美女看黄片| 99热这里只有精品一区| 亚洲av五月六月丁香网| 高清日韩中文字幕在线| 亚洲五月婷婷丁香| 可以在线观看毛片的网站| 国产一区二区三区视频了| 久久中文看片网| 18禁裸乳无遮挡免费网站照片| 欧美国产日韩亚洲一区| 999久久久精品免费观看国产| 国产精品,欧美在线| 午夜福利在线观看吧| 日本在线视频免费播放| 精品国产超薄肉色丝袜足j| 男女那种视频在线观看| 波多野结衣巨乳人妻| 亚洲 欧美 日韩 在线 免费| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全电影3| 久久久国产精品麻豆| 波多野结衣高清作品| 国产成人福利小说| 搡老妇女老女人老熟妇| 黄片大片在线免费观看| 可以在线观看的亚洲视频| 欧美3d第一页| 国产高清videossex| 国产av麻豆久久久久久久| 亚洲国产精品合色在线| 国产免费一级a男人的天堂| 国产男靠女视频免费网站| www.色视频.com| 极品教师在线免费播放| 三级毛片av免费| 在线国产一区二区在线| 美女大奶头视频| 亚洲av第一区精品v没综合| 老司机福利观看| 欧美日韩瑟瑟在线播放| 色吧在线观看| 一本一本综合久久| 欧美高清成人免费视频www| 国产男靠女视频免费网站| 一个人免费在线观看电影| 成年女人毛片免费观看观看9| 老司机深夜福利视频在线观看| 美女大奶头视频| 欧美在线一区亚洲| 最新美女视频免费是黄的| 国产精品亚洲一级av第二区| 十八禁网站免费在线| 最近在线观看免费完整版| 国产av麻豆久久久久久久| 国产精品久久久久久亚洲av鲁大| 亚洲av一区综合| 免费看日本二区| 岛国在线观看网站| 国产欧美日韩一区二区三| 国产成人影院久久av| 久久久精品大字幕| 久久欧美精品欧美久久欧美| 国内毛片毛片毛片毛片毛片| 国产精品永久免费网站| 欧美日韩综合久久久久久 | 成人精品一区二区免费| 亚洲成人精品中文字幕电影| 欧美zozozo另类| 精品久久久久久成人av| 动漫黄色视频在线观看| 久久久久久大精品| 狂野欧美激情性xxxx| 内射极品少妇av片p| 国产一区二区在线观看日韩 | 18+在线观看网站| 婷婷精品国产亚洲av在线| 最近视频中文字幕2019在线8| 在线看三级毛片| 又黄又粗又硬又大视频| 高清日韩中文字幕在线| 精品一区二区三区视频在线 | 日韩av在线大香蕉| 国产色爽女视频免费观看| 日本黄大片高清| 久久久国产精品麻豆| 久久精品国产99精品国产亚洲性色| 又黄又爽又免费观看的视频| 十八禁网站免费在线| 精品无人区乱码1区二区| 国产伦一二天堂av在线观看| 国产免费av片在线观看野外av| 男人的好看免费观看在线视频| 午夜精品一区二区三区免费看| 97碰自拍视频| av视频在线观看入口| 在线播放国产精品三级| 婷婷精品国产亚洲av在线| av视频在线观看入口| 亚洲天堂国产精品一区在线| 99久久久亚洲精品蜜臀av| 露出奶头的视频| 成熟少妇高潮喷水视频| 1000部很黄的大片| 日本a在线网址| 成人18禁在线播放| 大型黄色视频在线免费观看| 亚洲欧美一区二区三区黑人| 亚洲人成网站高清观看| 国产三级在线视频| www.www免费av| 免费在线观看日本一区| 一级a爱片免费观看的视频| 午夜老司机福利剧场| 怎么达到女性高潮| 国产不卡一卡二| 性色avwww在线观看| 中文在线观看免费www的网站| 欧美日韩精品网址| 最新中文字幕久久久久| 亚洲性夜色夜夜综合| 99热6这里只有精品| 久久精品国产亚洲av涩爱 | 在线观看午夜福利视频| 午夜日韩欧美国产| 久久精品国产清高在天天线| 色视频www国产| 免费无遮挡裸体视频| 欧美黄色片欧美黄色片| 欧美一级毛片孕妇| 三级国产精品欧美在线观看| 两个人视频免费观看高清| 真实男女啪啪啪动态图| 色在线成人网| 男人和女人高潮做爰伦理| 人人妻人人看人人澡| 中文字幕人成人乱码亚洲影| 午夜福利视频1000在线观看| 在线播放无遮挡| 国产中年淑女户外野战色| 日韩大尺度精品在线看网址| 3wmmmm亚洲av在线观看| 午夜精品在线福利| 亚洲天堂国产精品一区在线| 每晚都被弄得嗷嗷叫到高潮| 夜夜爽天天搞| 91在线观看av| 日本精品一区二区三区蜜桃| 色av中文字幕| 亚洲真实伦在线观看| 夜夜爽天天搞| 亚洲国产日韩欧美精品在线观看 | 99热6这里只有精品| 99热这里只有精品一区| 99久久99久久久精品蜜桃| 国产一级毛片七仙女欲春2| 亚洲av成人av| 国产成人欧美在线观看| 一进一出抽搐gif免费好疼| 搡老妇女老女人老熟妇| 精品久久久久久久人妻蜜臀av| 少妇丰满av| 成人性生交大片免费视频hd| 小说图片视频综合网站| av视频在线观看入口| 日日夜夜操网爽| 亚洲精品成人久久久久久| 久久九九热精品免费| 欧美+日韩+精品| 久久中文看片网| 成人鲁丝片一二三区免费| 国产午夜福利久久久久久| 日韩欧美免费精品| 国产精品电影一区二区三区| 亚洲片人在线观看| 国产视频内射| 亚洲成人免费电影在线观看| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久久久免费视频| 色噜噜av男人的天堂激情| 精品一区二区三区视频在线 | 国产精品嫩草影院av在线观看 | 亚洲国产中文字幕在线视频| 午夜影院日韩av| 久久午夜亚洲精品久久| 国语自产精品视频在线第100页| 国产三级黄色录像| 天堂影院成人在线观看| 国产av在哪里看| 日韩高清综合在线| 757午夜福利合集在线观看| 免费高清视频大片| 日本在线视频免费播放| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 热99re8久久精品国产| 又紧又爽又黄一区二区| 日日干狠狠操夜夜爽| 99精品在免费线老司机午夜| 午夜a级毛片| 中文字幕精品亚洲无线码一区| 欧美最新免费一区二区三区 | 嫩草影院入口| 国产午夜精品论理片| 深夜精品福利| 在线播放国产精品三级| 在线观看一区二区三区| 国产精品亚洲av一区麻豆| 999久久久精品免费观看国产| av视频在线观看入口| 中文字幕av在线有码专区| 国产免费av片在线观看野外av| 色视频www国产| 一区二区三区高清视频在线| 最近手机中文字幕大全| h日本视频在线播放| 99九九线精品视频在线观看视频| 少妇的逼水好多| 亚洲综合色惰| 熟妇人妻不卡中文字幕| 少妇人妻一区二区三区视频| 国产乱来视频区| 成人毛片60女人毛片免费| 91久久精品国产一区二区成人| 波多野结衣巨乳人妻| 国产精品国产三级国产av玫瑰| 女的被弄到高潮叫床怎么办| 2021天堂中文幕一二区在线观| 日本欧美国产在线视频| 欧美区成人在线视频| 日本免费a在线| 成人一区二区视频在线观看| 国产男人的电影天堂91| 国产成人a区在线观看| 国产不卡一卡二| 非洲黑人性xxxx精品又粗又长| 精品国产露脸久久av麻豆 | 中文资源天堂在线| 久久久久久久午夜电影| 秋霞伦理黄片| 国产亚洲91精品色在线| 国产亚洲av片在线观看秒播厂 | 少妇熟女欧美另类| 亚洲精品乱码久久久久久按摩| 男插女下体视频免费在线播放| 日韩,欧美,国产一区二区三区| 国产黄a三级三级三级人| 久久精品国产鲁丝片午夜精品| 亚洲成人精品中文字幕电影| 美女高潮的动态| 久久精品人妻少妇| 国产精品美女特级片免费视频播放器| 99视频精品全部免费 在线| av在线播放精品| 久久久欧美国产精品| 久久久国产一区二区| 国产大屁股一区二区在线视频| 免费播放大片免费观看视频在线观看| 亚洲一区高清亚洲精品| 久久99精品国语久久久| 国产精品一二三区在线看| 精品人妻偷拍中文字幕| 好男人在线观看高清免费视频| 人体艺术视频欧美日本| 午夜老司机福利剧场| 伦理电影大哥的女人| 欧美一区二区亚洲| 中文资源天堂在线| 男女国产视频网站| 久久久欧美国产精品| 2018国产大陆天天弄谢| 久久精品熟女亚洲av麻豆精品 | 亚洲av在线观看美女高潮| 亚洲成人中文字幕在线播放| 女的被弄到高潮叫床怎么办| 日韩中字成人| 欧美成人精品欧美一级黄| 2022亚洲国产成人精品| 国产视频内射| 日韩欧美精品v在线| 免费在线观看成人毛片| 只有这里有精品99| 啦啦啦啦在线视频资源| 国产真实伦视频高清在线观看| 国产中年淑女户外野战色| 国内少妇人妻偷人精品xxx网站| 春色校园在线视频观看| 欧美xxⅹ黑人| 国产片特级美女逼逼视频| 亚洲国产高清在线一区二区三| 欧美xxxx黑人xx丫x性爽| 岛国毛片在线播放| 在线观看免费高清a一片| 看免费成人av毛片| av在线播放精品| 十八禁网站网址无遮挡 | 你懂的网址亚洲精品在线观看| 日韩av在线大香蕉| 黄色日韩在线| 午夜福利在线观看吧| 极品少妇高潮喷水抽搐| 日韩不卡一区二区三区视频在线| 色尼玛亚洲综合影院| 3wmmmm亚洲av在线观看| 丝瓜视频免费看黄片| 日韩中字成人| 国产久久久一区二区三区| 黄片wwwwww| 亚洲国产成人一精品久久久| 日韩不卡一区二区三区视频在线| 高清在线视频一区二区三区| 97超视频在线观看视频| 99热这里只有是精品50| 亚洲国产欧美在线一区| 日韩电影二区| 成人特级av手机在线观看| 欧美变态另类bdsm刘玥| av在线天堂中文字幕| 国产精品人妻久久久久久| 国产视频首页在线观看| 国产午夜精品论理片| 国产精品女同一区二区软件| 久久久久久久久久成人| 久久99热这里只频精品6学生| 亚洲在久久综合| 国产黄片视频在线免费观看| 国产白丝娇喘喷水9色精品| 亚洲aⅴ乱码一区二区在线播放| 国产成人freesex在线| 最近中文字幕高清免费大全6| 国产成人freesex在线| 免费观看性生交大片5| 搞女人的毛片| 一个人观看的视频www高清免费观看| 亚洲,欧美,日韩| 偷拍熟女少妇极品色| 国产精品不卡视频一区二区| 超碰av人人做人人爽久久| 亚洲在久久综合| 久久99热这里只频精品6学生| 国产一区二区三区综合在线观看 | 偷拍熟女少妇极品色| 伦精品一区二区三区| 日韩成人av中文字幕在线观看| 国产伦在线观看视频一区| 免费看美女性在线毛片视频| 岛国毛片在线播放| 午夜老司机福利剧场| 国产有黄有色有爽视频| 国产精品熟女久久久久浪| 两个人的视频大全免费| 欧美xxxx性猛交bbbb| 亚洲精品乱码久久久久久按摩| 免费在线观看成人毛片| 精品久久久久久久久亚洲| 免费在线观看成人毛片| 麻豆久久精品国产亚洲av| 欧美激情在线99| 国产伦精品一区二区三区四那| av国产免费在线观看| 日本欧美国产在线视频| 黄色一级大片看看| 夫妻午夜视频| 欧美成人午夜免费资源| 国产乱人视频| 69人妻影院| 成人一区二区视频在线观看| 搡女人真爽免费视频火全软件| 一本一本综合久久| 国产探花极品一区二区| xxx大片免费视频| 亚洲av在线观看美女高潮| 国产精品福利在线免费观看| 亚洲av免费高清在线观看| 麻豆成人av视频| 亚洲av福利一区| 99久久人妻综合| 精品人妻视频免费看| 国产精品一区www在线观看| 国产精品久久久久久久电影| 色播亚洲综合网| 国产成人精品福利久久| 亚洲欧美清纯卡通| 久久热精品热| 蜜桃亚洲精品一区二区三区| 婷婷六月久久综合丁香| 日韩视频在线欧美| 成人漫画全彩无遮挡| 亚洲av中文字字幕乱码综合| 亚洲熟女精品中文字幕| 日韩人妻高清精品专区| 国产一区二区亚洲精品在线观看| 亚洲一区高清亚洲精品| 午夜福利在线在线| 日韩 亚洲 欧美在线| 亚洲av日韩在线播放| 五月天丁香电影| 成年女人看的毛片在线观看| 国产av国产精品国产| 中文资源天堂在线| 精品久久国产蜜桃| 99久国产av精品| 99视频精品全部免费 在线| 免费观看无遮挡的男女| 人人妻人人看人人澡| 亚洲精品国产成人久久av| 亚洲av不卡在线观看| 在线观看av片永久免费下载| 少妇猛男粗大的猛烈进出视频 | 国产伦精品一区二区三区四那| 婷婷色综合大香蕉| 三级毛片av免费| 大话2 男鬼变身卡| 久久韩国三级中文字幕| av在线亚洲专区| 2018国产大陆天天弄谢| 观看免费一级毛片| 亚洲美女视频黄频| 蜜桃久久精品国产亚洲av| 日本与韩国留学比较| 超碰97精品在线观看| 亚洲精品成人久久久久久| 女人被狂操c到高潮| 全区人妻精品视频| 在线观看免费高清a一片| 国产精品一区二区三区四区免费观看| 在线免费观看的www视频| 国产v大片淫在线免费观看| 国产免费福利视频在线观看| 亚洲最大成人手机在线| 夫妻性生交免费视频一级片| 一二三四中文在线观看免费高清| 国产成人a区在线观看| 美女大奶头视频| 久久精品国产亚洲网站| 人人妻人人澡欧美一区二区| 国产精品国产三级国产专区5o| 久久精品人妻少妇| 成年人午夜在线观看视频 | 亚洲国产欧美人成| 久久精品熟女亚洲av麻豆精品 | 干丝袜人妻中文字幕| 精品国产露脸久久av麻豆 | 亚洲国产日韩欧美精品在线观看| 深爱激情五月婷婷| 免费看光身美女| 69人妻影院| 精品酒店卫生间| 亚洲高清免费不卡视频| 久久久久久久久久人人人人人人| 丰满少妇做爰视频| 免费观看av网站的网址| 男女视频在线观看网站免费| 日本熟妇午夜| 超碰av人人做人人爽久久| a级一级毛片免费在线观看| 久久久精品94久久精品| 有码 亚洲区| 精品人妻一区二区三区麻豆| 亚洲精品视频女| 亚洲国产最新在线播放| 色5月婷婷丁香| 丰满乱子伦码专区| 国产色爽女视频免费观看| 日韩欧美三级三区| 特大巨黑吊av在线直播| 久久精品国产自在天天线| 亚洲人成网站在线播| 国产成人免费观看mmmm| 精品一区二区三卡| 少妇的逼好多水| 午夜免费男女啪啪视频观看| 99九九线精品视频在线观看视频| 欧美激情国产日韩精品一区| 日日摸夜夜添夜夜爱| 人人妻人人澡人人爽人人夜夜 | 老女人水多毛片| 亚洲熟妇中文字幕五十中出| 国产欧美日韩精品一区二区| 男女国产视频网站| 中国美白少妇内射xxxbb| 日本爱情动作片www.在线观看| 女人久久www免费人成看片| av卡一久久| 人人妻人人澡人人爽人人夜夜 | 蜜桃久久精品国产亚洲av| 久久精品久久久久久噜噜老黄| av黄色大香蕉| 日本猛色少妇xxxxx猛交久久| 又黄又爽又刺激的免费视频.| 高清av免费在线| 国产色婷婷99| 色哟哟·www| 国产在线男女| 欧美一区二区亚洲| 久久久久久久久久久免费av| 卡戴珊不雅视频在线播放| 伊人久久精品亚洲午夜| 国产av国产精品国产| 亚洲成人av在线免费| 精品欧美国产一区二区三| 91久久精品电影网| 国产精品一二三区在线看| 久久国产乱子免费精品| 干丝袜人妻中文字幕| 日韩制服骚丝袜av| 日本一本二区三区精品| 国国产精品蜜臀av免费| 伦精品一区二区三区| 亚洲婷婷狠狠爱综合网| 成人亚洲精品一区在线观看 | 亚洲性久久影院| 一区二区三区四区激情视频| 看十八女毛片水多多多| 毛片女人毛片| 国产高潮美女av| 国产精品嫩草影院av在线观看| 男人舔奶头视频| 亚洲精品成人av观看孕妇| 亚洲国产精品sss在线观看| 别揉我奶头 嗯啊视频| 99久久中文字幕三级久久日本| 国产探花在线观看一区二区| 夫妻午夜视频| 天天一区二区日本电影三级| 日韩欧美精品免费久久| 免费观看在线日韩| 午夜爱爱视频在线播放| 免费观看在线日韩| 免费黄色在线免费观看| 国产av在哪里看| 久久精品国产亚洲av天美| 亚洲av福利一区| 精品久久久久久久人妻蜜臀av| 免费看美女性在线毛片视频| 成人亚洲欧美一区二区av| 国产老妇伦熟女老妇高清| 午夜福利在线在线| 99久久人妻综合| 久久久久久久午夜电影| 国产成人一区二区在线| 性插视频无遮挡在线免费观看| 亚洲国产高清在线一区二区三| 五月天丁香电影| 国产精品国产三级专区第一集| 午夜福利成人在线免费观看| 国产 亚洲一区二区三区 | 晚上一个人看的免费电影| 我的女老师完整版在线观看| 午夜免费观看性视频| 欧美激情在线99| 免费观看av网站的网址| 久久久久性生活片| 日韩伦理黄色片| 免费观看精品视频网站| 免费观看性生交大片5| 精品不卡国产一区二区三区| 非洲黑人性xxxx精品又粗又长| 黄片wwwwww| 天堂网av新在线| 一级二级三级毛片免费看| 久久久久久久久久人人人人人人| 午夜爱爱视频在线播放| 深夜a级毛片| 午夜福利成人在线免费观看| 国产成年人精品一区二区| 久久亚洲国产成人精品v| 美女内射精品一级片tv| 纵有疾风起免费观看全集完整版 |