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

超聲橢圓振動車削三維形貌形成研究

2017-11-10 00:40:36張國華李咚咚李茂偉張德遠
兵工學報 2017年10期
關鍵詞:方向振動

張國華, 李咚咚, 李茂偉, 張德遠

(1.北京衛星制造廠, 北京 100094;2.北京航空航天大學 機械工程及自動化學院, 北京 100191)

超聲橢圓振動車削三維形貌形成研究

張國華1, 李咚咚1, 李茂偉1, 張德遠2

(1.北京衛星制造廠, 北京 100094;2.北京航空航天大學 機械工程及自動化學院, 北京 100191)

為研究橢圓振動車削表面三維幾何形貌形成規律及影響因素,采用金屬切削理論對橢圓振動切削過程和三維切削模型進行了分析。研究結果發現:切削過程中相鄰兩轉之間不同的相位差特征值對最終表面形貌有著重要的影響;通過已有的試驗結果表明,相位差特征值分別為0和0.5時所形成的表面形貌截然不同,與理論分析結果相吻合。

機械制造工藝與設備; 超聲橢圓振動; 表面形貌; 精密加工

0 引言

超聲橢圓振動車削是一種新型加工方法[1],具有降低切削力[2-4]、提高零件系統的穩定性[5-7]、提高表面質量的優勢,特別是在航空航天、兵器工業和難加工材料精密及超精密等加工領域內發揮了重要的作用。

然而,從目前絕大部分的研究成果和文獻可以看出,針對超聲橢圓振動車削表面形貌的研究還比較少,現有的研究主要基于實驗手段來驗證表面質量提高,或是基于振動往復熨壓的切削特性來分析研究表面成因[8-12]。對于橢圓振動車削表面微觀幾何形貌形成的研究領域幾乎處于空白,這束縛了超聲橢圓振動車削在精密及超精密等加工領域內的進一步推廣應用。

本文針對超聲橢圓振動精密加工過程中表面微觀幾何形貌形成進行了詳細研究,通過已有的實驗結果進行了驗證,驗證了分析的準確性。

1 橢圓振動車削三維形貌形成機理

1.1 表面幾何形貌理論公式的推導

對相位差不存在時橢圓振動車削表面形成機理前人已進行了分析,但是對相位差存在時表面幾何形貌形成機理還沒有系統的研究,故本文對表面幾何形貌形成機理進行全面的研究。

選擇帶圓弧的刀具,取工件一個區域進行分析。普通車削加工時,圓弧刀具在已加工表面形成了圖1所示的進給刀紋紋路。如果進給量fv滿足fv≤2Rεsink′r,k′r為刀具的副偏角,則進給刀紋最大高度由幾何關系[13]可知:

(1)

式中:Rε為刀尖圓弧半徑。

圖1 工件車削位置Fig.1 Cutting position of workpiece

由(1)式可知,進給刀紋的幾何高度僅與fv、Rε有關。精密加工時,由于兩個參數都是在同一數量級,因此fv對H的影響相比Rε的影響更大一些。

橢圓振動車削過程中,刀具同時存在兩個方向相互垂直振動:其中一個振動方向平行于切削速度,使得刀具產生周期性的分離運動;另一個振動方向垂直于切削速度,在切削方向上形成一系列具有一定高度的振紋,進給方向上的刀紋和切削方向上的振紋近似于相互垂直。在切削的過程中兩者相互影響,彼此有一定的削弱作用,最終形成了橢圓振動車削所特有的表面形貌,其中垂直于切削速度的振動對已加工面的表面幾何高度有著重要的影響。

超聲振動對系統顫振有著極強的抑制作用,故不考慮系統顫振的影響,從幾何角度來分析表面形貌的形成規律。如圖2所示,取圓弧刀刃與進給刀紋相切的一點G進行分析,任意時刻刀具橢圓振動控制方程[12]如(2)式所示:

(2)

式中:AXmax為X坐標進給方向最大振幅值;AX為X坐標進給方向任意振幅值;AYmax為Y坐標進給方向最大振幅值;AY為Y坐標進給方向任意振幅值;w為超聲橢圓振動加工頻率;θ為刀具振動相位角。

圖2 刀具- 工件位置關系Fig.2 Position of cutting tool and workpiece

刀具由橢圓軌跡不同振動點開始切入時,由于刀具振動形式不同,切入工件后所形成的表面形貌也有所不同。工件轉動一圈,刀具從A點到B點的切削行程為L,定義為刀具的一次切削行程,如圖1所示。圖2及圖3中的工件為沿周向的展開圖,對比之前的研究方法,引入頻率行程速度系數KL:

