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

    鋼筋混凝土結(jié)構(gòu)多尺度建模界面連接方法研究

    2022-08-18 08:32:22吳怡凡沈建金世杰夏晉
    低溫建筑技術(shù) 2022年7期
    關(guān)鍵詞:軸力彎矩尺度

    吳怡凡, 沈建, 金世杰, 夏晉

    (1.浙江大學(xué)結(jié)構(gòu)工程研究所,杭州 310058;2.浙江大學(xué)建筑設(shè)計(jì)研究院有限公司,杭州 310028)

    0 引言

    隨著有限元技術(shù)的不斷發(fā)展,在大型土木工程結(jié)構(gòu)的模擬分析中,為了滿(mǎn)足整體結(jié)構(gòu)的力學(xué)響應(yīng),并反應(yīng)關(guān)鍵節(jié)點(diǎn)的塑性發(fā)展,多尺度建模成為了一種常用策略,即在整體結(jié)構(gòu)的非關(guān)鍵部位使用計(jì)算效率較高的宏觀單元,而在需要精細(xì)分析的關(guān)鍵節(jié)點(diǎn)部位使用精度較高的微觀單元,并通過(guò)一定的連接方式將兩種類(lèi)型的單元耦合起來(lái),以達(dá)到既兼顧計(jì)算效率又保證精度的目的[1-3]。多點(diǎn)約束法是目前土木工程領(lǐng)域進(jìn)行多尺度界面連接最常用的一種方法,它的核心思想是通過(guò)建立不同尺度界面節(jié)點(diǎn)自由度間的協(xié)調(diào)方程使兩界面達(dá)到有效的耦合[4]。在過(guò)去的研究中,學(xué)者們往往采用三種理論來(lái)構(gòu)建多點(diǎn)約束方程,分別為:位移協(xié)調(diào)、能量協(xié)調(diào)以及位移協(xié)調(diào)與能量協(xié)調(diào)相結(jié)合。林旭川等運(yùn)用平截面假定列出位移協(xié)調(diào)方程,采用MSC.MARC模擬了鋼筋混凝土框架多尺度模型在地震下的受力情況,并通過(guò)試驗(yàn)驗(yàn)證[5]。Yue和Wu等[6,7]利用做功相等能量協(xié)調(diào)原理,分別通過(guò)ABAQUS和ANSYS建立了鋼筋混凝土框架多尺度模型,并與實(shí)驗(yàn)比較,結(jié)果表明該方法可以準(zhǔn)確模擬結(jié)構(gòu)局部破壞。但在鋼筋混凝土結(jié)構(gòu)的數(shù)值模擬分析中,由于材料非線性、幾何非線性等因素,基于單一理論的多點(diǎn)約束法存在一定的局限性[8-10]?;谖灰茀f(xié)調(diào)與能量協(xié)調(diào)相結(jié)合的多點(diǎn)約束法可以避開(kāi)兩種單一理論的缺點(diǎn)得到更好的模擬結(jié)果,但需要插入數(shù)量龐大的約束方程,建模前處理繁瑣復(fù)雜[11]。

    文中提出了一種基于連續(xù)分布耦合理論的多點(diǎn)約束法以實(shí)現(xiàn)不同尺度單元間的連接,針對(duì)鋼筋混凝土節(jié)點(diǎn),使用ABAQUS比較分析了連續(xù)分布耦合與位移協(xié)調(diào)、能量協(xié)調(diào)、位移協(xié)調(diào)與能量協(xié)調(diào)相結(jié)合4種多尺度界面連接方法的優(yōu)劣性。

    1 多尺度界面連接方式分析

    文中以梁-實(shí)體單元的連接為例介紹不同理論實(shí)現(xiàn)多尺度界面連接的基本原理和方法,梁-實(shí)體單元連接如圖1所示。

    圖1 梁-實(shí)體單元連接

    1.1 基于位移協(xié)調(diào)的界面連接法

    連接界面上實(shí)體單元共有n個(gè)節(jié)點(diǎn),每個(gè)節(jié)點(diǎn)有x、y、z3個(gè)自由度;梁?jiǎn)卧c實(shí)體單元通過(guò)節(jié)點(diǎn)B連接,節(jié)點(diǎn)B有xB、yB、zB、θxB、θyB、θzB6個(gè)自由度。根據(jù)梁、實(shí)體單元可能產(chǎn)生的4種位移模式,可知界面節(jié)點(diǎn)間產(chǎn)生的變形協(xié)調(diào)如圖2所示。將精細(xì)實(shí)體單元節(jié)點(diǎn)的3個(gè)自由度分別用宏觀梁?jiǎn)卧?jié)點(diǎn)的6個(gè)自由度表示,一共形成3n個(gè)位移協(xié)調(diào)方程如式(1)。

    圖2 梁-實(shí)體單元多尺度界面變形協(xié)調(diào)

    式中,bxi、byi為實(shí)體單元節(jié)點(diǎn)i在X、Y坐標(biāo)軸上的坐標(biāo)。

    1.2 基于能量協(xié)調(diào)的界面連接法

    在軸力、彎矩、剪力、扭矩的作用下界面節(jié)點(diǎn)變形示意如圖3所示,根據(jù)虛功原理可列出做功相等方程如式(2)。

    圖3 梁-實(shí)體單元多尺度界面受力變形

    式中,F(xiàn)j為作用在宏觀單元B節(jié)點(diǎn)方向的力;σi1,Fj為精細(xì)單元i節(jié)點(diǎn)處由Fj引起的軸向應(yīng)力;τi2,Fj、τi3,Fj為精細(xì)單元i節(jié)點(diǎn)處由Fj引起的切向應(yīng)力;u1S、u2S、u3S為精細(xì)單元i節(jié)點(diǎn)處軸向及切向的位移。

    結(jié)合材料力學(xué)公式,并用形函數(shù)來(lái)表示節(jié)點(diǎn)位移X={N}{X},Y={N}{Y},Z={N}{Z},代入四節(jié)點(diǎn)四邊形單元形函數(shù)對(duì)積分方程求解,可得到約束方程如式(3)。

    式中,Cxi、Cyi、Czi、Cθxi、Cθyi為與實(shí)體單元截面尺寸及連接界面點(diǎn)數(shù)個(gè)數(shù)相關(guān)的系數(shù)。

    1.3 基于位移協(xié)調(diào)與能量協(xié)調(diào)相結(jié)合的界面連接法

    基于位移協(xié)調(diào)的多點(diǎn)約束法在切向存在過(guò)約束的缺陷,而基于能量協(xié)調(diào)的多點(diǎn)約束法在轉(zhuǎn)角方向存在材料非線性和幾何非線性的缺陷。

    因此將兩種方法在不同自由度上進(jìn)行方程聯(lián)立,在切向上使用基于能量協(xié)調(diào)的多點(diǎn)約束方程,而在軸向和轉(zhuǎn)角方向上使用基于位移協(xié)調(diào)的多點(diǎn)約束方程,聯(lián)立得到的方程組如式(4)。

    1.4 基于連續(xù)分布耦合的界面連接法

    連續(xù)分布耦合示意如圖4所示,連續(xù)分布耦合允許約束面上各節(jié)點(diǎn)發(fā)生相對(duì)移動(dòng),通過(guò)對(duì)約束面上各節(jié)點(diǎn)的運(yùn)動(dòng)進(jìn)行加權(quán)處理,使在此區(qū)域上受到的合力和合力距與施加在參考點(diǎn)上的力和力矩相等效。因此,使用連續(xù)分布耦合理論建立的多尺度模型在理論上可以避免位移協(xié)調(diào)理論產(chǎn)生的切向過(guò)約束。ABAQUS連續(xù)分布耦合有4種加權(quán)系數(shù)計(jì)算方式,包括一致、線性、二次方和三次方,計(jì)算公式如式(5)~式(8)。經(jīng)過(guò)反復(fù)試算文中采用一致加權(quán)系數(shù)計(jì)算方式對(duì)界面節(jié)點(diǎn)進(jìn)行加權(quán)處理。

    圖4 連續(xù)分布耦合

    式中,ωi為界面節(jié)點(diǎn)i的加權(quán)系數(shù);ri為界面節(jié)點(diǎn)i到參考點(diǎn)的徑向距離;ω0為最遠(yuǎn)點(diǎn)的徑向距離。

    2 不同多尺度界面連接法算例分析

    分析模型如圖5所示,構(gòu)件高2m,截面為邊長(zhǎng)0.4m的正方形,多尺度界面距離底部1m。圖5(a)為精細(xì)化模型,圖5(b)~圖5(e)分別為基于位移協(xié)調(diào)、能量協(xié)調(diào)、位移協(xié)調(diào)與能量協(xié)調(diào)相結(jié)合、連續(xù)分布耦合理論建立的多尺度模型。分別命名為位移法模型、能量法模型、結(jié)合法模型和DC法模型。

    圖5 鋼筋混凝土節(jié)點(diǎn)梁-實(shí)體單元多尺度連接模型

    對(duì)于實(shí)體單元部分鋼筋混凝土采用分離式建模處理,其建模原理是將混凝土和鋼筋分別建模,混凝土采用C3D8R單元,鋼筋采用T3D2單元,并利用嵌固作用將鋼筋嵌入到混凝土中,實(shí)現(xiàn)鋼筋與混凝土的良好結(jié)合。對(duì)于梁?jiǎn)卧糠值匿摻罨炷敛捎谜w式模型,其建模原理是在同一位置建立兩個(gè)共節(jié)點(diǎn)B31梁?jiǎn)卧?,梁?jiǎn)卧?為混凝土,梁?jiǎn)卧?用箱形截面等效模擬鋼筋,其建模方式如圖6所示。混凝土采用塑性損傷本構(gòu)模型,鋼筋采用無(wú)屈服點(diǎn)的彈塑性加工硬化兩折線模型。模型網(wǎng)格尺寸為50mm。

    圖6 箱形截面等效配筋

    2.1 軸力工況分析

    在分析模型的底部設(shè)置固定約束,在頂部設(shè)置加載點(diǎn)通過(guò)施加位移控制加載,豎向位移從0逐漸增長(zhǎng)到10mm。各模型在塑性階段混凝土單元的應(yīng)力云圖如圖7所示,各模型在破壞階段混凝土受壓損傷云圖如圖8所示。通過(guò)觀察加載過(guò)程中應(yīng)力云圖的變化可以發(fā)現(xiàn),在軸力工況下,隨著豎向位移的增大,各模型的單元應(yīng)力狀態(tài)逐漸由彈性進(jìn)入塑性,其中構(gòu)件的兩端和中部最先進(jìn)入塑性,混凝土在達(dá)到極限抗壓強(qiáng)度后發(fā)生破壞,應(yīng)力下降。在加載過(guò)程中,構(gòu)件的中部逐漸形成損傷累積區(qū)并最終發(fā)生破壞。對(duì)比發(fā)現(xiàn),除了位移法模型,其余各多尺度模型均能較好地反應(yīng)構(gòu)件在彈塑性階段的應(yīng)力和損傷分布,印證了基于位移協(xié)調(diào)理論構(gòu)建的多尺度模型在僅受軸力作用時(shí),多尺度界面存在側(cè)向過(guò)約束。

    圖7 軸力工況下塑性階段各模型混凝土壓應(yīng)力云圖

    圖8 軸力工況下各模型混凝土受壓損傷累積云圖

    提取各模型在加載過(guò)程中加載點(diǎn)的位移和對(duì)應(yīng)的軸力,作出頂點(diǎn)位移-軸力曲線如圖9所示。5組模型加載結(jié)束時(shí)頂點(diǎn)軸力分別為1510、1550、1610、1500、1530kN。相比精細(xì)化模型,各多尺度模型相對(duì)誤差分別為2.65%、5.88%、0.2%、1.12%。

    圖9 頂點(diǎn)位移-軸力曲線

    按照象限編號(hào)法則,由內(nèi)到外,在各模型多尺度界面上選取16個(gè)特征單元,如圖10所示。當(dāng)頂點(diǎn)豎向位移加載到2.15mm時(shí)構(gòu)件處于塑性階段,分別提取各模型16個(gè)特征單元的Min.principle應(yīng)力,統(tǒng)計(jì)精細(xì)化模型與多尺度模型在特征單元處的應(yīng)力值并計(jì)算相對(duì)誤差。觀察發(fā)現(xiàn),平均相對(duì)誤差從大到小分別為位移法模型10.1%、能量法模型1.8%、結(jié)合法模型1.6%、DC法模型4.41%。位移法模型的13、14號(hào)特征單元相對(duì)誤差最大,達(dá)到32.1%。能量法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到2.5%。結(jié)合法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到2.1%。DC法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到5.1%。在軸力工況下,位移法模型的相對(duì)誤差率高于其他多尺度模型,印證了位移協(xié)調(diào)法在僅受軸力作用下存在側(cè)向過(guò)約束局限性的分析。DC法模型的平均相對(duì)誤差為4.4%,控制在5%之內(nèi),印證了使用連續(xù)分布耦合約束可以有效改善側(cè)向過(guò)約束降低誤差率。

    圖10 多尺度界面特征單元編號(hào)

    2.2 彎矩工況分析

    在分析模型的底部設(shè)置固定約束,在頂部設(shè)置加載點(diǎn)通過(guò)施加繞X軸的彎矩荷載控制加載,彎矩從0逐漸增長(zhǎng)到40kN·m。

    各模型在塑性階段混凝土單元的應(yīng)力云圖如圖11所示,各模型在破壞階段混凝土受拉損傷云圖如圖12所示。通過(guò)觀察加載過(guò)程中應(yīng)力云圖的變化可以發(fā)現(xiàn),混凝土最大拉、壓應(yīng)力分別沿受拉、受壓側(cè)邊緣分布。隨著不斷加載,構(gòu)件在受拉側(cè),沿柱身不斷形成損傷累積區(qū)并最終發(fā)生破壞。對(duì)比各模型可以發(fā)現(xiàn),能量法模型相比精細(xì)化模型應(yīng)力分布和損傷累積差別都較大。位移法模型和DC發(fā)模型與精細(xì)化模型在應(yīng)力分布上吻合度比較高,但在多尺度界面附近的損傷累積有所差別。結(jié)合法模型與精細(xì)化模型吻合度極高,能很好地反應(yīng)構(gòu)件的應(yīng)力分布及損傷累積。

    圖11 彎矩工況下塑性階段各模型混凝土軸向應(yīng)力云圖

    圖12 彎矩工況下各模型混凝土受拉損傷累積云圖

    提取各模型在加載過(guò)程中加載點(diǎn)繞X軸的彎矩和對(duì)應(yīng)繞X軸的轉(zhuǎn)角,作出頂點(diǎn)轉(zhuǎn)角-彎矩曲線如圖13所示。5組模型加載結(jié)束時(shí)頂點(diǎn)轉(zhuǎn)角分別為0.00718、0.00923、0.01768、0.00676、0.00493。相比精細(xì)化模型,各多尺度模型相對(duì)誤差分別為28.55%、91.55%、5.85%、31.94%。

    圖13 頂點(diǎn)轉(zhuǎn)角-彎矩曲線

    當(dāng)頂點(diǎn)彎矩加載到36.8kN·m時(shí)構(gòu)件處于塑性階段,分別提取各模型16個(gè)特征單元的S33應(yīng)力,統(tǒng)計(jì)精細(xì)化模型與多尺度模型在特征單元處的應(yīng)力值并計(jì)算相對(duì)誤差,觀察發(fā)現(xiàn),平均相對(duì)誤差從大到小分別為能量法模型78.5%、位移法模型27.7%、DC法模型6.6%、結(jié)合法模型0.7%。能量法模型的3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到239.2%。位移法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到50%。DC法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到7.1%。結(jié)合法模型的9、10、11、12號(hào)特征單元相對(duì)誤差最大,達(dá)到2.3%。在彎矩工況下,能量法模型的相對(duì)誤差率遠(yuǎn)大于其他多尺度模型,印證了能量協(xié)調(diào)法在彎矩作用下局限性的分析,既當(dāng)構(gòu)件進(jìn)入塑性階段后,彎矩產(chǎn)生的正應(yīng)力呈非線性分布,基于能量協(xié)調(diào)的多點(diǎn)約束方程不再適用。

    2.3 剪力工況分析

    在分析模型的底部設(shè)置固定約束,在頂部設(shè)置加載點(diǎn)通過(guò)施加沿X軸的剪力荷載控制加載,剪力從0逐漸增長(zhǎng)到20kN。

    各模型在塑性階段混凝土單元的應(yīng)力云圖如圖14所示,各模型在破壞階段混凝土受拉損傷云圖如圖15所示。通過(guò)觀察加載過(guò)程中應(yīng)力云圖的變化可以發(fā)現(xiàn),構(gòu)件在端部剪力的作用下,底部混凝土首先進(jìn)入塑性,并在拉壓兩側(cè)產(chǎn)生應(yīng)力集中。隨著不斷加載,構(gòu)件底部發(fā)生破壞,塑性變形快速增大。對(duì)比各模型可以發(fā)現(xiàn),位移法模型與結(jié)合法模型的應(yīng)力分布、損傷累積計(jì)算結(jié)果均較貼合精細(xì)化模型所得結(jié)果。能量法模型在多尺度界面出現(xiàn)了明顯應(yīng)力突變,損傷積累也與精細(xì)化模型有較大差別。

    圖14 剪力工況下塑性階段各模型混凝土軸向應(yīng)力云圖

    提取各模型在加載過(guò)程中加載點(diǎn)沿X軸的剪力和對(duì)應(yīng)沿X軸的位移,作出頂點(diǎn)位移-剪力曲線如圖16所示。5組模型加載結(jié)束時(shí)頂點(diǎn)位移分別為1.92、1.89、1.89、1.87、1.89mm。相比精細(xì)化模型,各多尺度模型相對(duì)誤差分別為1.22%、1.22%、2.25%。

    圖16 頂點(diǎn)位移-剪力曲線

    當(dāng)頂點(diǎn)剪力加載到20kN時(shí)構(gòu)件處于塑性階段,分別提取各模型16個(gè)特征單元的S33應(yīng)力,統(tǒng)計(jì)精細(xì)化模型與多尺度模型在特征單元處的應(yīng)力值并計(jì)算相對(duì)誤差見(jiàn)表1,觀察發(fā)現(xiàn),在剪力荷載工況下,平均相對(duì)誤差從大到小分別為能量法模型142.2%、位移法模型23.7%、DC法模型15.2%、結(jié)合法模型3.6%。能量法模型的1號(hào)特征單元相對(duì)誤差最大,達(dá)到341.6%。位移法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到33.3%。DC法模型的1、2、3、4號(hào)特征單元相對(duì)誤差最大,達(dá)到45.1%。結(jié)合法模型的6、7、8號(hào)特征單元相對(duì)誤差最大,達(dá)到5.4%。在剪力工況下,能量法模型的相對(duì)誤差率遠(yuǎn)大于其他多尺度模型,是因?yàn)樵诙瞬考袅Φ淖饔孟聵?gòu)件底部產(chǎn)生彎矩,混凝土進(jìn)入塑性后彎矩產(chǎn)生的正應(yīng)力呈非線性分布,基于能量協(xié)調(diào)的多點(diǎn)約束方程不再適用。

    表1 鋼筋混凝土節(jié)點(diǎn)平均應(yīng)力相對(duì)誤差統(tǒng)計(jì)

    3 結(jié)語(yǔ)

    文中針對(duì)鋼筋混凝土節(jié)點(diǎn),基于位移協(xié)調(diào)、能量協(xié)調(diào)、位移協(xié)調(diào)與能量協(xié)調(diào)相結(jié)合以及基于連續(xù)分布耦合理論實(shí)現(xiàn)多尺度界面的連接并進(jìn)行算例分析,所得具體研究成果如下:

    (1) 基于位移協(xié)調(diào)實(shí)現(xiàn)多尺度模型界面連接時(shí),除了在僅受軸力這種工況外,可以較為合理地反應(yīng)構(gòu)件的應(yīng)力分布和損傷累積。

    (2) 基于能量協(xié)調(diào)實(shí)現(xiàn)多尺度模型界面連接時(shí),除了在軸力工況下得到的結(jié)果較為精確外,在剪力工況、彎矩工況下,由于材料進(jìn)入塑性,彎矩產(chǎn)生的正應(yīng)力分布公式不再成立,約束方程不再適用,塑性區(qū)域單元迭代運(yùn)算逐漸出現(xiàn)誤差累積。

    (3) 基于位移協(xié)調(diào)和能量協(xié)調(diào)相結(jié)合實(shí)現(xiàn)多尺度模型界面連接時(shí),在各工況下均可合理地反應(yīng)構(gòu)件應(yīng)力分布和損傷累積,得到的結(jié)果也最接近精細(xì)化模型,但需要輸入大量約束方程,建模前處理較為復(fù)雜。

    (4) 基于連續(xù)分布耦合理論實(shí)現(xiàn)多尺度模型界面連接時(shí),在各工況下均可合理地反應(yīng)構(gòu)件的應(yīng)力分布和損傷累積,相比位移協(xié)調(diào)理論不僅可以有效改善側(cè)向過(guò)約束,在軸力工況下得到較為精確的結(jié)果,還能避免輸入大量的約束方程等繁瑣的建模前處理。

    猜你喜歡
    軸力彎矩尺度
    軸力作用下軟土基坑力學(xué)場(chǎng)演化規(guī)律的原位試驗(yàn)研究*
    淺談支撐軸力伺服系統(tǒng)在深基坑中的應(yīng)用
    財(cái)產(chǎn)的五大尺度和五重應(yīng)對(duì)
    基坑鋼支撐軸力監(jiān)測(cè)的優(yōu)化及實(shí)踐
    水道港口(2020年2期)2020-06-08 04:41:12
    零彎矩設(shè)計(jì)理論在連續(xù)梁橋中的應(yīng)用研究
    鋼支撐軸力在安裝過(guò)程中的監(jiān)測(cè)分析
    鐵道勘察(2016年3期)2016-08-01 07:20:14
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    CFRP-PCPs復(fù)合筋連續(xù)梁開(kāi)裂截面彎矩計(jì)算方法研究
    鋼-混疊合連續(xù)梁負(fù)彎矩區(qū)計(jì)算分析
    板孔式有彎矩平衡梁應(yīng)用技術(shù)及研究
    亚洲成国产人片在线观看| 欧美乱码精品一区二区三区| 精品欧美一区二区三区在线| 国产97色在线日韩免费| 女人精品久久久久毛片| 操美女的视频在线观看| 一本—道久久a久久精品蜜桃钙片| 成人三级做爰电影| 欧美中文综合在线视频| 亚洲情色 制服丝袜| 亚洲精品久久午夜乱码| 又黄又粗又硬又大视频| 欧美精品亚洲一区二区| 性色av一级| av线在线观看网站| 久热这里只有精品99| 欧美黄色片欧美黄色片| 99国产综合亚洲精品| 亚洲欧美一区二区三区黑人| 最近最新免费中文字幕在线| 制服诱惑二区| 美女主播在线视频| 国产成人av教育| 成人国产av品久久久| 国产在线视频一区二区| 午夜福利在线免费观看网站| 天堂中文最新版在线下载| 1024视频免费在线观看| 久9热在线精品视频| 亚洲免费av在线视频| 夫妻午夜视频| 久久国产亚洲av麻豆专区| 老司机深夜福利视频在线观看 | 老司机午夜福利在线观看视频 | 手机成人av网站| 无遮挡黄片免费观看| 黑丝袜美女国产一区| 女警被强在线播放| 日日夜夜操网爽| 久久久精品区二区三区| 99精品欧美一区二区三区四区| 久久久久国内视频| svipshipincom国产片| 一边摸一边抽搐一进一出视频| 无遮挡黄片免费观看| 国内毛片毛片毛片毛片毛片| 亚洲综合色网址| 欧美 亚洲 国产 日韩一| 精品亚洲乱码少妇综合久久| 51午夜福利影视在线观看| 桃花免费在线播放| www日本在线高清视频| 国产亚洲精品第一综合不卡| 国产伦人伦偷精品视频| 久久久国产一区二区| 欧美性长视频在线观看| 国产主播在线观看一区二区| 亚洲五月婷婷丁香| 久久国产亚洲av麻豆专区| 伊人久久大香线蕉亚洲五| 女人久久www免费人成看片| 久久香蕉激情| 美国免费a级毛片| 波多野结衣av一区二区av| 精品国产一区二区三区久久久樱花| 菩萨蛮人人尽说江南好唐韦庄| 久久国产亚洲av麻豆专区| 男人舔女人的私密视频| 女性生殖器流出的白浆| 精品欧美一区二区三区在线| 欧美成人午夜精品| 啦啦啦 在线观看视频| 麻豆国产av国片精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品国产av蜜桃| 十八禁网站免费在线| 欧美日韩视频精品一区| 精品卡一卡二卡四卡免费| 日韩欧美国产一区二区入口| 久久精品国产亚洲av高清一级| 久久久久国产一级毛片高清牌| 免费日韩欧美在线观看| 91麻豆精品激情在线观看国产 | 乱人伦中国视频| 国产一区二区三区在线臀色熟女 | 人妻 亚洲 视频| 亚洲国产欧美一区二区综合| 亚洲av电影在线观看一区二区三区| 国产精品国产三级国产专区5o| 亚洲成人国产一区在线观看| 国产欧美日韩综合在线一区二区| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲精品国产一区二区精华液| 国产老妇伦熟女老妇高清| 精品国产一区二区三区四区第35| a在线观看视频网站| 精品一区二区三区av网在线观看 | 波多野结衣av一区二区av| 亚洲欧美精品自产自拍| 欧美日韩亚洲国产一区二区在线观看 | 91成人精品电影| 免费观看人在逋| 91麻豆av在线| 热re99久久精品国产66热6| 91精品三级在线观看| 最黄视频免费看| 亚洲国产欧美网| av在线app专区| 日韩欧美免费精品| 欧美精品av麻豆av| 亚洲伊人色综图| 欧美激情高清一区二区三区| 中文字幕人妻熟女乱码| 夜夜夜夜夜久久久久| 宅男免费午夜| 国产精品99久久99久久久不卡| 精品亚洲成a人片在线观看| 韩国高清视频一区二区三区| 日本91视频免费播放| 国产不卡av网站在线观看| 在线观看免费日韩欧美大片| 王馨瑶露胸无遮挡在线观看| 两个人免费观看高清视频| 嫁个100分男人电影在线观看| 久久综合国产亚洲精品| 男女午夜视频在线观看| 男女免费视频国产| 精品国产乱子伦一区二区三区 | 十八禁人妻一区二区| 久久久久久久精品精品| 国产精品久久久av美女十八| 午夜福利在线观看吧| 两个人看的免费小视频| 亚洲全国av大片| 99国产综合亚洲精品| 亚洲avbb在线观看| 咕卡用的链子| 岛国在线观看网站| 女人精品久久久久毛片| 黄色a级毛片大全视频| 欧美精品亚洲一区二区| 亚洲国产毛片av蜜桃av| 日本精品一区二区三区蜜桃| 老司机午夜十八禁免费视频| 久久久久网色| 一二三四社区在线视频社区8| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看 | 精品亚洲乱码少妇综合久久| 大型av网站在线播放| 国产精品一区二区在线不卡| 久久人人97超碰香蕉20202| 精品国产乱码久久久久久男人| 国产男女内射视频| 999精品在线视频| 一个人免费在线观看的高清视频 | 十八禁高潮呻吟视频| 亚洲精品中文字幕一二三四区 | 天堂俺去俺来也www色官网| 国产亚洲欧美在线一区二区| 一级毛片电影观看| 亚洲国产欧美日韩在线播放| 天天操日日干夜夜撸| 久久热在线av| 亚洲欧美色中文字幕在线| 中亚洲国语对白在线视频| 精品一区二区三区av网在线观看 | 久久精品国产亚洲av香蕉五月 | 亚洲精品日韩在线中文字幕| 国产成人精品久久二区二区免费| 考比视频在线观看| 女警被强在线播放| 一边摸一边做爽爽视频免费| 国产免费现黄频在线看| av一本久久久久| 亚洲专区中文字幕在线| 久久免费观看电影| 黑人操中国人逼视频| 久久亚洲精品不卡| 两人在一起打扑克的视频| 一级a爱视频在线免费观看| 亚洲久久久国产精品| 亚洲av男天堂| 另类精品久久| 纵有疾风起免费观看全集完整版| 亚洲中文日韩欧美视频| 精品视频人人做人人爽| 一个人免费在线观看的高清视频 | 超碰97精品在线观看| 日韩电影二区| 黄频高清免费视频| 后天国语完整版免费观看| 欧美黄色淫秽网站| 一区二区三区激情视频| 最近最新中文字幕大全免费视频| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区激情| 欧美日韩视频精品一区| 丝袜脚勾引网站| 精品人妻熟女毛片av久久网站| 可以免费在线观看a视频的电影网站| 97精品久久久久久久久久精品| 亚洲精品粉嫩美女一区| 精品一区在线观看国产| 少妇精品久久久久久久| 99香蕉大伊视频| 欧美激情久久久久久爽电影 | 国产av又大| 日韩免费高清中文字幕av| 欧美老熟妇乱子伦牲交| 777米奇影视久久| 999精品在线视频| 狠狠狠狠99中文字幕| 国产精品一区二区免费欧美 | 男女国产视频网站| 少妇精品久久久久久久| 99久久国产精品久久久| 亚洲中文字幕日韩| 一级片免费观看大全| 中文字幕高清在线视频| 1024香蕉在线观看| 久久久国产一区二区| 男人舔女人的私密视频| 欧美成狂野欧美在线观看| 久久精品亚洲熟妇少妇任你| 午夜福利在线免费观看网站| 啦啦啦中文免费视频观看日本| 日韩电影二区| 国产成人欧美在线观看 | 一级黄色大片毛片| 精品免费久久久久久久清纯 | 多毛熟女@视频| 好男人电影高清在线观看| 伊人亚洲综合成人网| 精品高清国产在线一区| 男女床上黄色一级片免费看| 丰满人妻熟妇乱又伦精品不卡| 精品亚洲乱码少妇综合久久| 精品国产乱码久久久久久小说| 777久久人妻少妇嫩草av网站| 国产成人精品在线电影| 淫妇啪啪啪对白视频 | 精品一区二区三区四区五区乱码| 999久久久精品免费观看国产| 51午夜福利影视在线观看| 日日爽夜夜爽网站| 久久av网站| 亚洲国产欧美在线一区| 亚洲国产欧美在线一区| 午夜日韩欧美国产| 国产野战对白在线观看| 久久人妻熟女aⅴ| 色综合欧美亚洲国产小说| 五月天丁香电影| 日韩欧美一区二区三区在线观看 | 一边摸一边抽搐一进一出视频| 久久国产精品男人的天堂亚洲| 免费观看av网站的网址| 欧美日韩亚洲国产一区二区在线观看 | 久久人人97超碰香蕉20202| 欧美xxⅹ黑人| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 亚洲国产精品一区二区三区在线| 国产高清videossex| 丝袜在线中文字幕| 俄罗斯特黄特色一大片| 午夜成年电影在线免费观看| 999久久久精品免费观看国产| 国产免费现黄频在线看| 一级片'在线观看视频| 亚洲精品在线美女| 亚洲av国产av综合av卡| 曰老女人黄片| 在线看a的网站| 99国产精品99久久久久| 18禁国产床啪视频网站| 热99re8久久精品国产| 啦啦啦中文免费视频观看日本| 欧美日韩视频精品一区| 美女扒开内裤让男人捅视频| 国产亚洲av高清不卡| 90打野战视频偷拍视频| 精品久久久久久电影网| netflix在线观看网站| 国产日韩欧美亚洲二区| 纵有疾风起免费观看全集完整版| 久久久久精品国产欧美久久久 | 精品一区二区三卡| 久久精品成人免费网站| 久久影院123| 欧美日韩黄片免| 亚洲av电影在线进入| 色综合欧美亚洲国产小说| 国产亚洲一区二区精品| 国产成人精品在线电影| 91麻豆av在线| 啦啦啦啦在线视频资源| 大陆偷拍与自拍| 国产精品免费大片| 欧美亚洲日本最大视频资源| 在线精品无人区一区二区三| 日韩制服骚丝袜av| 欧美日韩精品网址| 亚洲第一青青草原| 女警被强在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 日韩一卡2卡3卡4卡2021年| 美女国产高潮福利片在线看| 亚洲久久久国产精品| 日韩中文字幕视频在线看片| 欧美日韩精品网址| 黄色视频不卡| 一级毛片电影观看| 日本五十路高清| 国产成人精品久久二区二区91| 久久ye,这里只有精品| 秋霞在线观看毛片| 成在线人永久免费视频| √禁漫天堂资源中文www| netflix在线观看网站| 男女无遮挡免费网站观看| 最黄视频免费看| 老司机在亚洲福利影院| h视频一区二区三区| 黄色片一级片一级黄色片| 狂野欧美激情性bbbbbb| 在线天堂中文资源库| 各种免费的搞黄视频| 亚洲精品一卡2卡三卡4卡5卡 | 欧美日韩亚洲综合一区二区三区_| 亚洲综合色网址| 人妻人人澡人人爽人人| 久久狼人影院| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影 | 性色av一级| 国产精品自产拍在线观看55亚洲 | 久久久久精品人妻al黑| 亚洲精品乱久久久久久| 99国产精品一区二区三区| 99国产精品免费福利视频| 国产色视频综合| 亚洲国产精品成人久久小说| 成年人免费黄色播放视频| 极品人妻少妇av视频| 一级黄色大片毛片| 十八禁网站网址无遮挡| 桃花免费在线播放| 国产精品1区2区在线观看. | 成人影院久久| 亚洲 国产 在线| 丝袜美腿诱惑在线| 高清av免费在线| 国产老妇伦熟女老妇高清| 国产淫语在线视频| 麻豆av在线久日| 一级毛片女人18水好多| 日韩制服丝袜自拍偷拍| 亚洲精品中文字幕一二三四区 | 美女扒开内裤让男人捅视频| av免费在线观看网站| 捣出白浆h1v1| 欧美少妇被猛烈插入视频| 丝袜在线中文字幕| 男女免费视频国产| 成人影院久久| 在线亚洲精品国产二区图片欧美| 国产福利在线免费观看视频| 久久人人爽av亚洲精品天堂| 成年av动漫网址| 国产成人影院久久av| 老司机亚洲免费影院| 十八禁网站网址无遮挡| 一本色道久久久久久精品综合| 久久人人爽人人片av| 亚洲人成77777在线视频| 在线观看免费视频网站a站| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| 亚洲av欧美aⅴ国产| 欧美精品av麻豆av| 老熟女久久久| 麻豆av在线久日| 日韩免费高清中文字幕av| 日韩欧美一区视频在线观看| 在线观看免费午夜福利视频| 亚洲午夜精品一区,二区,三区| h视频一区二区三区| 日韩有码中文字幕| 亚洲午夜精品一区,二区,三区| 亚洲精品美女久久久久99蜜臀| 国产成人一区二区三区免费视频网站| 国产精品熟女久久久久浪| 99re6热这里在线精品视频| 日本猛色少妇xxxxx猛交久久| 国产精品影院久久| 欧美激情高清一区二区三区| 脱女人内裤的视频| 人人妻人人爽人人添夜夜欢视频| 自拍欧美九色日韩亚洲蝌蚪91| 少妇粗大呻吟视频| 啪啪无遮挡十八禁网站| 亚洲国产精品一区三区| www.自偷自拍.com| 这个男人来自地球电影免费观看| 一本久久精品| 成年动漫av网址| 亚洲美女黄色视频免费看| 亚洲专区国产一区二区| 丝袜美腿诱惑在线| 两个人看的免费小视频| 在线永久观看黄色视频| 国产精品欧美亚洲77777| www.自偷自拍.com| 又大又爽又粗| 午夜激情av网站| 中文字幕另类日韩欧美亚洲嫩草| 久久人人爽av亚洲精品天堂| 日韩电影二区| 99热网站在线观看| 丝袜在线中文字幕| 亚洲人成电影观看| 色婷婷久久久亚洲欧美| 美女高潮到喷水免费观看| 日韩中文字幕视频在线看片| 欧美亚洲日本最大视频资源| 欧美激情高清一区二区三区| 黄色怎么调成土黄色| 一本一本久久a久久精品综合妖精| 免费久久久久久久精品成人欧美视频| 999久久久国产精品视频| 成人黄色视频免费在线看| 国产日韩欧美视频二区| 亚洲av欧美aⅴ国产| 99九九在线精品视频| 欧美日韩中文字幕国产精品一区二区三区 | 丁香六月欧美| 亚洲视频免费观看视频| 精品久久久久久电影网| 在线亚洲精品国产二区图片欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品一区二区在线不卡| 欧美精品高潮呻吟av久久| 日本一区二区免费在线视频| 黄网站色视频无遮挡免费观看| av超薄肉色丝袜交足视频| 一边摸一边抽搐一进一出视频| 亚洲九九香蕉| 丰满迷人的少妇在线观看| 欧美黑人精品巨大| av片东京热男人的天堂| 国产精品一二三区在线看| 国产成人精品久久二区二区91| 亚洲自偷自拍图片 自拍| 国产伦人伦偷精品视频| 日日爽夜夜爽网站| 午夜福利,免费看| 两个人看的免费小视频| 亚洲第一青青草原| 久久久精品94久久精品| 男人添女人高潮全过程视频| 国产区一区二久久| 国产深夜福利视频在线观看| 大香蕉久久成人网| 久久久久网色| 欧美日韩亚洲综合一区二区三区_| 免费一级毛片在线播放高清视频 | 一进一出抽搐动态| 欧美成人午夜精品| 国产精品 国内视频| 亚洲一区中文字幕在线| 久久中文字幕一级| 欧美日韩国产mv在线观看视频| 免费高清在线观看日韩| 黄色毛片三级朝国网站| 深夜精品福利| 9色porny在线观看| 欧美午夜高清在线| 午夜久久久在线观看| 久久毛片免费看一区二区三区| 国产亚洲午夜精品一区二区久久| 精品国产一区二区三区四区第35| av片东京热男人的天堂| 免费黄频网站在线观看国产| 久久久久久亚洲精品国产蜜桃av| 久久久久久久国产电影| 成人黄色视频免费在线看| 老熟妇仑乱视频hdxx| 久久精品亚洲av国产电影网| 日韩 亚洲 欧美在线| 亚洲七黄色美女视频| 9191精品国产免费久久| 90打野战视频偷拍视频| 99久久综合免费| 国产精品九九99| 黄色毛片三级朝国网站| 国产精品.久久久| 99久久99久久久精品蜜桃| 99re6热这里在线精品视频| 久久人人爽av亚洲精品天堂| 人人妻人人爽人人添夜夜欢视频| 午夜免费鲁丝| 亚洲av美国av| 国产精品一区二区免费欧美 | 亚洲av男天堂| 亚洲欧美色中文字幕在线| 老司机午夜福利在线观看视频 | 国产亚洲av高清不卡| 丝袜人妻中文字幕| 另类精品久久| 久9热在线精品视频| 久久久久久免费高清国产稀缺| 久久天躁狠狠躁夜夜2o2o| 岛国在线观看网站| 亚洲性夜色夜夜综合| 午夜福利视频精品| 免费高清在线观看日韩| 91精品伊人久久大香线蕉| 成年人免费黄色播放视频| 91老司机精品| 亚洲久久久国产精品| 天堂8中文在线网| 精品人妻熟女毛片av久久网站| 黄色a级毛片大全视频| 在线天堂中文资源库| 黑人操中国人逼视频| 香蕉丝袜av| 老司机深夜福利视频在线观看 | 天天躁日日躁夜夜躁夜夜| 久久毛片免费看一区二区三区| 久久久久精品国产欧美久久久 | 亚洲精品自拍成人| 欧美激情久久久久久爽电影 | 国产亚洲精品第一综合不卡| 操出白浆在线播放| av天堂在线播放| 高清av免费在线| 午夜福利乱码中文字幕| 国产精品 欧美亚洲| 少妇被粗大的猛进出69影院| 伦理电影免费视频| 黄色视频不卡| av电影中文网址| 性色av乱码一区二区三区2| 欧美激情久久久久久爽电影 | 国产91精品成人一区二区三区 | 亚洲欧美清纯卡通| 国产精品av久久久久免费| 高清欧美精品videossex| 国内毛片毛片毛片毛片毛片| 王馨瑶露胸无遮挡在线观看| 日韩视频一区二区在线观看| 亚洲欧美一区二区三区久久| 免费观看av网站的网址| 欧美精品亚洲一区二区| 国产精品久久久久久人妻精品电影 | 视频区图区小说| 日韩欧美国产一区二区入口| 最新在线观看一区二区三区| 美女午夜性视频免费| 伦理电影免费视频| 欧美黄色片欧美黄色片| 亚洲 欧美一区二区三区| 天天添夜夜摸| 老司机深夜福利视频在线观看 | 97人妻天天添夜夜摸| 男女边摸边吃奶| 亚洲欧美日韩另类电影网站| 后天国语完整版免费观看| 三上悠亚av全集在线观看| 可以免费在线观看a视频的电影网站| 午夜福利一区二区在线看| 亚洲 欧美一区二区三区| 黄频高清免费视频| 男人爽女人下面视频在线观看| 1024视频免费在线观看| 动漫黄色视频在线观看| 欧美日韩一级在线毛片| 不卡一级毛片| 69精品国产乱码久久久| 999久久久国产精品视频| 99久久精品国产亚洲精品| 十八禁高潮呻吟视频| 狂野欧美激情性xxxx| 精品一品国产午夜福利视频| 欧美少妇被猛烈插入视频| 亚洲国产精品999| 国产99久久九九免费精品| 亚洲av电影在线进入| 我要看黄色一级片免费的| 国产精品久久久久久人妻精品电影 | 亚洲欧洲精品一区二区精品久久久| 亚洲成人国产一区在线观看| 精品少妇久久久久久888优播| 国产在线免费精品| 女人爽到高潮嗷嗷叫在线视频| 亚洲av电影在线进入| 女性被躁到高潮视频| 美女大奶头黄色视频| 国产亚洲一区二区精品| 成人国产av品久久久| 69av精品久久久久久 | 美女脱内裤让男人舔精品视频| 热re99久久精品国产66热6| 人妻久久中文字幕网| 久久精品人人爽人人爽视色| 久久久精品区二区三区| 亚洲精品国产一区二区精华液| 天天添夜夜摸| 国产成人欧美| 亚洲精品国产色婷婷电影| 精品免费久久久久久久清纯 | 欧美亚洲日本最大视频资源| 午夜成年电影在线免费观看|