999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

金屬材料液態與固態結構相關性的研究進展

2017-12-04 10:35:39孫永軍
山東冶金 2017年4期
關鍵詞:結構研究

孫永軍

(濟鋼集團山東建設工程有限公司,山東 濟南250101)

專論綜述

金屬材料液態與固態結構相關性的研究進展

孫永軍

(濟鋼集團山東建設工程有限公司,山東 濟南250101)

綜述了最近幾十年國內外液態結構的研究方法、研究成果以及應用情況。通過高溫X-ray衍射、高溫黏度、計算機模擬等方法已取得了一些有關高溫液態結構的信息,為液態-固態結構轉變、合金優異組織性能的獲取提供理論指導。為獲得更多的熔體結構信息,建議研制高精尖設備測試液態結構,尋找盡可能精確、形式簡單的勢函數對液態結構進行模擬。

金屬材料;液態結構;固態結構;相關性;研究進展

1 前言

金屬材料的加工如鑄造、焊接、粉末冶金等大多需要經歷液-固態轉變過程,液態的結構特征無疑會對固態的組織和性能起到重要的影響[1]。早在20世紀20年代,法國科學家就發現金屬液態結構中存在著和固態結構相同的微小顆粒,認為鑄造中存在“遺傳性”[2];20世紀30年代初及50年代,前蘇聯著名冶金學家A.C.庫什尼爾斯基和C.M.沃羅諾夫亦提出了組織的遺傳性問題[3]。但由于液態結構具有“長程有序”、“短程無序”的結構特征,又像氣體一樣具有流動性,研究其結構和性質較艱難和復雜,人們對液態的認識僅僅局限于氣體、夾雜、微量元素等異質組成對最終組織的影響[4]。直到近幾年才逐漸認識到,即使在純凈的合金液態體系中,液態結構的變化對固態結構的組織性能和質量也存在直接和重要的影響。而從熔體結構控制的角度來改善和控制凝固尚是經驗性的,沒有形成比較系統的理論。實驗發現,液態和固態結構之間的聯系不僅僅是在原子近鄰存在相似形,而且對于某些熔體來說在中程序尺度上也存在關聯[5],認識這種關聯對改善材料的組織性能等具有深遠的意義。本研究綜述了最近幾十年國內外液態結構的研究方法及研究成果,希望通過獲得液態結構的信息以及溫度、元素添加等對液態結構的影響,控制金屬液態預結晶狀態和冷卻速度,達到顯著改善金屬材料的組織、性能和質量的目的。

2 液態結構研究方法

目前主要有三大類方法被采用到金屬液態結構的研究中[6],分別是直接測量法、間接測量法和計算機模擬法。

2.1 直接測量法[6]

應用較多的有中子衍射、電子衍射、X射線衍射和X射線吸收精細結構(X-ray absorption fine structure,XAFS),其基本原理為:根據運動粒子束或電磁波等與熔體內的原子相互作用產生的相干散射來獲取熔體結構內部原子排列信息。一般中子衍射適合測定物質中輕原子的分布;電子衍射適用于測定材料的微細結構;而X射線衍射適用于測定和分析團簇結構和結構有序度,可在很寬的波譜范圍提供大量的信息,如結構因子S(Q)、配位數、徑向分布函數RDF、雙體相關函數g(r)等,此外,相對于電子和中子衍射,X射線衍射所需設備小且操作方便,因而在測量熔體結構的研究中應用更為廣泛。

XAFS是一種同步福射特有的結構分析方法,可提供小范圍內原子團簇的結構信息。與X射線和中子散射相比,XAFS吸收邊由于具有原子特征,調節X射線的能量可獲得合金中指定原子的局域結構,這是對X射線衍射所獲得的團簇信息的有效補充。但由于XAFS對實驗條件要求高,特別是對于液體樣品,在樣品制備與密封方面要求茍刻,很難成功采集到高質量的XAFS譜,并且在數據處理方面十分復雜。因此,國內外進行金屬液態結構XAFS實驗研究的團隊相對較少,主要的幾個研究團隊包括中科院物理所陸坤權課題組和意大Camerino大學Di Cicco課題組。

2.2 間接測量法

