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

    水下滑翔器整體外形設(shè)計(jì)及水動(dòng)力性能分析

    2014-10-11 06:20:20孟凡豪嚴(yán)天宏徐新勝
    海洋工程 2014年2期
    關(guān)鍵詞:后掠角滑翔攻角

    孟凡豪,嚴(yán)天宏,何 波,徐新勝

    (1.中國(guó)計(jì)量學(xué)院機(jī)電工程學(xué)院,浙江杭州 310018;2.中國(guó)海洋大學(xué)信息科學(xué)與工程學(xué)院,山東青島025506)

    水下滑翔器(AUG,Autonomous Underwater Glider)是20世紀(jì)90年代以浮標(biāo)技術(shù)為基礎(chǔ)發(fā)展起來(lái)的新型水下測(cè)量?jī)x器搭載平臺(tái)[1]。與傳統(tǒng)的水下自航行平臺(tái)相比較,其具有航程大,低噪聲,低成本等優(yōu)勢(shì)。近年來(lái),隨著海洋監(jiān)測(cè)和海洋開(kāi)發(fā)的需求增長(zhǎng),水下滑翔器的研究與應(yīng)用前景得到了海洋研究領(lǐng)域的充分認(rèn)識(shí),并廣泛地開(kāi)展起來(lái)[2-4]。目前可以成功進(jìn)行實(shí)際運(yùn)用的水下滑翔器主要有三大類(lèi)型:Scripps海洋研究所研制的Spray滑翔器,華盛頓大學(xué)設(shè)計(jì)的Seaglider滑翔器,Webb公司生產(chǎn)的Slocum滑翔器。這3種滑翔器的航程都超過(guò)了2 000 km,水下續(xù)航時(shí)間達(dá)到百余天,甚至1年[5-9]。

    現(xiàn)在水下滑翔器的驅(qū)動(dòng)方式有依靠浮力驅(qū)動(dòng)、電能驅(qū)動(dòng)、溫差驅(qū)動(dòng),還出現(xiàn)利用太陽(yáng)能驅(qū)動(dòng)的水下滑翔器[10-11]。根據(jù)水下滑翔器工作環(huán)境的不同采用不同的驅(qū)動(dòng)方式,這些都是為了提高水下滑翔器的水下速度、延長(zhǎng)續(xù)航能力等特性。通過(guò)查閱資料了解到,現(xiàn)在的水下滑翔器很難同時(shí)滿足高速、長(zhǎng)續(xù)航時(shí)間這兩個(gè)主要性能指標(biāo)。水下滑翔器的經(jīng)濟(jì)性和穩(wěn)定性是目前研發(fā)人員最關(guān)注的問(wèn)題之一。減小水下滑翔器的俯仰力矩,提高升阻比,是現(xiàn)階段提高水下滑翔器經(jīng)濟(jì)性和穩(wěn)定性的主要手段。

    基于CFD方法,對(duì)水下滑翔器的主體進(jìn)行了設(shè)計(jì)仿真,通過(guò)不同構(gòu)型主體的水動(dòng)力仿真對(duì)比,獲得最優(yōu)外形;同時(shí),還對(duì)水下滑翔器的附體機(jī)翼和尾翼進(jìn)行了翼型的選擇及設(shè)計(jì),進(jìn)一步提高了水下滑翔器的水動(dòng)力特性;對(duì)設(shè)計(jì)后的水下滑翔器進(jìn)行了整體的水動(dòng)力分析,并與Slocum等幾種經(jīng)典水下滑翔器樣機(jī)模型進(jìn)行了性能比較,獲得在特定攻角范圍內(nèi)的最佳升阻比,同時(shí)驗(yàn)證了新的水下滑翔器模型的水動(dòng)力性能有了很大的提升。

    1 主體設(shè)計(jì)及選擇

    1.1 主體模型的建立

    水下滑翔器的基本功能相差不大,已有水下滑翔器的總體質(zhì)量大約為50~70 kg左右。設(shè)計(jì)的目標(biāo)水下滑翔器總質(zhì)量定在60 kg左右。

    根據(jù)水下滑翔器的質(zhì)量M,初步估算滑翔器在水中的排水量:

    滑翔器的主體外形大體定為長(zhǎng)2 m,主體直徑0.2 m。下面有五種主體外形設(shè)計(jì)構(gòu)型,在選擇時(shí)要充分考慮主體的水動(dòng)力性能、可加工性以及內(nèi)部空間的可利用性。如圖1~圖3所示,構(gòu)型Ⅰ到構(gòu)型Ⅲ是基于Slocum水下滑翔器樣機(jī)在已有的尺寸上進(jìn)行修改得到的[12-13]。

    圖1 AUG1主體構(gòu)型IFig.1 Body Type I

    圖2 AUG2主體構(gòu)型IIFig.2 Body Type II

    圖3 AUG3主體構(gòu)型IIIFig.3 Body Type III

    構(gòu)型Ⅰ:主體兩端采用直徑為0.2 m的半球體,中間部分是直徑為0.2 m,長(zhǎng)為1.8 m的圓柱體,如圖1所示。該構(gòu)型的主要特點(diǎn)是結(jié)構(gòu)加工簡(jiǎn)單,成本低,空間有效利用率較高,但是水動(dòng)力特性很差。

    構(gòu)型Ⅱ:主體的首端依然是直徑為0.2 m的半球體,尾端修正為偏心率e為3.5的半橢圓體,中間部分為直徑0.2 m的圓柱體,如圖2所示。該構(gòu)型的主要特點(diǎn)是相對(duì)構(gòu)型1其水動(dòng)力特性得到了一定的提高,但是加工復(fù)雜性增加,內(nèi)部空間有效利用率下降。

    構(gòu)型Ⅲ:主體首端修正為偏心率e為3的半橢圓體,尾端為偏心率e為3.5的半橢圓體,中間部分為直徑0.2 m的圓柱體,如圖3所示。該構(gòu)型進(jìn)一步提升了水下滑翔器的水動(dòng)力特性,但在加工、有效空間利用等方面的性能被進(jìn)一步的削弱。前3種構(gòu)型是水下滑翔器現(xiàn)在應(yīng)用比較廣泛的幾種主流主體外形,綜合其特點(diǎn)為加工簡(jiǎn)單、成本低廉、空間利用率較高,有較好的水動(dòng)力性能。

    構(gòu)型Ⅳ:依據(jù)標(biāo)準(zhǔn)橢圓線性方程設(shè)計(jì)的頭尾對(duì)稱的雙參數(shù)橢圓回轉(zhuǎn)體作為水下滑翔器的主體,其長(zhǎng)短軸分為1.8和0.24 m,如圖4所示。該構(gòu)型在滑翔器的水動(dòng)力性能和有效空間利用上都有很大的提升空間。

    構(gòu)型Ⅴ:采用NACA0012流線型線性方程設(shè)計(jì)的“雨滴”狀回轉(zhuǎn)體,其長(zhǎng)度為2 m,最大直徑為0.24 m,如圖5所示。該構(gòu)型與前三種構(gòu)型相比較,其外形可以推遲邊界層分離、減少漩渦阻力,使水下滑翔器的水動(dòng)力性能有很大的提高,相對(duì)不足之處是有效空間利用有所下降。用“雨滴”狀回轉(zhuǎn)體作為水下滑翔器的主體是一個(gè)積極而有益的探索。

    圖4 AUG4主體構(gòu)型IVFig.4 Body Type IV

    圖5 AUG5主體構(gòu)型VFig.5 Body Type V

    下文采用計(jì)算流體力學(xué)軟件Fluent?分別對(duì)上述五種主體構(gòu)型進(jìn)行水動(dòng)力計(jì)算與分析,通過(guò)綜合對(duì)比分析得到最優(yōu)結(jié)果。

    1.2 網(wǎng)格的劃分及計(jì)算條件的設(shè)定

    滑翔器在水下航行時(shí),以雷諾數(shù)來(lái)判別流動(dòng)狀態(tài),雷諾數(shù)的定義公式為

    式中:ρ為流體密度,v為入流速度,L為流體特征長(zhǎng)度,μ為流體粘性系數(shù)。對(duì)于外流場(chǎng)判斷湍流的標(biāo)準(zhǔn)為[14]:外掠平板時(shí)雷諾數(shù)Re大于5×105,外掠障礙物時(shí)大于2×104。這五種構(gòu)型中,雷諾數(shù)Re在1.6×106左右,故為湍流。

    網(wǎng)格劃分采用非結(jié)構(gòu)化網(wǎng)格,數(shù)量在100萬(wàn)到110萬(wàn)。在靠近水下滑翔器殼體的壁面處劃分出一個(gè)小的區(qū)域并進(jìn)行網(wǎng)格加密,這樣可以保證網(wǎng)格疏密合理分布,又能減少總體網(wǎng)格數(shù)目,在不降低準(zhǔn)確度的條件下可以加快計(jì)算速度。

    計(jì)算條件的設(shè)定

    1)計(jì)算區(qū)域?yàn)殚L(zhǎng)12 m,半徑3 m的圓柱體區(qū)域,邊界條件為速度入口,壓力出口。計(jì)算區(qū)域是模型的六倍,可以充分模擬流場(chǎng)環(huán)境,提高數(shù)值模擬的準(zhǔn)確性,同時(shí)采用壓力出口條件,這樣可以考慮回流情況,使計(jì)算更容易收斂。

    2)水下滑翔器的所有壁面均設(shè)定為靜止無(wú)滑移壁面,壁面粗糙度為0,其它壁面均為無(wú)滑移壁面。流體材料設(shè)定為水密度998.2 kg/m3,運(yùn)動(dòng)粘性0.001 003 kg/m·s。

    3)采用三維單精度壓力基求解器,壓力基求解器主要用于低速不可壓縮流動(dòng)的求解,密度求解器收斂速度快,但需要更多的內(nèi)存和計(jì)算量,這里考慮到計(jì)算的是一個(gè)低速不可壓縮流場(chǎng),網(wǎng)格的數(shù)量很大,故采用壓力基求解器。湍流模型選取為RNG k-ε模型,RNG k-ε模型考慮了湍流漩渦,RNG理論提供了一個(gè)考慮低雷諾數(shù)流動(dòng)粘性的解析公式,這些特點(diǎn)使得RNG k-ε模型比標(biāo)準(zhǔn)k-ε模型等在更廣泛的流動(dòng)中有更高的可信度和精度[4]。

    4)壓力-速度耦合采用SIMPLE方法;離散格式中壓力采用PRESTO離散,動(dòng)量方程中的對(duì)流項(xiàng)采用二階迎風(fēng)差分格式,其余為一階迎風(fēng)差分格式。PRESTO離散可以保證計(jì)算的穩(wěn)定性,二階迎風(fēng)差分格式可以提高計(jì)算精度,但是會(huì)增加計(jì)算量。所有參數(shù)的收斂殘差標(biāo)準(zhǔn)為0.000 1。

    1.3 模型數(shù)值計(jì)算準(zhǔn)確性的驗(yàn)證

    為了驗(yàn)證數(shù)值計(jì)算的合理性和準(zhǔn)確性,按照上述方法計(jì)算的Slocum水下滑翔器在不同攻角下的阻力系數(shù)和文獻(xiàn)[15-16]中的試驗(yàn)數(shù)據(jù),對(duì)比如表1所示。由表1可知,數(shù)值計(jì)算結(jié)果和試驗(yàn)結(jié)果相近,文獻(xiàn)[16]與試驗(yàn)結(jié)果最大誤差為9.35%,本文最大誤差為7.2%,完全滿足計(jì)算的精度,驗(yàn)證了本文CFD模型的正確性與可靠性。

    表1 針對(duì)阻力系數(shù)[15]CD本文數(shù)值模擬結(jié)果對(duì)比驗(yàn)證Tab.1 The comparison of the CDbetween the proposed numerical simulation and Reference[15]

    1.4 模型數(shù)值計(jì)算收斂性的驗(yàn)證

    為了驗(yàn)證數(shù)值計(jì)算的收斂性,計(jì)算模型2在速度為0.5 kn、攻角0°下的水動(dòng)力值。所有參數(shù)的收斂殘差標(biāo)準(zhǔn)為0.000 1。改變模型2的網(wǎng)格疏密程度,圖6和圖7分別是在80萬(wàn)網(wǎng)格和120萬(wàn)網(wǎng)格下,模型2的升力系數(shù)收斂曲線。從圖6到圖7可以看出,網(wǎng)格疏密不影響計(jì)算結(jié)果,迭代1 000次后升力系數(shù)不再變化,說(shuō)明數(shù)值計(jì)算已收斂。

    圖6 80萬(wàn)網(wǎng)格時(shí)的模型2升力系數(shù)CL收斂曲線Fig.6 Lift CLconvergence curve for Model 2

    圖7 120萬(wàn)網(wǎng)格時(shí)模型2的升力系數(shù)CL收斂曲線Fig.7 Lift CLconvergence curve for Model 2

    1.5 水動(dòng)力性能的計(jì)算與分析

    利用Fluent?的后處理功能可以計(jì)算出主體受到的阻力,無(wú)量綱化得到阻力系數(shù)。阻力系數(shù)和升力系數(shù)的計(jì)算公式分別為:

    式中:ρ為水的密度;V為入流速度;Sd為主體的迎流截面面積或表面面積(選用不同面積時(shí),阻力系數(shù)不同,這里選迎流截面面積計(jì)算);當(dāng)計(jì)算的是機(jī)翼時(shí),Sl為翼面的最大面積;L為升力;D為阻力。表2是五種方案主體在0.5 kn速度下的主要水動(dòng)力參數(shù)。

    表2 五種主體主要水動(dòng)力參數(shù)Tab.2 The main hydraulic parameters of five body types

    表2的計(jì)算結(jié)果表明:在速度為0.5 kn,攻角為0°的情況下,從構(gòu)型AUG1到AUG5水下滑翔器主體構(gòu)型的水動(dòng)力性能不斷提高,AUG4橢圓回轉(zhuǎn)體與前三種主體構(gòu)型相比,阻力分別減少了27%、10%、4%;AUG5“雨滴”狀回轉(zhuǎn)體與其他四種主體構(gòu)型相比,阻力分別減少了33%、17.8%、12.2%和8.5%;這兩種構(gòu)型使水下滑翔器的阻力都有很大幅度的降低。圖8是這五種主體構(gòu)型在攻角為0°、不同速度條件下繪制的阻力曲線對(duì)比圖。圖8表明:隨著水下滑翔器速度的增大,主體所受到的阻力也不斷增大。通過(guò)對(duì)比可以看出AUG4、AUG5所受阻力的增長(zhǎng)幅度始終低于其余三種主體,并且隨著速度的增加,AUG4、AUG5的阻力增長(zhǎng)幅度越來(lái)越小于其余三種主體的阻力增長(zhǎng)幅度。AUG5“雨滴”狀回轉(zhuǎn)體的尾端比較狹長(zhǎng),使有效空間的利用率很低,加工難度較大;AUG4水動(dòng)力性能良好,空間利用率較高,加工難度不大,綜合考慮,這里選取AUG4作為水下滑翔器的主體。從圖9、圖10,可以分別看出AUG4主體構(gòu)型的動(dòng)壓力分布和速度分布情況。

    圖8 攻角為0°不同速度下主體阻力曲線對(duì)比Fig.8 The drag force vs.the velocity with attack angle at 0°

    圖9 AUG4的動(dòng)壓(Dynamic Pressure)分布Fig.9 The distribution of dynamic pressure of AUG Type IV

    圖10 AUG4的二維速度矢量圖Fig.10 The 2D velocity vector of AUG Body Type IV

    2 附體設(shè)計(jì)及優(yōu)選

    水下滑翔器的附體包括機(jī)翼和方向尾翼。水下滑翔器要求滑翔機(jī)翼及方向尾翼在水下滑翔機(jī)的攻角、漂角范圍內(nèi),盡量避免在翼上發(fā)生流動(dòng)分離現(xiàn)象,以減少航行阻力和能耗。同時(shí)機(jī)翼應(yīng)具備較大的升阻比,以提高兩翼的工作效率。

    2.1 機(jī)翼設(shè)計(jì)

    2.1.1 翼型選擇

    因?yàn)樗禄铏C(jī)要在水中完成上浮和下潛運(yùn)動(dòng),攻角會(huì)出現(xiàn)正負(fù)交替的現(xiàn)象,要保證在兩種狀態(tài)下翼型均有良好表現(xiàn),選用對(duì)稱翼型。低速的對(duì)稱翼型可選用平板翼和曲面翼兩種。

    相關(guān)資料表明[17-18]:平板翼型的特點(diǎn)是型面簡(jiǎn)單、易于加工、成本低廉但相應(yīng)的流體動(dòng)力性能比曲面翼型差得多,易發(fā)生流動(dòng)分離使阻力增大。而曲面翼型則相反,不易發(fā)生流動(dòng)分離現(xiàn)象。因此,翼型的選擇要在效率和加工工藝之間加以權(quán)衡從而決定采用哪種。為了研究?jī)煞N翼型對(duì)升阻比的影響,采用Fluent?軟件對(duì)面積相同的平板翼和曲面翼進(jìn)行水動(dòng)力模擬。圖11是模擬水流速度為0.3 m/s,不同攻角下的兩組翼型的升阻比變化曲線,顯而易見(jiàn),曲面翼的升阻比性能優(yōu)于平板翼。通過(guò)圖11對(duì)比可知,兩種翼型的升阻比都會(huì)隨著攻角的增大而提高,曲面翼的升阻比增大趨勢(shì)明顯大于平板翼,可以提高50%左右,其水動(dòng)力性能更加優(yōu)越,雖然加工相對(duì)復(fù)雜,但曲面翼仍會(huì)成為水下滑翔器的主流翼型。

    2.1.2 機(jī)翼優(yōu)選

    目前水下滑翔器選擇的曲面翼型基本為NACA系列翼型,本文將采用NACA0005翼型。在主體參數(shù)確定的情況下改變NACA0005機(jī)翼的翼展長(zhǎng)l,梢弦長(zhǎng)b1,后掠角λ,如圖12所示。

    圖11 兩種機(jī)翼升阻比隨攻角變化關(guān)系曲線圖Fig.11 The lift and drag ratio curves of 2 types of airfoils vs.the attack angle

    圖12 翼的平面形狀[19]Fig.12 The planar shape of airfoil[19]

    在攻角為6°,航速為0.3 m/s,根梢比η=2.5的情況下進(jìn)行3因素3水平正交試驗(yàn)。正交試驗(yàn)法[16]可以得到設(shè)計(jì)參數(shù)的多種組合,使實(shí)驗(yàn)結(jié)果更加豐富,可以分析相應(yīng)規(guī)律得到最優(yōu)結(jié)果。正交參數(shù)仿真組合情況如表3所示。

    表3 正交參數(shù)仿真組合表Tab.3 The orthogonal parameters combination for simulation

    機(jī)翼對(duì)水下滑翔器經(jīng)濟(jì)性的影響,可以由滑翔器的升阻比L/D來(lái)評(píng)價(jià)[16],升阻比越大,滑翔經(jīng)濟(jì)性越高。正交試驗(yàn)的組合模型,采用CFD軟件進(jìn)行水動(dòng)力模擬仿真實(shí)驗(yàn),得到每個(gè)模型的升阻比。正交試驗(yàn)結(jié)果如表4所示。

    表4 正交參數(shù)表及仿真試驗(yàn)結(jié)果Tab.4 The orthogonal parameters and corresponding simulation results

    正交試驗(yàn)結(jié)果表明,第4,7,8,9組翼型的實(shí)驗(yàn)結(jié)果與其他幾組相比水動(dòng)力性能提升很高。圖13與圖14是第9組機(jī)翼表面速度分布和壓力分布情況??煽闯龅?組實(shí)驗(yàn)翼型能夠有效避免翼上發(fā)生流動(dòng)分離現(xiàn)象,以減少航行阻力、降低能耗;同時(shí)具有較大的升阻比,可提高兩翼的工作效率。

    圖13 機(jī)翼表面速度分布Fig.13 The velocity distribution for airfoil

    圖14 機(jī)翼表面動(dòng)壓力分布Fig.14 The dynamic pressure distribution for airfoil

    2.2 尾翼設(shè)計(jì)

    尾翼的主要作用是保證水下滑翔器在最大速度及帶攻角運(yùn)動(dòng)時(shí)具有足夠的效率,通過(guò)與機(jī)翼配合產(chǎn)生的流體動(dòng)力來(lái)提供水下滑翔器航向的穩(wěn)定性。

    根據(jù)文獻(xiàn)[19],尾翼選用NACA0024翼型,可以保證足夠的效率和穩(wěn)定性。尾翼的阻力效率受展弦比及后掠角的影響,選用較大展弦比及小后掠角盡量減小尾翼的阻力。

    表5 尾翼后掠角性能比較Tab.5 The comparison of performances of different angles

    翼展l=0.2 m,展弦比λ=2,根梢比η=1.6,根弦長(zhǎng) b0=0.16 m,梢弦長(zhǎng) b1=0.1 m。設(shè)定如表5中設(shè)計(jì)4個(gè)尾翼的后掠角方案,攻角為 0°,速度 0.3 m/s。利用Matlabc的曲線擬合工具,將阻力系數(shù)和后掠角的關(guān)系擬合為二次函數(shù),并得到兩者的經(jīng)驗(yàn)公式,如圖15所示。

    由經(jīng)驗(yàn)公式和圖形以及表5得出結(jié)論:尾翼的后掠角在15°到20°之間阻力系數(shù)最小,考慮到加工和材料的因素,選擇15°和20°作為較優(yōu)值。

    2.3 具體模型分析

    圖15 阻力系數(shù)和后掠角的關(guān)系擬合曲線Fig.15 Relationship between CDand sweep angle

    根據(jù)飛機(jī)的外形設(shè)計(jì)經(jīng)驗(yàn)可知,從減小阻力的角度考慮,機(jī)翼最大橫截面位置在機(jī)身長(zhǎng)度的60%處最好。水下滑翔器機(jī)翼的最大橫截面布置在主載體長(zhǎng)度60%處,可獲得最好的水動(dòng)力性能。根據(jù)上述研究和分析結(jié)果,為進(jìn)一步研究在不同航行攻角下附體對(duì)水下滑翔器的升阻比L/D和最大俯仰力矩 Mpitch的影響,對(duì)表5所示的5個(gè)具有不同機(jī)翼參數(shù)的具體模型進(jìn)行航速為0.5 kn、攻角為0°~20°的數(shù)值模擬。升阻比L/D和最大俯仰力矩 Mpitch決定水下滑翔器的經(jīng)濟(jì)性和穩(wěn)定性[16]。依據(jù)正交試驗(yàn)結(jié)果和尾翼后掠角性能比較結(jié)果,本文設(shè)計(jì)了5種水下滑翔器模型,這5種模型是由正交試驗(yàn)第4,7,8,9組實(shí)驗(yàn)翼型和尾翼模型中的方案2和3組合而成,具體模型參數(shù)如表5所示。圖16和圖17分別是模型2和模型3的外形效果圖。圖18為8°攻角下模型2的表面速度矢量分布圖,圖19為8°攻角下模型3的表面動(dòng)壓力分布云圖。各模型其他參數(shù)的模擬結(jié)果比較如圖20和圖21所示。

    圖16 模型2的外形效果圖Fig.16 The shape of Model 2

    圖17 模型3的外形效果圖Fig.17 The shape of Model 3

    圖18 模型2的表面速度矢量分布圖Fig.18 The velocity distribution for Model 2

    圖19 模型3的表面動(dòng)壓力分布云圖Fig.19 The dynamic pressure distribution for Model 3

    表6 具體模型參數(shù)Tab.6 Parameters of specific models

    由表6和圖19中的模型1、模型2、模型3和模型5可知,水下滑翔器的升阻比L/D隨著翼展的增加而增加,而隨著后掠角的增大而下降;在上面機(jī)翼的正交試驗(yàn)中,機(jī)翼的升阻比L/D隨著后掠角的增大而增大,由此可推知,單獨(dú)進(jìn)行水動(dòng)力性能分析最優(yōu)的機(jī)翼不能確保水下滑翔器整體的水動(dòng)力性能達(dá)到最優(yōu)。由圖21中的模型1、模型2和模型5可知,水下滑翔器的最大俯仰力矩 Mpitch隨著后掠角的增加而減小,即穩(wěn)定性增加,機(jī)翼的重心位置越靠后其穩(wěn)定性越高。由圖20和圖21中的模型4和模型5可知,在低攻角范圍內(nèi),水下滑翔器的升阻比L/D隨著尾翼后掠角的增大而減小,當(dāng)攻角超過(guò)6°后,隨著尾翼后掠角的增大而略微增大;水下滑翔器的最大俯仰力矩 Mpitch隨著尾翼后掠角的增大而增大,即穩(wěn)定性下降。

    由圖20可知,模型3的升阻比明顯大于其余模型,模型1和模型2的升阻比值相近,模型2整體上略微優(yōu)于模型1。模型3的升阻比與模型2、模型4和模型5相比提高了4.4%、9%和10%,模型2與模型4和模型5相比提高了4.5%和4.8%。而由圖21可以看出,模型3的最大俯仰力矩大于其余模型,其穩(wěn)定性下降明顯。模型3的穩(wěn)定性與模型2、模型4和模型5相比下降了19.8%、23%和24%。模型2與模型1最大俯仰力矩曲線基本相近,比模型4和模型5增大了4.3%和4.7%。

    圖20 升阻比L/D隨攻角的變化關(guān)系Fig.20 Relationship between L/D and attack angle

    圖21 最大俯仰力矩 Mpitch隨攻角的變化關(guān)系Fig.21 Relationship between Mpitch and attack angle

    由上述分析可知,在滑翔狀態(tài)下,為提高滑翔經(jīng)濟(jì)性應(yīng)盡量增加機(jī)翼的翼展,同時(shí)減小機(jī)翼的后掠角,但為提高航行器的穩(wěn)定性,又需要增大機(jī)翼的后掠角來(lái)實(shí)現(xiàn)??梢?jiàn),提高滑翔器滑翔經(jīng)濟(jì)性和增加穩(wěn)定性這兩方面是矛盾的,設(shè)計(jì)需要根據(jù)其他性能要求來(lái)選擇。

    3 具體模型與樣機(jī)水動(dòng)力性能對(duì)比分析

    為進(jìn)一步選取最優(yōu)模型,將模型2和模型3與Spray、Seaglider、Slocum這三種經(jīng)典水下滑翔機(jī)相比較,它們機(jī)體尺寸和重量相差不多,體長(zhǎng)2~3 m左右,重量約為50~70 kg,根據(jù)Spray、Seaglider、Slocum水下滑翔機(jī)的外形設(shè)計(jì)參數(shù)[20-21],可得到這三種樣機(jī)模型,如圖22~圖24所示。這里通過(guò)CFD模擬仿真,獲得三種經(jīng)典樣機(jī)的水動(dòng)力性能參數(shù),其數(shù)據(jù)可靠性已經(jīng)在上文中驗(yàn)證。模型2和模型3與這三種樣機(jī)在航速0.5 kn,不同的攻角下進(jìn)行水動(dòng)力性能分析對(duì)比,具體仿真參數(shù)和指標(biāo)見(jiàn)表7。

    圖22 模型 I——Spray樣機(jī)Fig.22 The classic type I-Spray Model

    圖23 模型II——Seaglider樣機(jī)Fig.23 The classic type II-Seaglider Model

    圖24 模型III——Slocum樣機(jī)Fig.24 The classic model III-Slocum Model

    表7 模型2和模型3與樣機(jī)主要參數(shù)對(duì)比Tab.7 The comparison of main parameters between the specific models and the classic types

    圖25 航速為0.5 kn不同攻角下具體模型和樣機(jī)升阻比曲線對(duì)照Fig.25 The L/D ratio comparison between the specific models and the classic type under different attack angles

    通過(guò)表7可知,模型2、模型3和三種主流樣機(jī)相比較,排水量相對(duì)增加明顯,增加了機(jī)體的內(nèi)部空間;最大升阻比提升很大,使水下滑翔器具有更好的經(jīng)濟(jì)性,但是最大俯仰力矩增大,降低了水下滑翔器的穩(wěn)定性。圖25是航速0.5 kn,不同攻角下具體模型和樣機(jī)升阻比曲線對(duì)照?qǐng)D。由圖25中可知,模型2和模型3的升阻比在0°~40°的攻角范圍內(nèi)一直大于樣機(jī),兩個(gè)模型的水動(dòng)力性能提升非常明顯,均滿足設(shè)計(jì)要求。同時(shí),模型2和模型3在攻角8°左右時(shí)具有最大的升阻比,故在8°左右攻角時(shí)具有最好的滑翔經(jīng)濟(jì)性。表8是模型2、模型3和樣機(jī)達(dá)到最大升阻比時(shí)的水動(dòng)力參數(shù)對(duì)照表。

    表8 模型2和模型3與樣機(jī)最大升阻比時(shí)水動(dòng)力參數(shù)對(duì)照表Tab.8 The hydraulic performance comparison of specific models and the classic model with maximal L/D ratio

    為了綜合考慮升阻比和俯仰力矩這兩個(gè)影響水下滑翔器水動(dòng)力性能的因素,本文引用文獻(xiàn)[22]中的加權(quán)公式:

    式中:FL(X)、FD(X)和 Mpitch( X )分別是水下滑翔器在攻角8°受到的升力、阻力和最大俯仰力矩,ω1和ω2是加權(quán)系數(shù)。

    為了能夠保證持續(xù)高效航行,一般會(huì)讓水下滑翔器的攻角長(zhǎng)時(shí)間保持在最大升阻比時(shí)的角度,這時(shí)水下滑翔器具有最佳滑翔經(jīng)濟(jì)性,因此水下滑翔器在這一角度的穩(wěn)定性也尤為重要。利用加權(quán)公式綜合考慮這兩大因素,獲得的加權(quán)值,如表8所示。其中模型2的加權(quán)值最高,即模型2的綜合水動(dòng)力性能最優(yōu)。以上數(shù)據(jù)分析表明,模型2在保持長(zhǎng)時(shí)間航行時(shí)具有更高的航速以及工作效率。

    4 結(jié)語(yǔ)

    在Slocum等幾種典型水下滑翔器樣機(jī)的基礎(chǔ)上,對(duì)主體及其附體機(jī)翼和尾翼進(jìn)行了設(shè)計(jì),得到了一種全新的水下滑翔器外形,它的水動(dòng)力性能相比于Spray、Seaglider、Slocum三大主流樣機(jī)有了大幅度的提升。當(dāng)滑翔器在8°左右攻角航行時(shí),具有最大的升阻比,可以保證其在水下長(zhǎng)時(shí)間高效率航行,并且能滿足工程設(shè)計(jì)和研究的要求。

    [1] Daneil L Rudnick,Russ E Davis,Charles C Eriksen,et al.Underwater gliders for ocean research[J].Marine Technology,Society Journal,2004,38(1):48-59.

    [2] 張宏偉,李智,王延輝,等.水下滑翔器浮力驅(qū)動(dòng)機(jī)構(gòu)布局分析[J].海洋技術(shù),2009,28(1):1-4.(ZHANG Hong-wei,LI Zhi,WANG Yan-hui,et al.Research on the layout of the buoyancy engine of underwater glider[J].Ocean Technology,2009,28(1):1-4.(in Chinese))

    [3] 王樹(shù)新,王延輝,張大濤,等.溫差能驅(qū)動(dòng)的水下滑翔器設(shè)計(jì)與實(shí)驗(yàn)研究[J].海洋技術(shù),2006,25(1):1-5.(WANG Shu-xin,WANG Yan-hui,ZHANG Da-tao,et al.Design and trial on an underwater glider propelled by thermal engine[J].Ocean Technology,2006,25(1):1-5.(in Chinese))

    [4] 張奇峰,俞建成,唐元貴,等.水下滑翔機(jī)器人運(yùn)動(dòng)分析與載體設(shè)計(jì)[J].海洋工程,2006,24(1):74-78.(ZHANG Qi-feng,YU Jian-cheng,TANG Yuan-gui,et al.Motion analysis and body design of an underwater glider[J].The Ocean Engineering,2006,24(1):74-78.(in Chinese))

    [5] Douglas C W,Paul J S,Clayton P J.SLOCUM:An underwater glider propelled by environmental energy[J].IEEE Journal of Oceanic Engineering,2001,26(4):447-452.

    [6] Eriksen C C,Osse T J,Light R D,et al.Seaglider:A long range autonomous underwater vehicle for oceanographic research[J].IEEE Journal of Oceanic Engineering,2001,26(4):424-436.

    [7] Sherman J,Davis R E,Owens W B,et al.The autonomous underwater glider“Spray”[J].IEEE Journal of Oceanic Engineering,2001,26(4):437-446.

    [8] Davis R E,Eriksen C E,Jones C P.The technology and applications of autonomous underwater vehicles[M].Griffiths Ged.London:Taylor and Francis,2002.

    [9] Dr Stephen Wood,P E ToddAllen,Sean Kuhn Joseph Caldwell.The Development of an autonomous underwater powered glider for deep-sea biological,chemical and physical oceanography[C]//OCEANS 2007-Europe.IEEE Journal of Oceanic Engineering,2007:1-6.

    [10] Chappell S G,Komerska R J,Blidberg D R,et al.Recent field experience with multiple cooperating solar-powered AUVs[C]//The Fifteenth International Symposium on Unmanned Untethered Submersible Technology.IEEE Journal of Oceanic Engineering,2007:1-8.

    [11] Arima M,Okashima T,Yamada T.Development of a solar-powered underwater glider[C]//Underwater Technology(UT),2011 IEEE Symposium on and 2011 Workshop on Scientific Use of Submarine Cables and Related Technologies(SSC).IEEE Journal of Oceanic Engineering,2011:1-5.

    [12] Alvarez A.Redesigning the SLOCUM glider for torpedo tube launching[J].IEEE Journal of Oceanic Engineering,2010,35(4):984-991.

    [13]張禹,李雀屏,楊國(guó)哲,等.水下滑翔機(jī)器人載體外形設(shè)計(jì)與優(yōu)化[J].機(jī)械設(shè)計(jì)與制造,2010,10:148-150.(ZHANG Yu,LI Que-ping,YANG Guo-zhe,et al.Design and optimum shape of underwater glider[J].Machinery Design & Manufacture,2010,10:148-150.(in Chinese))

    [14]李鵬飛,徐敏義,王飛飛.精通CFD工程仿真與案例實(shí)戰(zhàn)[M].北京:人民郵電出版社,2011.(LI Peng-fei,XU Min-yi,WANG Fei-fei.Fluent gambit icem CFD tecplot[M].Beijing:Posts & Telecom Press,2011.(in Chinese))

    [15]Joshua G G,Ralf B,Naomi E L.Underwater glider model parameter identification[C]//13th Unmanned Untethered Submersible Technology(UUST).2003:1-12.

    [16]武建國(guó),陳超英,王樹(shù)新,等.混合驅(qū)動(dòng)水下滑翔器滑翔狀態(tài)機(jī)翼水動(dòng)力特性[J].天津大學(xué)學(xué)報(bào),2010,43(1):84-89.(WU Jian-guo,CHEN Chao-ying,WANG Shu-xin,et al.Hydrodynamic characteristics of the wings of hybrid-driven underwater glider in glide mode[J].Journal of Tianjin University,2010,43(1):84-89.(in Chinese))

    [17]胡克,俞建成,張奇峰.水下滑翔機(jī)器人載體外形設(shè)計(jì)與優(yōu)化[J].機(jī)器人,2005,27(2):108-113.(HU Ke,YU Jian-cheng,ZHANG Qi-feng.Design and optimization of underwater glider shape[J].Robot,2005,27(2):108-113.(in Chinese))

    [18]王強(qiáng),葛彤,吳超,等.結(jié)合固定翼飛機(jī)原理重水潛器設(shè)計(jì)[J].海洋工程,2012,30(2):143-149.(WANG Qiang,GE Tong,WU Chao,et al.Design of the HUV based on the airplane's principles[J].The Ocean Engineering,2012,30(2):143-149.(in Chinese))

    [19]諸敏.水下滑翔機(jī)設(shè)計(jì)優(yōu)化與運(yùn)動(dòng)分析[D].杭州:浙江大學(xué),2007.(ZHU Min.Optimizing&kinetic study of underwater glide[D].Hangzhou:Zhejiang University,2007.(in Chinese))

    [20] Daniel L R,Russ E D,Charles C E,et al.Underwater gliders for ocean research[J].Marine Technology Society Journal,2004,38(1):48-59.

    [21] Davis R E,Eriksen C C,Jones C P.Autonomous buoyancy-driven underwater gliders[C]//Technology and Applications of Autonomous Underwater Vehicles.London:Taylor& Francis,2002:37-58.

    [22]谷海濤,林揚(yáng),胡志強(qiáng),等.基于代理模型的水下滑翔機(jī)機(jī)翼設(shè)計(jì)優(yōu)化方法[J].機(jī)械工程學(xué)報(bào),2009,45(12):7-14.(GU Hai-tao,LIN Yang,HU Zhi-qiang,et al.Surrogate models based optimization methods for the design of underwater glider wing[J].Journal of Mechanical Engineering,2009,45(12):7-14.(in Chinese))

    猜你喜歡
    后掠角滑翔攻角
    蘭花螳螂會(huì)滑翔
    不同后掠角大展弦比復(fù)合材料機(jī)翼氣動(dòng)特性
    攻天掠地的先鋒武器——滑翔導(dǎo)彈
    風(fēng)標(biāo)式攻角傳感器在超聲速飛行運(yùn)載火箭中的應(yīng)用研究
    大攻角狀態(tài)壓氣機(jī)分離流及葉片動(dòng)力響應(yīng)特性
    讓你心跳加速的“滑翔飛板”
    基于CFD的最優(yōu)變后掠規(guī)律研究
    不同后掠角三角翼的靜態(tài)地面效應(yīng)數(shù)值模擬
    空中滑翔大比拼(下)——滑翔傘
    附加攻角效應(yīng)對(duì)顫振穩(wěn)定性能影響
    嫩草影院新地址| 国产综合懂色| 午夜激情欧美在线| 边亲边吃奶的免费视频| 久久国内精品自在自线图片| 国产日韩欧美在线精品| 99热6这里只有精品| 内地一区二区视频在线| 乱人视频在线观看| 国产一区亚洲一区在线观看| 国产极品天堂在线| 国产av码专区亚洲av| 久久精品国产鲁丝片午夜精品| 亚洲av免费高清在线观看| 亚洲国产精品sss在线观看| 亚洲精品456在线播放app| 变态另类丝袜制服| 一级爰片在线观看| 插逼视频在线观看| 最后的刺客免费高清国语| 哪个播放器可以免费观看大片| 亚洲欧美日韩卡通动漫| 国产国拍精品亚洲av在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av专区在线播放| av国产久精品久网站免费入址| 秋霞在线观看毛片| 国产 一区 欧美 日韩| 亚洲国产日韩欧美精品在线观看| 在线天堂最新版资源| 村上凉子中文字幕在线| 久久精品综合一区二区三区| 亚洲婷婷狠狠爱综合网| 一级av片app| 亚洲三级黄色毛片| 国产毛片a区久久久久| av天堂中文字幕网| 97人妻精品一区二区三区麻豆| 我的女老师完整版在线观看| 一个人免费在线观看电影| 久99久视频精品免费| 免费黄网站久久成人精品| 久久这里有精品视频免费| 国产毛片a区久久久久| 精品人妻偷拍中文字幕| 黑人高潮一二区| 国产视频首页在线观看| 免费观看a级毛片全部| 成人性生交大片免费视频hd| 十八禁国产超污无遮挡网站| 国产精品一及| 国产亚洲一区二区精品| 午夜福利在线观看免费完整高清在| 26uuu在线亚洲综合色| 亚洲精品456在线播放app| 美女内射精品一级片tv| 日韩中字成人| 国产av在哪里看| 一个人看的www免费观看视频| 欧美日韩国产亚洲二区| av专区在线播放| 国产一级毛片在线| 高清av免费在线| 全区人妻精品视频| 亚洲怡红院男人天堂| 中文字幕人妻熟人妻熟丝袜美| 日韩制服骚丝袜av| 一个人看视频在线观看www免费| 国产精品嫩草影院av在线观看| 精品久久国产蜜桃| 国产精品1区2区在线观看.| 久久欧美精品欧美久久欧美| 精品99又大又爽又粗少妇毛片| 超碰av人人做人人爽久久| 国产真实伦视频高清在线观看| 久久久午夜欧美精品| 亚洲av.av天堂| 日本免费a在线| 一级黄片播放器| 国产在线一区二区三区精 | 日本三级黄在线观看| 国产精品野战在线观看| 小说图片视频综合网站| 草草在线视频免费看| 搡老妇女老女人老熟妇| 男人狂女人下面高潮的视频| 久久久国产成人免费| 能在线免费观看的黄片| 亚洲人成网站在线观看播放| 搞女人的毛片| 精品酒店卫生间| 色尼玛亚洲综合影院| 看黄色毛片网站| 综合色丁香网| 国产免费一级a男人的天堂| 免费黄色在线免费观看| 床上黄色一级片| av在线亚洲专区| eeuss影院久久| 国产老妇女一区| 亚洲av不卡在线观看| 熟女人妻精品中文字幕| 一级黄片播放器| 国产精品国产三级国产专区5o | 久久久久久久国产电影| 精品人妻一区二区三区麻豆| 我要搜黄色片| 国产精品1区2区在线观看.| 亚洲国产精品久久男人天堂| 日本与韩国留学比较| 欧美一级a爱片免费观看看| 男人舔女人下体高潮全视频| 小蜜桃在线观看免费完整版高清| 精品人妻一区二区三区麻豆| a级一级毛片免费在线观看| 少妇熟女aⅴ在线视频| 91午夜精品亚洲一区二区三区| 最近最新中文字幕免费大全7| 色播亚洲综合网| 免费看日本二区| 亚洲无线观看免费| av线在线观看网站| 亚洲三级黄色毛片| 纵有疾风起免费观看全集完整版 | 成人性生交大片免费视频hd| 国产成人精品一,二区| 一个人看的www免费观看视频| 2022亚洲国产成人精品| 国产精品女同一区二区软件| 久久久久九九精品影院| 国产单亲对白刺激| 国产淫语在线视频| 国内精品美女久久久久久| 亚洲高清免费不卡视频| 三级国产精品片| 真实男女啪啪啪动态图| 波野结衣二区三区在线| 国产毛片a区久久久久| 啦啦啦啦在线视频资源| 免费一级毛片在线播放高清视频| 欧美性猛交╳xxx乱大交人| 国产亚洲精品av在线| 2021天堂中文幕一二区在线观| 久99久视频精品免费| АⅤ资源中文在线天堂| 欧美日韩精品成人综合77777| 一区二区三区高清视频在线| 久久鲁丝午夜福利片| 欧美日韩综合久久久久久| 国产日韩欧美在线精品| 欧美高清性xxxxhd video| 狠狠狠狠99中文字幕| 久久精品91蜜桃| 色尼玛亚洲综合影院| 99久国产av精品国产电影| 日本-黄色视频高清免费观看| 亚洲人成网站在线播| 国产精华一区二区三区| 舔av片在线| 精品熟女少妇av免费看| 亚洲丝袜综合中文字幕| 精品欧美国产一区二区三| 欧美一区二区国产精品久久精品| 少妇丰满av| 一个人看视频在线观看www免费| 亚洲成人av在线免费| 男人舔女人下体高潮全视频| av在线播放精品| a级毛色黄片| 日韩欧美在线乱码| 国产精品综合久久久久久久免费| 在线免费观看的www视频| 久久韩国三级中文字幕| 国产精品伦人一区二区| 中国美白少妇内射xxxbb| 夜夜看夜夜爽夜夜摸| 亚洲激情五月婷婷啪啪| 99久久精品热视频| 国产亚洲av嫩草精品影院| 亚洲图色成人| 黄色一级大片看看| 日韩精品有码人妻一区| 禁无遮挡网站| 国产成人午夜福利电影在线观看| 日本与韩国留学比较| 在线观看美女被高潮喷水网站| 久久午夜福利片| 精品久久久久久久久av| 亚洲国产精品专区欧美| 一个人看的www免费观看视频| 麻豆av噜噜一区二区三区| 亚洲自偷自拍三级| 中文字幕免费在线视频6| 综合色丁香网| 国产不卡一卡二| 亚洲在久久综合| 日韩高清综合在线| 99国产精品一区二区蜜桃av| 简卡轻食公司| 青青草视频在线视频观看| 日韩成人av中文字幕在线观看| 一级二级三级毛片免费看| 亚洲精品乱久久久久久| 在线a可以看的网站| 亚洲欧美日韩东京热| 国产精品1区2区在线观看.| 精品少妇黑人巨大在线播放 | 久久久精品94久久精品| 久久亚洲国产成人精品v| 亚洲天堂国产精品一区在线| 中文精品一卡2卡3卡4更新| 精品久久国产蜜桃| 亚洲精品国产av成人精品| 国产国拍精品亚洲av在线观看| 深夜a级毛片| 国产在线男女| 极品教师在线视频| 最近2019中文字幕mv第一页| 男女视频在线观看网站免费| 高清在线视频一区二区三区 | 边亲边吃奶的免费视频| 国产一区二区在线av高清观看| 亚洲国产精品成人久久小说| 亚洲综合精品二区| 午夜福利成人在线免费观看| h日本视频在线播放| 国产精品久久久久久久电影| www.色视频.com| 最后的刺客免费高清国语| 成人综合一区亚洲| 亚洲欧美精品专区久久| 日日干狠狠操夜夜爽| 欧美成人免费av一区二区三区| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 精品久久久久久久末码| 国产免费又黄又爽又色| 91aial.com中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 日本三级黄在线观看| 国产高清不卡午夜福利| 丰满人妻一区二区三区视频av| 天堂√8在线中文| 国产成年人精品一区二区| 91aial.com中文字幕在线观看| 欧美不卡视频在线免费观看| 超碰97精品在线观看| 日日啪夜夜撸| 国产免费又黄又爽又色| 99热全是精品| 最近视频中文字幕2019在线8| 精品久久久久久久久久久久久| 免费黄色在线免费观看| 一本一本综合久久| 床上黄色一级片| 亚洲av电影在线观看一区二区三区 | 国产黄色小视频在线观看| 国产真实伦视频高清在线观看| 午夜免费男女啪啪视频观看| 精品无人区乱码1区二区| 日日干狠狠操夜夜爽| 久久99精品国语久久久| 一级黄色大片毛片| 中文在线观看免费www的网站| 国内精品宾馆在线| 亚洲自偷自拍三级| 亚洲五月天丁香| 国产精品伦人一区二区| 亚洲国产成人一精品久久久| 夜夜爽夜夜爽视频| 久久人妻av系列| 一区二区三区高清视频在线| 男女啪啪激烈高潮av片| www.色视频.com| 美女黄网站色视频| 菩萨蛮人人尽说江南好唐韦庄 | 少妇被粗大猛烈的视频| 亚洲精品自拍成人| 亚洲av一区综合| 尤物成人国产欧美一区二区三区| 在线免费观看不下载黄p国产| av黄色大香蕉| 91午夜精品亚洲一区二区三区| 色吧在线观看| 久久久国产成人精品二区| 最近最新中文字幕大全电影3| 一级黄片播放器| 国产精品久久视频播放| 久久6这里有精品| 亚洲aⅴ乱码一区二区在线播放| .国产精品久久| 亚洲av电影不卡..在线观看| 日日撸夜夜添| 欧美不卡视频在线免费观看| av.在线天堂| 高清毛片免费看| 精品少妇黑人巨大在线播放 | 成人鲁丝片一二三区免费| 狠狠狠狠99中文字幕| 变态另类丝袜制服| 亚洲精品影视一区二区三区av| 中文字幕制服av| 又黄又爽又刺激的免费视频.| 51国产日韩欧美| 国产亚洲av片在线观看秒播厂 | 国产精品不卡视频一区二区| 九草在线视频观看| 欧美日本视频| 天堂影院成人在线观看| 99热这里只有精品一区| av线在线观看网站| 久久久精品大字幕| 国产熟女欧美一区二区| 国产免费福利视频在线观看| 欧美性猛交黑人性爽| 成人性生交大片免费视频hd| 91狼人影院| 亚洲国产精品sss在线观看| 毛片一级片免费看久久久久| 午夜福利成人在线免费观看| 成人亚洲欧美一区二区av| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 日本一本二区三区精品| 亚洲乱码一区二区免费版| 亚洲aⅴ乱码一区二区在线播放| 国产探花在线观看一区二区| 欧美精品一区二区大全| 最近中文字幕高清免费大全6| 成人美女网站在线观看视频| 搞女人的毛片| 日韩三级伦理在线观看| 极品教师在线视频| 22中文网久久字幕| 欧美成人一区二区免费高清观看| 别揉我奶头 嗯啊视频| 中文亚洲av片在线观看爽| 国产亚洲91精品色在线| 亚洲av不卡在线观看| 两个人的视频大全免费| 国产午夜精品久久久久久一区二区三区| 少妇猛男粗大的猛烈进出视频 | 激情 狠狠 欧美| 欧美xxxx黑人xx丫x性爽| 91久久精品国产一区二区成人| 高清在线视频一区二区三区 | 国产成人精品婷婷| 一区二区三区高清视频在线| 男女下面进入的视频免费午夜| 色噜噜av男人的天堂激情| 搡女人真爽免费视频火全软件| 天天躁夜夜躁狠狠久久av| 久久精品国产自在天天线| 免费看美女性在线毛片视频| 国产一区有黄有色的免费视频 | 国产成人精品婷婷| 1024手机看黄色片| 精品酒店卫生间| 日本黄大片高清| 最近视频中文字幕2019在线8| 亚洲人成网站在线播| 边亲边吃奶的免费视频| 99久久无色码亚洲精品果冻| 免费观看精品视频网站| 久久精品夜夜夜夜夜久久蜜豆| 三级毛片av免费| 又爽又黄无遮挡网站| 亚洲色图av天堂| 久久久久久久午夜电影| 国产高清三级在线| 久久久亚洲精品成人影院| 美女国产视频在线观看| 国产高清有码在线观看视频| 欧美三级亚洲精品| 久久精品国产自在天天线| 99久久无色码亚洲精品果冻| 99在线人妻在线中文字幕| 精品欧美国产一区二区三| 国产日韩欧美在线精品| 干丝袜人妻中文字幕| 国产精品人妻久久久久久| 国产精品不卡视频一区二区| 久久久久免费精品人妻一区二区| av.在线天堂| 一级毛片久久久久久久久女| 国产色爽女视频免费观看| 国产极品天堂在线| 国产亚洲av嫩草精品影院| 成人高潮视频无遮挡免费网站| 久久精品影院6| 能在线免费看毛片的网站| 在线播放无遮挡| 亚洲自偷自拍三级| 国产v大片淫在线免费观看| 在线观看av片永久免费下载| 18禁在线无遮挡免费观看视频| 久久久久免费精品人妻一区二区| 国产精品日韩av在线免费观看| 色尼玛亚洲综合影院| 精品人妻视频免费看| 水蜜桃什么品种好| 亚洲欧美成人精品一区二区| 天堂√8在线中文| 老司机影院成人| 99久久九九国产精品国产免费| 一夜夜www| 国产探花在线观看一区二区| 久久久久久伊人网av| 美女内射精品一级片tv| 亚洲在久久综合| 男女国产视频网站| 国产黄片视频在线免费观看| 亚洲丝袜综合中文字幕| 九九久久精品国产亚洲av麻豆| 久久久久性生活片| 99久久精品一区二区三区| 国产探花在线观看一区二区| 亚洲成人中文字幕在线播放| 久久久久性生活片| 免费观看精品视频网站| www.色视频.com| 精品免费久久久久久久清纯| 午夜视频国产福利| 亚洲国产欧洲综合997久久,| 日韩欧美国产在线观看| 只有这里有精品99| 99久久精品国产国产毛片| 三级经典国产精品| 91久久精品国产一区二区成人| av专区在线播放| 亚洲欧美日韩无卡精品| 精品人妻熟女av久视频| 国产精品一区二区性色av| 床上黄色一级片| 久久精品国产鲁丝片午夜精品| 亚洲成人中文字幕在线播放| 岛国在线免费视频观看| 秋霞伦理黄片| 亚洲国产精品sss在线观看| 精品人妻熟女av久视频| 午夜福利在线观看免费完整高清在| 一个人观看的视频www高清免费观看| 国产精品不卡视频一区二区| 中文精品一卡2卡3卡4更新| 精品人妻熟女av久视频| 男女国产视频网站| 日本三级黄在线观看| 99热网站在线观看| 国产大屁股一区二区在线视频| 国产男人的电影天堂91| 天堂网av新在线| 99在线人妻在线中文字幕| 又黄又爽又刺激的免费视频.| 国产高清视频在线观看网站| 亚洲国产成人一精品久久久| 亚洲伊人久久精品综合 | 国产av一区在线观看免费| 只有这里有精品99| 欧美成人精品欧美一级黄| 久久婷婷人人爽人人干人人爱| 中文字幕制服av| av免费观看日本| 免费黄色在线免费观看| 午夜福利在线在线| 亚洲精品456在线播放app| 我的老师免费观看完整版| 久久国内精品自在自线图片| 日本黄大片高清| 午夜激情福利司机影院| 人人妻人人澡欧美一区二区| 欧美极品一区二区三区四区| 日韩精品有码人妻一区| 天堂网av新在线| 亚洲av福利一区| 久久久欧美国产精品| 18禁在线播放成人免费| 伦精品一区二区三区| 精品久久久久久久末码| h日本视频在线播放| 欧美色视频一区免费| 欧美日本亚洲视频在线播放| 免费看日本二区| 我的女老师完整版在线观看| 18+在线观看网站| 国产视频首页在线观看| 精品酒店卫生间| 国产精华一区二区三区| 国产一级毛片七仙女欲春2| 青春草国产在线视频| 草草在线视频免费看| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 午夜免费男女啪啪视频观看| 国产免费视频播放在线视频 | 欧美一区二区精品小视频在线| 久久国内精品自在自线图片| 久久久午夜欧美精品| 好男人在线观看高清免费视频| 国产女主播在线喷水免费视频网站 | 欧美性猛交黑人性爽| 赤兔流量卡办理| 嫩草影院精品99| 精品久久久久久久人妻蜜臀av| 搡老妇女老女人老熟妇| 中文亚洲av片在线观看爽| 长腿黑丝高跟| 国产淫语在线视频| 在线观看av片永久免费下载| 51国产日韩欧美| 日韩av在线免费看完整版不卡| 亚洲欧美精品综合久久99| 美女国产视频在线观看| 国产色爽女视频免费观看| АⅤ资源中文在线天堂| 欧美人与善性xxx| 久久久精品欧美日韩精品| 九草在线视频观看| 毛片女人毛片| 国产免费一级a男人的天堂| 美女高潮的动态| 在线观看66精品国产| 久久韩国三级中文字幕| 久久人人爽人人片av| 国产伦理片在线播放av一区| 日日摸夜夜添夜夜添av毛片| 亚洲真实伦在线观看| 麻豆成人午夜福利视频| 国产免费男女视频| 久久久久久大精品| 天堂影院成人在线观看| 26uuu在线亚洲综合色| 国产真实乱freesex| 国产女主播在线喷水免费视频网站 | 久久久精品大字幕| 久久综合国产亚洲精品| 午夜免费男女啪啪视频观看| 国产真实伦视频高清在线观看| 亚洲一区高清亚洲精品| 亚洲国产精品sss在线观看| 国产片特级美女逼逼视频| 久久久久久九九精品二区国产| 免费看av在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av成人精品一区久久| 一级毛片久久久久久久久女| 国产免费男女视频| 国产亚洲精品久久久com| 国产精品无大码| 男人狂女人下面高潮的视频| 国产欧美另类精品又又久久亚洲欧美| 秋霞伦理黄片| 亚洲av熟女| 国产精品一区二区三区四区免费观看| 日本爱情动作片www.在线观看| 日本黄色视频三级网站网址| 乱系列少妇在线播放| 日韩中字成人| 国产在线男女| 一个人免费在线观看电影| 三级国产精品欧美在线观看| 男插女下体视频免费在线播放| 亚洲成人av在线免费| 免费在线观看成人毛片| 亚洲精品成人久久久久久| 国产在视频线在精品| 噜噜噜噜噜久久久久久91| 波野结衣二区三区在线| 国产亚洲av嫩草精品影院| 亚洲精品,欧美精品| 亚洲人成网站在线播| 国国产精品蜜臀av免费| 十八禁国产超污无遮挡网站| 狠狠狠狠99中文字幕| 一区二区三区高清视频在线| 51国产日韩欧美| 色吧在线观看| 久久久久免费精品人妻一区二区| 免费观看性生交大片5| 1000部很黄的大片| 亚洲精品456在线播放app| 精品国产三级普通话版| 成人一区二区视频在线观看| 久久韩国三级中文字幕| 国产精品久久久久久久电影| 国产精品国产三级国产专区5o | 简卡轻食公司| 91av网一区二区| 青春草视频在线免费观看| 一二三四中文在线观看免费高清| 久久精品91蜜桃| 亚洲av电影不卡..在线观看| 国产一区二区在线av高清观看| 国产精品.久久久| 日本欧美国产在线视频| 91av网一区二区| 可以在线观看毛片的网站| 我要看日韩黄色一级片| av在线蜜桃| 男人狂女人下面高潮的视频| 国产成人免费观看mmmm| 特级一级黄色大片| 亚洲婷婷狠狠爱综合网| 亚洲成人精品中文字幕电影| 国产免费又黄又爽又色| 一个人免费在线观看电影| 日韩欧美 国产精品| 久久精品国产99精品国产亚洲性色| 国产精品,欧美在线| 精品少妇黑人巨大在线播放 | 在线观看av片永久免费下载| 婷婷色av中文字幕| 长腿黑丝高跟| 两性午夜刺激爽爽歪歪视频在线观看| 插阴视频在线观看视频| 亚洲av中文av极速乱|