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

光束指向裝置高分辨率傳動機構(gòu)

2023-03-10 08:18:28李寶宇陳雨康范大鵬
光學精密工程 2023年3期
關(guān)鍵詞:變形

李寶宇, 陳 寧, 陳雨康, 謝 馨, 范大鵬

(國防科技大學 智能科學學院,湖南 長沙 410073)

1 引 言

光束指向裝置用于光束指向的調(diào)整或者改變成像視軸,實現(xiàn)對目標的捕獲、跟蹤、瞄準和定位,廣泛應用于空間觀測、偵察監(jiān)視、紅外對抗、顯微探測和機器視覺等領(lǐng)域[1]。隨著裝備自動化、智能化程度的提高,對光束指向裝置的指向精度及分辨率、結(jié)構(gòu)緊湊性和輕量化等方面提出了更高的要求。目前,普遍使用的光束指向機構(gòu)存在著機構(gòu)體積大、占用空間大、分辨率不高等問題[2],難以滿足裝備的實際應用需求。

在光束或視軸的精密指向調(diào)整系統(tǒng)中,旋轉(zhuǎn)雙棱鏡系統(tǒng)通過兩棱鏡的共軸獨立旋轉(zhuǎn)改變光的傳播方向,具有結(jié)構(gòu)緊湊,指向精度高,可實現(xiàn)大角度偏轉(zhuǎn)、整體造價低的特點[2]。精密光束指向裝置在光軸的初始校準及零位微調(diào)等應用中,旋轉(zhuǎn)雙棱鏡系統(tǒng)的高分辨率微步進運動可實現(xiàn)光軸的精密微調(diào),與傳統(tǒng)的手動微調(diào)工作臺相比,具有實時自動調(diào)整的優(yōu)勢,同時滿足裝備智能化的發(fā)展趨勢。

扭輪摩擦傳動機構(gòu)利用主動輪與從動輪之間的摩擦來傳遞運動和動力,主要由驅(qū)動軸和扭輪螺母組成,3個均布扭輪通過一定預緊力壓在驅(qū)動軸的光桿上,通過摩擦力使扭輪螺母實現(xiàn)螺旋進給運動[3-9]。該機構(gòu)具有傳動精度及分辨率高、傳動平穩(wěn)、無空程和反向間隙以及結(jié)構(gòu)簡單的特點,可實現(xiàn)傳動平臺的超精密定位,定位分辨率達到納米級水平[4],已在精密機床、超精密工作臺等領(lǐng)域得到了廣泛應用。Mizumoto[4]等設(shè)計了一種采用扭輪摩擦傳動的埃級定位系統(tǒng),實現(xiàn)了0.2 nm的定位分辨率。羅兵[6]深入分析了扭輪摩擦傳動的運動學、動力學問題,并采用扭輪摩擦傳動機構(gòu)設(shè)計了一維工作臺,實測運動分辨率達到10 nm。田軍委[9]等認為機構(gòu)外部載荷很小時可忽略預位移和微觀滑移對機構(gòu)的影響,對扭輪摩擦傳動進行動力學分析,討論了扭輪純滾動的動力學條件,并指出外部載荷對傳動軸切向驅(qū)動力的影響較小。上述研究表明,在外部載荷較小的低速運動場景下,扭輪摩擦傳動特別適合實現(xiàn)超精密定位。

