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

    基于人工魚(yú)群算法的隨機(jī)結(jié)構(gòu)AMD控制系統(tǒng)LQR權(quán)矩陣優(yōu)化

    2016-05-20 03:53:47李森萍廣州大學(xué)工程抗震研究中心廣州510405
    振動(dòng)與沖擊 2016年8期

    王 磊, 譚 平, 李森萍(廣州大學(xué) 工程抗震研究中心,廣州 510405)

    ?

    基于人工魚(yú)群算法的隨機(jī)結(jié)構(gòu)AMD控制系統(tǒng)LQR權(quán)矩陣優(yōu)化

    王磊, 譚平, 李森萍(廣州大學(xué) 工程抗震研究中心,廣州510405)

    摘要:針對(duì)隨機(jī)結(jié)構(gòu)在平穩(wěn)隨機(jī)激勵(lì)下線性二次最優(yōu)控制的權(quán)參數(shù)選取問(wèn)題,提出了基于人工魚(yú)群算法的隨機(jī)結(jié)構(gòu)AMD控制系統(tǒng)權(quán)矩陣優(yōu)化設(shè)計(jì)方法。該方法以結(jié)構(gòu)隨機(jī)響應(yīng)和控制力均方值的加權(quán)組合為目標(biāo)函數(shù),考慮了結(jié)構(gòu)和外激勵(lì)的隨機(jī)性對(duì)于控制效果的影響。其優(yōu)化結(jié)果不僅保證了控制器的控制效果,而且降低了控制效果對(duì)于隨機(jī)參數(shù)的敏感性,增強(qiáng)了控制器的魯棒性。最后結(jié)合數(shù)值算例,驗(yàn)證了所述方法的有效性和正確性。

    關(guān)鍵詞:隨機(jī)結(jié)構(gòu);權(quán)矩陣;主動(dòng)質(zhì)量阻尼器;魚(yú)群算法

    線性二次型經(jīng)典最優(yōu)控制算法(Linear Quadratic Optimal Control,LQR)因其性能指標(biāo)物理意義明確,最優(yōu)控制解可以解析表達(dá),易于實(shí)現(xiàn)閉環(huán)反饋控制,而被廣泛運(yùn)用于結(jié)構(gòu)振動(dòng)控制之中[1-4]。LQR控制器的控制效果直接取決于加權(quán)矩陣的選取,如果加權(quán)矩陣選取不當(dāng),則可能造成控制裝置不能發(fā)揮最優(yōu)性能,造成資源浪費(fèi),甚至對(duì)結(jié)構(gòu)響應(yīng)產(chǎn)生不利的影響。對(duì)于權(quán)矩陣的選取較常用的為試探法,即選擇不同的權(quán)矩陣比較計(jì)算選取最優(yōu)值。此種方法一般僅能得到局部最優(yōu)解。鑒于此,發(fā)展了多種優(yōu)化權(quán)矩陣的方法。郭一峰等[5-6]基于遺傳算法采用不同的目標(biāo)函數(shù)對(duì)權(quán)矩陣進(jìn)行了優(yōu)化,其優(yōu)化后的控制效果都滿足設(shè)計(jì)和結(jié)構(gòu)性能的要求;張葛祥等[7]利用參數(shù)滿意度對(duì)權(quán)參數(shù)進(jìn)行了優(yōu)化,綜合考慮了控制性能和能量間的關(guān)系;祁皚等[8]提出了一種新的權(quán)矩陣選擇方法,并詳細(xì)分析了權(quán)矩陣參數(shù)對(duì)控制效果的影響;楊剛等[9]提出了僅包含一個(gè)可調(diào)參數(shù)的新目標(biāo)函數(shù),并利用非線性優(yōu)化理論對(duì)該目標(biāo)函數(shù)進(jìn)行了優(yōu)化分析。在以往的研究中,對(duì)于權(quán)矩陣的優(yōu)化大多集中于確定性結(jié)構(gòu)在特定激勵(lì)下的分析。然而在實(shí)際環(huán)境之中,存在著大量的不確定性,僅對(duì)特定激勵(lì)下確定性結(jié)構(gòu)進(jìn)行權(quán)矩陣優(yōu)化,忽略結(jié)構(gòu)和外部環(huán)境激勵(lì)的隨機(jī)性對(duì)于控制效果的影響,其優(yōu)化結(jié)果仍具有一定的局限性,難以全面表現(xiàn)控制裝置對(duì)于結(jié)構(gòu)響應(yīng)的控制效果。

    本文基于人工魚(yú)群算法,對(duì)隨機(jī)地震下隨機(jī)結(jié)構(gòu)AMD控制系統(tǒng)的權(quán)矩陣進(jìn)行了參數(shù)優(yōu)化。其中對(duì)于結(jié)構(gòu)參數(shù)和環(huán)境激勵(lì)的隨機(jī)性造成的復(fù)合隨機(jī)問(wèn)題,首先通過(guò)結(jié)構(gòu)響應(yīng)對(duì)隨機(jī)參數(shù)靈敏度分析,合理的確定結(jié)構(gòu)隨機(jī)參數(shù),進(jìn)一步利用Gegenbauer多項(xiàng)式將含λ-PDF分布隨機(jī)變量的隨機(jī)結(jié)構(gòu)化為等效擴(kuò)階系統(tǒng),采用Kanai-Tajimi地震動(dòng)模型為地震動(dòng)輸入,計(jì)算結(jié)構(gòu)隨機(jī)響應(yīng)。最后結(jié)合數(shù)值算例驗(yàn)證了所提方法的有效性及所優(yōu)化參數(shù)的魯棒性。

    1隨機(jī)結(jié)構(gòu)控制系統(tǒng)正交分解

    對(duì)于隨機(jī)結(jié)構(gòu)AMD振動(dòng)控制系統(tǒng),其動(dòng)力平衡方程為

    Ks(θ)Xs=-MsIf(t)+HfTMD

    (1)

    (2)

    fTMD=kT(xT-xn)+cT(xT-xn)-u(t,θ)

    (3)

    控制力u(t)可采用LQR控制算法進(jìn)行設(shè)計(jì),以層間位移為控制目標(biāo)。首先定義二次性能泛函為:

    (4)

    式中:Z(t)=[d1,d2,…,dn]T,其中di為結(jié)構(gòu)的第i層層間位移,Q為半正定矩陣,R為正定矩陣。

    不失一般性,假設(shè)式(1)中結(jié)構(gòu)控制系統(tǒng)的質(zhì)量、阻尼和剛度有下列形式[12]

    (5)

    (6)

    (7)

    (8)

    把式(5)~式(8)代入式(1)利用正交多項(xiàng)式逼近理論,得到其確定性系統(tǒng)如下

    ui=0,1,…,nii=1,…,N

    (9)

    2隨機(jī)結(jié)構(gòu)響應(yīng)靈敏度

    對(duì)于式(5)~式(7)所示的結(jié)構(gòu)隨機(jī)質(zhì)量、阻尼和剛度矩陣可分解為均值分量和零均值隨機(jī)分量之和:

    (10)

    (11)

    (12)

    隨機(jī)結(jié)構(gòu)中的隨機(jī)變量也能分解為

    (13)

    則式(1)所示隨機(jī)結(jié)構(gòu)響應(yīng)亦可做同樣分解

    (14)

    (15)

    若已知隨機(jī)結(jié)構(gòu)各個(gè)參數(shù)的分布參數(shù),利用攝動(dòng)法通過(guò)式(15)求解結(jié)構(gòu)隨機(jī)響應(yīng)對(duì)于各個(gè)隨機(jī)參數(shù)的靈敏度,從而以靈敏度為依據(jù)來(lái)選取結(jié)構(gòu)隨機(jī)變量,可有效降低隨機(jī)變量選取的盲目性和不確定性。

    3基于人工魚(yú)群算法的權(quán)矩陣優(yōu)化

    人工魚(yú)群算法是李曉磊等[13]提出的一類基于動(dòng)物行為的群體智能優(yōu)化算法。該算法通過(guò)模擬魚(yú)類的覓食、聚群、追尾、隨機(jī)等行為在搜索區(qū)域中進(jìn)行尋優(yōu),是一種群智能算法。其具有并行性、簡(jiǎn)單性、能很快跳出局部極優(yōu)值、尋優(yōu)速度快、對(duì)初值不敏感等特點(diǎn)。本文利用人工魚(yú)群算法,對(duì)權(quán)矩陣進(jìn)行參數(shù)優(yōu)化,其優(yōu)化流程如圖1所示。

    圖1 人工魚(yú)群算法優(yōu)化流程圖Fig.1 Flow chart for weighting matrices using artificial fish algorithm

    結(jié)構(gòu)在地震作用下的性能很大程度上取決于結(jié)構(gòu)各層的層間位移,因此在本文中所設(shè)計(jì)AMD系統(tǒng)主要的控制對(duì)象選取為結(jié)構(gòu)層間位移,權(quán)矩陣設(shè)定為Q=10βI,R=1,其中β為待優(yōu)化參數(shù)。對(duì)于實(shí)際工程,控制系統(tǒng)的輸入能量一定是有界的,為了實(shí)現(xiàn)前述控制目標(biāo)并兼顧能量消耗,本文定義食物濃度函數(shù)(目標(biāo)函數(shù))如下

    (16)

    4仿真分析

    4.1算例模型

    三層剪切型隨機(jī)結(jié)構(gòu),其結(jié)構(gòu)參數(shù)如表1所示,其中對(duì)于結(jié)構(gòu)質(zhì)量參數(shù),由于在實(shí)際工程中其變異系數(shù)常常較小,所以被考慮為確定性參數(shù)。

    表1 結(jié)構(gòu)隨機(jī)參數(shù)

    注:質(zhì)量、剛度、阻尼單位分別為:t, kN/m, kN·s/m.

    4.2隨機(jī)參數(shù)靈敏度分析

    表1中結(jié)構(gòu)隨機(jī)參數(shù)對(duì)于結(jié)構(gòu)層間位移響應(yīng)均方值的一階靈敏度為式(13)計(jì)算結(jié)果的絕對(duì)值,對(duì)于式(13)利用攝動(dòng)法求解(攝動(dòng)步長(zhǎng)為0.01σ,σ為隨機(jī)參數(shù)的標(biāo)準(zhǔn)差)。由各隨機(jī)參數(shù)靈敏度分析可知:① 剛度對(duì)于結(jié)構(gòu)各層層間位移均方值的影響較阻尼的影響大;② 各層剛度的隨機(jī)變化都對(duì)本層的層間位移均方值有較大影響,對(duì)其他各層影響程度較小且基本相當(dāng);③ 各層阻尼的隨機(jī)變化對(duì)層間位移均方值影響規(guī)律相同:對(duì)一層影響最大,三層次之,二層最小。依據(jù)以上分析,后文在對(duì)于控制系統(tǒng)權(quán)矩陣優(yōu)化時(shí)僅考慮一、二層剛度為隨機(jī)參數(shù),其分布如圖2所示,其中一層剛度變異系數(shù)為8%,二層為7%,其概率密度函數(shù)如下所示

    (17)

    圖2 概率密度曲線Fig.2 PDF curve

    4.3權(quán)矩陣優(yōu)化

    人工魚(yú)群算法最終優(yōu)化結(jié)果為β=15.013,其優(yōu)化過(guò)程和優(yōu)化前后魚(yú)群在 [10,16,-3,0]區(qū)域內(nèi)分布分別如圖3、4所示,由圖可知魚(yú)群算法對(duì)本文優(yōu)化目標(biāo)收斂快、計(jì)算效率高,在優(yōu)化的初始,魚(yú)群的分布較為均勻,隨著優(yōu)化的進(jìn)行,魚(yú)群通過(guò)覓食、聚群、追尾、隨機(jī)等行為逐漸向最優(yōu)值靠攏,尋優(yōu)過(guò)程相對(duì)簡(jiǎn)單、快捷,特別是由于考慮了魚(yú)群的擁擠度因子,將會(huì)有少數(shù)人工魚(yú)分布在次優(yōu)值附近,所以人工魚(yú)群算法較其他優(yōu)化算法容易跳出優(yōu)化函數(shù)局部最優(yōu)值的影響從而得到函數(shù)的全局最優(yōu)值。

    圖3 優(yōu)化過(guò)程Fig.3 Optimization process

    圖4 人工魚(yú)分布Fig.4 Artificial fish distribution

    參數(shù)量值參數(shù)量值ωg/(rad·s-1)15.6mt(t)1.44ξg0.6kt/(kN·m-1)23.7S0/(cm2·s-3)15.74ct/(kN·s·m-1)0.7

    表4對(duì)比了無(wú)控、TMD控制、AMD控制結(jié)構(gòu)的層間位移均方根。由表4分析可知,相對(duì)于無(wú)控結(jié)構(gòu),主、被動(dòng)控制系統(tǒng)對(duì)于隨機(jī)結(jié)構(gòu)層間位移響應(yīng)均方根都具有一定的控制效果。但是由于參數(shù)的不確定性,特別是結(jié)構(gòu)參數(shù)隨機(jī)性的影響,TMD對(duì)于隨機(jī)結(jié)構(gòu)響應(yīng)的控制效果不明顯。但是由于主動(dòng)控制力的引入,AMD控制效果較TMD有較大的提高,AMD能更有效地降低結(jié)構(gòu)各層層間位移響應(yīng),減小結(jié)構(gòu)在外激勵(lì)下的損傷。因此,采用本文所述魚(yú)群算法優(yōu)化的權(quán)矩陣,在綜合考慮了控制能量的前提下,取得了較好的控制效果,增大了控制系統(tǒng)對(duì)于外界不確定因素的魯棒性,降低了控制效果對(duì)于隨機(jī)參數(shù)的敏感性,增強(qiáng)了結(jié)構(gòu)的安全性與魯棒性。

    表3 人工魚(yú)群算法參數(shù)

    表4 結(jié)構(gòu)響應(yīng)均方根

    注:控制效果中括號(hào)內(nèi)外分別為AMD和TMD的控制效果。

    5結(jié)論

    本文在充分考慮了外激勵(lì)和結(jié)構(gòu)參數(shù)隨機(jī)性的影響下,利用魚(yú)群算法對(duì)隨機(jī)結(jié)構(gòu)AMD控制系統(tǒng)的LQR權(quán)矩陣進(jìn)行了優(yōu)化分析,得到了如下一些結(jié)論:

    (1) 實(shí)際工程之中參數(shù)的不確定性是不可避免的,在合理設(shè)計(jì)控制系統(tǒng)時(shí)應(yīng)充分考慮結(jié)構(gòu)參數(shù)和外激勵(lì)的隨機(jī)性。本文提出的基于人工魚(yú)群算法的LQR權(quán)矩陣優(yōu)化設(shè)計(jì)方法,不依賴于特定的外激勵(lì)時(shí)程,能有效減小結(jié)構(gòu)層間位移響應(yīng)的均方根,降低AMD系統(tǒng)控制效果對(duì)于隨機(jī)參數(shù)的敏感性,提高AMD控制系統(tǒng)的魯棒性,具有重要的理論意義和工程應(yīng)用前景;

    (2) 人工魚(yú)群算法對(duì)LQR權(quán)矩陣優(yōu)化過(guò)程表明其具有簡(jiǎn)單高效、不易陷入局部最優(yōu)、對(duì)初值不敏感的特點(diǎn);

    (3) 對(duì)于結(jié)構(gòu)隨機(jī)參數(shù)敏感性分析表明,相較于阻尼,剛度的隨機(jī)變化對(duì)結(jié)構(gòu)響應(yīng)的影響更大;各層剛度隨機(jī)變化對(duì)各層層間位移均方根的影響規(guī)律為:對(duì)本層影響最大,其他各層影響次之且程度基本相當(dāng);對(duì)于各層阻尼的隨機(jī)變化,其對(duì)于結(jié)構(gòu)層間位移影響規(guī)律基本相同:對(duì)一層影響最大,三層次之,二層最小;

    (4) 仿真分析結(jié)果表明,權(quán)矩陣優(yōu)化后AMD對(duì)于隨機(jī)結(jié)構(gòu)響應(yīng)的控制效果優(yōu)于TMD,在綜合考慮能量消耗的條件下,AMD能有效降低結(jié)構(gòu)的層間位移響應(yīng),增大結(jié)構(gòu)的安全性,使結(jié)構(gòu)性能滿足設(shè)計(jì)要求。

    參 考 文 獻(xiàn)

    [ 1 ] Ghaffarzadeh H, Younespour A.Active tendons control of structures using block pulse functions[J].Structural Control and Health Monitoring, 2014,21:1453-1464.

    [ 2 ] Ambrosio P, Cazzulani G, Resta F, et al.An optimal vibration control logic for minimizing fatigue damage in flexible structures[J].Journal of Sound and Vibration,2014, 333:1269-1280.

    [ 3 ] Gagge F L A,Da Conceic?o S M, Vasques C H,et al.Experimental identification and control of a cantilever beam using ERA/OKID with a LQR Controller[J].Journal of Control, Automation and Electrical Systems, 2014,25:161-173.

    [ 4 ] 符川,屈鐵軍,孫世國(guó).動(dòng)調(diào)頻液柱阻尼器基于遺傳算法的LQR控制優(yōu)化設(shè)計(jì)[J].振動(dòng)與沖擊,2015,34(2):210-214.

    FU Chuan,QU Tie-jun,SUN Shi-guo.Optimal design of ATLCD with LQR control based on genetic algorithm[J].Journal of Vibration and Shock, 2015, 34(2):210-214.

    [ 5 ] 郭一峰,徐趙東,涂青,等.基于遺傳算法的LQR算法中權(quán)矩陣的優(yōu)化分析[J].振動(dòng)與沖擊,2010,29(11):217-220.

    GUO Yi-feng, XU Zhao-dong, TU Qing,et al.Optimal analysis for weight matrices in LQR algorithm based on genetic algorithm[J].Journal of Vibration and Shock, 2010,29(11):217-220.

    [ 6 ] Jiang Bo, Wei Xin-jiang, Guo Yan-ying.Linear quadratic optimal control in active control of structural vibration systems[C]//2010 Chinese Control and Decision Conference.Xuzhou,China,2010.

    [ 7 ] 張葛祥,金煒東,胡來(lái)招.多變量系統(tǒng)控制器的參數(shù)滿意優(yōu)化設(shè)計(jì)[J].控制理論與應(yīng)用,2004,21(3):362-366.

    ZHANG Ge-xiang, JIN Wei-dong, HU Lai-zhao.Parameters satisfactory optimization of controller in multivariable control system[J].Control Theory and Applications, 2004,21(3):362-366.

    [ 8 ] 祁皚,李惠.主動(dòng)質(zhì)量控制系統(tǒng)的權(quán)矩陣選擇及其參數(shù)的敏感性分析[J].土木工程學(xué)報(bào),2003,36(11):93-98.

    QI Ai, LI Hui.Selection of weight matrix and sensitivity analysis of parameters in active mass damper system[J].China Civil Engineering Journal, 2003,36(11):93-98.

    [ 9 ] 楊剛,姚華.實(shí)用航空發(fā)動(dòng)機(jī)LQR權(quán)陣選取方法[J].南京航空航天大學(xué)學(xué)報(bào),2006,38(4):403-407.

    YANG Gang, YAO Hua.Choosing method for aeroengine LQR weighting[J].Journal of Nanjing University of Aeronautics and Astronautics, 2006,38(4):403-407.

    [10] Li Jie, Chen Jian-bing, Fan Wen-liang.The equivalent extreme-value event and evaluation of the structural system reliability[J].Structural Safety, 2007, 29: 112-131.

    [11] Spencer Jr B F, Dyke S J,Deoskar H S.Benchmark problems in structural control: Part I—active mass driven system[J].Earthquake Engineering and Structural Dynamics, 27(11), 1998, 1127-1139.

    [12] Wu Cun-li, Ma Xiao-ping, Fang Tong.A complementary note Gegenbauer polynomial approximation for random response problem of stochastic structure[J].Probabilistic Engineering Mechanics, 2006, 21:410-419.

    [13] 李曉磊,邵之江,錢(qián)積新.一種基于動(dòng)物自治體的尋優(yōu)模式:魚(yú)群算法[J].系統(tǒng)工程理論與實(shí)踐,2002,22(11):32-38.

    LI Xiao-lei, SHAO Zhi-jiang, QIAN Ji-xin.An optimizing method based on autonomous animats: fish-swarm algorithm[J].Systems Engineering—Theory and Practice,2002,22(11):32-38.

    [14] Soong T T, Dargush G F, Passive energy dissipation systems in structural engineering[M].Chichester :Wiley, 1997.

    Optimal analysis of weight matrices of LQR algorithm for stochastic structure-AMD system based on artificial fish algorithm

    WANGLei,TANPing,LISen-ping(Earthquake Engineering Research and Test Center, Guangzhou University, Guangzhou 510405, China)

    Abstract:In allusion to the choice of the weight matrices Q and R in LQR algorithm in the stochastic structure-AMD system under stationary random excitation, the optimization method of the weight matrices is proposed based on the artificial fish algorithm.The method resolves the weight matrices optimization via the objective function with the weighted root mean square of the stochastic structural responses and active control force, taking into account the effect of the randomness in the structure and excitation on the control effect.The effectiveness and validity of the procedures of the optimal designs of the weight matrices are proven by a numerical example.Numerical results preliminarily show that the optimized weight matrices can effectively suppress undesirable vibration, reduce the sensitivity of the control effectiveness to the random parameters and improve the robustness of AMD.

    Key words:stochastic structure; weight matrix; active mass damper(AMD); artificial fish algorithm

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

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

    DOI:10.13465/j.cnki.jvs.2016.08.024

    通信作者譚平 男,博士,研究員,博士生導(dǎo)師,1973年生

    收稿日期:2014-11-06修改稿收到日期:2015-04-21

    基金項(xiàng)目:國(guó)家基金項(xiàng)目(91315301-07;51408142);國(guó)家教育部長(zhǎng)江學(xué)者創(chuàng)新團(tuán)隊(duì)項(xiàng)目(IRT13057)

    第一作者 王磊 男,博士生,1987年生

    h视频一区二区三区| 人人妻人人看人人澡| 久久久精品免费免费高清| av福利片在线观看| 国产伦精品一区二区三区视频9| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 中文字幕亚洲精品专区| 免费大片黄手机在线观看| 亚洲精品日韩在线中文字幕| 黑人猛操日本美女一级片| 欧美成人一区二区免费高清观看| 国产精品久久久久久久久免| 国产在线视频一区二区| 一区在线观看完整版| 91久久精品国产一区二区三区| 亚洲精品国产av成人精品| 久久久久久久国产电影| 黄色一级大片看看| 欧美日韩精品成人综合77777| 麻豆国产97在线/欧美| 成人无遮挡网站| 国产v大片淫在线免费观看| 欧美变态另类bdsm刘玥| 午夜日本视频在线| 激情 狠狠 欧美| 大香蕉97超碰在线| 国产免费又黄又爽又色| 午夜福利网站1000一区二区三区| 国产国拍精品亚洲av在线观看| 大香蕉久久网| 99久久精品热视频| 一区二区三区免费毛片| 观看免费一级毛片| 国产爱豆传媒在线观看| 国产精品三级大全| a级毛色黄片| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 欧美国产精品一级二级三级 | 熟女人妻精品中文字幕| 赤兔流量卡办理| 久久久久久九九精品二区国产| 中文精品一卡2卡3卡4更新| 中文字幕精品免费在线观看视频 | 精品国产乱码久久久久久小说| 18禁在线播放成人免费| 午夜激情福利司机影院| 中国美白少妇内射xxxbb| 这个男人来自地球电影免费观看 | 成人毛片a级毛片在线播放| 国产黄色免费在线视频| 日本欧美国产在线视频| 3wmmmm亚洲av在线观看| 蜜桃在线观看..| 国产在视频线精品| 黄片无遮挡物在线观看| 国产探花极品一区二区| 久久精品国产亚洲网站| 97在线视频观看| 欧美日韩在线观看h| 伊人久久精品亚洲午夜| av黄色大香蕉| 亚洲国产最新在线播放| 国产精品欧美亚洲77777| 午夜福利在线观看免费完整高清在| av专区在线播放| 青春草国产在线视频| 日本av手机在线免费观看| 亚洲欧美成人综合另类久久久| 久久鲁丝午夜福利片| 精品人妻一区二区三区麻豆| 在现免费观看毛片| 2018国产大陆天天弄谢| 91午夜精品亚洲一区二区三区| 看免费成人av毛片| 国产精品麻豆人妻色哟哟久久| 国产精品99久久99久久久不卡 | 新久久久久国产一级毛片| 六月丁香七月| 国产精品欧美亚洲77777| 精品人妻视频免费看| 欧美日韩综合久久久久久| 日韩电影二区| 国产精品久久久久久精品电影小说 | 黑丝袜美女国产一区| 免费黄频网站在线观看国产| 久久久久精品性色| 美女福利国产在线 | 久久人人爽av亚洲精品天堂 | 亚洲va在线va天堂va国产| 97在线视频观看| 菩萨蛮人人尽说江南好唐韦庄| 大话2 男鬼变身卡| 看免费成人av毛片| 夜夜爽夜夜爽视频| 久久久色成人| 噜噜噜噜噜久久久久久91| 国产伦精品一区二区三区视频9| 免费大片黄手机在线观看| 亚洲国产精品成人久久小说| 精品久久国产蜜桃| 欧美日韩精品成人综合77777| .国产精品久久| a级一级毛片免费在线观看| 久久精品国产a三级三级三级| 国产爽快片一区二区三区| 日本一二三区视频观看| 一级片'在线观看视频| 久久99蜜桃精品久久| 国产av精品麻豆| 一个人看视频在线观看www免费| 高清在线视频一区二区三区| av女优亚洲男人天堂| av在线播放精品| 五月玫瑰六月丁香| www.色视频.com| 国产黄频视频在线观看| 久久久久久久久久人人人人人人| 精华霜和精华液先用哪个| 夫妻性生交免费视频一级片| 亚洲无线观看免费| videos熟女内射| 美女高潮的动态| 三级国产精品片| 国产成人免费无遮挡视频| 97超视频在线观看视频| 在线观看一区二区三区| 亚洲精品亚洲一区二区| 成人免费观看视频高清| 99久久精品一区二区三区| 一级毛片aaaaaa免费看小| 看十八女毛片水多多多| 欧美三级亚洲精品| 大话2 男鬼变身卡| 国产精品成人在线| 99久久精品热视频| 多毛熟女@视频| 欧美性感艳星| av女优亚洲男人天堂| 国产久久久一区二区三区| 国产又色又爽无遮挡免| 人妻制服诱惑在线中文字幕| 1000部很黄的大片| 国产精品福利在线免费观看| 99久久中文字幕三级久久日本| 欧美国产精品一级二级三级 | a级毛片免费高清观看在线播放| 国产成人免费无遮挡视频| 久久久久国产精品人妻一区二区| 国产亚洲av片在线观看秒播厂| 久久亚洲国产成人精品v| 肉色欧美久久久久久久蜜桃| 只有这里有精品99| 99久久中文字幕三级久久日本| 久久影院123| 亚洲综合精品二区| 欧美丝袜亚洲另类| 美女高潮的动态| 美女中出高潮动态图| 国产成人freesex在线| 亚洲欧美一区二区三区黑人 | 午夜免费男女啪啪视频观看| 亚洲成色77777| 亚洲精品aⅴ在线观看| 国产成人免费观看mmmm| 国产高清国产精品国产三级 | 高清午夜精品一区二区三区| 日韩一本色道免费dvd| 精华霜和精华液先用哪个| 精品久久久久久久久av| 91在线精品国自产拍蜜月| 99久久精品国产国产毛片| 国产精品蜜桃在线观看| 久热这里只有精品99| 亚洲久久久国产精品| 最近中文字幕2019免费版| 亚洲欧美日韩无卡精品| 国产精品人妻久久久影院| 国产综合精华液| 婷婷色综合www| 777米奇影视久久| 精品一区在线观看国产| 国产黄片视频在线免费观看| 下体分泌物呈黄色| 新久久久久国产一级毛片| 免费看光身美女| 国产永久视频网站| 国产精品av视频在线免费观看| 国产精品国产三级国产av玫瑰| 美女xxoo啪啪120秒动态图| 国产老妇伦熟女老妇高清| 久久99热这里只频精品6学生| 菩萨蛮人人尽说江南好唐韦庄| 久久综合国产亚洲精品| 国产精品一区www在线观看| 国产深夜福利视频在线观看| 日韩在线高清观看一区二区三区| 插逼视频在线观看| 只有这里有精品99| 大码成人一级视频| 久久6这里有精品| 51国产日韩欧美| 97在线人人人人妻| a级一级毛片免费在线观看| 国产av码专区亚洲av| 永久网站在线| 亚洲成人中文字幕在线播放| 国产精品免费大片| 日日摸夜夜添夜夜爱| 国产 精品1| 久久精品国产a三级三级三级| 中国美白少妇内射xxxbb| 美女脱内裤让男人舔精品视频| 亚洲精品亚洲一区二区| 老师上课跳d突然被开到最大视频| 亚洲av成人精品一区久久| 一个人看视频在线观看www免费| 制服丝袜香蕉在线| 男女下面进入的视频免费午夜| 99九九线精品视频在线观看视频| 亚洲欧美日韩另类电影网站 | 大香蕉久久网| 久久99精品国语久久久| 各种免费的搞黄视频| 天天躁日日操中文字幕| 国产无遮挡羞羞视频在线观看| 夜夜爽夜夜爽视频| 成人亚洲欧美一区二区av| 欧美日韩视频精品一区| 精品一区在线观看国产| 国语对白做爰xxxⅹ性视频网站| 狠狠精品人妻久久久久久综合| 国产探花极品一区二区| 一级毛片黄色毛片免费观看视频| 男人舔奶头视频| xxx大片免费视频| 欧美另类一区| 一本久久精品| 久久久久久久大尺度免费视频| 久热这里只有精品99| 大话2 男鬼变身卡| 日本vs欧美在线观看视频 | 日本vs欧美在线观看视频 | 激情 狠狠 欧美| 51国产日韩欧美| 国产乱人偷精品视频| 一个人看视频在线观看www免费| 日产精品乱码卡一卡2卡三| 黑人猛操日本美女一级片| 色婷婷久久久亚洲欧美| 高清在线视频一区二区三区| 亚洲,欧美,日韩| 亚洲精品乱码久久久久久按摩| 大片电影免费在线观看免费| 激情 狠狠 欧美| 久久久久国产精品人妻一区二区| av专区在线播放| 亚洲av欧美aⅴ国产| 99久久中文字幕三级久久日本| 卡戴珊不雅视频在线播放| 国产无遮挡羞羞视频在线观看| 男人添女人高潮全过程视频| 午夜福利在线观看免费完整高清在| 亚洲av.av天堂| 18禁在线无遮挡免费观看视频| 国产视频内射| 日韩在线高清观看一区二区三区| 亚洲av欧美aⅴ国产| 亚洲内射少妇av| 国产片特级美女逼逼视频| 亚洲图色成人| 亚洲中文av在线| 成人黄色视频免费在线看| 成人18禁高潮啪啪吃奶动态图 | 美女中出高潮动态图| 国产免费一级a男人的天堂| 一边亲一边摸免费视频| 美女cb高潮喷水在线观看| 国产 一区 欧美 日韩| 黑丝袜美女国产一区| 最新中文字幕久久久久| 女人十人毛片免费观看3o分钟| 韩国av在线不卡| 国产成人一区二区在线| 国产精品偷伦视频观看了| 夫妻午夜视频| 亚洲精品第二区| 亚洲无线观看免费| 日产精品乱码卡一卡2卡三| 男女国产视频网站| 亚洲国产精品成人久久小说| 男人添女人高潮全过程视频| 身体一侧抽搐| 十分钟在线观看高清视频www | 国产一区二区三区av在线| 欧美一区二区亚洲| 18禁动态无遮挡网站| 人人妻人人澡人人爽人人夜夜| 少妇人妻久久综合中文| 国模一区二区三区四区视频| 成人综合一区亚洲| 在线看a的网站| 最近的中文字幕免费完整| 欧美97在线视频| 交换朋友夫妻互换小说| 国产精品三级大全| 五月天丁香电影| 99久久精品一区二区三区| 国产永久视频网站| 成人特级av手机在线观看| 天天躁日日操中文字幕| 最近中文字幕高清免费大全6| 亚洲一区二区三区欧美精品| 老师上课跳d突然被开到最大视频| 97热精品久久久久久| 99热这里只有是精品在线观看| 国产欧美日韩一区二区三区在线 | 国产成人精品一,二区| 久久国产精品男人的天堂亚洲 | 日韩,欧美,国产一区二区三区| 亚洲精品久久久久久婷婷小说| 亚洲精品乱久久久久久| 99久久综合免费| 国产日韩欧美亚洲二区| 欧美区成人在线视频| 啦啦啦在线观看免费高清www| 午夜激情福利司机影院| 蜜桃在线观看..| 欧美精品一区二区免费开放| 成人影院久久| 狂野欧美白嫩少妇大欣赏| 久久毛片免费看一区二区三区| 午夜日本视频在线| 国产精品精品国产色婷婷| 麻豆精品久久久久久蜜桃| 国产又色又爽无遮挡免| 日韩av在线免费看完整版不卡| 日韩强制内射视频| 午夜福利影视在线免费观看| 精品一区二区三卡| 亚洲av中文av极速乱| 欧美日韩视频高清一区二区三区二| 在线亚洲精品国产二区图片欧美 | 久久影院123| 国产视频内射| 麻豆成人av视频| 深夜a级毛片| 亚洲高清免费不卡视频| 观看免费一级毛片| 亚洲第一区二区三区不卡| 色吧在线观看| 国产精品嫩草影院av在线观看| 精品一区二区三卡| 国产永久视频网站| 舔av片在线| 亚洲国产成人一精品久久久| 亚洲最大成人中文| 黄片无遮挡物在线观看| 国产高清有码在线观看视频| 又爽又黄a免费视频| 1000部很黄的大片| 中文精品一卡2卡3卡4更新| a 毛片基地| 一级毛片aaaaaa免费看小| 久久久午夜欧美精品| 久久人人爽av亚洲精品天堂 | 国产永久视频网站| 国产视频首页在线观看| 国产无遮挡羞羞视频在线观看| 人妻制服诱惑在线中文字幕| 精品国产三级普通话版| 能在线免费看毛片的网站| 中文在线观看免费www的网站| 亚洲欧美精品自产自拍| 欧美少妇被猛烈插入视频| 在线亚洲精品国产二区图片欧美 | 男人和女人高潮做爰伦理| 男女无遮挡免费网站观看| 日日啪夜夜撸| 视频中文字幕在线观看| 特大巨黑吊av在线直播| 99精国产麻豆久久婷婷| 亚洲精品色激情综合| 少妇熟女欧美另类| 伊人久久精品亚洲午夜| 日韩制服骚丝袜av| 在现免费观看毛片| 精品视频人人做人人爽| 亚洲熟女精品中文字幕| 午夜老司机福利剧场| 大又大粗又爽又黄少妇毛片口| 一区二区三区精品91| 好男人视频免费观看在线| 在线观看免费视频网站a站| 1000部很黄的大片| 在线观看美女被高潮喷水网站| av线在线观看网站| 亚洲av欧美aⅴ国产| 久久久精品94久久精品| 纯流量卡能插随身wifi吗| 纵有疾风起免费观看全集完整版| 免费大片18禁| 久久久久久人妻| 免费观看av网站的网址| 男女边吃奶边做爰视频| 欧美3d第一页| 大陆偷拍与自拍| 一本色道久久久久久精品综合| 国产 精品1| 亚洲最大成人中文| 亚州av有码| 免费看av在线观看网站| 亚洲成色77777| 久久久国产一区二区| 亚洲av不卡在线观看| 久久人妻熟女aⅴ| 人人妻人人添人人爽欧美一区卜 | 伊人久久精品亚洲午夜| 高清日韩中文字幕在线| 国产无遮挡羞羞视频在线观看| 国产av码专区亚洲av| 激情五月婷婷亚洲| 老熟女久久久| 午夜日本视频在线| 亚洲色图综合在线观看| 精品久久久噜噜| 欧美精品一区二区免费开放| 一本一本综合久久| 赤兔流量卡办理| 免费黄色在线免费观看| 中文资源天堂在线| 日日撸夜夜添| a级毛色黄片| 精品国产露脸久久av麻豆| 美女xxoo啪啪120秒动态图| 丝袜喷水一区| 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| 另类亚洲欧美激情| 久久99热这里只有精品18| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| av国产精品久久久久影院| 久久婷婷青草| 日韩制服骚丝袜av| a级毛色黄片| 欧美精品人与动牲交sv欧美| 91精品伊人久久大香线蕉| 99国产精品免费福利视频| 我的女老师完整版在线观看| 久久人人爽人人片av| 亚洲欧美日韩卡通动漫| 亚洲第一av免费看| 免费观看a级毛片全部| 国产片特级美女逼逼视频| 伊人久久精品亚洲午夜| 国产国拍精品亚洲av在线观看| 国产v大片淫在线免费观看| 欧美精品人与动牲交sv欧美| 国精品久久久久久国模美| 午夜免费男女啪啪视频观看| 色婷婷av一区二区三区视频| 成年免费大片在线观看| 国内少妇人妻偷人精品xxx网站| 美女福利国产在线 | 九草在线视频观看| 亚洲欧洲日产国产| 又爽又黄a免费视频| 少妇人妻一区二区三区视频| a 毛片基地| 色网站视频免费| 亚洲色图综合在线观看| 秋霞伦理黄片| 99热网站在线观看| 欧美极品一区二区三区四区| 久久久久久九九精品二区国产| 亚洲综合精品二区| 伦理电影大哥的女人| 噜噜噜噜噜久久久久久91| 黄色欧美视频在线观看| 久久99热这里只有精品18| 天堂俺去俺来也www色官网| 国产女主播在线喷水免费视频网站| 一区二区三区四区激情视频| 亚洲国产精品999| 欧美精品国产亚洲| 久久精品国产鲁丝片午夜精品| 久久99蜜桃精品久久| 国产精品99久久99久久久不卡 | 日产精品乱码卡一卡2卡三| 色婷婷久久久亚洲欧美| av天堂中文字幕网| 亚洲,一卡二卡三卡| 97超碰精品成人国产| 有码 亚洲区| 精品亚洲成国产av| 日韩中字成人| 在线免费十八禁| .国产精品久久| av国产精品久久久久影院| 亚洲,欧美,日韩| 国产一区亚洲一区在线观看| 99久久人妻综合| 亚洲av成人精品一二三区| 干丝袜人妻中文字幕| 亚洲精品乱码久久久久久按摩| 亚洲国产色片| 99精国产麻豆久久婷婷| 一二三四中文在线观看免费高清| 精品一品国产午夜福利视频| 99久久精品热视频| 亚洲精品456在线播放app| 91aial.com中文字幕在线观看| 99热全是精品| 日韩av不卡免费在线播放| 欧美成人一区二区免费高清观看| 免费看光身美女| 噜噜噜噜噜久久久久久91| 国产av一区二区精品久久 | 久久久精品94久久精品| 久久久久久久亚洲中文字幕| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 久久青草综合色| 蜜桃在线观看..| a级毛片免费高清观看在线播放| 久久久久久久亚洲中文字幕| 美女主播在线视频| 三级国产精品片| 亚洲av中文av极速乱| 男人爽女人下面视频在线观看| 国产精品久久久久久久电影| 天天躁日日操中文字幕| 午夜视频国产福利| 欧美精品国产亚洲| 联通29元200g的流量卡| 777米奇影视久久| 成人美女网站在线观看视频| 亚洲自偷自拍三级| 久久久久久久亚洲中文字幕| 各种免费的搞黄视频| 精品人妻一区二区三区麻豆| 一本—道久久a久久精品蜜桃钙片| 久久久久久久亚洲中文字幕| 亚洲欧美成人综合另类久久久| 亚洲精品456在线播放app| 深夜a级毛片| 男人狂女人下面高潮的视频| 在线观看免费日韩欧美大片 | 精品人妻偷拍中文字幕| 我要看黄色一级片免费的| 男女免费视频国产| 国产美女午夜福利| 亚洲av欧美aⅴ国产| 国产 精品1| 日产精品乱码卡一卡2卡三| 欧美精品亚洲一区二区| 久久精品人妻少妇| 国产中年淑女户外野战色| 汤姆久久久久久久影院中文字幕| 成人无遮挡网站| 亚洲经典国产精华液单| 欧美日韩国产mv在线观看视频 | 大片免费播放器 马上看| 成人亚洲欧美一区二区av| 最近手机中文字幕大全| 日本黄大片高清| 久久久色成人| 成人毛片a级毛片在线播放| 少妇 在线观看| 日产精品乱码卡一卡2卡三| 亚洲欧美日韩另类电影网站 | 精品一区在线观看国产| 91狼人影院| 激情五月婷婷亚洲| 国产精品不卡视频一区二区| 丰满迷人的少妇在线观看| 亚洲性久久影院| 亚洲欧美日韩卡通动漫| 夫妻性生交免费视频一级片| 久久国内精品自在自线图片| .国产精品久久| 久久精品久久精品一区二区三区| 看免费成人av毛片| 在线观看免费高清a一片| 人体艺术视频欧美日本| 男人和女人高潮做爰伦理| 国产永久视频网站| 国产精品国产三级国产专区5o| 黄片无遮挡物在线观看| 欧美日韩亚洲高清精品| 菩萨蛮人人尽说江南好唐韦庄| 夜夜爽夜夜爽视频| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久久久| 精品少妇久久久久久888优播| av福利片在线观看| 毛片女人毛片| 亚洲精品乱码久久久久久按摩| 日本黄色片子视频| av国产免费在线观看| av在线app专区| 国产一区亚洲一区在线观看| 联通29元200g的流量卡| 国产精品爽爽va在线观看网站| 人人妻人人爽人人添夜夜欢视频 | 菩萨蛮人人尽说江南好唐韦庄| tube8黄色片| 免费久久久久久久精品成人欧美视频 | 在线观看av片永久免费下载| 少妇人妻精品综合一区二区| 成人黄色视频免费在线看| 亚洲国产精品国产精品|