999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

高精度非回轉(zhuǎn)對稱非球面加工方法研究

2016-11-09 05:38:47張春雷
中國光學(xué) 2016年3期
關(guān)鍵詞:檢測

徐 樂,張春雷,代 雷,張 健

(中國科學(xué)院 長春光學(xué)精密機(jī)械與物理研究所 應(yīng)用光學(xué)國家重點(diǎn)實(shí)驗(yàn)室超精密光學(xué)工程研究中心,吉林 長春130033)

?

高精度非回轉(zhuǎn)對稱非球面加工方法研究

徐樂*,張春雷,代雷,張健

(中國科學(xué)院 長春光學(xué)精密機(jī)械與物理研究所 應(yīng)用光學(xué)國家重點(diǎn)實(shí)驗(yàn)室超精密光學(xué)工程研究中心,吉林 長春130033)

本文提出一種高精度非回轉(zhuǎn)對稱非球面加工方法。首先,通過范成法銑磨出非回轉(zhuǎn)對稱非球面的最佳擬合球;然后,利用古典拋光修正小磨頭確定拋光難以修正的中頻誤差;最后,利用高精度氣囊拋光設(shè)備(IRP)精確對位精修面形,在不引入額外中頻誤差條件下,通過高精度對位檢測技術(shù)實(shí)現(xiàn)非回轉(zhuǎn)對稱非球面高精度加工。將該方法應(yīng)用于定點(diǎn)曲率半徑為970.737 mm、k=-1、口徑為106 mm三次非球面加工,降低了加工難度,提高了加工精度,面形誤差收斂到1/30λ(RMS)。實(shí)驗(yàn)結(jié)果驗(yàn)證了本文加工方法的正確性和可行性,對高精度非回轉(zhuǎn)對稱非球面加工具有一定的指導(dǎo)意義。

非回轉(zhuǎn)對稱;非球面;氣囊拋光;IRP拋光

1 引 言

為了校正像差、改善像質(zhì)、擴(kuò)大視場[1]并簡化系統(tǒng)結(jié)構(gòu)以達(dá)到減輕設(shè)備質(zhì)量、降低系統(tǒng)成本的目的,在光學(xué)系統(tǒng)中開始大量采用非球面,其在X射線光學(xué)系統(tǒng)、功能光電器件以及信息與微電子等尖端科技領(lǐng)域應(yīng)用較為廣泛[2]。由于不同系統(tǒng)的需要,非球面參數(shù)的選擇多種多樣。系統(tǒng)中大量采用了三次非球面,通過波前編碼技術(shù)實(shí)現(xiàn)對景深的延拓[3],而非回轉(zhuǎn)對稱非球面的加工更為困難。

非球面元件需要經(jīng)過車削、磨削及研磨等方式成型,這些加工方法會在其表面留有車痕、磨痕以及裂紋、殘留應(yīng)力等亞表面損傷[4-6],且加工過程中軸承噪聲/振動、控制精度及加工工具磨損等因素均會造成工件面形中頻誤差;工件表面/亞表面損傷及面形誤差會都會對系統(tǒng)的光學(xué)性能造成很大影響,因此需要通過拋光技術(shù)對工件表面進(jìn)行修正。通常對表面中頻誤差的修正方式為剛性大尺寸磨頭與隨機(jī)路徑結(jié)合加工,但都需要工件是旋轉(zhuǎn)對稱結(jié)構(gòu),也就是說這些普遍應(yīng)用于偶次非球面的方法只能對對稱面形中頻誤差進(jìn)行修正[7-11]。

由于非回轉(zhuǎn)對稱非球面的非對稱特性,普遍應(yīng)用于對稱非球面的加工方法不再適用,因此非回轉(zhuǎn)對稱非球面的加工極為困難。本文對非回轉(zhuǎn)對稱非球面加工方法進(jìn)行了研究,先解決對稱項面形誤差,再通過面形補(bǔ)償?shù)姆绞浇鉀Q了非回轉(zhuǎn)對稱非球面加工的困難。此方法對非回轉(zhuǎn)對稱非球面普遍適用,同樣也是一種自由曲面加工的新方法。