鑒于扭輪摩擦傳動機構(gòu)具有低速高分辨率的運動特性,結(jié)合某型精密光束指向裝置實現(xiàn)光軸精密微調(diào)及空間緊湊的特點,本文提出了一種扭輪偏置的傳動構(gòu)型,即偏置型扭輪摩擦傳動機構(gòu)。與普通扭輪摩擦傳動機構(gòu)不同,該機構(gòu)中扭輪呈圓臺狀結(jié)構(gòu),其回轉(zhuǎn)平面不再通過從動輪的回轉(zhuǎn)中心,而是存在一定偏移量。扭輪偏置最大化地利用了裝置內(nèi)部空間,有效解決了棱鏡透光口徑要求較大而裝置外部尺寸要求較小的矛盾。與齒輪傳動相比,偏置型扭輪摩擦傳動機構(gòu)在有限空間內(nèi)實現(xiàn)了更大的傳動比,具有轉(zhuǎn)動分辨率高、空間緊湊、低成本、無需潤滑和維護簡單的優(yōu)勢。為研究偏置型扭輪摩擦傳動機構(gòu)的回轉(zhuǎn)輸出特性,分析了機構(gòu)傳動機理與運動學特性,并對接觸區(qū)的受力情況和擠壓變形量進行了理論分析;使用RecurDyn多體動力學仿真軟件建立扭輪摩擦傳動的剛?cè)狁詈蟿恿W模型,并對扭角對系統(tǒng)回轉(zhuǎn)位移特性的影響進行仿真分析和物理樣機實驗,為這種傳動機構(gòu)在光束指向裝置傳動系統(tǒng)中應用提供了理論支持。

2 原 理

2.1 機構(gòu)工作原理

為了在有限空間內(nèi)實現(xiàn)更高的傳動精度,本文設(shè)計的偏置型扭輪摩擦傳動機構(gòu)如圖1所示。與傳統(tǒng)的扭輪摩擦傳動直線輸出方式不同,偏置型扭輪摩擦傳動機構(gòu)采用回轉(zhuǎn)輸出形式,以扭輪為主動輪,通過深溝球軸承限制從動輪的軸向平動,使從動輪只沿周向旋轉(zhuǎn)。其中,主動輪與從動輪之間的夾角α定義為扭角,兩輪接觸點形成的偏移量為e。

圖1 偏置型扭輪摩擦傳動機構(gòu)Fig.1 Offset twist-roller frictional transmission mecha?nism

主動輪和從動輪之間受到正壓力作用,使兩摩擦輪在工作時產(chǎn)生微小變形,在接觸部位形成一個微小的接觸區(qū),設(shè)該接觸區(qū)為A,主動輪上的接觸區(qū)為A1,從動輪上的接觸區(qū)為A2。A1中心到主動輪旋轉(zhuǎn)軸線的垂線和A2中心到從動輪旋轉(zhuǎn)軸線的垂線成一定夾角,定義為偏置角β,即有:

因摩擦副材料在受載后發(fā)生微小變形,主動輪和從動輪在接觸區(qū)內(nèi)會產(chǎn)生相對滑動,這種微觀滑移情況是由摩擦副材料彈性造成的,在正常工作時不可避免。常用滑動率[10]定量表征這種微觀滑移,將主、從動輪在接觸點處圓周線速度的降低量與主動輪圓周線速度的比值定義為滑動率ε,即:

其中v1,v2分別為主、從動輪的圓周線速度。特別是對于橡膠/金屬摩擦副材料構(gòu)成的干式摩擦傳動機構(gòu),滑動率可達4%~10%,在計算機構(gòu)傳動比時已不能忽略。

在偏置型扭輪摩擦傳動機構(gòu)中,扭輪的旋轉(zhuǎn)運動帶動從動輪旋轉(zhuǎn),使得扭輪和從動輪之間持續(xù)存在宏觀滑移,為避免摩擦副磨損引起的傳動失效,此類傳動方式適用于低速、小負載工況。主動輪以聚氨酯、聚酰亞胺等耐磨橡膠為主,從動輪常用鋁合金等輕金屬材料制成,主、從動輪之間為干式摩擦傳動。

2.2 運動學分析

在偏置型扭輪摩擦傳動機構(gòu)中,主動輪的旋轉(zhuǎn)軸線與從動輪的旋轉(zhuǎn)軸線為空間異面直線,其夾角為扭角α。同時扭輪偏置安裝,其偏轉(zhuǎn)角為β。偏置型扭輪摩擦傳動機構(gòu)的運動分析如圖2所示。

圖2 偏置型扭輪摩擦傳動機構(gòu)運動學分析Fig.2 Kinematic analysis of offset twist-roller frictional transmission mechanism

