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

        ?

        高速HSK熱裝工具系統(tǒng)結(jié)合部的參數(shù)辨識(shí)

        2015-12-29 02:40:54鄭登升,陳建,陶德飛
        中國(guó)機(jī)械工程 2015年18期
        關(guān)鍵詞:參數(shù)辨識(shí)結(jié)合部可靠性

        高速HSK熱裝工具系統(tǒng)結(jié)合部的參數(shù)辨識(shí)

        鄭登升1陳建1陶德飛1王貴成1,2

        1.江蘇大學(xué),鎮(zhèn)江,2120132.南通理工學(xué)院,南通,226002

        摘要:HSK工具系統(tǒng)的動(dòng)態(tài)特性直接影響機(jī)床的加工效率與質(zhì)量,獲得準(zhǔn)確的工具系統(tǒng)結(jié)合部參數(shù)又是精確預(yù)測(cè)HSK工具系統(tǒng)動(dòng)態(tài)特性的基礎(chǔ)。將熱裝刀柄-刀具結(jié)合部簡(jiǎn)化為彈簧-阻尼模型,基于Timoshenko梁理論建立了HSK熱裝刀柄-刀具結(jié)合部的有限元模型,采用彈塑性理論辨識(shí)結(jié)合部的剛度和阻尼,獲得了HSK熱裝刀柄-刀具結(jié)合部?jī)?nèi)沿軸向不同位置的剛度和阻尼。進(jìn)而建立了HSK熱裝工具系統(tǒng)的有限元?jiǎng)恿W(xué)模型,分析工具系統(tǒng)刀尖點(diǎn)的頻響函數(shù)。最后進(jìn)行了HSK熱裝工具系統(tǒng)頻響函數(shù)的實(shí)驗(yàn)測(cè)量,與理論頻響函數(shù)相比較,驗(yàn)證了將結(jié)合部等效為彈簧-阻尼的辨識(shí)方法的合理性和有效性。

        關(guān)鍵詞:高速加工;可靠性;HSK熱裝工具系統(tǒng);結(jié)合部;參數(shù)辨識(shí)

        中圖分類(lèi)號(hào):TG502;TH113

        收稿日期:2014-12-05

        基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(51075192);國(guó)家重大科技專(zhuān)項(xiàng)(ZX201304009031);南通市重點(diǎn)實(shí)驗(yàn)室建設(shè)資助項(xiàng)目(CP12014002)

        作者簡(jiǎn)介:鄭登升,男,1989年生。江蘇大學(xué)機(jī)械工程學(xué)院碩士研究生。主要研究方向?yàn)楦咚偌庸さ?。陳建,男?987年生。江蘇大學(xué)機(jī)械工程學(xué)院博士研究生。陶德飛,男,1990年生。江蘇大學(xué)機(jī)械工程學(xué)院碩士研究生。王貴成(通信作者),男,1955年生。江蘇大學(xué)機(jī)械工程學(xué)院教授、南通理工學(xué)院機(jī)電工程學(xué)院教授。

        ParameterIdentificationofHigh-speedHSKShrinkToolingSystemCouplingPortion

        ZhengDengsheng1Chen Jian1Tao Defei1Wang Guicheng1,2

        1.JiangsuUniversity,Zhenjiang,Jiangsu,212013

        2.NantongPolytechnicCollege,Nantong,Jiangsu,226002

        Abstract:The dynamic characteristics of HSK tool system affected the machine processing efficiency and the workpiece surface quality,The basis of accurate prediction of the dynamic characteristics of HSK tool system was for obtaining coupling portion precise parameters of tool system.This paper simplified coupling portion of HSK shrink toolholder as spring-damping model,based on Timoshenko beam theory a finite element model of the HSK shrink-tool coupling portion was established,the stiffness and damping of coupling portion was identified by using elastoplastic theory,stiffness and damping of HSK shrink-tool coupling portion was got at different axial locations.And dynamics finite element model of HSK shrink tool system was established,frequency response function(FRF) of system tool tip point was analyzed.And experimental measurements of HSK shrink tool system FRF were made,which was compared with the theoretical FRF.Equivalent to a spring-damping the reasonable and effective identification method of coupling portion was validated.

        Keywords:highspeedmachining;reliability;HSKshrinktoolingsystem;couplingportion;parameteridentification

        0引言

        HSK熱裝工具系統(tǒng)由HSK熱裝刀柄與刀具組成[1],包含HSK熱裝刀柄-刀具結(jié)合面。由工具系統(tǒng)結(jié)合部剛度引起的變形量占刀尖點(diǎn)變形量的26%~50%[2],是影響切削穩(wěn)定性最重要的因素之一。高速加工中刀尖顫振一直受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注,顫振對(duì)加工質(zhì)量和加工效率影響非常大。刀尖顫振主要是由主軸系統(tǒng)結(jié)合部的柔性連接引起的,精確辨識(shí)主軸系統(tǒng)結(jié)合部剛度在預(yù)測(cè)機(jī)床主軸系統(tǒng)穩(wěn)定性的過(guò)程中起著關(guān)鍵的作用。國(guó)內(nèi)外學(xué)者對(duì)主軸系統(tǒng)結(jié)合部作了大量的研究:Erturk等[3]采用有限差分法和實(shí)驗(yàn)測(cè)量相結(jié)合的方式對(duì)彈簧刀柄-刀具結(jié)合部進(jìn)行了參數(shù)識(shí)別;Altintas等[4]建立主軸-刀柄-刀具系統(tǒng)的運(yùn)動(dòng)方程,通過(guò)測(cè)試刀尖點(diǎn)頻響函數(shù),辨識(shí)出了結(jié)合面的動(dòng)力學(xué)參數(shù);汪博等[5]基于并聯(lián)轉(zhuǎn)子的建模思想,將主軸-刀柄-刀具結(jié)合面簡(jiǎn)化為分布彈簧,并完成了主軸系統(tǒng)的有限元建模。這些研究對(duì)深入揭示高速加工工具系統(tǒng)結(jié)合面的參數(shù)變化規(guī)律具有重要作用。

        本文基于有限元法識(shí)別了HSK熱裝工具系統(tǒng)結(jié)合部的4個(gè)剛度參數(shù)和2個(gè)阻尼參數(shù),并基于Timoshenko梁理論建立了HSK熱裝工具系統(tǒng)梁?jiǎn)卧P停瑢⒆R(shí)別出的結(jié)合部參數(shù)代入熱裝工具系統(tǒng)有限元?jiǎng)恿W(xué)模型中,計(jì)算HSK熱裝工具系統(tǒng)的頻響函數(shù),并進(jìn)行實(shí)驗(yàn)驗(yàn)證。

        1參數(shù)辨識(shí)的有限元模型及理論基礎(chǔ)

        1.1熱裝刀柄-刀具結(jié)合部參數(shù)辨識(shí)的有限元模型

        本文主要使用有限元法來(lái)識(shí)別結(jié)合部的剛度和阻尼參數(shù)。HSK熱裝工具系統(tǒng)的三維模型如圖1所示。

        圖1 HSK熱裝工具系統(tǒng)的三維模型

        在進(jìn)行參數(shù)識(shí)別時(shí),只需要將熱裝刀柄-刀具結(jié)合部代入ANSYS軟件中進(jìn)行計(jì)算。有限元模型的熱裝刀柄夾持直徑為16mm,刀具直徑為16mm,總長(zhǎng)度為150mm,刀柄與刀具夾持長(zhǎng)度為40mm。建立模型以后,設(shè)置熱裝刀柄與刀具的材料,HSK熱裝所采用的材料是熱作磨具鋼4Cr5MoSiVA(H13),刀具采用的材料為硬質(zhì)合金,具體參數(shù)如表1所示;再分別對(duì)刀柄和刀具進(jìn)行網(wǎng)格劃分,并在結(jié)合部處設(shè)置接觸對(duì),目標(biāo)面為刀具外圓面,接觸面為夾頭內(nèi)孔面,接觸摩擦因數(shù)為0.2,刀具是通過(guò)熱脹冷縮的原理裝進(jìn)刀柄中,裝配好的刀柄-刀具等效過(guò)盈量為10μm[6];將刀柄末端設(shè)置為固定,刀具末端上施加相應(yīng)的載荷,便可進(jìn)行熱裝刀柄-刀具結(jié)合部的參數(shù)識(shí)別。熱裝刀柄-刀具結(jié)合部有限元模型如圖2所示。

        表1 刀柄和刀具的材料屬性

        圖2 熱裝刀柄-刀具結(jié)合部參數(shù)辨識(shí)的有限元模型

        1.2熱裝刀柄-刀具結(jié)合部參數(shù)辨識(shí)的基礎(chǔ)知識(shí)

        HSK主軸系統(tǒng)軸向剛度與扭轉(zhuǎn)剛度都非常大,在此兩方向上的振動(dòng)影響很小,但徑向剛度有限,HSK主軸系統(tǒng)切削不平穩(wěn)主要受到徑向振動(dòng)影響;HSK主軸系統(tǒng)為軸對(duì)稱(chēng)結(jié)構(gòu),雖然Timoshenko梁有6個(gè)自由度,但本文只需考慮HSK主軸系統(tǒng)XY平面內(nèi)的徑向動(dòng)態(tài)特性。因此熱裝刀柄-刀具結(jié)合部的參數(shù)定義為4個(gè)剛度參數(shù)和2個(gè)阻尼參數(shù),分別為:kyf、kym、kθf(wàn)、kθm、cyf和cym。圖3為HSK熱裝刀柄-刀具結(jié)合部參數(shù)示意圖。根據(jù)彈性理論中的胡克定理,熱裝刀柄-刀具結(jié)合部的參數(shù)辨識(shí)原理如下。由圖3所示,在刀具外圓柱面上施加y方向力f(t),刀柄-刀具結(jié)合面上沿z方向上任意一點(diǎn)A產(chǎn)生的位移為Δy,A點(diǎn)的剛度為kyf,阻尼為cyf。同樣施加力f(t)時(shí),結(jié)合面上沿z方向上任意一點(diǎn)A產(chǎn)生的轉(zhuǎn)動(dòng)角度為Δθ,A點(diǎn)的剛度為kθf(wàn)。撤去力f(t)后,施加力矩m(t),A點(diǎn)產(chǎn)生的位移為Δy,A點(diǎn)的剛度為kym,阻尼為cym。同樣的施加力矩m(t),A點(diǎn)產(chǎn)生的轉(zhuǎn)動(dòng)角度為Δθ,A點(diǎn)的剛度為kθm,阻尼為cθm。

        圖3 熱裝刀柄-刀具結(jié)合部參數(shù)示意圖

        有限元仿真主要為了識(shí)別結(jié)合面的剛度和阻尼參數(shù),因此可以忽略刀柄與刀具的自身變形。有限元模擬主要為兩個(gè)時(shí)間步的計(jì)算,第一個(gè)時(shí)間步為加載刀柄-刀具結(jié)合部的過(guò)盈量,第二個(gè)時(shí)間步為加載依次遞增的力和力矩。kyf和kθf(wàn)為受力載荷下的剛度,kym和kθm為受力矩載荷下的剛度。在熱裝刀柄-刀具結(jié)合部模型上施加力/力矩,且力/力矩是加載在靠近夾頭/刀具配合處的3個(gè)節(jié)點(diǎn)上,如圖4所示。這樣加載是為了減小刀具的局部變形,提高仿真精度。求解以后,取出結(jié)合面上沿z軸方向變化的位移參數(shù),并繪制每一位置的位移-力曲線,斜率即為每一位置的剛度值。

        圖4 熱裝刀柄-刀具有限元模型的力加載示意圖

        熱裝刀柄-刀具結(jié)合部的阻尼參數(shù)為cyf和cym。cyf和cym是刀柄-刀具結(jié)合部在力/力矩的作用下沿z軸方向產(chǎn)生微小的移動(dòng)而形成的位移阻尼。由于轉(zhuǎn)動(dòng)角度變化很小,可以將每一接觸單元沿y軸滑動(dòng)位移的絕對(duì)值來(lái)代替轉(zhuǎn)動(dòng)角度。并將此阻尼等效為黏性阻尼,其計(jì)算公式為

        (1)

        Fd,n=μpA

        式中,Fd,n為每個(gè)單元的摩擦阻力;|ln|為結(jié)合面每一個(gè)接觸單元沿z軸或y軸方向上滑動(dòng)位移的絕對(duì)值;p、A分別為接觸單元的壓力和面積;ω為接觸單元的角速度。

        與上述剛度識(shí)別方法一樣,分別在熱裝刀柄-刀具結(jié)合部模型上施加力和力矩求解以后,選取每一接觸單元的壓力、面積和z軸或y軸方向上滑動(dòng)位移的絕對(duì)值,進(jìn)而代入到式(1)中可求出相應(yīng)的阻尼。由于靠近中軸處單元的滑動(dòng)幾乎為0,故使用式(1)計(jì)算等效阻尼時(shí)會(huì)得到很大的阻尼值,這是因?yàn)橹休S附近的微小滑動(dòng)由幾何約束造成,這些微小的滑動(dòng)并不能夠正確的表示阻尼[6]。所以,本研究將中軸±15°內(nèi)的接觸單元不納入等效阻尼計(jì)算。

        2HSK熱裝刀柄-刀具結(jié)合部的參數(shù)識(shí)別

        2.1剛度參數(shù)的識(shí)別

        分別在熱裝刀柄-刀具結(jié)合部模型上施加力/力矩,力的加載大小分別為250N、500N、750N、1000N和1250N,力矩的大小分別為2.5N·m、5.0N·m、7.5N·m、10.0N·m和12.5N·m。在計(jì)算kθf(wàn)和kθm時(shí),計(jì)算z軸方向相鄰兩點(diǎn)的位移,可得到刀柄-刀具結(jié)合部沿z軸方向各個(gè)位置的轉(zhuǎn)動(dòng)角度Δθ,繪制每一位置的角度-力/力矩曲線,其斜率即為剛度kθf(wàn)和kθm。圖5所示為刀柄-刀具結(jié)合部沿z軸方向變化的剛度。

        2.2阻尼參數(shù)的識(shí)別

        分別在熱裝刀柄-刀具結(jié)合部模型上施加力和力矩,力為250N、500N、750N、1000N和1250N,力矩為2.5N·m、5.0N·m、7.5N·m、10.0N·m、12.5N·m。圖6a為作用力250N和1250N下的阻尼值cyfω,圖6b為力矩2.5N·m和12.5N·m下的阻尼值cymω。

        3辨識(shí)結(jié)果與實(shí)驗(yàn)驗(yàn)證

        3.1刀尖的頻響函數(shù)的有限元計(jì)算

        (a)沿z軸方向變化的剛度k yf

        (b)沿z軸方向變化的剛度k θf(wàn)

        (c)沿z軸方向變化的剛度k ym

        (d)沿z軸方向變化的剛度k θm 圖5 熱裝刀柄-刀具結(jié)合部z軸方向的剛度

        根據(jù)實(shí)際HSK熱裝工具系統(tǒng)結(jié)構(gòu)建立考慮結(jié)合面的HSK熱裝刀柄-刀具有限元模型,HSKA63刀柄和刀具的具體尺寸如表2、表3所示。其中1~10段為HSK刀柄錐部,11~19段為刀柄中部,20~36段為熱裝夾頭部分,37~57段為刀具部分,結(jié)合部的配合部分為27~36刀柄段和37~46刀具段。相關(guān)刀柄和刀具的材料如表1所示。對(duì)于結(jié)合部的連接采用彈簧-阻尼連接,在ANSSYS中使用MATRIX27剛度矩陣模擬彈簧連接和MATRIX28阻尼矩陣模擬阻尼連接。參考文獻(xiàn)[3],結(jié)合部采用5段分布彈簧最為合適,沿軸向不同位置的具體剛度阻尼參數(shù)如表4所示。建立好的梁?jiǎn)卧狧SK熱裝工具系統(tǒng)的動(dòng)力學(xué)模型如圖7所示。

        3.2實(shí)驗(yàn)設(shè)備及檢測(cè)方法

        實(shí)驗(yàn)所使用的HSK熱裝工具系統(tǒng)與上述參數(shù)識(shí)別所使用的模型相同,刀柄采用HSKA63型刀柄,刀具為整體硬質(zhì)合金四刃立銑刀。

        高速切削動(dòng)力學(xué)仿真與切削參數(shù)優(yōu)化系統(tǒng)(CUTPRO)是由加拿大UBC大學(xué)制造自動(dòng)化實(shí)驗(yàn)室(MAL)開(kāi)發(fā)的,該測(cè)試系統(tǒng)可以用來(lái)仿真加工和機(jī)床檢測(cè)。本實(shí)驗(yàn)使用CUTPRO系統(tǒng)的切削分析模塊,通過(guò)沖擊實(shí)驗(yàn),該模塊可以很好地測(cè)量刀尖點(diǎn)的FRF。實(shí)驗(yàn)中采用力錘對(duì)HSK熱裝工具系統(tǒng)的刀尖點(diǎn)進(jìn)行激勵(lì),選用金屬錘頭,力錘型號(hào)為9722A500。采集力信號(hào)的傳感器和放大器都集成在力錘中。采集響應(yīng)信號(hào)的傳感器為8778A500型加速度傳感器。搭建的HSK熱裝工具系統(tǒng)的測(cè)試系統(tǒng)如圖8所示。

        (a)沿z軸變化的阻尼c yfω

        (b)沿z軸變化的阻尼c ymω 圖6 熱裝刀柄-刀具結(jié)合部z軸方向的阻尼

        mm

        表3 刀具的尺寸  mm

        表4 結(jié)合部的剛度和阻尼參數(shù)

        圖7 HSK熱裝工具系統(tǒng)的有限元?jiǎng)恿W(xué)模型

        圖8 HSK熱裝工具系統(tǒng)的測(cè)試系統(tǒng)

        3.3理論分析結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比

        利用Timoshenko梁?jiǎn)卧ATRIX27單元和MATRIX28單元建立了HSK熱裝刀柄-刀具結(jié)合面的模型,并將識(shí)別的結(jié)合面的每個(gè)位置的參數(shù)代入模型中,隨后計(jì)算熱裝工具系統(tǒng)的響應(yīng),最后得到刀尖點(diǎn)的頻響函數(shù)。分別計(jì)算了實(shí)驗(yàn)、結(jié)合部剛性連接和彈性連接的頻響函數(shù),結(jié)果如圖9所示。

        圖9 理論結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比

        從圖9中還可以看出,以剛性連接處理結(jié)合面的仿真結(jié)果偏高于實(shí)驗(yàn)結(jié)果。前3階彈性連接模型的準(zhǔn)確性明顯優(yōu)越于剛性連接模型的準(zhǔn)確性,特別是第2階剛性連接模型的結(jié)果與實(shí)驗(yàn)結(jié)果偏差很大。第4階彈性連接模型的誤差有點(diǎn)偏大,但總體上來(lái)看,彈性連接模型與實(shí)驗(yàn)的吻合度還是很高的。將結(jié)合面等效為彈簧-阻尼模型的Timoshenko梁?jiǎn)卧P偷姆抡娼Y(jié)果更接近實(shí)驗(yàn)結(jié)果。驗(yàn)證了本研究提出的將HSK熱裝工具系統(tǒng)接合面等效為彈簧-阻尼參數(shù)辨識(shí)方法的可行性和有效性。同時(shí)相比實(shí)驗(yàn)測(cè)試的方法,Timoshenko梁?jiǎn)卧P筒粌H可以準(zhǔn)確地預(yù)測(cè)刀尖點(diǎn)的頻響函數(shù),還可以省去大量的測(cè)量,節(jié)省時(shí)間。Timoshenko梁?jiǎn)卧P驮诜抡孢\(yùn)算的速度上也明顯快于剛性連接的模型。

        4結(jié)論

        (1)將高速熱裝刀柄-刀具結(jié)合部簡(jiǎn)化為彈簧-阻尼模型,并基于有限元法辨識(shí)出結(jié)合部?jī)?nèi)各位置的剛度和阻尼參數(shù),為深入研究高速加工工具系統(tǒng)動(dòng)態(tài)特性奠定了基礎(chǔ)。

        (2)基于Timoshenko梁理論建立了HSK熱裝工具系統(tǒng)的有限元?jiǎng)恿W(xué)模型,計(jì)算系統(tǒng)刀尖點(diǎn)頻響函數(shù)。驗(yàn)證了將HSK熱裝工具系統(tǒng)接合面等效為彈簧-阻尼參數(shù)辨識(shí)方法的可行性和有效性,為進(jìn)一步工具系統(tǒng)結(jié)構(gòu)設(shè)計(jì)及工具系統(tǒng)優(yōu)化提供了理論依據(jù)。

        參考文獻(xiàn):

        [1]王貴成,王樹(shù)林,董廣強(qiáng).高速加工工具系統(tǒng)[M].北京:國(guó)防工業(yè)出版社,2005.

        [2]AgapiouJ,RivinE,XieC.Toolholder/SpindleInterfacesforCNCMachineTools[J].AnnalsofIRP,1995,41(1):383-387.

        [3]ErturkA,OzguvenHN,BudakE.EffectAnalysisofBearingandInterfaceDynamicsonToolPointFRFforChatterStabilityinMachineToolsbyUsingaNewAnalyticalModelforSpindle-toolAssemblies[J].InternationalJournalofMachineToolsandManufacture,2007,47:23-32.

        [4]AltintasY,CaoYuzhong.VirtualDesignandOptimizationofMachineToolSpindles[J].CIRPAnnals,ManufacturingTechnology,2005,54:379-382.

        [5]汪博,孫偉,聞邦椿.考慮主軸-刀柄-刀具結(jié)合面的高速主軸系統(tǒng)動(dòng)力學(xué)特性有限元建模[J].機(jī)械工程學(xué)報(bào),2012,48(15),83-89.

        WangBo,SunWei,WenBangchun.TheFiniteElementModelingofHigh-speedSpindleSystemDynamicswithSpindle-holder-toolJoints[J].JournalofMechanicalEngineering,2012,48(15):83-89.

        [6]SchmitzT,PowellK.ShrinkFitToolHolderConnectionStiffness/DampingModelingforFrequencyResponsePredictioninMilling[J].InternationalJournalofMachineTools&Manufacture,2007,47:1368-1380.

        (編輯袁興玲)

        猜你喜歡
        參數(shù)辨識(shí)結(jié)合部可靠性
        可靠性管理體系創(chuàng)建與實(shí)踐
        工電道岔結(jié)合部聯(lián)合管理機(jī)制的探討
        CTCS-3級(jí)列控系統(tǒng)RBC與ATP結(jié)合部異常信息處理
        電子制作(2017年2期)2017-05-17 03:55:06
        基于偏差補(bǔ)償遞推最小二乘法的熒光補(bǔ)償方法
        基于壓電陶瓷的納米定位與掃描平臺(tái)模型辨識(shí)算法研究
        電纜導(dǎo)體溫度間接測(cè)量的方法研究
        近端胃切除和全胃切除術(shù)治療食管胃結(jié)合部癌的臨床對(duì)照研究
        基于相關(guān)度的忙時(shí)話務(wù)量加權(quán)一階局域預(yù)測(cè)模型
        基于可靠性跟蹤的薄弱環(huán)節(jié)辨識(shí)方法在省級(jí)電網(wǎng)可靠性改善中的應(yīng)用研究
        亚洲精品成人无百码中文毛片| 亚洲色大成在线观看| 国产成社区在线视频观看| 亚洲成在人网站天堂日本| 黄片大全视频在线播放| 亚洲综合欧美在线一区在线播放| 日本中文字幕在线播放第1页| 久草精品手机视频在线观看| 海外华人在线免费观看| 亚洲色成人网站www永久| 比比资源先锋影音网| 特级毛片全部免费播放a一级| 日韩精品国产精品亚洲毛片| 免费无码又爽又刺激网站直播| 久久人人爽人人爽人人片亞洲| 99re国产电影精品| 日韩不卡一区二区三区色图| 国产欧美精品一区二区三区四区| a级毛片毛片免费观看久潮喷| 国产思思久99久精品| 日本一区二区三区四区在线视频 | 中文字幕人妻第一区| 狠狠人妻久久久久久综合| 国产在线观看网址不卡一区| 亚洲中文字幕剧情类别| 粉嫩被粗大进进出出视频| 人妻无码视频| 亚洲精品乱码久久麻豆| 蜜桃视频一区二区在线观看| 久久99精品国产99久久6男男| 国产呦系列呦交| 日韩一区二区av极品| 亚洲一区 日韩精品 中文字幕 | 中文不卡视频| 一区二区三区在线乱码| 免费无遮挡无码永久视频| 国产人碰人摸人爱视频| 亚洲成人激情在线影院| 在线免费观看一区二区| 国产成人av免费观看| 亚洲国产成人aⅴ毛片大全|