KL= 2wL/v=m+n,

(3)

式中:v為刀具本身切削速度(不含振動速度);m為整數;n為小數或0.

圖3 刀具任意振動狀態切入工件Fig.3 Cutting tool cutting into workpiece in any vibrational state

規定刀具一個完整的振動形態是走完一個完整的橢圓弧長s,橢圓軌跡走向沿進給方向看從起始點A點順時針再回到A點,如圖3所示。將n定義為相位差特征值,表達公式為

(4)

式中:s′(θ)是從起始點A點順時針振動到任一點的弧長。

對圖1、圖3進行分析,假如在B點完成車削時,n=0,則KL只有整數部分m,那么以A點瞬時振動形態作為參考,可知在切削行程L距離內刀具運動到B點時完成了m個的完整橢圓軌跡,即在A點與B點刀具的瞬時振動形態完全一致,由(2)式可知,刀具橢圓振動控制方程完全一致,故相鄰連轉不存在相位差;如果在B點完成車削時,n≠0,則KL由整數部分m和小數部分n構成的,那么以A點瞬時振動形態作為參考,可知在切削行程L距離內刀具運動到B點時完成了m個完整橢圓軌跡和n個橢圓軌跡,即A點與B點刀具的瞬時振動形態不一致,由(2)式可知,刀具橢圓振動控制方程不一致,必然造成A點的相位θ值與B點的相位θ值不等,故相鄰連轉存在相位差,因此其中n可以表征反映出相鄰兩轉之間的相位差。

1.2 相位差不存在時幾何形貌形成分析

刀具起振點位于橢圓軌跡上A點時,刀具開始切入工件,A點軌跡參數為

(5)

切完第1刀后,第2刀車削時,刀具由橢圓軌跡上B點開始切入工件,A點與B點刀具的瞬時振動形態完全一致,如(5)式所示,即相位差特征值n=0.A點與B點不存在相位差,在平行于進給方向的任意切平面的位置上,相鄰兩刀的橢圓最低點A′點和B′點必然在同一水平位置,如圖4所示。已加工表面的形貌如圖5(a)所示,其俯視圖投影后得到圖5(b)。

圖4 在橢圓最低點切削后截面形狀Fig.4 Shape at the lowest elliptical point

圖5 橢圓振動車削三維形貌(n=0)Fig.5 3-D micro-surface shapes in ultrasonic elliptical vibration cutting(n=0)

該類型的表面形成機理基于工件相鄰兩刀車削時刀具以相同振動狀態切入,李勛等[14]對此已經作了相應的研究,并對成型機理及影響因素進行了詳細分析,表面幾何高度值為車削形成的最高點至最低點的距離:H=RCC+REVC+R′CC,其中RCC為普通切削表面幾何高度,REVC為超聲切削表面幾何高度,R′CC為普通切削與超聲切削相互削弱的表面幾何高度值。

1.3 相位差存在時幾何形貌形成分析

但是大多數橢圓振動車削過程中,不可避免地存在相位差,這使得表面幾何形貌形成過程更加復雜。假定第1刀車削時,刀具由橢圓軌跡上任意形態A點車削,第2刀車削時,刀具由橢圓軌跡上任意形態B點開始切入工件,即相位差特征值n≠0,如圖3所示。

這種車削狀態下形成如圖6所示的三維表面形貌,由刀具Y方向振動和進給量的影響而形成的表面形貌幾何高度H滿足:H=RCC+REVC,圖6中給出了第1刀沿橢圓1車削分別形成波谷1和波峰3,第2刀沿橢圓2車削分別形成波谷2、波峰1和波峰2. 形成過程可以描述如下:刀具沿橢圓1切入后,受到垂直于切削速度方向的Y方向振動影響,形成一系列的振紋高點;刀具沿橢圓2切入后,受到相位差及進給量的影響,相鄰兩刀的刀具軌跡形成重疊,第2刀切削就會對第1刀切削所形成振紋高點的重疊區域進行重復切削,該高點被切去一部分,從而振紋高度有所降低。

圖6 橢圓振動車削三維形貌(n≠0)Fig.6 3-D micro-surface shape in ultrasonic elliptical vibration cutting(n≠0)