主動輪以恒角速度ω1旋轉(zhuǎn)時,將主動輪在接觸區(qū)A1的圓周線速度v1沿水平方向和豎直方向分解為v1h和v1u,則有:

式中r1為接觸區(qū)A1中心所在截面的半徑。將主動輪的豎直分速度v1u沿從動輪轉(zhuǎn)動的切向方向分解為v1t:

從動輪運動的角速度為ω2,則從動輪在接觸區(qū)A2的圓周線速度v2為:

因主、從動輪受載后的微小變形引起微觀滑移,在接觸區(qū)內(nèi)圓周線速度不相等。則偏置型扭輪摩擦傳動機構(gòu)的滑動率εb為:

綜合式(3)~式(6)可得,偏置型扭輪摩擦傳動的傳動比ib為:

2.3 接觸區(qū)擠壓變形量理論分析

扭輪和從動輪在正壓力作用下發(fā)生擠壓變形,接觸區(qū)擠壓變形是扭輪和從動輪間運動傳遞時產(chǎn)生微觀滑移的主要原因。為準確獲得接觸區(qū)擠壓變形量的計算值,應用赫茲接觸理論研究了正壓力與擠壓變形量的關(guān)系。

將扭輪簡化為半徑r1的圓柱體,扭輪和從動輪軸線的夾角為α,在正壓力P的作用下發(fā)生擠壓變形,設(shè)總的擠壓變形量為δ。假設(shè)滿足赫茲接觸理論計算局部變形的簡化條件,則扭輪接觸區(qū)A1形狀為橢圓形斑塊,橢圓斑的尺寸和位置及擠壓變形量δ可通過以下赫茲接觸公式[12]計算:

式中:G=,E1,υ1和E2,υ2分別為扭輪和從動輪的彈性模量和泊松比;a,b分別為橢圓斑的長半軸和短半軸為橢圓斑的離心率;中間參數(shù)m,n,K(e)通過查表確定。

通過赫茲接觸公式計算接觸區(qū)A1橢圓斑的尺寸與位置。如圖3所示,扭輪軸線與Y軸重合,在一定的正壓力作用下,扭輪和從動輪軸線的夾角α從30°到90°的變化過程中,接觸區(qū)A1的橢圓形長軸逐漸變小,短軸逐漸變大,同時橢圓長軸與Y軸的夾角逐漸減小到0,與實際情況相符。

圖3 橢圓斑與扭角的關(guān)系Fig.3 Relationship between elliptical spot and twist angle

擠壓變形量δ的變化曲線如圖4所示,擠壓變形量δ與正壓力P存在非線性正相關(guān)關(guān)系。由式(8)可知,擠壓變形量δ與P2/3成正比關(guān)系。

圖4 擠壓深度與正壓力的關(guān)系Fig.4 Relationship between extrusion depth and contact force

3 扭輪摩擦傳動剛?cè)狁詈蟿恿W模型

接觸是兩個物體在接觸界面上相互作用的復雜力學現(xiàn)象,在固體力學領(lǐng)域中普遍存在。接觸問題的復雜性不僅體現(xiàn)在接觸界面的變形和受力狀況,更主要在于接觸界面的邊界條件非線性,此外,接觸過程還常常涉及材料非線性和幾何非線性。因此,極少數(shù)的接觸問題可以解析處理,絕大多數(shù)只能采用有限元、離散元和邊界元等數(shù)值方法進行模擬,其中有限元法的應用最為廣泛。

RecurDyn是一款基于遞歸算法的多體系統(tǒng)仿真軟件,非常適合求解復雜接觸的多體動力學問題,同時借助其特有的MFBD(Multi Flexible Body Dynamics)多柔體動力學分析技術(shù),完成有限元柔性體和多體動力學的耦合計算,可以更加真實地仿真柔性體的非線性問題。RecurDyn強化了以運動中的接觸、大變形和其他非線性為仿真目的的FFlex有限元柔性體算法,它利用柔性體節(jié)點之間的相對位移和旋轉(zhuǎn)作為節(jié)點坐標描述柔性體的變形,能夠精準預測柔性體和剛性體之間的接觸問題,大幅提升仿真計算的穩(wěn)定性和精確性[11]。

