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

    基于IBA-RBF的水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)預(yù)測(cè)

    2020-08-22 07:27:30劉丹彤杜志達(dá)
    中國(guó)農(nóng)村水利水電 2020年8期
    關(guān)鍵詞:蝙蝠廠房水電站

    劉丹彤,杜志達(dá)

    (大連理工大學(xué)建設(shè)工程學(xué)部,遼寧 大連 116024)

    0 引 言

    水電站廠房結(jié)構(gòu)具有支撐發(fā)電機(jī)組和水體通過的雙重作用。隨著水電事業(yè)的飛速發(fā)展,在裝機(jī)容量日益擴(kuò)大的同時(shí),由多種因素引發(fā)的廠房結(jié)構(gòu)振動(dòng)問題也日漸突出,威脅著電站的安全運(yùn)行。目前各大型水電站對(duì)機(jī)組振動(dòng)及水流脈動(dòng)情況已建立了較為完善的監(jiān)測(cè)體系,然而受現(xiàn)場(chǎng)條件限制,較少在廠房結(jié)構(gòu)上布設(shè)測(cè)點(diǎn)[1]。為全面地掌握和控制電站的運(yùn)行情況,如何利用有限的監(jiān)測(cè)數(shù)據(jù)預(yù)測(cè)廠房結(jié)構(gòu)的振動(dòng)狀態(tài)和振動(dòng)響應(yīng)成為一項(xiàng)重要課題[2,3]。

    在水力、機(jī)械、電磁等因素的共同作用下,水電站廠房結(jié)構(gòu)的振動(dòng)特性呈現(xiàn)典型的高維非線性和多源不確定性特點(diǎn),為廠房結(jié)構(gòu)振動(dòng)響應(yīng)的計(jì)算分析帶來了不小的挑戰(zhàn)[4,5]。由于各類振源的產(chǎn)生機(jī)理與作用方式各不相同,其動(dòng)荷載大小及幅頻特性亦難以確定,機(jī)組與廠房結(jié)構(gòu)間耦聯(lián)振動(dòng)的存在進(jìn)一步增大了利用有限元及結(jié)構(gòu)動(dòng)力學(xué)方法求解的難度,因此廠房結(jié)構(gòu)的振動(dòng)問題一直是國(guó)內(nèi)外學(xué)術(shù)和工程領(lǐng)域的關(guān)注重點(diǎn)[6]。許多研究者做了大量的工作,練繼建等[7]提出基于BP神經(jīng)網(wǎng)絡(luò)的水電站廠房結(jié)構(gòu)振動(dòng)位移預(yù)測(cè)模型,結(jié)合三峽水電站這一實(shí)例驗(yàn)證了模型的適用性。徐國(guó)賓等[8]利用SSPSO-GRNN模型對(duì)廠頂溢流式水電站廠房結(jié)構(gòu)的振動(dòng)位移進(jìn)行預(yù)測(cè),獲得了在0.004左右波動(dòng)的均方誤差。宋志強(qiáng)等[9]提出了改進(jìn)的螢火蟲算法(IFA),并利用其優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的權(quán)值和閾值,提出水電站廠房結(jié)構(gòu)振動(dòng)加速度預(yù)測(cè)模型IFA-BPNN。苗宗偉等[10]通過對(duì)景洪水電站現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)的分析,證明了廠房結(jié)構(gòu)振動(dòng)與機(jī)組振動(dòng)間的耦聯(lián)作用及相關(guān)性。練繼建等[11]以二灘水電站原型觀測(cè)數(shù)據(jù)為基礎(chǔ),利用粒子群算法(PSO)優(yōu)化最小二乘支持向量機(jī)(LS-SVM)的參數(shù),對(duì)廠房結(jié)構(gòu)振動(dòng)響應(yīng)進(jìn)行了預(yù)測(cè),預(yù)測(cè)誤差在20%以內(nèi)。以上研究說明了通過機(jī)組振動(dòng)和尾水脈動(dòng)的監(jiān)測(cè)數(shù)據(jù)預(yù)測(cè)廠房結(jié)構(gòu)振動(dòng)響應(yīng)的可行性,對(duì)于保障水電站的安全經(jīng)濟(jì)運(yùn)行起到了積極作用,但預(yù)測(cè)精度上的不足,有待于進(jìn)一步完善。

    本文基于對(duì)蝙蝠算法尋優(yōu)策略分析,在全局尋優(yōu)和局部搜索階段分別引入動(dòng)態(tài)擾動(dòng)因子和正切隨機(jī)探索機(jī)制,采用慣性權(quán)重遞減策略對(duì)兩種搜索機(jī)制進(jìn)行平衡。以線性變化的方式更新脈沖音量和脈沖頻率,同時(shí)進(jìn)行種群刷新操作加快尋優(yōu)效率。提出改進(jìn)蝙蝠算法(IBA)并利用其對(duì)徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)(RBF)的擴(kuò)展參數(shù)P進(jìn)行優(yōu)選,基于機(jī)組振動(dòng)和尾水脈動(dòng)的監(jiān)測(cè)數(shù)據(jù)建立水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)預(yù)測(cè)的IBA-RBF模型,對(duì)水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)進(jìn)行了高精度預(yù)測(cè)。

    1 基本蝙蝠算法及其改進(jìn)

    1.1 基本蝙蝠算法

    蝙蝠算法自2010年Yang[12]提出以來,以其收斂快、參數(shù)少、通用性強(qiáng)的優(yōu)勢(shì),廣泛應(yīng)用于數(shù)據(jù)挖掘[13,14]、圖像處理[15,16]、模式識(shí)別[17]等諸多領(lǐng)域,并取得了較好的效果。但作為一種群體智能優(yōu)化算法,蝙蝠算法也存在著易陷入局部極值、后期收斂慢的缺陷。通過迭代過程中脈沖音量和脈沖頻率的適應(yīng)性改變,蝙蝠種群實(shí)現(xiàn)了全局尋優(yōu)與局部搜索兩部分間的動(dòng)態(tài)轉(zhuǎn)換,具體實(shí)現(xiàn)步驟如下:

    fi=fmin+(fmax-fmin)×β

    (1)

    (2)

    (3)

    式中:β是[0,1]上服從均勻分布的隨機(jī)向量;fmin、fmax分別為種群發(fā)射超聲波的最小和最大頻率。

    Step2:局部搜索—接近全局最優(yōu)解的蝙蝠(以隨機(jī)生成的均勻分布的隨機(jī)數(shù)rand1為判斷閾值,若rand1>第i只蝙蝠的脈沖頻率ri)采用以下方式在最優(yōu)個(gè)體附近產(chǎn)生新解xnew:

    xnew=xold+εAt

    (4)

    式中:ε為[-1,1]上均勻分布的隨機(jī)數(shù);At為第t次迭代時(shí)種群脈沖音量的平均值。

    在個(gè)體發(fā)現(xiàn)獵物后,脈沖音量Ai和脈沖頻率ri按下式更新:

    Ait+1=αAit

    (5)

    rit+1=ri0[1-exp(-γt)]

    (6)

    1.2 動(dòng)態(tài)擾動(dòng)因子

    基本蝙蝠算法中蝙蝠個(gè)體每次全局搜索時(shí)都受上一次迭代時(shí)所在位置與當(dāng)前最優(yōu)位置間距離的制約,捕食空間有限,易陷入局部最優(yōu)。借鑒文獻(xiàn)[18]的改進(jìn)思路,引入動(dòng)態(tài)擾動(dòng)因子δ對(duì)上一次迭代時(shí)所在位置進(jìn)行擾動(dòng),以提高種群捕食的靈活性。搜索前期擾動(dòng)力度較大,全局搜索的開拓能力較強(qiáng),隨著捕食過程的進(jìn)行,δ逐漸減小以保證算法收斂。動(dòng)態(tài)擾動(dòng)因子δ和改進(jìn)后的位置更新公式如下:

    (7)

    (8)

    式中:Tmax為最大迭代次數(shù);κ為擾動(dòng)偏離因子,κ∈[0.1,0.9],經(jīng)多次試驗(yàn)確定其取值,本文取κ=0.5;betarnd( )為服從貝塔分布的隨機(jī)數(shù),對(duì)自適應(yīng)減小的δ值進(jìn)行擾動(dòng)以提高全局搜索能力。

    1.3 慣性權(quán)重遞減策略

    蝙蝠算法和粒子群算法都是耦合個(gè)體在解空間的全局搜索與局部搜索實(shí)現(xiàn)尋優(yōu),可參考粒子群算法中的慣性權(quán)重w遞減策略平衡兩種搜索機(jī)制,改進(jìn)后的速度更新公式如式(11)所示:

    c=1/Tmax×ln(wmax/wmin)

    (9)

    w=wmax-(wmax-wmin)×ect

    (10)

    (11)

    式中:wmax、wmin分別為慣性權(quán)重w的最大值和最小值,本文取wmax=0.9,wmin=0.4。

    1.4 正切隨機(jī)探索機(jī)制

    針對(duì)基本蝙蝠算法在局部搜索階段以當(dāng)前最優(yōu)位置為基準(zhǔn)均勻向外探索、局部挖掘能力較弱的缺陷,利用柯西分布函數(shù)[19]按下式對(duì)局部搜索機(jī)制進(jìn)行改進(jìn):

    xnew=xold+Attan[π(λ-0.5)]

    (12)

    1.5 脈沖音量和頻率線性漸變

    基本蝙蝠算法產(chǎn)生新解的條件為:rand1>ri,接受新解的條件為:新解的適應(yīng)度優(yōu)于當(dāng)前最優(yōu)解的適應(yīng)度f(xnew)

    (13)

    (14)

    式中:Amax、Amin、rmax、rmin分別為整個(gè)蝙蝠種群的最大、最小脈沖音量及最大、最小脈沖頻率。

    1.6 種群刷新

    借鑒人工免疫算法(Immune Algorithm)的種群刷新策略,每次迭代用隨機(jī)產(chǎn)生的新個(gè)體替換適應(yīng)度排名后10%的個(gè)體,以加快尋優(yōu)效率。

    2 函數(shù)測(cè)試

    為比較改進(jìn)蝙蝠算法(IBA)與基本蝙蝠算法(BA)的尋優(yōu)性能,以Schaffer、Griewank、Ackley、Sphere 4個(gè)典型函數(shù)的最小化問題進(jìn)行測(cè)試,測(cè)試過程中參數(shù)設(shè)置如表1所示。

    表1 測(cè)試過程算法參數(shù)設(shè)置Tab.1 Test process algorithm parameter settings

    (1)Schaffer函數(shù)。

    f(xi)=0.5+

    (15)

    函數(shù)的測(cè)試區(qū)間為xi∈[-2.048,2.048],在xi=0處取最小值0。

    (2)Griewank函數(shù)。

    (16)

    測(cè)試區(qū)間為xi∈[-600,600],在xi=0處取全局最小值0,該函數(shù)是典型的非線性多模態(tài)函數(shù),難以找到全局最小值。

    (3)Ackley函數(shù)。

    (17)

    函數(shù)的測(cè)試區(qū)間為xi∈[-30,30],在xi=0處取全局最小值0。

    (4)Sphere函數(shù)。

    (18)

    函數(shù)的測(cè)試區(qū)間為xi∈[-100,100],該函數(shù)是典型的單峰函數(shù),在xi=0處取全局最小值0。

    設(shè)置測(cè)試函數(shù)的維數(shù)D=10,利用BA、IBA算法分別對(duì)上述4個(gè)函數(shù)最小值問題進(jìn)行尋優(yōu),將每個(gè)函數(shù)獨(dú)立運(yùn)行20次的結(jié)果進(jìn)行統(tǒng)計(jì)分析,如表2所示。

    表2 函數(shù)測(cè)試結(jié)果Tab.2 Function test results

    由表2可知,改進(jìn)后的蝙蝠算法對(duì)于多維的單峰、多峰函數(shù)均具有較好的尋優(yōu)性能和穩(wěn)定性,與基本蝙蝠算法相比最優(yōu)精度、最差精度均得到了較大提高,平均精度至少提高了1個(gè)數(shù)量級(jí)。

    3 基于IBA-RBF的水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)預(yù)測(cè)模型

    RBF神經(jīng)網(wǎng)絡(luò)于1988年提出[21],以非線性回歸理論為基礎(chǔ),結(jié)構(gòu)簡(jiǎn)單、不受初始權(quán)值影響,可以任意精度逼近復(fù)雜非線性函數(shù)[22],更適合于水電站機(jī)組與廠房結(jié)構(gòu)間耦聯(lián)振動(dòng)的預(yù)測(cè)分析。通過調(diào)用MATLAB工具箱中的newrbe(·)函數(shù)生成一個(gè)精確的RBF網(wǎng)絡(luò),其格式如下:

    net=newrbe(input,output,P)

    (19)

    式中:input、output分別為訓(xùn)練樣本的輸入和輸出向量;P為網(wǎng)絡(luò)的擴(kuò)展參數(shù),其值的選取對(duì)網(wǎng)絡(luò)性能有較大影響。

    為提高網(wǎng)絡(luò)的預(yù)測(cè)精度,本文采用IBA優(yōu)化RBF神經(jīng)網(wǎng)絡(luò)的擴(kuò)展參數(shù),將蝙蝠個(gè)體位置編碼為RBF網(wǎng)絡(luò)的擴(kuò)展參數(shù)P,以RBF網(wǎng)絡(luò)的均方誤差MSE為適應(yīng)度函數(shù),模型計(jì)算流程如圖1所示。

    圖1 基于IBA-RBF的水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)預(yù)測(cè)模型流程圖Fig.1 Flow chart of vibration response prediction model for hydropower house based on IBA-RBF

    4 仿真實(shí)例分析

    4.1 數(shù)據(jù)來源

    采用文獻(xiàn)[23]中的水電站原型觀測(cè)數(shù)據(jù)進(jìn)行分析,利用機(jī)組振動(dòng)和尾水脈動(dòng)壓力數(shù)據(jù)預(yù)測(cè)廠房結(jié)構(gòu)的振動(dòng)響應(yīng)。以機(jī)組頂蓋、上、下機(jī)架3個(gè)部位X、Y、Z三向振動(dòng)位移及尾水脈動(dòng)壓力10個(gè)影響因素95%的雙幅值為輸入因子,以廠房結(jié)構(gòu)定子基礎(chǔ)、下機(jī)架基礎(chǔ)、樓板、風(fēng)罩、左側(cè)柱5個(gè)測(cè)點(diǎn)Z向振動(dòng)位移標(biāo)準(zhǔn)差為輸出因子。為充分檢驗(yàn)?zāi)P托阅?,選擇與文獻(xiàn)[23]相同的訓(xùn)練與測(cè)試數(shù)據(jù),在訓(xùn)練與測(cè)試前,對(duì)數(shù)據(jù)進(jìn)行歸一化處理。

    4.2 模型的創(chuàng)建與訓(xùn)練

    確定模型的輸入層神經(jīng)元數(shù)目為10,與廠房結(jié)構(gòu)振動(dòng)的10個(gè)影響因子相對(duì)應(yīng),輸出層神經(jīng)元數(shù)目為5,對(duì)應(yīng)于廠房結(jié)構(gòu)5個(gè)測(cè)點(diǎn)的振動(dòng)位移標(biāo)準(zhǔn)差。設(shè)置IBA-RBF模型參數(shù)如下:蝙蝠種群規(guī)模M為50,最大迭代次數(shù)Tmax為50,蝙蝠個(gè)體頻率范圍[fmin,fmax]=[0,2],脈沖音量范圍[Amin,Amax]=[0,1],脈沖頻率范圍[rmin,rmax]=[0,0.5],蝙蝠個(gè)體位置范圍[xmin,xmax]=[10-5,5]。

    4.3 模型驗(yàn)證與誤差分析

    將5組不同負(fù)荷下的測(cè)試數(shù)據(jù)代入已訓(xùn)練好的IBA-RBF模型中進(jìn)行廠房結(jié)構(gòu)振動(dòng)響應(yīng)仿真,預(yù)測(cè)結(jié)果如表3所示,并與文獻(xiàn)[23]中的RVM模型進(jìn)行對(duì)比見圖2,繪制各負(fù)荷下5個(gè)測(cè)點(diǎn)振動(dòng)響應(yīng)預(yù)測(cè)結(jié)果的平均相對(duì)誤差如圖3所示。

    表3 廠房結(jié)構(gòu)振動(dòng)響應(yīng)預(yù)測(cè)結(jié)果Tab.3 Prediction results of vibration response of hydropower house structure

    由表3、圖2可以看出,IBA-RBF模型預(yù)測(cè)精度較高,除以負(fù)荷50 MW工況為測(cè)試數(shù)據(jù)誤差略大外,200、350、480、600 MW負(fù)荷下模型的預(yù)測(cè)效果均優(yōu)于RVM模型的預(yù)測(cè)效果,與實(shí)測(cè)值非常接近。由表3、圖3可知,IBA-RBF模型對(duì)于廠房結(jié)構(gòu)各測(cè)點(diǎn)、各工況的預(yù)測(cè)誤差大多位于5%以下,5組負(fù)荷下的平均相對(duì)誤差為5.6741%。而文獻(xiàn)[23]所提出的相關(guān)向量機(jī)(RVM)模型預(yù)測(cè)的相對(duì)誤差大多位于8%以下,平均相對(duì)誤差為7.476 5%。與RVM模型相比,IBA-RBF模型的泛化能力有較大提升,更適用于水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)的高維非線性預(yù)測(cè)。

    5 結(jié) 語

    本文在基本蝙蝠算法的基礎(chǔ)上,對(duì)蝙蝠種群的全局尋優(yōu)、局部搜索以及脈沖音量和脈沖頻率的更新機(jī)制進(jìn)行改進(jìn),同時(shí)引入慣性權(quán)重遞減和種群刷新策略,提出改進(jìn)蝙蝠算法IBA。并利用該算法優(yōu)化RBF網(wǎng)絡(luò)的擴(kuò)展參數(shù)P,基于機(jī)組振動(dòng)和尾水脈動(dòng)的監(jiān)測(cè)數(shù)據(jù)建立水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)預(yù)測(cè)的IBA-RBF模型,主要結(jié)論如下。

    圖2 各測(cè)點(diǎn)振動(dòng)位移實(shí)測(cè)值與預(yù)測(cè)值對(duì)比Fig. 2 Comparison of measured and predicted values of vibration displacement of each measuring point

    圖3 不同工況下平均相對(duì)誤差對(duì)比Fig.3 Comparison of average relative error under different working conditions

    (1)改進(jìn)后的蝙蝠算法一定程度上彌補(bǔ)了基本蝙蝠算法易陷入局部極值的不足,與基本蝙蝠算法相比最優(yōu)精度、最差精度、平均精度均得到了較大提高,對(duì)于多維的單峰及多峰函數(shù)都具有較好的尋優(yōu)性能。

    (2)與RVM模型相比,IBA-RBF模型預(yù)測(cè)精度和泛化能力得到了顯著提高,可為水電站廠房結(jié)構(gòu)振動(dòng)響應(yīng)監(jiān)測(cè)與控制提供一種新方法,為電站的安全運(yùn)行提供一定技術(shù)保障。

    猜你喜歡
    蝙蝠廠房水電站
    耍惹石主和他的三萬平舊廠房“畫布”
    里底水電站工程
    西北水電(2022年1期)2022-04-16 05:20:06
    白鶴灘水電站有多“?!?
    金橋(2021年8期)2021-08-23 01:06:54
    讓老廠房變文創(chuàng)新地標(biāo)
    舊廠房改造之庭院微氣候創(chuàng)建
    蝙蝠
    小型水電站技術(shù)改造探討
    蝙蝠女
    蝙蝠在黑暗處如何捕食
    電子廠房多種接地連接方式的探討與實(shí)踐
    国产成人一区二区在线| 高清不卡的av网站| 色精品久久人妻99蜜桃| 国产伦理片在线播放av一区| 九九爱精品视频在线观看| 精品少妇一区二区三区视频日本电影 | 97精品久久久久久久久久精品| 久久ye,这里只有精品| 大码成人一级视频| 极品人妻少妇av视频| 成人黄色视频免费在线看| 亚洲七黄色美女视频| 在线观看免费高清a一片| 在线 av 中文字幕| 亚洲成人一二三区av| 国产成人欧美| 波多野结衣av一区二区av| 午夜免费男女啪啪视频观看| 亚洲,欧美,日韩| 一区二区三区乱码不卡18| 免费不卡黄色视频| 又大又爽又粗| 亚洲av日韩在线播放| 欧美日韩亚洲综合一区二区三区_| 捣出白浆h1v1| 欧美激情极品国产一区二区三区| 99re6热这里在线精品视频| 国产不卡av网站在线观看| 一区二区日韩欧美中文字幕| 无遮挡黄片免费观看| 日本欧美视频一区| 国产乱人偷精品视频| 女的被弄到高潮叫床怎么办| 久久久精品区二区三区| 久久精品亚洲熟妇少妇任你| 老司机深夜福利视频在线观看 | 啦啦啦中文免费视频观看日本| 国产精品人妻久久久影院| 最近最新中文字幕免费大全7| 大香蕉久久成人网| 国产av码专区亚洲av| 日本黄色日本黄色录像| 国产精品二区激情视频| 久久久欧美国产精品| 欧美在线黄色| 国产亚洲精品第一综合不卡| 伊人久久国产一区二区| 亚洲av在线观看美女高潮| 亚洲av成人不卡在线观看播放网 | 久久久国产精品麻豆| 国产成人一区二区在线| 国产精品久久久人人做人人爽| 男人添女人高潮全过程视频| 午夜免费观看性视频| 亚洲人成网站在线观看播放| 午夜免费男女啪啪视频观看| 亚洲成人国产一区在线观看 | 韩国精品一区二区三区| 在线亚洲精品国产二区图片欧美| 男人舔女人的私密视频| 晚上一个人看的免费电影| 国产爽快片一区二区三区| 18禁动态无遮挡网站| 天天躁夜夜躁狠狠躁躁| av线在线观看网站| 天天影视国产精品| 热99久久久久精品小说推荐| 如何舔出高潮| 黄色 视频免费看| 自线自在国产av| av线在线观看网站| 日韩 欧美 亚洲 中文字幕| 免费少妇av软件| tube8黄色片| 男人添女人高潮全过程视频| 精品人妻在线不人妻| 日韩 欧美 亚洲 中文字幕| 巨乳人妻的诱惑在线观看| 午夜福利一区二区在线看| 久久精品熟女亚洲av麻豆精品| 纯流量卡能插随身wifi吗| 侵犯人妻中文字幕一二三四区| 看免费av毛片| www.精华液| 日本wwww免费看| 中文欧美无线码| 90打野战视频偷拍视频| 久久精品aⅴ一区二区三区四区| 国产色婷婷99| 国产日韩欧美亚洲二区| 亚洲av在线观看美女高潮| 精品国产超薄肉色丝袜足j| av电影中文网址| 一级毛片我不卡| 91精品伊人久久大香线蕉| 国产精品秋霞免费鲁丝片| 人人妻人人澡人人看| 欧美日韩成人在线一区二区| 一级,二级,三级黄色视频| 男女边摸边吃奶| 日韩一本色道免费dvd| 国产福利在线免费观看视频| 亚洲七黄色美女视频| 精品国产乱码久久久久久男人| www.av在线官网国产| 高清视频免费观看一区二区| 成人午夜精彩视频在线观看| 欧美av亚洲av综合av国产av | 日韩,欧美,国产一区二区三区| 日韩欧美一区视频在线观看| 国产日韩欧美亚洲二区| 免费黄网站久久成人精品| 飞空精品影院首页| 婷婷色av中文字幕| 久久久久国产精品人妻一区二区| 国产片特级美女逼逼视频| 下体分泌物呈黄色| 婷婷色综合www| 国产1区2区3区精品| 毛片一级片免费看久久久久| 美女福利国产在线| av卡一久久| 在线观看免费高清a一片| 一二三四中文在线观看免费高清| 国产精品女同一区二区软件| 国产视频首页在线观看| 中文精品一卡2卡3卡4更新| 天天操日日干夜夜撸| 亚洲人成77777在线视频| 黄色一级大片看看| 国产亚洲精品第一综合不卡| 如何舔出高潮| 中文字幕av电影在线播放| 中国国产av一级| 免费久久久久久久精品成人欧美视频| 欧美国产精品一级二级三级| 国产精品亚洲av一区麻豆 | 99热国产这里只有精品6| 亚洲色图 男人天堂 中文字幕| 日韩人妻精品一区2区三区| 午夜福利网站1000一区二区三区| 日韩大片免费观看网站| 国产精品无大码| 中文字幕人妻熟女乱码| 操美女的视频在线观看| 亚洲av男天堂| 爱豆传媒免费全集在线观看| 国产精品亚洲av一区麻豆 | 亚洲国产精品成人久久小说| 欧美国产精品一级二级三级| 在线观看国产h片| 男女高潮啪啪啪动态图| 国产精品一区二区在线不卡| 综合色丁香网| 久久精品久久精品一区二区三区| 91老司机精品| 亚洲色图综合在线观看| 国产精品女同一区二区软件| 精品一区二区三区av网在线观看 | 卡戴珊不雅视频在线播放| 高清不卡的av网站| 亚洲精品国产av成人精品| 国产精品久久久久久精品电影小说| 中文乱码字字幕精品一区二区三区| a 毛片基地| 久久久精品区二区三区| 女的被弄到高潮叫床怎么办| 国产在线一区二区三区精| 人人妻,人人澡人人爽秒播 | 亚洲美女视频黄频| 少妇人妻 视频| 久久久久视频综合| 欧美国产精品一级二级三级| 久久久久久久久免费视频了| 久久 成人 亚洲| 一级毛片我不卡| 国产精品99久久99久久久不卡 | 亚洲,欧美,日韩| 日本黄色日本黄色录像| 波野结衣二区三区在线| 满18在线观看网站| 国产黄色免费在线视频| 大香蕉久久成人网| 欧美日韩福利视频一区二区| 久久久久精品久久久久真实原创| 日本av免费视频播放| 一级毛片 在线播放| 亚洲av男天堂| 国产男人的电影天堂91| 国产精品 欧美亚洲| av国产精品久久久久影院| 国产亚洲午夜精品一区二区久久| 天堂8中文在线网| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩国产mv在线观看视频| 国产精品嫩草影院av在线观看| a级毛片黄视频| 99热全是精品| 晚上一个人看的免费电影| 欧美日韩亚洲高清精品| 狠狠精品人妻久久久久久综合| 999久久久国产精品视频| 亚洲美女搞黄在线观看| 亚洲国产成人一精品久久久| 建设人人有责人人尽责人人享有的| 亚洲熟女毛片儿| 久久久久久人妻| 18禁动态无遮挡网站| 三上悠亚av全集在线观看| 欧美av亚洲av综合av国产av | 国产精品 国内视频| 亚洲精品久久久久久婷婷小说| 最黄视频免费看| e午夜精品久久久久久久| 18禁裸乳无遮挡动漫免费视频| 国产精品 欧美亚洲| 人人澡人人妻人| 国语对白做爰xxxⅹ性视频网站| 秋霞在线观看毛片| 久久久久久人妻| 天天躁狠狠躁夜夜躁狠狠躁| h视频一区二区三区| 少妇被粗大的猛进出69影院| 啦啦啦 在线观看视频| 精品酒店卫生间| 日韩人妻精品一区2区三区| bbb黄色大片| 2018国产大陆天天弄谢| av卡一久久| a 毛片基地| 一级爰片在线观看| 久久久国产精品麻豆| 两性夫妻黄色片| 男女边吃奶边做爰视频| 一区二区三区激情视频| 亚洲精品久久成人aⅴ小说| 亚洲人成网站在线观看播放| 五月开心婷婷网| 国产成人精品福利久久| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品一区三区| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人一二三区av| 亚洲成人免费av在线播放| 19禁男女啪啪无遮挡网站| 亚洲av国产av综合av卡| 免费不卡黄色视频| 国产精品三级大全| 国产精品女同一区二区软件| 美女国产高潮福利片在线看| 嫩草影视91久久| 精品少妇黑人巨大在线播放| 人成视频在线观看免费观看| 国产国语露脸激情在线看| √禁漫天堂资源中文www| 久久97久久精品| www日本在线高清视频| 18禁观看日本| 国产av一区二区精品久久| 捣出白浆h1v1| 日韩人妻精品一区2区三区| 99热国产这里只有精品6| √禁漫天堂资源中文www| 精品少妇久久久久久888优播| 一本大道久久a久久精品| 纯流量卡能插随身wifi吗| 免费看av在线观看网站| 少妇精品久久久久久久| 国产又色又爽无遮挡免| 国产av码专区亚洲av| 亚洲专区中文字幕在线 | 国产免费视频播放在线视频| 成人影院久久| 咕卡用的链子| 欧美亚洲日本最大视频资源| 丁香六月天网| 亚洲久久久国产精品| 欧美激情高清一区二区三区 | 午夜免费鲁丝| 热re99久久国产66热| 老司机亚洲免费影院| 欧美成人精品欧美一级黄| 国产色婷婷99| 欧美人与善性xxx| 又大又爽又粗| 十分钟在线观看高清视频www| 9191精品国产免费久久| 亚洲欧美精品综合一区二区三区| 久久人人97超碰香蕉20202| 1024香蕉在线观看| 满18在线观看网站| 免费不卡黄色视频| xxxhd国产人妻xxx| 国产亚洲av高清不卡| 免费少妇av软件| 国产免费现黄频在线看| 婷婷色麻豆天堂久久| 男女边吃奶边做爰视频| 一二三四在线观看免费中文在| 亚洲图色成人| 亚洲成av片中文字幕在线观看| 欧美精品一区二区大全| 精品一区二区免费观看| 欧美日韩视频高清一区二区三区二| 色婷婷久久久亚洲欧美| 亚洲精品一区蜜桃| 亚洲,欧美精品.| a级毛片在线看网站| 一级毛片 在线播放| 美女福利国产在线| 亚洲av综合色区一区| 99精品久久久久人妻精品| 男女高潮啪啪啪动态图| 国产日韩欧美视频二区| 亚洲男人天堂网一区| 两性夫妻黄色片| 91精品国产国语对白视频| 国产精品国产三级国产专区5o| 赤兔流量卡办理| 男女之事视频高清在线观看 | 久久久久精品性色| videos熟女内射| 色94色欧美一区二区| 午夜av观看不卡| 亚洲av电影在线观看一区二区三区| 色视频在线一区二区三区| av在线播放精品| 亚洲欧美激情在线| 少妇精品久久久久久久| 国产免费又黄又爽又色| 久久亚洲国产成人精品v| 国产精品二区激情视频| 久久精品国产亚洲av高清一级| 日韩制服骚丝袜av| 亚洲,一卡二卡三卡| 婷婷色综合www| 久久久久视频综合| 热re99久久国产66热| 久久久久久免费高清国产稀缺| 日韩电影二区| 在线精品无人区一区二区三| 街头女战士在线观看网站| 久久久久视频综合| 亚洲成人手机| av线在线观看网站| avwww免费| 黄频高清免费视频| 99香蕉大伊视频| 日日爽夜夜爽网站| 蜜桃国产av成人99| 赤兔流量卡办理| 91精品三级在线观看| 国产精品国产三级专区第一集| 婷婷色综合大香蕉| 国产一区二区 视频在线| 青春草国产在线视频| 日韩一区二区视频免费看| 永久免费av网站大全| 欧美日韩一级在线毛片| 久久久久精品国产欧美久久久 | 欧美日韩av久久| 国产一级毛片在线| 99re6热这里在线精品视频| 国产精品无大码| 天堂8中文在线网| 女的被弄到高潮叫床怎么办| kizo精华| videosex国产| 99热国产这里只有精品6| 亚洲色图综合在线观看| 亚洲精品国产区一区二| 深夜精品福利| 中文欧美无线码| 久久久久久久久免费视频了| 777久久人妻少妇嫩草av网站| 国产色婷婷99| 亚洲欧美精品自产自拍| 亚洲一级一片aⅴ在线观看| 国产在视频线精品| 久久久久久久大尺度免费视频| 国产伦理片在线播放av一区| 久久99一区二区三区| 18禁裸乳无遮挡动漫免费视频| 91精品伊人久久大香线蕉| 老熟女久久久| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 久久国产亚洲av麻豆专区| 久久久久久人妻| 黑丝袜美女国产一区| 亚洲国产av影院在线观看| 热re99久久国产66热| av在线老鸭窝| 成人毛片60女人毛片免费| 欧美精品av麻豆av| www日本在线高清视频| 国产免费又黄又爽又色| 午夜免费男女啪啪视频观看| 亚洲精品日韩在线中文字幕| 精品一区在线观看国产| 夜夜骑夜夜射夜夜干| 波多野结衣av一区二区av| 人体艺术视频欧美日本| 日韩制服骚丝袜av| 王馨瑶露胸无遮挡在线观看| 悠悠久久av| 成年人免费黄色播放视频| 国产精品偷伦视频观看了| 久久久精品94久久精品| 九色亚洲精品在线播放| 日本黄色日本黄色录像| 爱豆传媒免费全集在线观看| 亚洲精品国产区一区二| 大陆偷拍与自拍| 伦理电影免费视频| 国产在线一区二区三区精| 亚洲国产欧美一区二区综合| 久久韩国三级中文字幕| 亚洲精品一二三| 久久久精品国产亚洲av高清涩受| 一级a爱视频在线免费观看| 亚洲人成77777在线视频| 成人影院久久| videos熟女内射| 天堂中文最新版在线下载| 一本—道久久a久久精品蜜桃钙片| 亚洲国产最新在线播放| 在线观看人妻少妇| 满18在线观看网站| 国产精品一区二区在线不卡| 九色亚洲精品在线播放| 亚洲视频免费观看视频| 亚洲色图 男人天堂 中文字幕| 亚洲情色 制服丝袜| 女的被弄到高潮叫床怎么办| 90打野战视频偷拍视频| 成人国产av品久久久| 下体分泌物呈黄色| 男女床上黄色一级片免费看| 99精品久久久久人妻精品| 免费日韩欧美在线观看| 欧美激情 高清一区二区三区| av福利片在线| 午夜福利影视在线免费观看| 亚洲av综合色区一区| 国产精品一区二区在线观看99| 亚洲av福利一区| 麻豆乱淫一区二区| 在线精品无人区一区二区三| 亚洲精品乱久久久久久| 一二三四在线观看免费中文在| 精品国产一区二区三区四区第35| 汤姆久久久久久久影院中文字幕| 亚洲成人手机| 80岁老熟妇乱子伦牲交| 大香蕉久久成人网| 久久精品久久精品一区二区三区| 欧美日韩视频精品一区| 卡戴珊不雅视频在线播放| 免费少妇av软件| 免费观看a级毛片全部| 精品久久久精品久久久| 免费观看性生交大片5| 99久久人妻综合| 又大又爽又粗| netflix在线观看网站| 亚洲国产毛片av蜜桃av| 免费高清在线观看视频在线观看| 国产色婷婷99| 少妇被粗大猛烈的视频| 亚洲精品第二区| 亚洲中文av在线| 一本—道久久a久久精品蜜桃钙片| 日韩中文字幕欧美一区二区 | 亚洲精品久久午夜乱码| 国产淫语在线视频| 日韩精品有码人妻一区| 黄网站色视频无遮挡免费观看| 99国产综合亚洲精品| 啦啦啦中文免费视频观看日本| 午夜福利一区二区在线看| 大码成人一级视频| 18禁观看日本| 夫妻性生交免费视频一级片| 男女高潮啪啪啪动态图| 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区| 99久国产av精品国产电影| 女人久久www免费人成看片| www.熟女人妻精品国产| 国产爽快片一区二区三区| 这个男人来自地球电影免费观看 | 亚洲精品成人av观看孕妇| 国产在线免费精品| 免费av中文字幕在线| 国产日韩欧美亚洲二区| 在线观看免费日韩欧美大片| 又黄又粗又硬又大视频| 十八禁网站网址无遮挡| 午夜久久久在线观看| 精品人妻一区二区三区麻豆| 只有这里有精品99| 五月天丁香电影| 国产精品一区二区在线观看99| 久久综合国产亚洲精品| 日本欧美视频一区| 少妇人妻 视频| 最近2019中文字幕mv第一页| 国产不卡av网站在线观看| 日韩中文字幕视频在线看片| 97精品久久久久久久久久精品| 久久婷婷青草| 多毛熟女@视频| 精品亚洲成a人片在线观看| 日本vs欧美在线观看视频| 18禁动态无遮挡网站| 中文精品一卡2卡3卡4更新| 欧美黄色片欧美黄色片| 又黄又粗又硬又大视频| 国产精品一区二区在线不卡| 老司机靠b影院| xxx大片免费视频| 中文字幕人妻丝袜一区二区 | 香蕉丝袜av| 亚洲精品成人av观看孕妇| 精品一区二区三区四区五区乱码 | www.av在线官网国产| 中文乱码字字幕精品一区二区三区| 波多野结衣一区麻豆| 巨乳人妻的诱惑在线观看| 久久热在线av| 男女之事视频高清在线观看 | av在线老鸭窝| 91老司机精品| 你懂的网址亚洲精品在线观看| 一区二区三区激情视频| 少妇的丰满在线观看| 99re6热这里在线精品视频| 亚洲四区av| 国产 一区精品| 久久国产精品男人的天堂亚洲| 精品少妇久久久久久888优播| 香蕉丝袜av| 超碰97精品在线观看| 国产精品亚洲av一区麻豆 | 亚洲综合精品二区| 国产精品无大码| 国产又爽黄色视频| 久久久久久久精品精品| 交换朋友夫妻互换小说| www.自偷自拍.com| 一级毛片电影观看| 美女脱内裤让男人舔精品视频| 国产乱来视频区| 亚洲成国产人片在线观看| 中文字幕人妻丝袜一区二区 | 亚洲熟女精品中文字幕| 美女中出高潮动态图| 国产1区2区3区精品| 日韩av在线免费看完整版不卡| 十八禁高潮呻吟视频| 9色porny在线观看| 国产免费福利视频在线观看| 纯流量卡能插随身wifi吗| 黑人猛操日本美女一级片| 亚洲成人av在线免费| 性色av一级| 中文精品一卡2卡3卡4更新| 国产成人精品福利久久| 久热爱精品视频在线9| e午夜精品久久久久久久| 人妻一区二区av| 欧美国产精品va在线观看不卡| 国产精品成人在线| 搡老岳熟女国产| 丁香六月欧美| 又大又爽又粗| 欧美av亚洲av综合av国产av | 亚洲av日韩在线播放| 亚洲专区中文字幕在线 | 亚洲成色77777| 黄网站色视频无遮挡免费观看| av又黄又爽大尺度在线免费看| 国产成人精品福利久久| 99热全是精品| 精品午夜福利在线看| 性高湖久久久久久久久免费观看| 一边摸一边做爽爽视频免费| 最近最新中文字幕免费大全7| 久久久久精品人妻al黑| 最近2019中文字幕mv第一页| 成人国产av品久久久| 校园人妻丝袜中文字幕| 亚洲av欧美aⅴ国产| 极品少妇高潮喷水抽搐| 久久99热这里只频精品6学生| 欧美在线一区亚洲| 99久久人妻综合| 精品一区二区三区av网在线观看 | 夜夜骑夜夜射夜夜干| 中文字幕人妻熟女乱码| 久久天堂一区二区三区四区| 国产在视频线精品| 看非洲黑人一级黄片| 考比视频在线观看| 国产精品久久久久成人av| 久久久久精品久久久久真实原创| 欧美老熟妇乱子伦牲交| a级毛片黄视频| 日韩制服骚丝袜av| 十八禁网站网址无遮挡| 国产精品久久久久成人av| av在线观看视频网站免费| 乱人伦中国视频| 黑人欧美特级aaaaaa片|