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

    高性能銅合金熱力學(xué)數(shù)據(jù)庫的開發(fā)及其在材料設(shè)計中的應(yīng)用

    2011-11-24 08:37:56劉興軍王翠萍甘世溪大沼郁雄貝沼亮介石田清仁
    中國有色金屬學(xué)報 2011年10期
    關(guān)鍵詞:銅合金基合金熱力學(xué)

    劉興軍,王翠萍,甘世溪,大沼郁雄,貝沼亮介,石田清仁

    (1. 廈門大學(xué) 材料學(xué)院,廈門 361005;2. 廈門大學(xué) 圣諾有色金屬研究院,廈門 361005;3. 日本東北大學(xué) 工學(xué)部 材料科學(xué)系,仙臺 980-8579)

    高性能銅合金熱力學(xué)數(shù)據(jù)庫的開發(fā)及其在材料設(shè)計中的應(yīng)用

    劉興軍1,2,王翠萍1,2,甘世溪1,2,大沼郁雄3,貝沼亮介3,石田清仁3

    (1. 廈門大學(xué) 材料學(xué)院,廈門 361005;2. 廈門大學(xué) 圣諾有色金屬研究院,廈門 361005;3. 日本東北大學(xué) 工學(xué)部 材料科學(xué)系,仙臺 980-8579)

    利用相圖計算(CALPHAD)方法,采用亞規(guī)則溶體模型描述溶體相的吉布斯自由能,采用亞點陣模型描述金屬間化合物和有序相的吉布斯自由能,并結(jié)合相平衡和熱力學(xué)性質(zhì)的實驗數(shù)據(jù),優(yōu)化與計算Cu-X二元系以及Cu-Fe、Cu-Ni、Cu-Cr、Cu-Co、Cu-Mo和Cu-W基各三元系的相圖,獲得自洽性良好的熱力學(xué)參數(shù),并建立銅合金熱力學(xué)數(shù)據(jù)庫。該數(shù)據(jù)庫可以提供穩(wěn)定和亞穩(wěn)的相圖計算、相分數(shù)計算、液相面計算、熱力學(xué)性質(zhì)的計算等多種信息,為外推計算銅基多元合金系的相平衡提供理論基礎(chǔ),并為高性能銅合金材料的設(shè)計及制備提供重要的理論指導(dǎo)。

    銅合金;相圖;熱力學(xué)計算;材料設(shè)計

    銅合金由于具有優(yōu)良的導(dǎo)電和導(dǎo)熱性能,以及較高的強度和良好的塑性等優(yōu)良綜合性能,作為功能材料和結(jié)構(gòu)材料已廣泛應(yīng)用于電子、電器和電子封裝等眾多領(lǐng)域,如電阻焊電極、開關(guān)觸橋、電子電路中的基板及導(dǎo)線材料等。近年來,隨著電子產(chǎn)品的小型化和集成化,將對銅合金性能提出更高的要求,以滿足其產(chǎn)品的設(shè)計要求。通常新型銅合金的開發(fā)采用“試錯法”(Trial and error)來實現(xiàn),即材料開發(fā)人員通過大量的實驗和經(jīng)驗來選擇材料的成分、穩(wěn)定工藝參數(shù)等。這樣既消耗了大量的人力和物力,又不利于系統(tǒng)地探討材料改性的機理。因此,20世紀后期,材料科學(xué)家提出了帶有預(yù)測性的材料設(shè)計理念,其目的在于“通過理論與計算來預(yù)測新材料的組分、結(jié)構(gòu)和性能”以及“通過理論來‘訂做’具有特定功能的新材料”。

    相圖作為材料設(shè)計的“地圖”,對新型銅合金材料的開發(fā)具有重要的指導(dǎo)意義。相圖計算的 CALPHAD技術(shù)[1]是通過建立熱力學(xué)模型來計算體系的相圖和熱力學(xué)性質(zhì),其特點是通過二元和三元等低組元系的實驗數(shù)據(jù)為主建立的熱力學(xué)模型和多元系的少量關(guān)鍵實驗數(shù)據(jù)相結(jié)合,可以預(yù)測實用多元合金的相平衡性質(zhì)、熱力學(xué)性質(zhì)、組元的活度和蒸汽壓、相變驅(qū)動力等合金的性質(zhì)等。CALPHAD 方法已經(jīng)成為多元合金設(shè)計的有效手段,并且已廣泛應(yīng)用于新型材料的設(shè)計和開發(fā)[2],為了開發(fā)高性能銅基合金材料,建立銅基合金相圖的熱力學(xué)數(shù)據(jù)庫實現(xiàn)銅基合金的成分與組織的精確設(shè)計,將是一項具有重要理論價值的研究工作。

    本研究利用相圖計算(CALPHAD)方法,結(jié)合各種相平衡和熱力學(xué)性質(zhì)的實驗數(shù)據(jù),建立 Cu-X二元系以及Cu-Fe、Cu-Ni、Cu-Cr、Cu-Co、Cu-Mo和Cu-W基各三元系的熱力學(xué)數(shù)據(jù)庫,為新型銅基合金的成分和組織設(shè)計提供重要的理論指導(dǎo),同時也為銅基合金的動力學(xué)計算提供基礎(chǔ)的熱力學(xué)參數(shù)。

    1 相圖計算熱力學(xué)模型

    在進行相圖計算和熱力學(xué)評估時,首先要選擇合理的熱力學(xué)模型。下面分別就本研究中所涉及到的純組元、液相和端際固溶體相、化學(xué)計量比化合物、金屬間化合物溶體相、有序相等所選用的熱力學(xué)模型進行闡述。

    1.1 純組元

    對于純組元,由于Gibbs自由能的絕對值無法確定,但選定參考態(tài)(通常以298.15 K和0.1 MPa大氣壓下純固態(tài)元素穩(wěn)定態(tài)的摩爾始值為參考態(tài)),即穩(wěn)定元素參考態(tài)(STableelement reference, SER),通常以(298.15 K)表示。其相對于參考態(tài)的差值(ΔG)具有重要的物理意義,則任一元素 i以φ相存在時的Gibbs自由能就可以表示為

    通常,恒壓熱容Cp與溫度之間的關(guān)系可表示為

    恒壓下,根據(jù)熱力學(xué)函數(shù)之間的關(guān)系式,Gibbs自由能和溫度的關(guān)系可用下式表示:

    事實上,式(2)僅適用于一個有限的溫度范圍,為避免式(3)中的參數(shù)過多,通常采用劃分溫度區(qū)間的方法。各溫度區(qū)間內(nèi)Cp和Gibbs自由能的表達式的形式相同,但參數(shù)值不同。而在溫度間斷處,Gibbs自由能和Cp值都被規(guī)定為連續(xù)的,并且溫度區(qū)間的劃分須有理論根據(jù),不能隨意劃分。

    本研究中,涉及到的純組元的 Gibbs自由能函數(shù)均取自Dinsdale[3]評估和優(yōu)化的SGTE元素數(shù)據(jù)庫。

    1.2 液相和端際固溶體相

    液相和端際固溶體相都是采用亞正規(guī)溶體模型描述的,其摩爾Gibbs自由能表示為

    式中:xi為液相中組分i的摩爾分數(shù);為純組分i的液相摩爾Gibbs自由能;為液相的摩爾過剩自由能。

    在三元系中,通過其3個子二元系計算三元系的計算,并基于三元系的實驗數(shù)據(jù),估出三元相互作用參數(shù)(),三元系的過剩自由能()如下所示:

    magGφ是磁性對φ相自由能的貢獻部分,本研究采用 HILLERT和 JARL[5]提出的磁性轉(zhuǎn)變自由能模型,可用如下多項式來表示:

    式中: βφ是與總磁熵有關(guān)的量,大多數(shù)情況下,它等于每摩爾原子的玻爾磁矩; τφ定義為 T /是φ相磁有序的臨界溫度,對于鐵磁性的相為居里溫度Tc,對于反鐵磁性的相為奈爾溫度TN。它們可分別表示為

    函數(shù) f (τφ) 可用下面的式子表示:

    本研究中,液相的和端際固溶體相的摩爾自由能均采用上述亞正規(guī)溶體模型描述。

    1.3 化學(xué)計量比化合物

    對于A-B 二元系中化學(xué)計量比化合物(AaBb),本研究中將它劃分為兩個亞點陣,A和B元素各占據(jù)其一。假設(shè)純組元A和B分別以純元素穩(wěn)定相為參考態(tài),則AaBb的摩爾生成Gibbs自由能可以表示為

    式中的a′′、b′′和c′′的值可以基于實驗數(shù)據(jù)優(yōu)化得到。

    1.4 金屬間化合物溶體相

    在某些體系中,部分化合物具有一定的固溶度,這類相的熱力學(xué)模型可以用雙亞點陣模型來描述。亞點陣模型首先是由 HILLERT和 STAFANSSON[6]于1970年提出的,該模型是把溶體看成是由多個亞點陣組成,固溶體的混合熵等于各亞點陣的混合熵之和。它主要用于處理合金體系中的化學(xué)計量比化合物、置換式固溶體、間隙式固溶體和離子溶體等。

    經(jīng)過30多年的發(fā)展,亞點陣模型已經(jīng)成為物理意義明確、普適性好且含有多個亞點陣的多組分溶體理論模型。其中由Hillert發(fā)展的雙亞點陣(Two sublattice)模型具有較強的普適性。即把溶體看成由兩類亞點陣組成的一種“化合物”,用(M1, M2, …)a(N1, N2, …)c表示,這里(M1, M2, …)a是一個亞點陣,(N1, N2, …)c是另一個亞點陣,a和c是兩個亞點陣的結(jié)點數(shù)的比例,并認為M型原子與N型原子在各自亞點陣結(jié)點的分布是隨機的。以模型(M1, M2)a(N1, N2)c為例,其摩爾Gibbs自由能表達式如下:

    G?:?表示當(dāng)該相中第一個亞點陣被一種組元 M1或 M2占據(jù),第二個亞點陣被組元N1或N2占據(jù)時化合物的Gibbs自由能。 LM1,M2:?, L?:N1,N2中的冒號用來分隔不同點陣中的元素,逗號分隔同一亞點陣中相互作用的元素。 LM1,M2:?表示當(dāng)?shù)谝粋€亞點陣被組元 M1和 M2占據(jù),第二個亞點陣被組元N1或N2占據(jù)時第一個亞點陣中組元M1和M2之間的相互作用參數(shù);L?:N1,N2表示當(dāng)?shù)谝粋€亞點陣被組元M1或M2占據(jù),第二個亞點陣被組元 N1和N2占據(jù)時,第二個亞點陣中N1和 N2組元之間的相互作用參數(shù)。其中,xi表示合金元素的摩爾分數(shù),分別表示組元M1或M2在第一個亞點陣及N1或N2在第二個亞點陣中所占的摩爾分數(shù)。合金元素的摩爾分數(shù)與點陣分數(shù)滿足: xM1+ xM2+= 1 且= 1 。

    1.5 bcc(B2)有序相

    本研究中的 bcc相的無序-有序轉(zhuǎn)變,采用了HILLERT[6]發(fā)展的雙亞點陣(Two sublattice)模型來描述bcc(B2)的Gibbs自由能,其摩爾Gibbs自由能表達式如下:

    2 銅基合金熱力學(xué)數(shù)據(jù)庫的建立

    近年來, 本研究組在銅基合金相圖的優(yōu)化和熱力學(xué)計算方面做了大量的研究工作[7-10],利用相圖計算的CALPHAD方法,結(jié)合相平衡和熱力學(xué)性質(zhì)的相關(guān)實驗數(shù)據(jù),對Cu-X二元系以及Cu-Fe、Cu-Ni、Cu-Cr、Cu-Co、Cu-Mo和Cu-W基各三元系及銅基多元合金相圖進行了熱力學(xué)優(yōu)化與計算,建立了銅基合金的熱力學(xué)數(shù)據(jù)庫,所包含的合金體系如表1所列。從表1可見,該數(shù)據(jù)庫包含了 Cu-X二元系、重要的銅基三元合金系以及多元體系,所含信息量龐大。本研究還采用合金法、擴散偶法、XRD衍射技術(shù)、電子探針微區(qū)成分分析(EPMA)以及DSC差示掃描量熱儀等方法獲得包括一級相變和一些重要體系的有序-無序轉(zhuǎn)變等的相平衡實驗信息[11-17],例如,在Cu-Al[15]、Cu-Sn[16]和 Cu-Fe-Al[17]體系中的 BCC有序相的確定。從而保證了該熱力學(xué)數(shù)據(jù)庫建立在大量的實驗信息的基礎(chǔ)上。該數(shù)據(jù)庫的熱力學(xué)優(yōu)化與計算工作是在 SUNDMAN等[1]開發(fā)的 Thermo-Calc軟件上完成的。其中液相和端際固溶體相的Gibbs自由能采用亞規(guī)則溶體模型來描述,金屬間化合物相和BCC有序相的Gibbs自由能采用亞點陣模型來描述,而氣相的Gibbs自由能采用理想氣體模型來描述。液相熱力學(xué)參數(shù)的評估主要依據(jù)相平衡的實驗數(shù)據(jù)和液相混合焓、液相自由能、活度等熱力學(xué)性質(zhì)的實驗數(shù)據(jù),而固相熱力學(xué)參數(shù)的評估則主要依據(jù)于相平衡的實驗信息和熵、焓等熱力學(xué)性質(zhì)的實驗結(jié)果。計算結(jié)果與實驗值取得了良好的一致性,獲得了合理地描述 Cu基二元、三元系中各相自由能的熱力學(xué)參數(shù)。目前該數(shù)據(jù)庫可以為銅基合金的設(shè)計提供重要的基礎(chǔ)理論參考,如圖1所示。

    表1 銅基合金熱力學(xué)數(shù)據(jù)庫Table1 Thermodynamic database for Cu-base Alloy systems

    圖1 Cu基合金熱力學(xué)數(shù)據(jù)庫的評估和應(yīng)用Fig.1 Outline of assessment and application of thermodynamic database of Cu-base alloy systems

    3 銅基合金熱力學(xué)數(shù)據(jù)庫的應(yīng)用

    本研究組的最終目標(biāo)是建立完善的銅基合金的熱力學(xué)設(shè)計系統(tǒng),達到高效地設(shè)計新型銅基合金的目的。目前所建立的銅基合金的熱力學(xué)數(shù)據(jù)庫可以提供相平衡及熱力學(xué)性質(zhì)等多種信息,例如平衡和亞穩(wěn)相圖的計算、熱力學(xué)性質(zhì)(生成焓、Gibbs自由能、活度、熵等)的計算、相分數(shù)與體積分數(shù)的計算等。該熱力學(xué)數(shù)據(jù)庫將作為銅基合金的熱力學(xué)設(shè)計系統(tǒng)的重要組成部分,為新型高性能銅基合金的設(shè)計提供高效的理論指導(dǎo)。

    圖2 Cu-X二元系的計算相圖Fig.2 Calculated phase diagram in Cu-X binary system

    3.1 相圖的熱力學(xué)計算

    3.1.1 二元相圖的計算

    本研究已經(jīng)建立了較完整的銅基二元合金系相圖的熱力學(xué)數(shù)據(jù)庫。利用該數(shù)據(jù)庫,計算的部分銅基二元系相圖如圖2(a)~(d)所示。圖2(a)所示為Cu-Fe二元系相圖的計算結(jié)果。從計算相圖中可以看出,在低溫區(qū)域,Cu和Fe的相互溶解度很小。而在高溫區(qū)域,存在著亞穩(wěn)態(tài)液相的溶解度間隙(或兩相分離)。圖2(b)所示為Cu-Ni二元系的計算相圖。從圖2(b)可以看出,在全成分范圍內(nèi),Cu和Ni能形成連續(xù)固溶體,并在低溫區(qū)域存在著FCC相的溶解度間隙。圖2(c)所示為Cu-Sn二元系的計算相圖,該體系中存在BCC相(β相)的有序-無序轉(zhuǎn)變以及多個線性化合物相,這些復(fù)雜的相平衡關(guān)系同樣可以實現(xiàn)精確計算。圖2(d)所示為Cu–Be二元系相圖的計算結(jié)果。從圖2(d)中可以看出,隨著溫度的降低,Be在FCC(Cu)中的固溶度減小,富銅側(cè)的 γ(B2)相為 BCC有序相,在高溫區(qū)存在β(A2)+γ(B2)相的兩相平衡區(qū)域,在低溫區(qū)存在著FCC(Cu)+γ(B2)相的兩相區(qū)。

    3.1.2 三元相圖的計算

    利用銅基合金熱力學(xué)數(shù)據(jù)庫,計算的部分銅基三元系相圖如圖 3所示。圖 3(a)~(d)所示為計算的Cu-Fe-(Cr, Co, V, Ni) 各三元系在1 473 K時的等溫截面相圖與實驗數(shù)據(jù)的比較??梢钥闯?,計算結(jié)果與實驗結(jié)果符合的很好。圖4(a)~(b)所示為Cu-Fe-(Co, Mo)各三元系縱截面相圖的計算結(jié)果,計算結(jié)果很好地再現(xiàn)了亞穩(wěn)液相的溶解度間隙的實驗結(jié)果。圖5所示為Cu-Fe-(Co, Mn, V)各三元系中穩(wěn)定和亞穩(wěn)液相的溶解度間隙的計算結(jié)果。這些信息為利用液相的溶解度間隙設(shè)計和制備新型的銅基合金提供重要的理論指導(dǎo)。

    圖3 計算的Cu-Fe-Cr, Cu-Fe-Co, Cu-Fe-V和Cu-Fe-Ni三元系在1 473 K時的等溫截面相圖Fig.3 Calculated isothermal section diagrams at 1 473 K for Cu-Fe-Cr(a), Cu-Fe-Co(b), Cu-Fe-V(c) and Cu-Fe-Ni(d) systems

    3.2 銅基合金熱力學(xué)數(shù)據(jù)庫的應(yīng)用

    3.2.1 銅基合金的固溶度的計算

    圖4 Fe-Cu-Mo三元系縱截面相圖的計算結(jié)果Fig.4 Calculated results of vertical section diagrams including sTableand metasTablemiscibility gaps of liquid phase in Fe-Cu-Mo and Fe-Cu-Co systems: (a) 10%Mo; (b) 10%Co

    圖5 計算的 Fe-Cu-V, Fe-Cu-Co和Fe-Cu-Mn三元系中液相的溶解度間隙Fig.5 Calculated miscibility gap of liquid phase in Fe-Cu-V(a), Fe-Cu-Co(b)and Fe-Cu-Mn(c) systems

    鈹青銅合金是典型的析出強化型銅合金,具有良好的導(dǎo)電導(dǎo)熱性能、較高的強度以及較好的機械加工性能,可以廣泛應(yīng)用于電子材料、控制軸承、電器材料零件以及焊接材料設(shè)備等領(lǐng)域。由于鈹?shù)难趸锞哂卸拘?,為了減少工業(yè)鈹青銅合金中鈹?shù)暮?,通過少量添加合金元素來減小Be在FCC(Cu)中的固溶度,增加 γ(B2)相的析出量以提高鈹青銅合金的強度是鈹青銅合金設(shè)計的主要方向。圖6所示為Ni添加量對Be在 FCC(Cu)相中的固溶度的影響。通過與 Cu-Be二元相圖中的固溶度線相比可知,在FCC(Cu)相中,隨著Ni含量的增加,Be在Cu中的固溶度急劇減少,表明少量Be和Ni的復(fù)合添加可能在Cu-Be-Ni三元系中獲得更顯著的析出強化效果。

    3.2.2 Cu-Ni基合金中FCC相的溶解度間隙的計算

    FCC相的失穩(wěn)分解型銅合金也是新型高強銅合金實用材料的重點研究對象,因此,正確地把握合金元素對溶解度間隙的影響對于失穩(wěn)分解型銅合金的組織設(shè)計是非常必要的。目前在工業(yè)中應(yīng)用的有Cu-30Ni-2.8Cr合金[18],將該合金在 900~1 000 ℃保溫,然后在450~760 ℃溫度區(qū)間慢冷,可以得到FCC相的失穩(wěn)分解組織及最佳的綜合性能。此外,Cu-Ni-Sn三元系中也存在著FCC相的失穩(wěn)分解反應(yīng),但由于反應(yīng)溫度低,基于實驗測定很難準確地確定其相邊界。基于本研究中的熱力學(xué)數(shù)據(jù)庫,可以很容易地計算出FCC相的溶解度間隙,其計算結(jié)果如圖7所示。由圖7可見,隨著Sn的添加量的增加,F(xiàn)CC相的溶解度間隙的溫度表示出有明顯的上升,表明可以在更高溫度下獲得FCC相的失穩(wěn)分解組織。

    3.2.3 核/殼結(jié)構(gòu)銅鐵復(fù)合粉體的設(shè)計與制備

    目前,銅及銅合金粉末的產(chǎn)量僅次于鐵粉、鋁粉,主要用于粉末冶金制品、金剛石工具、電碳、顏料等行業(yè),在工業(yè)生產(chǎn)中發(fā)揮著重要作用。鐵、銅粉的燒結(jié)體和鐵、銅、石墨的混合粉末燒結(jié)體被廣泛應(yīng)用于汽車、電子和工業(yè)機械領(lǐng)域。但是由于石墨和銅粉的膨脹而帶來體積空間的變化以及鐵和銅的潤濕性不好,導(dǎo)致了燒結(jié)體的成分偏析和組織不均勻性等問題[19]。近年來,本文作者巧妙地利用液相溶解度間隙的性質(zhì),利用霧化法制粉工藝(粉體內(nèi)外存在溫度梯度而導(dǎo)致液相界面能的差別作為驅(qū)動力),在通常重力場的條件下,不需要任何復(fù)合工藝,在世界上率先制備出具有核/殼結(jié)構(gòu)的銅鐵復(fù)合粉體[20]。

    圖 8所示為通過利用熱力學(xué)數(shù)據(jù)庫計算的(Fe-7.2Cr)-(Cu-7.2Cr)的縱截面相圖,圖中可見存在著穩(wěn)定的液相的溶解度間隙。圖 9所示為 Cu-32.4Fe-7.2Cr-0.4C(質(zhì)量分數(shù),%)合金的核/殼結(jié)構(gòu)粉末的微觀組織。該復(fù)合粉體具有以富(Fe-Cr)為核,以富銅為殼的復(fù)合結(jié)構(gòu),可以有效地避免鐵銅粉末燒結(jié)體和鐵、銅、石墨混合粉燒結(jié)體的組織不均勻性。由于富銅相具有較高的導(dǎo)電和導(dǎo)熱性能,而富Fe相具有較高的強度和硬度,因此,該類復(fù)合粉的燒結(jié)體可兼有高強高導(dǎo)性能。

    圖6 計算的Cu-Ni-Be三元體系中FCC(Cu)的固溶度曲線Fig.6 Calculated solubility line of FCC(Cu) in Cu-Ni-Be system

    圖7 計算的Cu-Ni-Sn體系中FCC相的溶解度間隙Fig.7 Calculated metasTablemiscibility gap of FCC phase in Cu-Ni-Sn system

    3.2.4 高性能銅基合金的相組成設(shè)計

    銅合金的相組成與合金的性能密切相關(guān),因此,通過計算銅基合金在溫度變化過程中相分數(shù)的變化以實現(xiàn)銅合金的組織設(shè)計。代表性實用合金Cu-9Ni-2.3Sn(質(zhì)量分數(shù),%)的相分數(shù)與溫度的關(guān)系的計算結(jié)果如圖10所示。由圖10可看出,隨著溫度的降低,Ni3Sn化合物相開始析出,在更低的溫度條件下,將出現(xiàn) FCC相的溶解度間隙。圖 11所示為Cu-3Ni-0.7Si(質(zhì)量分數(shù),%) 合金的相分數(shù)隨溫度變化的計算結(jié)果,其低溫組織由 FCC相和少量的 Ni5Si2,Ni2Si化合物組成。這些信息對析出強化型合金的組織設(shè)計有著重要的參考價值。

    3.2.5 Cu基板與焊料間溶解反應(yīng)的速度預(yù)測

    在半導(dǎo)體電子封裝中,基板與焊料之間發(fā)生的溶解反應(yīng),對于器件的可靠性來說是非常重要的。而該問題可以通過 DICTRA軟件[21],基于擴散動力學(xué)理論,預(yù)測其溶解的過程。例如,將Cu/Sn擴散偶在250℃下保溫,通過動力學(xué)計算,可以模擬銅在溶化的Sn中的溶解擴散行為。圖12(a)所示為銅在液相Sn中隨時間增加濃度變化的計算結(jié)果。從圖12(a)中可看出,熱處理時間為500 s時,液相濃度幾乎達到平衡位置。圖 12(b)所示為 FCC(Cu)/液相(Sn)的界面在不同時間的位置和移動速度。結(jié)果表明:在擴散的初始階段,銅快速溶解到液相(Sn)中,銅分解的速度隨時間的增加而下降,Cu/Sn界面快速向銅側(cè)移動。這些信息為電子封裝中焊接工藝的制定提供理論指導(dǎo)。

    圖8 計算的Cu-Fe-Cr三元系縱截面相圖Fig.8 Calculated vertical section diagram of Cu-Fe-Cr system

    圖9 Cu-32.4Fe-7.2Cr-0.4C合金粉末的微觀組織Fig.9 Microstructure of Cu-32.4Fe-7.2Cr-0.4C alloy powder

    圖10 Cu-9Ni-2.3Sn合金的相分數(shù)的計算結(jié)果Fig.10 Calculated results of phase mole fraction of Cu-9Ni-2.3Sn alloy

    圖11 Cu-3Ni-0.7Si合金相分數(shù)的計算結(jié)果Fig.11 Calculated results of phase mole fraction of Cu-3Ni-0.7Si system in Cu-Ni-Si system

    3.2.6 核/殼結(jié)構(gòu)無鉛焊接球的設(shè)計

    隨著電子產(chǎn)品和半導(dǎo)體器件向小型化和高密度化發(fā)展,電子行業(yè)對封裝技術(shù)的要求也越來越高。近年來,(Ball grid array, BAG)封裝方法已經(jīng)成為封裝互連技術(shù)的主流。與一般的封裝技術(shù)相比,BAG封裝技術(shù)具有更高的輸入/輸出終端密度、更小的臺面面積和更高的可靠性等優(yōu)點[22],這對于焊接材料本身提出了更高的要求。一般 BAG用的焊接球是通過霧化制粉的方法制備的 Sn基材料,雖然具有較好的焊接性能,但導(dǎo)電、導(dǎo)熱等性能還有待進一步的提高。而且目前無鉛焊接球大部分依賴國外進口,沒有我國自主知識產(chǎn)權(quán)。為此,研發(fā)新的復(fù)合粉體并制備成焊接性能優(yōu)良的焊接材料是電子封裝領(lǐng)域急需解決的問題。

    近年來,作者通過利用銅基合金的熱力學(xué)數(shù)據(jù)庫來設(shè)計和制備核/殼結(jié)構(gòu)低熔點無鉛焊接球[20]。圖 13所示為 Cu-15Sn-50Bi(質(zhì)量分數(shù),%)合金粉末的微觀組織。該粉體的組織具有以富Cu-Sn相為核心,以富Sn-Bi相為殼的結(jié)構(gòu)特點。目前工業(yè)上應(yīng)用的低熔點焊接球的尺寸通常為700 μm,但是工業(yè)上需要的尺寸在100 μm以下的微焊接球卻很難用傳統(tǒng)的制備方法獲得。本研究中(見圖13)的復(fù)合粉體的尺寸在80 μm左右,可以預(yù)見該類復(fù)合粉體在未來的微焊接領(lǐng)域具有良好的應(yīng)用前景。

    圖12 銅在錫中擴散的計算結(jié)果Fig.12 Dissolution behavior of Cu in molten Sn: (a) Concentration change of liquid phase; (b) Movement velocity and position of fcc and liquid boundary

    圖13 Cu-15Sn-50Bi合金核/殼復(fù)合粉體的組織Fig.13 Microstructure of core and shell of Cu-15Sn-50Bi alloy powder

    3.2.7 Sn-Bi-Cu復(fù)合粉焊膏

    通過調(diào)節(jié)霧化制粉工藝,在Sn-Bi-Cu合金系中也可以制備出高熔點的富(Cu, Sn)相彌散型顆粒分布在低熔點的富(Bi, Sn)相的基體上的彌散型復(fù)合粉體。該類復(fù)合粉體中的低熔點富(Bi, Sn)相,主要分布在復(fù)合粉的外部,熔點在138.1 ℃附近,可在較低的焊接溫度下熔化將元器件焊接起來。高熔點的富(Cu, Sn)相,主要分布在復(fù)合粉的內(nèi)部,熔點在587.1 ℃附近,在焊接溫度下不熔化,起導(dǎo)電和導(dǎo)熱作用。本研究組制備了Sn-Bi-Cu復(fù)合粉體(核/殼結(jié)構(gòu)的復(fù)合粉(見圖13)和彌散型復(fù)合粉的比例控制在一定范圍內(nèi))并研發(fā)出與其相匹配的助焊劑制成焊膏,焊膏在較低的焊接溫度下流動性良好。圖14所示為Sn-Bi-Cu復(fù)合粉焊膏的大焊盤焊點、小焊盤焊點和 BGA焊點的形貌。通過中試研究結(jié)果表明,Sn-Bi-Cu復(fù)合粉焊膏在 LED封裝、BGA封裝、Flip Chip封裝中,焊后效果良好,有望實現(xiàn)實用化。

    圖14 Sn-Bi-Cu復(fù)合粉焊膏的大焊盤焊點、小焊盤焊點和BGA焊點形貌圖Fig.14 Morphologies of big pad(a), small pad(b) and BGA welding spot(c)

    4 結(jié)論

    本研究利用CALPHAD方法,建立了Cu-X二元系以及 Cu-Fe、Cu-Ni、Cu-Cr、Cu-Co、Cu-Mo 和 Cu-W基各三元系相圖的熱力學(xué)數(shù)據(jù)庫。該數(shù)據(jù)庫是新型銅基合金成分和組織設(shè)計的重要工具,將為新型銅基合金的研發(fā)提供高效的理論指導(dǎo)。

    REFERENCES

    [1] SUNDMAN B, JANSSON B, ANDERSON J O. The Thermo-Calc databank system[J]. CALPHAD, 1985, 9(2):153-190.

    [2] SAUNDERS N, MIODOWNIK A P. CALPHAD (calculation of phase diagrams): A comprehensive guide[M]. Oxford, UK:Pergamon, 1998.

    [3] DINSDSLE A T. SGTE data for pure elemrnts[J]. CALPHAD,1991, 15(4): 317-425.

    [4] REDLICH O, KISTER A T. Algebraic representation of thermodynamic properties and the classification of solutions[J].Industr Eng Chem, 1948, (40): 345-348.

    [5] HILLERT M, JARL M. A model for alloying effects in ferromagnetic metals[J]. CALPHAD, 1978, 2: 227-238.

    [6] HILLERT M, STAFANSSON L I. The regular-solution model for stoichometric phases and ionic melts[J]. Acta Chem Scand,1970, 24: 3618-3626.

    [7] OHNUMA I, LIU X J, OHTANI H, ISHIDA K. Thermodynamic database for phase diagrams in micro-soldering alloys[J]. Journal of Electronic Materials Mater, 1999, 28(11): 1163-1171.

    [8] LIU X J, WANG C P, OHNUMA I, KAINUMA R, ISHIDA K.Thermodynamic assessment of the phase diagrams of the Cu-Sb and Sb-Zn systems[J]. Journal of Phase Equilibria, 2000, 21(5):432-442.

    [9] WANG C P, LIU X J, OHNUMA I, KAINUMA R, ISHIDA K.Thermodynamic database of the phase diagrams in Cu-Fe base ternary systems[J]. Journal of Phase Equilibria and Diffusion,2004, 25(4): 320-328.

    [10] WANG C P, LIU X J, JIANG M, OHNUMA I, KAINUMA R,ISHIDA K. Thermodynamic database of the phase diagrams in copper base alloy systems[J]. Journal of Physics and Chemistry of Solids, 2005, 66(2/4): 256-260.

    [11] WANG C P, LIU X J, OHNUMA I, KAINUMA R, ISHIDA K.Phase equilibria in Fe-Cu-X (X: Co, Cr, Si, V) ternary systems[J].Journal of Phase Equilibria, 2002, 23(3): 236-245.

    [12] LIU X J, WANG C P, OHNUMA I, KAINUMA R, ISHIDA K.Phase stability among the α(A1), β(A2), and γ(D83) phases in the Cu-Al-X system[J]. Journal of Phase Equilibria, 2001, 22(4):431-438.

    [13] WANG C P, LIU X J, OHNUMA I, KAINUMA R, ISHIDA K.Thermodynamic assessment of the Cu-Ni-Pb system[J].CALPHAD, 2000, 24(2): 149-167.

    [14] WANG C P, LIU X J, OHNUMA I, KAINUMA R, ISHIDA K,HAO S M. Phase equilibria in the Cu-Fe-Mo and Cu-Fe-Nb systems[J]. Journal of Phase Equilibria, 2000, 21(1): 54-62.

    [15] LIU X J, OHNUMA I, KAINUMA R, ISHIDA K. Phase equilibria in the Cu-rich portion of the Cu-Al binary system[J].Journal of Alloys and Compounds, 1998, 264(1/2):201-208.

    [16] LIU X J, WANG C P, OHNUMA I, KAINUMA R, ISHIDA K.Experimental investigation and thermodynamic calculation of the phase equilibria in the Cu-Sn and Cu-Sn-Mn systems[J]Metallurgical and Materials Transactions A, 2004,35(6):1641-1654.

    [17] WANG C P, LIU X J, OHNUMA I, HAO S M ,ISHIDA K.Ordering and phase separation of the bcc phase in the Fe-Cu-Al system[J]. Zeitschrift fuer Metallkunde, 1998, 89(12): 828-835.

    [18] 劉 平, 任風(fēng)章, 賈淑果. 銅合金及其應(yīng)用[M]. 北京: 化學(xué)出版社, 2007: 24-25.

    LIU Ping, REN Feng-zhang, JIA Shu-guo. Copper base alloys and its applications[M]. Beijing: Chemical Industry Press, 2007:24-25.

    [19] JAMIL S J, CHADWICK G A. Investigation and analysis of liquid phase sintering of Fe-Cu and Fe-Cu-C compacts[J].Powder Metall., 1985, 28(2): 65-71.

    [20] WANG C P, LIU X J, OHNUMA I, KAINUMA R, ISHIDA K.Formation of immiscible alloy powders with egg-type microstructure[J]. Science, 2002, 297(5583): 990-993.

    [21] ANDESSON J O, HOGLUND L, JONSSON B, AGREN J.Fundamental and application of ternary diffusion[M]. New York:Pergamon Press, 1990: 153-163.

    [22] BRADLEY E, BAERJI K. Effect of PCB finish on the reliability and wettability of ball grid array packages[J]. IEEE Trans Comp Packaging Manuf Technol B, 1996, 19(2): 320-330.

    Development of thermodynamic database for copper base alloy systems and its application in material design

    LIU Xing-jun1,2, WANG Cui-ping1,2, GAN Shi-xi1,2, I. Ohnuma3, R. Kainuma3, K. Ishida3
    (1. College of Materials, Xiamen University, Xiamen 361005, China;2. Institute of Advanced Nonferrous Metal Research, Xiamen University, Xiamen 361005, China;2. Department of Materials Science, Graduate School of Engineering, Tohoku University,Aoba-yama 6-6-02, Sendai 980-8579, Japan)

    A thermodynamic database on the phase equilibria of the copper base alloys, including Cu-X binary system and Cu-Fe, Cu-Ni, Cu-Cr, Cu-Co, Cu-Mo and Cu-W base ternary systems were developed by the calculation of phase diagrams (CALPHAD) method. The Gibbs free energies of the solution phases were described by subregular solution model with the Redlich-Kister equation, and the Gibbs free energies of the intermetallic compounds were described by sublattice model. The thermodynamic parameters describing the Gibbs energies of different phases were evaluated by optimizing the experimental data of phase equilibria and thermodynamic properties. The present thermodynamic database can provide much information, such as sTableand metasTablephase equilibria, phase fraction, liquidus projection and various thermodynamic quantities, which is expected to play an important role in the design of copper base alloys.

    Cu alloys; phase diagram; thermodynamic calculation; materials design

    TG113.14

    A

    1004-0609(2011)10-2511-12

    國家高新技術(shù)研究發(fā)展計劃資助項目(2009AA03Z101);國家科技部國際科技合作項目(2009DFA52170);國家自然科學(xué)基金資助項目(51031003)

    2011-05-10;

    2011-07-20

    劉興軍,教授,博士;電話:0592-2187888;E-mail: mylermm@163.com

    (編輯 李艷紅)

    猜你喜歡
    銅合金基合金熱力學(xué)
    超細鉬銅復(fù)合粉體及細晶鉬銅合金的制備
    珩磨技術(shù)在銅合金零件中的應(yīng)用
    銅及銅合金無硝酸光亮酸洗工藝應(yīng)用
    新型銅合金化學(xué)去毛刺劑的研制
    lncoloy 800鎳基合金加熱管腐蝕性能研究
    上海金屬(2016年3期)2016-11-23 05:19:43
    Fe-C-Mn-Si-Cr的馬氏體開始轉(zhuǎn)變點的熱力學(xué)計算
    上海金屬(2016年1期)2016-11-23 05:17:24
    活塞的靜力學(xué)與熱力學(xué)仿真分析
    電子制作(2016年19期)2016-08-24 07:49:54
    Alloy20鐵鎳基合金焊接
    焊接(2016年2期)2016-02-27 13:01:16
    高鉻鎳基合金熔敷金屬的組織和性能
    焊接(2015年9期)2015-07-18 11:03:53
    一類非奇異黑洞的熱力學(xué)穩(wěn)定性
    中文字幕最新亚洲高清| 日本午夜av视频| av片东京热男人的天堂| 中文字幕av电影在线播放| 人妻少妇偷人精品九色| 黑人巨大精品欧美一区二区蜜桃 | 欧美日韩精品成人综合77777| 欧美精品人与动牲交sv欧美| 最近中文字幕2019免费版| 免费观看性生交大片5| 日本与韩国留学比较| 波多野结衣一区麻豆| 国产精品成人在线| 日韩伦理黄色片| 亚洲av成人精品一二三区| 寂寞人妻少妇视频99o| 亚洲经典国产精华液单| 不卡视频在线观看欧美| 九草在线视频观看| 成年人免费黄色播放视频| 高清视频免费观看一区二区| 欧美+日韩+精品| 精品少妇黑人巨大在线播放| www.av在线官网国产| 亚洲av电影在线进入| 久久国产精品大桥未久av| 日本-黄色视频高清免费观看| 欧美xxⅹ黑人| 1024视频免费在线观看| 国产永久视频网站| 少妇 在线观看| 国产精品久久久久久av不卡| 在现免费观看毛片| 亚洲精品第二区| 欧美日韩亚洲高清精品| 97在线人人人人妻| 精品国产一区二区三区四区第35| 精品一区二区三卡| 另类亚洲欧美激情| 亚洲性久久影院| 国产一区二区三区综合在线观看 | 十分钟在线观看高清视频www| 在线观看三级黄色| 又粗又硬又长又爽又黄的视频| √禁漫天堂资源中文www| 亚洲国产精品一区三区| 久久99热这里只频精品6学生| 久久精品国产综合久久久 | 日本黄大片高清| 蜜臀久久99精品久久宅男| 亚洲内射少妇av| 国产精品熟女久久久久浪| 在线观看人妻少妇| 久久99一区二区三区| 一级毛片 在线播放| 午夜福利视频精品| 国产日韩欧美亚洲二区| 在线天堂最新版资源| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕av电影在线播放| 欧美成人精品欧美一级黄| 亚洲国产精品999| 亚洲一区二区三区欧美精品| 亚洲国产精品一区三区| 欧美丝袜亚洲另类| 国产毛片在线视频| 两性夫妻黄色片 | 亚洲第一av免费看| 日韩人妻精品一区2区三区| 一级片'在线观看视频| 最近手机中文字幕大全| 国产色婷婷99| 亚洲精品日韩在线中文字幕| 国产激情久久老熟女| 在线 av 中文字幕| 欧美变态另类bdsm刘玥| 美女国产视频在线观看| 日韩欧美精品免费久久| 亚洲欧美色中文字幕在线| 欧美日韩一区二区视频在线观看视频在线| 少妇的逼好多水| 国产精品久久久久久精品电影小说| 美女中出高潮动态图| 成年美女黄网站色视频大全免费| 国产男女内射视频| 国产黄色视频一区二区在线观看| 秋霞伦理黄片| 超色免费av| av天堂久久9| 秋霞在线观看毛片| 免费大片黄手机在线观看| 国产一区二区在线观看日韩| 91精品国产国语对白视频| 男人舔女人的私密视频| 香蕉丝袜av| 一本—道久久a久久精品蜜桃钙片| 咕卡用的链子| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 青春草亚洲视频在线观看| 久久99热6这里只有精品| 免费在线观看完整版高清| 国产精品国产三级专区第一集| 两个人看的免费小视频| 久久久久久久国产电影| 制服丝袜香蕉在线| 国产亚洲午夜精品一区二区久久| 一边摸一边做爽爽视频免费| 国产在线免费精品| 一本—道久久a久久精品蜜桃钙片| 99久久精品国产国产毛片| 午夜福利视频在线观看免费| 母亲3免费完整高清在线观看 | 在线观看美女被高潮喷水网站| 亚洲av综合色区一区| 99re6热这里在线精品视频| 午夜av观看不卡| 日日摸夜夜添夜夜爱| 18禁观看日本| 日本黄大片高清| 国产精品三级大全| 秋霞在线观看毛片| 久久精品人人爽人人爽视色| 高清av免费在线| 51国产日韩欧美| 亚洲精华国产精华液的使用体验| 国产极品粉嫩免费观看在线| 春色校园在线视频观看| 久久狼人影院| 欧美日韩成人在线一区二区| 精品少妇久久久久久888优播| 久热久热在线精品观看| 国产精品一区二区在线观看99| 日韩人妻精品一区2区三区| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av在线观看视频网站免费| 天堂8中文在线网| 黄色 视频免费看| 亚洲中文av在线| 日韩制服骚丝袜av| 日本欧美国产在线视频| 中国美白少妇内射xxxbb| 哪个播放器可以免费观看大片| 在线观看免费高清a一片| 高清在线视频一区二区三区| 日韩av在线免费看完整版不卡| 9色porny在线观看| 在线观看人妻少妇| 母亲3免费完整高清在线观看 | 久久韩国三级中文字幕| 乱码一卡2卡4卡精品| 在线看a的网站| 激情五月婷婷亚洲| 精品亚洲乱码少妇综合久久| 男女啪啪激烈高潮av片| 热re99久久国产66热| 国产成人精品一,二区| 亚洲精品日韩在线中文字幕| 精品国产国语对白av| 午夜影院在线不卡| 69精品国产乱码久久久| 欧美 亚洲 国产 日韩一| 赤兔流量卡办理| 欧美人与性动交α欧美精品济南到 | 男女免费视频国产| 男女高潮啪啪啪动态图| 久久精品国产亚洲av天美| 亚洲第一区二区三区不卡| 色哟哟·www| 另类精品久久| 亚洲第一区二区三区不卡| 91久久精品国产一区二区三区| 婷婷色av中文字幕| 最近中文字幕2019免费版| 有码 亚洲区| 美女视频免费永久观看网站| 欧美日韩一区二区视频在线观看视频在线| 国产综合精华液| av黄色大香蕉| 91国产中文字幕| videossex国产| 精品视频人人做人人爽| av福利片在线| 国产精品久久久久久精品电影小说| 色5月婷婷丁香| 色婷婷av一区二区三区视频| 天堂俺去俺来也www色官网| 成人免费观看视频高清| 又粗又硬又长又爽又黄的视频| 久久99热6这里只有精品| 久久女婷五月综合色啪小说| 午夜久久久在线观看| 亚洲av综合色区一区| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级| 丁香六月天网| 亚洲色图 男人天堂 中文字幕 | 亚洲精品美女久久av网站| 少妇人妻 视频| 久久久久网色| www.av在线官网国产| 亚洲综合色网址| 另类精品久久| 母亲3免费完整高清在线观看 | 最近手机中文字幕大全| 亚洲国产精品999| 99热网站在线观看| 欧美3d第一页| 一级片免费观看大全| 一区二区三区四区激情视频| 日韩一区二区视频免费看| 赤兔流量卡办理| 欧美亚洲 丝袜 人妻 在线| av卡一久久| 久久精品aⅴ一区二区三区四区 | 各种免费的搞黄视频| 在线免费观看不下载黄p国产| 日韩熟女老妇一区二区性免费视频| 在线观看三级黄色| 国产在线一区二区三区精| 久久久久久久久久久久大奶| 亚洲综合精品二区| 亚洲精品美女久久av网站| 欧美精品国产亚洲| 黑人巨大精品欧美一区二区蜜桃 | 女的被弄到高潮叫床怎么办| 日本黄大片高清| 校园人妻丝袜中文字幕| 久久毛片免费看一区二区三区| 日韩不卡一区二区三区视频在线| 久久精品夜色国产| 免费女性裸体啪啪无遮挡网站| 久久ye,这里只有精品| 亚洲图色成人| 午夜av观看不卡| www日本在线高清视频| 不卡视频在线观看欧美| 九色成人免费人妻av| 久久久久久人人人人人| 最新的欧美精品一区二区| 精品卡一卡二卡四卡免费| 少妇的逼水好多| 久久精品久久久久久噜噜老黄| 高清欧美精品videossex| 免费在线观看完整版高清| 亚洲成国产人片在线观看| 亚洲精品中文字幕在线视频| 国产女主播在线喷水免费视频网站| xxxhd国产人妻xxx| 少妇被粗大的猛进出69影院 | 国产成人欧美| 亚洲精品自拍成人| 国产精品嫩草影院av在线观看| 在线天堂中文资源库| 黑人猛操日本美女一级片| 国产乱人偷精品视频| 亚洲av国产av综合av卡| 中文乱码字字幕精品一区二区三区| 国产在线免费精品| 综合色丁香网| 欧美老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 亚洲国产精品国产精品| 国产欧美日韩综合在线一区二区| 麻豆精品久久久久久蜜桃| 亚洲欧美精品自产自拍| 99久久综合免费| 欧美激情 高清一区二区三区| 黑人高潮一二区| 黄色 视频免费看| 日韩av不卡免费在线播放| 99九九在线精品视频| 国产伦理片在线播放av一区| 精品久久国产蜜桃| 午夜影院在线不卡| av免费在线看不卡| 久久人人97超碰香蕉20202| 国产黄色免费在线视频| 国产精品一区www在线观看| 色94色欧美一区二区| 成人影院久久| 在线观看免费视频网站a站| 色婷婷av一区二区三区视频| 久久人人97超碰香蕉20202| 欧美日韩成人在线一区二区| a 毛片基地| 国产熟女欧美一区二区| 热re99久久国产66热| 深夜精品福利| 亚洲精品国产色婷婷电影| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 韩国精品一区二区三区 | 高清在线视频一区二区三区| 日本av免费视频播放| 男人爽女人下面视频在线观看| 色94色欧美一区二区| 久久人人爽av亚洲精品天堂| 热99久久久久精品小说推荐| 亚洲综合色网址| 久久精品国产鲁丝片午夜精品| 国产精品成人在线| 国产日韩欧美视频二区| 91午夜精品亚洲一区二区三区| 亚洲精品国产av成人精品| 五月开心婷婷网| 青春草亚洲视频在线观看| 免费人成在线观看视频色| 亚洲国产av新网站| 韩国精品一区二区三区 | 涩涩av久久男人的天堂| 亚洲伊人久久精品综合| 黑人巨大精品欧美一区二区蜜桃 | 久久久国产欧美日韩av| 久久久久久伊人网av| 纯流量卡能插随身wifi吗| 汤姆久久久久久久影院中文字幕| 亚洲五月色婷婷综合| 亚洲欧美一区二区三区黑人 | 亚洲美女搞黄在线观看| 久久精品久久久久久噜噜老黄| 激情视频va一区二区三区| 男女边摸边吃奶| 日本欧美国产在线视频| 寂寞人妻少妇视频99o| av国产久精品久网站免费入址| 99国产精品免费福利视频| 欧美少妇被猛烈插入视频| 国产精品麻豆人妻色哟哟久久| 日韩精品有码人妻一区| 亚洲国产成人一精品久久久| 日韩中字成人| 晚上一个人看的免费电影| 丝袜人妻中文字幕| 人妻系列 视频| 人成视频在线观看免费观看| a级毛片黄视频| 日本-黄色视频高清免费观看| 美女内射精品一级片tv| 国产男女超爽视频在线观看| 男男h啪啪无遮挡| 欧美日本中文国产一区发布| 看免费av毛片| 国产探花极品一区二区| av国产精品久久久久影院| 亚洲性久久影院| 久久久久人妻精品一区果冻| 欧美激情 高清一区二区三区| www.熟女人妻精品国产 | 在线观看美女被高潮喷水网站| 国产又色又爽无遮挡免| 国产精品一区二区在线不卡| kizo精华| 亚洲av电影在线观看一区二区三区| 少妇被粗大的猛进出69影院 | videos熟女内射| 日韩精品有码人妻一区| 99久久人妻综合| 日韩不卡一区二区三区视频在线| av线在线观看网站| 亚洲国产精品专区欧美| 涩涩av久久男人的天堂| 18禁裸乳无遮挡动漫免费视频| 极品少妇高潮喷水抽搐| 亚洲av男天堂| av黄色大香蕉| 日本av免费视频播放| 欧美日韩国产mv在线观看视频| 下体分泌物呈黄色| 女性被躁到高潮视频| 在线天堂中文资源库| 欧美另类一区| 成年av动漫网址| 亚洲精品aⅴ在线观看| 好男人视频免费观看在线| 色5月婷婷丁香| 国产国拍精品亚洲av在线观看| 日本欧美国产在线视频| 晚上一个人看的免费电影| 亚洲精品久久久久久婷婷小说| 久久99热这里只频精品6学生| 亚洲丝袜综合中文字幕| 最黄视频免费看| 日本欧美国产在线视频| 欧美变态另类bdsm刘玥| 日本wwww免费看| 亚洲国产欧美日韩在线播放| 国产福利在线免费观看视频| 18禁在线无遮挡免费观看视频| 国产探花极品一区二区| 少妇的逼好多水| 欧美 日韩 精品 国产| 免费高清在线观看视频在线观看| 亚洲国产精品成人久久小说| 波野结衣二区三区在线| 久久午夜综合久久蜜桃| 成年女人在线观看亚洲视频| www.色视频.com| 欧美成人午夜精品| 国产熟女欧美一区二区| 大片免费播放器 马上看| 国产免费又黄又爽又色| 五月伊人婷婷丁香| 亚洲av免费高清在线观看| 亚洲第一av免费看| 18禁裸乳无遮挡动漫免费视频| 久久精品国产a三级三级三级| 日韩伦理黄色片| 晚上一个人看的免费电影| 国产永久视频网站| 王馨瑶露胸无遮挡在线观看| 秋霞伦理黄片| 国产一区亚洲一区在线观看| 丰满饥渴人妻一区二区三| 你懂的网址亚洲精品在线观看| 国产成人av激情在线播放| h视频一区二区三区| 精品国产一区二区久久| 在线免费观看不下载黄p国产| 中文字幕另类日韩欧美亚洲嫩草| 黄色怎么调成土黄色| 精品一品国产午夜福利视频| 日本vs欧美在线观看视频| 99热这里只有是精品在线观看| 寂寞人妻少妇视频99o| 国产精品久久久久久av不卡| 欧美老熟妇乱子伦牲交| 欧美精品国产亚洲| 亚洲av在线观看美女高潮| 国产精品无大码| 精品一区二区三区视频在线| 国产精品一区二区在线观看99| 欧美精品av麻豆av| 国产xxxxx性猛交| 免费播放大片免费观看视频在线观看| 精品卡一卡二卡四卡免费| 久久久久久久亚洲中文字幕| 欧美bdsm另类| 韩国av在线不卡| 久久99一区二区三区| 欧美人与性动交α欧美软件 | 9191精品国产免费久久| 日本-黄色视频高清免费观看| av女优亚洲男人天堂| 视频在线观看一区二区三区| 香蕉精品网在线| 国产国语露脸激情在线看| 久久久久久人人人人人| 亚洲经典国产精华液单| 又黄又粗又硬又大视频| 人成视频在线观看免费观看| 国产精品一二三区在线看| 日本午夜av视频| 久久精品国产鲁丝片午夜精品| 色吧在线观看| 在线观看免费视频网站a站| 在线观看国产h片| 亚洲精品美女久久av网站| av国产久精品久网站免费入址| 制服丝袜香蕉在线| 欧美成人午夜精品| 视频区图区小说| 少妇精品久久久久久久| 亚洲欧洲精品一区二区精品久久久 | 欧美少妇被猛烈插入视频| 日韩成人伦理影院| 国产免费视频播放在线视频| 午夜福利在线观看免费完整高清在| 日韩 亚洲 欧美在线| 在线观看免费日韩欧美大片| 国产深夜福利视频在线观看| 国产一区亚洲一区在线观看| 日韩在线高清观看一区二区三区| 一二三四中文在线观看免费高清| 国产精品久久久久久久电影| 亚洲国产精品成人久久小说| 国产精品一国产av| 91成人精品电影| 国国产精品蜜臀av免费| 亚洲成av片中文字幕在线观看 | 2021少妇久久久久久久久久久| 精品人妻偷拍中文字幕| 女人被躁到高潮嗷嗷叫费观| 亚洲精品色激情综合| 久久这里有精品视频免费| 欧美日韩视频高清一区二区三区二| 成人免费观看视频高清| 成年人午夜在线观看视频| 免费高清在线观看视频在线观看| 高清毛片免费看| 亚洲av福利一区| 免费在线观看完整版高清| 国产精品久久久久久av不卡| 交换朋友夫妻互换小说| 久热久热在线精品观看| 少妇的丰满在线观看| 男女无遮挡免费网站观看| 免费大片18禁| 国产熟女欧美一区二区| 久久人人爽人人片av| 日韩精品有码人妻一区| 久久精品熟女亚洲av麻豆精品| 久久影院123| 国产成人精品在线电影| av女优亚洲男人天堂| 国产精品一国产av| 少妇精品久久久久久久| 久久热在线av| 欧美最新免费一区二区三区| 久久久久精品性色| 亚洲熟女精品中文字幕| 午夜福利视频在线观看免费| 久久影院123| 欧美最新免费一区二区三区| 国产黄色视频一区二区在线观看| 看非洲黑人一级黄片| 免费黄频网站在线观看国产| 亚洲精品国产色婷婷电影| 日韩熟女老妇一区二区性免费视频| 久久久久网色| 五月伊人婷婷丁香| 亚洲国产看品久久| 成人18禁高潮啪啪吃奶动态图| 国产免费现黄频在线看| 亚洲人与动物交配视频| 少妇人妻 视频| 久久这里有精品视频免费| 亚洲精品av麻豆狂野| 国产成人精品无人区| 热re99久久精品国产66热6| 亚洲av.av天堂| 国产午夜精品一二区理论片| 99久久综合免费| 久久免费观看电影| 丝瓜视频免费看黄片| 亚洲内射少妇av| 免费大片黄手机在线观看| 五月天丁香电影| 国产精品一国产av| 巨乳人妻的诱惑在线观看| 久久人人爽av亚洲精品天堂| av.在线天堂| 妹子高潮喷水视频| 五月玫瑰六月丁香| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 最后的刺客免费高清国语| 美女大奶头黄色视频| 国产爽快片一区二区三区| 精品久久国产蜜桃| 久久精品人人爽人人爽视色| 97精品久久久久久久久久精品| 一区二区av电影网| 赤兔流量卡办理| 9191精品国产免费久久| 超色免费av| 午夜福利视频在线观看免费| 综合色丁香网| 久久99一区二区三区| 国产欧美日韩综合在线一区二区| 一区二区日韩欧美中文字幕 | av视频免费观看在线观看| av卡一久久| 国产黄色免费在线视频| 国产亚洲精品第一综合不卡 | 国产男女超爽视频在线观看| 国产日韩欧美视频二区| 亚洲一级一片aⅴ在线观看| 欧美人与善性xxx| 欧美人与性动交α欧美软件 | 免费黄频网站在线观看国产| 视频中文字幕在线观看| 在线看a的网站| 狠狠精品人妻久久久久久综合| 日本猛色少妇xxxxx猛交久久| videosex国产| 在线免费观看不下载黄p国产| 美女国产视频在线观看| 国产毛片在线视频| 色婷婷久久久亚洲欧美| 99热国产这里只有精品6| 精品一品国产午夜福利视频| 美女内射精品一级片tv| 大香蕉久久成人网| 国产高清三级在线| xxxhd国产人妻xxx| 22中文网久久字幕| 亚洲av在线观看美女高潮| 国产日韩一区二区三区精品不卡| 日本-黄色视频高清免费观看| 欧美 日韩 精品 国产| 夫妻午夜视频| 欧美日韩综合久久久久久| 激情视频va一区二区三区| 在线观看三级黄色| 在线精品无人区一区二区三| 国内精品宾馆在线| av免费观看日本| 人妻 亚洲 视频| 99九九在线精品视频| 制服诱惑二区| h视频一区二区三区| 亚洲国产最新在线播放| 国产 一区精品| 在线观看免费视频网站a站| 9191精品国产免费久久| 国产一区二区三区av在线| 在线观看免费视频网站a站| 99九九在线精品视频| 国产 一区精品| 久热这里只有精品99| 国产在视频线精品| 看非洲黑人一级黄片| 97超碰精品成人国产|