• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮時(shí)變磨損的機(jī)器人磨拋參數(shù)優(yōu)化研究

    2023-02-15 01:25:04鄭乾健肖聚亮劉海濤
    關(guān)鍵詞:砂紙磨損建模

    王?超,鄭乾健,肖聚亮,周?宇,劉海濤,黃?田

    考慮時(shí)變磨損的機(jī)器人磨拋參數(shù)優(yōu)化研究

    王?超1, 2,鄭乾健1,肖聚亮1,周?宇1,劉海濤1,黃?田1

    (1. 天津大學(xué)機(jī)械工程學(xué)院,天津 300354;2. 濰柴動(dòng)力股份有限公司工藝工匠研究院,濰坊 261061)

    針對(duì)自動(dòng)化磨拋過程中存在的磨拋工具磨損現(xiàn)象,使得難以建立精準(zhǔn)的磨拋去除模型進(jìn)而達(dá)到穩(wěn)定一致的磨拋效果這一問題,提出了一種考慮時(shí)變磨損的機(jī)器人磨拋參數(shù)優(yōu)化方法.首先,將物理建模方法中的豐富信息和回歸建模方法的跟蹤學(xué)習(xí)特性相結(jié)合,提出了一種融合先驗(yàn)知識(shí)的磨拋材料去除回歸模型,極大地減少了原本回歸模型所需的大量實(shí)驗(yàn)數(shù)據(jù),使得提出的模型能夠跟蹤砂紙的磨損變化;其次,將所使用模型的預(yù)測(cè)結(jié)果與其他文獻(xiàn)所提出的理論模型以及基于支持向量回歸的預(yù)測(cè)模型進(jìn)行了對(duì)比分析;最后,使用所得的預(yù)測(cè)模型對(duì)機(jī)器人磨拋參數(shù)進(jìn)行優(yōu)化,在混聯(lián)機(jī)器人拋光實(shí)驗(yàn)臺(tái)上進(jìn)行實(shí)驗(yàn)驗(yàn)證,結(jié)果表明,所提出的融合先驗(yàn)知識(shí)的磨拋材料去除模型相比于常用模型具有更高的預(yù)測(cè)精度和更少的實(shí)驗(yàn)數(shù)據(jù),所提出的機(jī)器人磨拋參數(shù)優(yōu)化方法能夠有效地補(bǔ)償磨拋工具的磨損,保證磨拋過程中材料去除的一致性,同時(shí)提高了磨拋加工效率.

    時(shí)變磨損;支持向量回歸;材料去除模型;參數(shù)優(yōu)化

    模具磨拋精加工在整個(gè)制造所費(fèi)時(shí)間中的占比較大,約占37%[1].其目的是為了降低工件表面的粗糙度,同時(shí)在一定程度上去除模具成型時(shí)形成的凹凸層及裂紋層,即刀紋,并修復(fù)工件表面的面型誤差,以期獲得所要求的形狀、尺寸、位置精度和表面粗糙度.傳統(tǒng)的模具打磨和拋光主要是依賴于熟練工人的手工加工,是一種低效率重復(fù)性的勞動(dòng),且隨機(jī)性較大,很難保證加工質(zhì)量的穩(wěn)定性和一致性[2].此外,磨拋?zhàn)鳂I(yè)的環(huán)境噪音大,工人暴露在磨拋加工的粉塵中,對(duì)其健康也會(huì)造成嚴(yán)重的影響[3].因此開展磨拋?zhàn)詣?dòng)化的研究是具有重要意義的.

    在自動(dòng)化磨拋過程中,機(jī)器人可以使用多種拋光工具對(duì)不同類型的工件進(jìn)行拋光,通過調(diào)控拋光力、主軸轉(zhuǎn)速、進(jìn)給速度等參數(shù)進(jìn)行穩(wěn)定和均勻的去除,但是如何正確綜合地處理這些因素以達(dá)到在保證加工效果的情況下效率最大化始終是一個(gè)難點(diǎn),這需要建立準(zhǔn)確的磨拋去除模型來作為綜合調(diào)控的依據(jù).此外,由于加工過程中磨拋工具存在不可避免的磨損現(xiàn)象,其磨損程度對(duì)磨拋去除量的影響很大,其動(dòng)態(tài)變化的特性對(duì)建立準(zhǔn)確的磨拋去除模型進(jìn)而達(dá)到恒定均勻的磨拋效果造成了很大的阻礙.目前主要有兩類磨拋建模方法:一為物理建模方法;二為回歸建模方法.

    回歸建模方法指的是通過選擇可觀測(cè)的關(guān)鍵影響因素作為輸入,選擇合適的建模指標(biāo)作為輸出,利用數(shù)學(xué)統(tǒng)計(jì)和機(jī)器學(xué)習(xí)的方法訓(xùn)練磨拋的數(shù)學(xué)模型.國(guó)內(nèi)外在這方面的研究較少.浙江大學(xué)馮海濤[8]以砂帶速度、機(jī)器人進(jìn)給速度和工件曲率為輸入,以加工后的粗糙度為因變量,建立了基于支持向量回歸的粗糙度預(yù)測(cè)模型,根據(jù)不同的拋光參數(shù)對(duì)加工后的表面粗糙度進(jìn)行分析和預(yù)測(cè),但是該模型并沒有考慮磨拋工具磨損的影響,且沒有從材料去除的角度考慮.Yue等[9]提出了一種基于支持向量回歸的磨削深度預(yù)測(cè)模型,其模型根據(jù)現(xiàn)場(chǎng)測(cè)量數(shù)據(jù)對(duì)模型進(jìn)行修正,能夠在一定程度上適應(yīng)環(huán)境的變化如砂帶的逐漸磨損、新的待加工工件等,但是為了得到較佳的模型,實(shí)驗(yàn)需要收集大量訓(xùn)練用數(shù)據(jù),實(shí)驗(yàn)所需時(shí)間和花費(fèi)較大.Pandiyan等[10]使用田口實(shí)驗(yàn)確定的磨拋參數(shù)對(duì)不同的統(tǒng)計(jì)回歸方法進(jìn)行評(píng)估,得出砂粒尺寸對(duì)材料去除的影響最大,但是沒有考慮到磨拋過程中砂粒尺寸會(huì)因磨損而變化.相較于物理建模的方法而言,使用回歸建模是一種更能有效評(píng)估磨損變化的建模方法,其能夠利用實(shí)驗(yàn)數(shù)據(jù),抓住主要變量.但是現(xiàn)有的回歸建模方法所采用的“黑箱子”式的算法,基本沒有利用已有物理理論模型的物理、幾何知識(shí),這導(dǎo)致回歸模型輸入的維度較大(包括主軸轉(zhuǎn)速、進(jìn)給速度、磨拋壓力、磨拋角度、磨拋工具半徑、工件表面曲率、磨拋時(shí)間等),為了得到較為精準(zhǔn)的去除模型,每個(gè)維度的數(shù)據(jù)必須要有一定的密度,因此回歸所需的數(shù)據(jù)量呈指數(shù)級(jí)上升.而在實(shí)際的工作過程中,由于需要測(cè)得的數(shù)據(jù)點(diǎn)和數(shù)據(jù)類型較多,所需的實(shí)驗(yàn)時(shí)間是非常巨大的,這顯然與實(shí)際制造加工的要求不符.清華大學(xué)呂洪波等[11]創(chuàng)新性地提出了融合先驗(yàn)知識(shí)的建模方法,基于半經(jīng)驗(yàn)公式生成虛擬樣本從而解決回歸建模中所需樣本數(shù)量過多的問題,但是在該文中將半經(jīng)驗(yàn)公式中的磨削系數(shù)以及磨削壓力、工件進(jìn)給速度、砂帶速度的次數(shù)均作為全局系數(shù),即在整個(gè)磨削過程中保持不變,其適應(yīng)性較差.

    針對(duì)上述兩種磨拋模型建模方法的優(yōu)缺點(diǎn),并且受文獻(xiàn)[11]的啟發(fā),本文提出了一種融合先驗(yàn)知識(shí)的材料去除回歸模型,將基于機(jī)理的物理建模方法包含的大量關(guān)于磨拋問題的信息利用起來,極大地減少了回歸建模方法所需要的實(shí)驗(yàn)數(shù)據(jù)量;將回歸建模方法中使用的機(jī)器學(xué)習(xí)方法跟蹤磨拋工具磨損程度變化的優(yōu)勢(shì)利用起來,對(duì)理論模型中的Preston參數(shù)進(jìn)行跟蹤學(xué)習(xí),使得最終的模型能夠反映磨損的變化過程.最后,根據(jù)所提出的融合先驗(yàn)知識(shí)的材料去除回歸模型,建立了一種機(jī)器人磨拋參數(shù)優(yōu)化框架,對(duì)磨拋參數(shù)進(jìn)行了優(yōu)化,使得磨拋工具徹底失效之前盡可能加工更多的表面,并達(dá)到較為一致的拋光效果和材料去除深度.

    1?機(jī)器人磨拋參數(shù)優(yōu)化框架

    本文構(gòu)建的機(jī)器人磨拋參數(shù)優(yōu)化框架主要包括3部分:數(shù)據(jù)收集、模型建立、參數(shù)優(yōu)化實(shí)驗(yàn).框架如圖1所示.

    圖1?機(jī)器人磨拋參數(shù)優(yōu)化框架

    在數(shù)據(jù)收集部分,由于加工工件的磨拋質(zhì)量受到多種因素影響,因此需要對(duì)磨拋模型中可控可觀測(cè)的參數(shù)進(jìn)行獲取,如主軸轉(zhuǎn)速、機(jī)器人加工進(jìn)給速度、磨拋力等.此外為了評(píng)估磨拋效果,需要對(duì)工件表面材料去除深度、砂紙表面形貌信息進(jìn)行采集,對(duì)被加工表面的表面粗糙度進(jìn)行測(cè)量.

    在模型建立部分,把磨拋去除理論模型中磨拋參數(shù)間的關(guān)系作為先驗(yàn)知識(shí),使用支持向量回歸的方法對(duì)理論模型中的Preston系數(shù)進(jìn)行回歸,實(shí)現(xiàn)在較少實(shí)驗(yàn)數(shù)據(jù)量下能夠精準(zhǔn)地預(yù)測(cè)磨拋的材料去除量.

    在參數(shù)優(yōu)化部分,本文針對(duì)進(jìn)給速度進(jìn)行優(yōu)化,保持磨拋壓力、主軸轉(zhuǎn)速等參數(shù)為恒定,結(jié)合混聯(lián)機(jī)器人連桿模具拋光實(shí)驗(yàn)臺(tái)進(jìn)行實(shí)驗(yàn),將其與參數(shù)優(yōu)化前的樣件進(jìn)行對(duì)比,驗(yàn)證了本文所提出的機(jī)器人磨拋參數(shù)優(yōu)化方法的有效性.

    2?融合先驗(yàn)知識(shí)的材料去除回歸模型

    考慮時(shí)變磨損的材料去除回歸模型將物理理論模型與回歸模型的優(yōu)點(diǎn)相結(jié)合,將理論材料去除模型中的常數(shù)項(xiàng)Preston系數(shù)轉(zhuǎn)變?yōu)殡S磨拋過程進(jìn)行而變化的時(shí)變函數(shù),把理論模型中磨拋力、進(jìn)給速度、主軸轉(zhuǎn)速等內(nèi)在關(guān)系作為回歸模型的先驗(yàn)知識(shí),使用回歸的方法量化在磨拋進(jìn)程中磨拋工具磨損的影響,進(jìn)而準(zhǔn)確計(jì)算磨拋進(jìn)程中的材料去除量,為機(jī)器人磨拋參數(shù)優(yōu)化打下基礎(chǔ).

    2.1?模型建立基礎(chǔ)

    Preston方程于1927年提出[12],被廣泛應(yīng)用在磨拋加工中,它將磨拋去除深度描述為

    由于實(shí)際機(jī)器人控制拋光實(shí)驗(yàn)采用控制拋光盤沿著拋光軌跡連續(xù)進(jìn)給的方式,因此采用速度駐留形式的普雷斯頓方程,將式(1)改寫為沿拋光路徑的微分形式,即

    式中f表示拋光盤沿軌跡的進(jìn)給速度.

    根據(jù)文獻(xiàn)[13]所做的研究,如圖2所示,當(dāng)拋光盤沿直線在平面上運(yùn)動(dòng),進(jìn)給方向?yàn)榉较颍瑨伖鈪^(qū)域內(nèi)一點(diǎn)(,)的壓強(qiáng)為

    圖2?拋光盤沿直線進(jìn)給時(shí)的速度分布

    式中:n表示拋光盤的法向正壓力;p表示拋光盤的半徑.

    (,)的相對(duì)線速度分布為

    由于機(jī)器人的磨拋加工是一個(gè)累積過程,拋光盤在拋光接觸區(qū)域某一點(diǎn)(,)的磨拋加工可以看作是從2點(diǎn)到1點(diǎn)的累加,因此可以得到拋光盤沿軌跡的普雷斯頓方程,即

    通過實(shí)驗(yàn)對(duì)Preston系數(shù)p的測(cè)定,能夠得到新磨拋工具較為精準(zhǔn)的磨拋模型.但是在實(shí)際磨拋進(jìn)程中,磨拋工具會(huì)快速磨損直至失效,隨著磨拋工具表面形貌的迅速變化,Preston系數(shù)p也會(huì)隨之變化,使得磨拋模型失效.故需要從方法的角度跟蹤學(xué)習(xí)磨拋磨損這一變化過程,量化磨拋工具隨時(shí)間磨損的程度.本文采用支持向量回歸的方法對(duì)其進(jìn)行量化和建模.

    2.2?模型建立

    2.2.1?支持向量回歸算法原理

    支持向量機(jī)(support vector machine,SVM)算法最早是由Cortes等[14]提出的最大化間隔的分類算法,當(dāng)其用于回歸問題時(shí)稱作支持向量回歸(support vector regression,SVR),由Drucker等[15]提出.機(jī)器人磨拋去除模型建模是一個(gè)非線性回歸問題,對(duì)于非線性回歸,SVR采用的方式是先把樣本點(diǎn)映射到高維空間,然后用高維空間中的線性超平面進(jìn)行近似.對(duì)于給定數(shù)據(jù)

    SVR問題的優(yōu)化目標(biāo)為

    其限制條件為

    2.2.2?量化磨損及磨拋深度預(yù)測(cè)

    為了獲取模型構(gòu)建所需要的實(shí)驗(yàn)數(shù)據(jù),在混聯(lián)機(jī)器人拋光實(shí)驗(yàn)臺(tái)進(jìn)行實(shí)驗(yàn),采用Nanovea 白光干涉三維輪廓儀對(duì)試件輪廓進(jìn)行采集.如圖3(a)所示,由于受到輪廓儀的采集限制,實(shí)驗(yàn)試件設(shè)計(jì)為長(zhǎng)寬均為80mm、材質(zhì)為5CrNiMo的平面鋼板.如圖3(b)所示,在實(shí)際實(shí)驗(yàn)過程中,為保證磨拋過程的連續(xù)性,將3塊試件并排放置進(jìn)行實(shí)驗(yàn).

    圖3?實(shí)驗(yàn)試件與安裝

    如圖4所示為實(shí)驗(yàn)流程框圖,具體包括以下5個(gè)步驟.

    步驟1 針對(duì)并排固定的3塊試件做如圖5所示的軌跡規(guī)劃,圖5中1為磨拋軌跡的起始點(diǎn),2為磨拋軌跡的終點(diǎn),箭頭方向即為磨拋進(jìn)給方向.

    步驟2 使用Nanovea 白光干涉三維輪廓儀獲取砂紙初始表面反射光強(qiáng)度分布,該分布可以作為評(píng)價(jià)砂紙磨損程度的一個(gè)指標(biāo)[16],同時(shí)使用三維輪廓儀掃描獲取工件表面的初始輪廓高度.3塊鋼板的測(cè)量范圍分別為圖5中的點(diǎn)至點(diǎn)、點(diǎn)至點(diǎn)、點(diǎn)至點(diǎn).

    步驟3 在單次實(shí)驗(yàn)中,同一次磨拋保持法向磨拋力和進(jìn)給速度恒定.表1表示各組實(shí)驗(yàn)數(shù)據(jù)中相同的參數(shù)或設(shè)置.表2為實(shí)驗(yàn)數(shù)據(jù)記錄樣例.

    圖4?實(shí)驗(yàn)流程

    圖5?實(shí)驗(yàn)軌跡規(guī)劃

    表1?各組實(shí)驗(yàn)相同的參數(shù)

    Tab.1?Same parameters in each group of experiments

    表2?實(shí)驗(yàn)數(shù)據(jù)樣例

    Tab.2?Sample experimental data

    圖6?試件磨拋前后表面形貌

    在訓(xùn)練SVR回歸模型前,將式(10)作為先驗(yàn)知識(shí)輸入至模型中,使得磨拋去除深度hp與時(shí)間的變化關(guān)系轉(zhuǎn)化為Preston系數(shù)與時(shí)間的變化關(guān)系,進(jìn)而簡(jiǎn)化了問題.同時(shí),由于該回歸問題為關(guān)于時(shí)間序列的預(yù)測(cè)問題,因此不采用常用的K折交叉驗(yàn)證方法,而采用擴(kuò)展窗口交叉驗(yàn)證方法尋找最優(yōu)的超參數(shù)組合[17].當(dāng)實(shí)際測(cè)量計(jì)算得到的磨拋去除深度與SVR模型回歸預(yù)測(cè)得到的去除深度兩者之間的差值在區(qū)間內(nèi),則認(rèn)為SVR模型在該點(diǎn)的預(yù)測(cè)是正確的.圖7中的訓(xùn)練準(zhǔn)確率指的是數(shù)據(jù)集中SVR預(yù)測(cè)正確的點(diǎn)數(shù)量在數(shù)據(jù)集中所占的比重.考慮到由于測(cè)量時(shí)試件存在污漬或測(cè)量設(shè)備振動(dòng)等一些其他情況,使得實(shí)驗(yàn)會(huì)得到明顯偏離理想值的一些“壞點(diǎn)”,因此當(dāng)準(zhǔn)確率大于90%時(shí),即認(rèn)為訓(xùn)練得到的模型是可行的.

    圖8?Preston系數(shù)變化曲線

    表3?SVR參數(shù)

    Tab.3?SVR parameters

    圖9?一定工況下的材料去除深度預(yù)測(cè)

    此外,通過多組磨拋參數(shù)組合實(shí)驗(yàn)發(fā)現(xiàn),當(dāng)磨拋壓力較低(10~25N)時(shí),砂紙的磨損速度很慢,但是磨拋的效果也很差,當(dāng)磨拋壓力提高到一定程度(30~50N)時(shí),由于在高磨拋力下砂紙的表面磨??焖倌p,砂紙快速進(jìn)入磨損穩(wěn)定期,并且磨損速率趨于一致,與文獻(xiàn)[18]所得出的結(jié)論是一致的.進(jìn)給速度在1~4mm/s的范圍內(nèi)時(shí),其變化對(duì)砂紙磨損速率的影響非常小,與文獻(xiàn)[19]通過對(duì)不同磨拋參數(shù)的相對(duì)靈敏度分析得出砂帶磨損速率對(duì)進(jìn)給速度變化的敏感度最低的結(jié)論相一致.主軸轉(zhuǎn)速的變化對(duì)砂紙的磨損速度有較大的影響,當(dāng)主軸轉(zhuǎn)速過大時(shí),砂紙將快速磨損至失效,難以持久作業(yè);當(dāng)主軸轉(zhuǎn)速過小時(shí),雖然磨拋盤的磨損減緩,但是磨拋效果較差.綜上,可以通過對(duì)進(jìn)給速度與磨拋壓力在一定范圍內(nèi)的調(diào)整優(yōu)化來補(bǔ)償磨拋盤的磨損.同時(shí),為達(dá)到更佳的磨拋效果,也需要確定較佳的主軸轉(zhuǎn)速.

    2.2.3?模型對(duì)比分析

    為了確定本文方法在預(yù)測(cè)效果良好的前提下所需要的訓(xùn)練樣本數(shù)量,使用不同訓(xùn)練樣本數(shù)量對(duì)本文方法,即融合先驗(yàn)知識(shí)的SVR方法進(jìn)行訓(xùn)練,每50個(gè)數(shù)據(jù)點(diǎn)為一組,訓(xùn)練樣本數(shù)分別為150、300、450、600.從除訓(xùn)練集外的5組實(shí)驗(yàn)組中隨機(jī)抽取40個(gè)數(shù)據(jù)點(diǎn)作為測(cè)試集(每組抽取8個(gè)數(shù)據(jù)點(diǎn)),使用訓(xùn)練出來的模型進(jìn)行預(yù)測(cè),在結(jié)果輸出中,先將測(cè)試集按照實(shí)測(cè)值從小到大進(jìn)行排列.如圖10所示為使用不同訓(xùn)練樣本數(shù)量訓(xùn)練得到的模型在測(cè)試樣本中的預(yù)測(cè)值與實(shí)測(cè)值的對(duì)比.從圖10中可以看出,訓(xùn)練樣本數(shù)量為150的模型其預(yù)測(cè)效果已經(jīng)較為優(yōu)秀,隨著訓(xùn)練樣本數(shù)量的增加,幾種模型的預(yù)測(cè)性能相差無幾.考慮到減少實(shí)驗(yàn)成本和時(shí)間,使用樣本數(shù)量為150的樣本集即可滿足預(yù)測(cè)效果,之后使用該模型與其他算法進(jìn)行比對(duì).

    圖10  本文方法在不同訓(xùn)練樣本下的預(yù)測(cè)值與實(shí)測(cè)值的對(duì)比

    圖11?預(yù)測(cè)結(jié)果與實(shí)測(cè)值結(jié)果對(duì)比

    表4?誤差分析結(jié)果對(duì)比

    Tab.4?Error analysis and comparison results

    根據(jù)圖11與表4可以看出,理論模型由于其Preston計(jì)算方法以及未考慮磨損的原因,其大部分預(yù)測(cè)值相較于實(shí)測(cè)值都較大,偏離程度取決于磨拋加工時(shí)間的長(zhǎng)短;訓(xùn)練樣本數(shù)量為600的原SVR方法在一定程度上也能夠預(yù)測(cè)出磨拋去除深度,但是在部分?jǐn)?shù)據(jù)點(diǎn)上的預(yù)測(cè)值與實(shí)測(cè)值有較大偏離,本文方法的SSE、RMSE、MAPE、MAE相較于使用訓(xùn)練樣本數(shù)量為600的原SVR方法分別降低了0.6652μm2、0.0712μm、1.36%、0.0424μm,均是本文所提出的融合先驗(yàn)知識(shí)的SVR方法更為優(yōu)秀,分析其原因在于原SVR方法的輸入矩陣維度較大,而訓(xùn)練樣本的數(shù)量?jī)H有600,樣本對(duì)于磨拋參數(shù)的分布而言是較為稀疏的,訓(xùn)練樣本數(shù)不夠,導(dǎo)致原SVR方法對(duì)于未知樣本的預(yù)測(cè)能力較差.而當(dāng)訓(xùn)練樣本數(shù)增至1200時(shí),通過表4可以發(fā)現(xiàn)原SVR方法與本文所提出的融合先驗(yàn)知識(shí)的SVR方法的誤差評(píng)價(jià)指標(biāo)相差不大,綜合考慮測(cè)量誤差、加工誤差等其他因素,可以認(rèn)為兩者的預(yù)測(cè)能力基本一致,但是本文方法訓(xùn)練樣本數(shù)量?jī)H為150,樣本數(shù)量遠(yuǎn)少于前者.綜上,融合先驗(yàn)知識(shí)的磨拋材料去除模型能夠較好地量化評(píng)估磨拋工具的磨損并得到較佳的磨拋深度預(yù)測(cè)效果,在小樣本預(yù)測(cè)方面具有優(yōu)越性,大大減少了實(shí)驗(yàn)所需要的時(shí)間和經(jīng)濟(jì)成本.

    3?考慮時(shí)變磨損的機(jī)器人磨拋參數(shù)優(yōu)化實(shí)驗(yàn)

    為保證材料去除的一致性以及提高磨拋的效率,需要通過改進(jìn)磨拋參數(shù)對(duì)磨拋工具的磨損進(jìn)行補(bǔ)償.較為常用的可控可觀測(cè)磨拋參數(shù)為磨拋進(jìn)給速度f和磨拋壓力n,可以只調(diào)控某一磨拋參數(shù)如進(jìn)給速度f(或磨拋壓力n),也可以為了實(shí)現(xiàn)磨拋時(shí)間最短、磨拋參數(shù)的調(diào)控盡可能平滑等目標(biāo),同時(shí)調(diào)控進(jìn)給速度f和磨拋壓力n.為了簡(jiǎn)單起見,將磨拋壓力n設(shè)為常數(shù),只調(diào)控磨拋進(jìn)給速度f來保證較好的磨拋去除效果.設(shè)定磨拋深度為3μm,給定法向力為40N.為了便于工程實(shí)際應(yīng)用和計(jì)算,使用圖8中折線2所示的p變換關(guān)系,通過反算進(jìn)給速度f可以得到如圖12所示的進(jìn)給速度優(yōu)化曲線.

    圖12?進(jìn)給速度優(yōu)化曲線

    使用優(yōu)化后的磨拋參數(shù)在混聯(lián)機(jī)器人拋光實(shí)驗(yàn)臺(tái)進(jìn)行實(shí)驗(yàn),如圖13(a)所示.待加工工件表面的微觀形貌,如圖13(b)所示,表面可見有明顯密集的刀紋. 經(jīng)過參數(shù)調(diào)控磨拋后工件表面的微觀形貌,如圖13(c)所示,幾乎不可辨加工方向.

    同樣使用白光干涉三維輪廓儀對(duì)試件表面進(jìn)行檢測(cè)獲取磨拋過程中的去除深度.若使用文獻(xiàn)[12]所述理論模型,當(dāng)設(shè)定磨拋深度為3μm,給定法向力為40N,則可計(jì)算得到磨拋過程中的進(jìn)給速度為2.4mm/s,為一定值.使用理論模型未經(jīng)過優(yōu)化的進(jìn)給速度所測(cè)得的磨拋去除深度以及參數(shù)優(yōu)化后測(cè)得的磨拋去除深度與期望值的對(duì)比,如圖14所示,從圖14中可以看出,未經(jīng)過優(yōu)化的進(jìn)給速度所測(cè)得的磨拋去除深度在磨拋初期略大于3μm,而隨著磨拋工具的磨損,其去除深度逐漸減小并且始終小于3μm,無法在全過程中保持材料去除的一致性.而經(jīng)過參數(shù)優(yōu)化后的實(shí)驗(yàn)組其磨拋去除深度在全過程中與期望值較為一致,能夠較好地補(bǔ)償砂紙的磨損,保證材料去除量的一致性.

    圖13?磨拋前后工件表面微觀形貌對(duì)比

    圖14?參數(shù)優(yōu)化前后磨拋去除深度與期望值對(duì)比

    此外使用北京時(shí)代TR200高精度表面粗糙度儀分別測(cè)量磨拋時(shí)間為10s、20s、30s、120s、130s、140s、230s、240s、250s處試件的表面粗糙度,圖15所示為磨拋參數(shù)優(yōu)化前后在這9個(gè)采樣點(diǎn)處表面粗糙度的大小,從圖15中可以看出,經(jīng)過參數(shù)優(yōu)化后的實(shí)驗(yàn)組雖然在磨拋初期其表面粗糙度比未經(jīng)過參數(shù)優(yōu)化的實(shí)驗(yàn)組更大,但是其磨拋效果在全過程中較為一致,整體上磨拋效果更佳,在滿足表面粗糙度要求的情況下能夠加工更長(zhǎng)的時(shí)間,而未經(jīng)參數(shù)優(yōu)化的實(shí)驗(yàn)組其磨拋效果由于砂紙磨損而快速下降,在大部分時(shí)間的磨拋效果均劣于經(jīng)過參數(shù)優(yōu)化后的實(shí)驗(yàn)組,在230s后的表面粗糙度明顯大于參數(shù)優(yōu)化后的實(shí)驗(yàn)組.

    圖15?參數(shù)優(yōu)化前后表面粗糙度對(duì)比

    綜上,從磨拋去除深度和表面粗糙度兩個(gè)角度證明了本文所提出的磨拋參數(shù)優(yōu)化方法能夠較好地補(bǔ)償砂紙磨損,進(jìn)而提高砂紙的利用率和提高磨拋效率,是一種切實(shí)有效的機(jī)器人磨拋參數(shù)優(yōu)化方法.

    4?結(jié)?語

    本文提出了一種考慮時(shí)變磨損的機(jī)器人磨拋參數(shù)優(yōu)化方法.首先將理論模型中材料去除深度與磨拋參數(shù)之間的關(guān)系作為先驗(yàn)知識(shí),融合回歸模型能夠跟蹤學(xué)習(xí)磨拋磨損變化情況的特性,提出了一種融合先驗(yàn)知識(shí)的磨拋材料去除模型,極大地減少了回歸模型所需的數(shù)據(jù)維度和數(shù)據(jù)數(shù)量,方便工程實(shí)際應(yīng)用.通過與傳統(tǒng)SVR方法與傳統(tǒng)理論模型進(jìn)行對(duì)比分析證明了融合先驗(yàn)知識(shí)的SVR模型能夠用來預(yù)測(cè)磨拋去除深度,并且在小樣本預(yù)測(cè)上具有顯著優(yōu)勢(shì).根據(jù)提出的去除模型優(yōu)化機(jī)器人磨拋參數(shù),并在混聯(lián)機(jī)器人拋光實(shí)驗(yàn)臺(tái)進(jìn)行實(shí)驗(yàn),實(shí)驗(yàn)證明所提出的考慮時(shí)變磨損的機(jī)器人磨拋參數(shù)優(yōu)化方法能夠較好地補(bǔ)償磨拋工具的磨損,使得磨拋去除量在全加工過程中較為一致,磨拋加工效果更佳,此外在滿足表面粗糙度要求的情況下能夠加工更長(zhǎng)的時(shí)間,進(jìn)而提高了磨拋工具的利用率和磨拋效率.

    [1] Saito K. Finishing and polishing of free-form surface[J]. Bulletin of the Japan Society of Precision Engineering,1984,18(2):104-109.

    [2] Kamijyo K,Sakamoto M,Asao T,et al. Development of an automatic polishing system for metal molds and dies(1st report,automation of polishing operation based on CAD/CAM data)[J]. Transactions of the Japan Society of Mechanical Engineers Series C,1989,55(509):193-196.

    [3] Mohammad A,Wang D. Electrochemical mechanical polishing technology:Recent developments and future research and industrial needs[J]. International Journal of Advanced Manufacturing Technology,2016,86:1-16.

    [4] 曾?晰. 軟固結(jié)磨粒氣壓砂輪設(shè)計(jì)方法及材料去除特性研究[D]. 杭州:浙江工業(yè)大學(xué)機(jī)械工程學(xué)院,2013.

    Zeng Xi. Research on the Design Method and Material Removal Characteristics of the Soft-Consolidated Abra-sive Grain Pneumatic Grinding Wheel[D]. Hangzhou:School of Mechanical Engineering,Zhejiang University of Technology,2013(in Chinese).

    [5] Cao Z C,Cheung C F. Multi-scale modeling and simulation of material removal characteristics in computer-controlled bonnet polishing[J]. International Journal of Mechanical Sciences,2016,106:147-156.

    [6] 朱天寶,金曉怡,陳志鵬,等. 復(fù)雜曲面柔性拋光材料去除模型研究[J]. 組合機(jī)床與自動(dòng)化加工技術(shù),2020(9):76-79.

    Zhu Tianbao,Jin Xiaoyi,Chen Zhipeng,et al. Study on material removal model for flexible polishing of complex surfaces[J]. Modular Machine Tools and Automatic Machining Technology,2020(9):76-79(in Chinese).

    [7] Huang H,Gong Z M,Chen X Q,et al. Robotic grinding and polishing for turbine-vane overhaul[J]. Journal of Materials Processing Technology,2002,127(2):140-145.

    [8] 馮海濤. 機(jī)器人自動(dòng)化拋光系統(tǒng)關(guān)鍵技術(shù)的研究[D]. 杭州:浙江大學(xué)機(jī)械工程學(xué)院,2015.

    Feng Haitao. Research on Key Technologies of Robot Automatic Polishing System[D]. Hangzhou:School of Mechanical Engineering,Zhejiang University,2015(in Chinese).

    [9] Yue Y,Zhang J,Zhou Y,et al. Inverse input predic-tion model for robotic belt grinding[J]. International Journal of Intelligent Robotics and Applications,2021,5(4):465-476.

    [10]Pandiyan V,Caesarendra W,Glowacz A,et al. Modelling of material removal in abrasive belt grinding process:A regression approach[J]. Symmetry,2020,12(1):363-380.

    [11]呂洪波,宋亦旭,賈培發(fā). 機(jī)器人修磨中融合先驗(yàn)知識(shí)的適應(yīng)學(xué)習(xí)建模方法[J]. 機(jī)器人,2011,33(6):641-648.

    Lü Hongbo,Song Yixu,Jia Peifa. Adaptive learning modeling method incorporating prior knowledge in robot grinding[J]. Robot,2011,33(6):641-648(in Chi-nese).

    [12] Preston F W. The theory and design of glass plate polishing machine[J]. Journal of the Society of Glass Technology,1927,11:214-256.

    [13]樊?成,趙?繼,張?雷,等. 移動(dòng)拋光自由曲面材料去除的理論建模與試驗(yàn)研究[J]. 機(jī)械工程學(xué)報(bào),2014,50(5):173-181.

    Fan Cheng,Zhao Ji,Zhang Lei,et al. Theoretical modeling and experimental research on material removal of free-form surface in mobile polishing[J]. Journal of Mechanical Engineering,2014,50(5):173-181(in Chinese).

    [14] Cortes C,Vapnik V. Support vector networks[J]. Machine Learning,1995,20(3):273-297.

    [15] Drucker H,Burges C,Kaufman L,et al. Support vector regression machines[J]. Advances in Neural Information Processing Systems,1997,28(7):779-784.

    [16] He Z,Li J Y,Wu Y F,et al. Quantitative study on abrasive belt wear based on geometric parameters[J]. Key Engineering Materials,2018,764:156-163.

    [17]Bergmeir C,José M B. On the use of cross-validation for time series predictor evaluation[J]. Information Sciences,2012,191:192-213.

    [18]Yang C,Huang Y,Huang Z,et al. Experimental research on the abrasive belt grinding 1Cr13 stainless steel[J]. Mechanical ence and Technology for Aerospace Engineering,2011,30(6):1007-1010,1015.

    [19] Lin X,Xin X,Shan X,et al. Optimal parameter ranges of material removal depth of abrasive cloth wheel polishing based on sensitivity analysis[J]. The International Journal of Advanced Manufacturing Technology,2019,105(12):5165-5179.

    Research on Optimization of Robot Grinding and Polishing Parameters Considering Time-Varying Wear

    Wang Chao1, 2,Zheng Qianjian1,Xiao Juliang1,Zhou Yu1,Liu Haitao1,Huang Tian1

    (1. School of Mechanical Engineering,Tianjin University,Tianjin 300354,China;2. Craftsman Research Institute,Weichai Power Co.,Ltd.,Weifang 261061,China)

    Addressing the problem of grinding and polishing tool wear in the automatic grinding and polishing process,an issue that makes it difficult to establish an accurate grinding and polishing removal model to achieve a stable and consistent grinding and polishing effect,a robot grinding and polishing parameter optimization method considering time-varying wear is proposed. First of all,combining the rich information in the physical modeling method with the tracking learning characteristics of the regression modeling method,a regression model for grinding and polishing material removal that integrates prior knowledge is proposed. Such a model will greatly reduce the large amount of experimental data required. The proposed model can track the changes in grinding and polishing wear. The prediction results of the model used are compared with theoretical models and the support vector regression prediction model proposed in other research. Finally,the prediction model obtained is used to optimize the robot grinding and polishing parameters,and experimental verification is carried out on the polishing experiment platform of the hybrid robot. The results show that the proposed grinding and polishing material removal model fused with prior knowledge has higher prediction accuracy and less experimental data. The proposed optimization method for robot grinding and polishing parameters can effectively compensate for the wear of grinding and polishing tools,ensure the consistency of material removal in the grinding and polishing process,and improve the efficiency of grinding and polishing.

    time-varying wear;support vector regression;material removal model;parameter optimization

    10.11784/tdxbz202110006

    TH162

    A

    0493-2137(2023)02-0127-10

    2021-10-11;

    2021-12-02.

    王?超(1985—??),男,博士研究生,wangchao4585@163.com.

    鄭乾健,2020201185@tju.edu.cn.

    國(guó)家科技重大專項(xiàng)資助項(xiàng)目(2017ZX04021001-004).

    Supported by the National Science and Technology Major Project(No. 2017ZX04021001-004).

    (責(zé)任編輯:王曉燕)

    猜你喜歡
    砂紙磨損建模
    聯(lián)想等效,拓展建模——以“帶電小球在等效場(chǎng)中做圓周運(yùn)動(dòng)”為例
    一張砂紙
    天津詩人(2020年2期)2020-11-18 06:24:28
    青海湖(2020年6期)2020-09-07 08:08:22
    砂紙式復(fù)合片的防刺性能研究
    基于PSS/E的風(fēng)電場(chǎng)建模與動(dòng)態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    套管磨損機(jī)理研究
    不對(duì)稱半橋變換器的建模與仿真
    鉆頭磨損后怎樣更換
    普通砂紙有妙用
    三元組輻射場(chǎng)的建模與仿真
    在线观看午夜福利视频| 亚洲aⅴ乱码一区二区在线播放 | www.自偷自拍.com| 国产一区在线观看成人免费| 中文亚洲av片在线观看爽| 黄色片一级片一级黄色片| 纯流量卡能插随身wifi吗| 国产人伦9x9x在线观看| 91麻豆av在线| 十八禁人妻一区二区| 亚洲一码二码三码区别大吗| 免费观看人在逋| 十八禁网站免费在线| 99精品欧美一区二区三区四区| 男人舔女人下体高潮全视频| 亚洲美女黄片视频| 自线自在国产av| 久久婷婷成人综合色麻豆| 女性被躁到高潮视频| 我的亚洲天堂| 国产人伦9x9x在线观看| 久久伊人香网站| 无人区码免费观看不卡| 人人妻人人爽人人添夜夜欢视频| 亚洲熟女毛片儿| 黄色丝袜av网址大全| 免费不卡黄色视频| 黄片大片在线免费观看| 亚洲国产看品久久| 国产精品综合久久久久久久免费 | 欧美国产日韩亚洲一区| 高清毛片免费观看视频网站| 变态另类丝袜制服| 午夜福利免费观看在线| 日韩欧美三级三区| 成在线人永久免费视频| 桃红色精品国产亚洲av| 女人精品久久久久毛片| 九色国产91popny在线| 精品日产1卡2卡| 国产高清激情床上av| 美女大奶头视频| 精品熟女少妇八av免费久了| 99久久精品国产亚洲精品| 18禁观看日本| 18禁美女被吸乳视频| 日韩国内少妇激情av| 老鸭窝网址在线观看| 国产午夜精品久久久久久| 久久久久久久午夜电影| 亚洲人成电影免费在线| 国产成人啪精品午夜网站| 一级毛片高清免费大全| 男人舔女人的私密视频| 日韩欧美三级三区| 国产成人影院久久av| 非洲黑人性xxxx精品又粗又长| 熟女少妇亚洲综合色aaa.| 国产激情欧美一区二区| 久久香蕉精品热| 精品国产一区二区三区四区第35| 日本a在线网址| 精品国产一区二区久久| 久久 成人 亚洲| 亚洲国产高清在线一区二区三 | 大型av网站在线播放| 精品熟女少妇八av免费久了| 一级毛片精品| 12—13女人毛片做爰片一| 欧美午夜高清在线| 99riav亚洲国产免费| 成人亚洲精品av一区二区| 欧美日本视频| 久久久久亚洲av毛片大全| 天堂动漫精品| 国产欧美日韩一区二区精品| 精品国内亚洲2022精品成人| 国产高清视频在线播放一区| 亚洲成av片中文字幕在线观看| 亚洲五月天丁香| 亚洲欧美精品综合久久99| 亚洲国产欧美网| 一区二区三区高清视频在线| 99在线人妻在线中文字幕| 视频在线观看一区二区三区| 国产av一区在线观看免费| 亚洲午夜精品一区,二区,三区| 欧美日本亚洲视频在线播放| 久久人妻福利社区极品人妻图片| 成人国语在线视频| 免费在线观看视频国产中文字幕亚洲| 国产在线观看jvid| 日韩av在线大香蕉| 成人18禁在线播放| 亚洲人成网站在线播放欧美日韩| 亚洲第一青青草原| 亚洲色图 男人天堂 中文字幕| 亚洲欧洲精品一区二区精品久久久| www日本在线高清视频| 亚洲美女黄片视频| 日韩av在线大香蕉| 欧美大码av| 变态另类成人亚洲欧美熟女 | 国产精品 欧美亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 大型av网站在线播放| 老司机午夜福利在线观看视频| 咕卡用的链子| 国产又色又爽无遮挡免费看| 午夜福利在线观看吧| 多毛熟女@视频| 曰老女人黄片| 丰满的人妻完整版| 青草久久国产| 窝窝影院91人妻| 亚洲成人久久性| 中文字幕人妻丝袜一区二区| 久99久视频精品免费| 国产精品电影一区二区三区| 波多野结衣高清无吗| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 精品一区二区三区av网在线观看| 12—13女人毛片做爰片一| 在线天堂中文资源库| 精品欧美国产一区二区三| 九色国产91popny在线| 在线观看www视频免费| 亚洲第一青青草原| 女警被强在线播放| 国产伦一二天堂av在线观看| 美女午夜性视频免费| 人成视频在线观看免费观看| 国产片内射在线| 又紧又爽又黄一区二区| 亚洲中文日韩欧美视频| 在线观看66精品国产| 精品国产一区二区三区四区第35| av网站免费在线观看视频| 免费不卡黄色视频| 欧美国产精品va在线观看不卡| 美女大奶头视频| 99国产精品一区二区蜜桃av| 国产伦一二天堂av在线观看| 电影成人av| 变态另类成人亚洲欧美熟女 | 亚洲avbb在线观看| 极品教师在线免费播放| 国产一区二区三区视频了| 色尼玛亚洲综合影院| 国产亚洲精品久久久久5区| 亚洲伊人色综图| 成人精品一区二区免费| 免费在线观看视频国产中文字幕亚洲| 国产av一区二区精品久久| 色精品久久人妻99蜜桃| 99精品在免费线老司机午夜| 91成人精品电影| 日本免费a在线| 精品欧美一区二区三区在线| 日本黄色视频三级网站网址| 久久久久久国产a免费观看| 色综合欧美亚洲国产小说| 亚洲一区中文字幕在线| 国产91精品成人一区二区三区| 免费看十八禁软件| 亚洲 欧美一区二区三区| www.自偷自拍.com| 日日干狠狠操夜夜爽| www.精华液| 亚洲人成网站在线播放欧美日韩| 视频区欧美日本亚洲| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品一区二区www| 亚洲第一av免费看| 欧美乱妇无乱码| 夜夜躁狠狠躁天天躁| 亚洲一区二区三区不卡视频| 午夜福利欧美成人| 久久人妻熟女aⅴ| av福利片在线| 99国产精品一区二区三区| 精品一品国产午夜福利视频| 欧美乱妇无乱码| 男男h啪啪无遮挡| 女人高潮潮喷娇喘18禁视频| 啪啪无遮挡十八禁网站| 脱女人内裤的视频| 97人妻精品一区二区三区麻豆 | 久9热在线精品视频| 91在线观看av| 国产精品,欧美在线| 欧美日韩福利视频一区二区| 欧美日韩精品网址| 午夜福利免费观看在线| 久久国产精品人妻蜜桃| 中文字幕人成人乱码亚洲影| 三级毛片av免费| 老司机在亚洲福利影院| 国产精品自产拍在线观看55亚洲| 久久久久久久久免费视频了| 丁香欧美五月| 女警被强在线播放| 亚洲三区欧美一区| а√天堂www在线а√下载| 亚洲情色 制服丝袜| 久久欧美精品欧美久久欧美| 十八禁人妻一区二区| 一级毛片高清免费大全| 亚洲第一青青草原| 熟女少妇亚洲综合色aaa.| 亚洲成人久久性| 最近最新中文字幕大全免费视频| 看黄色毛片网站| 久久精品国产亚洲av香蕉五月| 99在线视频只有这里精品首页| 99国产精品一区二区三区| 男女床上黄色一级片免费看| 久久久久国产一级毛片高清牌| 999久久久精品免费观看国产| 国产精品自产拍在线观看55亚洲| 国内精品久久久久精免费| 国产成人系列免费观看| 一进一出抽搐动态| 亚洲熟女毛片儿| 久久久精品欧美日韩精品| 日本三级黄在线观看| 男人舔女人的私密视频| 桃红色精品国产亚洲av| 少妇熟女aⅴ在线视频| 岛国视频午夜一区免费看| 桃红色精品国产亚洲av| 久久久国产成人精品二区| 亚洲视频免费观看视频| 在线视频色国产色| 美女大奶头视频| 久久精品aⅴ一区二区三区四区| 欧美日本亚洲视频在线播放| 国产又色又爽无遮挡免费看| 村上凉子中文字幕在线| 精品乱码久久久久久99久播| 日韩欧美免费精品| 欧美日韩福利视频一区二区| 午夜成年电影在线免费观看| 热re99久久国产66热| 一级黄色大片毛片| 欧美日本亚洲视频在线播放| 国产私拍福利视频在线观看| 性色av乱码一区二区三区2| 99在线视频只有这里精品首页| 成人亚洲精品av一区二区| 国产熟女xx| 1024香蕉在线观看| 久久久久国产一级毛片高清牌| 不卡一级毛片| 麻豆av在线久日| 日日摸夜夜添夜夜添小说| 日韩 欧美 亚洲 中文字幕| 国产熟女午夜一区二区三区| 久久久久久免费高清国产稀缺| 亚洲七黄色美女视频| 一级毛片高清免费大全| 精品乱码久久久久久99久播| 日韩欧美在线二视频| 99国产极品粉嫩在线观看| 国产成人影院久久av| 91老司机精品| 亚洲国产毛片av蜜桃av| av片东京热男人的天堂| 此物有八面人人有两片| 久99久视频精品免费| 国产视频一区二区在线看| 免费观看人在逋| 精品国内亚洲2022精品成人| av在线播放免费不卡| 久99久视频精品免费| 亚洲欧美激情在线| 亚洲精品久久国产高清桃花| 九色亚洲精品在线播放| 99国产精品一区二区蜜桃av| 色婷婷久久久亚洲欧美| 亚洲欧美一区二区三区黑人| 精品久久久精品久久久| 在线天堂中文资源库| 中文字幕最新亚洲高清| 久久精品aⅴ一区二区三区四区| 夜夜夜夜夜久久久久| 亚洲成人久久性| www.www免费av| 久久精品亚洲熟妇少妇任你| 国产精品av久久久久免费| 精品国产乱子伦一区二区三区| 亚洲精品久久成人aⅴ小说| 日本撒尿小便嘘嘘汇集6| 不卡一级毛片| 制服诱惑二区| 国产精品亚洲一级av第二区| 亚洲国产欧美一区二区综合| 精品福利观看| 69av精品久久久久久| 精品国产国语对白av| 亚洲熟妇熟女久久| 亚洲色图综合在线观看| 天堂√8在线中文| 最近最新中文字幕大全电影3 | 国产男靠女视频免费网站| 别揉我奶头~嗯~啊~动态视频| 亚洲熟女毛片儿| 国产精品 国内视频| 欧美乱色亚洲激情| 国产av精品麻豆| 亚洲 欧美一区二区三区| 欧美黄色片欧美黄色片| 久久国产精品人妻蜜桃| 亚洲男人天堂网一区| 高清在线国产一区| 亚洲激情在线av| 国产亚洲欧美精品永久| 国产又爽黄色视频| 最新在线观看一区二区三区| 免费看a级黄色片| 久99久视频精品免费| 色哟哟哟哟哟哟| 国产成人影院久久av| 国产精品98久久久久久宅男小说| 亚洲五月天丁香| 欧美不卡视频在线免费观看 | 一个人免费在线观看的高清视频| 国产精品野战在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产精品影院| av超薄肉色丝袜交足视频| 色婷婷久久久亚洲欧美| 亚洲国产精品成人综合色| 亚洲一区二区三区不卡视频| 免费女性裸体啪啪无遮挡网站| 在线国产一区二区在线| 天天添夜夜摸| 国产成+人综合+亚洲专区| 久久久精品欧美日韩精品| 成人精品一区二区免费| 久久久精品欧美日韩精品| 看黄色毛片网站| 国产亚洲欧美在线一区二区| 看黄色毛片网站| 色综合站精品国产| 亚洲中文日韩欧美视频| 亚洲狠狠婷婷综合久久图片| 香蕉久久夜色| 18禁美女被吸乳视频| 老司机深夜福利视频在线观看| 丝袜美腿诱惑在线| 中文字幕精品免费在线观看视频| 日日干狠狠操夜夜爽| 国产成+人综合+亚洲专区| 亚洲九九香蕉| 欧美精品啪啪一区二区三区| 久久久久久久精品吃奶| 国产乱人伦免费视频| 一个人观看的视频www高清免费观看 | 精品人妻1区二区| 非洲黑人性xxxx精品又粗又长| 亚洲 国产 在线| 人人妻人人澡人人看| 国产精品综合久久久久久久免费 | 黄色成人免费大全| 黄片大片在线免费观看| 亚洲精品在线美女| 亚洲色图 男人天堂 中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 国产精品野战在线观看| 久久婷婷人人爽人人干人人爱 | 色哟哟哟哟哟哟| 欧美在线一区亚洲| 一二三四在线观看免费中文在| 国产成人精品无人区| 欧美 亚洲 国产 日韩一| 久久人妻福利社区极品人妻图片| 欧美精品啪啪一区二区三区| 最近最新免费中文字幕在线| 亚洲欧美激情在线| 国产熟女午夜一区二区三区| 久久欧美精品欧美久久欧美| www国产在线视频色| 国产又爽黄色视频| 亚洲第一电影网av| 国产精品免费一区二区三区在线| 成人18禁高潮啪啪吃奶动态图| 亚洲激情在线av| 最新美女视频免费是黄的| 91麻豆精品激情在线观看国产| 成人18禁在线播放| 日本精品一区二区三区蜜桃| 啦啦啦韩国在线观看视频| 久久久水蜜桃国产精品网| 免费一级毛片在线播放高清视频 | 国产乱人伦免费视频| 国产精品久久视频播放| 亚洲国产高清在线一区二区三 | 女人精品久久久久毛片| 日韩 欧美 亚洲 中文字幕| www.精华液| 色综合亚洲欧美另类图片| 久久国产精品人妻蜜桃| 国产一区二区激情短视频| 亚洲天堂国产精品一区在线| 精品久久久久久,| 欧美国产精品va在线观看不卡| 婷婷六月久久综合丁香| 男人操女人黄网站| 老司机深夜福利视频在线观看| 人人妻人人爽人人添夜夜欢视频| 99精品久久久久人妻精品| 国产成+人综合+亚洲专区| 变态另类丝袜制服| 高清在线国产一区| 在线视频色国产色| 国产欧美日韩综合在线一区二区| 午夜激情av网站| 叶爱在线成人免费视频播放| 久久久久久久久免费视频了| 亚洲国产看品久久| 国产麻豆69| 国产成人欧美| 亚洲五月天丁香| 免费在线观看视频国产中文字幕亚洲| 女生性感内裤真人,穿戴方法视频| 色尼玛亚洲综合影院| 亚洲国产毛片av蜜桃av| 欧美色视频一区免费| 12—13女人毛片做爰片一| 精品无人区乱码1区二区| 亚洲国产毛片av蜜桃av| 国产精品久久电影中文字幕| 国产精品二区激情视频| 黑人操中国人逼视频| 免费观看人在逋| 午夜激情av网站| 禁无遮挡网站| 老司机靠b影院| 亚洲精品美女久久久久99蜜臀| 中文字幕av电影在线播放| 90打野战视频偷拍视频| 在线十欧美十亚洲十日本专区| 午夜福利在线观看吧| 亚洲一区二区三区色噜噜| 午夜免费激情av| 黄色 视频免费看| 99久久国产精品久久久| 亚洲国产精品sss在线观看| 一级片免费观看大全| 国语自产精品视频在线第100页| 国产成人影院久久av| 成人av一区二区三区在线看| 女人被狂操c到高潮| 午夜a级毛片| 欧美久久黑人一区二区| 国产黄a三级三级三级人| 好看av亚洲va欧美ⅴa在| 亚洲国产精品成人综合色| 女同久久另类99精品国产91| 日韩有码中文字幕| 无遮挡黄片免费观看| 色尼玛亚洲综合影院| 国产欧美日韩精品亚洲av| 久久精品国产综合久久久| 亚洲国产精品成人综合色| 动漫黄色视频在线观看| 免费不卡黄色视频| 国产蜜桃级精品一区二区三区| 免费在线观看日本一区| 国产伦人伦偷精品视频| 欧美精品亚洲一区二区| 别揉我奶头~嗯~啊~动态视频| 91成人精品电影| www.熟女人妻精品国产| 女性生殖器流出的白浆| 国产精品国产高清国产av| 国产成人影院久久av| 亚洲人成网站在线播放欧美日韩| 国产极品粉嫩免费观看在线| 亚洲国产看品久久| 制服人妻中文乱码| 天天添夜夜摸| 亚洲精品一卡2卡三卡4卡5卡| 男人的好看免费观看在线视频 | 亚洲国产欧美网| 高清黄色对白视频在线免费看| 九色亚洲精品在线播放| 午夜精品久久久久久毛片777| 精品人妻1区二区| 国产激情久久老熟女| 欧美老熟妇乱子伦牲交| 久久精品国产综合久久久| 日本 欧美在线| 欧美日韩黄片免| 99久久精品国产亚洲精品| 男女做爰动态图高潮gif福利片 | 免费在线观看完整版高清| 波多野结衣巨乳人妻| 亚洲aⅴ乱码一区二区在线播放 | 黄网站色视频无遮挡免费观看| 国产精品美女特级片免费视频播放器 | 国产精华一区二区三区| 国产精品影院久久| 亚洲精品国产一区二区精华液| 日本精品一区二区三区蜜桃| 女人爽到高潮嗷嗷叫在线视频| 亚洲精华国产精华精| 日韩免费av在线播放| 一卡2卡三卡四卡精品乱码亚洲| www.999成人在线观看| 韩国av一区二区三区四区| 欧美成人一区二区免费高清观看 | 亚洲男人的天堂狠狠| 久久久水蜜桃国产精品网| 天天一区二区日本电影三级 | 97碰自拍视频| 国产精品永久免费网站| 露出奶头的视频| 免费无遮挡裸体视频| 99精品久久久久人妻精品| 中文字幕人成人乱码亚洲影| 啦啦啦韩国在线观看视频| 欧美成人性av电影在线观看| 日韩一卡2卡3卡4卡2021年| 午夜影院日韩av| 亚洲av美国av| 欧美国产日韩亚洲一区| 岛国视频午夜一区免费看| 人妻丰满熟妇av一区二区三区| www.999成人在线观看| 亚洲在线自拍视频| 99国产精品免费福利视频| 91老司机精品| 免费高清视频大片| 怎么达到女性高潮| 女生性感内裤真人,穿戴方法视频| 日日爽夜夜爽网站| 国产精品爽爽va在线观看网站 | 无人区码免费观看不卡| 亚洲少妇的诱惑av| 国产97色在线日韩免费| 久久国产乱子伦精品免费另类| 欧美一级毛片孕妇| 精品人妻在线不人妻| 国产av一区在线观看免费| 欧美老熟妇乱子伦牲交| 精品国产乱子伦一区二区三区| 国产一区二区激情短视频| 99久久99久久久精品蜜桃| 精品国产乱码久久久久久男人| 欧美精品啪啪一区二区三区| 电影成人av| 中文亚洲av片在线观看爽| 国产亚洲av嫩草精品影院| 国产视频一区二区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产精品男人的天堂亚洲| 亚洲 欧美一区二区三区| 免费搜索国产男女视频| 亚洲av成人不卡在线观看播放网| 色尼玛亚洲综合影院| 大陆偷拍与自拍| 一边摸一边做爽爽视频免费| 首页视频小说图片口味搜索| 午夜福利,免费看| 不卡av一区二区三区| 波多野结衣一区麻豆| 男人舔女人下体高潮全视频| 91老司机精品| 免费人成视频x8x8入口观看| 神马国产精品三级电影在线观看 | 色综合欧美亚洲国产小说| 国产一区二区激情短视频| 久久国产精品影院| 日本vs欧美在线观看视频| 国产精品,欧美在线| 99国产精品一区二区三区| 精品熟女少妇八av免费久了| 午夜福利一区二区在线看| 91九色精品人成在线观看| 大香蕉久久成人网| 日韩成人在线观看一区二区三区| 日日夜夜操网爽| 18美女黄网站色大片免费观看| 91成人精品电影| 啦啦啦韩国在线观看视频| 亚洲五月天丁香| 亚洲自拍偷在线| 成人欧美大片| 久久天躁狠狠躁夜夜2o2o| 最好的美女福利视频网| 国产精品一区二区免费欧美| 日韩av在线大香蕉| 人人妻,人人澡人人爽秒播| 人人澡人人妻人| 久久久久久久精品吃奶| av天堂在线播放| 怎么达到女性高潮| 久久天躁狠狠躁夜夜2o2o| 欧美av亚洲av综合av国产av| 看免费av毛片| 亚洲国产欧美一区二区综合| 村上凉子中文字幕在线| 日日爽夜夜爽网站| 十八禁人妻一区二区| 国产又色又爽无遮挡免费看| 99久久99久久久精品蜜桃| 好男人在线观看高清免费视频 | 国产成人影院久久av| 日本欧美视频一区| 亚洲国产欧美网| 久9热在线精品视频| 久久天堂一区二区三区四区|