2 非回轉(zhuǎn)對稱非球面加工方法

通常情況下,非回轉(zhuǎn)對稱非球面可以表示為如下方程:

(1)

式中,x、y、z分別為x、y、z軸坐標(biāo),c為頂點(diǎn)球曲率,k為曲線的形狀參數(shù),D為歸一化半徑α,β為高次項常數(shù)[12]。如果β=0,那么此非球面方程可以改寫為:

(2)

式中,r2=x2+y2,從方程可以看出,這個非球面方程為旋轉(zhuǎn)對稱結(jié)構(gòu),因此可以通過螺旋線加工方式對表面面形誤差進(jìn)行修正。然而式(1)中,如果β≠0,那么面形為非旋轉(zhuǎn)對稱結(jié)構(gòu),將無法進(jìn)行螺旋線加工,如果采用IRP氣囊磨頭柵格路徑加工將無法消除銑磨過程中帶來的中頻誤差,因此將無法實(shí)現(xiàn)高精度光學(xué)元件面形加工(1/30λ,RMS)。

本文將非回轉(zhuǎn)對稱非球面分為兩個部分進(jìn)行加工,加工思路如下:第一步,將非球面的對稱部分加工完畢,由于面形為偶次項(式(3)中z1),是旋轉(zhuǎn)對稱結(jié)構(gòu),因此可以先通過IRP氣囊拋光(此處采用聚氨酯拋光膜)去除銑磨過程中帶來的亞表面損傷(約10 μm),聚氨酯拋光膜加工過程中雖然會引入中頻誤差,但具有較高加工效率,再采用古典拋光對工件中頻誤差進(jìn)行修正達(dá)到理想的對稱面形;第二步通過面形誤差補(bǔ)償?shù)姆绞剑跈z測得到的絕對誤差圖基礎(chǔ)上加入非回轉(zhuǎn)對稱項面形誤差,通過IRP準(zhǔn)對心氣囊磨頭(拋光膜選用Uninap)加工的方式得到最終非對稱非球面面形。

為了保證加工效率,此方法只對非對稱非球面有較好的加工效果,因?yàn)樵谏鲜龅诙街胁捎肬ninap拋光膜,此拋光膜的特點(diǎn)是不會引入中頻誤差(尤其對4 mm以下頻段控制良好),但是去除效率較低,因此如果非回轉(zhuǎn)對稱項于對稱項偏離量過大則會使得加工周期變長,影響加工效率。

本文選用的非球面參數(shù)如表1所示。

表1 本文選用的非球面參數(shù)

其中,R為工件半徑,Φ為通光口徑。非球面方程為方程(3):

(3)

根據(jù)方程可知對稱項為拋物面,同樣為非球面。按照上一節(jié)的敘述,可以將此非球面分為對稱項和非對稱項分別進(jìn)行面形修正。通過方程可以計算出對稱項的干涉圖形,如圖1所示。同樣可以計算出非回轉(zhuǎn)對稱項的面形誤差,如圖2所示。由此可見,非回轉(zhuǎn)對稱項的面形誤差為±1.1 μm,即加工時應(yīng)去除最大量為2.2 μm。非球面偶次項加工難度最大,因?yàn)榕即雾椡瑯邮欠乔蛎妗2⑶遥跈z測過程中,由于面形為非回轉(zhuǎn)對稱非球面,并且非回轉(zhuǎn)對稱項為一個45°方向的傾斜,因此在檢測時沒有基準(zhǔn)面進(jìn)行參考,也就無法進(jìn)行絕對測量。在測量中,只能通過調(diào)整光圈個數(shù)為最少作為每次測量的基準(zhǔn),因此這種方法既增大了檢測難度又增加了加工難度。

圖1 非球面偶次項面形干涉圖 Fig.1 Surface shape interferogram of even-order aspheric surface

