• <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
    好男人在线观看高清免费视频| 亚洲美女视频黄频| 日韩av在线大香蕉| 少妇熟女aⅴ在线视频| 在现免费观看毛片| 一级二级三级毛片免费看| 91精品国产九色| 夫妻性生交免费视频一级片| 久久久欧美国产精品| 国产一区有黄有色的免费视频 | a级毛片免费高清观看在线播放| 九色成人免费人妻av| 视频中文字幕在线观看| videossex国产| 日本欧美国产在线视频| 中文字幕免费在线视频6| 精品久久久久久成人av| 中国美白少妇内射xxxbb| 免费无遮挡裸体视频| 欧美区成人在线视频| 十八禁国产超污无遮挡网站| 午夜日本视频在线| 国产淫语在线视频| 高清毛片免费看| 中文字幕制服av| 日韩三级伦理在线观看| 亚洲av中文av极速乱| 国产伦在线观看视频一区| 亚洲美女搞黄在线观看| 在线免费十八禁| 免费看美女性在线毛片视频| 精品国产一区二区三区久久久樱花 | 成年免费大片在线观看| 高清av免费在线| 国产一区二区三区综合在线观看 | 国产精品蜜桃在线观看| 可以在线观看毛片的网站| 精品酒店卫生间| 91午夜精品亚洲一区二区三区| 亚洲自拍偷在线| 国产午夜精品久久久久久一区二区三区| 亚洲国产精品专区欧美| 91在线精品国自产拍蜜月| 美女内射精品一级片tv| 韩国av在线不卡| 国产一区有黄有色的免费视频 | 天天躁夜夜躁狠狠久久av| 岛国毛片在线播放| 毛片女人毛片| 亚洲一级一片aⅴ在线观看| 91av网一区二区| kizo精华| 欧美高清性xxxxhd video| 欧美高清性xxxxhd video| 色尼玛亚洲综合影院| 高清毛片免费看| 99热这里只有精品一区| 亚洲丝袜综合中文字幕| 肉色欧美久久久久久久蜜桃 | 神马国产精品三级电影在线观看| 国产精品一区二区三区四区久久| 欧美日韩国产mv在线观看视频 | 丰满人妻一区二区三区视频av| 久久综合国产亚洲精品| 51国产日韩欧美| 极品教师在线视频| 秋霞伦理黄片| 国产成人精品一,二区| 欧美成人一区二区免费高清观看| 在线观看一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 九九久久精品国产亚洲av麻豆| 国产精品福利在线免费观看| 国产极品天堂在线| 美女脱内裤让男人舔精品视频| 日韩精品青青久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精华霜和精华液先用哪个| 街头女战士在线观看网站| 亚洲国产欧美在线一区| 国模一区二区三区四区视频| 在线 av 中文字幕| 真实男女啪啪啪动态图| 九九爱精品视频在线观看| 看非洲黑人一级黄片| 免费av观看视频| 91精品伊人久久大香线蕉| 中文精品一卡2卡3卡4更新| 91av网一区二区| 久久久久久九九精品二区国产| 插逼视频在线观看| 99久久精品热视频| 激情 狠狠 欧美| 一区二区三区免费毛片| 观看美女的网站| 久久久久网色| 欧美一级a爱片免费观看看| 天堂中文最新版在线下载 | 亚洲熟妇中文字幕五十中出| 99久久精品国产国产毛片| 丝袜喷水一区| 亚洲精品视频女| 国产视频首页在线观看| 亚洲国产成人一精品久久久| 精品不卡国产一区二区三区| 久久久久久伊人网av| 丰满少妇做爰视频| 免费人成在线观看视频色| 亚洲精华国产精华液的使用体验| 亚洲精品日韩在线中文字幕| 日本猛色少妇xxxxx猛交久久| xxx大片免费视频| 插逼视频在线观看| 尤物成人国产欧美一区二区三区| 一本久久精品| 亚洲国产av新网站| freevideosex欧美| 国产精品一区二区在线观看99 | 天美传媒精品一区二区| 国产真实伦视频高清在线观看| www.av在线官网国产| 日韩国内少妇激情av| 十八禁网站网址无遮挡 | 亚洲精品一二三| 激情 狠狠 欧美| 久久久久久国产a免费观看| 国产高清有码在线观看视频| 卡戴珊不雅视频在线播放| or卡值多少钱| 18禁在线无遮挡免费观看视频| 老司机影院成人| 91精品国产九色| a级一级毛片免费在线观看| 大陆偷拍与自拍| 国产黄片美女视频| 天堂中文最新版在线下载 | 可以在线观看毛片的网站| 亚洲无线观看免费| 国产精品.久久久| 亚洲激情五月婷婷啪啪| 久久99热这里只有精品18| 国产视频首页在线观看| 99久国产av精品国产电影| 欧美高清性xxxxhd video| 国产精品不卡视频一区二区| 免费在线观看成人毛片| 男人爽女人下面视频在线观看| 免费大片18禁| 偷拍熟女少妇极品色| 久久99热这里只频精品6学生| 亚洲精品影视一区二区三区av| 日本欧美国产在线视频| 亚洲欧美中文字幕日韩二区| 中文字幕制服av| 国产av国产精品国产| 欧美一级a爱片免费观看看| 99热这里只有是精品在线观看| 大片免费播放器 马上看| 777米奇影视久久| 成人毛片a级毛片在线播放| 欧美97在线视频| 亚洲精品第二区| 国产高清国产精品国产三级 | 国产黄频视频在线观看| 国产精品福利在线免费观看| 日韩精品有码人妻一区| 插阴视频在线观看视频| 又爽又黄a免费视频| h日本视频在线播放| 少妇人妻一区二区三区视频| 91在线精品国自产拍蜜月| 在线 av 中文字幕| 日韩在线高清观看一区二区三区| 国产高潮美女av| 免费av毛片视频| 亚洲精品国产成人久久av| 一级毛片我不卡| 免费黄网站久久成人精品| 欧美激情久久久久久爽电影| 亚洲一级一片aⅴ在线观看| 国产v大片淫在线免费观看| 卡戴珊不雅视频在线播放| 久久久久久久亚洲中文字幕| 日本色播在线视频| 99久久精品国产国产毛片| 少妇的逼水好多| 国产精品无大码| 精品久久久精品久久久| 尾随美女入室| 七月丁香在线播放| 国产精品伦人一区二区| 免费观看av网站的网址| 淫秽高清视频在线观看| 久久精品久久久久久久性| 欧美最新免费一区二区三区| 天堂√8在线中文| 超碰97精品在线观看| 黄色日韩在线| 中国国产av一级| 国产久久久一区二区三区| 日韩欧美 国产精品| 青春草亚洲视频在线观看| 精品久久国产蜜桃| 国内精品美女久久久久久| 99久国产av精品| 日本黄色片子视频| 国产精品伦人一区二区| 秋霞伦理黄片| 小蜜桃在线观看免费完整版高清| 婷婷色av中文字幕| 91av网一区二区| 亚洲av日韩在线播放| 欧美最新免费一区二区三区| 国产永久视频网站| 精品不卡国产一区二区三区| 国产乱人视频| 精品久久国产蜜桃| 久久久久精品性色| 老女人水多毛片| 国产精品一区二区在线观看99 | 一个人免费在线观看电影| 久久99热6这里只有精品| 国产亚洲精品av在线| av播播在线观看一区| 搡老妇女老女人老熟妇| 丰满少妇做爰视频| h日本视频在线播放| 亚洲熟女精品中文字幕| 校园人妻丝袜中文字幕| 国产成人a∨麻豆精品| 精品人妻熟女av久视频| 一本一本综合久久| 欧美丝袜亚洲另类| 精品国内亚洲2022精品成人| 成年女人看的毛片在线观看| 高清日韩中文字幕在线| 亚洲综合精品二区| av在线老鸭窝| 韩国av在线不卡| 亚洲精品,欧美精品| 婷婷色麻豆天堂久久| 中文天堂在线官网| av福利片在线观看| 久久久久久九九精品二区国产| 天堂√8在线中文| 久99久视频精品免费| 欧美一级a爱片免费观看看| 99re6热这里在线精品视频| 日本黄大片高清| 少妇被粗大猛烈的视频| 国产视频首页在线观看| 日日干狠狠操夜夜爽| 国产淫语在线视频| 免费观看a级毛片全部| 一级毛片我不卡| 天堂av国产一区二区熟女人妻| 精品熟女少妇av免费看| 亚洲乱码一区二区免费版| 青春草视频在线免费观看| 精品少妇黑人巨大在线播放| 岛国毛片在线播放| 波野结衣二区三区在线| av在线观看视频网站免费| 国产精品女同一区二区软件| 日本黄色片子视频| 九草在线视频观看| 日本一本二区三区精品| 人妻系列 视频| 亚洲精品乱码久久久v下载方式| 91精品一卡2卡3卡4卡| 免费高清在线观看视频在线观看| 搡老乐熟女国产| 午夜激情欧美在线| 中文字幕av在线有码专区| 国产免费福利视频在线观看| 男女边吃奶边做爰视频| 丝袜喷水一区| www.色视频.com| 国产成人aa在线观看| 看免费成人av毛片| 国产亚洲最大av| 99热这里只有是精品50| 搡女人真爽免费视频火全软件| 在线免费观看不下载黄p国产| 在线天堂最新版资源| 91在线精品国自产拍蜜月| 国产在线男女| 有码 亚洲区| 日韩欧美国产在线观看| 亚洲av成人精品一二三区| 三级国产精品片| 成人毛片60女人毛片免费| 一区二区三区免费毛片| 国产黄片美女视频| 卡戴珊不雅视频在线播放| 欧美日本视频| 久久久久久久久大av| 婷婷色av中文字幕| 亚洲av中文av极速乱| 亚洲综合精品二区| 最近中文字幕2019免费版| 国产黄片视频在线免费观看| 亚洲国产日韩欧美精品在线观看| av一本久久久久| 国产色爽女视频免费观看| 日韩精品有码人妻一区| 亚洲国产成人一精品久久久| 亚洲av电影在线观看一区二区三区 | 超碰av人人做人人爽久久| 亚洲av不卡在线观看| 国产极品天堂在线| 91久久精品国产一区二区成人| 国产伦理片在线播放av一区| 国产高清国产精品国产三级 | 伊人久久国产一区二区| 搡女人真爽免费视频火全软件| 成人鲁丝片一二三区免费| 国产精品综合久久久久久久免费| 大香蕉97超碰在线| 国产高清不卡午夜福利| 免费看美女性在线毛片视频| 日韩成人伦理影院| 搡女人真爽免费视频火全软件| 我要看日韩黄色一级片| av福利片在线观看| 边亲边吃奶的免费视频| 精品久久久久久久久亚洲| 超碰av人人做人人爽久久| 日韩 亚洲 欧美在线| 干丝袜人妻中文字幕| 欧美变态另类bdsm刘玥| 三级经典国产精品| 久久久久久久久大av| 国产片特级美女逼逼视频| 久久久久性生活片| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 女的被弄到高潮叫床怎么办| 国产精品不卡视频一区二区| av在线蜜桃| 高清av免费在线| 色视频www国产| 十八禁网站网址无遮挡 | 亚洲一级一片aⅴ在线观看| 人妻系列 视频| h日本视频在线播放| 日日干狠狠操夜夜爽| 久久久久久久久久成人| 国产高清有码在线观看视频| 99久久九九国产精品国产免费| 日日干狠狠操夜夜爽| 国产探花极品一区二区| 久久99热这里只频精品6学生| 亚洲,欧美,日韩| 少妇的逼好多水| 联通29元200g的流量卡| 日韩三级伦理在线观看| 青春草视频在线免费观看| 国产精品一区二区三区四区久久| 久久精品熟女亚洲av麻豆精品 | 国产一级毛片在线| eeuss影院久久| 国产精品一二三区在线看| 国产精品国产三级专区第一集| 99热全是精品| 婷婷色麻豆天堂久久| 秋霞在线观看毛片| 国产av在哪里看| 麻豆精品久久久久久蜜桃| 男女边摸边吃奶| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 91在线精品国自产拍蜜月| 日韩一区二区视频免费看| 2022亚洲国产成人精品| 一边亲一边摸免费视频| 久久久精品94久久精品| 免费在线观看成人毛片| 又黄又爽又刺激的免费视频.| 色网站视频免费| 日日啪夜夜撸| av卡一久久| 久久久久免费精品人妻一区二区| 亚洲av.av天堂| 蜜桃亚洲精品一区二区三区| 国产免费福利视频在线观看| 色播亚洲综合网| 久久久久久久亚洲中文字幕| 日韩欧美三级三区| 亚洲av成人精品一区久久| 在线免费观看的www视频| 国产 一区 欧美 日韩| 久久久久久久久久黄片| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 久久人人爽人人爽人人片va| 一区二区三区四区激情视频| 麻豆成人午夜福利视频| 亚洲在线自拍视频| 极品少妇高潮喷水抽搐| 成人高潮视频无遮挡免费网站| 亚洲熟妇中文字幕五十中出| 欧美激情在线99| 国产亚洲精品av在线| 99久久中文字幕三级久久日本| 亚洲精品,欧美精品| 亚洲最大成人av| 亚洲av免费高清在线观看| 国产免费又黄又爽又色| 菩萨蛮人人尽说江南好唐韦庄| 午夜久久久久精精品| 国产一区亚洲一区在线观看| 亚洲人成网站在线播| 国内揄拍国产精品人妻在线| 老师上课跳d突然被开到最大视频| 天天躁夜夜躁狠狠久久av| 久久久久性生活片| 两个人的视频大全免费| a级毛色黄片| 91狼人影院| 国产av国产精品国产| 亚洲性久久影院| or卡值多少钱| www.av在线官网国产| 一级片'在线观看视频| 小蜜桃在线观看免费完整版高清| 国内揄拍国产精品人妻在线| 国产一区有黄有色的免费视频 | 国产精品一区二区三区四区久久| 国产精品熟女久久久久浪| 22中文网久久字幕| av在线天堂中文字幕| 人妻制服诱惑在线中文字幕| 一级二级三级毛片免费看| 国产探花在线观看一区二区| 免费电影在线观看免费观看| 国产精品一二三区在线看| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 久久综合国产亚洲精品| 伦理电影大哥的女人| 青春草亚洲视频在线观看| 蜜桃亚洲精品一区二区三区| 美女国产视频在线观看| 欧美成人精品欧美一级黄| 搡女人真爽免费视频火全软件| 国产成人aa在线观看| 最近中文字幕2019免费版| av国产久精品久网站免费入址| 成人综合一区亚洲| 精品一区在线观看国产| 男女那种视频在线观看| 免费看a级黄色片| 一级av片app| 一区二区三区高清视频在线| 久久精品久久精品一区二区三区| 一个人免费在线观看电影| 天堂网av新在线| 亚洲欧美一区二区三区黑人 | 美女主播在线视频| 简卡轻食公司| 久久久久免费精品人妻一区二区| 免费播放大片免费观看视频在线观看| 国内精品宾馆在线| 国产在线男女| 亚洲自拍偷在线| 亚洲av成人av| 非洲黑人性xxxx精品又粗又长| 天堂影院成人在线观看| 欧美激情国产日韩精品一区| 少妇高潮的动态图| 国产午夜精品久久久久久一区二区三区| 一个人观看的视频www高清免费观看| 尾随美女入室| 久久久久久久午夜电影| 久久精品久久久久久久性| 亚洲国产高清在线一区二区三| 69人妻影院| 大陆偷拍与自拍| 国产免费视频播放在线视频 | 欧美精品国产亚洲| 国产伦理片在线播放av一区| 亚洲av成人av| 欧美激情国产日韩精品一区| 亚洲精品国产av成人精品| 国产毛片a区久久久久| 日本色播在线视频| 国产视频内射| 九九久久精品国产亚洲av麻豆| 日本黄色片子视频| 高清日韩中文字幕在线| 波多野结衣巨乳人妻| 一级毛片 在线播放| 小蜜桃在线观看免费完整版高清| 成人午夜高清在线视频| 黄色配什么色好看| 国内精品宾馆在线| 寂寞人妻少妇视频99o| 久久99热6这里只有精品| 亚洲av中文av极速乱| 一级黄片播放器| 人人妻人人澡欧美一区二区| 日韩一区二区三区影片| 草草在线视频免费看| 成人av在线播放网站| 伦理电影大哥的女人| 欧美zozozo另类| 免费少妇av软件| 久久综合国产亚洲精品| 免费av观看视频| 精品久久久久久久末码| 亚洲经典国产精华液单| 在线 av 中文字幕| 两个人的视频大全免费| 99久国产av精品| 亚洲国产精品sss在线观看| 久久久久久久国产电影| 亚洲熟妇中文字幕五十中出| 国产成人精品一,二区| 99九九线精品视频在线观看视频| av天堂中文字幕网| 久久韩国三级中文字幕| 你懂的网址亚洲精品在线观看| 欧美日韩视频高清一区二区三区二| 大香蕉久久网| 91精品伊人久久大香线蕉| 大陆偷拍与自拍| 国产黄色小视频在线观看| 久久久久久久久久久免费av| 欧美一区二区亚洲| 美女脱内裤让男人舔精品视频| 秋霞在线观看毛片| av在线观看视频网站免费| 色尼玛亚洲综合影院| 少妇被粗大猛烈的视频| 成年av动漫网址| 免费av不卡在线播放| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品在线观看| av在线蜜桃| 青春草亚洲视频在线观看| 女的被弄到高潮叫床怎么办| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 老司机影院成人| 人人妻人人澡人人爽人人夜夜 | 99久久精品热视频| 久久久久久久国产电影| 1000部很黄的大片| 日韩 亚洲 欧美在线| 男女啪啪激烈高潮av片| 亚洲最大成人手机在线| 国产探花极品一区二区| 美女主播在线视频| 2021少妇久久久久久久久久久| 亚洲激情五月婷婷啪啪| av黄色大香蕉| av又黄又爽大尺度在线免费看| 亚洲在久久综合| 青春草亚洲视频在线观看| 十八禁网站网址无遮挡 | 极品教师在线视频| 久久人人爽人人爽人人片va| 亚洲最大成人中文| 免费观看无遮挡的男女| 亚洲精品乱码久久久久久按摩| 欧美xxxx性猛交bbbb| 国产精品久久久久久精品电影| 日韩成人伦理影院| 日本爱情动作片www.在线观看| 麻豆久久精品国产亚洲av| 女的被弄到高潮叫床怎么办| 美女国产视频在线观看| 欧美激情国产日韩精品一区| 69av精品久久久久久| 日本-黄色视频高清免费观看| 色视频www国产| 国产一区二区在线观看日韩| 波多野结衣巨乳人妻| 超碰av人人做人人爽久久| 免费无遮挡裸体视频| 少妇熟女欧美另类| 91av网一区二区| 99久久精品一区二区三区| 久久精品综合一区二区三区| 波野结衣二区三区在线| 久久国内精品自在自线图片| 日韩av在线免费看完整版不卡| 寂寞人妻少妇视频99o| 一个人免费在线观看电影| 亚洲欧美精品自产自拍| 亚洲欧美一区二区三区国产| 永久网站在线| 内地一区二区视频在线| 2018国产大陆天天弄谢| 人妻一区二区av| 男人和女人高潮做爰伦理| 日韩av不卡免费在线播放| 观看免费一级毛片| 又粗又硬又长又爽又黄的视频| 国产午夜精品一二区理论片| 一级毛片久久久久久久久女| 日韩欧美 国产精品| a级毛色黄片| 天堂√8在线中文| 午夜视频国产福利| 午夜福利视频1000在线观看| 亚洲精品日韩av片在线观看| 直男gayav资源| 超碰av人人做人人爽久久| 国产黄片美女视频| av国产免费在线观看| 精品99又大又爽又粗少妇毛片| 一级a做视频免费观看| 老司机影院毛片|