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

    基于比例風(fēng)險(xiǎn)模型的地鐵車門可靠性評(píng)估方法研究

    2015-12-28 07:50:28何廣堅(jiān)任金寶邢宗義
    機(jī)械制造與自動(dòng)化 2015年3期

    何廣堅(jiān),任金寶,邢宗義

    (1. 廣州市地下鐵道總公司,廣東 廣州 510030; 2. 南京理工大學(xué),江蘇 南京 210094)

    ?

    基于比例風(fēng)險(xiǎn)模型的地鐵車門可靠性評(píng)估方法研究

    何廣堅(jiān)1,任金寶2,邢宗義2

    (1. 廣州市地下鐵道總公司,廣東 廣州 510030; 2. 南京理工大學(xué),江蘇 南京 210094)

    摘要:車門系統(tǒng)作為地鐵車輛的關(guān)鍵系統(tǒng),對(duì)車輛的安全性、舒適性有著重要的影響,故引入比例風(fēng)險(xiǎn)模型對(duì)車門系統(tǒng)的可靠性進(jìn)行評(píng)估。對(duì)車門系統(tǒng)的歷史檢修數(shù)據(jù)進(jìn)行篩選得到協(xié)變量因子;應(yīng)用主成分分析法對(duì)協(xié)變量進(jìn)行處理得到協(xié)變量矩陣;采用極大似然估計(jì)法和牛頓迭代法求取比例風(fēng)險(xiǎn)模型參數(shù)估計(jì)值;得到車門系統(tǒng)的比例風(fēng)險(xiǎn)模型。通過比例風(fēng)險(xiǎn)模型計(jì)算車門系統(tǒng)可靠度下降為0.95時(shí),車門壽命為170天。比例風(fēng)險(xiǎn)模型所得結(jié)果與Weibull分布模型所得結(jié)果比較,可知比例風(fēng)險(xiǎn)模型所得結(jié)果更符合現(xiàn)場(chǎng)工程師經(jīng)驗(yàn)。

    關(guān)鍵詞:車門可靠性;比例風(fēng)險(xiǎn)模型;似然估計(jì);牛頓迭代

    0引言

    車門系統(tǒng)作為地鐵車輛的關(guān)鍵系統(tǒng),對(duì)車輛的安全性、舒適性、維護(hù)方便性以及整體美觀性,有著重要的影響。根據(jù)國(guó)內(nèi)外統(tǒng)計(jì),在地鐵車輛系統(tǒng)中車門系統(tǒng)的故障占車輛系統(tǒng)總故障的30%以上,對(duì)客車運(yùn)行的安全性構(gòu)成了嚴(yán)重的威脅[1]。因此對(duì)地鐵車門系統(tǒng)進(jìn)行可靠性分析研究,對(duì)保證地鐵車輛運(yùn)行安全具有重要的理論意義和實(shí)用價(jià)值。

    針對(duì)車門系統(tǒng)這一高故障系統(tǒng),國(guó)內(nèi)外學(xué)者做了大量的可靠性分析研究。朱小娟等[2]采用故障樹分析法對(duì)車門系統(tǒng)的可靠性進(jìn)行了分析,得出導(dǎo)致車門故障的所有可能因素和薄弱環(huán)節(jié);蔡國(guó)強(qiáng)等[3]將GO法的操作符特征量計(jì)算方法應(yīng)用于城市軌道交通車輛車門系統(tǒng)可靠性分析,得出了車門系統(tǒng)的可靠度指標(biāo);Thierry Lecomte等[4]對(duì)車門控制系統(tǒng)進(jìn)行可靠性分析,得出控制系統(tǒng)的薄弱環(huán)節(jié),用來對(duì)控制系統(tǒng)提出改進(jìn)建議。以上學(xué)者的研究為地鐵車輛的安全運(yùn)行做出了貢獻(xiàn),但沒有分析車門系統(tǒng)中的日常檢修對(duì)車門系統(tǒng)壽命的影響。

    由于目前尚未檢索到檢修車門系統(tǒng)壽命預(yù)測(cè)的相關(guān)文獻(xiàn),現(xiàn)引入比例風(fēng)險(xiǎn)模型,將檢修數(shù)據(jù)作為協(xié)變量,對(duì)車門系統(tǒng)進(jìn)行可靠性分析研究。首先,將得到的車門系統(tǒng)檢修數(shù)據(jù)進(jìn)行初步篩選,得到協(xié)變量因子;其次,利用主成分分析法對(duì)協(xié)變量因子進(jìn)行預(yù)處理得到協(xié)變量矩陣;然后,利用極大似然估計(jì)和牛頓迭代法得到比例風(fēng)險(xiǎn)模型中的參數(shù)估計(jì)量;最后,得出車門系統(tǒng)的生存函數(shù),預(yù)測(cè)車門系統(tǒng)壽命。

    1比例風(fēng)險(xiǎn)模型

    比例風(fēng)險(xiǎn)模型PHM(proportional hazard model)是基于統(tǒng)計(jì)回歸的壽命預(yù)測(cè)方法中的一種,最早由Cox在1972年提出[5],之后很快成為一種統(tǒng)計(jì)數(shù)據(jù)分析工具,并逐步在生物醫(yī)學(xué)領(lǐng)域得到廣泛應(yīng)用。在可靠性工程領(lǐng)域,比例風(fēng)險(xiǎn)模型主要應(yīng)用于壽命預(yù)測(cè)和維修時(shí)序規(guī)劃。使用比例風(fēng)險(xiǎn)模型進(jìn)行設(shè)備的剩余壽命預(yù)測(cè)時(shí),需要完成樣本數(shù)據(jù)處理、參數(shù)估計(jì)、壽命預(yù)測(cè)3個(gè)基本步驟。首先介紹比例風(fēng)險(xiǎn)模型的基本形式,然后分別在樣本數(shù)據(jù)、參數(shù)估計(jì)、壽命預(yù)測(cè)3個(gè)方面介紹比例風(fēng)險(xiǎn)模型進(jìn)行壽命預(yù)測(cè)的過程。

    1.1基本形式

    一般地帶時(shí)變協(xié)變量的比例風(fēng)險(xiǎn)模型的函數(shù)形式[6]:

    (1)

    (2)

    將式(2)代入式(1),可得:

    (3)

    式(3)所示即為地鐵車門系統(tǒng)所用的比例風(fēng)險(xiǎn)模型基本形式。

    1.2樣本數(shù)據(jù)處理

    利用比例風(fēng)險(xiǎn)模型進(jìn)行壽命預(yù)測(cè)需要從檢修數(shù)據(jù)中得到以下數(shù)據(jù):

    a) 設(shè)備從開始運(yùn)行至失效或截尾的時(shí)間;

    b) 設(shè)備從開始運(yùn)行至失效或截尾過程中的協(xié)變量矢量;

    協(xié)辦量是影響比例風(fēng)險(xiǎn)模型的重要因素,協(xié)變量選取的結(jié)果直接關(guān)系比例風(fēng)險(xiǎn)模型的壽命預(yù)測(cè)準(zhǔn)確性,因此對(duì)協(xié)變量作如下處理。

    1) 協(xié)變量的預(yù)處理

    比例風(fēng)險(xiǎn)模型要求各協(xié)變量間的偏相關(guān)系數(shù)盡量小,主成分分析法是一種有效消除相關(guān)性的方法,所以可用主成分分析法得到車門系統(tǒng)相互獨(dú)立的協(xié)變量參數(shù)。假設(shè)車門系統(tǒng)協(xié)變量為X1,X2,,Xm,則對(duì)協(xié)變量數(shù)據(jù)作標(biāo)準(zhǔn)化處理,得到相關(guān)系數(shù)矩陣[8]:

    (4)

    式中,矩陣元素rij(i,j=1,2,,m)為協(xié)變量Xi于Xj的相關(guān)系數(shù)。

    利用主成分分析法得到協(xié)變量矩陣為:

    (5)

    2) 協(xié)變量綜合

    由于比例風(fēng)險(xiǎn)模型中的協(xié)變量參數(shù)有多個(gè),應(yīng)把原始數(shù)據(jù)中的多個(gè)協(xié)變量進(jìn)行綜合,綜合協(xié)變量函數(shù)如式(6)。

    (6)

    對(duì)原始協(xié)變量進(jìn)行線性回歸擬合[12],即:

    (7)

    可得協(xié)變量的綜合值如式(8)。

    (8)

    1.3模型參數(shù)估計(jì)

    極大似然估計(jì)具有優(yōu)良的統(tǒng)計(jì)性質(zhì)和較好的近似分布,同時(shí)考慮樣本數(shù)據(jù)中含有截尾數(shù)據(jù),故采用極大似然方法來得到模型中各有關(guān)參數(shù)的估計(jì)值[8-12]。

    設(shè)有n個(gè)樣本數(shù)據(jù),則模型的似然函數(shù)為:

    lnL(β,η,γ)=

    式中:q為失效的樣本總數(shù),p為協(xié)變量的總個(gè)數(shù),k=1, …,p。

    采用牛頓-拉夫森(N-R)迭代法求模型參數(shù),假設(shè)X=[,,]=[,,1,2, …,m],對(duì)式(9)分別求各參數(shù)的一階偏導(dǎo)數(shù)、二階偏導(dǎo)數(shù)和混合偏導(dǎo)數(shù)得到矩陣F(X)和矩陣G(X)。

    把矩陣F(X)和矩陣G(X)作為迭代因子,利用牛頓迭代公式可得:

    (10)

    1.4壽命預(yù)測(cè)

    由Weibull分布為基底函數(shù)的比例風(fēng)險(xiǎn)模型的失效率函數(shù)可得生存函數(shù)為:

    (11)

    由生存函數(shù)可得特定可靠度閾值下設(shè)備的無故障運(yùn)行時(shí)間。

    (12)

    由式(11)和式(12)可計(jì)算車門系統(tǒng)運(yùn)行到某一時(shí)間段的可靠度和在某一可靠度下車門無故障運(yùn)行時(shí)間。

    2實(shí)例分析

    選取某地鐵公司2號(hào)線的檢修數(shù)據(jù),假設(shè)年檢或架修后則認(rèn)為車門系統(tǒng)修復(fù)如新。通過分析2號(hào)線2011~2012年的檢修數(shù)據(jù),可知在車輛一次壽命周期內(nèi)正線故障、車輛的維修次數(shù)是影響車輛壽命的主要因素,人為因素導(dǎo)致車門無法正常工作也是出現(xiàn)次數(shù)較多的因素,因此選取正線故障次數(shù)、車輛維修次數(shù)以及由于乘客導(dǎo)致的車門故障次數(shù)作為車門系統(tǒng)的協(xié)變量,如表1。

    表1 車門系統(tǒng)檢修歷史記錄

    續(xù)表1

    對(duì)表1中原始協(xié)變量數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理,得到相關(guān)系數(shù)矩陣

    求出矩陣R的特征值和特征向量,得到協(xié)變量矩陣:

    通過線性回歸模型將多個(gè)協(xié)變量綜合成綜合協(xié)變量影響因子,得到協(xié)變量線性回歸的擬合結(jié)果如表2。

    表2 協(xié)變量參數(shù)的線性擬合結(jié)果

    結(jié)合表1中的數(shù)據(jù)以及公式(6~10),選擇初始變量X(0)=[2.3,500,-0.12,-0.289,0.03],可得比例風(fēng)險(xiǎn)模型的參數(shù)估計(jì)結(jié)果如表3。

    表3 比例風(fēng)險(xiǎn)模型的參數(shù)估計(jì)結(jié)果

    根據(jù)表2和表3中參數(shù)可得以weibull分布為基底函數(shù)的比例風(fēng)險(xiǎn)模型為:

    根據(jù)失效率函數(shù)得到車門系統(tǒng)的生存函數(shù)為:

    由車門系統(tǒng)的生存函數(shù)可得車門系統(tǒng)的可靠度函數(shù)圖如圖1。

    圖1 車門系統(tǒng)的可靠度函數(shù)圖

    通過車門系統(tǒng)的可靠度函數(shù)計(jì)算其可靠度降為0.95時(shí),可靠壽命為t0.95=170天。

    對(duì)比weibull分布模型和比例風(fēng)險(xiǎn)模型所得車門系統(tǒng)的可靠度函數(shù),可得圖2。

    圖2 PHM和weibull的可靠度函數(shù)圖

    由圖2可知,在初始投入使用的階段檢修對(duì)可靠性的影響較小,但兩個(gè)月后檢修可延長(zhǎng)車門系統(tǒng)的可靠壽命,所以比例風(fēng)險(xiǎn)模型預(yù)測(cè)車門系統(tǒng)的可靠壽命比weibull分布模型更加合理。

    3結(jié)語

    通過分析某地鐵公司2號(hào)線車門系統(tǒng)檢修歷史數(shù)據(jù),引入比例風(fēng)險(xiǎn)模型對(duì)數(shù)據(jù)進(jìn)行處理,得到了車門系統(tǒng)的可靠壽命t0.95=170天。對(duì)比weibull分布可靠性函數(shù)模型可知,把檢修記錄引入可靠性性壽命預(yù)測(cè),可使車門系統(tǒng)的壽命預(yù)測(cè)更符合現(xiàn)場(chǎng)實(shí)際,可為地鐵維修部門提供技術(shù)和管理支持。

    參考文獻(xiàn):

    [1] 時(shí)旭. 地鐵車門系統(tǒng)故障診斷與維修決策的方法研究[D]. 北京:北京交通大學(xué),2009.

    [2] 朱小娟,王建兵,印禎民. 上海地鐵車輛客室車門可靠性技術(shù)研究[J]. 城市軌道交通研究,2006,03:31-34.

    [3] 蔡國(guó)強(qiáng),周莉茗,李熙,等. 基于GO法的城市軌道交通車門系統(tǒng)可靠性分析[J]. 西南交通大學(xué)學(xué)報(bào),2011,02:264-270.

    [4] Thierry Lecomte. Safe and reliable metro platform screen doors control/command systems[J]. Springer-Verlag, 2008: 430-434.

    [5] Cox D R. Regression models and life-tables (with discussion) [J]. Journal of Royal Statistical Society. Series B (Methodological), 1972, 34(2), 187-220.

    [6] Mendes, A.C.; Fard, N., Reliability modeling for appliances using the Proportional Hazard Model, Reliability and Maintainability Symposium (RAMS), 2013: 28-31.

    [7] E.A.Elsayed. 可靠性工程[M]. 第二版. 北京:電子工業(yè)出版社, 2013.

    [8] 戎翔. 民航發(fā)動(dòng)機(jī)健康管理中的壽命預(yù)測(cè)與維修決策方法研究[D]. 南京:南京航空航天大學(xué),2008.

    [9] 張鵬. 基于主成分分析的綜合評(píng)價(jià)研究[D]. 南京:南京理工大學(xué),2004.

    [10] 王文,孟光,劉芳,等. 基于比例風(fēng)險(xiǎn)模型的板級(jí)無鉛焊點(diǎn)跌落壽命分析[J]. 振動(dòng)與沖擊,2011,(3).

    [11] 徐弘博,馮虎田,歐屹. 滾動(dòng)直線導(dǎo)軌副壽命試驗(yàn)方法[J]. 機(jī)械制造與自動(dòng)化,2014,05:65-68.

    [12] 洪東跑,馬小兵,趙宇. 基于比例風(fēng)險(xiǎn)模型的可靠性綜合評(píng)估[J]. 系統(tǒng)工程與電子技術(shù),2010,(10).

    歡迎投稿歡迎刊登廣告歡迎訂閱

    Assessment Method of Metro Door Fault Criticality Based

    on Proportional Hazard Model

    HE Guang-jian1,REN Jin-bao2,XING Zong-yi2

    (1. Guangzhou Metro Corporation, Guangzhou 510030,China;

    2. Nanjing University of Science and Technology, Nanjing 210094, China)

    Abstract:The reliability of Metro Door is directly related to the safety of passengers. The proportional hazard model (PHM) is proposed to evaluate metro door fault criticality. the proportional hazard model is introduced to the description of the relations between reliability and covariates. relevant historical failure data of metro door system, is sifted to obtain the function of the covariates of the model. the log-likelihood function is used to estimate the parameters in the PHM and the Newton-Raphson method is applied to calculating the parameters from the likelihood function. the reliability model is applied to evaluating the lifetime of the metro door. when R is 0.95, the reliable lifetime t0.95is 170. Compared with the Weibull distribution, the PHM is more adaptive and robust and the impact of many covariates is considered simultaneously.

    Keywords:metro door reliability; proportional hazard model; likelihood function; newton-raphson

    基金項(xiàng)目:蘇州市科技發(fā)展 (SYG201258)

    收稿日期:2015-03-16

    中圖分類號(hào):U231

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

    文章編號(hào):1671-5276(2015)03-0187-04

    作者簡(jiǎn)介:何廣堅(jiān)(1960-),男,廣東順德人,副總工程師,研究方向:自動(dòng)控制。

    国产一区有黄有色的免费视频| 国产综合精华液| 满18在线观看网站| 精品第一国产精品| 国产免费现黄频在线看| 免费黄频网站在线观看国产| www日本在线高清视频| 国产亚洲午夜精品一区二区久久| 精品福利永久在线观看| 亚洲精品久久成人aⅴ小说| 久久久精品区二区三区| 中文字幕av电影在线播放| 国产男人的电影天堂91| 日本av免费视频播放| 午夜影院在线不卡| 亚洲欧美中文字幕日韩二区| 国产亚洲午夜精品一区二区久久| av福利片在线| av.在线天堂| 成年人免费黄色播放视频| 国产毛片在线视频| 国产精品久久久av美女十八| a级毛色黄片| av不卡在线播放| av卡一久久| 亚洲色图综合在线观看| 中文精品一卡2卡3卡4更新| 亚洲伊人久久精品综合| 亚洲欧洲国产日韩| 久久久精品免费免费高清| 熟女人妻精品中文字幕| 我要看黄色一级片免费的| 精品国产乱码久久久久久小说| 亚洲第一区二区三区不卡| 不卡视频在线观看欧美| 亚洲欧美一区二区三区国产| 亚洲一级一片aⅴ在线观看| 视频区图区小说| 免费高清在线观看日韩| tube8黄色片| 18禁动态无遮挡网站| 成人二区视频| 狂野欧美激情性bbbbbb| 曰老女人黄片| 热re99久久国产66热| 亚洲色图综合在线观看| 观看美女的网站| 国产视频首页在线观看| 伦理电影免费视频| 午夜激情久久久久久久| 国产av精品麻豆| 老女人水多毛片| 中文字幕精品免费在线观看视频 | 在线 av 中文字幕| 欧美日韩视频高清一区二区三区二| 一级毛片 在线播放| 国产极品天堂在线| 香蕉丝袜av| av视频免费观看在线观看| 久久久久久久大尺度免费视频| 久久久久国产精品人妻一区二区| 欧美人与性动交α欧美软件 | 国产成人精品婷婷| 亚洲av日韩在线播放| 欧美97在线视频| 久久久久久人妻| 亚洲欧美清纯卡通| 久久综合国产亚洲精品| 一级毛片黄色毛片免费观看视频| 国产免费又黄又爽又色| 亚洲国产精品国产精品| 午夜激情av网站| 国产 精品1| 91精品伊人久久大香线蕉| 夫妻午夜视频| 亚洲精品色激情综合| 色吧在线观看| 亚洲国产精品专区欧美| 国产成人91sexporn| 日本vs欧美在线观看视频| 国产麻豆69| 桃花免费在线播放| 精品酒店卫生间| 国产探花极品一区二区| 在线免费观看不下载黄p国产| 久久午夜综合久久蜜桃| 欧美性感艳星| 日本与韩国留学比较| 91精品三级在线观看| 亚洲婷婷狠狠爱综合网| 黑人高潮一二区| www.色视频.com| 亚洲美女黄色视频免费看| 精品福利永久在线观看| 高清黄色对白视频在线免费看| 久久人妻熟女aⅴ| 寂寞人妻少妇视频99o| 国产成人精品久久久久久| 满18在线观看网站| 美女大奶头黄色视频| 国产精品 国内视频| 久久这里只有精品19| 国产熟女欧美一区二区| 午夜视频国产福利| 婷婷成人精品国产| 欧美国产精品va在线观看不卡| 18禁在线无遮挡免费观看视频| 亚洲成国产人片在线观看| 日本与韩国留学比较| 亚洲第一区二区三区不卡| 亚洲中文av在线| 在线精品无人区一区二区三| 久热这里只有精品99| 国产高清三级在线| 高清毛片免费看| 内地一区二区视频在线| 在线观看www视频免费| 国产男人的电影天堂91| 永久网站在线| 日韩,欧美,国产一区二区三区| 午夜福利,免费看| 少妇的逼水好多| 高清视频免费观看一区二区| 久久久久精品人妻al黑| 五月开心婷婷网| 伦理电影大哥的女人| 女人久久www免费人成看片| 久热这里只有精品99| 久久精品国产亚洲av天美| 男人舔女人的私密视频| 欧美日韩亚洲高清精品| 晚上一个人看的免费电影| 国产精品一二三区在线看| 国产69精品久久久久777片| 久久人人97超碰香蕉20202| av黄色大香蕉| av福利片在线| 男女啪啪激烈高潮av片| 免费人妻精品一区二区三区视频| 超色免费av| 日本午夜av视频| 国产男人的电影天堂91| 国产黄色免费在线视频| 久久久久精品久久久久真实原创| 人体艺术视频欧美日本| 国产免费一级a男人的天堂| 国产精品国产av在线观看| 丝袜喷水一区| 精品视频人人做人人爽| 国产成人欧美| 日韩欧美一区视频在线观看| 激情视频va一区二区三区| 亚洲av免费高清在线观看| 热99国产精品久久久久久7| 免费观看性生交大片5| 国产精品女同一区二区软件| 黄色毛片三级朝国网站| 热re99久久国产66热| 人妻少妇偷人精品九色| 在线观看www视频免费| 久久久久久人妻| 日韩欧美一区视频在线观看| 久久久久久久久久久免费av| 国产黄色免费在线视频| 中文精品一卡2卡3卡4更新| 91精品国产国语对白视频| av在线播放精品| 日韩中文字幕视频在线看片| 如日韩欧美国产精品一区二区三区| 制服丝袜香蕉在线| 国产av精品麻豆| 国产男女超爽视频在线观看| 免费高清在线观看日韩| 女人久久www免费人成看片| 看十八女毛片水多多多| 国产成人精品福利久久| 大香蕉久久网| 女性生殖器流出的白浆| 看十八女毛片水多多多| 中文字幕另类日韩欧美亚洲嫩草| 五月开心婷婷网| 久久久久久久久久成人| 国产精品不卡视频一区二区| 国产精品不卡视频一区二区| 中国美白少妇内射xxxbb| 久久这里有精品视频免费| 十八禁网站网址无遮挡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品自拍成人| 国产欧美日韩综合在线一区二区| 一本—道久久a久久精品蜜桃钙片| 极品人妻少妇av视频| 成人国语在线视频| av又黄又爽大尺度在线免费看| 91久久精品国产一区二区三区| 精品一区在线观看国产| 黄色怎么调成土黄色| av国产久精品久网站免费入址| 久久久久网色| 蜜桃国产av成人99| 亚洲国产精品国产精品| 欧美老熟妇乱子伦牲交| 亚洲av电影在线进入| 亚洲精品久久久久久婷婷小说| 色婷婷av一区二区三区视频| 五月玫瑰六月丁香| 天美传媒精品一区二区| 咕卡用的链子| 亚洲欧美日韩另类电影网站| 久久精品熟女亚洲av麻豆精品| 韩国精品一区二区三区 | 久久鲁丝午夜福利片| 日韩人妻精品一区2区三区| 九九在线视频观看精品| 午夜av观看不卡| 日韩av不卡免费在线播放| 亚洲伊人久久精品综合| 国产一区二区在线观看日韩| 日韩大片免费观看网站| 亚洲欧洲国产日韩| 国产麻豆69| 久久午夜综合久久蜜桃| 黄色怎么调成土黄色| 最后的刺客免费高清国语| 最近的中文字幕免费完整| 国产成人欧美| 久热这里只有精品99| 亚洲经典国产精华液单| 亚洲精品乱码久久久久久按摩| 成人国产麻豆网| xxx大片免费视频| 国语对白做爰xxxⅹ性视频网站| 亚洲国产成人一精品久久久| 亚洲精品一区蜜桃| 国产乱人偷精品视频| 熟妇人妻不卡中文字幕| 九色亚洲精品在线播放| 青春草国产在线视频| 天天躁夜夜躁狠狠躁躁| 欧美xxⅹ黑人| 国产男女超爽视频在线观看| 交换朋友夫妻互换小说| 国产 一区精品| 18禁裸乳无遮挡动漫免费视频| av一本久久久久| 国产av国产精品国产| 夜夜骑夜夜射夜夜干| 51国产日韩欧美| 中文字幕亚洲精品专区| 午夜免费男女啪啪视频观看| 女性生殖器流出的白浆| 成人毛片a级毛片在线播放| 亚洲av中文av极速乱| 最黄视频免费看| 欧美国产精品一级二级三级| 亚洲欧洲精品一区二区精品久久久 | 国产不卡av网站在线观看| 亚洲av电影在线进入| a级毛片在线看网站| 亚洲欧美一区二区三区黑人 | 另类精品久久| 久久国产精品男人的天堂亚洲 | av又黄又爽大尺度在线免费看| 9色porny在线观看| www日本在线高清视频| 久久久久久久大尺度免费视频| 久久毛片免费看一区二区三区| 免费高清在线观看日韩| 大香蕉久久成人网| 欧美另类一区| 一级毛片我不卡| a级毛色黄片| 国产av一区二区精品久久| 天天躁夜夜躁狠狠久久av| 在线看a的网站| 啦啦啦在线观看免费高清www| 一边亲一边摸免费视频| 国产精品99久久99久久久不卡 | 最近中文字幕2019免费版| 91国产中文字幕| 人人妻人人澡人人看| 中文字幕最新亚洲高清| 久久97久久精品| 9热在线视频观看99| 亚洲成人一二三区av| 精品一品国产午夜福利视频| 欧美老熟妇乱子伦牲交| 亚洲欧美一区二区三区黑人 | 亚洲熟女精品中文字幕| 国产高清国产精品国产三级| 色婷婷久久久亚洲欧美| 涩涩av久久男人的天堂| 日韩熟女老妇一区二区性免费视频| 一区二区三区精品91| √禁漫天堂资源中文www| videossex国产| 亚洲国产欧美日韩在线播放| 天天影视国产精品| 久久国产精品男人的天堂亚洲 | 五月玫瑰六月丁香| 777米奇影视久久| 国产精品一区二区在线不卡| 国产亚洲精品久久久com| 如何舔出高潮| 80岁老熟妇乱子伦牲交| 人人妻人人澡人人爽人人夜夜| 少妇被粗大猛烈的视频| 午夜福利影视在线免费观看| 日本欧美视频一区| 18禁国产床啪视频网站| av.在线天堂| 草草在线视频免费看| 9191精品国产免费久久| 久久精品久久久久久久性| 欧美97在线视频| 九草在线视频观看| 18在线观看网站| 一本色道久久久久久精品综合| 精品国产国语对白av| 成人二区视频| 久久久欧美国产精品| 纵有疾风起免费观看全集完整版| 色婷婷av一区二区三区视频| 99久久综合免费| 在线亚洲精品国产二区图片欧美| 日韩,欧美,国产一区二区三区| 亚洲色图 男人天堂 中文字幕 | 99久久精品国产国产毛片| 日韩制服丝袜自拍偷拍| 亚洲成av片中文字幕在线观看 | 99久久人妻综合| 日本黄色日本黄色录像| 欧美日韩视频精品一区| 观看美女的网站| 少妇 在线观看| 亚洲三级黄色毛片| 国产免费福利视频在线观看| 日韩大片免费观看网站| 午夜91福利影院| 亚洲国产欧美日韩在线播放| 午夜福利,免费看| 91aial.com中文字幕在线观看| 男女边吃奶边做爰视频| 51国产日韩欧美| 国产午夜精品一二区理论片| 午夜av观看不卡| 一级黄片播放器| 国产毛片在线视频| 插逼视频在线观看| 亚洲精品色激情综合| 日本av手机在线免费观看| 女性被躁到高潮视频| 欧美日韩亚洲高清精品| 精品一品国产午夜福利视频| 精品熟女少妇av免费看| 观看av在线不卡| 看免费av毛片| a 毛片基地| 黄色一级大片看看| 成人毛片a级毛片在线播放| av一本久久久久| 人妻系列 视频| 久久精品熟女亚洲av麻豆精品| 大陆偷拍与自拍| 精品一区二区三区四区五区乱码 | freevideosex欧美| 久久午夜综合久久蜜桃| 天堂俺去俺来也www色官网| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人免费无遮挡视频| 久久久久人妻精品一区果冻| 亚洲国产最新在线播放| 精品99又大又爽又粗少妇毛片| 国产精品人妻久久久久久| 亚洲国产日韩一区二区| 欧美日韩视频精品一区| 综合色丁香网| 欧美3d第一页| 日韩成人伦理影院| freevideosex欧美| 国产日韩一区二区三区精品不卡| 久久精品夜色国产| 亚洲精华国产精华液的使用体验| 久久久亚洲精品成人影院| 免费久久久久久久精品成人欧美视频 | 午夜久久久在线观看| a级毛片黄视频| 夜夜爽夜夜爽视频| 国产女主播在线喷水免费视频网站| 欧美精品av麻豆av| 国产亚洲一区二区精品| 国产成人精品福利久久| 最新中文字幕久久久久| 国产成人a∨麻豆精品| 亚洲av在线观看美女高潮| 亚洲精品美女久久av网站| 精品国产一区二区久久| 丰满少妇做爰视频| 亚洲欧美色中文字幕在线| 国产永久视频网站| 黄色怎么调成土黄色| 亚洲国产精品国产精品| 国产福利在线免费观看视频| 亚洲欧洲日产国产| 青春草国产在线视频| 国产毛片在线视频| 免费播放大片免费观看视频在线观看| 卡戴珊不雅视频在线播放| 国产一区有黄有色的免费视频| 亚洲激情五月婷婷啪啪| 日韩成人av中文字幕在线观看| 女人被躁到高潮嗷嗷叫费观| 色网站视频免费| √禁漫天堂资源中文www| 欧美 日韩 精品 国产| 精品国产露脸久久av麻豆| 欧美另类一区| 18禁观看日本| 欧美xxⅹ黑人| 两个人看的免费小视频| 老熟女久久久| 久久久亚洲精品成人影院| 免费观看无遮挡的男女| 又黄又爽又刺激的免费视频.| 久久99精品国语久久久| 亚洲四区av| av在线老鸭窝| 午夜福利在线观看免费完整高清在| 亚洲av.av天堂| 女人被躁到高潮嗷嗷叫费观| 男女边摸边吃奶| 老司机亚洲免费影院| 国产精品.久久久| 美女xxoo啪啪120秒动态图| 自线自在国产av| 90打野战视频偷拍视频| 亚洲欧洲国产日韩| av卡一久久| 久久久久久久久久成人| 制服诱惑二区| 日本爱情动作片www.在线观看| 国产日韩欧美亚洲二区| 香蕉丝袜av| 国产av国产精品国产| 高清毛片免费看| 久久国产精品大桥未久av| √禁漫天堂资源中文www| 51国产日韩欧美| 夫妻性生交免费视频一级片| 少妇的逼好多水| 国产精品久久久久久久久免| 一级毛片黄色毛片免费观看视频| 男女边吃奶边做爰视频| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕大全免费视频 | 精品久久久久久电影网| 人妻少妇偷人精品九色| 久久女婷五月综合色啪小说| 狂野欧美激情性xxxx在线观看| 日本91视频免费播放| 日韩中字成人| 一边摸一边做爽爽视频免费| 天天操日日干夜夜撸| 国产精品久久久久久av不卡| 内地一区二区视频在线| 国产高清三级在线| 亚洲精品,欧美精品| 国语对白做爰xxxⅹ性视频网站| 午夜免费男女啪啪视频观看| 成人国产麻豆网| 亚洲欧美一区二区三区国产| 夫妻性生交免费视频一级片| 青青草视频在线视频观看| 久久久a久久爽久久v久久| 国产极品天堂在线| 在线观看美女被高潮喷水网站| 午夜福利视频在线观看免费| 精品国产一区二区三区久久久樱花| 2018国产大陆天天弄谢| 国产精品久久久久久久久免| 交换朋友夫妻互换小说| 一二三四中文在线观看免费高清| 日韩欧美一区视频在线观看| 久久这里只有精品19| 国产日韩欧美视频二区| 一级毛片黄色毛片免费观看视频| 亚洲,一卡二卡三卡| 男女无遮挡免费网站观看| 久久久久精品久久久久真实原创| 人妻少妇偷人精品九色| 亚洲在久久综合| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕亚洲精品专区| 2022亚洲国产成人精品| 春色校园在线视频观看| 亚洲精品美女久久av网站| 国产精品免费大片| 老熟女久久久| 黄片无遮挡物在线观看| 国产老妇伦熟女老妇高清| 国产白丝娇喘喷水9色精品| 欧美精品av麻豆av| 国产一区有黄有色的免费视频| 久久精品国产自在天天线| 亚洲精品第二区| av.在线天堂| 在线观看免费高清a一片| 中国美白少妇内射xxxbb| 大码成人一级视频| 国产深夜福利视频在线观看| 成人手机av| 下体分泌物呈黄色| 男女免费视频国产| 国产精品久久久久久av不卡| 久久狼人影院| 欧美bdsm另类| 97精品久久久久久久久久精品| 少妇高潮的动态图| 亚洲欧洲精品一区二区精品久久久 | 91在线精品国自产拍蜜月| 亚洲欧美色中文字幕在线| 最黄视频免费看| 美女国产高潮福利片在线看| 性高湖久久久久久久久免费观看| 亚洲欧洲精品一区二区精品久久久 | 美女国产高潮福利片在线看| 欧美最新免费一区二区三区| 亚洲情色 制服丝袜| 精品一区二区三卡| 美女国产视频在线观看| 欧美成人精品欧美一级黄| 亚洲精品乱码久久久久久按摩| 久久国内精品自在自线图片| 国产精品一区二区在线观看99| 看免费av毛片| 国产精品人妻久久久久久| 大片电影免费在线观看免费| 亚洲国产精品国产精品| 波多野结衣一区麻豆| 少妇熟女欧美另类| 亚洲一码二码三码区别大吗| 又黄又粗又硬又大视频| 香蕉丝袜av| 日韩免费高清中文字幕av| 香蕉丝袜av| 国产精品久久久久成人av| 亚洲经典国产精华液单| 色视频在线一区二区三区| 日韩欧美一区视频在线观看| 国产成人免费观看mmmm| 天天躁夜夜躁狠狠躁躁| 国产精品一国产av| 欧美成人午夜精品| 男女下面插进去视频免费观看 | a级片在线免费高清观看视频| 亚洲成国产人片在线观看| 九九爱精品视频在线观看| 男人操女人黄网站| 国产成人aa在线观看| 男女免费视频国产| 婷婷色综合www| 最近中文字幕2019免费版| 啦啦啦视频在线资源免费观看| 欧美激情国产日韩精品一区| 伦理电影免费视频| 亚洲第一av免费看| 中国三级夫妇交换| 久久久国产欧美日韩av| 啦啦啦在线观看免费高清www| 亚洲第一av免费看| 久久人人爽人人爽人人片va| 99久久人妻综合| 亚洲欧洲国产日韩| 亚洲国产成人一精品久久久| 成人黄色视频免费在线看| 久久99热6这里只有精品| freevideosex欧美| 黄色视频在线播放观看不卡| 男的添女的下面高潮视频| 在现免费观看毛片| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久av不卡| 国产在线免费精品| 亚洲,一卡二卡三卡| 不卡视频在线观看欧美| 成人漫画全彩无遮挡| 欧美激情 高清一区二区三区| 搡老乐熟女国产| 国产精品不卡视频一区二区| av一本久久久久| 国产午夜精品一二区理论片| 久久久精品94久久精品| 亚洲丝袜综合中文字幕| 日本猛色少妇xxxxx猛交久久| 在线观看一区二区三区激情| 成人国产麻豆网| 久久精品国产亚洲av天美| 精品一区二区免费观看| 91成人精品电影| 精品一区二区免费观看| 91成人精品电影| 久久久久久伊人网av| 涩涩av久久男人的天堂| 亚洲av男天堂| 丝袜喷水一区| 亚洲精品色激情综合| 熟妇人妻不卡中文字幕| 亚洲av中文av极速乱| av不卡在线播放| 亚洲欧美一区二区三区黑人 | 2018国产大陆天天弄谢| 国产精品蜜桃在线观看| 日韩av不卡免费在线播放| 人妻少妇偷人精品九色| 人体艺术视频欧美日本|