圖2 非球面非回轉(zhuǎn)對稱項面形干涉圖 Fig.2 Surface shape interferogram of non-rotationally symmetrical aspheric surface

圖3 非球面與最佳擬合球面形差 Fig.3 Difference of surface shape interferogram between asphere and optimum fitting sphere

為了降低加工難度,我們將偶次項改為此偶次非球面的最佳擬合球。最佳擬合球的方法擁有更大優(yōu)勢,在加工方面,由圖可知此拋物面最大矢高為-1.432 9 mm,因此其最佳擬合球曲率半徑為971.460 mm,如圖3所示,最大面形誤差為2.2 μm,與前面所述方法去除量基本相同,球面的加工比非球面的加工更加容易,并且可以通過古典拋光方式修正表面中頻誤差,從而得到更好的面形誤差。在檢測方面,我們可以將此非球面當(dāng)做是一個擁有2.2 μm PV面形誤差、曲率半徑為971.460 mm的球面工件來進(jìn)行檢測。通過上述分析可知,這種方法降低了加工和檢測的難度,大大提高了非回轉(zhuǎn)對稱非球面加工的可行性,并且由于引入古典拋光從而提高了工件面形質(zhì)量。

3 加工過程

由上述分析可知,本實(shí)驗(yàn)在銑磨過程參數(shù)為:曲率半徑為971.460 mm,口徑為106 mm。在磨削加工過程中玻璃工件會產(chǎn)生大量的亞表面損傷[13],這些損傷將直接降低工件的強(qiáng)度、長期穩(wěn)定性、成像質(zhì)量、鍍膜質(zhì)量和抗激光損傷閾值等重要性能指標(biāo)[14]。因此,本實(shí)驗(yàn)的第一步就是將銑磨后的工件進(jìn)行快速拋光,既要對面形有一定的修正,使得拋光后可以進(jìn)行干涉儀檢測,從而進(jìn)行進(jìn)一步面形的精修,又要去除一定的厚度以達(dá)到去除亞表面損傷層的目的。本實(shí)驗(yàn)通過IRP200對工件進(jìn)行快速拋光。IRP采用CCOS技術(shù)控制氣囊磨頭進(jìn)行拋光,過Preston 方程描述氣囊拋光光學(xué)表面材料去除量與加工參數(shù)的關(guān)系[15]:

(4)

(5)

通過卷積定理可以算出駐留時間D(x,y),從而確定IRP氣囊磨頭加工過程每個點(diǎn)的加工速率,進(jìn)而實(shí)現(xiàn)向理想面形的收斂。拋光后結(jié)果如圖4所示。

圖4 快速拋光干涉圖 Fig.4 Interferogram of fast polishing

圖5 光順拋光干涉圖和1~3 mm帶通濾波結(jié)果 Fig.5 Inferferogram after smoothing polishing and filter result of 1-3 mm bandpass

從圖4可看出存在明顯柵格狀紋理,這種紋理將會形成中頻誤差,而IRP氣囊拋光對于這種中頻誤差沒有去除能力。由于此工件為非對稱非球面,加工時對位精度存在一定誤差,因此面形誤差較大。為了使面形誤差降低到1/30λ則中頻誤差必須?1/30λ,加工過程中需要引入古典拋光來對中頻誤差進(jìn)行修正。古典拋光后面形誤差為67.274 nm RMS,1~3 mm帶通濾波結(jié)果為0.903 nm RMS(如圖5(a)、(b)所示)。最終面形誤差目標(biāo)為1/30λ,即21 nm RMS。而1~3 mm頻段為氣囊拋光無法去除的頻段,即該頻段會一直殘留于最終面型,而此時的1~3 mm帶通濾波結(jié)果為0.903 nm RMS,已經(jīng)遠(yuǎn)遠(yuǎn)小于21 nm RMS,因此,此時1~3 mm頻段對最終結(jié)果影響甚微,可以進(jìn)入下一道工序。

