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

    前緣尖化對飛翼布局飛行器氣動隱身性能影響

    2022-04-27 01:46:18洪葦江李小婷楊道寧
    關(guān)鍵詞:飛翼迎角前緣

    洪葦江,李小婷,薛 暉,楊 樂,楊道寧

    (1.軍委裝備發(fā)展部某中心,北京 100034;2.北京機電工程研究所,北京 100074;3.航天工程大學(xué)宇航科學(xué)與技術(shù)系,北京 101416)

    飛翼布局氣動效率與隱身性能俱佳[1],是下一代作戰(zhàn)飛機主流外形布局設(shè)計方案之一。對于巡航速度設(shè)計點為亞聲速的飛行器,出于和緩前緣流動加速、降低逆壓梯度及在較大的迎角范圍內(nèi)保持附著流動的考慮,翼型一般采用頭部較為飽滿的超臨界翼型或者層流翼型。其中,層流翼型的低頭力矩較超臨界翼型小,更適用于俯仰通道操縱能力相對正常式布局弱的飛翼布局飛行器。由于翼身高度融合,無垂尾、平尾等部件,飛翼布局擁有優(yōu)秀的外形隱身能力?;诖?,圍繞翼身融合體開展精細化修形是進一步提高飛翼布局隱身性能的重要設(shè)計途徑之一。

    目前,美國B?2 轟炸機與X?47B 兩型飛機的機頭附近區(qū)域采取減小前緣半徑、下表面向內(nèi)凹陷等前緣尖化設(shè)計手段提升整機隱身性能,實現(xiàn)了工程化應(yīng)用[2]。近年來,國內(nèi)學(xué)者對飛翼布局飛行器前緣尖化開展了相關(guān)研究。張彬乾等[3]采用數(shù)值模擬仿真手段,針對翼型開展了氣動隱身優(yōu)化設(shè)計,給出了前緣半徑、彎度、厚度等翼型設(shè)計參數(shù)選取原則;張樂[4]提出一種尖前緣飛翼布局構(gòu)型,并圍繞巡航狀態(tài)開展了氣動隱身計算分析。樊華羽等[5]采用多目標粒子群優(yōu)化算法對飛翼布局無人機開展了氣動隱身多目標優(yōu)化設(shè)計。

    飛翼布局飛行器由于取消了平尾,使得氣動焦點靠前,縱向靜穩(wěn)定度較正常式布局低,特別是對于中/大展弦比飛翼構(gòu)型,縱向尺寸較短,俯仰操縱能力設(shè)計裕度較低[6]。因此在設(shè)計飛翼布局飛行器時,除了關(guān)注巡航狀態(tài)下的氣動性能外,設(shè)計人員也非常關(guān)注起降性能。本文以亞聲速飛翼布局飛行器為研究對象,在飛行包線內(nèi)選取起降狀態(tài)和巡航狀態(tài)進行分析。采用CFD 方法開展前緣尖化范圍對飛翼布局飛行器氣動性能影響的研究工作。采用低頻算法對典型電磁波入射頻率下,飛行器電磁散射特征開展研究。重點分析前緣尖化范圍對隱身性能帶來的影響,為此類飛行器氣動隱身綜合優(yōu)化設(shè)計提供參考。

    1 飛翼布局模型

    本文所研究的飛翼布局構(gòu)型與控制參數(shù)如圖1 所示。飛翼布局構(gòu)型翼身融合程度高,飛行器前緣鏡面反射是其頭向角域雷達散射截面(Radar cross section,RCS)重要的貢獻源,而較大的前緣半徑可能會帶來更為顯著的鏡面反射效應(yīng)[7?8]。因此,采用合理適當?shù)耐庑涡扌问侄?,將鏡面反射削減為較弱的尖頂散射和邊緣繞射,是一種較為可行的隱身修形設(shè)計思路。

    圖1 飛翼布局飛行器模型與控制參數(shù)Fig.1 Flying-wing aircraft model and control parameters

    本文以大幅縮減翼型前緣半徑作為尖化方法,鈍前緣翼型采用r1=0.035l1的層流翼型,尖前緣翼型則基于該層流翼型,縮減翼型前緣半徑r2至0.002 4l1(約為6.9%r1)。并以上述兩種翼型為基礎(chǔ),構(gòu)建3 種沿展向尖化范圍不同的飛翼飛行器三維外形,分別是鈍前緣外形、尖前緣A 外形和尖前緣B 外形,以研究前緣尖化對飛翼布局飛行器氣動隱身性能影響,具體構(gòu)建方法如下:

    對于鈍前緣外形,6 個控制剖面a~f翼型基于鈍前緣翼型設(shè)計,具有明顯的鈍頭部特征。對于尖前緣外形,由機身對稱面a向外一定范圍內(nèi),用尖前緣翼型替換鈍前緣翼型。為研究前緣尖化外形的氣動、隱身性能,本文設(shè)置了前緣尖化范圍約占整個內(nèi)翼段展向1/3 的尖前緣A 外形。該外形a~b站位使用尖前緣翼型,b~c為過渡區(qū),c~f站位使用鈍前緣翼型。為研究比對前緣尖化范圍對氣動、隱身性能的影響,本文設(shè)置尖化范圍更大的尖前緣B 外形。該外形a~c站位使用尖前緣翼型,c~d為過渡區(qū),d~f站位使用鈍前緣外形,前緣尖化范圍較尖前緣A外形擴大一倍。上述3 種外形的側(cè)視圖對比如圖2 所示,表1 給出了本文研究飛翼布局的具體參數(shù)。

    表1 飛翼布局飛行器具體尺寸參數(shù)Table 1 Specific size parameters of flying?wing aircraft

    圖2 3 種外形側(cè)視圖Fig.2 Side view of three shapes

    2 計算方法

    2.1 氣動性能計算方法及驗證

    本文的氣動性能數(shù)值模擬計算以RANS 方程作為控制方程,運用有限體積法離散,選取k?ωSST 兩方程湍流模型[9],該模型充分考慮了邊界層內(nèi)的流動,在近壁面求解原始k?ω方程,通過過渡轉(zhuǎn)換方程,逐漸向遠壁面求解k?e方程[10]。相比標準k?ω模型,其擁有更高的流場求解準確度。

    為分析網(wǎng)格數(shù)量對數(shù)值模擬計算結(jié)果的影響,針對尖前緣A 外形,在保證壁面函數(shù)中y+≈1 的情況下,采用H 型拓撲方式,對半模型進行網(wǎng)格劃分,分別構(gòu)建生成粗網(wǎng)格(168 萬個單元)、中網(wǎng)格(262 萬個單元)和細網(wǎng)格(493 萬個單元)3 套數(shù)量不同的結(jié)構(gòu)化網(wǎng)格。圍繞Ma=0.6,H=20 km 狀態(tài)開展數(shù)值模擬仿真。圖3 給出了阻力系數(shù)、升力系數(shù)和俯仰力矩系數(shù)對比。

    圖3 不同網(wǎng)格數(shù)量下全機氣動力/力矩曲線對比Fig.3 Comparison of longitudinal aerodynamic characteristics under different grids conditions

    由圖3 可以看出:中網(wǎng)格和細網(wǎng)格阻力系數(shù)基本一致,粗網(wǎng)格阻力系數(shù)略?。簧ο禂?shù)三者基本一致;網(wǎng)格數(shù)量對俯仰力矩系數(shù)斜率略有影響,網(wǎng)格量越大,隨著迎角增大,斜率略有降低,但總體上差異不大。為兼顧計算精度和效率,本文采用中網(wǎng)格,如圖4 所示。

    圖4 數(shù)值計算網(wǎng)格Fig.4 Grid for numerical simulation

    為驗證本文所選用數(shù)值模擬方法的可靠性,選取DLR?F6 翼身組合體[11]作為驗證算例。計算狀態(tài)為Ma=0.75、α=0.49°、Re=3.0×106,劃分生成半模結(jié)構(gòu)化網(wǎng)格,網(wǎng)格數(shù)量689 萬個(y+≈1),翼身組合體表面網(wǎng)格分布情況如圖5 所示。

    圖5 校驗計算所用網(wǎng)格Fig.5 Grid for checking computation

    圖6 給出了DLR?F6 翼身組合體4 個典型位置截面示意,y為沿展向坐標,b為全展長。圖7 為各截面壓力系數(shù)Cp計算結(jié)果和試驗值[12]的對比。由圖7 可以看出,計算值與試驗值吻合良好,激波位置與強度捕捉較準確,所采用的CFD 方法可以較有效地反映流場特性,說明了本文采用的CFD 方法和流場求解器合理,滿足研究所需。

    式中:為開口管樁在土塞不完全閉塞條件下的單樁承載力;為開口管樁在土塞完全閉塞條件下的單樁承載力;為樁外側(cè)摩阻力;為樁內(nèi)側(cè)摩阻力;為樁環(huán)底端阻力;為樁土塞阻力;為土塞自重。

    圖6 展向典型截面位置示意圖Fig.6 Arrange of typical section in spanwise direction

    圖7 DLR-F6 翼身組合體壓力系數(shù)分布與試驗數(shù)據(jù)[12]對比Fig.7 Pressure coefficient distribution comparison between computational results and experimental results[12] of DLR?F6 wing?body configuration

    2.2 隱身性能計算方法與校驗

    對于飛行器類目標,當前主流計算方法有多層快速多極子算法(Multi?level fast multiple method,MLFMM)、物理光學(xué)法(Physical optics,PO)以及矩量法。本文研究對象為電大尺寸目標,考慮前緣尖化部位尺寸較小,對于計算精確度有著較高的需求。本文電磁散射仿真采用MLFMM 算法。該算法是一種基于矩量法的快速算法,對于復(fù)雜電大目標計算擁有較高精確度,同時兼顧了計算效率。

    為驗證MLFMM 算法的有效性,選取圓錐為驗證對象。圓錐錐底半徑80.3 mm,錐角22°36',貼體生成三角形表面網(wǎng)格,網(wǎng)格單元最大尺寸為λ/10。雷達工作頻率6 GHz,水平極化。圖8 給出了角域0°~180°狀態(tài)下的計算與試驗值[13]對比??梢钥闯?,計算值與試驗結(jié)果具有較好的重合度,雖然計算結(jié)果的波谷點出現(xiàn)較少,但是對于RCS 曲線趨勢預(yù)測正確,針對波峰的幅值和出現(xiàn)方位角作出了較為精確的預(yù)測。說明本文采用的計算方法對于具有尖銳頭部特征的金屬體RCS 有效。

    圖8 圓錐體RCS 試驗值[13]與MLFMM 算法計算對比Fig.8 Comparison between experimental results[13] and com?putational results of MLFMM of cone model

    本文研究的飛翼布局飛行器模型俯仰角與滾轉(zhuǎn)角均設(shè)定為0°,計算方位角范圍為0°~360°,雷達波入射頻率為2,4 GHz。視模型為金屬目標,生成三角形表面網(wǎng)格,表面網(wǎng)格最大尺寸為λ/10,為確保前緣尖化附近部位的網(wǎng)格適應(yīng)性,對前后緣進行適當加密,模型的網(wǎng)格總量為20 萬個。重點分析±30°附近角域內(nèi)的RCS 變化情況。通過對比鈍前緣外形和兩種尖化范圍不同的尖前緣外形的電磁散射特征,研究分析前緣尖化對飛翼布局飛行器隱身性能的影響。

    3 氣動特性計算結(jié)果與分析

    3.1 起降狀態(tài)

    起降狀態(tài)下計算條件為:Ma=0.2、高度H=0 km,溫度、壓強等大氣參數(shù)由標準大氣表查表得到,計算迎角范圍為-2°~16°。

    圖9 給出了起降狀態(tài)下全機氣動力/力矩曲線對比。三外形阻力系數(shù)在中小迎角狀態(tài)下相差不大,在大迎角區(qū)域略有差異。升力方面,三外形在迎角為-2°~10°范圍內(nèi)均呈現(xiàn)線性增長,升力系數(shù)曲線斜率基本相當。迎角大于10°后非線性增長,鈍前緣外形失速迎角在14°附近,尖前緣B 外形升力系數(shù)拐點出現(xiàn)較早,失速迎角提前至12°附近。對于尖化范圍小于B 外形的A 外形,失速迎角和鈍前緣外形基本相當,但迎角進一步增加后升力系數(shù)下降更為迅速。俯仰力矩方面,在線性段,尖前緣外形的俯仰力矩系數(shù)曲線斜率較鈍前緣外形更為陡峭,氣動焦點無量綱位置相對鈍前緣外形略有后移。進入非線性段后,隨著迎角進一步增加,三外形的俯仰力矩系數(shù)曲線均表現(xiàn)出明顯的上揚趨勢,尖前緣外形在氣動非線性區(qū)域的俯仰力矩較鈍前緣外形更為和緩,在失速特性上有一定的改善。

    圖9 起降狀態(tài)下全機氣動力/力矩曲線對比(Ma=0.2, H=0 km)Fig.9 Comparison of longitudinal aerodynamic characteristics under takeoff and landing condition(Ma=0.2, H=0 km)

    大迎角狀態(tài)下,鈍前緣外形與尖前緣A 外形表面極限流線曲線如圖10 所示。α=10°時,由于機翼后掠效應(yīng),原始外形與尖前緣A 外形的機翼后緣部位展向流動顯著。隨著附面層向翼梢堆積,原始外形翼梢出現(xiàn)分離渦,而尖前緣A 外形翼梢僅前緣部分出現(xiàn)氣流附著差的情況,說明原始外形機翼展向流動更為劇烈。表現(xiàn)在全機俯仰力矩上,原始外形較尖前緣A 外形的俯仰力矩系數(shù)曲線上揚更為顯著,縱向靜穩(wěn)定度急劇變小,甚至反號。隨著迎角進一步增大至16°時,兩外形的上表面均出現(xiàn)較大面積的分離渦以及回流區(qū)域。由于尖前緣A 外形頭部附近前緣半徑非常小,氣流經(jīng)過頭部后附著能力差,產(chǎn)生較為明顯的展向流動趨勢,引起內(nèi)翼段后緣產(chǎn)生較大范圍的分離渦,隨著迎角增大,分離渦區(qū)域進一步增大,外翼段出現(xiàn)了較為嚴重回流以及分離渦。上述情形綜合作用下,兩外形的氣動性能進一步惡化,俯仰力矩曲線上揚,呈現(xiàn)縱向靜不穩(wěn)定。

    圖10 大迎角狀態(tài)鈍前緣與尖前緣A 外形上表面極限流線曲線對比Fig.10 Upper surface limit streamline comparison between shape with blunt leading edge and shape A with sharp leading edge at high angle of attack

    A、B 外形大迎角狀態(tài)下上表面極限流線曲線如圖11 所示。α=12°時,尖前緣B 外形外翼段形成分離渦且后緣出現(xiàn)較大區(qū)域回流,內(nèi)翼段后緣也已形成小范圍分離渦。隨著迎角進一步增加,外翼段和內(nèi)翼段的分離渦區(qū)域逐漸擴大的同時,相互靠近。α=16°時,前緣尖化部位來流無法附著,沿展向流動至內(nèi)外翼轉(zhuǎn)折處,促使原本在內(nèi)外翼的兩個分離渦相遇合并形成一組旋轉(zhuǎn)方向相同的分離渦,使得機翼中段出現(xiàn)大范圍分離區(qū)域。與此相比,尖前緣A 外形由于前緣尖化范圍小,其表面氣流附著情況好于B 外形。該兩外形大迎角狀態(tài)下的流場特征與其宏觀上表現(xiàn)出的氣動性能相符。

    圖11 大迎角狀態(tài)尖前緣A、B 外形上表面極限流線曲線對比Fig.11 Upper surface limit streamline comparison between shapes A and B at high angle of attack

    3.2 巡航飛行

    巡航飛行狀態(tài)下計算條件為:Ma=0.6、高度H=20 km,溫度、壓強等大氣參數(shù)由標準大氣表查表得到,計算迎角范圍為-2°~8°。

    圖12 給出了3 種外形的升阻力、升阻比以及俯仰力矩曲線??梢钥闯?,尖前緣外形的阻力系數(shù)要高于鈍前緣外形,但是隨著迎角增加,阻力系數(shù)量值上的差異逐漸減小。從升力系數(shù)量值上看,三者沒有明顯差別,前緣尖化修形沒有對升力產(chǎn)生顯著影響。升阻比方面,鈍前緣外形較經(jīng)過前緣尖化修形的外形略高,且隨著前緣尖化范圍的擴大,升阻比略有降低。三外形均在α=4°附近達到最大升阻比,鈍前緣外形、尖前緣A 外形和尖前緣B 外形的(CL/CD)max分別為17.5、17.2 和17.0。在相同正迎角下配平本文所研究的飛翼布局飛行器,應(yīng)當向上偏轉(zhuǎn)布置于機翼后緣的副翼產(chǎn)生抬頭力矩,以抵消干凈構(gòu)型本身的低頭力矩,建立起俯仰方向上的力矩平衡。從圖12(d)可以看出,隨著前緣尖化范圍擴大,俯仰力矩系數(shù)偏線逐漸向上平移,意味著干凈構(gòu)型低頭力矩量值的降低。此時副翼上偏更小的角度即可產(chǎn)生足夠的配平力矩??v向配平舵偏角的減小對于全機氣動力特性有兩方面影響。一方面由舵面上偏引起的升力損失變小,另一方面降低了舵面偏轉(zhuǎn)帶來的阻力增量,上述兩因素共同疊加,有利于巡航狀態(tài)下配平升阻比的提升。

    圖12 巡航狀態(tài)下縱向氣動特性對比Fig.12 Comparison of longitudinal aerodynamic characteristics under cruise condition(Ma=0.6,H=20 km)

    圖13 給出了在迎角4°時,不同站位處的壓力系數(shù)分布對比情況。由b站位處壓力分布可知,前緣尖化使得上表面氣流難以附著,導(dǎo)致上表面吸力較鈍前緣外形弱。下表面凹陷起到阻滯氣流作用,使得頭部下表面附近壓力系數(shù)明顯增加,由于遠離重心,產(chǎn)生了一定量的抬頭力矩。對于c站位,尖前緣A 外形已由尖前緣外形變?yōu)殁g前緣外形,因此,其頭部Cp曲線介于鈍前緣外形和尖前緣B 外形之間,解釋了宏觀氣動力矩上,隨著前緣尖化范圍擴大,Cm曲線逐漸向上移動的現(xiàn)象。3 種外形在d站位處均為鈍前緣翼型,Cp曲線基本重合。

    圖13 不同站位下壓力系數(shù)分布對比Fig.13 Pressure coefficient comparison at different stations

    4 隱身特性計算結(jié)果與分析

    對空警戒雷達通常采用主天線產(chǎn)生水平極化波束減少地面雜波干擾[14],本節(jié)開展入射雷達波頻率分別為2,4 GHz,水平極化條件下的尖前緣飛翼布局飛行器單站隱身特性仿真分析(圖14)。

    圖14 RCS 仿真計算結(jié)果對比Fig.14 Comparison of RCS computational results

    從圖14 可以看出:3 種外形的RCS 曲線趨勢較為一致,均在±35°、±90°、±145°附近出現(xiàn)6 個波峰。分析頭向±30°附近范圍內(nèi)RCS,可觀察到,采取前緣尖化的外形,頭向±30°范圍內(nèi)RCS 較鈍前緣外形明顯降低。原因是前緣尖化修形措施有效降低了來自飛行器頭部前緣的鏡面反射貢獻。同時,尖前緣外形在±35°的波峰較鈍前緣略有降低。對于亞聲速飛翼布局隱身飛機,機翼后掠角的設(shè)計與隱身設(shè)計思路高度相關(guān),一般應(yīng)大于頭向威脅角域,以達到將前向電磁波反射波峰“堆積”在前向雷達威脅區(qū)域范圍外的目的。從仿真結(jié)果看,前緣尖化是一種成功的隱身修形手段。

    從表2,3 可看出,在入射頻率2,4 GHz 條件下,尖前緣A、B 外形在±30°、±45°、±360°角域范圍內(nèi)的RCS 均值與峰值小于鈍前緣外形,特別是在前向±30°和±45°角域范圍內(nèi)尤為顯著,進一步說明了前緣尖化在隱身性能提高方面的設(shè)計有效性。

    表2 2 GHz 下RCS 統(tǒng)計對比Table 2 RCS comparison at frequency of 2 GHz dBsm

    表3 4 GHz 下RCS 統(tǒng)計對比Table 3 RCS comparison at frequency of 4 GHz dBsm

    分析前緣尖化范圍對頭向隱身性能的影響。入射波頻率為2 GHz 時,尖前緣A 外形和尖前緣B外形在±45°和±360°角域范圍內(nèi)的RCS 峰值均為4.78 dBsm 和-3.42 dBsm,尖前緣B 外形的RCS峰值較尖前緣A 外形更小。尖前緣A、B 兩外形在±30° 角域范圍內(nèi)的RCS 平均值分別為-29.60 dBsm 和-30.32dBsm,±45°角域范圍內(nèi)的RCS平均值分別為-24.91 dBsm和-25.04 dBsm,360°全向角域范圍內(nèi)的RCS 平均值分別為-23.40 dBsm 和-24.82 dBsm,B 外形略小于A 外形。隨著頻率增大到4 GHz,仍能得到相似的結(jié)論。綜上,前緣尖化范圍的擴大,有利于縮減前向區(qū)域RCS,對隱身性能的提升有一定的幫助。

    5 結(jié)論

    通過對不同前緣尖化范圍的飛翼布局飛行器開展起降、巡航飛行狀態(tài)下的氣動性能仿真分析,結(jié)合典型突防條件下的電磁散射仿真結(jié)果,得出以下主要結(jié)論:

    (1)氣動性能仿真結(jié)果表明,沿展向大范圍尖化前緣,使得失速迎角提前,減小飛行器可用迎角范圍,削弱起降性能。選擇適當?shù)那熬壖饣秶捎行Ы档颓熬壖饣瘜τ陲w翼飛行器起降狀態(tài)下的氣動性能不利影響,和緩氣動非線性段區(qū)域力矩,改善失速特性。對于巡航狀態(tài),隨著前緣尖化范圍的擴大,阻力增大,最大升阻比略有下降。但是下表面“凹陷”帶來的抬頭力矩有利于減小舵面配平帶來的升力損失與阻力增量,有利于配平升阻比的提升。

    (2)隱身性能仿真結(jié)果表明,前緣尖化修形可降低前向角域雷達回波強度,沿展向向外擴大前緣尖化范圍可進一步縮減前向附近范圍內(nèi)RCS,從而有效提升隱身突防性能。

    (3)前緣尖化修形與氣動/隱身性能高度關(guān)聯(lián)。對于亞聲速飛行器,前緣尖化作為一種有效的隱身性能優(yōu)化修形措施,為氣動性能帶來一定的不利影響。因此在飛翼飛行器氣動隱身一體化設(shè)計工作中應(yīng)當充分考慮尖化修形范圍對氣動與隱身性能的影響,追求高外形隱身性能的同時,考慮對氣動性能的影響。統(tǒng)籌考量飛行性能需求和突防能力需求,開展氣動與隱身專業(yè)聯(lián)合設(shè)計,確定適當?shù)那熬壖饣秶?,達到提升飛翼布局飛行器綜合戰(zhàn)技指標的目的。

    在平衡計算精度與計算能力的前提下,深入研究基于氣動/隱身緊耦合優(yōu)化思想的飛行器外形多目標設(shè)計方法,建立一套高效的優(yōu)化流程與算法,是下一步的研究重點,該工作對于提升優(yōu)化設(shè)計效率及縮短總體方案閉合周期具有重要意義。

    猜你喜歡
    飛翼迎角前緣
    連續(xù)變迎角試驗數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    一種飛機尾翼前緣除冰套安裝方式
    飛翼無人機機動飛行非線性魯棒控制方法
    飛翼布局飛機側(cè)風(fēng)起降特性
    飛翼無人機嗡鳴氣動彈性響應(yīng)分析
    飛翼隱身特性數(shù)值模擬
    深水沉積研究進展及前緣問題
    前緣
    小說月刊(2015年2期)2015-04-23 08:49:47
    失速保護系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    久久久久久久大尺度免费视频| 免费黄色在线免费观看| tube8黄色片| 国产一区二区三区综合在线观看| 亚洲精品aⅴ在线观看| 日韩不卡一区二区三区视频在线| 青青草视频在线视频观看| 国产精品久久久久成人av| 亚洲图色成人| 日韩,欧美,国产一区二区三区| 日韩大码丰满熟妇| 成人影院久久| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 叶爱在线成人免费视频播放| 中国三级夫妇交换| 中文字幕精品免费在线观看视频| 国产乱人偷精品视频| 深夜精品福利| 一区在线观看完整版| √禁漫天堂资源中文www| 亚洲欧美一区二区三区久久| 丁香六月天网| 精品人妻在线不人妻| 欧美日韩综合久久久久久| 又大又黄又爽视频免费| 男男h啪啪无遮挡| 老汉色∧v一级毛片| 亚洲欧美成人精品一区二区| 日韩av不卡免费在线播放| 亚洲综合色网址| 欧美人与性动交α欧美软件| 亚洲国产毛片av蜜桃av| 久久99一区二区三区| 日韩大片免费观看网站| 亚洲欧美清纯卡通| 精品国产一区二区三区四区第35| 老鸭窝网址在线观看| 色婷婷av一区二区三区视频| 夜夜骑夜夜射夜夜干| 侵犯人妻中文字幕一二三四区| 亚洲专区中文字幕在线 | 一本一本久久a久久精品综合妖精| 国产又色又爽无遮挡免| 丝袜美足系列| www日本在线高清视频| 国产免费一区二区三区四区乱码| 一本大道久久a久久精品| 男人爽女人下面视频在线观看| 久久精品亚洲熟妇少妇任你| 午夜福利影视在线免费观看| 又黄又粗又硬又大视频| 欧美 日韩 精品 国产| 少妇精品久久久久久久| 中文字幕亚洲精品专区| 考比视频在线观看| 国产av码专区亚洲av| 国产在线免费精品| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美一区二区三区黑人| 免费观看a级毛片全部| 国产黄色免费在线视频| 日日撸夜夜添| 日日爽夜夜爽网站| 欧美成人午夜精品| 午夜免费观看性视频| videosex国产| 美女大奶头黄色视频| 一区在线观看完整版| 成人三级做爰电影| 久久久久久久精品精品| 一级a爱视频在线免费观看| 成年av动漫网址| 精品亚洲成a人片在线观看| 欧美精品高潮呻吟av久久| 国产一卡二卡三卡精品 | 熟女av电影| 亚洲中文av在线| 天堂中文最新版在线下载| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品三级大全| 黑人猛操日本美女一级片| 天堂8中文在线网| 我要看黄色一级片免费的| 999精品在线视频| 国产精品久久久av美女十八| 人人妻人人爽人人添夜夜欢视频| 在线天堂中文资源库| 69精品国产乱码久久久| 日日爽夜夜爽网站| 老熟女久久久| 久久精品国产亚洲av高清一级| 日日撸夜夜添| 久久综合国产亚洲精品| 人体艺术视频欧美日本| 黄色视频不卡| 亚洲av日韩精品久久久久久密 | 乱人伦中国视频| 欧美黄色片欧美黄色片| 中文字幕人妻丝袜制服| 国产一区亚洲一区在线观看| 在线 av 中文字幕| 亚洲,欧美精品.| 成人黄色视频免费在线看| 久久97久久精品| 男人舔女人的私密视频| 男女下面插进去视频免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 777久久人妻少妇嫩草av网站| 另类亚洲欧美激情| 街头女战士在线观看网站| 看免费av毛片| av线在线观看网站| 香蕉丝袜av| 又粗又硬又长又爽又黄的视频| 一边摸一边抽搐一进一出视频| 久热这里只有精品99| 十八禁人妻一区二区| 人人澡人人妻人| 999久久久国产精品视频| 夫妻性生交免费视频一级片| 熟女少妇亚洲综合色aaa.| 久久久久国产精品人妻一区二区| 亚洲第一青青草原| 日韩制服骚丝袜av| 欧美 亚洲 国产 日韩一| 精品一品国产午夜福利视频| 美女主播在线视频| 午夜91福利影院| 日本欧美视频一区| 黄片无遮挡物在线观看| 无限看片的www在线观看| 日韩精品免费视频一区二区三区| 色精品久久人妻99蜜桃| 久久 成人 亚洲| 精品亚洲成国产av| 亚洲精品,欧美精品| 黄色毛片三级朝国网站| 精品一品国产午夜福利视频| 又大又黄又爽视频免费| 一二三四在线观看免费中文在| 欧美黄色片欧美黄色片| 久久久久视频综合| 人体艺术视频欧美日本| 精品一区在线观看国产| 亚洲av综合色区一区| 成人免费观看视频高清| 国产精品亚洲av一区麻豆 | 电影成人av| 久久久久人妻精品一区果冻| 国产日韩欧美亚洲二区| 久久精品亚洲av国产电影网| 夜夜骑夜夜射夜夜干| 国产av一区二区精品久久| 久久性视频一级片| av不卡在线播放| 国产成人午夜福利电影在线观看| 99久久精品国产亚洲精品| www日本在线高清视频| 熟女av电影| 中文字幕人妻熟女乱码| 欧美国产精品一级二级三级| 91精品伊人久久大香线蕉| 最近的中文字幕免费完整| 伊人久久大香线蕉亚洲五| 两个人免费观看高清视频| 亚洲精品aⅴ在线观看| 色婷婷av一区二区三区视频| 国产一区二区 视频在线| 精品少妇内射三级| 国产日韩欧美视频二区| 丝瓜视频免费看黄片| 久久av网站| 日韩人妻精品一区2区三区| 国产男人的电影天堂91| 婷婷色综合www| av在线老鸭窝| 午夜福利免费观看在线| 国产精品免费大片| 国产精品蜜桃在线观看| 一级,二级,三级黄色视频| 国产av国产精品国产| 精品亚洲成a人片在线观看| 亚洲精品自拍成人| 一本久久精品| 考比视频在线观看| 亚洲国产欧美一区二区综合| 悠悠久久av| 在线观看一区二区三区激情| 国产一区二区三区综合在线观看| 色精品久久人妻99蜜桃| 日韩av免费高清视频| 日日爽夜夜爽网站| 亚洲欧美激情在线| 亚洲欧美色中文字幕在线| 免费看不卡的av| 亚洲av日韩精品久久久久久密 | 亚洲av在线观看美女高潮| 人妻一区二区av| 十八禁人妻一区二区| 精品卡一卡二卡四卡免费| 日韩一区二区视频免费看| 69精品国产乱码久久久| 少妇人妻精品综合一区二区| 中文字幕另类日韩欧美亚洲嫩草| 天美传媒精品一区二区| 欧美97在线视频| 亚洲色图 男人天堂 中文字幕| 无限看片的www在线观看| 亚洲,欧美精品.| 色吧在线观看| 久久久久久久大尺度免费视频| 亚洲欧美精品自产自拍| 亚洲 欧美一区二区三区| 国产精品久久久av美女十八| 视频区图区小说| 亚洲欧美中文字幕日韩二区| 国产成人精品福利久久| 成人漫画全彩无遮挡| 欧美亚洲日本最大视频资源| 不卡视频在线观看欧美| 纯流量卡能插随身wifi吗| 国产爽快片一区二区三区| 日韩一区二区视频免费看| 人体艺术视频欧美日本| 99精品久久久久人妻精品| 国产亚洲欧美精品永久| 一级毛片我不卡| 999久久久国产精品视频| 国产日韩欧美在线精品| 男女免费视频国产| 久久久久国产一级毛片高清牌| 热re99久久国产66热| 涩涩av久久男人的天堂| 国产不卡av网站在线观看| 国产av码专区亚洲av| 婷婷色综合大香蕉| 91老司机精品| 国产精品久久久久久久久免| 久久鲁丝午夜福利片| 精品国产乱码久久久久久小说| 黄色怎么调成土黄色| 亚洲少妇的诱惑av| 最黄视频免费看| 欧美日韩精品网址| 国产免费福利视频在线观看| 久久精品国产综合久久久| 亚洲第一青青草原| 丰满乱子伦码专区| 欧美在线黄色| av片东京热男人的天堂| 亚洲精品国产区一区二| 免费高清在线观看视频在线观看| 欧美 亚洲 国产 日韩一| 19禁男女啪啪无遮挡网站| www.精华液| 国产在线一区二区三区精| 91国产中文字幕| 一边摸一边抽搐一进一出视频| 国产免费现黄频在线看| 黑人猛操日本美女一级片| 国产亚洲精品第一综合不卡| 久久久精品国产亚洲av高清涩受| 18禁动态无遮挡网站| 亚洲激情五月婷婷啪啪| 色婷婷久久久亚洲欧美| 国产免费现黄频在线看| 日日爽夜夜爽网站| 久久天堂一区二区三区四区| 男女边吃奶边做爰视频| 国产精品麻豆人妻色哟哟久久| 日本色播在线视频| 另类精品久久| 国产精品偷伦视频观看了| 成年人午夜在线观看视频| 90打野战视频偷拍视频| 亚洲色图综合在线观看| 亚洲国产精品999| 色综合欧美亚洲国产小说| 伊人久久大香线蕉亚洲五| 亚洲精品,欧美精品| 亚洲一区中文字幕在线| 波野结衣二区三区在线| 成人三级做爰电影| 亚洲av电影在线进入| 久久精品熟女亚洲av麻豆精品| 亚洲成人av在线免费| 七月丁香在线播放| 99九九在线精品视频| 国产日韩欧美亚洲二区| 天天操日日干夜夜撸| 精品第一国产精品| 悠悠久久av| 精品国产露脸久久av麻豆| 国产精品人妻久久久影院| 韩国高清视频一区二区三区| 亚洲人成电影观看| 国产乱来视频区| 永久免费av网站大全| 咕卡用的链子| 天天添夜夜摸| 中文字幕高清在线视频| 黄色视频不卡| 一区二区三区四区激情视频| 久久婷婷青草| 国产精品免费大片| 在线观看免费视频网站a站| 中国三级夫妇交换| 午夜福利视频在线观看免费| 桃花免费在线播放| 精品亚洲乱码少妇综合久久| 久久这里只有精品19| 如何舔出高潮| 别揉我奶头~嗯~啊~动态视频 | 国产成人啪精品午夜网站| 涩涩av久久男人的天堂| e午夜精品久久久久久久| 中文字幕亚洲精品专区| 看非洲黑人一级黄片| 精品国产超薄肉色丝袜足j| 天天躁夜夜躁狠狠久久av| 热99国产精品久久久久久7| 最近中文字幕2019免费版| 一区二区三区四区激情视频| 黄色一级大片看看| 亚洲,一卡二卡三卡| 免费在线观看完整版高清| 国产精品欧美亚洲77777| 午夜影院在线不卡| 色综合欧美亚洲国产小说| 丁香六月天网| 爱豆传媒免费全集在线观看| 久久久国产欧美日韩av| 国产97色在线日韩免费| 欧美精品高潮呻吟av久久| e午夜精品久久久久久久| 国产一级毛片在线| 亚洲欧美一区二区三区黑人| 一级毛片 在线播放| 亚洲人成电影观看| 国产亚洲午夜精品一区二区久久| 丝袜脚勾引网站| 欧美黑人欧美精品刺激| 99久久精品国产亚洲精品| 老司机在亚洲福利影院| 久久鲁丝午夜福利片| 狂野欧美激情性xxxx| 最近的中文字幕免费完整| 国产精品一区二区在线观看99| 午夜激情久久久久久久| 两个人看的免费小视频| 搡老乐熟女国产| 久热爱精品视频在线9| 一本一本久久a久久精品综合妖精| 少妇人妻 视频| 国产精品嫩草影院av在线观看| 成年人免费黄色播放视频| 国产片内射在线| 欧美在线黄色| 亚洲熟女精品中文字幕| 一本大道久久a久久精品| 尾随美女入室| 精品久久蜜臀av无| 一级黄片播放器| 国产欧美亚洲国产| 亚洲av男天堂| 成年人午夜在线观看视频| 国产欧美亚洲国产| 一级片'在线观看视频| 飞空精品影院首页| 亚洲七黄色美女视频| 观看美女的网站| 十八禁高潮呻吟视频| 成人国产麻豆网| 亚洲精品国产av蜜桃| 操出白浆在线播放| 99久久人妻综合| 国产麻豆69| 99久久人妻综合| 极品人妻少妇av视频| 啦啦啦啦在线视频资源| 久久亚洲国产成人精品v| 亚洲av在线观看美女高潮| 女人久久www免费人成看片| √禁漫天堂资源中文www| 久久久久网色| 热99久久久久精品小说推荐| 美女主播在线视频| 少妇人妻精品综合一区二区| 欧美激情高清一区二区三区 | 一边摸一边做爽爽视频免费| 日韩,欧美,国产一区二区三区| 大香蕉久久成人网| 精品国产国语对白av| 另类亚洲欧美激情| 操美女的视频在线观看| 狂野欧美激情性bbbbbb| 男女免费视频国产| 97人妻天天添夜夜摸| 日韩av不卡免费在线播放| 精品少妇一区二区三区视频日本电影 | 国产成人精品福利久久| 久久精品国产综合久久久| 欧美激情高清一区二区三区 | 青草久久国产| 亚洲国产成人一精品久久久| 日本vs欧美在线观看视频| 一个人免费看片子| 国产精品三级大全| 人人妻人人爽人人添夜夜欢视频| 啦啦啦视频在线资源免费观看| 久久女婷五月综合色啪小说| 夜夜骑夜夜射夜夜干| 国产福利在线免费观看视频| 亚洲精品成人av观看孕妇| 免费高清在线观看视频在线观看| 女的被弄到高潮叫床怎么办| 少妇人妻久久综合中文| 青春草视频在线免费观看| 国产精品三级大全| 亚洲,欧美,日韩| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 老司机亚洲免费影院| 一级片免费观看大全| 国产乱来视频区| 日本色播在线视频| 国产乱人偷精品视频| 国产国语露脸激情在线看| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 免费观看性生交大片5| 日韩欧美一区视频在线观看| 成人国产麻豆网| 狂野欧美激情性bbbbbb| 最近中文字幕2019免费版| 亚洲精品久久午夜乱码| 久久狼人影院| 操美女的视频在线观看| 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 国产成人精品福利久久| 最近中文字幕高清免费大全6| 久久久久久久久久久免费av| 99久国产av精品国产电影| av.在线天堂| 一级毛片我不卡| 色94色欧美一区二区| 欧美日韩亚洲高清精品| 日韩电影二区| 国产极品天堂在线| 国产成人a∨麻豆精品| 新久久久久国产一级毛片| 又大又爽又粗| 男人爽女人下面视频在线观看| 在线天堂最新版资源| 男男h啪啪无遮挡| 丝袜美腿诱惑在线| 国产黄色视频一区二区在线观看| 欧美日韩视频精品一区| 精品亚洲成国产av| 国产亚洲最大av| 99re6热这里在线精品视频| 亚洲国产欧美网| 国产福利在线免费观看视频| 亚洲精品在线美女| 性少妇av在线| 成人午夜精彩视频在线观看| 久久青草综合色| 亚洲成人av在线免费| 曰老女人黄片| 亚洲天堂av无毛| av在线app专区| 亚洲欧洲日产国产| 日韩成人av中文字幕在线观看| 精品酒店卫生间| 亚洲自偷自拍图片 自拍| 老鸭窝网址在线观看| 999精品在线视频| av片东京热男人的天堂| 午夜福利乱码中文字幕| 久久精品久久精品一区二区三区| 狠狠婷婷综合久久久久久88av| 国产 精品1| 91aial.com中文字幕在线观看| 亚洲精品一区蜜桃| 男女边吃奶边做爰视频| 欧美激情极品国产一区二区三区| 中文字幕精品免费在线观看视频| 日韩大码丰满熟妇| 亚洲伊人色综图| 大片电影免费在线观看免费| 精品人妻熟女毛片av久久网站| 国语对白做爰xxxⅹ性视频网站| 亚洲av欧美aⅴ国产| 久久久久久免费高清国产稀缺| 亚洲伊人久久精品综合| 欧美变态另类bdsm刘玥| 久久久精品94久久精品| 国产午夜精品一二区理论片| 男女高潮啪啪啪动态图| 一本大道久久a久久精品| 少妇精品久久久久久久| 亚洲欧洲国产日韩| 久久精品久久久久久久性| 精品福利永久在线观看| 欧美国产精品va在线观看不卡| 你懂的网址亚洲精品在线观看| 又黄又粗又硬又大视频| av.在线天堂| av不卡在线播放| 国产成人精品久久久久久| 亚洲四区av| 亚洲七黄色美女视频| 高清欧美精品videossex| 久久精品aⅴ一区二区三区四区| 高清视频免费观看一区二区| 国产一区亚洲一区在线观看| 免费黄频网站在线观看国产| 国产一区二区在线观看av| 久久久久久人人人人人| 黄片无遮挡物在线观看| 天堂俺去俺来也www色官网| 亚洲中文av在线| 美女福利国产在线| videosex国产| 久久精品久久精品一区二区三区| 九草在线视频观看| 亚洲,一卡二卡三卡| 午夜免费男女啪啪视频观看| 少妇人妻久久综合中文| 少妇精品久久久久久久| 欧美日韩视频精品一区| 熟妇人妻不卡中文字幕| 99热网站在线观看| 韩国高清视频一区二区三区| 国产熟女午夜一区二区三区| 不卡av一区二区三区| 美女大奶头黄色视频| 麻豆精品久久久久久蜜桃| a级毛片在线看网站| 在线免费观看不下载黄p国产| 国产成人精品无人区| 97精品久久久久久久久久精品| 一级,二级,三级黄色视频| 97精品久久久久久久久久精品| 国产精品一区二区在线观看99| 欧美日韩国产mv在线观看视频| 久久国产精品大桥未久av| 伦理电影免费视频| 99国产精品免费福利视频| 日本爱情动作片www.在线观看| 一本一本久久a久久精品综合妖精| 一区二区三区四区激情视频| 亚洲在久久综合| 中文字幕av电影在线播放| 欧美人与善性xxx| 亚洲色图综合在线观看| 91精品三级在线观看| 中文字幕高清在线视频| 在线天堂最新版资源| 热re99久久国产66热| 久久av网站| 国产免费视频播放在线视频| 90打野战视频偷拍视频| 一区二区三区四区激情视频| 黄色一级大片看看| 国产成人a∨麻豆精品| 国产一卡二卡三卡精品 | 精品一区在线观看国产| a级片在线免费高清观看视频| 亚洲美女视频黄频| 免费黄频网站在线观看国产| 极品人妻少妇av视频| 久久国产精品男人的天堂亚洲| 久久性视频一级片| 韩国高清视频一区二区三区| 免费在线观看黄色视频的| 亚洲三区欧美一区| 午夜av观看不卡| 丝瓜视频免费看黄片| 如日韩欧美国产精品一区二区三区| 精品少妇一区二区三区视频日本电影 | 人成视频在线观看免费观看| 中文乱码字字幕精品一区二区三区| 满18在线观看网站| 日韩欧美精品免费久久| 久久精品国产亚洲av高清一级| 下体分泌物呈黄色| 亚洲一区二区三区欧美精品| 黄色 视频免费看| 99香蕉大伊视频| 涩涩av久久男人的天堂| 亚洲精品乱久久久久久| 成年动漫av网址| 精品一区二区三区av网在线观看 | 欧美在线黄色| 又粗又硬又长又爽又黄的视频| 国产一卡二卡三卡精品 | 激情视频va一区二区三区| 观看av在线不卡| 国产亚洲精品第一综合不卡| 亚洲国产中文字幕在线视频| 80岁老熟妇乱子伦牲交| 久久综合国产亚洲精品| 亚洲av电影在线进入| 亚洲一区中文字幕在线| 亚洲欧洲国产日韩| 一本一本久久a久久精品综合妖精| 男人爽女人下面视频在线观看| 极品人妻少妇av视频| 精品人妻熟女毛片av久久网站| 一二三四在线观看免费中文在| 国产成人欧美在线观看 |