如果相位差特征值n由0逐漸增大時,波谷1就會在波峰1與波峰2之間變動游走,很明顯地發現波谷1與波谷2重合時,波峰3與波峰1重合,此時重復切削量最少,表面形貌幾何高度H最大,振動形態相同時,相位差特征值n為0;隨著相位差特征值n逐漸增大到0.5時,表面形貌幾何高度H值逐漸降低,當波谷1與波峰2重合,波峰3與波谷2重合,此時重復切削量最多,H值達到最??;隨著相位差特征值n繼續增大到1時,表面形貌幾何高度H值又開始逐漸增大,當波峰3與波峰2重合,此時重復切削量最少,表面形貌幾何高度H最大,此時相位差特征值n=1. 據此分析,相位差特征值n=0時,表面形貌幾何高度H值最大;相位差特征值n=0.5時,表面形貌幾何高度H值最??;相位差特征值n為其他值時,表面形貌幾何高度H值介于這兩者之間。

相位差特征值n=0.5時刀具位于兩個點的相對位置如圖7(a)所示,由橢圓軌跡上A點開始切入工件,A點軌跡參數為

(6)

切完第1刀后,第2刀車削時,刀具由橢圓軌跡上B點開始切入工件,B點參數表達同(5)式。

圖7 橢圓振動車削三維表面形貌(n=0.5)Fig.7 3-D micro-surface shapes in ultrasonic elliptical vibration cutting(n=0.5)

圖7(a)所處的加工位置開始加工,所形成的最終表面形貌類似圖6,將圖6在俯視圖方向投影得到如圖7(b)所示圖形。

此外,不同的相位差特征值n所得到最終三維形貌基本類似圖6所示,只是第1刀波峰和第2刀所形成的波谷相對位置不同,相對位置沿著X方向發生了錯動。針對普通狀況,列出了4種不同的相位差特征值n車削后所形成的表面形貌,俯視平面圖如圖8所示。

圖8 橢圓振動車削三維表面形貌Fig.8 3-D micro-surface shapes in ultrasonic elliptical vibration cutting

(7)

圖9 相鄰兩轉刀具軌跡幾何關系Fig.9 Geometrical relationship of cutting tool between two adjacent turns

將上述幾何關系代入(5)式,整理可得:

(8)

當振動頻率w=0時,REVC=0,有H=RCC,(8)式退化為

(9)

(9)式為普通車削表面形貌幾何高度公式。故(8)式將超聲振動車削與普通車削所形成的表面形貌幾何高度理論公式統一起來。橢圓振動車削刀具運動軌跡控制方程式如下:

(10)

分析圖10所示軌跡,為方便研究刀具在一個周期內的車削過程,不妨令刀具切入點位于Ps點,刀具位于P0點時,刀具Y負方向位于最大振幅位置,X正方向的振動速度達到大。此后X方向的振動速度開始逐漸減小,在Y方向的振幅為0時,振動速度減小至0. 此后刀具X方向振動速度開始增大,振動方向發生改變與工件運動的方向一致。當刀具運動到臨界點P1點時,刀具X方向的速度達到臨界速度,刀具開始與工件分離,刀具在下一次切入之前一直保持空切狀態。此后刀具X方向的振動速度繼續增大,當位于Y正方向最大振幅位置時,X負方向的振動速度達到最大。此后振動速度開始逐漸減小,在Y方向的振幅為0時,X負方向振動速度減小至0. 此后刀具X方向振動速度開始增大,振動方向發生改變與工件運動方向相反。由于對稱性,刀具運動到第2個臨界點P2點時,刀具X正方向的速度達到臨界速度。此后刀具X正方向的振動速度繼續增大,刀具運動到P3點時,刀具再次切入工件。至此刀具完成一個周期的車削,其中分離點和切入點如圖10所示。取P0點刀具軌跡時間和位置作相對起點,則此刻P1點、P2點、P3點所對應的時間及位置分別為t1、t2、t3及X1、X2、X3,圖10中P1點和P2點滿足對稱性,則有

(11)

圖10 刀具切入與分離過程Fig.10 Cutting-in of cutting tool and its separation from workpiece

兩個分離點速度方向滿足以下條件:

(12)

P1點和P2點滿足對稱性:

(13)

聯立(11)式、(12)式、(13)式,可求H.

2 橢圓振動切削參數對表面形貌的影響

2.1 頻率行程速度系數的影響

利用數值解析法求KL與H的關系,如圖11所示。

圖11 頻率行程速度系數- 表面幾何高度曲線Fig.11 Relation between frequency travel speed coefficient and surface geometric height