本文采用基于完全柔性體技術(shù)的RecurDyn/FFlex剛?cè)狁詈蟿恿W分析模塊,對扭輪摩擦傳動機構(gòu)進行動力學建模及分析,可較好地反應扭輪的受力變形及變形對傳動精度的影響。建模時,扭輪為基材為鋁合金2A12,覆材為聚酰亞胺(PI)制成的膠輪,從動輪為表面陽極氧化處理的鋁合金2A12。扭輪與從動輪的尺寸及材料參數(shù)如表1所示。

表1 扭輪與從動輪的尺寸及材料參數(shù)Tab.1 Size and material parameters of transmission mechanism

因扭輪的彈性模量遠低于從動輪,將扭輪看作柔性體,從動輪視為剛性體,建立扭輪摩擦傳動機構(gòu)的剛?cè)狁詈夏P停阋哉鎸嵎从衬Σ硫?qū)動的實際情況,同時簡化仿真計算量。

3.1 剛性體模型

在CAD軟件中建立扭輪摩擦傳動機構(gòu)模型,導入到RecurDyn動力學仿真軟件中,定義重力加速度方向沿-Z方向。首先,在扭輪摩擦傳動機構(gòu)的多剛體模型中,分別定義電機和支架之間固定鉸連接、電機及支架相對于大地在Z方向上的平移副、從動輪相對于大地的旋轉(zhuǎn)副、扭輪和相對于電機及支架的旋轉(zhuǎn)副;其次,為扭輪的旋轉(zhuǎn)運動添加驅(qū)動,通過函數(shù)表達式定義扭輪的轉(zhuǎn)動是與時間相關(guān)性的關(guān)于位移、速度、加速度的函數(shù),為電機及支架的平移運動施加載荷,定義大地對支架沿+Z方向的彈簧力,扭輪和從動輪在彈簧力的作用下相互接觸;最后,創(chuàng)建扭輪和從動輪兩剛性體之間的面接觸,方便后續(xù)柔性接觸建模時識別接觸面。對建立好的剛性體模型進行多剛體動力學分析,以驗證仿真模型的正確性。

3.2 扭輪有限元柔性體模型

利用RecurDyn/Mesher將多剛體仿真模型中的扭輪進行網(wǎng)格劃分,選擇Solid4實體單元并設(shè)置單元的材料屬性,劃分網(wǎng)格后選擇扭輪接觸面上的所有節(jié)點,創(chuàng)建輸出節(jié)點集。

退出Mesher建模環(huán)境后,扭輪成為有限元柔性體,同時在劃分網(wǎng)格后能夠保持原模型中的裝配關(guān)系,并在旋轉(zhuǎn)副處生成一個剛化區(qū)域,在已設(shè)定接觸的位置產(chǎn)生創(chuàng)建柔性接觸所必需的片集。偏置型扭輪摩擦傳動機構(gòu)的剛?cè)狁詈夏P腿鐖D5所示。

圖5 扭輪摩擦傳動剛?cè)狁詈夏P虵ig.5 Rigid flexible coupling model of twist-roller fric?tional transmission mechanism

3.3 有限元柔性體接觸模型

扭輪摩擦傳動機構(gòu)中,從動輪旋轉(zhuǎn)運動的動力由扭輪對從動輪摩擦力的切向分力提供。因此,對有限元柔性體和剛體的接觸建模,是扭輪摩擦傳動剛?cè)狁詈戏抡婺P椭凶铌P(guān)鍵的部分,而在柔性體和剛性體發(fā)生接觸時,接觸力計算是建模時最為重要的一個問題。

