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

        ?

        紫銅粗糙表面的光譜雙向反射分布函數(shù)測(cè)量研究

        2020-01-10 02:32:20楊博文馬王杰慧劉彥磊劉玉芳
        中國(guó)光學(xué) 2019年6期
        關(guān)鍵詞:鏡面反射紫銅頂角

        楊博文,馬王杰慧,劉彥磊,劉玉芳

        (河南師范大學(xué)物理與材料科學(xué)學(xué)院,河南 新鄉(xiāng) 453007)

        1 引 言

        雙向反射分布函數(shù)(Bidirectional Reflectance Distribution Function,BRDF)是表征目標(biāo)表面空間光學(xué)反射特性的物理量,能夠有效描述目標(biāo)表面的光學(xué)反射特性[1-2]。BRDF在目標(biāo)光散射計(jì)算、雜光分析、計(jì)算機(jī)圖像處理、地物遙感、目標(biāo)識(shí)別和材料的光學(xué)性能研究等領(lǐng)域有著廣泛應(yīng)用[3-7]。因此,近幾十年來關(guān)于BRDF的研究工作從未間斷過,這些工作主要集中在BRDF測(cè)量裝置的研制[6,8-12]、材料BRDF的測(cè)量、BRDF模型的建立和BRDF的應(yīng)用等方面[13-15]。準(zhǔn)確描述目標(biāo)表面的視覺特性是BRDF的重要應(yīng)用之一。眾所周知,目標(biāo)的顏色與光在介質(zhì)中的吸收和散射特性有關(guān),光澤則與光的反射方向和強(qiáng)度有關(guān)[16],光譜BRDF能夠很好地描述目標(biāo)的視覺特性,對(duì)工藝品視覺效果方面的設(shè)計(jì)具有重要意義。

        紫銅是一種常見的金屬材料,由于其具有良好的導(dǎo)電性能被廣泛應(yīng)用于電氣工業(yè)。紫銅具有暗紅的金屬光澤,不僅能夠呈現(xiàn)出現(xiàn)代感,還能體現(xiàn)出沉穩(wěn)、高貴的品質(zhì),是銅飾中應(yīng)用最廣泛的材料。測(cè)量紫銅表面的光譜BRDF能夠提供紫銅的光學(xué)反射特性,為紫銅工藝品加工設(shè)計(jì)提供準(zhǔn)確的視覺特性,有助于降低加工成本和優(yōu)化設(shè)計(jì)效果。

        關(guān)于銅表面BRDF的測(cè)量研究很少見到報(bào)道,戴景民等[17]測(cè)量了25~500 ℃的溫度范圍內(nèi)銅表面的BRDF,分析了溫度對(duì)測(cè)量結(jié)果的影響。結(jié)果表明:溫度對(duì)銅表面BRDF的測(cè)量結(jié)果有明顯的影響。然而,他們測(cè)量中采用632.8 nm的He-Ne激光器作為光源,主要研究單波長(zhǎng)下溫度對(duì)測(cè)量結(jié)果的影響,不能夠體現(xiàn)銅表面的光譜反射特性。

        本文的主要工作是基于自行研制的BRDF測(cè)量裝置測(cè)量了可見光波段紫銅表面的面內(nèi)光譜BRDF,并對(duì)影響測(cè)量結(jié)果的因素做了分析。論文的結(jié)構(gòu)如下:第二部分主要介紹本文所采用的BRDF測(cè)量方法與測(cè)量裝置,第三部分介紹樣品加工、測(cè)試和實(shí)驗(yàn)過程;第四部分給出測(cè)量結(jié)果與討論;第五部分給出基于Gauss和Lorentz分布的擬合結(jié)果與分析;最后對(duì)全文工作進(jìn)行了總結(jié)。

        2 測(cè)量原理與裝置

        2.1 BRDF測(cè)量原理

        BRDF的定義最早由Nicodemus給出[1],其將雙向反射分布函數(shù)fr定義為光輻射的反射輻亮度dLr與入射輻照度dEi的比值,即:

        fr(θi,φi,θr,φr,λ)=

        (1)

        式中,θ和φ分別表示天頂角和方位角,下標(biāo)i和r分別表示入射和反射,λ為波長(zhǎng)。BRDF的物理含義是:來自入射方向的表面輻照度的微增量與其所引起的反射方向上反射輻射亮度增量之間的比值。BRDF定義的幾何關(guān)系如圖1所示。

        圖1 雙向反射分布函數(shù)幾何關(guān)系圖 Fig.1 Geometric relationship diagram of incident and reflected beams for BRDF

        BRDF的測(cè)量方法有絕對(duì)測(cè)量和相對(duì)測(cè)量?jī)煞N。相對(duì)測(cè)量方法可以減小系統(tǒng)誤差,對(duì)雜散光也能起到較好的抑制作用。與絕對(duì)測(cè)量相比,相對(duì)測(cè)量對(duì)設(shè)備精度的要求較低,能夠節(jié)約測(cè)量成本。但是,相對(duì)測(cè)量方法的缺點(diǎn)也比較明顯,即必須依賴標(biāo)準(zhǔn)參考試樣。首先,目前尚未見到在所有波長(zhǎng)和所有環(huán)境中均適用的標(biāo)準(zhǔn)參考試樣。其次,并不是所有研究單位都具備精密測(cè)量標(biāo)準(zhǔn)參考試樣的條件,而計(jì)量部門提供的往往只是典型的入射和接收幾何條件下的數(shù)據(jù),與實(shí)際測(cè)量中的幾何條件不一定一致。因此,本文采用絕對(duì)測(cè)量方法測(cè)量紫銅表面的BRDF。

        本文采用的BRDF絕對(duì)測(cè)量方法在文獻(xiàn)[18]中有詳細(xì)介紹,這里僅給出測(cè)量表達(dá)式:

        fr(θi,φi,θr,φr,λ)=

        (2)

        對(duì)于均勻性、穩(wěn)定性和準(zhǔn)直性較好的光源,Lo(0,0,0,0,λ)是與波長(zhǎng)有關(guān)的函數(shù),完全由光源決定,可以在安放樣品前對(duì)光源進(jìn)行測(cè)量得到。對(duì)于固定的測(cè)量裝置,探測(cè)立體角Ωr是常數(shù),可以根據(jù)探測(cè)器的探測(cè)孔徑和探測(cè)距離計(jì)算得到;角度因子cosθi由入射天頂角得到。根據(jù)測(cè)量原理知,只需測(cè)量反射輻亮度Lr(θi,φi,θr,φr,λ)即可得到待測(cè)目標(biāo)的BRDF。

        2.2 測(cè)量裝置

        本文的測(cè)量所用裝置為一套自主設(shè)計(jì)的典型傳統(tǒng)BRDF測(cè)量裝置,其通過樣品、探測(cè)器和光源的配合轉(zhuǎn)動(dòng)實(shí)現(xiàn)不同入射和反射角度下BRDF的測(cè)量。測(cè)量裝置采用氙燈(GLORIA-X500A)作為光源,能夠提供入射光的波長(zhǎng)范圍為300~2 500 nm;采用光柵光譜儀作為探測(cè)系統(tǒng),可探測(cè)波長(zhǎng)范圍為250~2 500nm,光譜分辨率為0.4 nm。轉(zhuǎn)角裝置主要包括兩個(gè)轉(zhuǎn)臺(tái)和兩個(gè)轉(zhuǎn)臂,通過它們之間的配合轉(zhuǎn)動(dòng)可實(shí)現(xiàn)不同角度下BRDF的測(cè)量。測(cè)量裝置的設(shè)計(jì)圖如圖2所示,轉(zhuǎn)臺(tái)和轉(zhuǎn)臂的轉(zhuǎn)動(dòng)范圍見表1。

        表1 轉(zhuǎn)臂和轉(zhuǎn)臺(tái)的轉(zhuǎn)動(dòng)范圍

        本文實(shí)驗(yàn)系統(tǒng)的入射光路和測(cè)量光路如圖3(彩圖見期刊電子版)所示,分別用藍(lán)色和紅色箭頭表示。光路中M0~ M6是同一型號(hào)的保護(hù)銀膜平面反射鏡。虛線框中所示是本文所設(shè)計(jì)的探測(cè)聚焦附件,主要包括3個(gè)高透過會(huì)聚透鏡和一個(gè)通光孔徑為2 mm的光闌,這樣設(shè)計(jì)的主要目的是為了減小探測(cè)系統(tǒng)在樣品表面的有效探測(cè)面積。理論上講,安裝附件后,探測(cè)器只能夠探測(cè)到來自于透鏡L2焦點(diǎn)位置的反射光,能夠有效消除反射天頂角對(duì)有效探測(cè)面積的影響,這是絕對(duì)測(cè)量方法得以實(shí)現(xiàn)的基礎(chǔ)。

        3 樣品制備與測(cè)量

        本文準(zhǔn)備了4個(gè)半徑為25 mm、厚度為2 mm的圓片狀紫銅樣品。采用金相磨拋機(jī)對(duì)樣品進(jìn)行加工以得到不同的表面粗糙度。分別采用400目、600目、800目和1000目的碳化硅砂紙對(duì)樣品進(jìn)行打磨,依次標(biāo)記為1#、2#、3#和4#樣品。分別采用丙酮和無水乙醇對(duì)砂紙打磨后的樣品進(jìn)行清洗,之后,將樣品放入超聲波清洗機(jī)進(jìn)行進(jìn)一步的清洗,以確保樣品表面無污垢殘留。使用粗糙度儀TR220(時(shí)代公司)對(duì)處理完成后的樣品表面的粗糙度進(jìn)行測(cè)量,測(cè)量結(jié)果如表2所示。其中Ra、Rq和Rz分別為算術(shù)平均粗糙度、均方根粗糙度、輪廓最大高度。

        表2 4個(gè)樣品表面的粗糙度

        本文所采用的測(cè)量裝置是典型的傳統(tǒng)BRDF測(cè)量裝置,需要依靠轉(zhuǎn)角裝置實(shí)現(xiàn)不同角度下BRDF的測(cè)量。這導(dǎo)致測(cè)量過程十分耗時(shí),很難實(shí)現(xiàn)所有角度下目標(biāo)表面光譜BRDF的測(cè)量。本文的測(cè)量中,入射天頂角分別設(shè)置為30°、45°和60°,入射和反射方位角分別設(shè)定為0°和180°,反射天頂角的范圍為-80°~80°,測(cè)量間隔在鏡面反射方向附近設(shè)定為1°,遠(yuǎn)離鏡面反射方向的區(qū)域設(shè)定為2°;波長(zhǎng)范圍設(shè)置為380~780 nm,間隔為5 nm。如前文所述,本文采用絕對(duì)測(cè)量方法實(shí)現(xiàn)BRDF的測(cè)量,測(cè)量過程主要包括光源光譜亮度的測(cè)量和樣品表面反射亮度的測(cè)量。首先,取下樣品臺(tái),將入射天頂角、反射天頂角、入射方位角、反射方位角分別設(shè)置為90°、90°、0°、180°,測(cè)量氙燈的光譜亮度。然后,安裝樣品臺(tái)并放置樣品,按照設(shè)定的目標(biāo)依次測(cè)量4個(gè)樣品表面的面內(nèi)光譜BRDF。在分析測(cè)量結(jié)果之前需要指出的是:對(duì)于粗糙的表面而言,BRDF在一般情況下不可能等于零,然而在本文的測(cè)量中出現(xiàn)了等于零的情況,這主要?dú)w因于探測(cè)器的靈敏度有限。在本文的結(jié)果分析中,不再顯示BRDF等于零的部分。

        4 測(cè)量結(jié)果與討論

        為了研究紫銅表面的BRDF隨波長(zhǎng)和粗糙度的變化情況,本文測(cè)量了入射天頂角為30°時(shí)的面內(nèi)光譜BRDF,1#~4#樣品的測(cè)量結(jié)果在圖4(彩圖見期刊電子版)中給出。從圖4能夠看出,在入射波長(zhǎng)的情況下,BRDF隨著反射天頂角的增加,先增大后減小,在鏡面反射方向達(dá)到最大值,這一變化趨勢(shì)幾乎不受波長(zhǎng)變化的影響。另外,從圖中能夠看出,紫銅在380~550 nm波段存在較強(qiáng)的吸收,波長(zhǎng)大于550 nm時(shí),隨著波長(zhǎng)的增加吸收迅速減弱,并且這一變化趨勢(shì)不受反射天頂角,即觀察方向的影響,這與紫銅表面所呈現(xiàn)出的顏色是一致的。在可見光波段,目標(biāo)的光譜吸收和反射特性很大程度上取決于目標(biāo)表面顏色,紫銅的顏色是色素色不隨觀察方向的改變而發(fā)生變化。測(cè)量結(jié)果與紫銅所呈現(xiàn)出的視覺特性是一致的。對(duì)比圖4中的4幅子圖能夠發(fā)現(xiàn),不同粗糙度下紫銅表面的BRDF隨波長(zhǎng)和反射天頂角的變化趨勢(shì)比較相似。在給定波長(zhǎng)的情況下,鏡面反射方向的BRDF隨著粗糙度的減小而逐漸增大,這是因?yàn)殡S著表面粗糙度的減小,樣品表面越來越接近于光學(xué)光滑表面,導(dǎo)致鏡面反射方向的反射增強(qiáng)。根據(jù)圖4的測(cè)量結(jié)果能夠得到如下規(guī)律:紫銅表面的BRDF隨反射天頂角的變化不受波長(zhǎng)的影響,隨波長(zhǎng)的變化趨勢(shì)不受反射天頂角的影響。如果這個(gè)規(guī)律正確,那么可以將紫銅的空間光譜反射特性分為光譜反射特性和空間反射特性,即在單波長(zhǎng)下研究紫銅表面的空間反射特性,在單個(gè)反射方向下研究其光譜反射特性。這樣可以大幅度減小BRDF的測(cè)量時(shí)間和測(cè)量成本。

        圖4 θi =30°時(shí),樣品的面內(nèi)光譜BRDF,(a)~(d)依次為1#~4#樣品的測(cè)量結(jié)果 Fig.4 In-plane spectral BRDFs of samples at θi=30°:(a)1# sample, (b)2# sample, (c)3# sample and (d)4# sample

        為了驗(yàn)證上文所述規(guī)律的正確性,文中測(cè)量了入射天頂角為45°和60°時(shí)紫銅表面的光譜BRDF,測(cè)量結(jié)果分別在圖5(彩圖見期刊電子版)和圖6(彩圖見期刊電子版)中給出。從圖中可以看出,BRDF隨反射天頂角的變化趨勢(shì)不受波長(zhǎng)影響,隨波長(zhǎng)的變化趨勢(shì)不受反射天頂角的影響。這很好地驗(yàn)證了上述結(jié)規(guī)律的正確性,也說明將紫銅表面的空間光譜反射特性分為單入射波長(zhǎng)下的空間反射特性和單反射方向下的光譜反射特性是可行的。

        圖5 θi =45°時(shí),樣品的面內(nèi)光譜BRDF,(a)~(d)依次為1#~4#樣品的測(cè)量結(jié)果 Fig.5 In-plane spectral BRDFs of samples at θi =45°: (a) 1# sample, (b) 2# sample , (c) 3# sample and (d) 4# sample

        圖6 θi=60°時(shí),樣品的面內(nèi)光譜BRDF,(a)~(d)依次為1#~4#樣品的測(cè)量結(jié)果 Fig.6 In-plane spectral BRDFs of samples at θi=60°:(a)1# sample, (b)2# sample , (c)3# sample and (d)4# sample

        圖7 650 nm下,BRDF隨反射天頂角的變化曲線,(a)~(d)依次為1#~4#樣品的測(cè)量結(jié)果 Fig.7 Curves of BRDF varying with reflectance zenith at 650 nm:(a) 1# sample, (b) 2# sample , (c) 3# sample and (d) 4# sample

        為了更清晰地展示紫銅表面BRDF隨反射天頂角的變化,本文提取了650 nm下,入射天頂角為30°、45°和60°時(shí)4個(gè)樣品的測(cè)量結(jié)果,分別在圖7(a)、7(b)、7(c)和7(d)中給出。從圖7(a)中能夠看出, BRDF隨著反射天頂角的增加呈現(xiàn)出先增大后減小的趨勢(shì),在鏡面反射反向附近達(dá)到最大值。θi=45°和θi=60°時(shí)BRDF隨反射天頂角的變化趨勢(shì)與θi=30°時(shí)的相同,這說明入射天頂角對(duì)BRDF隨反射天頂角的變化趨勢(shì)沒有明顯的影響。另外,從圖7(a)中能夠看出,隨著反射天頂角的增加,BRDF在鏡面反射方向附近的值逐漸增加。從圖7(b)、7(c)和7(d)中能夠看出,2#~4#樣品的BRDF隨反射天頂角的變化趨勢(shì)與1#樣品相同,這說明表面粗糙度對(duì)BRDF隨反射天頂角的變化趨勢(shì)沒有明顯的影響。

        為了更清楚地展示粗糙度對(duì)測(cè)量結(jié)果的影響,圖8給出了θi=45°時(shí),4個(gè)樣品在λ=650 nm時(shí)的測(cè)量結(jié)果。

        圖8 650 nm下,θi=45°時(shí)BRDF隨反射天頂角的變化曲線 Fig.8 Curves of BRDF varying with refletance zenith at 650 nm, θi=45°

        很明顯,樣品的表面粗糙度對(duì)BRDF隨反射天頂角變化的趨勢(shì)沒有明顯的影響,4個(gè)樣品的BRDF隨反射天頂角的變化曲線均呈“鐘形”,曲線的最大值在鏡面反射方向。從圖8能夠看出,粗糙度對(duì)BRDF的大小有很大的影響,尤其在鏡面反射方向附近。以鏡面反射方向?yàn)槔?,隨著Ra由0.219 μm減小為0.072 μm,BRDF值由1.81 sr-1增大為16.80 sr-1。另外,隨著粗糙度的減小,BRDF隨θr變化的曲線變得越來越“高”、越來越“瘦”。這說明,隨著粗糙度的減小,樣品表面在鏡面反射方向附近的反射逐漸增強(qiáng),偏離鏡面反射方向的強(qiáng)度逐漸減小。

        從圖7的測(cè)量結(jié)果能夠看出,4個(gè)樣品的BRDF尤其是在鏡面反射方向的值隨入射天頂角的增大而增大。產(chǎn)生這種影響的原因可以歸結(jié)為兩個(gè)方面:(1)隨著θi的增大,樣品在鏡面反射方向的反射增強(qiáng);(2)角度因子的影響,測(cè)量原理式中有cosθi項(xiàng),隨著θi的增大,角度因子逐漸減小,從而導(dǎo)致BRDF增大,在θi較大時(shí)更加明顯。測(cè)量原理中的探測(cè)立體角Ωr由探測(cè)器的入光孔徑和探測(cè)器到樣品之間的距離決定,測(cè)量裝置設(shè)計(jì)完成后,探測(cè)立體角為常數(shù)。分析測(cè)量原理能夠發(fā)現(xiàn),fr· cosθi能夠很好地描述目標(biāo)表面的反射強(qiáng)度。為了進(jìn)一步分析反射強(qiáng)度和角度因子對(duì)BRDF的影響,本文給出了不同入射天頂角下鏡面反射方向的BRDF與fr· cosθi隨角度的變化曲線,如圖9所示。測(cè)量波長(zhǎng)為650 nm,θi的設(shè)置范圍為20°≤θi≤70°,測(cè)量間隔為5°。

        圖9 650 nm時(shí),樣品在鏡面反射方向的BRDF(a)和fr·cosθi(b) Fig.9 BRDF(a) and fr·cosθi(b) at specular direction for the samples at 650 nm

        由圖9(a)可以看出,在給定天頂角的情況下,BRDF隨著粗糙度的增加而減小。對(duì)于每一個(gè)樣品,當(dāng)θi≤45°時(shí),BRDF隨著θi的增大而緩慢增加;當(dāng)θi>45°時(shí),BRDF隨著θi的增大而快速增加。為了消除角度因子cosθi對(duì)測(cè)量結(jié)果的影響,圖9(b)給出了不同入射天頂角下鏡面反射方向的fr·cosθi。從圖中可以看出,fr·cosθi隨θi的變化趨勢(shì)與BRDF的變化趨勢(shì)相同。圖9的測(cè)量結(jié)果說明,隨著入射天頂角的增大,鏡面反射方向的反射強(qiáng)度與角度因子對(duì)BRDF均有積極作用。

        5 實(shí)驗(yàn)數(shù)據(jù)擬合

        這部分的主要工作是對(duì)測(cè)量結(jié)果進(jìn)行擬合,以比較不同分布函數(shù)對(duì)測(cè)量值的擬合結(jié)果?,F(xiàn)階段常用的BRDF模型,如Phong模型、Cook-Torrance模型和Ward模型等均采用Gauss分布函數(shù)作為基礎(chǔ)函數(shù)。這是因?yàn)榇蠖鄶?shù)BRDF模型將空間反射特性近似看做理想朗伯體與鏡面反射的疊加,而Gauss分布比較適用于這種近似。然而,目標(biāo)表面在遠(yuǎn)離鏡面反射的方向一般都不滿足理想朗伯體近似,即使是在BRDF相對(duì)測(cè)量方法中作為參照的的白板也不是理想朗伯體。尋找更合適的基礎(chǔ)函數(shù)對(duì)BRDF經(jīng)驗(yàn)或者半經(jīng)驗(yàn)?zāi)P偷慕⒂兄匾饬x。文中選取多種不同的分布函數(shù)對(duì)測(cè)量結(jié)果進(jìn)行擬合,發(fā)現(xiàn)采用Lorentz分布得到的擬合結(jié)果與測(cè)量值最接近,文中給出了入射天頂角為45°時(shí),采用Lorentz分布的擬合結(jié)果與采用Gauss分布的擬合結(jié)果進(jìn)行了對(duì)比,如圖10所示。

        圖10 θi=45°,λ=650 nm 4個(gè)樣品的測(cè)量與擬合結(jié)果 Fig.10 Measuring and fitting results of the four samples at θi =45° for λ=650 nm

        從圖10能夠看出,采用Lorentz分布的擬合結(jié)果比采用Gauss分布得到的結(jié)果更好。對(duì)比這兩種擬合結(jié)果能夠看出,在鏡面反射方向附近,Lorentz和Gauss分布的擬合結(jié)果與測(cè)量值符合的較好。在遠(yuǎn)離鏡面反射方向的角度,采用Lorentz分布的擬合結(jié)果優(yōu)于Gauss分布的擬合結(jié)果。這說明采用Lorentz分布作為基礎(chǔ)函數(shù)建立BRDF經(jīng)驗(yàn)或者半經(jīng)驗(yàn)?zāi)P?,可能得到更好的結(jié)果。

        6 結(jié) 論

        本文通過自主設(shè)計(jì)的BRDF測(cè)量裝置測(cè)量了紫銅表面的BRDF,詳細(xì)分析了入射波長(zhǎng),入射天頂角,表面粗糙度對(duì)測(cè)量結(jié)果的影響。結(jié)果表明:(1)BRDF隨波長(zhǎng)的變化趨勢(shì)不受入射和反射角度的影響;(2)鏡面反射方向附近的BRDF值隨著入射天頂角的增大而逐漸增大,同時(shí)受角度因子和樣品表面反射強(qiáng)度的影響;(3)測(cè)量結(jié)果容易受樣品表面粗糙度的影響,隨著表面粗糙度的增加,BRDF隨反射天頂角的變化曲線呈現(xiàn)出越來越“矮”,越來越“胖”的趨勢(shì)。另外,本文采用不同的分布函數(shù)對(duì)測(cè)量結(jié)果進(jìn)行了擬合,并將得到的最優(yōu)擬合結(jié)果與采用Guass分布的擬合結(jié)果進(jìn)行了比較,結(jié)果表明:采用Lorentz擬合能夠得到更好的結(jié)果。

        猜你喜歡
        鏡面反射紫銅頂角
        一般三棱鏡最大頂角與折射率的關(guān)系
        光滑物體表面反射光偏振特征分析及反射光分離技術(shù)*
        紫銅紅灰蝶
        涼亭中的數(shù)學(xué)
        基于最短路徑的GNSS-R鏡面反射點(diǎn)算法
        頂角為100°的等腰三角形性質(zhì)的應(yīng)用
        低合金鋼管板與紫銅換熱管焊接工藝的探索
        焊接(2015年2期)2015-07-18 11:02:38
        紫銅零件小尺寸螺紋孔的攻絲分析
        磷及稀土鑭元素對(duì)紫銅焊絲熔敷金屬組織和性能的影響
        雙元計(jì)算法求體積本征值
        亚洲国产精品久久久久秋霞小说 | 九九久久精品无码专区| 亚洲中文字幕乱码| 日本韩国黄色三级三级| 大香蕉国产av一区二区三区| 国产精品无码久久综合| 国产mv在线天堂mv免费观看| 亚洲第一区二区快射影院| 中文字幕乱码一区在线观看| 国产乡下妇女做爰| 国内a∨免费播放| 国产成人精品免费视频大全| 国产在线观看女主播户外| 国产av久久久久精东av| 国产自偷亚洲精品页65页| 国产传媒在线视频| 中文字幕综合一区二区| 97久久精品人妻人人搡人人玩| 亚洲av无码片一区二区三区| 日韩不卡无码三区| 中文字幕中文字幕在线中二区| 日韩欧美成人免费观看| 国产久视频国内精品999| 冲田杏梨av天堂一区二区三区| 亚洲啪啪视频一区二区| 97夜夜澡人人爽人人喊中国片| 久久夜色精品国产亚洲噜噜| 熟女一区二区国产精品| 免费人成激情视频在线观看冫| 久久精品国产夜色| 久久精品一区二区三区不卡牛牛| 久久久精品人妻一区二区三区四区| 人妻少妇精品中文字幕av蜜桃| 中出高潮了中文字幕| 尤物国产一区二区三区在线观看| 亚洲国产成人片在线观看| 国产美女在线一区二区三区| 国产av大片久久中文字幕| 亚洲精品成人无限看| 亚洲日韩精品国产一区二区三区| 人妻少妇喷水意淫诱惑|