高辣慎入hh出轨,男男黄gay片免费网站www,亚洲国产成人久久精品大牛影视,少妇性活bbbbbbbbb四川

預(yù)存
Document
當(dāng)前位置:文庫百科 ? 文章詳情
同步輻射全散射PDF從基礎(chǔ)到實例分析
來源:測試GO 時間:2022-08-09 09:14:32 瀏覽:9334次

一、PDF以及實驗基礎(chǔ)理論

1.1什么是PDF和RDF

眾所周知,結(jié)構(gòu)的測試方法有多種,如核磁共振(NMR)、穆斯堡爾譜學(xué)和衍射法,其中的NMR和穆斯堡爾譜學(xué)是用于檢測局域有序的核物理方法。然而,作為更大范圍使用的衍射法(包括有XRD、中子衍射、電子衍射以及同步輻射等),其關(guān)鍵的分析技術(shù)就是原子的分布函數(shù)法。

分布函數(shù)是一個統(tǒng)計物理學(xué)的概念,表示在一個原子附近的徑向原子分布狀態(tài),一般有三種形式:徑向分布函數(shù)(Radial Distribution Function,RDF)、雙體分布函數(shù)g(r)(Pair Distribution Function,PDF又稱對分布函數(shù))以及約化徑向分布函數(shù)G(r)。

圖1.1 原子結(jié)構(gòu)與徑向密度函數(shù)的關(guān)系示意圖


我們使用最多的是對分布函數(shù)PDF,對于N元體系來說,就有N(N-1)個偏雙體分布函數(shù): 

對于PDF的函數(shù)圖的分析,可以得到非晶合金結(jié)構(gòu)的重要參數(shù)[1, 2]:如圖1.2所示兩種不同的材料狀態(tài)的對分布函數(shù)。

A 最近鄰原子平均距離:雙體分布函數(shù)第一峰位值即為最近鄰原子平均值;

B 最近鄰配位數(shù):CN等于PDF函數(shù)的第一個半峰面積的兩倍,即:

C 原子的平均位移:位移等于RDF函數(shù)一個半峰寬度的1/2.36倍;

D 有序疇尺寸:在雙體分布函數(shù)約等于1,即為有序疇邊界,考慮實驗誤差一般取g(r)=1.02時的r對應(yīng)有序疇尺寸;

E 第一峰的高度和形狀可以反應(yīng)原子之間的結(jié)合強度,如峰表現(xiàn)為尖銳,表明在此半徑范圍內(nèi)的原子數(shù)密度比平均密度要高很多,中心原子與最近鄰原子的相互結(jié)合強度也比較大;

F 第一峰與第一谷的關(guān)系通過經(jīng)驗性的方法,可用來確定玻璃轉(zhuǎn)變點。


圖1.2 對分布函數(shù)



1.2靜態(tài)結(jié)構(gòu)因子

靜態(tài)結(jié)構(gòu)因子(Static Structure Factor,SSF) 是一個連接實驗分析與模擬分析的重要參量。對于衍射分析來說,它表征的是材料對射線的散射能力,反映結(jié)構(gòu)的平均信息。結(jié)構(gòu)因子S(Q)可以和上述的PDF相互轉(zhuǎn)化(互為傅里葉變換):

同理,如為中子散射,則將原子散射因子換成中子散射長度;

若為電子散射,則替換成相應(yīng)的散射因子。


1.3衍射實驗

X射線衍射技術(shù)(X-ray Diffraction, XRD)是材料分析的最常用手段之一。對于晶體材料,XRD測量布拉格衍射峰,提供樣品的物相和晶體結(jié)構(gòu)信息。

X射線衍射是屬于散射的一個特例了,是相干散射,強調(diào)入射波與目標(biāo)物質(zhì)發(fā)生相互作用后傳播方向變化的過程,且衍射偏向強調(diào)波被周期性結(jié)構(gòu)散射后的結(jié)果。對納米材料、半晶體、非晶、液體等材料,衍射的相干性降低,同樣的進(jìn)行XRD過程,我們稱為X射線散射。

所以,散射衍射本質(zhì)上是同樣的東西作用于不同物體上的不同稱呼,沒有本質(zhì)的差別。從散射的空間分布的剖面看,不強調(diào)某個特定方向的時候泛稱散射,而要強調(diào)某個特定的模式、方向的分布時則說衍射。

