999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于狀態(tài)的裝備故障預(yù)測技術(shù)綜述

2019-12-01 10:18:00徐宗昌張文俊
關(guān)鍵詞:故障方法模型

郭 建,徐宗昌,張文俊

(1.中國人民解放軍96853部隊,遼寧 沈陽 110000;2.陸軍裝甲兵學(xué)院 裝備保障與再制造系,北京 100072)

隨著監(jiān)測技術(shù)和信號處理技術(shù)的不斷發(fā)展,獲取裝備健康狀態(tài)技術(shù)的逐漸提高,加之對復(fù)雜裝備和重要裝備可靠性、安全性和經(jīng)濟性的要求不斷提高,視情維修(Condition-based Maintenance,CBM)已逐漸成為可能,其成功應(yīng)用不僅提高了裝備可靠性,減少意外停機,而且大大地減少了維修費用。

視情維修作為預(yù)防性維修的高級形式,其成功應(yīng)用改變了以往僅僅依靠以事件為主導(dǎo)的修復(fù)性維修和以時間為主導(dǎo)的定時維修的方式,使獲得裝備狀態(tài)實施準確的維修活動成為可能。視情維修是在裝備狀態(tài)監(jiān)測的基礎(chǔ)上,對裝備進行故障預(yù)測并依據(jù)結(jié)果制定合理的維修計劃,并主動實施的維修。這樣可以最大化地減少裝備維修的費用、增加裝備可用度、提高裝備的可靠性。

1 裝備故障預(yù)測的主要方法

復(fù)雜技術(shù)的應(yīng)用帶來了裝備的結(jié)構(gòu)復(fù)雜化、功能集成化以及技術(shù)綜合化,與此同時,也突顯出裝備故障的多樣化、致命化、隨機化等特點,人與裝備故障的長期斗爭中積累了大量裝備維修方面的實踐和經(jīng)驗。隨著裝備維修水平的不斷提高,人們不僅希望能夠在裝備發(fā)生故障之后能夠進行快速的診斷和制定合理的預(yù)防性維修計劃,還更加迫切地需要監(jiān)測裝備的狀態(tài),預(yù)測裝備的剩余壽命,掌握裝備的健康情況。裝備故障預(yù)測技術(shù)能夠準確的預(yù)測裝備的剩余壽命,既是對狀態(tài)監(jiān)控數(shù)據(jù)有效的使用,又可以作為制定維修決策的重要依據(jù),通過采取及時有效的維修措施來提高裝備的可靠性,減少發(fā)生故障導(dǎo)致的意外停機,是提高裝備使用可用度和減少壽命周期費用的有效手段,因此大力開展裝備故障預(yù)測研究有著非常重要的意義。可用于故障預(yù)測的方法有很多,具體來說可以分為基于傳統(tǒng)可靠性的預(yù)測方法、基于數(shù)據(jù)驅(qū)動的預(yù)測方法、基于統(tǒng)計的預(yù)測方法和基于失效物理的預(yù)測方法。

1.1 基于傳統(tǒng)可靠性的故障預(yù)測方法

基于傳統(tǒng)可靠性的預(yù)測方法主要可分為基于故障樹分析和基于壽命分布模型兩種。

基于故障樹分析的故障預(yù)測方法是通過建立故障樹模型,進行故障樹定性和定量分析來達到故障預(yù)測的目的,對于無故障征兆的情況,根據(jù)總體壽命期望值進行預(yù)測;對于有故障征兆的情況,通過分析系統(tǒng)內(nèi)部關(guān)系來確定影響因素和容易導(dǎo)致的故障。陳濤等[1]通過研究風電齒輪箱傳動系統(tǒng),將故障樹模型分解為幾個單元,并利用灰色關(guān)聯(lián)理論建立可靠性聯(lián)系,提高了故障樹分析在預(yù)測中的準確度。

