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

    基于時(shí)間序列分析的橋梁長期撓度分離與預(yù)測

    2016-07-22 07:40:44陳國良林訓(xùn)根
    關(guān)鍵詞:ARIMA模型分離撓度

    陳國良, 林訓(xùn)根, 岳 青, 劉 華

    (1. 中國礦業(yè)大學(xué) 環(huán)境與測繪學(xué)院, 江蘇 徐州 221116; 2. 中鐵大橋勘測設(shè)計(jì)院集團(tuán)有限公司, 湖北 武漢 430050)

    ?

    基于時(shí)間序列分析的橋梁長期撓度分離與預(yù)測

    陳國良1, 林訓(xùn)根1, 岳青2, 劉華2

    (1. 中國礦業(yè)大學(xué) 環(huán)境與測繪學(xué)院, 江蘇 徐州 221116; 2. 中鐵大橋勘測設(shè)計(jì)院集團(tuán)有限公司, 湖北 武漢 430050)

    摘要:針對撓度長期監(jiān)測信號的多尺度特性及溫度荷載下?lián)隙戎芷谛宰兓奶攸c(diǎn),利用修正平均周期圖法獲取信號分量周期長度,結(jié)合基于中心移動(dòng)平均法的撓度溫度效應(yīng)分離策略,提取結(jié)構(gòu)變形特征,運(yùn)用ARIMA模型建立結(jié)構(gòu)變形趨勢的預(yù)測模型.以黃岡長江大橋?yàn)槔懻摿藴囟扰c撓度的相關(guān)性,通過模擬數(shù)值試驗(yàn)驗(yàn)證分離算法的有效性.實(shí)測數(shù)據(jù)的處理與分析結(jié)果表明,本文方法具有較高的可靠性和精度,可為橋梁進(jìn)行損傷識別提供數(shù)據(jù)支持.

    關(guān)鍵詞:撓度; 分離; 長期趨勢; ARIMA模型; 預(yù)測

    撓度是橋梁結(jié)構(gòu)健康監(jiān)測中最能直觀反應(yīng)結(jié)構(gòu)安全狀態(tài)的參數(shù),直接反映了結(jié)構(gòu)的整體豎向剛度、承載能力.運(yùn)營期監(jiān)測的橋梁撓度信號是環(huán)境荷載、車輛荷載、材料劣化、噪聲等因素共同作用下的綜合響應(yīng)[1-2].準(zhǔn)確分離橋梁在各個(gè)因素作用下的撓度特性、預(yù)測結(jié)構(gòu)長期變形趨勢,有助于橋梁工作性能評估和結(jié)構(gòu)安全預(yù)警.結(jié)構(gòu)運(yùn)行狀態(tài)與環(huán)境參數(shù)密切相關(guān),因而針對撓度信號的分離與預(yù)測是近年來研究的熱點(diǎn)之一.文獻(xiàn)[3-5]在分析長期監(jiān)測數(shù)據(jù)的基礎(chǔ)上得出撓度與溫度并非一般的線性關(guān)系,采用主成分分析法討論溫度的作用規(guī)律,卻未從監(jiān)測數(shù)據(jù)中提取溫度影響成分,缺乏定量評價(jià);Tang[6]應(yīng)用經(jīng)驗(yàn)?zāi)B(tài)分解對信號進(jìn)行通道擴(kuò)展構(gòu)建基于獨(dú)立分量分析的撓度分離模型,但在實(shí)際應(yīng)用中區(qū)別特征信號與虛假信號的固有模態(tài)分量存在一定難度;梁宗保[7]通過建立相位差為半個(gè)日周期的撓度-溫度效應(yīng)回歸方程,從而剔除結(jié)構(gòu)信息中的溫度效應(yīng),但對于長期監(jiān)測而言,相位差并非均勻的半個(gè)日周期,因而無法確?;貧w方程的準(zhǔn)確性;劉夏平等[8]研究了基于奇異值分解的盲信號分離方法,然而其忽略了構(gòu)建分解矩陣時(shí)對分離信號周期長度的討論;文獻(xiàn)[9-11]分別提出了自適應(yīng)帶寬濾波法、最小二乘支持向量機(jī)、多最小二乘支持向量機(jī)進(jìn)行撓度-溫度效應(yīng)分離的方法,但算法較為復(fù)雜,實(shí)用性不足.

    本文利用撓度各響應(yīng)分量在時(shí)間尺度上不耦合的特點(diǎn),提出了基于時(shí)間序列分析的撓度分離與預(yù)測方法,實(shí)現(xiàn)撓度-溫度效應(yīng)的有效分離,從而定量評價(jià)結(jié)構(gòu)變形特征,再以求和自回歸滑動(dòng)平均(autoregressive integrated and moving average,ARIMA)模型進(jìn)行精確擬合,準(zhǔn)確預(yù)測橋梁變形趨勢,為橋梁撓度-溫度效應(yīng)的分離及長期趨勢的預(yù)測提供了一種思路.

    1撓度響應(yīng)分析

    1.1撓度響應(yīng)加法模型

    隨著橋梁服役期限的增加,預(yù)應(yīng)力逐漸損傷、混凝土收縮徐變、溫度變化等原因,主跨跨中均存在不同程度的下?lián)?,影響行車安全,因而長期撓度變形是橋梁結(jié)構(gòu)監(jiān)測的重要一環(huán).健康監(jiān)測系統(tǒng)獲得的撓度數(shù)據(jù)包含了溫度荷載、車輛荷載、結(jié)構(gòu)劣化效應(yīng)等信息.溫度效應(yīng)是結(jié)構(gòu)在外界環(huán)境影響下產(chǎn)生溫度變化而引起結(jié)構(gòu)不受構(gòu)件約束的均勻變形,一般認(rèn)為自然環(huán)境下的溫度荷載(Tt,t為時(shí)間.)包括日溫差荷載(Td)、年溫差荷載(Ty)及驟然升降溫荷載(Tu).車輛荷載(Vt)屬于隨機(jī)影響因素,大跨度橋梁結(jié)構(gòu)趨于柔性,對于鐵路橋梁而言,列車荷載造成的瞬時(shí)下?lián)陷^為突出,文獻(xiàn)[7]提出利用小波理論,根據(jù)撓度信號中溫度效應(yīng)與活載效應(yīng)時(shí)頻特征的不同,將監(jiān)測信號進(jìn)行多尺度分解,提取實(shí)際測量信息中的活載效應(yīng).撓度的變化趨勢(Lt)具有長期不穩(wěn)定的特點(diǎn),通常只能對變形的影響程度做出定性分析,而準(zhǔn)確預(yù)測結(jié)構(gòu)長期變形存在一定難度.

    結(jié)構(gòu)響應(yīng)的各個(gè)影響因素具有多尺度的特征,例如溫度效應(yīng)的周期包含日周期和年周期等,而結(jié)構(gòu)的長期變形有可能長達(dá)數(shù)月甚至數(shù)年,通常視各類效應(yīng)在時(shí)間尺度上并不耦合,即某一荷載所產(chǎn)生的變形不受同時(shí)作用的其他荷載的影響,因此選擇加法模型[8]作為撓度響應(yīng)時(shí)間序列Xt的表達(dá).

    (1)

    式中:Nt為環(huán)境噪聲.

    1.2分離策略

    撓度分離策略采用逐一剝離的方法.首先,利用小波分析剔除活載效應(yīng)與環(huán)境噪聲,再利用時(shí)間序列分析的中心移動(dòng)平均法分離日溫差效應(yīng)和年溫差效應(yīng),最后獲得撓度變形的長期趨勢項(xiàng).

    設(shè)撓度響應(yīng)時(shí)間序列觀測值為{Xt,t=1,2,…,n},其中n為觀測值總數(shù).在準(zhǔn)確確定周期長度(T)的情況下,移動(dòng)平均法可以在很大程度上削弱或消除原有序列中的溫度效應(yīng)周期項(xiàng)[12],使時(shí)間序列趨于平滑,從而得到監(jiān)測信號中的長期趨勢特征.環(huán)境噪聲與車輛荷載都屬于高頻信號,采用修正平均周期圖法(Welch法)對經(jīng)過小波濾波后的撓度信號進(jìn)行頻譜分析可以準(zhǔn)確獲得信號的頻域分布信息,轉(zhuǎn)化為周期長度,作為中心移動(dòng)平均的移動(dòng)步長參數(shù).以序列起點(diǎn)進(jìn)行移動(dòng)平均會(huì)導(dǎo)致分離結(jié)果存在滯后性,因此本文提出對移動(dòng)平均的時(shí)間相位進(jìn)行半周期延遲,即以周期中點(diǎn)為起點(diǎn)進(jìn)行移動(dòng)平均,以規(guī)避這種缺陷,時(shí)域相位延遲長度需根據(jù)相應(yīng)分離效應(yīng)的周期進(jìn)行調(diào)整.時(shí)序中含有多種周期成分,應(yīng)針對各個(gè)周期進(jìn)行數(shù)次移動(dòng)平均,由式(2)計(jì)算得到各個(gè)頻率所對應(yīng)中心移動(dòng)平均產(chǎn)生的各撓度效應(yīng)分量值,依次從原始序列減去求得的各個(gè)中心移動(dòng)平均序列,最后得到趨勢項(xiàng).移動(dòng)平均流程如圖1所示.

    圖1 移動(dòng)平均流程圖

    其中,第t期的移動(dòng)平均值為

    (2)

    2撓度變形趨勢預(yù)測模型

    時(shí)間序列分析從動(dòng)態(tài)角度描述某一現(xiàn)象的內(nèi)在數(shù)量關(guān)系及變化規(guī)律,由于建模過程簡單、預(yù)報(bào)精度高,因而廣泛應(yīng)用于變形監(jiān)測的數(shù)據(jù)處理中.橋梁撓度變形是以時(shí)間為特征的數(shù)據(jù)集,影響因素多樣,形變機(jī)理復(fù)雜,難以用確定性的函數(shù)進(jìn)行表達(dá),但長期監(jiān)測數(shù)據(jù)在統(tǒng)計(jì)上具有一定的相關(guān)關(guān)系,因此可采用時(shí)間序列分析方法來表征橋梁結(jié)構(gòu)的特性.設(shè){Xt,t=1,2,…,n}是零均值的平穩(wěn)序列,其值與前p階的值和前q階的干擾項(xiàng)有關(guān),根據(jù)多元線性回歸的思想,可得ARMA(autoregressive and moving average)模型[13].

    (3)

    ARMA模型要求建模的基礎(chǔ)數(shù)據(jù)具有平穩(wěn)性,即序列的統(tǒng)計(jì)特征不隨時(shí)間的推移而發(fā)生改變.而一般的工程數(shù)據(jù),如橋梁撓度變形,由于受到自身結(jié)構(gòu)特點(diǎn)及外界環(huán)境因素影響,無法滿足條件,需要進(jìn)行平穩(wěn)化處理,本文采用差分法,設(shè)Xt經(jīng)過d次差分后得到平穩(wěn)的新序列Yt,則

    (4)

    式中:B為后移算子,BdXt=Xt-d,d為差分階次.將Yt代入式(3)中進(jìn)行建模,即得到ARIMA模型.ARIMA模型的建立主要包括以下4個(gè)步驟.

    (1) 平穩(wěn)性檢驗(yàn).平穩(wěn)性檢驗(yàn)可通過時(shí)序圖直接判讀,也可利用單位根檢驗(yàn)來判斷,常用的單位根檢驗(yàn)方法有DF(Dickey-Fuller)檢驗(yàn)法和ADF(augmented Dickey-Fuller)檢驗(yàn)法.

    (5)

    (6)

    (7)

    (8)

    式中:ε*為殘差向量;β為待估參數(shù);Q為殘差平方和;βk為第k次參數(shù)計(jì)算結(jié)果.

    A中的各元素可由下式計(jì)算得到:

    (9)

    實(shí)際計(jì)算時(shí),需要設(shè)置一個(gè)合適的迭代閾值δ,βk+1-βk<δ,例如可取δ=10-4.

    (3) 模型定階AIC.AIC(Akaike information criterion)準(zhǔn)則由日本統(tǒng)計(jì)學(xué)家Akaike于1974年提出,該準(zhǔn)則體現(xiàn)了模型擬合的優(yōu)劣程度.

    (10)

    (4) 模型χ2檢驗(yàn).對已確定的模型進(jìn)行白噪聲(模型殘差)檢驗(yàn)[15].構(gòu)造χ2檢驗(yàn)的統(tǒng)計(jì)量

    (11)

    3工程應(yīng)用

    3.1數(shù)據(jù)來源

    黃岡長江大橋?yàn)橹骺?67 m的鋼桁梁斜拉橋,上層為公路橋面,下層為雙線鐵路,2014年6月建成通車后結(jié)構(gòu)健康監(jiān)測系統(tǒng)也正式開始運(yùn)行,監(jiān)測包含環(huán)境參數(shù)、結(jié)構(gòu)動(dòng)力特性、主梁撓度、拉索索力等各項(xiàng)內(nèi)容.主梁撓度監(jiān)測采用基于差壓變送器的方法,儀器參考精度為0.065%(與所設(shè)置量程相關(guān)),傳感器安裝如圖2所示.測點(diǎn)布置(圖3)于邊跨、輔助跨跨中,主跨撓度較大,受荷載作用各點(diǎn)撓度變換差異明顯,因此在1/4跨、跨中、3/4跨布設(shè)測點(diǎn),3號塔處為基準(zhǔn)點(diǎn).

    圖2 傳感器安裝

    圖3 測點(diǎn)布置

    3.2撓度信號的模擬與分離

    鋼桁梁為超靜定結(jié)構(gòu),雖然熱傳導(dǎo)性能較好,但是考慮日照輻射、結(jié)構(gòu)方位、周圍地理環(huán)境等因素,結(jié)構(gòu)內(nèi)部各構(gòu)件間存在局部溫差,構(gòu)件間的變形相互約束,產(chǎn)生溫度次應(yīng)力,加之溫度分布不均在結(jié)構(gòu)各層之間產(chǎn)生的溫度自應(yīng)力,二者對結(jié)構(gòu)受力的影響十分顯著[16].鋼桁梁斜拉橋構(gòu)造復(fù)雜,準(zhǔn)確計(jì)算溫度效應(yīng)存在較大的困難,利用實(shí)測數(shù)據(jù)通過相關(guān)性分析可以得到溫度荷載下鋼桁梁的撓度行為.圖4為黃岡長江大橋典型撓度日變化趨勢圖,可以看出實(shí)測數(shù)據(jù)與上文撓度加法模型的描述較為符合,日變化趨勢表現(xiàn)為規(guī)則的正弦曲線,而列車荷載在跨中產(chǎn)生的瞬時(shí)下?lián)峡蛇_(dá)40 mm.

    圖4 典型撓度日變化時(shí)程

    對于溫度與撓度的相關(guān)性(如圖5所示),雖然二者存在一定的線性關(guān)系,但其區(qū)間很寬,回歸估計(jì)的精確性較低,因而通過建立溫度與撓度回歸模型以剔除溫度荷載的做法并不適用.但回歸分析的結(jié)果可以用來估計(jì)單位溫度影響下?lián)隙鹊拇笾赂淖兞浚鳛閿?shù)值模擬的依據(jù),結(jié)合溫度變化規(guī)律,進(jìn)而得到模擬撓度值.擬合線性模型為

    f=-2.68w+79.73

    (12)

    式中:f為撓度;w為溫度.由模型可知,當(dāng)升溫1℃時(shí),主跨跨中將產(chǎn)生2.68 mm(設(shè)為d)的下?lián)?

    圖5 溫度與撓度相關(guān)性散點(diǎn)圖

    考察黃岡地區(qū)氣候條件,結(jié)合年度溫度變化的規(guī)律,設(shè)1 d中整體日照溫度變化wd在0~16℃,年溫差wy為30~40℃,溫度均為正弦變化的周期函數(shù)(以1 s為1 h計(jì)),溫度與橋梁結(jié)構(gòu)變形之間關(guān)系按回歸分析結(jié)果來考慮.根據(jù)撓度加法模型,總撓度f0由日溫差效應(yīng)(f1)、年溫差效應(yīng)(f2)及長期下?lián)?f3)構(gòu)成,假設(shè)隨著時(shí)間的推移,長期撓度將趨于平穩(wěn),因此利用冪函數(shù)來擬合長期變形趨勢.綜上所述,可得模擬撓度各分量的計(jì)算公式為

    (13)

    式中:t為時(shí)間,采樣頻率為1 Hz,模擬的總撓度信號f0如圖6所示.

    圖6 總撓度模擬值

    采用Welch法對f0進(jìn)行頻譜分析,要準(zhǔn)確確定周期長度,要求頻率分辨率盡量要高,應(yīng)增加進(jìn)行快速傅里葉變換(FFT)的點(diǎn)數(shù),窗函數(shù)選擇Hanning窗,窗口長度與FFT計(jì)算點(diǎn)數(shù)相等,分段序列重疊長度取窗口長度一半,結(jié)果如圖7.按照1.2節(jié)中的策略對f0進(jìn)行分離,各撓度分量與模擬值比較如圖8~10,可見分離值與模擬值的曲線波動(dòng)規(guī)律基本一致.但考慮到不同隨機(jī)幅值的波動(dòng)影響到分離結(jié)果的優(yōu)劣,為檢驗(yàn)算法的可靠性,進(jìn)行100次試驗(yàn),統(tǒng)計(jì)試驗(yàn)中模擬值與分離值的相關(guān)系數(shù),結(jié)果如表1.日溫差效應(yīng)、年溫差效應(yīng)及趨勢項(xiàng)模擬值與分離值的相關(guān)系數(shù)在0.8以上的占了100%,表明分離值與模擬值存在強(qiáng)相關(guān),算法的分離結(jié)果驗(yàn)證了中心移動(dòng)平均算法用于撓度分離的可靠性與有效性.

    3.3實(shí)測數(shù)據(jù)的分離

    以主跨跨中測點(diǎn)為例,由于系統(tǒng)連續(xù)、實(shí)時(shí)運(yùn)行,采集的數(shù)據(jù)量較大,先利用小波分析理論剔除列車荷載和噪聲的影響,再以1 h平均值作為該時(shí)段的特征值.圖11為2014年7月7日至9月10日的監(jiān)測數(shù)據(jù),其總體呈現(xiàn)連續(xù)正弦曲線的波動(dòng)特征,短期內(nèi)受溫差影響幅值有所變化,但周期長度較為穩(wěn)定,這為溫度效應(yīng)的分離提供了基礎(chǔ).實(shí)測數(shù)據(jù)的時(shí)間跨度為數(shù)月,從頻域信息中獲得日溫差效應(yīng)撓度的周期仍然近似于1“天”,而長期溫度效應(yīng)的周期長度則接近于“月”,因此以“月”周期代替“年”周期進(jìn)行長期溫度效應(yīng)的分離,分離后的各撓度分量如圖12a至12c所示.

    圖7 模擬信號的頻譜分析

    圖8 日溫差效應(yīng)分離前后對比

    圖9 年溫差效應(yīng)分離前后對比

    圖10 長期趨勢分離前后對比

    表1 試驗(yàn)相關(guān)系數(shù)統(tǒng)計(jì)

    圖11 撓度長期觀測數(shù)據(jù)

    由圖12可見,分離的日溫差效應(yīng)較平穩(wěn),相鄰時(shí)刻的幅值也十分接近,月溫差效應(yīng)的頻率明顯大于日溫差效應(yīng),與實(shí)際相符,分離出的趨勢項(xiàng)總體表現(xiàn)為上拱,這由于以“月”代替“年”作為長期溫度效應(yīng)周期而導(dǎo)致短期溫差無法消除,主梁表現(xiàn)為受降溫影響而產(chǎn)生上拱.采用中心移動(dòng)平均法對實(shí)測溫度數(shù)據(jù)進(jìn)行分離,這里不考慮驟然溫差的影響,分離出日溫差變化和月溫差變化,結(jié)果如圖12d至12f,剔除日溫度變化、年溫度變化后的殘余溫度整體為降溫趨勢,符合溫度的季節(jié)特性,同時(shí)與回歸分析及撓度的分離結(jié)果相吻合,這印證了主梁長期趨勢變化的根源,從另一面也驗(yàn)證了分離算法的正確性.

    3.4長期趨勢的預(yù)測

    將提取的撓度長期趨勢進(jìn)行模型擬合,經(jīng)過試驗(yàn)計(jì)算,選取500個(gè)數(shù)據(jù)作為建模基礎(chǔ)數(shù)據(jù)擬合效果較好.由圖12c可知原始數(shù)據(jù)具有明顯的趨勢項(xiàng),經(jīng)過二階差分后,時(shí)序的ADF統(tǒng)計(jì)檢驗(yàn)量為-4.073,小于1%時(shí)顯著水平的臨界值,表明二階差分后時(shí)序?yàn)槠椒€(wěn)序列,滿足建模條件.設(shè)定模型初始階數(shù),由低階至高階擬合,由于低階模型就能以較高精度實(shí)現(xiàn)擬合,因此限定擬合最高階次.分別計(jì)算各模型AIC值,部分模型AIC值如表2.搜索模型中AIC最小值,當(dāng)p=4,q=5時(shí),AIC取到最小,則模型應(yīng)為ARIMA(4,2,5),進(jìn)行卡方檢驗(yàn),殘差服從正態(tài)分布,模型檢驗(yàn)通過.從圖13可知,模型的擬合效果良好.對后50期的數(shù)據(jù)進(jìn)行預(yù)測,由圖14可知其短期預(yù)測殘差未超過1 mm,最大相對誤差為8.61%,而前20期中相對誤差均小于5%,但隨預(yù)測期數(shù)的增加預(yù)測殘差也隨之增大,實(shí)際操作中可采用加窗的方法,以20期為間隔,更新模型的起算數(shù)據(jù),可有效提高預(yù)測精度.

    a 日溫差效應(yīng)撓度

    b 月溫差效應(yīng)撓度

    c 撓度長期趨勢

    d 日溫差變化

    e 月溫差變化

    f 殘余溫度

    表2 部分模型AIC值統(tǒng)計(jì)

    4結(jié)論

    針對現(xiàn)有撓度溫度效應(yīng)分離算法研究不足、實(shí)用性差的問題,提出基于中心移動(dòng)平均法的撓度-溫度效應(yīng)分離策略,將Welch法獲得的信號分量周期長度作為移動(dòng)平均半周期延遲的參數(shù),避免了移動(dòng)平均的滯后性,實(shí)現(xiàn)移動(dòng)步長的自動(dòng)確定.通過多次數(shù)值模擬試驗(yàn),得到模擬值與分離值存在強(qiáng)相關(guān)的結(jié)論,驗(yàn)證了算法的有效性.同時(shí)根據(jù)橋梁結(jié)構(gòu)長期變形趨勢不穩(wěn)定的特點(diǎn),以實(shí)測數(shù)據(jù)為研究對象,采用時(shí)間序列分析建立了撓度趨勢的ARIMA預(yù)測模型,模型擬合效果較好,預(yù)測精度高,可為橋梁的結(jié)構(gòu)安全提供預(yù)警.

    圖13 模型擬合效果

    圖14 模型預(yù)測殘差值

    參考文獻(xiàn):

    [1]陳宏玉,黃騰,李桂華. 基于灰關(guān)聯(lián)GM(1,N)模型的特大跨徑橋梁撓度變化分析[J]. 測繪通報(bào),2009,55(1):31.

    CHEN Hongyu,HUANG Teng,LI Guihua. Analysis of oversize span bridge deflection based on grey relation GM(1,N)[J]. Bulletin of Surveying and Mapping,2009,55(1):31.

    [2]宗周紅,鐘儒勉,鄭沛娟,等. 基于健康監(jiān)測的橋梁結(jié)構(gòu)損傷預(yù)后和安全預(yù)后研究進(jìn)展及挑戰(zhàn)[J]. 中國公路學(xué)報(bào),2014,27(12):46.

    ZONG Zhouhong,ZHONG Rumian,ZHENG Peijuan,etal. Damage and safety prognosis of bridge structures based on structural health monitoring[J]. China Journal of Highway and Transport,2014,27(12):46.

    [3]James M W, Brownjohn,Ki-Young Koo,etal. Operational deformations in long-span bridges[J]. Structure and Infrastructure Engineering, 2015, 11(4): 556.

    [4]Ghali A,Elbadry M,Megally S. Two-year deflections of the Confederation Bridge[J]. Canadian Journal of Civil Engineering, 2000, 27(6): 1139.

    [5]E J Cross,K Y Koo,J M W Brownjohn,etal. Long-term monitoring and data analysis of the Tamar Bridge[J]. Mechanical Systems and Signal Processing, 2013, 35(1): 16.

    [6]Tang C H, Yang H, Liu X P,etal. Separation of bridge deflection signals based on ICA[J].Advanced Materials Research, 2011, 374-377(10): 2090.

    [7]梁宗保.基于監(jiān)測信息統(tǒng)計(jì)分析的橋梁結(jié)構(gòu)安全評價(jià)研究[D].重慶:重慶大學(xué),2006.

    LIANG Zongbao. Research based on the monitoring information statistics and analysis of the bridge structure safety evaluation[D]. Chongqing: Chongqing University,2006.

    [8]劉夏平,楊紅,孫卓,等. 基于奇異值分解的橋梁撓度分離研究[J]. 中山大學(xué)學(xué)報(bào):自然科學(xué)版,2013,52(3):11.

    LIU Xiaping,YANG Hong,SUN Zhuo,etal. Separation study of bridge deflection based on singular value decomposition[J]. Acta Scientiarum Naturalium Universitatis Sunyatseni,2013,52(3):11.

    [9]劉綱,邵毅敏,黃宗明,等. 長期監(jiān)測中結(jié)構(gòu)溫度效應(yīng)分離的一種新方法[J].工程力學(xué),2010,27(3):55.

    LIU Gang,SHAO Yimin,HUANG Zongming,etal. A new method to separate temperature effect from long-term structural health monitoring data[J]. Engineering Mechanics,2010,27(3):55.

    [10]劉夏平,楊紅,孫卓,等. 基于LS-SVM的橋梁撓度監(jiān)測中溫度效應(yīng)分離[J]. 鐵道學(xué)報(bào),2012,34(10):91.

    LIU Xiaping,YANG Hong,SUN Zhuo,etal. Study on separation of bridge deflection temperature effect based on LS-SVM[J]. Journal of the China Railway Society,2012,34(10):91.

    [11]楊紅,孫卓,劉夏平,等. 基于多最小二乘支持向量機(jī)的橋梁溫度撓度效應(yīng)的分離[J]. 振動(dòng)與沖擊,2014,33(1):71.

    YANG Hong,SUN zhuo,LIU Xiaping,etal. Separation of bridge temperature deflection effect based on M-S-SVM[J]. Journal of Vibration and Shock,2014,33(1):71.

    [12]彭令,牛瑞卿,吳婷. 時(shí)間序列分析與支持向量機(jī)的滑坡位移預(yù)測[J]. 浙江大學(xué)學(xué)報(bào):工學(xué)版,2013,47(9):1672.

    PENG Ling,NIU Ruiqing,WU Ting. Time series analysis and support vector machine for landslide displacement prediction[J]. Journal of Zhejiang University:Engineering Science,2013,47(9):1672.

    [13]王燕. 應(yīng)用時(shí)間序列分析[M]. 北京:中國人民大學(xué)出版社,2008.

    WANG Yan. Applied time series analysis[M]. Beijing: China Renmin University Press,2008.

    [14]曾浩. 基于時(shí)間序列分析的結(jié)構(gòu)損傷診斷[D]. 武漢:武漢理工大學(xué),2006.

    ZENG Hao. Structural damage detection using time series analysis[D]. Wuhan:Wuhan University of Technology,2006.

    [15]劉寒冰,李國恒,譚國金,等. 基于時(shí)間序列的邊坡位移實(shí)時(shí)預(yù)測方法[J]. 吉林大學(xué)學(xué)報(bào):工學(xué)版,2012,42(S1):193.

    LIU Hanbing,LI Guoheng,TAN Guojin,etal. Real-time prediction method of slop displacement based on time series analysis[J]. Journal of Jilin University:Engineering and Technology Edition,2012,42(S1):193.

    [16]陳海鋒. 收縮徐變及溫度作用下懸索橋鋼-混組合橋面系受力性能分析[D].長沙:長沙理工大學(xué),2013.

    CHEN Haifeng. Performance analysis on the steel-concrete composite bridge deck of suspension bridge under the action of shrinkage and creep as well as temperature [D]. Changsha:Changsha University of Science & Technology,2013.

    收稿日期:2015-06-29

    基金項(xiàng)目:國家自然科學(xué)基金面上項(xiàng)目(41371423);國家“八六三”高技術(shù)研究發(fā)展計(jì)劃(2013AA12A201);江蘇省高校優(yōu)勢學(xué)科建設(shè)工程(SZBF2011-6-B35)

    中圖分類號:TU196+.1

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

    Study on Separation and Forecast of Long-term Deflection Based on Time Series Analysis

    CHEN Guoliang1,LIN Xungen1,YUE Qing2,LIU Hua2

    (1. School of Environment Science and Spatial Informatics, China University of Mining and Technology, Xuzhou 221116, China; 2. China Railway Major Bridge Reconnaissance & Design Group Co., Ltd., Wuhan 430050, China)

    Abstract:In view of the multi-scale property of deflection and the periodicity under temperature load, a strategy was proposed to separate the deflection and extract the long-term trend, which combines Welch method with Center Moving Averages Method. A prediction model of the deformation was established based on ARIMA model. Taking Bridge of Huanggang Yangtze River as an example to issue the correlation between temperature and deflection, and test of numerical simulation was made to validate the effectiveness of the algorithm. Processing and analysis of the measured data indicate that the approach of this article has high precision and reliability which could provide data support for damage identification of bridge.

    Key words:deflection; separation; long-term trend; ARIMA model; forecast

    第一作者: 陳國良(1977—),男,副教授,工學(xué)博士,主要研究方向?yàn)樽冃伪O(jiān)測、導(dǎo)航位置服務(wù).E-mail:chglcumt@163.com

    猜你喜歡
    ARIMA模型分離撓度
    Spontaneous multivessel coronary artery spasm diagnosed with intravascular ultrasound imaging:A case report
    基于時(shí)間序列模型的中國出口總額分析及預(yù)測
    基于R軟件的金融時(shí)間序列的預(yù)測分析
    基于Eviews上證綜合指數(shù)預(yù)測
    轉(zhuǎn)型背景下的民辦高校管理理念創(chuàng)新
    融合與分離:作為一種再現(xiàn)的巫術(shù)、文字與影像世界的構(gòu)成
    高效液相色譜技術(shù)在石油化工中的應(yīng)用分析
    基于ARIMA模型的滬銅期貨價(jià)格預(yù)測研究
    商(2016年27期)2016-10-17 06:06:02
    氣體分離提純應(yīng)用變壓吸附技術(shù)的分析
    科技視界(2016年20期)2016-09-29 11:02:20
    懸高測量在橋梁撓度快速檢測中的應(yīng)用
    高潮久久久久久久久久久不卡| 久久午夜亚洲精品久久| 香蕉丝袜av| 美女高潮到喷水免费观看| 精品国产一区二区三区久久久樱花| 韩国av一区二区三区四区| 欧美一级毛片孕妇| 精品久久蜜臀av无| 久久国产亚洲av麻豆专区| 日韩中文字幕欧美一区二区| 日韩欧美三级三区| 精品福利永久在线观看| 黄色丝袜av网址大全| 色综合婷婷激情| 飞空精品影院首页| videosex国产| 色婷婷av一区二区三区视频| 久久精品亚洲精品国产色婷小说| 亚洲国产精品sss在线观看 | 精品人妻在线不人妻| 性色av乱码一区二区三区2| 亚洲欧美激情综合另类| 在线观看舔阴道视频| 日韩免费av在线播放| 99精品欧美一区二区三区四区| 久久久久国内视频| 久久人妻av系列| 亚洲精品一卡2卡三卡4卡5卡| 搡老岳熟女国产| 亚洲人成77777在线视频| 欧美日韩福利视频一区二区| 在线观看免费高清a一片| 超碰97精品在线观看| 日韩欧美一区二区三区在线观看 | 咕卡用的链子| av中文乱码字幕在线| 天堂动漫精品| 叶爱在线成人免费视频播放| 免费av中文字幕在线| 高清毛片免费观看视频网站 | 成人特级黄色片久久久久久久| 色老头精品视频在线观看| 在线十欧美十亚洲十日本专区| 亚洲欧美精品综合一区二区三区| 欧美国产精品va在线观看不卡| 老熟妇乱子伦视频在线观看| 免费日韩欧美在线观看| 久久精品亚洲精品国产色婷小说| tube8黄色片| 国产精品自产拍在线观看55亚洲 | 亚洲第一av免费看| 久久亚洲真实| 国产激情久久老熟女| 女性被躁到高潮视频| 免费少妇av软件| 51午夜福利影视在线观看| 国产区一区二久久| 一本综合久久免费| 日韩 欧美 亚洲 中文字幕| 免费女性裸体啪啪无遮挡网站| 亚洲免费av在线视频| 久久精品国产综合久久久| 自线自在国产av| 成人免费观看视频高清| 欧美丝袜亚洲另类 | 老司机福利观看| 精品久久久精品久久久| 午夜日韩欧美国产| 国产高清videossex| 女性被躁到高潮视频| 好男人电影高清在线观看| 国产一区二区三区综合在线观看| 亚洲一码二码三码区别大吗| 国产精品美女特级片免费视频播放器 | 日韩视频一区二区在线观看| 亚洲av电影在线进入| 亚洲av第一区精品v没综合| 久久午夜亚洲精品久久| 99久久综合精品五月天人人| 国产在线精品亚洲第一网站| 咕卡用的链子| 成人特级黄色片久久久久久久| 亚洲美女黄片视频| 丰满的人妻完整版| 国产一区二区激情短视频| 久久久久久久久久久久大奶| 99精品久久久久人妻精品| 色尼玛亚洲综合影院| 丝袜美足系列| 少妇 在线观看| 国产精品影院久久| 国产又色又爽无遮挡免费看| 欧美av亚洲av综合av国产av| 在线观看免费日韩欧美大片| 精品福利永久在线观看| 免费观看a级毛片全部| 欧美另类亚洲清纯唯美| 亚洲色图综合在线观看| 老司机亚洲免费影院| 悠悠久久av| 在线观看免费视频网站a站| 精品久久久久久久毛片微露脸| 18在线观看网站| 亚洲成国产人片在线观看| 高潮久久久久久久久久久不卡| 国产精品国产高清国产av | 亚洲精品美女久久久久99蜜臀| 国产精品综合久久久久久久免费 | 一二三四在线观看免费中文在| 国产又爽黄色视频| 久久精品91无色码中文字幕| 亚洲色图综合在线观看| 丁香欧美五月| 桃红色精品国产亚洲av| 人人妻人人添人人爽欧美一区卜| 午夜影院日韩av| а√天堂www在线а√下载 | 久久精品人人爽人人爽视色| 成人三级做爰电影| 精品国产亚洲在线| 日韩熟女老妇一区二区性免费视频| 亚洲人成伊人成综合网2020| 欧美性长视频在线观看| 欧美黄色片欧美黄色片| 免费不卡黄色视频| 在线观看免费视频日本深夜| 亚洲性夜色夜夜综合| 中出人妻视频一区二区| av免费在线观看网站| 久久精品国产清高在天天线| 老司机午夜十八禁免费视频| 成年人午夜在线观看视频| av视频免费观看在线观看| 欧美日韩成人在线一区二区| 男男h啪啪无遮挡| 精品亚洲成国产av| 桃红色精品国产亚洲av| 国产成+人综合+亚洲专区| 国产不卡一卡二| 免费看十八禁软件| 国产99久久九九免费精品| 黑人巨大精品欧美一区二区蜜桃| 精品免费久久久久久久清纯 | 岛国毛片在线播放| 成年女人毛片免费观看观看9 | 亚洲人成77777在线视频| 99久久综合精品五月天人人| 丝袜在线中文字幕| 99国产精品一区二区蜜桃av | 亚洲人成77777在线视频| 777米奇影视久久| 久久香蕉精品热| 国产精品九九99| 女人被狂操c到高潮| 下体分泌物呈黄色| 国产三级黄色录像| 欧美另类亚洲清纯唯美| 亚洲 国产 在线| 久久中文字幕人妻熟女| 欧美大码av| 一进一出抽搐gif免费好疼 | 一级片'在线观看视频| 又紧又爽又黄一区二区| 欧美日韩中文字幕国产精品一区二区三区 | 新久久久久国产一级毛片| 高清黄色对白视频在线免费看| 日韩欧美在线二视频 | 性少妇av在线| 精品国产亚洲在线| 亚洲国产精品一区二区三区在线| 人人妻人人添人人爽欧美一区卜| 日韩 欧美 亚洲 中文字幕| 91国产中文字幕| 国产成人精品久久二区二区免费| 日韩欧美一区二区三区在线观看 | 国产精品亚洲一级av第二区| 国产男女超爽视频在线观看| 久久午夜综合久久蜜桃| 欧美av亚洲av综合av国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 水蜜桃什么品种好| 露出奶头的视频| 美女扒开内裤让男人捅视频| 超碰97精品在线观看| 欧美日韩亚洲综合一区二区三区_| 一级片'在线观看视频| 身体一侧抽搐| 狠狠婷婷综合久久久久久88av| 伦理电影免费视频| 午夜影院日韩av| 丝袜人妻中文字幕| 亚洲午夜理论影院| 深夜精品福利| 国产激情久久老熟女| 丝袜在线中文字幕| 亚洲精品自拍成人| 老司机深夜福利视频在线观看| 18禁国产床啪视频网站| 十分钟在线观看高清视频www| 在线观看66精品国产| 成人亚洲精品一区在线观看| 国产精品一区二区免费欧美| 欧美av亚洲av综合av国产av| 国产精品美女特级片免费视频播放器 | 精品久久久久久电影网| 日韩中文字幕欧美一区二区| 在线看a的网站| 淫妇啪啪啪对白视频| 日韩成人在线观看一区二区三区| 久99久视频精品免费| 女性被躁到高潮视频| 天堂√8在线中文| 久久久国产成人精品二区 | 亚洲第一青青草原| 大码成人一级视频| 午夜精品国产一区二区电影| 一边摸一边抽搐一进一小说 | 亚洲精品粉嫩美女一区| 欧美精品一区二区免费开放| 国产成人啪精品午夜网站| 欧美最黄视频在线播放免费 | 国内久久婷婷六月综合欲色啪| 下体分泌物呈黄色| 人妻一区二区av| 午夜福利在线观看吧| 国产精品 欧美亚洲| 9热在线视频观看99| 黄色女人牲交| 亚洲精品自拍成人| 久久ye,这里只有精品| 精品少妇一区二区三区视频日本电影| a级片在线免费高清观看视频| 69精品国产乱码久久久| 免费观看精品视频网站| 两个人免费观看高清视频| 热99re8久久精品国产| 少妇 在线观看| 国产精品 欧美亚洲| 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 窝窝影院91人妻| 日韩精品免费视频一区二区三区| 国产欧美日韩综合在线一区二区| 手机成人av网站| 女性生殖器流出的白浆| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 国产麻豆69| 国产亚洲精品久久久久久毛片 | 成年版毛片免费区| 一个人免费在线观看的高清视频| 成人18禁高潮啪啪吃奶动态图| 午夜影院日韩av| 如日韩欧美国产精品一区二区三区| 女性被躁到高潮视频| 女性生殖器流出的白浆| 免费在线观看日本一区| 午夜福利欧美成人| 国产97色在线日韩免费| 亚洲av成人不卡在线观看播放网| 香蕉国产在线看| 在线十欧美十亚洲十日本专区| 欧美激情极品国产一区二区三区| 久久久久久久久免费视频了| 国产蜜桃级精品一区二区三区 | 久久精品国产亚洲av香蕉五月 | av欧美777| 午夜老司机福利片| ponron亚洲| 麻豆av在线久日| 久久精品亚洲精品国产色婷小说| 国产成人精品在线电影| 国产1区2区3区精品| 精品国产美女av久久久久小说| 夜夜夜夜夜久久久久| 老司机靠b影院| 母亲3免费完整高清在线观看| 91九色精品人成在线观看| 中文字幕人妻丝袜一区二区| 久久人妻av系列| 国产淫语在线视频| 国产成人欧美在线观看 | 国产单亲对白刺激| 国产在线一区二区三区精| 国产av又大| 一级a爱视频在线免费观看| 久久中文字幕人妻熟女| 精品国产国语对白av| 一二三四社区在线视频社区8| √禁漫天堂资源中文www| 精品亚洲成国产av| 两人在一起打扑克的视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月天丁香| 成人永久免费在线观看视频| 岛国毛片在线播放| 窝窝影院91人妻| 欧美激情 高清一区二区三区| 校园春色视频在线观看| 欧美av亚洲av综合av国产av| 欧美丝袜亚洲另类 | 欧美性长视频在线观看| 午夜福利在线免费观看网站| 黄片播放在线免费| 欧美老熟妇乱子伦牲交| 日韩欧美在线二视频 | 国产精品一区二区在线不卡| 一区二区日韩欧美中文字幕| 亚洲avbb在线观看| 久久国产精品男人的天堂亚洲| 免费在线观看亚洲国产| 好男人电影高清在线观看| 香蕉国产在线看| 国产亚洲av高清不卡| 交换朋友夫妻互换小说| 免费在线观看视频国产中文字幕亚洲| 欧美日韩亚洲综合一区二区三区_| 国产成人免费无遮挡视频| 亚洲男人天堂网一区| 国产精品久久电影中文字幕 | 日本一区二区免费在线视频| 国产欧美日韩一区二区三区在线| av福利片在线| 三上悠亚av全集在线观看| 黄片大片在线免费观看| 侵犯人妻中文字幕一二三四区| а√天堂www在线а√下载 | 如日韩欧美国产精品一区二区三区| 变态另类成人亚洲欧美熟女 | 大香蕉久久网| 亚洲欧美日韩高清在线视频| 中文字幕人妻丝袜一区二区| 久久午夜综合久久蜜桃| 在线免费观看的www视频| 亚洲综合色网址| 日韩成人在线观看一区二区三区| 久久精品国产99精品国产亚洲性色 | 很黄的视频免费| 三上悠亚av全集在线观看| 男人操女人黄网站| 美女视频免费永久观看网站| 中文字幕人妻熟女乱码| 欧美日韩亚洲高清精品| 女警被强在线播放| 亚洲 国产 在线| 欧美日韩福利视频一区二区| 丝袜在线中文字幕| 欧美乱妇无乱码| 欧美乱码精品一区二区三区| 亚洲午夜精品一区,二区,三区| 色播在线永久视频| 欧美日韩黄片免| 亚洲精品一二三| 中文字幕人妻熟女乱码| 精品卡一卡二卡四卡免费| 老汉色av国产亚洲站长工具| 我的亚洲天堂| 亚洲精品美女久久av网站| 人妻丰满熟妇av一区二区三区 | 欧美精品一区二区免费开放| 悠悠久久av| 波多野结衣一区麻豆| 成人特级黄色片久久久久久久| 国产成人影院久久av| 国产精品一区二区精品视频观看| 久久久久国内视频| 亚洲精品在线美女| 久9热在线精品视频| 999久久久精品免费观看国产| 高清毛片免费观看视频网站 | 精品国产国语对白av| 久久久久久人人人人人| 国产精品乱码一区二三区的特点 | 男人的好看免费观看在线视频 | 久久久精品国产亚洲av高清涩受| 欧美精品亚洲一区二区| 女性生殖器流出的白浆| 国产精品免费大片| ponron亚洲| 亚洲精华国产精华精| 91九色精品人成在线观看| 另类亚洲欧美激情| 熟女少妇亚洲综合色aaa.| 啦啦啦免费观看视频1| 成人三级做爰电影| 热99久久久久精品小说推荐| 亚洲在线自拍视频| 久久精品国产清高在天天线| 丰满迷人的少妇在线观看| 国产黄色免费在线视频| 久99久视频精品免费| 脱女人内裤的视频| 乱人伦中国视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美激情高清一区二区三区| xxx96com| 国产又爽黄色视频| 国产精品美女特级片免费视频播放器 | 亚洲五月色婷婷综合| 日韩熟女老妇一区二区性免费视频| 精品少妇久久久久久888优播| 啦啦啦 在线观看视频| 欧美最黄视频在线播放免费 | 18禁美女被吸乳视频| 欧美日韩瑟瑟在线播放| 久久精品91无色码中文字幕| 国产无遮挡羞羞视频在线观看| 欧美乱妇无乱码| 午夜视频精品福利| 两人在一起打扑克的视频| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美在线二视频 | 精品高清国产在线一区| 欧美日韩乱码在线| 国产精品1区2区在线观看. | 欧美精品亚洲一区二区| 一区二区三区精品91| 久久婷婷成人综合色麻豆| 操出白浆在线播放| 国产日韩欧美亚洲二区| 欧洲精品卡2卡3卡4卡5卡区| 久久精品aⅴ一区二区三区四区| av一本久久久久| 国产91精品成人一区二区三区| 精品一区二区三区av网在线观看| 国产一区二区激情短视频| 9191精品国产免费久久| 在线看a的网站| 国产精品免费视频内射| 波多野结衣一区麻豆| 嫩草影视91久久| 亚洲成国产人片在线观看| 久久狼人影院| 成人18禁在线播放| 天堂√8在线中文| 精品卡一卡二卡四卡免费| 国产伦人伦偷精品视频| 亚洲欧美日韩另类电影网站| 麻豆国产av国片精品| 嫩草影视91久久| 久久精品国产99精品国产亚洲性色 | 国产精品欧美亚洲77777| 亚洲精品一二三| 久久青草综合色| 国产一区二区三区在线臀色熟女 | 国产xxxxx性猛交| 国产精品永久免费网站| 99在线人妻在线中文字幕 | 欧美精品人与动牲交sv欧美| 国产99白浆流出| 热99re8久久精品国产| 亚洲国产精品sss在线观看 | 少妇猛男粗大的猛烈进出视频| 国产一区二区三区在线臀色熟女 | 亚洲成人免费电影在线观看| a在线观看视频网站| 欧美亚洲日本最大视频资源| 后天国语完整版免费观看| 亚洲专区国产一区二区| 国产成人影院久久av| 午夜视频精品福利| 狂野欧美激情性xxxx| 国产精品久久久人人做人人爽| 黄色怎么调成土黄色| 午夜福利在线免费观看网站| 老熟妇乱子伦视频在线观看| tocl精华| 国产在线一区二区三区精| 成熟少妇高潮喷水视频| 在线免费观看的www视频| 真人做人爱边吃奶动态| 国产成人精品久久二区二区91| 一级黄色大片毛片| 操出白浆在线播放| 美女午夜性视频免费| 搡老乐熟女国产| 国产精品.久久久| 一本大道久久a久久精品| 母亲3免费完整高清在线观看| 妹子高潮喷水视频| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 夜夜躁狠狠躁天天躁| 99精国产麻豆久久婷婷| 又紧又爽又黄一区二区| 热99re8久久精品国产| 欧美在线黄色| 我的亚洲天堂| 少妇裸体淫交视频免费看高清 | 国产成人欧美| 国产区一区二久久| 91麻豆精品激情在线观看国产 | 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| a级片在线免费高清观看视频| 精品高清国产在线一区| 夜夜爽天天搞| 999久久久国产精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产91精品成人一区二区三区| 男男h啪啪无遮挡| 精品福利永久在线观看| 久久久久久人人人人人| 大片电影免费在线观看免费| 日本黄色日本黄色录像| 老司机亚洲免费影院| 青草久久国产| 波多野结衣av一区二区av| av一本久久久久| 女性生殖器流出的白浆| 午夜激情av网站| 一区在线观看完整版| 久久国产乱子伦精品免费另类| 久久午夜综合久久蜜桃| 国产av精品麻豆| 九色亚洲精品在线播放| 日韩有码中文字幕| 精品一品国产午夜福利视频| 脱女人内裤的视频| 男女下面插进去视频免费观看| 国产成人欧美| 一级片'在线观看视频| 99国产精品免费福利视频| 中文字幕精品免费在线观看视频| 国产精品成人在线| 久久国产乱子伦精品免费另类| 最新美女视频免费是黄的| 黄色视频,在线免费观看| 亚洲五月色婷婷综合| 很黄的视频免费| 精品亚洲成a人片在线观看| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 亚洲成人免费av在线播放| 国产国语露脸激情在线看| 成人av一区二区三区在线看| 国产色视频综合| 国产精品亚洲一级av第二区| 精品久久久久久电影网| 久久天堂一区二区三区四区| 在线观看舔阴道视频| 18禁裸乳无遮挡免费网站照片 | 黑丝袜美女国产一区| 黑人巨大精品欧美一区二区mp4| 人人妻人人澡人人看| 一区二区日韩欧美中文字幕| 亚洲中文日韩欧美视频| 久久精品成人免费网站| videos熟女内射| 中文字幕最新亚洲高清| 精品国产乱码久久久久久男人| 久久久精品区二区三区| 精品久久蜜臀av无| 国产午夜精品久久久久久| a级毛片黄视频| bbb黄色大片| 午夜福利乱码中文字幕| 午夜日韩欧美国产| 嫁个100分男人电影在线观看| 亚洲成人免费av在线播放| 淫妇啪啪啪对白视频| 亚洲精品中文字幕一二三四区| 黄色 视频免费看| 色综合欧美亚洲国产小说| 久久精品aⅴ一区二区三区四区| 啦啦啦视频在线资源免费观看| 母亲3免费完整高清在线观看| 日韩制服丝袜自拍偷拍| 中文字幕最新亚洲高清| 狠狠婷婷综合久久久久久88av| 婷婷成人精品国产| 99久久综合精品五月天人人| 亚洲精品国产精品久久久不卡| 国产精品一区二区精品视频观看| 亚洲欧美激情综合另类| 一本一本久久a久久精品综合妖精| 99国产精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 一区福利在线观看| 精品国产国语对白av| 两性夫妻黄色片| 国产一区在线观看成人免费| 少妇猛男粗大的猛烈进出视频| 久久久久国内视频| 日韩成人在线观看一区二区三区| 国产99久久九九免费精品| 电影成人av| 涩涩av久久男人的天堂| xxxhd国产人妻xxx| a级毛片在线看网站| videos熟女内射| 亚洲精品久久午夜乱码| 日韩人妻精品一区2区三区| 国产日韩欧美亚洲二区| xxxhd国产人妻xxx| 中亚洲国语对白在线视频| 人人妻,人人澡人人爽秒播| 国产日韩一区二区三区精品不卡| 亚洲欧美一区二区三区黑人| 久久人人97超碰香蕉20202| 欧美另类亚洲清纯唯美| 免费高清在线观看日韩| 妹子高潮喷水视频| 亚洲国产中文字幕在线视频| 99精国产麻豆久久婷婷| 老司机靠b影院| 黄色丝袜av网址大全| 男女之事视频高清在线观看| 18禁观看日本| 国产高清国产精品国产三级| 亚洲七黄色美女视频| 亚洲 国产 在线| 国产在线精品亚洲第一网站| 亚洲黑人精品在线| 精品亚洲成a人片在线观看|