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

    聚脲噴涂鋼制罐體抗爆性能試驗(yàn)及數(shù)值模擬研究

    2023-02-24 02:16:30楊光張博一韋建樹(shù)蔣月新徐世林
    關(guān)鍵詞:抗爆聚脲模擬計(jì)算

    楊光,張博一,韋建樹(shù),蔣月新,徐世林

    (1.上海爵格工業(yè)工程有限公司,上海 200082;2.哈爾濱工業(yè)大學(xué) 土木工程學(xué)院,哈爾濱 150001;3.哈爾濱工業(yè)大學(xué)(深圳)土木與環(huán)境工程學(xué)院,廣東 深圳 518055)

    近年來(lái),石油化工罐區(qū)逐漸成為恐怖襲擊的目標(biāo)。2020年1月,敘利亞巴尼亞斯儲(chǔ)油罐區(qū)就曾遭到恐怖分子爆炸襲擊,造成超過(guò)10億美元的直接經(jīng)濟(jì)損失。有效降低爆炸造成的損失、控制災(zāi)害影響范圍一直是業(yè)內(nèi)關(guān)注的焦點(diǎn),其中,大型鋼制儲(chǔ)油罐爆炸因具有巨大的危害性而受到重點(diǎn)關(guān)注。當(dāng)罐區(qū)某一位置遭到爆炸襲擊時(shí),產(chǎn)生的沖擊波會(huì)迅速波及附近儲(chǔ)油罐,如果儲(chǔ)油罐抗爆能力不足,則可能發(fā)生結(jié)構(gòu)破壞,進(jìn)而導(dǎo)致燃油泄漏,并引發(fā)二次爆炸乃至連環(huán)爆炸,從而造成巨大的經(jīng)濟(jì)損失和人員傷亡。

    噴涂聚脲彈性體(PSUA)是近20年來(lái)研制開(kāi)發(fā)的一種新型綠色施工技術(shù)[1],其原料——聚脲是一種高分子彈性體,具有耐磨、防水、抗疲勞等優(yōu)點(diǎn),廣泛應(yīng)用于工業(yè)涂裝領(lǐng)域。大量研究表明,聚脲涂層在結(jié)構(gòu)抗沖擊方面也有著優(yōu)異的性能[2-4]。Davidson等[5]對(duì)聚脲噴涂的砌體墻進(jìn)行了試驗(yàn)及數(shù)值模擬分析,研究表明,在墻體背面噴涂聚脲能更有效地提升墻體的抗爆性能。Raman等[6]則對(duì)聚脲涂覆鋼筋混凝土板的抗爆能力進(jìn)行了數(shù)值模擬研究,結(jié)果表明,聚脲涂層能夠?qū)⒈ㄗ饔孟落摻罨炷涟宓奈灰祈憫?yīng)減小40%。將高分子彈性體用于金屬結(jié)構(gòu)抗爆能力提升的相關(guān)研究也取得了較為顯著的成果[7-12]。Xue等[13]通過(guò)試驗(yàn)與數(shù)值分析發(fā)現(xiàn),噴涂高分子彈性體的雙層金屬板結(jié)構(gòu)存在頸縮遲滯效應(yīng)(neck retardation),使雙層板結(jié)構(gòu)的頸縮極限顯著提高,從而增強(qiáng)了其承受沖擊荷載的能力;Mcshane等[14]針對(duì)Xue等[13]的理論進(jìn)行聚合物—金屬雙層板動(dòng)態(tài)斷裂研究發(fā)現(xiàn),在拉伸能量的吸收方面,足夠堅(jiān)硬和堅(jiān)固的聚合物涂層對(duì)結(jié)構(gòu)存在一定的增強(qiáng)作用;Jiang等[15]對(duì)聚脲復(fù)合鋼板在低速?zèng)_擊下的性能進(jìn)行試驗(yàn)研究發(fā)現(xiàn),在低速?zèng)_擊作用下,聚脲涂層能夠有效地減小鋼板的塑性變形,且當(dāng)涂層厚度較薄時(shí)效果顯著,而厚度增加帶來(lái)的強(qiáng)度提升并不明顯。目前,現(xiàn)有研究多局限于單個(gè)板構(gòu)件(如聚脲—金屬雙層板)的局部承載研究,而對(duì)于完整金屬結(jié)構(gòu)物的爆炸試驗(yàn)及數(shù)值研究較少,對(duì)于聚脲應(yīng)用于整體金屬結(jié)構(gòu)的抗爆性能提升也缺乏較明確的論據(jù)。

    筆者針對(duì)某大型鋼制儲(chǔ)罐的抗爆性能進(jìn)行相似結(jié)構(gòu)模型的爆炸試驗(yàn),通過(guò)試驗(yàn)數(shù)據(jù)的比較,研究聚脲涂層對(duì)薄壁鋼制罐體結(jié)構(gòu)抗爆能力的提升作用,并結(jié)合試驗(yàn)?zāi)P瓦M(jìn)行有限元分析與對(duì)比。

    1 爆炸試驗(yàn)

    1.1 試驗(yàn)設(shè)計(jì)

    爆炸試驗(yàn)在陸軍工程大學(xué)室外爆炸試驗(yàn)場(chǎng)地進(jìn)行,共計(jì)2次,炸藥為固體TNT炸藥,質(zhì)量2 kg,形狀為圓柱體,如圖1所示。

    圖1 試驗(yàn)用TNT炸藥Fig.1 TNT explosive in test

    試驗(yàn)?zāi)康氖强疾榫垭逋繉訉?duì)大型儲(chǔ)罐建造項(xiàng)目中雙層鋼制儲(chǔ)罐抗爆能力的提升,所以,試驗(yàn)所用模型罐是依據(jù)相似原理對(duì)實(shí)際儲(chǔ)罐外罐進(jìn)行比例縮放得到的,比例為1:26.7??s放后的試驗(yàn)罐體模型結(jié)構(gòu)分為拱頂、圓柱形罐壁和法蘭3部分,其中,拱頂及罐壁厚度為1 mm,法蘭厚度為2 mm,罐體材料為Q235低碳鋼,模型罐尺寸如圖2所示。

    圖2 模型罐體幾何尺寸簡(jiǎn)圖(單位:mm)Fig.2 Geometric dimension sketch of tank model(Unit:mm)

    罐體不同構(gòu)件之間采用焊接,并在法蘭中線(xiàn)處均勻布置32個(gè)螺栓孔,通過(guò)打入地面的膨脹螺栓實(shí)現(xiàn)對(duì)罐體的約束,如圖3(a)所示。聚脲涂層均勻噴涂在罐體內(nèi)表面,噴涂區(qū)域涵蓋整個(gè)罐壁及拱頂,不包括法蘭,如圖3(b)所示。在罐體迎爆側(cè)距底面高度275 mm處布置位移測(cè)點(diǎn),測(cè)點(diǎn)處安裝拉桿式直線(xiàn)位移傳感器,以采集罐體表面位移時(shí)程數(shù)據(jù),如圖3(c)、(d)所示。罐體殘余位移則采用卷尺手動(dòng)測(cè)量,如圖3(e)所示。

    圖3 試驗(yàn)設(shè)施Fig.3 Test facility

    試驗(yàn)共進(jìn)行2次爆炸,分別考查在不同爆炸距離下無(wú)涂層鋼罐及噴涂聚脲涂層鋼罐的動(dòng)力響應(yīng),炸藥及鋼罐布置如圖4所示。

    圖4 試驗(yàn)炸藥及鋼罐布置Fig.4 Layout of TNT and steel tanks in test

    1)第1次爆炸,無(wú)涂層鋼罐與噴涂層鋼罐分立炸藥兩側(cè),炸藥中心與儲(chǔ)罐迎爆面外壁距離為4 m;

    2)第2次爆炸,無(wú)涂層鋼罐與噴涂層鋼罐分立炸藥兩側(cè),炸藥中心與儲(chǔ)罐迎爆面外壁距離為3.2 m。

    1.2 試驗(yàn)結(jié)果

    1)第1次爆炸,罐體表面殘余位移很小,絕大部分變形處于彈性范圍內(nèi),如圖5所示。

    從圖5可以看出,炸藥距離罐體4 m時(shí),爆炸造成的破壞較小,罐體結(jié)構(gòu)響應(yīng)以彈性變形為主,很難從爆炸后罐體形態(tài)上直觀(guān)地體現(xiàn)出涂層的作用。無(wú)涂層罐體和噴涂層罐體的測(cè)點(diǎn)位移如圖6所示,噴涂層罐體測(cè)點(diǎn)位置最大位移為26.24 mm,而無(wú)涂層罐體測(cè)點(diǎn)位置最大位移為39.13 mm,噴涂1.5 mm聚脲涂層可使罐體測(cè)點(diǎn)位置的最大位移降低32.9%。另外,通過(guò)對(duì)罐體表面殘余位移的測(cè)量發(fā)現(xiàn),噴涂層罐體的殘余位移為0.42 mm,而無(wú)涂層罐體的殘余位移為4.33 mm,內(nèi)部噴涂聚脲涂層對(duì)于降低結(jié)構(gòu)殘余位移也有很明顯的效果。

    圖5 第1次爆炸后的罐體變形Fig.5 Tank deformation after first explosion

    圖6 第1次爆炸時(shí)罐體測(cè)點(diǎn)位移時(shí)程曲線(xiàn)Fig.6 Displacement time-history curve of tank test point during first explosion

    2)第2次爆炸,罐體表面有明顯的殘余位移,迎爆面鋼板已進(jìn)入塑性范圍,如圖7所示。

    圖7 第2次爆炸后的罐體變形Fig.7 Tank deformation after second explosion

    從圖7可以看出,炸藥距離罐體3.2 m時(shí),爆炸造成的破壞較大,此時(shí)噴涂層罐體爆炸后的形態(tài)明顯好于無(wú)涂層罐體。無(wú)涂層和噴涂層罐體的測(cè)點(diǎn)位移如圖8所示,噴涂層罐體測(cè)點(diǎn)位置最大位移為44.27 mm,而無(wú)涂層罐體測(cè)點(diǎn)位置最大位移為63.6 mm,噴涂聚脲涂層可使罐體測(cè)點(diǎn)位置的最大位移降低30.4%。另外,對(duì)罐體表面殘余位移進(jìn)行測(cè)量,噴涂層罐體的殘余位移為16.51 mm,而無(wú)涂層罐體的殘余位移為24.47 mm,說(shuō)明噴涂聚脲涂層對(duì)結(jié)構(gòu)殘余位移的降低也有很明顯的效果。

    圖8 第2次爆炸時(shí)罐體測(cè)點(diǎn)位移時(shí)程曲線(xiàn)Fig.8 Displacement time-history curve of tank test point during second explosion

    從兩次試驗(yàn)的結(jié)果可以看出,在炸藥不同距離的爆炸作用下,噴涂層罐體的測(cè)點(diǎn)最大位移與殘余位移均較無(wú)涂層罐體有明顯的降低,降幅約為30%。說(shuō)明無(wú)論結(jié)構(gòu)響應(yīng)處于彈性區(qū)間還是塑性區(qū)間,聚脲涂層對(duì)結(jié)構(gòu)強(qiáng)度的提升都有明顯效果。

    2 數(shù)值模擬

    ALE法(任意拉格朗日—?dú)W拉方法)作為一種流固耦合計(jì)算方法,能實(shí)現(xiàn)空氣與罐體網(wǎng)格間力的傳遞,相較于直接加載壓力時(shí)程曲線(xiàn)或者LBE法,其優(yōu)點(diǎn)在于整個(gè)模擬過(guò)程滿(mǎn)足能量守恒,且能考慮罐體與空氣域的相互作用,缺點(diǎn)在于空氣網(wǎng)格數(shù)量龐大,計(jì)算耗時(shí)較長(zhǎng)。應(yīng)用LS-DYNA軟件,采用ALE法分析罐體在炸藥自由爆炸沖擊波作用下的響應(yīng)。

    2.1 有限元模型

    與試驗(yàn)過(guò)程相同,數(shù)值計(jì)算也分爆炸距離為3.2、4 m兩組工況,并分別與試驗(yàn)結(jié)果進(jìn)行比較。由于鋼罐體屬于薄壁結(jié)構(gòu),采用2D薄殼單元建模較為合適,既可以控制網(wǎng)格數(shù)量,又能保證計(jì)算精度。采用Belytschko-Tsay shell單元對(duì)罐體、法蘭、涂層及地面建模,空氣及炸藥選用六面體Solid單元,為保證流體單元與結(jié)構(gòu)單元良好耦合,兩者尺寸應(yīng)盡量接近。ALE算法的計(jì)算結(jié)果受網(wǎng)格尺寸影響較大,網(wǎng)格尺寸越小,峰值超壓越大,計(jì)算結(jié)果越精確。師燕超[16]和都浩等[17]通過(guò)研究發(fā)現(xiàn),當(dāng)空氣網(wǎng)格小于100 mm時(shí),得到的沖擊波超壓已經(jīng)可以滿(mǎn)足精度要求。計(jì)算模型的基本網(wǎng)格尺寸選為30 mm,以4 m無(wú)涂層爆炸工況進(jìn)行敏感性分析。相對(duì)于20 mm網(wǎng)格尺寸(40.34 mm),該尺寸網(wǎng)格得到的測(cè)點(diǎn)位移響應(yīng)峰值(38.6 mm)的誤差約為4.5%,說(shuō)明30 mm網(wǎng)格已能夠滿(mǎn)足基本精度要求,并可以節(jié)省大量計(jì)算時(shí)間。由于計(jì)算工況中的模型及荷載均滿(mǎn)足對(duì)稱(chēng)特征,為節(jié)省計(jì)算資源,采用1/2模型,并在對(duì)稱(chēng)面施加法向約束。罐體與法蘭單元采用共節(jié)點(diǎn)連接,地面采用剛性地面,法蘭與地面間采用ATUOMATIC_SURFACE_TO_SURFACE接觸,忽略摩擦力,僅考慮接觸面間的法向力,在螺栓實(shí)際位置處約束相鄰4個(gè)節(jié)點(diǎn),以實(shí)現(xiàn)螺栓的緊固作用??諝庥蜻吔绮捎肗ON_REFLECTING無(wú)反射邊界條件。有限元模型如圖9所示(空氣網(wǎng)格不顯示)。

    圖9 有限元模型示意圖Fig.9 Sketch of FE model

    數(shù)值模擬采用mm-ton-s單位制,TNT炸藥選用MAT_HIGH_EXPLOSIVE_BURN材料模型以及JWL狀態(tài)方程,空氣選用MAT_NULL材料模型以及線(xiàn)性多項(xiàng)式狀態(tài)方程,兩種狀態(tài)方程表述為

    JWL:

    式中:p為壓強(qiáng),MPa;E為炸藥內(nèi)能,J·m-3;V為當(dāng)前相對(duì)體積,其余各參數(shù)見(jiàn)表1。

    表1 炸藥材料及狀態(tài)方程參數(shù)[18]Table 1 Parameters of TNT material and EOS

    線(xiàn)性多項(xiàng)式:

    表2 線(xiàn)性多項(xiàng)式狀態(tài)方程參數(shù)[19]Table 2 Parameters of linear polynomial EOS

    罐體鋼材采用MAT_SIMPLIFIED_JOHNSON_COOK本構(gòu)模型,該模型能夠考慮到材料在高應(yīng)變率下的應(yīng)力強(qiáng)化及后屈服狀態(tài)下的塑形應(yīng)變,廣泛應(yīng)用于沖擊與爆炸領(lǐng)域的分析計(jì)算。為得到試驗(yàn)中罐體所用鋼材準(zhǔn)確的材料參數(shù),對(duì)其進(jìn)行了準(zhǔn)靜態(tài)拉伸、壓縮和動(dòng)態(tài)壓縮試驗(yàn)[20],通過(guò)試驗(yàn)數(shù)據(jù)擬合得到的各材料參數(shù)見(jiàn)表3。

    表3 罐體材料參數(shù)Table 3 Parameters of steel tank material

    聚脲材料采用MAT_PLASTICITY_POLYMER本構(gòu)模型,該模型依據(jù)Cowper-Symonds理論,通過(guò)動(dòng)力放大系數(shù)(DIF)的形式來(lái)表現(xiàn)材料的應(yīng)變率強(qiáng)化效應(yīng)。同樣地,為得到聚脲材料不同應(yīng)變率條件下的應(yīng)力—應(yīng)變關(guān)系,進(jìn)行準(zhǔn)靜態(tài)拉伸試驗(yàn)及霍普金斯壓桿試驗(yàn)[7]。通過(guò)對(duì)試驗(yàn)數(shù)據(jù)的擬合,得到該聚脲材料的各材料參數(shù),見(jiàn)表4。

    表4 聚脲材料參數(shù)Table 4 Parameters of ployurea coating

    對(duì)于聚脲涂層與鋼罐之間的連接,采用CONTACT_TIEBREAK_SURFACE_TO_SURFACE接 觸,它是一種帶有失效模式的綁定接觸,其失效判定主要由拉伸失效應(yīng)力(NFLS)和剪切失效應(yīng)力(SFLS)決定,具體判定準(zhǔn)則為

    式中:σn為單元拉伸應(yīng)力;σs為單元剪切應(yīng)力。

    試驗(yàn)中所用聚脲材料與鋼板連接的NFLS和SFLS由準(zhǔn)靜態(tài)試驗(yàn)得到,均為7.94 MPa。

    2.2 計(jì)算結(jié)果

    無(wú)涂層及有涂層罐體在兩次爆炸后的狀態(tài)分別見(jiàn)圖10、圖11。

    圖10 3.2 m距離爆炸在200 ms時(shí)刻的爆炸計(jì)算位移云圖Fig.10 Cloud image of caltulated explosion displacement of 3.2 m blast at 200 ms

    當(dāng)計(jì)算進(jìn)行到200 ms時(shí),沖擊波作用已結(jié)束,罐體振蕩逐漸衰減,該時(shí)刻的位移云圖接近于爆炸結(jié)束后罐體的最終狀態(tài)。從圖10、圖11中的位移云圖來(lái)看,模擬計(jì)算能夠較為直觀(guān)地反映涂層對(duì)降低罐體變形幅度的影響,與試驗(yàn)照片對(duì)比發(fā)現(xiàn),迎爆面凹陷程度和褶皺分布也較為接近,但無(wú)涂層罐體的頂面凹陷幅度較大,與試驗(yàn)結(jié)果有一定的差距。分析原因,可能是由于拱頂與罐壁連接處屬于結(jié)構(gòu)突變位置,在受力時(shí)發(fā)生折疊而使材料加速進(jìn)入塑性階段,轉(zhuǎn)角處無(wú)法繼續(xù)提供反向彎矩,導(dǎo)致位移響應(yīng)過(guò)大。試驗(yàn)中的罐體拱頂與罐壁采用焊接連接,焊角在一定程度上緩解了該位置處的折疊變形,從而降低了位移響應(yīng)。

    圖11 4 m距離爆炸在200 ms時(shí)刻的爆炸計(jì)算位移云圖Fig.11 Cloud image of calculated explosion displacement of 4 m blast at 200 ms

    對(duì)于測(cè)點(diǎn)位置處的位移,將模擬計(jì)算時(shí)程曲線(xiàn)與試驗(yàn)數(shù)據(jù)進(jìn)行比較,如圖12所示。從測(cè)點(diǎn)位移曲線(xiàn)的對(duì)比可知,無(wú)論是有涂層還是無(wú)涂層,ALE計(jì)算方法得到的兩次爆炸罐體測(cè)點(diǎn)位移峰值與試驗(yàn)值均較為接近,具體數(shù)據(jù)對(duì)照見(jiàn)表5。

    表5 測(cè)點(diǎn)位移計(jì)算結(jié)果與試驗(yàn)數(shù)據(jù)對(duì)照Table 5 Comparation between calculating results and test data of test points

    從表5中可以看出,采用ALE法計(jì)算得到的測(cè)點(diǎn)位移響應(yīng)與試驗(yàn)值較為接近,最大誤差不超過(guò)20%,相比用LBE法模擬得到的結(jié)果[18],準(zhǔn)確性更高,具有一定的應(yīng)用價(jià)值。模擬計(jì)算結(jié)果中,有涂層在4 m時(shí)對(duì)于罐體測(cè)點(diǎn)位移的降低幅度為34.1%,3.2 m時(shí)降低幅度為21.2%,與試驗(yàn)得到的32.9%和30.4%較為符合,說(shuō)明模擬計(jì)算中采用的聚脲材料參數(shù)與實(shí)際值貼近度較好,計(jì)算得到的結(jié)果具有較高的參考價(jià)值。但從圖12中也能看到,模擬計(jì)算得到的測(cè)點(diǎn)殘余位移較大,遠(yuǎn)高于試驗(yàn)值,說(shuō)明雖然罐體位移大致相同,但計(jì)算模型中的鋼材塑性變形程度更深,材料的屈服強(qiáng)度較實(shí)際情況偏小,具體原因還需要通過(guò)進(jìn)一步的材料測(cè)試分析得到。

    圖12 爆炸模擬測(cè)點(diǎn)位移結(jié)果與試驗(yàn)數(shù)據(jù)對(duì)比Fig.12 Comparison between simulating results and test data of test point displacement

    進(jìn)一步比較聚脲涂層對(duì)罐體最大位移響應(yīng)的影響,兩次測(cè)試模擬結(jié)果見(jiàn)圖13及表6。與測(cè)點(diǎn)位移結(jié)果相似,聚脲涂層也能顯著地降低罐體最大位移,4 m時(shí)可降低20.5%,3.2 m時(shí)可降低24.0%。

    表6 罐體最大位移計(jì)算值對(duì)照Table 6 Comparison between the calculated maximum displacement of tank

    圖13 模擬計(jì)算罐體最大位移時(shí)程曲線(xiàn)Fig.13 Time-history curves of max displacement of tank in simulated calculation

    比較模擬計(jì)算結(jié)果和試驗(yàn)測(cè)值可知,采用ALE法對(duì)鋼罐體在TNT爆炸荷載作用下的響應(yīng)進(jìn)行模擬具有一定的精確度,模擬得到的測(cè)點(diǎn)位移峰值與試驗(yàn)結(jié)果基本一致,但殘余位移值大于試驗(yàn)結(jié)果。模擬計(jì)算結(jié)果表明,聚脲涂層能夠較為明顯地降低鋼罐體在爆炸荷載下的位移響應(yīng),降低幅度約為20%~35%。

    2.3 變參數(shù)分析

    2.3.1 粘接設(shè)置對(duì)涂層防護(hù)效果的影響 常規(guī)的涂層與鋼板連接方式模擬可以用帶有失效模式的TIEBREAK接觸,也可以不考慮連接失效,采用TIED綁定接觸,甚至直接將兩者單元共節(jié)點(diǎn)處理。區(qū)別在于,TIED綁定接觸無(wú)法考慮涂層與鋼板連接失效,而單元共節(jié)點(diǎn)既不能考慮連接失效,也不能體現(xiàn)涂層和鋼板的相對(duì)位置,它的優(yōu)點(diǎn)在于設(shè)置簡(jiǎn)單,同時(shí),節(jié)點(diǎn)數(shù)量的縮減也可以節(jié)省計(jì)算資源。

    前文根據(jù)涂層鋼板實(shí)測(cè)數(shù)據(jù),采用TIEBREAK接觸形式進(jìn)行模擬,接下來(lái)進(jìn)行兩組試算,比較3種不同接觸形式對(duì)模擬產(chǎn)生的影響。對(duì)照工況選用響應(yīng)較為明顯的3.2 m爆炸工況,分別考查測(cè)點(diǎn)位移曲線(xiàn)與最大位移曲線(xiàn),結(jié)果見(jiàn)圖14。

    從圖14中可以 看 出,TIEBREAK和TIED接觸對(duì)于位移結(jié)果幾乎沒(méi)有影響,兩條曲線(xiàn)峰值幾乎完全重合,在后期往復(fù)震蕩時(shí)才有少許偏差。相對(duì)而言,共節(jié)點(diǎn)方式和前兩者的位移峰值差距小于3%,尚可接受,但在隨后的振蕩中,差距不斷放大,最大可達(dá)30%。

    圖14 不同接觸形式下最大位移時(shí)程曲線(xiàn)對(duì)照Fig.14 Comparison between time-history curves of maximum displacement of tank under different contact types

    試驗(yàn)結(jié)果中,雖然鋼板存在較大的位移與殘余變形,但涂層與鋼板連接良好,并無(wú)脫膠情況發(fā)生,說(shuō)明以現(xiàn)有的噴涂技術(shù),可以保證即使鋼結(jié)構(gòu)有輕微破壞,也很難發(fā)生連接失效。所以,在缺少實(shí)測(cè)數(shù)據(jù)的情況下,可以先假定兩者連接不會(huì)失效而采用TIED綁定接觸進(jìn)行模擬,這并不會(huì)影響結(jié)果的精確程度。共節(jié)點(diǎn)連接方式僅能用來(lái)考查峰值位移,并不推薦在其他情況下應(yīng)用。

    2.3.2 罐體鋼材對(duì)涂層防護(hù)效果的影響 雖然針對(duì)噴涂聚脲涂層對(duì)Q235低碳鋼制罐體抗爆能力的提升進(jìn)行了分析,而工業(yè)生產(chǎn)中常常會(huì)使用強(qiáng)度更高的鋼材,如Q345低碳合金鋼、45#碳素結(jié)構(gòu)鋼等。針對(duì)這3種常用鋼材進(jìn)行試算比較,考查聚脲涂層對(duì)于不同鋼材抗爆能力的影響。

    試驗(yàn)中僅用到Q235鋼材,沒(méi)有Q345和45#鋼材的實(shí)測(cè)數(shù)據(jù),在模擬計(jì)算中,參考于文靜等[21]和陳剛等[22]的數(shù)據(jù)對(duì)這兩種鋼材進(jìn)行材料參數(shù)設(shè)置,具體參數(shù)見(jiàn)表7。

    表7 鋼板材料參數(shù)Table 7 Parameters of steel plates

    同樣選取3.2 m爆炸距離作為對(duì)照工況,比較不同鋼材罐體在無(wú)涂層及噴涂層時(shí)的最大位移時(shí)程曲線(xiàn),結(jié)果如圖15及表8所示。

    表8 不同鋼材罐體最大位移計(jì)算值對(duì)照Table 8 Comparison between maximum displacements of tanks with different steels

    圖15 不同鋼材罐體最大位移時(shí)程曲線(xiàn)對(duì)照Fig.15 Comparison between time-history curves of maximum displacement of different steel tanks

    從無(wú)涂層罐體的位移曲線(xiàn)中看出,3種鋼材的強(qiáng)度順序是45#>Q345>Q235,而三者噴涂聚脲后的抗爆強(qiáng)度都有一定程度提升,但提升的幅度有所不同。Q235鋼板噴涂聚脲能夠減少24%的最大位移,而同樣的方式下Q345鋼僅減小11.37%,45#鋼僅減小6.45%。隨著鋼材強(qiáng)度的提高,聚脲涂層帶來(lái)的強(qiáng)度增益有所減小,但該趨勢(shì)是否適用于一般情況,還需進(jìn)一步試驗(yàn)研究。

    2.3.3 涂層厚度對(duì)防護(hù)效果的影響 沿用之前的計(jì)算模型及參數(shù),僅改變模型涂層厚度,將涂層厚度設(shè)為0.5、1.5、2、2.5 mm,分別計(jì)算,比較不同模型的最大位移值,并繪制涂層厚度增強(qiáng)效果曲線(xiàn),結(jié)果如表9及圖16所示。

    由表9和圖16可以看出,隨著涂層厚度的增大,罐體的抗爆效果雖不斷增強(qiáng),但曲線(xiàn)斜率逐漸平緩,意味著繼續(xù)增大涂層厚度,增益會(huì)越來(lái)越低。在工程應(yīng)用中,僅靠增大涂層厚度并不能解決一切問(wèn)題,反而會(huì)使成本大幅增加,根據(jù)實(shí)際情況選擇合適的涂層厚度才能最大化發(fā)揮涂層的作用。

    圖16 涂層厚度增強(qiáng)效果曲線(xiàn)Fig.16 Curve of coating thickness enhancement

    表9 不同涂層厚度罐體最大位移計(jì)算值對(duì)照Table 9 Comparison between maximum displacements of tank with different coating thickness

    3 結(jié)論

    針對(duì)聚脲涂層對(duì)鋼制罐體抗爆能力的影響進(jìn)行了試驗(yàn)研究與數(shù)值模擬,并將結(jié)果進(jìn)行對(duì)比分析,得到以下結(jié)論:

    1)在2 kg TNT炸藥不同距離的爆炸作用下,噴涂層罐體的測(cè)點(diǎn)最大位移與殘余位移均較無(wú)涂層罐體有明顯的降低,降幅約為30%,說(shuō)明聚脲涂層對(duì)于鋼罐體結(jié)構(gòu)強(qiáng)度的提升效果較為明顯。

    2)利用ALE流固耦合分析方法對(duì)鋼罐體及周?chē)諝庥蜻M(jìn)行計(jì)算,并將計(jì)算結(jié)果與試驗(yàn)值進(jìn)行比較,結(jié)果表明,ALE方法能夠較為準(zhǔn)確地反映罐體結(jié)構(gòu)響應(yīng),測(cè)點(diǎn)位移峰值與試驗(yàn)值較為接近,但該方法得到的殘余位移比試驗(yàn)值大。從模擬分析結(jié)果來(lái)看,噴涂聚脲涂層能夠有效地降低罐體測(cè)點(diǎn)位移與最大位移,降幅在20%~35%,與試驗(yàn)結(jié)論相符合。

    3)將聚脲涂層與鋼罐體采用不同連接形式進(jìn)行模擬計(jì)算,結(jié)果表明,TIED綁定接觸和TIEBREAK失效接觸兩種形式得到的結(jié)果基本一致,而共節(jié)點(diǎn)連接形式只有位移峰值與前兩者較為接近,其他結(jié)果均差距明顯。在缺少聚脲涂層與鋼板連接的試驗(yàn)數(shù)據(jù)時(shí),可以采用TIED綁定接觸形式,能夠滿(mǎn)足大多數(shù)情況,而共節(jié)點(diǎn)形式則不推薦采用。

    4)對(duì)3種常用鋼材制作的罐體進(jìn)行了模擬分析,結(jié)果表明,聚脲涂層對(duì)不同材料罐體的抗爆能力均有一定程度的提升,但隨著鋼材強(qiáng)度的提高,聚脲涂層帶來(lái)的強(qiáng)度增益有所降低。

    5)對(duì)不同涂層厚度的罐體進(jìn)行了模擬分析并比較其增強(qiáng)效果,結(jié)果表明,隨著聚脲涂層厚度的增加,罐體抗爆效果也不斷增強(qiáng),但增加厚度帶來(lái)的增益逐漸下降。

    猜你喜歡
    抗爆聚脲模擬計(jì)算
    復(fù)合夾芯纖維增強(qiáng)水泥板抗爆性能研究
    R1234ze PVTx熱物性模擬計(jì)算
    能源工程(2022年1期)2022-03-29 01:06:26
    石油化工控制室抗爆設(shè)計(jì)規(guī)范
    天津化工(2021年3期)2021-01-08 18:40:57
    淺談聚脲材料的性能分析及其在土木工程中的應(yīng)用
    水下爆炸載荷作用下聚脲材料對(duì)鋼結(jié)構(gòu)防護(hù)效果研究
    抗爆裝甲門(mén)的程序化設(shè)計(jì)
    山西建筑(2016年27期)2016-11-14 09:32:33
    擠出發(fā)泡片材褶皺分析及模擬計(jì)算
    SK單組分聚脲在混凝土缺陷處理中的應(yīng)用
    三維樣條小波單元構(gòu)造及其在地下箱型結(jié)構(gòu)抗爆數(shù)值模擬中的應(yīng)用
    乙二醇改性聚脲微膠囊相變材料的研制
    av在线老鸭窝| 欧美性猛交╳xxx乱大交人| 亚洲va在线va天堂va国产| 国产欧美日韩一区二区精品| 国产aⅴ精品一区二区三区波| 亚洲成av人片在线播放无| 国产视频一区二区在线看| 3wmmmm亚洲av在线观看| 舔av片在线| 日本欧美国产在线视频| 黄色欧美视频在线观看| 国产精品一区二区性色av| 99久久精品热视频| 亚洲精品久久国产高清桃花| 色av中文字幕| 亚洲色图av天堂| 久久国产乱子免费精品| 国产精品爽爽va在线观看网站| 一边摸一边抽搐一进一小说| 亚洲人成网站在线播| av在线亚洲专区| 成人毛片a级毛片在线播放| 久久精品夜夜夜夜夜久久蜜豆| 免费看美女性在线毛片视频| 最近的中文字幕免费完整| 久久久久久久亚洲中文字幕| 日本黄色视频三级网站网址| 色吧在线观看| 亚洲欧美清纯卡通| 日本一本二区三区精品| 精品不卡国产一区二区三区| 高清日韩中文字幕在线| 性色avwww在线观看| 久久精品国产亚洲av天美| 伊人久久精品亚洲午夜| 美女大奶头视频| 深爱激情五月婷婷| 免费看日本二区| 亚洲精华国产精华液的使用体验 | 亚洲国产高清在线一区二区三| 欧美日韩在线观看h| 有码 亚洲区| 高清毛片免费观看视频网站| 级片在线观看| 精品一区二区三区视频在线观看免费| 波多野结衣高清作品| 一级毛片我不卡| 直男gayav资源| 伦理电影大哥的女人| 亚洲精品一区av在线观看| 综合色av麻豆| 最新在线观看一区二区三区| 日本熟妇午夜| 日韩制服骚丝袜av| 在线观看午夜福利视频| 成人美女网站在线观看视频| 国产精品综合久久久久久久免费| 国产精品爽爽va在线观看网站| 51国产日韩欧美| 美女被艹到高潮喷水动态| 色哟哟哟哟哟哟| 黑人高潮一二区| 最近手机中文字幕大全| 99热精品在线国产| 午夜福利在线在线| 亚洲七黄色美女视频| 又爽又黄无遮挡网站| 伦精品一区二区三区| 国产精品福利在线免费观看| 少妇裸体淫交视频免费看高清| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添av毛片| 黄片wwwwww| 国产精品99久久久久久久久| 在线a可以看的网站| 亚洲成人中文字幕在线播放| 美女内射精品一级片tv| 成人美女网站在线观看视频| 久久久久久九九精品二区国产| 哪里可以看免费的av片| 麻豆精品久久久久久蜜桃| 中文资源天堂在线| 午夜a级毛片| 免费观看人在逋| 久久久国产成人精品二区| 精品少妇黑人巨大在线播放 | 永久网站在线| 狂野欧美激情性xxxx在线观看| 一本一本综合久久| 免费人成在线观看视频色| 国产爱豆传媒在线观看| 日本爱情动作片www.在线观看 | 亚洲av美国av| 亚洲国产色片| 亚洲,欧美,日韩| 亚洲成人中文字幕在线播放| 午夜福利视频1000在线观看| 国产在线精品亚洲第一网站| 国产精品久久久久久亚洲av鲁大| 一夜夜www| 校园春色视频在线观看| 日本黄色视频三级网站网址| 舔av片在线| 成年女人毛片免费观看观看9| 两性午夜刺激爽爽歪歪视频在线观看| www日本黄色视频网| 狠狠狠狠99中文字幕| 精品久久久久久成人av| 欧美区成人在线视频| 国产精品国产三级国产av玫瑰| 亚洲最大成人手机在线| 一级毛片aaaaaa免费看小| 此物有八面人人有两片| 一级毛片久久久久久久久女| 天堂影院成人在线观看| 人妻夜夜爽99麻豆av| 久久久久久国产a免费观看| 97超级碰碰碰精品色视频在线观看| 国产精品99久久久久久久久| 国产精品国产三级国产av玫瑰| 在线免费十八禁| 人人妻人人看人人澡| 午夜影院日韩av| av在线播放精品| 亚洲图色成人| 大型黄色视频在线免费观看| 搡老妇女老女人老熟妇| 色av中文字幕| 国产午夜福利久久久久久| a级毛片免费高清观看在线播放| 99热这里只有是精品在线观看| 免费看光身美女| 男女做爰动态图高潮gif福利片| 一个人观看的视频www高清免费观看| 久久久久精品国产欧美久久久| 国产一级毛片七仙女欲春2| 亚洲精品色激情综合| 麻豆成人午夜福利视频| 久久久久国产网址| 国产亚洲91精品色在线| 干丝袜人妻中文字幕| 婷婷精品国产亚洲av| 日韩欧美三级三区| 免费看日本二区| 小说图片视频综合网站| 噜噜噜噜噜久久久久久91| 老司机午夜福利在线观看视频| 国语自产精品视频在线第100页| 色吧在线观看| 中文字幕精品亚洲无线码一区| 欧美激情国产日韩精品一区| 久久婷婷人人爽人人干人人爱| 春色校园在线视频观看| 欧美日韩国产亚洲二区| 精品乱码久久久久久99久播| 欧美3d第一页| 三级国产精品欧美在线观看| 偷拍熟女少妇极品色| 成人二区视频| 亚洲成人久久性| 亚洲性夜色夜夜综合| 亚洲欧美精品综合久久99| 国产精品亚洲一级av第二区| 久久综合国产亚洲精品| 联通29元200g的流量卡| 精品国产三级普通话版| 99在线视频只有这里精品首页| 欧美日本视频| 国产精品av视频在线免费观看| 国内精品宾馆在线| 老熟妇乱子伦视频在线观看| 淫秽高清视频在线观看| 九九爱精品视频在线观看| 夜夜看夜夜爽夜夜摸| 欧美性猛交黑人性爽| 精品午夜福利视频在线观看一区| 成人亚洲欧美一区二区av| 国产一区二区在线av高清观看| 女人被狂操c到高潮| 欧美区成人在线视频| 欧美一区二区国产精品久久精品| 色尼玛亚洲综合影院| 午夜福利高清视频| av在线播放精品| 精品久久久噜噜| 人人妻,人人澡人人爽秒播| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻| 人人妻人人澡欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 联通29元200g的流量卡| 少妇丰满av| 精品一区二区三区视频在线| 国语自产精品视频在线第100页| 少妇的逼好多水| 成年女人毛片免费观看观看9| 久久精品91蜜桃| 欧美性猛交╳xxx乱大交人| 天堂影院成人在线观看| 精品久久久久久久久亚洲| 国产精华一区二区三区| 久久久久国产精品人妻aⅴ院| 国产男人的电影天堂91| 日本黄色片子视频| 免费在线观看影片大全网站| 99久久无色码亚洲精品果冻| 久久久久久久久中文| 国产高清激情床上av| 乱系列少妇在线播放| 看十八女毛片水多多多| 免费看美女性在线毛片视频| 日韩一区二区视频免费看| 精品久久久久久久久亚洲| 1024手机看黄色片| 一级黄色大片毛片| 欧美日韩乱码在线| 午夜日韩欧美国产| 免费大片18禁| 成人一区二区视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产美女午夜福利| av.在线天堂| 日韩制服骚丝袜av| 美女 人体艺术 gogo| 国产精品99久久久久久久久| 久久国产乱子免费精品| 99在线视频只有这里精品首页| 欧美激情久久久久久爽电影| 综合色av麻豆| 成人一区二区视频在线观看| 久久中文看片网| 欧美中文日本在线观看视频| 在线观看av片永久免费下载| 国产精品久久久久久av不卡| 久久久久久久久久成人| 亚洲国产欧美人成| 国产精品一区二区免费欧美| 精品人妻视频免费看| 男人和女人高潮做爰伦理| 欧美激情久久久久久爽电影| 亚洲性夜色夜夜综合| 久久久久精品国产欧美久久久| 蜜臀久久99精品久久宅男| 一级毛片我不卡| 国产免费男女视频| 亚洲精品国产成人久久av| 久久久久精品国产欧美久久久| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 99视频精品全部免费 在线| 欧美性感艳星| 国产精品人妻久久久久久| 草草在线视频免费看| 国产中年淑女户外野战色| 亚洲欧美成人精品一区二区| 床上黄色一级片| 在线播放国产精品三级| 一a级毛片在线观看| 少妇熟女欧美另类| 久久久久久九九精品二区国产| 亚洲高清免费不卡视频| 国产精品一区二区三区四区免费观看 | 国产精品人妻久久久影院| 成年女人毛片免费观看观看9| ponron亚洲| 日韩欧美精品免费久久| 黄色日韩在线| 岛国在线免费视频观看| 国产探花极品一区二区| 舔av片在线| 国产色婷婷99| 亚州av有码| 晚上一个人看的免费电影| 又粗又爽又猛毛片免费看| 一级毛片我不卡| 久久鲁丝午夜福利片| 91久久精品国产一区二区三区| 国产又黄又爽又无遮挡在线| 男人舔女人下体高潮全视频| 国产亚洲欧美98| 国产精品综合久久久久久久免费| 老司机影院成人| 亚洲av二区三区四区| 深爱激情五月婷婷| 91在线观看av| 久久久精品欧美日韩精品| 99热这里只有是精品50| 久久99热这里只有精品18| 蜜桃亚洲精品一区二区三区| 免费看a级黄色片| 欧美色视频一区免费| av天堂中文字幕网| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| av视频在线观看入口| 欧美高清成人免费视频www| 精品午夜福利在线看| 久久久久精品国产欧美久久久| 精华霜和精华液先用哪个| 三级毛片av免费| 日本-黄色视频高清免费观看| 一个人观看的视频www高清免费观看| 欧美日韩国产亚洲二区| 美女xxoo啪啪120秒动态图| 啦啦啦观看免费观看视频高清| 丰满的人妻完整版| 亚洲人成网站在线播| 国产真实伦视频高清在线观看| 狂野欧美白嫩少妇大欣赏| 麻豆久久精品国产亚洲av| 97热精品久久久久久| 成人高潮视频无遮挡免费网站| or卡值多少钱| 欧美日本视频| av在线天堂中文字幕| 日韩一区二区视频免费看| 日本精品一区二区三区蜜桃| 中文在线观看免费www的网站| 少妇人妻一区二区三区视频| 最新在线观看一区二区三区| 亚洲美女黄片视频| 最近2019中文字幕mv第一页| 欧美+亚洲+日韩+国产| 日韩亚洲欧美综合| 全区人妻精品视频| 给我免费播放毛片高清在线观看| 国产视频一区二区在线看| 亚洲五月天丁香| 亚洲国产精品合色在线| av视频在线观看入口| 欧美绝顶高潮抽搐喷水| av在线播放精品| 内射极品少妇av片p| 黄色视频,在线免费观看| 一级毛片aaaaaa免费看小| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 国产高清视频在线观看网站| 亚洲在线自拍视频| 亚洲av成人av| 亚洲欧美日韩高清在线视频| 少妇的逼好多水| 啦啦啦观看免费观看视频高清| 国产成年人精品一区二区| 国产av一区在线观看免费| 午夜a级毛片| 午夜福利成人在线免费观看| 日本 av在线| 99热这里只有精品一区| 成人国产麻豆网| 精品午夜福利视频在线观看一区| 97人妻精品一区二区三区麻豆| 亚洲成a人片在线一区二区| 日日啪夜夜撸| 最新在线观看一区二区三区| 美女免费视频网站| 成熟少妇高潮喷水视频| 久久精品国产自在天天线| 蜜桃亚洲精品一区二区三区| 少妇丰满av| 精品99又大又爽又粗少妇毛片| 搡老岳熟女国产| 欧美3d第一页| 国产不卡一卡二| 99久国产av精品| 亚洲成人中文字幕在线播放| 精品99又大又爽又粗少妇毛片| 色噜噜av男人的天堂激情| 欧美极品一区二区三区四区| 国产高清激情床上av| 国产久久久一区二区三区| 在线观看一区二区三区| 九九爱精品视频在线观看| 一区福利在线观看| 日韩欧美国产在线观看| 久久久久久久久久黄片| 日韩大尺度精品在线看网址| 男女边吃奶边做爰视频| 18+在线观看网站| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添小说| 变态另类丝袜制服| 亚洲成av人片在线播放无| 99九九线精品视频在线观看视频| 天天一区二区日本电影三级| 免费搜索国产男女视频| 天天一区二区日本电影三级| 国产男人的电影天堂91| 特大巨黑吊av在线直播| 国产伦在线观看视频一区| 伦精品一区二区三区| 香蕉av资源在线| 国产日本99.免费观看| 免费在线观看影片大全网站| 亚洲五月天丁香| 国产伦在线观看视频一区| 一区福利在线观看| 一级黄片播放器| 国产高清视频在线播放一区| 欧美3d第一页| 久久久国产成人免费| 成人二区视频| 亚洲性久久影院| 亚洲av美国av| 成熟少妇高潮喷水视频| 亚洲国产日韩欧美精品在线观看| 一区二区三区高清视频在线| 观看美女的网站| 日韩一区二区视频免费看| 色综合亚洲欧美另类图片| 国产日本99.免费观看| 欧美日韩一区二区视频在线观看视频在线 | 18禁在线无遮挡免费观看视频 | 超碰av人人做人人爽久久| 午夜精品一区二区三区免费看| 久久精品国产自在天天线| 久久精品夜色国产| 久久久精品欧美日韩精品| 欧美成人免费av一区二区三区| 精品国产三级普通话版| 尤物成人国产欧美一区二区三区| 亚洲av成人av| 国产亚洲精品av在线| 日本爱情动作片www.在线观看 | 国产淫片久久久久久久久| 91久久精品电影网| 三级毛片av免费| 亚洲成人精品中文字幕电影| 日本三级黄在线观看| 精品乱码久久久久久99久播| 联通29元200g的流量卡| 美女xxoo啪啪120秒动态图| 日韩,欧美,国产一区二区三区 | 国产精品女同一区二区软件| 尾随美女入室| 亚洲欧美日韩卡通动漫| 亚洲国产高清在线一区二区三| 亚洲国产色片| 亚洲自拍偷在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美一区二区三区在线观看| 淫妇啪啪啪对白视频| 久久久精品94久久精品| 国产成人freesex在线 | 亚洲精品在线观看二区| 老司机福利观看| 免费av毛片视频| 十八禁国产超污无遮挡网站| 久久精品影院6| 91在线观看av| 六月丁香七月| 一进一出抽搐gif免费好疼| 偷拍熟女少妇极品色| 国产免费一级a男人的天堂| 69av精品久久久久久| 国产午夜精品论理片| 91久久精品国产一区二区三区| 亚洲欧美日韩高清在线视频| 久久久久久久久久黄片| av黄色大香蕉| eeuss影院久久| 亚洲一区高清亚洲精品| 国产一级毛片七仙女欲春2| 俺也久久电影网| 成人亚洲欧美一区二区av| 免费人成视频x8x8入口观看| 国产精品乱码一区二三区的特点| 美女xxoo啪啪120秒动态图| 真人做人爱边吃奶动态| 亚洲精品色激情综合| 观看免费一级毛片| 丰满的人妻完整版| 欧美国产日韩亚洲一区| 3wmmmm亚洲av在线观看| 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 99热全是精品| 久久精品夜夜夜夜夜久久蜜豆| 国产在线精品亚洲第一网站| 国产精品精品国产色婷婷| 国产白丝娇喘喷水9色精品| 久久九九热精品免费| 婷婷精品国产亚洲av在线| 99热这里只有精品一区| 国产av不卡久久| 欧美人与善性xxx| 一边摸一边抽搐一进一小说| 亚洲精品亚洲一区二区| 老司机福利观看| 精品午夜福利在线看| 91狼人影院| 成人无遮挡网站| 大型黄色视频在线免费观看| 精品久久久久久成人av| 床上黄色一级片| 国产成人freesex在线 | 亚洲色图av天堂| 老司机影院成人| 国产午夜精品久久久久久一区二区三区 | 国产精品99久久久久久久久| 高清毛片免费观看视频网站| 日产精品乱码卡一卡2卡三| 亚洲国产精品久久男人天堂| 亚洲内射少妇av| 天天一区二区日本电影三级| 99久国产av精品国产电影| 人人妻人人澡欧美一区二区| 亚洲天堂国产精品一区在线| 亚洲无线在线观看| 免费av不卡在线播放| 亚洲国产精品成人久久小说 | а√天堂www在线а√下载| 成人特级黄色片久久久久久久| aaaaa片日本免费| 国产高清激情床上av| 可以在线观看毛片的网站| 午夜亚洲福利在线播放| 亚洲欧美精品自产自拍| 久久久国产成人精品二区| 精品乱码久久久久久99久播| 国产午夜精品论理片| 高清毛片免费观看视频网站| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全电影3| or卡值多少钱| 国产精品亚洲一级av第二区| 99热这里只有是精品在线观看| 一进一出好大好爽视频| 国产精品人妻久久久久久| 精品久久久久久久末码| 又粗又爽又猛毛片免费看| 内地一区二区视频在线| 色综合站精品国产| 在线观看66精品国产| 一级黄片播放器| 色综合亚洲欧美另类图片| 亚洲精华国产精华液的使用体验 | 中文字幕人妻熟人妻熟丝袜美| 97在线视频观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美xxxx性猛交bbbb| 久久婷婷人人爽人人干人人爱| 成人精品一区二区免费| 国产精品久久视频播放| 亚洲av中文av极速乱| 国产在视频线在精品| 听说在线观看完整版免费高清| 国产黄片美女视频| 午夜激情福利司机影院| 中文字幕精品亚洲无线码一区| 最新中文字幕久久久久| 麻豆一二三区av精品| 久久天躁狠狠躁夜夜2o2o| av天堂在线播放| 久久久a久久爽久久v久久| 国产精品女同一区二区软件| 国产极品精品免费视频能看的| 国产单亲对白刺激| 国产国拍精品亚洲av在线观看| 大香蕉久久网| 国产精品美女特级片免费视频播放器| 国产伦精品一区二区三区视频9| 亚洲欧美日韩高清在线视频| 成人午夜高清在线视频| 日韩一本色道免费dvd| 亚洲图色成人| 亚洲性夜色夜夜综合| 久久精品91蜜桃| 99久国产av精品| 男女边吃奶边做爰视频| 成人亚洲欧美一区二区av| 久久99热6这里只有精品| 国产精品人妻久久久影院| 亚洲国产精品合色在线| 国产麻豆成人av免费视频| 香蕉av资源在线| 精品少妇黑人巨大在线播放 | 一本久久中文字幕| 国产真实伦视频高清在线观看| 大又大粗又爽又黄少妇毛片口| 丝袜美腿在线中文| 国产乱人视频| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 久久热精品热| 尤物成人国产欧美一区二区三区| 成人漫画全彩无遮挡| 国产一级毛片七仙女欲春2| 精品99又大又爽又粗少妇毛片| 一级毛片电影观看 | 内地一区二区视频在线| 亚洲性夜色夜夜综合| 国产成人精品久久久久久| 亚洲av电影不卡..在线观看| 成人特级av手机在线观看| 插逼视频在线观看| 精品午夜福利视频在线观看一区| 色综合站精品国产| 欧美日韩一区二区视频在线观看视频在线 | 春色校园在线视频观看| eeuss影院久久| 亚洲av中文字字幕乱码综合| 一本久久中文字幕| 欧美日韩在线观看h| 国产探花在线观看一区二区| 18禁在线无遮挡免费观看视频 | 国产大屁股一区二区在线视频| 久久久久久大精品| 欧美日韩一区二区视频在线观看视频在线 | 中文字幕精品亚洲无线码一区| 久久久久久伊人网av| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 免费av毛片视频| 欧美日韩国产亚洲二区|