亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        人工免疫算法在橋梁健康監(jiān)測中的應用

        2016-02-09 09:28:28陳偉民
        重慶理工大學學報(自然科學) 2016年12期
        關(guān)鍵詞:活荷載撓度故障診斷

        吳 軍,陳偉民

        (1.西南醫(yī)科大學醫(yī)學信息與工程學院,四川瀘州646000;2.重慶大學光電技術(shù)與系統(tǒng)教育部重點實驗室,重慶400044)

        人工免疫算法在橋梁健康監(jiān)測中的應用

        吳 軍1,陳偉民2

        (1.西南醫(yī)科大學醫(yī)學信息與工程學院,四川瀘州646000;2.重慶大學光電技術(shù)與系統(tǒng)教育部重點實驗室,重慶400044)

        橋梁健康監(jiān)測中故障樣本不易獲取,模型試驗不僅成本高,而且通過模型試驗建立的故障樣本與實際橋梁相差甚遠。為此提出了一種基于人工免疫算法,以撓度、應變等橋梁靜態(tài)數(shù)據(jù)的角度進行分析,并從橋梁健康運行時的長期歷史數(shù)據(jù)中尋找有效信息的數(shù)據(jù)評估方法。經(jīng)過在菜園壩長江大橋一年期的試運行和實踐,證明了基于人工免疫算法的橋梁狀態(tài)評估方法的有效性,其有效率達到了40%。

        橋梁健康監(jiān)測;狀態(tài)評估;人工免疫

        在現(xiàn)代公路、鐵路交通中,橋梁是最重要和常見的基礎(chǔ)設施之一。如果橋梁發(fā)生安全問題,將會造成重大的財產(chǎn)損失及人員傷亡。傳統(tǒng)土木工程行業(yè)對橋梁故障診斷采取的方法是測試動態(tài)數(shù)據(jù),即模態(tài)、振型等[1]。低階模態(tài)對故障反應不夠靈敏,需要采用高階模態(tài)來進行橋梁故障診斷,但高階模態(tài)往往難以提取。所以,從另一個角度來看,可從靜態(tài)數(shù)據(jù)的角度來進行橋梁健康監(jiān)測的數(shù)據(jù)評估及故障診斷。所謂靜態(tài)數(shù)據(jù)是指應變[2-3]、撓度、索力等。靜態(tài)數(shù)據(jù)具有采集方便、處理較為容易等優(yōu)點,因此可采用智能算法、數(shù)據(jù)挖掘、機器學習等進行數(shù)據(jù)評估和故障模式識別。另一個重要的實際問題就是橋梁故障診斷不像機械故障診斷,比較容易獲得故障樣本,橋梁的故障診斷和模式識別中不能從實際橋梁中獲取故障樣本,如果做模型試驗,成本太高,而且模型很難做到與實際橋梁的情況相一致。設計正常的橋梁在大部分時間應該處于正常運行,所以本文提出了一種從橋梁長期運行的數(shù)據(jù)中獲取橋梁正常運行的數(shù)據(jù)模式,從而發(fā)現(xiàn)異常模式的橋梁故障識別方法。

        橋梁健康診斷技術(shù)發(fā)展至今已經(jīng)歷了3個發(fā)展階段:第1階段是以土木工程領(lǐng)域?qū)<业母泄俸蛯I(yè)經(jīng)驗為基礎(chǔ)的經(jīng)驗診斷技術(shù),對診斷信息只能作簡單的數(shù)據(jù)處理;第2階段是以傳感器技術(shù)和動態(tài)測試技術(shù)為手段,以信號處理和建模技術(shù)等為基礎(chǔ)的現(xiàn)代診斷技術(shù),在工程中已得到了廣泛的應用;第3階段為智能診斷技術(shù)階段。近年來,為了滿足大型復雜結(jié)構(gòu)的健康診斷要求,診斷技術(shù)進入了以知識處理為核心,數(shù)據(jù)處理、信號處理與知識處理相融合的智能化數(shù)據(jù)處理階段,智能化故障模式識別和遠程橋梁健康監(jiān)測已經(jīng)成為橋梁故障診斷的發(fā)展趨勢和研究熱點[4]。因此,目前國內(nèi)外已展開大量基于橋梁故障診斷的傳感器技術(shù)、物聯(lián)網(wǎng)和故障識別方法等的研究。

        橋梁健康監(jiān)測和故障診斷的主要特點有:①橋梁結(jié)構(gòu)復雜,所處環(huán)境惡劣;②橋梁健康監(jiān)測系統(tǒng)是一套長期的、在線的監(jiān)測系統(tǒng),包含多種傳感器[5-9],每種傳感器數(shù)量也很多,分別布設在不同的自由度;③必須實現(xiàn)無損檢測,即橋梁健康監(jiān)測系統(tǒng)不能對結(jié)構(gòu)本身性能產(chǎn)生任何不良影響;④屬于多學科領(lǐng)域交叉的研究。橋梁故障診斷靜態(tài)數(shù)據(jù)識別方法必須考慮臺風、地震、車載及溫度等環(huán)境參數(shù)變化對橋梁結(jié)構(gòu)效應的影響,因此在進行橋梁數(shù)據(jù)評估之前需要先剔除這些噪聲和干擾信息。本研究從長期的歷史信息和數(shù)據(jù)中,挖掘橋梁正常運行時的健康模式樣本,即把橋梁健康監(jiān)測中采集到的靜態(tài)數(shù)據(jù)視為時間序列,從長期正常運行的數(shù)據(jù)中進行統(tǒng)計學習和分析,不需要建立故障模式樣本或?qū)<規(guī)斓?。如果出現(xiàn)異常數(shù)據(jù),即給出預警信息。故本方法屬于上述的第3階段橋梁故障診斷方法。根據(jù)所提出的需求,最近十幾年來出現(xiàn)的一種新的智能仿生算法——人工免疫算法比較符合本文思想和方法的實際需要[10-12]。因此,提出了一種基于人工免疫算法的新的橋梁健康監(jiān)測數(shù)據(jù)評估方法。

        1 橋梁健康監(jiān)測中的兩類錯誤

        在橋梁健康監(jiān)測中存在兩類錯誤,如下所示:

        第1類錯誤:當橋梁出現(xiàn)劣化的時候,應該預警的時候而沒有給出預警,它的概率如式(1)所示。

        式中:符號“|”表示條件概率;A表示橋梁健康監(jiān)測系統(tǒng)給出預警;B則表示橋梁已經(jīng)發(fā)生了劣化。

        第2類錯誤:當橋梁在正常運行的狀態(tài)下,健康監(jiān)測系統(tǒng)給出虛假預警信號的概率,如式(2)所示。

        式(1)表示橋梁已經(jīng)出現(xiàn)故障,監(jiān)測系統(tǒng)卻未給出預警信號的概率;式(2)則表示橋梁還未出現(xiàn)故障,監(jiān)測系統(tǒng)卻給出預警信號的概率。由此可知,第2類錯誤可適當容忍,第1類錯誤卻需要盡量避免。兩類錯誤屬于此消彼長的關(guān)系,即降低第1類錯誤的概率,就會增加第2類錯誤的概率;反之亦然。因此,本文橋梁健康監(jiān)測系統(tǒng)故障診斷的目標,就是在控制第1類錯誤的同時,盡量降低第2類錯誤的概率。

        2 橋梁健康監(jiān)測數(shù)據(jù)的預處理

        在橋梁長期運行的監(jiān)測數(shù)據(jù)中,包含了活荷載效應、溫差效應、風載效應、劣化效應和噪聲?;詈奢d效應是指由車流所造成的效應。對于溫差效應,本文只考慮日溫差變化所造成的效應,而年溫差不予考慮(因為對橋梁來說,未達到滿載運行的情況下,日溫差將是主要的荷載因素)。噪聲包含傳感器的測量噪聲、電噪聲等。要對橋梁的靜態(tài)數(shù)據(jù)進行健康評估,就需要事先剔除無用的活荷載信息(隨機信息與車流有關(guān))和周期性的溫度信息(隨機信息與氣溫有關(guān))。

        本文以重慶菜園壩長江大橋的健康監(jiān)測和故障診斷為例進行分析。菜園壩長江大橋全橋跨徑為800 m,是全國最早采用健康監(jiān)測系統(tǒng)的橋梁之一,在橋上布置有撓度傳感器、索力傳感器、應變傳感器等。菜園壩長江大橋?qū)儆诠皹颍謨蓪?,上層為公路,下層為重慶軌道交通3號線。圖1中以其中跨跨中節(jié)點為例,給出了撓度數(shù)據(jù)的活荷載效應提取結(jié)果。因為活荷載作用的時間尺度較短,因此本文采用了多尺度分析的方法剔除橋梁活荷載效應[13]。

        圖1采用了集合經(jīng)驗模態(tài)分解法(Ensemble Empirical Mode Decomposition,EEMD)剔除橋梁活荷載效應,經(jīng)過對比分析和實際驗證發(fā)現(xiàn),采用EEMD法比采用小波分析法效果更好,活荷載效應剔除得更“干凈”。圖2則給出了中跨跨中節(jié)點的日溫差效應提取及剔除日溫差效應后的結(jié)果。由于橋梁的日溫差效應具有與日氣溫變化同樣的規(guī)律,因此本文采用了以一天24 h的時間為窗口,進行加窗傅里葉變換,剔除掉與日溫差變化頻率相同的信息。

        圖2中給出的剩余荷載效應即為剔除了活荷載效應、日溫差效應后的預處理結(jié)果,即“殘余信息”。這個“殘余信息”中主要包含了劣化效應和各種噪聲信號,噪聲屬于白噪聲,從長期運行數(shù)據(jù)的結(jié)果來看無關(guān)緊要,因此最重要的是劣化效應。劣化效應包含了橋梁健康狀況劣化及承載力下降的信息,因此本文將以該預處理后的“殘余信息”作為人工免疫算法分析和處理的對象。

        圖1 菜園壩長江大橋跨中節(jié)點的撓度活荷載效應提取

        圖2 菜園壩長江大橋跨中節(jié)點的撓度日溫差效應提取

        3 人工免疫算法及其在橋梁健康監(jiān)測中的應用

        3.1 人工免疫算法簡介

        自從1991年比利時學者Bersini和法國學者Varela首次使用人工免疫算法來解決問題[14-15]以來,人工免疫算法已經(jīng)得到了飛速發(fā)展,國外大批學者投入了該領(lǐng)域的研究。在最近的十年里,人工免疫算法在國內(nèi)也逐漸成為研究熱點之一[16]。

        人工免疫算法是受人體免疫系統(tǒng)的體細胞理論和網(wǎng)絡理論啟發(fā)而提出的一種仿生智能計算方法,它通過模擬人體免疫系統(tǒng)的自然防御機理和學習機制,提供了噪聲容忍、無導師學習、自組織、記憶等進化學習機理,融合了分類器、神經(jīng)網(wǎng)絡和機器推理等許多智能算法的優(yōu)點。生命體的免疫系統(tǒng)是一個自組織、自調(diào)節(jié)的控制系統(tǒng),人工免疫算法就是模擬其功能的一種智能仿生算法。人工免疫算法是一種自適應算法,對異常數(shù)據(jù)的識別能力較強,識別率也較高。因此,在橋梁健康監(jiān)測中,人工免疫算法可以保證在控制第1類錯誤的同時盡量降低第2類錯誤。

        3.2 抗原和抗體

        抗原:是指能刺激和誘導機體的免疫系統(tǒng)使其產(chǎn)生免疫應答,并能與相應的免疫應答產(chǎn)物在體內(nèi)或體外發(fā)生特異性反應的物質(zhì)。在人工免疫系統(tǒng)中,一般指問題及其約束,與進化算法中與適應度函數(shù)類似。具體地是問題目標的函數(shù),是人工免疫系統(tǒng)算法的始動因子及重要的度量標準[17]。

        抗體:是指免疫系統(tǒng)受抗原刺激后,免疫細胞轉(zhuǎn)化為漿細胞并產(chǎn)生能與抗原發(fā)生特異性結(jié)合的免疫球蛋白,該免疫球蛋白即為抗體。在人工免疫系統(tǒng)中一般指問題的候選解,與進化算法中的個體相似,抗體的集合稱為抗體群。在實踐中,一般抗體以編碼的形式出現(xiàn),常用的編碼形式有二進制和十進制。

        抗體和抗原的親和力:為兩個字符串(或向量)之間的距離,如歐氏距離等[18]。

        3.3 免疫和免疫系統(tǒng)

        免疫(immune)是指生物體的一種特異性生理反應,它能夠識別并清除抗原性物質(zhì),“通過免疫可以使機體維持在健康的免疫平衡狀態(tài)”。

        免疫系統(tǒng)(immune system)是機體產(chǎn)生免疫應答反應的物質(zhì)基礎(chǔ),是生物體執(zhí)行免疫功能的組織機構(gòu)。各種淋巴細胞、T細胞和B細胞等各種免疫細胞共同組成了免疫系統(tǒng)。

        免疫應答的基本過程可分為3個階段。首先是感應階段,抗原信息通過提成細胞交給淋巴細胞,T細胞和B細胞通過其表面抗原識受體進行識別。然后是反應階段,淋巴細胞識別后開始分化、增殖,其中B細胞將產(chǎn)生記憶細胞,以便后期能夠迅速識別相同抗原。最后是效應階段,T細胞直接對抗體進行攻擊,B細胞合成、分泌大量抗體并清除抗原。從免疫系統(tǒng)應答的過程可以看出,免疫系統(tǒng)包括初次應答、再次應答和交叉應答等方式。

        3.4 特征向量的提取

        撓度傳感器采用光電液位傳感器,而應變則將光纖應變傳感器預埋進混凝土中獲得。撓度傳感器布置在主跨(中跨)跨中、主跨1/4節(jié)點處及邊跨跨中等位置,本文的數(shù)據(jù)分析均以主跨跨中處的傳感器為例。在每個監(jiān)測節(jié)點處均配置有撓度傳感器和應變傳感器,應變傳感器設置在上下表面處各兩個,上下應變均取這兩個應變傳感器的平均值。本文的監(jiān)測數(shù)據(jù)均為撓度、應變等靜態(tài)數(shù)據(jù),且2 h采集一次,采集數(shù)據(jù)按時間順序排列成隨機時間序列。本文取差分系數(shù)作為評價指標,差分系數(shù)定義如式(3)所示。ii-

        本文分別用字母s1、p、s2代表上應變、撓度、下應變。于是ti時刻的3個差分系數(shù)構(gòu)成向量為

        式(4)所示的向量就構(gòu)成了橋梁健康監(jiān)測中靜態(tài)數(shù)據(jù)的特征向量。

        3.5 人工免疫算法在橋梁健康監(jiān)測中的應用

        一般的機械故障診斷,需要提前建立故障樣本庫,然后進行模式匹配。因為橋梁造價較高,不能進行破壞試驗,所以橋梁的故障樣本不易獲取。而人工免疫算法可以在實時運行中自動識別橋梁撓度和應變的粗大值或異常值等。正常橋梁的撓度變化范圍一般與跨度有關(guān),跨度越大的橋梁撓度值也就越大。所以在土木工程中衡量撓度變化范圍一般不是一個具體的數(shù)值,而是一個百分數(shù)(撓度和跨度的比值)。本文對撓度和應變均采用了差分值來進行計算,在程序中所設定的閾值為10‰(該數(shù)值根據(jù)歷史數(shù)據(jù)的經(jīng)驗得出)。

        本文先采用K-鄰域均值法篩選出最初的核心抗體群(歷史信息中橋梁健康正常運行時的數(shù)據(jù)樣本),然后把采集到的新數(shù)據(jù)作為抗原加入到健康監(jiān)測系統(tǒng)中進行識別。算法設計如圖3所示。

        圖3中的“親和突變”是為了能更好地尋找全局最優(yōu)解,以免陷入局部最優(yōu)解。新采集到的數(shù)據(jù)經(jīng)過有限次迭代之后,如果能匹配將加入到新的抗體群,再用它們來識別新的抗原。如此循環(huán),直至不能匹配為止,將給出警報信號。

        圖3 橋梁健康監(jiān)測中的人工免疫算法設計

        3.6 監(jiān)測結(jié)果

        經(jīng)過一年時間的測試運行,診斷結(jié)果如表1所示。

        表1 人工免疫算法在橋梁健康監(jiān)測中的診斷結(jié)果

        4 結(jié)論和展望

        從表1中的結(jié)果可得:采用人工免疫算法,從橋梁長期正常運行的歷史數(shù)據(jù)中去尋找有用信息是可行的,而且是一種值得深入研究的橋梁健康監(jiān)測數(shù)據(jù)評估方法。但其有效性是否具有統(tǒng)計學意義,還需要在其他橋梁中做進一步的研究和觀察。主要問題是傳感器的壽命與橋梁設計使用年限相比相差太遠,橋梁遠未出現(xiàn)問題時,傳感器已經(jīng)發(fā)生了故障。因此,研究如何提高傳感器的使用壽命是一個非常重要的課題。

        橋梁健康監(jiān)測中的方法都是從機械故障診斷中移植過來的,雖現(xiàn)在機械故障診斷中已經(jīng)得到了一些有效的應用,但由于橋梁的復雜性高、結(jié)構(gòu)較大等特點和所處自然環(huán)境變化較多等實際情況,橋梁健康監(jiān)測也有許多自己的獨特之處,需要根據(jù)實際情況具體進行處理,并對各種診斷算法進行改進才能應用到橋梁健康監(jiān)測中。

        橋梁健康監(jiān)測系統(tǒng)作為人工檢查和維護的輔助手段是有效的,雖然目前橋梁健康監(jiān)測系統(tǒng)還不能達到完全不需人工檢查的地步,還僅僅是土木專業(yè)人工檢查的輔助手段,但橋梁健康監(jiān)測還是一個值得研究的課題和方向,不過還需進一步研究如何提高其有效性和準確性。

        因為人工免疫算法是一種仿生智能算法,所以其收斂性是可以保證的,該算法對異常數(shù)據(jù)也較為敏感。因此從理論上分析和一年的試運行性結(jié)果來看,確實能夠做到在控制第1類錯誤的同時降低第2類錯誤發(fā)生的概率,即漏報率比誤報率要低得多。但是,兩類錯誤概率的理論值還需要進一步進行分析和計算,并需長期觀察實際運行的結(jié)果與理論值是否能夠保持一致。

        [1] 葛玉梅,邵帥,潘輝,等.泰州大橋懸索橋南錨碇基礎(chǔ)的變形及超載安全度分析[J].科技導報,2013,31(26):35-39.

        [2] 孫陽陽,王源,章征林,等.表面粘貼式光纖布拉格光柵應變傳遞規(guī)律分析與實驗研究?[J].功能材料,2016,47(7):07046-07050.

        [3] 王源,章征林,孫陽陽,等.不同表面粘貼方式對光纖光柵應變傳遞的影響[J].壓電與聲光,2016(1):106-110.

        [4] 董聰,楊培捷.城市橋梁智能健康診斷研究與進展[J].科技導報,1999,19(10):57-59,60.

        [5] 盧凱,黃文,劉思祎,等.基于PVDF的柔性壓力傳感器陣列的制備及仿真研究[J].電子元件與材料,2016,35(3):40-43.

        [6] 母坤,童杏林,胡畔,等.氫氣傳感器的技術(shù)現(xiàn)狀及發(fā)展趨勢[J].激光雜志,2016(5):1-5.

        [7] 劉林鋒,郭平,趙娟,等.無線傳感器網(wǎng)絡中基于移動sink最優(yōu)路徑的數(shù)據(jù)收集策略[J].四川兵工學報,2015,36(1):118-121.

        [8] 李梅菊.水下無線傳感器網(wǎng)絡綜述[J].重慶理工大學學報(自然科學),2016(8):92-98,121.

        [9] 唐海松,陳亦新,倪立新.現(xiàn)代轉(zhuǎn)矩轉(zhuǎn)速傳感器在電機能效測試應用中的技術(shù)[J].電機與控制應用,2016,43(1):95-100.

        [10]DIAO Y X,PASSION K M.Fault diagnosis for a turbine engine[J].Control Engineering Practice,2004(12):1151-1165.

        [11]竇唯,劉樹林,孫明,等.生物免疫機理在往復壓縮機在線狀態(tài)監(jiān)測中的應用[J].流體機械,2004,32(5):16-19.

        [12]XIONG H,SUN C X.Artificial Immune Network Classification Algorithm for Fault Diagnosis of Power Transformer[J].IEEE Transactions on Power Delivery,2007,22(2):930-935.

        [13]WU J,CHENW M.An Extraction Method for Live-Load Effect of Bridge Based on EEMD[J].Applied Mechanics&Materials,2011,105/107:1144-1148.

        [14]BERSINIH,VARELA F J.Hints for adaptive problem solving gleaned from immune networks[M]//Parallel Problem Solving from Nature.2006:343-354.

        [15]BERSINIH.Immune network and adaptive control[C]//Proceedings of the 1stEuropean Conference on Artificial Life.1991:217-226.

        [16]欒林林.基于人工免疫算法的故障診斷方法研究[D].上海:華東理工大學,2011:86-88.

        [17]岑健.基于人工免疫系統(tǒng)的機組復合故障診斷技術(shù)研究[D].廣州:華南理工大學,2010:78-92.

        [18]段海濱,張祥銀,徐春芳.仿生智能計算[M].北京:科學出版社,2011:168-188.

        (責任編輯楊文青)

        Application of Artificial Immune Algorithm in Bridge Health Monitoring

        WU Jun1,CHENWei-min2
        (1.College of Medical Informatics Engineering,Southwest Medical University,Luzhou 646000,China;2.Key Lab of Ministry of Education for Optoelectronic Technology&System,Chongqing University,Chongqing 400044,China)

        The fault samples in bridge health monitoring are difficult for acquisition.The cost of model test of bridge is very high,and the fault samples which are got from model test are far from the actual bridge.For these factors,a method of data assessment based on artificial immune is proposed,which analyzes from the perspective of static data,such as deflection,strain,and so on.The historical information is very effective when bridge is running healthfully for a long time.In a year of commissioning and trial in Caiyuanba Yangtz River Bridge,it is proved that the method of bridge state assessment based on artificial immune algorithm is effective and feasible.And the effective rate reached 40%.

        bridge health monitoring;state assessment;artificial immunity

        TP39

        A

        1674-8425(2016)12-0129-06

        10.3969/j.issn.1674-8425(z).2016.12.020

        2016-05-25

        遵義醫(yī)學院博士啟動資金資助(F-741)

        吳軍(1977—),男,博士,副教授,主要從事智能儀器與測試、數(shù)據(jù)分析與大數(shù)據(jù)處理等方面研究,E-mail:wj2135187@126.com。

        吳軍,陳偉民.人工免疫算法在橋梁健康監(jiān)測中的應用[J].重慶理工大學學報(自然科學),2016(12):129-134.

        format:WU Jun,CHEN Wei-min.Application of Artificial Immune Algorithm in Bridge Health Monitoring[J].Journal of Chongqing University of Technology(Natural Science),2016(12):129-134.

        猜你喜歡
        活荷載撓度故障診斷
        活荷載
        北方建筑(2022年2期)2022-11-21 14:57:16
        Spontaneous multivessel coronary artery spasm diagnosed with intravascular ultrasound imaging:A case report
        火電廠主廠房活荷載折減系數(shù)取用分析
        因果圖定性分析法及其在故障診斷中的應用
        PKPM關(guān)于梁側(cè)面水平荷載及活荷載最不利布置的計算結(jié)果分析與思考
        河南建材(2014年5期)2014-03-20 03:29:58
        懸高測量在橋梁撓度快速檢測中的應用
        基于LCD和排列熵的滾動軸承故障診斷
        基于WPD-HHT的滾動軸承故障診斷
        機械與電子(2014年1期)2014-02-28 02:07:31
        高速泵的故障診斷
        河南科技(2014年3期)2014-02-27 14:05:48
        淺談研發(fā)廠房樓面均布活荷載及可變荷載地震組合值系數(shù)的取值
        97精品国产高清自在线看超| 亚洲精彩av大片在线观看| 亚洲一区二区三区码精品色| 国产av在线观看一区二区三区| 亚洲欧洲成人a∨在线观看| 久久伊人少妇熟女大香线蕉| 日韩精品区一区二区三vr| 老司机在线精品视频网站| 欧美 亚洲 国产 日韩 综AⅤ | 国产一区二区三区护士| 大肉大捧一进一出好爽视频动漫| 艳z门照片无码av| 欲妇荡岳丰满少妇岳| 日本一本久道| 三级黄色片一区二区三区 | 日本视频在线观看二区| 国产乱人伦av在线麻豆a| 男人j进女人j啪啪无遮挡| 国产丝袜在线精品丝袜不卡| 亚洲xx视频| 网站在线观看视频一区二区| 性色欲情网站| 国产亚洲日韩欧美一区二区三区 | 闺蜜张开腿让我爽了一夜| 国产精品麻花传媒二三区别 | 无码丰满熟妇一区二区| 越南女子杂交内射bbwbbw| 国内精品福利在线视频| 国产亚洲专区一区二区| 欧美成人看片一区二区三区尤物| 久久婷婷色综合一区二区| 亚洲免费成年女性毛视频| 精品少妇一区二区av免费观看| 高潮抽搐潮喷毛片在线播放| 香蕉久久人人97超碰caoproen| 在线视频一区二区亚洲| 一本久道竹内纱里奈中文字幕| 色费女人18毛片a级毛片视频| 肉体裸交丰满丰满少妇在线观看| 免费va国产高清不卡大片 | 亚洲欧美另类精品久久久|