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

    焊接加固熱作用對(duì)工形截面壓彎鋼構(gòu)件承載性能的影響

    2016-12-07 06:48:49蔣立王元清戴國欣張?zhí)焐?/span>石永久
    關(guān)鍵詞:壓彎熱源鋼結(jié)構(gòu)

    蔣立,王元清,戴國欣,張?zhí)焐?,石永?/p>

    (1.重慶市設(shè)計(jì)院,重慶 400015;2.清華大學(xué) 土木工程安全與耐久教育部重點(diǎn)實(shí)驗(yàn)室;土木工程系,北京 100084;3.重慶大學(xué) 土木工程學(xué)院,重慶 400045)

    ?

    焊接加固熱作用對(duì)工形截面壓彎鋼構(gòu)件承載性能的影響

    蔣立1,3,王元清2,戴國欣3,張?zhí)焐?,石永久2

    (1.重慶市設(shè)計(jì)院,重慶 400015;2.清華大學(xué) 土木工程安全與耐久教育部重點(diǎn)實(shí)驗(yàn)室;土木工程系,北京 100084;3.重慶大學(xué) 土木工程學(xué)院,重慶 400045)

    鋼結(jié)構(gòu)加固;負(fù)載下;焊接熱作用;壓彎構(gòu)件;承載性能;有限元分析

    中國工程建設(shè)發(fā)展正經(jīng)歷從大規(guī)模新建向新建與加固改造并重過渡的階段,中國用于對(duì)既有建筑的加固改造費(fèi)用逐年增加,特別是工業(yè)建筑。面對(duì)中國工業(yè)化對(duì)城市化推動(dòng)力嚴(yán)重不足,淘汰落后產(chǎn)能及推進(jìn)技術(shù)改造已成為當(dāng)前政府工作的重要部署和企業(yè)生產(chǎn)力的發(fā)展重心,因而對(duì)已有鋼結(jié)構(gòu)在加固技術(shù)方面提出了更高的要求。負(fù)載下焊接加固技術(shù)是傳統(tǒng)鋼結(jié)構(gòu)加固方法中最具普遍運(yùn)用價(jià)值的一種,已在大量加固工程中被采用,但焊接熱影響問題的存在已成為完善負(fù)載下焊接加固技術(shù)急需解決的關(guān)鍵問題之一[1]。

    目前的負(fù)載下焊接加固鋼結(jié)構(gòu)研究幾乎全部集中于軸壓構(gòu)件和受彎構(gòu)件,僅李少甫等[2]采用數(shù)值積分方法得出長程殘余應(yīng)力下偏壓桿極限承載力的計(jì)算公式,并作為《鋼結(jié)構(gòu)加固技術(shù)規(guī)范》(CECS77:96)[3]的基礎(chǔ),而此前尚未見關(guān)于負(fù)載下焊接加固壓彎構(gòu)件試驗(yàn)研究的報(bào)道[4],結(jié)合試驗(yàn)的相關(guān)分析和討論仍留存空白。實(shí)際工程中純粹的軸壓柱是不存在的,所有鋼柱都應(yīng)當(dāng)作壓彎構(gòu)件來處理[5]。因而,負(fù)載下焊接加固壓彎構(gòu)件及其焊接熱影響問題的研究對(duì)于鋼結(jié)構(gòu)負(fù)載下加固技術(shù)體系的完善極其必要。

    1 有限元建模

    1.1 試驗(yàn)簡況

    為研究負(fù)載焊接加固的初始負(fù)載和焊接加固過程熱影響對(duì)鋼構(gòu)件受力性能的影響,在清華大學(xué)土木工程安全與耐久教育部重點(diǎn)實(shí)驗(yàn)室進(jìn)行的工字形壓彎構(gòu)件負(fù)載加固試驗(yàn)是有限元驗(yàn)證的基礎(chǔ)[4]。如圖1所示,被加固鋼柱均為Q345B鋼焊接工字形截面,采用翼緣外對(duì)稱貼焊鋼板加固方案,試件在柱底固接、柱頂面內(nèi)自由面外無平動(dòng)的約束條件下,柱頂面內(nèi)固定100 mm偏心距加載。

    圖1 試驗(yàn)裝置及試件加固后截面尺寸示意

    試件的初始負(fù)載、實(shí)測(cè)尺寸及初始幾何缺陷見表1,常溫下拉伸試驗(yàn)所得材料力學(xué)性能見表2。

    表1 試件初始負(fù)載、實(shí)測(cè)尺寸和幾何初始缺陷

    表2 材性試驗(yàn)結(jié)果

    1.2 有限元模型

    圖2 間接熱結(jié)構(gòu)耦合法分析流程圖

    1.2.2 焊接熱輸入 焊接熱源具有瞬時(shí)性(加熱速度極快,電弧焊時(shí)可達(dá)1 500 ℃/s以上)、局部集中性和運(yùn)動(dòng)性的特點(diǎn),在時(shí)間及空間域內(nèi)均易形成較大梯度的非均勻溫度場,對(duì)于焊件,這種局部不均勻熱輸入是產(chǎn)生焊接應(yīng)力與變形的決定因素,因而,焊接數(shù)值模擬的關(guān)鍵是焊接熱源模型的合理選取[9]。

    移動(dòng)高斯熱源模型被認(rèn)為是對(duì)目前常用焊接方法進(jìn)行數(shù)值模擬的較理想模型,常用來進(jìn)行較精細(xì)的焊接殘余應(yīng)力分析,但相當(dāng)耗時(shí)[10]。Variyar等[11]采用9節(jié)點(diǎn)的退化殼單元按高斯熱源加載模擬焊接殘余應(yīng)力,結(jié)果表明,使用高斯熱源和殼單元對(duì)薄板和厚板均適用,與實(shí)驗(yàn)結(jié)果吻合較好。牛犇等[12]采用三維實(shí)體單元建模在表面施加熱流密度的熱源模型,對(duì)1∶5縮尺構(gòu)件進(jìn)行熱分析獲得溫度場結(jié)果。對(duì)于實(shí)際加固構(gòu)件的全長角焊縫情況,受彎和軸壓構(gòu)件的焊接加固研究采用移動(dòng)點(diǎn)熱源的簡化模型[6-7],但熱輸入以1 500 ℃為控制目標(biāo)尚未能體現(xiàn)實(shí)際熱輸入的影響。筆者基于與移動(dòng)高斯熱源模型熱量輸入等效原則,采用ANSYS殼單元模擬焊縫建立了簡化的移動(dòng)分段點(diǎn)熱源模型,即按試驗(yàn)分段采用只在每道焊縫中心施加熱輸入功率的特殊串熱源模型。

    相關(guān)數(shù)值算例表明,采用基于與移動(dòng)高斯熱源熱量輸入等效原則的串熱源模型,可以在保持精度的同時(shí)提高計(jì)算效率高達(dá)200倍,適于對(duì)構(gòu)件的焊接模擬[13]。高斯熱源的熱流密度分布可描述為[9]

    q(r)=qmexp(-Kr2)

    (1)

    式中:q(r)和qm分別為熱流密度和加熱斑點(diǎn)中心最大比熱流,J/(m2·s);K為能量集中系數(shù),m-2。

    串熱源模型的簡化效率得到證實(shí),而另一方面,實(shí)際試驗(yàn)?zāi)芸刂频氖请娏骱碗妷?,需要借助有限元工具建立二者的銜接。通用有限元軟件ANSYS為瞬態(tài)熱分析提供了BF或BFE命令[14],可對(duì)通過節(jié)點(diǎn)或單元施加熱生成載荷HGEN(體荷載HGEN表示單位體積上的生熱速率W/m3)實(shí)現(xiàn)熱輸入。考慮在ANSYS建模時(shí)對(duì)每一道焊縫劃分為4個(gè)單元并作為一組,采用只在焊縫單元組中心節(jié)點(diǎn)用BF命令施加熱輸入功率的特殊串熱源(或點(diǎn)熱源)模型,仍基于熱量輸入等效原則,建立如下關(guān)系進(jìn)行銜接

    HGEN=P/V

    (2)

    P=ηUI?HGEN=ηUI/V

    (3)

    式中:P表示有效熱輸入功率;η表示功率有效系數(shù),為簡化討論取η=1.0;V表示所施加熱輸入的焊縫單元組的體積或節(jié)點(diǎn)的從屬體積。

    1.2.4 初始?xì)堄鄳?yīng)力 班慧勇等[17]匯總了普通鋼材及高強(qiáng)鋼材的殘余應(yīng)力試驗(yàn)數(shù)據(jù),提出了可以考慮不同強(qiáng)度等級(jí)鋼材、不同截面尺寸的焊接工字形截面殘余應(yīng)力統(tǒng)一分布模型,歐洲鋼結(jié)構(gòu)協(xié)會(huì)(ECCS)《鋼結(jié)構(gòu)穩(wěn)定設(shè)計(jì)手冊(cè)》[18]對(duì)單獨(dú)焰切割加固板的殘余應(yīng)力采用矩形分布模型。考慮到鋼結(jié)構(gòu)設(shè)計(jì)規(guī)范所用典型截面[19]尺寸與筆者所用尺寸差別較大,中國沒有相關(guān)的單獨(dú)焰切割加固板殘余應(yīng)力研究報(bào)道,故對(duì)工字形截面和加固板分別采用班慧勇模型和ECCS模型,得到的初始?xì)堄鄳?yīng)力分布分別如圖3(a)和圖4(a)所示,有限元考慮網(wǎng)格劃分后所取用的焊接殘余應(yīng)力分布簡化模型如圖3(b)和圖4(b)所示。

    圖3 焊接工字形截面的初始?xì)堄鄳?yīng)力分布

    圖4 加固板截面的初始?xì)堄鄳?yīng)力分布

    1.2.5 加固焊接順序及焊接模擬 試驗(yàn)中加固板先焊偏心受壓遠(yuǎn)側(cè),再焊受壓近側(cè)。圖5為鋼柱負(fù)載焊接加固順序,遠(yuǎn)側(cè)和近側(cè)均按A-B-C-…-H-I的區(qū)段順序,對(duì)長度為200 mm的A區(qū)段分每道焊長50 mm的8道焊;對(duì)區(qū)段長度均為350 mm的B~I(xiàn)在每區(qū)段分焊長為70 mm的10道焊。每道焊縫現(xiàn)場焊接約20 s,焊后停歇40~60 s。

    圖5 焊接加固施焊順序

    分析時(shí)為盡量模擬試驗(yàn)真實(shí)的焊接過程,點(diǎn)熱源作用的坐標(biāo)定位按試驗(yàn)的施焊次序隨時(shí)間變化,熱輸入功率則由試驗(yàn)實(shí)測(cè)有效功率確定。

    2 有限元分析結(jié)果

    2.1 焊接溫度場

    有限元熱分析的瞬態(tài)溫度場模擬結(jié)果見圖6,模擬的環(huán)境溫度取為試驗(yàn)實(shí)測(cè)室內(nèi)溫度28.6 ℃。焊接時(shí),焊點(diǎn)附近高溫區(qū)云圖形狀呈以焊接方向?yàn)殚L軸的近似紡錘形或橢圓形,這與焊接工藝參數(shù)相近的多層多道焊試驗(yàn)的首層焊溫度分布結(jié)果相似[20],對(duì)稱施焊的兩焊點(diǎn)溫度場相互疊加在加固板中間形成鞍部。每焊道焊接時(shí)高溫區(qū)溫度最高分別可達(dá)1 024.9 ℃(對(duì)A區(qū)段)和725.7 ℃(對(duì)B~I(xiàn)區(qū)段),但溫度梯度大,距離焊點(diǎn)越遠(yuǎn),溫度迅速降低,如圖6(a)和(b)所示。每道焊縫焊接操作間歇時(shí),最高溫度迅速降低到250 ℃以下,溫度場分布如圖6(c)所示,在加固板平面內(nèi)溫度場呈鉞形,在腹板平面內(nèi)呈半紡錘形略滯后于加固板溫度場,隨焊點(diǎn)移動(dòng)而溫度場分布動(dòng)態(tài)穩(wěn)定地整體向前推移。

    每道焊縫均采用相同熱功率輸入,以及相同的焊接20 s和焊后停歇60 s的時(shí)間,而每道焊縫長A區(qū)段比其他區(qū)段小,故A區(qū)段升溫更高。

    圖6 試件BCS3的溫度場云圖

    試驗(yàn)采用手持式紅外線測(cè)溫儀對(duì)試件BCS3進(jìn)行焊接及冷卻過程溫度測(cè)試,獲得了溫度場特征點(diǎn)變化數(shù)據(jù),其中,焊接受壓近側(cè)各區(qū)段焊接完畢對(duì)應(yīng)的區(qū)段加固板中心和腹板中心溫度實(shí)測(cè)值列于表3。有限元與試驗(yàn)結(jié)果的相應(yīng)對(duì)比可知,加固板處結(jié)果吻合良好。注意到,腹板中心有限元結(jié)果比試驗(yàn)結(jié)果略小,原因是加固板與翼緣之間實(shí)際存在一定程度的接觸,可以傳遞熱流,而有限元未能考慮。而加固板中心溫度有限元結(jié)果比試驗(yàn)結(jié)果在開始焊接時(shí)(A區(qū)段)略高,原因是柱底與剛性基座之間實(shí)際存在溫度傳遞,有限元模型熱分析未建立基座。

    表3 各區(qū)段焊接完畢時(shí)對(duì)應(yīng)區(qū)段加固板中心和腹板中心溫度的有限元結(jié)果與試驗(yàn)結(jié)果對(duì)比(BCS3)

    續(xù)表3

    2.2 極限變形狀態(tài)

    采用考慮熱影響的有限元模型得到3個(gè)負(fù)載下焊接加固壓彎構(gòu)件的相應(yīng)分析結(jié)果,發(fā)現(xiàn)鋼柱均發(fā)生同試驗(yàn)中破壞模態(tài)一致的空間的彎扭失穩(wěn),典型試件BCS2的極限變形狀態(tài)如圖7所示。

    圖7 典型試件BCS2的極限變形

    2.3 焊接位移時(shí)程曲線

    試驗(yàn)中對(duì)負(fù)載下焊接加固過程進(jìn)行了位移監(jiān)測(cè),得到焊接時(shí)程曲線如圖8所示,考慮熱影響的有限元分析得到的位移時(shí)程曲線同時(shí)作于圖上,對(duì)比可見,兩者在趨勢(shì)上吻合較好,有限元可以揭示與試驗(yàn)結(jié)果一致的規(guī)律。

    圖8 焊接位移時(shí)程曲線有限元結(jié)果與試驗(yàn)結(jié)果對(duì)比

    有限元的時(shí)程曲線呈局部和總體波動(dòng)變化。局部的波動(dòng)是由于每側(cè)焊接時(shí)各區(qū)段采用對(duì)稱交替的焊接道次順序;而“一峰一谷”的總體波動(dòng)變化,是由先焊接受壓近側(cè),再焊接受壓遠(yuǎn)側(cè)的焊接次序決定,每焊接一側(cè)時(shí),柱頂面內(nèi)位移先向與焊接側(cè)相反的方向變動(dòng),然后向與焊接側(cè)相同的方向變動(dòng)??梢?,焊接次序決定焊接變形的發(fā)展規(guī)律。注意到試驗(yàn)的位移時(shí)程曲線的峰谷處趨平且對(duì)稱交替焊接的局部波動(dòng)不明顯,是由于試驗(yàn)?zāi)Σ劣绊懀邢拊捎谒捎媚Σ聊P团c生死單元法之間的兼容性問題,未能反映持載焊接加固過程的摩擦影響。

    試驗(yàn)與有限元的焊接時(shí)程曲線相應(yīng)位移變化的特征值見表4。由表4可知,初始負(fù)載越大,焊接殘余變形越大;初始負(fù)載越大,焊接熱影響的峰谷波動(dòng)位移幅值與焊接殘余位移的比值越小。此外,由圖8也可看出,初始負(fù)載越大,焊接殘余變形將趨于超過焊接過程由焊接高溫?zé)釕?yīng)力引起的波峰幅值,提醒關(guān)注焊接過程的熱輸入。因此,焊接熱輸入大小和初始負(fù)載大小共同決定了持載焊接的位移變化范圍和焊接殘余變形的大小。

    表4 焊接加固過程的柱頂面內(nèi)位移變化特征值

    2.4 截面應(yīng)力應(yīng)變

    不考慮熱影響的有限元由于激活單元后僅是恢復(fù)先前被殺死單元的物理力學(xué)參數(shù)[8],不能模擬持載焊接過程的相互作用,因而,得不到截面焊接應(yīng)力應(yīng)變重分布,截面應(yīng)變與激活單元前一致。

    試驗(yàn)時(shí)獲得了負(fù)載下焊接加固試件腹板個(gè)別位置的應(yīng)變時(shí)程數(shù)據(jù),試件BCS3柱中截面和柱三分點(diǎn)截面腹板的考慮熱影響的有限元應(yīng)變分布與試驗(yàn)實(shí)測(cè)結(jié)果對(duì)比如圖9所示,包括焊接加固前后及極限承載力時(shí)的情況。圖9中有限元應(yīng)變結(jié)果為分析獲得的彈性應(yīng)變EPEL、塑性應(yīng)變EPPL和熱應(yīng)變EPTH之和。對(duì)比可知考慮熱影響的有限元得到與試驗(yàn)結(jié)果一致的規(guī)律:焊接加固后,腹板壓應(yīng)變整體相對(duì)加固前增大,即焊接加固過程引起截面應(yīng)力應(yīng)變重分布現(xiàn)象。此外,分析結(jié)果還表明,初始負(fù)載越大,應(yīng)力應(yīng)變重分布往偏心受壓近側(cè)發(fā)展越多。由于試驗(yàn)比未考慮摩擦的有限元極限承載力高,極限承載力對(duì)應(yīng)的應(yīng)變?cè)囼?yàn)值比有限元值偏大較多。

    圖9 柱中截面和柱三分點(diǎn)截面腹板有限元應(yīng)變分布與試驗(yàn)實(shí)測(cè)結(jié)果對(duì)比

    試驗(yàn)時(shí)由于焊接高溫使得貼于翼緣的應(yīng)變片燒損,未能得到有效數(shù)據(jù),故通過有限元分析來探討焊接加固過程的截面主要受力部分的應(yīng)力應(yīng)變重分布情況。考慮熱影響的有限元給出了受力最大的柱底截面的焊接加固前后柱底截面應(yīng)力云圖如圖10,提取受壓近側(cè)翼緣和加固板的單元應(yīng)力結(jié)果,分別作出應(yīng)力分布如圖11、12所示。

    圖10 焊接加固前后柱底截面(軸向)應(yīng)力云圖(單位:MPa)

    圖11 焊接加固前后柱底截面翼緣(軸向)單元應(yīng)力

    圖12 焊接加固前后柱底截面加固板(軸向)

    據(jù)圖可知,焊接前加固板無初始應(yīng)力,焊接后偏壓近側(cè)加固板應(yīng)力并不自相平衡而趨于出現(xiàn)合壓力,而偏壓近側(cè)翼緣的合壓力趨于減小,且初始負(fù)載越大,應(yīng)力重分布越明顯。注意到偏壓遠(yuǎn)側(cè)這種現(xiàn)象并不明顯。所以,截面應(yīng)力應(yīng)變重分布往偏心受力的方向發(fā)展,這與腹板情況的結(jié)論相一致。

    綜上可知,經(jīng)歷負(fù)載下焊接加固的試件,其焊接熱影響下的截面應(yīng)力變化情況非常復(fù)雜。由于逐漸形成的加固焊縫帶動(dòng)截面新舊部分同時(shí)受力,焊接殘余應(yīng)力形成過程中同時(shí)伴隨初始負(fù)載在新舊部分應(yīng)力重分布。此外,負(fù)載下焊接使得構(gòu)件產(chǎn)生殘余變形,由此造成的荷載效應(yīng)放大也將帶來新增應(yīng)力疊加效果。

    圖13 有限元與試驗(yàn)荷載位移曲線結(jié)果對(duì)比

    對(duì)于負(fù)載下加固構(gòu)件BCS1~BCS3,其考慮摩擦的有限元分析過程中由于過盈漸變的摩擦設(shè)置不適用于生死單元模型[8],故不能模擬出加載初期的剛度偏高現(xiàn)象;而混合接觸設(shè)置的摩擦影響主要體現(xiàn)在激活后的繼續(xù)加載階段,后期曲線剛度及趨勢(shì)與試驗(yàn)較好吻合??紤]熱影響和摩擦的BCS1在達(dá)承載力附近模擬出了受摩擦影響的承載力相對(duì)位移變化的不規(guī)則波動(dòng)上升現(xiàn)象,這與試驗(yàn)現(xiàn)象相吻合,可能是因?yàn)槟Σ潦锹窂较嚓P(guān)的影響因素,BCS1的初始幾何缺陷相比BCS2和BCS3較大造成柱頂偏向約束梁的趨勢(shì)更大而引起。

    2.6 邊緣屈服承載力

    CECS 77:96對(duì)接近軸壓受力的壓彎構(gòu)件采用同軸心受力的加固折減系數(shù)公式,可以考慮初始負(fù)載的影響,而對(duì)大部分偏心受力構(gòu)件的加固計(jì)算統(tǒng)一取等同受彎構(gòu)件的加固折減系數(shù)ηEM為0.9,尚未考慮初始負(fù)載的影響。新編國家標(biāo)準(zhǔn)《鋼結(jié)構(gòu)加固設(shè)計(jì)規(guī)范》征求意見稿按0.2、0.4、0.65三檔初始應(yīng)力比分級(jí)分別取定折減系數(shù),對(duì)壓彎構(gòu)件接近軸心受壓時(shí)仍取同軸心受壓構(gòu)件,取值在分級(jí)處保持原CECS 77:96公式水平,而對(duì)壓彎構(gòu)件偏心受力時(shí)所取用系數(shù)介于軸壓構(gòu)件和受彎構(gòu)件之間來實(shí)現(xiàn)過渡,均考慮初始負(fù)載的影響。

    現(xiàn)有規(guī)范對(duì)于平面外穩(wěn)定計(jì)算均采用的是基于穩(wěn)定屈曲理論的簡化直線公式,只是公式中套用軸壓構(gòu)件和受彎構(gòu)件的穩(wěn)定系數(shù)來進(jìn)行協(xié)調(diào)。此外,折減系數(shù)的來源也主要是平面內(nèi)承載力的數(shù)值計(jì)算結(jié)果,基于加固后全截面邊緣準(zhǔn)則的修正[21]。

    表5給出了分別按CECS 77:96規(guī)范和新編《鋼結(jié)構(gòu)加固設(shè)計(jì)規(guī)范》征求意見稿的設(shè)計(jì)承載力與有限元邊緣屈服時(shí)對(duì)應(yīng)承載力結(jié)果對(duì)比??梢娬髑笠庖姼逯械脑O(shè)計(jì)承載力水平大致是不考慮熱影響的有限元的邊緣屈服承載力的水平,而考慮熱影響的有限元的邊緣屈服承載力結(jié)果則過于保守。

    表5 不考慮摩擦的有限元邊緣屈服承載力結(jié)果對(duì)比

    2.7 極限承載力

    BCS1~BCS3試件在熱作用和摩擦兩種影響因素不同組合下分析得到的有限元極限承載力與試驗(yàn)結(jié)果對(duì)比如表6所示。結(jié)合圖13也可看出,僅就極限承載力而言,初始負(fù)載較大的BCS3考慮摩擦和熱影響的有限元結(jié)果與試驗(yàn)結(jié)果較吻合,初始負(fù)載較小的BCS1和BCS2考慮熱影響的承載力結(jié)果整體比不考慮熱影響的有限元結(jié)果低,也比試驗(yàn)結(jié)果更低。

    表6 有限元極限承載力結(jié)果與試驗(yàn)結(jié)果對(duì)比

    雖然采用考慮熱影響的有限元能夠模擬出焊接時(shí)程曲線,獲得試驗(yàn)類似的平臺(tái)段,但在平臺(tái)段長度和承載力及其對(duì)應(yīng)變形方面與試驗(yàn)相比仍有一些差異,產(chǎn)生差異的原因主要有:1)試驗(yàn)實(shí)際存在的摩擦對(duì)受力及變形的影響復(fù)雜,有限元難于準(zhǔn)確考慮;2)實(shí)際焊接過程的復(fù)雜多變,有限元難于完全真實(shí)模擬;3)有限元采用的高溫下材性與試驗(yàn)鋼材真實(shí)的高溫本構(gòu)存在一定差異;4)有限元同時(shí)考慮摩擦和熱影響分析使得問題的非線性程度更大,增加了模擬的復(fù)雜性和收斂難度。

    排除有限元中邊界摩擦的復(fù)雜影響,僅考察熱影響對(duì)承載力的降低程度,對(duì)比不考慮熱影響和考慮熱影響的有限元結(jié)果(見表7),發(fā)現(xiàn)初始負(fù)載越大,考慮熱影響的結(jié)果相對(duì)不考慮熱影響的結(jié)果有折減越大的趨勢(shì),從BCS1到BCS3依次降低約7%、8%、9%,即初始負(fù)載越大,越應(yīng)在承載力折減中將熱影響考慮進(jìn)去。再進(jìn)而考察初始負(fù)載對(duì)承載力的總體折減程度,發(fā)現(xiàn)BCS1和BCS2不考慮熱影響的有限元結(jié)果以及BCS3考慮熱影響的有限元結(jié)果的折減比值分別為0.99、0.99和0.86,這與試驗(yàn)的相應(yīng)折減比值0.97、1.02和0.89較為吻合。即在初始負(fù)載較小時(shí)采用不考慮熱影響的有限元,在初始負(fù)載較大時(shí)采用考慮熱影響的有限元計(jì)算結(jié)果可以獲得較好的結(jié)果。

    表7 排除摩擦的有限元承載力結(jié)果與試驗(yàn)結(jié)果對(duì)比

    至此,對(duì)有限元承載力結(jié)果的分析尚未涉及初始?xì)堄鄳?yīng)力,下面將給出考慮初始?xì)堄鄳?yīng)力和不考慮初始?xì)堄鄳?yīng)力對(duì)極限承載力結(jié)果對(duì)比。如表8所示,除不考慮熱影響的BCS2之外,考慮初始?xì)堄鄳?yīng)力后承載力均只略降低,可見不論考慮熱影響與否,初始?xì)堄鄳?yīng)力對(duì)極限承載力都幾乎沒有影響。

    然而與規(guī)范設(shè)計(jì)承載力結(jié)果對(duì)比可見,雖然不考慮摩擦的有限元結(jié)果比試驗(yàn)結(jié)果偏低,但相當(dāng)于邊緣屈服承載力水平的規(guī)范計(jì)算值偏低更多,因此,規(guī)范設(shè)計(jì)方法有可以提升的空間,有待進(jìn)一步擴(kuò)大參數(shù)化分析進(jìn)行驗(yàn)證。

    表8 不考慮摩擦的有限元承載力結(jié)果與試驗(yàn)結(jié)果對(duì)比

    3 結(jié) 論

    1)建立了能夠考慮實(shí)際熱輸入的方法,簡化串熱源模型的熱分析,獲得了負(fù)載焊接溫度場分布,與試驗(yàn)結(jié)果吻合較好,表明有限元熱分析的有效性。

    3)采用考慮熱影響的有限元方法模擬負(fù)載下焊接加固壓彎構(gòu)件的受力性能具有一定可行性和總體安全性,為后續(xù)開展參數(shù)化研究提供了方法和依據(jù)。

    4)初始負(fù)載較低時(shí)可采用不考慮熱影響的有限元,初始負(fù)載較大時(shí)應(yīng)采用考慮熱影響的有限元。

    5)考慮熱影響的有限元分析結(jié)果表明,初始?xì)堄鄳?yīng)力幾乎不影響負(fù)載下焊接加固壓彎構(gòu)件的整體穩(wěn)定極限承載力。

    6)新編《鋼結(jié)構(gòu)加固設(shè)計(jì)規(guī)范》征求意見稿的設(shè)計(jì)承載力水平大致是不考慮熱影響的有限元的邊緣屈服承載力的水平,比考慮熱影響的有限元的極限承載力更低,比較保守,設(shè)計(jì)方法有提升的空間。

    [1] 蔣立, 王元清, 戴國欣, 等. 焊接熱作用對(duì)鋼構(gòu)件負(fù)載下加固的影響分析[C]//第十四屆全國現(xiàn)代結(jié)構(gòu)工程學(xué)術(shù)研討會(huì)論文集, 天津, 2014: 755-763.

    JIANG L, WANG Y Q, DAI G X, et al. Analysis of welding heat effect on steel members strengthened by weld-ing while under load [C]// Symposium of 14th National Modern Structure Conference, Tianjin, 2014: 755-763. (in Chinese)

    [2] 李少甫, 張納新. 長程殘余應(yīng)力下偏壓桿的極限承載力[J]. 建筑結(jié)構(gòu), 1994, 6: 13-16.

    LI S F, ZHANG N X. Ultimate load-carrying capacity of beam-column under long range residual stress [J]. Building Structure, 1994, 6: 13-16. (in Chinese)

    [3] 中國工程建設(shè)標(biāo)準(zhǔn)化協(xié)會(huì). 鋼結(jié)構(gòu)加固技術(shù)規(guī)范: CECS 77:96 [S]. 北京: 中國計(jì)劃出版社, 2005.

    China Association for Engineering Construction Standardization. Technical code for strengthening steel structures: CECS 77:96 [S]. Beijing: China Planning Press, 2005. (in Chinese)

    [4] 王元清, 蔣立, 戴國欣, 等. 負(fù)載下鋼結(jié)構(gòu)工字形壓彎構(gòu)件焊接加固試驗(yàn)[J]. 哈爾濱工業(yè)大學(xué)學(xué)報(bào), 2016,48(6):30-37.

    WANG Y Q, JIANG L, DAI G X, et al. Experimental study on I section steel beam-columns strengthened by welding while under load [J]. Journal of Harbin Institute of Technology, 2016,48(6):30-37.(in Chinese)

    [5] Chen W F, Atsuta T. Theory of beam-columns, Volume 1: In-plane behavior and design [M]. United States: J Ross Publishing, 2007.

    [6] 王元清, 祝瑞祥, 戴國欣, 等. 負(fù)載下焊接加固受彎工形鋼梁的受力特性分析[J]. 建筑科學(xué)與工程學(xué)報(bào), 2013, 30(4): 112-120.

    WANG Y Q, ZHU R X, DAI G X, et al. Analysis on load-carrying behavior of I section steel beams strengthened with weld while under load by ANSYS [J]. Journal of Architecture and Civil Engineering, 2013, 30(4): 112-120. (in Chinese)

    [7] 王元清, 祝瑞祥, 戴國欣, 等. 工形鋼柱負(fù)載下焊接加固的受力特性[J]. 沈陽建筑大學(xué)學(xué)報(bào) (自然科學(xué)版), 2014, 30(1): 25-33.

    WANG Y Q, ZHU R X, DAI G X, et al. Analysis on load-carrying behavior of I section steel columns reinforced by welding with initial compressive load [J]. Journal of Shenyang Jianzhu University (Natural Science), 2014, 30(1): 25-33. (in Chinese)

    [8] 蔣立, 王元清, 戴國欣, 等. 負(fù)載下工形截面壓彎鋼構(gòu)件焊接加固的承載性能分析[J]. 天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版), 2015(Sup1): 60-66.

    JIANG L, WANG Y Q, DAI G X, et al. Analysis on load-carrying behavior of I section steel beam-columns strengthened by welding while under load [J]. Journal of Tianjin University (Science and Technology), 2015(Sup1): 60-66. (in Chinese)

    [9] 張文鉞. 焊接冶金學(xué): 基本原理[M]. 北京: 機(jī)械工業(yè)出版社, 1993.

    [10] 莫春立, 錢百年, 國旭明, 等. 焊接熱源計(jì)算模式的研究進(jìn)展[J]. 焊接學(xué)報(bào), 2001, 22(3): 93-96.

    MO C L, QIAN B N, GUO X M, et al. The development of models about welding heat source’s calculation [J]. Transactions of the China Welding Institution, 2001, 22(3): 93-96. (in Chinese)

    [11] Variyar S M, Prasad N S. Estimation of residual stresses in weldments using a nine-noded degenerated shell element [J]. Journal of Materials Processing Technology, 1999, 91(1): 150-160.

    [12] 牛犇, 劉曉珂, 陳志華, 等. 負(fù)載下焊接加固過程對(duì)鋼梁及鋼支撐受力性能的影響研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào), 2014, 35(7): 87-95.

    NIU B, LIU X K, CHEN Z H, et al. Study on influence of welding reinforcement on mechanical behavior of steel beam and braces in load status [J]. Journal of Building Structures, 2014, 35(7): 87-95. (in Chinese)

    [13] 蔡志鵬, 趙海燕, 吳甦, 等. 串熱源模型及其在焊接數(shù)值模擬中的應(yīng)用[J]. 機(jī)械工程學(xué)報(bào), 2001, 37(4): 25-28.

    CAI Z P, ZHAO H Y, WU S, et al. Model of string heat source in welding numerical simulations [J]. Chinese Journal of Mechanical Engineering, 2001, 37(4): 25-28. (in Chinese)

    [14] ANSYS Theory Reference Release 13.0 [M]. SAS IP, Inc., 2010.

    [15] 李國強(qiáng), 韓林海, 婁國彪, 等. 鋼結(jié)構(gòu)及鋼-混凝土組合結(jié)構(gòu)抗火設(shè)計(jì)[M]. 北京: 中國建筑工業(yè)出版社, 2006.

    [16] European Committee for Standardization. Eurocode 3: Design of steel structures-Part 1-2: general rules-structural fire design: EN 1993-1-2 [S]. 2005.

    [17] 班慧勇, 施剛, 石永久. 高強(qiáng)鋼焊接構(gòu)件工字形橫截面殘余應(yīng)力試驗(yàn)及統(tǒng)一分布模型研究[J]. 工程力學(xué), 2012, 31(8): 83-91.

    BAN H Y, SHI G, SHI Y J. Experimental and unified model investigations on residual stress within high strength steel welded I-sections [J]. Engineering Mechanics, 2012, 31(8): 83-91. (in Chinese)

    [18] European Convention for Constructional Steelworks. Manual on stability of steel structures [S]. 2nd ed. Bruxelles: ECCS Publ., 1976.

    [19] 李開禧, 饒曉峰. 單向偏心鋼壓桿彎扭屈曲臨界力[J]. 重慶建筑工程學(xué)院學(xué)報(bào), 1983, 4: 1-21.

    LI K X, RAO X F. Torsional-flexural buckling force of one-way eccentric steel strut [J]. Journal of Chongqing Institute of Architecture and Engineering, 1983, 4: 1-21. (in Chinese)

    [20] 李慧娟, 黃振華, 張京燾. 厚板多層多道焊的數(shù)值模擬分析[J]. 宇航材料工藝, 2007, 37(5): 30-33.

    LI H J, HUANG Z H, ZHANG J T. Numerical simulation of multilayer welding of thick plates [J]. Aerospace Materials & Technology, 2007, 37(5): 30-33. (in Chinese)

    [21] 蔣立, 王元清, 戴國欣, 等. 鋼結(jié)構(gòu)受壓柱負(fù)載下焊接加固設(shè)計(jì)方法比較分析[J]. 工業(yè)建筑, 2015, 2: 159-164.

    JIANG L, WANG Y Q, DAI G X, et al. Comparative study on design methods for steel columns strengthened by welding while under axial and eccentric load in different codes [J]. Industrial Construction, 2015, 2: 159-164. (in Chinese)

    (編輯 胡英奎)

    Analysis of welding heat effect on load-carrying behavior of I section steel beam-columns strengthened by welding

    Jiang Li1,3,Wang Yuanqing2,Dai Guoxin3,Zhang Tianshen2,Shi Yongjiu2

    (1. Chongqing Architectural Design Institute, Chongqing 400015, P. R. China; 2. Key Laboratory of Civil Engineering Safety and Durability of China Education Ministry; Department of Civil Engineering, Tsinghua University, Beijing 100084, P. R. China; 2. School of Civil Engineering, Chongqing University, Chongqing 400045, P. R. China)

    In order to investigate the influence of reinforcement welding heat effect and different initial load on load-carrying behavior of I section steel beam-columns, three specimens strengthened by welding while under different load were simulated. Based on the thermal-structural coupling method considering heat effect, the analysis improved thermal input of heat source model and also considered initial geometric imperfection, initial residual stress and friction. Time history of displacement during welding, web stress-strain redistribution and load-displacement relationship were discussed, and meanwhile the effectiveness of the finite element analyses(FEA) were verified by comparing the FEA results and test results. Furthermore, welding temperature field, stress-strain redistribution between flange and strengthening plate and flange margin yield capacity, which could not be measured by test, were acquired. And by comparison of bearing capacity results and code calculation results, the present design methods were discussed. The results showed that welding strengthening procedure decided the development of welding residual deformation, while thermal input and initial load effected the range of displacement change during welding under load and magnitude of welding residual deformation. Higher initial load propelled the development of stress-strain redistribution toward the eccentric loading direction, thus leading to lower bearing capacity. However, initial residual stress had little effect on ultimate capacity. FEA method considering heat effect showed certain practicability and overall security, and code design method remained space for improvement still.

    strengthening steel structure; under load; welding heat effect; beam-column; load-carrying behavior; finite element analysis

    2016-02-18

    《鋼結(jié)構(gòu)加固設(shè)計(jì)規(guī)范》國家標(biāo)準(zhǔn)管理組科研專項(xiàng)課題(No.2013-1)

    蔣立(1989-),男,主要從事鋼結(jié)構(gòu)研究, (E-mail)johnlee@cqu.edu.cn。

    王元清(通信作者),男,教授,博士生導(dǎo)師,(E-mail)wang-yq@mail.tsinghua.edu.cn。

    Foundation item:Specialized Funding Research of National Standards Management Group onCodeforDesignofStrengtheningSteelStructure(No.2013-1)

    TU391

    A

    1674-4764(2016)04-0096-12

    10.11835/j.issn.1674-4764.2016.04.015

    Received:2016-02-18

    Author brief:Jiang Li (1989-), main research interest:steel structure, (E-mail) johnlee@cqu.edu.cn.

    Wang Yuanqing(corresponding author), professor, doctorial supervisor, (E-mail)wang-yq@mail.tsinghua.edu.cn.

    猜你喜歡
    壓彎熱源鋼結(jié)構(gòu)
    頂推滑移法在鋼結(jié)構(gòu)橋梁施工中的應(yīng)用
    鋼結(jié)構(gòu)防火設(shè)計(jì)規(guī)范及要點(diǎn)探討
    首件鋼結(jié)構(gòu)構(gòu)件確認(rèn)機(jī)制在鋼結(jié)構(gòu)施工中的應(yīng)用
    山東冶金(2022年1期)2022-04-19 13:40:48
    一種“Z”形壓彎的級(jí)進(jìn)模設(shè)計(jì)
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    椰子樹為什么斜向海邊成長?
    6 45 度壓彎摩托車
    決策探索(2019年19期)2019-10-21 09:16:38
    基于啟發(fā)式動(dòng)態(tài)規(guī)劃的冷熱源優(yōu)化控制
    電子制作(2017年19期)2017-02-02 07:08:31
    中部槽激光-MAG復(fù)合熱源打底焊焊接工藝研究
    焊接(2015年8期)2015-07-18 10:59:13
    多類型熱源點(diǎn)共存下的區(qū)域熱力統(tǒng)籌供應(yīng)探討
    免费搜索国产男女视频| 国产一区二区在线观看日韩| 能在线免费观看的黄片| 日韩欧美一区二区三区在线观看| 婷婷亚洲欧美| 国产成年人精品一区二区| 亚洲成人免费电影在线观看| 成人精品一区二区免费| 夜夜夜夜夜久久久久| 亚洲无线在线观看| 日韩中字成人| 久久久久久久久中文| 亚洲美女搞黄在线观看 | 三级男女做爰猛烈吃奶摸视频| 最近最新免费中文字幕在线| 欧美潮喷喷水| 韩国av在线不卡| 日韩av在线大香蕉| 美女免费视频网站| 日韩精品中文字幕看吧| 亚洲欧美日韩高清在线视频| 久久久久性生活片| 日本免费a在线| 97人妻精品一区二区三区麻豆| 亚洲在线观看片| 男插女下体视频免费在线播放| 国产一区二区亚洲精品在线观看| 久久国产乱子免费精品| 真人做人爱边吃奶动态| bbb黄色大片| 国产亚洲91精品色在线| 美女大奶头视频| 看免费成人av毛片| 观看免费一级毛片| 在线看三级毛片| 国产黄a三级三级三级人| 亚洲自拍偷在线| 身体一侧抽搐| 精品99又大又爽又粗少妇毛片 | 国内精品一区二区在线观看| 超碰av人人做人人爽久久| 中文字幕av成人在线电影| 99久国产av精品| 亚洲av第一区精品v没综合| 深夜a级毛片| 国产毛片a区久久久久| 乱系列少妇在线播放| 18禁在线播放成人免费| 国产不卡一卡二| 12—13女人毛片做爰片一| 熟妇人妻久久中文字幕3abv| 亚洲精品国产成人久久av| 免费不卡的大黄色大毛片视频在线观看 | 亚洲成av人片在线播放无| 国产黄a三级三级三级人| 久久久久久久久久久丰满 | 久久久久久久久中文| 搡老熟女国产l中国老女人| 成人综合一区亚洲| 日韩一区二区视频免费看| 国产亚洲精品综合一区在线观看| 精品久久久久久久人妻蜜臀av| 精品久久久久久久人妻蜜臀av| 人人妻人人看人人澡| 国产精品乱码一区二三区的特点| 欧美人与善性xxx| 人妻制服诱惑在线中文字幕| 欧美日韩国产亚洲二区| 欧美精品啪啪一区二区三区| 亚洲av五月六月丁香网| 一级黄色大片毛片| 日韩中文字幕欧美一区二区| 久久热精品热| 国内少妇人妻偷人精品xxx网站| 男女啪啪激烈高潮av片| 久久久久久九九精品二区国产| 男女啪啪激烈高潮av片| 日日啪夜夜撸| 91av网一区二区| 丰满的人妻完整版| 精品午夜福利在线看| 亚洲精品国产成人久久av| 深爱激情五月婷婷| 日韩欧美免费精品| 岛国在线免费视频观看| 日本 欧美在线| 久9热在线精品视频| 午夜激情福利司机影院| 无人区码免费观看不卡| aaaaa片日本免费| 久久精品国产99精品国产亚洲性色| 人妻久久中文字幕网| 欧美国产日韩亚洲一区| 最近最新中文字幕大全电影3| 国产精品电影一区二区三区| 久久久久免费精品人妻一区二区| 亚洲国产精品合色在线| 日本-黄色视频高清免费观看| 一区二区三区高清视频在线| 欧美日本亚洲视频在线播放| 日韩欧美一区二区三区在线观看| 国产成人福利小说| 国产伦人伦偷精品视频| 精品久久久久久久末码| 国产淫片久久久久久久久| x7x7x7水蜜桃| 欧美成人免费av一区二区三区| 乱码一卡2卡4卡精品| 啦啦啦韩国在线观看视频| 联通29元200g的流量卡| 天天一区二区日本电影三级| 亚洲人成网站在线播放欧美日韩| 久久热精品热| 人人妻,人人澡人人爽秒播| 亚洲自拍偷在线| 久久精品夜夜夜夜夜久久蜜豆| 高清日韩中文字幕在线| 国产精品日韩av在线免费观看| 超碰av人人做人人爽久久| 99久久无色码亚洲精品果冻| 男插女下体视频免费在线播放| 精品一区二区三区视频在线| 一级av片app| 两人在一起打扑克的视频| 国内精品宾馆在线| 在线免费十八禁| 亚洲精品久久国产高清桃花| 国产精品女同一区二区软件 | 日日摸夜夜添夜夜添av毛片 | 网址你懂的国产日韩在线| 亚洲欧美日韩东京热| av在线蜜桃| 婷婷精品国产亚洲av在线| 久久婷婷人人爽人人干人人爱| 九九久久精品国产亚洲av麻豆| 成人永久免费在线观看视频| 白带黄色成豆腐渣| 精品一区二区三区人妻视频| 国产精品综合久久久久久久免费| av国产免费在线观看| 97超视频在线观看视频| 国产高清视频在线观看网站| 国产精品爽爽va在线观看网站| 久久久久久久精品吃奶| 亚洲欧美清纯卡通| x7x7x7水蜜桃| 成人美女网站在线观看视频| 黄色一级大片看看| 日韩中文字幕欧美一区二区| 嫩草影视91久久| 亚洲精品日韩av片在线观看| 亚洲美女视频黄频| 亚洲av成人av| 露出奶头的视频| 日本免费a在线| 日韩亚洲欧美综合| 久久精品国产亚洲网站| 亚洲久久久久久中文字幕| 最近最新中文字幕大全电影3| 少妇的逼水好多| 色5月婷婷丁香| 欧美日韩综合久久久久久 | 女生性感内裤真人,穿戴方法视频| 中文字幕熟女人妻在线| 久久午夜福利片| 欧美丝袜亚洲另类 | 精品一区二区免费观看| 淫秽高清视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 简卡轻食公司| 人妻久久中文字幕网| 久久久精品大字幕| 91在线观看av| 99久久九九国产精品国产免费| 我的老师免费观看完整版| 美女大奶头视频| 久久精品国产清高在天天线| 日韩强制内射视频| 99riav亚洲国产免费| 熟女人妻精品中文字幕| 成人亚洲精品av一区二区| 啦啦啦观看免费观看视频高清| 波多野结衣巨乳人妻| 国产探花极品一区二区| 久久久久国产精品人妻aⅴ院| h日本视频在线播放| 12—13女人毛片做爰片一| 大型黄色视频在线免费观看| ponron亚洲| 亚洲专区中文字幕在线| 神马国产精品三级电影在线观看| 欧美日本亚洲视频在线播放| 精品国产三级普通话版| 日本精品一区二区三区蜜桃| 国产精品精品国产色婷婷| 久久精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 校园人妻丝袜中文字幕| 97热精品久久久久久| 少妇人妻精品综合一区二区 | 国产欧美日韩精品一区二区| 国产精品人妻久久久久久| 色综合色国产| 欧美区成人在线视频| 色在线成人网| 国产男人的电影天堂91| 99九九线精品视频在线观看视频| 男女视频在线观看网站免费| 麻豆国产av国片精品| 欧美高清性xxxxhd video| 美女被艹到高潮喷水动态| 成人国产综合亚洲| 国产精品综合久久久久久久免费| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区亚洲精品在线观看| 亚洲五月天丁香| 非洲黑人性xxxx精品又粗又长| 在线播放国产精品三级| 少妇的逼水好多| h日本视频在线播放| 欧美又色又爽又黄视频| 久久国产乱子免费精品| 日韩精品有码人妻一区| 免费高清视频大片| 亚洲内射少妇av| 亚洲最大成人av| 日韩av在线大香蕉| 日本撒尿小便嘘嘘汇集6| 人妻久久中文字幕网| 中文字幕免费在线视频6| 亚洲在线观看片| 在线a可以看的网站| 国产一区二区三区视频了| 又粗又爽又猛毛片免费看| 91麻豆精品激情在线观看国产| 美女xxoo啪啪120秒动态图| 欧美成人a在线观看| 干丝袜人妻中文字幕| 亚洲人成网站高清观看| 在线国产一区二区在线| 国产亚洲精品久久久久久毛片| 在线观看av片永久免费下载| 黄片wwwwww| 国国产精品蜜臀av免费| www.www免费av| 中出人妻视频一区二区| 亚洲成av人片在线播放无| 国产精品一区二区性色av| 搡女人真爽免费视频火全软件 | 不卡一级毛片| 精品国产三级普通话版| 国产黄色小视频在线观看| 少妇人妻精品综合一区二区 | 成人国产一区最新在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产在视频线在精品| 国产精品不卡视频一区二区| 1000部很黄的大片| 免费av毛片视频| 桃红色精品国产亚洲av| 成年免费大片在线观看| 精品国内亚洲2022精品成人| 不卡一级毛片| 成人三级黄色视频| 日本免费a在线| 国产亚洲精品久久久com| 91麻豆av在线| 亚洲无线观看免费| 99久久久亚洲精品蜜臀av| 亚洲美女视频黄频| 夜夜爽天天搞| 别揉我奶头 嗯啊视频| 日韩国内少妇激情av| 色噜噜av男人的天堂激情| 一个人观看的视频www高清免费观看| 人妻久久中文字幕网| 亚洲av五月六月丁香网| 草草在线视频免费看| 哪里可以看免费的av片| 99热6这里只有精品| 亚洲欧美日韩东京热| 男女做爰动态图高潮gif福利片| 亚洲av成人av| 久99久视频精品免费| 亚洲在线观看片| 国产精品久久久久久久电影| 国产精品,欧美在线| 午夜老司机福利剧场| 国产成人a区在线观看| 国产一级毛片七仙女欲春2| 免费电影在线观看免费观看| videossex国产| 亚洲美女搞黄在线观看 | 18+在线观看网站| 91久久精品电影网| 欧美日韩中文字幕国产精品一区二区三区| 禁无遮挡网站| 精品久久久噜噜| 99久久久亚洲精品蜜臀av| 偷拍熟女少妇极品色| 美女高潮的动态| 搡老熟女国产l中国老女人| 婷婷色综合大香蕉| 亚洲精品影视一区二区三区av| 国产av一区在线观看免费| 永久网站在线| 少妇的逼好多水| 午夜激情欧美在线| 亚洲,欧美,日韩| 一a级毛片在线观看| 欧美xxxx黑人xx丫x性爽| 国产高清三级在线| 中文字幕av成人在线电影| 我要看日韩黄色一级片| 久久草成人影院| 又紧又爽又黄一区二区| 日韩精品有码人妻一区| 久久午夜亚洲精品久久| 中文字幕高清在线视频| 3wmmmm亚洲av在线观看| 久久午夜亚洲精品久久| 在线国产一区二区在线| 日本 av在线| 人人妻人人澡欧美一区二区| 他把我摸到了高潮在线观看| 一进一出抽搐动态| 亚洲av成人精品一区久久| 99国产极品粉嫩在线观看| 久久6这里有精品| 国产色爽女视频免费观看| 色哟哟哟哟哟哟| or卡值多少钱| 国内久久婷婷六月综合欲色啪| 成人国产一区最新在线观看| 国产黄色小视频在线观看| 最近最新中文字幕大全电影3| av中文乱码字幕在线| 国内毛片毛片毛片毛片毛片| 夜夜夜夜夜久久久久| 春色校园在线视频观看| 黄色丝袜av网址大全| 国产精品久久久久久久电影| 亚洲狠狠婷婷综合久久图片| 日本熟妇午夜| 久久久久久久久大av| 欧美色欧美亚洲另类二区| 精品人妻熟女av久视频| 窝窝影院91人妻| 婷婷精品国产亚洲av在线| 一级黄色大片毛片| 在线观看一区二区三区| 免费观看的影片在线观看| 性欧美人与动物交配| 久久久久久久亚洲中文字幕| 精品久久久久久,| 最近在线观看免费完整版| 日韩精品青青久久久久久| 午夜免费成人在线视频| 日韩欧美三级三区| 深爱激情五月婷婷| 简卡轻食公司| 国产高清三级在线| 亚洲真实伦在线观看| 亚洲真实伦在线观看| 久久午夜福利片| 18禁裸乳无遮挡免费网站照片| 亚洲av电影不卡..在线观看| 国内精品久久久久精免费| 欧洲精品卡2卡3卡4卡5卡区| 亚洲 国产 在线| 观看美女的网站| 日韩国内少妇激情av| 色综合婷婷激情| 国产三级中文精品| 亚洲精品在线观看二区| 嫩草影院新地址| 亚洲一区高清亚洲精品| 99国产精品一区二区蜜桃av| 久久人人精品亚洲av| 国产精品98久久久久久宅男小说| 亚洲成av人片在线播放无| 又紧又爽又黄一区二区| 成人毛片a级毛片在线播放| 亚洲专区国产一区二区| 我要看日韩黄色一级片| 日日干狠狠操夜夜爽| 成人鲁丝片一二三区免费| 欧美另类亚洲清纯唯美| 国产一区二区三区在线臀色熟女| 日本欧美国产在线视频| 精品久久久久久久久久免费视频| 亚洲va日本ⅴa欧美va伊人久久| 中国美白少妇内射xxxbb| 永久网站在线| 2021天堂中文幕一二区在线观| 亚洲第一区二区三区不卡| 极品教师在线视频| 国产精品亚洲一级av第二区| 国产一区二区三区视频了| 夜夜爽天天搞| 麻豆一二三区av精品| 俄罗斯特黄特色一大片| 日韩欧美精品v在线| 毛片一级片免费看久久久久 | 精品久久久久久成人av| 免费观看在线日韩| 欧美三级亚洲精品| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 久久国产精品人妻蜜桃| 亚洲av五月六月丁香网| 精品久久久久久久久亚洲 | 欧美xxxx性猛交bbbb| 欧美不卡视频在线免费观看| 免费av不卡在线播放| 精品一区二区免费观看| 啦啦啦啦在线视频资源| 99热这里只有是精品在线观看| 欧美高清成人免费视频www| 国产综合懂色| 国产午夜福利久久久久久| 精品久久久久久久久av| 精品久久久久久久末码| 精品99又大又爽又粗少妇毛片 | av在线老鸭窝| 久久久久免费精品人妻一区二区| 欧美一区二区精品小视频在线| a级毛片a级免费在线| 国产精品永久免费网站| 国产白丝娇喘喷水9色精品| av国产免费在线观看| 亚洲国产色片| 亚洲av中文字字幕乱码综合| 欧美日韩精品成人综合77777| 中国美女看黄片| 欧美人与善性xxx| 精品一区二区免费观看| 国产中年淑女户外野战色| 精品久久久久久久久久久久久| 波多野结衣巨乳人妻| 婷婷六月久久综合丁香| 欧美日韩中文字幕国产精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 国产av麻豆久久久久久久| 亚洲人成网站在线播放欧美日韩| 久久精品国产清高在天天线| 色播亚洲综合网| 少妇的逼水好多| 国产亚洲精品久久久久久毛片| 长腿黑丝高跟| 亚洲欧美激情综合另类| 国产老妇女一区| 久久精品国产清高在天天线| 日韩 亚洲 欧美在线| av福利片在线观看| 日韩欧美三级三区| 欧美另类亚洲清纯唯美| 国产爱豆传媒在线观看| 国产伦一二天堂av在线观看| 国产精品乱码一区二三区的特点| 久久久久久伊人网av| 精品一区二区三区视频在线| 成人美女网站在线观看视频| 欧美三级亚洲精品| 免费在线观看影片大全网站| 日本五十路高清| 噜噜噜噜噜久久久久久91| 亚洲五月天丁香| 十八禁国产超污无遮挡网站| 夜夜看夜夜爽夜夜摸| 欧美一区二区国产精品久久精品| 人妻少妇偷人精品九色| 三级男女做爰猛烈吃奶摸视频| 亚洲av熟女| 又黄又爽又刺激的免费视频.| 亚洲精品日韩av片在线观看| 欧美丝袜亚洲另类 | 成年免费大片在线观看| 女同久久另类99精品国产91| 免费搜索国产男女视频| 欧美最黄视频在线播放免费| 搡老岳熟女国产| 可以在线观看的亚洲视频| 中文资源天堂在线| 日韩人妻高清精品专区| 深爱激情五月婷婷| 最近在线观看免费完整版| 成人欧美大片| 国产精品98久久久久久宅男小说| 免费不卡的大黄色大毛片视频在线观看 | 18禁裸乳无遮挡免费网站照片| 精品久久久久久成人av| 成年女人看的毛片在线观看| 一进一出好大好爽视频| 国模一区二区三区四区视频| 成人永久免费在线观看视频| 日韩精品有码人妻一区| 久久久国产成人免费| 亚洲精华国产精华液的使用体验 | 91麻豆精品激情在线观看国产| 最近最新免费中文字幕在线| 中文字幕免费在线视频6| 熟妇人妻久久中文字幕3abv| 精品不卡国产一区二区三区| 国产欧美日韩精品亚洲av| 国产精华一区二区三区| 欧美绝顶高潮抽搐喷水| 免费电影在线观看免费观看| av女优亚洲男人天堂| 亚洲在线自拍视频| 国产精品1区2区在线观看.| 制服丝袜大香蕉在线| 久久人妻av系列| 国产91精品成人一区二区三区| 亚洲av熟女| 欧美最黄视频在线播放免费| 自拍偷自拍亚洲精品老妇| 久久九九热精品免费| 日本黄色视频三级网站网址| 国内毛片毛片毛片毛片毛片| 国产综合懂色| 一级a爱片免费观看的视频| 亚洲av熟女| 老师上课跳d突然被开到最大视频| 国产一级毛片七仙女欲春2| av在线天堂中文字幕| 免费在线观看日本一区| 日韩欧美三级三区| 一级av片app| 亚洲成人免费电影在线观看| 狂野欧美激情性xxxx在线观看| av在线老鸭窝| 久久欧美精品欧美久久欧美| 两个人视频免费观看高清| 精品国内亚洲2022精品成人| 日日撸夜夜添| 亚洲国产精品久久男人天堂| 国产精品久久久久久亚洲av鲁大| 嫩草影院新地址| 欧美一级a爱片免费观看看| 国产高清激情床上av| av女优亚洲男人天堂| 天天一区二区日本电影三级| 91av网一区二区| 亚洲五月天丁香| 国产aⅴ精品一区二区三区波| 国产精品不卡视频一区二区| 亚洲精品成人久久久久久| 麻豆av噜噜一区二区三区| 国产视频一区二区在线看| 91久久精品国产一区二区三区| 日韩欧美在线二视频| 国产主播在线观看一区二区| 18禁黄网站禁片午夜丰满| 国产熟女欧美一区二区| www.www免费av| 中文字幕久久专区| 成人三级黄色视频| 国产成人aa在线观看| 国产高潮美女av| 亚洲国产欧美人成| 亚洲人成网站在线播| 精品一区二区三区人妻视频| 一区二区三区四区激情视频 | 老女人水多毛片| 两个人的视频大全免费| 久久亚洲真实| 婷婷精品国产亚洲av| 久久精品综合一区二区三区| 日韩欧美国产在线观看| 97超视频在线观看视频| 午夜福利成人在线免费观看| 国内精品美女久久久久久| 国语自产精品视频在线第100页| 国内精品美女久久久久久| 一区二区三区高清视频在线| 免费观看人在逋| 内地一区二区视频在线| 精品一区二区三区人妻视频| 露出奶头的视频| 国产亚洲精品综合一区在线观看| 在线免费观看的www视频| 亚洲avbb在线观看| 99久久中文字幕三级久久日本| 亚洲av电影不卡..在线观看| 亚洲自偷自拍三级| av在线老鸭窝| 国国产精品蜜臀av免费| 成人精品一区二区免费| 国产成人福利小说| 国产亚洲精品av在线| 国产精品爽爽va在线观看网站| 久久精品国产亚洲av涩爱 | 身体一侧抽搐| 久久精品国产亚洲av香蕉五月| 我要看日韩黄色一级片| 久久久久久久亚洲中文字幕| 一级a爱片免费观看的视频| 国内精品久久久久久久电影| 成人av一区二区三区在线看| 人妻少妇偷人精品九色| 成人高潮视频无遮挡免费网站| 国产一区二区三区在线臀色熟女| 欧美日韩国产亚洲二区| 少妇的逼好多水| 91在线观看av| 日本精品一区二区三区蜜桃| 国内精品久久久久久久电影| 国产 一区精品| 亚洲成a人片在线一区二区| 国产不卡一卡二| 国产视频内射| 欧美bdsm另类| 久久精品国产99精品国产亚洲性色| 久久精品国产亚洲av香蕉五月|