• <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| 午夜视频精品福利| 禁无遮挡网站| 亚洲中文字幕一区二区三区有码在线看 | 桃红色精品国产亚洲av| 中文字幕人成人乱码亚洲影| 黄色成人免费大全| 国产成人av激情在线播放| 最新在线观看一区二区三区| 国产精品av久久久久免费| 国产伦精品一区二区三区视频9 | 国产日本99.免费观看| 桃红色精品国产亚洲av| 国产精品九九99| 日本 av在线| 老司机深夜福利视频在线观看| 精品国内亚洲2022精品成人| 日韩欧美国产在线观看| 丰满人妻熟妇乱又伦精品不卡| 这个男人来自地球电影免费观看| 日日干狠狠操夜夜爽| 午夜福利在线在线| 91九色精品人成在线观看| 国内久久婷婷六月综合欲色啪| 亚洲国产色片| 在线观看美女被高潮喷水网站 | 99热这里只有精品一区 | 色av中文字幕| 一本久久中文字幕| 我要搜黄色片| 免费电影在线观看免费观看| 久久久久久人人人人人| 亚洲最大成人中文| 日本一本二区三区精品| 欧美色欧美亚洲另类二区| 午夜福利在线观看吧| 亚洲国产欧美人成| 男女下面进入的视频免费午夜| 免费在线观看影片大全网站| 免费看光身美女| 精品久久久久久久久久免费视频| 夜夜夜夜夜久久久久| 免费av毛片视频| 97碰自拍视频| 亚洲专区国产一区二区| 麻豆一二三区av精品| 国产av一区在线观看免费| 亚洲一区高清亚洲精品| 非洲黑人性xxxx精品又粗又长| 国产成人影院久久av| 久久精品91无色码中文字幕| 日韩av在线大香蕉| 国产成人福利小说| 国产蜜桃级精品一区二区三区| 亚洲精品一区av在线观看| 国产精品综合久久久久久久免费| 无限看片的www在线观看| 淫妇啪啪啪对白视频| 久久草成人影院| 青草久久国产| 欧美日韩乱码在线| 久久久水蜜桃国产精品网| 国产一区二区在线观看日韩 | 国产成人一区二区三区免费视频网站| a级毛片在线看网站| 免费搜索国产男女视频| 免费大片18禁| 国产欧美日韩精品亚洲av| 久久99热这里只有精品18| 嫩草影视91久久| av福利片在线观看| 国产熟女xx| 亚洲欧美日韩高清在线视频| 99热这里只有精品一区 | 在线播放国产精品三级| 国产日本99.免费观看| 久久性视频一级片| 最近最新中文字幕大全免费视频| 天堂√8在线中文| 国产主播在线观看一区二区| 嫩草影视91久久| 动漫黄色视频在线观看| 国内精品一区二区在线观看| 人人妻人人看人人澡| 九色国产91popny在线| 亚洲精品456在线播放app | 成年人黄色毛片网站| 免费看十八禁软件| 午夜亚洲福利在线播放| 91av网站免费观看| 国产精品日韩av在线免费观看| 在线永久观看黄色视频| 日韩成人在线观看一区二区三区| 国产精品1区2区在线观看.| 久9热在线精品视频| 精品久久久久久久末码| 黄色视频,在线免费观看| 亚洲国产欧洲综合997久久,| 男女做爰动态图高潮gif福利片| 亚洲乱码一区二区免费版| 国内精品美女久久久久久| 巨乳人妻的诱惑在线观看| 久久精品国产清高在天天线| 午夜亚洲福利在线播放| 99久久国产精品久久久| 精品欧美国产一区二区三| 日韩欧美在线乱码| 久久久国产欧美日韩av| 999久久久国产精品视频| 999久久久精品免费观看国产| 亚洲精品一卡2卡三卡4卡5卡| 日韩精品中文字幕看吧| 熟妇人妻久久中文字幕3abv| 在线国产一区二区在线| 精品久久久久久久人妻蜜臀av| 波多野结衣巨乳人妻| 午夜日韩欧美国产| 精品久久久久久久久久久久久| 欧美黄色淫秽网站| 色视频www国产| 真人一进一出gif抽搐免费| 一边摸一边抽搐一进一小说| 日本a在线网址| 免费大片18禁| 久久中文字幕一级| 舔av片在线| 色综合站精品国产| 九色成人免费人妻av| www.999成人在线观看| 国产av不卡久久| 中文字幕人妻丝袜一区二区| 综合色av麻豆| 性色av乱码一区二区三区2| 香蕉国产在线看| 色视频www国产| 日本在线视频免费播放| 啦啦啦韩国在线观看视频| 真人做人爱边吃奶动态| 9191精品国产免费久久| 国产精品电影一区二区三区| 美女被艹到高潮喷水动态| 在线观看舔阴道视频| 曰老女人黄片| 国产一级毛片七仙女欲春2| av福利片在线观看| 久久香蕉精品热| 黄色成人免费大全| 在线观看免费视频日本深夜| 一进一出抽搐gif免费好疼| 国产三级在线视频| 18禁黄网站禁片免费观看直播| 男女做爰动态图高潮gif福利片| 99久久久亚洲精品蜜臀av| 999精品在线视频| 我的老师免费观看完整版| 我的老师免费观看完整版| 99久久久亚洲精品蜜臀av| 一进一出抽搐gif免费好疼| 国产精品综合久久久久久久免费| 国产精品综合久久久久久久免费| 一进一出抽搐gif免费好疼| 日韩免费av在线播放| 久久久久久九九精品二区国产| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产高清在线一区二区三| 成人特级黄色片久久久久久久| 日本黄色视频三级网站网址| 国产99白浆流出| 日韩欧美 国产精品| 人人妻,人人澡人人爽秒播| 最好的美女福利视频网| 日本 av在线| 国产精品亚洲一级av第二区| 夜夜看夜夜爽夜夜摸| 久久欧美精品欧美久久欧美| 伊人久久大香线蕉亚洲五| 欧美又色又爽又黄视频| 2021天堂中文幕一二区在线观| 精品人妻1区二区| 亚洲欧美日韩无卡精品| 五月伊人婷婷丁香| 亚洲片人在线观看| 国产精品久久久久久精品电影| 亚洲激情在线av| 亚洲av成人av| 亚洲av电影在线进入| 亚洲国产精品sss在线观看| 亚洲电影在线观看av| 12—13女人毛片做爰片一| 97碰自拍视频| 一级a爱片免费观看的视频| 久久这里只有精品中国| 丰满的人妻完整版| 亚洲精品美女久久久久99蜜臀| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| 九色国产91popny在线| 国产成人av激情在线播放| 国产伦精品一区二区三区四那| 757午夜福利合集在线观看| 亚洲成人免费电影在线观看| 欧美又色又爽又黄视频| 一进一出好大好爽视频| 极品教师在线免费播放| 亚洲色图 男人天堂 中文字幕| 99久久精品国产亚洲精品| www国产在线视频色| 国产精品 欧美亚洲| 亚洲人成网站高清观看| 日本熟妇午夜| 一区二区三区激情视频| 悠悠久久av| 999精品在线视频| 国产男靠女视频免费网站| 久久人人精品亚洲av| 久久婷婷人人爽人人干人人爱| 十八禁网站免费在线| 亚洲无线在线观看| 性色av乱码一区二区三区2| 中国美女看黄片| 在线观看免费午夜福利视频| 欧美3d第一页| 伦理电影免费视频| 一区二区三区国产精品乱码| 亚洲av五月六月丁香网| 国内少妇人妻偷人精品xxx网站 | 啪啪无遮挡十八禁网站| 99久久精品国产亚洲精品| www国产在线视频色| 欧美另类亚洲清纯唯美| 国产精品野战在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品98久久久久久宅男小说| 麻豆久久精品国产亚洲av| 国产三级在线视频| 一进一出好大好爽视频| 九九热线精品视视频播放| 91久久精品国产一区二区成人 | 黄色女人牲交| 欧美日韩中文字幕国产精品一区二区三区| 日本a在线网址| 欧美色视频一区免费| 国产视频内射| 级片在线观看| 亚洲成av人片在线播放无| 国产精品98久久久久久宅男小说| 国产高潮美女av| 国产又色又爽无遮挡免费看| 亚洲精品美女久久久久99蜜臀| 欧美性猛交╳xxx乱大交人| 一边摸一边抽搐一进一小说| 美女 人体艺术 gogo| 观看免费一级毛片| 亚洲av五月六月丁香网| av天堂中文字幕网| 日本免费一区二区三区高清不卡| 精华霜和精华液先用哪个| 国产高清视频在线观看网站| cao死你这个sao货| 亚洲色图 男人天堂 中文字幕| 成人鲁丝片一二三区免费| 国产在线精品亚洲第一网站| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| 国产三级在线视频| 久久久久久久久久黄片| 九九在线视频观看精品| 中文字幕av在线有码专区| 亚洲人成网站在线播放欧美日韩| 婷婷精品国产亚洲av在线| 99精品欧美一区二区三区四区| 最近最新免费中文字幕在线| 欧美成狂野欧美在线观看| 国产单亲对白刺激| 99久久99久久久精品蜜桃| 老司机在亚洲福利影院| 日韩人妻高清精品专区| 男女做爰动态图高潮gif福利片| 成人特级av手机在线观看| 免费搜索国产男女视频| 女警被强在线播放| 久久伊人香网站| 99久久99久久久精品蜜桃| ponron亚洲| 岛国视频午夜一区免费看| 日韩高清综合在线| 成年免费大片在线观看| 亚洲欧美日韩高清在线视频| 曰老女人黄片| 国产精品久久久av美女十八| 舔av片在线| 真人一进一出gif抽搐免费| 亚洲天堂国产精品一区在线| 国产精品1区2区在线观看.| 成人三级黄色视频| 国产精品爽爽va在线观看网站| 国产精品99久久99久久久不卡| 欧美日本亚洲视频在线播放| 亚洲国产欧洲综合997久久,| 成人无遮挡网站| 法律面前人人平等表现在哪些方面| 国产毛片a区久久久久| 午夜免费成人在线视频| 亚洲国产欧美人成| 一本一本综合久久| 国产激情欧美一区二区| 97碰自拍视频| 观看美女的网站| 很黄的视频免费| 国产av麻豆久久久久久久| 国产极品精品免费视频能看的| 久久久成人免费电影| 亚洲国产精品成人综合色| 成人亚洲精品av一区二区| 香蕉久久夜色| 亚洲成人久久性| 最新中文字幕久久久久 | 中亚洲国语对白在线视频| 成年免费大片在线观看| 免费一级毛片在线播放高清视频| 成人特级黄色片久久久久久久| 男人舔女人下体高潮全视频| 亚洲欧美一区二区三区黑人| 成年免费大片在线观看| 国产av麻豆久久久久久久| 国产一级毛片七仙女欲春2| 久久天堂一区二区三区四区| 丝袜人妻中文字幕| 色尼玛亚洲综合影院| 亚洲自偷自拍图片 自拍| 精品电影一区二区在线| 九九热线精品视视频播放| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 一级毛片精品| 两个人的视频大全免费| 国产91精品成人一区二区三区| av视频在线观看入口| 久久精品国产99精品国产亚洲性色| 99久久99久久久精品蜜桃| 999久久久国产精品视频| 亚洲国产欧洲综合997久久,| a在线观看视频网站| 91麻豆av在线| 在线观看舔阴道视频| 男女午夜视频在线观看| 我要搜黄色片| 色视频www国产| 国产在线精品亚洲第一网站| 黑人巨大精品欧美一区二区mp4| 少妇丰满av| 国产伦人伦偷精品视频| 成年女人毛片免费观看观看9| 波多野结衣巨乳人妻| 我的老师免费观看完整版| 五月玫瑰六月丁香| 亚洲成av人片免费观看| 九色国产91popny在线| 在线观看免费视频日本深夜| 男女做爰动态图高潮gif福利片| 99国产精品一区二区蜜桃av| 久久久久国产一级毛片高清牌| 国产成+人综合+亚洲专区| 精品人妻1区二区| 国产一区二区三区在线臀色熟女| 黄色片一级片一级黄色片| 国产精品野战在线观看| 午夜福利高清视频| 亚洲国产欧美人成| 黄色片一级片一级黄色片| 最近在线观看免费完整版| 天堂影院成人在线观看| 日本成人三级电影网站| 欧美午夜高清在线| 久久久久九九精品影院| 欧美一级a爱片免费观看看| 亚洲电影在线观看av| 欧美一区二区国产精品久久精品| 日本撒尿小便嘘嘘汇集6| 国产精品一区二区三区四区免费观看 | 国内精品一区二区在线观看| 亚洲乱码一区二区免费版| 我要搜黄色片| 亚洲第一欧美日韩一区二区三区| 亚洲一区二区三区不卡视频| 亚洲美女黄片视频| 亚洲精品中文字幕一二三四区| 成人亚洲精品av一区二区| 网址你懂的国产日韩在线| 色综合亚洲欧美另类图片| 一本久久中文字幕| 午夜精品一区二区三区免费看| 一区二区三区高清视频在线| 色综合站精品国产| 97超级碰碰碰精品色视频在线观看| 国产伦在线观看视频一区| 成人三级黄色视频| 日韩av在线大香蕉| 伦理电影免费视频| 亚洲国产高清在线一区二区三| 国产高清视频在线观看网站| 99精品久久久久人妻精品| 日韩欧美在线乱码| 国产又色又爽无遮挡免费看| 久久精品91无色码中文字幕| 国产精品一区二区免费欧美| 亚洲精品一区av在线观看| 久久久久久久精品吃奶| 两个人看的免费小视频| 男女午夜视频在线观看| 最近最新中文字幕大全免费视频| 欧美国产日韩亚洲一区| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美人成| 国产97色在线日韩免费| 国产精品久久久久久久电影 | 女人高潮潮喷娇喘18禁视频| 精品免费久久久久久久清纯| 国产精品久久电影中文字幕| 黄色丝袜av网址大全| 久久精品亚洲精品国产色婷小说| avwww免费| 两性夫妻黄色片| 99久久精品一区二区三区| 午夜福利在线在线| 国产成人精品久久二区二区91| 国产综合懂色| ponron亚洲| 成人午夜高清在线视频| 国产高清三级在线| av中文乱码字幕在线| 亚洲五月婷婷丁香| 亚洲av成人一区二区三| 人人妻,人人澡人人爽秒播| 国产精品亚洲一级av第二区| 欧美在线黄色| 亚洲精品久久国产高清桃花| 一级毛片女人18水好多| 2021天堂中文幕一二区在线观| 免费人成视频x8x8入口观看| 国产精品一区二区免费欧美| 欧美+亚洲+日韩+国产| 国产高清有码在线观看视频| 亚洲欧洲精品一区二区精品久久久| 免费在线观看影片大全网站| 我要搜黄色片| 91老司机精品| 日本与韩国留学比较| 成人无遮挡网站| 色综合婷婷激情| 美女午夜性视频免费| 亚洲中文av在线| 精品日产1卡2卡| 国产99白浆流出| 91在线精品国自产拍蜜月 | 女人被狂操c到高潮| 熟女电影av网| 精品国产乱码久久久久久男人| 99久久成人亚洲精品观看| 国产一级毛片七仙女欲春2| 免费大片18禁| 亚洲人成网站高清观看| 欧美日本亚洲视频在线播放| 免费人成视频x8x8入口观看| 日日夜夜操网爽| 丁香欧美五月| 一进一出抽搐动态| 在线观看免费视频日本深夜| 欧美性猛交╳xxx乱大交人| 亚洲天堂国产精品一区在线| 日韩免费av在线播放| 久久久久久久精品吃奶| 亚洲欧美日韩高清在线视频| 成人三级黄色视频| 男人和女人高潮做爰伦理| 国产免费男女视频| 精品乱码久久久久久99久播| 国产69精品久久久久777片 | 搡老妇女老女人老熟妇| 精品一区二区三区av网在线观看| 色在线成人网| 午夜精品在线福利| 无限看片的www在线观看| 久久久久久久久久黄片| 又爽又黄无遮挡网站| 不卡一级毛片| 亚洲精品中文字幕一二三四区| 中文字幕久久专区| 久久亚洲精品不卡| 亚洲成a人片在线一区二区| 久久精品91蜜桃| 少妇裸体淫交视频免费看高清| 蜜桃久久精品国产亚洲av| 精品电影一区二区在线| 成人无遮挡网站| 在线十欧美十亚洲十日本专区| 国产黄a三级三级三级人| 精品久久久久久成人av| 亚洲真实伦在线观看| 免费观看的影片在线观看| 天天一区二区日本电影三级| 午夜福利在线在线| 午夜免费激情av| 在线十欧美十亚洲十日本专区| 亚洲精品粉嫩美女一区| 日韩欧美国产在线观看| 国产一区在线观看成人免费| 精品熟女少妇八av免费久了| 国产精品98久久久久久宅男小说| 香蕉丝袜av| 中文字幕av在线有码专区| 亚洲精品美女久久久久99蜜臀| 亚洲欧洲精品一区二区精品久久久| 啦啦啦免费观看视频1| 在线免费观看不下载黄p国产 | 黑人操中国人逼视频| 男女床上黄色一级片免费看| 国内久久婷婷六月综合欲色啪| 男女午夜视频在线观看| tocl精华| av福利片在线观看| 国内毛片毛片毛片毛片毛片| 亚洲 欧美一区二区三区| 在线永久观看黄色视频| 亚洲精品456在线播放app | 欧美日韩综合久久久久久 | 国产野战对白在线观看| 美女cb高潮喷水在线观看 | 久久久精品欧美日韩精品| 中文字幕久久专区| 少妇丰满av| 精品一区二区三区视频在线 | 成人高潮视频无遮挡免费网站| 99在线人妻在线中文字幕| 精品久久久久久久久久久久久| 欧美在线黄色| 日本黄色片子视频| 日韩人妻高清精品专区| 两性夫妻黄色片| 在线观看66精品国产| 两个人视频免费观看高清| 午夜日韩欧美国产| 长腿黑丝高跟| av黄色大香蕉| 久久精品aⅴ一区二区三区四区| 国产成人欧美在线观看| 成人国产一区最新在线观看| 欧美色视频一区免费| 欧美三级亚洲精品| 这个男人来自地球电影免费观看| 手机成人av网站| 最近最新中文字幕大全免费视频| 国产亚洲精品av在线| 一级黄色大片毛片| 日日干狠狠操夜夜爽| 熟妇人妻久久中文字幕3abv| 亚洲性夜色夜夜综合| 亚洲午夜理论影院| 国产午夜精品久久久久久| 看免费av毛片| 变态另类丝袜制服| 国产在线精品亚洲第一网站| 精品国产美女av久久久久小说| 国产精品 欧美亚洲| 九九久久精品国产亚洲av麻豆 | 少妇的逼水好多| 一二三四社区在线视频社区8| 观看免费一级毛片| 精品国产美女av久久久久小说| 国产激情久久老熟女| 99re在线观看精品视频| 成年女人毛片免费观看观看9| 2021天堂中文幕一二区在线观| 午夜激情福利司机影院| 久久精品人妻少妇| 操出白浆在线播放| 午夜视频精品福利| 别揉我奶头~嗯~啊~动态视频| 国产一区二区三区视频了| 国产亚洲av嫩草精品影院| 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说| 黄色日韩在线| 18美女黄网站色大片免费观看| 一个人看的www免费观看视频| 日韩av在线大香蕉| 日韩精品中文字幕看吧| 成人午夜高清在线视频| 岛国在线观看网站| 亚洲欧美日韩卡通动漫| 岛国在线观看网站| 91久久精品国产一区二区成人 | 99久国产av精品| av国产免费在线观看| 亚洲男人的天堂狠狠| 久久亚洲精品不卡| 色尼玛亚洲综合影院| 在线永久观看黄色视频| 真实男女啪啪啪动态图| 亚洲av免费在线观看| 怎么达到女性高潮| 床上黄色一级片| 天堂av国产一区二区熟女人妻| or卡值多少钱| 欧美日韩国产亚洲二区| 亚洲九九香蕉| 成人高潮视频无遮挡免费网站| 91在线精品国自产拍蜜月 | 午夜日韩欧美国产| 欧美性猛交╳xxx乱大交人| 国产亚洲精品一区二区www| 香蕉丝袜av| 岛国在线观看网站| a级毛片在线看网站|