圖6 需要面形精修的面形誤差 Fig.6 Error map of surface shape to be fine polished

圖7 最終面形誤差 Fig.7 Final error map of surface shape

為了得到最終面形,采用IRP200的Uninap氣囊拋光,此拋光膜特點(diǎn)為不會引入明顯的中頻誤差,并且加工精度較高。在加工過程中,由于去除量較大,因此在加工過程中如果固定加工軌跡將會帶來中頻誤差。我們可以給工件定義一個正方向,將工件以正方向?yàn)榛鶞?zhǔn)順時針分為多個方向,這些方向都加工并依次稱之為一輪加工。在加工過程中采取從工件多個方向依次走柵格路徑的方式進(jìn)行加工,每次加工采用的間隔不同,但是去除材料的量相同,這樣合成的材料去除量與預(yù)期相同,并且既避免了由于機(jī)床精度誤差帶來的象散,又可以屏蔽掉相同加工路徑大量去除過程中帶來的軌跡中頻誤差。將圖3與圖5(a)面形進(jìn)行疊加就會得到IRP200進(jìn)行修正的誤差圖,如圖6所示。利用IRP200進(jìn)行面形精修后的結(jié)果如圖7所示。圖7所示的面形誤差為18.224 nm(<1/30λ),已經(jīng)達(dá)到最終要求。

4 結(jié)果與討論

圖8 不同轉(zhuǎn)角時的面型誤差 Fig.8 Error map of surface shape rotated by different angles

非回轉(zhuǎn)對稱非球面檢測的角度對準(zhǔn)是非常重要的,由于面形本身的非對稱性,且IRP高精度拋光機(jī)床又為準(zhǔn)對心近似點(diǎn)對點(diǎn)加工,因此檢測的旋轉(zhuǎn)誤差將會對加工帶來錯誤的指導(dǎo)。本實(shí)驗(yàn)采用ZYGO的12寸VWS干涉儀進(jìn)行面形檢測,干涉儀的重復(fù)精度為0.02 nm,因此完全滿足檢測要求。根據(jù)干涉儀探測器像元尺寸和檢測鏡片口徑,測量旋轉(zhuǎn)角度偏差不超過±0.2°,因此只要在這一范圍內(nèi)的面形誤差遠(yuǎn)小于最終非球面面形誤差就可認(rèn)為基本滿足測量要求。如圖8(a)、(b)所示,經(jīng)過±0.2°旋轉(zhuǎn)后與原面形誤差相差PV值與RMS值一致,并且RMS值為0.968 nm?1/30λ,因此滿足測量過程中由于角度對準(zhǔn)帶來的面形誤差要求。

在角度偏轉(zhuǎn)給加工帶來的誤差的同時,工件的中高頻同樣需要關(guān)注。光學(xué)元件的表面信息是一個非常復(fù)雜的混合信號。傳統(tǒng)的峰谷值(PV)和均方根(RMS)等所表征的表面形態(tài)是有限的功率譜密度函數(shù)(PSD),是具有綜合分析意義的概率統(tǒng)計函數(shù),它可以從頻域上研究隨機(jī)振動的各頻率成分的統(tǒng)計含量[16]。這一方式可以定量地給出表面輪廓的空間頻率分布,從而確定中高頻段對整個面形的影響。圖9中給出了利用Uninap拋光膜加工前后的PSD曲線。Uninap加工后的曲線整體在加工前的下方,這說明Uninap對于各個頻段都沒有明顯的惡化作用,抑制了中高頻的同時將低頻面型修正到18.244 nm RMS。通過在氣囊磨頭上使用Uninap拋光膜走多方向柵格路徑的方式可以規(guī)避掉柵格路徑本身產(chǎn)生的中高頻,并且柵格路徑可以修正低頻面形。而進(jìn)入Uninap氣囊拋光階段之前為古典光順方式拋光,該拋光方法為螺旋路徑加工,兩種加工方式迭代時會出現(xiàn)彼此路徑相互垂直的階段,這一階段可以去除不同路徑產(chǎn)生的固有頻率,因此圖9中顯示為全頻段PSD曲線下降。

