• <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ù)及研究
    国产一卡二卡三卡精品| 一边摸一边抽搐一进一小说| 中文字幕色久视频| netflix在线观看网站| 成人特级黄色片久久久久久久| 国产精品香港三级国产av潘金莲| 最近最新中文字幕大全免费视频| 成人亚洲精品av一区二区| 黑人操中国人逼视频| 12—13女人毛片做爰片一| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 亚洲精品国产色婷婷电影| 亚洲一区中文字幕在线| 国产主播在线观看一区二区| 免费人成视频x8x8入口观看| 两人在一起打扑克的视频| 国产成人精品久久二区二区免费| 97碰自拍视频| 俄罗斯特黄特色一大片| 少妇的丰满在线观看| 欧美日本视频| 亚洲国产毛片av蜜桃av| 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看| 一个人免费在线观看的高清视频| 无人区码免费观看不卡| netflix在线观看网站| 韩国精品一区二区三区| 国产亚洲精品久久久久5区| 欧美色欧美亚洲另类二区 | 国产在线观看jvid| 国产精品自产拍在线观看55亚洲| 国产区一区二久久| 99国产精品一区二区三区| 国产野战对白在线观看| 日日干狠狠操夜夜爽| 性欧美人与动物交配| 在线国产一区二区在线| 亚洲自拍偷在线| 免费高清在线观看日韩| 狂野欧美激情性xxxx| 一本综合久久免费| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 女人被狂操c到高潮| 国产亚洲精品一区二区www| 久久欧美精品欧美久久欧美| 亚洲九九香蕉| 精品午夜福利视频在线观看一区| 岛国在线观看网站| 精品国产一区二区三区四区第35| 人成视频在线观看免费观看| 欧美大码av| 亚洲三区欧美一区| 亚洲第一青青草原| 给我免费播放毛片高清在线观看| 99精品在免费线老司机午夜| 亚洲伊人色综图| 淫秽高清视频在线观看| 日韩一卡2卡3卡4卡2021年| 免费观看精品视频网站| av天堂在线播放| 嫩草影院精品99| 久久这里只有精品19| 丝袜美足系列| 亚洲伊人色综图| 99国产综合亚洲精品| 桃色一区二区三区在线观看| 超碰成人久久| 久久精品国产清高在天天线| 一区在线观看完整版| 正在播放国产对白刺激| 黄色a级毛片大全视频| 成人永久免费在线观看视频| 一级a爱视频在线免费观看| 久久国产乱子伦精品免费另类| 在线天堂中文资源库| 亚洲狠狠婷婷综合久久图片| 叶爱在线成人免费视频播放| 免费人成视频x8x8入口观看| 操出白浆在线播放| 午夜成年电影在线免费观看| 看片在线看免费视频| 9热在线视频观看99| 老熟妇仑乱视频hdxx| av欧美777| 亚洲av电影在线进入| 亚洲五月天丁香| 国产麻豆69| 国产aⅴ精品一区二区三区波| 日韩国内少妇激情av| 亚洲成人久久性| 亚洲第一青青草原| 午夜老司机福利片| 欧美日韩瑟瑟在线播放| 精品乱码久久久久久99久播| 久久久精品国产亚洲av高清涩受| 亚洲自偷自拍图片 自拍| 国产又爽黄色视频| 亚洲人成电影观看| 亚洲中文日韩欧美视频| 国产xxxxx性猛交| 69av精品久久久久久| 法律面前人人平等表现在哪些方面| 亚洲欧美精品综合一区二区三区| 国产欧美日韩精品亚洲av| 人人妻人人澡人人看| 丝袜美足系列| 久久婷婷人人爽人人干人人爱 | 欧美一区二区精品小视频在线| 亚洲成a人片在线一区二区| 欧美日本视频| 精品国产美女av久久久久小说| 久久性视频一级片| 精品国产国语对白av| 无遮挡黄片免费观看| cao死你这个sao货| 黄色丝袜av网址大全| 在线观看www视频免费| 亚洲精品美女久久av网站| 一进一出抽搐动态| 久久精品亚洲熟妇少妇任你| 久久精品aⅴ一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | netflix在线观看网站| 丁香欧美五月| 国产成人影院久久av| 一级毛片女人18水好多| 亚洲国产毛片av蜜桃av| 亚洲成av人片免费观看| 国产免费男女视频| 国产精品,欧美在线| 99在线视频只有这里精品首页| 精品久久久久久久久久免费视频| 九色亚洲精品在线播放| 亚洲第一av免费看| 久久午夜亚洲精品久久| 麻豆一二三区av精品| 午夜福利18| 亚洲aⅴ乱码一区二区在线播放 | 国产日韩一区二区三区精品不卡| 日本免费一区二区三区高清不卡 | 黄色成人免费大全| 久久久久久大精品| 国产片内射在线| 女人被躁到高潮嗷嗷叫费观| 亚洲性夜色夜夜综合| 老司机福利观看| 女同久久另类99精品国产91| 亚洲熟妇熟女久久| 精品久久蜜臀av无| 老司机靠b影院| 天天添夜夜摸| 91大片在线观看| 18禁观看日本| 久久久久久久久久久久大奶| 黄片播放在线免费| 日本五十路高清| 亚洲aⅴ乱码一区二区在线播放 | 50天的宝宝边吃奶边哭怎么回事| 免费少妇av软件| 高清毛片免费观看视频网站| 熟妇人妻久久中文字幕3abv| 国产精品二区激情视频| 中出人妻视频一区二区| 亚洲国产欧美日韩在线播放| 久久久久久久精品吃奶| 色综合婷婷激情| avwww免费| 久久午夜亚洲精品久久| 欧美日韩瑟瑟在线播放| 成人亚洲精品一区在线观看| 久久午夜综合久久蜜桃| 搡老岳熟女国产| 一本久久中文字幕| 色综合婷婷激情| 日本 av在线| 一级,二级,三级黄色视频| 成人亚洲精品av一区二区| 国产成人免费无遮挡视频| 午夜视频精品福利| 日韩视频一区二区在线观看| 极品人妻少妇av视频| 精品久久久久久久人妻蜜臀av | 69av精品久久久久久| 最新在线观看一区二区三区| 一级毛片精品| 亚洲伊人色综图| 亚洲一码二码三码区别大吗| 欧美国产日韩亚洲一区| 9热在线视频观看99| 久久亚洲精品不卡| 国产精品九九99| 欧美乱妇无乱码| xxx96com| 久久国产亚洲av麻豆专区| 国产一区二区三区视频了| 丰满的人妻完整版| 亚洲国产欧美一区二区综合| 精品第一国产精品| 成人国产一区最新在线观看| 琪琪午夜伦伦电影理论片6080| 又大又爽又粗| xxx96com| 国产成人欧美| 免费在线观看完整版高清| 色精品久久人妻99蜜桃| 校园春色视频在线观看| 美国免费a级毛片| 久久国产乱子伦精品免费另类| 欧美日本亚洲视频在线播放| 国产精品日韩av在线免费观看 | 一级a爱视频在线免费观看| 少妇被粗大的猛进出69影院| 欧美成人一区二区免费高清观看 | 超碰成人久久| 国产xxxxx性猛交| 国产成人精品无人区| 色综合欧美亚洲国产小说| 亚洲午夜精品一区,二区,三区| 亚洲国产看品久久| 高清黄色对白视频在线免费看| 欧美一级毛片孕妇| 国产精品99久久99久久久不卡| 亚洲视频免费观看视频| 久久天躁狠狠躁夜夜2o2o| 久久国产亚洲av麻豆专区| 日韩欧美一区视频在线观看| 美女 人体艺术 gogo| 亚洲欧美日韩另类电影网站| 精品免费久久久久久久清纯| 真人做人爱边吃奶动态| 亚洲av美国av| 91麻豆av在线| 91成人精品电影| 男人的好看免费观看在线视频 | 免费一级毛片在线播放高清视频 | 国产精品野战在线观看| 亚洲男人天堂网一区| 国产成人欧美在线观看| 久久久久亚洲av毛片大全| www.www免费av| 久久精品亚洲精品国产色婷小说| 国产成年人精品一区二区| 精品日产1卡2卡| 日本欧美视频一区| 美女 人体艺术 gogo| 国产精品亚洲一级av第二区| 黄频高清免费视频| 热re99久久国产66热| 国产欧美日韩一区二区三| 亚洲免费av在线视频| 日本黄色视频三级网站网址| 亚洲国产精品成人综合色| 精品卡一卡二卡四卡免费| 欧美黄色片欧美黄色片| 国产高清videossex| 亚洲精品在线观看二区| 女警被强在线播放| 国产精品亚洲一级av第二区| 国产精品 国内视频| 在线国产一区二区在线| 十八禁网站免费在线| 熟女少妇亚洲综合色aaa.| 亚洲伊人色综图| 久久午夜亚洲精品久久| 夜夜看夜夜爽夜夜摸| 欧美成人性av电影在线观看| 美国免费a级毛片| 免费高清视频大片| 91大片在线观看| 美女午夜性视频免费| 一进一出抽搐动态| 一区二区三区高清视频在线| 成年人黄色毛片网站| 免费无遮挡裸体视频| 国产精品爽爽va在线观看网站 | 精品卡一卡二卡四卡免费| 深夜精品福利| 日本在线视频免费播放| 久久国产精品男人的天堂亚洲| 美女 人体艺术 gogo| 欧美日韩黄片免| 最近最新免费中文字幕在线| 日韩国内少妇激情av| 精品久久久久久久毛片微露脸| 久久人妻av系列| 一a级毛片在线观看| 国产亚洲欧美98| 少妇粗大呻吟视频| 久9热在线精品视频| 久久国产精品影院| 18美女黄网站色大片免费观看| 日韩国内少妇激情av| 亚洲人成77777在线视频| 亚洲在线自拍视频| 免费在线观看完整版高清| 欧洲精品卡2卡3卡4卡5卡区| 美女免费视频网站| 精品久久蜜臀av无| 国产成人系列免费观看| 中文字幕色久视频| 国产男靠女视频免费网站| 99在线视频只有这里精品首页| 侵犯人妻中文字幕一二三四区| 性色av乱码一区二区三区2| 一级毛片高清免费大全| 视频区欧美日本亚洲| 午夜免费鲁丝| 国产成年人精品一区二区| 黄色丝袜av网址大全| 精品国产一区二区三区四区第35| 精品日产1卡2卡| 午夜福利18| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲国产欧美日韩在线播放| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲一级av第二区| 一本久久中文字幕| 亚洲国产高清在线一区二区三 | 在线观看免费午夜福利视频| 精品电影一区二区在线| 国产成人av教育| 老司机靠b影院| 国产亚洲av嫩草精品影院| 黑人巨大精品欧美一区二区mp4| 免费搜索国产男女视频| 九色亚洲精品在线播放| 麻豆av在线久日| 精品福利观看| 大型av网站在线播放| 欧美 亚洲 国产 日韩一| 极品教师在线免费播放| netflix在线观看网站| 琪琪午夜伦伦电影理论片6080| 精品国内亚洲2022精品成人| 国产一区二区在线av高清观看| 国产成人av激情在线播放| 欧美老熟妇乱子伦牲交| 黄色丝袜av网址大全| 久久精品国产亚洲av香蕉五月| 国产亚洲欧美精品永久| 日本 av在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩成人在线观看一区二区三区| 精品久久蜜臀av无| 一区福利在线观看| 亚洲国产精品久久男人天堂| 他把我摸到了高潮在线观看| 久久热在线av| www.999成人在线观看| 久久影院123| 1024香蕉在线观看| 亚洲情色 制服丝袜| 亚洲,欧美精品.| 欧美大码av| 1024香蕉在线观看| 久久人人爽av亚洲精品天堂| 十八禁网站免费在线| 黄色a级毛片大全视频| 日本精品一区二区三区蜜桃| 可以在线观看毛片的网站| 在线天堂中文资源库| 亚洲 欧美一区二区三区| 亚洲熟妇熟女久久| 男男h啪啪无遮挡| 久久久国产精品麻豆| 看片在线看免费视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产片内射在线| 日韩欧美一区二区三区在线观看| 级片在线观看| 精品乱码久久久久久99久播| 免费在线观看日本一区| 成人免费观看视频高清| av电影中文网址| 叶爱在线成人免费视频播放| 欧美绝顶高潮抽搐喷水| 亚洲国产精品成人综合色| 国产乱人伦免费视频| 亚洲精华国产精华精| 亚洲欧美精品综合一区二区三区| 大香蕉久久成人网| 身体一侧抽搐| 精品国产一区二区久久| 日韩成人在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| 多毛熟女@视频| 久久亚洲精品不卡| xxx96com| 国产三级黄色录像| 国产精品亚洲美女久久久| 黄频高清免费视频| 国产精品1区2区在线观看.| 视频在线观看一区二区三区| 成人特级黄色片久久久久久久| 欧美成狂野欧美在线观看| 少妇粗大呻吟视频| 亚洲成a人片在线一区二区| 久久久国产欧美日韩av| 精品午夜福利视频在线观看一区| 在线观看免费视频网站a站| 色尼玛亚洲综合影院| 色综合亚洲欧美另类图片| 精品第一国产精品| 女性生殖器流出的白浆| 中文字幕精品免费在线观看视频| 久久久久久久久久久久大奶| 成人特级黄色片久久久久久久| 男人操女人黄网站| 久久久久亚洲av毛片大全| 最新在线观看一区二区三区| 亚洲片人在线观看| 91大片在线观看| 亚洲av成人av| 欧美日本视频| 免费人成视频x8x8入口观看| 日本黄色视频三级网站网址| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| 欧美一级毛片孕妇| 国产私拍福利视频在线观看| 亚洲精品av麻豆狂野| 成在线人永久免费视频| 曰老女人黄片| 黄片大片在线免费观看| 一进一出抽搐gif免费好疼| bbb黄色大片| 色在线成人网| 亚洲男人的天堂狠狠| 在线永久观看黄色视频| 99精品在免费线老司机午夜| 亚洲九九香蕉| 午夜激情av网站| 亚洲五月婷婷丁香| www.自偷自拍.com| 嫁个100分男人电影在线观看| 久久久久国产精品人妻aⅴ院| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 男女之事视频高清在线观看| 精品国产国语对白av| 一个人免费在线观看的高清视频| 又黄又爽又免费观看的视频| 国产精华一区二区三区| 神马国产精品三级电影在线观看 | 18禁黄网站禁片午夜丰满| 91成人精品电影| 法律面前人人平等表现在哪些方面| 99热只有精品国产| 国产成人精品久久二区二区免费| 日日干狠狠操夜夜爽| 亚洲欧美日韩另类电影网站| 国产av一区二区精品久久| 色综合欧美亚洲国产小说| 国产欧美日韩一区二区三| 中文字幕久久专区| 中文字幕另类日韩欧美亚洲嫩草| 又紧又爽又黄一区二区| 国产国语露脸激情在线看| 97人妻精品一区二区三区麻豆 | 成人精品一区二区免费| 日本a在线网址| 精品不卡国产一区二区三区| bbb黄色大片| 久久国产精品男人的天堂亚洲| 国产人伦9x9x在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 此物有八面人人有两片| 精品欧美一区二区三区在线| 宅男免费午夜| 亚洲午夜理论影院| 亚洲国产欧美一区二区综合| 一进一出抽搐gif免费好疼| 女人被狂操c到高潮| 丰满的人妻完整版| 中文字幕高清在线视频| 国产成人精品久久二区二区91| 在线观看免费视频日本深夜| 久久久精品国产亚洲av高清涩受| 日韩精品免费视频一区二区三区| 怎么达到女性高潮| 欧美日韩乱码在线| 亚洲 欧美 日韩 在线 免费| 亚洲成a人片在线一区二区| 婷婷六月久久综合丁香| 国产精品野战在线观看| 男女下面插进去视频免费观看| 亚洲精品在线美女| 亚洲一区高清亚洲精品| 男女下面进入的视频免费午夜 | 免费在线观看视频国产中文字幕亚洲| 国产成人影院久久av| 97碰自拍视频| 久久伊人香网站| 亚洲精品粉嫩美女一区| 又黄又爽又免费观看的视频| 波多野结衣一区麻豆| 久久人人精品亚洲av| 国内久久婷婷六月综合欲色啪| 黄色视频不卡| 亚洲午夜精品一区,二区,三区| 亚洲色图 男人天堂 中文字幕| 色在线成人网| 黑人欧美特级aaaaaa片| 国产精品精品国产色婷婷| 久久国产精品影院| 国内精品久久久久久久电影| 国产精品久久久久久精品电影 | 自线自在国产av| 成人亚洲精品av一区二区| 亚洲成人国产一区在线观看| 精品国产乱子伦一区二区三区| 国产一区二区三区综合在线观看| 看免费av毛片| 国产亚洲精品综合一区在线观看 | 两性夫妻黄色片| 久久婷婷成人综合色麻豆| 极品人妻少妇av视频| 国产精品一区二区精品视频观看| 久久精品91无色码中文字幕| 精品国产超薄肉色丝袜足j| 欧美人与性动交α欧美精品济南到| 久久人人97超碰香蕉20202| 免费不卡黄色视频| 丁香欧美五月| 欧美国产日韩亚洲一区| 亚洲男人天堂网一区| 窝窝影院91人妻| 俄罗斯特黄特色一大片| 欧美绝顶高潮抽搐喷水| av视频在线观看入口| 日韩国内少妇激情av| 日韩视频一区二区在线观看| 色播在线永久视频| 90打野战视频偷拍视频| 色播在线永久视频| 国产午夜精品久久久久久| 久久久久精品国产欧美久久久| 国产高清激情床上av| 校园春色视频在线观看| 日本欧美视频一区| 大码成人一级视频| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 亚洲精品一区av在线观看| 国产熟女午夜一区二区三区| 色哟哟哟哟哟哟| 久久久久国内视频| 国产亚洲av高清不卡| 老鸭窝网址在线观看| 亚洲欧美精品综合一区二区三区| 亚洲国产欧美一区二区综合| 欧美精品亚洲一区二区| 大型av网站在线播放| 亚洲欧美日韩另类电影网站| 青草久久国产| 亚洲熟女毛片儿| 成人手机av| 激情视频va一区二区三区| 嫁个100分男人电影在线观看| 国产一级毛片七仙女欲春2 | 一区在线观看完整版| 12—13女人毛片做爰片一| 久99久视频精品免费| 自线自在国产av| 在线av久久热| 欧美在线黄色| 欧美日韩亚洲国产一区二区在线观看| 欧美黑人精品巨大| 欧美国产日韩亚洲一区| 一区福利在线观看| 国产精品永久免费网站| 亚洲激情在线av| 中文字幕最新亚洲高清| 在线免费观看的www视频| 久久热在线av| 国产欧美日韩精品亚洲av| 99精品在免费线老司机午夜| 国产一区二区在线av高清观看| 99久久99久久久精品蜜桃| 国产精品精品国产色婷婷| 两性午夜刺激爽爽歪歪视频在线观看 | 女警被强在线播放| 日本黄色视频三级网站网址| 欧美人与性动交α欧美精品济南到| 欧美成人午夜精品| 欧美乱妇无乱码| 免费在线观看影片大全网站| 国产欧美日韩综合在线一区二区| 嫩草影视91久久| 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 久9热在线精品视频| 亚洲 国产 在线| 欧美中文日本在线观看视频| 国产高清激情床上av| 国产日韩一区二区三区精品不卡| 亚洲,欧美精品.| 免费观看精品视频网站| 亚洲精品粉嫩美女一区| 亚洲国产日韩欧美精品在线观看 | 国产精品一区二区精品视频观看| 黄片大片在线免费观看| 18禁观看日本| www国产在线视频色| 国产免费av片在线观看野外av| 精品国内亚洲2022精品成人| 99国产综合亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 欧美乱妇无乱码|