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

    細(xì)晶銅材料力學(xué)性能及藥型罩領(lǐng)域應(yīng)用研究

    2017-03-28 01:16:44黃正祥祖旭東肖強(qiáng)強(qiáng)
    彈道學(xué)報(bào) 2017年1期
    關(guān)鍵詞:粗晶細(xì)晶純銅

    閆 超,黃正祥,祖旭東,肖強(qiáng)強(qiáng),賈 鑫

    (南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094)

    細(xì)晶銅材料力學(xué)性能及藥型罩領(lǐng)域應(yīng)用研究

    閆 超,黃正祥,祖旭東,肖強(qiáng)強(qiáng),賈 鑫

    (南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094)

    為研究細(xì)晶純銅材料藥型罩的侵徹性能,對細(xì)晶純銅材料進(jìn)行力學(xué)性能研究。利用二級輕氣炮加速飛片的方法對細(xì)晶純銅材料進(jìn)行了沖擊加載實(shí)驗(yàn),給出了細(xì)晶純銅材料在70~320 GPa范圍內(nèi)的沖擊波速度和質(zhì)點(diǎn)速度的關(guān)系。計(jì)算得到了適用于LS-DYNA的Gruneisen狀態(tài)方程參數(shù)。利用擬合得到的參數(shù)分析了粗晶材料及細(xì)晶材料所形成的射流在斷裂時(shí)間、侵徹威力方面的差異。根據(jù)仿真,在160 mm最佳炸高進(jìn)行靜破甲實(shí)驗(yàn)。結(jié)果表明,細(xì)晶材料藥型罩的侵徹威力較粗晶材料明顯提高,說明將藥型罩材料的晶粒組織細(xì)化是提高破甲威力的重要手段。

    細(xì)晶純銅;動(dòng)態(tài)力學(xué)性能;聚能射流;侵徹

    在現(xiàn)代戰(zhàn)場上,間隔裝甲、復(fù)合裝甲、爆炸反應(yīng)裝甲在坦克裝甲車輛上的成熟運(yùn)用,使得現(xiàn)階段各國主戰(zhàn)坦克的防御能力大幅提高,這就使破甲戰(zhàn)斗部在與裝甲防護(hù)技術(shù)和防護(hù)性能的博弈中處于明顯的劣勢。因此,提升現(xiàn)有破甲戰(zhàn)斗部的毀傷威力,已經(jīng)成為急需解決的關(guān)鍵性問題。在評價(jià)破甲戰(zhàn)斗部威力時(shí),射流的有效長度和斷裂時(shí)間是衡量破甲能力的重要指標(biāo)。根據(jù)國內(nèi)外關(guān)于射流的理論研究[1],破甲金屬射流是一種超塑性金屬固體流,而金屬的超塑性與其晶粒尺寸的平方成反比。另外,大量的動(dòng)靜態(tài)射流測量和破甲實(shí)驗(yàn)也表明,金屬藥型罩材料的破甲威力與其初始材料晶粒尺寸有著密切關(guān)系,隨著藥型罩材料初始晶粒尺寸的降低,射流的侵徹威力顯著增加。但是傳統(tǒng)加工手段的藥型罩晶粒尺寸一般在幾十微米以上,至多細(xì)化到10 μm。若能將藥型罩材料的晶粒尺寸進(jìn)一步細(xì)化到亞微米甚至納米級別,則有望大幅度提升傳統(tǒng)材料破甲彈的威力。因此,將超細(xì)晶及納米晶材料運(yùn)用到藥型罩的加工生產(chǎn)當(dāng)中,則可為藥型罩侵徹能力的提升開拓新途徑、新方法。

    目前,國內(nèi)外科研單位及學(xué)者針對藥型罩材料晶粒尺寸及織構(gòu)對射流長度和斷裂時(shí)間的影響,開展了大量的研究工作。Bourne等[2]研究了初始藥型罩材料晶粒尺寸和結(jié)晶織構(gòu)對聚能裝藥藥型罩侵徹威力的共同影響,結(jié)果表明藥型罩材料晶粒組織的細(xì)化以及組織均勻性的改變,對聚能裝藥侵徹威力的提升意義重大。Lichtenberger[3]針對純銅材料藥型罩射流的斷裂時(shí)間和晶粒尺寸之間的關(guān)系進(jìn)行了研究,指出當(dāng)藥型罩材料晶粒尺寸在20~90 μm范圍內(nèi)變化時(shí),藥型罩形成金屬射流的斷裂時(shí)間隨著晶粒的細(xì)化而呈現(xiàn)增長趨勢。Walters[4]對晶粒尺寸范圍在10~120 μm內(nèi)的無氧純銅藥型罩侵徹結(jié)果進(jìn)行了研究,得到了與Lichtenberger相類似的結(jié)論。Golaski[5]通過實(shí)驗(yàn)指明射流的斷裂時(shí)間直接受藥型罩材料晶粒尺寸的影響。王鐵福等[6]利用靜破甲威力實(shí)驗(yàn)和X光射線照相技術(shù)對25~70 μm范圍內(nèi)變化的4種紫銅藥型罩進(jìn)行了研究,得出藥型罩材料晶粒度對射流動(dòng)態(tài)延展性有重要影響,材料晶粒細(xì)化可增加射流長度、延緩斷裂時(shí)間的結(jié)論。Hirsh[7]提出了一種能夠解釋細(xì)晶材料破甲射流形成與補(bǔ)充的模型,研究指出當(dāng)金屬晶粒小于5 μm時(shí),其可獲得類似液體表面張力的性質(zhì)。

    本文以提升傳統(tǒng)紫銅藥型罩侵徹威力為目標(biāo),通過等位角擠壓(equal channel angular pressing,ECAP)手段將紫銅材料的晶粒組織進(jìn)行細(xì)化,來實(shí)現(xiàn)破甲能力的提升。針對細(xì)晶純銅材料,采用二級輕氣炮作為加載手段,對細(xì)晶純銅的Gruneisen狀態(tài)方程參數(shù)進(jìn)行計(jì)算,得到了適用于LS-DYNA的細(xì)晶純銅材料參數(shù)。計(jì)算求解了2種材料藥型罩在不同炸高條件下的侵徹威力,最終選取最佳炸高條件進(jìn)行了靜破甲實(shí)驗(yàn)研究,驗(yàn)證仿真參數(shù)的正確性。

    1 細(xì)晶純銅材料的沖擊壓縮實(shí)驗(yàn)

    1.1 細(xì)晶純銅材料的制備

    本次試驗(yàn)所選用的純銅原材料為工業(yè)無氧高導(dǎo)電性純銅TU1,初始的外形尺寸為32 mm厚的熱軋態(tài)板料??紤]到純銅板料的制造過程會(huì)對材料顯微結(jié)構(gòu)造成破壞,因此在后續(xù)的晶粒細(xì)化之前,要對切割好的銅塊(32 mm×32 mm×160 mm)進(jìn)行真空退火處理。經(jīng)過熱處理之后,純銅材料的晶粒尺寸約為100 μm。

    在粗晶以及細(xì)晶純銅藥型罩的加工制造過程中,純銅材料主要經(jīng)歷了如圖1所示的過程,最終形成2種晶粒尺寸的純銅藥型罩。由粗晶純銅材料到超細(xì)晶純銅藥型罩主要經(jīng)歷劇烈的塑性變形晶粒細(xì)化、軋制板材、旋壓成形這三大過程。

    圖1 粗晶及細(xì)晶純銅藥型罩加工過程示意圖

    本文根據(jù)細(xì)晶純銅藥型罩實(shí)際加工制造流程,選取其中3種狀態(tài)的超細(xì)晶純銅材料作為研究對象,開展細(xì)晶純銅材料的動(dòng)高壓加載的沖擊壓縮實(shí)驗(yàn)。選取的材料分別為:①經(jīng)過劇烈塑性變形(ECAP8)之后的純銅,②劇烈塑性變形加交叉軋制(ECAP8+CR)之后的純銅,③劇烈塑性變形加軋制和熱處理(ECAP8+CR+Treated)之后的純銅。上述3種細(xì)晶純銅材料和常規(guī)粗晶純銅在準(zhǔn)靜態(tài)拉伸條件下的力學(xué)性能參數(shù),以及通過EBSD觀察的平均晶粒尺寸見表1,表中,s為晶粒尺寸,ρ為材料密度,σs為屈服強(qiáng)度,σb為抗拉強(qiáng)度,δEB為斷裂延伸率,δVE為均勻延伸率。

    表1 細(xì)晶純銅材料力學(xué)性能

    1.2 實(shí)驗(yàn)方法及測量原理

    針對3種狀態(tài)的細(xì)晶純銅材料,采用二級輕氣炮加速發(fā)射飛片的方法,使得飛片與靶板撞擊,對靶板內(nèi)部的材料進(jìn)行沖擊壓縮。實(shí)驗(yàn)整體的方法原理如圖2所示,圖中,vw為飛片速度。

    圖2 沖擊壓縮實(shí)驗(yàn)原理圖

    在碰撞面前后的物理量滿足如下的質(zhì)量、動(dòng)量、能量守恒方程:

    ρ0(vD-vu0)=ρ1(vD-vu1)

    (1)

    p1-p0=ρ0(vD-vu0)(vu1-vu0)

    (2)

    (3)

    式中:p為壓強(qiáng),ρ為密度,vu為質(zhì)點(diǎn)速度,e為比內(nèi)能,vD為樣品中的沖擊波速度,V為比體積,下標(biāo)0表示沖擊壓縮前狀態(tài),1表示沖擊壓縮后狀態(tài)。上述5個(gè)待測物理量,需要測量出其中的任意2個(gè),再根據(jù)式(1)~式(3)就可以計(jì)算出其余的參數(shù)??紤]到現(xiàn)有的技術(shù)手段及測量準(zhǔn)確性等因素,本文選擇測量沖擊波速度vD及粒子速度vu。沖擊波速度vD采用直接測量的方法,粒子速度則是根據(jù)阻抗匹配法進(jìn)行測量??紤]到本文研究的沖擊壓力范圍較高,因此選用2種標(biāo)材:無氧銅和Ta,其中選擇Ta是為了實(shí)現(xiàn)更高的沖擊壓力。當(dāng)標(biāo)材為無氧銅時(shí),屬于對稱碰撞,此時(shí)飛片的速度vw等于樣品中粒子速度的2倍,即vw=2vu。當(dāng)選擇Ta為標(biāo)材時(shí),待測材料中的沖擊波后粒子速度可以表示為

    (4)

    式中:

    F=ρ0tλt

    G=-ρ0tC0f-2λtvwρ0t-ρ0svD

    (5)

    H=ρ0tvw(λtvw+C0f)

    式中:ρ0f為標(biāo)材的密度;C0f,λf分別為標(biāo)材的Hugoniot參數(shù);ρ0s為待測樣品的密度。vw,vD通過實(shí)驗(yàn)測量得到。因此就可以通過式(4)、式(5)求解出待測樣品中的沖擊波后粒子速度vu。通過上述手段就可以將vu的測量轉(zhuǎn)變?yōu)閺椝賤w的測量。

    本次實(shí)驗(yàn)是在中國工程物理研究院φ28 mm口徑二級輕氣炮試驗(yàn)系統(tǒng)上進(jìn)行的。彈丸的擊靶速度采用激光遮斷式測試技術(shù)(OBB)測量,通過安裝在發(fā)射管末端的2組光源,當(dāng)彈丸通過時(shí)會(huì)依次遮擋光束,此時(shí)示波器會(huì)顯示2個(gè)脈沖的響應(yīng)時(shí)間差,再結(jié)合兩光源的距離,就可以得到彈丸的速度vw。樣品中的沖擊波速度采用成熟的電探針測試技術(shù)測量,電探針的布局如圖3所示。

    圖3 電探針布局圖

    樣品前表面16根探針,后表面7根探針,另外還有4根觸發(fā)探針(T1,T2,T3,T4),4根接地探針(G1,G2,G3,G4)。利用波陣面還原法計(jì)算樣品中的沖擊波傳播時(shí)間,將樣品前表面的探針響應(yīng)時(shí)間及空間坐標(biāo)帶入下式:

    t=A1x+A2y+A3x2+A4y2+A5xy+A6

    (6)

    式中:t為探針的響應(yīng)時(shí)間;x,y為探針的平面坐標(biāo);A1,A2,…,A6為方程系數(shù)。

    對式(6)進(jìn)行擬合,計(jì)算曲面方程系數(shù),再將樣品后界面探針的空間位置坐標(biāo)帶入到曲面方程中,可以得到樣品前界面沖擊波所對應(yīng)的到達(dá)時(shí)間tb,i,樣品前界面的沖擊波響應(yīng)時(shí)間為tf,i,則樣品中沖擊波經(jīng)過的時(shí)間ta計(jì)算如下:

    (7)

    式中:k表示樣品后界面探針的個(gè)數(shù)。利用測量得到的樣品厚度就可以計(jì)算出樣品中的沖擊波速vD。

    1.3 實(shí)驗(yàn)結(jié)果與分析

    根據(jù)實(shí)驗(yàn)具體測量的彈丸速度vw、樣品中的沖擊波速度vD、通過阻抗匹配法計(jì)算得到的粒子速度vu及其他相關(guān)參數(shù)的計(jì)算結(jié)果見表2。表中,φ為沖擊波速度計(jì)算得出的相對測量不確定度,pH為沖擊壓力計(jì)算值,V為沖擊壓縮后的比體積,V0為沖擊壓縮前的比體積。

    表2 3種細(xì)晶純銅材料的沖擊壓縮實(shí)驗(yàn)數(shù)據(jù)

    表2中的1-1發(fā)實(shí)驗(yàn)數(shù)據(jù),由于飛片未精加工,實(shí)驗(yàn)并不滿足精密物理實(shí)驗(yàn)的要求,因此可信度差,實(shí)驗(yàn)數(shù)據(jù)無效。其余3種材料的vD-vu關(guān)系如下:

    vD1=4.058(±0.029)+1.476(±0.049)vu1

    (8)

    vD2=3.934(±0.054)+1.506(±0.058)vu2

    (9)

    vD3=4.000(±0.075)+1.487(±0.026)vu3

    (10)

    式中:各量下標(biāo)的序號和表2中的材料號一一對應(yīng)。3種材料的vD-vu關(guān)系顯示測量結(jié)果較為接近,說明晶粒細(xì)化并不能夠影響材料在超高壓條件下的沖擊壓縮性能。前處理手段的不同,改變的僅僅是材料顯微結(jié)構(gòu)組織均勻性等因素,并未徹底改變材料的組成成分。而對某種材料的沖擊壓縮性能而言,影響的主要因素卻是初始密度及組成成分。針對本次實(shí)驗(yàn)所進(jìn)行的10發(fā)有效數(shù)據(jù),由于3種材料的測量結(jié)果都比較相近,因此可以近似地看成是針對破甲彈藥型罩所用細(xì)晶純銅材料在70~320 GPa壓力范圍內(nèi)所測量得到的沖擊壓縮數(shù)據(jù)。將上述所有數(shù)據(jù)進(jìn)行線性擬合,圖4為所有數(shù)據(jù)擬合得到的vD-vu關(guān)系,可以得到vD-vu關(guān)系式:

    vD=4.007(±0.06)+1.485(±0.02)vu

    (11)

    圖4 數(shù)據(jù)擬合得到的細(xì)晶銅vD-vu關(guān)系

    若要擬合細(xì)晶純銅材料的Gruneisen狀態(tài)方程系數(shù),就需要對其格林愛森系數(shù)γ隨比體積變化的關(guān)系式進(jìn)行求解,目前主要采用經(jīng)驗(yàn)公式的方法進(jìn)行求解[8]:

    (12)

    式中:常態(tài)格林愛森系數(shù)γ0采用:

    (13)

    式中:c為流體聲速,可以近似地等于零壓下的材料聲速c0=4.007km/s;αv為體積膨脹系數(shù),cV為定容比熱,KT為等溫體積模量。由于細(xì)晶純銅材料是一種新型材料,許多材料參數(shù)仍然是未知的,因此在計(jì)算求解過程中,對于一些材料的物性常量,利用純銅材料的數(shù)值進(jìn)行代替。αv,cV等材料參數(shù)的取值參考文獻(xiàn)[8]。將上述數(shù)據(jù)帶入式(12),計(jì)算求得常態(tài)下的格林愛森系數(shù)γ0=2.07。因此格林愛森系數(shù)隨著比體積的變化關(guān)系可以表示為

    (14)

    細(xì)晶純銅材料的Gruneisen狀態(tài)方程參數(shù)見表3,表中參數(shù)C,S1,S2,S3是利用飛片沖擊壓縮得到的vD-vu關(guān)系獲得,γ0為Gruneisen系數(shù);a為γ0的一階體積修正。

    表3 細(xì)晶純銅材料的Gruneisen狀態(tài)方程參數(shù)

    2 數(shù)值仿真計(jì)算及實(shí)驗(yàn)研究

    2.1 數(shù)值計(jì)算模型及材料參數(shù)

    采用有限元軟件LS-DYNA建立全尺寸靜破甲實(shí)驗(yàn)的數(shù)值計(jì)算模型。聚能裝藥采用φ56 mm無殼體裝藥,大量的彈靶實(shí)驗(yàn)證明該種成型裝藥結(jié)構(gòu)合理,破甲性能穩(wěn)定。基準(zhǔn)裝藥的基本尺寸見圖5。藥型罩材料都是純銅,兩者僅在微觀晶粒尺寸上存在差異。后效靶的材料為45#鋼。

    圖5 基準(zhǔn)裝藥的尺寸

    炸藥爆炸壓垮藥型罩形成射流及射流與靶板的相互作用,屬于典型的多物質(zhì)相互作用的大變形運(yùn)動(dòng)。采用Van Leer ALE算法進(jìn)行模擬。該算法為二階精度算法,可較精確地模擬炸藥爆炸過程中動(dòng)量與能量之間的轉(zhuǎn)化。

    針對上述的實(shí)際物理模型,采用有限元軟件進(jìn)行建模。建立的模型如圖6所示。模型整體由四部分區(qū)域組成,其中,炸藥、藥型罩、空氣域采用三維Euler網(wǎng)格建模,后效靶采用三維Lagrange建模,這2種網(wǎng)格采用ALE算法耦合在一起。考慮到模型屬于典型的回轉(zhuǎn)體結(jié)構(gòu),因此選擇1/4模型進(jìn)行建模來減少計(jì)算量,提高計(jì)算效率。在模型的對稱面上施加對稱滑移約束,將空氣域的邊界設(shè)置為非反射邊界。

    圖6 有限元模型

    炸藥裝藥為JH-2炸藥,材料模型采用高能炸藥爆轟模型和JWL狀態(tài)方程。空氣域的材料模型采用Null模型和LINEAR_POLYNOMIAL狀態(tài)方程[9]。藥型罩材料為粗晶純銅和超細(xì)晶純銅,后效靶材料為45#鋼,材料模型都采用Johnson-Cook模型和Gruneisen狀態(tài)方程。其中粗晶純銅和45#鋼的材料參數(shù)參考文獻(xiàn)[10]。細(xì)晶純銅材料的材料參數(shù)根據(jù)準(zhǔn)靜態(tài)及動(dòng)態(tài)拉伸實(shí)驗(yàn)得到的結(jié)果進(jìn)行擬合,見表4,表中,A為考慮應(yīng)變率及溫度下的屈服強(qiáng)度;B,n為應(yīng)變率強(qiáng)化系數(shù);C為應(yīng)變率敏感系數(shù);m為溫度軟化系數(shù)。

    表4 細(xì)晶純銅材料的J-C模型相關(guān)參數(shù)

    2.2 仿真計(jì)算結(jié)果與分析

    2.2.1 晶粒尺寸對射流斷裂時(shí)間的影響

    根據(jù)3.1節(jié)建立的有限元模型,在不添加后效靶的情況下研究材料晶粒尺寸對射流斷裂時(shí)間的影響。分別截取粗晶純銅射流出現(xiàn)斷裂的時(shí)刻56 μs,細(xì)晶純銅射流出現(xiàn)斷裂的時(shí)刻74 μs,以及2種射流都出現(xiàn)斷裂后80 μs的射流形態(tài),計(jì)算結(jié)果如圖7所示。從圖中不同時(shí)刻的射流形態(tài)對比可直觀地發(fā)現(xiàn),粗晶純銅藥型罩形成的射流在56 μs時(shí)射流頭部就已經(jīng)出現(xiàn)頸縮斷裂,而同時(shí)刻細(xì)晶純銅藥型罩形成的射流形態(tài)完好,仍然處于穩(wěn)定的長度增加階段。直至74 μs,細(xì)晶純銅藥型罩形成的射流才出現(xiàn)頭部斷裂,相比于粗晶藥型罩,射流開始出現(xiàn)斷裂的時(shí)間延長18 μs。更晚的射流斷裂時(shí)間則意味著形成射流的有效長度更長。對比同時(shí)刻的粗晶純銅藥型罩的射流,其射流的頭部已經(jīng)可以觀測到較為明顯的連續(xù)斷裂。隨著時(shí)間的繼續(xù)推移,細(xì)晶純銅藥型罩的射流頭部才出現(xiàn)斷斷續(xù)續(xù)的射流斷裂,但斷裂的程度并不是很明顯,只有頭部一小部分發(fā)生與射流主體的分離,總的來說新增加的斷裂數(shù)量有限,射流形態(tài)仍然具有很好的一致性。粗晶純銅在80 μs時(shí)形成射流,射流斷裂的范圍繼續(xù)增加,整體顯現(xiàn)出較長范圍的顆粒狀斷裂射流。

    圖7 不同時(shí)刻超細(xì)晶純銅與粗晶純銅射流形態(tài)對比

    從仿真結(jié)果的對比情況來看,細(xì)晶藥型罩形成的射流具有較長的累積射流和較遲的射流斷裂時(shí)間,且明顯優(yōu)于粗晶藥型罩形成的射流,這說明藥型罩材料初始晶粒細(xì)化有助于射流性能的增強(qiáng)。同時(shí),由于聚能射流的超塑性是在高應(yīng)變率下由動(dòng)態(tài)再結(jié)晶引起的,而動(dòng)態(tài)再結(jié)晶的難易、快慢,又與材料的晶粒度密切相關(guān),因此材料的初始晶粒越細(xì)其動(dòng)態(tài)再結(jié)晶也就越容易,再結(jié)晶的速度也就越快。最終導(dǎo)致的結(jié)果就是連續(xù)射流的長度增加。

    2.2.2 炸高對藥型罩侵徹威力的影響

    相同裝藥結(jié)構(gòu)的藥型罩形成的射流在飛行過程中受到慣性力的作用而不斷拉長,頭部速度也會(huì)降低,故藥型罩具有最佳炸高。細(xì)晶材料由于具有更好的動(dòng)態(tài)延展性,因此在不同炸高條件下,其表現(xiàn)出的侵徹威力可能與粗晶材料不一致。因此針對相同裝藥結(jié)構(gòu)不同晶粒尺寸的2種材料開展炸高對比研究,具有一定的科學(xué)意義。

    利用3.1節(jié)中數(shù)值計(jì)算模型,計(jì)算2種材料藥型罩在不同炸高時(shí)侵徹性能的差異。炸高選取80 mm,120 mm,160 mm,200 mm,240 mm,280 mm,320 mm,360 mm,炸高的變化范圍為1.5~6.5倍裝藥口徑。計(jì)算后效靶的侵徹深度對比見圖8。圖中,h為炸高,P為侵徹深度。

    從上述的計(jì)算結(jié)果可以看出,細(xì)晶藥型罩和粗晶藥型罩的侵徹深度都隨著炸高的增加呈現(xiàn)先增加后減小的趨勢,且都在3倍裝藥口徑附近達(dá)到最佳侵徹威力,但侵徹威力始終是細(xì)晶材料優(yōu)于粗晶材料。隨著炸高的繼續(xù)增加,在超過4倍裝藥口徑時(shí),2種材料藥型罩的侵徹威力開始出現(xiàn)降低,但降幅不盡相同。粗晶純銅藥型罩的侵徹深度隨著炸高的增加降幅明顯,而細(xì)晶藥型罩則降幅平緩。2種材料在相同炸高條件下的侵徹深度增幅δ如表5所示。

    圖8 不同炸高條件下2種藥型罩侵徹深度對比

    表5 2種材料藥型罩不同炸高條件下侵徹深度增幅對比

    伴隨著炸高的增加,在5倍裝藥口徑炸高范圍內(nèi),侵徹深度增幅主要呈現(xiàn)出先增加后減小的趨勢,且在3倍裝藥口徑附近取得最大值。在高于5倍裝藥口徑的情況下,粗晶藥型罩侵徹深度明顯減小,而細(xì)晶藥型罩的侵徹深度減小不明顯,使得在大炸高的情況下,2種材料的侵徹深度增幅又逐漸增大。造成這種情況的原因:在80 mm,120 mm的小炸高情況下,由于射流并未得到完全伸長,因此超細(xì)晶材料較粗晶材料侵徹深度的增幅有限,不能體現(xiàn)出細(xì)晶材料在超塑性方面的優(yōu)勢。隨著炸高的增加,粗晶純銅藥型罩在炸高160 mm時(shí)侵徹深度最大,達(dá)到236 mm,約為4倍裝藥口徑。細(xì)晶純銅藥型罩也在炸高160 mm時(shí)侵徹深度最大,達(dá)到273 mm,約為5倍裝藥口徑。細(xì)晶純銅藥型罩的侵徹深度隨著炸高的增加,上升速度快且幅度大。隨著炸高的進(jìn)一步增加,在200~360 mm范圍時(shí),細(xì)晶純銅藥型罩侵徹深度隨著炸高的增加呈現(xiàn)緩慢減小的趨勢,而粗晶純銅藥型罩的侵徹深度隨著炸高的增加而降幅明顯。特別是在320 mm,360 mm炸高情況下,侵深增幅增加,達(dá)到23.9%和23.4%,說明在大炸高情況下,雖然2種材料形成射流的整體侵徹威力較低炸高而言減小,但是細(xì)晶材料射流由于具有更好的延展性,在大炸高的情況下仍然可以保持一定的侵徹威力。這主要是因?yàn)樵诖笳ǜ邨l件下,射流頭部在未接觸到靶板之前,由于速度梯度的存在,頭部產(chǎn)生拉伸斷裂,形成顆粒狀的不連續(xù)射流,最終形成斷裂侵徹,使得侵徹性能下降。

    2.3 實(shí)驗(yàn)驗(yàn)證

    仿真結(jié)果顯示,粗晶及細(xì)晶材料藥型罩的最佳炸高都在160 mm,因此對2種材料在此炸高條件下進(jìn)行靜破甲實(shí)驗(yàn)研究。試驗(yàn)現(xiàn)場的布置如圖9所示,主要包括雷管、成型裝藥、炸高筒、靶板四部分,其中靶板由不同高度的圓柱形45#鋼疊加而成。

    圖9 靜破甲實(shí)驗(yàn)現(xiàn)場布置

    表6 實(shí)驗(yàn)結(jié)果對比

    采用細(xì)晶純銅藥型罩的裝藥,在160 mm炸高條件下侵徹深度分別達(dá)到了282 mm,280 mm,278 mm,與相同裝藥結(jié)構(gòu)的粗晶藥型罩在最佳炸高條件下的平均侵徹深度244 mm相比,侵徹深度增加均在14%以上。同在160 mm炸高下,實(shí)驗(yàn)與仿真誤差為5.6%,這說明根據(jù)力學(xué)實(shí)驗(yàn)擬合的J-C材料模型參數(shù)及輕氣炮實(shí)驗(yàn)測量的細(xì)晶純銅狀態(tài)方程參數(shù)真實(shí)可靠,可以為細(xì)晶純銅藥型罩的研究提供參考。

    相同炸高條件下,粗晶純銅與細(xì)晶純銅藥型罩的侵徹孔形貌見圖10。粗晶銅藥型罩的平均侵徹深度為244 mm,侵徹孔截面均為規(guī)整圓形,射流在軸線方向上保持很好的一致性。進(jìn)一步從剖面侵徹孔通道形貌分析可知,侵徹孔前半段過渡均勻,未見凹凸不均變化,而侵徹孔底部則見明顯“波浪形”變化。分析可知:射流在侵徹到末端時(shí),射流尾部由于存在速度差,出現(xiàn)明顯的分段斷裂,但是射流的方向并沒有太大的變化,故斷裂的射流不斷形成新的侵徹孔。相比于粗晶純銅藥型罩,細(xì)晶純銅藥型罩平均侵徹深度為280 mm,較粗晶純銅藥型罩平均侵徹深度增加了14.7%,表明細(xì)晶純銅藥型罩的侵徹威力明顯優(yōu)于粗晶態(tài)。從侵徹孔形貌分析可知,細(xì)晶純銅形成的射流侵徹孔不僅侵徹深度增加明顯,侵徹孔內(nèi)壁也更加光滑平直,同時(shí)孔底部出現(xiàn)“波浪形”的范圍也更長。由此可以看出,細(xì)晶純銅藥型罩形成的射流長度要比粗晶純銅藥型罩長得多,但超細(xì)晶純銅藥型罩形成的射流尾部速度可能較低,低于45#鋼的臨界侵徹速度,可以通過選擇高能炸藥來進(jìn)一步提高超細(xì)晶純銅藥型罩的侵徹深度。

    圖10 炸高160 mm侵徹孔形貌對比

    3 結(jié)論

    本文通過力學(xué)性能實(shí)驗(yàn),研究了細(xì)晶純銅材料的動(dòng)態(tài)力學(xué)性能,得到了適用于LS-DYNA的材料模型參數(shù)及狀態(tài)方程參數(shù)。利用仿真模型分析了藥型罩晶粒尺寸對射流斷裂時(shí)間及侵徹威力的影響,最后根據(jù)仿真結(jié)果選擇最佳炸高進(jìn)行了靜破甲實(shí)驗(yàn)驗(yàn)證,得到的主要結(jié)論如下。

    ①準(zhǔn)靜態(tài)拉伸實(shí)驗(yàn)結(jié)果顯示,經(jīng)過ECAP8及ECAP8+CR之后的純銅材料,屈服強(qiáng)度分別達(dá)到360 MPa和400 MPa,較粗晶狀態(tài)時(shí)增加70 MPa。同時(shí)根據(jù)EBSD微觀組織觀察,材料的晶粒尺寸被細(xì)化到500 nm左右,在經(jīng)過軋制處理后,晶粒尺寸進(jìn)一步減小至300 nm。

    ②根據(jù)實(shí)驗(yàn)測量數(shù)據(jù),通過阻抗匹配法計(jì)算得到?jīng)_擊波后粒子速度vu。并擬合得到了3種材料在70~320 GPa范圍內(nèi)的vD-vu關(guān)系。由于測量結(jié)果較為接近,可將3種細(xì)晶純銅材料近似地認(rèn)為是一種應(yīng)用于破甲戰(zhàn)斗部的細(xì)晶純銅材料,得到?jīng)_擊波速度vD和波后粒子速度的關(guān)系為vD=4.007(±0.006)+1.485(±0.02)vu,并計(jì)算得到了細(xì)晶純銅的Gruneisen系數(shù)隨比體積變化的關(guān)系:γ=2.07(V/V0)。

    ③仿真結(jié)果表明,藥型罩材料的初始晶粒尺寸對聚能射流有著重要的影響。在相同裝藥結(jié)構(gòu)下,細(xì)晶材料形成的射流有效長度及斷裂時(shí)間更長。在160 mm最佳炸高情況下,細(xì)晶材料平均侵徹深度 280 mm,粗晶材料平均侵徹深度244 mm,增幅達(dá)14.7%。

    ④由實(shí)驗(yàn)分析可知,細(xì)晶純銅藥型罩侵徹性能較粗晶純銅藥型罩而言提升明顯,因此可以利用細(xì)晶純銅替代現(xiàn)有的粗晶純銅作為藥型罩材料。

    [1] PETIT J,JEANCLAUDE V,FRESSENGEAS C.Breakup of copper shaped-charge jets:experiment,numerical simulations,and analytical modeling[J].Journal of Applied Physics,2005,98(12):123521.

    [2] BOURNE B,JONES P N,WARREN R H.Grain size and crystallographic texture effects on the performance of shaped charge[C]//14th Internatioanl Symposium on Ballistics.Quebec,Canada:IBC,1993:19-124.

    [3] LICHTENBERGER A,SCHARF M,BOHMANN A.Influence of the structural and metallurgical state of a liner on the performance of a shaped charge:CO 218/8[R].Saint-Louis:French-German Research Institute,1981.

    [4] WALTERS W.Introduction to shaped charges:MD21005-5069[R].Aberdeen,MD:Army Research Laboratory Aberdeen Proving Ground,2007.

    [5] GOLASKI S K.Effect of liner grain size on shaped charge jet performance and characteristics:BRL-TR-2800[R].Aberdeen,MD:US Army Ballistic Research Laboratory,1987.

    [6] 王鐵福,王雷.藥型罩材料的晶粒度對射流性能的影響[J].高壓物理學(xué)報(bào),1996,10(4):291-298. WANG Tie-fu,WANG Lei.The effects of the grain size of a liner on the performance of shaped charge jets[J].Chinese Journal of High Pressure Physics,1996,10(4):291-298.(in Chinese)

    [7] 董文健.大口徑純鎢藥型罩CVD制備工藝的優(yōu)化與破甲性能的研究[D].北京:北京理工大學(xué),2015. DONG Wen-jian.Study of CVD process optimization for large diameter pure tungsten liner and penetration performance[D].Beijing:Beijing Institute of Technology,2015.(in Chinese)

    [8] 湯文輝,張若棋.物態(tài)方程理論及計(jì)算概論[M].北京:高等教育出版社,2008. TANG Wen-hui,ZHANG Ruo-qi.Introduction theory and computation of equations of state[M].Beijing:High Education Press,2008.(in Chinese)

    [9] 尚曉江,蘇建宇,王化鋒.ANSYS/LS-DYNA動(dòng)力分析方法與工程實(shí)例[M].北京:中國水利水電出版社,2006. SHANG Xiao-jiang,SU Jian-yu,WANG Hua-feng.ANSYS/LS-DYNA dynamic analysis method and the engineering examples[M].Beijing:China Water and Power Press,2006.(in Chinese)

    [10] JOHNSON G R,COOK W H.A constitutive model and data for metals subjected to large strains,high strain rates and high temperatures[C]//Proceedings of the 7th International Symposium on Ballistics.[S.l.]:IBC,1983:541-547.

    Mechanical Property of Ultra-fine Grained Copper Material and Study on Application of Liner Field

    YAN Chao,HUANG Zheng-xiang,ZU Xu-dong,XIAO Qiang-qiang,JIA Xin

    (School of Mechanical Engineering,Nanjing University of Science and Technology,Nanjing 210094,China)

    In order to study the penetration performance of ultra-fine grained(UFG)copper liner,the mechanical property of UFG copper was studied.The impact experiment of UFG copper were conducted by two-stage light-gas-gun planar impact technique.The relation between shock velocity and particle velocity was described in the pressure range from 70 GPa to 320 GPa.The parameters of Gruneisen EOS applied to LS-DYNA were calculated.The differences of break-up time and penetration power between the coarse grained copper and UFG copper were analyzed by fitting parameters.The penetration experiment was carried out by choosing the best burst-height 160 mm.The results show that the penetration power of UFG copper liner significantly improves compared with the coarse grained copper liner,and it is the effective way to improve the penetration power to reduce the grain size of liner material.

    ultra-fine grained copper;dynamic mechanical property;shaped charge jet;penetration

    2016-04-27

    閆超(1991- ),男,碩士研究生,研究方向?yàn)椴牧系膭?dòng)態(tài)力學(xué)行為。E-mail:18601408248@163.com。

    TJ413.2

    A

    1004-499X(2017)01-0085-08

    猜你喜歡
    粗晶細(xì)晶純銅
    抑制6xxx系制品拉毛缺陷和粗晶暴露的平衡氧化工藝研究
    超細(xì)鉬銅復(fù)合粉體及細(xì)晶鉬銅合金的制備
    淺析陽極氧化工藝對粗晶缺陷的影響
    遇上了行家
    制備參數(shù)對HCCM水平連鑄純銅板坯組織與力學(xué)性能的影響
    多硫化鈉法表征純銅帶表面清潔度
    滑移爆轟條件下高純銅的層裂行為
    高性能細(xì)晶鎢及鎢合金的研究進(jìn)展
    退火時(shí)間對細(xì)晶高強(qiáng)IF 鋼的織構(gòu)和晶界特征分布的影響
    上海金屬(2015年4期)2015-11-29 01:12:26
    超低碳鋼熱軋盤條表面環(huán)形粗晶組織及控制措施
    上海金屬(2015年1期)2015-11-28 06:01:07
    a级毛色黄片| 精品酒店卫生间| 日韩欧美一区视频在线观看 | 国产黄色视频一区二区在线观看| 啦啦啦中文免费视频观看日本| 中文精品一卡2卡3卡4更新| 99久久中文字幕三级久久日本| 亚洲激情五月婷婷啪啪| 国产高清不卡午夜福利| 亚洲经典国产精华液单| 黄色配什么色好看| 色5月婷婷丁香| 91狼人影院| 亚洲丝袜综合中文字幕| 国产综合懂色| 九九爱精品视频在线观看| 中文字幕制服av| 秋霞在线观看毛片| 老司机影院成人| 最近手机中文字幕大全| 久久午夜福利片| 久久精品久久精品一区二区三区| 中国国产av一级| 久久久精品94久久精品| 亚洲av成人av| 国产精品av视频在线免费观看| 少妇丰满av| 亚洲怡红院男人天堂| 亚洲国产精品专区欧美| 亚洲成人精品中文字幕电影| 精品人妻视频免费看| 久久热精品热| 1000部很黄的大片| 美女脱内裤让男人舔精品视频| 久久久久久久久久久免费av| 国产精品一区二区三区四区免费观看| 亚洲欧美日韩卡通动漫| 99热网站在线观看| 国产精品国产三级国产av玫瑰| 高清日韩中文字幕在线| 亚洲av男天堂| av免费在线看不卡| 美女主播在线视频| 91av网一区二区| 欧美bdsm另类| 亚洲伊人久久精品综合| 日韩精品有码人妻一区| 99久久精品热视频| 国产 亚洲一区二区三区 | 97人妻精品一区二区三区麻豆| 免费观看a级毛片全部| 国产不卡一卡二| 美女主播在线视频| 日韩中字成人| 婷婷色综合大香蕉| 精品一区二区免费观看| 免费电影在线观看免费观看| 国产精品人妻久久久影院| 麻豆国产97在线/欧美| 欧美另类一区| 搡老乐熟女国产| 深爱激情五月婷婷| 日日干狠狠操夜夜爽| 久久97久久精品| 大片免费播放器 马上看| 少妇猛男粗大的猛烈进出视频 | 欧美性猛交╳xxx乱大交人| 日韩强制内射视频| 日本av手机在线免费观看| 免费观看性生交大片5| 午夜福利视频精品| 亚洲精品第二区| 高清在线视频一区二区三区| 联通29元200g的流量卡| kizo精华| 黄色配什么色好看| 国产黄片美女视频| av天堂中文字幕网| 寂寞人妻少妇视频99o| 纵有疾风起免费观看全集完整版 | 七月丁香在线播放| 亚洲欧洲日产国产| 日日摸夜夜添夜夜爱| 成年人午夜在线观看视频 | 97精品久久久久久久久久精品| 国精品久久久久久国模美| 国产成人福利小说| 国产一区亚洲一区在线观看| 精品久久国产蜜桃| 天美传媒精品一区二区| 一级毛片我不卡| 能在线免费看毛片的网站| 男插女下体视频免费在线播放| 久久久久网色| 国产中年淑女户外野战色| 亚洲精品456在线播放app| 内地一区二区视频在线| 亚洲精品乱码久久久v下载方式| 国产一区有黄有色的免费视频 | 熟女电影av网| 欧美日韩综合久久久久久| 午夜福利网站1000一区二区三区| 只有这里有精品99| 久久久久精品久久久久真实原创| 国产精品久久久久久精品电影| 又粗又硬又长又爽又黄的视频| 日韩强制内射视频| 69人妻影院| 毛片女人毛片| 亚洲精品成人av观看孕妇| 国产乱来视频区| 国产在线一区二区三区精| av黄色大香蕉| 久久这里有精品视频免费| 日韩国内少妇激情av| 高清午夜精品一区二区三区| 91午夜精品亚洲一区二区三区| 3wmmmm亚洲av在线观看| 亚洲伊人久久精品综合| 日韩av免费高清视频| 欧美日韩视频高清一区二区三区二| 欧美xxxx黑人xx丫x性爽| 男人爽女人下面视频在线观看| 嘟嘟电影网在线观看| 国产一区有黄有色的免费视频 | 国产免费一级a男人的天堂| 亚洲精品视频女| 国产不卡一卡二| 精品久久久久久久久av| 日本av手机在线免费观看| 最后的刺客免费高清国语| 国产人妻一区二区三区在| 国产亚洲一区二区精品| 久久精品国产鲁丝片午夜精品| 少妇的逼水好多| 免费无遮挡裸体视频| 色综合亚洲欧美另类图片| 日日啪夜夜撸| 亚洲国产精品专区欧美| av一本久久久久| 一二三四中文在线观看免费高清| 日韩亚洲欧美综合| 久久99热6这里只有精品| 色尼玛亚洲综合影院| 欧美潮喷喷水| 国产探花极品一区二区| 天堂中文最新版在线下载 | 国产伦精品一区二区三区四那| 日韩在线高清观看一区二区三区| 三级经典国产精品| 男的添女的下面高潮视频| 五月伊人婷婷丁香| 国产av不卡久久| av在线天堂中文字幕| 亚洲综合精品二区| 欧美丝袜亚洲另类| 亚洲精品456在线播放app| av专区在线播放| 国产毛片a区久久久久| 亚洲精品乱码久久久久久按摩| 日本-黄色视频高清免费观看| 2022亚洲国产成人精品| 菩萨蛮人人尽说江南好唐韦庄| 久久久国产一区二区| 69人妻影院| 最近的中文字幕免费完整| 免费黄网站久久成人精品| 秋霞伦理黄片| 亚洲欧美日韩东京热| 夜夜看夜夜爽夜夜摸| 国产精品无大码| 亚洲av福利一区| 校园人妻丝袜中文字幕| 亚洲国产最新在线播放| 简卡轻食公司| 亚洲美女搞黄在线观看| 午夜免费男女啪啪视频观看| 大陆偷拍与自拍| av网站免费在线观看视频 | 校园人妻丝袜中文字幕| 国内精品宾馆在线| 蜜桃亚洲精品一区二区三区| 国产成人精品一,二区| 久久精品国产亚洲av涩爱| 免费观看av网站的网址| 亚洲熟女精品中文字幕| 久久久久久国产a免费观看| 精品一区二区三区视频在线| 免费观看精品视频网站| 久久久久久久久中文| 小蜜桃在线观看免费完整版高清| 成人亚洲欧美一区二区av| 国产精品福利在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 婷婷六月久久综合丁香| 久久久久久国产a免费观看| 性色avwww在线观看| 国产乱人偷精品视频| 国产精品.久久久| 久久99热这里只有精品18| 国产成人a∨麻豆精品| 国产久久久一区二区三区| 午夜免费观看性视频| 寂寞人妻少妇视频99o| 亚洲欧美成人综合另类久久久| 日韩强制内射视频| 性色avwww在线观看| 99热全是精品| 久久久久久久久久久丰满| 日韩电影二区| 人人妻人人澡人人爽人人夜夜 | 亚洲国产欧美在线一区| 亚洲av一区综合| 亚洲人与动物交配视频| 十八禁国产超污无遮挡网站| 免费大片黄手机在线观看| 22中文网久久字幕| 九九久久精品国产亚洲av麻豆| 日韩强制内射视频| 亚洲欧美日韩卡通动漫| 能在线免费观看的黄片| 日韩欧美国产在线观看| 高清视频免费观看一区二区 | 禁无遮挡网站| 免费看日本二区| 欧美高清成人免费视频www| 熟女电影av网| 少妇被粗大猛烈的视频| a级毛片免费高清观看在线播放| 六月丁香七月| 国产69精品久久久久777片| 免费观看精品视频网站| 九九在线视频观看精品| 日韩中字成人| 乱人视频在线观看| 午夜福利视频1000在线观看| 晚上一个人看的免费电影| 成人午夜精彩视频在线观看| 午夜福利在线观看免费完整高清在| 3wmmmm亚洲av在线观看| 淫秽高清视频在线观看| 久久久色成人| 亚洲精品一区蜜桃| 成人高潮视频无遮挡免费网站| 日本与韩国留学比较| av一本久久久久| 亚洲国产精品成人久久小说| 国产亚洲av片在线观看秒播厂 | 成年女人在线观看亚洲视频 | 午夜精品一区二区三区免费看| a级一级毛片免费在线观看| 精品久久久久久电影网| 狠狠精品人妻久久久久久综合| 丝袜喷水一区| 亚洲成色77777| 男的添女的下面高潮视频| 亚洲欧洲国产日韩| 国国产精品蜜臀av免费| 天天躁日日操中文字幕| 久久人人爽人人片av| 精品人妻一区二区三区麻豆| 色5月婷婷丁香| 国产极品天堂在线| 国产精品爽爽va在线观看网站| 午夜激情福利司机影院| 99热这里只有是精品50| 不卡视频在线观看欧美| 纵有疾风起免费观看全集完整版 | 日本-黄色视频高清免费观看| 亚洲国产精品sss在线观看| 一边亲一边摸免费视频| 亚洲国产最新在线播放| 亚洲国产成人一精品久久久| 久久久欧美国产精品| 久久久久久久久久久免费av| 日本与韩国留学比较| 久久久精品欧美日韩精品| 午夜激情福利司机影院| 免费观看精品视频网站| 日韩人妻高清精品专区| 欧美另类一区| 非洲黑人性xxxx精品又粗又长| 国产探花在线观看一区二区| 麻豆久久精品国产亚洲av| 国产国拍精品亚洲av在线观看| 老司机影院成人| 久久久a久久爽久久v久久| 全区人妻精品视频| 国产成人免费观看mmmm| 91精品伊人久久大香线蕉| 99久国产av精品| 亚洲色图av天堂| 亚洲伊人久久精品综合| 国产在视频线精品| 精品亚洲乱码少妇综合久久| 亚洲国产精品成人综合色| 中文字幕免费在线视频6| 国产又色又爽无遮挡免| 天天躁夜夜躁狠狠久久av| 亚洲美女视频黄频| 欧美日韩精品成人综合77777| 3wmmmm亚洲av在线观看| 欧美激情久久久久久爽电影| 高清在线视频一区二区三区| 国产在视频线精品| 国产精品日韩av在线免费观看| 一级毛片久久久久久久久女| 久99久视频精品免费| 免费看不卡的av| 男人舔奶头视频| 欧美一区二区亚洲| 亚洲国产成人一精品久久久| 少妇高潮的动态图| 亚洲av.av天堂| 非洲黑人性xxxx精品又粗又长| 天堂俺去俺来也www色官网 | 亚洲精品自拍成人| 超碰av人人做人人爽久久| 国产亚洲一区二区精品| 亚洲最大成人中文| 精品国产三级普通话版| 午夜老司机福利剧场| 国产一区有黄有色的免费视频 | 亚洲伊人久久精品综合| 高清欧美精品videossex| 免费少妇av软件| 亚洲精品aⅴ在线观看| av一本久久久久| 特大巨黑吊av在线直播| 黄色一级大片看看| 菩萨蛮人人尽说江南好唐韦庄| 18禁裸乳无遮挡免费网站照片| 乱码一卡2卡4卡精品| 日本-黄色视频高清免费观看| 五月玫瑰六月丁香| 国产精品一区www在线观看| 97超碰精品成人国产| 国产午夜精品久久久久久一区二区三区| 免费看美女性在线毛片视频| 免费看av在线观看网站| 天天一区二区日本电影三级| 日本wwww免费看| 亚洲av免费高清在线观看| 午夜免费观看性视频| 汤姆久久久久久久影院中文字幕 | 日本黄色片子视频| 欧美激情久久久久久爽电影| 久久久成人免费电影| 最近手机中文字幕大全| 欧美潮喷喷水| 街头女战士在线观看网站| 在线观看免费高清a一片| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 欧美性感艳星| 亚洲精品国产av蜜桃| 熟女电影av网| 夫妻午夜视频| 国产单亲对白刺激| 国产毛片a区久久久久| 久久久成人免费电影| 国产精品久久视频播放| 国产大屁股一区二区在线视频| 国产成人精品久久久久久| 国产精品三级大全| 麻豆精品久久久久久蜜桃| 免费大片黄手机在线观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧洲日产国产| 久久久久网色| 亚洲精品一区蜜桃| 国产精品1区2区在线观看.| 亚洲18禁久久av| 亚洲国产色片| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久| 色播亚洲综合网| 高清在线视频一区二区三区| 人人妻人人澡人人爽人人夜夜 | 亚洲国产精品成人综合色| 又黄又爽又刺激的免费视频.| 1000部很黄的大片| 精品久久久久久久久av| 国产美女午夜福利| 亚洲性久久影院| 男女下面进入的视频免费午夜| 欧美成人一区二区免费高清观看| 亚洲精品一区蜜桃| 肉色欧美久久久久久久蜜桃 | 亚州av有码| 秋霞在线观看毛片| 日日摸夜夜添夜夜添av毛片| 草草在线视频免费看| 插阴视频在线观看视频| 啦啦啦韩国在线观看视频| 精品久久久久久久久亚洲| av国产久精品久网站免费入址| 成人欧美大片| 美女主播在线视频| 亚洲欧美清纯卡通| 欧美变态另类bdsm刘玥| 欧美bdsm另类| 国内精品宾馆在线| 欧美日本视频| 26uuu在线亚洲综合色| 日日摸夜夜添夜夜爱| 黄色欧美视频在线观看| 国产av在哪里看| 丰满少妇做爰视频| 欧美潮喷喷水| 精华霜和精华液先用哪个| 男插女下体视频免费在线播放| 久久亚洲国产成人精品v| 男女下面进入的视频免费午夜| 亚洲人与动物交配视频| 大陆偷拍与自拍| 视频中文字幕在线观看| 天堂网av新在线| 久久国产乱子免费精品| 禁无遮挡网站| 日本熟妇午夜| 精品久久久久久久久av| 高清毛片免费看| 中文字幕av成人在线电影| 国产成人a∨麻豆精品| 赤兔流量卡办理| 99久久中文字幕三级久久日本| 久久精品国产亚洲网站| 夜夜看夜夜爽夜夜摸| 日本免费a在线| 成人综合一区亚洲| 日韩精品有码人妻一区| 久久人人爽人人片av| 亚洲电影在线观看av| 欧美zozozo另类| 亚洲av中文字字幕乱码综合| 精品国产三级普通话版| 在线免费观看不下载黄p国产| 国产男女超爽视频在线观看| 狂野欧美激情性xxxx在线观看| 熟妇人妻不卡中文字幕| 久久久久久九九精品二区国产| 午夜福利视频精品| 久久久欧美国产精品| 最近最新中文字幕免费大全7| 欧美丝袜亚洲另类| 全区人妻精品视频| 日日干狠狠操夜夜爽| 一级毛片aaaaaa免费看小| 亚洲伊人久久精品综合| 麻豆精品久久久久久蜜桃| 国产精品一区二区性色av| 亚洲国产精品成人综合色| 欧美高清成人免费视频www| 乱码一卡2卡4卡精品| 网址你懂的国产日韩在线| 一级毛片电影观看| 在线观看人妻少妇| 3wmmmm亚洲av在线观看| 51国产日韩欧美| 一级毛片黄色毛片免费观看视频| 最近中文字幕高清免费大全6| 国产亚洲91精品色在线| 插阴视频在线观看视频| 亚洲最大成人中文| 狂野欧美激情性xxxx在线观看| 一级爰片在线观看| 国产一区二区三区av在线| 午夜爱爱视频在线播放| 夜夜看夜夜爽夜夜摸| 在线观看免费高清a一片| 国产精品不卡视频一区二区| 又粗又硬又长又爽又黄的视频| 久久精品国产鲁丝片午夜精品| 精品人妻视频免费看| 男人爽女人下面视频在线观看| 波野结衣二区三区在线| 青青草视频在线视频观看| 国产在线男女| 国产精品人妻久久久影院| 亚洲国产成人一精品久久久| 日产精品乱码卡一卡2卡三| av国产久精品久网站免费入址| 看免费成人av毛片| 舔av片在线| 青春草视频在线免费观看| 国产精品一及| 亚洲久久久久久中文字幕| 国产精品无大码| 久久人人爽人人片av| 边亲边吃奶的免费视频| 午夜福利在线观看吧| 高清视频免费观看一区二区 | 99热这里只有是精品50| 国产不卡一卡二| 亚洲精品456在线播放app| 免费高清在线观看视频在线观看| 在线免费十八禁| 免费观看在线日韩| 国产免费视频播放在线视频 | 深爱激情五月婷婷| 国产淫语在线视频| 亚州av有码| 女的被弄到高潮叫床怎么办| freevideosex欧美| 国产黄色视频一区二区在线观看| 亚洲av日韩在线播放| 高清欧美精品videossex| 精品国产露脸久久av麻豆 | 精品久久久久久久久久久久久| 精华霜和精华液先用哪个| 麻豆av噜噜一区二区三区| 亚洲精品日韩在线中文字幕| 中文乱码字字幕精品一区二区三区 | 日韩av在线免费看完整版不卡| 在现免费观看毛片| 国产av码专区亚洲av| 午夜福利在线观看吧| 九草在线视频观看| 欧美性感艳星| 波多野结衣巨乳人妻| 亚洲内射少妇av| av专区在线播放| 国产亚洲91精品色在线| 观看免费一级毛片| 亚洲伊人久久精品综合| 亚洲av在线观看美女高潮| 在线观看一区二区三区| 中文字幕制服av| 看非洲黑人一级黄片| 欧美丝袜亚洲另类| 午夜福利在线在线| 日韩 亚洲 欧美在线| 午夜福利高清视频| 好男人在线观看高清免费视频| 国产精品久久久久久久电影| 色播亚洲综合网| 岛国毛片在线播放| 国产av国产精品国产| 欧美激情国产日韩精品一区| 欧美性感艳星| 国产亚洲精品av在线| 国产精品日韩av在线免费观看| 亚洲第一区二区三区不卡| 成人性生交大片免费视频hd| 日韩人妻高清精品专区| 99热6这里只有精品| 中文欧美无线码| 九九爱精品视频在线观看| 99久久中文字幕三级久久日本| 亚洲在线观看片| 亚洲自偷自拍三级| 国产精品麻豆人妻色哟哟久久 | 中文字幕制服av| 亚洲成人精品中文字幕电影| 99热这里只有是精品在线观看| 久久99蜜桃精品久久| 三级毛片av免费| 久久精品人妻少妇| 少妇人妻精品综合一区二区| 日韩av不卡免费在线播放| 午夜激情福利司机影院| 男女边摸边吃奶| 青春草国产在线视频| 亚洲欧美日韩无卡精品| 超碰97精品在线观看| 国产一区二区在线观看日韩| 精品久久久久久成人av| 精品久久久久久电影网| 国产久久久一区二区三区| 99久久精品热视频| 高清在线视频一区二区三区| 日韩大片免费观看网站| 亚洲在线自拍视频| 欧美精品国产亚洲| 精品国内亚洲2022精品成人| 一级毛片电影观看| 国产成人一区二区在线| 日韩成人av中文字幕在线观看| 九草在线视频观看| videossex国产| 日本免费在线观看一区| 国产亚洲av嫩草精品影院| 亚洲av日韩在线播放| 精品人妻熟女av久视频| 亚洲av中文av极速乱| 女人被狂操c到高潮| 亚洲成人中文字幕在线播放| 亚洲av电影不卡..在线观看| 亚洲av成人精品一二三区| 精品久久久久久久久久久久久| 国产欧美另类精品又又久久亚洲欧美| 联通29元200g的流量卡| a级毛片免费高清观看在线播放| 97超碰精品成人国产| 国产综合精华液| 日韩国内少妇激情av| 久久韩国三级中文字幕| 精品国产露脸久久av麻豆 | 亚洲国产精品国产精品| 搡女人真爽免费视频火全软件| 别揉我奶头 嗯啊视频| 日韩大片免费观看网站| 午夜免费观看性视频| 欧美一区二区亚洲| 成人亚洲精品一区在线观看 | 美女被艹到高潮喷水动态| 久久久久久国产a免费观看| 在线播放无遮挡| 国产黄片视频在线免费观看| 人人妻人人澡人人爽人人夜夜 | 成人美女网站在线观看视频| 国产成人福利小说| 91狼人影院| 国产成人a∨麻豆精品|