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

    新概念高速穿梭艇系列船型及其直航性能

    2017-01-02 13:56:29魏成柱易宏李英輝
    中國(guó)艦船研究 2017年2期
    關(guān)鍵詞:干舷片體雙體船

    魏成柱,易宏,李英輝

    1上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200240

    2高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240

    新概念高速穿梭艇系列船型及其直航性能

    魏成柱1,2,易宏1,2,李英輝1

    1上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200240

    2高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240

    [目的]為了將高速與高耐波性能相結(jié)合,上海交通大學(xué)開發(fā)了高速穿梭艇系列復(fù)合船型,目前已發(fā)展出單體、雙體和三體船型。[方法]介紹高速穿梭艇系列船型的新進(jìn)展和船型設(shè)計(jì)特點(diǎn),并通過數(shù)值水池實(shí)驗(yàn)對(duì)高速穿梭艇系列船型在靜水中的直航性能進(jìn)行研究。數(shù)值水池通過求解URANS方程和采用重疊網(wǎng)格技術(shù)來預(yù)報(bào)船體受力和運(yùn)動(dòng)。[結(jié)果]數(shù)值水池實(shí)驗(yàn)結(jié)果表明,高速穿梭艇系列船型具有優(yōu)良的快速性和航行姿態(tài),船體興波和飛濺隨船型的不同而有所差異。[結(jié)論]展示了多種創(chuàng)新的船舶設(shè)計(jì)方案,可為研究人員提供定量和定性的參考。

    高速穿梭艇;船型開發(fā);高性能船舶;直航性能;數(shù)值仿真

    0 引 言

    高性能船舶的發(fā)展圍繞著不斷提高快速性和耐波性展開。自上世紀(jì)40年代開始,出現(xiàn)了多種經(jīng)典的高性能船型,其中的典型代表有水翼艇、氣墊船、深V型船、滑行艇、多體船、穿浪雙體船、槽道滑行艇、小水線面雙體船等,大體上可以分為雙(多)體船型、動(dòng)力增升船型和復(fù)合式船型。這些船型的出現(xiàn)促進(jìn)了艦船的發(fā)展,在船舶發(fā)展史上占有重要的位置。包括上海交通大學(xué)在內(nèi)的眾多科研機(jī)構(gòu)正致力于開發(fā)更多高性能船舶以滿足各領(lǐng)域的需求。

    高速穿梭艇(High speed shuttle vessel of SJTU,SV-SJTU)是上海交通大學(xué)開發(fā)的高速高耐波系列船型,屬于中小型快艇,目前已發(fā)展出高速穿梭艇單體船型(High speed monohull shuttle vessel of SJTU,SV-SJTU-M)、高速穿梭艇雙體船型(High speed catamaran shuttle vessel of SJTU,SV-SJTU-C)和高速穿梭艇三體船型(High speed trimaran shuttle vessel of SJTU,SV-SJTU-T)。高速穿梭艇屬于復(fù)合船型,設(shè)計(jì)上吸收了多種船型的特點(diǎn)。例如:高速穿梭艇單體船型融合了穿浪船型、深V船型、滑行(半滑行)船型的船型特征;雙體船型融合了雙體滑行艇、雙體穿浪船、單體穿浪船的船型特征;三體船型融合了穿浪船型、深V船型、滑行(半滑行)船型、三體船型的船型特征。不同的船型方案具有不同的靜水力和水動(dòng)力特性。穿梭艇追求船體良好的快速性、耐波性和良好的航行姿態(tài),并兼顧船體的隱身特性。高速穿梭艇船型設(shè)計(jì)家族化,船型方案具有繼承性,造型設(shè)計(jì)臉譜化。

    本文是高速穿梭艇系列船型研究中的一部分。擬介紹高速穿梭艇不同系列船型的最新研究進(jìn)展和船型設(shè)計(jì)特點(diǎn),并基于數(shù)值水池實(shí)驗(yàn)對(duì)高速穿梭艇系列船型靜水中的直航性能進(jìn)行研究,揭示不同穿梭艇船型靜水中的船體直航水動(dòng)力特性。本文將展示和對(duì)比多種形式高性能船舶的船型設(shè)計(jì)和水動(dòng)力特點(diǎn),以便為研究人員提供定量和定性的參考。

    1 穿梭艇系列船型介紹

    高速穿梭艇屬于中小型快艇,其主要設(shè)計(jì)思想為:采取適當(dāng)設(shè)計(jì),在獲得良好的縱向姿態(tài)的前提下獲得船體平穩(wěn)抬升,減小船體靜水阻力;采取適當(dāng)設(shè)計(jì),保證良好的船體破浪姿態(tài),實(shí)現(xiàn)波浪中的緩和運(yùn)動(dòng)和優(yōu)良的快速性。高速穿梭艇在單體船型的基礎(chǔ)上,吸收新的船型特征,衍生出了雙體、三體船型。雙體和三體船型相比單體船型,船型特征既有繼承性又融合了新的設(shè)計(jì)元素。表1為穿梭艇典型參考方案的部分船體主尺度。

    1.1 高速穿梭艇單體船型

    高速穿梭艇單體船型最先得到開發(fā),先后進(jìn)行了概念設(shè)計(jì)[1]和船型特征研究[2]。圖1所示為2個(gè)早期的高速穿梭艇單體船型原型設(shè)計(jì)。在原型設(shè)計(jì)的基礎(chǔ)上,進(jìn)一步開發(fā)了新一代高速穿梭艇單體船型。原型方案主船體采用多折角式設(shè)計(jì),新一代的高速穿梭艇單體船型包括尖舭(Hard chine)方案和圓舭(Soft chine)方案。圖2所示為單體船型的船體輪廓和3D渲染圖。

    新一代高速穿梭艇單體船型融合了穿浪船型、深V船型、滑行(半滑行)船型的船型特點(diǎn)。

    穿梭艇單體船型尖舭方案使用瘦削的穿浪船艏設(shè)計(jì),以減小波浪對(duì)船體的擾動(dòng),減小波浪砰擊、船體加速度和波浪中的阻力。干舷和船艏內(nèi)傾,具有隱身特性。船艏和上層建筑采用一體融合式設(shè)計(jì),以減小波浪的砰擊。主船體水下部分采用深V型設(shè)計(jì)。艇艏相對(duì)尖削的快艇在高速時(shí)艏部常常會(huì)插入波浪之中導(dǎo)致“失速”,通常采用加高船艏干舷和增加船艏外飄來克服,但穿浪船艏比較瘦削,不能采用該設(shè)計(jì)。折角型艇與波浪遭遇時(shí),艇艏迅速抬起,可改善船體埋入波浪的現(xiàn)象,中和瘦削穿浪船艏帶來的一些不利影響。為了兼顧水動(dòng)力和耐波性的要求,將底部橫向斜升角沿船長(zhǎng)方向變化,使滑行面發(fā)生扭曲,艇體前部有較大的橫向斜升,從而減小波浪的拍擊。船底橫向斜升角較小可以獲得較大的虛長(zhǎng)度和動(dòng)升力,適合航速高、排水量小的快艇,同時(shí)有利于布置噴水推進(jìn)器。艏部舭龍骨采用下沉式設(shè)計(jì),吃水由艉部向船艏線性加深,形成流暢的船體,并保證船體始終有足夠的浸濕長(zhǎng)度。艉部采用方艉設(shè)計(jì)。船體在中部向船艉收縮,以減小船體剩余阻力。新一代高速穿梭艇吸收了前期半滑行穿浪船和壓浪干舷[3]等研究成果,利用壓浪干舷來控制船體淹濕、增加船體儲(chǔ)備浮力和動(dòng)穩(wěn)性,通過干舷折角設(shè)計(jì)來彌補(bǔ)因干舷內(nèi)傾帶來的容積損失。通過以上措施,快艇在很小的艉傾角下能夠獲得平穩(wěn)的抬升,減小船體濕表面積,減小船體阻力,提高靜水中的快速性,同時(shí)可避免船底暴露于空氣和波浪中,保證船體良好的縱向姿態(tài)來使快艇順利切開波浪,提高耐波性。

    穿梭艇單體船型圓舭方案是由尖舭方案衍生而來,將尖舭船型船底舭部倒角,主要目的是兼顧中低速段的船體阻力。船底曲面采用凸殼式設(shè)計(jì)(橫剖面曲線向船外凸出)可增加船體強(qiáng)度,如圖2(b)所示。

    1.2 高速穿梭艇雙體船型

    高速穿梭艇雙體船型繼承了新一代高速穿梭艇單體船型的設(shè)計(jì)特點(diǎn),并融合了高速雙體滑行艇、高速雙體穿浪船的船型特征,是一種創(chuàng)新的復(fù)合船型,如圖3所示。

    該船型以高速穿梭艇單體船型尖舭方案為母型船,將主船體從中縱剖面處分開形成雙體式設(shè)計(jì),中間加以連接橋連接。船型設(shè)計(jì)同高速雙體槽道滑行艇有一定的相似性,但也存在明顯的差別。同競(jìng)速雙體滑行艇相比,高速穿梭艇雙體船型移除了前段的連接橋,并將上層建筑修改為流線型,連接橋底部設(shè)計(jì)有分水舯來減小波浪中的船體運(yùn)動(dòng)和砰擊。不同于傳統(tǒng)高速槽道滑行艇利用槽道提供動(dòng)升力,雙體高速穿梭艇將槽道加寬加深,減小和避免了波浪對(duì)連接橋的砰擊。高速穿梭艇雙體船型艏部采用穿浪設(shè)計(jì),以減小波浪對(duì)船體的擾動(dòng),上層建筑采用流線型設(shè)計(jì)以便上浪緩和地流過船體。同單體船型相比,雙體船型修改了舭部折角的設(shè)計(jì),外側(cè)干舷沒有采用壓浪式設(shè)計(jì),但為了控制槽道內(nèi)的船體淹濕和上浪,內(nèi)部干舷采用了壓浪式設(shè)計(jì)。與單體船型相比,雙體高速穿梭艇甲板面積增加,初穩(wěn)性增加顯著(表1)。同時(shí),由片體距中縱剖面的距離更遠(yuǎn)所帶來的橫傾回復(fù)力矩的增加可知,雙體船型的動(dòng)穩(wěn)性必將有所增加,但雙體船型的船體初始濕表面積也有所增加。

    1.3 高速穿梭艇三體船型

    研究人員對(duì)三體船研究后認(rèn)為:三體船主體細(xì)長(zhǎng),當(dāng)長(zhǎng)寬比大約在12~18之間時(shí)既能有效降低高速航行時(shí)的興波阻力又不會(huì)由于增加太大的濕表面積而導(dǎo)致摩擦阻力增加;但主體細(xì)長(zhǎng)會(huì)損失一定的初穩(wěn)性,故需要增加左右對(duì)稱的片體。片體排水量不超過主體的10%,其主要作用是彌補(bǔ)初穩(wěn)性損失,增加橫搖阻尼,改善耐波性。但這些結(jié)論多是針對(duì)大型三體船。即使是高速三體船,其對(duì)應(yīng)的傅汝德數(shù)也要小于快艇的傅汝德數(shù)。因此,上述結(jié)論是否適用于高速快艇還有待驗(yàn)證,這里僅作為設(shè)計(jì)參考。

    高速穿梭艇三體船型由單體船型衍生而來,以單體船型為母型船衍生出主船體,并在船體兩側(cè)增加片體來提高船體穩(wěn)性。主船體是由單體船型沿船長(zhǎng)方向拉伸1.5倍而來,并進(jìn)行了局部修改,目的在于利用細(xì)長(zhǎng)體的優(yōu)勢(shì),進(jìn)一步改善航行姿態(tài)、穿浪能力和耐波性。主船體水線處長(zhǎng)寬比約為9.4。未采用更大的長(zhǎng)寬比是考慮到高速快艇體積和總體布置的限制,過分細(xì)長(zhǎng)的船體并不合適。對(duì)應(yīng)于單體船型的尖舭和圓舭方案,主船體分別設(shè)計(jì)了尖舭和圓舭方案。尖舭主船體相比于尖舭單體船型,其底部橫剖面形狀做了一定的修改。如圖4所示,橫剖線外側(cè)折角由原直線平行式設(shè)計(jì)改為了圓弧式設(shè)計(jì)以減小船體飛濺,折角區(qū)域進(jìn)一步向前延伸。

    高速穿梭艇三體船型片體采用對(duì)稱割劃式設(shè)計(jì)來減小阻力、提供橫搖和橫蕩阻尼力,并采用類似壓浪干舷的設(shè)計(jì)來限制高速時(shí)的船體上浪和淹濕。

    通常有2種方式來安排片體和主船體的橫向位置。第1種是通過大量的優(yōu)化計(jì)算尋找片體和主船體能夠形成有利興波干擾的方案,這是學(xué)者通常比較感興趣的學(xué)術(shù)研究?jī)?nèi)容;第2種是將片體遠(yuǎn)離主船體,避免片體和主船體之間發(fā)生干擾。文獻(xiàn)[4]認(rèn)為,在高傅汝德數(shù)區(qū)域片體前置會(huì)減小船體阻力。作者在研究過程中發(fā)現(xiàn),在高速下,如果片體離主船體過近,會(huì)發(fā)生嚴(yán)重的堵水現(xiàn)象,船體阻力也會(huì)隨之增加。本文認(rèn)為,片體前置使其處于主船體的興波范圍之外所帶來的流暢的船體周圍流場(chǎng)是阻力減小的主要原因之一。高速穿梭艇三體船型將片體遠(yuǎn)離主船體,使片體處在船體興波和飛濺的劇烈區(qū)域之外。該措施在避免主船體和片體發(fā)生干擾及高速下堵水現(xiàn)象的同時(shí),也進(jìn)一步改善了船體的穩(wěn)性和橫搖。

    文獻(xiàn)[5]指出,三體船側(cè)體縱向位置位于主船體中部略向后一點(diǎn)時(shí),在大部分航速下耐波性能比較好。根據(jù)文獻(xiàn)[5]中的相關(guān)結(jié)論及作者的研究經(jīng)驗(yàn),將左右對(duì)稱的片體設(shè)置在船體中部,片體浮心和主船體浮心在船長(zhǎng)方向上的位置相同,而并未采用許多三體船片體后置的設(shè)計(jì)。船體布局也參考了高速戰(zhàn)機(jī)的布置。

    因此,從阻力出發(fā),三體船型采用片體遠(yuǎn)離主船體和主船體興波的措施;從耐波性出發(fā),三體船型將片體放置在船體中后部,并將主船體的浮心和片體的浮心進(jìn)行了在船長(zhǎng)方向上的坐標(biāo)對(duì)齊處理。

    2 穿梭艇直航性能研究

    高速穿梭艇直航性能研究基于數(shù)值水池實(shí)驗(yàn)完成。為了完成本次研究,論文主要基于CFD,通過求解URANS方程來研究高速穿梭艇直航水動(dòng)力特征。隨著技術(shù)和學(xué)術(shù)的積累,數(shù)值計(jì)算的精度和可靠性得到很大的提升,在科學(xué)研究和工程應(yīng)用上得到了越來越廣泛的使用和認(rèn)可[6],成為除模型實(shí)驗(yàn)外船型優(yōu)化設(shè)計(jì)非常有效的手段。當(dāng)傅汝德數(shù)較高時(shí),船體存在著較大的縱向和垂向耦合運(yùn)動(dòng)。因此,在數(shù)值計(jì)算中為了更準(zhǔn)確地預(yù)報(bào)船體性能,需要考慮船舶航行運(yùn)動(dòng)改變對(duì)預(yù)報(bào)結(jié)果帶來的影響。獲取和模擬船體運(yùn)動(dòng)最有效的方法為使用動(dòng)網(wǎng)格。為方便對(duì)比,所有船型均縮尺為具有近似尺度和船體重量的模型。高速穿梭艇系列船型直航性能研究基于模型尺度進(jìn)行。計(jì)算所用模型的簡(jiǎn)略主尺度如表2所示。

    2.1 計(jì)算方法及驗(yàn)證

    通過使用CAD建模和CFD技術(shù)搭建了數(shù)值水池。使用重疊網(wǎng)格技術(shù)和六自由度模型來模擬和預(yù)報(bào)船體運(yùn)動(dòng),背景網(wǎng)格使用切割體網(wǎng)格來控制網(wǎng)格密度,減小網(wǎng)格數(shù)量和提高網(wǎng)格的分層性與正交性;隨體網(wǎng)格使用了通過合理布置而獲得良好分層性和正交性的多面體網(wǎng)格。合理布置網(wǎng)格密度及改善網(wǎng)格的正交性和分層性對(duì)消除自由面鋸齒現(xiàn)象及捕捉興波與飛濺非常關(guān)鍵[7]。文獻(xiàn)[7]中的圖7~圖10對(duì)此進(jìn)行了形象的表述。在利用CFD進(jìn)行相關(guān)研究時(shí),通常需要預(yù)先進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證計(jì)算。研究人員進(jìn)行了很多的網(wǎng)格獨(dú)立性研究,研究點(diǎn)主要集中于網(wǎng)格單元數(shù)量對(duì)數(shù)值計(jì)算結(jié)果的影響。由文獻(xiàn)[8]可知,當(dāng)網(wǎng)格單元數(shù)量超過150萬時(shí),其對(duì)計(jì)算結(jié)果的影響就很小了。但是,網(wǎng)格對(duì)計(jì)算結(jié)果的影響不僅來自于網(wǎng)格單元數(shù)量。比如,在網(wǎng)格單元數(shù)量足夠的前提下,網(wǎng)格單元的分布、近壁面網(wǎng)格的劃分和壁面函數(shù)的選擇對(duì)計(jì)算結(jié)果的影響也非常顯著。嚴(yán)格的網(wǎng)格獨(dú)立性研究是非常巨大的工程。為解決網(wǎng)格獨(dú)立性分析問題,作者參考相關(guān)研究結(jié)果直接填充了足夠數(shù)量的網(wǎng)格單元,網(wǎng)格單元數(shù)量大于200萬,根據(jù)預(yù)計(jì)算合理地分布網(wǎng)格密度,并進(jìn)行了y+對(duì)計(jì)算結(jié)果影響的研究。高y+壁面函數(shù)配合較為粗糙的近壁面網(wǎng)格具有足夠的精度。對(duì)壁面函數(shù)的選擇除了依據(jù)作者自己的驗(yàn)證計(jì)算外,文獻(xiàn)[9]也提供了對(duì)壁面函數(shù)選擇的參考。計(jì)算機(jī)硬件的發(fā)展使網(wǎng)格單元數(shù)量不再成為限制數(shù)值仿真計(jì)算的瓶頸。而且較少數(shù)量的網(wǎng)格單元也不能很好地捕捉船體飛濺?;谏鲜鲈?,作者沒有進(jìn)行專門的網(wǎng)格獨(dú)立性計(jì)算驗(yàn)證,而是參考已有的研究成果和網(wǎng)格劃分經(jīng)驗(yàn),一次性填充足夠的網(wǎng)格單元,對(duì)網(wǎng)格進(jìn)行了保守處理。本部分研究所用到的網(wǎng)格單元數(shù)量在200萬~280萬左右,滿足計(jì)算精度對(duì)網(wǎng)格單元數(shù)量的要求。

    求解非穩(wěn)態(tài)RANS方程時(shí)選用結(jié)合了k-ω和k-ε湍流模型特點(diǎn)的SST k-ω[10]湍流模型。研究使用流體體積(Volume ofFluid,VOF)模型捕捉自由面。

    為進(jìn)一步驗(yàn)證計(jì)算精度及保證計(jì)算結(jié)果的可靠性,本文對(duì)數(shù)值仿真計(jì)算方法進(jìn)行了水池實(shí)驗(yàn)驗(yàn)證。該部分研究使用文獻(xiàn)[11]中公開的NACA滑行楔形體和文獻(xiàn)[12]中公開的滑行艇幾何模型及相關(guān)的實(shí)驗(yàn)設(shè)置與實(shí)驗(yàn)數(shù)據(jù)。限于文章篇幅,本節(jié)主要展示基于文獻(xiàn)[12]中的滑行艇的驗(yàn)證結(jié)果,并簡(jiǎn)要介紹基于文獻(xiàn)[11]中的滑行楔形體的驗(yàn)證結(jié)果。文獻(xiàn)[12]中滑行艇的幾何輪廓如圖5所示,模型幾何尺寸如表3所示?;谖墨I(xiàn)[12]中滑行艇的數(shù)值仿真計(jì)算采用自由模式,釋放垂蕩和縱搖2個(gè)自由度;基于NACA滑行楔形體的數(shù)值仿真計(jì)算采用固定模式。計(jì)算設(shè)置、計(jì)算域和網(wǎng)格劃分如圖6所示。

    文獻(xiàn)[11]中的滑行楔形體的水池實(shí)驗(yàn)所測(cè)得的升力平均相對(duì)誤差為3%,最大相對(duì)誤差為8%;測(cè)得的阻力平均相對(duì)誤差為20%,最大相對(duì)誤差為50%。本文基于NACA滑行楔形體的驗(yàn)證計(jì)算的升力誤差平均為5.11%,與實(shí)驗(yàn)值3%的固有平均誤差相比稍微偏大,但誤差均小于最大相對(duì)誤差;而阻力計(jì)算誤差(平均為6.15%)則遠(yuǎn)小于其20%的參考值。文獻(xiàn)[13]和文獻(xiàn)[14]也基于NACA滑行楔形體進(jìn)行了CFD計(jì)算精度的驗(yàn)證,得出了數(shù)值仿真計(jì)算在滑行艇性能預(yù)報(bào)上具有足夠精度的結(jié)論。文獻(xiàn)[13]中阻力最大誤差為20.28%,平均誤差為10.06%;文獻(xiàn)[14]中阻力最大誤差為13.26%,平均誤差為9.65%。與文獻(xiàn)[13]和文獻(xiàn)[14]的計(jì)算結(jié)果相比,本文采用的數(shù)值仿真方法的阻力預(yù)報(bào)精度稍高,而升力預(yù)報(bào)的誤差則稍微偏大。

    由圖7可以看出,基于文獻(xiàn)[12]中公開的滑行艇數(shù)值仿真計(jì)算結(jié)果(CFD)和實(shí)驗(yàn)結(jié)果(EXP)在大小和變化趨勢(shì)上符合得很好。計(jì)算所得阻力的平均相對(duì)誤差為8.98%,最大相對(duì)誤差為14.57%。由數(shù)值水池實(shí)驗(yàn)得到的艉傾角與實(shí)驗(yàn)值相比略微偏大,但最大偏差僅為0.45°。若對(duì)參考文獻(xiàn)[12]中升沉值小數(shù)點(diǎn)后的位數(shù)進(jìn)行數(shù)值仿真計(jì)算值的舍入處理,船體升沉計(jì)算值幾乎和實(shí)驗(yàn)值吻合。由于文獻(xiàn)[11]只給出了船體型線圖,沒有給出船體型值表,考慮到CAD逆向建模的誤差,計(jì)算所用的模型與其實(shí)驗(yàn)所用模型存在一定誤差,再考慮到物理模型、網(wǎng)格離散產(chǎn)生的誤差,可以認(rèn)為本文研究所采用的方法對(duì)阻力和運(yùn)動(dòng)的數(shù)值預(yù)報(bào)精度滿足工程需要。

    同時(shí),考慮基于兩個(gè)模型的驗(yàn)證結(jié)果,并對(duì)比文獻(xiàn)[13]和文獻(xiàn)[14]中NACA滑行楔形體水動(dòng)力預(yù)報(bào)的計(jì)算誤差,可以認(rèn)為本文中數(shù)值水池對(duì)船體阻力的預(yù)報(bào)達(dá)到了很高的精度。因此,可以在同一網(wǎng)格和計(jì)算設(shè)置的前提下進(jìn)行船型水動(dòng)力性能預(yù)報(bào)和對(duì)比。

    2.2 穿梭艇系列船型船體阻力

    本節(jié)展示和對(duì)比高速穿梭艇不同船型的靜水阻力,并分析不同船型的阻力成分,如圖8所示。文獻(xiàn)[15]中公開了M型船(M-HULL)的阻力實(shí)驗(yàn)數(shù)據(jù)。由于其模型重量和尺度與本研究中所用的模型重量和尺度相近,故被用作對(duì)比船型。M型船是美國(guó)M船舶公司(M Ship Company)研發(fā)的新型高性能多用途船型,在快速性、隱身性等方面獨(dú)具特色。

    如圖8所示,采用圓舭設(shè)計(jì)的穿梭艇單體船型和三體船型在低速段的總阻力較小(低速段是相對(duì)于小型高速快艇而言,對(duì)于大型船舶,本文中的傅汝德數(shù)區(qū)間位于高速及超高速段)。在高速段,穿梭艇單體系列船型阻力最小,三體系列阻力最大。當(dāng)體積傅汝德數(shù)F▽=5.5時(shí),阻力大小從高到低依次為 SV-SJTU-M(Hard chine),SVSJTU-M(Soft chine),SV-SJTU-C,SV-SJTU-T(Hard chine),SV-SJTU-T(Soft chine)。尖舭船體方案的阻力在高速段表現(xiàn)出了較大的優(yōu)勢(shì)。當(dāng)F▽=5.5時(shí),SV-SJTU-M(Hard chine)的阻力比SV-SJTU-M(Soft chine)減小了7.1%,SV-SJTU-T(Hard chine)的阻力比SV-SJTU-T(Soft chine)減小了11.4%。與M型船相比,高速穿梭艇系列船型在所計(jì)算的航速段均具有更小的船體阻力。

    由圖8可知,摩擦阻力是高速穿梭艇阻力的主要組成部分。以F▽=5.5為例,摩擦阻力占總阻力的77%~87%。剩余阻力隨速度的變化遠(yuǎn)沒有摩擦阻力劇烈,且并未一直隨速度的增加而增加。對(duì)比5種穿梭艇船型阻力及其成分可知,采用細(xì)長(zhǎng)設(shè)計(jì)的三體穿梭艇的剩余阻力較小,但摩擦阻力較大。雙體船型和三體船型與單體船型相比,阻力增加是由于摩擦阻力增加所導(dǎo)致,而摩擦阻力增加是由于雙體船型和三體船型較單體船型濕表面積增加所致。

    2.3 穿梭艇系列船型船體航行姿態(tài)

    航行姿態(tài)是穿梭艇設(shè)計(jì)時(shí)重點(diǎn)控制的參數(shù)。計(jì)算結(jié)果表明,文中涉及的多種船型的航行姿態(tài)均達(dá)到了在利用設(shè)計(jì)獲得良好縱向姿態(tài)的前提下船體獲得平穩(wěn)抬升的設(shè)計(jì)目標(biāo)(圖9),保證了船體在破波時(shí)能夠使用船艏順利地穿透波浪而不會(huì)用船底面沖擊波面。穿梭艇系列船型在全航速段的艉傾均很小,最大艉傾不超過3°。總體上,三體船型由于船體比較細(xì)長(zhǎng),其船體抬升和艉傾值均小于單體船型。高速穿梭艇系列船型的艉傾角變化趨勢(shì)與滑行艇和M型船不同。如圖9所示,穿梭艇并沒有在低速段出現(xiàn)艉傾峰值,峰值向高速段偏移。三體船型需要更高的航速達(dá)到峰值,在本文的速度范圍內(nèi)沒有出現(xiàn)艉傾峰值。

    2.4 穿梭艇系列船型的船體淹濕和興波

    圖10示出了F▽=3.3及F▽=5.5時(shí)穿梭艇系列船型的船體淹濕分布。圖10中藍(lán)色區(qū)域代表船體濕表面。由于船型之間幾何特征的繼承和衍生關(guān)系,不同穿梭艇船型之間的船體淹濕表現(xiàn)出了相似性和差異性。由圖10可知,穿梭艇均提供了足夠的浸濕長(zhǎng)度來保證耐波性的需求。尖舭方案和圓舭方案在船體淹濕上表現(xiàn)出了差異。在F▽=3.3時(shí),船體還未獲得充分的抬升,此時(shí)壓浪干舷發(fā)揮作用,將船體淹濕限制在了壓浪干舷折角線以下;在F▽=5.5時(shí),船體獲得充分的抬升,尖舭式設(shè)計(jì)使船體淹濕控制在舭部折角線以下,而圓舭船體的船體淹濕仍控制在壓浪干舷折角線以下。雙體船型和三體船型由于其艉傾角小于單體船型,故其舭龍骨浸沒比例要大。

    由于穿梭艇系列船型之間的船型特征差異,其飛濺和興波也存在一定的差異,如圖11所示。

    對(duì)于穿梭艇單體船型圓舭方案,飛濺是由壓浪干舷壓制船體上浪產(chǎn)生的。對(duì)于穿梭艇單體船型尖舭方案,在F▽=3.3時(shí),飛濺是由舭部折角線和壓浪干舷壓制船體上浪產(chǎn)生的,在F▽=5.5時(shí),飛濺是由舭部折角線產(chǎn)生的。在F▽=3.3時(shí),飛濺分布在穿梭艇單體船型船長(zhǎng)范圍內(nèi);隨著航速的提高,興波向船中后部發(fā)展,此時(shí)飛濺分布在船體中后部。

    由于未采用壓浪干舷設(shè)計(jì),穿梭艇雙體船型的飛濺是由舭部的折角產(chǎn)生的。其飛濺特征同穿梭艇單體船型尖舭方案類似。在F▽=3.3時(shí),飛濺分布在船長(zhǎng)范圍內(nèi);隨著航速的提高,興波向船中后部發(fā)展,此時(shí)飛濺分布在船體中后部。觀察槽道內(nèi)的興波可知,槽道內(nèi)的自由面很平坦,槽道內(nèi)流場(chǎng)基本未受到片體的影響。與單體船型不同,雙體船型在船尾處有兩股雞尾流并在船尾遠(yuǎn)場(chǎng)處合并成一股。

    與單體船型和雙體船型不同,2種三體船型在F▽=3.3時(shí)船體上浪和飛濺很小,壓浪干舷沒有發(fā)揮有效的作用。隨著航速的提高,飛濺開始出現(xiàn)。與單體船型相似,在F▽=5.5時(shí),對(duì)于穿梭艇三體船型圓舭方案,飛濺是由壓浪干舷壓制船體上浪產(chǎn)生的;對(duì)于穿梭艇三體船型尖舭方案,飛濺是由舭部折角線產(chǎn)生的。觀察圖11可知,片體處于主船體興波和飛濺區(qū)域之外,片體和主船體之間沒有發(fā)生堵水現(xiàn)象。采用割劃式設(shè)計(jì)的片體對(duì)其周圍流場(chǎng)擾動(dòng)非常小,片體的興波非常小。由于片體的興波微小,因此三體船型尾部興波主要是由主船體造成,故其雞尾流特征同單體船型類似。

    3 結(jié) 論

    作為高速穿梭艇系列船型研究的一部分,本文介紹了穿梭艇系列船型的船型設(shè)計(jì)特點(diǎn)并通過數(shù)值水池實(shí)驗(yàn)研究了系列船型靜水中的直航水動(dòng)力特性。通過以上研究,得到以下主要結(jié)論。

    1)穿梭艇系列船型達(dá)到了在保證船體良好航行姿態(tài)的前提下具有優(yōu)良的快速性的設(shè)計(jì)初衷。

    2)由單體船型衍生出來的雙體船型和三體船型相比于母型船,其橫穩(wěn)性得到了很大提升,橫搖也將得到進(jìn)一步的改善。

    3)相比于穿梭艇雙體船型和三體船型,單體船型在高速段具有阻力小的優(yōu)勢(shì)。穿梭艇雙體船型和三體船型的阻力增加是由濕表面積增加帶來的船體摩擦阻力增加所致。尖舭船型在高速段其阻力要小于圓舭船型,但在低速段阻力要大于圓舭船型。

    4)穿梭艇單體船型興波與穿梭艇雙體船型興波特征存在較大的差異;由于采用割劃式設(shè)計(jì)的片體興波微小,三體船型尾部興波主要由主船體造成,故其雞尾流特征同單體船型類似。

    5)穿梭艇尖舭船型飛濺與圓舭船型飛濺產(chǎn)生的原因不同,尖舭船型的飛濺主要是由舭部折角產(chǎn)生,而圓舭船型的飛濺是由壓浪干舷產(chǎn)生。不同于單體和雙體船型,細(xì)長(zhǎng)的三體船型需要更高的航速才會(huì)產(chǎn)生劇烈的飛濺。

    更多的技術(shù)(如船體斷階、氣泡減阻、可調(diào)式壓浪干舷等)還有待在穿梭艇上進(jìn)行嘗試,穿梭艇系列船型的其他水動(dòng)力性能(如耐波性、波浪中的穩(wěn)性和橫搖、操作性系數(shù)等)還有待研究、整理和發(fā)布,這些工作將在不久以后完成。

    [1] 魏成柱.穿梭艇性能特征與船型優(yōu)化[D].上海:上海交通大學(xué),2013.

    WEI C Z.The shuttle vessel performance characteris?tics and hull form optimization[D].Shanghai:Shang?hai Jiao Tong University,2013(in Chinese).

    [2] 魏成柱,李英輝,易宏.基于CAD與CFD的穿梭艇局部船型特征分析[J].船舶工程,2014,36(3):28-32.

    WEI C Z,LI Y H,YI H.Analysis of shuttle vessel's lo?cal hull form characteristics based on CAD and CFD[J].Ship Engineering,2014,36(3):28-32(in Chinese).

    [3] 魏成柱,毛立夫,李英輝,等.單體半滑行穿浪船船型與靜水航行性能[J].中國(guó)艦船研究,2015,10(5):16-21.

    WEI C Z,MAO L F,LI Y H,et al.Analysis of the hull form and sailing characters in calm water of a semi-planing wave-piercing boat[J].Chinese Journal of Ship Research,2015,10(5):16-21(in Chinese).

    [4] 閆蕾,宗智,倪少玲,等.半滑行前三體船型模型阻力對(duì)比試驗(yàn)研究[J].船海工程,2013 42(1):41-44,49.

    YAN L,ZONG Z,NI S L,et al.Resistance model test?ing of semi-planning forward placed outrigger trimaran hull forms[J].Ship&Ocean Engineering,2013,42(1):41-44,49(in Chinese).

    [5] 張文鵬,宗智,倪少玲,等.三體船耐波性的模型試驗(yàn)研究[J].水動(dòng)力學(xué)研究與進(jìn)展,2007,22(5):619-624.

    ZHANG W P,ZONG Z,NI S L,et al.Model testing of seakeeping performance of trimaran[J].Journal of Hydrodynamics,2007,22(5):619-624(in Chinese).

    [6] LARSSON L,STERN F,VISONNEAU M.Numerical ship hydrodynamics:an assessment of the Gothenburg 2010 Workshop[M].Netherlands:Springer,2013.

    [7] 魏成柱,李英輝,易宏.多面體網(wǎng)格在船體繞流計(jì)算中的應(yīng)用[J].上海交通大學(xué)學(xué)報(bào),2016,50(8):1144-1151.

    WEI C Z,LI Y H,YI H.Application of polyhedral mesh in numerical simulations of flow around hulls[J]. Journal of Shanghai Jiaotong University,2016,50(8):1144-1151(in Chinese).

    [8] YOUSEFIR, SHAFAGHAT R, SHAKERIM. High-speed planing hull drag reduction using tunnels[J].Ocean Engineering,2014,84:54-60.

    [9] 孫華偉,馬偉佳,朱江波.影響滑行艇阻力數(shù)值計(jì)算的網(wǎng)格因素研究[J].中國(guó)造船,2015,56(2):170-178.

    SUN H W,MA W J,ZHU J B.Research on grid factor in numerical calculation of planing craft resistance[J]. Shipbuilding of China,2015,56(2):170-178(in Chinese).

    [10] MENTER F R.Two-equation eddy-viscosity turbu?lence models for engineering applications[J].AIAA Journal,1994,32(8):1598-1605.

    [11] CHAMBLISS D B,BOYD JR G M.The planing char?acteristics of two v-shaped prismatic surfaces having angles of dead rise of 20°and 40°[R].NACA Techni?cal Note 2876.Washington:NACA,1953.

    [12] TAUNTON D J,HUDSON D A,SHENOI R A.Char?acteristics of a series of high speed hard chine planing hulls-part 1:performance in calm water[J].Interna?tional Journal of Small Craft Technology,2010,152:55-75.

    [13] BRIZZOLARA S,SERRA F.Accuracy of CFD codes in the prediction of planing surfaces hydrodynamic characteristics[C]//Proceedings of the 2nd Interna?tional Conference on Marine Research and Transporta?tion.Naples:ICMRT,2007:147-159.

    [14] 王碩,蘇玉民,龐永杰,等.高速滑行艇CFD精度研究[J].船舶力學(xué),2013,17(10):1107-1114.

    WANG S,SU Y M,PANG Y J,et al.Study on the accuracy in the hydrodynamic prediction of high-speed planing crafts of CFD method[J].Journal of Ship Mechanics,2013,17(10):1107-1114(in Chinese).

    [15] 陳輝.M船型水氣兩相流場(chǎng)特性研究[D].武漢:武漢理工大學(xué),2011.

    CHEN H.The research of M-ship's two-phrase flow field characteristic[D].Wuhan:Wuhan University of Technology,2011(in Chinese).

    Hull forms and straight forward CFD free running trials of high-speed shuttle vessels

    WEI Chengzhu1,2,YI Hong1,2,LI Yinghui1
    1 State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China
    2 Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai 200240,China

    SV-SJTU high-speed shuttle vessels are developed by Shanghai Jiao Tong University for fine seakeeping and high speed purposes.A series of SV-SJTUs have been developed,and are introduced in this paper.Straight forward CFD free running trials were conducted and the results are also presented.Hull resistance and motions are predicted by solving URANS equations and adopting the overset mesh method. The results of the straight forward CFD free running trials prove that SV-SJTUs have little resistance and fine hull motion in calm water,and their wave-making and splashing differ with different hull forms.This paper presents the designs of a variety of high performance ships,thereby providing quantitative and quali?tative references for researchers.

    high-speed shuttle vessel;hull form development;high-performance ships;straight forward CFD free running trials;numerical simulations

    U661.3,U662.2

    A

    10.3969/j.issn.1673-3185.2017.02.002

    http://kns.cnki.net/kcms/detail/42.1755.TJ.20170313.1617.034.html

    魏成柱,易宏,李英輝.新概念高速穿梭艇系列船型及其直航性能[J].中國(guó)艦船研究,2017,12(2):12-21.

    WEI C Z,YI H,LI Y H.Hull forms and straight forward CFD free running trials of high-speed shuttle vessels[J]. Chinese Journal of Ship Research,2017,12(2):12-21.

    2016-09-07 < class="emphasis_bold"> 網(wǎng)絡(luò)出版時(shí)間:

    時(shí)間:2017-3-13 16:17

    上海交通大學(xué)海洋工程國(guó)家重點(diǎn)實(shí)驗(yàn)室自主研究課題(GKZD010061)

    魏成柱,男,1987年生,博士生。研究方向:新船型開發(fā),智能水面無人平臺(tái)研發(fā)與數(shù)值仿真計(jì)算。

    E-mail:weichengzhu@sjtu.edu.cn

    易宏(通信作者),男,1962年生,教授,博士生導(dǎo)師。研究方向:潛器與特種船舶開發(fā),海上裝置與系統(tǒng)開發(fā)設(shè)計(jì),系統(tǒng)可靠性與人因工程。E-mail:yihong@sjtu.edu.cn

    李英輝,男,1973年生,博士,講師。研究方向:新型船舶開發(fā)和數(shù)值計(jì)算

    期刊網(wǎng)址:www.ship-research.com

    猜你喜歡
    干舷片體雙體船
    基于有限元方法的滾裝船片體分段疊層放置方案設(shè)計(jì)
    船海工程(2023年4期)2023-08-18 05:42:14
    基于STAR-CCM+的雙體船阻力預(yù)報(bào)
    關(guān)于冬季和北大西洋冬季干舷的探討
    型深建造誤差對(duì)勘劃載重線標(biāo)志的影響
    某大型碼頭模塊的合攏方法及力學(xué)計(jì)算
    圓鋸片受熱時(shí)的應(yīng)力變化規(guī)律
    天津冶金(2018年1期)2018-06-13 08:47:06
    628客位珠江雙體游船的設(shè)計(jì)
    載重線公約的規(guī)定、釋義及應(yīng)用(三)
    片體形式對(duì)WPCat快速性和耐波性的影響
    內(nèi)河雙尾船甲板上浪統(tǒng)計(jì)特性影響因素分析
    廣東造船(2014年4期)2014-05-30 11:47:57
    精品国产国语对白av| 亚洲精品久久久久久婷婷小说| 色网站视频免费| 国产免费福利视频在线观看| 超色免费av| 成人黄色视频免费在线看| 精品人妻偷拍中文字幕| 成人亚洲欧美一区二区av| 九草在线视频观看| 国产精品二区激情视频| 亚洲第一青青草原| 国产成人一区二区在线| 18+在线观看网站| 日韩电影二区| 婷婷成人精品国产| 热99国产精品久久久久久7| 高清av免费在线| 人人妻人人添人人爽欧美一区卜| 如何舔出高潮| 成年人免费黄色播放视频| 国产极品天堂在线| 久久人人爽av亚洲精品天堂| 婷婷成人精品国产| 色94色欧美一区二区| av免费观看日本| 国产人伦9x9x在线观看 | 精品人妻偷拍中文字幕| 又大又黄又爽视频免费| 免费日韩欧美在线观看| 最近最新中文字幕大全免费视频 | 午夜激情久久久久久久| 女人精品久久久久毛片| 午夜免费鲁丝| 色婷婷久久久亚洲欧美| 男人爽女人下面视频在线观看| 成年人午夜在线观看视频| 国产在视频线精品| 免费日韩欧美在线观看| 少妇 在线观看| 咕卡用的链子| 亚洲情色 制服丝袜| 可以免费在线观看a视频的电影网站 | 久久久久久人妻| 国产精品女同一区二区软件| 国产免费一区二区三区四区乱码| 国产成人a∨麻豆精品| 99久国产av精品国产电影| 下体分泌物呈黄色| 日韩一区二区三区影片| 一区二区三区激情视频| 久久久久久人妻| 涩涩av久久男人的天堂| av视频免费观看在线观看| 99re6热这里在线精品视频| 免费大片黄手机在线观看| 日本av免费视频播放| 天美传媒精品一区二区| 日韩制服丝袜自拍偷拍| 亚洲人成77777在线视频| 免费在线观看黄色视频的| 日日摸夜夜添夜夜爱| 欧美日韩亚洲高清精品| 久久国产精品大桥未久av| 日本vs欧美在线观看视频| 男男h啪啪无遮挡| 男人舔女人的私密视频| 在线天堂最新版资源| 日日撸夜夜添| 99久国产av精品国产电影| 久久国产精品男人的天堂亚洲| 中文字幕另类日韩欧美亚洲嫩草| 各种免费的搞黄视频| 一本久久精品| 久久久精品免费免费高清| 伊人亚洲综合成人网| 亚洲欧洲日产国产| 国产乱来视频区| 1024香蕉在线观看| 少妇被粗大猛烈的视频| 久久精品国产综合久久久| 亚洲伊人色综图| 久久国产精品大桥未久av| av片东京热男人的天堂| 曰老女人黄片| 久久毛片免费看一区二区三区| 亚洲av中文av极速乱| 午夜影院在线不卡| 制服丝袜香蕉在线| 国产亚洲最大av| 韩国av在线不卡| 欧美老熟妇乱子伦牲交| 欧美日韩精品成人综合77777| 久久97久久精品| 久久精品国产a三级三级三级| 亚洲av电影在线进入| 精品一区二区免费观看| 老司机影院成人| 午夜福利乱码中文字幕| 在线观看免费日韩欧美大片| 丝袜美足系列| 午夜福利视频在线观看免费| 日本免费在线观看一区| 欧美少妇被猛烈插入视频| 亚洲精品第二区| 国产野战对白在线观看| 99热全是精品| 满18在线观看网站| 男女下面插进去视频免费观看| 日韩人妻精品一区2区三区| 国产黄频视频在线观看| 久久久久视频综合| 国产高清国产精品国产三级| 校园人妻丝袜中文字幕| 热re99久久国产66热| 国产 精品1| 国产爽快片一区二区三区| 80岁老熟妇乱子伦牲交| 精品人妻熟女毛片av久久网站| 国产精品.久久久| 人人妻人人澡人人看| 亚洲国产最新在线播放| 欧美变态另类bdsm刘玥| 一区二区av电影网| 亚洲国产精品成人久久小说| 国产探花极品一区二区| 中国三级夫妇交换| 一边摸一边做爽爽视频免费| 国产成人精品在线电影| 可以免费在线观看a视频的电影网站 | 国产一区二区三区av在线| 黑人欧美特级aaaaaa片| 天天躁夜夜躁狠狠躁躁| 叶爱在线成人免费视频播放| 国产男女内射视频| 日本免费在线观看一区| 日韩中文字幕视频在线看片| 国产精品久久久久成人av| 欧美精品高潮呻吟av久久| 亚洲av日韩在线播放| xxx大片免费视频| 久久久久网色| 黄片小视频在线播放| 久久精品亚洲av国产电影网| 国产成人a∨麻豆精品| 久久精品国产鲁丝片午夜精品| 看十八女毛片水多多多| 亚洲av.av天堂| 不卡av一区二区三区| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区大全| 国产精品女同一区二区软件| 十八禁高潮呻吟视频| 丝袜美腿诱惑在线| √禁漫天堂资源中文www| 一二三四在线观看免费中文在| 久久毛片免费看一区二区三区| 国产精品偷伦视频观看了| 天堂8中文在线网| 丝袜在线中文字幕| 国产高清不卡午夜福利| 少妇人妻 视频| 久久这里只有精品19| 97在线人人人人妻| 欧美成人精品欧美一级黄| 精品99又大又爽又粗少妇毛片| av线在线观看网站| 国产精品麻豆人妻色哟哟久久| 日韩欧美一区视频在线观看| 久久国产精品大桥未久av| 成年美女黄网站色视频大全免费| 亚洲精品aⅴ在线观看| 成人手机av| 极品人妻少妇av视频| 搡老乐熟女国产| 99九九在线精品视频| 一边摸一边做爽爽视频免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久久久久婷婷小说| 国产人伦9x9x在线观看 | 大陆偷拍与自拍| 亚洲一级一片aⅴ在线观看| a级毛片黄视频| 熟女av电影| 一区二区av电影网| 久久久久国产精品人妻一区二区| 岛国毛片在线播放| 亚洲国产欧美在线一区| 自线自在国产av| 国产福利在线免费观看视频| 黄色毛片三级朝国网站| 午夜免费男女啪啪视频观看| 日韩中文字幕视频在线看片| 亚洲美女搞黄在线观看| 日韩av在线免费看完整版不卡| 美女福利国产在线| 男女无遮挡免费网站观看| 亚洲美女黄色视频免费看| 日日摸夜夜添夜夜爱| 精品视频人人做人人爽| 亚洲国产av新网站| 精品国产超薄肉色丝袜足j| 成人亚洲欧美一区二区av| 国产视频首页在线观看| 欧美av亚洲av综合av国产av | 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 精品福利永久在线观看| xxxhd国产人妻xxx| 成人午夜精彩视频在线观看| 人人妻人人澡人人看| 97人妻天天添夜夜摸| 黄色配什么色好看| √禁漫天堂资源中文www| 亚洲美女搞黄在线观看| 街头女战士在线观看网站| 美国免费a级毛片| 欧美亚洲 丝袜 人妻 在线| 日本av手机在线免费观看| 一个人免费看片子| 久久99一区二区三区| 黄色怎么调成土黄色| 乱人伦中国视频| 久久久亚洲精品成人影院| 黄色毛片三级朝国网站| 深夜精品福利| 久久久久国产一级毛片高清牌| 性色av一级| 日本黄色日本黄色录像| av一本久久久久| 亚洲图色成人| 欧美日韩视频精品一区| 99热全是精品| 欧美精品人与动牲交sv欧美| 2021少妇久久久久久久久久久| 黄片小视频在线播放| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 最近手机中文字幕大全| 性少妇av在线| 天美传媒精品一区二区| 蜜桃国产av成人99| 欧美黄色片欧美黄色片| 天天躁夜夜躁狠狠久久av| 国产精品 欧美亚洲| 国产免费视频播放在线视频| 啦啦啦啦在线视频资源| 在线观看免费视频网站a站| 欧美xxⅹ黑人| 亚洲国产av影院在线观看| 9色porny在线观看| 一边摸一边做爽爽视频免费| 日本欧美视频一区| av.在线天堂| 桃花免费在线播放| 午夜久久久在线观看| 亚洲精品日本国产第一区| 国产精品二区激情视频| 丰满乱子伦码专区| 国产一区有黄有色的免费视频| 午夜av观看不卡| 高清欧美精品videossex| 日韩,欧美,国产一区二区三区| 2018国产大陆天天弄谢| 丁香六月天网| 色婷婷av一区二区三区视频| 丰满饥渴人妻一区二区三| 亚洲欧美一区二区三区黑人 | 亚洲视频免费观看视频| 2018国产大陆天天弄谢| 免费少妇av软件| 2022亚洲国产成人精品| 欧美 亚洲 国产 日韩一| 国产一区二区三区综合在线观看| 青青草视频在线视频观看| 一区福利在线观看| 成年女人毛片免费观看观看9 | 99九九在线精品视频| 国产探花极品一区二区| 国产综合精华液| 男人添女人高潮全过程视频| 日日撸夜夜添| 亚洲国产精品国产精品| 亚洲欧美成人综合另类久久久| 韩国av在线不卡| 国产成人精品久久久久久| 观看美女的网站| 午夜影院在线不卡| 亚洲欧洲国产日韩| 又黄又粗又硬又大视频| 国产无遮挡羞羞视频在线观看| 又粗又硬又长又爽又黄的视频| 看免费av毛片| 一本色道久久久久久精品综合| 激情五月婷婷亚洲| 性色av一级| 男女无遮挡免费网站观看| 只有这里有精品99| 中文字幕最新亚洲高清| 水蜜桃什么品种好| 1024香蕉在线观看| 国产精品 国内视频| 亚洲精品国产色婷婷电影| 97精品久久久久久久久久精品| 18禁国产床啪视频网站| 两性夫妻黄色片| 啦啦啦中文免费视频观看日本| 亚洲av在线观看美女高潮| 在线天堂最新版资源| 亚洲 欧美一区二区三区| 亚洲,欧美精品.| 2018国产大陆天天弄谢| 久久久a久久爽久久v久久| 亚洲国产av新网站| 亚洲伊人色综图| 精品国产乱码久久久久久小说| 麻豆精品久久久久久蜜桃| 精品人妻偷拍中文字幕| 亚洲欧美精品综合一区二区三区 | 日韩av免费高清视频| 国产一区有黄有色的免费视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产成人一精品久久久| 欧美人与性动交α欧美精品济南到 | 精品视频人人做人人爽| av又黄又爽大尺度在线免费看| 乱人伦中国视频| 久久精品久久久久久噜噜老黄| 久久人人爽人人片av| 大话2 男鬼变身卡| 亚洲av免费高清在线观看| 欧美最新免费一区二区三区| 午夜老司机福利剧场| 人人妻人人澡人人爽人人夜夜| 日日摸夜夜添夜夜爱| 国产一区二区三区综合在线观看| 日本av手机在线免费观看| 亚洲男人天堂网一区| av天堂久久9| 午夜福利一区二区在线看| 亚洲内射少妇av| 少妇精品久久久久久久| 在线观看三级黄色| 一区福利在线观看| av网站在线播放免费| 一个人免费看片子| 各种免费的搞黄视频| 国产野战对白在线观看| 爱豆传媒免费全集在线观看| 欧美av亚洲av综合av国产av | 免费在线观看视频国产中文字幕亚洲 | 成人18禁高潮啪啪吃奶动态图| 国产又爽黄色视频| 制服诱惑二区| 我的亚洲天堂| 国产国语露脸激情在线看| 免费av中文字幕在线| 麻豆精品久久久久久蜜桃| 精品少妇黑人巨大在线播放| 少妇人妻 视频| 国产精品熟女久久久久浪| 少妇猛男粗大的猛烈进出视频| 午夜激情av网站| 美女中出高潮动态图| 在线观看免费日韩欧美大片| 人人妻人人添人人爽欧美一区卜| 国产黄色视频一区二区在线观看| 亚洲精品久久成人aⅴ小说| 国产极品粉嫩免费观看在线| 久久 成人 亚洲| 午夜影院在线不卡| 一级毛片 在线播放| 国产在线一区二区三区精| 一区在线观看完整版| 欧美日韩亚洲高清精品| 精品一区二区三区四区五区乱码 | 亚洲国产欧美在线一区| 日韩制服丝袜自拍偷拍| 寂寞人妻少妇视频99o| 久久鲁丝午夜福利片| 狠狠婷婷综合久久久久久88av| 亚洲av福利一区| 久久久久久人人人人人| 少妇被粗大猛烈的视频| 老熟女久久久| 亚洲精品乱久久久久久| 精品一品国产午夜福利视频| 国产精品一区二区在线不卡| 老汉色av国产亚洲站长工具| 国产精品久久久久久av不卡| 99精国产麻豆久久婷婷| 我要看黄色一级片免费的| 一级毛片黄色毛片免费观看视频| 亚洲一区中文字幕在线| av.在线天堂| 国产精品久久久久成人av| 免费观看在线日韩| 成年人午夜在线观看视频| 伊人亚洲综合成人网| 亚洲精品第二区| 国产精品偷伦视频观看了| 久久精品久久久久久久性| 哪个播放器可以免费观看大片| 亚洲精品中文字幕在线视频| 国产欧美日韩综合在线一区二区| 久久精品国产鲁丝片午夜精品| 最近最新中文字幕免费大全7| 在线观看一区二区三区激情| 久久久精品94久久精品| 久久av网站| 看免费av毛片| 国产av精品麻豆| 国产在视频线精品| 日韩人妻精品一区2区三区| 最近2019中文字幕mv第一页| 在线观看免费日韩欧美大片| 精品99又大又爽又粗少妇毛片| 少妇人妻久久综合中文| 亚洲欧美色中文字幕在线| 有码 亚洲区| av视频免费观看在线观看| 两个人看的免费小视频| 免费大片黄手机在线观看| 嫩草影院入口| 国产福利在线免费观看视频| 青春草国产在线视频| 婷婷色综合大香蕉| 精品亚洲成国产av| 2018国产大陆天天弄谢| 欧美黄色片欧美黄色片| 精品少妇内射三级| 国产成人精品久久二区二区91 | 欧美日韩av久久| 亚洲激情五月婷婷啪啪| 久久久久久免费高清国产稀缺| 看非洲黑人一级黄片| 新久久久久国产一级毛片| 日韩中文字幕视频在线看片| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 女性生殖器流出的白浆| 亚洲av免费高清在线观看| 日韩av在线免费看完整版不卡| 另类精品久久| 久久这里只有精品19| 午夜影院在线不卡| 久久影院123| 国产一区二区三区综合在线观看| 一区二区三区乱码不卡18| 大话2 男鬼变身卡| 亚洲欧美日韩另类电影网站| 日本色播在线视频| 日韩制服丝袜自拍偷拍| 大片电影免费在线观看免费| 深夜精品福利| 日韩欧美精品免费久久| 男女啪啪激烈高潮av片| 美女视频免费永久观看网站| 亚洲一级一片aⅴ在线观看| 成年动漫av网址| 黄片播放在线免费| 最近中文字幕2019免费版| 亚洲欧美中文字幕日韩二区| 一边摸一边做爽爽视频免费| 熟女电影av网| 亚洲国产色片| 国产精品成人在线| 熟妇人妻不卡中文字幕| 久久久久视频综合| 熟女av电影| 日本av手机在线免费观看| 黑丝袜美女国产一区| 婷婷色综合www| 黑丝袜美女国产一区| 亚洲精品日本国产第一区| 一区二区三区精品91| 国产精品久久久久久久久免| 18禁裸乳无遮挡动漫免费视频| 在线观看免费视频网站a站| 色视频在线一区二区三区| 成人亚洲精品一区在线观看| 亚洲四区av| 麻豆av在线久日| 丰满饥渴人妻一区二区三| 精品一区二区三区四区五区乱码 | 欧美国产精品va在线观看不卡| 男的添女的下面高潮视频| 老汉色∧v一级毛片| 国产97色在线日韩免费| 秋霞伦理黄片| 国产综合精华液| 人人妻人人澡人人爽人人夜夜| 亚洲一区中文字幕在线| 免费观看a级毛片全部| 伊人久久国产一区二区| 亚洲欧洲日产国产| 日本黄色日本黄色录像| 亚洲欧美清纯卡通| 久久久久久久久久久久大奶| 亚洲国产欧美日韩在线播放| av片东京热男人的天堂| 亚洲精华国产精华液的使用体验| 午夜福利视频精品| 国产成人精品无人区| 亚洲五月色婷婷综合| 亚洲欧美成人综合另类久久久| 日本欧美视频一区| 久久99热这里只频精品6学生| 中文字幕人妻丝袜制服| 麻豆乱淫一区二区| 日韩制服骚丝袜av| 国产精品国产av在线观看| 高清不卡的av网站| 亚洲av日韩在线播放| 激情视频va一区二区三区| 久久久精品国产亚洲av高清涩受| 久久婷婷青草| 一区二区三区四区激情视频| 国产成人精品福利久久| 超碰成人久久| 两性夫妻黄色片| tube8黄色片| 久久综合国产亚洲精品| 最近中文字幕2019免费版| 欧美成人午夜精品| av片东京热男人的天堂| 尾随美女入室| 人妻人人澡人人爽人人| 久久99蜜桃精品久久| 免费黄频网站在线观看国产| 青青草视频在线视频观看| av.在线天堂| 国产无遮挡羞羞视频在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲男人天堂网一区| 精品久久久精品久久久| 成人手机av| 涩涩av久久男人的天堂| 黑人欧美特级aaaaaa片| 午夜福利在线观看免费完整高清在| 建设人人有责人人尽责人人享有的| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站在线观看播放| 国产成人a∨麻豆精品| 91午夜精品亚洲一区二区三区| 国产爽快片一区二区三区| 国产精品 欧美亚洲| 90打野战视频偷拍视频| 久久精品夜色国产| 亚洲综合色惰| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 少妇的逼水好多| 91精品伊人久久大香线蕉| 1024视频免费在线观看| 丁香六月天网| 天美传媒精品一区二区| 人成视频在线观看免费观看| 欧美成人午夜免费资源| 国产精品一区二区在线不卡| 免费观看在线日韩| 久热久热在线精品观看| 一级片'在线观看视频| 亚洲中文av在线| 欧美日韩精品网址| 日韩精品有码人妻一区| 久久精品aⅴ一区二区三区四区 | 国产在线视频一区二区| 黄色视频在线播放观看不卡| 人人妻人人澡人人看| 国产成人午夜福利电影在线观看| 精品第一国产精品| 寂寞人妻少妇视频99o| 三上悠亚av全集在线观看| 免费av中文字幕在线| 国产av精品麻豆| 婷婷色av中文字幕| 国产激情久久老熟女| 夫妻性生交免费视频一级片| 最近的中文字幕免费完整| 亚洲av男天堂| 不卡视频在线观看欧美| 国产成人aa在线观看| 久久精品久久久久久噜噜老黄| 久久精品aⅴ一区二区三区四区 | 黄色一级大片看看| 免费黄色在线免费观看| 爱豆传媒免费全集在线观看| 欧美激情极品国产一区二区三区| 97精品久久久久久久久久精品| 成人国语在线视频| 伦精品一区二区三区| 一级毛片我不卡| 最近最新中文字幕大全免费视频 | 黄片无遮挡物在线观看| 国产欧美日韩一区二区三区在线| 99久久中文字幕三级久久日本| 久久久久精品人妻al黑| 亚洲欧美精品综合一区二区三区 | 熟妇人妻不卡中文字幕| 久久影院123| 高清黄色对白视频在线免费看| 黑人巨大精品欧美一区二区蜜桃| 久久ye,这里只有精品| 在线天堂最新版资源| 欧美亚洲日本最大视频资源| 免费av中文字幕在线| 中文字幕人妻丝袜制服| 精品99又大又爽又粗少妇毛片| 成人18禁高潮啪啪吃奶动态图| 亚洲一码二码三码区别大吗| 免费黄网站久久成人精品| 美女xxoo啪啪120秒动态图| 永久免费av网站大全| 国产日韩欧美亚洲二区| 午夜免费鲁丝| 亚洲精品视频女|