圖9 Uninap拋光膜加工前后PSD曲線 Fig.9 PSD curves before and after processing by Uninap

5 結(jié) 論

高精度非回轉(zhuǎn)對稱非球面加工不僅需要對元件面形精確加工,還需對檢測角度精度對準(zhǔn)。本文根據(jù)非回轉(zhuǎn)對稱非球面特點(diǎn),提出該類非球面優(yōu)化加工方案:銑磨出最佳擬合球,而后古典拋光去除中頻,最后通過面形精修的方式調(diào)整面形誤差,這種方法細(xì)化了加工過程,規(guī)避了加工風(fēng)險,降低了各個階段的加工難度,從而滿足最終面形誤差要求。采用該方式對曲率半徑為970.737 mm、k=-1、口徑為106 mm三次非回轉(zhuǎn)對稱非球面加工,檢測中精確對準(zhǔn)獲得并分析面形誤差,加工中精確對準(zhǔn)標(biāo)記并采用Uninap拋光膜不會引入新的中頻誤差,最終實(shí)現(xiàn)1/30λ(RMS)面形誤差加工,驗(yàn)證了本文加工工藝的正確性和可行性,可用于其它高精度非回轉(zhuǎn)對稱非球面加工。

[1]王貴林,李圣怡,戴一帆.光學(xué)非球面復(fù)合加工機(jī)床的設(shè)計與精度分析[J].中國機(jī)械工程,2004,15(2):99-102.

WANG G L,LI SH Y,DAI Y F. Design method and accuracy analysis of aspherical optical compound machine tool[J]. 中國機(jī)械工程,2004,15(2):99-102.(in Chinese)

[2]袁巨龍,吳喆,呂冰海,等.非球面超精密拋光技術(shù)研究現(xiàn)狀[J].機(jī)械工程學(xué)報,2012,48(23):167-177.

YUAN J L,WU ZH,LV B H,etal.. Review on ultra-precision polishing technology of aspheric surface[J].J.MechanicalEngineering,201248(23):167-177.(in Chinese)

[3]潘超,陳家璧,張榮福,等.用位相板實(shí)現(xiàn)景深延拓的原理與模擬實(shí)驗(yàn)研究[J].光子學(xué)報,2009,38(7):1679-1682.

PAN CH,CHEN J B,ZHANG R F,etal.. Study on principle and simulation of extending of depth of field with phase plate[J].ActaPhotonicaSinica,2009,38(7):1679-1682.(in Chinese)

[4]BEAUCAMP A,F(xiàn)REEMAN R,MORTON R,etal.. Removal of diamond-turning signatures on X-ray mandrels and metal optics by fluid-jet polishing[J].SPIE,2008,7018:351-359.

[5]GAO SH,KANG R K,GUO D M,etal.. Study on the subsurface damage distribution of the silicon wafer ground by diamond wheel[J].AdvancedMaterialsResearch,2010,126-128:113-118.

[6]SURATWALA T I,WONG L L,MILLER P E,etal.. Sub-surface mechanical damage distributions during grinding of fused silica[J].J.Non-CrystallineSolids,2006,352(52-54):5601-5617.

[7]張健,代雷,王飛,等.小磨頭自適應(yīng)拋光抑制高精度非球面中頻誤差[J].光學(xué)學(xué)報,2013,33(8):0822022.

ZHANG J,DAI L,WNAG F,etal.. Restraint of mid-spatial-frequency aspheric surface by small-tool aptive polishing[J].ActaOpticaSinica,2013,33(8):0822022.(in Chinese)

[8]高松濤,王高文,張健,等.用計算全息圖校正非球面的畸變[J].光學(xué) 精密工程,2013,21(8):1929-1935.

