• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    TC4鈦合金動態(tài)II型裂紋擴展速度與加載速度的關系研究

    2017-01-07 03:01:47苗潤王偉力宋楊
    兵工學報 2016年12期
    關鍵詞:尖端裂紋試件

    苗潤, 王偉力, 宋楊

    (1.海軍航空工程學院 兵器科學與技術系, 山東 煙臺 264001; 2.海軍裝備部, 北京 100071)

    TC4鈦合金動態(tài)II型裂紋擴展速度與加載速度的關系研究

    苗潤1, 王偉力1, 宋楊2

    (1.海軍航空工程學院 兵器科學與技術系, 山東 煙臺 264001; 2.海軍裝備部, 北京 100071)

    為了確定TC4鈦合金動態(tài)II型斷裂實驗中的起裂時間,計算TC4鈦合金在發(fā)生II型斷裂時塑性區(qū)能量吸收和裂紋擴展的速度,使用分離式霍普金森壓桿裝置對TC4鈦合金進行動態(tài)II型裂紋擴展實驗。采用高速數(shù)字散斑相關方法獲取觀測點處位移- 時間歷程曲線,確定了動態(tài)II型裂紋擴展的起裂時間,并計算出裂紋擴展速度。使用ABAQUS/Explicit軟件建立有限元模型與實驗結(jié)果進行同工況比對,在確認模型準確性后采用不同的加載速度研究II型加載條件下的裂紋擴展形式差異,計算出不同加載速度下的II型裂紋擴展速度,并得出加載速度與II型裂紋擴展速度的關系式。

    金屬材料; II型裂紋; 裂紋擴展速度; 分離式霍普金森壓桿; 高速數(shù)字散斑相關方法

    0 引言

    由于大量動態(tài)載荷問題的產(chǎn)生,人們發(fā)現(xiàn)在低應變率條件下的許多常用材料力學性能參數(shù)無法滿足實際工程與實踐需要,例如在高加載率條件下材料斷裂韌性的率敏感性[1]、絕熱剪切帶的產(chǎn)生[2]以及變形溫度對材料塑性變形行為的影響[3]等。使用分離式霍普金森壓桿(SHPB)裝置配合用電測法和數(shù)字散斑相關技術方法,可以對材料的基本力學性能進行計算和觀察,但由于實驗條件和實驗成本的原因,一些條件無法得到滿足:

    1)由于高速數(shù)字散斑相關方法(HS-DSCM)需架設高速攝像機,相機距離試樣很近,導致子彈加載速度不宜過快,否則容易損壞相機,無法實現(xiàn)更高應變率下的動態(tài)測試;

    2)試樣數(shù)目較少,數(shù)據(jù)采集成功率較低,無法采用更多組的樣本來進行加載從而驗證加載速度與裂紋擴展速度的關系;

    3)由于裂紋擴展范圍較大,當速度過快或疲勞裂紋過長時,裂紋擴展將會脫離高速攝影視場。無法對脫離視場的部分進行全面分析。

    本文通過數(shù)值模擬與實驗的對照組來驗證模型的準確性,再通過增加試樣的樣本容量,來研究加載速度與II型裂紋擴展的關系。

    1 動態(tài)II型裂紋擴展的實驗分析

    1.1 實驗材料

    材料選用TC4鈦合金,鈦合金具有密度低、比強度高、耐腐蝕性好等特點,已在航空工業(yè)、化工產(chǎn)業(yè)等領域得到廣泛應用[4]。

    TC4鈦合金是α+β兩相鈦合金,采用中溫變形進行制造加工[11](溫度區(qū)間為0.3~0.5Tm,Tm是材料熔點),由于TC4合金的導熱性能較差,易產(chǎn)生絕熱剪切帶,此現(xiàn)象已經(jīng)在Timothy等[5]、Rogers等[6]的實驗中得到證實。對TC4材料進行II型裂紋加載將比較容易產(chǎn)生熱塑失穩(wěn)現(xiàn)象,從而對絕熱剪切帶形成和擴展的研究產(chǎn)生影響。該實驗用TC4材料的成分組成與基本物理參數(shù)分別見表1和表2.

    1.2 實驗方案

    試樣制作成單邊裂紋的形式并含有線切割,其中試樣寬W=60 mm,高L=100 mm,厚度B=6.5 mm,在試樣的中部設置一條線切割裂紋,該裂紋上表面距離上端50 mm,線切割長度a=22 mm,寬度為0.3 mm,將該試樣設定為標準試樣。為實現(xiàn)實驗目的,除標準試樣外,在線切割裂紋的基礎上,采用高頻疲勞實驗機預制長度b分別為0.2 mm、0.3 mm、0.5 mm、1 mm、2 mm、3 mm的疲勞裂紋。試樣圖示和實際樣件如圖1所示。將試件固定于特制夾具上,如圖2所示。采用發(fā)射槍擊發(fā)橫截面直徑為16 mm的鋼制柱狀子彈撞擊橫截面直徑為16 mm、長為1 400 mm的入射桿,通過間接撞擊試件的加載方式來實現(xiàn)動態(tài)II型裂紋加載。實驗整體裝置如圖3所示。

    表1 TC4的成分組成Tab.1 Composition of TC4 %

    表2 TC4材料的基本物理參數(shù)Tab.2 Basic physical parameters of TC4

    圖1 試樣二維幾何圖樣Fig.1 2D geometric pattern of sample

    圖2 實驗現(xiàn)場照片F(xiàn)ig.2 Photo of experimental site

    圖3 實驗裝置圖Fig.3 Schematic diagram of experimental apparatus

    1.3 測試方法

    首先采用傳統(tǒng)電測法對裂紋尖端附近應變場進行測試,根據(jù)Kalthoff等[7]對于II型裂紋的失效方式研究中得到結(jié)論,測試裂紋尖端應變的應變片粘貼位置如圖4所示,選取φ≈60°,α≈10°,r=5 mm.

    圖4 試樣上應變片位置示意圖Fig.4 Position of strain gauge on sample

    但傳統(tǒng)電測法無法獲得形象的細節(jié)特征,對材料斷裂過程破壞模式的分析,需要采取后效觀察方式,因此,通過高速數(shù)字散斑相關技術達到的實時觀測,結(jié)合計算機軟件對于試件上噴涂的物質(zhì)點的追蹤,能夠較為準確地描述高應變率條件下試樣變形、裂紋擴展的全過程。該方法不需要均勻應力- 應變狀態(tài),也不需要力的平衡條件, 憑借著隨機噴在試件上黑- 白相間的物質(zhì)點,實時動態(tài)地記錄試件變形的每一個瞬間。通過校準板的標定,可準確地給出試件表面的應變場- 時間歷程,它給出的實驗結(jié)果說明均勻性假設是有一定局限性的。HS-DSCM測試技術的這一特色,使得Ⅱ型裂紋的起裂時間以及沖擊載荷下裂紋失穩(wěn)擴展的速度測量成為可能。

    為使散斑圖像清晰并節(jié)約計算成本,散斑大小設定約為300×300像素/in,試件測試面面積為60 cm2,散斑總數(shù)約為5 800左右。本實驗所使用的高速攝影設備德國DANTEC DYNAMICS公司的Q450,當發(fā)射槍擊發(fā)柱狀子彈時發(fā)出巨大聲響,通過聲控開關激發(fā)Q450運行,該設備可滿足每1 μs輸出一幅圖像。因設計總時長較長,為節(jié)約計算成本,在滿足追蹤裂紋擴展全過程的條件下,將輸出頻率改為每10 μs輸出一幅圖像。再配備ISTRA 4D軟件,對試件變形乃至破壞的全過程進行跟蹤,可直接給出觀測點x方向位移。觀測點的選取主要考慮兩個方面:首先本實驗后續(xù)希望計算得到II型裂紋動態(tài)斷裂韌性,因此需要選擇x方向位移占總位移主要成分的點,因此選擇過裂紋尖端平行于水平面的直線上的點;其次由于在個別區(qū)域物質(zhì)點噴涂不均勻,可能導致觀測區(qū)域數(shù)據(jù)獲取失敗,因此需要選擇噴涂較為均勻的點作為參考。ISTRA 4D軟件可判定觀測點周邊區(qū)域數(shù)據(jù)獲取是否完整,通過多次比對,設定觀測點距離裂紋尖端0.5 mm時獲取數(shù)據(jù)最為完整。圖5為0.2 mm預制裂紋試樣觀測區(qū)域和觀測點的位置圖。

    圖5 高速數(shù)字散斑相關方法觀測點位置Fig.5 Position of observation point by HS-DSCM

    1.4 實驗結(jié)果分析

    1.4.1 II型裂紋起裂時間的確定

    當試件承受載荷作用時,應變片上所測得的應變隨載荷的增大而增大,當裂紋起裂時,產(chǎn)生卸載波,使得應變劇烈減小。因此,應變片上所測應變的最大值對應的時間減去應力波從裂尖傳播至應變片所需要的時間就是起裂時間。

    計算試樣裂紋擴展速度和材料動態(tài)斷裂韌性都需要得到材料在II型加載時的起裂時間[8]。對于Ⅰ型裂紋,裂紋起裂的瞬間,在其尖端處,彈性能釋放使得應變迅速下降。由于沒有剪應力作用,應力波傳至裂紋尖端到裂紋擴展的過程異常迅速,因此對于起裂時間的計算也較為明顯。此方法在II裂紋研究中就出現(xiàn)了問題。如圖6所示,在應力波傳播至裂紋尖端到裂紋開始擴展的過程中,應變出現(xiàn)了幾次明顯的波動。僅僅通過電測法無法確定起裂時間的起點,可將由HS-DSCM獲取的觀測點(見圖5)的總位移- 時間曲線(見圖7)與高速攝影獲得的間斷點處實時照片進行對照,從而確定動態(tài)II型裂紋擴展過程中的起裂時間[9]。

    圖6 II型裂紋動態(tài)斷裂韌性實驗應變片信號Fig.6 Strain gauge signal of Mode II crack dynamic fracture toughness

    圖8(a)、圖8(b)、圖8(c)中每一步長為10 μs,當90 μs時,應力波傳播至裂紋尖端,而此時裂紋并未開裂,能量在塑性區(qū)累積,如圖8(a)所示;在之后30 μs的時間內(nèi),預制裂紋張開,此時裂紋尖端并未擴展,裂尖塑性區(qū)產(chǎn)生;120 μs時,預制裂紋的張開行為停止,動能迅速轉(zhuǎn)化為裂紋尖端的塑性能,積累時間約30 μs,如圖8(b)所示;當再次產(chǎn)生位移增量時即為起裂起點,持續(xù)90 μs,裂紋開始擴展,因此需要的起裂時間即為從150 μs至240 μs段,共計90 μs,如圖8(c)所示。將高速攝影的拍攝圖片與HS-DSCM得到的位移- 時間曲線(見圖7)進行對照得到了動態(tài)II型裂紋擴展從應力到達裂紋尖端,塑性區(qū)能量積累到裂紋失穩(wěn)擴展的全過程??梢源_定該試樣動態(tài)II型裂紋擴展起裂時間為90 μs.

    圖8 90 μs、150 μs、240 μs實況圖Fig.8 Crack initiations for 90 μs,150 μs and 240 μs

    1.4.2 II型裂紋擴展速度的確定

    通過金相照片(見圖9)得出裂紋尖端位移和裂紋擴展長度。

    圖9 試樣裂紋尖端金相照片(放大500倍)Fig.9 Metallograph of sample’s crack tip(500×)

    隨著裂紋擴展長度的增加,裂紋張開位移在減小。裂紋擴展停止時裂紋尖端張開位移為16.96 μm. 根據(jù)Kachanov等[10]提出的一維蠕變損傷理論,在較大載荷下,斷裂行為與小載荷下明顯不同,該實驗說明了這個現(xiàn)象:裂紋擴展并未沿金屬晶粒的邊緣擴展,而是發(fā)生了與其近乎垂直方向的穿晶脆性斷裂。在忽略變形的情況下,裂紋擴展長度L=1.03 mm,起裂時間ΔtL=90 μs. 擴展速度為

    (1)

    將裂紋擴展速度近似為勻速擴展,可以得到脆性斷裂的時間約為150 μs. 與Kachanov等[10]得到的有損傷無變形脆性斷裂時間吻合。計算公式為

    (2)

    式中:ν為材料泊松比;c為該材料波速;σ0為臨界應力。

    2 動態(tài)II型裂紋擴展的有限元對照實驗

    為了驗證有限元模型的準確性,首先設定與實驗相同的工況進行計算和對比。有限元建模過程中子彈和試樣的尺寸都完全按照具體實驗中加工的試件尺寸定義。子彈的橫截面直徑為16 mm,入射桿截面直徑為16 mm,長1.4 m,試樣寬60 mm,高100 mm,厚度為6.5 mm. 在試樣中部沿水平位置設置一條長22 mm、寬0.3 mm的線切割裂縫,為了更加準確地模擬線切割效果,在裂縫內(nèi)端部設置兩個0.05°的圓形倒角。預制疲勞裂紋的模擬是在線切割裂縫端部中間位置設置一條一定長度,開口10 μm的尖端裂紋。由于試樣模型符合平面應變條件,因此模擬中采用CPE4R和CPE3平面應變單元。裂紋區(qū)域是發(fā)生變形局部化的敏感區(qū)域,所以該區(qū)域的單元劃分必須非常細密。為了使該區(qū)域劃分的網(wǎng)格比較規(guī)整,首先通過分區(qū)將試樣預制裂紋部分分離出來,然后采用一個過線切割端部倒角兩端點的橢圓將模型的裂紋部分和其他區(qū)域分割開來進行網(wǎng)格加密,橢圓的長軸為4.15 mm,短軸為2.64 mm. 裂紋區(qū)域和遠離裂紋的兩邊部分采用結(jié)構(gòu)化網(wǎng)格劃分方法,其間的過渡區(qū)域采用自由網(wǎng)格劃分方法。如果計算中網(wǎng)格劃分過粗則不能真實地反映裂紋起裂過程,而網(wǎng)格劃分過細,則增大計算成本,造成不必要的資源浪費。因此選擇裂紋周邊及裂紋面方向的單元最小尺寸為20 μm×40 μm,遠離裂紋面處網(wǎng)格尺寸為1 mm. 以預制裂紋長度為0.2 mm的模型為例,網(wǎng)格劃分完成后,試樣上的節(jié)點數(shù)目為46 450個,單元數(shù)目為46 738個。網(wǎng)格劃分后試樣的模型圖和局部放大圖如圖10所示。

    圖10 網(wǎng)格劃分局部放大圖Fig.10 Partially enlarged view of sample after mesh generation

    材料采用Johnson-Cook熱黏塑性本構(gòu)模型、Mie-Gruneisen狀態(tài)方程并通過Johnson-Cook失效模型定義材料失效。

    通過紅外測速裝置測得該組實驗子彈的加載速度為34.88 m/s,采用相同的速度進行數(shù)值模擬,最終加載效果對比如圖11所示。

    圖11 實驗與數(shù)值模擬結(jié)果對比Fig.11 Comparison of experimental and simulated results

    預制裂紋長度為0.2 mm,在模擬中直接輸入速度信號,分別選取90 μs、150 μs、300 μs的應變場圖像,這3個時間點分別代表應力波傳播至裂紋尖端、裂紋開裂瞬間、裂紋開裂后。將其與數(shù)字散斑相關方法得到的高速攝影圖進行對比,如圖12所示。

    圖12 子彈打擊2 mm預制裂紋試樣時的有限元分析與高速攝影實時對比Fig.12 Real-time comparison of high-speed photography method and finite element method when bullet strikes the sample with 2 mm pre-production crack

    對數(shù)值模擬所得到的裂紋擴展速度進行統(tǒng)計。裂紋擴展速度隨時間的變化曲線如圖13所示。

    圖13 有限元法裂紋擴展速度時間曲線Fig.13 FEM crack velocity-time curve

    根據(jù)圖13得到的速度- 時間曲線,去除干擾點,將得到的裂紋擴展速度求平均值,得到打擊速度為34.88 m/s時裂紋擴展速度為12.74 m/s,實驗得到裂紋擴展速度為11.45 m/s.

    表3 不同預制裂紋長度的實驗數(shù)據(jù)Tab.3 Experimental data about different precrack lengths

    通過以上各組實驗與數(shù)值模擬對比,計算誤差較小,說明數(shù)值模擬方法計算準確,模型設計合理。

    3 不同加載速度下的數(shù)值模擬

    數(shù)值模擬可以在更廣泛的速度范圍內(nèi),選取更多裂紋種類進行Ⅱ型裂紋擴展的實驗模擬,這樣可以更詳細地了解試件的變形破壞情況。選取6種不同長度的預制疲勞裂紋進行了計算,6種裂紋長度分別為0.2 mm、0.3 mm、0.5 mm、1.0 mm和2.0 mm和3.0 mm,對每種長度的裂紋以不同的加載速度進行沖擊加載,加載速度范圍為20 m/s到100 m/s. 以0.2 mm裂紋為例。在加載到60 m/s時斷裂方式發(fā)生明顯變化。選擇加載速度從20~70 m/s的相應塑性應變分布如圖14所示。

    圖14 不同加載速度下的Ⅱ型裂紋試樣最終變形及相應塑性應變分布Fig.14 Final deformation of Mode II crack specimen and plastic strain distribution at different loading speeds

    Ⅱ型裂紋試樣在沖擊加載狀態(tài)下,塑性應變逐漸集中到裂紋擴展區(qū)域,并在裂尖處開始發(fā)生破壞。計算過程顯示,當加載速度低于12 m/s時,裂尖處只發(fā)生輕微的塑性變形,并沒有起裂;加載速度為12 m/s時,裂尖處單元開始發(fā)生失效,失效單元被刪除,裂紋開始沿裂紋面方向擴展。當加載速度低于50 m/s時,裂紋首先沿水平方向迅速傳播,然后傳播速度逐漸下降,之后裂紋沿著右上方約45°方向失穩(wěn)擴展。這是由于子彈和試樣作用的初期能夠滿足Ⅱ型裂紋擴展的條件,而隨著子彈和試樣的相互作用,試樣的上半部分發(fā)生傾斜,子彈右端和試樣左表面的接觸不再滿足Ⅱ型裂紋的擴展條件,此時,裂紋發(fā)展成為混合型裂紋。隨著加載速度不斷增大,裂紋沿著水平方向擴展的距離也不斷增大;當加載速度達到60 m/s時,整個試樣被裂紋截斷為兩部分。即加載速度越高,裂紋最終的擴展長度越大,試樣的形變率越高??梢钥闯觯敿虞d速度高于60 m/s時,整個計算過程完全滿足純剪切條件,裂紋擴展形式發(fā)生明顯變化。下面就以60 m/s打擊速度的數(shù)值模擬為例進行分析。

    取裂紋擴展區(qū)域單元45 795和單元45 813的應力- 時間曲線和應變- 時間曲線如圖15所示。

    圖15 60 m/s加載速度下的等效應力- 時間與等效應變- 時間曲線Fig.15 Equivalent stress-time curves and equivalent strain-time curves at 60 m/s

    由圖15可以看出,加載速度為60 m/s時,單元45 795和單元45 813的應力在下降階段都下降到0,說明裂紋擴展區(qū)域上的這兩個單元都發(fā)生了失效。由等效應變- 時間曲線可以看出,單元45 795的等效應變首先達到最大失效應變,大約16 μs時單元45 813才達到最大失效應變,這說明裂紋由裂尖處開始沿裂紋擴展區(qū)域傳播,擴展區(qū)域沿線的單元是先后發(fā)生實效的。在這一加載速度下最后試樣被裂紋截為兩半,試樣完全破壞。加載速度為60 m/s時,不同時刻裂紋擴展區(qū)域的等效應力云圖如圖16所示。

    圖16 加載速度為60 m/s時不同時刻的等效應力云圖Fig.16 Equivalent stress contour plots at different moments at 60 m/s

    通過數(shù)值計算可以得出60 m/s加載速度下,裂紋擴展速度對時間的變化曲線如圖17所示。

    圖17 初始裂紋為0.2 mm時60 m/s加載速度下裂紋傳播速度- 時間曲線Fig.17 Crack propagation speed-time curve at 60 m/s whith crack of 0.2 mm in length

    對于同樣裂紋長度的試樣,計算不同打擊速度下裂紋擴展速度可得到圖18的曲線圖,通過對點信息進行擬合,可以得出加載速度與裂紋擴展速度的近似方程為

    y=ae-e-k(x-xc),

    (3)

    式中:擬合出的系數(shù)a=778.05,xc=70.88,k=0.12.

    將每一種打擊速度下的裂紋擴展速度取平均值,見表4,匯總?cè)鐖D19所示。

    圖18 裂紋平均傳播速度- 加載速度關系擬合曲線Fig.18 Fitted curve of average crack propagation speed and loading speed

    圖19 不同裂紋長度對應不同加載速度時的裂紋擴展速度Fig.19 Crack propagation speeds at different loading speeds relative to different crack lengths

    4 結(jié)論

    1)對TC4材料進行動態(tài)Ⅱ型裂紋擴展實驗,采用電測法和數(shù)字散斑相關方法,結(jié)合金相照片,確定了起裂時間,得到在高加載速率下TC4材料的裂紋擴度速度等重要參數(shù)。

    2)通過對照實驗驗證了有限元模型的準確性,從而可以通過數(shù)值模擬方法擴充不同加載速度下的實驗樣本。發(fā)現(xiàn)當加載速度達到一定值時,材料的破壞形式將發(fā)生變化,從向斜上方約45°的裂紋擴展形式轉(zhuǎn)變?yōu)檠厮椒较驍U展的破壞形式。

    表4 不同裂紋長度對應不同加載速度時的裂紋擴展速度表Tab.4 Crack propagation speeds at different loading speeds relative to different crack lengths

    3)通過不同加載速度下裂紋擴展速度的計算,得到了裂紋擴展速度與加載速度的曲線,并擬合出其對應關系,得出關系方程。該方程可在得到某種材料某一加載速度下的裂紋擴展速度后,預測其他加載速度下的裂紋擴展速度,可減少計算時間,提高計算效率。

    References)

    [1] 許澤建, 李玉龍, 劉元庸, 等. 兩種高強鋼在高加載速率下的II型動態(tài)斷裂韌性[J]. 金屬學報, 2006, 42(6): 635-640. XU Ze-jian, LI Yu-long, LIU Yuan-yong, et al. Mode II dynamic fracture toughness of two high strength steels under high loading rate[J]. Acta Metallurgica Sinica, 2006, 42(6): 635-640. (in Chinese)

    [2] 蘇冠龍, 龔煦, 李玉龍, 等. TC4在動態(tài)載荷下的剪切行為研究[J]. 爆炸與沖擊, 2015, 35(4): 527-529. SU Guan-long, GONG Xu, LI Yu-long, et al. Shear behavior of TC4 alloy under dynamic loading[J]. Explosion and Shock Waves, 2015, 35(4): 527-529. (in Chinese)

    [3] 楊新亮, 葛長健, 楊智. 變形溫度對TC4合金中溫塑性變形行為的影響[J]. 熱加工工藝, 2012, 41(16): 71-73,76. YANG Xin-liang, GE Chang-jian, YANG Zhi.The influence of stress state on the development of adiabatic shear for TC4 alloys[J]. Hot Working Technology, 2012, 41(16): 71-73,76. (in Chinese)

    [4] 張長清, 謝蘭生, 陳明和, 等. 高應變率下TC4-DT鈦合金的動態(tài)力學性能與塑性本構(gòu)關系[J]. 中國有色金屬學報. 2015, 25(2): 323-324. ZHANG Chang-qing, XIE Lan-sheng, CHEN Ming-he, et al. Dynamic mechanical property and plastic constitutive relation of TC4-DT Ti alloy under high strain rate[J] The Chinese Journal of Nonferrous, 2015, 25(2): 323-324. (in Chinese)

    [5] Timothy S P, Hutchings I M. The structrue of adiabatic shear band in a titanium alloy[J]. Acta Metallurgica, 1985, 33(4): 667-676.

    [6] Roger H C. Adiabatic plastic deformation[J]. Annual Reciew of Materials Science, 1979, 9(1):283-311.

    [7] Kalthoff J F, Winkler S. Modes of dynamic shear failure in solids[J].International Journal of Fracture, 2000, 101(1):1-31.

    [8] 董新龍, 虞吉林, 胡時勝, 等. 高加載率下Ⅱ型裂紋試樣的動態(tài)應力強度因子及斷裂行為[J]. 爆炸與沖擊, 1998, 18(1): 62-68. DONG Xin-long, YU Ji-lin, HU Shi-sheng, et al. Stress intensity factor and fracture behavior for mode II crack specimen under high shear loading rate[J].Explosion and Shock Waves,1998, 18(1): 62-68. (in Chinese)

    [9] 苗潤, 范亞夫, 王偉力. 動態(tài)II型加載條件下兩種金屬材料起裂時間判定方法的實驗研究[C]∥第十四屆戰(zhàn)斗部與毀傷技術學術交流會. 重慶: 中國宇航學會無人飛行器分會, 2015: 1080-1081. MIAO Run, FAN Ya-fu, WANG Wei-li. Research on the judgement method of crack time about two metal mertials in dynamic mode II loading condition[C]∥The Fourth Symposium on Warheads and Damcege Technology. Chongqing: Branch of Unmanned Air Vehicles, Chinese Society of Astronautics, 2015: 1080-1081.(in Chinese)

    [10] Kachanov L M, Rabotnov Y N. Time of the rupture process under creep conditons[J]. Izvestiya Akademii Nauk SSSR Otdelenie Tekniches,1958,8: 26-31.

    [11] 郝露菡, 叢陽陽, 彭艷. 考慮相變的TC4鈦合金流動應力研究[J]. 材料科學與工藝, 2015, 23(5): 12-14. HAO Lu-han, CONG Yang-yang, PENG Yan. Research of flow stress of TC4 titanium considering phase change[J]. Materials Science & Technology, 2015, 23(5): 12-14. (in Chinese)

    Research on Relationship between Crack Propagation Speed of TC4 Alloy and Loading Speed under Mode II Dynamic Loading

    MIAO Run1, WANG Wei-li1, SONG Yang2

    (1.Department of ordnance Science and Technology, Naval Aeronautical and Astronautical University, Yantai 264001,Shandong, China; 2.Naval Equipment Department, Beijing 100071, China)

    In order to ascertain the crack initiation time of TC4 in Mode II crack experiment, the energy absorption in plastic zone and the crack propagation speed are calculated. Split Hopkinson pressure bar (SHPB) is used to experiment upon Mode II crack propagation of TC4. High-speed digital speckle correlation method (HS-DSCM) is used to get displacement-time curve at a observation point, the Mode II crack initiation time of TC4 is determined, and the propagation speed of Mode II crack is calculated. A finite element model is established using ABAQUS/Explicit. The simulated results is compared to the experimental results. The propagation forms of crack at different impact velocities under loading condition are researched. Mode II crack propagation speeds at different loading speeds are calculated. Relation between loading speed and Mode II crack propagation speed is obtained.

    metallic material; Mode II crack; crack propagation speed; split Hopkinson pressure bar; HS-DSCM

    2016-05-05

    苗潤(1989—),男,博士研究生。E-mail: miaorun1769@163.com; 王偉力(1962—),男,教授,博士生導師。E-mail: w.l.wang@tom.com

    10.3969/j.issn.1000-1093.2016.12.020

    O346.1+2; O347.3

    A

    1000-1093(2016)12-2331-09

    猜你喜歡
    尖端裂紋試件
    復材管纖維纏繞角度對約束混凝土軸壓性能的影響研究
    廣東建材(2022年5期)2022-06-10 08:48:20
    Epidermal growth factor receptor rs17337023 polymorphism in hypertensive gestational diabetic women: A pilot study
    不同參數(shù)對開圓洞鋼板剪力墻抗震性能的影響
    科學中國人(2018年8期)2018-07-23 02:26:56
    微裂紋區(qū)對主裂紋擴展的影響
    國際標準攻堅克難“S試件”美玉漸成
    鏡頭看展
    預裂紋混凝土拉壓疲勞荷載下裂紋擴展速率
    基于位移相關法的重復壓裂裂縫尖端應力場研究
    斷塊油氣田(2014年5期)2014-03-11 15:33:49
    加速尖端機床國產(chǎn)化
    亚洲av熟女| av电影中文网址| 久久久久国产精品人妻aⅴ院| 亚洲中文字幕日韩| 亚洲激情在线av| 国产av精品麻豆| a级毛片在线看网站| 国产区一区二久久| 波多野结衣巨乳人妻| 日韩中文字幕欧美一区二区| 一进一出好大好爽视频| 在线观看午夜福利视频| 欧美另类亚洲清纯唯美| 熟妇人妻久久中文字幕3abv| 神马国产精品三级电影在线观看 | 亚洲少妇的诱惑av| 国产成人精品在线电影| 国产av又大| 色尼玛亚洲综合影院| 国产精品乱码一区二三区的特点 | 午夜久久久在线观看| 国产熟女xx| 999久久久精品免费观看国产| 久久狼人影院| 国产成人精品久久二区二区91| 国产亚洲欧美98| 国产激情欧美一区二区| 18禁观看日本| 国产午夜精品久久久久久| 视频在线观看一区二区三区| √禁漫天堂资源中文www| 免费在线观看视频国产中文字幕亚洲| 国产野战对白在线观看| 亚洲欧美一区二区三区黑人| 国产亚洲欧美精品永久| 777久久人妻少妇嫩草av网站| 亚洲国产毛片av蜜桃av| 高清黄色对白视频在线免费看| 99久久精品国产亚洲精品| 成人欧美大片| ponron亚洲| 99国产精品免费福利视频| 18禁裸乳无遮挡免费网站照片 | 欧美最黄视频在线播放免费| 日本免费一区二区三区高清不卡 | 一本综合久久免费| 美女大奶头视频| 日韩视频一区二区在线观看| 女人被狂操c到高潮| 亚洲国产欧美网| 午夜精品国产一区二区电影| 国产三级黄色录像| 国内精品久久久久久久电影| 国产熟女xx| 999久久久国产精品视频| 91精品三级在线观看| 欧美成人午夜精品| 美女免费视频网站| 国产精品,欧美在线| 88av欧美| 国产私拍福利视频在线观看| АⅤ资源中文在线天堂| 十八禁网站免费在线| 日本精品一区二区三区蜜桃| 欧美黑人欧美精品刺激| 国产91精品成人一区二区三区| 岛国视频午夜一区免费看| 国产精品自产拍在线观看55亚洲| 亚洲av熟女| 好男人在线观看高清免费视频 | 99久久综合精品五月天人人| 天堂影院成人在线观看| 国产成人免费无遮挡视频| 亚洲最大成人中文| 国产亚洲av高清不卡| 日韩精品青青久久久久久| 国产精品香港三级国产av潘金莲| 国产一区二区三区综合在线观看| 无人区码免费观看不卡| 制服人妻中文乱码| 国产精品国产高清国产av| 91九色精品人成在线观看| 成熟少妇高潮喷水视频| xxx96com| 老司机靠b影院| 人人妻人人爽人人添夜夜欢视频| 久99久视频精品免费| 人妻丰满熟妇av一区二区三区| 日本黄色视频三级网站网址| 日本精品一区二区三区蜜桃| 人人妻人人澡人人看| 国产1区2区3区精品| 国产午夜福利久久久久久| 宅男免费午夜| 欧美人与性动交α欧美精品济南到| 亚洲欧美日韩高清在线视频| 12—13女人毛片做爰片一| 十八禁网站免费在线| 欧美日韩福利视频一区二区| 一区二区三区国产精品乱码| 91字幕亚洲| 美女高潮到喷水免费观看| 中文字幕人妻熟女乱码| 看免费av毛片| 精品人妻1区二区| 大陆偷拍与自拍| 极品人妻少妇av视频| 亚洲国产日韩欧美精品在线观看 | 亚洲avbb在线观看| 国产精品98久久久久久宅男小说| 99国产精品免费福利视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品国产综合久久久| 免费看a级黄色片| 少妇裸体淫交视频免费看高清 | 母亲3免费完整高清在线观看| 亚洲国产精品久久男人天堂| 国产一级毛片七仙女欲春2 | 桃色一区二区三区在线观看| 老司机在亚洲福利影院| 中文字幕人妻丝袜一区二区| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 日韩精品免费视频一区二区三区| 中国美女看黄片| 最近最新中文字幕大全电影3 | 男女午夜视频在线观看| 国产高清有码在线观看视频 | www.www免费av| 久久热在线av| 亚洲精品一区av在线观看| 制服丝袜大香蕉在线| 99久久综合精品五月天人人| 变态另类成人亚洲欧美熟女 | 韩国av一区二区三区四区| 国产三级在线视频| 一个人观看的视频www高清免费观看 | 91精品国产国语对白视频| 欧美+亚洲+日韩+国产| 国产精品精品国产色婷婷| 国产单亲对白刺激| 激情在线观看视频在线高清| 国产黄a三级三级三级人| 国产成人一区二区三区免费视频网站| 激情在线观看视频在线高清| 国产蜜桃级精品一区二区三区| 青草久久国产| 国产亚洲精品av在线| 亚洲国产毛片av蜜桃av| 欧美成狂野欧美在线观看| 国产精品野战在线观看| 国产真人三级小视频在线观看| 电影成人av| 亚洲在线自拍视频| 涩涩av久久男人的天堂| 国产成人影院久久av| 国内精品久久久久精免费| 久久人人97超碰香蕉20202| 久久中文字幕人妻熟女| 久久人妻av系列| 午夜视频精品福利| 亚洲av熟女| xxx96com| 国产精品久久久av美女十八| 成人手机av| 日日摸夜夜添夜夜添小说| 亚洲成人精品中文字幕电影| 国产精品久久视频播放| 热re99久久国产66热| 亚洲国产毛片av蜜桃av| 岛国在线观看网站| 亚洲精品一区av在线观看| 国产高清有码在线观看视频 | www.www免费av| 日日摸夜夜添夜夜添小说| 亚洲欧美日韩无卡精品| av超薄肉色丝袜交足视频| 9色porny在线观看| 久久国产乱子伦精品免费另类| 免费人成视频x8x8入口观看| 国产成人免费无遮挡视频| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频 | 国产成人系列免费观看| 18禁裸乳无遮挡免费网站照片 | 免费无遮挡裸体视频| 国产精品一区二区三区四区久久 | 色尼玛亚洲综合影院| 女人被躁到高潮嗷嗷叫费观| 一级片免费观看大全| 一级毛片女人18水好多| 午夜影院日韩av| 国产精品久久视频播放| 美国免费a级毛片| 婷婷六月久久综合丁香| 国产麻豆69| 亚洲国产精品合色在线| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| 神马国产精品三级电影在线观看 | av在线天堂中文字幕| 成年女人毛片免费观看观看9| 免费观看精品视频网站| 搡老熟女国产l中国老女人| 亚洲人成网站在线播放欧美日韩| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| 不卡一级毛片| 免费看美女性在线毛片视频| 亚洲欧美日韩高清在线视频| 一a级毛片在线观看| 亚洲成av人片免费观看| 少妇熟女aⅴ在线视频| 中文字幕人妻丝袜一区二区| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| 91av网站免费观看| 黄色毛片三级朝国网站| 黄片小视频在线播放| 非洲黑人性xxxx精品又粗又长| 色老头精品视频在线观看| 免费高清视频大片| 亚洲国产中文字幕在线视频| 涩涩av久久男人的天堂| 欧美成狂野欧美在线观看| 熟妇人妻久久中文字幕3abv| 美女高潮到喷水免费观看| 亚洲精品国产一区二区精华液| 色哟哟哟哟哟哟| 久久人妻熟女aⅴ| 两个人看的免费小视频| 欧美在线黄色| 美女 人体艺术 gogo| 老司机在亚洲福利影院| 日韩欧美免费精品| 窝窝影院91人妻| 亚洲片人在线观看| 亚洲免费av在线视频| 国产精品亚洲一级av第二区| 久久久精品欧美日韩精品| 欧美中文综合在线视频| 美国免费a级毛片| av片东京热男人的天堂| 国产日韩一区二区三区精品不卡| 欧美性长视频在线观看| 欧美av亚洲av综合av国产av| 免费女性裸体啪啪无遮挡网站| 亚洲成a人片在线一区二区| 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 精品国产乱码久久久久久男人| 国产精华一区二区三区| 国产片内射在线| 亚洲人成网站在线播放欧美日韩| 免费观看人在逋| 国产精品一区二区免费欧美| 男女下面进入的视频免费午夜 | 99riav亚洲国产免费| 一级毛片精品| 国产欧美日韩一区二区三| e午夜精品久久久久久久| 成人特级黄色片久久久久久久| 91av网站免费观看| 久久人人爽av亚洲精品天堂| 久热这里只有精品99| 极品人妻少妇av视频| 日韩中文字幕欧美一区二区| 免费搜索国产男女视频| 成人精品一区二区免费| 亚洲精品美女久久久久99蜜臀| 成年女人毛片免费观看观看9| 丝袜在线中文字幕| 日韩高清综合在线| 午夜福利影视在线免费观看| 美女高潮到喷水免费观看| 两人在一起打扑克的视频| 日韩国内少妇激情av| 国产一区在线观看成人免费| 精品免费久久久久久久清纯| 性欧美人与动物交配| 亚洲国产精品合色在线| 亚洲av日韩精品久久久久久密| 国产精品永久免费网站| 99re在线观看精品视频| bbb黄色大片| 欧美成人午夜精品| 一级片免费观看大全| 91在线观看av| 亚洲熟妇熟女久久| 最近最新中文字幕大全电影3 | 国产高清激情床上av| 亚洲人成伊人成综合网2020| 黄色 视频免费看| 啪啪无遮挡十八禁网站| 好看av亚洲va欧美ⅴa在| 亚洲av美国av| 亚洲成国产人片在线观看| 亚洲国产欧美日韩在线播放| 一边摸一边抽搐一进一小说| 亚洲国产精品合色在线| 免费看美女性在线毛片视频| 久久 成人 亚洲| 亚洲第一青青草原| 国产成人啪精品午夜网站| 亚洲av成人不卡在线观看播放网| 麻豆国产av国片精品| 美女国产高潮福利片在线看| 久久人妻福利社区极品人妻图片| 免费一级毛片在线播放高清视频 | 一级毛片精品| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品久久男人天堂| 国产又色又爽无遮挡免费看| 久久这里只有精品19| 夜夜看夜夜爽夜夜摸| 国内毛片毛片毛片毛片毛片| 亚洲中文日韩欧美视频| 99精品欧美一区二区三区四区| 人妻久久中文字幕网| 91字幕亚洲| 欧美成狂野欧美在线观看| 国语自产精品视频在线第100页| 成人亚洲精品一区在线观看| 国产精品电影一区二区三区| 久久天堂一区二区三区四区| 亚洲精品国产色婷婷电影| 亚洲第一av免费看| 首页视频小说图片口味搜索| 精品人妻1区二区| 欧美色视频一区免费| 一二三四在线观看免费中文在| 免费观看人在逋| 麻豆一二三区av精品| 欧美午夜高清在线| 午夜成年电影在线免费观看| 亚洲自拍偷在线| 国产成人一区二区三区免费视频网站| 精品电影一区二区在线| 久久久国产成人免费| 欧美成人午夜精品| 亚洲精华国产精华精| 又黄又粗又硬又大视频| 免费观看人在逋| 色播在线永久视频| 欧美日韩中文字幕国产精品一区二区三区 | 男男h啪啪无遮挡| 亚洲精品美女久久久久99蜜臀| 亚洲欧洲精品一区二区精品久久久| 欧美av亚洲av综合av国产av| 久久久国产成人精品二区| 日本免费一区二区三区高清不卡 | 黄色 视频免费看| 中国美女看黄片| 精品日产1卡2卡| 97碰自拍视频| 久久这里只有精品19| 国产高清视频在线播放一区| 老鸭窝网址在线观看| 久久国产精品人妻蜜桃| 日本vs欧美在线观看视频| 成熟少妇高潮喷水视频| av电影中文网址| 国产在线观看jvid| 国产成人影院久久av| 丝袜美足系列| 淫秽高清视频在线观看| 午夜福利影视在线免费观看| 最近最新中文字幕大全免费视频| 又黄又粗又硬又大视频| 日本a在线网址| 久久精品亚洲熟妇少妇任你| 男女之事视频高清在线观看| 精品福利观看| 亚洲精品久久国产高清桃花| 亚洲中文字幕日韩| 欧美人与性动交α欧美精品济南到| 黄色 视频免费看| 午夜精品久久久久久毛片777| 身体一侧抽搐| 国产欧美日韩综合在线一区二区| 日本五十路高清| 99久久精品国产亚洲精品| 最好的美女福利视频网| 可以免费在线观看a视频的电影网站| 国产精品久久久av美女十八| 欧洲精品卡2卡3卡4卡5卡区| 亚洲视频免费观看视频| 人妻丰满熟妇av一区二区三区| 在线天堂中文资源库| 高潮久久久久久久久久久不卡| 99精品欧美一区二区三区四区| 高潮久久久久久久久久久不卡| 99在线人妻在线中文字幕| 中文字幕久久专区| 国产精品爽爽va在线观看网站 | 国产精品电影一区二区三区| 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 12—13女人毛片做爰片一| bbb黄色大片| 韩国精品一区二区三区| 久久草成人影院| 日韩 欧美 亚洲 中文字幕| 亚洲国产中文字幕在线视频| 久久中文看片网| 精品国产乱子伦一区二区三区| 国产成年人精品一区二区| 欧美最黄视频在线播放免费| 成在线人永久免费视频| 国产精品乱码一区二三区的特点 | 国产高清视频在线播放一区| 免费在线观看亚洲国产| 国产精品久久久av美女十八| 在线播放国产精品三级| 亚洲精品美女久久av网站| 黄色毛片三级朝国网站| 香蕉丝袜av| 欧美国产精品va在线观看不卡| 男女下面插进去视频免费观看| 免费无遮挡裸体视频| 日韩高清综合在线| 午夜精品国产一区二区电影| 国产精品亚洲av一区麻豆| 国产又色又爽无遮挡免费看| 国产日韩一区二区三区精品不卡| 国产精品日韩av在线免费观看 | 夜夜爽天天搞| 亚洲国产高清在线一区二区三 | 亚洲熟妇中文字幕五十中出| tocl精华| 久久欧美精品欧美久久欧美| 久久香蕉激情| 波多野结衣高清无吗| 欧美人与性动交α欧美精品济南到| 此物有八面人人有两片| 三级毛片av免费| 欧美av亚洲av综合av国产av| 成年女人毛片免费观看观看9| 久久精品aⅴ一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 亚洲男人天堂网一区| 午夜老司机福利片| 国产野战对白在线观看| 日本在线视频免费播放| 欧美日本中文国产一区发布| 亚洲精品av麻豆狂野| 叶爱在线成人免费视频播放| 国产精品免费视频内射| 久久天躁狠狠躁夜夜2o2o| 亚洲国产高清在线一区二区三 | 久久久久久久久久久久大奶| 亚洲欧美激情在线| 嫁个100分男人电影在线观看| 成人av一区二区三区在线看| 中文字幕人妻熟女乱码| 国产精品久久久av美女十八| 少妇的丰满在线观看| av视频免费观看在线观看| 欧美日韩亚洲国产一区二区在线观看| 91麻豆精品激情在线观看国产| 久久精品国产99精品国产亚洲性色 | 国产亚洲精品第一综合不卡| 亚洲国产中文字幕在线视频| 亚洲国产看品久久| 午夜免费鲁丝| 欧美日韩乱码在线| 女性被躁到高潮视频| 桃红色精品国产亚洲av| 国产亚洲精品久久久久久毛片| 在线观看日韩欧美| 午夜精品国产一区二区电影| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 可以在线观看的亚洲视频| 日韩精品青青久久久久久| aaaaa片日本免费| 欧美黑人精品巨大| www.熟女人妻精品国产| 久久精品影院6| 可以在线观看毛片的网站| 久久久精品欧美日韩精品| 国产熟女xx| 1024视频免费在线观看| 在线观看免费视频日本深夜| 欧美在线一区亚洲| 成人av一区二区三区在线看| 免费在线观看黄色视频的| 日本欧美视频一区| 又黄又爽又免费观看的视频| 俄罗斯特黄特色一大片| 久久精品亚洲熟妇少妇任你| 欧美激情高清一区二区三区| 亚洲五月色婷婷综合| 又大又爽又粗| 此物有八面人人有两片| 成人国语在线视频| 午夜福利,免费看| 高潮久久久久久久久久久不卡| 俄罗斯特黄特色一大片| 欧美精品啪啪一区二区三区| 禁无遮挡网站| 国产欧美日韩一区二区三区在线| 精品久久蜜臀av无| 精品人妻在线不人妻| 精品国内亚洲2022精品成人| 欧美日韩中文字幕国产精品一区二区三区 | 一边摸一边做爽爽视频免费| 国产av又大| 欧美av亚洲av综合av国产av| 亚洲国产中文字幕在线视频| 国产精品九九99| 一区二区三区国产精品乱码| 母亲3免费完整高清在线观看| 国产亚洲av高清不卡| 免费在线观看黄色视频的| 欧美日本亚洲视频在线播放| 久久热在线av| av有码第一页| av免费在线观看网站| 亚洲人成电影观看| 免费看a级黄色片| 大码成人一级视频| 老汉色av国产亚洲站长工具| 一级毛片高清免费大全| 久久亚洲真实| 国产午夜福利久久久久久| 女人被狂操c到高潮| 午夜a级毛片| 国内精品久久久久精免费| 琪琪午夜伦伦电影理论片6080| 中文字幕高清在线视频| av超薄肉色丝袜交足视频| 激情视频va一区二区三区| 女人被狂操c到高潮| av在线播放免费不卡| 国产成人免费无遮挡视频| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 少妇粗大呻吟视频| 亚洲avbb在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 在线十欧美十亚洲十日本专区| 国产精品一区二区三区四区久久 | 国产精品爽爽va在线观看网站 | 亚洲av片天天在线观看| 午夜亚洲福利在线播放| 久久精品亚洲精品国产色婷小说| 欧美乱妇无乱码| 身体一侧抽搐| 大码成人一级视频| 国产91精品成人一区二区三区| 国产亚洲精品一区二区www| 国产一区二区在线av高清观看| 又黄又粗又硬又大视频| 久久午夜亚洲精品久久| 午夜影院日韩av| 欧美一级a爱片免费观看看 | 男女下面进入的视频免费午夜 | 精品国内亚洲2022精品成人| or卡值多少钱| 精品国产亚洲在线| 久久精品国产99精品国产亚洲性色 | 一区二区日韩欧美中文字幕| 一级毛片高清免费大全| 欧美日韩乱码在线| 国产高清视频在线播放一区| 首页视频小说图片口味搜索| 国产欧美日韩精品亚洲av| 国产亚洲精品av在线| 青草久久国产| 欧美+亚洲+日韩+国产| 中文字幕最新亚洲高清| 国产欧美日韩精品亚洲av| 搡老熟女国产l中国老女人| 香蕉国产在线看| 97碰自拍视频| xxx96com| 黑丝袜美女国产一区| 黄网站色视频无遮挡免费观看| 一进一出抽搐动态| 欧美老熟妇乱子伦牲交| 丁香六月欧美| 1024视频免费在线观看| 天天躁夜夜躁狠狠躁躁| 久久久久精品国产欧美久久久| 午夜福利,免费看| 色老头精品视频在线观看| 人成视频在线观看免费观看| av天堂在线播放| 国产单亲对白刺激| 国产一区二区在线av高清观看| 少妇 在线观看| 欧美日韩福利视频一区二区| 欧美午夜高清在线| 中文字幕av电影在线播放| 日韩高清综合在线| 好男人电影高清在线观看| 久久人妻av系列| 亚洲狠狠婷婷综合久久图片| 不卡一级毛片| 操美女的视频在线观看| АⅤ资源中文在线天堂| 精品久久久精品久久久| 别揉我奶头~嗯~啊~动态视频| 在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 在线播放国产精品三级| 日本免费a在线| 女人被躁到高潮嗷嗷叫费观| 国产成人av教育| 国产蜜桃级精品一区二区三区| 欧美绝顶高潮抽搐喷水|