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

加載速率對(duì)40Cr 鋼Ⅱ型動(dòng)態(tài)斷裂特性的影響*

2021-09-10 02:54:14范昌增許澤建何曉東黃風(fēng)雷
爆炸與沖擊 2021年8期
關(guān)鍵詞:裂紋信號(hào)

范昌增,許澤建,何曉東,黃風(fēng)雷

(北京理工大學(xué)爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100081)

40Cr 鋼是航空航天以及國(guó)防領(lǐng)域中使用最廣泛的高強(qiáng)鋼之一,因其淬透性良好,經(jīng)熱處理后有良好的綜合力學(xué)性能和低溫沖擊韌性,常被用于制作機(jī)翼大梁、坦克變速器扭軸以及裝甲車輛的軸承類零件等[1],但這些構(gòu)件在服役過(guò)程中容易受到?jīng)_擊剪切載荷的作用而發(fā)生Ⅱ型斷裂失效。因此,工程材料在Ⅱ型加載下的動(dòng)態(tài)斷裂問(wèn)題日益受到人們的關(guān)注,Ⅱ型裂紋在動(dòng)態(tài)載荷下的起裂及擴(kuò)展行為也成為力學(xué)和材料學(xué)科的交叉前沿課題。20 世紀(jì)80 年代,Kalthoff 等[2-5]發(fā)現(xiàn)42CrMo4 鋼在低加載速率下Ⅱ型裂紋發(fā)生拉伸型失效,裂紋擴(kuò)展方向與韌帶方向成?70°夾角,而在高加載速率下裂紋沿韌帶方向以絕熱剪切帶形式發(fā)生失效。Zhou 等[6-8]的研究表明,材料的失效模式轉(zhuǎn)變行為更多地依賴于材料特性,而不是試樣的幾何形狀。Ravi-Chandar 等[9-10]發(fā)現(xiàn),隨著沖擊速度的增大,聚碳酸酯(PC)的失效模式存在從延性到脆性再到延性的轉(zhuǎn)變,且第二次模式轉(zhuǎn)變會(huì)導(dǎo)致裂紋尖端附近形成剪切帶。Mason 等[11]對(duì)C-300 鋼單缺口試件進(jìn)行動(dòng)態(tài)剪切試驗(yàn),研究了不同失效形式之間的轉(zhuǎn)換。Chu 等[12]對(duì)金屬材料的韌脆模式轉(zhuǎn)變進(jìn)行了相場(chǎng)損傷演化,發(fā)現(xiàn)應(yīng)力波的傳播在脆性-韌性破壞模式轉(zhuǎn)變中起著重要作用。在實(shí)驗(yàn)測(cè)試方面,莫爾云紋法、全息法和光彈性法等光測(cè)力學(xué)方法[13-15]已被應(yīng)用到斷裂試驗(yàn)中[16-20],并取得了滿意的效果。但是,目前對(duì)于高強(qiáng)鋼Ⅱ型裂紋在沖擊載荷下失效行為的研究大多仍采用類似Kalthoff 的加載方法或單桿加載方法對(duì)單邊裂紋試樣進(jìn)行測(cè)試[21],這類加載方法不易于對(duì)斷口形貌特征的轉(zhuǎn)變行為進(jìn)行定量控制和精確測(cè)量。近期,許澤建等[22-24]提出了一種新型Ⅱ型動(dòng)態(tài)斷裂實(shí)驗(yàn)方法,該方法可以實(shí)現(xiàn)對(duì)加載速率的精確控制,從而可以對(duì)材料的失效模式轉(zhuǎn)變行為進(jìn)行系統(tǒng)的量化研究。他們通過(guò)對(duì)TC4 材料的研究,發(fā)現(xiàn)失效模式轉(zhuǎn)變(FMT)不是突變的結(jié)果,而是受不同微觀機(jī)制主導(dǎo)的漸變過(guò)程。隨著加載速率的增加,材料的失效模式及斷口特征均發(fā)生改變,是多種失效機(jī)制協(xié)同作用的結(jié)果。以上結(jié)果為Ⅱ斷裂的失效行為研究提供了新的方法和思路。

