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

    基于整群抽樣和支持向量回歸模型的高功率半導(dǎo)體激光器剩余使用壽命預(yù)測

    2021-07-15 02:55:04嚴(yán)建文鐘小虎郭三敏1
    中國機(jī)械工程 2021年13期
    關(guān)鍵詞:整群激光器階段

    嚴(yán)建文 鐘小虎 范 煜 郭三敏1,

    1.合肥工業(yè)大學(xué)管理學(xué)院,合肥,2300092.合肥工業(yè)大學(xué)航空結(jié)構(gòu)件成形制造與裝備安徽省重點(diǎn)實(shí)驗(yàn)室,合肥,2300093.浙江大學(xué)機(jī)械工程學(xué)院,杭州,3100584.安徽皖維集團(tuán)有限責(zé)任公司,合肥,238002

    0 引言

    高功率半導(dǎo)體激光器及其陣列具有體積小、質(zhì)量小、能耗低、光斑易調(diào)節(jié)、光電轉(zhuǎn)換效率較高的優(yōu)點(diǎn),已廣泛應(yīng)用于金屬材料加工,該種激光器可用于連續(xù)性焊接不同型號的合金鋼,獲得大面積深度均勻的相變硬化層,也能夠精確地控制熔覆層結(jié)構(gòu)及其幾何形狀[1]。為保證材料加工過程的安全性和可靠性,需利用預(yù)測與健康管理(prognostics and health management, PHM)技術(shù)對高功率半導(dǎo)體激光器的狀態(tài)進(jìn)行監(jiān)測和分析[2]。PHM技術(shù)的主要目的在于預(yù)測故障發(fā)生的時(shí)間和地點(diǎn),估計(jì)剩余使用壽命(remaining useful life, RUL)的分布或期望[3],提高激光器的可靠性。準(zhǔn)確的RUL預(yù)測可為決策者提前制定維修計(jì)劃和優(yōu)化供應(yīng)鏈管理,減少不必要的維修或更換成本,提供有效的信息。

    相比于傳統(tǒng)的基于可靠性的方法[4-5],現(xiàn)有的數(shù)據(jù)驅(qū)動(dòng)方法主要通過對狀態(tài)監(jiān)測或檢測數(shù)據(jù)進(jìn)行分析,捕捉數(shù)據(jù)中潛在的退化規(guī)律和故障失效信息,通??梢苑譃閮深?。①基于退化模型的方法[6](如維納過程模型[7]、伽馬過程模型[8]、逆高斯過程模型[9]、隱(半)馬爾可夫模型[10-11]等),旨在對設(shè)備的退化演化過程進(jìn)行建模,如采用逆高斯過程模型對退化過程進(jìn)行分析,以預(yù)估GaAs激光器的RUL分布[12-14]。這類方法均假定退化模型是事先知道的,模型參數(shù)利用狀態(tài)監(jiān)測數(shù)據(jù)進(jìn)行估計(jì),然而實(shí)際工程中很難事先確定一個(gè)適當(dāng)?shù)耐嘶P停彝嘶P瓦x擇不當(dāng)將會嚴(yán)重影響RUL的預(yù)測精度[15]。②基于機(jī)器學(xué)習(xí)的方法[16](如人工神經(jīng)網(wǎng)絡(luò)(artificial neutral network, ANN)[17]、基于相似性的學(xué)習(xí)算法(SbRP)[18]、支持向量回歸(support vector regression, SVR)方法[19-20]等),試圖將退化數(shù)據(jù)與RUL之間的關(guān)系進(jìn)行直接映射。如文獻(xiàn)[21]利用SbRP方法對GaAs激光器的RUL進(jìn)行了預(yù)測;文獻(xiàn)[22]基于自回歸模型和SbRP方法,提出了一種基于相似性的差值分析(SbDA)方法用于GaAs激光器的RUL預(yù)測。

    基于機(jī)器學(xué)習(xí)的方法能夠有效克服退化模型未知的問題,近年來已得到了廣泛的研究和應(yīng)用。然而,如何可靠地在失效歷史數(shù)據(jù)有限的情況下實(shí)現(xiàn)對目標(biāo)設(shè)備的實(shí)時(shí)預(yù)測仍然具有極大的挑戰(zhàn)。SVR方法是支持向量機(jī)(support vector machine,SVM)算法在RUL預(yù)測領(lǐng)域最常見的應(yīng)用方式,可有效地解決小樣本情況下的預(yù)測問題。BENKEDJOUH等[19-20]利用SVR方法將退化時(shí)間序列映射成非線性回歸,然后將得到的回歸擬合到功率模型中,并用于機(jī)械設(shè)備的RUL預(yù)測。LIU等[23]建立了一種改進(jìn)的概率SVR模型用于預(yù)測核電站裝備部件的RUL。FUMEO等[24]通過優(yōu)化精度與計(jì)算效率之間的權(quán)衡,開發(fā)了一種在線SVR模型用于軸承的RUL預(yù)測。上述研究雖然拓展了SVR模型在RUL預(yù)測問題中的應(yīng)用,但均未考慮關(guān)鍵預(yù)警階段特別是臨近故障失效前的預(yù)測可靠性問題,即使SVR方法所訓(xùn)練模型的回歸曲線的整體誤差最小,也有可能因?yàn)殛P(guān)鍵預(yù)警階段數(shù)據(jù)擬合得不好而導(dǎo)致實(shí)時(shí)預(yù)測的結(jié)果未必能可靠地支持維護(hù)決策。關(guān)鍵預(yù)警階段RUL預(yù)測的準(zhǔn)確性是決策者在運(yùn)行可靠性和成本之間權(quán)衡的決策基礎(chǔ),特別是臨近故障失效前。HUYNH等[25]提出了一種考慮RUL預(yù)測精度的預(yù)測維護(hù)決策框架,其決策框架的性能優(yōu)劣取決于臨近失效階段的RUL預(yù)測精度。ZHAO等[26]根據(jù)不同的可靠性水平將產(chǎn)品退化過程簡化為5種狀態(tài),其中,狀態(tài)0為正常運(yùn)行狀態(tài),狀態(tài)1為退化加速狀態(tài),狀態(tài)2和狀態(tài)3為關(guān)鍵預(yù)警狀態(tài),狀態(tài)4為失效狀態(tài)。相對于整個(gè)退化過程,狀態(tài)2和狀態(tài)3階段的RUL預(yù)測結(jié)果對維護(hù)決策更具有實(shí)際意義。

    基于上述情況,本文提出了一種基于整群抽樣和SVR模型的方法用于高功率半導(dǎo)體激光器的RUL預(yù)測。對測試樣本狀態(tài)2和狀態(tài)3階段的數(shù)據(jù)進(jìn)行多次整群抽樣后用于SVR模型測試,SVR模型中的參數(shù)使得SVR模型對訓(xùn)練樣本的后期數(shù)據(jù)擬合得更好,更能滿足維護(hù)決策需要。最后,所提方法的有效性、準(zhǔn)確性和穩(wěn)健性通過GaAs激光器數(shù)據(jù)集得到了驗(yàn)證。

    1 基于整群抽樣和SVR模型的剩余壽命預(yù)測框架

    SVR方法基于結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則,保證了模型最大的泛化能力,特別適合處理小樣本、非線性和高維數(shù)據(jù),它在RUL預(yù)測領(lǐng)域最常見的應(yīng)用方式為最小二乘支持向量回歸(least square support vector regression, LS-SVR)方法。給定一組訓(xùn)練集{(xp,yp)|p=1,2,…,n},其中xp∈RN為N維輸入向量(N≥1),yp∈R為輸出值,n為訓(xùn)練樣本總數(shù),則LS-SVR模型可用一個(gè)非線性映射函數(shù)φ(·)來表示,即

    y=wTφ(x)+b

    (1)

    式中,w為權(quán)向量;b為偏差項(xiàng)。

    通過將原始輸入數(shù)據(jù)映射到高維空間,非線性可分問題在該空間上變?yōu)榫€性可分問題,則LS-SVR問題可變?yōu)榍蠼鉃槿缦聝?yōu)化問題:

    (2)

    式中,ep為隨機(jī)誤差;γ為調(diào)節(jié)參數(shù);m為支持向量的個(gè)數(shù),且m≤n。

    此外,根據(jù)拉格朗日函數(shù)和KKT定理,采用非線性核函數(shù)時(shí)的LS-SVR可表示為

    (3)

    式中,ap為拉格朗日乘子;K(x,xp)表示非線性核函數(shù)。

    圖1給出了某個(gè)測試樣本采用現(xiàn)有SVR方法的測試結(jié)果,可以看出,即使訓(xùn)練得到了相對于整個(gè)數(shù)據(jù)集最優(yōu)的模型參數(shù),但狀態(tài)2和狀態(tài)3這兩個(gè)階段的擬合曲線與實(shí)際曲線相差甚遠(yuǎn),遠(yuǎn)遠(yuǎn)達(dá)不到制定最佳維護(hù)策略的要求,甚至?xí)斐商崆巴C(jī)或失效風(fēng)險(xiǎn)。如何提高算法預(yù)測結(jié)果的可靠性是當(dāng)前工業(yè)應(yīng)用的實(shí)際需求,然而直接對式(2)進(jìn)行改良是個(gè)極大的挑戰(zhàn),可能會使算法失去其一般性。鑒于此,本文提出了基于整群抽樣的方法,通過對狀態(tài)2和狀態(tài)3兩個(gè)階段的測試數(shù)據(jù)進(jìn)行多次整群抽樣來間接增大上述兩個(gè)階段數(shù)據(jù)在算法實(shí)現(xiàn)過程中的誤差權(quán)重,使得算法更側(cè)重于這兩個(gè)階段的預(yù)測準(zhǔn)確度。假設(shè)Z={z1,z2,…,zp,…,zq,…,zn}為數(shù)據(jù)集,其中z1,z2,…,zn為數(shù)據(jù)向量,則整群抽樣方法在本文方法中的應(yīng)用形式如圖2所示,k(k=1,2,…,K)為整個(gè)退化過程中的某一階段,其中K為整個(gè)退化過程劃分的階段數(shù)。使用整群抽樣的本文方法與傳統(tǒng)SVR方法主要區(qū)別在于:本文方法將測試樣本k階段的數(shù)據(jù)Zk={zp,…,zq}進(jìn)行整體抽樣h次并重新放回,從而得到新的測試樣本,對應(yīng)的測試數(shù)據(jù)較原始測試數(shù)據(jù)多出了h組k階段的數(shù)據(jù),

    圖1 傳統(tǒng)SVR模型某樣本測試結(jié)果Fig.1 Testing result of a sample by traditional SVR model

    圖2 整群抽樣在本文方法中的應(yīng)用形式Fig.2 Usage mode of cluster sampling in the proposed method

    用該測試數(shù)據(jù)進(jìn)行測試得到的最優(yōu)模型參數(shù)即為k階段的模型參數(shù),對應(yīng)模型對整個(gè)退化過程中的k階段也是擬合最優(yōu)的。整群抽樣后,Zk中每個(gè)數(shù)據(jù)點(diǎn)的權(quán)重從1/n增大到wZk,其中,wZk可由下式計(jì)算:

    (4)

    圖3 基于整群抽樣和SVR模型的剩余使用壽命預(yù)測流程Fig.3 RUL prediction flow chart based on cluster sampling and SVR model

    2 實(shí)例分析

    2.1 數(shù)據(jù)來源

    本文采用的GaAs激光器是一種廣泛應(yīng)用的高功率半導(dǎo)體激光器,其相關(guān)數(shù)據(jù)參見文獻(xiàn)[22]。本研究的所有記錄數(shù)據(jù)如表1所示,均為GaAs激光器在80 ℃恒定熱應(yīng)力下工作電流隨時(shí)間變化而增加的百分比,已被廣泛用作性能退化指標(biāo)進(jìn)行RUL估計(jì)[12-14,21-22]。表1中共有15個(gè)實(shí)驗(yàn)樣本,每250 h測試一次數(shù)據(jù),至4000 h為止。假設(shè)產(chǎn)品的失效閾值為6(即工作電流增加6%,GaAs激光器功能失效),狀態(tài)2和狀態(tài)3下對應(yīng)的性能退化數(shù)據(jù)分別為[4,5)和[5,6),則本實(shí)驗(yàn)注重性能退化數(shù)據(jù)在[4,6)之間的RUL預(yù)測準(zhǔn)確性,實(shí)驗(yàn)數(shù)據(jù)退化曲線見圖4。由于GaAs激光器的退化曲線均近似于線性回歸模型,因此其各自的失效時(shí)間可通過非參數(shù)局部線性估計(jì)方法來確定[27]。選取樣本1~5為測試樣本,樣本6~15為訓(xùn)練樣本,采用本文方法對測試樣本在狀態(tài)2和狀態(tài)3下的RUL進(jìn)行估計(jì),并與傳統(tǒng)SVR方法、SbRP方法和SbDA方法這三種小樣本方法進(jìn)行比較分析。性能評價(jià)準(zhǔn)則依據(jù)平均絕對誤差(mean absolute error,MAE),MAE值越小,表明性能越佳。

    圖4 GaAs激光器實(shí)驗(yàn)數(shù)據(jù)退化曲線Fig.4 Degradation curves of GaAs laser experimental data

    表1 GaAs激光器在80 ℃熱應(yīng)力下的性能退化數(shù)據(jù)(工作電流增加百分比)

    2.2 結(jié)果與分析

    本實(shí)驗(yàn)中,設(shè)定最大抽樣次數(shù)H=10,分別對測試樣本狀態(tài)2和狀態(tài)3下的數(shù)據(jù)進(jìn)行整群抽樣,則可獲得121種測試結(jié)果。狀態(tài)2和狀態(tài)3下的抽樣次數(shù)均為0時(shí),本文方法等同于傳統(tǒng)SVR方法。以樣本1為例,在狀態(tài)2和狀態(tài)3兩個(gè)狀態(tài)下所有測試結(jié)果的誤差如圖5所示,其中,最大MAE為473.76 h,最小MAE為4.27 h。測試結(jié)果表明,采用本文方法可以有效地減小離失效閾值更為接近的關(guān)鍵預(yù)警階段的預(yù)測誤差,可以使得設(shè)備在更接近失效階段預(yù)測RUL更加準(zhǔn)確,更有利于維護(hù)決策。在本實(shí)驗(yàn)中,假定性能退化狀態(tài)值取4時(shí)為關(guān)鍵預(yù)警值,當(dāng)性能退化狀態(tài)值超過4但不超過6時(shí)則認(rèn)為設(shè)備性能狀態(tài)處于關(guān)鍵預(yù)警階段,即狀態(tài)2和狀態(tài)3為關(guān)鍵預(yù)警階段。對于樣本1,在1500~2250 h時(shí)的性能退化數(shù)據(jù)超過關(guān)鍵預(yù)警值且處于狀態(tài)2和狀態(tài)3下,期間的預(yù)測結(jié)果直接關(guān)系到維護(hù)決策任務(wù)的經(jīng)濟(jì)性和可靠性。測試結(jié)果如表2所示,假設(shè)決策者要求的預(yù)測結(jié)果置信區(qū)間為預(yù)測值±10%,則在樣本1的測試結(jié)果中,僅本文方法在1500 h、1750 h和2000 h處達(dá)到要求,傳統(tǒng)SVR方法、SbRP方法和SbDA方法的預(yù)測結(jié)果均不符合要求。這主要是因?yàn)樵谟?xùn)練樣本有限的情況下,訓(xùn)練樣本的性能退化模式不能很好地概括產(chǎn)品總體的所有退化模式,當(dāng)要進(jìn)行預(yù)測的目標(biāo)樣本的退化曲線與訓(xùn)練樣本的群體退化曲線差異較大時(shí),有限訓(xùn)練樣本的不足就很明顯,這會造成傳統(tǒng)SVR方法、SbRP方法和SbDA方法在實(shí)際應(yīng)用時(shí)很容易受到數(shù)據(jù)不確定性的影響。相比之下,采用本文方法則更為穩(wěn)健,其預(yù)測結(jié)果的MAE最小且遠(yuǎn)小于其他幾種方法對應(yīng)的MAE,由圖6可以看出,采用本文方法得到的預(yù)測結(jié)果更接近實(shí)際值。

    圖5 H=10時(shí),樣本1采用本文方法的全部測試誤差Fig.5 While H=10,all the testing errors of sample 1 by the proposed method

    表2 樣本1采用幾種常用方法的預(yù)測結(jié)果

    圖6 樣本1的預(yù)測結(jié)果Fig.6 Predicted results of sample 1

    同樣,樣本2~5在狀態(tài)2和狀態(tài)3下的測試結(jié)果見圖7。圖7a中,樣本2與樣本1類似,其退化模式與訓(xùn)練樣本的群體退化模式差異較大,采用傳統(tǒng)SVR方法和SbRP方法得到的預(yù)測結(jié)果與實(shí)際值相差較大,采用SbDA方法得到的預(yù)測結(jié)果較傳統(tǒng)SVR方法和SbRP方法的預(yù)測結(jié)果相對更準(zhǔn)確一些,但這三種方法的預(yù)測結(jié)果均不如采用本文方法準(zhǔn)確。圖7b~圖7d中,采用傳統(tǒng)SVR方法得到的預(yù)測結(jié)果比采用SbRP方法和SbDA方法準(zhǔn)確,更為接近實(shí)際值,但不如采用本文方法準(zhǔn)確。此外,表3詳細(xì)列出了所有測試樣本采用這幾種方法的預(yù)測誤差,傳統(tǒng)SVR方法、SbRP方法和SbDA方法互有優(yōu)劣,相比之下,本文方法更加準(zhǔn)確有效,且更為穩(wěn)健。

    (a) 樣本2

    表3 本文方法與幾種常用方法的預(yù)測誤差對比結(jié)果

    3 結(jié)論

    本文提出了一種基于整群抽樣和支持向量回歸(SVR)模型的剩余使用壽命(RUL)預(yù)測方法,對測試樣本關(guān)鍵預(yù)警階段觀測數(shù)據(jù)進(jìn)行多次整群抽樣后用于SVR模型測試,以尋找對后期數(shù)據(jù)擬合更好的模型參數(shù),并構(gòu)建了符合維修決策要求的預(yù)測模型。

    (1)相比于傳統(tǒng)SVR方法在模型參數(shù)尋優(yōu)過程中側(cè)重于使得測試樣本的整體誤差最小,本文方法在參數(shù)尋優(yōu)過程中更注重對關(guān)鍵預(yù)警階段數(shù)據(jù)的擬合,可保證所得預(yù)測模型在關(guān)鍵預(yù)警階段的輸出最優(yōu),進(jìn)而有效地提高了關(guān)鍵預(yù)警階段剩余壽命預(yù)測結(jié)果的準(zhǔn)確性和可靠性。

    (2)通過實(shí)驗(yàn)發(fā)現(xiàn),本文方法相比傳統(tǒng)SVR方法及兩種基于相似性的方法具有更好的準(zhǔn)確性和可靠性,能有效用于支持激光器的維護(hù)決策,也為提高RUL預(yù)測方法的性能提供了一條新的技術(shù)途徑。

    (3)相比于重新設(shè)計(jì)機(jī)器學(xué)習(xí)算法中的函數(shù),通過整群抽樣來優(yōu)化提高模型預(yù)測效果更便于在實(shí)際故障預(yù)測與健康管理中實(shí)施,值得進(jìn)一步研究整群抽樣方法在其他數(shù)據(jù)驅(qū)動(dòng)的RUL預(yù)測方法中的應(yīng)用。

    猜你喜歡
    整群激光器階段
    肥皂泡制成微小激光器?
    軍事文摘(2024年4期)2024-03-19 09:40:02
    階梯整群隨機(jī)試驗(yàn)
    基于整群隨機(jī)樣本評估的簡單隨機(jī)抽樣精度探討
    關(guān)于基礎(chǔ)教育階段實(shí)驗(yàn)教學(xué)的幾點(diǎn)看法
    激光器發(fā)明60周年
    科學(xué)(2020年6期)2020-02-06 09:00:06
    在學(xué)前教育階段,提前搶跑,只能跑得快一時(shí),卻跑不快一生。
    莫愁(2019年36期)2019-11-13 20:26:16
    為什么整群豬圍攻一頭豬
    一體化半導(dǎo)體激光器的ANSYS熱仿真及結(jié)構(gòu)設(shè)計(jì)
    基于注入鎖定法激光器的研究
    大熱的O2O三個(gè)階段,你在哪?
    營銷界(2015年22期)2015-02-28 22:05:18
    成人精品一区二区免费| 国产97色在线日韩免费| 在线免费观看的www视频| 亚洲午夜精品一区,二区,三区| 午夜福利成人在线免费观看| 国产一区二区三区在线臀色熟女| 亚洲 国产 在线| 亚洲欧美日韩高清在线视频| 亚洲成av人片免费观看| www.999成人在线观看| 亚洲 欧美 日韩 在线 免费| 国产av一区在线观看免费| 亚洲av熟女| 搡老岳熟女国产| 亚洲色图av天堂| 极品教师在线免费播放| 精品国产超薄肉色丝袜足j| 精品一区二区三区av网在线观看| 欧美日韩一级在线毛片| 狠狠狠狠99中文字幕| 亚洲激情在线av| 一区二区三区激情视频| 人人澡人人妻人| 青草久久国产| 午夜a级毛片| 国产精品乱码一区二三区的特点| 亚洲一卡2卡3卡4卡5卡精品中文| 91麻豆精品激情在线观看国产| 757午夜福利合集在线观看| 国产精品美女特级片免费视频播放器 | 在线观看舔阴道视频| 亚洲一码二码三码区别大吗| 亚洲三区欧美一区| 中文字幕久久专区| 亚洲一卡2卡3卡4卡5卡精品中文| 91成人精品电影| 50天的宝宝边吃奶边哭怎么回事| 高清毛片免费观看视频网站| 久久久久久久久久黄片| 精品福利观看| 男女那种视频在线观看| 禁无遮挡网站| 国产av一区在线观看免费| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产综合亚洲| 久久久久国产一级毛片高清牌| 日韩av在线大香蕉| 亚洲片人在线观看| 亚洲一区二区三区不卡视频| 桃色一区二区三区在线观看| 国产精品日韩av在线免费观看| 国产爱豆传媒在线观看 | 国产色视频综合| 成人免费观看视频高清| 此物有八面人人有两片| 好看av亚洲va欧美ⅴa在| 日本撒尿小便嘘嘘汇集6| 亚洲国产欧美一区二区综合| 少妇被粗大的猛进出69影院| 在线观看www视频免费| netflix在线观看网站| 成人三级做爰电影| 99riav亚洲国产免费| 欧美黑人欧美精品刺激| 97碰自拍视频| 香蕉久久夜色| 老司机午夜福利在线观看视频| 国产国语露脸激情在线看| 久久久久久久精品吃奶| 久久国产精品男人的天堂亚洲| 亚洲人成网站在线播放欧美日韩| 亚洲精品av麻豆狂野| 法律面前人人平等表现在哪些方面| 在线观看日韩欧美| 久久久精品欧美日韩精品| 久久久国产欧美日韩av| 亚洲午夜精品一区,二区,三区| 欧美久久黑人一区二区| 操出白浆在线播放| 亚洲九九香蕉| avwww免费| 又黄又爽又免费观看的视频| 法律面前人人平等表现在哪些方面| 亚洲av第一区精品v没综合| 免费在线观看完整版高清| 伊人久久大香线蕉亚洲五| 国产精品1区2区在线观看.| 日韩大尺度精品在线看网址| 免费在线观看黄色视频的| 亚洲一区二区三区色噜噜| 法律面前人人平等表现在哪些方面| 国产精华一区二区三区| 搡老妇女老女人老熟妇| 91麻豆精品激情在线观看国产| 精品国产亚洲在线| 国产熟女xx| 国产在线精品亚洲第一网站| 99热这里只有精品一区 | 99精品久久久久人妻精品| or卡值多少钱| 男女床上黄色一级片免费看| 婷婷六月久久综合丁香| 欧美又色又爽又黄视频| 麻豆久久精品国产亚洲av| 少妇 在线观看| 天天一区二区日本电影三级| 99国产精品99久久久久| 精品久久久久久久久久免费视频| 国产视频内射| 久久天躁狠狠躁夜夜2o2o| 久久中文字幕人妻熟女| 国产一区二区激情短视频| 亚洲第一青青草原| 男人舔女人下体高潮全视频| 老司机靠b影院| 又紧又爽又黄一区二区| 超碰成人久久| 丁香六月欧美| 国产亚洲精品一区二区www| 夜夜躁狠狠躁天天躁| 国内精品久久久久精免费| 级片在线观看| 久久久久国产一级毛片高清牌| 在线看三级毛片| 国产精品久久视频播放| 成人亚洲精品一区在线观看| 免费高清视频大片| 日韩欧美国产一区二区入口| 亚洲专区国产一区二区| 国产伦一二天堂av在线观看| 特大巨黑吊av在线直播 | 搡老熟女国产l中国老女人| 亚洲色图 男人天堂 中文字幕| 免费看美女性在线毛片视频| 久久精品影院6| 国产精品二区激情视频| 巨乳人妻的诱惑在线观看| 国产精品综合久久久久久久免费| 久热爱精品视频在线9| 老司机深夜福利视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一区中文字幕在线| 在线免费观看的www视频| 天天躁狠狠躁夜夜躁狠狠躁| 嫩草影视91久久| 亚洲全国av大片| 熟女电影av网| 婷婷精品国产亚洲av| 不卡av一区二区三区| 真人做人爱边吃奶动态| 最新美女视频免费是黄的| svipshipincom国产片| 亚洲在线自拍视频| 无人区码免费观看不卡| 久久天躁狠狠躁夜夜2o2o| 中文字幕另类日韩欧美亚洲嫩草| 变态另类丝袜制服| 国产精品一区二区三区四区久久 | 国产午夜精品久久久久久| 岛国视频午夜一区免费看| 他把我摸到了高潮在线观看| 色综合亚洲欧美另类图片| 在线观看一区二区三区| av中文乱码字幕在线| 国产高清视频在线播放一区| 亚洲中文字幕日韩| 国语自产精品视频在线第100页| 久久国产精品影院| 成人永久免费在线观看视频| 国产99久久九九免费精品| 日韩大尺度精品在线看网址| 亚洲第一电影网av| 女警被强在线播放| 欧美色视频一区免费| 成年人黄色毛片网站| 久久精品91无色码中文字幕| 少妇熟女aⅴ在线视频| 国产v大片淫在线免费观看| 日韩有码中文字幕| 久久香蕉激情| 国产真实乱freesex| 制服诱惑二区| 在线观看免费午夜福利视频| 给我免费播放毛片高清在线观看| 中国美女看黄片| 18美女黄网站色大片免费观看| cao死你这个sao货| 搡老妇女老女人老熟妇| 日日干狠狠操夜夜爽| 两个人视频免费观看高清| 香蕉丝袜av| 日本免费a在线| 免费女性裸体啪啪无遮挡网站| 少妇裸体淫交视频免费看高清 | √禁漫天堂资源中文www| 首页视频小说图片口味搜索| e午夜精品久久久久久久| tocl精华| 满18在线观看网站| 精品免费久久久久久久清纯| 日韩欧美免费精品| 欧美在线一区亚洲| 午夜福利在线观看吧| 国产成人精品久久二区二区91| 亚洲一码二码三码区别大吗| 男女之事视频高清在线观看| 99国产综合亚洲精品| 亚洲av第一区精品v没综合| 最近最新中文字幕大全免费视频| 午夜亚洲福利在线播放| 精品久久久久久久毛片微露脸| 黑丝袜美女国产一区| 成人免费观看视频高清| 午夜亚洲福利在线播放| 夜夜爽天天搞| 桃红色精品国产亚洲av| 欧美成狂野欧美在线观看| 久久久久九九精品影院| 在线观看免费日韩欧美大片| 国产精品美女特级片免费视频播放器 | 国产精品 国内视频| 中文字幕另类日韩欧美亚洲嫩草| avwww免费| 韩国精品一区二区三区| 波多野结衣高清无吗| 亚洲欧美日韩无卡精品| 99热只有精品国产| 国产精品爽爽va在线观看网站 | 久久久久九九精品影院| 在线观看免费午夜福利视频| 成人亚洲精品av一区二区| 国产精品电影一区二区三区| 欧美国产精品va在线观看不卡| 一边摸一边做爽爽视频免费| 亚洲精品国产区一区二| 12—13女人毛片做爰片一| 国产一区二区三区视频了| 人人妻人人澡欧美一区二区| 国产精品久久久久久精品电影 | www.精华液| xxx96com| netflix在线观看网站| 亚洲专区国产一区二区| 午夜免费成人在线视频| av在线天堂中文字幕| 国产精品av久久久久免费| 啪啪无遮挡十八禁网站| 国产亚洲av嫩草精品影院| 日本撒尿小便嘘嘘汇集6| 高潮久久久久久久久久久不卡| 亚洲在线自拍视频| 悠悠久久av| 久久青草综合色| 99riav亚洲国产免费| 国产精品久久久久久人妻精品电影| 亚洲av电影在线进入| 亚洲aⅴ乱码一区二区在线播放 | 国产高清有码在线观看视频 | 精品国产一区二区三区四区第35| 国产精品日韩av在线免费观看| 久久精品91无色码中文字幕| 久久久久久亚洲精品国产蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 色老头精品视频在线观看| 97碰自拍视频| 亚洲中文字幕一区二区三区有码在线看 | 国产亚洲av高清不卡| 日韩免费av在线播放| 宅男免费午夜| 久久人妻福利社区极品人妻图片| 又大又爽又粗| 欧美丝袜亚洲另类 | 精品午夜福利视频在线观看一区| 午夜免费鲁丝| 在线观看66精品国产| 巨乳人妻的诱惑在线观看| 国产精品爽爽va在线观看网站 | 国产亚洲欧美98| 啦啦啦观看免费观看视频高清| www日本在线高清视频| 亚洲一区二区三区色噜噜| 精华霜和精华液先用哪个| 99精品欧美一区二区三区四区| 国产成人精品久久二区二区91| 青草久久国产| 欧美一级毛片孕妇| 国内精品久久久久精免费| 老熟妇乱子伦视频在线观看| 久久欧美精品欧美久久欧美| 免费高清视频大片| 一边摸一边抽搐一进一小说| 日本撒尿小便嘘嘘汇集6| 国产又黄又爽又无遮挡在线| 男男h啪啪无遮挡| 真人做人爱边吃奶动态| 日本熟妇午夜| www日本在线高清视频| 国产精品1区2区在线观看.| 999精品在线视频| 久久婷婷人人爽人人干人人爱| 国产精品自产拍在线观看55亚洲| 亚洲五月天丁香| 国产单亲对白刺激| 麻豆一二三区av精品| 动漫黄色视频在线观看| 色在线成人网| 97人妻精品一区二区三区麻豆 | av在线天堂中文字幕| 一个人观看的视频www高清免费观看 | 又紧又爽又黄一区二区| 2021天堂中文幕一二区在线观 | 亚洲激情在线av| 亚洲人成77777在线视频| 99久久无色码亚洲精品果冻| 99久久国产精品久久久| 日本三级黄在线观看| 韩国av一区二区三区四区| 每晚都被弄得嗷嗷叫到高潮| 国产精品香港三级国产av潘金莲| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美一区二区三区黑人| 日韩高清综合在线| 精品久久久久久成人av| 国产熟女xx| 国产爱豆传媒在线观看 | 美女免费视频网站| 在线天堂中文资源库| 男人的好看免费观看在线视频 | 禁无遮挡网站| 国产久久久一区二区三区| 精品卡一卡二卡四卡免费| 精品国产亚洲在线| 国产精品香港三级国产av潘金莲| 日韩欧美免费精品| 可以在线观看的亚洲视频| 国产精品永久免费网站| 欧美乱码精品一区二区三区| 嫩草影视91久久| 久久久精品欧美日韩精品| 亚洲成国产人片在线观看| 怎么达到女性高潮| 国产精品一区二区精品视频观看| 亚洲国产精品sss在线观看| 国产一区二区在线av高清观看| 国产精品美女特级片免费视频播放器 | 久久国产亚洲av麻豆专区| 麻豆久久精品国产亚洲av| 天堂√8在线中文| 午夜精品久久久久久毛片777| 男女那种视频在线观看| 亚洲精品av麻豆狂野| 午夜两性在线视频| 欧美午夜高清在线| 免费人成视频x8x8入口观看| 亚洲一码二码三码区别大吗| 国内精品久久久久久久电影| 狂野欧美激情性xxxx| 可以在线观看毛片的网站| 午夜福利高清视频| 日本精品一区二区三区蜜桃| www.www免费av| 一本精品99久久精品77| 精品少妇一区二区三区视频日本电影| 少妇被粗大的猛进出69影院| 欧美成人性av电影在线观看| 在线av久久热| 美女国产高潮福利片在线看| 国产成人啪精品午夜网站| 成人精品一区二区免费| 精品电影一区二区在线| 日韩欧美三级三区| 久久性视频一级片| 黄片播放在线免费| 又紧又爽又黄一区二区| 一a级毛片在线观看| 日日干狠狠操夜夜爽| 国产亚洲欧美精品永久| 757午夜福利合集在线观看| 一区二区日韩欧美中文字幕| 欧美黑人巨大hd| 日韩欧美在线二视频| 变态另类成人亚洲欧美熟女| 精品国产美女av久久久久小说| 久久人妻福利社区极品人妻图片| 精品一区二区三区四区五区乱码| 91九色精品人成在线观看| 国产精品影院久久| 国产欧美日韩精品亚洲av| 黄色片一级片一级黄色片| netflix在线观看网站| 亚洲国产看品久久| a在线观看视频网站| 精品久久久久久久久久免费视频| 在线观看www视频免费| 色综合亚洲欧美另类图片| 听说在线观看完整版免费高清| 国产精华一区二区三区| 国产精品永久免费网站| 欧美 亚洲 国产 日韩一| 日韩成人在线观看一区二区三区| 国产精品爽爽va在线观看网站 | 精品日产1卡2卡| 老司机在亚洲福利影院| 免费av毛片视频| 别揉我奶头~嗯~啊~动态视频| 亚洲av中文字字幕乱码综合 | 亚洲五月天丁香| 九色国产91popny在线| 亚洲专区中文字幕在线| or卡值多少钱| 精品欧美一区二区三区在线| 又紧又爽又黄一区二区| 国产精品,欧美在线| 大香蕉久久成人网| 国产片内射在线| 欧美不卡视频在线免费观看 | 欧美国产精品va在线观看不卡| 亚洲精品久久成人aⅴ小说| av超薄肉色丝袜交足视频| 91大片在线观看| 亚洲av片天天在线观看| 巨乳人妻的诱惑在线观看| 亚洲五月天丁香| 淫妇啪啪啪对白视频| 日本五十路高清| 99久久无色码亚洲精品果冻| 最新美女视频免费是黄的| 亚洲专区中文字幕在线| 国产亚洲精品久久久久5区| 国产精品久久电影中文字幕| 老汉色∧v一级毛片| 淫秽高清视频在线观看| 午夜福利成人在线免费观看| 熟女电影av网| 亚洲无线在线观看| 精品欧美一区二区三区在线| 亚洲熟女毛片儿| 欧美性猛交黑人性爽| 免费看日本二区| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 欧美中文日本在线观看视频| 日本成人三级电影网站| 女性生殖器流出的白浆| 久久久久国产一级毛片高清牌| 亚洲专区国产一区二区| 88av欧美| 夜夜夜夜夜久久久久| 国产精品1区2区在线观看.| 亚洲精品美女久久久久99蜜臀| svipshipincom国产片| av天堂在线播放| 在线永久观看黄色视频| 51午夜福利影视在线观看| 精品第一国产精品| 精品人妻1区二区| 亚洲一码二码三码区别大吗| 亚洲欧美日韩高清在线视频| 亚洲自拍偷在线| 日韩高清综合在线| 热re99久久国产66热| 日韩欧美一区视频在线观看| 欧美色视频一区免费| 亚洲av成人一区二区三| 国产蜜桃级精品一区二区三区| 日本一区二区免费在线视频| 国语自产精品视频在线第100页| 亚洲精品美女久久久久99蜜臀| 免费看日本二区| 一区二区三区激情视频| a级毛片a级免费在线| 精品欧美一区二区三区在线| 欧美三级亚洲精品| 欧美成人性av电影在线观看| 99精品在免费线老司机午夜| 国产又色又爽无遮挡免费看| 99精品久久久久人妻精品| 两性午夜刺激爽爽歪歪视频在线观看 | 国内精品久久久久精免费| 亚洲成人精品中文字幕电影| 成年人黄色毛片网站| 欧美性长视频在线观看| 久久精品影院6| 亚洲男人的天堂狠狠| 91字幕亚洲| 午夜免费激情av| 好看av亚洲va欧美ⅴa在| 国产91精品成人一区二区三区| 男人操女人黄网站| 成人永久免费在线观看视频| 美女高潮喷水抽搐中文字幕| 国产精品亚洲av一区麻豆| 欧美日韩黄片免| 亚洲七黄色美女视频| 久久国产亚洲av麻豆专区| 日韩欧美一区视频在线观看| 美女 人体艺术 gogo| 欧美性长视频在线观看| 变态另类丝袜制服| 日本在线视频免费播放| 女人高潮潮喷娇喘18禁视频| 成人手机av| 岛国在线观看网站| 久久久久久久午夜电影| 法律面前人人平等表现在哪些方面| 一区二区三区高清视频在线| 亚洲av熟女| 色播在线永久视频| 88av欧美| 欧洲精品卡2卡3卡4卡5卡区| 韩国精品一区二区三区| www.www免费av| 韩国av一区二区三区四区| 欧美黑人欧美精品刺激| 成人亚洲精品av一区二区| 亚洲 国产 在线| 1024视频免费在线观看| 一本大道久久a久久精品| 国产一区二区在线av高清观看| 国内久久婷婷六月综合欲色啪| 婷婷精品国产亚洲av在线| 这个男人来自地球电影免费观看| 国产免费男女视频| a级毛片a级免费在线| 夜夜躁狠狠躁天天躁| 久久精品国产亚洲av香蕉五月| 97超级碰碰碰精品色视频在线观看| 欧美精品啪啪一区二区三区| 美女大奶头视频| 一夜夜www| 久久草成人影院| 亚洲一区二区三区不卡视频| 成人国产一区最新在线观看| www.999成人在线观看| a在线观看视频网站| √禁漫天堂资源中文www| 国产精品影院久久| 后天国语完整版免费观看| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 欧美最黄视频在线播放免费| 日韩欧美三级三区| 黄色丝袜av网址大全| 国产又爽黄色视频| www国产在线视频色| 日韩欧美国产一区二区入口| 中文字幕人妻熟女乱码| 熟女电影av网| 久久精品91无色码中文字幕| 欧美不卡视频在线免费观看 | 国产精品二区激情视频| 一本精品99久久精品77| 男女做爰动态图高潮gif福利片| 婷婷丁香在线五月| 亚洲成人精品中文字幕电影| 欧美色欧美亚洲另类二区| 午夜视频精品福利| 欧美色视频一区免费| 日本成人三级电影网站| 人成视频在线观看免费观看| 一本大道久久a久久精品| 国产精品精品国产色婷婷| 亚洲中文日韩欧美视频| 国产av不卡久久| 最近最新中文字幕大全电影3 | 国产精品爽爽va在线观看网站 | 久热爱精品视频在线9| 成人午夜高清在线视频 | 欧美激情极品国产一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 国产一卡二卡三卡精品| 亚洲国产中文字幕在线视频| 久久久国产欧美日韩av| 精品一区二区三区四区五区乱码| 日本免费一区二区三区高清不卡| 欧美色视频一区免费| 身体一侧抽搐| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲一级av第二区| 久久亚洲真实| 久热爱精品视频在线9| 久久天躁狠狠躁夜夜2o2o| 国产精品1区2区在线观看.| 亚洲五月天丁香| 一级黄色大片毛片| 免费看十八禁软件| 亚洲精品一卡2卡三卡4卡5卡| 久久精品成人免费网站| 欧美日韩瑟瑟在线播放| 国产精品av久久久久免费| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久, | 亚洲国产精品sss在线观看| 欧美丝袜亚洲另类 | 亚洲 欧美 日韩 在线 免费| 这个男人来自地球电影免费观看| 手机成人av网站| 亚洲七黄色美女视频| 又黄又爽又免费观看的视频| 亚洲一码二码三码区别大吗| 婷婷亚洲欧美| 1024视频免费在线观看| 亚洲熟女毛片儿| 成年人黄色毛片网站| 成人av一区二区三区在线看| 99久久国产精品久久久| 久久久久亚洲av毛片大全| 亚洲成国产人片在线观看| 久久人人精品亚洲av| e午夜精品久久久久久久| 国产野战对白在线观看| 女人被狂操c到高潮| 美女大奶头视频| 精品久久久久久成人av|