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

        ?

        磨床用自定心中心架的機構(gòu)設(shè)計與優(yōu)化

        2016-10-10 05:05:12趙永強王耀鋒侯紅玲
        中國機械工程 2016年17期
        關(guān)鍵詞:定心滾輪凸輪

        趙永強 王耀鋒 侯紅玲

        陜西理工大學(xué),漢中,723001

        ?

        磨床用自定心中心架的機構(gòu)設(shè)計與優(yōu)化

        趙永強王耀鋒侯紅玲

        陜西理工大學(xué),漢中,723001

        設(shè)計了一種磨床用自定心中心架的原型機構(gòu),建立了自定心中心架的參數(shù)組合和平動凸輪的參數(shù)化模型。通過對機構(gòu)的設(shè)計參數(shù)組合篩選,提出了機構(gòu)的優(yōu)化設(shè)計目標和優(yōu)化函數(shù),設(shè)計了優(yōu)化的計算程序。在保證中心架的傳力性能、夾持穩(wěn)定性和夾持精度的前提下,對機構(gòu)中平動凸輪的兩個重要參數(shù)進行了優(yōu)化,為磨床用高精度自定心中心架的設(shè)計提供了完整的設(shè)計方案。

        自定心中心架;細長軸磨削;平動凸輪;參數(shù)化設(shè)計

        0 引言

        細長軸類零件的加工一直以來都是機械工程師的關(guān)注焦點,為了提高軸類零件的加工精度,各種不同的加工工藝被引入實際加工中,其中,采用中心架作為輔助支承提高零件的剛度是一種理想的工藝方法。尤其是在超長絲杠和螺桿的磨削過程中,中心架的定心定位精度成為提高零件加工精度的制約因素[1]。

        國內(nèi)對中心架的研究大多停留在傳統(tǒng)的設(shè)計基礎(chǔ)之上,如彭克立等[2]設(shè)計了一種高速、超高速磨削的自動跟蹤式中心架;孫金海等[3]設(shè)計了一種采用V形塊自動定心的絲杠自鎖夾緊的自動定心中心架;王偉順等[4]采用靜壓原理設(shè)計了一種重型機床加工大型軸類零件用的高精度重載靜壓中心架;牛福林等[5]設(shè)計了一種采用V形滾輪式夾緊機構(gòu)、T形絲杠自鎖和液壓剎緊機構(gòu)的手動高剛性自定心中心架。以上的研究都存在中心架自動化程度不高、夾緊力調(diào)節(jié)繁瑣、適應(yīng)性較差等缺點。

        早在1960年前后,德國某夾具廠商率先推出了一種自定心中心架。隨后瑞典、美國等歐美國家也研制出了多種自定心中心架。此類中心架利用凸輪無窮多精確點的特性,在定位精度和夾持可靠性等方面具有突出優(yōu)勢,因而該類型自定心中心架迄今仍是市場上的主流產(chǎn)品[6-12]。

        近些年來,國內(nèi)針對自定心中心架的研究也取得了一些成果。李京魯[13]對自定心中心架的楔形凸輪的輪廓曲線進行了理論推導(dǎo);涂庭波[14]對液壓自定心中心架的凸輪機構(gòu)進行了研究,給出了凸輪設(shè)計的方法;盧學(xué)玉等[15]對液壓自定心中心架凸輪曲線進行了理論推導(dǎo);肖凱[16]從加工的角度給出了凸輪曲線的求解方法;Wu等[17]和Wang等[18]從液壓角度分析了影響中心架的工作壓力的因素,并對參數(shù)進行了優(yōu)化。由此可見:國內(nèi)學(xué)者針對自定心中心架的研究主要是圍繞凸輪機構(gòu)的理論計算和液壓夾緊力的調(diào)節(jié)兩個方面進行的,而這兩個因素正是保證夾持中心恒定和夾持力穩(wěn)定的兩個關(guān)鍵點,其中凸輪輪廓的精確設(shè)計與加工是實現(xiàn)夾持中心恒定的前提。

        以上對凸輪輪廓的設(shè)計多采用包絡(luò)法、瞬心法、復(fù)變量法等。包絡(luò)法利用凸輪和從動件的幾何關(guān)系推導(dǎo)出接觸點的軌跡方程;瞬心法可同時確定壓力角、平面凸輪輪廓和凸輪曲率半徑;復(fù)變量法可同時求得壓力角和曲率半徑,但包絡(luò)法和瞬心法只是近似求解方法,求解過程繁瑣,精度不高;復(fù)變量法的求解精度高,但是求解方程式比較復(fù)雜,關(guān)系推導(dǎo)困難[19-20]。

        近年來,基于參數(shù)化數(shù)值求解的解析法在復(fù)雜輪廓凸輪設(shè)計中的優(yōu)越性越來越明顯,多軸數(shù)控銑削加工和CNC成形磨削技術(shù)的發(fā)展,使復(fù)雜輪廓凸輪的設(shè)計與加工逐步成熟。筆者設(shè)計了一種磨床用自定心中心架,建立了關(guān)鍵零件平動凸輪的數(shù)值計算模型,利用參數(shù)化設(shè)計手段,將平動凸輪與滾輪從動件的壓力角、夾持力的周向間隔角度以及擺桿桿長作為約束條件,對自定心中心架機構(gòu)進行優(yōu)化設(shè)計??紤]自定心中心架系對稱結(jié)構(gòu),故文中涉及的計算及優(yōu)化設(shè)計除特別說明外,均只針對機構(gòu)中的上半部分進行。

        1 自定心中心架的工作原理

        自定心中心架機構(gòu)原理如圖1所示。具有對稱工作曲線的平動凸輪做左右平移,滾輪1、2、3的結(jié)構(gòu)尺寸相同,滾輪3的中心D固定在平動凸輪的右端,工件的夾持中心始終在O3點。中心架夾緊工件的過程為:在外部驅(qū)動力作用下,平動凸輪及滾輪3向右移動,與凸輪接觸的兩個滾輪4、5受平動凸輪的作用,帶動擺桿1、2繞著各自的鉸點O1、O2轉(zhuǎn)動,從而帶動滾輪1和滾輪2轉(zhuǎn)動,實現(xiàn)對工件的夾緊。在夾持過程中,滾輪1、2、3始終同步運動,同時夾緊或松開工件,從而使被加工工件軸心始終保持在O3位置,以實現(xiàn)中心架對工件的夾持和定心。

        圖1 自定心中心架的機構(gòu)原理

        在自定心中心架設(shè)計中,為了推算平動凸輪的廓形,首先保持機架和被加工工件軸心O3不變,初步設(shè)計擺桿1、2的長度,通過改變被夾持工件的半徑使擺桿1、2轉(zhuǎn)動,滾輪4、5的中心擺動和滾輪3的移動形成平動凸輪的廓形。中心架的外加驅(qū)動力為F0,滾輪1、2、3對工件的夾持力依次為F4、F5、F3。

        2 平動凸輪廓形的解析計算

        2.1坐標系的建立

        自定心中心架的參數(shù)定義如圖2所示。以擺桿1的軸心O1為原點,以平動凸輪的平移軸為X1軸建立靜坐標系O1X1Y1。以平動凸輪的左端面為OYZ平面,以過被夾持工件中心O3的直線為X軸,建立動坐標系OXY,O1X1Y1與OXY之間的Y向距離為L1cosθ,X向距離為L3,則

        (1)

        L4=L0cos(α+β)

        其中,L5為平動凸輪的預(yù)設(shè)值;Lu為平動凸輪沿X方向的移動距離。

        圖2 自定心中心架的參數(shù)定義

        2.2平動凸輪廓形計算

        設(shè)定擺桿的左擺臂O1A的長度為L0,右擺臂O1C的長度為L2,鉸點O1與工件軸心O3的距離O1O3為L1,O1O3與右擺臂O1C的夾角∠O3O1C=α;左右擺臂之間的夾角∠AO1C=φ;滾輪4、5的半徑為r,r0為滾輪1、2、3的半徑;R為被夾持工件的半徑,被夾持工件半徑R的最大值記為最大夾持半徑Rmax,最小值記為最小夾持半徑Rmin,則R∈[Rmin,Rmax];Rmax對應(yīng)的右擺臂O1C與X軸之間的夾角記為ε。

        最大夾持半徑Rmax和最小夾持半徑Rmin對應(yīng)的右擺臂O1C與O1O3之間的夾角分別記為αmax和αmin,則

        (2)

        (3)

        滾輪1的中心C與工件中心O3之間的距離

        (4)

        工件中心O3與擺桿1、2的鉸點O1、O2的連線之間的夾角∠O3O1O2=θ,其值為

        (5)

        左擺臂O1A在最高位置時與X軸之間的夾角記為β,其值為

        (6)

        如圖3所示,被夾持工件的半徑從Rmin變化至Rmax時,滾輪4的中心A從A1移動至A2,同時平動凸輪向左平移的距離

        Lu=Lx-(L1-L2)

        (7)

        圖3 凸輪上A點坐標計算

        當(dāng)被夾持工件的半徑R最大時,β=βmax,滾輪4的中心A處于A1點,可推算

        L4=L0cosβmax

        (8)

        在O1X1Y1坐標系中,滾輪4的中心點A的坐標xA、yA為

        (9)

        根據(jù)平動凸輪的理論廓線,基于包絡(luò)法可求解得到平動凸輪的實際廓形坐標:

        (10)

        y1=L0cos(α+β)

        將x0、y0從坐標系O1X1Y1變換到坐標系OXY中即可得到凸輪的最終廓線方程:

        (11)

        3 主要參數(shù)計算

        在自定心中心架中,初始參數(shù)選擇不當(dāng)會導(dǎo)致凸輪輪廓設(shè)計不合理,機構(gòu)不能實現(xiàn)預(yù)期的功能,或者運行中的動力學(xué)特性欠佳。根據(jù)平動凸輪廓形計算過程推導(dǎo)和前期的預(yù)驗算,選擇自定心中心架中平動凸輪的壓力角和影響夾持質(zhì)量的夾持力間隔角度兩個參數(shù)分別進行分析計算[21-22]。

        3.1平動凸輪的壓力角

        如圖1和圖2所示,設(shè)α1為凸輪實際廓線的法線n-n與水平線之間的夾角,則

        (12)

        設(shè)β1為凸輪對滾輪4的力作用線n-n與左擺臂之間的夾角:

        β1=α1-π/2-α+ψ-θ

        (13)

        滾輪1對工件的壓力作用線與水平軸線之間的夾角

        (14)

        滾輪1對工件的力作用線與右擺臂之間的夾角

        (15)

        平動凸輪的壓力角

        α1=π/2-β1

        (16)

        當(dāng)被夾持工件的半徑變化時,可以得出平動凸輪壓力角的變化趨勢,如圖4所示。從圖4可以看出,平動凸輪的壓力角隨被夾持工件半徑的增大而增大,近似成線性比例關(guān)系。

        圖4 工件半徑與平動凸輪壓力角的關(guān)系

        3.2夾持力周向間隔角

        夾持力周向間隔角定義為擺動夾持滾輪1(或2)對工件的夾持力F4(或F5)與滾輪3對工件的夾持力F3之間的夾角,其值為

        φ=π-δ1

        (17)

        在自定心中心架夾持工件時,滾輪1、2、3的三個中心C、E、D,滾輪與工件的切點和工件的夾持中心O3之間都是相互共線的,并且與各夾持力的作用線相重合。改變工件半徑之后夾持力周向間隔角變化如圖5所示,三個夾持力隨被夾持工件半徑和夾持力軸向間隔角的變化如圖6和圖7所示。在圖7中,工件未受切削力作用,φ越趨近120°,三個夾持滾輪對工件的夾持力在數(shù)值上越接近,中心架的夾持越穩(wěn)定,夾持質(zhì)量越好。

        圖5 工件半徑與夾持力周向間隔角的關(guān)系

        圖6 三個夾持力與工件半徑的關(guān)系

        圖7 三個夾持力與周向間隔角的關(guān)系

        由圖6、圖7可見:在當(dāng)前尺寸組合條件及工況下,隨著被夾持工件半徑的增大,三個夾持力的差值先減小后增大,即夾持質(zhì)量先提高后降低。圖6、圖7中曲線的交點代表該中心架夾持質(zhì)量最理想、最穩(wěn)定的夾持情況。

        4 平動凸輪廓形的優(yōu)化設(shè)計

        (1)設(shè)計變量。該中心架設(shè)計的夾持工件的半徑R的范圍為20~40 mm,預(yù)取L1=81 mm,r=5 mm,ψ=108°,θ=45°,L0、L2為獨立變量,則設(shè)計變量為

        x=(x1,x2)T=(L0,L2)T

        (18)

        (2)目標函數(shù)。為了適應(yīng)夾持不同半徑工件的情況,隨著夾持工件半徑R的變化,平動凸輪壓力角α與夾持力周向間隔角φ都會變化,α的變化會使驅(qū)動力F0顯著變化,φ值的變化是中心架穩(wěn)定的一個重要保證,因此,將α和φ作為目標函數(shù)對于中心架的夾持質(zhì)量和穩(wěn)定性具有重要意義。目標函數(shù)為

        (19)

        (20)

        目標函數(shù)式(19)和式(20)是兩個獨立的函數(shù),根據(jù)約束條件與前述的參數(shù)計算方法,采用一維搜索法可以進行求解優(yōu)化。

        (3)約束條件。各擺桿長度的臨界值設(shè)計為

        s.t.L0-L2≤0

        50mm≤L0≤70mm

        50mm≤L2≤70mm

        (4)優(yōu)化程序設(shè)計。利用MATLAB優(yōu)化求解,首先在變量的可變范圍內(nèi)尋求目標的最優(yōu)值,在最優(yōu)值處檢驗條件是否滿足,再根據(jù)最優(yōu)值求出各變量的最優(yōu)組合,具體方法如圖8所示。

        圖8 優(yōu)化設(shè)計流程

        在給L1賦值時,為了使平動推桿的長度合理,需要對θ進行限制,其初始取值范圍限定在40°~50°之間,具體取值還需考慮其他設(shè)計因素。

        具體算法如下:

        clc;clearall;

        symsL2%將L2作為未知參變量

        r=7;%夾持滾輪的半徑

        Rmax=45+r;%被夾持工件半徑最大時,被夾持工件中心與滾輪中心的距離

        L1=81; %假想的固定桿桿長的初設(shè)值L1

        epsi=62*pi/180;%被夾持工件半徑最大時,上夾持滾輪中心、夾持中心點連線與水平軸線之間的夾角

        format long

        alpha=acos((L1^2+L2^2-Rmax^2)/2/L1/L2);

        %夾持不同半徑工件對應(yīng)的夾持角度

        delta=asin(L2*sin(alpha)/Rmax);

        %夾持最大工件時對應(yīng)的∠O1OC角

        theta=delta+epsi-pi/2;

        %假想的L1桿與豎直軸線的夾角

        theta=subs(theta);

        L2=solve(theta-pi/4,L2);

        %預(yù)設(shè)給定的θ為π/4

        L2=subs(L2);

        fprintf('L2=%f ',L2);

        自定心中心架的設(shè)計流程如圖8所示。腳本1、腳本2也可以合并為一個腳本文件(此處為了流程圖的明晰,未加以合并),腳本文件2也包含有靜力學(xué)、動力學(xué)分析等反映動態(tài)性能的數(shù)據(jù)輸出;腳本文件3用作凸輪廓線的圖形輸出及點陣(點云)輸出。對于流程中涉及的腳本文件代碼,限于篇幅此處不作詳述。

        由圖9可知,優(yōu)化后不論是壓力角還是夾持力周向間隔角度都得到了很好的控制:壓力角的最小值較優(yōu)化前減小了50%,壓力角的最大值較優(yōu)化前也減小了50%;夾持力周向間隔角度較優(yōu)化前分布區(qū)間也更加集中,可以很好地滿足機構(gòu)對傳力性能及穩(wěn)定性的要求。

        圖9 優(yōu)化后平動凸輪壓力角、夾持力周向間隔角與工件半徑的關(guān)系

        5 結(jié)論

        本文給出了一種磨床用自定心中心架的機構(gòu)原型,建立了自定心中心架的參數(shù)組合和平動凸輪的參數(shù)化模型。對機構(gòu)中的平動凸輪的兩個重要參數(shù)進行優(yōu)化。設(shè)計計算結(jié)果表明:

        (1)平動凸輪的輪廓曲線滿足了自定心中心架所要求的復(fù)雜運動規(guī)律;

        (2)建立的平動凸輪廓形優(yōu)化設(shè)計數(shù)學(xué)模型合理,通過優(yōu)化設(shè)計使得自定心中心架的工作壓力角的最小值減小了50%,壓力角最大值也減小了50%,夾持力周向間隔角度的分布區(qū)間得到了改善,實現(xiàn)了傳力性能與夾持質(zhì)量的綜合優(yōu)化。

        本研究的設(shè)計方法和思路可以為同類機構(gòu)的設(shè)計提供參考。

        [1]趙永強,王耀鋒.自定心中心架的力學(xué)分析與仿真[J].制造業(yè)自動化,2015,37(7):99-102.Zhao Yongqiang,Wang Yaofeng. Mechanical Analysis and Simulation of Self-centering Center Rest [J]. Manufacturing Automation,2015,37(7):99-102.

        [2]彭克立, 崔中, 劉運泉. 一種自動跟蹤式中心架的研制[J]. 制造技術(shù)與機床, 2009(11):123-125.

        Peng Keli,Cui Zhong,Liu Yunquan.Design of an Auto-tracking Central Shelves[J]. Manufacturing Technology & Machine Tool,2009(11):123-125.

        [3]孫金海,楊曉磊. 自定心、滾動支撐中心架的設(shè)計[J]. 制造技術(shù)與機床,2013(8):145-146.

        Sun Jinhai,Yang Xiaolei. Design of Self- centering and Rolling Supporting Steady Rest[J]. Manufacturing Technology & Machine Tool,2013(8):145-146.

        [4]王偉順,運同樹. 重載機床靜壓刀架的設(shè)計與研究[J]. 金屬加工(冷加工), 2014(14):45-46.

        Wang Weishun,Yun Tongshu.The Design and Research of Heavy-duty Machine Tool Hydrostatic Carriage[J].Metal Working(Metal Cutting),2014(14):45-46.

        [5] 牛福林. 手動高剛性自定心中心架[J]. 制造技術(shù)與機床,2015(10):28-29.

        Niu Fuxin. Handing Good Rigidity Self-centring Centre Frame[J]. Manufacturing Technology & Machine Tool,2015(10):28-29.

        [6]魏嵩.凸輪連桿自定心中心架的設(shè)計[D].新竹:清華大學(xué),2014.

        [7]Takano Y, Kimura Y, Misu A. Center Frame in Crawl-er Travel Device and Method for Manufacturing Center Frame: US,8066302[P]. 2011-11-29.

        [8]Rehm F. Steady Rest:US,8226992B2[P]. 2012-09-18.

        [9]Huang H C, Wei S J, Su H C. Auto-centering Appar-atus of a Steady Rest: US, 14/091066[P]. 2013-11-26.

        [10]Helm P, Maurer E. Self-centering Steady Rest: US, 8726772[P]. 2014-05-20.

        [11]Tsai P T H. Center Frame of CNC Machine Tool: US,14/272564[P]. 2014-05-08.

        [12]Hangleiter E,Schenk P. Electric-drive Steady Rest:US, 8813616B2[P].2014-08-26.

        [13]李京魯.自定心中心架工作楔形曲線的計算[J].機械工程師,1991(11):27-28.

        Li Jinglu. Calculation of Wedge-shape Curces in Self-centering Steady Rest[J]. Machinery Manufacturing Engineer, 1991(11):27-28.

        [14]涂庭波. 液壓自定心中心架的機構(gòu)設(shè)計研究[D]. 合肥:合肥工業(yè)大學(xué), 2008.

        [15]盧學(xué)玉,謝能斌. 液壓自定心中心架凸輪曲線的推導(dǎo)研究[J].機械設(shè)計與制造, 2008(4):24-25.

        Lu Xueyu,Xie Nengbin.The Cam Curve Deduce Research of the Hydraulic Self-centred Steady Rest[J]. Machinery Design & Manufacture, 2008(4):24-25.

        [16]肖凱. 自定心中心架內(nèi)部進程凸輪曲線的計算[J]. 裝備制造技術(shù), 2013(3):149-150.

        XiaoKai.Self-centeringCenterRackInternalProcessofCalculationoftheCamCurve[J].EquipmentManufacturingTechnology, 2013(3):149-150.

        [17]WuB,YuXD,ChangXM,etal.InfluenceofWorkingParametersonDynamicPressureEffectofHeavyConstantFlowHydrostaticCenterRest[J].AppliedMechanicsandMaterials, 2013, 274: 82-86.

        [18]WangWS,DongXZ,WuY,etal.OptimalDesignofHeavyLoadHydrostaticCenterRest[J].AppliedMechanicsandMaterials, 2014, 494: 611-615.

        [19]石永剛. 凸輪機構(gòu)設(shè)計[M].上海: 上??茖W(xué)技術(shù)出版社, 1995.

        [20]黃尚兵, 金光,安源,等.滾子擺動從動件盤形凸輪的設(shè)計及優(yōu)化[J]. 工程設(shè)計學(xué)報,2012,19(6):449-453.

        HuangShangbing,JinGuang,AnYuan,etal.DesignandOptimizationofOscillatingRollerFollowerCam[J].ChineseJournalofEngineeringDesign, 2012,19(6):449-453.

        [21]WangZ,DuanR,XuX.ModelIdentificationofHydrostaticCenterFrameControlSystembasedonMATLAB[J].JournalofNetworks, 2013, 8(6): 1315-1321.

        [22]李軍.基于MATLAB的平面盤形凸輪機構(gòu)參數(shù)化設(shè)計[J]. 煤炭技術(shù),2011,30(3):22-24.

        LiJun.ParametricDesignofPlaneDiscMechanismBasedonMATLAB[J].CoalTechnology,2011,30(3):22-24.

        (編輯盧湘帆)

        MechanismDesignandOptimizationofSelf-centeringSteadyRestUsedinGrindingMachine

        ZhaoYongqiangWangYaofengHouHongling

        ShaanxiSci-TechUniversity,Hanzhong,Shaanxi,723001

        Aprototypeofself-centeringsteadyrestforgrindingmachinewasdesigned;theparametercombinationsofsteadyrestandparametricmodeloftranslationalcamwereestablished.Optimizationdesigngoalandoptimizationfunctionwerepresentedbycombiningscreeningofdesignparameters.TheoptimizationcalculationprograminMATLABwasaccomplished.Twoparametersofthetranslationalcamwereoptimizedbasedonguaranteeingtransmissionperformance,clampingstabilityandespeciallyclampingaccuracy.So,acompletesolutionforhighaccuracy-centeringsteadyrestdesignwasprovided.

        self-centeringsteadyrest;slenderaxisgrinding;translationalcam;parameteriza-tiondesign

        2015-10-19

        陜西省科技統(tǒng)籌創(chuàng)新工程重點實驗室項目(2014SZS16-K03);陜西省教育廳專項科研計劃資助項目(16JK1166);陜西理工大學(xué)研究生創(chuàng)新基金資助項目(SLGYCX1610)

        TH161

        10.3969/j.issn.1004-132X.2016.17.009

        趙永強,男, 1976年生。陜西理工大學(xué)機械工程學(xué)院副教授。主要研究方向為先進成形工藝與精密加工技術(shù)設(shè)備及控制、電動汽車及節(jié)能技術(shù)。發(fā)表論文40余篇。王耀鋒,男,1988年生。陜西理工大學(xué)機械工程學(xué)院碩士研究生。侯紅玲,女,1976年生。陜西理工大學(xué)機械工程學(xué)院副教授。

        猜你喜歡
        定心滾輪凸輪
        定心以應(yīng)世
        凸輪零件的內(nèi)花鍵拉削工藝的自動化生產(chǎn)線
        滾輪滑雪運動多元價值與推廣策略研究
        冰雪運動(2019年2期)2019-09-02 01:41:22
        用三爪自定心卡盤夾持車削偏心工件
        基于UG&VERICUT的弧面凸輪多軸數(shù)控加工仿真實現(xiàn)
        基于MATLAB的盤形凸輪逆向工程
        馬一?。阂x書,先定心
        1+1=150
        凸輪機構(gòu)在“S”型無碳小車中應(yīng)用的可行性
        滾輪移動式鑄鐵機改滾輪固定式鑄鐵機設(shè)計改造
        亚洲人免费| 亚洲自偷精品视频自拍| 东北老女人高潮大喊舒服死了| 蜜桃精品免费久久久久影院| 色噜噜狠狠色综合欧洲| 久久精品视频日本免费| 中文有码无码人妻在线| 精品国产aⅴ无码一区二区| 国产精品99精品一区二区三区∴| 久久精品成人一区二区三区蜜臀| 久久精品熟女亚洲av麻| а√资源新版在线天堂| 精品国产亚洲一区二区在线3d| 日本一区二区三区一级免费| 爆操丝袜美女在线观看| 欧美丰满熟妇xxxx性| 99国产超薄丝袜足j在线观看| 精品久久日产国产一区| 蜜桃视频网站在线观看一区| 亚洲日产精品一二三四区| 国产在线不卡AV观看| 男女啪啪免费视频网址| 无人区乱码一区二区三区| 国产成人啪精品视频免费软件| 中文字幕无码免费久久99| 久久精品视频日本免费| 草草地址线路①屁屁影院成人| 婷婷亚洲综合五月天小说| 国产一区二区在线观看我不卡| 青青河边草免费在线看的视频| 国产一区二区三区免费观看在线 | 国产精品久久综合桃花网| 国产av一区麻豆精品久久| 欧美性受xxxx黑人猛交| 无码人妻少妇色欲av一区二区 | 女同性恋一区二区三区四区| 日韩精品综合一本久道在线视频| 风间由美性色一区二区三区 | 日韩精品免费一区二区三区观看| 六月婷婷久香在线视频| 精品人妻丰满久久久a|