本文中采用該新型Ⅱ型動(dòng)態(tài)斷裂實(shí)驗(yàn)方法,基于分離式霍普金森壓桿(SHPB)技術(shù)[25],對(duì)40Cr 高強(qiáng)鋼開(kāi)展不同加載速率下的動(dòng)態(tài)斷裂實(shí)驗(yàn);采用實(shí)驗(yàn)-數(shù)值方法確定材料的動(dòng)態(tài)斷裂韌性值,并通過(guò)高速攝影觀察裂紋的起裂及擴(kuò)展過(guò)程;結(jié)合試樣斷口的宏、微觀形貌特征分析,對(duì)材料的失效模式及微觀機(jī)理進(jìn)行研究。

1 實(shí)驗(yàn)方法

采用SHPB 系統(tǒng)[26]對(duì)Ⅱ型裂紋進(jìn)行動(dòng)態(tài)加載和測(cè)量,并采用實(shí)驗(yàn)-數(shù)值方法對(duì)裂尖的應(yīng)力強(qiáng)度因子和斷裂韌性進(jìn)行確定。實(shí)驗(yàn)時(shí),通過(guò)入射桿對(duì)特殊設(shè)計(jì)的Ⅱ型斷裂試樣進(jìn)行加載,并在裂尖形成Ⅱ型載荷,試樣的幾何外形及尺寸如圖1 所示。入射桿、透射桿和試樣上均貼有應(yīng)變片,分別記錄入射波 εi(t) 、透射波 εt(t) 信號(hào)和試樣的起裂信號(hào)。為避開(kāi)裂尖的塑性區(qū),試樣上應(yīng)變片的粘貼位置位于裂尖的一側(cè),應(yīng)變片中心點(diǎn)位于距離裂尖約3 mm 處(見(jiàn)圖1)。根據(jù)一維應(yīng)力波理論,可獲得壓桿與試樣接觸面處的載荷和位移曲線,如下式:

圖1 Ⅱ型斷裂試樣幾何尺寸(單位:mm)Fig. 1 Geometric dimensions of the mode Ⅱfracture specimen (unit: mm)

把實(shí)驗(yàn)所測(cè)入射桿的應(yīng)力波波形與實(shí)際位移邊界條件等作為輸入量,對(duì)試樣的受載過(guò)程進(jìn)行數(shù)值模擬,并通過(guò)位移法[22]計(jì)算裂尖的應(yīng)力強(qiáng)度因子曲線KⅡ(t),如下式:

式中:Fin、Fout為入射桿與透射桿端部的載荷,Uin、Uout為入射桿與透射桿端部的位移,A為橫截面積,E為楊氏模量,c0為桿的縱波波速,r為距裂尖的距離,G為剪切模量,μ為泊松比,u(r,t)為裂尖位移場(chǎng)。

在此基礎(chǔ)上,結(jié)合試樣的起裂信號(hào)可以確定材料的Ⅱ型動(dòng)態(tài)斷裂韌性。

另外,裂尖的加載速率K˙Ⅱd可由下式確定:

式中:tf為試樣的起裂時(shí)間,KⅡd為動(dòng)態(tài)斷裂韌性。

實(shí)驗(yàn)材料為40Cr 高強(qiáng)鋼,元素成分見(jiàn)表1,熱處理工藝為:40Cr 正火850 ℃,冷卻方式為空冷,40Cr 淬火溫度為780 ℃,冷卻方式為水冷,回火溫度為200 ℃,冷卻方式為空冷。針對(duì)同批次、同爐熱處理的材料進(jìn)行了力學(xué)性能測(cè)試,得到40Cr 材料的力學(xué)性能測(cè)試參數(shù),見(jiàn)表2,表中ρ 為密度,E為彈性模量,μ為泊松比,σb為抗拉強(qiáng)度。實(shí)驗(yàn)中通過(guò)控制氣壓大小和子彈長(zhǎng)度來(lái)控制加載速率。

