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

    基于卡爾曼濾波的翼型表面壓力實時重構(gòu)方法

    2023-05-09 08:42:48劉余丹周楷文劉應(yīng)征
    空氣動力學(xué)學(xué)報 2023年4期
    關(guān)鍵詞:升力卡爾曼濾波流場

    劉余丹,周楷文,劉應(yīng)征,溫 新,*

    (1.上海交通大學(xué) 機械與動力工程學(xué)院,上海 200240;2.空氣動力學(xué)國家重點實驗室,綿陽 621000)

    0 引 言

    飛行器實時控制是實現(xiàn)飛行器姿態(tài)調(diào)整、提高飛行器機動性能的有效方法之一,而飛行器環(huán)境參數(shù)的感知是實現(xiàn)其飛行控制的基礎(chǔ)。飛行器在飛行過程中往往面臨著復(fù)雜的非定常流動問題,但傳統(tǒng)的傳感器設(shè)置方案,如慣性測量裝置、磁力儀、風(fēng)速儀等,受測量技術(shù)、安裝空間等限制,僅能感知到局部的、離散的參數(shù),難以為飛行操縱系統(tǒng)及時提供飛行控制所需的全部信息。在此情況下,難以對飛行器的關(guān)鍵氣動參數(shù)進(jìn)行精準(zhǔn)識別,進(jìn)而難以實現(xiàn)運行監(jiān)測與主動控制。因此,利用有限的觀測數(shù)據(jù)實時感知飛行器周圍流場的全貌具有重要意義。

    許多學(xué)者研究了如何通過有限的測量數(shù)據(jù)重建整個流場的問題。其中,一部分學(xué)者的研究主要基于Everson和Sirovich[1]在1995年提出的一種恢復(fù)缺失數(shù)據(jù)的gappy POD方法展開。gappy POD方法可用于將風(fēng)洞實驗所獲得的有限表面壓力測量數(shù)據(jù)與CFD計算數(shù)據(jù)融合,以感知流場的全貌[2],但gappy POD所采用的求解流場重構(gòu)系數(shù)的最小二乘法對噪聲較為敏感。還有一些學(xué)者基于機器學(xué)習(xí)方法進(jìn)行了研究。Deng等[3]采用長短時間序列(LSTM)神經(jīng)網(wǎng)絡(luò)建立離散測量點數(shù)據(jù)與POD模態(tài)系數(shù)之間的關(guān)系,進(jìn)而實現(xiàn)流場全重構(gòu)。但機器學(xué)習(xí)的方法存在缺乏物理解釋性、過擬合等問題。

    Candes等[4]2006年提出的壓縮感知理論指出,如果一個信號在某一特定基上具有稀疏表示,則該信號可以從相對較少的隨機采樣中得到恢復(fù),該結(jié)論為從有限的測量數(shù)據(jù)實現(xiàn)流場實時重構(gòu)指明了新的方向。此前已有部分學(xué)者嘗試將壓縮感知理論應(yīng)用于流場重構(gòu)。Huang[5]對壓縮感知在航空航天應(yīng)用中的潛力進(jìn)行了探討。Zhao等[6]提出了一種基于壓縮感知重構(gòu)翼型表面壓力的精細(xì)方法,以獲得翼型表面完整的壓力分布信息以及精確的升力和力矩特性。Bright等[7]利用壓縮感知技術(shù),通過圓柱體表面的有限離散壓力測量數(shù)據(jù)重構(gòu)圓柱體周圍的流動特性。Wen等[8]開發(fā)了一種基于壓縮感知的數(shù)據(jù)挖掘方法,從低速圓柱尾流誘導(dǎo)的非定常壓力場的高噪聲測量中對流場數(shù)據(jù)進(jìn)行重構(gòu)。但基于壓縮感知的流場重構(gòu)方法同樣不具魯棒性,易受噪聲影響。

    本文所提出的基于壓縮感知與卡爾曼濾波器的實時流場重構(gòu)方法,采用卡爾曼濾波器,針對基于壓縮感知的流場重構(gòu)方法進(jìn)行改進(jìn)。該重構(gòu)方法具體可分為線上和線下兩個部分:線下全場采樣,通過本征正交分解建立數(shù)據(jù)庫;線上稀疏采樣,通過壓縮感知優(yōu)化重構(gòu)系數(shù),然后采用基于內(nèi)嵌動態(tài)模態(tài)分解的卡爾曼濾波器對所求解時間系數(shù)進(jìn)行實時去噪,提高重構(gòu)方法的魯棒性。

    1 重構(gòu)方法

    基于壓縮感知與卡爾曼濾波器的實時流場重構(gòu)方法流程圖如圖1所示。

    1.1 本征正交分解

    本征正交分解(POD)是提取離散數(shù)據(jù)主導(dǎo)特征的有效方法[9-10]。本文使用POD對帶有噪聲的全場測量數(shù)據(jù)Pfull進(jìn)行處理,構(gòu)建主導(dǎo)空間模態(tài) ψ 與對應(yīng)時間系數(shù)afull組成的線下數(shù)據(jù)庫,該過程對應(yīng)于圖1中紅色A所標(biāo)記的過程。采用POD可將流場Pfull分解為以下形式:

    圖1 重構(gòu)方法流程Fig.1 Outline of reconstruction method

    其中:k表示流場快照序號;M表示流場快照總數(shù);表示流場快照均值; ψi表示POD主導(dǎo)空間模態(tài);表示主導(dǎo)模態(tài)對應(yīng)的POD時間系數(shù)。該分解過程可利用奇異值分解(SVD)實現(xiàn):

    其中:λ表示各POD模態(tài)能量貢獻(xiàn)。

    以能量貢獻(xiàn)排序的前幾階POD模態(tài)捕獲了流場的主導(dǎo)流動結(jié)構(gòu),可用于構(gòu)建原始流場的低階近似模型:

    其中:m表示低階模型維度。構(gòu)建低階模型可簡化對復(fù)雜流場的描述,同時,通過選擇主導(dǎo)空間模態(tài)并忽略其他不重要的模態(tài),可在一定程度上起到濾除數(shù)據(jù)噪聲的效果。

    1.2 壓縮感知

    壓縮感知是一種尋找欠定線性系統(tǒng)的稀疏解的技術(shù),已在各個領(lǐng)域得到廣泛應(yīng)用[11-13]。壓縮感知理論指出,對于稀疏信號或可壓縮信號,可通過遠(yuǎn)低于奈奎斯特采樣率的采樣數(shù)據(jù)實現(xiàn)其精確重構(gòu)。

    本文通過求解L1范數(shù)下的最優(yōu)化問題求解流場重構(gòu)系數(shù)acs:

    其中:Pscattered表示線上稀疏采樣數(shù)據(jù);ψ表示線下所得POD主導(dǎo)模態(tài);矩陣 ?存儲離散傳感器的位置信息。詳情參見文獻(xiàn)[14]。最終,可通過流場重構(gòu)系數(shù)與POD主導(dǎo)模態(tài)的線性組合實現(xiàn)全場數(shù)據(jù)的重構(gòu):

    其中:m表示用于實現(xiàn)流場重構(gòu)的主導(dǎo)模態(tài)數(shù)量。

    1.3 卡爾曼濾波器

    卡爾曼濾波器提供了一種基于線性系統(tǒng)模型與測量數(shù)據(jù)對系統(tǒng)狀態(tài)進(jìn)行最優(yōu)估計的方法[15-16]。本文引入卡爾曼濾波器對壓縮感知求解所得流場重構(gòu)系數(shù)acs進(jìn)行實時去噪,進(jìn)一步提高流場重構(gòu)方法的魯棒性。卡爾曼濾波器中系統(tǒng)的線性模型可表示為:

    其中:x表示系統(tǒng)狀態(tài);A表示系統(tǒng)的線性狀態(tài)轉(zhuǎn)換矩陣;w表示協(xié)方差為Q的高斯分布的系統(tǒng)誤差??柭鼮V波器中的測量過程可表示為:

    其中:y表示測量結(jié)果;H表示觀測矩陣;v表示協(xié)方差為R的高斯分布的測量誤差。

    卡爾曼濾波器包含預(yù)估、矯正兩個步驟。預(yù)估步為:

    矯正步為:

    其中:上標(biāo) p表 示預(yù)估值;上標(biāo) c 表示校正值;F表示預(yù)估協(xié)方差矩陣;K表示卡爾曼增益。本文中,所要估測的系統(tǒng)狀態(tài)為POD時間系數(shù)a。

    1.4 動態(tài)模態(tài)分解

    采用動態(tài)模態(tài)分解構(gòu)建卡爾曼濾波器的線性系統(tǒng)模型,對非線性流動問題進(jìn)行線性化近似[17]:

    本文中,

    矩陣A的最佳形式可以表示為:

    其中:上標(biāo)+表示Moore-Penrose偽逆。

    2 數(shù)據(jù)準(zhǔn)備

    采用商業(yè)軟件ANSYS模擬處于圓柱尾渦影響下的翼型表面壓力變化,CFD模型示意圖如圖2所示,選用k-ω SST作為湍流模型,來流速度設(shè)置為30 m/s,翼型表面隨機設(shè)置12個離散壓力傳感器,傳感器布置位置如圖中紅色·所示。

    圖2 CFD模擬模型示意圖及離散傳感器布置位置Fig.2 Schematic diagram of CFD simulation model and deployment of discrete sensors

    為模擬線下全場采樣中可能產(chǎn)生的噪聲問題,對所模擬翼型表面壓力數(shù)據(jù)加入標(biāo)準(zhǔn)差為流場參數(shù)均值10%的高斯噪聲,所得單周期無噪聲與含噪聲流場的演化過程對比如圖3所示。同時,為模擬線上稀疏采樣過程中可能產(chǎn)生的噪聲問題,對所合成離散數(shù)據(jù)加入標(biāo)準(zhǔn)差為流場參數(shù)均值10%的高斯噪聲,所得單周期無噪聲和有噪聲的離散數(shù)據(jù)隨時間變化的曲線圖如圖4所示。

    3 方法驗證

    采用POD對線下所獲得含有噪聲的全場數(shù)據(jù)Pfull進(jìn)行處理,計算其主導(dǎo)空間模態(tài) ψ 與其對應(yīng)的時間系數(shù)afull。前四階POD空間模態(tài)和時間系數(shù)如圖5所示。POD主導(dǎo)空間模態(tài)捕獲了流場的主要流動特征,與之對應(yīng)的時間系數(shù)描述了其在時間上的動態(tài)演變。

    圖5 前4階POD空間模態(tài)及其對應(yīng)的時間系數(shù)Fig.5 The first four POD spatial modes and corresponding time coefficients

    選取前十階POD模態(tài)作為壓縮感知的稀疏基,以含有噪聲的線上采樣稀疏數(shù)據(jù)Pscattered作為觀測值,在L1范數(shù)最小化的條件下求解最優(yōu)時間系數(shù)acs。壓縮感知計算所得前四階模態(tài)對應(yīng)的時間系數(shù)如圖6所示,同樣計算所得時間系數(shù)嚴(yán)重扭曲,含有大量噪聲。直接通過這些含有大量噪聲的時間系數(shù)對流場進(jìn)行重構(gòu),重構(gòu)所得結(jié)果如圖7所示,與真實翼型表面壓力分布情況存在較大不符,無法滿足對翼型表面壓力進(jìn)行全場感知的要求。

    圖6 壓縮感知所得前4階模態(tài)時間系數(shù)Fig.6 The first four modal time coefficients obtained by compressed sensing

    圖7 壓縮感知所重構(gòu)翼型表面壓力Fig.7 Airfoil surface pressure reconstructed by compressed sensing

    為了解決上述噪聲干擾問題,在壓縮感知流場重構(gòu)方法的基礎(chǔ)上引入了卡爾曼濾波器。采用DMD對線下所得流場時間系數(shù)進(jìn)行處理,以DMD所得線性近似模型x′=Ax作為卡爾曼濾波器的系統(tǒng)模型;以線上壓縮感知計算得到的含有噪聲的時間系數(shù)作為觀測值進(jìn)行濾波??柭鼮V波器中過程噪聲方差矩陣Q和 測量噪聲方差矩陣R應(yīng)基于系統(tǒng)模型誤差與傳感器測量誤差的統(tǒng)計學(xué)信息設(shè)置,本文此處分別設(shè)置為0.0075Eeye(m)和 2Eeye(m),其中Eeye(m)表示m階單位矩陣。經(jīng)卡爾曼濾波器實時去噪后所得的時間系數(shù)得到顯著平滑,如圖8所示。利用濾波后的時間系數(shù)進(jìn)行重構(gòu),重構(gòu)所得結(jié)果如圖9所示,與各時刻翼型表面真實壓力分布吻合良好,滿足對翼型表面壓力進(jìn)行全場感知的要求。

    圖8 經(jīng)卡爾曼濾波器處理后前4階模態(tài)時間系數(shù)Fig.8 The first four modal time coefficients processed by Kalman filter

    圖9 經(jīng)卡爾曼濾波器所重構(gòu)翼型表面壓力Fig.9 Airfoil surface pressure reconstructed by Kalman filter

    基于翼型表面壓力數(shù)據(jù)積分可計算其升力l。經(jīng)壓縮感知所重構(gòu)升力和卡爾曼濾波器處理所重構(gòu)升力與真實升力變化曲線的對比如圖10所示??梢钥闯觯?jīng)壓縮感知所重構(gòu)升力在總體變化趨勢上與真實情況基本一致,但存在較嚴(yán)重的噪聲問題,無法滿足飛行器飛行過程中精準(zhǔn)感知的要求。在采用卡爾曼濾波器進(jìn)行處理之后,重構(gòu)結(jié)果與真實值之間的差異顯著降低,重構(gòu)結(jié)果基本滿足對升力進(jìn)行準(zhǔn)確感知的要求。

    圖10 重構(gòu)升力系數(shù)對比Fig.10 Comparison of lift coefficients reconstructed by compressed sensing and Kalman filter

    為量化卡爾曼濾波器對升力重構(gòu)準(zhǔn)確性的提升效果,定義升力的重構(gòu)誤差為:

    其中:M表示流場快照總數(shù);lrecon表示重構(gòu)所得升力;ltruth表示真實升力。經(jīng)壓縮感知所重構(gòu)流場與經(jīng)卡爾曼濾波器處理后所重構(gòu)升力誤差分別為6.39%和18.15%,誤差相較下降64.79%

    4 結(jié) 論

    本文提出了一種魯棒的基于數(shù)據(jù)融合的實時流場重構(gòu)方法,可以從有限的測量數(shù)據(jù)中實現(xiàn)全場數(shù)據(jù)的重構(gòu)。通過壓縮感知和DMD線性模型內(nèi)嵌的卡爾曼濾波器,顯著增強了其魯棒性,有效抑制了噪聲的影響。

    以翼型模擬數(shù)據(jù)為例測試了該方法的性能,并對是否采用卡爾曼濾波器進(jìn)行重構(gòu)的結(jié)果進(jìn)行了量化分析。對于翼型表面壓力數(shù)據(jù),升力重構(gòu)誤差從18.15%減小到6.39%,相較下降64.79%。

    另外值得注意的是,卡爾曼濾波器的成功應(yīng)用在很大程度上依賴于模型誤差方差矩陣及傳感器的誤差方差矩陣的準(zhǔn)確性,而在實際應(yīng)用中方差矩陣很難確定。為了對此進(jìn)行改進(jìn),我們目前正在開發(fā)基于DMD線性模型內(nèi)嵌的自適應(yīng)卡爾曼濾波器,引入對噪聲統(tǒng)計特性的實時動態(tài)估測,解決實際應(yīng)用中誤差方差矩陣難以準(zhǔn)確估測的問題。

    此外,本文目前所采用POD、DMD、卡爾曼濾波等流場分解與建模方法均為線性化方法,且DMD線性建模方法僅對較為簡單的周期性流動問題具有較好的適用性。后期考慮采用卷積神經(jīng)網(wǎng)絡(luò)、擴(kuò)展DMD、集合卡爾曼濾波等非線性建模方法,進(jìn)一步提高該方法的普適性。

    猜你喜歡
    升力卡爾曼濾波流場
    高速列車車頂–升力翼組合體氣動特性
    大型空冷汽輪發(fā)電機轉(zhuǎn)子三維流場計算
    無人機升力測試裝置設(shè)計及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計分析
    基于模糊卡爾曼濾波算法的動力電池SOC估計
    升力式再入飛行器體襟翼姿態(tài)控制方法
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    亚洲av综合色区一区| 日韩 欧美 亚洲 中文字幕| 国产精品香港三级国产av潘金莲 | 日韩,欧美,国产一区二区三区| 成年av动漫网址| 免费黄色在线免费观看| 亚洲国产欧美日韩在线播放| 亚洲人成网站在线观看播放| 最近最新中文字幕免费大全7| 考比视频在线观看| 黄片无遮挡物在线观看| 国产精品国产三级国产专区5o| 人人妻人人澡人人看| 青春草视频在线免费观看| 婷婷色综合www| 一本色道久久久久久精品综合| 人人妻人人澡人人爽人人夜夜| 丁香六月天网| 侵犯人妻中文字幕一二三四区| 精品一区二区三区av网在线观看 | 99九九在线精品视频| 亚洲一区二区三区欧美精品| 国产午夜精品一二区理论片| 妹子高潮喷水视频| 嫩草影院入口| 一本色道久久久久久精品综合| 免费女性裸体啪啪无遮挡网站| 精品一区二区三卡| 亚洲精品第二区| 咕卡用的链子| 丰满乱子伦码专区| bbb黄色大片| 黄色视频不卡| 久久国产精品男人的天堂亚洲| av女优亚洲男人天堂| 亚洲国产精品一区三区| 99久久精品国产亚洲精品| 国产有黄有色有爽视频| 国产又色又爽无遮挡免| 三上悠亚av全集在线观看| 老司机影院成人| 久久久精品94久久精品| 视频区图区小说| 久久99一区二区三区| 成人午夜精彩视频在线观看| 国产无遮挡羞羞视频在线观看| 成人漫画全彩无遮挡| 99热国产这里只有精品6| 国产高清不卡午夜福利| 国产日韩欧美在线精品| 亚洲精品日本国产第一区| 亚洲精品国产色婷婷电影| 街头女战士在线观看网站| 黄色一级大片看看| 国产色婷婷99| 国产精品一区二区精品视频观看| 亚洲欧美精品自产自拍| 咕卡用的链子| 两个人免费观看高清视频| 国产亚洲精品第一综合不卡| 在线观看一区二区三区激情| 19禁男女啪啪无遮挡网站| 国产99久久九九免费精品| 美女脱内裤让男人舔精品视频| 午夜激情av网站| 国产 一区精品| 一级爰片在线观看| 国产1区2区3区精品| 免费观看a级毛片全部| 欧美日本中文国产一区发布| 亚洲精品国产av蜜桃| 天堂俺去俺来也www色官网| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品视频女| 捣出白浆h1v1| 日韩 欧美 亚洲 中文字幕| 99久久人妻综合| 欧美亚洲日本最大视频资源| 欧美乱码精品一区二区三区| 成年女人毛片免费观看观看9 | 人人妻人人添人人爽欧美一区卜| 黄频高清免费视频| 精品人妻熟女毛片av久久网站| 日韩电影二区| 自线自在国产av| 国产乱来视频区| 国产精品嫩草影院av在线观看| 嫩草影院入口| xxxhd国产人妻xxx| 熟女少妇亚洲综合色aaa.| 日日撸夜夜添| 69精品国产乱码久久久| 美女大奶头黄色视频| 精品国产乱码久久久久久男人| 女人高潮潮喷娇喘18禁视频| 亚洲五月色婷婷综合| 精品少妇内射三级| 在线观看三级黄色| 国产黄色免费在线视频| 欧美日韩亚洲综合一区二区三区_| 国产不卡av网站在线观看| 免费黄网站久久成人精品| 中文乱码字字幕精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 久久精品久久久久久噜噜老黄| 久久99精品国语久久久| 电影成人av| 美女大奶头黄色视频| 亚洲国产精品一区二区三区在线| 在线天堂最新版资源| 久久久国产精品麻豆| 男女床上黄色一级片免费看| 国产伦理片在线播放av一区| 一区二区日韩欧美中文字幕| 秋霞伦理黄片| 久久狼人影院| 男人添女人高潮全过程视频| 久久精品熟女亚洲av麻豆精品| 自线自在国产av| 精品亚洲乱码少妇综合久久| 中文字幕亚洲精品专区| 久久久久久免费高清国产稀缺| 黄网站色视频无遮挡免费观看| 在线 av 中文字幕| 9191精品国产免费久久| 王馨瑶露胸无遮挡在线观看| 精品一区二区三区av网在线观看 | 黄片小视频在线播放| 日韩 欧美 亚洲 中文字幕| 国产亚洲av片在线观看秒播厂| 精品人妻一区二区三区麻豆| 国产亚洲av片在线观看秒播厂| 国产日韩欧美在线精品| 亚洲国产毛片av蜜桃av| 欧美av亚洲av综合av国产av | 最新的欧美精品一区二区| 亚洲精品日韩在线中文字幕| 亚洲人成电影观看| 两性夫妻黄色片| 国产精品偷伦视频观看了| 久久天躁狠狠躁夜夜2o2o | av电影中文网址| 校园人妻丝袜中文字幕| 成人漫画全彩无遮挡| 夫妻午夜视频| 2018国产大陆天天弄谢| 少妇人妻久久综合中文| 免费观看a级毛片全部| av福利片在线| 日韩一卡2卡3卡4卡2021年| 亚洲国产av影院在线观看| 日本vs欧美在线观看视频| 亚洲成国产人片在线观看| 午夜久久久在线观看| 国产精品 欧美亚洲| 又粗又硬又长又爽又黄的视频| 国产精品免费大片| 成人18禁高潮啪啪吃奶动态图| 丁香六月欧美| 国产高清国产精品国产三级| 成人18禁高潮啪啪吃奶动态图| 一级毛片电影观看| 国产激情久久老熟女| av又黄又爽大尺度在线免费看| bbb黄色大片| av免费观看日本| 国产精品一国产av| 观看av在线不卡| 欧美精品一区二区免费开放| 国产精品一国产av| 国产极品天堂在线| 国产精品无大码| 美女大奶头黄色视频| 午夜日本视频在线| 汤姆久久久久久久影院中文字幕| 国产精品久久久久久人妻精品电影 | 国产精品一国产av| 永久免费av网站大全| 成年美女黄网站色视频大全免费| 女性生殖器流出的白浆| 99re6热这里在线精品视频| 久久国产精品男人的天堂亚洲| 成年av动漫网址| 女性被躁到高潮视频| 超碰成人久久| 久热爱精品视频在线9| 秋霞在线观看毛片| 国产精品一区二区精品视频观看| 欧美日韩国产mv在线观看视频| 欧美日韩视频高清一区二区三区二| 天堂8中文在线网| 成人18禁高潮啪啪吃奶动态图| 成人三级做爰电影| 日韩av免费高清视频| 制服人妻中文乱码| 天堂8中文在线网| 免费观看性生交大片5| 午夜日韩欧美国产| 日韩免费高清中文字幕av| 热99国产精品久久久久久7| 亚洲一区中文字幕在线| 国产99久久九九免费精品| 中国国产av一级| 精品少妇久久久久久888优播| 亚洲国产欧美网| 国产精品久久久久成人av| 国产精品久久久久久精品电影小说| 美女国产高潮福利片在线看| 久久久久精品性色| 9热在线视频观看99| 一边摸一边做爽爽视频免费| 欧美日韩亚洲高清精品| 乱人伦中国视频| 日韩 欧美 亚洲 中文字幕| 男女无遮挡免费网站观看| 亚洲人成网站在线观看播放| 国产淫语在线视频| h视频一区二区三区| 在线精品无人区一区二区三| 亚洲久久久国产精品| 另类亚洲欧美激情| 黄片小视频在线播放| 男人添女人高潮全过程视频| 久热这里只有精品99| 黄片小视频在线播放| 美女高潮到喷水免费观看| 国产又色又爽无遮挡免| 午夜免费男女啪啪视频观看| 久久鲁丝午夜福利片| 九草在线视频观看| 成人亚洲欧美一区二区av| 夫妻性生交免费视频一级片| 国产av码专区亚洲av| 亚洲精品久久久久久婷婷小说| 亚洲激情五月婷婷啪啪| 最黄视频免费看| 久久99热这里只频精品6学生| 免费日韩欧美在线观看| 五月天丁香电影| 精品第一国产精品| 国产成人免费观看mmmm| 老鸭窝网址在线观看| 十八禁人妻一区二区| 精品一区二区免费观看| 啦啦啦在线观看免费高清www| av一本久久久久| 王馨瑶露胸无遮挡在线观看| 九色亚洲精品在线播放| 纯流量卡能插随身wifi吗| 久久人人97超碰香蕉20202| 十八禁高潮呻吟视频| 精品一品国产午夜福利视频| 亚洲男人天堂网一区| 国产欧美亚洲国产| 制服丝袜香蕉在线| 日韩大码丰满熟妇| 精品一区二区免费观看| 熟女少妇亚洲综合色aaa.| 日本av手机在线免费观看| av国产久精品久网站免费入址| 少妇人妻精品综合一区二区| 亚洲精品国产av蜜桃| 成人午夜精彩视频在线观看| 久久精品人人爽人人爽视色| 51午夜福利影视在线观看| 操出白浆在线播放| 国产精品欧美亚洲77777| 亚洲国产欧美网| 最近的中文字幕免费完整| 久久久久久久国产电影| 日本欧美视频一区| 老汉色∧v一级毛片| 亚洲成人免费av在线播放| 久久毛片免费看一区二区三区| 亚洲成av片中文字幕在线观看| 欧美日韩av久久| 啦啦啦在线免费观看视频4| 毛片一级片免费看久久久久| 午夜福利视频精品| 亚洲,欧美,日韩| 免费女性裸体啪啪无遮挡网站| 成人午夜精彩视频在线观看| 大香蕉久久成人网| √禁漫天堂资源中文www| 2021少妇久久久久久久久久久| 精品人妻在线不人妻| 亚洲一区二区三区欧美精品| 综合色丁香网| 国产精品亚洲av一区麻豆 | 永久免费av网站大全| 日韩免费高清中文字幕av| 视频区图区小说| 一级毛片我不卡| 国产一区有黄有色的免费视频| 男人操女人黄网站| 99热国产这里只有精品6| 欧美日韩精品网址| 国产av精品麻豆| 精品一品国产午夜福利视频| 精品国产乱码久久久久久小说| 国产伦人伦偷精品视频| av在线老鸭窝| 丝袜在线中文字幕| 免费av中文字幕在线| 操美女的视频在线观看| 成年人午夜在线观看视频| a级片在线免费高清观看视频| 桃花免费在线播放| 最近2019中文字幕mv第一页| 国产伦人伦偷精品视频| 黑人猛操日本美女一级片| 久久久国产欧美日韩av| 国产精品国产三级专区第一集| 日日爽夜夜爽网站| 亚洲图色成人| 亚洲av欧美aⅴ国产| 久久久久人妻精品一区果冻| 日韩视频在线欧美| 最黄视频免费看| 777久久人妻少妇嫩草av网站| 青春草视频在线免费观看| 亚洲av成人精品一二三区| 欧美国产精品va在线观看不卡| 美女高潮到喷水免费观看| 成人国产av品久久久| 嫩草影视91久久| 国产av精品麻豆| 中文字幕色久视频| 乱人伦中国视频| 精品免费久久久久久久清纯 | 国产黄频视频在线观看| 少妇 在线观看| 另类亚洲欧美激情| 日本猛色少妇xxxxx猛交久久| 侵犯人妻中文字幕一二三四区| 欧美人与善性xxx| 久久精品久久久久久噜噜老黄| 日韩中文字幕欧美一区二区 | 国产成人精品久久二区二区91 | 一个人免费看片子| 欧美日韩成人在线一区二区| 久久久久国产精品人妻一区二区| 成人亚洲精品一区在线观看| 午夜免费鲁丝| 日韩 欧美 亚洲 中文字幕| 亚洲国产av新网站| 97精品久久久久久久久久精品| 制服丝袜香蕉在线| 午夜日韩欧美国产| 亚洲精品一二三| 大片电影免费在线观看免费| 成年av动漫网址| 日本色播在线视频| 亚洲国产精品成人久久小说| 狠狠婷婷综合久久久久久88av| 免费av中文字幕在线| 男女无遮挡免费网站观看| av免费观看日本| 日韩一区二区三区影片| 精品免费久久久久久久清纯 | 色播在线永久视频| 午夜福利免费观看在线| 国产激情久久老熟女| 国产成人系列免费观看| 黄片无遮挡物在线观看| 一边摸一边抽搐一进一出视频| 人体艺术视频欧美日本| 久久人人爽av亚洲精品天堂| 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www| 日韩制服丝袜自拍偷拍| 天天躁日日躁夜夜躁夜夜| 日韩精品免费视频一区二区三区| 亚洲av成人精品一二三区| 久久久久久久国产电影| 国产女主播在线喷水免费视频网站| 观看美女的网站| 男女午夜视频在线观看| 黄片无遮挡物在线观看| 免费在线观看视频国产中文字幕亚洲 | 国产成人精品在线电影| 啦啦啦在线免费观看视频4| 男女无遮挡免费网站观看| xxx大片免费视频| 亚洲五月色婷婷综合| 国产精品 国内视频| 王馨瑶露胸无遮挡在线观看| 波野结衣二区三区在线| 亚洲国产精品一区三区| 国产成人一区二区在线| 如日韩欧美国产精品一区二区三区| 少妇 在线观看| 最近2019中文字幕mv第一页| 亚洲熟女毛片儿| 国产亚洲一区二区精品| 国产亚洲精品第一综合不卡| 国产又色又爽无遮挡免| 一级爰片在线观看| 在线观看人妻少妇| av线在线观看网站| www.av在线官网国产| 又粗又硬又长又爽又黄的视频| 亚洲成人免费av在线播放| 91国产中文字幕| 日韩中文字幕欧美一区二区 | 久久久久国产一级毛片高清牌| 一区二区三区激情视频| 只有这里有精品99| 国产无遮挡羞羞视频在线观看| 又大又爽又粗| 天堂俺去俺来也www色官网| 亚洲精品日韩在线中文字幕| 亚洲四区av| 亚洲欧美一区二区三区黑人| 欧美中文综合在线视频| 国产淫语在线视频| 69精品国产乱码久久久| 亚洲成人手机| 在线精品无人区一区二区三| 99热网站在线观看| 午夜福利视频在线观看免费| 亚洲欧美精品自产自拍| 热re99久久国产66热| 精品一品国产午夜福利视频| 亚洲欧洲日产国产| 日韩欧美一区视频在线观看| av又黄又爽大尺度在线免费看| 久久狼人影院| 爱豆传媒免费全集在线观看| 青春草视频在线免费观看| 美女主播在线视频| 中文字幕亚洲精品专区| 纵有疾风起免费观看全集完整版| av卡一久久| 色94色欧美一区二区| 综合色丁香网| 精品酒店卫生间| av天堂久久9| 宅男免费午夜| 桃花免费在线播放| 黄色 视频免费看| 哪个播放器可以免费观看大片| 国产日韩一区二区三区精品不卡| 美女午夜性视频免费| 高清av免费在线| 涩涩av久久男人的天堂| 国产精品亚洲av一区麻豆 | 中文字幕人妻熟女乱码| 欧美日韩精品网址| 久久天躁狠狠躁夜夜2o2o | 亚洲伊人色综图| 免费人妻精品一区二区三区视频| 999精品在线视频| 成年美女黄网站色视频大全免费| av女优亚洲男人天堂| 天天躁夜夜躁狠狠久久av| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 卡戴珊不雅视频在线播放| 午夜福利免费观看在线| 最近中文字幕高清免费大全6| 国产一区二区激情短视频 | 久久人妻熟女aⅴ| 亚洲自偷自拍图片 自拍| 久久99一区二区三区| 精品卡一卡二卡四卡免费| 精品国产露脸久久av麻豆| 成人18禁高潮啪啪吃奶动态图| 男女床上黄色一级片免费看| 日日爽夜夜爽网站| 99国产精品免费福利视频| 亚洲 欧美一区二区三区| 免费在线观看完整版高清| xxxhd国产人妻xxx| 亚洲成色77777| 啦啦啦视频在线资源免费观看| 老司机影院毛片| 在线观看免费高清a一片| 天天添夜夜摸| 2021少妇久久久久久久久久久| 亚洲成色77777| 高清欧美精品videossex| 人人妻人人澡人人爽人人夜夜| 一级毛片我不卡| 18禁国产床啪视频网站| 免费久久久久久久精品成人欧美视频| 1024香蕉在线观看| 视频区图区小说| 亚洲成人免费av在线播放| 男女下面插进去视频免费观看| 精品国产国语对白av| 99国产精品免费福利视频| 777米奇影视久久| 日本91视频免费播放| 青草久久国产| 国产精品一二三区在线看| 日韩人妻精品一区2区三区| 男男h啪啪无遮挡| 亚洲精品美女久久av网站| 性少妇av在线| 伊人久久国产一区二区| 国产麻豆69| 免费黄色在线免费观看| 亚洲一区二区三区欧美精品| 日韩 亚洲 欧美在线| 国产精品.久久久| 亚洲 欧美一区二区三区| 日日啪夜夜爽| 只有这里有精品99| 天堂俺去俺来也www色官网| 精品国产一区二区三区久久久樱花| 啦啦啦视频在线资源免费观看| 99国产精品免费福利视频| 成年女人毛片免费观看观看9 | 日本91视频免费播放| 免费在线观看黄色视频的| 欧美 日韩 精品 国产| 免费观看a级毛片全部| 亚洲自偷自拍图片 自拍| 国产av精品麻豆| 亚洲欧美一区二区三区久久| 高清视频免费观看一区二区| 中国国产av一级| 赤兔流量卡办理| 亚洲国产av新网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲一码二码三码区别大吗| 9热在线视频观看99| 国产熟女午夜一区二区三区| 国产无遮挡羞羞视频在线观看| 精品一区二区免费观看| 国产淫语在线视频| 男女国产视频网站| 欧美激情极品国产一区二区三区| 综合色丁香网| 免费观看人在逋| 最近的中文字幕免费完整| 成人毛片60女人毛片免费| 老司机深夜福利视频在线观看 | 伊人久久大香线蕉亚洲五| 制服诱惑二区| 纯流量卡能插随身wifi吗| 哪个播放器可以免费观看大片| 午夜影院在线不卡| 久久青草综合色| 国产又爽黄色视频| 国产99久久九九免费精品| 欧美日韩亚洲综合一区二区三区_| 欧美乱码精品一区二区三区| 亚洲精品日韩在线中文字幕| 97在线人人人人妻| e午夜精品久久久久久久| 18禁国产床啪视频网站| 欧美日韩视频高清一区二区三区二| 亚洲人成77777在线视频| 一级a爱视频在线免费观看| 亚洲国产欧美一区二区综合| 天天躁夜夜躁狠狠躁躁| 日本猛色少妇xxxxx猛交久久| 成人午夜精彩视频在线观看| 亚洲婷婷狠狠爱综合网| 超色免费av| 亚洲三区欧美一区| 国产精品无大码| 免费日韩欧美在线观看| 老司机影院毛片| 久久久久久久精品精品| 亚洲久久久国产精品| av片东京热男人的天堂| 国产又爽黄色视频| 久久这里只有精品19| av在线老鸭窝| 99久久精品国产亚洲精品| 免费看av在线观看网站| bbb黄色大片| 9热在线视频观看99| 久久久亚洲精品成人影院| 天天躁日日躁夜夜躁夜夜| 五月开心婷婷网| www.熟女人妻精品国产| 国产在线视频一区二区| 久久97久久精品| 久久久久久久久久久久大奶| 狂野欧美激情性bbbbbb| 最新的欧美精品一区二区| 亚洲少妇的诱惑av| e午夜精品久久久久久久| 黄片小视频在线播放| 我要看黄色一级片免费的| 久久久亚洲精品成人影院| 亚洲少妇的诱惑av| 91老司机精品| 伦理电影大哥的女人| 97人妻天天添夜夜摸| 免费不卡黄色视频| 巨乳人妻的诱惑在线观看| 免费观看性生交大片5| 母亲3免费完整高清在线观看| 丰满饥渴人妻一区二区三| 亚洲精品乱久久久久久| 高清视频免费观看一区二区| 黄色毛片三级朝国网站| 欧美变态另类bdsm刘玥| 日韩一卡2卡3卡4卡2021年| av在线老鸭窝| 激情五月婷婷亚洲| 欧美av亚洲av综合av国产av | 少妇猛男粗大的猛烈进出视频| 色网站视频免费| 亚洲国产精品成人久久小说| 国产精品成人在线| 欧美97在线视频| 嫩草影院入口| 欧美精品一区二区免费开放| 18禁动态无遮挡网站|