圖11中分離點左側區域為分離區,右側區域為不分離區。分析圖11,在分離區:其他參數確定后,表面形貌幾何高度H值隨著KL的增大而逐漸降低,并且降低的趨勢逐漸變緩,最終接近某一個數值并不會一直減小至0. 因此在滿足車削加工的各項指標時,可以采用盡可能大的頻率行程速度系數KL. 但是隨著KL值增大到一定范圍后,試圖通過繼續增加KL值來改善表面加工質量,效果會變的很微弱。隨著KL值的增大,形成的表面形貌幾何H值越來越低,這符合圖11中的趨勢。

2.2 切削參數的影響

為了進一步研究其他切削參數對表面形貌幾何高度的影響,根據(8)式,作出了進給量—表面幾何高度關系曲線、Y方向振幅—表面幾何高度關系曲線、刀尖圓弧半徑—表面幾何高度關系曲線,如圖12、圖13及圖14所示。

圖12 進給量—表面幾何高度曲線Fig.12 Relation between feed and surface geometric height

圖13 Y方向振幅—表面幾何高度曲線Fig.13 Relation between Y amplitude and surface geometric height

圖14 刀尖圓弧半徑—表面幾何高度曲線Fig.14 Relation between cutting tool’s circular arc radius and surface geometric height

從圖12可以看出:頻率行程速度系數KL確定后,隨著進給量fv值的增加,表面形貌幾何高度H逐漸增大,并且當fv接近2Rε時,H值逐漸接近于刀尖圓弧半徑Rε. 從圖12中還可以知道,KL值逐漸增大可以在抑制由于進給量fv值增加而造成的已加工表面量惡化。但是頻率行程速度系數KL增到一定范圍后,這種抑制作用就變得不明顯了。同時也可看出進給量fv值對H的影響要大于頻率行程速度系數KL值。

圖13中分離點左側區域為相應參數下的分離區,右側區域為不分離區。從圖13可以看出:在分離區,頻率行程速度系數KL值確定后,隨著Y方向振幅AY的增加,表面形貌幾何高度H逐漸減小;當AY值增大到幾十微米后,H值逐漸接近某一固定數值。事實上對于超聲振動車削,由于受到現有材料及換能器設計原理的限制,刀具單邊振幅很難達到很大,因此很難作到已加工表面的H值接近圖13中比較理想的狀態。對比KL=400.5、KL=700.5和KL=1 000.5這3條曲線,可以看出,KL值越大,H值對Y方向振幅變化越不敏感。KL值比較大的時候,試圖通過調整Y方向振幅大小來改善表面形貌,效果并不是很明顯。在不分離切削區,Y方向振幅增大,表面形貌幾何高度H隨之增大。因此當刀具處于不分離區進行切削時,為獲得比較好的表面形貌,Y方向振幅應盡可能的小。

從圖14可以看出,表面形貌幾何高度H隨著刀尖圓弧半徑增大而減小,頻率行程速度系數KL越大,刀尖圓弧半徑的這種影響也就越弱。

3 橢圓振動車削試驗驗證

圖15 雙彎復合橢圓振動系統原理圖Fig.15 Double bending composite elliptical vibration system

為了更清楚地驗證第2節相位差對表面三維形貌影響的準確性,以李勛等[14]所做的車削試驗結果來驗證前面分析的正確性,車削試驗系統為雙彎橢圓振動系統及試驗系統,振動原理如圖15、圖16所示。

圖16 車削試驗系統圖Fig.16 Vibration cutting test system

橢圓振動車削過程中主要的切削條件及參數如表1所示。

表1 超聲橢圓振動車削試驗條件

依據李勛等[14]試驗所得到的表面形貌照片如圖17所示。在頻率行程速度系數KL中的相位差特征值n=0時,對比理論分析所得到圖5(a)與試驗所得到的圖17(a),可以發現兩者的表面形貌一致性非常好;相位差特征值n=0.5時,對比理論分析圖5(b)與試驗所得到的圖17(b),可以發現兩者的表面形貌一致性非常好。

圖17 不同頻率行程速度系數的微觀幾何形貌Fig.17 Micro-surface shapes at different frequency travel speeds

通過相位差特征值n取兩種不同值的對比分析,驗證了第2節分析的頻率行程速度系數KL對表面幾何形貌影響的正確性,也說明了在實際加工過程中相位差特征值對微觀表面幾何形貌的形成有著重要的影響。

4 結論

1)通過理論分析,推導出了超聲橢圓振動表面形貌理論公式,將普通車削與超聲振動車削理論公式統一起來;橢圓振動車削時,切削表面幾何形貌是由進給刀紋和Y方向振動產生的振紋彼此削弱、共同作用而形成的。