表1 高強(qiáng)鋼40Cr 的元素成分及質(zhì)量分?jǐn)?shù)Table 1 Composition and mass fraction of high-strength steel 40Cr

表2 高強(qiáng)鋼40Cr 的力學(xué)性能參數(shù)Table 2 Mechanical properties of high strength steel 40Cr

在1.08~5.53 TPa·m1/2/s 加載速率范圍內(nèi)對(duì)40Cr 試樣進(jìn)行動(dòng)態(tài)加載,為避免壓桿中反射的應(yīng)力波對(duì)試樣造成重復(fù)加載,實(shí)驗(yàn)中采用了單脈沖加載技術(shù)[27]。圖2 為在4.86 TPa·m1/2/s 加載速率下獲得的原始波形。由圖2 可知,反射波在入射桿端部經(jīng)單脈沖系統(tǒng)吸收后,剩余的反射波成分很少,幅值接近于零,因此可以實(shí)現(xiàn)對(duì)試樣的單次加載。通過(guò)該技術(shù),可以避免壓桿中因應(yīng)力波的反射給試樣造成重復(fù)加載,有利于對(duì)材料的失效模式進(jìn)行確定和研究。另外,試樣的起裂時(shí)間采用應(yīng)變片法[22]進(jìn)行測(cè)定,原理為當(dāng)應(yīng)力波作用于試樣時(shí),試樣先發(fā)生彈性變形,此時(shí)應(yīng)變信號(hào)從零開(kāi)始線性增加;當(dāng)裂尖出現(xiàn)塑性區(qū)并開(kāi)始起裂時(shí),模擬中試樣應(yīng)變信號(hào)與實(shí)測(cè)信號(hào)發(fā)生分離,其分離點(diǎn)即為起裂時(shí)間。這部分時(shí)間包括了裂尖起裂時(shí)的應(yīng)力波傳到應(yīng)變片的時(shí)間[23],見(jiàn)下式:

圖2 典型實(shí)驗(yàn)信號(hào)(4.86 TPa·m1/2/s)Fig. 2 Typical experimental signals (4.86 TPa·m1/2/s)

式中:D為應(yīng)變片中心到裂尖的距離,B為試樣的厚度,c0為卸載波的波速, 一般為彈性波縱波波速。

2 實(shí)驗(yàn)結(jié)果

在1.25 TPa·m1/2/s 加載速率下的實(shí)驗(yàn)波形如圖3 所示,試樣發(fā)生了拉伸型起裂。由圖3 可知,反射波電壓先上升到0.2 V,然后驟降;與此同時(shí),透射波信號(hào)也不斷降低。在138 μs 左右,反射波恢復(fù)到0.2 V并保持穩(wěn)定,透射波對(duì)應(yīng)的信號(hào)上升至?0.1 V。結(jié)合高速攝影觀察分析知,反射波與透射波信號(hào)下降階段對(duì)應(yīng)于應(yīng)力波對(duì)試樣的加載過(guò)程,此過(guò)程中試樣始終與壓桿保持接觸,且試樣發(fā)生起裂但尚未完全斷裂。當(dāng)反射波信號(hào)上升至0.2 V、透射波信號(hào)上升至?0.1 V 時(shí),試樣完全起裂,且通過(guò)高速攝影觀察可知,試樣單邊支撐端先發(fā)生斷裂。此后透射端載荷未降至0 V,說(shuō)明該階段透射桿與試樣支撐端仍保持接觸狀態(tài)。

圖3 試樣拉伸起裂的典型信號(hào)(1.25 TPa·m1/2/s)Fig. 3 Typical signals of tensile fracture initiation of a specimen (1.25 TPa·m1/2/s)

