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

    矩形管輥彎成形角部充滿度分析及試驗(yàn)驗(yàn)證

    2023-12-01 07:15:08邢夢(mèng)龍杜鳳山付玉濤
    中國(guó)機(jī)械工程 2023年7期

    邢夢(mèng)龍 杜鳳山 付玉濤

    摘要:為了提高矩形管輥彎成形角部的成形精度,根據(jù)體積不變定律和金屬秒流量相等原則對(duì)矩形管角部欠充滿原因進(jìn)行分析;結(jié)合工廠實(shí)際生產(chǎn)參數(shù)建立輥彎成形有限元模型,分析圓管管徑和壁厚對(duì)成形過程中增厚系數(shù)、角部充滿度的影響;根據(jù)角部欠充滿區(qū)域中性層弧長(zhǎng)線性對(duì)應(yīng)關(guān)系進(jìn)行圓管管徑修正,確定系數(shù)修正關(guān)系式;對(duì)角部成形過程進(jìn)行損傷判定,篩選出合適的圓管管徑,將其仿真結(jié)果與試驗(yàn)進(jìn)行對(duì)比。結(jié)果表明:輥彎成形過程中壁厚的增大使得角部處于欠充滿狀態(tài);壁厚越小,成形過程中的增厚系數(shù)越小,所需圓管管徑越??;圓管管徑越大,角部充滿度越好,同時(shí)管徑的增大會(huì)增加輥彎出口處角部?jī)?nèi)側(cè)損傷的可能性;仿真與試驗(yàn)結(jié)果的誤差較小,驗(yàn)證了仿真模型的準(zhǔn)確性。

    關(guān)鍵詞:輥彎成形;角部充滿度;系數(shù)修正;角部損傷

    中圖分類號(hào):TG306

    DOI:10.3969/j.issn.1004132X.2023.07.013

    Analysis and Experimental Verification of Corner Filling Degree of

    Rectangular Tube Roll Forming

    XING Menglong1 DU Fengshan1 FU Yutao2

    1.National Cold Rolling Strip Equipment and Process Engineering Technology Research Center of

    Yanshan University,Qinhuangdao,Hebei,066004

    2.Qinhuangdao Chen-Ming Special-shaped Pipe Manufacturing Co.,Ltd.,Qinhuangdao,Hebei,066000

    Abstract: In order to improve the forming accuracy of square rectangular tube corners formed by roller forming, the causes of underfilling rectangular tube corners were analyzed according to the law of constant volume and the principle of equal metal flow per second. The finite element model of roll forming was established according to the actual production parameters in the factories, and the influences of pipe diameter and wall thickness on the thickening coefficient and corner filling degree in the forming processes were analyzed. According to the linear relation of arc length of neutral layer in underfilled regions of the corner, the tube diameter was corrected and the coefficient correction relation was determined. The damage of the corner forming processes was judged, and the appropriate diameter of the round pipe was selected. The simulation results were compared with the experimental ones. The results show that the corners are underfilled with the increase of wall thickness during roll forming. The decrease of wall thickness leads to the decrease of thickening coefficient, and the required diameters of circular pipes are smaller. With the increase of pipe diameters, the corner filling degree is better, at the same time, the increase of pipe diameters will increase the possibility of damage in the inner corner at the exit of roll forming. The errors between the test results and the experimental results are small, which verifies the accuracy of the simulation model.

    Key words: roll forming; corner filling degree; coefficient modification; corner damage

    0 引言

    輥彎成形的斷面鋼材廣泛應(yīng)用于建筑、汽車、交通運(yùn)輸、石油化工等領(lǐng)域[1]。連續(xù)輥彎成形是將板帶或管材通過輥彎?rùn)C(jī)成形輸出各種異型管的工藝。輥彎成形根據(jù)原材料幾何類型的不同,可以分為開式輥彎成形和閉式輥彎成形。其中,圓管成形方矩形管為閉式輥彎成形,優(yōu)點(diǎn)是產(chǎn)品質(zhì)量更好,機(jī)組成形速度更快,效率更高,所以得到了更廣泛的應(yīng)用。

    在閉式輥彎成形的相關(guān)研究中,NAGAMACHI等[2-3]通過仿真與試驗(yàn)研究了立輥偏置成形對(duì)方管角部充滿度的影響,并對(duì)方管管端變形進(jìn)行了分析。ABRINI等[4]和FARAHMAND等[5]根據(jù)輥彎過程中的金屬流動(dòng)規(guī)則,對(duì)管坯的變形區(qū)和自由區(qū)進(jìn)行劃分,基于上限法推導(dǎo)了輥彎過程的解析解,并結(jié)合試驗(yàn)驗(yàn)證了解析模型的準(zhǔn)確性。BAYOUMI[6-7]給出了圓管冷拔正多邊形管的解析解,研究了截面形狀、摩擦因數(shù)對(duì)力能參數(shù)的影響,并通過與試驗(yàn)結(jié)果及有限元模擬對(duì)比,驗(yàn)證了解析模型的準(zhǔn)確性,該方案可應(yīng)用于低碳鋼標(biāo)準(zhǔn)圓管。LEU等[8]基于更新的拉格朗日彈塑性有限元模型研究了多因素對(duì)圓管縮方過程中不對(duì)稱程度的影響。HOSSEINZADEH 等[9]基于有限元法(FEM)和響應(yīng)面法(RMS)對(duì)方管拉拔工藝參數(shù)進(jìn)行了優(yōu)化,通過與有限元結(jié)果對(duì)比可知,該方法能夠高度預(yù)測(cè)最優(yōu)工藝參數(shù)組合。YAO 等[10-12]基于有限元法提出了一種預(yù)測(cè)殘余應(yīng)力和等效塑性應(yīng)變的方法,通過獲得不同參數(shù)的數(shù)據(jù)建立了殘余應(yīng)力和等效塑性應(yīng)變的解析模型。彭雪鋒等[13]對(duì)高強(qiáng)鋼方管角部進(jìn)行在線感應(yīng)加熱,并對(duì)不同溫度角部?jī)?nèi)側(cè)裂紋的演化進(jìn)行了研究。韓飛等[14]對(duì)不同材料和厚度的方管進(jìn)行不同部位的取樣,通過拉伸測(cè)試發(fā)現(xiàn)R角區(qū)域的屈服強(qiáng)度大于平面區(qū)域;曾國(guó)等[15]對(duì)輥彎方矩形管采用X射線衍射的方法對(duì)其內(nèi)外表面進(jìn)行了殘余應(yīng)力的測(cè)試,發(fā)現(xiàn)采用圓成方的縱向殘余應(yīng)力水平比采用方成方的殘余應(yīng)力大;徐樹成等[16-17]采用三維大變形彈塑性有限元法,考慮材料和幾何的非線性,描述了圓管成形方矩形管的非穩(wěn)定軋制到穩(wěn)定軋制的過程;杜鳳山等[18-20]采用三維大變形理論模擬了四輥輥彎成形過程,并對(duì)管壁壁厚進(jìn)行了研究。

    盡管大多數(shù)學(xué)者都對(duì)閉式輥彎過程金屬流動(dòng)進(jìn)行了研究,但在圓管輥彎成形方矩形管過程中會(huì)出現(xiàn)管壁增厚現(xiàn)象[21],尤其是角部更加明顯。本文基于體積不變定律和金屬秒流量相等原則對(duì)圓管輥彎成形方矩形管過程中的參數(shù)變化進(jìn)行分析,建立輥彎成形有限元模型,分析圓管管徑和壁厚對(duì)輥彎過程中增厚系數(shù)和角部充滿度的影響,基于角部欠充滿狀態(tài),根據(jù)角部中性層弧長(zhǎng)對(duì)圓管管徑進(jìn)行系數(shù)修正,并對(duì)可能存在的角部損傷問題進(jìn)行損傷判定,對(duì)確定的圓管管徑進(jìn)行輥彎試驗(yàn),將仿真結(jié)果與試驗(yàn)進(jìn)行對(duì)比,驗(yàn)證了仿真模型的準(zhǔn)確性。

    1 理論方法

    1.1 體積不變定律

    在圓管輥彎成形方矩形管過程中滿足體積不變定律[22]:

    式中,K為總擴(kuò)展系數(shù);S0為圓管截面面積;Sn為矩形管截面面積;tn為方矩形管壁厚;t0為圓管壁厚;μ為延伸系數(shù);ε為增厚系數(shù)。

    在孔型設(shè)計(jì)中,ε=1,即K=μ,但在實(shí)際生產(chǎn)中,由于金屬的流動(dòng)tn會(huì)出現(xiàn)不同程度的增大,從而導(dǎo)致輥彎成形后Sn增大,容易造成角部處于欠充滿狀態(tài)。

    1.2 金屬秒流量相等原則

    在輥彎成形過程中,軋輥整體由一臺(tái)電機(jī)驅(qū)動(dòng),各個(gè)道次軋輥角速度相同,假設(shè)機(jī)架間滿足金屬秒流量相等原則[22]:

    式中,Rn為軋輥?zhàn)詈笞冃蔚来纹骄霃?;R0為軋輥第一次變形道次平均半徑;δ為平均延伸系數(shù)。

    理論上δ=μ為常數(shù),但在實(shí)際生產(chǎn)中,由于軋輥孔型上接觸點(diǎn)的變化使得質(zhì)點(diǎn)速度不同,δ會(huì)產(chǎn)生變化,δ的變化會(huì)產(chǎn)生一定的張力來保持輥彎成形的穩(wěn)定。

    1.3 輥彎成形參數(shù)變化分析

    由于在實(shí)際生產(chǎn)過程中存在管壁增厚的現(xiàn)象,同時(shí)角部容易處于欠充滿狀態(tài),故需要考慮增大圓管管徑對(duì)角部進(jìn)行補(bǔ)償填充,但這樣又會(huì)導(dǎo)致壁厚的繼續(xù)增大,即圓管管徑的增大一部分用于縱向延伸,一部分用于管壁的增厚,其中參數(shù)變化如表1所示。

    2 有限元仿真

    2.1 有限元建模

    圖1為閉式輥彎成形示意圖,其中第一道次為定徑道次不參與變形,從第二道次開始進(jìn)行壓下分配,經(jīng)過多道次輥彎成形后輸出異型管。

    對(duì)圖1過程進(jìn)行有限元建模來分析增厚系數(shù)分布及角部充滿度情況,結(jié)合工廠生產(chǎn)的40 mm×27.5 mm×3 mmSUS304不銹鋼矩形管,在MSC.MARC中根據(jù)輥花圖建立輥彎成形仿真模型,由于成形具有對(duì)稱性,所以采用1/4有限元模型,如圖2所示。

    圓管管徑根據(jù)孔型設(shè)計(jì)計(jì)算獲取。對(duì)圓管角部成形區(qū)域進(jìn)行網(wǎng)格細(xì)化,其中平輥為驅(qū)動(dòng)輥,立輥為從動(dòng)輥,各個(gè)道次角速度相同均為1 rad/s,平輥摩擦因數(shù)為0.2,立輥摩擦因數(shù)為0.001。由于第一道次不參與變形,所以模型從第二道次開始建立,第二道次變形軋輥直徑為160 mm,最后道次軋輥直徑為164 mm,中間道次中后一個(gè)道次比前一個(gè)道次按比例逐漸遞增。咬入階段由推板將圓管推入孔型,咬入后撤去推板,根據(jù)金屬秒流量相等原則推板的速度為79.7 mm/s。

    在材料模型中,彈性階段設(shè)定泊松比為0.3,彈性模量為210 GPa;塑性階段流動(dòng)應(yīng)力采用Swift模型對(duì)SUS304不銹鋼進(jìn)行擬合,如圖3所示,可以看出試驗(yàn)結(jié)果與流動(dòng)模型擬合度較高,能夠應(yīng)用于仿真模型中。

    2.2 仿真結(jié)果分析

    根據(jù)輥花圖計(jì)算圓管外半徑為20.87 mm,設(shè)計(jì)孔型與仿真結(jié)果對(duì)比如圖4所示??梢钥闯觯和獗砻骈L(zhǎng)邊與短邊重合度很高,內(nèi)表面存在一定的內(nèi)移,說明管壁的增厚是通過內(nèi)壁的內(nèi)移實(shí)現(xiàn)的,這是因?yàn)橥獗砻嬗捎谲堓伒南拗谱饔茫瑢?dǎo)致金屬只有通過向內(nèi)流動(dòng)實(shí)現(xiàn)管壁的增厚;在角部成形方面,內(nèi)側(cè)與外側(cè)均出現(xiàn)了欠充滿的現(xiàn)象,且存在一定誤差,其中金屬在內(nèi)側(cè)及外側(cè)剪應(yīng)力矢量分布如圖5所示。根據(jù)變形區(qū)(Ⅰ區(qū))可以看出,角部?jī)?nèi)側(cè)剪應(yīng)力整體高于外側(cè),內(nèi)側(cè)剪應(yīng)力在角部與長(zhǎng)短邊過渡區(qū)域附近最高,由于金屬的流動(dòng)阻力造成了金屬堆積;同時(shí)根據(jù)角部外側(cè)剪應(yīng)力分布可知,角部外側(cè)處于向角部填充狀態(tài),從而引起了角部增厚,隨著成形的結(jié)束(圖中Ⅱ區(qū)和Ⅲ區(qū))剪應(yīng)力逐漸減小,并逐漸改變方向,趨向于殘余應(yīng)力分布,即外拉內(nèi)壓[15]。根據(jù)體積不變定律可知平衡增厚系數(shù)與角部充滿度對(duì)矩形管成形有著至關(guān)重要的影響。

    2.3 圓管半徑和厚度對(duì)增厚系數(shù)的影響

    對(duì)1/4圓管網(wǎng)格節(jié)點(diǎn)進(jìn)行編號(hào),如圖6所示,分別采用圓管壁厚為1,2,3 mm和圓管外半徑R為20.87,22,22.5 mm進(jìn)行仿真模擬,在保持40 mm×27.5 mm×3 mm矩形管孔型不變的前提下,分析厚度為3 mm不同圓管管徑對(duì)增厚系數(shù)的影響。由圖7a可以看出,隨著管徑的增大,增厚系數(shù)也呈現(xiàn)增大趨勢(shì),在節(jié)點(diǎn)編號(hào)12和18附近增厚系數(shù)最大。這是因?yàn)榻遣刻幱谧杂沙尚螀^(qū),角部與長(zhǎng)邊及短邊交界區(qū)域金屬流動(dòng)劇烈,但由于管壁相對(duì)較厚,阻礙了其金屬流動(dòng),便在該區(qū)域堆積,造成增厚。分析在相同圓管管徑下不同厚度對(duì)增厚系數(shù)的影響,如圖7b所示??梢钥闯?,隨著厚度的增大,增厚系數(shù)逐漸增大,但當(dāng)t=1 mm時(shí),增厚系數(shù)變化不大,最大增厚系數(shù)出現(xiàn)在節(jié)點(diǎn)編號(hào)15附近,此時(shí)管壁厚度較小,阻礙金屬流動(dòng)的能力變?nèi)?,造成?5號(hào)節(jié)點(diǎn)附近金屬堆積,這與圖5的剪應(yīng)力矢量分布情況一致。

    2.4 圓管半徑和厚度對(duì)充滿度的影響

    圖8所示為圓管半徑和壁厚對(duì)角部充滿度的影響,可以看出,在厚度t為3 mm時(shí),采用圓管外半徑R=22.5 mm效果較好,R=20.87 mm和R=22 mm時(shí)角部處于欠充滿狀態(tài);當(dāng)厚度為2 mm時(shí),采用R=22 mm效果較好,R=22.5 mm時(shí)會(huì)使角部處于過充狀態(tài),R=20.87 mm時(shí)出現(xiàn)欠充滿;當(dāng)厚度為1 mm時(shí),采用R=22.5 mm和R=22 mm時(shí)會(huì)出現(xiàn)角部過充滿的情況,R=20.87 mm時(shí)效果較好??梢钥闯鲈谕环N孔型下,不同厚度應(yīng)匹配不同圓管管徑。

    3 線性補(bǔ)償與損傷判定

    3.1 圓管的線性補(bǔ)償

    理論上,延伸系數(shù)應(yīng)根據(jù)式(2)選取,但由于孔型上每個(gè)質(zhì)點(diǎn)的速度不同以及張力的作用,所以一般取值為1.02~1.04之間,由于輥彎過程中壁厚的增加及角部的欠充滿狀態(tài),所以需要根據(jù)角部區(qū)域中性層弧長(zhǎng)對(duì)圓管進(jìn)行系數(shù)修正,以補(bǔ)償角部欠充滿區(qū)域。

    以40 mm×27.5 mm×3 mm的矩形管孔型為例,設(shè)計(jì)圓管外半徑為20.87mm,根據(jù)圖4角部欠充滿狀態(tài)進(jìn)行補(bǔ)償計(jì)算,對(duì)設(shè)計(jì)孔型與仿真孔型進(jìn)行系數(shù)修正。充滿度采用成形后角部平均弧長(zhǎng)與設(shè)計(jì)平均弧長(zhǎng)比值作為評(píng)價(jià)指標(biāo),即

    式中,C3、C1分別為成形后角部區(qū)域內(nèi)側(cè)、外側(cè)弧長(zhǎng);C2、C4分別為設(shè)計(jì)角部區(qū)域外側(cè)、內(nèi)側(cè)弧長(zhǎng)。

    補(bǔ)償量計(jì)算如下:

    式中,yc為1/4段圓管對(duì)應(yīng)補(bǔ)償量。

    系數(shù)修正后圓管外半徑為

    式中,γ為修正系數(shù);t-n為矩形管平均壁厚;Cn為圓管中性層長(zhǎng)度。

    γ=1時(shí),計(jì)算結(jié)果如表2所示。

    3.2 角部?jī)?nèi)側(cè)損傷判定

    由有限元結(jié)果可知,在矩形管角部區(qū)域內(nèi)側(cè)出現(xiàn)了應(yīng)力集中,增加圓管管徑進(jìn)行補(bǔ)償?shù)耐瑫r(shí)會(huì)增加角部?jī)?nèi)側(cè)損傷的可能性,所以采用全解耦的韌性損傷準(zhǔn)則作為判定依據(jù)對(duì)角部區(qū)域進(jìn)行損傷傾向判定。

    應(yīng)力三軸度定義為

    式中,σm為靜水壓力;σe為等效應(yīng)力;σ1為第一主應(yīng)力,即最大主應(yīng)力;σ2為第二主應(yīng)力;σ3為第三主應(yīng)力。

    根據(jù)應(yīng)力三軸度的大小采用LeRoy準(zhǔn)則和Oh準(zhǔn)則分別進(jìn)行損傷判定。LeRoy準(zhǔn)則能夠準(zhǔn)確預(yù)測(cè)低應(yīng)力三軸度下的損傷,其定義為

    式中,εf為斷裂應(yīng)變;εp為等效塑性應(yīng)變;C為斷裂閾值。

    采用剪切缺口試樣對(duì)LeRoy準(zhǔn)則進(jìn)行標(biāo)定,損傷區(qū)域進(jìn)行局部網(wǎng)格加密,具體標(biāo)定結(jié)果如圖9所示。為了更好地說明輥彎過程中角部損傷位置,取圓管R=22.5 mm,建立輥彎損傷模型,對(duì)其角部及圓管尾部進(jìn)行網(wǎng)格加密。在MSC.MARC中采用PLOTV子程序輸出應(yīng)力三軸度和LODE角參數(shù),損傷值可以采用PLOTV或UDAMAGE_INDICATOR進(jìn)行輸出。由于兩種方法迭代計(jì)算方式不同,會(huì)產(chǎn)生一定的偏差,本文采用計(jì)算結(jié)果偏保守的PLOTV子程序,在子程序中三個(gè)主應(yīng)力需要根據(jù)應(yīng)力分量進(jìn)行轉(zhuǎn)換計(jì)算求得。損傷值計(jì)算

    結(jié)果如圖10a所示,可以看出其損傷值未達(dá)到損傷閾值,初步判定不會(huì)發(fā)生損傷。但在實(shí)際成形過程中,產(chǎn)品需要在成形截面尺寸滿足要求的情況下,確保角部不發(fā)生損傷,由于成形時(shí)角部復(fù)雜的應(yīng)力狀態(tài)、角部成形時(shí)的塑性變化以及準(zhǔn)則標(biāo)定誤差,所以需要考慮一個(gè)合理的“安全因子”進(jìn)行評(píng)價(jià)。采用工程應(yīng)變作為“安全因子”進(jìn)行損傷判定,結(jié)果如圖10b所示,可以看出輥彎結(jié)束狀態(tài)角部?jī)?nèi)側(cè)損傷值超過損傷閾值,即圖10b中A區(qū)域,此時(shí)應(yīng)力三軸度處于低水平,LODE角參數(shù)θ-在[-1,1]區(qū)間內(nèi)的分布情況如圖10b所示,可以判斷該區(qū)域?yàn)槲kU(xiǎn)區(qū)域。

    Oh準(zhǔn)則對(duì)高應(yīng)力三軸度和負(fù)應(yīng)力三軸度的預(yù)測(cè)結(jié)果較好[23],所以采用Oh準(zhǔn)則對(duì)高應(yīng)力三軸度和負(fù)應(yīng)力三軸度區(qū)域進(jìn)行損傷判定,表達(dá)式為

    采用圓角缺口試樣對(duì)Oh準(zhǔn)則進(jìn)行標(biāo)定,同樣在損傷區(qū)域進(jìn)行局部網(wǎng)格加密,標(biāo)定結(jié)果如圖11所示。取圓管R=22.5 mm進(jìn)行輥彎仿真,損傷值分布如圖12a所示,可以看出仿真結(jié)果未達(dá)到損傷閾值,初步判定不發(fā)生損傷。采用工程應(yīng)變的計(jì)算結(jié)果如圖12b所示,當(dāng)應(yīng)力三軸度小于-1/3時(shí),認(rèn)定斷裂應(yīng)變?yōu)闊o(wú)窮大,即材料永遠(yuǎn)不會(huì)發(fā)生韌性斷裂[24],可以看出在負(fù)應(yīng)力三軸度條件下,會(huì)出現(xiàn)少量的損傷區(qū)域。綜上,角部的損傷區(qū)域處于低應(yīng)力三軸度或負(fù)應(yīng)力三軸度狀態(tài)。

    由圖8可以看出,采用圓管R=22.5 mm、t=3 mm時(shí)輥彎成形角部充滿度較好,但在工程應(yīng)變的預(yù)測(cè)下,角部?jī)?nèi)側(cè)存在損傷危險(xiǎn)區(qū)域,故不能采用。采用圓管R=22 mm時(shí),提取工程應(yīng)變計(jì)算結(jié)果,如圖13所示,由角部?jī)?nèi)側(cè)損傷云圖可以看出,在低應(yīng)力三軸度下,未出現(xiàn)損傷區(qū)域,但在負(fù)應(yīng)力三軸度下,仍然存在少部分的區(qū)域超過損傷閾值,故不能采用。圖14所示為R=21.5 mm時(shí),在工程應(yīng)變條件下LeRoy準(zhǔn)則和Oh準(zhǔn)則的損傷值分布,可以看出并未超過損傷閾值。由于在R=22.5 mm時(shí)采用張量應(yīng)變計(jì)算損傷值是安全的,所以當(dāng)R=22 mm和R=21.5 mm時(shí)采用張量應(yīng)變計(jì)算也不會(huì)發(fā)生損傷,即當(dāng)R=21.5 mm時(shí),采用張量應(yīng)變和工程應(yīng)變計(jì)算都滿足損傷判據(jù),所以采用R=21.5 mm可行。同時(shí),根據(jù)GB/T6728—2002可知,采用R=22.5 mm和R=23 mm時(shí),在金屬堆積節(jié)點(diǎn)位置會(huì)出現(xiàn)壁厚超差的情況(圖7);采用R=21.5 mm時(shí)滿足壁厚偏差要求及R角要求。

    4 試驗(yàn)驗(yàn)證

    采用修正后圓管外半徑21.5 mm進(jìn)行40 mm×27.5 mm×3 mm的SUS304不銹鋼矩形管輥彎成形試驗(yàn)。將矩形管的仿真、試驗(yàn)、設(shè)計(jì)結(jié)果進(jìn)行對(duì)比,由圖15可以看出,在外表面,長(zhǎng)邊與短邊仿真結(jié)果、設(shè)計(jì)結(jié)果與試驗(yàn)結(jié)果吻合度很高;在內(nèi)表面,長(zhǎng)邊與短邊出現(xiàn)了一定程度的內(nèi)移,導(dǎo)致壁厚出現(xiàn)少量增加,仿真、設(shè)計(jì)結(jié)果與試驗(yàn)結(jié)果吻合度較高;在角部區(qū)域,外表面與內(nèi)表面均出現(xiàn)了欠充滿現(xiàn)象,同時(shí)厚度增加明顯,仿真結(jié)果與試驗(yàn)結(jié)果吻合度較高,但與設(shè)計(jì)結(jié)果存在一定差距,在可接受范圍內(nèi)。這與上述分析結(jié)果一致,從而驗(yàn)證了仿真模型的準(zhǔn)確性。輥彎圓管管徑尺寸整體確定具體流程如圖16所示。

    5 結(jié)論

    本文根據(jù)體積不變定律和金屬秒流量相等原則對(duì)圓管輥彎成形方矩形管進(jìn)行了參數(shù)分析,針對(duì)40 mm×27.5 mm×3 mm的SUS304不銹鋼矩形管建立了輥彎成形有限元模型,分析了管徑和壁厚對(duì)增厚系數(shù)和角部充滿度的影響。根據(jù)線性對(duì)應(yīng)關(guān)系對(duì)圓管管徑進(jìn)行了系數(shù)修正,并基于韌性損傷準(zhǔn)則對(duì)角部?jī)?nèi)側(cè)可能出現(xiàn)的損傷區(qū)域進(jìn)行判定,將確定的圓管仿真結(jié)果與試驗(yàn)結(jié)果進(jìn)行了對(duì)比,得到以下結(jié)論:

    (1)圓管成形方矩形管過程中由于內(nèi)側(cè)無(wú)約束存在,外側(cè)受到軋輥的限制,會(huì)使金屬沿壁厚方向產(chǎn)生流動(dòng),使得管壁厚度增加,厚度的增加會(huì)使方矩形管角部處于欠充滿狀態(tài)。

    (2)通過增加圓管管徑能夠有效緩解角部欠充滿問題,但管徑的增加會(huì)使壁厚再次增大,同時(shí)角部?jī)?nèi)側(cè)區(qū)域在輥彎成形結(jié)束階段發(fā)生金屬堆積,有引起角部損傷的可能。

    (3)在工程應(yīng)變條件下,采用R≥22 mm的圓管在輥彎結(jié)束階段會(huì)在角部?jī)?nèi)側(cè)出現(xiàn)低應(yīng)力三軸度或負(fù)應(yīng)力三軸度的損傷;采用R=21.5 mm圓管的試驗(yàn)與仿真結(jié)果吻合度較高,但與設(shè)計(jì)結(jié)果存在一定差距,在可接受范圍內(nèi),驗(yàn)證了有限元模型的準(zhǔn)確性。

    參考文獻(xiàn):

    [1] HALMOS G T. Roll Forming Handbook[M]. Boca Raton:CRC Press,2006:2-5.

    [2] NAGAMACHI T, NAKAKO T, NAKAMURA D. Effects of Forming Conditions of Roll Offset Method on Sectional Shape at the Corner of Square Steel Pipe[J]. Materials Transactions, 2013, 54:1703-1708.

    [3] NAGAMACHI T, NAKAKO T, NAKAMURA D. Mechanism of Pipe End Deformation after Cutting of Square Steel Pipe Formed by Roll Forming[J]. Materials Transactions, 2013, 54:1910-1915.

    [4] ABRINI K, FARAHMAND H R. An Upper Bound Analysis for the Reshaping of Thick Tubes with Experimental Verification[J]. International Journal of Mechanical Sciences, 2008,50:342-358.

    [5] FARAHMAND H R, ABRINI K. An Upper Bound Analysis for Reshaping Thick Tubes to Polygonal Cross-section Tubes through Multistage Roll Forming Process[J]. International Journal of Mechanical Sciences, 2015,100:90-98.

    [6] BAYOUMI L S. Cold Drawing of Regular Polygonal Tubular Sections from Round Tubes[J]. International Journal of Mechanical Sciences, 2001,43:2541-2553.

    [7] BAYOUMI L S. Analysis of Flow and Stresses in Flattening a Circular Tube by Rolling[J]. Journal of Materials Processing Technology, 2002, 128:130-135.

    [8] LEU D K, WU J Y. Finite Element Simulation of the Squaring of Circular Tube[J]. International Journal of Advanced Manufacturing Technology, 2005, 25:691-699.

    [9] HOSSEINZADEH M, MOUZIRAJI M G. An Analysis of Tube Drawing Process Used to Produce Squared Sections from Round Tubes through FE Simulation and Response Surface Methodology[J].International Journal of Advanced Manufacturing Technology, 2016,87:2179-2194.

    [10] YAO Y, QUACH W M, YOUNG B. Finite Element-based Method for Residual Stresses in Cold-formed Steel Hollow Sections[J]. Engineering Structures, 2019, 188:24-42.

    [11] YAO Y, QUACH W M, YOUNG B. Cross-section Behavior of Cold-formed Steel Elliptical Hollow Sections—a Numerical Study[J]. Engineering Structures, 2019,201:109797.

    [12] YAO Y, QUACH W M, YOUNG B. Simplified Models for Residual Stresses and Equivalent Plastic Strains in Cold-formed Steel Elliptical Hollow Sections[J]. Thin-Walled Structures,2020,154:106835.

    [13] 彭雪鋒,陳強(qiáng),王宇,等.局部感應(yīng)加熱對(duì)高強(qiáng)鋼方管二次輥壓成形的影響[J].工程科學(xué)學(xué)報(bào),2016,38(7):973-978.

    PENG Xuefeng, CHEN Qiang, WANG Yu, et al. Effect of Local Induction Heating on Secondary Roll-forming of High Strength Square Tubes[J]. Chinese Journal of Engineering, 2016,38(7):973-978.

    [14] 韓飛,張濤.中厚壁方管冷彎成形效應(yīng)研究[J]. 焊管,2016,39(7):61-64.

    HAN Fei, ZHANG Tao. Experimental Investigation on Cold-forming Effect of Middle Thick-walled Square Tube[J]. Welded Pipe and Pipe, 2016,39(7):61-64.

    [15] 曾國(guó),郭永進(jìn),于忠奇,等.輥彎成形方形型鋼的殘余應(yīng)力[J]. 鋼鐵研究學(xué)報(bào),2008,20(11):29-32.

    ZENG Guo, GUO Yongjin, YU Zhongqi, et al. Residual Stresses of Roll-formed Shape with Square Section[J]. Journal of Iron and Steel Research, 2008,20(11):29-32.

    [16] 徐樹成,劉才.方矩形管第一道次輥彎成形過程的變形和應(yīng)力場(chǎng)[J]. 燕山大學(xué)學(xué)報(bào),2000,24(1):15-18.

    XU Shucheng, LIU Cai. The Deformation and Stress Field for Square Rectangle Tube Roll Forming Process in First Pass[J]. Journal of Yanshan University, 2000,24(1):15-18.

    [17] 徐樹成,劉才.方矩形管連續(xù)輥彎成形過程的變形和應(yīng)力場(chǎng)[J]. 機(jī)械工程學(xué)報(bào),2000,36(3):93-96.

    XU Shucheng, LIU Cai. Deformation and Stress Field of Square Rectangle Tube in Continue Roll Forming Process[J]. Chinese Journal of Mechanical Engineering, 2000,36(3):93-96.

    [18] 杜鳳山,陸衛(wèi)平,張建宇,等.方管閉式輥彎成形過程的變形及流場(chǎng)[J]. 東北重型機(jī)械學(xué)院學(xué)報(bào),1997,21(2):100-103.

    DU Fengshan,LU Weiping,ZHANG Jianyu, et al. Deformation and Flow Fields of the Square Tube Forming of 4-roll Bending[J]. Journal of Northeast Heavy Machinery Institute, 1997,21(2):100-103.

    [19] 杜鳳山,付志強(qiáng),于輝.矩形管冷彎成形的輥花設(shè)計(jì)及成形規(guī)律[J]. 中國(guó)機(jī)械工程,2016,27(21):2841-2845.

    DU Fengshan,F(xiàn)U Zhiqiang,YU Hui.Roll Flower Design and Forming Mechanism of Rectangular Pipes in Cold Roll Forming Processes[J]. China Mechanical Engineering, 2016,27(21):2841-2845.

    [20] 付志強(qiáng),安子軍,杜鳳山,等.方矩形管連續(xù)輥彎成形延伸規(guī)律[J]. 吉林大學(xué)學(xué)報(bào)(工學(xué)版),2015, 45(2):487-493.

    FU Zhiqiang, AN Zijun, DU Fengshan, et al. Elongation of Rectangular Tube in Continuous Roll Forming[J]. Journal of Jilin University (Engineering and Technology Edition), 2015, 45(2):487-493.

    [21] 曹國(guó)富,盧啟威.方矩管管壁增厚的原因分析[J]. 焊管,2010,33(7):55-59.

    CAO Guofu, LU Qiwei. Discussion on Thickening Reason of Rectangular Pipe Wall Thickness[J]. Welded Pipe and Tube, 2010,33(7):55-59.

    [22] 杜鳳山,邢夢(mèng)龍,裴新元,等.一種考慮增厚系數(shù)的圓管成型方矩形管孔型設(shè)計(jì)方法:CN, 202110057360.5[P]. 2022-03-15.

    DU Fengshan,XING Menglong,PEI Xinyuan,et al. Circular Tube Forming Square Rectangular Tube Pass Design Method Based on Thickening Coefficient:CN,202110057360.5[P]. 2022-03-15.

    [23] MA H, XU W C,JIN B C, et al. Damage Evaluation in Tube Spinnability Test with Ductile Fracture Criteria[J]. International Journal of Mechanical Sciences, 2015,100:99-111.

    [24] BAO Y B, WIERZBICKI T. On the Cut-off Value of Negative Triaxiality for Fracture[J]. Engineering Fracture Mechanics, 2005, 72:1049-1069.

    国产精品久久久久成人av| 精品欧美一区二区三区在线| 亚洲黑人精品在线| 肉色欧美久久久久久久蜜桃| 黑丝袜美女国产一区| 成在线人永久免费视频| 国产精品av久久久久免费| 亚洲中文av在线| 精品一品国产午夜福利视频| 国产xxxxx性猛交| 蜜桃在线观看..| 亚洲成人国产一区在线观看| 国产一区二区三区综合在线观看| 午夜福利在线观看吧| 中文字幕色久视频| 精品高清国产在线一区| 91字幕亚洲| 欧美在线黄色| 黄色视频不卡| 69精品国产乱码久久久| 国产成人免费观看mmmm| 国精品久久久久久国模美| 婷婷成人精品国产| 国产一区二区激情短视频 | 99国产精品一区二区三区| 午夜免费成人在线视频| 欧美激情高清一区二区三区| 欧美另类亚洲清纯唯美| 老司机影院成人| 国产一区有黄有色的免费视频| 国产精品自产拍在线观看55亚洲 | 欧美精品人与动牲交sv欧美| 亚洲精品国产色婷婷电影| 欧美日韩精品网址| 无遮挡黄片免费观看| 黄色视频,在线免费观看| 国产极品粉嫩免费观看在线| 黑人操中国人逼视频| 欧美+亚洲+日韩+国产| 日韩三级视频一区二区三区| 国产麻豆69| 日日摸夜夜添夜夜添小说| 免费在线观看影片大全网站| 国产成人精品在线电影| 久久久久视频综合| 一级a爱视频在线免费观看| 国产精品香港三级国产av潘金莲| 永久免费av网站大全| 老司机福利观看| 国产老妇伦熟女老妇高清| 国产日韩欧美在线精品| 亚洲国产欧美一区二区综合| 成人18禁高潮啪啪吃奶动态图| 国产在视频线精品| 精品国内亚洲2022精品成人 | 亚洲精品久久久久久婷婷小说| 老鸭窝网址在线观看| a在线观看视频网站| 大香蕉久久成人网| 91国产中文字幕| 欧美日韩视频精品一区| 亚洲人成电影免费在线| 国产成人av教育| 欧美另类亚洲清纯唯美| 男女高潮啪啪啪动态图| 精品国产乱子伦一区二区三区 | 日韩熟女老妇一区二区性免费视频| 中文字幕人妻丝袜一区二区| 性高湖久久久久久久久免费观看| 亚洲七黄色美女视频| 操出白浆在线播放| 50天的宝宝边吃奶边哭怎么回事| 久久人人97超碰香蕉20202| 女人久久www免费人成看片| 女警被强在线播放| 老司机亚洲免费影院| 搡老岳熟女国产| 亚洲欧美精品自产自拍| 美女扒开内裤让男人捅视频| 亚洲精品国产区一区二| 黑丝袜美女国产一区| 在线av久久热| 亚洲五月色婷婷综合| 久久亚洲国产成人精品v| 国产精品成人在线| 99国产精品99久久久久| 精品一区二区三卡| 亚洲美女黄色视频免费看| 高潮久久久久久久久久久不卡| 久久久久久人人人人人| 国产伦人伦偷精品视频| 中文字幕最新亚洲高清| 丝瓜视频免费看黄片| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 亚洲五月色婷婷综合| 人人妻人人爽人人添夜夜欢视频| 日本wwww免费看| 精品国产国语对白av| 丁香六月天网| 男女边摸边吃奶| 中国国产av一级| 女人被躁到高潮嗷嗷叫费观| avwww免费| 1024视频免费在线观看| 久久国产亚洲av麻豆专区| 亚洲少妇的诱惑av| 天天操日日干夜夜撸| 久久久精品94久久精品| 91字幕亚洲| 香蕉丝袜av| 少妇的丰满在线观看| 亚洲午夜精品一区,二区,三区| 国产精品久久久人人做人人爽| a在线观看视频网站| 脱女人内裤的视频| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| netflix在线观看网站| av在线老鸭窝| 午夜福利在线观看吧| 色视频在线一区二区三区| 伊人久久大香线蕉亚洲五| 大陆偷拍与自拍| 午夜激情久久久久久久| 飞空精品影院首页| 丝袜美腿诱惑在线| 国产亚洲欧美精品永久| 国产精品成人在线| 国产免费一区二区三区四区乱码| 欧美午夜高清在线| 色婷婷av一区二区三区视频| 99热全是精品| 老汉色av国产亚洲站长工具| 日韩制服骚丝袜av| 久久人人爽人人片av| 捣出白浆h1v1| 狂野欧美激情性xxxx| 下体分泌物呈黄色| 91精品三级在线观看| 亚洲国产av新网站| 秋霞在线观看毛片| 国产av精品麻豆| 亚洲avbb在线观看| 狠狠精品人妻久久久久久综合| 黑丝袜美女国产一区| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 欧美日本中文国产一区发布| 亚洲国产欧美在线一区| 悠悠久久av| 大陆偷拍与自拍| 法律面前人人平等表现在哪些方面 | 人妻 亚洲 视频| 制服诱惑二区| 五月天丁香电影| 欧美黄色片欧美黄色片| 国产成人免费观看mmmm| 999久久久精品免费观看国产| 人人妻人人澡人人爽人人夜夜| 大片电影免费在线观看免费| 天天躁夜夜躁狠狠躁躁| 最新的欧美精品一区二区| 日日爽夜夜爽网站| 精品人妻在线不人妻| 日本av免费视频播放| 9191精品国产免费久久| 1024香蕉在线观看| 色婷婷av一区二区三区视频| 欧美日韩黄片免| 大陆偷拍与自拍| 久久久精品94久久精品| 丁香六月天网| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品日韩在线中文字幕| 女人爽到高潮嗷嗷叫在线视频| avwww免费| 一区在线观看完整版| 精品久久久精品久久久| 一二三四社区在线视频社区8| 亚洲伊人久久精品综合| 久热爱精品视频在线9| www.av在线官网国产| 少妇的丰满在线观看| 午夜两性在线视频| 精品免费久久久久久久清纯 | 久久久国产一区二区| 亚洲精品国产区一区二| 电影成人av| 久久99热这里只频精品6学生| 在线观看人妻少妇| 91国产中文字幕| av电影中文网址| a级片在线免费高清观看视频| 精品国产一区二区久久| 一二三四在线观看免费中文在| 国产一区二区激情短视频 | 捣出白浆h1v1| 色综合欧美亚洲国产小说| 超碰97精品在线观看| 久久人妻福利社区极品人妻图片| 久久综合国产亚洲精品| 日韩视频在线欧美| 亚洲欧美成人综合另类久久久| 黄色视频在线播放观看不卡| 天天添夜夜摸| 亚洲色图综合在线观看| 在线永久观看黄色视频| 欧美日韩福利视频一区二区| 欧美成人午夜精品| 亚洲七黄色美女视频| 久久久久久久精品精品| 久久精品国产亚洲av香蕉五月 | 如日韩欧美国产精品一区二区三区| 搡老熟女国产l中国老女人| 亚洲精品一二三| 99热网站在线观看| 午夜视频精品福利| 视频在线观看一区二区三区| 不卡av一区二区三区| 亚洲精品美女久久久久99蜜臀| 亚洲成国产人片在线观看| 国产亚洲一区二区精品| a级毛片在线看网站| 午夜日韩欧美国产| 捣出白浆h1v1| 国产一区二区 视频在线| tube8黄色片| 麻豆国产av国片精品| 在线十欧美十亚洲十日本专区| 中文字幕人妻丝袜一区二区| 高清在线国产一区| 欧美在线黄色| 国产97色在线日韩免费| 亚洲 国产 在线| 一级,二级,三级黄色视频| 午夜91福利影院| 亚洲国产精品999| 各种免费的搞黄视频| 另类精品久久| 中亚洲国语对白在线视频| 久久 成人 亚洲| 两个人免费观看高清视频| 国产成人一区二区三区免费视频网站| 欧美成人午夜精品| 亚洲一码二码三码区别大吗| 国产亚洲精品一区二区www | 亚洲久久久国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 美女视频免费永久观看网站| 国产成人影院久久av| 一本久久精品| 国产成人精品在线电影| 人人妻人人添人人爽欧美一区卜| 两个人免费观看高清视频| 国产区一区二久久| 黑人巨大精品欧美一区二区蜜桃| 99精品久久久久人妻精品| 又紧又爽又黄一区二区| 亚洲av欧美aⅴ国产| netflix在线观看网站| 99国产精品免费福利视频| 熟女少妇亚洲综合色aaa.| av天堂在线播放| 成年人黄色毛片网站| 欧美日韩精品网址| 国产成人精品久久二区二区91| 91精品国产国语对白视频| 久久中文看片网| 中文精品一卡2卡3卡4更新| 嫩草影视91久久| 亚洲av国产av综合av卡| 亚洲欧美精品综合一区二区三区| 最近最新中文字幕大全免费视频| 王馨瑶露胸无遮挡在线观看| 视频区欧美日本亚洲| 91大片在线观看| 国产在线观看jvid| 首页视频小说图片口味搜索| av在线app专区| 麻豆国产av国片精品| 一级黄色大片毛片| 中亚洲国语对白在线视频| 欧美人与性动交α欧美软件| 天天添夜夜摸| 啦啦啦在线免费观看视频4| 老汉色av国产亚洲站长工具| 亚洲av电影在线观看一区二区三区| 国产一区有黄有色的免费视频| 国产成人av教育| 超碰97精品在线观看| 色精品久久人妻99蜜桃| 老熟妇乱子伦视频在线观看 | 久久 成人 亚洲| 欧美黄色片欧美黄色片| 精品第一国产精品| 国产伦理片在线播放av一区| 亚洲av电影在线观看一区二区三区| 久久狼人影院| 久久久久久亚洲精品国产蜜桃av| 精品人妻1区二区| 美女主播在线视频| 韩国高清视频一区二区三区| 黄片大片在线免费观看| 一边摸一边做爽爽视频免费| 中国美女看黄片| 亚洲午夜精品一区,二区,三区| 久久中文字幕一级| 91九色精品人成在线观看| 黄色视频不卡| 国产成人免费无遮挡视频| 亚洲成av片中文字幕在线观看| 高清在线国产一区| 国产不卡av网站在线观看| 中文精品一卡2卡3卡4更新| 国产亚洲精品第一综合不卡| 桃花免费在线播放| 一个人免费看片子| 成年女人毛片免费观看观看9 | 一本综合久久免费| 99精国产麻豆久久婷婷| 欧美乱码精品一区二区三区| 国产av精品麻豆| 国产成人a∨麻豆精品| 老汉色av国产亚洲站长工具| 久久国产精品人妻蜜桃| 亚洲av美国av| 精品一区二区三区av网在线观看 | 人人妻人人澡人人爽人人夜夜| 亚洲 欧美一区二区三区| 亚洲视频免费观看视频| 丰满迷人的少妇在线观看| 久久毛片免费看一区二区三区| av网站在线播放免费| 少妇精品久久久久久久| 各种免费的搞黄视频| 国产野战对白在线观看| 狂野欧美激情性xxxx| 国产野战对白在线观看| 中文字幕色久视频| 男女无遮挡免费网站观看| 精品久久蜜臀av无| 国产真人三级小视频在线观看| 一本—道久久a久久精品蜜桃钙片| 飞空精品影院首页| bbb黄色大片| 少妇粗大呻吟视频| av福利片在线| 黑人巨大精品欧美一区二区蜜桃| 少妇人妻久久综合中文| 新久久久久国产一级毛片| 日本wwww免费看| 韩国精品一区二区三区| 99热网站在线观看| 亚洲欧美色中文字幕在线| 国产区一区二久久| 巨乳人妻的诱惑在线观看| 成人国语在线视频| 国产又爽黄色视频| 多毛熟女@视频| 考比视频在线观看| 久久国产精品人妻蜜桃| 一二三四社区在线视频社区8| 国产成人av激情在线播放| 亚洲精品一二三| 国产精品一区二区在线不卡| 美女高潮喷水抽搐中文字幕| 一区二区三区四区激情视频| 19禁男女啪啪无遮挡网站| 亚洲色图综合在线观看| 久久av网站| 精品久久蜜臀av无| 国产精品久久久久成人av| 天天躁日日躁夜夜躁夜夜| 久久狼人影院| 最黄视频免费看| 久久精品国产综合久久久| 国产精品香港三级国产av潘金莲| 中文字幕人妻熟女乱码| 777久久人妻少妇嫩草av网站| 国产精品一区二区在线观看99| √禁漫天堂资源中文www| 老熟女久久久| 欧美亚洲 丝袜 人妻 在线| 一本大道久久a久久精品| 日本黄色日本黄色录像| 大片电影免费在线观看免费| 91字幕亚洲| 国产成人影院久久av| 老司机深夜福利视频在线观看 | 成人国语在线视频| 黄色 视频免费看| 交换朋友夫妻互换小说| 99久久99久久久精品蜜桃| 精品乱码久久久久久99久播| av网站在线播放免费| 99国产极品粉嫩在线观看| 香蕉国产在线看| 97精品久久久久久久久久精品| 黄频高清免费视频| 9热在线视频观看99| 老鸭窝网址在线观看| 精品少妇一区二区三区视频日本电影| 欧美精品一区二区免费开放| 男人舔女人的私密视频| 女人精品久久久久毛片| 久久国产精品男人的天堂亚洲| 亚洲avbb在线观看| 黄片小视频在线播放| 亚洲欧美日韩高清在线视频 | 国产极品粉嫩免费观看在线| 亚洲成人免费av在线播放| 国产男人的电影天堂91| 1024香蕉在线观看| 日日爽夜夜爽网站| 少妇裸体淫交视频免费看高清 | 精品国产乱码久久久久久男人| 精品久久久精品久久久| 亚洲黑人精品在线| 亚洲精品在线美女| 亚洲欧美成人综合另类久久久| 新久久久久国产一级毛片| 欧美一级毛片孕妇| 18禁国产床啪视频网站| 啦啦啦视频在线资源免费观看| 91成年电影在线观看| 国产亚洲精品一区二区www | 狠狠婷婷综合久久久久久88av| 热99国产精品久久久久久7| 一级片免费观看大全| 久久人人爽人人片av| 成人黄色视频免费在线看| 伊人久久大香线蕉亚洲五| 在线av久久热| 少妇被粗大的猛进出69影院| 色视频在线一区二区三区| 成人免费观看视频高清| 亚洲自偷自拍图片 自拍| 母亲3免费完整高清在线观看| 国产精品国产三级国产专区5o| 后天国语完整版免费观看| 在线亚洲精品国产二区图片欧美| 日韩大片免费观看网站| 日韩欧美一区二区三区在线观看 | 欧美日韩福利视频一区二区| 午夜福利视频在线观看免费| 男女国产视频网站| 国产精品九九99| 美国免费a级毛片| 啪啪无遮挡十八禁网站| 久久久久久免费高清国产稀缺| 99久久国产精品久久久| 美女午夜性视频免费| 国产高清国产精品国产三级| 女人久久www免费人成看片| 夜夜骑夜夜射夜夜干| 久久久久久久久久久久大奶| 美女视频免费永久观看网站| 女人精品久久久久毛片| 在线永久观看黄色视频| 免费一级毛片在线播放高清视频 | 操美女的视频在线观看| 久久精品国产亚洲av香蕉五月 | 精品一区在线观看国产| 大香蕉久久网| 国产成人啪精品午夜网站| 亚洲第一av免费看| 国产高清videossex| 国产欧美亚洲国产| 欧美 亚洲 国产 日韩一| 视频在线观看一区二区三区| 天天躁夜夜躁狠狠躁躁| 大片电影免费在线观看免费| 国产精品二区激情视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲av欧美aⅴ国产| 亚洲精品国产区一区二| 国产精品亚洲av一区麻豆| 欧美黑人精品巨大| 免费在线观看日本一区| 狠狠婷婷综合久久久久久88av| 久久国产精品影院| 国产成人av教育| 青春草视频在线免费观看| 91老司机精品| 久久中文字幕一级| 欧美成人午夜精品| 亚洲伊人久久精品综合| 少妇被粗大的猛进出69影院| 午夜福利视频精品| 亚洲黑人精品在线| 纯流量卡能插随身wifi吗| 各种免费的搞黄视频| 18禁国产床啪视频网站| 99久久99久久久精品蜜桃| 老鸭窝网址在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美一级毛片孕妇| 十八禁网站免费在线| 午夜福利影视在线免费观看| 国产成人影院久久av| 丝瓜视频免费看黄片| 国产男女内射视频| 日韩三级视频一区二区三区| 岛国在线观看网站| 免费观看a级毛片全部| 国产一区有黄有色的免费视频| 丰满迷人的少妇在线观看| 老司机影院毛片| 777久久人妻少妇嫩草av网站| 午夜精品久久久久久毛片777| tube8黄色片| 精品第一国产精品| 91av网站免费观看| 丰满饥渴人妻一区二区三| 免费不卡黄色视频| 久久天堂一区二区三区四区| 午夜日韩欧美国产| 又大又爽又粗| 午夜免费成人在线视频| 制服人妻中文乱码| 精品国产乱码久久久久久男人| 麻豆乱淫一区二区| 极品人妻少妇av视频| 国产精品久久久人人做人人爽| 午夜福利视频精品| 久久天躁狠狠躁夜夜2o2o| 成人国产一区最新在线观看| 欧美另类一区| 国产成人精品在线电影| 亚洲精品成人av观看孕妇| 亚洲国产av新网站| 亚洲精品一二三| 一区二区三区乱码不卡18| 欧美日韩成人在线一区二区| 亚洲成人免费av在线播放| 两个人看的免费小视频| 黄频高清免费视频| 国产成人欧美| av网站在线播放免费| 欧美激情极品国产一区二区三区| 成年人午夜在线观看视频| 亚洲色图 男人天堂 中文字幕| 一个人免费看片子| 日本av手机在线免费观看| 国产免费福利视频在线观看| 国产成人欧美| 国产亚洲欧美精品永久| 国产亚洲欧美在线一区二区| 水蜜桃什么品种好| 极品人妻少妇av视频| 欧美老熟妇乱子伦牲交| 一本综合久久免费| 亚洲熟女精品中文字幕| 日韩精品免费视频一区二区三区| 国产成人a∨麻豆精品| 免费看十八禁软件| 日本91视频免费播放| 国产成人av激情在线播放| 男女免费视频国产| 午夜视频精品福利| 黑人猛操日本美女一级片| 午夜福利在线免费观看网站| 一边摸一边做爽爽视频免费| 亚洲成国产人片在线观看| 捣出白浆h1v1| 99久久综合免费| 亚洲欧美一区二区三区黑人| 国产免费一区二区三区四区乱码| 动漫黄色视频在线观看| 十八禁人妻一区二区| 不卡一级毛片| 国产欧美日韩综合在线一区二区| 99国产极品粉嫩在线观看| 亚洲欧洲日产国产| 国产在线视频一区二区| 国产精品久久久久成人av| a在线观看视频网站| 成年av动漫网址| 热99re8久久精品国产| 国产在线观看jvid| 国产日韩一区二区三区精品不卡| 女性生殖器流出的白浆| 亚洲午夜精品一区,二区,三区| 国产欧美日韩一区二区精品| 久久久国产欧美日韩av| 亚洲avbb在线观看| 91精品三级在线观看| 国产高清国产精品国产三级| 亚洲,欧美精品.| 久久久久久久精品精品| 19禁男女啪啪无遮挡网站| 窝窝影院91人妻| 天天添夜夜摸| 老司机在亚洲福利影院| 91国产中文字幕| 国产深夜福利视频在线观看| 欧美一级毛片孕妇| 欧美激情 高清一区二区三区| xxxhd国产人妻xxx| 久久女婷五月综合色啪小说| 欧美性长视频在线观看| 日韩大码丰满熟妇| 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| 欧美黑人精品巨大| 亚洲 欧美一区二区三区| 久久人妻熟女aⅴ| 欧美成人午夜精品| 国产国语露脸激情在线看| 亚洲第一av免费看| 欧美久久黑人一区二区| 国产一区二区三区在线臀色熟女 | 97在线人人人人妻| 制服人妻中文乱码|