金屬液態的物理性質,包括黏度、擴散系數、密度、熱容、磁化率等,與液態的結構密切相關,是液態微觀結構變化的宏觀表現[7]。液態物理性質的變化可間接反映其結構的變換,因此可以用測量金屬液態黏度、密度、電阻等方法來間接研究液態的結構。如黏度是被人們廣泛應用的物理性質,是體系中原子遷移性質的一種表現,反映了原子間結合力的大小,是金屬液態最敏感的物理性質之一。黏度隨溫度的變化規律反映了液態結構的穩定性[8-9];密度與液態中原子的配位數及原子間距有關;電阻率則從更微觀的電子角度來探究金屬的液態結構,當液態結構發生突變時,其值會在相應位置發生突變[10]。

2.3 計算機模擬[11-12]

在實驗無法觸及或一些極端的條件下,計算機模擬無疑成為實驗研究與理論研究的橋梁。計算機模擬可全面給出反映熔體結構信息的參數,如雙體分布函數、結構因子、原子間相互作用勢函數和鍵取向序。因此從20世紀初開始,計算機模擬得到了廣泛的應用,目前原子、電子尺度上應用最廣泛的研究方法是分子動力學模擬(MD)和蒙特卡羅模擬(MC),圖1顯示了MD的研究過程,而分子尺度上的動力模擬法有:等溫等壓分子動力學,非平衡分子動力學,從頭算分子動力學方法和能量最小化方法。分子動力學法是在給定的粒子間相互作用勢下,通過數值求解系統的運動方程組,計算在各瞬時原胞中N個粒子的速度和坐標,由此分析系統的性質,如結構、動力學、熱力學等[12]。

圖1 分子動力學模擬研究液體結構的技術路線[11]

3 理論成果

程素娟等人[13]采用高溫液態X-ray衍射儀研究了 In-1%Cu、In-1%Al和 In-1%Sn(質量分數)合金熔體在280~750℃溫度區間的原子團簇微觀結構的變化。研究發現,在液相線以上較大范圍內隨溫度的升高,熔體的原子團簇仍然存在著微觀熱收縮現象,Cu和Al的加入加強了團簇的微觀熱收縮程度,而Sn的加入降低了團簇的微觀熱收縮程度。

Han等人[14]對三元Co-Gd-Ti系相分離合金的相平衡及凝固行為進行了研究。液態XRD結果表明Co30Gd35Ti35合金在1417~1487℃溫度區間內有兩種液相同時存在;1327~1507℃溫度區間內,除了之前的兩種液相,還有部分CoTi固相析出;繼續降溫,相分離現象消失,只存在一種液相和繼續析出的CoTi相;最終凝固形成Gd、CoTi和CoTi2三種相。同樣以CdTi二元不混溶合金為基體,Co25Gd37.5Ti37.5、Co10Gd45Ti45和Co30Gd30Ti50合金都發生了液相分離,而Co30Gd50T20合金卻沒有,說明是否出現兩種液相共存與焰體中形成的結構單元能否互溶密切相關。通過對相分離合金溶體局域結構的研究有望進一步揭示相分離機理,從而實現對相分離材料制備過程的控制。

孫益民等人[15]采用差示掃描量熱儀,將Al-1.1%Fe合金溶體過熱至液相線以上不同溫度進行冷卻,發現溶體溫度低于900℃時,隨著過熱溫度的升高,其過冷度急劇增大,固態組織中的FeAl3相顯著細化,但超過900℃過熱,FeAl3相則細化不明顯。Li等人[16]通過分子動力學模擬方法研究了常壓、不同溫度條件下液體Ni3Al的結構特征,結果在雙體分布函數曲線上發現了肩峰的存在,經過選取不同的鍵對類型對該曲線進行分析,發現肩峰的強度不受所選系統大小的影響,而是液體中具有較強結合力團簇的體現。這項工作首次為表征液體中的中程有序結構提供了依據。

