• <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)定性能影響
    婷婷成人精品国产| 永久网站在线| 亚洲熟女精品中文字幕| 欧美精品国产亚洲| 日韩免费高清中文字幕av| 久热久热在线精品观看| 国产一区二区激情短视频 | 色播在线永久视频| 亚洲少妇的诱惑av| 天堂中文最新版在线下载| 亚洲国产av新网站| 99久久综合免费| 女人被躁到高潮嗷嗷叫费观| 久久久久精品人妻al黑| 69精品国产乱码久久久| 91在线精品国自产拍蜜月| 韩国av在线不卡| 久热这里只有精品99| 久久人妻熟女aⅴ| 少妇人妻久久综合中文| 高清不卡的av网站| 欧美成人午夜免费资源| 成人亚洲欧美一区二区av| 国产极品粉嫩免费观看在线| 亚洲欧美一区二区三区久久| 国产日韩欧美亚洲二区| 国产成人精品久久久久久| 色吧在线观看| 女性被躁到高潮视频| 丝袜喷水一区| 成年美女黄网站色视频大全免费| 高清av免费在线| 99久国产av精品国产电影| 亚洲av在线观看美女高潮| 国产精品.久久久| 观看av在线不卡| 欧美日韩视频高清一区二区三区二| 侵犯人妻中文字幕一二三四区| 日韩电影二区| 亚洲精品久久久久久婷婷小说| 欧美另类一区| 国产日韩欧美视频二区| 十八禁网站网址无遮挡| 日韩精品有码人妻一区| 亚洲精品国产av成人精品| 精品人妻一区二区三区麻豆| 香蕉丝袜av| 亚洲av.av天堂| 精品卡一卡二卡四卡免费| 狠狠精品人妻久久久久久综合| 老汉色av国产亚洲站长工具| 色婷婷av一区二区三区视频| 女性被躁到高潮视频| 亚洲三级黄色毛片| 久久99一区二区三区| 国产男人的电影天堂91| 久久精品夜色国产| a级片在线免费高清观看视频| 18在线观看网站| 中文字幕最新亚洲高清| 久久久久久人人人人人| 精品国产一区二区三区久久久樱花| 99热网站在线观看| 热99国产精品久久久久久7| 少妇精品久久久久久久| a级毛片在线看网站| 亚洲欧洲精品一区二区精品久久久 | 99热国产这里只有精品6| 亚洲色图综合在线观看| 国产亚洲av片在线观看秒播厂| 精品国产露脸久久av麻豆| 亚洲综合色惰| 欧美日韩精品网址| 亚洲在久久综合| 99香蕉大伊视频| 啦啦啦在线免费观看视频4| 午夜精品国产一区二区电影| 国产有黄有色有爽视频| 欧美国产精品一级二级三级| 国产精品国产三级国产专区5o| 欧美黄色片欧美黄色片| 久久久精品国产亚洲av高清涩受| 亚洲国产精品成人久久小说| 国产不卡av网站在线观看| 亚洲精品久久久久久婷婷小说| 国产一区二区 视频在线| 亚洲成人一二三区av| 成年人免费黄色播放视频| 69精品国产乱码久久久| 欧美 亚洲 国产 日韩一| 日本欧美视频一区| 亚洲精品日韩在线中文字幕| 最近中文字幕高清免费大全6| 在线亚洲精品国产二区图片欧美| 欧美日韩一区二区视频在线观看视频在线| 在线观看三级黄色| 一级a爱视频在线免费观看| kizo精华| 搡女人真爽免费视频火全软件| 欧美日韩国产mv在线观看视频| 丁香六月天网| 亚洲精品日韩在线中文字幕| 亚洲精品久久午夜乱码| 亚洲精品一区蜜桃| 美女大奶头黄色视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲综合精品二区| 精品视频人人做人人爽| 免费人妻精品一区二区三区视频| 日日摸夜夜添夜夜爱| 亚洲国产精品国产精品| 日韩av不卡免费在线播放| av天堂久久9| 黄色 视频免费看| 成人亚洲精品一区在线观看| 国产激情久久老熟女| 亚洲五月色婷婷综合| 欧美黄色片欧美黄色片| 日韩大片免费观看网站| 少妇熟女欧美另类| 黄色一级大片看看| 中文精品一卡2卡3卡4更新| a 毛片基地| 亚洲欧美成人精品一区二区| 久久久精品国产亚洲av高清涩受| 亚洲欧美精品综合一区二区三区 | 大片电影免费在线观看免费| 国产一区二区 视频在线| 飞空精品影院首页| 寂寞人妻少妇视频99o| 欧美国产精品一级二级三级| 三上悠亚av全集在线观看| 国产高清不卡午夜福利| 久久ye,这里只有精品| 国产精品一二三区在线看| 婷婷色综合大香蕉| 久久精品久久精品一区二区三区| 精品亚洲成a人片在线观看| 熟女av电影| 欧美人与性动交α欧美精品济南到 | 纯流量卡能插随身wifi吗| 欧美日韩综合久久久久久| xxx大片免费视频| 秋霞在线观看毛片| 久久精品人人爽人人爽视色| 多毛熟女@视频| 青春草国产在线视频| 国产成人精品福利久久| 亚洲国产欧美在线一区| 丝袜在线中文字幕| 高清黄色对白视频在线免费看| 一级毛片电影观看| 国产乱人偷精品视频| 永久免费av网站大全| 久久人人爽av亚洲精品天堂| 欧美av亚洲av综合av国产av | 中文字幕最新亚洲高清| 老熟女久久久| a级毛片在线看网站| 精品少妇久久久久久888优播| 国产在线一区二区三区精| 亚洲精品aⅴ在线观看| 久久影院123| 狠狠精品人妻久久久久久综合| 黄频高清免费视频| 日本av手机在线免费观看| 精品亚洲成a人片在线观看| 日韩av不卡免费在线播放| 啦啦啦中文免费视频观看日本| 日本91视频免费播放| 伊人亚洲综合成人网| 精品少妇久久久久久888优播| 国产野战对白在线观看| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 国产野战对白在线观看| 成人黄色视频免费在线看| 亚洲综合精品二区| 1024香蕉在线观看| 一本—道久久a久久精品蜜桃钙片| 狂野欧美激情性bbbbbb| 最近中文字幕2019免费版| 国产精品欧美亚洲77777| 超碰97精品在线观看| 亚洲成人av在线免费| 欧美国产精品一级二级三级| 国产亚洲精品第一综合不卡| 亚洲色图 男人天堂 中文字幕| 午夜老司机福利剧场| 国产av一区二区精品久久| 亚洲熟女精品中文字幕| 国产黄色视频一区二区在线观看| 久久国产精品大桥未久av| 欧美另类一区| 免费黄网站久久成人精品| 高清黄色对白视频在线免费看| 日韩成人av中文字幕在线观看| 亚洲国产毛片av蜜桃av| 国产黄频视频在线观看| a级毛片黄视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产av影院在线观看| 亚洲第一区二区三区不卡| 一区二区av电影网| 天堂8中文在线网| 在线观看www视频免费| 久久精品aⅴ一区二区三区四区 | 亚洲精品,欧美精品| 9191精品国产免费久久| 欧美+日韩+精品| 亚洲一级一片aⅴ在线观看| 亚洲五月色婷婷综合| 国产精品一区二区在线观看99| 免费av中文字幕在线| 纯流量卡能插随身wifi吗| 欧美精品一区二区大全| 最近中文字幕高清免费大全6| 国产xxxxx性猛交| 亚洲精品日韩在线中文字幕| 观看美女的网站| 性色av一级| 少妇的逼水好多| av网站免费在线观看视频| 国产一区亚洲一区在线观看| 91精品伊人久久大香线蕉| 欧美日韩亚洲高清精品| 看免费成人av毛片| 大香蕉久久成人网| 国产视频首页在线观看| 成年人午夜在线观看视频| 国产乱人偷精品视频| 欧美最新免费一区二区三区| 亚洲伊人久久精品综合| 91精品伊人久久大香线蕉| 乱人伦中国视频| 人体艺术视频欧美日本| 亚洲成国产人片在线观看| 丝袜美足系列| 免费在线观看完整版高清| 日韩中文字幕视频在线看片| 久久久久久久国产电影| 日本av免费视频播放| 久久久久网色| 色网站视频免费| 成人漫画全彩无遮挡| 午夜影院在线不卡| 色婷婷av一区二区三区视频| 中文字幕亚洲精品专区| 最新的欧美精品一区二区| 亚洲欧美清纯卡通| 久久精品国产亚洲av涩爱| 国产激情久久老熟女| 国产在线一区二区三区精| 久久99蜜桃精品久久| 在线 av 中文字幕| 你懂的网址亚洲精品在线观看| 最近中文字幕高清免费大全6| 日韩视频在线欧美| 成人国语在线视频| 丝瓜视频免费看黄片| 日产精品乱码卡一卡2卡三| 久久国内精品自在自线图片| www.av在线官网国产| 国产亚洲最大av| 亚洲欧美色中文字幕在线| 99热国产这里只有精品6| 一边摸一边做爽爽视频免费| av网站免费在线观看视频| 9191精品国产免费久久| videosex国产| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 国产精品麻豆人妻色哟哟久久| 日韩 亚洲 欧美在线| 日日啪夜夜爽| 久久国产精品大桥未久av| 精品一区二区免费观看| 人妻人人澡人人爽人人| 精品亚洲成a人片在线观看| 寂寞人妻少妇视频99o| 777米奇影视久久| 色婷婷久久久亚洲欧美| 日韩视频在线欧美| 交换朋友夫妻互换小说| 免费av中文字幕在线| 国产精品久久久久久久久免| 亚洲人成77777在线视频| 一边亲一边摸免费视频| 亚洲图色成人| 中文字幕人妻熟女乱码| 国产精品久久久久久精品古装| 久久久久久久久久久久大奶| 欧美日韩av久久| 熟女电影av网| 亚洲精品国产av成人精品| 亚洲中文av在线| 在线观看美女被高潮喷水网站| 精品视频人人做人人爽| 亚洲国产欧美在线一区| 日本欧美视频一区| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久国产电影| 26uuu在线亚洲综合色| 亚洲精品在线美女| 一区二区三区精品91| 9热在线视频观看99| 人体艺术视频欧美日本| 性色avwww在线观看| 五月开心婷婷网| 精品少妇久久久久久888优播| 亚洲精品一二三| 日韩av免费高清视频| 91成人精品电影| 国产欧美日韩综合在线一区二区| 纯流量卡能插随身wifi吗| 两性夫妻黄色片| 久久精品夜色国产| 亚洲激情五月婷婷啪啪| 水蜜桃什么品种好| 日韩av不卡免费在线播放| 亚洲精品日本国产第一区| 欧美激情高清一区二区三区 | 一级a爱视频在线免费观看| 男女下面插进去视频免费观看| 人妻人人澡人人爽人人| 一级毛片黄色毛片免费观看视频| 国产在视频线精品| 中文字幕另类日韩欧美亚洲嫩草| 我要看黄色一级片免费的| 夫妻性生交免费视频一级片| 久久精品国产综合久久久| 岛国毛片在线播放| 亚洲精品av麻豆狂野| 高清视频免费观看一区二区| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 国产精品久久久久久av不卡| 亚洲经典国产精华液单| 欧美日韩国产mv在线观看视频| 久久久久久久久久久免费av| 男女边吃奶边做爰视频| 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线| 精品午夜福利在线看| 久久久久久久国产电影| 久久久精品国产亚洲av高清涩受| 十八禁高潮呻吟视频| 少妇人妻 视频| 中文字幕制服av| av免费观看日本| 国产成人欧美| 人成视频在线观看免费观看| 久久精品久久精品一区二区三区| 男人添女人高潮全过程视频| 母亲3免费完整高清在线观看 | 一区二区三区激情视频| 精品酒店卫生间| 免费播放大片免费观看视频在线观看| 日韩中文字幕欧美一区二区 | 亚洲精品aⅴ在线观看| www日本在线高清视频| 免费女性裸体啪啪无遮挡网站| 亚洲欧美一区二区三区久久| 色网站视频免费| 午夜老司机福利剧场| 老汉色av国产亚洲站长工具| 欧美精品av麻豆av| 亚洲,欧美精品.| 久久午夜综合久久蜜桃| 99re6热这里在线精品视频| 观看美女的网站| 十八禁网站网址无遮挡| 国产成人a∨麻豆精品| 久久久久久久久久久免费av| 一本久久精品| 满18在线观看网站| 国产黄色视频一区二区在线观看| 一区二区三区精品91| 天堂俺去俺来也www色官网| 国产极品天堂在线| 国产精品久久久久久精品古装| 国产 一区精品| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品国产精品| 在线观看三级黄色| 久久毛片免费看一区二区三区| 欧美精品av麻豆av| 精品人妻熟女毛片av久久网站| 欧美日韩视频高清一区二区三区二| 91精品三级在线观看| 久久精品久久精品一区二区三区| 伦理电影大哥的女人| 久久国产亚洲av麻豆专区| 国产精品.久久久| 9191精品国产免费久久| 大话2 男鬼变身卡| 亚洲精品第二区| 大香蕉久久成人网| 亚洲精品国产av蜜桃| 日韩精品有码人妻一区| tube8黄色片| 大码成人一级视频| h视频一区二区三区| 欧美少妇被猛烈插入视频| 国产不卡av网站在线观看| 美女xxoo啪啪120秒动态图| 欧美+日韩+精品| 国精品久久久久久国模美| 久久久国产精品麻豆| 亚洲国产最新在线播放| 国产免费福利视频在线观看| 一级a爱视频在线免费观看| 免费人妻精品一区二区三区视频| 在线亚洲精品国产二区图片欧美| av免费在线看不卡| 国产精品香港三级国产av潘金莲 | 亚洲精品乱久久久久久| 精品人妻在线不人妻| 成人18禁高潮啪啪吃奶动态图| 99re6热这里在线精品视频| 久久 成人 亚洲| 亚洲经典国产精华液单| www日本在线高清视频| 久久久亚洲精品成人影院| 麻豆乱淫一区二区| 免费在线观看黄色视频的| 极品少妇高潮喷水抽搐| 有码 亚洲区| 欧美激情 高清一区二区三区| 精品少妇久久久久久888优播| 成人免费观看视频高清| 婷婷色综合www| 欧美精品人与动牲交sv欧美| 亚洲成人一二三区av| 视频区图区小说| 一个人免费看片子| 91久久精品国产一区二区三区| 国产成人一区二区在线| 精品少妇内射三级| 国产成人精品无人区| 免费黄色在线免费观看| 少妇被粗大猛烈的视频| av在线老鸭窝| 91久久精品国产一区二区三区| 26uuu在线亚洲综合色| 丰满乱子伦码专区| 波多野结衣av一区二区av| 久久久久精品久久久久真实原创| 日韩 亚洲 欧美在线| 免费女性裸体啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 午夜免费观看性视频| 日韩不卡一区二区三区视频在线| 亚洲成人一二三区av| 一级,二级,三级黄色视频| 热re99久久国产66热| 欧美+日韩+精品| av一本久久久久| 国产精品香港三级国产av潘金莲 | 99国产精品免费福利视频| 97在线视频观看| 美女中出高潮动态图| √禁漫天堂资源中文www| 国产伦理片在线播放av一区| 午夜日韩欧美国产| 男的添女的下面高潮视频| 性色av一级| 亚洲国产av新网站| 久久精品aⅴ一区二区三区四区 | 蜜桃在线观看..| 日本爱情动作片www.在线观看| 亚洲成人av在线免费| 欧美日韩av久久| 亚洲色图 男人天堂 中文字幕| 97在线人人人人妻| 丰满乱子伦码专区| 成人免费观看视频高清| 亚洲国产色片| 五月伊人婷婷丁香| 精品亚洲成a人片在线观看| 亚洲欧洲日产国产| 黄色一级大片看看| a级毛片黄视频| 如日韩欧美国产精品一区二区三区| 欧美成人午夜免费资源| 午夜久久久在线观看| 精品人妻熟女毛片av久久网站| av电影中文网址| 视频在线观看一区二区三区| 狠狠精品人妻久久久久久综合| xxx大片免费视频| 亚洲男人天堂网一区| videossex国产| a级毛片在线看网站| 女人被躁到高潮嗷嗷叫费观| av又黄又爽大尺度在线免费看| 欧美激情 高清一区二区三区| 好男人视频免费观看在线| 美女大奶头黄色视频| 欧美日韩亚洲高清精品| 欧美精品国产亚洲| 国产在线视频一区二区| 色婷婷久久久亚洲欧美| 亚洲熟女精品中文字幕| 免费少妇av软件| 国产成人av激情在线播放| 黑丝袜美女国产一区| 亚洲欧美中文字幕日韩二区| av福利片在线| 欧美在线黄色| 黄色怎么调成土黄色| 日本免费在线观看一区| 天堂俺去俺来也www色官网| 日本vs欧美在线观看视频| 亚洲人成网站在线观看播放| 亚洲成人手机| 亚洲av中文av极速乱| 黑人巨大精品欧美一区二区蜜桃| 免费av中文字幕在线| 性高湖久久久久久久久免费观看| 亚洲一区中文字幕在线| 哪个播放器可以免费观看大片| 久久久国产欧美日韩av| 国产综合精华液| 狠狠婷婷综合久久久久久88av| 99热全是精品| 一级毛片 在线播放| 国产午夜精品一二区理论片| 亚洲婷婷狠狠爱综合网| 丝袜人妻中文字幕| 亚洲少妇的诱惑av| 亚洲国产精品一区三区| 国产有黄有色有爽视频| 久久久久久免费高清国产稀缺| 桃花免费在线播放| 人人妻人人爽人人添夜夜欢视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲成人手机| 亚洲av国产av综合av卡| 新久久久久国产一级毛片| 制服人妻中文乱码| 日本免费在线观看一区| 婷婷成人精品国产| 2018国产大陆天天弄谢| 女性生殖器流出的白浆| 久久久久国产精品人妻一区二区| 国产熟女午夜一区二区三区| 国产欧美日韩一区二区三区在线| 亚洲国产av新网站| 欧美av亚洲av综合av国产av | 热re99久久国产66热| 男女边摸边吃奶| 亚洲国产精品国产精品| 伦理电影免费视频| 色网站视频免费| 国产亚洲一区二区精品| 久久久久国产一级毛片高清牌| 日本黄色日本黄色录像| 午夜激情av网站| 国产一区二区在线观看av| 18禁国产床啪视频网站| av电影中文网址| 波多野结衣一区麻豆| 国产淫语在线视频| 精品卡一卡二卡四卡免费| 国产极品天堂在线| 亚洲欧美色中文字幕在线| 亚洲 欧美一区二区三区| 久久人妻熟女aⅴ| 18禁动态无遮挡网站| 国产男女内射视频| 亚洲成色77777| 国产精品国产三级国产专区5o| 亚洲精品乱久久久久久| 久久鲁丝午夜福利片| 国产成人aa在线观看| 亚洲综合色惰| 亚洲在久久综合| 中文乱码字字幕精品一区二区三区| 免费播放大片免费观看视频在线观看| 久久毛片免费看一区二区三区| 亚洲综合精品二区| 婷婷色综合大香蕉| av视频免费观看在线观看| av电影中文网址| 欧美人与性动交α欧美精品济南到 | 边亲边吃奶的免费视频| 日韩制服丝袜自拍偷拍| 亚洲av电影在线观看一区二区三区| 日本-黄色视频高清免费观看| 中文字幕制服av| 国产人伦9x9x在线观看 | 黄网站色视频无遮挡免费观看| 日韩三级伦理在线观看| 亚洲中文av在线| 三级国产精品片| 在线观看三级黄色| 国产精品一区二区在线观看99| 久久精品久久久久久久性| 一边亲一边摸免费视频| 老司机影院毛片| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 熟女电影av网| 中文字幕制服av| 青春草亚洲视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 女的被弄到高潮叫床怎么办| 亚洲欧洲精品一区二区精品久久久 | 美女高潮到喷水免费观看| 国产日韩欧美视频二区| av网站在线播放免费| 欧美日韩综合久久久久久| 如何舔出高潮|