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

    自適應(yīng)粒子群算法在成橋索力優(yōu)化中的應(yīng)用

    2022-02-22 05:04:52張宇航王林柯楊雅勛高宏偉
    公路交通科技 2022年1期
    關(guān)鍵詞:系桿索力吊桿

    蔡 儆,張宇航,王林柯,楊雅勛,高宏偉

    (1.陜西科技大學(xué),陜西 西安 710021;2.長安大學(xué) 公路學(xué)院,陜西 西安 710064;3.西安市軌道交通集團(tuán)有限公司,陜西 西安 710016)

    0 引言

    隨著中國公路建設(shè)的發(fā)展,系桿拱橋以其較大的整體剛度、較好的力學(xué)性能、便捷的施工方法、較低的造價等特點被國內(nèi)外廣泛應(yīng)用。系桿拱橋是一種基于拱肋、吊桿、系桿的三元結(jié)構(gòu),其外部靜定內(nèi)部超靜定,系桿拱橋的設(shè)計過程是以合理成橋狀態(tài)為目標(biāo),即在恒載作用下系桿、拱肋受力合理、主橋線形平順。在拱橋結(jié)構(gòu)已經(jīng)確定的情況下,通過改變吊桿的索力可以顯著影響結(jié)構(gòu)受力。因此為了使結(jié)構(gòu)達(dá)到目標(biāo)狀態(tài)就需要吊桿具有合適的成橋索力,系桿拱橋的優(yōu)化問題便轉(zhuǎn)化成吊桿索力的優(yōu)化問題。

    戴杰[1]等人對當(dāng)下斜拉橋索力求解方法進(jìn)行了綜述,主要方法有:指定結(jié)構(gòu)狀態(tài)的優(yōu)化方法、數(shù)學(xué)優(yōu)化法、分布優(yōu)化法等。而目前系桿拱橋吊桿的索力求解方法主要是沿用了一些斜拉橋的索力求解方法,Leonhardt率先提出了剛性支撐連續(xù)梁法并最早應(yīng)用于美國P-K橋的分析中;肖汝城[2-3]結(jié)合影響矩陣推導(dǎo)了彎矩最小能量法的公式;梁鵬[4]提出通過調(diào)整結(jié)構(gòu)剛度來近似的求解最小彎曲能量法。傳統(tǒng)的調(diào)索方法因受限于各自的局限性而不具有普遍適用性,如剛性支撐連續(xù)梁法不能考慮多條件的約束問題,因而用其計算的索力值往往不合理;彎曲最小能量法在結(jié)合影響矩陣情況下,想要精確求解計算會十分繁瑣,而采用近似求解方法也需要不斷調(diào)節(jié)結(jié)構(gòu)的剛度才能使其達(dá)到預(yù)期結(jié)果。

    為了避免傳統(tǒng)調(diào)索方法的缺點,且能夠快速有效的求解出索力,出現(xiàn)了將計算機(jī)智能優(yōu)化算法與數(shù)學(xué)方法相結(jié)合的調(diào)索方法,Ha[5]等人提出微遺傳算法并對斜拉橋成橋索力進(jìn)行了優(yōu)化,陳志軍[6]、吳霄[7]等人分別采用傳統(tǒng)的粒子群和遺傳算法對斜拉橋索力進(jìn)行了優(yōu)化,但傳統(tǒng)的智能優(yōu)化算法因其自身的局限性導(dǎo)致結(jié)果往往不夠理想。為了使結(jié)果更加準(zhǔn)確,本研究以實際工程項目為依托,通過對粒子群優(yōu)化算法加以改進(jìn)并與影響矩陣相結(jié)合,在同時兼顧計算效率和多目標(biāo)約束的條件下求解吊桿索力,并與其他多種方法進(jìn)行分析比較,從而驗證該方法的可行性。

    1 基于改進(jìn)粒子群算法的索力優(yōu)化

    1.1 粒子群算法原理

    粒子群算法最早是由Kennedy和Eberhart[8]提出的一種基于群鳥覓食的仿生智能算法。其原理是群鳥在覓食的過程中,每只鳥的初始位置和飛行方向都處于隨機(jī)的狀態(tài),而且不知道最佳的覓食點在何處。鳥群間通過相互學(xué)習(xí)、信息共享,并在每次覓食過程中結(jié)合自身經(jīng)驗和種群之間的信息傳遞不斷調(diào)整速度和位置,最終步步逼近食物,示意圖如圖1所示。

    圖1 粒子尋優(yōu)示意圖Fig.1 Schematic diagram of particle optimization

    如果將該方法運用到實際應(yīng)用中,則是把鳥群抽象為m個沒有質(zhì)量和體積的粒子,再將其延伸至N維空間,粒子在N維空間中的位置表示為Xi=(xi1,xi2,…,xij),速度表示為Vi=(vi1,vi2,…,vij),且其位置和速度需限制在[Xmin,Xmax]和[Vmin,Vmax]范圍內(nèi)。每個粒子都有一個由目標(biāo)函數(shù)決定的適應(yīng)度值fitness,并且知道自己迄今為止發(fā)現(xiàn)的最好位置pb=(pi1,pi2,…,pij)和群體發(fā)現(xiàn)的最好位置gb=(g1,g2,…,gj)。粒子通過追蹤pb和gb來更新自己的速度和位置,且速度和位置更新公式為:

    (1)

    (2)

    式中,i為粒子的個數(shù),i=1,2,…,m;j為例子的維度,j=1,2,…,N;k為迭代步數(shù);w是慣性因子;c1,c2為學(xué)習(xí)因子;r1,r2為[0,1]范圍內(nèi)的均勻隨機(jī)數(shù);vmax為用戶自定義常數(shù),一般取vmax=aXmax,vmin=aXmin,通常a為0.1~0.2。

    1.2 調(diào)索目標(biāo)函數(shù)

    目標(biāo)函數(shù)下每一個粒子都有一個適應(yīng)度值,它作為每次迭代的臨時儲存變量,并不具有記憶性,優(yōu)化算法可通過粒子每次迭代的適應(yīng)度值來評價結(jié)果的好壞,進(jìn)而不斷地尋優(yōu)最終靠近目標(biāo)。實踐證明[9-10],通過限制拱橋結(jié)構(gòu)的彎曲能量所求解的吊桿索力能使橋梁結(jié)構(gòu)處于一個較好的受力狀態(tài),然而拱橋在受彎的同時也受壓,因此本研究以系桿拱橋的系梁彎曲應(yīng)變能和拱肋壓縮應(yīng)變能之和為目標(biāo),并將結(jié)構(gòu)離散化,離散后的應(yīng)變能公式為:

    (3)

    式中,U為結(jié)構(gòu)的應(yīng)變能之和;m為離散單元的數(shù)量;Li,Ai,Ei,Ii,Mi,Ni分別為單元的長度、面積、材料的彈性模量、截面慣性矩、彎矩和軸力。

    為了能夠建立起設(shè)計變量(拱橋的吊桿索力)與優(yōu)化目標(biāo)間的函數(shù)關(guān)系,可采用影響矩陣原理作為連接二者的橋梁,進(jìn)而達(dá)到索力優(yōu)化的目的。影響矩陣公式具體表示為:

    AY=D,

    (4)

    式中,Y為施調(diào)向量,在本研究中是指吊桿索力組成的列向量;D為受調(diào)向量,是指結(jié)構(gòu)中關(guān)心截面上若干獨立元素所組成的列向量。這些元素一般是指關(guān)心截面的內(nèi)力、位移,通過將這些元素調(diào)整到期望狀態(tài)來求解施調(diào)向量Y。A為影響矩陣,是指當(dāng)施調(diào)向量的某一向量發(fā)生單位變化時,引起受調(diào)向量D的變化量。

    從式(3)可以看出,當(dāng)結(jié)構(gòu)確定時,其單元及材料的具體特性已經(jīng)確定,因而只能通過改變截面彎矩和軸力來調(diào)整應(yīng)變能的大小。因此不妨令受調(diào)向量為截面的彎矩和軸力,且調(diào)索前彎矩向量為M0、軸力向量為N0,施調(diào)向量X=[S1,S2,…Sn]T,其中S為吊桿索力,則調(diào)索后的彎矩M和軸力N為:

    (5)

    將式(5)代入式(3)可得結(jié)構(gòu)的目標(biāo)函數(shù)為:

    U=M0TBM0+M0TBC1X+

    XTC1TBM0+XTC1TBC1X+

    N0TQN0+N0TQC2X+XT

    C2TQN0+XTC2TQC2X,

    (6)

    式中,矩陣B和Q均為系數(shù)矩陣,是對角陣,分別表示單元柔度對單元彎矩和軸力的加權(quán)矩陣,bii=li/2EiIi、qii=li/2EiAi,C1、C2分別表示單位索力引起彎矩和軸力變化的影響矩陣。

    1.3 索力約束條件

    對于系桿拱橋,其系梁主要承受彎矩,拱肋受軸力作用,因此就必需對其彎矩、軸力加以約束以使其落在允許的范圍內(nèi)。彎矩的約束條件可表示為Mmin≤M≤Mmax,軸力的約束條件可表示為Nmin≤N≤Nmax。根據(jù)影響矩陣?yán)碚?,在對彎矩、軸力進(jìn)行限制的同時也相當(dāng)于對索力進(jìn)行了約束,從而使調(diào)索后的最優(yōu)解具有了合理性。即當(dāng)Mmin≤M≤Mmax、Nmin≤N≤Nmax時,可得索力約束條件為:

    式中,Si為第i根吊桿的索力;Mjmax,Mjmin分別為系梁上第j個單元所允許的最大、最小彎矩;M0j為系梁第j個單元初始彎矩;Njmax,Njmin分別為拱肋上第j個單元所允許的最大、最小軸力;N0j為拱肋第j個單元初始軸力。同時為了滿足索力均勻原則,需將相鄰索力的差值限定在可控范圍內(nèi),另對索力做以下約束:

    (9)

    式中Δ為索力均勻度的評價閾值。

    1.4 粒子群算法改進(jìn)

    由于傳統(tǒng)粒子群算法的程序?qū)崿F(xiàn)過程十分簡潔,且需要調(diào)整的參數(shù)較少,因此在隨機(jī)優(yōu)化算法中具有強(qiáng)勁的優(yōu)勢,并被廣泛應(yīng)用于工程領(lǐng)域。但同時它也存在收斂速度慢、局部搜索能力差等缺點。為了提高算法的適應(yīng)范圍,使其結(jié)果更加準(zhǔn)確可靠,為此我們需要對傳統(tǒng)的算法進(jìn)行改進(jìn)[11-13]。研究結(jié)果表明[14-18],當(dāng)w較大時,算法的全局搜素能力較強(qiáng),w較小時算法的局部搜索能力較強(qiáng),因此算法的優(yōu)化問題便可轉(zhuǎn)化成一個w值的選取問題,即當(dāng)粒子目標(biāo)值趨于局部最優(yōu)時,需要增大慣性因子;當(dāng)粒子目標(biāo)值比較分散時,則需減小慣性因子。

    為了合理選取w,本研究將權(quán)重進(jìn)行自適應(yīng)調(diào)整,通過判斷粒子當(dāng)前的目標(biāo)函數(shù)值U的好壞來自動調(diào)整w,即當(dāng)U>Uavg,此時粒子的目標(biāo)函數(shù)值要比平均目標(biāo)值差,為了使其向較好的搜索區(qū)靠攏,對應(yīng)的慣性因子要較大;當(dāng)U′avg

    (10)

    式中,wmax,wmin分別為w的最大和最小值;Umin,Uavg分別表示當(dāng)前所有粒子的最小適應(yīng)值和平均適應(yīng)值。由于本研究是以彎曲能量的最小值為目標(biāo),因此U′avg是將所有小于Uavg的適應(yīng)值取平均所得。

    1.5 調(diào)索路線優(yōu)化

    本研究將索力調(diào)整和計算機(jī)智能優(yōu)化算法結(jié)合,并考慮了慣性權(quán)重的實時變化,使索力調(diào)整過程自動化,從而實現(xiàn)索力的高效尋優(yōu)。其具體步驟如下:

    (1)利用有限元軟件建立拱橋模型,并提取相關(guān)結(jié)構(gòu)數(shù)據(jù)及影響矩陣作為基礎(chǔ)數(shù)據(jù)庫。

    (2)建立目標(biāo)函數(shù),利用Matlab軟件編程粒子群算法并搜索在約束范圍內(nèi)的最優(yōu)解,如果滿足停止條件,則輸出當(dāng)前解,否則更新權(quán)重繼續(xù)搜索。

    (3)將搜索值回代有限元模型,校核結(jié)果的有效性。

    本研究優(yōu)化后的調(diào)索路線具體計算流程如圖2所示。

    圖2 計算流程圖Fig.2 Flowchart of calculation

    2 工程應(yīng)用

    2.1 工程實例

    本研究以某剛性系桿剛性拱橋位工程實例,主橋總長160 m,計算跨徑156.28 m,橋梁總寬13.9 m,拱軸線為二次拋物線,矢跨比為1/5。拱肋采用啞鈴型鋼管混凝土,采用Q345q鋼材,內(nèi)充C50微膨脹混凝土。系梁采用箱形斷面,系梁和橫梁為預(yù)應(yīng)力混凝土結(jié)構(gòu),橋面2%橫坡通過橫梁高度的變化進(jìn)行調(diào)整,吊桿間距為7.8 m,每片拱肋設(shè)吊桿19根,吊桿編號由小里程到大里程依次為1#~19#,如圖3所示。

    圖3 吊桿編號示意圖Fig.3 Schematic diagram of suspender numbering

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

    本研究系桿拱橋采用有限元軟件Abaqus建模,如圖4所示。除吊桿采用桁架單元外,其余均采用梁單元。恒載作用下主梁的初始彎矩值M0和拱肋軸力值N0通過軟件的數(shù)據(jù)提取功能來提取。通過分別改變每根吊桿的初拉力,使其發(fā)生單位1的變化,得到關(guān)鍵節(jié)點的彎矩和軸力變化值從而組成影響矩陣。模型的主要材料特性見表1。

    圖4 系桿拱有限元模型Fig.4 Finite element model of tied arch

    表1 材料特性Tab.1 Material properties

    利用Matlab工具編寫改進(jìn)前和改進(jìn)后的粒子群算法程序并進(jìn)行迭代求解,改進(jìn)前后除慣性權(quán)重外,其他系數(shù)取值相同。目前w較典型的取值范圍是0.4~0.9,算法在此范圍內(nèi)取值能獲得較好的局部和全局收斂性能。由1.4可知,為使算法能在初期計算時找到較好的解,需要取較大的慣性權(quán)重值,因此本研究對改進(jìn)前的慣性權(quán)重w恒取0.9,改進(jìn)后慣性權(quán)重系數(shù)按式(9)計算,學(xué)習(xí)因子取c1、c2取Bergh F[19]的推薦值并精確至小數(shù)點后兩位,種群規(guī)模m取50,索力范圍取[500,1 500],速度范圍取索力范圍的15%,具體數(shù)據(jù)見表2。由圖5、圖6可以看出,在迭代次數(shù)都為300次的條件下,改進(jìn)后的算法收斂速度得到明顯提高,在收斂精度方面也有一定的改善。

    表2 改進(jìn)后粒子群算法參數(shù)Tab.2 Improved PSO parameters

    圖5 算法改進(jìn)前后效果對比Fig.5 Comparison of effects before and after algorithm improvement

    圖6 算法改進(jìn)前后索力對比Fig.6 Comparison of cable forces before and after algorithm improvement

    圖7 索力對比Fig.7 Comparison of cable forces

    2.3 優(yōu)化結(jié)果分析

    如何控制拱橋的內(nèi)力和線形使其在安全的范圍內(nèi)是施工監(jiān)控過程中的重要環(huán)節(jié)。為了能夠清楚地展示粒子群優(yōu)化算法在實際工程中的效果,本節(jié)除使用研究中提出的改進(jìn)粒子群算法外,另使用近似的彎曲能量法、剛性支撐連續(xù)梁法及未知荷載系數(shù)法對索力進(jìn)行求解,4種方法求得的索力見圖7。從圖7中可以看出,利用改進(jìn)后的粒子群優(yōu)化算法求得的索力值較其他3種方法更加均勻適中。

    將4種方法求得的索力再分別代入有限元模型中,得到各索力下系梁彎矩和豎向位移,如圖8、圖9所示。從圖8可以看出,用不同方法約束的系梁彎矩在數(shù)值和均勻度上都有較大的差異,其中未知荷載系數(shù)法約束的彎矩上下峰值差距較大,剛性支撐連續(xù)梁法約束的彎矩出現(xiàn)驟增現(xiàn)象,這與該方法在吊桿數(shù)量有限情況下,可能會導(dǎo)致結(jié)構(gòu)其余部位出現(xiàn)狀態(tài)異常有關(guān),而本研究的改進(jìn)粒子群優(yōu)化法和近似最小能量法優(yōu)化后的系梁彎矩都比較均勻,但本研究方法得到的結(jié)果比近似最小彎曲能量法得到的結(jié)果在數(shù)值上更小。從圖9可以看出,剛性支撐連續(xù)梁法和近似最小彎曲能量法優(yōu)化后的豎向位移較大,本研究的改進(jìn)粒子群優(yōu)化法和未知荷載系數(shù)法優(yōu)化后的豎向位移相近,但本研究的改進(jìn)粒子群優(yōu)化法較未知荷載系數(shù)法優(yōu)化后的豎向位移更加均勻。4種方法對結(jié)構(gòu)具體影響見表3。

    圖8 系梁彎矩圖Fig.8 Bending moments of tie beam

    圖9 豎向位移Fig.9 Vertical displacements

    表3 四種優(yōu)化方法對結(jié)構(gòu)的影響Tab.3 Influence of 4 optimization methods on structure

    3 結(jié)論

    本研究在粒子群優(yōu)化算法的基礎(chǔ)上對其慣性權(quán)重進(jìn)行改進(jìn),同時考慮了系桿拱橋在索力優(yōu)化過程中的主要約束條件,提出了一種新的系桿拱橋索力優(yōu)化的方法,并以某系桿拱橋為工程背景進(jìn)行應(yīng)用,主要結(jié)論如下:

    (1)基于粒子群算法的調(diào)索方法可以借助有限元軟件和數(shù)值分析軟件實現(xiàn)索力的自動化、智能化求解。在改進(jìn)后算法的尋優(yōu)能力得到提高,從而避免了繁瑣的試算過程,進(jìn)而可以高效、準(zhǔn)確地求解出索力,極大提高了調(diào)索效率。

    (2)粒子群算法改進(jìn)后較改進(jìn)前在收斂速度和精度上有一定的改善。

    (3)4種索力優(yōu)化方法的結(jié)果表明,改進(jìn)后的粒子群算法較其他3種方法能更加全面地考慮影響因素的作用,進(jìn)而規(guī)避不利因素對優(yōu)化結(jié)果的影響,以達(dá)到更好的優(yōu)化效果。

    (4)本研究提出的方法可進(jìn)一步推廣于各施工階段過程中吊桿初拉力的確定問題。

    猜你喜歡
    系桿索力吊桿
    江蘇索力得新材料集團(tuán)有限公司
    反對稱異型系桿拱橋設(shè)計與研究
    上海公路(2019年2期)2019-10-08 09:05:44
    某下承式系桿拱橋設(shè)計分析
    淺談下承式系桿拱結(jié)構(gòu)設(shè)計及計算分析
    飛燕式鋼箱系桿拱橋吊桿疲勞性能影響因素研究
    12CrlMoV鍋爐吊桿用鋼的開發(fā)生產(chǎn)實踐
    山東冶金(2018年5期)2018-11-22 05:12:06
    基于差值法的系桿拱橋更換吊桿方案優(yōu)化研究
    鐵道建筑(2018年7期)2018-08-01 06:29:10
    系桿拱橋穩(wěn)定性探索——考慮拱梁共同作用
    上海公路(2018年3期)2018-03-21 05:55:42
    預(yù)應(yīng)力鋼絞線網(wǎng)加固混凝土橋梁的索力分布試驗研究
    基于拉索振動特征的索力檢測與評估方法
    亚洲av片天天在线观看| 夫妻午夜视频| 91成人精品电影| 高清av免费在线| 国产精品香港三级国产av潘金莲 | 午夜久久久在线观看| 日韩精品免费视频一区二区三区| 少妇猛男粗大的猛烈进出视频| 岛国毛片在线播放| 悠悠久久av| 国产麻豆69| 大香蕉久久网| 亚洲av在线观看美女高潮| 国产激情久久老熟女| 交换朋友夫妻互换小说| 午夜福利乱码中文字幕| 欧美另类一区| kizo精华| 9191精品国产免费久久| 精品一区在线观看国产| svipshipincom国产片| 午夜老司机福利片| 久久这里只有精品19| 男女免费视频国产| 性色av乱码一区二区三区2| 国产在线视频一区二区| avwww免费| 老司机靠b影院| 久久国产精品影院| 亚洲精品av麻豆狂野| 啦啦啦在线免费观看视频4| 黑人猛操日本美女一级片| 又大又黄又爽视频免费| 69精品国产乱码久久久| 国产麻豆69| 午夜老司机福利片| 欧美日韩成人在线一区二区| 久9热在线精品视频| 精品国产乱码久久久久久小说| 久久精品国产a三级三级三级| 亚洲精品一卡2卡三卡4卡5卡 | 老汉色av国产亚洲站长工具| 国产日韩欧美在线精品| 91字幕亚洲| 久久精品久久久久久久性| 性色av一级| 一级黄色大片毛片| 中文乱码字字幕精品一区二区三区| 高清欧美精品videossex| 老司机午夜十八禁免费视频| 国产精品久久久久久人妻精品电影 | 一级毛片我不卡| av一本久久久久| 国产成人啪精品午夜网站| 亚洲伊人久久精品综合| 美女扒开内裤让男人捅视频| 国产一区亚洲一区在线观看| 美女扒开内裤让男人捅视频| 男女国产视频网站| 亚洲av成人精品一二三区| bbb黄色大片| 91国产中文字幕| 人妻 亚洲 视频| 国产精品免费视频内射| 成人国语在线视频| 国产精品.久久久| 老司机深夜福利视频在线观看 | 免费少妇av软件| 男男h啪啪无遮挡| 亚洲人成电影免费在线| 国产三级黄色录像| 99久久人妻综合| 国产一区二区 视频在线| 嫁个100分男人电影在线观看 | 国产欧美日韩综合在线一区二区| 一级片免费观看大全| 丝袜美腿诱惑在线| 香蕉丝袜av| 精品免费久久久久久久清纯 | 啦啦啦在线观看免费高清www| 国产免费一区二区三区四区乱码| www.自偷自拍.com| 亚洲人成77777在线视频| 日本av免费视频播放| 成人三级做爰电影| 国产精品久久久久成人av| 99久久人妻综合| 国产在线视频一区二区| 青草久久国产| 久久青草综合色| 亚洲成人免费电影在线观看 | 天堂8中文在线网| 亚洲欧美中文字幕日韩二区| 80岁老熟妇乱子伦牲交| 成年人黄色毛片网站| 久久99热这里只频精品6学生| 蜜桃国产av成人99| 波野结衣二区三区在线| 桃花免费在线播放| 日本欧美国产在线视频| 国产主播在线观看一区二区 | 亚洲精品国产区一区二| 国产精品二区激情视频| av福利片在线| netflix在线观看网站| 男人舔女人的私密视频| 老汉色av国产亚洲站长工具| 日本vs欧美在线观看视频| 大片免费播放器 马上看| 老熟女久久久| 最近中文字幕2019免费版| 桃花免费在线播放| 五月天丁香电影| 国产xxxxx性猛交| 最近最新中文字幕大全免费视频 | 免费在线观看完整版高清| 最新的欧美精品一区二区| 男女高潮啪啪啪动态图| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩综合久久久久久| 女人精品久久久久毛片| 成年人午夜在线观看视频| 成年动漫av网址| a级毛片在线看网站| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 日本黄色日本黄色录像| 精品熟女少妇八av免费久了| 精品久久久精品久久久| 日本欧美视频一区| 18在线观看网站| 欧美久久黑人一区二区| 日韩制服丝袜自拍偷拍| 亚洲七黄色美女视频| 国产男人的电影天堂91| 国产视频一区二区在线看| 1024香蕉在线观看| 亚洲av成人精品一二三区| 丰满少妇做爰视频| 久久久精品免费免费高清| 1024视频免费在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av国产av综合av卡| 美女福利国产在线| 国产真人三级小视频在线观看| 永久免费av网站大全| 国产精品欧美亚洲77777| 国产亚洲av高清不卡| 亚洲,欧美精品.| 成人手机av| 日本av手机在线免费观看| 欧美 日韩 精品 国产| 性高湖久久久久久久久免费观看| 大码成人一级视频| 建设人人有责人人尽责人人享有的| 亚洲国产精品999| 成人黄色视频免费在线看| 免费在线观看日本一区| 免费黄频网站在线观看国产| 国产成人精品无人区| 夫妻午夜视频| 亚洲 国产 在线| 大陆偷拍与自拍| 18禁裸乳无遮挡动漫免费视频| 高清欧美精品videossex| 在线观看免费视频网站a站| 这个男人来自地球电影免费观看| 日日夜夜操网爽| 午夜免费成人在线视频| a级毛片在线看网站| 久久久久久久久免费视频了| 国产在线一区二区三区精| 精品免费久久久久久久清纯 | 久久亚洲国产成人精品v| 日本色播在线视频| 99久久综合免费| 国产精品麻豆人妻色哟哟久久| 国产精品久久久久久人妻精品电影 | 国产熟女欧美一区二区| 中文字幕最新亚洲高清| 日本五十路高清| 最近最新中文字幕大全免费视频 | 国产精品一区二区在线观看99| 狂野欧美激情性xxxx| 亚洲精品第二区| 十八禁网站网址无遮挡| 国产男人的电影天堂91| 狠狠精品人妻久久久久久综合| 国产免费福利视频在线观看| 久久99热这里只频精品6学生| 新久久久久国产一级毛片| 伊人亚洲综合成人网| 国产日韩一区二区三区精品不卡| 精品亚洲成国产av| 免费看十八禁软件| 日日夜夜操网爽| 久久人人爽人人片av| 亚洲精品美女久久av网站| 少妇的丰满在线观看| 欧美成人午夜精品| 国产av精品麻豆| 国产男女内射视频| 乱人伦中国视频| 亚洲一码二码三码区别大吗| 国产成人精品久久二区二区免费| 91麻豆精品激情在线观看国产 | 欧美日韩视频高清一区二区三区二| 久久精品国产a三级三级三级| 亚洲一码二码三码区别大吗| 好男人视频免费观看在线| 日韩视频在线欧美| 老司机午夜十八禁免费视频| 婷婷成人精品国产| 精品免费久久久久久久清纯 | 五月开心婷婷网| av又黄又爽大尺度在线免费看| 秋霞在线观看毛片| 久久中文字幕一级| 国产av国产精品国产| 国产精品麻豆人妻色哟哟久久| 久久性视频一级片| 老司机亚洲免费影院| 人人妻人人爽人人添夜夜欢视频| 性少妇av在线| 久久精品久久久久久久性| 欧美激情极品国产一区二区三区| 免费观看a级毛片全部| 国产免费视频播放在线视频| 午夜两性在线视频| 亚洲第一av免费看| 亚洲av成人精品一二三区| 日韩,欧美,国产一区二区三区| 国产精品秋霞免费鲁丝片| 好男人视频免费观看在线| 国产av一区二区精品久久| 久久久久久久久久久久大奶| 狂野欧美激情性xxxx| www日本在线高清视频| 美国免费a级毛片| 婷婷色综合www| av在线app专区| 欧美日韩成人在线一区二区| 在现免费观看毛片| 黄片小视频在线播放| 两人在一起打扑克的视频| 亚洲欧美一区二区三区久久| 午夜福利,免费看| 国产精品秋霞免费鲁丝片| 午夜日韩欧美国产| 国产不卡av网站在线观看| 91国产中文字幕| 无限看片的www在线观看| 精品一品国产午夜福利视频| 亚洲精品自拍成人| 欧美性长视频在线观看| 亚洲精品一二三| av片东京热男人的天堂| 美女国产高潮福利片在线看| 男女无遮挡免费网站观看| 亚洲人成电影免费在线| 日本欧美国产在线视频| 性高湖久久久久久久久免费观看| 精品久久蜜臀av无| 国产高清视频在线播放一区 | 九草在线视频观看| 精品国产乱码久久久久久小说| 一二三四在线观看免费中文在| 另类亚洲欧美激情| 永久免费av网站大全| 中文字幕色久视频| 国产精品一二三区在线看| 黄频高清免费视频| 色综合欧美亚洲国产小说| 久久久久久久大尺度免费视频| 19禁男女啪啪无遮挡网站| 嫁个100分男人电影在线观看 | 人妻一区二区av| 久久精品亚洲av国产电影网| 色婷婷久久久亚洲欧美| av有码第一页| 亚洲综合色网址| 日韩中文字幕欧美一区二区 | 91精品伊人久久大香线蕉| 乱人伦中国视频| 日韩大码丰满熟妇| 中文字幕av电影在线播放| 婷婷色av中文字幕| 国产免费福利视频在线观看| 国产亚洲欧美精品永久| 国产片特级美女逼逼视频| 黄色视频不卡| 成年女人毛片免费观看观看9 | 乱人伦中国视频| 亚洲人成77777在线视频| 国产一区二区 视频在线| 亚洲精品在线美女| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡动漫免费视频| a 毛片基地| 国产成人精品在线电影| 首页视频小说图片口味搜索 | 深夜精品福利| 一区二区三区精品91| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播 | 国产激情久久老熟女| 亚洲欧美中文字幕日韩二区| 男女午夜视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 美女扒开内裤让男人捅视频| 亚洲成人免费电影在线观看 | 欧美少妇被猛烈插入视频| 久久 成人 亚洲| 中文字幕精品免费在线观看视频| 亚洲午夜精品一区,二区,三区| 婷婷色av中文字幕| 亚洲av美国av| 嫁个100分男人电影在线观看 | 欧美97在线视频| 日本a在线网址| 成人18禁高潮啪啪吃奶动态图| 国产成人a∨麻豆精品| 成人国产一区最新在线观看 | 18禁裸乳无遮挡动漫免费视频| 久久热在线av| 亚洲av电影在线观看一区二区三区| 校园人妻丝袜中文字幕| 乱人伦中国视频| 一二三四社区在线视频社区8| 最近最新中文字幕大全免费视频 | 亚洲欧美精品综合一区二区三区| 99久久人妻综合| 亚洲av日韩在线播放| 日韩一卡2卡3卡4卡2021年| 啦啦啦中文免费视频观看日本| 男人爽女人下面视频在线观看| 好男人电影高清在线观看| 久久久精品免费免费高清| 青草久久国产| 久久国产亚洲av麻豆专区| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 精品久久蜜臀av无| 亚洲欧美日韩高清在线视频 | 成人亚洲欧美一区二区av| 国产成人欧美在线观看 | 一二三四社区在线视频社区8| av一本久久久久| 亚洲国产成人一精品久久久| 亚洲精品久久久久久婷婷小说| 丝袜美足系列| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 国产精品久久久久久精品电影小说| 丝袜美腿诱惑在线| 丰满迷人的少妇在线观看| 久久狼人影院| 精品国产一区二区久久| 久久精品久久久久久久性| 一级毛片黄色毛片免费观看视频| 久久国产亚洲av麻豆专区| 在线观看人妻少妇| 亚洲成国产人片在线观看| 十八禁高潮呻吟视频| 母亲3免费完整高清在线观看| 最近手机中文字幕大全| 国产欧美日韩精品亚洲av| √禁漫天堂资源中文www| 亚洲av在线观看美女高潮| 日本午夜av视频| 老司机影院成人| 人体艺术视频欧美日本| 久久热在线av| 亚洲av美国av| 欧美日韩亚洲国产一区二区在线观看 | www日本在线高清视频| 叶爱在线成人免费视频播放| 女人爽到高潮嗷嗷叫在线视频| 久久天堂一区二区三区四区| 国产高清videossex| 亚洲伊人色综图| 中文字幕人妻熟女乱码| 精品福利观看| 国产成人啪精品午夜网站| 亚洲av电影在线观看一区二区三区| 久久人人97超碰香蕉20202| 国产精品久久久久久人妻精品电影 | 欧美日韩av久久| 你懂的网址亚洲精品在线观看| 国产日韩一区二区三区精品不卡| 午夜福利视频在线观看免费| 久久精品国产亚洲av高清一级| 欧美精品啪啪一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| 自线自在国产av| 大型av网站在线播放| 人人澡人人妻人| 国产91精品成人一区二区三区 | 国产又爽黄色视频| 国产成人av教育| 国产爽快片一区二区三区| 亚洲黑人精品在线| 高清不卡的av网站| 亚洲精品久久成人aⅴ小说| 日本vs欧美在线观看视频| 丝袜脚勾引网站| 美女午夜性视频免费| 汤姆久久久久久久影院中文字幕| 91精品国产国语对白视频| 19禁男女啪啪无遮挡网站| 久9热在线精品视频| 久久人人爽av亚洲精品天堂| 国产日韩欧美视频二区| 亚洲久久久国产精品| 少妇的丰满在线观看| 少妇精品久久久久久久| 69精品国产乱码久久久| 国产精品一二三区在线看| 成年人黄色毛片网站| 黄色 视频免费看| 久9热在线精品视频| 国产欧美日韩一区二区三区在线| 午夜老司机福利片| 亚洲精品国产区一区二| 99热国产这里只有精品6| 999久久久国产精品视频| 一级片免费观看大全| 亚洲欧美精品综合一区二区三区| 青草久久国产| 少妇人妻 视频| 国产精品久久久av美女十八| 欧美日韩精品网址| 国产欧美日韩一区二区三区在线| 亚洲图色成人| 国产午夜精品一二区理论片| 男女高潮啪啪啪动态图| 伦理电影免费视频| 青春草亚洲视频在线观看| 国产黄频视频在线观看| cao死你这个sao货| 狠狠婷婷综合久久久久久88av| 久久久精品区二区三区| av电影中文网址| 欧美久久黑人一区二区| 婷婷丁香在线五月| 国产成人av教育| 国产一区二区三区av在线| 侵犯人妻中文字幕一二三四区| 久久99一区二区三区| 老汉色av国产亚洲站长工具| 欧美xxⅹ黑人| bbb黄色大片| 丰满少妇做爰视频| 人人妻,人人澡人人爽秒播 | 老司机在亚洲福利影院| 80岁老熟妇乱子伦牲交| 久9热在线精品视频| 大片免费播放器 马上看| 极品少妇高潮喷水抽搐| 免费观看av网站的网址| 亚洲欧洲日产国产| 一区二区三区激情视频| 亚洲,欧美,日韩| 精品国产乱码久久久久久小说| 欧美日韩视频高清一区二区三区二| 飞空精品影院首页| 看十八女毛片水多多多| 国产精品秋霞免费鲁丝片| 无限看片的www在线观看| 男人添女人高潮全过程视频| 色精品久久人妻99蜜桃| 丁香六月欧美| 精品亚洲成a人片在线观看| 精品久久久久久久毛片微露脸 | 久久久欧美国产精品| 日韩大码丰满熟妇| 精品国产国语对白av| 国产淫语在线视频| 久久精品成人免费网站| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 婷婷成人精品国产| 巨乳人妻的诱惑在线观看| 国产精品人妻久久久影院| www.熟女人妻精品国产| 久久 成人 亚洲| 我的亚洲天堂| 亚洲午夜精品一区,二区,三区| 免费日韩欧美在线观看| 一个人免费看片子| 亚洲七黄色美女视频| 日韩av免费高清视频| 又黄又粗又硬又大视频| 天天躁夜夜躁狠狠躁躁| 男女床上黄色一级片免费看| 丝袜人妻中文字幕| 亚洲国产av影院在线观看| 国产精品欧美亚洲77777| 热re99久久国产66热| 成年人午夜在线观看视频| 嫩草影视91久久| 国产成人一区二区在线| 视频区欧美日本亚洲| 欧美人与性动交α欧美精品济南到| 啦啦啦在线观看免费高清www| 久久精品久久精品一区二区三区| 亚洲国产成人一精品久久久| 在线精品无人区一区二区三| 在线av久久热| 老司机在亚洲福利影院| 捣出白浆h1v1| bbb黄色大片| 一级黄片播放器| av欧美777| 视频在线观看一区二区三区| xxx大片免费视频| 蜜桃国产av成人99| 久久99一区二区三区| 伦理电影免费视频| 自线自在国产av| 大陆偷拍与自拍| 久久人人爽av亚洲精品天堂| 老汉色∧v一级毛片| 国产免费福利视频在线观看| 一边亲一边摸免费视频| 午夜视频精品福利| 欧美日韩亚洲综合一区二区三区_| 亚洲伊人久久精品综合| 亚洲精品第二区| 欧美日本中文国产一区发布| 操美女的视频在线观看| 亚洲精品第二区| 啦啦啦中文免费视频观看日本| 国产又爽黄色视频| 麻豆av在线久日| 女人爽到高潮嗷嗷叫在线视频| 一个人免费看片子| 69精品国产乱码久久久| 亚洲精品成人av观看孕妇| 捣出白浆h1v1| 欧美日韩亚洲国产一区二区在线观看 | 久久精品亚洲av国产电影网| 91精品国产国语对白视频| 欧美黑人精品巨大| 两性夫妻黄色片| 五月天丁香电影| 人妻 亚洲 视频| 18禁黄网站禁片午夜丰满| 美女福利国产在线| 91麻豆精品激情在线观看国产 | 国产一区二区三区av在线| 久久精品国产综合久久久| 亚洲久久久国产精品| 国产一区二区在线观看av| 久久久久久免费高清国产稀缺| 国产男女超爽视频在线观看| 在线观看国产h片| 精品人妻一区二区三区麻豆| 亚洲一码二码三码区别大吗| 9191精品国产免费久久| 久久精品aⅴ一区二区三区四区| 国产极品粉嫩免费观看在线| 在线 av 中文字幕| 波野结衣二区三区在线| 免费av中文字幕在线| 精品亚洲乱码少妇综合久久| 99国产精品免费福利视频| 黑人猛操日本美女一级片| 亚洲av在线观看美女高潮| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性bbbbbb| 国产精品av久久久久免费| 久久久久久久精品精品| av视频免费观看在线观看| 可以免费在线观看a视频的电影网站| 999久久久国产精品视频| 男女之事视频高清在线观看 | av网站在线播放免费| 免费黄频网站在线观看国产| 欧美人与性动交α欧美精品济南到| 免费日韩欧美在线观看| 国产高清不卡午夜福利| 亚洲欧洲精品一区二区精品久久久| 日韩视频在线欧美| 两性夫妻黄色片| a 毛片基地| 国产熟女欧美一区二区| 国产精品一区二区在线不卡| 午夜免费男女啪啪视频观看| 国产成人av教育| 一边摸一边做爽爽视频免费| 91精品三级在线观看| 操美女的视频在线观看| 亚洲 国产 在线| 一级片免费观看大全| 97人妻天天添夜夜摸| 人人妻人人澡人人爽人人夜夜| 看免费av毛片| 日本黄色日本黄色录像| 亚洲精品国产一区二区精华液| 午夜免费成人在线视频| 久久精品成人免费网站| 日韩 亚洲 欧美在线| 黑丝袜美女国产一区| 好男人视频免费观看在线| 国产三级黄色录像| 国产亚洲午夜精品一区二区久久| 中文欧美无线码| 成在线人永久免费视频| 制服人妻中文乱码| 亚洲精品一区蜜桃| 亚洲av成人精品一二三区| 国产成人欧美| 男女之事视频高清在线观看 | 日韩大片免费观看网站|