GAO S T,WANG G W,ZHANG J,etal.. Correction of distortion in asphere testing with computer-generated hologram[J].Opt.PrecisionEng.,2013,21(8):1929-1935.(in Chinese)

[9]曲興田,王宏一,樊成,等.重疊率螺旋線的非球面拋光軌跡規(guī)劃[J].西安交通大學(xué)學(xué)報,2015,49(6):y1-y7.

QU X T,WANG H Y,FAN CH,etal.. Uniform-overlap-rate spiral path for aspheric polishing[J].J.Xi′anJiaotongUniversity,2015,49(6):y1-y7.(in Chinese)

[10]周旭升,李圣怡,戴一帆,等.光學(xué)表面中頻誤差的控制方法—確定區(qū)域修正法[J].光學(xué) 精密工程,2007,15(11):1668 -1673.

ZHOU X SH,LI SH Y,DAI Y F,etal.. Correcting errors indefinite area: a new method for controlling mid-spatial-frequency errors in optical surface[J].Opt.PrecisionEng.,2007,15(11):1668-1673.(in Chinese)

[11]鄧偉杰,張峰,鄭立功.邊緣效應(yīng)的去除函數(shù)模型及實(shí)驗(yàn)[J].紅外與激光工程,2011,40(9):1473-1478.

DENG W J,ZHANG F,ZHENG L G. Removal functionmodel and experiment of edge effect[J].InfraredandLaserEngineering,2011,40(9):1743-1748.(in Chinese)

[12]孔繁星,于正林,樸承鎬,等.基于UMAC加工非球面光學(xué)零件PVT控制軌跡研究[J].長春理工大學(xué)學(xué)報,2009,32(2):196-206.

KONG F X,YU ZH L,PIAO CH H,etal.. J. Research on PVT control trajectory of processing aspheric optical parts based on UMAC[J].J.ChangchunUniversityofScienceandTechnology,2009,32(2):196-206.(in Chinese)

[13]胡陳林,畢果,葉卉,等.光學(xué)元件磨削加工亞表面損傷檢測研究[J].人工晶體學(xué)報2014,43(11):2929-2934.

HU CH L,BI G,YE H,etal.. Research on detection of subsurface damage on grinding optical elements[J].J.SyntheticCrystals,2014,43(11):2929-2934.(in Chinese)

[14]GENIN F Y,SALLEO A,PISTOR T V,etal.. Role of light intensification by cracks in optical breakdown on surfaces[J].J.AmericanCeramicSociety,2001,18(10):2607-2616.

[15]王貴林.光學(xué)材料超精密研拋關(guān)鍵技術(shù)研究[D].長沙:國防科學(xué)技術(shù)大學(xué),2002.

WANG G L. A study on key techniques in ultraprecision lapping and polishing for optical SIC materials[D]. Changsha:National University of Defense Technology,2002.(in Chinese)

[16]馬占龍,王君林.超高精度光學(xué)元件加工技術(shù)[J].紅外與激光工程,42(6):1485-1490.

MA ZH L,WANG J L. Ultra-precision optical fabrication technology[J].InfraredandLaserEngineering,42(6):1485-1490.(in Chinese)

Research on manufacturing method of non-rotationally symmetrical aspheric surface with high accuracy

XU Le*, ZHANG Chun-lei, DAI Lei, ZHANG Jian

(Engineering Research Center of Extreme Precision Optics,State KeyLaboratoryofAppliedOptics,ChangchunInstituteofOptics,FineMechanicsandPhysics,ChineseAcademyofSciences,Changchun130033,China)*Correspondingauthor,E-mail:xule198711@163.com

