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

    波浪推進艇的集成電機推進器設計及航速預報

    2021-06-03 08:18:04鄭躍洲王冬姣葉家瑋
    船舶力學 2021年5期
    關鍵詞:水翼推進器航速

    鄭躍洲,王冬姣,葉家瑋

    (1.華南理工大學土木與交通學院,廣州510640;2.中集海洋工程有限公司,廣東深圳518000)

    0 引 言

    波浪推進無人艇應用了無人智能技術,在海上依靠波浪能推進可實現(xiàn)遠距離的低速航行,替代有人艇進行海洋環(huán)境觀測。無人艇在靠離碼頭,或航行中遇到障礙物需要緊急避碰時,需開啟輔助動力系統(tǒng)。集成電機推進器IMP(Integrated Motor Propeller)是集電機、螺旋槳和導管于一身的新型電力推進裝置,其吊艙式設計可實現(xiàn)自由轉動,同時具有結構緊湊、安裝布置方便、可靠性高、易維護、噪音和振動小的優(yōu)點[1],是無人水面艇和AUV的理想推進器。

    集成電機推進器的設計涉及多個學科,包括電機設計和槳葉、導管的設計等。挪威科技大學Kr?vel 等人[2]研制了表貼式永磁體、吊艙式安裝的100 kW 船用集成電機推進器,通過實驗證明了其比傳統(tǒng)導管螺旋槳具有更好的水動力性能;武漢理工大學蘭加芬等[3]分析了無槳軸型集成電機推進器的槳葉在不同葉切面處厚度分布對推進器水動力性能和槳葉強度的影響。在船舶自航的數(shù)值模擬研究方面,鄭洋[4]考慮船-槳干擾問題,使用勢流渦格法計算了螺旋槳的性能,應用迭代型體積力法對集裝箱船模的自航性能進行了數(shù)值模擬,通過調節(jié)螺旋槳轉速使船模達到預計航速;吳浩[5]等通過描述型體積力法替代螺旋槳力模擬了某雙尾鰭型船模的自航,探究了虛擬盤的流入半徑、厚度和偏移量等對自航航速的影響。

    1 電機設計和空載分析

    有槳軸型集成電機推進器由永磁無刷直流電機、導管、螺旋槳和撐桿等部件組成?,F(xiàn)有波浪推進無人艇實艇的設計吃水為0.4 m,為了使推進器位于基線以上水線以下并考慮導管的尺寸,預估螺旋槳的直徑為300 mm。鑒于KA 螺旋槳的葉梢呈圓弧狀,易于與電機轉子內表面的連接,將其作為使用槳。四葉的槳葉數(shù)有利于通過CNC 整體加工得到螺旋槳和轉子,無需焊接。無人艇的目標航速為5 kn,通過螺旋槳圖譜法預取電機的轉速為1 000 r/min。由于無人艇攜帶的鋰電池為48 V,電機的功率預取4 kW。電機尺寸不大,因此采用星形繞組接法,兩相導通,三相六狀態(tài)的導通方式,雙層繞組,56 極、63 槽的分數(shù)槽和集中整距繞組方式,以及表貼式永磁體設計。通過Ansys RMxprt 模塊,使用參數(shù)化設計法確定電機定子和轉子的具體尺寸,得到電機的性能。電機主要性能為:效率89.14%,額定轉速994 r/min,額定轉矩38.4 N·m。電機安裝在導管的平行中段處,參照JD-7704 導管漿,依據(jù)轉子和定子的徑向和軸向尺寸,修改了導管的形狀,使導管內能夠容納電機。

    螺旋槳和導管共同決定集成電機推進器的水動力性能,其中螺旋槳的旋轉產生集成電機推進器的推力,導管的內表面形狀影響推力的分布和產生的流阻。推進器工作時,電機的轉矩需要大于螺旋槳和轉子的轉矩之和,轉子的轉矩可以由公式估算,依據(jù)導管槳的圖譜,選擇螺距比P/D=1.0 的KA4-55螺旋槳作為使用槳。圖1為定子線圈實物圖,圖2為安裝了螺旋槳、軸、電機轉子鐵心和N35永磁體的實物圖,圖3 為IMP 裝配圖,其中定子和轉子鐵心是由DW310-35 材質的硅鋼片通過線切割和激光焊接得到。

    圖1 定子線圈實物圖Fig.1 Stator coil

    圖2 轉子和螺旋槳實物圖Fig.2 Rotor and propeller

    圖3 IMP裝配圖Fig.3 IMP assembly drawing

    2 槳葉結構強度分析

    為了無人艇航行安全,必須保證輔助動力裝置中的螺旋槳具有足夠的強度,為此在設計集成電機推進器時須進行強度校核計算。流固耦合分析是流體力學分析和固體力學分析交叉耦合而生成的一門力學分析方法,它可分為兩種:單向流固耦合分析和雙向流固耦合分析。單向流固耦合分析是指耦合交界面處的數(shù)據(jù)傳遞是單向的,將CFD計算得到的結果傳遞給固體結構進行分析,但是沒有固體結構分析結果傳遞給流體分析的過程。雙向流固耦合分析的數(shù)據(jù)交換是雙向的,既有流體分析結果傳遞給固體結構分析,又有固體結構分析的結果反向傳遞給流體分析。

    2.1 計算方法驗證

    為驗證計算方法的有效性,通過Star ccm+軟件中雙向流固耦合的瞬態(tài)計算方法,計算與本文設計的集成電機推進器具有相同螺旋槳的導管槳槳葉的結構強度,其槳葉材質為耐海水腐蝕的5083鋁合金。設置的計算域包括螺旋槳和轉子的固體區(qū)域、旋轉域和靜態(tài)流體計算域。在旋轉域與固體區(qū)域之間創(chuàng)建流體-結構界面,在旋轉域與靜態(tài)計算域之間創(chuàng)建流體-流體界面。對旋轉設備中的流動問題可以通過創(chuàng)建一個與旋轉設備一起運動的坐標系來建模,即使用移動參考系法來處理,將轉動參考系設在與螺旋槳一起旋轉的流體區(qū)域(旋轉域)上,螺旋槳的轉速設置為994 r/min。計算時先凍結流體結構耦合求解器和固體應力求解器,待導管槳表面壓力趨于穩(wěn)定時,再解凍上述兩個求解器進行雙向流固耦合計算,最終得到導管槳在不同進速J 對應的水動力性能和槳葉結構強度,結果如圖4 和圖5所示。由圖4可知:導管槳的敞水性征曲線與圖譜試驗值[6]接近,其中在系柱(J=0)狀態(tài)下,導管槳水動力的平均誤差為0.376%;當進速系數(shù)J小于0.6時,平均誤差不到3%(平均誤差為推力系數(shù)、轉矩系數(shù)和效率誤差的平均值)。由圖5 可知:導管槳的槳葉最大應變和應力值均隨著進速系數(shù)J 的增大而減小,符合實際情況,說明本文采用的雙向流固耦合的計算方法合理可行。本文對槳葉強度校核均在系柱狀態(tài)下進行。

    圖4 導管槳敞水性征曲線Fig.4 Open water characteristic curve of ducted propeller

    圖5 槳葉最大應變和應力隨進速系數(shù)變化曲線Fig.5 Maximum strain and stress of blades as function ofadvance coefficient

    2.2 IMP槳葉強度分析

    采用與2.1節(jié)相同的方法對有槳軸型集成電機推進器在系柱狀態(tài)時進行槳葉強度計算,其中IMP的面網(wǎng)格、固體區(qū)域網(wǎng)格和網(wǎng)格剖面圖如圖6 所示,其中固體區(qū)域的網(wǎng)格數(shù)為50.8 萬,旋轉區(qū)域網(wǎng)格數(shù)為213.7萬,靜態(tài)計算域網(wǎng)格數(shù)為87.1萬。

    設置的轉速為電機額定轉速,計算得到在系柱狀態(tài)時IMP的推力約為1 113.4 N,而使用相同螺旋槳的導管槳對應的推力為1 170 N。相同進速系數(shù)和轉速時,IMP 的推力小于導管槳推力的主要原因是由于IMP的螺旋槳與電機轉子連接,轉子在導管內的旋轉產生了額外的流阻,同時螺旋槳的軸向誘導速度降低。由圖7 可知:有槳軸型集成電機推進器槳葉的導邊和隨邊與輪緣內表面連接處有應力集中,最大應力為42.38 MPa,最大應變出現(xiàn)在轉子處,最大值為0.231 mm,是槳葉直徑的0.077%。5083鋁合金的抗拉應力為270 MPa,依據(jù)公式(1),得到系柱狀態(tài)下槳葉的安全系數(shù)約為6.36。

    式中,kb為安全系數(shù),σb為強度極限應力,[ σ ]為許用應力。

    圖6 IMP相關網(wǎng)格Fig.6 IMP related grids

    圖7 系柱狀態(tài)下IMP的壓力、應力和應變云圖Fig.7 Contours of IMP pressure,stress and strain in mooring state

    在計算螺旋槳的強度時通常都是以船在全速航行時螺旋槳發(fā)出的推力及吸收的轉矩為依據(jù),考慮到船舶開航時,進速系數(shù)J 極小,從圖5 可知,此時對應的螺旋槳應力大于全速航行時的數(shù)值。另外,常規(guī)螺旋槳在工作中可能會遭遇漂浮物的碰擊而遭受突然負荷,因此常規(guī)螺旋槳所取用的安全系數(shù)較大(10 左右)。近年來國內外許多學者認為將安全系數(shù)降至8 左右為宜,對于內河船可取6 左右[7]??紤]到本文所研究的集成電機推進器,其導管可保護螺旋槳不與異物相碰,且以系柱狀態(tài)進行計算,取安全系數(shù)6 作為螺旋槳強度校核判別標準,因此可認為本文設計的IMP 槳葉強度能夠滿足螺旋槳安全系數(shù)要求。

    3 波浪推進無人艇航速預報

    3.1 方法驗證

    本文研究的波浪推進無人艇是一艘雙體船,在風浪中關閉輔助動力系統(tǒng)僅依靠安裝在船首尾位于水面以下的水翼隨艇體作縱搖和垂蕩運動產生的推力實現(xiàn)無人艇向前低速行駛[8]。在靜水中航行或遇到障礙物時才需要開啟輔助動力裝置。在無人艇進/出港遇到波浪較小的情況啟用輔助動力裝置,此時可將水翼保持在水下航行,也可以將水翼抬離水面以減小靜水阻力。通過Star ccm+重疊網(wǎng)格法計算水翼抬離水面工況下船模的靜水阻力,并與船模拖曳試驗結果進行比較,以驗證計算方法的可行性。無人艇模型的縮尺比為3,船模水線長為2.146 m,船寬為1.01 m,片體寬度為0.286 m,吃水為0.133 m,型深為0.233 m。數(shù)值模擬計算的坐標原點位于船體對稱面首部靜水面處,沿船首方向為x軸,z軸垂直向上為正。圖8所示為船模阻力計算域和船體表面及y=1.0 m剖面處的網(wǎng)格圖。

    圖8 船模靜水阻力數(shù)值模擬計算域及網(wǎng)格圖Fig.8 Numerical simulation calculation domain and grid diagram of ship model

    圖9 為船模靜水阻力及縱傾隨航速變化曲線,從圖9 可以看出:數(shù)值模擬計算結果與船模拖曳試驗數(shù)據(jù)曲線基本一致,船體總阻力和縱傾角的絕對值大致隨著航速的增加而增大(船體縱傾角負值表示尾傾)。與拖曳試驗結果比較,船模速度范圍在0.594~2.079 m/s 時,數(shù)值計算的阻力平均誤差為4.51%,最大誤差為10.28%。進一步分析可知:在1.188~1.485 m/s 和1.782~2.079 m/s 航速范圍內阻力增加較快,此時尾傾隨航速增加而增大;在航速為1.485~1.782 m/s時船體總阻力和船模尾縱傾角變化不大,說明船體的航態(tài)會影響總阻力的大小。

    圖9 船模靜水阻力及縱傾隨航速變化曲線Fig.9 Variation curve of ship model still water resistance and trim as a function of speed

    3.2 航速預報方法

    體積力法是對螺旋槳的推力和扭矩等效建模,在不實際求解螺旋槳幾何的情況下創(chuàng)建螺旋槳的推力和扭矩模型。該方法將體積力均勻地分布在圓柱形的虛擬盤體上,通過對螺旋槳盤面區(qū)域的計算單元表面施加軸向力和切向力,來模擬螺旋槳產生的推力和扭矩對流場的作用[9]。體積力fb在徑向上會發(fā)生變化,其分力的徑向分布遵循Goldstein的最佳分布,由公式(2)~(5)給出[10]:

    式中,fbx為軸向體積分力,fbθ為切向體積分力,r為徑向坐標,RH為輪轂半徑,RP為螺旋槳半徑。常數(shù)Ax和Aθ計算公式如式(6)和(7)所示:

    式中,T為推力,Q為扭矩,Δ為虛擬盤體厚度。

    體積力法有描述型和迭代型。描述型只需要輸入轉速、推力系數(shù)、扭矩系數(shù)、輪轂半徑、螺旋槳半徑、虛擬盤厚度等參數(shù),通過經驗公式計算槳葉的載荷分布,本文采用式(2)~(7)所描述的經驗公式確定體積力,適用于螺旋槳設計的初期;迭代型則通過計算實際螺旋槳的敞水性能來獲得體積力分布[9]。本文采用的是描述型,首先采用移動參考系法計算IMP 的敞水性能,即不同進速系數(shù)對應的推力系數(shù)和扭矩系數(shù),如表1 所示,再將敞水性征曲線輸入到虛擬盤計算體積力。

    3.3 實艇在靜水中航行時的航速預測

    第一代波浪推進無人艇可通過連接橋兩端的升降裝置來升降首尾水翼。水翼形狀為梯形,左右各一個,水翼本身重力與靜浮力平衡,兩端的翼型分別是NCACA-0018和NACA-0024,靜浮時水翼處于水平狀態(tài)。本節(jié)只分析實艇首部安裝固定水平翼和船后安裝一個IMP時的無人艇在靜水中自航時的航速。水翼保持在水面以下,波浪推進無人艇的面網(wǎng)格和在航速為4.4 kn 時對應的船體周圍水面分布分別如圖10和圖11所示。

    表1 集成電機推進器敞水性能Tab.1 Open water performance of IMP

    圖10 波浪推進無人艇的面網(wǎng)格Fig.10 Surface grid of waves propelling unmanned boat

    圖11 航速為4.4 kn時船體周圍水面分布圖Fig.11 Distribution of the water surface around the hull at 4.4 kn

    通過計算x方向(船長方向)凈力來獲得平衡入口流速,式(8)為x方向凈力Fx的計算公式。

    式中,θ為船體縱傾值,T為虛擬盤的總推力,Rt為船體總阻力。

    本文安裝的首固定翼在靜浮時處于水平狀態(tài),即水翼的安裝角為0°。式(8)中對應的船體縱傾角、虛擬盤總推力及船體總阻力皆為瞬時值。IMP 轉速為994 r/min,對于靜水中航行工況,取計算穩(wěn)定后的一段數(shù)據(jù)進行平均,得到圖12 和圖13 中各自數(shù)據(jù)的平均值。對于規(guī)則波中航行的船舶,其Fx值是隨遭遇周期波動的,可取計算穩(wěn)定后若干個遭遇周期內的瞬時Fx值求其平均。從圖12 可以看出:無人艇在靜水中航行時,水翼位于水下1 m 深度處,主要受粘性阻力作用,其阻力隨航速增大的曲線斜率明顯小于無人艇總阻力曲線,其中航速為4 kn 時水翼阻力約為102 N,占總阻力的26.7%。由圖13可知,縱坐標值為0的水平線與x方向凈力曲線的交點對應的橫坐標即為航速,安裝首水平固定翼和一個IMP時,無人艇在靜水中航行時的航速可達4.27 kn。

    圖12 靜水總阻力和水翼阻力隨航速變化曲線Fig.12 Total still water resistance and hydrofoil resistance versus speed

    圖13 x方向凈力隨航速變化曲線Fig.13 Net force in x direction versus speed

    3.4 實艇在迎浪規(guī)則波中航行時的運動響應及航速預測

    在3.3 節(jié)基礎上,使用五階斯托克斯波,預報無人艇在迎浪規(guī)則波中航行時的航速。模擬五階斯托克斯波的波長為15 m,水深為100 m,波高為0.4 m。坐標系x 軸從船尾指向船首方向,z 軸垂直向上,采用右手坐標系,坐標原點位于縱中剖面船首靜水面處,x=-3.39 m 為無人艇重心所在位置的縱向坐標值。在x=-3.39 m,y=10 m 處設置波面位移監(jiān)測點,目的是通過分析垂蕩、縱搖運動和波面位移的時歷曲線,得到運動垂蕩和縱搖運動的頻率響應函數(shù)及相位角。圖14 為無人艇航速是4 kn時,與船舶重心縱向位置相同點的波面位移時歷曲線。對x=-3.39 m,y=10 m處,t=80 s前八個周期的波進行統(tǒng)計,得到的平均波高為0.358 2 m,約為設置波高的89.55%,遭遇周期為2.17 s。

    圖14 航速4 kn時波面位移時歷曲線(x=-3.39 m,y=10 m)Fig.14 Time history of wave elevation at 4 kn(x=-3.39 m,y=10 m)

    圖15 無人艇在迎浪規(guī)則波中航行時的時歷曲線Fig.15 Time-history curve of unmanned boat travelling in regular head waves

    圖15結果顯示無人艇在迎浪規(guī)則波中,以4 kn的航速航行時的水翼受力、總阻力、虛擬盤推力的時歷曲線和縱搖運動響應均呈周期性變化。在一個波浪周期內,作用在水翼上的水平力有2 個峰值和2 個谷值,即在一個波浪周期內出現(xiàn)了二次周期性循環(huán),阻力和推力交替出現(xiàn),兩個峰值和谷值大小不等,水翼平均值為-13.2 N,負值表示推力;艇體、水翼及立柱在內的總阻力平均值是556.8 N,虛擬盤的平均推力為521.8 N,平均縱傾約為-1.46°,x方向的平均凈力約為-35.2 N。因此基于首水平固定翼的波浪推進無人艇,在波長為15 m、波高約為0.36 m 的迎浪中航行時,航速預計略小于4 kn。船體縱搖雙幅值為0.16°,可得縱搖頻率響應函數(shù)為0.454°/m,縱搖與波面位移間的相位差為-152.6°。

    圖16 和圖17 分別為波浪推進艇以4 kn 航速在迎浪規(guī)則波中航行時垂蕩運動響應和t=80 s 時刻的波形圖??芍故庍\動雙幅值為0.294 m,垂蕩頻率響應函數(shù)為0.82,與波面曲線對比可知,垂蕩運動與波面位移之間的相位差約為21.6°。

    圖16 垂蕩運動響應Fig.16 Heave motion response

    圖17 以4 kn航速航行時的波形圖(t=80 s)Fig.17 Free surface wave pattern at 4 kn(t=80 s)

    4 結 論

    本文通過等效磁路法設計了集成電機推進器所用的無刷直流電機,使用雙向流固耦合方法校核了槳葉的結構強度,預測了含首固定水平翼之波浪推進無人艇在安裝設計的集成電機推進器后分別在靜水中和迎浪規(guī)則波中航行時的航速。研究結果表明:有槳軸型集成電機推進器的槳葉最大應力發(fā)生在葉梢與輪緣內表面的連接處,最大應變發(fā)生在轉子外表面;描述型體積力法替代螺旋槳的作用能便捷且較為快速地預測航行器的航速,便于設計推進器的安裝位置;波浪推進無人艇在靜水中航行時水翼會產生較大阻力,而在波浪作用下,水翼隨艇體搖蕩,可產生推力。

    猜你喜歡
    水翼推進器航速
    VLCC在波浪中的航速優(yōu)化與能效優(yōu)化分析
    提升全回轉港作拖輪航速的有效途徑
    水上消防(2022年1期)2022-06-16 08:06:56
    波浪滑翔機橢圓形后緣水翼動力特性研究
    基于CFD扇翼推進器敞水性能預報分析
    袖珍水翼突防潛艇的設計構想及運用研究
    發(fā)揮考核“指揮棒”“推進器”作用
    當代陜西(2018年12期)2018-08-04 05:49:08
    低速水面目標航速精度分析及精確解算
    三維扭曲水翼空化現(xiàn)象CFD模擬
    基于CFD的波浪滑翔機航速預測
    讓黨建成為信仰播種機和工作推進器
    學習月刊(2015年10期)2015-07-09 03:35:12
    免费av毛片视频| 日韩制服骚丝袜av| 亚洲精品亚洲一区二区| 精品一区二区三区人妻视频| 国产熟女欧美一区二区| 全区人妻精品视频| 99视频精品全部免费 在线| 我的老师免费观看完整版| 免费观看在线日韩| 麻豆久久精品国产亚洲av| freevideosex欧美| 亚洲怡红院男人天堂| 免费大片18禁| 三级国产精品片| 国产激情偷乱视频一区二区| 男人舔奶头视频| 2018国产大陆天天弄谢| 久99久视频精品免费| 在线观看人妻少妇| 亚洲欧美精品自产自拍| 精品久久久久久久久av| 免费黄色在线免费观看| 日韩成人伦理影院| 国产极品天堂在线| 校园人妻丝袜中文字幕| 中文字幕免费在线视频6| 欧美激情国产日韩精品一区| 男的添女的下面高潮视频| 国产成年人精品一区二区| 亚洲欧美清纯卡通| 久久精品熟女亚洲av麻豆精品 | 一级a做视频免费观看| 亚洲国产最新在线播放| 十八禁国产超污无遮挡网站| 一级毛片黄色毛片免费观看视频| 亚洲18禁久久av| 一本久久精品| 中文字幕久久专区| 国产乱人视频| 青春草亚洲视频在线观看| 男人爽女人下面视频在线观看| 国产高清三级在线| 欧美三级亚洲精品| 亚洲国产欧美人成| 亚洲精品中文字幕在线视频 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品嫩草影院av在线观看| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| 国产成年人精品一区二区| 国产亚洲av嫩草精品影院| 精品久久久久久久久久久久久| 亚洲av二区三区四区| 热99在线观看视频| 免费观看的影片在线观看| 丝袜喷水一区| 日本黄大片高清| 国产成年人精品一区二区| 丰满少妇做爰视频| 99热6这里只有精品| 99久久精品一区二区三区| 国产精品综合久久久久久久免费| 我要看日韩黄色一级片| 中文字幕免费在线视频6| 久久久久久久久久久免费av| 免费电影在线观看免费观看| 乱系列少妇在线播放| 日韩一区二区视频免费看| av网站免费在线观看视频 | 午夜久久久久精精品| 国产男女超爽视频在线观看| 久久久久精品久久久久真实原创| 精品久久久久久久末码| 精品一区二区三卡| 2022亚洲国产成人精品| 免费观看性生交大片5| 亚洲国产成人一精品久久久| 成年女人看的毛片在线观看| 综合色av麻豆| 久久久久久久国产电影| 永久网站在线| or卡值多少钱| 蜜桃久久精品国产亚洲av| 三级经典国产精品| 精品午夜福利在线看| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验| av免费观看日本| 久久97久久精品| 国产91av在线免费观看| 寂寞人妻少妇视频99o| 国产精品人妻久久久影院| 成人午夜精彩视频在线观看| av播播在线观看一区| 少妇裸体淫交视频免费看高清| 久久草成人影院| 欧美不卡视频在线免费观看| 五月玫瑰六月丁香| 在线免费十八禁| 午夜免费男女啪啪视频观看| 亚洲国产精品sss在线观看| 一个人免费在线观看电影| 午夜福利成人在线免费观看| 99视频精品全部免费 在线| 亚洲av免费在线观看| 亚洲国产色片| 国产成人freesex在线| 少妇被粗大猛烈的视频| 最近手机中文字幕大全| 午夜视频国产福利| 91av网一区二区| 丝袜美腿在线中文| 国产免费又黄又爽又色| 中文欧美无线码| 狂野欧美激情性xxxx在线观看| 天堂俺去俺来也www色官网 | 五月伊人婷婷丁香| 欧美精品国产亚洲| 高清午夜精品一区二区三区| 欧美极品一区二区三区四区| 伊人久久国产一区二区| 91精品一卡2卡3卡4卡| 麻豆成人av视频| 欧美xxⅹ黑人| 日韩伦理黄色片| 最近视频中文字幕2019在线8| 春色校园在线视频观看| 又爽又黄无遮挡网站| 久久久精品免费免费高清| 亚洲成色77777| 国产一区二区三区av在线| 久久人人爽人人爽人人片va| 插阴视频在线观看视频| 99久国产av精品| 中国美白少妇内射xxxbb| 身体一侧抽搐| 久久久久久久久大av| 美女xxoo啪啪120秒动态图| 国产老妇伦熟女老妇高清| 伊人久久国产一区二区| www.av在线官网国产| 亚洲精品久久久久久婷婷小说| 国产精品日韩av在线免费观看| 国产在视频线精品| 欧美成人午夜免费资源| 高清视频免费观看一区二区 | 日韩一区二区三区影片| 国内精品一区二区在线观看| 国产精品一及| av国产久精品久网站免费入址| 伊人久久国产一区二区| 中国美白少妇内射xxxbb| av在线观看视频网站免费| 一级片'在线观看视频| 国产精品久久视频播放| 男人狂女人下面高潮的视频| 日韩成人伦理影院| 亚洲精品aⅴ在线观看| 午夜激情久久久久久久| 国产人妻一区二区三区在| 亚洲丝袜综合中文字幕| 麻豆久久精品国产亚洲av| 亚洲自偷自拍三级| 国模一区二区三区四区视频| 99久久人妻综合| 伊人久久精品亚洲午夜| 免费观看a级毛片全部| 国产成人精品婷婷| 18禁在线无遮挡免费观看视频| 久久草成人影院| 欧美日本视频| 91午夜精品亚洲一区二区三区| 国产成人免费观看mmmm| 丰满少妇做爰视频| 中文字幕制服av| 免费看光身美女| 亚洲精品,欧美精品| a级毛片免费高清观看在线播放| 午夜精品在线福利| 国产日韩欧美在线精品| 亚洲第一区二区三区不卡| 在线免费观看的www视频| 插阴视频在线观看视频| 亚洲精品一区蜜桃| 少妇高潮的动态图| 国产 一区精品| 亚洲精品久久午夜乱码| 干丝袜人妻中文字幕| 少妇高潮的动态图| 国产黄a三级三级三级人| 亚洲一区高清亚洲精品| 国产免费一级a男人的天堂| 在线观看美女被高潮喷水网站| 国产久久久一区二区三区| 亚洲精品aⅴ在线观看| 一夜夜www| 99久国产av精品| 国内精品美女久久久久久| 国产综合精华液| 蜜桃久久精品国产亚洲av| 成人毛片a级毛片在线播放| 国产精品女同一区二区软件| 成人综合一区亚洲| 日韩不卡一区二区三区视频在线| 美女高潮的动态| 精品熟女少妇av免费看| 免费大片18禁| 日韩电影二区| 亚洲av男天堂| 国产精品伦人一区二区| 免费观看在线日韩| 啦啦啦中文免费视频观看日本| 亚洲av成人精品一区久久| av在线观看视频网站免费| 一区二区三区乱码不卡18| 国产精品不卡视频一区二区| 亚洲伊人久久精品综合| 麻豆乱淫一区二区| 插逼视频在线观看| 国精品久久久久久国模美| 联通29元200g的流量卡| 亚洲精品成人久久久久久| 天天躁日日操中文字幕| 国产 一区 欧美 日韩| 精品不卡国产一区二区三区| 在线观看一区二区三区| 久久99热这里只频精品6学生| 99久国产av精品| 久久久久久久久久久免费av| 日日撸夜夜添| 插阴视频在线观看视频| 亚洲av成人精品一区久久| 国产精品一区二区性色av| 国国产精品蜜臀av免费| 国产精品美女特级片免费视频播放器| 国产伦精品一区二区三区四那| 日韩精品青青久久久久久| 搡老乐熟女国产| 国产精品久久久久久久久免| 99久久九九国产精品国产免费| 99热网站在线观看| 日本一二三区视频观看| 免费人成在线观看视频色| 欧美三级亚洲精品| 夫妻午夜视频| 亚洲一区高清亚洲精品| 三级经典国产精品| 久久久精品欧美日韩精品| 中文在线观看免费www的网站| 日韩欧美精品v在线| 精品欧美国产一区二区三| 白带黄色成豆腐渣| 欧美一区二区亚洲| 啦啦啦韩国在线观看视频| 99久久精品一区二区三区| 在线观看一区二区三区| 亚洲成人一二三区av| 永久免费av网站大全| 亚洲婷婷狠狠爱综合网| 成人高潮视频无遮挡免费网站| 亚洲国产欧美人成| 精品国内亚洲2022精品成人| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成年免费大片在线观看| 久久精品人妻少妇| 乱系列少妇在线播放| 美女脱内裤让男人舔精品视频| 免费大片18禁| 免费看不卡的av| 美女黄网站色视频| av线在线观看网站| 欧美bdsm另类| 亚洲精品456在线播放app| 日韩一本色道免费dvd| 亚洲精品国产av成人精品| 久久鲁丝午夜福利片| 综合色av麻豆| 国产在线男女| 春色校园在线视频观看| 日韩视频在线欧美| 精品国产一区二区三区久久久樱花 | 久99久视频精品免费| 人妻系列 视频| 亚洲色图av天堂| 美女黄网站色视频| 日本黄色片子视频| 肉色欧美久久久久久久蜜桃 | 亚洲国产欧美在线一区| 干丝袜人妻中文字幕| 亚洲国产高清在线一区二区三| 噜噜噜噜噜久久久久久91| 久久久久九九精品影院| 国产又色又爽无遮挡免| 日本与韩国留学比较| 一级毛片久久久久久久久女| 亚洲美女搞黄在线观看| 免费观看av网站的网址| 精品一区在线观看国产| 久久久久性生活片| 伦理电影大哥的女人| 在线免费观看的www视频| 国产成人aa在线观看| 国产精品99久久久久久久久| av免费在线看不卡| 成人性生交大片免费视频hd| 国产成人一区二区在线| 菩萨蛮人人尽说江南好唐韦庄| 在线免费十八禁| 成人高潮视频无遮挡免费网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美性猛交╳xxx乱大交人| 成人亚洲精品av一区二区| 插逼视频在线观看| 综合色av麻豆| 午夜久久久久精精品| 亚洲精品亚洲一区二区| 97热精品久久久久久| 97超碰精品成人国产| 国内揄拍国产精品人妻在线| 水蜜桃什么品种好| 国产精品.久久久| 日韩强制内射视频| 九九爱精品视频在线观看| 国产午夜精品久久久久久一区二区三区| 能在线免费观看的黄片| 国产乱人视频| 亚洲精品成人av观看孕妇| 午夜福利高清视频| 亚洲人与动物交配视频| 日日啪夜夜撸| 麻豆精品久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 精品一区二区三区视频在线| 1000部很黄的大片| 久久久色成人| 成年女人看的毛片在线观看| 国产精品嫩草影院av在线观看| 日本免费在线观看一区| 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 午夜福利视频精品| 三级男女做爰猛烈吃奶摸视频| 少妇猛男粗大的猛烈进出视频 | 国内揄拍国产精品人妻在线| 国产久久久一区二区三区| 男人舔女人下体高潮全视频| 亚洲不卡免费看| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 精品久久久久久电影网| 亚洲av电影不卡..在线观看| 日韩伦理黄色片| 18禁动态无遮挡网站| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲网站| 免费看a级黄色片| 少妇熟女欧美另类| 国产亚洲5aaaaa淫片| 免费在线观看成人毛片| 中文欧美无线码| 日本黄色片子视频| 五月伊人婷婷丁香| 欧美高清性xxxxhd video| 26uuu在线亚洲综合色| 成年免费大片在线观看| av国产免费在线观看| 蜜桃亚洲精品一区二区三区| 日韩精品有码人妻一区| 精品亚洲乱码少妇综合久久| 色吧在线观看| 午夜精品一区二区三区免费看| 亚洲自偷自拍三级| 秋霞伦理黄片| 亚洲最大成人av| 国产精品人妻久久久久久| or卡值多少钱| 伦精品一区二区三区| 日本免费a在线| 欧美人与善性xxx| 午夜亚洲福利在线播放| 看免费成人av毛片| 69av精品久久久久久| 欧美不卡视频在线免费观看| 日韩欧美精品免费久久| 国产精品一区二区性色av| 亚洲av免费高清在线观看| 精品一区二区三区人妻视频| 51国产日韩欧美| 午夜日本视频在线| 亚洲精品中文字幕在线视频 | 欧美成人一区二区免费高清观看| 国产老妇女一区| 能在线免费看毛片的网站| 成人亚洲欧美一区二区av| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久人人人人人人| 午夜免费激情av| 人体艺术视频欧美日本| 18禁裸乳无遮挡免费网站照片| 全区人妻精品视频| 久久午夜福利片| 最近中文字幕2019免费版| 男插女下体视频免费在线播放| 欧美三级亚洲精品| 亚洲人与动物交配视频| 日韩欧美精品v在线| 女的被弄到高潮叫床怎么办| 精品一区二区免费观看| 国产精品麻豆人妻色哟哟久久 | 亚洲欧美成人精品一区二区| 校园人妻丝袜中文字幕| av福利片在线观看| 中文乱码字字幕精品一区二区三区 | 91在线精品国自产拍蜜月| 精品酒店卫生间| 一区二区三区免费毛片| 日本黄大片高清| 国产美女午夜福利| 日日摸夜夜添夜夜添av毛片| or卡值多少钱| 黄色一级大片看看| 一级毛片我不卡| 寂寞人妻少妇视频99o| 久99久视频精品免费| 亚洲精品久久久久久婷婷小说| 免费看av在线观看网站| 亚洲怡红院男人天堂| 久久久a久久爽久久v久久| 午夜激情久久久久久久| 久久亚洲国产成人精品v| 久久99热这里只频精品6学生| 国产伦在线观看视频一区| 久久精品久久精品一区二区三区| 三级男女做爰猛烈吃奶摸视频| 大香蕉97超碰在线| 日韩av不卡免费在线播放| 午夜激情福利司机影院| 自拍偷自拍亚洲精品老妇| 女的被弄到高潮叫床怎么办| 国产 亚洲一区二区三区 | 又黄又爽又刺激的免费视频.| 久久人人爽人人爽人人片va| 尾随美女入室| 嫩草影院入口| 国产一级毛片七仙女欲春2| 国产黄色视频一区二区在线观看| 国产亚洲5aaaaa淫片| 99久国产av精品| 99热网站在线观看| 极品少妇高潮喷水抽搐| 久久精品夜夜夜夜夜久久蜜豆| 免费看日本二区| 亚洲精品,欧美精品| 国产黄色视频一区二区在线观看| 亚洲最大成人手机在线| 99久久精品一区二区三区| 淫秽高清视频在线观看| 尾随美女入室| 亚洲婷婷狠狠爱综合网| 91久久精品电影网| 99热这里只有是精品在线观看| 欧美bdsm另类| 成人美女网站在线观看视频| 国产精品精品国产色婷婷| 日韩,欧美,国产一区二区三区| 性插视频无遮挡在线免费观看| 又粗又硬又长又爽又黄的视频| 国产激情偷乱视频一区二区| 啦啦啦中文免费视频观看日本| h日本视频在线播放| 国内精品宾馆在线| 国产乱人视频| 能在线免费观看的黄片| 亚洲精品aⅴ在线观看| 午夜精品在线福利| 一级毛片黄色毛片免费观看视频| 午夜老司机福利剧场| av福利片在线观看| 男女视频在线观看网站免费| 国产精品一区二区三区四区免费观看| 亚洲精品久久午夜乱码| 国产伦在线观看视频一区| 超碰av人人做人人爽久久| 亚洲性久久影院| 搡老妇女老女人老熟妇| 精品少妇黑人巨大在线播放| 乱人视频在线观看| 天堂av国产一区二区熟女人妻| 亚洲18禁久久av| 99热这里只有是精品50| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产亚洲精品久久久com| 亚洲成色77777| 亚洲精品影视一区二区三区av| 国产黄片美女视频| 色播亚洲综合网| www.av在线官网国产| 69人妻影院| 人妻制服诱惑在线中文字幕| 午夜亚洲福利在线播放| 夫妻性生交免费视频一级片| 最近2019中文字幕mv第一页| 26uuu在线亚洲综合色| 99热这里只有是精品50| 青青草视频在线视频观看| 老司机影院毛片| 亚洲成人久久爱视频| 3wmmmm亚洲av在线观看| 国产探花在线观看一区二区| 日日干狠狠操夜夜爽| 亚洲欧美精品自产自拍| 国产单亲对白刺激| 在线免费观看的www视频| 中文资源天堂在线| 好男人视频免费观看在线| 在线观看av片永久免费下载| 亚洲久久久久久中文字幕| 少妇的逼好多水| 免费大片黄手机在线观看| 免费观看的影片在线观看| 人人妻人人澡欧美一区二区| 极品少妇高潮喷水抽搐| 青春草亚洲视频在线观看| 大陆偷拍与自拍| 51国产日韩欧美| 少妇人妻一区二区三区视频| 国产日韩欧美在线精品| 搡老妇女老女人老熟妇| 我的老师免费观看完整版| 蜜桃亚洲精品一区二区三区| 国产亚洲91精品色在线| 亚洲欧美日韩东京热| 免费观看的影片在线观看| 黄片无遮挡物在线观看| 老师上课跳d突然被开到最大视频| 男女边吃奶边做爰视频| 亚洲经典国产精华液单| 少妇人妻精品综合一区二区| 婷婷六月久久综合丁香| 国产精品人妻久久久久久| 欧美成人一区二区免费高清观看| 2022亚洲国产成人精品| 中文在线观看免费www的网站| 乱人视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 欧美成人一区二区免费高清观看| 国产黄色免费在线视频| videos熟女内射| 乱人视频在线观看| 免费观看av网站的网址| 嘟嘟电影网在线观看| 中文精品一卡2卡3卡4更新| 亚洲在久久综合| 好男人在线观看高清免费视频| 欧美日韩精品成人综合77777| 欧美97在线视频| 69av精品久久久久久| 日韩一本色道免费dvd| 色尼玛亚洲综合影院| 嫩草影院精品99| 你懂的网址亚洲精品在线观看| 99热这里只有是精品在线观看| 国产 亚洲一区二区三区 | 亚洲精品一二三| 亚洲欧美清纯卡通| 91av网一区二区| 欧美潮喷喷水| videossex国产| 久久99蜜桃精品久久| 久久久久久久亚洲中文字幕| 观看美女的网站| 国产免费一级a男人的天堂| 亚洲精品国产av蜜桃| 日本猛色少妇xxxxx猛交久久| 亚洲成色77777| 建设人人有责人人尽责人人享有的 | 国产精品人妻久久久久久| 18禁裸乳无遮挡免费网站照片| 国产精品无大码| 国产淫语在线视频| 亚洲精品影视一区二区三区av| 国产亚洲5aaaaa淫片| 亚洲在线自拍视频| av在线播放精品| 国产永久视频网站| 亚洲av在线观看美女高潮| 嫩草影院精品99| 寂寞人妻少妇视频99o| 国产91av在线免费观看| 激情 狠狠 欧美| 亚洲欧美精品自产自拍| 熟妇人妻不卡中文字幕| 亚洲av电影在线观看一区二区三区 | 永久网站在线| 少妇熟女欧美另类| 日韩精品有码人妻一区| 久久久精品94久久精品| 99久久精品一区二区三区| 有码 亚洲区| 国产成人一区二区在线| 欧美激情在线99| 亚洲国产高清在线一区二区三| 亚洲精品第二区| 亚洲自拍偷在线| 99久国产av精品| 国产欧美另类精品又又久久亚洲欧美| 亚洲人成网站在线观看播放| 少妇裸体淫交视频免费看高清| 九草在线视频观看| 国产亚洲精品久久久com| 久久99精品国语久久久| 久久精品国产亚洲av天美| 成人国产麻豆网| 欧美一级a爱片免费观看看| 最近视频中文字幕2019在线8|