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

    基于雙策略改進(jìn)遺傳算法的斜塔斜拉橋索力優(yōu)化分析

    2021-03-31 01:01:42鄭凱鋒曾瓊瑤
    四川建筑 2021年1期
    關(guān)鍵詞:橋塔索力斜拉橋

    潘 偉, 鄭凱鋒, 曾瓊瑤

    (西南交通大學(xué)土木工程學(xué)院,四川成都 610031)

    斜拉橋是一種組合體系橋梁,由橋塔、主梁和拉索三個(gè)部分組成,拉索對(duì)主梁提供彈性支撐,使得主梁在截面較小的情況下,具備較大的跨越能力[1]。同時(shí),斜拉橋美學(xué)特點(diǎn)突出,被廣泛應(yīng)用于大跨度橋梁建設(shè)中。斜塔斜拉橋通過(guò)將橋塔傾斜,利用其不平衡力矩改善整體結(jié)構(gòu)的受力性能[2]。

    斜拉橋的合理成橋狀態(tài)是指在該狀態(tài)下,主梁和橋塔的線形符合設(shè)計(jì)要求,結(jié)構(gòu)彎矩較小,拉索受力均勻,結(jié)構(gòu)整體性能優(yōu)良。斜拉橋的合理成橋狀態(tài)一般可以通過(guò)調(diào)整拉索內(nèi)力來(lái)實(shí)現(xiàn)[3]。針對(duì)索力優(yōu)化問(wèn)題,相關(guān)學(xué)者提出多種不同方法進(jìn)行求解。例如指定結(jié)構(gòu)狀態(tài)的優(yōu)化方法,它包括:零位移法,剛性支撐連續(xù)法,零支反力法等。但該方法針對(duì)主、邊跨非對(duì)稱斜拉橋不適用,易造成主梁局部彎矩過(guò)大的問(wèn)題,且無(wú)法考慮橋塔的內(nèi)力和變形[4]。另一種思路是通過(guò)將索力優(yōu)化問(wèn)題轉(zhuǎn)化為數(shù)學(xué)模型,運(yùn)用數(shù)學(xué)優(yōu)化分析求解。于玲等[5]采用顯示梯度法針對(duì)混合梁斜拉橋合理索力進(jìn)行求解,戴杰等[6]采用有效約束集法優(yōu)化了某斜拉橋成橋索力。楊景瑜[7]將索力優(yōu)化問(wèn)題轉(zhuǎn)化為非線性規(guī)劃問(wèn)題并利用罰函數(shù)法進(jìn)行求解。一些學(xué)者將數(shù)學(xué)優(yōu)化領(lǐng)域的成熟算法引入索力調(diào)整中,取得了較好的效果。吳霄[8]運(yùn)用遺傳算法針對(duì)大跨度混合梁斜拉橋進(jìn)行了索力優(yōu)化,陳驍[9]結(jié)合神經(jīng)網(wǎng)絡(luò)與遺傳算法運(yùn)用于獨(dú)塔斜拉橋索力優(yōu)化和施工控制中,張玉平[10]提出了一種粒子群智能優(yōu)化算法的索力優(yōu)化方法,用于斜拉橋索力和梁拱體系橋梁吊桿力確定。針對(duì)斜塔斜拉橋非對(duì)稱結(jié)構(gòu),無(wú)法運(yùn)用傳統(tǒng)指定狀態(tài)方法進(jìn)行求解合理成橋索力,本文運(yùn)用遺傳算法,以某斜塔有背索斜拉橋?yàn)槔?,將全橋結(jié)構(gòu)彎曲應(yīng)變能最小化作為目標(biāo),結(jié)合python編程語(yǔ)言對(duì)通用有限元軟件ABAQUS進(jìn)行二次開(kāi)發(fā),求解合理成橋狀態(tài)索力,并采用較優(yōu)種群初始化策略與精英留存策略,對(duì)遺傳算法進(jìn)行改進(jìn)。

    1 優(yōu)化算法

    1.1 基礎(chǔ)理論

    遺傳算法是一種模擬生物界繁衍遺傳和自然選擇機(jī)制的搜索型算法,通過(guò)種群交叉、遺傳、變異及選優(yōu)產(chǎn)生子代種群,重復(fù)該過(guò)程直到種群滿足某種收斂要求[11]。針對(duì)拉索索力和結(jié)構(gòu)力學(xué)指標(biāo)之間的非線性映射關(guān)系,遺傳算法可在解空間進(jìn)行高效率尋優(yōu),避免一般調(diào)索過(guò)程中對(duì)經(jīng)驗(yàn)的依賴。

    遺傳算法的編碼策略一般可分為實(shí)數(shù)編碼和二進(jìn)制編碼,實(shí)數(shù)編碼的優(yōu)勢(shì)在于避免對(duì)種群基因進(jìn)行編碼和解碼操作,操作簡(jiǎn)單,提高程序運(yùn)算效率,但由于搜索區(qū)間連續(xù),容易陷入局部最優(yōu)解。二進(jìn)制編碼種群基因變化性較大,可以在更大空間內(nèi)進(jìn)行解尋優(yōu),但其收斂速度較慢,可以通過(guò)一些策略提升算法效率。

    選擇操作體現(xiàn)優(yōu)勝略汰的原理,根據(jù)評(píng)價(jià)函數(shù)判斷個(gè)體優(yōu)劣,對(duì)優(yōu)秀個(gè)體進(jìn)行選擇培育子代。采用錦標(biāo)賽策略,首先在種群中隨機(jī)抽取一定個(gè)體,之后中選出最優(yōu)個(gè)體,重復(fù)此過(guò)程,即可得到育種族群。

    交叉模擬種群然基因進(jìn)行交換的過(guò)程。采用均勻交叉算子,即父代基因均有概率交換,避免陷入局部最優(yōu)。

    變異可以擴(kuò)大種群多樣性,避免過(guò)早收斂。針對(duì)二進(jìn)制編碼,位翻轉(zhuǎn)突變是一種有效突變,可以充分發(fā)揮二進(jìn)制編碼尋求全局最優(yōu)解的特性。

    1.2 評(píng)價(jià)函數(shù)

    索力優(yōu)化的主要目標(biāo)即降低橋塔和主梁的彎矩和位移,可采用結(jié)構(gòu)的彎矩應(yīng)變能作為種群評(píng)價(jià)函數(shù),目標(biāo)函數(shù)為:

    (1)

    式中:U為結(jié)構(gòu)的彎曲應(yīng)變能,Mij、Mik表示第i號(hào)梁?jiǎn)卧猨、k端的彎矩,E為材料彈性模量,l為單元長(zhǎng)度,I為單元截面抗彎慣性矩。

    針對(duì)斜塔斜拉橋,橋塔的彎曲應(yīng)變能不可忽略,由于主梁、橋塔的目標(biāo)值存在數(shù)量級(jí)差異,可對(duì)橋塔項(xiàng)增加權(quán)重,同等的考慮兩者的效應(yīng),公式如下:

    minU=min(UB+φ·UT)

    (2)

    式中:UB、UT分別表示主梁和橋塔的彎曲應(yīng)變能值。

    1.3 雙策略改進(jìn)算法

    一般遺傳算法采用隨機(jī)方式生成初始種群,尋優(yōu)效率較低,通過(guò)指定一組較優(yōu)解作為生成初始種群的基準(zhǔn),可以較大幅度提升求解效率,同時(shí)采用二進(jìn)制編碼和合理的變異策略,可以避免陷入局部最優(yōu)解的情況。

    在簡(jiǎn)單遺傳算法中一般用子代替代全部父代,可能造成父代中優(yōu)秀個(gè)體流失,通過(guò)應(yīng)用精英保留策略,使優(yōu)秀的個(gè)體保留,可以使得算法收斂速度提升。

    本文通過(guò)引入上述精英保留策略和優(yōu)質(zhì)初始化種群策略(簡(jiǎn)稱雙策略)進(jìn)行算法改進(jìn),以提升算法的計(jì)算效率。算法流程如圖1所示。

    圖1 改進(jìn)算法流程

    2 計(jì)算實(shí)例

    2.1 工程概況

    某斜拉橋主結(jié)構(gòu)為(180+90+75)m單斜塔單索面斜拉橋,主梁采用單箱四室鋼箱梁截面,橋塔采用箱形截面混凝土橋塔,塔墩梁固結(jié),全橋共12對(duì)斜拉索呈不對(duì)成布置。邊跨主梁的斜拉索集中布置在輔助墩附近,標(biāo)準(zhǔn)索距3 m,主跨主梁上則均勻布置,標(biāo)準(zhǔn)索距12 m。結(jié)構(gòu)立面圖如圖2所示,結(jié)構(gòu)仿真計(jì)算參數(shù)如表1所示。主跨一側(cè)斜拉索由橋塔向邊墩方向依次編號(hào)為S01~S12,邊跨一側(cè)斜拉索由橋塔向邊墩方向依次編號(hào)為H01~H12。

    2.2 建立空間梁?jiǎn)卧P颓蠼庖淮温淞籂顟B(tài)

    首先運(yùn)用Midas軟件建立全橋一次落架梁?jiǎn)卧P?。橋塔、橋墩和主梁采用梁?jiǎn)卧M,斜拉索采用桁架單元模擬,分析模型如圖3所示,利用該模型采用一次落梁方法獲得橋梁狀態(tài)初始解。

    圖2 結(jié)構(gòu)立面(單位:m)

    表1 結(jié)構(gòu)計(jì)算參數(shù)

    圖3 Midas計(jì)算模型

    2.3 通過(guò)未知荷載系數(shù)法求解橋梁狀態(tài)較優(yōu)解

    結(jié)合斜塔斜拉橋的特點(diǎn),通過(guò)多次試算,指定約束條件:主梁位移小于200 mm;主梁彎矩小于15×105kN·m;橋塔彎矩小于4×104kN·m,運(yùn)用未知荷載系數(shù)法求解橋梁狀態(tài)較優(yōu)解,用于生成遺傳算法初始種群。一次落梁索力和較優(yōu)解索力對(duì)比如圖4所示。

    3 索力優(yōu)化

    3.1 建立ABAQUS模型

    運(yùn)用ABAQUS軟件建立分析模型如圖5所示。主梁、橋塔采用B33梁?jiǎn)卧?,斜拉索采用T3D2桁架單元模擬,主梁和橋塔與斜拉索之間通過(guò)MPC模擬主從約束,可以使之間剛性連桿只出現(xiàn)剛體轉(zhuǎn)動(dòng)[12],索力可采用初始狀態(tài)功能進(jìn)行設(shè)置,斜拉索垂度效應(yīng)采用等效彈性模量法考慮。

    圖5 ABAQUS計(jì)算模型

    3.2 索力優(yōu)化流程

    利用Python編程語(yǔ)言對(duì)ABAQUS進(jìn)行二次開(kāi)發(fā),使之能在不同初始索力條件下進(jìn)行多次迭代尋優(yōu),即首先將初始種群分別寫(xiě)入索力文本文檔,將索力信息作為初拉力對(duì)有限元模型進(jìn)行更新,調(diào)用ABAQUS求解器進(jìn)行求解,在后處理過(guò)程中提取主梁和橋塔單元的彎矩以及長(zhǎng)度用以計(jì)算種群適應(yīng)度函數(shù),即結(jié)構(gòu)彎曲應(yīng)變能函數(shù)值,根據(jù)種群內(nèi)個(gè)體的優(yōu)劣程度,進(jìn)行一系列選擇、交叉、變異操作,生成子代索力種群,對(duì)新的種群再次調(diào)用求解器求解,進(jìn)行多次迭代求解,最終得到索力最優(yōu)值。

    超參數(shù)設(shè)置:種群迭代次數(shù):40;初始種群數(shù)目:50;編碼方式:二進(jìn)制編碼;變異概率:0.5;交叉概率:0.4。

    迭代大約30次后,適應(yīng)度函數(shù)值趨于收斂,最終得到一組最優(yōu)索力,較優(yōu)解、最優(yōu)解索力對(duì)比如圖6所示。

    圖6 較優(yōu)解、最優(yōu)解索力計(jì)算

    3.3 計(jì)算結(jié)果分析

    取一次落架模型、較優(yōu)解索力模型和優(yōu)化結(jié)束的最優(yōu)解索力模型的全橋結(jié)構(gòu)位移及內(nèi)力狀態(tài)進(jìn)行比較,主梁彎矩及豎向位移、橋塔彎矩及縱向位移計(jì)算結(jié)果如圖7所示。結(jié)果顯示,結(jié)構(gòu)內(nèi)力及位移情況均得到較大改善。一次落架狀態(tài)主跨最大彎矩為3.81×105kN·m,經(jīng)過(guò)優(yōu)化后降低為4.58×104kN·m,降低了87.9 %;橋塔最大彎矩1.72×104kN·m優(yōu)化后降低為5.54×104kN·m,降低了67.8 %。主跨最大位移豎向位移由453.1 mm降低為44.7 mm,降低90.1 %;橋塔縱向位移由101.7 mm降低為3.7 mm,降低96.3 %。優(yōu)化后結(jié)構(gòu)線形符合設(shè)計(jì)要求,內(nèi)力有效降低,提高了整體性能。

    圖7 一次落架、較優(yōu)解、最優(yōu)解主跨及橋塔彎矩、位移對(duì)比

    3.3 算法改進(jìn)效果對(duì)比

    為對(duì)比改進(jìn)算法的優(yōu)越性,采用普通遺傳算法與上述結(jié)果進(jìn)行比較,如圖8所示,可見(jiàn)改進(jìn)算法在具有優(yōu)秀初始種群的情況下,適應(yīng)度函數(shù)在迭代初期更加符合要求,在整個(gè)求解過(guò)程中,改進(jìn)后算法收斂速度更快,計(jì)算效率較大提高,具有一定優(yōu)勢(shì)。

    圖8 改進(jìn)前后算法收斂效果對(duì)比

    4 結(jié)論

    本文針對(duì)斜塔斜拉橋的特殊受力情況,為使主梁和橋塔的線形符合設(shè)計(jì)要求,改善結(jié)構(gòu)整體力學(xué)性能,需求解合理成橋狀態(tài)索力,因此提出了一種基于遺傳算法的索力優(yōu)化方法,并引入精英保留策略和優(yōu)質(zhì)初始化種群策略對(duì)算法進(jìn)行改進(jìn)。首先建立Midas全橋梁?jiǎn)卧P停⑦\(yùn)用未知荷載系數(shù)法功能求解出一組符合主梁和橋塔彎矩、位移限定條件的合理索力較優(yōu)解,將該較優(yōu)解作為基準(zhǔn),生成了一組有較強(qiáng)競(jìng)爭(zhēng)優(yōu)勢(shì)的初始種群。之后利用Python編程語(yǔ)言對(duì)ABAQUS進(jìn)行二次開(kāi)發(fā),通過(guò)將初始種群分別寫(xiě)入索力文本文檔對(duì)有限元模型進(jìn)行更新,調(diào)用求解器進(jìn)行求解,并計(jì)算結(jié)構(gòu)彎曲應(yīng)變能。以遺傳算法為索力更新策略,根據(jù)種群內(nèi)個(gè)體的優(yōu)劣程度,進(jìn)行一系列選擇、交叉、變異操作,進(jìn)行多次迭代求解,最終得到索力最優(yōu)值。

    通過(guò)以某斜塔斜拉橋進(jìn)行實(shí)例計(jì)算,結(jié)果表明:與一次落架結(jié)構(gòu)狀態(tài)相比,整體結(jié)構(gòu)各項(xiàng)力學(xué)性能均有效改善,其中主跨最大彎矩降低了87.9 %,橋塔最大彎矩降低了67.8 %,主跨最大位移豎向位移由453.1 mm降低為44.7 mm,降低90.1 %;橋塔縱向位移由101.7 mm降低為3.7 mm,降低96.3 %。由此證明了算法的有效性。當(dāng)結(jié)構(gòu)較為復(fù)雜或種群基數(shù)較大時(shí),改進(jìn)算法較普通算法計(jì)算效率更高,在收斂速度上具有一定優(yōu)勢(shì)。

    猜你喜歡
    橋塔索力斜拉橋
    外懸掛附著式塔機(jī)對(duì)橋塔線形影響分析
    斜拉橋風(fēng)致振動(dòng)控制及其未來(lái)發(fā)展趨勢(shì)
    江蘇索力得新材料集團(tuán)有限公司
    矮塔斜拉橋彈塑性地震響應(yīng)分析
    (112+216+112)m部分斜拉橋設(shè)計(jì)
    超高異形橋塔及支架自振特性研究
    斜拉橋橋塔開(kāi)裂計(jì)算與分析
    基于索面布置的扇形索面斜拉橋橋塔穩(wěn)定性分析
    上地斜拉橋:天上滴落的水珠
    預(yù)應(yīng)力鋼絞線網(wǎng)加固混凝土橋梁的索力分布試驗(yàn)研究
    久久天躁狠狠躁夜夜2o2o| 三级国产精品欧美在线观看 | 国产精品99久久99久久久不卡| 白带黄色成豆腐渣| 色老头精品视频在线观看| 日本黄大片高清| 色在线成人网| 美女大奶头视频| 美女午夜性视频免费| 午夜福利免费观看在线| 亚洲五月婷婷丁香| 久久久久性生活片| 中出人妻视频一区二区| 亚洲人成网站高清观看| 岛国在线观看网站| 伊人久久大香线蕉亚洲五| 亚洲美女视频黄频| 午夜免费观看网址| 免费看美女性在线毛片视频| 美女cb高潮喷水在线观看 | 91字幕亚洲| 午夜免费观看网址| 国产精品爽爽va在线观看网站| 久久天躁狠狠躁夜夜2o2o| 在线看三级毛片| 欧美日韩亚洲国产一区二区在线观看| 日本一本二区三区精品| 欧美乱色亚洲激情| 在线观看66精品国产| 一区二区三区高清视频在线| 久久国产乱子伦精品免费另类| 色噜噜av男人的天堂激情| av天堂在线播放| 国产成人系列免费观看| 免费看a级黄色片| 成人三级黄色视频| 美女免费视频网站| 国产高清激情床上av| 亚洲国产欧洲综合997久久,| 国产美女午夜福利| 国产激情欧美一区二区| 免费一级毛片在线播放高清视频| 九九热线精品视视频播放| 日韩欧美免费精品| 一个人免费在线观看的高清视频| 91av网站免费观看| 欧美激情在线99| 99国产精品一区二区蜜桃av| 黑人巨大精品欧美一区二区mp4| 成人午夜高清在线视频| 五月伊人婷婷丁香| 亚洲av美国av| 久久中文字幕一级| 欧美成人一区二区免费高清观看 | www.精华液| 国产又黄又爽又无遮挡在线| 天堂网av新在线| 亚洲一区二区三区不卡视频| 搡老岳熟女国产| 一夜夜www| 91麻豆av在线| 国产精品国产高清国产av| 人人妻人人澡欧美一区二区| 精品国产美女av久久久久小说| 国产伦精品一区二区三区视频9 | 国产精品一及| 少妇人妻一区二区三区视频| 无遮挡黄片免费观看| 国产91精品成人一区二区三区| 久久久国产欧美日韩av| 在线视频色国产色| 无遮挡黄片免费观看| 一进一出好大好爽视频| 国产一区二区激情短视频| 啪啪无遮挡十八禁网站| 欧美国产日韩亚洲一区| 亚洲最大成人中文| 老汉色∧v一级毛片| 日韩欧美国产一区二区入口| 精品久久久久久,| 亚洲精品在线观看二区| 2021天堂中文幕一二区在线观| 精品久久久久久久久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 狂野欧美白嫩少妇大欣赏| 亚洲欧洲精品一区二区精品久久久| 精品电影一区二区在线| 精品国产美女av久久久久小说| a级毛片a级免费在线| 久久中文字幕人妻熟女| 国产伦人伦偷精品视频| 在线观看日韩欧美| 人妻久久中文字幕网| 国产 一区 欧美 日韩| 日日夜夜操网爽| 中文资源天堂在线| 老汉色∧v一级毛片| 三级男女做爰猛烈吃奶摸视频| 国产精品 欧美亚洲| 欧美黄色淫秽网站| 午夜福利成人在线免费观看| 国产高清三级在线| 中文字幕久久专区| 久久久久免费精品人妻一区二区| 国产精品99久久99久久久不卡| 色哟哟哟哟哟哟| 久久国产乱子伦精品免费另类| 一个人免费在线观看电影 | av在线蜜桃| 亚洲精品美女久久av网站| 黄色片一级片一级黄色片| 一级毛片女人18水好多| 成人高潮视频无遮挡免费网站| av福利片在线观看| 精品国产亚洲在线| 男女午夜视频在线观看| 亚洲精品美女久久av网站| 欧美成狂野欧美在线观看| 色av中文字幕| 99热这里只有精品一区 | 久久精品国产清高在天天线| 免费观看人在逋| 精品熟女少妇八av免费久了| 国产综合懂色| 小说图片视频综合网站| 久久99热这里只有精品18| 91久久精品国产一区二区成人 | 亚洲九九香蕉| 亚洲av熟女| 可以在线观看的亚洲视频| 欧美色欧美亚洲另类二区| 中文字幕高清在线视频| 成人国产一区最新在线观看| 国产精品久久久久久精品电影| 91久久精品国产一区二区成人 | 操出白浆在线播放| 九色国产91popny在线| 日本成人三级电影网站| 国产精品久久久人人做人人爽| 搡老妇女老女人老熟妇| 久久精品国产99精品国产亚洲性色| 欧美高清成人免费视频www| 成人特级黄色片久久久久久久| 男女午夜视频在线观看| 久久久久久国产a免费观看| 一级作爱视频免费观看| 一级毛片女人18水好多| 欧美日韩精品网址| 久久久国产欧美日韩av| 美女高潮喷水抽搐中文字幕| 中亚洲国语对白在线视频| 久久欧美精品欧美久久欧美| 国产野战对白在线观看| 久久婷婷人人爽人人干人人爱| 日韩 欧美 亚洲 中文字幕| 亚洲,欧美精品.| 亚洲狠狠婷婷综合久久图片| 日韩高清综合在线| 午夜精品在线福利| 国产精华一区二区三区| 狂野欧美激情性xxxx| 国产aⅴ精品一区二区三区波| 欧美高清成人免费视频www| 国产99白浆流出| 午夜福利欧美成人| 成人特级av手机在线观看| 国产淫片久久久久久久久 | 岛国在线免费视频观看| 亚洲男人的天堂狠狠| 国产1区2区3区精品| 欧美激情久久久久久爽电影| 日韩欧美国产一区二区入口| 一区福利在线观看| 波多野结衣高清作品| 国产真实乱freesex| 国产午夜精品论理片| 久久99热这里只有精品18| 国产真实乱freesex| 特级一级黄色大片| 亚洲精品色激情综合| 成年免费大片在线观看| 亚洲一区二区三区色噜噜| 最近视频中文字幕2019在线8| 99热这里只有精品一区 | 亚洲九九香蕉| 精品午夜福利视频在线观看一区| 国产探花在线观看一区二区| 欧美高清成人免费视频www| 香蕉久久夜色| 精品电影一区二区在线| a级毛片在线看网站| 99久久精品一区二区三区| 97人妻精品一区二区三区麻豆| 99久久精品国产亚洲精品| 99国产精品一区二区三区| 国产欧美日韩一区二区三| 国产麻豆成人av免费视频| 在线十欧美十亚洲十日本专区| 国产91精品成人一区二区三区| 老司机午夜十八禁免费视频| 国产精品国产高清国产av| 男女做爰动态图高潮gif福利片| 久久久久久国产a免费观看| 久久久久久久精品吃奶| x7x7x7水蜜桃| 成人精品一区二区免费| av中文乱码字幕在线| 欧美xxxx黑人xx丫x性爽| 免费一级毛片在线播放高清视频| 国产免费av片在线观看野外av| 国产精品美女特级片免费视频播放器 | 网址你懂的国产日韩在线| 嫩草影视91久久| 亚洲真实伦在线观看| 99热只有精品国产| 高清毛片免费观看视频网站| 国产成人一区二区三区免费视频网站| 日韩欧美国产在线观看| 超碰成人久久| 久久精品人妻少妇| 欧美黑人巨大hd| 美女被艹到高潮喷水动态| 三级国产精品欧美在线观看 | 国产激情久久老熟女| 免费一级毛片在线播放高清视频| 91老司机精品| 国产免费男女视频| 黄色 视频免费看| 少妇丰满av| 日本 欧美在线| 天堂影院成人在线观看| 色噜噜av男人的天堂激情| 久99久视频精品免费| 国产精品99久久99久久久不卡| 成人特级av手机在线观看| 国产高清videossex| 免费看a级黄色片| 禁无遮挡网站| 亚洲欧美精品综合一区二区三区| 日本免费a在线| 啦啦啦观看免费观看视频高清| 久久人人精品亚洲av| 亚洲国产精品成人综合色| 黄色片一级片一级黄色片| 精品久久久久久久久久久久久| 日本熟妇午夜| 男女下面进入的视频免费午夜| 日韩有码中文字幕| 亚洲欧美精品综合一区二区三区| 欧美日韩乱码在线| 欧美高清成人免费视频www| 在线免费观看的www视频| 久久精品国产99精品国产亚洲性色| 999精品在线视频| 欧美在线一区亚洲| 国产在线精品亚洲第一网站| 色哟哟哟哟哟哟| 亚洲18禁久久av| 99riav亚洲国产免费| 性欧美人与动物交配| 两性午夜刺激爽爽歪歪视频在线观看| 1024香蕉在线观看| 日韩欧美在线二视频| 国产男靠女视频免费网站| 男人的好看免费观看在线视频| 日韩欧美在线乱码| 国产精品香港三级国产av潘金莲| 日本免费a在线| 一级作爱视频免费观看| 国产av一区在线观看免费| 国产熟女xx| 欧美又色又爽又黄视频| www.www免费av| 在线十欧美十亚洲十日本专区| 90打野战视频偷拍视频| 国产激情欧美一区二区| 亚洲熟妇中文字幕五十中出| xxx96com| 久久天堂一区二区三区四区| 欧美国产日韩亚洲一区| 久久人人精品亚洲av| 欧美日韩一级在线毛片| 国产一区二区三区视频了| 国产一区二区在线av高清观看| 久久久久国产一级毛片高清牌| 黄色成人免费大全| 黄色女人牲交| 一本久久中文字幕| 色综合婷婷激情| 巨乳人妻的诱惑在线观看| bbb黄色大片| 亚洲成人免费电影在线观看| 在线观看午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 精品免费久久久久久久清纯| 色综合婷婷激情| 国产欧美日韩精品亚洲av| 精品久久久久久久久久久久久| 五月伊人婷婷丁香| 欧美一级毛片孕妇| av天堂在线播放| 精品熟女少妇八av免费久了| 国产高潮美女av| 黄片大片在线免费观看| 欧美日韩综合久久久久久 | 91老司机精品| 欧美黑人欧美精品刺激| 美女高潮的动态| 欧美日韩精品网址| 一级a爱片免费观看的视频| 国产精品精品国产色婷婷| 99久久精品国产亚洲精品| netflix在线观看网站| 我的老师免费观看完整版| 99热只有精品国产| 51午夜福利影视在线观看| 久久伊人香网站| 三级男女做爰猛烈吃奶摸视频| 变态另类丝袜制服| 成人三级做爰电影| 一区二区三区激情视频| 国产精品一及| 97超级碰碰碰精品色视频在线观看| 国产黄片美女视频| 亚洲黑人精品在线| 又黄又爽又免费观看的视频| 亚洲av熟女| 又黄又粗又硬又大视频| 日本免费a在线| 久久精品国产综合久久久| 免费电影在线观看免费观看| 国产视频内射| 99久久99久久久精品蜜桃| 老汉色∧v一级毛片| 深夜精品福利| 国模一区二区三区四区视频 | 性欧美人与动物交配| 国产精品98久久久久久宅男小说| 亚洲激情在线av| 两性夫妻黄色片| 黑人操中国人逼视频| 十八禁网站免费在线| 日本一本二区三区精品| 国产高清视频在线观看网站| 人妻丰满熟妇av一区二区三区| 宅男免费午夜| 99热这里只有精品一区 | 久久久久久久久久黄片| 免费搜索国产男女视频| 亚洲国产中文字幕在线视频| 国产精品98久久久久久宅男小说| 亚洲美女视频黄频| tocl精华| 精品一区二区三区视频在线观看免费| 男人和女人高潮做爰伦理| 久久国产精品影院| 男插女下体视频免费在线播放| 日本成人三级电影网站| 禁无遮挡网站| 悠悠久久av| 久9热在线精品视频| 国产99白浆流出| 成年人黄色毛片网站| 亚洲国产高清在线一区二区三| 欧美xxxx黑人xx丫x性爽| 熟妇人妻久久中文字幕3abv| 长腿黑丝高跟| 老司机午夜十八禁免费视频| 熟妇人妻久久中文字幕3abv| 色噜噜av男人的天堂激情| 桃红色精品国产亚洲av| 亚洲 欧美一区二区三区| 搡老岳熟女国产| e午夜精品久久久久久久| 午夜福利在线在线| 成人鲁丝片一二三区免费| 国产野战对白在线观看| 婷婷丁香在线五月| 婷婷六月久久综合丁香| 久久热在线av| 亚洲精品粉嫩美女一区| 9191精品国产免费久久| 精品免费久久久久久久清纯| 精品福利观看| 亚洲五月婷婷丁香| 午夜视频精品福利| 三级毛片av免费| 成人特级黄色片久久久久久久| 国产精品女同一区二区软件 | 精品日产1卡2卡| 一区二区三区激情视频| 久久香蕉国产精品| 亚洲av中文字字幕乱码综合| 叶爱在线成人免费视频播放| 在线十欧美十亚洲十日本专区| 欧美zozozo另类| 精品福利观看| 久久性视频一级片| 国产精品综合久久久久久久免费| 老司机午夜十八禁免费视频| 巨乳人妻的诱惑在线观看| 成人亚洲精品av一区二区| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久电影中文字幕| 国产精品,欧美在线| 亚洲第一欧美日韩一区二区三区| 日韩欧美国产一区二区入口| 欧美色视频一区免费| 欧美高清成人免费视频www| 亚洲精品一卡2卡三卡4卡5卡| 精品国产美女av久久久久小说| 99精品久久久久人妻精品| 欧美绝顶高潮抽搐喷水| 床上黄色一级片| 国产1区2区3区精品| 久久久精品欧美日韩精品| 国产伦在线观看视频一区| 两个人视频免费观看高清| 网址你懂的国产日韩在线| 亚洲欧美一区二区三区黑人| 小蜜桃在线观看免费完整版高清| 日本五十路高清| 国产精品久久视频播放| 99热这里只有是精品50| 搡老岳熟女国产| 黄片大片在线免费观看| 成人高潮视频无遮挡免费网站| 亚洲av日韩精品久久久久久密| 91在线观看av| 欧美zozozo另类| 在线a可以看的网站| 两性夫妻黄色片| 日本免费一区二区三区高清不卡| 99久久精品国产亚洲精品| 亚洲成人中文字幕在线播放| av黄色大香蕉| 长腿黑丝高跟| 国产亚洲精品一区二区www| 欧美成人性av电影在线观看| 女同久久另类99精品国产91| av女优亚洲男人天堂 | 夜夜夜夜夜久久久久| 97人妻精品一区二区三区麻豆| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产欧洲综合997久久,| 青草久久国产| 黄色成人免费大全| 男人和女人高潮做爰伦理| 91字幕亚洲| 国产成人欧美在线观看| 国产免费av片在线观看野外av| 亚洲在线观看片| 视频区欧美日本亚洲| 欧美绝顶高潮抽搐喷水| 99riav亚洲国产免费| 国产高清视频在线播放一区| 又爽又黄无遮挡网站| 久久精品国产亚洲av香蕉五月| 精品久久久久久久久久免费视频| 亚洲av电影在线进入| www日本在线高清视频| 18美女黄网站色大片免费观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲va日本ⅴa欧美va伊人久久| 日本 欧美在线| 亚洲一区二区三区色噜噜| 老司机午夜十八禁免费视频| 国产主播在线观看一区二区| 高潮久久久久久久久久久不卡| 亚洲狠狠婷婷综合久久图片| 国产黄片美女视频| av女优亚洲男人天堂 | 国产高清videossex| 亚洲精品美女久久av网站| 好看av亚洲va欧美ⅴa在| 男女视频在线观看网站免费| 国内精品美女久久久久久| 在线看三级毛片| 久久亚洲精品不卡| 日日夜夜操网爽| 少妇裸体淫交视频免费看高清| 变态另类成人亚洲欧美熟女| 国产美女午夜福利| 色噜噜av男人的天堂激情| 国产精品一区二区免费欧美| 欧美性猛交╳xxx乱大交人| 精品福利观看| 在线观看免费午夜福利视频| 国产一区在线观看成人免费| 夜夜躁狠狠躁天天躁| 国产精品自产拍在线观看55亚洲| 免费无遮挡裸体视频| 精品欧美国产一区二区三| 日本与韩国留学比较| 国产精品 欧美亚洲| 床上黄色一级片| 两性午夜刺激爽爽歪歪视频在线观看| 黄片大片在线免费观看| 人人妻人人澡欧美一区二区| 国产人伦9x9x在线观看| 亚洲av电影在线进入| 亚洲色图av天堂| 国产欧美日韩一区二区精品| 国产精品久久久av美女十八| 国产精品一区二区免费欧美| 精品国产乱码久久久久久男人| 久久久久国产精品人妻aⅴ院| 日韩成人在线观看一区二区三区| 成年免费大片在线观看| www.自偷自拍.com| 在线免费观看不下载黄p国产 | 日韩免费av在线播放| 香蕉国产在线看| 国内精品美女久久久久久| 给我免费播放毛片高清在线观看| 亚洲 欧美 日韩 在线 免费| 欧美成人一区二区免费高清观看 | 亚洲精品久久国产高清桃花| 久久精品91蜜桃| 99久久精品国产亚洲精品| 国产精品久久久人人做人人爽| 真人做人爱边吃奶动态| 久久亚洲真实| 日韩免费av在线播放| 欧美极品一区二区三区四区| 哪里可以看免费的av片| 久久久久国产精品人妻aⅴ院| 精品欧美国产一区二区三| 亚洲国产精品999在线| 久久久久久久久免费视频了| 制服丝袜大香蕉在线| 动漫黄色视频在线观看| 免费看日本二区| 亚洲18禁久久av| 欧美绝顶高潮抽搐喷水| 后天国语完整版免费观看| 97人妻精品一区二区三区麻豆| 精品久久久久久久久久免费视频| 性色av乱码一区二区三区2| 成年女人看的毛片在线观看| 成人鲁丝片一二三区免费| 97超视频在线观看视频| 国产伦在线观看视频一区| 久久性视频一级片| 91av网一区二区| 一级毛片女人18水好多| 亚洲欧美日韩卡通动漫| 日韩免费av在线播放| 日韩欧美在线二视频| 成人特级av手机在线观看| av福利片在线观看| 99久久精品国产亚洲精品| 美女黄网站色视频| 一个人看视频在线观看www免费 | 欧美大码av| 亚洲九九香蕉| 不卡一级毛片| 欧美极品一区二区三区四区| 成年女人永久免费观看视频| 国产精品一区二区精品视频观看| 非洲黑人性xxxx精品又粗又长| 久久天堂一区二区三区四区| 成人永久免费在线观看视频| 久久性视频一级片| av国产免费在线观看| 国产欧美日韩精品一区二区| 日本 欧美在线| 人人妻,人人澡人人爽秒播| 观看免费一级毛片| 亚洲欧美一区二区三区黑人| ponron亚洲| 久久伊人香网站| 国产精品自产拍在线观看55亚洲| 亚洲电影在线观看av| 99久久综合精品五月天人人| 19禁男女啪啪无遮挡网站| 国产精品1区2区在线观看.| 午夜激情欧美在线| 国产精品,欧美在线| 不卡av一区二区三区| 久99久视频精品免费| 国产精品久久久久久精品电影| 香蕉久久夜色| 麻豆国产av国片精品| 男女做爰动态图高潮gif福利片| netflix在线观看网站| 成年人黄色毛片网站| 一a级毛片在线观看| 中文字幕av在线有码专区| 免费看十八禁软件| 亚洲成人精品中文字幕电影| 看免费av毛片| 久久天躁狠狠躁夜夜2o2o| 2021天堂中文幕一二区在线观| 欧美一级a爱片免费观看看| 一个人观看的视频www高清免费观看 | 日韩三级视频一区二区三区| 国产黄片美女视频| 最新在线观看一区二区三区| 国产av不卡久久| 亚洲欧美精品综合一区二区三区| e午夜精品久久久久久久| 长腿黑丝高跟| 三级男女做爰猛烈吃奶摸视频| 女人被狂操c到高潮| 亚洲精品粉嫩美女一区| 国产视频一区二区在线看| 精品午夜福利视频在线观看一区| 久久精品综合一区二区三区| 亚洲熟女毛片儿| 特大巨黑吊av在线直播| 一个人免费在线观看的高清视频| 国产极品精品免费视频能看的| 熟女少妇亚洲综合色aaa.| 黄色视频,在线免费观看|