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

        ?

        基于切深模型的汽車曲軸軸頸巴厘線磨削精度控制

        2019-01-24 05:18:24上海大學(xué)高華鈺沈南燕方明倫
        關(guān)鍵詞:巴厘廓形修整

        上海大學(xué) 李 靜 高華鈺 沈南燕 方明倫

        上海機(jī)床廠有限公司 黃海濤

        汽車曲軸作為發(fā)動機(jī)關(guān)鍵零件,其主軸頸、連桿頸巴厘線的磨削精度對發(fā)動機(jī)性能有較大影響。切入磨削過程中,成形砂輪輪廓將直接反映在曲軸軸頸上,因此研究了基于砂輪修整切深模型的成形砂輪廓形修整誤差在線測量及補(bǔ)償方法,通過控制砂輪廓形修整精度來保證曲軸軸頸巴厘線的磨削精度。采用聲發(fā)射傳感器搭建了砂輪修整過程監(jiān)測系統(tǒng),建立了砂輪圓弧修整過程中任意位置的聲發(fā)射信號均方根值與修整切深的數(shù)學(xué)模型。通過構(gòu)建切深模型參數(shù)與修整進(jìn)給速度和修整圓弧中凸量的函數(shù)關(guān)系,得到了砂輪圓弧修整變參數(shù)切深模型,提高了切深模型的準(zhǔn)確性。在此模型的基礎(chǔ)上,提出了砂輪廓形修整誤差補(bǔ)償策略及實現(xiàn)流程。實驗證明:采用變參數(shù)切深模型可以準(zhǔn)確獲得砂輪實際修整廓形,定量評定砂輪廓形修整誤差,并且修整誤差補(bǔ)償策略能夠有效地保證巴厘線廓形磨削精度。

        一、 前言

        曲軸是汽車發(fā)動機(jī)上的關(guān)鍵零件之一,需求量大,組成要素多,輪廓形狀復(fù)雜,精度要求高。曲軸軸頸圓柱表面母線為具有一定公差要求的圓弧形,主軸頸呈凸形,連桿頸呈凹形,稱為巴厘線(源自德文Ballig)。主軸頸巴厘線能避免高扭矩狀態(tài)下的曲軸軸頸與軸瓦出現(xiàn)邊角接觸;而連桿頸巴厘線則對動壓油膜的建立有利,從而提升發(fā)動機(jī)可靠性[1]。因此曲軸主軸頸、連桿軸頸巴厘線的磨削精度對發(fā)動機(jī)性能有較大影響。

        巴厘線通常是通過對砂輪進(jìn)行成形修整,再由切入磨來實現(xiàn)的。國內(nèi)的李小新等[2]提出采用矢量切點(diǎn)式跟蹤磨削的方法,通過改造砂輪基體結(jié)構(gòu)以及磨削時調(diào)整砂輪線速度實現(xiàn)曲軸巴厘線的磨削。該方法不能夠很好地控制巴厘線形狀及尺寸,且不易操作,因此實用性不高。目前廣泛應(yīng)用的曲軸隨動磨削技術(shù),通過控制C軸和X軸的聯(lián)動實現(xiàn)一次裝夾對所有主軸頸和連桿頸的磨削[3-4],并利用成形砂輪磨削出曲軸軸頸巴厘線。因此,高效率大批量生產(chǎn)中汽車曲軸軸頸巴厘線的中凸或中凹量主要靠砂輪廓形的修整精度保證。

        砂輪廓形在機(jī)修整通過Z-X軸聯(lián)動的圓弧插補(bǔ)實現(xiàn)[5],但由于曲軸磨床加工需求和配置特點(diǎn),通常采用伺服電機(jī)驅(qū)動配合滾珠絲杠傳動實現(xiàn)沿工件軸向的Z軸移動,其運(yùn)動定位精度相對采用直線電機(jī)直驅(qū)的砂輪架進(jìn)給方向的X軸明顯較低,導(dǎo)致磨床Z-X聯(lián)動插補(bǔ)精度較低,由此降低了砂輪輪廓的修整精度,產(chǎn)生曲軸軸徑巴厘線加工誤差。

        因此,在生產(chǎn)節(jié)拍快、連續(xù)工作周期長的汽車發(fā)動機(jī)生產(chǎn)線用曲軸磨床上,為保證曲軸軸頸巴厘線的加工質(zhì)量及穩(wěn)定性,需要監(jiān)測砂輪修形過程,以便及時補(bǔ)償砂輪輪廓的修整誤差,保證修整精度和修整效率。聲發(fā)射(Acoustic Emission, AE)信號因其高頻的特點(diǎn),可避開低頻噪聲的干擾且靈敏度極高,已被廣泛應(yīng)用于對工件磨削及砂輪修整過程的監(jiān)控[6-11]。但此方面研究大多集中在利用聲發(fā)射信號進(jìn)行定性分析,如砂輪磨損狀態(tài)判別[8-9]以及砂輪修整時與修整器的接觸狀態(tài)判別[10-11],而沒有定量分析聲發(fā)射信號與砂輪修整廓型的對應(yīng)關(guān)系。

        本文通過構(gòu)建砂輪修整監(jiān)測平臺,定量分析不同修整參數(shù)下聲發(fā)射信號的特征值與修整切深之間的關(guān)系,從而通過修整過程中的聲發(fā)射信號廓形得到砂輪實際修整廓形的估計值,進(jìn)而獲得廓形中凸量修整誤差。以此為基礎(chǔ),研究修整廓形中凸量補(bǔ)償策略,達(dá)到通過控制砂輪廓形修整精度保證曲軸軸頸巴厘線磨削精度的目的。

        二、巴厘線磨削精度控制系統(tǒng)方案

        曲軸軸頸巴厘線的磨削精度主要靠砂輪廓形精度保證,切入磨削過程中砂輪廓形將直接反映在曲軸軸頸上。因此可以通過控制砂輪廓形的修整精度,來滿足曲軸軸頸巴厘線磨削精度。為此,本文設(shè)計了如圖1所示的砂輪廓形修整精度控制系統(tǒng)方案:通過布置在修整器附近的聲發(fā)射傳感器及相關(guān)采集設(shè)備實時獲取砂輪修整過程中的聲發(fā)射信號;在濾波等信號預(yù)處理的基礎(chǔ)上,基于某位置附近微小時段聲發(fā)射信號均方根值與該位置實際修整切深之間的數(shù)學(xué)模型,對砂輪實際修整廓形進(jìn)行在線估計;根據(jù)砂輪實際修整廓形與理論修整廓形,計算補(bǔ)償后的砂輪修整圓弧半徑,并更新具有中凸量誤差補(bǔ)償能力的砂輪修整NC程序;最后通過CNC系統(tǒng)執(zhí)行此修整程序,以此保證砂輪廓形修整精度,從而實現(xiàn)曲軸軸頸巴厘線磨削精度在線控制。

        圖1 砂輪廓形修整精度的控制系統(tǒng)

        三、基于聲發(fā)射信號的砂輪圓弧廓形修整切深建模

        軸頸巴厘線廓形通常為圓弧,該圓弧的弦高即為巴厘線的中凸或中凹量。砂輪的理論修整廓形與工件巴厘線廓形相同。通過在機(jī)直接測量獲得砂輪廓形的修整質(zhì)量十分困難[12],因此采用聲發(fā)射(AE)信號,在砂輪修平的基礎(chǔ)上,通過研究AE信號與修整切深間的數(shù)學(xué)模型,用某位置附近微小時段的AE信號均方根值來表征該位置的修整切深,進(jìn)而獲得砂輪實際修整圓弧廓形,實現(xiàn)其在機(jī)間接測量。

        1. 砂輪平修切深模型

        磨削過程AE信號的均方根值與法向或切向磨削力的變化曲線相近[13-14],在外圓切入磨削過程中AE信號均方根值與磨削力成一定的比例關(guān)系[15-16]。

        而修整過程中修整器與砂輪之間的相互作用與磨削過程中砂輪與工件的關(guān)系相似,因此修整過程中AE信號的均方根值與修整力之間也有類似的關(guān)系。借鑒磨削力的經(jīng)驗公式[17],砂輪修整過程中AE信號的均方根值與砂輪修整參數(shù)具有類似的函數(shù)關(guān)系,其表達(dá)形式如下:

        式中,

        Xae——AE信號均方根值/V;

        Xgr——修整切深/mm;

        Fz——Z方向進(jìn)給速度/mm·min-1;

        Vs——砂輪轉(zhuǎn)速/m·s-1;

        a0 、a1 、a2 、a3——待定參數(shù)。

        將式(1)變形,則修整切深可表示為

        經(jīng)簡化得砂輪平修時的切深模型:

        在圓弧修整過程中,圓弧中凸量相較于圓弧半徑r非常小,基于微元法的思想,將每一微小時段看作是沿Z軸方向的平修。如圖2所示, 為修整圓弧軌跡,C為圓弧最高點(diǎn),O為圓弧圓心,D為CO與AB的交點(diǎn)。已知修整圓弧中凸量=ap,修整圓弧寬度在直角△OAD中,滿足如下關(guān)系:

        得圓弧半徑:

        若以A為圓弧起點(diǎn),某時刻修整到P點(diǎn),修整器沿Z方向運(yùn)動了zi距離,此時對應(yīng)的圓心角θ滿足如下關(guān)系:

        根據(jù)圖中的幾何關(guān)系,進(jìn)給速度F與Z軸也呈θ角,其Z軸分量:

        圖2 圓弧修整原理幾何示意

        由式(5)~(7)可得圓弧軌跡中任意位置zi處的Z方向進(jìn)給速度

        將式(8)代入式(3),即可得到該位置砂輪圓弧修整的切深模型。

        四、砂輪圓弧廓形修整切深模型的改進(jìn)

        1 .砂輪圓弧廓形修整切深模型參數(shù)辨識

        砂輪圓弧廓形修整過程中,AE信號均方根值及Z方向進(jìn)給速度均在變化,而砂輪轉(zhuǎn)速始終為恒定值,因此切深模型中與砂輪轉(zhuǎn)速有關(guān)的一項也為常系數(shù)項,可將該項與常系數(shù)項 1合并處理,由此式(3)的切深模型將簡化為:

        式中,b0為合并后的常系數(shù)項指數(shù),為原b0與b3lgVs之和。

        令上式兩邊分別取對數(shù),將該非線性模型轉(zhuǎn)換為線性模型,得到:

        則各參數(shù)的求解可轉(zhuǎn)換為多元線性回歸方程系數(shù)的求解。根據(jù)修整進(jìn)給速度F、修整圓弧中凸量ap、修整圓弧寬度d,利用式(8)計算出各位置Z方向進(jìn)給速度Fz。同時聯(lián)立修整過程中各位置附近微小時段內(nèi)的AE信號均方根值Xae及實際修整廓形測量值Xgr作為樣本輸入模型,基于最小二乘的原則,將求得的式(10)中系數(shù)最優(yōu)值作為砂輪圓弧廓形修整切深模型參數(shù)辨識結(jié)果。

        但由于圓弧修整過程的復(fù)雜性,如不同修整參數(shù)產(chǎn)生的系統(tǒng)伺服誤差不同等,使得不同修整參數(shù)求解得到的模型參數(shù)不盡相同。若直接將模型參數(shù)定為常數(shù),則模型準(zhǔn)確性不佳。因此本文將通過分析并建立模型參數(shù)與修整參數(shù)的關(guān)系,來改進(jìn)砂輪圓弧廓形修整切深模型。

        2.修整參數(shù)對模型參數(shù)的影響

        為了探究不同修整參數(shù)對模型參數(shù)的影響,并對切深模型進(jìn)行改進(jìn),本文設(shè)計了單因素試驗進(jìn)行分析。對標(biāo)實際生產(chǎn)中的主軸頸巴厘線磨削,即欲將工件磨成中凸廓形,砂輪相應(yīng)地需修整成中凹廓形。為了保證修整切深與廓形的準(zhǔn)確映射,每次砂輪圓弧廓形修整之前應(yīng)保證砂輪已修平。

        實際應(yīng)用中,砂輪轉(zhuǎn)速通常根據(jù)砂輪類型和磨床能力選定后無需調(diào)整。因此本文不考慮砂輪轉(zhuǎn)速變化的影響,將模型參數(shù)b0、b1、b2均看作是進(jìn)給速度F及修整圓弧中凸量ap的函數(shù)。為此,本文分別設(shè)計了中凸量為6μm、進(jìn)給速度在40mm/min~150mm/min范圍變化的單因素試驗及進(jìn)給速度為100mm/min、中凸量在1~12μm范圍變化的單因素試驗。試驗過程中砂輪轉(zhuǎn)速恒定為30m/s?;谏闲」?jié)中的參數(shù)辨識方法,根據(jù)試驗結(jié)果求解得到的參數(shù)及其變化趨勢分別如圖3、圖4所示。其中圖3為僅進(jìn)給速度變化時模型參數(shù)變化趨勢,圖4為僅圓弧中凸量變化時模型參數(shù)變化趨勢。

        圖3 進(jìn)給速度變化時參數(shù)變化趨勢

        圖4 圓弧中凸量變化時參數(shù)變化趨勢

        圖3結(jié)果表明:進(jìn)給速度在40mm/min至150mm/min之間變化時,參數(shù)b1與b0變化趨勢相同而與b2變化趨勢相反,三個參數(shù)的變化規(guī)律基本相同:進(jìn)給速度在50mm/min至80mm/min之間變化時,參數(shù)變化劇烈,在60mm/min至70mm/min間出現(xiàn)最值;其余位置變化平緩,基本呈高斯函數(shù)變化趨勢。圖4結(jié)果表明:圓弧中凸量在1~12μm之間變化時,參數(shù)b1與b2變化趨勢相同而與b0變化趨勢相反,三個參數(shù)的變化規(guī)律基本相同:中凸量在2μm以下時,參數(shù)有明顯下降(上升),在2μm以上時則變化平緩,基本呈冪函數(shù)變化趨勢。

        3. 砂輪圓弧修整變參數(shù)切深模型

        基于上節(jié)分析,分別以高斯函數(shù)及冪函數(shù)為基礎(chǔ),以二次多項式補(bǔ)充修正,對兩組單因素試驗的參數(shù)變化曲線進(jìn)行擬合,進(jìn)而建立砂輪圓弧修整變參數(shù)切深模型。

        若記:

        則可將參數(shù)曲線擬合函數(shù)表示為以下形式:

        為提高切深模型的準(zhǔn)確性,設(shè)計了如圖5所示的砂輪圓弧修整變參數(shù)切深模型參數(shù)辨識流程:對模型估計相對誤差進(jìn)行補(bǔ)償,計算產(chǎn)生新的參數(shù),再進(jìn)行參數(shù)曲線擬合,直至模型估計相對誤差小于誤差限10%,可實現(xiàn)基于模型估計誤差的參數(shù)識別結(jié)果自調(diào)整。其中新的參數(shù)基于以下規(guī)則產(chǎn)生:

        由于模型中包含參數(shù)b0的一項為常系數(shù)項,為方便計算,將模型估計誤差全補(bǔ)償至b0,則產(chǎn)生的新參數(shù)值為

        式中j表示參數(shù)調(diào)整次數(shù),j=0、1、2、……。其中b0(0)、b1(0)、b2(0)為上節(jié)中得到的參數(shù)。

        圖5 砂輪圓弧修整變參數(shù)切深模型參數(shù)辨識流程

        最終求得擬合后結(jié)果為

        依據(jù)參數(shù)擬合函數(shù)的基本函數(shù)構(gòu)成,可將上述單因素試驗結(jié)果進(jìn)行整合,假定切深模型各參數(shù)具有如下數(shù)學(xué)表達(dá)形式

        式中Bi為 的系數(shù)矩陣,i=0、1、2…

        根據(jù)式(15)~(21)即可求得參數(shù)b0、b1、b2的系數(shù)矩陣B0、B1、B2,從而可得到b0、b1、b2關(guān)于進(jìn)給速度F及修整圓弧中凸量ap的函數(shù)表達(dá)式,進(jìn)而代入式(9)可得到砂輪圓弧修整任意位置變參數(shù)切深模型:

        五、砂輪修整廓形精度控制策略

        1.砂輪圓弧修整中凸量誤差補(bǔ)償規(guī)則

        根據(jù)修整過程中砂輪寬度上各位置得到的AE信號均方根值,由建立的砂輪圓弧修整變參數(shù)切深模型,可以獲得連續(xù)的砂輪實際修整廓形。與砂輪理論修整廓形對比,可以發(fā)現(xiàn)實際修整廓形仍近似為圓弧,砂輪修整廓形誤差整體反映在圓弧的中凸量上。由于存在伺服滯后誤差以及修整筆磨損等原因,實際修整廓形的中凸量總是小于理論修整廓形,中凸量絕對誤差(即模型估計得到的實際中凸量與理論中凸量的絕對誤差) 。因此保持修整圓弧起點(diǎn)和終點(diǎn)不變的情況下,本文通過增大修整程序中的圓弧軌跡半徑來實現(xiàn)砂輪廓形中凸量的補(bǔ)償。

        如圖6所示,虛線表示補(bǔ)償前后的理論修整廓形,而實線表示補(bǔ)償前后的實際修整廓形;為補(bǔ)償前的理論修整廓形,其圓弧圓心位于O點(diǎn),圓弧中凸量=a;為補(bǔ)償后的理論修整廓形,其圓弧圓心位于O′點(diǎn),圓弧中凸量=ap′。根據(jù)幾何關(guān)系得到的補(bǔ)償前的理論圓弧半徑為:

        將補(bǔ)償前后理論圓弧廓形中凸量的增量設(shè)定為與補(bǔ)償前的中凸量絕對誤差值ε相等如圖所示,補(bǔ)償后的理論圓弧中凸量變?yōu)?,得因此補(bǔ)償后的理論圓弧半徑:

        圖6 中凸量誤差補(bǔ)償原理幾何示意

        基于上節(jié)單因素試驗,進(jìn)行中凸量誤差補(bǔ)償試驗,補(bǔ)償前后試驗結(jié)果如表1所示:

        表1 中凸量誤差補(bǔ)償試驗方案及結(jié)果

        試驗結(jié)果表明:補(bǔ)償前,實際中凸量與理論中凸量相對誤差較大,多數(shù)都在20%以上;而補(bǔ)償后,實際中凸量與理論中凸量的相對誤差均控制在5%以內(nèi),補(bǔ)償效果良好。

        2.砂輪廓形修整誤差補(bǔ)償流程

        根據(jù)汽車曲軸巴厘線磨削精度控制方案,本文設(shè)計了如圖7所示的砂輪圓弧修整中凸量誤差補(bǔ)償流程。其中中凸量相對誤差δ的范圍設(shè)定由以下規(guī)則確定:

        已知實測廓形的相對誤差δms(即實測廓形中凸量ap測量量與理論中凸量ap的相對誤差)滿足:

        模型估計相對誤差δpr(e即估計廓形中凸量與實際測量得廓形中凸量的相對誤差)滿足:

        中凸量相對誤差δ(即由模型估計得到的實際修整廓形中凸量與理論中凸量ap的相對誤差)滿足

        結(jié)合式(27)~(28)得

        圖7 汽車曲軸巴厘線磨削精度控制流程

        六、巴厘線磨削精度控制實驗驗證

        為了完整驗證軸頸巴厘線磨削精度控制方案的有效性,在偏心軸非圓磨床H367上進(jìn)行了驗證實驗。如圖8所示,實驗以圓棒代替實際生產(chǎn)中的主軸頸,聲發(fā)射傳感器利用磁座吸附在修整點(diǎn)附近進(jìn)行修整過程中AE信號的獲??;根據(jù)砂輪修整廓形中凸量精度控制流程進(jìn)行砂輪圓弧廓形修整后,對工件進(jìn)行外圓切入磨削,將砂輪廓形印刻在工件上;磨削完成后,工件轉(zhuǎn)移至SE13-J10曲軸測量機(jī)上進(jìn)行實際修整廓形的測量。實驗隨機(jī)選取了七組修整參數(shù),實驗條件見表2。

        圖8 巴厘線磨削精度控制驗證實驗

        表2 中凸量誤差補(bǔ)償試驗方案及結(jié)果

        在砂輪修整過程中,傳感器采集的信號并不只包含由于砂輪修整而產(chǎn)生的有效信號,而是各種噪聲與試驗有效信號的交疊,對試驗結(jié)果會造成影響。鑒于噪聲和有效信號往往具有不同的頻率分布特性,通過快速傅里葉(FFT)頻譜分析[18],對比修整前與修整時聲發(fā)射信號頻譜的變化,確定砂輪修整時的有效信號的頻段分布范圍集中在100kHz到225kHz之間。

        聲發(fā)射的信號強(qiáng)度一般選用其均方根值來表征,見式(31)。式中V(t)為任意時刻采集的聲發(fā)射信號值,T為時間常數(shù)。本文取時間常數(shù)為100ms。

        驗證實驗結(jié)果如表3所示,試驗結(jié)果表明,最終廓形誤差均控制在10%以內(nèi),驗證了砂輪圓弧修整切深模型的正確性,以及軸頸巴厘線磨削精度控制方案的可行性。

        表3 巴厘線磨削精度控制驗證實驗

        七、結(jié)論

        為實現(xiàn)曲軸軸徑巴厘線磨削精度控制,研究了成形砂輪廓形修整誤差在線監(jiān)測及補(bǔ)償方法,具體包括:

        (1)通過類比磨削力模型,建立了基于聲發(fā)射信號的砂輪圓弧修整切深模型,從而獲得了修整過程中聲發(fā)射信號廓形與砂輪實際修整廓形之間的定量關(guān)系,利用聲發(fā)射信號估計出砂輪實際修整廓形。

        (2)為提高切深模型的準(zhǔn)確性,基于單因素試驗結(jié)果,分析模型參數(shù)隨修整進(jìn)給速度和修整圓弧中凸量的變化規(guī)律,發(fā)現(xiàn)模型參數(shù)隨前者的變化基本呈高斯函數(shù)形式,隨后者變化基本呈冪函數(shù)形式。以此為基礎(chǔ),構(gòu)建了砂輪圓弧修整變參數(shù)切深模型,并設(shè)計了模型參數(shù)辨識流程及參數(shù)產(chǎn)生規(guī)則,實現(xiàn)了基于模型誤差的參數(shù)自調(diào)整。

        (3)通過分析實際修整廓形,確定廓形誤差主要體現(xiàn)在圓弧中凸量上,由此提出通過改變修整圓弧半徑進(jìn)行補(bǔ)償?shù)囊?guī)則。并設(shè)計了砂輪圓弧修整中凸量誤差自動補(bǔ)償流程。

        上述砂輪圓弧修整變參數(shù)切實模型及中凸量誤差補(bǔ)償策略,通過實際修整及磨削實驗,驗證了其在巴厘線磨削精度控制中的有效性。在具體應(yīng)用中,由于聲發(fā)射傳感器安裝位置對采集信號有一定影響,因此下一步將對不同安裝位置采得信號對切深模型參數(shù)的影響進(jìn)行研究,以增強(qiáng)模型的魯棒性。

        猜你喜歡
        巴厘廓形修整
        懂經(jīng)濟(jì)學(xué)的巴厘獼猴
        懂經(jīng)濟(jì)學(xué)的巴厘獼猴
        意林(2021年17期)2021-09-26 02:05:50
        懂經(jīng)濟(jì)學(xué)的巴厘獼猴
        知識窗(2021年7期)2021-08-28 18:09:33
        意在輕閑
        力量廓形
        單點(diǎn)金剛石修整筆的現(xiàn)狀分析及發(fā)展趨勢
        窄氣感廓形
        ‘巴厘’菠蘿離體培養(yǎng)技術(shù)
        大衣+牛仔褲 走路得帶風(fēng)
        日本5R型曲軸連桿頸磨床修整器的維修
        日本边添边摸边做边爱| av少妇偷窃癖在线观看| av狠狠色丁香婷婷综合久久| 日韩内射美女片在线观看网站| 国产av综合网站不卡| a黄片在线视频免费播放| 视频一区视频二区自拍偷拍 | 国产性感丝袜在线观看| 国产三级精品和三级男人| 国产三级不卡视频在线观看 | 精品国产亚洲级一区二区| 亚洲一区二区三区日本久久九| 婷婷久久精品国产色蜜蜜麻豆| 精人妻无码一区二区三区| 最爽无遮挡行房视频| 人妻 色综合网站| 公粗挺进了我的密道在线播放贝壳 | 亚洲综合小综合中文字幕| 国产一区二区三区免费在线播放| 国产内射视频在线观看| 国产男女猛烈无遮挡免费视频网址| 久久精品有码中文字幕1| 久久精品国产亚洲av热明星| 久久亚洲综合亚洲综合| 久久精品av在线观看| 色五月丁香五月综合五月| 无码人妻精品一区二区三区9厂| 日本大片免费观看视频| 久久成人国产精品免费软件| 亚洲国产成人久久一区| 国产精品免费观看久久| 久久伊人色av天堂九九| 色一情一区二| 99久久人妻无码精品系列蜜桃| 精品国产av无码一道| 亚洲乱精品中文字字幕| 精品黑人一区二区三区| 一区二区三区视频在线免费观看 | 亚洲无码观看a| 国产精品自产拍av在线| 日韩av在线亚洲女同|