999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

船舶在風(fēng)浪中的加速能力評估方法

2023-12-27 01:48:40王國強孫海曉蔣曙暉
船海工程 2023年6期
關(guān)鍵詞:船舶能力

王國強,孫海曉,蔣曙暉

(上海船舶研究設(shè)計院,上海 201203)

內(nèi)貿(mào)集裝箱船由于航程短,航速普遍不高,裝機功率儲備不會太大,惡劣海況下用于加速的功率不多,加速能力問題就更為凸顯,所以,有必要對船舶在海上風(fēng)浪中的加速能力進行分析和評估,以保證船舶在更為復(fù)雜的海況下正常營運。為此,考慮從肥大型內(nèi)貿(mào)集裝箱船特點出發(fā),結(jié)合螺旋槳水動力特性,對海上風(fēng)浪中的船舶加速能力的相關(guān)內(nèi)容進行分析,形成能被實船應(yīng)用的評估方法。

1 螺旋槳水動力特性與實船模型

螺旋槳的水動力性能是指:一定幾何形體的螺旋槳在水中運動時所產(chǎn)生的推力、消耗的轉(zhuǎn)矩和效率與其運動(進速Va和轉(zhuǎn)速n)間的關(guān)系[1]。

推力系數(shù)為

(1)

轉(zhuǎn)矩系數(shù)為

(2)

式中:T為推力,kN;ρ為水密度,t/m3;n為螺旋槳轉(zhuǎn)速,r/s;D為螺旋槳直徑,m;Q為轉(zhuǎn)矩,kN/m。

當(dāng)螺旋槳以進速Va和轉(zhuǎn)速n工作時,必須吸收主機所供給的轉(zhuǎn)矩Q才能發(fā)出推力T,其所吸收的主機功率為

PME=2πnQ

(3)

螺旋槳進速系數(shù)J為

(4)

對于幾何形狀一定的螺旋槳而言,kt、kq及效率η0僅與J有關(guān),船舶航行加速能力與螺旋槳密切相關(guān)。根據(jù)螺旋槳敞水特性,kt、kq在J=0時達到最大,所以,螺旋槳推力極限即船舶推力極限的研究一般把重點放在J=0的狀態(tài)下。

1)正常航行狀態(tài)螺旋槳某半徑處葉原體受力見圖1[2],當(dāng)推力為正值時,螺旋槳尚能克服船舶部分或全部阻力,為推動船舶前行提供推力。見式(4),對槳徑D已確定的螺旋槳,其J隨Va和n變化。見圖1,當(dāng)轉(zhuǎn)速保持不變,隨著Va(亦即J值)的增加,攻角αk隨之減小,從而推力和轉(zhuǎn)矩也相應(yīng)減小。

圖1 槳葉原體受力示意

2)假定一種極端情況,當(dāng)船舶在沒有洋流的高海況下航行,主機保持服務(wù)功率輸出不變,海況高到船舶服務(wù)航速無法維持,航速不斷降低直至為0,主機輸出的功率全部用來抵抗風(fēng)浪造成的增阻,螺旋槳有效推力與船舶遭受的阻力達到平衡,此時船舶就相當(dāng)于受到一個系柱力作用而原地不動,這種極端情況可以簡化為螺旋槳系柱狀態(tài)。

系柱狀態(tài)下,螺旋槳進速系數(shù)J=0,螺旋槳只旋轉(zhuǎn)而不前進,其速度和力的關(guān)系就由船舶航行狀態(tài)下螺旋槳受力圖1轉(zhuǎn)變?yōu)槭芰D2,升力dL將與推力dT重合,槳葉上各葉原體具有最大的攻角αk,所以推力T和轉(zhuǎn)矩Q都達到某個轉(zhuǎn)速下的最大值。

圖2 槳葉原體受力示意

