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

    混凝土斜拉橋的優(yōu)化設(shè)計

    2016-08-26 07:53:01彭琳琳中交二航局第四工程有限公司
    珠江水運 2016年13期
    關(guān)鍵詞:橋塔斜拉橋拉索

    ◎ 彭琳琳 中交二航局第四工程有限公司

    ?

    混凝土斜拉橋的優(yōu)化設(shè)計

    ◎ 彭琳琳 中交二航局第四工程有限公司

    斜拉橋設(shè)計包括大量設(shè)計變量和設(shè)計目標?;炷列崩瓨騼?yōu)化在此被歸結(jié)為多目標優(yōu)化問題,優(yōu)化目的是成本最低,偏移最小及應(yīng)力最小,開發(fā)數(shù)值方法以獲得此類結(jié)構(gòu)的優(yōu)化設(shè)計。結(jié)構(gòu)分析能說明所有相關(guān)的效應(yīng)(混凝土?xí)r間依存性效應(yīng),施工階段和幾何非線性效應(yīng))。本文通過對真實大小混凝土斜拉橋優(yōu)化的數(shù)值例子進行闡述,證明了該方法的特點和適用性。

    結(jié)構(gòu)優(yōu)化 混凝土斜拉橋 橫斷面設(shè)計變量

    混凝土斜拉橋優(yōu)化被歸結(jié)為多目標優(yōu)化問題,其目的是獲得最低成本,最小偏移和應(yīng)力,并尋求帕累托解決方案。通過凸標量函數(shù)的最小化發(fā)現(xiàn)了該解決方案。數(shù)值方法允許纜索安裝力、纜索調(diào)整力和斜拉索區(qū)域、橋面和橋塔橫截面的計算,以滿足整個結(jié)構(gòu)位移和應(yīng)力,無論是在架設(shè)期間或竣工。因此本文所提出的優(yōu)化算法可以被用來研究拓撲和幾何設(shè)計變量的效應(yīng)。

    在MATLAB開發(fā)的計算機程序被用于結(jié)構(gòu)分析、靈敏性分析和優(yōu)化。結(jié)構(gòu)分析采用一種有限元模型,其包括由于混凝土施工順序、幾何非線性和時間依存性效應(yīng)造成的負載歷史和幾何變化。施工過程和混凝土流變行為顯著影響了混凝土斜拉橋的應(yīng)力和變形。在分析中,同樣應(yīng)考慮處理纜索或大型靈活結(jié)構(gòu)時產(chǎn)生的幾何非線性情況。

    1.時間依存性效應(yīng)建模

    在這項研究中,根據(jù)歐洲規(guī)范2公式進行評估混凝土的蠕變、收縮和老化的時間依存賴性效應(yīng)。蠕變模型基于線性粘彈性及考慮老化影響。收縮應(yīng)變依存于時間,但與應(yīng)力無關(guān)。由作者在最近的研究工作中提出了有關(guān)時間依存性效應(yīng)建模的詳細考慮。

    1.1混凝土老化

    由于固化的結(jié)果,混凝土的強度和彈性系數(shù)隨時間增加。在早期,強度和彈性系數(shù)迅速增加,然后增加逐漸停滯,但不會完全停止。根據(jù)EN1992-1-1(2010),天數(shù)t,混凝土的彈性系數(shù)由下式給出(1)。

    其中,Ecm是混凝土彈性系數(shù)在28天內(nèi)的平均值,t是混凝土固化天數(shù),s是取決于水泥類型的系數(shù)。

    1.2混凝土蠕變

    在時刻t0施加單軸應(yīng)力σc的混凝土試件,其在時刻t的總應(yīng)變可以寫成相關(guān)應(yīng)力?ca(t,t0)和與無關(guān)應(yīng)力?cn(t)的總和,應(yīng)變?yōu)椋?)。

    其中J(t,t0)是蠕變功能;如果應(yīng)力小于45%的混凝土抗壓強度標準值,則疊加原理是有效的,并且蠕變應(yīng)變隨著所施加的應(yīng)力而線性變化。

    在施工階段和結(jié)構(gòu)使用壽命期間,斜拉橋的應(yīng)力不斷變化。根據(jù)變量的應(yīng)力以及使用疊加原理,方程(2)可改寫為(3)。

    已經(jīng)提出了幾種方法來解決這個方程,包括簡化方法,逐步數(shù)值積分以及蠕變函數(shù)的粗略估算。在這篇文章中,由Dirichlet級數(shù)(Bazant,1988)概略估算該蠕變函數(shù),從而得到(4)。

    其中,n是Dirichlet級數(shù)的項數(shù),系數(shù)aj通過使用最小二乘方法從曲線擬合得出。系數(shù)1/aj被稱為延遲時間,并被選擇來覆蓋用于計算蠕變系數(shù)的時間值范圍以。

    1.3混凝土收縮

    根據(jù)EN1992-1-1(2010),在時刻t的總收縮應(yīng)變?cs(t),是自收縮(?ca)和干燥收縮(?cd)的總和。時刻t的干燥收縮被定義為(5)。

    其中,系數(shù)βds(t,ts)和kh取決于構(gòu)件理論尺寸和混凝土干燥收縮開始時刻。參數(shù)εcd(t)取決于環(huán)境相對濕度、水泥類型和混凝土抗壓強度。

    自收縮是混凝土早期硬化期間化學(xué)反應(yīng)的結(jié)果,它可以通過時刻t由以下公式表示(6)。

    其中,εca(∞)是自收縮應(yīng)力的長期值,βas(t)是自收縮隨時間演變的函數(shù)。

    1.4時間依存性效應(yīng)模擬

    在結(jié)構(gòu)分析中,通過產(chǎn)生相同位移場的等效節(jié)點力模擬時間依存效應(yīng),并將其作為時間依存性效應(yīng)及從中計算出實際變形狀態(tài)。每個時間間隔,使用有限元公式和相應(yīng)的蠕變值以及根據(jù)先前提出公式計算出的收縮應(yīng)變,將這些計算的應(yīng)力作為初始變形。于是,僅用應(yīng)力和機械原點變形之間的彈性本構(gòu)關(guān)系計算該應(yīng)力。

    1.5幾何非線性效應(yīng)

    斜拉橋中幾何非線性有三個主要來源:由于自身重量造成下垂的非線性軸向力-延伸關(guān)聯(lián)斜拉索;根據(jù)聯(lián)合彎曲和軸向力,非線性軸向力和彎矩-變形關(guān)聯(lián)橋塔和橋面;以及大位移引起的幾何變化。在這篇文章中,通過二階彈性分析方法考慮幾何非線性效應(yīng)。

    考慮斜拉索中幾何非線性的一種廣泛使用方法是用于考慮能描述纜索懸索效應(yīng)并具有等效彈性系數(shù)(Ernst 1965年)的等效直弦構(gòu)件。

    纜索等效彈性系數(shù)值由下式給出(7)。

    其中,Eeq是等效的纜索彈性系數(shù),E是等效的纜索材料彈性系數(shù),y是纜索材料的比重,L是弦的長度,α是索弦與水平方向之間的角度,σ是纜索的張力應(yīng)力。

    二階效應(yīng)被認為是采用等效側(cè)力法,也被稱為虛擬側(cè)向荷載法或P-A迭代方法,通常用于建筑結(jié)構(gòu)二階分析。

    1.6包括架設(shè)階段的結(jié)構(gòu)分析

    繼施工階段后,使用正向分析程序,進行架設(shè)階段的建模和分析,然后獲得整個施工的應(yīng)力和位移,使其能直接考慮時間依存性效應(yīng)。這適合于結(jié)構(gòu)優(yōu)化目的,因為在各階段結(jié)束時和要求優(yōu)化模塊之前,可同時使用所有有關(guān)應(yīng)力、位移及其靈敏度的信息。

    在斜拉橋施工中,采用最普遍和廣泛使用的方法——平衡懸臂施工法進行架設(shè)。施工一開始,要建造橋塔,并在橋塔的兩側(cè)開始懸臂施工。在隨后的階段,架設(shè)其他橋面節(jié)段和斜拉索,直到橋面合攏。

    使用一個有限元計算機程序進行結(jié)構(gòu)分析,該橋梁被構(gòu)建成一個二維構(gòu)架結(jié)構(gòu)。該模型不考慮橋面扭轉(zhuǎn)和結(jié)構(gòu)的三維特性。然而,由于該研究以架設(shè)和使用情況為中心,因而利用二維模型是足夠用的。

    2.優(yōu)化設(shè)計公式

    在斜拉橋的優(yōu)化設(shè)計中,由于設(shè)計變量數(shù)量高且設(shè)計目標非線性特征也會產(chǎn)生沖突,因此搜索空間很復(fù)雜。這里作為一個多目標優(yōu)化問題被提出,從中獲得帕累托最優(yōu)解向量。這意味著在不增加至少一個目標的情況下,不存在其他可行的能降低一個目標的向量。它涉及設(shè)計變量、設(shè)計目標和目標函數(shù)的定義。

    2.1設(shè)計變量

    考慮的設(shè)計變量是斜拉索區(qū)域、預(yù)應(yīng)力、橋面和橋塔的橫截面尺寸。設(shè)計變量用xi來表示,而全球設(shè)計變量向量是(8)。

    斜拉索由0.6英寸直徑的鋼索制成(15.7mm的公稱直徑和1.5cm2的橫截面積)。對于橋塔的橫截面,考慮使用矩形空心型鋼,對于橋面考慮使用三種橫截面類型:梁板、單室箱和三室箱。橋面和橋塔的橫截面設(shè)計變量對減重(或降低成本)有直接的影響。斜拉索區(qū)域和斜拉索索力在整個結(jié)構(gòu)的應(yīng)力分布中發(fā)揮極為重要的作用,因為斜拉索區(qū)域和斜拉索索力對橋面的梁型特性的范圍作出了定義。此外,斜拉索區(qū)域和斜拉索索力是調(diào)整橋梁幾何形狀和擾度控制的基礎(chǔ),否則可以只通過橋面的嚴格加勁來實現(xiàn),這與預(yù)期的減少材料相反。

    2.2設(shè)計目標

    斜拉橋的設(shè)計涉及到實現(xiàn)一些設(shè)計目標,以檢驗使用和強度標準。該目標應(yīng)以規(guī)范化的形式來配置。這些目標會因處于施工階段及恒載和活荷載下的整個橋梁位移和應(yīng)力極限值而產(chǎn)生。該設(shè)計還應(yīng)力求將該結(jié)構(gòu)的成本降至最低??紤]到這一點,第一個目標可表示為(9)。

    其中,C是該結(jié)構(gòu)的現(xiàn)時成本,C0為參比成本,這對應(yīng)于每個分析和優(yōu)化周期的初始成本。這確保了在每個周期,成本始終是優(yōu)化算法的主要目標之一。該結(jié)構(gòu)的成本被當做材料成本(混凝土、加強鋼和預(yù)應(yīng)力鋼)。這些材料的單價可以通過咨詢葡萄牙供應(yīng)商公司而獲得。第二組目標產(chǎn)生于限制橋面的垂直位移和橋塔的水平位移,以獲得要求的最終橋面縱剖面,并使橋塔彎曲擾度最小化(10)。

    其中,δ和δ0分別是位移值和控制下的位移限值。

    第三組的目標會因處于施工階段和恒載下的整個橋梁橋面位移和橋塔應(yīng)力極限值而產(chǎn)生。這些目標都與使用情況有關(guān)。

    根據(jù)歐洲規(guī)范2(EN1992-1-1 2010)的建議,混凝土的抗壓應(yīng)力限定為45%的混凝土抗壓強度(fck)的特征值,從而使混凝土保持在線性蠕變范圍內(nèi),并預(yù)防縱向開裂?;炷晾瓚?yīng)力限定為5%的混凝土典型軸向抗拉強度的分位數(shù)(fctk,0.05),以避免開裂,確保耐久性,得(11)、(12)。

    其中,σc是在混凝土構(gòu)件的作用應(yīng)力。使用抗拉和抗壓應(yīng)力及相應(yīng)的容許應(yīng)力時,要考慮相應(yīng)的信號。

    也應(yīng)對混凝土構(gòu)件的最大應(yīng)力進行檢查。該目標可以表示為(13)。

    其中,σc是在混凝土構(gòu)件的作用應(yīng)力,σallow是抗拉或抗壓中相應(yīng)的容許應(yīng)力。作用應(yīng)力通過作用軸向力(NEd)和彎矩(MEd)計算得出。容許值被確定為一個應(yīng)力,等同于橫截面組合的軸向力-彎矩設(shè)計阻力(NRd;MRd)。此值通過混凝土構(gòu)件相應(yīng)橫截面產(chǎn)生的無量綱交互圖獲得。對于這些計算結(jié)果,加強鋼區(qū)域不是一個設(shè)計變量,被認為是混凝土截面積的2%。這被用來作為混凝土構(gòu)件的加強鋼區(qū)域的平均值,表示普遍的實用價值。

    斜拉索中應(yīng)力的剩余目標為(14)、(15)。

    其中,σ和fpk分別是斜拉索中的作用應(yīng)力和預(yù)應(yīng)力鋼抗拉強度的特征值。如果作用應(yīng)力大于0.1fpk,則方程式(14)適用;如果作用應(yīng)力小于或等于0.1fpk,則方程式(15)適用。斜拉索拉伸應(yīng)力的0.1fpk下限值被考慮,以確保其結(jié)構(gòu)效能。方程式(14)中的k值等于架設(shè)期間的0.55,使用情況下的0.50和強度校核的0.74。

    2.3目標函數(shù)

    多目標優(yōu)化的目標是最大限度地減少該組所有目標的設(shè)計變量。這通過最小最大值優(yōu)化來實現(xiàn),得(16)、(17)。

    該問題是不連續(xù)、不可微的問題,因此不好解決。然而,如Simoes 和Templeman所述(1989),可以通過一個不受約束的凸標量最小函數(shù)來間接獲得該求解方法,結(jié)果是Kreisselmeier-Steinhauser函數(shù)(17)。這是連續(xù)、可微的問題,因此很容易解決。該函數(shù)只取決于一個控制參數(shù)p,它不能被減小。

    我們的目標函數(shù)gj(x)沒有一個明確的代數(shù)形式,只能通過一個特定的設(shè)計變量向量的結(jié)構(gòu)分析結(jié)果得出數(shù)值。采取的策略是通過明確的近似模型迭代序列手段,來解方程式(17)。一個明確的近似值可以通過線性項后截斷的所有目標函數(shù)gj(x)的泰勒級數(shù)展開式制定。該算法得(18)。

    其中,N和M分別是設(shè)計變量數(shù)量和目標數(shù)量。goj(x)和dg0j(x)/dxi是根據(jù)當前設(shè)計變量向量(x0)計算的目標及其靈敏度,其中制定了泰勒級數(shù)展開式。g0j(x)特定數(shù)值求解方程(18)僅構(gòu)成該問題完整求解方法的一個迭代法。使用MATLAB函數(shù)fmincon,進行目標函數(shù)的最小值計算,采用一系列二次問題,使處于約束范圍內(nèi)的各變量標量函數(shù)最小化。

    該算法與多個起點相關(guān)聯(lián),提供一個優(yōu)化的工程求解方法,通過重新排列斜拉索的支承剛度及橋面和橋塔之間的質(zhì)量分布變得合理。目標函數(shù)使過程從一個可行或不可行的求解方法開始。

    2.4敏感性分析

    考慮到源代碼的可用性和大量的目標(應(yīng)力和位移)與設(shè)計變量數(shù)量相對,則采用分析直接離散方法,進行靈敏度分析。相對于設(shè)計變量的位移靈敏度通過區(qū)分平衡方程獲得(19)~(21)。

    可在形式上被重寫,其中,Qvi是相對于ith設(shè)計變量的系統(tǒng)虛擬假載向量。位移靈敏度可表示為(22)。其中需要存儲剛度矩陣,給剛度矩陣和右側(cè)導(dǎo)數(shù)預(yù)編程序,以使位移導(dǎo)數(shù)可以通過N假載右側(cè)求解方法來計算。應(yīng)力靈敏度通過有限元應(yīng)力位移關(guān)系的鏈求導(dǎo)來確定(22)、(23)。

    在對導(dǎo)數(shù)表達式預(yù)編程序并請求這一階段的條件下,右側(cè)的第一項可在全球系統(tǒng)元素貢獻計算期間直接計算出來。由于位移導(dǎo)數(shù)已知,則很容易計算出來右側(cè)第二項。

    圖1 橋的幾何形狀和斜拉索的數(shù)目

    3.數(shù)值例子和結(jié)果

    3.1數(shù)值模型的說明

    數(shù)值模型涉及一個總跨徑為284m,中跨距為148m的對稱混凝土斜拉橋。橋塔的總高度為52m,橋面置于基礎(chǔ)以上15m的位置。圖1說明了作為示例橋的幾何結(jié)構(gòu)。

    使用平衡懸臂施工法的按工序劃分施工階段。有關(guān)完整的橋,橋面到橋塔連接僅僅是垂直位置的連接,但要保證結(jié)構(gòu)穩(wěn)定,在假設(shè)階段,該連接固定。

    在施工階段考慮采取的行動包括橋面自重和1.0kN/m2(這取決于人員及手動工具),0.5kN/m2(這取決于臨時設(shè)備),和懸臂邊緣上400kN的集中載荷(EN1991-1-6 2005)。

    在第1階段,在橋塔各側(cè)執(zhí)行第一橋面節(jié)段,并且安裝第一對斜拉索。在以后的每個階段,在橋塔各側(cè)上對稱地增加一個橋面節(jié)段和一個斜拉索。在第9階段,在橫向跨距末端處增加一個簡單的支座;而第10階段,在中跨處增加控制彎矩和軸向力的連接件,以模擬橋梁合攏。

    在兩個步驟(安裝施力和調(diào)整力)中計算斜拉索施力。使用第3節(jié)中所描述的優(yōu)化方法,來確定安裝施力,以確保在第9階段中橋面要求的幾何結(jié)構(gòu)(在橫向跨距中增加支座并在主跨距中實現(xiàn)橋梁合攏),并且施工階段期間該結(jié)構(gòu)的應(yīng)力要保持在容許限值內(nèi)。要做到這一點,設(shè)計目標gj包括在第9階段的位移及第9階段和施工階段期間的應(yīng)力。在第10階段,對索力進行了調(diào)整,以獲得在自重和和額外恒載為2.5kN/m2(橋面、人行道、安全防護欄和護欄)的條件下,橋梁施工結(jié)束時要求的橋面縱剖面。在第10階段后,對永久狀態(tài)且沒有時間效應(yīng)影響下的橋梁進行了三種荷載情況的分析,要考慮最大應(yīng)力的設(shè)計目標。荷載情況包括施加于整個橋面或僅施加于中跨或邊跨上的恒載加上4 kN/m2(道路交通)的活荷載。

    3.2設(shè)計變量

    在所分析例子的數(shù)值計算法中,橋的整體幾何形狀和橋面寬度(19.0m)為預(yù)先指定的恒定設(shè)計參數(shù)。為確定定徑設(shè)計變量,考慮了橋塔的兩個區(qū)域和橋面的三個區(qū)域。橋面橫截面的幾何形狀設(shè)計變量分析有三個選項:梁板、單室箱和三室箱截面。橋塔的設(shè)計變量,考慮了矩形空心型材。

    3.3優(yōu)化設(shè)計結(jié)果

    實施例1:梁板橫截面,例2:單室箱截面,實施例3:三室箱截面。索力從橋塔到中跨增加。對稱分布的斜拉索力相似。最大力發(fā)生在平衡與邊跨相比最大的中央跨度負荷支索上。支索力基本控制橋塔頂部的水平位移,從而控制橋塔彎曲變形和應(yīng)力。這是斜拉橋充分結(jié)構(gòu)行為的必要特征。實施例1、2、3的優(yōu)化結(jié)果是成本分別降低17.8%、15.7%和19.1%。這是由于在新設(shè)計中,減少了橋面、斜拉索和橋塔的橫截面。

    永久荷載作用下的全橋橋塔頂部的水平位移約為1cm。結(jié)果表明,在結(jié)構(gòu)收尾橋面的正應(yīng)力仍在使用條件允許的范圍之內(nèi),并且橋面幾乎完全受壓。總結(jié)了總成本和橋面成本、橋塔和拉索成本??梢哉f,實施例2的成本最小。橋面占橋成本最大,為總成本的68.8%-71.1%。拉索約占橋梁成本的12%,而橋塔占其余成本。

    4.結(jié)論

    綜上所述,混凝土斜拉橋的設(shè)計作為一個多目標優(yōu)化問題提出,涉及到拉索預(yù)應(yīng)力和設(shè)計變量即拉索、橋塔和橋面的橫截面尺寸。優(yōu)化程序為整個結(jié)構(gòu)的足夠位移值和應(yīng)力值提供了結(jié)構(gòu)效率解決方案,證明離散直接法靈敏度分析(每次迭代只需要進行一次結(jié)構(gòu)分析),在預(yù)測橋梁結(jié)構(gòu)性能上是有效的。

    鑒于相當數(shù)量的目標和設(shè)計變量,解決方案只由少量迭代后獲得,因此本文中使用的優(yōu)化算法是有效的。

    猜你喜歡
    橋塔斜拉橋拉索
    “拉索”精確測量最亮伽馬暴
    軍事文摘(2024年4期)2024-03-19 09:40:02
    外懸掛附著式塔機對橋塔線形影響分析
    建筑機械化(2022年7期)2022-07-29 13:25:34
    斜拉橋風致振動控制及其未來發(fā)展趨勢
    手緩解拉索優(yōu)化設(shè)計
    矮塔斜拉橋彈塑性地震響應(yīng)分析
    (112+216+112)m部分斜拉橋設(shè)計
    超高異形橋塔及支架自振特性研究
    斜拉橋橋塔開裂計算與分析
    基于索面布置的扇形索面斜拉橋橋塔穩(wěn)定性分析
    上地斜拉橋:天上滴落的水珠
    99精国产麻豆久久婷婷| 久久久久久人人人人人| 老司机在亚洲福利影院| 成年人午夜在线观看视频| 国产精品98久久久久久宅男小说| 激情在线观看视频在线高清 | av天堂久久9| 欧美久久黑人一区二区| 两个人看的免费小视频| 丁香六月欧美| 国产精品久久久久成人av| 99国产精品一区二区三区| 手机成人av网站| 在线观看免费视频日本深夜| 国产视频一区二区在线看| 美女视频免费永久观看网站| 婷婷成人精品国产| 99久久人妻综合| videos熟女内射| 国产淫语在线视频| 日韩视频在线欧美| 欧美老熟妇乱子伦牲交| 另类亚洲欧美激情| 99国产精品免费福利视频| 亚洲男人天堂网一区| 精品人妻在线不人妻| 亚洲欧美日韩另类电影网站| 国产亚洲精品一区二区www | 亚洲三区欧美一区| 悠悠久久av| 国产精品麻豆人妻色哟哟久久| 久久精品aⅴ一区二区三区四区| 一本一本久久a久久精品综合妖精| 久久久久精品人妻al黑| 叶爱在线成人免费视频播放| 中国美女看黄片| 国产熟女午夜一区二区三区| 欧美黄色淫秽网站| svipshipincom国产片| 久久青草综合色| 我要看黄色一级片免费的| 侵犯人妻中文字幕一二三四区| 免费在线观看影片大全网站| 99久久人妻综合| 亚洲伊人色综图| 国产一区有黄有色的免费视频| 久久毛片免费看一区二区三区| www.999成人在线观看| 欧美日韩亚洲综合一区二区三区_| 国产精品久久久久久人妻精品电影 | 国产精品久久久久久人妻精品电影 | 国产日韩一区二区三区精品不卡| 亚洲成人手机| 777久久人妻少妇嫩草av网站| 欧美亚洲日本最大视频资源| 精品亚洲成a人片在线观看| 在线永久观看黄色视频| 视频区欧美日本亚洲| 免费女性裸体啪啪无遮挡网站| 中亚洲国语对白在线视频| 午夜精品久久久久久毛片777| 成人国产av品久久久| 欧美精品啪啪一区二区三区| 亚洲三区欧美一区| 国产野战对白在线观看| 国产又色又爽无遮挡免费看| 精品人妻熟女毛片av久久网站| 香蕉丝袜av| 日韩大码丰满熟妇| 日本撒尿小便嘘嘘汇集6| 黑人猛操日本美女一级片| 免费一级毛片在线播放高清视频 | 天天躁狠狠躁夜夜躁狠狠躁| 宅男免费午夜| 黄色视频,在线免费观看| 午夜福利在线观看吧| 黄色片一级片一级黄色片| 亚洲精品国产一区二区精华液| 一级毛片精品| 激情在线观看视频在线高清 | 国产真人三级小视频在线观看| 中文字幕人妻丝袜一区二区| 亚洲av第一区精品v没综合| 成年人黄色毛片网站| 国产单亲对白刺激| 欧美久久黑人一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产在线观看jvid| 欧美午夜高清在线| 黄色视频,在线免费观看| 桃红色精品国产亚洲av| 国产97色在线日韩免费| 日本av手机在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 在线亚洲精品国产二区图片欧美| 97在线人人人人妻| 午夜福利一区二区在线看| 精品久久久久久久毛片微露脸| 视频在线观看一区二区三区| 老司机午夜十八禁免费视频| 久久久久精品国产欧美久久久| 色精品久久人妻99蜜桃| 日本av免费视频播放| 欧美中文综合在线视频| 最新美女视频免费是黄的| 50天的宝宝边吃奶边哭怎么回事| 熟女少妇亚洲综合色aaa.| 老司机福利观看| 青青草视频在线视频观看| 又黄又粗又硬又大视频| 麻豆成人av在线观看| 女人高潮潮喷娇喘18禁视频| 一区福利在线观看| 91精品三级在线观看| 一夜夜www| 亚洲国产av新网站| 成人18禁在线播放| 脱女人内裤的视频| 一级毛片精品| 国产亚洲av高清不卡| 欧美av亚洲av综合av国产av| 美女高潮喷水抽搐中文字幕| 国产精品av久久久久免费| 欧美日韩福利视频一区二区| av一本久久久久| 午夜福利影视在线免费观看| 男人舔女人的私密视频| 蜜桃在线观看..| 满18在线观看网站| 国产色视频综合| 伦理电影免费视频| 亚洲一区中文字幕在线| 极品少妇高潮喷水抽搐| 一本一本久久a久久精品综合妖精| 男女免费视频国产| 欧美日韩福利视频一区二区| 天堂8中文在线网| 老司机在亚洲福利影院| 国产高清激情床上av| 色精品久久人妻99蜜桃| 一级片'在线观看视频| 丰满少妇做爰视频| 手机成人av网站| 欧美亚洲日本最大视频资源| 免费观看av网站的网址| 天堂8中文在线网| 久久久久久久国产电影| 性色av乱码一区二区三区2| 国产激情久久老熟女| 国产亚洲一区二区精品| 国产伦人伦偷精品视频| 国产成人av激情在线播放| 好男人电影高清在线观看| 18禁美女被吸乳视频| 国产片内射在线| 12—13女人毛片做爰片一| 精品国产一区二区三区久久久樱花| 亚洲国产中文字幕在线视频| 黑人巨大精品欧美一区二区蜜桃| 国产精品 国内视频| 一级毛片女人18水好多| 国产淫语在线视频| 热99re8久久精品国产| 色婷婷久久久亚洲欧美| 老汉色av国产亚洲站长工具| 9191精品国产免费久久| 一个人免费看片子| 欧美亚洲 丝袜 人妻 在线| 在线观看免费高清a一片| 天天影视国产精品| 欧美成狂野欧美在线观看| 亚洲欧美日韩高清在线视频 | 国产av精品麻豆| 99热国产这里只有精品6| 日日摸夜夜添夜夜添小说| 午夜视频精品福利| 亚洲熟妇熟女久久| a级毛片在线看网站| 午夜久久久在线观看| 欧美国产精品va在线观看不卡| 亚洲美女黄片视频| 91成年电影在线观看| 视频区图区小说| 国产成人av激情在线播放| 免费在线观看影片大全网站| 黄色视频不卡| 成年人黄色毛片网站| 人人澡人人妻人| 欧美 日韩 精品 国产| 十八禁网站网址无遮挡| 两性夫妻黄色片| 久久ye,这里只有精品| 97人妻天天添夜夜摸| 亚洲国产看品久久| 狠狠狠狠99中文字幕| 久久久久久亚洲精品国产蜜桃av| 99国产精品99久久久久| 午夜精品国产一区二区电影| 99re6热这里在线精品视频| 黄片大片在线免费观看| 亚洲欧洲精品一区二区精品久久久| 欧美日韩av久久| 国产又爽黄色视频| 欧美激情高清一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产成人免费无遮挡视频| 美女视频免费永久观看网站| 中文字幕人妻丝袜一区二区| 久久性视频一级片| 欧美久久黑人一区二区| 亚洲熟女毛片儿| 久久天堂一区二区三区四区| 69av精品久久久久久 | 啦啦啦免费观看视频1| 久久精品国产亚洲av高清一级| 自拍欧美九色日韩亚洲蝌蚪91| 高清欧美精品videossex| 亚洲成人免费av在线播放| 天天添夜夜摸| 99热网站在线观看| 亚洲av片天天在线观看| 精品人妻1区二区| 极品教师在线免费播放| 国产成人免费观看mmmm| 搡老熟女国产l中国老女人| 丁香六月天网| 精品免费久久久久久久清纯 | 亚洲情色 制服丝袜| 建设人人有责人人尽责人人享有的| 久久天堂一区二区三区四区| 亚洲人成电影观看| 熟女少妇亚洲综合色aaa.| 欧美日韩av久久| 另类亚洲欧美激情| a在线观看视频网站| 久久影院123| 日韩中文字幕欧美一区二区| 久久精品aⅴ一区二区三区四区| 两性夫妻黄色片| 热re99久久国产66热| 久久久国产欧美日韩av| 51午夜福利影视在线观看| 精品午夜福利视频在线观看一区 | kizo精华| tocl精华| 日日夜夜操网爽| 丝瓜视频免费看黄片| 成年人免费黄色播放视频| 国产国语露脸激情在线看| 9191精品国产免费久久| 少妇精品久久久久久久| 亚洲欧美一区二区三区黑人| 久久久国产一区二区| 国产精品98久久久久久宅男小说| 久久毛片免费看一区二区三区| 高清毛片免费观看视频网站 | 无遮挡黄片免费观看| 18禁美女被吸乳视频| svipshipincom国产片| 一区二区三区激情视频| 亚洲精品美女久久av网站| 欧美日韩一级在线毛片| 午夜日韩欧美国产| 天天躁狠狠躁夜夜躁狠狠躁| 最近最新中文字幕大全免费视频| 欧美中文综合在线视频| 亚洲三区欧美一区| 欧美日韩中文字幕国产精品一区二区三区 | 好男人电影高清在线观看| 90打野战视频偷拍视频| 在线av久久热| 曰老女人黄片| 国产精品久久久人人做人人爽| 18在线观看网站| 色视频在线一区二区三区| 久久久久久人人人人人| 大陆偷拍与自拍| 国产亚洲欧美精品永久| 9热在线视频观看99| 国产不卡av网站在线观看| 久久九九热精品免费| 啦啦啦中文免费视频观看日本| 99久久人妻综合| 国产成人精品久久二区二区91| 多毛熟女@视频| 日韩一区二区三区影片| 女人久久www免费人成看片| 欧美黄色淫秽网站| 一级a爱视频在线免费观看| 99国产精品免费福利视频| 日韩一卡2卡3卡4卡2021年| 亚洲一区中文字幕在线| 人妻一区二区av| 午夜老司机福利片| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美亚洲二区| 欧美精品人与动牲交sv欧美| 午夜免费鲁丝| 亚洲美女黄片视频| 色精品久久人妻99蜜桃| 国产av国产精品国产| 欧美精品高潮呻吟av久久| 日韩人妻精品一区2区三区| 国产精品一区二区免费欧美| 嫩草影视91久久| 男女下面插进去视频免费观看| 国产aⅴ精品一区二区三区波| 老熟妇仑乱视频hdxx| 99精国产麻豆久久婷婷| 国产男女内射视频| 999久久久精品免费观看国产| 美女高潮到喷水免费观看| 亚洲成国产人片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 一二三四社区在线视频社区8| 99热网站在线观看| 黄色视频在线播放观看不卡| 欧美黑人欧美精品刺激| 99久久精品国产亚洲精品| 97人妻天天添夜夜摸| 欧美日本中文国产一区发布| 夜夜夜夜夜久久久久| 国产精品99久久99久久久不卡| 国产精品免费一区二区三区在线 | 高清黄色对白视频在线免费看| 精品第一国产精品| 午夜福利视频在线观看免费| av网站免费在线观看视频| 高清在线国产一区| 久久久久久久大尺度免费视频| 中文字幕另类日韩欧美亚洲嫩草| www.熟女人妻精品国产| 纯流量卡能插随身wifi吗| 亚洲av第一区精品v没综合| 久久热在线av| 乱人伦中国视频| 久久亚洲真实| 久久青草综合色| 精品一区二区三卡| 国产精品偷伦视频观看了| 91九色精品人成在线观看| 他把我摸到了高潮在线观看 | www.自偷自拍.com| cao死你这个sao货| 日日夜夜操网爽| 亚洲精品久久成人aⅴ小说| 一区在线观看完整版| 欧美乱码精品一区二区三区| 丝袜喷水一区| 91国产中文字幕| 色综合婷婷激情| 久久久国产欧美日韩av| 国产成人影院久久av| 欧美大码av| 亚洲伊人色综图| 久久精品国产亚洲av高清一级| 久久热在线av| 国产在线一区二区三区精| 在线看a的网站| 久久ye,这里只有精品| 99久久人妻综合| 中亚洲国语对白在线视频| 午夜福利视频精品| 两个人免费观看高清视频| 麻豆国产av国片精品| 欧美性长视频在线观看| 十八禁网站免费在线| 首页视频小说图片口味搜索| 亚洲成人免费电影在线观看| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 国产男靠女视频免费网站| 99riav亚洲国产免费| 欧美乱码精品一区二区三区| kizo精华| 国产成人欧美在线观看 | 成人国产一区最新在线观看| a级毛片在线看网站| 亚洲中文日韩欧美视频| 80岁老熟妇乱子伦牲交| 视频区图区小说| 九色亚洲精品在线播放| 国产男女超爽视频在线观看| 男女高潮啪啪啪动态图| 国产深夜福利视频在线观看| 久久久久久久久久久久大奶| 91精品国产国语对白视频| 少妇被粗大的猛进出69影院| 一区二区av电影网| 黄色视频在线播放观看不卡| 成人免费观看视频高清| 男女高潮啪啪啪动态图| 国产男靠女视频免费网站| 一区二区三区精品91| 99国产精品一区二区三区| 国产一区二区三区综合在线观看| 欧美日韩亚洲国产一区二区在线观看 | 美女扒开内裤让男人捅视频| 国产亚洲一区二区精品| 天堂中文最新版在线下载| 亚洲欧美日韩高清在线视频 | 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| 黄片大片在线免费观看| 亚洲专区字幕在线| 国产在线视频一区二区| 亚洲欧洲日产国产| 亚洲av成人不卡在线观看播放网| 视频在线观看一区二区三区| 精品久久久精品久久久| 捣出白浆h1v1| 正在播放国产对白刺激| 国产日韩欧美亚洲二区| 视频在线观看一区二区三区| 久久精品亚洲精品国产色婷小说| 天堂动漫精品| 女人精品久久久久毛片| 国产成人欧美| 日韩人妻精品一区2区三区| 丰满少妇做爰视频| 99国产精品免费福利视频| 午夜成年电影在线免费观看| 国产深夜福利视频在线观看| 成人手机av| 一边摸一边抽搐一进一出视频| 久久精品人人爽人人爽视色| 久久久久久久久免费视频了| 国产男靠女视频免费网站| 久久精品亚洲av国产电影网| 免费在线观看视频国产中文字幕亚洲| 国产免费av片在线观看野外av| 99久久精品国产亚洲精品| 韩国精品一区二区三区| 亚洲精华国产精华精| 久久久久网色| 亚洲欧美一区二区三区黑人| netflix在线观看网站| 99精品欧美一区二区三区四区| 老司机午夜福利在线观看视频 | 亚洲av电影在线进入| 久久久久久久久免费视频了| 久久精品熟女亚洲av麻豆精品| 一进一出抽搐动态| 侵犯人妻中文字幕一二三四区| 18禁美女被吸乳视频| 视频在线观看一区二区三区| 在线观看舔阴道视频| 国产有黄有色有爽视频| 国产精品美女特级片免费视频播放器 | 国产片内射在线| 丁香六月天网| 最黄视频免费看| 国产亚洲一区二区精品| 欧美精品人与动牲交sv欧美| 两性午夜刺激爽爽歪歪视频在线观看 | tocl精华| 国产亚洲av高清不卡| 国产成人精品无人区| 国产亚洲精品一区二区www | 亚洲欧美一区二区三区黑人| 9191精品国产免费久久| 国产在线视频一区二区| 国产日韩一区二区三区精品不卡| avwww免费| 一本色道久久久久久精品综合| 亚洲精品国产区一区二| 久久精品91无色码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲黑人精品在线| 我要看黄色一级片免费的| 日韩视频一区二区在线观看| 国产av国产精品国产| 国产精品二区激情视频| 最新美女视频免费是黄的| 午夜福利一区二区在线看| 亚洲综合色网址| 国产成人系列免费观看| 国产单亲对白刺激| 人人妻,人人澡人人爽秒播| 黄色片一级片一级黄色片| 国产欧美日韩一区二区三区在线| 香蕉国产在线看| 精品国产超薄肉色丝袜足j| 免费看十八禁软件| 一二三四在线观看免费中文在| 69av精品久久久久久 | 亚洲av日韩在线播放| 丰满少妇做爰视频| 成人免费观看视频高清| 高潮久久久久久久久久久不卡| h视频一区二区三区| 成年女人毛片免费观看观看9 | 亚洲欧美激情在线| 色精品久久人妻99蜜桃| 久久国产精品影院| 电影成人av| 欧美日韩精品网址| 交换朋友夫妻互换小说| 亚洲自偷自拍图片 自拍| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美网| 国产免费现黄频在线看| 国产1区2区3区精品| 视频区欧美日本亚洲| 91精品国产国语对白视频| 黄色视频不卡| 国产av国产精品国产| 亚洲精品粉嫩美女一区| 午夜免费鲁丝| 成年人免费黄色播放视频| 啦啦啦在线免费观看视频4| videosex国产| 一本大道久久a久久精品| 亚洲欧美一区二区三区黑人| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美免费精品| 亚洲欧美一区二区三区久久| 美女高潮到喷水免费观看| 91精品国产国语对白视频| 大香蕉久久网| 国产av国产精品国产| 女性生殖器流出的白浆| 男男h啪啪无遮挡| 亚洲精品粉嫩美女一区| 怎么达到女性高潮| 十八禁网站免费在线| 国产精品一区二区免费欧美| 男女床上黄色一级片免费看| 最近最新中文字幕大全电影3 | 久久人人爽av亚洲精品天堂| 91精品三级在线观看| 91老司机精品| 国产不卡av网站在线观看| 国产精品一区二区在线不卡| 欧美变态另类bdsm刘玥| www.精华液| av视频免费观看在线观看| 国产日韩一区二区三区精品不卡| 美女高潮喷水抽搐中文字幕| 50天的宝宝边吃奶边哭怎么回事| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 国产精品九九99| 成人精品一区二区免费| 老司机影院毛片| 又黄又粗又硬又大视频| 男女下面插进去视频免费观看| 9热在线视频观看99| 这个男人来自地球电影免费观看| 国产视频一区二区在线看| 一区二区三区精品91| 美女视频免费永久观看网站| av福利片在线| 国产精品影院久久| 国产成人欧美在线观看 | 亚洲成人免费av在线播放| 中文字幕色久视频| 国产主播在线观看一区二区| 国产高清videossex| 欧美黑人欧美精品刺激| 亚洲精品av麻豆狂野| 国产精品一区二区在线观看99| 欧美乱码精品一区二区三区| 久久人妻熟女aⅴ| 色婷婷久久久亚洲欧美| 久久免费观看电影| 成人特级黄色片久久久久久久 | xxxhd国产人妻xxx| 岛国毛片在线播放| 精品国产国语对白av| 国产有黄有色有爽视频| 老熟妇乱子伦视频在线观看| 亚洲精华国产精华精| 亚洲三区欧美一区| 国产熟女午夜一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲免费av在线视频| 国产在线观看jvid| 亚洲精品乱久久久久久| 国产成人系列免费观看| 91麻豆av在线| 日韩欧美一区视频在线观看| 久久ye,这里只有精品| 天堂动漫精品| 亚洲av成人一区二区三| 日本a在线网址| 一级a爱视频在线免费观看| 欧美成狂野欧美在线观看| 大型黄色视频在线免费观看| 久久久久久亚洲精品国产蜜桃av| 国产真人三级小视频在线观看| 国产男女内射视频| 人人妻人人爽人人添夜夜欢视频| 国产1区2区3区精品| 久久人妻福利社区极品人妻图片| 国产日韩欧美视频二区| 69精品国产乱码久久久| 极品少妇高潮喷水抽搐| 人妻久久中文字幕网| 在线永久观看黄色视频| 脱女人内裤的视频| 美女高潮到喷水免费观看| 女人爽到高潮嗷嗷叫在线视频| 老鸭窝网址在线观看| 黄网站色视频无遮挡免费观看| 99香蕉大伊视频| 精品亚洲成国产av| 成人三级做爰电影| 美女高潮喷水抽搐中文字幕| 午夜福利视频精品| 精品免费久久久久久久清纯 | 中文亚洲av片在线观看爽 | 大片电影免费在线观看免费| 另类亚洲欧美激情| 中文字幕av电影在线播放|