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

    基于HHT能量和最大Lyapunov指數(shù)的蛇行分類方法

    2023-08-30 03:35:38李艷萍陳春俊
    計(jì)算機(jī)測量與控制 2023年8期
    關(guān)鍵詞:蛇行踏面構(gòu)架

    王 敏,寧 靜,趙 飛,李艷萍,陳春俊

    (西南交通大學(xué) 機(jī)械工程學(xué)院,成都 610031)

    0 引言

    蛇行運(yùn)動(dòng)是車輛動(dòng)力學(xué)系統(tǒng)的核心問題之一,車輛在運(yùn)行過程中發(fā)生收斂較慢的小幅蛇行或者劇烈的大幅蛇行都會(huì)嚴(yán)重地影響車輛運(yùn)行安全,因此對車輛小幅蛇行和大幅蛇行的在線監(jiān)測都至關(guān)重要。學(xué)者們通過大量的理論研究總結(jié)出輪軌參數(shù)和懸掛系統(tǒng)參數(shù)對車輛系統(tǒng)蛇行運(yùn)動(dòng)的影響規(guī)律[1-5]。在這些理論研究的成果上,更多的學(xué)者致力于探索如何將這些成果應(yīng)用到高速列車蛇行監(jiān)測領(lǐng)域。宋興武[6]提出橫向位移峰值方法(LMP),通過計(jì)算輪對和構(gòu)架的橫向加速度確定輪對橫移量的峰值來識(shí)別車輛的蛇行運(yùn)動(dòng)。樸明偉等[7]發(fā)現(xiàn)當(dāng)車輛系統(tǒng)服從超臨界分岔時(shí),隨著運(yùn)行速度的增加,轉(zhuǎn)向架會(huì)發(fā)生小幅值蛇行,且隨著速度的增加,轉(zhuǎn)向架的蛇行幅值變大,最終會(huì)出現(xiàn)大幅蛇行失穩(wěn)。蔡里軍[8]總結(jié)了大量的實(shí)測數(shù)據(jù)認(rèn)為當(dāng)構(gòu)架橫向加速度峰值連續(xù)6次達(dá)到或超過2 m/s2,車輛發(fā)生了蛇行運(yùn)動(dòng)。方明寬[9]等考慮到構(gòu)架橫向加速度數(shù)據(jù)的非平穩(wěn)性,提出NKJADE方法對多傳感器橫向加速度數(shù)據(jù)進(jìn)行特征融合,從而對車輛的正常運(yùn)行,小幅蛇行,大幅蛇行狀態(tài)進(jìn)行識(shí)別。Sun[10]發(fā)現(xiàn)當(dāng)車輛出現(xiàn)蛇行運(yùn)動(dòng)時(shí),構(gòu)架與車體的加速度信號(hào)存在相位延遲,通過對車體和構(gòu)架的橫向加速度信號(hào)進(jìn)行互相關(guān)分析,將互相關(guān)指標(biāo)作為判別車輛蛇行失穩(wěn)的閥值。曾元辰、趙飛[11-12]等發(fā)現(xiàn)車輛出現(xiàn)蛇行運(yùn)動(dòng)時(shí),輪對、構(gòu)架和車體的橫向加速度都具有較強(qiáng)的周期性,提出利用信號(hào)周期性特點(diǎn)來識(shí)別車輛的蛇行運(yùn)動(dòng)。崔萬里[13]通過樣本熵理論以及等距映射算法(ISOMAP)對高速列車構(gòu)架橫向加速度信號(hào)進(jìn)行特征提取,然后通過LS-SVM進(jìn)行車輛小幅蛇行異常識(shí)別,達(dá)到了較好的識(shí)別效果。葉運(yùn)廣[14]通過MEEMD進(jìn)行特征提取與使用LS-SVM對車輛的運(yùn)行狀態(tài)進(jìn)行特征識(shí)別亦取得很好的效果。冉偉[15]因而提出了一種基于EEMD-SVD-LTSA的特征提取框架,用于識(shí)別高速列車小幅蛇行運(yùn)動(dòng)的演變趨勢。陳楊[15]通過建立道岔模型,研究高速列車在小幅蛇行的狀態(tài)下通過道岔的演變規(guī)律。寧云志[15]提出基于1D-CNN和CGAN的預(yù)測方法對小幅蛇行進(jìn)行預(yù)警,實(shí)現(xiàn)了蛇行失穩(wěn)的提前預(yù)警。王曉東[18]利用1D CNN自適應(yīng)地對構(gòu)架橫向加速度數(shù)據(jù)進(jìn)行特征提取,使用LSTM對車輛的正常、小幅蛇行和大幅蛇行狀態(tài)進(jìn)行分類,取得了很好的分類效果。但是目前,列車實(shí)際運(yùn)行過程中對蛇行運(yùn)動(dòng)的監(jiān)測還存在以下不足:

    問題一:現(xiàn)有的蛇行監(jiān)測標(biāo)準(zhǔn)將構(gòu)架橫向加速度信號(hào)峰值連續(xù)6次達(dá)到8 m/s2時(shí)視為大幅蛇行,但實(shí)際過程中即使構(gòu)架橫向加速度沒有達(dá)到現(xiàn)有的監(jiān)測標(biāo)準(zhǔn)也會(huì)出現(xiàn)大幅的振動(dòng)。此外,現(xiàn)有標(biāo)準(zhǔn)沒有定量反映蛇行失穩(wěn)對車輛的影響程度,僅用一個(gè)固定值來評判是否產(chǎn)生蛇行,具有一定的局限性。

    問題二:某些車輛受到部分軌道不平順周期性的干擾,雖然發(fā)生了明顯的諧波振動(dòng),但只是運(yùn)行中的極少時(shí)刻,且會(huì)快速收斂,不會(huì)形成穩(wěn)定的周期蛇行運(yùn)動(dòng),對運(yùn)行安全性沒有較大的影響[19],將此類狀態(tài)認(rèn)為快速蛇行收斂。

    本文針對以上問題,提出HHT(Hilbert-Huang transform)能量法從信號(hào)頻域主頻的大小與頻譜的集中性以及頻率值等方面來判斷高速列車是否存在蛇行運(yùn)動(dòng),對兩個(gè)蛇行運(yùn)動(dòng)狀態(tài)(小幅蛇行和大幅蛇行)的蛇行程度的大小進(jìn)行定量的評估。利用最大Lyapunov指數(shù)進(jìn)一步表征信號(hào)的周期性從而將快速蛇行收斂和蛇行運(yùn)動(dòng)(小幅蛇行和大幅蛇行)區(qū)分開來。最后通過HHT能量法與最大Lyapunov指數(shù)相結(jié)合,從高速列車實(shí)際運(yùn)行的需求出發(fā),對車輛系統(tǒng)進(jìn)行正常運(yùn)行、小幅蛇行、快速蛇行收斂、大幅蛇行的定性識(shí)別和蛇行程度大小的定量分析,以實(shí)現(xiàn)對蛇行運(yùn)動(dòng)的具體監(jiān)測。

    1 仿真模型建立

    1.1 高速列車動(dòng)力模型

    本文使用SIMPACK軟件建立國內(nèi)某型號(hào)高速列車的整車動(dòng)力學(xué)模型,該車輛的基本參數(shù)如表1所示。該車輛包括1個(gè)車體、2個(gè)構(gòu)架、4個(gè)輪對、8個(gè)軸箱、一系懸掛系統(tǒng)(一系鋼簧、一系垂向減振器、轉(zhuǎn)臂式軸箱定位裝置)和二系懸掛系統(tǒng)(空氣彈簧、二系垂向減振器、二系橫向減振器、抗蛇行減振器、抗側(cè)滾扭桿、牽引拉桿、橫向止擋)。在建立仿真模型時(shí),忽略了車輛系統(tǒng)各部件的彈性變形,將車體、構(gòu)架、輪對、軸箱等部件假定為剛體。車體、構(gòu)架和輪對有6個(gè)自由度(縱向、橫向、垂向、側(cè)滾、點(diǎn)頭、搖頭),軸箱只有1個(gè)點(diǎn)頭自由度。因此,該車輛動(dòng)力學(xué)模型共有50個(gè)自由度。

    表1 高速列車動(dòng)力學(xué)參數(shù)表

    1.2 建立仿真工況類型

    關(guān)于輪對踏面的選擇,踏面類型對高速列車系統(tǒng)極限環(huán)分岔影響顯著,為了模擬高速列車的不同運(yùn)行狀態(tài),同時(shí)使高速列車的運(yùn)行狀態(tài)更符合實(shí)際,本文分別選用了在一個(gè)旋修周期(2.5×105km)內(nèi)3種不同行駛里程的S1002G實(shí)測踏面。

    踏面1至踏面3對應(yīng)行駛的里程依次增大,踏面的磨損程度依次增加。其中,踏面1為列車行駛5×104km的行駛里程左右時(shí),測量輪對外形輪廓獲得的實(shí)測磨損踏面;踏面2為列車行駛 10×104km 的行駛里程左右時(shí),測量輪對外形輪廓獲得的實(shí)測磨損踏面;踏面3為列車行駛15×104km 的行駛里程左右時(shí),測量輪對外形輪廓獲得的實(shí)測磨損踏面。

    關(guān)于抗蛇行減振器的參數(shù)設(shè)定,抗蛇行減振器的參數(shù)對高速列車系統(tǒng)的臨界速度和分岔類型也有很大影響。本文為充分模擬高速列車的不同運(yùn)行狀態(tài),使用3種不同阻尼特性的抗蛇行減振器。抗蛇行減振器的阻尼特性曲線如圖1所示,從抗蛇行減振器1到抗蛇行減振器3對應(yīng)的卸荷力和卸荷速度依次增加。使用不同行駛里程的S1002G實(shí)測踏面和不同阻尼特性的抗蛇行減振器模擬高速列車(速度:300~400 km/h)3種工況:

    圖1 抗蛇行減振器阻尼系數(shù)圖

    1)工況1:使用踏面1和抗蛇行減振器1仿真,模擬正常運(yùn)行;

    2)工況2:使用踏面2和抗蛇行減振器2仿真,模擬小幅蛇行;

    3)工況3:使用踏面3和抗蛇行減振器3仿真,模擬大幅蛇行和快速蛇行收斂。

    根據(jù)以上工況,計(jì)算得到各工況下的極限環(huán)分岔圖,如圖2所示。工況1的高速列車動(dòng)力學(xué)模型分岔類型為超臨界分岔,在此工況下車輛模型在300~400 km/h區(qū)間運(yùn)行時(shí),車輛處于正常運(yùn)行狀態(tài),如圖3所示。工況2的高速列車動(dòng)力學(xué)模型分岔類型同樣為超臨界分岔,在此工況下車輛模型在300~400 km/h區(qū)間運(yùn)行時(shí),車輛的輪對處于小幅振動(dòng)狀態(tài),此時(shí)車輛處于小幅周期蛇行狀態(tài),隨著速度增加,橫向振動(dòng)幅值也逐漸增大,如圖4所示。工況3下高速列車動(dòng)力學(xué)模型分岔類型為亞臨界分岔,在此工況下,車輛模型在300~400 km/h區(qū)間運(yùn)行時(shí),車輛系統(tǒng)極限環(huán)處于不穩(wěn)定區(qū)域,受軌道不平順的影響,軌道不平順激擾大時(shí),車輛系統(tǒng)平衡點(diǎn)收斂為穩(wěn)定極限環(huán)值產(chǎn)生大幅蛇行運(yùn)動(dòng),如圖5所示。當(dāng)軌道不平順激擾不足以讓車輛系統(tǒng)產(chǎn)生穩(wěn)定的極限環(huán)時(shí),車輛系統(tǒng)可能會(huì)出現(xiàn)短暫的蛇行運(yùn)動(dòng)隨后車輛系統(tǒng)平衡點(diǎn)收斂為0,無法形成穩(wěn)定的周期運(yùn)動(dòng),即快速蛇行收斂,如圖6所示。

    圖2 車輛系統(tǒng)極限環(huán)分岔圖

    圖3 正常運(yùn)行

    圖4 小幅蛇行

    圖5 大幅蛇行

    圖6 快速蛇行收斂

    通過建立上述工況,確定了高速列車動(dòng)力學(xué)模型4種運(yùn)行狀態(tài)的仿真條件。為構(gòu)建大量的仿真數(shù)據(jù)集,將車輛運(yùn)行速度分為(300 km/h、320 km/h、340 km/h、360 km/h、380 km/h、400 km/h)6個(gè)速度級(jí),對上文用到的軌道不平順分別設(shè)置了0.5、1.0、2.0三種幅值比例系數(shù),擴(kuò)充軌道不平順數(shù)據(jù)集。每種運(yùn)行狀態(tài)對應(yīng)6個(gè)速度級(jí),3種幅值比例軌道不平順,共18組仿真條件,4種運(yùn)行狀態(tài)共對應(yīng)72組仿真數(shù)據(jù)。對每組仿真得到的構(gòu)架橫向加速度信號(hào)隨機(jī)取20個(gè)樣本,每個(gè)樣本段取樣時(shí)長為4 s,采樣頻率為250 Hz,總計(jì)1 440個(gè)樣本段的仿真數(shù)據(jù)集,供后續(xù)實(shí)驗(yàn)使用。

    2 理論方法

    2.1 HHT能量理論

    大量的研究表明蛇行信號(hào)和正常信號(hào)的主頻不同,基于這一特點(diǎn),通過研究信號(hào)所有頻率成分在頻域段內(nèi)的變化,可以更加準(zhǔn)確地找到正常信號(hào)和蛇行信號(hào)在頻域內(nèi)的不同特征。再結(jié)合信號(hào)的時(shí)域特征進(jìn)行分析,可以準(zhǔn)確地得知信號(hào)的時(shí)頻信息,從而更加全面地分析信號(hào)特征。孫新建[20]通過HHT能量方法對非平穩(wěn)的爆炸震動(dòng)信號(hào)進(jìn)行分析,對希爾伯特邊際譜進(jìn)行頻域積分計(jì)算HHT能量,研究結(jié)構(gòu)體的爆炸震動(dòng)損傷。HHT能量法結(jié)合了信號(hào)頻域主頻的大小與頻譜的集中性以及頻率值和時(shí)域幅值的相關(guān)信息,因此本文采用此方法研究蛇行運(yùn)動(dòng),

    首先對蛇行信號(hào)和正常信號(hào)的主頻及頻域范圍進(jìn)行分析,確定蛇行頻域主要范圍。然后通過HHT邊際譜統(tǒng)計(jì)正常信號(hào)和蛇行信號(hào)在特定頻率段范圍內(nèi)的頻域信息,最后結(jié)合信號(hào)的時(shí)域幅值和頻域信息,通過HHT能量法全面反映蛇行信號(hào)和正常信號(hào)不同的時(shí)頻特征,從而對車輛運(yùn)行狀態(tài)準(zhǔn)確的識(shí)別。

    HHT具體變換原理如下:對原始信號(hào)x(t)進(jìn)行經(jīng)驗(yàn)?zāi)J椒纸獾玫絥個(gè)IMF(intrinsic mode function)信號(hào)分量:

    (1)

    式中,ci為信號(hào)的基本模式分量,rn稱為殘余函數(shù)。對ci進(jìn)行Hilbert變換得到:

    (2)

    式中,τ為時(shí)間延長間隔。構(gòu)造解析信號(hào)zi(t):

    zi(t)=ci(t)+jGi(t)=ai(t)ejΦi(t)

    (3)

    式中,幅值函數(shù)為ai(t),相位函數(shù)為Φi(t),ci分量信號(hào)的Hilbert譜為Hi(ω,t):

    Hi(ω,t)=RP[ai(t)ejΦi(t)dt]

    (4)

    式中,RP表示取實(shí)部。對Hi(ω,t)進(jìn)行時(shí)域上的積分得到邊際譜hi(ω):

    (5)

    式中,T為信號(hào)時(shí)間長度。對hi(ω)(i=1,2,…,n)邊際譜分析,對主頻在2 Hz以上的邊際譜進(jìn)行疊加求得最終的邊際譜h(ω):

    (6)

    式中,hk(ω)(k=1,2,...j)為主頻在2 Hz以上的邊際譜,在頻域內(nèi)對h(ω)進(jìn)行積分得到最終的HHT能量值。

    (7)

    2.2 最大李雅普諾夫指數(shù)理論

    蛇行信號(hào)和正常信號(hào)除了時(shí)頻上的差異,其周期性特征也有著明顯的區(qū)別,所以有必要對信號(hào)的周期性特征進(jìn)行分析。李雅普諾夫方法是研究運(yùn)動(dòng)穩(wěn)定性的重要方法,對于非線性動(dòng)力學(xué)系統(tǒng),李雅普諾夫指數(shù)可以從狀態(tài)變量的時(shí)間序列中計(jì)算出來。最大李雅普諾夫指數(shù)為相空間相鄰軌跡的平均指數(shù)發(fā)散率的數(shù)值特征,常用來判定系統(tǒng)的混沌特性,因此其可以用來量化信號(hào)的周期性。最大李雅普諾夫指數(shù)計(jì)算步驟如下。

    對于時(shí)間序列{x(t)}構(gòu)造n維空間Rn,T為延遲時(shí)間,T=k△t(k=1,2,…),△t為時(shí)間間隔即:

    X(t)=[x(t),x(t+T),...,x(t+(n-1)T)]T

    (8)

    N維相空間中的某一時(shí)刻,兩條鄰近軌跡之間的距離可以分解在n個(gè)不同的方向,這個(gè)n個(gè)不同方向的上的距離增長率是不同的,每一個(gè)增長率就是一個(gè)李雅普諾夫指數(shù)。系統(tǒng)可以寫成n個(gè)自治一階微分方程組的形式即:

    (9)

    取兩條鄰近軌跡L1和L2,起始點(diǎn)分別為x0和y0兩起始點(diǎn)的距離為d0=y0-x0,經(jīng)過△t時(shí)間后分別運(yùn)動(dòng)到x1和y1,此時(shí)距離為d1=y1-x1,如此循環(huán)下去經(jīng)過m△t后得到m個(gè)di(i=1,2,…,m),di=yi-xi,最大李雅普諾夫指數(shù)λ1即:

    (10)

    3 基于仿真數(shù)據(jù)的蛇行狀態(tài)研究

    3.1 研究步驟

    由于實(shí)測數(shù)據(jù)較少,本文利用仿真數(shù)據(jù)進(jìn)行分析,再通過實(shí)測數(shù)據(jù)進(jìn)行驗(yàn)證,通過HHT能量法和最大Lyapunov指數(shù)法對蛇行狀態(tài)進(jìn)行分類且定量反映蛇行程度的大小,其具體操作如下:

    1)仿真得到不同運(yùn)行狀態(tài)下的構(gòu)架橫向加速度信號(hào)數(shù)據(jù),采用0.5~10 Hz對構(gòu)架橫向加速度進(jìn)行濾波,提取時(shí)間長度為4 s的信號(hào)進(jìn)行分析。

    2)通過對大量的仿真數(shù)據(jù)進(jìn)行EMD(empirical mode decomposition)分析,得出蛇行信號(hào)主頻和能量頻帶基本在2 Hz以上。

    3)對主頻在2 Hz以上IMF(intrinsic mode function)分量信號(hào)的Hilbert邊際譜進(jìn)行疊加求和。

    4)在頻域內(nèi)對最終得到的Hilbert邊際譜進(jìn)行積分得到最終的HHT能量值。

    5)利用SVM(support vector machines)分類得到HHT能量閾值,根據(jù)HHT能量值將正常運(yùn)行、小幅蛇行、大幅蛇行區(qū)分開,且通過HHT能量值的大小反映蛇行程度的大小。

    6)利用最大Lyapunov指數(shù)對高速列車不同運(yùn)行狀態(tài)信號(hào)的周期性進(jìn)行分析。

    7)利用SVM分類得到最大李雅普諾夫指數(shù)閾值,根據(jù)最大Lyapunov指數(shù)值將快速蛇行收斂與蛇行運(yùn)動(dòng)(小幅蛇行和大幅蛇行)區(qū)分開來。

    3.2 HHT能量分析

    對圖3、4、5、6的構(gòu)架橫向加速度信號(hào)進(jìn)行示例分析。如表2所示,正常運(yùn)行信號(hào)的主頻和能量頻帶分布比較均勻,主頻2 Hz以上IMF信號(hào)分量能量占比62%,2 Hz以下IMF信號(hào)分量能量值占比38%。如表3和表4所示,蛇行信號(hào)(小幅蛇行、大幅蛇行)主頻高的IMF分量能量占比高,蛇行頻率主頻和能量頻帶基本在2 Hz以上,信號(hào)分量能量值均占比達(dá)到90%以上。如表5所示,快速蛇行收斂信號(hào)因?yàn)榘咝刑卣?,IMF分量主頻在2 Hz以上的占比也高達(dá)80%。對主頻在2 Hz以上IMF分量信號(hào)邊際譜進(jìn)行疊加,得到最終的Hilbert邊際譜,最后在頻域內(nèi)對Hilbert邊際譜進(jìn)行積分得到最終的HHT能量值。

    表2 正常運(yùn)行(2 Hz以上HHT能量所占百分比為62.1%)

    表3 小幅蛇行(2 Hz以上HHT能量所占百分比為95.87%)

    表4 大幅蛇行(2 Hz以上HHT能量所占百分比為99.5%)

    表5 快速蛇行收斂(2 Hz以上HHT能量所占百分比為82.9%)

    對仿真數(shù)據(jù)集(1 440個(gè)樣本)的正常運(yùn)行、小幅蛇行和大幅蛇行3類數(shù)據(jù)的HHT能量值進(jìn)行分析,利用SVM方法計(jì)算得出其分類閾值,分別為0.82和2.34,如圖7所示。當(dāng)信號(hào)HHT能量值小于0.82時(shí),將其視為正常運(yùn)行,當(dāng)信號(hào)HHT能量值大于0.82小于2.34時(shí)視為小幅蛇行,當(dāng)信號(hào)HHT能量值大于2.34時(shí)視為大幅蛇行。HHT能量值的大小反映了蛇行程度的大小,即HHT能量值越大蛇行程度越劇烈,HHT能量值越小蛇行程度越小。如圖8所示,在此3類數(shù)據(jù)上加入快速蛇行收斂數(shù)據(jù),發(fā)現(xiàn)快速蛇行收斂HHT能量值和小幅蛇行、大幅蛇行HHT能量值重疊,表明HHT能量值指標(biāo)無法區(qū)分這3種運(yùn)行狀態(tài),因此需要補(bǔ)充新的指標(biāo)將其區(qū)分。

    圖7 三類仿真數(shù)據(jù)HHT能量統(tǒng)計(jì)圖

    圖8 四類仿真數(shù)據(jù)HHT能量統(tǒng)計(jì)圖

    大幅蛇行和小幅蛇行的最大Lyapunov指數(shù)值信號(hào)周期性較強(qiáng),正常運(yùn)行和快速蛇行收斂的信號(hào)周期性較弱。因此可通過最大Lyapunov指數(shù)值將快速蛇行收斂和正常運(yùn)行歸為周期性較弱的一類,將小幅蛇行和大幅蛇行歸為周期性較強(qiáng)的一類。

    對仿真數(shù)據(jù)集(1 440個(gè)樣本)的正常運(yùn)行、小幅蛇行、大幅蛇行、快速蛇行收斂數(shù)據(jù)的最大Lyapunov指數(shù)進(jìn)行分析,如圖9所示。利用SVM方法計(jì)算得出其分類閾值10.26,李雅普諾夫指數(shù)值低于10.26時(shí),將其認(rèn)為屬于正常運(yùn)行和快速蛇行收斂一類,最大Lyapunov指數(shù)值大于10.26時(shí)將其認(rèn)為小幅蛇行和大幅蛇行一類。

    圖9 仿真數(shù)據(jù)最大Lyapunov指數(shù)統(tǒng)計(jì)圖

    通過最大Lyapunov指數(shù)法和HHT能量法結(jié)合,將正常運(yùn)行、小幅蛇行、大幅蛇行、快速蛇行收斂區(qū)分開來。當(dāng)信號(hào)HHT能量值小于0.82時(shí),將其視為正常運(yùn)行,當(dāng)信號(hào)HHT能量值大于0.82小于2.34且最大Lyapunov指數(shù)值小于10.26時(shí)視為小幅蛇行,當(dāng)信號(hào)HHT能量值大于2.34且最大Lyapunov指數(shù)值小于10.26時(shí)視為大幅蛇行,當(dāng)信號(hào)HHT能量值大于0.82且最大Lyapunov指數(shù)值大于10.26時(shí)視為快速蛇行收斂。此方法既考慮了蛇行的周期性特點(diǎn),也考慮了蛇行信號(hào)的頻域特性包括頻率主頻的大小、頻譜的集中性以及頻率值等特點(diǎn),從而可以定性識(shí)別車輛系統(tǒng)不同的運(yùn)行狀態(tài),且能夠定量表示蛇行信號(hào)的能量值從而反映蛇行程度的大小。

    4 實(shí)測數(shù)據(jù)驗(yàn)證

    為驗(yàn)證本文所提方法的正確性,采用實(shí)測數(shù)據(jù)對其進(jìn)行驗(yàn)證。圖10為某高速列車速度在300~400 km/h區(qū)間的部分實(shí)測構(gòu)架橫向加速度數(shù)據(jù),時(shí)長為1 220 s,單個(gè)樣本數(shù)據(jù)長度為4 s,樣本總數(shù)為305。

    圖10 實(shí)測構(gòu)架橫向加速度信號(hào)

    用HHT-最大Lyapunov指數(shù)方法對實(shí)測數(shù)據(jù)進(jìn)行分析,圖11為實(shí)測數(shù)據(jù)的HHT能量值,圖12為實(shí)測數(shù)據(jù)的最大Lyapunov指數(shù)圖。圖13~15為局部信號(hào)監(jiān)測放大圖。

    圖11 實(shí)測數(shù)據(jù)HHT能量圖

    圖12 實(shí)測數(shù)據(jù)Lyapunov指數(shù)圖

    圖13 小幅-大幅監(jiān)測圖

    從圖13的識(shí)別結(jié)果中可以觀察到,704~708 s這一段信號(hào),傳統(tǒng)的方法會(huì)將其識(shí)別為正常運(yùn)行,但這段信號(hào)出現(xiàn)了明顯的周期性特征,本文所提方法將其識(shí)別為小幅蛇行狀態(tài),這段小幅蛇行識(shí)別可以起到預(yù)警作用。708~712 s這一段信號(hào),其加速度峰值未達(dá)到目前所使用的蛇行報(bào)警標(biāo)準(zhǔn)(連續(xù)6個(gè)周期加速度峰值達(dá)到8 m/s2),傳統(tǒng)的方法會(huì)將其識(shí)別為正常運(yùn)行,但通過觀察可以看出,該信號(hào)峰值(6~8 m/s2之間)接近目前所使用的蛇行報(bào)警標(biāo)準(zhǔn),且出現(xiàn)了明顯的周期性,本文所提方法將其識(shí)別為大幅蛇行狀態(tài),故本文認(rèn)為其識(shí)別結(jié)果合理。

    從圖14的識(shí)別結(jié)果可以看到712~716 s這一段信號(hào),傳統(tǒng)的方法無法識(shí)別其小幅蛇行狀態(tài),但其具有明顯的周期性特征,716~720 s這一段信號(hào)本方法識(shí)別為快速蛇行收斂,從圖中可以看出其周期性特征逐漸消失,證明了本方法的合理性。

    圖14 小幅-快速蛇行收斂監(jiān)測圖

    從圖15的識(shí)別結(jié)果看到,1 084~1 092 s這一段信號(hào)從正常運(yùn)行變?yōu)樾》咝械臓顟B(tài)。其信號(hào)周期性特征逐漸變強(qiáng),幅值不斷上升,車輛系統(tǒng)出現(xiàn)小幅蛇行。傳統(tǒng)的方法會(huì)將這一段信號(hào)整體識(shí)別為正常運(yùn)行,本方法可以準(zhǔn)確地監(jiān)測這一變化過程。

    圖15 正常-小幅監(jiān)測圖

    本文所提方法準(zhǔn)確地識(shí)別了車輛蛇行運(yùn)動(dòng)的不同狀態(tài),能夠及時(shí)采取措施預(yù)防列車發(fā)生嚴(yán)重的蛇行運(yùn)動(dòng),保證列車運(yùn)行安全。綜上,本文利用實(shí)際數(shù)據(jù)驗(yàn)證了方法具有可行性,能夠在車輛運(yùn)行過程中實(shí)現(xiàn)更為準(zhǔn)確的在線監(jiān)測。

    5 結(jié)束語

    現(xiàn)有的高速列車蛇行失穩(wěn)標(biāo)準(zhǔn)僅針對大幅蛇行且單一的從信號(hào)時(shí)域的角度分析。本文將HHT能量值和最大Lyapunov指數(shù)值兩個(gè)指標(biāo)進(jìn)行了結(jié)合,綜合考慮了信號(hào)的時(shí)域、頻域和周期性特點(diǎn),以此作為區(qū)分高速列車不同蛇行狀態(tài)的參考指標(biāo)。本方法重點(diǎn)并非確定不同蛇行運(yùn)行狀態(tài)的準(zhǔn)確閾值,由于仿真數(shù)據(jù)的局限性,本文通過仿真數(shù)據(jù)所確立的閾值標(biāo)準(zhǔn)存在一定的誤差。閾值的準(zhǔn)確性問題可通過臺(tái)架試驗(yàn)或線路實(shí)測等方式對真實(shí)數(shù)據(jù)進(jìn)行補(bǔ)充,再結(jié)合遷移學(xué)習(xí)等機(jī)器學(xué)習(xí)算法,提高實(shí)際監(jiān)測中閾值的準(zhǔn)確性。

    猜你喜歡
    蛇行踏面構(gòu)架
    建筑安裝造價(jià)控制核心要點(diǎn)構(gòu)架
    踏面清掃器動(dòng)作邏輯與成本控制
    急診PCI治療急性心肌梗死的護(hù)理探索構(gòu)架
    高速列車抗蛇行減振器故障分析
    機(jī)車蛇行狀態(tài)橫向平穩(wěn)性仿真研究
    高速列車可變阻尼抗蛇行減振器適應(yīng)性研究
    基于機(jī)器視覺的礦車踏面磨耗檢測設(shè)計(jì)
    高可靠全平臺(tái)ICT超融合云構(gòu)架的設(shè)計(jì)與實(shí)現(xiàn)
    電子制作(2018年17期)2018-09-28 01:56:54
    基于小幅蛇行視角的地鐵車輛橫向穩(wěn)定性的研究
    略論意象間的主體構(gòu)架
    中華詩詞(2018年1期)2018-06-26 08:46:42
    色视频在线一区二区三区| 麻豆精品久久久久久蜜桃| 97在线人人人人妻| 看免费成人av毛片| 看非洲黑人一级黄片| 国产视频首页在线观看| 少妇人妻精品综合一区二区| 久久久久久久亚洲中文字幕| 大陆偷拍与自拍| 欧美精品一区二区免费开放| 国国产精品蜜臀av免费| 51国产日韩欧美| 99热这里只有是精品50| 三级国产精品欧美在线观看| 久久久国产一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲国产成人一精品久久久| 不卡视频在线观看欧美| 一区二区三区精品91| 亚洲欧洲日产国产| 欧美精品亚洲一区二区| 午夜免费观看性视频| 夜夜看夜夜爽夜夜摸| 亚洲欧美成人综合另类久久久| 精品少妇久久久久久888优播| 国产成人aa在线观看| 夜夜爽夜夜爽视频| 身体一侧抽搐| 精品久久久久久久久av| 水蜜桃什么品种好| 欧美高清成人免费视频www| 交换朋友夫妻互换小说| 一级毛片电影观看| 人妻系列 视频| 青春草视频在线免费观看| 国产av码专区亚洲av| 日韩制服骚丝袜av| 丰满少妇做爰视频| 久久精品久久久久久久性| 一级av片app| 久久久久国产精品人妻一区二区| 精品一区二区三区视频在线| 韩国av在线不卡| 国产美女午夜福利| 一级av片app| 久久国产精品男人的天堂亚洲 | 国产精品麻豆人妻色哟哟久久| 成人午夜精彩视频在线观看| 日本黄色日本黄色录像| 交换朋友夫妻互换小说| 久久久久精品久久久久真实原创| 精品久久久久久久久av| 最近的中文字幕免费完整| 亚洲丝袜综合中文字幕| 午夜福利在线在线| 校园人妻丝袜中文字幕| 国产91av在线免费观看| 国产av精品麻豆| 亚洲国产精品专区欧美| 一区二区三区免费毛片| 亚洲av在线观看美女高潮| 看十八女毛片水多多多| 亚州av有码| 中文字幕精品免费在线观看视频 | 免费av中文字幕在线| 久久久久精品性色| 97超碰精品成人国产| 男女边摸边吃奶| 亚洲国产欧美人成| 韩国高清视频一区二区三区| 亚洲,欧美,日韩| 亚洲欧美精品专区久久| 久久久久久久精品精品| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| 日韩中文字幕视频在线看片 | 天堂俺去俺来也www色官网| 国产成人a区在线观看| 伊人久久精品亚洲午夜| 国产伦精品一区二区三区视频9| 毛片一级片免费看久久久久| 国产精品麻豆人妻色哟哟久久| 日本黄大片高清| 777米奇影视久久| 五月天丁香电影| 国产 一区 欧美 日韩| 九九在线视频观看精品| 国产国拍精品亚洲av在线观看| 久久久久网色| 日韩成人伦理影院| 久久精品国产自在天天线| 久久 成人 亚洲| 免费黄网站久久成人精品| 少妇精品久久久久久久| 99精国产麻豆久久婷婷| 婷婷色综合大香蕉| 噜噜噜噜噜久久久久久91| 国产乱来视频区| 91狼人影院| 国产一区二区在线观看日韩| 国产精品久久久久久av不卡| 夜夜爽夜夜爽视频| 欧美日韩在线观看h| 天美传媒精品一区二区| 尤物成人国产欧美一区二区三区| 交换朋友夫妻互换小说| 最近中文字幕高清免费大全6| 熟妇人妻不卡中文字幕| 人体艺术视频欧美日本| 国产精品爽爽va在线观看网站| 久久久久性生活片| 欧美激情国产日韩精品一区| 国产成人91sexporn| 春色校园在线视频观看| 中文资源天堂在线| 少妇被粗大猛烈的视频| 丰满少妇做爰视频| av专区在线播放| 自拍偷自拍亚洲精品老妇| 日韩成人av中文字幕在线观看| 熟妇人妻不卡中文字幕| 国产精品一区二区性色av| 麻豆国产97在线/欧美| 国产中年淑女户外野战色| 精品少妇黑人巨大在线播放| 色吧在线观看| 亚洲精品成人av观看孕妇| 国产探花极品一区二区| 丰满人妻一区二区三区视频av| 一本一本综合久久| 国产在线视频一区二区| 国产精品秋霞免费鲁丝片| 王馨瑶露胸无遮挡在线观看| 99久久精品一区二区三区| 日产精品乱码卡一卡2卡三| 天美传媒精品一区二区| 日本欧美视频一区| 街头女战士在线观看网站| 久久99热这里只频精品6学生| 直男gayav资源| 久久人妻熟女aⅴ| 日本猛色少妇xxxxx猛交久久| 建设人人有责人人尽责人人享有的 | 亚洲图色成人| 男女免费视频国产| 国产精品女同一区二区软件| 国产精品麻豆人妻色哟哟久久| 中文字幕av成人在线电影| 国产精品精品国产色婷婷| 女的被弄到高潮叫床怎么办| 国产免费福利视频在线观看| 久久人人爽av亚洲精品天堂 | 又粗又硬又长又爽又黄的视频| 久久久久性生活片| 你懂的网址亚洲精品在线观看| 一边亲一边摸免费视频| 欧美丝袜亚洲另类| 国产精品一区二区在线观看99| 国产亚洲91精品色在线| 热re99久久精品国产66热6| 国产精品一区二区三区四区免费观看| 久久久国产一区二区| 国产亚洲一区二区精品| 国产免费福利视频在线观看| 亚洲欧美日韩卡通动漫| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 搡老乐熟女国产| 久久久久国产精品人妻一区二区| 日本黄色日本黄色录像| 3wmmmm亚洲av在线观看| 99久久精品国产国产毛片| 国产在线男女| 精品久久久久久电影网| 99热全是精品| 久久国产精品大桥未久av | 狂野欧美激情性bbbbbb| 亚洲经典国产精华液单| 日本猛色少妇xxxxx猛交久久| 嘟嘟电影网在线观看| 色综合色国产| 我要看日韩黄色一级片| 国产在线一区二区三区精| 高清黄色对白视频在线免费看 | 女的被弄到高潮叫床怎么办| 寂寞人妻少妇视频99o| 舔av片在线| 成年av动漫网址| 亚洲第一区二区三区不卡| 一级二级三级毛片免费看| 97在线人人人人妻| 欧美精品一区二区免费开放| 国产精品熟女久久久久浪| 在线播放无遮挡| 少妇人妻 视频| 国产淫片久久久久久久久| 伦精品一区二区三区| 国产精品秋霞免费鲁丝片| 熟女人妻精品中文字幕| 久久女婷五月综合色啪小说| 精品一区在线观看国产| 97在线人人人人妻| 日本色播在线视频| 国产老妇伦熟女老妇高清| 国产探花极品一区二区| 男女国产视频网站| 中文字幕人妻熟人妻熟丝袜美| 这个男人来自地球电影免费观看 | av国产精品久久久久影院| 国产日韩欧美在线精品| 亚洲天堂av无毛| 男人和女人高潮做爰伦理| 亚洲性久久影院| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美日韩另类电影网站 | 观看美女的网站| 亚洲第一av免费看| 搡女人真爽免费视频火全软件| 久久国内精品自在自线图片| 国产伦在线观看视频一区| 国产成人免费观看mmmm| 中文字幕制服av| 高清黄色对白视频在线免费看 | 成年av动漫网址| 午夜福利高清视频| 亚洲欧美日韩卡通动漫| 尾随美女入室| 国国产精品蜜臀av免费| 国产成人a区在线观看| av免费观看日本| 在现免费观看毛片| 不卡视频在线观看欧美| 在线观看国产h片| 97在线人人人人妻| 在线观看一区二区三区| 男人狂女人下面高潮的视频| 久久久午夜欧美精品| 插阴视频在线观看视频| 亚洲欧美日韩东京热| 国产男女内射视频| 乱码一卡2卡4卡精品| 日产精品乱码卡一卡2卡三| 日韩免费高清中文字幕av| 日日摸夜夜添夜夜添av毛片| 2018国产大陆天天弄谢| 日日啪夜夜爽| 精品久久久久久久久亚洲| 久久久久精品性色| 天堂中文最新版在线下载| 国产黄片视频在线免费观看| 免费看不卡的av| 久久久久久伊人网av| 免费看av在线观看网站| 亚洲精品久久久久久婷婷小说| 亚洲国产最新在线播放| 黄色配什么色好看| 午夜免费鲁丝| 日韩欧美一区视频在线观看 | 精品久久久精品久久久| 99re6热这里在线精品视频| 日韩中文字幕视频在线看片 | 久久久国产一区二区| 老熟女久久久| 国产真实伦视频高清在线观看| 丝袜脚勾引网站| 久久久久久久久久久免费av| 一区二区三区四区激情视频| 91精品一卡2卡3卡4卡| 亚洲国产欧美在线一区| av黄色大香蕉| 少妇被粗大猛烈的视频| 国产视频内射| 国产精品国产三级专区第一集| 婷婷色综合www| 不卡视频在线观看欧美| 青青草视频在线视频观看| 在线观看免费视频网站a站| 日日撸夜夜添| 国产深夜福利视频在线观看| 国产精品无大码| 麻豆成人午夜福利视频| 大话2 男鬼变身卡| 日韩免费高清中文字幕av| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 国产成人午夜福利电影在线观看| xxx大片免费视频| 免费黄网站久久成人精品| 国产精品成人在线| 亚洲精品成人av观看孕妇| 在线看a的网站| 久久午夜福利片| 亚洲最大成人中文| 在线观看美女被高潮喷水网站| 人妻 亚洲 视频| 高清日韩中文字幕在线| av在线app专区| 日本午夜av视频| 97精品久久久久久久久久精品| 亚洲欧美一区二区三区黑人 | 高清毛片免费看| 夫妻性生交免费视频一级片| 舔av片在线| 久久久久久久亚洲中文字幕| 黄片无遮挡物在线观看| 久久久久网色| 边亲边吃奶的免费视频| 国内精品宾馆在线| 欧美区成人在线视频| 少妇猛男粗大的猛烈进出视频| 欧美日本视频| 成人影院久久| 日韩视频在线欧美| 日本av免费视频播放| 久久久久人妻精品一区果冻| 亚洲国产精品专区欧美| 97精品久久久久久久久久精品| 国产伦在线观看视频一区| 色婷婷av一区二区三区视频| 男人爽女人下面视频在线观看| 国产成人a区在线观看| 日韩av不卡免费在线播放| 丰满乱子伦码专区| 免费看日本二区| 永久免费av网站大全| 欧美xxxx黑人xx丫x性爽| 亚洲中文av在线| 免费大片18禁| 日韩一区二区视频免费看| 日本-黄色视频高清免费观看| 日韩一区二区视频免费看| 只有这里有精品99| 人人妻人人添人人爽欧美一区卜 | 肉色欧美久久久久久久蜜桃| 久久国内精品自在自线图片| 在线天堂最新版资源| 18禁在线无遮挡免费观看视频| 日本黄色日本黄色录像| 亚洲色图综合在线观看| 伦理电影大哥的女人| 啦啦啦中文免费视频观看日本| 人人妻人人爽人人添夜夜欢视频 | 美女福利国产在线 | 偷拍熟女少妇极品色| 午夜视频国产福利| 丰满乱子伦码专区| 久久亚洲国产成人精品v| 伦理电影免费视频| 亚洲国产精品999| 久久影院123| av视频免费观看在线观看| 国产精品三级大全| 18+在线观看网站| 亚洲精品国产成人久久av| 欧美变态另类bdsm刘玥| 亚洲精品国产色婷婷电影| 午夜视频国产福利| 欧美性感艳星| 久久久久久久久久成人| 一区二区三区四区激情视频| 亚洲成色77777| 六月丁香七月| 欧美日韩综合久久久久久| 美女福利国产在线 | 观看美女的网站| 久久久久久久久久成人| 免费观看无遮挡的男女| 国产高清不卡午夜福利| 毛片女人毛片| 天天躁夜夜躁狠狠久久av| 国产片特级美女逼逼视频| 国产女主播在线喷水免费视频网站| 国产老妇伦熟女老妇高清| 国产免费一级a男人的天堂| 久久国内精品自在自线图片| 男的添女的下面高潮视频| 国产美女午夜福利| 日本与韩国留学比较| 国产精品国产三级国产专区5o| 汤姆久久久久久久影院中文字幕| 成年免费大片在线观看| 成人黄色视频免费在线看| 欧美3d第一页| 最黄视频免费看| 熟妇人妻不卡中文字幕| 十分钟在线观看高清视频www | 18禁动态无遮挡网站| 亚洲熟女精品中文字幕| 国产深夜福利视频在线观看| 久久久久视频综合| a级一级毛片免费在线观看| 一级a做视频免费观看| 一级毛片aaaaaa免费看小| 国产亚洲欧美精品永久| 成人午夜精彩视频在线观看| 美女xxoo啪啪120秒动态图| 在线精品无人区一区二区三 | 成人亚洲欧美一区二区av| 国产免费一区二区三区四区乱码| 午夜福利在线观看免费完整高清在| 我要看黄色一级片免费的| 在线看a的网站| 精品人妻熟女av久视频| 汤姆久久久久久久影院中文字幕| 亚洲精品456在线播放app| 3wmmmm亚洲av在线观看| av黄色大香蕉| 国产男女内射视频| 91狼人影院| 久久青草综合色| 岛国毛片在线播放| 18禁裸乳无遮挡免费网站照片| 在线天堂最新版资源| 日韩亚洲欧美综合| 男女边摸边吃奶| 成年免费大片在线观看| 国内揄拍国产精品人妻在线| 新久久久久国产一级毛片| 丰满人妻一区二区三区视频av| 日韩视频在线欧美| 在线观看国产h片| 一个人看视频在线观看www免费| 久久 成人 亚洲| 欧美区成人在线视频| 老司机影院成人| 亚洲精品亚洲一区二区| 如何舔出高潮| 十分钟在线观看高清视频www | 国产精品一区二区在线不卡| 成人美女网站在线观看视频| a级毛片免费高清观看在线播放| 色5月婷婷丁香| 嫩草影院新地址| 女人久久www免费人成看片| 精品久久久精品久久久| 22中文网久久字幕| 美女视频免费永久观看网站| 美女国产视频在线观看| 成人综合一区亚洲| 欧美一区二区亚洲| 亚洲色图av天堂| 中文字幕精品免费在线观看视频 | 欧美日韩亚洲高清精品| www.色视频.com| 国产色婷婷99| 久久人人爽av亚洲精品天堂 | 国产人妻一区二区三区在| 亚洲精品乱码久久久久久按摩| 成人二区视频| av免费在线看不卡| 国产欧美日韩精品一区二区| 深爱激情五月婷婷| 久久国产精品大桥未久av | 久久婷婷青草| 久久人人爽av亚洲精品天堂 | a 毛片基地| 看免费成人av毛片| 久久精品久久精品一区二区三区| 国产av码专区亚洲av| 精品人妻熟女av久视频| 亚洲成人中文字幕在线播放| 免费黄频网站在线观看国产| 日韩强制内射视频| 黄色一级大片看看| 久久国内精品自在自线图片| 黑丝袜美女国产一区| 久久6这里有精品| 亚洲国产欧美人成| 一本久久精品| 国产成人精品福利久久| 在线观看一区二区三区| 久久国产亚洲av麻豆专区| 97超碰精品成人国产| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品第二区| 久久人妻熟女aⅴ| 日本爱情动作片www.在线观看| 丝袜喷水一区| 久久久久精品久久久久真实原创| 欧美成人精品欧美一级黄| 男女免费视频国产| 一区二区三区免费毛片| 精品少妇黑人巨大在线播放| 精品午夜福利在线看| 国产深夜福利视频在线观看| 美女主播在线视频| 99热这里只有是精品在线观看| 国产高潮美女av| 欧美成人一区二区免费高清观看| 免费大片黄手机在线观看| 亚洲精品国产色婷婷电影| 国产色婷婷99| 久久久成人免费电影| 丝瓜视频免费看黄片| 亚洲成人中文字幕在线播放| 国产中年淑女户外野战色| 我的老师免费观看完整版| 美女cb高潮喷水在线观看| 精品国产露脸久久av麻豆| 久久午夜福利片| 色网站视频免费| 丰满少妇做爰视频| 少妇高潮的动态图| 久久ye,这里只有精品| 成人综合一区亚洲| 成人国产av品久久久| 国产一区有黄有色的免费视频| 免费播放大片免费观看视频在线观看| 欧美成人精品欧美一级黄| av播播在线观看一区| 亚洲成人一二三区av| 只有这里有精品99| a级毛色黄片| av播播在线观看一区| 成人二区视频| 精品人妻视频免费看| 国产无遮挡羞羞视频在线观看| 午夜福利高清视频| 国产精品一区www在线观看| 国产成人精品久久久久久| 亚洲精品日本国产第一区| 亚洲av欧美aⅴ国产| 国语对白做爰xxxⅹ性视频网站| 成人亚洲欧美一区二区av| 国产精品久久久久久av不卡| 五月开心婷婷网| 国产黄片视频在线免费观看| 五月开心婷婷网| 日韩亚洲欧美综合| 黄片wwwwww| 久久国产精品男人的天堂亚洲 | 成人亚洲欧美一区二区av| 免费黄色在线免费观看| 亚洲伊人久久精品综合| 亚洲欧美成人精品一区二区| 久久6这里有精品| 欧美最新免费一区二区三区| 好男人视频免费观看在线| 亚洲在久久综合| 久久综合国产亚洲精品| 人人妻人人爽人人添夜夜欢视频 | 一区二区三区乱码不卡18| 亚洲精品久久久久久婷婷小说| 欧美成人午夜免费资源| 国产精品爽爽va在线观看网站| 久久国产精品大桥未久av | 色5月婷婷丁香| 18禁裸乳无遮挡动漫免费视频| 久久精品熟女亚洲av麻豆精品| 日本av手机在线免费观看| av在线蜜桃| 国产精品久久久久久久久免| 亚洲内射少妇av| 国产白丝娇喘喷水9色精品| 国产 精品1| 草草在线视频免费看| 80岁老熟妇乱子伦牲交| 欧美xxxx黑人xx丫x性爽| 免费高清在线观看视频在线观看| 97在线人人人人妻| 青春草亚洲视频在线观看| 免费av不卡在线播放| 免费黄色在线免费观看| av免费在线看不卡| 国产探花极品一区二区| a 毛片基地| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美成人精品一区二区| 性色avwww在线观看| 婷婷色麻豆天堂久久| 国产 一区 欧美 日韩| a级毛片免费高清观看在线播放| 日本欧美视频一区| 亚洲人成网站在线观看播放| 有码 亚洲区| 国产伦在线观看视频一区| 久久精品久久精品一区二区三区| 亚洲在久久综合| 免费av中文字幕在线| 人人妻人人看人人澡| 成人特级av手机在线观看| av视频免费观看在线观看| 又粗又硬又长又爽又黄的视频| 一级爰片在线观看| 丰满人妻一区二区三区视频av| 黄色欧美视频在线观看| 一边亲一边摸免费视频| 岛国毛片在线播放| 亚洲精品成人av观看孕妇| 爱豆传媒免费全集在线观看| 久久毛片免费看一区二区三区| 嫩草影院入口| 亚洲美女黄色视频免费看| h视频一区二区三区| 天天躁日日操中文字幕| 日产精品乱码卡一卡2卡三| 亚洲第一av免费看| 免费观看无遮挡的男女| 日产精品乱码卡一卡2卡三| 夜夜爽夜夜爽视频| 精品久久久久久久末码| 自拍偷自拍亚洲精品老妇| 男女无遮挡免费网站观看| 一区二区三区免费毛片| 欧美日韩视频精品一区| 亚洲丝袜综合中文字幕| 99久久综合免费| 国产精品偷伦视频观看了| 免费人成在线观看视频色| 天天躁日日操中文字幕| 国产一区二区三区av在线| 久热这里只有精品99| 99久久精品国产国产毛片| 日韩不卡一区二区三区视频在线| 91狼人影院|