• <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
    一级毛片aaaaaa免费看小| 免费在线观看成人毛片| av黄色大香蕉| 免费大片黄手机在线观看| 日韩av不卡免费在线播放| 亚洲自偷自拍三级| av女优亚洲男人天堂| 亚洲美女搞黄在线观看| 日本黄色日本黄色录像| 国产精品福利在线免费观看| 激情五月婷婷亚洲| 51国产日韩欧美| 免费人成在线观看视频色| 中文字幕人妻丝袜制服| 亚洲av福利一区| 一本久久精品| 国产成人精品一,二区| 国产男女超爽视频在线观看| 亚洲国产精品专区欧美| 一本大道久久a久久精品| 高清在线视频一区二区三区| 中文字幕人妻丝袜制服| 国产亚洲91精品色在线| 高清不卡的av网站| 高清在线视频一区二区三区| 永久网站在线| 国产欧美另类精品又又久久亚洲欧美| 下体分泌物呈黄色| 免费看不卡的av| 久久免费观看电影| 国产日韩一区二区三区精品不卡 | 免费观看性生交大片5| 日韩一区二区三区影片| 国产精品99久久久久久久久| 国产一区二区在线观看日韩| 午夜91福利影院| 国产深夜福利视频在线观看| 夜夜看夜夜爽夜夜摸| 九九久久精品国产亚洲av麻豆| 中文字幕av电影在线播放| 丝袜在线中文字幕| 色94色欧美一区二区| 国产极品粉嫩免费观看在线 | 99久久精品国产国产毛片| 涩涩av久久男人的天堂| 永久网站在线| 午夜久久久在线观看| 中文乱码字字幕精品一区二区三区| 亚洲国产毛片av蜜桃av| 成人毛片a级毛片在线播放| 中文精品一卡2卡3卡4更新| 春色校园在线视频观看| 亚洲丝袜综合中文字幕| 亚洲真实伦在线观看| 久久影院123| 国产日韩欧美视频二区| 亚洲欧美日韩东京热| 国产真实伦视频高清在线观看| 国产日韩一区二区三区精品不卡 | 黑人巨大精品欧美一区二区蜜桃 | 99国产精品免费福利视频| 国产精品福利在线免费观看| 卡戴珊不雅视频在线播放| 午夜精品国产一区二区电影| 九九在线视频观看精品| 日韩av免费高清视频| 老司机影院成人| 日本欧美视频一区| 国产极品天堂在线| 国产av国产精品国产| 中文字幕亚洲精品专区| 在线天堂最新版资源| 18禁裸乳无遮挡动漫免费视频| 精品久久久精品久久久| 久久国产精品男人的天堂亚洲 | 成人亚洲欧美一区二区av| 日韩制服骚丝袜av| 99久国产av精品国产电影| 精品久久久噜噜| 人妻夜夜爽99麻豆av| 精品久久久久久久久av| 热re99久久精品国产66热6| 亚洲久久久国产精品| 日本色播在线视频| 少妇裸体淫交视频免费看高清| 高清av免费在线| 中文在线观看免费www的网站| 熟妇人妻不卡中文字幕| 日韩电影二区| 婷婷色麻豆天堂久久| 久久久久人妻精品一区果冻| 美女视频免费永久观看网站| 夫妻性生交免费视频一级片| 99久久中文字幕三级久久日本| 美女cb高潮喷水在线观看| 免费少妇av软件| 欧美精品一区二区大全| 9色porny在线观看| 久久精品国产鲁丝片午夜精品| 成人亚洲欧美一区二区av| 国产精品国产三级国产av玫瑰| 高清欧美精品videossex| 国产亚洲av片在线观看秒播厂| 日本免费在线观看一区| 丰满迷人的少妇在线观看| 亚洲成人一二三区av| 成人亚洲精品一区在线观看| 一级毛片 在线播放| 亚洲综合精品二区| kizo精华| 国产一区二区在线观看av| 国模一区二区三区四区视频| 国产在线男女| 久久久久精品久久久久真实原创| 免费观看的影片在线观看| 精品视频人人做人人爽| 哪个播放器可以免费观看大片| 大片免费播放器 马上看| 国产高清三级在线| 91aial.com中文字幕在线观看| 久久6这里有精品| 精品亚洲乱码少妇综合久久| 2022亚洲国产成人精品| 人人妻人人添人人爽欧美一区卜| 亚洲图色成人| 一级,二级,三级黄色视频| 国产精品女同一区二区软件| 如何舔出高潮| 国产成人精品无人区| 亚洲欧洲精品一区二区精品久久久 | 国产熟女欧美一区二区| 伦理电影免费视频| 亚洲精品一区蜜桃| 男人舔奶头视频| 免费观看性生交大片5| 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 日本猛色少妇xxxxx猛交久久| 国产精品一二三区在线看| www.av在线官网国产| 久久久久久伊人网av| 亚洲国产色片| 国产精品.久久久| 各种免费的搞黄视频| 精品一区二区三区视频在线| 成年人免费黄色播放视频 | 嘟嘟电影网在线观看| 国产黄片美女视频| 国产视频首页在线观看| 精品人妻一区二区三区麻豆| 日本与韩国留学比较| 久久韩国三级中文字幕| av视频免费观看在线观看| 亚洲精品日韩av片在线观看| 日韩免费高清中文字幕av| a级一级毛片免费在线观看| 久久99蜜桃精品久久| 91精品伊人久久大香线蕉| 一级爰片在线观看| 一区在线观看完整版| 男男h啪啪无遮挡| 日韩视频在线欧美| 久久精品国产自在天天线| 一级毛片我不卡| 日韩精品免费视频一区二区三区 | 人妻制服诱惑在线中文字幕| 中文精品一卡2卡3卡4更新| 韩国高清视频一区二区三区| av天堂中文字幕网| 久久久久人妻精品一区果冻| 亚洲国产毛片av蜜桃av| 少妇被粗大猛烈的视频| 国产精品99久久99久久久不卡 | 人人妻人人澡人人看| 热99国产精品久久久久久7| 蜜桃在线观看..| 少妇人妻 视频| 中文乱码字字幕精品一区二区三区| 久久av网站| 日本wwww免费看| 麻豆成人午夜福利视频| 亚洲成色77777| 男女边吃奶边做爰视频| 国产乱人偷精品视频| 91午夜精品亚洲一区二区三区| av在线播放精品| 一级爰片在线观看| 哪个播放器可以免费观看大片| 中文字幕精品免费在线观看视频 | 国产欧美日韩精品一区二区| 日韩中字成人| 大香蕉97超碰在线| 亚洲国产日韩一区二区| 国产亚洲午夜精品一区二区久久| 日日啪夜夜爽| 黄片无遮挡物在线观看| 国产片特级美女逼逼视频| 亚洲在久久综合| 中国三级夫妇交换| 18禁动态无遮挡网站| 两个人的视频大全免费| 能在线免费看毛片的网站| 午夜免费男女啪啪视频观看| 爱豆传媒免费全集在线观看| 99热这里只有精品一区| 免费人妻精品一区二区三区视频| 亚洲精品日本国产第一区| 国产在线视频一区二区| 成人影院久久| 欧美97在线视频| av免费观看日本| 成年人午夜在线观看视频| 久久精品国产亚洲av涩爱| 国产精品国产三级国产av玫瑰| 日本黄色片子视频| 久久免费观看电影| 日本欧美国产在线视频| 亚洲成色77777| 精品一区在线观看国产| 日韩视频在线欧美| 亚洲av.av天堂| 色视频www国产| 日本免费在线观看一区| 热re99久久国产66热| 两个人免费观看高清视频 | 黄色欧美视频在线观看| 国产免费一级a男人的天堂| 大码成人一级视频| 在线观看美女被高潮喷水网站| 国产精品久久久久成人av| 国产精品国产三级国产专区5o| 成人亚洲欧美一区二区av| 久久影院123| 又黄又爽又刺激的免费视频.| 久久99蜜桃精品久久| 美女cb高潮喷水在线观看| 免费看光身美女| 亚洲国产欧美在线一区| 久久精品国产鲁丝片午夜精品| 国产极品天堂在线| 久久国内精品自在自线图片| 国产伦在线观看视频一区| av又黄又爽大尺度在线免费看| 美女视频免费永久观看网站| 最新的欧美精品一区二区| 夜夜骑夜夜射夜夜干| 十八禁网站网址无遮挡 | 十八禁网站网址无遮挡 | 九草在线视频观看| 一级二级三级毛片免费看| 少妇精品久久久久久久| 搡女人真爽免费视频火全软件| 高清欧美精品videossex| 精品人妻熟女毛片av久久网站| xxx大片免费视频| 国产欧美日韩一区二区三区在线 | 搡老乐熟女国产| 久久久久久久亚洲中文字幕| 精品一区二区免费观看| 亚洲人成网站在线播| 婷婷色麻豆天堂久久| 久久国内精品自在自线图片| 最近中文字幕2019免费版| 日本av手机在线免费观看| 久热这里只有精品99| 99久久精品热视频| 三上悠亚av全集在线观看 | 丝瓜视频免费看黄片| 欧美亚洲 丝袜 人妻 在线| 中文字幕精品免费在线观看视频 | 日韩成人伦理影院| 欧美国产精品一级二级三级 | 我要看日韩黄色一级片| 欧美变态另类bdsm刘玥| 2018国产大陆天天弄谢| 十八禁高潮呻吟视频 | 国产一区有黄有色的免费视频| 免费观看a级毛片全部| 国产美女午夜福利| 22中文网久久字幕| 国产老妇伦熟女老妇高清| av天堂中文字幕网| 嫩草影院入口| 国产白丝娇喘喷水9色精品| 久久午夜综合久久蜜桃| 欧美变态另类bdsm刘玥| 亚洲欧洲日产国产| 日韩精品免费视频一区二区三区 | 永久网站在线| 高清欧美精品videossex| 女性生殖器流出的白浆| 国产av国产精品国产| 欧美老熟妇乱子伦牲交| 国产视频首页在线观看| 人妻 亚洲 视频| 色吧在线观看| 亚洲精品一区蜜桃| 精品久久久久久电影网| 黑人巨大精品欧美一区二区蜜桃 | freevideosex欧美| 狂野欧美白嫩少妇大欣赏| 午夜福利在线观看免费完整高清在| 午夜福利视频精品| 日韩熟女老妇一区二区性免费视频| 久久午夜综合久久蜜桃| 亚洲精华国产精华液的使用体验| 久久久精品94久久精品| h视频一区二区三区| 各种免费的搞黄视频| 国产在线一区二区三区精| 五月伊人婷婷丁香| 中文字幕制服av| 欧美日韩在线观看h| 狂野欧美激情性bbbbbb| a级毛片在线看网站| 国产亚洲精品久久久com| 老司机影院毛片| 国产av国产精品国产| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久亚洲中文字幕| 黑人高潮一二区| 97在线人人人人妻| 亚洲av男天堂| videos熟女内射| 搡女人真爽免费视频火全软件| 国产精品伦人一区二区| 免费观看a级毛片全部| 另类精品久久| 一区在线观看完整版| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| 99热6这里只有精品| 国产成人精品久久久久久| 99九九在线精品视频 | 曰老女人黄片| 日韩人妻高清精品专区| 熟女人妻精品中文字幕| 国产女主播在线喷水免费视频网站| 大香蕉97超碰在线| 国产免费一级a男人的天堂| 2021少妇久久久久久久久久久| 99视频精品全部免费 在线| 91精品伊人久久大香线蕉| av有码第一页| 黄色毛片三级朝国网站 | 久久国内精品自在自线图片| 制服丝袜香蕉在线| 99久国产av精品国产电影| 你懂的网址亚洲精品在线观看| 久久午夜福利片| 国产白丝娇喘喷水9色精品| 免费观看无遮挡的男女| 国产白丝娇喘喷水9色精品| 啦啦啦视频在线资源免费观看| 亚洲欧美日韩东京热| 亚洲精品久久午夜乱码| 九草在线视频观看| 18禁在线无遮挡免费观看视频| av有码第一页| 久久午夜福利片| av福利片在线观看| 插阴视频在线观看视频| 午夜福利在线观看免费完整高清在| 夜夜骑夜夜射夜夜干| 夜夜看夜夜爽夜夜摸| 最近手机中文字幕大全| 在线看a的网站| 女性被躁到高潮视频| 欧美xxⅹ黑人| 简卡轻食公司| 国产黄色免费在线视频| 性色av一级| 国产午夜精品久久久久久一区二区三区| 99热6这里只有精品| 久久久午夜欧美精品| 97精品久久久久久久久久精品| 日韩av在线免费看完整版不卡| 狂野欧美激情性bbbbbb| 国产精品福利在线免费观看| 另类精品久久| 久久精品久久久久久噜噜老黄| 99热这里只有精品一区| 女人精品久久久久毛片| 99热这里只有是精品50| 夫妻午夜视频| 亚洲性久久影院| 91久久精品国产一区二区三区| 国产日韩欧美亚洲二区| 深夜a级毛片| 国产精品一区www在线观看| 你懂的网址亚洲精品在线观看| 日本av手机在线免费观看| 日韩一区二区视频免费看| 国产成人91sexporn| 欧美 日韩 精品 国产| 国产精品国产三级专区第一集| 最近手机中文字幕大全| 麻豆成人av视频| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲欧美精品专区久久| 国产精品蜜桃在线观看| 精品国产乱码久久久久久小说| 亚洲av不卡在线观看| 免费播放大片免费观看视频在线观看| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 男女无遮挡免费网站观看| 国产国拍精品亚洲av在线观看| 99久久精品一区二区三区| 日韩强制内射视频| 日本黄色日本黄色录像| 午夜91福利影院| 亚洲第一区二区三区不卡| 在线 av 中文字幕| 蜜臀久久99精品久久宅男| 成人二区视频| 日韩不卡一区二区三区视频在线| 日本91视频免费播放| 777米奇影视久久| a级一级毛片免费在线观看| 国产精品一区二区在线不卡| 国产欧美亚洲国产| 午夜福利影视在线免费观看| 欧美精品人与动牲交sv欧美| 国产视频内射| 国产极品粉嫩免费观看在线 | 久久久欧美国产精品| 丰满迷人的少妇在线观看| 欧美日本中文国产一区发布| 亚洲国产最新在线播放| 久久99一区二区三区| 卡戴珊不雅视频在线播放| 久久国产精品大桥未久av | 国产精品偷伦视频观看了| 如何舔出高潮| 十分钟在线观看高清视频www | 久久久久久久久大av| 亚洲第一av免费看| 中国国产av一级| 免费大片黄手机在线观看| 成人无遮挡网站| 国产免费福利视频在线观看| 国产欧美日韩精品一区二区| 国产亚洲午夜精品一区二区久久| 久久婷婷青草| 视频中文字幕在线观看| 久久国产精品大桥未久av | 亚洲欧美精品专区久久| 黑人巨大精品欧美一区二区蜜桃 | 国内揄拍国产精品人妻在线| 国产一区有黄有色的免费视频| 午夜久久久在线观看| 国产毛片在线视频| 好男人视频免费观看在线| 男女免费视频国产| 日韩一区二区三区影片| 日韩不卡一区二区三区视频在线| 黄色欧美视频在线观看| 欧美性感艳星| 久久久欧美国产精品| 国产亚洲最大av| 少妇的逼水好多| a级毛片免费高清观看在线播放| 久久这里有精品视频免费| 中文精品一卡2卡3卡4更新| 欧美精品高潮呻吟av久久| 亚洲欧美成人综合另类久久久| 嘟嘟电影网在线观看| 日本av免费视频播放| 久久精品国产亚洲av天美| 国产男人的电影天堂91| 国产伦在线观看视频一区| 大片免费播放器 马上看| 在线看a的网站| 青春草国产在线视频| 又大又黄又爽视频免费| 美女中出高潮动态图| 久久人人爽av亚洲精品天堂| 欧美区成人在线视频| a 毛片基地| 亚洲天堂av无毛| 日本91视频免费播放| 一级毛片电影观看| 看免费成人av毛片| 男人和女人高潮做爰伦理| 只有这里有精品99| 精品亚洲乱码少妇综合久久| 一级毛片电影观看| 日韩熟女老妇一区二区性免费视频| 99久国产av精品国产电影| 国产亚洲5aaaaa淫片| 老熟女久久久| 国产亚洲午夜精品一区二区久久| 国产欧美亚洲国产| 成年人午夜在线观看视频| 91aial.com中文字幕在线观看| 高清在线视频一区二区三区| 一级毛片电影观看| 三级经典国产精品| 国产成人aa在线观看| 3wmmmm亚洲av在线观看| 欧美97在线视频| 少妇裸体淫交视频免费看高清| 欧美日韩综合久久久久久| 丝袜脚勾引网站| 一区在线观看完整版| 欧美精品高潮呻吟av久久| 99九九线精品视频在线观看视频| 精品国产一区二区三区久久久樱花| 欧美xxxx性猛交bbbb| 成人国产麻豆网| h视频一区二区三区| 免费观看性生交大片5| 久久人妻熟女aⅴ| 国产精品一二三区在线看| 国产成人91sexporn| 少妇的逼水好多| 毛片一级片免费看久久久久| 精品人妻偷拍中文字幕| 在线亚洲精品国产二区图片欧美 | 国产 精品1| 国产熟女午夜一区二区三区 | 纵有疾风起免费观看全集完整版| 国产一级毛片在线| 国产亚洲91精品色在线| 亚洲国产日韩一区二区| 中文乱码字字幕精品一区二区三区| 亚洲天堂av无毛| 免费少妇av软件| 精品少妇久久久久久888优播| 日本-黄色视频高清免费观看| 啦啦啦啦在线视频资源| 亚洲自偷自拍三级| 国产成人freesex在线| 2021少妇久久久久久久久久久| 欧美日韩在线观看h| 亚洲av成人精品一区久久| 一级毛片久久久久久久久女| 精品熟女少妇av免费看| 丝袜在线中文字幕| 国产男女内射视频| 久久久国产精品麻豆| 少妇 在线观看| 国产伦精品一区二区三区视频9| 欧美变态另类bdsm刘玥| 午夜91福利影院| 久久99精品国语久久久| 99热这里只有是精品50| 亚洲色图综合在线观看| 国产永久视频网站| 欧美成人午夜免费资源| 成人影院久久| 夫妻性生交免费视频一级片| 高清欧美精品videossex| 男的添女的下面高潮视频| 狂野欧美白嫩少妇大欣赏| av播播在线观看一区| 自线自在国产av| 日韩大片免费观看网站| 青春草视频在线免费观看| 国产精品人妻久久久影院| 夜夜爽夜夜爽视频| √禁漫天堂资源中文www| 午夜福利,免费看| 国内精品宾馆在线| 精华霜和精华液先用哪个| h视频一区二区三区| 99久久综合免费| 欧美亚洲 丝袜 人妻 在线| av在线老鸭窝| 亚洲色图综合在线观看| a级毛片在线看网站| 97在线人人人人妻| 国产综合精华液| 国产精品一区二区性色av| 三上悠亚av全集在线观看 | 亚洲中文av在线| 国产中年淑女户外野战色| 成人国产av品久久久| 久久99一区二区三区| 狠狠精品人妻久久久久久综合| 精品久久国产蜜桃| 国产av精品麻豆| 欧美丝袜亚洲另类| 一边亲一边摸免费视频| h日本视频在线播放| 麻豆精品久久久久久蜜桃| 五月玫瑰六月丁香| 女人久久www免费人成看片| 97精品久久久久久久久久精品| 我的女老师完整版在线观看| 久久精品国产亚洲av涩爱| 精品人妻熟女av久视频| 国产一区二区在线观看av| 人人妻人人爽人人添夜夜欢视频 | 一级毛片黄色毛片免费观看视频| 国产一级毛片在线| 欧美另类一区| 日本免费在线观看一区| 在线亚洲精品国产二区图片欧美 | 热re99久久精品国产66热6| 午夜激情久久久久久久| 久久女婷五月综合色啪小说| 妹子高潮喷水视频| 免费久久久久久久精品成人欧美视频 | 春色校园在线视频观看| 欧美日韩精品成人综合77777| 少妇被粗大的猛进出69影院 | 久久久国产一区二区| 国产男女超爽视频在线观看| 男人狂女人下面高潮的视频| 色哟哟·www| 夜夜看夜夜爽夜夜摸| 午夜精品国产一区二区电影| 中国国产av一级|