基于壽命分布模型的預(yù)測是根據(jù)同類裝備壽命服從的分布來對故障時間進行預(yù)測,常見壽命分布有指數(shù)分布、威布爾分布、正態(tài)分布和對數(shù)正態(tài)分布。這類方法的優(yōu)點是不需要對裝備進行狀態(tài)監(jiān)測,以及可以對同類設(shè)備的使用壽命進行大概的估計,缺點是只能針對同類裝備進行預(yù)測,對于單個裝備無法進行準確的預(yù)測。以往開展的很多關(guān)于維修策略的文獻中,假設(shè)存在退化趨勢的機電產(chǎn)品的壽命分布服從威布爾分布,電子產(chǎn)品的壽命服從指數(shù)分布。

1.2 基于數(shù)據(jù)驅(qū)動的故障預(yù)測方法

基于數(shù)據(jù)驅(qū)動的故障預(yù)測方法優(yōu)點是不需要建立物理模型,無需假設(shè)參數(shù)和經(jīng)驗性估計,通常來說復(fù)雜系統(tǒng)建立物理模型是相當困難的,建立數(shù)學(xué)模型進行故障預(yù)測相對簡單,但是通常需要大量準確的歷史數(shù)據(jù)。

典型的方法有基于人工神經(jīng)網(wǎng)絡(luò)的預(yù)測方法,它不僅能用于故障診斷,還能很好地應(yīng)用于故障預(yù)測,經(jīng)過長期的發(fā)展與應(yīng)用,人工神經(jīng)網(wǎng)絡(luò)的主要類型有BP神經(jīng)網(wǎng)絡(luò)、RBF神經(jīng)網(wǎng)絡(luò)和小波神經(jīng)網(wǎng)絡(luò)等。其主要的優(yōu)點在于它具有其他常規(guī)算法和專家系統(tǒng)不具備的自學(xué)習和自適應(yīng)功能,能夠處理多變量分析,提供非線性預(yù)測,并且不需要先驗知識。神經(jīng)網(wǎng)絡(luò)經(jīng)常與其他方法結(jié)合使用,可以取得更好的預(yù)測效果,黨香俊等[2]研究了基于BP神經(jīng)網(wǎng)絡(luò)和小波神經(jīng)網(wǎng)絡(luò)結(jié)合對產(chǎn)品退化軌跡進行預(yù)測,不僅實現(xiàn)了較長距離的預(yù)測,而且提高了預(yù)測的精度。

基于濾波器的預(yù)測方法也是一種基于數(shù)據(jù)驅(qū)動的故障預(yù)測方法,常見的有卡爾曼濾波器、擴展式卡爾曼濾波器、強跟蹤濾波器和粒子濾波器。卡爾曼濾波器是最佳的線性估計器,可實現(xiàn)極小化系統(tǒng)狀態(tài)估計誤差,存在的不足是需要精確的數(shù)學(xué)模型并且對模型不確定性的魯棒性較差。強跟蹤濾波器對模型的不確定性的魯棒性較強,且對緩變和突變狀態(tài)都有極強的跟蹤能力,但是只能用于解決高斯噪聲系統(tǒng)的故障預(yù)測。粒子濾波器與強跟蹤濾波器相反,其具有解決非高斯系統(tǒng)的濾波和預(yù)測的能力,但是對于突變狀態(tài)的跟蹤能力較差。杜黨波等[3]研究了在復(fù)雜系統(tǒng)存在缺失數(shù)據(jù)時,如何采用改進的小波-卡爾曼濾波解決故障預(yù)測問題,并進行了數(shù)值仿真和實例驗證,證明了算法的有效性和可行性。

1.3 基于統(tǒng)計的故障預(yù)測方法

基于統(tǒng)計的預(yù)測方法的最大優(yōu)點在于可以通過觀察到的統(tǒng)計數(shù)據(jù)建立所需的概率密度函數(shù),這些概率密度函數(shù)能夠給出足夠的預(yù)測結(jié)果置信區(qū)間,主要包括貝葉斯網(wǎng)絡(luò)、隱馬爾可夫和隱半馬爾可夫模型、回歸分析方法等。