在4.86 TPa·m1/2/s 加載速率下的實(shí)驗(yàn)波形見(jiàn)圖4,試樣發(fā)生絕熱剪切型斷裂。由圖4 可知,反射波電壓在上升到0.45 V 后開(kāi)始緩慢下降;與此同時(shí),透射波信號(hào)也較為平緩地下降。在38 μs 左右,反射波恢復(fù)到0.5 V 并保持穩(wěn)定,透射波對(duì)應(yīng)的信號(hào)上升到0 V。結(jié)合高速攝影(見(jiàn)圖5)觀察分析知,反射波與透射波信號(hào)下降階段對(duì)應(yīng)于應(yīng)力波對(duì)試樣的加載過(guò)程,此過(guò)程中試樣始終與壓桿保持接觸,且試樣發(fā)生起裂但尚未完全斷裂。當(dāng)反射波電壓值保持為0.5 V、透射波電壓保持為0 V 時(shí),試樣完全起裂,且通過(guò)高速攝影觀察可知試樣雙邊支撐端同時(shí)發(fā)生斷裂。此后透射波電壓降為0 V,說(shuō)明該階段透射桿與試樣兩支撐端已脫離接觸。

圖4 試樣絕熱剪切起裂的典型信號(hào)(4.86 TPa·m1/2/s)Fig. 4 Typical signals of ASB initiation of a specimen (4.86 TPa·m1/2/s)

以上兩次實(shí)驗(yàn)的加載過(guò)程如圖5 所示,其中圖5(a)、(e)分別為試樣初始受載時(shí)刻;圖5(b)、(f)為試樣裂尖開(kāi)始起裂的時(shí)刻,可以觀察到此時(shí)裂尖局部出現(xiàn)白亮區(qū)域,說(shuō)明已經(jīng)發(fā)生了小范圍屈服現(xiàn)象;圖5(c)、(g)為不同加載速率下試樣完全斷裂的失效形態(tài)。根據(jù)以上分析可將加載過(guò)程分為2 個(gè)階段:第1 階段為從試樣開(kāi)始受載到臨界起裂(即出現(xiàn)白點(diǎn))的階段,此階段表現(xiàn)為試樣吸收能量發(fā)生彈性變形,同時(shí)在裂尖區(qū)域形成小范圍塑性變形和微裂紋;第2 階段是從裂尖處的微裂紋集聚形成新的界面(即發(fā)生起裂)至試樣以拉伸型或絕熱剪切型失效模式發(fā)生完全斷裂。在第2 階段,隨著加載速率的增加,裂紋擴(kuò)展角及擴(kuò)展路徑也隨之發(fā)生變化。該現(xiàn)象說(shuō)明裂紋擴(kuò)展過(guò)程中裂尖的應(yīng)力狀態(tài)和材料特性會(huì)發(fā)生改變并決定著裂紋面的發(fā)展趨勢(shì)。

圖5 拉伸型起裂(1.25 TPa·m1/2/s)與ASB 型起裂(4.86 TPa·m1/2/s)的高速攝影圖Fig. 5 High-speed photographic images of tensile fracture initiation (1.25 TPa·m1/2/s) and ASB initiation (4.86 TPa·m1/2/s)

3 分析及討論

3.1 有限元模擬及驗(yàn)證

有限元模擬采用ABAQUS/Explicit,此軟件的顯示動(dòng)力學(xué)模塊在計(jì)算高速?zèng)_擊加載和接觸性問(wèn)題時(shí)具有顯著的優(yōu)勢(shì)。根據(jù)實(shí)際工況和尺寸進(jìn)行1∶1建模,并將實(shí)測(cè)的入射桿應(yīng)力波作為輸入量。入射桿、透射桿和試樣的材料均定義為線彈性模型,參數(shù)見(jiàn)表3;入射桿、透射桿與試樣的接觸端面均設(shè)為法向硬接觸,并在試樣支撐端的側(cè)面施加橫向位移約束以模擬實(shí)際工況中夾具的限位作用。入射桿與透射桿均采用八節(jié)點(diǎn)六面體單元;為提高計(jì)算精度,在裂尖區(qū)域采用六面體-楔形體-四面體的過(guò)渡方式劃分網(wǎng)格,并進(jìn)行網(wǎng)格細(xì)化,使裂紋端部與裂紋半長(zhǎng)處的網(wǎng)格尺寸比值控制在0.001~0.02 之間,此時(shí)裂尖彈性應(yīng)變場(chǎng)的計(jì)算誤差約為1%。模型裝配圖如圖6 所示。

