亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        不同對(duì)接工藝下焊接殘余應(yīng)力對(duì)典型縱骨極限強(qiáng)度的影響研究

        2014-06-22 05:08:08李永正岳亞霖韋朋余曾慶波
        船舶力學(xué) 2014年4期
        關(guān)鍵詞:縱骨外板溫度場(chǎng)

        李永正 , 岳亞霖 , 韋朋余, 曾慶波

        (1中國(guó)船舶科學(xué)研究中心,江蘇 無(wú)錫 214082;2江蘇科技大學(xué) 船舶與海洋工程學(xué)院,江蘇 鎮(zhèn)江 212003)

        不同對(duì)接工藝下焊接殘余應(yīng)力對(duì)典型縱骨極限強(qiáng)度的影響研究

        李永正1,2, 岳亞霖1, 韋朋余1, 曾慶波1

        (1中國(guó)船舶科學(xué)研究中心,江蘇 無(wú)錫 214082;2江蘇科技大學(xué) 船舶與海洋工程學(xué)院,江蘇 鎮(zhèn)江 212003)

        隨著模塊化造船的發(fā)展和焊接工藝水平的提高,總段合攏的方式正逐步從傳統(tǒng)的“階梯式”轉(zhuǎn)變?yōu)椤耙坏洱R”方式。文章運(yùn)用ANSYS軟件對(duì)焊接溫度場(chǎng)和應(yīng)力場(chǎng)進(jìn)行間接耦合,對(duì)不同對(duì)接工藝下的船舶縱骨的殘余應(yīng)力分布進(jìn)行數(shù)值模擬,得到相應(yīng)的殘余應(yīng)力的分布規(guī)律;并把數(shù)值模擬得到的殘余應(yīng)力作為初始荷載施加到典型縱骨結(jié)構(gòu)上進(jìn)行極限承載力計(jì)算。通過(guò)對(duì)不同對(duì)接形式的縱骨結(jié)構(gòu)的數(shù)值模擬和試驗(yàn)研究表明,不同對(duì)接形式下焊接殘余應(yīng)力的分布差別較大,對(duì)結(jié)構(gòu)的極限強(qiáng)度有一定的影響,為后續(xù)的相關(guān)研究提供了可靠的依據(jù)。

        一刀齊;對(duì)接焊工藝;焊接殘余應(yīng)力;數(shù)值模擬;極限強(qiáng)度

        1 引 言

        隨著模塊化造船的發(fā)展和焊接工藝水平的提高,總段合攏的方式正逐步從傳統(tǒng)的“階梯式”轉(zhuǎn)變?yōu)椤耙坏洱R”方式,這意味著不同船體分段中縱骨和外板的對(duì)接焊縫將出現(xiàn)在同一位置,從而存在復(fù)雜的焊接殘余應(yīng)力。殘余應(yīng)力對(duì)結(jié)構(gòu)的性能有較大的影響[1],在焊接部位容易產(chǎn)生疲勞裂紋,使得結(jié)構(gòu)的強(qiáng)度和韌性下降,因此對(duì)比這兩種對(duì)接方式下焊接殘余應(yīng)力的大小及分布規(guī)律,能更好地為結(jié)構(gòu)強(qiáng)度計(jì)算提供研究依據(jù)。

        然而現(xiàn)階段對(duì)船體焊接構(gòu)件的殘余應(yīng)力還沒(méi)有簡(jiǎn)單易行的測(cè)量方法[2-5],隨著計(jì)算機(jī)軟硬件性能的提高和數(shù)值模擬技術(shù)的完善,可以運(yùn)用有限元方法對(duì)某些簡(jiǎn)單的焊接過(guò)程進(jìn)行數(shù)值模擬計(jì)算,然后進(jìn)行驗(yàn)證試驗(yàn)來(lái)證明其準(zhǔn)確性及可行性,從而對(duì)其他焊接構(gòu)件不必再進(jìn)行大量的試驗(yàn)就可以達(dá)到研究目的。

        2 焊接數(shù)值模擬的理論基礎(chǔ)

        焊接過(guò)程包含的物理現(xiàn)象異常復(fù)雜,其中包括了焊接時(shí)的傳熱過(guò)程、金屬的熔化和凝固、金屬的相變問(wèn)題、焊接殘余應(yīng)力與變形等[6-7],在焊接數(shù)值模擬過(guò)程中,我們著重分析焊接結(jié)構(gòu)的溫度場(chǎng)和應(yīng)力場(chǎng)的瞬態(tài)變化情況。

        2.1 焊接溫度場(chǎng)的控制方程及邊界條件

        焊接溫度場(chǎng)分析屬于典型的非線性瞬態(tài)熱傳導(dǎo)問(wèn)題的分析,三維瞬態(tài)溫度場(chǎng)控制方程如下:

        式中:c為材料比熱容,ρ為材料密度,T為溫度場(chǎng)分布函數(shù),λ為導(dǎo)熱系數(shù),Q為內(nèi)熱源強(qiáng)度,t為傳熱時(shí)間。

        焊接溫度場(chǎng)熱交換邊界條件:

        有限元單元集成后的矩陣表示形式:

        2.2 焊接溫度場(chǎng)的控制方程及邊界條件

        焊接熱應(yīng)力場(chǎng)可以看成是材料非線性瞬態(tài)問(wèn)題,其中包含復(fù)雜的彈塑性變形。處于彈性或塑性狀態(tài)下的材料,其應(yīng)力應(yīng)變滿足以下本構(gòu)關(guān)系:

        式中:dσ為應(yīng)力增量,dε為應(yīng)變?cè)隽?,dt為溫度增量,D為彈性或彈塑性矩陣,c為材料比熱容。

        在ANSYS中,我們可以采用兩種方式來(lái)耦合焊接溫度場(chǎng)和應(yīng)力場(chǎng),從而求解焊接殘余應(yīng)力,即:間接耦合法和直接耦合法[8],本文采用前者,即先求解溫度場(chǎng),得到準(zhǔn)確的節(jié)點(diǎn)溫度后以體載荷的方式施加到應(yīng)力場(chǎng)求解的過(guò)程。

        3 典型縱骨對(duì)接焊數(shù)值模擬

        3.1 有限元模型

        為了研究縱骨及其帶板在不同對(duì)接工藝下焊接殘余應(yīng)力的大小及分布,本文選取了兩種不同對(duì)接焊縫布置方式的模型進(jìn)行研究計(jì)算,如表1所示。

        同位式:采用同一截面對(duì)接焊縫,面板1 500×250×8,骨材為10號(hào)球扁鋼,面板及骨材焊縫均位于距左端750 mm處的同一截面;有限元模型網(wǎng)格劃分如圖1所示。

        交錯(cuò)式:采用交錯(cuò)布置對(duì)接焊縫,面板1 500×250×8,骨材為10號(hào)球扁鋼,面板焊縫位于距左端600 mm處,縱骨焊縫位于距左端900 mm處;有限元模型網(wǎng)格劃分如圖2所示。

        表1 模型幾何尺寸Tab.1 Model geometry

        圖1 同位式有限元模型Fig.1 ‘In-plane welded connection’FEM model

        圖2 交錯(cuò)式有限元模型Fig.2 ‘Step join’ FEM model

        3.2 材料屬性

        模型材料采用某種高強(qiáng)度鋼[11],材料的物理性能參數(shù)及力學(xué)性能參數(shù)均隨溫度的變化而變化,表2和表3給出了各溫度下的材料熱物理性能參數(shù)和力學(xué)性能參數(shù)。

        表2 材料熱物理性能參數(shù)Tab.2 Material thermal physical property parameters

        表3 材料在不同溫度下的應(yīng)力應(yīng)變Tab.3 The stress and strain of materials under different temperatures

        3.3 縱骨焊接殘余應(yīng)力分布

        縱骨焊接殘余應(yīng)力分布如圖3所示,從圖3(a)中可以看出在同一截面對(duì)接焊縫處的焊接殘余應(yīng)力達(dá)到峰值,維持在600-700 MPa之間,隨著距離焊縫越來(lái)越遠(yuǎn),殘余應(yīng)力隨之降低;在骨材與外板的結(jié)合處也分布著較大的焊接殘余應(yīng)力,其值約為600 MPa左右。從圖3(b)中可以看出在交錯(cuò)式對(duì)接工藝下,骨材和外板相應(yīng)的對(duì)接焊縫處亦存在較大焊接殘余應(yīng)力,維持在600-700 MPa之間,且相互影響;但骨材和外板結(jié)合處不受其影響。

        圖3 不同對(duì)接工藝下焊接殘余應(yīng)力(mises)云圖Fig.3 The welding residual stress under different docking process

        圖4 同一截面焊接工藝下焊縫處殘余應(yīng)力分布曲線Fig.4 Weld residual stress distribution curve under the same section welding plate

        圖5 交錯(cuò)焊接工藝下外板焊縫處殘余應(yīng)力分布曲線Fig.5 Weld residual stress distribution curve under the step section welding plate

        為了更直觀地分析和比較焊縫附近的殘余應(yīng)力的分布情況,沿外板和縱骨焊縫的溶合線處的殘余應(yīng)力分布情況如圖4-5所示。

        從圖4中可以得出同位式對(duì)接結(jié)構(gòu)的外板和縱骨焊縫溶合線處的殘余應(yīng)力分布情況,外板焊縫處的縱向應(yīng)力是三個(gè)方向中最大的應(yīng)力,大小約在300 MPa之間,峰值約為560 MPa。其次為X方向的橫向殘余應(yīng)力,峰值大小約為330 MPa,在板兩端處存在壓應(yīng)力。Z向的應(yīng)力則基本為恒定的壓應(yīng)力,大小約在120 MPa左右。

        圖5為交錯(cuò)式對(duì)接結(jié)構(gòu)外板和縱骨焊縫溶合線處的殘余應(yīng)力分布情況,外板焊縫溶合線處的焊接殘余應(yīng)力穩(wěn)定分布在220-280 MPa之間,比焊縫內(nèi)的焊接殘余應(yīng)力600 MPa低。其中縱向殘余應(yīng)力(平行于焊縫方向的SY)大小為230-300 MPa之間,在板兩端處有所降低。橫向殘余應(yīng)力(垂直焊縫方向的SX)與平板對(duì)接焊接的殘余應(yīng)力分布相似,在焊縫中間處存在較低的拉應(yīng)力約為100 MPa,在板的端部則存在壓應(yīng)力,在與外板垂直的Z方向存在一定的拉應(yīng)力,大小約為110 MPa,分布較為均勻。

        4 焊接殘余應(yīng)力對(duì)典型縱骨極限承載能力的影響

        4.1 極限強(qiáng)度的有限元分析

        典型縱骨在不同對(duì)接工藝下的焊接殘余應(yīng)力,作為初始應(yīng)力施加在有限元模型上,再利用ANSYS進(jìn)行非線性極限強(qiáng)度分析,獲取結(jié)構(gòu)的極限承載能力。

        如圖6所示,不同結(jié)構(gòu)形式受拉伸載荷作用下的載荷因子—位移曲線,在加載初始階段,各單元材料處于彈性階段,載荷位移曲線成線性分布;隨著載荷不斷增加,各部分單元先后進(jìn)入塑性階段,由于材料塑性強(qiáng)化模量小于彈性模量,載荷位移曲線斜率逐漸變小,當(dāng)結(jié)構(gòu)達(dá)到極限狀態(tài)時(shí)曲線接近一條直線,此時(shí),位移不斷增加,載荷不再增加反而有所降低。圖6中的三條曲線分別為無(wú)對(duì)接焊縫、交錯(cuò)式對(duì)接焊縫和同一截面對(duì)接焊縫三種形式下結(jié)構(gòu)的載荷因子—位移曲線。

        圖6 不同對(duì)接形式下的載荷-位移曲線Fig.6 Load-stress curve under different welding techniques

        4.2 試驗(yàn)結(jié)果對(duì)比

        利用1 500噸靜態(tài)試驗(yàn)系統(tǒng)開展不同對(duì)接形式典型結(jié)構(gòu)的極限強(qiáng)度試驗(yàn),荷載從0 kN開始施加直至結(jié)構(gòu)斷裂失效,試驗(yàn)得到的最大荷載就是該結(jié)果的極限荷載。

        圖7 不同對(duì)接形式的極限強(qiáng)度試驗(yàn)Fig.7 Ultimate strength test of different joint forms

        表4列出不同對(duì)接形式下縱骨結(jié)構(gòu)的極限載荷:整體式縱骨結(jié)構(gòu)不存在對(duì)接焊縫,無(wú)殘余應(yīng)力影響,在受載過(guò)程中,截面的應(yīng)力分布均勻,逐漸趨于材料的極限強(qiáng)度,結(jié)構(gòu)極限承載能力的數(shù)值模擬值結(jié)果為1 917 kN,試驗(yàn)值為為1 964 kN,相對(duì)誤差為2.4%。交錯(cuò)式縱骨結(jié)構(gòu)的焊接殘余應(yīng)力分布在外板與骨材的焊縫處,在加載階段焊縫處部分單元中的的拉應(yīng)力緩慢增加,逐步達(dá)到塑性階段,其他區(qū)域的應(yīng)力隨著載荷的增加,應(yīng)力值增加,直到與焊縫區(qū)的應(yīng)力值均勻相等,結(jié)構(gòu)極限載荷的數(shù)值模擬結(jié)果為1 835 kN,與試驗(yàn)值1 874 kN相差在2.1%。同位式縱骨結(jié)構(gòu)中的焊接殘余應(yīng)力集中在外板與縱骨的焊縫交界處,其應(yīng)力變化趨勢(shì)與交錯(cuò)式相近,結(jié)構(gòu)極限載荷的數(shù)值模擬結(jié)果為1 819 kN,與試驗(yàn)值相差1.3%。

        表4 結(jié)構(gòu)極限載荷Tab.4 Structural ultimate load

        不同對(duì)接形式產(chǎn)生的殘余應(yīng)力對(duì)結(jié)構(gòu)極限強(qiáng)度有一定程度影響,其中交錯(cuò)式對(duì)接結(jié)構(gòu)的極限荷載與沒(méi)有對(duì)接焊縫的整體式相比,數(shù)值模擬結(jié)果相對(duì)差值為4.2%,試驗(yàn)結(jié)果相對(duì)差值為4.6%;同位式對(duì)接結(jié)構(gòu)的極限荷載與沒(méi)有對(duì)接焊縫的整體式相比,數(shù)值模擬結(jié)果相對(duì)差值為5.1%,試驗(yàn)結(jié)果相對(duì)差值為6.3%。通過(guò)上述結(jié)果表明,不同對(duì)接形式產(chǎn)生的殘余應(yīng)力對(duì)結(jié)構(gòu)極限強(qiáng)度均有不同程度的影響,其中交錯(cuò)式對(duì)接形式的影響相對(duì)較小,同位式對(duì)接形式的影響較大。

        5 結(jié) 論

        (1)交錯(cuò)式對(duì)接焊縫產(chǎn)生的殘余應(yīng)力分布存在縱向焊接殘余應(yīng)力(沿焊縫方向)和橫向殘余應(yīng)力(垂直焊縫方向),且縱向殘余應(yīng)力的值在焊縫及其附近具有最大值,大小接近材料的屈服點(diǎn)。而垂直焊縫方向的橫向殘余應(yīng)力值的大小比縱向殘余應(yīng)力值小,且向板兩端處遞減,在板兩端處存在較小的壓應(yīng)力,并且在溶合線以外的區(qū)域,殘余應(yīng)力隨著距離焊縫越來(lái)越遠(yuǎn),殘余應(yīng)力值也隨之迅速下降。

        (2)同一截面對(duì)接焊縫產(chǎn)生的焊接殘余應(yīng)力值大小較交錯(cuò)式結(jié)構(gòu)大,且在帶板和縱骨焊縫交界處的位置,存在應(yīng)力峰值,接近材料的屈服極限,在外板與縱骨連接處,也存在較大的應(yīng)力值,約為材料的屈服點(diǎn),并且在X、Y、Z方向均存在殘余應(yīng)力,且數(shù)值較大。

        (3)縱骨由于其結(jié)構(gòu)的特點(diǎn),產(chǎn)生了較大的焊接殘余應(yīng)力,在受載時(shí),殘余應(yīng)力與工作應(yīng)力將互相疊加,互相影響,對(duì)結(jié)構(gòu)極限強(qiáng)度有一定的影響,其中同位式對(duì)接形式的影響相比交錯(cuò)式對(duì)接而言影響更大。因此,在“一刀齊”合攏方式日趨主流的今天,研究船體結(jié)構(gòu)極限強(qiáng)度時(shí),考慮焊接殘余應(yīng)力的影響是非常必要的。

        [1]黃潔瓊,薛鴻祥,唐文勇,等.不同對(duì)接焊縫布置下艦體結(jié)構(gòu)的疲勞強(qiáng)度和極限承載分析[J].艦船科學(xué)技術(shù),2008(1):62-66.

        Huang Jieqiong,Xue Hongxiang,Tang Wenyong,et al.Fatigue and ultimate strength analysis of naval structures with different butt-welding lines distribution[J].Ship Science and Technology,2008(1):62-66.(in Chinese)

        [2]孫文婷,萬(wàn)正權(quán).對(duì)接焊殘余應(yīng)力的有限元分析[J].船舶力學(xué),2007,11(1):94-101.

        Sun Wenting,Wan Zhengquan.Finite element analysis of butt-welding residual stresses[J].Journal of Ship Mechanics,2007,11(1):94-101.(in Chinese)

        [3]張喜梅,張世聯(lián),胡鐵牛.T型材角接焊的變形和殘余應(yīng)力分析[J].艦船科學(xué)技術(shù),2008(3):72-77.

        Zhang Ximei,Zhang Shilian,Hu Tieniu.Analysis of the deformation and residual stress of T-joint fillt weld[J].Ship Science and Technology,2008(3):72-77.(in Chinese)

        [4]李良碧,潘廣善,萬(wàn)正權(quán),等.高強(qiáng)鋼錐柱結(jié)合殼焊接殘余應(yīng)力的數(shù)值模擬和試驗(yàn)研究[J].船舶力學(xué),2010,14(10):1143-1150.

        Li Liangbi,Pan Guangshan,Wan Zhengquan,et al.Numerical simulation and experiments study of welding residual stress of the cone-cylinder pressure hull of high tensile strength steel[J].Journal of Ship Mechanics,2010,14(10):1143-1150.(in Chinese)

        [5]潘廣善,王自力,李良碧,萬(wàn)正權(quán).高強(qiáng)鋼焊接結(jié)構(gòu)殘余熱應(yīng)力的有限元分析[J].船舶工程,2011(3):79-83.

        Pan Guangshan,Wang Zili,Li Liangbi,Wan Zhengquan.Numerical method of welding residual stress of the weld structure of high tensile steel[J].Ship Engineering,2011(3):79-83.(in Chinese)

        [6]陳 楚.數(shù)值分析在焊接中的應(yīng)用[M].上海:上海交通大學(xué)出版社,1985.

        Chen Chu.Numerical analysis in the application of the welding[M].Shanghai:Shanghai Jiao Tong University Press,1985.

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

        Dieter Radaj.Welding residual stress deformation heating temperature field[M].Xiong Dijing,transl.Beijing:Mechanical Industry Press,1997.

        [8]王新敏.ANSYS工程結(jié)構(gòu)數(shù)值分析[M].北京:人民交通出版社,2007.

        Wang Xinming.ANSYS numerical analysis of engineering structures[M].Beijing:People’s Traffic Press,2007.

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

        Ni Hongfang,Lin Xiang,Tu Shandong.Multi-channel welding residual stress field of 3 d finite element simulation[J].Mechanical Strength,2004,26(2):218-222.(in Chinese)

        [10]潘紅良,李培寧,李 煒.焊接殘余應(yīng)力測(cè)試方法的研究[J].壓力容器,1991(3):36-39.

        Pan Hongliang,Li Peining,Li Wei.Study of welding residual stress test method[J].Pressure Vessel Technology,1991(3):36-39.(in Chinese)

        [11]羅瑞鋒,唐文勇,宋友良.船體總段對(duì)接焊縫結(jié)構(gòu)疲勞壽命分析評(píng)估[J].艦船科學(xué)技術(shù),2009(6):27-30.

        Luo Ruifeng,Tang Wenyong,Song Youliang.Longitudinal section of the butt weld structure fatigue life analysis evaluation[J].Ship Science and Technology,2009(6):27-30.(in Chinese)

        Effect of welding residual stress under different welding procedures on ultimate strength of typical longitudinal

        LI Yong-zheng1,2,YUE Ya-lin1,WEI Peng-yu1,ZENG Qing-bo1

        (1 China Ship Scientific Research Center,Wuxi 214082,China;2 School of Naval Architecture and Ocean Engineering,Jiangsu University of Science and Technology,Zhenjiang 212003,China)

        With the development of modular technology and welding procedure,the fold way of ship block has gradually changed from the ‘step joint’ to the ‘in-plane welded connection’.In this paper,the coupling between the welding temperature field and stress field are simulated,and the distributions of the residual stress are obtained under different welding techniques.The results of the residual stress apply to the typical longitudinal to derive the ultimate strength.By analyzing the simulated results and the experiments,the distributions of residual stress are very different under different welding techniques,and the structural ultimate strength is affected.The result provides the reliable basis for follow-up study.

        in-plane welded connection;welding process;welding residual stress;numerical simulation;ultimate strength

        U661.43

        A

        10.3969/j.issn.1007-7294.2014.04.012

        1007-7294(2014)04-0441-07

        2013-12-26

        江蘇省船舶先進(jìn)設(shè)計(jì)制造技術(shù)重點(diǎn)實(shí)驗(yàn)室開放研究基金(CJ1202);江蘇省綠色船舶技術(shù)重點(diǎn)實(shí)驗(yàn)室資助

        李永正(1980-),男,中國(guó)船舶科學(xué)研究中心博士生,江蘇科技大學(xué)講師,E-mail:justyzli@163.com;

        岳亞霖(1964-),男,中國(guó)船舶科學(xué)研究中心研究員。

        猜你喜歡
        縱骨外板溫度場(chǎng)
        無(wú)頂?shù)什坌螜M艙壁附近甲板構(gòu)件布置及型式優(yōu)化
        船海工程(2021年5期)2021-10-25 08:50:34
        鋁合金加筋板焊接溫度場(chǎng)和殘余應(yīng)力數(shù)值模擬
        側(cè)圍外板轉(zhuǎn)角深拉伸起皺缺陷研究
        模具制造(2019年10期)2020-01-06 09:13:02
        側(cè)圍外板尾燈處排料困難的解決方案
        模具制造(2019年10期)2020-01-06 09:13:00
        汽車側(cè)圍外板A柱起皺缺陷分析及處理方法
        模具制造(2019年4期)2019-06-24 03:36:42
        基于紋影法的溫度場(chǎng)分布測(cè)量方法
        MJS工法與凍結(jié)法結(jié)合加固區(qū)溫度場(chǎng)研究
        建筑科技(2018年6期)2018-08-30 03:41:08
        超大型集裝箱船縱骨貫穿孔的形狀優(yōu)化
        江蘇船舶(2017年4期)2017-10-12 01:32:46
        某型快艇船體外板腐蝕點(diǎn)焊補(bǔ)工藝研究
        承受擱置腳載荷的艦船甲板縱骨設(shè)計(jì)方法
        99热免费观看| 最近高清中文在线字幕观看| 久久久精品国产亚洲av网深田| 中文字幕亚洲乱亚洲乱妇| 在线亚洲人成电影网站色www| 欧美日韩在线视频| 日本熟女视频一区二区三区| 四虎精品免费永久在线| 亚洲成av人片在线观看麦芽| 蜜桃一区二区三区视频| 国产偷闻隔壁人妻内裤av| 少妇高潮喷水正在播放| 亚洲国产精品无码中文字| 丝袜美腿在线播放一区二区| 久久久久综合一本久道| 亚洲色欲久久久久综合网| 国产综合精品一区二区三区| 最新69国产精品视频| 亚洲乱精品中文字字幕| 亚洲AV无码一区二区三区日日强| 精品国内在视频线2019| 中文字幕精品一区二区三区 | 在线一区二区三区国产精品| 成人av在线免费播放| 精品无码久久久久久久久粉色| av无码天堂一区二区三区| 国产精品久久国产精品99 gif| 精品极品一区二区三区| 亚洲国产精品一区亚洲国产| 久久综合视频网站| 亚洲中久无码永久在线观看软件| 精品www日韩熟女人妻| 色婷婷一区二区三区四区成人网| 男女动态91白浆视频| 国产黄色看三级三级三级| 亚洲国产福利成人一区二区| 久久综合第一页无码| 国产人与禽zoz0性伦| 精品无码国产一区二区三区av| 18禁在线永久免费观看| 国产av天堂亚洲av刚刚碰|