彈性散射是指在散射中沒有頻率位移、無能量損失,也就是相干散射。而非彈性散射,也就是說散射后不但能量改變,方向也有變化,如Compton散射。我們在獲得PDF的時候需要的就是相干散射數(shù)據(jù),所以稱為X射線全散射實驗(X-ray total scattering),其操作方式、原理和XRD相同。

二、實驗到PDF

2.1PDF與全散射

我們利用實驗室XRD可以得到衍射圖譜,進(jìn)而得到晶體材料的結(jié)構(gòu)信息。但是對于非晶等材料,常規(guī)XRD分析僅能獲得寬化的衍射峰或背景包,無法得到有意義的結(jié)構(gòu)信息,如圖2.1中右所示。

圖2.1 晶體(左),與非晶合金(右)XRD圖譜

根據(jù)上章節(jié)所說的原子對分布函數(shù)(PDF, Pair distribution function),描述了在所研究材料中發(fā)現(xiàn)距離為r的一對原子的概率。此方法以高能硬X射線測量樣品廣角全散射數(shù)據(jù),因此也稱為Total scattering全散射分析,同時對布拉格衍射峰和漫散射進(jìn)行歸一化和傅里葉變換等處理,不僅提供長程(>10 nm)原子有序性信息,還提供材料中短程結(jié)構(gòu)信息,如短程有序/無序排布、鍵長、局部缺陷等。

通過對不同狀態(tài)的同類樣品的 PDF 數(shù)據(jù)進(jìn)行差異化分析,還可以進(jìn)一步研究過程中材料精細(xì)結(jié)構(gòu)的變化,獲得材料物理性能或化學(xué)性能的變化與材料結(jié)構(gòu)變化之間的關(guān)系,深入研究變化/反應(yīng)過程機理。PDF極大拓展了X射線結(jié)構(gòu)表征的分析范圍,樣品不再局限于晶態(tài)材料,非晶、液體等均可測量。

2.2如何獲得的PDF

獲得PDF數(shù)據(jù)有三個步驟,如下圖2.2所示:

1. 對樣品的全散射數(shù)據(jù)進(jìn)行校正,包括背景(對空白樣品容器進(jìn)行同樣條件的單獨測量獲得)、康普頓散射、探測器死區(qū)時間、吸收、衍射幾何和偏振;

2. 對校正后的X射線散射數(shù)據(jù)進(jìn)行歸一化,計算簡化結(jié)構(gòu)函數(shù)S(Q);

3. 對結(jié)構(gòu)函數(shù)進(jìn)行傅里葉變換,得到原子對分布函數(shù)g(r)。

 

圖2.2 PDF計算過程

注意事項:

A 在結(jié)構(gòu)函數(shù)S(Q)中,低Q部分的數(shù)據(jù)主要為布拉格散射的尖銳強峰,它們對應(yīng)于材料中原子的長程排序;在高Q值部分,布拉格散射峰基本消失,強度極低的漫散射寬峰才可以識別出來[4]。因此,數(shù)據(jù)的Q值范圍足夠大,才能獲得盡可能多的漫散射信號。

B 由于Q=4πsinθ/λ,為獲得盡可能大的Q值(Qmax)的原始數(shù)據(jù),PDF實驗要求使用盡可能短波長的X射線,并且測量角度范圍(2θ)應(yīng)盡可能大。這樣在最終計算出PDF函數(shù)時,實空間的分辨率也會更好。此外,由于漫散射信號強度比布拉格散射峰信號強度低5-6個數(shù)量級,在測試時一般需要對全譜數(shù)據(jù)進(jìn)行分段測量,在高角度(高Q部分)使用足夠長的測量時間來獲得足夠強度和信噪比的清晰、平滑的高質(zhì)量原始數(shù)據(jù)。

總結(jié):實驗需要短波長和高強度的光源,且在實驗時需要保證光源波長遠(yuǎn)小于或稍微長于測試元素的吸收限。在實際工作中,同步輻射光源和加速器線站天然具有高強度多波長的射線源,因此經(jīng)常在粉末衍射線站搭建PDF光路,使用單色器選取短波長高能射線進(jìn)行PDF實驗,歐洲光源(ESRF)等還有專用的高能PDF 線站(ID11)用于實現(xiàn)更高Qmax和測量效率的PDF數(shù)據(jù)。