在實際應用中,一般綜合采用幾種方法對液態結構進行研究。秦敬玉等人[17]應用液態X-ray衍射儀和逆MonteCarlo(RMC)模擬的方法研究了Fe3Si和FeSi合金熔體的預峰并討論了其微觀機制。預峰對應的真實空間的距離與Fe3Si的DO3類型晶體中的Si-Si之間的距離相近,偏差≯2.4%。此外,在Fe3Si合金熔體中只在局部結構SSiSi(Q)中發現預峰,而在FeSi合金熔體中預峰在SSiSi(Q)和SFeFe(Q)出現。兩種合金中的Si原子周圍Fe原子的數目分布成高斯分布,在FeSi合金熔體中視為最主要的,結果表明了FeSi晶體熔化過程中存在Fe7Si和FeSi7型團簇。Fe3Si熔體中的中程有序是由Si-Si團簇引起的,而FeSi熔體中的中程有序是由于Fe-Fe和Si-Si團簇引起的。邊秀房等采用高溫黏度儀和液態X-ray衍射儀研究了 Cu75Al25、Cu87Sn13、Al-12.5%Si及純 Cu、純Sn熔體的黏度和液態結構,研究表明:純金屬和具有短程有序結構的共晶合金在液相線以上很大溫度范圍內遵循Arrihenius定律,而當熔體結構由中程有序向短程有序結構轉變時,Arrihenius曲線上出現裂點,表明黏度的變化是熔體結構改變的一個重要特征[8]。趙巖等人利用高溫液態X-ray衍射儀和高溫黏度儀研究了 In30Sn70[18]、Cu80Ag20[19]、Cu60Sn40[20]二元合金熔體。In30Sn70合金熔體(液相線172℃、固相線120℃)團簇的相關半徑rc,相關長度D和原子數目Nc在溫度降低到液相線時突然增大,表明形核過程中有結構突變,將高溫黏度儀測量的黏度和配位數Nmin結合比得出,溫度和成分引起的不連續結構轉變區域應在145~165℃[18]。Cu80Ag20二元合金熔體(液相線950℃、固相線750℃)團簇的相關半徑和原子數目在低于1100℃時隨著溫度的降低有較小的增長,而相關長度與溫度幾乎呈線性變化關系,表明該合金熔體在凝固過程中結構轉變穩定。由于熔體團簇在形核時起到至關重要的作用,所以Cu80Ag20合金熔體中的團簇可能成為形核時的晶胚,這與In30Sn70合金熔體是兩種不同的形核機制,可以認為前者是穩定型凝固模型,而后者是跳躍型凝固模型[19]。通過高溫X-ray衍射試驗、黏度試驗和差示掃描量熱儀研究表明,在Cu60Sn40二元熔體合金中存在液-液相轉變的現象,在加熱和冷卻該合金時存在著逆轉變[20]:

Kaban等人除了利用高溫液態X-ray衍射儀和高溫黏度儀,還測量了合金的熔體密度(方法見文獻[21]),從而研究了溫度對In20Sn80合金熔體結構的影響[22]。密度隨溫度的變化關系如圖2所示,從圖中可以看出,密度隨溫度的升高而降低,表明熔體的體積隨溫度的升高增加了,可能是原子間距的增加或者是第一峰位向低角度偏移的結果,但是從高溫X-ray衍射曲線中并沒有發現第一峰位的偏移,原子間距還稍微變小了,這就說明該熔體結構是不均勻的。經過綜合分析認為,In20Sn80合金熔體結構是由組成團簇的原子和填充團簇之間空隙的具有高遷移率的自由電子組成的,在較寬的溫度區間內原子密排和團簇的成分是穩定的[22]。這些發現都進一步揭示了金屬液態經結構的本質,豐富了液態結構方面的知識,有助于更進一步了解金屬熔化和凝固的本質。

圖2 液態In20Sn80合金密度隨溫度成線性變化[22]

4 實踐應用

根據金屬遺傳性和液態結構變化規律的原理,俄羅斯學者提出了一種稱為熔體高溫處理(high temperature treatment of themelt,HTTM)的新工藝[23]。通過這種工藝,使合金熔體達到完全均勻化的平衡狀態,從而改善材料和制品的凝固組織、結構和性能。熔體高溫處理工藝包括加熱金屬或合金到某一臨界溫度以上,保持一段時間,然后冷卻到澆注溫度進行澆注的過程。熔體高溫處理工藝與傳統熔煉工藝的比較表明,1850℃熔體高溫處理工藝與傳統1600℃熔煉工藝相比,可顯著提高ZhS3DK-VI合金的力學性能[24]。王海峰等[25]研究發現超溫處理會使新型鎳基單晶高溫熔體的結構發生變化,從而導致溶質分配系數發生變化。

M963是一種W、Mo、Nb含量較高的鑄造鎳基多晶高溫合金,具有較高的高溫強度,但塑性較低。用高溫液態X-ray射線衍射的方法證明該合金經1520℃熔體過熱后熔體中除存在短程有序外,仍然存在Ni3A1型中程有序原子團簇和殘余MC碳化物顆粒(或MC碳化物型中程有序原子團簇),熔體熱歷史不同,相同溫度下測試的衍射譜亦明顯不同,表明M963合金的液態結構與熔體熱歷史之間的密切相關性[26]。將熔體高溫處理工藝(1550~1650℃、常規熔體精煉)應用到M963合金后,975℃/225 MPa條件下的持久壽命和持久塑性明顯提高,并且發現熔體處理對M963合金凝固組織的影響主要表現在晶粒大小、MC碳化物形貌和γ’相的熱穩定性上[26]。