圖6 試樣裂尖網(wǎng)格細(xì)化Fig. 6 Mesh refinement of the specimen crack tips

表3 入射桿、透射桿的力學(xué)性能參數(shù)Table 3 Mechanical properties of the incident and transmission bars

將模擬結(jié)果中試樣粘貼應(yīng)變片區(qū)域的應(yīng)變信號(hào)與實(shí)測(cè)的應(yīng)變信號(hào)進(jìn)行對(duì)比,如圖7 所示。當(dāng)材料處于線彈性階段時(shí),試樣的模擬真實(shí)應(yīng)變與實(shí)測(cè)的試樣應(yīng)變信號(hào)基本吻合,佐證了材料參數(shù)的準(zhǔn)確性以及有限元模擬的可靠性。當(dāng)時(shí)材料處于塑性變形階段,且實(shí)測(cè)應(yīng)變?yōu)?.005 時(shí),模擬信號(hào)與實(shí)測(cè)信號(hào)發(fā)生分離,該分離點(diǎn)作為起裂時(shí)刻。由圖8 可知,隨著子彈速度的增加,試樣的起裂時(shí)間明顯提前,且基本保持線性關(guān)系。采用位移法確定的應(yīng)力強(qiáng)度因子(dynamic stress intensity factor, DSIF)曲線如圖9 所示。根據(jù)起裂時(shí)間,在DSIF 曲線上確定了試樣Ⅱ型斷裂韌性,結(jié)果見(jiàn)表4。圖10 中還給出了文獻(xiàn)[21]采用單桿加載單邊裂紋試樣方法所測(cè)出的40Cr 鋼動(dòng)態(tài)斷裂韌性結(jié)果。從圖10 中可看出,新型Ⅱ型加載方式測(cè)得的動(dòng)態(tài)斷裂韌性DFT(dynamic fracture toughness)明顯高于文獻(xiàn)結(jié)果,這可能是由于文獻(xiàn)中的單桿加載方式在對(duì)試樣加載的過(guò)程中,單邊裂紋在彎矩作用下出現(xiàn)了Ⅰ型斷裂現(xiàn)象(即張開(kāi)型斷裂),因此裂尖起裂時(shí)受到Ⅰ型斷裂的影響,導(dǎo)致所測(cè)得的Ⅱ型動(dòng)態(tài)斷裂韌性值偏低。

圖7 試樣的實(shí)測(cè)應(yīng)變與模擬應(yīng)變Fig. 7 The measured strain and simulated strain of the specimen

圖8 子彈速度與起裂時(shí)間的關(guān)系Fig. 8 Relationship between bullet velocity and crack initiation time

圖9 不同加載速率下的DSIF 曲線Fig. 9 DSIF curves at different loading rates

圖10 不同加載方式下的DFT 與加載速率關(guān)系圖Fig. 10 Relationship between DFT and loading rate under different loading methods

表4 高強(qiáng)鋼40Cr 的動(dòng)態(tài)斷裂韌性值Table 4 Dynamic fracture toughness values of high strength steel 40Cr

綜上所述,當(dāng)加載速率在1.08~5.53 TMPa·m1/2/s范圍時(shí),40Cr 材料Ⅱ型斷裂韌性與加載速率呈現(xiàn)正相關(guān)趨勢(shì)。其中當(dāng)≤2.92 TPa·m1/2/s時(shí),40Cr 試樣均屬于拉伸型失效模式,裂紋擴(kuò)展方向與韌帶方向夾角由?60°逐漸減小至?45°,動(dòng)態(tài)斷裂韌性(DFT)與加載速率的擬合表達(dá)式為y=10.45x+18.83;當(dāng)≥3.19 TPa·m1/2/s時(shí),試樣均表現(xiàn)為絕熱剪切型失效模式,裂紋擴(kuò)展與韌帶方向夾角約為?16°,擬合表達(dá)式為y=2.09x+43.39。

