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

    框架填充墻等效單桿斜撐簡化模型

    2017-03-22 02:58:06鄧宏宇孫柏濤
    關(guān)鍵詞:梁柱對(duì)角砌體

    鄧宏宇,孫柏濤

    ?

    框架填充墻等效單桿斜撐簡化模型

    鄧宏宇1, 3,孫柏濤1, 2

    (1. 哈爾濱工程大學(xué)航天與建筑工程學(xué)院,黑龍江哈爾濱,150001;2. 中國地震局工程力學(xué)研究所地震工程與工程振動(dòng)重點(diǎn)實(shí)驗(yàn)室,黑龍江哈爾濱,150080;3. 東北農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,黑龍江哈爾濱,150030)

    在已有的填充墻等效對(duì)角斜撐模型的基礎(chǔ)上,結(jié)合近年來地震現(xiàn)場(chǎng)普遍出現(xiàn)的填充墻頂部與框架梁柱節(jié)點(diǎn)接觸處的磚砌體在大震作用下通常會(huì)被壓碎這一實(shí)際震害現(xiàn)象,改變等效斜撐桿與框架柱的接觸位置,提出一種新的簡化計(jì)算模型。通過理論推導(dǎo)給出等效斜撐墻片寬度的計(jì)算公式,利用ABAQUS有限元模擬的方法,確定等效斜撐墻片與框架柱合理接觸位置的經(jīng)驗(yàn)計(jì)算公式。最后通過將簡化有限元模型計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比,驗(yàn)證該簡化建模方法的正確性。研究結(jié)果表明:傳統(tǒng)的對(duì)角斜撐模型低估框架填充墻構(gòu)件的水平承載力,新模型與蘆山地震框架結(jié)構(gòu)震害結(jié)果更相符。本模型可為框架結(jié)構(gòu)抗震分析與設(shè)計(jì)驗(yàn)算提供參考意義。

    填充墻;框架結(jié)構(gòu);等效斜撐模型;抗剪承載力;剛度

    鋼筋混凝土框架填充墻結(jié)構(gòu)是我國應(yīng)用最廣泛的結(jié)構(gòu)形式之一,近年來發(fā)生的地震中,由于填充墻與框架梁柱之間相互作用造成的震害十分嚴(yán)重[1]。為研究框架填充墻結(jié)構(gòu)的破壞機(jī)理,人們常常需要用數(shù)值模擬的方法,建立有限元模型進(jìn)行分析和計(jì)算。因此,研究與框架填充墻構(gòu)件實(shí)際受力最為接近的有限元模型的建模方法意義十分重大。在地震工程領(lǐng)域,大致有2類建模方法:一是微觀模型,即采用實(shí)體單元精細(xì)化模型,這種模型能夠精確的模擬填充墻體的變形、開裂和與框架梁柱的局部應(yīng)力作用,但也存在單元個(gè)數(shù)多、計(jì)算量大和不容易收斂的問題,無法應(yīng)用于整體模型的大規(guī)模計(jì)算。二是宏觀等效模型,即基于一定的假設(shè),將填充墻簡化為等效斜撐桿來模擬計(jì)算。等效斜撐的概念由POLYAJCOV[2]基于試驗(yàn)提出,他認(rèn)為填充墻可以簡化為在填充墻對(duì)角線方向上只受壓不受拉的斜撐桿。由于等效對(duì)角斜撐模型概念簡單、便于應(yīng)用、宏觀上描述了填充墻與框架之間的相互作用,因此在國內(nèi)外有廣泛的應(yīng)用。但該模型也存在一些明顯的不足,該計(jì)算模型等效斜撐桿的位置與填充墻對(duì)角線方向重合,上部桿端作用在梁柱節(jié)點(diǎn)處,這在彈性階段或者構(gòu)件剛進(jìn)入彈塑性階段是合理的。但是近年來的震害實(shí)際和框架填充墻構(gòu)件的擬靜力推覆試驗(yàn)表明[3]:在大震作用下,填充墻頂部與梁柱節(jié)點(diǎn)接觸處的角部磚砌體通常會(huì)被壓碎從而退出工作,這一破壞形式在地震現(xiàn)場(chǎng)比較普遍[4]。針對(duì)大震中框架結(jié)構(gòu)這種較常出現(xiàn)的破壞模式,在極限承載力條件下,仍按等效斜撐桿端支撐在梁柱節(jié)點(diǎn)處來處理是不盡合理的,用該方法建立的模型來模擬結(jié)構(gòu)在彈塑性階段的地震反應(yīng)與地震現(xiàn)場(chǎng)看到的實(shí)際的震害不符,因?yàn)樘畛鋲w頂部退出工作時(shí)墻體實(shí)際支撐位置下移,使得框架節(jié)點(diǎn)與支撐位置的部分框架承受很大的剪力,經(jīng)常造成框架根部剪切破壞或節(jié)點(diǎn)破壞。對(duì)于等效斜撐模型,最重要的參數(shù)就是確定等效斜撐墻片的有效寬度,國內(nèi)外的研究學(xué)者給出了多種計(jì)算方法。其中,HOLMES[5]建議等效斜撐墻片的有效寬度取填充墻對(duì)角線長度的1/3;SMITH等[6]通過理論分析引入相對(duì)剛度參數(shù)的概念,認(rèn)為等效墻片寬度與框架和填充墻的剛度比有關(guān);SANEINEJAD等[7]考慮了結(jié)構(gòu)的非線性,修正了單個(gè)等效支撐的面積公式,將等效墻片的有效寬度用與框架梁柱接觸長度的方式來表示;TUCKER[8]通過一系列的框架填充墻構(gòu)件的試驗(yàn),對(duì)比了18個(gè)等效斜撐寬度的計(jì)算公式,得出了1個(gè)較適合鋼筋混凝土框架砌體填充墻構(gòu)件的等效斜撐寬度公式。由此可見,不同的學(xué)者給出的計(jì)算公式差別很大,等效斜撐墻片寬度的取值有待于研究。基于以上的研究現(xiàn)狀,本文作者從實(shí)際工程應(yīng)用的角度出發(fā),在傳統(tǒng)的對(duì)角斜撐模型的基礎(chǔ)上,調(diào)整斜撐桿端與框架柱的接觸位置,將斜撐桿上端自梁柱節(jié)點(diǎn)處下移某一距離,提出一種新的等效單桿斜撐模型。結(jié)合理論分析和有限元模擬的方法,給出等效斜撐墻片寬度和等效斜撐桿與框架柱合理的接觸位置的計(jì)算公式。該簡化模型適用于大震作用下,在非線性階段,填充墻角部砌體被壓碎這一普遍出現(xiàn)的破壞模式,可應(yīng)用于整體結(jié)構(gòu)大規(guī)模的分析計(jì)算。

    1 震害實(shí)例

    1.1 結(jié)構(gòu)破壞情況

    2013?04?20T08:02,中國四川省雅安市蘆山縣發(fā)生7.0級(jí)地震,震源深度13 km,震中烈度為Ⅸ度。地震造成大量的人員傷亡和房屋建筑破壞。本文作者在第一時(shí)間進(jìn)入雨城區(qū)上里鎮(zhèn)進(jìn)行震害考察,發(fā)現(xiàn)該地區(qū)鋼筋混凝土框架結(jié)構(gòu)的震害十分嚴(yán)重。尤其是教學(xué)樓結(jié)構(gòu),多采用單跨外廊式的鋼筋混凝土框架結(jié)構(gòu)。圖1所示為蘆山地震雨城區(qū)上里中學(xué)教學(xué)樓震害情況。在地震作用下,框架柱柱頭受剪破壞,填充墻嚴(yán)重開裂破壞,填充墻與框架梁柱節(jié)點(diǎn)處的角部磚砌體被壓碎。這種破壞形式在震區(qū)框架結(jié)構(gòu)建筑中十分普遍[9]。為研究框架填充墻構(gòu)件的受力機(jī)制和破壞機(jī)理,下面用ABAQUS有限元模擬的方法,模擬計(jì)算圖1所示的框架填充墻構(gòu)件。

    1.2 實(shí)體單元模型

    首先,通過ABAQUS有限元軟件建立框架填充墻構(gòu)件的實(shí)體單元有限元模型。有研究表明,實(shí)體單元模型的計(jì)算結(jié)果與實(shí)際構(gòu)件最為接近[10]。為驗(yàn)證實(shí)體單元計(jì)算結(jié)果的準(zhǔn)確性,根據(jù)MEHRABI等[11]所做的一系列框架填充墻構(gòu)件模型擬靜力試驗(yàn),模擬其中的構(gòu)件1、構(gòu)件8和構(gòu)件9,并將有限元計(jì)算結(jié)果與通過試驗(yàn)所得實(shí)際數(shù)據(jù)進(jìn)行對(duì)比。其中構(gòu)件1為空框架,構(gòu)件8為空心砌塊填充墻,構(gòu)件9為實(shí)心磚填充墻,相關(guān)的試件尺寸和材料參數(shù)見文獻(xiàn)[11]。

    ABAQUS中建模過程如下:框架梁柱采用實(shí)體單元C3D20R,鋼筋采用桁架單元T3D2,本構(gòu)關(guān)系選擇對(duì)應(yīng)的混凝土和鋼筋塑性本構(gòu)[12],填充墻采用實(shí)體單元C3D20R,本構(gòu)關(guān)系采用砌體塑性本構(gòu)[13],將材料的本構(gòu)關(guān)系賦予到ABSQUS中的混凝土損傷塑性模型進(jìn)行彈塑性的靜力模擬分析。鋼筋與混凝土框架梁柱之間采用embedded連接,填充墻與框架梁柱之間采用contact連接模擬填充墻與框架梁柱之間的相互作用。采用隱式靜力分析的方法,在梁端施加側(cè)向荷載,最終能夠得到1條荷載?位移曲線見圖2。

    (a) 框架填充墻破壞;(b) 填充墻角部被壓碎

    從圖2中試件荷載?位移曲線的試驗(yàn)值與各自對(duì)應(yīng)的實(shí)體單元有限元模擬值的對(duì)比可以看出:2條曲線的最大水平承載力值比較吻合,驗(yàn)證了采用實(shí)體單元有限元模型靜力分析得到的構(gòu)件的最大水平承載力模擬結(jié)果與實(shí)際相符。因此,在后續(xù)的分析過程中,認(rèn)為實(shí)體單元模型的計(jì)算結(jié)果即為構(gòu)件的實(shí)際值。

    1—實(shí)心磚填充墻;2—空心砌塊填充墻;3—空框架。

    1.3 震害實(shí)例計(jì)算結(jié)果

    將圖1中的框架填充墻構(gòu)件單獨(dú)提取出來,采用已驗(yàn)證的實(shí)體單元有限元建模方法對(duì)其進(jìn)行模擬計(jì)算分析,并將計(jì)算結(jié)果與實(shí)際結(jié)果進(jìn)行對(duì)比。圖3所示為構(gòu)件達(dá)到極限水平承載力后某一時(shí)刻的Mises應(yīng)力云圖。

    從圖3(a)可以看出:隨著框架填充墻構(gòu)件側(cè)向變形的增加,填充墻右上角的磚砌體受到框架梁柱節(jié)點(diǎn)的擠壓,當(dāng)超過砌體材料的極限變形能力時(shí),磚砌體就會(huì)因?yàn)閴核槎顺龉ぷ?。從填充墻的Mises應(yīng)力云圖可以看出填充墻中應(yīng)力的分布和變化趨勢(shì),在初始階段,填充墻右上角的角部磚砌體應(yīng)力最大,對(duì)角斜撐的作用比較明顯。當(dāng)位移進(jìn)一步增大時(shí),角部的最大應(yīng)力會(huì)有一定程度的下移(圖3(b)),此時(shí)對(duì)角斜撐模型不再適用。圖3(c)中,從框架柱與填充墻接觸面的壓應(yīng)力變化也可以看出,最大應(yīng)力的位置隨著框架側(cè)向位移的增大而逐漸下移。以上計(jì)算結(jié)果最終狀態(tài)與實(shí)際震害比較相符,說明高烈度區(qū)的結(jié)構(gòu)遭到大震作用時(shí),填充墻對(duì)框架的對(duì)角斜撐位置有一定的下移,原來始終保持在對(duì)角線方向上的對(duì)角斜撐模型不盡合理,應(yīng)對(duì)模型進(jìn)行改進(jìn)和調(diào)整,才能更好地適用非線性階段的模擬分析。

    (a) 整體Mises應(yīng)力云圖;(b) 填充墻Mises應(yīng)力云圖;(c) 接觸面處壓應(yīng)力變化

    基于以上的實(shí)際震害和有限元模擬的初步分析,提出一種適用于普遍出現(xiàn)的填充墻角部砌體被壓碎的受力模型,根據(jù)該模型,從受力分析的角度,推導(dǎo)出等效斜撐墻片寬度的計(jì)算公式,并通過設(shè)計(jì)一系列構(gòu)件對(duì)其計(jì)算結(jié)果進(jìn)行統(tǒng)計(jì)回歸,給出等效斜撐桿與框架柱接觸位置的經(jīng)驗(yàn)計(jì)算公式。

    2 等效斜撐墻片寬度

    2.1 計(jì)算簡圖和基本假設(shè)

    單榀框架填充墻構(gòu)件的計(jì)算簡圖如圖4所示,在地震荷載作用下,認(rèn)為只有一定寬度的斜撐墻片參與受力,其與框架柱的接觸長度為c,墻片厚度與填充墻厚度相同,等效斜撐墻片的中軸線與水平方向的夾角為′。針對(duì)該模型,做出如下基本假設(shè):1) 普通磚砌體填充墻為各向同性的理想勻質(zhì)材料。2) 在以下的受力分析中,為簡便起見假設(shè)填充墻的邊界與框架梁柱的中軸線一致,即′≈,′≈,其中′和′為框架軸線處的高度和跨度,和為填充墻的高度和跨度。3) 簡化的斜撐填充墻頂部在構(gòu)件進(jìn)入極限承載力,角部磚砌體被壓碎之后認(rèn)為與梁沒有接觸,此時(shí)斜撐墻片的頂端與框架梁之間存在一定的距離,且與柱的接觸長度不超過0.3倍填充墻高,見式(1)。4) 框架柱、框架梁內(nèi)最大彈性彎矩與屈服彎矩的比值不超過0.2[14]。見下式:

    c≤0.3(1)

    ≤(2)

    式中:c為等效斜撐墻片與框架柱的接觸長度與框架柱高度的比值;為填充墻的高度;c和b為框架柱、框架梁的最大彈性彎矩;pc和pb分別為框架柱、框架梁的屈服彎矩;c和b分別為框架柱、框架梁內(nèi)彈性彎矩與屈服彎矩的比值;0為框架柱、框架梁內(nèi)彈性彎矩與屈服彎矩的比值最大值。

    圖4 框架填充墻構(gòu)件在角部磚砌體被壓碎后的計(jì)算簡圖及受力

    2.2 理論公式推導(dǎo)

    等效接觸面長度指的是等效斜撐墻片與框架柱的接觸長度。這里等效斜撐墻片的寬度以等效接觸面長度的形式給出。在極限水平承載力作用下,等效斜撐墻片在與框架柱接觸面處受剪、受壓綜合作用而導(dǎo)致破壞。根據(jù)特雷斯卡(Tresca)六棱柱體屈服準(zhǔn)則[14]得:

    (4)

    式中:為填充墻等效斜撐墻片與框架柱之間均布正應(yīng)力;為填充墻等效斜撐墻片與框架柱之間均布剪應(yīng)力;為框架梁柱的高跨比,即=′/′;為等效斜撐填充墻與框架柱之間的摩擦因數(shù),一般取0.45[16]。聯(lián)立式(3)和式(4),得出接觸面均布正應(yīng)力為

    式中:c0為填充墻與框架柱之間名義均布正應(yīng)力,c為填充墻材料的有效抗壓強(qiáng)度,可按計(jì)算;為填充墻材料的抗壓強(qiáng)度;為調(diào)整系數(shù),可取[16]。

    圖5所示為框架填充墻受力變形圖。由圖5可見:在水平荷載作用下,假設(shè)等效斜撐填充墻片與框架柱分離的點(diǎn)為,為簡化分析,在接觸長度范圍內(nèi)均布剪應(yīng)力忽略不計(jì),將柱段單獨(dú)提取出來作為分離體進(jìn)行受力分析,根據(jù)彎矩平衡得

    式中:c為框架柱彈性彎矩;pj為節(jié)點(diǎn)塑性抵抗彎矩,取節(jié)點(diǎn)處框架梁、柱的屈服彎矩及節(jié)點(diǎn)屈服彎矩的最小值。將式(2)代入式(6),并結(jié)合式(1),求解出等效接觸面長度c為

    (7)

    通過以上的推導(dǎo),將等效斜撐墻片寬度以斜撐墻片與框架柱接觸長度c的形式給出。當(dāng)1個(gè)框架填充墻構(gòu)件的尺寸、材料等級(jí)確定之后,除之外的其余參數(shù)就都已經(jīng)確定。即等效斜撐墻片自梁柱節(jié)點(diǎn)處下移的距離,不僅是等效斜撐墻片寬度計(jì)算公式的必要參數(shù)之一,也是將實(shí)體單元模型簡化為單桿斜撐模型,從而可以應(yīng)用于整體結(jié)構(gòu)大規(guī)模的分析計(jì)算的關(guān)鍵。以下介紹經(jīng)驗(yàn)計(jì)算公式的確定過程和簡化模型的建模方法。

    (a) 整體框架變形示意圖;(b) 柱EA段受力分析圖

    3 斜撐桿與框架柱的接觸位置

    3.1 構(gòu)件設(shè)計(jì)

    設(shè)計(jì)一系列層高′和跨度′不同,跨高比相同的框架填充墻構(gòu)件。其中框架梁、柱的截面尺寸和配筋相同,框架柱截面尺寸(長×寬)為350 mm×350 mm,框架梁截面尺寸(長×寬)為500 mm×250 mm,混凝土強(qiáng)度等級(jí)取C30,保護(hù)層厚度取35 mm,縱向受力鋼筋采用直徑為16 mm的HRB335鋼筋,箍筋采用直徑為8 mm的HPB235鋼筋,填充墻砌體強(qiáng)度等級(jí)為MU10/M5,填充墻厚度取120 mm,框架的抗震等級(jí)為3級(jí),設(shè)計(jì)軸壓比取0.2,跨高比均為1.5。具體設(shè)計(jì)參數(shù)見表1。

    表1 構(gòu)件設(shè)計(jì)參數(shù)

    3.2的確定

    首先采用實(shí)體單元模型的建模方法,在ABAQUS中對(duì)表1中的各個(gè)構(gòu)件進(jìn)行模擬計(jì)算。然后在實(shí)體單元有限元模型的基礎(chǔ)上,采用桁架單元T3D2,將實(shí)體單元的填充墻替換為等效斜撐桿。桿兩端與框架之間采用join連接,即節(jié)點(diǎn)之間允許有相對(duì)轉(zhuǎn)角,但不允許有相對(duì)位移,對(duì)應(yīng)于框架梁柱和填充墻之間只傳遞剪力,不傳遞彎矩的情況;材料屬性與填充墻使用的砌體本構(gòu)關(guān)系一樣,在設(shè)置材料屬性時(shí)規(guī)定只承受壓力不能承受拉力。簡化模型情況如圖6所示。

    圖6 等效單桿斜撐模型

    逐漸調(diào)整的取值即變化等效斜撐墻片與框架柱的接觸位置,直至等效單桿斜撐模型的荷載?位移曲線與實(shí)體單元有限元模型的荷載?位移曲線相吻合,記錄此時(shí)的,認(rèn)為此時(shí)的值即為該構(gòu)件的合理取值。各個(gè)構(gòu)件框架梁和框架柱的屈服彎矩按下式計(jì)算:

    (9)

    從圖7可以看出:當(dāng)3=350 mm時(shí),簡化斜撐模型與實(shí)體單元模型計(jì)算的荷載?位移曲線在平臺(tái)段最為符合,因此認(rèn)為對(duì)于構(gòu)件FR-1,的合理取值為350 mm。按照上述方法,對(duì)其余9個(gè)構(gòu)件采用同樣的分析過程,找出當(dāng)?shù)刃螚U斜撐模型的荷載?位移曲線的水平承載力最大值與實(shí)體單元有限元模型的計(jì)算結(jié)果相符時(shí)的取值。

    3.3 經(jīng)驗(yàn)公式擬合

    通過式(7)的推導(dǎo)看出, 等效壓桿的有效寬度可以通過含的表達(dá)式來表示,兩者之間成正比例關(guān)系。有研究認(rèn)為,等效壓桿的有效寬度與′cos′和′有關(guān)[17]。再結(jié)合觀察有限元實(shí)體單元模型的分析結(jié)果,本文認(rèn)為與′cos′有關(guān),/′cos′可以寫成含有′的函數(shù)表達(dá)式。其中,為相對(duì)剛度參數(shù),是由SMITH等[6]提出的,其計(jì)算公式為

    式中:m為填充墻材料的彈性模量;c為混凝土的彈性模量;為填充墻厚度;為填充墻對(duì)角線與水平方向的夾角;c為框架柱的截面慣性矩;為填充墻的高度。

    本文認(rèn)為/(′cos′)和′之間存在函數(shù)關(guān)系,分別計(jì)算上述設(shè)計(jì)的10個(gè)構(gòu)件中當(dāng)取合適的值時(shí)對(duì)應(yīng)的/(′cos′)和′,將/(′cos′)作為縱坐標(biāo),′作為橫坐標(biāo),并放在同一直角坐標(biāo)系下,得到一系列的分布點(diǎn),如圖8所示。

    對(duì)以上這些分布點(diǎn)擬合曲線得到的經(jīng)驗(yàn)計(jì)算式為

    通過式(11)可以計(jì)算,即等效斜撐墻片自框架梁柱節(jié)點(diǎn)下移的距離。當(dāng)1個(gè)框架填充墻構(gòu)件的截面尺寸和材料等級(jí)確定后,可通過式(11)和式(7)計(jì)算出等效斜撐桿與框架柱的接觸位置和等效斜撐墻片的寬度,即可建立對(duì)應(yīng)的等效簡化有限元模型。

    1—實(shí)體單元模型;2—斜撐墻片下降1=150 mm;3—斜撐墻片下降2=250 mm;4—斜撐墻片下降3=350 mm;5—斜撐墻片下降4=450 mm;6—斜撐墻片下降5=550 mm。

    圖7 構(gòu)件FR-1荷載?位移曲線

    Fig. 7 Load?displacement curve of member FR-1

    圖8 分布點(diǎn)和公式擬合曲線

    4 實(shí)例應(yīng)用

    4.1 試驗(yàn)及相關(guān)數(shù)據(jù)

    Shine和Mehrabi為研究填充墻對(duì)框架性能的影響,進(jìn)行了一系列共12個(gè)框架填充墻構(gòu)件的擬靜力試驗(yàn)。本文對(duì)其中的構(gòu)件8和構(gòu)件9進(jìn)行模擬驗(yàn)證。構(gòu)件8的填充墻為混凝土空心砌塊,構(gòu)件9的填充墻為實(shí)心磚砌體,構(gòu)件尺寸、材料參數(shù)和實(shí)驗(yàn)情況詳見文獻(xiàn)[11]。

    4.2 有限元計(jì)算與實(shí)驗(yàn)結(jié)果對(duì)比

    對(duì)以上2個(gè)構(gòu)件,分別采用實(shí)體單元接觸連接模型、本文提出的改進(jìn)的等效簡化斜撐模型和傳統(tǒng)的對(duì)角斜撐模型進(jìn)行有限元模擬。其中,傳統(tǒng)的對(duì)角斜撐簡化模型的等效墻片寬度的取值采用TUCKER[8]提出的計(jì)算公式:

    式中:為等效斜撐的寬度;′為填充墻對(duì)角線的長度;′為框架的高度;為相對(duì)剛度參數(shù)。TUCKER[8]提出的等效斜撐寬度的計(jì)算公式在對(duì)角斜撐模型中應(yīng)用最為廣泛,具有較普遍的意義[18]。構(gòu)件8、構(gòu)件9的3種有限元建模方法計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果的比較見圖9。

    從圖9可以看出:采用本文所述的簡化建模方法所得到的框架填充墻構(gòu)件的最大水平承載力與試驗(yàn)所得數(shù)據(jù)和采用實(shí)體單元接觸連接模型的計(jì)算結(jié)果三者之間比較吻合,驗(yàn)證了該等效簡化模型的合理性和準(zhǔn)確性。對(duì)于傳統(tǒng)的對(duì)角斜撐模型,始終保持在對(duì)角線方向上的對(duì)角斜撐模型不盡合理,等效墻片寬度采用TUCKER[8]公式所得的最大水平承載力低于上述3種結(jié)果,說明該對(duì)角斜撐模型低估了框架填充墻構(gòu)件的最大水平承載力,與構(gòu)件的實(shí)際受力誤差較大。

    (a) 構(gòu)件8;(b) 構(gòu)件9

    5 結(jié)論

    1) 在等效對(duì)角斜撐桿模型的基礎(chǔ)上,結(jié)合震害實(shí)例中填充墻頂部與框架梁柱節(jié)點(diǎn)接觸處的磚砌體在大震作用下通常會(huì)被壓碎的情況,將等效斜撐桿與框架柱的接觸位置下移一定的距離,提出一種新的簡化計(jì)算模型。在該模型的基礎(chǔ)上進(jìn)行受力分析和理論公式推導(dǎo),給出了等效斜撐墻片寬度的計(jì)算公式。等效斜撐桿與框架柱接觸位置的經(jīng)驗(yàn)計(jì)算公式則是通過設(shè)計(jì)一系列框架填充墻構(gòu)件進(jìn)行有限元模擬,并對(duì)結(jié)果進(jìn)行統(tǒng)計(jì)回歸的方法得到。由此可以建立一種改進(jìn)的等效單桿斜撐簡化模型,與震害實(shí)際更相符。

    2) 實(shí)體單元模型可以很好的模擬框架填充墻構(gòu)件的受力情況,但計(jì)算單元多,計(jì)算時(shí)間長且不容易收斂。按本文所介紹的等效簡化模型的數(shù)值模擬結(jié)果與實(shí)體單元模型計(jì)算結(jié)果和實(shí)驗(yàn)所得實(shí)測(cè)數(shù)據(jù)三者之間吻合得很好,采用等效簡化模型不僅大大節(jié)省了計(jì)算成本,同時(shí)也使得帶填充墻的框架結(jié)構(gòu)整體模型的大規(guī)模時(shí)程分析計(jì)算成為可能。

    3) 該簡化計(jì)算模型存在一定的局限性,即只適用于在大震作用下,填充墻角部磚砌體被壓碎這一破壞模式的模擬計(jì)算。在彈性階段,填充墻與框架梁柱之間的傳力方向是按對(duì)角線方向傳遞的,此時(shí)采用已有的對(duì)角斜撐模型即可。只有在構(gòu)件進(jìn)入非線性后,框架填充墻構(gòu)件才會(huì)普遍出現(xiàn)填充墻與框架柱接觸處最大應(yīng)力區(qū)域下移的情況。對(duì)于處于高烈度區(qū)、嚴(yán)重破壞且符合本文所述的破壞模式的框架填充墻結(jié)構(gòu),采用本文所介紹的簡化有限元模型進(jìn)行分析更符合震害實(shí)際,模擬結(jié)果也更準(zhǔn)確。

    [1] 葉列平, 曲哲, 馬千里, 等. 從汶川地震框架結(jié)構(gòu)震害談“強(qiáng)柱弱梁”屈服機(jī)制的實(shí)現(xiàn)[J]. 建筑結(jié)構(gòu), 2008, 38(11): 52?67. YE Lieping, QU Zhe, MA Qianli, et al. Study on ensuring the strong column?weak beam mechanism for RC frames based on the damage analysis in the Wenchuan earthquake[J]. Building Structure, 2008, 38(11): 52?67.

    [2] POLYAJCOV S V. On the interactions between masonry filler walls and enclosing frame when loaded in the plane of the wall[J]. Translations in Earthquake Engineering Research Institute, 1956, 32(2): 1?25.

    [3] 金煥.填充墻RC框架結(jié)構(gòu)地震破壞機(jī)理及關(guān)鍵抗震措施研究[D]. 哈爾濱: 中國地震局工程力學(xué)研究所, 2014: 35?47. JIN Huan. Research on seismic failure mechanisms and key aseismic measures for masonry-infilled RC frame structures[D]. Harbin: Institute of Engineering Mechanics, 2014: 35?47.

    [4] SUN Baitao, YAN Peilei, CHEN Xiangzhao. New characteristics of intensity assessment of Sichuan Lushan“4.20”Ms7.0 earthquake[J]. Earthquake Engineering and Engineering Vibration, 2014, S1(13): 123?139.

    [5] HOLMES M. Combined loading on infilled frames[J]. Proceedings of the Institution of Civil Engineers, 1963, 25(5): 31?38.

    [6] SMITH B S, CARTER C. A method of analysis for infilled frames[J]. Proceedings Institution of Civil Engineers, 1969, 44(1): 31?48.

    [7] SANEINEJAD A, HOBBS B. Inelastic design of infilled frames[J]. Journal of Structural Engineering, 1995, 121(4): 634?650.

    [8] TUCKER C J. Predicting the in-plane capacity of masonry infilled frames[D]. USA: The University of Tennessee. College of Civil Engineering, 2007: 57?69.

    [9] 孫柏濤, 閻培雷, 王明振, 等. 四川省蘆山“4.20”7.0級(jí)強(qiáng)烈地震建筑物震害圖集[M]. 北京: 地震出版社, 2014: 11?13. SUN Baitao, YAN Peilei, WANG Mingzhen, et al. Earthquake damage atlas of Sichuan Lushan“4.20”Ms7.0 earthquake[M]. Beijing: Earthquake Publishing House, 2014: 11?13.

    [10] KOUTROMANOS I, STAVRIDIS A, SHING P B, et al. Numerical modeling of masonry-infilled RC frames subjected to seismic loads[J]. Computers and Structures, 2011, 89(1): 1026?1037.

    [11] MEHRABI A B, SHING P B, SCHULLER M P, et al. Experimental evaluation of masonry-infilled RC frames[J]. Journal of Structural Engineering, 1996, 122(3): 228?237.

    [12] GB 50010—2010, 混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[S]. GB 50010—2010, Code for design of concrete structures[S].

    [13] 劉桂秋. 砌體結(jié)構(gòu)基本受力性能的研究[D]. 長沙: 湖南大學(xué)土木工程學(xué)院, 2005: 21?25. LIU Guiqiu. The research on the basic mechanical behavior of masonry structure[D]. Changsha: Hunan University. College of Civil Engineering, 2005: 21?25.

    [14] ACI 318M—05, Building code requirements for structural concrete[S].

    [15] CHEN W F. Plasticity in reinforced concrete[M]. New York:McGraw-Hill Book Co, 1982: 29?38.

    [16] ACI530—02, Building code requirements for masonry structures[S].

    [17] 金曉飛, 高松召, 楊曉杰, 等. 內(nèi)嵌墻板鋼框架等效單斜撐簡化模型[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2013, 45(6): 16?21. JIN Xiaofei, GAO Songzhao, YANG Xiaojie, et al. Equivalent single strut simplified model for steel frames with embedded retaining wall[J]. Journal of Harbin institute of technology, 2013, 45(6): 16?21.

    [18] 閻紅霞. 樓板和填充墻對(duì)RC框架結(jié)構(gòu)抗震性能的影響[D]. 北京: 北京交通大學(xué), 2012: 121?133. YAN Hongxia. Influences of slabs and infill walls on seismic performance of reinforced concrete frame structures[D]. Beijing: Beijin Jiaotong University, 2012: 121?133.

    (編輯 陳愛華)

    Equivalent diagonal bracing model of a frame-filled wall

    DENG Hongyu1, SUN Baitao1, 2

    (1. College of Space Technology and Civil Engineering, Harbin Engineering University, Harbin 150001, China;2. Key Laboratory of Earthquake Engineering and Engineering Vibration,Institute of Engineering Mechanics, China Earthquake Administration, Harbin 150080, China;3. School of Water Conservancy and Civil Engineering, Northeast Agricultural University, Harbin 150030, China)

    Based on the existing equivalent diagonal strut model, combined with the actual earthquake damage phenomena that the brick masonry between the top of the infill wall and the beam-column connections is often crushed during a strong earthquake, and changing the contact position between the equivalent bracing wall and column, a new calculation model was generated. Through theoretical derivation, the formulae for calculating the contact length between the equivalent bracing wall and column was given. Using the ABAQUS finite element simulation, an empirical formula for the reasonable contact position between the equivalent bracing wall and column was determined. Finally through comparing the finite element model calculation results with the experimental result, the validity of the simplified modeling method was verified. The results show that the traditional diagonal strut model underestimates the horizontal bearing capacity of the infill-frame, and the new model is more consistent with the seismic damage of frame structure in Lushan earthquake. The model has reference significance to seismic analysis of frame structure.

    masonry infill panel; reinforced concrete frame structure; equivalent strut mode; shear strength; stiffness

    10.11817/j.issn.1672?7207.2017.02.023

    TU31

    A

    1672?7207(2017)02?0440?08

    2016?03?15;

    2016?06?10

    國家自然科學(xué)基金資助項(xiàng)目(91315301);地震行業(yè)科研專項(xiàng)(201508026-01,201508010-06)(Project(91315301) supported by the National Natural Science Foundation of China; Projects(201508026-01, 201508010-06) supported by Earthquake Industry Research)

    鄧宏宇,博士,講師,從事結(jié)構(gòu)抗震研究;E-mail:dengblue@163.com

    猜你喜歡
    梁柱對(duì)角砌體
    砌體墻上安裝摩擦型阻尼器施工技術(shù)探討
    豎向開槽砌體墻燃?xì)獗▌?dòng)力響應(yīng)及加固
    大截面鋼管勁性柱梁柱節(jié)點(diǎn)施工技術(shù)的要點(diǎn)分析
    斜腿剛構(gòu)拱橋貝雷梁柱式支架的安全性驗(yàn)算
    擬對(duì)角擴(kuò)張Cuntz半群的某些性質(zhì)
    鋼框架腹板雙角鋼連接梁柱子結(jié)構(gòu)抗倒塌性能分析
    論建筑工程中砌體結(jié)構(gòu)現(xiàn)狀及前景
    采高對(duì)砌體梁關(guān)鍵層位置的影響分析
    低屈服點(diǎn)鋼板剪力墻梁柱剛接與梁柱鉸接對(duì)比分析
    非奇異塊α1對(duì)角占優(yōu)矩陣新的實(shí)用簡捷判據(jù)
    80岁老熟妇乱子伦牲交| 亚洲精品中文字幕在线视频 | 人体艺术视频欧美日本| 亚洲一级一片aⅴ在线观看| 六月丁香七月| 国产91av在线免费观看| 日韩,欧美,国产一区二区三区| 亚洲av欧美aⅴ国产| 18禁裸乳无遮挡动漫免费视频| 国产av国产精品国产| 妹子高潮喷水视频| 国产伦在线观看视频一区| 男女啪啪激烈高潮av片| 我要看日韩黄色一级片| 一级黄片播放器| 日韩视频在线欧美| 亚洲欧美一区二区三区黑人 | 一级爰片在线观看| 国产成人精品久久久久久| 我的女老师完整版在线观看| 热re99久久国产66热| 纯流量卡能插随身wifi吗| 国产伦在线观看视频一区| 汤姆久久久久久久影院中文字幕| 欧美国产精品一级二级三级 | av在线观看视频网站免费| 久久人人爽人人爽人人片va| 国产成人aa在线观看| 午夜影院在线不卡| 欧美精品一区二区大全| 成人影院久久| 99热全是精品| 一个人看视频在线观看www免费| 午夜日本视频在线| www.色视频.com| xxx大片免费视频| 精品久久久久久久久av| 一级a做视频免费观看| 国产精品三级大全| 高清av免费在线| 偷拍熟女少妇极品色| 性色avwww在线观看| 蜜桃久久精品国产亚洲av| 十八禁网站网址无遮挡 | 如何舔出高潮| 欧美老熟妇乱子伦牲交| 成人亚洲欧美一区二区av| 免费av不卡在线播放| 只有这里有精品99| 69精品国产乱码久久久| 成人国产麻豆网| 日韩不卡一区二区三区视频在线| 桃花免费在线播放| 亚洲精品色激情综合| 欧美另类一区| 精品少妇内射三级| 亚洲av免费高清在线观看| 能在线免费看毛片的网站| 乱系列少妇在线播放| 亚洲国产成人一精品久久久| 免费人成在线观看视频色| 性高湖久久久久久久久免费观看| 寂寞人妻少妇视频99o| 欧美日韩精品成人综合77777| 日韩不卡一区二区三区视频在线| 日韩精品有码人妻一区| 九色成人免费人妻av| 日本av手机在线免费观看| 中国国产av一级| 一本色道久久久久久精品综合| 亚洲激情五月婷婷啪啪| 大陆偷拍与自拍| 汤姆久久久久久久影院中文字幕| 视频区图区小说| 国产亚洲午夜精品一区二区久久| 成人特级av手机在线观看| 91久久精品国产一区二区成人| 国模一区二区三区四区视频| 国产高清国产精品国产三级| 亚洲一级一片aⅴ在线观看| 黄片无遮挡物在线观看| 交换朋友夫妻互换小说| 免费观看a级毛片全部| 亚洲激情五月婷婷啪啪| 一级爰片在线观看| 国产精品久久久久久精品电影小说| 在线免费观看不下载黄p国产| av女优亚洲男人天堂| 亚洲欧美精品专区久久| 在线观看一区二区三区激情| 亚洲欧美日韩卡通动漫| 亚洲精品aⅴ在线观看| 2018国产大陆天天弄谢| 老女人水多毛片| av天堂久久9| 亚洲国产精品一区二区三区在线| 欧美3d第一页| 亚州av有码| 国产在线视频一区二区| 亚洲国产色片| 少妇丰满av| h视频一区二区三区| 精品久久久精品久久久| 亚洲精品日韩av片在线观看| 一级爰片在线观看| 久久99热这里只频精品6学生| 99久久精品国产国产毛片| 精品午夜福利在线看| 国产毛片在线视频| 国产精品久久久久久精品电影小说| 精华霜和精华液先用哪个| av专区在线播放| 91精品国产国语对白视频| av在线观看视频网站免费| 涩涩av久久男人的天堂| 亚洲欧美日韩卡通动漫| 在线观看av片永久免费下载| av线在线观看网站| 嫩草影院新地址| 美女大奶头黄色视频| 伊人久久国产一区二区| 91在线精品国自产拍蜜月| 国产成人精品无人区| 国产日韩欧美亚洲二区| 亚洲,一卡二卡三卡| 国产精品福利在线免费观看| 日韩一区二区视频免费看| 中文字幕精品免费在线观看视频 | 卡戴珊不雅视频在线播放| 国产亚洲5aaaaa淫片| 男人添女人高潮全过程视频| 黄色配什么色好看| 制服丝袜香蕉在线| 又爽又黄a免费视频| 最近中文字幕高清免费大全6| 91久久精品电影网| 97在线人人人人妻| 美女内射精品一级片tv| 最黄视频免费看| 欧美人与善性xxx| 国产亚洲一区二区精品| 一级毛片电影观看| 成人毛片60女人毛片免费| 国产一区有黄有色的免费视频| 18禁在线无遮挡免费观看视频| 乱系列少妇在线播放| 国产乱人偷精品视频| 大片免费播放器 马上看| 26uuu在线亚洲综合色| 精品少妇黑人巨大在线播放| 久久精品久久久久久久性| 日韩强制内射视频| 久久久亚洲精品成人影院| 大又大粗又爽又黄少妇毛片口| 久久午夜福利片| 26uuu在线亚洲综合色| 久久久久国产网址| 午夜日本视频在线| 纵有疾风起免费观看全集完整版| 国产亚洲av片在线观看秒播厂| 一二三四中文在线观看免费高清| 国产精品成人在线| 午夜福利,免费看| 国产有黄有色有爽视频| 一区在线观看完整版| av一本久久久久| 国产片特级美女逼逼视频| 青春草国产在线视频| 久久精品熟女亚洲av麻豆精品| 亚洲中文av在线| 国产亚洲av片在线观看秒播厂| 午夜福利视频精品| 草草在线视频免费看| 美女脱内裤让男人舔精品视频| 黑人巨大精品欧美一区二区蜜桃 | 乱码一卡2卡4卡精品| 人人妻人人爽人人添夜夜欢视频 | 日韩一本色道免费dvd| 亚洲精品视频女| 婷婷色麻豆天堂久久| 国产精品福利在线免费观看| 欧美少妇被猛烈插入视频| 国产亚洲精品久久久com| 成人二区视频| 精品亚洲乱码少妇综合久久| 一本一本综合久久| 中国国产av一级| 亚洲欧美成人综合另类久久久| 亚洲成色77777| 国产日韩欧美在线精品| 在线免费观看不下载黄p国产| 精品一区在线观看国产| av免费观看日本| 久久久久久久久久人人人人人人| 国产乱人偷精品视频| 亚洲精品日本国产第一区| 国产91av在线免费观看| 国产伦理片在线播放av一区| 欧美日本中文国产一区发布| av天堂中文字幕网| 精品国产一区二区三区久久久樱花| 久久99一区二区三区| 国产成人精品婷婷| 色婷婷久久久亚洲欧美| 99热这里只有是精品50| 日韩欧美一区视频在线观看 | 午夜久久久在线观看| 欧美bdsm另类| 国产精品久久久久久久电影| 男女无遮挡免费网站观看| 免费看日本二区| 久久国内精品自在自线图片| 乱码一卡2卡4卡精品| 国产av一区二区精品久久| 中文字幕人妻熟人妻熟丝袜美| 99久国产av精品国产电影| 成年人免费黄色播放视频 | 亚洲成人一二三区av| 亚洲综合精品二区| 国产高清国产精品国产三级| 亚洲av成人精品一二三区| 中文天堂在线官网| 免费大片黄手机在线观看| 亚洲第一av免费看| 免费看av在线观看网站| 人妻人人澡人人爽人人| 久久ye,这里只有精品| 极品人妻少妇av视频| 日韩av在线免费看完整版不卡| 亚洲久久久国产精品| 国产精品麻豆人妻色哟哟久久| 国产黄色视频一区二区在线观看| 国产综合精华液| 亚洲自偷自拍三级| 亚洲精品久久久久久婷婷小说| 久久鲁丝午夜福利片| 亚洲第一av免费看| 嘟嘟电影网在线观看| 免费观看的影片在线观看| 亚洲精品国产色婷婷电影| 精品久久久久久电影网| 亚洲av中文av极速乱| 国产av精品麻豆| 色5月婷婷丁香| 中文字幕人妻熟人妻熟丝袜美| 欧美精品亚洲一区二区| 亚洲高清免费不卡视频| 伦理电影免费视频| 女人精品久久久久毛片| 国产精品成人在线| 人妻夜夜爽99麻豆av| 伦精品一区二区三区| 国产色婷婷99| 国产无遮挡羞羞视频在线观看| 最后的刺客免费高清国语| 91久久精品国产一区二区三区| 亚洲欧美成人综合另类久久久| 纵有疾风起免费观看全集完整版| 十八禁网站网址无遮挡 | 亚洲精品视频女| 夜夜看夜夜爽夜夜摸| 亚洲久久久国产精品| 精品人妻熟女av久视频| 国产男女内射视频| 国产精品偷伦视频观看了| 我的老师免费观看完整版| 国产在线视频一区二区| 午夜视频国产福利| 伊人久久国产一区二区| 如何舔出高潮| 国产 精品1| 大香蕉97超碰在线| 午夜福利视频精品| 天堂8中文在线网| 亚洲欧洲日产国产| 爱豆传媒免费全集在线观看| 日韩精品有码人妻一区| 久久99蜜桃精品久久| 伊人亚洲综合成人网| 99九九线精品视频在线观看视频| 在线亚洲精品国产二区图片欧美 | 久久精品久久久久久久性| 人妻少妇偷人精品九色| 美女cb高潮喷水在线观看| 日本爱情动作片www.在线观看| 九色成人免费人妻av| 日韩不卡一区二区三区视频在线| 亚洲精品视频女| 男的添女的下面高潮视频| 欧美三级亚洲精品| 波野结衣二区三区在线| 人人妻人人爽人人添夜夜欢视频 | 一个人看视频在线观看www免费| 成人综合一区亚洲| 在线观看一区二区三区激情| 精品亚洲乱码少妇综合久久| 亚洲国产成人一精品久久久| 热re99久久精品国产66热6| 免费黄色在线免费观看| 亚洲国产色片| 性色avwww在线观看| 中文字幕av电影在线播放| 久久久久久久精品精品| 亚洲情色 制服丝袜| 久久久精品94久久精品| 欧美精品国产亚洲| 日本色播在线视频| 国产精品偷伦视频观看了| 国产日韩一区二区三区精品不卡 | 日本vs欧美在线观看视频 | 久久久久久人妻| 久久av网站| 一区在线观看完整版| 成人国产麻豆网| 日韩精品免费视频一区二区三区 | 久久婷婷青草| 亚洲精品日韩av片在线观看| 噜噜噜噜噜久久久久久91| av国产精品久久久久影院| 高清不卡的av网站| 毛片一级片免费看久久久久| 亚洲高清免费不卡视频| 日韩欧美精品免费久久| 婷婷色综合www| 女性被躁到高潮视频| 国产淫片久久久久久久久| 精品人妻熟女av久视频| 久久久a久久爽久久v久久| 永久免费av网站大全| 91久久精品国产一区二区三区| 精品久久久精品久久久| 99热这里只有精品一区| av视频免费观看在线观看| 欧美xxⅹ黑人| av免费在线看不卡| 六月丁香七月| 久久精品国产a三级三级三级| 大陆偷拍与自拍| 精品人妻偷拍中文字幕| 汤姆久久久久久久影院中文字幕| 婷婷色麻豆天堂久久| 哪个播放器可以免费观看大片| 免费观看a级毛片全部| 日本wwww免费看| 精品午夜福利在线看| 国产亚洲欧美精品永久| 男女国产视频网站| 人人澡人人妻人| av卡一久久| 2021少妇久久久久久久久久久| 99久国产av精品国产电影| 少妇 在线观看| 如何舔出高潮| 亚洲美女搞黄在线观看| 亚洲av成人精品一区久久| 中文字幕久久专区| 亚洲无线观看免费| 精品少妇黑人巨大在线播放| h日本视频在线播放| 欧美最新免费一区二区三区| 日本-黄色视频高清免费观看| 国产精品一区二区性色av| 精品人妻偷拍中文字幕| 午夜日本视频在线| 最黄视频免费看| 亚州av有码| 男女边吃奶边做爰视频| 校园人妻丝袜中文字幕| 性色avwww在线观看| 国产精品嫩草影院av在线观看| 国产精品一区www在线观看| 免费看av在线观看网站| 夜夜看夜夜爽夜夜摸| 精品久久国产蜜桃| 亚洲欧美日韩卡通动漫| 丰满少妇做爰视频| 亚洲人与动物交配视频| 国产无遮挡羞羞视频在线观看| 亚洲国产欧美日韩在线播放 | 午夜日本视频在线| 六月丁香七月| 青春草视频在线免费观看| 欧美人与善性xxx| 一级毛片黄色毛片免费观看视频| 美女大奶头黄色视频| 最黄视频免费看| 久久久久视频综合| 综合色丁香网| 热99国产精品久久久久久7| 熟女人妻精品中文字幕| 欧美丝袜亚洲另类| 国产色爽女视频免费观看| 秋霞伦理黄片| 免费av不卡在线播放| 国产免费一级a男人的天堂| 亚洲av国产av综合av卡| 亚洲国产精品专区欧美| 久久ye,这里只有精品| av天堂久久9| 丰满少妇做爰视频| 国产精品一二三区在线看| 中文在线观看免费www的网站| 久久精品熟女亚洲av麻豆精品| 麻豆成人av视频| 久久99一区二区三区| 又大又黄又爽视频免费| 一级毛片电影观看| 人妻人人澡人人爽人人| 热re99久久精品国产66热6| 午夜影院在线不卡| 九草在线视频观看| 伦精品一区二区三区| √禁漫天堂资源中文www| 能在线免费看毛片的网站| 欧美丝袜亚洲另类| 久久国产精品男人的天堂亚洲 | 亚洲真实伦在线观看| 日韩中文字幕视频在线看片| 国精品久久久久久国模美| 熟女人妻精品中文字幕| 亚洲av不卡在线观看| 亚洲av成人精品一区久久| 高清视频免费观看一区二区| 成人毛片60女人毛片免费| 国产精品秋霞免费鲁丝片| 欧美性感艳星| 我要看黄色一级片免费的| 午夜精品国产一区二区电影| 蜜桃在线观看..| 久久99一区二区三区| 精品一区二区免费观看| 亚洲欧洲精品一区二区精品久久久 | 高清视频免费观看一区二区| 午夜福利网站1000一区二区三区| 精品久久久久久久久亚洲| 久久ye,这里只有精品| 99久国产av精品国产电影| 99久久人妻综合| 成年人免费黄色播放视频 | 天美传媒精品一区二区| 日日啪夜夜爽| 国内少妇人妻偷人精品xxx网站| 午夜91福利影院| 精品99又大又爽又粗少妇毛片| 久久久精品免费免费高清| 久久人妻熟女aⅴ| 亚洲av.av天堂| 最近中文字幕高清免费大全6| 人妻少妇偷人精品九色| 国产成人一区二区在线| 我的女老师完整版在线观看| 免费看av在线观看网站| 亚洲精华国产精华液的使用体验| 日本免费在线观看一区| 亚洲国产av新网站| 亚洲国产毛片av蜜桃av| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 亚洲成人一二三区av| 你懂的网址亚洲精品在线观看| 色哟哟·www| 国产精品熟女久久久久浪| 国产成人精品一,二区| 成年美女黄网站色视频大全免费 | 国产av码专区亚洲av| 亚洲国产av新网站| 嫩草影院入口| 亚洲国产成人一精品久久久| 亚洲国产日韩一区二区| 亚洲熟女精品中文字幕| 日韩成人av中文字幕在线观看| 亚洲精品第二区| 蜜桃在线观看..| 久久99精品国语久久久| 国产日韩一区二区三区精品不卡 | www.色视频.com| 啦啦啦啦在线视频资源| 亚洲欧美一区二区三区黑人 | 精品一区在线观看国产| 99热这里只有精品一区| 国产成人精品久久久久久| 嫩草影院新地址| 亚洲精品一二三| 免费观看a级毛片全部| 黄色毛片三级朝国网站 | 日本免费在线观看一区| 777米奇影视久久| 国产美女午夜福利| 日韩欧美一区视频在线观看 | 美女福利国产在线| 午夜精品国产一区二区电影| 人人妻人人澡人人爽人人夜夜| 国产极品粉嫩免费观看在线 | 日本与韩国留学比较| 国产精品久久久久久精品电影小说| 亚洲欧美一区二区三区黑人 | 亚洲国产成人一精品久久久| 免费观看av网站的网址| 国产男女内射视频| 少妇熟女欧美另类| 99国产精品免费福利视频| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| av网站免费在线观看视频| 麻豆成人午夜福利视频| 寂寞人妻少妇视频99o| av天堂久久9| 18禁动态无遮挡网站| av天堂中文字幕网| 亚洲综合色惰| 高清欧美精品videossex| 精品人妻一区二区三区麻豆| 成人二区视频| 波野结衣二区三区在线| 亚洲av中文av极速乱| 黄色怎么调成土黄色| 亚洲图色成人| 寂寞人妻少妇视频99o| 亚洲欧洲精品一区二区精品久久久 | 交换朋友夫妻互换小说| 美女xxoo啪啪120秒动态图| 99热这里只有精品一区| 久久女婷五月综合色啪小说| 久久久久视频综合| 一级黄片播放器| 亚洲在久久综合| 国产一区二区在线观看av| 女性被躁到高潮视频| 亚洲性久久影院| 亚洲国产精品专区欧美| av女优亚洲男人天堂| 国产又色又爽无遮挡免| 亚洲精品日韩av片在线观看| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 亚洲欧美精品专区久久| 成人无遮挡网站| 一级爰片在线观看| 国产精品偷伦视频观看了| 美女福利国产在线| 免费人成在线观看视频色| av免费在线看不卡| 国产成人精品一,二区| 日韩中文字幕视频在线看片| 赤兔流量卡办理| 国产成人精品婷婷| 99久久精品一区二区三区| 有码 亚洲区| 99久久精品一区二区三区| 精品久久国产蜜桃| 80岁老熟妇乱子伦牲交| 又大又黄又爽视频免费| 久久久久精品性色| 国产成人一区二区在线| 97超视频在线观看视频| 少妇人妻精品综合一区二区| 日韩熟女老妇一区二区性免费视频| 97超视频在线观看视频| 少妇 在线观看| 亚洲成人av在线免费| 搡老乐熟女国产| 少妇精品久久久久久久| 国产午夜精品久久久久久一区二区三区| 男男h啪啪无遮挡| 九草在线视频观看| 亚洲国产最新在线播放| 免费观看在线日韩| 国产成人aa在线观看| 伦理电影大哥的女人| 国产亚洲欧美精品永久| 蜜桃久久精品国产亚洲av| 亚洲精品aⅴ在线观看| 男男h啪啪无遮挡| a级片在线免费高清观看视频| 最新中文字幕久久久久| 亚洲成人一二三区av| 最新中文字幕久久久久| 高清在线视频一区二区三区| 一区二区av电影网| 欧美丝袜亚洲另类| 交换朋友夫妻互换小说| av国产久精品久网站免费入址| 赤兔流量卡办理| 国产成人91sexporn| 亚洲av综合色区一区| 天天躁夜夜躁狠狠久久av| 精品久久久精品久久久| 欧美区成人在线视频| 高清欧美精品videossex| 日本黄大片高清| 国产精品一区www在线观看| 一区二区三区精品91| 日韩一区二区三区影片| 久久99一区二区三区| 欧美日韩国产mv在线观看视频| 亚洲精华国产精华液的使用体验| 少妇人妻久久综合中文| 亚洲成人一二三区av| 亚洲一区二区三区欧美精品| 老司机影院成人| 九九在线视频观看精品| 人人妻人人添人人爽欧美一区卜| 男人添女人高潮全过程视频| 女人精品久久久久毛片| 自线自在国产av| 69精品国产乱码久久久| a级毛色黄片| 亚洲精品乱码久久久久久按摩| 一区二区三区乱码不卡18| 看十八女毛片水多多多| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧洲日产国产| 国产免费又黄又爽又色| 一本—道久久a久久精品蜜桃钙片| 特大巨黑吊av在线直播|