王震等人利用液態金屬X-ray射線衍射儀測試了一種Ni-Cr-Co-W-Mo-T-Al-Ti系鎳基合金熔體結構在1450~1600℃溫度范圍內的變化。發現隨著熔體處理溫度的升高,合金熔體的結構更加均勻;并研究分析了經過不同過熱處理后合金的凝固組織,找到了適合于處理該合金的臨界溫度:1700℃,在此溫度下,形成的凝固組織枝晶干、枝晶間γ’的尺寸均減小,形狀更規則,并且枝晶干、枝晶間γ’尺寸差別減小,枝晶偏析比趨近于1[27]。

非晶態合金具有非常優異和特殊的物理、化學、力學性能,至今為止,仍是科研人員研究的熱點[28]。液態是非晶材料的母態,想了解合金非晶形成能力的本質,可以從了解合金熔體結構開始。邊秀房等人提出了過熱熔體的脆性理論:即合金液相線以上的結構隨溫度的變化是由構型的動力學因素決定的,通過測量合金熔體液相線以上的黏度,經以下公式:

擬和可計算出液體的過熱脆性值(見圖3),M越小,非晶形成能力越好[29]。2007年,邊秀房等人[30]又運用高溫X-ray衍射的方法,通過研究金屬熔體團簇的穩定性提出了預測非晶形成能力的參數:B=dr/dT(r為團簇的相關半徑,T是絕對溫度)。當B<1×10-4nm/K時,金屬熔體可能形成大塊非晶。在這兩個參數的指導下,有望獲得大塊的非晶合金,使其優異的性能得到廣泛的應用,如利用非晶合金的低矯頑力、低損耗、強磁導率和良好的高頻性能等制作的高效節能非晶合金變壓器。

圖3 AlCoCe系統合金的黏度數據和擬和曲線[30]

5 結論

以上研究可見,國內外利用高溫X-ray衍射、高溫黏度、計算機模擬等方法對液態結構進行了大量的研究,取得了一些有關高溫液態結構性質的信息,為液態-固態結構轉變、合金優異組織性能的獲取提供了理論指導,具有重要的實際意義。

但目前獲得的有關液態結構性質的信息量還是很薄弱。國內高溫X-ray衍射儀、高溫黏度儀等測量液態結構性質的高精尖測試設備稀缺,大部分合金液態結構信息的獲得主要靠計算機模擬的方法,但模擬法和實際測量值往往存在著差異,在指導實踐的過程中,需要反復摸索。今后,研制高精尖測試液態結構性質的設備及尋找盡可能精確、形式簡單的勢函數是豐富液態結構信息的熱點。

[1] E KARAK?SE,M KESKIN.Structural investigations ofmechanical properties of Albased rapidly solidified alloys[J].Materials&Design,2011,32(10):4970-4979.

[2]A Levi.Heredity in cast iron[J].The iron age,6(1972)960-965.

[3] 陳光明,蔡英文,李建國,等.熔體熱處理研究及其應用[J].河北科技大學學報,1998,19(1):6.

[4] S R Ellitt.Medium-range structural ordering in covalent amorphous solids[J].Nature,1991,354:445.

[5] 邊秀房,潘學民,秦緒波,等.金屬熔體中程有序結構[J].中國科學(E 輯),2002,32(2):145.

[6]杜松昭.Al-Ni-La合金熔體結構演變規律及其對凝固組織的影響[D].濟南:濟南大學,2014.

[7] 白延文.合金熔體局域結構及其遺傳性研究[D].濟南:山東大學,2014.

[8]X F Bian,J X Zhang,Y B Jia,et al.Viscosity charateristic inmetallicmelts withmedium/short-range order structures[J].hinese Physics Letters,2005,22:644-647.

[9] Y Q Wang,Y Q W,J T Liu,et al.Discontinuous structural phase transition behaviour inmultiple component alloymelts[J].Chinese Physics Letters,2006,23:2513-2515.

[10] H S Chen,F Q Zu,J Chen,et.al.The effect ofmelt overheating on themelt structure transition and solidified structures of Sn-Bi40 alloy[J].Science in China(Series E),2008,51(9):1402-1408.