貝葉斯網(wǎng)絡(luò)是利用對象個體的狀態(tài)監(jiān)控數(shù)據(jù)代替事件數(shù)據(jù),進行可靠性評估,準確性強烈依賴對各種趨勢特征閾值的準確確定。蔡志強等[4]將故障模式、影響及危害性分析(FMECA)與貝葉斯網(wǎng)絡(luò)結(jié)合用于故障預(yù)測,經(jīng)過實驗表明其在不確定性和定量分析方面具有優(yōu)勢,能夠準確地預(yù)測復(fù)雜裝備的故障。

隱馬爾可夫模型(HMM)是由馬爾可夫模型發(fā)展而來的,是一種描述隨機過程統(tǒng)計特性的概率模型,與馬爾可夫模型不同,HMM是雙重隨機過程,一個隨機過程描述狀態(tài)轉(zhuǎn)移,另一個描述狀態(tài)與觀察值之間的統(tǒng)計對應(yīng)關(guān)系。隱馬爾可夫模型根據(jù)觀察信號性質(zhì)不同,可分為連續(xù)隱馬爾可夫模型和離散隱馬爾可夫模型,隱半馬爾可夫模型(HSMM)是HMM的一種擴展形式,介于連續(xù)和離散HMM的半連續(xù)HMM。馮輔周等[5]為解決微弱信號特征提取問題和故障預(yù)測精度不高的問題,提出一種將小波相關(guān)濾波和排列熵結(jié)合的特征提取方法,并將HMM用于故障預(yù)測,最后通過軸承數(shù)據(jù)驗證了算法的有效性和優(yōu)越性。徐宇亮等[6]為解決電子設(shè)備的故障預(yù)測問題,提出一種改進流形學(xué)習與HSMM結(jié)合的健康評估與故障預(yù)測方法,并應(yīng)用與某型導(dǎo)彈電子設(shè)備,驗證了算法的有效性。

回歸分析方法的主要步驟是根據(jù)預(yù)測目標確定自變量和因變量,建立回歸預(yù)測模型,進行相關(guān)分析,模型檢驗與修改,計算預(yù)測值。可以分為一元線性回歸分析預(yù)測法、多元線性回歸分析預(yù)測法和非線性回歸分析預(yù)測法。其主要特點是預(yù)測過程簡單、技術(shù)成熟,存在的不足之處是誤差較大,外推特性差,且要求大量樣本數(shù)據(jù)。鄧力等[7]提出一種將網(wǎng)絡(luò)中單組件的故障率多項式轉(zhuǎn)換成多元線性回歸模型進行求解,給出的預(yù)測區(qū)間具有較高的預(yù)測精度。另外,基于時間序列也是一種重要的預(yù)測方法,這種方法具有計算方便快捷,能夠提供線性平穩(wěn)預(yù)測的優(yōu)點,但是在解決非線性問題時受到局限,而且模型識別需要確定很多參數(shù),常用的時間序列模型有自回歸(AR)模型、滑動平均(MA)模型和自回歸滑動平均(ARMA)模型,李向前[8]提出基于自回歸滑動平均模型與人工神經(jīng)網(wǎng)絡(luò)的狀態(tài)預(yù)測方法,克服了自回歸滑動平均模型無法解決復(fù)雜非線性問題預(yù)測的局限。杜占龍等[9]將卡爾曼濾波與自相關(guān)模型融合用于故障預(yù)測,提出帶有次優(yōu)漸消因子的強跟蹤平方根容積卡爾曼濾波方法,解決了非線性復(fù)雜系統(tǒng)參數(shù)不易測量的問題。

1.4 基于失效物理的故障預(yù)測方法

基于失效物理的預(yù)測方法是利用裝備壽命周期的載荷和故障的失效機理來評估裝備的可靠性以及預(yù)測故障,其優(yōu)勢在于需要的數(shù)據(jù)量少。基于失效物理的預(yù)測方法主要有基于帕里斯法則裂紋擴展建模、基于Forman規(guī)律的裂紋擴展建模、疲勞剝落擴展模型和基于剛度的損傷規(guī)律模型等。

