日本在线www-日本在线播放一区-日本在线不卡免费视频一区-日本在线不卡视频-成人影院久久久久久影院-成人影院一区二区三区

最新廣告
關(guān)注中國自動(dòng)化產(chǎn)業(yè)發(fā)展的先行者!
工業(yè)智能邊緣計(jì)算2025年會(huì)
CAIAC 2025
2025工業(yè)安全大會(huì)
OICT公益講堂
當(dāng)前位置:首頁 >> 案例 >> 案例首頁

案例頻道

基于PCA和RBF網(wǎng)絡(luò)的減速箱運(yùn)行狀態(tài)診斷技術(shù)
  • 企業(yè):控制網(wǎng)     領(lǐng)域:工業(yè)以太網(wǎng)    
  • 點(diǎn)擊數(shù):2987     發(fā)布時(shí)間:2008-01-12 03:49:32
  • 分享到:



    楊帆(1973—)
男,陜西禮泉人,講師,碩士,主要研究領(lǐng)域?yàn)橹悄苄畔⑻幚砑夹g(shù)、計(jì)算機(jī)監(jiān)控技術(shù)。

基金項(xiàng)目:陜西科技大學(xué)自然科學(xué)基金資助(ZX05-37);陜西省工業(yè)攻關(guān)項(xiàng)目(2006K05-G18)

摘要:針對減速箱運(yùn)行狀態(tài)和特征參數(shù)之間存在的復(fù)雜非線性關(guān)系,提出了基于主成分分析的RBF神經(jīng)網(wǎng)絡(luò)減速箱運(yùn)行狀態(tài)診斷方法。該方法用主成分分析方法將高維相關(guān)特征參數(shù)轉(zhuǎn)化為低維相互獨(dú)立的特征參數(shù),在此基礎(chǔ)上建立了RBF網(wǎng)絡(luò)分類器,并用該網(wǎng)絡(luò)對某汽輪機(jī)減速箱的運(yùn)行狀態(tài)進(jìn)行識(shí)別。理論分析和實(shí)驗(yàn)結(jié)果表明,基于PCA和RBF網(wǎng)絡(luò)方法的減速箱運(yùn)行狀態(tài)診斷技術(shù)具有模型簡單,檢測速度快等優(yōu)點(diǎn),可以在實(shí)際應(yīng)用中發(fā)揮有效作用。

關(guān)鍵詞:主成分分析;RBF網(wǎng)絡(luò);特征提取;狀態(tài)診斷

Abstract: As to the complicated nonlinear relation existing between running status of gear reducer and characteristic parameters, PCA-based RBF neural network reducer running status diagnostics is put forward. High-dimensional correlated characteristic parameters are transformed into low-dimensional independent characteristic parameters, a RBF neural network classifier is built hereto and the network is used to identify a turbine reducer running status. Theory analysis and Experiment result indicate that PCA and RBF neural network based reducer running status diagnostics has advantages of simple model and quick detecting speed, and can play effective role in practical application.

Key Words: PCA;RBF network;Characteristics extraction;Status Diagnosis

    減速箱是各類機(jī)械傳動(dòng)部分的關(guān)鍵部件。根據(jù)研究和統(tǒng)計(jì)分析表明,減速箱內(nèi)各零件失效比最大,約占機(jī)械傳動(dòng)系統(tǒng)故障的 60%。因此進(jìn)行狀態(tài)監(jiān)測和故障診斷,識(shí)別設(shè)備的工作狀態(tài),及時(shí)發(fā)現(xiàn)異常情況,對故障進(jìn)行早期診斷和預(yù)報(bào),從而有針對性、有計(jì)劃性地采取檢修措施,確保設(shè)備安全運(yùn)行就顯得非常重要。

    在機(jī)械故障診斷中一般描述系統(tǒng)采用的參數(shù)越多,被認(rèn)為對系統(tǒng)的認(rèn)識(shí)越深刻。但是如果利用過多的系統(tǒng)參數(shù)作為診斷識(shí)別系統(tǒng)的數(shù)據(jù),一方面由于參數(shù)之間不可避免的相關(guān)性以及量測過程中噪聲的引入,就會(huì)使得常用的診斷算法性能下降;另外一方面,引入的參數(shù)越多,就會(huì)占用大量的機(jī)器處理時(shí)間和存貯空間,影響診斷的速度。所以在診斷過程中提取核心特征參數(shù)、抑制噪聲數(shù)據(jù)是簡化計(jì)算過程,提高診斷識(shí)別率的重要步驟。本文采用了主成分分析( PCA) 方法進(jìn)行特征再提取,通過創(chuàng)建一個(gè)替換的、較小的變量集來“組合”原始特征的精華,原始數(shù)據(jù)可以投影到該較小的集合中。在PCA特征提取的基礎(chǔ)上,用RBF神經(jīng)網(wǎng)絡(luò)分類器進(jìn)行了某汽輪機(jī)減速箱的運(yùn)行狀態(tài)診斷,取得了較好的診斷結(jié)果。