這種情況下,倘若螺旋槳仍然有條件增大轉(zhuǎn)矩運轉(zhuǎn),則船舶推力等于阻力的平衡就可以被打破,主機可以通過增加功率輸出提高螺旋槳轉(zhuǎn)矩,槳葉原體受力就由圖2逐漸轉(zhuǎn)變?yōu)閳D1,攻角αk縮小,水動力螺距角β增大,進程P逐步增加,船舶在這種極端海況下就具有了正值推力,也就具有了加速能力,從而可以在惡劣海況下逐步加速前進,進而脫困。換言之,倘若這種情況下主機再無功率余量提供給螺旋槳,已達功率輸出限制,則這種海上的“僵持”狀態(tài)就難以打破,船舶脫困就只有期待海況降低增阻減小,十分被動。所以,在螺旋槳某個轉(zhuǎn)速下的系柱狀態(tài),其旋轉(zhuǎn)功率相對主機在該轉(zhuǎn)速下的限制功率的余量大小就成為船舶加速能力大小的一個體現(xiàn)。

3)另一種情況是船舶正常航行中遭遇了3級、4級或更高海況,并未處于J=0的系柱狀態(tài),仍然能夠以一定的航速前進,倘若需要加速前進,就需要看主機在常用轉(zhuǎn)速下是否還有轉(zhuǎn)矩余量傳輸給螺旋槳,所以該種情況的分析就成為螺旋槳帶海況的輕載功率在常用轉(zhuǎn)速下距離主機轉(zhuǎn)矩限制是否仍然有余量的分析。

2 螺旋槳重載工況下加速能力的影響因素

船舶設(shè)計方案確定后,螺旋槳可獲得數(shù)據(jù)kt、kq、n、D。主機可獲得數(shù)據(jù)有額定功率MCR、額定轉(zhuǎn)速N、轉(zhuǎn)速禁區(qū)BSR。通過這些已知數(shù)據(jù)計算當(dāng)螺旋槳進速系數(shù)J=0時的系柱狀態(tài)下,螺旋槳旋轉(zhuǎn)產(chǎn)生的最大轉(zhuǎn)矩Qmax。根據(jù)螺旋槳特性式(2)中的參數(shù)變化如下。

(5)

(6)

由此得到系柱螺旋槳最大轉(zhuǎn)速為

(7)

式中:kq為螺旋槳轉(zhuǎn)矩系數(shù);N為主機額定轉(zhuǎn)速,r/min。

系柱拖力用螺旋槳收到功率P表示,則系柱拖力曲線成為螺旋槳收到功率P關(guān)于螺旋槳轉(zhuǎn)速n的三次冪函數(shù):

