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

    水下隱身動目標(biāo)激光探測水面能量散射特性研究*

    2016-08-10 03:32:49孫曉峰陶永勤
    艦船電子工程 2016年7期
    關(guān)鍵詞:海面潛艇激光

    姚 瑤 孫曉峰 陶永勤

    (海軍702廠 上海 200434)

    ?

    水下隱身動目標(biāo)激光探測水面能量散射特性研究*

    姚瑤孫曉峰陶永勤

    (海軍702廠上海200434)

    摘要為了提高衛(wèi)星對水下隱身動目標(biāo)的探測效率,對不同海況下、不同入射角度時的三維激光海面漫反射特性進(jìn)行了建模仿真計算,并在實驗室條件下對激光入射波動水面后產(chǎn)生的散射場的能量分布特性進(jìn)行了的研究。首先,根據(jù)麥克斯韋爾方程和邊界條件對散射場在各個方向上的分量之間的耦合關(guān)系進(jìn)行了研究,并列出矩陣方程;然后,利用三維錐形波對入射激光束進(jìn)行了模擬;最后,通過稀疏矩陣規(guī)則網(wǎng)格法求解矩陣方程,得出三維雙站散射系數(shù),并通過實驗對仿真結(jié)果加以驗證。結(jié)果表明該方法能夠準(zhǔn)確地表示波動水面的三維激光漫反射特性,為進(jìn)一步研究三維激光海面漫反射特性奠定了基礎(chǔ)。

    關(guān)鍵詞激光; 潛艇; 無人水下自航器; 漫反射; 雙站散射系數(shù); 海面

    Class NumberTN219

    1引言

    潛艇、UUV等水下隱身動目標(biāo)具有極大的戰(zhàn)略威懾力和極強(qiáng)的戰(zhàn)場信息偵察能力,使得傳統(tǒng)的聲吶探測手段遇到了新的挑戰(zhàn)。因此,利用綠激光對海水的穿透特性以及綠激光對水下隱身動目標(biāo)特征的Mie散射效應(yīng),可以作為探潛體系的一種補(bǔ)充手段,以進(jìn)一步提高對水下隱身動目標(biāo)的探測效能。由于不同海況、不同角度入射的激光散射能量分布具有很大的差異,對于探測效率有著非常大的影響。因此,針對不同海況,得出準(zhǔn)確的三維激光海面漫反射特性,對于提高衛(wèi)星對水下隱身動目標(biāo)的探測范圍和準(zhǔn)確率具有重要的意義。

    由于入射激光束經(jīng)過粗糙界面漫反射后,不同偏振狀態(tài)的激光在各個方向上有著復(fù)雜的耦合關(guān)系。因此為了方便計算,以往求解方法大多基于一定的近似條件,如基爾霍夫近似(Kirchhoff Approach, KA)方法[1],這些近似條件不但對粗糙面的參數(shù)有一定的限制,由于其固有的物理近似勢必帶來誤差,使得計算結(jié)果達(dá)不到理想的效果。二維前后向迭代法(Forward-Backward Method, FBM)[2]雖然是數(shù)值計算方法,但是它是將三維的粗糙界面簡化為二維界面后得出的數(shù)值模擬。由于其所建立的二維粗糙界面與現(xiàn)實中的粗糙界面差別很大,因此仿真計算結(jié)果的準(zhǔn)確性較差。

    因此,為了得到準(zhǔn)確的激光在波動水面上的漫反射能量分布特性,本文采用稀疏矩陣規(guī)則網(wǎng)格法(Sparse Matrix Canonical Grid Method, SMCG)來計算激光在波動水面上的雙站散射系數(shù)(Bistatic Scattering Coefficient, BSC),該算法不僅減小了計算機(jī)內(nèi)存的需求,而且由于在每步迭代過程中矩陣與向量的運算又采用快速傅里葉變換(Fast Fourier Transition, FFT)算法,進(jìn)一步減小了CPU的計算時間。

    2激光漫反射原理

    假設(shè)一束激光入射到波動水面S上,則三維激光漫反射的剖面示意圖如圖1所示。

    圖1 激光在波動水面漫反射的幾何關(guān)系

    (1)

    那么,根據(jù)Stratton-Chu公式[3],結(jié)合由麥克斯韋爾方程推導(dǎo)得到散射場的電場和磁場的積分方程和邊界條件[4],通過對粗糙海面散射場切向和法向進(jìn)行分析,可以得出以下方程:

    (2a)

    (2b)

    (2c)

    (2d)

    (3)

    (4)

    (5)

    3矩陣離散

    將光斑在x-y平面上的投影區(qū)域Sxy均勻劃分成Nx×Ny個網(wǎng)格,該區(qū)域的每個采樣點位于每個網(wǎng)格的中心處。在x和y方向上,每個網(wǎng)格的邊長分別為Δlx=Lx/Nx和Δly=Ly/Ny。令m,n分別代表該點在網(wǎng)格中的橫向和豎向序號,則該區(qū)域上的積分順序如圖2所示。

    從圖2中可以看出,采樣點的排序是按照從下到上,再從左到右的順序。其中,rd的定義將在下節(jié)進(jìn)行介紹。

    (6a)

    (6b)

    (6c)

    圖2 Sxy的離散示意圖

    (7a)

    (7b)

    (7c)

    經(jīng)過矢量分解后,根據(jù)式(2)可以得到以下六個關(guān)于標(biāo)量場的方程。為了便于數(shù)值計算,根據(jù)矩量法(Method of Moment, MOM)對這六個標(biāo)量方程進(jìn)行離散,離散后的矩陣方程表示為

    (8)

    4稀疏矩陣規(guī)則網(wǎng)格法(SMCG)

    由于源點和接收點的水平距離越近,源點對接收點的貢獻(xiàn)就越大,因此,離接收點近的區(qū)域為強(qiáng)區(qū),反之為弱區(qū)。SMCG[6]是一種效率很高的數(shù)值計算方法,它對強(qiáng)區(qū)源點的貢獻(xiàn)進(jìn)行詳細(xì)計算,而對弱區(qū)源點的貢獻(xiàn)則采取近似計算,在保證準(zhǔn)確度的同時,極大地節(jié)約了計算成本。

    為了便于說明,假設(shè)存在如下方程:

    (9)

    按照上述的強(qiáng)弱區(qū)域劃分原則,可以將式(9)中的系數(shù)矩陣分解為如下形式:

    Z=Z(s)+Z(w)

    (10)

    用稀疏矩陣來表示強(qiáng)區(qū)矩陣Z(s),對于弱區(qū)矩陣Z(w)的元素,將其中的格林函數(shù)用泰勒級數(shù)展開求和的方式來近似表示,即:

    (11)

    (12)

    其中,am,bm分別代表利用泰勒級數(shù)對G(R)和g(R)進(jìn)行展開時的第m級系數(shù),下標(biāo)或上標(biāo)中的0代表大氣中的相應(yīng)參數(shù),1代表水中的相應(yīng)參數(shù)。

    在數(shù)值計算中,為了達(dá)到理想的精度,將泰勒級數(shù)展開至第6項,即M=5。根據(jù)泰勒級數(shù)的展開式,可以將弱區(qū)矩陣Z(w)表示為級數(shù)求和的形式,即:

    (13)

    令式(13)中的第1項為x-y平面的貢獻(xiàn)Z(FS),即:

    (14)

    則SMCG的迭代過程可以簡單表述如下:

    (15a)

    (15b)

    (15c)

    第n次迭代過程結(jié)束后的相對誤差Er定義為

    (16)

    在具體的求解過程中,采用共軛梯度法(Conjugate Gradient Method, CGM)[7]來進(jìn)行計算。由于Z(FS)是一個滿矩陣,因此CGM在每步的迭代過程中,要消耗巨大的計算量來計算它與向量的乘積。通過Z(FS)的表達(dá)式發(fā)現(xiàn),ZFS是一個分塊托普利茲矩陣。它與向量的乘積滿足卷積運算,因而可以用FFT來實現(xiàn)。令T為Z(FS)的一個分塊矩陣,T的第1列為T1,則具體的計算方法為:

    Tx=T1?x=ifft(fft(T1)·fft(x))

    (17)

    其中,?表示卷積運算符,fft表示快速傅里葉變換,ifft表示快速傅里葉反變換。

    5錐形入射波

    (18)

    (19)

    (20)

    (21)

    (22)

    t=tx+ty

    (23)

    (24)

    (25)

    (26)

    入射激光束的總能量為

    (27)

    6數(shù)值計算與分析

    假設(shè)波動水面是一種分布式的隨機(jī)過程,則采用歸一化微分散射截面,即散射系數(shù)來表征粗糙面角度性的散射特征。三維激光漫反射在不同偏振狀態(tài)下的雙站散射系數(shù)可以被定義為[9]

    (28)

    通過對雙站散射系數(shù)進(jìn)行立體角的積分求和,可以得到不同偏振狀態(tài)的入射激光以角度θi入射時所產(chǎn)生的散射場的總能量與入射能量的比值Prβ(θi),具體的表達(dá)式為

    (29)

    本文采用基于高斯譜的高斯隨機(jī)過程[10]來模擬二維隨機(jī)粗糙海面。為了模擬粗糙海面不同的粗糙度,建立了A、B兩種海面模型來模擬不同粗糙度的波動水面。其中模型A的均方高度hx=hy=0.5cm,相關(guān)長度lx=ly=6cm;模型B的均方高度hx=hy=1cm,相關(guān)長度lx=ly=5cm。

    為了說明三維激光海面漫反射的能量分布特性,以水平偏振為例,仿真計算了三維激光海面漫反射的雙站散射系數(shù),即γ=γhh。在仿真計算實例中,分別取入射角θi1=20°和θi2=60°,入射方位角φi1=φi2=0°。由于激光的能量相當(dāng)集中,為了保證仿真結(jié)果的準(zhǔn)確性,將計算了光斑的全部面積,參與散射的粗糙面的面積取為光斑面積的16倍,則相應(yīng)的錐形波照射區(qū)域參數(shù)和粗糙水面的參數(shù)分別為gx1=gy1=0.01cm,Lx1=Ly1=0.04cm和gx2=0.02cm,gy2=0.04cm,Lx2=0.08cm,Ly2=0.16cm。激光波長為λ=0.532μm,為了保證激光的電磁特性,在x和y方向上的采樣間隔為Δlx=Δly=λ/10。

    采用SMCG進(jìn)行數(shù)值計算,當(dāng)取入射角θi1=20°時,用來區(qū)分強(qiáng)弱區(qū)域的臨界距離rd1=40λ,當(dāng)取入射角θi2=60°時,用來區(qū)分強(qiáng)弱區(qū)域的臨界距離rd2=100λ。為了得到散射場能量分布的統(tǒng)計特性,采用蒙特卡洛方法進(jìn)行仿真,蒙特卡洛仿真次數(shù)為500次。具體的三維激光海面漫反射雙站散射系數(shù)的仿真結(jié)果如圖3和圖4所示。

    圖3 當(dāng)入射角θi1=20°入射不同粗糙度的界面時的雙站散射系數(shù)

    由于Matlab中沒有直接繪制屬于三維極坐標(biāo)系下數(shù)據(jù)圖形的命令,如果直接采用散射角θs和散射方位角φs作為直角坐標(biāo)系的變量的話,當(dāng)散射方位角θs相同時,φs=0°和φs=360°的雙站散射系數(shù)是無法重合的。因此,為了形象、直觀地描述散射場三維能量分布的角度特性,將散射角θs和散射方位角φs轉(zhuǎn)換為直角坐標(biāo)系下的坐標(biāo)值Px=sinθscosφs,Py=sinθssinφs。

    從圖3和圖4中不難看出,當(dāng)界面的粗糙度相對較小時,散射場的能量分布相對集中,鏡面反射方向的峰值能量相對較大,鄰近區(qū)域的能量相對較小;當(dāng)界面的粗糙度適度增加時,鏡面反射方向的峰值能量減少,鄰近區(qū)域的能量有所增長,散射場的能量分布更加均勻,即散射場能量的覆蓋范圍變得更大。但是在散射場的邊緣處,能量衰減速度仍然非??臁?/p>

    圖4 當(dāng)入射角θi2=60°入射不同粗糙度的界面時的雙站散射系數(shù)

    隨著入射角度的增大,在粗糙面上的光斑面積和粗糙面的截斷區(qū)域面積都會隨之增加,因此就同一粗糙面而言,入射角度不同,粗糙面對激光的漫反射程度也不同,散射場的能量分布也存在這較大差別。通過圖3(b)和圖4(b)的對比可以看出,當(dāng)入射角度θi1=20°時,散射場的能量分布相對于入射角度θi2=60°時的散射場的能量分布而言,峰值能量相對較弱,總體能量分布變得更加均勻。這就說明對于相同粗糙度的粗糙面模型,入射角度越小,粗糙面對激光的漫反射程度就會越大,散射場的能量分布就會越均勻。

    7實驗驗證與分析

    為了對三維激光海面漫反射雙站散射系數(shù)的數(shù)值計算結(jié)果進(jìn)行驗證,在實驗室中進(jìn)行了實驗。在實驗室中,分別取入射角θi1=20°和θi2=60°,入射方位角φi1=φi2=0°,入射高度H=0.1m,采用功率為P=100mW,激光光束的發(fā)散角a=1mrad,波長λ=0.532μm的綠光水平偏振連續(xù)激光器。為了獲得較為準(zhǔn)確的實驗數(shù)據(jù),利用造波器生成兩種不同粗糙度的波動水面,其中波動水面A的均方高度hx=hy=0.5cm,相關(guān)長度lx=ly=6cm;波動水面B的均方高度hx=hy=1cm,相關(guān)長度lx=ly=5cm。利用147型激光功率計對散射場中不同方位的采樣點的能量進(jìn)行多次測量,測量結(jié)果經(jīng)過處理后,再與仿真計算結(jié)果進(jìn)行對比。

    實驗中,對激光在入射角度θi2=60°,入射高度H=0.1m的情況下入射在不同粗糙度的波動水面上時,將散射場的能量變化呈現(xiàn)在另一端的背景墻上,并對散射場的漫反射能量分布情況進(jìn)行了拍攝,具體的拍攝結(jié)果如圖5所示。

    圖5 不同粗糙度的水面產(chǎn)生的散射場

    通過圖5(a)與圖5(b)的對比,可以觀察到圖5(b)中背景墻上的光斑面積明顯增大,整個光斑的亮度分布仍然前向鏡面反射方向附近的亮度最高,越靠近邊緣,亮度越低。但是,相對于圖5(a)而言,光斑中心的亮度明顯降低,而且在邊緣處的能量衰減速度也相對較慢。這說明粗糙度較大的水面使得散射場的能量分布更加均勻,散射場的覆蓋面積明顯增大。

    以上通過直接觀察得到的激光漫反射能量分布特性,與對仿真結(jié)果進(jìn)行分析后得出的激光漫反射能量分布特性是一致的。接下來,為了獲得更加客觀、量化的激光漫反射能量分布特性,將利用通過多次實驗測得的數(shù)據(jù),對三種不同方法(KA、2DFBM和3DSMCG)得出的散射場能量分布特性進(jìn)行驗證、對比和分析。

    為了更加清楚地表示散射場能量在前后向上的分布,令散射方位角φs=0°時的散射角為正,即為前向。φs=180°時的散射角為負(fù),即為后向,則采用的三種不同方法所得出的漫反射率以及實驗數(shù)據(jù)之間的對比結(jié)果如圖6、圖7所示。

    圖6 θi=20°時的漫反射率

    圖7 θi=60°時的漫反射率

    通過圖6和圖7中二維和三維仿真結(jié)果與實驗數(shù)據(jù)的比對,發(fā)現(xiàn)通過三維數(shù)值計算仿真得出的漫反射率與實驗數(shù)據(jù)更為接近。由于KA是一種解析算法,只能大致地估算散射場的能量分布特性,因此它是三種方法中最不準(zhǔn)確的一個,甚至在圖6(b)中,并未體現(xiàn)出散射場的后向增強(qiáng)現(xiàn)象。而二維數(shù)值計算的仿真結(jié)果相對于實驗數(shù)據(jù)要大一些。從理論上分析,這主要是因為二維數(shù)值計算的仿真結(jié)果所代表的能量之和就是粗糙面上所有源點產(chǎn)生的貢獻(xiàn)之和,即二維散射場的總能量。而三維數(shù)值計算的仿真結(jié)果與真實情況更為接近,圖5和圖6中代表三維漫反射率的曲線所代表的能量之和只是粗糙面上所有源點產(chǎn)生的貢獻(xiàn)的一部分,并不是三維散射場的總能量。雖然三種計算方法在準(zhǔn)確度上存在差異,但是所反映出來的散射場的能量分布特性上來講還是相近的。

    關(guān)于實驗數(shù)據(jù)與三維雙站散射系數(shù)之間存在的偏差,經(jīng)過認(rèn)真分析之后,認(rèn)為主要是由于以下原因造成的:

    1) 實驗條件的客觀因素,不可避免的會出現(xiàn)一些誤差;

    2) 粗糙海面的理論模型和真實水面還有一定差距。

    8結(jié)語

    本文采用數(shù)值計算方法對激光在波動水面上的三維漫反射特性進(jìn)行了研究和實驗驗證。并且通過與其他方法得出的計算結(jié)果進(jìn)行對比之后,認(rèn)為三維數(shù)值計算方法得出的三維BSC是最能準(zhǔn)確反映散射場能量分布特性的。通過分析研究,發(fā)現(xiàn)散射場的能量分布與入射角度和水面的粗糙度有著密切的關(guān)系,區(qū)分散射場能量前后向分布大小的入射角度隨著粗糙度的變化而變化,粗糙度越大,入射角的分界點就越大。并且較為粗糙的水面會使得散射場的能量分布更加均勻,在散射場的邊緣

    處能量迅速衰減。

    參 考 文 獻(xiàn)

    [1] 劉春香,程傳福,任曉榮,等.隨機(jī)表面散射光場的格林函數(shù)法與基爾霍夫近似的比較[J].物理學(xué)報,2004,53(2):427-435.

    [2] 梁華鋒,賴建軍,彭中良,等.用于氣體傳感的一維光子晶體窄帶紅外發(fā)射光源[J].紅外與毫米波學(xué)報,2009,28(6):414-417.

    [3] 楊雨川,譚吉春,盛定儀,等.導(dǎo)體板上單極天線對電磁脈沖響應(yīng)特性的矩量法分析[J].強(qiáng)激光與粒子束,2008,20(4):649-652.

    [4] 崔志偉,韓平,李明磊.磁化等離子覆蓋導(dǎo)體散射問題的FE/BI方法分析[J].電波科學(xué)學(xué)報,2010(5):973-978.

    [5] 吳國榮,鐘偉芳,吳永東.半空間彈性波強(qiáng)散射雙參數(shù)非線性反演[J].華中科技大學(xué)(自然科學(xué)版),2003,31(4):105-107.

    [6] H T Chou and J T Johnson. A novel acceleration algorithm for the computation of scattering from rough surfaces with the forward-backward method[J]. Radio Sci.,1998,33(4):1277-1287.

    [7] Tsang L, Kong J A. Scattering of Electromagnetic Waves-Numerical Simulation[M]. New York: John Wiley,2001:201-203.

    [8] 葉紅霞,金亞秋.隨機(jī)粗糙面上介質(zhì)目標(biāo)差場散射的快速計算方法[J].電子學(xué)報,2006,34(6):1074-1079.

    [9] Tsang L, Kong J A. Scattering of Electromagnetic Waves-Numerical Simulation[M]. New York: John Wiley,2001:56-60.

    [10] Lei Kuang, Yaqiu Jin. Bistatic Scattering From a Three-Dimensional Object Over a Randomly Rough Surface Using the FDTD Algorithm[J]. Antennas and Propagation, IEEE Transactions on,2007,55(8):2302-2312.

    收稿日期:2016年1月14日,修回日期:2016年3月2日

    作者簡介:姚瑤,女,博士,工程師,研究方向:水下目標(biāo)探測。

    中圖分類號TN219

    DOI:10.3969/j.issn.1672-9730.2016.07.037

    Characteristics of Laser Scattering from Rough Water Surface Based on Detection of Underwater Stealth Moving Target

    YAO YaoSUN XiaofengTAO Yongqin

    (No.702 Factory of Navy, Shanghai200434)

    AbstractTo enhance the detection efficiency of satellite against underwater stealth moving targets, the characteristics of the scattered field of laser from the water surface in different roughness and incident angles was studied. First, the coupling relationships among components in different directions were analyzed and the corresponding matrix equations from Maxwell equations and boundary conditions were solved with SMCG. Second, the 3D model of incident laser beam was established with 3D tapered wave. Finally, the 3D BSC was get and verified by experiments. The results showed that the 3D BSC was enough accurate to express the characteristics of laser diffusion from water surface so that some contributions were made for the further study about the characteristics of laser diffusion from sea surface.

    Key Wordslaser, submarine, UUV, diffusion, bistatic scattering coefficient, sea surface

    猜你喜歡
    海面潛艇激光
    十分鐘讀懂潛艇史(下)
    潛艇哥別撞我
    Er:YAG激光聯(lián)合Nd:YAG激光在口腔臨床醫(yī)學(xué)中的應(yīng)用
    海面床,輕輕搖
    十分鐘讀懂潛艇史(上)
    潛艇躍進(jìn)之黃金時代
    激光誕生60周年
    科學(xué)(2020年5期)2020-11-26 08:19:24
    第六章 邂逅“胖胖號”
    請給激光點個贊
    国产97色在线日韩免费| 亚洲人成电影观看| 国产精品99久久99久久久不卡| 日韩大尺度精品在线看网址 | 亚洲国产精品久久男人天堂| 电影成人av| 免费女性裸体啪啪无遮挡网站| 亚洲中文av在线| 女生性感内裤真人,穿戴方法视频| 久久婷婷人人爽人人干人人爱 | 国产高清有码在线观看视频 | 久久国产亚洲av麻豆专区| 国产单亲对白刺激| 热99re8久久精品国产| 叶爱在线成人免费视频播放| 国产激情欧美一区二区| 国内久久婷婷六月综合欲色啪| 亚洲国产中文字幕在线视频| 免费观看精品视频网站| 搡老熟女国产l中国老女人| 精品国产乱子伦一区二区三区| 亚洲国产精品久久男人天堂| 亚洲av片天天在线观看| 亚洲专区中文字幕在线| 天天添夜夜摸| 国产亚洲欧美98| 成人精品一区二区免费| 免费看a级黄色片| 亚洲精品国产精品久久久不卡| 国产亚洲av嫩草精品影院| 在线播放国产精品三级| 一区二区三区激情视频| 在线免费观看的www视频| 国产精品 国内视频| 国产亚洲欧美精品永久| 欧美精品啪啪一区二区三区| 国产高清激情床上av| 免费搜索国产男女视频| 正在播放国产对白刺激| 亚洲一区二区三区色噜噜| 欧美激情 高清一区二区三区| 最近最新免费中文字幕在线| 黄色丝袜av网址大全| 国产精品一区二区精品视频观看| 久久婷婷人人爽人人干人人爱 | 黄色毛片三级朝国网站| 自线自在国产av| 美女高潮到喷水免费观看| 一区在线观看完整版| 日韩欧美国产在线观看| 精品国产亚洲在线| 久久久国产精品麻豆| 国产精品美女特级片免费视频播放器 | 一二三四在线观看免费中文在| 美女免费视频网站| 免费看美女性在线毛片视频| 成人亚洲精品av一区二区| 久久精品亚洲熟妇少妇任你| 一级a爱片免费观看的视频| 欧美日韩乱码在线| 91麻豆av在线| 麻豆成人av在线观看| 男女做爰动态图高潮gif福利片 | 色婷婷久久久亚洲欧美| 欧美av亚洲av综合av国产av| 色精品久久人妻99蜜桃| 一级毛片女人18水好多| 国产成人欧美| 亚洲成av片中文字幕在线观看| 国产色视频综合| 91老司机精品| 欧美一区二区精品小视频在线| 91精品三级在线观看| 黄片大片在线免费观看| 法律面前人人平等表现在哪些方面| 曰老女人黄片| 久久国产精品男人的天堂亚洲| 搡老妇女老女人老熟妇| 校园春色视频在线观看| 一边摸一边抽搐一进一出视频| 国产精品香港三级国产av潘金莲| 桃红色精品国产亚洲av| 看片在线看免费视频| 97碰自拍视频| 亚洲九九香蕉| 亚洲av电影不卡..在线观看| 久久久久国产精品人妻aⅴ院| 国产真人三级小视频在线观看| 亚洲国产精品成人综合色| 法律面前人人平等表现在哪些方面| 一二三四在线观看免费中文在| 色综合亚洲欧美另类图片| 久久久久九九精品影院| 亚洲国产欧美网| 国产亚洲精品第一综合不卡| 99国产精品一区二区蜜桃av| 日韩欧美在线二视频| 麻豆成人av在线观看| 亚洲成av片中文字幕在线观看| 在线观看免费午夜福利视频| 搡老熟女国产l中国老女人| 国产精品久久久av美女十八| 在线观看午夜福利视频| 可以在线观看毛片的网站| 大型黄色视频在线免费观看| 叶爱在线成人免费视频播放| 国产aⅴ精品一区二区三区波| 妹子高潮喷水视频| 一个人免费在线观看的高清视频| 婷婷丁香在线五月| 成人特级黄色片久久久久久久| 十八禁网站免费在线| 在线观看免费午夜福利视频| 成人国产综合亚洲| 日韩三级视频一区二区三区| www国产在线视频色| 老熟妇乱子伦视频在线观看| 在线天堂中文资源库| 三级毛片av免费| 99久久99久久久精品蜜桃| 国产亚洲精品第一综合不卡| 日韩av在线大香蕉| 中文字幕av电影在线播放| 日韩欧美国产一区二区入口| 大型av网站在线播放| 日本三级黄在线观看| 午夜两性在线视频| videosex国产| 欧洲精品卡2卡3卡4卡5卡区| 成人特级黄色片久久久久久久| 国产一区二区在线av高清观看| e午夜精品久久久久久久| 成人国产综合亚洲| 人人妻人人爽人人添夜夜欢视频| 非洲黑人性xxxx精品又粗又长| 欧美乱色亚洲激情| 免费av毛片视频| 麻豆国产av国片精品| 国产成人精品久久二区二区免费| 在线观看一区二区三区| 欧美日韩福利视频一区二区| 免费无遮挡裸体视频| 91麻豆精品激情在线观看国产| 国产熟女午夜一区二区三区| 亚洲国产欧美网| 人人妻,人人澡人人爽秒播| 国产精品二区激情视频| 亚洲一码二码三码区别大吗| 香蕉国产在线看| 婷婷丁香在线五月| 啦啦啦韩国在线观看视频| 777久久人妻少妇嫩草av网站| ponron亚洲| 国产主播在线观看一区二区| 成人亚洲精品一区在线观看| 国产色视频综合| 欧美国产日韩亚洲一区| 精品卡一卡二卡四卡免费| 国产乱人伦免费视频| 亚洲成人精品中文字幕电影| 亚洲精品国产精品久久久不卡| 大型av网站在线播放| 国产97色在线日韩免费| 在线国产一区二区在线| 动漫黄色视频在线观看| 丰满的人妻完整版| 国产成人精品久久二区二区91| 两人在一起打扑克的视频| 国产精品野战在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲久久久国产精品| 丝袜人妻中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 亚洲一区中文字幕在线| 亚洲国产精品999在线| 亚洲人成77777在线视频| 99精品久久久久人妻精品| 日韩大尺度精品在线看网址 | 91国产中文字幕| 中文字幕色久视频| 精品一区二区三区av网在线观看| 97人妻天天添夜夜摸| 涩涩av久久男人的天堂| 国产一区二区激情短视频| 高清毛片免费观看视频网站| 黄色丝袜av网址大全| 国产精品99久久99久久久不卡| 亚洲精品美女久久av网站| 亚洲第一欧美日韩一区二区三区| 身体一侧抽搐| 亚洲国产精品成人综合色| 一级毛片女人18水好多| 精品国产一区二区三区四区第35| 国产成人av教育| 色老头精品视频在线观看| 久久精品影院6| 成年女人毛片免费观看观看9| 国产精品一区二区精品视频观看| 国产av精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品国产区一区二| 久久人妻熟女aⅴ| 正在播放国产对白刺激| 国产又爽黄色视频| 99久久国产精品久久久| 国产精品1区2区在线观看.| 久久久水蜜桃国产精品网| 精品欧美一区二区三区在线| 999久久久国产精品视频| 高清黄色对白视频在线免费看| 人人澡人人妻人| 日韩欧美一区二区三区在线观看| 久久久国产成人免费| 欧美性长视频在线观看| 色综合站精品国产| 国产精品一区二区免费欧美| 免费少妇av软件| 国产欧美日韩一区二区精品| 搞女人的毛片| 午夜福利免费观看在线| 两个人看的免费小视频| 亚洲精品久久成人aⅴ小说| 中文字幕人妻熟女乱码| 亚洲三区欧美一区| 亚洲aⅴ乱码一区二区在线播放 | 夜夜躁狠狠躁天天躁| 久久国产乱子伦精品免费另类| 丝袜人妻中文字幕| 9热在线视频观看99| 色综合站精品国产| 亚洲欧美日韩高清在线视频| 久久国产精品男人的天堂亚洲| 99在线视频只有这里精品首页| 成人国产综合亚洲| 制服丝袜大香蕉在线| www国产在线视频色| 欧美 亚洲 国产 日韩一| 两个人免费观看高清视频| cao死你这个sao货| 波多野结衣巨乳人妻| 两个人免费观看高清视频| 国产精品秋霞免费鲁丝片| 在线观看66精品国产| 国产一区二区三区在线臀色熟女| 两个人看的免费小视频| 男人舔女人下体高潮全视频| e午夜精品久久久久久久| 亚洲国产高清在线一区二区三 | 一级,二级,三级黄色视频| 麻豆一二三区av精品| 免费少妇av软件| 在线观看日韩欧美| 精品一区二区三区视频在线观看免费| 搡老妇女老女人老熟妇| 精品久久久久久久人妻蜜臀av | 亚洲免费av在线视频| 亚洲欧美激情在线| 老鸭窝网址在线观看| 色综合婷婷激情| 亚洲视频免费观看视频| 校园春色视频在线观看| 亚洲国产精品sss在线观看| 午夜精品久久久久久毛片777| 又黄又爽又免费观看的视频| 丝袜美腿诱惑在线| 成人精品一区二区免费| 十八禁网站免费在线| 亚洲国产中文字幕在线视频| 国产在线观看jvid| 国产熟女午夜一区二区三区| 午夜两性在线视频| 啪啪无遮挡十八禁网站| 久久精品亚洲精品国产色婷小说| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲av成人一区二区三| АⅤ资源中文在线天堂| 正在播放国产对白刺激| 国产精品香港三级国产av潘金莲| 亚洲成a人片在线一区二区| 丝袜在线中文字幕| 国产色视频综合| 激情在线观看视频在线高清| 欧美日韩亚洲综合一区二区三区_| 国产99久久九九免费精品| 欧美黑人精品巨大| 成年女人毛片免费观看观看9| 999久久久国产精品视频| 一本久久中文字幕| 一级a爱视频在线免费观看| 免费观看人在逋| 我的亚洲天堂| 久久久久久人人人人人| av有码第一页| 禁无遮挡网站| 国产日韩一区二区三区精品不卡| 一级a爱片免费观看的视频| 国产极品粉嫩免费观看在线| 欧美一区二区精品小视频在线| 最近最新免费中文字幕在线| 老熟妇乱子伦视频在线观看| 男女之事视频高清在线观看| 国产亚洲精品一区二区www| 国产成年人精品一区二区| 成人亚洲精品av一区二区| 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 亚洲国产欧美日韩在线播放| 给我免费播放毛片高清在线观看| 一边摸一边做爽爽视频免费| 麻豆av在线久日| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久免费高清国产稀缺| 亚洲aⅴ乱码一区二区在线播放 | 欧美乱码精品一区二区三区| 国产精品日韩av在线免费观看 | 欧美午夜高清在线| 亚洲一卡2卡3卡4卡5卡精品中文| 99精品久久久久人妻精品| 熟女少妇亚洲综合色aaa.| 999久久久精品免费观看国产| 午夜福利免费观看在线| 国产人伦9x9x在线观看| 日韩欧美国产一区二区入口| 又大又爽又粗| 亚洲午夜理论影院| 国产精品香港三级国产av潘金莲| 日本一区二区免费在线视频| 日韩精品中文字幕看吧| 亚洲成人国产一区在线观看| 国产成人一区二区三区免费视频网站| 一二三四在线观看免费中文在| 国产精品自产拍在线观看55亚洲| 久久人妻av系列| 免费少妇av软件| 免费在线观看黄色视频的| 国产精品1区2区在线观看.| 久久精品人人爽人人爽视色| 老司机福利观看| 91字幕亚洲| 久久久久久久精品吃奶| 日韩大尺度精品在线看网址 | 91国产中文字幕| 精品国产国语对白av| 一级a爱片免费观看的视频| 一本久久中文字幕| 怎么达到女性高潮| 亚洲精品一区av在线观看| 久久久久久久精品吃奶| 国内精品久久久久精免费| 欧美+亚洲+日韩+国产| 成人手机av| 亚洲 欧美 日韩 在线 免费| 久久九九热精品免费| 一边摸一边做爽爽视频免费| 大香蕉久久成人网| 在线观看www视频免费| 丝袜人妻中文字幕| 999久久久精品免费观看国产| 亚洲五月婷婷丁香| 不卡一级毛片| 69精品国产乱码久久久| 精品电影一区二区在线| 国产精品秋霞免费鲁丝片| 亚洲av成人不卡在线观看播放网| 操出白浆在线播放| av天堂在线播放| 精品无人区乱码1区二区| 此物有八面人人有两片| 天天添夜夜摸| 午夜福利影视在线免费观看| 国产高清有码在线观看视频 | 人妻久久中文字幕网| 美女大奶头视频| 成人国语在线视频| 美女午夜性视频免费| 成人亚洲精品av一区二区| 男人舔女人下体高潮全视频| 国产片内射在线| 亚洲成人国产一区在线观看| 免费在线观看日本一区| 午夜免费观看网址| 人妻丰满熟妇av一区二区三区| 天堂√8在线中文| 国产精品电影一区二区三区| 精品久久久久久久毛片微露脸| 亚洲午夜理论影院| 免费女性裸体啪啪无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 亚洲第一欧美日韩一区二区三区| 午夜a级毛片| 亚洲国产中文字幕在线视频| 欧美老熟妇乱子伦牲交| 中文字幕久久专区| or卡值多少钱| 一本大道久久a久久精品| 18禁美女被吸乳视频| 欧美日韩瑟瑟在线播放| a级毛片在线看网站| 黄色成人免费大全| 在线观看免费日韩欧美大片| 纯流量卡能插随身wifi吗| 日韩有码中文字幕| av天堂在线播放| 亚洲一区高清亚洲精品| 精品一品国产午夜福利视频| 九色国产91popny在线| 亚洲熟妇中文字幕五十中出| 一区二区三区国产精品乱码| 国产av精品麻豆| 不卡av一区二区三区| 亚洲一码二码三码区别大吗| 国产视频一区二区在线看| 欧美成狂野欧美在线观看| 身体一侧抽搐| 日本撒尿小便嘘嘘汇集6| 91成人精品电影| 精品欧美一区二区三区在线| 免费av毛片视频| 精品一区二区三区av网在线观看| 色综合站精品国产| 99国产精品一区二区三区| 久久人妻福利社区极品人妻图片| 午夜精品在线福利| 亚洲 国产 在线| 国内精品久久久久久久电影| 亚洲欧美日韩无卡精品| svipshipincom国产片| 国产91精品成人一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 久久久久久久午夜电影| 免费在线观看视频国产中文字幕亚洲| 日本 欧美在线| 亚洲熟女毛片儿| 午夜精品久久久久久毛片777| 此物有八面人人有两片| 自拍欧美九色日韩亚洲蝌蚪91| 日本 av在线| 亚洲欧洲精品一区二区精品久久久| 黄片小视频在线播放| 中文字幕久久专区| 亚洲九九香蕉| 亚洲精品一卡2卡三卡4卡5卡| 99riav亚洲国产免费| 亚洲国产欧美日韩在线播放| 99香蕉大伊视频| 国产亚洲精品久久久久久毛片| 日本a在线网址| 亚洲九九香蕉| 成人18禁高潮啪啪吃奶动态图| 欧美日本亚洲视频在线播放| 国产精华一区二区三区| 国产精品99久久99久久久不卡| 精品免费久久久久久久清纯| 9191精品国产免费久久| 啦啦啦 在线观看视频| 九色国产91popny在线| av视频免费观看在线观看| 成人三级做爰电影| 国产在线观看jvid| 久久人妻av系列| 久久久久久国产a免费观看| 欧美日本视频| 免费搜索国产男女视频| 不卡一级毛片| 国产97色在线日韩免费| 色综合亚洲欧美另类图片| 久久青草综合色| 天天添夜夜摸| 一二三四在线观看免费中文在| 精品欧美国产一区二区三| 国产精品影院久久| √禁漫天堂资源中文www| 88av欧美| 女人爽到高潮嗷嗷叫在线视频| 亚洲一码二码三码区别大吗| 亚洲中文字幕一区二区三区有码在线看 | 国产1区2区3区精品| 精品国产乱码久久久久久男人| 国产精品98久久久久久宅男小说| 午夜精品久久久久久毛片777| 国产亚洲精品第一综合不卡| 成人18禁在线播放| 国产精品日韩av在线免费观看 | 亚洲成av人片免费观看| 此物有八面人人有两片| √禁漫天堂资源中文www| 亚洲精品美女久久av网站| 国产伦一二天堂av在线观看| 亚洲精品久久国产高清桃花| www.自偷自拍.com| www.精华液| 亚洲精品美女久久av网站| 无限看片的www在线观看| 午夜亚洲福利在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 波多野结衣av一区二区av| 国内精品久久久久久久电影| 村上凉子中文字幕在线| 久久精品国产99精品国产亚洲性色 | 91麻豆精品激情在线观看国产| 九色国产91popny在线| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 国产一区二区激情短视频| 国产成人免费无遮挡视频| 色综合站精品国产| 亚洲av片天天在线观看| 嫩草影院精品99| 好看av亚洲va欧美ⅴa在| 亚洲欧美一区二区三区黑人| 亚洲一区中文字幕在线| 亚洲欧美一区二区三区黑人| 免费无遮挡裸体视频| 免费看十八禁软件| 不卡一级毛片| 国产欧美日韩综合在线一区二区| 97碰自拍视频| 国产区一区二久久| 亚洲精品在线观看二区| 少妇的丰满在线观看| 色在线成人网| 亚洲精品国产一区二区精华液| 久久精品aⅴ一区二区三区四区| 高清在线国产一区| 亚洲成av片中文字幕在线观看| 午夜免费鲁丝| 久久久久久国产a免费观看| 亚洲一区二区三区不卡视频| 国产高清有码在线观看视频 | 啪啪无遮挡十八禁网站| 女同久久另类99精品国产91| 亚洲人成电影观看| 夜夜看夜夜爽夜夜摸| 国产精华一区二区三区| 一区二区三区高清视频在线| 丁香六月欧美| 露出奶头的视频| 精品久久久精品久久久| 满18在线观看网站| 曰老女人黄片| 桃红色精品国产亚洲av| 亚洲成人久久性| 性少妇av在线| 中文字幕人成人乱码亚洲影| √禁漫天堂资源中文www| 国产91精品成人一区二区三区| 淫妇啪啪啪对白视频| av有码第一页| 亚洲午夜理论影院| 久久人妻福利社区极品人妻图片| 亚洲最大成人中文| 1024视频免费在线观看| 丝袜在线中文字幕| 免费在线观看亚洲国产| 精品国产一区二区三区四区第35| 黄色丝袜av网址大全| 久久精品影院6| 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| 欧美激情久久久久久爽电影 | av免费在线观看网站| 两性夫妻黄色片| 18禁观看日本| 成人国语在线视频| 日韩欧美一区视频在线观看| www.精华液| 亚洲久久久国产精品| 午夜日韩欧美国产| 免费女性裸体啪啪无遮挡网站| 欧美国产精品va在线观看不卡| 久久久久久久精品吃奶| 亚洲三区欧美一区| 国产麻豆成人av免费视频| 日韩成人在线观看一区二区三区| 校园春色视频在线观看| 黄色女人牲交| xxx96com| 黄色视频,在线免费观看| 搞女人的毛片| 色播亚洲综合网| 亚洲人成网站在线播放欧美日韩| 国产成人av教育| 在线观看一区二区三区| 黄色a级毛片大全视频| 91成年电影在线观看| 欧美日本视频| av视频在线观看入口| 久久精品91蜜桃| 热99re8久久精品国产| 亚洲第一av免费看| 欧美国产日韩亚洲一区| 99热只有精品国产| 久久 成人 亚洲| 国产精品av久久久久免费| 色婷婷久久久亚洲欧美| 久久精品成人免费网站| 麻豆久久精品国产亚洲av| 欧美日韩福利视频一区二区| 日韩精品青青久久久久久| 性少妇av在线| av在线播放免费不卡| 亚洲天堂国产精品一区在线| 亚洲成av片中文字幕在线观看| 日日爽夜夜爽网站| 国产精品av久久久久免费| 一级a爱视频在线免费观看| 男人操女人黄网站| 亚洲自拍偷在线| 国产成人啪精品午夜网站| 国产私拍福利视频在线观看| 免费在线观看日本一区| 国产伦人伦偷精品视频| 少妇 在线观看| 在线观看免费午夜福利视频| 国产精品久久久av美女十八|