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

    基于偏振光的干涉測長法

    2021-06-01 08:51:22胡婷婷杜玉嫻白在橋
    物理實(shí)驗(yàn) 2021年5期
    關(guān)鍵詞:金屬桿光程干涉儀

    胡婷婷,杜玉嫻,白在橋

    (北京師范大學(xué) 物理學(xué)系,北京 100875)

    利用光的干涉測量微小位移或光程差是一類重要的大學(xué)物理實(shí)驗(yàn). 常見的實(shí)驗(yàn)包括利用邁克耳孫干涉儀測量空氣折射率、薄膜折射率或厚度、金屬彈性模量/線脹系數(shù)、材料的電/磁致伸長系數(shù)等[1-7]. 實(shí)驗(yàn)中通常是通過數(shù)條紋的數(shù)量來確定動鏡位移或其他因素導(dǎo)致的光程差. 數(shù)條紋雖然操作簡便,但也存在一些問題,比如:

    1)只適合位移緩慢變化的情況. 如果位移變化過快,人眼會看不清條紋的結(jié)構(gòu);

    2)只適合位移(光程差)單方向運(yùn)動的情況. 如果位移方向不固定,條紋既有涌出也有陷入,人工計(jì)數(shù)將變得十分困難;

    3)讀數(shù)限定為半波長的整數(shù)倍. 為保證測量精度,一般需要幾十甚至上百個(gè)條紋,不適合測量波長或亞波長范圍的位移量.

    前2個(gè)問題可以用(單向或可逆)自動條紋計(jì)數(shù)器解決[8-12]. 但為了解決第3個(gè)問題,必須采用比數(shù)條紋更精細(xì)的方法.

    1 實(shí)驗(yàn)原理

    光強(qiáng)分別為I1和I2的2束光的干涉光強(qiáng)為

    (1)

    其中相位差Δφ直接與位移相關(guān). 測出I1,I2和I,就可以根據(jù)(1)式算出cos Δφ. 不過在1個(gè)周期內(nèi)cos Δφ與Δφ是一對二的關(guān)系,除非已知位移的運(yùn)動規(guī)律(比如勻速運(yùn)動),一般情況下并不能由干涉光強(qiáng)推知位移的變化情況. 本文針對此問題,利用光的偏振自由度,無需對位移的變化方式做任何假設(shè),即可以測出干涉相位,從而可以穩(wěn)定地實(shí)現(xiàn)位移的連續(xù)測量.

    1.1 原理圖與基本思路

    光路如圖1所示,在傳統(tǒng)邁克耳孫干涉儀光路(虛線框)的基礎(chǔ)上做了改進(jìn):

    1)激光入射分束鏡(BS1)之前先經(jīng)過1塊偏振片(P)和1塊分束鏡(BS0);

    2)激光從干涉儀出射后經(jīng)過1塊偏振分束鏡(PBS);

    3)使用3個(gè)光電探測器分別(D0,Dp,Ds)測量從BS0分出的光強(qiáng)I0和從PBS出射的2路光強(qiáng)Ip0和Is0.

    偏振片的透振方向與水平方向大致成45°,這樣入射光中p分量和s分量比例固定,基本為1∶1. BS0和光探測器D0的作用是監(jiān)視入射光強(qiáng). 相對光強(qiáng)為

    (2)

    結(jié)果將不受激光器輸出狀態(tài)漂移的影響.

    圖1 測量原理圖

    由于p光和s光之間不會產(chǎn)生干涉,所以干涉儀中同時(shí)存在2套干涉光路,它們在空間上完全重合,直至PBS把它們分開,進(jìn)入對應(yīng)的探測器. 因此有

    (3)

    相位差來自干涉儀兩臂的光程差和光學(xué)元件導(dǎo)致的光程差,即

    (4)

    式中λ為激光波長,ΔL為干涉儀兩臂的光程差,Δφp0與Δφs0分別為光學(xué)元件對p光和s光引入的相位差. 由于分束鏡鍍膜對光的反射與透射都依賴于光的偏振狀態(tài),一般Δφp0≠Δφs0. 因此,當(dāng)ΔL改變時(shí),Ip和Is的變化有固定的相位差δ0≡Δφp0-Δφs0. 只要δ0不等于π的整數(shù)倍,(Ip,Is)就會在固定的橢圓上移動. 記錄測量點(diǎn)在橢圓上的運(yùn)動軌跡,就可推知ΔL如何隨時(shí)間變化,這就是本文測量方法的基本思路. 雖然原則上只需要δ0≠kπ,但如果δ0過于接近kπ,橢圓就變得很扁長,測量結(jié)果就會因噪聲的影響變得不穩(wěn)定. 在本文的實(shí)驗(yàn)裝置中δ0大致為45°,為降低噪聲的影響. 如果條件允許,可在干涉光路的1條臂中插入可調(diào)相位延遲片,即可對δ0進(jìn)行調(diào)節(jié).

    1.2 計(jì)算方法

    根據(jù)(3)式和(4)式,不妨設(shè)

    (5)

    經(jīng)過仿射(平移加線性)變換:

    (6)

    (7)

    (Ip′,Is′)在單位圓上,而且輻角隨光程差均勻變化,光程差改變λ,(Ip′,Is′)在單位圓上運(yùn)動1周. 如果采樣率足夠高或者光程差變化足夠慢,可以保證前后2個(gè)采樣點(diǎn)之間的光程差改變量小于λ/2,那么就可以通過連續(xù)跟蹤,得到光程差隨時(shí)間的變化曲線.

    ΔL很可能不隨時(shí)間勻速變化,因此一般不能通過正弦曲線擬合得到Ap,Bp,As,Bs與δ0. 下面給出根據(jù)數(shù)據(jù)點(diǎn)的幾何性質(zhì)(即它們位于個(gè)橢圓上),來擬合這5個(gè)參量的方法. 平面上橢圓的一般方程為

    (8)

    其中a,b,c,d,e,f為常數(shù),而且矩陣

    (9)

    是正定的. 方程(8)可寫成XTP=0,其中向量X和P定義為

    (10)

    設(shè)測量點(diǎn)的集合為{(Ipi,Isi)|i=1,2,…,N}. 由于存在測量誤差,一般并不存在1個(gè)向量P與所有的Xi都正交. 對所有實(shí)測數(shù)據(jù)點(diǎn)求和,得到對稱矩陣

    (11)

    然后用M的最小特征值對應(yīng)的特征向量作為對P的估計(jì). 由于特征向量可以相差常數(shù)因子,這里要求a>0. 由于D是正定的,同時(shí)一定有b>0和ab-c2>0.

    a,b,c,d,e,f確定后,就可以計(jì)算Ap,Bp,As,Bs與δ0了. 首先,(Ap,As)為橢圓的中心,可表示為

    (12)

    (Bp,Bs)可根據(jù)下式計(jì)算

    (13)

    其中

    (14)

    最后,相位差δ0可根據(jù)下式計(jì)算

    (15)

    由于三角函數(shù)的特性,可以得到4種不同的仿射變換(分別對應(yīng)δ0→±δ0,2π±δ0). 實(shí)驗(yàn)中固定選擇1種即可.

    利用上述方法得到仿射變換的參量,將數(shù)據(jù)點(diǎn)映射到單位圓上. 計(jì)算輻角(即相位),并根據(jù)前一時(shí)刻的值,通過適當(dāng)平移2kπ,把輻角連續(xù)地映射在實(shí)軸上. 最后乘以λ/2π(或者λ/4π)得到光程差ΔL(或反射鏡相對位移S)的變化曲線.

    2 裝置與方法

    實(shí)驗(yàn)裝置以Thorlabs公司的邁克耳孫干涉儀教學(xué)套裝為基礎(chǔ)搭建,實(shí)物如圖2所示. 光源采用氦氖激光器(λ=632.8 nm,功率1.5 mW). 干涉儀的反射鏡M1固定在二維調(diào)整架上,反射鏡M2固定在待測金屬桿的前端,金屬桿的尾部通過連接桿用螺絲固定在光學(xué)面包板上. 金屬桿外面包裹加熱箔,可通過電流改變金屬桿的溫度. 金屬桿的對稱軸上鉆有1個(gè)深孔,可以將溫度傳感器(K型溫差電偶)插入,測量金屬桿的溫度. 偏振分束鏡使用Thorlabs PBS201. 3個(gè)光電探測器均為硅光電二極管(Thorlabs DET36A2). 探測器輸出至數(shù)據(jù)采集卡(NI USB-6212). 由于D0探測的光信號變化緩慢,所以采用了較大的耦合電阻(10 kΩ),為了提高響應(yīng)速度,Ds和Dp使用了較小的耦合電阻(1 kΩ). 因此Is和Ip要比I0小1個(gè)數(shù)量級. 所有元件,除了激光器和偏振片,都固定在30 cm×45 cm光學(xué)面包板上,整體放在光學(xué)平臺上.

    裝置的調(diào)節(jié)方法與通常的邁克耳孫干涉儀相同. 為觀察調(diào)節(jié)效果,需要在激光器輸出位置插入擴(kuò)束鏡,并用白屏觀察干涉圖樣. 要求兩臂的光程盡量相等,使干涉圖樣的中心斑點(diǎn)遠(yuǎn)大于探測器的有效接收區(qū)域(2.2 mm×2.2 mm). 測量時(shí)取走擴(kuò)束鏡,讓光直接打在探測器上.

    1.He-Ne激光器 2.偏振片 3.分束鏡(非偏振) 4.反射鏡M1 5.反射鏡M2 6.PBS 7.探測器圖2 裝置照片

    利用此裝置可完成以下2個(gè)實(shí)驗(yàn):

    a.測量系統(tǒng)的沖擊響應(yīng). 輕微敲擊光學(xué)面包板,干涉條紋出現(xiàn)晃動,然后逐步平穩(wěn). 編寫LabVIEW程序,當(dāng)監(jiān)測到振動(干涉光強(qiáng)發(fā)生劇烈變化)后,記錄1 s的數(shù)據(jù)(采樣率105s-1). 用全部的測量數(shù)據(jù)擬合橢圓參量,確定仿射變換,對光強(qiáng)進(jìn)行變換,計(jì)算輻角,得到光程差的變化曲線.

    b.測量金屬的線脹系數(shù). 加熱金屬桿使其長度發(fā)生改變,M2發(fā)生位移,引起干涉條紋移動. 實(shí)驗(yàn)所用數(shù)字溫度計(jì)可將溫度按1 mV/℃的格式轉(zhuǎn)換成電壓輸出. 利用數(shù)字萬用表(Rigol DM3058E)測量輸出電壓,并且編程通過VISA自動讀取測量結(jié)果. 直流電源(Rigol DP1116A)也可以通過VISA控制. 程序控制加熱溫度超過TH=36 ℃就斷電,冷卻到溫度低于TL=31 ℃就通電加熱(加熱電流0.3 A). 可盡量減小人為操作失誤和讀數(shù)誤差,便于驗(yàn)證測量結(jié)果的可重復(fù)性. 金屬桿升溫和降溫的時(shí)間都較長(大約5 min),光程差變化較平緩. 實(shí)驗(yàn)需要接近實(shí)時(shí)地顯示金屬桿的長度變化,實(shí)驗(yàn)中每秒采集1 000組光強(qiáng)數(shù)據(jù)(采樣率103s-1),這些數(shù)據(jù)對應(yīng)的光程差變化不大,不能形成完整的橢圓. 因此保存30 s的光強(qiáng)數(shù)據(jù)(每次采集到新數(shù)據(jù)后都需要更新),用所有數(shù)據(jù)擬合橢圓,確定仿射變換參量,并且把當(dāng)前1 s的數(shù)據(jù)轉(zhuǎn)換成桿的伸長量.

    3 結(jié)果與分析

    3.1 沖擊響應(yīng)測量

    圖3為1次采集到的光強(qiáng)波形. 可以看出光強(qiáng)變化相當(dāng)劇烈. 在前面一段時(shí)間(大致0~0.6 s),光強(qiáng)變化的最大值與最小值基本固定,對應(yīng)光程差的變化幅度大于λ,有超過1個(gè)條紋的吞吐. 之后光強(qiáng)的極值隨時(shí)間改變,對應(yīng)光程差的變化幅度小于λ,變化小于1個(gè)完整的條紋.

    圖3 不同偏振相對干涉光強(qiáng)隨時(shí)間變化曲線

    圖4給出了其中一部分的放大. 注意圖中標(biāo)出的9個(gè)光強(qiáng)的極大值,其中1,3,5,7,9號Ip和Is看起來是同相位的,而2,4,6,8號Ip和Is有明顯的相位差. 前一組極大值對應(yīng)光程差的極值(即反射鏡振動位移到極值處),這對p光和s光沒有差別. 而后一組極值對應(yīng)光程差達(dá)到為λ的整數(shù)倍,與偏振方向有關(guān),所以Ip和Is有明顯的相位差.

    圖4 圖3波形的細(xì)節(jié)放大

    將(Ip,Is)數(shù)據(jù)按xy方式繪圖,得到圖5. 雖然存在較大的噪聲,但數(shù)據(jù)點(diǎn)形成了很明顯的橢圓(圖中紅色為擬合橢圓). 注意Ip和Is的數(shù)值越大,噪聲也越大. 用數(shù)據(jù)擬合橢圓,并用所得的仿射變換作用于原始數(shù)據(jù),結(jié)果見圖6. 可以看出(Ip′,Is′)都在單位圓(紅色曲線)附近. 雖然結(jié)果同樣存在較大的噪聲,但數(shù)據(jù)距離原點(diǎn)都足夠遠(yuǎn),輻角的計(jì)算是可靠的.

    將輻角折合成為反射鏡的相對位移,得到振動波形,結(jié)果見圖7. 圖中可以看到明顯的振動衰減過程. 放大對應(yīng)圖4的部分(見圖8),可見清晰的正弦振動. 注意圖中的各位移極值點(diǎn),它們與圖4中的1,3,5,7,9號極值對應(yīng). 而且噪聲水平明顯小于圖4中的光強(qiáng)信號. 噪聲減小的原因是Ip和Is的噪聲存在顯著的正向關(guān)聯(lián),因此噪聲主要影響數(shù)據(jù)點(diǎn)到原點(diǎn)的距離,而對輻角的影響相對較小.

    圖5 在(Ip,Is)平面顯示數(shù)據(jù)

    圖6 經(jīng)過仿射變換后的數(shù)據(jù)

    圖7 系統(tǒng)振動衰減曲線

    圖7中的衰減曲線存在明顯的拍,說明系統(tǒng)存在2個(gè)相差不大的固有頻率. 對衰減曲線做傅里葉變換,結(jié)果見圖9. 從圖中可得這2個(gè)固有頻率分別為206 Hz和224 Hz. 此外,系統(tǒng)還存在260 Hz,398 Hz,484 Hz等固有頻率.

    圖8 衰減曲線細(xì)節(jié)放大

    圖9 衰減曲線的頻譜

    進(jìn)行多次測量,每次敲擊光學(xué)面包板的位置或力度不相同. 表1為21次測量所得各個(gè)參量擬合值的平均值值與標(biāo)準(zhǔn)差. 可以發(fā)現(xiàn)Ap,Bp,As和Bs擬合值的相對標(biāo)準(zhǔn)偏差都小于0.1%,δ0的標(biāo)準(zhǔn)偏差只有0.002 drad≈0.12°,這說明實(shí)驗(yàn)光路、探測系統(tǒng)和數(shù)據(jù)處理算法的穩(wěn)定性很好.

    表1 擬合參量統(tǒng)計(jì)

    3.2 金屬線脹系數(shù)測量

    圖10顯示了連續(xù)4 000 s的測量數(shù)據(jù),圖10(a)每個(gè)點(diǎn)是1 s內(nèi)1 000次桿伸長量的平均值,圖10(b)是從數(shù)字萬用表讀出的溫度值. 圖中可見有7次完整的升溫和冷卻過程. 伸長量δL和溫度T的變化形式都類似電容充放電. 這是因?yàn)闇囟仍礁?,散熱越快,?dǎo)致加熱段升溫和冷卻段降溫都越來越慢.

    (a)

    (b)圖10 金屬桿伸長量δL與溫度T隨時(shí)間變化曲線

    在(T,δL)平面中畫出測量數(shù)據(jù),結(jié)果見圖11. 圖中可以看出在升溫段和冷卻段,δL都與T呈很好的直線關(guān)系. 但這兩段明顯不重合,T的變化落后于δL,表現(xiàn)為滯回性. 可理解為:δL反映了桿平均溫度的改變,溫度探頭測量的是金屬桿內(nèi)部的溫度,但加熱和散熱都發(fā)生在金屬桿的表面,因此T的變化要滯后于δL.

    圖11 在(T,δL)平面顯示測量數(shù)據(jù)

    取T在31.5~35.5 ℃之間連續(xù)的升溫或冷卻階段的數(shù)據(jù). 這14段數(shù)據(jù)的線性相關(guān)性系數(shù)都大于0.999 93. 用最小二乘法計(jì)算δL對T變化的斜率κ=(1.880 5±0.007 8) μm/℃,結(jié)果見表2. 可看出升溫段(奇數(shù)次)的斜率總體上略小于冷卻段(偶數(shù)次)的斜率. 原因可能是升溫段時(shí)間略短于冷卻段,因此受滯回效應(yīng)的影響會更大. 滯回效應(yīng)總會導(dǎo)致測量斜率偏小,所以升溫段的斜率比冷卻段要小. 對擬合斜率取平均,并取等效桿長L=79 mm,計(jì)算鋁的線脹系數(shù),結(jié)果在參考值(2.1×10-5~2.4×10-5℃-1)范圍之內(nèi).

    表2 δL對T的斜率κ擬合值

    上面測量結(jié)果最大的誤差來源于等效桿長L的測量. 因?yàn)榻饘贄U是通過連接桿用螺絲緊固在光學(xué)面包板上,不易確定后端的位置. 79 mm是從緊固螺絲與金屬桿的接觸點(diǎn)算起的,本身不便測量,誤差也比較大. 不過這個(gè)問題與本文主旨無關(guān). 重要的是采取本文的方法,可以對微小位移進(jìn)行連續(xù)測量. 在持續(xù)超過1 h的時(shí)間內(nèi),結(jié)果具有很好的可重復(fù)性. 由于位移分辨率遠(yuǎn)大于數(shù)條紋的方式,這種方法不僅提高了探測靈敏度,同時(shí)也清晰地展示了一些有趣的物理過程或現(xiàn)象,比如桿的散熱和溫度非均勻分布.

    4 討 論

    干涉法測量微小位移(光程差)本質(zhì)上是測量2維向量(即干涉復(fù)振幅)的角度. 干涉光強(qiáng)對應(yīng)的是該向量在1個(gè)方向上的投影. 有2個(gè)方向的投影就可以唯一地恢復(fù)出向量的角度. 可逆條紋計(jì)數(shù)器就是利用了此思想[10-12]:用光纖探頭在不同空間位置測量干涉光強(qiáng),得到有一定相位差(最好是π/2)的2個(gè)信號. 根據(jù)2個(gè)信號達(dá)到極值的時(shí)序可以判斷條紋是吞進(jìn)或吐出. 這個(gè)方法的缺點(diǎn)是需要使用擴(kuò)束鏡,而且2個(gè)光強(qiáng)取樣點(diǎn)的相位很容易受到擾動,導(dǎo)致信號弱且重復(fù)性不好. 本文使用2個(gè)不同偏振狀態(tài)的干涉光強(qiáng),它們在空間完全重合,測量時(shí)也不需要擴(kuò)束. 因此數(shù)據(jù)信噪比高,結(jié)果可重復(fù)性好. 經(jīng)過適當(dāng)?shù)臄?shù)學(xué)處理,可以實(shí)現(xiàn)位移(光程差)的連續(xù)測量.

    根據(jù)圖8中峰-峰值小于1個(gè)波長的正弦振動曲線估計(jì),單次位移測量的誤差應(yīng)該比波長小2個(gè)數(shù)量級. 為了估計(jì)位移測量的分辨率,保持干涉儀不動,連續(xù)測量反射鏡的相對位移(固定變換參量). 圖12中展示了典型的結(jié)果. 可以看見測量結(jié)果存在持續(xù)的漲落. 漲落的低頻分量來自不可控制的環(huán)境擾動導(dǎo)致的光程差漂移. 這部分反映的是反射鏡實(shí)際的相對位移. 漲落的高頻分量(即短時(shí)間的漲落)是由探測器和采集卡的噪聲造成的,也與數(shù)據(jù)處理算法有關(guān). 即便光程差嚴(yán)格固定不動,這部分高頻噪聲仍然存在,因此它決定了實(shí)際的位移測量分辨率. 把測量結(jié)果通過截止頻率為20 Hz低通數(shù)字濾波器,得到其低頻分量(圖12中紅色曲線). 用原始信號減去低頻分量,得到高頻分量. 圖13為高頻分量的統(tǒng)計(jì)直方圖. 可以看出形式很接近高斯分布,而且絕大部分值在±0.004 μm之間. 事實(shí)上,實(shí)驗(yàn)發(fā)現(xiàn)高頻分量的標(biāo)準(zhǔn)差在0.001~0.005 μm之間,與光強(qiáng)的噪聲水平正相關(guān). 如果用3倍的標(biāo)準(zhǔn)差作為位移測量分辨率,則結(jié)果在0.003~0.015 μm之間,比波長小2個(gè)數(shù)量級. 此處估算的是單次測量的分辨率. 如果位移變化緩慢,提高采樣率,用多次采樣取平均值,可以進(jìn)一步提高位移測量的分辨率(不過這樣會降低時(shí)間分辨率).

    圖12 一段反射鏡相對位移測量波形(黑色)及其經(jīng)過低通濾波后的波形(紅色)

    利用本文方法,還進(jìn)行了音叉共振曲線測量、剛體擺的隨機(jī)振動測量等實(shí)驗(yàn). 該方法結(jié)合了光的干涉和偏振,并利用LabVIEW進(jìn)行數(shù)據(jù)采集與處理(以及儀器控制),可以研究的物理問題也比較豐富,因此是很好的大學(xué)物理綜合設(shè)計(jì)性實(shí)驗(yàn)的選題. 本文的實(shí)驗(yàn)裝置可以根據(jù)實(shí)際情況進(jìn)行簡化. 比如用普通的分束鏡和2塊偏振片代替PBS;用示波器代替數(shù)據(jù)采集卡,在XY模式下觀察測量點(diǎn)在橢圓上的運(yùn)動實(shí)現(xiàn)條紋的可逆計(jì)數(shù). 數(shù)據(jù)處理也可以根據(jù)具體情況做適當(dāng)?shù)暮喕? 比如在測量短時(shí)間過程時(shí),可以忽略激光器輸出狀態(tài)的漂移,不做光強(qiáng)標(biāo)度,從而省去第1個(gè)分束鏡和探測器.

    圖13 位移測量噪聲的統(tǒng)計(jì)直方圖

    致謝:感謝北京師范大學(xué)物理學(xué)系熊俊教授在本文寫作過程提供的建設(shè)性意見.

    猜你喜歡
    金屬桿光程干涉儀
    基于改進(jìn)的邁克爾遜干涉儀對熱變形特性的研究
    用于原子干涉儀的光學(xué)鎖相環(huán)系統(tǒng)
    用微元法解含金屬桿的電磁感應(yīng)電路問題
    基于傅里葉變換色散條紋法的實(shí)驗(yàn)研究*
    非對稱干涉儀技術(shù)及工程實(shí)現(xiàn)
    金屬桿“勻加速直線”切割磁感線的幾類問題解析
    基于最優(yōu)模糊的均勻圓陣干涉儀測向算法
    擺鏡掃描傅里葉變換光譜儀光程差計(jì)算
    嵌在墻上的自行車鎖
    激光全息照相景深擴(kuò)展實(shí)驗(yàn)探討
    一二三四在线观看免费中文在 | 久久99蜜桃精品久久| 中文天堂在线官网| 精品人妻偷拍中文字幕| 亚洲精品日韩在线中文字幕| 日韩一本色道免费dvd| 亚洲人成77777在线视频| 亚洲在久久综合| 久久久国产精品麻豆| 日韩熟女老妇一区二区性免费视频| 日本免费在线观看一区| 免费观看在线日韩| 国产免费福利视频在线观看| 18在线观看网站| 婷婷色综合大香蕉| 美女国产视频在线观看| 国产高清国产精品国产三级| 热re99久久国产66热| 久久国内精品自在自线图片| 亚洲国产成人一精品久久久| 国产毛片在线视频| 国产精品国产三级专区第一集| √禁漫天堂资源中文www| 免费日韩欧美在线观看| 午夜影院在线不卡| 欧美精品人与动牲交sv欧美| 搡老乐熟女国产| 狂野欧美激情性bbbbbb| 国产成人午夜福利电影在线观看| 91午夜精品亚洲一区二区三区| 国产精品免费大片| 亚洲五月色婷婷综合| 久久狼人影院| 亚洲精品视频女| 夫妻午夜视频| 天堂俺去俺来也www色官网| 少妇的逼水好多| 亚洲av欧美aⅴ国产| 免费播放大片免费观看视频在线观看| 草草在线视频免费看| 久久99蜜桃精品久久| 男女啪啪激烈高潮av片| 日韩中字成人| 日韩不卡一区二区三区视频在线| 亚洲精品色激情综合| 欧美性感艳星| 9191精品国产免费久久| 日韩 亚洲 欧美在线| 午夜精品国产一区二区电影| 少妇精品久久久久久久| kizo精华| 香蕉丝袜av| av在线观看视频网站免费| 精品一区二区免费观看| 激情五月婷婷亚洲| 成人无遮挡网站| 日韩一区二区三区影片| 久久久a久久爽久久v久久| 亚洲av电影在线进入| 国产麻豆69| 99久久综合免费| av线在线观看网站| 亚洲av电影在线进入| 色婷婷av一区二区三区视频| 亚洲国产色片| 高清视频免费观看一区二区| 国产欧美日韩综合在线一区二区| 夜夜骑夜夜射夜夜干| 在线观看一区二区三区激情| 纵有疾风起免费观看全集完整版| 丰满迷人的少妇在线观看| 熟女人妻精品中文字幕| 日韩制服丝袜自拍偷拍| 精品亚洲成a人片在线观看| 少妇的逼水好多| 波野结衣二区三区在线| 亚洲激情五月婷婷啪啪| 亚洲伊人久久精品综合| 99热国产这里只有精品6| 少妇的丰满在线观看| 午夜免费鲁丝| 热99国产精品久久久久久7| 精品99又大又爽又粗少妇毛片| 伦理电影大哥的女人| 国产精品秋霞免费鲁丝片| 丝袜美足系列| 久久婷婷青草| 波多野结衣一区麻豆| 日本爱情动作片www.在线观看| 另类精品久久| 最近中文字幕2019免费版| 18禁裸乳无遮挡动漫免费视频| 春色校园在线视频观看| kizo精华| 视频区图区小说| 咕卡用的链子| 久久久精品免费免费高清| 日韩av在线免费看完整版不卡| 狂野欧美激情性xxxx在线观看| 啦啦啦啦在线视频资源| 免费高清在线观看日韩| √禁漫天堂资源中文www| 日韩一区二区三区影片| 波野结衣二区三区在线| 欧美人与性动交α欧美精品济南到 | 国产色婷婷99| 色网站视频免费| 日本午夜av视频| 欧美国产精品一级二级三级| 高清欧美精品videossex| 日韩 亚洲 欧美在线| 免费黄网站久久成人精品| 热99久久久久精品小说推荐| 国产av精品麻豆| 亚洲av国产av综合av卡| 国产熟女午夜一区二区三区| 18禁裸乳无遮挡动漫免费视频| 黑人猛操日本美女一级片| 精品少妇久久久久久888优播| videossex国产| 国产精品久久久久久精品古装| 国产xxxxx性猛交| 国产免费现黄频在线看| 男人操女人黄网站| 国产亚洲欧美精品永久| 看免费av毛片| 亚洲人与动物交配视频| tube8黄色片| 亚洲经典国产精华液单| 免费观看a级毛片全部| 毛片一级片免费看久久久久| 久久精品aⅴ一区二区三区四区 | 桃花免费在线播放| 欧美精品亚洲一区二区| 免费av不卡在线播放| 国产一区亚洲一区在线观看| 韩国av在线不卡| www.av在线官网国产| 色婷婷av一区二区三区视频| 国产精品国产三级国产专区5o| 国产在线视频一区二区| 大片电影免费在线观看免费| 天天躁夜夜躁狠狠躁躁| 91午夜精品亚洲一区二区三区| 精品少妇久久久久久888优播| 伊人久久国产一区二区| 汤姆久久久久久久影院中文字幕| 国产免费一区二区三区四区乱码| 亚洲丝袜综合中文字幕| 精品国产一区二区久久| 亚洲欧美中文字幕日韩二区| 少妇高潮的动态图| 亚洲av福利一区| 麻豆精品久久久久久蜜桃| av有码第一页| 狂野欧美激情性xxxx在线观看| 18禁动态无遮挡网站| 欧美日韩综合久久久久久| 大片电影免费在线观看免费| 九九爱精品视频在线观看| 满18在线观看网站| 国产成人免费无遮挡视频| 免费黄频网站在线观看国产| 亚洲av中文av极速乱| 久久 成人 亚洲| 欧美97在线视频| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频 | 久久韩国三级中文字幕| 国产精品国产av在线观看| 一区在线观看完整版| av电影中文网址| 久久久欧美国产精品| 国产国语露脸激情在线看| 女性被躁到高潮视频| 亚洲国产精品999| 少妇人妻 视频| av女优亚洲男人天堂| 在线精品无人区一区二区三| 国产国拍精品亚洲av在线观看| 在线观看www视频免费| 久久国产精品大桥未久av| 亚洲欧洲日产国产| 丁香六月天网| 日韩中字成人| 秋霞在线观看毛片| 亚洲欧美一区二区三区黑人 | 免费看av在线观看网站| 亚洲人成77777在线视频| 国产精品欧美亚洲77777| 久久久久久久大尺度免费视频| 中文字幕人妻丝袜制服| 免费观看无遮挡的男女| 赤兔流量卡办理| 一区二区av电影网| 妹子高潮喷水视频| 亚洲av欧美aⅴ国产| 免费观看无遮挡的男女| 亚洲欧美成人精品一区二区| 国产男人的电影天堂91| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 欧美日韩国产mv在线观看视频| 国产午夜精品一二区理论片| 七月丁香在线播放| 女的被弄到高潮叫床怎么办| 91午夜精品亚洲一区二区三区| 嫩草影院入口| 18+在线观看网站| 日本vs欧美在线观看视频| 亚洲精品第二区| 国产男女超爽视频在线观看| 日韩成人av中文字幕在线观看| 国产成人精品一,二区| 国产精品无大码| 国产亚洲精品久久久com| 亚洲国产精品999| 18禁裸乳无遮挡动漫免费视频| 在线观看www视频免费| 免费观看在线日韩| 丝袜美足系列| 捣出白浆h1v1| 久久久久人妻精品一区果冻| 男人添女人高潮全过程视频| 成人国产麻豆网| 久久精品夜色国产| 日本免费在线观看一区| 精品少妇内射三级| 欧美人与善性xxx| 777米奇影视久久| 欧美日韩亚洲高清精品| 亚洲一码二码三码区别大吗| 热re99久久国产66热| 18+在线观看网站| 日韩大片免费观看网站| 天天躁夜夜躁狠狠躁躁| 亚洲色图 男人天堂 中文字幕 | 欧美日韩综合久久久久久| 精品第一国产精品| 久久精品aⅴ一区二区三区四区 | 成人18禁高潮啪啪吃奶动态图| 亚洲少妇的诱惑av| 久久精品国产自在天天线| 欧美日韩av久久| 90打野战视频偷拍视频| 久久99热这里只频精品6学生| 国产日韩欧美在线精品| 女人精品久久久久毛片| 免费黄网站久久成人精品| 亚洲一区二区三区欧美精品| 国产永久视频网站| 国产亚洲最大av| tube8黄色片| 26uuu在线亚洲综合色| 日韩 亚洲 欧美在线| 最近中文字幕2019免费版| 色视频在线一区二区三区| 男女无遮挡免费网站观看| 国产精品秋霞免费鲁丝片| 国产亚洲一区二区精品| 国产av码专区亚洲av| xxxhd国产人妻xxx| 亚洲,欧美精品.| 国产欧美另类精品又又久久亚洲欧美| 久久人人97超碰香蕉20202| 国产极品天堂在线| 色婷婷久久久亚洲欧美| 国产av码专区亚洲av| 亚洲 欧美一区二区三区| 精品人妻在线不人妻| 成人国产av品久久久| 中文字幕av电影在线播放| 欧美变态另类bdsm刘玥| 各种免费的搞黄视频| 国产爽快片一区二区三区| 蜜臀久久99精品久久宅男| 国产精品久久久久成人av| 成人免费观看视频高清| 美女xxoo啪啪120秒动态图| 精品视频人人做人人爽| 大码成人一级视频| tube8黄色片| 少妇人妻 视频| 免费看av在线观看网站| 91国产中文字幕| av女优亚洲男人天堂| 最近中文字幕2019免费版| av在线播放精品| 18+在线观看网站| 午夜激情久久久久久久| 丝袜喷水一区| 成人国语在线视频| 国产深夜福利视频在线观看| 免费高清在线观看视频在线观看| 免费观看性生交大片5| 三上悠亚av全集在线观看| 久久精品人人爽人人爽视色| 制服诱惑二区| av天堂久久9| 久久精品夜色国产| 免费在线观看完整版高清| 少妇被粗大猛烈的视频| 高清不卡的av网站| 国产一区有黄有色的免费视频| 天天躁夜夜躁狠狠久久av| 伦理电影免费视频| 99热6这里只有精品| 日韩成人伦理影院| 午夜av观看不卡| 久久久精品94久久精品| 下体分泌物呈黄色| 久久精品熟女亚洲av麻豆精品| 熟女电影av网| 日本vs欧美在线观看视频| 黄色视频在线播放观看不卡| 人人妻人人爽人人添夜夜欢视频| 日韩制服骚丝袜av| 99久国产av精品国产电影| 色5月婷婷丁香| 亚洲国产av新网站| 国产精品久久久久久精品电影小说| 人人妻人人爽人人添夜夜欢视频| 日韩制服骚丝袜av| 国产黄色视频一区二区在线观看| 久久精品国产亚洲av涩爱| 成人国语在线视频| 久久99热6这里只有精品| 国产男女超爽视频在线观看| 欧美少妇被猛烈插入视频| 久久久久久久国产电影| 国产成人精品无人区| 久久综合国产亚洲精品| 在线观看免费视频网站a站| av又黄又爽大尺度在线免费看| 亚洲少妇的诱惑av| 我的女老师完整版在线观看| 久久久a久久爽久久v久久| 亚洲综合色惰| tube8黄色片| 国产国语露脸激情在线看| 欧美亚洲日本最大视频资源| 中文字幕人妻熟女乱码| 久久久久久久国产电影| 草草在线视频免费看| 国产爽快片一区二区三区| 五月伊人婷婷丁香| 欧美日韩亚洲高清精品| 波多野结衣一区麻豆| 三级国产精品片| 国产永久视频网站| 欧美日韩av久久| 日本-黄色视频高清免费观看| 国产精品久久久久久av不卡| 久久久国产精品麻豆| 亚洲色图综合在线观看| 国产在线免费精品| 一二三四中文在线观看免费高清| 大香蕉久久网| 啦啦啦视频在线资源免费观看| 亚洲成av片中文字幕在线观看 | xxx大片免费视频| 亚洲精品一区蜜桃| 久久久久久伊人网av| 久久久久久人妻| 五月伊人婷婷丁香| 午夜精品国产一区二区电影| 欧美bdsm另类| 国产欧美日韩一区二区三区在线| 91国产中文字幕| 韩国av在线不卡| 午夜日本视频在线| 国产日韩欧美视频二区| 日产精品乱码卡一卡2卡三| 最近的中文字幕免费完整| 国产女主播在线喷水免费视频网站| 国产亚洲精品久久久com| 国产有黄有色有爽视频| 春色校园在线视频观看| 国产免费福利视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲精品日韩在线中文字幕| 国产免费一区二区三区四区乱码| 全区人妻精品视频| 好男人视频免费观看在线| 国产免费现黄频在线看| 咕卡用的链子| 日本-黄色视频高清免费观看| 性色avwww在线观看| 国产av精品麻豆| 亚洲国产最新在线播放| 精品熟女少妇av免费看| 亚洲精品日本国产第一区| 免费观看av网站的网址| 五月天丁香电影| 最近最新中文字幕大全免费视频 | 大香蕉久久成人网| 女的被弄到高潮叫床怎么办| 国产成人91sexporn| 乱码一卡2卡4卡精品| 日韩中文字幕视频在线看片| 久久久精品免费免费高清| 狂野欧美激情性xxxx在线观看| 久久亚洲国产成人精品v| 国产成人av激情在线播放| 91精品三级在线观看| 高清不卡的av网站| 插逼视频在线观看| 国产老妇伦熟女老妇高清| 亚洲av欧美aⅴ国产| 国产成人aa在线观看| freevideosex欧美| 18禁动态无遮挡网站| 国产麻豆69| 丝袜喷水一区| 欧美另类一区| 亚洲av成人av| 午夜日韩欧美国产| 久久人妻熟女aⅴ| 亚洲一区中文字幕在线| 黄片大片在线免费观看| 中文字幕制服av| 性少妇av在线| 久久精品熟女亚洲av麻豆精品| a级片在线免费高清观看视频| 国产日韩一区二区三区精品不卡| 99在线人妻在线中文字幕 | 老司机影院毛片| 精品卡一卡二卡四卡免费| 少妇裸体淫交视频免费看高清 | 久久香蕉精品热| 成年人黄色毛片网站| 日韩免费av在线播放| 免费不卡黄色视频| 久久亚洲精品不卡| 午夜免费鲁丝| 国产av精品麻豆| 久久久久国产精品人妻aⅴ院 | 大型黄色视频在线免费观看| 首页视频小说图片口味搜索| 少妇被粗大的猛进出69影院| 午夜免费鲁丝| 俄罗斯特黄特色一大片| 欧美精品亚洲一区二区| 大香蕉久久成人网| 国产高清视频在线播放一区| 久久午夜亚洲精品久久| 在线观看免费视频日本深夜| 亚洲精华国产精华精| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜制服| 免费一级毛片在线播放高清视频 | 热99国产精品久久久久久7| 丝袜美腿诱惑在线| 色综合婷婷激情| 国产欧美亚洲国产| xxx96com| 国产不卡av网站在线观看| 欧美成人免费av一区二区三区 | 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 在线天堂中文资源库| 精品久久久久久久久久免费视频 | 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 好看av亚洲va欧美ⅴa在| 久久香蕉国产精品| 最近最新中文字幕大全电影3 | 热99久久久久精品小说推荐| 久久久国产精品麻豆| 最近最新中文字幕大全电影3 | 日韩熟女老妇一区二区性免费视频| 亚洲午夜理论影院| 夜夜爽天天搞| 女人被躁到高潮嗷嗷叫费观| 狂野欧美激情性xxxx| 黄色女人牲交| 国产成人免费无遮挡视频| 一夜夜www| 在线观看免费视频日本深夜| 国产成人一区二区三区免费视频网站| 精品视频人人做人人爽| 精品少妇久久久久久888优播| 国产av精品麻豆| 九色亚洲精品在线播放| 午夜精品在线福利| 免费黄频网站在线观看国产| 91大片在线观看| 五月开心婷婷网| www日本在线高清视频| 黄片大片在线免费观看| 天堂俺去俺来也www色官网| 国产精品成人在线| 动漫黄色视频在线观看| 国产精品久久视频播放| 欧美性长视频在线观看| 黄频高清免费视频| 欧美乱色亚洲激情| 香蕉丝袜av| 美国免费a级毛片| 亚洲性夜色夜夜综合| 色综合婷婷激情| 老司机影院毛片| 欧美乱码精品一区二区三区| 久久中文字幕人妻熟女| 男女免费视频国产| 成人影院久久| 亚洲视频免费观看视频| 精品欧美一区二区三区在线| 国产一区有黄有色的免费视频| 在线观看免费午夜福利视频| 国产99白浆流出| 人人妻人人澡人人看| 欧美国产精品va在线观看不卡| 老鸭窝网址在线观看| 欧美日韩黄片免| 欧美日韩福利视频一区二区| 波多野结衣av一区二区av| 777米奇影视久久| 亚洲精品一二三| 国产精品永久免费网站| 精品乱码久久久久久99久播| 一级毛片女人18水好多| 日韩视频一区二区在线观看| 一边摸一边抽搐一进一出视频| 麻豆av在线久日| 国产一区二区三区视频了| 精品无人区乱码1区二区| 亚洲第一av免费看| 中出人妻视频一区二区| 国产亚洲精品第一综合不卡| 99热网站在线观看| 国产成人系列免费观看| 日韩欧美在线二视频 | 国产av又大| 悠悠久久av| 在线观看免费日韩欧美大片| 亚洲精品成人av观看孕妇| 国产成人影院久久av| 午夜视频精品福利| 欧美精品一区二区免费开放| 亚洲三区欧美一区| 久久九九热精品免费| 丰满饥渴人妻一区二区三| 亚洲五月婷婷丁香| 亚洲精品中文字幕在线视频| av福利片在线| 国产精品影院久久| 首页视频小说图片口味搜索| 美女午夜性视频免费| 多毛熟女@视频| 亚洲视频免费观看视频| 亚洲情色 制服丝袜| 天堂√8在线中文| 亚洲成人手机| 91av网站免费观看| 国产精品乱码一区二三区的特点 | 女人高潮潮喷娇喘18禁视频| 在线av久久热| 成年动漫av网址| 制服人妻中文乱码| 日本黄色日本黄色录像| netflix在线观看网站| 51午夜福利影视在线观看| 老司机在亚洲福利影院| 一边摸一边抽搐一进一小说 | 国产色视频综合| 欧美日韩黄片免| 99久久国产精品久久久| 涩涩av久久男人的天堂| 大片电影免费在线观看免费| 国产欧美日韩一区二区三| 看免费av毛片| 99精国产麻豆久久婷婷| 色老头精品视频在线观看| 高清毛片免费观看视频网站 | 色婷婷久久久亚洲欧美| 亚洲国产欧美一区二区综合| 在线观看www视频免费| 国产欧美日韩综合在线一区二区| 国产亚洲精品第一综合不卡| 亚洲avbb在线观看| 美女高潮喷水抽搐中文字幕| 一区二区三区国产精品乱码| 97人妻天天添夜夜摸| 青草久久国产| 欧美在线黄色| 中文字幕高清在线视频| 男人的好看免费观看在线视频 | 99国产精品一区二区蜜桃av | √禁漫天堂资源中文www| 国产精品自产拍在线观看55亚洲 | 老司机午夜十八禁免费视频| 一本一本久久a久久精品综合妖精| 在线十欧美十亚洲十日本专区| 午夜精品在线福利| 欧美黑人精品巨大| 好看av亚洲va欧美ⅴa在| 一区在线观看完整版| av视频免费观看在线观看| 香蕉久久夜色| 亚洲人成电影免费在线| 宅男免费午夜| a级毛片黄视频| 色婷婷av一区二区三区视频| 十八禁网站免费在线| 一二三四在线观看免费中文在| 首页视频小说图片口味搜索| 岛国毛片在线播放| 亚洲成人免费电影在线观看| 成人三级做爰电影| 成人亚洲精品一区在线观看| 国产99白浆流出| 国内久久婷婷六月综合欲色啪| 久久人妻熟女aⅴ| 黑人巨大精品欧美一区二区蜜桃|