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

    過(guò)冷度對(duì)Ni-Cu合金枝晶偏析影響的相場(chǎng)法模擬

    2016-08-05 00:57:04侯超杰靳玉春趙宇宏佟樂(lè)樂(lè)中北大學(xué)材料科學(xué)與工程學(xué)院太原030051
    關(guān)鍵詞:相場(chǎng)枝晶偏析

    侯超杰,靳玉春,趙宇宏,侯 華,佟樂(lè)樂(lè)(中北大學(xué) 材料科學(xué)與工程學(xué)院,太原 030051)

    過(guò)冷度對(duì)Ni-Cu合金枝晶偏析影響的相場(chǎng)法模擬

    侯超杰,靳玉春,趙宇宏,侯 華,佟樂(lè)樂(lè)
    (中北大學(xué) 材料科學(xué)與工程學(xué)院,太原 030051)

    采用耦合熱擾動(dòng)的相場(chǎng)模型,對(duì)Ni-Cu合金枝晶生長(zhǎng)中的枝晶形貌和微觀偏析進(jìn)行計(jì)算和分析,并研究過(guò)冷度對(duì)其凝固過(guò)程中溶質(zhì)分布的影響。結(jié)果表明:隨著過(guò)冷度的增加,二次枝晶變的更發(fā)達(dá),枝晶主干變細(xì),一次枝晶干軸對(duì)稱(chēng)中心處的溶質(zhì)濃度升高,一次枝晶干和二次枝晶臂間富集的溶質(zhì)更多。過(guò)冷度越大,固液界面前沿溶質(zhì)擴(kuò)散層越薄,枝晶的微觀偏析越嚴(yán)重,即溶質(zhì)微觀偏析程度隨過(guò)冷度的增大而增大。

    Ni-Cu合金;相場(chǎng)法;過(guò)冷度;枝晶偏析

    Ni-Cu合金在工業(yè)領(lǐng)域有著非常重要的應(yīng)用,其結(jié)晶方式為樹(shù)枝狀生長(zhǎng),一次枝晶臂兩側(cè)存在大量二次枝晶[1]。在其實(shí)際凝固中難以得到完全均衡態(tài)的組織,只能進(jìn)行非平衡凝固。凝固過(guò)程中溶質(zhì)元素的重新分配會(huì)造成合金元素在晶粒內(nèi)部和晶界的分布不均勻,從而導(dǎo)致成分偏析[2]。溶質(zhì)的微觀偏析會(huì)降低金屬產(chǎn)品的力學(xué)性能,因此,對(duì)金屬產(chǎn)品進(jìn)行定量預(yù)測(cè)就顯得尤為重要[3]。采用實(shí)驗(yàn)方法研究微觀偏析時(shí)由于觀察及測(cè)定困難不能準(zhǔn)確的進(jìn)行分析。隨著計(jì)算模擬技術(shù)在凝固領(lǐng)域的廣泛應(yīng)用,發(fā)展出許多凝固微觀組織的數(shù)值模擬方法。其中相場(chǎng)法[4](Phase field method)由于可以對(duì)凝固過(guò)程微觀組織的演化進(jìn)行更為貼近真實(shí)數(shù)據(jù)的模擬,已成為預(yù)測(cè)凝固微觀組織的重要手段[5-7]。以Ginzburg-Landau理論為物理基礎(chǔ)[8],通過(guò)微分方程來(lái)體現(xiàn)具有特定物理機(jī)制的擴(kuò)散、有序化勢(shì)和熱力學(xué)驅(qū)動(dòng)的綜合作用[9],再通過(guò)計(jì)算機(jī)編程求解上述方程,從而獲枝晶的凝固微觀形貌。

    金屬材料的性能不僅會(huì)受到其晶粒大小的影響,更取決于枝晶的細(xì)化程度及顯微偏析的分布。因此,掌握凝固過(guò)程的組織生長(zhǎng)是獲得理想產(chǎn)品的關(guān)鍵。本文作者采用KKS模型,通過(guò)耦合熱擾動(dòng)的相場(chǎng)模型模擬Ni-Cu二元合金枝晶生長(zhǎng)過(guò)程中過(guò)冷度對(duì)枝晶偏析的影響。對(duì)枝晶偏析情況進(jìn)行了預(yù)測(cè),可據(jù)此優(yōu)化工藝,得到力學(xué)性能更優(yōu)良的鑄件。

    1 相場(chǎng)模型及理論求解

    1.1 相場(chǎng)方程

    式中:相場(chǎng)變量Φ=0代表液態(tài),Φ=1代表固態(tài); ?t為Φ對(duì)時(shí)間求導(dǎo)數(shù);為弛豫時(shí)間;?為拉普拉斯算子;為界面能量梯度系數(shù); ?m代表分別對(duì)于x、y求偏導(dǎo)為無(wú)量綱溫度;Tm、L和cp分別是金屬熔點(diǎn)、潛熱和定壓比熱容;t為時(shí)間;D為熱擴(kuò)散系數(shù),λ為跟界面動(dòng)力學(xué)相關(guān)的參數(shù);P(Φ)為界面區(qū)域固相率;和是跟界面法相矢量相關(guān)的參數(shù)。

    1.2 溶質(zhì)擴(kuò)散方程

    在合金模擬時(shí)耦合了溶質(zhì)場(chǎng)方程,溶質(zhì)擴(kuò)散方程

    用自由能密度表示,即

    式中:D(Φ)為溶質(zhì)擴(kuò)散系數(shù);?c和?cc是自由能密度?對(duì)濃度的一階和二階偏微分;Δ為無(wú)量綱過(guò)冷度。

    1.3 溫度場(chǎng)控制方程

    在研究中,潛熱影響通過(guò)在相應(yīng)的溫度場(chǎng)網(wǎng)格中相場(chǎng)變量的變化量之和來(lái)估計(jì)[10]。這樣熱擴(kuò)散方程可表示為

    式中:D是熱擴(kuò)散系數(shù),D=λ/(ρcp)(λ為導(dǎo)熱系數(shù));T表示溫度;L表示潛熱。

    1.4 擾動(dòng)模型

    微觀上的擾動(dòng)主要源于固液兩相的熱流擾動(dòng),而界面噪聲是兩相間原子交換的結(jié)果。擾動(dòng)分為能量守恒擾動(dòng)和非能量守恒擾動(dòng),其中起決定作用的是能量守恒擾動(dòng)[11]。以下為加入能量守恒擾動(dòng)后的無(wú)量綱溫度場(chǎng)控制方程:

    式中:?tu為無(wú)量綱溫度為時(shí)間的偏導(dǎo);?tφ為相場(chǎng)對(duì)時(shí)間的偏導(dǎo)。

    式中:Tm為熔點(diǎn)溫度;L為凝固潛熱;W為界面厚度;kB為Boltzmann常數(shù);W0為界面厚度。

    2 數(shù)值求解

    2.1 初始條件和邊界條件

    假設(shè)初始晶核半徑為r0,則

    式中:x和y分別是橫坐標(biāo)和縱坐標(biāo);Δ是無(wú)量綱過(guò)冷度;C0是合金初始成分,相場(chǎng)計(jì)算區(qū)域邊界采用Zero-Neumann邊界條件[12]。

    2.2 合金熱物性參數(shù)

    模擬所用鎳銅合金物性參數(shù)如表1所列。

    表1 Ni-Cu合金物性參數(shù)[13]Table 1 Physical Parameter of Ni-Cu alloy

    2.3 數(shù)值計(jì)算方法

    為了差分方便,對(duì)于相場(chǎng)和溫度場(chǎng)的控制方程都采用均勻網(wǎng)格顯示算法,即向前Euler法。引入空間步長(zhǎng) Δ x、Δy,且 Δ x=Δy,時(shí)間步長(zhǎng) Δ t。為了獲得可靠和穩(wěn)定的計(jì)算結(jié)果,空間步長(zhǎng)Δx和時(shí)間步長(zhǎng)Δt應(yīng)滿(mǎn)足Δx≤W0和Δt≤(τ0Δx2)/(5W02)[14]。

    為了計(jì)算二次枝晶臂間距λ(z)和幅值A(chǔ)(z),首先要計(jì)算一個(gè)穩(wěn)態(tài)的枝晶主干尖端形狀x0(z),然后再計(jì)算耦合強(qiáng)度系數(shù)為λ的枝晶主干尖端形狀x0(z,t),則

    式中:z是距枝晶主干尖端的距離;v是枝晶主干尖端速度;N(z)表示在距枝晶主干尖端z處在(t2-t1)時(shí)間內(nèi)曲線(xiàn)x(z,t)極值點(diǎn)的個(gè)數(shù)。

    3 模擬結(jié)果與分析

    3.1 枝晶形貌及溶質(zhì)分布

    分別模擬過(guò)冷度系數(shù)Δ為0.195、0.295、0.395時(shí)鎳銅合金的枝晶生長(zhǎng),模擬中加入了熱噪聲擾動(dòng),并著重研究過(guò)冷度對(duì)枝晶生長(zhǎng)的影響。

    圖1所示為在1267 K下對(duì)鎳銅合金過(guò)冷熔體的凝固過(guò)程進(jìn)行模擬而獲得的相場(chǎng)分布、溶質(zhì)分布和溫度分布情況。由圖1可以看出,溶質(zhì)場(chǎng)及枝晶形貌的分布與實(shí)際情況相吻合。在基面{001}上,4條枝晶分別沿著〈2 0 0〉、〈2 0 0〉、〈02 0〉、〈0 2 0〉 4個(gè)晶向?qū)ΨQ(chēng)生長(zhǎng),一次枝晶主干間互成直角,一次枝晶與其上生長(zhǎng)的二次枝晶臂間也約成90°的夾角。

    圖 1(b1)所示為溶質(zhì)的具體分布圖。由圖 1可看出,固相中的溶質(zhì)濃度偏低,貧銅區(qū)基本集中在一次枝晶和二次枝晶臂的軸對(duì)稱(chēng)中心處。這是由于凝固時(shí)枝晶尖端曲率效應(yīng)引起過(guò)冷,使固相線(xiàn)向下移動(dòng),而固相中溶質(zhì)的擴(kuò)散速度又遠(yuǎn)遠(yuǎn)小于枝晶的生長(zhǎng)速度[15]。

    溶質(zhì)的富集主要集中在液固界面前沿,在凝固界面前沿附近出現(xiàn)了溶質(zhì)Cu的富集。其中一次枝晶根部液固界面前沿的溶質(zhì)富集程度最高,原因有以下兩點(diǎn):首先,鎳銅合金的凝固屬于非平衡結(jié)晶過(guò)程,而且合金在凝固過(guò)程中需要對(duì)溶質(zhì)進(jìn)行再分配。在非平衡凝固條件下,液相中溶質(zhì)的擴(kuò)散速度小于枝晶生長(zhǎng)速度,導(dǎo)致凝固析出的溶質(zhì)不能及時(shí)充分地?cái)U(kuò)散到液相中,從而在枝晶前沿富集。其次,一次枝晶根部的溶質(zhì)擴(kuò)散通道被二次枝晶臂包圍,溶質(zhì)擴(kuò)散較困難,嚴(yán)重阻礙了該處溶質(zhì)的擴(kuò)散,因此,該處也就成為溶質(zhì)富集程度最高的區(qū)域。

    圖1 不同過(guò)冷度系數(shù)時(shí)枝晶的相場(chǎng)、溶質(zhì)場(chǎng)和溫度場(chǎng)的形貌Fig. 1 Morphologies of dendritic ((a1), (a2), (a3)), solute field ((b1), (b2), (b3)) and thermal field ((c1), (c2), (c3)) at different supercoolings: (a1), (b1), (c1) Δ=0.195; (a2), (b2), (c2) Δ=0.295; (a3), (b3), (c3) Δ=0.395

    圖 1(a1)、(a2)、(a3)所示分別為同一時(shí)刻不同過(guò)冷度下合金凝固的枝晶形貌。由圖 1(a1)、(a2)和(a3)可知,隨著過(guò)冷度的加大,枝晶的側(cè)向分支變發(fā)達(dá),并且枝晶主干也相對(duì)要細(xì)一點(diǎn)。過(guò)冷度較小時(shí),只有少量的二次枝晶出現(xiàn),二次枝晶也比較小,且枝晶主干相對(duì)其他兩個(gè)也要更粗一點(diǎn)。這表明,過(guò)冷度的增大會(huì)促進(jìn)枝晶的生長(zhǎng),且隨著過(guò)冷度的變大,枝晶的生長(zhǎng)速度增加,二次枝晶更加發(fā)達(dá),枝晶主干變細(xì)。

    圖 1(b1)、(b2)和(b3)所示分別為不同初始溫度下枝晶的溶質(zhì)場(chǎng)分布。由圖1(b1)、(b2)和(b3)可以看出,溶質(zhì)的分布和枝晶的生長(zhǎng)保持一致。隨著過(guò)冷度的增加,枝晶尖端生長(zhǎng)速率增大,會(huì)析出更多的溶質(zhì),所以在大過(guò)冷度下時(shí),固液界面前的溶質(zhì)濃度梯度和最高溶質(zhì)濃度都要比小過(guò)冷度大。

    合金在不同初始溫度下凝固時(shí),相同時(shí)間下枝晶生長(zhǎng)的溫度場(chǎng)分布如圖 1(c1)、(c2)和(c3)所示。從圖1(c1)、(c2)和(c3)可以看出,溫度分布與枝晶生長(zhǎng)形貌相符,由于凝固潛熱的影響,液相中的溫度低于固相的,固液界面處的溫度最高。當(dāng)過(guò)冷度較大時(shí),初始溫度較低,所以枝晶生長(zhǎng)的快,釋放的潛熱更多,固相溫度升高的就更快,從而導(dǎo)致熱量來(lái)不及充分?jǐn)U散,所以最高溫度也就越大,溫度差也變的更大。

    3.2 過(guò)冷度對(duì)一次枝晶對(duì)稱(chēng)軸各處溶質(zhì)分布的影響

    圖2所示為凝固時(shí)間t=12000Δt時(shí)不同過(guò)冷度下一次枝晶對(duì)稱(chēng)軸處沿〈2 0 0〉晶向生長(zhǎng)的溶質(zhì)濃度分布曲線(xiàn)圖。曲線(xiàn)1、2、3分別對(duì)應(yīng)過(guò)冷度系數(shù)Δ為0.195、0.295、0.395時(shí),枝晶軸對(duì)稱(chēng)中心方向各處的溶質(zhì)分布曲線(xiàn)。cL1、cL2、cL3分別為 3種不同過(guò)冷下枝晶尖端前沿固液界面處液相一側(cè)的最高溶質(zhì)濃度;cS1、cS2、cS3則分別為〈2 0 0〉晶向一次枝晶干軸對(duì)稱(chēng)中心溶質(zhì)Cu的濃度。從圖2中可以看出,固相中溶質(zhì)濃度比較低。在液固界面固相一側(cè)附近,溶質(zhì)濃度急劇增大,在界面處達(dá)到最大值,在界面前沿的液相區(qū)存在一個(gè)擴(kuò)散層,在擴(kuò)散層中濃度以指數(shù)關(guān)系朝著液相初始濃度遞減,一定距離后濃度值又等于初始濃度。這是典型的由擴(kuò)散產(chǎn)生的溶質(zhì)再分配圖。

    圖2 不同過(guò)冷度下一次枝晶干軸對(duì)稱(chēng)中心沿〈2 0 0〉晶向的濃度分布Fig. 2 Concentration distribution of center of single crystal dry axis along 〈2 0 0〉 crystal orientation under different supercoolings

    由圖2可看出,隨著過(guò)冷度的不斷增加,一次枝晶軸對(duì)稱(chēng)中心處的溶質(zhì)濃度也持續(xù)提高。原因有以下兩點(diǎn):1) 隨著凝固溫度的降低,凝固點(diǎn)在相圖上的液相線(xiàn)與固相線(xiàn)向右移動(dòng),因而凝固時(shí)的平衡濃度增大,從而導(dǎo)致在凝固過(guò)程中固液界面溶質(zhì)再分配時(shí)固相中的濃度增大,過(guò)冷度越大,固相中溶質(zhì)濃度也就越高;2) 無(wú)論過(guò)冷度大小與否,液固界面處總是處于局部平衡狀態(tài)。因此,即使在非平衡結(jié)晶過(guò)程中,液固界面處固相側(cè)溶質(zhì)濃度 cS和液相側(cè)溶質(zhì)濃度 cL的比值依然滿(mǎn)足溶質(zhì)平衡擴(kuò)散系數(shù)k0,即k0=cS/cL,根據(jù)前文得出的 cL1<cL2<cL3即得到 cS1>cS1>cS1,即過(guò)冷度越大固相中溶質(zhì)濃度越高。

    從圖2中還可以看到,過(guò)冷度會(huì)影響枝晶尖端前沿的溶質(zhì)擴(kuò)散層厚度 δ。隨著過(guò)冷度增大,枝晶尖端前沿的溶質(zhì)擴(kuò)散層變薄。這是由于合金凝固時(shí),從固相中析出的溶質(zhì)需通過(guò)溶質(zhì)擴(kuò)散層釋放到液相中去,溶質(zhì)擴(kuò)散的長(zhǎng)度λs=DL/v取決于液相溶質(zhì)擴(kuò)散系數(shù)DL和枝晶生長(zhǎng)速度v的比值。由于枝晶生長(zhǎng)速度v隨過(guò)冷度的增加而增大,故熔體中溶質(zhì)擴(kuò)散的長(zhǎng)度λs隨過(guò)冷度的增大而減小,從而導(dǎo)致枝晶尖端的溶質(zhì)擴(kuò)散層厚度δ隨過(guò)冷度的增大減小。

    3.3 過(guò)冷度對(duì)合金元素顯微偏析的影響

    當(dāng)合金以樹(shù)枝狀方式凝固時(shí),由于溶質(zhì)的重新分布引起顯微偏析,從而導(dǎo)致枝晶的晶內(nèi)偏析。溶質(zhì)元素在枝晶干和枝晶間的濃度差反映出溶質(zhì)元素的偏析程度,因此通常利用枝晶間最大溶質(zhì)濃度cmax與枝晶干最小溶質(zhì)濃度cmin的比值即偏析比SR=cmax/cmin來(lái)表示微觀偏析的大小。合金元素的顯微偏析對(duì)合金的均勻化程度及力學(xué)性能具有重要影響,因此,理論和實(shí)踐方面對(duì)元素偏析進(jìn)行研究具有重要意義。

    凝固過(guò)程中的偏析取決于元素的有效分配系數(shù),在不同的過(guò)冷度條件下凝固速度、凝固溫度范圍、溶質(zhì)擴(kuò)散速率等有明顯的差別,因此不同的過(guò)冷度會(huì)導(dǎo)致不同程度的微觀偏析。圖3所示為不同過(guò)冷度下溶質(zhì)的偏析比曲線(xiàn)。從圖3中可以看出,過(guò)冷度越大,微觀偏析就越嚴(yán)重。

    圖3 不同過(guò)冷度下元素Cu的偏析比曲線(xiàn)Fig. 3 Segregation ratio curves of Cu under different supercoolings

    由圖3可知,Cu元素在晶間的濃度要高于晶干的濃度(SR>1),為正偏析。即鎳銅合金凝固時(shí)會(huì)向液相中析出溶質(zhì)Cu,而Cu在液相中的擴(kuò)散速度小于枝晶生長(zhǎng)速度,析出的溶質(zhì)擴(kuò)散不充分,從而在液固界面前沿富集,形成圖3所示的濃度分布情況。隨著過(guò)冷度增大,枝晶生長(zhǎng)速度加快,兩者成正比。隨著枝晶生長(zhǎng)速度的加快,析出到液相中的溶質(zhì)變多,二次枝晶也變的越發(fā)達(dá),對(duì)界面前沿富集溶質(zhì)擴(kuò)散的阻礙更大,溶質(zhì)越難擴(kuò)散[16]。凝固后的枝晶軸富集了高熔點(diǎn)的鎳元素,而在枝晶間區(qū)內(nèi)低熔點(diǎn)的銅元素濃度較高。最終導(dǎo)致溶質(zhì)的偏析比變大,即微觀偏析隨著過(guò)冷度的增加變得嚴(yán)重。

    4 結(jié)論

    1) 鎳銅合金凝固時(shí)高熔點(diǎn)的鎳集中于一次枝晶干和二次枝晶臂的中心,低熔點(diǎn)的溶質(zhì)銅富集于液固界面前沿。

    2) 在大過(guò)冷度下時(shí),枝晶的生長(zhǎng)速度增加,二次枝晶更加發(fā)達(dá),枝晶主干變細(xì)。

    3) 隨著過(guò)冷度增大,枝晶對(duì)稱(chēng)軸中心溶質(zhì)濃度增大,枝晶尖端前沿溶質(zhì)擴(kuò)散層厚度減小,固液界面前沿的溶質(zhì)濃度梯度和最高溶質(zhì)濃度都要比小過(guò)冷度下的大。

    4) 過(guò)冷度越大,SR越大,即微觀偏析越嚴(yán)重。

    REFERENCES

    [1]朱昌盛, 王智平, 荊 濤, 肖榮振. 二元合金微觀偏析的相場(chǎng)法數(shù)值模擬[J]. 物理學(xué)報(bào), 2006, 55(3): 1502-1507. ZHU Chang-sheng, WANG Zhi-ping, JING Tao, XIAO Rong-zhen. Numerical simulation of solute segregation patterns for a binary alloy using phase-field approach[J]. Acta Physica Sinica, 2006, 55(3): 1502-1507.

    [2]劉明翔, 王 開(kāi), 夏 丹, 蔣 濤, 朱子宗. Al-Si二元合金微觀偏析的相場(chǎng)法模擬[J]. 材料熱處理學(xué)報(bào), 2014, 35(3):213-217. LIU Ming-xiang, WANG Kai, XIA Dan, JIANG Tao, ZHU Zi-zong. Numerical simulation of the microsegregation of Al-Si binary alloy using phase-field method[J]. Transactions of Materials and Heat Treatment, 2014, 35(3): 213-217.

    [3]XIAO Rong-zhen, AN Guo-sheng, ZHU Chang-sheng, WANG Zhi-ping, YANG Shi-yin. Comparative analysis of isothermal and non-isothermal solidification of binary alloys using phase-field model[J]. Transactions of Nonferrous Metals Society of China, 2014, 24(11): 3639-3644.

    [4]陳 志, 宋慶軍, 陳安琪, 李 陽(yáng), 孔佑超, 李 峰. 相場(chǎng)法模擬凝固微觀組織演化研究進(jìn)展[J]. 鑄造技術(shù), 2011, 32(3):384-387. CHEN Zhi, SONG Qing-jun, CHEN An-qi, LI Yang, KONG You-chao, LI Feng. Recent research and progress on solidification simulated microstructures by phase field method[J]. Foundry Technology, 2011, 32(3): 384-387.

    [5]KUNDIN J, SIQUIERI R, EMMERICH H. A quantitative multi-phase-field modeling of the microstructure evolution in a peritectic Al-Ni alloy[J]. Physica D, 2013, 243(1): 116-127.

    [6]ZHANG Shu-zhou, ZHANG Rui-jie, QU Xuan-hui, FANG Wei,LIU Ming-zhi. Phase field simulation for non-isothermal solidification of multicomponent alloys coupled with thermodynamics database[J]. Transactions of Nonferrous Metals Society of China, 2013, 23(8): 2361-2367.

    [7]王明光, 趙宇宏, 任娟娜, 穆彥青, 王 偉, 楊偉明, 李?lèi)?ài)紅,葛洪浩, 侯 華. 相場(chǎng)法模擬NiCu合金非等溫凝固枝晶生長(zhǎng)[J]. 物理學(xué)報(bào), 2011, 60(4): 93-98. WANG Ming-guang, ZHAO Yu-hong, REN Juan-na, MU Yan-qing, WANG Wei, YANG Wei-ming, LI Ai-hong, GE Hong-hao, HOU Hua. Phase-field simulation of Non-Isothermal dendritic growth of NiCu alloy[J]. Acta Physica Sinica, 2011,60(4): 93-98.

    [8]龍文元, 蔡啟舟, 陳立亮, 魏伯康. 二元合金等溫凝固過(guò)程的相場(chǎng)模型[J]. 物理學(xué)報(bào), 2005, 54(1): 256-262. LONG Wen-yuan, CAI Qi-zhou, CHEN Li-liang, WEI Bo-kang. Phase-field modeling of isothermal solidification in binary alloy[J]. Acta Physica Sinica, 2005, 54(1): 256-262.

    [9]袁訓(xùn)鋒, 丁雨田. 強(qiáng)界面能各向異性下二元 Ni-Cu合金枝晶生長(zhǎng)過(guò)程的相場(chǎng)法模擬[J]. 中國(guó)有色金屬學(xué)報(bào), 2011, 21(7):1656-1663. YUAN Xun-feng, DING Yu-tian. Phase-field simulation of dendrite growth process for binary Ni-Cu alloy with anisotropy of strong interface energy[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(7): 1656-1663.

    [10]LONG Wen-yuan, LV Dong-lan, XIA Chun, PAN Mei-man, CAI Qi-zhou, CHEN Li-liang. Phase-field simulation of non-isothermal solidification dendrite growth of binary alloy under the force flow [J]. Acta Physica Sinica, 2009, 58(11):7802-7808.

    [11]張玉妥, 李東輝, 王承志, 李依依. 用相場(chǎng)方法模擬 Fe-C合金枝晶生長(zhǎng)[J]. 材料研究學(xué)報(bào), 2009, 23(3): 317-322. ZHANG Yu-tuo, LI Dong-hui, WANG Cheng-zhi, LI Yi-yi. Simulation of dendrite growth of Fe-C alloy using Phase field method[J]. Chinese Journal of Materials Research, 2009, 23(3):317-322.

    [12]CHEN Yu-juan, CHEN Chang-le. Simulation of the influence of convection velocity on upstream dendritic growth using phase-field method[J]. Acta Physica Sinica, 2008, 57(7):4585-4589.

    [13]劉 強(qiáng), 楊湘杰, 劉志凌. 強(qiáng)迫對(duì)流對(duì)合金凝固枝晶微觀偏析影響的相場(chǎng)法模擬[J]. 鑄造技術(shù), 2013, 34(12): 1698-1701. LIU Qiang, YANG Xiang-jie, LIU Zhi-ling. Simulation on effect of forced flow on microsegregation in binaryalloy based on phase field method[J]. Foundry Technology, 2013, 34(12):1698-1701.

    [14]ZHAO Yu-hong, LIU Wei-jin, HOU Hua, ZHAO Yu-hui. Impact on solidification dendrite growth by interfacial atomic motion time with phase-field method[J]. Materials Science Forum, 2013,749: 660-667.

    [15]WANG Jun-wei, WANG Zhi-ping, LU Yang, ZHU Chang-sheng,F(xiàn)ENG Li, XIAO Rong-zhen. Effect of forced lamina flow on microsegregation simulated by phase field method quantitatively[J]Transactions of Nonferrous Metals Society of China, 2012, 22(2): 391-397.

    [16]堯軍平, 李翔光, 龍文元, 張 磊, 王法軍, 李 文. 鎂合金微觀偏析相場(chǎng)法定量數(shù)值模擬[J]. 中國(guó)有色金屬學(xué)報(bào), 2014,24(1): 36-43. YAO Jun-ping, LI Xiang-guang, LONG Wen-yuan, ZHANG Lei,WANG Fa-jun, LI Wen. Numerical simulation of microsegregation of magnesium alloys quantitatively using phase-field method[J]. The Chinese Journal of Nonferrous Metals, 2014, 24(1): 36-43.

    (編輯 李艷紅)

    Numerical simulation for influence of super-cooling degree by phase-field method on microsegregation of nickel-copper alloy

    HOU Chao-jie, JIN Yu-chun, ZHAO Yu-hong, HOU Hua, TONG Le-le
    (College of Materials Science and Engineering, North University of China, Taiyuan 030051, China)

    Based on KKS model coupling the phase field and concentration field, the dendritic growth process model for Ni-Cu alloy during isothermal solidifications was developed. The solute distribution of nickel-copper alloy during the metal solidification was simulated, the effects of supercooling (ΔT) on the highest solute concentration in front of the dendrite tip and solute segregation ratio were studied emphatically. The results indicate that, the larger the supercooling ΔT is, the higher the concentration maximum value in front of the dendrite tip is, the severer the microsegregation is. The supercooling can trigger the growth of the secondary arms and increase the growth rate of the primary dendrite. As the supercooling increases, the microsegregation increases.

    Ni-Cu alloy; phase-field method; super-cooling degree; microsegregation

    Project (2014DFA50320) supported by International Cooperation of Ministry of Science and Technology of China; Projects (51574207, 51574206, 51204147, 51274175) supported by the National Natural Science Foundation of China; Projects (2013081017, 2012081013) supported by International Science and Technology Cooperation of Shanxi Province, China

    date: 2015-02-09; Accepted data: 2015-05-07

    ZHAO Yu-hong; Tel: +86-15035172958; E-mail: zyh388@sina.com

    1004-0609(2016)-01-0060-06

    TG146.1

    A

    科技部國(guó)際科技合作項(xiàng)目(2014DFA50320);國(guó)家自然科學(xué)基金資助項(xiàng)目(51574207,51574206,51204147,51274175);山西省國(guó)際科技合作項(xiàng)目(2013081017,2012081013)

    2015-02-09;

    2015-05-07

    趙宇宏,教授,博士,電話(huà):15035172958,E-mail:zyh388@sina.com

    猜你喜歡
    相場(chǎng)枝晶偏析
    基于子單元光滑有限元的混凝土相場(chǎng)損傷模型研究
    鑄件凝固微觀組織仿真程序開(kāi)發(fā)
    基于相場(chǎng)理論的瀝青自愈合微觀進(jìn)程與機(jī)理研究進(jìn)展
    石油瀝青(2018年1期)2018-04-12 07:31:51
    基于COMSOL的相場(chǎng)模擬研究
    科技視界(2017年8期)2017-07-31 10:31:17
    A356合金近液相線(xiàn)半固態(tài)鑄造非枝晶組織與模鍛成形性能
    改善60Si2MnA彈簧鋼小方坯中心碳偏析的研究
    上海金屬(2016年4期)2016-11-23 05:39:12
    原鋁旋轉(zhuǎn)偏析法凈化工藝
    改善82B 小方坯碳偏析的連鑄工藝研究
    新疆鋼鐵(2015年1期)2015-11-07 03:11:03
    SWRH82B小方坯中心偏析的改善
    上海金屬(2014年1期)2014-12-18 06:52:12
    不同形狀橫向限制對(duì)枝晶間距影響作用的相場(chǎng)法模擬
    美女大奶头视频| 国产乱人伦免费视频| 中文亚洲av片在线观看爽| 大型黄色视频在线免费观看| 91字幕亚洲| 亚洲一区二区三区色噜噜| 国产亚洲精品久久久久久毛片| 黄色成人免费大全| 精品国内亚洲2022精品成人| 中出人妻视频一区二区| 女人被狂操c到高潮| 亚洲成国产人片在线观看| 可以免费在线观看a视频的电影网站| 欧美中文日本在线观看视频| 久久久国产成人精品二区| 亚洲欧美一区二区三区黑人| 中亚洲国语对白在线视频| 久久天躁狠狠躁夜夜2o2o| av网站免费在线观看视频| 亚洲欧美一区二区三区黑人| 免费看十八禁软件| 久久精品91无色码中文字幕| 午夜精品在线福利| 久久 成人 亚洲| 侵犯人妻中文字幕一二三四区| 无遮挡黄片免费观看| 露出奶头的视频| 别揉我奶头~嗯~啊~动态视频| 天堂影院成人在线观看| 精品国产亚洲在线| 亚洲av美国av| 麻豆国产av国片精品| 多毛熟女@视频| 欧美绝顶高潮抽搐喷水| 美女大奶头视频| 天天添夜夜摸| 欧美黄色淫秽网站| 欧美日本视频| 村上凉子中文字幕在线| 黄频高清免费视频| 黄片大片在线免费观看| 国产av又大| 别揉我奶头~嗯~啊~动态视频| 极品人妻少妇av视频| 精品熟女少妇八av免费久了| 午夜老司机福利片| 精品久久蜜臀av无| 两个人看的免费小视频| 久久久久久久久中文| netflix在线观看网站| 亚洲专区中文字幕在线| 午夜影院日韩av| 少妇被粗大的猛进出69影院| 亚洲av美国av| 免费人成视频x8x8入口观看| 女人精品久久久久毛片| 国产高清videossex| 日日摸夜夜添夜夜添小说| 欧美 亚洲 国产 日韩一| 亚洲最大成人中文| 久久性视频一级片| 操美女的视频在线观看| 欧美中文日本在线观看视频| 人妻久久中文字幕网| 在线观看午夜福利视频| tocl精华| 美女 人体艺术 gogo| 精品午夜福利视频在线观看一区| 少妇裸体淫交视频免费看高清 | 深夜精品福利| 国产黄a三级三级三级人| 欧美日韩乱码在线| 18禁美女被吸乳视频| 欧美成人一区二区免费高清观看 | 国产日韩一区二区三区精品不卡| 人人澡人人妻人| 日日摸夜夜添夜夜添小说| 成人国产综合亚洲| 女人被躁到高潮嗷嗷叫费观| 日韩 欧美 亚洲 中文字幕| 久久久精品欧美日韩精品| bbb黄色大片| 男男h啪啪无遮挡| 自线自在国产av| 琪琪午夜伦伦电影理论片6080| 免费在线观看黄色视频的| 窝窝影院91人妻| 久久人人97超碰香蕉20202| 亚洲无线在线观看| e午夜精品久久久久久久| 亚洲专区字幕在线| 久久亚洲精品不卡| 日韩高清综合在线| 国产成人av激情在线播放| 亚洲国产欧美日韩在线播放| 国产一区二区三区在线臀色熟女| 精品久久蜜臀av无| 久久香蕉精品热| 1024香蕉在线观看| 无人区码免费观看不卡| 老司机深夜福利视频在线观看| 无遮挡黄片免费观看| 久久精品国产亚洲av高清一级| 老司机福利观看| 色老头精品视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲国产欧美一区二区综合| 我的亚洲天堂| 国产野战对白在线观看| 最近最新免费中文字幕在线| 可以在线观看毛片的网站| 午夜福利高清视频| 亚洲欧美一区二区三区黑人| 最近最新中文字幕大全电影3 | 亚洲精品在线观看二区| 狠狠狠狠99中文字幕| 日韩欧美在线二视频| 一级黄色大片毛片| 国产欧美日韩精品亚洲av| 黄色片一级片一级黄色片| 黑人操中国人逼视频| 亚洲美女黄片视频| 欧美成人性av电影在线观看| 9色porny在线观看| 色在线成人网| 免费人成视频x8x8入口观看| 91大片在线观看| 一区二区三区激情视频| 久久香蕉国产精品| 国产区一区二久久| 50天的宝宝边吃奶边哭怎么回事| 日本免费一区二区三区高清不卡 | 国产精品一区二区精品视频观看| 黄片小视频在线播放| 69精品国产乱码久久久| 中文字幕精品免费在线观看视频| 国产精品免费一区二区三区在线| 99在线视频只有这里精品首页| 亚洲av成人一区二区三| 少妇被粗大的猛进出69影院| 久久久精品国产亚洲av高清涩受| 少妇熟女aⅴ在线视频| 日韩欧美一区视频在线观看| 欧美日韩一级在线毛片| 757午夜福利合集在线观看| 久久婷婷人人爽人人干人人爱 | 人人妻,人人澡人人爽秒播| 99国产综合亚洲精品| 午夜福利成人在线免费观看| 国产一级毛片七仙女欲春2 | 十分钟在线观看高清视频www| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久精品电影 | 国内精品久久久久精免费| 午夜精品国产一区二区电影| 两个人视频免费观看高清| 亚洲精品国产色婷婷电影| 女同久久另类99精品国产91| 男女做爰动态图高潮gif福利片 | 欧美精品啪啪一区二区三区| 亚洲专区国产一区二区| 亚洲精品国产区一区二| 人人妻人人爽人人添夜夜欢视频| 午夜精品在线福利| 久久久久久久久免费视频了| 午夜久久久久精精品| 黑人操中国人逼视频| 亚洲电影在线观看av| 欧美亚洲日本最大视频资源| 国产单亲对白刺激| 伦理电影免费视频| 中文字幕高清在线视频| 免费在线观看影片大全网站| 国产不卡一卡二| 精品久久久久久久毛片微露脸| 国产精品 国内视频| 久久草成人影院| 亚洲av片天天在线观看| 色哟哟哟哟哟哟| 婷婷六月久久综合丁香| 性色av乱码一区二区三区2| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| 久久久久久免费高清国产稀缺| 久久婷婷人人爽人人干人人爱 | 久久国产精品人妻蜜桃| 亚洲专区国产一区二区| 精品人妻在线不人妻| 亚洲一码二码三码区别大吗| 91麻豆av在线| 757午夜福利合集在线观看| 巨乳人妻的诱惑在线观看| 亚洲欧美一区二区三区黑人| 熟女少妇亚洲综合色aaa.| 曰老女人黄片| 久久久国产精品麻豆| 中出人妻视频一区二区| 在线观看www视频免费| 国产精品自产拍在线观看55亚洲| 9色porny在线观看| 美女午夜性视频免费| 青草久久国产| 免费无遮挡裸体视频| 天堂√8在线中文| 国产精品亚洲av一区麻豆| 精品少妇一区二区三区视频日本电影| 丝袜美腿诱惑在线| 亚洲情色 制服丝袜| 亚洲美女黄片视频| 桃红色精品国产亚洲av| 天天一区二区日本电影三级 | 国产亚洲精品一区二区www| 亚洲激情在线av| 无人区码免费观看不卡| 成人精品一区二区免费| 美女免费视频网站| 国产欧美日韩一区二区精品| 变态另类成人亚洲欧美熟女 | 一级毛片精品| 日韩精品青青久久久久久| 久久亚洲真实| 精品日产1卡2卡| 亚洲av五月六月丁香网| 国产亚洲欧美在线一区二区| 免费久久久久久久精品成人欧美视频| 免费看a级黄色片| 十八禁人妻一区二区| 成人精品一区二区免费| 午夜福利一区二区在线看| 在线观看免费午夜福利视频| 中文亚洲av片在线观看爽| 国内毛片毛片毛片毛片毛片| 成人欧美大片| 婷婷六月久久综合丁香| 精品久久蜜臀av无| 精品人妻1区二区| 中国美女看黄片| 91在线观看av| av有码第一页| 黄色a级毛片大全视频| 99国产极品粉嫩在线观看| 亚洲av第一区精品v没综合| 国产精品,欧美在线| 精品一品国产午夜福利视频| 欧美成人性av电影在线观看| 欧美中文日本在线观看视频| 美国免费a级毛片| 日本欧美视频一区| 91国产中文字幕| 亚洲免费av在线视频| 九色国产91popny在线| 欧美精品亚洲一区二区| 男女做爰动态图高潮gif福利片 | 国产亚洲欧美在线一区二区| 成年人黄色毛片网站| 啦啦啦韩国在线观看视频| 自线自在国产av| 免费女性裸体啪啪无遮挡网站| www国产在线视频色| 9热在线视频观看99| 欧美黑人欧美精品刺激| 亚洲精品久久成人aⅴ小说| 国产国语露脸激情在线看| 一二三四在线观看免费中文在| 亚洲熟妇熟女久久| 日韩欧美国产在线观看| 日本欧美视频一区| 首页视频小说图片口味搜索| 久久欧美精品欧美久久欧美| 国产欧美日韩综合在线一区二区| 久久婷婷成人综合色麻豆| 免费人成视频x8x8入口观看| 99国产精品一区二区三区| 久久香蕉国产精品| 免费观看精品视频网站| 日韩欧美三级三区| 欧美大码av| 中文字幕人成人乱码亚洲影| 亚洲成国产人片在线观看| 日韩免费av在线播放| 一本久久中文字幕| 女同久久另类99精品国产91| a级毛片在线看网站| 亚洲一区二区三区不卡视频| 淫秽高清视频在线观看| 99国产精品99久久久久| 国产午夜精品久久久久久| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| 在线av久久热| 国产亚洲精品第一综合不卡| 欧美色欧美亚洲另类二区 | 操美女的视频在线观看| 日韩欧美国产一区二区入口| 1024视频免费在线观看| 欧美中文综合在线视频| 亚洲国产欧美一区二区综合| 精品久久久久久成人av| 精品久久久精品久久久| 欧美日本亚洲视频在线播放| 美女高潮到喷水免费观看| 夜夜夜夜夜久久久久| 一本大道久久a久久精品| 亚洲欧美日韩无卡精品| 亚洲aⅴ乱码一区二区在线播放 | 曰老女人黄片| 美女高潮到喷水免费观看| 黄片播放在线免费| 亚洲视频免费观看视频| 精品无人区乱码1区二区| www日本在线高清视频| 中文字幕精品免费在线观看视频| 高清在线国产一区| 国产亚洲精品一区二区www| 午夜影院日韩av| 亚洲国产精品成人综合色| 亚洲成国产人片在线观看| 黑人巨大精品欧美一区二区mp4| 黄片小视频在线播放| 女人被狂操c到高潮| 欧美不卡视频在线免费观看 | 午夜免费观看网址| 成人手机av| 男男h啪啪无遮挡| 国产高清视频在线播放一区| 国产精品av久久久久免费| 免费不卡黄色视频| 色婷婷久久久亚洲欧美| 人人妻人人爽人人添夜夜欢视频| 操出白浆在线播放| 制服人妻中文乱码| 最好的美女福利视频网| 久久久久久久午夜电影| 久久精品91无色码中文字幕| 亚洲av电影不卡..在线观看| 黑人操中国人逼视频| 国产乱人伦免费视频| 免费高清视频大片| 91国产中文字幕| 级片在线观看| 成人国产一区最新在线观看| 51午夜福利影视在线观看| 精品第一国产精品| 999久久久国产精品视频| 十八禁网站免费在线| 国产熟女午夜一区二区三区| 中文字幕久久专区| 亚洲国产精品sss在线观看| 日韩av在线大香蕉| 制服人妻中文乱码| 亚洲成人免费电影在线观看| 99精品欧美一区二区三区四区| 日本撒尿小便嘘嘘汇集6| 久久精品亚洲精品国产色婷小说| 午夜影院日韩av| 中国美女看黄片| 少妇的丰满在线观看| 岛国在线观看网站| 一二三四社区在线视频社区8| 亚洲五月婷婷丁香| 欧美大码av| 久久欧美精品欧美久久欧美| 亚洲国产看品久久| 天堂动漫精品| 欧美黄色淫秽网站| 无限看片的www在线观看| 女人高潮潮喷娇喘18禁视频| а√天堂www在线а√下载| 首页视频小说图片口味搜索| 搞女人的毛片| 午夜福利18| 亚洲色图 男人天堂 中文字幕| 91老司机精品| av免费在线观看网站| videosex国产| 国产区一区二久久| 精品一区二区三区视频在线观看免费| 在线免费观看的www视频| 香蕉国产在线看| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 午夜精品国产一区二区电影| 国产三级黄色录像| 人成视频在线观看免费观看| 脱女人内裤的视频| 亚洲自拍偷在线| 免费看十八禁软件| 亚洲情色 制服丝袜| 午夜免费成人在线视频| 午夜亚洲福利在线播放| 怎么达到女性高潮| 亚洲一区二区三区色噜噜| 久久久水蜜桃国产精品网| 色综合欧美亚洲国产小说| 亚洲国产欧美日韩在线播放| 亚洲熟妇熟女久久| www.熟女人妻精品国产| 18禁观看日本| 深夜精品福利| 欧美中文日本在线观看视频| 国产99久久九九免费精品| 涩涩av久久男人的天堂| 欧美日韩亚洲综合一区二区三区_| 久久人人爽av亚洲精品天堂| 欧美性长视频在线观看| 十八禁网站免费在线| 国产精品美女特级片免费视频播放器 | 欧美精品啪啪一区二区三区| 国产精品影院久久| 黄频高清免费视频| 日本三级黄在线观看| 亚洲美女黄片视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲午夜理论影院| 久久久久久人人人人人| 免费女性裸体啪啪无遮挡网站| 中文字幕人妻熟女乱码| 亚洲av电影不卡..在线观看| 国产欧美日韩一区二区三区在线| 午夜福利18| 免费高清视频大片| bbb黄色大片| 69精品国产乱码久久久| 99久久综合精品五月天人人| 天天一区二区日本电影三级 | 色精品久久人妻99蜜桃| 亚洲精品久久成人aⅴ小说| 真人一进一出gif抽搐免费| 制服人妻中文乱码| 99久久久亚洲精品蜜臀av| 电影成人av| 亚洲精品av麻豆狂野| 18禁国产床啪视频网站| 极品人妻少妇av视频| 午夜福利高清视频| 久久久久久久午夜电影| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩一区二区三区在线| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 啦啦啦免费观看视频1| a级毛片在线看网站| 国产精品乱码一区二三区的特点 | 可以免费在线观看a视频的电影网站| 一区二区三区激情视频| 一级片免费观看大全| 国产激情欧美一区二区| 黑人巨大精品欧美一区二区蜜桃| 国产精华一区二区三区| 国产1区2区3区精品| tocl精华| 日本a在线网址| 久久久国产成人精品二区| 黄色a级毛片大全视频| 久久人人97超碰香蕉20202| 日本精品一区二区三区蜜桃| 国产熟女午夜一区二区三区| √禁漫天堂资源中文www| 麻豆成人av在线观看| 国产野战对白在线观看| 国产精品爽爽va在线观看网站 | www国产在线视频色| 91九色精品人成在线观看| 亚洲五月色婷婷综合| 久久国产乱子伦精品免费另类| 亚洲欧美一区二区三区黑人| 极品人妻少妇av视频| 亚洲最大成人中文| 精品人妻1区二区| 十八禁网站免费在线| 国产国语露脸激情在线看| 黑人欧美特级aaaaaa片| 欧美亚洲日本最大视频资源| 久久婷婷人人爽人人干人人爱 | av片东京热男人的天堂| 成人av一区二区三区在线看| 又大又爽又粗| 三级毛片av免费| 岛国在线观看网站| 日韩视频一区二区在线观看| 久久午夜亚洲精品久久| 在线国产一区二区在线| АⅤ资源中文在线天堂| 人人妻人人澡人人看| 欧美最黄视频在线播放免费| 午夜福利18| 极品教师在线免费播放| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| 午夜福利免费观看在线| 可以在线观看的亚洲视频| 国产精品一区二区精品视频观看| 变态另类丝袜制服| 精品熟女少妇八av免费久了| 看黄色毛片网站| 女生性感内裤真人,穿戴方法视频| 男女下面插进去视频免费观看| 欧美乱码精品一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲成av人片免费观看| av免费在线观看网站| 琪琪午夜伦伦电影理论片6080| 亚洲情色 制服丝袜| 亚洲自拍偷在线| 久久九九热精品免费| 又黄又爽又免费观看的视频| 国产精品二区激情视频| 女警被强在线播放| 好男人电影高清在线观看| av网站免费在线观看视频| 久久精品亚洲熟妇少妇任你| 亚洲男人天堂网一区| 国产精品秋霞免费鲁丝片| 人人妻,人人澡人人爽秒播| 免费搜索国产男女视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人免费无遮挡视频| 日本免费一区二区三区高清不卡 | 99国产精品99久久久久| 天天躁夜夜躁狠狠躁躁| 亚洲成av人片免费观看| 老司机深夜福利视频在线观看| 少妇被粗大的猛进出69影院| 熟女少妇亚洲综合色aaa.| 成人av一区二区三区在线看| 女性生殖器流出的白浆| 国产97色在线日韩免费| 亚洲欧美激情在线| 一区二区三区精品91| 欧美精品啪啪一区二区三区| 搡老妇女老女人老熟妇| 黑人欧美特级aaaaaa片| 亚洲av成人不卡在线观看播放网| av福利片在线| 欧美成人午夜精品| 亚洲七黄色美女视频| 久久国产精品人妻蜜桃| 人人妻,人人澡人人爽秒播| aaaaa片日本免费| 国产视频一区二区在线看| 久久九九热精品免费| 亚洲色图av天堂| 99国产精品一区二区蜜桃av| 脱女人内裤的视频| АⅤ资源中文在线天堂| √禁漫天堂资源中文www| 少妇粗大呻吟视频| 精品熟女少妇八av免费久了| 一级黄色大片毛片| 国产精品久久视频播放| 亚洲专区字幕在线| 亚洲av第一区精品v没综合| 国产日韩一区二区三区精品不卡| 自线自在国产av| 久久这里只有精品19| 日韩中文字幕欧美一区二区| 一级,二级,三级黄色视频| 免费观看精品视频网站| 亚洲少妇的诱惑av| www国产在线视频色| 午夜福利高清视频| 99久久精品国产亚洲精品| 9热在线视频观看99| 18美女黄网站色大片免费观看| 两个人看的免费小视频| 亚洲国产高清在线一区二区三 | 免费一级毛片在线播放高清视频 | 国产精品爽爽va在线观看网站 | 黄片小视频在线播放| 免费在线观看黄色视频的| 老汉色∧v一级毛片| 18禁美女被吸乳视频| www.精华液| 乱人伦中国视频| 亚洲黑人精品在线| 亚洲av成人一区二区三| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区免费欧美| 久久香蕉国产精品| 91老司机精品| 在线天堂中文资源库| 国产av一区在线观看免费| 嫩草影视91久久| 性少妇av在线| 精品电影一区二区在线| 国产97色在线日韩免费| 大型黄色视频在线免费观看| 午夜老司机福利片| 亚洲成人免费电影在线观看| 真人一进一出gif抽搐免费| 久久香蕉精品热| 国产熟女午夜一区二区三区| 日韩精品中文字幕看吧| 女性生殖器流出的白浆| 国产亚洲精品一区二区www| 黄色视频不卡| 久久精品国产综合久久久| 免费高清视频大片| 色综合婷婷激情| 国产午夜精品久久久久久| 在线观看舔阴道视频| 变态另类丝袜制服| 国产亚洲精品av在线| 久久香蕉国产精品| 免费观看精品视频网站| 国产蜜桃级精品一区二区三区| 久久久精品国产亚洲av高清涩受| 夜夜看夜夜爽夜夜摸| 久久久国产成人免费| 午夜两性在线视频| 香蕉丝袜av| 亚洲国产精品合色在线| 18禁美女被吸乳视频| 国产亚洲av高清不卡| 母亲3免费完整高清在线观看|