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

    火災(zāi)條件下冷彎薄壁型鋼樓板體系的耐火性能

    2021-09-10 01:35楊成羅浪宋謙益
    土木建筑與環(huán)境工程 2021年5期

    楊成 羅浪 宋謙益

    摘 要:冷彎薄壁型鋼樓板的耐火性能是決定其作為裝配式結(jié)構(gòu)能否推廣的重要因素,目前可供借鑒的足尺火災(zāi)試驗(yàn)數(shù)量有限,且大部分?jǐn)?shù)值模擬局限于用二維傳熱模型開展溫度場分析,或基于簡化的線性溫度分布開展三維模型的熱力耦合分析。冷彎薄壁型鋼樓板構(gòu)造復(fù)雜,利用二維傳熱模型以及簡化線性溫度分布假定開展的熱力耦合分析,由于未考慮三維結(jié)構(gòu)在受火時(shí)真實(shí)的非線性傳熱機(jī)制,其所反映的結(jié)構(gòu)抗火性能與真實(shí)情況可能存在顯著差異。為更準(zhǔn)確地描述這類結(jié)構(gòu)的耐火性能,針對由冷彎薄壁型鋼骨架、結(jié)構(gòu)膠合板和石膏板、巖棉等材料組成的樓板體系建立三維傳熱模型,按照ISO-834標(biāo)準(zhǔn)升溫環(huán)境模擬火災(zāi)情境,進(jìn)行非線性傳熱過程和熱力耦合分析,描述三維非線性傳熱機(jī)制和具有熱力耦合特征的結(jié)構(gòu)行為;提出更為準(zhǔn)確的接觸行為建模方法,以反映自攻螺釘、板材、輕鋼構(gòu)件之間接觸傳力關(guān)系和變形協(xié)調(diào)特征。在與兩組足尺樓板體系火災(zāi)試驗(yàn)數(shù)據(jù)對比的基礎(chǔ)上,給出由數(shù)值模擬得到的結(jié)構(gòu)行為機(jī)理解釋。研究表明:提出的三維傳熱和三維熱力耦合分析數(shù)值模型能準(zhǔn)確反映火災(zāi)條件下結(jié)構(gòu)的力學(xué)行為,模擬的結(jié)構(gòu)非線性溫度分布特征、托梁的變形演化機(jī)制以及最終破壞形態(tài)與試驗(yàn)結(jié)果吻合較好。

    關(guān)鍵詞:冷彎薄壁型鋼;樓板體系;耐火性能;火災(zāi)情境;傳熱過程;力學(xué)行為

    Abstract: Cold-formed light gauge steel framing (LSF) structure is a new type of fabricated steel structure. Fire resisting performance of LSF structure is the key factor that determines whether it can be promoted as a fabricated structure. However, for the LSF floor as one of the main load-bearing systems in LSF structures, the number of full-scale fire tests available for reference is limited. Most numerical simulations of LSF floors are limited to using two-dimensional heat transfer models to conduct temperature field analysis, or to conduct thermal-mechanical coupling analysis of three-dimensional models based on simplified linear temperature distributions. The difference between the analysis results of two-dimension and three-dimension structural members could accumulate and evolve into the essential discrepancy during force transfer when the structural scale is large and the system is complex. In order to describe the fire resistance of this type of structure more accurately, a three-dimensional heat transfer model has been established for the floor system composed of cold-formed thin-walled steel framework, structural plywood, gypsum board, rock wool and other materials.According to the ISO-834 standard,the heating environment is used to simulate the fire situation, the nonlinear heat transfer process and the thermal-mechanical coupling analysis have been performed, and the three-dimensional nonlinear heat transfer mechanism and the structural behavior with thermal-mechanical coupling characteristics have been described. In addition, a more accurate modeling method of contact behavior has been proposed to reflect the contact force relationship and deformation coordination characteristics among self-tapping screws, structural plates, and light steel members.Based on the comparison with two sets of fire testing data for the full-scale floor systems, the explanation of the structure behavior mechanism obtained by numerical simulation has been given. In general, this model based on three-dimensional heat conduction and thermal-mechanical coupling can effectively describe the LSF spatial mechanical behavior under fire conditions. The results show that the proposed three-dimensional heat transfer and three-dimensional thermal-mechanical coupling analysis numerical model can be utilized to predict the fire resistance of LSF floor systems and evaluate the effective factors on structural high-temperature behaviors via parametric analysis.The simulated non-linear temperature distribution characteristics of the structure, the deformation evolution mechanism of the joists and the final failure form are in good agreement with the test results.

    Keywords: cold-formed light gauge steel; floor systems; fire resistance; fire conditions; heat transfer process; mechanical behavior

    冷彎薄壁型鋼結(jié)構(gòu)作為一種裝配式鋼結(jié)構(gòu),主要通過自攻螺釘將冷彎薄壁型鋼骨架連接成框架,并在框架一側(cè)或兩側(cè)覆以結(jié)構(gòu)板材形成墻體、樓板等受力構(gòu)件。并在結(jié)構(gòu)構(gòu)件基體外再覆以石膏板、玻鎂板等具備防火、隔聲、裝飾等性能的建筑板材,以滿足建筑結(jié)構(gòu)的功能性需求。鋼結(jié)構(gòu)在高溫下的強(qiáng)度和剛度衰退非常顯著[1-2],冷彎薄壁型鋼樓板作為該結(jié)構(gòu)體系中重要的承重結(jié)構(gòu)構(gòu)件,其耐火性能得到日益關(guān)注[3-5]。

    其他國家開展冷彎薄壁型鋼結(jié)構(gòu)火災(zāi)試驗(yàn)研究[6-8]已積累了一定的經(jīng)驗(yàn)和數(shù)據(jù),而中國的試驗(yàn)研究[9-11]相對有限。相比火災(zāi)試驗(yàn)耗資大、耗時(shí)長,數(shù)值分析提供了更便捷的方式以了解其在火災(zāi)下的溫度分布[12-15]和結(jié)構(gòu)性能[16]。Keerthan等[17]利用SAFIR進(jìn)行了冷彎薄壁型鋼墻體的二維瞬態(tài)傳熱模擬。Feng等[12]利用ABAQUS對冷彎薄壁型鋼墻體截面建立了二維傳熱模型。Baleshan等[18]和Jatheeshan等[19]使用簡化的線性溫度分布對火災(zāi)下冷彎薄壁型鋼托梁的結(jié)構(gòu)行為進(jìn)行了模擬。

    近年來,冷彎薄壁型鋼樓板的耐火性能數(shù)值模擬主要存在以下局限性:1)基于結(jié)構(gòu)體系截面建立的二維傳熱模型僅能考慮單個(gè)截面的溫度分布和升溫規(guī)律。但三維空間中冷彎薄壁型鋼樓板體系的板材組合情況、托梁腹板開孔、空腔和托梁截面形狀的改變以及結(jié)構(gòu)部分受火時(shí)的情況,均會(huì)導(dǎo)致結(jié)構(gòu)體系溫度發(fā)展不均勻,因此,二維模型適用性非常有限,需采用三維模型進(jìn)行傳熱模擬[20]。同時(shí),二維傳熱模型的計(jì)算結(jié)果無法作為空間溫度場數(shù)據(jù)用于火災(zāi)下結(jié)構(gòu)熱力耦合的研究。2)一些三維熱力耦合分析模型使用的溫度場數(shù)據(jù)采用簡化的線性溫度分布假定[18],但無法考慮樓板體系由于復(fù)雜構(gòu)造等因素導(dǎo)致的非線性溫度分布及其影響。3)既有的三維熱力耦合模型往往是一個(gè)獨(dú)立的鋼托梁構(gòu)件[19],該模型考慮的樓面板材、鋼托梁的約束情況以及荷載施加方式對鋼托梁性能造成的影響過于簡化。因此,有必要提出一種能夠更充分考慮冷彎薄壁型鋼樓板結(jié)構(gòu)體系耦合熱力行為的模擬方法:首先建立三維傳熱模型,計(jì)算并得出結(jié)構(gòu)的非線性溫度分布和升溫規(guī)律,繼而分析熱力耦合效應(yīng)下冷彎薄壁型鋼結(jié)構(gòu)體系的結(jié)構(gòu)性能。

    參考既有試驗(yàn)的參數(shù)信息,筆者基于非線性有限元平臺(tái)ABAQUS建立了冷彎薄壁型鋼樓板的三維傳熱模型;根據(jù)既有研究經(jīng)驗(yàn)和試算對比,選取材料的熱工和力學(xué)性能參數(shù),計(jì)算分析樓板在ISO-834標(biāo)準(zhǔn)火災(zāi)條件[21]下的傳熱行為,并利用兩組火災(zāi)試驗(yàn)結(jié)果[8]對傳熱分析結(jié)果進(jìn)行驗(yàn)證;基于驗(yàn)證的傳熱分析結(jié)果,建立三維熱力耦合模型,對樓板體系在火災(zāi)下的結(jié)構(gòu)行為開展模擬分析,進(jìn)一步驗(yàn)證耐火極限試驗(yàn)結(jié)果[8]。

    1 三維傳熱模型

    1.1 模型的建立

    以Baleshan等[8]的試件為分析對象,利用有限元分析平臺(tái)ABAQUS建立兩個(gè)2 400 mm×2 100 mm(長×寬)的冷彎薄壁型鋼樓板模型,包括4根鋼托梁和兩根導(dǎo)軌組成的鋼骨架以及受火側(cè)、背火側(cè)板材。所用冷彎薄壁型鋼的鋼材等級為G500,采用C型截面鋼(托梁C180 mm×40 mm×15 mm×1.15 mm,腹板高度h×翼緣寬度b×卷邊寬度a×厚度t)及U型截面鋼(導(dǎo)軌U182 mm×50 mm×1.15 mm,h×b×t),鋼托梁(J1、J2、J3、J4)之間間距S均為600 mm。試件沿截面厚度方向的板材組合和鋼托梁布置如圖1所示。

    值得注意的是,試件在受火側(cè)空腔面設(shè)有3個(gè)200 mm寬的石膏板背擋塊,如圖2(a)所示。此處空腔形狀已發(fā)生改變,結(jié)構(gòu)沿跨度方向溫度分布不均勻,故二維截面?zhèn)鳠崮P筒辉龠m用[22]。

    傳熱分析數(shù)值模擬中,石膏板、巖棉和膠合板采用8節(jié)點(diǎn)三維實(shí)體傳熱單元(DC3D8),冷彎薄壁型鋼骨架采用4節(jié)點(diǎn)殼體傳熱單元(DS4)。根據(jù)敏感性分析結(jié)果,選擇模型平面網(wǎng)格密度為80 mm,在鋼骨架與板材相交區(qū)域網(wǎng)格細(xì)化為10 mm;厚度方向網(wǎng)格密度為5 mm,它決定了板材厚度方向的傳熱情況,因此,在此方向采用更細(xì)的網(wǎng)格劃分,如圖2(b)所示。

    采用Baleshan等[8]的火災(zāi)試驗(yàn)結(jié)果,以對比驗(yàn)證冷彎薄壁型鋼樓板的傳熱和熱力耦合分析結(jié)果。

    在試驗(yàn)開始前,以增量的方式給樓板試件背火側(cè)施壓至荷載比0.4,然后受火側(cè)開始升溫,直至結(jié)構(gòu)喪失耐火承載力。試件1的耐火極限時(shí)間為107 min,試件2為139 min。試驗(yàn)后,觀察結(jié)果發(fā)現(xiàn):試件1受火側(cè)石膏板有燒蝕、脫落的現(xiàn)象,其背火側(cè)膠合板在空腔側(cè)形成燒焦層,其他表面狀況良好,基本保持了背火側(cè)板材的完整性。試件2受火側(cè)巖棉發(fā)生收縮,導(dǎo)致接縫的開放和內(nèi)側(cè)石膏板層的暴露、燒蝕及脫落。

    在后處理過程中,提取與Baleshan等[8]試驗(yàn)相對應(yīng)的試件溫度測點(diǎn)(見圖3)的平均溫度值(如:P1-P2間測點(diǎn)的平均值),輸出時(shí)間溫度曲線。為監(jiān)測試件在火災(zāi)下的溫度場分布和各部件的升溫規(guī)律,試件背火面、受火面以及各板材之間各設(shè)置了5個(gè)熱電偶,如圖3所示。其中,各層的其中3個(gè)熱電偶分別位于中間腔體中垂線的0.25L、0.5L、0.75L處,另外2個(gè)熱電偶分別位于兩側(cè)腔體中垂線的0.5L處。內(nèi)側(cè)兩根托梁分別在其冷、熱翼緣和腹板的0.25L、0.5L和0.75L處設(shè)置熱電偶。外側(cè)兩根托梁只在0.5L處的冷、熱翼緣以及腹板設(shè)置熱電偶。

    1.2 邊界條件

    由于實(shí)際樓板體系中鋼骨架與兩側(cè)板材以及各板材之間接觸非常緊密,故對其相交部分定義為綁定連接,以確保它們之間的固體傳熱。

    在分析過程中,受火側(cè)的火災(zāi)升溫曲線是根據(jù)ISO-834標(biāo)準(zhǔn)升溫曲線[21]確定的,同時(shí)考慮環(huán)境溫度為20 ℃。Stefan-Boltzmann常數(shù)取5.668×10-8 W/m2/℃4,既有研究對于熱對流系數(shù)和熱輻射系數(shù)的取值見表1,彼此之間沒有顯著差異。參考Song等[23-24]的研究,受火側(cè)對流系數(shù)與輻射系數(shù)分別取25 W/(m2·℃)和0.9;背火側(cè)分別取10 W/(m2·℃)和0.7。模型空腔內(nèi)均未填充隔熱材料,因此,腔體內(nèi)部表面設(shè)置空腔輻射系數(shù),取值0.7,如圖4所示。

    1.3 材料熱工參數(shù)的確定

    選取高溫下石膏板、巖棉、膠合板與冷彎薄壁型鋼合適的熱工參數(shù)(如:導(dǎo)熱系數(shù)、比熱容和相對密度)是傳熱分析的關(guān)鍵。石膏板的熱工參數(shù)主要受到內(nèi)部化學(xué)成分、燒蝕以及含水率的影響,為此分別采用Keerthan等[17]提出的高溫下石膏板的導(dǎo)熱系數(shù)、相對密度及Feng等[12]提出的比熱容。巖棉比熱容和密度在高溫下變化不大,采用Jatheeshan等[26]的取值,其導(dǎo)熱系數(shù)則采用Sterner等[27]的數(shù)據(jù)。膠合板和鋼材的熱工參數(shù)分別按照J(rèn)atheeshan等[26]和BS EN 1993-1-2[28]的建議取值。

    1.3.1 石膏板

    根據(jù)已有試驗(yàn)結(jié)果[12-14, 17],選取合適的石膏板熱工參數(shù)進(jìn)行傳熱模擬。石膏板在火災(zāi)影響下煅燒收縮形成小裂縫后部分脫落。這將導(dǎo)致更多的熱量通過樓板腔體,故選用參數(shù)[17]考慮了1 000 ℃后石膏板開裂造成導(dǎo)熱系數(shù)增大的影響,如圖5所示。

    圖6顯示了使用的石膏板比熱變化曲線。在100 ℃時(shí),石膏板內(nèi)部結(jié)合水吸熱蒸發(fā),溫度上升將會(huì)有一個(gè)延遲,這會(huì)顯著影響石膏板的比熱。在125 ℃時(shí),石膏板比熱達(dá)到峰值24 572.32 J/(kg·℃)[12]。

    圖7為使用的石膏板的相對密度曲線。125~175 ℃時(shí),石膏板因高溫產(chǎn)生的質(zhì)量損失為10%,在接下來的火災(zāi)過程中,其密度保持不變[17]。

    1.3.2 巖棉

    圖8給出了幾種既有的巖棉導(dǎo)熱系數(shù)曲線[15, 17, 27, 29]。各地巖棉的化學(xué)成分有所不同,熱性能略有差異,但其導(dǎo)熱系數(shù)均隨溫度的升高而上升。在分析中,巖棉比熱容采用定值840 J/(kg·℃),密度恒為100 kg/m3[26]。

    1.3.3 鋼材

    鋼材的導(dǎo)熱系數(shù)和比熱容依據(jù)歐洲規(guī)范BS EN 1993-1-2[28]的建議公式取值,該公式考慮了鋼材在735 ℃時(shí)比熱容出現(xiàn)峰值的現(xiàn)象,如圖9、圖10所示。在整個(gè)火災(zāi)過程中,鋼材的密度保持7 850 kg/m3不變[26]。

    1.3.4 膠合板 在試件1樓板體系中,膠合板被用作背火側(cè)結(jié)構(gòu)板,仍需考慮其在火災(zāi)下的燒蝕行為。根據(jù)Jatheeshan等[26]的建議對高溫下的膠合板導(dǎo)熱系數(shù)進(jìn)行取值,其導(dǎo)熱系數(shù)逐漸增大,如圖11所示。

    由于試件1膠合板未直接接觸靠火側(cè)石膏板,所以忽略100~105 ℃的比熱峰值,取恒定值1 500 J/(kg·℃)[26]。膠合板的含水量約為9%~11%,在100 ℃時(shí)板材內(nèi)部水分以水蒸氣的形式釋放出來。同時(shí),其在高溫下燃燒的現(xiàn)象,均會(huì)導(dǎo)致相對密度隨著溫度的升高而降低。膠合板環(huán)境溫度下密度為500 kg/m3[26],采用BS EN 1995-1-2[30]的建議進(jìn)行高溫下相對密度的取值,如圖12所示。

    1.4 傳熱分析與驗(yàn)證

    對上述說明的三維傳熱模型進(jìn)行升溫計(jì)算,在后處理中提取模型對應(yīng)試驗(yàn)測溫位置的節(jié)點(diǎn)溫度,輸出節(jié)點(diǎn)的平均時(shí)間溫度曲線,然后與Baleshan等[8]的試驗(yàn)結(jié)果進(jìn)行對比分析。

    兩個(gè)試件在數(shù)值模擬和試驗(yàn)中不同板材表面(受火面、背火面、板材間的面、受火側(cè)空腔表面以及背火側(cè)空腔表面)上的平均溫度發(fā)展曲線,如圖13、圖14所示。

    在圖13中,數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果吻合較好,僅P1-P2面的計(jì)算結(jié)果在65 min后略高于試驗(yàn)值,原因是采用的石膏板導(dǎo)熱系數(shù)在1 000 ℃后快速增大,熱量傳遞增加。試件2傳熱分析的計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)整體吻合較好,見圖14。

    兩個(gè)試件在數(shù)值模擬和試驗(yàn)中兩根中間托梁J2、J3(熱翼緣、腹板、冷翼緣)的平均溫度發(fā)展曲線見圖15~圖16。

    在圖15~圖16中,托梁模擬計(jì)算的溫度在60 min前均略低于試驗(yàn)值,原因是板材內(nèi)部水分在高溫下蒸發(fā)成水蒸氣,導(dǎo)致空腔內(nèi)存在水分子的遷移、相變過程,這會(huì)使托梁上的溫度測點(diǎn)在前期所測溫度值偏高。

    圖15中,火災(zāi)發(fā)生60 min后,兩根托梁熱翼緣的升溫與試驗(yàn)結(jié)果非常接近,托梁腹板、冷翼緣的計(jì)算值與試驗(yàn)值吻合較好,最大溫度差值為31 ℃。

    在圖16中,鋼托梁腹板與冷翼緣的計(jì)算升溫曲線與試驗(yàn)數(shù)據(jù)吻合較好。托梁J2和J3熱翼緣的計(jì)算結(jié)果略小于試驗(yàn)所得數(shù)據(jù),可給出的解釋是:模型未考慮托梁中自攻螺釘?shù)臒針蛐?yīng),略微減少了傳遞至熱翼緣的熱量;使用的材料熱工參數(shù)與試驗(yàn)材料有一定的差異。

    整個(gè)火災(zāi)過程中,試件1和試件2背火側(cè)石膏板的最高溫度均未超過160 ℃,試件未喪失隔熱性。

    1.5 與二維傳熱模型的對比分析

    二維傳熱模型一般采用結(jié)構(gòu)中部截面建立模型,僅能考慮單個(gè)截面的溫度發(fā)展,不能考慮到結(jié)構(gòu)沿跨度方向截面的變化及結(jié)構(gòu)部分受火等因素造成的不均勻溫度發(fā)展[22],有使用局限性。

    以試件1中0.25和0.5長度處截面建立二維傳熱模型,對比三維模型分析托梁J2熱翼緣的溫度發(fā)展情況,具體位置及編號見圖17。二維傳熱模型中,相交部分定義綁定連接;材料熱工參數(shù)和熱邊界條件同三維傳熱模型一致。

    圖18顯示了以A和B點(diǎn)截面建立的二維傳熱模型在ISO-834火災(zāi)下107 min時(shí)的溫度分布。與三維傳熱模型結(jié)果對比如圖19(a)所示,可以發(fā)現(xiàn),二維模型A點(diǎn)的升溫遠(yuǎn)低于三維模型A點(diǎn),是因?yàn)楦鶕?jù)A點(diǎn)截面建立的二維傳熱模型默認(rèn)結(jié)構(gòu)沿跨度所有截面形狀均一致,然而受火側(cè)空腔面存在200 mm寬的背擋塊石膏板(圖17),空腔截面形狀已改變。由圖19(b)可見,B點(diǎn)的升溫存在一定差異,是因?yàn)槿S模型中一部分熱量被背擋塊石膏板吸收,B點(diǎn)建立的二維模型不能考慮到背擋塊的存在,因此,三維模型B點(diǎn)處的升溫略小于二維模型B點(diǎn)。同時(shí),圖19中二維模型的A、B截面之間J2熱翼緣的升溫曲線差異較大。

    整體來看,三維傳熱模型的計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)整體吻合程度較好,且能考慮到結(jié)構(gòu)構(gòu)造等因素造成的不均勻溫度發(fā)展。

    2 三維熱力耦合分析

    以傳熱分析得到的不均勻溫度分布為溫度場基礎(chǔ),建立三維熱力耦合模型,對火災(zāi)環(huán)境下冷彎薄壁型鋼樓板體系的結(jié)構(gòu)性能進(jìn)行分析。

    2.1 模型的建立

    由于受火側(cè)石膏板抗拉強(qiáng)度和剛度很低,且火災(zāi)下很快發(fā)生開裂脫落等現(xiàn)象,因此,未考慮其對鋼骨架的約束作用[31-32]。建立的數(shù)值模型僅包含鋼托梁、鋼導(dǎo)軌和背火側(cè)板材。熱力耦合分析中,背火側(cè)板采用8節(jié)點(diǎn)三維實(shí)體單元(C3D8R),冷彎薄壁型鋼骨架采用4節(jié)點(diǎn)殼體單元(S4R)。根據(jù)試驗(yàn)的加荷方式以及托梁翼緣自攻螺釘?shù)奈恢?,對整個(gè)模型施加了適當(dāng)?shù)倪吔鐥l件,如圖20所示。邊界條件包括:結(jié)構(gòu)兩側(cè)鋼導(dǎo)軌的簡支約束,以及通過耦合連接(模擬自攻螺釘作用)約束托梁冷翼緣的位移。在4根托梁的0.25和0.75長度處分別以增量的方式在背火側(cè)板材上施加豎向荷載,逐漸增大至4.5 kN(荷載比0.4)后保持不變,然后導(dǎo)入傳熱分析結(jié)果進(jìn)行熱力耦合分析。

    此外,需要補(bǔ)充說明導(dǎo)軌、托梁、背火側(cè)板材之間的連接關(guān)系。導(dǎo)軌與托梁在相交部分定義為綁定連接。為了更好地模擬自攻螺釘作用,鋼骨架與背火側(cè)板材兩者間定義為“點(diǎn)環(huán)面”的耦合連接(圖20),以更準(zhǔn)確約束托梁冷翼緣x、y、z方向的位移。同時(shí),為體現(xiàn)兩者間接觸傳力關(guān)系,在法線方向設(shè)置硬接觸,切線方向設(shè)置摩擦系數(shù)0.3。

    2.2 材料力學(xué)性能

    所用1.15 mm厚的G500冷彎薄壁型鋼屈服強(qiáng)度和彈性模量分別為612、210 260 MPa,泊松比為0.3。根據(jù)文獻(xiàn)[16],采用彈塑性材料模型模擬火災(zāi)環(huán)境下冷彎薄壁型鋼的非線性行為,并且鋼材的彈性模量、屈服強(qiáng)度等力學(xué)性能隨著溫度的升高而下降。

    表2給出了Kankanamge等[33]提出的冷彎型鋼高溫下力學(xué)性能折減系數(shù),以精確模擬高溫下冷彎薄壁型鋼力學(xué)性能降低的影響。

    背火側(cè)膠合板與石膏板的力學(xué)性能見表3。膠合板根據(jù)BS EN 1995-1-2[30]對高溫下木材的力學(xué)性能進(jìn)行調(diào)整;石膏板由于其常溫強(qiáng)度較小,未予以高溫下力學(xué)性能的調(diào)整。

    為確定建立的樓板模型與試驗(yàn)試件力學(xué)性能符合,計(jì)算整個(gè)樓板體系在20 ℃下的極限承載能力為94 kN,與Baleshan等[8]試驗(yàn)的90 kN極限荷載工況較為相近,可做參考。

    2.3 初始幾何缺陷和殘余應(yīng)力

    初始幾何缺陷會(huì)影響結(jié)構(gòu)性能。根據(jù)文獻(xiàn)[18],模型中設(shè)置了初始幾何缺陷,初始幾何缺陷的大小為0.000 6b(b為腹板高度),然后進(jìn)行熱力耦合計(jì)算。

    殘余應(yīng)力影響托梁抗彎強(qiáng)度,導(dǎo)致其抗彎強(qiáng)度降低。隨著溫度的升高,托梁內(nèi)部殘余應(yīng)力迅速減小。根據(jù)Baleshan等[18]的研究結(jié)論以及既有研究常用方法,火災(zāi)條件下可以不考慮殘余應(yīng)力。

    2.4 耐火失效準(zhǔn)則

    根據(jù)《建筑構(gòu)件耐火試驗(yàn)第1部分:通用要求》(GB/T 9978.1—2008)[34]中相關(guān)規(guī)定,冷彎薄壁型鋼承重樓板體系的耐火極限模擬分析以穩(wěn)定性和隔熱性為判別標(biāo)準(zhǔn)。

    1)穩(wěn)定性。對于抗彎構(gòu)件,試件超過以下任一判定準(zhǔn)則時(shí),均認(rèn)為試件喪失承載能力。

    2)隔熱性。試件背火面溫度溫升發(fā)生超過以下任一限定情況均認(rèn)為試件喪失隔熱性:平均溫度溫升超過初始平均溫度140 ℃;任一點(diǎn)的溫度溫升超過初始溫度180 ℃。

    3 三維熱力耦合分析的驗(yàn)證

    火災(zāi)條件下,鋼托梁強(qiáng)度下降,且非均勻升溫膨脹,試件撓度逐漸增大。整個(gè)試驗(yàn)過程中,試件背火側(cè)溫度均低于耐火失效溫度160 ℃(圖13、圖14),滿足隔熱性要求?;馂?zāi)升溫環(huán)境下,托梁撓度迅速增大,從而喪失穩(wěn)定性。

    火災(zāi)試驗(yàn)和數(shù)值分析中,托梁J2的冷、熱翼緣和腹板溫度均略高于托梁J3,因此,熱力耦合分析驗(yàn)證主要以火災(zāi)下承載力最低的托梁J2作為研究對象。

    3.1 試件1的模擬分析

    3.1.1 托梁的撓度發(fā)展

    圖21(a)為環(huán)境溫度下加荷步驟后的托梁J2,其最大垂直撓度為4.56 mm。圖21(b)顯示了火災(zāi)下升溫114 min時(shí),由于鋼材受熱膨脹和高溫強(qiáng)度降低的影響,托梁J2的跨中撓度顯著增大,達(dá)到32.48 mm。

    火災(zāi)環(huán)境下,托梁J2跨中撓度在110 min時(shí)有顯著增大,構(gòu)件逐漸喪失穩(wěn)定性,如圖22所示。最終得到的耐火失效時(shí)間為114 min,文獻(xiàn)[18]數(shù)值模擬預(yù)測的失效時(shí)間為120 min。從撓度分析結(jié)果看,相比文獻(xiàn)[18],模型的精度改善作用明顯,原因?yàn)椋菏褂玫臏囟葓鰹槿S傳熱模型計(jì)算的非線性結(jié)果,而文獻(xiàn)[18]施加的溫度場為沿跨度方向均相同的簡化線性溫度分布;對背火側(cè)板材、輕鋼構(gòu)件之間的連接方式和接觸傳力關(guān)系定義更為精細(xì)。

    數(shù)值模擬的托梁冷翼緣處的膠合板約束始終存在,直到結(jié)構(gòu)耐火失效。而事實(shí)上,膠合板部分燒焦以及冷翼緣自攻螺釘屈曲等將對托梁撓度的發(fā)展造成影響。因此,數(shù)值模擬的耐火時(shí)間比試驗(yàn)結(jié)果略長。

    3.1.2 托梁的應(yīng)力發(fā)展

    起初托梁受力產(chǎn)生應(yīng)力,隨著火災(zāi)下托梁內(nèi)部溫度升高,強(qiáng)度下降,最后達(dá)到與溫度相關(guān)的屈服應(yīng)力,托梁耐火失效時(shí)的von Mises應(yīng)力分布如圖23所示。在圖23中,托梁J2支座和加荷處的冷翼緣出現(xiàn)最大von Mises應(yīng)力。圖24顯示了這兩處von Mises應(yīng)力與屈服應(yīng)力隨時(shí)間變化的曲線。

    從圖24可看出,支座處冷翼緣與加荷處冷翼緣達(dá)到屈服強(qiáng)度的時(shí)間有所差異,在114 min時(shí)托梁喪失穩(wěn)定性。圖24中的虛線有所差異,是因?yàn)槭褂玫臏囟葓鰹榉蔷€性溫度分布,這兩處同為托梁冷翼緣,但溫度發(fā)展并不完全一致,導(dǎo)致與溫度有關(guān)的屈服應(yīng)力值略有差異。

    3.1.3 失效模式

    數(shù)值模擬結(jié)果顯示,托梁未發(fā)生繞弱軸的彎扭屈曲,與試驗(yàn)結(jié)果[8]一致。試驗(yàn)和數(shù)值模擬中,托梁腹板沿長度方向均呈現(xiàn)波浪狀屈曲的失效模式,如圖25(a)所示。托梁J2在支座處均發(fā)生局部屈曲破壞,如圖25(c)所示。圖25證明了經(jīng)數(shù)值模擬的失效模式與試驗(yàn)結(jié)果一致。

    3.2 試件2的模擬分析

    圖26為火災(zāi)下數(shù)值模擬157 min耐火失效時(shí),試件2托梁J2的垂直撓度達(dá)到最大值35.69 mm,此時(shí)支座處和跨中冷翼緣達(dá)到屈服應(yīng)力,結(jié)構(gòu)喪失穩(wěn)定性。失效模式與試件1相同(圖25),均為支座處局部屈曲和沿跨長的腹板屈曲。圖27為數(shù)值模擬、火災(zāi)試驗(yàn)及Baleshan等[18]的托梁撓度隨時(shí)間的變化曲線。

    在熱力耦合分析過程中,全程考慮了背火側(cè)石膏板的約束。然而,試驗(yàn)中背火側(cè)石膏板發(fā)生多處開裂的現(xiàn)象,如圖28所示。該現(xiàn)象影響了托梁冷翼緣處的約束條件,因此,數(shù)值模擬得到的耐火失效時(shí)間長于火災(zāi)試驗(yàn)結(jié)果。

    本文模型在ISO-834標(biāo)準(zhǔn)火災(zāi)下的結(jié)構(gòu)性能與試驗(yàn)結(jié)果整體吻合較好,說明所建立的三維熱力耦合模型可用于分析冷彎薄壁型鋼樓板體系在標(biāo)準(zhǔn)火災(zāi)下的結(jié)構(gòu)性能。

    4 結(jié)論

    在ISO-834標(biāo)準(zhǔn)火災(zāi)條件下,針對兩種不同的冷彎薄壁型鋼承重樓板體系,改進(jìn)了熱力耦合分析中的二維數(shù)值模擬方法,建立了較為精細(xì)的三維傳熱模型,并根據(jù)時(shí)間溫度試驗(yàn)數(shù)據(jù)對模擬結(jié)果進(jìn)行驗(yàn)證。除此之外,為模擬該種結(jié)構(gòu)在標(biāo)準(zhǔn)火災(zāi)條件下背火側(cè)受壓的結(jié)構(gòu)行為,建立了三維熱力耦合模型,數(shù)值模擬經(jīng)與火災(zāi)試驗(yàn)數(shù)據(jù)對比,吻合程度較好。得到以下結(jié)論:

    1)相比二維傳熱模型對空間分析和構(gòu)造影響分析缺失造成的局限性,對冷彎薄壁型鋼樓板體系建立三維傳熱模型,對組成材料石膏板、膠合板、巖棉以及鋼材,合理選取其熱工性能參數(shù),充分考慮了三維空間影響及實(shí)際構(gòu)造等影響因素,更準(zhǔn)確地得到了結(jié)構(gòu)體系溫度分布和升溫規(guī)律。

    2)基于三維數(shù)值模型的熱力耦合分析,在使用經(jīng)驗(yàn)證的傳熱分析數(shù)據(jù)基礎(chǔ)上,顯著改善了結(jié)構(gòu)熱力耦合行為的預(yù)測結(jié)果,包括火災(zāi)條件下托梁的受力彎曲、屈服、失穩(wěn)和考慮非線性溫度分布影響的托梁熱膨脹、熱彎曲。

    3)相比使用較多的二維模型,提出的三維模型更能反映真實(shí)的結(jié)構(gòu)行為,適用范圍更大??梢愿鶕?jù)實(shí)際工程需要,考慮其他種類的板材、組合方式和鋼材截面形狀,拓展到更多的組合樓板結(jié)構(gòu)體系熱力耦合分析中。

    參考文獻(xiàn):

    [1] 韓林海, 宋天詣, 譚清華. 鋼混凝土組合結(jié)構(gòu)抗火設(shè)計(jì)原理研究[J]. 工程力學(xué), 2011, 28(Sup2): 54-66.

    HAN L H, SONG T Y, TAN Q H. Fire resistance design of steel-concrete composite structures [J]. Engineering Mechanics, 2011, 28(Sup2): 54-66. (in Chinese)

    [2] 王廣勇, 韓林海. 局部火災(zāi)下鋼筋混凝土平面框架結(jié)構(gòu)的耐火性能研究[J]. 工程力學(xué), 2010, 27(10): 81-89.

    WANG G Y, HAN L H. Fire performance of reinforced concrete plane frame exposed to local fire [J]. Engineering Mechanics, 2010, 27(10): 81-89. (in Chinese)

    [3] KOLARKAR P N. Structural and thermal performance of cold-formed steel stud wall systems under fire conditions [D]. Queensland: Queensland University of Technology, 2011.

    [4] 葉繼紅, 陳偉, 尹亮. C形冷彎薄壁型鋼承重組合墻體足尺耐火試驗(yàn)研究[J]. 土木工程學(xué)報(bào), 2013, 46(8): 1-10.

    YE J H, CHEN W, YIN L. Full-scale fire resistance tests on load-bearing C-shape cold-formed steel wall systems [J]. China Civil Engineering Journal, 2013, 46(8): 1-10. (in Chinese)

    [5] 余志武, 王中強(qiáng), 蔣麗忠. 火災(zāi)下鋼筋混凝土板的溫度場分析[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2004, 1(1):58-61.

    YU Z W, WANG Z Q, JIANG L Z. Analysis of temperature field in reinforced concrete slabs exposed to fire [J]. Journal of Railway Science and Engineering, 2004, 1(1): 58-61.(in Chinese)

    [6] FENG M, WANG Y C. An experimental study of loaded full-scale cold-formed thin-walled steel structural panels under fire conditions [J]. Fire Safety Journal, 2005, 40(1): 43-63.

    [7] JATHEESHAN V, MAHENDRAN M. Experimental study of cold-formed steel floors made of hollow flange channel section joists under fire conditions [J]. Journal of Structural Engineering, 2016, 142(2): 04015134.

    [8] BALESHAN B, MAHENDRAN M. Experimental study of light gauge steel framing floor systems under fire conditions [J]. Advances in Structural Engineering, 2017, 20(3): 426-445.

    [9] 葉繼紅, 陳偉, 汪正流, 等. 冷彎薄壁C型鋼托梁加氣混凝土板組合樓蓋足尺耐火試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2015, 36(8): 108-115.

    YE J H, CHEN W, WANG Z L, et al. Full scale fire experiments on cold-formed steel C-shape joists-ALC composite floors [J]. Journal of Building Structures, 2015, 36(8): 108-115. (in Chinese)

    [10] 蔣首超, 李國強(qiáng), 周宏宇, 等. 鋼混凝土組合樓蓋抗火性能的試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2004, 25(3): 45-50.

    JIANG S C, LI G Q, ZHOU H Y, et al. Experimental study of behavior of steel-concrete composite slabs subjected to fire [J]. Journal of Building Structures, 2004, 25(3): 45-50.(in Chinese)

    [11] 王衛(wèi)永, 李國強(qiáng), 陳玲珠, 等. 鋼筋桁架樓承板鋼組合梁抗火性能試驗(yàn)研究[J]. 土木工程學(xué)報(bào), 2015, 48(9): 67-75.

    WANG W Y, LI G Q, CHEN L Z, et al. Experimental study on fire resistance of steel bar truss slab and steel composite beams [J]. China Civil Engineering Journal, 2015, 48(9): 67-75. (in Chinese)

    [12] FENG M, WANG Y C, DAVIES J M. Thermal performance of cold-formed thin-walled steel panel systems in fire [J]. Fire Safety Journal, 2003, 38(4): 365-394.

    [13] SULTAN M A. Fire resistance of steel C-joist floor assemblies [J]. Fire Technology, 2010, 46(2): 375-405.

    [14] THOMAS G. Modelling thermal performance of gypsum plasterboard-lined light timber frame walls using SAFIR and TASEF [J]. Fire and Materials, 2010, 34(8): 385-406.

    [15] WANG H B. Heat transfer analysis of components of construction exposed to fire [D]. Salford: University of Salford, 1995: 84-85.

    [16] GUNALAN S, MAHENDRAN M. Finite element modelling of load bearing cold-formed steel wall systems under fire conditions [J]. Engineering Structures, 2013, 56: 1007-1027.

    [17] KEERTHAN P, MAHENDRAN M. Numerical modelling of non-load-bearing light gauge cold-formed steel frame walls under fire conditions [J]. Journal of Fire Sciences, 2012, 30(5): 375-403.

    [18] BALESHAN B, MAHENDRAN M. Numerical study of high strength LSF floor systems in fire [J]. Thin-Walled Structures, 2016, 101: 85-99.

    [19] JATHEESHAN V, MAHENDRAN M. Numerical study of LSF floors made of hollow flange channels in fire [J]. Journal of Constructional Steel Research, 2015, 115: 236-251.

    [20] RUSTHI M, KEERTHAN P, ARIYANAYAGAM A, et al. Numerical studies of gypsum plasterboard and MgO board lined LSF walls exposed to fire [C]//Proceedings of the Second International Conference on Performance-based and Life-cycle Structural Engineering (PLSE 2015). Brisbane, QLD, Australia. School of Civil Engineering, The University of Queensland, 2015: 1077-1084.

    [21] International Standard Organization. Fire-resistance tests-elements of building construction, Part 1: General requirements: ISO-834 [S]. Geneva: International Standard Organization, 1999.

    [22] RUSTHI M, KEERTHAN P, MAHENDRAN M, et al. Investigating the fire performance of LSF wall systems using finite element analyses [J]. Journal of Structural Fire Engineering, 2017, 8(4): 354-376.

    [23] SONG Q Y, WU X B, HAN L H. Investigation on temperature distribution of light gauge steel frame (LSF) wall [C]//9th International Conference on Steel and Aluminium Structures. Bradford: Routledge, 2019: 1447-1458.

    [24] SONG Q Y, WU X B, HAN L H. Modelling of temperature rise of light steel frame walls in fire [C]//2019 International Conference on Advances in Civil Engineering and Materials (ACEM2019) and 2nd World Symposium on Sustainable Bio-composite Materials and Structures (SBMS2). Fuzhou, 2019: 25-27.

    [25] British Standards Institution. Eurocode 1: Design of concrete structures: part 1-2: General actions-actions on structures exposed to fire: BS EN 1991-1-2 [S]. London: British Standards Institution, 2002.

    [26] JATHEESHAN V, MAHENDRAN M. Thermal performance of LSF floors made of hollow flange channel section joists under fire conditions [J]. Fire Safety Journal, 2016, 84: 25-39.

    [27] STERNER E, WICKSTRM U. TASEF-temperature analysis of structures exposed to fire [M]. 1990: 98-99.

    [28] British Standards Institution. Eurocode 3: Design of steel structures: part 1-2: General rules-structural fire design: BS EN 1993-1-2 [S]. London: British Standards Institution, 2005.

    [29] CHEN W, YE J H. Simplified prediction of the thermal and mechanical behavior of a cold-formed steel composite floor at room and elevated temperatures [J]. Journal of Structural Engineering, 2016, 142(6): 04016028.

    [30] British Standards Institution. Eurocode 5: Design of timber structures: part 1-2: General -structural fire design: BS EN 1995-1-2 [S]. London: British Standards Institution, 2004.

    [31] YOUNG S A, CLANCY P. Structural modelling of light-timber framed walls in fire [J]. Fire Safety Journal, 2001, 36(3): 241-268.

    [32] NASSIF A Y, YOSHITAKE I, ALLAM A. Full-scale fire testing and numerical modelling of the transient thermo-mechanical behaviour of steel-stud gypsum board partition walls [J]. Construction and Building Materials, 2014, 59: 51-61.

    [33] KANKANAMGE N D, MAHENDRAN M. Mechanical properties of cold-formed steels at elevated temperatures [J]. Thin-Walled Structures, 2011, 49(1): 26-44.

    [34] 建筑構(gòu)件耐火試驗(yàn)方法: GB/T 9978—1988 [S]. 北京: 中國標(biāo)準(zhǔn)出版社, 2009.

    Fire-resistance tests: Elements of building construction: GB/T 9978-1988 [S]. Beijing: Standards Press of China, 2009. (in Chinese)

    (編輯 王秀玲)

    少妇的逼水好多| 日本vs欧美在线观看视频| 校园人妻丝袜中文字幕| 狂野欧美激情性bbbbbb| 在线观看免费视频网站a站| 在线天堂最新版资源| 亚洲色图 男人天堂 中文字幕| 丰满乱子伦码专区| 国产成人aa在线观看| 老熟女久久久| 久久久国产精品麻豆| 日本午夜av视频| 国产精品欧美亚洲77777| 热99国产精品久久久久久7| 久久久久久人人人人人| 国产成人午夜福利电影在线观看| 日韩人妻精品一区2区三区| av免费在线看不卡| 我要看黄色一级片免费的| 久热久热在线精品观看| 国产伦理片在线播放av一区| 青草久久国产| 两个人看的免费小视频| 999精品在线视频| 久久精品国产亚洲av天美| 亚洲三级黄色毛片| 国产成人精品在线电影| 天堂中文最新版在线下载| 国产精品秋霞免费鲁丝片| 久久精品久久精品一区二区三区| 中文字幕亚洲精品专区| 久久久国产精品麻豆| 国产1区2区3区精品| 蜜桃国产av成人99| 大香蕉久久网| 国产极品天堂在线| 久久精品久久久久久噜噜老黄| 日韩伦理黄色片| 五月伊人婷婷丁香| 亚洲成国产人片在线观看| 国产高清国产精品国产三级| 91精品伊人久久大香线蕉| 美女国产高潮福利片在线看| 黄色 视频免费看| 欧美精品高潮呻吟av久久| 熟女av电影| 国产av精品麻豆| 2022亚洲国产成人精品| 极品人妻少妇av视频| 精品一区在线观看国产| 一二三四中文在线观看免费高清| 五月天丁香电影| 在线观看www视频免费| 在线观看人妻少妇| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人精品久久久久久| 国产精品无大码| 国产精品无大码| 熟女少妇亚洲综合色aaa.| 国产高清国产精品国产三级| 男女边吃奶边做爰视频| 人人澡人人妻人| 搡老乐熟女国产| 制服人妻中文乱码| 好男人视频免费观看在线| 熟妇人妻不卡中文字幕| 美女脱内裤让男人舔精品视频| 99热网站在线观看| av电影中文网址| 免费黄频网站在线观看国产| 老鸭窝网址在线观看| 国产精品 欧美亚洲| 国产成人aa在线观看| 久久毛片免费看一区二区三区| 亚洲一区二区三区欧美精品| 久久青草综合色| 国产欧美日韩一区二区三区在线| 久久久久久久久久久免费av| 各种免费的搞黄视频| 超碰成人久久| 少妇被粗大的猛进出69影院| 中文字幕人妻熟女乱码| 亚洲综合精品二区| 亚洲国产看品久久| 亚洲av电影在线进入| av在线观看视频网站免费| 国产精品秋霞免费鲁丝片| 久久久国产一区二区| 啦啦啦中文免费视频观看日本| 一区二区三区乱码不卡18| 精品少妇一区二区三区视频日本电影 | 亚洲国产欧美在线一区| 女性被躁到高潮视频| 日韩电影二区| 日本色播在线视频| 国产免费一区二区三区四区乱码| 男人添女人高潮全过程视频| 久久99蜜桃精品久久| 天堂8中文在线网| 久热久热在线精品观看| 黄频高清免费视频| 免费av中文字幕在线| 一级黄片播放器| 在线精品无人区一区二区三| 亚洲av男天堂| 一级a爱视频在线免费观看| 久久久欧美国产精品| 国产黄频视频在线观看| 性色avwww在线观看| 国产色婷婷99| 国产精品一国产av| 夜夜骑夜夜射夜夜干| 日韩视频在线欧美| 青草久久国产| 国产在视频线精品| 久久久欧美国产精品| 国产 一区精品| 91精品伊人久久大香线蕉| 久久久国产欧美日韩av| 成人毛片a级毛片在线播放| 国产人伦9x9x在线观看 | 一本—道久久a久久精品蜜桃钙片| 不卡av一区二区三区| 人妻 亚洲 视频| 香蕉国产在线看| 男男h啪啪无遮挡| 波多野结衣av一区二区av| 精品福利永久在线观看| 啦啦啦在线免费观看视频4| 国产日韩欧美亚洲二区| 欧美亚洲 丝袜 人妻 在线| 啦啦啦中文免费视频观看日本| 乱人伦中国视频| 亚洲第一区二区三区不卡| 激情五月婷婷亚洲| av线在线观看网站| 91午夜精品亚洲一区二区三区| 午夜免费鲁丝| 天天影视国产精品| 国产一区二区三区av在线| 91精品伊人久久大香线蕉| 欧美精品高潮呻吟av久久| 亚洲精品自拍成人| 毛片一级片免费看久久久久| 成年女人在线观看亚洲视频| 欧美日韩一级在线毛片| av女优亚洲男人天堂| 深夜精品福利| 菩萨蛮人人尽说江南好唐韦庄| 高清欧美精品videossex| 少妇的逼水好多| 日韩中字成人| 亚洲欧美精品综合一区二区三区 | 国产激情久久老熟女| 国产精品一二三区在线看| www.av在线官网国产| 一级毛片 在线播放| 丝袜脚勾引网站| 国产毛片在线视频| 亚洲成国产人片在线观看| 精品少妇一区二区三区视频日本电影 | 妹子高潮喷水视频| 大香蕉久久网| 伦理电影免费视频| av在线app专区| 亚洲av欧美aⅴ国产| 在线观看美女被高潮喷水网站| 欧美成人午夜免费资源| 亚洲激情五月婷婷啪啪| 丰满少妇做爰视频| 国产毛片在线视频| 18禁裸乳无遮挡动漫免费视频| 精品少妇一区二区三区视频日本电影 | 欧美精品高潮呻吟av久久| 午夜日本视频在线| 婷婷色综合www| 午夜免费男女啪啪视频观看| 一级毛片 在线播放| 秋霞在线观看毛片| 肉色欧美久久久久久久蜜桃| 久久99热这里只频精品6学生| 丁香六月天网| 制服诱惑二区| 久久国内精品自在自线图片| av有码第一页| 欧美精品av麻豆av| av在线老鸭窝| 99热国产这里只有精品6| 2021少妇久久久久久久久久久| 欧美老熟妇乱子伦牲交| 少妇的丰满在线观看| 韩国av在线不卡| 国产男女超爽视频在线观看| 亚洲成国产人片在线观看| 亚洲精品久久午夜乱码| 国产精品无大码| 亚洲精品自拍成人| 最近最新中文字幕免费大全7| 中文字幕av电影在线播放| 国产av国产精品国产| 国产精品国产三级国产专区5o| 国产xxxxx性猛交| 久久久久久久久久人人人人人人| 亚洲人成电影观看| 欧美中文综合在线视频| 99久久综合免费| 婷婷成人精品国产| 国产毛片在线视频| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 人体艺术视频欧美日本| 久久青草综合色| 青春草视频在线免费观看| 欧美人与性动交α欧美精品济南到 | 最近2019中文字幕mv第一页| 叶爱在线成人免费视频播放| 久久久久久久久免费视频了| 两个人免费观看高清视频| 精品国产乱码久久久久久小说| 美女国产高潮福利片在线看| 在线天堂最新版资源| 久久久久久久久免费视频了| 久久精品aⅴ一区二区三区四区 | 91精品国产国语对白视频| 一级黄片播放器| 伦理电影免费视频| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 叶爱在线成人免费视频播放| 一区二区三区精品91| 精品一区在线观看国产| 秋霞在线观看毛片| av视频免费观看在线观看| 美女大奶头黄色视频| 青春草亚洲视频在线观看| 91成人精品电影| 制服诱惑二区| 韩国av在线不卡| 日本免费在线观看一区| 亚洲精品av麻豆狂野| 久久精品熟女亚洲av麻豆精品| 观看av在线不卡| 狂野欧美激情性bbbbbb| 大话2 男鬼变身卡| 丰满少妇做爰视频| 成人免费观看视频高清| 人人澡人人妻人| 男人添女人高潮全过程视频| 中国三级夫妇交换| 999久久久国产精品视频| 亚洲激情五月婷婷啪啪| a级毛片黄视频| 午夜福利在线免费观看网站| 91在线精品国自产拍蜜月| 80岁老熟妇乱子伦牲交| 啦啦啦在线观看免费高清www| 自线自在国产av| 欧美亚洲 丝袜 人妻 在线| 亚洲四区av| 国产成人精品福利久久| 国产一区二区 视频在线| 精品国产一区二区三区四区第35| 桃花免费在线播放| 天天躁夜夜躁狠狠躁躁| 观看av在线不卡| 精品少妇一区二区三区视频日本电影 | 免费黄频网站在线观看国产| 精品一区二区免费观看| 男人舔女人的私密视频| 欧美亚洲 丝袜 人妻 在线| 在线 av 中文字幕| 熟女电影av网| 好男人视频免费观看在线| 伦精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦在线观看免费高清www| www.自偷自拍.com| 最黄视频免费看| 麻豆乱淫一区二区| 视频在线观看一区二区三区| 亚洲欧美中文字幕日韩二区| 久久99一区二区三区| 99九九在线精品视频| 不卡视频在线观看欧美| 人人妻人人添人人爽欧美一区卜| 国产精品国产av在线观看| 精品午夜福利在线看| 男人舔女人的私密视频| 最黄视频免费看| 天堂俺去俺来也www色官网| 成人毛片60女人毛片免费| 咕卡用的链子| 九草在线视频观看| 中文字幕色久视频| 精品国产国语对白av| 久久热在线av| 精品99又大又爽又粗少妇毛片| 亚洲经典国产精华液单| 少妇人妻 视频| 久久国内精品自在自线图片| 国产日韩欧美在线精品| 国产成人精品久久二区二区91 | 色网站视频免费| 精品午夜福利在线看| 免费黄色在线免费观看| 9色porny在线观看| 老汉色∧v一级毛片| 久久久久久久久久人人人人人人| 嫩草影院入口| 90打野战视频偷拍视频| 国产免费视频播放在线视频| 亚洲在久久综合| 久久久精品94久久精品| 久久精品熟女亚洲av麻豆精品| 国产精品亚洲av一区麻豆 | 午夜福利在线免费观看网站| av卡一久久| 国产一区二区三区综合在线观看| 日本av手机在线免费观看| 国产成人午夜福利电影在线观看| 亚洲欧洲精品一区二区精品久久久 | 一个人免费看片子| 亚洲伊人久久精品综合| 伊人久久国产一区二区| 日韩一卡2卡3卡4卡2021年| 中文字幕亚洲精品专区| 视频区图区小说| 久久久久久久大尺度免费视频| 91国产中文字幕| 国产黄色视频一区二区在线观看| 1024视频免费在线观看| 最近2019中文字幕mv第一页| 亚洲av中文av极速乱| 久久久久精品性色| 伊人亚洲综合成人网| 女性被躁到高潮视频| 肉色欧美久久久久久久蜜桃| 日本欧美国产在线视频| 国产女主播在线喷水免费视频网站| 超碰成人久久| 国产成人精品无人区| 日韩免费高清中文字幕av| 性高湖久久久久久久久免费观看| 亚洲国产欧美在线一区| 蜜桃国产av成人99| 啦啦啦啦在线视频资源| 国产xxxxx性猛交| 亚洲av中文av极速乱| videossex国产| 国产精品不卡视频一区二区| 一区二区av电影网| 97精品久久久久久久久久精品| 国产成人精品在线电影| 国产精品女同一区二区软件| 亚洲国产看品久久| 视频区图区小说| 九色亚洲精品在线播放| 看十八女毛片水多多多| 亚洲美女视频黄频| 亚洲欧美清纯卡通| 欧美人与善性xxx| 18禁裸乳无遮挡动漫免费视频| 最近的中文字幕免费完整| 亚洲欧美成人精品一区二区| 亚洲av成人精品一二三区| 下体分泌物呈黄色| 熟女av电影| 菩萨蛮人人尽说江南好唐韦庄| 在线 av 中文字幕| 亚洲成av片中文字幕在线观看 | 69精品国产乱码久久久| 精品卡一卡二卡四卡免费| 日韩一区二区视频免费看| 成人黄色视频免费在线看| 99久久精品国产国产毛片| 综合色丁香网| 欧美成人精品欧美一级黄| 叶爱在线成人免费视频播放| 欧美少妇被猛烈插入视频| av电影中文网址| 国产av一区二区精品久久| 少妇人妻 视频| 国产精品99久久99久久久不卡 | 91精品三级在线观看| 青青草视频在线视频观看| 久久av网站| 青青草视频在线视频观看| 亚洲精品美女久久久久99蜜臀 | 搡老乐熟女国产| 日韩在线高清观看一区二区三区| 黄色怎么调成土黄色| 亚洲精品美女久久av网站| 久久久久人妻精品一区果冻| av在线老鸭窝| 久久久久久久久久人人人人人人| 免费高清在线观看视频在线观看| 亚洲国产最新在线播放| 成年女人在线观看亚洲视频| 欧美日韩精品成人综合77777| 侵犯人妻中文字幕一二三四区| 九草在线视频观看| 国产高清国产精品国产三级| 男人添女人高潮全过程视频| 一区在线观看完整版| 亚洲成人手机| 亚洲,欧美精品.| 精品午夜福利在线看| 纯流量卡能插随身wifi吗| 男女边摸边吃奶| 一区二区三区四区激情视频| 亚洲视频免费观看视频| 激情五月婷婷亚洲| 午夜福利乱码中文字幕| 久久久久久伊人网av| 三上悠亚av全集在线观看| 2022亚洲国产成人精品| 下体分泌物呈黄色| 人体艺术视频欧美日本| 国产成人免费无遮挡视频| 看免费成人av毛片| 国产亚洲最大av| 国产亚洲av片在线观看秒播厂| 香蕉国产在线看| 午夜日韩欧美国产| 成人亚洲精品一区在线观看| 黑人欧美特级aaaaaa片| 亚洲国产毛片av蜜桃av| 亚洲欧美成人综合另类久久久| 亚洲三区欧美一区| 免费高清在线观看视频在线观看| a 毛片基地| 欧美日韩综合久久久久久| av国产精品久久久久影院| 赤兔流量卡办理| 国产成人精品婷婷| 亚洲欧美色中文字幕在线| 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 男人操女人黄网站| 三级国产精品片| 1024视频免费在线观看| 亚洲精品第二区| 欧美亚洲 丝袜 人妻 在线| 国产精品一二三区在线看| 久久精品亚洲av国产电影网| 青春草国产在线视频| 欧美亚洲 丝袜 人妻 在线| 国产精品一二三区在线看| 久久免费观看电影| 久久久久精品性色| 一级,二级,三级黄色视频| 国产片特级美女逼逼视频| 免费高清在线观看日韩| 一级毛片我不卡| 中文字幕另类日韩欧美亚洲嫩草| 少妇猛男粗大的猛烈进出视频| 国产一级毛片在线| 久久久欧美国产精品| 啦啦啦在线观看免费高清www| 精品亚洲成国产av| av片东京热男人的天堂| 天天躁夜夜躁狠狠躁躁| 免费观看a级毛片全部| 国产乱来视频区| 最新的欧美精品一区二区| 人妻少妇偷人精品九色| 亚洲国产精品一区二区三区在线| 一级毛片我不卡| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 国产精品.久久久| 国产成人免费无遮挡视频| 岛国毛片在线播放| 欧美日韩视频高清一区二区三区二| 一区二区日韩欧美中文字幕| 一本—道久久a久久精品蜜桃钙片| 日本91视频免费播放| 性少妇av在线| 久久久久久人人人人人| 香蕉国产在线看| 国产精品一二三区在线看| 午夜日本视频在线| 一二三四在线观看免费中文在| 99精国产麻豆久久婷婷| 国产精品久久久久久av不卡| 亚洲中文av在线| 高清不卡的av网站| 免费大片黄手机在线观看| 校园人妻丝袜中文字幕| 一级毛片黄色毛片免费观看视频| 一级黄片播放器| 免费播放大片免费观看视频在线观看| 欧美97在线视频| 我的亚洲天堂| 肉色欧美久久久久久久蜜桃| 国产精品一二三区在线看| 极品人妻少妇av视频| 亚洲国产毛片av蜜桃av| av网站在线播放免费| 亚洲精品久久成人aⅴ小说| 免费久久久久久久精品成人欧美视频| 精品国产一区二区三区四区第35| 欧美日韩亚洲国产一区二区在线观看 | 久久青草综合色| 免费在线观看完整版高清| 一级毛片电影观看| h视频一区二区三区| 国产亚洲欧美精品永久| 青春草视频在线免费观看| 国产成人精品婷婷| 亚洲国产看品久久| 日韩av免费高清视频| 成年动漫av网址| 男女高潮啪啪啪动态图| 午夜免费男女啪啪视频观看| 成人18禁高潮啪啪吃奶动态图| 深夜精品福利| 欧美人与性动交α欧美软件| 久久99精品国语久久久| 老女人水多毛片| 成人漫画全彩无遮挡| 老司机影院毛片| 超碰成人久久| 波多野结衣一区麻豆| 午夜影院在线不卡| 久久青草综合色| 久久99一区二区三区| 如何舔出高潮| 免费久久久久久久精品成人欧美视频| 黄色一级大片看看| 国产 精品1| 国产国语露脸激情在线看| 在线观看免费视频网站a站| 久久 成人 亚洲| 伊人亚洲综合成人网| 国产成人欧美| 日韩不卡一区二区三区视频在线| 精品人妻一区二区三区麻豆| 青草久久国产| 在线观看三级黄色| 新久久久久国产一级毛片| 性少妇av在线| 国产一区有黄有色的免费视频| 亚洲av福利一区| 日韩视频在线欧美| 黑人猛操日本美女一级片| 亚洲第一青青草原| av国产精品久久久久影院| 女人久久www免费人成看片| 久久精品国产综合久久久| 欧美日韩一区二区视频在线观看视频在线| av不卡在线播放| 麻豆av在线久日| 蜜桃国产av成人99| 日本wwww免费看| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品一区二区在线不卡| 香蕉丝袜av| 性色av一级| 国产一区亚洲一区在线观看| 看免费成人av毛片| 最近2019中文字幕mv第一页| 极品人妻少妇av视频| 亚洲,欧美精品.| 亚洲四区av| 欧美bdsm另类| 热99久久久久精品小说推荐| av国产久精品久网站免费入址| 2022亚洲国产成人精品| 欧美av亚洲av综合av国产av | 99久国产av精品国产电影| 亚洲熟女精品中文字幕| 久久婷婷青草| 精品国产国语对白av| 国产精品偷伦视频观看了| 波多野结衣av一区二区av| 久久99热这里只频精品6学生| 成人毛片a级毛片在线播放| 久久精品国产a三级三级三级| 久久ye,这里只有精品| av.在线天堂| 中文字幕制服av| 伦理电影免费视频| 免费女性裸体啪啪无遮挡网站| 欧美成人精品欧美一级黄| 免费av中文字幕在线| 国产熟女欧美一区二区| 成人二区视频| 国产精品三级大全| 99国产综合亚洲精品| 婷婷色av中文字幕| 熟女少妇亚洲综合色aaa.| 欧美最新免费一区二区三区| 久久久久久免费高清国产稀缺| 久久久国产一区二区| 国产男女超爽视频在线观看| 亚洲成人一二三区av| 国产亚洲午夜精品一区二区久久| 最近中文字幕2019免费版| 久热久热在线精品观看| 男女高潮啪啪啪动态图| 免费观看在线日韩| 国产欧美日韩一区二区三区在线| 成人影院久久| 久久99一区二区三区| 久久精品国产鲁丝片午夜精品| 国产综合精华液| 国产精品 国内视频| 一个人免费看片子| 国产1区2区3区精品| 大香蕉久久成人网| 在线观看免费视频网站a站| 美国免费a级毛片| 26uuu在线亚洲综合色| 精品人妻在线不人妻|