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

    爆炸與電磁加載下無氧銅環(huán)、柱殼的斷裂模式轉(zhuǎn)變*

    2017-12-21 11:15:56郭昭亮劉明濤任國武湯鐵鋼劉倉理
    爆炸與沖擊 2017年6期
    關(guān)鍵詞:柱殼局域收縮率

    郭昭亮,范 誠,劉明濤,任國武,湯鐵鋼,劉倉理

    爆炸與電磁加載下無氧銅環(huán)、柱殼的斷裂模式轉(zhuǎn)變*

    郭昭亮,范 誠,劉明濤,任國武,湯鐵鋼,劉倉理

    (中國工程物理研究院流體物理研究所,四川 綿陽621999)

    考慮斷面收縮率、局域斷裂應(yīng)變以及平均斷裂應(yīng)變,并基于電磁膨脹環(huán)、爆炸膨脹環(huán)(柱殼)實(shí)驗(yàn)平臺,研究了高純無氧銅(TU1)環(huán)及柱殼在高應(yīng)變率載荷下的膨脹斷裂行為。采用高速攝影技術(shù)拍攝柱殼外壁的膨脹斷裂形貌演化過程,用于確定柱殼平均斷裂應(yīng)變;利用激光干涉測速技術(shù)獲得樣品徑向膨脹速度歷史,用于確定加載應(yīng)變率;利用樣品的全回收測量及微觀表征,確定了無氧銅環(huán)、柱殼的局域斷裂應(yīng)變及斷裂模式。實(shí)驗(yàn)表明,隨著應(yīng)變率的增加,TU1材料的平均斷裂應(yīng)變增加,斷面的收縮程度加劇,并在應(yīng)變率約為1.0×104s-1附近會出現(xiàn)明顯的斷裂模式轉(zhuǎn)變,斷面收縮率出現(xiàn)量級上的跳躍,從100變化至約103,局域斷裂應(yīng)變呈現(xiàn)明顯的分區(qū)現(xiàn)象。

    頸縮;斷裂;膨脹環(huán);高應(yīng)變率;斷裂應(yīng)變

    材料及結(jié)構(gòu)在高應(yīng)變率載荷加載下的動力學(xué)響應(yīng),涉及材料的快速變形及應(yīng)變局域化、動態(tài)斷裂等相關(guān)過程。材料在高應(yīng)變率下的力學(xué)行為與準(zhǔn)靜態(tài)加載下不同,需要考慮慣性效應(yīng)與應(yīng)變率效應(yīng)帶來的影響[1],特別是應(yīng)變率敏感材料,材料的屈服應(yīng)力,斷裂應(yīng)變等都會隨著應(yīng)變率的變化而發(fā)生變化[2]。此外,斷裂現(xiàn)象的發(fā)生是一個(gè)動態(tài)過程,涉及應(yīng)力波的傳播與能量的快速釋放,研究表明,能量總是會選擇最快的釋放路徑[3]。因此,不同載荷作用下的動態(tài)斷裂將展現(xiàn)出不同的特征形貌。為了研究各種類型的斷裂現(xiàn)象,人們采用了不同的比較準(zhǔn)則,其中斷裂應(yīng)變作為可直接測量的物理量被廣泛使用。然而,針對不同的斷裂模式,與之相關(guān)的斷裂應(yīng)變的選擇標(biāo)準(zhǔn)則會出現(xiàn)較大的差別。在厚壁圓筒膨脹斷裂研究中,湯鐵鋼等[4]采用裂紋張開后爆轟產(chǎn)物泄漏時(shí)刻的工程應(yīng)變作為衡量柱殼承受動態(tài)載荷程度的標(biāo)準(zhǔn)。隨著診斷技術(shù)的進(jìn)步,任國武等[5-6]在實(shí)驗(yàn)中利用高速攝影、多點(diǎn)密排激光干涉測速等技術(shù)可精密診斷外壁裂紋萌生時(shí)刻的應(yīng)變。然而需要注意的是,上述兩類應(yīng)變都是采用殼體環(huán)向的平均變形來計(jì)算,雖然可以作為衡量結(jié)構(gòu)承載程度的標(biāo)準(zhǔn),但這種基于平均假定的計(jì)算方法,會抹平局域頸縮帶來的應(yīng)變局域化特征,它只是一種近似,并不能準(zhǔn)確描述局域變形特征。另一方面,從定性描述的角度出發(fā),人們將材料的斷裂分為脆性斷裂與韌性斷裂,但這種定性描述并不能準(zhǔn)確描述脆性斷裂與韌性斷裂的程度。比如對于金屬這類韌性較好的材料,不同的應(yīng)變率加載下,宏觀上可能都展現(xiàn)出韌性斷裂的特征,然而其變形局域化及斷口形貌可能會出現(xiàn)較大的差別,甚至?xí)霈F(xiàn)斷裂模式的轉(zhuǎn)變。H.Zhang等[7]在研究中發(fā)現(xiàn),金屬環(huán)在膨脹斷裂過程中存在3種斷裂模式:(1)diffuse necking;(2)sheet necking;(3)shear banding。雖然斷面收縮率較小的diffuse necking與斷面收縮率較大的sheet necking都屬于韌性斷裂,然而顯著的塑性局域化差異已經(jīng)不能單純的使用韌脆性這一偏宏觀唯象的描述來表征。

    為了準(zhǔn)確描述金屬環(huán)、柱殼膨脹斷裂過程中結(jié)構(gòu)的局域化特征,本文中嘗試引入描述局域化變形與斷裂程度的斷面收縮率及局域斷裂應(yīng)變;利用電磁膨脹環(huán)實(shí)驗(yàn)平臺開展一維環(huán)的拉伸斷裂研究,并利用新型的線起爆加載平臺,開展金屬環(huán)、柱殼膨脹斷裂實(shí)驗(yàn)研究;通過各類實(shí)驗(yàn)研究,重點(diǎn)關(guān)注平均斷裂應(yīng)變、斷面收縮率及局域斷裂應(yīng)變隨應(yīng)變率的變化關(guān)系,用以闡述平均斷裂應(yīng)變、局域斷裂應(yīng)變在描述結(jié)構(gòu)綜合抗變形斷裂能力以及材料韌性方面的應(yīng)用。

    1 斷面收縮率與局域斷裂應(yīng)變

    H.Zhang等[8]發(fā)現(xiàn),回收破片中的應(yīng)變并不均勻,因此認(rèn)為平均斷裂應(yīng)變并不是一個(gè)很好的描述試樣環(huán)變形的物理量。事實(shí)上,試樣斷裂之后的環(huán)向平均應(yīng)變是結(jié)構(gòu)宏觀抗變形能力的綜合表征,無法準(zhǔn)確展現(xiàn)局域的斷裂特征,這由其定義所決定。為了能夠比較不同載荷強(qiáng)度下的材料與結(jié)構(gòu)的斷裂程度,本文中引入斷面的收縮率來定量給出塑性局域化程度。針對本文實(shí)驗(yàn)中發(fā)現(xiàn)的頸縮導(dǎo)致的兩種斷裂形式diffuse necking、sheet necking,定義斷面收縮率為初始尺寸與斷面截面之比,對于環(huán)結(jié)構(gòu):

    式中:S0為初始截面面積,Sf為斷面面積。對于柱殼,由于軸向變化相對較小,故斷面收縮率近似為:

    式中:h0為柱殼的初始厚度,hf為斷面厚度。由上述定義可以看出,斷面收縮率越大,頸縮程度越嚴(yán)重。

    若材料在對應(yīng)的載荷作用下體積變化不明顯,則可視為不可壓材料,由dV=d(Sl)=Sdl+ldS=0可得

    式中:l為與斷裂面相垂直方向的尺度,在本文討論的兩類情況下,表示環(huán)向尺度。于是可以通過測量斷面的面積變化得到樣品局域的斷裂應(yīng)變

    因此,為了準(zhǔn)確描述環(huán)、柱殼膨脹斷裂特征,本文中選取的斷裂應(yīng)變分為兩類:第一,平均斷裂應(yīng)變,為斷面/斷口完全張開時(shí)刻的環(huán)向應(yīng)變;第二,局域斷裂應(yīng)變,由斷面收縮率決定,如公式(4)所示。

    2 實(shí)驗(yàn)結(jié)果及其分析

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

    采用爆炸膨脹環(huán)實(shí)驗(yàn)平臺(圖1(a))以及電磁膨脹環(huán)實(shí)驗(yàn)平臺(圖1(b))開展實(shí)驗(yàn)研究。樣品材料均為高純無氧銅TU1,樣品從同一根棒材截取加工而成。爆炸膨脹環(huán)(外徑42mm、內(nèi)徑40mm、高2mm)、柱殼(外徑42mm、內(nèi)徑40mm、高80mm)實(shí)驗(yàn)中,驅(qū)動器為20鋼,炸藥為黑索金(RDX)粉末或泰安(PETN)粉末填充而成,通過改變裝藥的直徑獲得不同應(yīng)變率的加載狀態(tài);電磁膨脹環(huán)(外徑42mm、內(nèi)徑40mm、高1mm)實(shí)驗(yàn)中,通過改變加載電壓獲得不同應(yīng)變率的加載狀態(tài)。TU1試樣環(huán)、柱殼實(shí)物圖如圖2所示。

    圖1 爆炸膨脹環(huán)實(shí)驗(yàn)平臺[9]和電磁膨脹環(huán)實(shí)驗(yàn)平臺示意圖[10]Fig.1Sketch of experimental setup of explosive expanding ring and electromagnetic expanding ring

    實(shí)驗(yàn)采用中國工程物理研究院流體物理研究所自主研制的激光干涉測速系統(tǒng)(Doppler pins system,DPS),獲得TU1樣品膨脹斷裂過程的徑向速度歷史,以此計(jì)算加載的應(yīng)變率。樣品的斷裂模式通過斷面的形貌進(jìn)行判斷,平均斷裂應(yīng)變采用樣品完全拉斷時(shí)刻的環(huán)向膨脹量進(jìn)行估計(jì)。膨脹環(huán)實(shí)驗(yàn)中,采用全回收裝置,可以利用回收樣品幾何測量獲得最終的平均斷裂應(yīng)變;而在柱殼的膨脹斷裂中,很難做到全回收,于是我們同時(shí)使用了高速攝影技術(shù),獲得柱殼膨脹斷裂的演化圖像,利用圖像判讀柱殼在不同應(yīng)變率加載下的平均斷裂應(yīng)變;斷面收縮率及局域斷裂應(yīng)變利用SEM進(jìn)行斷口的微觀表征幾何測量獲得。

    圖2 高純無氧銅TU1環(huán)和柱殼樣品實(shí)物圖Fig.2Specimen of oxygen-free high-conductivity copper(TU1)

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

    實(shí)驗(yàn)采用DPS測速技術(shù)、高速攝影及破片回收分析,獲得電磁膨脹環(huán)、爆炸膨脹環(huán)(柱殼)實(shí)驗(yàn)數(shù)據(jù)如表1所示,表中ε·為加載應(yīng)變率,εf-average為平均斷裂應(yīng)變,ψ為斷面收縮率,εf-local為局域斷裂應(yīng)變其中應(yīng)變率為加載峰值應(yīng)變率。

    表1 電磁膨脹環(huán)、爆炸膨脹環(huán)(柱殼)實(shí)驗(yàn)數(shù)據(jù)Table 1Strain rate and fracture strain of experiments

    從表1可以看出,隨著應(yīng)變率的增加,電磁膨脹環(huán)的平均斷裂應(yīng)變也隨之而增加,然而,僅僅通過試樣完全拉斷時(shí)刻的環(huán)向平均應(yīng)變,無法描述與斷裂模式相關(guān)的信息。為了給出斷裂模式對應(yīng)變率的依賴關(guān)系,我們在SEM上對膨脹環(huán)樣品的斷面進(jìn)行了微觀表征。圖3中分別給出了電磁膨脹環(huán)與爆炸膨脹環(huán)的斷面形貌,從圖中可以看出,在爆炸驅(qū)動加載與電磁驅(qū)動加載下,樣品的斷裂模式均為局域的頸縮拉伸導(dǎo)致的韌性斷裂。不同之處在于,爆炸膨脹環(huán)斷面收縮率較小,在本文的實(shí)驗(yàn)狀態(tài)下約為4,且斷面存在明顯的韌窩(圖3(d));而電磁膨脹環(huán)由于試樣環(huán)中存在大電流,溫升較高,在環(huán)向失穩(wěn)之后,由于局域的熱軟化存在,頸縮程度相比而言更為嚴(yán)重。加載電壓為3.0kV時(shí)(圖3(a)),樣品局域的表觀斷面收縮率為6,然而我們注意到樣品在斷口局域出現(xiàn)了較為明顯的燒蝕熔化現(xiàn)象,表觀的斷面收縮率并非純粹由熱軟化導(dǎo)致,而是耦合了其它物理(局部放電)效應(yīng),不適合用來判斷局域斷裂應(yīng)變。而加載電壓為3.5kV(圖3(b))、4.0kV(圖3(c))時(shí),斷口局域并未出現(xiàn)明顯燒蝕,且斷面存在一定程度的韌窩,所以可以較為準(zhǔn)確反映材料自身的特征,其斷面收縮率分別為約為78、13。

    實(shí)驗(yàn)中,利用高速分幅相機(jī)拍攝了TU1柱殼膨脹斷裂過程,典型的破壞圖像如圖4所示,其中圖4(a)為實(shí)驗(yàn)5結(jié)果,其應(yīng)變率約為8.6×103s-1,圖4(b)為實(shí)驗(yàn)6結(jié)果,應(yīng)變率約為1.2×104s-1。通過高速攝影圖像可以看出,應(yīng)變率較低時(shí),斷裂模式為diffuse necking,斷面收縮率較??;而當(dāng)應(yīng)變率較大時(shí),斷裂模式為sheet necking,斷面收縮率較大。

    為了更清晰地展現(xiàn)不同斷裂模式下的斷口特征,并對斷面收縮率進(jìn)行定量測量,利用SEM對回收樣品斷口形貌進(jìn)行了表征,圖5中分別給出了 TU1柱殼在應(yīng)變率為8.6×103s-1(圖5(a))與12.0×103s-1(圖5(b))加載下的斷面形貌。其中圖5(a))的斷面出現(xiàn)大小不同的韌窩,是典型的韌性斷裂;但圖5(b))卻呈現(xiàn)出類似刀刃狀的斷口。通過斷口幾何測量,發(fā)現(xiàn)應(yīng)變率較低時(shí),柱殼徑向的收縮率約為2.5;而當(dāng)應(yīng)變率較高時(shí),由于斷面過于鋒利,斷面收縮率很難精確測量,約為500~1 000的量級。雖然加載應(yīng)變率變化不大,卻出現(xiàn)了完全不同的斷裂模式,斷面收縮率出現(xiàn)了顯著的量級增加。

    圖3 電磁膨脹環(huán)、爆炸膨脹環(huán)典型斷面形貌Fig.3Fracture feature of expanding rings under electromagnetic and explosive loading

    圖4 爆炸膨脹柱殼高速攝影圖像Fig.4High speed camera images of expanding cylindrical shell

    圖5 膨脹柱殼拉伸破片及斷面微觀表征圖片F(xiàn)ig.5Fragments and microscopic fracture surface of expanding cylindrical shell

    2.3 平均斷裂應(yīng)變、局域斷裂應(yīng)變對應(yīng)變率的依賴關(guān)系

    通過SEM的分析已經(jīng)發(fā)現(xiàn),隨著應(yīng)變率的增加,斷裂模式發(fā)生了轉(zhuǎn)變,為了對這一過程進(jìn)行量化,本節(jié)將嘗試尋找斷裂模式轉(zhuǎn)變與斷裂應(yīng)變之間的半定量關(guān)系。研究發(fā)現(xiàn)應(yīng)變率對材料與結(jié)構(gòu)的變形與斷裂起著非常重要的影響,D.Grady[2]指出,斷裂應(yīng)變εf與應(yīng)變率ε·之間滿足關(guān)系εf∝ε·2/3。為了更好地理解TU1斷裂應(yīng)變與應(yīng)變率之間的關(guān)系,我們利用文獻(xiàn)[11]中M態(tài)TU1膨脹環(huán)碎裂數(shù)據(jù)計(jì)算其平均斷裂應(yīng)變,如表2所示。在膨脹環(huán)實(shí)驗(yàn)研究中,加載應(yīng)變率尚未達(dá)到104s-1,為了與膨脹柱殼的應(yīng)變率進(jìn)行對比,選擇文獻(xiàn)[12]中關(guān)于TU1環(huán)數(shù)值模擬的數(shù)據(jù),并將這些數(shù)據(jù)與本實(shí)驗(yàn)結(jié)果繪制于同一張圖中,如圖6所示。

    表2 不同應(yīng)變率下M態(tài)TU1環(huán)碎裂數(shù)據(jù)[11]Table 2Fragmentations data of M state TU1rings at varying strain rates

    圖6 電磁膨脹環(huán)、爆炸膨脹環(huán)/柱殼中平均斷裂應(yīng)變與應(yīng)變率的關(guān)系Fig.6Relationship between strain rate and average fracture strain in electromagnetic ring and explosive expanding ring/cylinder

    結(jié)果表明,文獻(xiàn)[11]中的TU1膨脹環(huán)應(yīng)變率與斷裂應(yīng)變之間的關(guān)系,與本文電磁膨脹環(huán)的數(shù)據(jù)落在同一個(gè)區(qū)域,電磁膨脹環(huán)的斷裂主要以sheet necking模式為主,斷面的收縮率約在101~102的量級。不同載荷形式下的TU1環(huán)、柱殼中平均斷裂應(yīng)變隨著應(yīng)變率的增加,呈現(xiàn)出先增加,然后在達(dá)到一定值之后趨于穩(wěn)定,這一趨勢的過渡恰好位于應(yīng)變率104s-1附近,斷裂模式發(fā)生轉(zhuǎn)變之處。這一改變對于爆炸加載下的TU1樣品而言,意味著應(yīng)變率誘導(dǎo)的斷裂模式轉(zhuǎn)變,與之對應(yīng)的斷面收縮率則發(fā)生了量級的跳躍。爆炸膨脹環(huán)、柱殼結(jié)構(gòu),在低于此應(yīng)變率時(shí),斷面的收縮率(~100)較小,以diffuse necking為主,而高于此應(yīng)變率時(shí),斷面的收縮率(102~103)急劇變大,以sheet necking為主。需要特別指出的是,在應(yīng)變率-平均斷裂應(yīng)變圖上,需要特殊的標(biāo)記才能給出不同斷裂模式之間的分區(qū),在圖6中利用虛線進(jìn)行區(qū)分。除非事先知道實(shí)驗(yàn)結(jié)果,否則單從應(yīng)變率-平均斷裂應(yīng)變圖上,很難進(jìn)行區(qū)分。這一現(xiàn)象從另個(gè)一側(cè)面說明,采用平均思想獲得的斷裂應(yīng)變,本質(zhì)上而言,很難準(zhǔn)確抓住不同斷裂模式的特征。

    圖7 斷面收縮率隨應(yīng)變率的變化Fig.7Variation of area reduction with strain rate

    為了更清晰地定量描述斷裂模式的轉(zhuǎn)變,圖7中給出了電磁膨脹環(huán)、爆炸膨脹環(huán)/柱殼斷面收縮率隨著應(yīng)變率的變化關(guān)系。

    電磁膨脹環(huán)在3.0kV加載下,由于無氧銅韌性較好,應(yīng)變率較低時(shí)局域斷裂較慢,于是頸縮區(qū)域會有充裕的時(shí)間進(jìn)行局部放電,由此導(dǎo)致的燒蝕熔化影響了對斷面收縮率的解讀。為此,在圖7中只給出了加載3.5、4.0kV的數(shù)據(jù)。在材料體積變化不顯著的情況下,利用公式(4)可以在斷裂收縮率與局域斷裂應(yīng)變之間建立起定量關(guān)系,圖8中給出了電磁與爆炸加載下平均斷裂應(yīng)變與局域斷裂應(yīng)變隨應(yīng)變率的變化。

    圖8 爆炸與電磁加載下膨脹環(huán)、柱殼局域斷裂應(yīng)變與平均斷裂應(yīng)變Fig.8Local fracture strain and average fracture strain of explosive expanding ring/cylinder(EE-ring/EE-Cylinder)and electromagnetic ring(EM-ring)

    結(jié)果發(fā)現(xiàn),高應(yīng)變率加載下無氧銅環(huán)的平均斷裂應(yīng)變(0.333)高于低應(yīng)變率加載下的平均斷裂應(yīng)變(0.224);但高應(yīng)變率加載下的斷面收縮率(13)反而要低于低應(yīng)變率加載下的斷面收縮率(78)。這種反轉(zhuǎn)現(xiàn)象產(chǎn)生的原因在于,平均斷裂應(yīng)變是通過回收樣品的環(huán)向伸長量去計(jì)算的,既包含了均勻變形區(qū)域在環(huán)向的變形貢獻(xiàn),又包含了頸縮點(diǎn)、斷裂點(diǎn)對變形區(qū)域的貢獻(xiàn),高應(yīng)變率加載下頸縮點(diǎn)、斷裂點(diǎn)的數(shù)目多于低應(yīng)變加載,因此高應(yīng)變率加載下的平均斷裂應(yīng)變要高于低應(yīng)變率加載下的平均斷裂應(yīng)變,因此平均斷裂應(yīng)變是樣品環(huán)向整體的變形表征,材料與結(jié)構(gòu)抗變形能力的體現(xiàn);而局域斷裂應(yīng)變是通過斷口的收縮率確定的,是斷裂點(diǎn)處的變形表征,是材料的抗變形能力的體現(xiàn),隨著應(yīng)變率的增加,局域斷裂應(yīng)變的變小可能與高應(yīng)變率脆斷特征相關(guān)。這一發(fā)現(xiàn)告訴我們,利用平均斷裂應(yīng)變評估材料的韌性,可能會極大地低估材料的抗變形能力,而斷裂點(diǎn)的局域斷裂應(yīng)變則能更好地評估材料的抗變形能力。此外,在所有的加載狀態(tài)下,平均斷裂應(yīng)變均顯著低于局域斷裂應(yīng)變。這一特征與平均斷裂應(yīng)變、局域斷裂應(yīng)變的計(jì)算相關(guān)。

    同時(shí),我們通過塑性功的計(jì)算發(fā)現(xiàn),爆炸加載下膨脹環(huán)在初期受到?jīng)_擊載荷之后自由膨脹,在結(jié)構(gòu)失穩(wěn)之前溫升較小,約為20~30℃,應(yīng)變率大約在104s-1以下時(shí),斷面的收縮率在100的量級,這一區(qū)域?qū)?yīng)較低應(yīng)變率加載下的爆炸膨脹環(huán)/柱殼;而當(dāng)應(yīng)變率在104s-1之上時(shí),斷裂模式轉(zhuǎn)變導(dǎo)致斷面收縮率急劇增加,因此在斷面收縮率圖(見圖7)上展現(xiàn)出較為明顯的分區(qū)現(xiàn)象。另一方面,由于電磁膨脹環(huán)在加載階段環(huán)內(nèi)存在感應(yīng)電流,將導(dǎo)致約200~300℃的溫升,溫度的升高相當(dāng)于附加給試樣環(huán)較高的環(huán)境溫度,帶來比爆炸加載更大的軟化效應(yīng)。溫度的軟化效應(yīng)使得金屬環(huán)更容易發(fā)生塑性流動,從而導(dǎo)致斷裂應(yīng)變的增加,斷面收縮率約101~102。所以與爆炸膨脹環(huán)相比,相同應(yīng)變率加載下的電磁膨脹環(huán)斷裂應(yīng)變更大。因此在圖7中,電磁加載下的斷面收縮率要明顯高于爆炸加載下的斷面收縮率。

    通過上述的分析我們發(fā)現(xiàn),爆炸膨脹環(huán)、柱殼斷裂模式的轉(zhuǎn)變,主要以應(yīng)變率效應(yīng)為主導(dǎo),隨著應(yīng)變率的增加,超過一定的臨界值(TU1材料約為104s-1)之后,斷面的收縮率會出現(xiàn)量級的跳躍,從100變化至約103;而爆炸膨脹環(huán)與電磁膨脹環(huán)之間的斷裂模式轉(zhuǎn)變則是溫度效應(yīng)與應(yīng)變率效應(yīng)共同主導(dǎo),相同應(yīng)變率加載下,電磁膨脹環(huán)的局域變形程度要明顯高于爆炸膨脹環(huán)。同時(shí),我們通過圖6~8之間的對比可以發(fā)現(xiàn),在較低應(yīng)變率的狀態(tài)時(shí),爆炸膨脹環(huán)與電磁膨脹環(huán)的環(huán)向平均斷裂應(yīng)變相對比較接近,但斷面的收縮率卻落在完全不同的區(qū)域,這直接導(dǎo)致了平均斷裂應(yīng)變與局域斷裂應(yīng)變的顯著差異。這一特征告訴我們,描述斷裂現(xiàn)象尤其是涉及斷裂模式轉(zhuǎn)變時(shí),除了平均斷裂應(yīng)變之外,還需要斷面收縮率/局域斷裂應(yīng)變來確定不同應(yīng)力狀態(tài)下的斷裂特征差異。

    3 結(jié) 論

    本研究利用電磁膨脹環(huán),爆炸膨脹環(huán)、柱殼實(shí)驗(yàn),研究了高純無氧銅TU1的動態(tài)拉伸斷裂行為,為了更好地刻畫斷裂模式的轉(zhuǎn)變,在平均斷裂應(yīng)變的基礎(chǔ)上引入了斷面收縮率以及局域斷裂應(yīng)變,初步得到如下結(jié)論:

    (1)高純無氧銅TU1環(huán)、柱殼的平均斷裂應(yīng)變隨著應(yīng)變率的增加而增加,且在應(yīng)變率約為104s-1附近會出現(xiàn)斷裂模式轉(zhuǎn)變,平均斷裂應(yīng)變增長率驟然變緩,但斷面收縮率從約100變化至102~103,發(fā)生數(shù)量級的跳躍;低于此應(yīng)變率,高純無氧銅TU1環(huán)、柱殼的斷裂模式以diffuse necking為主,高于此應(yīng)變率,斷裂模式以sheet necking為主。

    (2)在較低應(yīng)變率加載下,雖然爆炸膨脹環(huán)與電磁膨脹環(huán)的平均斷裂應(yīng)變相對接近,但其斷面收縮率出現(xiàn)約1~2個(gè)量級的差異,斷面收縮率-應(yīng)變率圖上,呈現(xiàn)出相對孤立的區(qū)域,與之相應(yīng)的局域斷裂應(yīng)變同樣差異顯著;其原因在于電磁膨脹環(huán)實(shí)驗(yàn)中,樣品環(huán)在加載初期即存在大電流,導(dǎo)致了環(huán)的整體軟化,使之更易于進(jìn)行塑性流動,隨后在頸縮局域產(chǎn)生更大的塑性變形,增加了斷面收縮率。

    (3)所有的加載狀態(tài)下,平均斷裂應(yīng)變均顯著低于局域斷裂應(yīng)變,原因在于平均斷裂應(yīng)變的計(jì)算既包含了均勻變形區(qū)域在環(huán)向的變形貢獻(xiàn),又包含了頸縮點(diǎn)、斷裂點(diǎn)對變形區(qū)域的貢獻(xiàn),而局域斷裂應(yīng)變只體現(xiàn)斷裂點(diǎn)處的變形特征;所以可以認(rèn)為平均斷裂應(yīng)變是材料與結(jié)構(gòu)抗變形能力的綜合體現(xiàn),局域斷裂應(yīng)變是材料抗變形能力的體現(xiàn)。

    從更為廣泛的角度而言,在材料與結(jié)構(gòu)的動態(tài)斷裂現(xiàn)象描述中,斷裂應(yīng)變使用時(shí)需要重點(diǎn)區(qū)分平均斷裂應(yīng)變與局域斷裂應(yīng)變,平均斷裂應(yīng)變是材料與結(jié)構(gòu)抗變形能力的表征,與結(jié)構(gòu)的失穩(wěn)相關(guān),是材料動態(tài)性能、樣品的宏觀幾何結(jié)構(gòu)以及加載歷史耦合所致;而局域斷裂應(yīng)變則主要是材料動態(tài)性能與加載歷史相互作用所致。由于結(jié)構(gòu)中尚存在較大比例的均勻變形區(qū)域,利用平均斷裂應(yīng)變評估材料的韌性,會極大低估材料的抗變形能力。因此,除了環(huán)向平均斷裂應(yīng)變之外,還需要引入描述頸縮程度的斷面收縮率、局域斷裂應(yīng)變及頸縮點(diǎn)斷裂點(diǎn)的數(shù)目,這樣才能更完整地展現(xiàn)材料與結(jié)構(gòu)的斷裂特征。

    感謝張世文在文章撰寫過程中對本文的指導(dǎo)幫助。

    [1] 王禮立.應(yīng)力波基礎(chǔ)[M].2版.北京:國防工業(yè)出版社,2005.

    [2] Grady D.Fragmentation of rings and shells-The legacy of N.F.Mott[M].Berlin:Springer,2006.

    [3] 周風(fēng)華,郭麗娜,王禮立.脆性固體碎裂過程中的最快卸載特性[J].固體力學(xué)學(xué)報(bào),2010,31(3):286-295.Zhou Fenghua,Guo Lina,Wang Lili.The rapidest unloading characteristics in the fragmentation process of brittle solids[J].Chinese Journal of Solid Mechanics,2010,31(3):286-295.

    [4] 湯鐵鋼,李慶忠,孫學(xué)林,等.45鋼柱殼膨脹斷裂的應(yīng)變率效應(yīng)[J].爆炸與沖擊,2006,26(2):129-133.Tang Tiegang,Li Qingzhong,Sun Xuelin,et al.Strain-rate effects of expanding fracture of 45steel cylinder shells driven by detonation[J].Explosion and Shock Waves,2006,26(2):129-133.

    [5] 任國武,郭昭亮,湯鐵鋼,等.高應(yīng)變率加載下金屬柱殼斷裂的實(shí)驗(yàn)研究[J].兵工學(xué)報(bào),2016,37(1):77-82.Ren Guowu,Guo Zhaoliang,Tang Tiegang,et al.Experimental research on fracture of metal case under loading at high strain rate[J].Acta Armamentarii,2016,37(1):77-82.

    [6] Ren Guowu,Guo Zhaoliang,F(xiàn)an Cheng,et al.Dynamic shear fracture of an explosively-driven metal cylindrical shell[J].International Journal of Impact Engineering,2016,95:35-39.

    [7] Zhang H,Ravi-Chandar K.Dynamic fragmentation of ductile materials[J].Journal of Physics D:Applied Physics,2009,42(21):214010.

    [8] Zhang H,Ravi-Chandar K.On the dynamics of necking and fragmentation-I.Real-time and post-mortem observations in Al 6061-O[J].International Journal of Fracture,2006,142:183-217.

    [9] 湯鐵鋼,李慶忠,陳永濤,等.實(shí)現(xiàn)材料高應(yīng)變率拉伸加載的爆炸膨脹環(huán)技術(shù)[J].爆炸與沖擊,2009,29(5):546-549.Tang Tiegang,Li Qingzhong,Chen Yongtao,et al.An improved technique for dynamic tension of metal ring by explosive loading[J].Explosion and Shock Waves,2009,29(5):546-549.

    [10] 桂毓林,孫承緯,李強(qiáng),等.實(shí)現(xiàn)金屬環(huán)動態(tài)拉伸的電磁加載技術(shù)研究[J].爆炸與沖擊,2006,26(6):481-485.Gui Yulin,Sun Chengwei,Li Qiang,et al.Experimental studies on dynamic tension of metal ring by electromagnetic loading[J].Explosion and Shock Waves,2006,26(6):481-485.

    [11] 桂毓林.電磁加載下金屬膨脹環(huán)的動態(tài)斷裂與碎裂研究[D].四川綿陽:中國工程物理研究院,2007.

    [12] 陳磊,周風(fēng)華,湯鐵鋼.韌性金屬環(huán)高速膨脹碎裂過程的有限元模擬[J].力學(xué)學(xué)報(bào),2011,43(5):861-870.Chen Lei,Zhou Fenghua,Tang Tiegang.Finite element simulation of the high velocity expansion and fragmentation of ductile metallic rings[J].Chinese Journal of Theoretical and Applied Mechanics,2011,43(5):861-870.

    Fracture mode transition in expanding ring and cylindrical shell under electromagnetic and explosive loadings

    Guo Zhaoliang,F(xiàn)an Cheng,Liu Mingtao,Ren Guowu,Tang Tiegang,Liu Cangli
    (Institute of Fluid Physics,China Academy of Engineering Physics,Mianyang621999,Sichuan,China)

    In the present study,we designed the electromagnetic and explosive driving expanding ring/cylinder experiments and investigated the expanding fracture characteristics of oxygen-free high-conductivity copper(OFHC)in consideration of the conception of the reduction of area,the local fracture strain and the average fracture strain.We used a high speed camera to record the fracture process and obtain the fracture strain of the copper cylinder and the Doppler pins system (DPS)to obtain the radial velocity of the specimen in order to achieve the strain rate of the loading.We verified the local fracture strain and the fracture mode by analyzing the soft-recovered fragments of the expanding ring and the cylinder.Based on the experimental results,we found that the average fracture strain and the reduction of the area increases as does the strain rate.Moreover,the fracture mode transition may occur at the strain rate of about 1.0×104s-1,and the reduction of the area may increase by an order of magnitude,i.e.from the order of 100to that of 103,and the local fracture strain exhibits an obvious subarea.

    necking;fracture;expanding ring;high strain rate;fracture strain

    O347.3 國標(biāo)學(xué)科代碼:13015

    A

    10.11883/1001-1455(2017)06-1072-08

    2016-04-21;

    2016-09-21

    國家自然科學(xué)基金項(xiàng)目(11172279)

    郭昭亮(1984— ),男,博士研究生,助理研究員;通信作者:范 誠,fancheng@caep.cn。

    (責(zé)任編輯 曾月蓉)

    猜你喜歡
    柱殼局域收縮率
    滌綸短纖維干熱收縮率測量不確定度分析
    影響滌綸短纖維干熱收縮率因素探討
    聚丙烯材料收縮率影響因素研究
    上海塑料(2021年3期)2022-01-06 14:05:02
    電磁驅(qū)動薄壁鋁合金柱殼的動態(tài)膨脹變形特征
    母線生產(chǎn)中端面線收縮率與寬厚比間關(guān)系的統(tǒng)計(jì)規(guī)律
    山東冶金(2018年5期)2018-11-22 05:12:20
    局域積分散列最近鄰查找算法
    電子測試(2018年18期)2018-11-14 02:30:34
    偏心柱殼自由振動的級數(shù)變換求解方法
    1/3含口蓋復(fù)合材料柱殼后屈曲性能
    PET成像的高分辨率快速局域重建算法的建立
    基于局域波法和LSSVM的短期負(fù)荷預(yù)測
    電測與儀表(2015年7期)2015-04-09 11:39:50
    亚洲成a人片在线一区二区| 日本撒尿小便嘘嘘汇集6| 最近视频中文字幕2019在线8| 久久久久国产精品人妻aⅴ院| 一卡2卡三卡四卡精品乱码亚洲| 国内毛片毛片毛片毛片毛片| 老司机在亚洲福利影院| 精品99又大又爽又粗少妇毛片 | 国产精品亚洲av一区麻豆| 国产精品女同一区二区软件 | 两人在一起打扑克的视频| 观看美女的网站| 99热只有精品国产| 好男人电影高清在线观看| www日本在线高清视频| 欧美色欧美亚洲另类二区| 久久中文字幕人妻熟女| 免费高清视频大片| 久久人人精品亚洲av| 18禁黄网站禁片免费观看直播| 少妇丰满av| 国产一区在线观看成人免费| 精品国内亚洲2022精品成人| 国产成人欧美在线观看| 女警被强在线播放| 偷拍熟女少妇极品色| 亚洲国产欧洲综合997久久,| 日韩精品青青久久久久久| 熟女电影av网| 小蜜桃在线观看免费完整版高清| 在线观看免费午夜福利视频| 亚洲国产欧美人成| 岛国视频午夜一区免费看| 国产亚洲av高清不卡| 91老司机精品| 给我免费播放毛片高清在线观看| 哪里可以看免费的av片| 亚洲欧美日韩卡通动漫| 亚洲av日韩精品久久久久久密| 天堂√8在线中文| 黄片小视频在线播放| avwww免费| 无遮挡黄片免费观看| 美女午夜性视频免费| 亚洲美女视频黄频| 久久中文看片网| 99在线人妻在线中文字幕| 麻豆国产97在线/欧美| 欧美性猛交黑人性爽| 淫秽高清视频在线观看| 亚洲九九香蕉| 成人无遮挡网站| 欧美极品一区二区三区四区| 婷婷亚洲欧美| 1024香蕉在线观看| 最近在线观看免费完整版| 真实男女啪啪啪动态图| 小说图片视频综合网站| 欧美黄色片欧美黄色片| 久99久视频精品免费| 国产成人av激情在线播放| 久久香蕉精品热| 欧美午夜高清在线| 亚洲国产色片| 亚洲无线观看免费| 在线观看一区二区三区| 精品久久久久久久末码| 日本一本二区三区精品| 国产高潮美女av| 一二三四社区在线视频社区8| 青草久久国产| 国产伦精品一区二区三区四那| 一区二区三区高清视频在线| 中文字幕熟女人妻在线| 精品久久蜜臀av无| 一a级毛片在线观看| 日韩欧美免费精品| 色av中文字幕| 国内精品一区二区在线观看| 黄色女人牲交| 精品欧美国产一区二区三| 精品久久久久久久久久免费视频| 成人av一区二区三区在线看| 嫩草影院精品99| 久久久久性生活片| 在线观看舔阴道视频| 18禁观看日本| 我要搜黄色片| 欧美国产日韩亚洲一区| 成人特级黄色片久久久久久久| 一进一出好大好爽视频| 后天国语完整版免费观看| 精品国内亚洲2022精品成人| 男女午夜视频在线观看| 国产精品久久久久久人妻精品电影| 成人三级做爰电影| 久久精品91蜜桃| 国产精品一区二区三区四区久久| 热99re8久久精品国产| 午夜日韩欧美国产| av国产免费在线观看| 成人午夜高清在线视频| 国产 一区 欧美 日韩| 久久国产精品影院| 床上黄色一级片| 欧美日韩一级在线毛片| 国产精品99久久久久久久久| 九九久久精品国产亚洲av麻豆 | 91在线观看av| 男女下面进入的视频免费午夜| 免费人成视频x8x8入口观看| 男女那种视频在线观看| av天堂在线播放| 亚洲九九香蕉| 久久久久精品国产欧美久久久| 精品不卡国产一区二区三区| 亚洲精品国产精品久久久不卡| 欧美日韩福利视频一区二区| 搞女人的毛片| 久久中文看片网| 最近在线观看免费完整版| 国产1区2区3区精品| 97碰自拍视频| cao死你这个sao货| 免费观看精品视频网站| 精品久久久久久,| 亚洲欧美日韩卡通动漫| 亚洲欧美日韩高清在线视频| 久久人人精品亚洲av| 国产高清有码在线观看视频| 久久人妻av系列| 一区二区三区激情视频| 国产精品爽爽va在线观看网站| 国产视频内射| 小说图片视频综合网站| 露出奶头的视频| 又爽又黄无遮挡网站| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久久电影 | 亚洲国产色片| 亚洲欧美日韩东京热| www日本黄色视频网| 最近最新免费中文字幕在线| 成人精品一区二区免费| 国内揄拍国产精品人妻在线| 亚洲国产色片| 色老头精品视频在线观看| 国产精品一及| 免费在线观看成人毛片| 亚洲精品一卡2卡三卡4卡5卡| 岛国视频午夜一区免费看| 可以在线观看的亚洲视频| 久久精品国产综合久久久| 国产真人三级小视频在线观看| 嫩草影院入口| 99久久精品一区二区三区| 男女午夜视频在线观看| 亚洲五月婷婷丁香| 国产伦在线观看视频一区| 一a级毛片在线观看| 男插女下体视频免费在线播放| 脱女人内裤的视频| 国产一区二区在线观看日韩 | 18禁黄网站禁片午夜丰满| 成人永久免费在线观看视频| 久久久久久久久久黄片| 国产亚洲欧美98| 亚洲精品一卡2卡三卡4卡5卡| 一二三四社区在线视频社区8| 欧美色欧美亚洲另类二区| 手机成人av网站| 精品乱码久久久久久99久播| 欧美黄色淫秽网站| 亚洲一区二区三区色噜噜| or卡值多少钱| 国产综合懂色| 国产精品久久久av美女十八| 久久午夜亚洲精品久久| 免费在线观看视频国产中文字幕亚洲| 国产成年人精品一区二区| 99热6这里只有精品| 欧美精品啪啪一区二区三区| 久久久水蜜桃国产精品网| av片东京热男人的天堂| 国产精品久久久久久久电影 | 99久久精品国产亚洲精品| 亚洲色图av天堂| 五月玫瑰六月丁香| 成熟少妇高潮喷水视频| 好男人电影高清在线观看| 最近最新中文字幕大全免费视频| 亚洲在线观看片| 人人妻,人人澡人人爽秒播| 亚洲国产精品成人综合色| 国产高清三级在线| 一级a爱片免费观看的视频| 日韩高清综合在线| 老司机深夜福利视频在线观看| 国产免费av片在线观看野外av| 欧美日韩福利视频一区二区| 在线a可以看的网站| 视频区欧美日本亚洲| 成熟少妇高潮喷水视频| 精品人妻1区二区| 日本 欧美在线| 美女被艹到高潮喷水动态| 精华霜和精华液先用哪个| 一本一本综合久久| 亚洲va日本ⅴa欧美va伊人久久| 亚洲,欧美精品.| 国产97色在线日韩免费| 午夜视频精品福利| 一边摸一边抽搐一进一小说| 成人高潮视频无遮挡免费网站| 亚洲av免费在线观看| 国产亚洲欧美98| 日本黄色视频三级网站网址| 久久久精品欧美日韩精品| 亚洲天堂国产精品一区在线| 亚洲国产精品sss在线观看| 久久午夜亚洲精品久久| 一二三四在线观看免费中文在| 欧美日韩中文字幕国产精品一区二区三区| 午夜久久久久精精品| 人人妻,人人澡人人爽秒播| 又大又爽又粗| 天天一区二区日本电影三级| 国产免费av片在线观看野外av| 久久这里只有精品19| 亚洲专区国产一区二区| 男女之事视频高清在线观看| 久久久久久久久免费视频了| 欧美av亚洲av综合av国产av| 可以在线观看的亚洲视频| 国产精品98久久久久久宅男小说| 国产黄a三级三级三级人| 18禁黄网站禁片午夜丰满| 一级毛片精品| 国产爱豆传媒在线观看| 欧美一区二区国产精品久久精品| 中文字幕最新亚洲高清| 国产一区二区在线观看日韩 | 成人亚洲精品av一区二区| 国产精品综合久久久久久久免费| 夜夜爽天天搞| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇熟女久久| 熟女少妇亚洲综合色aaa.| 国产爱豆传媒在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美大码av| 又黄又爽又免费观看的视频| 此物有八面人人有两片| 露出奶头的视频| 99久久综合精品五月天人人| 国产毛片a区久久久久| 丰满人妻一区二区三区视频av | 香蕉av资源在线| 好男人在线观看高清免费视频| 国产激情欧美一区二区| 日本五十路高清| 免费看美女性在线毛片视频| 精品一区二区三区视频在线 | 国产精品亚洲美女久久久| 国产精品香港三级国产av潘金莲| 久久精品综合一区二区三区| 麻豆国产av国片精品| 日本一本二区三区精品| 午夜福利在线观看免费完整高清在 | 欧美一区二区国产精品久久精品| aaaaa片日本免费| 国产免费男女视频| 国产真人三级小视频在线观看| 成人无遮挡网站| 精品国产乱子伦一区二区三区| 精品电影一区二区在线| 久久中文字幕人妻熟女| 国产精品野战在线观看| АⅤ资源中文在线天堂| 国产在线精品亚洲第一网站| 国产精品久久久久久精品电影| 国产97色在线日韩免费| 日本五十路高清| 久久久国产成人免费| 变态另类成人亚洲欧美熟女| 99国产精品一区二区蜜桃av| 超碰成人久久| 日本免费a在线| 美女被艹到高潮喷水动态| 宅男免费午夜| 国产精品 欧美亚洲| 亚洲成人精品中文字幕电影| 日本撒尿小便嘘嘘汇集6| 免费在线观看影片大全网站| 给我免费播放毛片高清在线观看| 国产精品久久久久久久电影 | 法律面前人人平等表现在哪些方面| 久久久久亚洲av毛片大全| 午夜精品在线福利| 一区二区三区激情视频| 亚洲五月婷婷丁香| 黄色成人免费大全| 嫩草影院入口| 免费av不卡在线播放| 亚洲欧美一区二区三区黑人| 久久精品亚洲精品国产色婷小说| 露出奶头的视频| 国产激情欧美一区二区| 国产午夜精品久久久久久| 久久国产精品人妻蜜桃| 精品乱码久久久久久99久播| 久久人妻av系列| 久久婷婷人人爽人人干人人爱| 少妇的丰满在线观看| 成人国产一区最新在线观看| 操出白浆在线播放| 日韩免费av在线播放| 久久性视频一级片| 日韩中文字幕欧美一区二区| www.www免费av| 曰老女人黄片| 日韩av在线大香蕉| 又黄又爽又免费观看的视频| 久久久久久久精品吃奶| 色老头精品视频在线观看| 免费观看人在逋| netflix在线观看网站| 国产午夜精品久久久久久| 国产精品一区二区免费欧美| 88av欧美| 国产aⅴ精品一区二区三区波| 午夜a级毛片| 欧美中文日本在线观看视频| 97人妻精品一区二区三区麻豆| 非洲黑人性xxxx精品又粗又长| 久久久久久人人人人人| 国产成人aa在线观看| 999久久久精品免费观看国产| 亚洲精品久久国产高清桃花| 国产精品1区2区在线观看.| 国产成人精品久久二区二区91| 国产成+人综合+亚洲专区| 久久久久久大精品| 色综合婷婷激情| 19禁男女啪啪无遮挡网站| 看免费av毛片| 国产一区二区在线av高清观看| 国产在线精品亚洲第一网站| 一二三四社区在线视频社区8| 久久久久久久久免费视频了| 十八禁网站免费在线| 男女那种视频在线观看| 十八禁网站免费在线| 黑人欧美特级aaaaaa片| 十八禁网站免费在线| 精品一区二区三区av网在线观看| 哪里可以看免费的av片| 亚洲av中文字字幕乱码综合| 成人特级av手机在线观看| 国产伦精品一区二区三区四那| 日韩三级视频一区二区三区| 成人欧美大片| 成人特级av手机在线观看| 真人做人爱边吃奶动态| 天天添夜夜摸| 国产伦在线观看视频一区| 俄罗斯特黄特色一大片| 久久精品国产综合久久久| 久久久精品欧美日韩精品| 亚洲av片天天在线观看| 在线观看美女被高潮喷水网站 | 成年女人毛片免费观看观看9| 免费看a级黄色片| 在线看三级毛片| 亚洲精品一区av在线观看| 在线看三级毛片| av中文乱码字幕在线| 午夜激情福利司机影院| 亚洲真实伦在线观看| 精品一区二区三区四区五区乱码| 午夜福利在线在线| 少妇的丰满在线观看| 桃红色精品国产亚洲av| 两个人的视频大全免费| 操出白浆在线播放| 国产男靠女视频免费网站| 国产视频内射| 巨乳人妻的诱惑在线观看| 99久久综合精品五月天人人| 久久精品国产亚洲av香蕉五月| 成人18禁在线播放| 宅男免费午夜| 脱女人内裤的视频| 国产一区二区三区在线臀色熟女| 欧美乱色亚洲激情| 99国产精品一区二区蜜桃av| 麻豆成人av在线观看| 欧美成人一区二区免费高清观看 | 男人和女人高潮做爰伦理| 九九在线视频观看精品| 欧美乱妇无乱码| 黄色 视频免费看| 久99久视频精品免费| 日韩精品青青久久久久久| 国产日本99.免费观看| 村上凉子中文字幕在线| 91麻豆av在线| 久久久成人免费电影| 制服人妻中文乱码| 波多野结衣高清作品| 成人欧美大片| 国产精品女同一区二区软件 | 美女cb高潮喷水在线观看 | 日本黄色片子视频| 午夜福利欧美成人| 级片在线观看| 国产黄a三级三级三级人| 少妇的逼水好多| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| 亚洲七黄色美女视频| h日本视频在线播放| 成人无遮挡网站| 亚洲精品在线美女| 好男人在线观看高清免费视频| 国产爱豆传媒在线观看| 成人鲁丝片一二三区免费| 日韩 欧美 亚洲 中文字幕| 亚洲男人的天堂狠狠| 国产伦精品一区二区三区视频9 | 两人在一起打扑克的视频| 亚洲美女黄片视频| 人妻久久中文字幕网| 日韩欧美国产一区二区入口| 午夜激情欧美在线| 久久欧美精品欧美久久欧美| 99久久综合精品五月天人人| 久久精品国产清高在天天线| 免费在线观看日本一区| 国产精品自产拍在线观看55亚洲| 又黄又爽又免费观看的视频| 久久久久九九精品影院| 好男人电影高清在线观看| 最新在线观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 最新美女视频免费是黄的| 变态另类丝袜制服| 国产伦一二天堂av在线观看| 免费搜索国产男女视频| 国产一级毛片七仙女欲春2| 久99久视频精品免费| 午夜影院日韩av| 国产v大片淫在线免费观看| 黄色女人牲交| 窝窝影院91人妻| 国产精品av久久久久免费| 少妇熟女aⅴ在线视频| 欧美成人一区二区免费高清观看 | e午夜精品久久久久久久| 日韩 欧美 亚洲 中文字幕| 亚洲天堂国产精品一区在线| 国产精品九九99| 无限看片的www在线观看| 嫩草影院入口| 国产精品综合久久久久久久免费| 女人被狂操c到高潮| 国产精品香港三级国产av潘金莲| 亚洲天堂国产精品一区在线| 久久久水蜜桃国产精品网| 久久久久免费精品人妻一区二区| 男女视频在线观看网站免费| 不卡一级毛片| 桃色一区二区三区在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品一区二区免费欧美| 国产一区二区三区在线臀色熟女| 亚洲国产精品合色在线| 日本一二三区视频观看| 午夜福利在线在线| 国产久久久一区二区三区| 亚洲国产欧洲综合997久久,| 一本综合久久免费| 午夜福利在线观看吧| 身体一侧抽搐| 久久精品国产99精品国产亚洲性色| 日本 av在线| 天堂网av新在线| 亚洲性夜色夜夜综合| 老司机在亚洲福利影院| 国产真实乱freesex| 色在线成人网| 最新美女视频免费是黄的| 欧美性猛交黑人性爽| 欧美激情久久久久久爽电影| 后天国语完整版免费观看| 2021天堂中文幕一二区在线观| 欧美国产日韩亚洲一区| 国产亚洲欧美在线一区二区| 激情在线观看视频在线高清| 麻豆国产av国片精品| 日韩成人在线观看一区二区三区| 国产综合懂色| av在线天堂中文字幕| 网址你懂的国产日韩在线| 亚洲欧美日韩高清专用| 男女午夜视频在线观看| 中文在线观看免费www的网站| 亚洲av日韩精品久久久久久密| 岛国在线观看网站| 久9热在线精品视频| 日本与韩国留学比较| 天堂影院成人在线观看| 亚洲熟妇熟女久久| 亚洲欧美日韩东京热| 给我免费播放毛片高清在线观看| 国产野战对白在线观看| 最近最新中文字幕大全电影3| 噜噜噜噜噜久久久久久91| 白带黄色成豆腐渣| 日韩有码中文字幕| 日本 av在线| 美女被艹到高潮喷水动态| 最近最新中文字幕大全免费视频| 老汉色∧v一级毛片| 欧美黄色淫秽网站| 亚洲av免费在线观看| 国产伦一二天堂av在线观看| 男人的好看免费观看在线视频| 久久欧美精品欧美久久欧美| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| 黑人巨大精品欧美一区二区mp4| 非洲黑人性xxxx精品又粗又长| 久久这里只有精品中国| 国产精品一区二区免费欧美| 极品教师在线免费播放| 日韩中文字幕欧美一区二区| 日本黄色片子视频| 亚洲精品色激情综合| 国产免费男女视频| 给我免费播放毛片高清在线观看| 少妇人妻一区二区三区视频| 亚洲av成人精品一区久久| 午夜亚洲福利在线播放| 国内精品美女久久久久久| 亚洲第一欧美日韩一区二区三区| 波多野结衣高清作品| www日本在线高清视频| 香蕉国产在线看| 一级a爱片免费观看的视频| 怎么达到女性高潮| 日韩av在线大香蕉| www.自偷自拍.com| 国产av一区在线观看免费| 国产欧美日韩精品亚洲av| 又黄又爽又免费观看的视频| 一进一出抽搐动态| 国产精品av视频在线免费观看| 看黄色毛片网站| 中国美女看黄片| 人人妻,人人澡人人爽秒播| 中文在线观看免费www的网站| 国产精品av久久久久免费| 精品国产超薄肉色丝袜足j| 国产v大片淫在线免费观看| 久久久水蜜桃国产精品网| 成人国产综合亚洲| 美女 人体艺术 gogo| 国产91精品成人一区二区三区| netflix在线观看网站| 国产成人影院久久av| 波多野结衣高清无吗| 欧美大码av| 91av网站免费观看| 999精品在线视频| 亚洲欧美精品综合久久99| 女生性感内裤真人,穿戴方法视频| 婷婷精品国产亚洲av| 中文字幕精品亚洲无线码一区| 91麻豆精品激情在线观看国产| 久久久久久久久中文| 日韩高清综合在线| 97人妻精品一区二区三区麻豆| 99热精品在线国产| 久久亚洲真实| 又粗又爽又猛毛片免费看| 在线播放国产精品三级| 亚洲性夜色夜夜综合| 日本黄大片高清| www.精华液| 欧美一区二区国产精品久久精品| 成年版毛片免费区| xxxwww97欧美| 国产成年人精品一区二区| 国产毛片a区久久久久| 99热6这里只有精品| 深夜精品福利| 亚洲色图av天堂| 国产亚洲欧美98| 五月玫瑰六月丁香| 亚洲欧美日韩东京热| 久99久视频精品免费| 午夜免费观看网址| 欧美精品啪啪一区二区三区| 午夜激情欧美在线| 免费看美女性在线毛片视频| 免费在线观看影片大全网站| 色尼玛亚洲综合影院| 欧美日韩乱码在线| 12—13女人毛片做爰片一| 国内久久婷婷六月综合欲色啪| 日韩欧美免费精品| 18禁观看日本| xxx96com| 久久亚洲真实|