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

    冰載荷電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)研究

    2018-06-05 07:06:03肖能齊周瑞平侯之超
    船舶力學(xué) 2018年5期
    關(guān)鍵詞:冰區(qū)階次軸系

    肖能齊 , 周瑞平 , 徐 翔 , 侯之超

    (1.水電機(jī)械設(shè)備設(shè)計(jì)與維護(hù)湖北省重點(diǎn)實(shí)驗(yàn)室(三峽大學(xué)),湖北 宜昌 443002;2.汽車安全與能國(guó)家重點(diǎn)實(shí)驗(yàn)室(清華大學(xué)),北京 100084;3.武漢理工大學(xué) 能源與動(dòng)力工程學(xué)院,武漢 430063)

    0 引 言

    為爭(zhēng)奪北極地區(qū)漁業(yè)資源和礦產(chǎn)資源,現(xiàn)今各國(guó)已經(jīng)開始設(shè)計(jì)和制造冰區(qū)航行的科學(xué)考察船和鉆井船等[1-2];同時(shí),為縮短貨物的航程、降低航運(yùn)成本從而提高競(jìng)爭(zhēng)力,航運(yùn)界開始運(yùn)用冰區(qū)貨船、冰區(qū)油船和破冰LNG船舶[3]等航行于北極航道中,目前柴電混合動(dòng)力推進(jìn)系統(tǒng)和純電力推進(jìn)系統(tǒng)被廣泛用于冰區(qū)航行船舶中。為了冰區(qū)航行船舶電力推進(jìn)軸系的安全性,各國(guó)船級(jí)社對(duì)冰區(qū)航行船舶軸系扭轉(zhuǎn)振動(dòng)提出了相應(yīng)的要求。

    目前國(guó)內(nèi)外學(xué)者從冰載荷下的螺旋槳激勵(lì)力[4]、電機(jī)振動(dòng)激勵(lì)力和冰載荷下的軸系扭轉(zhuǎn)振動(dòng)算法等三個(gè)方面進(jìn)行了一定的研究。例如,Dan和Bose等[5]以加拿大海岸警衛(wèi)隊(duì)破冰船的螺旋槳為對(duì)象建立螺旋槳模型,借助理論計(jì)算與實(shí)驗(yàn)測(cè)試手段研究了螺旋槳與冰塊相互作用過程中螺旋槳的水動(dòng)力特性以及冰塊對(duì)螺旋槳推力和轉(zhuǎn)矩的變化。胡志寬[6]提出利用ANSYS軟件對(duì)所建立的螺旋槳進(jìn)行了冰載荷作用下的靜力分析;最后,借助冰塊與槳葉的碰撞模型,分別進(jìn)行了單個(gè)冰塊、兩個(gè)冰塊和三個(gè)冰塊與螺旋槳槳葉的碰撞分析。Barro[7]提出冰塊作用在螺旋槳上負(fù)載與冰的強(qiáng)度指數(shù)、原動(dòng)機(jī)的類型以及推進(jìn)系統(tǒng)剛?cè)崽匦杂嘘P(guān);同時(shí)考慮冰載荷下的螺旋槳激勵(lì)力,對(duì)冰區(qū)柴油機(jī)推進(jìn)軸系扭轉(zhuǎn)振動(dòng)為例進(jìn)行了理論計(jì)算與測(cè)試分析。江攀[8]對(duì)電機(jī)電磁振動(dòng)和螺旋槳的激勵(lì)力進(jìn)行了研究,利用ANSYS軟件在考慮激勵(lì)力作用下對(duì)電力推進(jìn)軸系進(jìn)行了縱向振動(dòng)和橫向振動(dòng)分析。肖能齊等[9]通過建立電機(jī)的機(jī)電耦合振動(dòng)模型,分析了電機(jī)在起動(dòng)過程中的扭轉(zhuǎn)振動(dòng)特性;同時(shí)以1000噸潤(rùn)滑油船舶電力推進(jìn)軸系為研究對(duì)象,在不考慮電磁勵(lì)磁轉(zhuǎn)矩和考慮電磁激振力矩2種情況下,進(jìn)行了強(qiáng)迫振動(dòng)特性對(duì)比分析分析。楊紅軍等[10]闡述了冰區(qū)航行的柴油機(jī)推進(jìn)軸系螺旋槳與冰塊相互作用的螺旋槳激勵(lì)力,使用Newmark法對(duì)軸系進(jìn)行扭轉(zhuǎn)振動(dòng)計(jì)算,其瞬態(tài)計(jì)算結(jié)果與傳遞矩陣法計(jì)算結(jié)果和扭轉(zhuǎn)振動(dòng)測(cè)試結(jié)果一致。

    本文以冰載荷混合動(dòng)力推進(jìn)系統(tǒng)的PTI推進(jìn)模式為研究對(duì)象,首先對(duì)冰載荷與螺旋槳相互作用過程進(jìn)行分析,從而對(duì)冰載荷激勵(lì)力矩的時(shí)域和頻域特性進(jìn)行研究;其次研究推進(jìn)電機(jī)氣隙分布不均勻?qū)е碌碾姶呸D(zhuǎn)矩波動(dòng)而產(chǎn)生的階次振動(dòng)激勵(lì)力;最后運(yùn)用系統(tǒng)矩陣法對(duì)僅考慮電機(jī)的階次激勵(lì)力和冰載荷激勵(lì)力與電機(jī)階次激勵(lì)力共同作用二種情況對(duì)軸系的扭轉(zhuǎn)振動(dòng)特性展開研究,并采用應(yīng)變法進(jìn)行實(shí)船測(cè)試與分析,驗(yàn)證相關(guān)理論和算法的正確性。

    1 冰載荷作用下的螺旋槳激勵(lì)力

    1.1 螺旋槳與冰塊相互作用過程

    船舶推進(jìn)軸系在冰區(qū)航行時(shí),槳葉與冰塊間將產(chǎn)生接觸負(fù)荷,稱為螺旋槳與冰塊間的相互作用,圖1所示為冰區(qū)航行中的電力推進(jìn)軸系示意圖。螺旋槳與冰塊間的相互作用過程可以分為三個(gè)階段:(1)冰塊與槳葉接觸-加載:該階段葉片與冰塊接觸,水與冰塊之間的相對(duì)運(yùn)動(dòng)將增加,同時(shí)載荷也增加;(2)槳葉穿透:隨著載荷的增加達(dá)到一定值時(shí),冰塊破碎后葉片將穿透冰塊破碎區(qū);(3)冰塊與槳葉分離:葉片與冰塊逐漸分離,載荷逐漸減小直至槳葉與冰塊之間分離,此時(shí)螺旋槳與冰塊之間的相互作用負(fù)荷為零。螺旋槳與冰塊間的相互作用如圖2所示。

    根據(jù)如圖2所示冰區(qū)航行船舶的螺旋槳槳葉與冰塊的相互作用過程,將螺旋槳單個(gè)葉片對(duì)冰塊產(chǎn)生的扭矩激勵(lì)分為3個(gè)工況,建立單個(gè)葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力函數(shù):

    圖1 冰區(qū)航行船舶電力推進(jìn)軸系示意圖Fig.1 Ice ships electric propulsion shafting diagram

    圖2 螺旋槳槳葉與冰塊的相互作用過程Fig.2 Interaction of the propeller blades and ice

    其中:Qmax為螺旋槳與冰相互作用在螺旋槳上產(chǎn)生的最大冰塊扭矩;θ為螺旋槳旋轉(zhuǎn)角度;Cq和am與螺旋槳和冰塊相互作用的扭矩激勵(lì)工況有關(guān),如表1所示。表中冰區(qū)加強(qiáng)船舶Cq和am的取值中,Z為螺旋槳的葉片數(shù)。

    表1 冰區(qū)加強(qiáng)船舶Cq和am的取值Tab.1 Ice strengthened vessel Cq and am

    1.2 冰載荷激勵(lì)力

    以5葉螺旋槳為例,在葉片1至葉片5將依次銑削冰塊過程中,根據(jù)葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力函數(shù)(1)式和(2)式,可以得到如圖3所示的工況1-單葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力曲線。

    根據(jù)圖3所示的扭矩激勵(lì)力變化曲線可知:由于5葉螺旋槳相鄰槳葉夾角為72°,葉片與冰塊從進(jìn)入冰塊與槳葉接觸-加載階段到葉片與冰塊逐漸分離的角度是72°,因此葉片1-5依次銑削冰塊,曲線不存在重疊區(qū)域。

    同理可以得到如圖4所示的工況2-單葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力的變化曲線。

    圖3 工況1-單葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力變化Fig.3 Case 1-change of torque excitation force produced by single blade and ice cube

    圖4 工況2-單葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力變化Fig.4 Case 2-change of torque excitation force produced by single blade and ice cube

    根據(jù)圖4所示的曲線圖可知:螺旋槳相鄰槳葉夾角為72°,而葉片與冰塊從進(jìn)入冰塊與槳葉接觸-加載階段到葉片與冰塊分離的角度是135°,曲線重疊區(qū)域是由于葉片與冰塊之間未分離時(shí),而相鄰葉片已經(jīng)開始切削冰塊。

    同理可以得到如圖5所示的工況3-單葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力的變化曲線。

    根據(jù)圖5所示曲線圖可知:由于兩個(gè)冰塊相位差36°,而螺旋槳相鄰槳葉夾角為72°,因此在0°到36°和在36°到72°范圍內(nèi)葉片1分別銑削第1個(gè)冰塊和第2個(gè)冰塊,其扭矩激勵(lì)力曲線分別為圖中的藍(lán)色實(shí)線和紅色虛線。

    圖5 工況3-葉片與冰塊沖擊產(chǎn)生的扭矩激勵(lì)力變化Fig.5 Case 3-change of torque excitation force produced by single blade and ice cube

    圖6 工況1-工況3螺旋槳與冰塊沖擊產(chǎn)生的總扭矩激勵(lì)力與旋轉(zhuǎn)角的關(guān)系Fig.6 The relationship between the total torque excitation force and the rotation angle of the propeller and the ice cube in case1 to 3

    為了得到槳葉與冰塊相互作用的總扭矩激勵(lì)力,將圖3-5所示的工況1-工況3中的螺旋槳各槳葉的扭矩激勵(lì)力進(jìn)行疊加后可以得到如圖9所示的在工況1-工況3下螺旋槳與冰塊沖擊產(chǎn)生的總扭矩激勵(lì)力與旋轉(zhuǎn)角的關(guān)系。

    將其進(jìn)行FFT變換,對(duì)螺旋槳與冰塊相互作用產(chǎn)生的總扭矩激勵(lì)力進(jìn)行簡(jiǎn)諧分析,可以得到如圖7-9所示的工況1-3的前32諧次的螺旋槳冰載荷總扭矩激勵(lì)力幅值。

    根據(jù)圖7-9所示的螺旋槳冰載荷總扭矩激勵(lì)力幅值分析可知:對(duì)于工況1和工況2的冰載荷下總扭矩激勵(lì)力主諧次為5諧次、10諧次和15諧次;而工況3中冰載荷下的總扭矩激勵(lì)力主諧次為10諧次、20諧次和30諧次。

    2 電機(jī)振動(dòng)階次特性分析

    電機(jī)在生產(chǎn)制造過程中存在加工誤差,使得電機(jī)磁場(chǎng)必然是非正弦分布。為研究電機(jī)的階次振動(dòng)機(jī)理以及獲得電磁轉(zhuǎn)矩的解析計(jì)算公式,建立電機(jī)A、B、C相電機(jī)物理模型,如圖10所示。

    圖7 工況1中1-32諧次的螺旋槳冰載荷總扭矩激勵(lì)力幅值Fig.7 The total torque excitation force amplitude of the propeller ice load of 1-32 harmonic times in case 1

    圖8 工況2中1-32諧次的螺旋槳冰載荷總扭矩激勵(lì)力幅值Fig.8 The total torque excitation force amplitude of the propeller ice load of 1-32 harmonic times in case 2

    圖9 工況2中1-32諧次的螺旋槳冰載荷總扭矩激勵(lì)力幅值Fig.9 The total torque excitation force amplitude of the propeller ice load of 1-32 harmonic times in case 3

    圖10 電機(jī)物理模型Fig.10 The physical model of the electric machine

    在ABC坐標(biāo)系下,根據(jù)A相、B相和C相磁鏈之間的相關(guān)關(guān)系,可以建立如下所示的三相磁鏈矩陣表達(dá)式:

    其中:τm和τ1為磁場(chǎng)分布系數(shù);θ為主磁極與A相之間的夾角;B1和Br為磁密度;為氣隙磁通密度(2i- 1 )階諧波幅值。

    根據(jù))和磁通密度)計(jì)算公式(4)進(jìn)行聯(lián)合求解,可以得到 ψa(θ)

    其中:rs為定子外圓半徑;ls為定子軸向長(zhǎng)度;

    將ABC坐標(biāo)系下的三相磁鏈矩陣表達(dá)式(3)經(jīng)過坐標(biāo)變換,可以得到在dqo坐標(biāo)系下的磁鏈方程:

    若不考慮繞組中的時(shí)間諧波電流,則在dqo坐標(biāo)系下三相對(duì)稱定子繞組中由基波正弦電流所產(chǎn)生的磁鏈為:

    其中:Ld、Lq、Lo分 別為 d、q、o 軸定子電感;id、iq、io分別為 d、q、o 軸電流。

    根據(jù)dqo坐標(biāo)系下的磁鏈方程(7)式和基波正弦電流所產(chǎn)生的磁鏈(8)式,可以得到總磁鏈為:

    在dqo坐標(biāo)系下對(duì)總磁鏈?zhǔn)剑?)進(jìn)行求導(dǎo),可以得到感應(yīng)電動(dòng)勢(shì)為:

    其中:ωr為轉(zhuǎn)子電角速度。

    根據(jù)電機(jī)的原理,電磁轉(zhuǎn)矩為:

    將(10)式和(11)式進(jìn)行聯(lián)立求解,電磁轉(zhuǎn)矩的解析式解為:

    上述電磁轉(zhuǎn)矩的解析式(12)中的第一項(xiàng)轉(zhuǎn)矩是電機(jī)在理想狀態(tài)下的電磁轉(zhuǎn)矩;但是第二項(xiàng)轉(zhuǎn)矩和第三項(xiàng)轉(zhuǎn)矩在均由磁場(chǎng)諧波引起的6k階轉(zhuǎn)矩。因此在冰載荷電力推進(jìn)軸系進(jìn)行扭轉(zhuǎn)振動(dòng)計(jì)算過程中,在電機(jī)質(zhì)量點(diǎn)處需要考慮由場(chǎng)諧波引起的6k階轉(zhuǎn)矩。

    3 電力推進(jìn)軸系計(jì)算與測(cè)試

    3.1 電力推進(jìn)軸系實(shí)例計(jì)算

    本文以雙機(jī)單槳混合動(dòng)力系統(tǒng)的PTI模式推進(jìn)軸系為研究對(duì)象,分析冰塊-螺旋槳相互作用激勵(lì)力與電機(jī)電磁激勵(lì)力對(duì)推進(jìn)軸系扭轉(zhuǎn)振動(dòng)的影響,如圖11所示為電力推進(jìn)模式示意圖。在PTI模式下,減速齒輪箱離合器H和J脫開而離合器J合排,1#大柴油機(jī)和2#小柴油機(jī)均不工作,電機(jī)經(jīng)過齒輪箱驅(qū)動(dòng)船舶運(yùn)行。采用集總參數(shù)對(duì)電力推進(jìn)軸系的部件進(jìn)行分塊建模,再將各部件扭轉(zhuǎn)振動(dòng)數(shù)學(xué)模型進(jìn)行組裝,可得如圖12所示的電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)模型和表2所示的系統(tǒng)當(dāng)量參數(shù)。

    圖11 電力推進(jìn)模式示意圖Fig.11 Schematic diagram of electric propulsion model

    圖12 電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)模型Fig.12 Torsional vibration model of electric propulsion shafting

    表2 電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)當(dāng)量參數(shù)Tab.2 Equivalent parameters of torsional vibration of electric propulsion shafting

    表3 電力模式的軸系自由振動(dòng)固有頻率Tab.3 Natural frequency of free vibration of shaft system in power mode

    采用系統(tǒng)矩陣法以及利用MATLAB軟件進(jìn)行編程,從而可以得到無(wú)阻尼自由振動(dòng)系統(tǒng)固有頻率值,如表3所示電力推進(jìn)模式的軸系自由振動(dòng)固有頻率。

    如果船舶推進(jìn)軸系的某諧次激勵(lì)力矩的變化頻率N等于軸系的某階自振頻率Nn時(shí),即產(chǎn)生第ν諧次共振;同時(shí)電磁諧波的影響將產(chǎn)生6k諧次電磁激勵(lì)力矩,因此可以得到電機(jī)各諧次激勵(lì)力矩的臨界轉(zhuǎn)速nc,各參數(shù)之間的相互關(guān)系見公式(13)所示。由于篇幅原因僅列出前9階次固有頻率下各諧次的臨界轉(zhuǎn)速,如表4所示。

    表4 1-9階次固有頻率下各諧次的臨界轉(zhuǎn)速Tab.4 Critical speed of the harmonic times of 1-9 order natural frequency

    采用系統(tǒng)矩陣法對(duì)圖15的軸系扭轉(zhuǎn)振動(dòng)數(shù)學(xué)模型進(jìn)行強(qiáng)迫振動(dòng)計(jì)算,本文僅考慮電機(jī)的階次激勵(lì)力和冰載荷激勵(lì)力與電機(jī)階次激勵(lì)力共同作用二種情況進(jìn)行理論計(jì)算,其計(jì)算結(jié)果如圖13-16所示。

    圖13 第11-12質(zhì)量點(diǎn)間(中間軸2#)扭振應(yīng)力曲線圖(僅考慮電機(jī)激勵(lì)力)Fig.13 Torsional vibration stress curves between 11st and 12nd mass points(considering only the motor excitation force)

    圖14 第14-15質(zhì)量點(diǎn)間(螺旋槳軸)扭振應(yīng)力曲線圖(僅考慮電機(jī)激勵(lì)力)Fig.14 Torsional vibration stress curves between 14th and 15th mass points(considering only the motor excitation force)

    圖15 第11-12質(zhì)量點(diǎn)(中間軸2#)應(yīng)力曲線圖(考慮電機(jī)和冰載荷工況1激勵(lì)力)Fig.15 Torsional vibration stress curves between 11st and 12nd mass points(considering the load of motor and ice load in case 1)

    圖16 第14-15質(zhì)量點(diǎn)(螺旋槳軸)應(yīng)力曲線圖(考慮電機(jī)和冰載荷工況3激勵(lì)力)Fig.16 Torsional vibration stress curves between 14th and 15th mass points(considering the load of motor and ice load in case3)

    根據(jù)圖13-16所示的電力推進(jìn)軸系強(qiáng)迫振動(dòng)計(jì)算結(jié)果可知:

    (1)在電機(jī)的階次激勵(lì)力作用下,對(duì)軸系扭轉(zhuǎn)振動(dòng)影響最大的前三諧次分別為6諧次、12諧次和18諧次。

    (2)在電機(jī)6諧次激勵(lì)力下,電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)在speed=245 r/min和speed=514 r/min處共振較為明顯,與表4所示的前9階固有頻率下各諧次的臨界轉(zhuǎn)速進(jìn)行比較,分別對(duì)應(yīng)自由振動(dòng)的3階6諧次、4階6諧次和6階6諧次。

    (3)電機(jī)6諧次激勵(lì)力下,電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)在speed=122 r/min和speed=257 r/min處共振較為明顯,與表4所示的前9階固有頻率下各諧次的臨界轉(zhuǎn)速進(jìn)行比較,分別對(duì)應(yīng)自由振動(dòng)的3階12諧次和4階12諧次。

    (4)冰載荷激勵(lì)力對(duì)軸系扭轉(zhuǎn)振動(dòng)的影響比電機(jī)階次激勵(lì)力對(duì)軸系扭轉(zhuǎn)振動(dòng)影響更大;冰載荷工況1對(duì)軸系應(yīng)力影響最大,冰載荷工況3對(duì)軸系應(yīng)力影響最小,與本文1.2節(jié)中所分析的冰載荷激勵(lì)力理論一致。

    (5)在冰載荷工況1作用下,冰載荷激勵(lì)力的主諧次5諧次、10諧次和15諧次;而在冰載荷工況3的激勵(lì)力作用下,冰載荷激勵(lì)力的主諧次10諧次、20諧次和30諧次,與本文1.2節(jié)中所分析的冰載荷激勵(lì)力理論一致。

    3.2 電力推進(jìn)軸系測(cè)試

    為了驗(yàn)證電力推進(jìn)軸系扭轉(zhuǎn)振動(dòng)計(jì)算結(jié)果的正確性,采用應(yīng)變法進(jìn)行軸系振動(dòng)測(cè)試,其測(cè)點(diǎn)分別如圖17所示的測(cè)點(diǎn)布置示意圖和圖18所示的在實(shí)船中測(cè)試儀器布置圖。測(cè)試系統(tǒng)的主要由粘貼在中間軸處的應(yīng)變片、無(wú)線接收裝置、無(wú)線應(yīng)變采集節(jié)點(diǎn)以及上位機(jī)分析系統(tǒng)。

    本次測(cè)試過程中上位機(jī)分析軟件是DHDAS動(dòng)態(tài)信號(hào)采集分析系統(tǒng),如圖19所示。應(yīng)變法軸系扭轉(zhuǎn)振動(dòng)的測(cè)試原理:

    (1)軸系在運(yùn)轉(zhuǎn)過程中軸系產(chǎn)生的扭轉(zhuǎn)振動(dòng)使得軸系產(chǎn)生微小變形,粘貼在軸系表面上的應(yīng)變片發(fā)生與之等量的變形,并通過自身電阻值的變化將軸系的形變表現(xiàn)出來(lái);

    (2)通過惠斯通電橋?qū)?yīng)變片電阻值的變化轉(zhuǎn)化為電壓信號(hào)輸出,無(wú)線動(dòng)態(tài)應(yīng)變采集節(jié)點(diǎn)通過檢測(cè)電橋輸出電壓信號(hào)的變化,獲取應(yīng)變片電阻值的變化量,從而得到軸系在運(yùn)轉(zhuǎn)過程中所產(chǎn)生的軸系應(yīng)變值;

    (3)采用無(wú)線傳輸技術(shù),將無(wú)線動(dòng)態(tài)應(yīng)變采集節(jié)點(diǎn)所采集的數(shù)據(jù)輸出至計(jì)算機(jī)系統(tǒng)中,然后通過上位機(jī)分析軟件進(jìn)行數(shù)據(jù)處理與計(jì)算,從而得到軸系的實(shí)時(shí)振動(dòng)應(yīng)力值。

    圖17 測(cè)點(diǎn)布置示意圖Fig.17 Schematic diagram of measuring point layout

    圖18 在實(shí)船中測(cè)試儀器布置圖Fig.18 Test instrument layout in a ship

    圖19 DHDAS動(dòng)態(tài)信號(hào)采集分析系統(tǒng)Fig.19 DHDASdynamic signal acquisition and analysis system

    圖20 中間軸扭轉(zhuǎn)振動(dòng)應(yīng)力曲線測(cè)試值Fig.20 Test values of intermediate shaft torsional vibration stress curve

    本次測(cè)試過程中對(duì)雙機(jī)單槳混合動(dòng)力系統(tǒng)的PTI模式推進(jìn)軸系工況進(jìn)行了測(cè)試,由于受條件約束未進(jìn)行冰區(qū)航行,因此僅獲得無(wú)冰載荷時(shí)軸系扭轉(zhuǎn)測(cè)試數(shù)據(jù)。根據(jù)上述應(yīng)變法測(cè)試的原理對(duì)測(cè)試數(shù)據(jù)進(jìn)行處理后,可以得到如圖20所示的中間軸扭轉(zhuǎn)振動(dòng)應(yīng)力曲線測(cè)試值。

    對(duì)圖13所示的中間軸扭轉(zhuǎn)振動(dòng)應(yīng)力理論計(jì)算值與圖20所示的中間軸扭轉(zhuǎn)振動(dòng)應(yīng)力測(cè)試值進(jìn)行對(duì)比分析:

    在3.1節(jié)中僅考慮電機(jī)階次激勵(lì)力的推進(jìn)軸系強(qiáng)迫振動(dòng)理論計(jì)算可知,中間軸2#在speed=245 r/min和speed=514 r/min處共振較為明顯,其應(yīng)力值分別為15.6 N/mm2和12.7 N/mm2;而與圖20所示的中間軸2#在共振點(diǎn)處的應(yīng)力值分別為14.2 N/mm2和12.1 N/mm2相比較,其誤差分別為9.2%和4.7%,理論值與測(cè)試值基本一致。

    4 結(jié) 論

    (1)本文對(duì)螺旋槳槳葉與冰塊的相互作用過程及激勵(lì)力特性進(jìn)行了研究,得到了槳葉與冰塊相互作用產(chǎn)生的激勵(lì)力時(shí)域曲線;同時(shí)對(duì)其進(jìn)行頻域分析時(shí),對(duì)于5葉螺旋槳而言工況1和工況2中冰載荷下的總扭矩激勵(lì)力主諧次為5諧次、10諧次和15諧次;而工況3中冰載荷下的總扭矩激勵(lì)力主諧次為10諧次、20諧次和30諧次。

    (2)對(duì)電機(jī)振動(dòng)機(jī)理以及電機(jī)振動(dòng)過程中闡述的階次激勵(lì)力進(jìn)行了分析,提出了電機(jī)存在階次振動(dòng)以及通過階次振動(dòng)機(jī)理和公式推導(dǎo)給出電磁轉(zhuǎn)矩的解析計(jì)算公式,為研究電力推進(jìn)軸系的扭轉(zhuǎn)振動(dòng)特性奠定了理論基礎(chǔ)。

    (3)以混合動(dòng)力推進(jìn)軸系的電動(dòng)機(jī)模式下的軸系為實(shí)例,在僅考慮電機(jī)階次激勵(lì)力作用和考慮電機(jī)階次激勵(lì)力與冰載荷激勵(lì)力作用二種情況下,運(yùn)用系統(tǒng)矩陣法對(duì)軸系的扭轉(zhuǎn)振動(dòng)特性進(jìn)行分析,驗(yàn)證了冰載荷作用下的螺旋槳激勵(lì)頻域特性和電機(jī)階次振動(dòng)相關(guān)理論的正確性,具有十分重要的工程應(yīng)用價(jià)值。

    [1]Cho SR.The resistance of ice breaking ship using the regression analysis of model tests[J].Ocean Engineering,2015,108(2015):692-703.

    [2]Kwon Y H,Lee T K,Choi K.A study on measurements of local ice pressure for ice breaking research vessel‘ARAON’at the Amundsen Sea[J].International Journal of Naval Architecture&Ocean Engineering,2015,7(3):490-499.

    [3]王傳榮.冰區(qū)LNG船未來(lái)發(fā)展方向[J].船舶物資與市場(chǎng),2010(5):15-17.Wang Chuanrong.Ice LNGtanker development direction in the future[J].Marine Equipment/Materials&Marketing,2010(5):15-17.

    [4]Tian Ximin,Zou Zaojian,Yu Jijun,Wang Fuhua.Review on advances in research of ice loads on ice-going ships[J].Journal of Ship Mechanics,2015,19(3):337-348.

    [5]Dan W,Bose N,Yamaguchi H,et al.Hydrodynamic loads on ice-class propellers during propeller-ice interaction[J].Journal of Marine Science&Technology,1997,2(1):12-20.

    [6]胡志寬.冰載荷下螺旋槳靜力分析及冰與槳碰撞動(dòng)力響應(yīng)研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2014.Hu Zhikuan.Static analysis of propeller under ice loads and dynamic response of the collision between ice and blade[D].Harbin:Harbin Institute of Technology,2014.

    [7]Barro R D,Lee D C.Excitation response estimation of polar class vessel propulsion shafting system[J].Transactions of the Korean Society for Noise and Vibration Engineering,2011,21(12):1166-1176.

    [8]江 攀.電力推進(jìn)船舶軸系振動(dòng)特性及其仿真研究[D].武漢:武漢理工大學(xué),2011.Jiang pan.Research on shafting vibration and simulation electrical marine[D].Wuhan:Wuhan University of Technology,2011.

    [9]肖能齊,周瑞平,林晞晨.冰區(qū)航行船舶電力推進(jìn)軸系機(jī)電耦合的扭振分析[J].船舶工程,2015(4):45-48.Xiao Nengqi,Zhou Ruiping,Lin Xichen.Electromechanical coupling torsional vibration analysis of electrical propulsion in icy conditions[J].Ship Engineering,2015(4):45-48.

    [10]楊紅軍,車馳東,張維競(jìng),等.冰載荷沖擊下的船舶推進(jìn)軸系瞬態(tài)扭轉(zhuǎn)振動(dòng)響應(yīng)分析[J].船舶力學(xué),2015,19(1):176-181.Yang Hongjun,Che Chidong,Zhang Weijing,Qiu Ting.Transient torsional vibration analysis for ice impact of ship propulsion shaft[J].Journal of Ship Mechanics,2015,19(1):176-181.

    猜你喜歡
    冰區(qū)階次軸系
    照亮回家的路
    我國(guó)高校首艘破冰船“中山大學(xué)極地”號(hào)成功開展冰區(qū)試航
    重覆冰區(qū)220kV雙回路窄基鋼管塔設(shè)計(jì)及試驗(yàn)研究
    吉林電力(2022年1期)2022-11-10 09:20:48
    臥式異步電機(jī)軸系支撐載荷研究
    冰區(qū)船舶壓載艙防凍方案研究
    能源工程(2022年2期)2022-05-23 13:51:44
    階次分析在驅(qū)動(dòng)橋異響中的應(yīng)用
    基于Vold-Kalman濾波的階次分析系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)*
    雙機(jī)、雙槳軸系下水前的安裝工藝
    基于齒輪階次密度優(yōu)化的變速器降噪研究
    軸系校中參數(shù)與軸系振動(dòng)特性相關(guān)性仿真研究
    色综合婷婷激情| 精品国内亚洲2022精品成人| 久久久久国内视频| 亚洲av日韩精品久久久久久密| www.色视频.com| 国产男靠女视频免费网站| 亚洲中文字幕一区二区三区有码在线看| 亚洲va在线va天堂va国产| 少妇裸体淫交视频免费看高清| 亚洲自偷自拍三级| 久久精品国产亚洲av香蕉五月| 日日啪夜夜撸| 91久久精品电影网| ponron亚洲| 欧美区成人在线视频| 在线观看免费视频日本深夜| 九九在线视频观看精品| 国产亚洲精品久久久久久毛片| 一级黄色大片毛片| 国产乱人伦免费视频| 成人国产综合亚洲| 黄色一级大片看看| 黄色日韩在线| 熟女电影av网| 99视频精品全部免费 在线| 看免费成人av毛片| 91狼人影院| 特级一级黄色大片| 搡老熟女国产l中国老女人| 日日摸夜夜添夜夜添av毛片 | 搡女人真爽免费视频火全软件 | 国产精品99久久久久久久久| 国产亚洲精品综合一区在线观看| 色综合婷婷激情| 搡老熟女国产l中国老女人| 亚洲va日本ⅴa欧美va伊人久久| 少妇熟女aⅴ在线视频| 欧美精品啪啪一区二区三区| 国产精品人妻久久久久久| 成人永久免费在线观看视频| 成人国产综合亚洲| 一区二区三区激情视频| 成年女人看的毛片在线观看| 波野结衣二区三区在线| 小说图片视频综合网站| 国产老妇女一区| 久久香蕉精品热| 在线a可以看的网站| 网址你懂的国产日韩在线| 三级男女做爰猛烈吃奶摸视频| 国产精品1区2区在线观看.| 伦理电影大哥的女人| 悠悠久久av| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩高清在线视频| 淫妇啪啪啪对白视频| 啦啦啦韩国在线观看视频| 小蜜桃在线观看免费完整版高清| 国产精品不卡视频一区二区| 成人永久免费在线观看视频| 亚洲av二区三区四区| 成年女人看的毛片在线观看| 波多野结衣高清无吗| av中文乱码字幕在线| 日韩欧美精品v在线| 久久精品国产自在天天线| 国产精品av视频在线免费观看| 亚州av有码| 国产亚洲91精品色在线| 欧美性猛交黑人性爽| 免费看美女性在线毛片视频| 国产精品女同一区二区软件 | 国产精品不卡视频一区二区| netflix在线观看网站| 日韩一区二区视频免费看| 又黄又爽又刺激的免费视频.| 人妻丰满熟妇av一区二区三区| 日韩大尺度精品在线看网址| 午夜日韩欧美国产| 国产精品一区二区三区四区免费观看 | 日本a在线网址| 亚洲五月天丁香| 好男人在线观看高清免费视频| 18禁黄网站禁片午夜丰满| 久久精品国产99精品国产亚洲性色| 在线观看舔阴道视频| 亚洲天堂国产精品一区在线| 国产精品一区二区免费欧美| 黄色欧美视频在线观看| 国产精品99久久久久久久久| 国语自产精品视频在线第100页| 日韩大尺度精品在线看网址| 中文字幕免费在线视频6| 日韩欧美在线二视频| 亚洲avbb在线观看| 国产在视频线在精品| 精品久久久久久久久久免费视频| 窝窝影院91人妻| 韩国av一区二区三区四区| 国产在线男女| 看免费成人av毛片| 国产午夜福利久久久久久| 午夜福利在线在线| 国产一区二区三区av在线 | 国产精品永久免费网站| 观看免费一级毛片| 午夜精品一区二区三区免费看| 赤兔流量卡办理| 国产探花极品一区二区| 亚洲专区国产一区二区| 999久久久精品免费观看国产| 免费看a级黄色片| 婷婷六月久久综合丁香| 国产一区二区在线av高清观看| 黄色视频,在线免费观看| 久久亚洲真实| 国产美女午夜福利| 韩国av在线不卡| 最新在线观看一区二区三区| 国产久久久一区二区三区| 99热精品在线国产| 日本五十路高清| 在线a可以看的网站| 日本免费一区二区三区高清不卡| 亚洲国产精品sss在线观看| 天堂动漫精品| 免费看美女性在线毛片视频| 99久久无色码亚洲精品果冻| 在线免费十八禁| 中文字幕免费在线视频6| 久久人妻av系列| 成人精品一区二区免费| 国产精品人妻久久久久久| 亚洲最大成人中文| 在线天堂最新版资源| 男女啪啪激烈高潮av片| 久久这里只有精品中国| 国产精品一区二区三区四区久久| 免费看美女性在线毛片视频| 天天躁日日操中文字幕| 亚洲国产色片| 亚洲人成伊人成综合网2020| 国产精品久久久久久久久免| 真实男女啪啪啪动态图| 午夜日韩欧美国产| 看十八女毛片水多多多| 亚洲中文字幕一区二区三区有码在线看| 网址你懂的国产日韩在线| 午夜福利高清视频| 美女免费视频网站| 国产国拍精品亚洲av在线观看| 偷拍熟女少妇极品色| 99久久精品国产国产毛片| 一级黄片播放器| 高清在线国产一区| 国产欧美日韩精品一区二区| 欧美激情久久久久久爽电影| 欧美一区二区精品小视频在线| 日韩精品有码人妻一区| 国产精品三级大全| 最后的刺客免费高清国语| 性欧美人与动物交配| 九色国产91popny在线| 国产伦在线观看视频一区| 熟女电影av网| 国产三级在线视频| 日韩欧美一区二区三区在线观看| 长腿黑丝高跟| 观看美女的网站| av女优亚洲男人天堂| 婷婷丁香在线五月| 99久久精品一区二区三区| 婷婷精品国产亚洲av| 99久久久亚洲精品蜜臀av| 久久精品国产自在天天线| 国产av不卡久久| 丝袜美腿在线中文| 国产精品一区www在线观看 | 神马国产精品三级电影在线观看| 国产精品av视频在线免费观看| 国产精品精品国产色婷婷| 日本一二三区视频观看| 老司机福利观看| 中文字幕av在线有码专区| 亚洲成人免费电影在线观看| 国产一区二区在线观看日韩| 最近视频中文字幕2019在线8| 日本欧美国产在线视频| 在线播放国产精品三级| 日日摸夜夜添夜夜添av毛片 | 男女下面进入的视频免费午夜| 国产黄a三级三级三级人| 性插视频无遮挡在线免费观看| 亚洲人成伊人成综合网2020| 国产高清三级在线| 搡老熟女国产l中国老女人| 动漫黄色视频在线观看| 在线观看一区二区三区| 91av网一区二区| av在线亚洲专区| 亚洲av.av天堂| 精品久久久久久久久久久久久| 亚洲成人免费电影在线观看| 久99久视频精品免费| 1024手机看黄色片| 婷婷丁香在线五月| av在线亚洲专区| 国内少妇人妻偷人精品xxx网站| 九九久久精品国产亚洲av麻豆| 99九九线精品视频在线观看视频| 99在线视频只有这里精品首页| 别揉我奶头 嗯啊视频| 欧美不卡视频在线免费观看| 日韩欧美国产一区二区入口| 88av欧美| 性插视频无遮挡在线免费观看| 婷婷丁香在线五月| 天堂av国产一区二区熟女人妻| 日韩中文字幕欧美一区二区| 91av网一区二区| 91久久精品国产一区二区三区| 日韩亚洲欧美综合| 成人二区视频| 成人综合一区亚洲| 国产精品98久久久久久宅男小说| 亚洲成人中文字幕在线播放| 欧美日韩中文字幕国产精品一区二区三区| 午夜亚洲福利在线播放| 久9热在线精品视频| 九色国产91popny在线| 亚洲最大成人中文| 美女免费视频网站| 欧美性猛交黑人性爽| 午夜爱爱视频在线播放| 乱码一卡2卡4卡精品| 亚洲av五月六月丁香网| 亚洲乱码一区二区免费版| 亚洲国产欧美人成| 久9热在线精品视频| 九色国产91popny在线| 成人国产综合亚洲| 五月伊人婷婷丁香| 乱人视频在线观看| 精品人妻熟女av久视频| 淫妇啪啪啪对白视频| 亚洲成人免费电影在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品女同一区二区软件 | 中文字幕av成人在线电影| 久久精品综合一区二区三区| 男人的好看免费观看在线视频| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产亚洲av涩爱 | 两人在一起打扑克的视频| 中文在线观看免费www的网站| 久久香蕉精品热| 99久久中文字幕三级久久日本| 在线观看午夜福利视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲,欧美,日韩| 神马国产精品三级电影在线观看| 成人特级av手机在线观看| 最后的刺客免费高清国语| 亚洲第一电影网av| 国产激情偷乱视频一区二区| 两个人的视频大全免费| 长腿黑丝高跟| 真人一进一出gif抽搐免费| 成人美女网站在线观看视频| av在线蜜桃| 麻豆国产97在线/欧美| 国产精品美女特级片免费视频播放器| 天堂av国产一区二区熟女人妻| 国产精品久久久久久亚洲av鲁大| 极品教师在线免费播放| 中出人妻视频一区二区| 久久久久国产精品人妻aⅴ院| 国产免费av片在线观看野外av| 日韩欧美精品免费久久| 午夜久久久久精精品| 人妻夜夜爽99麻豆av| 听说在线观看完整版免费高清| 美女 人体艺术 gogo| 亚洲成av人片在线播放无| 亚洲精品一区av在线观看| 欧美性感艳星| 一级av片app| 精品福利观看| 久久久久久伊人网av| 久久九九热精品免费| 中文字幕久久专区| 国产免费一级a男人的天堂| 国国产精品蜜臀av免费| 美女 人体艺术 gogo| 国内精品一区二区在线观看| 他把我摸到了高潮在线观看| 给我免费播放毛片高清在线观看| 成人一区二区视频在线观看| 国产乱人视频| 免费看美女性在线毛片视频| 成熟少妇高潮喷水视频| 在线a可以看的网站| 长腿黑丝高跟| 精品无人区乱码1区二区| 日韩欧美在线乱码| 少妇人妻一区二区三区视频| 69av精品久久久久久| 久久精品国产鲁丝片午夜精品 | 国产精品精品国产色婷婷| 午夜福利成人在线免费观看| 亚洲电影在线观看av| a级毛片a级免费在线| 蜜桃亚洲精品一区二区三区| 99视频精品全部免费 在线| 99热只有精品国产| 亚洲四区av| 看黄色毛片网站| eeuss影院久久| 国模一区二区三区四区视频| 亚洲在线自拍视频| 国产高清视频在线播放一区| 亚洲精品粉嫩美女一区| 欧美丝袜亚洲另类 | 成人毛片a级毛片在线播放| 亚洲18禁久久av| 一区二区三区四区激情视频 | 亚洲国产精品sss在线观看| 91在线精品国自产拍蜜月| 一本精品99久久精品77| 免费观看在线日韩| 国产精品三级大全| 亚洲一区高清亚洲精品| 黄色日韩在线| 变态另类丝袜制服| 久久精品国产99精品国产亚洲性色| 国内精品一区二区在线观看| 色综合亚洲欧美另类图片| 高清毛片免费观看视频网站| 大又大粗又爽又黄少妇毛片口| 无遮挡黄片免费观看| 国产精品久久视频播放| 淫妇啪啪啪对白视频| av天堂在线播放| а√天堂www在线а√下载| 亚洲avbb在线观看| 亚洲成人免费电影在线观看| 精品久久久久久久久av| 国产成人aa在线观看| 久久久久免费精品人妻一区二区| 日韩精品有码人妻一区| 国产精品嫩草影院av在线观看 | 九九在线视频观看精品| 欧美黑人欧美精品刺激| 在线看三级毛片| 波野结衣二区三区在线| 亚洲av成人av| aaaaa片日本免费| 老熟妇乱子伦视频在线观看| 亚洲av美国av| 干丝袜人妻中文字幕| 国产老妇女一区| a级毛片a级免费在线| 国产v大片淫在线免费观看| 亚洲av二区三区四区| 欧美一区二区国产精品久久精品| 国产单亲对白刺激| 国产精品亚洲一级av第二区| 国产精品日韩av在线免费观看| 久久精品国产99精品国产亚洲性色| 国产白丝娇喘喷水9色精品| 少妇的逼好多水| 夜夜夜夜夜久久久久| 亚洲经典国产精华液单| 亚洲国产欧洲综合997久久,| 一级av片app| 亚洲国产高清在线一区二区三| a级毛片a级免费在线| 亚洲不卡免费看| 午夜福利在线观看免费完整高清在 | 日本一二三区视频观看| 看片在线看免费视频| 欧美bdsm另类| 在线免费十八禁| 天天一区二区日本电影三级| 久久久久久久久大av| 国产麻豆成人av免费视频| 久久精品国产亚洲av涩爱 | 国产亚洲精品综合一区在线观看| 黄片wwwwww| av女优亚洲男人天堂| 尾随美女入室| 88av欧美| 中文字幕精品亚洲无线码一区| 麻豆成人av在线观看| 国国产精品蜜臀av免费| 高清在线国产一区| 国产男靠女视频免费网站| 久久午夜福利片| 少妇的逼水好多| 一区福利在线观看| 国产精华一区二区三区| 亚洲,欧美,日韩| 国产一区二区在线av高清观看| 极品教师在线视频| 国产一级毛片七仙女欲春2| 丝袜美腿在线中文| 少妇的逼水好多| 国产亚洲91精品色在线| 内地一区二区视频在线| 三级国产精品欧美在线观看| 免费看美女性在线毛片视频| 国产在视频线在精品| 2021天堂中文幕一二区在线观| 春色校园在线视频观看| 亚洲中文日韩欧美视频| 国产主播在线观看一区二区| 99精品久久久久人妻精品| 欧美不卡视频在线免费观看| 国产国拍精品亚洲av在线观看| 国产精品伦人一区二区| 国内少妇人妻偷人精品xxx网站| 69av精品久久久久久| 国产色婷婷99| 一区福利在线观看| bbb黄色大片| 成人特级黄色片久久久久久久| 99九九线精品视频在线观看视频| 观看免费一级毛片| 精品人妻视频免费看| 久久久色成人| 欧美激情国产日韩精品一区| 九色成人免费人妻av| 国产精品永久免费网站| 亚洲精品在线观看二区| 亚洲最大成人手机在线| 成人国产综合亚洲| 色哟哟·www| 又爽又黄a免费视频| 精品一区二区免费观看| 看十八女毛片水多多多| 最近视频中文字幕2019在线8| 国产精品一区二区免费欧美| eeuss影院久久| 精品无人区乱码1区二区| 一区二区三区激情视频| 亚洲七黄色美女视频| 人妻丰满熟妇av一区二区三区| 成人国产麻豆网| 国产真实乱freesex| .国产精品久久| 日本爱情动作片www.在线观看 | 赤兔流量卡办理| 久久久久久久久大av| 免费看日本二区| 亚洲第一区二区三区不卡| 久久久久久久精品吃奶| 91在线观看av| 欧美精品国产亚洲| 精品人妻一区二区三区麻豆 | 亚洲色图av天堂| a级毛片a级免费在线| 国内精品久久久久精免费| 高清日韩中文字幕在线| 久久久国产成人精品二区| 国内少妇人妻偷人精品xxx网站| www.www免费av| 国产毛片a区久久久久| 嫩草影院精品99| 国产精品电影一区二区三区| 99久久精品热视频| 蜜桃亚洲精品一区二区三区| 久久热精品热| 男人和女人高潮做爰伦理| 看十八女毛片水多多多| 淫秽高清视频在线观看| 免费在线观看日本一区| 亚洲av.av天堂| 极品教师在线视频| 午夜福利在线在线| 亚洲成人久久爱视频| 日本黄大片高清| 男人的好看免费观看在线视频| 成人国产综合亚洲| 免费大片18禁| 国产精品人妻久久久久久| 国产av不卡久久| 亚洲在线观看片| 国产精品免费一区二区三区在线| 99热这里只有精品一区| eeuss影院久久| 性欧美人与动物交配| 久久久午夜欧美精品| 熟女人妻精品中文字幕| 精品人妻偷拍中文字幕| a在线观看视频网站| 日韩欧美在线乱码| 国产国拍精品亚洲av在线观看| 十八禁国产超污无遮挡网站| 久久久成人免费电影| 中文字幕免费在线视频6| 亚洲精品成人久久久久久| 日本撒尿小便嘘嘘汇集6| 男人舔奶头视频| 可以在线观看的亚洲视频| 1000部很黄的大片| 日本一二三区视频观看| 婷婷亚洲欧美| 日日干狠狠操夜夜爽| 天天一区二区日本电影三级| 精品免费久久久久久久清纯| 国产伦在线观看视频一区| 欧美性猛交╳xxx乱大交人| 在线国产一区二区在线| 精品久久国产蜜桃| 精品久久久噜噜| 我的女老师完整版在线观看| 精品久久久久久久久久久久久| 久久婷婷人人爽人人干人人爱| 特大巨黑吊av在线直播| 亚洲美女黄片视频| 两人在一起打扑克的视频| 自拍偷自拍亚洲精品老妇| aaaaa片日本免费| 欧美又色又爽又黄视频| 免费在线观看成人毛片| 午夜免费成人在线视频| aaaaa片日本免费| 自拍偷自拍亚洲精品老妇| 亚洲成a人片在线一区二区| 久久久午夜欧美精品| 久久99热6这里只有精品| 日本五十路高清| 99热这里只有是精品在线观看| 国产男靠女视频免费网站| aaaaa片日本免费| 午夜影院日韩av| 欧美丝袜亚洲另类 | 一个人观看的视频www高清免费观看| 日本熟妇午夜| 一个人免费在线观看电影| 两性午夜刺激爽爽歪歪视频在线观看| 超碰av人人做人人爽久久| 内射极品少妇av片p| 日本a在线网址| 亚洲最大成人中文| 欧美激情在线99| 黄色一级大片看看| 日本 欧美在线| 欧美日韩综合久久久久久 | 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久精品电影| 内地一区二区视频在线| 国产高清激情床上av| 免费黄网站久久成人精品| 搡女人真爽免费视频火全软件 | 亚洲四区av| 人妻夜夜爽99麻豆av| 国产精品一区www在线观看 | 国产美女午夜福利| 又粗又爽又猛毛片免费看| 亚洲成人久久爱视频| 日韩欧美国产一区二区入口| 精品无人区乱码1区二区| 免费看a级黄色片| 国产单亲对白刺激| 少妇被粗大猛烈的视频| 日本在线视频免费播放| 婷婷六月久久综合丁香| 亚洲第一电影网av| 深夜精品福利| 亚洲三级黄色毛片| 欧美性猛交黑人性爽| 韩国av一区二区三区四区| 免费观看人在逋| 少妇熟女aⅴ在线视频| 男女视频在线观看网站免费| 日本成人三级电影网站| 国内精品宾馆在线| ponron亚洲| 国产中年淑女户外野战色| 精品人妻1区二区| 亚洲一区高清亚洲精品| 3wmmmm亚洲av在线观看| 别揉我奶头 嗯啊视频| 精品乱码久久久久久99久播| 免费在线观看日本一区| 成年女人永久免费观看视频| 精品久久久久久久人妻蜜臀av| 极品教师在线免费播放| 国产精品乱码一区二三区的特点| 色播亚洲综合网| 国产高清不卡午夜福利| 男人和女人高潮做爰伦理| 又粗又爽又猛毛片免费看| 夜夜夜夜夜久久久久| 国产不卡一卡二| 一本精品99久久精品77| 国产av不卡久久| 亚洲成av人片在线播放无| 狠狠狠狠99中文字幕| 国产乱人伦免费视频| 国产亚洲av嫩草精品影院| 久久热精品热| 一级黄片播放器| 男女下面进入的视频免费午夜| 天堂影院成人在线观看| 精品久久久久久久久av| 在线观看舔阴道视频| 女人十人毛片免费观看3o分钟| 精品不卡国产一区二区三区| 最近在线观看免费完整版| 波多野结衣高清无吗| 99久国产av精品| 精品一区二区免费观看| 欧美最黄视频在线播放免费|