• <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
    梯狀光晶格中自旋軌道耦合的排斥費米氣體
    三維復式晶格的熱容研究
    岛国毛片在线播放| 日韩精品免费视频一区二区三区| 亚洲午夜精品一区,二区,三区| 免费在线观看视频国产中文字幕亚洲 | 建设人人有责人人尽责人人享有的| 亚洲伊人色综图| 欧美黑人欧美精品刺激| 色精品久久人妻99蜜桃| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 免费黄频网站在线观看国产| 免费少妇av软件| 丝袜脚勾引网站| 国产色视频综合| 男女下面插进去视频免费观看| a 毛片基地| 热re99久久精品国产66热6| 99九九在线精品视频| 欧美乱码精品一区二区三区| 99国产精品一区二区三区| 下体分泌物呈黄色| 侵犯人妻中文字幕一二三四区| 日韩 亚洲 欧美在线| 日韩一区二区三区影片| 免费在线观看影片大全网站 | 国产免费又黄又爽又色| 好男人电影高清在线观看| 丝袜喷水一区| 天堂中文最新版在线下载| 视频区图区小说| 黄色毛片三级朝国网站| 亚洲伊人色综图| 人人妻人人澡人人爽人人夜夜| 99热国产这里只有精品6| 高潮久久久久久久久久久不卡| 亚洲欧美激情在线| 十八禁网站网址无遮挡| 在线观看免费午夜福利视频| 2018国产大陆天天弄谢| 日日夜夜操网爽| 妹子高潮喷水视频| 久久鲁丝午夜福利片| 国产av精品麻豆| 久久99精品国语久久久| 18禁裸乳无遮挡动漫免费视频| 精品一区在线观看国产| 国产精品久久久久久精品古装| 亚洲国产精品999| 亚洲av国产av综合av卡| 一级毛片女人18水好多 | 精品国产一区二区久久| 国产成人精品久久二区二区免费| 午夜精品国产一区二区电影| 老司机影院毛片| 日本一区二区免费在线视频| 91九色精品人成在线观看| 国产黄色视频一区二区在线观看| 另类亚洲欧美激情| 国产精品国产三级专区第一集| 亚洲精品国产色婷婷电影| 久久国产精品影院| 国产欧美日韩精品亚洲av| 亚洲美女黄色视频免费看| 亚洲情色 制服丝袜| 日韩一区二区三区影片| 高清视频免费观看一区二区| 免费在线观看黄色视频的| 久久青草综合色| 丝袜在线中文字幕| 成人亚洲欧美一区二区av| 极品人妻少妇av视频| 精品国产国语对白av| 午夜免费成人在线视频| 国产福利在线免费观看视频| 搡老乐熟女国产| 国产亚洲av高清不卡| 蜜桃在线观看..| 汤姆久久久久久久影院中文字幕| 青春草亚洲视频在线观看| 国产视频一区二区在线看| 亚洲视频免费观看视频| 又大又爽又粗| 国产精品亚洲av一区麻豆| 亚洲国产av新网站| 成人亚洲精品一区在线观看| 亚洲精品久久久久久婷婷小说| 色综合欧美亚洲国产小说| 国产99久久九九免费精品| 女性被躁到高潮视频| 午夜精品国产一区二区电影| 少妇精品久久久久久久| 高清视频免费观看一区二区| 51午夜福利影视在线观看| 日韩av在线免费看完整版不卡| 后天国语完整版免费观看| 日韩中文字幕视频在线看片| 丰满人妻熟妇乱又伦精品不卡| 精品少妇内射三级| 亚洲五月色婷婷综合| 三上悠亚av全集在线观看| 在线观看免费视频网站a站| 亚洲中文av在线| a级毛片在线看网站| 亚洲国产av新网站| 美国免费a级毛片| 亚洲欧美色中文字幕在线| 久久99一区二区三区| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 亚洲精品久久久久久婷婷小说| 人人澡人人妻人| 欧美日韩亚洲高清精品| 亚洲人成网站在线观看播放| 80岁老熟妇乱子伦牲交| 国产99久久九九免费精品| 亚洲精品av麻豆狂野| 少妇猛男粗大的猛烈进出视频| www日本在线高清视频| 啦啦啦啦在线视频资源| 亚洲国产毛片av蜜桃av| 日韩欧美一区视频在线观看| 午夜福利一区二区在线看| 老汉色av国产亚洲站长工具| 精品国产乱码久久久久久男人| 一区二区三区精品91| 久久久精品国产亚洲av高清涩受| 欧美国产精品一级二级三级| 看免费av毛片| 侵犯人妻中文字幕一二三四区| 女人高潮潮喷娇喘18禁视频| kizo精华| 国产熟女欧美一区二区| 日韩伦理黄色片| 欧美国产精品一级二级三级| 亚洲九九香蕉| 自线自在国产av| 另类亚洲欧美激情| 亚洲精品美女久久久久99蜜臀 | 国产一区二区在线观看av| 成人影院久久| av电影中文网址| 欧美日韩亚洲高清精品| 80岁老熟妇乱子伦牲交| 老司机影院毛片| 精品熟女少妇八av免费久了| 夫妻午夜视频| 91九色精品人成在线观看| 蜜桃国产av成人99| 午夜两性在线视频| 欧美精品高潮呻吟av久久| 国产xxxxx性猛交| 亚洲av美国av| 欧美日韩福利视频一区二区| 手机成人av网站| 一区福利在线观看| 国产成人精品久久久久久| 狂野欧美激情性bbbbbb| 看免费av毛片| 色婷婷久久久亚洲欧美| 国产爽快片一区二区三区| 国产免费又黄又爽又色| 十八禁网站网址无遮挡| 欧美精品一区二区大全| 青草久久国产| 各种免费的搞黄视频| 91老司机精品| 男人爽女人下面视频在线观看| 欧美+亚洲+日韩+国产| 一本一本久久a久久精品综合妖精| 中文字幕另类日韩欧美亚洲嫩草| 亚洲美女黄色视频免费看| 久久午夜综合久久蜜桃| 99国产精品一区二区蜜桃av | av又黄又爽大尺度在线免费看| av在线播放精品| 一级黄色大片毛片| 国产又色又爽无遮挡免| 国产97色在线日韩免费| 国产成人欧美| 久久天躁狠狠躁夜夜2o2o | 99国产综合亚洲精品| 一边亲一边摸免费视频| 亚洲图色成人| 青草久久国产| 高清欧美精品videossex| 蜜桃在线观看..| 国产成人一区二区三区免费视频网站 | 天天躁夜夜躁狠狠躁躁| 日韩一本色道免费dvd| 一区二区三区四区激情视频| 国产av一区二区精品久久| 精品熟女少妇八av免费久了| 久久国产精品影院| 久久人妻熟女aⅴ| 搡老乐熟女国产| 亚洲伊人色综图| 亚洲精品自拍成人| 午夜福利影视在线免费观看| 亚洲成人国产一区在线观看 | 一区二区日韩欧美中文字幕| 免费在线观看黄色视频的| 亚洲精品日本国产第一区| 亚洲男人天堂网一区| 色婷婷av一区二区三区视频| 婷婷色av中文字幕| 在线看a的网站| 久久青草综合色| 巨乳人妻的诱惑在线观看| 国产精品成人在线| 在线 av 中文字幕| 美女高潮到喷水免费观看| 一边亲一边摸免费视频| 男女午夜视频在线观看| 国产黄频视频在线观看| 手机成人av网站| 国产在线免费精品| 777米奇影视久久| 老汉色∧v一级毛片| 51午夜福利影视在线观看| 高潮久久久久久久久久久不卡| 国产在线观看jvid| 亚洲欧洲国产日韩| 妹子高潮喷水视频| 男女边摸边吃奶| 王馨瑶露胸无遮挡在线观看| 日本黄色日本黄色录像| www.熟女人妻精品国产| 国产爽快片一区二区三区| 伊人亚洲综合成人网| 亚洲,一卡二卡三卡| 国产精品香港三级国产av潘金莲 | 男人添女人高潮全过程视频| 亚洲国产av新网站| 欧美 日韩 精品 国产| 国产成人精品无人区| 国精品久久久久久国模美| 午夜免费观看性视频| 久久久精品国产亚洲av高清涩受| 亚洲五月色婷婷综合| 精品国产一区二区三区久久久樱花| 欧美大码av| 天天躁狠狠躁夜夜躁狠狠躁| 少妇人妻 视频| 欧美成人午夜精品| 男女下面插进去视频免费观看| 高清不卡的av网站| 尾随美女入室| 爱豆传媒免费全集在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲成色77777| 久久性视频一级片| www.熟女人妻精品国产| 又黄又粗又硬又大视频| 80岁老熟妇乱子伦牲交| 国产精品 国内视频| 大码成人一级视频| 亚洲国产精品成人久久小说| 免费看不卡的av| 后天国语完整版免费观看| 好男人视频免费观看在线| 少妇精品久久久久久久| 欧美黄色片欧美黄色片| 亚洲熟女精品中文字幕| 午夜日韩欧美国产| 精品人妻1区二区| 亚洲九九香蕉| 日韩制服骚丝袜av| 国产1区2区3区精品| 妹子高潮喷水视频| 女性生殖器流出的白浆| 免费不卡黄色视频| 国产精品99久久99久久久不卡| 欧美激情极品国产一区二区三区| 国产精品人妻久久久影院| 成人影院久久| 97精品久久久久久久久久精品| 国产一卡二卡三卡精品| 久久久久久久久久久久大奶| 日韩av免费高清视频| 9191精品国产免费久久| 纯流量卡能插随身wifi吗| h视频一区二区三区| 人人妻人人澡人人看| 亚洲专区国产一区二区| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 久久精品亚洲熟妇少妇任你| 欧美老熟妇乱子伦牲交| 麻豆av在线久日| 丰满迷人的少妇在线观看| 欧美日韩黄片免| av在线老鸭窝| 亚洲七黄色美女视频| 中文字幕色久视频| 欧美国产精品va在线观看不卡| 青春草视频在线免费观看| 麻豆av在线久日| 久久综合国产亚洲精品| 黑人巨大精品欧美一区二区蜜桃| 90打野战视频偷拍视频| 日韩中文字幕视频在线看片| 777久久人妻少妇嫩草av网站| 成人午夜精彩视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲欧美日韩高清在线视频 | 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 亚洲欧洲日产国产| h视频一区二区三区| 国产精品亚洲av一区麻豆| 久久久久久久国产电影| av福利片在线| 超碰97精品在线观看| 久久天躁狠狠躁夜夜2o2o | 国产一区二区三区av在线| 一本综合久久免费| 久久久久精品人妻al黑| 我的亚洲天堂| 精品少妇黑人巨大在线播放| 亚洲欧洲日产国产| 久热这里只有精品99| 亚洲免费av在线视频| 纯流量卡能插随身wifi吗| 中文欧美无线码| 五月天丁香电影| 国产有黄有色有爽视频| av电影中文网址| 国产精品一区二区在线观看99| 一级毛片女人18水好多 | 99九九在线精品视频| 国产一区二区激情短视频 | 亚洲精品中文字幕在线视频| 亚洲天堂av无毛| 久久人人爽人人片av| 电影成人av| 手机成人av网站| 久久99一区二区三区| 99国产精品一区二区三区| 久久人妻福利社区极品人妻图片 | 天天躁夜夜躁狠狠久久av| av不卡在线播放| 国产高清不卡午夜福利| 一本综合久久免费| 天堂中文最新版在线下载| 天天躁日日躁夜夜躁夜夜| av有码第一页| 色播在线永久视频| 丝袜美腿诱惑在线| 欧美人与性动交α欧美软件| 大话2 男鬼变身卡| 2021少妇久久久久久久久久久| 精品国产乱码久久久久久男人| 精品卡一卡二卡四卡免费| av有码第一页| 王馨瑶露胸无遮挡在线观看| 欧美日韩综合久久久久久| 日日摸夜夜添夜夜爱| 19禁男女啪啪无遮挡网站| 久久99热这里只频精品6学生| 久久av网站| 成人影院久久| 中文字幕人妻丝袜一区二区| 国产99久久九九免费精品| 国产欧美日韩精品亚洲av| 人人妻人人添人人爽欧美一区卜| 91精品三级在线观看| 欧美日韩视频高清一区二区三区二| 国产不卡av网站在线观看| 久久久精品94久久精品| 国产欧美日韩一区二区三 | 少妇猛男粗大的猛烈进出视频| 国产欧美亚洲国产| 热re99久久国产66热| 久久女婷五月综合色啪小说| 日韩免费高清中文字幕av| 丝袜美腿诱惑在线| www日本在线高清视频| 国产成人精品久久二区二区91| 久久久久网色| 日韩人妻精品一区2区三区| 午夜激情av网站| 日韩一本色道免费dvd| 狂野欧美激情性bbbbbb| 超色免费av| 久久久精品区二区三区| 日本色播在线视频| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| 少妇粗大呻吟视频| 视频区图区小说| av天堂在线播放| 成人国产av品久久久| 一本大道久久a久久精品| 香蕉丝袜av| 一本大道久久a久久精品| 精品免费久久久久久久清纯 | 久久性视频一级片| 亚洲国产成人一精品久久久| 免费不卡黄色视频| 男人添女人高潮全过程视频| 国产视频首页在线观看| 国产成人一区二区在线| 亚洲精品av麻豆狂野| 亚洲av日韩精品久久久久久密 | 国产在线视频一区二区| 精品卡一卡二卡四卡免费| 亚洲九九香蕉| 亚洲一区中文字幕在线| 汤姆久久久久久久影院中文字幕| 老司机靠b影院| 精品熟女少妇八av免费久了| av在线app专区| 爱豆传媒免费全集在线观看| 人妻 亚洲 视频| 一级毛片电影观看| 亚洲第一青青草原| 啦啦啦视频在线资源免费观看| 成人亚洲精品一区在线观看| 老司机深夜福利视频在线观看 | 亚洲欧洲日产国产| 尾随美女入室| 一级毛片电影观看| 国产在线观看jvid| 亚洲国产av新网站| 亚洲精品第二区| 叶爱在线成人免费视频播放| 久久影院123| 99久久人妻综合| 色精品久久人妻99蜜桃| 欧美日韩国产mv在线观看视频| 色94色欧美一区二区| 看免费av毛片| av有码第一页| 高潮久久久久久久久久久不卡| 人人妻人人爽人人添夜夜欢视频| 国产高清videossex| 男人添女人高潮全过程视频| 国产一区亚洲一区在线观看| 男男h啪啪无遮挡| av有码第一页| 精品一区二区三区av网在线观看 | 久久免费观看电影| 精品卡一卡二卡四卡免费| 午夜免费男女啪啪视频观看| 国产成人欧美在线观看 | 天天影视国产精品| 麻豆国产av国片精品| 免费看av在线观看网站| cao死你这个sao货| 久久精品国产综合久久久| 亚洲第一青青草原| 免费不卡黄色视频| 日韩免费高清中文字幕av| 亚洲久久久国产精品| 午夜激情久久久久久久| 曰老女人黄片| 成人亚洲欧美一区二区av| 亚洲,欧美,日韩| 国产黄色免费在线视频| 国产人伦9x9x在线观看| 精品人妻一区二区三区麻豆| 中文欧美无线码| 不卡av一区二区三区| 午夜福利在线免费观看网站| 中文字幕av电影在线播放| 亚洲国产精品999| 欧美少妇被猛烈插入视频| 精品亚洲乱码少妇综合久久| 国产精品久久久久久人妻精品电影 | 成在线人永久免费视频| 91精品伊人久久大香线蕉| 久久精品人人爽人人爽视色| 成人18禁高潮啪啪吃奶动态图| 精品人妻一区二区三区麻豆| 欧美日韩成人在线一区二区| 婷婷色综合www| 国产免费现黄频在线看| 天堂8中文在线网| 欧美老熟妇乱子伦牲交| 激情五月婷婷亚洲| 国产精品久久久人人做人人爽| 免费av中文字幕在线| 精品人妻在线不人妻| 一级毛片电影观看| 久久亚洲精品不卡| 777米奇影视久久| h视频一区二区三区| 伦理电影免费视频| 成年动漫av网址| 国产精品三级大全| 午夜免费男女啪啪视频观看| 欧美大码av| 男女床上黄色一级片免费看| 天堂8中文在线网| 久久鲁丝午夜福利片| 国精品久久久久久国模美| 首页视频小说图片口味搜索 | www.av在线官网国产| 中国国产av一级| 国产午夜精品一二区理论片| 国产成人免费观看mmmm| 国产精品国产三级专区第一集| 日韩av在线免费看完整版不卡| 国产色视频综合| 天天躁夜夜躁狠狠躁躁| 一级毛片黄色毛片免费观看视频| 亚洲国产精品一区三区| svipshipincom国产片| 国产亚洲一区二区精品| 国产亚洲精品久久久久5区| 国产成人精品久久二区二区免费| 午夜免费观看性视频| 亚洲av成人精品一二三区| 一级a爱视频在线免费观看| 成人免费观看视频高清| 色综合欧美亚洲国产小说| www.自偷自拍.com| 一级,二级,三级黄色视频| 咕卡用的链子| 亚洲,欧美,日韩| 多毛熟女@视频| 精品亚洲成国产av| 这个男人来自地球电影免费观看| 9191精品国产免费久久| 91精品三级在线观看| 男女边吃奶边做爰视频| 亚洲av成人不卡在线观看播放网 | 国产色视频综合| 水蜜桃什么品种好| 中文字幕人妻丝袜制服| 日日摸夜夜添夜夜爱| 亚洲国产精品999| 成年女人毛片免费观看观看9 | 我的亚洲天堂| 日韩欧美一区视频在线观看| 操美女的视频在线观看| 国产一级毛片在线| 国产成人影院久久av| a级毛片黄视频| 成年美女黄网站色视频大全免费| 日日摸夜夜添夜夜爱| www.999成人在线观看| 久久午夜综合久久蜜桃| 男人添女人高潮全过程视频| 精品久久久久久电影网| 老司机靠b影院| 久久99精品国语久久久| 中文乱码字字幕精品一区二区三区| 色94色欧美一区二区| 久久久久久久精品精品| 热99久久久久精品小说推荐| 精品久久蜜臀av无| 日本五十路高清| 久久久精品国产亚洲av高清涩受| 妹子高潮喷水视频| 日韩av在线免费看完整版不卡| 久久久久网色| 欧美激情高清一区二区三区| 在线观看免费高清a一片| 国产欧美亚洲国产| 久久综合国产亚洲精品| 亚洲精品久久成人aⅴ小说| 国产精品人妻久久久影院| 97在线人人人人妻| 一区二区av电影网| www.自偷自拍.com| 亚洲国产av影院在线观看| 晚上一个人看的免费电影| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 日韩av不卡免费在线播放| 五月天丁香电影| 精品免费久久久久久久清纯 | 国产成人精品久久二区二区免费| 久久鲁丝午夜福利片| 99国产精品一区二区三区| 国产精品九九99| 一区二区三区四区激情视频| 国产精品二区激情视频| 在线观看国产h片| kizo精华| 中文字幕最新亚洲高清| 另类亚洲欧美激情| 黑人猛操日本美女一级片| 午夜免费男女啪啪视频观看| 国产成人啪精品午夜网站| 精品一区二区三卡| 免费久久久久久久精品成人欧美视频| 欧美精品高潮呻吟av久久| 久久久久久亚洲精品国产蜜桃av| 精品亚洲乱码少妇综合久久| 免费不卡黄色视频| 极品少妇高潮喷水抽搐| 国产一区二区在线观看av| 欧美+亚洲+日韩+国产| 国产精品秋霞免费鲁丝片| 自线自在国产av| 日韩一卡2卡3卡4卡2021年| 精品亚洲乱码少妇综合久久| 国产成人av教育| 男人舔女人的私密视频| 国产精品秋霞免费鲁丝片| 新久久久久国产一级毛片| 久久久久久亚洲精品国产蜜桃av| 新久久久久国产一级毛片| 国产精品亚洲av一区麻豆| 无遮挡黄片免费观看| 18禁裸乳无遮挡动漫免费视频| 91成人精品电影| 中文字幕色久视频| 美女午夜性视频免费| 免费av中文字幕在线| 久久亚洲国产成人精品v| 黑人欧美特级aaaaaa片|