P(n)=2πkqρD5·n3(0

(8)

繪制函數(shù)P(n)曲線見圖3,為系柱拖力曲線。

圖3 系柱拖力對比

當(dāng)系柱狀態(tài)下螺旋槳收到最大功率Pmax的時候,主機軸上輸出最大轉(zhuǎn)矩Qmax,對應(yīng)的螺旋槳重載轉(zhuǎn)速也為最大nmax。從曲線變化可得到如下特征:①相同的螺旋槳收到功率,對應(yīng)的轉(zhuǎn)速越小,則P(n)曲線越陡峭;②同一轉(zhuǎn)速下螺旋槳遭受的阻力矩越大,P(n)曲線越陡峭,表明船舶在極端風(fēng)浪中的加速能力就越差。

如圖3所示,當(dāng)螺旋槳收到功率Pmax的時候,如果對應(yīng)主機轉(zhuǎn)速n1max,則系柱拖力曲線為P1(n),該曲線對應(yīng)一個確定的設(shè)計槳。但如果對設(shè)計槳進行一些修改,增加輕轉(zhuǎn)裕度,則會得到系柱拖力曲線P2(n)。很顯然曲線P2(n)比曲線P1(n)更平緩一些,輕轉(zhuǎn)裕度越大,拖力曲線越平緩,則在同一系柱轉(zhuǎn)速下所對應(yīng)的主機功率P就越小。從式(8)P(n)函數(shù)可以看到, 影響P值大小的參數(shù)有kq、D、n。在主機不變、轉(zhuǎn)速n一定的情況下,通過調(diào)整螺旋槳設(shè)計,改變kq、D這2個參數(shù)能達到減小功率曲線P(n)斜率增幅的目標(biāo),使得功率P距離主機轉(zhuǎn)矩限制線更遠,從而獲得更大的轉(zhuǎn)矩余量,船舶在風(fēng)浪中的加速能力就更好。

同時,從圖3也可以看到,當(dāng)主機一定,螺旋槳收到功率P一定時,通過提高對應(yīng)的螺旋槳轉(zhuǎn)速n1至n2,則系柱拖力曲線也會變得更為平緩,整個曲線距離主機轉(zhuǎn)矩限制線的距離也就越遠,在螺旋槳同一轉(zhuǎn)速下船舶的加速能力就越強,遇到大風(fēng)大浪脫困能力也就越強。

綜上所述,船舶加速能力可以通過螺旋槳在同一系柱轉(zhuǎn)速下收到功率距離主機轉(zhuǎn)矩限制線之間的距離來體現(xiàn),距離越大,加速能力越強,距離越小,加速能力越弱。所以在船舶設(shè)計階段,可以通過調(diào)節(jié)這一功率差值改變船舶加速能力,方式如下:①螺旋槳收到功率不變的情況下直接增加螺旋槳設(shè)計轉(zhuǎn)速,可通過提高主機選型轉(zhuǎn)速達成;②通過加大螺旋槳的輕轉(zhuǎn)裕度來提高船舶的加速性能是一種非常有效的方法[3],也是保證螺旋槳收到功率不變的情況下改變系柱曲線最簡單直接的方式。國際拖曳水池會議(International Towing Tank Cofference,ITTC)在其2008年頒發(fā)的指導(dǎo)文件中建議輕轉(zhuǎn)裕度取5%~7%[4];③螺旋槳收到功率不變的情況下,減小螺旋槳轉(zhuǎn)矩系數(shù)kq值;④螺旋槳收到功率不變的情況下,減小螺旋槳直徑。這些修改方案都涉及到螺旋槳的修改,對敞水效率有影響,設(shè)計槳時應(yīng)視具體情況而定。

3 主機越過轉(zhuǎn)速禁區(qū)的能力

3.1 轉(zhuǎn)速禁區(qū)上界的功率余量

如圖4所示,主機轉(zhuǎn)速禁區(qū)最好分布在系柱拖力曲線上穿主機轉(zhuǎn)矩限制線之前,這樣,當(dāng)主機轉(zhuǎn)速越過轉(zhuǎn)速禁區(qū)上界的時候,系柱拖力處于主機轉(zhuǎn)矩范圍之內(nèi),船舶在跨越主機轉(zhuǎn)速禁區(qū)后遭遇惡劣工況,主機不會因為輸出的功率小于螺旋槳系柱拖力而無法提供更多的轉(zhuǎn)矩,主機轉(zhuǎn)速禁區(qū)的設(shè)置要充分考慮這一點。所以轉(zhuǎn)速禁區(qū)的位置與螺旋槳系柱曲線、主機轉(zhuǎn)矩限制線之間要有合理的關(guān)系。通過計算螺旋槳拖力曲線、主機轉(zhuǎn)矩限制線,確保主機轉(zhuǎn)速禁區(qū)上界低于拖力曲線與主機轉(zhuǎn)矩限制線的交點,并且距離主機轉(zhuǎn)矩限制線有足夠的余量,該余量就是判斷主機能否快速越過轉(zhuǎn)速禁區(qū)的標(biāo)尺,是判斷船舶加速能力的關(guān)鍵,可作為量化指標(biāo)要求確保船舶的抗風(fēng)浪能力。當(dāng)轉(zhuǎn)速禁區(qū)功率裕度(轉(zhuǎn)矩余量)較高時,通過轉(zhuǎn)速禁區(qū)時間較短,但如果轉(zhuǎn)速禁區(qū)功率(轉(zhuǎn)矩余量)較低時,通過時間將快速增加[5]。

圖4 拖力、轉(zhuǎn)矩、主機功率線(轉(zhuǎn)矩換算為功率后繪制)

DNVGL船級社對轉(zhuǎn)速禁區(qū)上界的主機功率給出了量化建議:轉(zhuǎn)速禁區(qū)上界的功率余量應(yīng)足以確保轉(zhuǎn)速禁區(qū)在所有預(yù)期條件下的合理通過時間;10%的功率余量被認為是足夠的,但是,如果可以接受較低的功率裕度,則需要經(jīng)驗和/或擴展分析的支持;無需在船上測試驗證功率余量[6]。

3.2 提高主機轉(zhuǎn)矩限制提供轉(zhuǎn)矩余量

避免低速通過轉(zhuǎn)速禁區(qū)最基本的指導(dǎo)意見是避免轉(zhuǎn)速禁區(qū)擴大到主機轉(zhuǎn)速的60%以上,倘若高于該值,則需要格外注意禁區(qū)上界的主機扭矩限制值是否已經(jīng)低于螺旋槳系柱拖力值,若如此,該船在航行過程中是無法快速越過轉(zhuǎn)速禁區(qū)的,加速能力就很差。

如圖5所示,轉(zhuǎn)速禁區(qū)正好包含了拖力曲線與主機轉(zhuǎn)矩限制線的交點,意味著螺旋槳在重載工況下,主機不能有效越過轉(zhuǎn)速禁區(qū),所以此時就需要主機開啟主機動態(tài)限制器(dynamic limiter function,DLF)。在轉(zhuǎn)速禁區(qū)上界,系柱拖力與主機轉(zhuǎn)矩之間只有很小的轉(zhuǎn)矩余量甚至無余量時,可以采用DLF增大主機某一轉(zhuǎn)速區(qū)間的功率。

圖5 轉(zhuǎn)速禁區(qū)主機轉(zhuǎn)矩限制增大示意

DLF功能說明:目前MAN新主機的主機控制系統(tǒng)(ECS)都自帶一套快速通過轉(zhuǎn)速禁區(qū)的控制模塊,簡稱DLF,主機調(diào)試工程師可以通過調(diào)整DLF模塊的設(shè)定值,使轉(zhuǎn)速禁區(qū)內(nèi)的主機的轉(zhuǎn)距限制線局部提高10%~20%至機器本身的極限,從而提供局部轉(zhuǎn)速下更大的轉(zhuǎn)矩。

3.3 重載轉(zhuǎn)速與輕載轉(zhuǎn)速之間的差值

根據(jù)經(jīng)驗及主機廠家的推薦,相對于螺旋槳輕載曲線,重載曲線(系柱拖力曲線)處于距離輕載轉(zhuǎn)速15%~20%轉(zhuǎn)速的位置能保證主機快速越過轉(zhuǎn)速禁區(qū)。倘若沒有更準(zhǔn)確的數(shù)據(jù)可用時,可以參考該條建議。通常取距離輕載轉(zhuǎn)速的17.5%轉(zhuǎn)速值作為螺旋槳系柱拖力曲線轉(zhuǎn)速的預(yù)設(shè)值,見圖4。從另一個角度看,螺旋槳設(shè)計好后所得的拖力曲線轉(zhuǎn)速是否距離輕載轉(zhuǎn)速15%~20%區(qū)間也是檢驗船舶加速能力的一個參考。如果根據(jù)設(shè)計槳所得到的拖力曲線超過了該推薦區(qū)間,則該船的加速能力需要引起關(guān)注,需要計算評估高海況下的螺旋槳推力是否仍然在主機轉(zhuǎn)矩限制線范圍以內(nèi)。雖然如此,該檢驗數(shù)值僅僅作為參考,不能作為評判,不能作為評判螺旋槳性能的決定性指標(biāo)。

4 加速能力評估方法

4.1 第一層衡準(zhǔn)

4.1.1 計算方法

設(shè)計槳確定的情況下,將主機轉(zhuǎn)矩限制曲線與螺旋槳系柱拖力曲線在主機轉(zhuǎn)速禁區(qū)上界進行比較,得到轉(zhuǎn)矩余量大于10%。

需要已知的數(shù)據(jù):①螺旋槳kt、kq、n、D;②主機裝機功率MCR、額定轉(zhuǎn)速N、轉(zhuǎn)速禁區(qū)BSR。

由此可通過式(8)計算系柱拖力,對比轉(zhuǎn)速禁區(qū)上界對應(yīng)的主機轉(zhuǎn)矩限制值,計算轉(zhuǎn)矩余量。

4.1.2 預(yù)測方法

無設(shè)計槳情況下,默認由輕載曲線減小17.5%的轉(zhuǎn)速得到螺旋槳系柱拖力,由此預(yù)測螺旋槳系柱拖力,再計算轉(zhuǎn)速禁區(qū)上界轉(zhuǎn)矩余量,判斷該余量是否大于10%。

4.2 第二層衡準(zhǔn)

在第一層衡準(zhǔn)不能滿足的情況下,或者輕、重轉(zhuǎn)速差超出推薦范圍時,應(yīng)判斷船舶在實際可能遭遇的高海況下及常用轉(zhuǎn)速(一般為CSR轉(zhuǎn)速)下是否具有加速能力,即3級、4級或者更高海況下轉(zhuǎn)矩余量是否大于0。

4.3 計算須輸入的數(shù)據(jù)

1)螺旋槳kt、kq、n、D。