1 PCA方法基本原理

    成分分析是用來在數(shù)據(jù)中尋找恰當(dāng)?shù)奶卣鞅硎镜姆椒āV鞒煞址治觯╬rinciple component analysis)PCA一般用來對輸入樣本集進(jìn)行預(yù)處理,可把多個(gè)存在復(fù)雜非線性關(guān)系的特征參數(shù)轉(zhuǎn)化為較少的彼此不相關(guān)的綜合變量,接著以這些較少的向量形成新的網(wǎng)絡(luò)輸入變量,從而降低神經(jīng)網(wǎng)絡(luò)輸入的變量數(shù),優(yōu)化網(wǎng)絡(luò)結(jié)構(gòu)。PCA的對象是過程變量的樣本數(shù)據(jù)矩陣。數(shù)據(jù)矩陣的行表示采樣值或觀察值,列表示變量。PCA產(chǎn)生一個(gè)壓縮的統(tǒng)計(jì)模型——主元模型,模型給出了變量的線性組合,描述了數(shù)據(jù)變化的主要趨勢。主元模型使原標(biāo)準(zhǔn)差的平方重新分布,大多數(shù)標(biāo)準(zhǔn)差平方會(huì)分布在第一主元上,其次分布在第二主元上,依此類推。按某種準(zhǔn)則將最后幾個(gè)主元視為分解殘差予以忽略,則有可能利用最少主元來說明最多的信息。主元模型舍棄了部分殘差而保留體現(xiàn)數(shù)據(jù)變異的主要方向,從而達(dá)到抽取系統(tǒng)信息,清除系統(tǒng)干擾的目的。

    實(shí)際應(yīng)用中,合理確定主元個(gè)數(shù)進(jìn)而建立主元模型非常重要。通常采用方差累積貢獻(xiàn)率百分比(CPV)原則,選擇百分比大于85%的主元個(gè)數(shù)。在基于PCA的過程狀態(tài)、故障的診斷中,一般用來描述正常運(yùn)行過程的主元不超過3個(gè)(Kresta等,1991)。主成分分析的主要計(jì)算步驟如下:

    (1) 原始樣本標(biāo)準(zhǔn)化

    為了消除量綱和數(shù)量級不同的影響,采用均值標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化方法處理原始樣本數(shù)據(jù)。

    (2) 建立標(biāo)準(zhǔn)化變量的協(xié)方差矩陣,求解矩陣的特征值和特征向量

    利用標(biāo)準(zhǔn)化值計(jì)算變量之間的相關(guān)系數(shù),有k個(gè)特征參數(shù)可建立k階相關(guān)矩陣。由此矩陣可獲得由大到小排列的特征值A(chǔ)i(i=1,2,L,k),k個(gè)特征值對應(yīng)k個(gè)特征向量,每一特征向量包含k個(gè)分量。

    (3) 根據(jù)要求的累計(jì)貢獻(xiàn)率, 選取主成分

    計(jì)算第i個(gè)主成分對總方差的貢獻(xiàn)率,即方差貢獻(xiàn)率:。通常取累計(jì)方差貢獻(xiàn)率大于85% 所需的主成分,能夠代表k 個(gè)原始變量所能提供的絕大部分信息。

    (4) 建立主成分方程、計(jì)算各主成分值

    各主成分值方程為:,其中aj為對應(yīng)于j 的特征向量的分量, xj為各變量的標(biāo)準(zhǔn)化數(shù)值。計(jì)算出所需要的各主成分值, 形成新的訓(xùn)練樣本集和測試樣本集。

2 基于PCA的RBF網(wǎng)絡(luò)診斷方法

    2.1 RBF神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

    神經(jīng)網(wǎng)絡(luò)對復(fù)雜問題具有自適應(yīng)和自學(xué)習(xí)能力,為解決復(fù)雜系統(tǒng)的信息處理和控制等問題提供了新的思想和方法。作為一種前饋型神經(jīng)網(wǎng)絡(luò),RBF網(wǎng)絡(luò)避免了BP神經(jīng)網(wǎng)絡(luò)的冗長繁瑣計(jì)算,學(xué)習(xí)速度較通常的BP方法快的多,具有良好的泛化能力,能以任意精度逼近非線性函數(shù)。RBF神經(jīng)網(wǎng)絡(luò)一般組成結(jié)構(gòu)如圖1所示:


圖1  RBF神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖

    RBF網(wǎng)絡(luò)結(jié)構(gòu)中,輸入層由信號(hào)源節(jié)點(diǎn)組成,傳遞信號(hào)到隱層;輸出層一般是簡單的線性函數(shù),對輸入模式作出響應(yīng)。隱層節(jié)點(diǎn)的變換函數(shù)(基函數(shù))是對中心點(diǎn)徑向?qū)ΨQ且衰減的非負(fù)非線性函數(shù),對輸入信號(hào)將在局部產(chǎn)生響應(yīng),也就是說,當(dāng)輸入信號(hào)靠近基函數(shù)的中央范圍時(shí),隱層節(jié)點(diǎn)將產(chǎn)生較大的輸出,由此看出這種網(wǎng)絡(luò)具有局部逼近能力,所以徑向基函數(shù)網(wǎng)絡(luò)也成為局部感知場網(wǎng)絡(luò)。作為基函數(shù)的形式,最常用的一般是高斯函數(shù):

   

    其中x是n維輸入向量;Ci是第i個(gè)基函數(shù)的中心,與x具有相同維數(shù)的向量,i是i個(gè)感知的變量,決定該基函數(shù)圍繞中心點(diǎn)的寬度;m是感知單元的個(gè)數(shù);是向X-Ci量的范數(shù),表示兩個(gè)向量之間的距離,Ri(x)在ci處有唯一的最大值,隨著的增大,Ri(x)迅速衰減到0。對于給定的輸入,只有一小部分靠近x的中心被激活。

    從RBF網(wǎng)絡(luò)結(jié)構(gòu)可以看出,輸入層實(shí)現(xiàn)從的非線性映射,輸出層實(shí)現(xiàn)的線性映射,即:,其中是輸出節(jié)點(diǎn)的個(gè)數(shù)。

    RBF 神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程主要分為兩個(gè)階段。首先,根據(jù)輸入樣本求出各隱層節(jié)點(diǎn)高斯核函數(shù)的中心ci和每個(gè)中心的半徑Ri。隱層參數(shù)確定以后, 還需要求出隱層和輸出層之間的權(quán)值wik。

    2.2 基于PCA的RBF網(wǎng)絡(luò)診斷模型

    基于PCA和RBF神經(jīng)網(wǎng)絡(luò)的減速箱運(yùn)行狀態(tài)診斷系統(tǒng)實(shí)現(xiàn)框圖如圖2所示。



圖2  基于PCA和RBF神經(jīng)網(wǎng)絡(luò)的減速箱運(yùn)行狀態(tài)診斷系統(tǒng)實(shí)現(xiàn)框圖

    如前所述,基于PCA的RBF診斷就是利用PCA的方法進(jìn)行原始數(shù)據(jù)的預(yù)處理,從而獲取系統(tǒng)描述主成分,然后利用主成分構(gòu)造出診斷算法需要的訓(xùn)練數(shù)據(jù)集,再基于RBF神經(jīng)網(wǎng)絡(luò)進(jìn)行學(xué)習(xí)和訓(xùn)練。對于采集到的未知狀態(tài)的運(yùn)行數(shù)據(jù),首先將其變換到主成分空間,然后利用訓(xùn)練好的RBF網(wǎng)絡(luò)進(jìn)行運(yùn)行狀態(tài)的診斷和識(shí)別。其具體處理步驟如下:

    (1)獲取描述某汽輪機(jī)減速箱運(yùn)行特征參數(shù)以及運(yùn)行狀態(tài)分類決策結(jié)果,形成原始數(shù)據(jù)集。

    (2)基于PCA方法抽取訓(xùn)練數(shù)據(jù)集的主成分;

    (3)構(gòu)造基于主成分的訓(xùn)練數(shù)據(jù)集;

    (4)基于主成分訓(xùn)練數(shù)據(jù)集構(gòu)造并訓(xùn)練RBF神經(jīng)網(wǎng)絡(luò)模式分類器;

    (5)將新的觀測向量(與原始數(shù)據(jù)維數(shù)相同)變換到主成分空間,得到主成分空間的觀測對象變量。

    在主成分的訓(xùn)練數(shù)據(jù)集上基于訓(xùn)練好的RBF神經(jīng)網(wǎng)絡(luò)進(jìn)行推理分類。