創(chuàng)建扭輪柔性面和從動輪剛性面的接觸,扭輪柔性面設(shè)置為作用面,根據(jù)網(wǎng)格的形狀用三角形或者四邊形的片集替代,從動輪剛性面設(shè)置為基礎(chǔ)面,程序會根據(jù)幾何形狀將作用面和基礎(chǔ)面離散為三角形碎片,用于接觸分析。為改善三角形曲面的逼真程度,通過降低平面容許誤差系數(shù)和平面最小容許誤差來解決。

扭輪柔性面和從動輪剛性面在上述彈簧力的作用下擠壓變形,RecurDyn通過接觸穿透深度計算法向接觸力fn,其計算公式為:

其中:k為剛度接觸系數(shù),c為阻尼系數(shù),δ為接觸穿透深度為接觸點的相對速度,m1,m2,m3分別為剛度指數(shù)、阻尼指數(shù)和凹痕指數(shù)。

在獲得接觸力后,設(shè)置動、靜摩擦系數(shù)及各自發(fā)生作用時的相應速度,使用滑動摩擦模型計算接觸摩擦力。為避免速度換向過程中的摩擦力瞬變,在正負靜滑動速度間的靜摩擦系數(shù)通過連續(xù)插值過渡[11]。

4 扭角對輸出特性的影響

在扭輪摩擦傳動中,從動輪旋轉(zhuǎn)的驅(qū)動力為它受到摩擦力的余弦分量,因此扭輪和從動輪之間的扭角直接影響機構(gòu)的傳動性能。當扭角為0°時,即為圓柱平摩擦輪傳動,此時從動輪受到的摩擦力最大;扭角逐漸增大時,從動輪的驅(qū)動力矩逐漸降低,當從動輪的驅(qū)動力矩不足以克服阻力矩時,從動輪停止轉(zhuǎn)動,此時扭輪和從動輪之間產(chǎn)生持續(xù)滑移,扭輪磨損量急劇增大。為避免出現(xiàn)此類傳動失效現(xiàn)象,需要對扭輪和從動輪間的扭角與機構(gòu)輸出特性的關(guān)系進行研究。

搭建如圖6所示的實驗系統(tǒng)。其中,從動機構(gòu)固定在Z向平移臺上,和Z向平移臺一起固定在X向平移臺上;驅(qū)動機構(gòu)、三向力傳感器和手動旋轉(zhuǎn)臺依次固定連接,三者一起固定在Y向平移臺上。三向平移臺用于確定扭輪和從動輪摩擦副的接觸位置,即偏移量e,Z向平移臺和三向力傳感器配合,確定并調(diào)整扭輪和從動輪之間的正壓力N,手動旋轉(zhuǎn)臺用于確定和調(diào)整扭角α。

圖6 物理樣機實驗系統(tǒng)Fig.6 Experimental system of prototype

電機控制和信號采集系統(tǒng)采用dSPACE半實物仿真系統(tǒng),實驗系統(tǒng)實物如圖7所示,實驗儀器與設(shè)備如表2所示。

圖7 扭角與傳動機構(gòu)輸出特性分析實驗系統(tǒng)Fig.7 Experiment system for relation analysis between twist angle and output characteristics of transmis?sion system

表2 實驗儀器與設(shè)備Tab.2 Experimental instruments and equipment

4.1 扭角對輸出角速度穩(wěn)定性的影響

從動輪轉(zhuǎn)動時存在摩擦力矩作用,當扭角達到一定值時,扭輪提供給從動輪的驅(qū)動力矩不足以克服摩擦力矩,此時傳動失效。經(jīng)多次實驗測試,現(xiàn)有實驗條件下扭角α的極限值為85°。因此,在扭角α為80°,82°和84°三種條件下,通過仿真分析和實驗測試,研究了扭角對機構(gòu)的輸出速度穩(wěn)定性、傳動比和回轉(zhuǎn)分辨率的影響。

扭輪以恒角速度運轉(zhuǎn)時的仿真分析和測試結(jié)果如圖8所示。可以看出,隨著扭角的增大,從動輪運轉(zhuǎn)一周的時間明顯延長,仿真結(jié)果和實驗結(jié)果吻合度較好。