2)螺旋槳收到功率PL的輕載值。

3)主機裝機功率MCR、額定轉(zhuǎn)速N、轉(zhuǎn)速禁區(qū)BSR。

4)船舶在3級、4級或更高海況下的風(fēng)、浪增阻,得到螺旋槳帶海況的PL輕載功率值。

4.4 加速能力評估計算輸出的數(shù)據(jù)

1)主機在轉(zhuǎn)速禁區(qū)上界、下界的轉(zhuǎn)矩(功率)余量X1,見圖6。

圖6 計算輸出示意

2)系柱狀態(tài)下(螺旋槳進速系數(shù)J=0 時),螺旋槳旋轉(zhuǎn)產(chǎn)生的最大轉(zhuǎn)矩Qmax、最大功率Pmax,及對應(yīng)最大系柱轉(zhuǎn)速nmax,計算輕載線與系柱重載線在這一功率下的轉(zhuǎn)速差X2,見圖6。

3)計算在主機CSR轉(zhuǎn)速下,螺旋槳的輕載值、3級海況值、4級海況值或更高海況對應(yīng)的功率距離主機轉(zhuǎn)矩限制線的轉(zhuǎn)矩余量X3、X4、X5(見圖6),判斷在最高海況下X5是否大于0,再參考船舶以后正常航行可能面臨的海況決定保留余量值。