基于帕里斯法則裂紋擴展建模,利用有限元分析使得基于部件幾何特征、缺陷、載荷與速度等材料應(yīng)力計算成為可能,但是計算的成本過高。鄭戰(zhàn)光等[10]推導(dǎo)了一種基于損傷力學(xué)的一般損傷演化過程,并結(jié)合帕里斯公式,以有效當量應(yīng)力幅和材料損傷特性參數(shù)為控制變量得到了新的疲勞損傷演化模型。

基于Forman規(guī)律的裂紋擴展建模,能夠?qū)⒉考臓顟B(tài)監(jiān)測數(shù)據(jù)和缺陷擴展物理特性與壽命關(guān)聯(lián)起來,但是該方法的簡化假設(shè)需要驗證,對于復(fù)雜狀況,模型中的參數(shù)仍需確定。劉文光等[11]將廣義的Forman方程模擬疲勞裂紋擴展應(yīng)用到懸梁臂表面裂紋上,研究了裂紋面摩擦效應(yīng)對裂紋疲勞擴展的影響。薛飛等[12]使用簡化的四參數(shù)全范圍Forman模型表征國產(chǎn)壓水堆核電站主管道母材與焊材全范圍的疲勞擴展規(guī)律。王愛紅等[13]人運用Forman公式估算橋式起重機的疲勞剩余壽命,為橋式起重機的工作安全提供了保證。

疲勞剝落模型和疲勞剝落擴展模型可計算剝落發(fā)生時間以及失效時間,主要用于研究循環(huán)應(yīng)力作用下的非金屬夾雜物、組織不均勻、碳化物等金屬內(nèi)部發(fā)生的裂紋,在接觸表面下發(fā)生并不斷擴展就導(dǎo)致了剝離。王艷等[14]在確定材料在受到反復(fù)沖擊條件下的疲勞磨損表現(xiàn)為疲勞剝離后,采用膜層來提高材料表面硬度,增加耐磨性。趙新澤等[15]研究了鋼絲鋁絞線在干燥和酸性兩種典型的氣候條件下的微動磨損行為,在觀察試樣的表面外貌和磨損微粒的外貌后,發(fā)現(xiàn)磨損機理為黏著和疲勞剝落。

基于剛度的損傷規(guī)律模型可將軸承的頻率、加速與運行和失效時間關(guān)聯(lián),計算時也需要確定材料常數(shù)。馮培鋒等[16]研究了復(fù)合材料在疲勞作用下的剩余強度衰退模型和剩余剛度衰退模型以及二者之間的關(guān)系,并用實驗對模型進行了驗證。顧怡[17]研究了纖維樹脂復(fù)合材料的剩余強度隨疲勞載荷加載次數(shù)的退化規(guī)律,提出了纖維樹脂復(fù)合材料剩余強度退化模型,并得到了大量實驗結(jié)果的支持。蘇小萍[18]針對復(fù)合材料損傷機理復(fù)雜性和損傷形式多樣性,采用試驗方法來評價損傷困難的問題,研究了采用電阻變化來反映層板結(jié)構(gòu)疲勞損傷的方法。

以上文獻綜述了基于失效物理的故障預(yù)測研究,對于復(fù)雜裝備來說,考慮到裝備的復(fù)雜性和故障隨機性很難建模,因此對于大部分復(fù)雜裝備來說,基于失效物理并不是首選的故障預(yù)測方法。

2 綜合集成故障預(yù)測方法