2)最終的表面形貌幾何高度H值是由頻率行程速度系數KL、進給量、Y方向振幅和刀尖圓弧半徑共同決定的,每個參數對H值的影響各有不同。

3)通過試驗驗證了相位差特征值n是否為0對表面幾何形貌有著重要的影響,兩者呈現出截然不同的表面形貌。

References)

[1] 隈部淳一郎. 精密加工- 振動切削基礎與應用[M]. 韓一昆,譯. 北京: 機械工業出版社, 1985.

Kumabe J.Fundamentals and application of vibration cutting[M].HAN Yi-kun,translated.Beijing:China Machine Press,1985.(in Chinese)

[2] Ma C X,Shamoto E,Moriwaki T.Study on the thrust cutting force in ultrasonic elliptical vibration cutting [J].Materials Science Forum, 2004: 396-400, 467-470.

[3] 張德遠, 張成茂. 飛機交點孔超聲橢圓振動精密加工技術[J]. 中國機械工程, 2012, 23(1):39-41.

ZHANG De-yuan, ZHANG Cheng-mao.Ultrasonic elliptical vibration precision machining for aircraft intersection holes [J].China Mechanical Engineering, 2012, 23(1):39-41.(in Chinese)

[4] Zhou X Q, Zuo C M, Liu Q, et al. Surface generation of freeform surfaces in diamond turning by applying double-frequency elliptical vibration cutting[J]. International Journal of Machine Tools & Manufacture, 2016, 104:45-57.

[5] Zhang J, Suzuki N, Wang Y, et al. Ultra-precision nano-structure fabrication by amplitude control sculpturing method in elliptical vibration cutting[J]. Precision Engineering, 2015, 39:86-99.

[6] Kim G D,Loh B G.An ultrasonic elliptical vibration cutting device for micro V-groove machining:kinematical analysis and micro V-groove machining characteristics[J].Journal of Materials Processing Technology,2007, 190(1):181-188.

[7] 馬春翔,社本英二,森肋俊道,等.超聲波橢圓振動切削提高加工系統穩定性的研究[J].兵工學報,2004,25(6):752-756.

MA Chun-xiang,Shamoto E,Moriwaki T,et a1.Study on the improvement of machining system stability by applying ultrasonic elliptical vibration cutting[J].Acta Armamentarii,2004,25(6):752-756.(in Chinese)

[8] Jung H, Hayasaka T, Shamoto E. Mechanism and suppression of frictional chatter in high-efficiency elliptical vibration cutting[J]. CIRP Annals-Manufacturing Technology, 2016, 65(1):369-372.

[9] Lin J, Han J, Zhou X, et al. Study on predictive model of cutting force and geometry parameters for oblique elliptical vibration cutting[J]. International Journal of Mechanical Sciences, 2016, 117:43-52.

[10] Nath C,Rahman M,Neo K S.A study on the effect of tool nose radius in ultrasonic elliptical vibration cutting of tungsten carbide [J].Journal of Materials Processing Technology,2009,209(17):5830-5836.

[11] Kurniawan R, Kiswanto G, Ko T J. Micro-dimple pattern process and orthogonal cutting force analysis of elliptical vibration texturing[J]. International Journal of Machine Tools & Manufacture, 2016, 106:127-140.

[12] 張成茂,張德遠.導套式橢圓超聲鏜削模擬實驗研究[J].兵工學報,2013,34(4):465-470.

ZHANG Cheng-mao, ZHANG De-yuan.Simulation experiment study on elliptical ultrasonic vibration boring of leading bushing [J].Acta Armamentarii,2013,34(4):465-470.(in Chinese)

[13] 陳日曜. 金屬切削原理[M]. 第2版.北京: 機械工業出版社,2012.

CHEN Ri-yao.Metal-cutting principie [M].2nd ed. Beijing:China Machine Press,2012.(in Chinese)

[14] 李勛, 張德遠. 超聲橢圓振動切削表面形貌形成機理的研究[J]. 中國機械工程, 2009, 2(7): 807-811.

LI Xun,ZHANG De-yuan.Research on micro-surface formation mechanism of ultrasonic elliptical vibration cutting [J].China Mechanical Engineering,2009, 2(7): 807-811.(in Chinese)

ResearchonFormationof3-DMicro-surfaceinUltrasonicEllipticalVibrationCutting