計算流程見圖7~9。

圖7 轉(zhuǎn)矩余量計算流程

圖8 重載轉(zhuǎn)速與輕載轉(zhuǎn)速差計算流程

圖9 高海況下船舶加速能力判斷流程

4.5 計算實例

以2020年上海船舶研究設(shè)計院為大連信風(fēng)海運公司設(shè)計的某集裝箱船為例。

1)已知數(shù)據(jù)見表1。

表1 實例計算輸入信息

2)加速能力評估。計算主機轉(zhuǎn)速禁區(qū)(BSR)轉(zhuǎn)矩余量計算,見表2。

表2 主機轉(zhuǎn)速禁區(qū)(BSR上界)轉(zhuǎn)矩余量計算

表2中計算系柱拖力由P(n)=2πkqρD5·n3

計算得到;主機轉(zhuǎn)矩限制由主機廠家提供計算數(shù)值,同時由主機廠家提供轉(zhuǎn)速禁區(qū)范圍。計算得到在主機轉(zhuǎn)速禁區(qū)上屆轉(zhuǎn)矩余量X1為16.5%,大于要求的10%,符合要求,認為該船在風(fēng)浪中的加速能力可被接受,主機也有能力快速越過轉(zhuǎn)速禁區(qū)。

另外,從圖10可見系柱拖力曲線上穿主機轉(zhuǎn)矩限制線的位置高于60%N,且遠離轉(zhuǎn)速禁區(qū),表明拖力曲線、轉(zhuǎn)速禁區(qū)與主機轉(zhuǎn)矩限制線相對位置合理。