[11] 邊秀房,王偉民,李輝,等.金屬熔體結構[M].上海:上海交通大學出版社,2003.

[12] 張景祥,李輝,邊秀房.Fe3Al合金熔體分子動力學計算機模擬[J].鑄造,1999(12):22-25.

[13] S J Cheng,X F Bian,W M Wang,et al.Effect of copper,aluminum and tin addition on thermal contraction of indiummelt clusters[J].Physics B,2005,366:67-73.

[14] Han J H,Mattern N,Kaban I,et al.Phase separation in ternary Co-Gd-Tiliquids[J].JournalofPhysicsCondensedMatteranInstitute of Physics Journal,2013,25(24):245104.

[15] 孫益民,秦敬玉.過熱處理對Al-1.1%Fe合金固態組織的影響[J].山東冶金,2008,30(4):35-36.

[16] HLi,F Pederiva.Anomalies in liquid structureof Ni3Al alloys during a rapid cooling process[J].Physical Review B,2003,68:054210.

[17] J Y Qin,T K Gu,X F Bian.A study on themedium range order inmolten Fe3Si and FeSi alloys[J].J.phys:Condense Matter,2004,16:4753-4760.

[18] Y Zhao,X F Bian,X B Qin,et al.X-ray diffraction experiments on In30Sn70 from normal liquid to solidus[J].Physics Letters A,2006,356:385-391.

[19] Y Zhao,J Y Qin,X B Qin.X-ray diffraction experiments on the solidification process of Cu80Ag20alloy[J].Physics Letters A,2006,357:479-484.

[20] Y Zhao,X X Hou,X F Bian.Phase transition in liquid Cu60Sn40alloy[J].Materials letters,2008,62:3542-3544.

[21] I Kaban,S Mhiaoui,W Hoyer,et al.Surface tension and density of binary lead and lead-free Sn-based solders[J].J.Phys.:Condens.Mat.,2005,17:7867-7873.

[22] I Kaban,S Gruner,W Hoyer,et al.Effect of temperature on the structure of liquid In20Sn80[J].J.Non-crystalline Solids,2007,353:1979-1984.

[23] 桂忠樓.鎳基高溫合金BTOP工藝的發展[J].航空工程與維修,1995(4):12.

[25] 王海峰,蘇海軍,張軍,等.熔體超溫處理溫度對新型鎳基單晶高溫合金溶質分配行為的影響[J].金屬學報,2016,52(4):419-425.

[26] 殷鳳仕,管恒榮,孫曉峰,等.鑄態鎳基高溫合金M963的液態結構和熔體處理[J].金屬學報,2005,41:1190-1194.

[27] 王震,李金國,趙乃仁,等.熔體處理溫度對鎳基單晶高溫合金熔體結構和凝固組織的影響[J].金屬學報,2002,38:920-924.

[28] 閆春蕾,沈建興,馬元.塊體非晶合金制備方法的研究進展[J].山東輕工業學院學報,2009,23:8-12.

[29] X F Bian,B A Sun,Li A Hu,et al.Fragility of superheatedmelts and glass-forming ability in Al-based alloys[J].Physics Letters A,2005,335:61-67.

[30]X F Bian,J Guo,X Q Lv,et al.Prediction of glass-forming ability ofmetallic liquids[J].Applied Physics Letters,2007,91:221910.

·信息園地·

科技論文摘要寫作要求

摘要是對文章內容準確、概括而不加評論的簡短陳述。摘要應具有獨立性和完整性。一般以第三人稱的語氣,避免用“我們”、“作者”、“本文”、“本研究”等作為文摘的開頭。

1)報道性摘要。報道性摘要即資料性摘要或情報性摘要,用來報道論文所反映的作者的主要研究成果,向讀者提供論文中全部創新內容和盡可能多的定量或定性的信息。尤其適用于試驗研究和專題研究類論文。篇幅100~300字。其內容比較詳細,應包括研究的目的、方法、結果和結論,應是具體內容,盡量反映文章的主要信息。重點是結果和結論,將論文創新的、具有特色的內容表達出來。一般情況下,向學術性期刊投稿應選用報道性摘要形式。

2)指示性摘要。僅僅指出文章用什么方法研究了什么問題,而不涉及結果及結論。以數學解析為主的論文、綜述性文章常用這種摘要。篇幅一般50~100字。

