• <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
    国产成人系列免费观看| 欧美国产精品一级二级三级| 国产在视频线精品| 免费看av在线观看网站| 国产伦理片在线播放av一区| 久久性视频一级片| 久久午夜综合久久蜜桃| 啦啦啦啦在线视频资源| 十八禁网站网址无遮挡| 久久99精品国语久久久| 中国国产av一级| 黄色视频在线播放观看不卡| 国产精品久久久人人做人人爽| 性色av乱码一区二区三区2| 两性夫妻黄色片| av天堂久久9| 精品福利观看| 国产91精品成人一区二区三区 | 亚洲精品一区蜜桃| 亚洲av日韩精品久久久久久密 | 亚洲欧美精品综合一区二区三区| 成人午夜精彩视频在线观看| 日韩大码丰满熟妇| 久久影院123| 国产黄色视频一区二区在线观看| av欧美777| 天堂中文最新版在线下载| e午夜精品久久久久久久| 国产高清videossex| 国产片特级美女逼逼视频| 新久久久久国产一级毛片| 国产成人精品在线电影| 少妇的丰满在线观看| 99九九在线精品视频| www.自偷自拍.com| 肉色欧美久久久久久久蜜桃| 又粗又硬又长又爽又黄的视频| 国产成人精品久久二区二区免费| 精品人妻1区二区| 成人免费观看视频高清| 黄色片一级片一级黄色片| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品古装| 国产91精品成人一区二区三区 | 19禁男女啪啪无遮挡网站| 99精国产麻豆久久婷婷| 极品人妻少妇av视频| 一区二区三区四区激情视频| 久热爱精品视频在线9| 欧美日韩亚洲高清精品| 中文字幕人妻丝袜一区二区| 男女之事视频高清在线观看 | 欧美成人精品欧美一级黄| 亚洲久久久国产精品| 青草久久国产| 国产成人免费无遮挡视频| 国产欧美日韩一区二区三区在线| 黄色一级大片看看| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕av电影在线播放| 女人被躁到高潮嗷嗷叫费观| 成人亚洲欧美一区二区av| 国产一区二区在线观看av| 国产免费视频播放在线视频| 午夜免费成人在线视频| 亚洲欧美一区二区三区黑人| 一级毛片电影观看| 久久久精品区二区三区| 久久精品人人爽人人爽视色| 国产精品秋霞免费鲁丝片| 大片免费播放器 马上看| 亚洲 国产 在线| 国产在线免费精品| 精品国产一区二区久久| 女性被躁到高潮视频| 久久ye,这里只有精品| 狂野欧美激情性bbbbbb| 午夜福利一区二区在线看| 女人久久www免费人成看片| 91精品三级在线观看| a级片在线免费高清观看视频| 国产真人三级小视频在线观看| 精品少妇黑人巨大在线播放| 母亲3免费完整高清在线观看| 久久久精品国产亚洲av高清涩受| 亚洲国产av影院在线观看| 午夜老司机福利片| 久久久久视频综合| 色综合欧美亚洲国产小说| svipshipincom国产片| 久久人人爽av亚洲精品天堂| 国产成人啪精品午夜网站| 成人国产av品久久久| av欧美777| 老司机影院毛片| 久久久精品区二区三区| 国产精品99久久99久久久不卡| 热99国产精品久久久久久7| 狠狠精品人妻久久久久久综合| 亚洲成人国产一区在线观看 | 99九九在线精品视频| 女人被躁到高潮嗷嗷叫费观| 亚洲av美国av| 国产成人精品久久久久久| 18禁黄网站禁片午夜丰满| 久久精品久久精品一区二区三区| 亚洲三区欧美一区| h视频一区二区三区| 久久鲁丝午夜福利片| 男女之事视频高清在线观看 | 欧美乱码精品一区二区三区| 欧美成人精品欧美一级黄| 亚洲人成网站在线观看播放| 一级毛片黄色毛片免费观看视频| 亚洲国产欧美一区二区综合| 久9热在线精品视频| 两性夫妻黄色片| 可以免费在线观看a视频的电影网站| 欧美国产精品一级二级三级| 亚洲精品久久午夜乱码| www日本在线高清视频| 日韩一区二区三区影片| 亚洲成人免费av在线播放| 91精品三级在线观看| 国产成人精品久久久久久| 亚洲色图 男人天堂 中文字幕| 国产男女内射视频| 在现免费观看毛片| 欧美激情 高清一区二区三区| 成年人黄色毛片网站| 精品少妇久久久久久888优播| 好男人电影高清在线观看| 天天添夜夜摸| 中文字幕最新亚洲高清| 成人国产一区最新在线观看 | 日韩 亚洲 欧美在线| 亚洲精品一二三| 观看av在线不卡| 亚洲国产欧美网| 国产精品麻豆人妻色哟哟久久| 欧美精品一区二区大全| 国精品久久久久久国模美| 悠悠久久av| 黄色一级大片看看| 嫁个100分男人电影在线观看 | 久久国产精品大桥未久av| 亚洲国产欧美在线一区| 美女视频免费永久观看网站| 亚洲av日韩在线播放| 欧美国产精品一级二级三级| 丰满饥渴人妻一区二区三| 国产亚洲av高清不卡| 亚洲国产欧美在线一区| 两人在一起打扑克的视频| 99热网站在线观看| 国产一区亚洲一区在线观看| 国产精品一区二区在线不卡| 亚洲中文日韩欧美视频| 天堂俺去俺来也www色官网| 好男人电影高清在线观看| 啦啦啦 在线观看视频| 天天操日日干夜夜撸| 久久久久视频综合| 精品卡一卡二卡四卡免费| 国产成人欧美| 亚洲成人国产一区在线观看 | 黑人欧美特级aaaaaa片| 女性被躁到高潮视频| 亚洲av综合色区一区| 国产又色又爽无遮挡免| 日韩av在线免费看完整版不卡| 精品久久蜜臀av无| 精品第一国产精品| 深夜精品福利| 欧美日韩成人在线一区二区| 国产高清videossex| 精品少妇一区二区三区视频日本电影| 日本欧美视频一区| 欧美黑人精品巨大| 在线天堂中文资源库| 欧美精品一区二区免费开放| 国产高清videossex| 美女国产高潮福利片在线看| 午夜福利乱码中文字幕| 欧美激情高清一区二区三区| 99热国产这里只有精品6| 搡老乐熟女国产| 久久精品熟女亚洲av麻豆精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产最新在线播放| 三上悠亚av全集在线观看| 国产成人一区二区在线| 国产黄频视频在线观看| 欧美精品av麻豆av| 狂野欧美激情性bbbbbb| 国产欧美日韩一区二区三 | 欧美精品亚洲一区二区| 99久久99久久久精品蜜桃| 最近手机中文字幕大全| 女人高潮潮喷娇喘18禁视频| 欧美97在线视频| 这个男人来自地球电影免费观看| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 午夜福利乱码中文字幕| 亚洲一码二码三码区别大吗| 久久久久国产一级毛片高清牌| 高清欧美精品videossex| 精品人妻在线不人妻| 国产精品一区二区在线观看99| 电影成人av| 欧美精品av麻豆av| 五月天丁香电影| 一区二区三区激情视频| 在线亚洲精品国产二区图片欧美| 亚洲视频免费观看视频| 国产极品粉嫩免费观看在线| 亚洲中文日韩欧美视频| 老司机靠b影院| 亚洲中文av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美中文综合在线视频| 免费日韩欧美在线观看| 女人高潮潮喷娇喘18禁视频| av有码第一页| 99国产精品一区二区蜜桃av | 精品一区二区三区av网在线观看 | 99精国产麻豆久久婷婷| 18在线观看网站| 欧美中文综合在线视频| 国产精品熟女久久久久浪| 91字幕亚洲| 精品一品国产午夜福利视频| 午夜av观看不卡| 丁香六月欧美| 久久精品国产a三级三级三级| 90打野战视频偷拍视频| 精品熟女少妇八av免费久了| 成人三级做爰电影| 不卡av一区二区三区| 搡老岳熟女国产| 免费观看人在逋| 18禁黄网站禁片午夜丰满| 精品卡一卡二卡四卡免费| 搡老岳熟女国产| 丝袜人妻中文字幕| 亚洲 欧美一区二区三区| 国产高清视频在线播放一区 | 国产精品.久久久| 国产福利在线免费观看视频| 久久久久国产一级毛片高清牌| 久久精品国产综合久久久| 少妇精品久久久久久久| 美女高潮到喷水免费观看| 最新的欧美精品一区二区| 久久久国产一区二区| 黄色毛片三级朝国网站| 又黄又粗又硬又大视频| 亚洲精品一区蜜桃| 大陆偷拍与自拍| 国产麻豆69| 国产成人免费无遮挡视频| 曰老女人黄片| 精品亚洲成a人片在线观看| 国产免费福利视频在线观看| 又粗又硬又长又爽又黄的视频| 热99国产精品久久久久久7| 欧美+亚洲+日韩+国产| 国产男女超爽视频在线观看| 国产成人免费观看mmmm| 亚洲国产欧美网| 国产免费又黄又爽又色| 视频区欧美日本亚洲| 国产熟女欧美一区二区| 精品一区二区三区四区五区乱码 | 18禁观看日本| 成人免费观看视频高清| 国产97色在线日韩免费| 丰满少妇做爰视频| 欧美精品高潮呻吟av久久| 久久精品国产亚洲av高清一级| 性少妇av在线| 中文字幕av电影在线播放| 国产无遮挡羞羞视频在线观看| 国产一区二区三区综合在线观看| 日本一区二区免费在线视频| 少妇猛男粗大的猛烈进出视频| 亚洲精品一区蜜桃| 欧美激情高清一区二区三区| av网站在线播放免费| 在现免费观看毛片| 国产亚洲av高清不卡| 精品亚洲成a人片在线观看| 天堂中文最新版在线下载| 精品久久蜜臀av无| 中文字幕高清在线视频| av线在线观看网站| 岛国毛片在线播放| 午夜老司机福利片| 欧美人与性动交α欧美软件| 国产免费又黄又爽又色| 久久久久久久久久久久大奶| 精品福利永久在线观看| 色综合欧美亚洲国产小说| 婷婷色av中文字幕| 高清黄色对白视频在线免费看| 久久99热这里只频精品6学生| 精品人妻一区二区三区麻豆| 99国产精品免费福利视频| 少妇裸体淫交视频免费看高清 | 亚洲熟女毛片儿| 尾随美女入室| 国产精品一国产av| 亚洲 欧美一区二区三区| 亚洲精品国产av蜜桃| 亚洲av电影在线进入| 各种免费的搞黄视频| 成年人午夜在线观看视频| 久久久精品免费免费高清| 亚洲美女黄色视频免费看| 国产一区二区 视频在线| 亚洲中文字幕日韩| 在线精品无人区一区二区三| 别揉我奶头~嗯~啊~动态视频 | 在线 av 中文字幕| 成年动漫av网址| 亚洲精品日韩在线中文字幕| 日日夜夜操网爽| 男人舔女人的私密视频| 国产av精品麻豆| 爱豆传媒免费全集在线观看| 午夜福利影视在线免费观看| 精品国产一区二区三区四区第35| 麻豆av在线久日| 国产片内射在线| 啦啦啦啦在线视频资源| av有码第一页| 成人国语在线视频| 欧美日韩精品网址| 男女床上黄色一级片免费看| 精品熟女少妇八av免费久了| 午夜91福利影院| 欧美精品高潮呻吟av久久| 老司机影院毛片| av网站在线播放免费| 久久精品久久久久久久性| 欧美人与善性xxx| 久久久精品国产亚洲av高清涩受| 亚洲av片天天在线观看| 久久精品国产a三级三级三级| 精品免费久久久久久久清纯 | 国产精品久久久久成人av| 久久毛片免费看一区二区三区| 久久人妻熟女aⅴ| 99香蕉大伊视频| 国产片内射在线| 无限看片的www在线观看| 国产精品一二三区在线看| 久久久精品免费免费高清| 麻豆国产av国片精品| 国产一区二区 视频在线| 97在线人人人人妻| 欧美在线一区亚洲| 国产日韩一区二区三区精品不卡| 纯流量卡能插随身wifi吗| 香蕉丝袜av| 亚洲av片天天在线观看| 国产有黄有色有爽视频| 五月天丁香电影| 最近最新中文字幕大全免费视频 | 只有这里有精品99| 婷婷成人精品国产| 精品国产乱码久久久久久小说| 久久久久久人人人人人| 99精品久久久久人妻精品| 免费在线观看影片大全网站 | 国产精品国产三级国产专区5o| 日韩av在线免费看完整版不卡| 久久久久国产精品人妻一区二区| 国产黄频视频在线观看| 久久久久久久大尺度免费视频| 观看av在线不卡| 久久久久久久大尺度免费视频| 中文字幕人妻熟女乱码| www.自偷自拍.com| 在线av久久热| 成年动漫av网址| 国产精品欧美亚洲77777| 精品人妻在线不人妻| 久久久精品94久久精品| 亚洲国产精品999| 欧美久久黑人一区二区| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 欧美日本中文国产一区发布| 免费在线观看日本一区| 国产欧美亚洲国产| 亚洲欧美一区二区三区久久| 亚洲国产精品一区二区三区在线| 国产黄色视频一区二区在线观看| 91国产中文字幕| 亚洲精品国产区一区二| 日本一区二区免费在线视频| 国产野战对白在线观看| 又粗又硬又长又爽又黄的视频| 午夜免费观看性视频| 亚洲精品国产色婷婷电影| 亚洲,欧美精品.| 亚洲欧美日韩高清在线视频 | 18禁国产床啪视频网站| 久久久久久久精品精品| 午夜福利免费观看在线| 2021少妇久久久久久久久久久| 国产成人精品无人区| 在线 av 中文字幕| 久久久亚洲精品成人影院| 97精品久久久久久久久久精品| 只有这里有精品99| 国产爽快片一区二区三区| 99热国产这里只有精品6| 飞空精品影院首页| 两人在一起打扑克的视频| 亚洲伊人色综图| 一区二区av电影网| www.自偷自拍.com| xxx大片免费视频| 久热这里只有精品99| 久久热在线av| 大香蕉久久网| 午夜福利影视在线免费观看| 美女大奶头黄色视频| 国产日韩欧美视频二区| 久久精品久久久久久噜噜老黄| 十八禁高潮呻吟视频| 国产人伦9x9x在线观看| 蜜桃在线观看..| 黑人欧美特级aaaaaa片| 国产免费现黄频在线看| 日本欧美国产在线视频| 汤姆久久久久久久影院中文字幕| 蜜桃国产av成人99| 七月丁香在线播放| 一区二区日韩欧美中文字幕| 黄片播放在线免费| 成人亚洲欧美一区二区av| 欧美精品啪啪一区二区三区 | 亚洲国产av新网站| 日韩 欧美 亚洲 中文字幕| 欧美精品亚洲一区二区| 国产精品偷伦视频观看了| 999久久久国产精品视频| 女人精品久久久久毛片| 成人午夜精彩视频在线观看| 脱女人内裤的视频| 免费观看人在逋| 国产黄色视频一区二区在线观看| 麻豆国产av国片精品| 午夜免费成人在线视频| 国产亚洲午夜精品一区二区久久| 国产欧美日韩精品亚洲av| 精品卡一卡二卡四卡免费| 香蕉丝袜av| 欧美精品亚洲一区二区| 观看av在线不卡| 在线观看人妻少妇| 男女高潮啪啪啪动态图| 国产一区二区激情短视频 | 日韩视频在线欧美| 久久人人爽av亚洲精品天堂| 天天添夜夜摸| 一级a爱视频在线免费观看| 亚洲天堂av无毛| 韩国精品一区二区三区| 又粗又硬又长又爽又黄的视频| a级片在线免费高清观看视频| 免费一级毛片在线播放高清视频 | 国产男女超爽视频在线观看| 黄色视频在线播放观看不卡| 久久九九热精品免费| 久久精品国产a三级三级三级| av网站在线播放免费| 亚洲av美国av| 人妻人人澡人人爽人人| 国产黄色免费在线视频| 久久久久久久国产电影| 一级片'在线观看视频| 午夜福利,免费看| 亚洲国产欧美网| 交换朋友夫妻互换小说| 美女视频免费永久观看网站| 久久久久视频综合| 深夜精品福利| 最近手机中文字幕大全| 久久99精品国语久久久| 十八禁高潮呻吟视频| svipshipincom国产片| 国产欧美亚洲国产| 国产片特级美女逼逼视频| 制服诱惑二区| 波多野结衣av一区二区av| 久久ye,这里只有精品| 国产亚洲欧美在线一区二区| 在线 av 中文字幕| 欧美人与性动交α欧美精品济南到| 丝袜喷水一区| 国产精品一区二区在线不卡| 成人亚洲欧美一区二区av| 九色亚洲精品在线播放| 99精品久久久久人妻精品| 国产日韩欧美在线精品| 久久综合国产亚洲精品| 中文字幕人妻丝袜制服| 精品国产超薄肉色丝袜足j| 操美女的视频在线观看| 久久毛片免费看一区二区三区| 少妇猛男粗大的猛烈进出视频| 99国产精品一区二区蜜桃av | 日韩中文字幕视频在线看片| 欧美人与善性xxx| 亚洲欧美清纯卡通| 国产色视频综合| 在线观看人妻少妇| 久久久久久久久免费视频了| 狠狠婷婷综合久久久久久88av| 亚洲欧美激情在线| av片东京热男人的天堂| 国产又色又爽无遮挡免| 咕卡用的链子| 国产精品偷伦视频观看了| 91成人精品电影| 久久午夜综合久久蜜桃| 久久人妻福利社区极品人妻图片 | 亚洲av男天堂| 尾随美女入室| 老司机深夜福利视频在线观看 | 国产av国产精品国产| 大型av网站在线播放| 国产高清视频在线播放一区 | 亚洲欧美一区二区三区国产| 午夜老司机福利片| 免费看不卡的av| 王馨瑶露胸无遮挡在线观看| 中国美女看黄片| 三上悠亚av全集在线观看| www.精华液| 久久久国产精品麻豆| 九色亚洲精品在线播放| 菩萨蛮人人尽说江南好唐韦庄| 免费看不卡的av| 一区二区三区乱码不卡18| 欧美日韩成人在线一区二区| tube8黄色片| 午夜免费鲁丝| 国产成人精品久久二区二区91| 午夜视频精品福利| 欧美国产精品一级二级三级| 久久国产精品大桥未久av| 亚洲熟女毛片儿| 各种免费的搞黄视频| 久久久久精品人妻al黑| 亚洲欧美中文字幕日韩二区| 一区二区av电影网| 欧美老熟妇乱子伦牲交| 免费在线观看视频国产中文字幕亚洲 | 国产成人精品久久久久久| 每晚都被弄得嗷嗷叫到高潮| 午夜久久久在线观看| 波野结衣二区三区在线| 日韩中文字幕欧美一区二区 | 国产欧美日韩一区二区三 | 建设人人有责人人尽责人人享有的| 一本色道久久久久久精品综合| 十八禁网站网址无遮挡| a 毛片基地| 国产麻豆69| 国产伦理片在线播放av一区| 色婷婷久久久亚洲欧美| 日韩伦理黄色片| 免费看不卡的av| 亚洲自偷自拍图片 自拍| 99香蕉大伊视频| 一区福利在线观看| 欧美久久黑人一区二区| 女人久久www免费人成看片| 99热国产这里只有精品6| 欧美性长视频在线观看| 大话2 男鬼变身卡| 国产精品国产三级国产专区5o| 欧美日韩一级在线毛片| 少妇精品久久久久久久| 久久久精品94久久精品| 久久久欧美国产精品| 搡老乐熟女国产| 日韩,欧美,国产一区二区三区| 男男h啪啪无遮挡| 国产熟女欧美一区二区| 欧美日韩av久久| 国产成人影院久久av| 悠悠久久av| 爱豆传媒免费全集在线观看| 国产精品国产av在线观看| 国产av国产精品国产| 免费观看人在逋| 少妇 在线观看| 日本午夜av视频| 久久久久久久国产电影| 极品人妻少妇av视频| 在线观看一区二区三区激情| 久久久国产一区二区| 一本一本久久a久久精品综合妖精| 黄色怎么调成土黄色| 国产成人精品久久二区二区91| 人人妻,人人澡人人爽秒播 | 久久人人97超碰香蕉20202| 久久热在线av|