亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        亞穩(wěn)β 鈦合金的變形機制的研究進展

        2021-04-10 09:08:52吳松全黃愛軍
        有色金屬材料與工程 2021年1期
        關鍵詞:孿晶馬氏體鈦合金

        徐 迪,楊 義,吳松全,黃愛軍

        (上海理工大學 材料科學與工程學院,上海 200093)

        亞穩(wěn)β 鈦合金具有高強度、低彈性模量和優(yōu)異的生物相容性,因此被廣泛應用于航空航天工業(yè)和生物醫(yī)療器械等多個領域[1-6],吸引了世界范圍內(nèi)眾多研究人員的注意。

        亞穩(wěn)β 鈦合金的變形機制與β 鈦合金的有所不同。β 鈦合金自β 相區(qū)淬火到室溫的相組成為穩(wěn)定的β 相,其在室溫變形時只發(fā)生位錯滑移,不會產(chǎn)生應力誘發(fā)相變,只有在加熱或者低溫變形時才會發(fā)生相變。而亞穩(wěn)β 鈦合金自β 相區(qū)淬火到室溫的相組成為不穩(wěn)定的β 相,對其進行室溫變形時,除發(fā)生位錯滑移外,還可能發(fā)生應力誘發(fā)α"馬氏體相變(SIMα")、應力誘發(fā)ω 相變(SIω),以及{112}<111>孿生、{332}<113>孿生等[7-9]。這些變形機制通常受到β 相穩(wěn)定性、變形溫度、變形速率和變形量等因素的影響,導致亞穩(wěn)β 鈦合金的變形機制十分復雜。

        本文介紹了亞穩(wěn)β 鈦合金的多種變形機制,以及影響變形機制的因素,并提出了該領域可能的研究方向。

        1 亞穩(wěn)β 鈦合金常見的變形機制

        1.1 滑移

        滑移在金屬的塑性變形過程中是一種主要的變形方式。體心立方結(jié)構(gòu)是一種非密排結(jié)構(gòu),因此體心立方金屬的滑移面不太穩(wěn)定,通常在低溫時為{112},在中溫時為{110},在高溫時為{123},不過它的滑移方向很穩(wěn)定,總是<111>?;谱冃卧趤喎€(wěn)β 型鈦合金中是一個普遍現(xiàn)象,無論相穩(wěn)定性系數(shù)的大小,在變形時或多或少都有1/2<111>位錯滑移的參與。

        亞穩(wěn)β 鈦合金體系在變形中存在多種變形機制,其中位錯滑移會在塑性變形階段與其他的變形產(chǎn)物產(chǎn)生相互作用。由于細晶強化是位錯與界面的交互作用造成的,因此位錯與新形成的變形產(chǎn)物的交互作用會形成動態(tài)細晶強化效果。在變形量較大的樣品中,位錯滑移與孿生或應力誘發(fā)相變之間產(chǎn)生了某種競爭機制,主導變形產(chǎn)物會隨著變形量的增大而變化[10]。

        也有研究認為,除了競爭機制,在應變量較小的范圍內(nèi),位錯的產(chǎn)生會通過降低孿晶形核能,增加孿晶的形核位置,從而促進孿晶的形核[11],因此在一定程度上位錯滑移會增加變形機制的多樣性。

        1.2 孿生

        孿生是金屬塑性變形過程中另一種重要的變形機制,通常認為孿生會在滑移不易進行時發(fā)生。在亞穩(wěn)β 鈦合金中最常見的孿晶為{332}<113>孿晶和{112}<111>孿晶,它們與基體的晶體學取向差為以<110>為軸分別旋轉(zhuǎn)50.5°和70.3°,這兩種孿晶多見于室溫變形后的合金中,前者的尺寸比后者要大得多。在傳統(tǒng)的體心立方晶體中,最常見的孿晶是{112}<111>孿晶,其孿晶面對應著晶體的密排面[12]。{332}<113>孿晶是亞穩(wěn)β 鈦合金中一種特殊的孿晶,這種高指數(shù)孿晶最早是由Blackburn 等[13]于1971 年在Ti-11.5Mo-6.0Zr-4.5Sn 合金中發(fā)現(xiàn),其孿晶面與切變方向和密排面及密排方向并不對應。理論上孿晶很難生成,有研究認為這與β 相的不穩(wěn)定性有關[12]。圖1 為變形后亞穩(wěn)β 鈦合金的電子背向散射衍射(electron back-scattered diffraction,EBSD)分析,其中的透鏡狀層片狀結(jié)構(gòu)為{332}<113>孿晶[14]。

        圖1 變形亞穩(wěn)β 鈦合金的電子背向散射衍射分析圖[14]Fig.1 EBSD analysis images of the deformed metastable β titanium alloy[14]

        {332}<113>孿晶常見于二元Ti-Nb、Ti-V、Ti-Mo、Ti-Cr、Ti-Fe 合金和Ti-Nb 基、Ti-V 基、Ti-Ta 基的合金體系中。由于{332}<113>孿晶特殊的晶體學特征,眾多學者對其形成機制提出了假說,比如Takemoto等[15]提出的β ?α"馬氏體機制,認為β 相基體中先產(chǎn)生α"相,α"相再逆轉(zhuǎn)變?yōu)閷\晶。Kawabata 等[16]提出位錯機制,認為此孿晶來源于1/2[11] →1/11[55]+1/22[11]的位錯分解和不全位錯的滑移。Tobe等[12]提出晶格不穩(wěn)定機制,認為通過晶格調(diào)制{332}<113>孿晶在亞穩(wěn)β 鈦合金中比{112}<111>孿晶更容易生成。這些理論都能在一定程度上解釋此種孿晶的形成機制,但都缺少充足的實驗支持[8]。

        變形過程中的孿生行為會影響合金的力學性能,根據(jù)動態(tài)Hall-Petch 方程,由于產(chǎn)生孿晶而產(chǎn)生的晶粒細化效應可以提高材料的強度,當孿晶的密度變大時,由于晶界增多而導致位錯堆積,從而導致更高的應變硬化率,使得材料的強度和均勻伸長率同時提高,形成孿生增塑效應(twinning induced plasticity,TWIP)。Ren 等[17]和Zhang 等[18]通過在拉伸變形過程中誘導出納米級孿晶,先后研發(fā)出了抗拉強度達1.1 GPa 并具有高伸長率(> 25%)的亞穩(wěn)β 鈦合金。

        1.3 應力誘發(fā)α"馬氏體相變

        亞穩(wěn)β 型鈦合金的α"馬氏體轉(zhuǎn)變開始溫度(Ms)低于室溫,因此從β 相區(qū)淬火至室溫時不會發(fā)生α"馬氏體轉(zhuǎn)變。但是如對合金進行塑性變形,α"馬氏體轉(zhuǎn)變將會發(fā)生,此時形成的馬氏體稱為應力誘發(fā)α"馬氏體(SIMα")。

        α"馬氏體屬于斜方結(jié)構(gòu)。有理論認為[19]應力誘發(fā)α"馬氏體相變是通過一個不變平面應變(invariant plane strain, IPS)形成的。在塑性變形的過程中,由于某些晶粒中原子的位置和距離被改變,晶體結(jié)構(gòu)也會變化,當β 相轉(zhuǎn)變?yōu)棣?相的時候,β 基體相中原子發(fā)生了協(xié)同運動,由于其慣習面是不變的,這個面作為β 基體相與α"馬氏體之間的界面不會產(chǎn)生旋轉(zhuǎn)和畸變,因此產(chǎn)生了與β 基體相高共格的α"馬氏體。由于β 相與α"相之間存在[001]β∥[100]α"和(110)β∥(001)α"的 取向 關 系,使由形變產(chǎn)生相變成為可能。

        應力誘發(fā)α"馬氏體相變對鈦合金力學性能的影響有著大量研究[1,6,19-24]。最近的研究表明,可以通過調(diào)控SIMα"的產(chǎn)生來形成相變誘發(fā)塑性效應(transformation induced plasticity,TRIP),從而改善亞穩(wěn)β 鈦合金的力學性能[25],控制SIMα"為優(yōu)化鈦合金的力學性能提供了一條新途徑。在早期的鈦合金相變與塑性的研究中,如對Ti-10V-2Fe-3Al[26],β-Cez[20]和Ti-8.0Mo-3.9Nb-2.0V-3.1Al[27]等合金的研究,發(fā)現(xiàn)通過控制亞穩(wěn)態(tài)β 相中α"馬氏體的含量可以平衡材料的強度和延展性。但通常情況下,α"馬氏體轉(zhuǎn)變的發(fā)生會導致材料的屈服強度下降,不利于材料的綜合性能[23,28]。

        1.4 應力誘發(fā)ω 相變

        ω 相是一種六方結(jié)構(gòu)的相,其結(jié)構(gòu)如圖2 所示。ω 相自被發(fā)現(xiàn)以來一直受到廣泛的關注,對于其形成機制、力學和物理性能,眾多學者都有系統(tǒng)的研究[3,29-36]。在亞穩(wěn)β 鈦合金體系中,ω 相會有多種形成方式,比如通過熱處理,或通過外加應變誘發(fā)。其中熱處理方式又可分為從β 單相區(qū)快速冷卻(淬火)和對材料進行等溫時效,而對于外加應變誘發(fā)的方式,ω 相則是通過高應變速率壓縮(沖擊載荷),以機械方式形成的[37],其形貌呈層片狀。

        圖2 ω 相的晶體結(jié)構(gòu)及其與β 相的位向關系Fig. 2 Crystal structure of ω phase and the orientation relationship with β phase

        應力誘發(fā)ω 相的形成機制在Samiee 等[38]眾多學者的研究中都有討論。ω 相與β 基體之間的取向為 [111]β//[0001]ω 和(1 1ˉ0 )β//(110)ω。Hsiung 等[39]認為應力誘發(fā)ω 相形成的機制是由1/2<111>全位錯分解形成的1/3<111>和2 個1/12<111>分位錯的移動形成的。

        由于應力誘發(fā)ω 相變往往伴隨著其他多種變形機制發(fā)生,因此,有此相變的合金會有很高的應變硬化率。

        1.5 剪切

        近年來的研究發(fā)現(xiàn),一些亞穩(wěn)β 型鈦合金[40-45]能通過剪切進行塑性變形。例如,Xu 等[40]在室溫壓縮變形的Ti-Nb-Ta-In/Cr 合金中、Wang 等[41]在冷軋的Ti-35Nb-2Ta-3Zr 合金中、Hao 等[42-45]在冷軋的Ti-24.0Nb-4.0Zr-7.9Sn 合金中、Yang 等[24]在大變形的Ti-22.40Nb-0.73Ta-2.00Zr-1.34O 合金中均發(fā)現(xiàn)了剪切帶的存在。當晶體處于某一特殊取向,并有大量的局部切應力集中在某一滑移面上時,晶體便會沿著該面進行剪切變形,并形成一個薄片狀變形微區(qū),即剪切帶。剪切帶內(nèi)的晶體發(fā)生轉(zhuǎn)動和嚴重變形,在Ti-24.0Nb-4.0Zr-7.9Sn 合金中甚至碎化成了納米晶,而與之相鄰的兩側(cè)晶體沒有明顯的變化,剪切變形通常發(fā)生在變形后期或嚴重變形區(qū)域。

        1.6 扭折

        扭折一般在合金經(jīng)歷較大程度的塑性變形時產(chǎn)生。Yang 等[10,46]發(fā)現(xiàn)在大變形時,應力誘發(fā)相變已經(jīng)停止,合金中的位錯塞積變得非常嚴重,靠應力誘發(fā)相變和位錯滑移已不能使其產(chǎn)生進一步的變形。為對抗材料破環(huán),晶體發(fā)生了局部彎曲,即扭折。扭折作為一種協(xié)調(diào)機制,其形成有助于協(xié)調(diào)相鄰晶?;蛘咄痪Я2煌课坏淖冃?,使應力松弛,進而提高合金的變形能力。

        2 亞穩(wěn)β 鈦合金變形機制的影響因素

        2.1 β 相穩(wěn)定性對變形機制的影響

        β 相的穩(wěn)定性通常與合金元素有關。一般來講,β 相穩(wěn)定元素含量的增加會提高β 相的穩(wěn)定性,抑制應力誘發(fā)相變,塑性變形機制會逐漸轉(zhuǎn)向以滑移為主導[47]。將能夠自β 相區(qū)淬火時把β 相保持到室溫而不發(fā)生馬氏體相變的β 相穩(wěn)定元素最小濃度稱為臨界密度,常見β 相穩(wěn)定元素的臨界質(zhì)量分數(shù)分別為Mo=10 %、V=15 %、Nb=36 %、Ta=40 %、Cr=7 %、Fe=5 %、Mn=6.4 %、Co=7 %、Ni=9 %,規(guī)定此時的鉬當量([Mo])均為10。通常,將[Mo]在10~25 時的合金歸為亞穩(wěn)β 鈦合金,但是對于不同的合金體系,這一規(guī)定并不是很嚴格。

        但β 相的穩(wěn)定性對變形機制的影響還存在更深層的原因。2003 年,Saito 等[48]設計出了一類具有眾多優(yōu)異性能的亞穩(wěn)β 鈦合金,這些鈦合金同時具有高強度、低彈性模量、超彈性和優(yōu)異的冷加工性能,并被命名為“Gum metal”(即“橡膠金屬”)。經(jīng)過第一性原理計算分析,這類材料同時滿足3 個參數(shù):(1) 平均價電子數(shù)e/a=4.24;(2) 鍵級Bo=2.87;(3) d 電子軌道能級Md=2.45。這3 個參數(shù)被認為與β 相的穩(wěn)定性有關。其中最具有代表性的Ti-23.0Nb-0.7Ta-2.0Zr-1.2O 合金在經(jīng)過變形量為90%的冷變形后抗拉強度由1 000 MPa 提升至1 200 MPa,彈性模量由100 GPa 下降至55 GPa 左右,并且在-173~227 ℃的大溫度區(qū)間內(nèi)都保持恒定。在對“橡膠金屬”的冷變形行為進行研究之后,Saito等[49-51]開展了一系列的實驗驗證,認為此系列合金之所以具有如此優(yōu)越的力學性能,是因為其獨特的“無位錯”變形機制,其變形是通過“橡膠金屬”大尺寸缺陷實現(xiàn)的,沒有位錯滑移、應力誘發(fā)相變和孿生的參與,這些結(jié)論引起了其他研究人員極大的興趣。隨著研究的深入,越來越多的研究者在合金中發(fā)現(xiàn)了傳統(tǒng)的變形產(chǎn)物。

        盡管“無位錯”塑性變形機制存在巨大的爭議,但其設計高強度低彈性模量亞穩(wěn)β 鈦合金的方法為開發(fā)新合金和研究變形機制提供了新的思路。由于Bo 與Md 值決定了β 相的穩(wěn)定性,根據(jù)這兩個值可以預測鈦合金材料的變形產(chǎn)物。圖3 為描述相穩(wěn)定性的Bo-Md 圖[52]。許多研究開始利用Bo-Md圖對鈦合金進行設計。在圖3 中,合金的Bo 與Md值分布在不同區(qū)域,意味著在變形時有可能發(fā)生此區(qū)域?qū)淖冃螜C制。

        圖3 描述相穩(wěn)定性的Bo-Md 圖[52]Fig.3 Bo-Md diagram for describing the phase stability[52]

        2.2 變形條件對變形機制的影響

        在傳統(tǒng)的α+β 雙相鈦合金和β 鈦合金領域,Khan 等[53]、Nemat-Nasser 等[54]、Hokka 等[55]和Lee等[56]進行了動態(tài)變形實驗,發(fā)現(xiàn)流變應力對高應變速率和高溫敏感。由于亞穩(wěn)定β 鈦合金在形變時易發(fā)生相變,各相對于材料力學性能的影響也不盡相同。再加上不同生成相之間還會互相影響,變形溫度、應變速率和變形量的改變會令合金的變形機制更加復雜。

        變形溫度對亞穩(wěn)β 鈦合金變形機制的影響是由于變形溫度能改變β 相的穩(wěn)定性。例如Zhan 等[57]研究了20~600 ℃內(nèi)Ti-Nb-Zr-Mo-Sn 合金在高應變速率下的力學性能,發(fā)現(xiàn)由于β 相的穩(wěn)定性隨著溫度升高而升高,變形機制也由低溫時的孿生、位錯滑移和應力誘發(fā)相變共存變?yōu)楦邷貢r的純位錯滑移。同樣,Samiee 等[38]研究了亞穩(wěn)β 鈦合金Ti-10V-3Fe-3Al 在725 ℃下以10-3s-1的應變速率的壓縮行為。發(fā)現(xiàn)在725 ℃的高溫下β 相過于穩(wěn)定,很難生成α"相。眾多研究指出合金在低溫下β 相的穩(wěn)定性較低,這造成了變形機制的多樣性。

        變形量也會對變形機制產(chǎn)生影響,Tsuchiya 等[58]發(fā)現(xiàn),TiNi 形狀記憶合金變形量超過25%時馬氏體相變就會受到抑制。Yang 等[10]發(fā)現(xiàn)Ti-Nb-Ta-Zr-O 系合金的變形機制不僅與β 相的穩(wěn)定性有關,還與材料的變形量有關。當真應變在35%以下時,材料變形會產(chǎn)生{332}<113>孿晶、{112}<111>孿晶、應力誘發(fā)α"和ω 相以及位錯滑移,當真應變?yōu)?5%~77%時則發(fā)生扭折,當真應變?yōu)?7%~149%時材料的形變機制則為剪切,這種主要變形機制的變化是由于在材料的變形過程中缺陷密度和晶體取向發(fā)生了轉(zhuǎn)變,并且由于主導的變形機制改變,材料還表現(xiàn)出了多次應變硬化現(xiàn)象。

        應變速率對亞穩(wěn)β 鈦合金的變形機制也有明顯的影響,有研究[59-62]提出,高速壓縮會導致孿晶密度高。例如,Ahmed 等[62]在室溫下通過改變應變速率(10-3、10-1、10 s-1和102s-1)測試了亞穩(wěn)β 鈦合金Ti-10V-3Fe-3Al-0.27O 的壓縮性能。發(fā)現(xiàn)材料的變形機制會隨著應變速率的提高而變化,在應變速率低于10-3s-1時,主要的變形機制為應力誘發(fā)α"馬氏體相變,而應變速率處于10-3s-1至10-1s-1的區(qū)間時則出現(xiàn)了SIMα"與{332}<113>孿生的競爭,在應變速率高于10 s-1的條件下,主要的變形機制轉(zhuǎn)變?yōu)閧332}<113>孿生。Ji 等[61]研究了 Ti-10Mo-1Fe 合金在2.8×10-5s-1至2.8×10-1s-1應變速率下的拉伸性能,發(fā)現(xiàn)位錯密度隨應變速率提高而增大,在應變量較低時孿晶密度也隨應變速率提高而有明顯增大。應變速率與合金β 相的穩(wěn)定性對變形機制的影響如圖4 所示。Ahmed 等[62]認為,之所以孿晶的含量隨著應變速率的提高而增大,是因為在高的應變速率下,位錯密度和位錯移動速度也會上升,而高的位錯密度有利于孿晶的成核。這種理論在其他研究[11]中也有涉及,較高的位錯密度會降低孿晶成核所需的活化能。

        圖4 應變速率與變形機制變化示意圖,虛線為應變速率臨界值的變化趨勢[62]Fig. 4 Schematic diagram for the change of deformation mechanism with respect to strain rate. dashed lines indicate the change trend of the strain rate critical values[62]

        Sadeghpour 等[1]研究了經(jīng)過0.7×10-4s-1至0.7×10-1s-1的應變速率壓縮后的Ti-4Al-7Mo-3V-3Cr 合金,在所有條件下都觀察到了應力誘發(fā)的α"馬氏體,并且還發(fā)現(xiàn)材料的馬氏體觸發(fā)應力和屈服強度都與應變速率的對數(shù)呈線性增加的關系。眾多研究[1,63-64]結(jié)果表明,應變速率過高往往會抑制應力誘發(fā)α"馬氏體的生成,這種作用是通過生成α"馬氏體的自由能變化來影響的,此關系如下式[63]所示:

        式中:ΔG為SIMα"轉(zhuǎn)變的自由能變化;ΔH為焓變;T為測試溫度;ΔS為熵變; δEirr為相界在運動中克服摩擦力所做的功;ΔEel為在基體中產(chǎn)生的內(nèi)部彈性勢能,由相變時的體積變化造成;ρ為密度; σ為軸向應力; εtr是由馬氏體在應力方向上的晶格畸變造成的應變;1 /2ρ(σ2/EM-σ2/Eβ)為在彈性變形期間由于β 相和a"相的彈性模量的差異而產(chǎn)生的彈性勢能;EM為 SIMα"的彈性模量;Eβ為β 相的彈性模量。式中,ΔH和ΔS均與應變速率無關;由于β 相和a′′相的密度幾乎相同,ΔEel可以忽略不計; εtr、EM和Eβ不受應變速率的影響;只有 δEirr隨著應變速率的增加而增加[59],是唯一與應變速率顯著相關的變量,并且由于應變速率攀升,ΔG將隨 δEirr的增加而增加,所以生成α"的能量壁壘升高,生成α"的難度上升。

        上述研究充分說明了變形條件對亞穩(wěn)β 鈦合金變形機制的影響,作為此合金的主要變形機制之一,{332}<113>孿生對變形溫度和應變速率的變化都比較敏感,應力誘發(fā)α"相也會受到變形條件的影響。變形機制在材料變形中的轉(zhuǎn)變則會進一步造成材料力學性能的變化。

        3 結(jié) 論

        亞穩(wěn)β 鈦合金雖然具有優(yōu)異的力學性能,但其變形機制復雜,近年來許多研究者就亞穩(wěn)β 鈦合金的變形機制開展了深入的研究。其變形機制與力學性能的關系逐步厘清,最新的研究成果甚至實現(xiàn)了將此類合金的高強度和高塑性同時保留。但這些研究的實驗條件大多集中于低速準靜態(tài)變形,而對高速動態(tài)變形的研究才剛剛起步??紤]到亞穩(wěn)β 鈦合金在航空航天領域的應用中經(jīng)常處在較為極端的服役條件,亞穩(wěn)β 鈦合金在動態(tài)變形的領域還需開展進一步研究。

        猜你喜歡
        孿晶馬氏體鈦合金
        中低碳系列馬氏體不銹鋼開發(fā)與生產(chǎn)
        山東冶金(2022年1期)2022-04-19 13:40:20
        鎂金屬孿晶變形的實驗和理論模型研究進展
        激光制備預壓應力超高強韌馬氏體層的組織與性能
        “神的金屬”鈦合金SHINE YOUR LIFE
        中國自行車(2018年8期)2018-09-26 06:53:32
        硅量子點中的孿晶對其電子結(jié)構(gòu)和光學性能的影響
        鈦合金板鍛造的工藝實踐
        四川冶金(2017年6期)2017-09-21 00:52:30
        Fe-C-Mn-Si-Cr的馬氏體開始轉(zhuǎn)變點的熱力學計算
        上海金屬(2016年1期)2016-11-23 05:17:24
        關于Fe-1.4C合金馬氏體相變的研究
        醫(yī)用鈦合金的研究與應用
        TC17鈦合金超大規(guī)格棒材的制備
        成人在线视频亚洲国产| 亚洲字幕av一区二区三区四区| 亚洲熟妇无码八av在线播放| 日韩专区欧美专区| 国产麻豆放荡av激情演绎| 日本一区二区三区视频免费在线| 乱色精品无码一区二区国产盗| av人摸人人人澡人人超碰小说| 美女黄频视频免费国产大全| 黄片国产一区二区三区| 国产欧美一区二区精品久久久| 伊人久久无码中文字幕| 亚洲性无码av在线| 日本综合视频一区二区| 午夜精品射精入后重之免费观看| 国产午夜影视大全免费观看| 成人免费丝袜美腿视频| 人妻少妇精品专区性色anvn| 国产精品99久久久久久猫咪| 亚洲欧美日韩在线观看一区二区三区| 国产精品狼人久久久影院| 蜜桃臀av一区二区三区| 亚洲日韩av无码| 亚洲AV无码成人品爱| 丰满熟女人妻一区二区三区| 国产无套乱子伦精彩是白视频| 99精品热这里只有精品| 亚洲无AV码一区二区三区| 亚洲国产国语对白在线观看| 亚洲人成无码区在线观看| 精品国产高清a毛片无毒不卡 | mm在线精品视频| 男女激情视频网站免费在线| 国模雨珍浓密毛大尺度150p| 在线免费观看国产精品| av免费一区在线播放| 亚洲国产精品一区二区成人片国内| 亚洲高潮喷水无码av电影| 一区在线播放| 一区二区三区四区草逼福利视频| 国产真实强被迫伦姧女在线观看 |