隨著裝備結(jié)構(gòu)的復(fù)雜化和故障模式的多樣化,使用單一的故障預(yù)測方法的預(yù)測精度很難滿足需求,另外不同的故障預(yù)測方法都有其局限性和適用條件,因此,若能融合多種預(yù)測方法共同進行故障預(yù)測,則能在克服傳統(tǒng)預(yù)測方法不足的同時提高預(yù)測的精確度。趙建忠等[19-20]針對故障預(yù)測時中數(shù)據(jù)采樣間隔不均勻、采樣難度大、數(shù)據(jù)不足的問題,提出了一種基于數(shù)據(jù)融合和改進的灰色模型結(jié)合的故障預(yù)測方法,并且用實例證明了方法的有效性。高甜容等[21]針對現(xiàn)有的數(shù)控機床主軸傳動系統(tǒng)故障預(yù)測的不足,以及故障的多發(fā)性,提出了一種基于故障先兆判定模型和動態(tài)置信度匹配的主軸潤滑故障預(yù)測方法,該方法能夠?qū)?shù)控機床的歷史故障數(shù)據(jù)用于定義故障先兆狀態(tài)序列,并利用小波分析和概率神經(jīng)網(wǎng)絡(luò)構(gòu)建預(yù)測模型,實驗證明,該方法能夠?qū)?shù)控機床主軸傳動系統(tǒng)的故障進行預(yù)測。曾慶虎等[22]為了消除數(shù)據(jù)冗余,獲得高維的故障數(shù)據(jù),對信號進行小波相關(guān)濾波處理,并構(gòu)建了小波相關(guān)特征尺度熵向量,結(jié)合隱半馬爾可夫模型(HSMM)對設(shè)備的退化狀態(tài)進行識別并預(yù)測故障。費勝巍等[23]為了對變壓器進行故障預(yù)測,提出了一種結(jié)合粗糙集和灰色模型的故障預(yù)測方法,這種方法可以結(jié)合預(yù)測對象的狀態(tài)特征和診斷規(guī)則,盡早發(fā)現(xiàn)潛在的故障,及時地采取應(yīng)對措施,避免造成經(jīng)濟損失。文獻[24]提出了一種基于神經(jīng)網(wǎng)絡(luò)和威布爾分布的剩余壽命預(yù)測方法,將軸承的振動信號用于預(yù)測并將軸承的狀態(tài)分為健康狀態(tài)和6個退化狀態(tài)來表征軸承的健康程度,實驗結(jié)果表明此方法可以準確地對剩余壽命進行預(yù)測,并且可以應(yīng)用于除軸承外的其他機械設(shè)備中。文獻[25]針對離散事件后勤系統(tǒng)仿真模型對剩余壽命預(yù)測不可靠的問題,提出了兩步剩余壽命預(yù)測評估離散事件后勤系統(tǒng)仿真模型的規(guī)則和操作方法,并且通過實驗證明兩步剩余壽命預(yù)測方法比傳統(tǒng)的單步預(yù)測方法更加有效。文獻[26]描述了一種結(jié)合粒子群算法與支持向量機的飛機發(fā)動機剩余壽命預(yù)測方法,這種基于粒子群算法與支持向量機的混合模型的優(yōu)點在于其不需要發(fā)動機的歷史狀態(tài)數(shù)據(jù),并且用實驗證明了這種方法的實用性。文獻[27]運用關(guān)聯(lián)向量機對鋰電池進行剩余壽命預(yù)測,但是預(yù)測過程中樣本量少和精確度低的問題限制了關(guān)聯(lián)向量機在鋰電池故障預(yù)測中的應(yīng)用,要提高預(yù)測的精度,就要提高動態(tài)訓(xùn)練和在線學(xué)習能力;并提出了通過增加靈活性和增值優(yōu)化關(guān)聯(lián)向量機算法來對鋰電池進行更加準確的剩余壽命預(yù)測。

3 裝備故障預(yù)測中的不確定性

裝備故障預(yù)測中的不確定性不僅會直接導(dǎo)致預(yù)測準確度下降,還會間接地影響到維修策略的制定。造成故障預(yù)測不確定的因素很多,主要的因素有獲取的退化狀態(tài)參數(shù)不準確、考慮導(dǎo)致故障的因素不全面、模型參數(shù)的不確定性、模型結(jié)構(gòu)的不確定性等。孫強等[28]對近年來的基于不確定性的故障預(yù)測方法進行了綜述,并且通過不確定屬性的不同將故障預(yù)測分為基于隨機性、模糊性、灰性及混合不確定性4類,綜述了每一類方法的研究現(xiàn)狀和不足。