圖8 扭角與輸出角位置的關(guān)系Fig.8 Relationship between twist angle and angle posi?tion

從動輪回轉(zhuǎn)角速度的對比結(jié)果如圖9所示。從動輪運轉(zhuǎn)一周的過程中,不同扭角條件下的從動輪角速度仿真結(jié)果和實驗結(jié)果出現(xiàn)不同程度的波動,在實驗結(jié)果中,角速度波動范圍較大,且隨著扭角的增大,角速度的波動程度越大,在20°~80°位置時,不同扭角條件下的角速度出現(xiàn)相似的波動情況,與從動輪絕對角位置明顯相關(guān)。

圖9 扭角對輸出角速度的影響Fig.9 Effect of twist angle on angular velocity

根據(jù)文獻[12]可知,在扭輪和從動輪的接觸過程中,扭輪接觸區(qū)的基體與從動輪金屬產(chǎn)生黏著,在驅(qū)動力的作用下,黏著結(jié)點被剪切而產(chǎn)生滑動,同時產(chǎn)生新的黏著結(jié)點,扭輪和從動輪間的滑動摩擦就是黏著結(jié)點的形成和剪切交替發(fā)生的過程。因此在仿真分析結(jié)果中,扭輪的受壓變形是從動輪角速度出現(xiàn)微小波動的主要原因。實際測試時,受扭輪和從動輪的加工精度和表面質(zhì)量的影響,摩擦條件發(fā)生改變,使得扭輪和從動輪之間的摩擦力發(fā)生變化。從動輪在變化的驅(qū)動力作用下做加速、減速轉(zhuǎn)動,角速度表現(xiàn)出波動性,且在扭角較大時,從動輪做加、減速轉(zhuǎn)動的時間越長,角速度波動越劇烈。

4.2 扭角對輸出角位置分辨率的影響

精密光束指向裝置中,旋轉(zhuǎn)分辨率是一個重要的性能指標,直接影響著光束的指向精度。為研究扭角對輸出角位置分辨率的影響,根據(jù)現(xiàn)有實驗條件,分別在扭角α為80°,82°和84°三種條件下,對步進電機輸入固定脈沖數(shù)的階梯信號,階梯信號間隔1 s,采集從動輪位置編碼器測量值,測試結(jié)果如圖10所示。

圖10 輸出角位置與脈沖數(shù)的關(guān)系Fig.10 Relationship between angle position and pulse number

由圖10可知,從動輪角位置輸出信號同樣為階梯形狀,扭角相同時,輸入信號脈沖數(shù)越大,經(jīng)歷相同時間后從動輪輸出的絕對角位置越大;輸入信號脈沖數(shù)相同的情況下,扭角越大,經(jīng)歷相同時間后從動輪輸出的絕對角位置越小。

根據(jù)上述實驗結(jié)果,扭角α為80°時,取輸入信號脈沖數(shù)為15,扭角α為82°和84°時,取輸入信號脈沖數(shù)為20,對偏置型扭輪摩擦傳動機構(gòu)在階梯信號下的響應進行仿真分析,并計算從動輪輸出角位置的分辨率,實驗和仿真結(jié)果對比如圖11所示。

圖11 角位置分辨率實驗和仿真結(jié)果Fig. 11 Experimental and simulation results of angular resolution

實驗和仿真結(jié)果表明,扭角越大,機構(gòu)輸出角位置分辨率越小,三種扭角條件下的機構(gòu)分辨率基本低于0.009°;在扭角α為84°時,偏置型扭輪摩擦傳動機構(gòu)能實現(xiàn)0.005°(18″)的低速微步進運動。而同樣的安裝空間內(nèi)采用單級圓柱齒輪傳動,最大傳動比為5.5,理論上最高分辨率為0.022°,遠大于偏置型扭輪摩擦傳動機構(gòu)的輸出分辨率。

5 結(jié) 論

