• <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
    最新的欧美精品一区二区| 欧美少妇被猛烈插入视频| 国产淫语在线视频| 国产人伦9x9x在线观看 | 日本av手机在线免费观看| 亚洲av电影在线进入| 久久精品久久久久久久性| 亚洲美女黄色视频免费看| 久久人人97超碰香蕉20202| 国产片内射在线| 免费观看av网站的网址| 久久久久久人人人人人| av片东京热男人的天堂| 午夜免费男女啪啪视频观看| 一本久久精品| 亚洲综合色惰| 不卡视频在线观看欧美| 在线免费观看不下载黄p国产| 中国三级夫妇交换| 美女高潮到喷水免费观看| 麻豆乱淫一区二区| 日韩不卡一区二区三区视频在线| 人妻系列 视频| 日韩人妻精品一区2区三区| 精品人妻一区二区三区麻豆| 日韩伦理黄色片| 日韩视频在线欧美| 啦啦啦在线免费观看视频4| 久久人妻熟女aⅴ| 国产成人精品一,二区| 伦理电影免费视频| 国产一区二区三区综合在线观看| 日韩 亚洲 欧美在线| 亚洲精品,欧美精品| 亚洲图色成人| av国产精品久久久久影院| 最近中文字幕2019免费版| 国产乱来视频区| 午夜激情久久久久久久| 最黄视频免费看| 中文天堂在线官网| 国产xxxxx性猛交| 大片免费播放器 马上看| 男女边摸边吃奶| 亚洲精品久久成人aⅴ小说| 国产精品嫩草影院av在线观看| 国产成人91sexporn| 国产男女超爽视频在线观看| 亚洲成人一二三区av| 夜夜骑夜夜射夜夜干| 久久久久久伊人网av| 1024香蕉在线观看| 狂野欧美激情性bbbbbb| 精品一区在线观看国产| 国产成人午夜福利电影在线观看| 飞空精品影院首页| 免费日韩欧美在线观看| 国产一区二区 视频在线| 精品一区二区三区四区五区乱码 | 日本av手机在线免费观看| 亚洲精品,欧美精品| 欧美激情 高清一区二区三区| 国产精品国产三级国产专区5o| 看免费av毛片| 亚洲一级一片aⅴ在线观看| 日本免费在线观看一区| 晚上一个人看的免费电影| 久久久久久久久免费视频了| 最新的欧美精品一区二区| 高清视频免费观看一区二区| 久久久久久免费高清国产稀缺| 色婷婷av一区二区三区视频| 老司机影院毛片| 女人高潮潮喷娇喘18禁视频| 老司机亚洲免费影院| 国产一区二区三区av在线| 老汉色av国产亚洲站长工具| 999精品在线视频| 国产精品女同一区二区软件| 免费观看a级毛片全部| 在线观看免费视频网站a站| 国产一级毛片在线| 一级a爱视频在线免费观看| 性色avwww在线观看| av一本久久久久| 亚洲av中文av极速乱| 国产免费又黄又爽又色| 老司机亚洲免费影院| 久久久精品94久久精品| 亚洲熟女精品中文字幕| 日本黄色日本黄色录像| 熟妇人妻不卡中文字幕| 久久久久久久国产电影| 国产亚洲一区二区精品| 婷婷色综合www| 日本爱情动作片www.在线观看| 国产激情久久老熟女| 亚洲伊人久久精品综合| 日韩av不卡免费在线播放| 国产熟女欧美一区二区| 精品国产一区二区三区久久久樱花| 精品酒店卫生间| 日韩一区二区视频免费看| tube8黄色片| 一边摸一边做爽爽视频免费| 春色校园在线视频观看| 老司机亚洲免费影院| 亚洲三区欧美一区| 欧美日韩av久久| 国产日韩欧美视频二区| 欧美成人午夜精品| 免费看av在线观看网站| 免费黄网站久久成人精品| 男女午夜视频在线观看| 久久国产精品男人的天堂亚洲| 一级爰片在线观看| 乱人伦中国视频| 亚洲少妇的诱惑av| 亚洲第一av免费看| 另类亚洲欧美激情| 十八禁网站网址无遮挡| 精品国产露脸久久av麻豆| 国产成人精品无人区| 桃花免费在线播放| 最近手机中文字幕大全| 亚洲精品一区蜜桃| 国产精品av久久久久免费| 国产黄色免费在线视频| 18禁裸乳无遮挡动漫免费视频| 王馨瑶露胸无遮挡在线观看| 久久久久久人妻| 18禁国产床啪视频网站| 欧美日韩成人在线一区二区| 男女国产视频网站| 在线观看免费高清a一片| 国产在视频线精品| 精品亚洲成a人片在线观看| 秋霞伦理黄片| 国产精品 国内视频| 在线观看三级黄色| 国产精品嫩草影院av在线观看| 涩涩av久久男人的天堂| 久久久久久久久久久免费av| 视频区图区小说| 又粗又硬又长又爽又黄的视频| 国产麻豆69| 香蕉丝袜av| 久久人人爽av亚洲精品天堂| 丰满少妇做爰视频| 欧美 日韩 精品 国产| 另类亚洲欧美激情| 日产精品乱码卡一卡2卡三| 国产亚洲欧美精品永久| 热99国产精品久久久久久7| 97在线人人人人妻| 国产高清不卡午夜福利| 精品久久蜜臀av无| www日本在线高清视频| 国产亚洲av片在线观看秒播厂| 国产精品.久久久| 精品国产乱码久久久久久小说| 久久99热这里只频精品6学生| 建设人人有责人人尽责人人享有的| 成年动漫av网址| 午夜福利在线观看免费完整高清在| 亚洲av成人精品一二三区| 9191精品国产免费久久| 18禁观看日本| 黄频高清免费视频| 日韩三级伦理在线观看| 一二三四中文在线观看免费高清| 国产精品久久久久久久久免| 欧美日韩综合久久久久久| 亚洲欧洲国产日韩| 夫妻午夜视频| 国产精品国产三级国产专区5o| 日韩制服丝袜自拍偷拍| 久久女婷五月综合色啪小说| 七月丁香在线播放| 在现免费观看毛片| 国产精品国产av在线观看| 麻豆av在线久日| 午夜久久久在线观看| 男人爽女人下面视频在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲在久久综合| 国产精品偷伦视频观看了| 黑人欧美特级aaaaaa片| 啦啦啦啦在线视频资源| 伦精品一区二区三区| 欧美日韩视频精品一区| 精品卡一卡二卡四卡免费| 99久久人妻综合| 亚洲人成77777在线视频| 极品人妻少妇av视频| 久久午夜福利片| 乱人伦中国视频| 精品国产一区二区三区久久久樱花| 欧美日韩成人在线一区二区| 免费久久久久久久精品成人欧美视频| 激情视频va一区二区三区| 午夜福利乱码中文字幕| 亚洲成av片中文字幕在线观看 | 久久 成人 亚洲| www日本在线高清视频| 在线 av 中文字幕| 久久久久国产精品人妻一区二区| 欧美日韩av久久| www.自偷自拍.com| 国产精品久久久久久精品古装| 国产精品麻豆人妻色哟哟久久| 欧美+日韩+精品| 亚洲国产日韩一区二区| 一区二区三区精品91| 久久久a久久爽久久v久久| 久久女婷五月综合色啪小说| 亚洲综合色惰| 97人妻天天添夜夜摸| 午夜福利一区二区在线看| 美女大奶头黄色视频| 午夜久久久在线观看| 侵犯人妻中文字幕一二三四区| 久久久欧美国产精品| 久久99精品国语久久久| 免费观看无遮挡的男女| 三上悠亚av全集在线观看| 成人午夜精彩视频在线观看| 亚洲,欧美,日韩| 国产白丝娇喘喷水9色精品| 电影成人av| 精品久久蜜臀av无| 亚洲国产精品成人久久小说| 美女午夜性视频免费| 人人妻人人爽人人添夜夜欢视频| 久久人妻熟女aⅴ| 久久精品国产自在天天线| 有码 亚洲区| 久久 成人 亚洲| 卡戴珊不雅视频在线播放| av电影中文网址| 久久久久久久大尺度免费视频| 十八禁高潮呻吟视频| 激情视频va一区二区三区| 亚洲视频免费观看视频| 大香蕉久久成人网| 久久 成人 亚洲| 波多野结衣av一区二区av| 夫妻性生交免费视频一级片| 两性夫妻黄色片| av一本久久久久| 老司机亚洲免费影院| 日韩中文字幕欧美一区二区 | 亚洲欧洲日产国产| 丝袜脚勾引网站| 1024香蕉在线观看| 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 国产亚洲av片在线观看秒播厂| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| 国产亚洲一区二区精品| 伦理电影大哥的女人| 一本色道久久久久久精品综合| 99久久综合免费| 校园人妻丝袜中文字幕| 超碰97精品在线观看| av.在线天堂| 美女主播在线视频| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美清纯卡通| 女的被弄到高潮叫床怎么办| 亚洲成人av在线免费| 亚洲av成人精品一二三区| 日本欧美国产在线视频| 新久久久久国产一级毛片| 成人毛片60女人毛片免费| 久久热在线av| 韩国精品一区二区三区| freevideosex欧美| 久久午夜福利片| 精品一区二区三卡| 在线观看免费日韩欧美大片| 国产免费视频播放在线视频| 精品一区二区三卡| 亚洲天堂av无毛| 午夜91福利影院| 丝袜喷水一区| 中文字幕另类日韩欧美亚洲嫩草| 多毛熟女@视频| 男女午夜视频在线观看| 看非洲黑人一级黄片| 欧美成人午夜精品| 视频在线观看一区二区三区| 9热在线视频观看99| 一级毛片黄色毛片免费观看视频| 亚洲国产欧美网| 久久久久精品人妻al黑| 国产片内射在线| 日韩视频在线欧美| 亚洲av男天堂| 男女午夜视频在线观看| 色婷婷av一区二区三区视频| 亚洲伊人色综图| 美女主播在线视频| 女人被躁到高潮嗷嗷叫费观| 日本av免费视频播放| 久久久久国产精品人妻一区二区| 久久久亚洲精品成人影院| 亚洲综合精品二区| 夫妻性生交免费视频一级片| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| 综合色丁香网| 久久精品国产亚洲av涩爱| 国产精品欧美亚洲77777| 亚洲欧美色中文字幕在线| 人人妻人人澡人人看| 婷婷色av中文字幕| 国产一区二区三区av在线| 亚洲欧美日韩另类电影网站| 美女大奶头黄色视频| 不卡av一区二区三区| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 熟女av电影| 又粗又硬又长又爽又黄的视频| 成人二区视频| 日本-黄色视频高清免费观看| 色视频在线一区二区三区| 亚洲精品国产av蜜桃| 天堂中文最新版在线下载| 欧美精品av麻豆av| 九九爱精品视频在线观看| 中文字幕人妻丝袜制服| 精品第一国产精品| 亚洲精品一区蜜桃| 婷婷色麻豆天堂久久| 亚洲国产欧美网| 亚洲国产看品久久| 69精品国产乱码久久久| 国产精品久久久久久精品电影小说| 国产探花极品一区二区| 日韩精品有码人妻一区| 97人妻天天添夜夜摸| 五月开心婷婷网| 国产亚洲午夜精品一区二区久久| 狠狠精品人妻久久久久久综合| 欧美精品一区二区大全| 亚洲色图综合在线观看| 美女中出高潮动态图| 美女主播在线视频| 女性被躁到高潮视频| 99re6热这里在线精品视频| 亚洲国产欧美在线一区| 美女中出高潮动态图| 国产一区亚洲一区在线观看| 男人添女人高潮全过程视频| √禁漫天堂资源中文www| 国产一区有黄有色的免费视频| 日韩中字成人| 最近的中文字幕免费完整| 熟女电影av网| av一本久久久久| 亚洲av在线观看美女高潮| 欧美精品人与动牲交sv欧美| 一区二区三区四区激情视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成电影观看| 一区二区三区乱码不卡18| 国产伦理片在线播放av一区| 日韩大片免费观看网站| 婷婷色av中文字幕| 满18在线观看网站| 黄色怎么调成土黄色| 国产xxxxx性猛交| 欧美日韩一级在线毛片| 日韩 亚洲 欧美在线| 久久久国产精品麻豆| 欧美国产精品一级二级三级| 亚洲美女搞黄在线观看| 亚洲精品久久久久久婷婷小说| 一区二区三区乱码不卡18| 人妻一区二区av| 中文字幕色久视频| 晚上一个人看的免费电影| 日韩一本色道免费dvd| 下体分泌物呈黄色| 一级,二级,三级黄色视频| 免费黄色在线免费观看| 亚洲美女黄色视频免费看| 性色av一级| 亚洲熟女精品中文字幕| 亚洲美女搞黄在线观看| 亚洲一级一片aⅴ在线观看| 乱人伦中国视频| 少妇人妻 视频| 精品99又大又爽又粗少妇毛片| 中文字幕人妻熟女乱码| 高清在线视频一区二区三区| 桃花免费在线播放| 亚洲第一av免费看| 啦啦啦在线观看免费高清www| 精品视频人人做人人爽| 日韩av免费高清视频| 成年动漫av网址| 激情五月婷婷亚洲| 91精品三级在线观看| 国精品久久久久久国模美| 女人久久www免费人成看片| 男男h啪啪无遮挡| 欧美日韩综合久久久久久| 又黄又粗又硬又大视频| 国产精品 国内视频| 韩国av在线不卡| 久久99精品国语久久久| 三上悠亚av全集在线观看| 精品久久久久久电影网| 国产乱人偷精品视频| 各种免费的搞黄视频| 国产成人精品一,二区| 有码 亚洲区| 国产不卡av网站在线观看| 青青草视频在线视频观看| 久久鲁丝午夜福利片| 啦啦啦在线免费观看视频4| 波野结衣二区三区在线| 色婷婷久久久亚洲欧美| 国产精品国产三级专区第一集| 免费看av在线观看网站| 亚洲国产av新网站| 亚洲国产毛片av蜜桃av| 欧美中文综合在线视频| 欧美+日韩+精品| 卡戴珊不雅视频在线播放| 人体艺术视频欧美日本| 97人妻天天添夜夜摸| 欧美日韩成人在线一区二区| 午夜精品国产一区二区电影| 国产成人免费观看mmmm| h视频一区二区三区| 成人18禁高潮啪啪吃奶动态图| 国产探花极品一区二区| 亚洲国产最新在线播放| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 免费观看a级毛片全部| 国产在线一区二区三区精| 久久国产精品大桥未久av| 99九九在线精品视频| 国产白丝娇喘喷水9色精品| av国产精品久久久久影院| 男女午夜视频在线观看| 亚洲人成网站在线观看播放| 免费在线观看黄色视频的| 97人妻天天添夜夜摸| 香蕉丝袜av| 国产免费视频播放在线视频| 深夜精品福利| 久久久久人妻精品一区果冻| 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区| 黄色视频在线播放观看不卡| 精品国产乱码久久久久久男人| 精品少妇内射三级| 国产精品不卡视频一区二区| 男女边吃奶边做爰视频| 女性被躁到高潮视频| 久久鲁丝午夜福利片| 人妻人人澡人人爽人人| 久久久亚洲精品成人影院| 999精品在线视频| 久久久久精品久久久久真实原创| 国产日韩欧美在线精品| 两个人看的免费小视频| 成人亚洲欧美一区二区av| 国产片特级美女逼逼视频| 在线观看免费日韩欧美大片| 乱人伦中国视频| 2018国产大陆天天弄谢| 91精品三级在线观看| 欧美日韩国产mv在线观看视频| 人成视频在线观看免费观看| 在线观看www视频免费| 亚洲国产av影院在线观看| 免费少妇av软件| 青春草亚洲视频在线观看| 五月开心婷婷网| 国产精品熟女久久久久浪| 欧美亚洲 丝袜 人妻 在线| 18禁动态无遮挡网站| 亚洲伊人久久精品综合| 午夜福利一区二区在线看| 久久97久久精品| 成年美女黄网站色视频大全免费| 国产日韩欧美视频二区| 一区在线观看完整版| 国产乱人偷精品视频| 午夜激情久久久久久久| 国产毛片在线视频| 久久精品国产综合久久久| a级毛片在线看网站| 久久久久久人人人人人| 亚洲一级一片aⅴ在线观看| 久久婷婷青草| 中文字幕av电影在线播放| 免费观看在线日韩| 性色av一级| 国产成人免费观看mmmm| 亚洲欧美成人精品一区二区| 久久精品国产综合久久久| 在线看a的网站| 国产精品久久久久久精品古装| 十分钟在线观看高清视频www| av视频免费观看在线观看| 亚洲 欧美一区二区三区| 肉色欧美久久久久久久蜜桃| 久久精品久久久久久久性| 国产成人免费观看mmmm| 亚洲第一青青草原| 人妻人人澡人人爽人人| 免费人妻精品一区二区三区视频| 亚洲精品一区蜜桃| 国产亚洲最大av| 人人妻人人澡人人爽人人夜夜| 国产片特级美女逼逼视频| 免费观看a级毛片全部| 国产日韩欧美在线精品| 久久久久久久精品精品| 最近最新中文字幕大全免费视频 | 久久影院123| 黑丝袜美女国产一区| 侵犯人妻中文字幕一二三四区| 亚洲av成人精品一二三区| 久久精品aⅴ一区二区三区四区 | 久久久欧美国产精品| 亚洲国产欧美网| 亚洲精品日韩在线中文字幕| 亚洲伊人色综图| 欧美精品av麻豆av| 99国产综合亚洲精品| 亚洲精品成人av观看孕妇| 亚洲美女视频黄频| av视频免费观看在线观看| 99久久综合免费| 午夜日本视频在线| 亚洲少妇的诱惑av| 精品第一国产精品| 天堂中文最新版在线下载| 国产精品久久久久久精品古装| 日韩中文字幕视频在线看片| 精品少妇黑人巨大在线播放| 亚洲四区av| 国产一区二区三区av在线| 天天影视国产精品| 黄色怎么调成土黄色| 久久99一区二区三区| 国产欧美日韩一区二区三区在线| 国产亚洲午夜精品一区二区久久| 亚洲经典国产精华液单| 一边摸一边做爽爽视频免费| 一本—道久久a久久精品蜜桃钙片| 国产成人一区二区在线| 色婷婷久久久亚洲欧美| 免费在线观看黄色视频的| 人妻 亚洲 视频| 精品卡一卡二卡四卡免费| 18在线观看网站| 久久久a久久爽久久v久久| 我的亚洲天堂| 亚洲精品aⅴ在线观看| 久久人人97超碰香蕉20202| 国产片内射在线| 熟女电影av网| 丝袜在线中文字幕| 国产精品三级大全| 夫妻性生交免费视频一级片| 激情视频va一区二区三区| 亚洲av福利一区| 亚洲中文av在线| 亚洲精品成人av观看孕妇| 中国三级夫妇交换| 亚洲精品美女久久av网站| 亚洲国产色片| 精品久久久久久电影网| 人人妻人人爽人人添夜夜欢视频| 女性被躁到高潮视频| 视频区图区小说| 国产精品不卡视频一区二区| 99久国产av精品国产电影| 日本午夜av视频| 两个人免费观看高清视频| 精品卡一卡二卡四卡免费| 国产精品偷伦视频观看了| 欧美97在线视频| 99久久综合免费| 黄片无遮挡物在线观看| 国产xxxxx性猛交| 国产成人a∨麻豆精品| 午夜日本视频在线| 国产精品麻豆人妻色哟哟久久| 高清不卡的av网站| 2022亚洲国产成人精品| 亚洲国产毛片av蜜桃av| 黄色毛片三级朝国网站| 国产极品天堂在线| 巨乳人妻的诱惑在线观看| 91精品伊人久久大香线蕉| 国产精品偷伦视频观看了| 青春草视频在线免费观看| 嫩草影院入口| 国产一区二区三区av在线| 久久精品aⅴ一区二区三区四区 | 欧美激情极品国产一区二区三区| 在线免费观看不下载黄p国产| 爱豆传媒免费全集在线观看|