模型參數(shù)和性能參數(shù)的不確定是最常見的引起不確定性的因素,許喆平等[29]針對由于飛機性能參數(shù)包含噪聲和非線性而造成故障預(yù)測不確定性的問題,采用增加線性約束的方式來解決噪聲對故障預(yù)測精度的影響,仿真結(jié)果證明,這種方法能夠降低噪聲,提高故障預(yù)測的準確性,從而增加飛機的安全性。鄧宏琛等[30]將動態(tài)不確定因果圖理論模型應(yīng)用于核電站的故障預(yù)測中,該方法能夠得到關(guān)注的假設(shè)事件及狀態(tài)概率表達,并實現(xiàn)了推理機軟件和算法優(yōu)化。張磊等[31]提出了一種基于關(guān)聯(lián)向量機回歸的故障預(yù)測算法,預(yù)測過程采用多步時間序列預(yù)測中的遞推計算的思路,將每一步預(yù)測中的不確定性作為下一步預(yù)測迭代的輸入因素充分考慮,并用仿真實驗證明改進的算法優(yōu)于傳統(tǒng)的故障預(yù)測算法。

除參數(shù)造成的不確定性,環(huán)境因素也可能對故障預(yù)測的結(jié)果帶來不確定性,由于裝備通常需要在多變的復(fù)雜環(huán)境中作業(yè),因此考慮環(huán)境因素對裝備壽命產(chǎn)生的影響是非常必要的。蔡志強等[32]針對不確定環(huán)境下設(shè)備故障的傳播推理問題,提出了一種基于貝葉斯網(wǎng)絡(luò)的故障推理模型,利用網(wǎng)絡(luò)結(jié)構(gòu)與概率分布有效表達裝備中各部件故障狀態(tài)、關(guān)聯(lián)關(guān)系及傳播方式。徐宇亮等[33]針對導(dǎo)彈電子設(shè)備故障預(yù)測問題,提出了一種基于綜合環(huán)境加速退化試驗和粒子濾波的故障預(yù)測的新方法,并且考慮到環(huán)境應(yīng)力帶來的不確定因素對故障預(yù)測帶來的影響,定義了電子設(shè)備退化度的概念。文獻[34]將貝葉斯框架運用于復(fù)雜多變的隨機環(huán)境中飛機發(fā)動機壽命預(yù)測,提出用4種疲勞模型來驗證貝葉斯框架,通過更新數(shù)據(jù)參數(shù),貝葉斯框架能夠提供有用的性能信息和不確定約束,結(jié)果顯示預(yù)測的疲勞壽命分布符合實驗結(jié)果。

另外,同時考慮隨機和動態(tài)因素產(chǎn)生的不確定性,以及主觀不確定性和客觀不確定性,采用綜合量化的方法來盡量消除不確定性對故障預(yù)測產(chǎn)生的影響也是常見的一類方法。文獻[35]研究了發(fā)動機故障和可靠性預(yù)測中的隨機和動態(tài)不確定性,采用最小二乘法支持向量機和非線性迭代濾波來更新可靠性數(shù)據(jù)的準確性,并且用兩個例子來比較先前文獻中的神經(jīng)網(wǎng)絡(luò)和支持向量機模型與此方法,實驗結(jié)果證明提出的新方法在可靠性預(yù)測有更高的準確性。文獻[36]研究了在不確定條件下飛機渦輪的疲勞壽命預(yù)測,預(yù)測包含了出現(xiàn)在框架結(jié)構(gòu)完整性評估中的不確定因素,采用綜合不確定量化方法來量化多種類型的不確定因素,通過實驗結(jié)果證明預(yù)測模型具有很高的預(yù)測精度。

4 結(jié)束語

