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

    多層多道焊接殘余應(yīng)力與變形三維數(shù)值模擬

    2016-12-13 06:53:28蹤雪梅張貴芝張立平占小紅
    電焊機(jī) 2016年3期
    關(guān)鍵詞:平板加密焊縫

    蹤雪梅,張貴芝,黃 松,張立平,占小紅

    (1.徐工集團(tuán)江蘇徐州工程機(jī)械研究院,江蘇徐州221004;2.南京航空航天大學(xué),江蘇南京211106)

    多層多道焊接殘余應(yīng)力與變形三維數(shù)值模擬

    蹤雪梅1,張貴芝1,黃松1,張立平1,占小紅2

    (1.徐工集團(tuán)江蘇徐州工程機(jī)械研究院,江蘇徐州221004;2.南京航空航天大學(xué),江蘇南京211106)

    大型復(fù)雜焊接結(jié)構(gòu)件的主要接頭形式為平板對(duì)接,開(kāi)展對(duì)平板對(duì)接多層多道焊接三維數(shù)值模擬研究十分必要。通過(guò)控制網(wǎng)格尺寸和邊界條件進(jìn)行優(yōu)化來(lái)平衡模擬精度和計(jì)算效率的問(wèn)題,并采用試驗(yàn)測(cè)量和MSC.MARC有限元模擬相結(jié)合方法分析焊接殘余應(yīng)力與變形趨勢(shì)。結(jié)果表明,該多層多道焊接數(shù)值模擬采用位移約束和彈簧約束混和邊界條件,在焊縫最大網(wǎng)格尺寸為2 mm時(shí),計(jì)算效率和精度匹配效果最佳,有限元計(jì)算結(jié)果與試驗(yàn)測(cè)量結(jié)果吻合良好,證明該有限元模型的合理性。

    數(shù)值模擬;多層多道焊;殘余應(yīng)力;變形;彈簧約束

    0 前言

    在中厚板焊接中,考慮到熔寬及熔深大小的局限性,一般采用多層多道焊,而多層多道焊的模擬涉及到金屬的逐步填充問(wèn)題,這在數(shù)學(xué)建模上比較難以實(shí)現(xiàn),這種情況下一般采用將模型簡(jiǎn)化為二維或者不考慮逐步填充,材料一開(kāi)始就全部填入焊道的方法[1]。這些方法從不同程度上降低了模擬精度,都不太理想。

    有限元計(jì)算中的單元“生死”技術(shù)提供了解決這種問(wèn)題的一種方法。其原理是:首先在前處理時(shí)將這些單元“殺死”,即將它們的熱物性參數(shù)變?yōu)楹苄〉闹担诒WC剛度矩陣穩(wěn)定的前提下,將這部分單元的作用減至最小,這樣在實(shí)際計(jì)算中就相當(dāng)于這部分單元不存在。然后在焊接過(guò)程中再根據(jù)焊接速度逐步將這部分單元激活,即將它們的熱物性參數(shù)恢復(fù)原值,相當(dāng)于這部分材料被逐步填充入構(gòu)件[2]。這種技術(shù)在中厚板焊接數(shù)值模擬中得到了廣泛應(yīng)用[3-6]。

    在實(shí)際生產(chǎn)過(guò)程中,涉及中厚板多層多道焊的結(jié)構(gòu)多為復(fù)雜的大型結(jié)構(gòu),采用單元“生死”技術(shù)時(shí),對(duì)數(shù)值模擬的計(jì)算效率要求較高。在模擬精度滿足需求的前提下,如何提高計(jì)算效率是國(guó)內(nèi)外焊接工作者急需解決的研究課題。

    為解決某大型復(fù)雜焊接結(jié)構(gòu)件的焊接應(yīng)力及焊接變形的數(shù)值模擬精度和效率的問(wèn)題,采用平板對(duì)接多層多道焊接接頭來(lái)模擬該結(jié)構(gòu)件的主要接頭形式。通過(guò)控制網(wǎng)格尺寸劃分和邊界條件來(lái)平衡模擬精度和效率的問(wèn)題,并通過(guò)局部網(wǎng)格加密和試驗(yàn)測(cè)試的方法驗(yàn)證仿真結(jié)果,為該大型復(fù)雜焊接結(jié)構(gòu)件的合理模擬提供參考依據(jù)。

    1 焊接試驗(yàn)及殘余應(yīng)力測(cè)試

    試驗(yàn)對(duì)象是開(kāi)60°V型坡口的兩塊平板,材料選用Q235。試驗(yàn)所用焊絲為ER50-6實(shí)心焊絲;采用MAG焊方法,在室溫下進(jìn)行打底焊和填充蓋面焊,保護(hù)氣體為φ(Ar)80%+φ(CO2)20%混合氣體。

    試件平板尺寸為200 mm×200 mm×10 mm,兩層滿焊。試件表面經(jīng)過(guò)手砂輪打磨拋光,砂紙打磨,酒精清洗等表面處理,然后對(duì)測(cè)量表面畫線,分別沿焊縫方向(即縱向近焊縫處)和垂直焊縫方向(即橫向)粘貼應(yīng)變片,采用盲孔法殘余應(yīng)力測(cè)試儀鉆孔,測(cè)量出應(yīng)變值,計(jì)算出焊接殘余應(yīng)力,獲得最大主應(yīng)力σ1。試件物理模型及實(shí)測(cè)點(diǎn)位置如圖1所示。

    圖1 物理模型及殘余應(yīng)力測(cè)試點(diǎn)位置Fig.1Physical model and locations of measurements for residual stress

    2 模型建立

    2.1熱源模型

    MSC.Marc焊接模塊中有兩種缺省的熱源模型:一種是Goldark圓形面熱源,用來(lái)定義二維面熱源。另一種是Goldark雙橢球體熱源,用來(lái)定義三維體熱源[7]。雙橢球熱源模型中熱流密度沿長(zhǎng)軸呈高斯分布前半部分是1/4橢球,后半部分是3/4橢球。前、后橢球的熱分布函數(shù)分別為

    式中Q為輸入熱源功率;ff、fr為熱流密度分布系數(shù);a、b、c1、c2為定義橢球形狀的參數(shù);c1、c2表示前、后半部橢球的長(zhǎng)度;a為影響熔寬,b為影響熔深[8]。

    2.2材料參數(shù)

    材料的熱物理及力學(xué)參數(shù)均隨溫度而變化,焊接是一個(gè)瞬時(shí)高溫過(guò)程,材料參數(shù)變化會(huì)影響模擬結(jié)果的準(zhǔn)確性。模擬材料為Q235,隨溫度變化的部分物性參數(shù)和力學(xué)參數(shù)見(jiàn)圖2。

    2.3網(wǎng)格模型

    對(duì)模型進(jìn)行實(shí)體單元網(wǎng)格劃分,為保證計(jì)算精度,母材網(wǎng)格劃分為3層,為提高計(jì)算速度,將焊縫和熱影響區(qū)的單元網(wǎng)格細(xì)化,而遠(yuǎn)離焊縫區(qū)域的網(wǎng)格較為稀疏,如圖3所示。

    2.4力學(xué)邊界條件

    一般計(jì)算中邊界采用位移約束,對(duì)于平板焊接,如果三個(gè)方向均采用0位移約束,雖消除了模型的剛體位移但約束過(guò)大,經(jīng)過(guò)多次調(diào)試,采用力學(xué)邊界條件包括:位移約束和彈簧約束兩部分,寬度方向在中截面下表面選擇兩個(gè)節(jié)點(diǎn)位置,用來(lái)限制Z向位移,而不影響縱向收縮變形;通過(guò)焊縫背面沿焊縫長(zhǎng)度方向節(jié)點(diǎn)來(lái)限制X向位移,而不影響橫向收縮;工件焊縫背面中心選擇單元表面設(shè)置Y向彈簧約束,如圖3所示。

    2.5局部網(wǎng)格加密算法

    焊接是一個(gè)高溫瞬態(tài)的熱過(guò)程,在焊縫及熱影響位置處,溫度梯度及應(yīng)力梯度均較大。為了得到較為精確的結(jié)果,在焊接及熱影響區(qū)的位置進(jìn)行局部網(wǎng)格加密,在此采用盒形算法加密。

    對(duì)于焊接及熱影響區(qū)中的第i個(gè)單元,如果其變形后的任意一個(gè)節(jié)點(diǎn)落在指定的箱盒范圍內(nèi),則認(rèn)為不滿足箱盒要求,需細(xì)分網(wǎng)格。定義的箱盒按照總體的坐標(biāo)系下的X、Y、Z值給定

    式中X、Y、Z為初始節(jié)點(diǎn)的坐標(biāo);U、V、W為總位移,nodej為單元i的任意一個(gè)節(jié)點(diǎn);f1~f6為定義的箱盒的邊界值[9]。

    圖2 Q235熱-力參量與溫度的關(guān)系Fig.2Relationship between thermo-mechanical parameters and temperature of Q235

    圖3 有限元模型及約束狀態(tài)Fig.3FE model and boundary condition

    2.6解的可靠性

    為使仿真分析的數(shù)據(jù)更具可靠性,必須重新劃分網(wǎng)格進(jìn)行精度收斂分析。由于有限元分析的精度和準(zhǔn)確性與網(wǎng)格的疏密程度有關(guān),在焊道及熱影響區(qū)部位,溫度梯度與應(yīng)力梯度變化較大,為測(cè)試焊縫位置處的網(wǎng)格疏密程度對(duì)解的準(zhǔn)確性影響,對(duì)有限元分析的準(zhǔn)確性進(jìn)行收斂驗(yàn)證。對(duì)焊縫和熱影響區(qū)的網(wǎng)格進(jìn)行劃分,最大網(wǎng)格尺寸分別為5 mm、3 mm和2 mm,并采用MSC.Marc的局部網(wǎng)格加密功能對(duì)焊縫最大尺寸為2 mm的網(wǎng)格進(jìn)行進(jìn)一步加密。圖4c為最大網(wǎng)格尺寸為2 mm的網(wǎng)格模型,圖4d為加密8倍網(wǎng)格的模型,即最大網(wǎng)格尺寸為1 mm的網(wǎng)格模型。

    焊接網(wǎng)格細(xì)化精度收斂情況如表1所示。由表1可知,焊縫最大網(wǎng)格尺寸由5 mm減至1 mm時(shí),計(jì)算時(shí)間呈幾何倍數(shù)增加,第一主應(yīng)力和角變形趨近于某一確定值。最大網(wǎng)格尺寸為5 mm和3 mm的角變形與網(wǎng)格加密后相比,誤差在18.71%。最大網(wǎng)格尺寸為2 mm的角變形與網(wǎng)格加密后相比,誤差僅1.60%,并且計(jì)算時(shí)間僅為加密后的1/6。說(shuō)明網(wǎng)格加密后,網(wǎng)格數(shù)量、節(jié)點(diǎn)數(shù)和計(jì)算時(shí)間大大增加,而加密前后誤差在5%之內(nèi),工程上可以接受。

    圖4 焊縫網(wǎng)格細(xì)化Fig.4Mesh refinement graph of weld bead

    表1 焊縫網(wǎng)格細(xì)化精度收斂分析Tab.1Convergence analysis of mesh refinement accuracy of weld

    3 焊接變形與殘余應(yīng)力分析

    3.1變形結(jié)果分析

    最大網(wǎng)格尺寸為2 mm的仿真變形結(jié)果如圖5所示。對(duì)接平板橫向收縮的主要原因是母材在焊接過(guò)程中首先受熱膨脹,當(dāng)焊縫金屬凝固時(shí),已膨脹的母材金屬必然收縮,這就是對(duì)接接頭橫向收縮的主要組成部分。沿焊縫方向發(fā)生縱向收縮,主要出現(xiàn)在焊縫首尾處。平板垂直焊縫方向上都呈收縮狀態(tài),遠(yuǎn)離焊縫兩側(cè)的母材邊緣向上翹曲,造成角變形,變形趨勢(shì)與實(shí)際情況相吻合。焊縫最大網(wǎng)格尺寸為2 mm時(shí)模擬計(jì)算結(jié)果與實(shí)際焊接變形對(duì)比情況見(jiàn)表2,可以看出,橫向收縮量誤差為8.7%,角變形誤差僅為5.6%。

    圖5 焊后變形云圖Fig.5Contour of welding deformation

    表2 模擬與實(shí)際焊接變形比較Tab.2Comparison of deformation between simulation and measurements

    3.2應(yīng)力結(jié)果分析

    冷至室溫時(shí),對(duì)照試驗(yàn)測(cè)點(diǎn)位置,提取母材面板上表面距焊縫單元邊緣約2 mm的縱向路徑上的縱向殘余應(yīng)力,并與實(shí)際結(jié)果對(duì)比,如圖6所示。可以看出,數(shù)值模擬結(jié)果與實(shí)際測(cè)量獲得的應(yīng)力分布趨勢(shì)基本一致,平行焊縫方向中段處應(yīng)力趨于穩(wěn)定,高于兩端起弧、熄弧段。其值均超過(guò)材料的屈服極限。

    圖6 沿焊縫方向焊趾部位的殘余應(yīng)力分布Fig.6Residual stress distribution of weld toe along the length of weld

    在焊縫中段選取一條垂直于焊縫的路徑,提取表面應(yīng)力分布曲線,并與實(shí)際測(cè)量對(duì)比如圖7所示??梢钥闯觯厔?shì)基本一致,隨距離焊縫中心距離減少,殘余拉應(yīng)力呈遞增趨勢(shì),且焊縫兩側(cè)母材處應(yīng)力呈對(duì)稱分布。不同的是,模擬計(jì)算在焊縫中心的應(yīng)力值有所回落,而實(shí)測(cè)值在中心最大。主要是由于焊縫寬度小,應(yīng)變片有一定尺寸,在焊縫上只能貼一片應(yīng)變片,測(cè)試點(diǎn)比較少,沒(méi)有完全反應(yīng)焊縫上的應(yīng)力分布情況。

    4 結(jié)論

    基于MARC軟件中焊接有限元分析功能,考慮網(wǎng)格尺寸大小、約束、材料等對(duì)模型的影響,對(duì)平板對(duì)接多層多道焊接接頭進(jìn)行三維數(shù)值模擬,給出焊接應(yīng)力場(chǎng)及變形計(jì)算結(jié)果云圖,并通過(guò)局部網(wǎng)格加密的方法和試驗(yàn)測(cè)試的方法對(duì)仿真結(jié)果進(jìn)行了驗(yàn)證,結(jié)果表明:

    (1)通過(guò)網(wǎng)格細(xì)化分析測(cè)試解的可靠性,驗(yàn)證有限元分析收斂可靠,從模擬計(jì)算和試驗(yàn)測(cè)量結(jié)果的對(duì)比分析來(lái)看,兩者吻合較好,說(shuō)明所建立的Marc有限元模型是可行的。

    圖7 工件垂直焊縫方向上的表面殘余應(yīng)力分布Fig.7Surface residual stresses distribution of workpiece in vertical direction of weld

    (2)平板垂直焊縫方向上呈收縮狀態(tài),遠(yuǎn)離焊縫兩側(cè)的母材邊緣向上翹曲,造成角變形,變形趨勢(shì)與實(shí)際情況相一致。

    (3)實(shí)際測(cè)量結(jié)果與仿真模擬獲得應(yīng)力分布趨勢(shì)基本一致:最大焊接殘余應(yīng)力都出現(xiàn)在焊縫中心和靠近焊縫熔合線的地方。平行焊縫方向上中段處的殘余應(yīng)力趨于穩(wěn)定,高于兩端,垂直焊縫方向上的殘余應(yīng)力呈對(duì)稱分布。

    (4)采用位移約束和彈簧約束混合邊界條件,當(dāng)焊縫最大網(wǎng)格尺寸為2 mm時(shí),該多層多道焊接三維數(shù)值模擬的計(jì)算效率和精度匹配效果最佳。

    [1]李國(guó)成,王靖濤,曾靜.大型箱型梁對(duì)接焊的三維數(shù)值模擬[J].工業(yè)建筑,2002,32(11):61-62.

    [2]孟慶國(guó),方洪淵,徐文立,等.考慮金屬逐步填充的多道焊溫度場(chǎng)數(shù)值模擬[J].焊接學(xué)報(bào),2004,25(5):53-59.

    [3]邢淑清,張曉燕,麻永林,等.16Mn特厚鋼板多道焊溫度場(chǎng)數(shù)值模擬[J].電焊機(jī),2012,42(1):18-21.

    [4]黎超文,王勇,韓濤.焊接順序?qū)形接頭殘余應(yīng)力和變形的影響[J].焊接學(xué)報(bào),2011,32(10):37-40.

    [5]陳章蘭,熊云峰,李宗民.船用低溫高強(qiáng)鋼三維多層焊接變形有限元模擬[J].焊接學(xué)報(bào),2008,29(8):109-112.

    [6]羅金華,梁曉燕,王春明,等.中厚板多道焊接溫度場(chǎng)和應(yīng)力場(chǎng)三維數(shù)值模擬[J].電焊機(jī),2006,36(8):32-35.

    [7]MSC.Marc User’s Manual,Volume A,Chapter6.

    [8]John Goldak,Aditya Chakravarti,Malcolm Bibby.A new finite element model for welding heat sources[J].Metallurgical Transactions B,1984(15):299-305.

    [9]陳火紅.MSC.Marc接觸分析培訓(xùn)教程[M].北京:科學(xué)出版社,2004.

    3-D simulation of welding residual stresses and distortion of the multi-pass welding

    ZONG Xuemei1,ZHANG Guizhi1,HUANG Song1,ZHANG Liping1,ZHAN Xiaohong2
    (1.Xuzhou Construction Machinery Group,Jiangsu Xuzhou Construction Machinery Research Institute,Xuzhou 221004 China;2.Najing University of Aeronautics and Astronautics,Nanjing 211106,China)

    Butt joint is the main joint form of large welding structure,so it is essential to carry out 3-D simulation of multi-pass welding of plate welded joints.The paper optimized and balanced the simulation accuracy and computational efficiency by controlling the mesh size and boundary conditions,and then experimental measurement and MSC.Marc finite element simulation were combined to analyze the welding residual stress distribution and distortion trend.The results show that computation efficiency of numerical simulation of multi-pass welding matches well with calculation accuracy,when the mixed boundary condition of displacement constraint and spring constraint are applied to the model,and the weld has maximum element size of 2 mm.Simulation analysis results are in good agreement with actual results,which verify the reasonability of finite element analysis of this model.

    numerical simulation;multi-pass weld;welding residual stress;deformation;spring constraint

    TG404

    A

    1001-2303(2016)03-0117-05

    10.7512/j.issn.1001-2303.2016.03.25

    2015-07-06

    蹤雪梅(1975—),女,江蘇徐州人,學(xué)士,高級(jí)工程師,主要從事工程機(jī)械方面的焊接工藝研究。

    猜你喜歡
    平板加密焊縫
    基于焊縫余高對(duì)超聲波探傷的影響分析
    屬于你的平板電腦
    TP347制氫轉(zhuǎn)油線焊縫裂紋返修
    出彩的立體聲及豐富的畫面層次 華為|平板M6
    一種基于熵的混沌加密小波變換水印算法
    機(jī)器人在輪輞焊縫打磨工藝中的應(yīng)用
    認(rèn)證加密的研究進(jìn)展
    光譜分析在檢驗(yàn)焊縫缺陷中的應(yīng)用
    基于ECC加密的電子商務(wù)系統(tǒng)
    基于格的公鑰加密與證書基加密
    久久久国产精品麻豆| 国产亚洲精品第一综合不卡 | 成人亚洲欧美一区二区av| 成人亚洲精品一区在线观看| 国产成人精品无人区| 交换朋友夫妻互换小说| 精品久久久久久久久亚洲| 亚洲av.av天堂| 精品一区在线观看国产| 久久青草综合色| 丝袜脚勾引网站| 亚洲精品中文字幕在线视频| 中文字幕亚洲精品专区| 九色亚洲精品在线播放| 在线免费观看不下载黄p国产| 亚洲国产精品成人久久小说| 一边摸一边做爽爽视频免费| 九色亚洲精品在线播放| 欧美亚洲日本最大视频资源| 精品亚洲乱码少妇综合久久| 精品一区二区三卡| 秋霞伦理黄片| 熟女av电影| videosex国产| 国产女主播在线喷水免费视频网站| 一区二区三区精品91| av免费观看日本| 国产国拍精品亚洲av在线观看| 男女国产视频网站| 丰满迷人的少妇在线观看| 91在线精品国自产拍蜜月| 天堂8中文在线网| videossex国产| 欧美 亚洲 国产 日韩一| 精品一区二区三区视频在线| 免费日韩欧美在线观看| 国产国拍精品亚洲av在线观看| 久久韩国三级中文字幕| 久久ye,这里只有精品| 9色porny在线观看| 国产精品一区二区在线不卡| 久久久久久久精品精品| 能在线免费看毛片的网站| 亚洲人与动物交配视频| 曰老女人黄片| 内地一区二区视频在线| 国产伦理片在线播放av一区| 国产精品.久久久| 亚洲第一区二区三区不卡| 视频在线观看一区二区三区| 日韩中字成人| 国产精品人妻久久久影院| 五月玫瑰六月丁香| 成人影院久久| 男男h啪啪无遮挡| 国产老妇伦熟女老妇高清| 久久精品久久久久久久性| 国产深夜福利视频在线观看| 国产色婷婷99| 91精品一卡2卡3卡4卡| 亚洲成色77777| 高清黄色对白视频在线免费看| 日韩av不卡免费在线播放| 亚洲精品乱久久久久久| 中文字幕精品免费在线观看视频 | 亚洲性久久影院| 亚洲精品亚洲一区二区| 97在线视频观看| 最近的中文字幕免费完整| 少妇精品久久久久久久| 在线看a的网站| 国产男女内射视频| 日本黄色片子视频| 国产在线视频一区二区| 色婷婷av一区二区三区视频| 女性被躁到高潮视频| 国产成人精品一,二区| 大话2 男鬼变身卡| 中文字幕最新亚洲高清| 18禁动态无遮挡网站| 免费大片黄手机在线观看| 欧美日韩精品成人综合77777| 91在线精品国自产拍蜜月| 午夜精品国产一区二区电影| 国产av精品麻豆| 成人国语在线视频| 欧美xxxx性猛交bbbb| 另类亚洲欧美激情| 欧美日韩视频精品一区| 国产成人免费无遮挡视频| 亚洲精品国产av蜜桃| 五月开心婷婷网| 搡女人真爽免费视频火全软件| 久久国产亚洲av麻豆专区| 高清毛片免费看| 久久久久久久久久人人人人人人| 日韩大片免费观看网站| 啦啦啦视频在线资源免费观看| 久久久久国产精品人妻一区二区| √禁漫天堂资源中文www| 国产av国产精品国产| 久热久热在线精品观看| a级片在线免费高清观看视频| 一个人看视频在线观看www免费| 91午夜精品亚洲一区二区三区| 我的女老师完整版在线观看| 日本av手机在线免费观看| 秋霞在线观看毛片| 亚洲三级黄色毛片| 国产有黄有色有爽视频| 久久久久久久久大av| 国产成人av激情在线播放 | 欧美成人精品欧美一级黄| 国产精品99久久久久久久久| 一边亲一边摸免费视频| 国产免费一区二区三区四区乱码| 亚洲图色成人| 女性生殖器流出的白浆| 欧美亚洲 丝袜 人妻 在线| av又黄又爽大尺度在线免费看| 大片免费播放器 马上看| 亚洲av免费高清在线观看| 久久久久国产精品人妻一区二区| 国产高清三级在线| 亚洲经典国产精华液单| 日产精品乱码卡一卡2卡三| 日韩中文字幕视频在线看片| 人人妻人人添人人爽欧美一区卜| 最后的刺客免费高清国语| 热re99久久精品国产66热6| av在线app专区| 国产毛片在线视频| 午夜福利网站1000一区二区三区| 搡女人真爽免费视频火全软件| 青春草亚洲视频在线观看| 99久久精品国产国产毛片| 国产精品蜜桃在线观看| 一级毛片aaaaaa免费看小| 老司机影院毛片| 日本与韩国留学比较| 一二三四中文在线观看免费高清| 中文字幕久久专区| 少妇被粗大猛烈的视频| 99久国产av精品国产电影| 欧美变态另类bdsm刘玥| 26uuu在线亚洲综合色| av电影中文网址| 久久久精品免费免费高清| 只有这里有精品99| 精品少妇内射三级| 亚洲人成网站在线观看播放| 搡女人真爽免费视频火全软件| 免费日韩欧美在线观看| 国产午夜精品一二区理论片| 肉色欧美久久久久久久蜜桃| 女性生殖器流出的白浆| 十八禁高潮呻吟视频| 国产高清三级在线| 中国美白少妇内射xxxbb| 母亲3免费完整高清在线观看 | 一本大道久久a久久精品| 三上悠亚av全集在线观看| 亚洲人成网站在线观看播放| 肉色欧美久久久久久久蜜桃| 亚洲国产av影院在线观看| 亚洲欧美日韩卡通动漫| 久久国内精品自在自线图片| 久久久久久久国产电影| 不卡视频在线观看欧美| 男女边摸边吃奶| 亚洲av.av天堂| 国语对白做爰xxxⅹ性视频网站| 美女视频免费永久观看网站| 九九在线视频观看精品| 校园人妻丝袜中文字幕| tube8黄色片| 一区二区三区精品91| 精品卡一卡二卡四卡免费| 久久精品熟女亚洲av麻豆精品| 国产精品成人在线| 嫩草影院入口| 日韩精品免费视频一区二区三区 | 国产精品国产三级国产av玫瑰| 在线免费观看不下载黄p国产| 伦理电影大哥的女人| 久久久久视频综合| 日本黄色片子视频| 肉色欧美久久久久久久蜜桃| a级毛片免费高清观看在线播放| 在线亚洲精品国产二区图片欧美 | 一边摸一边做爽爽视频免费| 女性生殖器流出的白浆| 黑人欧美特级aaaaaa片| 久久精品熟女亚洲av麻豆精品| 日本色播在线视频| 精品一区二区三区视频在线| 久久久国产精品麻豆| 多毛熟女@视频| 18禁观看日本| 2022亚洲国产成人精品| 在线观看免费高清a一片| 飞空精品影院首页| 老女人水多毛片| 精品国产国语对白av| av不卡在线播放| 999精品在线视频| 亚洲丝袜综合中文字幕| 精品国产一区二区三区久久久樱花| 国产一区二区三区av在线| 伦理电影免费视频| 黄色视频在线播放观看不卡| 一级爰片在线观看| 久久精品国产亚洲av涩爱| 国产熟女欧美一区二区| av免费在线看不卡| 大片免费播放器 马上看| 少妇熟女欧美另类| av专区在线播放| 三级国产精品欧美在线观看| 日韩不卡一区二区三区视频在线| 国产精品久久久久久av不卡| 麻豆乱淫一区二区| 久久久欧美国产精品| 成人毛片60女人毛片免费| 成人18禁高潮啪啪吃奶动态图 | 欧美三级亚洲精品| 亚洲高清免费不卡视频| 久久久久久久久久成人| 欧美另类一区| 国产免费又黄又爽又色| 汤姆久久久久久久影院中文字幕| 亚洲欧洲精品一区二区精品久久久 | 国产成人精品久久久久久| 少妇人妻久久综合中文| 久久韩国三级中文字幕| 黄色视频在线播放观看不卡| 亚洲国产欧美日韩在线播放| 考比视频在线观看| 一个人看视频在线观看www免费| 日日撸夜夜添| 欧美 亚洲 国产 日韩一| 久久久久久久亚洲中文字幕| 在线看a的网站| 99热网站在线观看| 亚洲人与动物交配视频| 考比视频在线观看| 婷婷色麻豆天堂久久| 亚洲精品日韩在线中文字幕| 欧美日韩成人在线一区二区| 日本爱情动作片www.在线观看| 国精品久久久久久国模美| 老司机影院成人| 中文字幕久久专区| 大香蕉久久成人网| 激情五月婷婷亚洲| 黄色视频在线播放观看不卡| 免费看光身美女| 99热网站在线观看| 一本一本综合久久| 久久av网站| 国产精品成人在线| 啦啦啦中文免费视频观看日本| 亚洲人成网站在线播| 久久精品人人爽人人爽视色| 久久精品久久久久久噜噜老黄| 午夜免费鲁丝| 国产一区二区在线观看av| 能在线免费看毛片的网站| 91精品三级在线观看| 2018国产大陆天天弄谢| 91成人精品电影| 最黄视频免费看| 在线免费观看不下载黄p国产| 你懂的网址亚洲精品在线观看| 亚洲av电影在线观看一区二区三区| 国产黄色视频一区二区在线观看| 午夜91福利影院| 日韩成人av中文字幕在线观看| 在线免费观看不下载黄p国产| 午夜免费男女啪啪视频观看| 在线看a的网站| av一本久久久久| 一边摸一边做爽爽视频免费| 色5月婷婷丁香| 日本-黄色视频高清免费观看| 免费av中文字幕在线| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 国产午夜精品久久久久久一区二区三区| 中文字幕av电影在线播放| 国产在线一区二区三区精| 狠狠婷婷综合久久久久久88av| 18禁裸乳无遮挡动漫免费视频| 免费日韩欧美在线观看| h视频一区二区三区| 纵有疾风起免费观看全集完整版| 9色porny在线观看| 亚洲精品一二三| 777米奇影视久久| 欧美成人午夜免费资源| 日韩精品有码人妻一区| 三上悠亚av全集在线观看| 国产精品偷伦视频观看了| 精品人妻熟女av久视频| 久久久久久久久久久久大奶| 欧美日韩av久久| 欧美三级亚洲精品| 婷婷色综合www| 久久精品人人爽人人爽视色| 热re99久久国产66热| 不卡视频在线观看欧美| 人人妻人人添人人爽欧美一区卜| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡 | 亚洲精品日韩在线中文字幕| 男人爽女人下面视频在线观看| 99国产综合亚洲精品| 啦啦啦视频在线资源免费观看| 尾随美女入室| 国产综合精华液| av视频免费观看在线观看| 久久久国产一区二区| 涩涩av久久男人的天堂| 久久精品久久久久久噜噜老黄| 婷婷成人精品国产| 一区二区三区精品91| 97精品久久久久久久久久精品| 国产有黄有色有爽视频| 一本大道久久a久久精品| 夫妻性生交免费视频一级片| 毛片一级片免费看久久久久| 欧美最新免费一区二区三区| 91精品伊人久久大香线蕉| 久久久久国产网址| 一区二区三区四区激情视频| www.av在线官网国产| 免费大片18禁| 国产毛片在线视频| 天美传媒精品一区二区| 人妻制服诱惑在线中文字幕| 国产av国产精品国产| 在线精品无人区一区二区三| 久久精品久久精品一区二区三区| 伊人久久精品亚洲午夜| 中文字幕av电影在线播放| 精品人妻熟女毛片av久久网站| 狂野欧美激情性xxxx在线观看| 成年女人在线观看亚洲视频| 有码 亚洲区| 黄色欧美视频在线观看| 久久久国产一区二区| 又黄又爽又刺激的免费视频.| 七月丁香在线播放| 一个人免费看片子| 久久久久视频综合| 男女无遮挡免费网站观看| 啦啦啦在线观看免费高清www| 新久久久久国产一级毛片| www.色视频.com| 乱人伦中国视频| 亚洲欧洲国产日韩| av线在线观看网站| 我的女老师完整版在线观看| 国产日韩欧美亚洲二区| av免费在线看不卡| 午夜影院在线不卡| 又粗又硬又长又爽又黄的视频| 国产成人精品一,二区| 国产欧美另类精品又又久久亚洲欧美| 新久久久久国产一级毛片| 美女国产视频在线观看| 亚洲少妇的诱惑av| 精品久久蜜臀av无| 青春草视频在线免费观看| 亚洲av福利一区| 欧美精品高潮呻吟av久久| 99久久中文字幕三级久久日本| 亚洲少妇的诱惑av| 久久久久视频综合| 亚洲精品久久成人aⅴ小说 | 丝袜喷水一区| 高清av免费在线| 日韩欧美一区视频在线观看| 久久影院123| 在线播放无遮挡| 美女xxoo啪啪120秒动态图| 人妻一区二区av| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 免费观看性生交大片5| 国产在线免费精品| 午夜91福利影院| 一级毛片黄色毛片免费观看视频| 亚洲精品国产av成人精品| 永久网站在线| 精品一区二区三区视频在线| 国产 一区精品| 免费看av在线观看网站| 最新中文字幕久久久久| 日本91视频免费播放| 亚州av有码| 国产日韩欧美在线精品| √禁漫天堂资源中文www| 午夜免费男女啪啪视频观看| 18禁在线播放成人免费| 女人久久www免费人成看片| 91精品伊人久久大香线蕉| 黑丝袜美女国产一区| 日韩欧美精品免费久久| 国产女主播在线喷水免费视频网站| 午夜激情av网站| 亚洲精品久久久久久婷婷小说| 韩国高清视频一区二区三区| 久久精品国产自在天天线| 91精品国产九色| 精品国产乱码久久久久久小说| 国产深夜福利视频在线观看| 精品午夜福利在线看| 秋霞伦理黄片| 久久影院123| 精品人妻熟女毛片av久久网站| 91久久精品国产一区二区成人| 欧美日韩综合久久久久久| 成人无遮挡网站| 99视频精品全部免费 在线| 一本一本综合久久| 日韩人妻高清精品专区| 久久久精品区二区三区| 26uuu在线亚洲综合色| 97在线人人人人妻| 久久久国产一区二区| 国产精品国产av在线观看| 99九九线精品视频在线观看视频| 伦理电影大哥的女人| 国语对白做爰xxxⅹ性视频网站| 亚洲国产精品一区二区三区在线| 在线免费观看不下载黄p国产| 欧美日韩在线观看h| 自线自在国产av| 中国三级夫妇交换| 日韩欧美一区视频在线观看| 黄色怎么调成土黄色| 五月开心婷婷网| 亚洲欧美清纯卡通| 满18在线观看网站| 99热全是精品| 26uuu在线亚洲综合色| 精品亚洲成国产av| 一级a做视频免费观看| a级片在线免费高清观看视频| 一级a做视频免费观看| 亚洲国产成人一精品久久久| 色婷婷av一区二区三区视频| 精品一品国产午夜福利视频| 99九九线精品视频在线观看视频| 国产欧美亚洲国产| 美女cb高潮喷水在线观看| 伊人亚洲综合成人网| a级毛色黄片| av福利片在线| 午夜影院在线不卡| 成人漫画全彩无遮挡| 国产精品一区二区三区四区免费观看| 人妻人人澡人人爽人人| 中文天堂在线官网| 妹子高潮喷水视频| 亚洲综合色惰| 欧美变态另类bdsm刘玥| 免费观看无遮挡的男女| 大话2 男鬼变身卡| 国产69精品久久久久777片| 人人妻人人澡人人看| 日本欧美视频一区| 18禁在线播放成人免费| 黑丝袜美女国产一区| 一本一本综合久久| av播播在线观看一区| 老熟女久久久| √禁漫天堂资源中文www| 亚洲无线观看免费| 亚洲国产最新在线播放| 国产精品 国内视频| 另类精品久久| 成人国语在线视频| 久久亚洲国产成人精品v| 国产成人a∨麻豆精品| 在线看a的网站| 性色av一级| 亚洲伊人久久精品综合| 一级毛片aaaaaa免费看小| 观看美女的网站| 成人亚洲欧美一区二区av| 欧美精品亚洲一区二区| 精品国产乱码久久久久久小说| 亚洲国产毛片av蜜桃av| 久久午夜综合久久蜜桃| 亚洲成人av在线免费| 亚洲人与动物交配视频| 成人国产麻豆网| 国产淫语在线视频| 亚洲欧美日韩卡通动漫| 在线看a的网站| 国产无遮挡羞羞视频在线观看| 日本vs欧美在线观看视频| 国产av一区二区精品久久| 啦啦啦视频在线资源免费观看| 寂寞人妻少妇视频99o| 天堂俺去俺来也www色官网| 国产伦精品一区二区三区视频9| 久久国产精品男人的天堂亚洲 | 麻豆精品久久久久久蜜桃| 国产一区二区在线观看av| 99久久精品一区二区三区| 高清不卡的av网站| 成人毛片a级毛片在线播放| 免费看av在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 伦精品一区二区三区| 午夜老司机福利剧场| 亚洲精品自拍成人| 亚洲精品一区蜜桃| 国产高清国产精品国产三级| 色网站视频免费| 秋霞在线观看毛片| 99热网站在线观看| 中国美白少妇内射xxxbb| 最后的刺客免费高清国语| 2022亚洲国产成人精品| 涩涩av久久男人的天堂| 一级二级三级毛片免费看| 最近中文字幕2019免费版| 久久人人爽人人片av| 亚洲国产av影院在线观看| 少妇被粗大的猛进出69影院 | 熟女人妻精品中文字幕| 久久久久精品久久久久真实原创| 久久这里有精品视频免费| 少妇猛男粗大的猛烈进出视频| 母亲3免费完整高清在线观看 | 嫩草影院入口| 久久久亚洲精品成人影院| 国产69精品久久久久777片| 七月丁香在线播放| 色5月婷婷丁香| 高清欧美精品videossex| 91aial.com中文字幕在线观看| 又黄又爽又刺激的免费视频.| 综合色丁香网| 亚洲国产精品一区二区三区在线| 欧美 日韩 精品 国产| 一边亲一边摸免费视频| 岛国毛片在线播放| 51国产日韩欧美| 国产成人91sexporn| 欧美激情 高清一区二区三区| 亚洲欧美清纯卡通| 日本色播在线视频| 久久这里有精品视频免费| 嫩草影院入口| 超色免费av| 日本91视频免费播放| 观看美女的网站| 国产免费一级a男人的天堂| 狠狠精品人妻久久久久久综合| 97在线视频观看| 成人亚洲欧美一区二区av| 天堂俺去俺来也www色官网| 国产欧美另类精品又又久久亚洲欧美| 日韩,欧美,国产一区二区三区| 丝袜脚勾引网站| 精品久久久久久电影网| 午夜福利,免费看| 色婷婷久久久亚洲欧美| 精品久久久久久久久亚洲| 一区二区日韩欧美中文字幕 | 中文乱码字字幕精品一区二区三区| 免费观看a级毛片全部| 秋霞伦理黄片| 丰满少妇做爰视频| 亚洲美女黄色视频免费看| 男人添女人高潮全过程视频| av卡一久久| 在线观看免费日韩欧美大片 | 午夜免费鲁丝| 美女大奶头黄色视频| 乱人伦中国视频| 黄片播放在线免费| 寂寞人妻少妇视频99o| 熟妇人妻不卡中文字幕| 中国国产av一级| 满18在线观看网站| 亚洲图色成人| 在线观看美女被高潮喷水网站| 欧美精品一区二区大全| 欧美成人午夜免费资源| 9色porny在线观看| 日本黄大片高清| 亚洲美女搞黄在线观看| 久久ye,这里只有精品| 亚洲欧美精品自产自拍| 蜜桃久久精品国产亚洲av| 成人18禁高潮啪啪吃奶动态图 | 能在线免费看毛片的网站| 亚洲熟女精品中文字幕| 考比视频在线观看| 午夜老司机福利剧场| 欧美日韩在线观看h| 一区二区日韩欧美中文字幕 | 国产精品一二三区在线看| 啦啦啦在线观看免费高清www| 欧美激情国产日韩精品一区| 丝袜脚勾引网站| av免费在线看不卡| 国产在线免费精品| 久久99热6这里只有精品|