3.2 斷口形貌分析

采用Phenom Pro掃描電鏡對(duì)試樣斷口形貌進(jìn)行觀察,根據(jù)形貌特征分為脆性、韌性以及絕熱剪切型3種斷口。由圖11可知,拉伸型失效模式下存在脆性、韌性兩種斷口形貌,而絕熱剪切型失效模式下僅存在絕熱剪切型斷口形貌特征。在不同的斷口形貌下,DFT與的線性擬合參數(shù)見(jiàn)表5。

表5 不同斷口形貌下DFT與加載速率的線性擬合參數(shù)Table 5 Linear fitting parameters of DFT and the loading rates under different fracture morphology

圖11 不同斷口形貌特征的DFT 與加載速率關(guān)系圖Fig.11 Relationship between DFT and loading rate of different fracture morphology characteristics

在中加載速率(2.12 TPa·m1/2/s≤≤2.92 TPa·m1/2/s)下,40Cr 試樣斷口宏觀上呈亮灰色,斷面有平滑的剪切滑移帶且有金屬光澤。將圖13局部區(qū)域放大后可觀察到明顯的韌性斷口形貌特征:拋物線型的剪切韌窩以及平滑的長(zhǎng)條狀形貌。其原因一般為微孔洞在剪切應(yīng)力作用下萌生、成核、長(zhǎng)大、集聚形成剪切韌窩,隨著加載速率的增加,形變不斷加劇,由此產(chǎn)生的絕熱溫升使剪切韌窩軟化,同時(shí)在剪切應(yīng)力作用下被快速剪切為平滑的長(zhǎng)條狀形貌,甚至形成無(wú)特征的平坦面,該區(qū)域面積占比伴隨著加載速率的增大而不斷增大。因此,在中加載速率下,40Cr 材料斷口屬于由應(yīng)變/應(yīng)變率硬化以及熱軟化機(jī)制共同主導(dǎo)的韌性斷裂模式。

圖13 40Cr 試樣韌性斷裂特征(2.92 TPa·m1/2/s)Fig.13 Ductile fracture characteristics of 40Cr specimen (2.92 TPa·m1/2/s)

圖14 40Cr 試樣絕熱剪切型斷裂特征(3.87 TPa·m1/2/s)Fig.14 Adiabatic shear fracture characteristics of 40Cr specimen (3.87 TPa·m1/2/s)

由以上分析可知,不同加載速率下,材料的斷口形貌特征也不同。低加載速率下,40Cr 材料斷口屬于由應(yīng)變率硬化機(jī)制主導(dǎo)的脆性斷裂模式;中加載速率下,屬于由應(yīng)變/應(yīng)變率硬化以及熱軟化機(jī)制共同主導(dǎo)的韌性斷裂模式;高加載速率下,屬于由熱軟化機(jī)制主導(dǎo)的絕熱剪切型斷裂模式。

4 結(jié) 論

采用新型Ⅱ型動(dòng)態(tài)斷裂測(cè)試技術(shù),對(duì)40Cr 鋼進(jìn)行了不同加載速率下的Ⅱ型動(dòng)態(tài)斷裂特性研究,得到如下結(jié)論:

(1)在加載速率范圍(1.08~5.53 TPa·m1/2/s)內(nèi),40Cr 材料的Ⅱ型斷裂韌性值與加載速率基本呈正相關(guān)趨勢(shì)。

(2)在1.08~5.53 TPa·m1/2/s的加載速率范圍內(nèi),40Cr 存在拉伸型失效模式向絕熱剪切型失效模式轉(zhuǎn)變,且臨界加載速率為2.92 TPa·m1/2/s。