圖10 實船計算輸出示意

3)螺旋槳轉(zhuǎn)速差評估(參考)。螺旋槳輕/重載轉(zhuǎn)速差(LRM)見表3。

表3 螺旋槳輕/重載轉(zhuǎn)速差(LRM)計算

通過式(7)計算,螺旋槳最大系柱轉(zhuǎn)速為67.1 r/min,對應(yīng)收到功率8 208.5 kW;船模試驗結(jié)果:該功率對應(yīng)的輕載轉(zhuǎn)速為83.1 r/min,則可計算得到轉(zhuǎn)速差X2為23.9%,超出了船級社及主機廠家推薦的15%~20%的范圍。加之此前計算得到轉(zhuǎn)速禁區(qū)上界轉(zhuǎn)矩余量(16.5%)大于10%,滿足加速能力的要求,所以轉(zhuǎn)速差X2的計算結(jié)果僅作為參考。

但倘若沒有X1計算作為支撐,而X2結(jié)果又不理想,就有必要對風(fēng)浪中的加速能力作進一步評估,即計算高海況下的轉(zhuǎn)矩余量計算。

4)4級海況下加速能力評估。該船將來營運航線最高只需考慮4級海況,所以本船就以4級海況作為高海況計算條件。要評估船舶在4級海況中的加速能力,首先要計算風(fēng)浪增阻,具體計算依據(jù)ITTC(國際拖曳水池會議)標(biāo)準(zhǔn)計算公式,不再詳述。

表5中PL_SS4為4級海況下螺旋槳的收到功率,包含了通過ITTC方法計算的風(fēng)浪增阻。通過表5計算可以看到,4級海況下主機在CSR轉(zhuǎn)速下,依然具有24.3%(X5)的轉(zhuǎn)矩余量,結(jié)果遠大于0,表明該船在CSR轉(zhuǎn)速下可承受更高的海況。

表5 4級海況下CSR轉(zhuǎn)速下轉(zhuǎn)矩余量計算

通過以上計算可以判斷,該集裝箱船在風(fēng)浪中的加速能力良好。雖然X2超出推薦范圍,當(dāng)僅作為參考,通過X5的計算可以證明本船在實際海況下加速性能良好。

5 結(jié)論

1)基于螺旋槳水動力特性將微觀槳葉受力分析與實船宏觀加速工況相結(jié)合,將船舶在極端風(fēng)浪中的加速工況簡化為螺旋槳系柱工況的處理方法合理可行。

2)可以通過螺旋槳在同一系柱轉(zhuǎn)速下收到功率距離主機轉(zhuǎn)矩限制線之間的距離來體現(xiàn)加速能力。距離越大,加速能力越強;距離越小,加速能力越弱。通過計算船舶系柱狀態(tài)下在轉(zhuǎn)速禁區(qū)上界的轉(zhuǎn)矩余量作為第一層衡準(zhǔn),對船舶在風(fēng)浪中的加速能力做出判斷。在不能滿足第一層衡準(zhǔn)的情況下,通過計算船舶高海況下在螺旋槳常用轉(zhuǎn)速下的轉(zhuǎn)矩余量進入第二層衡準(zhǔn),判斷實船在實際可預(yù)見的帶風(fēng)浪增阻的高海況下是否具有轉(zhuǎn)矩余量。

