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

    船用高速發(fā)電機(jī)基波磁場(chǎng)下機(jī)座電磁振動(dòng)

    2016-07-26 09:10:48張賢彪顧太平劉德志
    振動(dòng)與沖擊 2016年2期
    關(guān)鍵詞:有限元法

    張賢彪, 王 東, 顧太平, 劉德志

    (海軍工程大學(xué) 艦船綜合電力技術(shù)國防科技重點(diǎn)實(shí)驗(yàn)室,武漢 430033)

    ?

    船用高速發(fā)電機(jī)基波磁場(chǎng)下機(jī)座電磁振動(dòng)

    張賢彪, 王東, 顧太平, 劉德志

    (海軍工程大學(xué) 艦船綜合電力技術(shù)國防科技重點(diǎn)實(shí)驗(yàn)室,武漢430033)

    摘要:以船用高速同步發(fā)電機(jī)為對(duì)象進(jìn)行電磁振動(dòng)特性分析。用有限元法計(jì)算空載下電機(jī)基波氣隙磁場(chǎng)引起的徑向電磁力波在機(jī)座機(jī)腳處垂向振動(dòng)加速度。通過對(duì)不同機(jī)座安裝方式、材料剛度及結(jié)構(gòu)阻尼的振動(dòng)仿真對(duì)比,分析影響振動(dòng)加速度主因及有限元模型簡化處理對(duì)計(jì)算結(jié)果影響。用三維諧波場(chǎng)及瞬態(tài)場(chǎng)有限元法計(jì)算機(jī)腳振動(dòng)加速度,用電動(dòng)機(jī)運(yùn)行法測(cè)試發(fā)電機(jī)空載下機(jī)座機(jī)腳振動(dòng)加速度。結(jié)果表明,仿真與試驗(yàn)結(jié)果基本一致,機(jī)座及機(jī)腳振動(dòng)受機(jī)座安裝方式影響較大,定子振幅受材料彈性模量響應(yīng)較大,在彈性安裝/自由狀態(tài)下機(jī)座形變分布主要受結(jié)構(gòu)阻尼影響。該結(jié)果可為船用高速電機(jī)減振降噪設(shè)計(jì)、研究提供依據(jù)及參考。

    關(guān)鍵詞:電磁振動(dòng);基波磁場(chǎng);船用發(fā)電機(jī);有限元法

    由于對(duì)船舶聲隱身性能要求提高,振動(dòng)、噪聲已成為衡量船用發(fā)電機(jī)的主要技術(shù)指標(biāo)之一。電機(jī)振動(dòng)按產(chǎn)生機(jī)理可分為電磁、機(jī)械及空氣動(dòng)力三類[1]。其中,電磁振動(dòng)主要為氣隙磁場(chǎng)引起電磁力振動(dòng)(暫不考慮磁致伸縮引起的振動(dòng)[2-3]及定子端部磁場(chǎng)引起的繞組振動(dòng)[4]),作為電機(jī)的寄生效應(yīng)[5]廣受關(guān)注。電磁力按氣隙磁場(chǎng)作用方向可分為徑向、切向及軸向力,徑向力中一般考慮基波磁場(chǎng)產(chǎn)生的倍頻力波及由定轉(zhuǎn)子諧波磁場(chǎng)相互作用產(chǎn)生的徑向力波,為電磁振動(dòng)主要激勵(lì)源。各徑向力波作用于定轉(zhuǎn)子鐵心,使定子鐵心、機(jī)座及轉(zhuǎn)子出現(xiàn)隨時(shí)間周期性變化的徑向形變(振幅)。

    目前對(duì)各類電機(jī)電磁振動(dòng)研究多針對(duì)定子鐵心及機(jī)座。對(duì)振動(dòng)激勵(lì)源及傳遞途徑,通過分析徑向電磁力波模型及機(jī)座結(jié)構(gòu)模態(tài),可揭示氣隙長度、槽配合、轉(zhuǎn)子槽形等主要電磁參數(shù)對(duì)電磁力波影響規(guī)律及結(jié)構(gòu)固有振動(dòng)特性。以電磁力波避開結(jié)構(gòu)固有頻率及模態(tài)階次為主要原則,整體優(yōu)化電機(jī)電磁、結(jié)構(gòu)參數(shù)實(shí)現(xiàn)電機(jī)低電磁振動(dòng)[6-12]。對(duì)振動(dòng)響應(yīng),由于氣隙電磁場(chǎng)與定子振動(dòng)為弱耦合關(guān)系[13-14],主要采用電磁場(chǎng)、應(yīng)力場(chǎng)相結(jié)合的有限元法[15];據(jù)電機(jī)實(shí)際安裝方式及分析方便,計(jì)算工況多為電機(jī)空載運(yùn)行,且機(jī)座與基礎(chǔ)之間為剛性約束或自由狀態(tài)[16]。

    由于基波磁場(chǎng)產(chǎn)生的倍頻力波存在階次低、幅值大且難以削弱等特點(diǎn),是電磁振動(dòng)激勵(lì)源的主要分量。實(shí)現(xiàn)船用高速發(fā)電機(jī)低電磁振動(dòng)優(yōu)化設(shè)計(jì),需分析基波磁場(chǎng)產(chǎn)生的倍頻力波振動(dòng)。為此,本文以船用高速同步發(fā)電機(jī)為對(duì)象,用機(jī)腳振動(dòng)加速度參數(shù)作為振動(dòng)響應(yīng)表征,采用磁場(chǎng)、應(yīng)力場(chǎng)相結(jié)合的有限元法計(jì)算空載下電機(jī)基波氣隙磁場(chǎng)產(chǎn)生的倍頻力波振動(dòng)。通過對(duì)不同機(jī)座安裝方式、材料剛度及結(jié)構(gòu)阻尼的振動(dòng)仿真對(duì)比,分析影響振動(dòng)加速度主因及有限元模型簡化對(duì)計(jì)算結(jié)果影響,在此基礎(chǔ)上分別采用三維諧波場(chǎng)及瞬態(tài)場(chǎng)有限元法計(jì)算機(jī)腳振動(dòng)加速度,并實(shí)驗(yàn)驗(yàn)證。

    1電磁振動(dòng)數(shù)學(xué)模型

    1.1電磁振動(dòng)響應(yīng)模型

    多自由度系統(tǒng)由彈性恢復(fù)力、阻尼力、慣性力及電磁激振力形成的平衡系統(tǒng),據(jù)達(dá)朗貝爾原理得運(yùn)動(dòng)方程為

    (1)

    運(yùn)動(dòng)方程在模態(tài)坐標(biāo)系下分解為n個(gè)獨(dú)立方程,令{u}=[s]{η},則有

    (2)

    電機(jī)氣隙基波磁密為B(θ,t)=B0cos(pθ-wet)。其中B0為基波磁密幅值;p為電機(jī)極對(duì)數(shù);θ為定子圓周相位;we為電頻率。作用于定子鐵心內(nèi)表面單位面積的徑向磁拉力正比于磁通密度平方,即:f(θ,t)=B2(θ,t)/(2μ0),其中μ0為氣隙長度。不考慮常數(shù)項(xiàng),得倍頻力波為

    f(θ,t)=f0cos(2wet-2pθ)

    (3)

    據(jù)式(3),設(shè){f}=(f0)cos(2wet),則{h}=[s]T[M]-1(f0)cos(2wet)={H0}cos(2wet),代入式(2)求解得穩(wěn)態(tài)電磁振動(dòng)響應(yīng)為

    (4)

    式中:H0i為模態(tài)坐標(biāo)下第i階廣義電磁力波幅值;φi為模態(tài)坐標(biāo)下第i階節(jié)點(diǎn)位移向量相位。

    已知{u}=[s]{η},代入式(4)得

    (5)

    式中:s(i)為坐標(biāo)系變換矩陣第i階元素。

    由式(5)可知,振動(dòng)響應(yīng)幅度與電磁力波幅值H0i成正比,阻尼比ξi越大振幅ui越小,并與電頻率we及模態(tài)頻率wi密切相關(guān)。以定子電磁振動(dòng)為例,倍頻力波頻率與定子模態(tài)頻率相等且階次相同時(shí)(wi=2we,i=2p),電機(jī)定子振幅最大。

    1.2定子結(jié)構(gòu)模型

    結(jié)構(gòu)固有頻率隨剛度增大而增大,隨質(zhì)量增大而減小。對(duì)形狀及尺寸一定的均勻材料結(jié)構(gòu),增加材料彈性模量或減小材料密度能提高固有頻率,反之則降低。定子繞組由多匝電磁線經(jīng)絕緣包扎、浸漆處理后形成,與實(shí)心導(dǎo)電用銅結(jié)構(gòu)不同,其彈性模量遠(yuǎn)小于實(shí)心結(jié)構(gòu)(E=1.08×1011N/m2),但質(zhì)量與實(shí)心結(jié)構(gòu)相當(dāng)。有限元建模中,繞組多被作為附加質(zhì)量計(jì)入鐵心或鐵心磁極[17]。而對(duì)較大型定子結(jié)構(gòu),有效部分繞組對(duì)固有頻率影響較大,不能簡單以附加質(zhì)量形式計(jì)算[18]。故分析繞組對(duì)振動(dòng)影響時(shí)忽略端部結(jié)構(gòu),有效部分按均勻?qū)嶓w建模,選不同彈性模量進(jìn)行計(jì)算。

    定子鐵心由硅鋼片疊壓而成,與實(shí)心結(jié)構(gòu)不同。在各向異性模型中有一各向同性面(徑向平面xoy),z向?yàn)榀B壓方向,其彈性模量Ez相對(duì)較低。定子鐵心材料參數(shù)模型為

    (6)

    式中:Er,vr,Gxy為徑向彈性模量、泊松比及剪切模量。參數(shù)值由硅鋼材料決定;Ez與疊壓工藝有關(guān),須通過試驗(yàn)測(cè)試。鐵心彈性模量一般取Er=2.0×1011N/m2,Ez=1.5×1011N/m2。

    為分析定子鐵心對(duì)振動(dòng)影響,分別按各向同性、異性對(duì)比計(jì)算。

    1.3機(jī)座結(jié)構(gòu)模型

    船用發(fā)電機(jī)座通過隔振器安裝在基礎(chǔ)上,為彈性安裝。有限元建模中忽略隔振器質(zhì)量只計(jì)其連接剛度,采用彈性模擬。據(jù)隔振器隔振頻率fL及額定載荷ML計(jì)算垂向連接剛度為

    KL=ML(2πfL)2

    (7)

    連接剛度越大機(jī)座固有頻率越高。為分析安裝方式對(duì)振動(dòng)影響,對(duì)比計(jì)算機(jī)座自由狀態(tài)、剛性及彈性安裝三種方式的振動(dòng)響應(yīng)。機(jī)座主要為普通鋼結(jié)構(gòu),阻尼比ξ一般直接按經(jīng)驗(yàn)值0.02計(jì)算,或分解為質(zhì)量阻尼及剛度阻尼系數(shù)后計(jì)算,即

    (8)

    式中:α,β為質(zhì)量、剛度阻尼系數(shù),經(jīng)驗(yàn)值分別取3.0及0.0001。為分析結(jié)構(gòu)阻尼對(duì)振動(dòng)影響,對(duì)比計(jì)算有、無阻尼振動(dòng)響應(yīng)。

    2三維有限元模型及分析

    船用發(fā)電機(jī)系統(tǒng)結(jié)構(gòu)主要包括電機(jī)定子、機(jī)座(隔板、底板、側(cè)板、機(jī)腳)、轉(zhuǎn)子與端蓋式軸承及定子以上部件(整流裝置、勵(lì)磁裝置、風(fēng)機(jī)等)。定子及機(jī)座見圖1,其中1為定子,2為側(cè)板,3為機(jī)腳,4為底板,5為隔板。

    圖1 定子及機(jī)座三維虛擬樣機(jī)圖Fig.1 Structure of stator core and frame

    為便于有限元建模,對(duì)定子及機(jī)座作常規(guī)簡化:① 忽略定子通風(fēng)槽、繞組端部等結(jié)構(gòu)對(duì)計(jì)算影響,將其簡化為圓筒;繞組與鐵心固連,鐵心與隔板固連。② 忽略結(jié)構(gòu)焊接部分,將隔板、側(cè)板、底板各部分視為固連。③ 忽略對(duì)計(jì)算影響不大的孔、槽及倒角等。整機(jī)采用體單元SOLID186,定子模型采用六面體網(wǎng)格劃分。該發(fā)電機(jī)極對(duì)數(shù)為2??蛰d、額定轉(zhuǎn)速下定子徑向電磁力波作用在定子內(nèi)表面,有限元模型見圖2。

    圖2 定子及機(jī)座有限元模型Fig.2 Finite element model of stator core and frame

    分別采用三維諧波場(chǎng)、瞬態(tài)場(chǎng)有限元法計(jì)算機(jī)腳處4個(gè)振動(dòng)測(cè)點(diǎn)的平均加速度。瞬態(tài)場(chǎng)分析采用Full完全法,選加載時(shí)間為5個(gè)周期,每周期分40個(gè)載荷步。為獲得穩(wěn)定振動(dòng),只分析最后一個(gè)周期數(shù)據(jù)結(jié)果。

    3振動(dòng)特性分析

    3.1機(jī)座自由狀態(tài)下振動(dòng)影響因素分析

    3.1.1定子、機(jī)座形變

    機(jī)座為自由狀態(tài)時(shí)機(jī)腳無自由度約束。定子分兩種處理情況:① 低模量:繞組按均勻?qū)嶓w銅建模,不考慮彈性模量,鐵心疊壓結(jié)構(gòu)按各向異性建模,徑向彈性模量取2.0×1011N/m2,疊加方向取1.5×1011N/m2。② 高模量:繞組按均勻?qū)嶓w銅建模,彈性模量取1.08×1011N/m2;鐵心疊壓結(jié)構(gòu)按各向同性建模,彈性模量取2.0×1011N/m2。機(jī)構(gòu)結(jié)構(gòu)阻尼處理為有阻尼時(shí)阻尼比取0.02,無阻尼時(shí)阻尼比取0。定子為低模量下機(jī)座形變見圖3。由圖3可知,有結(jié)構(gòu)阻尼機(jī)座最大形變主要集中在側(cè)板上半部及尾部隔板處;無結(jié)構(gòu)阻尼機(jī)座最大形變量主要集中在側(cè)板、底板頭部。有結(jié)構(gòu)阻尼時(shí)定子形變見圖4。

    圖3 機(jī)座自由狀態(tài)形變?cè)茍DFig.3 Deformation distribution of frame in free installation

    圖4 機(jī)座自由狀態(tài)下的定子形變?cè)茍DFig.4 Deformation distribution of stator core with free installed frame

    機(jī)座自由狀態(tài)下定子、機(jī)座形變見表1。由表1知,有、無結(jié)構(gòu)阻尼下定子最大形變值基本相同。無阻尼下機(jī)座最大形變量較大,主要分布在機(jī)座頭部;相反,有阻尼的機(jī)座最大形變量較小,且主要分布在機(jī)座尾部。機(jī)座及定子最大形變量隨材料彈性模量增大而減小。因此,對(duì)自由狀態(tài),機(jī)座最大形變量及分布主要受阻尼影響,而最大形變量隨阻尼及定子材料的彈性模量增大而減小。定子最大形變量主要受材料彈性模量影響。

    表1 自由狀態(tài)下的定子及機(jī)座形變

    3.1.2機(jī)腳垂向加速度

    自由狀態(tài)下機(jī)腳振動(dòng)加速度見表2。由表2知,較結(jié)構(gòu)阻尼的機(jī)腳垂向振動(dòng),有結(jié)構(gòu)阻尼最大加速度值較小,平均加速度值較大。定子為高模量時(shí)加速度小于低模量。由此得機(jī)座自由狀態(tài)的機(jī)腳垂向平均加速度級(jí)為111~112 dB。

    表2 自由狀態(tài)下機(jī)腳垂向振動(dòng)加速度

    3.2機(jī)座剛性安裝的振動(dòng)影響因素分析

    3.2.1定子、機(jī)座形變

    機(jī)座剛性安裝時(shí),將機(jī)腳與隔振器螺栓連接視為與基礎(chǔ)固連,約束垂向位移為0。定子低模量下機(jī)座形變見圖5。由圖5知,有、無結(jié)構(gòu)阻尼的機(jī)座最大形變量及分布基本一致,形變主要集中在機(jī)座中部及頭部底板。

    圖5 機(jī)座剛性安裝的形變?cè)茍DFig.5 Deformation distribution of frame in fixed installation

    剛性安裝的定子、機(jī)座形變見表3。由表3看出,剛性安裝時(shí)機(jī)座、定子形變量及分布穩(wěn)定,受結(jié)構(gòu)阻尼及材料彈性模量影響較小。

    表3 剛性安裝定子、機(jī)座形變

    3.2.2機(jī)腳垂向振動(dòng)加速度

    剛性安裝的機(jī)腳振動(dòng)加速度見表4。由表4可知,此時(shí)振動(dòng)加速度受結(jié)構(gòu)阻尼影響可忽略;定子為高模量時(shí)機(jī)腳振動(dòng)加速度略小于低模量。

    表4 剛性安裝機(jī)腳垂向振動(dòng)加速度

    3.3機(jī)座彈性安裝的振動(dòng)影響因素分析

    機(jī)座彈性安裝時(shí)隔振器按彈簧單元建模。據(jù)計(jì)算結(jié)果分析知,此時(shí)形變及機(jī)腳振動(dòng)加速度與自由狀態(tài)的值基本一致。定子為低模量的機(jī)座形變見圖6。

    圖6 彈性安裝的機(jī)座形變?cè)茍DFig.6 Deformation distribution of frame in flexible installation

    彈性安裝時(shí)機(jī)腳振動(dòng)加速度見表5。

    表5 彈性安裝機(jī)腳垂向加速度

    3.4振動(dòng)特性分析

    據(jù)機(jī)座不同安裝方式的振動(dòng)影響因素分析知,① 機(jī)座形變主要受安裝方式影響,安裝剛度越大機(jī)座形變量越小。剛性安裝時(shí)機(jī)座形變基本不受結(jié)構(gòu)阻尼及定子材料彈性模量影響;彈性安裝或自由狀態(tài)時(shí)形變分布主要受結(jié)構(gòu)阻尼影響,形變量隨阻尼及彈性模量增大而減小。② 定子形變量主要受彈性模量影響,其次為機(jī)座安裝方式,基本不受結(jié)構(gòu)阻尼影響。彈性模量越高或機(jī)座安裝剛度越大定子最大形變量越小。③ 機(jī)腳垂向加速度受機(jī)座安裝方式影響最大,其次為材料彈性模量,基本不受結(jié)構(gòu)阻尼影響。安裝剛度或彈性模量越大振動(dòng)越小。剛性安裝的彈性模量對(duì)振動(dòng)影響相較彈性安裝小,因安裝方式影響占主導(dǎo)。④ 已知電機(jī)彈性安裝,機(jī)腳加速度分布式測(cè)量,取平均加速度值作為測(cè)試結(jié)果。定子繞組實(shí)際結(jié)構(gòu)中因絕緣、槽楔等其彈性模量較低,繞組按均勻?qū)嶓w銅建??刹挥?jì)彈性模量;鐵心疊壓結(jié)構(gòu)應(yīng)按各向異性建模,因此機(jī)座彈性安裝且定子低模量時(shí)機(jī)腳振動(dòng)平均值加速度更符合電機(jī)實(shí)際安裝及定子鐵心結(jié)構(gòu)。仿真計(jì)算值為111~112 dB。

    4實(shí)驗(yàn)驗(yàn)證及分析

    電磁振動(dòng)實(shí)驗(yàn)中,為避免軸系耦合對(duì)振動(dòng)測(cè)量的影響,將發(fā)電機(jī)作為同步電動(dòng)機(jī)空載運(yùn)行,通過三相變頻器拖動(dòng)至額定轉(zhuǎn)速。機(jī)座通過4個(gè)隔振器與試驗(yàn)底板連接,振動(dòng)測(cè)點(diǎn)布置在機(jī)座與減振器的連接螺栓處,測(cè)量機(jī)腳垂向加速度。電磁振動(dòng)測(cè)試方案見圖7。圖中,1為振動(dòng)測(cè)點(diǎn),2為隔振器,3為試驗(yàn)底板。

    圖7 電磁振動(dòng)測(cè)試方案Fig.7 Test scheme of electromagnetic vibration

    采用同一電機(jī)定子、機(jī)座,分別裝配3種不同結(jié)構(gòu)軸承進(jìn)行振動(dòng)測(cè)試。4個(gè)測(cè)點(diǎn)主要倍頻(對(duì)轉(zhuǎn)頻而言)振動(dòng)分量見表6。據(jù)電機(jī)軸承轉(zhuǎn)速、載荷及穩(wěn)定性要求等選擇可傾瓦軸承,常用結(jié)構(gòu)形式有4片、5片瓦兩種。表6中A為4片瓦軸承;B、C為5片瓦軸承。據(jù)上、下側(cè)瓦片數(shù)不同分別為上側(cè)兩片、下側(cè)三片及上側(cè)三片、下側(cè)兩片?;ù艌?chǎng)產(chǎn)生的倍頻力波為電頻率的2倍、轉(zhuǎn)頻的4倍。由表6知, 4倍轉(zhuǎn)頻平均垂向加速度為0.51~ 0.99 m/s2,轉(zhuǎn)換成分貝為114.2~119.9 dB。仿真結(jié)果為111~112 dB,與實(shí)測(cè)值基本一致,相對(duì)誤差為2%~7%。仿真誤差主要源于有限元求解及模型簡化處理。

    表6 主要倍頻振動(dòng)分量

    為排除有限元求解誤差,對(duì)比三維諧波場(chǎng)及瞬態(tài)場(chǎng)有限元法,所求機(jī)腳垂向加速度相對(duì)誤差見表7。由表7知,用兩種方法計(jì)算的差值較小,相對(duì)誤差在3%以內(nèi),由此知,有限元求解誤差較小。在有限元模型常規(guī)簡化中,由于忽略過盈裝配連接及焊接部分,將機(jī)座、定子鐵心及繞組視為整體,使整個(gè)結(jié)構(gòu)剛性增加,造成計(jì)算結(jié)果偏小。鐵心疊片結(jié)構(gòu)取各向異性,徑向彈性模量按經(jīng)驗(yàn)取值2.0×1011N/m2,該值偏大而使計(jì)算結(jié)果偏小。

    表7 加速度相對(duì)誤差

    5結(jié)論

    通過對(duì)船用高速發(fā)電機(jī)基波磁場(chǎng)下機(jī)座電磁振動(dòng)特性有限元仿真分析,結(jié)論如下:

    (1) 機(jī)腳垂向加速度受機(jī)座安裝方式影響最大,定子材料彈性模量次之,結(jié)構(gòu)阻尼無影響。機(jī)座安裝剛度或定子彈性模量越大振動(dòng)越小。

    (2) 機(jī)座電磁振幅(形變量)及分布主要受機(jī)座安裝方式影響。剛性安裝時(shí)機(jī)座形變基本不受結(jié)構(gòu)阻尼及彈性模量影響;彈性安裝或自由狀態(tài)時(shí)形變分布主要受結(jié)構(gòu)阻尼影響,形變量隨結(jié)構(gòu)阻尼及材料彈性模量增大而減小。

    (3) 定子形變量主要受彈性模量影響,其次為機(jī)座安裝方式,結(jié)構(gòu)阻尼無影響。彈性模量越高或機(jī)座安裝剛度越大定子最大形變量越小。

    (4) 機(jī)腳垂向加速度仿真值與實(shí)測(cè)值基本一致,相對(duì)誤差為2%~7%;有限元求解誤差較小,用兩種有限元求解的機(jī)腳垂向加速度級(jí)相對(duì)誤差在3%以內(nèi);有限元建模中對(duì)各結(jié)構(gòu)簡化處理會(huì)導(dǎo)致仿真誤差。

    參 考 文 獻(xiàn)

    [1] Lisner R P, Timar P L. A new approach to electric motor acoustic noise standards and test Proeedures [J]. IEEE Transactions on Energy Conversion,1999,14(3):692-697.

    [2] 郝清亮,胡義軍,朱少林. 磁致伸縮在電機(jī)電磁振動(dòng)中的貢獻(xiàn)分析[J]. 電機(jī)與控制應(yīng)用,2011, 38(10):31-35.

    HAO Qing-liang, HU Yi-jun,ZHU Shao-lin. Contribution analysis of magnetostriction in motor electromagnetic vibration [J]. Electric Machines & Control Application, 2011,38(10): 31-35.

    [3] Delaere K, Heylen W, Belmans R,et al. Comparison of induction machine stator vibration spectra induced by reluctance forces and magnetostriction[J]. Magnetics, IEEE Transactions on, 2002, 38(2): 969-972.

    [4] 鮑曉華,方勇,程曉巍,等.大型充水式潛水電機(jī)定子端部繞組電磁振動(dòng)研究[J].電工技術(shù)學(xué)報(bào),2013, 28(11):138-146.

    BAO Xiao-hua, FANG Yong, CHENG Xiao-wei, et al. Research on electromagnetic vibration of stator end winding of large water filling submersible motor [J]. Transactions of China Electrotechnical Society, 2013, 28(11): 138-146.

    [5] 楊浩東.永磁同步電機(jī)電磁振動(dòng)分析[D].杭州:浙江大學(xué),2011.

    [6] Sakamoto S,Hirata T. Vibration analysis considering higher harmonics of electromagnetic forces for rotating electricmachines[J]. IEEE Transactions on Magnetics,1999, 35(3): 1662-1665.

    [7] Schlensok C, Schmülling B, van der Giet M, et al. Electromagnetically excited andible noise-evaluation and optimization of electrical machines by numerical simulation[J]. COMPEL: The International Journal for Computation and Mathematics inElectrical and Electronic Engineering, 2007,26(3): 727-742.

    [8] Dai Ying, Cui Shu-mei, Zhang Qian-fan. Analysis on electromagnetic vibration/noise of induction motors for EV drives [J]. Proceedings of the CSEE, 2012, 33(32): 89-97.

    [9] 張磊,溫旭輝.車用永磁同步電機(jī)徑向電磁振動(dòng)特性[J].電機(jī)與控制學(xué)報(bào),2012,16(5):33-39.

    ZHANG Lei, WEN Xu-hui. Radial electromagnetic vibration model characteristics of PMSMs for electric vehicles[J]. Electric Machines and Control, 2012, 16(5): 33-39.

    [10] 鐘焱,丁妍,馮海軍,等.異步電動(dòng)機(jī)電磁振動(dòng)試驗(yàn)研究[J].噪聲與振動(dòng)控制,2011,5:165-167.

    ZHONG Yan, DING Yan, FENG Hai-jun,et al. Experimental study on electromagnetic vibration of asynchronous motor [J].Noise and Vibration Control, 2011, 5:165-167.

    [11] 楊浩東,陳陽生.分?jǐn)?shù)槽永磁同步電機(jī)電磁振動(dòng)的分析與抑制[J].中國電機(jī)工程學(xué)報(bào),2011,24(31):83-89.

    YANG Hao-dong, CHEN Yang-sheng. Electromagnetic vibration analysis and suppression ofpermanent magnet synchronous motor with fractional slot combination[J]. Proceedings of the CSEE, 2011, 24(31): 83-89.

    [12] 左曙光,林福,孫慶,等. 極槽配合和繞組層數(shù)對(duì)永磁同步電機(jī)振動(dòng)的影響分析[J].振動(dòng)與沖擊,2014,33(13):130-134.

    ZUO Shu-guang, LIN Fu, SUN Qing, et al. Vibration of permanent magnet synchronous motor with different slot-pole combinations and winding layers[J].Journal of Vibration and Shock,2014,33(13):130-134.

    [13] Rainer S, Biro O, Weilharter B, et al. Weak coupling between electromagnetic and structural models for electrical machines [J]. Magnetics, IEEE Transactions, 2010,46(8): 2807-2810.

    [14] Lee J H, Lee Y H, Kim D H,et al.Dynamic vibration analysis of switched reluctance motor using magnetic charge force density and mechanical analysis[J]. Applied Superconductivity, IEEE Transactions, 2002, 12(1): 1511-1514.

    [15] 王玎,祝長生,符嘉靖.基于有限元的異步電機(jī)電磁振動(dòng)分析[J].振動(dòng)與沖擊, 2012, 31(2):140-144.

    WANG Ding, ZHU Chang-sheng, FU Jia-jing. Electromagnetically excited vibration analysis for an asynchronous electrical machine with finite element method[J]. Journal of Vibration and Shock, 2012, 31(2): 140-144.

    [16] 韋春健,劉德志,王東.空載情況下三相異步電動(dòng)機(jī)徑向電磁振動(dòng)特性[J].電機(jī)與控制應(yīng)用,2010, 37(12):13-17.

    WEI Chun-jian, LIU De-zhi, WANG Dong. Three-phase asynchronous motor radial electromagnetic vibration analysis without loads[J]. Electric Machines & Control Application, 2010, 37(12): 13-18.

    [17] 吳建華. 基于物理模型開關(guān)磁阻電機(jī)定子模態(tài)和固有頻率的研究[J]. 中國電機(jī)工程學(xué)報(bào), 2004, 24(8): 110-114.

    WU Jian-hua. Study on the stator mode shapes and resonant frequencies of switched reluctance motor based on real structural model[J]. Proceedings of the CSEE, 2004, 24(8): 110-114.

    [18] 王天煜, 王鳳翔. 大型異步電動(dòng)機(jī)定子振動(dòng)與模態(tài)分析[J]. 中國電機(jī)工程學(xué)報(bào), 2007, 27(12): 41-45.

    WANG Tian-yu, WANG Feng-xiang. Vibration and modal analysis of stator of large induction motors[J]. Proceedings of the CSEE, 2007, 27(12):41-45.

    基金項(xiàng)目:國家自然科學(xué)重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃973項(xiàng)目(2013CB035601);國家自然科學(xué)基金資助項(xiàng)目(51137005,51222705)

    收稿日期:2014-09-18修改稿收到日期:2014-11-06

    通信作者王東 男,博士,教授,博士生導(dǎo)師,1978年生

    中圖分類號(hào):TH113

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

    DOI:10.13465/j.cnki.jvs.2016.02.033

    Electromagnetic vibration characteristics of marine high-speed generator framederived from fundamental magnetic field

    ZHANG Xian-biao, WANG Dong, GU Tai-ping, LIU De-zhi

    (National Key Laboratory for Vessel Integrated Power System Technology,Naval University of Engineering, Wuhan 430033, China)

    Abstract:With the increase of speed and power-density, the vibration problem of marine high-speed generator has become much more severe. Aiming at this problem, the vertical vibration acceleration at the machine feet of a marine high-speed synchronous generator was studied by finite element method (FEM), which is derived from the radial electromagnetic force due to the fundamental air-gap magnetic field in no-load condition. The vibration characteristicsis were analysed under different type of frame installations, different material stiffness and different structural damping. The main influential factors were classified clearly by comparative analysis. Moreover, the simplification of finite element model was introduced. The vibration acceleration at the machine feet was calculated by using the three dimensional harmonic field FEM and the transient field FEM respectively. The tested vertical vibration acceleration at machine feet of marine high-speed synchronous generator when running as a no-load motor was presented. The simulation results were compared with the test results, where a substantial agreement was achieved. The results show that frame and feet’s vibrations are impacted by installations greatly and the material stiffness is the main influential factor of stator’s vibration amplitude. Moreover, the frame’s damping has a significant effect on the vibration amplitude distribution under free or flexible installation. The results provide a theoretical reference to the analysis of marine high-speed generator vibration and noise reduction.

    Key words:electromagnetic vibration; fundamental magnetic field; marine generator; finite element method

    第一作者 張賢彪 男,博士生,1986年10月生

    猜你喜歡
    有限元法
    正交各向異性材料裂紋疲勞擴(kuò)展的擴(kuò)展有限元法研究
    基于有限元法的高頻變壓器繞組損耗研究
    基于有限元法副發(fā)動(dòng)機(jī)托架輕量化設(shè)計(jì)
    專用汽車(2016年8期)2016-03-01 04:16:43
    傳遞矩陣法與有限元法計(jì)算電機(jī)轉(zhuǎn)子臨界轉(zhuǎn)速的對(duì)比分析
    Sine-Gordon方程H1-Galerkin非協(xié)調(diào)混合有限元法的誤差分析
    三維有限元法在口腔正畸生物力學(xué)研究中發(fā)揮的作用
    RKDG有限元法求解一維拉格朗日形式的Euler方程
    集成對(duì)稱模糊數(shù)及有限元法的切削力預(yù)測(cè)
    有限元法在機(jī)械設(shè)計(jì)方向中的教學(xué)實(shí)踐
    基于HCSR和CSR-OT的油船疲勞有限元法對(duì)比分析
    船海工程(2013年6期)2013-03-11 18:57:25
    午夜久久久久精精品| 成年版毛片免费区| 国产乱人视频| 色吧在线观看| 能在线免费看毛片的网站| 精品无人区乱码1区二区| 国产又黄又爽又无遮挡在线| 欧美日韩国产亚洲二区| 91午夜精品亚洲一区二区三区| 少妇丰满av| 色哟哟·www| 欧美成人一区二区免费高清观看| 国产v大片淫在线免费观看| 亚洲丝袜综合中文字幕| 亚洲国产高清在线一区二区三| av在线观看视频网站免费| 国产av不卡久久| 观看美女的网站| 成人av在线播放网站| 精品久久久久久久久久久久久| 久久久久久久久久成人| 亚洲欧美精品专区久久| 一级毛片久久久久久久久女| 成人鲁丝片一二三区免费| 波多野结衣巨乳人妻| 亚洲电影在线观看av| 午夜a级毛片| 搡老妇女老女人老熟妇| 中文字幕制服av| 久久亚洲国产成人精品v| ponron亚洲| 日韩三级伦理在线观看| 一本久久精品| 国产在视频线在精品| 两性午夜刺激爽爽歪歪视频在线观看| 日产精品乱码卡一卡2卡三| 九色成人免费人妻av| 久久久久九九精品影院| 亚洲在线观看片| 国产高清不卡午夜福利| 我的女老师完整版在线观看| 日本一二三区视频观看| 高清在线视频一区二区三区 | 久久精品人妻少妇| 男女边吃奶边做爰视频| 久久久色成人| 天堂网av新在线| 老司机影院毛片| 91av网一区二区| 亚洲精品国产成人久久av| 日韩亚洲欧美综合| 中文字幕人妻熟人妻熟丝袜美| 日韩,欧美,国产一区二区三区 | 久热久热在线精品观看| www.av在线官网国产| 波多野结衣高清无吗| 国产一区二区在线av高清观看| 国产在视频线精品| 中文字幕久久专区| 22中文网久久字幕| 我要看日韩黄色一级片| 麻豆av噜噜一区二区三区| 免费av观看视频| 一区二区三区高清视频在线| 在线免费观看的www视频| 久久午夜福利片| 亚洲久久久久久中文字幕| 免费看光身美女| 日本五十路高清| 久久久久久久午夜电影| 日韩一区二区三区影片| 欧美日本视频| 深夜a级毛片| 中文字幕熟女人妻在线| 成人鲁丝片一二三区免费| 尾随美女入室| 亚洲精品久久久久久婷婷小说 | 高清毛片免费看| 少妇熟女aⅴ在线视频| 欧美日韩在线观看h| 成人漫画全彩无遮挡| 韩国av在线不卡| 精品国内亚洲2022精品成人| 国产乱人视频| 久久久久久久午夜电影| 女人被狂操c到高潮| 欧美潮喷喷水| 神马国产精品三级电影在线观看| or卡值多少钱| 如何舔出高潮| 欧美最新免费一区二区三区| 亚洲一级一片aⅴ在线观看| 能在线免费看毛片的网站| 亚洲国产精品国产精品| 99久久九九国产精品国产免费| 在线观看一区二区三区| 性色avwww在线观看| 99久久人妻综合| 国产精品麻豆人妻色哟哟久久 | 丰满人妻一区二区三区视频av| 国语对白做爰xxxⅹ性视频网站| 国产在线一区二区三区精 | 一级毛片久久久久久久久女| 内地一区二区视频在线| 国产精品一区二区性色av| 亚洲av福利一区| 99久久精品国产国产毛片| 黄色欧美视频在线观看| 国产精品国产三级专区第一集| 久久久久久大精品| 1024手机看黄色片| 午夜免费男女啪啪视频观看| 亚洲精品亚洲一区二区| 69人妻影院| 99久国产av精品| 国产精品乱码一区二三区的特点| 美女国产视频在线观看| 久久久国产成人免费| 国产淫片久久久久久久久| 十八禁国产超污无遮挡网站| 纵有疾风起免费观看全集完整版 | 中国国产av一级| 一级黄片播放器| 国产午夜精品一二区理论片| 国产 一区精品| 日本av手机在线免费观看| 嘟嘟电影网在线观看| 伊人久久精品亚洲午夜| 色综合亚洲欧美另类图片| 一区二区三区免费毛片| 夜夜爽夜夜爽视频| av在线观看视频网站免费| 午夜亚洲福利在线播放| 亚洲图色成人| 成人毛片60女人毛片免费| 久久精品影院6| 亚洲国产色片| 少妇熟女aⅴ在线视频| 国内揄拍国产精品人妻在线| 免费观看的影片在线观看| 久久久久久久亚洲中文字幕| 日产精品乱码卡一卡2卡三| 青青草视频在线视频观看| 国产精品蜜桃在线观看| 黄色一级大片看看| 少妇被粗大猛烈的视频| 亚洲四区av| 精品国内亚洲2022精品成人| 国产真实伦视频高清在线观看| 久久99热这里只频精品6学生 | 欧美日本亚洲视频在线播放| 嫩草影院入口| 成人综合一区亚洲| 美女大奶头视频| 一夜夜www| 99久久成人亚洲精品观看| 国产免费福利视频在线观看| av在线天堂中文字幕| 蜜桃亚洲精品一区二区三区| 一个人看的www免费观看视频| 亚洲精品色激情综合| 成人av在线播放网站| 亚洲成av人片在线播放无| 欧美又色又爽又黄视频| 免费搜索国产男女视频| 亚洲av成人精品一区久久| 国产高潮美女av| 国产精品三级大全| 美女高潮的动态| 欧美不卡视频在线免费观看| 国产亚洲91精品色在线| 国产单亲对白刺激| 在线观看66精品国产| 高清日韩中文字幕在线| 毛片女人毛片| 亚洲欧美成人精品一区二区| 精品久久久久久电影网 | 日本五十路高清| 久久精品久久精品一区二区三区| 乱码一卡2卡4卡精品| 我的老师免费观看完整版| 综合色av麻豆| 国产在线一区二区三区精 | 99国产精品一区二区蜜桃av| 国产免费又黄又爽又色| 精品人妻视频免费看| 免费搜索国产男女视频| 国产久久久一区二区三区| 插逼视频在线观看| 欧美性猛交黑人性爽| 亚洲一区高清亚洲精品| 大香蕉久久网| 色综合站精品国产| av线在线观看网站| 亚洲精品久久久久久婷婷小说 | 欧美日韩国产亚洲二区| 免费电影在线观看免费观看| 成年免费大片在线观看| 天堂中文最新版在线下载 | 中文字幕av成人在线电影| 日韩成人av中文字幕在线观看| 级片在线观看| 男女啪啪激烈高潮av片| 午夜免费激情av| 免费av观看视频| 亚洲综合色惰| 国产乱来视频区| 九草在线视频观看| 国产国拍精品亚洲av在线观看| 久久久国产成人精品二区| 插逼视频在线观看| 亚洲国产精品成人综合色| 免费观看精品视频网站| 美女国产视频在线观看| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 国产黄色小视频在线观看| 亚洲在线观看片| 婷婷色综合大香蕉| 国产精品久久久久久精品电影小说 | 久久精品人妻少妇| 国产精华一区二区三区| 亚洲经典国产精华液单| 最近中文字幕高清免费大全6| 一级爰片在线观看| 国产伦精品一区二区三区视频9| 成人性生交大片免费视频hd| 日本免费在线观看一区| 69av精品久久久久久| 国产精品电影一区二区三区| 国产成人免费观看mmmm| 老司机福利观看| 午夜福利在线观看免费完整高清在| 搡老妇女老女人老熟妇| 日韩在线高清观看一区二区三区| 国产高清国产精品国产三级 | 日韩强制内射视频| 免费观看精品视频网站| 高清日韩中文字幕在线| 欧美高清性xxxxhd video| 亚洲激情五月婷婷啪啪| 男插女下体视频免费在线播放| av在线亚洲专区| 九九久久精品国产亚洲av麻豆| 久久精品国产亚洲网站| 午夜福利网站1000一区二区三区| 亚洲精品影视一区二区三区av| 国产午夜精品一二区理论片| 国产成人精品婷婷| 亚洲国产成人一精品久久久| 日产精品乱码卡一卡2卡三| 久久久久九九精品影院| 女人久久www免费人成看片 | 国产毛片a区久久久久| 99久久中文字幕三级久久日本| 免费看av在线观看网站| 亚洲中文字幕日韩| 久久久精品欧美日韩精品| 久久99热6这里只有精品| 日韩一区二区视频免费看| 国产不卡一卡二| av免费观看日本| 美女内射精品一级片tv| 日韩欧美国产在线观看| 老司机影院毛片| 超碰av人人做人人爽久久| 啦啦啦啦在线视频资源| 中文亚洲av片在线观看爽| 中文欧美无线码| 搡老妇女老女人老熟妇| 国产成人精品一,二区| av国产免费在线观看| 国产精品久久久久久av不卡| 久久久久久久亚洲中文字幕| 人妻制服诱惑在线中文字幕| 激情 狠狠 欧美| 成人鲁丝片一二三区免费| 不卡视频在线观看欧美| av免费在线看不卡| 晚上一个人看的免费电影| 伦理电影大哥的女人| 国内揄拍国产精品人妻在线| 能在线免费观看的黄片| 国产精品乱码一区二三区的特点| 菩萨蛮人人尽说江南好唐韦庄 | 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 日本猛色少妇xxxxx猛交久久| 国产久久久一区二区三区| 全区人妻精品视频| 国产一区二区亚洲精品在线观看| 99久国产av精品| 汤姆久久久久久久影院中文字幕 | 久久久精品欧美日韩精品| 亚洲一区高清亚洲精品| 欧美区成人在线视频| 一夜夜www| 99久久人妻综合| 成年女人看的毛片在线观看| 日本一二三区视频观看| 久久久久久久久中文| 国产三级在线视频| 精品无人区乱码1区二区| 国产探花极品一区二区| 日韩成人av中文字幕在线观看| 别揉我奶头 嗯啊视频| 国产极品天堂在线| 18禁在线无遮挡免费观看视频| 久久久久久久久久成人| 亚洲av福利一区| 级片在线观看| 亚洲国产精品sss在线观看| 亚洲高清免费不卡视频| 免费人成在线观看视频色| 自拍偷自拍亚洲精品老妇| 国产美女午夜福利| 国产精品.久久久| 国产成年人精品一区二区| 精品欧美国产一区二区三| .国产精品久久| 国产一区亚洲一区在线观看| 日韩一区二区三区影片| 国产亚洲5aaaaa淫片| 久久午夜福利片| 联通29元200g的流量卡| 国产av在哪里看| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 国产av不卡久久| 国产 一区 欧美 日韩| 日韩欧美 国产精品| 欧美性猛交╳xxx乱大交人| 欧美高清性xxxxhd video| 亚洲国产精品sss在线观看| 能在线免费看毛片的网站| 大又大粗又爽又黄少妇毛片口| 日本与韩国留学比较| 国产人妻一区二区三区在| 亚洲在久久综合| 免费电影在线观看免费观看| av黄色大香蕉| 汤姆久久久久久久影院中文字幕 | 精品久久久久久久久av| 亚洲欧美精品综合久久99| 黄色日韩在线| 卡戴珊不雅视频在线播放| 精品久久久噜噜| 久久精品国产自在天天线| 欧美成人精品欧美一级黄| 精品国产露脸久久av麻豆 | 国产精品av视频在线免费观看| 最后的刺客免费高清国语| 老司机影院毛片| 精品午夜福利在线看| 国产精品一二三区在线看| 国产精品乱码一区二三区的特点| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| 日本午夜av视频| 欧美3d第一页| av在线蜜桃| 18禁在线播放成人免费| 久久亚洲国产成人精品v| 亚洲综合色惰| 如何舔出高潮| www日本黄色视频网| 久久久久网色| 亚洲国产精品合色在线| 91在线精品国自产拍蜜月| 观看免费一级毛片| av播播在线观看一区| 久久久久性生活片| 亚洲第一区二区三区不卡| 99九九线精品视频在线观看视频| 色5月婷婷丁香| 99久久中文字幕三级久久日本| 国产黄色视频一区二区在线观看 | 久久久久久伊人网av| 亚洲av免费高清在线观看| 寂寞人妻少妇视频99o| 日本免费a在线| 十八禁国产超污无遮挡网站| 97人妻精品一区二区三区麻豆| 国产69精品久久久久777片| 亚洲国产高清在线一区二区三| 成人欧美大片| 成年女人永久免费观看视频| 亚洲三级黄色毛片| 色吧在线观看| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 最近的中文字幕免费完整| 久久这里有精品视频免费| 看非洲黑人一级黄片| 老司机影院成人| 欧美不卡视频在线免费观看| 亚洲婷婷狠狠爱综合网| 97超碰精品成人国产| 99久久中文字幕三级久久日本| 久久久久久久亚洲中文字幕| 最新中文字幕久久久久| 国产老妇伦熟女老妇高清| 天堂中文最新版在线下载 | 欧美丝袜亚洲另类| 色网站视频免费| АⅤ资源中文在线天堂| 国产精品电影一区二区三区| 亚洲国产精品专区欧美| 国产精品伦人一区二区| 精品久久久久久久久亚洲| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99 | 国产成年人精品一区二区| 亚洲在线观看片| 亚洲国产精品sss在线观看| 亚洲av免费高清在线观看| 国产精品乱码一区二三区的特点| 欧美xxxx黑人xx丫x性爽| 亚洲成av人片在线播放无| www.色视频.com| 国产高清三级在线| 久久这里只有精品中国| 观看免费一级毛片| 国产不卡一卡二| 乱系列少妇在线播放| 日韩av在线免费看完整版不卡| 亚洲婷婷狠狠爱综合网| 日韩国内少妇激情av| 性插视频无遮挡在线免费观看| 国产又色又爽无遮挡免| 久久久久久久久大av| 又黄又爽又刺激的免费视频.| 五月伊人婷婷丁香| 狂野欧美白嫩少妇大欣赏| 中文欧美无线码| 久久久色成人| 一级黄色大片毛片| 精品欧美国产一区二区三| 少妇熟女aⅴ在线视频| 观看美女的网站| 国产激情偷乱视频一区二区| 最近最新中文字幕免费大全7| 韩国av在线不卡| 日韩欧美精品v在线| 能在线免费观看的黄片| 欧美三级亚洲精品| 日韩亚洲欧美综合| 精品熟女少妇av免费看| 国产美女午夜福利| 日韩欧美 国产精品| 亚洲国产最新在线播放| 纵有疾风起免费观看全集完整版 | 成人漫画全彩无遮挡| 国产 一区 欧美 日韩| 少妇被粗大猛烈的视频| 国产 一区 欧美 日韩| 久久精品久久久久久噜噜老黄 | 亚洲真实伦在线观看| 国内精品宾馆在线| 亚洲av中文av极速乱| 国产精品熟女久久久久浪| 黄色日韩在线| 中文乱码字字幕精品一区二区三区 | 少妇熟女欧美另类| 网址你懂的国产日韩在线| 久久久精品欧美日韩精品| av专区在线播放| av在线天堂中文字幕| 秋霞伦理黄片| or卡值多少钱| 看黄色毛片网站| 国产精品一区二区三区四区久久| 97超碰精品成人国产| 久久精品影院6| 精品人妻一区二区三区麻豆| 天堂√8在线中文| 国产综合懂色| 免费播放大片免费观看视频在线观看 | 亚洲欧美日韩高清专用| av在线观看视频网站免费| 好男人在线观看高清免费视频| 只有这里有精品99| 国产视频内射| 男的添女的下面高潮视频| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 99久国产av精品国产电影| 日韩一本色道免费dvd| 美女内射精品一级片tv| 97超碰精品成人国产| 日本猛色少妇xxxxx猛交久久| 高清在线视频一区二区三区 | 中文字幕熟女人妻在线| 99久久精品一区二区三区| 免费电影在线观看免费观看| 日韩欧美在线乱码| 国产国拍精品亚洲av在线观看| 国产精品久久电影中文字幕| 成人国产麻豆网| 99久国产av精品| 爱豆传媒免费全集在线观看| 亚洲天堂国产精品一区在线| 国内揄拍国产精品人妻在线| 秋霞伦理黄片| 欧美一区二区国产精品久久精品| 日韩国内少妇激情av| 婷婷六月久久综合丁香| 成人av在线播放网站| 永久免费av网站大全| 亚洲一区高清亚洲精品| 少妇人妻一区二区三区视频| 少妇熟女欧美另类| 免费av毛片视频| 日本三级黄在线观看| 99在线人妻在线中文字幕| 精品久久久噜噜| 日韩制服骚丝袜av| 国产成人午夜福利电影在线观看| 18禁在线无遮挡免费观看视频| 中文字幕av成人在线电影| 久久久国产成人免费| 日本午夜av视频| 啦啦啦啦在线视频资源| 亚洲一级一片aⅴ在线观看| 亚洲天堂国产精品一区在线| 2022亚洲国产成人精品| 色尼玛亚洲综合影院| 色综合亚洲欧美另类图片| 人人妻人人看人人澡| 欧美性猛交黑人性爽| 国产精品野战在线观看| 欧美成人a在线观看| 伦理电影大哥的女人| 久99久视频精品免费| 国产成人精品久久久久久| 免费看日本二区| 日韩欧美精品v在线| 国产精品国产三级国产专区5o | 久久久久九九精品影院| 日日啪夜夜撸| 青春草国产在线视频| 午夜福利在线在线| 老女人水多毛片| av免费在线看不卡| 伦精品一区二区三区| 免费观看性生交大片5| 国产又色又爽无遮挡免| 国产精品爽爽va在线观看网站| 欧美变态另类bdsm刘玥| 午夜福利在线观看吧| 一级黄色大片毛片| 能在线免费看毛片的网站| 亚洲真实伦在线观看| 婷婷色av中文字幕| 长腿黑丝高跟| 成人欧美大片| 好男人视频免费观看在线| 一个人看的www免费观看视频| 在线免费观看不下载黄p国产| 九九在线视频观看精品| 禁无遮挡网站| 国产真实伦视频高清在线观看| 成人二区视频| 高清av免费在线| 国产淫语在线视频| 亚洲av福利一区| a级一级毛片免费在线观看| 少妇的逼好多水| 久久久欧美国产精品| a级毛片免费高清观看在线播放| 小说图片视频综合网站| 午夜福利视频1000在线观看| 欧美日韩在线观看h| 91精品国产九色| 亚洲乱码一区二区免费版| 中文字幕人妻熟人妻熟丝袜美| 精品午夜福利在线看| 99久久九九国产精品国产免费| 成人无遮挡网站| 欧美日韩在线观看h| 亚洲人成网站在线播| 国产精品人妻久久久影院| 大又大粗又爽又黄少妇毛片口| 色综合站精品国产| 中文亚洲av片在线观看爽| 亚洲国产色片| 菩萨蛮人人尽说江南好唐韦庄 | 日本熟妇午夜| 在线免费十八禁| 欧美精品一区二区大全| 男人狂女人下面高潮的视频| 精品一区二区免费观看| 99久久成人亚洲精品观看| 在线观看美女被高潮喷水网站| 亚洲,欧美,日韩| 日韩亚洲欧美综合| 婷婷色麻豆天堂久久 | 免费无遮挡裸体视频| 又粗又硬又长又爽又黄的视频| 99久久九九国产精品国产免费| 精品久久久久久久末码| 欧美日韩在线观看h| 成人鲁丝片一二三区免费| 中文欧美无线码| 久久这里有精品视频免费| 中文字幕熟女人妻在线| 国内精品美女久久久久久| 国产美女午夜福利| 亚洲国产精品专区欧美| 亚洲三级黄色毛片| 午夜福利在线观看吧| 永久免费av网站大全| 精品人妻一区二区三区麻豆| 亚洲国产欧洲综合997久久,| 国产精品.久久久|