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

    魚(yú)眼鏡頭光學(xué)系統(tǒng)的非球面優(yōu)化設(shè)計(jì)

    2019-10-15 03:13:24呂麗軍
    應(yīng)用光學(xué) 2019年5期
    關(guān)鍵詞:優(yōu)化評(píng)價(jià)

    張 瀟,呂麗軍

    (上海大學(xué) 精密儀器系,上海 200444)

    引言

    魚(yú)眼鏡頭的視場(chǎng)角一般可以達(dá)到180°甚至270°[1],能同時(shí)獲取很大視場(chǎng)范圍的場(chǎng)景圖像。由于魚(yú)眼鏡頭系統(tǒng)大視場(chǎng)成像光束以很大斜入角照射在光學(xué)元件表面,偏離近軸成像,導(dǎo)致嚴(yán)重的像差。非球面光學(xué)元件可以有效地消除像差,提高成像質(zhì)量[2]。常用的確定非球面的方法是無(wú)像差法,但該方法為了實(shí)現(xiàn)穩(wěn)定的高精度測(cè)量,對(duì)人員要求和操作環(huán)境要求較高。利用三坐標(biāo)測(cè)量?jī)x獲得坐標(biāo),然后進(jìn)行最小二乘曲面擬合[3],從而降低了對(duì)人員和環(huán)境的要求,但計(jì)算量很大。研究如何選擇系統(tǒng)中哪個(gè)光學(xué)面為非球面及確定非球面面型是十分有意義的。

    對(duì)于魚(yú)眼鏡頭光學(xué)系統(tǒng),大視場(chǎng)物點(diǎn)發(fā)出的光束以很大的入射角度入射到前組負(fù)彎月形透鏡的表面,其成像具有平面對(duì)稱(chēng)的特性。賽德?tīng)栂癫罾碚摬辉龠m合此類(lèi)系統(tǒng)的分析[4-8]。非球面的光路計(jì)算考慮通過(guò)迭代逼近真實(shí)光路情況[9],缺陷在于很難提出一個(gè)適合于各種場(chǎng)合逼近精度的數(shù)值判據(jù)。目前,LU等應(yīng)用波像差的方法發(fā)展了平面對(duì)稱(chēng)光學(xué)系統(tǒng)的像差理論[10],并基于此理論發(fā)展了超大視場(chǎng)光學(xué)系統(tǒng)的優(yōu)化設(shè)計(jì)方法,最近,進(jìn)一步發(fā)展了適用于超大視場(chǎng)光學(xué)系統(tǒng)的六階波像差理論[11],可對(duì)大孔徑魚(yú)眼鏡頭系統(tǒng)進(jìn)行像差分析計(jì)算。本文在這些工作的基礎(chǔ)上,推導(dǎo)了二次圓錐曲面的面型系數(shù),建立了五階像差的評(píng)價(jià)函數(shù),針對(duì)一個(gè)全球面魚(yú)眼鏡頭,應(yīng)用六階波像差理論,分析各光學(xué)面的波像差分布以及它們隨非球面系數(shù)的變化趨勢(shì),從而選定應(yīng)用非球面的光學(xué)表面,并使用自適應(yīng)歸一化實(shí)數(shù)編碼遺傳算法(ANRCGA)[12]對(duì)其進(jìn)行優(yōu)化。

    本文簡(jiǎn)要介紹了六階波像差理論,推導(dǎo)了適用于該理論的面型系數(shù)公式,分析其與波像差的關(guān)系,并建立基于五階像差的評(píng)價(jià)函數(shù),以便以數(shù)值方式評(píng)價(jià)優(yōu)化設(shè)計(jì)。提出了如何為魚(yú)眼鏡頭光學(xué)系統(tǒng)選擇合適的光學(xué)面為非球面的方法,對(duì)所提出的方法進(jìn)行數(shù)值驗(yàn)證,并對(duì)結(jié)果進(jìn)行比較分析。

    1 波像差理論

    1.1 平面對(duì)稱(chēng)光學(xué)系統(tǒng)的像差計(jì)算

    對(duì)由g個(gè)光學(xué)面組成的平面對(duì)稱(chēng)光學(xué)系統(tǒng),系統(tǒng)的波像差是各光學(xué)面的波像差之和,即:

    (1)

    式中:W(k)、nk和wij0(k)分別表示第k個(gè)光學(xué)面的波像差、折射率和波像差系數(shù);xk和yk表示光線(xiàn)入射在第k個(gè)光學(xué)面上的孔徑坐標(biāo)。若將最后一個(gè)光學(xué)面上光線(xiàn)的孔徑坐標(biāo)(xg,yg)作為參考坐標(biāo),上式也可寫(xiě)成[11]:

    (2)

    式中:nk/g=nk/ng;Ak/g和gk/g分別為在第k和第g光學(xué)面之間x、y坐標(biāo)的線(xiàn)性傳遞系數(shù)。

    波像差理論研究超大視場(chǎng)光學(xué)系統(tǒng)時(shí),在全視場(chǎng)范圍內(nèi)取多個(gè)視場(chǎng)物點(diǎn),分別計(jì)算各視場(chǎng)物點(diǎn)的波像差,然后求和來(lái)反映光學(xué)系統(tǒng)在整個(gè)視場(chǎng)范圍的成像性能[10]。因此,對(duì)像差的貢獻(xiàn)占主導(dǎo)地位的波像差是球差波像差Wsph、彗差波像差Wcoma和場(chǎng)曲波像差Wcurv。其中場(chǎng)曲是通過(guò)光學(xué)面的出射波陣面和參考波陣面的幾何關(guān)系得出的[13],計(jì)算公式如下:

    Wsph=W400x4+W220x2y2+W040y4+W600x6+

    W420x4y2+W240x2y4+W060y6

    (3)

    Wcoma=W300x3+W120xy2+W500x5+

    W320x3y2+W140xy4

    (4)

    (5)

    一般魚(yú)眼鏡頭后光組的光瞳尺寸比前光組的大很多,因此,Wsph和Wcoma的主要貢獻(xiàn)來(lái)自于后光組光學(xué)面。另外,大視場(chǎng)物點(diǎn)發(fā)出的光束以較大的斜入射角入射到前組光學(xué)面上,所以場(chǎng)曲的主要貢獻(xiàn)來(lái)自于前光組的光學(xué)面。

    通過(guò)對(duì)波像差關(guān)于子午和弧矢方向分別求導(dǎo),得到光線(xiàn)在像面上的五階像差(x′,y′)表達(dá)式為

    (6)

    式中dij0和hij0為像差系數(shù)[11]。

    從上述表達(dá)式可以發(fā)現(xiàn),波像差和像差都受到面型系數(shù)ci,j的影響,因此可以通過(guò)cij控制波像差,提高成像質(zhì)量。魚(yú)眼鏡頭光學(xué)系統(tǒng)的孔徑一般較大,參考文獻(xiàn)[10]所述四階精度不夠精確,還需推導(dǎo)適用于六階波像差理論的面型系數(shù)。

    1.2 二次圓錐曲面的面型系數(shù)

    圖1 主光線(xiàn)經(jīng)過(guò)二次圓錐曲面反射時(shí)坐標(biāo)系統(tǒng)xyz的原點(diǎn)(主光線(xiàn)和光學(xué)面的交點(diǎn))坐標(biāo)Fig.1 Origin (intersection of chief ray and optical surface) coordinates of xyz coordinate system when chief ray is reflected by quadratic conical surface

    如圖1所示,在以光學(xué)面和光軸的交點(diǎn)O′為坐標(biāo)原點(diǎn)的坐標(biāo)系x′y′z′,繞z′軸旋轉(zhuǎn)的二次圓錐曲面方程為

    x′2+y′2=a1z′+a2z′2,a1=2R0

    (7)

    式中:R0是二次圓錐曲面在O′點(diǎn)處的曲率半徑;a2是二次圓錐曲面的非球面系數(shù),當(dāng)a2=-1時(shí)該曲面為球面,當(dāng)a2≠-1時(shí)為非球面。

    波像差理論中描述面型時(shí)采用以光線(xiàn)與光學(xué)面交點(diǎn)O為原點(diǎn)的坐標(biāo)系xyz,因此將常規(guī)的二次圓錐曲面方程式(1)轉(zhuǎn)換到坐標(biāo)系xyz中,x軸為原點(diǎn)O沿光學(xué)面的切線(xiàn),y軸為原點(diǎn)O沿光學(xué)面的法線(xiàn),y軸的斜率為tanθ。根據(jù)圖1所示,坐標(biāo)系x′y′z′和坐標(biāo)系xyz的轉(zhuǎn)化關(guān)系為

    (8)

    將(8)式帶入(7)式,并作泰勒級(jí)數(shù)展開(kāi),就可以得到適用于六階波像差理論的二次圓錐曲面方程:

    z=c2,0x2+c0,2y2+c3.0x3+c1,2xy2+c4,0x4+

    c0,4y4+c2,2x2y2+c5,0x5+c3,2x3y2+c1,4xy4+

    c6,0x6+c4,2x4y2+c2,4x2y4+c0,6y6

    (9)

    式中ci,j為各階對(duì)應(yīng)的面型系數(shù),具體表達(dá)式如下:

    c2,0=-a21B3c0,2=-1Bc3,0=-4Aa21x?0B6c1,2=-4Ax?0B4c4,0=-a21(a21C2+16A2X?20)B9c0,4=-C3B5c2,2=-2(a21C2+8A2x?20)B7c5,0=-4Aa21x?0(3a21C2+16A2x?20)B12c1,4=-12AC2x?0B8c3,2=-8Ax?0(3a21C2+8A2x?20)B10c6,0=-2a21[a41C4+16A2x?20(3a21C2+8A2x?20)]B12c4,2=-2[3a41C4+32A2x?20(3a21C2+4A2x?20)]B13C2,4=-6C2(a21C2+16A2x?20)B11c0,6=-2C2B9

    雖然波像差表達(dá)式中各項(xiàng)參數(shù)都受到面型系數(shù)ci,j的影響,但是每一個(gè)光學(xué)面的面型影響成像質(zhì)量的程度不一樣。從以上公式可知,ci,j受a1和非球面系數(shù)a2影響,a1=2R0,所以a2是影響ci,j的主要因素,通過(guò)a2改變ci,j,進(jìn)而影響波像差。

    2 利用非球面優(yōu)化魚(yú)眼鏡頭

    2.1 評(píng)價(jià)函數(shù)的定義

    為了對(duì)魚(yú)眼鏡頭光學(xué)系統(tǒng)進(jìn)行數(shù)值評(píng)價(jià),在視場(chǎng)范圍內(nèi)取k個(gè)視場(chǎng)物點(diǎn),基于它們?cè)谙衩嫔系南癫铧c(diǎn)列圖分布,定義魚(yú)眼鏡頭光學(xué)系統(tǒng)的像差評(píng)價(jià)函數(shù)為

    (10)

    其中,

    式中:εi表示不同視場(chǎng)角的權(quán)重因子;Wq和L分別為光束在最后光學(xué)面上沿子午(x)和弧矢(y)方向的投影長(zhǎng)度;x′和y′為光線(xiàn)在像面上的像差,可由(6)式計(jì)算得到。(10)式中積分域采用橢圓是因?yàn)榭讖焦怅@一般采用圓形孔,但光束在子午和弧矢平面內(nèi)的焦點(diǎn)位置不一樣,從而造成光束截面一般呈橢圓形。

    2.2 應(yīng)用非球面的優(yōu)化思路

    利用非球面優(yōu)化魚(yú)眼鏡頭的基本思路是,找到系統(tǒng)中對(duì)成像質(zhì)量影響最大的光學(xué)面,利用非球面系數(shù)a2減小其波像差。對(duì)于一個(gè)光學(xué)系統(tǒng)來(lái)說(shuō),光學(xué)面的波像差越大,其各類(lèi)光學(xué)參量對(duì)加工制造的精度要求越高[14]。因此選擇波像差最大的光學(xué)面為非球面,在減小波像差、提高成像質(zhì)量的同時(shí),也有助于增大非球面的可允許公差范圍,降低非球面的制造難度。

    對(duì)魚(yú)眼鏡頭系統(tǒng)各視場(chǎng)的物點(diǎn)而言,主要影響成像質(zhì)量的是各光學(xué)面的球差波像差Wsph和彗差波像差Wcoma及前光組的場(chǎng)曲波像差Wcurv。大孔徑時(shí)各光學(xué)面的Wsph和Wcoma對(duì)像差的貢獻(xiàn)遠(yuǎn)大于Wcurv,分析波像差分布,選擇波像差較大的光學(xué)面為非球面。若幾個(gè)光學(xué)面的波像差均較大,則需進(jìn)一步結(jié)合系統(tǒng)總的波像差情況,分析各光學(xué)面Wsph和Wcoma隨非球面系數(shù)a2的變化趨勢(shì)。假設(shè)系統(tǒng)總的球差和彗差波像差符號(hào)相反,應(yīng)該選擇Wsph和Wcoma隨a2的變化趨勢(shì)相反的光學(xué)面;若系統(tǒng)總的球差和彗差波像差符號(hào)相同,則應(yīng)選擇Wsph和Wcoma隨a2變化相同的光學(xué)面。另外,如果有多個(gè)光學(xué)面符合上述條件,則選擇隨a2變化最迅速的光學(xué)面,以減小加工難度。小孔徑時(shí),Wcurv對(duì)像差的貢獻(xiàn)大于各光學(xué)面的Wsph和Wcoma。結(jié)合之前的研究[15],本文認(rèn)為選擇第1和第3光學(xué)面為非球面可以有效校正場(chǎng)曲。

    確定非球面后,利用評(píng)價(jià)函數(shù)式(9)和自適應(yīng)歸一化實(shí)數(shù)編碼遺傳算法(ANRCGA)[12],得到評(píng)價(jià)函數(shù)最優(yōu)時(shí)的非球面系數(shù)值及相應(yīng)的光學(xué)系統(tǒng)參數(shù)。

    3 數(shù)值驗(yàn)證

    3.1 魚(yú)眼鏡頭實(shí)例

    應(yīng)用上述方法,對(duì)專(zhuān)利文獻(xiàn)[16]中的魚(yú)眼鏡頭進(jìn)行優(yōu)化設(shè)計(jì),其光路如圖2所示,工作波段為可見(jiàn)光,最大半視場(chǎng)角ω0=80°。分別進(jìn)行如下5種情況的分析比對(duì):1) 對(duì)專(zhuān)利設(shè)計(jì)進(jìn)行全球面優(yōu)化設(shè)計(jì);2) 取較大孔徑(F#3.3)時(shí)后光組的非球面設(shè)計(jì);3) 取較小孔徑(F#7.6)時(shí)后光組的非球面設(shè)計(jì);4) 取較大孔徑(F#3.3)時(shí),前光組第1、3光學(xué)面的非球面優(yōu)化;5)取較小孔徑(F#7.6)時(shí),前光組第1、3光學(xué)面的非球面優(yōu)化。表1為上述5種情況優(yōu)化后的光學(xué)參數(shù),其中:rk表示光學(xué)面的曲率半徑;dk表示鏡片間距。表2為5種設(shè)計(jì)對(duì)應(yīng)的非球面所在光學(xué)面及系統(tǒng)接收孔徑。

    圖2 魚(yú)眼鏡頭光學(xué)系統(tǒng)光路圖Fig. 2 Optical path of fisheye lens

    表1 魚(yú)眼鏡頭Design x的光學(xué)參量 mm

    表2 Design x的非球面所在光學(xué)面及a2

    根據(jù)以上所述,由(3)式計(jì)算專(zhuān)利原始設(shè)計(jì)0°,10°,…,80°共9個(gè)視場(chǎng)的波像差Wcurv、Wsph和Wcoma,取ω0=80°時(shí)波像差的分布,如圖3所示。其中圖3 (a)和圖3(c)分別為F#3.3和F#7.6時(shí)各光學(xué)面波像差分布情況,圖3(b)為系統(tǒng)的前組場(chǎng)曲波像差?!苉為系統(tǒng)總的球差波像差∑Wsph和彗差波像差∑Wcoma。由于系統(tǒng)前5個(gè)光學(xué)面的Wsph和Wcoma均小于0.01,圖3中未列出。

    圖4(a)、(b)和(c)分別為當(dāng)?shù)?7、18和21面為非球面時(shí),該光學(xué)面的Wsph和Wcoma隨非球面系數(shù)a2的變化趨勢(shì)。

    圖3 魚(yú)眼鏡頭的波像差分布Fig. 3 Wave aberration distribution of fisheye lens

    圖4 Wsph和Wcoma隨第k個(gè)光學(xué)面a2的變化趨勢(shì)圖Fig. 4 Changing trend of the k-th surface’s Wsph和 Wcoma with aspherical coefficient

    根據(jù)圖3和圖4可知:

    1) 第8、9、17、18和21光學(xué)面的波像差比較突出,但是第8和第9面曲率半徑相似且空氣間隔很小,兩者的波像差大小基本相同,符號(hào)相反,可以相互抵消,不納入非球面的考慮范圍。

    2) 由于系統(tǒng)的∑Wsph和∑Wcoma同號(hào),從優(yōu)化角度考慮,要求光學(xué)面Wsph和Wcoma隨 變化趨勢(shì)相同。在第17、18和19面中,只有第17面符合條件。

    3) 當(dāng)大孔徑F#3.3時(shí),Wcurv相對(duì)較小,應(yīng)優(yōu)先校正Wsph和Wcoma,選擇第17面為非球面;當(dāng)小孔徑F#7.6時(shí),Wcurv的作用顯著,應(yīng)優(yōu)先校正前光組場(chǎng)曲,選擇第1和3面為非球面。

    3.2 優(yōu)化與討論

    確定待優(yōu)化的光學(xué)面后,對(duì)5種情況依次應(yīng)用(9)式計(jì)算出對(duì)應(yīng)的評(píng)價(jià)函數(shù),令εi=1,將該光學(xué)面的非球面系數(shù)a2遍歷[-4,4](涵蓋了常見(jiàn)的易于加工的非球面類(lèi)型[17]),變化步距為0.01。選擇使評(píng)價(jià)函數(shù)達(dá)到最小時(shí)的a2為相應(yīng)的非球面系數(shù)初值。

    基于評(píng)價(jià)函數(shù),將(9)式中60°、70°和80°這3個(gè)視場(chǎng)角的權(quán)重因子εi取為5,其他權(quán)重因子εi取為1。然后應(yīng)用自適應(yīng)歸一化實(shí)數(shù)編碼遺傳算法(ANRCGA)[12]分別對(duì)5種情況進(jìn)行全局優(yōu)化。種群數(shù)為300,每一代種群規(guī)模為100,交叉概率為0.85。采取粗、精兩輪優(yōu)化,粗優(yōu)20次,各光學(xué)參量的優(yōu)化搜索范圍為其初值的±50%;從中選擇評(píng)價(jià)函數(shù)最小的5次作為精優(yōu)的初值,搜索范圍取初值的±20%,所有視場(chǎng)的物點(diǎn)均為無(wú)限遠(yuǎn)。圖5為優(yōu)化后5種設(shè)計(jì)在ω0=0°,ω0=40°和ω0=80°時(shí)的點(diǎn)列圖,圖中坐標(biāo)為Image width,單位為μm。表3為它們的評(píng)價(jià)函數(shù)Q,圖6為Design 2和Design5的波像差分布。

    圖5 不同視場(chǎng)角時(shí)最后成像面上的孔徑光線(xiàn)像差(最右側(cè)一列參數(shù)表示視場(chǎng)角)Fig. 5 Aperture ray aberrations on final image plane at different field angles (The far right column parameter means field angles)

    圖6 Design 2和Design 5波像差分布圖Fig. 6 Wave aberration distribution of Design 2 and 5

    表3 魚(yú)眼鏡頭的評(píng)價(jià)函數(shù)值

    分析圖5和表3可以得出如下結(jié)論:

    1) 4種非球面優(yōu)化設(shè)計(jì)比全球面優(yōu)化設(shè)計(jì)Design 1的點(diǎn)列圖的范圍減小1/2,評(píng)價(jià)函數(shù)值提高了約1個(gè)數(shù)量級(jí),說(shuō)明ANRCGA算法可用于魚(yú)眼鏡頭的非球面優(yōu)化,且效果良好。

    2) 小孔徑時(shí)的Design 4和Design 5比大孔徑時(shí)的Design 2和Design 3的評(píng)價(jià)函數(shù)值小。說(shuō)明本文建立的評(píng)價(jià)函數(shù)可以準(zhǔn)確反映魚(yú)眼鏡頭光學(xué)系統(tǒng)的成像質(zhì)量。

    3) Design 2的結(jié)果優(yōu)于Design 3,說(shuō)明在大孔徑系統(tǒng)中,按照各光學(xué)面的波像差分布,選擇波像差較大且對(duì)非球面系數(shù)敏感的光學(xué)面為非球面,可以有效校正系統(tǒng)的波像差,并改善鏡頭成像質(zhì)量。

    4) Design 5的結(jié)果優(yōu)于Design 4,說(shuō)明隨著孔徑的減小,各光學(xué)面的Wsph和Wcoma減小,Wcurv的作用顯著。此時(shí)Design 4優(yōu)化Wsph和Wcoma的效果不如Design 5優(yōu)化前光組Wcurv,因此在小孔徑光學(xué)系統(tǒng)中應(yīng)該優(yōu)先校正場(chǎng)曲波像差。

    5) 當(dāng)孔徑較大時(shí),Design 2優(yōu)化效果顯著;孔徑較小時(shí),Design 5優(yōu)化效果顯著。比較圖6與圖4可以看出,經(jīng)非球面優(yōu)化后Design 2的Wsph和Wcoma大幅降低,Design 5的Wcurv大幅降低,說(shuō)明本文提出的優(yōu)化方法有效。

    4 結(jié)論

    推導(dǎo)了適用于六階波像差理論的面型方程,分析了波像差與光學(xué)面面型的關(guān)系,提出了一種有效選擇非球面的方法。對(duì)驗(yàn)證實(shí)例系統(tǒng)來(lái)說(shuō),當(dāng)系統(tǒng)的F#7.6,前組的場(chǎng)曲波像差較大時(shí),選擇對(duì)場(chǎng)曲波像差影響較大的第1和3光學(xué)面為非球面;當(dāng)系統(tǒng)的F#3.3,各光學(xué)面的球差和彗差波像差較大,根據(jù)系統(tǒng)的波像差分布情況,選擇波像差較大,且隨a2的變化趨勢(shì)有助于系統(tǒng)兩類(lèi)總的波像差減小的光學(xué)面(第17面)為非球面。通過(guò)實(shí)例證明,本文提出的方法可以快速選出非球面,對(duì)利用非球面優(yōu)化魚(yú)眼鏡頭有借鑒意義。

    猜你喜歡
    優(yōu)化評(píng)價(jià)
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    SBR改性瀝青的穩(wěn)定性評(píng)價(jià)
    石油瀝青(2021年4期)2021-10-14 08:50:44
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    中藥治療室性早搏系統(tǒng)評(píng)價(jià)再評(píng)價(jià)
    基于Moodle的學(xué)習(xí)評(píng)價(jià)
    基于低碳物流的公路運(yùn)輸優(yōu)化
    關(guān)于項(xiàng)目后評(píng)價(jià)中“專(zhuān)項(xiàng)”后評(píng)價(jià)的探討
    国产一区二区在线观看av| 免费在线观看影片大全网站 | 午夜福利免费观看在线| 婷婷色综合大香蕉| 一级毛片女人18水好多 | 如日韩欧美国产精品一区二区三区| 精品欧美一区二区三区在线| 午夜福利影视在线免费观看| 在线观看一区二区三区激情| 免费高清在线观看日韩| 视频区图区小说| 久久99热这里只频精品6学生| 一级a爱视频在线免费观看| 日本欧美视频一区| 男的添女的下面高潮视频| 最新在线观看一区二区三区 | 午夜福利,免费看| 日韩av在线免费看完整版不卡| av网站在线播放免费| netflix在线观看网站| 亚洲精品国产av蜜桃| 欧美黑人欧美精品刺激| 亚洲五月婷婷丁香| 啦啦啦中文免费视频观看日本| 国产精品国产三级国产专区5o| 一区福利在线观看| 免费看十八禁软件| 欧美日韩视频高清一区二区三区二| 国产视频一区二区在线看| av欧美777| 亚洲av日韩在线播放| 亚洲成国产人片在线观看| 亚洲国产欧美网| 国产爽快片一区二区三区| 一二三四社区在线视频社区8| 男女边摸边吃奶| 亚洲欧洲日产国产| 日韩中文字幕欧美一区二区 | 午夜福利视频精品| 人人妻,人人澡人人爽秒播 | 色播在线永久视频| 国产亚洲av片在线观看秒播厂| 亚洲av日韩精品久久久久久密 | 国产亚洲一区二区精品| 久久鲁丝午夜福利片| 如日韩欧美国产精品一区二区三区| 精品第一国产精品| 悠悠久久av| 老汉色av国产亚洲站长工具| 国产成人一区二区三区免费视频网站 | 97精品久久久久久久久久精品| 亚洲精品久久成人aⅴ小说| 国产女主播在线喷水免费视频网站| 大码成人一级视频| a级毛片黄视频| 亚洲国产看品久久| 国产精品国产三级专区第一集| 国产视频首页在线观看| 美女视频免费永久观看网站| 国产高清不卡午夜福利| 啦啦啦在线观看免费高清www| 又黄又粗又硬又大视频| 欧美黑人精品巨大| 欧美精品一区二区免费开放| 真人做人爱边吃奶动态| 日本91视频免费播放| 美女脱内裤让男人舔精品视频| 少妇人妻久久综合中文| 色婷婷久久久亚洲欧美| 国产精品久久久人人做人人爽| 久久久久久久久久久久大奶| 亚洲情色 制服丝袜| 啦啦啦啦在线视频资源| 国产免费现黄频在线看| 亚洲国产欧美在线一区| 黄色视频在线播放观看不卡| 国产精品二区激情视频| 性少妇av在线| 国产片特级美女逼逼视频| 国产淫语在线视频| 亚洲五月婷婷丁香| 丝袜美足系列| 国产xxxxx性猛交| 国产女主播在线喷水免费视频网站| 国产激情久久老熟女| 亚洲欧美中文字幕日韩二区| 人妻一区二区av| 国产黄色免费在线视频| 啦啦啦在线免费观看视频4| 国产爽快片一区二区三区| 老司机影院毛片| 女人精品久久久久毛片| 亚洲精品一卡2卡三卡4卡5卡 | 成人免费观看视频高清| 女人精品久久久久毛片| 一边摸一边抽搐一进一出视频| 日韩欧美一区视频在线观看| 欧美日韩视频精品一区| 欧美乱码精品一区二区三区| 久久精品aⅴ一区二区三区四区| 麻豆av在线久日| 国产高清国产精品国产三级| 伊人久久大香线蕉亚洲五| 天天添夜夜摸| 男女边吃奶边做爰视频| 黄色 视频免费看| 国产精品九九99| 日本wwww免费看| 咕卡用的链子| 亚洲国产欧美日韩在线播放| av天堂在线播放| 九草在线视频观看| 黄频高清免费视频| 久久精品亚洲av国产电影网| 国产精品一国产av| 欧美日韩视频高清一区二区三区二| 不卡av一区二区三区| 欧美日韩黄片免| 精品一品国产午夜福利视频| 国产福利在线免费观看视频| 亚洲国产av影院在线观看| 久久ye,这里只有精品| netflix在线观看网站| 亚洲成人免费电影在线观看 | 欧美黄色淫秽网站| 日日爽夜夜爽网站| 老司机靠b影院| kizo精华| 亚洲av国产av综合av卡| av网站在线播放免费| 美国免费a级毛片| www.自偷自拍.com| 91麻豆精品激情在线观看国产 | 亚洲国产欧美网| 久久久久久免费高清国产稀缺| 国产精品久久久av美女十八| 色94色欧美一区二区| 少妇人妻 视频| 免费在线观看日本一区| 狂野欧美激情性bbbbbb| 精品亚洲成a人片在线观看| 国产一区二区三区综合在线观看| 亚洲黑人精品在线| 最新在线观看一区二区三区 | 色播在线永久视频| 九色亚洲精品在线播放| 国产在线视频一区二区| 777米奇影视久久| 一级a爱视频在线免费观看| 亚洲,欧美,日韩| 国产精品秋霞免费鲁丝片| 最近手机中文字幕大全| 国产精品二区激情视频| 午夜激情久久久久久久| 三上悠亚av全集在线观看| 99国产综合亚洲精品| 一边亲一边摸免费视频| 日韩制服丝袜自拍偷拍| 欧美亚洲 丝袜 人妻 在线| 日日摸夜夜添夜夜爱| 久久久国产欧美日韩av| 中文欧美无线码| 中文字幕高清在线视频| 少妇被粗大的猛进出69影院| 极品人妻少妇av视频| 欧美日韩av久久| 99国产精品免费福利视频| 亚洲伊人久久精品综合| 亚洲欧美激情在线| 国产精品.久久久| 国产成人欧美在线观看 | 黄色片一级片一级黄色片| 日本欧美国产在线视频| 黄色视频不卡| 国产精品 欧美亚洲| 欧美 日韩 精品 国产| 肉色欧美久久久久久久蜜桃| 久久午夜综合久久蜜桃| 搡老岳熟女国产| 在线观看免费日韩欧美大片| 热99国产精品久久久久久7| 91麻豆精品激情在线观看国产 | 一级毛片我不卡| 午夜免费观看性视频| 欧美人与性动交α欧美精品济南到| 高清欧美精品videossex| 久久天堂一区二区三区四区| 大码成人一级视频| 男人舔女人的私密视频| 永久免费av网站大全| 日本av手机在线免费观看| 亚洲国产毛片av蜜桃av| 最新的欧美精品一区二区| 99国产精品99久久久久| 亚洲成人国产一区在线观看 | 午夜福利乱码中文字幕| 日本vs欧美在线观看视频| 水蜜桃什么品种好| 叶爱在线成人免费视频播放| 人人妻人人添人人爽欧美一区卜| 精品人妻在线不人妻| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| 午夜福利,免费看| 十分钟在线观看高清视频www| 一级黄片播放器| 黑丝袜美女国产一区| 多毛熟女@视频| 18禁观看日本| 亚洲人成电影观看| 在线 av 中文字幕| 女性生殖器流出的白浆| 久久ye,这里只有精品| 人人妻,人人澡人人爽秒播 | 国产亚洲av高清不卡| 黑丝袜美女国产一区| 国产精品一区二区在线不卡| 国产真人三级小视频在线观看| 91国产中文字幕| 又黄又粗又硬又大视频| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 久久久久久久国产电影| 又大又爽又粗| 欧美精品啪啪一区二区三区 | 满18在线观看网站| 亚洲av片天天在线观看| 国产欧美日韩精品亚洲av| 亚洲国产欧美在线一区| 女人高潮潮喷娇喘18禁视频| 如日韩欧美国产精品一区二区三区| 一级毛片女人18水好多 | 亚洲欧美清纯卡通| 国产欧美日韩一区二区三区在线| 久久国产精品大桥未久av| 国产成人欧美在线观看 | 免费看不卡的av| 别揉我奶头~嗯~啊~动态视频 | 另类精品久久| 欧美激情 高清一区二区三区| 免费在线观看黄色视频的| 亚洲欧美一区二区三区久久| 精品亚洲成a人片在线观看| 国产精品久久久久久人妻精品电影 | 97在线人人人人妻| 叶爱在线成人免费视频播放| 久久精品国产亚洲av高清一级| 中文精品一卡2卡3卡4更新| 超碰成人久久| 国产又爽黄色视频| 国产成人精品无人区| 欧美精品av麻豆av| 亚洲三区欧美一区| 一级片免费观看大全| netflix在线观看网站| 亚洲男人天堂网一区| 777米奇影视久久| 韩国精品一区二区三区| 美女主播在线视频| 另类精品久久| a级毛片在线看网站| 国产在视频线精品| 19禁男女啪啪无遮挡网站| 精品一区在线观看国产| h视频一区二区三区| 亚洲av男天堂| 老熟女久久久| bbb黄色大片| 久久久精品区二区三区| 日韩大码丰满熟妇| 亚洲欧美一区二区三区国产| 大话2 男鬼变身卡| 婷婷色综合www| 女性生殖器流出的白浆| 日日夜夜操网爽| 18禁国产床啪视频网站| 黑人猛操日本美女一级片| 日本91视频免费播放| 国产欧美日韩综合在线一区二区| 欧美成人午夜精品| 国产黄色免费在线视频| 九色亚洲精品在线播放| 中文字幕另类日韩欧美亚洲嫩草| av国产久精品久网站免费入址| 一二三四社区在线视频社区8| 欧美另类一区| 80岁老熟妇乱子伦牲交| 大片电影免费在线观看免费| 国产又爽黄色视频| 精品久久蜜臀av无| 搡老岳熟女国产| 日韩制服骚丝袜av| 天天躁夜夜躁狠狠躁躁| 午夜av观看不卡| 狂野欧美激情性bbbbbb| a级毛片在线看网站| 热99国产精品久久久久久7| 亚洲av电影在线进入| 一级毛片女人18水好多 | 亚洲精品av麻豆狂野| 欧美黄色片欧美黄色片| 久久热在线av| 欧美黄色片欧美黄色片| 精品第一国产精品| 日本五十路高清| 久久免费观看电影| 9热在线视频观看99| 黄频高清免费视频| 免费在线观看视频国产中文字幕亚洲 | 波野结衣二区三区在线| 久久精品久久精品一区二区三区| 亚洲久久久国产精品| 国产高清国产精品国产三级| 在线观看免费视频网站a站| 狂野欧美激情性xxxx| 亚洲av成人不卡在线观看播放网 | 亚洲欧美中文字幕日韩二区| cao死你这个sao货| 手机成人av网站| 黄色视频不卡| 免费一级毛片在线播放高清视频 | 18禁国产床啪视频网站| 伦理电影免费视频| 国产精品久久久人人做人人爽| 精品久久蜜臀av无| 精品国产一区二区三区四区第35| 成年女人毛片免费观看观看9 | 精品熟女少妇八av免费久了| 中文字幕制服av| 欧美日韩视频精品一区| 免费在线观看影片大全网站 | 国产精品99久久99久久久不卡| 黄频高清免费视频| 在现免费观看毛片| 国产成人精品久久久久久| 精品久久久精品久久久| 中文字幕另类日韩欧美亚洲嫩草| 2021少妇久久久久久久久久久| 国产精品一区二区在线观看99| 亚洲精品av麻豆狂野| 老汉色av国产亚洲站长工具| 高清不卡的av网站| 久久久国产精品麻豆| 国精品久久久久久国模美| 国产不卡av网站在线观看| 国产一区亚洲一区在线观看| 久久国产精品男人的天堂亚洲| 成人国产av品久久久| 欧美精品啪啪一区二区三区 | 最新的欧美精品一区二区| 亚洲精品一卡2卡三卡4卡5卡 | 热re99久久国产66热| 啦啦啦 在线观看视频| 国产激情久久老熟女| 久久狼人影院| 欧美国产精品一级二级三级| 1024香蕉在线观看| 黑人欧美特级aaaaaa片| 日韩,欧美,国产一区二区三区| 校园人妻丝袜中文字幕| 欧美黄色片欧美黄色片| 久久久久久久精品精品| 看免费成人av毛片| 国产1区2区3区精品| 国产三级黄色录像| 久久精品国产综合久久久| 黑丝袜美女国产一区| 久久精品国产综合久久久| 老司机影院毛片| av国产久精品久网站免费入址| 国产成人精品在线电影| 最近手机中文字幕大全| 亚洲国产精品一区二区三区在线| 我要看黄色一级片免费的| 狠狠婷婷综合久久久久久88av| 一本大道久久a久久精品| 免费观看人在逋| 一区二区三区乱码不卡18| 蜜桃在线观看..| 99热网站在线观看| 国产又爽黄色视频| 亚洲成国产人片在线观看| 超碰97精品在线观看| 亚洲精品自拍成人| 欧美日韩黄片免| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| 又大又黄又爽视频免费| 97在线人人人人妻| 婷婷丁香在线五月| 另类亚洲欧美激情| 黑人欧美特级aaaaaa片| 美女高潮到喷水免费观看| 黄色毛片三级朝国网站| 叶爱在线成人免费视频播放| 不卡av一区二区三区| 叶爱在线成人免费视频播放| 亚洲国产最新在线播放| 久久人妻熟女aⅴ| 久久久久网色| 国产精品国产三级专区第一集| 性色av一级| 后天国语完整版免费观看| 丝袜在线中文字幕| 极品人妻少妇av视频| 在线观看人妻少妇| 一级片'在线观看视频| 日韩精品免费视频一区二区三区| 欧美精品av麻豆av| 欧美日本中文国产一区发布| av国产久精品久网站免费入址| 精品国产乱码久久久久久男人| 欧美精品啪啪一区二区三区 | 成人手机av| 每晚都被弄得嗷嗷叫到高潮| 性少妇av在线| 久久狼人影院| 亚洲情色 制服丝袜| 亚洲成av片中文字幕在线观看| 国产av国产精品国产| 老汉色av国产亚洲站长工具| 99久久99久久久精品蜜桃| 人体艺术视频欧美日本| 丝袜脚勾引网站| 亚洲一码二码三码区别大吗| 曰老女人黄片| 成人手机av| 三上悠亚av全集在线观看| 黄色片一级片一级黄色片| 久久精品国产a三级三级三级| 首页视频小说图片口味搜索 | 久久热在线av| 日韩熟女老妇一区二区性免费视频| 9191精品国产免费久久| 免费观看av网站的网址| 精品少妇一区二区三区视频日本电影| 青草久久国产| 狂野欧美激情性bbbbbb| 性色av乱码一区二区三区2| 成人18禁高潮啪啪吃奶动态图| 啦啦啦中文免费视频观看日本| 乱人伦中国视频| 欧美亚洲 丝袜 人妻 在线| 咕卡用的链子| 国产一区二区三区综合在线观看| 国产精品人妻久久久影院| 伊人久久大香线蕉亚洲五| 国产黄频视频在线观看| 黄色视频在线播放观看不卡| 香蕉丝袜av| 中文字幕高清在线视频| 国产精品秋霞免费鲁丝片| 侵犯人妻中文字幕一二三四区| 母亲3免费完整高清在线观看| 一边摸一边做爽爽视频免费| 日本av免费视频播放| 国产深夜福利视频在线观看| 2018国产大陆天天弄谢| 亚洲成av片中文字幕在线观看| 国产成人精品久久二区二区免费| 天堂中文最新版在线下载| 97精品久久久久久久久久精品| 日本色播在线视频| 新久久久久国产一级毛片| 午夜久久久在线观看| 久久av网站| 中文字幕高清在线视频| 交换朋友夫妻互换小说| 热re99久久精品国产66热6| 日本一区二区免费在线视频| 国产真人三级小视频在线观看| 国产精品一区二区在线不卡| 午夜福利免费观看在线| 手机成人av网站| 久久九九热精品免费| av线在线观看网站| 国产精品三级大全| 国产欧美亚洲国产| 日本欧美国产在线视频| 久久九九热精品免费| av天堂在线播放| 亚洲欧洲日产国产| www.av在线官网国产| 免费在线观看视频国产中文字幕亚洲 | 亚洲av日韩精品久久久久久密 | 成人亚洲欧美一区二区av| 一边摸一边抽搐一进一出视频| 久久九九热精品免费| 日本五十路高清| 国产精品欧美亚洲77777| 国产欧美日韩精品亚洲av| 最近手机中文字幕大全| 久久99精品国语久久久| 国产日韩欧美亚洲二区| av天堂久久9| 国产成人精品在线电影| 九色亚洲精品在线播放| 18禁黄网站禁片午夜丰满| 搡老岳熟女国产| 黄片播放在线免费| 精品国产一区二区三区久久久樱花| 国产精品久久久久久精品电影小说| 国产一区二区三区av在线| 日本五十路高清| 男女国产视频网站| 久久久国产精品麻豆| 丝袜脚勾引网站| 中文字幕人妻熟女乱码| 美女中出高潮动态图| 一级毛片女人18水好多 | 国产亚洲一区二区精品| 亚洲免费av在线视频| 亚洲视频免费观看视频| 亚洲欧美一区二区三区久久| 18在线观看网站| 亚洲伊人色综图| 一级毛片电影观看| 午夜福利免费观看在线| 97精品久久久久久久久久精品| bbb黄色大片| 午夜免费观看性视频| 久久 成人 亚洲| 精品福利永久在线观看| 欧美日韩亚洲高清精品| 美女扒开内裤让男人捅视频| 亚洲国产日韩一区二区| 亚洲午夜精品一区,二区,三区| 久9热在线精品视频| 国产xxxxx性猛交| 一级毛片女人18水好多 | 蜜桃国产av成人99| 久久久久国产一级毛片高清牌| 高清视频免费观看一区二区| 一本大道久久a久久精品| 久久九九热精品免费| 欧美av亚洲av综合av国产av| 高潮久久久久久久久久久不卡| 精品免费久久久久久久清纯 | 大码成人一级视频| 捣出白浆h1v1| 亚洲欧美中文字幕日韩二区| 黄色片一级片一级黄色片| 久久久精品免费免费高清| 精品少妇内射三级| 一级,二级,三级黄色视频| 亚洲av成人不卡在线观看播放网 | 黄网站色视频无遮挡免费观看| 日韩制服骚丝袜av| 亚洲成人国产一区在线观看 | 免费人妻精品一区二区三区视频| 丁香六月天网| 色视频在线一区二区三区| 最近最新中文字幕大全免费视频 | 91麻豆av在线| 国产精品九九99| 亚洲欧美精品综合一区二区三区| 日韩制服骚丝袜av| 欧美激情 高清一区二区三区| 无限看片的www在线观看| 制服诱惑二区| avwww免费| 国产成人a∨麻豆精品| 一级片免费观看大全| 久久精品国产a三级三级三级| av天堂久久9| 国产亚洲精品久久久久5区| 老司机靠b影院| www.精华液| 国产男女超爽视频在线观看| 久久久精品区二区三区| 色综合欧美亚洲国产小说| 在线天堂中文资源库| 下体分泌物呈黄色| 老司机午夜十八禁免费视频| 亚洲av美国av| 曰老女人黄片| 日韩一卡2卡3卡4卡2021年| 免费看不卡的av| 中文字幕精品免费在线观看视频| 久9热在线精品视频| 精品人妻在线不人妻| 欧美乱码精品一区二区三区| 午夜av观看不卡| 午夜免费观看性视频| 天天躁日日躁夜夜躁夜夜| 两性夫妻黄色片| 久久狼人影院| netflix在线观看网站| 国产一区二区三区av在线| avwww免费| 丰满少妇做爰视频| 后天国语完整版免费观看| 亚洲精品久久成人aⅴ小说| 亚洲欧美一区二区三区久久| 亚洲色图 男人天堂 中文字幕| av欧美777| 丰满少妇做爰视频| 免费观看a级毛片全部| 成人午夜精彩视频在线观看| 国产在视频线精品| 亚洲人成电影免费在线| 另类亚洲欧美激情| 不卡av一区二区三区| 午夜老司机福利片| 人人妻人人爽人人添夜夜欢视频| 久久久久网色| 性少妇av在线| 日本黄色日本黄色录像| 侵犯人妻中文字幕一二三四区| 狂野欧美激情性xxxx| 欧美日韩成人在线一区二区| 波多野结衣av一区二区av| 免费一级毛片在线播放高清视频 | 国产精品 国内视频| 日本91视频免费播放|