3 實(shí)驗(yàn)過程以及結(jié)果分析

    為了驗(yàn)證本文提出方法的正確性和有效性,選取一組實(shí)際采集的某汽輪機(jī)減速箱運(yùn)行實(shí)驗(yàn)數(shù)據(jù)(表1) 進(jìn)行運(yùn)算驗(yàn)證。

    表1  汽輪機(jī)減速箱運(yùn)行狀態(tài)特征數(shù)據(jù):

    表1中共有100組實(shí)測減速箱運(yùn)行數(shù)據(jù),a1~a8為采集到的8個(gè)特征參數(shù),表中狀態(tài)值1表示正常,2表示故障。將樣本數(shù)據(jù)分為2組,選前80個(gè)數(shù)據(jù)進(jìn)行訓(xùn)練網(wǎng)絡(luò),后20個(gè)數(shù)據(jù)進(jìn)行檢驗(yàn)。為了消除由于數(shù)據(jù)量綱對于最后判決帶來的誤差,首先進(jìn)行歸一化處理,然后進(jìn)行主成分分析,計(jì)算得到主成分的方差累計(jì)貢獻(xiàn)率為:76.2117、86.0981、93.3289、98.1578、   99.3166、99.8645、99.9573、100.0000,按照方差累計(jì)貢獻(xiàn)率>90%獲取主成分,最終生成3個(gè)主成分x1,x2,x3,其計(jì)算方法如下所示。
u1=(-0.3986  -0.3657  -0.3570  -0.3944  -0.3307  -0.3760  -0.2335  -0.3452)’;

    u2=(0.0189  -0.3638  -0.3765  -0.07312  -0.0469  0.0691  0.8024  0.2633 )’;

    u3=(-0.1959  0.2471  0.3770  0.1793  0.0052   -0.4213  0.5197  -0.5279)’;

    使U=(u1 u2 u3)’,與標(biāo)準(zhǔn)化處理后得到的數(shù)據(jù)集X按照式計(jì)算,即可得到原始數(shù)據(jù)的主要成分表示,如表2所示。

    表2  主成分訓(xùn)練數(shù)據(jù)集:

    利用主成分訓(xùn)練數(shù)據(jù)集的前80個(gè)數(shù)據(jù)進(jìn)行RBF網(wǎng)絡(luò)的構(gòu)造和訓(xùn)練,后20個(gè)作為測試數(shù)據(jù),試驗(yàn)結(jié)果表明:基于PCA和RBF神經(jīng)網(wǎng)絡(luò)的減速箱運(yùn)行狀態(tài)診斷方法的診斷結(jié)果與實(shí)際完全一致。可見,基于PCA方法將觀測向量投影到低維空間,構(gòu)造出相互獨(dú)立的主成分觀測向量后,可以簡化RBF神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)和計(jì)算,對于提高診斷速度的目標(biāo)作用十分明顯。

4 結(jié)束語

    在減速器狀態(tài)監(jiān)測和故障診斷中,將主成分分析與RBF神經(jīng)網(wǎng)絡(luò)技術(shù)相結(jié)合,通過對原始特征參數(shù)的降維處理,得到互不相關(guān)低維主成分空間特征參數(shù),以此為RBF神經(jīng)網(wǎng)絡(luò)的輸入向量,進(jìn)行網(wǎng)絡(luò)的構(gòu)造和訓(xùn)練,最后利用訓(xùn)練好的網(wǎng)絡(luò)進(jìn)行減速箱的運(yùn)行狀態(tài)診斷。相比傳統(tǒng)的單純神經(jīng)網(wǎng)絡(luò)狀態(tài)診斷方法,本文提出的診斷方法由于PCA技術(shù)的引入,減少了輸入向量的維度,消除了相關(guān)信息和噪聲數(shù)據(jù)的影響,簡化了神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu),提高了分析運(yùn)算的速度,可以在減速箱運(yùn)行狀態(tài)的實(shí)時(shí)診斷中發(fā)揮有效作用。

其它作者:

    張玉杰(1967-),男,陜西武功人,副教授,碩士,主要從事嵌入式智能系統(tǒng)研究;

    張彩麗(1973-),女,陜西合陽人,碩士,主要從事機(jī)電系統(tǒng)故障診斷技術(shù)研究。

