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

    Al-7Si-0.36Mg合金定向凝固一次枝晶臂間距實(shí)驗(yàn)和模擬

    2015-11-19 09:40:34許慶彥柳百成清華大學(xué)材料學(xué)院先進(jìn)成形制造教育部重點(diǎn)實(shí)驗(yàn)室北京100084
    中國有色金屬學(xué)報 2015年10期
    關(guān)鍵詞:枝晶限值定向

    陳 瑞,許慶彥,柳百成(清華大學(xué) 材料學(xué)院 先進(jìn)成形制造教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)

    Al-7Si-0.36Mg合金定向凝固一次枝晶臂間距實(shí)驗(yàn)和模擬

    陳瑞,許慶彥,柳百成
    (清華大學(xué) 材料學(xué)院 先進(jìn)成形制造教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)

    通過Al-7Si-0.36Mg合金定向凝固實(shí)驗(yàn)和元胞自動機(jī)模型,開展定向凝固枝晶形貌演化和一次枝晶臂間距選擇過程的實(shí)驗(yàn)和模擬。結(jié)果表明:在給定的凝固條件下,一次枝晶臂間距范圍是一個連續(xù)的變化區(qū)間。在恒定溫度梯度和不同凝固速度條件下,測得Al-7Si-0.36Mg合金一次枝晶臂間距上限值(λmax)、下限值(λmin)和平均值(λave)以及生長速率之間的關(guān)系,且上限值和下限值的比值接近3。模擬結(jié)果與實(shí)驗(yàn)結(jié)果的吻合程度明顯優(yōu)于Hunt-Lu等解析模型的預(yù)測結(jié)果,表明CA模型在枝晶定向凝固過程枝晶形貌演化模擬和枝晶臂間距預(yù)測等方面的準(zhǔn)確性。結(jié)合模擬研究和文獻(xiàn)調(diào)研分析影響定向凝固一次枝晶臂選擇的因素,包括抽拉速度v、溫度梯度G、界面能大小、溶質(zhì)擴(kuò)散系數(shù)DL、枝晶生長取向與熱流方向的偏離角度θ等。

    定向凝固;枝晶臂間距;生長速度;元胞自動機(jī);數(shù)值模擬

    定向凝固鑄件中枝晶形貌及特征直接影響著鑄件最終力學(xué)性能,因此,通過深入了解和掌握定向凝固過程中枝晶的演化規(guī)律,有助于控制枝晶的形成過程,進(jìn)而達(dá)到對鑄件性能進(jìn)行有效預(yù)測和控制的目的[1]。枝晶臂間距(尤其是一次枝晶臂間距),是用于表征定向凝固枝晶細(xì)化程度的重要參數(shù),枝晶臂間距的大小直接影響著溶質(zhì)成分分布、中間相、共晶相、雜質(zhì)相的尺寸和分布以及裂紋和縮松等缺陷的形成[2-4],因此,研究定向凝固不同工藝條件下枝晶臂間距的演化過程并進(jìn)行準(zhǔn)確的定量預(yù)測,對控制凝固組織、提高鑄件性能具有十分重要的意義。

    過去的幾十年中,許多學(xué)者基于一些假設(shè)條件建立了枝晶臂間距與凝固參數(shù)及合金特性相關(guān)聯(lián)的解析模型[5-7]。HUNT[5]、KURZ等[6]、TRIVEDI[7]分別通過假設(shè)不同的枝晶尖端形狀(球形或是橢圓形)和完全一致的枝晶列排布,并利用臨界穩(wěn)定性判據(jù)(Marginal stability criterion)建立了穩(wěn)態(tài)條件下定向凝固枝晶一次枝晶臂間距λ1與固液界面溫度梯度G和生長速度v的單值函數(shù)關(guān)系。這些模型可以用式(1)統(tǒng)一表示:

    式中:A表示與材料特性參數(shù)相關(guān)的函數(shù);C0是合金成分。從式(1)可以看到,在材料和成分固定的情況下,即A和C0是定值,一次枝晶臂間距λ1只由固液界面的溫度梯度G和生長速度v決定,并且在同一種條件下,λ1具有唯一值。隨著枝晶生長界面穩(wěn)定性、枝晶側(cè)向分枝機(jī)制等因素對枝晶臂間距影響研究的深化[8-10]以及一些實(shí)驗(yàn)數(shù)據(jù)的驗(yàn)證[11],研究人員逐漸意識到枝晶的演化不僅受實(shí)驗(yàn)條件參數(shù)的影響,而且與演化歷史密切相關(guān),所以在同一凝固參數(shù)下,一次枝晶臂間距λ1并不是一個固定值,而是一個連續(xù)變化區(qū)間。HUNT等[12]通過分析凝固界面前沿溶質(zhì)橫向擴(kuò)散所引起的交互作用影響以及基于枝晶陣列穩(wěn)定性分析,推導(dǎo)出了一次枝晶臂間距與溫度梯度、凝固速度等因素的關(guān)系,確定了穩(wěn)態(tài)和非穩(wěn)態(tài)條件下一次枝晶臂間距的下限值[12]:

    為了能更準(zhǔn)確描述枝晶臂間距的演化過程,模型中必須要考慮枝晶臂的淹沒和側(cè)向分枝的作用,而這兩種機(jī)制主要受枝晶凝固界面穩(wěn)定性(界面張力和界面各向異性)、溶質(zhì)交互作用以及枝晶復(fù)雜的空間擴(kuò)展結(jié)構(gòu)等因素的影響。顯微組織數(shù)值模擬技術(shù)作為一種能準(zhǔn)確描述這些因素變化的方法,已經(jīng)成功應(yīng)用在定向凝固過程枝晶形貌和枝晶臂演化等方面的預(yù)測[13-18],并且證明了數(shù)值模擬預(yù)測的結(jié)果要比上述解析模型預(yù)測的結(jié)果更為準(zhǔn)確[13, 15]。目前,顯微組織數(shù)值模擬的主要方法有元胞自動機(jī)(CA)法和相場(PF)法,其中,CA方法具有較高的計算效率,通過耦合熱擴(kuò)散、溶質(zhì)擴(kuò)散以及枝晶生長動力學(xué)等,建立較復(fù)雜的描述枝晶固/液界面形態(tài)以及溶質(zhì)分布的數(shù)學(xué)方程,從而在更大尺度上實(shí)現(xiàn)對枝晶形貌演化過程的模擬,近些年來發(fā)展較快,在實(shí)際中得到了廣泛的應(yīng)用。WANG等[13]建立了溶質(zhì)擴(kuò)散控制的CA-FD模型,該模型考慮了成分過冷和曲率過冷,并改進(jìn)了GANDIN 等[19]提出的用于模擬晶粒生長的偏心算法(Decentered Square Algorithm),實(shí)現(xiàn)了任意取向枝晶生長的模擬,他們采用該模型詳細(xì)研究二元鎳基高溫合金定向凝固過程中初始枝晶臂數(shù)目對穩(wěn)態(tài)枝晶整列一次枝晶臂間距的影響,同時也分析了凝固速度對枝晶生長的影響,結(jié)果表明初始枝晶臂數(shù)目對穩(wěn)態(tài)時的枝晶臂間距影響不大,而凝固速度對枝晶臂間距的影響很大,同時也給出了枝晶臂間距上限值與下限值呈接近3倍的關(guān)系。隨后,單博煒等[15]也建立了相應(yīng)的CA模型模擬分析了定向凝固一次枝晶臂間距的選擇過程,分析了初始晶核數(shù)目和變抽拉速度對于一次枝晶臂間距允許范圍的影響,結(jié)果表明一次枝晶臂間距與這兩種數(shù)值實(shí)驗(yàn)方法無關(guān),同時給出了一次枝晶臂間距的上線值和下限值相對于凝固速度的冪函數(shù)關(guān)系式,并與實(shí)驗(yàn)結(jié)果進(jìn)行了良好對比。ZHANG等[16]采用3D-CA模型模擬了Al-11.6Cu-0.85Mg(質(zhì)量分?jǐn)?shù),%)合金定向凝固枝晶演變過程,指出三次枝晶臂能否發(fā)展成為穩(wěn)態(tài)一次枝晶臂主要由固液界面穩(wěn)定性以及相鄰一次枝晶臂間溶質(zhì)交互作用程度共同決定,同時,模擬得到的平均一次枝晶臂間距與實(shí)驗(yàn)結(jié)果相吻合。目前,大多數(shù)的研究局限于分析抽拉速度和溫度梯度等工藝參數(shù)的影響,對于影響一次枝晶臂間距的其他因素研究的較少。

    本文作者首先通過實(shí)驗(yàn)獲得不同凝固速度下Al-7Si-0.36Mg合金的定向凝固組織和枝晶臂間距的統(tǒng)計數(shù)據(jù),然后基于之前建立的CA模型基礎(chǔ)之上[20-21],預(yù)測了一次枝晶臂間距的上限值和下限值,并模擬了該合金定向凝固過程中不同凝固速度下枝晶臂間距演化過程,同時與實(shí)驗(yàn)結(jié)果進(jìn)行對比。此外,還對影響定向凝固一次枝晶臂選擇的諸多因素進(jìn)行了詳細(xì)分析。

    1 實(shí)驗(yàn)

    1.1實(shí)驗(yàn)方法

    本實(shí)驗(yàn)用于定向凝固實(shí)驗(yàn)合金的成分為Al-7Si-0.36Mg(質(zhì)量分?jǐn)?shù),%)。由于該合金不同組元之間的密度差異小,自然對流對枝晶生長的影響相對較小。定向凝固實(shí)驗(yàn)在Bridgman 定向凝固爐中完成。首先,把直徑為10mm,長度為120mm的試棒放入到高純石墨坩堝中,在含有氬氣保護(hù)的氣氛中把試樣加熱到840℃,保溫0.5 h,然后把試樣按設(shè)定的速度往下抽拉到裝有Ga-In-Sn液態(tài)金屬的冷卻罐中,使合金從下往上長出柱狀晶組織。實(shí)驗(yàn)過程中溫度梯度為15 K/mm[22],完成了4種不同抽拉速度下的定向凝固試樣,分別為50、100、150、200μm/s。截取高度在75~85mm區(qū)間的試棒用于金相試樣的制取并在偏振光顯微鏡下觀察縱截面和橫截面的枝晶組織。通過在橫截面上采用截線法來統(tǒng)計各種凝固條件下柱狀晶的一次枝晶臂間距λ1,每種凝固條件都統(tǒng)計200個數(shù)據(jù)點(diǎn),從而可以通過統(tǒng)計的方式得到一次枝晶臂間距λ1的概率分布圖以及一次枝晶臂間距最大值λmax、平均值λave、最小值λmin等特征參數(shù)。

    1.2實(shí)驗(yàn)結(jié)果

    圖1所示為4種不同抽拉速度下凝固試樣橫截面的定向凝固組織。從圖1中可以看出,在一次枝晶臂上發(fā)展出了明顯的二次枝晶臂,在二次枝晶臂上出現(xiàn)了三次枝晶臂。對于Al-7Si-0.36Mg合金來說,該合金的初生相凝固溫度區(qū)間相對較窄,在共晶形成之前,熔體中的枝晶組織沒有足夠的時間進(jìn)行枝晶臂的粗化、合并,所以該合金最終的初生相呈現(xiàn)為比較發(fā)達(dá)的樹枝晶形貌。比較不同抽拉速度下的枝晶組織可以明顯看到,隨著抽拉速度從50增加到150μm/s,枝晶發(fā)生明顯的細(xì)化(見圖1(a)~(c)),而當(dāng)抽拉速度從150增大到200μm/s時,枝晶的細(xì)化效果相對不是很明顯。此外,從圖1中也可以看到在同一凝固條件下,相鄰一次枝晶臂間距也并不是固定的。為了定量表征不同凝固速度下枝晶間距的變化,圖2所示為一次枝晶臂間距的概率分布圖。由圖2可見,一次枝晶臂間距存在較大的穩(wěn)態(tài)區(qū)間,且隨著抽拉速度的增大,一次枝晶臂穩(wěn)態(tài)區(qū)間減小,概率分布更加趨于正態(tài)分布。表1所列為不同凝固條件下一次枝晶臂間距最大值λmax、平均值λave、最小值λmin以及最大值與最小值的比值λmax/λmin。由表1可以看出,λmax/λmin的比值大致分布在3.0左右。許多文獻(xiàn)從實(shí)驗(yàn)和模擬的角度對一次枝晶臂的穩(wěn)態(tài)區(qū)間進(jìn)行了相關(guān)研究。對于高溫合金,QUESTED等[23]、黃太文等[24]、WANG等[13]的研究表明穩(wěn)態(tài)時一次枝晶臂上限值與下限值的比值λmax/λmin接近3。對于丁二腈-2.5%乙醇(質(zhì)量分?jǐn)?shù))透明合金,HUANG等[11]和單博煒等[15]研究結(jié)果顯示,在不同凝固速度下,λmax/λmin的比值在2左右波動,而對于丁二腈-0.055.%丙酮溶液,HAN等[25]研究發(fā)現(xiàn),λmax/λmin的比值約為1.3。由此可見,一次枝晶臂間距上下限的比值與合金體系密切相關(guān)。

    圖1 Al-7Si-0.36Mg合金定向凝固不同凝固速度下柱狀晶橫截面的枝晶形貌Fig.1 Dendrite morphologies on transverse sections of Al-7Si-0.36Mg alloy directional solidified at different growth velocities:(a)50μm/s;(b)100μm/s;(c)150μm/s;(d)200μm/s

    圖2 不同凝固速度下一次枝晶臂間距λ1的概率分布Fig.2 Size distribution of primary dendrite spacing λ1solidified at different growth velocities:(a)50μm/s;(b)100μm/s;(c)150μm/s;(d)200μm/s

    表1 不同凝固速度下Al-7Si-0.36Mg合金一次枝晶臂間距最大值λmax、平均值λave、最小值λmin以及λmax/λmin比值Table 1 λmax, λave, λminand λmax/λminof characteristic primary dendrite spacing at different growth velocities in directional solidification of Al-7Si-0.36Mg alloy

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

    利用CA模型對Al-7Si-0.36Mg合金定向凝固條件下枝晶生長過程開展相應(yīng)數(shù)值模擬研究,模擬采用的算法在文獻(xiàn)[20-21]中進(jìn)行了詳細(xì)描述,合金的熱物性參數(shù)請參見文獻(xiàn)[22]。本文模擬內(nèi)容主要為:首先通過設(shè)置不同初始一次枝晶臂間距來考察不同凝固速度下穩(wěn)態(tài)一次枝晶臂間距的上下限值,然后模擬了圖1中不同凝固條件下的枝晶選擇過程,預(yù)測平均一次枝晶臂間距,并與實(shí)驗(yàn)結(jié)果進(jìn)行對比,最后利用CA模型對影響枝晶臂選擇過程的因素進(jìn)行了分析。

    2.1一次枝晶臂間距λmax和λmin預(yù)測

    由于一次枝晶臂間距的分布區(qū)間決定了定向凝固組織的均勻程度,而一次枝晶臂間距上限值λmax和下限值λmin是用于體現(xiàn)該分布區(qū)間的重要參數(shù),因此準(zhǔn)確預(yù)測不同凝固條件下一次枝晶臂間距的上下限值顯得尤為重要。有研究表明[13, 26],定向凝固一次枝晶臂間距的選擇是一個自組織過程,在固定的溫度梯度G和凝固速度v作用下,柱狀晶一次枝晶臂間距存在一個穩(wěn)態(tài)區(qū)間λmin~λmax(見表1)。當(dāng)枝晶臂間距λ<λmin,由于相鄰枝晶臂之間強(qiáng)烈的溶質(zhì)交互作用,會出現(xiàn)競爭生長,通過淹沒機(jī)制來調(diào)整枝晶臂間距至穩(wěn)態(tài)區(qū)間;當(dāng)枝晶臂間距λ>λmax,枝晶間的溶質(zhì)成分富集較低,過冷度大,有利于通過側(cè)向分叉機(jī)制產(chǎn)生新的枝晶臂來降低枝晶臂間距。為了預(yù)測穩(wěn)態(tài)枝晶臂間距的上下限值,在固定的計算區(qū)域(2.6mm×5mm,網(wǎng)格尺寸5μm)中通過設(shè)置不同的初始晶粒密度來調(diào)整初始一次枝晶臂間距,同時采用與圖1中實(shí)驗(yàn)一致的凝固條件進(jìn)行模擬。

    圖3所示為初始枝晶臂間距在λmax~λmin區(qū)間的枝晶生長情況。以圖3(b1)~(b4)為例,當(dāng)初始枝晶臂間距λini=170μm時(見圖3(b1)),初始一次枝晶臂間距發(fā)生淹沒現(xiàn)象,淘汰了部分的一次枝晶臂,使得一次枝晶臂間距增大。而當(dāng)λini超過180μm時(見圖3(b2)),穩(wěn)態(tài)一次枝晶臂間距等于初始值λini,表明在凝固速度v=100μm/s時,最小一次枝晶臂間距λmin約為180μm。當(dāng)λini=575μm時(見圖3(b3)),凝固過程中二次枝晶臂上形成大量的三次枝晶臂,但是由于相鄰枝晶臂空間上的限制和自身二次枝晶臂的阻礙,這些三次枝晶臂很難發(fā)展成為穩(wěn)定的一次枝晶臂。當(dāng)λini增大到590μm 時(見圖3(b4)),部分三次枝晶臂發(fā)展成了一次枝晶臂,從而降低了相鄰枝晶臂間距,這表明上限值λmax處于575~590μm。

    穩(wěn)態(tài)一次枝晶臂間距的下限值λmin也可以通過Hunt-Lu模型來預(yù)測,即式(2)。該模型并沒有給出上限值的函數(shù)表達(dá)式,而是認(rèn)為上限值λmax約為下限值λmin的2倍。圖4所示為模擬得到的一次枝晶臂間距上下限值、表1中實(shí)驗(yàn)值以及Hunt-Lu模型預(yù)測結(jié)果(λmax=2λmin)之間的比較,從圖4中可以看到,CA模型模擬得到的結(jié)果與實(shí)驗(yàn)結(jié)果非常接近,吻合程度明顯優(yōu)于Hunt-Lu模型的預(yù)測結(jié)果。DING等[26]基于透明合金丁二腈-乙醇溶液定向凝固實(shí)驗(yàn)結(jié)果總結(jié)出一次枝晶臂上限值λmax和下限值λmin相對于凝固速度v呈冪函數(shù)關(guān)系式:λmax=avb、λmin=cvd。將CA模擬結(jié)果進(jìn)行數(shù)據(jù)擬合,得到的方程為:λmax=8837v-0.61、λmin=3495v-0.65,可以看到上限值λmax和下限值λmin基

    本呈3倍的關(guān)系,與實(shí)驗(yàn)結(jié)果吻合良好。

    圖3 不同凝固速度和不同初始一次枝晶臂間距條件下定向生長枝晶陣列Fig.3 Steady state dendrite arrays in directional solidification at different initial primary dendrite arm spacings (λini)and different growth velocities:(a1)-(a4)50μm/s;(b1)-(b4)100μm/s;(c1)-(c4)150μm/s;(d1)-(d4)200μm/s

    圖4 不同凝固速度下一次枝晶臂間距上下限λmax、λmin預(yù)測結(jié)果和實(shí)驗(yàn)結(jié)果比較Fig.4 Comparison of allowable range of primary dendrite spacing between CA simulations, analytical predictions and experimental results (λmaxpredicted by Hunt-Lu model is twice of λmin;λmax,simand λmin, simpredicted by CA model)

    2.2不同凝固速度下枝晶定向生長過程模擬

    圖5 定向凝固條件下模擬得到的不同凝固速度下枝晶穩(wěn)態(tài)生長時的枝晶陣列Fig.5 Simulated steady state dendrite arrays at different growth velocities in directional solidification:(a)200μm/s;(b)150μm/s;(c)100μm/s;(d)50μm/s

    文獻(xiàn)[26]指出,一次枝晶臂間距的選擇不僅取決于穩(wěn)態(tài)凝固參數(shù),而且還與枝晶生長的歷史過程相關(guān)。本研究中采用的是單一凝固條件,即枝晶在整個生長過程中采用恒定的凝固速度和溫度梯度,因此,可以認(rèn)為在本研究中枝晶生長歷史對一次枝晶臂間距的影響較小。為了模擬圖1中4種凝固條件下的枝晶生長過程,選取計算域大小為6.25mm×15mm,CA單元尺寸5μm。初始在計算域底部放置一定數(shù)量的晶核。圖5所示為模擬得到的4種凝固條件下的枝晶定向凝固演化過程。從圖5中可以看到,在凝固初期,從晶核的二次枝晶臂上發(fā)展出了大量的三次枝晶臂,并且部分三次枝晶臂發(fā)展成了穩(wěn)定的一次枝晶臂。這主要是因?yàn)槌跏嫉木Ш碎g距λini遠(yuǎn)大于一次枝晶臂間距上限值λmax,從而出現(xiàn)側(cè)向分枝而形成新的一次枝晶臂,但由于大部分枝晶臂相鄰間距小于一次枝晶臂間距下限值λmin,從而出現(xiàn)淹沒,只有少數(shù)能成為穩(wěn)定的一次枝晶臂。此外,從圖5中也可以看到,凝固不久后系統(tǒng)就可以達(dá)到穩(wěn)定狀態(tài),若沒有外界的擾動,則后續(xù)凝固過程中枝晶可以一直穩(wěn)定生長,枝晶臂間距不變。圖6所示為CA模擬結(jié)果、解析模型預(yù)測結(jié)果以及實(shí)驗(yàn)結(jié)果的比較。從圖6中可以看到,Kurz-Fisher模型預(yù)測的結(jié)果普遍高于實(shí)驗(yàn)結(jié)果,而Hunt模型預(yù)測的結(jié)果則小于實(shí)驗(yàn)結(jié)果,并且兩個模型的預(yù)測結(jié)果與實(shí)驗(yàn)值偏離的程度都較大。Trivedi模型相對較為準(zhǔn)確,但是與CA模擬的結(jié)果相比,誤差還是較大,這表明CA模型在枝晶定向凝固過程枝晶形貌演化模擬和枝晶臂間距預(yù)測等方面的準(zhǔn)確性。對于上述的解析模型,界面能各向異性的影響是沒有進(jìn)行考慮的,而事實(shí)證明,界面能各向異性對于枝晶生長具有較大影響,而CA模型可以很好地把界面能各向異性的影響耦合進(jìn)去,這可能是CA模型預(yù)測準(zhǔn)確性較高的原因之一。此外,Hunt等模型認(rèn)為一次枝晶臂均勻分布,顯然與實(shí)際不相符,這也可能導(dǎo)致預(yù)測結(jié)果不準(zhǔn)確。類似的可以對CA模擬的結(jié)果進(jìn)行數(shù)據(jù)擬合,得到的關(guān)系式為λave=4788v-0.59。

    圖6 不同凝固速度下平均一次枝晶臂間距λave預(yù)測結(jié)果和實(shí)驗(yàn)結(jié)果比較Fig.6 Comparison of analytical predictions and experimental results of average primary dendrite arm spacing λavebetween CA simulations

    2.3影響一次枝晶臂間距因素分析

    對于固定成分的合金,定向凝固枝晶臂間距的選擇過程受外因和內(nèi)因的共同作用。外因主要是抽拉速度v、溫度梯度G等可控因素,而內(nèi)因包括界面能大小、溶質(zhì)擴(kuò)散系數(shù)DL、枝晶生長方向與熱流的偏離角度θ等。外應(yīng)和內(nèi)因的共同作用會改變溫度場和溶質(zhì)場,使枝晶凝固前沿的過冷度ΔT發(fā)生變化,改變枝晶尖端生長速度 (νn∝ΔT2)和界面穩(wěn)定性,從而影響枝晶的生長和一次枝晶臂的選擇過程。

    2.3.1界面能和擴(kuò)散系數(shù)DL的影響

    界面能大小是影響枝晶尖端生長動力學(xué)的一個重要因素,體現(xiàn)了枝晶的分枝能力。LANGER等[27]通過枝晶尖端的線性分析,提出了尖端界面穩(wěn)定性理論,即LMK理論。該理論認(rèn)為枝晶尖端半徑r和尖端生長速度v之間滿足關(guān)系:vr2=2ΓDL/(σ*ΔT0)。式中:Γ為表征界面能的Gibbs-Thomson系數(shù);DL為液相擴(kuò)散系數(shù);ΔT0為結(jié)晶溫度區(qū)間;σ*是尖端穩(wěn)定性參數(shù),約為1/(4π2)。從該模型中可以看到,Γ越小,會促使尖端半徑減小,當(dāng)尖端半徑減小到不能保證界面穩(wěn)定時,就會通過界面分叉的形式來保證生長過程中的界面穩(wěn)定性,從而減小一次枝晶臂間距。枝晶凝固過程中,起到穩(wěn)定固液界面作用的界面能效應(yīng)通過Γζ來體現(xiàn),其中ζ為表征界面能各向異性的函數(shù),表示為ζ=1-15εcos[4(φ-θ)],ε為界面能各向異性系數(shù)。隨著界面各向異性系數(shù)ε的增大,ζ將減小,導(dǎo)致界面能穩(wěn)定作用減弱,造成枝晶臂容易出現(xiàn)側(cè)向分枝,減小一次枝晶臂間距[28-29]。此外,由文獻(xiàn)[30-31]中可以看出,隨著各向異性系數(shù)ε增大,尖端穩(wěn)定性參數(shù)σ*呈冪指數(shù)增大,這也有利于枝晶臂的分枝形成。YIN 等[32]通過CA方法研究了Gibbs-Thomson系數(shù)Γ和界面能各向異性系數(shù)ε對于枝晶形貌的影響,模擬得到的結(jié)果與上述的解釋相一致。從LMK理論還可以看到,溶質(zhì)擴(kuò)散系數(shù)DL減小也會降低枝晶界面穩(wěn)定性,從而影響枝晶生長。DL越小,溶質(zhì)擴(kuò)散能力減弱,相鄰枝晶間的溶質(zhì)場交疊作用降低,促進(jìn)固液界面前沿的過冷度(ΔT=Tliq(CL)-T*,式中:Tliq為成分為CL的液相線溫度;T*為實(shí)際界面溫度)增大,降低枝晶界面穩(wěn)定性,有利于枝晶臂的分枝,降低枝晶臂間距以保證穩(wěn)定生長。表2所列為兩組不同液相溶質(zhì)擴(kuò)散系數(shù)條件下的平均一次枝晶臂間距的模擬結(jié)果。綜上所述,界面能和溶質(zhì)擴(kuò)散系數(shù)的減小都有利于定向凝固過程一次枝晶臂間距的減小。

    表2 液相溶質(zhì)擴(kuò)散系數(shù)DL對一次枝晶臂間距的影響Table 2 Effects of diffusion coefficients on primary dendrite spacing during directional solidification

    2.3.2凝固參數(shù)(v,G)的影響

    定向凝固一次枝晶臂間距主要取決于溫度梯度G和生長速度v這兩個可控的外部因素,許多研究者[24,26]針對這兩個參數(shù)對定向凝固一次枝晶臂間距的影響進(jìn)行了大量實(shí)驗(yàn)研究,表明一次枝晶臂間距隨著溫度梯度和生長速度的增大而減小,這一結(jié)論與式(1)相一致。溫度梯度和生長速度對枝晶臂間距的影響機(jī)制類似,兩者的變化會改變固液界面前沿的散熱效果和溶質(zhì)場分布,從而影響界面前沿的過冷度ΔT,導(dǎo)致枝晶臂間距發(fā)生變化。為了說明這一點(diǎn),圖7所示為初始枝晶臂間距λini=320μm、溫度梯度G=15 K/mm、凝固速度v為100和200μm/s兩種條件下凝固前沿的過冷度分布。從圖7中可以看出,由于相鄰枝晶臂間的“V”形區(qū)域落后于枝晶尖端,使得“V”形區(qū)域的溫度值T*較低,熱過冷度較大,從而使得該處的局部冷度比枝晶尖端大,因此,該處容易發(fā)展出穩(wěn)定的一次枝晶臂來填充該區(qū)域(見圖3)。比較圖7中兩組不同凝固速度下的模擬結(jié)果可以看到,當(dāng)凝固速度從100增加至200μm/s時,由于冷卻速度增大(CR=G·v),所以在凝固至同一高度時,在v=200μm/s條件下的溫度值T*要比在v=100μm/s條件下的低,從而枝晶凝固前沿的過冷度ΔT在v=200μm/s條件下明顯比v=100μm/s時的大。過冷度的增大促進(jìn)枝晶生長速度增大,從而導(dǎo)致一次枝晶臂的形成和一次枝晶臂間距的減小。

    圖7 初始枝晶臂間距λini=320μm、凝固速度v=100μm/s 和v=200μm/s下凝固前沿的過冷度分布Fig.7 Undercooling distributions of solidification front with same initial dendrite spacing (λini)of 320μm at different growth velocities:(a)v=100μm/s;(b)v=200μm/s

    2.3.3生長方向偏離熱流方向角度θ的影響

    實(shí)際定向凝固過程中,枝晶生長方向并不完全平行于熱流方向,兩者往往存在一定角度θ,如圖8所示。從圖8中可以看到,晶粒A和晶粒B的〈100〉方向并不完全沿著溫度梯度G的方向,而是存在一定的夾角,且晶粒B的夾角θB(15°)要大于θA(2°)。通過測定A和B晶粒的一次枝晶臂間距發(fā)現(xiàn),B晶粒的平均一次枝晶臂間距(234μm)要大于A晶粒的平均一次枝晶間距(206μm),這與文獻(xiàn)[33-34]中的實(shí)驗(yàn)結(jié)果是相一致的,即枝晶臂間距會隨著枝晶〈100〉生長方向與溫度梯度方向之間夾角θ的增大而增大。從式(1)和(2)可以看到,HUNT等[5]提出的預(yù)測枝晶臂間距的理論模型中并沒有考慮枝晶生長取向?qū)χчg距的影響,因此,需要進(jìn)行修正才能解釋晶體取向?qū)σ淮沃П坶g距的影響。若一次枝晶臂的生長方向與熱流方向的夾角為θ,則枝晶沿〈100〉方向的生長速度和溫度梯度分別為v/cosθ、Gcosθ,代入式(1)可得

    根據(jù)式(3)可以看到,對于固定成分C0,在恒定溫度梯度G和生長速度v作用下,當(dāng)θ增大時,cosθ減小,從而使得一次臂間距λ1增大。GANDIN等[35]在研究晶體取向?qū)χL的影響時,指出一次枝晶臂間距對θ的依賴關(guān)系主要是由三次枝晶臂和二次枝晶臂相互間的競爭機(jī)制決定的,并基于分枝機(jī)制建立了一次枝晶臂間距預(yù)測模型:

    式中:λ(C0,v,G)、e、d均為正數(shù),從而也可以說明枝晶臂間距隨著擇優(yōu)取向偏離角度θ的增大而增大。

    圖8 生長方向偏離熱流方向的角度θ大小對枝晶生長和枝晶臂間距的影響(G=15 K/mm, v=200μm/s)Fig.8 Effects of deviation angle θ between dendrite growth direction and heat flux direction on dendrite growth and dendrite spacing (G=15 K/mm, v=200μm/s)

    通過模擬研究了θ對一次枝晶臂間距的影響。初始在計算域底部設(shè)置不同生長取向的晶核(5°~30°),且G=15 K/mm、v=200μm/s,并測定穩(wěn)態(tài)生長時的枝晶間距。圖9所示為平均枝晶臂間距隨θ的變化曲線,可以看到λave隨著θ的增大而增大,這與式(3)和(4)的解釋相一致。

    圖9 平均一次枝晶臂間距隨生長取向角度θ的變化曲線(G=15 K/mm, v=200μm/s)Fig.9 Predicted dendrite arm spacing varying with deviation angle θ (G=15 K/mm, v=200μm/s)

    3 結(jié)論

    1)通過開展Al-7Si-0.36Mg合金在不同凝固速度下的定向凝固實(shí)驗(yàn),測得了一次枝晶臂間距最大值λmax、平均值λave、最小值λmin以及λmax/λmin比值。

    2)采用CA模型預(yù)測了Al-7Si-0.36Mg合金不同凝固速度下穩(wěn)態(tài)一次枝晶臂間距的上限值、下限值及平均值,分別可以表示為:λmax=8837v-0.61、λmin=3495v-0.65、λave=4788v-0.59。模擬結(jié)果與實(shí)驗(yàn)結(jié)果的吻合程度明顯優(yōu)于Hunt-Lu等模型的預(yù)測結(jié)果,表明本模型在枝晶定向凝固過程枝晶形貌演化模擬和枝晶臂間距預(yù)測等方面的準(zhǔn)確性。

    3)結(jié)合模擬研究和文獻(xiàn)調(diào)研分析了影響定向凝固一次枝晶臂選擇的因素,包括抽拉速度v、溫度梯度G、界面能大小、溶質(zhì)擴(kuò)散系數(shù)DL、枝晶生長方向與熱流的偏離角度θ等。結(jié)果表明,界面能和溶質(zhì)擴(kuò)散系數(shù)的減小都有利于定向凝固過程一次枝晶臂間距的減小;溫度梯度G和生長速度v對一次枝晶臂間距的影響主要通過影響凝固前沿的過冷度來體現(xiàn);枝晶生長取向偏離熱流方向越大,一次枝晶臂間距越大。

    REFERENCES

    [1] BOETTINGER W J, CORELL S R, GREERA L, KARAM A, KURA W, RAPPAZ M, TRIVEDI R.Solidification microstructures:Recent developments, further directions[J].Acta Materialia, 2000, 48:43-70.

    [2] MELO M L N M, RIZZO E M S, SANTOS R G.Predicting dendrite arm spacing and their effect on microporosity formation in directionally solidified Al-Cu alloy[J].Journal of Materials Science, 2005, 40:1599-1609.

    [3] ROCHA O L, SIQUEIRA C A, GARCIA A.Heat flow parameters affecting dendrite spacing during unsteady-state solidification of Sn-Pb and Al-Cu alloy[J].Metallurgical and Materials Transactions A, 2003, 34:995-1006.

    [4] QUAREMA J M V, SANTOS C A, GARCIA A.Correlation between unsteady-state solidification condition, dendrite spacing and mechanical properties of Al-Cu alloys[J].Metallurgical and Materials Transactions A, 2000, 31:3167-3178.

    [5] HUNT J D.Solidification and casting of metals[M].London:The Metal Society, 1979.

    [6] KURZ W, FISHER D J.Dendrite growth at the limit of stability:tip radius and spacing[J].Acta Metallurgica, 1981, 29:11-20.

    [7] TRIVEDI R.Interdendritic spacing:partⅡ.A comparison of theory and experiment[J].Metallurgical Transactions A, 1984, 15:977-982.

    [8] WARREN J A, LANGER J S.Stability of dendritic spacings[J].Physics Review A, 1990, 42:3518-3525.

    [9] WARREN J A, LANGER J S.Prediction of dendritic spacings in a directional solidification experiment[J].Physics Review E, 1993, 47:2702-2712.

    [10] MA D.Modeling of primary spacing selection in dendrite arrays during directional solidification[J].Metallurgical and Materials Transactions B, 2002, 33:223-233.

    [11] HUANG W D, GENG X G, ZHOU Y H.Primary spacing selection of constrained dendritic growth[J].Journal of Crystal Growth, 1993, 134:105-115.

    [12] HUNT J D, LU S Z.Modeling of cellular-dendritic array growth:spacing and structure predictions[J].Metallurgical and Materials Transactions A, 1996, 27:611-623.

    [13] WANG W, LEE P D, MELEAN M.A model of solidification microstructures in nickel based superalloys:Predicting primary dendrite spacing selection[J].Acta Materialia, 2003, 51:2971-2987.

    [14] DIEPERS H J, MA D, STEINBACH I.History effects during the selection of primary dendrite spacing:Comparison of phase-field simulation with experimental observations[J].Journal of Crystal Growth, 2002, 237/239:149-153.

    [15] 單博煒, 黃衛(wèi)東, 林鑫, 魏雷.元胞自動機(jī)模型模擬枝晶一次間距的選擇[J].金屬學(xué)報, 2008, 44:1042-1050.SHAN Bo-wei, HUANG Wei-dong, LIN Xin, WEI Lei.Dendrite primary spacing selection simulation by the cellular automaton model[J].Acta Metallurgica Sinica, 2008, 44:1042-1050.

    [16] ZHANG Xian-fei, ZHAO Jiu-zhou.Dendritic microstructure formation in a directionally solidified Al-11.6Cu-0.85Mgalloy[J].Journal of Crystal Growth, 2014, 391:52-58.

    [17] 張航, 許慶彥, 史振學(xué), 柳百成.DD6高溫合金定向凝固枝晶生長的數(shù)值模擬研[J].金屬學(xué)報, 2014, 50:345-354.ZHANG Hang, XU Qing-yan, SHI Zhen-xue, LIU Bai-cheng.Numerical simulation of dendrite grain growth of DD6 superalloy during directional solidification process[J].Acta Metallurgica Sinica, 2014, 50:345-354.

    [18] 石玉峰, 許慶彥, 龔銘, 柳白成.定向凝固過程中NH4Cl-H2O 枝晶生長的數(shù)值模[J].金屬學(xué)報, 2011, 47:620-627.SHI Yu-feng, XU Qing-yan, GONG Ming, LIU Bai-cheng.Simulation of NH4Cl-H2O dendritic growth in directional solidification[J].Acta Metallurgica Sinica, 2011, 47:620-627.

    [19] GANDIN C A, RAPPAZ M.A 3D cellular automaton algorithm for the prediction of dendritic grain growth[J].Acta Materials, 1997, 45(5):2187-2195.

    [20] CHEN Rui, XU Qing-yan, LIU Bai-cheng.Cellular automaton simulation of three-dimensional dendrite growth in Al-7Si-Mg ternary aluminum alloys[J].Computational Materials Science, 2015, 105:90-100.

    [21] CHEN Rui, XU Qing-yan, LIU Bai-cheng.A modified cellular automaton model for the quantitative prediction of equiaxed and columnar dendritic growth[J].Journal of Materials Science and Technology, 2014, 30:1311-1320.

    [22] 石玉峰.鋁合金砂型鑄件微觀組織及力學(xué)性能的數(shù)值模擬[D].北京:清華大學(xué), 2013.SHI Yu-feng.Numerical simulation of microstructure and mechanical properties of aluminum alloy sand casting[D].Beijing:Tsinghua University, 2013.

    [23] QUESTED P N, MCLEAN M.Solidification morphologies in directionally solidified superalloys[J].Material Science and Engineering, 1984, 65:171-180.

    [24] 黃太文, 劉林, 張衛(wèi)國, 張軍, 傅恒志.抽拉速度躍遷對定向凝固單晶高溫合金 DD3一次枝晶間距和微觀偏析的影響[J].金屬學(xué)報, 2009, 45(10):1225-1231.HUANG Tai-wen, LIU Lin, ZHANG Wei-guo, ZHANG Jun, FU Heng-zhi.The primary dendrite arm spacing and microsegregation of directionally solidified single crystal superalloy DD3[J].Acta Metallurgica Sinica, 2009, 45(10):1225-1231.

    [25] HAN S H, TRIVEDI R.Primary spacing selection in directionally solidified alloys[J].Acta Metallurgica and Materialia, 1994, 42:25-41.

    [26] DING G L, HUANG W D, HUANG X, LIN X, ZHOU Y H.On primary dendritic spacing during unidirectional solidification[J].Acta Metallurgica, 1996, 44:3705-3709.

    [27] LANGER J S, MüLLER-KRUMBHAAR H.Theory of dendritic growth-Ⅰ.Elements of a stability analysis[J].Acta Metallurgica, 1978, 26:1681-1687.

    [28] 張顯飛, 趙九州.合金定向凝固一次枝晶間距模擬[J].中國有色金屬學(xué)報, 2012, 22(10):2868-2874.ZHANG Xian-fei, ZHAO Jiu-zhou.Simulation of primary dendrite spacing in unidirectionally solidified alloy[J].The Chinese Journal of Nonferrous Metals, 2012, 22(10):2868-2874.

    [29] 張?jiān)迄i, 林鑫, 魏雷, 王猛, 彭東劍, 黃衛(wèi)東.用CA方法模擬界面能各向異性對胞晶生長形態(tài)的影響[J].物理學(xué)報, 2012, 61:228106.ZHANG Yun-peng, LIN Xin, WEI Lei, WANG Meng, PENG Dong-jian, HUANG Wei-dong.Effect of surface tension anisotropy on the growth patterns of cellulars in directional solidification[J].Acta Physics Sinica, 2012, 61:228106.

    [30] KESSLER D A, LEVINE H.Velocity selection in dendritic growth[J].Physics Review B, 1986, 33:7867-7870.

    [31] XU J J.Generalized needle solutions, interfacial instabilities and pattern formations[J].Physics Review E, 1996, 53:5051-5062.

    [32] YIN H, FELICELLI S D, WANG L.Simulation of a dendrite microstructure with the lattice Boltzmann and cellular automaton methods[J].Acta Materialia, 2011, 59:3124-3136.

    [33] 王賢斌, 林鑫, 王理林, 白貝貝, 王猛, 黃衛(wèi)東.晶體取向?qū)Χㄏ蚰讨L的影響[J].物理學(xué)報, 2013, 62(10):108103.WANG Xian-bin, LIN Xin, WANG Li-lin, BAI Bei-bei, WANG-Meng, HUANG Wei-dong.Effect of crystallographic orientation on dendrite growth in directional solidification[J].Acta Physics Sinica, 2013, 62(10):108103

    [34] 楊初斌, 劉林, 趙新寶, 劉剛, 張軍, 傅恒志.〈001〉和〈011〉取向DD407單晶高溫合金枝晶間距和微觀偏析[J].金屬學(xué)報, 2011, 47:1246-1250.YANG Chu-bin, LIU Lin, ZHAO Xin-bao, LIU Gang, ZHANG Jun, FU Heng-zhi.Dendrite arm spacings and microsegregations in 〈001〉 and 〈011〉 orientated single crystal superalloys DD407[J].Acta Metallurgica Sinica, 2011, 47:1246-1250.

    [35] GANDIN C H A, ESHELMAN M, TRIVEDI R.Orientation dependence of primary dendrite spacing[J].Metallurgical and Materials Transactions A, 1996, 27:2727-2739.

    (編輯王超)

    Experimental and simulation of primary dendrite spacing in directional solidification of Al-7Si-0.36Mg alloy

    CHEN Rui, XU Qing-yan, LIU Bai-cheng
    (Key Laboratory for Advanced Materials Processing Technology, Ministry of Education, School of Materials Science and Engineering, Tsinghua University, Beijing 100084, China)

    The directional solidification dendritic morphology evolution and primary dendrite arm spacing selection process of Al-7Si-0.36Mg alloy were studied by directional solidification experiments and cellular automaton model.The results indicate that there is a continuous range of primary dendrite arm spacing under the given solidification condition.Under the condition of the given temperature gradient and various solidification velocities, the relationships between the primary dendrite arm spacing parameters (λmax, λave, λmin)of Al-7Si-0.36Mg alloy and growth velocity were expressed, and λmax/λmin≈3.The simulated results show a quite good agreement with the experimental results, which is better than predicted results of Hunt-Lu model.The comparisons reveal that the present CA model has a high accuracy in simulating the evolution of dendrite morphology and predicting primary dendrite arm spacing in directional solidification.Based on the predictions and related literatures, the factors influencing the selection of primary dendrite arm spacing, such as growth velocity v, temperature gradient G, interfacial energy effect, solute diffusion coefficient DL, as well as the deviation angle θ between the dendrite growth direction, and heat flux direction were analyzed.

    directional solidification;primary dendrite spacing;growth velocity;cellular automaton;numerical simulation

    TG290

    A

    1004-0609(2015)10-2613-10

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃項(xiàng)目(2011CB706801);國家自然科學(xué)基金資助項(xiàng)目(51374137, 51171089)

    2015-01-12;

    2015-05-16

    許慶彥,教授,博士;電話:010-62795482;E-mail:scjxqy@mail.tsinghua.edu.cn

    猜你喜歡
    枝晶限值定向
    關(guān)于廢水排放特別限值的思考
    偏序集上的相對定向集及其應(yīng)用
    定向越野
    遼寧省遼河流域石油煉制排放限值的制定
    A356合金近液相線半固態(tài)鑄造非枝晶組織與模鍛成形性能
    中美煉鋼行業(yè)污染物排放限值研究
    不同形狀橫向限制對枝晶間距影響作用的相場法模擬
    定向馴化篩選耐毒酵母
    環(huán)境保護(hù)部解讀新發(fā)布的大氣污染物特別排放限值
    合金定向凝固一次枝晶間距模擬
    亚洲在线自拍视频| 欧美精品啪啪一区二区三区| 又黄又粗又硬又大视频| 乱人伦中国视频| 国产91精品成人一区二区三区| 精品久久久久久久毛片微露脸| 日韩 欧美 亚洲 中文字幕| 亚洲男人的天堂狠狠| 国产精品亚洲av一区麻豆| 亚洲欧美激情综合另类| 丰满人妻熟妇乱又伦精品不卡| 波多野结衣巨乳人妻| 日韩精品中文字幕看吧| 午夜激情av网站| 午夜精品国产一区二区电影| 精品一品国产午夜福利视频| 亚洲激情在线av| 国产亚洲精品第一综合不卡| 亚洲专区中文字幕在线| 亚洲伊人色综图| 波多野结衣av一区二区av| АⅤ资源中文在线天堂| 色综合站精品国产| 天天添夜夜摸| 深夜精品福利| 日本精品一区二区三区蜜桃| 长腿黑丝高跟| 69精品国产乱码久久久| 伦理电影免费视频| svipshipincom国产片| √禁漫天堂资源中文www| 国产精品秋霞免费鲁丝片| 国产成人一区二区三区免费视频网站| www.熟女人妻精品国产| 国产精品98久久久久久宅男小说| 1024香蕉在线观看| 国产欧美日韩一区二区精品| 一级毛片精品| 久久欧美精品欧美久久欧美| 久久婷婷人人爽人人干人人爱 | 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品一区二区www| 嫩草影视91久久| 国产午夜精品久久久久久| 99在线视频只有这里精品首页| 亚洲欧美日韩高清在线视频| 日韩免费av在线播放| 变态另类丝袜制服| 极品人妻少妇av视频| 久久国产亚洲av麻豆专区| 日韩免费av在线播放| 国产三级黄色录像| 人人妻人人澡人人看| 日日摸夜夜添夜夜添小说| 亚洲 欧美 日韩 在线 免费| 人人妻人人爽人人添夜夜欢视频| 可以在线观看的亚洲视频| 国产激情欧美一区二区| 久久久久亚洲av毛片大全| videosex国产| 国产极品粉嫩免费观看在线| 亚洲一区二区三区不卡视频| 精品一品国产午夜福利视频| a级毛片在线看网站| √禁漫天堂资源中文www| av免费在线观看网站| 男女做爰动态图高潮gif福利片 | 午夜福利影视在线免费观看| av在线播放免费不卡| 亚洲国产精品成人综合色| 亚洲av电影不卡..在线观看| 日韩中文字幕欧美一区二区| 成人免费观看视频高清| 可以在线观看的亚洲视频| 波多野结衣巨乳人妻| 乱人伦中国视频| av在线天堂中文字幕| 国产精品亚洲av一区麻豆| 日韩国内少妇激情av| 亚洲国产精品sss在线观看| 日本免费a在线| 可以免费在线观看a视频的电影网站| 日韩国内少妇激情av| 啦啦啦免费观看视频1| 色播亚洲综合网| 法律面前人人平等表现在哪些方面| 欧美精品亚洲一区二区| 亚洲,欧美精品.| 给我免费播放毛片高清在线观看| 久久国产亚洲av麻豆专区| 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 久久热在线av| 日本五十路高清| 90打野战视频偷拍视频| 亚洲avbb在线观看| 亚洲电影在线观看av| 长腿黑丝高跟| 性欧美人与动物交配| 在线十欧美十亚洲十日本专区| 精品国产一区二区三区四区第35| 91国产中文字幕| 中文字幕久久专区| 一级a爱片免费观看的视频| 日韩av在线大香蕉| 淫妇啪啪啪对白视频| 每晚都被弄得嗷嗷叫到高潮| 免费看美女性在线毛片视频| 中文字幕最新亚洲高清| 亚洲电影在线观看av| 啦啦啦 在线观看视频| 热99re8久久精品国产| 亚洲国产日韩欧美精品在线观看 | 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 国产精品综合久久久久久久免费 | 50天的宝宝边吃奶边哭怎么回事| 男女午夜视频在线观看| 国产成人影院久久av| 精品国产乱子伦一区二区三区| 12—13女人毛片做爰片一| 欧美日韩中文字幕国产精品一区二区三区 | 成人三级做爰电影| 精品国产亚洲在线| 久久中文字幕一级| 天堂影院成人在线观看| 最好的美女福利视频网| 一本久久中文字幕| 国产成人av激情在线播放| 男人舔女人的私密视频| 丝袜美腿诱惑在线| av天堂久久9| 黄色 视频免费看| 国产精品自产拍在线观看55亚洲| 法律面前人人平等表现在哪些方面| 日本黄色视频三级网站网址| 女生性感内裤真人,穿戴方法视频| 性少妇av在线| 在线观看免费视频日本深夜| 极品教师在线免费播放| 91大片在线观看| 999精品在线视频| 亚洲人成77777在线视频| 国产成人欧美在线观看| 不卡一级毛片| 亚洲欧美一区二区三区黑人| 国内精品久久久久精免费| 久久欧美精品欧美久久欧美| 1024视频免费在线观看| 一级毛片高清免费大全| 国产又爽黄色视频| 桃红色精品国产亚洲av| 欧美乱码精品一区二区三区| 亚洲片人在线观看| 真人一进一出gif抽搐免费| 精品卡一卡二卡四卡免费| 国产一区在线观看成人免费| 日本一区二区免费在线视频| 欧美日韩瑟瑟在线播放| 国产精品九九99| 女人精品久久久久毛片| 身体一侧抽搐| 亚洲aⅴ乱码一区二区在线播放 | 一卡2卡三卡四卡精品乱码亚洲| 国产精品免费视频内射| 一区二区日韩欧美中文字幕| 18禁美女被吸乳视频| 久久香蕉激情| 国产亚洲欧美精品永久| 脱女人内裤的视频| 精品国产美女av久久久久小说| 国产亚洲欧美98| 国产精品久久电影中文字幕| 欧美久久黑人一区二区| 久久精品亚洲熟妇少妇任你| 老汉色∧v一级毛片| 99re在线观看精品视频| 久久中文字幕人妻熟女| 一边摸一边抽搐一进一小说| 久久伊人香网站| 国产精品乱码一区二三区的特点 | 可以在线观看的亚洲视频| 久久伊人香网站| 欧美日韩瑟瑟在线播放| 国产成人系列免费观看| 国产精品日韩av在线免费观看 | 亚洲欧美精品综合一区二区三区| 91成人精品电影| 国产精品久久久久久精品电影 | 午夜福利在线观看吧| 国产aⅴ精品一区二区三区波| 九色亚洲精品在线播放| 久久亚洲精品不卡| 久久久久九九精品影院| 欧洲精品卡2卡3卡4卡5卡区| 国产精品1区2区在线观看.| 日韩 欧美 亚洲 中文字幕| 男女下面插进去视频免费观看| 午夜福利成人在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品免费视频一区二区三区| 18禁裸乳无遮挡免费网站照片 | www国产在线视频色| 一本大道久久a久久精品| 好男人电影高清在线观看| 岛国视频午夜一区免费看| 少妇裸体淫交视频免费看高清 | 啦啦啦韩国在线观看视频| 欧美绝顶高潮抽搐喷水| 天天躁夜夜躁狠狠躁躁| 国产精品一区二区三区四区久久 | av在线天堂中文字幕| 桃色一区二区三区在线观看| 欧美日韩黄片免| 免费观看人在逋| 黄片播放在线免费| 国产高清激情床上av| 在线观看日韩欧美| 精品少妇一区二区三区视频日本电影| or卡值多少钱| 国内精品久久久久精免费| 免费少妇av软件| 69精品国产乱码久久久| 欧美日韩精品网址| 亚洲中文av在线| √禁漫天堂资源中文www| 91字幕亚洲| 一级,二级,三级黄色视频| 亚洲中文日韩欧美视频| 国产欧美日韩一区二区三| 亚洲欧美激情综合另类| 一进一出抽搐gif免费好疼| 久久精品国产亚洲av高清一级| 国产99久久九九免费精品| 男人舔女人下体高潮全视频| 日本撒尿小便嘘嘘汇集6| 国产精品亚洲一级av第二区| 亚洲av熟女| 成人欧美大片| 色综合亚洲欧美另类图片| 免费在线观看亚洲国产| 色哟哟哟哟哟哟| 国产成人精品久久二区二区91| 亚洲色图综合在线观看| 欧美黄色片欧美黄色片| 国产精品永久免费网站| 欧美在线黄色| 亚洲免费av在线视频| 国产一级毛片七仙女欲春2 | 成人18禁高潮啪啪吃奶动态图| av中文乱码字幕在线| 成年人黄色毛片网站| 欧美日韩瑟瑟在线播放| videosex国产| 窝窝影院91人妻| 一区二区日韩欧美中文字幕| 精品日产1卡2卡| 午夜福利高清视频| 国产精华一区二区三区| 国语自产精品视频在线第100页| 成人亚洲精品一区在线观看| 99精品在免费线老司机午夜| 亚洲男人天堂网一区| 午夜福利成人在线免费观看| 波多野结衣av一区二区av| 午夜精品久久久久久毛片777| 亚洲色图av天堂| 黄色视频,在线免费观看| 黄色a级毛片大全视频| 国产不卡一卡二| 精品一区二区三区四区五区乱码| 亚洲精品国产区一区二| 久久性视频一级片| 变态另类成人亚洲欧美熟女 | 女同久久另类99精品国产91| 母亲3免费完整高清在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲熟妇中文字幕五十中出| 亚洲精品国产色婷婷电影| 色在线成人网| e午夜精品久久久久久久| 成人av一区二区三区在线看| 一级,二级,三级黄色视频| 国产免费男女视频| 国内毛片毛片毛片毛片毛片| 黄色毛片三级朝国网站| 黄色视频不卡| 可以在线观看毛片的网站| 亚洲成人国产一区在线观看| 黄色丝袜av网址大全| 日本欧美视频一区| 老司机午夜福利在线观看视频| 亚洲av成人不卡在线观看播放网| 免费在线观看黄色视频的| 成人特级黄色片久久久久久久| 99国产精品一区二区蜜桃av| 免费在线观看视频国产中文字幕亚洲| 久久精品人人爽人人爽视色| 麻豆一二三区av精品| 人人妻人人爽人人添夜夜欢视频| 午夜久久久久精精品| 国产在线观看jvid| 久久久国产成人精品二区| 久久人人97超碰香蕉20202| 国产色视频综合| 国产精品国产高清国产av| 国产精品av久久久久免费| 免费在线观看日本一区| 美女高潮喷水抽搐中文字幕| xxx96com| 深夜精品福利| 欧美一级a爱片免费观看看 | 久久精品国产亚洲av香蕉五月| av片东京热男人的天堂| 男女午夜视频在线观看| 国产亚洲欧美98| 男女下面进入的视频免费午夜 | 伦理电影免费视频| 久久精品国产亚洲av香蕉五月| 亚洲第一欧美日韩一区二区三区| 亚洲美女黄片视频| 成人国语在线视频| 人人妻人人澡欧美一区二区 | 在线观看一区二区三区| 亚洲avbb在线观看| 国产高清激情床上av| 91麻豆av在线| 在线观看一区二区三区| 久久久久国产一级毛片高清牌| 97碰自拍视频| 一级作爱视频免费观看| 国产高清视频在线播放一区| 给我免费播放毛片高清在线观看| 不卡一级毛片| 午夜精品久久久久久毛片777| 动漫黄色视频在线观看| 欧美乱妇无乱码| 国产精品免费视频内射| 色播亚洲综合网| www日本在线高清视频| 大陆偷拍与自拍| 欧美最黄视频在线播放免费| 女生性感内裤真人,穿戴方法视频| 精品国产乱子伦一区二区三区| 黄片小视频在线播放| 欧美国产日韩亚洲一区| 少妇被粗大的猛进出69影院| 国产片内射在线| 一区二区日韩欧美中文字幕| 亚洲精品国产色婷婷电影| 婷婷六月久久综合丁香| 欧洲精品卡2卡3卡4卡5卡区| 国产精品秋霞免费鲁丝片| 国产成人欧美| 电影成人av| 亚洲成国产人片在线观看| 性少妇av在线| av电影中文网址| 午夜福利一区二区在线看| 在线观看免费日韩欧美大片| 咕卡用的链子| 久久狼人影院| 日日干狠狠操夜夜爽| 50天的宝宝边吃奶边哭怎么回事| 国产高清视频在线播放一区| 嫩草影院精品99| 正在播放国产对白刺激| 久久婷婷成人综合色麻豆| 国产精品自产拍在线观看55亚洲| 午夜福利成人在线免费观看| 亚洲欧美激情在线| 90打野战视频偷拍视频| 日本三级黄在线观看| 久久人人97超碰香蕉20202| av超薄肉色丝袜交足视频| 我的亚洲天堂| 黑人巨大精品欧美一区二区蜜桃| 中国美女看黄片| 91在线观看av| 国产91精品成人一区二区三区| 69精品国产乱码久久久| 日日摸夜夜添夜夜添小说| 久久欧美精品欧美久久欧美| 亚洲av日韩精品久久久久久密| 久久这里只有精品19| netflix在线观看网站| 精品一品国产午夜福利视频| 亚洲成av人片免费观看| 亚洲午夜理论影院| 日韩欧美一区二区三区在线观看| 热99re8久久精品国产| av免费在线观看网站| 亚洲av电影在线进入| 最新在线观看一区二区三区| 国产精品国产高清国产av| 婷婷六月久久综合丁香| 视频在线观看一区二区三区| 啦啦啦韩国在线观看视频| 国产野战对白在线观看| 嫁个100分男人电影在线观看| 看黄色毛片网站| 国产91精品成人一区二区三区| 国产成人欧美在线观看| 免费在线观看影片大全网站| 久久天躁狠狠躁夜夜2o2o| 搡老熟女国产l中国老女人| 久久久水蜜桃国产精品网| 欧洲精品卡2卡3卡4卡5卡区| av中文乱码字幕在线| 好男人在线观看高清免费视频 | 久久久久国产一级毛片高清牌| 黄色成人免费大全| 在线观看一区二区三区| 久久精品91无色码中文字幕| 窝窝影院91人妻| 在线观看免费视频网站a站| 国产精品美女特级片免费视频播放器 | 久久久久精品国产欧美久久久| 欧美日本亚洲视频在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲专区国产一区二区| 日韩欧美免费精品| 波多野结衣巨乳人妻| av视频免费观看在线观看| 如日韩欧美国产精品一区二区三区| 国产激情欧美一区二区| 国产又爽黄色视频| 在线观看午夜福利视频| 久久 成人 亚洲| 亚洲中文日韩欧美视频| 韩国av一区二区三区四区| 日本精品一区二区三区蜜桃| 国产成人啪精品午夜网站| 男女午夜视频在线观看| 人人妻人人澡人人看| 麻豆久久精品国产亚洲av| av网站免费在线观看视频| 男男h啪啪无遮挡| 亚洲片人在线观看| 国产成人av教育| 一二三四在线观看免费中文在| 午夜福利,免费看| 亚洲一区高清亚洲精品| 一边摸一边抽搐一进一出视频| 90打野战视频偷拍视频| 国产精品一区二区免费欧美| 国产精品香港三级国产av潘金莲| 亚洲一区高清亚洲精品| 好男人电影高清在线观看| 久久精品成人免费网站| 免费在线观看日本一区| av天堂在线播放| 日韩av在线大香蕉| 精品国内亚洲2022精品成人| 波多野结衣av一区二区av| 午夜免费激情av| 757午夜福利合集在线观看| 国产三级黄色录像| 成人亚洲精品一区在线观看| av天堂久久9| 美女免费视频网站| 久久九九热精品免费| 中文字幕av电影在线播放| 久久国产精品男人的天堂亚洲| 久久久久国内视频| 精品久久久精品久久久| 90打野战视频偷拍视频| 国产成人免费无遮挡视频| 国产不卡一卡二| 精品午夜福利视频在线观看一区| 波多野结衣av一区二区av| 国产亚洲欧美98| www.自偷自拍.com| 在线观看免费午夜福利视频| 色哟哟哟哟哟哟| 后天国语完整版免费观看| 精品一区二区三区视频在线观看免费| 午夜福利在线观看吧| 亚洲精品一卡2卡三卡4卡5卡| 满18在线观看网站| 免费在线观看视频国产中文字幕亚洲| 19禁男女啪啪无遮挡网站| 欧美大码av| 精品国产美女av久久久久小说| 熟女少妇亚洲综合色aaa.| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 久久国产精品男人的天堂亚洲| www.精华液| 久久精品aⅴ一区二区三区四区| 成人亚洲精品av一区二区| 男女午夜视频在线观看| 久久久国产成人免费| 一边摸一边做爽爽视频免费| a级毛片在线看网站| 亚洲一区二区三区不卡视频| 啦啦啦 在线观看视频| 午夜成年电影在线免费观看| 女人高潮潮喷娇喘18禁视频| 波多野结衣一区麻豆| 老熟妇乱子伦视频在线观看| av在线天堂中文字幕| 亚洲色图综合在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产色婷婷电影| 久热这里只有精品99| 老司机午夜福利在线观看视频| 少妇被粗大的猛进出69影院| 精品高清国产在线一区| 欧美日韩中文字幕国产精品一区二区三区 | 午夜精品在线福利| 欧美精品亚洲一区二区| 俄罗斯特黄特色一大片| 久久国产亚洲av麻豆专区| 亚洲精品美女久久久久99蜜臀| 亚洲精品国产一区二区精华液| 成人三级黄色视频| 十分钟在线观看高清视频www| 男女午夜视频在线观看| 操美女的视频在线观看| 国产av在哪里看| 99久久国产精品久久久| 亚洲电影在线观看av| 国产精品 国内视频| 久久久久国产精品人妻aⅴ院| 天天添夜夜摸| 免费一级毛片在线播放高清视频 | 久久久久国产一级毛片高清牌| 国产精品久久久av美女十八| 成年女人毛片免费观看观看9| 丁香欧美五月| 国产日韩一区二区三区精品不卡| 久久久久久大精品| 国产亚洲精品一区二区www| 久久中文字幕一级| 又黄又粗又硬又大视频| 嫩草影院精品99| 亚洲午夜理论影院| 亚洲中文av在线| 亚洲精品av麻豆狂野| 美女高潮喷水抽搐中文字幕| 伦理电影免费视频| 国产精品久久久人人做人人爽| 午夜久久久在线观看| 巨乳人妻的诱惑在线观看| 欧美丝袜亚洲另类 | 99精品久久久久人妻精品| 亚洲国产看品久久| 窝窝影院91人妻| 夜夜爽天天搞| 男人舔女人的私密视频| 亚洲第一电影网av| 亚洲狠狠婷婷综合久久图片| 又黄又粗又硬又大视频| а√天堂www在线а√下载| 久久久久国产精品人妻aⅴ院| 亚洲av第一区精品v没综合| 国产亚洲精品av在线| 亚洲av日韩精品久久久久久密| 一级毛片高清免费大全| www.999成人在线观看| 我的亚洲天堂| 亚洲精品粉嫩美女一区| 国产亚洲精品综合一区在线观看 | 中文字幕高清在线视频| 男女下面插进去视频免费观看| 国产精品综合久久久久久久免费 | 欧美成人午夜精品| 后天国语完整版免费观看| 免费在线观看影片大全网站| 久久九九热精品免费| 亚洲电影在线观看av| 国产精华一区二区三区| 亚洲精品中文字幕在线视频| 成人手机av| 久久婷婷成人综合色麻豆| 黄色片一级片一级黄色片| www.自偷自拍.com| 亚洲欧美激情综合另类| 亚洲成人精品中文字幕电影| 精品久久久久久成人av| 国产三级黄色录像| 在线观看66精品国产| 久久久国产成人免费| 黑人巨大精品欧美一区二区蜜桃| 99久久精品国产亚洲精品| 无遮挡黄片免费观看| 少妇的丰满在线观看| 免费少妇av软件| 欧洲精品卡2卡3卡4卡5卡区| 国产在线观看jvid| 欧美日韩亚洲国产一区二区在线观看| 一区二区日韩欧美中文字幕| 亚洲精品国产色婷婷电影| 国产精品二区激情视频| 男男h啪啪无遮挡| 性色av乱码一区二区三区2| 国产亚洲精品久久久久久毛片| 精品久久蜜臀av无| 在线免费观看的www视频| 国产精品电影一区二区三区| 亚洲狠狠婷婷综合久久图片| 91在线观看av| 亚洲精品中文字幕一二三四区| 9热在线视频观看99| 高清黄色对白视频在线免费看| 色老头精品视频在线观看| 大型av网站在线播放| 麻豆成人av在线观看| 99riav亚洲国产免费| 亚洲精品中文字幕在线视频| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品sss在线观看| 黄色毛片三级朝国网站| 午夜免费激情av|