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

    厚板焊接接頭殘余應(yīng)力和損傷分布的研究

    2015-04-02 13:37楊娜等
    關(guān)鍵詞:數(shù)值分析試驗(yàn)研究

    楊娜等

    摘要:對(duì)典型焊接方法下的厚鋼板構(gòu)件焊接性能進(jìn)行實(shí)驗(yàn)和數(shù)值模擬研究.首先采用磁彈儀獲得MPσ標(biāo)定曲線(xiàn),通過(guò)磁彈性法測(cè)量得到厚鋼板焊接接頭殘余應(yīng)力;然后基于ABAQUS及其子程序DFLUX,模擬了厚鋼板對(duì)接焊接全過(guò)程.對(duì)比數(shù)值模擬和實(shí)驗(yàn)結(jié)果表明:構(gòu)件焊接接頭縱向應(yīng)力大于橫向應(yīng)力;沿Y軸的縱向焊接殘余應(yīng)力隨著板厚度的增加而增加,厚度小于60 mm焊件增加尤為明顯.此外,本文耦合了焊接過(guò)程中的溫度場(chǎng)、應(yīng)力場(chǎng),并考慮了積累損傷的影響,完成了焊接構(gòu)件在低周往復(fù)荷載作用下的計(jì)算分析,結(jié)果表明,在低周往復(fù)荷載作用下焊接構(gòu)件中的殘余應(yīng)力對(duì)累積損傷分布有較大影響.

    關(guān)鍵詞:厚鋼板;焊接殘余應(yīng)力;試驗(yàn)研究;積累損傷;數(shù)值分析

    近年來(lái),焊接厚鋼板構(gòu)件廣泛地應(yīng)用于高層和大跨結(jié)構(gòu).例如,國(guó)家體育場(chǎng)“鳥(niǎo)巢”、央視新臺(tái)址大樓、上海環(huán)球金融中心、北京新保利大廈等典型工程都采用了厚度在100 mm以上的鋼板,作為柱腳、主桁架等主要承力構(gòu)件1-2.同時(shí),厚鋼板在焊接過(guò)程中產(chǎn)生的應(yīng)力集中、強(qiáng)度降低、分層等現(xiàn)象也引起了工程界對(duì)焊接鋼構(gòu)件殘余應(yīng)力和損傷分布的廣泛關(guān)注.相對(duì)于普通板件的焊接,厚鋼板焊接構(gòu)件由于受焊接過(guò)程中焊道增多、自身約束增大等因素影響,焊接完成后的殘余應(yīng)力往往會(huì)更大更復(fù)雜,并容易在構(gòu)件中產(chǎn)生初始裂紋和缺陷.在鋼結(jié)構(gòu)中,厚鋼板構(gòu)件在地震狀態(tài)下要承受很大的循環(huán)荷載作用;船舶結(jié)構(gòu)在航行中要承受著連續(xù)的波浪載荷.這樣結(jié)構(gòu)所受的外載荷和殘余應(yīng)力往往都會(huì)疊加,使結(jié)構(gòu)分析變得復(fù)雜3.因此,需要對(duì)厚鋼板構(gòu)件的焊接殘余應(yīng)力分布和外荷載作用下厚板損傷分布進(jìn)行更深入的研究.

    國(guó)內(nèi)外學(xué)者對(duì)焊接殘余應(yīng)力已經(jīng)進(jìn)行了不少研究.美國(guó)的Dean Deng4,ATTARHA M J5使用ABAQUS對(duì)平板和環(huán)形鋼板的焊接進(jìn)行三維數(shù)值模擬研究,并取得相應(yīng)鋼構(gòu)件的殘余應(yīng)力和變形特征;澳大利亞的E.A. FloresJohnson6采用ABAQUS軟件內(nèi)置的JC本構(gòu)模型對(duì)鋼板焊接殘余應(yīng)力及其對(duì)沖擊性能的影響進(jìn)行了研究;張國(guó)棟等人7也用ABAQUS軟件并編制子程序?qū)︿摌?gòu)件的殘余應(yīng)力及其蠕變損傷進(jìn)行了研究.王元清等人8通過(guò)拉伸試驗(yàn)發(fā)現(xiàn)焊縫對(duì)厚鋼板Z向性能有較大影響.但是,目前關(guān)于厚鋼板構(gòu)件焊接殘余應(yīng)力的實(shí)驗(yàn)研究還比較少,關(guān)于循環(huán)荷載作用下殘余應(yīng)力對(duì)構(gòu)件損傷分布影響的闡述也很少.

    本文采用磁彈儀對(duì)厚鋼板對(duì)接焊縫進(jìn)行實(shí)驗(yàn)研究,并借助ABAQUS及其子程序DFLUX,對(duì)焊縫的殘余應(yīng)力進(jìn)行數(shù)值模擬研究和驗(yàn)證.同時(shí),采用子程序UMAT引入鋼積累損傷模型,探討在循環(huán)荷載作用下厚鋼板焊接接頭損傷的分布規(guī)律.

    1厚板焊接殘余應(yīng)力實(shí)驗(yàn)測(cè)試

    對(duì)厚鋼板焊接的現(xiàn)有研究多數(shù)建立在數(shù)值模擬上,但是數(shù)值模擬有很多理論都是建立在假設(shè)簡(jiǎn)化基礎(chǔ)上,而且實(shí)際焊接過(guò)程中受到很多工序影響,因此數(shù)值模擬的全面性和準(zhǔn)確性需要實(shí)驗(yàn)來(lái)驗(yàn)證.

    本次厚鋼板焊接殘余應(yīng)力測(cè)試實(shí)驗(yàn)共分為3部分:焊接過(guò)程,標(biāo)定實(shí)驗(yàn),磁彈法測(cè)定殘余應(yīng)力實(shí)驗(yàn).通過(guò)實(shí)測(cè)對(duì)厚鋼板焊接殘余應(yīng)力分布進(jìn)行初步研究,驗(yàn)證數(shù)值模擬結(jié)果,總結(jié)出能夠適用于實(shí)際工程的殘余應(yīng)力規(guī)律和應(yīng)力模型.

    1.1實(shí)驗(yàn)原理

    實(shí)驗(yàn)選用板材厚度分別為25 mm,40 mm,60 mm和80 mm,試件參數(shù)如表1所示,每個(gè)厚度尺寸取2塊試件,試件坡口尺寸符合《GB 986-88埋弧焊焊縫坡口的基本形式和尺寸》.

    表1焊接參數(shù)

    鋼板選用Q345B鋼,焊接全程都在江蘇中建鋼構(gòu)有限公司進(jìn)行.焊接試件需要預(yù)留6 mm左右的鈍邊,如圖1所示,防止焊接過(guò)程中焊穿;焊接之前切除固定用的馬板,減小焊接過(guò)程中的約束;構(gòu)件需要預(yù)留大小坡口,小坡口深度大約是大坡口的一半,焊接過(guò)程中每焊接3~5層焊道需要進(jìn)行一次翻面,以減小焊接變形.本文采用磁彈法測(cè)量厚板焊接接頭殘余應(yīng)力值.磁彈法是發(fā)展較快、應(yīng)用比較成熟的一種殘余應(yīng)力測(cè)試方法,具有無(wú)損、方便、快速、準(zhǔn)確的特點(diǎn).磁彈法是完全無(wú)損的測(cè)量,本次焊接殘余應(yīng)力的測(cè)試所用儀器為美國(guó)AST公司生產(chǎn)的STRESSCAN500C型磁彈性殘余應(yīng)力測(cè)試儀,儀器本身有預(yù)設(shè)的測(cè)試深度,但基本只能測(cè)量表面應(yīng)力.試件和儀器如圖2所示.

    由于焊接厚鋼板殘余應(yīng)力不僅與工件中所受的應(yīng)力有關(guān),而且還與材料的組織狀態(tài)、硬度、晶粒尺寸及表面狀況有關(guān),所以采用磁彈性方法測(cè)量殘余應(yīng)力時(shí),必須對(duì)不同材料構(gòu)件進(jìn)行標(biāo)定試驗(yàn).標(biāo)定試樣的材質(zhì)、組織狀態(tài)以及表面硬度必須與被測(cè)鋼板相同.試件在拉伸壓縮之前退火消除初始應(yīng)力.

    磁彈儀工作原理在于磁致伸縮效應(yīng),即鐵磁材料在外加磁場(chǎng)的作用下將發(fā)生晶格的彈性變形.通常用巴克豪森噪聲的強(qiáng)度記作MP來(lái)反映材料內(nèi)部的應(yīng)力水平,當(dāng)構(gòu)件應(yīng)力表現(xiàn)為拉應(yīng)力增加時(shí)磁滯回線(xiàn)的斜率和MP值增高,反之,壓應(yīng)力增加則磁滯回線(xiàn)的斜率和MP值降低.

    標(biāo)定實(shí)驗(yàn)主要是測(cè)量巴克豪森噪聲強(qiáng)度記作MP與施加的應(yīng)力值關(guān)系,繪制出標(biāo)定曲線(xiàn)MPσ,該曲線(xiàn)是測(cè)量焊接板件的殘余應(yīng)力的基礎(chǔ).

    實(shí)驗(yàn)測(cè)試通過(guò)磁彈儀在構(gòu)件上測(cè)得MP值,再通過(guò)實(shí)驗(yàn)室得到MP與應(yīng)力的對(duì)應(yīng)關(guān)系標(biāo)定曲線(xiàn),對(duì)照標(biāo)定曲線(xiàn)便可得出焊接構(gòu)件的應(yīng)力值.測(cè)試殘余應(yīng)力的方法是:在工件上放置磁彈性?xún)x,儀器中的探測(cè)線(xiàn)圈檢測(cè)出在磁彈性?xún)x的交變磁場(chǎng)作用下鐵磁材料產(chǎn)生的巴克豪森噪聲MP,從而間接測(cè)出構(gòu)件表面殘余應(yīng)力.

    1.2實(shí)驗(yàn)步驟

    本次實(shí)驗(yàn)由焊接技術(shù)人員制定實(shí)驗(yàn)焊接作業(yè)指導(dǎo)書(shū),將試件的組裝、定位焊、焊前預(yù)熱、焊接過(guò)程控制及要求、焊后保溫等要求做出相應(yīng)規(guī)定,全過(guò)程經(jīng)歷“切割馬板-火焰預(yù)熱→氣保焊打底→清根→自動(dòng)埋弧焊圖3→清根→蓋面焊→保溫圖4→引弧板切割”等,焊接完成試件如圖5所示.測(cè)量應(yīng)力時(shí)的操作如圖6所示,實(shí)驗(yàn)過(guò)程中需要不斷變更方向獲得不同方向的應(yīng)力.

    1.3標(biāo)定實(shí)驗(yàn)結(jié)果

    通過(guò)拉伸和壓縮實(shí)驗(yàn),得到拉壓應(yīng)力與巴克豪森噪聲強(qiáng)度代表值MP的關(guān)系,即標(biāo)定曲線(xiàn),如圖7所示.可取比較有代表性的點(diǎn)輸入儀器,在測(cè)量構(gòu)件的時(shí)候直接輸出殘余應(yīng)力值.

    2厚板焊接接頭殘余應(yīng)力數(shù)值模擬及與

    測(cè)試結(jié)果對(duì)比

    本文基于商用有限元軟件ABAQUS建立焊接鋼厚板構(gòu)件的數(shù)值分析模型,并在模型中利用DFLUX子程序?qū)隚oldak9橢球熱源模型來(lái)模擬焊接熱源的移動(dòng).該熱源模型能夠方便應(yīng)用于薄板、厚板、鋼管等構(gòu)件的焊接過(guò)程熱輸入,其模型方程如式1所示:

    qx,y,z=63fQabcπexp -3x2a2+y2b2+z2c2.1

    式中:x,y,z為模型的整體坐標(biāo)值;Q為輸入功率,f,a,b,c為熱源形狀參數(shù),a,b,c為實(shí)際焊接熔池尺寸縮小20%得到9.本文將重點(diǎn)考慮4個(gè)不同厚度Q345鋼焊接試件的殘余應(yīng)力分布,模型如圖8所示.

    4個(gè)構(gòu)件由相同的材料Q345B鋼組成,焊接參數(shù)與表1相同.材料將參照參考文獻(xiàn)10-11]的16Mn鋼進(jìn)行設(shè)置,隨溫度變化的熱物理性能比熱、熱導(dǎo)率、密度和熱力學(xué)性能屈服強(qiáng)度、彈性模量、泊松比,如圖9所示10-11.熱邊界條件計(jì)入熱輻射和熱交換邊界條件,初始溫度取20°C.已有的研究表明,對(duì)接焊縫的尺寸變化和焊接冶金性能對(duì)焊接力學(xué)性能的影響很小12,因此在有限元建模過(guò)程中予以忽略.

    溫度℃

    本文采用順序耦合法對(duì)厚鋼板模型進(jìn)行殘余應(yīng)力耦合計(jì)算,順序耦合法是先求解溫度場(chǎng),再用溫度場(chǎng)耦合求解應(yīng)力場(chǎng).為了模擬焊料的填充過(guò)程,采用在INP模型文件中加入“model change, addremove”程序段來(lái)實(shí)現(xiàn)生死單元法,對(duì)焊縫單元集合進(jìn)行“殺死”和“激活”設(shè)置.

    由于實(shí)際工程焊接厚鋼板都是平放在地面進(jìn)行的,而且研究表明最符合實(shí)際厚板焊接過(guò)程的邊界條件需要約束焊縫區(qū)13,所以殘余應(yīng)力場(chǎng)的求解過(guò)程邊界條件如圖10所示.底面所有點(diǎn)都需要控制Z方向的豎向移動(dòng),底面點(diǎn)1固定各個(gè)方向的移動(dòng),點(diǎn)2和點(diǎn)4控制Y向和Z向的移動(dòng).

    殘余應(yīng)力可以分為橫向殘余應(yīng)力、縱向殘余應(yīng)力和沿厚度方向的殘余應(yīng)力3種;在ABAQUS有限元分析軟件中,定義平行焊縫方向Y軸的S33為縱向應(yīng)力,垂直焊縫方向X軸的S11為橫向應(yīng)力,沿厚度方向Z軸的應(yīng)力為S22.

    取對(duì)稱(chēng)一半構(gòu)件可看到熱源移動(dòng),從圖11的熱源移動(dòng)中可以看出,整個(gè)焊接階段的熱源最高溫度在1 500~1 600 ℃之間,與實(shí)際情況相符9.冷卻階段大約經(jīng)歷1 h,當(dāng)厚板構(gòu)件整體溫度到達(dá)100 ℃以下可以認(rèn)為冷卻放熱完畢.

    在得到焊接溫度場(chǎng)后,改變模型單元屬性和材料參數(shù),并導(dǎo)入之前模擬的溫度場(chǎng),采用順序耦合法求出最終的殘余應(yīng)力場(chǎng).

    3厚板焊接接頭殘余應(yīng)力測(cè)試結(jié)果對(duì)比

    由于磁彈性?xún)x測(cè)量深度最大只有0.2 mm,而且研究表明,由于最后一道焊對(duì)構(gòu)件應(yīng)力分布的影響比較大,厚鋼板構(gòu)件表面應(yīng)力往往比內(nèi)部應(yīng)力大14.所以本文主要對(duì)焊接厚鋼板構(gòu)件的表面殘余應(yīng)力進(jìn)行測(cè)量,通過(guò)對(duì)構(gòu)件表面的焊接殘余應(yīng)力測(cè)量反映整個(gè)厚鋼板構(gòu)件焊接殘余應(yīng)力.測(cè)點(diǎn)沿長(zhǎng)度方向均勻分布;垂直長(zhǎng)度方向,距離焊縫越近,布置越多.測(cè)量結(jié)果分為垂直焊縫的X軸向應(yīng)力和平行焊縫的Y軸向應(yīng)力.

    3.1沿Y軸應(yīng)力

    首先取平行于焊接方向的沿Y軸縱向分布?xì)堄鄳?yīng)力進(jìn)行研究,以40 mm厚構(gòu)件為例.實(shí)驗(yàn)過(guò)程中測(cè)量了平行焊縫分布的縱向應(yīng)力S33和橫向應(yīng)力S11,它們的分布如圖12所示.

    應(yīng)力MPa

    可以看出,在焊縫周?chē)?,沿Y軸縱向應(yīng)力總體要大于橫向應(yīng)力,縱向應(yīng)力值維持在300 MPa左右;橫向應(yīng)力值比較小,大小為50 MPa左右.焊接接頭邊緣與中間部位應(yīng)力有差異,主要原因是經(jīng)過(guò)火焰切割引弧板的工序?qū)е潞缚p邊緣應(yīng)力出現(xiàn)突變.

    由于縱向應(yīng)力值比較大,將重點(diǎn)討論縱向應(yīng)力S33的分布,取厚鋼板板件的殘余應(yīng)力實(shí)驗(yàn)和數(shù)值模擬結(jié)果對(duì)比,如圖13所示.

    可以看出,縱向應(yīng)力的實(shí)驗(yàn)測(cè)量數(shù)據(jù)點(diǎn)基本在數(shù)值模擬結(jié)果曲線(xiàn)周?chē)梢哉J(rèn)為數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果吻合.沿Y軸縱向應(yīng)力主要呈現(xiàn)中部高,邊緣低的趨勢(shì),和數(shù)值模擬分析結(jié)果一致.

    取厚鋼板大小坡口兩表面的殘余應(yīng)力實(shí)驗(yàn)測(cè)試結(jié)果對(duì)比,如圖14.可以看出,大坡口面縱向應(yīng)力比小坡口面的應(yīng)力大,去除邊緣的應(yīng)力取平均值,可得出40 mm厚鋼板小坡口面應(yīng)力與大坡口面應(yīng)力比值為0.7,而80 mm厚鋼板的比值為0.93.

    沿Y軸距離mm

    圖1440 mm厚構(gòu)件應(yīng)力對(duì)比

    Fig.14 Comparison of stresses in 40 mmthick member

    本文將4個(gè)厚度的構(gòu)件數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比,如圖15和16所示.

    可以看出,25 mm厚構(gòu)件應(yīng)力值比較小,縱向應(yīng)力大都小于270 MPa,焊縫及其熱影響區(qū)HAZ內(nèi)縱向殘余應(yīng)力整體上隨著板厚度的增加而增加,這種趨勢(shì)在25~60 mm之間表現(xiàn)得很明顯,60~80 mm時(shí)有所減弱.主要原因是厚鋼板焊接過(guò)程中板厚度越大,其受到自身的約束越多,所以焊接熱應(yīng)力也越大;但是焊接過(guò)程中每一層先焊的焊道對(duì)周?chē)哪覆挠蓄A(yù)熱作用,而預(yù)熱會(huì)對(duì)焊接熱殘余應(yīng)力有顯著的削弱.

    沿Y軸距離mm

    可以看出,縱向應(yīng)力實(shí)驗(yàn)值在距離焊縫20 mm左右基本達(dá)到峰值,之后在30~80 mm之間數(shù)值迅速減小,在遠(yuǎn)離焊縫的區(qū)域縱向應(yīng)力數(shù)值比較小,甚至出現(xiàn)負(fù)值,但是變化都不大.測(cè)試結(jié)果基本在數(shù)值模擬結(jié)果曲線(xiàn)周?chē)f(shuō)明實(shí)驗(yàn)結(jié)果和數(shù)值模擬結(jié)果比較吻合.

    對(duì)40~80 mm厚的鋼板沿X軸分布的縱向殘余應(yīng)力求均值,并取比較有代表性的5個(gè)點(diǎn)繪成沿X軸分布的縱向殘余應(yīng)力一般分布圖,如圖19所示.對(duì)厚度小于40 mm的厚鋼板縱向殘余應(yīng)力需要將模型數(shù)值適當(dāng)調(diào)低.模型基本體現(xiàn)了沿X軸分布的縱向應(yīng)力的分布特征,在近焊縫區(qū)域應(yīng)力比較大,在遠(yuǎn)離焊縫區(qū)域應(yīng)力比較小.

    沿X軸距離mm

    4構(gòu)件在低周往復(fù)荷載作用下?lián)p傷分布

    損傷是描述材料性能弱化的一個(gè)指標(biāo),通常用一個(gè)變量D來(lái)描述損傷大小,當(dāng)某個(gè)單元的D達(dá)到1時(shí)該單元失效了.本文將采用積累損傷來(lái)評(píng)價(jià)在外部循環(huán)荷載作用下殘余應(yīng)力對(duì)厚板構(gòu)件的影響.

    本文的損傷累積模型基于塑性應(yīng)變進(jìn)行計(jì)算,將在低周往復(fù)荷載作用下的塑性應(yīng)變累積描述成一個(gè)損傷變量D,受力狀態(tài)下的損傷變量可以表示為式215

    D=1-βεpmεpu+∑Ni=1βiεpiεpu.2

    損傷變量的值為N次循環(huán)中的每次循環(huán)中塑性應(yīng)變?chǔ)舙i與極限塑性應(yīng)變?chǔ)舙u的比值之和,β為一個(gè)權(quán)值,對(duì)Q345鋼取值0.011 816.材料彈性模量、屈服強(qiáng)度和強(qiáng)化系數(shù)因損傷累積而產(chǎn)生,隨荷載循環(huán)次數(shù)增加而退化15.

    該損傷模型是基于常用的鋼材雙線(xiàn)型模型,由沈祖炎等15提出適合鋼材的考慮損傷累積的滯回模型.北京交通大學(xué)的楊娜17等基于用戶(hù)子程序UMAT,將積累損傷模型應(yīng)用于ABAQUS的相應(yīng)數(shù)值模擬中.本文首先基于ABAQUS子程序DFLUX對(duì)厚鋼板構(gòu)件焊接殘余應(yīng)力進(jìn)行了順序耦合分析,然后對(duì)存在殘余應(yīng)力的模型施加低周往復(fù)荷載,最后采用UMAT引入累積損傷模型對(duì)焊接殘余應(yīng)力影響下的構(gòu)件累積損傷分布進(jìn)行研究,模擬流程如圖20所示.

    本節(jié)所選焊接構(gòu)件為本文表1所示的兩塊300×250×t mm3鋼板對(duì)接焊接,非焊接構(gòu)件為600×250×t mm3模型,其中t為板厚,各取25 mm,40 mm,0 mm 3試件共6塊板件模型進(jìn)行數(shù)值模擬分析,同厚度構(gòu)件模型網(wǎng)格劃分也相同.試件損傷分析前,提取本文第二部分模擬得到的焊接殘余應(yīng)力場(chǎng),然后導(dǎo)入相同網(wǎng)格的構(gòu)件模型內(nèi),嵌入已驗(yàn)證的子程序UMAT分析構(gòu)件累積損傷

    可以看出,滯回曲線(xiàn)呈飽滿(mǎn)的仿錐形,隨著加載次數(shù)增加,各級(jí)循環(huán)的彈性模量逐漸降低,體現(xiàn)了循環(huán)累積損傷對(duì)滯回性能的影響.通過(guò)對(duì)比可以看出,有殘余應(yīng)力影響的構(gòu)件承載力與剛度較無(wú)殘余應(yīng)力影響構(gòu)件小,說(shuō)明焊接殘余應(yīng)力對(duì)構(gòu)件剛度和承載力有一定的削弱作用.

    圖23為不同厚度構(gòu)件在循環(huán)加載制度下第一象限位移荷載骨架曲線(xiàn)及其對(duì)比.可以看出,骨架曲線(xiàn)表現(xiàn)為隨著位移的增大,荷載平穩(wěn)增大.在循環(huán)加載過(guò)程中有明顯的強(qiáng)化現(xiàn)象,但在加載到10 mm后強(qiáng)化現(xiàn)象不明顯.焊接殘余應(yīng)力對(duì)構(gòu)件承載力有顯著影響,并且隨著厚度增加,影響增大.3組構(gòu)件中,60 mm厚板的殘余應(yīng)力對(duì)構(gòu)件往復(fù)荷載作用下的承載力削弱影響最大.

    由于式2中的損傷值D在ABAQUS中尚未開(kāi)發(fā),本文專(zhuān)門(mén)將損傷值D作為一個(gè)獨(dú)立變量Depvar寫(xiě)入U(xiǎn)MAT中,在ABAQUS模擬完成后采用SDV變量輸出的方法得到構(gòu)件的損傷分布.位移mm

    由本文第3部分可知,殘余應(yīng)力主要集中在焊縫及其熱影響區(qū),因此取構(gòu)件焊縫表面的損傷分布進(jìn)行研究.以起焊點(diǎn)為原點(diǎn),沿焊縫全長(zhǎng)的損傷分布如圖24所示,可知損傷分布也沿焊縫中間大、兩邊小分布.對(duì)比圖中不同厚度構(gòu)件的模擬結(jié)果,可知在相同大小的低周往復(fù)荷載作用下,構(gòu)件厚度越小,承載力相對(duì)越小,屈服后產(chǎn)生的塑性應(yīng)變?cè)蕉?,焊縫表面的損傷值越大.可知構(gòu)件承載力仍然是決定其累積損傷值大小的主要因素.對(duì)比相同構(gòu)件有無(wú)殘余應(yīng)力影響的結(jié)果可知,殘余應(yīng)力對(duì)構(gòu)件焊縫區(qū)有較大的影響,但沒(méi)有達(dá)到破壞損傷值1,可知?dú)堄鄳?yīng)力增大了焊縫表面的累積損傷.

    距離mm

    5結(jié)論

    采用磁彈法對(duì)不同厚度厚板焊接接頭進(jìn)行殘余應(yīng)力測(cè)量,得到了Q345鋼的標(biāo)定曲線(xiàn),也得到厚板焊接接頭殘余應(yīng)力值.同時(shí),基于ABAQUS及其子程序DFLUX對(duì)殘余應(yīng)力進(jìn)行了數(shù)值模擬驗(yàn)證,并修改子程序UMAT引入積累損傷模型,把求解溫度場(chǎng)、殘余應(yīng)力、損傷評(píng)價(jià)連成一體,輸出損傷值,為今后實(shí)際工程中鋼構(gòu)件焊接殘余應(yīng)力的估算和損傷評(píng)價(jià)打下了理論基礎(chǔ).

    1 采用埋弧自動(dòng)焊進(jìn)行焊接的實(shí)驗(yàn)研究和數(shù)值模擬驗(yàn)證,主要研究厚鋼板表面應(yīng)力分布.焊縫處的縱向殘余應(yīng)力總體大于橫向殘余應(yīng)力.焊縫及其熱影響區(qū)縱向殘余應(yīng)力的拉應(yīng)力較大.

    2 沿Y軸向呈中間大邊緣小分布.小坡口面應(yīng)力小于大坡口面應(yīng)力,兩者比值處于0.65~0.95之間.沿焊縫長(zhǎng)度分布的橫向應(yīng)力在邊緣表現(xiàn)為比較大的壓應(yīng)力.沿垂直焊縫方向分布的縱向應(yīng)力在焊縫及其熱影響區(qū)內(nèi)出現(xiàn)峰值,遠(yuǎn)離焊縫變小,并可能出現(xiàn)壓應(yīng)力.

    3 由于隨著板厚度的增加,厚鋼板在焊接過(guò)程中自身約束的不斷增多,沿Y軸分布的縱向殘余應(yīng)力不斷增大,這種趨勢(shì)在板厚處于25~60 mm之間表現(xiàn)得很明顯,但在構(gòu)件厚度大于60 mm后不顯著.橫向殘余應(yīng)力隨厚度變化不顯著.

    4 焊接構(gòu)件在低周往復(fù)剪切位移荷載作用下,計(jì)算得到的骨架曲線(xiàn)有明顯的強(qiáng)化現(xiàn)象;構(gòu)件承載力仍然是決定累積損傷值大小的主要因素;焊接殘余應(yīng)力對(duì)厚鋼板焊縫區(qū)積累損傷分布有一定影響,并隨著厚度增加,影響增大.

    參考文獻(xiàn)

    1]柴昶.厚板鋼材在鋼結(jié)構(gòu)工程中的應(yīng)用及其材性選用J].鋼結(jié)構(gòu),2004,195:47-53.

    CHAI Chang.Application of steel thick plates to steel work constructions and selection of their character J].Steel Construction,2004,195:47-53.In Chinese

    2]王元清,周暉,石永久,等.鋼結(jié)構(gòu)厚板層狀撕裂及其防治措施的研究現(xiàn)狀J].建筑鋼結(jié)構(gòu)進(jìn)展,2010,125: 26-34.

    WANG Yuanqing,ZHOU Hui,SHI yongjiu,et al.A research review on mechanism and prevention of lamellar tearing in thick plate of steel structuresJ].Progress in Steel Building Structures,2010,125:26-34. In Chinese

    3]LEGGAT R H.Residual stresses in welded structuresJ]. International Journal of Pressure Vessels and Piping. 2008,853:144-151.

    4]DENG D,MURAKAWA H.Numerical simulation of temperature field and residual stress in multipass welds in stainless steel pipe and comparison with experimental measurementsJ].Computational Materials Science,2006, 37:269-277.

    5]ATTARHA M J,SATTARIFAR I.Study on welding temperature distribution in thin welded plates through experimental measurements and finite element simulationJ].Journal of Materials Processing Technology,2011211:688-694.

    6]FLORESJOHNSON E A, MURNSKY O.Numerical analysis of the effect of weldinduced residual stress and plastic damage on the ballistic performance of welded steel plateJ].Computational Materials Science,2012,58:131-139.

    7]張國(guó)棟,周昌玉.焊接接頭殘余應(yīng)力及蠕變損傷的有限元模擬J].金屬學(xué)報(bào),2008,447:848-852.

    ZHANG Guodong,ZHOU Changyu.Finite element simulations of welding residual stress and creep damage J].Acta Metallurgica Sinica,2008,447:848-852.In Chinese

    8]王元清,張?jiān)?,石永?鋼厚板母材及其焊接影響區(qū)的Z向拉伸試驗(yàn)J].湖南大學(xué)學(xué)報(bào):自然科學(xué)版,2014,412:26-31.

    WANG Yuanqing,ZHANG Yuanyuan,SHI Yongjiu. Experimental research of the Zdirection tensile properties of thick plate steel and its heat affected zoneJ].Journal of Hunan University:Natural Sciences,2014,412:26-31.In Chinese

    9]郭曉凱.模式搜索法反演多絲埋弧焊雙橢球熱源模型參數(shù)D].上海交通大學(xué)材料科學(xué)與工程學(xué)院,2009.

    GUO Xiaokai.Inversing parameter values of double ellipsoid source model during multiple wires submerged arc welding by using Pattern Search Method D].Shanghai: School of Materials Science and Engineering, Shanghai Jiaotong University,2009.In Chinese

    10]譚真,郭廣文.工程合金熱物性M].北京:冶金工業(yè)出版社,1994:75.

    TAN Zhen,GUO Guangwen.Thermophysical properties of engineering alloys M].Beijing:Metallurgical Industry Press,1994:75.In Chinese

    11]張家榮,趙廷元.工程常用物質(zhì)的熱物理性質(zhì)手冊(cè)M].北京:新時(shí)代出版社,1987:143-145.

    ZHANG Jiarong,ZHAO Tingyuan. Handbook of engineering thermophysical properties of ordinary material M].Beijing:New Age Press,1987:143-145. In Chinese

    12]WU T,CORET M,COMBESCURE A.Numerical simulation of welding induced damage and residual stress of martensitic steel 155PH J].International Journal of Solids and Structures,2008,451819:4973-4989.

    13]李萌盛,嚴(yán)紅丹.不同約束下試板中心堆焊焊接變形的模擬J].熱加工工藝,2009,383:119-123.

    LI Mengsheng,YAN Hongdan.Simulation of resurfacing welding deformation in different constraint J]. Hot Working Technology,2009,383:119-123.In Chinese

    14]TAKAHASH E, LWAI K, SATOH K.A method of measuring triaxial residual stress in heavy section butt weldments J].Transactions of the Japan Welding Society,1979,101:36-45.

    15]董寶,沈祖炎.空間鋼構(gòu)件考慮損傷累積效應(yīng)的恢復(fù)力模型及試驗(yàn)驗(yàn)證J].上海力學(xué),1999,204:341-345.

    DONG Bao,SHEN Zuyan.A hysteretic model and test verification for space steel member with consideration of the damage accumulation effects J].Shanghai Journal of Mechanics,1999,204:341-345.In Chinese

    16]吳艾輝.鋼結(jié)構(gòu)焊接節(jié)點(diǎn)焊材的損傷累積模型及分析D].上海:同濟(jì)大學(xué)土木工程學(xué)院,2007:36-38.

    WU Aihui. Damage accumulation model and analysis of weld material in steel structure welded joint D]. Shanghai: College of Civil Engineering,Tongji University,2007:36-38.In Chinese

    17]張菁菁,楊娜.損傷和殘余應(yīng)力對(duì)鋼構(gòu)件性能的影響J].低溫建筑技術(shù),2010,238:43-45.

    ZHANG Jingjing,YANG Na.The effect of the damage model and residual stress on the performance of the steel componentsJ].Low Temperature Architecture Technology, 2010,238:43-45

    猜你喜歡
    數(shù)值分析試驗(yàn)研究
    單電機(jī)調(diào)節(jié)無(wú)級(jí)變速器試驗(yàn)研究
    加筋土技術(shù)的研究現(xiàn)狀與展望
    彎槽段冰塞形成臨界條件的試驗(yàn)研究
    壓力溶腔對(duì)巖溶隧道施工安全影響的數(shù)值分析
    土與支護(hù)結(jié)構(gòu)相互作用及邊坡穩(wěn)定性分析
    探討補(bǔ)償回彈沖壓件模具設(shè)計(jì)的方法
    控制棒驅(qū)動(dòng)機(jī)構(gòu)傳動(dòng)件耐磨性研究
    南水北調(diào)中線(xiàn)工程輝縣段水泥改性膨脹土工程特性的試驗(yàn)研究
    基于創(chuàng)新和應(yīng)用能力的數(shù)值分析課程教學(xué)研究與實(shí)踐
    一边摸一边抽搐一进一小说| 成人精品一区二区免费| av天堂中文字幕网| 婷婷六月久久综合丁香| 国产精品女同一区二区软件 | 18禁裸乳无遮挡免费网站照片| 久久精品国产99精品国产亚洲性色| 人妻丰满熟妇av一区二区三区| 在线a可以看的网站| 在线观看一区二区三区| 女人十人毛片免费观看3o分钟| 国内精品久久久久久久电影| 91狼人影院| 久久久久久九九精品二区国产| 亚洲精品色激情综合| 日本 欧美在线| 色视频www国产| 国产伦精品一区二区三区四那| 久久精品91蜜桃| 99久久精品热视频| 日韩高清综合在线| 日本一本二区三区精品| 每晚都被弄得嗷嗷叫到高潮| 琪琪午夜伦伦电影理论片6080| 久久午夜亚洲精品久久| 美女被艹到高潮喷水动态| 亚洲中文字幕一区二区三区有码在线看| 他把我摸到了高潮在线观看| a在线观看视频网站| 97超级碰碰碰精品色视频在线观看| 久久午夜亚洲精品久久| 亚洲第一区二区三区不卡| 国产黄片美女视频| 亚洲精品在线美女| 久久久久久久久久成人| 日本免费一区二区三区高清不卡| 婷婷六月久久综合丁香| 天堂√8在线中文| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 又爽又黄无遮挡网站| 蜜桃久久精品国产亚洲av| 99在线视频只有这里精品首页| 国产麻豆成人av免费视频| 日本撒尿小便嘘嘘汇集6| 国产淫片久久久久久久久 | 国产成人aa在线观看| 亚洲第一欧美日韩一区二区三区| 欧美色欧美亚洲另类二区| 精品久久久久久久末码| 久久久久久久精品吃奶| 91久久精品国产一区二区成人| 一本久久中文字幕| 91久久精品国产一区二区成人| 一本久久中文字幕| 国产色婷婷99| 嫩草影视91久久| 国产精品三级大全| 美女高潮的动态| 91狼人影院| 精品人妻视频免费看| 一个人免费在线观看电影| 美女高潮的动态| 成人永久免费在线观看视频| 男女做爰动态图高潮gif福利片| 夜夜爽天天搞| 免费在线观看日本一区| 最好的美女福利视频网| 色播亚洲综合网| 内地一区二区视频在线| 久久精品91蜜桃| 精品午夜福利在线看| 精品欧美国产一区二区三| 国产色婷婷99| 亚洲成人精品中文字幕电影| 亚洲精品影视一区二区三区av| 免费黄网站久久成人精品 | 麻豆国产av国片精品| 欧美潮喷喷水| 日本五十路高清| 淫秽高清视频在线观看| 精品无人区乱码1区二区| 在线国产一区二区在线| 久久精品久久久久久噜噜老黄 | 亚洲国产精品久久男人天堂| 亚洲 欧美 日韩 在线 免费| 国产av不卡久久| 亚洲成人久久爱视频| 免费av观看视频| 2021天堂中文幕一二区在线观| 一级毛片久久久久久久久女| 欧美高清性xxxxhd video| 国产白丝娇喘喷水9色精品| 国产精品av视频在线免费观看| 中国美女看黄片| 中文在线观看免费www的网站| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 午夜两性在线视频| 亚洲av日韩精品久久久久久密| xxxwww97欧美| 搡老熟女国产l中国老女人| 深夜a级毛片| 国产三级在线视频| 国产精品女同一区二区软件 | 国产探花极品一区二区| 国产91精品成人一区二区三区| 97热精品久久久久久| 国产又黄又爽又无遮挡在线| 此物有八面人人有两片| 两个人视频免费观看高清| 欧美另类亚洲清纯唯美| 12—13女人毛片做爰片一| 久久久久久大精品| 久久久精品欧美日韩精品| 一个人免费在线观看的高清视频| 精品99又大又爽又粗少妇毛片 | 宅男免费午夜| bbb黄色大片| 国产午夜精品论理片| 国产一区二区在线观看日韩| 99riav亚洲国产免费| 淫妇啪啪啪对白视频| 91久久精品电影网| 国产高清有码在线观看视频| 99热这里只有是精品50| 国产三级在线视频| 国产精华一区二区三区| 中文资源天堂在线| 成人毛片a级毛片在线播放| 午夜免费成人在线视频| 身体一侧抽搐| 如何舔出高潮| 亚洲精品乱码久久久v下载方式| 99久久九九国产精品国产免费| 亚洲乱码一区二区免费版| 偷拍熟女少妇极品色| 在线看三级毛片| 免费一级毛片在线播放高清视频| 亚洲最大成人中文| 成人国产一区最新在线观看| 婷婷色综合大香蕉| 国产黄a三级三级三级人| 男人的好看免费观看在线视频| 精品一区二区三区av网在线观看| 国产精品亚洲一级av第二区| 美女高潮喷水抽搐中文字幕| 久久亚洲真实| 欧美3d第一页| 夜夜躁狠狠躁天天躁| 老司机深夜福利视频在线观看| 国产高潮美女av| 欧美性感艳星| 精品一区二区三区av网在线观看| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久电影中文字幕| 成人一区二区视频在线观看| 淫秽高清视频在线观看| 国产欧美日韩精品亚洲av| 国产高清有码在线观看视频| 亚洲av美国av| 久久久久久久精品吃奶| av视频在线观看入口| 久久九九热精品免费| 精品熟女少妇八av免费久了| 桃色一区二区三区在线观看| 午夜精品一区二区三区免费看| 亚洲无线观看免费| 成人特级黄色片久久久久久久| av中文乱码字幕在线| 天堂√8在线中文| 69av精品久久久久久| 一进一出抽搐动态| 久久伊人香网站| 久久久久久久久久成人| 亚洲国产高清在线一区二区三| 久久久久亚洲av毛片大全| 国产欧美日韩精品亚洲av| 久久久久国内视频| 国产伦在线观看视频一区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av不卡在线观看| 久久久精品大字幕| 怎么达到女性高潮| 天天一区二区日本电影三级| 欧美一区二区精品小视频在线| 日韩中字成人| 日韩欧美三级三区| 亚洲av熟女| 18禁黄网站禁片午夜丰满| 亚洲久久久久久中文字幕| 一个人看视频在线观看www免费| 亚洲第一区二区三区不卡| 97人妻精品一区二区三区麻豆| 一二三四社区在线视频社区8| 日本 欧美在线| 一级av片app| 亚洲成av人片在线播放无| 最好的美女福利视频网| 色精品久久人妻99蜜桃| 日日摸夜夜添夜夜添av毛片 | 91麻豆精品激情在线观看国产| 日韩有码中文字幕| 欧美最新免费一区二区三区 | 精品人妻一区二区三区麻豆 | 欧美最黄视频在线播放免费| 男女做爰动态图高潮gif福利片| 热99re8久久精品国产| 白带黄色成豆腐渣| 最近视频中文字幕2019在线8| 国内久久婷婷六月综合欲色啪| 日韩欧美国产一区二区入口| 国产亚洲精品久久久com| 91在线观看av| 国产91精品成人一区二区三区| 欧美色欧美亚洲另类二区| 白带黄色成豆腐渣| 一个人免费在线观看电影| 精品一区二区三区av网在线观看| 国产精品av视频在线免费观看| 亚洲成人久久爱视频| 亚洲经典国产精华液单 | 中文字幕免费在线视频6| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片| 日日摸夜夜添夜夜添av毛片 | 日韩av在线大香蕉| 18+在线观看网站| 男女下面进入的视频免费午夜| 午夜精品久久久久久毛片777| 欧美日韩瑟瑟在线播放| 久久久精品欧美日韩精品| 观看免费一级毛片| 国产主播在线观看一区二区| 国产亚洲欧美98| 亚洲电影在线观看av| 我要搜黄色片| 18+在线观看网站| 在线免费观看的www视频| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 久久天躁狠狠躁夜夜2o2o| 国产熟女xx| 色综合亚洲欧美另类图片| 老司机午夜十八禁免费视频| 99久久久亚洲精品蜜臀av| 久久久久久久精品吃奶| 欧美绝顶高潮抽搐喷水| 国产精品一区二区三区四区久久| 久久久久精品国产欧美久久久| 国产午夜精品久久久久久一区二区三区 | 国产精品久久久久久久久免 | 欧美xxxx黑人xx丫x性爽| 国产精品不卡视频一区二区 | 亚洲成av人片在线播放无| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 午夜影院日韩av| 在现免费观看毛片| 亚洲成人免费电影在线观看| 日韩av在线大香蕉| 亚洲最大成人中文| 一个人免费在线观看电影| 窝窝影院91人妻| 他把我摸到了高潮在线观看| 国产成人啪精品午夜网站| 人人妻,人人澡人人爽秒播| 青草久久国产| 在线观看66精品国产| 国产黄色小视频在线观看| 欧美最黄视频在线播放免费| 亚洲18禁久久av| 99热这里只有精品一区| 91av网一区二区| bbb黄色大片| 动漫黄色视频在线观看| 成人无遮挡网站| 9191精品国产免费久久| 成年女人永久免费观看视频| 国内精品久久久久精免费| 一区福利在线观看| 欧美极品一区二区三区四区| 亚洲自拍偷在线| 18禁黄网站禁片午夜丰满| 国产黄片美女视频| 丰满人妻一区二区三区视频av| 久久久久国内视频| 国产精品久久久久久久电影| 亚洲av美国av| www.www免费av| a级一级毛片免费在线观看| 欧美bdsm另类| 午夜免费激情av| 欧美日本视频| 成人国产一区最新在线观看| 免费搜索国产男女视频| 狠狠狠狠99中文字幕| 精品久久久久久久久久免费视频| 村上凉子中文字幕在线| 性色avwww在线观看| 在线播放无遮挡| 五月玫瑰六月丁香| 亚洲欧美精品综合久久99| 欧美日韩黄片免| 99国产精品一区二区蜜桃av| 淫妇啪啪啪对白视频| 国产野战对白在线观看| 亚洲精品影视一区二区三区av| 欧美最黄视频在线播放免费| 亚洲美女黄片视频| 高清毛片免费观看视频网站| 人人妻人人看人人澡| 在线观看一区二区三区| 国产欧美日韩一区二区精品| 熟妇人妻久久中文字幕3abv| 美女黄网站色视频| 老司机午夜十八禁免费视频| 国产精品1区2区在线观看.| 精品久久久久久久久av| 在线播放国产精品三级| 国产精品影院久久| 天天躁日日操中文字幕| 国产精品野战在线观看| 最新中文字幕久久久久| 亚洲av免费高清在线观看| 亚洲,欧美,日韩| 国产主播在线观看一区二区| 在线免费观看不下载黄p国产 | 蜜桃亚洲精品一区二区三区| 欧美不卡视频在线免费观看| 九色成人免费人妻av| 国产中年淑女户外野战色| 午夜视频国产福利| 久久精品91蜜桃| 国产精品爽爽va在线观看网站| 在线看三级毛片| 网址你懂的国产日韩在线| 国产私拍福利视频在线观看| 午夜激情福利司机影院| 亚洲专区国产一区二区| 高清日韩中文字幕在线| av女优亚洲男人天堂| 91狼人影院| 久久欧美精品欧美久久欧美| 日韩欧美免费精品| 国产色爽女视频免费观看| 成人特级黄色片久久久久久久| 天堂√8在线中文| 国产精品美女特级片免费视频播放器| 成人精品一区二区免费| 最新中文字幕久久久久| 熟女人妻精品中文字幕| 成年版毛片免费区| ponron亚洲| 男人舔奶头视频| 舔av片在线| 夜夜夜夜夜久久久久| 国产精品一区二区三区四区免费观看 | 97超视频在线观看视频| 色视频www国产| 亚洲专区国产一区二区| 亚洲国产欧洲综合997久久,| 在线看三级毛片| 黄色日韩在线| 亚洲精品粉嫩美女一区| 国产亚洲欧美98| 久久香蕉精品热| 99久久精品国产亚洲精品| a级一级毛片免费在线观看| 亚洲精品在线观看二区| 丝袜美腿在线中文| 国产激情偷乱视频一区二区| 久久久久久久久大av| 国产成人欧美在线观看| 两个人的视频大全免费| 欧美成人a在线观看| 欧美乱色亚洲激情| 亚洲熟妇中文字幕五十中出| 日韩欧美在线二视频| 99久久久亚洲精品蜜臀av| 91久久精品国产一区二区成人| 伊人久久精品亚洲午夜| 日日干狠狠操夜夜爽| 亚洲va日本ⅴa欧美va伊人久久| 搡老妇女老女人老熟妇| 亚洲经典国产精华液单 | 国产爱豆传媒在线观看| 在线免费观看不下载黄p国产 | 国产精品女同一区二区软件 | 少妇丰满av| 国产精品久久久久久久久免 | 日本黄色片子视频| 日韩中字成人| 级片在线观看| 变态另类丝袜制服| 国产久久久一区二区三区| 身体一侧抽搐| 熟妇人妻久久中文字幕3abv| 91av网一区二区| 日韩欧美精品免费久久 | 国产精品免费一区二区三区在线| 99热6这里只有精品| 午夜两性在线视频| 国产伦人伦偷精品视频| 国产在线男女| 99热只有精品国产| 欧美3d第一页| 欧美日韩黄片免| 桃红色精品国产亚洲av| 在线播放无遮挡| 哪里可以看免费的av片| 日本撒尿小便嘘嘘汇集6| 国语自产精品视频在线第100页| 有码 亚洲区| 国产aⅴ精品一区二区三区波| 午夜福利在线在线| 久久久色成人| 变态另类成人亚洲欧美熟女| 国产黄色小视频在线观看| 国产在线男女| 天堂动漫精品| 黄色丝袜av网址大全| 97人妻精品一区二区三区麻豆| 可以在线观看的亚洲视频| 成年女人看的毛片在线观看| 夜夜爽天天搞| 人妻制服诱惑在线中文字幕| 淫秽高清视频在线观看| 怎么达到女性高潮| 久久精品国产清高在天天线| 国产精品久久久久久久久免 | 免费大片18禁| 成人精品一区二区免费| 国产欧美日韩精品亚洲av| 日韩欧美国产在线观看| 国产成人福利小说| 91狼人影院| 91久久精品国产一区二区成人| 亚洲人成网站在线播放欧美日韩| 熟女人妻精品中文字幕| 久久国产乱子免费精品| 中文字幕人妻熟人妻熟丝袜美| 天堂动漫精品| 日韩欧美在线乱码| 欧美在线一区亚洲| 国产极品精品免费视频能看的| 简卡轻食公司| 精品人妻视频免费看| 久久精品国产亚洲av天美| 深夜精品福利| 亚洲成av人片在线播放无| avwww免费| 亚洲五月天丁香| 亚洲精品在线美女| 99久久精品国产亚洲精品| 色噜噜av男人的天堂激情| 国产 一区 欧美 日韩| 麻豆国产97在线/欧美| 欧美日韩国产亚洲二区| 97热精品久久久久久| 12—13女人毛片做爰片一| 久久99热6这里只有精品| 国产一区二区在线观看日韩| 国产成人a区在线观看| 露出奶头的视频| 99国产精品一区二区三区| 国产精品久久视频播放| 久久热精品热| 一卡2卡三卡四卡精品乱码亚洲| 精品乱码久久久久久99久播| 午夜日韩欧美国产| 色播亚洲综合网| 国产一级毛片七仙女欲春2| 很黄的视频免费| 一区福利在线观看| 免费人成在线观看视频色| 亚洲av一区综合| 麻豆成人av在线观看| av女优亚洲男人天堂| 日韩欧美精品v在线| 99久久成人亚洲精品观看| 床上黄色一级片| 欧美+日韩+精品| 久久午夜福利片| 国产爱豆传媒在线观看| 欧美极品一区二区三区四区| 一二三四社区在线视频社区8| 欧美日韩国产亚洲二区| 国产成人影院久久av| 亚洲精品亚洲一区二区| 国产久久久一区二区三区| 亚洲内射少妇av| 一个人看视频在线观看www免费| 十八禁人妻一区二区| 亚洲av成人不卡在线观看播放网| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久久av| 又黄又爽又免费观看的视频| 久9热在线精品视频| 国产真实伦视频高清在线观看 | 全区人妻精品视频| 免费在线观看亚洲国产| 色综合亚洲欧美另类图片| 麻豆一二三区av精品| 亚洲av熟女| 麻豆久久精品国产亚洲av| 欧美中文日本在线观看视频| 国产精品免费一区二区三区在线| 韩国av一区二区三区四区| 国模一区二区三区四区视频| 日韩精品青青久久久久久| 欧美日韩亚洲国产一区二区在线观看| 国产欧美日韩精品亚洲av| 欧美绝顶高潮抽搐喷水| 一本一本综合久久| 精品人妻视频免费看| 国产在线精品亚洲第一网站| 午夜亚洲福利在线播放| 不卡一级毛片| 亚洲欧美精品综合久久99| 亚洲熟妇中文字幕五十中出| 亚洲最大成人中文| 久久久国产成人免费| 精品人妻视频免费看| 亚洲在线自拍视频| 51午夜福利影视在线观看| 亚洲成人免费电影在线观看| 日本黄色片子视频| 亚洲精品色激情综合| 熟妇人妻久久中文字幕3abv| 高清毛片免费观看视频网站| 欧美日韩黄片免| 一级作爱视频免费观看| 国产亚洲精品av在线| 午夜免费成人在线视频| 国内精品久久久久久久电影| 久久久久亚洲av毛片大全| 99久久精品国产亚洲精品| 亚洲人与动物交配视频| 熟女电影av网| 嫁个100分男人电影在线观看| 日韩中文字幕欧美一区二区| 国产淫片久久久久久久久 | 夜夜爽天天搞| 神马国产精品三级电影在线观看| 国内久久婷婷六月综合欲色啪| 亚洲av电影在线进入| 在线观看免费视频日本深夜| 桃色一区二区三区在线观看| 亚洲精品色激情综合| 伊人久久精品亚洲午夜| 午夜精品久久久久久毛片777| avwww免费| 成人精品一区二区免费| 99riav亚洲国产免费| 无遮挡黄片免费观看| 99热这里只有是精品50| 91在线精品国自产拍蜜月| 每晚都被弄得嗷嗷叫到高潮| 亚洲美女搞黄在线观看 | 久久精品影院6| 成人国产综合亚洲| 无人区码免费观看不卡| 91麻豆精品激情在线观看国产| 国产精品亚洲美女久久久| 在现免费观看毛片| 最新中文字幕久久久久| 18+在线观看网站| 俺也久久电影网| 成人三级黄色视频| 免费在线观看日本一区| 毛片一级片免费看久久久久 | 久久婷婷人人爽人人干人人爱| 免费av观看视频| 美女cb高潮喷水在线观看| av国产免费在线观看| 天堂av国产一区二区熟女人妻| 午夜福利在线观看吧| 亚洲美女黄片视频| 精品国内亚洲2022精品成人| 欧美黑人巨大hd| 日本五十路高清| 亚洲欧美日韩高清在线视频| 夜夜爽天天搞| 国产男靠女视频免费网站| 成人美女网站在线观看视频| 免费人成视频x8x8入口观看| 波多野结衣高清无吗| 一本精品99久久精品77| 天美传媒精品一区二区| 久久热精品热| 老女人水多毛片| 51午夜福利影视在线观看| 成人av一区二区三区在线看| aaaaa片日本免费| 精品久久久久久久久亚洲 | 少妇丰满av| 精品99又大又爽又粗少妇毛片 | 99热只有精品国产| 免费一级毛片在线播放高清视频| 在线观看美女被高潮喷水网站 | 97超视频在线观看视频| 赤兔流量卡办理| 欧美zozozo另类| 成年女人看的毛片在线观看| 国内精品久久久久久久电影| 精品国产三级普通话版| 国产色婷婷99| 免费观看人在逋| 午夜两性在线视频| 无人区码免费观看不卡| 久久精品久久久久久噜噜老黄 | 一级av片app|