3)報道-指示性摘要。以報道性文摘的形式表述論文中最重要的內容,即信息價值較高的部分,而以指示性文摘的形式表述其余作者認為不太重要的內容,以避免摘要寫得太長且重點不突出。

(燕明宇)

Research Progress on the Relationship between Liquid and Solid Metal Material Structure

SUN Yongjun
(Shandong Construction Engineering Co.,Ltd.,Jinan Iron and Steel Group,Jinan 250101,China)

This paper summarizes the researchmethods,research results and application of liquid structure at home and abroad in recent decades.Some information about high temperature liquid structure has been obtained by high temperature X-ray diffraction,high temperature viscosity and computer simulation,which provides theoretical guidance for the transformation of liquid-solid structure and the acquisition of alloys with excellent properties.In order to obtainmore information on themelt structure,it is suggested to develop a high-precision equipment to test the liquid structure,and to find the potential function which is as accurate and simple as possible to simulate the liquid structure.

metalmaterial;liquid structure;solid structure;correlation;research progress

TG111.4

A

1004-4620(2017)04-0001-04

2016-08-05;

2017-07-12

孫永軍,男,1982年生,2005年畢業于山東建筑工程學院金屬材料工程專業。現為濟鋼集團山東建設工程有限公司工程師,研究方向為焊接和鑄造。

猜你喜歡
結構研究
FMS與YBT相關性的實證研究
《形而上學》△卷的結構和位置
哲學評論(2021年2期)2021-08-22 01:53:34
2020年國內翻譯研究述評
遼代千人邑研究述論
視錯覺在平面設計中的應用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
論結構
中華詩詞(2019年7期)2019-11-25 01:43:04
新型平衡塊結構的應用
模具制造(2019年3期)2019-06-06 02:10:54
EMA伺服控制系統研究
新版C-NCAP側面碰撞假人損傷研究
論《日出》的結構
主站蜘蛛池模板: m男亚洲一区中文字幕| 99九九成人免费视频精品| 色综合中文字幕| 久爱午夜精品免费视频| 亚洲性视频网站| 五月天香蕉视频国产亚| 日韩人妻无码制服丝袜视频| 日韩视频免费| 9999在线视频| 欧美一区二区精品久久久| 亚洲精品视频网| 国产欧美成人不卡视频| 亚洲一级毛片在线播放| 97视频在线精品国自产拍| 网友自拍视频精品区| 亚洲人成色在线观看| 午夜视频在线观看区二区| 亚洲性网站| 国产在线精品美女观看| 呦女亚洲一区精品| 日韩欧美中文在线| 精品夜恋影院亚洲欧洲| 在线精品视频成人网| 全免费a级毛片免费看不卡| 99热国产这里只有精品无卡顿" | 亚洲视频色图| 欧美性色综合网| 98超碰在线观看| 国产三级毛片| 国产高清无码麻豆精品| 国产福利小视频高清在线观看| 区国产精品搜索视频| 国产美女人喷水在线观看| 91精品国产丝袜| 人妻精品全国免费视频| 午夜日b视频| 97se综合| 高清无码手机在线观看| 91丝袜乱伦| 狼友视频一区二区三区| 天天色天天操综合网| 91高清在线视频| 精品91自产拍在线| 日日碰狠狠添天天爽| 日韩精品无码免费一区二区三区 | 亚洲成在线观看| 日韩经典精品无码一区二区| 99精品影院| 久久女人网| 狠狠色狠狠色综合久久第一次| 国产精品亚洲专区一区| 免费xxxxx在线观看网站| 国内精品视频在线| 国产一区在线观看无码| m男亚洲一区中文字幕| 在线观看国产精品第一区免费| 久久精品国产91久久综合麻豆自制| 51国产偷自视频区视频手机观看| 中国特黄美女一级视频| 国产一级在线观看www色| 在线观看国产小视频| 亚洲天堂.com| 亚洲丝袜中文字幕| 亚洲精选高清无码| 欧美激情第一区| 午夜毛片免费观看视频 | 91免费观看视频| 色有码无码视频| 国模沟沟一区二区三区| 无码高潮喷水专区久久| 久草视频精品| 免费av一区二区三区在线| 日本尹人综合香蕉在线观看 | 色婷婷狠狠干| 精品五夜婷香蕉国产线看观看| 欧美性天天| 最新精品国偷自产在线| 久久网欧美| 国产精品白浆在线播放| 国产新AV天堂| 色综合激情网| h网址在线观看|