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

    考慮波長(zhǎng)因素的軌道不平順預(yù)測(cè)研究

    2016-01-25 06:50:34呂宏李再幃何越磊
    關(guān)鍵詞:預(yù)測(cè)方法

    呂宏 ,李再幃,2,3, 何越磊

    (1.上海工程技術(shù)大學(xué) 城市軌道交通學(xué)院, 上海 201620;

    2.華東交通大學(xué) 鐵路環(huán)境振動(dòng)與噪聲教育部工程研究中心,江西 南昌 330013;

    3.北京交通大學(xué) 土木建筑工程學(xué)院, 北京100044)

    ?

    考慮波長(zhǎng)因素的軌道不平順預(yù)測(cè)研究

    呂宏1,李再幃1,2,3, 何越磊1

    (1.上海工程技術(shù)大學(xué) 城市軌道交通學(xué)院, 上海 201620;

    2.華東交通大學(xué) 鐵路環(huán)境振動(dòng)與噪聲教育部工程研究中心,江西 南昌 330013;

    3.北京交通大學(xué) 土木建筑工程學(xué)院, 北京100044)

    摘要:準(zhǔn)確預(yù)測(cè)軌道不平順的變化規(guī)律是提高軌道養(yǎng)護(hù)維修效率、降低維修成本的關(guān)鍵問(wèn)題之一。鑒于軌道不平順波長(zhǎng)對(duì)列車(chē)運(yùn)行狀態(tài)的影響,基于小波理論,通過(guò)選取合適的小波基函數(shù)及分解層數(shù),對(duì)原始不平順時(shí)域序列進(jìn)行小波分解及重構(gòu),并采用不平順狀態(tài)推移方法對(duì)重構(gòu)后不平順序列進(jìn)行分析處理,利用不同時(shí)刻下軌道不平順的狀態(tài)分布,從而實(shí)現(xiàn)對(duì)軌道不平順預(yù)測(cè);通過(guò)2個(gè)實(shí)際算例對(duì)本方法正確性進(jìn)行驗(yàn)證,并對(duì)預(yù)測(cè)結(jié)果進(jìn)行誤差分析。研究結(jié)果表明:利用本文提出方法得出的預(yù)測(cè)結(jié)果與實(shí)際結(jié)果誤差較小,可以有效地預(yù)測(cè)軌道不平順的變化趨勢(shì),且能夠?qū)崿F(xiàn)不同波長(zhǎng)段軌道不平順的預(yù)測(cè),可為我國(guó)軌道養(yǎng)護(hù)維修方案的制定提供一定的技術(shù)支持。

    關(guān)鍵詞:小波理論;軌道不平順;預(yù)測(cè)方法;狀態(tài)推移法

    軌道不平順是軌道結(jié)構(gòu)形變的外部綜合反映,是產(chǎn)生車(chē)輛系統(tǒng)振動(dòng)的主要根源,同時(shí)也是軌道系統(tǒng)動(dòng)力分析與軌道養(yǎng)護(hù)維修及質(zhì)量評(píng)估中的重要基礎(chǔ)資料,其對(duì)輪軌系統(tǒng)動(dòng)力傳播、輪軌噪聲和線路養(yǎng)護(hù)維修等方面的研究都有十分重要的作用。由于不同波長(zhǎng)的軌道不平順對(duì)于軌道線路的影響不盡相同,對(duì)于高速鐵路而言,在列車(chē)高速運(yùn)行區(qū)段,其長(zhǎng)波不平順對(duì)列車(chē)運(yùn)行狀態(tài)影響較為明顯,而在列車(chē)低速運(yùn)行區(qū)段,其局部的短波不平順影響較為嚴(yán)重,易多激擾起較大輪軌力,對(duì)列車(chē)運(yùn)行安全產(chǎn)生較大影響。因此,利用現(xiàn)有檢測(cè)數(shù)據(jù)實(shí)現(xiàn)不同波長(zhǎng)區(qū)段軌道不平順的預(yù)測(cè)具有十分重要的意義。影響軌道不平順變化的因素較多,各因素的影響程度也各不相同,且由于環(huán)境、路基狀態(tài)等部分因素難以具體描述。因此,國(guó)內(nèi)外學(xué)者針對(duì)軌道不平順變化問(wèn)題展開(kāi)了大量研究:日本學(xué)者[1-3]先后通過(guò)多參數(shù)分析法、回歸分析法、指數(shù)平滑法和頻域分析法等方法,利用軌檢車(chē)數(shù)據(jù)及軌道基本參數(shù)從幅值角度對(duì)不平順預(yù)測(cè)模型進(jìn)行了探索;我國(guó)在軌道不平順預(yù)測(cè)方面也取得了部分研究成果:許玉德等[4]利用線性預(yù)測(cè)模型,對(duì)鐵路軌道質(zhì)量指數(shù)(TQI)的發(fā)展進(jìn)行了預(yù)測(cè)研究;陳憲麥等[5]提出了一種用于軌道不平順發(fā)展預(yù)測(cè)的綜合因子法,對(duì)軌道不平順發(fā)展趨勢(shì)進(jìn)行預(yù)測(cè);曲建軍等[6]提出利用灰色預(yù)測(cè)模型,實(shí)現(xiàn)了對(duì)軌道TQI指數(shù)變化趨勢(shì)的預(yù)測(cè)。高建敏等[7]提出了一種基于不平順?lè)植己瘮?shù)的不平順發(fā)展統(tǒng)計(jì)預(yù)測(cè)方法,實(shí)現(xiàn)軌道不平順的預(yù)測(cè)。上述方法均能有效地預(yù)測(cè)出軌道不平順的變化趨勢(shì),對(duì)鐵路線路養(yǎng)護(hù)維修周期制定起到了積極作用,但值得注意的是,這些研究多是局限于軌道不平順?lè)到嵌冗M(jìn)行預(yù)測(cè),并未涉及到不同波長(zhǎng)軌道不平順的變化趨勢(shì)。而隨著高速鐵路及城市軌道交通的不斷發(fā)展,不同里程段上軌道不平順呈現(xiàn)出的波長(zhǎng)分布特性存在一定的差異。同時(shí),對(duì)于不同里程段上所關(guān)注的軌道不平順波長(zhǎng)范圍也不盡相同。因此,需要一種可以實(shí)現(xiàn)不同波長(zhǎng)特征的軌道不平順預(yù)測(cè)的方法。鑒于小波分解理論具有多頻段分析和適應(yīng)性分析的特點(diǎn)[8],且在不平順局部波型辨識(shí)、軌道病害識(shí)別等方面都有一定地應(yīng)用[9]。本文在已有的研究基礎(chǔ)上,引入小波理論,進(jìn)而提出一種考慮波長(zhǎng)因素的軌道不平順預(yù)測(cè)方法,從而為我國(guó)鐵路養(yǎng)護(hù)維修及工務(wù)部門(mén)的精細(xì)化和高效管理方面提供相應(yīng)技術(shù)支撐。

    1基本原理

    1.1 小波分解及重構(gòu)原理

    對(duì)于函數(shù)x(t)∈L2(R),L2(R)表示一個(gè)平方可積的實(shí)數(shù)空間,x(t)的離散小波變換定義為:

    (1)

    ψj,k=α-j/2Ψ(α-jn-kβ)

    (2)

    其中:α和β為標(biāo)度因子。令α=2,β=1,則式(2)表示二進(jìn)正交小波基函數(shù),再通過(guò)此式對(duì)x(t)進(jìn)行小波分解,即可得相應(yīng)小波系數(shù):

    W(j,k)={U,di,di-1,…,d1}

    (3)

    其中:U為低頻小波系數(shù);i為小波分解層數(shù);di為第i層高頻小波系數(shù)。對(duì)離散小波的逆變換重構(gòu)定義為:

    (4)

    1.2 軌道不平順狀態(tài)推移模型

    基于軌道不平順隨時(shí)間變化特性及不同時(shí)刻下軌道不平順的概率分布,這里采用一種狀態(tài)推移的模型來(lái)對(duì)軌道不平順變化進(jìn)行描述,模型建立基本流程為:根據(jù)不平順數(shù)據(jù)的范圍,將軌道不平順劃分為n個(gè)離散區(qū)間,這些區(qū)間可視為軌道不平順的不同狀態(tài),則屬于這些狀態(tài)的檢測(cè)點(diǎn)的數(shù)量為S1,S2,···,Sm;再計(jì)算由某時(shí)刻η-1到時(shí)刻η的各狀態(tài)轉(zhuǎn)換的概率,記為Pab,則有以下表達(dá)式:

    (5)

    當(dāng)事件狀態(tài)的概率分布未知時(shí),若樣本的容量較大,則可以通過(guò)對(duì)樣本計(jì)算近似地描述理論的狀態(tài)概率分布。本文利用實(shí)測(cè)的軌道不平順數(shù)據(jù),選取不同時(shí)刻下固定線路里程中的測(cè)量點(diǎn)作為樣本。通過(guò)對(duì)計(jì)算樣本的推移頻率近似得到理論的推移概率。其計(jì)算公式為:

    Pab=aab/ab(a,b=1,2,…,m)

    (6)

    2基于小波分解的軌道不平順預(yù)測(cè)方法

    本文在狀態(tài)推移預(yù)測(cè)模型基礎(chǔ)上,引入小波分解理論,通過(guò)對(duì)原始軌道不平順時(shí)域序列進(jìn)行小波分解及重構(gòu),得到新的不平順時(shí)域序列,并對(duì)其應(yīng)用狀態(tài)推移預(yù)測(cè)模型,從而實(shí)現(xiàn)不同波長(zhǎng)的不平順預(yù)測(cè)功能。其具體步驟如下:

    流程圖如圖1所示。

    圖1 預(yù)測(cè)方法流程圖Fig.1 Flow chart of theprediction method

    3算例驗(yàn)證

    為驗(yàn)證上述方法的正確性,本文通過(guò)選取2個(gè)算例,利用2013年合武線及2011年京廣線軌檢車(chē)實(shí)測(cè)數(shù)據(jù),對(duì)軌道不平順狀態(tài)的變化規(guī)律進(jìn)行研究與分析。

    3.1 算例1

    3.1.1數(shù)據(jù)選擇

    根據(jù)軌檢車(chē)實(shí)測(cè)數(shù)據(jù)的保存形式,本文選取了連續(xù)3次同一位置的合武線上行2 km里程范圍內(nèi)的實(shí)測(cè)數(shù)據(jù)。需說(shuō)明的是:所選里程段均為有砟軌道,設(shè)計(jì)時(shí)速為250 km/h,測(cè)試車(chē)輛為0號(hào)綜合軌檢車(chē),采樣間隔為0.25 m。且通過(guò)計(jì)算可知所選取的樣本均值,方差與據(jù)檢車(chē)全段檢測(cè)數(shù)據(jù)的均值,方差大致相等,因此說(shuō)明本文所選取的樣本容量具有代表性,可以表征線路整體不平順的分布情況。同時(shí),由于預(yù)測(cè)效果與線路維修息息相關(guān),因此本文通過(guò)查閱相關(guān)維修記錄,從而確保所選的3次檢測(cè)時(shí)間范圍內(nèi)進(jìn)行的維修作業(yè)量為最少。此外,由于軌檢車(chē)檢測(cè)對(duì)象包括軌距、左高低、右高低、左軌向、右軌向、水平和三角坑等,對(duì)于不同的檢測(cè)對(duì)象,它們對(duì)軌道質(zhì)量的影響程度也不盡相同。因此針對(duì)不同種類(lèi)的軌道不平順,本文提出的預(yù)測(cè)方法在檢測(cè)不平順狀態(tài)區(qū)間劃分及檢測(cè)周期的選取上都不盡相同。限于篇幅本文僅以高低不平順為例,選擇實(shí)測(cè)的左高低數(shù)據(jù)作為樣本數(shù)據(jù)。

    3.1.2小波基與分解層數(shù)選取

    本文選取具有緊支集及較好正則性的Daubechies小波(dbN)作為小波基函數(shù)。Daubechies小波(dbN)系列隨著序號(hào)N的增加,小波基函數(shù)頻帶劃分效果和時(shí)域波形的光滑性越好,表現(xiàn)出較好的局域性能。但同時(shí)會(huì)使時(shí)域緊支撐性減弱,同時(shí)計(jì)算量增加,實(shí)時(shí)性變差,因此綜合考慮本文選取db6小波作為小波基函數(shù)[10]。

    由于軌檢車(chē)的采樣間隔為0.25 m,且由于軌檢車(chē)檢測(cè)波長(zhǎng)的范圍為1.5~120 m,而小波分解后第i層對(duì)應(yīng)波長(zhǎng)區(qū)間大致為[λs·2i,λs·2i+1](λs為采樣間隔),因此當(dāng)選取小波分解層數(shù)為8時(shí),最后一層的波長(zhǎng)區(qū)間為64, 128 m,已滿(mǎn)足所需波長(zhǎng)范圍的要求,因此本文選用小波分解的層數(shù)為8層。

    3.1.3預(yù)測(cè)結(jié)果分析

    本文選取前兩個(gè)時(shí)刻下的高低不平順序數(shù)據(jù)為參考序列,通過(guò)本文提出的方法對(duì)下一時(shí)刻下的高低不平順序進(jìn)行預(yù)測(cè),并將預(yù)測(cè)結(jié)果與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比。為減小不同時(shí)刻下檢測(cè)數(shù)據(jù)存在的里程偏差,本文采用灰色關(guān)聯(lián)法進(jìn)行里程匹配[11]。通過(guò)計(jì)算分別將前兩次的檢測(cè)數(shù)據(jù)平移55.25和59.5 m作為標(biāo)準(zhǔn)里程。得到新不平順時(shí)域序列后,再對(duì)3組高低不平順序列進(jìn)行小波分解及重構(gòu),選取上述小波基及分解層數(shù),限于篇幅,圖2僅為第1次檢測(cè)數(shù)據(jù)小波重構(gòu)后的軌道不平順預(yù)測(cè)序列。

    圖2 各波長(zhǎng)段小波重構(gòu)后的不平順序列Fig.2 Wavelet reconstructed irregularity of different wavelength

    在得到重構(gòu)不平順序列的基礎(chǔ)上,根據(jù)實(shí)測(cè)數(shù)據(jù)的極差,進(jìn)行高低不平順的狀態(tài)劃分,由于狀態(tài)區(qū)間個(gè)數(shù)的確定與線路實(shí)際養(yǎng)護(hù)維修具體要求有關(guān),通常情況下區(qū)間個(gè)數(shù)越小,即區(qū)間長(zhǎng)度越大,則頻數(shù)分布較為集中,預(yù)測(cè)效果越好。而在實(shí)際運(yùn)用中由于不同里程及波長(zhǎng)段的不平順要求不盡相同,因此狀態(tài)區(qū)間的劃分需按線路的實(shí)際要求來(lái)確定,從而滿(mǎn)足現(xiàn)場(chǎng)實(shí)際需求。所以,在滿(mǎn)足現(xiàn)場(chǎng)預(yù)測(cè)精度的前提下,應(yīng)盡可能地多劃分狀態(tài)區(qū)間,從而達(dá)到更加精確地預(yù)測(cè)效果;通過(guò)前期大量的試算和分析,并結(jié)合相關(guān)文獻(xiàn)[7]的研究結(jié)果,這里選取了15個(gè)狀態(tài),即劃分16個(gè)區(qū)間。此外,由于不同波長(zhǎng)段小波重構(gòu)后軌道不平順的幅值有所變化,所以需確保3組不平順序列進(jìn)行小波重構(gòu)后對(duì)應(yīng)波長(zhǎng)段的區(qū)間劃分一致。

    狀態(tài)區(qū)間劃分后,即可通過(guò)狀態(tài)推移法計(jì)算各區(qū)間頻數(shù)的推移概率Pab,具體方法如下:設(shè)3組不平順序列分別對(duì)應(yīng)時(shí)刻1,2和3所采用的軌檢車(chē)檢測(cè)數(shù)據(jù),先統(tǒng)計(jì)每一次不平順檢測(cè)數(shù)據(jù)各區(qū)間段檢測(cè)點(diǎn)頻數(shù)ab(b=1,2,…,15),再統(tǒng)計(jì)時(shí)刻1檢測(cè)數(shù)據(jù)某一區(qū)間段檢測(cè)點(diǎn)頻數(shù)ab變化到時(shí)刻2各區(qū)間的頻數(shù)aab(a,b=1,2,…,15),即表示原來(lái)滿(mǎn)足時(shí)刻1下的某一狀態(tài)區(qū)間的檢測(cè)點(diǎn)由于軌道惡化從而不平順值發(fā)生了變化,導(dǎo)致這些檢測(cè)點(diǎn)在時(shí)刻2劃分區(qū)間時(shí)被統(tǒng)計(jì)到了其他各個(gè)區(qū)間。統(tǒng)計(jì)完區(qū)間頻數(shù)ab以及區(qū)間變化頻數(shù)aab后,即可通過(guò)式(6)計(jì)算出相鄰兩時(shí)刻下不平順區(qū)間頻數(shù)的推移概率Pab。最后即可在時(shí)刻2檢測(cè)數(shù)據(jù)的基礎(chǔ)上,由式(5)預(yù)測(cè)出時(shí)刻3的不平順區(qū)間頻數(shù)分布。第3次不平順檢測(cè)數(shù)據(jù)各波長(zhǎng)段的預(yù)測(cè)結(jié)果與實(shí)際結(jié)果的對(duì)比如圖3所示:

    由圖3可知,利用本文所提方法得出的預(yù)測(cè)值與軌檢車(chē)的實(shí)際值之間較為符合。對(duì)于預(yù)測(cè)結(jié)果和實(shí)際值之間存在部分偏差,可能是由于線路的不斷維修,導(dǎo)致在部分位置上的檢測(cè)數(shù)據(jù)發(fā)生較大變化,從而產(chǎn)生誤差。

    3.1.4誤差分析

    為驗(yàn)證本文所提方法準(zhǔn)確性,這里采用傳統(tǒng)的狀態(tài)推移法與本文方法進(jìn)行對(duì)比。通過(guò)計(jì)算預(yù)測(cè)值和實(shí)際值之間的相關(guān)系數(shù),分析2種方法預(yù)測(cè)結(jié)果的精度,從而確定方法準(zhǔn)確性,計(jì)算結(jié)果如表1所示。

    表1 預(yù)測(cè)方法評(píng)價(jià)指標(biāo)

    (a)波長(zhǎng)0.5~1 m;(b)波長(zhǎng)1~2 m;(c)波長(zhǎng)2~4 m;(d)波長(zhǎng)4~8 m;(e)波長(zhǎng)8~16 m;(f)波長(zhǎng)16~32 m;(g)波長(zhǎng)32~64 m;(h)波長(zhǎng)64~128 m圖3 高速鐵路第3次檢測(cè)數(shù)據(jù)各波長(zhǎng)段預(yù)測(cè)值與實(shí)際值分布對(duì)比Fig.3 Contrast between the third detection data’s predicted value and actual value of different wavelength in high-speed railway

    由表1可得,本文方法預(yù)測(cè)值與實(shí)際值擬合情況較好,且與其他方法相比,其相關(guān)系數(shù)大致相等,說(shuō)明本文方法與其他方法的預(yù)測(cè)精度相近,可以滿(mǎn)足實(shí)際養(yǎng)護(hù)維修工作預(yù)測(cè)精度要求。

    3.2 算例2

    數(shù)據(jù)來(lái)源為京廣既有提速線上行有砟軌道線路,3次同一里程區(qū)間檢測(cè)數(shù)據(jù),長(zhǎng)度為2 km,軌檢車(chē)類(lèi)型及相關(guān)參數(shù)與算例一相同,同樣針對(duì)軌道高低不平順進(jìn)行分析。通過(guò)計(jì)算分別將前兩次的檢測(cè)數(shù)據(jù)平移62.5和67.75 m作為標(biāo)準(zhǔn)里程。則通過(guò)上述方法得到的第3次不平順檢測(cè)數(shù)據(jù)各波長(zhǎng)段的預(yù)測(cè)結(jié)果與實(shí)際結(jié)果的對(duì)比如圖4所示。

    由圖4可知,利用本文所提方法得出的預(yù)測(cè)值與軌檢車(chē)的實(shí)際值之間較為符合。通過(guò)計(jì)算各波長(zhǎng)段預(yù)測(cè)值與實(shí)際值之間的相關(guān)系數(shù),得到方法的預(yù)測(cè)精度,計(jì)算結(jié)果與算例一相近。結(jié)果表明本文所提方法可以滿(mǎn)足實(shí)際養(yǎng)護(hù)維修工作預(yù)測(cè)精度要求。

    (a)波長(zhǎng)0.5~1 m;(b)波長(zhǎng)1~2 m;(c)波長(zhǎng)2~4 m;(d)波長(zhǎng)4~8 m;(e)波長(zhǎng)8~16 m;(f)波長(zhǎng)16~32 m;(g)波長(zhǎng)32~64 m;(h)波長(zhǎng)64~128 m圖4 既有線第3次檢測(cè)數(shù)據(jù)各波長(zhǎng)段預(yù)測(cè)值與實(shí)際值分布對(duì)比Fig.4 Contrast between the third detection data’s predicted value and actual value of different wavelength in existing speed line

    綜上所述,本文通過(guò)對(duì)2種不同運(yùn)行速度下的軌道不平順檢測(cè)數(shù)據(jù)進(jìn)行方法的驗(yàn)證和分析,其預(yù)測(cè)結(jié)果與實(shí)際值較為符合。因此,說(shuō)明本文提出的方法可以正確預(yù)測(cè)不平順的變化趨勢(shì),且與其他方法相比,本文方法還可以實(shí)現(xiàn)不同波長(zhǎng)段軌道不平順的預(yù)測(cè)。

    4結(jié)論

    1)利用小波分解的預(yù)測(cè)方法,可以實(shí)現(xiàn)軌道不平順有效預(yù)測(cè);且與現(xiàn)有的軌道不平順預(yù)測(cè)方法相比,可以實(shí)現(xiàn)不同波長(zhǎng)段的軌道不平順預(yù)測(cè)。

    2)預(yù)測(cè)所需要的原始數(shù)據(jù)可通過(guò)軌檢車(chē)數(shù)據(jù)處理后直接得到,數(shù)據(jù)源容易獲取,使得該方法簡(jiǎn)單且實(shí)用。

    3)該方法可以通過(guò)各狀態(tài)的變化推移計(jì)算出軌道不平順由某一范圍變化到另一范圍的數(shù)量及概率,可以反映出在不同時(shí)刻下軌道整體質(zhì)量的變化趨勢(shì),在評(píng)估軌道質(zhì)量,制定養(yǎng)護(hù)維修計(jì)劃等方面都有一定應(yīng)用價(jià)值。

    4)由于線路的不斷維修,在部分位置上的檢測(cè)數(shù)據(jù)可能發(fā)生較大變化,因此還需進(jìn)一步完善該方法,對(duì)其進(jìn)行相關(guān)處理,進(jìn)而實(shí)現(xiàn)更為精準(zhǔn)的預(yù)測(cè)。

    參考文獻(xiàn):

    [1] 佐藤吉彥.日本高速鐵路的全壽命維護(hù)[J].中國(guó)鐵道科學(xué), 2001, 22(1): 6-14.

    Misaki Sato Yoshihiko.The life cycle ma-intenance of high speed railways [J].China Railway Science, 2001, 22 (1): 6-14.

    [2] Hideo Miyamoto.日本東海道新干線三十年軌道維修經(jīng)驗(yàn)[J].上海鐵道科技, 1998(2):43-45.

    Hideo Miyamoto.The thirty years of trackmaintenance experience of Japanese Tokai-do Shinkansen[J].Shanghai Railway Science & Technology, 1998(2): 43-45.

    [3] 內(nèi)田雅夫.高速鐵路軌道的維修管理[J].鐵道建筑, 1992(10): 6-9.

    Uchida Masao.The high speed railway track repair management[J].Railway Engineering, 1992(10): 6-9.

    [4] 許玉德, 吳紀(jì)才.利用線性預(yù)測(cè)模型分析軌道不平順發(fā)展[J].石家莊鐵道學(xué)院學(xué)報(bào), 2005, 18(1): 6-9.

    XU Yude, WU Jicai.Analysis on development of track irregularities with linear forecast Model[J].Journal of Shijiazhuang Railway Institute, 2005, 18 (1): 6-9.

    [5] 陳憲麥.軌道不平順時(shí)頻域分析及預(yù)測(cè)方法的研究[D].北京:鐵道科學(xué)研究院, 2006: 99-101.

    CHEN Xianmai.Time frequence analysis and prediction study on track irregularities[D].Beijing: Chinese Railway Science Research Institute, 2006: 99-101.

    [6] 曲建軍, 高亮, 張新奎, 等.基于灰色GM(1,1)非等時(shí)距修正模型的軌道質(zhì)量預(yù)測(cè)[J].中國(guó)鐵道科學(xué), 2009, 30(3): 5-8.

    QU Jianjun, GAO Liang, ZHANG Xinkui,et al.A prediction method for track quality based on grey GM(1,1) nonequal time interval correction model[J].China Railway Science.2009, 30(3): 5-8.

    [7] 高建敏.基于狀態(tài)轉(zhuǎn)移概率矩陣的軌道不平順發(fā)展預(yù)測(cè)研究[J].鐵路建筑, 2011(7): 140-143.

    GAO Jianmin.The development of the track irregularity prediction research based on The state transition probability matrix[J].Railway Engineering, 2011(7): 140-143.

    [8] 陳憲麥, 王瀾, 陶夏新, 等.基于小波分析理論的軌道不平順?lè)治鯷J].鐵道工程學(xué)報(bào), 2008, 1(1): 57-61.

    CHEN Xianmai, WANG Lian,TAO Xiaxin,et al.The analysis of track irregularity based on the theory of wavelet analysis[J].Journal of Railway Engineering Society, 2008, 1(1): 57-61.

    [9] 徐磊, 陳憲麥, 徐偉昌, 等.基于小波和Wigner-Ville分布的軌道不平順特征識(shí)別[J].中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 44(8): 3344-3350.

    XU Lei, CHEN Xianmai, XU Weichang,et al.Explored of track irregularity’s characteristic identification based on wavelet method and WignerVille distribution[J].Journalof Central South University(Science and Technology), 2013, 44(8): 3344-3350.

    [10] 曲銘, 許玉德.基于小波的列車(chē)加速度和軌道不平順關(guān)系分析[J].華東交通大學(xué)學(xué)報(bào), 2008(5): 36-41.

    QU Ming, XU Yude.An analysis of the relationship of acceleration and track irregularity based on wavelet[J].Journal of East China Jiaotong University, 2008(5): 36-41.

    [11] 李再幃, 練松良.改進(jìn)HHT方法在軌道不平順信號(hào)分析中的應(yīng)用[J].深圳大學(xué)學(xué)報(bào)(理工版), 2012, 29(3): 270-275.

    LI Zaiwei, LIAN Songliang.Analysisof track irregularity signal by improved Hilbert Huang Transform[J].Journal of Shenzhen University Science and Engineering, 2012, 29(3): 270-275.

    (編輯蔣學(xué)東)

    Theprediction method considering the factors of track irregularity wavelength

    Lü Hong1,LI Zaiwei1,2,3, HE Yuelei1

    (1.School of Urban Rail Transportation, Shanghai University of Engineering Science, Shanghai 201620, China;

    2.Engineering Research Center of Railway Environmental Vibration and Noise of Ministry

    of Education, East China Jiaotong University, Nanchang 330013, China;

    3.School of Civil Engineering and Architecture, Beijing Jiaotong University, Beijing 100044, China)

    Abstract:The accurate prediction of the change law of track irregularity is one of the key problems to improve the efficiency of track maintenance and reduce maintenance cost.According to the change trend of the track irregularity of different wavelengths and the influence on operation, and basing on the theory of wavelet decomposition, this paper selected the appropriate wavelet basis function and decomposition levels to complete the wavelet decomposition and its reconstruction for the original irregularity time domain sequence.By using the state shifting method to analyze the irregularity sequence after reconstruction, and making use of the different times of orbit the irregularity of the state distribution, the track irregularity was predicted.The feasibility of this mothod was verified by a practical example and the error between the outual results and the prediction results was analyzed.The results show that the error between the predicted results and the actual results using the method proposed in this paper is small.It can effectively predict the change tendency of track irregularity, and this method can realize the prediction of track irregularity of different wavelengths.The method can provide technical support for the research of track irregularity prediction in our country.

    Key words:wavelet theory; track irregularity; prediction method; state shifting method

    通訊作者:李再幃(1983-),男,吉林大安人,副教授,博士,從事軌道結(jié)構(gòu)振動(dòng)研究;E-mail:lzw5220964@163.com

    基金項(xiàng)目:上海市科學(xué)技術(shù)委員會(huì)地方院校能力建設(shè)資助項(xiàng)目(14110501300);上海市研究生教育創(chuàng)新計(jì)劃學(xué)位點(diǎn)引導(dǎo)布局與建設(shè)培育資助項(xiàng)目(13sc002);上海市教育委員會(huì)科研創(chuàng)新資助項(xiàng)目(14YZ137);上海工程技術(shù)大學(xué)研究生科研創(chuàng)新資助項(xiàng)目(14KY1005)

    收稿日期:2015-04-09

    中圖分類(lèi)號(hào):U216

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1672-7029(2015)06-1312-07

    猜你喜歡
    預(yù)測(cè)方法
    電網(wǎng)規(guī)劃負(fù)荷預(yù)測(cè)方法探討
    區(qū)域經(jīng)濟(jì)知識(shí)需求演化預(yù)測(cè)方法探討
    風(fēng)電短期發(fā)電功率預(yù)測(cè)方法探討
    臨汾區(qū)塊煤系地層天然氣水合物危害與預(yù)防措施
    水驅(qū)油效率影響因素研究進(jìn)展
    預(yù)測(cè)策略教學(xué)技巧在初中英語(yǔ)閱讀課的實(shí)踐探索
    船舶交通流量及幾種預(yù)測(cè)方法分析
    科技視界(2016年12期)2016-05-25 08:58:42
    油氣管道腐蝕剩余壽命預(yù)測(cè)方法的研究進(jìn)展
    再論科技期刊出版時(shí)滯
    湖北省49.5MW風(fēng)電場(chǎng)水土流失分析及保持方案研究
    綠色科技(2014年1期)2014-11-28 23:55:13
    国产av麻豆久久久久久久| 中国美女看黄片| av欧美777| 午夜福利高清视频| 精品国产超薄肉色丝袜足j| 亚洲 国产 在线| 久久人人精品亚洲av| 精品人妻1区二区| 亚洲精品美女久久久久99蜜臀| 久久精品91无色码中文字幕| 熟女少妇亚洲综合色aaa.| 一本大道久久a久久精品| 亚洲av成人不卡在线观看播放网| 久久久国产精品麻豆| 俺也久久电影网| 国产片内射在线| 亚洲色图av天堂| 18禁观看日本| 天天躁狠狠躁夜夜躁狠狠躁| 色在线成人网| 久久午夜综合久久蜜桃| 日韩欧美在线乱码| 成人高潮视频无遮挡免费网站| 亚洲va日本ⅴa欧美va伊人久久| 变态另类成人亚洲欧美熟女| 不卡一级毛片| 久久欧美精品欧美久久欧美| 一区福利在线观看| 十八禁网站免费在线| 欧美zozozo另类| 国产视频内射| 午夜久久久久精精品| 国产精品亚洲美女久久久| 在线观看66精品国产| 亚洲精品中文字幕在线视频| 一边摸一边做爽爽视频免费| 男人舔女人的私密视频| 国产一级毛片七仙女欲春2| 女同久久另类99精品国产91| 免费在线观看影片大全网站| 99国产极品粉嫩在线观看| 啦啦啦免费观看视频1| 国产蜜桃级精品一区二区三区| 五月玫瑰六月丁香| 叶爱在线成人免费视频播放| 18美女黄网站色大片免费观看| 精品熟女少妇八av免费久了| 国产成人欧美在线观看| 久久人妻av系列| 亚洲国产欧美网| 又粗又爽又猛毛片免费看| 亚洲国产欧美人成| 国产亚洲欧美98| 日日摸夜夜添夜夜添小说| 999精品在线视频| 久久亚洲精品不卡| 国产黄色小视频在线观看| 国产成人一区二区三区免费视频网站| 中文字幕高清在线视频| 村上凉子中文字幕在线| 国产精品精品国产色婷婷| 黄色片一级片一级黄色片| 91大片在线观看| 国产亚洲av高清不卡| 国产又色又爽无遮挡免费看| a级毛片在线看网站| 校园春色视频在线观看| 757午夜福利合集在线观看| 亚洲真实伦在线观看| 琪琪午夜伦伦电影理论片6080| 高清毛片免费观看视频网站| 50天的宝宝边吃奶边哭怎么回事| 亚洲全国av大片| 久久香蕉国产精品| 日本一区二区免费在线视频| 亚洲av成人一区二区三| 精品人妻1区二区| 欧美性猛交黑人性爽| 亚洲五月婷婷丁香| 91老司机精品| 一区二区三区高清视频在线| 国产精品自产拍在线观看55亚洲| 非洲黑人性xxxx精品又粗又长| 日日夜夜操网爽| 色哟哟哟哟哟哟| 亚洲精品色激情综合| 青草久久国产| 日日爽夜夜爽网站| 99国产极品粉嫩在线观看| 长腿黑丝高跟| 免费在线观看黄色视频的| 亚洲在线自拍视频| 一二三四在线观看免费中文在| 美女大奶头视频| 亚洲九九香蕉| 可以免费在线观看a视频的电影网站| 亚洲一区二区三区不卡视频| 国产精品香港三级国产av潘金莲| 国产成人影院久久av| 亚洲欧洲精品一区二区精品久久久| 欧美精品啪啪一区二区三区| 国产精品久久久久久精品电影| 母亲3免费完整高清在线观看| 国产成人欧美在线观看| 国产精品99久久99久久久不卡| 男女下面进入的视频免费午夜| 制服丝袜大香蕉在线| 国产69精品久久久久777片 | 九色国产91popny在线| 婷婷精品国产亚洲av| 亚洲精品久久成人aⅴ小说| 亚洲国产精品成人综合色| 欧美中文综合在线视频| 亚洲 欧美 日韩 在线 免费| 欧美一区二区国产精品久久精品 | 麻豆一二三区av精品| 久久久久亚洲av毛片大全| 欧美性猛交╳xxx乱大交人| 久久草成人影院| 亚洲av五月六月丁香网| 国产黄色小视频在线观看| 校园春色视频在线观看| 少妇粗大呻吟视频| 黄色视频,在线免费观看| 亚洲欧美激情综合另类| 亚洲色图av天堂| 免费在线观看完整版高清| 精品欧美国产一区二区三| 99国产精品99久久久久| 美女午夜性视频免费| 亚洲aⅴ乱码一区二区在线播放 | 欧美日韩乱码在线| 两个人看的免费小视频| 日韩免费av在线播放| or卡值多少钱| 午夜福利成人在线免费观看| 麻豆久久精品国产亚洲av| 淫秽高清视频在线观看| 午夜福利成人在线免费观看| 久久久久久九九精品二区国产 | 老司机在亚洲福利影院| 五月玫瑰六月丁香| 国产区一区二久久| 免费无遮挡裸体视频| 男女床上黄色一级片免费看| 少妇的丰满在线观看| 久热爱精品视频在线9| 亚洲熟妇中文字幕五十中出| 精品午夜福利视频在线观看一区| 日韩欧美 国产精品| 舔av片在线| 亚洲国产中文字幕在线视频| 人妻夜夜爽99麻豆av| 国产精品自产拍在线观看55亚洲| 国产又黄又爽又无遮挡在线| 久久国产乱子伦精品免费另类| 国产精品久久视频播放| 搡老熟女国产l中国老女人| 国产又色又爽无遮挡免费看| 操出白浆在线播放| 亚洲一区中文字幕在线| 美女高潮喷水抽搐中文字幕| 欧美一级毛片孕妇| 日韩欧美 国产精品| 国产精品,欧美在线| 久久久久久久精品吃奶| 欧美中文综合在线视频| 日韩中文字幕欧美一区二区| 两个人免费观看高清视频| 亚洲人成77777在线视频| 给我免费播放毛片高清在线观看| 高清在线国产一区| 日韩国内少妇激情av| 亚洲全国av大片| 99在线人妻在线中文字幕| 午夜福利视频1000在线观看| svipshipincom国产片| 757午夜福利合集在线观看| 无遮挡黄片免费观看| 最近视频中文字幕2019在线8| 亚洲人与动物交配视频| 精品国产乱子伦一区二区三区| 动漫黄色视频在线观看| 99久久精品热视频| 国产成人影院久久av| 国产精品一区二区免费欧美| 免费在线观看影片大全网站| 九色成人免费人妻av| 1024手机看黄色片| 亚洲av五月六月丁香网| 99国产综合亚洲精品| 亚洲精品国产一区二区精华液| 一边摸一边做爽爽视频免费| 美女黄网站色视频| 国产精品爽爽va在线观看网站| 午夜免费成人在线视频| 日韩欧美精品v在线| 免费高清视频大片| www日本黄色视频网| 国产精品精品国产色婷婷| 一级作爱视频免费观看| 免费在线观看日本一区| 好看av亚洲va欧美ⅴa在| 人妻夜夜爽99麻豆av| 90打野战视频偷拍视频| 精品第一国产精品| 日韩 欧美 亚洲 中文字幕| 亚洲欧美日韩高清专用| 老司机在亚洲福利影院| 一个人免费在线观看的高清视频| 国产在线精品亚洲第一网站| 999久久久精品免费观看国产| 亚洲免费av在线视频| 久久久国产精品麻豆| 99精品久久久久人妻精品| 天堂av国产一区二区熟女人妻 | 久久精品国产亚洲av香蕉五月| 免费搜索国产男女视频| 亚洲精品国产一区二区精华液| 正在播放国产对白刺激| 亚洲真实伦在线观看| 国产伦人伦偷精品视频| 黄色毛片三级朝国网站| 91麻豆精品激情在线观看国产| x7x7x7水蜜桃| 日本成人三级电影网站| 亚洲九九香蕉| 麻豆一二三区av精品| 搡老岳熟女国产| 麻豆成人午夜福利视频| 一级黄色大片毛片| 久久中文字幕一级| 精品久久久久久久人妻蜜臀av| 在线观看一区二区三区| 国产97色在线日韩免费| 黄色 视频免费看| 99国产精品一区二区蜜桃av| 中文资源天堂在线| 久久香蕉国产精品| 夜夜躁狠狠躁天天躁| 五月伊人婷婷丁香| 日本撒尿小便嘘嘘汇集6| 亚洲熟女毛片儿| 在线观看美女被高潮喷水网站 | 国产野战对白在线观看| 在线观看日韩欧美| 亚洲精品在线美女| 久久久久久久精品吃奶| 亚洲一区二区三区不卡视频| 视频区欧美日本亚洲| 欧美极品一区二区三区四区| 亚洲真实伦在线观看| 少妇裸体淫交视频免费看高清 | 欧美高清成人免费视频www| 精品国产超薄肉色丝袜足j| 美女黄网站色视频| 亚洲五月婷婷丁香| 亚洲av成人精品一区久久| 亚洲中文字幕一区二区三区有码在线看 | 国产69精品久久久久777片 | 啦啦啦观看免费观看视频高清| 无限看片的www在线观看| 国产精品综合久久久久久久免费| 亚洲精品粉嫩美女一区| 欧美成狂野欧美在线观看| 中文字幕高清在线视频| 国产精品久久久久久亚洲av鲁大| 变态另类丝袜制服| 18禁黄网站禁片午夜丰满| 久久久久免费精品人妻一区二区| 99久久久亚洲精品蜜臀av| 亚洲男人天堂网一区| 午夜免费激情av| 天堂影院成人在线观看| 国内少妇人妻偷人精品xxx网站 | 男人舔奶头视频| av有码第一页| 毛片女人毛片| 亚洲人成网站高清观看| 我的老师免费观看完整版| 美女大奶头视频| 久久中文字幕一级| 亚洲人与动物交配视频| 精品国内亚洲2022精品成人| 精品一区二区三区视频在线观看免费| 国产av不卡久久| 中文字幕精品亚洲无线码一区| 小说图片视频综合网站| 男女视频在线观看网站免费 | 嫩草影院精品99| 欧美另类亚洲清纯唯美| 在线观看免费视频日本深夜| 国产精品一区二区三区四区久久| 后天国语完整版免费观看| 久久国产精品人妻蜜桃| 亚洲国产高清在线一区二区三| 亚洲激情在线av| 十八禁人妻一区二区| 亚洲人成网站高清观看| 妹子高潮喷水视频| 欧美丝袜亚洲另类 | 舔av片在线| 免费在线观看亚洲国产| 黑人操中国人逼视频| 好男人电影高清在线观看| 欧美日韩亚洲综合一区二区三区_| 老司机午夜福利在线观看视频| 老司机在亚洲福利影院| 岛国视频午夜一区免费看| tocl精华| 亚洲av成人一区二区三| 999久久久国产精品视频| 黄色视频不卡| 国产单亲对白刺激| 亚洲成av人片在线播放无| 国产不卡一卡二| 精品国产乱码久久久久久男人| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 两个人免费观看高清视频| 少妇的丰满在线观看| av天堂在线播放| 欧美日韩福利视频一区二区| 日韩欧美 国产精品| 国产亚洲av嫩草精品影院| 小说图片视频综合网站| 99久久无色码亚洲精品果冻| 91老司机精品| 免费在线观看影片大全网站| 最近视频中文字幕2019在线8| www国产在线视频色| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| xxxwww97欧美| 国产精品乱码一区二三区的特点| 啦啦啦免费观看视频1| 久久精品91无色码中文字幕| 夜夜躁狠狠躁天天躁| 神马国产精品三级电影在线观看 | 欧美高清成人免费视频www| 久久人人精品亚洲av| 黄色视频,在线免费观看| 欧美3d第一页| 国产精品综合久久久久久久免费| 99精品久久久久人妻精品| 精品一区二区三区视频在线观看免费| 一本一本综合久久| 神马国产精品三级电影在线观看 | 国产精品久久久av美女十八| 又爽又黄无遮挡网站| 国产精品影院久久| 国内精品久久久久久久电影| 亚洲美女黄片视频| 99久久精品国产亚洲精品| 欧美又色又爽又黄视频| 国产成人一区二区三区免费视频网站| av免费在线观看网站| 19禁男女啪啪无遮挡网站| 国产成人系列免费观看| 久久亚洲真实| 久久久久久久久中文| 巨乳人妻的诱惑在线观看| 午夜精品久久久久久毛片777| 亚洲va日本ⅴa欧美va伊人久久| 深夜精品福利| 国产亚洲欧美98| 亚洲中文字幕一区二区三区有码在线看 | 亚洲欧美激情综合另类| 日本黄色视频三级网站网址| 国产成人精品无人区| 看片在线看免费视频| 精华霜和精华液先用哪个| 一级作爱视频免费观看| 国产乱人伦免费视频| 黄色视频不卡| 欧美国产日韩亚洲一区| 午夜福利在线在线| 日韩大码丰满熟妇| 黑人操中国人逼视频| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 国产精品一及| 日韩大尺度精品在线看网址| 九九热线精品视视频播放| 亚洲成av人片在线播放无| 俄罗斯特黄特色一大片| 欧美激情久久久久久爽电影| 麻豆av在线久日| 精品少妇一区二区三区视频日本电影| 午夜精品在线福利| 亚洲成人久久性| 日本一区二区免费在线视频| 国产野战对白在线观看| 亚洲精品国产一区二区精华液| 婷婷六月久久综合丁香| 99热这里只有是精品50| 琪琪午夜伦伦电影理论片6080| 久久久久久久久中文| av在线天堂中文字幕| 免费在线观看视频国产中文字幕亚洲| 免费人成视频x8x8入口观看| 麻豆成人av在线观看| 国产精华一区二区三区| 久久久久久九九精品二区国产 | 亚洲国产看品久久| 国产精品永久免费网站| 亚洲av美国av| 国产1区2区3区精品| 最近视频中文字幕2019在线8| 久久精品亚洲精品国产色婷小说| 69av精品久久久久久| 国语自产精品视频在线第100页| 中国美女看黄片| 欧美三级亚洲精品| 亚洲成人久久爱视频| av中文乱码字幕在线| 狂野欧美激情性xxxx| av在线播放免费不卡| 久久伊人香网站| 99热只有精品国产| 中文字幕高清在线视频| 久久精品国产亚洲av高清一级| 色老头精品视频在线观看| 欧美黑人巨大hd| 一边摸一边做爽爽视频免费| 久久香蕉激情| 2021天堂中文幕一二区在线观| 狂野欧美激情性xxxx| 欧美黑人精品巨大| 制服人妻中文乱码| 国产亚洲欧美在线一区二区| 一本综合久久免费| 男女视频在线观看网站免费 | 在线看三级毛片| 日韩成人在线观看一区二区三区| 久久伊人香网站| 99精品欧美一区二区三区四区| 看片在线看免费视频| 久久伊人香网站| 欧美午夜高清在线| 精品国产乱子伦一区二区三区| 哪里可以看免费的av片| 看片在线看免费视频| 精品国产乱码久久久久久男人| 亚洲午夜精品一区,二区,三区| 日韩大码丰满熟妇| 久久久久久九九精品二区国产 | 久久精品国产清高在天天线| 免费在线观看日本一区| 最近视频中文字幕2019在线8| 成人高潮视频无遮挡免费网站| 亚洲精品久久国产高清桃花| 12—13女人毛片做爰片一| 精品无人区乱码1区二区| 亚洲成人中文字幕在线播放| 亚洲欧洲精品一区二区精品久久久| 丰满人妻一区二区三区视频av | 国内精品久久久久久久电影| www日本黄色视频网| 欧美大码av| 亚洲国产欧美人成| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利在线观看吧| 日本黄大片高清| 国产成人av教育| 亚洲国产精品999在线| 国产片内射在线| 19禁男女啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品成人综合色| 亚洲成人免费电影在线观看| 国产成人av教育| 啦啦啦免费观看视频1| 成人18禁高潮啪啪吃奶动态图| 女生性感内裤真人,穿戴方法视频| 国产午夜精品论理片| 中亚洲国语对白在线视频| 色哟哟哟哟哟哟| 欧美色欧美亚洲另类二区| 国内揄拍国产精品人妻在线| 人妻丰满熟妇av一区二区三区| 极品教师在线免费播放| 国产精品99久久99久久久不卡| 制服诱惑二区| 最好的美女福利视频网| www.熟女人妻精品国产| 变态另类丝袜制服| 欧美 亚洲 国产 日韩一| 久久久久国内视频| 97碰自拍视频| 欧美日韩亚洲国产一区二区在线观看| 99精品欧美一区二区三区四区| 韩国av一区二区三区四区| 日本 欧美在线| 国产爱豆传媒在线观看 | 国产成年人精品一区二区| 91在线观看av| 99久久无色码亚洲精品果冻| 成人av一区二区三区在线看| 国产单亲对白刺激| 精品乱码久久久久久99久播| 国产成人精品久久二区二区免费| 日韩av在线大香蕉| 我的老师免费观看完整版| 国产成人欧美在线观看| 免费在线观看完整版高清| 欧美日本视频| 精品久久久久久久末码| 又粗又爽又猛毛片免费看| 欧美极品一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 国产亚洲精品一区二区www| 久久精品成人免费网站| 91字幕亚洲| 精品久久久久久成人av| 久久久水蜜桃国产精品网| 日本熟妇午夜| 国产欧美日韩一区二区精品| 国产精品久久久av美女十八| 淫妇啪啪啪对白视频| 制服丝袜大香蕉在线| 欧美三级亚洲精品| 欧美中文日本在线观看视频| 欧美黑人欧美精品刺激| 搞女人的毛片| 给我免费播放毛片高清在线观看| 日韩成人在线观看一区二区三区| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 动漫黄色视频在线观看| 国产97色在线日韩免费| 宅男免费午夜| 在线观看舔阴道视频| 18禁裸乳无遮挡免费网站照片| 此物有八面人人有两片| 色尼玛亚洲综合影院| 亚洲av日韩精品久久久久久密| 婷婷亚洲欧美| 在线观看免费午夜福利视频| 欧美黑人精品巨大| 男男h啪啪无遮挡| 午夜福利在线观看吧| 午夜福利在线在线| 国产aⅴ精品一区二区三区波| 极品教师在线免费播放| 婷婷丁香在线五月| 九色成人免费人妻av| 激情在线观看视频在线高清| av福利片在线| 欧美黑人精品巨大| 亚洲成av人片在线播放无| 午夜免费激情av| 亚洲国产精品sss在线观看| 夜夜躁狠狠躁天天躁| 国产高清有码在线观看视频 | 日韩高清综合在线| 久久久国产欧美日韩av| 看免费av毛片| 国产精品 国内视频| 亚洲精品在线观看二区| 法律面前人人平等表现在哪些方面| 亚洲国产欧美人成| 日本成人三级电影网站| 欧美黄色淫秽网站| 天堂√8在线中文| 免费无遮挡裸体视频| 久久这里只有精品19| 宅男免费午夜| 亚洲欧美日韩东京热| 国产精品美女特级片免费视频播放器 | 国产一区在线观看成人免费| 久久 成人 亚洲| 国产精品久久视频播放| 久久人人精品亚洲av| 精品国产乱子伦一区二区三区| av视频在线观看入口| 妹子高潮喷水视频| 精品久久久久久成人av| 舔av片在线| 欧美日韩黄片免| 91麻豆av在线| 久久久久久大精品| 亚洲熟女毛片儿| av片东京热男人的天堂| 99精品欧美一区二区三区四区| 色老头精品视频在线观看| 中文字幕精品亚洲无线码一区| 亚洲,欧美精品.| 久久久久国产精品人妻aⅴ院| 女警被强在线播放| 91麻豆精品激情在线观看国产| 免费在线观看视频国产中文字幕亚洲| 在线观看日韩欧美| 香蕉av资源在线| 欧美 亚洲 国产 日韩一| 欧美日韩福利视频一区二区| 91麻豆精品激情在线观看国产| 欧美性长视频在线观看| 欧美日韩福利视频一区二区| 亚洲激情在线av| 日本熟妇午夜| 精品一区二区三区四区五区乱码| 久久中文看片网| 成人国产综合亚洲| 97超级碰碰碰精品色视频在线观看| 日本精品一区二区三区蜜桃| 成年版毛片免费区| 国产97色在线日韩免费| 亚洲精品国产精品久久久不卡| 欧美成人午夜精品| 亚洲欧美精品综合久久99| 亚洲av电影在线进入| 妹子高潮喷水视频| 欧美国产日韩亚洲一区| 国产黄色小视频在线观看|