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

    基于溫度撓度識別斜拉橋拉索損傷的研究

    2021-12-20 08:48:42譚冬梅陳方望
    振動與沖擊 2021年23期
    關鍵詞:斜拉橋拉索撓度

    譚冬梅,姚 歡,陳方望,吳 浩

    (1.武漢理工大學 道路橋梁與結(jié)構工程湖北省重點實驗室,武漢 430070;2.華中師范大學 城市與環(huán)境科學學院,武漢 430079)

    斜拉索作為斜拉橋的主要受力構件,了解其安全狀態(tài)尤為重要,故近年來斜拉索的損傷識別成為了研究的熱點,目前一大趨勢為采用實際監(jiān)測結(jié)構響應(如實測撓度)對拉索進行損傷識別。但橋梁結(jié)構在運營過程中承受多種荷載作用,因此得到的實時監(jiān)測撓度是多項因素共同作用的結(jié)果。據(jù)已有研究表明,這些影響因素主要包括車載、溫度作用及混凝土的收縮徐變等[1],若直接對多種荷載作用下的響應進行分析無法判斷橋梁所處的真實狀態(tài)。損傷信號的提取會受到較大干擾,故從分離出單項因素作用出發(fā)可很好避免這一點。其中溫度荷載在時間尺度上與動荷載存在很大差異,為從總體響應中分離出溫度響應提供了便利,同時分離出來的溫度響應也作為單項因素作用下的響應為判斷橋梁的運營狀態(tài)提供了可能。近年來,國內(nèi)外學者針對斜拉橋的溫度效應的提取和損傷診斷進行了一系列研究。

    丁幼亮等[2]研究了黃岡長江大橋在溫度荷載和車載分別作用下的變化規(guī)律,結(jié)果表明溫度撓度與主梁溫度之間存在較強的線性相關性,輔助跨跨中、主跨 1/4 跨、主跨 3/4跨與主跨跨中的溫度撓度具有較強的空間線性相關性;Cao等[3]借助湛江跨海灣大橋的健康監(jiān)測系統(tǒng)對該斜拉橋的主梁、拉索及橋塔的熱效應時延和梯度進行了分析,結(jié)果表明混凝土的溫度明顯滯后于空氣溫度;Yang等[4]研究了實測溫度與橋塔位移和主梁位移之間的關系,發(fā)現(xiàn)溫度與橋塔位移和主梁位移具有很強的線性關系,并構建了溫度與位移之間的回歸模型,擬合結(jié)果較好;劉夏平等[5-6]以溫度為輸入,溫度撓度為輸出分別構建了溫度與溫度撓度的LS-SVM(least squares support vector machines)及M-LS-SVM(multiple least squares support vector machines)模型;陳夏春等[7-8]構建了溫度與溫度撓度的多元線性回歸模型。其構建的溫度與溫度撓度間的模型,可通過損傷前后同樣溫度的輸入帶來輸出溫度撓度的變化來識別結(jié)構的損傷。吳向男等[9]總結(jié)了已有橋梁損傷診斷方法的優(yōu)缺點,并提出考慮環(huán)境影響和不確定因素是未來損傷診斷的發(fā)展方向,損傷診斷方法應更多的結(jié)合橋梁結(jié)構運營過程中實測監(jiān)測數(shù)據(jù);Wang等[10]在枝江斜拉橋上布置溫度傳感器和位移傳感器獲取了多個測點的溫度和撓度變化曲線,分析了橋梁的溫度分布規(guī)律以及結(jié)合溫度和撓度的變化規(guī)律,得出每個測點的撓度隨溫度變化的變化率(LSD)保持恒定的觀點,并基于此提出利用固定測點LSD的變化來識別橋梁的拉索損傷;Huang等[11]提出了一種基于卡爾曼濾波與協(xié)整的損傷識別新方法,該方法能夠有效剔除環(huán)境噪聲和溫度變化的影響并對結(jié)構進行損傷識別,并成功識別出了天津永和橋兩種損壞方案。

    針對拉索的損傷識別,考慮拉索損傷僅會影響局部測點在溫度作用下變化規(guī)律的特點,且相鄰測點間具有較強相關性的特點,本文提出一種只需要撓度數(shù)據(jù)即可對拉索進行損傷識別的方法。本方法解決了在運營中斜拉橋溫度分布規(guī)律不易測量的難題;并且由于實際溫度撓度變化緩慢的特點,便可降低采樣頻率,減少了數(shù)據(jù)的需求量,在一定程度上提升了計算效率,降低了操作的成本。

    1 實時在線斜拉橋的監(jiān)測撓度分析

    1.1 相鄰測點間溫度撓度的相關性分析

    基于武漢市某過江大跨斜拉橋在一年內(nèi)的長期健康監(jiān)測數(shù)據(jù),對沿橋梁縱向各測點的撓度進行相關性分析。如圖1所示為武漢某斜拉橋,為雙塔雙索面混合型斜拉橋,該橋在運營中結(jié)合北斗定位導航技術,對其主體結(jié)構的三維變形進行實時監(jiān)測。北斗監(jiān)測系統(tǒng)較傳統(tǒng)監(jiān)測方式,具有精度較高、不受氣候條件影響且實時高效率等特點,橋面測點布置如圖2所示。

    圖1 武漢某斜拉橋Fig.1 A cable-stayed bridge in Wuhan

    監(jiān)測橋梁的動態(tài)測點總共36個,采樣頻率為每秒一次。如圖3所示為選取的跨中測點(BD11和BD34)的實測撓度數(shù)據(jù),采樣時間為2018年06月17日—2018年06月18日。從圖3中可以看出,監(jiān)測橋梁主要受到車輛荷載和溫度作用的影響,為了能準確對溫度撓度進行分析,依據(jù)車輛荷載作用下響應的高頻特征及溫度作用下響應的低頻特征[12],本文采用集成經(jīng)驗模態(tài)分解[13-15]將溫度撓度從測點的總體響應中提取出來,分離得到BD11和BD34的溫度撓度,如圖4所示。

    (a)BD11

    (a)BD11

    對分解得到的相鄰測點間的溫度撓度做相關性分析,相關系數(shù)定義為

    (1)

    式中:cov(·)為序列間的協(xié)方差;si,sj為對應測點的溫度撓度;ρ的值越趨近于1,表明相關性越高。任意選取了2018年6月份3個時間段各測點的溫度撓度進行了相關系數(shù)計算,結(jié)果如表1所示。由表1可知:各測點溫度撓度均具有較高相似性,且隨著測點間距的增加,相關系數(shù)相對減小,但仍具有較高相關性。根據(jù)以上分析可近似認為相鄰測點間的溫度撓度存在較為理想的線性關系。

    表1 相鄰測點間溫度撓度的相關系數(shù)Tab.1 Correlation coefficient of temperature deflection between adjacent measuring points

    1.2 相鄰測點間線性回歸模型的建立

    根據(jù)對實際斜拉橋相鄰測點間溫度效應的相似性分析可知:斜拉橋相鄰測點間溫度撓度近似為一次線性關系,因此可采用線性回歸模型對之進行分析。一元線性回歸模型定義為

    Di(t)=ai+1,i×Di+1(t)+bi+1,i

    (2)

    Q(ai+1,i,bi)=Di(t)-ai+1,i×Di+1(t)-bi+1,i

    (3)

    式中:i為測點編號;ai+1,i和bi+1,i為待定系數(shù);Di為相應測點的溫度撓度;Di+1為相應測點相鄰測點的溫度撓度;Q(ai+1,i,bi)為預測殘差;對于模型參數(shù)ai+1,i和bi+1,i的求取,利用最小二乘法進行估值,即給定任意(Di(t),Di+1(t))的值,使式Q(ai+1,i,bi)的值最小即可。

    上述回歸模型的準確度通過回歸系數(shù)R2來評價,R2的取值越接近于1則回歸擬合程度越高,定義為

    (4)

    選取2018年06月01日—06月02日測點BD33、BD34、BD35、BD36及BD37的溫度撓度作為分析對象,建立RM(BD33,BD34),RM(BD34,BD35),RM(BD35,BD36)以及RM(BD36,BD37)4個回歸模型,其中RM(BD33,BD34)表示為D34(t)=a33,34×D33(t)+b33,34其他以此類推。各個回歸模型的擬合結(jié)果如圖5所示,各個回歸模型的回歸系數(shù),如表2所示。從圖5中可以看出,回歸擬合結(jié)果與實際值較接近,同時表2顯示回歸系數(shù)均在0.98以上,說明擬合程度較好。

    (a)RM(BD33,BD34)

    表2 部分測點回歸模型的回歸系數(shù)Tab.2 Regression coefficient of regression RM for several measuring points

    2 拉索損傷識別的數(shù)值仿真模擬

    2.1 溫度作用下主梁撓度對拉索損傷的敏感性分析

    探究溫度作用下主梁撓度對拉索損傷的敏感性。在Midas Civil軟件上建立武漢某過江斜拉橋有限元模型并對結(jié)構進行仿真分析,模型如圖6所示。全橋模型共有1 141個節(jié)點,1 126個單元。如圖7所示為拉索編號的設置及測點的布置示意圖。對于溫度的模擬采用整體升溫的形式,對模型整體升溫1 ℃,橋梁未受損傷情況下的變形圖,如圖8所示。

    圖6 武漢某斜拉橋模型Fig.6 A cable-stayed bridge model in Wuhan

    圖7 拉索編號的設置及測點的布置示意圖Fig.7 Schematic diagram of cable number setting and measuring point arrangement

    圖8 斜拉橋在溫度作用下的變形圖Fig.8 Deformation diagram of cable-stayed bridge under temperature

    斜拉橋拉索損傷的模擬可通過拉索截面面積的縮減或彈性模量的降低,本文采用對拉索截面面積的縮減來模擬拉索出現(xiàn)的損傷[16],即損傷率表示為

    (5)

    表3 拉索的損傷工況Tab.3 Damage condition of the cable

    (6)

    式中:DVTi,dmage為損傷后單位溫度產(chǎn)生的位移;DVTi,original為損傷前單位溫度產(chǎn)生的位移。

    表3中#SN01D10為#SN01拉索截面面積縮減10%,以此類推。

    各測點在不同種損傷工況下的φ(i)如圖9所示。從圖9(a)中可以看出:對于僅#SN01拉索出現(xiàn)損傷的情況,測點編號為D11、D12和D13的φ隨著#SN01拉索損傷程度的增加而增大,距離損傷拉索較遠測點的φ變化不明顯,且基本等于零。從圖9(b)中可以看出:隨著#SN17拉索損傷程度增加,測點編號為D15、D16和D17的其附近測點的φ也隨損傷程度而變大,且#SN01也因為出現(xiàn)損傷而使拉索附近測點D11、D12和D13的φ相比未損傷較大。由上述分析可知,出現(xiàn)損傷的拉索附近的主梁在單位溫度下的撓度會隨著損傷程度的增加而增大,未損傷拉索附近的主梁在單位溫度作用下基本無變化。

    (a)

    (a)單根拉索出現(xiàn)損傷

    2.2 溫度作用下拉索的損傷識別

    由1.2節(jié)可知相鄰測點間的溫度撓度可近似通過線性回歸表示,即

    Di(t)=ai+1,i×Di+1(t)+bi+1,i

    (7)

    依據(jù)文獻[18-19]對溫度效應與溫度之間的擬合,可做如下定義

    Di(t)=DVTi×T×sin(2×π×1/24)

    (8)

    Di+1(t)=DVTi+1×T×sin(2×π×1/24)

    (9)

    式中,T為溫度,當拉索未損傷時,由式(7)~式(9)得

    (10)

    當拉索損傷出現(xiàn)時,DVTi和DVTi+1都發(fā)生改變,但僅在拉索損傷處,其附近主梁撓度的DVTi變化較明顯,故定義拉索損傷后

    (11)

    則損傷前后ai+1,i的變化率為

    (12)

    式中:ai+1,i,damage為損傷后得到的系數(shù);ai+1,i,original為損傷前得到的系數(shù)。

    因此可以通過計算斜拉橋拉索在運營中ai+1,i是否變化來判定拉索損傷的位置,γ(i)確定拉索損傷的程度。ai+1,i的值發(fā)生變化,則測點i附近的拉索出現(xiàn)損傷,γ(i)值越大,則測點i附近的拉索出現(xiàn)的損傷越嚴重。

    為更好模擬實際情況,結(jié)合武漢地區(qū)天氣,擬定一年內(nèi)的日溫差為5~20 ℃,采用Matlab軟件生成365個5~20隨機數(shù)作為每天的日溫差,對模型施加時變溫度荷載,如圖10所示。其中為了顯示方便只顯示前10天的時變溫度荷載圖。

    圖10 模擬溫度Fig.10 Simulated temperature

    得到溫度作用下測點的溫度撓度后,以相鄰測點的溫度撓度分別作為回歸模型輸入和輸出,不同損傷工況下相鄰測點間回歸方程的回歸模型參數(shù)a的變化率γ,如圖11所示。從圖11可以看出:對于單根拉索出現(xiàn)損傷的情況,圖11(a)顯示,測點距離#SN01損傷拉索越近,回歸模型RM(D11,D12)和RM(D12,D13)的γ顯著大于其他模型,且隨著拉索損傷程度的增加,兩個模型的γ逐步增大;對于兩根拉索出現(xiàn)損傷的情況,圖11(b)顯示,測點距離#SN01和#SN17損傷拉索越近,回歸模型RM(D11,D12)、RM(D12,D13)、RM(D15,D16)和RM(D16,D17)的γ顯著大于其他模型,且隨著#SS17拉索損傷程度的增加,回歸模型RM(D15,D16)和RM(D16,D17)的γ逐步增大。據(jù)此可知,在拉索出現(xiàn)損傷時,通過其附近測點溫度撓度構建的回歸模型的參數(shù)的變化率相比未損傷的地方較大,且隨著損傷程度的增加,其回歸模型參數(shù)的變化率也會更大。

    (a)單根拉索出現(xiàn)損傷

    上述分析是從回歸模型參數(shù)變化率角度分析的,以下直觀分析回歸模型參數(shù)的值,其中每天計算一個值,共365個數(shù)值。圖12為#SN01拉索出現(xiàn)不同程度損傷時,RM(D11,D12)和RM(D12,D13)回歸模型參數(shù)a。從圖12可以看出:在理想狀態(tài)下,求得的回歸模型參數(shù)值會在同一種工況下保持恒定,但考慮實際中監(jiān)測值受環(huán)境噪聲及測量誤差的影響,有必要分析本文方法對噪聲的敏感性。圖13為信號信噪比為30 dB時,不同工況下拉索附近測點的回歸模型參數(shù)值。

    從圖13中可以看出:相鄰測點間回歸模型的參數(shù)ai+1,i均隨著時間發(fā)生改變,雖然存在波動,但波動幅度并不大,對此可引入損傷概率的概念[20],對損傷概率的定義如圖14所示。定義如下

    圖14 損傷概率的定義的示意圖Fig.14 Schematic diagram of the definition of damage probability

    (13)

    (14)

    (15)

    式中:xl為健康狀況下?lián)p傷指標概率分布95%置信區(qū)間的上限值;xu為未知工況下?lián)p傷指標概率分布95%置信區(qū)間的上限值;ftest為未知工況損傷指標的概率密度函數(shù);fhealth為健康工況損傷指標的概率密度函數(shù)。

    (a)RM(D11,D12)

    (a)RM(D11,D12)

    據(jù)上述對損傷概率的定義,不同工況下RM(D11,D12)和RM(D12,D13)的回歸模型參數(shù)a的概率密度分布如圖15所示。計算得到對應回歸模型測點附近拉索損傷的概率,如表4所示(將表中每列最大值加粗)。

    (a)RM(D11,D12)

    從表4中可以很明顯看出:在#SN01拉索出現(xiàn)損傷時,回歸模型RM(D11,D12)和RM(D12,D13)中測點附近出現(xiàn)損傷拉索概率遠大于其他模型,且隨著#SN01拉索出現(xiàn)損傷程度的增加,回歸模型RM(D11,D12)和RM(D12,D13)中測點附近出現(xiàn)損傷拉索概率也逐漸增大,測點D12也正是距離#SN01拉索最近的測點,依據(jù)本文方法可以有效識別出拉索損傷的位置和程度。

    表4 不同工況下回歸模型中測點附近的拉索出現(xiàn)損傷的概率Tab.4 Probability of damage to the cable near the measuring point in the regression model under different working conditions

    3 實時監(jiān)測橋梁的拉索損傷診斷

    為能使提出的拉索損傷識別方法應用到在役斜拉橋的拉索損傷識別中,就必須滿足分析的溫度撓度數(shù)據(jù)是斜拉橋在整體升降溫作用下產(chǎn)生的,但主梁截面日溫差撓度和整體日溫差的變化頻率均為一天,即從時間尺度的差異性來分離二者是很難實現(xiàn)的,特別是當結(jié)構的溫度分布難以測量的情況下,分離二者單獨作用的效應值更是不易實現(xiàn)。參考Wang等的研究對枝江斜拉橋主梁溫度撓度與溫度的關系分析可知,在1:00—8:00由于受到太陽光照微弱,枝江斜拉橋主梁溫度撓度受到主梁截面溫差的影響較小,因此其值主要受整體溫差的影響。故本文通過取1:00—8:00撓度數(shù)據(jù),驗證是否只受整體溫差的影響,即驗證是否滿足式(7)。驗證過程的主要內(nèi)容如下:

    (1)首先取背景橋梁主跨跨中附近6個測點(BD11、BD12、BD13、BD34、BD35、BD36)作為分析對象,監(jiān)測時間區(qū)間為2018年01月01日—2018年12月31日的數(shù)據(jù),采樣頻率為1次/h,對于各測點的溫度撓度按照本文所提辦法分離得到。

    (2)對于分離得到的日溫差撓度,按照假定取各個測點在每天1:00—8:00的溫度撓度數(shù)據(jù),如圖16所示為以兩個相鄰測點的溫度撓度分別作為橫坐標和縱坐標繪制的散點圖,溫度撓度的時間范圍為2018年5月—2018年6月,從圖16中可以看出:兩個相鄰測點的溫度撓度具有較好的線性關系。證明了實際橋梁拉索相鄰測點在1:00—8:00的溫度撓度應滿足式(7)。

    (a)測點BD11和測點BD12

    為驗證本文提出拉索損傷診斷方法在實際中應用的有效性,選取實際橋梁測點編號為上游BD11、BD12、BD13及BD14和下游BD34、BD35、BD36及BD37。構建的回歸模型為:RM(BD11,BD12),RM(BD12,BD13),RM(BD13,BD14),RM(BD34,BD35),RM(BD35,BD36),RM(BD36,BD37),以每天1:00—8:00主梁日溫差撓度做一次計算得到每個回歸模型對應的回歸參數(shù)ai+1,i,如圖17所示為每個回歸模型參數(shù)ai+1,i在2018年01月01日—2018年12月31日時間段內(nèi)隨時間的變化情況。

    (a)RM(BD11,BD12)

    各模型參數(shù)ai+1,i概率分布,如圖18所示。2018年人工巡檢及荷載試驗表明該斜拉橋拉索整體運營狀況良好,未出現(xiàn)損傷,所以將前6個月作為損傷識別的健康工況,2018年07月—2018年12月作為損傷識別的未知工況計算得到每個回歸模型的損傷概率,如表5所示。從表5可以看出,損傷概率較小,說明結(jié)構整體完好,與2018年全年檢修結(jié)果相符合,證明了本文方法的正確性。

    (a)RM(BD11,BD12)

    表5 每個回歸模型中測點的附近拉索出現(xiàn)損傷的概率Tab.5 Probability of damage to the cable near the point in each regression model

    4 結(jié) 論

    基于武漢某斜拉橋的撓度監(jiān)測數(shù)據(jù),提出一元線性回歸模型結(jié)合損傷概率的拉索損傷識別方法,具體分析結(jié)果如下:

    (1)計算了沿縱向相鄰測點間主梁溫度撓度的相關系數(shù),可知相鄰測點有極高的相關性;對相鄰測點間進行線性回歸擬合,線性擬合程度較好。

    (2)建立斜拉橋有限元模型,分析了在溫度作用下主梁撓度對拉索損傷的敏感性,分析了相鄰測點相關系數(shù),并構造一元線性回歸模型,從相鄰測點間溫度撓度的回歸參數(shù)ai+1,i、γ(i)對單根拉索出現(xiàn)損傷以及對兩個不同位置的拉索出現(xiàn)損傷進行了準確識別;引入損傷概率的定義考慮噪聲的影響,證明此方法受噪聲影響較小。

    (3)證明了1:00—8:00溫度撓度僅受整體溫差的影響,并對武漢某斜拉橋進行了損傷識別,結(jié)果表明結(jié)構整體完好,與2018年全年檢修結(jié)果相符合,證明了本文方法的正確性。

    猜你喜歡
    斜拉橋拉索撓度
    “拉索”精確測量最亮伽馬暴
    軍事文摘(2024年4期)2024-03-19 09:40:02
    斜拉橋風致振動控制及其未來發(fā)展趨勢
    手緩解拉索優(yōu)化設計
    Spontaneous multivessel coronary artery spasm diagnosed with intravascular ultrasound imaging:A case report
    矮塔斜拉橋彈塑性地震響應分析
    (112+216+112)m部分斜拉橋設計
    上地斜拉橋:天上滴落的水珠
    VOF法在斜拉索風雨激振數(shù)值模擬中的應用
    懸高測量在橋梁撓度快速檢測中的應用
    采用向量式有限元的斜拉索振動控制仿真
    九色国产91popny在线| 91成人精品电影| 欧美最黄视频在线播放免费| 午夜福利视频1000在线观看| 亚洲男人天堂网一区| 天天躁夜夜躁狠狠躁躁| 日日摸夜夜添夜夜添小说| 国内精品久久久久久久电影| 91av网站免费观看| 久久精品aⅴ一区二区三区四区| 国产视频一区二区在线看| 日韩大尺度精品在线看网址| av福利片在线| 亚洲欧美一区二区三区黑人| 男人舔女人下体高潮全视频| 欧美三级亚洲精品| 啦啦啦韩国在线观看视频| 两人在一起打扑克的视频| 香蕉丝袜av| 久久99热这里只有精品18| 免费看a级黄色片| 韩国av一区二区三区四区| 久久 成人 亚洲| www日本黄色视频网| 曰老女人黄片| 手机成人av网站| 少妇熟女aⅴ在线视频| 人人妻人人看人人澡| www日本在线高清视频| 伦理电影免费视频| 美女 人体艺术 gogo| 日韩视频一区二区在线观看| 欧美黄色片欧美黄色片| 亚洲美女黄片视频| 男男h啪啪无遮挡| 国产97色在线日韩免费| 日本撒尿小便嘘嘘汇集6| 日韩一卡2卡3卡4卡2021年| 成人国语在线视频| 日韩一卡2卡3卡4卡2021年| 桃色一区二区三区在线观看| 亚洲人成伊人成综合网2020| 黄色视频,在线免费观看| 亚洲人成77777在线视频| 亚洲精品色激情综合| 成人国产一区最新在线观看| 婷婷精品国产亚洲av| 一本大道久久a久久精品| 看免费av毛片| 精品日产1卡2卡| 视频区欧美日本亚洲| 国产主播在线观看一区二区| 一区二区日韩欧美中文字幕| 大型黄色视频在线免费观看| 熟女电影av网| 一级a爱视频在线免费观看| 国产精品一区二区精品视频观看| 成人欧美大片| 亚洲欧美日韩高清在线视频| 在线播放国产精品三级| 久久这里只有精品19| 亚洲欧美激情综合另类| 久久草成人影院| 国产极品粉嫩免费观看在线| 黄色视频不卡| 亚洲一区二区三区色噜噜| 老熟妇乱子伦视频在线观看| 美国免费a级毛片| 国产伦在线观看视频一区| 亚洲成国产人片在线观看| 国产精品久久久久久人妻精品电影| 国产欧美日韩精品亚洲av| 亚洲国产毛片av蜜桃av| 免费在线观看影片大全网站| 成年版毛片免费区| 在线观看舔阴道视频| 啦啦啦观看免费观看视频高清| 日本在线视频免费播放| 欧美成人免费av一区二区三区| 最新在线观看一区二区三区| 免费看十八禁软件| 免费女性裸体啪啪无遮挡网站| 欧美黑人欧美精品刺激| 国产精品久久久av美女十八| 不卡一级毛片| 欧美性长视频在线观看| 美国免费a级毛片| 看片在线看免费视频| 国产人伦9x9x在线观看| 男女那种视频在线观看| 看免费av毛片| 亚洲精品av麻豆狂野| 日韩大尺度精品在线看网址| 一级作爱视频免费观看| 一级片免费观看大全| 视频在线观看一区二区三区| 俺也久久电影网| 久久性视频一级片| 丝袜美腿诱惑在线| 久久久水蜜桃国产精品网| 日韩精品青青久久久久久| 国产99白浆流出| 不卡一级毛片| av中文乱码字幕在线| 欧美+亚洲+日韩+国产| 一本综合久久免费| 成人欧美大片| 日本在线视频免费播放| 老鸭窝网址在线观看| 国产精品精品国产色婷婷| 久久精品国产亚洲av高清一级| 国产成人影院久久av| 亚洲国产精品成人综合色| 18禁美女被吸乳视频| 日本 欧美在线| cao死你这个sao货| 国产一区二区在线av高清观看| 他把我摸到了高潮在线观看| 国产精品免费一区二区三区在线| 波多野结衣高清作品| 午夜免费观看网址| 久久天堂一区二区三区四区| 国产欧美日韩精品亚洲av| 欧美午夜高清在线| 久久久久久国产a免费观看| 精华霜和精华液先用哪个| av天堂在线播放| 国产一级毛片七仙女欲春2 | ponron亚洲| 精品一区二区三区视频在线观看免费| 日韩中文字幕欧美一区二区| 精品电影一区二区在线| 亚洲三区欧美一区| 黑人操中国人逼视频| 国产精品乱码一区二三区的特点| 国产亚洲精品av在线| 久久久久免费精品人妻一区二区 | 在线十欧美十亚洲十日本专区| 欧美乱码精品一区二区三区| 久久国产亚洲av麻豆专区| 脱女人内裤的视频| 亚洲在线自拍视频| 男女那种视频在线观看| 中文在线观看免费www的网站 | 日韩一卡2卡3卡4卡2021年| bbb黄色大片| 操出白浆在线播放| 真人一进一出gif抽搐免费| 高清毛片免费观看视频网站| 一级片免费观看大全| 一个人观看的视频www高清免费观看 | 一区二区三区高清视频在线| 身体一侧抽搐| 免费在线观看影片大全网站| 成人国产一区最新在线观看| 午夜福利高清视频| 波多野结衣高清无吗| 午夜a级毛片| 午夜福利一区二区在线看| 亚洲精品色激情综合| 日韩视频一区二区在线观看| 国产成人影院久久av| a在线观看视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 一进一出抽搐动态| 国产成人欧美在线观看| 高潮久久久久久久久久久不卡| 超碰成人久久| 一级a爱视频在线免费观看| 日本黄色视频三级网站网址| 欧美日韩黄片免| 老熟妇仑乱视频hdxx| 亚洲欧美日韩无卡精品| 一区二区日韩欧美中文字幕| www.999成人在线观看| 久久人人精品亚洲av| 亚洲国产毛片av蜜桃av| 波多野结衣av一区二区av| 国产av在哪里看| 国产成人系列免费观看| 免费av毛片视频| 亚洲av美国av| 免费看美女性在线毛片视频| 成人免费观看视频高清| 国产精品日韩av在线免费观看| 久久久水蜜桃国产精品网| 国产欧美日韩精品亚洲av| 极品教师在线免费播放| 免费人成视频x8x8入口观看| 后天国语完整版免费观看| 久99久视频精品免费| 啪啪无遮挡十八禁网站| 99国产综合亚洲精品| а√天堂www在线а√下载| 国产成人精品无人区| 国产精品久久久人人做人人爽| 久久99热这里只有精品18| 黄片小视频在线播放| АⅤ资源中文在线天堂| 亚洲av日韩精品久久久久久密| 久久中文字幕一级| 亚洲黑人精品在线| 91国产中文字幕| 欧美色欧美亚洲另类二区| aaaaa片日本免费| 欧美成人免费av一区二区三区| 中出人妻视频一区二区| √禁漫天堂资源中文www| 国产片内射在线| 人妻丰满熟妇av一区二区三区| 真人做人爱边吃奶动态| 女人爽到高潮嗷嗷叫在线视频| 日本三级黄在线观看| 亚洲avbb在线观看| 观看免费一级毛片| 久久 成人 亚洲| 免费在线观看视频国产中文字幕亚洲| 又紧又爽又黄一区二区| 999精品在线视频| 日韩成人在线观看一区二区三区| 中文字幕人成人乱码亚洲影| 久久香蕉国产精品| 色精品久久人妻99蜜桃| av福利片在线| 国产精品永久免费网站| 国产免费男女视频| 欧美激情久久久久久爽电影| 亚洲午夜精品一区,二区,三区| 12—13女人毛片做爰片一| 久久天堂一区二区三区四区| 制服丝袜大香蕉在线| 国产国语露脸激情在线看| 日本在线视频免费播放| 中文资源天堂在线| 欧美乱码精品一区二区三区| 婷婷亚洲欧美| 精品国产亚洲在线| 日本五十路高清| 精品国产美女av久久久久小说| 久久中文字幕人妻熟女| 久久香蕉激情| www.熟女人妻精品国产| 伦理电影免费视频| 观看免费一级毛片| 91在线观看av| 国产精品一区二区精品视频观看| 久久精品人妻少妇| 亚洲第一青青草原| 亚洲七黄色美女视频| 国产成人影院久久av| 亚洲成a人片在线一区二区| 色尼玛亚洲综合影院| 欧美乱码精品一区二区三区| 久久香蕉精品热| 91字幕亚洲| 日韩精品免费视频一区二区三区| 久久香蕉国产精品| av超薄肉色丝袜交足视频| 精品无人区乱码1区二区| 精品日产1卡2卡| 欧美一级毛片孕妇| 丰满人妻熟妇乱又伦精品不卡| 最近最新中文字幕大全电影3 | 久久草成人影院| 最近最新中文字幕大全电影3 | 不卡一级毛片| 国产精品久久久av美女十八| 男女床上黄色一级片免费看| 免费女性裸体啪啪无遮挡网站| 制服诱惑二区| 91麻豆av在线| 看片在线看免费视频| 久久国产精品影院| 巨乳人妻的诱惑在线观看| 99热只有精品国产| 亚洲一码二码三码区别大吗| 国产亚洲欧美98| 91大片在线观看| 欧美 亚洲 国产 日韩一| 国产精品亚洲一级av第二区| 日韩欧美国产一区二区入口| 黄网站色视频无遮挡免费观看| 999久久久国产精品视频| 在线观看免费视频日本深夜| av天堂在线播放| 欧美黑人欧美精品刺激| 最近最新中文字幕大全免费视频| 欧美人与性动交α欧美精品济南到| www国产在线视频色| 国产99久久九九免费精品| 欧美不卡视频在线免费观看 | 中文字幕精品免费在线观看视频| 日韩中文字幕欧美一区二区| 日韩 欧美 亚洲 中文字幕| 亚洲av成人av| 制服诱惑二区| 精品国产乱子伦一区二区三区| 久久国产精品影院| 大型黄色视频在线免费观看| 黄网站色视频无遮挡免费观看| 黄色 视频免费看| 国产主播在线观看一区二区| 欧美zozozo另类| 亚洲欧美精品综合久久99| 男男h啪啪无遮挡| 亚洲国产高清在线一区二区三 | 国产不卡一卡二| 露出奶头的视频| av超薄肉色丝袜交足视频| 国产亚洲欧美98| 精品卡一卡二卡四卡免费| 久久婷婷成人综合色麻豆| 国产免费av片在线观看野外av| 亚洲欧美一区二区三区黑人| 一进一出抽搐动态| 成人18禁在线播放| 亚洲色图av天堂| 无遮挡黄片免费观看| 日本精品一区二区三区蜜桃| 一本久久中文字幕| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 日本五十路高清| 国产亚洲欧美在线一区二区| 18禁裸乳无遮挡免费网站照片 | 亚洲人成网站高清观看| 亚洲专区国产一区二区| 精品人妻1区二区| 久久久国产精品麻豆| 久99久视频精品免费| 国产97色在线日韩免费| 精品卡一卡二卡四卡免费| 日韩免费av在线播放| 欧美日韩中文字幕国产精品一区二区三区| 后天国语完整版免费观看| 久久热在线av| 黑人操中国人逼视频| 麻豆成人av在线观看| 久久久国产欧美日韩av| 成人亚洲精品一区在线观看| 欧美人与性动交α欧美精品济南到| 久久人妻av系列| 国产黄片美女视频| 久久狼人影院| 这个男人来自地球电影免费观看| 亚洲精品国产一区二区精华液| 欧美日韩瑟瑟在线播放| 午夜福利18| 亚洲成av片中文字幕在线观看| 成人手机av| 午夜久久久久精精品| 国产一级毛片七仙女欲春2 | 禁无遮挡网站| 午夜久久久久精精品| 亚洲欧美精品综合一区二区三区| 精品国产一区二区三区四区第35| 欧美激情高清一区二区三区| 黄色片一级片一级黄色片| 啦啦啦韩国在线观看视频| 亚洲国产高清在线一区二区三 | 亚洲精品中文字幕一二三四区| 国产97色在线日韩免费| 国产亚洲精品综合一区在线观看 | 动漫黄色视频在线观看| 免费电影在线观看免费观看| 久久精品成人免费网站| 日本三级黄在线观看| 久久久久国产精品人妻aⅴ院| 亚洲自拍偷在线| 极品教师在线免费播放| 在线看三级毛片| 久久精品国产综合久久久| 亚洲专区中文字幕在线| 日韩三级视频一区二区三区| 高清毛片免费观看视频网站| 制服人妻中文乱码| 又紧又爽又黄一区二区| 国产欧美日韩一区二区三| 色精品久久人妻99蜜桃| 色尼玛亚洲综合影院| 少妇 在线观看| 他把我摸到了高潮在线观看| 视频在线观看一区二区三区| 欧美一级a爱片免费观看看 | 欧美 亚洲 国产 日韩一| 啦啦啦免费观看视频1| 99久久无色码亚洲精品果冻| 国产av一区在线观看免费| 老汉色∧v一级毛片| 亚洲人成77777在线视频| 熟女少妇亚洲综合色aaa.| 国产视频一区二区在线看| 国产精品日韩av在线免费观看| 老司机午夜福利在线观看视频| 看黄色毛片网站| 国产成人系列免费观看| 母亲3免费完整高清在线观看| 国产在线精品亚洲第一网站| 国产精品影院久久| av超薄肉色丝袜交足视频| 啦啦啦观看免费观看视频高清| 亚洲专区中文字幕在线| 午夜亚洲福利在线播放| 色综合欧美亚洲国产小说| 日本熟妇午夜| 久久久国产欧美日韩av| 亚洲精品中文字幕一二三四区| 天堂影院成人在线观看| 少妇粗大呻吟视频| 久久中文字幕一级| 国产片内射在线| 午夜两性在线视频| 久久久久久久久久黄片| 操出白浆在线播放| 在线国产一区二区在线| 波多野结衣av一区二区av| 日韩欧美免费精品| 久久伊人香网站| 免费无遮挡裸体视频| 少妇被粗大的猛进出69影院| 欧美人与性动交α欧美精品济南到| 亚洲av片天天在线观看| 亚洲国产毛片av蜜桃av| 久久精品成人免费网站| 50天的宝宝边吃奶边哭怎么回事| 国产爱豆传媒在线观看 | 国产又爽黄色视频| 母亲3免费完整高清在线观看| 午夜久久久久精精品| 一区二区三区高清视频在线| 淫秽高清视频在线观看| 精品午夜福利视频在线观看一区| 久久香蕉精品热| 宅男免费午夜| 国产精华一区二区三区| 成人精品一区二区免费| 亚洲精品在线观看二区| 1024视频免费在线观看| 亚洲av片天天在线观看| 色在线成人网| 欧美成人一区二区免费高清观看 | av欧美777| 1024手机看黄色片| 国产伦在线观看视频一区| 久久精品国产亚洲av香蕉五月| 真人一进一出gif抽搐免费| 国产色视频综合| 欧美黑人巨大hd| 免费在线观看亚洲国产| 欧美精品啪啪一区二区三区| 91麻豆精品激情在线观看国产| 99在线视频只有这里精品首页| 色综合婷婷激情| 日韩精品免费视频一区二区三区| 一区二区三区精品91| 男男h啪啪无遮挡| 午夜福利欧美成人| 午夜免费激情av| 在线国产一区二区在线| 国产精品香港三级国产av潘金莲| 亚洲五月婷婷丁香| 99久久99久久久精品蜜桃| 美女高潮喷水抽搐中文字幕| 黄色丝袜av网址大全| 91成人精品电影| 精品一区二区三区四区五区乱码| 12—13女人毛片做爰片一| 搡老妇女老女人老熟妇| 高清在线国产一区| 久久精品人妻少妇| 免费搜索国产男女视频| 99精品在免费线老司机午夜| 国产精品久久久av美女十八| 在线观看免费日韩欧美大片| 久久 成人 亚洲| 男女下面进入的视频免费午夜 | 色综合欧美亚洲国产小说| 动漫黄色视频在线观看| 叶爱在线成人免费视频播放| 亚洲欧美日韩高清在线视频| 国产精品香港三级国产av潘金莲| 一级作爱视频免费观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品粉嫩美女一区| 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| 香蕉国产在线看| 手机成人av网站| 叶爱在线成人免费视频播放| 国产精品电影一区二区三区| 色综合亚洲欧美另类图片| 桃色一区二区三区在线观看| 色哟哟哟哟哟哟| 免费看十八禁软件| 国产伦在线观看视频一区| 欧美黄色片欧美黄色片| √禁漫天堂资源中文www| www.自偷自拍.com| 后天国语完整版免费观看| 久久久久国产精品人妻aⅴ院| av有码第一页| 最近最新中文字幕大全电影3 | 日韩欧美一区二区三区在线观看| 99久久国产精品久久久| 欧美不卡视频在线免费观看 | 久久天堂一区二区三区四区| 18禁黄网站禁片午夜丰满| 日本成人三级电影网站| 大型av网站在线播放| 亚洲国产日韩欧美精品在线观看 | 草草在线视频免费看| 麻豆成人午夜福利视频| 在线观看66精品国产| 亚洲七黄色美女视频| 中文资源天堂在线| 19禁男女啪啪无遮挡网站| 午夜成年电影在线免费观看| 女人高潮潮喷娇喘18禁视频| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇 在线观看| 日日干狠狠操夜夜爽| 最近最新中文字幕大全电影3 | 男人舔女人下体高潮全视频| 免费在线观看黄色视频的| 黄色毛片三级朝国网站| 午夜福利欧美成人| 在线观看免费日韩欧美大片| 欧美性长视频在线观看| 亚洲精品国产精品久久久不卡| 午夜成年电影在线免费观看| 亚洲自拍偷在线| 成年人黄色毛片网站| 欧美日韩瑟瑟在线播放| 亚洲无线在线观看| 国产精品一区二区免费欧美| 午夜福利欧美成人| 一进一出抽搐动态| 色哟哟哟哟哟哟| 99久久无色码亚洲精品果冻| 国产亚洲精品av在线| 精品国产乱码久久久久久男人| 欧美性长视频在线观看| 欧美日韩亚洲综合一区二区三区_| 人妻丰满熟妇av一区二区三区| 久热这里只有精品99| 久久这里只有精品19| 久久精品91蜜桃| 视频在线观看一区二区三区| 又黄又爽又免费观看的视频| 久久久久久国产a免费观看| 精品一区二区三区av网在线观看| 色综合亚洲欧美另类图片| 亚洲avbb在线观看| 亚洲第一av免费看| 成年免费大片在线观看| 此物有八面人人有两片| 亚洲人成电影免费在线| 成人特级黄色片久久久久久久| 男女视频在线观看网站免费 | 精品国内亚洲2022精品成人| 欧美黑人精品巨大| 欧美中文综合在线视频| 欧美av亚洲av综合av国产av| 很黄的视频免费| 自线自在国产av| 欧美成狂野欧美在线观看| 一级a爱视频在线免费观看| 叶爱在线成人免费视频播放| 久久 成人 亚洲| 亚洲第一青青草原| 国产精品永久免费网站| 欧美在线黄色| www国产在线视频色| 91成人精品电影| 国产日本99.免费观看| 欧美一区二区精品小视频在线| 午夜精品在线福利| 欧美成人午夜精品| 亚洲精品中文字幕一二三四区| netflix在线观看网站| 国产欧美日韩精品亚洲av| 国产成人av激情在线播放| 欧美黄色片欧美黄色片| 国产aⅴ精品一区二区三区波| 国产精品久久久久久精品电影 | 精品电影一区二区在线| 99精品久久久久人妻精品| 美女免费视频网站| 日韩一卡2卡3卡4卡2021年| 精品人妻1区二区| 十分钟在线观看高清视频www| 人人妻,人人澡人人爽秒播| 成年免费大片在线观看| av天堂在线播放| 十分钟在线观看高清视频www| 波多野结衣av一区二区av| 少妇的丰满在线观看| 久久国产乱子伦精品免费另类| 国产精品亚洲av一区麻豆| 美女免费视频网站| 校园春色视频在线观看| 男人舔女人的私密视频| 黄色片一级片一级黄色片| 好男人电影高清在线观看| 99热6这里只有精品| 国产亚洲精品第一综合不卡| 亚洲第一欧美日韩一区二区三区| 国产又色又爽无遮挡免费看| 在线播放国产精品三级| 深夜精品福利| 亚洲 欧美一区二区三区| 视频区欧美日本亚洲| 精品国产国语对白av| 亚洲无线在线观看| 啦啦啦韩国在线观看视频|