亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        遠(yuǎn)程迎頭探測(cè)飛機(jī)目標(biāo)的紅外輻射特性研究

        2017-03-22 10:52:19康麗珠趙勁松
        紅外技術(shù) 2017年4期
        關(guān)鍵詞:頂角輻射強(qiáng)度紅外光

        康麗珠,趙勁松,周 倩,倪 凱,唐 晗,趙 強(qiáng),陶 亮

        ?

        遠(yuǎn)程迎頭探測(cè)飛機(jī)目標(biāo)的紅外輻射特性研究

        康麗珠1,趙勁松1,周 倩2,倪 凱2,唐 晗1,趙 強(qiáng)1,陶 亮1

        (1. 昆明物理研究所,云南 昆明 650223;2. 清華大學(xué)深圳研究生院,深圳 518055)

        研究了遠(yuǎn)程迎頭探測(cè)飛機(jī)目標(biāo)的紅外輻射特性。分別介紹了飛機(jī)迎頭的主要紅外輻射源,飛機(jī)紅外輻射特性的計(jì)算流程,以及紅外光電系統(tǒng)入瞳接收到的蒙皮紅外輻射能通量的計(jì)算模型,最后對(duì)遠(yuǎn)程探測(cè)水平路徑大氣透過率和飛機(jī)目標(biāo)紅外輻射強(qiáng)度進(jìn)行仿真計(jì)算。

        紅外輻射特性;飛機(jī)目標(biāo);紅外輻射強(qiáng)度

        0 引言

        隨著隱身戰(zhàn)機(jī)在現(xiàn)代戰(zhàn)場(chǎng)上的廣泛應(yīng)用,特別是F-22隱身戰(zhàn)機(jī)的出現(xiàn),給現(xiàn)代雷達(dá)系統(tǒng)帶來了巨大的挑戰(zhàn),紅外光電系統(tǒng)能顯著彌補(bǔ)雷達(dá)對(duì)隱身空中目標(biāo)探測(cè)能力的不足,因而遠(yuǎn)程探測(cè)紅外光電系統(tǒng)的研制恰逢其時(shí)[1-4]。在此背景下,國(guó)內(nèi)外持續(xù)20多年對(duì)飛機(jī)的紅外輻射特性進(jìn)行研究,為紅外光電系統(tǒng)的設(shè)計(jì)與評(píng)價(jià)提供理論依據(jù),對(duì)于遠(yuǎn)程探測(cè)和預(yù)警敵方戰(zhàn)機(jī)具有重要的現(xiàn)實(shí)意義。

        但是,紅外隱身技術(shù)的發(fā)展使飛機(jī)的紅外輻射特性呈現(xiàn)顯著減小的趨勢(shì),極大地增大了紅外探測(cè)的困難。與此形成鮮明對(duì)比的是,飛機(jī)超音速飛行時(shí),尾焰、尾噴管和蒙皮均會(huì)產(chǎn)生較大的紅外輻射強(qiáng)度[5]。根據(jù)參考文獻(xiàn)[6]中圖9分析,F(xiàn)-22在左右機(jī)動(dòng)、俯沖、拉升、翻滾時(shí)的紅外輻射特征都較強(qiáng),尤其在突破音障時(shí),高溫尾氣流包裹著機(jī)身,形成強(qiáng)烈的紅外輻射[6]。F-22紅外輻射最小值發(fā)生在亞音速飛行時(shí)迎頭方向,這是紅外光電系統(tǒng)預(yù)警和探測(cè)的難點(diǎn),也是研究的重點(diǎn)。本文對(duì)飛機(jī)亞音速飛行時(shí)迎頭方向探測(cè)的紅外輻射特性進(jìn)行深入研究和仿真分析。

        1 飛機(jī)紅外輻射源及計(jì)算流程

        1.1 紅外輻射源

        本文參考北約空中目標(biāo)紅外輻射特性計(jì)算軟件NIRATAM[7-9],將亞音速飛行時(shí)飛機(jī)迎頭紅外輻射源分為兩類,如圖1所示。

        1)飛機(jī)蒙皮輻射:氣動(dòng)加熱、內(nèi)部熱源產(chǎn)生的紅外輻射。飛機(jī)在飛行過程中與空氣摩擦產(chǎn)生氣動(dòng)加熱的熱量,發(fā)動(dòng)機(jī)、尾噴管、噴嘴等熱部件向外傳導(dǎo)的熱量對(duì)蒙皮進(jìn)行的再加溫,這些因素對(duì)飛機(jī)蒙皮紅外輻射量的貢獻(xiàn)是不可忽視的。

        2)飛機(jī)蒙皮受到背景的輻射:飛機(jī)蒙皮受到太陽、天空、大地的背景輻射也會(huì)產(chǎn)生反射紅外輻射量。

        圖1 飛機(jī)的典型紅外輻射源

        上述各部分輻射源的比重依賴于飛機(jī)飛行參數(shù)、飛機(jī)蒙皮特性參數(shù)、飛機(jī)結(jié)構(gòu)性能參數(shù)和背景條件等因素。飛機(jī)飛行參數(shù)包括飛機(jī)姿態(tài)、飛行速度、飛行高度等。飛機(jī)蒙皮特性參數(shù)包括蒙皮吸收率、蒙皮發(fā)射率、蒙皮反射率,以及機(jī)身材料密度、比熱容、導(dǎo)熱系數(shù)等物理參數(shù)。飛機(jī)結(jié)構(gòu)性能參數(shù)包括飛機(jī)零部件的位置、尺寸與形狀等參數(shù)。背景條件包括天氣狀況、地面狀況以及飛機(jī)與太陽地球的相對(duì)位置等。

        1.2 計(jì)算流程

        飛機(jī)目標(biāo)紅外輻射特性的計(jì)算是本研究的重要內(nèi)容。在計(jì)算飛機(jī)紅外輻射量時(shí),需考慮機(jī)頭與內(nèi)部熱部件及外部環(huán)境進(jìn)行的能量交換,因此,飛機(jī)紅外輻射特性的求解是非常復(fù)雜的,本文將在一定合理簡(jiǎn)化的基礎(chǔ)上開展研究。另外,飛機(jī)紅外輻射通過大氣傳輸?shù)郊t外光電系統(tǒng)的過程中,必定會(huì)伴隨著輻射的衰減,因此計(jì)算中還需要考慮紅外大氣透過率等參數(shù)。圖2給出了飛機(jī)亞音速飛行時(shí)迎頭紅外輻射特性的計(jì)算流程。在計(jì)算中,我們依據(jù)紅外輻射源建立了相應(yīng)的計(jì)算模型。

        2 飛機(jī)輻射能通量計(jì)算模型

        2.1 位置關(guān)系的確定

        受空間位置關(guān)系的影響,入射到飛機(jī)不同零部件表面的紅外輻射能量存在一定的差異,不僅如此,甚至飛機(jī)同一部件的不同局部表面也不相同。隨著時(shí)間的變化,太陽光與飛機(jī)蒙皮各部分夾角也不停變化,具體情況根據(jù)飛機(jī)所處在的實(shí)際位置決定。因此,計(jì)算中需要弄清飛機(jī)、太陽和地球的位置關(guān)系。為此,本文建立3個(gè)坐標(biāo)系:球坐標(biāo)系(--)、地表坐標(biāo)系(--)、飛機(jī)坐標(biāo)系(--),如圖3所示。

        圖2 飛機(jī)迎頭紅外輻射特性計(jì)算流程

        同樣,當(dāng)飛機(jī)處于不同的空間位置和飛行姿態(tài)時(shí),入射到紅外光電系統(tǒng)的紅外輻射量也是不相同的,為此需要確定飛機(jī)和紅外光電系統(tǒng)兩者之間的空間幾何位置關(guān)系,如圖4所示。以機(jī)頭為觀測(cè)的起點(diǎn),觀測(cè)角以逆時(shí)針方向?yàn)檎瑒t紅外光電系統(tǒng)正迎頭探測(cè)的角度(觀測(cè)角)為0°。

        2.2 熱部件傳導(dǎo)部分

        紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮發(fā)射熱部件傳導(dǎo)的輻射能通量unit()如公式(1)所示:

        式中:A1為飛機(jī)蒙皮接受內(nèi)部熱部件傳熱的有效面積;假設(shè)熱部件發(fā)射面元dAu和飛機(jī)吸收面元dA1之間進(jìn)行熱交換;l為兩面元中心連線的距離;a1和a2為面元法線與面元連線夾角;eu為熱部件發(fā)射率;Eunit(l)為飛機(jī)熱部件工作時(shí)T1產(chǎn)生的紅外輻照度;e為蒙皮發(fā)射率;t(l)為飛機(jī)到紅外光電探測(cè)系統(tǒng)的大氣透過率;u為紅外光電探測(cè)系統(tǒng)對(duì)飛機(jī)的觀測(cè)角。

        圖4 紅外光電系統(tǒng)與飛機(jī)的位置關(guān)系

        2.3 由氣動(dòng)加熱產(chǎn)生部分

        當(dāng)氣體迎面流向鈍頭物體時(shí),在頭部附近受到阻滯,其前端有一個(gè)滯止點(diǎn),流動(dòng)速度為零,稱為駐點(diǎn),其溫度為滯止溫度0,用公式(2)計(jì)算[10-12]:

        式中:為比熱容比;amb為周圍環(huán)境大氣溫度;為飛機(jī)飛行時(shí)的馬赫數(shù)。

        恢復(fù)溫度不等于氣流滯止溫度,這是由于氣體內(nèi)部有熱傳導(dǎo)而引起的差異。飛機(jī)恢復(fù)溫度可用公式(3)計(jì)算[13-14]:

        式中:R為飛機(jī)的恢復(fù)溫度;為恢復(fù)系數(shù)。

        紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮發(fā)射由氣動(dòng)加熱產(chǎn)生的輻射能通量flow()如公式(4)所示:

        總而言之,農(nóng)村氣象服務(wù)建設(shè)作為一項(xiàng)長(zhǎng)期性、系統(tǒng)性、復(fù)雜性工程,必須投入更多人力、物力、財(cái)力,全面提升農(nóng)村氣象的整體服務(wù)水平。加大氣象服務(wù)工程建設(shè)力度,為有效預(yù)防和處理農(nóng)業(yè)災(zāi)害、增加農(nóng)民收入和農(nóng)作物產(chǎn)量發(fā)揮出不可取代的作用。此外,農(nóng)村氣象服務(wù)是新農(nóng)村重點(diǎn)建設(shè)內(nèi)容,對(duì)于社會(huì)進(jìn)步和國(guó)家繁榮發(fā)展等目標(biāo)實(shí)現(xiàn)提供了有力條件。

        式中:2為蒙皮發(fā)射由氣動(dòng)加熱產(chǎn)生的輻射的有效面積;flow()為由氣動(dòng)加熱導(dǎo)致蒙皮升溫產(chǎn)生的輻照度。

        2.4 太陽直接照射部分

        式中:3為蒙皮發(fā)射太陽直接照射的有效面積;s為飛機(jī)蒙皮對(duì)太陽輻射的吸收率;s()為太陽到飛機(jī)的大氣透過率;太陽天頂角度取決于目標(biāo)與太陽的相對(duì)位置關(guān)系;sun()為太陽輻照度。

        2.5 反射太陽部分

        太陽天頂角的計(jì)算式如公式(6)所示[15]:

        cos=sinsin+coscoscos (6)

        式中:是觀察者所在地理緯度;是太陽赤緯角;是太陽時(shí)角。

        將反射分為兩部分:一是蒙皮直接對(duì)太陽光的反射,二是蒙皮反射來自地球的太陽反照光。紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮反射太陽的輻射能通量sun-r()如公式(7)所示:

        2.6 地球自身輻射部分

        紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮反射地球自身輻射部分的輻射能通量earth()如公式(8)所示:

        式中:5為飛機(jī)蒙皮反射地球輻射的有效面積;e為地球發(fā)射率;e()為地球到飛機(jī)的大氣透過率;為飛機(jī)蒙皮法線與地球之間的夾角。

        2.7 天空自身輻射部分

        紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮反射天空的輻射能通量sky()如公式(9)所示:

        式中:6為飛機(jī)蒙皮反射天空輻射的有效面積;s為天空發(fā)射率;sky()為天空自身輻射部分對(duì)飛機(jī)的輻照度;為飛機(jī)蒙皮法線與天空之間的夾角。

        2.8 紅外光電系統(tǒng)入瞳接收到的輻射總量

        綜上所述,紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮發(fā)射或反射的總輻射能通量radi()如公式(10)所示:

        radi()=unit()+flow()+sun-d()+

        sun-r()+earth()+sky() (10)

        3 飛機(jī)紅外輻射量仿真

        本文根據(jù)上述“飛機(jī)輻射能通量計(jì)算模型”的推導(dǎo)過程,編制“飛機(jī)目標(biāo)紅外輻射量計(jì)算程序”APIR[16],并計(jì)算紅外光電系統(tǒng)入瞳接收到的F-22迎頭探測(cè)時(shí)的蒙皮輻射量。首先,計(jì)算水平路徑海拔10km處遠(yuǎn)距離條件下的大氣透過率值,為后續(xù)紅外輻射強(qiáng)度的計(jì)算提供參考;其次,分析不同天頂角對(duì)蒙皮反射太陽紅外輻射強(qiáng)度的影響;最后,分析不同距離下紅外光電系統(tǒng)入瞳接收到的紅外輻射強(qiáng)度總量。

        3.1 遠(yuǎn)距水平路徑大氣透過率的仿真

        針對(duì)飛機(jī)目標(biāo)遠(yuǎn)程迎頭探測(cè)的特點(diǎn),利用“紅外大氣透過率計(jì)算程序”(IRAT)計(jì)算水平路徑海拔10km處遠(yuǎn)距離條件下的大氣透過率值[17]。距離分別取250km、300km、350km、400km,假設(shè)地面溫度為27℃,地面相對(duì)濕度為65%,能見度為50km。水平路徑3~5mm大氣透過率(海拔:10km,距離:150~400km)如圖5所示,水平路徑8~12mm大氣透過率(海拔:10km,距離:150~400km)如圖6所示。

        圖5 水平路徑3~5mm大氣透過率(海拔:10km,距離:150~400km)

        圖6 水平路徑8~12mm大氣透過率(海拔:10km,距離:150~400km)

        3.2 天頂角影響的仿真

        根據(jù)公式(6)編制APIR中天頂角的仿真計(jì)算程序,假設(shè)研究地址為北京(經(jīng)度116.41667°,緯度39.91667°),日期為2012年7月1日,天頂角計(jì)算程序界面如圖7所示。根據(jù)公式(7)編制APIR中反射太陽紅外輻射量的計(jì)算程序,分析不同天頂角對(duì)蒙皮反射太陽的紅外輻射強(qiáng)度的影響。假設(shè)海拔高度為10km,地球發(fā)射率為0.3,馬赫數(shù)為0.6,觀測(cè)角度為0°,在不同的太陽天頂角影響下,紅外光電系統(tǒng)入瞳接收到的蒙皮反射太陽的輻射強(qiáng)度值如表1所示,天頂角和紅外輻射強(qiáng)度的關(guān)系如圖8所示,橫坐標(biāo)為不同時(shí)刻的天頂角,縱坐標(biāo)為紅外光電系統(tǒng)入瞳接收到的飛機(jī)蒙皮反射太陽的輻射強(qiáng)度。

        由圖8可知,正午12點(diǎn)太陽直射時(shí)天頂角最小,此時(shí)紅外光電系統(tǒng)入瞳接收到的蒙皮反射太陽的輻射強(qiáng)度最強(qiáng),反射太陽輻射強(qiáng)度隨著天頂角的增加而較小,當(dāng)天頂角q>90°時(shí),太陽輻射強(qiáng)度為0,這說明飛機(jī)輻射量隨時(shí)間(或太陽天頂角)是呈現(xiàn)晝夜周期變化的。蒙皮反射太陽輻射部分主要以中波為主。以上結(jié)論與文獻(xiàn)[18]給出的結(jié)論基本一致。

        圖7 天頂角計(jì)算程序界面

        表1 天頂角及蒙皮反射太陽輻射強(qiáng)度值

        圖8 不同時(shí)刻蒙皮反射太陽的紅外輻射強(qiáng)度

        3.3 紅外光電系統(tǒng)入瞳接收到的蒙皮輻射總量

        利用APIR做仿真計(jì)算,分析紅外光電系統(tǒng)入瞳接收到F-22迎頭探測(cè)時(shí)的總蒙皮輻射量。假設(shè)地址為北京(經(jīng)度116.41667°,緯度39.91667°),日期為2012年7月1日,時(shí)間為10:0:0,海拔高度為10km,觀測(cè)角度為0°,蒙皮發(fā)射率為0.2,地球發(fā)射率為0.3,馬赫數(shù)為0.6,紅外光電系統(tǒng)到飛機(jī)的距離為250km、300km、350km、400km,采用IRAT計(jì)算的遠(yuǎn)距水平路徑大氣透過率值。在不同距離下,紅外光電系統(tǒng)入瞳接收到的紅外輻射強(qiáng)度總量如表2所示,距離與紅外輻射強(qiáng)度總量的關(guān)系如圖9所示,橫坐標(biāo)為紅外光電系統(tǒng)到飛機(jī)的距離,縱坐標(biāo)為紅外光電系統(tǒng)入瞳接收到F-22迎頭探測(cè)時(shí)的蒙皮輻射強(qiáng)度總量。

        圖9 距離對(duì)蒙皮輻射強(qiáng)度總量的影響

        由圖9可知,隨著距離增加,紅外光電系統(tǒng)入瞳接收到蒙皮輻射強(qiáng)度總量呈現(xiàn)逐漸減小的趨勢(shì),主要原因是大氣透過率隨著距離的增加而減小,從而導(dǎo)致蒙皮輻射強(qiáng)度總量隨之減小。該結(jié)論與文獻(xiàn)[18]給出的結(jié)論基本一致。

        表2 不同距離下的總輻射強(qiáng)度值

        4 結(jié)束語

        本文主要研究了遠(yuǎn)程迎頭探測(cè)飛機(jī)目標(biāo)的紅外輻射特性。計(jì)算過程中綜合考慮了飛機(jī)蒙皮所吸收的由氣動(dòng)加熱、內(nèi)部熱源產(chǎn)生的紅外輻射,以及飛機(jī)蒙皮反射太陽、地球、天空的紅外輻射。依據(jù)輻射源的分類,本文建立了紅外光電系統(tǒng)入瞳接收到飛機(jī)輻射能通量的計(jì)算模型。并根據(jù)該計(jì)算模型編制了計(jì)算程序,重點(diǎn)仿真計(jì)算了不同天頂角的飛機(jī)蒙皮反射太陽輻射強(qiáng)度和不同距離下的蒙皮輻射強(qiáng)度總量,得到的結(jié)論與文獻(xiàn)[18]給出的結(jié)論基本一致,由此證明了本文研究結(jié)果的正確性。

        [1] 王慧君, 趙寅, 倪樹新. 機(jī)載紅外預(yù)警探測(cè)系統(tǒng)[J]. 激光與紅外, 2007, 37(11): 1133-1136.

        WANG Huijun, ZHAO Yin, NI Shuxin. Airborne IR early warn detect system[J]., 2007, 37(11): 1133-1136.

        [2] Ohammad Rahamat Kalam, Akhtar Akella Ravikiran. A cost-effective early warning system for improving the reliability of power system[J]., 2013, 9(3): 1-6.

        [3] Jadi Amr, Zedan Hussein, Alghamdi Turki. Risk management based early warning system for healthcare industry[C]//, 2013: 1-6(DOI: 1001109//ICCMA.2013.6506181).

        [4] Allison Fusswinkel. The advanced hawkeye as a system of systems[C]//, 2007, 4: 1-12.

        [5] TONG Zhongxiang, LU Yanlong, LI Chuanliang , et al. Optimal method for interfering with infrared homing missile by airborne infrared smoke[C]//, 2008: 588-592( DOI: 10.1109/ASCI-ICSC.2008.4576428).

        [6] 康麗珠, 趙勁松, 李振華, 等. 飛機(jī)目標(biāo)紅外輻射特性研究現(xiàn)狀概述[J]. 紅外技術(shù), 2017, 39(2): 105-115.

        KANG Lizhu, ZHAO Jinsong, LI Zhenhua, et al. The overview of the research work developments on infrared signature of aircrafts[J]., 2017, 39(2): 105-115

        [7] BEIER K. Infrared radiation model for aircraft and reentry vehicle[C]//, 1988, 972: 363- 374.

        [8] Bakker E J. Modeling multi-spectral imagery data with NIRATAM v3.1 and NPLUME v1.6[C]//, 1991, 3699: 275-282.

        [9] Sidonie Lefebvre, Antoine Roblin, Suzanne Varet, et al. A methodological approach for statistical evaluation of aircraft infrared signature[J]., 2010, 95: 484-493.

        [10] YU Yanwen, LIU Yiyun. Generation of realistic infrared image for moving objects[J]., 2004, 25(7): 1087-1097.

        [11] John D McGlynn, Steven P Auerbach. IR signature prediction errors for skin-heated aerial targets[C]//, 3062: 22-30.

        [12]Hudson R D.[M]. New York: Wiley-Interscience, 1968.

        [13] 苗瑞生, 居賢銘, 吳甲生. 導(dǎo)彈空氣動(dòng)力學(xué)[M]. 北京: 國(guó)防工業(yè)出版, 2006: 15-39.

        MIAO Ruisheng, JU Xianming, WU Jiasheng.[M]. Beijing: National Defense Industry Press, 2006: 15-39.

        [14] LIU Juan. Researches on modeling and simulation of airplane infrared characteristic[C]//2010, 2010, 5: 47-51(DOI: 10.1109//ICCASM.2010.5619361).

        [15] 劉鑒民. 太陽能利用原理技術(shù)工程[M]. 北京: 電子工業(yè)出版社, 2010.

        LIU Jianmin.[M]. Beijing: Electronic Industry Press. 2010.

        [16] KANG Lizhu. The computation program of Infrared radiation intensity of aircrafts (APIR)[Z]. Kunming institute of physics, 2016. (internal sources)

        [17] KANG Lizhu. The computation program of infrared atmospheric transmittance(IRAT)[Z]. Kunming institute of physics, 2016. (internal sources)

        [18] Campana S B.:[M]. Bellingham: SPIE-The International Society for Optical Engineering, 1993.

        Research on Infrared Signature for Remotely Detection from the Nose of Aircrafts

        KANG Lizhu1,ZHAO Jinsong1,ZHOU Qian2,NI Kai2,TANG Han1,ZHAO Qiang1,TAO Liang1

        (1.,650223,; 2.,,518055,)

        The infrared signature for remotely detection from the nose of aircrafts is studied in this paper. The major infrared sources, the calculation process of infrared signature, and the calculation model of infrared radiation flux which reach into the entrance pupil of infrared imaging systems is introduced respectively. Finally, the atmospheric transmittance from horizontal remote path and the infrared radiation intensity are simulated.

        infrared,infrared signature,aircraft target,Infrared radiation intensity

        TN214

        A

        1001-8891(2017)04-0365-07

        2017-02-18;

        2017-03-17.

        康麗珠,女(1983-),博士研究生,主要從事紅外光學(xué)技術(shù)的研究。E-mail:kangbingxin@126.com。

        猜你喜歡
        頂角輻射強(qiáng)度紅外光
        神奇窗戶讓室內(nèi)四季如春
        一般三棱鏡最大頂角與折射率的關(guān)系
        甲硝唑配合紅外光治療慢性宮頸炎的有效性及對(duì)復(fù)發(fā)率的影響
        石墨烯纖維可應(yīng)用于中紅外光電子器件通訊
        涼亭中的數(shù)學(xué)
        氧氣A(O,O)波段氣輝體發(fā)射率和臨邊輻射強(qiáng)度模擬與分析
        頂角為100°的等腰三角形性質(zhì)的應(yīng)用
        有關(guān)紅外光電子物理研究的幾個(gè)問題
        基于模擬太陽輻射強(qiáng)度對(duì)自然循環(huán)式PV/T系統(tǒng)的實(shí)驗(yàn)研究
        紅外地球敏感器輻射強(qiáng)度在軌數(shù)據(jù)分析*
        国产午夜福利在线观看红一片| 国产一极毛片| 天天躁人人躁人人躁狂躁| 国产午夜视频免费观看| 久久久久久国产福利网站| 国产美女胸大一区二区三区| 成人av在线免费播放| 久久九九精品国产不卡一区| 蜜桃av精品一区二区三区| 柠檬福利第一导航在线| 欧美操逼视频| 亚洲色AV性色在线观看| 激情在线视频一区二区三区| 蜜桃视频羞羞在线观看| 狠狠躁18三区二区一区| 久久免费的精品国产v∧| 日本专区一区二区三区| 国产又大大紧一区二区三区| 一区在线视频免费播放| 中国孕妇变态孕交xxxx| 国产精品久久久久av福利动漫| 久久久AV无码精品免费| 久久精品国产亚洲不卡| 中文字幕免费人成在线网站| 99精品国产成人一区二区| 国产精品伦一区二区三级视频| 曰本女人牲交全视频免费播放| 亚洲午夜无码久久yy6080 | 精品黄色一区二区三区| 丰满少妇被啪啪到高潮迷轩| 亚洲av永久无码精品一福利| 亚洲精品乱码久久久久久日本蜜臀| 亚洲产在线精品亚洲第一站一| 亚洲AV无码一区二区水蜜桃| 少妇呻吟一区二区三区| 少妇无码av无码专线区大牛影院| 丰满老熟妇好大bbbbb| 久久精品这里只有精品| 免费在线观看蜜桃视频| 亚洲精品视频中文字幕| 少妇人妻在线无码天堂视频网|