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

    翼尖形狀對小展弦比擺翼水動力性能的影響

    2022-05-31 09:20:18張生浩閆文輝周軍偉
    船舶力學(xué) 2022年5期
    關(guān)鍵詞:翼尖水翼步長

    張生浩,閆文輝,周軍偉,梅 蕾,于 東

    (1.哈爾濱工業(yè)大學(xué)(威海)海洋工程學(xué)院,山東 威海 264209;2.北方工業(yè)大學(xué)機械與材料工程學(xué)院,北京 100144)

    0 引 言

    談及仿生推進,多數(shù)人會聯(lián)想到那些掌握著神秘的游動機理,擁有超高的推進效率,能在水中來去自如的海洋生物。近年來,越來越多的科研工作者致力于揭開這層神秘機理的“面紗”,不斷解析海洋生物游動中暗藏的物理規(guī)律,并逐步應(yīng)用于工程實踐中。其中,仿魚類波狀游動的推進器設(shè)計制造已成為研究熱點之一,如仿金槍魚月牙形尾鰭的機器人[1-2]、三關(guān)節(jié)仿鯊魚機器人[3]、機電一體化擺動機器魚[4]等。

    魚類綜合利用水動力的機動方式給人們以靈感,人們在觀察和歸納魚類推進作用機理的同時,也在逐步擺脫對魚類生物構(gòu)造的局限性,涌現(xiàn)新的設(shè)計理念,擺翼推進器就是其中之一。Belibassakis等[5-6]研究了將擺動的水翼置于船體下方作為助推裝置,通過與海浪耦合,利用船體的升沉運動產(chǎn)生附加推力,達到節(jié)能的目的。Mao[7]研究了擺翼作為海洋結(jié)構(gòu)物主動式穩(wěn)定鰭,通過調(diào)節(jié)擺動角度與擺動幅度來適應(yīng)周圍環(huán)境變化,能夠較好地維持航行器低航速下的穩(wěn)定性。此外,也有研究者將擺翼設(shè)計成能量收集裝置,Kumar等[8]在擺翼前方設(shè)計了彈性板,使得擺翼可以在海浪中被動地發(fā)生擺動,達到能量收集的目的。

    在擺翼的應(yīng)用中,其運動模型一般是前進、轉(zhuǎn)動以及周期性平動的多自由度耦合運動。為了得到更高的推進效率,Lin 等[9]研究了二維擺翼在純縱搖、純升沉以及耦合運動的水動力性能,發(fā)現(xiàn)純升沉運動比純縱搖擺動的擺翼速度更高,而耦合運動的擺翼效率更高。大多數(shù)研究中,擺翼的運動被認為符合正弦模型,Das 等[10]及劉煥興等[11]研究了非正弦式擺動對擺翼水動力性能的影響,指出不同形狀的運動波形會影響前進速度及效率,在水動力性能大幅降低的高斯特勞哈爾數(shù)工況,使用改進的非正弦擺動可以使尾流區(qū)的旋渦恢復(fù)反卡門渦街的形式,提高擺翼性能。為了建立擺翼性能與流場結(jié)構(gòu)之間的關(guān)系,Andersen 等[12]采用數(shù)值與試驗相結(jié)合的方法研究了擺翼后方的尾跡,用皂膜法獲得了更加清晰的尾流場,更加深入研究了尾流場與推力及阻力的變化規(guī)律。苑宗敬等[13]用三維矩形及橢圓形水翼模仿魚尾鰭的周期性擺動,研究其水動力性能,展示了不同展弦比水翼的尾渦結(jié)構(gòu),直觀地看出三維效應(yīng)對尾渦場的影響。Abbaspour 等[14]在研究振蕩翼時,以獨特的視角顯示了三維水翼兩邊扭曲中間平坦的尾渦層結(jié)構(gòu)形式,尾渦在水翼后方發(fā)展,形成“環(huán)環(huán)相扣”的三維渦街。與二維水翼不同,有限展長的三維水翼在前進的過程中會在每個翼尖處形成一個拖拽的旋渦,如同一個“小型龍卷風(fēng)”,將對周圍流場起強烈的速度誘導(dǎo)作用,在升力下降的同時產(chǎn)生誘導(dǎo)阻力[15],而且翼梢還可能存在空化[16]。Birch等[17]研究指出,翼尖渦強度在翼梢后緣立即達到最大值,并在下游流場的兩倍弦長距離范圍內(nèi),強度幾乎保持恒定。翼尖渦在尾流場與水翼非定常運動甩出的首緣渦相互作用,產(chǎn)生了更加復(fù)雜的尾流場[18],這種相互作用將影響水翼性能。在擺翼推進器非定常運動中,翼尖渦將引起巨大的能量損耗,設(shè)法降低翼尖損失是工程實踐中提高推進器水動力性能的重要途經(jīng)之一。

    前期研究發(fā)現(xiàn),較大展弦比的擺翼推進效率接近理想二維擺翼,但隨著展弦比的減小,效率明顯降低??紤]到實際應(yīng)用中的擺翼在結(jié)構(gòu)上以小展弦比為主,以及參考飛機機翼設(shè)計中對翼尖的處理辦法,本文以改變擺翼的翼尖形狀來降低翼尖損失為出發(fā)點,對不同翼尖形狀的小展弦比擺翼水動力性能進行研究。文中以NACA 0012 翼型為基礎(chǔ),生成了兩組共6 種不同翼尖形狀的水翼模型;而后,從網(wǎng)格拓撲結(jié)構(gòu)和尺寸、求解時間步長、湍流模型及計算收斂特性四個方面對數(shù)值方法進行驗證。最后,對不同翼尖形狀的擺翼性能和流場特征進行對比分析。

    1 水翼模型

    本文采用平直水翼作為研究基準,由NACA 0012翼型沿展向等弦長拉伸得到??紤]到實際作為推進器的擺翼在結(jié)構(gòu)上很難做到非常大的展長,本文針對展弦比AR=2的水翼進行分析,如圖1(a)所示??紤]到平直水翼在端部存在直角,可能引起端面流動分離,因此設(shè)計了倒圓的端面形狀,如圖1(b)所示。Jung等[19]在研究帶端板的三維翼時指出,端板能夠改善尖端渦流的形成,減小誘導(dǎo)阻力,從而使總阻力減小。參考這種思路,本文在水翼兩端設(shè)計了由兩個不同尺寸扁橢球構(gòu)造的端板,其中心位置在翼型最大厚度的弦長處,以水翼端面為橢球的赤道平面,前半個扁橢球在弦向的長半軸為0.034 m,后半個扁橢球在弦向的半軸長為0.074 m,二者的短軸為0.01 m,展向最大厚度為0.008 m(極半徑為0.004 m)。水翼均以弦長c=0.1 m,展弦比AR=2的矩形水翼為參考,水翼外形及投影形狀如圖1所示。

    圖1 不同翼尖結(jié)構(gòu)三維圖及投影圖Fig.1 3D and projection drawings of hydrofoil with different wing tip structures

    在后面的研究中我們發(fā)現(xiàn),翼尖倒圓和端板都沒有很好地抑制翼尖損失,因而又仿照螺旋槳的槳葉,設(shè)計了三種展弦比AR=2 的不同翼尖收縮形式。其中,橢圓形水翼幾何參數(shù)簡單,便于重復(fù)驗證;翼尖加寬的橢圓形用于分析翼尖收縮程度的影響;帶后掠的翼尖形狀是在加寬橢圓形的基礎(chǔ)上仿照側(cè)斜螺旋槳的一個嘗試。三種翼尖收縮形式的水翼如圖2所示。

    圖2 不同翼尖收縮示意圖Fig.2 Schematic diagram of hydrofoil with different wing tip contractions

    2 運動描述

    水翼運動分析采用絕對坐標系xoy,如圖3 所示。水翼三自由度耦合運動分別為勻速前進、橫蕩、艏搖,其中橫蕩運動和艏搖運動方程分別為

    其中旋轉(zhuǎn)參考點為翼型導(dǎo)緣的頂點。

    在前期研究中發(fā)現(xiàn),較大擺幅比有利于提高效率[20],但由于過大的擺幅不易實現(xiàn),本文折衷選擇擺幅比Y0/c=2.5,即橫蕩擺幅為Y0=0.25 m,運動頻率初步選用f=1 Hz??紤]到水翼運動的最大擺角θ0將影響翼尖流動分離,本文對0.3 rad、0.5 rad及0.7 rad進行了對照研究,橫蕩與艏搖運動相位差為φ=π/2。一周期內(nèi)水翼運動、受力情況如圖3所示。

    圖3 水翼運動和受力情況示意圖Fig.3 Schematic diagram of foil motion and force

    3 性能參數(shù)

    對比螺旋槳的進速系數(shù)定義,本文選用橫蕩擺幅B作為特征長度,定義進程hP=TVx與橫蕩擺幅B的比值為進速系數(shù),以J表示

    式中,進程與橫蕩擺幅如圖3所示,T為擺動周期,Vx為前進速度。由上式可知,在擺動頻率f和橫蕩擺幅B一定的條件下,改變進速Vx值即可得到不同的J。本文水翼運動的不同工況是指不同進速系數(shù)的工況。許多文獻中,采用St數(shù)代替進速系數(shù)來描述擺翼的不同工況,兩者關(guān)系表示為

    定義絕對坐標系下,擺翼穩(wěn)定工作后一個周期內(nèi)沿前進方向力的積分均值為平均推力Fˉx,表示為

    參照螺旋槳性能參數(shù)的定義,以水翼掃掠面積BH為參考面積,即橫蕩擺幅B與水翼展長H所圍面積,以fB為參考速度,定義推力系數(shù)KT,表示為

    按相對速度計算,擺翼推進器的工作雷諾數(shù)Re可表示為

    式中,ν為運動粘度系數(shù)。本文各工況雷諾數(shù)Re在1.6 × 105~5.0 × 105范圍內(nèi)。

    定義擺翼穩(wěn)定工作后一個周期上橫蕩運動與艏搖運動瞬時功率的積分均值為平均輸入功率,表示為

    則推進效率可表示為

    同時,與推力類似,將側(cè)向力Fy和擺翼所受轉(zhuǎn)矩Mz無量綱化為

    式中,c為水翼弦長。

    4 數(shù)值方法

    4.1 控制方程

    本文采用彈性變形網(wǎng)格以適應(yīng)擺翼的多自由度耦合運動的模擬要求,其采用的連續(xù)方程和動量方程分別為

    式中,Ω為控制體體積,A→為控制面,V→為流場速度,Vi為V→的分量,V→d為控制體表面運動速度,τij為黏性應(yīng)力張量分量,p為壓強。

    本文采用FINE/Marine 不可壓粘流求解器實現(xiàn)上述方程的求解。由于擺翼工作中伴隨著明顯的分離流,文中對比了常見的Menter[21-22]的k-ω模型和基于SST 模型的DES 模型[23]。其中k-ω模型是工程中應(yīng)用最為廣泛的兩方程渦粘性模型之一,該模型主要求解湍動能k及其比耗散率ω的對流輸運方程,對附著邊界層湍流及適度分離湍流都有著較高的精度。DES模型在1997年由Spalart[24]提出,屬于RANS/LES耦合算法之一,在分離流動的情況下預(yù)測準確性良好。

    4.2 網(wǎng)格劃分及邊界條件

    本文使用NUMECA/HEXPRESS 生成計算網(wǎng)格,仿真模型計算域是一個4 m×2 m×1 m的長方體,圖4為計算域的三個加密區(qū)域,其網(wǎng)格尺寸如表1所示。

    圖4 計算域加密區(qū)域示意圖Fig.4 Schematic diagram of calculation domain refinement area

    表1 網(wǎng)格劃分情況Tab.1 Details of mesh refinement

    翼型表面網(wǎng)格單元大小參考Nordanger K[25]和魏良[26]的工作,初步選定為7 次細分,即翼型表面的網(wǎng)格長度為c/27。由于水翼運動關(guān)于xoy面對稱,因此使用半模計算。

    計算中采用的全六面體非結(jié)構(gòu)網(wǎng)格,能精確捕捉物體的外形,而且計算域絕大部分網(wǎng)格單元都接近于長方體,網(wǎng)格的正交性十分優(yōu)異,翼型表面及其周圍的網(wǎng)格分布如圖5所示。

    圖5 水翼表面及周圍網(wǎng)格分布示意圖Fig.5 Grid distribution of foil surface and surrounding area

    水翼壁面為無滑移壁面條件,水翼前方邊界為遠場,左、右側(cè)及后方邊界均為零壓梯度,上方為描述壓力,下方邊界為鏡像邊界,如圖6所示。

    圖6 邊界條件Fig.6 Boundary conditions

    4.3 數(shù)值可靠性驗證

    為驗證NUMECA 軟件數(shù)值模擬擺翼運動的可靠性,本文驗證了Read[27]和Schouveiler[28]實驗結(jié)果,同樣采用NACA 0012翼型,升沉幅度比為0.75,旋轉(zhuǎn)參考點在距離翼型前緣1/3弦長處,選取最大攻角為20°與25°的兩種工況對比,如圖7所示。

    圖7 數(shù)值驗證Fig.7 Results of numerical verification

    由圖可知,軟件求解結(jié)果與Read[27]和Schouveiler[28]的工作較吻合,本文所研究的水翼攻角與驗證攻角相近,由此來看,數(shù)值方法對擺翼性能的計算結(jié)果較為可靠。

    4.4 網(wǎng)格無關(guān)性驗證

    為了保證網(wǎng)格質(zhì)量不會對計算結(jié)果產(chǎn)生影響,本文進行了大量的驗證。這里列出其中比較關(guān)鍵的、以及較為有意義的一些影響因素的結(jié)果。

    第一輪驗證中,為研究沿弦向分布的水翼壁面最大網(wǎng)格大小對計算結(jié)果的影響,本文以弦長c為基準,通過改變水翼壁面的弦向網(wǎng)格大小,設(shè)計了四種網(wǎng)格細化方案,如表2 所示,驗證工作的艏搖擺角均取為0.3 rad,湍流模型選用k-ω模型,時間步長為0.002 s。

    表2 弦向網(wǎng)格細化方案Tab.2 Scheme of mesh refinement

    計算結(jié)果如圖8 所示,壁面網(wǎng)格大小對計算結(jié)果有較為明顯的影響。方案1 的網(wǎng)格在一些分離較大的工況得到的結(jié)果曲線明顯向右下方發(fā)生了偏移,峰值偏小。方案2 的網(wǎng)格求解結(jié)果較更細密網(wǎng)格的曲線有較小的偏差,部分工況數(shù)值仍然偏小。細分程度更高的方案3及方案4的網(wǎng)格求解結(jié)果趨于重合,可見,方案3的網(wǎng)格已經(jīng)能滿足網(wǎng)格無關(guān)性要求。

    圖8 不同壁面網(wǎng)格大小的結(jié)果性征曲線Fig.8 Performance curves under different wall mesh sizes

    本文分析了壁面網(wǎng)格長度引起性能差異的原因。圖9為不同壁面網(wǎng)格長度算例的首尾端網(wǎng)格節(jié)點分布,顯然,隨著壁面網(wǎng)格數(shù)量增加,翼型細節(jié)逐漸能夠被準確捕捉。

    圖9 首尾端壁面網(wǎng)格分布示意圖Fig.9 Schematic diagram of wall grid distribution at leading edge and trailing edge

    圖10還列出了四種網(wǎng)格長度在J=4時的流場,這是水翼攻角較大的一個工況,所選取的時間為尾渦充分發(fā)展后的某個整周期時刻。從圖10(a)可以看出,采用方案1 的算例中,首緣已經(jīng)出現(xiàn)明顯的分離渦,而細化程度較高的算例結(jié)果中,尾渦依然附著在水翼表面,根據(jù)機翼的工作特點,可以推測,這是導(dǎo)致效率與推力系數(shù)下降的主要原因。

    圖10 J=4工況下某整周期時刻不同網(wǎng)格算例的尾流場Fig.10 Wake of different grid simulations at the whole period under J=4 condition

    在第二輪驗證中,考慮到水翼周圍網(wǎng)格的加密程度會影響尾渦的分辨率,進而可能影響其性能。本文以細化方案3 的網(wǎng)格為基礎(chǔ),在加密區(qū)域II 和III(見圖4),以及翼型附近的區(qū)域設(shè)計了進一步加密的方案(見圖11),組合得到5種加密方案,如表3所示。

    圖11 翼型周圍網(wǎng)格加密示意圖Fig.11 Mesh refinement around foil

    表3 尾流區(qū)域網(wǎng)格劃分情況Tab.3 Details of mesh generation in wake region

    計算結(jié)果如圖12所示,無論怎樣改變水翼尾流區(qū)域的網(wǎng)格大小,求解得到的性能曲線基本重合。相對來講,對水翼周圍網(wǎng)格加密后,性能得到了微幅提升,特別是較高進速系數(shù)工況的效率,提高了約1%~2%。但是整體來看,加密水翼周圍網(wǎng)格,計算結(jié)果并沒有出現(xiàn)大幅度改變。因此,從節(jié)省計算機資源的角度考慮,水翼周圍網(wǎng)格大小重要性低于壁面網(wǎng)格大小,可在資源充沛時適當加密??梢姡谘芯克硇阅苓^程中,尾流區(qū)域的網(wǎng)格大小并不是重點考慮的因素,但尾渦脫落后流場發(fā)展的模擬仍需要足夠的網(wǎng)格精度,要研究尾渦發(fā)展后的形態(tài)時,尾流區(qū)域有必要加密。

    圖12 翼型周圍網(wǎng)格加密算例的效率及推力系數(shù)曲線Fig.12 Performance curves of simulations in wake region with different mesh sizes

    通過本節(jié)研究發(fā)現(xiàn),研究水翼性能時,必須確保水翼壁面網(wǎng)格滿足網(wǎng)格無關(guān)性要求,而尾流區(qū)域的網(wǎng)格則可以適當調(diào)整,以控制網(wǎng)格總數(shù),加快求解速度。后文中將選用方案3 和組合1 的網(wǎng)格進行分析。

    5 結(jié)果與分析

    5.1 湍流模型及時間步長的影響分析

    由于非定常運動的水翼存在大量分離渦,不同湍流模型很可能使得分離渦狀態(tài)不同,進而引起水翼性能求解差異,本文對比研究了k-ω(SST)以及DES兩種常見的分離渦仿真湍流模型的結(jié)果。

    圖13 為不同湍流模型求解的性能曲線,整體來看湍流模型的影響較小,兩種湍流模型得到的曲線基本相同。但在分離較大(進速系數(shù)較小)的工況下,性能出現(xiàn)了微幅差異,湍流模型有略微影響,而在高效運行工況,兩種湍流模型得到的結(jié)果基本一致。

    圖13 不同湍流模型算例的效率及推力系數(shù)曲線Fig.13 Performance curves of different turbulence models

    在較高進速系數(shù)工況,水翼工作攻角總小于失速角,尾渦呈現(xiàn)為“彗星式拖尾”狀,分離較小,湍流模型影響較小。本文選取了分離較大工況下不同時刻的水翼尾渦場進行對比。圖14 為J=4.4 時,水翼在T/8、T/4、T時刻的渦量場。對比圖14(a)發(fā)現(xiàn),兩種湍流模型求解得到的水翼跨中分離渦形狀不同,而翼尖渦形狀差異很小,k-ω模型中,旋渦在壁面附近較為完整,隨即在后方變得斷斷續(xù)續(xù)。而DES 模型在水翼跨中后方產(chǎn)生了多個小渦環(huán),水翼壁面的旋渦也更為復(fù)雜。對比圖14(b)發(fā)現(xiàn),k-ω模型得到的首緣渦為單個脫落渦,而DES模型得到的脫落渦由多個較小的旋渦組成,這種多旋渦脫落將引起水翼壓力發(fā)生波動,由此推測,這很可能就是與k-ω模型求解結(jié)果出現(xiàn)差異的原因。不過,在水翼性能研究中,湍流模型有一定影響,但并沒有引起劇烈的性能變化,差異在可接受范圍內(nèi),因而在隨后性能研究中,本文初步選用k-ω模型。

    圖14 J=4.4時不同時刻水翼的尾渦場Fig.14 Wake with massive separation at different times

    而后,本文研究了求解器時間步長對結(jié)果的影響,設(shè)計了時間步長為0.01 s、0.005 s、0.002 s、0.001 s 的算例進行對照。圖15 為不同時間步長的計算結(jié)果,時間步長對各工況性能均存在較大影響,其中,較大的時間步長(如0.01 s與0.005 s)在非定常求解過程中存在較大加速度的時刻,造成數(shù)值解偏離真解,使得性能曲線形狀發(fā)生改變,特別是圖15(b)中,部分工況區(qū)間的推力系數(shù)曲線呈現(xiàn)明顯凹凸性不一致的情況。雖然0.005 s 時求解得到的效率曲線已經(jīng)基本趨于收斂,但推力系數(shù)曲線仍不能呈現(xiàn)正確的凹凸走勢,時間步長需進一步減小。時間步長達到足夠小(如0.002 s與0.001 s),算例性能曲線基本穩(wěn)定,結(jié)果向真解收斂,本文選取0.002 s的時間步長,已能夠滿足求解精度。

    圖15 不同時間步長算例的性征曲線Fig.15 Performance curves of simulations with different time steps

    水翼性能與渦分離狀態(tài)息息相關(guān),為了研究時間步長對首緣渦分離的影響,本文選取了性能差距較大工況下的尾渦場進行研究。圖16 為J=3.6 時不同時間步長的算例在整周期時刻的尾流場,對比發(fā)現(xiàn),時間步更小的算例,首緣渦更早發(fā)生分離,而較大時間步長對首緣渦分離較為遲鈍。圖16(a)中首緣渦仍附著在翼型表面,與其他結(jié)果對比,該時刻的首緣渦原本應(yīng)當已經(jīng)分離,而較大時間步長使得分離滯后,造成水翼兩側(cè)仍存在較大壓力差,表現(xiàn)為該時刻仍然保持了較高的推力系數(shù),因而推測這也是圖15(b)中推力系數(shù)曲線發(fā)生明顯凹凸不一致的原因。隨著時間步長減小,旋渦分離差異越來越小,時間步長為0.002 s與0.001 s的流場基本相似,本文合理地選用時間步長0.002 s。

    圖16 J=3.6工況不同時間步長的尾流場Fig.16 Wake field with different time steps under J=3.6 condition

    以下分析中,本文將選用k-ω模型,時間步長選用0.002 s。

    5.2 性能的收斂性分析

    為監(jiān)控求解的性能曲線是否達到收斂,本文定義某時刻的均值效率η0為該時刻在一個周期時窗上的積分均值,表達式為

    某時刻均值推力系數(shù)KT0的表達式與之類似,為

    根據(jù)式(13)及式(14)對求解器產(chǎn)生的瞬時數(shù)據(jù)進行處理,得到不同工況下的性能收斂性曲線,如圖17所示,對比發(fā)現(xiàn),J較小的工況下推力系數(shù)收斂較慢,但求解過程總能在兩個運動周期后達到穩(wěn)定??梢姡硇阅苎芯恐辽傩枰蠼?個以上運動周期,在舍去前兩個周期數(shù)據(jù)的基礎(chǔ)上,進行性能分析。

    圖17 性能收斂特性曲線Fig.17 Convergence characteristic curve of performance

    5.3 端部影響

    根據(jù)前文多次計算結(jié)果,對本文設(shè)計的三種翼尖結(jié)構(gòu)的水翼進行求解。圖18為不同翼尖結(jié)構(gòu)的水翼性能曲線計算結(jié)果,圖18(a)顯示,改變翼尖結(jié)構(gòu)并沒有對效率產(chǎn)生明顯的提高,三種結(jié)構(gòu)在各工況效率差距很小。不過,進速系數(shù)越高的工況,倒圓及端板翼尖呈現(xiàn)出越來越高的效率優(yōu)勢,比平直翼尖高約1%~3%。圖18(b)顯示,端板在啟動階段推力系數(shù)微幅下降,但在高效運行的工況區(qū)間推力系數(shù)又微高于倒圓和平面翼尖的水翼,本文所研究的三種翼尖結(jié)構(gòu)并沒有像預(yù)期那樣大幅改善推進性能。綜合比較發(fā)現(xiàn),雖然倒圓翼尖的水翼推力系數(shù)略微降低,但綜合性能良好,一方面,倒圓水翼在進速系數(shù)較高時仍能保持較高的推進效率;另一方面,倒圓水翼結(jié)構(gòu)復(fù)雜度比端板水翼小,更容易設(shè)計制造。因此,選用倒圓翼尖為佳。

    圖18 不同端部的水翼性征曲線Fig.18 Performance curves of different wing tip shapes

    圖19 是三種水翼峰值效率工況(J=6)側(cè)向力及艏搖力矩系數(shù)時歷曲線。對比發(fā)現(xiàn),三種水翼驅(qū)動力在運動的半周期位置差異最大,此時水翼橫蕩速度及工作攻角均在峰值附近,翼尖形狀對分離的影響最為明顯。倒圓水翼的側(cè)向力及艏搖力矩均處于較低水平。輸入功率的主要貢獻來自于側(cè)向力,端板水翼所需驅(qū)動力最大,雖然端板降低了翼尖損失,但其驅(qū)動力隨之增加,性能提升效果不如倒圓水翼。

    圖19 不同端部的水翼峰值效率工況側(cè)向力及艏搖力矩系數(shù)時歷曲線Fig.19 Time history of lateral force and yaw moment coefficient of foil with different ends under peak efficiency condition

    為了研究不同翼尖形式對水翼表面壓力分布的影響,本文選取了最高效率點工況的水翼。圖20為J=6.0、t=4 s 時刻的跨中翼面壓力分布曲線,對比看出,翼尖結(jié)構(gòu)能夠改變翼面壓力分布,平直翼尖與倒圓翼尖正負壓翼面分布大致相同,因此推力及效率十分相近。不同的是,端板水翼正壓面的分布曲線更為尖銳,正壓力有所減小,而負壓面分布與其他兩種翼尖基本相同。由此推測,這種正壓翼面壓力微弱減小是推力系數(shù)較另外兩種結(jié)構(gòu)略微增大的主要原因。

    圖20 不同翼尖結(jié)構(gòu)的翼面壓力分布曲線Fig.20 Pressure curves of different wing tip structures

    考慮到攻角對端壁繞流流動可能具有較大的影響,本文選用綜合性能較優(yōu)的倒圓水翼做進一步研究,分別計算了0.5 rad、0.7 rad的工況,并與0.3 rad工況進行對照,得到的性征曲線如圖21所示。

    圖21 不同擺角的水翼敞水性征曲線Fig.21 Performance curves of foil with different yaw angles

    圖21(a)顯示,倒圓翼尖水翼還可以通過改變擺角的方式進一步提高效率,但是擺角增大后,導(dǎo)致擺翼推進器高效工作的區(qū)間減小。僅從峰值效率來看,擺角為0.5 rad的效率最高,高效工作的工況區(qū)間為J=3~5。此外,圖21(b)顯示,擺角越大,啟動狀態(tài)的推力系數(shù)越大,但在最高效率點附近的推力系數(shù)減小。因此,這提供了一種推進器啟動過程的設(shè)計思路,即變擺角啟動,以較大擺角啟動,可在較低進速系數(shù)的工況達到高效狀態(tài),然后隨速度增加逐漸減小擺角。

    本文推測擺角不同的水翼峰值效率出現(xiàn)差異的原因可能是由翼尖渦強度不同引起的。圖22 是不同擺角的倒圓水翼在峰值效率點工況不同時刻的尾渦場。對比同時刻在水翼對稱切面內(nèi)的尾渦強度可以發(fā)現(xiàn),擺角越大,水翼產(chǎn)生的翼尖渦強度越低。由圖21 可知,擺角越大,最高效率點工況的進速系數(shù)越低,意味著機翼來流速度降低很可能是翼尖渦強度減小的主要原因,由于翼尖渦能量損失減小,最高效率有所提高。當進一步增大擺角,雖然翼尖渦強度降低,但水翼各工況旋渦分離程度較大,流動的不穩(wěn)定性對翼型水動力特性產(chǎn)生較大影響,最高效率反而降低。

    圖22 不同擺角水翼的尾渦場Fig.22 Wake field of foil with different yaw angles

    為了提供擺翼推進器設(shè)計中所需要的數(shù)據(jù),這里還列舉了不同擺角下的最大側(cè)向力系數(shù)和最大艏搖力矩系數(shù)曲線,如圖23所示。

    圖23 不同擺角水翼在各工況下的最大側(cè)向力及最大艏搖力矩系數(shù)曲線Fig.23 Maximum lateral force and maximum yaw moment coefficient of foil with different yaw angles

    5.4 翼尖收縮形式的影響

    最后,本文對四種翼尖收縮形式的水翼性能進行了研究,設(shè)置擺角為0.5 rad,得到的結(jié)果如圖24所示。

    圖24 不同翼尖收縮的水翼敞水性征曲線Fig.24 Performance curves of foils with different wing tip contractions

    可以看出,標準橢圓形收縮的翼尖再一次提高了效率,圖24(a)中標準橢圓形水翼的峰值效率達到了76.5%,相比矩形水翼在相同工況提高了約2.8%,加寬橢圓形水翼在高效工況區(qū)間效率略低于標準橢圓,后掠形水翼推進性能與矩形水翼相當。圖24(b)顯示,后掠形的水翼有較高的推力系數(shù),但是兩種橢圓水翼推力系數(shù)損失較明顯,特別是起動階段,矩形水翼較其他形狀反而有較大的推力系數(shù)。總之,翼尖收縮接近橢圓,效率將提高,但相同展長下的投影面積減小,會損失推力,工程中可通過調(diào)整水翼尺度來達到指定推力。

    本文推測水翼運行中橢圓形收縮的水翼輸入功率較小、效率較高的原因可能是翼尖渦強度有所降低。圖25為不同翼尖收縮形式水翼處于峰值效率工況(J=4.0)工作時,T、T/8及T/4時刻的尾渦切片圖,對比證實了以上猜測。圖中顯示,各個時刻橢圓形、加寬橢圓形及后掠形收縮的翼尖渦強度依次增大。標準橢圓形收縮的翼尖有最低水平的翼尖渦強度,誘導(dǎo)阻力減小,翼尖渦損耗減小,因而效率提高。從圖25(b)可看出,加寬橢圓水翼的翼尖渦強度略微增大,因此效率略低于標準橢圓水翼。后掠形收縮方式與矩形倒圓水翼的翼尖渦強度相當,因而性能曲線與倒圓水翼幾乎重合。翼尖渦損耗是三維水翼不可避免的問題,通過改變翼尖收縮形式,可以有效降低翼尖渦強度??傮w來看,翼尖按橢圓形收縮的翼尖形狀最好,翼尖收縮越快,翼尖渦強度越大,效率越低。圖26顯示了標準橢圓、加寬橢圓及后掠形水翼的最大側(cè)向力及艏搖力矩系數(shù)曲線。

    圖25 不同翼尖收縮形式的水翼尾渦切片圖Fig.25 Wake vortices of hydrofoils with different wing tip contractions

    圖26 不同收縮形式的水翼最大側(cè)向力及最大艏搖力矩系數(shù)曲線Fig.26 Maximum lateral force and maximum yaw moment coefficient of foil with different wing tip contractions

    6 結(jié) 論

    本文詳細驗證了擺翼的水動力性能數(shù)值模擬方法,分析了6 種不同翼尖形狀小展弦比擺翼(AR=2)的水動力性能,并探討了引起性能差異的原因。經(jīng)過大量數(shù)值計算及分析,得到以下結(jié)論:

    (1)水翼表面弦向網(wǎng)格大小對數(shù)值結(jié)果有顯著影響,當弦向網(wǎng)格長度小于c/256時,與更精密的網(wǎng)格求解結(jié)果基本相同。而在尾流場區(qū)域,網(wǎng)格尺寸對性能求解結(jié)果的影響較小。

    (2)由于擺翼運動中的流場是明顯的非定常過程,求解時間步長會顯著影響水翼的水動力性能求解結(jié)果。此外,從計算收斂性來看,至少需要求解2 個周期以上的流場,求解得到的各項性能參數(shù)才能穩(wěn)定。

    (3)不同翼尖結(jié)構(gòu)能夠改變水翼表面的壓力分布,但對水翼性能改善未能到達預(yù)期,后掠形擺翼也未能提高性能;橢圓形翼尖明顯提高了推進效率,但推力系數(shù)有所降低。倒圓形翼尖處于居中水平,效率較高,推力系數(shù)也較大。

    (4)存在最優(yōu)的艏搖擺角,以實現(xiàn)最高效率。低進速系數(shù)情況下,選擇較大的艏搖擺角有較好的性能,而在高進速時,較小的艏搖擺角有一定的優(yōu)勢。

    (5)本階段僅粗略探索了翼尖形狀對性能的影響,現(xiàn)有結(jié)果還不能直接應(yīng)用于工程實踐,需要更細致的參數(shù)對比分析,而后經(jīng)過試驗驗證,才能實現(xiàn)水翼參數(shù)化設(shè)計。

    猜你喜歡
    翼尖水翼步長
    基于Armijo搜索步長的BFGS與DFP擬牛頓法的比較研究
    中高速條件下不同翼尖小翼的數(shù)值模擬分析
    波浪滑翔機橢圓形后緣水翼動力特性研究
    袖珍水翼突防潛艇的設(shè)計構(gòu)想及運用研究
    基于翼尖渦物理特征的誘導(dǎo)阻力減阻機制實驗研究
    三維扭曲水翼空化現(xiàn)象CFD模擬
    基于流動顯示的翼尖渦不穩(wěn)定頻率測量
    基于 FFD 技術(shù)的民用運輸機翼尖裝置設(shè)計研究
    基于逐維改進的自適應(yīng)步長布谷鳥搜索算法
    湍流進流誘發(fā)的二維水翼振動噪聲特性研究
    亚洲aⅴ乱码一区二区在线播放 | 久99久视频精品免费| 高清av免费在线| 老司机福利观看| 久久天躁狠狠躁夜夜2o2o| 亚洲av成人不卡在线观看播放网| 亚洲 欧美 日韩 在线 免费| 日本欧美视频一区| 黄频高清免费视频| 欧美中文日本在线观看视频| 精品久久久精品久久久| 在线视频色国产色| 99国产精品一区二区蜜桃av| 午夜成年电影在线免费观看| 精品少妇一区二区三区视频日本电影| 一区二区日韩欧美中文字幕| 亚洲黑人精品在线| 国内久久婷婷六月综合欲色啪| 岛国在线观看网站| 日韩成人在线观看一区二区三区| 国产伦人伦偷精品视频| 色综合欧美亚洲国产小说| 女人爽到高潮嗷嗷叫在线视频| www.自偷自拍.com| 亚洲九九香蕉| 色综合欧美亚洲国产小说| 日韩国内少妇激情av| av中文乱码字幕在线| 老汉色av国产亚洲站长工具| 国产熟女午夜一区二区三区| 国产成人精品久久二区二区免费| 最好的美女福利视频网| 久久人妻福利社区极品人妻图片| 亚洲成国产人片在线观看| 美国免费a级毛片| 97人妻天天添夜夜摸| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全免费视频| 亚洲免费av在线视频| 色综合婷婷激情| 一区二区三区激情视频| 亚洲男人天堂网一区| 亚洲中文av在线| 中亚洲国语对白在线视频| a在线观看视频网站| 热99re8久久精品国产| 久久国产亚洲av麻豆专区| 12—13女人毛片做爰片一| 男女床上黄色一级片免费看| 亚洲人成电影免费在线| 麻豆一二三区av精品| 欧美一区二区精品小视频在线| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 黄色a级毛片大全视频| 妹子高潮喷水视频| 亚洲,欧美精品.| 十八禁人妻一区二区| 久久影院123| 国产精品亚洲一级av第二区| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 黄色怎么调成土黄色| 咕卡用的链子| 黄色成人免费大全| 亚洲激情在线av| 在线观看午夜福利视频| 久久热在线av| 亚洲 国产 在线| 久久久久久大精品| 91麻豆av在线| 激情视频va一区二区三区| 亚洲,欧美精品.| 欧美精品一区二区免费开放| 看黄色毛片网站| 波多野结衣高清无吗| 大陆偷拍与自拍| 两人在一起打扑克的视频| 好看av亚洲va欧美ⅴa在| 涩涩av久久男人的天堂| a级毛片在线看网站| 天堂影院成人在线观看| 999久久久国产精品视频| 久久久久久大精品| 高清av免费在线| 亚洲专区国产一区二区| 99热只有精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人免费无遮挡视频| www.自偷自拍.com| 国产精品久久久久久人妻精品电影| 成熟少妇高潮喷水视频| 亚洲成人免费电影在线观看| 欧美黄色片欧美黄色片| 一进一出抽搐动态| 欧美乱码精品一区二区三区| 婷婷六月久久综合丁香| 99香蕉大伊视频| 美女福利国产在线| 亚洲国产精品sss在线观看 | 高潮久久久久久久久久久不卡| 亚洲成人久久性| 亚洲中文av在线| 99久久精品国产亚洲精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲第一av免费看| 国产精品 国内视频| 中文字幕最新亚洲高清| 女人爽到高潮嗷嗷叫在线视频| 村上凉子中文字幕在线| 中国美女看黄片| 91老司机精品| 新久久久久国产一级毛片| 亚洲精品中文字幕一二三四区| 国产精品久久电影中文字幕| 亚洲精品国产区一区二| 香蕉丝袜av| 精品午夜福利视频在线观看一区| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| 如日韩欧美国产精品一区二区三区| 日韩大尺度精品在线看网址 | 亚洲中文av在线| 亚洲免费av在线视频| 亚洲午夜精品一区,二区,三区| 欧美老熟妇乱子伦牲交| 波多野结衣一区麻豆| 丁香六月欧美| 日韩成人在线观看一区二区三区| 最近最新免费中文字幕在线| 亚洲中文日韩欧美视频| 午夜影院日韩av| 人成视频在线观看免费观看| 亚洲欧美激情综合另类| 黄色a级毛片大全视频| 国产成人精品久久二区二区免费| 女性被躁到高潮视频| 一级黄色大片毛片| 日韩精品青青久久久久久| 午夜福利,免费看| 色精品久久人妻99蜜桃| 成年女人毛片免费观看观看9| 久久久久久久久中文| 久久精品91蜜桃| avwww免费| 久久人妻福利社区极品人妻图片| 99久久国产精品久久久| 人人澡人人妻人| 欧美不卡视频在线免费观看 | 国产午夜精品久久久久久| 大香蕉久久成人网| 丰满迷人的少妇在线观看| 精品电影一区二区在线| 香蕉丝袜av| 黄色丝袜av网址大全| 午夜视频精品福利| 亚洲精品美女久久av网站| 国产成人av教育| 亚洲专区字幕在线| 丰满的人妻完整版| 免费不卡黄色视频| 国产精品国产高清国产av| 黄片小视频在线播放| 欧美日韩福利视频一区二区| 这个男人来自地球电影免费观看| 国产精品野战在线观看 | 午夜久久久在线观看| 在线观看www视频免费| 成人免费观看视频高清| 伦理电影免费视频| 波多野结衣高清无吗| 交换朋友夫妻互换小说| 一级毛片高清免费大全| 国产高清激情床上av| 别揉我奶头~嗯~啊~动态视频| 无限看片的www在线观看| 久久精品国产清高在天天线| 欧美中文日本在线观看视频| 午夜视频精品福利| 免费看十八禁软件| 成人永久免费在线观看视频| 国产欧美日韩精品亚洲av| 国产亚洲欧美98| 亚洲精品成人av观看孕妇| 国产三级黄色录像| 精品午夜福利视频在线观看一区| 自线自在国产av| 国产精品久久久久久人妻精品电影| 免费在线观看日本一区| 午夜福利一区二区在线看| 国产极品粉嫩免费观看在线| 国产精品99久久99久久久不卡| 在线天堂中文资源库| 国产成人欧美在线观看| 老熟妇仑乱视频hdxx| 国产黄a三级三级三级人| 亚洲七黄色美女视频| 国产亚洲欧美98| 成人永久免费在线观看视频| 精品人妻1区二区| 叶爱在线成人免费视频播放| 91老司机精品| 男女做爰动态图高潮gif福利片 | 成人三级做爰电影| 人人妻,人人澡人人爽秒播| 1024视频免费在线观看| 色播在线永久视频| 亚洲熟妇中文字幕五十中出 | 在线免费观看的www视频| 97超级碰碰碰精品色视频在线观看| 日韩国内少妇激情av| 欧美黑人欧美精品刺激| 久久草成人影院| 午夜两性在线视频| 亚洲av日韩精品久久久久久密| 黑人猛操日本美女一级片| 男人舔女人的私密视频| 一区二区三区精品91| 久久精品亚洲精品国产色婷小说| 久久久水蜜桃国产精品网| 国产亚洲欧美精品永久| 国产精品影院久久| 久久久国产成人免费| 激情在线观看视频在线高清| 亚洲性夜色夜夜综合| 中文字幕最新亚洲高清| avwww免费| 国产麻豆69| 午夜两性在线视频| av超薄肉色丝袜交足视频| 日日摸夜夜添夜夜添小说| 亚洲第一av免费看| 99久久国产精品久久久| 动漫黄色视频在线观看| 啪啪无遮挡十八禁网站| av电影中文网址| 国产男靠女视频免费网站| 99精国产麻豆久久婷婷| 99精品在免费线老司机午夜| 亚洲片人在线观看| 男人的好看免费观看在线视频 | 18禁裸乳无遮挡免费网站照片 | 桃色一区二区三区在线观看| 又紧又爽又黄一区二区| 在线看a的网站| 欧美最黄视频在线播放免费 | 成年女人毛片免费观看观看9| 丝袜在线中文字幕| 久久久久国内视频| 亚洲一区二区三区不卡视频| 亚洲中文av在线| 色尼玛亚洲综合影院| 亚洲专区中文字幕在线| 国产精品成人在线| 变态另类成人亚洲欧美熟女 | 99久久综合精品五月天人人| 新久久久久国产一级毛片| 在线观看免费日韩欧美大片| 国产乱人伦免费视频| 久久中文字幕一级| 免费观看精品视频网站| 两个人免费观看高清视频| 日本vs欧美在线观看视频| 亚洲 欧美 日韩 在线 免费| 亚洲一区二区三区欧美精品| 欧美日韩精品网址| 久久99一区二区三区| 一级毛片精品| 国产成年人精品一区二区 | 亚洲精品美女久久久久99蜜臀| 国产单亲对白刺激| 涩涩av久久男人的天堂| 国产麻豆69| 999精品在线视频| 母亲3免费完整高清在线观看| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 国产一区在线观看成人免费| 欧美日韩乱码在线| 久久精品国产99精品国产亚洲性色 | 亚洲国产精品sss在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 曰老女人黄片| 久久中文看片网| 真人一进一出gif抽搐免费| 热re99久久国产66热| 乱人伦中国视频| а√天堂www在线а√下载| 精品久久久久久成人av| 精品福利观看| 国产欧美日韩一区二区三| 男人操女人黄网站| 亚洲一区二区三区不卡视频| 国产又色又爽无遮挡免费看| 日韩有码中文字幕| 老熟妇仑乱视频hdxx| 黄片大片在线免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲欧美一区二区三区黑人| 中文字幕精品免费在线观看视频| 国产av精品麻豆| 黄色毛片三级朝国网站| 国产黄a三级三级三级人| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 久久人妻av系列| 久久精品国产亚洲av香蕉五月| 久久久久国产一级毛片高清牌| 99久久人妻综合| 在线观看日韩欧美| 99久久久亚洲精品蜜臀av| 男人的好看免费观看在线视频 | 精品国产亚洲在线| 在线观看免费日韩欧美大片| 亚洲专区国产一区二区| 久久精品91无色码中文字幕| 麻豆成人av在线观看| 亚洲黑人精品在线| 欧美精品亚洲一区二区| 日本a在线网址| 精品久久久精品久久久| 免费高清在线观看日韩| 精品国产国语对白av| 国产亚洲欧美在线一区二区| 日日摸夜夜添夜夜添小说| 91精品三级在线观看| 国产aⅴ精品一区二区三区波| 日韩精品免费视频一区二区三区| 亚洲国产中文字幕在线视频| 久久欧美精品欧美久久欧美| 99热国产这里只有精品6| 不卡一级毛片| 99香蕉大伊视频| 热re99久久精品国产66热6| 在线观看午夜福利视频| 亚洲五月色婷婷综合| 亚洲国产欧美一区二区综合| 变态另类成人亚洲欧美熟女 | 国产片内射在线| 国产不卡一卡二| 中文字幕色久视频| 巨乳人妻的诱惑在线观看| 国产亚洲欧美98| 日本五十路高清| 人人妻人人澡人人看| 午夜精品国产一区二区电影| 精品国产一区二区久久| 亚洲精品在线观看二区| 免费不卡黄色视频| 在线视频色国产色| 国产亚洲欧美精品永久| 黄色女人牲交| 丰满迷人的少妇在线观看| 久久精品国产清高在天天线| 国产精品一区二区免费欧美| 老汉色av国产亚洲站长工具| 久久国产亚洲av麻豆专区| 女生性感内裤真人,穿戴方法视频| 久久青草综合色| 国产欧美日韩精品亚洲av| 黄色a级毛片大全视频| 国产精品二区激情视频| 日日摸夜夜添夜夜添小说| 欧美激情极品国产一区二区三区| 热re99久久国产66热| 久久久久久久午夜电影 | 精品福利永久在线观看| 少妇 在线观看| 一级,二级,三级黄色视频| 精品福利观看| 十八禁人妻一区二区| 午夜两性在线视频| 波多野结衣一区麻豆| 亚洲成人久久性| 99香蕉大伊视频| 纯流量卡能插随身wifi吗| 欧美人与性动交α欧美精品济南到| 黄色 视频免费看| 在线观看日韩欧美| av天堂在线播放| 亚洲免费av在线视频| 久久国产乱子伦精品免费另类| av片东京热男人的天堂| 女人高潮潮喷娇喘18禁视频| 欧美中文日本在线观看视频| 精品一品国产午夜福利视频| 天堂影院成人在线观看| 免费高清视频大片| 午夜日韩欧美国产| 久久人妻福利社区极品人妻图片| 亚洲一区中文字幕在线| 又大又爽又粗| 99久久久亚洲精品蜜臀av| 久99久视频精品免费| 日韩精品青青久久久久久| 操美女的视频在线观看| 美女午夜性视频免费| 人妻久久中文字幕网| 国产主播在线观看一区二区| 精品第一国产精品| 国产一区二区三区综合在线观看| 久久国产亚洲av麻豆专区| 淫秽高清视频在线观看| 高清欧美精品videossex| 最近最新中文字幕大全免费视频| 少妇 在线观看| 黄片大片在线免费观看| 日本免费一区二区三区高清不卡 | 国产亚洲精品综合一区在线观看 | 亚洲精品成人av观看孕妇| 久久久久九九精品影院| 一进一出抽搐gif免费好疼 | 亚洲激情在线av| 精品久久久久久久久久免费视频 | 丰满人妻熟妇乱又伦精品不卡| 精品国产美女av久久久久小说| 丰满迷人的少妇在线观看| 天堂√8在线中文| 丁香六月欧美| 免费观看人在逋| 757午夜福利合集在线观看| 18禁美女被吸乳视频| 99久久99久久久精品蜜桃| 欧美激情极品国产一区二区三区| 久久热在线av| 一区二区日韩欧美中文字幕| 很黄的视频免费| 亚洲人成网站在线播放欧美日韩| 午夜亚洲福利在线播放| 真人做人爱边吃奶动态| 十八禁网站免费在线| 国产精品 欧美亚洲| 久久国产乱子伦精品免费另类| 热99re8久久精品国产| 精品一区二区三区视频在线观看免费 | 欧美中文日本在线观看视频| 国产蜜桃级精品一区二区三区| 搡老乐熟女国产| 亚洲精品美女久久久久99蜜臀| 成人精品一区二区免费| 久久精品国产亚洲av高清一级| av免费在线观看网站| 国产色视频综合| 伊人久久大香线蕉亚洲五| 色婷婷av一区二区三区视频| 午夜a级毛片| 久久午夜综合久久蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品综合一区在线观看 | 日韩欧美在线二视频| 亚洲精品一二三| 人妻丰满熟妇av一区二区三区| 女性生殖器流出的白浆| 免费在线观看视频国产中文字幕亚洲| 中国美女看黄片| 麻豆成人av在线观看| 成人黄色视频免费在线看| 亚洲伊人色综图| 人妻久久中文字幕网| 宅男免费午夜| 在线国产一区二区在线| 亚洲黑人精品在线| 老汉色∧v一级毛片| 日本vs欧美在线观看视频| 国产不卡一卡二| 黑人欧美特级aaaaaa片| 久久久久久大精品| 手机成人av网站| 久久午夜亚洲精品久久| 午夜影院日韩av| svipshipincom国产片| 十八禁网站免费在线| 欧美乱码精品一区二区三区| 麻豆av在线久日| 黄色 视频免费看| 国产区一区二久久| 他把我摸到了高潮在线观看| 美女福利国产在线| av网站免费在线观看视频| 欧美日韩黄片免| av片东京热男人的天堂| 免费av毛片视频| 久久国产亚洲av麻豆专区| 51午夜福利影视在线观看| 亚洲欧美激情在线| 免费看十八禁软件| 999久久久国产精品视频| 久久久久国产一级毛片高清牌| 国产91精品成人一区二区三区| 麻豆一二三区av精品| 亚洲国产精品一区二区三区在线| 久久精品91蜜桃| aaaaa片日本免费| avwww免费| 免费观看精品视频网站| 亚洲av第一区精品v没综合| 亚洲中文字幕日韩| 精品一区二区三区四区五区乱码| 一区二区日韩欧美中文字幕| 一级毛片精品| 99国产精品一区二区蜜桃av| 亚洲成人久久性| 成年女人毛片免费观看观看9| 国产精品免费一区二区三区在线| 丝袜美足系列| 麻豆av在线久日| 伊人久久大香线蕉亚洲五| 熟女少妇亚洲综合色aaa.| 亚洲av成人av| 欧美国产精品va在线观看不卡| 中文欧美无线码| av中文乱码字幕在线| 黄色丝袜av网址大全| 国产精品日韩av在线免费观看 | 日韩 欧美 亚洲 中文字幕| 我的亚洲天堂| 99精国产麻豆久久婷婷| 欧美人与性动交α欧美软件| 91精品三级在线观看| 亚洲一码二码三码区别大吗| 搡老熟女国产l中国老女人| 黄色视频,在线免费观看| 咕卡用的链子| 黑丝袜美女国产一区| 久久精品亚洲精品国产色婷小说| 午夜福利在线免费观看网站| 一级毛片精品| av有码第一页| 麻豆一二三区av精品| 国产亚洲欧美在线一区二区| 亚洲色图综合在线观看| 国产精品日韩av在线免费观看 | 91老司机精品| 在线天堂中文资源库| 中文字幕高清在线视频| 亚洲性夜色夜夜综合| 亚洲自偷自拍图片 自拍| 久久99一区二区三区| 国产乱人伦免费视频| 久久精品成人免费网站| 一边摸一边做爽爽视频免费| 国产精品偷伦视频观看了| 12—13女人毛片做爰片一| 午夜老司机福利片| 久久人人精品亚洲av| 老司机午夜福利在线观看视频| 国产精品香港三级国产av潘金莲| 亚洲国产中文字幕在线视频| 宅男免费午夜| 看免费av毛片| 午夜福利在线免费观看网站| 老司机深夜福利视频在线观看| cao死你这个sao货| 一边摸一边做爽爽视频免费| 99精品久久久久人妻精品| 午夜影院日韩av| 午夜久久久在线观看| 欧美黄色片欧美黄色片| 在线观看免费高清a一片| 成在线人永久免费视频| 日日爽夜夜爽网站| 国产99白浆流出| 国产三级黄色录像| 欧美激情久久久久久爽电影 | 午夜福利在线免费观看网站| 老司机深夜福利视频在线观看| 久久久久久久久久久久大奶| 精品国产亚洲在线| 视频区图区小说| 在线观看一区二区三区激情| 久久久国产一区二区| 两个人免费观看高清视频| 黄色a级毛片大全视频| 国产一区二区三区在线臀色熟女 | 亚洲av第一区精品v没综合| 国产精品99久久99久久久不卡| 成年人黄色毛片网站| 婷婷精品国产亚洲av在线| 美女扒开内裤让男人捅视频| 亚洲男人天堂网一区| 国产成人免费无遮挡视频| 欧美成人性av电影在线观看| 在线观看免费视频日本深夜| 欧美日韩亚洲高清精品| 国产aⅴ精品一区二区三区波| 黄片小视频在线播放| 正在播放国产对白刺激| www.999成人在线观看| 色综合婷婷激情| www.999成人在线观看| 美女福利国产在线| 国内毛片毛片毛片毛片毛片| 一本综合久久免费| 最近最新中文字幕大全免费视频| 亚洲欧美激情在线| 亚洲人成伊人成综合网2020| 午夜福利欧美成人| 国产高清videossex| 日本黄色日本黄色录像| 亚洲国产欧美网| 国产成人精品久久二区二区91| 中文亚洲av片在线观看爽| 波多野结衣一区麻豆| 男女午夜视频在线观看| √禁漫天堂资源中文www| 91精品三级在线观看| 精品久久久精品久久久| 中国美女看黄片| 精品电影一区二区在线| 成年女人毛片免费观看观看9| 精品国产美女av久久久久小说| 黄色毛片三级朝国网站| 免费观看人在逋| 亚洲中文字幕日韩|