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

    Al-Ti-B合金中AlB2、TiB2和TiAl3的第一性原理研究

    2018-10-11 06:07:20黃元春邵虹榜肖政兵任賢魏郭曉芳
    中國有色金屬學報 2018年8期
    關鍵詞:熱容費米結合能

    黃元春,邵虹榜,肖政兵,任賢魏,郭曉芳

    ?

    Al-Ti-B合金中AlB2、TiB2和TiAl3的第一性原理研究

    黃元春1, 2, 3,邵虹榜1,肖政兵1, 2, 3,任賢魏1,郭曉芳1

    (1. 中南大學 輕合金研究院, 長沙 410083;2. 中南大學 高性能復雜制造國家重點實驗室, 長沙 410083;3. 中南大學 機電工程學院, 長沙 410083)

    采用基于密度泛函理論的第一性原理方法,研究0 K條件下Al-Ti-B體系中AlB2、TiB2和TiAl3的平衡晶格常數(shù)、形成熱、結合能及電子特性,并結合準簡諧德拜模型對3種物相在高溫條件下的熱力學基本性質(zhì)進行了分析。計算結果表明:AlB2、TiB2和TiAl3的形成熱與結合能均為負值,表明3種相均具有結構穩(wěn)定性,且穩(wěn)定性由大到小依次為TiB2、TiAl3、AlB2;3種物相的離子性與金屬性由大到小依次為TiAl3、AlB2、TiB2,TiAl3與AlB2顯示離子鍵特性,而TiB2鍵合特性以共價鍵為主;3種物相的熱穩(wěn)定性由大到小依次為TiB2、AlB2、TiAl3。計算所得TiB2、AlB2和TiAl3的熱力學性質(zhì)可以用于分析氟鹽法制備Al-Ti-B時Al熔體與氟鹽之間反應發(fā)生情況以及中間合金最終的物相組成。

    Al-Ti-B合金;第一性原理;準簡諧德拜模型;結構特性;熱力學性質(zhì)

    商業(yè)用Al-Ti-B的常用制備工藝是將高純鋁錠熔化,通過加入不同比例的混合氟鹽(K2TiF6、KBF4)至鋁熔體中并保溫反應一段時間后澆鑄,從而得到不同Ti/B比例的Al-Ti-B中間合金[1?2]。Al-Ti-B中的第二相粒子(如TiB2、TiAl3和(Al,Ti)B2等)有助于鋁發(fā)生異質(zhì)形核而達到細化晶粒的效果。例如,Al-5Ti-1B加入至Al熔體后釋放出的高熔點TiB2粒子可以作為(Al)異質(zhì)形核的核心,低熔點的TiAl3在鋁熔體中能被迅速溶解并釋放出Ti元素,為(Al)發(fā)生異質(zhì)形核創(chuàng)造了便利條件[3?4];而KORI等[5]研究發(fā)現(xiàn)Al-1Ti-3B中AlB2含量較高,適于Al-Si合金的晶粒細化。盡管Al-Ti-B中間合金具有細化效率高、成本低廉等優(yōu)點而廣受青睞,然而在實際應用過程中仍存在如內(nèi)部第二相粒子聚集沉淀、尺寸寬化等諸多缺陷[6?7],嚴重限制Al-Ti-B在高質(zhì)量鋁材生產(chǎn)中的應用。因此,深入了解這些異質(zhì)形核粒子的結構與性質(zhì)就顯得尤為重要。

    到目前為止,不少研究者[8?10]試圖解釋制備過程中氟鹽在鋁熔體中的反應機制及情況,以改良Al-Ti-B的制備工藝,從而進一步開發(fā)Al-Ti-B的細化潛能。然而,K2TiF6、KBF4與Al的反應生成TiB2、TiAl3等中間相的過程受化學反應機制、元素擴散等眾多因素影響,難以精確分析熔體中氟鹽反應機理。本文作者采用基于密度泛函理論(Density functional theory,DFT)[11?12]的第一性原理方法,研究AlB2、TiB2和TiAl3的結構穩(wěn)定性以及電子特性,并結合準簡諧德拜模型對這3種物相在高溫條件下的熱力學基本性質(zhì)進行分析,期望從原子角度揭示3種物相在Al-Ti-B體系下的基本性質(zhì),為開發(fā)新型高品質(zhì)Al-Ti-B制備工藝提供理論支撐。

    1 計算方法與模型

    在結構穩(wěn)定性與電子結構的研究方面,本文作者采用基于密度泛函理論的第一性原理平面波贗勢方法,并在計算過程中使用廣義梯度近似 (Generalized gradient approximation, GGA)[13]處理交換關聯(lián)能,交換關聯(lián)勢采取Perdew-Burke-Ernzerhof (PBE)[14]形式 進行。各元素價電子的選取情況如下:B-2s22p1、Al-3s23p1、Ti-3s23p63d24s2。AlB2、TiB2和TiAl3晶胞模型簡約布里淵區(qū)K-point采取Monkhorst-Pack[15]方法劃分,單胞網(wǎng)格數(shù)分別取10×10×8、10×10×8和8×8×8。自洽循環(huán)計算的能量收斂精度設置為1×10?6eV/atom,各個原子之間的相互作用力低于0.01 eV/nm,公差偏移低于5.0×10?4nm,應力偏差為低于0.05 GPa。

    采用準簡諧德拜模型[16]對TiB2、AlB2以及TiAl3的等容摩爾熱容、熵值以及熱平衡體積與熱膨脹系數(shù)隨溫度的變化進行了研究。3種物相的非平衡Gibbs函數(shù)*(;,)為可描述為

    式中:()表示3種物相每個原胞體積下的總能量;、和分別為壓強、體積和溫度;()為德拜溫度;Vib為振動Helmholtz自由能。

    而基于準諧近似和使用聲子態(tài)密度的德拜模型,Vib可以用以下形式表示[17]:

    式中:代表每個原胞中所含原子數(shù)目;為波爾茲曼常數(shù);(/)為德拜積分。

    各向同性固體的德拜溫度表示如下[16]:

    式中:為每個原胞中分子的質(zhì)量;s是表征晶體壓縮率的絕熱體彈模量,可以近似寫為

    2 結果與討論

    2.1 晶體結構與穩(wěn)定性

    TiB2(AlB2)屬于六方晶系,空間點群為6/,每個晶胞含有1個Ti(Al)原子和2個B原子,Ti(Al)原子占據(jù)六棱柱的頂角和底心位置,而硼原子處于中心;TiAl3具有四方晶體結構,屬于4/空間點群,晶胞由兩個面心立方晶格堆垛構成,包含2個Ti原子和6個Al原子。通過充分馳豫晶胞結構以及晶內(nèi)原子位置以達到能量最低狀態(tài),對3種物相的結構穩(wěn)定性進行研究分析。

    圖1 TiB2、AlB2和TiAl3的晶體結構模型

    晶體形成熱(Heat of formation, Δ)用以表征物質(zhì)由單質(zhì)狀態(tài)轉變成晶體化合物狀態(tài)的反應過程中能量的變化情況:Δ為負值時,表示化合物在形成過程中放熱,反應可自發(fā)進行;反之則為吸熱,反應需要外界提供能量才能進行,且生成的化合物不穩(wěn)定。結合能(Cohesive energy,coh)是指原子由自由態(tài)結合成為晶體化合物的過程中所釋放的能量,用來表征化合物形成之后在結構上的穩(wěn)定性,結合能的絕對值越大,則形成的晶體化合物越穩(wěn)定[21]。對于化合物AB,其形成熱與結合能的計算公式分別如下[22]:

    AlB2、TiB2和TiAl3的晶體結構參數(shù)以及形成熱與結合能的計算結果如表1所示。通過比較可以發(fā)現(xiàn),所得計算結果與實驗測量值誤差較小,表明本研究所采用的計算參數(shù)設置合理。

    計算結果顯示:TiB2、AlB2以及TiAl3的形成熱均為負值,表明3種物相的形成過程無需外界提供熱量,均能自發(fā)進行;從結合能的計算結果來看,3種物相的結合能絕對值由大到小的依次為TiB2、TiAl3、AlB2,說明Ti原子與B原子結合能力最強,所形成的TiB2相熔點最高。此外,比較兩種硼化物可以發(fā)現(xiàn)TiB2的結合能低于AlB2,說明在Al-Ti-B中間合金中TiB2更加穩(wěn)定,與文獻[23]的研究結果相符。形成熱與結合能的結果表明,3種金屬間化合物的穩(wěn)定性由大到小依次為TiB2、TiAl3、AlB2。

    2.2 電子結構

    為了解TiB2、AlB2和TiAl3的成鍵特性,并進一步揭示出3種物相的結構穩(wěn)定性差別,本研究中分別計算了3種物相的總態(tài)密度(Total density of state, DOS)和分態(tài)密度(Partial density of state, PDOS),結果如圖2所示,圖中虛線表示費米能級(Fermi level)。

    表1 計算所得TiB2、AlB2以及TiAl3的晶體結構參數(shù)、形成熱與結合能以及實驗值

    圖2 物相的分態(tài)密度與總態(tài)密度

    TiB2、AlB2和TiAl3在費米能級處的態(tài)密度均大于0,顯示不同程度的金屬特性。圖2(a)顯示TiB2中B 2s與2p的態(tài)密度發(fā)生重疊而產(chǎn)生sp態(tài)雜化。Ti原子PDOS顯示3d軌道在費米能級附近呈現(xiàn)出兩個尖峰,而與3d軌道重合的s和p軌道在費米能級附近的態(tài)密度較弱,這是由于Ti 3d、4s和4p軌道之間的雜化而造成的[28]。B 2p與Ti 3d軌道態(tài)密度尖峰分列在費米能級左右而顯示出較強的雜化,并使得TiB2的總態(tài)密度在費米能級處出現(xiàn)贗能隙,反映了TiB2存在較強共價鍵特性;圖2(b)顯示AlB2中B 2s與2p態(tài)密度同樣存在重疊現(xiàn)象,但與TiB2相比,重疊的能量范圍更為寬泛,而且在費米能級以上(0~20 eV)空軌道中的態(tài)密度更大,這表明AlB2的反鍵合態(tài)更高。另一方面,Al原子在費米能級以下存在的sp態(tài)雜化較弱,而且在高能量區(qū)域的反鍵合態(tài)較強。因此,Al原子的價電子對AlB2的共價鍵幾乎沒有貢獻。AlB2的總態(tài)密度相對平滑,而且在費米能級處存在平坦的能隙,意味著成鍵軌道與反鍵軌道之間的能量差異較小,共價鍵的鍵合強度比TiB2弱。圖2(c)顯示TiAl3的總態(tài)密度集中在?10~15 eV的能量范圍,并且在費米能級附近存在贗能隙,表明TiAl3同樣具有一定的共價鍵特性。然而該能隙范圍較窄,說明TiAl3的共價鍵強度弱于兩種硼化物。分析TiAl3的PDOS可知其總態(tài)密主要來自Al 3s和3p以及Ti 3d軌道價電子貢獻(d軌道的PDOS為較大的尖峰,該軌道電子相對比較局域)。此外,Al 3p與Ti 3d的PDOS在費米能級附近發(fā)生重疊,存在Al 3s和Ti 3d態(tài)電子雜化。

    Mulliken電荷布居分析可以定量地描述3種物相內(nèi)部電荷轉移情況,其結果如表2所示。

    從表2中可以看出,TiB2中Ti原子轉移到B原子的電荷數(shù)為1.16,轉移的電荷主要來自Ti原子4s軌道;AlB2中Al原子轉移到B原子的電荷數(shù)為1.29;而TiAl3中除了Al 3s和Ti 3d軌道電子雜化產(chǎn)生電荷轉移(Ti原子轉移到Al原子的電荷數(shù)為0.06)外,Al的sp軌道雜化使得單位Al原子之間產(chǎn)生0.10的電荷轉移。以上結果表明,3種相均出現(xiàn)了電荷轉移,原子之間除了共價鍵之外,還存在著離子鍵特性,離子性按以下順序增加,由大到小依次為TiAl3、AlB2、TiB2。

    表2 TiB2、AlB2和TiAl3 的Mulliken電荷布居

    TiB2、AlB2以及TiAl3的金屬性(m)可以通過以下計算公式[29]得到

    式中:F為費米能級處態(tài)密度值;m、e分別為晶胞內(nèi)熱激活電子、價電子的密度,而e可由晶胞體積(cell)與總價電子數(shù)()計算得到:e=/cell。相關參數(shù)以及計算結果如表3所示。由表3可以看到,TiAl3的m相對較大,其次是AlB2,而TiB2的金屬性最小。

    表3 物相費米能級處態(tài)密度值、總價電子數(shù)、晶胞體積以及金屬性參數(shù)

    2.3 熱力學性質(zhì)

    圖3所示為零壓條件下計算所得TiB2、AlB2和TiAl3的定容摩爾熱容c與熵在0~1200 K隨溫度的變化關系。圖3(a)顯示在室溫(300 K)條件下3種物相的c分別為48.715、47.878和84.394 J/(mol?K),隨著溫度升高,3種物相的熱容均有不同程度的增大。0~300 K的低溫條件下,熱容主要受到晶內(nèi)原子熱振動的影響,故在此階段3種相的c隨著溫度的升高而迅速增大,服從德拜3定律;當溫度升至300 K以上時,熱容增速開始減小,但TiAl3的熱容增速仍然高于兩種硼化物;而與絕大多數(shù)固體在高溫下的情況相同,在>1000 K的高溫區(qū)域內(nèi),3種物相內(nèi)部非簡諧效應被抑制,c增速逐漸放緩并趨近于各自的Dulong-petit極限。在0~1200 K溫度范圍內(nèi),兩種硼化物粒子熱容變化情況幾乎一致,而這可能是由于兩種物相晶型相同,內(nèi)部原子的振動細節(jié)相似而造成的。進一步考察高溫區(qū)域(900~1200 K),TiB2和AlB2的熱容與溫度大致呈線性關系,且前者熱容始終大于后者。

    表4 TiB2、AlB2和TiAl3的體彈模量B0以及其對壓強的一階導數(shù)B'

    1) Cal: Calculated; Ref: Reference; Exp: Experimental.

    圖3 TiB2、AlB2與TiAl3的熱容cV與熵S的關系

    圖3(b)顯示出零壓條件下TiB2、AlB2與TiAl3的熵隨溫度的變化情況??梢钥闯鲈?00~500 K溫度范圍內(nèi),TiAl3的熵值隨溫度的升高而劇烈增大。當溫度繼續(xù)上升時,熵增速率放緩(可以預計當溫度高于某個值后,熵會逐漸趨于一個極限值);而TiB2和AlB2的熵值增大速率低于TiAl3,且當>200 K后,熵增速率比較恒定。從900~1200 K高溫區(qū)內(nèi)變化情況來看,TiB2的熵值大于AlB2,且兩種硼化物的熵值與溫度基本呈線性增加的關系。900 K時,TiB2、AlB2和TiAl3的熵分別為102.367、101.191和187.088 J/(mol?K)。

    圖4 TiB2、AlB2和TiAl3的熱平衡體積V與熱膨脹系數(shù)隨溫度的變化

    采用氟鹽法制備Al-Ti-B時,鋁熔體與混合氟鹽(K2TiF6和KBF4)之間發(fā)生的反應以及最終得到的物相組成可以從以上計算結果加以分析:低溫熔煉時,具有較高熱容與熵值的TiAl3需要一定的過熱度才能形成,而K2TiF6加入鋁熔體后立即與鋁發(fā)生鋁熱反應并大量放熱,使得熔體溫度迅速升高,有利于大量生成TiAl3[34]。隨著反應進行,熔體溫度逐漸升高,呈金屬鍵與離子鍵混合鍵合特性的TiAl3熱穩(wěn)定性較差而易熔解于鋁液中,并使鋁液中Ti的濃度升高,熔體中將存有大量的TiAl3塊狀小晶體;高溫熔煉時,TiAl3熔解速率加快,少量未徹底溶解的TiAl3受到鋁熔體中Ti元素的長程擴散的影響而擇優(yōu)生長,從而長大成為桿狀或針片狀。

    TiB2的熱容與熵值均高于AlB2,生成TiB2的自由焓更低,氟鹽與鋁液反應首先生成TiB2,剩余的KBF4與Al熔體才會生成AlB2,同時,高溫條件下AlB2的熱穩(wěn)定性差于TiB2而處于亞穩(wěn)定狀態(tài),由于Ti的存在,Ti將置換AlB2中的Al原子而逐漸形成(Al, Ti1?x)B2復合相直至TiB2[9]。因此,在高(Ti)/(B)的Al-Ti-B中間合金中AlB2含量遠低于TiB2含量;在制備低(Ti)/(B)的Al-Ti-B(如Al-1Ti-3B)時,由于Ti含量相對較少,且Al熔體中添加過多的KBF4后黏度增大,在局部區(qū)域內(nèi)Ti分布不均勻,因此,AlB2不能完全轉化TiB2而存在部分AlB2在凝固過程中保留下來。

    3 結論

    1) 形成熱與結合能的計算結果表明:3種金屬間化合物的穩(wěn)定性由大到小依次為TiB2、AlB2、TiAl3。TiB2的結合能低于AlB2,表明在Al-Ti-B中間合金中TiB2更加穩(wěn)定;

    2) 電子結構分析表明:3種物相的離子性與金屬性由大到小依次為TiAl3、AlB2、TiB2。TiAl3與AlB2的鍵合特性以離子鍵為主,而TiB2鍵合特性以B 2p與Ti 3d軌道態(tài)電子雜化的共價鍵為主;

    3) 準簡諧德拜模型計算顯示:TiB2、AlB2和TiAl3的定容摩爾熱容與熵隨溫度升高而不同程度地增大,通過計算熱膨脹系數(shù)發(fā)現(xiàn)3種物相的熱穩(wěn)定性由大到小依次為TiB2、AlB2、TiAl3。

    [1] AURADI V, KORI S A. Influence of reaction temperature and reaction time for the manufacturing of Al-Ti-B (Ti:B=5:1, 1:3) master alloys and their grain refining efficiency on Al-7Si alloys[J]. Transactions of the Indian Institute of Metals, 2012, 65(6): 637?645.

    [2] BIROL Y. An improved practice to manufacture Al-Ti-B master alloys by reacting halide salts with molten aluminium[J]. Journal of Alloys and Compounds, 2006, 420(1/2): 71?76.

    [3] JONES G P, PEARSON J. Factors affecting the grain-refinement of aluminum using titanium and boron additives[J]. Metallurgical and Materials Transactions B, 1976, 7(2): 223?234.

    [4] 李 賀, 柴麗華, 馬騰飛, 陳子勇. 高溫熔體反應法制備Al-5Ti-1B細化劑[J]. 材料工程, 2017, 45(2): 39?45. LI He, CHAI Li-hua, MA Teng-fei, CHEN Zi-yong. Synthesis of Al-5Ti-1B refiner by melt reaction method[J]. Journal of Materials Engineering, 2017, 45(2): 39?45.

    [5] KORI S A, MURTY B S, CHAKRABORTY M. Development of an efficient grain refiner for Al-7Si alloy and its modification with strontium[J]. Materials Science and Engineering A, 2000, 280(1): 58?61.

    [6] YU L, LIU X. The relationship between viscosity and refinement efficiency of pure aluminum by Al-Ti-B refiner[J]. Journal of Alloys and Compounds, 2006, 425(1/2): 245?250.

    [7] 肖政兵, 鄧運來, 唐建國, 陳 祺, 張新明. Al-Ti-C與Al-Ti-B晶粒細化劑的Zr中毒機理[J]. 中國有色金屬學報, 2012, 22(2): 371?378. XIAO Zheng-bing, DENG Yun-lai, TANG Jian-guo, CHEN Qi, ZHANG Xin-ming. Poisoning mechanism of Zr on grain refiner of Al-Ti-C and Al-Ti-B[J]. The Chinese Journal of Nonferrous Metals, 2012, 22(2): 371?378.

    [8] 廖成偉, 陳聞天, 陳 歡, 付建平, 潘春旭. 原料添加順序對氟鹽法制備Al-Ti-B中間合金的影響[J]. 中國有色金屬學報, 2016, 26(1): 204?211.LIAO Cheng-wei, CHEN Wen-tian, CHEN Huan, FU Chun-ping, PAN Chun-xu. Effect of feeding order on preparation of Al-Ti-B master alloy by fluoride salt method[J]. The Chinese Journal of Nonferrous Metals, 2016, 26(1): 204?211.

    [9] FJELLSTEDT J, JARFORS A. On the precipitation of TiB2in aluminum melts from the reaction with KBF4and K2TiF6[J]. Materials Science and Engineering A, 2005, 413(6): 527?532.

    [10] MIRKOVIC D, GROEBNER J, SCHMID-FETZER R, FABRICHNAYA O, LUKAS H L. Experimental study and thermodynamic re-assessment of the Al-B system[J]. Journal of Alloys and Compounds, 2004, 384(1/2): 168?174.

    [11] HERRING C. A new method for calculating wave functions in crystals[J]. Physical Review, 1940, 57(12): 1169?1177.

    [12] KOHN W, SHAM L J. Self-consistent equations including exchange and correlation effects[J]. Physical Review, 1965, 140(4A): A1133/A1138.

    [13] XIE Y P, WANG Z Y, HOU Z F. The phase stability and elastic properties of MgZn2and Mg4Zn7in Mg-Zn alloys[J]. Scripta Materialia, 2013, 68(7): 495?498.

    [14] PERDEW J P, BURKE K, ERNZERHOF M. ERRATA: Generalized gradient approximation made simple[J]. Physical Review Letters, 1996, 77(18): 3865?3868.

    [15] MONKHORST H J. Special points for Brillouin-zone integrations[J]. Physical Review B Condensed Matter, 1976, 16(4): 1748?1749.

    [16] BLANCO M A, FRANCISCO E, LUA?A V. GIBBS: Isothermal-isobaric thermodynamics of solids from energy curves using a quasi-harmonic Debye model[J]. Computer Physics Communications, 2004, 158(1): 57?72.

    [17] FLóREZ M, RECIO J M, FRANCISCO E, BLANCO M A, PENDáS A M. First-principles study of the rocksalt-cesium chloride relative phase stability in alkali halides[J]. Physical Review B, 2002, 66(14): 144112?144123.

    [18] WANG M. The investigation of dependences of mechanical and electronic properties of TiB2on pressure using the first-principles method[J]. Physica Scripta, 2014, 89(11): 115702?115714.

    [19] DUAN Y H, SUN Y, GUO Z Z, PENG M J, ZHU P X, HE J H. Elastic constants of AlB2-type compounds from first-principles calculations[J]. Computational Materials Science, 2012, 51(1): 112?116.

    [20] DUAN Y H, SUN Y, LU L. Thermodynamic properties and thermal conductivities of TiAl3-type intermetallics in Al-Pt-Ti system[J]. Computational Materials Science, 2013, 68(2): 229?233.

    [21] MEDVEDEVA N I, GORNOSTYREV Y N, NOVIKOV D L, MRYASOV O N, FREEMAN A J. Ternary site preference energies, size misfits and solid solution hardening in NiAl and FeAl[J]. Acta Materialia, 1998, 46(10): 3433?3442.

    [22] WU M M, JIANG Y, WANG J W, WU J, TANG B Y, PENG L M, DING W J. Structural, elastic and electronic properties of Mg(Cu1?xZnx)2alloys calculated by first-principles[J]. Journal of Alloys and Compounds, 2011, 509(6): 2885?2890.

    [23] FJELLSTEDT J, JARFORS A E W, SVENDSEN L. Experimental analysis of the intermediary phases AlB2, AlB12and TiB2in the Al-B and Al-Ti-B systems[J]. Journal of Alloys & Compounds, 1999, 283(1/2): 192?197.

    [24] KELLY P M, REN H P, QIU D, ZHANG M X. Identifying close-packed planes in complex crystal structures[J]. Acta Materialia, 2010, 58(8): 3091?3095.

    [25] VAJEESTON P, RAVINDRAN P, RAVI C, ASOKAMANI R. Electronic structure, bonding, and ground-state properties of AlB2-type transition-metal diborides[J]. Physical Review B, 2001, 63(4): 33?41.

    [26] BURKHARDT U, GURIN V, HAARMANN F, BORRMANN H, SCHNELLE W, YARESKO A, GRIN Y. On the electronic and structural properties of aluminum diboride Al0.9B2[J]. Journal of Solid State Chemistry, 2004, 177(2): 389?394.

    [27] ZHANG M X, KELLY P M, EASTON M A, TAYLOR J A. Crystallographic study of grain refinement in aluminum alloys using the edge-to-edge matching model[J]. Acta Materialia, 2005, 53(5): 1427?1438.

    [28] ZHANG H L, HAN Y F, WANG J, DAI Y B, SUN B D. An Ab initio molecular dynamics study on the structural and electronic properties of AlB2, TiB2and (Al,Ti(1?x))B2in Al-Ti-B master alloys[J]. Journal of Alloys and Compounds, 2014, 585(3): 529?534.

    [29] LI Y, GAO Y, XIAO B, MIN T, FAN Z, MA S, XU L. Theoretical study on the stability, elasticity, hardness and electronic structures of W-C binary compounds[J]. Journal of Alloys and Compounds, 2010, 502(1): 28?37.

    [30] SPOOR P S, MAYNARD J D, PAN M J, GREEN D J, HELLMANN J R, TANAKA T. Elastic constants and crystal anisotropy of titanium diboride[J]. Applied Physics Letters, 1997, 70(15): 1959?1961.

    [31] LOA I, KUNC K, SYASSEN K, BOUVIER P. Crystal structure and lattice dynamics of AlB2under pressure and implications for MgB2[J]. Physical Review B, 2002, 66(13): 134101.

    [32] GHOSH G, WALLE A V D, ASTA M. First-principles phase stability calculations of pseudobinary alloys of (Al,Zn)3Ti with L12, D022, and D023structures[J]. Journal of Phase Equilibria and Diffusion, 2007, 28(1): 9?22.

    [33] OKADA Y, TOKUMARU Y. Precise determination of lattice parameter and thermal expansion coefficient of silicon between 300 and 1500 K[J]. Journal of Applied Physics, 1984, 56(2): 314?320.

    [34] EL-MAHALLAWY N, TAHA M A, JARFORS A E W, FREDRIKSSON H. On the reaction between aluminium, K2TiF6and KBF4[J]. Journal of Alloys and Compounds, 1999, 292(1/2): 221?229.

    First principle study of AlB2, TiB2and TiAl3in Al-Ti-B alloy

    HUANG Yuan-chun1, 2, 3, SHAO Hong-bang1, XIAO Zheng-bing1, 2, 3, REN Xian-wei1, GUO Xiao-fang1

    (1. Light Alloy Research Institute, Central South University, Changsha 410083, China; 2. State Key Laboratory of High Performance and Complex Manufacturing, Central South University, Changsha 410083, China; 3. College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China)

    The crystal structure, heat of formation, cohesive energy and electronic properties of AlB2, TiB2and TiAl3in Al-Ti-B system under 0 K condition were investigated by using first principles based on density functional theory, and the quasi-harmonic Debye model was performed to calculate the thermodynamic properties of these phases. The results show that the heat of formation and cohesive energy of AlB2, TiB2and TiAl3are both negative, which indicate that these phases possess structural stability, and the structural stability decreases in the order from big to little of TiB2, TiAl3, AlB2. The iconicity and metallicity of these compounds increase in the following sequence from big to little of TiAl3, AlB2, TiB2, and the dominant bonding in TiAl3and AlB2is ionic, while the bonding in TiB2is mainly covalent. The thermodynamic calculation results show that the thermal stabilities of the three phases decrease in the order from big to little of TiB2, AlB2, TiAl3, which can be used to illustrate the reactions between Al melt and fluoride salts when preparing Al-Ti-B by fluoride salt process.

    Al-Ti-B alloy; first-principles; quasi-harmonic Debye model; structural properties; thermodynamics properties

    Project(2015CB057305) supported by the National Basic Research Development Program of China

    2017-06-05;

    2017-09-20

    HUANG Yuan-chun; Tel: +86-13507315123; E-mail: science@csu.edu.cn

    10.19476/j.ysxb.1004.0609.2018.08.01

    1004-0609(2018)-08-1491-08

    TG146.2

    A

    國家重點基礎研究發(fā)展計劃資助項目(2015CB057305)

    2017-06-05;

    2017-09-20

    黃元春,教授,博士;電話:13507315123;E-mail: science@csu.edu.cn

    (編輯 李艷紅)

    猜你喜歡
    熱容費米結合能
    費米能區(qū)重離子反應中對稱能系數(shù)的提取
    晶體結合能對晶格動力學性質(zhì)的影響
    大學物理(2022年9期)2022-09-28 01:10:52
    借鑒躍遷能級圖示助力比結合能理解*
    物理通報(2020年7期)2020-07-01 09:28:02
    費米氣體光晶格模型的漸近軌線
    費米:假裝是司機
    熱容式熱流密度傳感器的研制及應用
    中國測試(2018年10期)2018-11-17 01:58:50
    ε-CL-20/F2311 PBXs力學性能和結合能的分子動力學模擬
    對“結合能、比結合能、質(zhì)能方程、質(zhì)量虧損”的正確認識
    文理導航(2015年26期)2015-09-29 14:12:24
    梯狀光晶格中自旋軌道耦合的排斥費米氣體
    三維復式晶格的熱容研究
    1024视频免费在线观看| 久久久久久久国产电影| 丝袜美足系列| 国产人伦9x9x在线观看 | 天堂8中文在线网| 亚洲av男天堂| 久久精品国产自在天天线| 狠狠婷婷综合久久久久久88av| 高清视频免费观看一区二区| 1024香蕉在线观看| 午夜免费观看性视频| 九九爱精品视频在线观看| 午夜91福利影院| 91午夜精品亚洲一区二区三区| 久久国产亚洲av麻豆专区| 国产精品 欧美亚洲| 久久久久久久精品精品| 亚洲成人手机| 国产免费一区二区三区四区乱码| 亚洲一区中文字幕在线| 一区二区三区乱码不卡18| 各种免费的搞黄视频| 天天躁夜夜躁狠狠躁躁| 9热在线视频观看99| 亚洲四区av| 我要看黄色一级片免费的| av国产久精品久网站免费入址| 久久久精品94久久精品| 爱豆传媒免费全集在线观看| 日韩不卡一区二区三区视频在线| 少妇被粗大猛烈的视频| 国产精品无大码| 超色免费av| 国产极品天堂在线| 一区在线观看完整版| 2022亚洲国产成人精品| 伊人亚洲综合成人网| 亚洲欧美一区二区三区黑人 | 亚洲第一区二区三区不卡| 成人亚洲精品一区在线观看| 91精品三级在线观看| kizo精华| 亚洲男人天堂网一区| 热99久久久久精品小说推荐| 秋霞伦理黄片| 免费少妇av软件| 精品亚洲成国产av| 精品99又大又爽又粗少妇毛片| 亚洲精品中文字幕在线视频| 亚洲精品第二区| 免费黄色在线免费观看| 国产黄色视频一区二区在线观看| 亚洲国产色片| 色婷婷久久久亚洲欧美| 在线 av 中文字幕| 欧美人与性动交α欧美软件| 人妻人人澡人人爽人人| 日韩免费高清中文字幕av| 不卡av一区二区三区| 五月天丁香电影| 久久韩国三级中文字幕| 国产精品一区二区在线不卡| 制服人妻中文乱码| 久久午夜综合久久蜜桃| 国产片内射在线| 午夜久久久在线观看| 一级片'在线观看视频| 亚洲精品乱久久久久久| 一本—道久久a久久精品蜜桃钙片| 99热网站在线观看| 桃花免费在线播放| 日韩中文字幕视频在线看片| 18禁裸乳无遮挡动漫免费视频| 啦啦啦在线观看免费高清www| 亚洲综合色网址| 亚洲内射少妇av| 大码成人一级视频| 色婷婷av一区二区三区视频| xxx大片免费视频| 一边亲一边摸免费视频| 欧美日韩亚洲国产一区二区在线观看 | 日韩制服骚丝袜av| av免费观看日本| 久久精品国产自在天天线| 亚洲av成人精品一二三区| 菩萨蛮人人尽说江南好唐韦庄| 精品国产一区二区久久| 国产成人免费观看mmmm| 777米奇影视久久| 男女边吃奶边做爰视频| 日韩av免费高清视频| 国产精品 欧美亚洲| 日韩成人av中文字幕在线观看| 国产亚洲最大av| 亚洲国产精品一区二区三区在线| 黑丝袜美女国产一区| 中文字幕另类日韩欧美亚洲嫩草| 电影成人av| 日本av免费视频播放| 美女主播在线视频| 丝袜美腿诱惑在线| 波野结衣二区三区在线| 国产日韩欧美在线精品| 赤兔流量卡办理| 国产精品.久久久| 丁香六月天网| 精品亚洲成国产av| 午夜福利视频在线观看免费| 国产日韩欧美在线精品| 成年女人在线观看亚洲视频| 18禁国产床啪视频网站| 亚洲在久久综合| 看十八女毛片水多多多| 最近中文字幕2019免费版| 永久免费av网站大全| 午夜日本视频在线| 26uuu在线亚洲综合色| 日日撸夜夜添| 91久久精品国产一区二区三区| 日韩成人av中文字幕在线观看| 嫩草影院入口| 精品卡一卡二卡四卡免费| 美国免费a级毛片| 亚洲四区av| 亚洲精品成人av观看孕妇| 亚洲欧美日韩另类电影网站| 性色av一级| 丝袜在线中文字幕| 国产深夜福利视频在线观看| 男人添女人高潮全过程视频| h视频一区二区三区| 国产又色又爽无遮挡免| 亚洲久久久国产精品| www.av在线官网国产| 亚洲一码二码三码区别大吗| 欧美最新免费一区二区三区| 国产黄色视频一区二区在线观看| 一级黄片播放器| 999精品在线视频| 久久精品国产亚洲av天美| 天天躁狠狠躁夜夜躁狠狠躁| 日韩精品免费视频一区二区三区| 亚洲,欧美精品.| 久久精品久久久久久久性| 精品亚洲成a人片在线观看| 久久精品久久久久久久性| 精品国产一区二区久久| av在线老鸭窝| 亚洲精华国产精华液的使用体验| 黄频高清免费视频| 菩萨蛮人人尽说江南好唐韦庄| 青春草亚洲视频在线观看| 亚洲,欧美,日韩| 久久久a久久爽久久v久久| 最新的欧美精品一区二区| 色网站视频免费| 电影成人av| 青草久久国产| 国产成人免费观看mmmm| 97在线视频观看| 观看美女的网站| 搡老乐熟女国产| 精品卡一卡二卡四卡免费| 五月伊人婷婷丁香| 人妻 亚洲 视频| 亚洲综合精品二区| 一区二区三区乱码不卡18| 一区二区三区激情视频| 免费高清在线观看视频在线观看| 国产高清不卡午夜福利| 人人妻人人澡人人爽人人夜夜| 久久免费观看电影| 99热国产这里只有精品6| 国产男人的电影天堂91| 观看美女的网站| 成年人午夜在线观看视频| 啦啦啦啦在线视频资源| 日本91视频免费播放| 久久毛片免费看一区二区三区| 日本黄色日本黄色录像| 亚洲成人av在线免费| 一级片免费观看大全| 久久精品国产自在天天线| 亚洲av福利一区| 天美传媒精品一区二区| av不卡在线播放| 久久久久久久大尺度免费视频| 美女中出高潮动态图| 国产欧美亚洲国产| 边亲边吃奶的免费视频| 91aial.com中文字幕在线观看| 国产一区亚洲一区在线观看| 校园人妻丝袜中文字幕| 亚洲一区中文字幕在线| 国产精品麻豆人妻色哟哟久久| 一区福利在线观看| 2021少妇久久久久久久久久久| 亚洲三区欧美一区| 国产色婷婷99| 国产成人精品无人区| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品电影小说| 婷婷色av中文字幕| 宅男免费午夜| 国产精品偷伦视频观看了| 欧美国产精品一级二级三级| 下体分泌物呈黄色| 国产淫语在线视频| 交换朋友夫妻互换小说| 欧美日韩亚洲国产一区二区在线观看 | 国产精品国产av在线观看| 亚洲av日韩在线播放| 精品久久蜜臀av无| 久久久久网色| 国产成人精品一,二区| 欧美激情高清一区二区三区 | 国产男人的电影天堂91| 亚洲国产色片| 人妻 亚洲 视频| 亚洲av中文av极速乱| 亚洲欧洲精品一区二区精品久久久 | 亚洲av成人精品一二三区| 丝袜美足系列| 丁香六月天网| 男女免费视频国产| 中文精品一卡2卡3卡4更新| 亚洲一区中文字幕在线| 考比视频在线观看| 汤姆久久久久久久影院中文字幕| 成人影院久久| 蜜桃国产av成人99| 日韩,欧美,国产一区二区三区| 日韩中字成人| 啦啦啦在线免费观看视频4| 男人舔女人的私密视频| 亚洲成人av在线免费| 一本—道久久a久久精品蜜桃钙片| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品久久成人aⅴ小说| 国产精品 国内视频| 久久影院123| 亚洲精品中文字幕在线视频| 美女高潮到喷水免费观看| 亚洲伊人久久精品综合| 日韩免费高清中文字幕av| 欧美在线黄色| 天天躁夜夜躁狠狠久久av| 青春草亚洲视频在线观看| 久热久热在线精品观看| 亚洲经典国产精华液单| 亚洲成国产人片在线观看| 国产精品熟女久久久久浪| 秋霞伦理黄片| 黄片无遮挡物在线观看| 久久国产亚洲av麻豆专区| 国产乱人偷精品视频| 在线天堂中文资源库| 国产精品一国产av| 亚洲精品av麻豆狂野| 精品99又大又爽又粗少妇毛片| 狠狠婷婷综合久久久久久88av| 啦啦啦在线观看免费高清www| 青春草国产在线视频| 在线天堂中文资源库| 成人国语在线视频| 老司机亚洲免费影院| 亚洲少妇的诱惑av| 成年美女黄网站色视频大全免费| √禁漫天堂资源中文www| 久久久国产一区二区| 一区二区日韩欧美中文字幕| 91精品三级在线观看| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 人人妻人人澡人人爽人人夜夜| 天堂8中文在线网| 中文字幕人妻丝袜一区二区 | 在线观看免费日韩欧美大片| 亚洲精品第二区| 久久这里只有精品19| 99热国产这里只有精品6| 制服诱惑二区| a级毛片在线看网站| 久久99一区二区三区| 国产精品.久久久| 国产日韩欧美亚洲二区| 国产探花极品一区二区| 亚洲av男天堂| 国产精品久久久久久精品古装| 国产有黄有色有爽视频| 亚洲色图综合在线观看| 午夜免费男女啪啪视频观看| 爱豆传媒免费全集在线观看| 天天躁夜夜躁狠狠躁躁| 男女啪啪激烈高潮av片| 国产 一区精品| 中文字幕人妻丝袜制服| www.自偷自拍.com| 免费黄网站久久成人精品| 久久精品国产a三级三级三级| 亚洲欧美清纯卡通| 亚洲精品av麻豆狂野| 亚洲 欧美一区二区三区| 王馨瑶露胸无遮挡在线观看| av电影中文网址| 大片电影免费在线观看免费| 久久久久精品性色| 大陆偷拍与自拍| 国产精品女同一区二区软件| 久久婷婷青草| 亚洲av.av天堂| 成人二区视频| 国产成人精品无人区| av卡一久久| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 最近最新中文字幕免费大全7| 亚洲激情五月婷婷啪啪| 最近最新中文字幕大全免费视频 | 国产精品免费大片| www.熟女人妻精品国产| 午夜激情久久久久久久| 哪个播放器可以免费观看大片| 亚洲国产av影院在线观看| 少妇熟女欧美另类| 天天操日日干夜夜撸| 毛片一级片免费看久久久久| 在线精品无人区一区二区三| 中文字幕最新亚洲高清| 考比视频在线观看| xxx大片免费视频| 丝瓜视频免费看黄片| 免费人妻精品一区二区三区视频| 97人妻天天添夜夜摸| 黑丝袜美女国产一区| 观看美女的网站| 午夜免费观看性视频| 十分钟在线观看高清视频www| 丝袜在线中文字幕| 大香蕉久久成人网| 久久人人爽人人片av| 亚洲欧美一区二区三区黑人 | 肉色欧美久久久久久久蜜桃| 色播在线永久视频| 美女高潮到喷水免费观看| 美女大奶头黄色视频| 丝袜美腿诱惑在线| 国产精品国产av在线观看| av不卡在线播放| 9191精品国产免费久久| 久久97久久精品| 欧美日韩一级在线毛片| 亚洲综合色惰| 午夜福利视频精品| 欧美成人午夜免费资源| 国产成人精品福利久久| 午夜福利在线观看免费完整高清在| 亚洲精品自拍成人| av网站在线播放免费| 在现免费观看毛片| 丝袜喷水一区| 狠狠婷婷综合久久久久久88av| 美女国产视频在线观看| 啦啦啦在线免费观看视频4| 日韩一区二区三区影片| 国产精品香港三级国产av潘金莲 | 午夜日本视频在线| 丰满乱子伦码专区| 亚洲欧美日韩另类电影网站| 少妇人妻精品综合一区二区| 男女无遮挡免费网站观看| 丝袜人妻中文字幕| 天堂8中文在线网| 99久久人妻综合| 丝袜人妻中文字幕| 免费日韩欧美在线观看| 久久久久精品人妻al黑| 国产激情久久老熟女| 久热久热在线精品观看| 边亲边吃奶的免费视频| 我的亚洲天堂| 精品亚洲成a人片在线观看| 你懂的网址亚洲精品在线观看| 国产精品三级大全| 如何舔出高潮| 熟女电影av网| 欧美另类一区| 99热全是精品| 宅男免费午夜| 久久99蜜桃精品久久| 青春草亚洲视频在线观看| 飞空精品影院首页| 欧美最新免费一区二区三区| 国产一区亚洲一区在线观看| 看非洲黑人一级黄片| 亚洲婷婷狠狠爱综合网| 亚洲国产看品久久| av在线观看视频网站免费| 国产免费一区二区三区四区乱码| 久久亚洲国产成人精品v| 考比视频在线观看| 90打野战视频偷拍视频| 男女高潮啪啪啪动态图| 午夜精品国产一区二区电影| 欧美xxⅹ黑人| 天美传媒精品一区二区| 精品一区二区三卡| 日韩不卡一区二区三区视频在线| 99九九在线精品视频| 另类亚洲欧美激情| 最新的欧美精品一区二区| 久久97久久精品| 成年女人在线观看亚洲视频| 高清av免费在线| 精品视频人人做人人爽| 欧美国产精品一级二级三级| 国产精品不卡视频一区二区| 观看美女的网站| 91精品三级在线观看| xxxhd国产人妻xxx| 色94色欧美一区二区| 亚洲婷婷狠狠爱综合网| 天堂中文最新版在线下载| 亚洲一区中文字幕在线| 亚洲精品视频女| 久久人人爽人人片av| 成人手机av| 婷婷成人精品国产| 日韩精品有码人妻一区| 97在线人人人人妻| 精品久久蜜臀av无| 成人午夜精彩视频在线观看| 久久热在线av| 国产欧美亚洲国产| 国产免费一区二区三区四区乱码| 在线看a的网站| av女优亚洲男人天堂| 亚洲一区二区三区欧美精品| 国产精品久久久久久久久免| 亚洲国产av新网站| 欧美精品高潮呻吟av久久| 啦啦啦在线观看免费高清www| 国产淫语在线视频| 99久久中文字幕三级久久日本| 美女视频免费永久观看网站| 青春草视频在线免费观看| 人体艺术视频欧美日本| 国产av国产精品国产| 国产野战对白在线观看| 视频区图区小说| 亚洲欧美日韩另类电影网站| 国产免费福利视频在线观看| 免费观看性生交大片5| 青春草视频在线免费观看| 26uuu在线亚洲综合色| 狂野欧美激情性bbbbbb| 精品酒店卫生间| 亚洲内射少妇av| 免费黄网站久久成人精品| 久久99精品国语久久久| 夜夜骑夜夜射夜夜干| 欧美日本中文国产一区发布| 男女边吃奶边做爰视频| 久久久a久久爽久久v久久| 26uuu在线亚洲综合色| 男女午夜视频在线观看| 97在线视频观看| 男女啪啪激烈高潮av片| 国产日韩欧美视频二区| 五月伊人婷婷丁香| 丝袜喷水一区| 搡老乐熟女国产| 国产探花极品一区二区| 一本久久精品| 电影成人av| 婷婷成人精品国产| 国产97色在线日韩免费| 午夜免费男女啪啪视频观看| 国产1区2区3区精品| 精品国产一区二区三区久久久樱花| 成年人午夜在线观看视频| 午夜日本视频在线| 水蜜桃什么品种好| 欧美日韩精品成人综合77777| 亚洲伊人色综图| 久久久国产一区二区| kizo精华| 卡戴珊不雅视频在线播放| 永久网站在线| 黄网站色视频无遮挡免费观看| 欧美日韩精品成人综合77777| 在线观看人妻少妇| 男女高潮啪啪啪动态图| av有码第一页| 亚洲精品日本国产第一区| 最新的欧美精品一区二区| 18禁观看日本| xxx大片免费视频| 欧美激情 高清一区二区三区| 91精品伊人久久大香线蕉| 日本-黄色视频高清免费观看| 久久精品夜色国产| 久久精品国产自在天天线| 国产亚洲精品第一综合不卡| 久久久精品94久久精品| 免费大片黄手机在线观看| 美国免费a级毛片| 母亲3免费完整高清在线观看 | 国产成人精品一,二区| 成人漫画全彩无遮挡| 久久久久久久久久久免费av| 国产 精品1| 男人操女人黄网站| 在线观看人妻少妇| 丰满迷人的少妇在线观看| 老司机影院毛片| 人妻少妇偷人精品九色| 久久精品久久久久久噜噜老黄| 国产免费又黄又爽又色| 看免费av毛片| 久久99一区二区三区| 999精品在线视频| 一本色道久久久久久精品综合| 少妇的逼水好多| 少妇猛男粗大的猛烈进出视频| 亚洲欧美一区二区三区黑人 | 日韩制服骚丝袜av| av视频免费观看在线观看| av在线播放精品| 免费女性裸体啪啪无遮挡网站| 秋霞在线观看毛片| 国产成人91sexporn| 午夜免费男女啪啪视频观看| 久久午夜福利片| 一级片免费观看大全| 天天躁夜夜躁狠狠躁躁| 国产亚洲午夜精品一区二区久久| 成年人免费黄色播放视频| 午夜福利影视在线免费观看| 亚洲一区二区三区欧美精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲内射少妇av| 尾随美女入室| 欧美人与性动交α欧美软件| 国产精品成人在线| 在线观看人妻少妇| 黄片播放在线免费| 精品国产乱码久久久久久男人| 少妇 在线观看| 国产麻豆69| 精品人妻在线不人妻| 日日爽夜夜爽网站| 老鸭窝网址在线观看| 三上悠亚av全集在线观看| 国产成人精品婷婷| 欧美+日韩+精品| 我要看黄色一级片免费的| 青草久久国产| av在线老鸭窝| 欧美日韩视频高清一区二区三区二| 叶爱在线成人免费视频播放| 性色av一级| 久久久久国产一级毛片高清牌| 日本黄色日本黄色录像| 五月天丁香电影| 日韩伦理黄色片| 日韩电影二区| 啦啦啦视频在线资源免费观看| 伊人亚洲综合成人网| 亚洲av综合色区一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲在久久综合| 亚洲av免费高清在线观看| 成人手机av| 国产色婷婷99| 午夜福利视频精品| 亚洲精品久久久久久婷婷小说| 久久精品国产鲁丝片午夜精品| 成人漫画全彩无遮挡| 国产成人精品久久久久久| 亚洲国产成人一精品久久久| 日韩av免费高清视频| 人人妻人人澡人人爽人人夜夜| 99国产精品免费福利视频| 一二三四中文在线观看免费高清| 波野结衣二区三区在线| 国产精品女同一区二区软件| 久久精品夜色国产| 亚洲精品久久久久久婷婷小说| 国产欧美日韩一区二区三区在线| 一级毛片 在线播放| 国产精品免费视频内射| 国产精品一国产av| 黄色视频在线播放观看不卡| 亚洲精品日韩在线中文字幕| 如何舔出高潮| 成人国产av品久久久| 男女午夜视频在线观看| 中文字幕色久视频| 亚洲欧美精品综合一区二区三区 | www.熟女人妻精品国产| 少妇的逼水好多| 久久人人爽av亚洲精品天堂| 免费大片黄手机在线观看| 亚洲婷婷狠狠爱综合网| 成年美女黄网站色视频大全免费| 中文乱码字字幕精品一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲欧美精品综合一区二区三区 | 成人18禁高潮啪啪吃奶动态图| 伦理电影大哥的女人| 国产精品国产三级国产专区5o| 在线观看三级黄色| 久久精品人人爽人人爽视色| 激情视频va一区二区三区|