對于整個測試數(shù)據(jù)的歸一化處理數(shù)據(jù)可以參考黃勝濤老師的《非晶態(tài)材料的結(jié)構(gòu)和分析》[5]。目前多采用的軟件為PDFgetX2,或者PDFgetX3[6]。


2.3實驗操作過程演示 

同步輻射X射線衍射(XRD)實驗在澳洲的ANSTO光束線上進(jìn)行。由于實驗設(shè)備的限制選取了17KeV能量的X射線,以保證不會出現(xiàn)重金屬的強烈吸收導(dǎo)致的熒光,從而進(jìn)一步對背底進(jìn)行干擾。在相同的實驗條件下,對裝材料的毛細(xì)管進(jìn)行測試,獲得背景數(shù)據(jù)。利用jade軟件去除毛刺。然后,扣除背底、樣品吸收,極化以及集合系數(shù)校正,并對非相干散射(Compton散射)和多重散射進(jìn)行消除,得到最后的相干散射強度Icoh(Q)。

圖2.3 背景擬合曲線

第一步,如圖2.3,對背景進(jìn)行擬合,由于同步輻射伴隨加壓過程,所以在前端低角度區(qū)域的強度比較高。

第二步,將去除毛刺的數(shù)據(jù)和背景數(shù)據(jù)帶入PDFgetx3 PDFgetX3, PDFgetN3 and PDFgetS3 — DiffPy documentation(網(wǎng)址https://www.diffpy.org/products/pdfgetx.html#pdfgetx3-and-pdfgetn3)中進(jìn)行歸一化處理,這里PDFgetX3相對于上一代程序更為方便的有了python的交互操作,操作更加的簡潔。PDFgetX2需要虛擬機器才能運行,但是有獨立的界面。

利用PDFgetX3的處理方法有多種,詳情可以參考PDFgetX3的manul,這里展示的是一種簡單具體的方法。

安裝好軟件后,打開shell,輸入:

pdfgetx3 –-createconfig=“filename.cfg” 創(chuàng)建一個全新的輸入文件

打開文件后對內(nèi)部參數(shù)進(jìn)行修改調(diào)整:

Dataformat= "twotheta", "QA", "Qnm"

輸入的強度數(shù)據(jù)的橫坐標(biāo)單位

inputfile = xrayfile01.chi

輸入文件的名稱

backgroundfile =“filename”

背景文件名(若沒有可以去掉這一項)

bgscale = 1

背景強度的調(diào)整,默認(rèn)為1

output = 輸出的文件名字

outputtype = fq, gr,iq,sq

輸出的數(shù)據(jù)的類型

force = yes

建議把強制打開,可以強行執(zhí)行,覆蓋前面的數(shù)據(jù)

mode = xray

選擇實驗的類型,中子散射或者X射線散射

wavelength = 0.7288

實驗波長,單位是埃,需要在輸入文件的橫坐標(biāo)為角度時使用

composition = Cu 56 Zr44

材料成分

rpoly = 0.63

F(Q)校正多項式中最大頻率的r限值(這個后續(xù)還可以調(diào)整)

qmin = 0.8    qmax = 20  

傅里葉變化的范圍 單位埃

qmaxinst = 20

輸入文件的有效值范圍 單位埃

rmin = 0.8

rmax = 30.0

rstep = 0.05     

計算出的PDF r-grid的界限和間距

plot = fq, gr,iq,sq

計算完成后的繪圖

interact = yes

使用Ipython交互    

所有的數(shù)據(jù)都修改完整后運行此文件

pdfgetx3 -c cu56.cfg  

出現(xiàn)如下圖所示的處理結(jié)果,并繪圖得到了pdf數(shù)據(jù)

在得到數(shù)據(jù)后,如果對數(shù)據(jù)不滿意可以在Ipython交互界面輸入tuneconfig(),然后進(jìn)行下一步的修改微調(diào),最后得到對分布函數(shù)。值得一說的是,若數(shù)據(jù)來自于中子衍射,可以在.cfg文件中改變-mode這個選項或者將運行命令改變成 pdfgetn3 ……


三、深入PDF

根據(jù)前面兩章的學(xué)習(xí),我們獲得的pdf數(shù)據(jù)是總體的結(jié)構(gòu)因子和徑向分布函數(shù)等。而且,結(jié)構(gòu)因子并不是唯一的,這一點要知道。在利用中子散射以及X射線的散射兩種方法做出來的結(jié)構(gòu)因子是有所差別的,如圖3.1所示:

 

圖3.1 XRD和ND結(jié)構(gòu)因子[7, 8]

兩者的趨勢相同,但是細(xì)節(jié)上有所不同。為何會有這種差別呢?這是由于一種材料結(jié)構(gòu)確定時,其對分布函數(shù)也確定了,傅里葉變換后會有不同的偏結(jié)構(gòu)因子,再利用系數(shù)進(jìn)行組合的時候由于系數(shù)的不同所以組合出來的總結(jié)構(gòu)因子會不同。XRD利用的是散射因子計算系數(shù),而中子散射則利用的是中子散射長度。所以同樣的物品有不同的結(jié)構(gòu)因子,當(dāng)然如果材料為單一元素將不會存在這個問題。

通常我們獲取PDF后,主要是了解其具體的原子分布規(guī)律,而總的PDF是不能滿足我們的需求的,所以需要換算偏函數(shù)。當(dāng)體系存在N種元素時,就有N(N-1)種的組合,也就是這么多種類的偏函數(shù)。但我們只知道總結(jié)構(gòu)因子,所以一個方程解三個未知數(shù)是不可能的,需要三個不同的實驗才能解出偏結(jié)構(gòu)函數(shù)。當(dāng)然,我們還可以利用異常X射線散射得到偏函數(shù)。

以前在我們得到了總的SQ或者gr后,可能直接利用模擬數(shù)據(jù)與之比較,確定其真實性。但隨著模擬技術(shù)的發(fā)展,出現(xiàn)了逆蒙特卡洛算法,可以利用SQ數(shù)據(jù)建模得到三維結(jié)構(gòu),然后求得徑向分布函數(shù)等。然而對于RMC的計算需要利用到結(jié)構(gòu)因子以及EXAFS的實驗數(shù)據(jù),雙向擬合才讓數(shù)據(jù)更具有真實性。

RMC的基礎(chǔ)是建立在MC[9]方法上的,MC方法認(rèn)為求解一個問題是一個概率模型隨機過程,通過對模型的反復(fù)抽樣,對結(jié)果不斷地統(tǒng)計分析然后反饋誤差,最后得到所求的參量的近似解。而RMC是一種新的改進(jìn),利用RMC求解結(jié)構(gòu)有以下幾個步驟[10]

所以在對EXAFS模擬時,需要知道單個背散射原子的貢獻(xiàn),以便于模擬。求解背散射因子利用軟件feff[11]獲得。近年來,已有很多學(xué)者通過RMC方法建模,分析非晶合金的結(jié)構(gòu)特征,取得很多不錯的成績[12-20]。

雖然已有學(xué)者通過MD和RMC的模型比對,證實了RMC結(jié)果的合理性[17],然而逆蒙特卡羅是一種抽樣的數(shù)學(xué)統(tǒng)計方法,本身是沒有物理意義的,僅僅是在原子堆垛的角度來考慮模擬的正確性,而且是通過一維的數(shù)據(jù)進(jìn)行三維構(gòu)型的模擬,也就是說可能存在不同的構(gòu)型,但是其一維數(shù)據(jù)相同,針對這一問題我們就需要給定足夠多的輸入條件來提高模型的可行度。

實例分析:

利用SQ數(shù)據(jù)和EXAFS數(shù)據(jù)帶入RMC程序中進(jìn)行反向擬合,最后可以得到偏對分布函數(shù)。 

  

 目前PDF分析表征多用于非晶合金領(lǐng)域,尋找無序結(jié)構(gòu)的規(guī)律。已有很多的學(xué)者通過XRD等衍射實驗進(jìn)行PDF分析得到真實理想的三位原子結(jié)構(gòu)[21-25]。


參考文獻(xiàn)

[1] 徐勇, 于美杰. 非晶態(tài)合金原子結(jié)構(gòu)及結(jié)晶動力學(xué) [M]. 2015.

[2] 張世良, 戚力, 高偉, et al. 分子模擬中常用的結(jié)構(gòu)分析與表征方法綜述 [J]. 燕山大學(xué)學(xué)報, 2015, (03): 28-35.

[3] Ziman, M. J. A theory of the electrical properties of liquid metals. I: The monovalent metals [J]. 1961, 6(68): 1013-34.

[4] Egami T. Underneath the Bragg peaks. Structural analysis of complex materials [M]. Underneath the bragg peaks : Structural analysis of complex materials, 2003.

[5] 黃勝濤. 非晶態(tài)材料的結(jié)構(gòu)和結(jié)構(gòu)分析 [M]. 北京:科學(xué)出版社, 1987.

[6] Juhás P, Davis T, Farrow C L, et al. PDFgetX3: A rapid and highly automatable program for processing powder diffraction data into total scattering pair distribution functions [J]. Journal of Applied Crystallography, 2013, 46(2): 560-6.

[7] Kaban I, Jóvári P, Escher B, et al. Atomic structure and formation of CuZrAl bulk metallic glasses and composites [J]. Acta Materialia, 2015, 100: 369-76.

[8] Mattern N, Jóvári P, Kaban I, et al. Short-range order of Cu–Zr metallic glasses [J]. Journal of Alloys and Compounds, 2009, 485(1-2): 163-9.

[9] Metropolis N, Ulam S. The Monte Carlo Method [J]. Journal of the American Statistical Association, 1949, 44(247): 335-41.

[10] 孫永麗. Cu-Zr基非晶合金結(jié)構(gòu)及動力學(xué)數(shù)值模擬 [D]; 哈爾濱工業(yè)大學(xué), 2010.

[11] Rehr J J, Ankudinov A L, Conradson S D. Real-space multiple-scattering calculation and interpretation of x-ray-absorption near-edge structure [J]. Physrevb, 1998, 58(12): 7565-76.

[12] Yu Q, Wang X D, Lou H B, et al. Atomic packing in Fe-based metallic glasses [J]. Acta Materialia, 2016, 102: 116-24.

[13] Li K, Chou Y-J, Gao F-L, et al. Atomic Structure of Cu49Hf42Al9 Metallic Glass with High Glass-Forming Ability and Plasticity [J]. Acta Metallurgica Sinica (English Letters), 2019, 32(7): 803-7.

[14] Babilas R, Dariusz ?, Laszlo T. Atomic structure of Mg-based metallic glass investigated with neutron diffraction, reverse Monte Carlo modeling and electron microscopy [J]. Beilstein Journal of Nanotechnology, 2017, 8(1): 1174-82.

[15] Khanna A, Fábián M, Hirdesh, et al. Structural analysis of WO3-TeO2 glasses by neutron, high energy X-ray diffraction, reverse Monte Carlo simulations and XANES [J]. Journal of Non-Crystalline Solids, 2018, 495: 27-34.

[16] Wetherall K M, Pickup D M, Newport R J, et al. The structure of calcium metaphosphate glass obtained from x-ray and neutron diffraction and reverse Monte Carlo modelling [J]. J Phys Condens Matter, 2009, 21(3): 035109.

[17] Fang X W, Huang L, Wang C Z, et al. Structure of Cu64.5Zr35.5 metallic glass by reverse Monte Carlo simulations [J]. Journal of Applied Physics, 2014, 115(5).

[18] Kaur R, Khanna A, Hirdesh, et al. Structure of strontium tellurite glass, anti-glass and crystalline phases by high-energy X-ray diffraction, reverse Monte Carlo and Rietveld analysis [J]. Acta Crystallographica Section B Structural Science, Crystal Engineering and Materials, 2020, 76(1): 108-21.

[19] 杜金成, 張久明, 趙修建. 硫系玻璃結(jié)構(gòu)的RMC計算機模擬 [J]. 材料研究學(xué)報, 1998, (2): 191-4.

[20] 張久明, 杜金成, 趙修建, et al. Ge-As-Te玻璃結(jié)構(gòu)的X射線衍射和RMC模擬研究 [J]. 無機材料學(xué)報, 1999, 14(3).

[21] Itoh K, Yamada R, Saida J, et al. Atomic-level characterization of free volume in the structure of Cu67Zr33amorphous alloy [J]. J Phys Condens Matter, 2021, 33(27).

[22] Guo G-Q, Wu S-Y, Luo S, et al. Detecting Structural Features in Metallic Glass via Synchrotron Radiation Experiments Combined with Simulations [J]. Metals, 2015, 5(4): 2093-108.

[23] Ellersdorfer P, Petersen T C, Opletal G, et al. Extracting nanoscale structures from experimental and synthetic data with reverse Monte Carlo [J]. Nano Futures, 2021, 5(2).

[24] Wajhal S, Shinde A B, Krishna P S R. Reverse Monte Carlo modeling of LiNbO3-TeO2 glasses [Z]. Dae Solid State Physics Symposium 2019. 2020.10.1063/5.0016729

[25] Jóvári P, Kaban I, Escher B, et al. Structure of glassy Cu47.5Zr47.5Ag5 investigated with neutron diffraction with isotopic substitution, X-ray diffraction, EXAFS and reverse Monte Carlo simulation [J]. Journal of Non-Crystalline Solids, 2017, 459: 99-102.


測試狗

評論 / 文明上網(wǎng)理性發(fā)言
12條評論
全部評論 / 我的評論
最熱 /  最新
全部 3小時前 四川
文字是人類用符號記錄表達(dá)信息以傳之久遠(yuǎn)的方式和工具?,F(xiàn)代文字大多是記錄語言的工具。人類往往先有口頭的語言后產(chǎn)生書面文字,很多小語種,有語言但沒有文字。文字的不同體現(xiàn)了國家和民族的書面表達(dá)的方式和思維不同。文字使人類進(jìn)入有歷史記錄的文明社會。
點贊12
回復(fù)
全部
查看更多評論
相關(guān)文章

一文詳解掃描電子顯微鏡(SEM)的工作原理及應(yīng)用技術(shù)

2023-10-08

熱重分析(TG-DTG)曲線的幾種解析方法

2023-12-26

接觸角測試(CA)的原理、樣品制備要求及實際應(yīng)用

2023-11-16

一文詳細(xì)介紹he染色的基本原理、實驗步驟及注意事項

2023-11-23

恒電流間歇滴定法GITT的基本原理以及測試教程

2022-08-12

TMA技術(shù)的基本概念、應(yīng)用領(lǐng)域以及實際操作中的要點

2023-12-06

熱門文章/popular

基礎(chǔ)理論丨一文了解XPS(概念、定性定量分析、分析方法、譜線結(jié)構(gòu))

晶體結(jié)構(gòu)可視化軟件 VESTA使用教程(下篇)

手把手教你用ChemDraw 畫化學(xué)結(jié)構(gòu)式:基礎(chǔ)篇

【科研干貨】電化學(xué)表征:循環(huán)伏安法詳解(上)

電化學(xué)實驗基礎(chǔ)之電化學(xué)工作站篇 (二)三電極和兩電極體系的搭建 和測試

【科研干貨】電化學(xué)表征:循環(huán)伏安法詳解(下)

微信掃碼分享文章
国产亚洲av无码av男人的天堂| 国产精品久久久久精品日日| 粉嫩被黑人两根粗大猛烈进出视频| 美味人妻2016| 成人电影在线观看| 99re6热在线精品视频播放| 2023国精产品一二二线精华液| 律师事务所在线咨询免费| 国产精品三级在线观看无码| 少妇人妻互换不带套| 天天干天天日| 亚洲精品久久久久中文字幕二区| 人妻熟女一二三区夜夜爱| 国产成人无码综合亚洲日韩| 久久变态刺激另类sm按摩| 玩弄白嫩少妇xxxxx性| 大地资源电影中文在线观看| 蜜桃传媒一区二区亚洲av婷婷| 日本韩国男男作爱gaywww| 国产丰满老熟妇乱xxx1区| 国产aⅴ激情无码久久久无码| 艳肉观音性三级dvd| 国产精品素颜系列| 亚洲熟妇AV日韩熟妇在线| 国色天香久久久久久久小说| 日本无码小泬粉嫩有套在线| 欧美精品久久久久久久自慰| 晚上想看点小片片| 五十六十丰满老熟妇hd| 老妇女xxxxx性黑寡妇小说| 被几个人摁住玩弄隐私作文| 亚洲精品无码久久久久秋霞| 国精产品一品二品国精破解| 久久精品国产亚洲av高清热| 小洁和公h文翁17| 高h各种姿势调教1v1| 妺妺窝人体色www在线小说| 美国式禁忌2| 最近2019年中文字幕无吗| 舌头伸进去添的我爽高潮| 大又大粗又爽又黄少妇毛片|