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

    厚鋼板對接焊接三維有限元數(shù)值模擬與分析

    2014-08-08 20:26:18郭彥林陳航袁星
    關(guān)鍵詞:板件層數(shù)計(jì)算結(jié)果

    郭彥林+陳航+袁星

    建筑科學(xué)與工程學(xué)報(bào)2014年文章編號:16732049(2014)01009008

    收稿日期:20131116

    基金項(xiàng)目:清華大學(xué)自主科研計(jì)劃杰出貢獻(xiàn)專項(xiàng)項(xiàng)目(2012Z10134)

    作者簡介:郭彥林(1958),男,陜西富平人,教授,博士研究生導(dǎo)師,工學(xué)博士

    摘要:采用有限元數(shù)值模擬方法,建立考慮熱固耦合效應(yīng)的鋼板對接焊接的三維有限元模型,運(yùn)用生死單元技術(shù)模擬焊接全過程,從焊接殘余變形和殘余應(yīng)力的角度,驗(yàn)證了分析模型的有效性,并分析了板件邊界約束、焊接層數(shù)與道數(shù)、板件厚度、焊接順序以及施工工藝對計(jì)算結(jié)果的影響。探討了簡化計(jì)算的可能性與正確性,并考慮了保溫和設(shè)置支撐等施工工藝對焊接施工的影響。研究成果可為設(shè)計(jì)與施工提供一定的參考。

    關(guān)鍵詞:厚鋼板;對接焊接;熱結(jié)構(gòu)耦合;殘余應(yīng)力;殘余變形;有限元;數(shù)值模擬

    中圖分類號:TU758文獻(xiàn)標(biāo)志碼:A

    3D Finite Element Numerical Simulation and Analysis on Butt

    Welding of Thick Steel PlateGUO Yanlin, CHEN Hang, YUAN Xing

    (Department of Civil Engineering, Tsinghua University, Beijing 100084, China)Abstract: Finite element numerical simulation was adopted to erect 3D finite element model, considering heatsolid coupling effect. Meanwhile, livedead element technology was applied to simulate the whole welding process. The effectiveness of the model was verified in the aspects of residual deformation and residual stress. The influences of plate boundary constraint, welding layer and pass, plate thickness, welding sequence and construction process on the behavior of steel plates on calculation results were analyzed. The possibility and correctness of simplified calculation were discussed. The influences of construction techniques such as heat preservation and support setting on welding were considered. The conclusions can provide some references to welding design and construction.

    Key words: thick steel plate; butt welding; thermalstructural coupling; residual stress; residual deformation; finite element; numerical simulation

    0引言

    目前,隨著高層結(jié)構(gòu)的發(fā)展,鋼板剪力墻、巨形鋼結(jié)構(gòu)柱和復(fù)雜鋼結(jié)構(gòu)節(jié)點(diǎn)的使用越來越多。在鋼結(jié)構(gòu)施工過程中,焊接產(chǎn)生的殘余應(yīng)力和殘余變形如果得不到有效的控制,會對結(jié)構(gòu)的受力和外觀產(chǎn)生很大的影響。

    焊接過程的物理現(xiàn)象非常復(fù)雜,是一個涉及高溫電弧物理、傳熱、冶金和力學(xué)的復(fù)雜過程[1]。焊接產(chǎn)生的殘余應(yīng)力與殘余變形之間存在相互的協(xié)調(diào),如殘余變形得到了釋放,則殘余應(yīng)力變小;反之,殘余變形被約束,殘余應(yīng)力就變大。目前獲取焊接殘余應(yīng)力和殘余變形的方法主要有試驗(yàn)測量法和數(shù)值模擬法。對于殘余變形的測量比較方便,而測量殘余應(yīng)力主要有應(yīng)力松弛法(包括鉆孔法、切條法等)、X射線法和裂紋法等。然而這些方法僅停留在實(shí)驗(yàn)室階段,用于小尺寸構(gòu)件的科學(xué)研究,無法適用于大型復(fù)雜構(gòu)件焊接殘余應(yīng)力的測量[2],因此對于鋼結(jié)構(gòu)焊接殘余應(yīng)力和殘余變形的研究多采用有限元數(shù)值模擬方法。

    對于簡單構(gòu)件焊接的有限元數(shù)值模擬方法,各國許多學(xué)者已做了一定的研究,但是大部分成果都集中在機(jī)械工程領(lǐng)域,主要涉及焊接溫度場、殘余應(yīng)力場的模擬,并沒有將焊接問題擴(kuò)展到土木工程領(lǐng)域,考慮約束度、焊接順序、施工工藝等因素對焊接殘余應(yīng)力與殘余變形的影響,對施工實(shí)踐幾乎沒有指導(dǎo)作用。焊接過程中溫度梯度大,導(dǎo)致材料的嚴(yán)重非線性,進(jìn)而使得求解必須采用瞬態(tài)分析并且求解過程收斂困難,各國學(xué)者的研究大多采用的是簡化模型,如二維平面假設(shè)或軸對稱單元[39],這導(dǎo)致實(shí)際問題與假設(shè)存在出入,計(jì)算結(jié)果誤差很大。三維有限元焊接數(shù)值模擬的文獻(xiàn)相對較少,國外學(xué)者對單道焊和多道焊的分析結(jié)果也主要集中在焊接溫度場和殘余應(yīng)力分布方面[1011],缺乏參數(shù)變化,對結(jié)構(gòu)分析意義不大?,F(xiàn)有的研究基本都集中在小尺寸鋼板對接焊接的模擬上,不涉及大尺寸和大厚度鋼板對接焊接的模擬,無法與工程實(shí)際接軌。

    本文中筆者應(yīng)用通用有限元計(jì)算軟件ANSYS,采用三維實(shí)體單元對大尺寸、大厚度鋼板和焊縫進(jìn)行建模,運(yùn)用生死單元技術(shù)模擬焊接全過程,研究了不同的約束度、焊接層數(shù)與道數(shù)、焊接順序和施工工藝對焊接殘余變形和殘余應(yīng)力大小和分布的影響,為設(shè)計(jì)和施工提供參考。

    1有限元模型及參數(shù)選取

    1.1熱結(jié)構(gòu)耦合關(guān)系

    焊接是熱固耦合問題,因此結(jié)構(gòu)的力學(xué)平衡方程為[1]

    σij,j=0(1)

    式中:σij,j為熱應(yīng)力。

    式(1)中包括了熱應(yīng)力的影響,其中σij包括了結(jié)構(gòu)受外力產(chǎn)生的應(yīng)力和熱應(yīng)力。結(jié)構(gòu)內(nèi)部應(yīng)力、應(yīng)變間的本構(gòu)方程為

    dσij=Dijkl(dεkl-dεpkl-dεckl-dεTkl)(2)

    式中:Dijkl為彈性本構(gòu)張拉系數(shù);dεkl,dεpkl,dεckl,dεTkl分別為總應(yīng)變、塑性應(yīng)變、蠕變應(yīng)變和熱應(yīng)變。

    熱應(yīng)變εTij與溫度和材料熱膨脹系數(shù)有關(guān),其表達(dá)式為

    εTij=αij(T-Tr)δij(3)

    δij=1i=j

    0i≠j(4)

    式中:αij為材料熱膨脹系數(shù);T為溫度;Tr為參考溫度;δij為δ算子。

    在通用有限元軟件ANSYS中,有2種考慮熱結(jié)構(gòu)耦合場的單元。一種是直接耦合單元,該單元除了具有普通單元的位移自由度之外,同時具有溫度自由度,在受到局部溫度作用后首先按照熱傳導(dǎo)理論計(jì)算單元中的溫度場,得到溫度場之后將溫度作為外荷載施加在單元上,再結(jié)合單元受力和位移約束得到單元內(nèi)力分布。整個計(jì)算過程中,結(jié)構(gòu)的內(nèi)力、位移以及溫度分布是同時得到的。這種類型的代表單元是實(shí)體單元Solid5,然而Solid5單元為彈性單元,不能考慮塑性[12]。另一種是間接耦合單元,即首先采用計(jì)算溫度場的單元計(jì)算結(jié)構(gòu)在局部溫度作用下的溫度場分布,得到整個結(jié)構(gòu)的溫度場之后,將溫度場單元轉(zhuǎn)化為相應(yīng)的結(jié)構(gòu)場單元,再進(jìn)行普通的內(nèi)力計(jì)算。這種類型的代表單元是溫度場單元Solid70與結(jié)構(gòu)場單元Solid185,其結(jié)構(gòu)場單元Solid185可以考慮塑性。

    endprint

    采用間接耦合單元,即只進(jìn)行溫度場與應(yīng)力、應(yīng)變場的單向耦合,只考慮溫度對應(yīng)力應(yīng)變關(guān)系的影響,而忽略應(yīng)力、應(yīng)變場對溫度場的影響。實(shí)踐證明[12],采用單向耦合的方法是合理的,本文中采用間接耦合單元進(jìn)行鋼板對接焊接過程的模擬和分析。

    1.2材料參數(shù)

    選擇2塊尺寸為1 m×1 m×0.03 m的板進(jìn)行對接焊接,焊縫為V形坡口,焊縫斷面尺寸如圖1所示,選用的鋼材為Q345普通結(jié)構(gòu)鋼,材料密度為7 850 kg·m-1,熱脹系數(shù)為1.48×10-5 m·℃-1,泊松比為0.29,熱傳導(dǎo)系數(shù)為34 W·(m·℃)-1,比熱容為420 J·(kg·℃)-1[13]。假設(shè)材料為理想彈塑性模型,Q345鋼材的力學(xué)性能見表1。

    圖1焊縫斷面尺寸

    Fig.1Crosssection Dimension of

    Welded Seam1.3焊接熱源

    ANSYS軟件中模擬焊接熱輸入的模型主要有高斯熱源模型和內(nèi)生熱率熱源模型。當(dāng)板件厚度不是特別大時,高斯熱源的熱供應(yīng)量是足夠的,而當(dāng)板件厚度比較大時,高斯熱源模擬焊接就與實(shí)際情況不一致,此時應(yīng)該采用內(nèi)生熱率熱源模型。本文中的研究對象是不同厚度的鋼板,采用內(nèi)生熱率熱源模型進(jìn)行模擬,通過施加內(nèi)生熱模擬電弧對焊縫的加熱作用[14]。

    內(nèi)生熱率熱源模型的內(nèi)部熱生成以熱生成強(qiáng)度

    表1Q345鋼材的力學(xué)性能

    Tab.1Mechanical Behaviors of Steel Q345溫度/℃205008001 2001 500初始彈性模量E0/GPa21017550101屈服強(qiáng)度fy/MPa34518030105屈服后彈性模量E/GPa00000來表示,即將有效的焊接熱輸入量換算成每道焊縫單元在單位體積、單位時間內(nèi)熱生成強(qiáng)度的供給熱能,在所有焊接熱能都施加到焊縫上后,內(nèi)部熱停止生成。熱生成強(qiáng)度q的計(jì)算公式為

    q=ηUIAweldvdt(4)

    式中:U為電壓;I為電流;η為電弧熱效率,η值在0.7~0.9范圍內(nèi),本文中取0.8;Aweld為焊縫的橫截面積;v為焊接速度;dt為每個荷載步的時間步長。

    焊接參數(shù)見表2。

    表2焊接參數(shù)

    Tab.2Welding Parameters焊道電流/A電壓/V焊接速度/(mm·s-1)首道160205中間道180225末道1602051.4邊界條件

    在進(jìn)行溫度場計(jì)算時,需要施加熱力學(xué)的邊界條件,包括恒溫邊界條件和熱對流邊界條件。恒溫邊界條件模擬周邊的外環(huán)境溫度;熱對流邊界條件模擬結(jié)構(gòu)與空氣的接觸面上的熱交換,假設(shè)對流換熱系數(shù)為常量13.9 W·(m2·K)-1。

    在進(jìn)行結(jié)構(gòu)場計(jì)算時,需要施加力學(xué)的邊界條件。本文中選擇了3個邊界條件分別進(jìn)行焊接計(jì)算和分析,即鋼板短邊固接、短邊簡支和短邊自由,長邊均完全自由,如圖2所示。

    圖2板件邊界約束條件(單位:mm)

    Fig.2Boundary Constraint Conditions of Plate (Unit:mm)自由板的約束條件為選擇模型中3個節(jié)點(diǎn)共6個自由度進(jìn)行約束。首先約束圖2(c)中節(jié)點(diǎn)1的3個平動自由度,形成一個球鉸,即假設(shè)自由板在該點(diǎn)處不發(fā)生平動位移;其次選擇節(jié)點(diǎn)2和節(jié)點(diǎn)3約束平面外的位移(z方向),即假設(shè)自由板在自由變形的過程中,節(jié)點(diǎn)1,2,3始終在同一平面上,即3個不共線的點(diǎn)可以確定一個平面;最后約束節(jié)點(diǎn)2的x方向的平動位移,保證整個板不會繞著節(jié)點(diǎn)1做剛體轉(zhuǎn)動。經(jīng)有限元計(jì)算驗(yàn)證,在這樣的約束條件下,3個節(jié)點(diǎn)的約束反力均為0,即表明沒有外力產(chǎn)生,板件發(fā)生自由變形,約束條件成立。

    1.5焊接過程的模擬方法

    采用生死單元技術(shù),將整條焊縫分為很多小段,每段分為若干道。首先將焊縫單元全部殺死,當(dāng)焊接到某一小段時,將此處的焊縫單元激活,施加內(nèi)生熱模擬焊接時的局部高溫區(qū)域,然后刪除該內(nèi)生熱,激活下一單元并開始下一步的焊接,以此循環(huán)完成全部焊接工作。每道焊接完成后,根據(jù)焊接過程中溫度的控制要求,控制冷卻時間。

    焊接初始的預(yù)熱和焊接完成后的保溫通過改變環(huán)境溫度來實(shí)現(xiàn),其中預(yù)熱溫度為60 ℃,焊接完成后以250 ℃的恒溫保溫90 min。2計(jì)算結(jié)果與分析

    2.1焊接層數(shù)與道數(shù)

    在精確的焊接模擬中,焊接層數(shù)與道數(shù)的劃分應(yīng)該與實(shí)際施工工藝相同。然而,對于大型鋼結(jié)構(gòu),如大型鋼板剪力墻、巨型柱的焊接施工,如果按照實(shí)際的焊接層數(shù)與道數(shù)進(jìn)行網(wǎng)格劃分和計(jì)算,必然會產(chǎn)生巨大的計(jì)算量,耗費(fèi)非常多的時間。因此需要找到簡化算法(如減少焊接層數(shù)或道數(shù)),分析簡化算法對焊接模擬結(jié)果的具體影響,評判采用該簡化算法的可能性和正確性。

    對本文中30 mm厚的鋼板,選擇的焊接層數(shù)分別為1,2,5,焊接層數(shù)與道數(shù)的劃分如圖3所示。圖4為自由板面外變形示意。表3,4中分別給出了自由板在道數(shù)不同時V形焊縫最大殘余變形和最大殘余應(yīng)力的計(jì)算結(jié)果,其中殘余應(yīng)力未注明方向的均為Von Mises應(yīng)力。

    圖3焊接層數(shù)與道數(shù)的劃分

    Fig.3Divisions of Welding Layer and Pass圖4自由板面外變形示意

    Fig.4Graphs of Outofplane Deformation for Free Plate由表3可以看出,焊接層數(shù)為1層和2層時的面外變形與焊接層數(shù)為5層時的面外變形方向相

    表3自由板V形焊縫最大殘余變形計(jì)算結(jié)果

    Tab.3Calculation Results of Maximum Residual

    Deformation of V Welded Seam for Free Platemm層數(shù)125道數(shù)每層1道每層1道每層1道底層1道,

    其他層各2道面內(nèi)x

    方向收縮-0.992-1.164-2.453-2.571面外變形-4.022-3.4677.8237.329表4自由板V形焊縫最大殘余應(yīng)力計(jì)算結(jié)果

    Tab.4Calculation Results of Maximum Residual

    Stress of V Welded Seam for Free PlateMPa層數(shù)125道數(shù)每層1道每層1道每層1道底層1道,

    其他層各2道保溫完成284248249249最終冷卻294293293293反,這說明劃分不同的層數(shù)對焊接模擬計(jì)算結(jié)果的正確性具有重大影響。當(dāng)劃分層數(shù)為1層時,計(jì)算耗時短,但是假設(shè)整條焊縫是一次成型,其收縮過程為整個焊縫截面同時收縮。圖5為1層1道焊收縮示意。由圖5可以看出:對于V形焊縫,開口處的寬度大,冷卻收縮量大;底部寬度小,冷卻收縮量小,導(dǎo)致同時收縮時沿z軸的負(fù)向發(fā)生面外變形。

    在實(shí)際焊接中,V形焊縫的殘余變形是沿圖5中z方向的正向,即有限元模型劃分層數(shù)為1層和2層時計(jì)算的結(jié)果與實(shí)際不符,而劃分層數(shù)為5層時計(jì)算的結(jié)果從定性上來說是正確的。因此對于厚鋼板多層焊接來說,有限元模型的焊接層數(shù)不能取得太少,否則會導(dǎo)致錯誤的計(jì)算結(jié)果。

    圖51層1道焊收縮示意

    Fig.5Graph of Cooling Shrinkage of

    One Layer with One Pass對比5層1道焊與5層2道焊的計(jì)算結(jié)果可知,分道后的面外殘余變形減小,最大殘余應(yīng)力保持不變,這說明分道使得焊接熱量輸入得到分散,同一層的累積收縮量相比于不分道時減小。應(yīng)用5層1道焊模型的計(jì)算結(jié)果比多道焊模型的計(jì)算結(jié)果偏于保守,是可行的。

    endprint

    2.2板件邊界的約束

    焊縫冷卻收縮變形的過程中,板件的變形會受到周圍構(gòu)件的約束,進(jìn)而在結(jié)構(gòu)內(nèi)部產(chǎn)生殘余應(yīng)力。約束程度不同,殘余應(yīng)力和殘余變形的大小也不同。表5,6中分別給出了焊接層數(shù)為5層且每層焊接1道情況下,短邊固接板、短邊簡支板和自由板的對接焊接最大殘余變形和最大殘余應(yīng)力。

    表5不同約束下V形焊縫最大殘余變形計(jì)算結(jié)果

    Tab.5Calculation Results of Maximum Residual

    Deformation of V Welded Seam Under

    Different Constraintsmm邊界條件短邊固接板短邊簡支板自由板面內(nèi)x方向收縮-1.173-2.450-2.453面外變形4.1177.7317.823表6不同約束下V形焊縫最大殘余應(yīng)力計(jì)算結(jié)果

    Tab.6Calculation Results of Maximum Residual Stress of

    V Welded Seam Under Different ConstraintsMPa邊界條件短邊固接板短邊簡支板自由板保溫完成345304249最終冷卻345345293在不同的約束條件下,焊接殘余應(yīng)力和殘余變形的大小不同。短邊固接板的邊界約束最強(qiáng),面內(nèi)和面外的變形受約束,變形量最小,但殘余應(yīng)力最大。短邊簡支板與自由板的殘余變形量相近,大于短邊固接板的變形量,但殘余應(yīng)力較小。這說明約束削弱后,變形和殘余應(yīng)力得到釋放,且變形主要以x方向的收縮和面外的變形為主。自由板的各向變形都得到釋放,最大殘余應(yīng)力沒有達(dá)到屈服值,而短邊固接板的x方向和y方向變形都被約束,短邊簡支板短邊上的y方向變形被約束,最大殘余應(yīng)力達(dá)到屈服值。

    自由板靠近焊縫處平行于焊縫的橫截面上的平均殘余應(yīng)力分布如圖6所示,其中,σr為殘余應(yīng)力。自由板垂直于焊縫方向鋼板短跨跨中截面上的平均殘余應(yīng)力分布如圖7所示。

    圖6自由板跨中截面x方向殘余應(yīng)力分布

    Fig.6Distribution of xdirection Residual Stress in

    Midspan Section for Free Plate圖7自由板垂直焊縫方向中間截面處y方向

    殘余應(yīng)力分布

    Fig.7Distribution of ydirection Residual Stress in Middle

    Section Perpendicular to Welded Seam for Free Plate由圖6和圖7得到的自由板焊接殘余應(yīng)力的分布與已有的定性結(jié)論相符合。對整個截面積分,得到整個截面上殘余應(yīng)力的合力為0,符合自由板的條件。同時計(jì)算自由板焊接過程中所施加約束的6個約束的約束反力,得到約束反力均為0,這進(jìn)一步說明了自由板約束的合理性與正確性。

    2.3板件厚度

    在現(xiàn)階段的高層結(jié)構(gòu)設(shè)計(jì)中,厚鋼板剪力墻使用日益廣泛,如天津117大廈中使用的核心筒鋼板剪力墻的最大厚度達(dá)到120 mm。相比于一般的薄板,厚鋼板剪力墻在焊接過程中產(chǎn)生的殘余應(yīng)力和殘余變形具有的特點(diǎn)是值得研究的。

    選取板件厚度t0分別為30,40,50,60 mm的自由板進(jìn)行5層1道焊的對比研究。圖8中給出了焊接過程中自由板不同板件厚度條件下最大面外變形的變化過程。

    圖8自由板不同板件厚度下最大面外變形的變化過程

    Fig.8Changing Processes of Maximum Outofplane

    Deformation of Different Plate Thicknesses for

    Free Plate從圖8可以看出,不同板件厚度在焊接過程中面外變形的發(fā)展趨勢是一樣的,隨著板件厚度的增加,鋼板的面外剛度增加,因此最初正向面外變形的幅值逐漸減小。對于厚板,坡口形狀相同時,最外1層的焊縫寬度大,焊縫收縮量大,最終導(dǎo)致鋼板發(fā)生負(fù)向的面外位移。

    焊縫中心處x方向殘余應(yīng)力沿鋼板厚度方向的分布如圖9所示,其中,σr,x為x方向的殘余應(yīng)力;tz為沿厚度方向的位置。圖9中z方向殘余應(yīng)力表現(xiàn)為中間受拉,兩端受壓,與x方向和y方向的殘余應(yīng)力相比,沿厚度方向的殘余應(yīng)力幅值很小。而x方向殘余應(yīng)力表現(xiàn)為中間受壓,兩端受拉,應(yīng)力幅值隨著板件厚度的增加而變大。這說明對于厚板,垂直于焊縫方向(x方向)的殘余應(yīng)力沿厚度方向是變化的。

    圖9焊縫中心處x方向殘余應(yīng)力沿鋼板厚度方向的分布

    Fig.9Distributions of Residual Stress Along Thickness of

    Steel Plate at Middle of Welded Seam2.4焊接順序

    鋼板焊接的順序影響著鋼板的焊接殘余應(yīng)力和殘余變形的大小。對于本文中研究的鋼板對接焊接,其焊接順序可以分為從上到下、從中間到兩邊和從兩邊到中間3種,如圖10所示。

    圖10焊接順序(單位:mm)

    Fig.10Welding Sequences (Unit:mm)自由板不同焊接順序下V形焊縫焊接最大殘余變形和最大殘余應(yīng)力計(jì)算結(jié)果如表7,8所示。由表7,8可以看出,從兩邊到中間的焊接順序,其面外變形要稍小于其他2種情況,3種不同焊接順序的最終最大殘余應(yīng)力相差不大。

    因此,在具體施工過程中,選取何種焊接順序并沒有統(tǒng)一的標(biāo)準(zhǔn),要根據(jù)實(shí)際情況,以方便施工、縮短工期為原則來選取。

    2.5施工工藝

    焊接過程中的施工工藝有很多種,包括焊接完成后是否保溫及施工過程中是否設(shè)置支撐等。

    表9,10中分別給出了不同保溫條件下焊接層數(shù)為5層且每層1道時自由板進(jìn)行對接焊接的最大殘余變形與最大殘余應(yīng)力。

    對比焊接完成后有無保溫的計(jì)算結(jié)果,無論從焊接殘余變形還是從殘余應(yīng)力上來看,有保溫條件的性能都明顯優(yōu)于無保溫的情況。有保溫時,鋼板

    表7自由板不同焊接順序下V形焊縫最大殘余變形

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

    Tab.7Calculation Results of Maximum Residual

    Deformation of V Welded Seam with Different

    Welding Sequences for Free Platemm焊接順序從上到下從兩邊到中間從中間到兩邊面內(nèi)x方向收縮-2.453-2.028-1.887面外變形7.8237.1817.366表8自由板不同焊接順序下V形焊縫最大殘余應(yīng)力

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

    Tab.8Calculation Results of Maximum Residual

    Stress of V Welded Seam with Different

    Welding Sequences for Free PlateMPa焊接順序從上到下從兩邊到中間從中間到兩邊保溫完成249259256最終冷卻293292293表9自由板不同保溫條件下V形焊縫最大殘余變形

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

    Tab.9Calculation Results of Maximum Residual

    Deformation of V Welded Seam Under Different

    Thermal Conditions for Free Platemm保溫條件有保溫?zé)o保溫面內(nèi)x方向收縮-2.453-2.457面外變形7.8238.035表10自由板不同保溫條件下V形焊縫最大殘余應(yīng)力

    endprint

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

    Tab.10Calculation Results of Maximum Residual

    Stress of V Welded Seam Under Different

    Thermal Conditionsfor Free PlateMPa保溫條件有保溫?zé)o保溫保溫完成249最終冷卻293345在冷卻到保溫溫度后通過熱傳遞使得溫度場分布均勻,溫差梯度小,撤除保溫后整塊鋼板能均勻冷卻到現(xiàn)場環(huán)境溫度;與無保溫的情況相比,有保溫時鋼板不均勻冷卻程度小,從而使得殘余變形和殘余應(yīng)力都較小,說明保溫措施對于控制焊接殘余變形和殘余應(yīng)力來說是必不可少的。

    在鋼板剪力墻焊接施工過程中,往往在2塊需要焊接的墻板上設(shè)置側(cè)向支撐(圖11),用來限制焊接過程中鋼板的面外變形對焊接的影響。在焊接完成后,撤去支撐,鋼板的面外變形增大。

    圖11側(cè)向支撐的設(shè)置(單位:cm)

    Fig.11Arrangement of Lateral Support (Unit:cm)選擇兩短邊簡支的鋼板對接焊接模型,在距焊縫x方向±37 cm處設(shè)置2道側(cè)向支撐,約束該處的面外位移,在保溫完成后撤去側(cè)向支撐,在無支撐情況下冷卻至環(huán)境溫度,其最大殘余變形的計(jì)算結(jié)果見表11。焊接過程中支撐對最大面外殘余變形的影響如圖12所示。

    表11不同支撐條件下V形焊縫最大殘余變形計(jì)算結(jié)果

    Tab.11Calculation Results of Maximum Residual

    Deformations of V Welded Seam Under Different

    Support Conditionsmm支撐設(shè)置有支撐無支撐面內(nèi)x方向收縮-2.258-2.450面外變形-12.8937.731圖12支撐對最大面外殘余變形的影響

    Fig.12Influences of Support on Maximum

    Outofplane Residual Deformation由圖12可以看出,支撐的存在限制了鋼板焊接前期的面外變形,使得鋼板在整個焊接過程中的面外變形都很小,這樣利于焊接的實(shí)施。撤去側(cè)向支撐后,鋼板約束減少,變形恢復(fù)。圖12中在有無支撐情況下最終面外變形方向相反的原因是焊接過程中每層焊縫收縮對鋼板變形趨勢的影響不同,支撐的存在限制了前期鋼板向z軸正向的變形,而第4層與第5層焊縫的收縮會使鋼板向z軸負(fù)向變形,最終導(dǎo)致變形方向相反。3結(jié)語

    (1)對于實(shí)際結(jié)構(gòu)中厚板采用的多層多道對接焊接,不能通過將焊接層數(shù)簡化成單層或雙層以減少計(jì)算量的方法達(dá)到簡化的目的,否則不能保證結(jié)果的正確性,應(yīng)根據(jù)實(shí)際焊接所采用的層數(shù)與道數(shù)進(jìn)行劃分。

    (2)板件的約束強(qiáng)弱影響焊接殘余變形和殘余應(yīng)力的大小。本文中提出了對接焊接的自由板模型,并通過有限元計(jì)算驗(yàn)證了該模型的正確性。計(jì)算結(jié)果表明,厚板的殘余應(yīng)力沿厚度方向變化,其影響不可忽略。

    (3)不同的焊接順序和施工工藝對最終的殘余變形和殘余應(yīng)力大小都有影響。其中,從兩邊到中間焊接的順序略優(yōu)于其他2種焊接順序,在施工過程中,要根據(jù)實(shí)際情況,以方便施工、作業(yè)時間為原則來選取。

    (4)焊接完成后的保溫處理可明顯減小殘余變形和殘余應(yīng)力的幅值。焊接時施加側(cè)向約束或支撐可以很好地控制焊接過程中板件的面外變形,但需要注意支撐撤去后變形的增加。參考文獻(xiàn):

    References:[1]鹿安理,史清宇,趙海燕,等.厚板焊接過程溫度場、應(yīng)力場的三維有限元數(shù)值模擬[J].中國機(jī)械工程,2001,12(2):183186.

    LU Anli,SHI Qingyu,ZHAO Haiyan,et al.Three Dimensional Numerical Simulation of Temperature and Stress Distribution in Welding of Thick Plate[J].China Mechanical Engineering,2001,12(2):183186.

    [2]達(dá)拉伊D.焊接熱效應(yīng)溫度場、殘余應(yīng)力、變形[M].熊第京,譯.北京:機(jī)械工業(yè)出版社,1997.

    RADAJ D.Heat Effects of Welding:Temperature Field,Residual Stress,Deformation[M].Translated by XIONG Dijing.Beijing:China Machine Press,1997.

    [3]倪紅芳,凌祥,涂善東.多道焊三維殘余應(yīng)力場有限元模擬[J].機(jī)械強(qiáng)度,2004,26(2):218222.

    NI Hongfang,LING Xiang,TU Shandong.Simulation of 3D Residual Stresses Field of Multipass Welding by FEM[J].Journal of Mechanical Strength,2004,26(2):218222.

    [4]陳楚,汪建華,楊洪慶.非線性焊接熱傳導(dǎo)的有限元分析和計(jì)算[J].焊接學(xué)報(bào),1983,4(3):139148.

    CHEN Chu,WANG Jianhua,YANG Hongqing.Analysis and Calculation of Nonlinear Welding Heat Conduction by Finite Element Method[J].Transactions of the China Welding Institution,1983,4(3):139148.

    [5]DONG P,BRUST F W.Welding Residual Stresses and Effects of Fracture in Pressure Vessel and Piping Components:A Millennium Review and Beyond[J].Journal of Pressure Vessel Technology,2000,122(3):329338.

    [6]彭崇梅,張啟偉,李元兵.焊接殘余應(yīng)力對鋼桁斜拉橋整體穩(wěn)定性的影響[J].長安大學(xué)學(xué)報(bào):自然科學(xué)版,2011,31(6):3642.

    PENG Chongmei,ZHANG Qiwei,LI Yuanbing.Influence of Welding Residual Stress on Stability of Steel Truss Cablestayed Bridge[J].Journal of Changan University:Natural Science Edition,2011,31(6):3642.

    [7]郭彥林,姜子欽.純彎等截面焊接工字形梁整體穩(wěn)定系數(shù)研究[J].建筑科學(xué)與工程學(xué)報(bào),2012,29(2):8995.

    GUO Yanlin,JIANG Ziqin.Research on Lateral Stability Factor of Welded Isection Beams Under Pure Moment[J].Journal of Architecture and Civil Engineering,2012,29(2):8995.

    endprint

    [8]黃永輝,王榮輝,甘泉.鋼桁梁橋整體節(jié)點(diǎn)焊接殘余應(yīng)力試驗(yàn)[J].中國公路學(xué)報(bào),2011,24(1):8388.

    HUANG Yonghui,WANG Ronghui,GAN Quan.Experiment on Welding Residual Stress of Integral Joint for Steel Truss Bridge[J].China Journal of Highway and Transport,2011,24(1):8388.

    [9]程斌,錢沁,黃仁鋒.部分內(nèi)填混凝土鋼桁梁焊接整體節(jié)點(diǎn)試驗(yàn)[J].中國公路學(xué)報(bào),2013,26(3):117123,136.

    CHENG Bin,QIAN Qin,HUANG Renfeng.Experiment on Partially Concretefilled and Welded Integral Joints in Steel Truss Girders[J].China Journal of Highway and Transport,2013,26(3):117123,136.

    [10]KARLSSON R I,JOSEFSON B L.Three Dimensional Finite Element Analysis of Temperatures and Stresses in a Singlepass Buttwelded Pipe[J].Journal of Pressure Vessel Technology,1990,112(1):7684.

    [11]BROWN S B,SONG H.Implication of Threedimensional Numerical Simulations of Welding of Large Structures[J].Welding Journal,1992,71(2):5562.

    [12]張國智,胡仁喜,陳繼剛.ANSYS 10.0熱力學(xué)有限元分析實(shí)例指導(dǎo)教程[M].北京:機(jī)械工業(yè)出版社,2007.

    ZHANG Guozhi,HU Renxi,CHEN Jigang.Examples of ANSYS 10.0 Finite Element Analysis of Thermodynamics[M].Beijing:China Machine Press,2007.

    [13]郭彥林,趙思遠(yuǎn).深圳灣體育館焊接合攏對結(jié)構(gòu)受力性能的影響分析[J].施工技術(shù),2011,40(3):7175.

    GUO Yanlin,ZHAO Siyuan.Welding Effects on Roof Closure Process of Shenzhenwan Stadium[J].Construction Technology,2011,40(3):7175.

    [14]方洪淵.焊接結(jié)構(gòu)學(xué)[M].北京:機(jī)械工業(yè)出版社,2008.

    FANG Hongyuan.Welding Structure[M].Beijing:China Machine Press,2008.

    endprint

    [8]黃永輝,王榮輝,甘泉.鋼桁梁橋整體節(jié)點(diǎn)焊接殘余應(yīng)力試驗(yàn)[J].中國公路學(xué)報(bào),2011,24(1):8388.

    HUANG Yonghui,WANG Ronghui,GAN Quan.Experiment on Welding Residual Stress of Integral Joint for Steel Truss Bridge[J].China Journal of Highway and Transport,2011,24(1):8388.

    [9]程斌,錢沁,黃仁鋒.部分內(nèi)填混凝土鋼桁梁焊接整體節(jié)點(diǎn)試驗(yàn)[J].中國公路學(xué)報(bào),2013,26(3):117123,136.

    CHENG Bin,QIAN Qin,HUANG Renfeng.Experiment on Partially Concretefilled and Welded Integral Joints in Steel Truss Girders[J].China Journal of Highway and Transport,2013,26(3):117123,136.

    [10]KARLSSON R I,JOSEFSON B L.Three Dimensional Finite Element Analysis of Temperatures and Stresses in a Singlepass Buttwelded Pipe[J].Journal of Pressure Vessel Technology,1990,112(1):7684.

    [11]BROWN S B,SONG H.Implication of Threedimensional Numerical Simulations of Welding of Large Structures[J].Welding Journal,1992,71(2):5562.

    [12]張國智,胡仁喜,陳繼剛.ANSYS 10.0熱力學(xué)有限元分析實(shí)例指導(dǎo)教程[M].北京:機(jī)械工業(yè)出版社,2007.

    ZHANG Guozhi,HU Renxi,CHEN Jigang.Examples of ANSYS 10.0 Finite Element Analysis of Thermodynamics[M].Beijing:China Machine Press,2007.

    [13]郭彥林,趙思遠(yuǎn).深圳灣體育館焊接合攏對結(jié)構(gòu)受力性能的影響分析[J].施工技術(shù),2011,40(3):7175.

    GUO Yanlin,ZHAO Siyuan.Welding Effects on Roof Closure Process of Shenzhenwan Stadium[J].Construction Technology,2011,40(3):7175.

    [14]方洪淵.焊接結(jié)構(gòu)學(xué)[M].北京:機(jī)械工業(yè)出版社,2008.

    FANG Hongyuan.Welding Structure[M].Beijing:China Machine Press,2008.

    endprint

    [8]黃永輝,王榮輝,甘泉.鋼桁梁橋整體節(jié)點(diǎn)焊接殘余應(yīng)力試驗(yàn)[J].中國公路學(xué)報(bào),2011,24(1):8388.

    HUANG Yonghui,WANG Ronghui,GAN Quan.Experiment on Welding Residual Stress of Integral Joint for Steel Truss Bridge[J].China Journal of Highway and Transport,2011,24(1):8388.

    [9]程斌,錢沁,黃仁鋒.部分內(nèi)填混凝土鋼桁梁焊接整體節(jié)點(diǎn)試驗(yàn)[J].中國公路學(xué)報(bào),2013,26(3):117123,136.

    CHENG Bin,QIAN Qin,HUANG Renfeng.Experiment on Partially Concretefilled and Welded Integral Joints in Steel Truss Girders[J].China Journal of Highway and Transport,2013,26(3):117123,136.

    [10]KARLSSON R I,JOSEFSON B L.Three Dimensional Finite Element Analysis of Temperatures and Stresses in a Singlepass Buttwelded Pipe[J].Journal of Pressure Vessel Technology,1990,112(1):7684.

    [11]BROWN S B,SONG H.Implication of Threedimensional Numerical Simulations of Welding of Large Structures[J].Welding Journal,1992,71(2):5562.

    [12]張國智,胡仁喜,陳繼剛.ANSYS 10.0熱力學(xué)有限元分析實(shí)例指導(dǎo)教程[M].北京:機(jī)械工業(yè)出版社,2007.

    ZHANG Guozhi,HU Renxi,CHEN Jigang.Examples of ANSYS 10.0 Finite Element Analysis of Thermodynamics[M].Beijing:China Machine Press,2007.

    [13]郭彥林,趙思遠(yuǎn).深圳灣體育館焊接合攏對結(jié)構(gòu)受力性能的影響分析[J].施工技術(shù),2011,40(3):7175.

    GUO Yanlin,ZHAO Siyuan.Welding Effects on Roof Closure Process of Shenzhenwan Stadium[J].Construction Technology,2011,40(3):7175.

    [14]方洪淵.焊接結(jié)構(gòu)學(xué)[M].北京:機(jī)械工業(yè)出版社,2008.

    FANG Hongyuan.Welding Structure[M].Beijing:China Machine Press,2008.

    endprint

    猜你喜歡
    板件層數(shù)計(jì)算結(jié)果
    填筑層數(shù)對土石壩應(yīng)力變形的影響研究
    上海發(fā)布藥品包裝物減量指南
    康復(fù)(2022年31期)2022-03-23 20:39:56
    基于車身板件定位切割焊接裝置的設(shè)計(jì)
    基于動態(tài)擇優(yōu)組合的板材切割下料算法
    不等高軟橫跨橫向承力索計(jì)算及計(jì)算結(jié)果判斷研究
    甘肅科技(2020年20期)2020-04-13 00:30:40
    MoS2薄膜電子性質(zhì)隨層數(shù)變化的理論研究
    電子制作(2019年11期)2019-07-04 00:34:50
    矩形鋼管截面延性等級和板件寬厚比相關(guān)關(guān)系
    住在哪一層
    鋁合金板件損傷修復(fù)
    超壓測試方法對炸藥TNT當(dāng)量計(jì)算結(jié)果的影響
    一区二区三区乱码不卡18| 麻豆成人午夜福利视频| 极品教师在线视频| 亚洲欧美成人精品一区二区| 日韩av在线免费看完整版不卡| 国产成人免费观看mmmm| 嫩草影院入口| 欧美高清性xxxxhd video| 人人妻人人看人人澡| 身体一侧抽搐| 97超视频在线观看视频| 国产精华一区二区三区| 国产淫语在线视频| 日韩一本色道免费dvd| 亚洲av免费在线观看| 一区二区三区免费毛片| 国产精品麻豆人妻色哟哟久久 | 国产精品福利在线免费观看| 久久精品影院6| 免费搜索国产男女视频| 建设人人有责人人尽责人人享有的 | 狂野欧美白嫩少妇大欣赏| 中文字幕免费在线视频6| 国产免费男女视频| 久久精品久久久久久噜噜老黄 | 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 搡老妇女老女人老熟妇| 免费在线观看成人毛片| 久久久欧美国产精品| 永久免费av网站大全| 亚洲精品乱久久久久久| 最近手机中文字幕大全| 天堂网av新在线| 亚洲电影在线观看av| av免费在线看不卡| 国产淫片久久久久久久久| 国产免费又黄又爽又色| 久久久亚洲精品成人影院| 久久人人爽人人爽人人片va| 亚洲人与动物交配视频| 色5月婷婷丁香| 亚洲精品久久久久久婷婷小说 | 成人鲁丝片一二三区免费| 亚洲成av人片在线播放无| 国产高潮美女av| 久久亚洲精品不卡| 最近手机中文字幕大全| a级一级毛片免费在线观看| 身体一侧抽搐| 最近视频中文字幕2019在线8| 狂野欧美白嫩少妇大欣赏| 成人午夜精彩视频在线观看| av播播在线观看一区| 听说在线观看完整版免费高清| 久久精品夜色国产| 久久精品国产亚洲av涩爱| 国产三级在线视频| 少妇猛男粗大的猛烈进出视频 | 国产视频内射| 九草在线视频观看| 久久精品国产自在天天线| 国产高清有码在线观看视频| 午夜福利成人在线免费观看| 黄色日韩在线| 美女黄网站色视频| 国产麻豆成人av免费视频| 国产精品野战在线观看| 欧美一区二区精品小视频在线| 长腿黑丝高跟| 久久午夜福利片| 中文字幕制服av| 网址你懂的国产日韩在线| 免费播放大片免费观看视频在线观看 | 日本黄色片子视频| 不卡视频在线观看欧美| 久久婷婷人人爽人人干人人爱| av卡一久久| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区 | 亚洲人成网站高清观看| 亚洲激情五月婷婷啪啪| 国产高清国产精品国产三级 | 色综合亚洲欧美另类图片| 午夜老司机福利剧场| 中文字幕人妻熟人妻熟丝袜美| 在线天堂最新版资源| 在线免费观看的www视频| 午夜福利视频1000在线观看| 可以在线观看毛片的网站| 国产爱豆传媒在线观看| 观看美女的网站| 亚洲av不卡在线观看| 黄色一级大片看看| 日韩三级伦理在线观看| 久久久久性生活片| 伦理电影大哥的女人| 99久久精品热视频| 日韩亚洲欧美综合| 久久久久免费精品人妻一区二区| 在线观看一区二区三区| 亚洲最大成人手机在线| 亚洲丝袜综合中文字幕| 国产 一区精品| 欧美成人a在线观看| 99久久成人亚洲精品观看| 夜夜看夜夜爽夜夜摸| 亚洲av成人av| 国产麻豆成人av免费视频| 中文资源天堂在线| 国产伦精品一区二区三区视频9| 久久久久九九精品影院| 18+在线观看网站| 村上凉子中文字幕在线| 一边摸一边抽搐一进一小说| 成人午夜精彩视频在线观看| 最新中文字幕久久久久| 麻豆成人午夜福利视频| 日韩,欧美,国产一区二区三区 | av专区在线播放| 久久久午夜欧美精品| 免费大片18禁| 中文字幕亚洲精品专区| 国产成人精品久久久久久| 99在线人妻在线中文字幕| 国产免费视频播放在线视频 | 国产精品乱码一区二三区的特点| 日韩欧美三级三区| 啦啦啦啦在线视频资源| 三级男女做爰猛烈吃奶摸视频| 黑人高潮一二区| 老司机影院成人| 赤兔流量卡办理| 欧美一区二区精品小视频在线| 最近视频中文字幕2019在线8| 蜜桃亚洲精品一区二区三区| 亚洲电影在线观看av| 久久精品综合一区二区三区| 成年av动漫网址| 熟女电影av网| 99九九线精品视频在线观看视频| 在线播放国产精品三级| 亚洲精品日韩在线中文字幕| 性插视频无遮挡在线免费观看| 久久久久久久亚洲中文字幕| 只有这里有精品99| 一区二区三区乱码不卡18| 秋霞伦理黄片| 久久精品夜夜夜夜夜久久蜜豆| 18禁裸乳无遮挡免费网站照片| 免费观看精品视频网站| 亚洲精品日韩av片在线观看| 毛片女人毛片| 精品少妇黑人巨大在线播放 | 欧美潮喷喷水| 久久精品国产自在天天线| 禁无遮挡网站| 久久精品国产亚洲av涩爱| 日日撸夜夜添| 亚洲一区高清亚洲精品| 亚洲中文字幕日韩| 18禁在线播放成人免费| 直男gayav资源| 永久免费av网站大全| 亚洲av一区综合| 色尼玛亚洲综合影院| 国产免费一级a男人的天堂| 中国国产av一级| 欧美极品一区二区三区四区| 中文亚洲av片在线观看爽| 精品久久国产蜜桃| 精品久久久久久久久久久久久| 一区二区三区高清视频在线| 免费搜索国产男女视频| 日韩在线高清观看一区二区三区| 日韩欧美三级三区| 欧美极品一区二区三区四区| 亚洲精品亚洲一区二区| 国产精品久久视频播放| 国产午夜精品论理片| 美女国产视频在线观看| 色视频www国产| 国产高清不卡午夜福利| 久久综合国产亚洲精品| 精品国内亚洲2022精品成人| 日韩精品有码人妻一区| 国产中年淑女户外野战色| 国产黄片美女视频| 国产亚洲最大av| 亚洲国产成人一精品久久久| 在线观看av片永久免费下载| 欧美成人一区二区免费高清观看| 高清视频免费观看一区二区 | 美女cb高潮喷水在线观看| 亚洲精华国产精华液的使用体验| 97超视频在线观看视频| 免费看av在线观看网站| 国产免费福利视频在线观看| 日日摸夜夜添夜夜爱| 大又大粗又爽又黄少妇毛片口| 中文欧美无线码| 美女内射精品一级片tv| 日韩在线高清观看一区二区三区| 亚洲国产欧洲综合997久久,| 九九爱精品视频在线观看| 日本免费一区二区三区高清不卡| 久久久久九九精品影院| 白带黄色成豆腐渣| 99九九线精品视频在线观看视频| 免费人成在线观看视频色| 看非洲黑人一级黄片| 午夜精品在线福利| 一个人看的www免费观看视频| 国产精品日韩av在线免费观看| 亚洲av电影不卡..在线观看| 边亲边吃奶的免费视频| 18禁裸乳无遮挡免费网站照片| 久久久久免费精品人妻一区二区| 变态另类丝袜制服| 精品人妻一区二区三区麻豆| 哪个播放器可以免费观看大片| 99久国产av精品国产电影| 伦精品一区二区三区| 五月玫瑰六月丁香| 久久精品国产亚洲av天美| 在线观看av片永久免费下载| 国内精品宾馆在线| 夫妻性生交免费视频一级片| 精品免费久久久久久久清纯| 国产成人午夜福利电影在线观看| 国产私拍福利视频在线观看| 久久草成人影院| 纵有疾风起免费观看全集完整版 | 在现免费观看毛片| 欧美日韩综合久久久久久| 日韩强制内射视频| 成人综合一区亚洲| 美女国产视频在线观看| 国产视频内射| 中文字幕av成人在线电影| 日本免费a在线| 欧美丝袜亚洲另类| av在线蜜桃| 国国产精品蜜臀av免费| 免费看日本二区| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 高清视频免费观看一区二区 | 乱系列少妇在线播放| 网址你懂的国产日韩在线| 中文字幕熟女人妻在线| 国产精品野战在线观看| 亚洲国产欧洲综合997久久,| 国产 一区 欧美 日韩| 亚洲精品成人久久久久久| 精品人妻熟女av久视频| 精品国产露脸久久av麻豆 | 男人的好看免费观看在线视频| 国产伦一二天堂av在线观看| 国产亚洲最大av| 男女那种视频在线观看| 久久久久久伊人网av| 亚洲精品456在线播放app| 国产高潮美女av| 日韩欧美精品免费久久| 一个人看的www免费观看视频| 成人漫画全彩无遮挡| 亚洲人成网站高清观看| 免费电影在线观看免费观看| 99久久人妻综合| 国产精品无大码| 国产精品一区二区三区四区久久| 搡老妇女老女人老熟妇| 18禁动态无遮挡网站| 网址你懂的国产日韩在线| 亚洲综合色惰| 一本一本综合久久| 成人欧美大片| 亚洲真实伦在线观看| 两个人视频免费观看高清| 亚洲国产精品久久男人天堂| 国产熟女欧美一区二区| 国产黄a三级三级三级人| 三级国产精品片| 伊人久久精品亚洲午夜| 国产精品av视频在线免费观看| 免费观看的影片在线观看| 99久久精品一区二区三区| 汤姆久久久久久久影院中文字幕 | 国产精品一区二区性色av| 嫩草影院入口| 久久久久久久久久久免费av| 69av精品久久久久久| 91狼人影院| a级毛片免费高清观看在线播放| 天堂网av新在线| 亚洲国产精品专区欧美| 蜜桃久久精品国产亚洲av| 中文欧美无线码| 51国产日韩欧美| 精品熟女少妇av免费看| 最近手机中文字幕大全| 亚洲国产高清在线一区二区三| 一二三四中文在线观看免费高清| 亚洲最大成人中文| 一个人免费在线观看电影| 欧美变态另类bdsm刘玥| 亚洲国产欧美人成| 麻豆一二三区av精品| 欧美精品国产亚洲| 成人美女网站在线观看视频| 日韩三级伦理在线观看| 久久久精品大字幕| 国产高清视频在线观看网站| 国产免费又黄又爽又色| 亚洲av成人av| 可以在线观看毛片的网站| 亚洲自偷自拍三级| 亚洲av.av天堂| 国产人妻一区二区三区在| av.在线天堂| 嫩草影院入口| 欧美日韩综合久久久久久| 国产成人a∨麻豆精品| 亚洲欧美日韩东京热| 青春草亚洲视频在线观看| 一级毛片aaaaaa免费看小| 五月玫瑰六月丁香| 免费观看精品视频网站| 有码 亚洲区| 国产免费又黄又爽又色| 亚洲内射少妇av| 欧美97在线视频| 在线免费观看的www视频| 中文字幕人妻熟人妻熟丝袜美| 草草在线视频免费看| 国产精品久久久久久久久免| 色噜噜av男人的天堂激情| 如何舔出高潮| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人免费av一区二区三区| 亚洲欧美日韩无卡精品| 国产精品伦人一区二区| 婷婷色综合大香蕉| 亚洲美女搞黄在线观看| 亚洲综合色惰| 久久久久久久亚洲中文字幕| 久久人人爽人人爽人人片va| 久久这里有精品视频免费| 国产成人福利小说| 国产又色又爽无遮挡免| 美女xxoo啪啪120秒动态图| 午夜免费男女啪啪视频观看| 噜噜噜噜噜久久久久久91| 久久久欧美国产精品| 男女国产视频网站| 日日摸夜夜添夜夜添av毛片| 搡老妇女老女人老熟妇| 国产精品.久久久| 欧美zozozo另类| 国产伦一二天堂av在线观看| 国产高清有码在线观看视频| 国语自产精品视频在线第100页| 欧美zozozo另类| 少妇高潮的动态图| 久久午夜福利片| 色网站视频免费| 精品久久久久久久末码| 亚洲中文字幕日韩| 国产淫片久久久久久久久| 国产女主播在线喷水免费视频网站 | 免费观看的影片在线观看| av.在线天堂| 亚洲四区av| 精品一区二区三区人妻视频| 亚洲人与动物交配视频| 亚洲乱码一区二区免费版| 久久久国产成人精品二区| 久久国产乱子免费精品| 午夜精品国产一区二区电影 | 看十八女毛片水多多多| 久久午夜福利片| 插逼视频在线观看| 男插女下体视频免费在线播放| 插逼视频在线观看| 久久久国产成人免费| 中文天堂在线官网| 中文字幕精品亚洲无线码一区| 日韩av不卡免费在线播放| 非洲黑人性xxxx精品又粗又长| 国产av在哪里看| 成人鲁丝片一二三区免费| 亚洲成人中文字幕在线播放| 蜜桃久久精品国产亚洲av| 久久99热这里只有精品18| 国产在线一区二区三区精 | 99久国产av精品| 黄片wwwwww| 国产久久久一区二区三区| www.色视频.com| 十八禁国产超污无遮挡网站| 色哟哟·www| 日韩欧美在线乱码| 久久精品久久精品一区二区三区| 免费av不卡在线播放| 一边摸一边抽搐一进一小说| 成年女人看的毛片在线观看| 在现免费观看毛片| 欧美又色又爽又黄视频| 国产精品蜜桃在线观看| 久久久久久久久大av| 观看免费一级毛片| 少妇的逼好多水| 蜜臀久久99精品久久宅男| 建设人人有责人人尽责人人享有的 | 精品久久久久久久久av| 级片在线观看| 国产一区二区在线观看日韩| 久久精品国产自在天天线| 嫩草影院入口| 久久久欧美国产精品| 长腿黑丝高跟| 精品99又大又爽又粗少妇毛片| 哪个播放器可以免费观看大片| 免费看av在线观看网站| 日韩人妻高清精品专区| 男插女下体视频免费在线播放| 国产精品一区二区在线观看99 | 日产精品乱码卡一卡2卡三| 男插女下体视频免费在线播放| 三级男女做爰猛烈吃奶摸视频| 嫩草影院新地址| 亚洲人成网站在线观看播放| 日本猛色少妇xxxxx猛交久久| 又爽又黄a免费视频| 日本色播在线视频| 国产真实乱freesex| 三级国产精品欧美在线观看| 亚洲综合精品二区| 欧美不卡视频在线免费观看| 国产视频首页在线观看| 精品久久久噜噜| 亚洲精品色激情综合| 久久人妻av系列| 日日撸夜夜添| 午夜福利成人在线免费观看| 2022亚洲国产成人精品| 亚洲综合精品二区| 日韩,欧美,国产一区二区三区 | 成人一区二区视频在线观看| 亚洲av成人av| 国产成人a区在线观看| 亚洲国产精品专区欧美| 嘟嘟电影网在线观看| 成人一区二区视频在线观看| 免费电影在线观看免费观看| 欧美成人午夜免费资源| 91av网一区二区| 国产成人aa在线观看| 国产免费福利视频在线观看| 插阴视频在线观看视频| 看十八女毛片水多多多| av在线播放精品| 不卡视频在线观看欧美| 我要搜黄色片| 尾随美女入室| 成人特级av手机在线观看| 亚洲在线观看片| av免费在线看不卡| 久久久久久久久久久丰满| 少妇丰满av| 一级av片app| 国产精品无大码| 国产亚洲一区二区精品| 亚洲精品日韩在线中文字幕| 狂野欧美白嫩少妇大欣赏| 日韩人妻高清精品专区| 国产人妻一区二区三区在| 久久久久久九九精品二区国产| 大香蕉97超碰在线| 在线天堂最新版资源| 亚洲综合精品二区| 精品国产露脸久久av麻豆 | 人人妻人人澡欧美一区二区| 亚洲国产色片| 寂寞人妻少妇视频99o| 丰满少妇做爰视频| 精品不卡国产一区二区三区| 免费搜索国产男女视频| 色综合色国产| 色播亚洲综合网| 黄色一级大片看看| 岛国在线免费视频观看| 看免费成人av毛片| 一级毛片久久久久久久久女| 91精品一卡2卡3卡4卡| 国国产精品蜜臀av免费| 国产免费又黄又爽又色| 日本wwww免费看| 麻豆乱淫一区二区| 国产三级中文精品| 1024手机看黄色片| 欧美区成人在线视频| 欧美成人a在线观看| 欧美区成人在线视频| 18禁裸乳无遮挡免费网站照片| 黄色一级大片看看| 日韩成人av中文字幕在线观看| 黄色配什么色好看| 亚洲五月天丁香| 在线a可以看的网站| 毛片一级片免费看久久久久| 日韩,欧美,国产一区二区三区 | 一级毛片aaaaaa免费看小| 亚洲18禁久久av| 亚洲在线自拍视频| 97超视频在线观看视频| 99久久人妻综合| 少妇丰满av| 亚洲美女搞黄在线观看| 国产一区二区在线观看日韩| 欧美极品一区二区三区四区| 国内揄拍国产精品人妻在线| 国产精品国产三级国产专区5o | 国模一区二区三区四区视频| 午夜精品一区二区三区免费看| 亚洲伊人久久精品综合 | 桃色一区二区三区在线观看| 国产乱人视频| 国产老妇伦熟女老妇高清| 能在线免费观看的黄片| 中文资源天堂在线| 久久精品综合一区二区三区| 成人欧美大片| 久久久久久久午夜电影| 国产精品野战在线观看| 男女边吃奶边做爰视频| 免费一级毛片在线播放高清视频| 又粗又硬又长又爽又黄的视频| 国产大屁股一区二区在线视频| 七月丁香在线播放| av在线老鸭窝| 国产精品,欧美在线| 午夜精品国产一区二区电影 | 中文精品一卡2卡3卡4更新| 亚洲av电影不卡..在线观看| www日本黄色视频网| 又爽又黄无遮挡网站| 国产一区二区三区av在线| 国产69精品久久久久777片| 神马国产精品三级电影在线观看| 大香蕉97超碰在线| 日本黄大片高清| 麻豆成人午夜福利视频| 国内少妇人妻偷人精品xxx网站| 亚洲四区av| 亚洲精品aⅴ在线观看| 精品一区二区三区视频在线| 两个人的视频大全免费| 简卡轻食公司| 永久免费av网站大全| 1024手机看黄色片| 精品国产露脸久久av麻豆 | 国产伦理片在线播放av一区| 国产高清有码在线观看视频| 欧美最新免费一区二区三区| 大香蕉97超碰在线| 亚洲av中文av极速乱| 亚洲av不卡在线观看| 日韩制服骚丝袜av| 国产成人精品一,二区| 亚洲精品国产成人久久av| 高清视频免费观看一区二区 | 少妇裸体淫交视频免费看高清| 免费看光身美女| 狠狠狠狠99中文字幕| 色综合色国产| 国产精品久久视频播放| 久久鲁丝午夜福利片| 亚洲精品影视一区二区三区av| 网址你懂的国产日韩在线| 97热精品久久久久久| 老司机影院毛片| 亚洲欧美清纯卡通| 日韩欧美 国产精品| 久99久视频精品免费| 男的添女的下面高潮视频| 国产片特级美女逼逼视频| 久久久亚洲精品成人影院| 国产亚洲一区二区精品| 日韩中字成人| 国产精品永久免费网站| 亚洲精品色激情综合| 亚洲人成网站在线观看播放| 免费观看人在逋| 国产精品一区二区在线观看99 | 九九热线精品视视频播放| 伦精品一区二区三区| av卡一久久| 永久免费av网站大全| 成年免费大片在线观看| 男女下面进入的视频免费午夜| 真实男女啪啪啪动态图| 国产亚洲午夜精品一区二区久久 | 国产精品熟女久久久久浪| 欧美又色又爽又黄视频| 亚洲av不卡在线观看| 美女黄网站色视频| 看黄色毛片网站| 一区二区三区四区激情视频| 欧美区成人在线视频| 国产精品美女特级片免费视频播放器| 亚洲乱码一区二区免费版| 免费黄网站久久成人精品| 国产成人精品婷婷|