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

    軌道車輛一系垂向減振器劣化辨識(shí)方法研究

    2021-05-12 02:40:24姜良奎盧昌宏方柳川
    機(jī)械 2021年3期
    關(guān)鍵詞:劣化減振器粒子

    姜良奎,盧昌宏,方柳川

    軌道車輛一系垂向減振器劣化辨識(shí)方法研究

    姜良奎1,盧昌宏2,方柳川2

    (1.中車青島四方機(jī)車車輛股份有限公司,山東 青島 266111;2.西南交通大學(xué)機(jī)械工程學(xué)院,四川 成都 610031)

    減振器是軌道車輛的重要部件,在目前的研究中,缺少對軌道車輛一系垂向減振器劣化階段進(jìn)行辨識(shí)的研究。為了實(shí)現(xiàn)對減振器的劣化階段進(jìn)行辨識(shí),提出了一種基于粒子群優(yōu)化支持向量機(jī)(PSO-SVM)模型的減振器性能劣化辨識(shí)方法。首先通過仿真軟件獲得減振器的加速度信號,提取特征后利用SVM模型對特征進(jìn)行篩選融合,其次使用粒子群尋優(yōu)算法對SVM模型的參數(shù)進(jìn)行尋優(yōu)處理,最后實(shí)現(xiàn)了對減振器服役性能劣化狀態(tài)的辨識(shí)。通過實(shí)驗(yàn)驗(yàn)證了該方法的有效性,實(shí)驗(yàn)結(jié)果表明:SVM模型對減振器進(jìn)行劣化辨識(shí)具有可行性,采用尋優(yōu)算法能夠提升方法的辨識(shí)率。

    減振器;劣化辨識(shí);PSO;SVM

    減振器作為高速列車轉(zhuǎn)向架的重要部件,在列車行駛過程中起著保護(hù)車體、緩沖減震的作用。減振器性能直接影響列車行駛的平順性、穩(wěn)定性和安全性,對其性能退化評估可以更為有效地服務(wù)于設(shè)備的主動(dòng)維護(hù)[1]。因此研究一種高效、準(zhǔn)確的智能診斷方法具有重要意義。

    在目前的一些機(jī)械設(shè)備,如刀具等的性能劣化狀態(tài)識(shí)別研究過程中,支持向量機(jī)、人工神經(jīng)網(wǎng)絡(luò)、隱馬爾可夫和深度學(xué)習(xí)等模型被廣泛應(yīng)用。

    支持向量機(jī)(Support Vector Machine,SVM)是一種監(jiān)督式機(jī)器學(xué)習(xí)模型,泛化能力強(qiáng)、非線性映射能力強(qiáng),目前在分類和回歸分析問題上被大量采用。吳志丹等[2]采用支持向量機(jī)對橫向減振器的故障狀態(tài)進(jìn)行分類識(shí)別;關(guān)山等[3]在研究刀具磨損狀態(tài)的過程中,對聲發(fā)射信號進(jìn)行處理,采用LS-SVM模型對處理后的數(shù)據(jù)進(jìn)行狀態(tài)辨識(shí);Qian等[4]為了建立工件表面紋理特征與刀具磨損之間的耦合關(guān)系,選擇通過投影法和Gabor濾波法提取特征,提取后將特征輸入到SVM模型中得到辨識(shí)結(jié)果。

    人工神經(jīng)網(wǎng)絡(luò)具有魯棒性良好、學(xué)習(xí)特征信息能力強(qiáng)、模型辨識(shí)穩(wěn)定等多個(gè)優(yōu)點(diǎn)。聶鵬等[5]建立了隱含層節(jié)點(diǎn)為9個(gè)的Elman神經(jīng)網(wǎng)絡(luò)對刀具磨損狀態(tài)進(jìn)行辨識(shí)。

    隱馬爾科夫模型(Hidden Markov Model,HMM)具有推理性好、數(shù)據(jù)結(jié)構(gòu)嚴(yán)謹(jǐn)和穩(wěn)定性良好等優(yōu)點(diǎn)。Dong等[6-7]在進(jìn)行液壓泵健康狀態(tài)辨識(shí)時(shí)采用了隱馬爾可夫模型;李威霖等[8]在樣本數(shù)較少的情況下建立了一個(gè)無監(jiān)督的隱馬爾可夫刀具磨損狀態(tài)辨識(shí)模型。

    深度學(xué)習(xí)模型具有可以處理復(fù)雜數(shù)據(jù)、自動(dòng)提取特征和泛化能力強(qiáng)等優(yōu)點(diǎn)。石朝[9]提出了一種基于深度學(xué)習(xí)長短時(shí)記憶網(wǎng)絡(luò)的刀具磨損階段辨識(shí)方法;吳昀璞等[10]針對模型泛化能力不足與魯棒性不足的問題,提出一種基于多域融合卷積神經(jīng)網(wǎng)絡(luò)的轉(zhuǎn)向架故障檢測方法。

    在目前的性能退化辨識(shí)研究中,軌道車輛一系垂向減振器劣化辨識(shí)存在缺少辨識(shí)方法、辨識(shí)較為困難等問題。因此,提出了一種軌道車輛一系垂向減振器性能劣化狀態(tài)辨識(shí)方法,并通過實(shí)驗(yàn)驗(yàn)證了該方法的有效性。

    1 辨識(shí)模型的構(gòu)建

    1.1 辨識(shí)方法與實(shí)驗(yàn)流程

    支持向量機(jī)分類方法具有良好的性能,在狀態(tài)辨識(shí)領(lǐng)域被廣泛應(yīng)用。SVM模型對非線性樣本有著很好的分類效果,在小樣本時(shí)也能保持較好的分類效果,對比其他方法具有明顯的優(yōu)勢,所以選擇支持向量機(jī)作為軌道車輛一系垂向減振器劣化辨識(shí)模型。同時(shí)SVM模型中存在對辨識(shí)效果影響較大兩組參數(shù),懲罰因子與核函數(shù)參數(shù),然而這兩種參數(shù)的確定具有較大隨機(jī)性。因此,選擇參數(shù)尋優(yōu)算法克服參數(shù)選擇問題上存在的缺點(diǎn)。最終構(gòu)成一種基于參數(shù)尋優(yōu)算法優(yōu)化SVM的一系垂向減振器劣化辨識(shí)方法。

    提出的減振器服役性能劣化狀態(tài)辨識(shí)方法的主要步驟如下:

    (1)通過仿真軟件模擬仿真減振器不同劣化階段,并采集對應(yīng)的加速度信號;

    (2)提取信號的時(shí)域和頻域特征,利用SVM模型測試每個(gè)通道的辨識(shí)精度,選取精度較高通道的數(shù)據(jù)特征組合形成多特征聯(lián)合向量;

    (3)劃分訓(xùn)練樣本與測試樣本,將訓(xùn)練樣本的多特征聯(lián)合向量輸入SVM模型進(jìn)行辨識(shí)訓(xùn)練;

    (4)利用參數(shù)尋優(yōu)算法對關(guān)鍵參數(shù)和進(jìn)行尋優(yōu)處理,獲得最優(yōu)的辨識(shí)精度;

    (5)對測試樣本進(jìn)行劣化辨識(shí),獲得辨識(shí)結(jié)果。

    劣化狀態(tài)辨識(shí)流程如圖1所示。

    1.2 優(yōu)化算法的選擇

    在SVM模型中,懲罰因子和核函數(shù)參數(shù)對提高模型的辨識(shí)精度起著非常重要的作用,但是這兩種參數(shù)的確定卻具有較大的隨機(jī)性,對SVM模型的辨識(shí)精度具有不良影響。因此,需要對這兩個(gè)參數(shù)進(jìn)行尋優(yōu)處理。

    在現(xiàn)有的尋優(yōu)處理手段中,傳統(tǒng)的試湊法具有較強(qiáng)的隨機(jī)性,需要花費(fèi)研究人員大量的時(shí)間和精力;簡單的網(wǎng)格尋優(yōu)法,雖然有時(shí)候效果比較好,但需要長時(shí)間的等待;除此之外,還有蟻群算法(Ant Colony Optimization,ACO)、遺傳算法(Genetic Algorithm,GA)、模擬退火算法(Simulated Annealing,SA)和粒子群算法(Particle Swarm Optimization,PSO)被廣泛應(yīng)用,這些算法具有自動(dòng)尋優(yōu)的優(yōu)點(diǎn),能節(jié)省大量的人力和時(shí)間,并減少人為因素的影響。將四種尋優(yōu)算法的優(yōu)缺點(diǎn)列入表1中進(jìn)行對比,選擇合適尋優(yōu)算法優(yōu)化SVM模型。

    圖1 辨識(shí)流程圖

    表1 尋優(yōu)算法對比

    粒子群算法具有高效、結(jié)構(gòu)簡單等優(yōu)點(diǎn),同時(shí)能夠有效地對懲罰因子和核函數(shù)參數(shù)進(jìn)行優(yōu)化。因此,選擇粒子群優(yōu)化算法對SVM模型參數(shù)進(jìn)行尋優(yōu),最終構(gòu)建出基于PSO- SVM模型的減振器性能劣化狀態(tài)辨識(shí)方法。

    2 實(shí)驗(yàn)驗(yàn)證

    2.1 數(shù)據(jù)集情況

    原始數(shù)據(jù)來源于SIMPack仿真軟件對建立的高速動(dòng)車組車輛系統(tǒng)動(dòng)力學(xué)模型進(jìn)行仿真的結(jié)果。在某一工況下,將減振器分別設(shè)置為100%、95%、90%、85%、80%、75%和70%共七種劣化狀態(tài)。然后進(jìn)行高速動(dòng)車組的動(dòng)力學(xué)仿真,獲得在這七種劣化狀態(tài)下減振器的振動(dòng)響應(yīng)。

    共采集30個(gè)通道的加速度信號,每個(gè)通道在各狀態(tài)下采集的數(shù)據(jù)點(diǎn)數(shù)為54980個(gè)。以長度1024、步長100的滑窗對數(shù)據(jù)進(jìn)行滑窗采樣,將一個(gè)長度為54980的數(shù)據(jù)處理成540個(gè)樣本,每個(gè)樣本長度為1024。所以單個(gè)通道得到的所有狀態(tài)下的樣本量為7×540=3780個(gè)。總共七種狀態(tài),每種狀態(tài)540個(gè)樣本,每個(gè)樣本包含1024個(gè)數(shù)據(jù)點(diǎn)。每個(gè)通道的3780個(gè)樣本按照4:1的比例劃分訓(xùn)練樣本和測試樣本,訓(xùn)練樣本為3780×0.8=3024個(gè),測試樣本為3780×0.2=756個(gè)。各通道具體情況如表2所示。

    2.2 減振器劣化特征提取與分析

    對每個(gè)樣本提取15個(gè)時(shí)域特征、3個(gè)頻域特征,共18個(gè)特征。以通道4車體前部垂向加速度和通道11構(gòu)架1架中部橫向加速度為例,提取部分時(shí)域特征與部分頻域特征并將其歸一化處理。其部分時(shí)域特征、頻域特征隨時(shí)間的變化趨勢如圖2、圖3所示。

    通過對特征提取的結(jié)果進(jìn)行分析,可以發(fā)現(xiàn):部分通道采集的加速度信號的有量綱時(shí)域特征,如峰峰值、方根幅值等,整體隨時(shí)間的變化而變化的情況較為明顯,可以較好地反映減振器的劣化狀態(tài)。如通道11采集到的加速度信號的有量綱時(shí)域特征與減振器劣化有很強(qiáng)的關(guān)聯(lián)性,而通道4采集到的加速度信號的有量綱時(shí)域特征則不能直觀地反映出減振器的劣化狀態(tài)。

    另外,從圖3可以看出,信號的頻域特征整體均與減振器劣化的關(guān)聯(lián)性不強(qiáng)。

    表2 各通道傳感器所在位置

    圖2 時(shí)域特征隨時(shí)間的變化趨勢

    2.3 聯(lián)合多特征向量的建立

    根據(jù)對特征變化趨勢的分析,并不是所有采集信號的時(shí)域、頻域特征都與減振器劣化過程有較強(qiáng)的關(guān)聯(lián)性,所以需要對這些采集信號進(jìn)行篩選。首先,將各個(gè)通道采集信號的時(shí)域、頻域特征輸入SVM分類模型進(jìn)行辨識(shí),得到各個(gè)通道采集信號對應(yīng)的辨識(shí)精度;其次,依據(jù)辨識(shí)精度的高低篩選出辨識(shí)效果好的通道;最后,將其特征組合在一起形成多特征聯(lián)合向量,再次輸入并訓(xùn)練SVM模型以提高辨識(shí)精度。

    各通道的辨識(shí)結(jié)果如圖4所示,可以看出,對于不同通道采集信號提取的特征,其辨識(shí)精度很大的差別。除安裝在車體與構(gòu)架1架的通道的辨識(shí)精度較高之外,安裝在其余位置的通道采集到的信號與減振器劣化基本沒有關(guān)聯(lián)。并且對于車體和構(gòu)架1架,只有采集到的橫向加速度信號辨識(shí)精度較高,垂向加速度信號同樣與減振器劣化關(guān)聯(lián)度不高。

    根據(jù)圖4的辨識(shí)結(jié)果,選擇辨識(shí)精度高于50%的通道對應(yīng)采集信號的時(shí)域、頻域特征進(jìn)行組合,其中包括通道1、2、3、7、9、11共六個(gè)通道。將其依次相連,組成一個(gè)包含18×6=108個(gè)特征的聯(lián)合多特征向量。

    圖3 頻域特征隨時(shí)間的變化趨勢

    圖4 PSO-SVM模型各通道辨識(shí)精度

    2.4 粒子群算法的尋優(yōu)過程

    粒子群優(yōu)化算法的具體步驟如下:

    (1)確定粒子群算法的適應(yīng)度函數(shù)。利用粒子群算法對SVM模型的和參數(shù)進(jìn)行優(yōu)化是為了提高模型的辨識(shí)精度,因此將粒子群算法的適應(yīng)度函數(shù)設(shè)置為SVM模型的辨識(shí)精度,表達(dá)式為:

    式中:y為訓(xùn)練樣本集中被正確分類的樣本的個(gè)數(shù);y為表示被分類錯(cuò)誤的個(gè)數(shù)。

    (2)劃分樣本。將所有樣本按照4:1的比例劃分為訓(xùn)練樣本與測試樣本;

    (3)初始化粒子群算法的各個(gè)參數(shù)。具體數(shù)值如表3所示。

    表3 粒子群算法參數(shù)設(shè)置

    (4)設(shè)置好參數(shù)尋優(yōu)范圍后,將訓(xùn)練樣本輸入到SVM模型中,獲得兩個(gè)參數(shù)值,并將這兩個(gè)參數(shù)值設(shè)置到二維坐標(biāo)中,每一個(gè)粒子的二維坐標(biāo)便對應(yīng)這兩個(gè)參數(shù)值。通過計(jì)算首先獲得粒子當(dāng)前的適應(yīng)度,其次每一個(gè)粒子會(huì)不斷將當(dāng)前的適應(yīng)度與自身的最優(yōu)適應(yīng)度進(jìn)行對比,從而更新自身最優(yōu)適應(yīng)度,然后通過將所有粒子的自身適應(yīng)度進(jìn)行比較得出當(dāng)前全局的最優(yōu)適應(yīng)度;

    (5)在粒子經(jīng)過步驟(4)后,會(huì)改變自身的位置以尋找新的最優(yōu)適應(yīng)度,其位置更新的速度和位置公式為:

    式中:1、2為兩個(gè)隨機(jī)函數(shù),取值范圍[0, 1],用于增加搜索范圍的隨機(jī)性;P為粒子個(gè)體所經(jīng)歷過的最好位置;P為種群所經(jīng)歷過的最好位置。

    (6)檢查當(dāng)前是否達(dá)到設(shè)置的迭代次數(shù),如果已達(dá)到,則輸出參數(shù)尋優(yōu)的結(jié)果,給出訓(xùn)練模型的最高辨識(shí)精度,及對應(yīng)的參數(shù)和。

    粒子群優(yōu)化算法對SVM模型進(jìn)行尋優(yōu)處理的流程如圖5所示。

    使用粒子群算法對和進(jìn)行參數(shù)尋優(yōu)的精度曲線如圖6所示。可以看出模型的辨識(shí)精度隨著迭代次數(shù)的增加在不斷地提高,最終輸出的辨識(shí)精度可達(dá)到98.8%,和尋優(yōu)的結(jié)果為=15439.651、=0.026。將此參數(shù)設(shè)置為最終的PSO-SVM模型的參數(shù),輸入測試樣本進(jìn)行測試。

    2.5 PSO-SVM模型的辨識(shí)結(jié)果分析

    為了驗(yàn)證利用粒子群算法進(jìn)行參數(shù)尋優(yōu)的有效性,將SVM模型作為對比,將特征選擇并組合后的特征分別輸入SVM模型和PSO-SVM模型中進(jìn)行測試。兩種模型的辨識(shí)結(jié)果如圖7所示,對比可知,PSO-SVM模型的辨識(shí)精度相較于SVM模型有較大的提升。

    pbest.個(gè)體歷史最優(yōu) gbest.全局最優(yōu)

    圖6 粒子群算法精度曲線

    優(yōu)化前后模型對減振器服役性能各個(gè)劣化狀態(tài)正確分類的精度如表4所示??芍?,SVM模型的辨識(shí)精度為94.7%,而PSO-SVM模型的辨識(shí)精度為98.8%。辨識(shí)精度具有一定程度的提升,說明了對SVM模型采用粒子群算法進(jìn)行參數(shù)尋優(yōu)是有效的,可以通過該方法有效地解決SVM模型中參數(shù)和選擇具有較大隨機(jī)性的問題。同時(shí)SVM模型與PSO-SVM模型對減振器性能劣化辨識(shí)狀態(tài)都具有較高的辨識(shí)精度,說明了該方法的可行性。

    1.性能狀態(tài)100% 2.性能狀態(tài)95% 3.性能狀態(tài)90% 4.性能狀態(tài)85% 5.性能狀態(tài)80% 6.性能狀態(tài)75% 7.性能狀態(tài)70%

    3 結(jié)論

    對于軌道車輛一系垂向減振器劣化辨識(shí)存在的缺少辨識(shí)方法、辨識(shí)較為困難等問題,提出了一種基于PSO-SVM模型的辨識(shí)方法,通過實(shí)驗(yàn)驗(yàn)證了該方法的有效性,得出以下結(jié)論:

    (1)將SVM模型應(yīng)用于軌道車輛一系垂向減振器劣化辨識(shí)的研究具有一定可行性;

    (2)利用PSO尋優(yōu)算法對SVM模型中的兩個(gè)參數(shù)進(jìn)行尋優(yōu)處理,獲得的優(yōu)化模型能夠取得更好的辨識(shí)效果。

    基于PSO-SVM模型的減振器服役性能劣化狀態(tài)辨識(shí)方法為軌道車輛一系垂向減振器劣化辨識(shí)提供了一種新的策略。

    表4 參數(shù)尋優(yōu)前后不同狀態(tài)辨識(shí)精度

    [1]井波,金煒東,秦娜,等. 高速列車橫向減振器性能退化的特征提取[J]. 噪聲與振動(dòng)控制,2015,35(2):57-60.

    [2]吳志丹,秦娜,金煒東. 基于EEMD的高速列車橫向減振器故障的排列熵特征分析[J]. 計(jì)算機(jī)科學(xué),2016,43(5):304-307.

    [3]關(guān)山,康振興,彭昶. 基于云理論與LS-SVM的刀具磨損識(shí)別方法[J]. 振動(dòng)、測試與診斷,2017,37(5):996-1003,1068-1069.

    [4]Yiqiu Qian,Jia Tian,Libing Liu,et al. Yingshu Chen.A Tool Wear Predictive Model Based on SVM[C]. Proceedings of 2010 Chinese Control and Decision Conference,2010.

    [5]聶鵬,崔凱奇,何超. 基于Elman神經(jīng)網(wǎng)絡(luò)的刀具磨損狀態(tài)識(shí)別技術(shù)[J]. 現(xiàn)代制造工程,2015(12):78-81,85.

    [6]Ming Dong,David He,Prashant Banerjee,et al. Equipment health diagnosis and prognosis using hidden semi-Markov models[J]. The International Journal of Advanced Manufacturing Technology, 2006,30(7-8):738-749.

    [7]Ming Dong,David He. A segmental hidden semi-Markov model (HSMM)-based diagnostics and prognostics framework and methodology[J]. Mechanical Systems and Signal Processing,2007,21(5):2248-2266.

    [8]李威霖,傅攀,李曉暉. 無先驗(yàn)知識(shí)下基于CHMM的刀具磨損監(jiān)測技術(shù)[J]. 機(jī)床與液壓,2013,41(15):37-41.

    [9]石朝. 基于力和振動(dòng)融合的刀具磨損狀態(tài)辨識(shí)方法研究[D]. 成都:西南交通大學(xué),2017.

    [10]吳昀璞,金煒東,黃穎坤. 基于多域融合CNN的高速列車轉(zhuǎn)向架故障檢測[J]. 系統(tǒng)仿真學(xué)報(bào),2018,30(11):4492-4497.

    Research on Deterioration Identification Method of the Primary Vertical Shock Absorber of Rail Vehicles

    JIANG Liangkui1,LU Changhong2,F(xiàn)ANG Liuchuan2

    (1.CRRC Qingdao Sifang Co. Ltd., Qingdao 266111, China; 2.School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China)

    Shock absorber is an important part of rail vehicles, while the research on identifying the deterioration stage of the primary vertical shock absorbers of rail vehicles is insufficient. In order to realize the identification of the deterioration stage of the shock absorber, a method for identifying the performance degradation of shock absorbers based on the PSO-SVM model is proposed. First, the acceleration signal of the shock absorber is obtained through the simulation software. After extracting the features, the SVM model is used to screen and fuse the features. Then, the particle swarm optimization algorithm is used to optimize the parameters of the SVM model. Finally, the identification of the deterioration state of the service performance of the shock absorber is realized. The effectiveness of the method is verified by the experiment, and the results show that: The SVM model is feasible to identify the deterioration of the shock absorber. The optimization algorithm can improve the identification rate of the method.

    shock absorber;deterioration identification;PSO;SVM

    U279.323

    A

    10.3969/j.issn.1006-0316.2021.03.006

    1006-0316 (2021) 03-0032-07

    2020-08-26

    國家重點(diǎn)研發(fā)計(jì)劃(2017YFB1201201)

    姜良奎(1974-),男,山東臨沂人,高級工程師,主要研究方向?yàn)檐囕v人機(jī)工程、設(shè)計(jì)虛擬評價(jià)等,E-mail:jlk@cqsf.com。

    猜你喜歡
    劣化減振器粒子
    基于S形試件五軸數(shù)控機(jī)床動(dòng)態(tài)性能劣化評價(jià)
    基于粒子群優(yōu)化的橋式起重機(jī)模糊PID控制
    基于粒子群優(yōu)化極點(diǎn)配置的空燃比輸出反饋控制
    跨海橋梁混凝土劣化分析及預(yù)養(yǎng)護(hù)
    上海公路(2017年1期)2017-07-21 13:38:33
    八鋼2500m3高爐適應(yīng)焦炭質(zhì)量劣化的生產(chǎn)實(shí)踐
    新疆鋼鐵(2016年3期)2016-02-28 19:18:52
    汽車減振器與磁流變材料分析與運(yùn)算
    饋能式磁流變減振器自供電特性研究
    汽車科技(2014年6期)2014-03-11 17:45:36
    裂紋齒輪嚙合剛度的劣化特性
    裝多級扭轉(zhuǎn)減振器的發(fā)動(dòng)機(jī)曲軸系統(tǒng)扭振分析
    基于Matlab的α粒子的散射實(shí)驗(yàn)?zāi)M
    物理與工程(2014年4期)2014-02-27 11:23:08
    久久久久精品性色| av.在线天堂| 婷婷色综合大香蕉| 日本av手机在线免费观看| 春色校园在线视频观看| 国产亚洲av片在线观看秒播厂| 成人一区二区视频在线观看| 日韩,欧美,国产一区二区三区| 免费观看在线日韩| xxx大片免费视频| 高清不卡的av网站| 亚洲精品乱码久久久v下载方式| 欧美高清成人免费视频www| 高清午夜精品一区二区三区| 亚洲怡红院男人天堂| 国产黄片视频在线免费观看| 国产一级毛片在线| 国产成人一区二区在线| 亚洲精品久久久久久婷婷小说| 精品午夜福利在线看| 性色avwww在线观看| 欧美 日韩 精品 国产| 人人妻人人爽人人添夜夜欢视频 | 看十八女毛片水多多多| 国产精品av视频在线免费观看| 偷拍熟女少妇极品色| 国产片特级美女逼逼视频| 国产色爽女视频免费观看| 人妻制服诱惑在线中文字幕| 亚洲怡红院男人天堂| 男男h啪啪无遮挡| 久热久热在线精品观看| 亚洲色图av天堂| 看非洲黑人一级黄片| 日本av手机在线免费观看| 久久精品夜色国产| 亚洲人与动物交配视频| 日韩人妻高清精品专区| av.在线天堂| 欧美日韩在线观看h| 少妇人妻久久综合中文| 色哟哟·www| 亚洲婷婷狠狠爱综合网| 欧美成人一区二区免费高清观看| 亚洲欧美日韩卡通动漫| 少妇猛男粗大的猛烈进出视频| 高清欧美精品videossex| 国产成人aa在线观看| 99久久综合免费| 午夜福利影视在线免费观看| 亚洲欧美清纯卡通| 亚洲欧美中文字幕日韩二区| 99热这里只有是精品在线观看| 亚洲成人一二三区av| 老熟女久久久| 涩涩av久久男人的天堂| 成年女人在线观看亚洲视频| 日韩成人伦理影院| h日本视频在线播放| 两个人的视频大全免费| 2018国产大陆天天弄谢| 免费在线观看成人毛片| 高清视频免费观看一区二区| 国产精品秋霞免费鲁丝片| 女人十人毛片免费观看3o分钟| 国产亚洲一区二区精品| av不卡在线播放| 国产亚洲av片在线观看秒播厂| 免费av不卡在线播放| 久久精品国产自在天天线| 在线观看美女被高潮喷水网站| 亚洲精品第二区| 久久毛片免费看一区二区三区| 人妻制服诱惑在线中文字幕| 九九在线视频观看精品| av.在线天堂| 成人毛片60女人毛片免费| 观看av在线不卡| 大片电影免费在线观看免费| 一级毛片电影观看| 视频中文字幕在线观看| 好男人视频免费观看在线| 久久久久网色| 韩国av在线不卡| 丝袜喷水一区| 国产成人freesex在线| 欧美日韩国产mv在线观看视频 | 亚洲美女搞黄在线观看| 亚洲精品一区蜜桃| av专区在线播放| 国产精品99久久久久久久久| 亚洲久久久国产精品| 精品国产露脸久久av麻豆| 日韩av在线免费看完整版不卡| 国产真实伦视频高清在线观看| 久久久久精品久久久久真实原创| 日韩中文字幕视频在线看片 | 成人免费观看视频高清| 五月玫瑰六月丁香| 在线免费观看不下载黄p国产| 国产av一区二区精品久久 | 啦啦啦中文免费视频观看日本| 成年免费大片在线观看| 久久av网站| 精华霜和精华液先用哪个| 嫩草影院新地址| 高清毛片免费看| 人体艺术视频欧美日本| 丰满乱子伦码专区| 极品教师在线视频| 亚洲欧美一区二区三区黑人 | 大香蕉97超碰在线| 色5月婷婷丁香| 欧美日韩精品成人综合77777| 久久99热这里只频精品6学生| 亚洲精品456在线播放app| 大香蕉久久网| 97超视频在线观看视频| 日本欧美国产在线视频| 内射极品少妇av片p| 亚洲精品日韩av片在线观看| 久久国内精品自在自线图片| 人妻一区二区av| 26uuu在线亚洲综合色| 国模一区二区三区四区视频| 青春草视频在线免费观看| 这个男人来自地球电影免费观看 | 欧美精品国产亚洲| 亚洲精品456在线播放app| 黄色日韩在线| 又粗又硬又长又爽又黄的视频| 久久久久网色| 欧美激情极品国产一区二区三区 | 一级爰片在线观看| 精品人妻视频免费看| 国产成人免费无遮挡视频| 国产亚洲av片在线观看秒播厂| 成年美女黄网站色视频大全免费 | 高清日韩中文字幕在线| 少妇人妻精品综合一区二区| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| 秋霞伦理黄片| 一区二区av电影网| 国产中年淑女户外野战色| 国产毛片在线视频| 一级毛片 在线播放| 亚洲国产色片| 国产精品国产av在线观看| 亚洲欧美日韩无卡精品| 欧美成人一区二区免费高清观看| 日韩一区二区视频免费看| 国产精品偷伦视频观看了| 国产免费视频播放在线视频| 久久婷婷青草| 少妇人妻 视频| 老司机影院成人| 少妇的逼水好多| 深夜a级毛片| 精品久久久精品久久久| 日日摸夜夜添夜夜添av毛片| 亚洲欧洲国产日韩| 亚洲人与动物交配视频| 国产精品99久久久久久久久| 视频区图区小说| 久久久久人妻精品一区果冻| 亚洲美女搞黄在线观看| 欧美区成人在线视频| 免费观看a级毛片全部| 午夜老司机福利剧场| 免费少妇av软件| 国产精品成人在线| 国产精品久久久久成人av| 少妇的逼水好多| 亚洲精品日韩av片在线观看| 国产美女午夜福利| 国产精品.久久久| 中文精品一卡2卡3卡4更新| 久久99热6这里只有精品| 三级经典国产精品| 中国国产av一级| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 一区二区三区四区激情视频| 亚洲三级黄色毛片| 一区二区av电影网| 啦啦啦啦在线视频资源| 午夜福利网站1000一区二区三区| 嫩草影院新地址| 一级爰片在线观看| 少妇人妻 视频| 国产精品一区www在线观看| 18禁裸乳无遮挡动漫免费视频| 中文乱码字字幕精品一区二区三区| 欧美日本视频| 精品午夜福利在线看| 男人狂女人下面高潮的视频| 一级毛片黄色毛片免费观看视频| 久久综合国产亚洲精品| 欧美97在线视频| 亚洲天堂av无毛| 一本一本综合久久| 国产在线免费精品| 18禁在线播放成人免费| 老熟女久久久| 这个男人来自地球电影免费观看 | 男女国产视频网站| 亚洲国产毛片av蜜桃av| 久久精品国产自在天天线| 99久国产av精品国产电影| 高清视频免费观看一区二区| 亚洲国产欧美人成| 午夜精品国产一区二区电影| 久久热精品热| 欧美成人a在线观看| 黄色欧美视频在线观看| 搡女人真爽免费视频火全软件| 国产精品嫩草影院av在线观看| 18+在线观看网站| 国产精品av视频在线免费观看| 亚洲国产色片| 夫妻性生交免费视频一级片| 日韩成人av中文字幕在线观看| 97在线人人人人妻| 黄片无遮挡物在线观看| 五月天丁香电影| 欧美激情极品国产一区二区三区 | 日韩成人伦理影院| 夫妻午夜视频| 国产免费又黄又爽又色| www.色视频.com| 精品人妻熟女av久视频| 亚洲欧美精品自产自拍| 国产成人a∨麻豆精品| 久热这里只有精品99| 亚洲国产av新网站| 嘟嘟电影网在线观看| 午夜福利网站1000一区二区三区| 美女内射精品一级片tv| 麻豆成人av视频| 大码成人一级视频| videos熟女内射| 亚洲欧美日韩无卡精品| 国产乱来视频区| 制服丝袜香蕉在线| 亚洲综合精品二区| 国产男女超爽视频在线观看| 国产av国产精品国产| 日韩欧美一区视频在线观看 | 91在线精品国自产拍蜜月| 成年女人在线观看亚洲视频| 亚洲精品色激情综合| 国产精品久久久久久精品古装| 国产精品久久久久久久电影| 毛片一级片免费看久久久久| 黄色配什么色好看| 天堂8中文在线网| 男人添女人高潮全过程视频| 亚州av有码| 国产91av在线免费观看| 国产午夜精品久久久久久一区二区三区| 亚洲精品国产av成人精品| 不卡视频在线观看欧美| 国产精品女同一区二区软件| 国产在线男女| 十分钟在线观看高清视频www | 亚洲精品国产av蜜桃| 国产乱人视频| 春色校园在线视频观看| 人体艺术视频欧美日本| 一区二区三区免费毛片| 亚洲欧美日韩东京热| 一边亲一边摸免费视频| 亚洲国产毛片av蜜桃av| 我的女老师完整版在线观看| 国产乱人偷精品视频| 免费观看a级毛片全部| 少妇人妻精品综合一区二区| 久久热精品热| 国产精品久久久久久精品电影小说 | av播播在线观看一区| 日韩亚洲欧美综合| 欧美最新免费一区二区三区| 观看美女的网站| 欧美97在线视频| 美女国产视频在线观看| 在线亚洲精品国产二区图片欧美 | 亚洲第一区二区三区不卡| 一区二区三区四区激情视频| 免费看光身美女| 亚洲欧洲日产国产| 久久这里有精品视频免费| 爱豆传媒免费全集在线观看| 国产色爽女视频免费观看| 观看av在线不卡| 国产精品国产三级国产专区5o| 成年美女黄网站色视频大全免费 | 一级二级三级毛片免费看| 久久婷婷青草| 麻豆国产97在线/欧美| 久久精品国产a三级三级三级| 在线观看免费日韩欧美大片 | 美女福利国产在线 | xxx大片免费视频| 国产深夜福利视频在线观看| 久久久久精品久久久久真实原创| 亚洲精品456在线播放app| 亚洲精品久久久久久婷婷小说| 2021少妇久久久久久久久久久| 精品酒店卫生间| 国产精品女同一区二区软件| 麻豆成人av视频| 成人影院久久| 中文资源天堂在线| 国国产精品蜜臀av免费| 欧美97在线视频| 免费不卡的大黄色大毛片视频在线观看| 18禁动态无遮挡网站| 黄片无遮挡物在线观看| 亚洲成人一二三区av| 亚洲精品久久久久久婷婷小说| 热99国产精品久久久久久7| 国产成人精品福利久久| 国产免费视频播放在线视频| 伦理电影大哥的女人| 黄色怎么调成土黄色| 免费看不卡的av| 黑人猛操日本美女一级片| 男人添女人高潮全过程视频| 免费久久久久久久精品成人欧美视频 | 成人免费观看视频高清| 91精品伊人久久大香线蕉| 久久精品国产亚洲av涩爱| 亚洲av电影在线观看一区二区三区| 国产成人a区在线观看| 少妇人妻 视频| 纯流量卡能插随身wifi吗| 一区在线观看完整版| 国产无遮挡羞羞视频在线观看| 日韩人妻高清精品专区| 国产 精品1| 边亲边吃奶的免费视频| 精品人妻熟女av久视频| 国产永久视频网站| 狂野欧美激情性bbbbbb| 天天躁夜夜躁狠狠久久av| 简卡轻食公司| 天天躁夜夜躁狠狠久久av| 丝袜脚勾引网站| 人妻夜夜爽99麻豆av| 国产午夜精品久久久久久一区二区三区| 久久久久精品性色| 日韩av不卡免费在线播放| 国产 精品1| 国产爱豆传媒在线观看| a 毛片基地| 亚洲精品aⅴ在线观看| 男人和女人高潮做爰伦理| 伦理电影大哥的女人| 国产免费视频播放在线视频| 国产男女超爽视频在线观看| 婷婷色综合www| 国产熟女欧美一区二区| 一本一本综合久久| 永久网站在线| 丝袜脚勾引网站| 久久久久久九九精品二区国产| 国产精品熟女久久久久浪| 午夜福利在线观看免费完整高清在| 亚洲高清免费不卡视频| 免费av不卡在线播放| 色视频在线一区二区三区| 色5月婷婷丁香| 一级毛片久久久久久久久女| 日本黄色日本黄色录像| 一二三四中文在线观看免费高清| 成人免费观看视频高清| 精品久久久精品久久久| av在线蜜桃| 国内精品宾馆在线| 两个人的视频大全免费| 丝袜脚勾引网站| 日韩av不卡免费在线播放| 免费黄色在线免费观看| 亚洲欧美精品专区久久| 99热全是精品| 哪个播放器可以免费观看大片| 美女高潮的动态| 亚洲av中文av极速乱| 日韩一区二区三区影片| 成人黄色视频免费在线看| 干丝袜人妻中文字幕| 亚洲天堂av无毛| 纯流量卡能插随身wifi吗| 亚洲美女视频黄频| 伊人久久国产一区二区| 中文精品一卡2卡3卡4更新| 精品视频人人做人人爽| 狂野欧美激情性xxxx在线观看| 人妻 亚洲 视频| 亚洲无线观看免费| 九色成人免费人妻av| 日韩视频在线欧美| 国产精品爽爽va在线观看网站| 国产午夜精品一二区理论片| 欧美性感艳星| 久久久久久久久久成人| 嫩草影院新地址| 不卡视频在线观看欧美| 夫妻性生交免费视频一级片| 一区二区三区乱码不卡18| 国产在线免费精品| 国产视频首页在线观看| 亚洲第一av免费看| 又爽又黄a免费视频| 插逼视频在线观看| av一本久久久久| 久久久久久久亚洲中文字幕| 欧美日韩精品成人综合77777| 丰满乱子伦码专区| 国产精品伦人一区二区| 人妻夜夜爽99麻豆av| 国产高清国产精品国产三级 | 内射极品少妇av片p| 精品一区二区免费观看| 香蕉精品网在线| 久久久午夜欧美精品| av不卡在线播放| 久久女婷五月综合色啪小说| 国产综合精华液| 国产伦精品一区二区三区四那| 菩萨蛮人人尽说江南好唐韦庄| 99热6这里只有精品| av.在线天堂| 18禁在线无遮挡免费观看视频| 欧美日本视频| 国产精品三级大全| 伊人久久国产一区二区| 看非洲黑人一级黄片| 日日撸夜夜添| 国产精品久久久久久久久免| 中文字幕av成人在线电影| 国产精品久久久久久av不卡| 人人妻人人爽人人添夜夜欢视频 | 精品久久久久久久末码| 国产av码专区亚洲av| 国产片特级美女逼逼视频| 久久久久久九九精品二区国产| 亚洲人与动物交配视频| 欧美精品一区二区免费开放| 亚洲美女搞黄在线观看| 久久人人爽人人爽人人片va| 91狼人影院| 国产精品久久久久久av不卡| av在线蜜桃| 欧美xxxx性猛交bbbb| 91aial.com中文字幕在线观看| 韩国高清视频一区二区三区| 日本wwww免费看| 久久毛片免费看一区二区三区| 欧美另类一区| 人人妻人人看人人澡| 亚洲精品乱久久久久久| 久久久久久人妻| 日本色播在线视频| 少妇 在线观看| 亚洲国产高清在线一区二区三| 国产精品一区二区在线不卡| 亚洲人成网站高清观看| 26uuu在线亚洲综合色| 欧美 日韩 精品 国产| 国产精品成人在线| 亚洲精品成人av观看孕妇| 日韩欧美 国产精品| 女人十人毛片免费观看3o分钟| 国产精品人妻久久久久久| a级毛片免费高清观看在线播放| tube8黄色片| 国产精品一及| 成年av动漫网址| 婷婷色综合大香蕉| 美女视频免费永久观看网站| 只有这里有精品99| 99热6这里只有精品| 永久网站在线| av不卡在线播放| 国产成人精品久久久久久| 一区二区三区乱码不卡18| 99精国产麻豆久久婷婷| 97热精品久久久久久| 免费av不卡在线播放| 国产免费视频播放在线视频| 人妻夜夜爽99麻豆av| 一级黄片播放器| 欧美日韩视频精品一区| 精品亚洲成国产av| 中国国产av一级| av福利片在线观看| 成年人午夜在线观看视频| 免费高清在线观看视频在线观看| 男人狂女人下面高潮的视频| 18禁裸乳无遮挡动漫免费视频| 日韩一区二区视频免费看| 久久精品久久精品一区二区三区| 插逼视频在线观看| a级毛色黄片| 少妇人妻久久综合中文| 久久久欧美国产精品| 少妇猛男粗大的猛烈进出视频| 欧美三级亚洲精品| 久久久久国产精品人妻一区二区| 亚洲成人中文字幕在线播放| 老女人水多毛片| 搡老乐熟女国产| .国产精品久久| 亚洲人与动物交配视频| 国产午夜精品久久久久久一区二区三区| 内地一区二区视频在线| av又黄又爽大尺度在线免费看| 国产午夜精品一二区理论片| 18禁在线无遮挡免费观看视频| 午夜福利网站1000一区二区三区| 妹子高潮喷水视频| av又黄又爽大尺度在线免费看| kizo精华| 午夜老司机福利剧场| 日韩一本色道免费dvd| 熟妇人妻不卡中文字幕| 亚洲欧洲日产国产| h日本视频在线播放| 99久国产av精品国产电影| 国产真实伦视频高清在线观看| 极品少妇高潮喷水抽搐| 狂野欧美激情性xxxx在线观看| 大香蕉久久网| 国产黄色视频一区二区在线观看| 最黄视频免费看| 一级毛片aaaaaa免费看小| 一区二区三区四区激情视频| 男的添女的下面高潮视频| 欧美日韩亚洲高清精品| 九色成人免费人妻av| 国产伦在线观看视频一区| 国产精品久久久久久av不卡| 22中文网久久字幕| 直男gayav资源| 亚洲精品国产成人久久av| 国产成人精品福利久久| 最近2019中文字幕mv第一页| 精品一区二区免费观看| 亚洲精品一二三| 丰满少妇做爰视频| 亚洲欧美精品专区久久| 夜夜爽夜夜爽视频| 一二三四中文在线观看免费高清| 纵有疾风起免费观看全集完整版| 国产精品伦人一区二区| 少妇的逼好多水| 欧美激情国产日韩精品一区| 日韩亚洲欧美综合| 免费观看a级毛片全部| 日韩电影二区| 91精品伊人久久大香线蕉| 亚洲精品视频女| 日日啪夜夜爽| 久久人人爽av亚洲精品天堂 | 最近的中文字幕免费完整| 亚洲伊人久久精品综合| 91狼人影院| 女性被躁到高潮视频| 久久精品久久久久久噜噜老黄| 成人综合一区亚洲| 女人十人毛片免费观看3o分钟| 人妻夜夜爽99麻豆av| 国产精品久久久久久久久免| 精品久久久久久久久av| av福利片在线观看| 99热全是精品| 亚州av有码| 黄色一级大片看看| 99热6这里只有精品| 色婷婷av一区二区三区视频| 麻豆成人av视频| 亚洲国产色片| 精品国产乱码久久久久久小说| 精品视频人人做人人爽| 精品久久久久久久末码| 亚洲精品久久久久久婷婷小说| 高清在线视频一区二区三区| 欧美日韩视频高清一区二区三区二| 欧美日韩综合久久久久久| 亚洲成人一二三区av| 精品久久久久久久久av| 欧美日韩综合久久久久久| 免费看日本二区| 一级片'在线观看视频| 波野结衣二区三区在线| 妹子高潮喷水视频| 国产精品三级大全| .国产精品久久| 国产无遮挡羞羞视频在线观看| av又黄又爽大尺度在线免费看| 国产成人freesex在线| 亚洲欧美一区二区三区国产| 成人国产av品久久久| 精品一区在线观看国产| 国产av码专区亚洲av| 晚上一个人看的免费电影| 一级毛片电影观看| 日本免费在线观看一区| 亚洲精品国产av蜜桃| 久久久久人妻精品一区果冻| 男女免费视频国产| 内地一区二区视频在线|