• <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
    www日本在线高清视频| 亚洲精品中文字幕一二三四区| 久久久国产成人免费| www.熟女人妻精品国产| 国产精品99久久久久久久久| 久久久久亚洲av毛片大全| 五月玫瑰六月丁香| 免费看十八禁软件| 久久久久国产一级毛片高清牌| 少妇的逼水好多| 精品一区二区三区视频在线观看免费| 热99在线观看视频| 国产精品 国内视频| 精品无人区乱码1区二区| 久久久精品大字幕| 免费在线观看影片大全网站| 国产欧美日韩精品亚洲av| 免费看十八禁软件| 国产又黄又爽又无遮挡在线| av福利片在线观看| 一进一出好大好爽视频| 亚洲精品久久国产高清桃花| 丰满人妻熟妇乱又伦精品不卡| 国产毛片a区久久久久| 亚洲美女黄片视频| 国产伦人伦偷精品视频| 免费看日本二区| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播放欧美日韩| 国产私拍福利视频在线观看| 国产成人影院久久av| 色综合欧美亚洲国产小说| 亚洲av片天天在线观看| 国产av在哪里看| 国产三级黄色录像| 中文字幕精品亚洲无线码一区| 91老司机精品| 999精品在线视频| 午夜免费观看网址| 欧美乱妇无乱码| 一本精品99久久精品77| 亚洲熟女毛片儿| 久久99热这里只有精品18| 国产成人精品久久二区二区91| 欧美丝袜亚洲另类 | 亚洲成人精品中文字幕电影| 免费av毛片视频| 美女高潮喷水抽搐中文字幕| 嫩草影院入口| www.精华液| 又黄又爽又免费观看的视频| 久久精品国产亚洲av香蕉五月| 91在线精品国自产拍蜜月 | 又黄又粗又硬又大视频| 日韩 欧美 亚洲 中文字幕| 日本免费a在线| 亚洲av成人av| tocl精华| 精品久久久久久久人妻蜜臀av| 亚洲欧美日韩东京热| 国产免费男女视频| 成人无遮挡网站| 欧美激情在线99| 国产欧美日韩一区二区精品| 色综合欧美亚洲国产小说| 国产v大片淫在线免费观看| 亚洲人成伊人成综合网2020| 狂野欧美激情性xxxx| 亚洲精品乱码久久久v下载方式 | 亚洲成a人片在线一区二区| 国产精品久久电影中文字幕| 精品国产乱子伦一区二区三区| 亚洲国产欧美网| 欧美精品啪啪一区二区三区| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久亚洲av鲁大| 少妇丰满av| 五月玫瑰六月丁香| 国产精品自产拍在线观看55亚洲| 中国美女看黄片| 久久伊人香网站| 欧美日韩精品网址| 深夜精品福利| 18禁黄网站禁片午夜丰满| 国产1区2区3区精品| 一区二区三区高清视频在线| 国产一区二区三区在线臀色熟女| 老司机午夜福利在线观看视频| 国产精品乱码一区二三区的特点| 亚洲激情在线av| 日本一二三区视频观看| 天堂av国产一区二区熟女人妻| 18禁观看日本| 免费观看人在逋| 人妻久久中文字幕网| 岛国在线观看网站| 在线免费观看的www视频| 国产亚洲av嫩草精品影院| 99热这里只有是精品50| 婷婷精品国产亚洲av在线| 黄色视频,在线免费观看| 欧美丝袜亚洲另类 | 精品久久久久久久末码| 欧美最黄视频在线播放免费| 国产野战对白在线观看| 国产精品久久视频播放| 亚洲 欧美一区二区三区| 嫩草影视91久久| 在线观看美女被高潮喷水网站 | а√天堂www在线а√下载| av女优亚洲男人天堂 | 人人妻,人人澡人人爽秒播| 午夜福利免费观看在线| 法律面前人人平等表现在哪些方面| 国产精品99久久99久久久不卡| 色综合亚洲欧美另类图片| 99riav亚洲国产免费| 日本一本二区三区精品| av中文乱码字幕在线| 免费看美女性在线毛片视频| 黄色丝袜av网址大全| 亚洲成人免费电影在线观看| 精品久久久久久久毛片微露脸| 热99re8久久精品国产| 熟女少妇亚洲综合色aaa.| 桃色一区二区三区在线观看| 久久99热这里只有精品18| 日韩欧美免费精品| 老司机午夜十八禁免费视频| 亚洲av中文字字幕乱码综合| 亚洲av五月六月丁香网| 午夜影院日韩av| 久久久国产精品麻豆| 国产麻豆成人av免费视频| 精品国产乱子伦一区二区三区| 亚洲国产高清在线一区二区三| 1024香蕉在线观看| 久久久久国内视频| 国产又黄又爽又无遮挡在线| 又紧又爽又黄一区二区| 国产欧美日韩一区二区三| 日日干狠狠操夜夜爽| 免费看日本二区| 特级一级黄色大片| 99久久久亚洲精品蜜臀av| 国产欧美日韩精品亚洲av| 精品熟女少妇八av免费久了| 日韩大尺度精品在线看网址| 亚洲午夜精品一区,二区,三区| 亚洲美女黄片视频| 亚洲av熟女| 国产精品 欧美亚洲| 精品福利观看| 久久精品国产亚洲av香蕉五月| 亚洲av第一区精品v没综合| 成年人黄色毛片网站| 亚洲国产精品久久男人天堂| 男女那种视频在线观看| 99riav亚洲国产免费| 国产成人aa在线观看| 国产高清三级在线| 国产 一区 欧美 日韩| 窝窝影院91人妻| 国产一级毛片七仙女欲春2| 国产精品免费一区二区三区在线| 亚洲欧美日韩卡通动漫| 日本黄色片子视频| 一个人看视频在线观看www免费 | 99国产精品一区二区蜜桃av| 两人在一起打扑克的视频| 男女那种视频在线观看| 9191精品国产免费久久| 在线观看免费视频日本深夜| 美女午夜性视频免费| 亚洲欧美精品综合一区二区三区| 午夜福利高清视频| 国产成年人精品一区二区| 欧美色视频一区免费| 脱女人内裤的视频| 99国产精品一区二区三区| 国产一级毛片七仙女欲春2| 亚洲自拍偷在线| 欧美在线黄色| 久久中文字幕人妻熟女| 国产淫片久久久久久久久 | 国内揄拍国产精品人妻在线| 日本a在线网址| 97超级碰碰碰精品色视频在线观看| 91九色精品人成在线观看| 亚洲中文字幕一区二区三区有码在线看 | 黄色 视频免费看| 国产精品久久久人人做人人爽| 亚洲性夜色夜夜综合| 国产伦一二天堂av在线观看| 日本成人三级电影网站| 国产成人系列免费观看| 偷拍熟女少妇极品色| av女优亚洲男人天堂 | 搡老岳熟女国产| 51午夜福利影视在线观看| 精品无人区乱码1区二区| 精品99又大又爽又粗少妇毛片 | 又大又爽又粗| 亚洲七黄色美女视频| 99国产极品粉嫩在线观看| 国产高潮美女av| 国产精品久久久久久精品电影| 日韩中文字幕欧美一区二区| 国语自产精品视频在线第100页| 亚洲国产精品sss在线观看| 久久精品影院6| 免费搜索国产男女视频| 久久人妻av系列| 黄色 视频免费看| 国产精品一及| 高清毛片免费观看视频网站| 热99re8久久精品国产| 精品久久久久久久久久久久久| 在线播放国产精品三级| 亚洲aⅴ乱码一区二区在线播放| 在线观看日韩欧美| 精品99又大又爽又粗少妇毛片 | 欧美乱码精品一区二区三区| 日本与韩国留学比较| 免费观看人在逋| 精品午夜福利视频在线观看一区| 欧美精品啪啪一区二区三区| 国产一区在线观看成人免费| 狠狠狠狠99中文字幕| 蜜桃久久精品国产亚洲av| 国语自产精品视频在线第100页| 久久精品国产亚洲av香蕉五月| 色精品久久人妻99蜜桃| 天堂网av新在线| 欧美一区二区精品小视频在线| 国产成年人精品一区二区| 一本一本综合久久| 啦啦啦韩国在线观看视频| 久久精品人妻少妇| 制服丝袜大香蕉在线| 精品久久久久久成人av| 又紧又爽又黄一区二区| 天堂动漫精品| 99久久国产精品久久久| 久久久久国产精品人妻aⅴ院| 97人妻精品一区二区三区麻豆| 国产精品 国内视频| 亚洲天堂国产精品一区在线| 国产一区二区三区视频了| 久久精品综合一区二区三区| 亚洲人成伊人成综合网2020| 99久久99久久久精品蜜桃| 午夜福利在线在线| 免费看美女性在线毛片视频| 久久草成人影院| 激情在线观看视频在线高清| 久久这里只有精品中国| 国产成人影院久久av| 午夜精品一区二区三区免费看| 18禁观看日本| 亚洲性夜色夜夜综合| 久久久久国产一级毛片高清牌| h日本视频在线播放| 久久婷婷人人爽人人干人人爱| 青草久久国产| 欧美最黄视频在线播放免费| av中文乱码字幕在线| 亚洲av熟女| 久久久久久久久免费视频了| 亚洲中文字幕一区二区三区有码在线看 | 精品国产乱码久久久久久男人| 久久草成人影院| 黄色日韩在线| 99热只有精品国产| 女警被强在线播放| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 老汉色∧v一级毛片| 久久久久久久午夜电影| 国产69精品久久久久777片 | 亚洲精品中文字幕一二三四区| 国产激情久久老熟女| 亚洲电影在线观看av| 国产三级黄色录像| 久久草成人影院| 小说图片视频综合网站| 日韩三级视频一区二区三区| 桃色一区二区三区在线观看| 国产乱人伦免费视频| 变态另类丝袜制服| 日韩av在线大香蕉| 国产久久久一区二区三区| 在线观看免费视频日本深夜| 1024香蕉在线观看| 一进一出抽搐动态| 亚洲中文日韩欧美视频| 精品一区二区三区视频在线观看免费| 国产免费av片在线观看野外av| 欧美乱妇无乱码| 亚洲一区二区三区不卡视频| 婷婷六月久久综合丁香| www.www免费av| 在线国产一区二区在线| 午夜a级毛片| 国产亚洲av嫩草精品影院| 国产激情欧美一区二区| 国产成人系列免费观看| 18禁美女被吸乳视频| 婷婷丁香在线五月| 九九在线视频观看精品| 亚洲性夜色夜夜综合| 欧美一级毛片孕妇| 在线观看66精品国产| 麻豆国产97在线/欧美| 欧美3d第一页| 在线观看66精品国产| 国产精品九九99| 国产单亲对白刺激| 夜夜看夜夜爽夜夜摸| 国产精品精品国产色婷婷| 国产精品 欧美亚洲| ponron亚洲| 51午夜福利影视在线观看| 亚洲人成网站高清观看| 婷婷精品国产亚洲av在线| 欧美黑人巨大hd| 国产精品野战在线观看| 亚洲片人在线观看| 在线观看免费午夜福利视频| 一个人看的www免费观看视频| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 欧美精品啪啪一区二区三区| 久久精品国产清高在天天线| 中文字幕精品亚洲无线码一区| 国产精品亚洲美女久久久| 亚洲国产看品久久| 久久午夜亚洲精品久久| 国产精品亚洲美女久久久| 国产91精品成人一区二区三区| 97超级碰碰碰精品色视频在线观看| 亚洲专区国产一区二区| 国产亚洲欧美在线一区二区| 精品日产1卡2卡| 成人特级黄色片久久久久久久| av在线天堂中文字幕| www.精华液| ponron亚洲| 亚洲国产精品成人综合色| 99国产精品99久久久久| 男女做爰动态图高潮gif福利片| 可以在线观看的亚洲视频| 白带黄色成豆腐渣| 国产午夜精品久久久久久| 色综合欧美亚洲国产小说| av在线蜜桃| 亚洲欧美日韩卡通动漫| 日日夜夜操网爽| av欧美777| 99久久国产精品久久久| 黄色日韩在线| 99国产综合亚洲精品| 一本精品99久久精品77| 午夜免费成人在线视频| 欧美成狂野欧美在线观看| 亚洲国产精品合色在线| 国产精品av视频在线免费观看| 国产亚洲精品久久久com| 在线十欧美十亚洲十日本专区| 国产黄a三级三级三级人| 婷婷丁香在线五月| 亚洲成人免费电影在线观看| 精品国产乱码久久久久久男人| 中文字幕久久专区| 黄色女人牲交| 国产视频内射| 999精品在线视频| 搡老妇女老女人老熟妇| 国内精品久久久久精免费| 亚洲片人在线观看| 久久久久性生活片| 国产亚洲精品久久久久久毛片| 国产精品av久久久久免费| 日韩欧美国产在线观看| 日韩人妻高清精品专区| 久久亚洲真实| 国产av麻豆久久久久久久| 久久这里只有精品中国| 日本撒尿小便嘘嘘汇集6| 免费av不卡在线播放| 亚洲性夜色夜夜综合| 夜夜躁狠狠躁天天躁| 村上凉子中文字幕在线| 亚洲欧美日韩高清专用| 免费无遮挡裸体视频| 综合色av麻豆| xxx96com| 一区福利在线观看| 精品一区二区三区四区五区乱码| 亚洲av五月六月丁香网| 久久久国产精品麻豆| 精品久久久久久久久久久久久| 久久久国产欧美日韩av| 男插女下体视频免费在线播放| 亚洲欧美日韩东京热| 欧美一区二区国产精品久久精品| av黄色大香蕉| 国产男靠女视频免费网站| 高清毛片免费观看视频网站| 国产成+人综合+亚洲专区| 国产av不卡久久| 国产av一区在线观看免费| 精品久久久久久,| 成人av在线播放网站| 欧美+亚洲+日韩+国产| 桃色一区二区三区在线观看| 亚洲成人中文字幕在线播放| 欧美国产日韩亚洲一区| 精品久久久久久,| 亚洲人成电影免费在线| 黑人巨大精品欧美一区二区mp4| 久久中文字幕一级| 国产精品野战在线观看| 欧美午夜高清在线| 亚洲av电影不卡..在线观看| 国产高潮美女av| 欧美乱妇无乱码| 黄色成人免费大全| 狠狠狠狠99中文字幕| 黑人操中国人逼视频| 国产欧美日韩精品亚洲av| 欧美国产日韩亚洲一区| 18禁黄网站禁片免费观看直播| 国产日本99.免费观看| 成人国产综合亚洲| 成人精品一区二区免费| 村上凉子中文字幕在线| 国产成人aa在线观看| 一级黄色大片毛片| 99国产综合亚洲精品| 国产精品一区二区免费欧美| 国产蜜桃级精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人不卡在线观看播放网| 99在线视频只有这里精品首页| 欧美日韩精品网址| 狠狠狠狠99中文字幕| 国产亚洲精品av在线| 午夜福利在线在线| 国产精品一及| 香蕉国产在线看| 国产精品99久久99久久久不卡| 免费看美女性在线毛片视频| 久久久久久久久中文| 五月玫瑰六月丁香| 国内少妇人妻偷人精品xxx网站 | 午夜福利在线在线| 亚洲精品一区av在线观看| 亚洲av第一区精品v没综合| 黑人欧美特级aaaaaa片| 国产精品香港三级国产av潘金莲| 在线免费观看不下载黄p国产 | 国产精品,欧美在线| 一级作爱视频免费观看| 亚洲专区国产一区二区| 91久久精品国产一区二区成人 | 99在线视频只有这里精品首页| av天堂中文字幕网| 99国产精品一区二区蜜桃av| 欧美乱色亚洲激情| 男人舔女人下体高潮全视频| 成人鲁丝片一二三区免费| 国产精品久久久久久亚洲av鲁大| 18禁黄网站禁片午夜丰满| 我的老师免费观看完整版| 欧美xxxx黑人xx丫x性爽| av福利片在线观看| 色精品久久人妻99蜜桃| 天天一区二区日本电影三级| 久久香蕉国产精品| 淫妇啪啪啪对白视频| 国产欧美日韩精品亚洲av| 熟女少妇亚洲综合色aaa.| 国产一区在线观看成人免费| 精品日产1卡2卡| avwww免费| 欧美国产日韩亚洲一区| 久久这里只有精品中国| 亚洲男人的天堂狠狠| 精品久久久久久,| 好男人在线观看高清免费视频| 亚洲 欧美 日韩 在线 免费| 久久久久国内视频| 色老头精品视频在线观看| 长腿黑丝高跟| 国产97色在线日韩免费| 18禁观看日本| 成人性生交大片免费视频hd| 999久久久精品免费观看国产| 午夜免费激情av| 成人特级黄色片久久久久久久| 一本久久中文字幕| 夜夜看夜夜爽夜夜摸| 国产亚洲精品久久久com| 欧美日韩黄片免| 精品熟女少妇八av免费久了| 国产乱人视频| 欧美日韩亚洲国产一区二区在线观看| 国产欧美日韩精品亚洲av| 熟妇人妻久久中文字幕3abv| 免费观看的影片在线观看| 婷婷精品国产亚洲av在线| 丰满人妻熟妇乱又伦精品不卡| 国产精品日韩av在线免费观看| 国产精品国产高清国产av| 久9热在线精品视频| 一级a爱片免费观看的视频| 99热这里只有是精品50| 一个人免费在线观看电影 | 国产美女午夜福利| 久久这里只有精品中国| 小蜜桃在线观看免费完整版高清| 99久久无色码亚洲精品果冻| 国产乱人视频| 国内久久婷婷六月综合欲色啪| 视频区欧美日本亚洲| 国产精品98久久久久久宅男小说| tocl精华| 亚洲 国产 在线| 很黄的视频免费| 美女免费视频网站| 欧美一级毛片孕妇| 国产主播在线观看一区二区| 一个人观看的视频www高清免费观看 | 激情在线观看视频在线高清| 一级作爱视频免费观看| av欧美777| 麻豆一二三区av精品| a级毛片在线看网站| 日韩精品青青久久久久久| 在线播放国产精品三级| 欧美性猛交╳xxx乱大交人| 中文亚洲av片在线观看爽| 国产亚洲欧美98| 亚洲av美国av| 日韩国内少妇激情av| 一本一本综合久久| 国产一区二区在线观看日韩 | 亚洲自偷自拍图片 自拍| 久久九九热精品免费| 搡老岳熟女国产| 国产精品电影一区二区三区| 老鸭窝网址在线观看| 亚洲 国产 在线| 国产人伦9x9x在线观看| 免费看美女性在线毛片视频| 欧美日韩黄片免| 成人三级做爰电影| 国产成人福利小说| 国产精品1区2区在线观看.| x7x7x7水蜜桃| 精品乱码久久久久久99久播| 免费看光身美女| 精品熟女少妇八av免费久了| 久久99热这里只有精品18| 国产亚洲精品av在线| 国产精华一区二区三区| 18禁黄网站禁片免费观看直播| xxx96com| 国产精品九九99| 国产亚洲精品一区二区www| 免费搜索国产男女视频| 99久久久亚洲精品蜜臀av| 亚洲精品乱码久久久v下载方式 | 国产97色在线日韩免费| www.999成人在线观看| 岛国在线观看网站| 美女 人体艺术 gogo| 亚洲,欧美精品.| 757午夜福利合集在线观看| 亚洲人与动物交配视频| 1024香蕉在线观看| 不卡av一区二区三区| 日韩有码中文字幕| 91麻豆精品激情在线观看国产| 搞女人的毛片| 国产美女午夜福利| 国产精品免费一区二区三区在线| 天天一区二区日本电影三级| 久久性视频一级片| 嫩草影视91久久| 欧美日韩亚洲国产一区二区在线观看| 长腿黑丝高跟| 三级男女做爰猛烈吃奶摸视频| av片东京热男人的天堂| 长腿黑丝高跟| 三级男女做爰猛烈吃奶摸视频| 久久久国产成人精品二区| 变态另类丝袜制服| 在线观看免费视频日本深夜| 中出人妻视频一区二区| 在线免费观看的www视频| 亚洲成av人片在线播放无| 日韩av在线大香蕉| 国产99白浆流出| 亚洲精品456在线播放app | 精品福利观看| 日本黄色片子视频| 亚洲精品粉嫩美女一区| 久久香蕉精品热| 欧美另类亚洲清纯唯美| 一个人免费在线观看的高清视频| 日本黄大片高清|