A novel manufacturing method of non-rotationally symmetrical aspheric surface with high accuracy is proposed in this paper. First, a optimum fit sphere of non-rotationally symmetrical asphere is ground through generating cutting, which is then polished through triditional polishing method in order to remove mid-spatial-frequency errors that cannot be realized through sub-aperture polishing tool. Finally, bonnet polishing method(IRP) is used for corrective polishing without generating extra mid-spatial-frequency errors in condition of high-accuracy positioning measurement technology. This method is employed in the manufacuring of a third-order aspheric surface with parameters of ROC=970.737 mm,k=-1 andD=106 mm. After several polishing iterations, the figure accuracy can be converged to 1/30λ(RMS). The correctness and feasibility of the manufacturing method are proved. This method would be meaningful in manufacturing non-rotationally symmetrical aspheric surface with small deviation in the future.

non-rotationally symmetrical;aspheric surface;bonnet polishing;IRP polishing

2016-01-26;

2016-02-19

國家科技重大專項資助項目(No.2009ZX02205)

2095-1531(2016)03-0364-07

TQ171.684

A

10.3788/CO.20160903.0364

徐樂(1987—),男,吉林長春人,碩士,研究實(shí)習(xí)員,2007年、2014年于北京理工大學(xué)分別獲得學(xué)士、碩士學(xué)位,主要從事光學(xué)精密加工方面的研究。E-mail:xule198711@163.com

Supported by National Major S&T Special Project of China(No.2009ZX02205)

猜你喜歡
檢測
QC 檢測
“不等式”檢測題
“一元一次不等式”檢測題
“一元一次不等式組”檢測題
“幾何圖形”檢測題
“角”檢測題
“有理數(shù)的乘除法”檢測題
“有理數(shù)”檢測題
“角”檢測題
“幾何圖形”檢測題
主站蜘蛛池模板: 激情网址在线观看| 亚洲天堂网在线观看视频| 国产精品第| 久久精品亚洲专区| 午夜高清国产拍精品| 色综合天天娱乐综合网| 韩日免费小视频| av免费在线观看美女叉开腿| 日本少妇又色又爽又高潮| 欧美第一页在线| 精品国产91爱| 国产视频入口| 天天操精品| 中文字幕人妻无码系列第三区| 97成人在线视频| 91外围女在线观看| 99成人在线观看| 2021亚洲精品不卡a| 成人综合在线观看| 欧美97欧美综合色伦图| 欧美精品黑人粗大| 成年人国产视频| 青青草91视频| yjizz视频最新网站在线| 久久精品无码专区免费| 国产91精品调教在线播放| 久久精品国产精品一区二区| 国产女人水多毛片18| 在线国产毛片| 亚洲国产精品日韩av专区| 成人在线亚洲| 国产黄色片在线看| 中文字幕在线播放不卡| 国产激情无码一区二区APP| 亚洲天堂网视频| 亚洲三级电影在线播放| 国产成人三级| av无码久久精品| 免费一级无码在线网站| 国产精品夜夜嗨视频免费视频| 在线播放国产99re| 成人免费视频一区二区三区| 人妻无码中文字幕一区二区三区| 国产剧情一区二区| 精品国产一区二区三区在线观看| 日本精品αv中文字幕| 欧美成人影院亚洲综合图| 国产精品福利尤物youwu| 国产精品亚洲天堂| 2024av在线无码中文最新| 亚洲第一极品精品无码| 免费国产小视频在线观看| 激情综合网激情综合| 国产性猛交XXXX免费看| 无码综合天天久久综合网| 久久综合九色综合97婷婷| 色视频久久| 精品亚洲麻豆1区2区3区| 成人在线不卡视频| 日韩欧美在线观看| 国产人成在线视频| 午夜日b视频| 91综合色区亚洲熟妇p| 污视频日本| 久久黄色视频影| 91久久精品日日躁夜夜躁欧美| 国产av色站网站| 亚洲精品麻豆| 91青青视频| 香蕉网久久| 免费国产高清视频| 午夜国产不卡在线观看视频| 在线日韩日本国产亚洲| 蜜芽一区二区国产精品| 九九香蕉视频| 国产幂在线无码精品| 在线观看亚洲精品福利片| 日韩无码黄色| 精品视频一区在线观看| 好吊妞欧美视频免费| 欧美色伊人| 亚洲欧美日本国产专区一区|