• <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)的建模與仿真
    亚洲人成网站在线播放欧美日韩| av视频在线观看入口| 人妻制服诱惑在线中文字幕| 欧美一级a爱片免费观看看| 亚洲不卡免费看| 男女那种视频在线观看| 一级黄片播放器| 男人狂女人下面高潮的视频| 18美女黄网站色大片免费观看| 日韩av在线大香蕉| 亚洲成人中文字幕在线播放| 亚洲最大成人中文| 色综合婷婷激情| 直男gayav资源| 18美女黄网站色大片免费观看| 日韩欧美精品免费久久 | 日韩欧美在线乱码| netflix在线观看网站| 天天一区二区日本电影三级| 黄色日韩在线| 91久久精品电影网| 精品人妻1区二区| 最新中文字幕久久久久| 久久中文看片网| 成人特级av手机在线观看| 亚洲国产欧洲综合997久久,| 国产爱豆传媒在线观看| 亚洲无线在线观看| 国产精品亚洲av一区麻豆| 黄色丝袜av网址大全| 亚洲人成电影免费在线| 在现免费观看毛片| 成人av在线播放网站| 日韩欧美三级三区| 日韩欧美一区二区三区在线观看| 一二三四社区在线视频社区8| 深夜精品福利| 国产一区二区亚洲精品在线观看| 极品教师在线免费播放| 欧美成人性av电影在线观看| 丝袜美腿在线中文| 69av精品久久久久久| 性欧美人与动物交配| 18+在线观看网站| 在线观看av片永久免费下载| 村上凉子中文字幕在线| 看免费av毛片| 国产国拍精品亚洲av在线观看| 少妇人妻一区二区三区视频| www日本黄色视频网| 国产精品久久久久久精品电影| 精品乱码久久久久久99久播| 色综合站精品国产| 少妇裸体淫交视频免费看高清| 露出奶头的视频| 麻豆成人av在线观看| 国产乱人伦免费视频| 永久网站在线| 久久久久免费精品人妻一区二区| 亚洲,欧美,日韩| 看十八女毛片水多多多| 久久精品91蜜桃| 国产乱人视频| 国产一区二区在线观看日韩| 宅男免费午夜| 午夜福利免费观看在线| 免费大片18禁| 一边摸一边抽搐一进一小说| 精品人妻熟女av久视频| 亚洲成a人片在线一区二区| 亚洲久久久久久中文字幕| 日本一本二区三区精品| 亚洲中文字幕一区二区三区有码在线看| 99热只有精品国产| 国产精品美女特级片免费视频播放器| 看片在线看免费视频| 1000部很黄的大片| 黄色视频,在线免费观看| 熟女人妻精品中文字幕| 中文资源天堂在线| 18禁黄网站禁片免费观看直播| 国产伦人伦偷精品视频| 成人鲁丝片一二三区免费| 国产精品一区二区免费欧美| 欧美最黄视频在线播放免费| 亚洲黑人精品在线| 精品人妻视频免费看| 露出奶头的视频| 久久精品综合一区二区三区| 欧美黄色片欧美黄色片| 18禁在线播放成人免费| 别揉我奶头~嗯~啊~动态视频| 中文字幕免费在线视频6| 日本撒尿小便嘘嘘汇集6| 亚洲第一欧美日韩一区二区三区| 99热这里只有是精品50| 国产在线精品亚洲第一网站| 91九色精品人成在线观看| 女同久久另类99精品国产91| 日本与韩国留学比较| 亚洲自偷自拍三级| 毛片一级片免费看久久久久 | 久久国产乱子伦精品免费另类| 亚洲片人在线观看| 国产伦人伦偷精品视频| 18美女黄网站色大片免费观看| 精品免费久久久久久久清纯| 97热精品久久久久久| bbb黄色大片| 色综合欧美亚洲国产小说| 国产aⅴ精品一区二区三区波| 99热精品在线国产| 国产精品久久久久久久电影| 在线观看66精品国产| 亚洲中文字幕日韩| 小说图片视频综合网站| 亚洲av五月六月丁香网| 黄色配什么色好看| 亚洲成人中文字幕在线播放| 麻豆国产97在线/欧美| 在线观看午夜福利视频| 校园春色视频在线观看| 亚洲无线观看免费| 国产精品久久久久久久电影| 日韩中字成人| 国产黄a三级三级三级人| 欧美成人性av电影在线观看| 亚洲av.av天堂| 亚洲欧美精品综合久久99| 1024手机看黄色片| 亚洲欧美激情综合另类| 欧美丝袜亚洲另类 | 国产单亲对白刺激| 国产 一区 欧美 日韩| 免费在线观看日本一区| 特级一级黄色大片| 国产一级毛片七仙女欲春2| 少妇人妻一区二区三区视频| 午夜福利欧美成人| 亚洲精品成人久久久久久| 男女下面进入的视频免费午夜| 国产三级在线视频| 可以在线观看毛片的网站| 免费大片18禁| 亚洲欧美日韩卡通动漫| 少妇人妻精品综合一区二区 | 久久精品91蜜桃| 亚洲内射少妇av| 国产乱人伦免费视频| 国内少妇人妻偷人精品xxx网站| 亚洲成人中文字幕在线播放| 深夜精品福利| 91在线精品国自产拍蜜月| 亚洲国产色片| 91麻豆av在线| av专区在线播放| 亚洲精品在线美女| 色噜噜av男人的天堂激情| 我要看日韩黄色一级片| 国产精品久久久久久人妻精品电影| 亚洲欧美激情综合另类| 窝窝影院91人妻| 色哟哟哟哟哟哟| 91麻豆精品激情在线观看国产| 91久久精品国产一区二区成人| 国产亚洲av嫩草精品影院| 如何舔出高潮| 久久久久久久久久成人| 最好的美女福利视频网| 亚洲欧美清纯卡通| 少妇熟女aⅴ在线视频| 亚洲av美国av| 婷婷丁香在线五月| ponron亚洲| 搡老妇女老女人老熟妇| 成人一区二区视频在线观看| 亚洲美女黄片视频| 成年人黄色毛片网站| 国产精品不卡视频一区二区 | 精品人妻一区二区三区麻豆 | 国产色爽女视频免费观看| 91午夜精品亚洲一区二区三区 | 一卡2卡三卡四卡精品乱码亚洲| 精品不卡国产一区二区三区| 制服丝袜大香蕉在线| 精品99又大又爽又粗少妇毛片 | 国产av不卡久久| 欧美黄色淫秽网站| 国产高潮美女av| 69人妻影院| 亚洲在线观看片| 欧美又色又爽又黄视频| 日本一二三区视频观看| 亚洲国产精品合色在线| 男女之事视频高清在线观看| 乱码一卡2卡4卡精品| 久久精品人妻少妇| 18美女黄网站色大片免费观看| 看片在线看免费视频| 久久精品久久久久久噜噜老黄 | 精品久久久久久久久久免费视频| 久久久久久久久久成人| 美女黄网站色视频| 精品午夜福利视频在线观看一区| 欧美一区二区精品小视频在线| 男女之事视频高清在线观看| 久久久久久久亚洲中文字幕 | 欧美国产日韩亚洲一区| 亚洲美女黄片视频| 精品熟女少妇八av免费久了| 久久香蕉精品热| a级毛片a级免费在线| 日本在线视频免费播放| 日韩欧美精品v在线| 青草久久国产| 婷婷精品国产亚洲av| 久久久久久国产a免费观看| 精品99又大又爽又粗少妇毛片 | 女生性感内裤真人,穿戴方法视频| 欧美中文日本在线观看视频| 老熟妇仑乱视频hdxx| 国产精品女同一区二区软件 | 热99在线观看视频| 国产精品国产高清国产av| 亚洲 欧美 日韩 在线 免费| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 极品教师在线视频| 日韩国内少妇激情av| 无遮挡黄片免费观看| 观看美女的网站| 久久精品国产99精品国产亚洲性色| 国产亚洲精品av在线| 日本与韩国留学比较| 免费大片18禁| 亚洲精品456在线播放app | 丁香六月欧美| 男人舔奶头视频| 一个人免费在线观看电影| 中出人妻视频一区二区| 国内少妇人妻偷人精品xxx网站| 久久亚洲真实| 国产精品乱码一区二三区的特点| 亚洲最大成人av| 九色国产91popny在线| 99国产极品粉嫩在线观看| 亚洲五月婷婷丁香| 三级男女做爰猛烈吃奶摸视频| 亚洲av电影在线进入| 国产大屁股一区二区在线视频| 国产精品美女特级片免费视频播放器| 小说图片视频综合网站| 免费搜索国产男女视频| 免费看光身美女| 女生性感内裤真人,穿戴方法视频| 亚洲第一电影网av| 精品久久久久久久久亚洲 | 看黄色毛片网站| 色综合亚洲欧美另类图片| 97超级碰碰碰精品色视频在线观看| 精品国产亚洲在线| 亚洲五月天丁香| 国产亚洲精品综合一区在线观看| 免费高清视频大片| 久久伊人香网站| 欧美日本亚洲视频在线播放| 亚洲成人久久爱视频| 国产毛片a区久久久久| 精品久久国产蜜桃| 精品一区二区三区av网在线观看| 丰满的人妻完整版| 色av中文字幕| 丁香欧美五月| 午夜激情福利司机影院| 我的老师免费观看完整版| 欧美黑人巨大hd| 搡老熟女国产l中国老女人| 九色成人免费人妻av| 久久久久精品国产欧美久久久| 精品久久久久久,| 国产人妻一区二区三区在| 国产精品亚洲一级av第二区| 18美女黄网站色大片免费观看| netflix在线观看网站| 热99在线观看视频| 久久午夜福利片| av在线天堂中文字幕| 国产精品一区二区三区四区久久| 嫁个100分男人电影在线观看| 又爽又黄a免费视频| 久久精品人妻少妇| 亚洲av第一区精品v没综合| 99久久成人亚洲精品观看| 成人毛片a级毛片在线播放| 黄色视频,在线免费观看| 级片在线观看| 神马国产精品三级电影在线观看| 午夜福利在线观看吧| 午夜日韩欧美国产| 日韩成人在线观看一区二区三区| 免费搜索国产男女视频| 人人妻人人看人人澡| 搡女人真爽免费视频火全软件 | 老司机深夜福利视频在线观看| 我的女老师完整版在线观看| 给我免费播放毛片高清在线观看| 日日夜夜操网爽| 色噜噜av男人的天堂激情| 欧美xxxx性猛交bbbb| 国产精品一区二区三区四区免费观看 | 黄色配什么色好看| 国产主播在线观看一区二区| 久久午夜福利片| 国产精品电影一区二区三区| 村上凉子中文字幕在线| 久久天躁狠狠躁夜夜2o2o| 成人欧美大片| 国产真实伦视频高清在线观看 | 久久九九热精品免费| 亚洲av中文字字幕乱码综合| 亚洲精品日韩av片在线观看| 在线观看午夜福利视频| 午夜视频国产福利| 国产不卡一卡二| 欧美乱妇无乱码| 国产午夜精品论理片| 国产精品免费一区二区三区在线| 欧美日韩亚洲国产一区二区在线观看| 五月伊人婷婷丁香| 色综合欧美亚洲国产小说| 久久国产乱子伦精品免费另类| aaaaa片日本免费| 97碰自拍视频| 99久久精品国产亚洲精品| 久久精品人妻少妇| 久久久久性生活片| 日本五十路高清| 三级毛片av免费| 日韩免费av在线播放| 午夜精品久久久久久毛片777| 亚洲 欧美 日韩 在线 免费| 在现免费观看毛片| 国产精品综合久久久久久久免费| 久99久视频精品免费| 每晚都被弄得嗷嗷叫到高潮| 可以在线观看的亚洲视频| 我的老师免费观看完整版| 亚洲av免费高清在线观看| 亚洲电影在线观看av| 十八禁人妻一区二区| 久久国产精品影院| 最后的刺客免费高清国语| 国产精品98久久久久久宅男小说| 久久午夜福利片| 波野结衣二区三区在线| 久久久久九九精品影院| 国产极品精品免费视频能看的| www.www免费av| 在线观看免费视频日本深夜| 久久性视频一级片| av在线蜜桃| 精品一区二区三区视频在线| 成人国产一区最新在线观看| 久久99热6这里只有精品| av中文乱码字幕在线| aaaaa片日本免费| 亚洲真实伦在线观看| 亚洲成a人片在线一区二区| av在线观看视频网站免费| 桃色一区二区三区在线观看| 成人三级黄色视频| 色av中文字幕| av天堂在线播放| 男人的好看免费观看在线视频| www.熟女人妻精品国产| 极品教师在线视频| 看十八女毛片水多多多| 哪里可以看免费的av片| 亚洲欧美精品综合久久99| 欧美日韩黄片免| 日本免费a在线| 高潮久久久久久久久久久不卡| www日本黄色视频网| 能在线免费观看的黄片| 婷婷精品国产亚洲av| 我要看日韩黄色一级片| 久久久久久久亚洲中文字幕 | 国产精品嫩草影院av在线观看 | 国产色爽女视频免费观看| 深夜a级毛片| 国产美女午夜福利| 亚洲人成网站在线播放欧美日韩| 国产伦人伦偷精品视频| 日日夜夜操网爽| 国产伦精品一区二区三区视频9| 国内精品久久久久久久电影| 精品一区二区三区视频在线观看免费| 欧美日韩乱码在线| 欧美精品啪啪一区二区三区| 超碰av人人做人人爽久久| 精品一区二区三区av网在线观看| 亚洲国产欧洲综合997久久,| 天堂av国产一区二区熟女人妻| 美女免费视频网站| xxxwww97欧美| 男人舔女人下体高潮全视频| 日韩欧美在线二视频| 99国产极品粉嫩在线观看| 欧美日韩福利视频一区二区| 欧美日本亚洲视频在线播放| 中文字幕免费在线视频6| 国产蜜桃级精品一区二区三区| 国产大屁股一区二区在线视频| 深夜精品福利| 看黄色毛片网站| 国产精品女同一区二区软件 | 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产亚洲av天美| 丰满的人妻完整版| 免费在线观看日本一区| 99国产综合亚洲精品| 国产免费男女视频| 欧美xxxx黑人xx丫x性爽| 成人特级黄色片久久久久久久| 日韩亚洲欧美综合| av中文乱码字幕在线| 国产精品久久久久久久久免 | 欧美中文日本在线观看视频| 男女做爰动态图高潮gif福利片| 男人狂女人下面高潮的视频| 91字幕亚洲| 在线观看免费视频日本深夜| 亚洲,欧美精品.| 熟女人妻精品中文字幕| 身体一侧抽搐| 亚洲性夜色夜夜综合| 一个人看视频在线观看www免费| 69人妻影院| .国产精品久久| 亚洲成av人片免费观看| 国产三级中文精品| 国产单亲对白刺激| 久久久久久久久大av| 久久久精品大字幕| 九九热线精品视视频播放| 久久久久久久久久成人| 午夜福利在线观看吧| 欧美日韩黄片免| 非洲黑人性xxxx精品又粗又长| 精品人妻视频免费看| 国产久久久一区二区三区| 久久久久久久久久黄片| 99久久99久久久精品蜜桃| 欧美成人性av电影在线观看| 国产精品久久久久久久电影| 国产麻豆成人av免费视频| 精品人妻视频免费看| 丰满人妻熟妇乱又伦精品不卡| 欧美精品国产亚洲| 国产精品亚洲美女久久久| 免费搜索国产男女视频| 我的老师免费观看完整版| 日本撒尿小便嘘嘘汇集6| 在线观看一区二区三区| 淫妇啪啪啪对白视频| 国产精品爽爽va在线观看网站| 亚洲欧美日韩东京热| 国产成人影院久久av| 久久热精品热| 婷婷丁香在线五月| 一进一出好大好爽视频| 99久久精品一区二区三区| 欧美性感艳星| 国产乱人伦免费视频| 男插女下体视频免费在线播放| 午夜福利高清视频| 老司机午夜福利在线观看视频| 久久久成人免费电影| 怎么达到女性高潮| 日韩欧美免费精品| 亚洲成人免费电影在线观看| 欧美xxxx性猛交bbbb| 久久久久亚洲av毛片大全| 99久久精品国产亚洲精品| 亚洲精品一区av在线观看| 很黄的视频免费| 免费看a级黄色片| 国产单亲对白刺激| 欧美在线黄色| 蜜桃久久精品国产亚洲av| 午夜福利欧美成人| 亚洲美女搞黄在线观看 | 蜜桃久久精品国产亚洲av| 久久国产乱子免费精品| 成人鲁丝片一二三区免费| 国产男靠女视频免费网站| 精品无人区乱码1区二区| 此物有八面人人有两片| 亚洲人与动物交配视频| 大型黄色视频在线免费观看| 亚洲人成电影免费在线| 午夜福利高清视频| 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 一进一出抽搐动态| 国产三级中文精品| 国内精品久久久久久久电影| 神马国产精品三级电影在线观看| 国内精品久久久久久久电影| 婷婷精品国产亚洲av在线| 51国产日韩欧美| 91午夜精品亚洲一区二区三区 | 亚洲欧美日韩高清在线视频| 久久久精品欧美日韩精品| .国产精品久久| 国产麻豆成人av免费视频| 在线观看一区二区三区| 成人欧美大片| 免费人成视频x8x8入口观看| ponron亚洲| 校园春色视频在线观看| 一二三四社区在线视频社区8| 搞女人的毛片| 欧美乱妇无乱码| 无遮挡黄片免费观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美高清性xxxxhd video| 变态另类丝袜制服| 757午夜福利合集在线观看| 亚洲欧美日韩高清专用| eeuss影院久久| 日韩大尺度精品在线看网址| а√天堂www在线а√下载| 久久精品91蜜桃| 成人鲁丝片一二三区免费| 亚洲av免费在线观看| 一本综合久久免费| 久久九九热精品免费| 少妇人妻一区二区三区视频| 国产色婷婷99| 欧美在线一区亚洲| 国产乱人视频| 久久久久精品国产欧美久久久| 中文字幕熟女人妻在线| 麻豆成人av在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品久久久久久精品电影| av福利片在线观看| 禁无遮挡网站| 日韩欧美 国产精品| 久久久久久久久大av| 亚洲五月婷婷丁香| 欧美区成人在线视频| 夜夜爽天天搞| 欧美区成人在线视频| 十八禁人妻一区二区| 在线a可以看的网站| 一本久久中文字幕| 悠悠久久av| 非洲黑人性xxxx精品又粗又长| 亚洲美女黄片视频| 亚洲人成网站在线播| 欧美午夜高清在线| 天堂网av新在线| 色综合站精品国产| 搞女人的毛片| 在线观看午夜福利视频| 国产午夜福利久久久久久| 精品熟女少妇八av免费久了| av在线观看视频网站免费| 极品教师在线视频| 亚洲综合色惰| 国产探花在线观看一区二区| 我要看日韩黄色一级片| 亚洲国产精品合色在线| 99riav亚洲国产免费| 国产午夜福利久久久久久| 18禁黄网站禁片午夜丰满| 日韩成人在线观看一区二区三区| 亚洲不卡免费看| 一二三四社区在线视频社区8| 淫妇啪啪啪对白视频| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| 直男gayav资源| 日韩欧美免费精品| 一区二区三区高清视频在线| 午夜福利成人在线免费观看| 欧美日韩中文字幕国产精品一区二区三区| .国产精品久久| 成人鲁丝片一二三区免费| 哪里可以看免费的av片| 亚洲专区国产一区二区| 最新在线观看一区二区三区| 少妇的逼好多水| 三级男女做爰猛烈吃奶摸视频| 欧美3d第一页| 亚洲内射少妇av| 欧美日韩综合久久久久久 | 精品国产三级普通话版| 欧美日韩综合久久久久久 | 国产一区二区三区视频了| 久久久久久久精品吃奶| 国产午夜福利久久久久久| 免费在线观看亚洲国产| 美女大奶头视频| 国产大屁股一区二区在线视频| 久久精品综合一区二区三区| 日本熟妇午夜| 男人的好看免费观看在线视频| 国产乱人视频| 欧美另类亚洲清纯唯美| 国产高清视频在线观看网站| 日日摸夜夜添夜夜添av毛片 | 中文字幕熟女人妻在线|