• <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
    国产成人av激情在线播放| 免费无遮挡裸体视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲成av片中文字幕在线观看| 亚洲三区欧美一区| 一区二区三区国产精品乱码| 色av中文字幕| 欧美成人免费av一区二区三区| 国产精品av久久久久免费| 日韩欧美国产一区二区入口| 成人特级黄色片久久久久久久| 男男h啪啪无遮挡| 又大又爽又粗| 国产精品乱码一区二三区的特点| 老鸭窝网址在线观看| 亚洲国产看品久久| 午夜激情av网站| 99久久99久久久精品蜜桃| 国产久久久一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 老熟妇乱子伦视频在线观看| 国产成人av激情在线播放| 一夜夜www| 99精品欧美一区二区三区四区| 一进一出好大好爽视频| 91成年电影在线观看| 免费观看人在逋| 禁无遮挡网站| 日韩一卡2卡3卡4卡2021年| 亚洲精品美女久久av网站| 搡老妇女老女人老熟妇| 免费在线观看黄色视频的| 亚洲第一欧美日韩一区二区三区| 国产成人欧美在线观看| 淫妇啪啪啪对白视频| 无限看片的www在线观看| 精品不卡国产一区二区三区| 午夜影院日韩av| 在线永久观看黄色视频| 国产亚洲欧美在线一区二区| 国产成人av教育| 精品久久久久久成人av| 精品国产国语对白av| 国产精品1区2区在线观看.| 亚洲中文字幕一区二区三区有码在线看 | 欧美日本亚洲视频在线播放| 99精品在免费线老司机午夜| 在线观看一区二区三区| 97碰自拍视频| 亚洲一区二区三区色噜噜| 亚洲国产精品合色在线| 亚洲av片天天在线观看| 国产在线精品亚洲第一网站| 老司机福利观看| 又大又爽又粗| 最近最新中文字幕大全免费视频| 亚洲精品久久国产高清桃花| 久久久水蜜桃国产精品网| 欧美激情 高清一区二区三区| 黑人巨大精品欧美一区二区mp4| 成人国产综合亚洲| 久久久久久亚洲精品国产蜜桃av| 少妇的丰满在线观看| 两个人看的免费小视频| 777久久人妻少妇嫩草av网站| 亚洲国产高清在线一区二区三 | 国产精品影院久久| 国产免费av片在线观看野外av| 亚洲国产看品久久| 午夜视频精品福利| 宅男免费午夜| 51午夜福利影视在线观看| 制服诱惑二区| 日本一本二区三区精品| 国产亚洲精品综合一区在线观看 | www国产在线视频色| 精品第一国产精品| 国产黄片美女视频| 欧美又色又爽又黄视频| 国产91精品成人一区二区三区| 一夜夜www| 欧美人与性动交α欧美精品济南到| 两人在一起打扑克的视频| 级片在线观看| 国产亚洲欧美精品永久| 99久久99久久久精品蜜桃| 亚洲 欧美 日韩 在线 免费| 国产高清激情床上av| 精品乱码久久久久久99久播| 国产精品久久久久久人妻精品电影| av在线播放免费不卡| 亚洲熟妇熟女久久| 久久九九热精品免费| 精品午夜福利视频在线观看一区| 国产真实乱freesex| 欧美大码av| 女人爽到高潮嗷嗷叫在线视频| 亚洲午夜精品一区,二区,三区| 日日摸夜夜添夜夜添小说| 色播在线永久视频| 一区二区三区激情视频| 国产又色又爽无遮挡免费看| 免费女性裸体啪啪无遮挡网站| 久久精品91无色码中文字幕| 校园春色视频在线观看| 99国产精品一区二区三区| 亚洲av片天天在线观看| 男女午夜视频在线观看| 国产一区二区在线av高清观看| 18禁美女被吸乳视频| 免费搜索国产男女视频| svipshipincom国产片| 精品卡一卡二卡四卡免费| 白带黄色成豆腐渣| 欧美 亚洲 国产 日韩一| 亚洲 国产 在线| 国产高清激情床上av| 日韩av在线大香蕉| 亚洲免费av在线视频| 一区二区三区国产精品乱码| 国产男靠女视频免费网站| 久久天堂一区二区三区四区| 真人一进一出gif抽搐免费| 99热6这里只有精品| 国产精品电影一区二区三区| 日本五十路高清| 久热这里只有精品99| 两个人视频免费观看高清| 亚洲免费av在线视频| 日日爽夜夜爽网站| 12—13女人毛片做爰片一| 亚洲中文字幕日韩| 久热这里只有精品99| 午夜福利视频1000在线观看| 法律面前人人平等表现在哪些方面| 欧美大码av| 岛国在线观看网站| ponron亚洲| 成人av一区二区三区在线看| 男人舔女人下体高潮全视频| 国内揄拍国产精品人妻在线 | 非洲黑人性xxxx精品又粗又长| 老司机午夜福利在线观看视频| 黄色丝袜av网址大全| 一区二区三区国产精品乱码| 午夜视频精品福利| a级毛片a级免费在线| 91麻豆精品激情在线观看国产| 日韩欧美国产在线观看| 国产精品亚洲av一区麻豆| 在线观看午夜福利视频| 亚洲精品美女久久久久99蜜臀| 日韩成人在线观看一区二区三区| 久久精品国产综合久久久| 日日摸夜夜添夜夜添小说| 搡老熟女国产l中国老女人| 亚洲av片天天在线观看| 老司机在亚洲福利影院| 亚洲专区字幕在线| 此物有八面人人有两片| 久久久水蜜桃国产精品网| 欧美成狂野欧美在线观看| 国产av一区在线观看免费| 激情在线观看视频在线高清| 国产精品香港三级国产av潘金莲| 在线观看一区二区三区| 欧美在线一区亚洲| 男女那种视频在线观看| 日本黄色视频三级网站网址| 日韩欧美三级三区| 午夜福利视频1000在线观看| 精品一区二区三区av网在线观看| 香蕉久久夜色| 国产男靠女视频免费网站| 国产精品爽爽va在线观看网站 | 大型av网站在线播放| 国产区一区二久久| 99久久99久久久精品蜜桃| 日韩精品青青久久久久久| 国产一区二区在线av高清观看| 夜夜爽天天搞| 日本在线视频免费播放| 麻豆成人av在线观看| 午夜a级毛片| 制服丝袜大香蕉在线| 国内精品久久久久久久电影| 一本综合久久免费| 2021天堂中文幕一二区在线观 | 俺也久久电影网| 熟女少妇亚洲综合色aaa.| 99久久99久久久精品蜜桃| 天天躁夜夜躁狠狠躁躁| 女生性感内裤真人,穿戴方法视频| 日日干狠狠操夜夜爽| 精品一区二区三区四区五区乱码| 久久午夜综合久久蜜桃| 亚洲av五月六月丁香网| 午夜福利高清视频| 国产av在哪里看| 2021天堂中文幕一二区在线观 | 女人被狂操c到高潮| 色播在线永久视频| 亚洲 国产 在线| 久久婷婷人人爽人人干人人爱| 久久亚洲精品不卡| 日本黄色视频三级网站网址| 1024香蕉在线观看| 男女视频在线观看网站免费 | 欧美色视频一区免费| 亚洲天堂国产精品一区在线| 亚洲中文字幕日韩| 亚洲成av片中文字幕在线观看| 色精品久久人妻99蜜桃| 久久久久久久久免费视频了| 亚洲成人免费电影在线观看| 18禁国产床啪视频网站| 久久久久久亚洲精品国产蜜桃av| 精品不卡国产一区二区三区| av视频在线观看入口| 日本一本二区三区精品| 国产一区在线观看成人免费| 国产成人一区二区三区免费视频网站| 国产黄片美女视频| 午夜福利免费观看在线| 黄片播放在线免费| 国产又色又爽无遮挡免费看| 99久久精品国产亚洲精品| 久热爱精品视频在线9| 久久久国产精品麻豆| 国产av一区在线观看免费| 熟女电影av网| 国产精品二区激情视频| 黄频高清免费视频| 国内久久婷婷六月综合欲色啪| 精品国产国语对白av| 国产99白浆流出| 长腿黑丝高跟| 非洲黑人性xxxx精品又粗又长| 午夜久久久在线观看| 欧美乱色亚洲激情| 午夜免费成人在线视频| 亚洲一码二码三码区别大吗| 美女扒开内裤让男人捅视频| 精品福利观看| 亚洲av美国av| 国产麻豆成人av免费视频| 一级毛片精品| 老司机靠b影院| 90打野战视频偷拍视频| x7x7x7水蜜桃| 50天的宝宝边吃奶边哭怎么回事| 欧美精品亚洲一区二区| 在线观看一区二区三区| 久久久久久国产a免费观看| 久久久久九九精品影院| 1024手机看黄色片| 成人亚洲精品av一区二区| 欧美日韩精品网址| 在线视频色国产色| 99国产综合亚洲精品| 别揉我奶头~嗯~啊~动态视频| 国产高清有码在线观看视频 | 国产aⅴ精品一区二区三区波| 人人妻人人看人人澡| 黑人巨大精品欧美一区二区mp4| 国产av一区二区精品久久| 在线观看午夜福利视频| 中文字幕久久专区| 亚洲一区二区三区色噜噜| 久久中文看片网| 亚洲成av人片免费观看| 精品国产国语对白av| 国产精品综合久久久久久久免费| 成人三级黄色视频| 首页视频小说图片口味搜索| 国产av一区在线观看免费| 麻豆成人av在线观看| 国产精品,欧美在线| 亚洲一码二码三码区别大吗| 美女扒开内裤让男人捅视频| 久久精品91蜜桃| 一二三四社区在线视频社区8| 91国产中文字幕| 久久国产精品人妻蜜桃| 久久久久久九九精品二区国产 | 亚洲精品中文字幕一二三四区| 国内精品久久久久久久电影| 亚洲av片天天在线观看| 久久精品国产亚洲av香蕉五月| 黄频高清免费视频| 午夜成年电影在线免费观看| 精品不卡国产一区二区三区| 亚洲精品美女久久av网站| 精品一区二区三区四区五区乱码| 亚洲欧美日韩高清在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩中文字幕欧美一区二区| 亚洲美女黄片视频| 色综合站精品国产| 久久草成人影院| 日韩精品中文字幕看吧| 麻豆久久精品国产亚洲av| 丁香欧美五月| 国产午夜精品久久久久久| 丝袜人妻中文字幕| 色综合欧美亚洲国产小说| 久久久久精品国产欧美久久久| 女警被强在线播放| 亚洲欧美精品综合久久99| 99久久综合精品五月天人人| 91国产中文字幕| 午夜老司机福利片| 日本一区二区免费在线视频| 国产区一区二久久| 女人被狂操c到高潮| 欧美最黄视频在线播放免费| 久久国产精品影院| 亚洲国产精品sss在线观看| 亚洲国产欧洲综合997久久, | 岛国视频午夜一区免费看| 99在线人妻在线中文字幕| av欧美777| av超薄肉色丝袜交足视频| 日韩大码丰满熟妇| 国产精品亚洲av一区麻豆| 国产在线精品亚洲第一网站| 又黄又爽又免费观看的视频| 日本熟妇午夜| 悠悠久久av| 国产精品香港三级国产av潘金莲| 免费av毛片视频| 好男人电影高清在线观看| 欧美在线黄色| 在线永久观看黄色视频| 黄色视频,在线免费观看| 一夜夜www| 久久久国产成人免费| 欧美成人免费av一区二区三区| 久久人妻av系列| 国产av一区在线观看免费| 日日夜夜操网爽| 黑人巨大精品欧美一区二区mp4| 老熟妇仑乱视频hdxx| 久久久久亚洲av毛片大全| 18禁国产床啪视频网站| 成人免费观看视频高清| 极品教师在线免费播放| 后天国语完整版免费观看| 一级毛片精品| 亚洲熟妇熟女久久| 亚洲精华国产精华精| 一本精品99久久精品77| 一区二区日韩欧美中文字幕| 日韩欧美在线二视频| 国产一区在线观看成人免费| 国产成人啪精品午夜网站| 巨乳人妻的诱惑在线观看| 亚洲一区高清亚洲精品| 欧美av亚洲av综合av国产av| 免费看日本二区| 2021天堂中文幕一二区在线观 | 又黄又粗又硬又大视频| 日韩欧美免费精品| 欧美绝顶高潮抽搐喷水| 亚洲精品中文字幕一二三四区| 真人做人爱边吃奶动态| 亚洲狠狠婷婷综合久久图片| 女人爽到高潮嗷嗷叫在线视频| 国产精品日韩av在线免费观看| 国产av又大| 欧美亚洲日本最大视频资源| 真人一进一出gif抽搐免费| 亚洲中文字幕一区二区三区有码在线看 | 悠悠久久av| 最近在线观看免费完整版| 悠悠久久av| 在线观看日韩欧美| 嫩草影视91久久| 国产精品美女特级片免费视频播放器 | 两个人看的免费小视频| 欧美黑人欧美精品刺激| 免费观看精品视频网站| 国产三级黄色录像| 亚洲国产欧美网| 国产欧美日韩一区二区三| 久久久久精品国产欧美久久久| 一区福利在线观看| 国产亚洲欧美98| 亚洲自拍偷在线| 国产区一区二久久| 国产不卡一卡二| 日韩欧美在线二视频| 国产蜜桃级精品一区二区三区| 国产精品久久电影中文字幕| 国产91精品成人一区二区三区| 国产久久久一区二区三区| 色综合亚洲欧美另类图片| 国产aⅴ精品一区二区三区波| 日韩欧美一区二区三区在线观看| 午夜成年电影在线免费观看| 欧美激情 高清一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 无遮挡黄片免费观看| 久久久久国产一级毛片高清牌| 神马国产精品三级电影在线观看 | 国产片内射在线| 久久精品国产清高在天天线| 91大片在线观看| 亚洲av美国av| 男男h啪啪无遮挡| 99国产综合亚洲精品| 免费在线观看视频国产中文字幕亚洲| 午夜亚洲福利在线播放| 国产精品一区二区精品视频观看| 夜夜爽天天搞| 亚洲精品在线美女| 在线观看日韩欧美| 91成人精品电影| 特大巨黑吊av在线直播 | 九色国产91popny在线| 亚洲av美国av| 亚洲中文日韩欧美视频| 中文字幕另类日韩欧美亚洲嫩草| 一本综合久久免费| 人成视频在线观看免费观看| 欧美成人午夜精品| 麻豆国产av国片精品| 国产爱豆传媒在线观看 | 免费观看精品视频网站| 天天添夜夜摸| 欧美日韩瑟瑟在线播放| 午夜福利在线观看吧| 欧美成人午夜精品| 欧美日韩福利视频一区二区| 色在线成人网| 久久精品国产99精品国产亚洲性色| 亚洲精品在线观看二区| av片东京热男人的天堂| 国产精品香港三级国产av潘金莲| 亚洲av电影不卡..在线观看| 色综合亚洲欧美另类图片| 色av中文字幕| 国产成年人精品一区二区| 日韩欧美国产在线观看| av有码第一页| 欧美性长视频在线观看| 俺也久久电影网| 成人av一区二区三区在线看| 首页视频小说图片口味搜索| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 亚洲av电影在线进入| 日本熟妇午夜| 大香蕉久久成人网| 91字幕亚洲| 国产区一区二久久| 在线永久观看黄色视频| 国产精品野战在线观看| 一夜夜www| 亚洲五月色婷婷综合| 亚洲成av人片免费观看| 女人高潮潮喷娇喘18禁视频| √禁漫天堂资源中文www| 精品乱码久久久久久99久播| 色播在线永久视频| 一个人观看的视频www高清免费观看 | 91国产中文字幕| 亚洲专区中文字幕在线| 性欧美人与动物交配| 午夜福利一区二区在线看| 男女那种视频在线观看| 亚洲性夜色夜夜综合| 女性被躁到高潮视频| 视频在线观看一区二区三区| 男人舔女人的私密视频| 精品久久蜜臀av无| 亚洲精品色激情综合| 看片在线看免费视频| 91av网站免费观看| 久久国产精品影院| 无人区码免费观看不卡| 美女午夜性视频免费| 欧美大码av| 亚洲av日韩精品久久久久久密| 精品卡一卡二卡四卡免费| 一夜夜www| 国产精品1区2区在线观看.| 一个人观看的视频www高清免费观看 | 亚洲电影在线观看av| 久久精品国产亚洲av高清一级| 国产亚洲欧美精品永久| 亚洲成人精品中文字幕电影| 亚洲人成伊人成综合网2020| 国产一级毛片七仙女欲春2 | 夜夜躁狠狠躁天天躁| 在线永久观看黄色视频| √禁漫天堂资源中文www| 欧美在线黄色| 国内精品久久久久久久电影| 大型黄色视频在线免费观看| 中文字幕人妻熟女乱码| 身体一侧抽搐| 最近在线观看免费完整版| 999久久久精品免费观看国产| 日本 欧美在线| 欧美日韩精品网址| 亚洲aⅴ乱码一区二区在线播放 | 欧美精品啪啪一区二区三区| 香蕉国产在线看| 亚洲美女黄片视频| 黑人欧美特级aaaaaa片| 黄色视频,在线免费观看| 在线观看66精品国产| 午夜精品久久久久久毛片777| 国产三级在线视频| 亚洲第一电影网av| 国产黄a三级三级三级人| 欧美精品啪啪一区二区三区| 两性夫妻黄色片| 午夜福利在线在线| 午夜成年电影在线免费观看| 女性生殖器流出的白浆| 国产精品1区2区在线观看.| 日本一区二区免费在线视频| 欧美黑人精品巨大| 一区二区三区国产精品乱码| 色精品久久人妻99蜜桃| 久久久久亚洲av毛片大全| 丰满人妻熟妇乱又伦精品不卡| 老汉色av国产亚洲站长工具| www日本在线高清视频| 男人操女人黄网站| 国产又色又爽无遮挡免费看| 日日爽夜夜爽网站| 一级片免费观看大全| 一进一出抽搐gif免费好疼| 亚洲欧美精品综合久久99| 亚洲精品美女久久久久99蜜臀| 亚洲av成人一区二区三| 久久精品91无色码中文字幕| 91字幕亚洲| 丰满人妻熟妇乱又伦精品不卡| 变态另类丝袜制服| 757午夜福利合集在线观看| 久久久国产成人免费| 国产又色又爽无遮挡免费看| 国产1区2区3区精品| 老鸭窝网址在线观看| 国产精品综合久久久久久久免费| 国产v大片淫在线免费观看| 88av欧美| www日本在线高清视频| 久9热在线精品视频| 一级a爱片免费观看的视频| 久久久久久久久免费视频了| 亚洲欧美一区二区三区黑人| 一进一出抽搐gif免费好疼| 精品久久蜜臀av无| 亚洲精品美女久久久久99蜜臀| 国产精品亚洲美女久久久| 香蕉国产在线看| 啦啦啦免费观看视频1| 日日爽夜夜爽网站| 1024视频免费在线观看| 波多野结衣高清无吗| 中文在线观看免费www的网站 | 制服诱惑二区| 精品国产一区二区三区四区第35| 黄片播放在线免费| 国产黄色小视频在线观看| 悠悠久久av| 国产成人啪精品午夜网站| 国产黄a三级三级三级人| 在线十欧美十亚洲十日本专区| 精品欧美国产一区二区三| 一a级毛片在线观看| 一级毛片女人18水好多| 亚洲成av人片免费观看| 波多野结衣av一区二区av| 最新美女视频免费是黄的| 琪琪午夜伦伦电影理论片6080| 好看av亚洲va欧美ⅴa在| 又大又爽又粗| 午夜免费成人在线视频| 最新美女视频免费是黄的| 香蕉丝袜av| 婷婷亚洲欧美| 很黄的视频免费| 国产成年人精品一区二区| 亚洲一区中文字幕在线| x7x7x7水蜜桃| 国产精品,欧美在线| 亚洲va日本ⅴa欧美va伊人久久| 91国产中文字幕| 国产又色又爽无遮挡免费看| 午夜日韩欧美国产| 久久九九热精品免费| 久久精品夜夜夜夜夜久久蜜豆 | 在线观看舔阴道视频| tocl精华| 久久99热这里只有精品18| 久久精品亚洲精品国产色婷小说| 国产精品久久视频播放| 亚洲成人精品中文字幕电影| 亚洲精品av麻豆狂野| 免费看十八禁软件| 亚洲第一av免费看| 免费一级毛片在线播放高清视频| 亚洲精华国产精华精| 无遮挡黄片免费观看| 国产激情偷乱视频一区二区| 1024视频免费在线观看| 在线视频色国产色| 波多野结衣巨乳人妻| 给我免费播放毛片高清在线观看|