通過總結(jié)分析以上文獻可知,在故障預(yù)測領(lǐng)域,很多機構(gòu)已經(jīng)開展了大量研究并取得了豐碩的研究成果,并且很多研究成果已經(jīng)應(yīng)用到實際。故障預(yù)測技術(shù)在解決復(fù)雜裝備的維修問題中起到了越來越重要的作用,遺憾的是故障預(yù)測的準確度和預(yù)測方法的適用性還有待提高,急需發(fā)展新的裝備故障預(yù)測技術(shù)或綜合多種有效的預(yù)測手段來一起解決這個難題。在剩余壽命預(yù)測中,使用預(yù)測壽命區(qū)間代替預(yù)測壽命值能更加符合實際情況,給出預(yù)測壽命值的置信區(qū)間可能成為未來壽命預(yù)測的發(fā)展方向。

筆者綜述了在視情維修背景下的故障預(yù)測的相關(guān)文獻,重點介紹了故障預(yù)測的主要方法及研究現(xiàn)狀。并對所綜述的文獻進行了分析總結(jié),對故障預(yù)測與維修決策未來發(fā)展方向進行了展望。高新技術(shù)在裝備上的廣泛應(yīng)用在帶動裝備發(fā)展的同時也增加了維修保障的難度,裝備“買得起,用不起”的難題也日益突出,亟待提出更加有效的故障預(yù)測技術(shù),盡量發(fā)揮裝備的最大效能。

猜你喜歡
故障方法模型
一半模型
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
故障一點通
3D打印中的模型分割與打包
奔馳R320車ABS、ESP故障燈異常點亮
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
捕魚
故障一點通
主站蜘蛛池模板: 18禁高潮出水呻吟娇喘蜜芽| 999在线免费视频| 国产无码高清视频不卡| 国产本道久久一区二区三区| 日韩少妇激情一区二区| 亚洲精品777| 欧美亚洲欧美区| 55夜色66夜色国产精品视频| 色网站在线视频| 亚洲AV一二三区无码AV蜜桃| 国产人人乐人人爱| 亚洲AV免费一区二区三区| 97se亚洲综合在线天天| 日韩视频福利| 久久一本日韩精品中文字幕屁孩| 99精品伊人久久久大香线蕉 | 国产一区二区三区在线观看免费| 一本大道香蕉中文日本不卡高清二区| 又爽又大又光又色的午夜视频| 亚洲色图综合在线| 狠狠躁天天躁夜夜躁婷婷| 国产免费a级片| 日本在线欧美在线| 四虎永久在线精品国产免费| 亚洲看片网| 极品私人尤物在线精品首页 | 99在线视频网站| 久久精品亚洲专区| 色综合网址| 国产一级在线观看www色| 四虎永久在线| 精品一区二区三区水蜜桃| 欧美视频二区| 嫩草国产在线| 国产国语一级毛片| 亚洲国产成人自拍| 亚洲精品日产AⅤ| 九九久久精品国产av片囯产区| 青青操视频在线| 亚洲第一页在线观看| 免费全部高H视频无码无遮掩| 91 九色视频丝袜| 在线欧美日韩| 久久亚洲日本不卡一区二区| 亚洲精品无码高潮喷水A| 日韩在线观看网站| 国产精品视频久| 免费日韩在线视频| 国产美女无遮挡免费视频网站 | 亚洲综合中文字幕国产精品欧美 | 国产精品lululu在线观看| 天天综合网站| 日韩123欧美字幕| 视频国产精品丝袜第一页| 婷婷成人综合| 国产精品19p| 久久伊伊香蕉综合精品| 日本国产精品一区久久久| 国产交换配偶在线视频| 亚洲午夜福利精品无码不卡| 亚洲中文字幕手机在线第一页| 在线观看网站国产| 婷婷丁香在线观看| 久无码久无码av无码| 精品国产香蕉在线播出| 无遮挡一级毛片呦女视频| 国产裸舞福利在线视频合集| 99久久精品国产综合婷婷| 婷婷亚洲综合五月天在线| 国产又粗又猛又爽视频| 亚洲一区无码在线| 国产中文在线亚洲精品官网| 制服丝袜无码每日更新| 日韩国产综合精选| 国产成人a在线观看视频| 亚洲永久免费网站| 国产精品女主播| 无码一区二区波多野结衣播放搜索| 亚洲福利网址| 久夜色精品国产噜噜| 亚洲欧美色中文字幕| 1级黄色毛片|