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

    方形中空夾層鋼管超高性能鋼纖維混凝土柱抗爆性能數(shù)值模擬與實(shí)驗(yàn)驗(yàn)證

    2017-02-14 09:26:35徐慎春劉中憲吳成清
    振動(dòng)與沖擊 2017年1期
    關(guān)鍵詞:抗爆性內(nèi)層外層

    徐慎春, 劉中憲, 吳成清

    (1.天津大學(xué) 建工學(xué)院,天津 300372; 2.天津城建大學(xué) 天津市土木建筑結(jié)構(gòu)防護(hù)與加固重點(diǎn)實(shí)驗(yàn)室,天津 300384;3.天津城建大學(xué)與悉尼科技大學(xué)防災(zāi)減災(zāi)聯(lián)合研究中心,天津 300384)

    方形中空夾層鋼管超高性能鋼纖維混凝土柱抗爆性能數(shù)值模擬與實(shí)驗(yàn)驗(yàn)證

    徐慎春1, 劉中憲2,3, 吳成清2,3

    (1.天津大學(xué) 建工學(xué)院,天津 300372; 2.天津城建大學(xué) 天津市土木建筑結(jié)構(gòu)防護(hù)與加固重點(diǎn)實(shí)驗(yàn)室,天津 300384;3.天津城建大學(xué)與悉尼科技大學(xué)防災(zāi)減災(zāi)聯(lián)合研究中心,天津 300384)

    建立了爆炸荷載作用下方形中空夾層鋼管超高性能鋼纖維混凝土(Ultra-High Performance Steel Fiber Reinforced Concrete Filled Double Skin Steel Tube, UHPSFRCFDST)柱動(dòng)態(tài)響應(yīng)及其損傷破壞三維有限元數(shù)值模型。首先通過(guò)模擬結(jié)果與爆炸破壞試驗(yàn)結(jié)果的對(duì)比分析,驗(yàn)證了數(shù)值模型和計(jì)算方法的有效性。進(jìn)而運(yùn)用參數(shù)化分析方法,研究了空心率、含鋼率、內(nèi)、外層鋼管厚度及其強(qiáng)度等關(guān)鍵參數(shù)對(duì)UHPSFRCFDST柱抗爆性能的影響。研究結(jié)果表明,UHPSFRCFDST柱具有優(yōu)越的抗爆性能,所建立的三維有限元模型能夠有效地分析UHPSFRCFDST柱在爆炸荷載作用下的動(dòng)態(tài)響應(yīng)及其損傷破壞;在一定范圍內(nèi)減小空心率及提高外層鋼管強(qiáng)度可有效提升UHPSFRCFDST柱抗爆性能; 提高含鋼率、減小內(nèi)、外層鋼管高厚比均能夠顯著提升UHPSFRCDST柱抗爆性能;內(nèi)層鋼管強(qiáng)度對(duì)UHPSFRCFDST柱的抗爆性能影響并不明顯。

    中空夾層鋼管超高性能鋼纖維混凝土柱;爆炸荷載;動(dòng)態(tài)響應(yīng);破壞模式;有限元模型

    實(shí)心鋼管混凝土通常被用作軸心受壓或荷載偏心較小的壓彎構(gòu)件當(dāng)中[1-3],當(dāng)長(zhǎng)細(xì)比或荷載偏心較大時(shí),構(gòu)件抗彎剛度在很大程度上決定了其承載力,而靠近截面形心軸的混凝土僅提供了小部分的抗彎剛度,卻增大大加了構(gòu)件自重,中空夾層鋼管混凝(Concrete-filled Double Skin Steel Tube, CFDST)土可顯著改善上述缺陷。因此,該結(jié)構(gòu)形式在大長(zhǎng)細(xì)比或大偏心率構(gòu)件的應(yīng)用方面具有潛在優(yōu)勢(shì)[4]。

    然而,近些年來(lái),一方面世界范圍內(nèi)爆炸恐怖襲擊事件頻繁發(fā)生。另一方面,隨著我國(guó)經(jīng)濟(jì)的發(fā)展,人們對(duì)化工產(chǎn)品的需求不斷增加,加之人們?cè)谏a(chǎn)、生活中的疏忽等意外因素,造成了諸如化工品泄露爆炸、粉塵爆炸和燃油泄露爆炸等意外爆炸事件也時(shí)有發(fā)生。鑒于CFDST柱通常被用作橋墩、海洋平臺(tái)支架柱、建筑中大直徑柱或其他高聳結(jié)構(gòu)柱等大長(zhǎng)細(xì)比或大偏心率構(gòu)件中,一旦其遭受爆炸荷載發(fā)生破壞,將會(huì)導(dǎo)致極為嚴(yán)重的生命、財(cái)產(chǎn)損失。因此,提高CFDST構(gòu)件的抗爆性能便成為了一個(gè)值得研究的課題。

    中空夾層超高性能鋼纖維混凝土(Ultra-High Performance Steel Fiber Reinforced Concrete-filled Double Skin Steel Tubes, UHPSFRCFDST)是在兩個(gè)同心放置的鋼管之間灌注超高性能鋼纖維混凝土而形成的新型鋼管混凝土組合結(jié)構(gòu)形式。研究表明,超高性能鋼纖維混凝土(Ultra-High Performance Steel Fiber Reinforced Concrete,UHPSFRC)具有高強(qiáng)度、高韌性、高耐久性等卓越的力學(xué)性能,因而能夠更加有效地抵抗爆炸、沖擊等高應(yīng)變率荷載[5-6]。已有少量研究表明超高性能鋼纖維混凝土結(jié)構(gòu)的抗爆性能遠(yuǎn)遠(yuǎn)優(yōu)于普通混凝土結(jié)構(gòu)[7]。因此,可以預(yù)見(jiàn),較之于傳統(tǒng)中空夾層鋼管混凝土結(jié)構(gòu),中空夾層鋼管超高性能鋼纖維混凝土結(jié)構(gòu)的抗爆性能會(huì)有較大提升。

    國(guó)內(nèi)外已有針對(duì)普通鋼管混凝土及復(fù)式空心鋼管混凝土抗爆性能的研究。例如,國(guó)內(nèi),孫珊珊等[8]通過(guò)理論分析、試驗(yàn)研究與數(shù)值模擬相結(jié)合的方法對(duì)爆炸荷載作用下鋼管混凝土柱動(dòng)態(tài)響應(yīng)、破壞機(jī)理等進(jìn)行了研究;崔瑩等[9]對(duì)復(fù)式空心鋼管混凝土在爆炸荷載作用下的動(dòng)態(tài)響應(yīng)及其損傷破壞機(jī)理進(jìn)行了試驗(yàn)研究及數(shù)值模擬;李國(guó)強(qiáng)等通過(guò)12個(gè)鋼管混凝土柱試件的現(xiàn)場(chǎng)爆炸試驗(yàn),研究了沖擊壓力波在柱表面的分布形式以及相關(guān)關(guān)鍵參數(shù)對(duì)柱動(dòng)態(tài)響應(yīng)的影響。國(guó)外,ZHANG等[11]結(jié)合現(xiàn)場(chǎng)爆炸試驗(yàn)以及數(shù)值模擬技術(shù)研究了鋼管混凝土在爆炸荷載作用下的動(dòng)態(tài)響應(yīng)及其耗能機(jī)理;Jama等[12]采用數(shù)值模擬方法研究了方鋼管混凝土梁在橫向爆炸荷載作用下的動(dòng)態(tài)響應(yīng)。綜上所述,現(xiàn)有報(bào)道鮮有提及鋼管超高性能混凝土柱抗爆性能的研究,而針對(duì)中空夾層鋼管超高性能混凝土柱抗爆性能的研究更是未見(jiàn)報(bào)道。

    本文針對(duì)UHPSFRCFDST柱的抗爆性能開(kāi)展研究,運(yùn)用顯式有限元分析軟件LS-DYNA建立了爆炸荷載作用下UPHSFRCFDST柱動(dòng)態(tài)響應(yīng)及其損傷破壞的三維有限元模型,通過(guò)模擬結(jié)果與野外爆炸破壞試驗(yàn)數(shù)據(jù)對(duì)比驗(yàn)證了模型和分析方法的有效性。進(jìn)而研究了空心率、含鋼率、內(nèi)、外層鋼管厚度及鋼管強(qiáng)度等關(guān)鍵參數(shù)對(duì)UPHSFRCFDST柱抗爆性能的影響,以期能夠?yàn)閁PHSFRCFDST結(jié)構(gòu)的抗爆設(shè)計(jì)提供一定的理論依據(jù)。

    1 數(shù)值模擬方法

    通過(guò)顯式有限元軟件LS-DYNA建立了3D有限元模型,進(jìn)而運(yùn)用此模型對(duì)UHPSFRCFDST柱爆炸荷載作用下動(dòng)態(tài)響應(yīng)及損傷破壞進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)驗(yàn)證。

    1.1 有限元模型

    建立的3D有限元模型如圖1所示。通過(guò)網(wǎng)格收斂性分析確定網(wǎng)格尺寸為5 mm。試件兩端為簡(jiǎn)支邊界,支撐與試件之間采用自動(dòng)單面接觸算法。鋼管與超高性能混凝土單元在接觸面上節(jié)點(diǎn)固結(jié),不考慮兩者之間的粘結(jié)滑移。爆炸荷載采用LS-DYNA內(nèi)置的ConWep程序施加[13-14]。

    1.2 材料本構(gòu)模型

    內(nèi)、外層鋼管均采用隨動(dòng)強(qiáng)化本構(gòu)模型[13-14],此模型能夠高效、準(zhǔn)確地模擬鋼筋、鋼管等力學(xué)性能,且已被大量研究驗(yàn)證[15-16]。隨動(dòng)強(qiáng)化模型假設(shè)在塑性變形過(guò)程中,屈服面的大小和形狀不變,僅在塑性應(yīng)變方向上平動(dòng),亦考慮了包辛格效應(yīng)[17]和應(yīng)變率效應(yīng),其中應(yīng)變率效應(yīng)采用Cowper-Symonds模型,如下式:

    (1)

    UHPSFRC采用K&C本構(gòu)模型,此模型考慮了混凝土應(yīng)變硬化、軟化特性、應(yīng)變率效應(yīng)以及高靜水壓的影響,已被廣泛應(yīng)用于混凝土結(jié)構(gòu)抗爆分析[18-20]。其中,師燕超采用此模型對(duì)某典型鋼筋混凝土柱爆炸荷載作用下的動(dòng)態(tài)響應(yīng)及損傷破壞模式進(jìn)行了研究,MAO等則運(yùn)用此模型對(duì)某超高性能鋼纖維混凝土板在爆炸荷載作用下的動(dòng)態(tài)響應(yīng)及損傷破壞進(jìn)行了數(shù)值模擬,并給出了適應(yīng)于抗壓強(qiáng)度為170 MPa 的UHPSFRC K&C本構(gòu)模型參數(shù)。

    在有限元軟件LS-DYNA中,僅輸入混凝土單軸抗壓強(qiáng)度,程序便可自動(dòng)計(jì)算適應(yīng)于此強(qiáng)度的本構(gòu)參數(shù)。然而,這些參數(shù)是以大量普通混凝土材料性能試驗(yàn)為依據(jù)生成的,并不完全適用于UHPSFRC[21]。因此,為真實(shí)反應(yīng)UHPSFRC宏觀力學(xué)性能,本文在MAO所提供的K&C本構(gòu)模型參數(shù)的基礎(chǔ)上,對(duì)UHPSFRC單軸抗壓性能進(jìn)行了數(shù)值模擬,通過(guò)與試驗(yàn)數(shù)據(jù)的對(duì)比分析對(duì)合理地調(diào)整了文獻(xiàn)[20]的參數(shù)。圖2比較了試驗(yàn)和數(shù)值模擬所獲取的單軸抗壓應(yīng)力應(yīng)變曲線,從中可以看出二者吻合較好。最終所采用的K&C本構(gòu)模型關(guān)鍵參數(shù)列于表1。

    圖2 UHSFRC單軸抗壓應(yīng)力應(yīng)變曲線Fig.2 Uniaxial compressive stress-strain curve

    關(guān)鍵參數(shù)數(shù)值a0y3.3480E+007a1y6.250E-001a2y1.717E-009a0m4.4340E+007a1m4.463E-001a2m5.387E-010a1r4.417E-001a2r7.887E-010b2-2

    現(xiàn)有研究表明,混凝土的應(yīng)變率效應(yīng)隨其單軸抗壓強(qiáng)度的不同而不同。NGO[22]提出了適用于單軸抗壓強(qiáng)度為32~160 MPa的混凝土的動(dòng)態(tài)增大系數(shù)模型,如下式:

    (2)

    然而,拉伸狀態(tài)下UHPSFRC的應(yīng)變率效應(yīng)仍在研究當(dāng)中,尚未形成統(tǒng)一結(jié)論。因此在考慮UHSFRC拉伸動(dòng)態(tài)應(yīng)變率效應(yīng)時(shí)忽略了混凝土強(qiáng)度的影響,仍然選用CEB模型[24],如下式:

    (3)

    2 數(shù)值模型實(shí)驗(yàn)驗(yàn)證

    采用UHPSFRCFDST柱的野外爆炸破壞試驗(yàn),驗(yàn)證本文所建立的3D有限元模型的有效性。試件長(zhǎng)2 500 mm,截面外徑200 mm,內(nèi)徑100 mm,內(nèi)、外層鋼管均為5 mm厚Q235無(wú)縫鋼管,截面如圖3(a)所示,圖3(b)給出了柱結(jié)構(gòu)示意圖。

    本試驗(yàn)在中國(guó)人民解放軍理工大學(xué)某靶場(chǎng)內(nèi)進(jìn)行,試驗(yàn)設(shè)置如圖3(c),3(d)及圖4所示。為放置試件及試驗(yàn)設(shè)備,挖掘一個(gè)3 m×0.4 m×1.5 m的試驗(yàn)坑洞,并使用20 mm厚的鋼板覆蓋試驗(yàn)坑洞表面以保護(hù)試驗(yàn)器材。爆源為48 kg 圓形TNT炸藥,其與試件迎爆面距離1.5 m,試件兩端為簡(jiǎn)支邊界。為記錄試件豎向位移動(dòng)態(tài)響應(yīng),在試件背爆面放置3個(gè)位移傳感器(LVDTs),間距為380 mm(從柱端到柱中編號(hào)分別為L(zhǎng)T1~LT3),表2給出了位移傳感器具體參數(shù)。同時(shí),在與UHPSFRCFDST柱相鄰的支撐鋼架上固定安裝美國(guó)PCB公司M113A32型爆炸壓力傳感器以精確測(cè)量反射超壓Pr(t),具體參數(shù)列于表2,3。

    圖3 爆炸破壞試驗(yàn)設(shè)置Fig.3 Diagram of blast test

    下節(jié)會(huì)對(duì)爆炸破壞試驗(yàn)數(shù)據(jù)與數(shù)值模擬結(jié)果進(jìn)行對(duì)比分析,以驗(yàn)證3D有限元模型的有效性,故此處不再贅述試驗(yàn)結(jié)果。

    圖4 爆炸破壞試驗(yàn)示意圖Fig.4 Diagram of blast test

    型號(hào)MEERKTC通用型行程/mm125機(jī)械行程/mm行程+7mm線性精度/%FS0.08%電阻/kΩ5解析度Infinite無(wú)斷解析重復(fù)性精度/mm0.01使用溫度范圍/℃-60~+150

    表3 爆炸壓力傳感器參數(shù)

    2.1 反射超壓驗(yàn)證

    為準(zhǔn)確分析方形中空夾層鋼管超高性能鋼纖維混凝土柱在爆炸荷載作用下的動(dòng)態(tài)響應(yīng)及其損傷破壞,對(duì)ConWep程序計(jì)算的爆炸荷載進(jìn)行了驗(yàn)證,以調(diào)整相關(guān)爆炸參數(shù),從而確保在數(shù)值模擬過(guò)程中對(duì)試件施加可靠的爆炸荷載。經(jīng)多次試算后,調(diào)整模擬炸藥當(dāng)量為54 kg。

    爆炸波反射超壓時(shí)程曲線如圖5所示,圖5展示了爆炸波超壓ConWep程序計(jì)算值、文獻(xiàn)[25]計(jì)算值以及實(shí)測(cè)值。由圖可以看出,ConWep程序計(jì)算的反射超壓峰值為51.7 MPa,正壓持時(shí)為1.0 ms,文獻(xiàn)[25]計(jì)算值分別為49.6 MPa和0.88 ms,而實(shí)測(cè)值分別為46.4 MPa和0.89 ms,各峰值超壓計(jì)算值與實(shí)測(cè)值誤差很小,分別為11.4%和6.9%。因此,可以說(shuō)明,本次數(shù)值模擬所施加的爆炸荷載是合理有效的。

    此外,盡管理論上而言,文獻(xiàn)[25]和ConWep程序?qū)τ诒_擊波相關(guān)參數(shù)的計(jì)算值應(yīng)該相同,但本文在根據(jù)文獻(xiàn)[25]計(jì)算爆炸沖擊波相關(guān)參數(shù)時(shí),在相關(guān)圖表中采用了人工取值方法,造成了其與ConWep程序計(jì)算值存在微小的誤差。而鑒于文獻(xiàn)[25]計(jì)算值主要用于驗(yàn)證本文所施加荷載的有效性,因此本文忽略了這一微小誤差。

    圖5 爆炸波壓力時(shí)程曲線Fig.5 Time history of blast overpressure

    2.2 動(dòng)態(tài)響應(yīng)驗(yàn)證

    試驗(yàn)中位移計(jì)L3未能獲得試件跨中豎向位移時(shí)程曲線,但L2依然有效測(cè)得了距跨中380 mm處測(cè)點(diǎn)豎向位移時(shí)程曲線,故本文主要通過(guò)對(duì)比此點(diǎn)豎向位移時(shí)程曲線驗(yàn)證3D模型的有效性。L2測(cè)點(diǎn)豎向位移時(shí)程曲線繪于圖6,圖中標(biāo)明了實(shí)測(cè)曲線及模擬曲線。各曲線關(guān)鍵參數(shù)列于表4。由圖及表可以看出,數(shù)值模擬曲線與實(shí)測(cè)曲線吻合很好,驗(yàn)證了所建立3D有限元模型的有效性。

    圖6 試件柱中豎向位移時(shí)程曲線Fig.6 Time history of vertical displacement

    實(shí)測(cè)值模擬值誤差/%峰值位移/mm32.731.24.6殘余位移/mm14.014.10.7

    2.3 破壞模式驗(yàn)證

    本文以Von Misses等效應(yīng)力判定鋼管損傷程度,采用K&C本構(gòu)模型的損傷縮放系數(shù)δ評(píng)估UHPSFRC單元損傷程度,其為等效塑性應(yīng)變的函數(shù),如式(4)、(5)所示。

    δ=2λ/(λ+λm)

    (4)

    (5)

    圖7展示了不同時(shí)刻外層鋼管Von Misses等效應(yīng)力云圖。由圖7可以看出,各時(shí)刻外層鋼管的Von Misses等效應(yīng)力最大值均超過(guò)了Q253鋼材的屈服應(yīng)力,表明外層鋼管發(fā)生了屈服,產(chǎn)生了塑性應(yīng)變。

    圖7 不同時(shí)刻試件外層鋼管 Von Misses應(yīng)力云圖Fig.7 Von Misses stress of outer steel tube at different times

    試件外層鋼管及內(nèi)部UHPSFRC損傷破壞如圖8所示。由圖可以看出,無(wú)論是模擬結(jié)果還是試驗(yàn)結(jié)果,外層鋼管均未出現(xiàn)鼓包屈曲現(xiàn)象,二者在一定程度上較為吻合。盡管應(yīng)變數(shù)據(jù)能更好的說(shuō)明二者的一致性,然而需要注意的是由于試驗(yàn)中所采用的炸藥當(dāng)量較大,爆源距離較近,產(chǎn)生了很大的爆炸沖擊波,沖擊波進(jìn)入空隙后切斷了數(shù)據(jù)采集電纜,導(dǎo)致未能有效測(cè)得外層鋼管應(yīng)變數(shù)據(jù)。因此,未采用應(yīng)變片數(shù)據(jù)對(duì)外層鋼管的損傷破壞進(jìn)行對(duì)比分析。內(nèi)部UHPSFRC在宏觀裂縫開(kāi)展范圍和發(fā)展程度方面與試驗(yàn)結(jié)果能很好地吻合。

    圖8 試件外層鋼管及內(nèi)部UHPSFRC最終形態(tài)Fig.8 Outer steel tube and inner UHPSFRC failure modes

    3 參數(shù)化分析

    上節(jié)已對(duì)所建立的3D有限元模型進(jìn)行了驗(yàn)證,結(jié)果充分表明了此模型能夠有效地模擬爆炸荷載作用下UHPSFRCFDST柱動(dòng)態(tài)響應(yīng)及損傷破壞,動(dòng)態(tài)響應(yīng)主要以距柱中380 mm處豎向位移為指標(biāo),損傷破壞主要以填充UHPSFRC損傷縮放系數(shù)δ為指標(biāo)。本節(jié)將運(yùn)用參數(shù)化分析方法研究不同參數(shù)對(duì)UHPSFRCFDST柱動(dòng)態(tài)響應(yīng)及損傷破壞的影響,主要分析參數(shù)為空心率、含鋼率、內(nèi)、外層鋼管高厚比及其強(qiáng)度等。

    3.1 空心率的影響

    參照相關(guān)規(guī)范,定義空心率(Ψ)為構(gòu)件截面中空心部分的面積與混凝土加空心部分總面積之比。在保持其他條件(內(nèi)、外層鋼管強(qiáng)度及厚度、軸壓比等)不變的情況下,通過(guò)改變?cè)嚰孛婵招牟糠置娣e,分析空心率對(duì)UHPSFRCFDST柱動(dòng)態(tài)響應(yīng)及損傷破壞的影響。

    圖9給出了空心部分面積分別為2 00 mm2、10 000 mm2、22 500 mm2(空心率分比為6.25%、25%、56.25%)時(shí)UHPSFRCFDST柱在相同爆炸荷載作用下距柱中380 mm處位移時(shí)程曲線圖。由圖可以看出,Ψ=56.25%時(shí),試件峰值位移及殘余位移最大;當(dāng)Ψ=25%時(shí),試件峰值位移和殘余位移雖然比Ψ=56.25%小,但比Ψ=6.25%時(shí)要大,差異卻并不明顯。這說(shuō)明,在Ψ<25%范圍內(nèi)減小空心率并不一定能夠提升UHPSFRCFDST柱的抗爆性能,而在25%<Ψ<56.25%范圍內(nèi)減小空心率卻能夠顯著提升UHPSFRCFDST柱的抗爆性能。主要是因?yàn)椋讦?25%范圍內(nèi),雖然Ψ的減小增大了試件中UHPSFRC面積,卻減小了內(nèi)層鋼管的面積,而內(nèi)層鋼管面積的減小對(duì)于試件截面剛度的削弱作用強(qiáng)于UHPSFRC面積增大對(duì)于試件截面剛度的增大作用,而在25%<Ψ<56.25%范圍內(nèi)卻恰恰相反。

    圖9 不同空心率試件距柱中380 mm位移時(shí)程曲線Fig.9 Displacement time history of specimens with different hollow ratios

    圖10給出了不同空心率試件內(nèi)UHPSFRC在相同爆炸荷載作用下的最終破壞形態(tài)。由圖可以看出UHPSFRC損傷隨著空心率的增大而趨于嚴(yán)重,Ψ=6.25%與Ψ=25% UHPSFRC損傷破壞差異并不明顯,而當(dāng)Ψ=56.25%時(shí),試件形成貫通裂縫,損傷破壞最為嚴(yán)重。由此可以說(shuō)明,在Ψ<25%范圍內(nèi)減小空心率并不一定能夠降低試件的損傷破壞,而在25%<Ψ<56.25%范圍內(nèi)減小空心率卻能夠顯著降低試件的損傷破壞。

    圖10 不同空心率試件內(nèi)UHPSFRC損傷δ云圖Fig.10 δ contour of specimens with different hollow ratios

    3.2 含鋼率的影響

    參照相關(guān)規(guī)范,定義含鋼率(αsc)為構(gòu)件截面中鋼管的面積與混凝土的面積之比。在保持其他條件(內(nèi)、外層鋼管強(qiáng)度、空心率、軸壓比等)不變的情況下,通過(guò)改變?cè)嚰孛驿摴苊娣e,分析含鋼率對(duì)UHPSFRCFDST柱動(dòng)態(tài)響應(yīng)及損傷破壞的影響。

    圖11給出了鋼管面積分別為3 552 m2、6 000 mm2、8 512 m2(含鋼率分比為14.3%、25%、36.8%)時(shí)UHPSFRCFDST柱在相同爆炸荷載作用下距柱中380 mm處豎向位移時(shí)程曲線圖。由圖可看出,UHPSFRCFDST柱的峰值位移及殘余位移均隨著含鋼率的增大而減小,說(shuō)明含鋼率對(duì)爆炸荷載作用下UHPSFRCFDST柱的動(dòng)態(tài)響應(yīng)有著顯著影響,隨著含鋼率的增大,其抗爆性能得到了有效提升。

    圖11 不同含鋼率試件柱中位移時(shí)程曲線Fig.11 Displacement time history of specimens with different steel ratios

    圖12給出了不同含鋼率率試件內(nèi)UHPSFRC在相同爆炸荷載作用下的最終破壞形態(tài)。由圖可以看出,隨著含鋼率的增大,無(wú)論是UHPSFRC的損傷范圍還是其損傷程度均均有減小的趨勢(shì)。可以說(shuō)明,提高試件的含鋼率能夠有效降低其損傷破壞,提升其抗爆性能。

    圖12 不同含鋼率試件內(nèi)UHPSFRC損傷δ云圖Fig.12 δ contour of specimens with different steel ratios

    3.3 內(nèi)、外層鋼管高厚比的影響

    參照相關(guān)規(guī)范,定義高厚比為鋼管外徑與鋼管厚度的比值。即內(nèi)層鋼管高厚比(H/ti)為內(nèi)層鋼管高度(H)與厚度(ti)的比值,外層鋼管高厚比(H/to)為外層鋼管高度(H)與外層鋼管厚度(to)的比值。通過(guò)分別改變內(nèi)、外層鋼管的厚度(ti、to)及強(qiáng)度(fi、fo),研究其對(duì)爆炸荷載作用下UHPSFRCFDST柱動(dòng)態(tài)響應(yīng)及損傷破壞的影響。

    圖13給出了內(nèi)層鋼管高厚比分比為 17.1、11、8.1(內(nèi)層鋼管厚度分別為3 mm、5 mm、7 mm)時(shí)UHPSFRCFDST柱在相同爆炸荷載作用下距柱中380 mm處豎向位移時(shí)程曲線圖。由圖可以看出,試件峰值位移及殘余位移均隨內(nèi)層鋼管高厚比的減小而減小,這說(shuō)明減小內(nèi)層鋼管高厚比有助于提升UHPSFRCFDST柱的抗爆性能,主要是因?yàn)閮?nèi)層鋼管高厚比的減小意味著其厚度的增大,進(jìn)而提高了試件整體抗彎剛度,最終提升了試件的抗爆性能。

    圖13 不同內(nèi)層鋼管高厚比試件位移時(shí)程曲線Fig.13 Displacement time history of specimens with different thickness of inner steel tube

    圖14給出了不同內(nèi)層鋼管高厚比試件內(nèi)UHPSFRC在相同爆炸荷載作用下的最終破壞形態(tài)。由圖可以看出,隨著內(nèi)層鋼管高厚比的減小,填充UHPSFRC在損傷范圍上有所擴(kuò)大,但在損傷程度上明顯減小,說(shuō)明減小內(nèi)層鋼管高厚比后,試件損傷分布趨于均勻,體現(xiàn)了試件抗爆耗能能力的提升。

    圖14 不同內(nèi)層鋼管高厚比試件內(nèi)UHPSFRC損傷δ云圖Fig.14 δ contour of specimens with different thickness of inner steel tube

    圖15給出了外層鋼管高厚比分比為33.3、20、14.3(外層鋼管厚度分別為3 mm、5 mm、7 mm)時(shí)UHPSFRCFDST柱在相同爆炸荷載作用下距柱中380 mm處豎向位移時(shí)程曲線圖。由圖可以看出,試件峰值位移及殘余位移均隨外層高厚比的減小而減小,這說(shuō)明減小外層鋼管高厚比有助于提升UHPSFRCFDST柱的抗爆性能,主要是因?yàn)橥鈱愉摴芨吆癖鹊臏p小意味著其厚度的增大,進(jìn)而提高了試件整體抗彎剛度,最終提升了試件的抗爆性能。

    圖15 不同外層鋼管高厚比試件柱中位移時(shí)程曲線Fig.15 Displacement time history of specimens with different thickness of outer steel tube

    圖16給出了不同外層鋼管高厚比試件內(nèi)UHPSFRC在相同爆炸荷載作用下的最終破壞形態(tài)。由圖可以看出,隨著外層鋼管高厚比的減小,試件內(nèi)UHPSFRC損傷沿柱縱向擴(kuò)展更為均勻,截面高度方向損傷程度減弱。可以說(shuō)明,減小UHPSFRCFDST外層鋼管高厚比能夠有效降低試件的損傷破壞,提升試件整體抗爆性能。

    圖16 不同外層鋼管高厚比試件內(nèi)UHPSFRC損傷δ云圖Fig.16 δ contour of specimens with different thickness of outer steel tube

    3.4 內(nèi)、外層鋼管強(qiáng)度的影響

    圖17給出了內(nèi)層鋼管鋼材分別為Q235、Q345、Q390時(shí)UHPSFRCFDST柱在相同爆炸荷載作用下距柱中380 mm處豎向位移時(shí)程曲線圖。由圖可以看出,無(wú)論是峰值位移還是殘余位移,內(nèi)層鋼管強(qiáng)度的影響并不顯著??梢哉f(shuō)明,內(nèi)層鋼管強(qiáng)度對(duì)UHPSFRCFDST柱抗爆性能影響并不明顯。

    圖17 不同強(qiáng)度內(nèi)層鋼管試件柱中位移時(shí)程曲線Fig.17 Displacement time history of specimens with different strength of inner steel tubes

    圖18給出了不同強(qiáng)度內(nèi)層鋼管試件UHPSFRC在相同爆炸荷載作用下的最終破壞形態(tài)。由圖可以看出,各試件UHPSFRC損傷范圍及程度相差不大,說(shuō)明內(nèi)層鋼管強(qiáng)度對(duì)UHPSFRCFDST柱的損傷破壞影響并不顯著。

    圖18 不同強(qiáng)度內(nèi)層鋼管試件UHPSFRC損傷δ云圖Fig.18 δ contour of specimens with different strength of inner steel tube

    圖19給出了外層鋼管鋼材分別為Q235、Q345、Q390時(shí)UHPSFRCFDST柱在相同爆炸荷載作用下距柱中380 mm處豎向位移時(shí)程曲線圖。由圖可以看出,當(dāng)外層鋼管強(qiáng)度分別為Q235和Q345時(shí),試件的峰值位移和殘余位移相差不大,而當(dāng)外層鋼管強(qiáng)度為Q390時(shí),試件的峰值位移和殘余位移顯著減小,說(shuō)明在強(qiáng)度小于Q345范圍內(nèi)提升外層鋼管的強(qiáng)度并不能提升UHPSFRCFDST柱的抗爆性能,而將外層鋼管強(qiáng)度提升到Q390以上卻能夠顯著提升UHPSFRCFDST柱的抗爆性能。

    圖19 不同強(qiáng)度外層鋼管試件柱中位移時(shí)程曲線Fig.19 Displacement time history of specimens with different strength of outer steel tube

    圖20給出了不同強(qiáng)度外層鋼管試件UHPSFRC在相同爆炸荷載作用下的最終破壞形態(tài)。由圖可以看出隨著外層鋼管強(qiáng)度的增大,UHPSFRC損傷范圍沿柱縱向變大,而沿柱橫向減小,雖然差異并不顯著,但也足以說(shuō)明提高外層鋼管強(qiáng)度能夠使得內(nèi)部UHPSFRC裂縫開(kāi)展更為均勻,進(jìn)而提升了試件整體的抗爆耗能性能。

    圖20 不同強(qiáng)度外層鋼管試件UHPSFRC損傷δ云圖Fig.20 δ contour of specimens with different strength of outer steel tube

    4 結(jié) 論

    (1) 本文中試驗(yàn)測(cè)得的爆炸波反射超壓與 ConWep 模擬值及文獻(xiàn)[25]計(jì)算值吻合很好。當(dāng)折合距離為 0.37 m/kg1/3時(shí),爆炸波反射超壓峰值及持時(shí)的實(shí)測(cè)值、ConWep 模擬值及文獻(xiàn)[25]計(jì)算值分別為46.4 MPa和0.89 ms、51.7 MPa和1.0 ms以及49.6 MPa和 0.88 ms。

    (2) 本文建立的3D有限元數(shù)值分析模型能夠有效的分析UHPSFRCFDST 柱在爆炸荷載作用下的動(dòng)態(tài)響應(yīng)。柱中最大位移與殘余位移實(shí)測(cè)值與模擬值誤差在5%以內(nèi)。基于Von Misses應(yīng)力及損傷縮放系數(shù)δ的損傷評(píng)估能夠較好地估計(jì)UHPSFRCFDST柱的損傷破壞。

    (3) 簡(jiǎn)支條件下,當(dāng)空心率Ψ<25%時(shí),減小空心率并不一定能夠提升UHPSFRCFDST柱的抗爆性能;當(dāng)25%<Ψ<56.25%時(shí),減小空心率卻能夠顯著減小UHPSFRCFDST柱的峰值位移和殘余位移,降低其損傷破壞,有效提升其抗爆性能;提高含鋼率亦能夠顯著減小UHPSFRCDST柱的峰值位移及殘余位移,有效降低其損傷破壞,提升其抗爆性能。

    (4) 簡(jiǎn)支條件下,減小內(nèi)、外層鋼管高厚比均能夠顯著減小UHPSFRCFDST柱的峰值位移和殘余位移,同時(shí)使得試件損傷破壞開(kāi)展地更為均勻,有效提升其抗爆性能。

    (5) 簡(jiǎn)支條件下,內(nèi)層鋼管強(qiáng)度對(duì)UHPSFRCFDST柱的峰值位移、殘余位移及其損傷破壞影響并不顯著,故而其對(duì)UHPSFRCDST柱抗爆性能的影響并不明顯;在強(qiáng)度小于Q345范圍內(nèi)提升外層鋼管的強(qiáng)度并不能顯著減小UHPSFRCFDST柱的峰值位移和殘余位移,但卻能夠使得試件損傷破壞開(kāi)展地更為均勻;將外層鋼管強(qiáng)度提升到Q390后,UHPSFRCFDST柱的峰值位移和殘余位移明顯減小,其損傷破壞區(qū)分布更為均勻, UHPSFRCFDST柱的抗爆性能得到了有效提升。

    [ 1 ] 鐘善桐. 鋼管混凝土結(jié)構(gòu)[M]. 北京:清華大學(xué)出版社,2003.

    [ 2 ] 蔡紹懷. 現(xiàn)代鋼管混凝土結(jié)構(gòu)[M]. 北京:人民交通出版社,2003.

    [ 3 ] 韓林海. 鋼管混凝土結(jié)構(gòu)-理論與實(shí)踐[M].北京:科學(xué)出版社,2004.

    [ 4 ] 黃宏.中空夾層鋼管混凝土壓彎構(gòu)件的力學(xué)性能研究[D].福州:福州大學(xué)土木工程學(xué)院, 2006.

    [ 5 ] HABEL K, GAUVREAU P. Response of ultra-high performance fiber reinforced concrete (UHPFRC) impact and static loading [J]. Cement and Concrete Composites, 2008, 30(10):938-946.

    [ 6 ] 巫緒濤.鋼纖維高強(qiáng)混凝土動(dòng)態(tài)力學(xué)性質(zhì)的研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué)力學(xué)和機(jī)械工程系,2006.

    [ 7 ] WU C,OEHLERS D J, REBENTROST M, et al. Blast testing of ultra-high performance fibre and FRP-retrofitted concrete slabs [J]. Engineering Structures, 2009, 31(9): 2060-2069.

    [ 8 ] 孫珊珊.爆炸荷載下鋼管混凝土柱抗爆性能研究[D].西安:長(zhǎng)安大學(xué)土木工程學(xué)院,2013.

    [ 9 ] 崔瑩.爆炸荷載下復(fù)式空心鋼管混凝土柱的動(dòng)態(tài)響應(yīng)及損傷評(píng)估[D].西安:長(zhǎng)安大學(xué)土木工程學(xué)院,2013.

    [10] 李國(guó)強(qiáng),瞿海雁,楊濤春,等.鋼管混凝土柱抗爆性能試驗(yàn)研究[J].建筑結(jié)構(gòu)學(xué)報(bào),2013,13(12):69-76. LI Guoqiang, QU Haiyan, YANG Taochun, et al. Experimental study of concrete-filled steel tubular columns under blast loading [J]. Journal of Building Structures,2013, 13(12):69-76.

    [11] ZHANG Fangrui, WU Chengqing, WANG Hongwei, et al. Numerical simulation of concrete filled steel tube columns against BLAST loads [J]. Thin-Walled Structures, 2015, 92:82-92.

    [12] JAMA H H, BAMBACH M R, NURICK G N, et al. Numerical modeling of square tubular steel beams subjected to transverse blast loads [J]. Thin-Walled Structures, 2009, 47(12):1523-1534.

    [13] LS-DYNA Keyword User’s Manual[M]. Livermore, California: Livermore Software Technology Corporation, 2015.

    [14] LS-DYNA Theoretical Manual[M]. Livermore, California: Livermore Software Technology Corporation, 2015.

    [15] 李忠獻(xiàn),師燕超,史祥生.爆炸荷載作用下鋼筋混凝土板破壞評(píng)定方法[J].建筑結(jié)構(gòu)學(xué)報(bào),2009(6):60-66. LI Zhongxian, SHI Yanchao, SHI Xiangsheng. Damage analysis and assessment of RC slabs under blast load [J]. Journal of Building Structures,2009(6):60-66.

    [16] 師燕超,李忠獻(xiàn).爆炸荷載作用下鋼筋混凝土柱的動(dòng)力響應(yīng)與破壞模式[J].建筑結(jié)構(gòu)學(xué)報(bào),2008(4):112-117. SHI Yanchao, LI Zhongxian. Dynamic response and failure modes of RC columns under blast loading [J]. Journal of Building Structures,2008(4):112-117.

    [17] BAUSCHINGER J.On the change of position of the elastic limit of iron and steel under cyclic variations of stress [J]. Mitt Mech Tech Lab Mü nchen, 1886, 13:1-115.

    [18] MALVAR L J, CRAWFORD J E, WESEVICH J W, et al. A plasticity concrete material model for DYNA3D [J]. International Journal of Impact Engineering, 1997, 19(9): 847-873.

    [19] TU Z, LU Y. Evaluation of typical concrete material models used in hydrocodes for high dynamic response simulations [J]. International Journal of Impact Engineering, 2009, 36(1): 132-146.

    [20] MAO L, BARNETT S, BEGG D, et al. Numerical simulation of ultra high performance fibre reinforced concrete panel subjected to blast loading[J]. International Journal of Impact Engineering, 2014, 64: 91-100.

    [21] MALVAR L J, CRAWFORD J E, MORILL K B,et al. K&C concrete material model release III-automated generation of material model input [J]. K&C Technical Report TR-99-24-B1, 2009.

    [22] NGO T,MENDIS P, KRAUTHAMMER T. Behavior of ultrahigh-strength prestressed concrete panels subjected to blast loading[J]. Journal of Structural Engineering, 2007, 133(11): 1582-1590.

    [23] ISAACS J,MAGALLANES J, REBENTROST M, et al. Exploratory dynamic material characterization tests on ultra-high performance fibre reinforced concrete[C]//Proceedings of 8th International Conference on Shock and Impact Loads on Structures, Adelaide, Australia, 2009.

    [24] Comite Euro-International Du Beton. CEB-FIP Model Code 1990: Design Code[S]. Thomas Telford, London ,1993.

    [25] UFC-3-340-02. Structures to resist the effect of accidental explosions[S].US Department of the Army,Navy and Air Force Technical Manual, 2014.

    [26] WANG F W Y, CHONG Q Y K, LIM C H. Reinforced concrete slab subjected to close-in explosion[C]//The Second International Workshop on Performanceprocection and Strengthening of Structure Under Extreme loading, 2009.

    Numerical simulation and test validation for ultra-high performance steel fiber reinforced concrete-filled double skin steel tube column under blast loading

    XU Shenchun1, LIU Zhongxian2,3, WU Chengqing2,3

    (1. School of Civil Engineering, Tianjin University, Tianjing 300372, China;2. Tianjin Key Laboratory of Civil Structure Protection and Reinforcement Tianjin Pro Chengjian University, Tianjin 300384, China;3. Tianjing Chengjian Univerisity-Adelaide University Disaster Prevention and Mitigation Joint Research Center, Tianjing 300384, China)

    A 3D FE model was developed to analyze dynamic response and damage mechanism of a UHPSFRCFDST(Ultra-High Performance Steel Fiber Reinforced Concrete Filled Double Skin Steel Tube)column under blast loading, they were validated through comparison of simulated results and blast testing ones of the UHPSFRCFDST column. The effects of main variables including hollow ratios, steel ratios, and thicknesses and strengths of inner steel tube and outer steel one on the blast resistant performance of the UHPSFRCFDST column were investigated by utilizing the 3D FE model. The results indicated that the UHPSFRCFDST column has excellent anti-blast performances, and the proposed 3D FEM model can be used to analyze the dynamic response of the UHPSFRCFDST column under blast loading efficiently; reducing hollow ratio and increasing strength of outer steel tube in certain ranges can enhance effectively the blast resistant performance of the UHPSFRCFDST column; increasing steel ratio and reducing height to thickness ratios of inner and outer steel tubes can significantly improve the blast resistant performance of the UHPSFRCFDST column; the strength of inner steel tube has a little effect on the blast resistant performance of the UHPSFRCFDST column, but increasing the strength of outer steel tube in a certain range can improve the blast resistant performance of the UHPSFRCFDST column significantly.

    ultra-high performance steel fiber reinforced concrete-filled double skin steel tube colums;blast loading;dynamic response; failure model;finite element model

    國(guó)家自然科學(xué)基金項(xiàng)目面上項(xiàng)目(51278326);國(guó)家科技支撐計(jì)劃重點(diǎn)項(xiàng)目(2012BAJ07B05);天津市科技支撐計(jì)劃重點(diǎn)項(xiàng)目(14ZCZDSF00016)

    2015-09-23 修改稿收到日期:2015-12-13

    徐慎春 男,博士生,1986年生

    劉中憲 男,副教授,1982年生 E-mail:zhongxian1212@163.com

    TU375.3

    A

    10.13465/j.cnki.jvs.2017.01.007

    猜你喜歡
    抗爆性內(nèi)層外層
    ASTM-CFR辛烷值機(jī)常見(jiàn)問(wèn)題及處理方法
    地聚物超高性能混凝土復(fù)合板抗接觸爆炸試驗(yàn)與數(shù)值模擬
    爆炸與沖擊(2024年6期)2024-10-31 00:00:00
    ◆ 裝飾板材
    ◆ 裝飾板材
    裝飾板材
    ◆ 裝飾板材
    一種溶液探測(cè)傳感器
    傳感器世界(2022年4期)2022-11-24 21:23:50
    航空汽油抗爆性指標(biāo)及測(cè)試方法概述
    95號(hào)汽油切換到92號(hào)汽油對(duì)發(fā)動(dòng)機(jī)性能影響的差異
    一種購(gòu)物袋
    科技資訊(2016年6期)2016-05-14 13:09:55
    麻豆乱淫一区二区| 美女高潮到喷水免费观看| 亚洲精品久久成人aⅴ小说| 巨乳人妻的诱惑在线观看| 美女中出高潮动态图| 成年动漫av网址| 国产午夜精品一二区理论片| 美女脱内裤让男人舔精品视频| 国产亚洲精品久久久久5区| 免费女性裸体啪啪无遮挡网站| 亚洲,欧美精品.| 欧美97在线视频| 午夜免费男女啪啪视频观看| 成年人免费黄色播放视频| 免费女性裸体啪啪无遮挡网站| 丝袜人妻中文字幕| 十分钟在线观看高清视频www| 国产精品欧美亚洲77777| 日日夜夜操网爽| 又大又黄又爽视频免费| 中文精品一卡2卡3卡4更新| 国产精品免费大片| 精品高清国产在线一区| 国产日韩欧美视频二区| 一边摸一边做爽爽视频免费| 午夜影院在线不卡| a级片在线免费高清观看视频| 欧美在线一区亚洲| 女警被强在线播放| av天堂久久9| 精品高清国产在线一区| 亚洲av综合色区一区| 欧美日韩亚洲高清精品| 精品视频人人做人人爽| 亚洲第一青青草原| 欧美老熟妇乱子伦牲交| 精品免费久久久久久久清纯 | 免费女性裸体啪啪无遮挡网站| 国产一卡二卡三卡精品| 欧美黑人精品巨大| 中文字幕人妻丝袜一区二区| 日韩欧美一区视频在线观看| 日韩,欧美,国产一区二区三区| 久久精品人人爽人人爽视色| 亚洲成色77777| 国产av精品麻豆| 大陆偷拍与自拍| 国产免费现黄频在线看| 黄色一级大片看看| 国产伦理片在线播放av一区| 久久精品aⅴ一区二区三区四区| 韩国精品一区二区三区| 亚洲精品第二区| 国产成人a∨麻豆精品| 丝袜人妻中文字幕| 久久久亚洲精品成人影院| 亚洲七黄色美女视频| 国产熟女欧美一区二区| 老汉色av国产亚洲站长工具| 久久精品熟女亚洲av麻豆精品| 狠狠精品人妻久久久久久综合| 2021少妇久久久久久久久久久| 日韩av免费高清视频| 国产97色在线日韩免费| 亚洲一卡2卡3卡4卡5卡精品中文| 中文欧美无线码| 麻豆av在线久日| 亚洲欧美精品自产自拍| 肉色欧美久久久久久久蜜桃| 成年动漫av网址| 男女午夜视频在线观看| 精品少妇黑人巨大在线播放| 日韩av在线免费看完整版不卡| 99香蕉大伊视频| 日韩大片免费观看网站| 久久久久精品国产欧美久久久 | 丝袜人妻中文字幕| 黄频高清免费视频| 国产老妇伦熟女老妇高清| 久久ye,这里只有精品| 黄色片一级片一级黄色片| 日韩一本色道免费dvd| 成年美女黄网站色视频大全免费| 在线观看免费视频网站a站| 亚洲少妇的诱惑av| 晚上一个人看的免费电影| 麻豆乱淫一区二区| 国产日韩欧美视频二区| 麻豆乱淫一区二区| 亚洲欧美日韩另类电影网站| 一级a爱视频在线免费观看| 在线 av 中文字幕| 高清欧美精品videossex| 天天躁日日躁夜夜躁夜夜| 99久久综合免费| 男人爽女人下面视频在线观看| 久久精品国产综合久久久| 老司机靠b影院| 国产精品久久久人人做人人爽| 天天添夜夜摸| 菩萨蛮人人尽说江南好唐韦庄| 女警被强在线播放| 国产精品一区二区免费欧美 | 大码成人一级视频| 国产有黄有色有爽视频| 午夜av观看不卡| 一边亲一边摸免费视频| 亚洲成人手机| 狠狠精品人妻久久久久久综合| 美女午夜性视频免费| 人妻一区二区av| 国产有黄有色有爽视频| 十分钟在线观看高清视频www| 一级毛片女人18水好多 | 久久人妻熟女aⅴ| 女警被强在线播放| 国产免费视频播放在线视频| 水蜜桃什么品种好| 欧美精品高潮呻吟av久久| 叶爱在线成人免费视频播放| 制服诱惑二区| 亚洲情色 制服丝袜| 亚洲午夜精品一区,二区,三区| 亚洲少妇的诱惑av| 成人亚洲精品一区在线观看| 亚洲美女黄色视频免费看| 美女扒开内裤让男人捅视频| 欧美人与性动交α欧美软件| 男女午夜视频在线观看| 女人久久www免费人成看片| 精品国产一区二区三区久久久樱花| 中国美女看黄片| av片东京热男人的天堂| 美国免费a级毛片| 欧美 日韩 精品 国产| 国产精品国产av在线观看| 大话2 男鬼变身卡| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人手机| 999久久久国产精品视频| 在线亚洲精品国产二区图片欧美| 老汉色av国产亚洲站长工具| xxx大片免费视频| 亚洲av片天天在线观看| 亚洲 国产 在线| 免费在线观看完整版高清| 精品福利永久在线观看| 亚洲综合色网址| 午夜福利视频精品| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看| 伊人亚洲综合成人网| 国产精品亚洲av一区麻豆| 色播在线永久视频| 丁香六月天网| 久久中文字幕一级| 中文精品一卡2卡3卡4更新| 午夜91福利影院| 超碰成人久久| 女人高潮潮喷娇喘18禁视频| 国产在线免费精品| 亚洲国产av新网站| 大片免费播放器 马上看| 三上悠亚av全集在线观看| 欧美黑人欧美精品刺激| 中文字幕人妻熟女乱码| 中文字幕高清在线视频| 午夜精品国产一区二区电影| 久久国产精品男人的天堂亚洲| 在线天堂中文资源库| 中文欧美无线码| 天堂8中文在线网| 亚洲欧美日韩另类电影网站| 欧美精品高潮呻吟av久久| 香蕉丝袜av| 少妇人妻 视频| 国产99久久九九免费精品| 男人爽女人下面视频在线观看| 一边摸一边做爽爽视频免费| 女性生殖器流出的白浆| 桃花免费在线播放| 99九九在线精品视频| 好男人电影高清在线观看| 国产欧美亚洲国产| 亚洲国产av影院在线观看| 黄色怎么调成土黄色| 一级a爱视频在线免费观看| 人妻 亚洲 视频| 交换朋友夫妻互换小说| 又粗又硬又长又爽又黄的视频| 日本欧美国产在线视频| 亚洲精品久久成人aⅴ小说| 久久精品国产亚洲av高清一级| 久久精品久久久久久噜噜老黄| 精品一区二区三区四区五区乱码 | 国产国语露脸激情在线看| 国产黄频视频在线观看| 99香蕉大伊视频| 国产精品 欧美亚洲| 亚洲欧洲精品一区二区精品久久久| 一级毛片电影观看| 久久人人爽av亚洲精品天堂| 中文字幕高清在线视频| 久久精品亚洲熟妇少妇任你| 侵犯人妻中文字幕一二三四区| 999精品在线视频| 亚洲 国产 在线| 日韩 欧美 亚洲 中文字幕| 亚洲视频免费观看视频| 免费看不卡的av| 亚洲第一青青草原| 久久久国产精品麻豆| 亚洲中文av在线| 婷婷色综合大香蕉| 99国产精品99久久久久| 狂野欧美激情性bbbbbb| 日本五十路高清| 老司机影院毛片| 亚洲免费av在线视频| 人妻人人澡人人爽人人| 免费少妇av软件| 精品国产乱码久久久久久男人| 中文字幕人妻熟女乱码| 精品少妇内射三级| 亚洲国产毛片av蜜桃av| 男的添女的下面高潮视频| www.熟女人妻精品国产| 黄网站色视频无遮挡免费观看| 在线精品无人区一区二区三| 观看av在线不卡| 日日摸夜夜添夜夜爱| 国产爽快片一区二区三区| 久久鲁丝午夜福利片| 国产一区二区三区av在线| 久久久久久久精品精品| 韩国精品一区二区三区| 免费高清在线观看日韩| 亚洲成国产人片在线观看| 大片免费播放器 马上看| 超碰97精品在线观看| 日韩欧美一区视频在线观看| 十八禁人妻一区二区| 欧美另类一区| 亚洲精品国产区一区二| av有码第一页| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 嫩草影视91久久| 国产亚洲精品久久久久5区| 国产成人91sexporn| 亚洲欧美清纯卡通| 男人添女人高潮全过程视频| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久| 大片免费播放器 马上看| 欧美精品啪啪一区二区三区 | 一级毛片 在线播放| 亚洲色图综合在线观看| 亚洲专区国产一区二区| 成年动漫av网址| 麻豆av在线久日| 我的亚洲天堂| 中文字幕制服av| 欧美日本中文国产一区发布| 又大又黄又爽视频免费| 国产在视频线精品| 日韩电影二区| 男女边摸边吃奶| 亚洲国产av影院在线观看| 一本综合久久免费| 国产在线观看jvid| cao死你这个sao货| 丰满迷人的少妇在线观看| 在线亚洲精品国产二区图片欧美| 成人亚洲精品一区在线观看| 中文字幕高清在线视频| 久久人人爽人人片av| 亚洲av日韩在线播放| 丁香六月天网| 又黄又粗又硬又大视频| 国产日韩欧美视频二区| 国产亚洲精品第一综合不卡| 91麻豆精品激情在线观看国产 | 一区福利在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日韩熟女老妇一区二区性免费视频| 国产成人精品无人区| 国产成人精品久久二区二区91| 国产精品国产三级专区第一集| 国产av精品麻豆| 国产成人精品无人区| 观看av在线不卡| 观看av在线不卡| 中国国产av一级| 超碰97精品在线观看| 女人久久www免费人成看片| 十分钟在线观看高清视频www| 久久久久网色| 一区福利在线观看| 亚洲人成电影观看| 黄色a级毛片大全视频| 久久人人97超碰香蕉20202| 99国产综合亚洲精品| 亚洲成国产人片在线观看| 亚洲一区中文字幕在线| 国产亚洲欧美在线一区二区| 黄片播放在线免费| 熟女av电影| av网站在线播放免费| 极品少妇高潮喷水抽搐| 日韩制服骚丝袜av| 黄色 视频免费看| 亚洲av国产av综合av卡| 国产精品99久久99久久久不卡| av欧美777| 国产精品国产三级专区第一集| 国精品久久久久久国模美| 97在线人人人人妻| 久久久亚洲精品成人影院| 亚洲欧美日韩高清在线视频 | 大型av网站在线播放| 美女扒开内裤让男人捅视频| 欧美亚洲日本最大视频资源| 午夜福利在线免费观看网站| 国产熟女午夜一区二区三区| 欧美国产精品va在线观看不卡| avwww免费| 久久99精品国语久久久| 国产免费一区二区三区四区乱码| 在线观看免费视频网站a站| 成年人黄色毛片网站| 国产免费一区二区三区四区乱码| 国产又色又爽无遮挡免| 国产日韩欧美视频二区| 午夜福利乱码中文字幕| 精品少妇黑人巨大在线播放| 在线看a的网站| 久久精品久久久久久噜噜老黄| 夫妻午夜视频| 国产成人精品久久二区二区免费| 婷婷色av中文字幕| 国产免费又黄又爽又色| 五月开心婷婷网| 亚洲一区二区三区欧美精品| 一本色道久久久久久精品综合| 夫妻午夜视频| 九色亚洲精品在线播放| 国产在线视频一区二区| 亚洲精品国产一区二区精华液| 一级黄色大片毛片| 久久国产精品大桥未久av| 亚洲成国产人片在线观看| 丝袜美足系列| 一个人免费看片子| 久久人妻福利社区极品人妻图片 | 亚洲国产精品成人久久小说| 日本黄色日本黄色录像| 成年女人毛片免费观看观看9 | 国产精品成人在线| 亚洲av电影在线观看一区二区三区| 欧美日韩亚洲高清精品| 男女下面插进去视频免费观看| 青青草视频在线视频观看| 人人妻,人人澡人人爽秒播 | 国产极品粉嫩免费观看在线| 男人爽女人下面视频在线观看| 一二三四社区在线视频社区8| 少妇的丰满在线观看| 麻豆乱淫一区二区| 午夜福利在线免费观看网站| 亚洲欧洲精品一区二区精品久久久| 婷婷成人精品国产| 一级毛片黄色毛片免费观看视频| 极品人妻少妇av视频| 国产福利在线免费观看视频| 国产福利在线免费观看视频| 免费不卡黄色视频| 国产成人a∨麻豆精品| 亚洲少妇的诱惑av| 肉色欧美久久久久久久蜜桃| 可以免费在线观看a视频的电影网站| 我要看黄色一级片免费的| 亚洲专区中文字幕在线| 18禁国产床啪视频网站| 久久精品国产亚洲av高清一级| 亚洲伊人久久精品综合| 日日摸夜夜添夜夜爱| 国产欧美日韩一区二区三区在线| 久久久久精品人妻al黑| 在线观看免费日韩欧美大片| 久久影院123| 久热爱精品视频在线9| 国产欧美日韩一区二区三 | 大香蕉久久网| 亚洲一卡2卡3卡4卡5卡精品中文| 国产高清国产精品国产三级| 91成人精品电影| 国产色视频综合| 天天躁日日躁夜夜躁夜夜| 国产又色又爽无遮挡免| 十分钟在线观看高清视频www| 婷婷色av中文字幕| 亚洲一区二区三区欧美精品| 欧美日韩国产mv在线观看视频| 欧美激情 高清一区二区三区| 少妇人妻久久综合中文| 建设人人有责人人尽责人人享有的| 久久99精品国语久久久| 欧美日韩亚洲综合一区二区三区_| 亚洲av电影在线观看一区二区三区| 久热这里只有精品99| 色综合欧美亚洲国产小说| 欧美97在线视频| 久久久精品94久久精品| 亚洲国产精品999| 天堂8中文在线网| 精品人妻熟女毛片av久久网站| 看免费av毛片| 90打野战视频偷拍视频| 在现免费观看毛片| 天天躁夜夜躁狠狠久久av| 黄网站色视频无遮挡免费观看| 后天国语完整版免费观看| 18禁观看日本| 狂野欧美激情性xxxx| h视频一区二区三区| 黑人猛操日本美女一级片| 国产欧美日韩一区二区三 | 丁香六月欧美| 精品国产乱码久久久久久男人| 亚洲一区中文字幕在线| 国产精品成人在线| av线在线观看网站| 大码成人一级视频| 熟女av电影| 国精品久久久久久国模美| 亚洲人成77777在线视频| 欧美性长视频在线观看| 汤姆久久久久久久影院中文字幕| 999精品在线视频| 精品国产一区二区久久| 美女中出高潮动态图| 久久久精品94久久精品| 啦啦啦 在线观看视频| 大片免费播放器 马上看| 精品免费久久久久久久清纯 | 别揉我奶头~嗯~啊~动态视频 | 国产伦人伦偷精品视频| 一级黄色大片毛片| 99国产精品一区二区蜜桃av | 亚洲中文av在线| 一本综合久久免费| 视频在线观看一区二区三区| 久久久久久人人人人人| 欧美黄色淫秽网站| 精品人妻在线不人妻| 久久国产精品人妻蜜桃| 国产精品.久久久| 欧美人与性动交α欧美软件| 又大又爽又粗| 成人国产一区最新在线观看 | 一本一本久久a久久精品综合妖精| 日本av手机在线免费观看| 精品久久久久久久毛片微露脸 | 亚洲精品国产区一区二| 欧美日韩成人在线一区二区| 国产视频一区二区在线看| 电影成人av| 老司机影院成人| 国产免费福利视频在线观看| 咕卡用的链子| 秋霞在线观看毛片| 亚洲国产看品久久| 久久精品国产亚洲av高清一级| 亚洲天堂av无毛| 美女视频免费永久观看网站| 免费在线观看黄色视频的| 蜜桃国产av成人99| 大香蕉久久成人网| 国产视频一区二区在线看| 国产成人精品无人区| 一本久久精品| 久久亚洲国产成人精品v| 2018国产大陆天天弄谢| 在线观看人妻少妇| 少妇人妻 视频| 婷婷色av中文字幕| 亚洲九九香蕉| 亚洲综合色网址| 欧美在线黄色| 免费观看a级毛片全部| 真人做人爱边吃奶动态| 国产主播在线观看一区二区 | 国产精品一区二区精品视频观看| 最新的欧美精品一区二区| 一区二区三区激情视频| 又大又爽又粗| 老司机在亚洲福利影院| 黄频高清免费视频| 中文乱码字字幕精品一区二区三区| 你懂的网址亚洲精品在线观看| 国产在视频线精品| 90打野战视频偷拍视频| 国产成人一区二区三区免费视频网站 | 国产精品久久久人人做人人爽| 一本色道久久久久久精品综合| 久久久久久人人人人人| 久久精品亚洲av国产电影网| 亚洲av综合色区一区| 美女主播在线视频| a级毛片黄视频| 国产有黄有色有爽视频| 国产成人a∨麻豆精品| 啦啦啦 在线观看视频| 啦啦啦在线免费观看视频4| 国产又色又爽无遮挡免| 久久热在线av| 亚洲国产毛片av蜜桃av| 少妇的丰满在线观看| 亚洲第一av免费看| 国产精品免费大片| 人妻一区二区av| 高清av免费在线| 巨乳人妻的诱惑在线观看| 又粗又硬又长又爽又黄的视频| 欧美人与性动交α欧美软件| 香蕉丝袜av| 午夜福利视频精品| 视频区欧美日本亚洲| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成77777在线视频| 亚洲av片天天在线观看| 国产精品.久久久| 精品少妇内射三级| 肉色欧美久久久久久久蜜桃| 黄色视频在线播放观看不卡| 最新的欧美精品一区二区| 日韩中文字幕欧美一区二区 | 美女国产高潮福利片在线看| 亚洲av日韩在线播放| 侵犯人妻中文字幕一二三四区| 国产成人精品久久久久久| 又紧又爽又黄一区二区| 国产一区二区三区av在线| 飞空精品影院首页| 日韩伦理黄色片| 老熟女久久久| 成人国语在线视频| 亚洲欧美色中文字幕在线| 看免费成人av毛片| 国产亚洲欧美精品永久| 国产精品麻豆人妻色哟哟久久| 自拍欧美九色日韩亚洲蝌蚪91| videos熟女内射| 狂野欧美激情性xxxx| 亚洲成色77777| 精品国产一区二区三区久久久樱花| 男女免费视频国产| 黄色毛片三级朝国网站| 亚洲视频免费观看视频| 久久精品久久久久久噜噜老黄| 热re99久久精品国产66热6| 欧美日本中文国产一区发布| 高清av免费在线| 丝袜喷水一区| 国产97色在线日韩免费| 国产无遮挡羞羞视频在线观看| 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 波野结衣二区三区在线| 久久久久精品国产欧美久久久 | 久久精品国产a三级三级三级| 婷婷色综合大香蕉| 婷婷色av中文字幕| 91国产中文字幕| 丝袜脚勾引网站| 成年美女黄网站色视频大全免费| 精品卡一卡二卡四卡免费| 成人免费观看视频高清| xxxhd国产人妻xxx| 亚洲av美国av| 男女高潮啪啪啪动态图| 少妇被粗大的猛进出69影院| 一边亲一边摸免费视频| 搡老岳熟女国产| 一本—道久久a久久精品蜜桃钙片| 女人高潮潮喷娇喘18禁视频| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 国产一区二区三区av在线| 热re99久久精品国产66热6| 中文字幕人妻熟女乱码| 97在线人人人人妻| 午夜激情久久久久久久| 两个人免费观看高清视频| 女人精品久久久久毛片| 天天影视国产精品| 亚洲精品乱久久久久久| 91成人精品电影| 国产97色在线日韩免费| 国产野战对白在线观看| 久久久精品国产亚洲av高清涩受| 国产深夜福利视频在线观看| 一区二区av电影网| 国产深夜福利视频在线观看| 桃花免费在线播放| av在线app专区| 桃花免费在线播放| 国产免费又黄又爽又色| 午夜激情久久久久久久| avwww免费| 欧美少妇被猛烈插入视频| 欧美激情 高清一区二区三区| 天堂8中文在线网| 国产欧美亚洲国产|