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

    多槳船附體阻力的預(yù)報(bào)方法*

    2010-04-03 02:00:42
    船舶 2010年2期
    關(guān)鍵詞:附體實(shí)船船模

    馮 毅

    (七○八研究所 上海 200011)

    多槳船附體阻力的預(yù)報(bào)方法*

    馮 毅

    (七○八研究所 上海 200011)

    多槳船;多附體;附體阻力;預(yù)報(bào)方法

    多槳船是艦船的必要船種,而“多槳船附體阻力預(yù)報(bào)方法”是多槳船實(shí)船快速性預(yù)報(bào)的基礎(chǔ)。對(duì)現(xiàn)有附體阻力預(yù)報(bào)方法進(jìn)行歸納,介紹了經(jīng)驗(yàn)公式估算方法、船模試驗(yàn)估算方法和CFD數(shù)值計(jì)算方法的基本原理,討論現(xiàn)有方法對(duì)于多槳船附體阻力預(yù)報(bào)的可行性,具有實(shí)際意義。

    0 引 言

    一般來(lái)說(shuō),將船殼外承擔(dān)一定功能的附屬體稱(chēng)為附體,又將水線(xiàn)以下的附體引起的水阻力稱(chēng)之為附體阻力。按照附體阻力不同成因,可分為由流體重力引起的興波阻力和由流體粘性引起的粘性阻力。

    一般船舶包括基本的附體如:槳軸、軸支架、螺旋槳、舵;可能還包括改善航向穩(wěn)定性的裝置,如呆木;改善船體的橫搖運(yùn)動(dòng)的減搖裝置,如舭龍骨、減搖鰭;改善船體的縱搖運(yùn)動(dòng)的減搖裝置如壓浪板、首鰭;改善流場(chǎng)狀態(tài)及結(jié)構(gòu)等多重因素的軸包套、尾鰭等;以及一些節(jié)能裝置的附體,如轂帽鰭,整流隔板等。另外,還有專(zhuān)為特定用途而在船體表面開(kāi)孔或船體表面凸出體所引起的阻力,如側(cè)推孔、聲吶圍井、導(dǎo)流槽或?qū)Я髡帧⒙晠妊b置等。

    船舶附體對(duì)艦船的快速性能有不可忽視的影響,隨著船舶建造技術(shù)的發(fā)展,船舶附體的種類(lèi)也愈加豐富。所以,在設(shè)計(jì)附體時(shí)不但要充分發(fā)揮附體本身的功能,另外也不得不考慮盡可能降低其阻力。

    多槳船的附體數(shù)量較多且附體阻力占總阻力百分比較大。螺旋槳數(shù)量的增加會(huì)導(dǎo)致軸、軸包套、軸支架、舵等附體數(shù)量增加,這對(duì)其快速性能有很大影響。若將每個(gè)螺旋槳所對(duì)應(yīng)的軸、軸包套、軸支架、舵等附體看成一組推進(jìn)裝置,則各組推進(jìn)裝置之間的橫向間距會(huì)隨螺旋槳數(shù)目的增加而減小,附體與附體之間流場(chǎng)的相互干涉會(huì)更加明顯,船體與附體之間的相互影響亦會(huì)更加復(fù)雜。因此,尋求恰當(dāng)便捷的方法估算多槳船附體阻力十分重要。

    圖1 船舶附體示意圖

    多槳船(雙槳以上)在船舶總量中所占比例較小,以往國(guó)內(nèi)外關(guān)于多槳船附體阻力估算的詳細(xì)研究資料很少公開(kāi),這給多槳船附體阻力的估算帶來(lái)了困難。因此,需要借鑒單槳及雙槳船的估算方法來(lái)研究多槳船附體阻力。目前,確定船舶附體阻力的方法通常分為依經(jīng)驗(yàn)公式估算的方法、船模試驗(yàn)估算的方法和CFD數(shù)值計(jì)算的方法。本文分別闡述經(jīng)驗(yàn)公式法、船模試驗(yàn)估算法及CFD數(shù)值計(jì)算法的基本情況,并討論其對(duì)于多槳船附體阻力估算時(shí)的優(yōu)點(diǎn)及不足。

    20世紀(jì)70年代,日本人大隅三彥對(duì)于高速快艇的推進(jìn)性能作過(guò)一些研究,其中涉及到了多槳船附體阻力估算方法。他認(rèn)為高速艇由舵、尾軸支架、海水吸口等附加物造成的水阻力與裸船體阻力之比值隨速度增加而逐漸變小,每軸約5%~10%,若附體阻力按照平均每軸7%計(jì),附加阻力系數(shù)k=0.07×軸數(shù)。

    陳良權(quán)等[5]以模型試驗(yàn)的方法為主,分別用通常的傅汝德?lián)Q算方法、考慮尺度效應(yīng)的傅汝德?lián)Q算方法(β法)、三因次法(1+k法)、附體阻力系數(shù)法(△C法)、附體百分比計(jì)算方法,對(duì)雙槳船的附體阻力換算方法進(jìn)行比較,得出了關(guān)于雙槳船附體阻力估算的一些結(jié)論:誤差隨Fn的增大而相對(duì)減小,附體阻力占總阻力的比重隨Fn增加而減小,隨吃水增大而減少;1+k法與β法較接近,△C法與傅汝德法較接近,附體百分比法居中。Fn>0.4后傅汝德法較好,低速時(shí)1+k法較好。

    鑒于小船模的附體阻力試驗(yàn)及其計(jì)算結(jié)果的誤差相對(duì)較大,以及小船模附體存在較大的尺度效應(yīng)影響等因素,李定尊建議制作大船模進(jìn)行阻力試驗(yàn),還提出了對(duì)于船模試驗(yàn)的其他一些要求[6],如流線(xiàn)試驗(yàn)是設(shè)計(jì)、安裝附體的重要依據(jù);同一船體不同附體方案進(jìn)行阻力性能比較時(shí),以同船長(zhǎng)、同吃水為宜;船模試驗(yàn)時(shí)的排水量應(yīng)包括實(shí)船上為非水密而在船模上為水密的那一部分(如長(zhǎng)軸包套),以免船模阻力試驗(yàn)時(shí)因不能模擬實(shí)船非水密而造成試驗(yàn)時(shí)阻力偏低;做“分解”附體阻力試驗(yàn)時(shí),要注意各附體裝和拆的順序,防止附體設(shè)計(jì)、加工、安裝的不對(duì)稱(chēng)造成船模航態(tài)橫傾??紤]尺度效應(yīng)的傅汝德?lián)Q算方法(β法),陳紀(jì)平等[8]采用多條船模試驗(yàn)與實(shí)船試航結(jié)果比較得出β系數(shù)平均值可作為參考指標(biāo),他們還就1+k法中形狀因子的取法提出了新的觀點(diǎn),并提出了船?!獙?shí)船相關(guān)分析的方法。

    針對(duì)CFD數(shù)值計(jì)算方法,楊卓懿對(duì)于多附體船舶阻力進(jìn)行了研究。文獻(xiàn)[12]中分別計(jì)算了考慮主附體之間干擾情況下的阻力,鑒于計(jì)算機(jī)的性能,未直接進(jìn)行全附體模型的計(jì)算。而是將船體總阻力分解為主船體摩擦力(平板公式計(jì)算)、各簡(jiǎn)化附體阻力(結(jié)合已有的經(jīng)驗(yàn)公式計(jì)算)和附體干擾下的主體興波阻力三部分來(lái)計(jì)算。研究附體對(duì)主流場(chǎng)的影響時(shí),將附體考慮成一個(gè)與主體船同速運(yùn)動(dòng)的獨(dú)立部分,利用細(xì)長(zhǎng)船理論,分別求得附體的Kochin函數(shù),再與主體船的Kochin函數(shù)進(jìn)行線(xiàn)性迭加,得到一種在附體干擾下的主體船阻力的計(jì)算方法[12]。結(jié)果表明后一種考慮主附體之間干擾情況下計(jì)算出的總阻力更接近船模試驗(yàn)值,而未考慮主附體之間干擾的總阻力略低。

    將現(xiàn)有的附體阻力估算方法歸納為以下三類(lèi)。

    1 經(jīng)驗(yàn)公式法[2,3]

    經(jīng)驗(yàn)公式法是一種結(jié)合經(jīng)驗(yàn)數(shù)據(jù)和理論公式來(lái)確定附體阻力的方法[1]。用經(jīng)驗(yàn)公式法直接估算各附體阻力,總的阻力則將估算所得的各部分附體阻力迭加即可。

    (1)舭龍骨

    一般長(zhǎng)度約為船長(zhǎng)L的30%~50%,沿流線(xiàn)方向設(shè)置于船體舭部,其阻力大小可在計(jì)算主體摩擦阻力時(shí)在濕面積中加進(jìn)舭龍骨的濕面積一并計(jì)算,或取為裸體阻力的1%~3%,或取舭龍骨所增加的總阻力,rbf是它的平板摩擦阻力)。

    (2)舵

    對(duì)于不同尾型、舵型及舵的安裝位置,所產(chǎn)生的附加阻力不同。流線(xiàn)型舵的阻力可取本身摩擦阻力的1.5倍。計(jì)算舵的摩擦阻力亦可按相當(dāng)平板計(jì)算,其雷諾數(shù)特征長(zhǎng)度用本身弦長(zhǎng),來(lái)流速度應(yīng)對(duì)船體伴流及螺旋槳尾流加以修正。單槳船舵的阻力一般占裸船體阻力的1%~2%;雙槳船的雙邊舵的阻力一般占裸船體阻力的3~5%。但由于舵處于螺旋槳的尾流中,對(duì)尾流有整流作用,可以提高推進(jìn)效率;有的舵(如雙槳船的中舵、駁船的舵等)相當(dāng)于增加有效船長(zhǎng),對(duì)減小船體粘壓阻力有利。因此,若舵的有利影響可抵消本身的阻力,則舵的阻力可忽略不計(jì)。

    (3)塢座龍骨

    其阻力可取作自身摩擦阻力的4倍。

    (4)軸包套

    通常在多槳船上設(shè)置,沿流線(xiàn)方向布置,使之不產(chǎn)生渦流。設(shè)計(jì)好的軸包套產(chǎn)生的粘壓阻力極小。軸包套的阻力一般為裸體阻力的2%~12%。

    表1 軸包套占裸船體百分比

    (5)軸支架和軸

    軸支架主要受到的是粘壓阻力,因此,剖面應(yīng)采用流線(xiàn)型,弦長(zhǎng)方向應(yīng)沿流線(xiàn)布置。軸支架的總阻力可按下列經(jīng)驗(yàn)公式計(jì)算:

    式中 k——系數(shù),當(dāng)軸支架剖面最大厚度為1/3弦

    長(zhǎng)時(shí),k=6.416;當(dāng)軸支架剖面最大厚度為

    1/5弦長(zhǎng)時(shí),k=2.139;

    S——軸支架的濕表面積(m2);

    v——軸支架與水流的相對(duì)速度(m/s)。b.軸支架經(jīng)驗(yàn)公式2

    式中 b——螺旋槳槳軸中心線(xiàn)到中線(xiàn)面的距離(m);

    v——船速(m/s)。

    此式僅適用于支臂總長(zhǎng)為3b、平均厚度為0.054b的情況,其他尺度情況可按正比換算得到。

    式中 l——尾軸在船外的長(zhǎng)度(m);

    d——尾軸直徑(m);

    v——船速(m/s);

    θ——水流方向與尾軸向的夾角(°)。一般尾軸前端θ值對(duì)于豐滿(mǎn)船約為20°,對(duì)瘦削船約為10°;其后端約為2°。應(yīng)分段確定夾角的大小,但實(shí)際情況該夾角是很難確定的。

    (6)螺旋槳軸制動(dòng)產(chǎn)生的阻力可按下式估算:

    式中 P/D——螺旋槳的螺距比;

    AE——螺旋槳的伸張面積(m2);

    vA——螺旋槳的進(jìn)速(m/s)。

    (7)25thITTC推薦用局部雷諾數(shù)及形狀因子的方法,對(duì)每個(gè)附體阻力的阻力系數(shù)進(jìn)行單獨(dú)計(jì)算:

    式中 n——附體的件數(shù);

    wi——第i個(gè)附體對(duì)應(yīng)的伴流分?jǐn)?shù);ki——第i個(gè)附體所對(duì)應(yīng)的形狀因子;

    CFsi——第i個(gè)附體所對(duì)應(yīng)的摩擦阻力系數(shù);

    Si——第i個(gè)附體所對(duì)應(yīng)的濕表面積;

    SS——船體的濕表面積。

    wi是模型試驗(yàn)得到,ki可以按(Hoerner(1965)或Kirkman和Klostsli(1980))的公式計(jì)算。

    (8)根據(jù)統(tǒng)計(jì)數(shù)據(jù),將幾種附體結(jié)合考慮的粗略估算公式,如前面提到的大隅三彥使用的對(duì)于多槳高速艇的附體阻力估算方法。

    假設(shè)附體有效功率系數(shù)為

    則帶有附體的總有效功率為

    式中 PeA——安裝全部附體后的有效功率;

    Peb——裸船體的有效功率。

    (9)利用大量的母型船的附體阻力統(tǒng)計(jì)資料估算目標(biāo)船附體阻力的方法。

    一般地,對(duì)于數(shù)量較少且設(shè)計(jì)合理的附體來(lái)說(shuō),附體阻力可能低于主船體總阻力的10%,而數(shù)量較多或設(shè)計(jì)不合理的附體來(lái)說(shuō),附體阻力可能超過(guò)主船體總阻力的25%。

    表2 附體阻力占裸船體阻力百分比

    對(duì)于一些附體數(shù)量較少、常規(guī)結(jié)構(gòu)附體的船來(lái)說(shuō),依靠經(jīng)驗(yàn)公式方法能夠較快速地將其附體阻力估算出來(lái),以預(yù)報(bào)實(shí)船附體阻力。該方法比較適合于初步設(shè)計(jì)時(shí)粗略估算附體阻力,因?yàn)椴⑽纯紤]到附體與附體、附體與主船體之間的相互干擾,尤其對(duì)于多槳船來(lái)說(shuō),其附體數(shù)量較多、結(jié)構(gòu)相對(duì)復(fù)雜的特點(diǎn)使得這種依經(jīng)驗(yàn)公式估算得到的附體阻力累加的方法不甚合理。在采用相當(dāng)平板公式計(jì)算附體表面的摩擦阻力時(shí),對(duì)于雷諾數(shù)的特征長(zhǎng)度的認(rèn)定方法仍有存疑;又如軸的計(jì)算公式中的θ和v實(shí)際上都是難以準(zhǔn)確測(cè)定的。另外,現(xiàn)有的經(jīng)驗(yàn)公式估算公式覆蓋面有限,僅能對(duì)一些常規(guī)的附體情況給予粗略計(jì)算。隨著造船技術(shù)的發(fā)展,當(dāng)遇到某些新形式或特殊外形附體時(shí),經(jīng)驗(yàn)公式估算方法可能無(wú)法進(jìn)行。此時(shí),需依靠船模試驗(yàn)估算方法或CFD數(shù)值計(jì)算的方法來(lái)解決。

    經(jīng)驗(yàn)公式法是依照既往經(jīng)驗(yàn)、數(shù)據(jù)資料總結(jié)而來(lái)。對(duì)于常規(guī)形式附體阻力的估算還可以借助大量試驗(yàn)資料建立數(shù)據(jù)庫(kù)的方法。對(duì)船型及附體類(lèi)型的試驗(yàn)資料進(jìn)行分類(lèi),以便能夠快速得到目標(biāo)船附體阻力的大致量。對(duì)于多附體船型附體之間的相互影響可嘗試根據(jù)數(shù)據(jù)庫(kù)資料歸納不同類(lèi)型附體間的影響系數(shù)公式,以便得到相對(duì)合理的附體阻力。解決這些新的問(wèn)題需要相當(dāng)長(zhǎng)時(shí)間的資料積累,相對(duì)于船模試驗(yàn)法和CFD數(shù)值計(jì)算法有一定局限性。

    2 船模試驗(yàn)估算法

    船模試驗(yàn)估算方法一般是借助船模試驗(yàn)得到的裸船體的阻力試驗(yàn)和帶有附體的船模試驗(yàn)來(lái)估算出實(shí)船附體阻力的方法。

    2.1 通常的傅汝德?lián)Q算方法

    帶附體的船模阻力試驗(yàn)結(jié)果應(yīng)用傅汝德阻力換算方法直接計(jì)算得到帶全附體實(shí)船阻力

    用該方法換算沒(méi)有考慮到附體模型尺寸較小所引起的尺度效應(yīng)影響,故計(jì)算的實(shí)船附體阻力結(jié)果一般偏高。

    2.2 考慮尺度效應(yīng)的傅汝德?lián)Q算方法(β法)

    通常稱(chēng)英國(guó)法或附體尺度效應(yīng)修正β方法(簡(jiǎn)稱(chēng)β法)。β法中認(rèn)為帶附體的模型阻力與不帶附體的模型阻力之差ΔRm可以用縮尺比λ的三次方的關(guān)系換算到實(shí)船的附體阻力ΔRs。若再考慮水密度的影響則有

    式中:β附體尺度效應(yīng)修正系數(shù),0.6~1.0。具體取決于各個(gè)拖曳水池的經(jīng)驗(yàn)數(shù)據(jù)[8]。中對(duì)于與實(shí)船試航結(jié)果最符合的內(nèi)插β值進(jìn)行了統(tǒng)計(jì)(βopt=0.707~0.960),平均值為0.844,一個(gè)典型值0.75被推薦用于常規(guī)雙槳船型文獻(xiàn)。這個(gè)方法經(jīng)證實(shí)在許多情況下很起作用,但是尺度效應(yīng)修正因子原則上應(yīng)取決于模型和實(shí)船的雷諾數(shù);進(jìn)行帶附體的船模試驗(yàn)時(shí),由于模型的雷諾數(shù)相對(duì)較低,在附體上容易產(chǎn)生不受控制的難預(yù)測(cè)的層流和分離;某一附體的阻力及其尺度效應(yīng)受到一系列復(fù)雜因素的影響,不可能用一個(gè)簡(jiǎn)單的β便將它們準(zhǔn)確表達(dá)[4]。

    2.3 三因次法(1+k法)

    在帶附體的船模阻力試驗(yàn)換算中采用1+k方法。

    形狀因子1+k由低速試驗(yàn)法、普羅哈斯卡方法、15thITTC推薦法或幾何相似船模組試驗(yàn)法來(lái)確定,關(guān)于ΔCF求解方法可見(jiàn)文獻(xiàn)[7]。

    因?yàn)?+k方法將受Fn影響較小的粘壓阻力從剩余阻力中分離出來(lái),所以相對(duì)于二因次法更加合理。對(duì)帶有軸支架,軸包套的雙槳船用幾何相似船模組試驗(yàn)法確定形狀因子1+k,可使尺度效應(yīng)的影響大大降低,但這樣作的工作量過(guò)于繁重[9,10]。而應(yīng)用低速試驗(yàn)法、普羅哈斯卡方法、15thITTC推薦法較難正確求取形狀因子1+k,會(huì)導(dǎo)致?lián)Q算結(jié)果不穩(wěn)定。

    2.4 附體阻力系數(shù)法(ΔC法)

    由于附體一般都位于自由面以下且離開(kāi)自由面有一定距離,故可假定附體不興波,認(rèn)為附體阻力與傅汝德數(shù)無(wú)關(guān);當(dāng)雷諾數(shù)達(dá)到一臨界值后,流動(dòng)發(fā)生轉(zhuǎn)淚,流動(dòng)由層流流動(dòng)轉(zhuǎn)為湍流流動(dòng),此時(shí)附體阻力系數(shù)ΔC為一常數(shù),而與雷諾數(shù)無(wú)關(guān)。實(shí)船的雷諾數(shù)較大,認(rèn)為實(shí)船附體的雷諾數(shù)已經(jīng)足夠大,達(dá)到臨界值,實(shí)船附體阻力系數(shù)為常數(shù),那么,模型試驗(yàn)得到的ΔC就可以作為實(shí)船附體阻力系數(shù)。

    式中 RmA、Rm、ρm需修正到相同水溫。

    在船模附體阻力試驗(yàn),附體實(shí)際上還是存在一定興波影響,所以ΔC并不是一個(gè)常數(shù)[5,6]。

    2.5 附體百分比計(jì)算法

    附體百分比法計(jì)算相對(duì)比較簡(jiǎn)便,在船舶的初步設(shè)計(jì)階段通常采用這一方法。該方法是由模型試驗(yàn)、同類(lèi)型船試驗(yàn)結(jié)果、統(tǒng)計(jì)資料或經(jīng)驗(yàn)公式對(duì)目標(biāo)船附體阻力進(jìn)行預(yù)報(bào)的方法[8]。模型試驗(yàn)中附體阻力所占裸船體阻力百分比即為實(shí)船附體阻力所占裸船體阻力百分比:

    一般來(lái)說(shuō)模型試驗(yàn)估算法方法相對(duì)成熟、應(yīng)用廣泛,適用于預(yù)估大多數(shù)形式附體的附體阻力。各單位通常有其常用的換算方法。對(duì)于單槳船型,模型試驗(yàn)估算法已經(jīng)被廣泛采用,并能夠用得到比較理想的結(jié)果。鑒于多槳船附體多、附體阻力占總阻力的比重大,各種模型試驗(yàn)估算法是否還能夠得到合理的附體阻力值得探討。

    此外,模型試驗(yàn)換算方法會(huì)因各水池的試驗(yàn)條件不同,估算結(jié)果也有所差異。在試驗(yàn)條件允許的情況下,應(yīng)盡量制作較大的模型進(jìn)行試驗(yàn),以減小其尺度效應(yīng)帶來(lái)的影響。再者,船模加工的精度、試驗(yàn)儀器性能等因素也影響著附體阻力預(yù)報(bào)結(jié)果的精度。

    3 CFD數(shù)值計(jì)算法

    CFD數(shù)值計(jì)算法一般是借助流體力學(xué)專(zhuān)業(yè)軟件通過(guò)計(jì)算機(jī)數(shù)值計(jì)算對(duì)船舶附體阻力進(jìn)行預(yù)報(bào)的方法。這種方法最大的好處之一就是考慮到了附體與主船體間,附體與附體間及粘性阻力與興波阻力間的耦合影響,并描繪出流場(chǎng)內(nèi)的流動(dòng)細(xì)節(jié)和壓力分布情況。

    隨著計(jì)算機(jī)硬件性能的提升,目前通常的做法是將帶附體的船型直接劃分網(wǎng)格進(jìn)行計(jì)算。原則上來(lái)說(shuō)在實(shí)尺度條件下,應(yīng)用計(jì)算軟件分別對(duì)裸船體和全附體船分別進(jìn)行計(jì)算,其差值即為實(shí)船的附體阻力,并且避免了附體的尺度效應(yīng)問(wèn)題。然而,實(shí)尺度船型的網(wǎng)格數(shù)過(guò)于龐大,又鑒于軟件及硬件的計(jì)算能力的限制,大多數(shù)情況下對(duì)實(shí)尺度船型進(jìn)行計(jì)算并不可行。因此,在實(shí)際計(jì)算中常以一定的縮尺比的幾何相似模型為對(duì)象進(jìn)行計(jì)算,以便提高計(jì)算的效率。將計(jì)算得到的裸船體阻力及全附體總阻力再結(jié)合模型試驗(yàn)估算法中的各種換算方法進(jìn)行換算,最終得到實(shí)船的附體阻力。其實(shí),這種CFD數(shù)值計(jì)算方法中所需的建模及計(jì)算過(guò)程就相當(dāng)于模型試驗(yàn)估算法中制作模型及水池試驗(yàn)的過(guò)程,而由于存在縮尺比,尺度效應(yīng)是不可避免的。但CFD數(shù)值計(jì)算法建模速度快、成本低,是模型試驗(yàn)法所不及的優(yōu)點(diǎn)。

    與前兩法相比,CFD數(shù)值計(jì)算法是一種新興的方法,尚處于發(fā)展完善階段,其計(jì)算結(jié)果的可靠性還有待進(jìn)一步論證。隨著軟件及硬件的計(jì)算能力提升,實(shí)尺度的計(jì)算及尺度效應(yīng)估算將成為可能,當(dāng)積累相當(dāng)豐富經(jīng)驗(yàn)及試驗(yàn)資料后,CFD數(shù)值計(jì)算法將會(huì)成為附體阻力預(yù)估一種十分有效的方法。

    4 結(jié) 語(yǔ)

    三種計(jì)算方法的特點(diǎn)歸納如下:

    (1)經(jīng)驗(yàn)公式法可以較快地預(yù)報(bào)出多槳船附體阻力,比較適合于初步設(shè)計(jì)時(shí)使用。但實(shí)際計(jì)算時(shí),部分參數(shù)選取有困難,對(duì)新形式附體阻力預(yù)報(bào)有局限性,對(duì)于多附體的多槳船難以考慮到附體與主體之間的相互影響,因此,預(yù)報(bào)結(jié)果可能誤差偏大;

    (2)船模試驗(yàn)估算法方法相對(duì)成熟、適用情況廣泛。雖然換算方法眾多,但是多數(shù)換算方法還都難以避免尺度效應(yīng)的影響。對(duì)于多槳船附體阻力換算沒(méi)有一致的附體阻力換算方法,須對(duì)各種換算方法進(jìn)行進(jìn)一步的比較評(píng)估;

    (3)CFD數(shù)值計(jì)算法模型精準(zhǔn),且考慮到了附體與主船體間、附體與附體間及粘性阻力與興波阻力間的耦合影響及流動(dòng)細(xì)節(jié)情況,但其相對(duì)于經(jīng)驗(yàn)公式法和模型試驗(yàn)法還不夠成熟,提高其計(jì)算結(jié)果的可靠性有待進(jìn)一步研究。

    由此可見(jiàn),三種估算方法各有其優(yōu)點(diǎn)及不足,應(yīng)結(jié)合實(shí)際需求選擇適合的方法。初步設(shè)計(jì)采用經(jīng)驗(yàn)公式法粗略估算,可以節(jié)約大量時(shí)間;在附體方案選型時(shí)可以用CFD數(shù)值計(jì)算法進(jìn)行比較選優(yōu),再將選定方案進(jìn)行模型試驗(yàn)和附體阻力預(yù)報(bào);在實(shí)船試驗(yàn)后,將得到的試驗(yàn)數(shù)據(jù)與模型試驗(yàn)數(shù)據(jù)進(jìn)行相關(guān)分析,對(duì)模型試驗(yàn)換算方法進(jìn)行改進(jìn)。相信隨著CFD數(shù)值計(jì)算技術(shù)的發(fā)展,尺度效應(yīng)問(wèn)題將會(huì)得到更好的解決。

    [1] 邵世明,趙連恩,朱念昌.船舶阻力[M].哈爾濱:哈爾濱工程大學(xué)出版社,2003.79.

    [2] 李世謨,船舶阻力.[M].北京:人民交通出版社,1989.110~112.

    [3] 張繼猛.船模試驗(yàn)中附體阻力的換算方法[J].艦船科學(xué)技術(shù),1989,4:32~36.

    [4] 25thITTC自航試驗(yàn)規(guī)程[S].

    [5] 陳良權(quán),'92船舶阻力與性能學(xué)術(shù)討論會(huì)論文集[C].成都:中國(guó)造船,1992:151~160.

    [6] 李定尊.'95船舶阻力與性能學(xué)術(shù)討論會(huì)論文集[C].威海:中國(guó)造船,1995:123~131.

    [7] 沈玉林,沈奇心,趙漢魂等.三因次換算方法中相關(guān)補(bǔ)貼系數(shù)ΔCF計(jì)算方法的探討[J].船舶,2003,4:8~11.

    [8] 陳紀(jì)平,趙漢魂,沈玉林.船舶附體阻力換算方法研究[J].船舶,2009,2:13~17.

    [9] 姜次平.雙槳船Lucy Ashton軸系附體雙槳阻力的尺度效應(yīng).

    [10] 姜次平,陳良權(quán).雙螺旋槳船的美人架阻力的尺度效應(yīng)[J].上海交通大學(xué)學(xué)報(bào),1983,3:75~85.

    [11] 大隅三彥,木原和之.小型高速艇推進(jìn)性能的估算方法.(日)西部造船會(huì)會(huì)報(bào),(75):36~51.

    [12] 楊卓懿.船舶多附體興波阻力計(jì)算的研究[D].哈爾濱工程大學(xué),2007.

    Forecast on appendage resistance for multi-screw ship

    Feng Yi

    multi-screw ship;multi-appendage;appendage resistance;forecast methods

    The forecast on appendage resistance for multi-screw ship,an indispensable type of warship,is the base for rapidity forecast.This paper summarizes the current appendage resistance forecast methods,and introduces the principles of estimation methods by experimental for mu la and model test,and CFD numerical algorithm.It also discusses the feasibility of current method for multi-screw appendage resistance forecast.

    U 661.31+1

    A

    1001-9855(2010)02-0015-06

    2009-12-20

    馮 毅(1982.04-),男,漢族,黑龍江人,助理工程師,主要從事船舶水動(dòng)力性能研究。

    猜你喜歡
    附體實(shí)船船模
    船舶阻力試驗(yàn)實(shí)船換算的蒙特卡羅仿真
    基于模糊PID的船模航向控制研究
    基于多種組合算法的船附體結(jié)構(gòu)設(shè)計(jì)優(yōu)化
    開(kāi)運(yùn)年會(huì)
    女報(bào)(2020年2期)2020-06-12 11:37:49
    這屆雪人跑偏啦
    長(zhǎng)江1 000車(chē)滾裝船EEDI驗(yàn)證分析
    船海工程(2019年6期)2019-12-25 06:12:50
    5 668 TEU集裝箱船球艏改型與浮態(tài)節(jié)能的實(shí)船驗(yàn)證
    陸樹(shù)銘 關(guān)公“附體”之后
    船模靜水橫搖試驗(yàn)的不確定度分析
    西洋船模王——童鑑良
    航海(2014年6期)2014-12-12 10:36:03
    97人妻天天添夜夜摸| 国产成人免费观看mmmm| 人人妻人人澡人人爽人人夜夜| 考比视频在线观看| 免费女性裸体啪啪无遮挡网站| 婷婷色综合大香蕉| 青春草亚洲视频在线观看| 日本欧美国产在线视频| tube8黄色片| 国产在线免费精品| 国产精品成人在线| 国产视频首页在线观看| 侵犯人妻中文字幕一二三四区| 亚洲国产毛片av蜜桃av| 国产免费福利视频在线观看| 国产男人的电影天堂91| 欧美日韩av久久| 自线自在国产av| 99re6热这里在线精品视频| 高清黄色对白视频在线免费看| 日本vs欧美在线观看视频| 久久国产亚洲av麻豆专区| 丝袜在线中文字幕| 国产一区二区激情短视频 | 国产精品久久久久久av不卡| 天天躁夜夜躁狠狠躁躁| 深夜精品福利| 日本色播在线视频| 高清不卡的av网站| 久热久热在线精品观看| 国产一区二区激情短视频 | 女人高潮潮喷娇喘18禁视频| 免费高清在线观看日韩| 国产欧美日韩综合在线一区二区| 久久 成人 亚洲| 超碰97精品在线观看| 亚洲,欧美精品.| 男女边摸边吃奶| 一边亲一边摸免费视频| 少妇 在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 满18在线观看网站| 欧美中文综合在线视频| 亚洲精品成人av观看孕妇| 在线观看免费视频网站a站| 国产成人精品久久久久久| 最近最新中文字幕免费大全7| 亚洲成人手机| 少妇精品久久久久久久| 欧美日韩一区二区视频在线观看视频在线| 777米奇影视久久| 一本久久精品| 男人添女人高潮全过程视频| 亚洲精品视频女| 在线观看三级黄色| 一级黄片播放器| 国产av一区二区精品久久| 亚洲精品aⅴ在线观看| 国产男女内射视频| 下体分泌物呈黄色| 亚洲欧美成人精品一区二区| 国产精品无大码| 黄色怎么调成土黄色| 在线观看www视频免费| 欧美老熟妇乱子伦牲交| 波多野结衣一区麻豆| 天堂8中文在线网| 啦啦啦啦在线视频资源| 91精品三级在线观看| 91国产中文字幕| 在线观看国产h片| 日韩一区二区三区影片| 欧美日本中文国产一区发布| 精品一区二区三区四区五区乱码 | 亚洲av男天堂| 亚洲精品国产av成人精品| 新久久久久国产一级毛片| 久久久久久久久久人人人人人人| av网站在线播放免费| 亚洲精品日本国产第一区| 日韩视频在线欧美| 亚洲av欧美aⅴ国产| 高清av免费在线| 国产精品99久久99久久久不卡 | 男女啪啪激烈高潮av片| av天堂久久9| 亚洲三区欧美一区| 国产精品.久久久| 亚洲欧洲国产日韩| 最近手机中文字幕大全| 亚洲国产毛片av蜜桃av| 麻豆av在线久日| 国产视频首页在线观看| 国产爽快片一区二区三区| 成人免费观看视频高清| 波野结衣二区三区在线| 久久这里只有精品19| 午夜日韩欧美国产| 十八禁高潮呻吟视频| 精品国产露脸久久av麻豆| 成年女人毛片免费观看观看9 | 国产成人精品婷婷| 久久久精品区二区三区| 大香蕉久久网| 亚洲伊人久久精品综合| 欧美 亚洲 国产 日韩一| 免费少妇av软件| 成人毛片60女人毛片免费| www.熟女人妻精品国产| 久久久久久人妻| 国产成人精品在线电影| 一级片'在线观看视频| 欧美人与性动交α欧美软件| 亚洲欧美成人精品一区二区| 男人操女人黄网站| 9191精品国产免费久久| 秋霞伦理黄片| 国产黄色免费在线视频| 在线观看免费高清a一片| 新久久久久国产一级毛片| 满18在线观看网站| 久久久久人妻精品一区果冻| 国产爽快片一区二区三区| 大陆偷拍与自拍| 99热网站在线观看| 人人妻人人添人人爽欧美一区卜| 青青草视频在线视频观看| 我的亚洲天堂| 久久韩国三级中文字幕| 9热在线视频观看99| 伊人久久大香线蕉亚洲五| 国产又色又爽无遮挡免| 香蕉丝袜av| 99香蕉大伊视频| av一本久久久久| 99久久精品国产国产毛片| 午夜精品国产一区二区电影| 国产精品不卡视频一区二区| 狠狠婷婷综合久久久久久88av| 26uuu在线亚洲综合色| 美女xxoo啪啪120秒动态图| 亚洲国产精品成人久久小说| 久热这里只有精品99| 一级,二级,三级黄色视频| 黑人猛操日本美女一级片| 国产精品一二三区在线看| 亚洲精品美女久久av网站| 国产成人免费观看mmmm| 久久精品久久久久久久性| 国产一区二区三区综合在线观看| 一本色道久久久久久精品综合| 午夜激情av网站| 欧美人与善性xxx| 极品人妻少妇av视频| www日本在线高清视频| 久热这里只有精品99| 日本色播在线视频| 国产在视频线精品| 赤兔流量卡办理| 国产一区二区三区av在线| 久久午夜福利片| 欧美97在线视频| 国产黄色免费在线视频| av电影中文网址| 午夜免费观看性视频| 91国产中文字幕| 婷婷色综合www| 亚洲av免费高清在线观看| 丝袜美腿诱惑在线| 亚洲国产av影院在线观看| 精品国产国语对白av| 老熟女久久久| 美女高潮到喷水免费观看| 国产精品一国产av| 制服人妻中文乱码| 青春草国产在线视频| 亚洲精品在线美女| 九草在线视频观看| 中文欧美无线码| 成人影院久久| 搡老乐熟女国产| www日本在线高清视频| 丁香六月天网| 成人免费观看视频高清| 国产极品粉嫩免费观看在线| 日本av手机在线免费观看| 国产日韩欧美视频二区| 国产午夜精品一二区理论片| 亚洲精品aⅴ在线观看| 中文精品一卡2卡3卡4更新| 欧美 亚洲 国产 日韩一| 丝袜喷水一区| 久久久久久久久久人人人人人人| 国产精品秋霞免费鲁丝片| 美女福利国产在线| 巨乳人妻的诱惑在线观看| 香蕉国产在线看| 我要看黄色一级片免费的| 人妻一区二区av| 美女国产视频在线观看| 女人被躁到高潮嗷嗷叫费观| 麻豆av在线久日| 91精品伊人久久大香线蕉| 青春草视频在线免费观看| 波多野结衣一区麻豆| 热99国产精品久久久久久7| 1024香蕉在线观看| 少妇的丰满在线观看| 韩国av在线不卡| 国产成人精品在线电影| 日韩成人av中文字幕在线观看| 涩涩av久久男人的天堂| 亚洲美女黄色视频免费看| 精品久久久精品久久久| 黄片小视频在线播放| 少妇人妻久久综合中文| 香蕉精品网在线| 国产成人精品久久久久久| 日韩成人av中文字幕在线观看| 在线观看一区二区三区激情| 视频区图区小说| 男人爽女人下面视频在线观看| 亚洲在久久综合| 777米奇影视久久| 久久久久久人人人人人| av在线观看视频网站免费| 午夜福利影视在线免费观看| 成人免费观看视频高清| 午夜免费鲁丝| 18禁国产床啪视频网站| 只有这里有精品99| 一区二区三区乱码不卡18| 久久亚洲国产成人精品v| 色婷婷av一区二区三区视频| 黄色 视频免费看| 久久精品国产亚洲av高清一级| 久久久久久久久免费视频了| h视频一区二区三区| 欧美bdsm另类| 另类精品久久| 亚洲激情五月婷婷啪啪| 女性被躁到高潮视频| 亚洲,欧美精品.| 久久精品久久久久久噜噜老黄| 亚洲精品国产av成人精品| 亚洲视频免费观看视频| 可以免费在线观看a视频的电影网站 | av.在线天堂| 国产视频首页在线观看| 国产深夜福利视频在线观看| 成人毛片a级毛片在线播放| 欧美日韩视频高清一区二区三区二| 高清在线视频一区二区三区| 丝袜人妻中文字幕| 在现免费观看毛片| 久久久久久久久久久免费av| av天堂久久9| tube8黄色片| 91aial.com中文字幕在线观看| 欧美日韩国产mv在线观看视频| 两个人免费观看高清视频| 亚洲欧洲国产日韩| 精品国产超薄肉色丝袜足j| 免费大片黄手机在线观看| 人成视频在线观看免费观看| 久久人人97超碰香蕉20202| 亚洲伊人色综图| 久久久欧美国产精品| 色网站视频免费| 亚洲国产av新网站| 久久国产亚洲av麻豆专区| 亚洲国产欧美日韩在线播放| 久久av网站| 2022亚洲国产成人精品| 中国国产av一级| 久久精品久久久久久久性| 在线精品无人区一区二区三| 亚洲一区二区三区欧美精品| 国产视频首页在线观看| 热re99久久精品国产66热6| 欧美日本中文国产一区发布| 一区在线观看完整版| 久久97久久精品| 成人二区视频| 大香蕉久久网| 日日啪夜夜爽| 久久狼人影院| 中文欧美无线码| 黑人欧美特级aaaaaa片| 黄色 视频免费看| 午夜福利在线观看免费完整高清在| 精品久久蜜臀av无| 久久久精品区二区三区| 精品国产露脸久久av麻豆| 国产男人的电影天堂91| 亚洲精品国产av蜜桃| 成人亚洲精品一区在线观看| 天天躁夜夜躁狠狠久久av| 久久午夜综合久久蜜桃| 精品99又大又爽又粗少妇毛片| 久久久久国产网址| 十分钟在线观看高清视频www| a级毛片黄视频| 1024香蕉在线观看| 七月丁香在线播放| 一二三四中文在线观看免费高清| 亚洲国产欧美日韩在线播放| 高清黄色对白视频在线免费看| 丰满乱子伦码专区| 亚洲国产欧美日韩在线播放| 久久精品国产亚洲av高清一级| 9191精品国产免费久久| 韩国精品一区二区三区| 男的添女的下面高潮视频| 色视频在线一区二区三区| 国产日韩欧美亚洲二区| 美女大奶头黄色视频| 天美传媒精品一区二区| 老熟女久久久| 精品国产露脸久久av麻豆| 夫妻性生交免费视频一级片| 亚洲经典国产精华液单| 自拍欧美九色日韩亚洲蝌蚪91| 少妇 在线观看| 黄色配什么色好看| 18在线观看网站| 国产成人免费无遮挡视频| 日韩欧美一区视频在线观看| 女人精品久久久久毛片| 婷婷成人精品国产| 香蕉丝袜av| 两个人看的免费小视频| 免费少妇av软件| 久久这里只有精品19| 超碰97精品在线观看| 亚洲av欧美aⅴ国产| 欧美日韩国产mv在线观看视频| 亚洲,欧美,日韩| 26uuu在线亚洲综合色| 春色校园在线视频观看| 亚洲一区中文字幕在线| 看十八女毛片水多多多| 亚洲精品国产一区二区精华液| 婷婷色av中文字幕| a级毛片在线看网站| 国产精品久久久久成人av| 久久久a久久爽久久v久久| av免费观看日本| 久久国内精品自在自线图片| 午夜福利,免费看| 成人国产麻豆网| 国产黄色视频一区二区在线观看| av有码第一页| 久久久国产一区二区| 国产成人aa在线观看| 成年美女黄网站色视频大全免费| 国产一区二区 视频在线| 久久久久久久亚洲中文字幕| 免费观看a级毛片全部| 亚洲av免费高清在线观看| 国产成人精品在线电影| 十分钟在线观看高清视频www| 国产精品蜜桃在线观看| 国产成人精品久久二区二区91 | 国产成人精品久久久久久| 日本av手机在线免费观看| 亚洲精品一区蜜桃| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 久久久亚洲精品成人影院| 久久久久国产一级毛片高清牌| 日韩一区二区三区影片| 中文天堂在线官网| 好男人视频免费观看在线| 亚洲国产日韩一区二区| 精品第一国产精品| 男女无遮挡免费网站观看| 国产精品一国产av| 性色avwww在线观看| 国产白丝娇喘喷水9色精品| 一级毛片电影观看| 国产xxxxx性猛交| 亚洲伊人久久精品综合| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 国产精品久久久久久精品电影小说| 欧美日本中文国产一区发布| 大片电影免费在线观看免费| 日日啪夜夜爽| 看非洲黑人一级黄片| 春色校园在线视频观看| 18+在线观看网站| 国产成人精品无人区| 欧美日韩亚洲国产一区二区在线观看 | 亚洲婷婷狠狠爱综合网| 亚洲国产精品999| 大片电影免费在线观看免费| 国产深夜福利视频在线观看| 丝袜人妻中文字幕| 亚洲四区av| 久久久国产精品麻豆| 国产日韩一区二区三区精品不卡| 亚洲少妇的诱惑av| 午夜av观看不卡| 亚洲五月色婷婷综合| 制服丝袜香蕉在线| 丝袜美足系列| 18禁国产床啪视频网站| 欧美亚洲 丝袜 人妻 在线| 国产白丝娇喘喷水9色精品| 美女国产视频在线观看| 国产精品国产三级国产专区5o| 99国产综合亚洲精品| 9色porny在线观看| 成年人午夜在线观看视频| 999精品在线视频| 你懂的网址亚洲精品在线观看| 中国国产av一级| 久久久久久久亚洲中文字幕| 国产精品 国内视频| 最近手机中文字幕大全| 国产成人精品一,二区| 成人毛片60女人毛片免费| 亚洲婷婷狠狠爱综合网| 中文字幕人妻丝袜一区二区 | 亚洲av.av天堂| 免费看av在线观看网站| 精品酒店卫生间| 香蕉丝袜av| 97在线视频观看| 男女无遮挡免费网站观看| 亚洲综合色网址| 一个人免费看片子| 国产97色在线日韩免费| 日日啪夜夜爽| 色婷婷av一区二区三区视频| 国产深夜福利视频在线观看| 精品国产一区二区久久| 国产精品久久久久久精品电影小说| 国产av精品麻豆| 毛片一级片免费看久久久久| 91精品三级在线观看| 男女国产视频网站| av视频免费观看在线观看| 男女高潮啪啪啪动态图| 国产又色又爽无遮挡免| 亚洲欧洲精品一区二区精品久久久 | 久久人人爽av亚洲精品天堂| 亚洲视频免费观看视频| 色吧在线观看| 欧美av亚洲av综合av国产av | 亚洲av电影在线观看一区二区三区| 另类亚洲欧美激情| 中文精品一卡2卡3卡4更新| 国产日韩欧美在线精品| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久久久免| 青草久久国产| 999精品在线视频| 日本午夜av视频| 侵犯人妻中文字幕一二三四区| 亚洲第一区二区三区不卡| videossex国产| 国产成人精品福利久久| 欧美日韩视频高清一区二区三区二| 99国产精品免费福利视频| 男人操女人黄网站| 免费在线观看视频国产中文字幕亚洲 | 日本黄色日本黄色录像| 青春草亚洲视频在线观看| 亚洲人成网站在线观看播放| 国产激情久久老熟女| tube8黄色片| 一本—道久久a久久精品蜜桃钙片| 伊人久久国产一区二区| 国产极品粉嫩免费观看在线| 成人二区视频| 日韩不卡一区二区三区视频在线| 欧美97在线视频| 美女xxoo啪啪120秒动态图| 成人亚洲精品一区在线观看| 下体分泌物呈黄色| av免费观看日本| 日本91视频免费播放| av一本久久久久| 国产在线免费精品| 国产老妇伦熟女老妇高清| 成人手机av| 97在线人人人人妻| www.熟女人妻精品国产| 久久精品人人爽人人爽视色| 女性被躁到高潮视频| 亚洲人成电影观看| 国产乱人偷精品视频| 亚洲精品国产一区二区精华液| 亚洲成av片中文字幕在线观看 | 亚洲精品美女久久久久99蜜臀 | 国产成人a∨麻豆精品| 在线看a的网站| 亚洲三区欧美一区| 久久久久久伊人网av| 成年动漫av网址| 最近中文字幕高清免费大全6| 丝袜美腿诱惑在线| 青春草国产在线视频| 精品久久久久久电影网| 国产午夜精品一二区理论片| 亚洲成人手机| 麻豆av在线久日| 久久人人97超碰香蕉20202| 2018国产大陆天天弄谢| 免费大片黄手机在线观看| 免费日韩欧美在线观看| 成人黄色视频免费在线看| 欧美精品高潮呻吟av久久| 国产高清国产精品国产三级| 999久久久国产精品视频| 少妇被粗大猛烈的视频| 国产免费又黄又爽又色| 亚洲av电影在线观看一区二区三区| 国产白丝娇喘喷水9色精品| a 毛片基地| 婷婷色综合大香蕉| 男女无遮挡免费网站观看| av在线播放精品| 满18在线观看网站| 国产精品蜜桃在线观看| 日韩av免费高清视频| 欧美老熟妇乱子伦牲交| 午夜福利影视在线免费观看| 91国产中文字幕| 丝袜喷水一区| 国产成人av激情在线播放| 国产精品av久久久久免费| 侵犯人妻中文字幕一二三四区| 精品一区二区三区四区五区乱码 | 日韩中文字幕欧美一区二区 | 又大又黄又爽视频免费| 日韩欧美一区视频在线观看| 毛片一级片免费看久久久久| 亚洲精品视频女| 免费黄网站久久成人精品| 日韩av免费高清视频| 欧美激情极品国产一区二区三区| 中文乱码字字幕精品一区二区三区| 五月伊人婷婷丁香| 久久婷婷青草| 丝袜美足系列| 汤姆久久久久久久影院中文字幕| 亚洲少妇的诱惑av| 黑人欧美特级aaaaaa片| 永久网站在线| 亚洲国产成人一精品久久久| √禁漫天堂资源中文www| 90打野战视频偷拍视频| 久久精品夜色国产| 久久这里只有精品19| kizo精华| 免费在线观看视频国产中文字幕亚洲 | av女优亚洲男人天堂| 少妇被粗大猛烈的视频| 精品人妻在线不人妻| 国产免费一区二区三区四区乱码| 超色免费av| 久久久久久久亚洲中文字幕| 国产男女内射视频| 亚洲欧美成人精品一区二区| 欧美另类一区| videossex国产| 99国产综合亚洲精品| 毛片一级片免费看久久久久| 90打野战视频偷拍视频| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线观看免费完整高清在| 丝袜在线中文字幕| 男女边吃奶边做爰视频| 99精国产麻豆久久婷婷| 婷婷成人精品国产| 国产亚洲午夜精品一区二区久久| 色网站视频免费| 亚洲成人一二三区av| 日韩av免费高清视频| 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| xxx大片免费视频| 伦理电影大哥的女人| 狂野欧美激情性bbbbbb| 性少妇av在线| 欧美人与善性xxx| 曰老女人黄片| 最黄视频免费看| 国产又色又爽无遮挡免| 2021少妇久久久久久久久久久| 日韩中字成人| 久久久欧美国产精品| 欧美 亚洲 国产 日韩一| 99热国产这里只有精品6| 国产不卡av网站在线观看| 亚洲欧美清纯卡通| av免费观看日本| 久久毛片免费看一区二区三区| 中文字幕精品免费在线观看视频| 在线观看免费日韩欧美大片| 巨乳人妻的诱惑在线观看| 热re99久久国产66热| 两性夫妻黄色片| 午夜福利乱码中文字幕| 国语对白做爰xxxⅹ性视频网站| 哪个播放器可以免费观看大片| 午夜免费观看性视频| 国产在视频线精品| 免费在线观看黄色视频的| 1024香蕉在线观看| 中文乱码字字幕精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄|