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

    矩量法計算煙幕粒子的紅外消光特性

    2021-03-06 02:47:32劉清海彭文聯(lián)代曉東
    紅外技術(shù) 2021年2期
    關(guān)鍵詞:煙幕旋轉(zhuǎn)體消光

    劉清海,蔣 云,彭文聯(lián),張 彤,代曉東

    〈制導與對抗〉

    矩量法計算煙幕粒子的紅外消光特性

    劉清海,蔣 云,彭文聯(lián),張 彤,代曉東

    (防化研究院,北京 102205)

    基于矩量法建立了旋轉(zhuǎn)體煙幕粒子的紅外消光模型,對圓片、等效球體和等效柱體等形狀的石墨粒子進行了消光性能計算,研究了形狀、尺寸、厚度等參數(shù)對粒子紅外消光能力的影響。計算結(jié)果表明,小尺寸粒子以吸收消光為主,大尺寸粒子以散射消光為主;圓片形狀的粒子具有最佳的紅外消光能力,厚度變小時,圓片粒子對紅外的消光能力顯著升高;半徑在1.5~2.1mm范圍內(nèi)時,100nm厚度的圓片對1~10mm范圍的紅外消光系數(shù)均大于5.0m2/g,在整個紅外波段范圍內(nèi)都表現(xiàn)出良好的消光性能。

    矩量法;旋轉(zhuǎn)體;煙幕粒子;紅外消光

    0 引言

    煙幕材料對電磁波產(chǎn)生衰減的基本原理是電磁波經(jīng)過物質(zhì)顆粒時,與顆粒發(fā)生了散射和吸收等相互作用,電磁波沿原有方向的傳輸能量減弱[1-2]。煙幕粒子對紅外輻射消光能力的強弱取決于顆粒的成分、復折射率、帶電性等自身特性參量[3-5]。建立適用于各種粒子的消光計算模型,可以揭示材料產(chǎn)生遮蔽/干擾效能的內(nèi)在作用機制,加深對粒子與電磁波相互作用原理的理解。同時,如果能夠借助大量的理論模擬和計算,獲得材料的形狀、尺寸、物理特性等基礎參數(shù),并基于最優(yōu)參數(shù)集進行干擾材料設計,將大大加快光電干擾材料的研發(fā)進度,節(jié)約研發(fā)成本。

    當前,可用于煙幕粒子消光性能計算的方法主要有兩類,一類是適用于球形粒子精確求解的Lorenz-Mie理論[6],該理論是最早發(fā)展的光散射模型,能夠?qū)η蛐瘟W拥南馓匦赃M行精確求解。另一類是適用于球形和非球形粒子精確求解的數(shù)值計算方法,此類方法通過在一定邊界條件下直接數(shù)值求解電磁波傳播方程,如Maxwell方程組及Helmholtz方程等,獲得粒子的電磁散射特性。按照計算原理,常見的數(shù)值散射理論及模型大致可分為[7-8]:①基于場展開方式的散射模型,主要包括T矩陣法、擴展邊界條件法、分離變量法和點匹配法等;②基于體積積分方程的散射模型,主要包括矩量法和離散偶極近似法等;③基于微元法的散射模型,主要包括時域有限差分法和有限元法等。每種計算方法都有自己的優(yōu)點和局限性,目前還沒有一種方法能夠精確計算所有類型粒子的電磁波散射特性。因此,在實踐中通常是對一種方法盡可能地優(yōu)化近似以提高其計算準確性,同時也可采用多種計算方法相結(jié)合的方式,互為補充、互相驗證,以獲取更加精確的計算結(jié)果。

    碳基材料煙幕對多波段電磁波有很好的衰減效果,也是目前應用最為廣泛的光電干擾材料之一[9-10]。本文采用計算精度較高的矩量法建立旋轉(zhuǎn)體煙幕粒子的消光計算模型,研究基本參數(shù)對石墨粒子紅外消光性能的影響,對其形狀、尺寸等參數(shù)進行優(yōu)化設計,以進一步提升其消光性能。矩量法離散得到的矩陣方程為滿陣,存儲復雜度和計算復雜度高,是一種比較耗費計算內(nèi)存和時間的數(shù)值方法[11]。旋轉(zhuǎn)體基本囊括了主要的煙幕粒子形狀,如球體、柱體、圓片等,旋轉(zhuǎn)體具有軸對稱性,計算區(qū)域由傳統(tǒng)的三維區(qū)域轉(zhuǎn)換為由母線和旋轉(zhuǎn)軸構(gòu)成的二維區(qū)域,求解區(qū)域和未知量的數(shù)目大大減小,可以節(jié)約計算內(nèi)存和計算時間,實現(xiàn)高速的計算。

    1 消光計算模型

    1.1 計算方法

    旋轉(zhuǎn)體粒子矩量法消光計算模型的基本思想是將粒子表面的未知電流密度積分方程用基函數(shù)和權(quán)函數(shù)離散化成矩陣形式的代數(shù)方程,然后計算矩陣元素并求解該矩陣方程,最后得到表面電磁場和各種衰減截面積[12]。該模型是在Mautz和Harrington[13-14]用矩量法研究旋轉(zhuǎn)體散射特性的基礎上建立的,此后又經(jīng)過不斷的完善和發(fā)展[15-17],現(xiàn)在廣泛用于解決旋轉(zhuǎn)體粒子的光學散射問題,本文只簡要介紹一下其基本原理。

    同樣,權(quán)函數(shù)可定義為:

    于是可得到矩陣方程如下:

    各矩陣元素的內(nèi)積表示形式為:

    由矩陣方程可求得各系數(shù),于是其散射振幅函數(shù)可表示為:

    其中:

    散射體的微分散射截面積可表示為:

    散射截面積為:

    式中:i為入射角;s為散射角;s為方位角。

    當散射方向為入射方向的反方向時為后向散射,此時,i+s=p可由式(18)計算得到后向散射截面積b。

    由前向散射原理,消光截面積為:

    吸收截面積為消光截面積與散射截面積之差:

    對于隨機取向的粒子,其各種衰減截面積的平均值都可用下式求得:

    消光效率因子ext、散射效率因子sca、吸收效率因子abs和后向散射效率因子bac可由相應的截面積除以粒子在電磁波傳播方向的投影面積得到,即:

    ext=e/R(22)

    sca=s/R(23)

    abs=a/R(24)

    bac=b/R(25)

    式中:R為旋轉(zhuǎn)體粒子在電磁波傳播方向的投影面積。

    效率因子均為無量綱的量,消光效率因子、散射效率因子和吸收效率因子之間有如下關(guān)系:

    ext=sca+abs(26)

    質(zhì)量消光系數(shù)是單位質(zhì)量粒子的消光截面,對于單一尺寸的粒子,質(zhì)量消光系數(shù)為:

    e/(27)

    式中:為旋轉(zhuǎn)體煙幕粒子的質(zhì)量。

    對于粒子尺度不均的多分散體系,質(zhì)量消光系數(shù)可用以下公式求得:

    式中:為粒子的尺寸;()為粒子消光截面積;()為粒子數(shù)量分布函數(shù);()為粒子質(zhì)量。

    1.2 模型參數(shù)

    圖1 圓片、球體和柱體粒子示意圖

    利用建立的旋轉(zhuǎn)體消光計算模型對3種形狀粒子進行消光性能計算,計算參數(shù)設置為:

    ①半徑:0.25~10mm;

    ②同形狀所有粒子為單分散體系(尺寸均一),三角展開函數(shù)個數(shù)取值為60~120,方位角展開函數(shù)個數(shù)球體取值1,其他形狀取20~40;

    ③入射波長中波紅外取=4mm,長波紅外取=10mm;

    ④石墨粒子密度取值為2.25g/cm3;

    ⑤復折射率的值取自文獻[18]。

    對于非球形粒子,入射角隨機取向進行角度平均,并且存在電磁極化,其消光性能計算結(jié)果取垂直極化和平行極化的平均值。

    2 消光計算結(jié)果

    2.1 粒子形狀

    圖2是計算得到的不同尺寸的3種形狀粒子對4mm和10mm紅外的消光效率因子ext、散射效率因子sca、吸收效率因子abs和后向散射因子bac。對4mm紅外而言,當圓片半徑<0.75mm、球體半徑<0.46mm、柱體半徑<0.43mm時,3種粒子的abs>sca,即吸收衰減在紅外消光效應中占主導。對10mm紅外而言,當圓片半徑<1.5mm、球體半徑<0.53mm、柱體半徑<0.32mm時,3種粒子的abs>sca,同樣是吸收衰減占主導。因此,對兩個波長的紅外來說,3種粒子在尺寸較小時都以吸收消光為主,尺寸變大后都以散射消光為主。另外,計算結(jié)果表明3種形狀粒子對紅外的后向散射效應都相對較弱。

    圖2 不同形狀的粒子對4mm和10mm紅外的效率因子

    圖3為不同尺寸的粒子對4mm和10mm紅外的消光系數(shù)計算結(jié)果,3種粒子的消光系數(shù)都是先隨半徑增大而增大,到達峰值后變?yōu)殡S著半徑增大而減小,說明粒子對特定波長電磁波的衰減應該有相應的最佳尺寸,此時消光系數(shù)達到最大值。對比3種形狀粒子的消光系數(shù)可以發(fā)現(xiàn),對4mm紅外而言:在所有尺寸下,等效球體的消光性能最差;半徑大于0.5mm時,圓片的消光性能最佳;半徑小于0.5mm時,等效柱體表現(xiàn)出最好的消光性能;3種形狀粒子對4mm紅外消光性能最大時對應的尺寸分別為圓片半徑0.75mm,球體半徑0.53mm,柱體半徑0.08mm。對10mm紅外而言,消光性能基本類似:等效球體在所有尺寸下都表現(xiàn)出最差的消光性能;半徑大于1.5mm時,圓片的消光性能最佳;圓片半徑小于1.5mm時,等效柱體表現(xiàn)出最好的消光性能;3種粒子對10mm紅外消光性能最大時對應的尺寸分別為圓片半徑2.0mm,球體半徑1.6mm,柱體半徑0.19mm。

    圖3 不同尺寸的粒子對4mm和10mm紅外的消光系數(shù)

    總體來看,無論是對中波還是長波紅外,球體的消光系數(shù)在三者中都是最小的,對10mm紅外,幾乎所有尺寸的消光系數(shù)都小于0.5m2/g,消光效率很低。實際上煙幕材料制備加工過程中,球體是較為容易獲取的,但從材料消光效率的角度來看,應盡可能避免將煙幕粒子設計成球形。非球形的圓片和柱體粒子對中波和長波紅外都表現(xiàn)出了良好的消光能力,特別是圓片形粒子半徑處于0.25~2.0mm之間時,對4mm紅外的消光系數(shù)均大于2.0m2/g,半徑處于1.4~2.4mm之間時,對10mm紅外的消光系數(shù)均大于2.0m2/g。因此,調(diào)控設計半徑在1.4~2.0mm的區(qū)間時,粒子對中波和長波紅外的消光系數(shù)均可超過2.0m2/g。柱體粒子雖然在中波和長波紅外也表現(xiàn)出了較好的消光能力,尤其在尺寸較小的時候甚至優(yōu)于圓片,但柱體粒子加工的技術(shù)難度大、成本高,且半徑大于1mm后其消光能力將顯著下降,因此綜合考慮平均消光能力和制造加工成本,圓片狀結(jié)構(gòu)將是石墨類紅外消光粒子的最優(yōu)選擇。

    2.2 圓片厚度

    圖4是不同厚度的圓片粒子對4mm和10mm紅外的消光系數(shù)計算結(jié)果,計算參數(shù)取值為100,取值為20。隨著厚度的減小,5種半徑粒子的消光系數(shù)均顯著升高。半徑2.0~5.0mm的4種粒子在片層厚度100nm時,對兩個波長紅外的消光系數(shù)均可達到4.0m2/g以上,厚度為50nm時可進一步提升到6.0m2/g以上,在厚度小于50nm的納米尺度,消光系數(shù)會隨著厚度的變薄而迅速增加。這一結(jié)果表明,圓片粒子對紅外的消光系數(shù)和消光能力會隨著厚度的減小而迅速提升。

    圖4 不同厚度圓片粒子對4mm和10mm紅外的消光系數(shù)

    為進一步研究厚度對圓片消光性能的影響,將圓片粒子的半徑固定在2mm,控制厚度在20~200nm之間變化,計算不同厚度的圓片對1~10mm各波長紅外的消光系數(shù),計算參數(shù)取值為60,取值為10,計算結(jié)果見圖5。在計算波長范圍內(nèi),圓片粒子對紅外的消光系數(shù)都隨著厚度的減小而增加,這進一步說明厚度對圓片粒子的紅外消光性能具有極為重要的影響。表1為半徑2mm圓片粒子對1~10mm紅外的平均消光系數(shù),20nm厚度圓片的平均消光系數(shù)達到了13.2m2/g。但隨著厚度的增加,平均消光系數(shù)呈現(xiàn)下降趨勢,當厚度增加到200nm時,平均消光系數(shù)下降到3.5m2/g。厚度為20nm的圓片,對近紅外和中紅外的消光系數(shù)最大,但往遠紅外延伸時消光系數(shù)卻出現(xiàn)了下降趨勢,對10mm紅外的消光系數(shù)已經(jīng)與40nm厚度的圓片接近,說明過小的厚度將不利于圓片粒子對遠紅外的消光能力。因此,在進行圓片粒子設計時應盡可能減小厚度,但為了兼顧中遠紅外的消光能力,其厚度不能低于某一范圍。對于半徑2mm圓片粒子來說,其厚度不應低于20nm,否則其對遠紅外的消光能力將大幅降低。

    2.3 薄層圓片半徑

    圓片粒子的厚度達到納米級后,其消光性能會有顯著的提升,為進一步研究納米級的圓片半徑對其消光性能的影響,固定圓片的厚度為100nm,通過調(diào)整其半徑進行圓片紅外消光性能計算,計算結(jié)果仍然取垂直極化和平行極化的平均值,計算結(jié)果如圖6所示。對每個波長的紅外而言,圓片粒子都存在一個最佳半徑,尺寸在該半徑值時圓片粒子對特定波長紅外的消光系數(shù)最大,圓片對各個波長紅外的最大消光系數(shù)和對應的最佳半徑如表2所示。每個波長下圓片的最佳半徑隨著波長的增大而增大,同時最佳直徑略小于相對應的波長的一半。另外,最大消光系數(shù)隨著紅外波長的增大反而減小,說明圓片對遠紅外的消光能力隨著波長的增加而逐漸降低。圖6表明,當半徑大于1.2mm時圓片對所有波長紅外的消光系數(shù)都大于3.0m2/g,當半徑大于1.5mm時圓片對所有波長紅外的消光系數(shù)都大于5.0m2/g,但當半徑超過2.1mm以后,圓片對所有波長紅外的消光系數(shù)均有隨半徑增大而變小的趨勢。因此,對于100nm厚度的圓片而言,如果將其半徑設計在1.5~2.1mm之間,則可使其對1~10mm紅外均表現(xiàn)出良好的消光性能。

    表1 半徑2mm的圓片粒子對1~10mm紅外的平均消光系數(shù)

    圖5 半徑2mm圓片粒子對1~10mm紅外的消光系數(shù)

    圖6 100nm厚圓片半徑對消光系數(shù)的影響

    表2 厚度為100nm的圓片粒子對1~10mm紅外的最大消光系數(shù)和對應的最佳半徑

    3 結(jié)論

    基于矩量法建立了旋轉(zhuǎn)體煙幕粒子的紅外消光模型,并對圓片、球體和柱體等基本粒子進行了消光性能計算。計算結(jié)果表明,3種基本形狀的粒子,在尺寸較小時都以吸收消光為主,在尺寸較大時都以散射消光為主。相對來說,球體粒子的消光能力最差,圓片粒子具有最佳的紅外消光能力。當半徑在1.4~2.0mm區(qū)間時,圓片對4mm和10mm紅外的消光系數(shù)均可超過2.0m2/g。通過減小圓片的厚度,可以顯著提升圓片粒子對紅外的消光能力,但考慮到遠紅外波段的消光效果,其厚度不能低于20nm。對于100nm厚度的圓片而言,當其半徑在1.5~2.1mm之間時,粒子對1~10mm紅外的消光系數(shù)均大于5.0m2/g,可使該尺寸的圓片粒子在整個紅外波段范圍內(nèi)都表現(xiàn)出良好的消光性能。石墨是當前較為高效的煙幕干擾材料之一,本計算結(jié)果可為進一步提高石墨干擾能力提供指導,也可為其他先進光電干擾材料的設計和研制提供理論基礎。

    [1] 霸書紅, 陳永進, 沙育林, 等. 抗紅外煙幕中固體消光材料的研究進展[J]. 含能材料, 2018, 26(4): 364-372.

    BA Shuhong, CHEN Yongjin, SHA Yulin, et al. Research progress of solid extinction materials in anti-infrared smokescreen[J]., 2018, 26(4): 364-372.

    [2] 丁國振, 張占月, 周思引, 等. 空間煙幕的擴散機理及衰減性能研究[J]. 紅外技術(shù), 2014, 36(11):914-919.

    DING Guozhen, ZHANG Zhanyue, ZHOU Siyin, et al. Study on diffusion mechanism and attenuation performance of space smoke -screen[J]., 2014, 36(11): 914-919.

    [3] 陳芳芳, 耿蕊, 呂勇. 激光紅外大氣傳輸透過率模型研究[J]. 紅外技術(shù), 2015, 37(6): 496-501.

    CHEN Fangfang, GENG Rui, LV Yong. Research on the transmittance model of laser infrared atmospheric transmission[J]., 2015, 37(6): 496-501.

    [4] 王紅霞, 孫紅輝, 宋仔標, 等. 基于蒙特卡羅方法的煙幕透過率計算與分析[J]. 紅外與激光工程, 2013, 41(5):1200-1205.

    WANG Hongxia, SUN Honghui, SONG Zibiao, et al. Numerical calculation and analysis of transmittance of smoke screen based on Monte Carlo method[J]., 2013, 41(5): 1200-1205.

    [5] 胡以華, 黃寶錕, 顧有林, 等. 生物顆粒遠紅外波段平均消光效率因子模型構(gòu)建[J]. 紅外與激光工程, 2018, 47(10): 1004003.

    HU Yihua, HUANG Baokun, GU Youlin, et al. Model construction of biological particles average extinction efficiency factor in far infrared band [J]., 2018, 47(10): 1004003.

    [6] Bohren C F, Huffman D R.[M]. New York: John Wiley & Sons, 1983.

    [7] 許麗生, 陳洪濱, 丁繼烈. 非球形粒子光散射計算研究的進展綜述[J]. 地球科學進展, 2014, 29(8): 903-912.

    XU Lisheng, CHEN Hongbin, DING Jilie. An overview of the advances in computational studies on light scattering by nonspherical particles[J]., 2014, 29(8): 903-912.

    [8] 高太長, 胡帥, 李浩. 氣溶膠光散射數(shù)值模擬的研究現(xiàn)狀及進展[J]. 氣象科學, 2017, 37(5): 598-609.

    GAO Taichang, HU Shuai, LI Hao. Actuality and prospect of aerosol scattering numerical simulation techniques[J]., 2017, 37(5): 598-609.

    [9] 劉清海, 劉海鋒, 代曉東, 等. 石墨烯煙幕紅外干擾性能研究[J]. 紅外技術(shù), 2019, 41(11): 1071-1076.

    LIU Qinghai, LIU Haifeng, DAI Xiaodong, et al. Infrared interfering performance of graphene smoke screen[J]., 2019, 41(11): 1071-1076.

    [10] 郭曉鐺, 喬小晶, 李旺昌, 等. 鐵磁體/碳復合材料多頻干擾性能[J]. 紅外與激光工程, 2016, 47(10): 0321001.

    GUO Xiaodang, QIAO Xiaojing, LI Wangchang, et al. Multi-frequency jamming of ferromagnet/carbon composite[J]., 2016, 47(10): 0321001.

    [11] Mishchenko M I, Travis L D, Lacis A A.,,[M]. New York: Cambridge, 2002.

    [12] Mautz J, Harrington R. Electromagnetic scattering from a homogeneous body of revolution[J]., 1979, 33: 71-80.

    [13] Mautz J, Harrington R. H-field, E-field, and combined-field solutions for conducting bodies of revolution[J]., 1978, 32: 157-164.

    [14] Mautz J, Harrington R. A combined-source solution for radiation and scattering from a perfectly conducting body[J]., 1979, 27(4): 445-454.

    [15] Davidson D, Mc Namara D. Predicting radiation patterns from aperture antennas on structures using the method of moments body of revolution technique[C], 1987, 25: 25-30.

    [16] Glisson A W, Wilton D R. Simple and efficient numerical methods for problems of electromagnetic radiation and scattering from surfaces[J]., 1980, 28: 593-603.

    [17] Medgyesi-Mitschg L N, Putnam J M. Electromagnetic scattering from axially inhomogeneous bodies of revolution[J]., 1984, 32: 797-806.

    [18] Aleksandra B Djuri?i?, E Herbert Li. Optical properties of graphite[J]., 1999, 85(10): 7404-7410.

    Infrared Extinction Calculations of Smokescreen Particles by Moment Method

    LIU Qinghai,JIANG Yun,PENG Wenlian,ZHANG Tong,DAI Xiaodong

    (,102205,)

    A moment-method-based infrared extinction model of rotating smokescreen particles is applied to perform calculations entailing three graphite particles including flakes, spheres, and cylinders, mainly focusing on the relationship between extinction ability and particle parameters, such as shape, diameter, and thickness. The calculations suggest that extinction is attributed to absorption when the particle size is small and attributed to scattering when the particle size is large. Flakes exhibit the best infrared extinction performance. When flakes become thinner, their extinction abilities are enhanced. Flakes with 100nm thickness and 1.5–2.1mm radius exhibit outstanding extinction performance in the 1–10mm infrared wavelength range, with an average infrared extinction coefficient as high as 5.0m2/g.

    moment method, rotating body, smoke particles, infrared extinction

    TN21

    A

    1001-8891(2021)02-0138-07

    2020-04-14;

    2020-09-08.

    劉清海(1984-),男,山東濰坊人,助理研究員,博士,主要從事光電干擾理論與技術(shù)研究。E-mail:liuqh_pub@163.com。

    代曉東(1974-),男,重慶人,研究員,博士,主要從事光電功能材料與特種裝備研究。E-mail:dxd010@163.com。

    國防預研項目。

    猜你喜歡
    煙幕旋轉(zhuǎn)體消光
    低聚合度消光PVC樹脂的研制
    中國氯堿(2023年1期)2023-02-22 06:34:34
    固化促進劑對聚酯/環(huán)氧粉末涂料消光性能的影響研究
    煙幕遮蔽載荷的性能測試
    煙幕干擾對紅外成像導引頭命中概率影響
    激光與紅外(2021年3期)2021-04-07 12:22:08
    基于柱殼法及柱坐標系求解旋轉(zhuǎn)體的體積
    旋轉(zhuǎn)體體積的探討
    基于微元法的旋轉(zhuǎn)體體積問題的一題多解
    旋轉(zhuǎn)體容球的一個有趣共性再探究
    考試周刊(2016年82期)2016-11-01 13:00:24
    戰(zhàn)場上的魔術(shù)師
    BOPP消光膜及其生產(chǎn)工藝
    河南化工(2011年3期)2011-09-24 03:57:30
    亚洲国产毛片av蜜桃av| 国产亚洲精品第一综合不卡| 亚洲专区字幕在线| 亚洲 欧美一区二区三区| 在线观看免费高清a一片| 欧美乱码精品一区二区三区| 黄色a级毛片大全视频| 校园春色视频在线观看| 久久ye,这里只有精品| 午夜免费鲁丝| 精品少妇久久久久久888优播| 黑人操中国人逼视频| 99国产精品一区二区三区| 一级黄色大片毛片| 久久香蕉精品热| 久久久久久久国产电影| tube8黄色片| 精品国产美女av久久久久小说| a级片在线免费高清观看视频| 午夜福利免费观看在线| 夜夜爽天天搞| 欧美日韩一级在线毛片| 国产精品亚洲av一区麻豆| 久久国产精品男人的天堂亚洲| 美女午夜性视频免费| 人人妻人人爽人人添夜夜欢视频| 午夜福利在线免费观看网站| 日韩视频一区二区在线观看| 国产精品久久久久久人妻精品电影| 男人的好看免费观看在线视频 | 欧美激情 高清一区二区三区| 成年动漫av网址| 久久精品国产亚洲av香蕉五月 | 两性夫妻黄色片| 国产熟女午夜一区二区三区| 欧美丝袜亚洲另类 | 黄色a级毛片大全视频| av电影中文网址| 热re99久久精品国产66热6| 国产男靠女视频免费网站| 国产不卡av网站在线观看| 国产一区二区三区综合在线观看| 三级毛片av免费| 亚洲精品av麻豆狂野| 12—13女人毛片做爰片一| 99国产精品免费福利视频| 91国产中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产精品亚洲一级av第二区| 久久久久精品人妻al黑| 久久久国产成人免费| 午夜免费鲁丝| 国产在线一区二区三区精| 黄色视频不卡| 女人精品久久久久毛片| www日本在线高清视频| 日韩成人在线观看一区二区三区| 满18在线观看网站| 99久久综合精品五月天人人| 9色porny在线观看| 国产欧美日韩一区二区精品| 宅男免费午夜| 亚洲精品国产精品久久久不卡| 国产欧美日韩一区二区三| 国产精品九九99| 国产精品久久久久久人妻精品电影| 久久久国产成人免费| 91麻豆av在线| 久久 成人 亚洲| 国产乱人伦免费视频| www.自偷自拍.com| av视频免费观看在线观看| 国产精品一区二区在线观看99| a级毛片在线看网站| 免费在线观看亚洲国产| 国产aⅴ精品一区二区三区波| 亚洲午夜理论影院| 亚洲成人免费电影在线观看| 91成人精品电影| 在线观看一区二区三区激情| 国产精品免费视频内射| 老熟女久久久| 99热只有精品国产| 午夜亚洲福利在线播放| tocl精华| 亚洲熟女精品中文字幕| xxx96com| 国产精品98久久久久久宅男小说| 天天影视国产精品| 精品乱码久久久久久99久播| av网站免费在线观看视频| 性少妇av在线| 国产aⅴ精品一区二区三区波| netflix在线观看网站| 亚洲成人免费电影在线观看| 亚洲视频免费观看视频| 国产精品久久久久久人妻精品电影| 中出人妻视频一区二区| 王馨瑶露胸无遮挡在线观看| 国产精品综合久久久久久久免费 | 精品乱码久久久久久99久播| 丰满的人妻完整版| 视频区欧美日本亚洲| 中亚洲国语对白在线视频| 两人在一起打扑克的视频| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品合色在线| 最近最新中文字幕大全免费视频| 成人黄色视频免费在线看| 首页视频小说图片口味搜索| 成人精品一区二区免费| 妹子高潮喷水视频| 午夜免费成人在线视频| 少妇被粗大的猛进出69影院| 久久久久精品人妻al黑| 啪啪无遮挡十八禁网站| 女警被强在线播放| 久久午夜综合久久蜜桃| www日本在线高清视频| 国产欧美日韩一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| www.熟女人妻精品国产| 国产精品一区二区免费欧美| 91九色精品人成在线观看| 国产精品免费视频内射| 三上悠亚av全集在线观看| 日韩 欧美 亚洲 中文字幕| 91成年电影在线观看| 国产精品成人在线| 国产男女内射视频| 精品久久久精品久久久| 超碰97精品在线观看| 久久狼人影院| 在线免费观看的www视频| 国产不卡av网站在线观看| 国产一区二区三区综合在线观看| 啦啦啦 在线观看视频| 黄片大片在线免费观看| 国内久久婷婷六月综合欲色啪| 一级a爱视频在线免费观看| 精品国产美女av久久久久小说| 日韩大码丰满熟妇| 中出人妻视频一区二区| 美女视频免费永久观看网站| 亚洲av日韩在线播放| www.熟女人妻精品国产| 亚洲伊人色综图| 精品国产一区二区三区久久久樱花| 久久性视频一级片| 可以免费在线观看a视频的电影网站| 亚洲专区字幕在线| 深夜精品福利| 精品欧美一区二区三区在线| 欧美亚洲日本最大视频资源| 大型av网站在线播放| 久久国产精品人妻蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品亚洲av国产电影网| 国产精品久久久久成人av| 国产免费av片在线观看野外av| 一区在线观看完整版| 精品卡一卡二卡四卡免费| 国产又色又爽无遮挡免费看| 99香蕉大伊视频| 亚洲专区字幕在线| 99在线人妻在线中文字幕 | 悠悠久久av| 极品教师在线免费播放| 黄色丝袜av网址大全| 成年版毛片免费区| 亚洲一区二区三区不卡视频| 黄色片一级片一级黄色片| 国产日韩一区二区三区精品不卡| 怎么达到女性高潮| 美女视频免费永久观看网站| 亚洲精品美女久久av网站| 精品熟女少妇八av免费久了| 久久热在线av| 国产精品99久久99久久久不卡| 99re在线观看精品视频| a级片在线免费高清观看视频| 女性生殖器流出的白浆| 亚洲五月天丁香| 成熟少妇高潮喷水视频| av不卡在线播放| 精品久久久久久,| 极品教师在线免费播放| 精品无人区乱码1区二区| 91精品国产国语对白视频| 亚洲精品国产区一区二| 久久精品国产清高在天天线| 久久天躁狠狠躁夜夜2o2o| 校园春色视频在线观看| 97人妻天天添夜夜摸| 精品视频人人做人人爽| 夜夜爽天天搞| av一本久久久久| 变态另类成人亚洲欧美熟女 | 国产精品亚洲一级av第二区| 精品人妻在线不人妻| 人成视频在线观看免费观看| 国产精品永久免费网站| 一边摸一边做爽爽视频免费| av中文乱码字幕在线| 国产免费现黄频在线看| 日本撒尿小便嘘嘘汇集6| 亚洲色图综合在线观看| 国产精品1区2区在线观看. | 日韩欧美一区二区三区在线观看 | 中文字幕高清在线视频| 亚洲成人免费av在线播放| 交换朋友夫妻互换小说| 久久午夜综合久久蜜桃| 国产亚洲一区二区精品| 欧美黑人欧美精品刺激| 大香蕉久久成人网| 我的亚洲天堂| 国产亚洲一区二区精品| 国产男靠女视频免费网站| 亚洲七黄色美女视频| 国产精品免费大片| 亚洲成人免费电影在线观看| 日本黄色日本黄色录像| 人人妻人人澡人人看| 日韩免费高清中文字幕av| 亚洲第一av免费看| 欧美日韩成人在线一区二区| 啦啦啦在线免费观看视频4| 嫩草影视91久久| 亚洲中文av在线| 日韩大码丰满熟妇| 国产日韩一区二区三区精品不卡| 欧美性长视频在线观看| 亚洲美女黄片视频| 国产成人精品无人区| 免费久久久久久久精品成人欧美视频| 亚洲一区中文字幕在线| 美女午夜性视频免费| 18禁裸乳无遮挡免费网站照片 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美另类亚洲清纯唯美| 深夜精品福利| 国产精华一区二区三区| 黄色丝袜av网址大全| 999精品在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 韩国精品一区二区三区| 亚洲精品国产区一区二| 欧美老熟妇乱子伦牲交| 国产成人精品无人区| 两个人看的免费小视频| 大型黄色视频在线免费观看| 亚洲性夜色夜夜综合| 国产精品秋霞免费鲁丝片| 搡老乐熟女国产| 亚洲精品一卡2卡三卡4卡5卡| 国产一卡二卡三卡精品| 十八禁人妻一区二区| 淫妇啪啪啪对白视频| 国产成人欧美| 亚洲久久久国产精品| 亚洲人成电影观看| 高清av免费在线| 免费看十八禁软件| www.精华液| 国产亚洲精品久久久久久毛片 | 午夜亚洲福利在线播放| 麻豆国产av国片精品| 午夜免费鲁丝| 黑人巨大精品欧美一区二区蜜桃| 少妇裸体淫交视频免费看高清 | 国产日韩欧美亚洲二区| a级毛片黄视频| 搡老岳熟女国产| 亚洲五月天丁香| 多毛熟女@视频| 激情视频va一区二区三区| 免费高清在线观看日韩| 久久久久精品国产欧美久久久| 久久人人97超碰香蕉20202| 亚洲av欧美aⅴ国产| 欧洲精品卡2卡3卡4卡5卡区| 纯流量卡能插随身wifi吗| 在线看a的网站| 成熟少妇高潮喷水视频| 免费日韩欧美在线观看| 久久国产乱子伦精品免费另类| 亚洲aⅴ乱码一区二区在线播放 | 五月开心婷婷网| 久久ye,这里只有精品| 国产xxxxx性猛交| 国产成人一区二区三区免费视频网站| 日韩免费av在线播放| 91国产中文字幕| 亚洲欧美精品综合一区二区三区| 99国产精品一区二区蜜桃av | 久久精品成人免费网站| 国产精品1区2区在线观看. | 亚洲自偷自拍图片 自拍| 女性生殖器流出的白浆| 精品一区二区三区av网在线观看| 美女高潮喷水抽搐中文字幕| 国产在线精品亚洲第一网站| 一进一出好大好爽视频| 在线播放国产精品三级| 久久精品aⅴ一区二区三区四区| 女人被狂操c到高潮| 男人舔女人的私密视频| bbb黄色大片| 丝袜在线中文字幕| 久久精品亚洲av国产电影网| 动漫黄色视频在线观看| 久久午夜亚洲精品久久| 婷婷成人精品国产| 男人操女人黄网站| 飞空精品影院首页| 精品熟女少妇八av免费久了| 每晚都被弄得嗷嗷叫到高潮| 69精品国产乱码久久久| 中文亚洲av片在线观看爽 | 国产激情久久老熟女| 日韩有码中文字幕| 丁香六月欧美| 午夜免费成人在线视频| 亚洲aⅴ乱码一区二区在线播放 | 一二三四社区在线视频社区8| 国产成人影院久久av| xxxhd国产人妻xxx| 精品亚洲成a人片在线观看| 这个男人来自地球电影免费观看| 日日摸夜夜添夜夜添小说| 国产又色又爽无遮挡免费看| 成年人黄色毛片网站| 曰老女人黄片| а√天堂www在线а√下载 | 两性夫妻黄色片| 国产亚洲精品久久久久久毛片 | 午夜精品国产一区二区电影| 99精品久久久久人妻精品| 动漫黄色视频在线观看| 桃红色精品国产亚洲av| 超碰97精品在线观看| 欧美精品人与动牲交sv欧美| 91在线观看av| 国产有黄有色有爽视频| 热99国产精品久久久久久7| 叶爱在线成人免费视频播放| 黄色毛片三级朝国网站| 日韩欧美一区二区三区在线观看 | 国产亚洲精品第一综合不卡| 国产成人欧美在线观看 | 中文字幕最新亚洲高清| 久久中文看片网| www.精华液| 国产又色又爽无遮挡免费看| 亚洲成国产人片在线观看| 十八禁高潮呻吟视频| 极品人妻少妇av视频| 久久精品国产a三级三级三级| 国产成人一区二区三区免费视频网站| 十八禁高潮呻吟视频| 国产麻豆69| 美女视频免费永久观看网站| 精品免费久久久久久久清纯 | 成人18禁在线播放| 国产精品一区二区精品视频观看| 精品乱码久久久久久99久播| 日韩一卡2卡3卡4卡2021年| 99国产精品一区二区三区| 国产免费av片在线观看野外av| 国产成人免费观看mmmm| 午夜91福利影院| 一个人免费在线观看的高清视频| 国产成人精品久久二区二区免费| 国产精品自产拍在线观看55亚洲 | 免费女性裸体啪啪无遮挡网站| 中文字幕精品免费在线观看视频| 久久久久久久精品吃奶| 又紧又爽又黄一区二区| 中文欧美无线码| 国产精品久久久久久精品古装| 如日韩欧美国产精品一区二区三区| 丰满饥渴人妻一区二区三| 嫁个100分男人电影在线观看| 一边摸一边做爽爽视频免费| 男人操女人黄网站| 伊人久久大香线蕉亚洲五| 国产精品秋霞免费鲁丝片| 人人妻,人人澡人人爽秒播| 日本撒尿小便嘘嘘汇集6| 亚洲美女黄片视频| 免费久久久久久久精品成人欧美视频| 久久天堂一区二区三区四区| 亚洲精品久久午夜乱码| 女人久久www免费人成看片| 看片在线看免费视频| 国产精品.久久久| 久久 成人 亚洲| 一个人免费在线观看的高清视频| 亚洲自偷自拍图片 自拍| 久久99一区二区三区| 国产乱人伦免费视频| 亚洲精品自拍成人| 美女视频免费永久观看网站| 久久这里只有精品19| 日韩免费av在线播放| 精品一区二区三区四区五区乱码| 亚洲欧美色中文字幕在线| av网站在线播放免费| 久久亚洲精品不卡| 久久人妻av系列| 捣出白浆h1v1| 99久久99久久久精品蜜桃| 99国产综合亚洲精品| 99国产精品一区二区三区| 精品福利永久在线观看| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区精品| av天堂在线播放| 久久久国产成人精品二区 | 久久99一区二区三区| 精品福利观看| 大香蕉久久成人网| 男人舔女人的私密视频| 久久精品国产99精品国产亚洲性色 | 国产精品久久电影中文字幕 | 国产亚洲欧美精品永久| 宅男免费午夜| 成人18禁在线播放| 国产xxxxx性猛交| 91av网站免费观看| 99精品久久久久人妻精品| 亚洲中文字幕日韩| 国产成人精品在线电影| 好看av亚洲va欧美ⅴa在| 欧美亚洲 丝袜 人妻 在线| 大片电影免费在线观看免费| 久久人妻av系列| 99国产极品粉嫩在线观看| 桃红色精品国产亚洲av| 欧美乱码精品一区二区三区| 久久久精品免费免费高清| x7x7x7水蜜桃| 亚洲精品国产区一区二| 视频在线观看一区二区三区| 免费久久久久久久精品成人欧美视频| 欧美人与性动交α欧美软件| 俄罗斯特黄特色一大片| 精品卡一卡二卡四卡免费| 天堂动漫精品| 国产激情久久老熟女| 日韩成人在线观看一区二区三区| 国产成人av教育| 国产精品久久久久久精品古装| 久久久久久久久久久久大奶| 下体分泌物呈黄色| 亚洲国产欧美一区二区综合| 欧美不卡视频在线免费观看 | 波多野结衣一区麻豆| 妹子高潮喷水视频| 亚洲五月婷婷丁香| 欧美日韩精品网址| 精品久久久久久久毛片微露脸| 黄色a级毛片大全视频| 1024香蕉在线观看| 欧美精品av麻豆av| 国产欧美日韩一区二区三| 男女高潮啪啪啪动态图| 99国产极品粉嫩在线观看| 制服诱惑二区| 中文亚洲av片在线观看爽 | 国产又爽黄色视频| 亚洲一卡2卡3卡4卡5卡精品中文| 激情在线观看视频在线高清 | 久久国产精品人妻蜜桃| 欧美日韩精品网址| 中文字幕av电影在线播放| 日本欧美视频一区| 99香蕉大伊视频| 国产免费男女视频| 色94色欧美一区二区| 日日摸夜夜添夜夜添小说| 91精品三级在线观看| 大片电影免费在线观看免费| 大型黄色视频在线免费观看| 欧美日韩一级在线毛片| 成人特级黄色片久久久久久久| 久久久久久久精品吃奶| 精品卡一卡二卡四卡免费| 欧美日韩亚洲高清精品| av在线播放免费不卡| 精品一区二区三区四区五区乱码| 久久国产精品影院| 国产成人精品久久二区二区91| 欧美乱色亚洲激情| 日本黄色日本黄色录像| 国产激情久久老熟女| 搡老岳熟女国产| 久久久国产成人精品二区 | 一区二区三区激情视频| 99riav亚洲国产免费| 变态另类成人亚洲欧美熟女 | 999久久久国产精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 老汉色∧v一级毛片| 国产成人影院久久av| 午夜两性在线视频| 成人手机av| 不卡一级毛片| 欧美日韩一级在线毛片| ponron亚洲| 无人区码免费观看不卡| 国产在线一区二区三区精| 久99久视频精品免费| 999久久久国产精品视频| 别揉我奶头~嗯~啊~动态视频| 亚洲精品国产区一区二| 青草久久国产| 亚洲精品一卡2卡三卡4卡5卡| 男女下面插进去视频免费观看| 国产精品久久久久成人av| 黑丝袜美女国产一区| 亚洲精品av麻豆狂野| 久久久久久久午夜电影 | 亚洲av片天天在线观看| 老司机影院毛片| 中文字幕高清在线视频| 色94色欧美一区二区| 亚洲av成人一区二区三| 俄罗斯特黄特色一大片| 超碰成人久久| 亚洲国产精品一区二区三区在线| 日日爽夜夜爽网站| 精品少妇一区二区三区视频日本电影| 在线观看午夜福利视频| 日韩中文字幕欧美一区二区| 久久香蕉激情| 12—13女人毛片做爰片一| 一进一出抽搐动态| 纯流量卡能插随身wifi吗| 老鸭窝网址在线观看| 女警被强在线播放| 男人舔女人的私密视频| 日本vs欧美在线观看视频| 免费看a级黄色片| e午夜精品久久久久久久| 亚洲av成人不卡在线观看播放网| 中文字幕人妻熟女乱码| 成人国产一区最新在线观看| 国产国语露脸激情在线看| 国产区一区二久久| 91成人精品电影| 久久精品国产99精品国产亚洲性色 | 午夜精品国产一区二区电影| 国产亚洲欧美98| av欧美777| 久久久久精品国产欧美久久久| 在线永久观看黄色视频| 久久久国产成人精品二区 | 国产激情久久老熟女| 一区福利在线观看| 欧美国产精品一级二级三级| 另类亚洲欧美激情| 国产片内射在线| av免费在线观看网站| 91av网站免费观看| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| 下体分泌物呈黄色| 午夜福利在线观看吧| 一二三四在线观看免费中文在| 亚洲五月婷婷丁香| 欧美日韩亚洲高清精品| 日本vs欧美在线观看视频| 国产一区有黄有色的免费视频| 日韩欧美在线二视频 | 精品久久久久久电影网| 亚洲国产毛片av蜜桃av| 日韩一卡2卡3卡4卡2021年| 国产成人av激情在线播放| 日韩三级视频一区二区三区| 美女福利国产在线| 天天添夜夜摸| 亚洲自偷自拍图片 自拍| 亚洲黑人精品在线| 欧美成人午夜精品| 国产精品成人在线| 三上悠亚av全集在线观看| 精品福利永久在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产又爽黄色视频| 亚洲欧洲精品一区二区精品久久久| 国产不卡av网站在线观看| 久久午夜综合久久蜜桃| 捣出白浆h1v1| 在线观看免费高清a一片| 国产精品一区二区在线观看99| 日韩免费高清中文字幕av| 大型av网站在线播放| 亚洲精品在线美女| 久99久视频精品免费| 亚洲成人免费av在线播放| 国产精品免费大片| 欧美乱码精品一区二区三区| 精品人妻熟女毛片av久久网站| 看黄色毛片网站| 18禁美女被吸乳视频| 亚洲精品成人av观看孕妇| 欧美乱妇无乱码| 日韩制服丝袜自拍偷拍| 久久精品国产99精品国产亚洲性色 | 国产区一区二久久| avwww免费| 亚洲精品av麻豆狂野| 久久精品国产亚洲av高清一级|