參考文獻(xiàn)

    [1]范金城. 數(shù)據(jù)分析[M]. 北京:科學(xué)出版社,2002 141~153.

    [2]梁循. 數(shù)據(jù)挖掘算法與應(yīng)用[M]. 北京:北京大學(xué)出版社,2006 184~187.

    [3]任若恩. 多元統(tǒng)計(jì)數(shù)據(jù)分析[M]. 北京:國防工業(yè)出版社,1997 92~109.

    [4]熊和金. 智能信息處理[M]. 北京:國防工業(yè)出版社,2006 123~135.

    [5]殷勤業(yè). 模式識(shí)別與神經(jīng)網(wǎng)絡(luò)[M]. 北京:機(jī)械工業(yè)出版社,1992 98~104.

    [6]曹龍漢. 柴油機(jī)智能化故障診斷技術(shù)[M]. 北京:國防工業(yè)出版社,2005 131~135.

    [7] Luo, Rongfu; Misra, Manish; Himmelblau, David M. Sensor fault detection via multiscale analysis and dynamic PCA Industrial and Engineering Chemistry Research, v 38, n 4, Apr, 1999, p 1489-1495.

    [8] Abdel-Qader, Ikhlas, PCA-Based algorithm for unsupervised bridge crack detection,Advances in Engineering Software, v 37, n 12, December, 2006, p 771-778.

熱點(diǎn)新聞

推薦產(chǎn)品

x
  • 在線反饋
1.我有以下需求:



