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

    BFe10白銅管材熱冷組合鑄型水平連鑄凝固溫度場(chǎng)模擬

    2012-11-03 03:32:08劉新華謝建新

    梅 俊,劉新華,謝建新

    (1. 北京科技大學(xué) 材料先進(jìn)制備技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083;2. 北京科技大學(xué) 新材料技術(shù)研究院,北京 100083)

    BFe10白銅管材熱冷組合鑄型水平連鑄凝固溫度場(chǎng)模擬

    梅 俊,劉新華,謝建新

    (1. 北京科技大學(xué) 材料先進(jìn)制備技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083;2. 北京科技大學(xué) 新材料技術(shù)研究院,北京 100083)

    建立了熱冷組合鑄型(HCCM)水平連鑄管材溫度場(chǎng)模擬模型,采用實(shí)驗(yàn)與模擬相結(jié)合的方法修正界面的換熱系數(shù)條件。所建立的HCCM水平連鑄全尺寸模擬模型和所施加邊界條件的誤差小于6%,可較好地模擬實(shí)際傳熱過程的溫度場(chǎng)。模擬結(jié)果表明:當(dāng)拉坯速度由20 mm/min增加到110 mm/min時(shí),兩相區(qū)寬度由20 mm增加至30 mm;當(dāng)熱型段加熱溫度由1 150 ℃提高到1 300 ℃時(shí),兩相區(qū)寬度由30 mm減小至12 mm;當(dāng)冷型段冷卻水流量由300 L/h增加到900 L/h時(shí),兩相區(qū)寬度由30 mm減小至20 mm;當(dāng)采用增加熱阻的改進(jìn)鑄型結(jié)構(gòu)時(shí),兩相區(qū)寬度由25 mm減小至12 mm。d 50 mm×5 mm BFe10管材HCCM水平連鑄合理的制備參數(shù)為:熔體保溫溫度1 250 ℃,連鑄拉坯速度50~80 mm/min,熱型段加熱溫度1 200~1 300 ℃,冷型段冷卻水流量500~700 L/h。關(guān)鍵詞:BFe10合金;熱冷組合鑄型;水平連鑄;溫度場(chǎng)模擬

    BFe10白銅合金具有較高的強(qiáng)度、良好的導(dǎo)熱性能和優(yōu)良的耐蝕性能,其管材作為冷凝管、熱交換器管等廣泛應(yīng)用于火電、船舶、海水淡化等領(lǐng)域[1]。白銅管材傳統(tǒng)生產(chǎn)采用“半連鑄實(shí)心鑄錠—熱擠壓管材—冷軋—拉拔”的工藝(簡(jiǎn)稱擠軋拉拔工藝),由于冷軋/拉伸道次多,加上酸洗、中間退火等工序,總加工工序多達(dá)20多個(gè)道次[2],導(dǎo)致工藝流程長(zhǎng)、能耗大、成材率低、成本高等一系列問題。

    本文作者等開發(fā)了一種白銅管材熱冷組合鑄型(Heating-cooling combined mold, 簡(jiǎn)稱HCCM)水平連鑄工藝[3]。前期研究表明,采用該工藝可連鑄出表面光亮、具有高軸向取向組織、致密度高的白銅管材。這種管材無需銑面等處理,可直接進(jìn)行大變形冷軋或拉拔成形。開發(fā)HCCM水平連鑄新工藝對(duì)于縮短傳統(tǒng)白銅管材生產(chǎn)工藝流程,提高成材率、降低生產(chǎn)成本具有重要意義。

    HCCM 水平連鑄過程控制的關(guān)鍵是通過工藝參數(shù)的合理匹配,精確控制管坯凝固的固液界面位置和溫度梯度,從而控制晶體的生長(zhǎng)方向和管坯的表面質(zhì)量。影響HCCM水平連鑄的因素較多,如熱型段加熱溫度(熱型溫度)、冷型段冷卻水流量以及拉坯速度等。對(duì)于復(fù)雜邊界條件和多因素耦合作用條件下的溫度場(chǎng)分析,數(shù)值模擬分析是非常有效的方法[4?5]。

    本文作者以BFe10白銅合金(10% Ni, 1% Fe, 1%Mn, 余量Cu,質(zhì)量分?jǐn)?shù))為對(duì)象,采用Procast軟件[6]對(duì) HCCM 水平連鑄過程鑄型內(nèi)的溫度場(chǎng)進(jìn)行模擬分析,研究制備參數(shù)及鑄型結(jié)構(gòu)對(duì)固液界面位置、兩相區(qū)寬度的影響規(guī)律,為工藝優(yōu)化與缺陷分析提供理論依據(jù)。

    圖1 熱冷組合鑄型水平連鑄工藝原理示意圖Fig. 1 Process principle schematic diagram for heating-cooling combined mold (HCCM) casting: (a) Schematic diagram of process principle; (b) Structure schematic diagram of HCCM; 1—Temperature measuring device; 2—Stopper; 3—Alloy melt; 4—Melting crucible; 5—Holding crucible; 6—Diversion pipe; 7—Mold heating device; 8—Water-cooled copper sleeve; 9—Secondary cooling water; 10—Traction device; 11—Tube; 12—Temperature measuring device; 13—Mold; 14—Core rod;Ⅰ—Section of heating mold;Ⅱ—Section of cooling mold

    1 模型與模擬方法

    1.1 HCCM水平連鑄工藝原理

    HCCM水平連鑄工藝的基本原理如圖1所示。合金在熔煉坩堝4中熔化后,經(jīng)導(dǎo)流管6進(jìn)入保溫坩堝5中保溫;當(dāng)保溫坩堝的溫度與熱型段Ⅰ的溫度達(dá)到設(shè)定值后,開啟牽引機(jī)構(gòu)10進(jìn)行管材連鑄成形。實(shí)驗(yàn)時(shí),可通過調(diào)整熔體溫度、熱型段加熱溫度、冷型段冷卻水流量和連鑄拉坯速度等,研究制備參數(shù)對(duì)管材質(zhì)量的影響。

    HCCM 水平連鑄工藝的主要特點(diǎn)是鑄型由加熱段(簡(jiǎn)稱熱型段)和冷卻段(簡(jiǎn)稱冷型段)組成。對(duì)熱型段采用加熱裝置進(jìn)行強(qiáng)制加熱,而對(duì)冷型段采用水冷銅套進(jìn)行強(qiáng)制冷卻,將管材凝固界面控制在熱型段和冷型段過渡區(qū)域內(nèi),在固液界面前沿建立起較高的軸向溫度梯度,從而可抑制晶粒沿管材徑向生長(zhǎng),促進(jìn)晶粒沿軸向生長(zhǎng),獲得沿軸向取向的結(jié)晶組織,消除普通水平連鑄(冷型連鑄)時(shí)發(fā)達(dá)的徑向柱狀晶組織。同時(shí),由于冷型段的存在,管材固液界面位置可在熱型段和冷型段過渡區(qū)域較大范圍內(nèi)變動(dòng),從而有利于避免連續(xù)定向凝固法(熱型連鑄, OCC法)[7]容易出現(xiàn)的拉漏事故,采用較大的拉鑄速度實(shí)現(xiàn)較大直徑和壁厚管材的連鑄。

    圖2 凝固溫度場(chǎng)模型示意圖Fig. 2 Schematic diagram for simulation model of solidification temperature field: (a) Simplified model;(b) Meshing and boundary conditions; 1, 2 and 3—Isothermal interface; 4—Convective heat transfer interface; 5 and 6—Radiation heat transfer interface; 7—Insulation interface

    1.2 傳熱方程及溫度場(chǎng)模型的建立

    鑄造過程溫度場(chǎng)的計(jì)算采用基于能量守恒的熱傳導(dǎo)微分方程[8?10]

    式中:ρ為密度,kg/m3;cp為定壓比熱容,J·kg/K;T為溫度,K;t是時(shí)間,s;x、y、z是空間坐標(biāo)軸;L為潛熱,J/kg;fs凝固過程中的金屬固相分?jǐn)?shù)。

    在熱傳導(dǎo)分析中,根據(jù)物體邊界處傳熱的特點(diǎn)不同,邊界條件可分為以下三類:

    1) 第一類邊界條件,給定物體邊界上任何時(shí)刻的溫度分布T=T0;

    2) 第二類邊界條件,給定物體邊界上任何時(shí)刻的熱流密度分布?nq=q0;

    3) 第三類邊界條件,給定物體邊界與周圍環(huán)境間的表面?zhèn)鳠嵯禂?shù)?nq=h(Text?T)。

    管材固液界面位于石墨鑄型內(nèi),因此,以包含鑄型、芯棒、水冷銅套和管材的一段熱冷組合鑄型為對(duì)象,建立幾何模型,分析鑄型長(zhǎng)度范圍內(nèi)管材溫度場(chǎng)。由于水平連鑄過程中管材的收縮導(dǎo)致鑄坯和結(jié)晶器間產(chǎn)生氣隙,且由于結(jié)晶器水平放置,重力作用使氣隙沿四周分布不均勻[11],導(dǎo)致管材凝固時(shí)上下傳熱存在差異,因此,在模型建立上選擇全尺寸模型,以便對(duì)這種傳熱的非對(duì)稱性進(jìn)行模擬。簡(jiǎn)化的幾何模型如圖2(a)所示。

    1.3 基本假設(shè)及邊界條件

    在模擬之前,對(duì)BFe10管材HCCM水平連鑄過程做如下簡(jiǎn)化處理:

    1) 考慮穩(wěn)態(tài)過程,即只考慮在設(shè)定的邊界條件下達(dá)到穩(wěn)定時(shí)的溫度場(chǎng);

    2) 忽略金屬液內(nèi)部的對(duì)流換熱和自然對(duì)流;

    3) 所有部件各向同性,物性參數(shù)僅與溫度有關(guān);

    4) 除結(jié)晶潛熱外,不考慮其它相變潛熱。

    圖2(b)所示為模型網(wǎng)格劃分及邊界條件。實(shí)際模擬時(shí),管材采用邊長(zhǎng)為1 mm的四邊形,鑄型、芯棒和銅套采用邊長(zhǎng)為2 mm的三角形進(jìn)行網(wǎng)格劃分。由于邊界條件是對(duì)稱施加的,因此圖2(b)中只給出了模型上半部分的邊界條件。其中邊界1、2、3均為等溫界面,屬于第一類邊界條件,邊界1、2為金屬液流入鑄型的端面,溫度等同于金屬液的溫度,邊界3為熱型段加熱溫度;邊界4上發(fā)生銅套與冷卻水之間的對(duì)流傳熱,邊界5、6上進(jìn)行輻射傳熱,均屬于第三類邊界條件,其傳熱系數(shù)通過計(jì)算確定[12?13];邊界7由于在管材內(nèi)部,可被認(rèn)為是絕熱面。

    1.4 管材與模具間界面換熱系數(shù)的處理

    由于普通水平連鑄過程管材凝固收縮和自重的影響,管材與鑄型和芯棒之間的換熱較為復(fù)雜。如圖 3所示,換熱界面包括管材外表面與鑄型頂部?jī)?nèi)表面(見圖 3(a)區(qū))、管材外表面與鑄型底部?jī)?nèi)表面(見圖 3(d)區(qū))、管材內(nèi)表面與芯棒頂部外表面(見圖 3(b)區(qū))、管材內(nèi)表面與芯棒底部?jī)?nèi)表面(見圖3(c)區(qū))等典型情況。管材在鑄型內(nèi)凝固收縮時(shí),與芯棒緊密接觸,并且由于重力作用,管材底部與鑄型同樣緊密接觸,在處理這類界面換熱系數(shù)(見圖 3(b)、(c)和(d)區(qū))時(shí),可采用實(shí)測(cè)底部石墨鑄型內(nèi)上下兩排多點(diǎn)溫度(如圖3所示),通過反問題求解法推出石墨鑄型內(nèi)側(cè)換熱邊界條件,以此作為管材與芯棒、管材底部與鑄型間界面換熱系數(shù)[14?16]。在處理管材外表面與頂部鑄型內(nèi)表面間界面換熱系數(shù)(見圖3(a)區(qū))時(shí),由于金屬液在凝固之前與鑄型接觸緊密,這時(shí)它們之間的傳熱以熱傳導(dǎo)的方式進(jìn)行,界面換熱系數(shù)較大;當(dāng)銅管發(fā)生凝固后,由于凝固收縮和自身重力的影響,使得銅管表面與石墨套之間形成了氣隙,這時(shí)它們之間的傳熱主要是通過氣隙的熱傳導(dǎo)進(jìn)行,因此,本研究設(shè)定金屬液與石墨鑄型間的換熱系數(shù)為常數(shù),金屬液凝固成管材后與鑄型間的換熱系數(shù)按氣隙傳熱計(jì)算[17]。

    圖3 管材與模具間換熱界面示意圖Fig. 3 Schematic diagram of heat transfer interface between tube and dies: Zone (a) Tube top between mold; Zone (b) Tube top between core rod; Zone (c) Tube bottom between core rod;Zone (d) Tube bottom between mold

    圖4 熱冷組合鑄型結(jié)構(gòu)示意圖及實(shí)驗(yàn)后石墨模具沿軸線剖面圖Fig. 4 Schematic diagram of mold location (a) and profile along axis after casting experiment (b)

    圖5 固液界面位置示意圖Fig. 5 Schematic diagram of liquid-solid interface position

    1.5 模擬誤差分析

    在制備參數(shù)為熔化溫度 1 250 ℃、保溫溫度 1 250 ℃、熱型段加熱溫度(熱型溫度)1 200 ℃、冷型段冷卻水流量700 L/h、拉坯速度50 mm/min條件下進(jìn)行預(yù)備實(shí)驗(yàn),實(shí)驗(yàn)結(jié)束后,將圖 4(a)所示的石墨鑄型取出,沿軸線剖切,可觀察到如圖4(b)所示的現(xiàn)象。

    從圖4(b)中劃痕初始位置以及金屬薄膜的位置,通過分別測(cè)量 L1和 L2,可確定固液界面位置,即圖4(b)中白色虛線位置,虛線左側(cè)為未凝固區(qū)(液相+固液兩相),右側(cè)為凝固區(qū)(固相)。

    將測(cè)量的L1和L2在HCCM裝配圖中標(biāo)記,即可得到圖5所示固液界面位置示意圖。由圖5可知,管材凝固的固液界面位于熱型段和冷型段之間的區(qū)域,并且固液界面與水平方向成小于90°的角度。

    借助上述確定管材固液界面位置的方法,對(duì)預(yù)備實(shí)驗(yàn)不同制備參數(shù)條件下的管材固液界面位置進(jìn)行實(shí)測(cè),將實(shí)測(cè)結(jié)果與模擬結(jié)果進(jìn)行對(duì)比,對(duì)界面換熱系數(shù)進(jìn)行修正,分析模型誤差,以確定更接近實(shí)際傳熱狀態(tài)的傳熱模型。

    圖6所示為熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 200 ℃、冷型段冷卻水流量700 L/h、拉坯速度50 mm/min條件下的模擬結(jié)果示意圖。其中圖6(a)所示為溫度分布云圖,圖6(b)所示為固相分?jǐn)?shù)圖,結(jié)合溫度分布和固相分?jǐn)?shù)圖,可以確定該制備參數(shù)條件下的固液界面位置。

    利用Procast中自帶的連鑄速度參數(shù)設(shè)置,進(jìn)行不同拉坯速度條件下溫度場(chǎng)的模擬,并將模擬結(jié)果以裝配圖中固液界面位置的形式表示,通過與實(shí)驗(yàn)獲得的數(shù)據(jù)對(duì)比,可以判斷模型的可靠性。

    圖7所示為不同拉坯速度條件下固液界面位置示意圖。由圖7可知,模擬結(jié)果與實(shí)測(cè)結(jié)果吻合較好,鑄型上部固液界面位置吻合度較高,誤差較小,鑄型下部誤差較大,最大誤差為 6%,出現(xiàn)在圖中圓圈標(biāo)記的數(shù)據(jù)位置。由此可知,本研究所建立的模型以及施加的邊界條件能夠反映實(shí)際傳熱情況,可用于溫度場(chǎng)模擬。由于管材自身重力的影響,管材上部的固液界面位置比下部更靠近鑄型出口。只要控制管材上部的固液界面位置在熱型段和冷型段交界處,則下部和中部的固液界面位置必定會(huì)在這一區(qū)間或熱型段區(qū)間,而不會(huì)位于冷型段區(qū)間。因此,以下的溫度提取及模擬結(jié)果分析主要考察管材上部的情形。

    圖6 模擬結(jié)果示意圖Fig. 6 Schematic diagram of simulation results: (a) Temperature contours; (b) Solid fraction contours

    圖7 不同拉坯速度條件下固液界面位置示意圖Fig. 7 Schematic diagram of liquid-solid interface position under different casting speeds

    1.6 模擬條件

    溫度場(chǎng)模擬時(shí)的管材尺寸為d50 mm ×5 mm。熔化溫度1 250 ℃和保溫溫度1 250 ℃固定不變。其余制備參數(shù)變化范圍如下:拉坯速度20~110 mm/min,熱型段加熱溫度1 150~1 300 ℃,冷型段冷卻水流量300~900 L/h。

    2 計(jì)算結(jié)果與討論

    2.1 制備參數(shù)對(duì)固液界面位置的影響

    采用本研究建立的模型,計(jì)算不同制備參數(shù)條件下的鑄型溫度場(chǎng),研究制備參數(shù)對(duì)管材凝固界面的影響規(guī)律。

    2.1.1 拉坯速度的影響

    在熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 200 ℃、冷型段冷卻水流量700 L/h的條件下,計(jì)算了不同拉坯速度時(shí)管材的溫度分布,結(jié)果如

    圖8所示。其中鑄型長(zhǎng)度為237 mm,右側(cè)為鑄型出口端,熱型段和冷型段的位置和寬度分別如圖8所示。圖中雙箭頭寬度為液相線與固相線之間的距離,即固液兩相區(qū)寬度;兩相區(qū)越寬,管材凝固界面前沿的溫度梯度越小,越不利于軸向熱傳導(dǎo)的建立。由于兩相區(qū)內(nèi)管材固液相共存,兩相區(qū)越寬,熱裂紋、表面粗糙等缺陷產(chǎn)生的幾率越大。為了獲得沿軸向取向顯著,表面質(zhì)量高的管材,固液界面應(yīng)控制在熱型段和冷型段之間的過渡區(qū)域內(nèi),并減小兩相區(qū)寬度,提高管材凝固界面前沿的溫度梯度,從而獲得組織沿軸向取向顯著的管材。

    圖8 不同拉坯速度管材溫度分布Fig. 8 Temperature distributions of tube under different drawing speeds in HCCM continuous casting process

    圖9 不同熱型段加熱溫度管材溫度分布Fig. 9 Temperature distributions of tube under different mold heating temperatures in HCCM continuous casting process

    由圖8可知,當(dāng)拉坯速度為20 mm/min時(shí),固液界面位置位于熱型段,兩相區(qū)寬度約為20 mm;隨著拉坯速度的提高,固液界面向鑄型出口的方向移動(dòng),兩相區(qū)寬度增加。當(dāng)拉坯速度為50 mm/min時(shí),固液界面位置位于熱型段和冷型段之間的過渡區(qū)域內(nèi),兩相區(qū)寬度約為25 mm;當(dāng)拉坯速度為80 mm/min時(shí),固液界面及整個(gè)兩相區(qū)均位于過渡區(qū)域內(nèi),兩相區(qū)寬度約為25 mm;而當(dāng)拉坯速度為110 mm/min時(shí),固液界面進(jìn)入冷型段,兩相區(qū)寬度約為30 mm。通過分析可知拉坯速度對(duì)固液界面位置的影響較大。另外隨著拉速的的提高,鑄型出口處管材溫度升高,實(shí)際實(shí)驗(yàn)過程中需要增加噴水冷卻系統(tǒng),以防止管材出鑄型后的過度氧化。

    綜上所述,在考慮固液界面位置和兩相區(qū)寬度以及制備效率的情況下,在本研究模擬條件下,合理的連鑄拉坯速度范圍為50~80 mm/min。

    2.1.2 熱型段加熱溫度的影響

    在熔化溫度1 250 ℃、保溫溫度1 250 ℃、冷型段冷卻水流量700 L/h、拉坯速度50 mm/min的條件下,計(jì)算了不同熱型段加熱溫度(熱型溫度)時(shí)管材的溫度分布,如圖9所示。由圖9可知,隨著熱型溫度由1 150 ℃升高至1 300 ℃,固液界面向鑄型出口方向移動(dòng),但均在熱型段和冷型段之間的過渡區(qū)域內(nèi),未進(jìn)入冷型段;兩相區(qū)寬度減小,由1 150 ℃時(shí)的30 mm減小到1 300 ℃時(shí)的12 mm;而鑄型出口管材溫度變化不大。由此可知,熱型溫度對(duì)固液界面位置影響較大,但對(duì)鑄型出口處管材溫度影響較?。惶岣邿嵝蜏囟龋欣诖龠M(jìn)結(jié)晶沿軸向(拉坯方向)的生長(zhǎng),提高結(jié)晶的軸向取向度。

    綜上所述,在考慮固液界面位置和兩相區(qū)寬度,以及有利于進(jìn)一步提高拉坯速度的情況下,在本模擬條件下,合理的熱型段加熱溫度為范圍為 1 200~ 1 300 ℃。

    2.1.3 冷型段冷卻水流量的影響

    在熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 200 ℃、拉坯速度50 mm/min的條件下,計(jì)算了冷型段不同冷卻水流量對(duì)管材溫度分布的影響,如圖10所示。由圖10可知,當(dāng)冷型段冷卻水流量為300 L/h時(shí),兩相區(qū)寬度約為30 mm,固液界面位于冷型段,并且鑄型出口處管材溫度較高。隨著冷卻水流量由500 L/h增大到900 L/h,固液界面在熱型段和冷型段的過渡區(qū)域向鑄型入口方向移動(dòng),但兩相區(qū)寬度變化不大,均維持在20 mm左右。冷卻水流量超過700 L/h后,固液界面位置變化不明顯。繼續(xù)增加冷卻水流量,并不能顯著增加熱型段的溫度梯度。

    綜上所述,在本研究模擬條件下,合理的冷型冷卻水流量范圍為500~700 L/h。

    圖10 不同冷型段冷卻水流量管材溫度分布Fig. 10 Temperature distributions of tube under different cooling water flow rates in HCCM continuous casting process

    圖11 不同鑄型結(jié)構(gòu)模擬結(jié)果(溫度云圖)Fig. 11 Simulation results of different mold structures (temperature contours): (a), (b) Mold model of whole structure; (c), (d) Mold model of improved structure

    2.2 模具結(jié)構(gòu)對(duì)固液界面位置的影響

    由于石墨材料具有良好的導(dǎo)熱能力,當(dāng)熱型段和冷型段石墨模具為一體,或雖然為二段結(jié)構(gòu),但其端面緊密接觸時(shí),沿石墨模具長(zhǎng)度方向產(chǎn)生大的熱傳導(dǎo),不利于建立高的軸向溫度梯度。為了提高熱型段和冷型段過渡區(qū)域內(nèi)的溫度梯度,改善鑄坯表面質(zhì)量和促進(jìn)結(jié)晶沿軸向取向生長(zhǎng),可采取措施在熱型段和冷型段之間增加熱阻,減小石墨鑄型沿軸向的傳熱,從而增大軸向溫度梯度。在熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 200 ℃、冷型段冷卻水流量700 L/h、拉坯速度50 mm/min的條件下,模擬了不同鑄型結(jié)構(gòu)對(duì)連鑄溫度場(chǎng)的影響,結(jié)果如圖11所示。兩種鑄型結(jié)構(gòu)的差別在于:整體鑄型(見圖11(a))由于石墨的良好導(dǎo)熱,熱型段和冷型段之間的熱干擾影響較大,如圖11(b)和圖12所示,兩相區(qū)寬度較大(約為25 mm),固液界面前沿溫度梯度較??;而改進(jìn)后的鑄型結(jié)構(gòu)(見圖11(c))通過在熱型段和冷型段人為增加熱阻,能夠有效減小鑄型加熱和冷卻之間的干擾,如圖11(d)和圖12可知,兩相區(qū)寬度(約12 mm)明顯減小,固液界面前沿溫度梯度增加。

    上述結(jié)果表明,在不提高熔體溫度或熱型段加熱溫度的情況下,通過合理的鑄型設(shè)計(jì)可以有效提高管材凝固界面前沿的溫度梯度,進(jìn)而改善管材的組織取向。

    根據(jù)模擬結(jié)果,當(dāng)制備參數(shù)控制在如下范圍:熔體保溫溫度 1 250 ℃,連鑄拉坯速度控制在 50~80 mm/min,熱型段加熱溫度1 200~1 300 ℃,冷型段冷卻水流量500~700 L/h,管材的凝固界面可控制在熱型段和冷型段過渡區(qū)內(nèi),HCCM水平連鑄可穩(wěn)定進(jìn)行,獲得具有良好表面質(zhì)量和沿拉坯方向的取向組織。

    圖12 不同鑄型結(jié)構(gòu)管材溫度分布Fig. 12 Temperature distributions of tube under different mold structures in HCCM continuous casting process

    圖13 連鑄管材表面的縱向條紋和橫向裂紋Fig. 13 Morphologies of cast tube: (a) Scratch under condition of melting temperature 1 250 ℃, holding temperature 1 250 ℃,heating temperature of mold 1 150 ℃, water flow of cooling mold 700 L/h, drawing speed 50 mm/min; (b) Hot cracking under condition of melting temperature 1 250 ℃, holding temperature 1 250 ℃, heating temperature of mold 1 200 ℃, water flow of cooling mold 700 L/h, drawing speed 110 mm/min

    3 實(shí)驗(yàn)驗(yàn)證

    采用整體鑄型,在熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 150 ℃、冷型段冷卻水流量700 L/h、拉坯速度50 mm/min條件下進(jìn)行了HCCM連鑄實(shí)驗(yàn),連鑄管材形貌如圖 13(a)所示。由圖 13(a)可知,管材表面出現(xiàn)了明顯的縱向條紋。根據(jù)圖8所示的模擬結(jié)果可知,雖然此時(shí)管材凝固兩相區(qū)寬度較窄,凝固界面前沿具有較大的溫度梯度,但由于整個(gè)兩相區(qū)及固液界面位置均靠近熱型段,管材凝固過程中與熱型段,凝固后與過渡段和冷型段的接觸摩擦長(zhǎng)度大,造成連鑄管材表面產(chǎn)生縱向條紋,并影響鑄型的使用壽命。

    在熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 200 ℃、冷型段冷卻水流量700 L/h、拉坯速度110 mm/min的條件下,進(jìn)行了HCCM連鑄實(shí)驗(yàn),連鑄管材形貌如圖13(b)所示。由圖13(b)可知,管材表面出現(xiàn)了橫向裂紋。根據(jù)圖8所示的模擬結(jié)果可知,此時(shí)管材的固液界面進(jìn)入了冷型段,且兩相區(qū)較寬,導(dǎo)致了熱裂紋缺陷的出現(xiàn)。

    采用改進(jìn)后的鑄型結(jié)構(gòu),在熔化溫度1 250 ℃、保溫溫度1 250 ℃、熱型段加熱溫度1 200 ℃、冷型段冷卻水流量700 L/h,拉坯速度50 mm/min條件下,制備了高質(zhì)量的白銅管材,管材形貌及金相如圖 14所示。

    由圖14(a)可知,連鑄管材的內(nèi)外表面光亮,無橘皮、微裂紋等表面缺陷,可不進(jìn)行銑面直接進(jìn)行后續(xù)軋制[3]。由圖14(b)可知,管材橫截面上壁厚均勻。由圖14(c)可知,連鑄管材組織為明顯沿軸向取向生長(zhǎng)的柱狀晶組織,枝晶主干方向與管材軸向存在一定角度,并且沿生長(zhǎng)方向不完全連續(xù),一次枝晶主干間距100~150 μm。這是由于水平連鑄上下傳熱不均勻,固液界面與水平方向成一定角度,同時(shí)枝晶沿近乎垂直于固液界面的方向生長(zhǎng),因此,枝晶與軸向會(huì)存在一定角度,并且無法沿軸向連續(xù)生長(zhǎng)。

    模擬分析和連鑄實(shí)驗(yàn)結(jié)果表明,只要將管材凝固界面位置控制在熱型段和冷型段過渡區(qū)間內(nèi),并將兩相區(qū)寬度控制在10~25 mm范圍內(nèi),就能保證連鑄管材具有沿軸向取向組織和良好的表面質(zhì)量。當(dāng)采用改進(jìn)的鑄型結(jié)構(gòu)時(shí),d 50 mm×5 mm BFe10管材合理的HCCM 水平連鑄制備參數(shù)范圍為:熔體保溫溫度 1 250 ℃,拉坯速度50~ 80 mm/min,熱型段加熱溫度1 200~ 1 300 ℃,冷型段冷卻水流量500~700 L/h。

    圖14 連鑄管材形貌及金相Fig. 14 Morphology and Metallographic of cast tube: (a) Exterior and internal surface of tube; (b) Transversal section of tube; (c)Microstructure of tube, under condition of improved mold, melting temperature 1 250 ℃, holding temperature 1 250 ℃, heating temperature of mold 1 200 ℃, water flow of cooling mold 700 L/h, drawing speed 50 mm/min

    4 結(jié)論

    1) 所建立的HCCM水平連鑄全尺寸模擬模型和所施加邊界條件的誤差小于 6%,可較好地模擬實(shí)際傳熱過程的溫度場(chǎng)。

    2) 隨著拉坯速度的提高,管材凝固界面向鑄型出口方向移動(dòng),兩相區(qū)寬度增加,由拉坯速度20 mm/min時(shí)的20 mm增加到110 mm/min時(shí)的30 mm;隨著熱型段加熱溫度的提高,管材凝固界面向鑄型出口的方向移動(dòng),兩相區(qū)寬度減小,由1 150 ℃時(shí)的30 mm減小到1 300 ℃時(shí)的12 mm;隨著冷卻水流量的增大,固液界面向鑄型入口方向移動(dòng),兩相區(qū)寬度變化不大,保持為20 mm左右。

    3) 不同鑄型結(jié)構(gòu)對(duì)HCCM水平連鑄溫度場(chǎng)有明顯的影響,采用在熱型段和冷型段之間人為增加熱阻的改進(jìn)鑄型結(jié)構(gòu),能顯著提高管材凝固界面前沿的溫度梯度,兩相區(qū)寬度由25 mm減小為12 mm。

    4) 溫度場(chǎng)模擬及實(shí)驗(yàn)結(jié)果表明:d50 mm×5 mm BFe10管材 HCCM 水平連鑄合理的制備參數(shù)范圍如下:合金熔體保溫溫度1 250 ℃,連鑄拉坯速度50~80 mm/min,熱型段加熱溫度1 200~1 300 ℃,冷型段冷卻水流量500~700 L/h。

    REFERENCES

    [1] GLOVER T J. Copper-Nickel alloy for the construction of ship and boat hulls [J]. British Corrosion Journal, 1982, 17(4): 155?157.

    [2] 李宏磊,婁花芬,馬可定. 銅加工生產(chǎn)技術(shù)問答[M]. 北京: 冶金工業(yè)出版社, 2008: 89?92.LI Hong-lei, LOU Hua-fen, MA Ke-ding. Question and answer of copper production technology [M]. Beijing: Metallurgical Industry Press, 2008: 89?92.

    [3] 謝建新, 梅 俊, 劉新華, 劉雪峰. 一種白銅管材熱冷組合鑄型水平連鑄工藝與設(shè)備[P]. 中國(guó), ZL 201010501407.4, 2012.XIE Jian-xin, MEI Jun, LIU Xin-hua, LIU Xue-feng. A kind of process and equipment for fabricating cupronickel pipes with heating-cooling combined mold casting [P]. CN, ZL 201010501407.4, 2012.

    [4] TIEU A K, KIM I S. Simulation of the continuous casting process by a mathematical model [J]. International Journal of Mechanical Sciences, 1997, 39(2): 185?192.

    [5] CHABCHOUB F, ARGYROPOULOS S A, MOSTAGHIMI J.Mathematical modeling and experimental measurements on the horizontal Ohno continuous casting process for pure tin [J].Canadian Metallurgical Quarterly, 1994, 33(1): 73?88.

    [6] 胡紅軍, 楊明波, 羅 靜, 王春歡, 陳 康. ProCAST軟件在鑄造凝固模擬中的應(yīng)用[J]. 材料科學(xué)與工藝, 2006, 14(3):293?295.HU Hong-jun, YANG Ming-bo, LUO Jing, WANG Chun-huan,CHEN Kang. Application of the software ProCAST in the casting of solidification simulation [J]. Materials Science &Technology, 2006, 14(3): 293?295.

    [7] 甘春雷, 劉雪峰, 黃海友, 謝建新. BFe10-1-1合金管材連續(xù)定向凝固制備工藝及組織力學(xué)性能[J]. 金屬學(xué)報(bào), 2010,46(12): 1549?1556.GAN Chun-lei, LIU Xue-feng, HUANG Hai-you, XIE Jian-xin.Fabrication process, microstructure and mechanical properties of BFe10-1-1 alloy tubes by continuous unidirectional solidification [J]. Acta Metallurgica Sinica, 2010, 46(12):1549?1556.

    [8] JANIK M, DYJA H. Modelling of three-dimensional temperature field inside the mould during continuous casting of steel [J]. Journal of Materials Processing Technology, 2004,157/158: 177?182.

    [9] ALVAREZ B D, MARTMHN M L, QUINTANA J M.Mathematical modeling of the process of continuous casting of aluminium and its alloys [J]. Finite Elements in Analysis and Design, 1999, 33: 43?59.

    [10] XIE S S, XIE W H, HUANG S H. Numerical simulation of temperature field of copper and copper alloy in horizontal continuous casting [J]. Rare Metals, 1999, 18(3): 195?196.

    [11] NISHIDA Y, DROSTE W, ENGLER S. The air-gap formation process at the casting-mold interface and the heat transfer mechanism through the gap [J]. Metallurgical and Materials Transactions B, 1986, 17(4): 833?844.

    [12] 侯忠霖, 姚 山, 王廷利, 張興國(guó), 金俊澤.一種鋁合金水冷界面換熱系數(shù)反求方法的研究[J]. 材料熱處理學(xué)報(bào), 2008,29(1): 157?161.HOU Zhong-lin, YAO Shan, WANG Ting-li, ZHANG Xing-guo,JIN Jun-ze. A method of inverse evaluation for interface heat transfer coefficient between aluminium alloy and cooling water[J]. Transactions of Materials and Heat Treatment, 2008, 29(1):157?161.

    [13] LEE J E, HAN H N, OH K H, YOON J K. A fully coupled analysis of fluid flow, heat transfer and stress in continuous round billet casting [J]. ISIJ International, 1999, 39(5): 435?444.

    [14] 朱憲華, 董增章, 蘇俊義. 水平連鑄鑄鐵棒材過程數(shù)值模擬[J]. 金屬學(xué)報(bào), 1990, 26(3): 171?176.ZHU Xian-hua, DONG Zeng-zhang, SU Jun-yi. Numerical simulation of solidification process of iron bar by horizontal continuous casting [J]. Acta Metallurgica Sinica, 1990, 26(3):171?176.

    [15] 郭 鵬, 張興國(guó), 郝 海, 金俊澤. AZ31鎂合金圓錠連鑄過程溫度場(chǎng)的數(shù)值模擬[J]. 中國(guó)有色金屬學(xué)報(bào), 2006, 16(9): 1570?1576.GUO Peng, ZHANG Xing-guo, HAO Hai, JIN Jun-ze.Temperature simulation of direct chill casting of AZ31 magnesium alloy billets [J]. The Chinese Journal of Nonferrous Metals, 2006, 16(9): 1570?1576.

    [16] 杜鳳山, 張 沛, 許志強(qiáng), 趙玲玲. 鋁鑄錠凝固邊界熱交換規(guī)律及溫度場(chǎng)模擬[J]. 中國(guó)有色金屬學(xué)報(bào), 2007, 17(11):1750?1754.DU Feng-shan, ZHANG Pei, XU Zhi-qiang, ZHAO Ling-ling.Law of heat transfer and simulation of temperature field for aluminum ingot solidification [J]. The Chinese Journal of Nonferrous Metals, 2007, 17(11): 1750?1754.

    [17] 張?jiān)迄i, 梁海奇, 楊秉儉, 蘇俊義, 朱憲華. 鑄鐵水平連鑄中鑄坯與石墨套界面換熱的數(shù)學(xué)模型[J]. 西安交通大學(xué)學(xué)報(bào),1998, 32(7): 90?93.ZHANG Yun-peng, LIANG Hai-qi, YANG Bing-jian, SU Jun-yi,ZHU Xian-hua. Modeling of interface heat exchange between casting bar and graphite sleeve of horizontal continuous casting[J]. Journal of Xi’an Jiao Tong University, 1998, 32(7): 90?93.

    Solidification temperature field simulation of BFe10 cupronickel tube during heating-cooling combined mold continuous casting

    MEI Jun, LIU Xin-hua, XIE Jian-xin

    (1. Key Laboratory for Advanced Materials Processing, Ministry of Education,University of Science and Technology Beijing, Beijing 100083, China;2. Institute for Advanced Materials and Technologies, University of Science and Technology Beijing,Beijing 100083, China)

    The model of solidification temperature field simulation during Heating-Cooling Combined Mold (HCCM)continuous casting was established and the heat transfer coefficient of interface was revised by the method of experiment combined with simulation. The established full-size model of HCCM continuous casting and the imposed boundary conditions could reflect the actual transfer process better with the error less than 6%. The simulation results show that with increasing drawing speed from 20 mm/min to 110 mm/min, the width of two-phase regions increases from 20 mm to 30 mm; with raising heating temperature from 1 150 ℃ to 1 300 ℃, the width of two-phase regions decreases from 30 mm to 12 mm; with increasing cooling water flow from 300 L/h to 900 L/h, the width of two-phase regions decreases from 30 mm to 20 mm; when an improved mold with a structure to add thermal resistance is adopted, the width of two-phase regions decreases from 25 mm to 12 mm. The appropriate parameters of d 50 mm×5 mm BFe10 tube by HCCM continuous casting are as follows: holding temperature of melt 1 250 ℃, drawing speed 50?80 mm/min, heating temperature of heating section 1 200?1 250 ℃, cooling water flow of cooling section 500?700 L/h.

    BFe10 alloy; heating-cooling combined mold; horizontal continuous casting; temperature field simulation

    TG249.7;TG146.1

    A

    1004-0609(2012)05-1430-10

    國(guó)家科技支撐計(jì)劃資助項(xiàng)目(2011BAE23B00)

    2011-05-29;

    2012-04-09

    謝建新,教授,博士;電話:010-62332254;E-mail: jxxie@mater.ustb.edu.cn

    (編輯 何學(xué)鋒)

    欧美日韩一级在线毛片| 嫁个100分男人电影在线观看| 久久婷婷人人爽人人干人人爱 | 日本 av在线| 狠狠狠狠99中文字幕| 欧美 亚洲 国产 日韩一| av超薄肉色丝袜交足视频| 久久国产精品男人的天堂亚洲| av超薄肉色丝袜交足视频| 少妇粗大呻吟视频| 成熟少妇高潮喷水视频| 午夜免费观看网址| 日日摸夜夜添夜夜添小说| 琪琪午夜伦伦电影理论片6080| 欧美日韩瑟瑟在线播放| 在线观看日韩欧美| 在线免费观看的www视频| 国产成人精品久久二区二区免费| 少妇裸体淫交视频免费看高清 | 嫩草影院精品99| 欧美绝顶高潮抽搐喷水| 在线天堂中文资源库| 一级a爱视频在线免费观看| 色av中文字幕| 亚洲精品粉嫩美女一区| 黄网站色视频无遮挡免费观看| 美女 人体艺术 gogo| 久久天堂一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 熟女少妇亚洲综合色aaa.| 成人三级黄色视频| 亚洲欧美日韩高清在线视频| 亚洲av片天天在线观看| 琪琪午夜伦伦电影理论片6080| 国产免费av片在线观看野外av| netflix在线观看网站| 法律面前人人平等表现在哪些方面| 亚洲七黄色美女视频| 麻豆久久精品国产亚洲av| 99精品在免费线老司机午夜| 成人手机av| 国产高清激情床上av| 国产伦一二天堂av在线观看| 18禁黄网站禁片午夜丰满| 色综合站精品国产| 午夜精品国产一区二区电影| 亚洲 欧美一区二区三区| 最近最新免费中文字幕在线| 久久婷婷人人爽人人干人人爱 | 99re在线观看精品视频| 九色国产91popny在线| 精品国产乱码久久久久久男人| 身体一侧抽搐| 午夜成年电影在线免费观看| 亚洲精品粉嫩美女一区| 不卡一级毛片| 精品卡一卡二卡四卡免费| 最新在线观看一区二区三区| 国产成人精品无人区| 精品国产美女av久久久久小说| 午夜两性在线视频| 热re99久久国产66热| 午夜成年电影在线免费观看| 日韩av在线大香蕉| 国产一区二区激情短视频| 国产不卡一卡二| 纯流量卡能插随身wifi吗| 国产成人免费无遮挡视频| 久久久精品国产亚洲av高清涩受| 一级片免费观看大全| 欧美黑人欧美精品刺激| 午夜两性在线视频| 香蕉久久夜色| 法律面前人人平等表现在哪些方面| www.www免费av| 黄色视频,在线免费观看| 久久久久久国产a免费观看| 十八禁网站免费在线| av欧美777| 亚洲片人在线观看| 国产一区在线观看成人免费| 色在线成人网| 黄色 视频免费看| 九色国产91popny在线| x7x7x7水蜜桃| 男人操女人黄网站| 久久久久久久午夜电影| 亚洲专区中文字幕在线| 亚洲av五月六月丁香网| 人妻丰满熟妇av一区二区三区| 女警被强在线播放| 老鸭窝网址在线观看| 成在线人永久免费视频| 精品一品国产午夜福利视频| 老汉色∧v一级毛片| 中文字幕精品免费在线观看视频| 人人澡人人妻人| 夜夜夜夜夜久久久久| 久久精品国产亚洲av高清一级| 精品人妻在线不人妻| 国产午夜福利久久久久久| 精品卡一卡二卡四卡免费| √禁漫天堂资源中文www| 久久精品91无色码中文字幕| www.999成人在线观看| 日本五十路高清| 夜夜看夜夜爽夜夜摸| 99国产精品99久久久久| 亚洲国产中文字幕在线视频| 日本一区二区免费在线视频| 亚洲av第一区精品v没综合| 午夜福利一区二区在线看| 人人妻人人澡欧美一区二区 | 麻豆国产av国片精品| 国产精品自产拍在线观看55亚洲| 国产真人三级小视频在线观看| 91九色精品人成在线观看| x7x7x7水蜜桃| 成人手机av| 好男人在线观看高清免费视频 | 成熟少妇高潮喷水视频| 亚洲五月天丁香| 欧美激情高清一区二区三区| 国产主播在线观看一区二区| 免费无遮挡裸体视频| 成人欧美大片| 亚洲在线自拍视频| 精品欧美一区二区三区在线| 日韩高清综合在线| 亚洲狠狠婷婷综合久久图片| av欧美777| 日本五十路高清| 黄网站色视频无遮挡免费观看| 久久精品亚洲熟妇少妇任你| 看片在线看免费视频| 美女高潮喷水抽搐中文字幕| 亚洲 欧美 日韩 在线 免费| 日本欧美视频一区| 国产精品美女特级片免费视频播放器 | www.999成人在线观看| 男人舔女人下体高潮全视频| 国产不卡一卡二| 国产亚洲精品第一综合不卡| 一区二区三区激情视频| 久久精品亚洲熟妇少妇任你| 欧美久久黑人一区二区| 女人被狂操c到高潮| 精品国产一区二区久久| 禁无遮挡网站| 国产精华一区二区三区| 夜夜爽天天搞| 精品免费久久久久久久清纯| 午夜福利免费观看在线| 亚洲精品美女久久久久99蜜臀| 在线观看免费视频网站a站| 一本大道久久a久久精品| 亚洲av片天天在线观看| 午夜激情av网站| 午夜精品国产一区二区电影| 黑人巨大精品欧美一区二区蜜桃| 午夜福利18| 纯流量卡能插随身wifi吗| 老司机深夜福利视频在线观看| 国产精品乱码一区二三区的特点 | 欧美日韩中文字幕国产精品一区二区三区 | 女性生殖器流出的白浆| 搡老妇女老女人老熟妇| 制服诱惑二区| 97超级碰碰碰精品色视频在线观看| 老司机福利观看| 人妻丰满熟妇av一区二区三区| 欧美精品啪啪一区二区三区| 久久精品91蜜桃| 18禁观看日本| 午夜福利一区二区在线看| 国产精品自产拍在线观看55亚洲| 夜夜爽天天搞| 久久久国产成人精品二区| 怎么达到女性高潮| 99国产精品99久久久久| 国产乱人伦免费视频| 天天躁夜夜躁狠狠躁躁| 久久国产精品男人的天堂亚洲| avwww免费| 国产精品一区二区三区四区久久 | 12—13女人毛片做爰片一| 日韩 欧美 亚洲 中文字幕| 在线免费观看的www视频| 性色av乱码一区二区三区2| 久久久久久久精品吃奶| 亚洲色图 男人天堂 中文字幕| 欧美成人免费av一区二区三区| 最新美女视频免费是黄的| 亚洲av片天天在线观看| 黄色成人免费大全| 久久婷婷成人综合色麻豆| 国产av在哪里看| 在线永久观看黄色视频| 成人免费观看视频高清| 午夜福利在线观看吧| 黄色片一级片一级黄色片| 老司机午夜福利在线观看视频| 91大片在线观看| 脱女人内裤的视频| av免费在线观看网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品在线美女| 搡老岳熟女国产| 99精品欧美一区二区三区四区| 亚洲aⅴ乱码一区二区在线播放 | 一边摸一边抽搐一进一出视频| 国产精品1区2区在线观看.| 精品一区二区三区视频在线观看免费| 神马国产精品三级电影在线观看 | 大码成人一级视频| 天堂影院成人在线观看| 在线观看舔阴道视频| 久久香蕉精品热| 村上凉子中文字幕在线| 久久天堂一区二区三区四区| 国产精品爽爽va在线观看网站 | 999精品在线视频| av中文乱码字幕在线| 日韩欧美三级三区| 大型av网站在线播放| 男人操女人黄网站| 性欧美人与动物交配| 99国产精品99久久久久| 国产亚洲精品久久久久久毛片| 999久久久国产精品视频| 久久精品aⅴ一区二区三区四区| 黄色视频,在线免费观看| 午夜免费激情av| 亚洲av电影不卡..在线观看| 国产99久久九九免费精品| 曰老女人黄片| 亚洲国产欧美日韩在线播放| 村上凉子中文字幕在线| 在线天堂中文资源库| 一级片免费观看大全| 桃色一区二区三区在线观看| 我的亚洲天堂| 欧美精品亚洲一区二区| av中文乱码字幕在线| 一二三四社区在线视频社区8| 免费人成视频x8x8入口观看| av片东京热男人的天堂| 久久久久精品国产欧美久久久| 又黄又粗又硬又大视频| 午夜两性在线视频| 在线观看一区二区三区| 免费看美女性在线毛片视频| 欧美人与性动交α欧美精品济南到| 中出人妻视频一区二区| 成人欧美大片| 久久国产精品影院| 日韩成人在线观看一区二区三区| 国产黄a三级三级三级人| 亚洲黑人精品在线| 很黄的视频免费| 免费观看精品视频网站| 久热爱精品视频在线9| 亚洲精品在线观看二区| 国产精品日韩av在线免费观看 | 99久久综合精品五月天人人| 91精品国产国语对白视频| 亚洲三区欧美一区| 黄色女人牲交| 亚洲七黄色美女视频| 十分钟在线观看高清视频www| 国产99久久九九免费精品| 999久久久精品免费观看国产| av视频免费观看在线观看| 最近最新中文字幕大全免费视频| 99久久国产精品久久久| 精品国产乱子伦一区二区三区| 国产亚洲av高清不卡| 91麻豆精品激情在线观看国产| 国产成人精品无人区| 女人被狂操c到高潮| 热re99久久国产66热| 纯流量卡能插随身wifi吗| 变态另类丝袜制服| 制服诱惑二区| 91老司机精品| 午夜两性在线视频| 久久婷婷人人爽人人干人人爱 | 精品卡一卡二卡四卡免费| 伊人久久大香线蕉亚洲五| 99精品欧美一区二区三区四区| 久久九九热精品免费| 九色国产91popny在线| 国产亚洲av高清不卡| 国产精品久久久人人做人人爽| 人成视频在线观看免费观看| 久久久精品国产亚洲av高清涩受| 一级毛片高清免费大全| 国产成人欧美| 中文字幕人妻丝袜一区二区| 给我免费播放毛片高清在线观看| 国产aⅴ精品一区二区三区波| 在线观看www视频免费| 日本精品一区二区三区蜜桃| 日本 欧美在线| 法律面前人人平等表现在哪些方面| 又黄又爽又免费观看的视频| 男女下面进入的视频免费午夜 | 欧美在线黄色| 亚洲一码二码三码区别大吗| 亚洲精品国产一区二区精华液| 欧美日韩瑟瑟在线播放| 国产真人三级小视频在线观看| 亚洲av熟女| 国产单亲对白刺激| 欧美性长视频在线观看| 97人妻精品一区二区三区麻豆 | 久久国产精品男人的天堂亚洲| 搡老妇女老女人老熟妇| 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av | 国产黄a三级三级三级人| 国产精品美女特级片免费视频播放器 | 国产高清有码在线观看视频 | 一级a爱视频在线免费观看| 久久国产乱子伦精品免费另类| 色综合站精品国产| 欧美日韩一级在线毛片| 精品卡一卡二卡四卡免费| 国产一区二区三区视频了| 老汉色av国产亚洲站长工具| 久9热在线精品视频| 中文字幕高清在线视频| 波多野结衣高清无吗| 国产av精品麻豆| 亚洲成国产人片在线观看| 亚洲最大成人中文| avwww免费| 亚洲av熟女| 国产精品 欧美亚洲| 国产乱人伦免费视频| av视频免费观看在线观看| 国产精品久久视频播放| 高清在线国产一区| 妹子高潮喷水视频| 国产三级黄色录像| 国产区一区二久久| 欧美 亚洲 国产 日韩一| 啦啦啦免费观看视频1| 国产精品久久视频播放| 两个人看的免费小视频| 色哟哟哟哟哟哟| 国产午夜精品久久久久久| 女人精品久久久久毛片| 国产成+人综合+亚洲专区| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 国产亚洲欧美精品永久| 成人国产综合亚洲| 午夜福利一区二区在线看| 久久精品91无色码中文字幕| 久久久久久久久免费视频了| 久久久久久人人人人人| 国产精品爽爽va在线观看网站 | 国产精品日韩av在线免费观看 | 老司机靠b影院| 国产午夜精品久久久久久| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 波多野结衣高清无吗| 欧美日韩乱码在线| 久久久久国产一级毛片高清牌| 久久精品亚洲精品国产色婷小说| 啦啦啦免费观看视频1| 国产高清有码在线观看视频 | 好看av亚洲va欧美ⅴa在| 国产av在哪里看| 黄色视频,在线免费观看| 亚洲成a人片在线一区二区| 亚洲最大成人中文| av视频免费观看在线观看| 成人手机av| 亚洲欧美日韩高清在线视频| 国产成人欧美在线观看| 国产精品 欧美亚洲| 久久国产精品影院| 国产91精品成人一区二区三区| 高清毛片免费观看视频网站| 国产在线精品亚洲第一网站| 一区二区三区高清视频在线| 久久人妻福利社区极品人妻图片| www.精华液| 麻豆久久精品国产亚洲av| av在线天堂中文字幕| 狂野欧美激情性xxxx| 精品不卡国产一区二区三区| 美女高潮到喷水免费观看| 久久久国产成人精品二区| 国产精品国产高清国产av| 在线永久观看黄色视频| 这个男人来自地球电影免费观看| 我的亚洲天堂| 婷婷精品国产亚洲av在线| 制服丝袜大香蕉在线| 级片在线观看| 91字幕亚洲| cao死你这个sao货| 亚洲 欧美一区二区三区| 国产亚洲精品第一综合不卡| 久久亚洲真实| 正在播放国产对白刺激| 动漫黄色视频在线观看| 国产欧美日韩精品亚洲av| 国产麻豆69| 久久久久精品国产欧美久久久| 不卡一级毛片| 在线永久观看黄色视频| 99在线视频只有这里精品首页| 国产亚洲精品av在线| av视频免费观看在线观看| 久久久精品国产亚洲av高清涩受| 美女国产高潮福利片在线看| av欧美777| 99在线视频只有这里精品首页| 午夜福利影视在线免费观看| 久久精品aⅴ一区二区三区四区| 免费高清在线观看日韩| av中文乱码字幕在线| 老鸭窝网址在线观看| 免费观看精品视频网站| 色在线成人网| 制服诱惑二区| 一本综合久久免费| 亚洲,欧美精品.| 很黄的视频免费| 国产精品1区2区在线观看.| 国产一区二区激情短视频| 午夜久久久在线观看| 999久久久精品免费观看国产| 免费高清在线观看日韩| 久久人妻av系列| 精品久久久久久,| 91国产中文字幕| 久久精品亚洲熟妇少妇任你| 99精品欧美一区二区三区四区| 中出人妻视频一区二区| 老汉色∧v一级毛片| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 午夜视频精品福利| 啪啪无遮挡十八禁网站| 国产亚洲精品综合一区在线观看 | 多毛熟女@视频| 日韩免费av在线播放| 欧美成人性av电影在线观看| 精品卡一卡二卡四卡免费| 国产精品免费一区二区三区在线| 99久久精品国产亚洲精品| 又黄又爽又免费观看的视频| 国产私拍福利视频在线观看| 亚洲一区高清亚洲精品| 精品一区二区三区视频在线观看免费| 免费观看人在逋| 久久热在线av| 亚洲精品一区av在线观看| 桃红色精品国产亚洲av| 亚洲国产日韩欧美精品在线观看 | 精品高清国产在线一区| 国产不卡一卡二| 亚洲人成伊人成综合网2020| 激情在线观看视频在线高清| 久久久水蜜桃国产精品网| 亚洲成国产人片在线观看| 身体一侧抽搐| 欧美日本中文国产一区发布| 国产精品美女特级片免费视频播放器 | 热99re8久久精品国产| 777久久人妻少妇嫩草av网站| 国产精品久久久久久精品电影 | 一边摸一边抽搐一进一小说| 热re99久久国产66热| 丁香欧美五月| 亚洲成av人片免费观看| 国产成+人综合+亚洲专区| 午夜免费观看网址| 久久久久久亚洲精品国产蜜桃av| 在线观看www视频免费| 黄频高清免费视频| 夜夜躁狠狠躁天天躁| 国产真人三级小视频在线观看| 欧美日韩一级在线毛片| 欧美日韩精品网址| 国产av精品麻豆| 免费一级毛片在线播放高清视频 | 日本 欧美在线| 在线国产一区二区在线| 日韩欧美一区二区三区在线观看| 亚洲欧美激情综合另类| 国产午夜福利久久久久久| 满18在线观看网站| 精品高清国产在线一区| 亚洲激情在线av| 91麻豆av在线| 国产精品影院久久| 一区二区三区精品91| 国产成人精品久久二区二区免费| 欧美午夜高清在线| 老司机福利观看| 国产av在哪里看| 久久性视频一级片| 又黄又粗又硬又大视频| 国产一区二区在线av高清观看| 欧美大码av| 亚洲视频免费观看视频| 精品久久蜜臀av无| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| 久久久久国产精品人妻aⅴ院| 99国产综合亚洲精品| 波多野结衣巨乳人妻| 国产精品久久久人人做人人爽| 日韩高清综合在线| 成人精品一区二区免费| 国产av在哪里看| 亚洲熟妇中文字幕五十中出| av在线天堂中文字幕| 国产精品 国内视频| 99久久精品国产亚洲精品| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| 免费看十八禁软件| 欧美绝顶高潮抽搐喷水| 亚洲人成77777在线视频| xxx96com| 国产私拍福利视频在线观看| 国产单亲对白刺激| 伊人久久大香线蕉亚洲五| 后天国语完整版免费观看| 一边摸一边抽搐一进一出视频| 一个人免费在线观看的高清视频| 日韩大尺度精品在线看网址 | 老司机福利观看| 亚洲伊人色综图| 在线观看日韩欧美| 黑人巨大精品欧美一区二区mp4| 国产高清有码在线观看视频 | 男女做爰动态图高潮gif福利片 | 好看av亚洲va欧美ⅴa在| 亚洲无线在线观看| 欧美最黄视频在线播放免费| 亚洲精品美女久久久久99蜜臀| 制服人妻中文乱码| 日韩 欧美 亚洲 中文字幕| 成人三级做爰电影| 欧美在线黄色| 成人国语在线视频| 搡老岳熟女国产| 欧美一级毛片孕妇| 亚洲专区中文字幕在线| 日本 欧美在线| 欧美不卡视频在线免费观看 | 精品国产乱子伦一区二区三区| 日本 av在线| 精品午夜福利视频在线观看一区| 美女大奶头视频| 久久精品国产99精品国产亚洲性色 | 午夜精品国产一区二区电影| 在线观看日韩欧美| 丝袜在线中文字幕| 后天国语完整版免费观看| 国产亚洲av高清不卡| 日韩高清综合在线| 久热爱精品视频在线9| 99riav亚洲国产免费| 97人妻天天添夜夜摸| 亚洲av熟女| 亚洲自偷自拍图片 自拍| 免费看美女性在线毛片视频| 国产精品亚洲美女久久久| 国产精品精品国产色婷婷| 精品国产国语对白av| 美女高潮喷水抽搐中文字幕| 九色国产91popny在线| 国产精品,欧美在线| 国产精品永久免费网站| 精品一品国产午夜福利视频| 欧洲精品卡2卡3卡4卡5卡区| 日韩精品免费视频一区二区三区| 成人三级做爰电影| 99re在线观看精品视频| 国产精品亚洲av一区麻豆| 亚洲国产精品999在线| 精品久久蜜臀av无| 国产又色又爽无遮挡免费看| 亚洲男人的天堂狠狠| 国产不卡一卡二| 成人18禁在线播放| 亚洲第一青青草原| 久久热在线av| 色综合亚洲欧美另类图片| 天堂影院成人在线观看| 91麻豆精品激情在线观看国产| 午夜久久久久精精品| 美女高潮喷水抽搐中文字幕| 久久热在线av| 十八禁网站免费在线| 黄片播放在线免费| 他把我摸到了高潮在线观看| 亚洲精品国产色婷婷电影| 国产成年人精品一区二区| 成人欧美大片| 国产一区二区三区在线臀色熟女| 欧美乱色亚洲激情| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 嫁个100分男人电影在线观看|