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

        ?

        基于NSGA2算法的ZPW2000A軌道電路維修策略優(yōu)化研究

        2018-10-08 07:59:56王瑞峰陶榮杰
        關(guān)鍵詞:維修策略軌道電路可靠性

        王瑞峰,陶榮杰

        基于NSGA2算法的ZPW2000A軌道電路維修策略優(yōu)化研究

        王瑞峰,陶榮杰

        (蘭州交通大學(xué) 自動(dòng)化與電氣工程學(xué)院,甘肅 蘭州 730070)

        針對(duì)ZPW2000A軌道電路傳統(tǒng)維修的低可靠性和高維修費(fèi)用問題,提出ZPW2000A維修策略多目標(biāo)優(yōu)化模型,對(duì)其維修策略進(jìn)行優(yōu)化。該模型采用求不交化最小路集方法求ZPW2000A的可靠度;考慮維修活動(dòng)對(duì)系統(tǒng)可靠性和維修費(fèi)用的影響,推導(dǎo)出ZPW2000A系統(tǒng)可靠度和維修費(fèi)用表達(dá)式。將系統(tǒng)的可靠性和維修費(fèi)用作為優(yōu)化目標(biāo),采用非劣排序遺傳算法(NSGA-II)進(jìn)行優(yōu)化計(jì)算,為決策者提供多種選擇,并與傳統(tǒng)方法進(jìn)行對(duì)比,驗(yàn)證了該方法的有效性。

        ZPW2000A;維修策略;不交化最小路集;NSGA-II;Pareto非劣解集

        ZPW2000A軌道電路作為信號(hào)的重要基礎(chǔ)設(shè)備,其可靠性對(duì)于行車安全極為重要。長(zhǎng)期以來,軌道電路設(shè)備維修面對(duì)的是“期望設(shè)備高性能的運(yùn)轉(zhuǎn)與維護(hù)資源有限”的矛盾。而現(xiàn)階段,我國鐵路針對(duì)信號(hào)設(shè)備的計(jì)劃維修主要以其服役時(shí)間作為根據(jù),然而這會(huì)使部分設(shè)備組件沒有將其功能充分發(fā)揮,造成維修資源的浪費(fèi),維修過量則影響設(shè)備的可靠性,造成安全隱患。因此,對(duì)ZPW2000A軌道電路維修策略進(jìn)行優(yōu)化研究,平衡維修資源和設(shè)備可靠性之間的矛盾,以最低的維修資源來達(dá)到對(duì)設(shè)備預(yù)期的性能要求,達(dá)到效益最優(yōu),非常具有現(xiàn)實(shí)意義。近年來,很多學(xué)者致力于軌道電路設(shè)備的維修策略和可靠性的研究。馬濤[1]在使用故障模式和影響分析方法對(duì)設(shè)備故障進(jìn)行分類的基礎(chǔ)上結(jié)合RCM邏輯決斷圖對(duì)設(shè)備維修策略進(jìn)行決策分析研究。米根鎖等[2]在對(duì)軌道電路設(shè)備剩余壽命預(yù)測(cè)的基礎(chǔ)上來調(diào)整維修方法。王文斌等[3]將故障模式影響分析方法和故障樹分析方法相結(jié)合來研究ZPW2000A系統(tǒng)的可靠性。由此可以看出,對(duì)于ZPW2000A軌道電路設(shè)備維修策略和可靠性的研究大多停留在以可靠性作為單一優(yōu)化目標(biāo)或者研究對(duì)象,并沒有考慮不同維修活動(dòng)對(duì)設(shè)備的綜合影響和平衡設(shè)備可靠性和經(jīng)濟(jì)性之間的關(guān)系。本文考慮不同維修活動(dòng)對(duì)ZPW2000A可靠性和維修費(fèi)用的影響,建立以ZPW2000A系統(tǒng)平均可靠性和維修費(fèi)用為優(yōu)化目標(biāo)的維修策略優(yōu)化模型。采用非劣排序遺傳算法(NSGA-II)進(jìn)行優(yōu)化求解[4?8],該算法通過模仿生物遺傳進(jìn)化過程中選擇、雜交和變異來完成對(duì)問題最優(yōu)解的搜索過程,得到一系列Pareto非劣解集,為決策者提供更多的最優(yōu)選擇。

        1 ZPW2000A維修策略多目標(biāo)優(yōu)化模型建立

        1.1 優(yōu)化目標(biāo)

        ZPW2000A維修的核心問題是其可靠性,在此基礎(chǔ)上追求經(jīng)濟(jì)性,從而達(dá)到效益最優(yōu)。本文將ZPW2000A的可靠性和維修費(fèi)用作為2個(gè)優(yōu)化 目標(biāo)。

        在建立模型之前,首先對(duì)模型進(jìn)行如下假設(shè):

        1) 設(shè)備的失效率服從威布爾分布;

        2) 當(dāng)設(shè)備運(yùn)行到維修周期時(shí),必須進(jìn)行維修;

        3) 當(dāng)設(shè)備發(fā)生故障后進(jìn)行維修,修復(fù)后的設(shè)備只能恢復(fù)到故障前一刻的狀態(tài),對(duì)設(shè)備可靠度不產(chǎn)生影響。

        1.1.1 ZPW2000A可靠性

        求ZPW2000A系統(tǒng)可靠性,應(yīng)對(duì)其系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)關(guān)系進(jìn)行分析,采用鄰接矩陣法求其最小路集[9],再對(duì)最小路集進(jìn)行不交化處理得到不交化最小路集,將每個(gè)設(shè)備的可靠度代入到系統(tǒng)不交化最小路集中,從而得到系統(tǒng)的可靠性表達(dá)式。

        為了降低復(fù)雜系統(tǒng)的求解難度和計(jì)算量,通常要對(duì)最小路集進(jìn)行不交化處理。若要求系統(tǒng)正常工作,則至少保證有一個(gè)最小路集正常。規(guī)定系統(tǒng)正常工作事件為:

        其中:為最小路集個(gè)數(shù)。通過集合運(yùn)算法則不交化得:

        具體不交化方法為:

        1) 將最小路集按照所含元素個(gè)數(shù)進(jìn)行排序和編號(hào);

        2) 假定所有最小路集含個(gè)元素,有個(gè)最小路集,建立矩陣×,最小路集排序后,對(duì)每個(gè)元素位置進(jìn)行標(biāo)記,包含為1,不包含則為0;

        3) 選出其中一個(gè)最小路集最為參照路集,把排在其前面的所有最小路集當(dāng)作比較路集,刪除比較路集和參照路集中一樣的元素,并對(duì)每個(gè)比較路集剩余的元素取反。將取反后的元素相乘并簡(jiǎn)化,然后與參照路集相乘,所得為不交化后的結(jié)果,將所有不交化結(jié)果相加即為不交化最小路集。

        通過上述方法可求ZPW2000A系統(tǒng)不交化最小路集。為了簡(jiǎn)化計(jì)算,忽略如鋼軌等因素的影響,系統(tǒng)中發(fā)送器為+1冗余,取值為1,接收器為雙機(jī)熱備并聯(lián)。ZPW2000A系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)如圖1 所示。

        圖1 ZPW2000A系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)

        由圖1可得ZPW2000A系統(tǒng)的鄰接矩陣:

        對(duì)最小路集進(jìn)行不交化處理,將得到的每個(gè)最小路集不交化結(jié)果相加,得到最終的ZPW2000A系統(tǒng)不交化最小路集為:

        在ZPW2000A系統(tǒng)中,不同設(shè)備性能變化水平不同,因此,針對(duì)不同設(shè)備及設(shè)備所處不同質(zhì)量狀態(tài)選擇合適的維修方式極為重要。目前,ZPW2000A的維修活動(dòng)主要可分為修復(fù)性維修和定期預(yù)防性維修。修復(fù)性維修指設(shè)備故障發(fā)生后進(jìn)行的維修活動(dòng),具有不可預(yù)見性。預(yù)防性維修是在便于調(diào)度管理和有準(zhǔn)備的前提下進(jìn)行的一種維修活動(dòng)。而修復(fù)性維修和預(yù)防性維修都是通過具體的維修方式來實(shí)現(xiàn)的。本文主要考慮3種維修方式,如表1所示。

        表1 維修方式

        由于修復(fù)性維修只會(huì)恢復(fù)設(shè)備原有的功能,不會(huì)改變?cè)O(shè)備故障率,所以不考慮其對(duì)設(shè)備可靠性的影響。服從二參數(shù)的威布爾ZPW2000A設(shè)備組件可靠性分布函數(shù)為

        若采用1維修方式,對(duì)其性能的改善為

        若采用3維修方式,對(duì)其性能的改善為:

        式(12)表示設(shè)備變?yōu)槿隆?/p>

        將式(7)代入式(6)可得系統(tǒng)的可靠度為

        由式(14)可對(duì)某維修策略下系統(tǒng)平均可靠性進(jìn)行計(jì)算,將其作為維修策略優(yōu)化目標(biāo)之一。

        1.1.2 系統(tǒng)維修費(fèi)用

        則在這段時(shí)間內(nèi),整個(gè)ZPW2000A系統(tǒng)修復(fù)性維修所需費(fèi)用如式(16)所示

        則ZPW2000A系統(tǒng)維修費(fèi)用如式(18)所示。

        1.2 約束條件

        因此得到ZPW2000A系統(tǒng)維修策略多目標(biāo)優(yōu)化模型為:

        優(yōu)化目標(biāo)如式(20)所示

        約束條件如式(19)所示

        2 NSGA-Ⅱ算法

        NSGA-Ⅱ算法是Deb等學(xué)者于2000年在NSGA算法的基礎(chǔ)之上提出來的,該算法能夠得到一系列分布均勻、多樣性較好的最優(yōu)解集,適用于多個(gè)目標(biāo)的優(yōu)化問題,且相對(duì)成熟,廣泛應(yīng)用于橋梁、核電等領(lǐng)域的維修策略優(yōu)化中[9?10]。

        圖2 NSGA-Ⅱ流程圖

        3 實(shí)例驗(yàn)證

        3.1 編碼和解碼

        采用二進(jìn)制編碼,該方法編碼和解碼操作簡(jiǎn)單容易實(shí)現(xiàn)。每?jī)晌淮硪淮尉S修活動(dòng)信息,每次維修活動(dòng)可以選擇不維修、保養(yǎng)、局部更換和全部更換4種級(jí)別,對(duì)應(yīng)的編碼表示如表2所示。

        若在優(yōu)化年限內(nèi)進(jìn)行了N次預(yù)防性維修,那么一個(gè)設(shè)備維修信息編碼長(zhǎng)度為2N個(gè),則系統(tǒng)的每一個(gè)維修計(jì)劃長(zhǎng)度為22N,即每個(gè)維修計(jì)劃可用長(zhǎng)度為22N的二進(jìn)制數(shù)表示。維修策略編碼后參與遺傳算法的優(yōu)化計(jì)算,得到優(yōu)化結(jié)果后還需進(jìn)行編碼的逆操作?解碼,解碼是按照編碼的方法將染色體進(jìn)行對(duì)應(yīng)的解釋來完成的。

        表2 維修方式編碼

        3.2 參數(shù)設(shè)置

        為簡(jiǎn)化計(jì)算,規(guī)定區(qū)段長(zhǎng)度為600 m,補(bǔ)償電容個(gè)數(shù)為7。不同維修方式的影響因子1,2和3設(shè)定為0.4,0.7,1,交叉算子為0.8,變異算子為0.2,種群數(shù)目設(shè)置為200,進(jìn)化代數(shù)為500,優(yōu)化年限為12 a。其他參數(shù)設(shè)置如表3所示。

        表3 系統(tǒng)參數(shù)

        3.3 優(yōu)化過程

        將初始種群個(gè)體代入系統(tǒng)可靠性和維修費(fèi)用目標(biāo)函數(shù)中,得到每個(gè)個(gè)體對(duì)應(yīng)的可靠度和維修費(fèi)用。將帶有維修信息的種群放進(jìn)NSGA2算法中進(jìn)行優(yōu)化,最后得到Pareto最優(yōu)解集。

        圖3為12次和16次維修次數(shù)下,優(yōu)化結(jié)果對(duì)比??梢钥闯?,隨著系統(tǒng)平均可靠性的增加,維修費(fèi)用也在不斷增加,當(dāng)超過一定極值后,可靠性的增加會(huì)使維修費(fèi)用大幅增加。不同維修次數(shù)對(duì)可靠性和維修費(fèi)用產(chǎn)生不同影響,當(dāng)維修次數(shù)較少時(shí),要達(dá)到較高的可靠度,維修費(fèi)用會(huì)越高。維修計(jì)劃的決策者可根據(jù)自身的實(shí)際情況選擇合適的維修策略。以圖3中A點(diǎn)為例,該點(diǎn)所對(duì)應(yīng)的維修策略進(jìn)行了16次維修,在優(yōu)化年限12 a內(nèi),投入75萬元維修費(fèi)用,ZPW2000A的系統(tǒng)平均可靠性最低為0.785。該點(diǎn)所對(duì)應(yīng)維修計(jì)劃如表4所示。根據(jù)某鐵路局某段在優(yōu)化年限內(nèi)的天窗維修數(shù)據(jù)及故障維修數(shù)據(jù)對(duì)其維修費(fèi)用和可靠度采用本文所建模型進(jìn)行計(jì)算,A策略年度累計(jì)維修成本和傳統(tǒng)維修維修策略年度累計(jì)維修成本對(duì)比如圖4所示。

        圖3 不同維修次數(shù)的優(yōu)化結(jié)果對(duì)比

        表4 ZPW2000A維修計(jì)劃

        圖4 2種維修策略年度累計(jì)維修成本對(duì)比圖

        由圖4可以看出,在A維修策略下維修耗費(fèi)成本低于傳統(tǒng)維修策略。表5為在優(yōu)化年限內(nèi)A策略下系統(tǒng)可靠度和傳統(tǒng)維修策略下系統(tǒng)可靠度的 對(duì)比。

        表5 2種策略可靠度對(duì)比

        由表5可以看出,A策略下系統(tǒng)可靠度優(yōu)于傳統(tǒng)維修策略。通過以上對(duì)比分析,說明本文所用優(yōu)化方法得到的維修策略能夠合理控制維修成本,達(dá)到的維修效果也優(yōu)于傳統(tǒng)維修策略,證明了該方法的有效性。

        4 結(jié)論

        1) 采用威布爾分布來描述ZPW2000A設(shè)備的失效過程,在考慮不同維修活動(dòng)的條件下,利用求不交化最小路集方法對(duì)ZPW2000A系統(tǒng)的可靠性進(jìn)行推導(dǎo)。為ZPW2000A可靠性領(lǐng)域進(jìn)一步研究提供借鑒。

        2) 建立以可靠性和維修費(fèi)用為優(yōu)化目標(biāo),以每個(gè)設(shè)備的可靠度下限為約束條件的多目標(biāo)優(yōu)化模型,用非劣排序遺傳算法NSGA2該模型進(jìn)行優(yōu)化,得到的一系列最優(yōu)解能為決策制定者制定合理科學(xué)的維修計(jì)劃提供依據(jù)。

        3) 通過對(duì)仿真結(jié)果的分析可知,通過維修使系統(tǒng)的可靠性增加時(shí),所需維修成本也隨之增加,超過一定極值后維修成本會(huì)大幅增加。而過高的可靠度也會(huì)要求維修次數(shù)的增加,對(duì)鐵路運(yùn)營(yíng)和調(diào)度產(chǎn)生影響。因此,決策者應(yīng)根據(jù)運(yùn)營(yíng)時(shí)間段以及其他自身實(shí)際情況選擇適當(dāng)?shù)木S修策略。

        4) 通過與傳統(tǒng)維修策略相比較,本文所用方法在控制維修成本和維修效果方面優(yōu)于傳統(tǒng)維修 策略。

        [1] 馬濤. 以可靠性為中心的ZPW-2000A無絕緣軌道電路設(shè)備維修研究[D]. 蘭州: 蘭州交通大學(xué), 2013. MA Tao. Research on reliability-centered maintence of ZPW-2000A track circuit equipment[D]. Lanzhou: Lanzhou Jiaotong University, 2013.

        [2] 米根鎖, 張鳳霞, 魏蕾. 基于剩余壽命的鐵路軌道電路調(diào)整型維修方法研究[J]. 鐵道學(xué)報(bào), 2015, 37(4): 69?74. MI Gensuo, ZHANG Fengxia, WEI Lei. Research on method of modulation maintenance in railway track circuit based on residual life[J]. Journal of the China Railway Society, 2015, 37(4): 69?74.

        [3] 王文斌, 蘇宏升. 基于FMEA和FTA的ZPW-2000A軌道電路可靠性分析[J]. 計(jì)算機(jī)工程與應(yīng)用, 2013, 49(21): 245?249. WANG Wenbin, SU Hongsheng. Analysis on reliability of ZPW-2000A track circuit based on FMEA and FTA[J]. Computer Engineering and Applications, 2013, 49(21): 245?249.

        [4] 劉洪, 李榮, 葛少云, 等. 考慮出租車隨機(jī)行為特性及路網(wǎng)行程時(shí)間可靠性的充電站多目標(biāo)規(guī)劃[J]. 電網(wǎng)技術(shù), 2016, 40(2): 433?441. LIU Hong, LI Rong, GE Shaoyun, et al. Multiple objective planning for electric taxi charging station considering random-probability behavioral trait of taxi and travel time reliability[J]. Power System Technology, 2016, 40(2): 433?441.

        [5] 肖懷碩, 李清泉, 張巖, 等. 基于非支配排序遺傳改進(jìn)算法的故障限流器的優(yōu)化配置[J]. 電網(wǎng)技術(shù), 2016, 40(8): 2444?2450. XIAO Huaishuo, LI Qingquan, ZHANG Yan, et al. Optimization of fault current limiter configuration based on improved NSGA-2[J]. Power System Technology, 2016, 40(8): 2444?2450.

        [6] 姜惠蘭, 安星, 王亞微, 等. 基于改進(jìn)NSGA2算法的考慮風(fēng)機(jī)接入電能質(zhì)量的多目標(biāo)電網(wǎng)規(guī)劃[J]. 中國電機(jī)工程學(xué)報(bào), 2015, 35(21): 5405?5411. JIANG Huilan, AN Xing, WANG Yawei, et al. Improved NSGA2 algorithm based multi-objective planning of power grid with wind farm considering power quality[J]. Proceedings of the CSEE, 2015, 35(21): 5405?5411.

        [7] 王瓏, 王同光, 羅源. 改進(jìn)的NSGA-Ⅱ算法研究風(fēng)力機(jī)葉片多目標(biāo)優(yōu)化[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2011, 32(6): 693?701. WANG Long, WANG Tongguang, LUO Yuan. Improved NSGA-Ⅱ in multi-objective optimization studies of wind turbine blades[J]. Applied Mathematics and Mechanics, 2011, 32(6): 693?701.

        [8] 聶瑞, 章衛(wèi)國, 李廣文, 等. 基于Tent映射的自適應(yīng)混沌混合多目標(biāo)遺傳算法[J]. 北京航空航天大學(xué)學(xué)報(bào), 2012, 38(8): 1010?1016. NIE Rui, ZHANG Weiguo, LI Guangwen, et al. Adaptive chaos hybrid multi-objective genetic algorithm based on the Tent map[J]. Journal of Beijing University of Aeronautics and Astronautics, 2012, 38(8): 1010?1016.

        [9] 呂言, 劉井泉, 曾聿赟. 基于非劣排序遺傳算法的核電廠維修決策多目標(biāo)優(yōu)化方法研究[J]. 核動(dòng)力工程, 2017, 38(1): 120?125. Lü Yan, LIU Jingquan, ZENG Yuyun. Research on multi-objective optimization method for maintenance decision of nuclear power plants[J]. Nuclear Power Engineering, 2017, 38(1): 120?125.

        [10] 王曉明, 賀耀北, 李瑜, 等. 基于并行NSGA-Ⅱ算法的橋梁網(wǎng)絡(luò)養(yǎng)護(hù)策略優(yōu)化[J]. 土木工程學(xué)報(bào), 2012, 45(1): 86?91. WANG Xiaoming, HE Yaobei, LI Yu, et al. Optimization of maintenance strategy for bridge networks using parallel NSGA-Ⅱ algorithm[J]. China Civil Engineering Journal, 2012, 45(1): 86?91.

        [11] Kannan S, Baskar S, McCalley J D, et al. Application of NSGA-II algorithm to generation expansion planning[J]. IEEE Transactions on Power Systems, 2009, 24(1): 454? 461.

        [12] Yen G, Lu H. Dynamic multiobjective evolutionary algorithm:Adative cell-based rank and density estimation[J]. IEEE Transactions on Evolutionary Computation, 2003, 7(3): 253?274.

        [13] Lapa, Celso Marcelo F, Claudio Marcio NA Pereira, et al. A model for preventive maintenance planning by genetic algorithms based in cost and reliability[J]. Reliability Engineering & System Safety, 2006, 91(2): 233?240.

        [14] Zio E, Bazzo R. A clustering procedure for reducing the number of representative solutions in the Pareto Front of multi-objective optimization problems[J]. European Journal of Operational Research, 2011, 210(3): 624?634.

        (編輯 陽麗霞)

        Research on optimization of ZPW2000A track circuit maintenance strategy based on NSGA2 algorithm

        WANG Ruifeng, TAO Rongjie

        (College of Automatic & Electrical Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)

        Concerning the problem that the low reliability and high maintenance cost of traditional maintenance of ZPW2000A track circuit, the authors proposed a multi - objective optimization model of ZPW2000A maintenance strategy to optimize the maintenance strategy of ZPW2000A. In this model, the method of disjointed minimal paths was used to obtain the reliability of ZPW2000A. The influence of maintenance activities on system reliability and maintenance cost was taken into account, and the expression of ZPW2000A system reliability and maintenance cost was derived. The reliability and maintenance costs of the ZPW2000A system were considered as optimization objective and the non-inferior ranking genetic algorithm (NSGA-II) was used for calculation. It provides a variety of choices for the decision maker.Compared with the traditional method, the validity of this method was verified.

        ZPW2000A; maintenance strategy optimization; disjointed minimal paths; NSGA-II; Pareto non-inferior solution set

        10.19713/j.cnki.43?1423/u.2018.09.028

        TP207

        A

        1672 ? 7029(2018)09 ? 2394 ? 07

        2017?07?02

        甘肅省自然科學(xué)基金資助項(xiàng)目(1310RJZA046)

        王瑞峰(1966?),女,內(nèi)蒙古卓資人,教授,從事交通運(yùn)輸自動(dòng)化與控制研究;E?mail:137109066@qq.com

        猜你喜歡
        維修策略軌道電路可靠性
        考慮故障相關(guān)性的風(fēng)電機(jī)組維修策略
        裝備維修策略選擇和優(yōu)化研究綜述
        基于HHT及LCS的軌道電路傳輸變化識(shí)別探討
        可靠性管理體系創(chuàng)建與實(shí)踐
        JXG-50S型相敏軌道電路接收器自動(dòng)測(cè)試臺(tái)
        ZPW-2000客專軌道電路掉碼故障分析
        基于狀態(tài)維修的防空反導(dǎo)裝備預(yù)防性維修策略
        電子制作(2017年2期)2017-05-17 03:55:06
        25Hz相敏軌道電路的計(jì)算和仿真
        電氣化鐵道(2016年2期)2016-05-17 03:42:36
        基于可靠性跟蹤的薄弱環(huán)節(jié)辨識(shí)方法在省級(jí)電網(wǎng)可靠性改善中的應(yīng)用研究
        日韩精品中文一区二区三区在线 | 黑人巨大精品欧美一区二区免费| 亚洲综合精品成人| 国产未成女年一区二区| 国内精品少妇久久精品| 日日噜噜夜夜狠狠久久丁香五月| 国产精品成人国产乱| 国产综合精品久久亚洲| 国产精品视频白浆免费看| 欧美性生交大片免费看app麻豆| 人妻献身系列第54部| 无码一区二区三区在线在看| 国产白浆大屁股精品视频拍| av无码国产精品色午夜| 久久久久久免费毛片精品| 精品国产AⅤ一区二区三区V免费| 经典黄色一区二区三区| 日日摸天天摸97狠狠婷婷| 精品人体无码一区二区三区| 国产精品亚洲av网站| 邻居美少妇张开腿让我爽了一夜 | 中国免费av网| 免费国产一区二区视频| 四川发廊丰满老熟妇| 色老头一区二区三区| 男女视频网站免费精品播放| 欧美成人家庭影院| 欧美成人一区二区三区| WWW拍拍拍| 成人自拍一二在线观看| 国产揄拍国产精品| 无码不卡免费一级毛片视频| 国产伦精品一区二区三区| 色欲人妻aaaaaaa无码| 国产人碰人摸人爱视频| 国产天堂av手机在线| 免费一级淫片日本高清| 国产女主播喷水视频在线观看| 97色综合| 亚洲不卡一区二区视频| 曰韩人妻无码一区二区三区综合部 |