本文針對扭角對偏置型摩擦傳動機構(gòu)輸出特性的影響,通過物理樣機實驗和多體動力學仿真進行分析。結(jié)果表明:扭角越大,機構(gòu)傳動比越大,但從動輪輸出角速度波動程度越大,傳動越不平穩(wěn);在某型光束指向裝置中,使用微型減速步進電機驅(qū)動、偏置型扭輪摩擦傳動方案巧妙地解決了裝置透光口徑與外部尺寸的矛盾,在高分辨率和低成本的要求下,實現(xiàn)了0.005?(18″)的低速微步進運動,優(yōu)于單級圓柱齒輪的傳動方案。

本文充分利用低負載下扭輪摩擦傳動高分辨率的優(yōu)點,在空間狹小的前提下設(shè)計了一種偏置型扭輪摩擦傳動機構(gòu),解決了傳統(tǒng)驅(qū)動傳動方案無法滿足某型光束指向裝置高分辨率和空間緊湊的要求,對于高分辨率傳動方式的構(gòu)型設(shè)計與實際應用具有指導意義。

猜你喜歡
變形
變形記
談詩的變形
中華詩詞(2020年1期)2020-09-21 09:24:52
柯西不等式的變形及應用
“變形記”教你變形
不會變形的云
“我”的變形計
會變形的折紙
童話世界(2018年14期)2018-05-29 00:48:08
變形巧算
例談拼圖與整式變形
會變形的餅
主站蜘蛛池模板: 亚洲三级影院| 狠狠色丁香婷婷| 国产亚洲视频免费播放| 欧美在线综合视频| 久久久久青草大香线综合精品| 亚洲欧美精品日韩欧美| 中文字幕首页系列人妻| 91精品国产一区自在线拍| 免费国产一级 片内射老| 亚洲Aⅴ无码专区在线观看q| 国产精品免费福利久久播放| 亚洲三级网站| 午夜不卡视频| 成年人久久黄色网站| 亚洲欧美日韩精品专区| 欧美在线一二区| 无码av免费不卡在线观看| 国产精品手机在线播放| 日韩黄色大片免费看| 特级毛片免费视频| 国产成人乱无码视频| 秘书高跟黑色丝袜国产91在线| 亚洲欧美在线综合图区| 久久大香伊蕉在人线观看热2| 日韩性网站| 在线精品欧美日韩| 性视频久久| 亚洲91精品视频| 欧美成人手机在线观看网址| 亚洲永久色| 一边摸一边做爽的视频17国产| 日韩欧美高清视频| 日韩中文欧美| 国产精品9| 中文成人在线视频| 久久久久无码精品国产免费| 亚洲欧美极品| 亚洲天堂区| 国产成人亚洲精品蜜芽影院| 免费观看成人久久网免费观看| 国产噜噜噜视频在线观看| 久久99精品久久久大学生| www.狠狠| 呦系列视频一区二区三区| 亚洲成人在线网| 91久久精品国产| 国产精品网址在线观看你懂的| 69国产精品视频免费| 中文无码日韩精品| 手机在线免费不卡一区二| 美女被操黄色视频网站| 国产免费高清无需播放器| 亚洲av成人无码网站在线观看| 久久中文电影| 日韩不卡高清视频| 久久青青草原亚洲av无码| 免费看的一级毛片| 免费无码一区二区| 欧美中文字幕无线码视频| av大片在线无码免费| 日韩AV手机在线观看蜜芽| 天天综合天天综合| 日韩天堂网| 国产麻豆永久视频| 日韩在线播放欧美字幕| 9966国产精品视频| 日韩av电影一区二区三区四区| 日本亚洲最大的色成网站www| 免费看美女毛片| 国产在线专区| 一级毛片视频免费| 久久精品国产国语对白| www.亚洲国产| 欧美日本激情| 国产二级毛片| 91香蕉视频下载网站| 综合久久久久久久综合网| 欧美成人h精品网站| 中文字幕在线免费看| 成年女人a毛片免费视频| 青青草原国产| 亚洲国产中文综合专区在|