3)螺旋槳重載轉(zhuǎn)速與輕載轉(zhuǎn)速之間的差值介于15%~20%之間的推薦可以作為設(shè)計初期(無設(shè)計槳信息情況下)的參考值,再使用文中描述的預(yù)測方法用來計算轉(zhuǎn)速禁區(qū)上界轉(zhuǎn)矩余量,也可以作為是否需要計算高海況下加速能力的前提條件。

4)應(yīng)避免轉(zhuǎn)速禁區(qū)包含拖力曲線與主機轉(zhuǎn)矩限制線交點的情況,否則螺旋槳在重載工況下,主機不能有效越過轉(zhuǎn)速禁區(qū),此時就需要主機提高局部轉(zhuǎn)速區(qū)間的轉(zhuǎn)矩余量。

5)所設(shè)計的船舶加速能力評估流程,經(jīng)實船計算應(yīng)用,證明合理可行。

猜你喜歡
船舶能力
消防安全四個能力
計算流體力學(xué)在船舶操縱運動仿真中的應(yīng)用
基于改進譜分析法的船舶疲勞強度直接計算
幽默是一種能力
船舶!請加速
BOG壓縮機在小型LNG船舶上的應(yīng)用
大興學(xué)習(xí)之風(fēng) 提升履職能力
你的換位思考能力如何
努力拓展無人機飛行能力
無人機(2017年10期)2017-07-06 03:04:36
船舶壓載水管理系統(tǒng)
中國船檢(2017年3期)2017-05-18 11:33:09
主站蜘蛛池模板: 亚洲国产欧美中日韩成人综合视频| 99久久精品国产精品亚洲| 天天色天天综合| 中文字幕色站| 国产丝袜第一页| 在线观看亚洲成人| 91欧美在线| 91精品专区国产盗摄| 丝袜无码一区二区三区| 国产精品久久久久久影院| 在线免费看片a| 人妻一本久道久久综合久久鬼色| 亚洲天堂精品在线| 91国内在线视频| 国产精品999在线| 精品人妻无码中字系列| 亚洲欧洲日韩综合色天使| 欧美啪啪视频免码| 五月天福利视频| 国产性精品| 久视频免费精品6| 国产主播一区二区三区| 呦视频在线一区二区三区| 久久精品免费国产大片| 欧美国产菊爆免费观看| 日韩 欧美 国产 精品 综合| 久久久久中文字幕精品视频| 波多野结衣视频一区二区| 欧美综合激情| 亚洲综合极品香蕉久久网| 婷婷色婷婷| 国内视频精品| 国产91精品最新在线播放| 91丝袜美腿高跟国产极品老师| 成人午夜天| 亚洲一区精品视频在线| 9久久伊人精品综合| 污网站在线观看视频| 日本国产一区在线观看| 国产男女免费完整版视频| AV在线天堂进入| 国产成人无码AV在线播放动漫| 欧美国产日本高清不卡| 在线观看国产精品一区| 亚洲午夜久久久精品电影院| 看看一级毛片| 91人人妻人人做人人爽男同| 中文字幕精品一区二区三区视频| 精品国产网| 日韩黄色在线| 国产在线视频欧美亚综合| 成年人免费国产视频| 77777亚洲午夜久久多人| 超清无码一区二区三区| 亚洲一区网站| 久久国产精品电影| 日韩在线观看网站| 国产高潮流白浆视频| 国产成人a在线观看视频| 国产免费高清无需播放器| 欧美国产日韩一区二区三区精品影视| 亚洲欧美日韩高清综合678| 免费人成在线观看成人片| 毛片网站在线看| 欧美在线观看不卡| 亚洲第一网站男人都懂| 中文纯内无码H| 日韩免费毛片视频| 日本免费一级视频| 中文纯内无码H| 99激情网| 综合人妻久久一区二区精品| 日韩精品无码不卡无码| 亚洲中文字幕在线观看| 国产一级二级三级毛片| 欧美性爱精品一区二区三区| 国产三区二区| 国产视频自拍一区| 丁香五月婷婷激情基地| 影音先锋亚洲无码| 欧美日韩动态图| 精品国产成人国产在线|