(3)當(dāng)K˙Ⅱd≤1.70 TPa·m1/2/s時(shí),40Cr 材料的斷口表現(xiàn)為脆性斷裂;當(dāng)K˙Ⅱd介于2.12~2.92 TPa·m1/2/s之間時(shí),斷口表現(xiàn)為韌性斷裂;當(dāng)K˙Ⅱd≥3.19 TPa·m1/2/s時(shí),斷口表現(xiàn)為絕熱剪切型形貌特征。

(4)脆性斷裂由應(yīng)變率硬化機(jī)制主導(dǎo),韌性斷裂由應(yīng)變/應(yīng)變率硬化和熱軟化機(jī)制共同主導(dǎo),絕熱剪切型斷裂由熱軟化機(jī)制主導(dǎo)。

猜你喜歡
裂紋信號(hào)
裂紋長(zhǎng)度對(duì)焊接接頭裂紋擴(kuò)展驅(qū)動(dòng)力的影響
一種基于微帶天線的金屬表面裂紋的檢測(cè)
信號(hào)
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
孩子停止長(zhǎng)個(gè)的信號(hào)
微裂紋區(qū)對(duì)主裂紋擴(kuò)展的影響
基于LabVIEW的力加載信號(hào)采集與PID控制
一種基于極大似然估計(jì)的信號(hào)盲抽取算法
預(yù)裂紋混凝土拉壓疲勞荷載下裂紋擴(kuò)展速率
主站蜘蛛池模板: A级毛片无码久久精品免费| 精品中文字幕一区在线| 国产迷奸在线看| 国产91蝌蚪窝| 欧美精品成人一区二区在线观看| 欧美成人影院亚洲综合图| 国产成人毛片| 亚洲香蕉在线| 国产视频 第一页| 国产人妖视频一区在线观看| 日韩精品久久无码中文字幕色欲| 东京热高清无码精品| 国产十八禁在线观看免费| 亚洲综合18p| 亚洲AV无码不卡无码| 99久久精品免费观看国产| 久青草免费视频| 毛片手机在线看| 日韩在线播放中文字幕| 国产女主播一区| 四虎影院国产| 久久99久久无码毛片一区二区| 国产69囗曝护士吞精在线视频| 波多野结衣在线se| 亚洲香蕉久久| 国产AV毛片| 成人av专区精品无码国产| www精品久久| 成人国产一区二区三区| 亚洲第一天堂无码专区| 精品国产香蕉伊思人在线| 第九色区aⅴ天堂久久香| 亚洲人成网站色7777| 国产尤物jk自慰制服喷水| 国产精品天干天干在线观看| 日本免费福利视频| 99精品国产电影| 日韩精品亚洲一区中文字幕| 成人在线观看一区| 国产一线在线| 国产一级毛片yw| 丁香六月综合网| 尤物亚洲最大AV无码网站| 婷婷亚洲视频| 国产高清无码第一十页在线观看| 国产精品入口麻豆| 精品人妻一区二区三区蜜桃AⅤ| 精品第一国产综合精品Aⅴ| 视频一区视频二区日韩专区| 欧美黄网在线| 国产福利小视频在线播放观看| 国产福利大秀91| 久久九九热视频| 久久精品国产亚洲麻豆| 欧美色视频日本| 呦视频在线一区二区三区| 日韩经典精品无码一区二区| 亚洲一区二区黄色| 99精品福利视频| 69av免费视频| 91精品国产情侣高潮露脸| 中文字幕永久视频| 91国内外精品自在线播放| 国产免费精彩视频| 67194亚洲无码| 欧美成人综合在线| 国产一区二区三区免费| 国产精品福利一区二区久久| 51国产偷自视频区视频手机观看| 免费一极毛片| 国产成人免费手机在线观看视频| 人人爽人人爽人人片| 粉嫩国产白浆在线观看| 亚洲视频一区| 特级毛片8级毛片免费观看| 99精品国产高清一区二区| 国产精品99一区不卡| 婷婷亚洲最大| 久久久噜噜噜| 久久久精品久久久久三级| 免费观看三级毛片| 中文字幕乱码二三区免费|