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

    修正金屬本構(gòu)模型在超高速撞擊模擬中的應(yīng)用*

    2022-10-10 02:24:50李名銳陳春林沈子楷
    爆炸與沖擊 2022年9期
    關(guān)鍵詞:超高速破片靶板

    馬 坤,李名銳,陳春林,沈子楷,2,周 剛

    (1. 西北核技術(shù)研究所,陜西 西安 710024;2. 清華大學(xué)工程物理系,北京 100084)

    柱形彈體超高速撞擊靶板研究是以空間碎片防護(hù)、反裝甲等為背景的研究工作。數(shù)值模擬是相關(guān)研究手段之一,主要通過(guò)離散近似的方式求解帶有初邊值條件、材料廣義應(yīng)力應(yīng)變關(guān)系的偏微分方程組,進(jìn)而獲得求解域內(nèi)的力學(xué)參量和動(dòng)力學(xué)行為。結(jié)果的準(zhǔn)確性主要取決于離散求解的算法以及材料的廣義應(yīng)力應(yīng)變關(guān)系(即本構(gòu)模型)。在算法方面,粒子類無(wú)網(wǎng)格法在解決超高速撞擊下材料大變形問(wèn)題優(yōu)勢(shì)明顯,目前光滑粒子流體動(dòng)力學(xué)(smooth particle hydrodynamics,SPH)方法應(yīng)用較為廣泛。而在本構(gòu)模型方面,如何更加準(zhǔn)確地反映超高速撞擊過(guò)程中材料的動(dòng)力學(xué)行為,是力學(xué)工作者主要的研究方向,具有重要的學(xué)術(shù)意義和應(yīng)用價(jià)值。

    金屬材料本構(gòu)模型一般包括物態(tài)方程、強(qiáng)度模型和失效模型。在爆炸與沖擊動(dòng)力學(xué)問(wèn)題的數(shù)值模擬中,金屬材料普遍采用Grüneisen 狀態(tài)方程。對(duì)于超高速撞擊問(wèn)題,為了能清楚地獲得材料的相態(tài)情況,唐蜜構(gòu)建了GRAY 三相完全物態(tài)方程并開(kāi)展了相關(guān)研究。在強(qiáng)度模型方面,Johnson-Cook 模型和Steinberg 模型是超高速撞擊數(shù)值模擬中普遍使用的強(qiáng)度模型。兩個(gè)模型各有特點(diǎn),Johnson-Cook 模型基于一維應(yīng)力實(shí)驗(yàn)得到,Steinberg 模型基于一維應(yīng)變實(shí)驗(yàn)得到。在超高速撞擊后期的靶板慣性擴(kuò)孔階段,材料的壓力、應(yīng)變率開(kāi)始下降,Johnson-Cook 強(qiáng)度模型能很好地描述此時(shí)的材料動(dòng)力學(xué)行為。在材料失效模型方面,Johnson-Cook 失效模型在超高速撞擊過(guò)程中也有較為廣泛的應(yīng)用。但Johnson-Cook 失效模型在計(jì)算壓力遠(yuǎn)大于屈服強(qiáng)度的層裂問(wèn)題上,并不能準(zhǔn)確計(jì)算材料的層裂失效過(guò)程。Collé等在材料層裂問(wèn)題中采用當(dāng)壓力超過(guò)失效閾值后偏應(yīng)力歸零而壓力不變的方法建立失效模型,處理方法不符合物理實(shí)際。Grady 失效模型本質(zhì)上是拉伸失效模型,能夠反映材料層裂失效問(wèn)題,廣泛應(yīng)用于超高速撞擊數(shù)值模擬中。Eftis 等基于微孔洞的成核和長(zhǎng)大過(guò)程建立了損傷演化模型,并模擬了球形彈體超高速撞擊問(wèn)題,獲得的靶板的層裂以及斷口破壞形貌與實(shí)驗(yàn)結(jié)果一致。該模型主要考慮的是靜水壓力產(chǎn)生的層裂問(wèn)題,沒(méi)有考慮剪切失效機(jī)制。封加波損傷演化模型是具有明確物理意義的描述材料內(nèi)部微孔洞總體積演化導(dǎo)致層裂問(wèn)題的材料模型,一般用于模擬飛片撞擊、滑移內(nèi)爆(或外爆)加載柱殼破壞等問(wèn)題,在超高速撞擊問(wèn)題中的應(yīng)用較少。Seisson 等在數(shù)值模擬鋼球超高速撞擊石墨成坑問(wèn)題中,鋼球材料的失效模型同時(shí)考慮了最大拉伸主應(yīng)力和等效塑性變形兩種失效準(zhǔn)則。Corbett在計(jì)算鋁球超高速撞擊鋁靶問(wèn)題中,彈體采用了最大拉應(yīng)力以及最大剪應(yīng)力失效模型,靶板選用了Johnson-Cook 失效模型。Randles 等建立了拉伸損傷和Johnson-Cook 損傷的聯(lián)合損傷演化模型,并驗(yàn)證了數(shù)值模擬的準(zhǔn)確性。因此,同時(shí)考慮拉伸和剪切失效是全面且準(zhǔn)確的。但還需要注意的是:在大多數(shù)的數(shù)值模擬中,當(dāng)材料滿足失效準(zhǔn)則后,材料立刻失效,相應(yīng)的單元/粒子的屈服強(qiáng)度立刻為零并且不再承受拉伸應(yīng)力。這將產(chǎn)生大量互相無(wú)作用力的單獨(dú)單元/粒子,這和彈靶撞擊形成破片群的物理現(xiàn)象不符。在韌性金屬?zèng)_擊碎裂問(wèn)題方面,曹祥等以韌性金屬長(zhǎng)桿動(dòng)態(tài)拉伸碎裂為例,分析了非線性損傷演化過(guò)程對(duì)韌性碎裂過(guò)程的影響規(guī)律。該斷裂演化模型應(yīng)用于超高速撞擊數(shù)值模擬中,將有助于相對(duì)準(zhǔn)確地計(jì)算彈靶材料的碎裂過(guò)程。

    在現(xiàn)有的數(shù)值模擬研究中,失效模型還不能全面反映超高速撞擊下材料在層裂、剪切斷裂時(shí)的失效機(jī)制,并給出失效后的應(yīng)力歸零的斷裂演化過(guò)程。因此,本文中擬針對(duì)超高速撞擊問(wèn)題,在材料本構(gòu)模型方面,引入GRAY 三相物態(tài)方程以描述材料相態(tài)變化,采用Johnson-Cook 強(qiáng)度模型描述靶板后期動(dòng)態(tài)擴(kuò)孔,結(jié)合封加波損傷演化模型以及Johnson-Cook 失效模型,描述不同應(yīng)力三軸度下材料的層裂、剪切失效行為,并引入曹祥斷裂演化模型,描述材料失效后應(yīng)力歸零的過(guò)程。并將該修正金屬本構(gòu)模型應(yīng)用于柱形93 鎢超高速撞擊Q345 鋼板數(shù)值模擬中。通過(guò)與實(shí)驗(yàn)結(jié)果的對(duì)比,驗(yàn)證本構(gòu)模型的適用性和準(zhǔn)確性。初步分析典型彈靶撞擊條件下破片群的空間分布特征。

    1 修正金屬本構(gòu)模型簡(jiǎn)介

    1.1 GRAY 物態(tài)方程

    GRAY 物態(tài)方程通過(guò)初始熔化內(nèi)能、 完全熔化內(nèi)能、恒定容比熱初始內(nèi)能,將材料的固相-液相區(qū)分為固相、熔化相、液相以及熱液相。計(jì)算中可根據(jù)內(nèi)能判斷材料所處相態(tài),進(jìn)而給出相應(yīng)的修正壓力以及溫度,具體表達(dá)式如下。

    當(dāng)≤時(shí):

    1.2 Johnson-Cook 強(qiáng)度模型

    1.3 Johnson-Cook 失效模型

    Johnson-Cook 失效模型主要包含了應(yīng)力三軸度、應(yīng)變率以及溫度對(duì)材料失效的影響。失效應(yīng)變的表達(dá)式為:

    為了避免材料失效后應(yīng)力立刻歸零,引入曹祥提出的非線性斷裂演化模型。該模型以線性內(nèi)聚力斷裂假設(shè)Mott-Grady 模型為基礎(chǔ),描述了材料在滿足失效準(zhǔn)則后的斷裂(裂紋軟化)過(guò)程。即材料在達(dá)到失效閾值后,材料應(yīng)力存在一個(gè)逐漸卸載至零的過(guò)程。此時(shí)損傷屈服強(qiáng)度與內(nèi)聚力斷裂假設(shè)下的損傷因子、材料初始失效時(shí)的等效應(yīng)力的關(guān)系為:

    損傷因子假設(shè)為:

    式中:為材料斷裂過(guò)程中消耗的表面能,即材料的斷裂能。因此,材料的損傷演化過(guò)程主要取決于材料斷裂能、非線性系數(shù)以及單元的尺寸。

    1.4 封加波損傷演化模型

    封加波等將材料宏觀微元中的微損傷(微裂紋、微孔洞)總體積與宏觀微元體積之比定義為損傷度函數(shù)。基于Griffith-Orowan 提出的裂紋擴(kuò)展能量觀點(diǎn),假設(shè)微孔洞長(zhǎng)大過(guò)程中因表面積增加所需要的表面能由微孔洞周圍介質(zhì)所存儲(chǔ)的應(yīng)變能提供。此外,提供微孔洞動(dòng)態(tài)長(zhǎng)大的周圍介質(zhì)范圍由應(yīng)力波擾動(dòng)傳播的速度即材料聲速來(lái)確定。損傷度函數(shù)的演化方程為:

    式中:為拉伸狀態(tài)下基體材料的壓力,σ為孔洞長(zhǎng)大的初始應(yīng)力閾值,為基體材料的體積模量,ω 為表面功(斷裂比功)。

    在材料損傷的后期,微孔洞相互間開(kāi)始連通,材料的損傷快速增長(zhǎng),材料發(fā)生分離性斷裂,并且所承受的應(yīng)力在斷裂發(fā)生后下降到零。為近似描述材料損傷后期的應(yīng)力狀態(tài),仍采用曹祥等提出的斷裂演化模型,含損傷的本構(gòu)關(guān)系定義為:

    式中:為單元損傷啟動(dòng)時(shí)的靜水壓力。代表材料成核、長(zhǎng)大過(guò)程的損傷度的初始值為,代表材料微孔洞聚集并導(dǎo)致最終破壞的損傷函數(shù)的初始值為0。材料的損傷演化過(guò)程是:首先,為0 并且不變,由起始隨損傷演化逐漸增加,當(dāng)達(dá)到時(shí),值不再變化,開(kāi)始隨損傷演化而逐漸增加,當(dāng)增加到1 時(shí)表示材料完全層裂破壞。

    2 超高速撞擊數(shù)值模擬

    選用ANSYS AUTODYN 3D 中的光滑粒子流體動(dòng)力學(xué)(SPH)方法模擬超高速撞擊過(guò)程。根據(jù)軟件接口編寫修正金屬本構(gòu)模型的自定義子程序模塊,再通過(guò)數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果的對(duì)比來(lái)驗(yàn)證本構(gòu)模型的準(zhǔn)確性。選取柱形93 鎢彈體超高速撞擊Q345 靶板問(wèn)題進(jìn)行計(jì)算,材料模型參數(shù)見(jiàn)表1~3。GRAY 物態(tài)方程的其他材料參數(shù)可參考Royce的研究工作。

    表1 93 鎢合金、Q345 鋼GRAY 物態(tài)方程及損傷演化模型材料參數(shù)Table 1 The material parameters of the GRAY equation of state and the damage evolution model for 93 tungsten alloy and Q345 steel

    表2 93 鎢合金、Q345 鋼Johnson-Cook 強(qiáng)度模型材料參數(shù)Table 2 The material parameters of the Johnson-Cook yield criterion for 93 tungsten alloy and Q345 steel

    表3 93 鎢合金、Q345 鋼Johnson-Cook 失效模型材料參數(shù)Table 3 The material parameters of the Johnson-Cook failure model for 93 tungsten alloy and Q345 steel

    上述Johnson-Cook 強(qiáng)度模型、Johnson-Cook 失效模型以及封加波損傷演化模型對(duì)應(yīng)的材料參數(shù)中,除部分參數(shù)引用了相關(guān)文獻(xiàn)外,其他參數(shù)均為實(shí)驗(yàn)所得。包括SHPB 實(shí)驗(yàn)、缺口拉伸實(shí)驗(yàn)、飛片撞擊實(shí)驗(yàn)等。

    2.1 數(shù)值模擬與實(shí)驗(yàn)宏觀測(cè)量結(jié)果對(duì)比

    針對(duì)柱形93 鎢彈體超高速撞擊Q345 鋼板問(wèn)題,可通過(guò)實(shí)驗(yàn)后效以及高速攝像獲得靶板穿孔直徑、破片群擴(kuò)展速度、彈體剩余長(zhǎng)度的測(cè)量結(jié)果。因此,根據(jù)實(shí)驗(yàn)彈靶結(jié)構(gòu)和撞擊條件,建立數(shù)值模擬模型(光滑長(zhǎng)度為0.1 mm)并開(kāi)展計(jì)算。獲得的相應(yīng)參量的數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果一同列于表4~5 中。其中,v為破片群軸向最大速度,v為破片群橫向最大速度,Δ為彈體侵蝕長(zhǎng)度。需要說(shuō)明的是,實(shí)驗(yàn)1~6 中未使用高速攝像,實(shí)驗(yàn)8~9 中受彈托材料撞擊影響無(wú)法獲得靶板穿孔直徑數(shù)據(jù)。

    表4 靶板穿孔直徑的數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比Table 4 Comparison of numerical simulation results and experimental results of perforation diameter of the target plate

    表5 破片群擴(kuò)展速度、彈體侵蝕長(zhǎng)度數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果對(duì)比Table 5 Comparison of numerical simulation results and experimental results of fragment group expansion speed and erosion length of projectile

    對(duì)比表4~5 中數(shù)據(jù)可知,在靶板穿孔直徑、破片群擴(kuò)展速度、彈體剩余長(zhǎng)度方面,數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果一致性較好,大部分誤差在8%以內(nèi)。

    2.2 靶板層裂失效對(duì)破片群形貌的影響

    實(shí)驗(yàn)研究發(fā)現(xiàn),所有靶板厚度()小于0.72 倍彈體直徑()的實(shí)驗(yàn)中,彈靶撞擊后將在靶后破片群前部形成“尖端”。主要是因?yàn)閺棸凶矒舢a(chǎn)生向靶板內(nèi)部傳播的強(qiáng)沖擊波來(lái)不及被來(lái)自彈體頭部側(cè)向自由面、靶板正向自由面的卸載波及時(shí)追趕卸載,導(dǎo)致沖擊波在靶板背表面反射形成較強(qiáng)的拉伸波,進(jìn)而使靶板產(chǎn)生拉伸層裂破壞并形成速度較高的“尖端”破片。而在數(shù)值模擬中,彈靶材料模型包含反映材料層裂失效機(jī)制的損傷演化模型。因此,數(shù)值模擬應(yīng)能夠有效預(yù)測(cè)靶板是否產(chǎn)生與實(shí)驗(yàn)相同的層裂結(jié)果。根據(jù)實(shí)驗(yàn)條件建模并開(kāi)展數(shù)值模擬,對(duì)比靶后破片群形貌數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果,如表6所示。

    表6 不同靶厚與彈體直徑比的破片群實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果對(duì)比Table 6 Comparison of experimental results and numerical simulation results of fragment group with different target thickness and projectile diameter ratio

    從表6 可以看出,在靶后破片群前部是否形成“尖端”問(wèn)題上,數(shù)值模擬能夠獲得同實(shí)驗(yàn)一致的現(xiàn)象,這一定程度上驗(yàn)證了數(shù)值模擬中靶板材料模型所包含的反映材料層裂失效機(jī)制的損傷演化模型的準(zhǔn)確性。

    2.3 超高速撞擊下材料的熔化情況

    為驗(yàn)證修正金屬本構(gòu)模型中GRAY 物態(tài)方程計(jì)算材料相態(tài)的準(zhǔn)確性,數(shù)值模擬了彈體超高速撞擊雙層靶板問(wèn)題,分析了首層靶板穿孔及第2 層靶板破片撞擊坑的材料熔化情況。彈體尺寸為?3.5 mm×17.5 mm,撞擊速度為2.76 km/s,兩層靶板厚度均為1.5 mm,建立三維面對(duì)稱模型,光滑長(zhǎng)度為0.1 mm。

    彈體穿透首層靶板后,破片群形貌如圖1 所示,圖中標(biāo)記了材料所處相態(tài),藍(lán)色為固相,綠色為固液混合相,紅色為液相。彈體撞擊首層靶后,有部分粒子發(fā)生熔化,但總體上,破片群沒(méi)有發(fā)生大范圍的熔化。破片群“尖端”部分的破片均為固相。

    圖1 數(shù)值模擬首層靶后破片群形貌及熔化狀態(tài)Fig. 1 The morphology and melting state of the fragment group after the first layer of target obtained by numerical simulation

    對(duì)實(shí)驗(yàn)結(jié)果的分析表明,在該速度段下,首層靶后破片群沒(méi)有產(chǎn)生大范圍材料熔化。而當(dāng)剩余彈體及破片群撞擊第2 層靶板后,受二次加載作用影響,撞擊區(qū)附近將發(fā)生大范圍的材料熔化。典型實(shí)驗(yàn)獲得的首層、第2 層靶板中心穿孔斷口掃描電子顯微鏡(SEM)結(jié)果如圖2 所示,第1 層靶板斷口較為光滑,未發(fā)現(xiàn)大范圍熔化,而第2 層靶板斷口則有較多熔融后凝固組織。在圖1 的數(shù)值模擬中,能夠觀察到首層靶板斷口、整個(gè)破片群尚未發(fā)生大范圍的材料熔化。

    圖2 實(shí)驗(yàn)第1 層和第2 層靶板斷口的SEM 結(jié)果Fig. 2 SEM results of fracture surface of the first and second layers of the target plates

    在數(shù)值模擬中,當(dāng)剩余彈體穿透第2 層靶板時(shí),除靶板中心區(qū)域受剩余彈體撞擊外,其他區(qū)域均為破片群撞擊影響區(qū)。如圖3 所示,在破片群撞擊后的撞擊坑處,大量的SPH 粒子發(fā)生熔化。對(duì)比圖1 中撞擊前的破片群的熔化狀態(tài),可以判斷在圖3 中靶板破片撞擊坑中的SPH 粒子的熔化,大部分是破片在撞擊第2 層靶板過(guò)程中發(fā)生的。

    圖3 數(shù)值模擬彈體穿透第2 層靶板時(shí)彈靶位置及靶板撞擊坑熔化狀態(tài)Fig. 3 Numerical simulation of projectile target position and melting state of impact crater when projectile penetrates the second layer of the target plate

    上述關(guān)于材料熔化的數(shù)值模擬分析結(jié)果和實(shí)驗(yàn)中靶板破片撞擊坑微觀組織分析的結(jié)論是一致的,實(shí)驗(yàn)通過(guò)掃描電鏡獲得的破片撞擊坑的形貌如圖4 所示??梢钥闯觯矒艨颖砻鎴A潤(rùn)光滑,像液體澆筑后的表面形貌,說(shuō)明撞擊坑中的材料已經(jīng)發(fā)生了熔化。因此,認(rèn)為數(shù)值模擬基本能夠反映超高速撞擊過(guò)程中彈靶材料所處的相態(tài)。

    圖4 實(shí)驗(yàn)第2 層靶板破片撞擊坑SEM 形貌圖Fig. 4 SEM morphologies of impact crater of fragments in the second layer of target

    2.4 破片群對(duì)效應(yīng)靶的穿孔統(tǒng)計(jì)分析

    當(dāng)93 鎢合金彈體以2~3.5 km/s 速度撞擊鋼板時(shí),撞擊后彈靶材料破碎將形成平均尺寸較大的破片群。為了獲得破片群特征尺寸的分布規(guī)律,可在首層靶板后設(shè)置一層薄效應(yīng)靶。通過(guò)統(tǒng)計(jì)效應(yīng)靶穿孔面積與累積數(shù)量的關(guān)系,間接分析破片群特征尺寸分布規(guī)律以及材料的動(dòng)態(tài)破碎行為。因此,為了進(jìn)一步驗(yàn)證修正金屬本構(gòu)模型的計(jì)算結(jié)果,參照實(shí)驗(yàn)條件建模并開(kāi)展計(jì)算。彈體直徑15 mm、長(zhǎng)15 mm。首層靶板厚度10 mm,后效靶板厚度1 mm,兩層靶板間距1 014 mm。實(shí)驗(yàn)彈體速度為2.25 km/s,模型光滑長(zhǎng)度為0.25 mm。選取數(shù)值模擬彈體撞擊首層靶板后1 ms 時(shí)刻后效靶板的穿孔結(jié)果進(jìn)行分析,實(shí)驗(yàn)和數(shù)值模擬效應(yīng)靶穿孔形貌對(duì)比如圖5 所示。圖6 為效應(yīng)靶穿孔的累積數(shù)量與穿孔面積曲線的對(duì)比情況。其中,累積數(shù)量表示大于橫坐標(biāo)穿孔面積的穿孔總數(shù)。需要說(shuō)明的是,實(shí)驗(yàn)中,彈體撞擊首層靶板后形成的破片群經(jīng)歷了與空氣介質(zhì)短暫的熱量交換才撞擊到薄效應(yīng)靶上,而數(shù)值模擬中,則考慮短暫熱交換所引起的溫度變化對(duì)破片穿孔尺寸的影響相對(duì)破片密度以及速度對(duì)穿孔尺寸的影響是次要因素,因而忽略了空氣的影響。

    圖5 實(shí)驗(yàn)與數(shù)值模擬后效靶穿孔形貌對(duì)比Fig. 5 Comparison of experiment and numerical simulation on the perforation morphology of aftereffect target

    在效應(yīng)靶穿孔形貌上,由于實(shí)驗(yàn)難以保證絕對(duì)的軸對(duì)稱撞擊,并且材料破碎有一定的隨機(jī)性,導(dǎo)致實(shí)驗(yàn)效應(yīng)靶穿孔形貌對(duì)稱性不明顯,而數(shù)值模擬有明顯的對(duì)稱性。但總的來(lái)說(shuō),較大破片穿孔更趨近于效應(yīng)靶中心區(qū)域,效應(yīng)靶邊緣區(qū)域穿孔較小,這是實(shí)驗(yàn)以及數(shù)值模擬的共同特征。從圖6 可以看出,后效靶穿孔累積數(shù)量與穿孔面積統(tǒng)計(jì)曲線的數(shù)值模擬結(jié)果和實(shí)驗(yàn)結(jié)果一致性較好,僅在統(tǒng)計(jì)曲線兩端略有差距。那么從統(tǒng)計(jì)角度講,數(shù)值模擬獲得的超高速撞擊破片群符合實(shí)驗(yàn)破片群的特征尺寸分布規(guī)律。

    圖6 實(shí)驗(yàn)與數(shù)值模擬后效靶穿孔累積數(shù)量與穿孔面積統(tǒng)計(jì)曲線對(duì)比Fig. 6 Comparison of experiment and numerical simulation on the statistical curves of cumulative perforation number and perforation area

    2.5 典型彈靶撞擊下破片群的空間分布

    開(kāi)展典型93 鎢彈體超高速撞擊Q345 鋼板問(wèn)題的數(shù)值模擬,彈體直徑3.2 mm,長(zhǎng)徑比5,靶板厚度1.5 mm,撞擊速度為3.0 km/s,模型光滑長(zhǎng)度取0.1 mm。保存模擬結(jié)果中的SPH 粒子位置坐標(biāo),參照Liang 等所使用的廣度優(yōu)先搜索(BFS)算法,識(shí)別破片并統(tǒng)計(jì)分析相關(guān)物理量的空間分布情況。

    圖7~8 為破片群的數(shù)量和質(zhì)量在各軸向比速度段、橫向比速度段內(nèi)的累計(jì)情況。在每個(gè)柱形圖內(nèi),不同顏色代表不同質(zhì)量范圍破片的累積結(jié)果。v為破片軸向速度,為彈體初始速度,v為破片橫向速度,v為破片群的最大橫向速度。從圖7 可以看出,在軸向方向上破片群前端的破片數(shù)量高、總質(zhì)量占比高。而低速段的大質(zhì)量破片主要為靶板擴(kuò)孔末期剝落的低速破片。從圖8 中可以看出,在橫向方向上,各比速度段的數(shù)量基本相同。大質(zhì)量破片主要集中在低橫向速度段內(nèi),即破片群的前端和尾端。

    圖7 破片數(shù)量、質(zhì)量關(guān)于破片群軸向比速度的分布Fig. 7 Distribution of fragments number and mass in relation to fragment specific axial velocity

    圖8 破片數(shù)量、質(zhì)量關(guān)于破片群橫向比速度的分布Fig. 8 Distribution of fragments number and mass in relation to fragment specific lateral velocity

    圖9 為破片軸向動(dòng)量、橫向動(dòng)量(絕對(duì)值)在各軸向比速度段內(nèi)的累積結(jié)果。其中,不同顏色代表不同動(dòng)量范圍內(nèi)的累積情況??梢钥闯?,隨軸向比速度的增大,軸向動(dòng)量呈指數(shù)增長(zhǎng)形式分布,高軸向動(dòng)量破片集中分布于破片群前端。另一方面,高橫向動(dòng)量也集中于破片群的前端,結(jié)合破片群質(zhì)量分布特征認(rèn)為,這是由于該區(qū)域的破片群累積質(zhì)量較高。綜上,破片群的前端具有較高的質(zhì)量、軸向動(dòng)量以及橫向動(dòng)量(絕對(duì)值)。

    圖9 破片軸向動(dòng)量、橫向動(dòng)量關(guān)于破片群軸向比速度的分布Fig. 9 Distributions of fragments axial momentum and lateral momentum in relation to fragment specific axial velocity

    3 結(jié) 論

    為建立適用于超高速撞擊問(wèn)題的修正金屬本構(gòu)模型,引入GRAY 三相物態(tài)方程描述材料相態(tài)變化,采用Johnson-Cook 強(qiáng)度模型描述撞擊后期材料的力學(xué)行為,綜合封加波損傷演化模型以及Johnson-Cook 失效模型描述層裂、剪切失效機(jī)制,并引用曹祥斷裂演化模型,描述材料失效后應(yīng)力歸零的過(guò)程。通過(guò)柱形93 鎢超高速撞擊Q345 鋼板數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比,驗(yàn)證了本構(gòu)模型的適用性。分析了典型彈靶撞擊條件下破片群的空間分布特征,得出結(jié)論如下。

    (1)基于修正金屬本構(gòu)模型,開(kāi)展了不同撞擊條件的柱形93 鎢超高速撞擊Q345 鋼板數(shù)值模擬,獲得的靶板穿孔直徑、彈體侵蝕長(zhǎng)度、破片群擴(kuò)展速度結(jié)果與實(shí)驗(yàn)結(jié)果一致。

    (2)GRAY 三相物態(tài)方程能夠相對(duì)準(zhǔn)確地計(jì)算93 鎢合金撞擊Q345 鋼靶過(guò)程中彈靶材料所處的相態(tài)。數(shù)值模擬發(fā)現(xiàn),撞擊速度在2~3.5 km/s 時(shí),彈靶材料并不會(huì)發(fā)生大范圍熔化,但當(dāng)剩余彈體和破片群撞擊第2 層靶板時(shí),破片撞擊坑將產(chǎn)生大范圍熔化,該現(xiàn)象與實(shí)驗(yàn)結(jié)果是一致的。

    (3)封加波損傷演化模型能夠準(zhǔn)確得出超高速撞擊過(guò)程中靶板是否產(chǎn)生層裂破壞,對(duì)于不同靶板厚度與彈體直徑比值的撞擊問(wèn)題,數(shù)值模擬獲得的靶板是否產(chǎn)生層裂破壞的結(jié)果與實(shí)驗(yàn)結(jié)果一致。

    (4)綜合封加波損傷演化模型、Johnson-Cook 失效模型以及曹祥提出的斷裂演化模型后,典型條件下數(shù)值模擬獲得的破片群對(duì)后效靶板穿孔的面積與累積數(shù)量統(tǒng)計(jì)曲線和實(shí)驗(yàn)結(jié)果一致。

    (5)基于修正金屬本構(gòu)模型,獲得了典型條件下柱形93 鎢彈體超高速撞擊Q345 靶板破片群空間分布統(tǒng)計(jì)結(jié)果,破片群的前端具有較高的質(zhì)量、軸向動(dòng)量以及橫向動(dòng)量(絕對(duì)值)。

    猜你喜歡
    超高速破片靶板
    爆轟驅(qū)動(dòng)下軸向預(yù)制破片飛散特性研究
    艙壁組合結(jié)構(gòu)抗高速破片侵徹性能研究
    疊合雙層靶抗球形破片的侵徹能耗
    具有攻角的鎢合金彈侵徹運(yùn)動(dòng)靶板的數(shù)值模擬研究
    彈丸斜撞擊間隔靶板的數(shù)值模擬
    彈藥動(dòng)態(tài)加載下破片測(cè)試方法
    中國(guó)風(fēng)投行業(yè)迎來(lái)超高速發(fā)展
    金橋(2018年1期)2018-09-28 02:24:38
    采用超高速離心鑄造技術(shù)開(kāi)發(fā)耐磨薄壁鑄管
    四川冶金(2018年1期)2018-09-25 02:39:26
    半預(yù)制破片戰(zhàn)斗部破片威力分析
    基于ST-SRCKF的超高速?gòu)?qiáng)機(jī)動(dòng)目標(biāo)跟蹤算法
    岛国视频午夜一区免费看| 很黄的视频免费| 国产激情久久老熟女| 韩国av一区二区三区四区| 国产一区二区三区视频了| 国产高清激情床上av| 精品国产乱码久久久久久男人| 国产麻豆成人av免费视频| 91老司机精品| 欧美一区二区精品小视频在线| 亚洲,欧美精品.| 母亲3免费完整高清在线观看| 久久久国产成人精品二区| 欧美高清成人免费视频www| 精品国产亚洲在线| 18禁国产床啪视频网站| 国产精品 国内视频| 中文字幕人妻丝袜一区二区| 久9热在线精品视频| 男女那种视频在线观看| 妹子高潮喷水视频| 久久精品影院6| 久久久久国内视频| 天堂影院成人在线观看| 精品一区二区三区视频在线观看免费| 亚洲国产欧洲综合997久久,| 制服诱惑二区| 天堂av国产一区二区熟女人妻 | 巨乳人妻的诱惑在线观看| 90打野战视频偷拍视频| 最近在线观看免费完整版| 久久热在线av| 又紧又爽又黄一区二区| 亚洲激情在线av| 国产精品 国内视频| 欧美中文综合在线视频| 精品少妇一区二区三区视频日本电影| 深夜精品福利| 国产精华一区二区三区| 亚洲精品av麻豆狂野| 免费在线观看影片大全网站| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 亚洲中文日韩欧美视频| 久久国产乱子伦精品免费另类| 国产精品自产拍在线观看55亚洲| 可以免费在线观看a视频的电影网站| 91九色精品人成在线观看| 日韩欧美国产一区二区入口| 91av网站免费观看| 老司机午夜福利在线观看视频| 熟妇人妻久久中文字幕3abv| 午夜福利欧美成人| 18美女黄网站色大片免费观看| 国产精品 欧美亚洲| 亚洲国产精品合色在线| 国产蜜桃级精品一区二区三区| cao死你这个sao货| 观看免费一级毛片| 老熟妇仑乱视频hdxx| 哪里可以看免费的av片| 欧美人与性动交α欧美精品济南到| 久久久久久久久免费视频了| 欧美一区二区精品小视频在线| 国产三级黄色录像| 女人被狂操c到高潮| 香蕉久久夜色| 美女午夜性视频免费| 叶爱在线成人免费视频播放| 亚洲男人的天堂狠狠| 日韩精品中文字幕看吧| 18禁国产床啪视频网站| 18禁黄网站禁片午夜丰满| 波多野结衣巨乳人妻| 日韩免费av在线播放| 一二三四社区在线视频社区8| 欧美色视频一区免费| 久久热在线av| 999久久久国产精品视频| 波多野结衣高清作品| 午夜精品在线福利| 亚洲国产日韩欧美精品在线观看 | 动漫黄色视频在线观看| av有码第一页| 桃色一区二区三区在线观看| 婷婷亚洲欧美| 国产欧美日韩精品亚洲av| 亚洲一卡2卡3卡4卡5卡精品中文| 国产伦一二天堂av在线观看| 国产高清激情床上av| 又大又爽又粗| 麻豆一二三区av精品| 久久欧美精品欧美久久欧美| 亚洲七黄色美女视频| 亚洲av电影不卡..在线观看| 欧美不卡视频在线免费观看 | 啦啦啦韩国在线观看视频| 欧美一区二区精品小视频在线| 午夜福利免费观看在线| 成人三级做爰电影| 久久久久久免费高清国产稀缺| 国产高清视频在线观看网站| 最新在线观看一区二区三区| 777久久人妻少妇嫩草av网站| 国产av一区二区精品久久| 日韩高清综合在线| 精品久久久久久久人妻蜜臀av| 女人高潮潮喷娇喘18禁视频| 国产精品久久视频播放| 免费观看精品视频网站| 色在线成人网| 国产成人影院久久av| 91老司机精品| 成人三级黄色视频| 禁无遮挡网站| 国产精品亚洲美女久久久| 一a级毛片在线观看| 最近最新免费中文字幕在线| 国产精品野战在线观看| 男女下面进入的视频免费午夜| 成人午夜高清在线视频| 怎么达到女性高潮| 精华霜和精华液先用哪个| 啪啪无遮挡十八禁网站| 午夜久久久久精精品| 日本黄大片高清| 19禁男女啪啪无遮挡网站| xxxwww97欧美| 国产亚洲精品久久久久久毛片| 亚洲精品美女久久av网站| 精品一区二区三区四区五区乱码| 国语自产精品视频在线第100页| 久久精品国产清高在天天线| 精品一区二区三区四区五区乱码| 在线观看66精品国产| 深夜精品福利| 亚洲欧美精品综合久久99| 久热爱精品视频在线9| 天堂动漫精品| 99久久综合精品五月天人人| 99久久久亚洲精品蜜臀av| 最近视频中文字幕2019在线8| 欧美一级毛片孕妇| 两个人的视频大全免费| 两个人视频免费观看高清| 伦理电影免费视频| 欧美又色又爽又黄视频| 免费看十八禁软件| 精品高清国产在线一区| 色综合亚洲欧美另类图片| 麻豆国产av国片精品| 欧美高清成人免费视频www| 丁香欧美五月| 欧美人与性动交α欧美精品济南到| 亚洲国产日韩欧美精品在线观看 | 午夜福利18| 亚洲乱码一区二区免费版| 日韩av在线大香蕉| 精品久久久久久久人妻蜜臀av| 亚洲精品一区av在线观看| 亚洲九九香蕉| 国产在线精品亚洲第一网站| 97碰自拍视频| 亚洲av熟女| 国产在线精品亚洲第一网站| 中文字幕人成人乱码亚洲影| 午夜精品在线福利| 成人高潮视频无遮挡免费网站| 国产精品一及| 日本一区二区免费在线视频| av免费在线观看网站| 深夜精品福利| 欧美成人免费av一区二区三区| 国产亚洲精品久久久久久毛片| 久久久久久久午夜电影| 熟女少妇亚洲综合色aaa.| 麻豆成人av在线观看| 久久久久久国产a免费观看| 亚洲欧美精品综合久久99| 动漫黄色视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲一码二码三码区别大吗| 成人精品一区二区免费| 色综合亚洲欧美另类图片| 狂野欧美白嫩少妇大欣赏| 亚洲一区中文字幕在线| 久久精品综合一区二区三区| 久久久久精品国产欧美久久久| 国产精品综合久久久久久久免费| 亚洲精品久久国产高清桃花| 欧美性长视频在线观看| 亚洲中文av在线| 精品日产1卡2卡| 亚洲精品在线美女| 成人午夜高清在线视频| 99久久99久久久精品蜜桃| 午夜a级毛片| 亚洲欧美日韩高清专用| 国产精品1区2区在线观看.| 国产午夜精品久久久久久| 99久久综合精品五月天人人| 神马国产精品三级电影在线观看 | 久久久久久人人人人人| 成人午夜高清在线视频| 韩国av一区二区三区四区| 91麻豆精品激情在线观看国产| x7x7x7水蜜桃| 国产一区二区三区视频了| 一夜夜www| 免费观看精品视频网站| 亚洲激情在线av| 欧美色欧美亚洲另类二区| 亚洲人成网站在线播放欧美日韩| 欧美成人性av电影在线观看| 亚洲中文日韩欧美视频| 一二三四社区在线视频社区8| 婷婷丁香在线五月| 人人妻人人澡欧美一区二区| 成人av在线播放网站| 欧美激情久久久久久爽电影| 久久久精品欧美日韩精品| 午夜久久久久精精品| 亚洲乱码一区二区免费版| 毛片女人毛片| 国产精品99久久99久久久不卡| 午夜福利免费观看在线| 免费在线观看日本一区| 精品久久久久久,| 男人的好看免费观看在线视频 | 亚洲精品久久国产高清桃花| 亚洲一区中文字幕在线| 我要搜黄色片| 欧美最黄视频在线播放免费| 国产欧美日韩一区二区三| 97超级碰碰碰精品色视频在线观看| 全区人妻精品视频| av超薄肉色丝袜交足视频| 在线国产一区二区在线| 中文字幕人妻丝袜一区二区| 国产人伦9x9x在线观看| 男女下面进入的视频免费午夜| 国产麻豆成人av免费视频| 小说图片视频综合网站| 99国产精品一区二区蜜桃av| 国产亚洲精品综合一区在线观看 | 搡老熟女国产l中国老女人| 亚洲激情在线av| 久久久国产成人精品二区| 12—13女人毛片做爰片一| 色精品久久人妻99蜜桃| 99国产精品一区二区三区| 国产精品九九99| √禁漫天堂资源中文www| 欧美+亚洲+日韩+国产| 亚洲精品久久国产高清桃花| 国产乱人伦免费视频| 国产精品爽爽va在线观看网站| 午夜久久久久精精品| 18禁黄网站禁片免费观看直播| 国产真实乱freesex| 国产乱人伦免费视频| 久久中文字幕人妻熟女| 麻豆一二三区av精品| 18禁国产床啪视频网站| 亚洲第一电影网av| 日韩 欧美 亚洲 中文字幕| 国产成人aa在线观看| 在线a可以看的网站| 在线观看美女被高潮喷水网站 | 国产91精品成人一区二区三区| 又爽又黄无遮挡网站| 国产精品综合久久久久久久免费| 精品欧美国产一区二区三| 国产精品99久久99久久久不卡| 精品国产乱码久久久久久男人| 国产精品电影一区二区三区| 欧美黄色片欧美黄色片| 制服人妻中文乱码| 国产蜜桃级精品一区二区三区| 国产单亲对白刺激| 亚洲真实伦在线观看| 国产av又大| 欧美精品亚洲一区二区| 欧美成人性av电影在线观看| 99国产精品一区二区三区| 90打野战视频偷拍视频| 99国产精品99久久久久| 国产三级黄色录像| 日日爽夜夜爽网站| 天堂√8在线中文| 在线播放国产精品三级| 亚洲片人在线观看| 老司机福利观看| 日本五十路高清| 久久午夜亚洲精品久久| av福利片在线观看| 亚洲国产精品成人综合色| 香蕉久久夜色| 久久中文字幕一级| 好看av亚洲va欧美ⅴa在| 国产一区二区激情短视频| 一本精品99久久精品77| 久久婷婷人人爽人人干人人爱| 国产区一区二久久| av中文乱码字幕在线| 麻豆国产av国片精品| 熟女电影av网| 久久久久九九精品影院| 麻豆成人av在线观看| 人成视频在线观看免费观看| 岛国在线观看网站| 大型av网站在线播放| 特级一级黄色大片| 午夜亚洲福利在线播放| 香蕉久久夜色| 免费观看人在逋| 中文字幕高清在线视频| 久久久国产精品麻豆| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产欧洲综合997久久,| 亚洲在线自拍视频| 亚洲精品美女久久久久99蜜臀| 级片在线观看| 亚洲,欧美精品.| 亚洲国产看品久久| 国产69精品久久久久777片 | 国产精品av久久久久免费| 天堂av国产一区二区熟女人妻 | 天堂影院成人在线观看| 18禁国产床啪视频网站| xxx96com| www.999成人在线观看| 欧美极品一区二区三区四区| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类 | 国产亚洲欧美在线一区二区| 亚洲精品av麻豆狂野| 亚洲精品美女久久久久99蜜臀| 日本一本二区三区精品| 波多野结衣高清作品| 黑人欧美特级aaaaaa片| 五月玫瑰六月丁香| 国产99久久九九免费精品| 国产精品久久久久久精品电影| 欧洲精品卡2卡3卡4卡5卡区| 久久精品aⅴ一区二区三区四区| 婷婷丁香在线五月| av欧美777| 成人三级做爰电影| 国产精品av视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看| 日本 av在线| 精品久久久久久成人av| 男女做爰动态图高潮gif福利片| 人妻久久中文字幕网| 一级a爱片免费观看的视频| av片东京热男人的天堂| 无遮挡黄片免费观看| 两人在一起打扑克的视频| av天堂在线播放| 狂野欧美白嫩少妇大欣赏| 日本免费一区二区三区高清不卡| 热99re8久久精品国产| www.自偷自拍.com| 精品一区二区三区四区五区乱码| 国产亚洲精品第一综合不卡| 女人被狂操c到高潮| 无遮挡黄片免费观看| 又爽又黄无遮挡网站| 天天躁狠狠躁夜夜躁狠狠躁| 成人三级做爰电影| 久久亚洲精品不卡| 国产高清videossex| 他把我摸到了高潮在线观看| 高清在线国产一区| 免费电影在线观看免费观看| 久久久久久久久免费视频了| 九色成人免费人妻av| 欧美大码av| av中文乱码字幕在线| 久99久视频精品免费| 久久久久久免费高清国产稀缺| 欧美日韩乱码在线| 午夜日韩欧美国产| 91av网站免费观看| 欧美不卡视频在线免费观看 | 18禁观看日本| aaaaa片日本免费| 午夜影院日韩av| 国产精品亚洲av一区麻豆| 国产男靠女视频免费网站| 久久这里只有精品19| 久久久久亚洲av毛片大全| 精品一区二区三区四区五区乱码| 久久国产乱子伦精品免费另类| 狂野欧美激情性xxxx| 搡老岳熟女国产| 欧美极品一区二区三区四区| 欧美午夜高清在线| 一级毛片女人18水好多| 色综合站精品国产| 91九色精品人成在线观看| 国产亚洲精品久久久久5区| 久久精品亚洲精品国产色婷小说| 99国产精品一区二区蜜桃av| 国产亚洲av高清不卡| 给我免费播放毛片高清在线观看| 一进一出抽搐gif免费好疼| bbb黄色大片| 男人舔奶头视频| 99久久国产精品久久久| 久久国产乱子伦精品免费另类| 特大巨黑吊av在线直播| 国产蜜桃级精品一区二区三区| 91大片在线观看| 亚洲精品美女久久久久99蜜臀| 窝窝影院91人妻| 国产精品爽爽va在线观看网站| 亚洲一区二区三区不卡视频| 欧美成人性av电影在线观看| 国产视频一区二区在线看| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 99热6这里只有精品| 亚洲成人精品中文字幕电影| 精品久久久久久成人av| 亚洲av美国av| 欧美一级a爱片免费观看看 | 欧美精品亚洲一区二区| 制服诱惑二区| 国产男靠女视频免费网站| 久久久久国产一级毛片高清牌| 国产主播在线观看一区二区| 亚洲人成网站高清观看| videosex国产| 久久久久国产精品人妻aⅴ院| 国产精品 国内视频| 日本黄色视频三级网站网址| 在线观看日韩欧美| 亚洲欧美日韩无卡精品| 精品午夜福利视频在线观看一区| 草草在线视频免费看| 精品第一国产精品| 国产精品九九99| 国产午夜精品论理片| cao死你这个sao货| 在线观看舔阴道视频| 舔av片在线| 99国产精品一区二区三区| 亚洲五月天丁香| 亚洲中文字幕一区二区三区有码在线看 | 国产黄a三级三级三级人| 人人妻人人看人人澡| 免费观看人在逋| 国语自产精品视频在线第100页| 黄片大片在线免费观看| 日韩欧美国产一区二区入口| 欧美日韩亚洲国产一区二区在线观看| 欧美黄色片欧美黄色片| 99久久99久久久精品蜜桃| 久久精品国产亚洲av高清一级| 又紧又爽又黄一区二区| 亚洲成人久久爱视频| 天天躁狠狠躁夜夜躁狠狠躁| 中文在线观看免费www的网站 | 在线十欧美十亚洲十日本专区| 久久天躁狠狠躁夜夜2o2o| 国产一级毛片七仙女欲春2| 女警被强在线播放| 成人三级黄色视频| 最新美女视频免费是黄的| 亚洲精品国产精品久久久不卡| 一本综合久久免费| 最近最新免费中文字幕在线| 一区二区三区高清视频在线| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 久久久久久久久中文| 成人午夜高清在线视频| 日本黄大片高清| 亚洲成a人片在线一区二区| 国产又色又爽无遮挡免费看| 男女那种视频在线观看| 听说在线观看完整版免费高清| 国产精品国产高清国产av| 日本一二三区视频观看| 12—13女人毛片做爰片一| 欧美丝袜亚洲另类 | 久久性视频一级片| 亚洲精品在线观看二区| 日本在线视频免费播放| 亚洲国产欧美人成| 婷婷亚洲欧美| 别揉我奶头~嗯~啊~动态视频| 最好的美女福利视频网| 国产高清videossex| av有码第一页| 国产精品一区二区三区四区免费观看 | 18禁国产床啪视频网站| 欧美高清成人免费视频www| 99久久无色码亚洲精品果冻| 麻豆久久精品国产亚洲av| 女人被狂操c到高潮| 1024香蕉在线观看| 亚洲中文日韩欧美视频| 欧美日韩福利视频一区二区| 亚洲国产精品999在线| 国内久久婷婷六月综合欲色啪| 天堂√8在线中文| 免费在线观看视频国产中文字幕亚洲| 久久久久九九精品影院| 最好的美女福利视频网| 国产aⅴ精品一区二区三区波| 欧美性猛交黑人性爽| 成年版毛片免费区| 国产精品美女特级片免费视频播放器 | 午夜免费观看网址| 亚洲中文日韩欧美视频| 国产精品久久视频播放| 神马国产精品三级电影在线观看 | 黄色视频不卡| 婷婷亚洲欧美| 免费在线观看黄色视频的| 97人妻精品一区二区三区麻豆| 国产精品一区二区免费欧美| 久久人妻福利社区极品人妻图片| 亚洲 欧美 日韩 在线 免费| 亚洲成a人片在线一区二区| 丰满的人妻完整版| 欧美高清成人免费视频www| 精品久久久久久久末码| 波多野结衣高清无吗| 黄色丝袜av网址大全| 国产精品日韩av在线免费观看| 免费无遮挡裸体视频| 毛片女人毛片| 国产精品久久久av美女十八| 日韩大尺度精品在线看网址| 最近最新中文字幕大全免费视频| 色哟哟哟哟哟哟| 午夜激情福利司机影院| 欧美午夜高清在线| 搡老熟女国产l中国老女人| a在线观看视频网站| 男女午夜视频在线观看| 国产在线精品亚洲第一网站| 日本a在线网址| 99久久久亚洲精品蜜臀av| 亚洲av电影不卡..在线观看| 国产三级黄色录像| av福利片在线观看| 99久久精品国产亚洲精品| 久久精品国产清高在天天线| 搡老熟女国产l中国老女人| 男人的好看免费观看在线视频 | 最近最新中文字幕大全电影3| 国产精品av久久久久免费| 一区二区三区激情视频| 欧美丝袜亚洲另类 | 操出白浆在线播放| 婷婷六月久久综合丁香| 午夜日韩欧美国产| 一本久久中文字幕| 丁香欧美五月| 日日摸夜夜添夜夜添小说| 亚洲黑人精品在线| 国产精品影院久久| 一本精品99久久精品77| 一个人免费在线观看电影 | 99在线人妻在线中文字幕| 成年版毛片免费区| 中文字幕精品亚洲无线码一区| 午夜免费观看网址| 一本大道久久a久久精品| 精品欧美国产一区二区三| 男女下面进入的视频免费午夜| 中文字幕av在线有码专区| 十八禁人妻一区二区| 久久人人精品亚洲av| 久久久精品大字幕| 免费人成视频x8x8入口观看| www.自偷自拍.com| 成人18禁高潮啪啪吃奶动态图| 一夜夜www| 999久久久国产精品视频| 精品久久久久久久久久久久久| 香蕉丝袜av| 我的老师免费观看完整版| 国产精品久久电影中文字幕| 19禁男女啪啪无遮挡网站| 变态另类成人亚洲欧美熟女| 成年免费大片在线观看| 亚洲18禁久久av| 免费无遮挡裸体视频| 国产视频内射| 哪里可以看免费的av片| 三级男女做爰猛烈吃奶摸视频| 欧美丝袜亚洲另类 | 黄色片一级片一级黄色片| 香蕉av资源在线| 不卡av一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线永久观看黄色视频| 香蕉久久夜色| 欧美日韩福利视频一区二区| 国产精品一及| 欧美激情久久久久久爽电影| 免费在线观看成人毛片| 一二三四社区在线视频社区8| 婷婷精品国产亚洲av| 女人被狂操c到高潮| av有码第一页| 一区二区三区激情视频| 美女扒开内裤让男人捅视频| 最近视频中文字幕2019在线8|