ZHANG Guo-hua1, LI Dong-dong1, LI Mao-wei1, ZHANG De-yuan2
(1.Beijing Satellite Manufacturing Factory,Beijing 100094, China;2.School of Mechanical Engineering and Automation, Beihang University, Beijing 100191, China)

In order to find the influence factors and the influence law of the 3-D micro-surface formation of ultrasonic elliptical vibration cutting, the elliptical vibration cutting process and 3-D cutting model are analyzed based on the metal cutting theory. The results show that the different phase difference eigenvalue between two adjacent turns in vibration cutting process has an important influence on the final surface geometry. The existing cutting test results show that the final surface shapes with the eigenvalues of 0 and 0.5 are very different, and the results agree with the theoretically analyzed results.

manufacturing technology and equipment; ultrasonic elliptical vibration; surface morphology; precision machining

2017-04-06

張國華(1985—),男,工程師。E-mail:zgh650@163.com

張德遠(1963—),男,教授,博士生導師。E-mail:zhangdy@buaa.edu.cn

TG506.5

A

1000-1093(2017)10-2002-08

10.3969/j.issn.1000-1093.2017.10.017

猜你喜歡
方向振動
振動的思考
科學大眾(2023年17期)2023-10-26 07:39:14
噴水推進高速艇尾部振動響應分析
2022年組稿方向
計算機應用(2022年2期)2022-03-01 12:33:42
2022年組稿方向
計算機應用(2022年1期)2022-02-26 06:57:42
2021年組稿方向
計算機應用(2021年4期)2021-04-20 14:06:36
2021年組稿方向
計算機應用(2021年3期)2021-03-18 13:44:48
2021年組稿方向
計算機應用(2021年1期)2021-01-21 03:22:38
This “Singing Highway”plays music
振動攪拌 震動創新
中國公路(2017年18期)2018-01-23 03:00:38
中立型Emden-Fowler微分方程的振動性
主站蜘蛛池模板: 成人亚洲国产| 亚洲啪啪网| 国产传媒一区二区三区四区五区| 亚洲成年人片| 中国丰满人妻无码束缚啪啪| 天堂网亚洲系列亚洲系列| AV老司机AV天堂| 91在线免费公开视频| 精品视频一区二区观看| 一区二区欧美日韩高清免费| 亚洲天堂网在线播放| 国产综合色在线视频播放线视| 久久精品这里只有精99品| 国产玖玖玖精品视频| 亚瑟天堂久久一区二区影院| 国产精品视频观看裸模| 亚洲成av人无码综合在线观看| 中文字幕在线看视频一区二区三区| 欧美不卡在线视频| 五月婷婷丁香综合| 天天干天天色综合网| 色婷婷在线播放| 亚洲欧美日本国产综合在线| 自偷自拍三级全三级视频 | 欧美一级高清免费a| 色一情一乱一伦一区二区三区小说| 老色鬼欧美精品| 高h视频在线| 全裸无码专区| 国产午夜精品鲁丝片| 成人福利视频网| 四虎国产精品永久一区| 亚洲日韩精品伊甸| 九色国产在线| 国产精品精品视频| 久久精品人人做人人爽| 日韩成人高清无码| 久久精品国产国语对白| 亚洲精品国产日韩无码AV永久免费网 | 亚洲国产欧美国产综合久久| 国产制服丝袜91在线| 一本大道在线一本久道| 青青青视频91在线 | 亚洲日本一本dvd高清| 日韩毛片免费视频| 亚洲小视频网站| 国内精品91| 亚洲人成色77777在线观看| 精品撒尿视频一区二区三区| 小说区 亚洲 自拍 另类| 成年人视频一区二区| 国内自拍久第一页| 亚洲香蕉在线| 欧美午夜在线观看| 亚洲香蕉伊综合在人在线| 四虎精品黑人视频| 女人18一级毛片免费观看 | www亚洲精品| 福利一区在线| 国产尤物jk自慰制服喷水| 国产精品蜜芽在线观看| 亚洲人成日本在线观看| 四虎永久在线视频| 日本黄色a视频| 九色在线观看视频| 亚洲国产清纯| 国产a网站| 99久久国产精品无码| 国产在线第二页| 影音先锋亚洲无码| 在线观看国产小视频| 国产福利在线免费| 精品一区二区无码av| 亚洲精品天堂在线观看| 国产福利观看| 国产成人免费高清AⅤ| 亚洲欧美日韩色图| 又大又硬又爽免费视频| 欧美综合中文字幕久久| 国产青榴视频| 欧美亚洲香蕉| 国产精品短篇二区|