2.詳細(xì)的需求:
姓名:
單位:
電話:
郵件:
主站蜘蛛池模板: 精品国产无码一区二区-精品国产无套在线观看-精品国产无限资源免费观看-精品国产午夜福利在线观看-精品国产午夜久久久久九九-精品国产午夜肉伦伦影院 | 美女三级毛片-美女牲交视频一级毛片无遮挡-美女视频大全视频a免费九-美女视频黄a视频全免费网站色窝-美女视频黄的全是免费-美女视频秀色福利视频 | 视频区国产亚洲.欧美-视频黄色在线-视频国产激情-十九岁日本电影免费粤语高清-久久99九九精品免费-久久99九九国产免费看小说 | 亚洲最大色网-亚洲最大色图-亚洲最大情网站在线观看-亚洲最大免费视频网-九一自拍-九一制片厂制作果冻传媒网站 | 婷婷综合激情五月中文字幕,美国性网站,国产精品亚洲欧美日韩久久,每日更新在线观看,久久免费视屏,www.成年人视频 | 天天噜噜色-天天看天天射天天视频-天天看天天射天天碰-天天看天天碰-国产成人高清-国产成人爱情动作片在线观看 | 激情综合欧美-激情综合激情五月-激情综合丁香-激情综-四虎影视精品-四虎影视国产在线观看精品 | 免费的成人a视频在线观看-免费岛国小视频在线观看-免费岛国-免费大香伊蕉在人线国产-免费大片黄在线观看日本-免费大片黄手机在线观看 | 亚洲综合在线观看视频-亚洲综合在线观看一区-亚洲综合在线观看一区www-亚洲综合在线视频-亚洲综合在线网-亚洲综合在线一区 | 精品无码免费黄色网站-精品无码免费一区二区三区-精品无码免费在线播放-精品无码免费专区-精品无码免费专区午夜-精品无码欧美黑人又粗又 | 二区乱码综合无码一区二区三-二区女人观看chinese中国真实乱-二区欧美三-二区欧美无遮挡中文字幕人成人-二区啪视频-二区日本成人动漫电影 | 一区二区国产精品-一区二区高清在线观看-一区二区高清在线-一区二区成人国产精品-九九视频精品全部免费播放-九九视频国产 | 久久道-久久大香伊蕉在人线国产昨爱-久久大香香蕉国产免费网站-久久大-天天做天天爽天天谢-天天做天天射 | 国产成人精品一区二区视频免费-国产成人精品一区二区视频-国产成人精品一区二区三在线观看-国产成人精品一区二区三在-国产成人精品一区二区三区小说-国产成人精品一区二区三区无码 | 国产美女小视频-国产美女网站视频-国产美女网站-国产美女特级嫩嫩嫩bbb-天天干视频在线观看-天天干视频在线 | 一个人看的在线www高清视频-一个人看的小说在线阅读-一个人看的手机视频www-一个人看的视频在线观看免费播放动漫-久久99精品久久久久久秒播放器-久久99精品久久久久久秒播 | 视频免费观看网址-视频免费观看视频-视频免费观看入口大全-视频免费1区二区三区-国产成人尤物精品一区-国产成人悠悠影院 | 欧美成人精品动漫在线专区-欧美成人精品区综合A片-欧美成人精品三区综合A片-欧美成人精品一区二区综合A片-欧美成人另类-欧美成人猛片AAAAAAA | 99欧美在线-99碰碰-99青草-99青草青草久热精品视频-99青青-99青青草 | 欧美一区不卡二区不卡三区,欧美另类日韩,日韩中文字幕免费版,亚洲一区二区免费看,欧美天天,亚洲欧美另类专区 | 久热久热精品在线观看-久青草影院-久青草国产在线视频亚瑟影视-久青草国产在线视频-正在播放黑人杂交派对卧槽-正在播放国产尾随丝袜美女 | 欧美视频综合-欧美视频中文字幕-欧美视频在线观在线看-欧美视频在线观看网站-直接观看黄网站免费视频-正在播放久久 | 久久久毛片-久久久麻豆-久久久鲁-久久久噜噜噜久久久-午夜你懂得-午夜男人网 | 九九激情网,日韩色综合,成人小视频网站,国产永久在线观看,污黄视频在线观看,看国产一级片 | 青草国产-青草草在线视频-青草草在线观看免费视频-青草草在线-国产高清自偷自在线观看-国产高清自拍视频 | 日韩美在线-日韩美一区二区三区-日韩美女在线视频-日韩美女视频在线观看-亚洲无线一二三四区手机-亚洲无专砖码直接进入 | 亚洲免费在线观看-做羞羞的事情的免费视频-最终痴汉电车在线观看-最新综艺-最新自拍偷拍-最新在线精品国自拍视频 | 波多野结衣的av一区二区三区-波多野结衣的电影-波多野结衣的中文-波多野结衣第二页视频-波多野结衣电影网-波多野结衣电影一区二区 | 视频免费观看网址-视频免费观看视频-视频免费观看入口大全-视频免费1区二区三区-国产成人尤物精品一区-国产成人悠悠影院 | 丝袜情趣在线资源二区-丝袜啪啪-丝袜秘书ol连裤袜在线播放-丝袜美女被出水一区-国产国产人在线成免费视频69-国产国产人免费视频成69堂 | 国产成a人亚洲精v品无码性色-国产成a人亚洲精v品无码樱花-国产成a人亚洲精v品在线观看-国产成a人亚洲精品无码樱花-国产成版人视频网-国产成版人视频网站免费下 | 韩国三级一区-韩国三级香港三级日本三级la-韩国三级香港三级日本三级-韩国三级视频网站-日韩欧美一及在线播放-日韩欧美一二三区 久久久久久久久国产-久久久久久久久97-久久久久久久国产视频-久久久久久久国产精品影院-午夜精-午夜寂寞院 | 精品久久久久久蜜臂a∨-精品久久久久久蜜臀-精品久久久久久免费看-精品久久久久久免费影院-精品久久久久久清纯-精品久久久久久日韩字幕无 | 欧美成人国产一区二区-欧美成人黄色-欧美成人黄色片-欧美成人家庭影院-欧美成人精精品一区二区三区-欧美成人精品a8198v无码 | 夜夜操狠狠操,99在线观看精品免费99,亲爱的妈妈2,色婷婷亚洲十月十月色天,精品福利一区二区在线观看,久久国产精品99精品国产 | 爆乳无码一区二区三区-爆乳熟妇一区二区三区霸乳-爆乳熟妇一区-爆乳少妇在办公室在线观看-爆乳护士一区二区三区在线播放-白丝一区二区三区 | 无人 视频 在线观看免费高清-无吗视频-无码中文字幕日韩专区-无码中文字幕乱在线观看-国产区1-国产情侣一区二区 | 制服师生一区二区三区在线-制服师生av在线-制服 丝袜 亚洲 中文 综合-直接看毛片-欧美视频在线观看视频-欧美视频在线观看免费最新 | 欧美中出,国产欧美又粗又猛又爽老,日本啊v,欧美日韩欧美,国产va在线观看,国产一区二区三区在线视频 | 欧美日韩亚洲一区二区三区在线观看-欧美日韩亚洲第一区在线-欧美日韩亚洲成人-欧美日韩午夜群交多人轮换-bbwvideos欧美老妇-bbwvideoa欧美老妇 | 国产三级高清午夜羞羞视频-国产三级高清在线观看-国产三级观看久久-国产三级国产av品爱网-国产三级国产精品-国产三级国产精品国产国在线观看 |