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

    錐體彈頭的嵌入式大氣數(shù)據(jù)傳感系統(tǒng)組合算法*

    2014-07-07 09:14:52肖前貴陸宇平肖地波
    傳感器與微系統(tǒng) 2014年10期
    關(guān)鍵詞:動(dòng)壓靜壓大氣

    張 勇,肖前貴,陸宇平,肖地波

    (1.南京航空航天大學(xué) 無人機(jī)研究院,江蘇 南京 210016; 2.南京航空航天大學(xué) 自動(dòng)化學(xué)院,江蘇 南京 210016; 3.南京航空航天大學(xué) 航天學(xué)院,江蘇 南京 210016)

    錐體彈頭的嵌入式大氣數(shù)據(jù)傳感系統(tǒng)組合算法*

    張 勇1,肖前貴1,陸宇平2,肖地波3

    (1.南京航空航天大學(xué) 無人機(jī)研究院,江蘇 南京 210016; 2.南京航空航天大學(xué) 自動(dòng)化學(xué)院,江蘇 南京 210016; 3.南京航空航天大學(xué) 航天學(xué)院,江蘇 南京 210016)

    嵌入式大氣數(shù)據(jù)傳感(FADS)系統(tǒng)空氣動(dòng)力學(xué)模型基于鈍頭體推導(dǎo),該模型適用于錐頭體解算大氣數(shù)據(jù),但存在一定的模型誤差。為消除模型誤差,分析了錐頭體模型誤差產(chǎn)生的機(jī)理,提出了一種FADS系統(tǒng)組合算法,并對(duì)某一典型錐體彈頭進(jìn)行了仿真驗(yàn)證。仿真結(jié)果表明:動(dòng)靜壓相對(duì)誤差在0.05 %以內(nèi),每組大氣數(shù)據(jù)解算平均耗時(shí)在20 ms以內(nèi),該算法有效提高了大氣數(shù)據(jù)解算精度和實(shí)時(shí)性。

    嵌入式大氣數(shù)據(jù)傳感系統(tǒng);錐頭體;模型誤差;組合算法

    0 引 言

    嵌入式大氣數(shù)據(jù)傳感(flush air data sensing,FADS)系統(tǒng)依靠分布在飛行器前端的壓力傳感器陣列來測(cè)量飛行器表面的壓力,并由壓力分布間接獲得飛行參數(shù)。FADS系統(tǒng)在大迎角飛行,減小雷達(dá)反射面積,高超聲速飛行等方面有著與傳統(tǒng)大氣數(shù)據(jù)測(cè)量系統(tǒng)無法比擬的優(yōu)勢(shì)[1]。在FADS系統(tǒng)的研究和應(yīng)用上,美國進(jìn)行了大量的研究,并取得大量成果,形成了成熟技術(shù)并應(yīng)用于航空工業(yè)中。國內(nèi)的FADS系統(tǒng)研究起步較晚,主要有方習(xí)高[2]、鄭成軍[3]、鄭守鐸[4]、葉瑋[5]、楊雨等人[6]在FADS系統(tǒng)框架、求解算法、故障檢測(cè)、FADS/INS系統(tǒng)組合以及互補(bǔ)濾波技術(shù)在FADS系統(tǒng)中的應(yīng)用方面取得了一系列階段性成果,上述研究成果均以鈍頭體飛行器為研究對(duì)象。

    對(duì)圓錐體外形,要建立表面壓力分布與大氣數(shù)據(jù)參數(shù)的模型非常困難。已有的FADS系統(tǒng)空氣動(dòng)力學(xué)模型是在球體上建立的,采用形壓系數(shù)來適應(yīng)在非球體外形上的應(yīng)用,這使得該模型有可能適用于圓錐體外形。目前,只有文獻(xiàn)[7]驗(yàn)證了FADS系統(tǒng)空氣動(dòng)力學(xué)模型在錐頭體飛行器適用性問題,仿真結(jié)果表明,根據(jù)鈍頭體推導(dǎo)的FADS系統(tǒng)空氣動(dòng)力學(xué)模型應(yīng)用于錐頭體飛行器存在一定的模型誤差,但文獻(xiàn)[7]未進(jìn)一步對(duì)模型誤差機(jī)理進(jìn)行分析。為進(jìn)一步消除模型誤差,文獻(xiàn)[7]中采用多項(xiàng)式擬合校準(zhǔn)的方法,該方法雖然有效減小了模型誤差,但模型誤差改善效果會(huì)隨著多項(xiàng)式次數(shù)的增加而趨于飽和。因此,根據(jù)模型誤差產(chǎn)生的機(jī)理和文獻(xiàn)[7]中解算方法的局限性,本文提出了一種精度和實(shí)時(shí)性更好的FADS系統(tǒng)解算方法。

    1 FADS系統(tǒng)模型

    對(duì)于典型的鈍頭體,F(xiàn)ADS系統(tǒng)動(dòng)力學(xué)模型為[7]

    p(θi)=qc(cos2θi+εsin2θi)+p∞,

    (1)

    cosθi=cosαecosβecosλi+sinβesinφisinλi+ sinαecosβecosφisinλi,

    (2)

    (3)

    ε=f(αe,βe,Ma).

    (4)

    其中,p(θi),θi,φi,λi分別為第i個(gè)測(cè)壓孔的壓力值、入射角、圓周角、圓錐角;qc為動(dòng)壓,p∞為靜壓,ε為無量綱的形壓系數(shù),αe,βe分別為當(dāng)?shù)赜呛彤?dāng)?shù)貍?cè)滑角,Ma為馬赫數(shù)。

    2 FADS系統(tǒng)模型誤差分析

    上述模型中式(1)~式(3)是顯式的數(shù)學(xué)公式,而式(4)是需要根據(jù)實(shí)驗(yàn)數(shù)據(jù)飛行前確定。該模型是在球體上建立的,采用形壓系數(shù)來適應(yīng)在非球體外形上的應(yīng)用。

    根據(jù)壓力分布情況計(jì)算當(dāng)?shù)赜铅羍和當(dāng)?shù)貍?cè)滑角βe[8],利用式(2)計(jì)算各測(cè)壓孔處的氣流入射角θi。

    對(duì)于一組飛行參數(shù)(α,β,Ma),利用cos2θi=1-sin2θi,將式(1)改寫成

    (5)

    寫成矩陣形式

    工作區(qū)地下水主要接受大氣降水入滲補(bǔ)給和地表水側(cè)向補(bǔ)給。巖溶裂隙、孔隙是地下水的主要徑流通道,地下水總體徑流方向與地表水徑流方向基本一致,即由北西向南東徑流。地下水排泄形式以泉、潛流及井為主。

    P=AX.

    (6)

    其中

    P=[p(θ1) …p(θn)]T,

    (7)

    (8)

    (9)

    根據(jù)線性方程組的求解理論,得到X的最小二乘解

    X=A+P.

    (10)

    再根據(jù)式(3)計(jì)算出不同馬赫數(shù)時(shí)的qc/p∞,然后就可以算出該組飛行參數(shù)對(duì)應(yīng)的形壓系數(shù)

    (11)

    即得到了一組(α,β,Ma)對(duì)應(yīng)的ε。取一系列的(α,β,Ma)k組合,重復(fù)以上步驟,可以得到對(duì)應(yīng)的形壓系數(shù)序列εk,即得到ε隨(α,β,Ma)變化的關(guān)系。

    根據(jù)已知的馬赫數(shù),利用式(3),求出動(dòng)靜壓之比

    (12)

    由于式(3)中馬赫數(shù)與動(dòng)靜壓之比的關(guān)系對(duì)各種外形(包括錐頭體和鈍頭體)的飛行器均適用的,因此,可認(rèn)為式(12)得到的動(dòng)靜壓之比是足夠精確的,可以由此來確定真實(shí)的動(dòng)靜壓。

    簡(jiǎn)化分析,對(duì)動(dòng)靜壓進(jìn)行歸一化處理,即靜壓為1,動(dòng)壓就等于動(dòng)靜壓之比,即得到動(dòng)靜壓的真實(shí)值

    (13)

    采用式(5)~式(10)的過程求出X,再根據(jù)已知馬赫數(shù)求出動(dòng)靜壓之比g,可以得到根據(jù)模型(1)計(jì)算的動(dòng)壓和靜壓值

    (14)

    模型動(dòng)壓和靜壓的相對(duì)誤差

    (15)

    由式(14)和式(15)可見,模型動(dòng)壓誤差和模型靜壓誤差相當(dāng)。實(shí)際上,根據(jù)式(9) 和式(10)可以看出,調(diào)整形壓系數(shù)ε可以改變靜壓或者動(dòng)壓,但是動(dòng)靜壓之和X(1)與形壓系數(shù)無關(guān)。式(1)用形壓系數(shù)來滿足不同外形的前提是式(1)必須能夠測(cè)得準(zhǔn)確的總壓(動(dòng)靜壓之和),然后通過調(diào)整形壓系數(shù)可以改變動(dòng)靜壓在總壓中占的比重。而通過以上分析發(fā)現(xiàn),對(duì)文獻(xiàn)[7]中研究的錐頭體外形,式(1)并不能準(zhǔn)確算出總壓。因此,通過調(diào)整形壓系數(shù),可以使得動(dòng)壓或者靜壓中的一個(gè)準(zhǔn)確,但卻使得另外一個(gè)更加不準(zhǔn)確,即無法通過改變形壓系數(shù)使得動(dòng)壓和靜壓同時(shí)精確。

    3 FADS系統(tǒng)組合算法

    由上一節(jié)分析可知,通過調(diào)整形壓系數(shù)不能減小FADS系統(tǒng)模型誤差。式(1)~式(4)描述的模型誤差主要是動(dòng)壓和靜壓誤差,迎角和側(cè)滑角誤差基本由校準(zhǔn)方法決定。

    雖然單個(gè)神經(jīng)網(wǎng)絡(luò)的方法能滿足大氣數(shù)據(jù)的計(jì)算精度要求,但由于采用單一神經(jīng)網(wǎng)絡(luò)計(jì)算所有大氣數(shù)據(jù)增大了故障管理難度和訓(xùn)練復(fù)雜度[10]。因此,為有效消除模型誤差和提高工程實(shí)現(xiàn)性,采用分層計(jì)算的方式實(shí)現(xiàn)FADS系統(tǒng)解算。采用三點(diǎn)法和校準(zhǔn)算法求解迎角和側(cè)滑角,再采用神經(jīng)網(wǎng)絡(luò)校準(zhǔn)動(dòng)靜壓,降低網(wǎng)絡(luò)的輸出維數(shù),最后根據(jù)計(jì)算所得的動(dòng)靜壓值求解馬赫數(shù)。

    這樣可以將整個(gè)FADS系統(tǒng)算法分成三層:三點(diǎn)法和校準(zhǔn)算法求解迎角α和側(cè)滑角β;神經(jīng)網(wǎng)絡(luò)法求解動(dòng)壓qc和靜壓p∞;根據(jù)動(dòng)靜壓計(jì)算馬赫數(shù)Ma。組合算法流程圖如圖1。

    圖1 組合算法流程Fig 1 Flow chart of combinatorial algorithm

    RBF網(wǎng)絡(luò)具有強(qiáng)大的函數(shù)逼近功能,并且比BP網(wǎng)絡(luò)設(shè)計(jì)更加簡(jiǎn)單,它將自動(dòng)選擇神經(jīng)元個(gè)數(shù),傳輸函數(shù)為徑向基函數(shù)。組合算法中動(dòng)靜壓求解采用RBF神經(jīng)網(wǎng)絡(luò),RBF神經(jīng)網(wǎng)絡(luò)構(gòu)建方法參見文獻(xiàn)[11]。神經(jīng)網(wǎng)絡(luò)在訓(xùn)練中確定了權(quán)值和閾值,在線計(jì)算也是簡(jiǎn)單的數(shù)學(xué)運(yùn)算,因此,該算法進(jìn)一步提高了算法的實(shí)時(shí)性。

    4 仿真驗(yàn)證

    采用文獻(xiàn)[7]中的錐體彈頭外形和壓力數(shù)據(jù),采用組合算法進(jìn)行大氣數(shù)據(jù)解算。由于迎角和側(cè)滑角的解算與文獻(xiàn)[7]相同,并且錐頭體的模型誤差體現(xiàn)在動(dòng)靜壓誤差方面,因此,本文著重比較動(dòng)靜壓的求解精度。

    本次設(shè)計(jì)中,RBF神經(jīng)網(wǎng)絡(luò)求解動(dòng)靜壓,其徑向基函數(shù)的密度常數(shù)選取為spread=0.55。將迎角α=-8°~8°、側(cè)滑角β=0°,±2°,±4°時(shí)的壓力數(shù)據(jù)作為RBF神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本數(shù)據(jù),并將α=-8°~8°,β=±1°,±3°時(shí)壓力數(shù)據(jù)作為測(cè)試樣本。組合算法解算得到的動(dòng)壓與靜壓相對(duì)誤差如圖 2和圖 3。

    圖2 動(dòng)壓相對(duì)誤差Fig 2 Relative errors of dynamic pressure

    圖3 靜壓相對(duì)誤差Fig 3 Relative errors of static pressure

    由圖 2和圖 3可知,本文的組合算法相比文獻(xiàn)[7]中的校準(zhǔn)算法精度更高,動(dòng)壓相對(duì)誤差控制在0.05 %以內(nèi),靜壓相對(duì)誤差達(dá)到10-12量級(jí)。

    進(jìn)一步,仿真對(duì)比本文組合算法和文獻(xiàn)[7]中的校準(zhǔn)算法解算時(shí)間。在同一臺(tái)計(jì)算機(jī)對(duì)兩種算法進(jìn)行測(cè)試,測(cè)試128組大氣數(shù)據(jù),得到每組大氣數(shù)據(jù)平均耗時(shí)如表 1。由表1可知,本文組合算法的平均一組大氣數(shù)據(jù)解算時(shí)間小于20 ms,相比文獻(xiàn)[7]大大縮短了大氣數(shù)據(jù)的解算時(shí)間。

    表1 平均一組大氣數(shù)據(jù)解算時(shí)間Tab 1 Average resolving time of a set of air datas

    5 結(jié) 論

    本文針對(duì)FADS系統(tǒng)空氣動(dòng)力學(xué)模型應(yīng)用于錐頭體飛行器存在模型誤差的機(jī)理進(jìn)行了分析,分析表明:通過調(diào)整形壓系數(shù)不能減小FADS系統(tǒng)的模型誤差。針對(duì)模型誤差產(chǎn)生的機(jī)理,提出了一種FADS系統(tǒng)組合算法,并對(duì)其進(jìn)行了數(shù)值仿真驗(yàn)證。仿真結(jié)果表明:動(dòng)壓與靜壓的相對(duì)誤差小于0.05 %,每組大氣數(shù)據(jù)解算時(shí)間小于20 ms,該算法有效提高了算法解算精度和實(shí)時(shí)性。

    [1] Srivastava Ankur ,Meade Andrew J,Mokhtarzadeh Ali Arya.A hybrid data-model fusion approach to calibrate a flush air data sensing system[R].Atlanta:AIAA,2010.

    [2] 方習(xí)高,陸玉平.嵌入式大氣數(shù)據(jù)傳感系統(tǒng)求解算法研究[M].計(jì)算機(jī)測(cè)量與控制,2008,16(3):398-400.

    [3] Zheng Chenjun,Lu Yuping,He zhen.Improved algorithms for flush airdata sensing system[J].Chinese Journal of Aeronatics,2006,9(4):334-339.

    [4] 鄭守鐸,陸宇平,葉 瑋.基于χ2檢驗(yàn)的FADS系統(tǒng)故障檢測(cè)與管理技術(shù)研究[J].計(jì)算機(jī)測(cè)量與控制,2007,15(11):1449-1551.

    [5] 葉 瑋,鄭守鐸,溫瑞珩.FADS/INS組合法迎角、側(cè)滑角測(cè)量方法研究[J].飛機(jī)設(shè)計(jì),2007,27(6):14-17.

    [6] 楊 雨,陸宇平,吳在桂.嵌入式大氣數(shù)據(jù)傳感系統(tǒng)中的組合濾波技術(shù)[J].傳感器與微系統(tǒng),2009,28(5):117-120.

    [7] 張 勇,肖前貴,陸宇平,等.一種適用錐頭體的嵌入式大氣數(shù)據(jù)傳感系統(tǒng)改進(jìn)校準(zhǔn)算法[J].計(jì)算機(jī)測(cè)量與控制,2014,22(5):1513-1515.

    [8] Cobleigh Brent R,Whitmore Stephen A,Haering Edward A.Flush airdata sensing (FADS) system calibration procedures and results for blunt forebodies[C]∥Dryden Flight Research Center Edwards,California,1999.

    [9] Edward J Artz,Nicholas W Dona,Dr Thomas R.Yechout.NASA orion flush air data sensing system feasibility determination and development[C]∥The 52nd Aerospace Sciences Meeting,Maryland,2014.

    [10] 肖地波,陸宇平,吳在桂.嵌入式大氣數(shù)據(jù)傳感系統(tǒng)求解算法研究[M].北京:科學(xué)出版社,2010.

    [11] 魏海坤.神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)設(shè)計(jì)的理論與方法[M].北京:國防工業(yè)出版社,2005.

    Combinatorial algorithm of flush air data sensing system for cone forebodies*

    ZHANG Yong1, XIAO Qian-gui1, LU Yu-ping2, XIAO Di-bo3

    (1.Research Institute of Unmanned Aircraft Vehicle,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China; 2.College of Automatic Engineering,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China; 3.College of Astronautics,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China)

    The aerodynamic model of flush air data sensing (FADS) system,which is derived from blunt forebodies,is suitable for cone forebodies,however,the model error is subsistent when the aerodynamic model is applied to cone forebodies.To eliminate model error,the mechanism of cone forebody model error is analyzed and a suit of combinatorial algorithm for FADS is proposed,the corresponding validation of simulation is accomplished for a typical cone forebody.Simulation results show that relative errors of dynamic pressure and static pressure are less than 0.05 % and average resolving time consuming of a set of air data is less than 20 ms,respectively,thus solving precision and real-time performance are improved effectively by the algorithm.

    flush air data sensing(FADS) system; cone forebody; model error; combinatorial algorithm

    10.13873/J.1000—9787(2014)10—0148—03

    2014—07—28

    江蘇省自然科學(xué)基金資助項(xiàng)目(BK20130817);南京航空航天大學(xué)基本科研業(yè)務(wù)費(fèi)資助項(xiàng)目(NS2014404)

    V 241.6

    A

    1000—9787(2014)10—0148—03

    張 勇(1986-),男,浙江衢州人,工學(xué)碩士,研究實(shí)習(xí)員,主要從事飛行控制方向的研究。

    猜你喜歡
    動(dòng)壓靜壓大氣
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    靜壓法沉樁對(duì)周邊環(huán)境影響及質(zhì)量控制
    國內(nèi)首個(gè)現(xiàn)代箔片氣體動(dòng)壓軸承技術(shù)培訓(xùn)班在長(zhǎng)沙成功舉辦
    靜壓托換樁在某濕陷性黃土場(chǎng)地地基加固中的應(yīng)用
    超精密液體靜壓轉(zhuǎn)臺(tái)裝配技術(shù)
    一種基于空氣靜壓支承的自調(diào)心裝置
    大氣古樸揮灑自如
    大氣、水之后,土十條來了
    南屯煤礦深部泵房硐室群動(dòng)壓失穩(wěn)機(jī)理及控制對(duì)策
    強(qiáng)烈動(dòng)壓巷道支護(hù)技術(shù)探討
    亚洲精华国产精华精| 另类亚洲欧美激情| 在线国产一区二区在线| x7x7x7水蜜桃| 日韩人妻精品一区2区三区| 欧美不卡视频在线免费观看 | 国产无遮挡羞羞视频在线观看| 级片在线观看| 中亚洲国语对白在线视频| 色哟哟哟哟哟哟| 女警被强在线播放| 久久久精品欧美日韩精品| 女人被躁到高潮嗷嗷叫费观| 我的亚洲天堂| 一区在线观看完整版| 亚洲精品粉嫩美女一区| 亚洲 国产 在线| 91老司机精品| 亚洲人成电影免费在线| 天堂动漫精品| 亚洲专区中文字幕在线| 91在线观看av| 久久久久国产一级毛片高清牌| 天堂中文最新版在线下载| 91成人精品电影| 欧美日韩瑟瑟在线播放| 熟女少妇亚洲综合色aaa.| 亚洲五月天丁香| 99国产精品免费福利视频| 欧美日韩亚洲高清精品| 咕卡用的链子| 大香蕉久久成人网| av在线天堂中文字幕 | 欧美一级毛片孕妇| 女生性感内裤真人,穿戴方法视频| 日日夜夜操网爽| 俄罗斯特黄特色一大片| 精品一品国产午夜福利视频| 日韩欧美在线二视频| 中文字幕另类日韩欧美亚洲嫩草| 国产蜜桃级精品一区二区三区| 岛国在线观看网站| 首页视频小说图片口味搜索| 十八禁网站免费在线| 丝袜在线中文字幕| 亚洲情色 制服丝袜| 黄色视频,在线免费观看| 999久久久国产精品视频| 91精品国产国语对白视频| 99精国产麻豆久久婷婷| 亚洲欧美精品综合一区二区三区| 亚洲一区高清亚洲精品| 另类亚洲欧美激情| 99久久99久久久精品蜜桃| 亚洲免费av在线视频| 亚洲av日韩精品久久久久久密| 男女午夜视频在线观看| 麻豆成人av在线观看| 黄色怎么调成土黄色| 最新在线观看一区二区三区| 久久精品91无色码中文字幕| 久久久水蜜桃国产精品网| 欧美日本亚洲视频在线播放| 午夜激情av网站| 亚洲情色 制服丝袜| 欧美人与性动交α欧美精品济南到| 亚洲狠狠婷婷综合久久图片| 男人舔女人下体高潮全视频| 亚洲熟女毛片儿| 真人做人爱边吃奶动态| 亚洲国产中文字幕在线视频| 精品一区二区三区四区五区乱码| 亚洲七黄色美女视频| 国产又色又爽无遮挡免费看| 亚洲五月色婷婷综合| 国产精品一区二区免费欧美| 亚洲免费av在线视频| 免费看十八禁软件| 国产精品国产av在线观看| 国产99白浆流出| 一区二区三区国产精品乱码| 日韩一卡2卡3卡4卡2021年| 五月开心婷婷网| 搡老岳熟女国产| 欧美日韩av久久| 黄色a级毛片大全视频| 欧美日韩一级在线毛片| 国产精品av久久久久免费| 日本 av在线| 欧美亚洲日本最大视频资源| 欧美日韩精品网址| 亚洲精品av麻豆狂野| 精品国产一区二区久久| 男女之事视频高清在线观看| 日韩欧美一区二区三区在线观看| 国产乱人伦免费视频| 久久伊人香网站| 亚洲精华国产精华精| 欧美日韩亚洲综合一区二区三区_| 久久天堂一区二区三区四区| 婷婷精品国产亚洲av在线| 成年版毛片免费区| 成人精品一区二区免费| 男女高潮啪啪啪动态图| 日本一区二区免费在线视频| 日韩人妻精品一区2区三区| 91成年电影在线观看| 亚洲国产毛片av蜜桃av| 妹子高潮喷水视频| 欧美 亚洲 国产 日韩一| av网站在线播放免费| 亚洲精品中文字幕一二三四区| 精品国产美女av久久久久小说| 国产高清视频在线播放一区| 黑人操中国人逼视频| 伦理电影免费视频| 在线av久久热| 男人操女人黄网站| 日韩人妻精品一区2区三区| 丝袜美腿诱惑在线| 日本精品一区二区三区蜜桃| 男人舔女人下体高潮全视频| 99国产极品粉嫩在线观看| 777久久人妻少妇嫩草av网站| 丰满迷人的少妇在线观看| 久久久国产成人免费| 亚洲熟妇熟女久久| 两个人免费观看高清视频| 波多野结衣高清无吗| 精品少妇一区二区三区视频日本电影| 99精品在免费线老司机午夜| 久久久久国产一级毛片高清牌| 国产av又大| 免费搜索国产男女视频| 午夜福利免费观看在线| 岛国视频午夜一区免费看| 亚洲精品久久成人aⅴ小说| 国产av在哪里看| 12—13女人毛片做爰片一| 精品少妇一区二区三区视频日本电影| 精品久久久精品久久久| 看片在线看免费视频| 国产精品永久免费网站| 三级毛片av免费| 国产精品偷伦视频观看了| 少妇裸体淫交视频免费看高清 | 欧美大码av| 日韩一卡2卡3卡4卡2021年| 国产色视频综合| 国产成人精品久久二区二区91| av有码第一页| av欧美777| 婷婷精品国产亚洲av在线| 亚洲精品国产区一区二| 九色亚洲精品在线播放| 欧美激情久久久久久爽电影 | 成熟少妇高潮喷水视频| 欧美中文日本在线观看视频| aaaaa片日本免费| 国产国语露脸激情在线看| 成人黄色视频免费在线看| 免费人成视频x8x8入口观看| 午夜亚洲福利在线播放| 天天添夜夜摸| 亚洲熟妇熟女久久| 一进一出好大好爽视频| 亚洲成人免费av在线播放| 成人av一区二区三区在线看| 变态另类成人亚洲欧美熟女 | 亚洲欧美日韩高清在线视频| 搡老乐熟女国产| 免费日韩欧美在线观看| 极品人妻少妇av视频| 黄片大片在线免费观看| av网站免费在线观看视频| 亚洲av第一区精品v没综合| 久久久国产成人精品二区 | 亚洲欧美一区二区三区久久| 人人澡人人妻人| 99久久综合精品五月天人人| 国产欧美日韩综合在线一区二区| 亚洲一区二区三区欧美精品| 在线看a的网站| 一夜夜www| 日本a在线网址| 精品国产乱码久久久久久男人| 国产精品自产拍在线观看55亚洲| 色综合站精品国产| 啦啦啦 在线观看视频| 日本精品一区二区三区蜜桃| 精品人妻在线不人妻| 国产xxxxx性猛交| 久久人人97超碰香蕉20202| 麻豆久久精品国产亚洲av | 亚洲精品粉嫩美女一区| 宅男免费午夜| 1024香蕉在线观看| 日韩 欧美 亚洲 中文字幕| 在线观看一区二区三区激情| 亚洲少妇的诱惑av| 婷婷精品国产亚洲av在线| 99国产极品粉嫩在线观看| 波多野结衣高清无吗| 亚洲少妇的诱惑av| 999久久久国产精品视频| 伊人久久大香线蕉亚洲五| 日本免费一区二区三区高清不卡 | 男男h啪啪无遮挡| 国产在线精品亚洲第一网站| 男女高潮啪啪啪动态图| 亚洲人成电影观看| 18禁裸乳无遮挡免费网站照片 | 午夜a级毛片| 日韩免费av在线播放| 日韩精品青青久久久久久| 99国产极品粉嫩在线观看| 国产男靠女视频免费网站| 婷婷精品国产亚洲av在线| 欧美日本亚洲视频在线播放| 亚洲伊人色综图| 女人高潮潮喷娇喘18禁视频| 99久久人妻综合| 精品电影一区二区在线| 亚洲精品在线观看二区| 亚洲美女黄片视频| 日韩 欧美 亚洲 中文字幕| 80岁老熟妇乱子伦牲交| 亚洲精品一二三| www.999成人在线观看| 19禁男女啪啪无遮挡网站| 精品久久久久久久毛片微露脸| 视频区欧美日本亚洲| 免费高清在线观看日韩| 久久中文字幕一级| 99国产综合亚洲精品| 天天添夜夜摸| 一夜夜www| 一级片免费观看大全| 嫁个100分男人电影在线观看| 国产欧美日韩一区二区精品| 精品久久久久久电影网| 校园春色视频在线观看| 亚洲在线自拍视频| 国产亚洲精品综合一区在线观看 | 亚洲人成伊人成综合网2020| 97超级碰碰碰精品色视频在线观看| 久久久久国产一级毛片高清牌| 国产片内射在线| 黄色毛片三级朝国网站| 咕卡用的链子| 精品无人区乱码1区二区| 真人一进一出gif抽搐免费| 午夜老司机福利片| 久久久久国产一级毛片高清牌| 国产高清videossex| 50天的宝宝边吃奶边哭怎么回事| 窝窝影院91人妻| 麻豆成人av在线观看| 久久人人精品亚洲av| 国产野战对白在线观看| 国产精品久久久av美女十八| 宅男免费午夜| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 久久伊人香网站| 日韩大尺度精品在线看网址 | 在线av久久热| 欧美中文综合在线视频| 国产激情欧美一区二区| 日韩国内少妇激情av| 黑人巨大精品欧美一区二区mp4| 悠悠久久av| 无遮挡黄片免费观看| 一区二区日韩欧美中文字幕| 宅男免费午夜| 日韩三级视频一区二区三区| 女性生殖器流出的白浆| 大陆偷拍与自拍| www.www免费av| 99久久国产精品久久久| 精品久久久久久成人av| 日日爽夜夜爽网站| 欧美激情久久久久久爽电影 | 丝袜美足系列| 在线观看免费视频日本深夜| 免费看a级黄色片| 亚洲熟妇中文字幕五十中出 | 国产高清videossex| 在线国产一区二区在线| 又紧又爽又黄一区二区| 老熟妇仑乱视频hdxx| www.自偷自拍.com| 一区二区日韩欧美中文字幕| 国产精品一区二区三区四区久久 | 18禁美女被吸乳视频| 制服诱惑二区| 人妻久久中文字幕网| 男男h啪啪无遮挡| 亚洲熟女毛片儿| 国产在线观看jvid| 成年人免费黄色播放视频| cao死你这个sao货| 日韩有码中文字幕| svipshipincom国产片| 亚洲精品国产一区二区精华液| 亚洲欧美一区二区三区久久| av电影中文网址| 黑人巨大精品欧美一区二区mp4| 久久人人97超碰香蕉20202| 一区二区日韩欧美中文字幕| 淫秽高清视频在线观看| 在线观看午夜福利视频| 国产97色在线日韩免费| 精品国内亚洲2022精品成人| 国产精品1区2区在线观看.| 女人高潮潮喷娇喘18禁视频| 日本精品一区二区三区蜜桃| 免费观看人在逋| 欧美成狂野欧美在线观看| 中出人妻视频一区二区| 看免费av毛片| 成年版毛片免费区| 亚洲 国产 在线| 亚洲视频免费观看视频| 99国产综合亚洲精品| 真人一进一出gif抽搐免费| 久久久国产一区二区| 真人做人爱边吃奶动态| 精品国产超薄肉色丝袜足j| 国产在线观看jvid| 亚洲av日韩精品久久久久久密| 日韩大码丰满熟妇| 老司机福利观看| 女性被躁到高潮视频| www.自偷自拍.com| 免费高清视频大片| 亚洲一区二区三区欧美精品| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品av麻豆狂野| 又紧又爽又黄一区二区| 精品久久久久久,| 国产一区二区激情短视频| 久久精品91无色码中文字幕| 99久久精品国产亚洲精品| 国产成人av激情在线播放| 9热在线视频观看99| 好男人电影高清在线观看| 亚洲av美国av| 老鸭窝网址在线观看| 精品久久蜜臀av无| 精品久久久精品久久久| 韩国精品一区二区三区| 免费高清视频大片| 大香蕉久久成人网| 久久伊人香网站| 欧美人与性动交α欧美软件| 男人舔女人的私密视频| 久久天堂一区二区三区四区| 免费女性裸体啪啪无遮挡网站| 久久香蕉激情| 久久久久久久久免费视频了| 久热爱精品视频在线9| 亚洲国产看品久久| 男女之事视频高清在线观看| 夜夜夜夜夜久久久久| 男男h啪啪无遮挡| 三上悠亚av全集在线观看| 国产精品免费一区二区三区在线| 交换朋友夫妻互换小说| www国产在线视频色| 亚洲国产欧美网| 99久久久亚洲精品蜜臀av| 一进一出好大好爽视频| 午夜福利,免费看| 亚洲一码二码三码区别大吗| 亚洲一区高清亚洲精品| 韩国av一区二区三区四区| 三上悠亚av全集在线观看| 午夜福利,免费看| 国产精品 国内视频| 女人被狂操c到高潮| 午夜91福利影院| av中文乱码字幕在线| 露出奶头的视频| 国产极品粉嫩免费观看在线| 黄片播放在线免费| 精品一区二区三区视频在线观看免费 | 欧美一区二区精品小视频在线| 国产精品久久久人人做人人爽| 男女下面插进去视频免费观看| 成人国语在线视频| 女人爽到高潮嗷嗷叫在线视频| 久久婷婷成人综合色麻豆| 亚洲av美国av| 久久精品亚洲av国产电影网| 久久午夜综合久久蜜桃| 水蜜桃什么品种好| 另类亚洲欧美激情| a在线观看视频网站| 两个人看的免费小视频| 欧美精品亚洲一区二区| 国产高清激情床上av| 亚洲一码二码三码区别大吗| 久久草成人影院| 高清欧美精品videossex| av中文乱码字幕在线| 极品教师在线免费播放| e午夜精品久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 欧美日本亚洲视频在线播放| 免费女性裸体啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 91在线观看av| 午夜影院日韩av| 性色av乱码一区二区三区2| aaaaa片日本免费| 91九色精品人成在线观看| 亚洲精品在线观看二区| 国产有黄有色有爽视频| 可以免费在线观看a视频的电影网站| 超碰成人久久| 一区福利在线观看| 亚洲aⅴ乱码一区二区在线播放 | 男人的好看免费观看在线视频 | 日本一区二区免费在线视频| 亚洲av成人一区二区三| 在线视频色国产色| 交换朋友夫妻互换小说| 国产高清videossex| 在线观看www视频免费| 亚洲精品美女久久久久99蜜臀| 亚洲五月婷婷丁香| 欧美久久黑人一区二区| 99久久久亚洲精品蜜臀av| 视频区欧美日本亚洲| 欧美黄色淫秽网站| 老司机亚洲免费影院| 国产免费av片在线观看野外av| 免费看十八禁软件| 日韩大尺度精品在线看网址 | 亚洲色图av天堂| 亚洲成人国产一区在线观看| 18禁黄网站禁片午夜丰满| 亚洲欧美一区二区三区久久| 亚洲第一av免费看| 日韩人妻精品一区2区三区| 国产高清国产精品国产三级| 老汉色∧v一级毛片| 女性生殖器流出的白浆| 久久久久久亚洲精品国产蜜桃av| 色婷婷久久久亚洲欧美| 午夜福利一区二区在线看| 这个男人来自地球电影免费观看| 久久影院123| 18禁观看日本| 国产亚洲精品久久久久5区| 国产激情久久老熟女| 国产精品久久久人人做人人爽| 亚洲七黄色美女视频| 国产高清videossex| 老汉色∧v一级毛片| 国产免费现黄频在线看| 久久久国产成人精品二区 | 18禁国产床啪视频网站| 高清毛片免费观看视频网站 | 精品久久久精品久久久| 久久狼人影院| 成人黄色视频免费在线看| 国产亚洲精品一区二区www| 国产又色又爽无遮挡免费看| 大香蕉久久成人网| 麻豆一二三区av精品| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲五月天丁香| 亚洲精品av麻豆狂野| 天天影视国产精品| 久久九九热精品免费| 好看av亚洲va欧美ⅴa在| 一进一出好大好爽视频| 国产麻豆69| 九色亚洲精品在线播放| 日本黄色日本黄色录像| 免费在线观看黄色视频的| 欧美性长视频在线观看| 亚洲第一青青草原| aaaaa片日本免费| 亚洲精品成人av观看孕妇| videosex国产| 香蕉丝袜av| 国产高清国产精品国产三级| www.自偷自拍.com| 麻豆av在线久日| 在线国产一区二区在线| 99re在线观看精品视频| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看 | 成人亚洲精品av一区二区 | 满18在线观看网站| 国产精品秋霞免费鲁丝片| 精品国内亚洲2022精品成人| 欧美日韩亚洲高清精品| 在线免费观看的www视频| 国产极品粉嫩免费观看在线| 亚洲午夜精品一区,二区,三区| 91精品三级在线观看| 欧美黄色淫秽网站| 侵犯人妻中文字幕一二三四区| 国产精品九九99| 精品午夜福利视频在线观看一区| 亚洲熟女毛片儿| 香蕉丝袜av| 我的亚洲天堂| 国产真人三级小视频在线观看| 黄片大片在线免费观看| 99国产精品99久久久久| 一区福利在线观看| 女警被强在线播放| 国产精品免费视频内射| 国产一区二区三区在线臀色熟女 | 亚洲伊人色综图| 一边摸一边抽搐一进一小说| 亚洲午夜理论影院| 他把我摸到了高潮在线观看| 极品人妻少妇av视频| 欧美成人午夜精品| 欧美激情高清一区二区三区| 国产精品一区二区免费欧美| 天天影视国产精品| 操出白浆在线播放| 又紧又爽又黄一区二区| 黄色视频,在线免费观看| 国产人伦9x9x在线观看| 色在线成人网| 99香蕉大伊视频| 亚洲av五月六月丁香网| 国产精品爽爽va在线观看网站 | 淫妇啪啪啪对白视频| 成人三级黄色视频| 欧美色视频一区免费| 美女国产高潮福利片在线看| 久久精品人人爽人人爽视色| 国产日韩一区二区三区精品不卡| 亚洲国产毛片av蜜桃av| 99热国产这里只有精品6| 在线观看66精品国产| 国产av精品麻豆| 午夜免费鲁丝| 国产主播在线观看一区二区| 日本vs欧美在线观看视频| 国产无遮挡羞羞视频在线观看| 亚洲精品中文字幕一二三四区| 丁香六月欧美| 国产一卡二卡三卡精品| 久久久久国产一级毛片高清牌| 欧美日韩亚洲高清精品| 这个男人来自地球电影免费观看| 国产伦一二天堂av在线观看| 国产黄a三级三级三级人| svipshipincom国产片| 欧美日本亚洲视频在线播放| 国产欧美日韩一区二区精品| 久热爱精品视频在线9| av网站免费在线观看视频| 亚洲专区国产一区二区| 如日韩欧美国产精品一区二区三区| 一进一出抽搐gif免费好疼 | 精品一区二区三区av网在线观看| 欧美在线黄色| 亚洲一区高清亚洲精品| 久久人人97超碰香蕉20202| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 99国产精品一区二区蜜桃av| 美国免费a级毛片| 91成年电影在线观看| 可以在线观看毛片的网站| 深夜精品福利| 黄色丝袜av网址大全| 国产主播在线观看一区二区| 精品午夜福利视频在线观看一区| 女人高潮潮喷娇喘18禁视频| 黄色片一级片一级黄色片| 一进一出抽搐gif免费好疼 | 亚洲一区中文字幕在线| 国产亚洲精品久久久久5区| 久久精品国产综合久久久| 午夜福利免费观看在线| 可以在线观看毛片的网站| 十八禁网站免费在线| 9色porny在线观看| 婷婷精品国产亚洲av在线| 亚洲一区二区三区欧美精品| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 在线观看舔阴道视频| 亚洲精品美女久久av网站| 国产亚洲欧美98| 纯流量卡能插随身wifi吗| 亚洲五月婷婷丁香| 在线观看舔阴道视频| 88av欧美| 9色porny在线观看| 美女高潮喷水抽搐中文字幕| 成人18禁高潮啪啪吃奶动态图| 日韩欧美一区二区三区在线观看| 亚洲男人的天堂狠狠| 在线观看免费视频网站a站| 校园春色视频在线观看| av在线播放免费不卡| 欧美成人免费av一区二区三区| 波多野结衣av一区二区av| 国产成人欧美在线观看| 国产精品久久久av美女十八| 可以免费在线观看a视频的电影网站| 欧美日本亚洲视频在线播放|