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

    PBX 炸藥損傷本構(gòu)模型及其工程運(yùn)用

    2022-03-17 07:28:00黃垂藝金朋剛
    含能材料 2022年3期
    關(guān)鍵詞:破片子程序本構(gòu)

    黃垂藝,時 巖,金朋剛,陳 凱

    (1. 南京理工大學(xué)機(jī)械工程學(xué)院,江蘇 南京 210094;2. 西安近代化學(xué)研究所,陜西 西安 710065)

    1 引言

    高聚物黏結(jié)炸藥((PBX)是一種以高能炸藥為主體,與黏結(jié)劑、降感劑和增塑劑等輔助成分混合制成的高能鈍感的混合炸藥,其使用場景中受到多種載荷的作用,可能產(chǎn)生各種損傷,影響其力學(xué)性能和可靠性。PBX 材料是武器裝備系統(tǒng)中最薄弱的環(huán)節(jié)。因此,對PBX 材料的力學(xué)性能及損傷特性進(jìn)行研究,對武器裝備的設(shè)計、可靠運(yùn)用有著重要作用。

    近年,世界相關(guān)學(xué)者對PBX 材料的研究報道主要集中在準(zhǔn)靜態(tài)及動態(tài)加載下的力學(xué)性能的研究與分析。張子敏[1]、屈可朋[2]、孫文旭等[3]通過準(zhǔn)靜態(tài)實(shí)驗(yàn)和霍普金森壓桿(SHPB)試驗(yàn)獲得了PBX 材料在一定應(yīng)變率范圍內(nèi)的力學(xué)性能,并用相關(guān)本構(gòu)模型進(jìn)行描述。目前國內(nèi)外對于PBX 材料研究的重難點(diǎn)在于其受到?jīng)_擊載荷后的損傷特性及其演化規(guī)律,其中的困難在于PBX 材料內(nèi)部的損傷難以捕捉與表征,以及由于工藝、設(shè)備、批次不同導(dǎo)致的初始損傷的無規(guī)律分布;此外,由于PBX 材料的損傷特性及其演化規(guī)律與應(yīng)變率、溫度有關(guān),缺乏考慮全面的變量來表示其細(xì)觀破壞與其宏觀力學(xué)性能響應(yīng)之間的關(guān)系。因此,目前對PBX 材料損傷的相關(guān)研究主要是采用專業(yè)儀器捕捉其破壞圖像,并進(jìn)行分析總結(jié),增進(jìn)對其破壞過程的認(rèn)識[4]。陳鵬萬等[5]對低速氣炮加載后的PBX 炸藥的沖擊損傷形貌進(jìn)行了顯微觀測;劉本德等[6]用X 射線限位層析成像、CT 圖像序列結(jié)合數(shù)字圖像處理算法對SHPB 沖擊試驗(yàn)后的樣本損傷裂紋進(jìn)行提取、三維重構(gòu)與分析。

    目前對PBX 材料的研究與評估中,相比試驗(yàn),數(shù)值仿真能展現(xiàn)更多細(xì)節(jié)信息,是一種效果更優(yōu)的研究手段。同時,商業(yè)有限元軟件給PBX 材料的研究提供巨大的便利性,如材料本構(gòu)的二次開發(fā)。雖然PBX 材料的力學(xué)行為和損傷特性沒有相應(yīng)的模型可對其進(jìn)行準(zhǔn)確描述[7],但由于宏觀上PBX 材料的力學(xué)響應(yīng)和橡膠材料、高分子材料及推進(jìn)劑材料的力學(xué)響應(yīng)相似,且對上述材料的研究更為深入和完整,可從上述材料的相關(guān)研究中參考有意義的部分對PBX 材料的數(shù)值仿真開展進(jìn)一步研究。

    本工作依托PBX 材料的動態(tài)力學(xué)性能試驗(yàn)結(jié)果,在有限元軟件ABAQUS 中的材料子程序VUMAT 功能下,以相關(guān)學(xué)者運(yùn)用較多的含損傷變量的Z-W-T 本構(gòu)模型為框架對PBX 炸藥進(jìn)行研究,解決了前人研究存在的擬合參數(shù)非多應(yīng)變率共用的問題;并對PBX 材料的動態(tài)力學(xué)性能及其損傷與演化過程的損傷變量進(jìn)行了定義,得到了可以進(jìn)行可視化損傷云圖表示且具備應(yīng)變率相關(guān)性的材料模型,并對PBX 材料進(jìn)行帶殼裝藥破片撞擊仿真與分析,該材料子程序可為后續(xù)PBX 材料的力學(xué)性能數(shù)值實(shí)現(xiàn)、損傷表征數(shù)值實(shí)現(xiàn)及安全性數(shù)值評估、裝藥結(jié)構(gòu)改進(jìn)的研究提供基礎(chǔ)。

    2 含損傷變量的Z-W-T 本構(gòu)模型

    Z-W-T 本構(gòu)模型是由朱兆祥、王禮立和唐志平根據(jù)有限粘彈性本構(gòu)方程和Green-Rivlin 本構(gòu)理論提出的適用于描述熱塑性或熱固性高分子材料非線性粘彈性行為的本構(gòu)模型[8],已被廣泛運(yùn)用于推進(jìn)劑材料等粘彈性材料的力學(xué)性能表征。

    Z-W-T 非線性粘彈性本構(gòu)模型的積分形式為:

    式中,ε為應(yīng)變;fe(ε)為材料的非線性彈性響應(yīng)項(xiàng),E0、α和β是對應(yīng)的非線性彈性常數(shù),E0單位為Pa;式(1)第二項(xiàng)描述低應(yīng)變率下的粘彈性響應(yīng),E1和θ1分別是其對應(yīng)的低頻Maxwell 單元的彈性常數(shù)和松弛時間;第三項(xiàng)描述高應(yīng)變率下的粘彈性響應(yīng),E2和θ2分別是所對應(yīng)的高頻Maxwell 單元的彈性常數(shù)和松弛時間,彈性常數(shù)的單位為Pa,松弛時間的單位為s。低頻和高頻Maxwell 單元即為分別負(fù)責(zé)低應(yīng)變率和高應(yīng)變率的彈簧—粘壺組合模型。

    早期在對Z-W-T 本構(gòu)模型進(jìn)行研究時,所研究的材料在應(yīng)變達(dá)到一定值后,會出現(xiàn)微裂紋損傷。隨著應(yīng)變的繼續(xù)增大,微裂紋損傷的數(shù)量增多,導(dǎo)致材料的抗變形能力削弱,引起材料的弱化和軟化。因此對本構(gòu)模型引入Kachanov 損傷變量[8],得到:

    式中,σ為無損傷應(yīng)力,Pa;σr為含損傷應(yīng)力,Pa;D為損傷變量,無量綱。Z-W-T 本構(gòu)的模型的提出者與早期研究者對D的定義如下:

    式中,εth為應(yīng)變損傷閾值,當(dāng)應(yīng)變大于該值時出現(xiàn)微裂紋;D0為初始損傷因子;b為損傷應(yīng)變指數(shù)因子;δ為率相關(guān)的指數(shù)因子;b和δ是材料參數(shù)。應(yīng)變損傷閾值εth在以往的研究中一般定為固定值或0,但根據(jù)試驗(yàn)獲得的材料本構(gòu)曲線,該值應(yīng)該與應(yīng)變率相關(guān),根據(jù)參考文獻(xiàn)[9],應(yīng)變損傷閾值與應(yīng)變率的關(guān)系應(yīng)該滿足:

    式中,A和B為一般參數(shù);ε˙為應(yīng)變率。

    近年來,研究PBX 炸藥力學(xué)性能的學(xué)者在研究過程中將Z-W-T 非線性粘彈性本構(gòu)模型運(yùn)用在描述PBX炸藥的力學(xué)性能中獲得了比較好的結(jié)果[1-3]。

    3 動態(tài)力學(xué)性能試驗(yàn)

    3.1 試驗(yàn)裝置

    采用SHPB 作為PBX 炸藥的動態(tài)壓縮性能測試裝置,如圖1 所示。PBX 炸藥強(qiáng)度較低,同時由于內(nèi)部高聚物黏結(jié)劑的黏彈性使其具有吸能緩沖效果,導(dǎo)致應(yīng)力加載過程較為緩慢,不同應(yīng)變率下動態(tài)屈服強(qiáng)度相差不大,失效應(yīng)變反而相差較大。要有效地獲取其具有應(yīng)變率相關(guān)性、非線性黏彈性及應(yīng)變軟化效應(yīng)的力學(xué)性能,對試件應(yīng)力均勻性、常應(yīng)變率加載及信號測量準(zhǔn)確性提出了較高的要求[10]。因此,從整形器與試件的設(shè)計方面加以考慮,設(shè)計了不同材料與尺寸的整形器與試件,通過試驗(yàn)篩選合適的組合以滿足測試要求。因此,采用波阻抗較低的鋁合金作為桿的材料,并采用半導(dǎo)體應(yīng)變片測量透射信號及將橡膠作為整形器材料。所設(shè)計的整形器、試件、子彈及壓桿的材料與尺寸見表1。

    圖1 霍普金森壓桿(SHPB)動態(tài)壓縮試驗(yàn)裝置Fig.1 Hopkinson dynamic compression test device

    表1 SHPB 主要設(shè)計參數(shù)Table 1 Main design parameters of SHPB

    3.2 試驗(yàn)結(jié)果與分析

    利用SHPB 測量得到入射波、反射波及透射波信號后,通過專門的數(shù)據(jù)處理軟件計算應(yīng)力-應(yīng)變曲線及應(yīng)變率-時間曲線。應(yīng)力均勻性與常應(yīng)變率加載是SHPB 準(zhǔn)確獲取動態(tài)力學(xué)性能的前提條件,圖2 為試驗(yàn)測量得到的入射波、反射波及透射波信號,從圖2 可以看出:(1)根據(jù)應(yīng)力均勻性判別條件,試件已達(dá)到應(yīng)力均勻;(2)半導(dǎo)體應(yīng)變片準(zhǔn)確地捕獲到了強(qiáng)度較弱的透射信號,為應(yīng)力、應(yīng)變及應(yīng)變率等力學(xué)性能參數(shù)的正確計算提供了保證;(3)反射波實(shí)現(xiàn)了平臺波,標(biāo)志著試件獲得了常應(yīng)變率加載,在圖3 所示的應(yīng)變率-時間變化曲線中得到了印證。

    圖2 PBX 炸藥典型入射波、反射波、透射波信號Fig.2 Signal of incident wave,reflected wave and transmitted wave

    圖3 PBX 炸藥典型應(yīng)變率-時間曲線Fig.3 Strain rate vs. time curve

    圖4 所示為PBX 炸藥在1300~6100 s-1應(yīng)變率范圍內(nèi)的應(yīng)力-應(yīng)變變化曲線,從圖4 可以發(fā)現(xiàn):(1)隨著應(yīng)變率的增高,動態(tài)屈服強(qiáng)度從0.96 MPa(1300 s-1)增加到4.52 MPa(6100 s-1),反映了一定的應(yīng)變率效應(yīng);(2)曲線開始階段,表現(xiàn)出較強(qiáng)的非線性黏彈性,且應(yīng)力上升較為緩慢。到達(dá)屈服后表現(xiàn)為應(yīng)變軟化效應(yīng),這是因?yàn)椴牧蟽?nèi)部發(fā)生損傷及演化,承載能力降低。由于PBX 的非線性黏彈性與應(yīng)變軟化效應(yīng),導(dǎo)致其失效應(yīng)變在應(yīng)變率范圍內(nèi)變化較大,從1300 s-1沖擊作用下的0.18 增加到6100 s-1沖擊作用下的0.62,說明PBX 炸藥具有較強(qiáng)的韌性與抗沖擊破壞能力。

    圖4 PBX 炸藥不同應(yīng)變率下的應(yīng)力-應(yīng)變曲線Fig.4 Stress strain curves of cast PBX under different strain rates

    3.3 本構(gòu)模型參數(shù)獲取

    得到PBX 的力學(xué)性能后需要擬合合適的本構(gòu)參數(shù)對其進(jìn)行表征。Z-W-T 本構(gòu)模型的相關(guān)研究表明,低頻Maxwell 單元松弛時間θ1通常是10~102量級,高頻Maxwell 單元的松弛時間θ2通常是10-4~10-6量級。本研究所進(jìn)行的動態(tài)力學(xué)性能試驗(yàn)僅在10-4s 內(nèi)完成,代表低頻Maxwell 單元的部分沒有足夠的時間松弛,根據(jù)參考文獻(xiàn)[2]中的處理方法,忽略代表低頻Maxwell 部分的積分項(xiàng),因此,式(1)可化為式(6):

    根據(jù)參考文獻(xiàn)[11]中提及的內(nèi)容及仿真軟件的算例驗(yàn)證方法,作為有限元分析基礎(chǔ)的材料應(yīng)力-應(yīng)變關(guān)系需要是全曲線準(zhǔn)確的。試驗(yàn)獲得的工程應(yīng)力應(yīng)變不能準(zhǔn)確地反映材料的應(yīng)力應(yīng)變情況,需按照式(7)和式(8)將工程應(yīng)力應(yīng)變轉(zhuǎn)換為真實(shí)應(yīng)力應(yīng)變。圖5為轉(zhuǎn)換后的真實(shí)應(yīng)力-應(yīng)變曲線。

    圖5 PBX 炸藥在不同應(yīng)變率下的真實(shí)應(yīng)力-應(yīng)變曲線Fig.5 Real stress-strain curves of PBX explosives at different strain rates

    首先,在多個應(yīng)變率達(dá)到損傷之前的真實(shí)本構(gòu)數(shù)據(jù)基礎(chǔ)上,運(yùn)用Origin Pro 軟件與Levenberg-Marquardt優(yōu)化算法對式(6)中的5 個參數(shù),E0、α、β、E2和θ2,進(jìn)行多參數(shù)擬合,其中上述5 個參數(shù)為共用參數(shù),而應(yīng)變率作為獨(dú)立變量不參與共用,擬合得到的參數(shù)見表2。同時,本研究采用的擬合方法解決了參考文獻(xiàn)[12]中不同應(yīng)變率下ZWT 非線性粘彈性本構(gòu)模型的參數(shù)值難以恒定,無法使用統(tǒng)一本構(gòu)參數(shù)描述材料力學(xué)行為的問題。擬合得到的結(jié)果為多應(yīng)變率共用,可以描述PBX 炸藥在不同應(yīng)變率情況下的力學(xué)行為。

    表2 無損傷本構(gòu)參數(shù)Table 2 Undamaged constitutive parameters

    圖6 為運(yùn)用擬合參數(shù)的Z-W-T 本構(gòu)模型曲線與試驗(yàn)得到的PBX 炸藥本構(gòu)曲線。經(jīng)過誤差分析,Z-W-T本構(gòu)方程對PBX 炸藥力學(xué)性能擬合的相關(guān)指數(shù)可達(dá)0.95 以上。擬合參數(shù)可較好的表示未達(dá)到損傷前的PBX 炸藥力學(xué)性能。

    圖6 擬合Z-W-T 本構(gòu)曲線與試驗(yàn)本構(gòu)曲線對比Fig.6 Comparison between fitting Z-W-T constitutive curves and experimental constitutive curves

    其次,為了擬合損傷變量中的參數(shù),將應(yīng)變率在3300 s-1和6100 s-1條件下達(dá)到損傷的應(yīng)變范圍內(nèi)的損傷應(yīng)力σr與擬合得到的不含損傷的本構(gòu)方程在相同應(yīng)變范圍內(nèi)的無損應(yīng)力σ相除,消除擬合方程中的自變量,得到式(9),并對其進(jìn)行參數(shù)擬合:

    D0初始損傷因子,根據(jù)所研究材料的真實(shí)條件測量評定得到,但目前對于PBX 炸藥材料,該值的測量與評定方法不成熟,因此本研究通過限定施加合理取值范圍進(jìn)行擬合,得到損傷變量參數(shù)D0=0.2,δ=1.28031,b=1.01783。

    擬合結(jié)果與實(shí)際結(jié)果的對比如圖7 所示。由圖7可見,擬合結(jié)果與實(shí)際結(jié)果較為吻合。

    圖7 本構(gòu)損傷部分的擬合結(jié)果與試驗(yàn)結(jié)果對比Fig.7 The comparison between the fitting results of constitutive damage and the experimental results

    最后,對如式(5)所示的Z-W-T 本構(gòu)模型的損傷閾值與應(yīng)變率相關(guān)關(guān)系,從試驗(yàn)得到的詳細(xì)數(shù)據(jù)中,找到對應(yīng)應(yīng)變率條件下本構(gòu)數(shù)據(jù)中出現(xiàn)損傷的數(shù)據(jù)點(diǎn),即隨應(yīng)變增大對應(yīng)應(yīng)力開始減小的數(shù)據(jù)點(diǎn),將該數(shù)據(jù)點(diǎn)的應(yīng)變值作為損傷閾值。1300 s-1應(yīng)變率對應(yīng)的損傷閾值為0.18637;3300 s-1應(yīng)變率對應(yīng)的損傷閾值為0.29994;6100 s-1應(yīng)變率對應(yīng)的損傷閾值為0.62453;用以上三組數(shù)據(jù)對式(5)中的參數(shù)進(jìn)行擬合,獲得表示應(yīng)變率與損傷閾值的關(guān)系,A=-0.12055,B=9.52536×10-5。

    綜上所述,對采用的Z-W-T 本構(gòu)模型,所有的參數(shù)見表3。

    表3 帶損傷Z-W-T 本構(gòu)模型參數(shù)Table 3 Parameters of Z-W-T constitutive model with damage

    4 數(shù)值仿真

    4.1 含損傷變量的Z-W-T 本構(gòu)材料子程序

    獲得合適的本構(gòu)參數(shù)對PBX 材料的力學(xué)性能進(jìn)行表征后,對其材料子程序進(jìn)行編寫。在大型商業(yè)有限元軟件ABAQUS 中,為了滿足多樣化的本構(gòu)模型需求,留有材料子程序的接口供使用者使用。本研究在高應(yīng)變率條件下的PBX 炸藥材料力學(xué)性能更符合非線性力學(xué)、瞬態(tài)響應(yīng)分析的范疇與ABAQUS/Explicit模塊的材料子程序VUMAT 的使用條件更為契合。根據(jù)VUMAT 材料子程序的開發(fā)邏輯,將帶損傷變量的Z-W-T 本構(gòu)模型從一維形式轉(zhuǎn)變?yōu)槿S增量形式是開發(fā)過程的重點(diǎn)和難點(diǎn)。

    根據(jù)彈塑性力學(xué),運(yùn)用第二Poila-Kirchhoff 應(yīng)力張量和Green 應(yīng)變張量將式(6)變換為三維張量形式[13]:

    式中,μ為泊松比。根據(jù)參考文獻(xiàn)[14]和[15]對PBX炸藥材料泊松比的描述可知,其泊松比與多種因素有關(guān),不同情況下其值不同,本研究根據(jù)參考文獻(xiàn)取值0.19。

    將高頻Maxwell 粘壺部分用Pij表示,當(dāng)Δt足夠小時,將t和t+ Δt時刻的Pij相減,可得到其在時間增量Δt內(nèi)的增量式(12):

    另外,對代表非線性彈性的部分直接進(jìn)行差分,得到其在時間增量的增量形式:

    結(jié)合式(12)和式(13),并且其中代表應(yīng)力與應(yīng)變的符號換為常用的符號,可得到不帶損傷的Z-W-T 本構(gòu)模型的三維張量增量形式:

    根據(jù)擬合出來的關(guān)系計算損傷閾值和損傷變量時,需要將三維張量形式的結(jié)果轉(zhuǎn)換為等效形式進(jìn)行計算。同時,本研究對損傷變量計算是假設(shè)其演化過程是各向同性的。等效應(yīng)變的計算形式為:

    根據(jù)式(5)計算應(yīng)變閾值,當(dāng)?shù)刃?yīng)變大于當(dāng)前時間增量應(yīng)變率對應(yīng)的應(yīng)變閾值時,運(yùn)用式(16)計算損傷變量D的增量形式:

    則t時刻的損傷變量為:

    因此,損傷后的應(yīng)力可按照式(18)計算:

    嚴(yán)格按照ABAQUS/VUMAT 材料子程序的接口,按照以上邏輯完成材料子程序的編寫。

    4.2 材料子程序的驗(yàn)證

    完成材料子程序的編寫后,根據(jù)參考文獻(xiàn)[16]中對二次開發(fā)材料本構(gòu)模型的驗(yàn)證方法,驗(yàn)證其正確性。建立了圖8 所示的模型,該模型是一個1 mm×1 mm×1 mm 的立方體,并且僅劃分為一個單元,設(shè)置為C3D8R 單元類型,底面施加U2=UR1=UR3=0 的約束,即約束了y 方向的平移自由度、繞x 軸和z 軸的旋轉(zhuǎn)自由度。為了便于控制應(yīng)變量頂面施加指向底面的固定位移載荷,相當(dāng)于對該單元進(jìn)行了壓縮受力仿真??色@得該單元的應(yīng)力-應(yīng)變曲線,并將其與實(shí)驗(yàn)結(jié)果進(jìn)行對比。

    圖8 驗(yàn)證子程序的計算模型Fig.8 Calculation model of verification subroutine

    對模型進(jìn)行的壓縮工況仿真,并將其力學(xué)響應(yīng)與試驗(yàn)結(jié)果進(jìn)行對比。該仿真中通過定義不同大小的位移載荷控制應(yīng)變量,通過定義相應(yīng)的幅值時間控制應(yīng)變率,由此進(jìn)行不同應(yīng)變率下的壓縮工況仿真。仿真結(jié)果與試驗(yàn)結(jié)果的對比如圖9 所示。從圖9 可以看出,Z-W-T 本構(gòu)材料子程序在仿真程序中對PBX 炸藥力學(xué)性能的再現(xiàn)整體較為可靠,在低應(yīng)變率與高應(yīng)變率的損傷部分的誤差會相對較大,總體結(jié)果良好。

    圖9 驗(yàn)證子程序的仿真結(jié)果與試驗(yàn)結(jié)果對比Fig.9 The comparison between the simulation results of the verification subroutine and the experimental results

    4.3 帶殼裝藥的破片撞擊仿真

    通過試驗(yàn)與材料子程序開發(fā)將PBX 炸藥材料的含損傷力學(xué)性能嵌入到商業(yè)有限元軟件后,對PBX 炸藥材料進(jìn)行典型沖擊工況的仿真,分析其損傷演化。武器裝備戰(zhàn)斗部受到破片撞擊是其在勤務(wù)過程中常見的意外情況,也是國內(nèi)外學(xué)者研究的典型工況。目前大多數(shù)研究都集中在破片撞擊作用下發(fā)生的燃燒、爆燃及爆轟等反應(yīng)。而破片撞擊導(dǎo)致的過度形變也會導(dǎo)致裝備戰(zhàn)斗部中的PBX 炸藥材料產(chǎn)生微孔洞和微裂紋之類的損傷,對炸藥的宏觀力學(xué)性能、爆轟及安全性能產(chǎn)生影響。因此,利用本課題開發(fā)的材料子程序,進(jìn)行帶殼裝藥的破片撞擊仿真,對損傷變量場進(jìn)行分析。

    實(shí)際情況中,破片多為不規(guī)則形狀,由于成本、工作量的原因,根據(jù)參考文獻(xiàn)[17]對物理模型進(jìn)行簡化,建立如圖10 所示的物理模型。破片為典型的平頂圓柱,尺寸為Φ12 mm×20 mm,破片與殼體間有1 mm間隙;PBX 炸藥尺寸為Φ70 mm×70 mm;戰(zhàn)斗部殼體厚度采用飛魚戰(zhàn)斗部的殼體厚度為20 mm;為減少計算成本,進(jìn)行四分之一建模。

    圖10 帶殼裝藥破片撞擊仿真物理模型Fig.10 Simulation physical model of shell charge fragment impact

    破片采用鋼材料,戰(zhàn)斗部殼體采用鋁材料,用Johnson-cook 材料模型定義,破片和戰(zhàn)斗部殼體的材料參數(shù)如表4 所示。

    表4 鋼和鋁的Johnson-cook 材料模型參數(shù)[18Table4 Johnson-cook material model parameters of steel and aluminum

    劃分完網(wǎng)格的模型有55萬個單元,根據(jù)參考文獻(xiàn)[19]中的破片速度梯度,將破片分別賦予600,800,1000 m·s-1的速度,對這三個工況進(jìn)行仿真。仿真結(jié)果如圖11 所示。

    從圖11 可以看出,在破片撞擊后,在與破片速度方向垂直的兩個平面,由于外殼的形變與PBX 炸藥的形變,這兩個部位出現(xiàn)了較大的間隙。在沖擊起爆的情況下,越大的裝藥間隙,越不容易形成爆轟反應(yīng),會殘留越多的藥量[20],同時也會在一定程度內(nèi)削弱爆轟波的強(qiáng)度[21]。從圖13 戰(zhàn)斗部裝藥損傷變量云圖可以看出,PBX 炸藥在直接遭受破片撞擊的表面會出現(xiàn)較大范圍的損傷,在直接撞擊面的遠(yuǎn)端面的中心和邊緣,也出現(xiàn)了損傷。對比圖12 和圖13,可以發(fā)現(xiàn)碰撞沖擊產(chǎn)生的應(yīng)力波經(jīng)過傳遞及作用后,對比損傷變量值較高位置的應(yīng)力與周圍損傷變量值較低位置的應(yīng)力,損傷值較高位置的應(yīng)力值較小,說明損傷影響了材料承載能力。同時,600,800,1000 m·s-1工況的仿真輸出結(jié)果表明:破片分別在80,84,90 μs 失去沖擊力,發(fā)生回彈。結(jié)合圖14 撞擊側(cè)單元D-t曲線中可以發(fā)現(xiàn),直接遭受撞擊部位并沒有在一開始即出現(xiàn)損傷,而直接遭受撞擊部位的邊緣卻在一開始出現(xiàn)了損傷。分析原因可能是直接撞擊部位的應(yīng)變率較大,損傷應(yīng)變閾值較高在撞擊瞬時并未達(dá)到其值;而撞擊的邊緣部位有一定大小的應(yīng)變且應(yīng)變率不大,因此在撞擊后迅速出現(xiàn)了損傷。隨著應(yīng)力波在PBX 炸藥材料中的作用,直接撞擊部位的應(yīng)變率降低后,應(yīng)變并沒有變小,則在直接撞擊面就出現(xiàn)了大面積的損傷。

    圖11 帶殼裝藥破片撞擊仿真整體應(yīng)力結(jié)果Fig.11 Simulation results of global stress of shell charging fragment impact

    圖12 帶殼裝藥破片撞擊仿真戰(zhàn)斗部應(yīng)力結(jié)果Fig.12 Simulation results of warhead stress of shell charging fragment impact

    圖13 戰(zhàn)斗部裝藥損傷變量云圖Fig.13 Damage variable cloud diagram of warhead charge

    圖14 各工況不同撞擊部位單元D-t 曲線Fig.14 D-t curves of different impact positions under different working conditions

    PBX 炸藥帶殼裝藥破片撞擊變形圖和PBX 炸藥損傷變量云圖分別從宏觀與微觀兩個角度,反映了帶殼戰(zhàn)斗部裝藥在破片撞擊情況下可能出現(xiàn)的對其爆炸性能及安全性能產(chǎn)生影響的區(qū)域,可為相關(guān)學(xué)者研究PBX 炸藥損傷特性與其裝藥結(jié)構(gòu)改進(jìn)提供依據(jù)。

    5 結(jié)論

    對PBX 炸藥材料從SHPB 試驗(yàn)獲取本構(gòu)、本構(gòu)模型表征、有限元軟件材料子程序二次開發(fā)和帶殼裝藥破片沖擊仿真進(jìn)行了帶損傷的PBX 炸藥材料模型的研究,得到以下結(jié)論:

    (1)通過SHPB 試驗(yàn)獲得了PBX 炸藥具有應(yīng)變率效應(yīng)且達(dá)到損傷范圍的本構(gòu)數(shù)據(jù),可真實(shí)反映PBX 炸藥的力學(xué)性能。

    (2)基于試驗(yàn)獲得的PBX 炸藥本構(gòu)數(shù)據(jù)和帶損傷的Z-W-T 本構(gòu)模型,采用分段、分步擬合的方法,擬合出了誤差極小的本構(gòu)模型參數(shù),可準(zhǔn)確表征PBX 炸藥的力學(xué)性能;并基于此完成了有限元軟件ABAQUS 材料子程序VUMAT 的開發(fā)工作,經(jīng)過仿真驗(yàn)證滿足實(shí)際工程的精度要求。

    (3)調(diào)用編寫的材料子程序進(jìn)行的帶殼裝藥破片撞擊仿真的仿真結(jié)果從宏觀變形與微觀損傷兩方面,對可能出現(xiàn)間隙、損傷微孔洞微裂紋的地方進(jìn)行了預(yù)測;同時仿真結(jié)果也表明,PBX 炸藥材料并不是在撞擊最劇烈的時刻出現(xiàn)損傷;損傷常出現(xiàn)在后續(xù)材料內(nèi)的應(yīng)力波傳遞過程中。相關(guān)仿真結(jié)果和材料子程序?qū)罄m(xù)研究者觀測其損傷部位與改進(jìn)其裝藥結(jié)構(gòu)提供了依據(jù)與工具。

    猜你喜歡
    破片子程序本構(gòu)
    爆轟驅(qū)動下軸向預(yù)制破片飛散特性研究
    離心SC柱混凝土本構(gòu)模型比較研究
    彈藥動態(tài)加載下破片測試方法
    中國測試(2018年10期)2018-11-17 01:58:40
    鋸齒形結(jié)構(gòu)面剪切流變及非線性本構(gòu)模型分析
    半預(yù)制破片戰(zhàn)斗部破片威力分析
    一種新型超固結(jié)土三維本構(gòu)模型
    淺談子程序在數(shù)控車編程中的應(yīng)用
    子程序在數(shù)控車加工槽中的應(yīng)用探索
    前向戰(zhàn)斗部破片散布均勻性研究*
    西門子840D系統(tǒng)JOG模式下PLC調(diào)用并執(zhí)行NC程序
    大码成人一级视频| 久久精品国产亚洲av高清一级| 亚洲 国产 在线| 国产成人精品久久二区二区免费| 亚洲免费av在线视频| 国产亚洲精品久久久久5区| 日本欧美视频一区| 侵犯人妻中文字幕一二三四区| 欧美精品av麻豆av| 青草久久国产| 成年人免费黄色播放视频| 成人国产一区最新在线观看| 国产在视频线精品| www.精华液| 中文欧美无线码| 久久精品人人爽人人爽视色| tube8黄色片| 国产亚洲av高清不卡| 久久青草综合色| 99热网站在线观看| 国产成人精品无人区| 欧美激情 高清一区二区三区| 国产无遮挡羞羞视频在线观看| 国产1区2区3区精品| 18禁裸乳无遮挡动漫免费视频| 亚洲av欧美aⅴ国产| 免费在线观看完整版高清| 亚洲专区中文字幕在线| 亚洲一区二区三区不卡视频| 日韩欧美在线二视频 | netflix在线观看网站| 精品欧美一区二区三区在线| 视频在线观看一区二区三区| 国产欧美日韩精品亚洲av| 国产亚洲精品第一综合不卡| 久热爱精品视频在线9| 在线观看免费高清a一片| 老司机福利观看| tocl精华| 性色av乱码一区二区三区2| 1024香蕉在线观看| 亚洲少妇的诱惑av| 国产精品一区二区精品视频观看| 啦啦啦视频在线资源免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品亚洲av国产电影网| 久久人妻福利社区极品人妻图片| 一进一出好大好爽视频| 日韩制服丝袜自拍偷拍| 9热在线视频观看99| 亚洲男人天堂网一区| 天天躁夜夜躁狠狠躁躁| 国产精品99久久99久久久不卡| 国产99久久九九免费精品| 中文字幕色久视频| 国产精华一区二区三区| 少妇裸体淫交视频免费看高清 | 久99久视频精品免费| 国产成人精品久久二区二区91| 国产欧美亚洲国产| 国产精品香港三级国产av潘金莲| 国产精品综合久久久久久久免费 | 国产免费av片在线观看野外av| 成人国产一区最新在线观看| 亚洲欧美日韩另类电影网站| 色播在线永久视频| 国产人伦9x9x在线观看| 国产成人免费观看mmmm| 后天国语完整版免费观看| 精品免费久久久久久久清纯 | 色尼玛亚洲综合影院| 国产成人一区二区三区免费视频网站| 一级,二级,三级黄色视频| 黑丝袜美女国产一区| a级片在线免费高清观看视频| 国产野战对白在线观看| 国产无遮挡羞羞视频在线观看| 欧美精品啪啪一区二区三区| 国精品久久久久久国模美| 夫妻午夜视频| 极品人妻少妇av视频| 桃红色精品国产亚洲av| 亚洲九九香蕉| 久久久久精品人妻al黑| 亚洲精品久久成人aⅴ小说| 99国产极品粉嫩在线观看| 天堂俺去俺来也www色官网| 女性被躁到高潮视频| 亚洲午夜理论影院| 亚洲在线自拍视频| 美女扒开内裤让男人捅视频| 美女 人体艺术 gogo| 高清毛片免费观看视频网站 | 在线观看一区二区三区激情| 老熟妇仑乱视频hdxx| 麻豆乱淫一区二区| 色婷婷久久久亚洲欧美| 一区二区三区国产精品乱码| 精品第一国产精品| www.精华液| 国产成人欧美在线观看 | 天堂俺去俺来也www色官网| 美女国产高潮福利片在线看| 欧美乱码精品一区二区三区| 精品国产乱子伦一区二区三区| 国内久久婷婷六月综合欲色啪| 在线av久久热| bbb黄色大片| tocl精华| 男人的好看免费观看在线视频 | 久久久精品区二区三区| 欧美激情久久久久久爽电影 | 免费观看人在逋| 久久久国产精品麻豆| 99精品欧美一区二区三区四区| 久久精品国产亚洲av香蕉五月 | 精品国产一区二区久久| 两性夫妻黄色片| 麻豆av在线久日| 欧美在线一区亚洲| 天堂中文最新版在线下载| 日韩熟女老妇一区二区性免费视频| 精品乱码久久久久久99久播| 国产蜜桃级精品一区二区三区 | 夜夜夜夜夜久久久久| 久99久视频精品免费| 男人操女人黄网站| 99riav亚洲国产免费| av网站免费在线观看视频| 狠狠婷婷综合久久久久久88av| 午夜福利免费观看在线| 国产午夜精品久久久久久| 一级a爱片免费观看的视频| 激情视频va一区二区三区| 午夜激情av网站| 国产精品乱码一区二三区的特点 | 精品国产一区二区三区久久久樱花| 一二三四社区在线视频社区8| 国产91精品成人一区二区三区| 80岁老熟妇乱子伦牲交| 日韩精品免费视频一区二区三区| 久久久国产欧美日韩av| 天天躁夜夜躁狠狠躁躁| 亚洲专区中文字幕在线| 女人被躁到高潮嗷嗷叫费观| 久久婷婷成人综合色麻豆| 久久久精品区二区三区| 国产亚洲av高清不卡| 欧美日韩精品网址| 欧美日韩瑟瑟在线播放| 欧美人与性动交α欧美软件| 9191精品国产免费久久| 看免费av毛片| 91老司机精品| 亚洲情色 制服丝袜| 久久久久久久久免费视频了| 久久人人爽av亚洲精品天堂| 校园春色视频在线观看| 亚洲综合色网址| 美女扒开内裤让男人捅视频| av网站在线播放免费| 国产激情久久老熟女| 夫妻午夜视频| 正在播放国产对白刺激| 国产高清videossex| 国产高清videossex| av天堂久久9| av天堂久久9| 精品久久久久久久久久免费视频 | 大香蕉久久成人网| 国产色视频综合| 久久午夜亚洲精品久久| 亚洲一卡2卡3卡4卡5卡精品中文| 18在线观看网站| 男人的好看免费观看在线视频| 禁无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清专用| 成人av在线播放网站| 熟女人妻精品中文字幕| 女警被强在线播放| 非洲黑人性xxxx精品又粗又长| 亚洲一区高清亚洲精品| 亚洲久久久久久中文字幕| 欧美大码av| 一级毛片女人18水好多| 日日夜夜操网爽| 亚洲午夜理论影院| 欧美乱色亚洲激情| 色老头精品视频在线观看| 国产乱人伦免费视频| 在线a可以看的网站| 精品一区二区三区视频在线观看免费| 精品国产亚洲在线| 人妻夜夜爽99麻豆av| 国产精品影院久久| 一级a爱片免费观看的视频| 午夜a级毛片| 精品乱码久久久久久99久播| АⅤ资源中文在线天堂| or卡值多少钱| 69人妻影院| 久久久久国内视频| 亚洲av免费高清在线观看| 一进一出抽搐动态| 国产单亲对白刺激| 国产久久久一区二区三区| 国产高清videossex| 免费观看精品视频网站| 人人妻人人看人人澡| 国产成+人综合+亚洲专区| 白带黄色成豆腐渣| 国产亚洲精品久久久久久毛片| 国产免费一级a男人的天堂| 此物有八面人人有两片| 亚洲真实伦在线观看| 少妇丰满av| 日本精品一区二区三区蜜桃| 久久精品影院6| 偷拍熟女少妇极品色| 久久6这里有精品| 俺也久久电影网| 亚洲人与动物交配视频| 国产成人欧美在线观看| 色尼玛亚洲综合影院| 看免费av毛片| 久久久久免费精品人妻一区二区| 岛国在线观看网站| 亚洲 欧美 日韩 在线 免费| 久久精品人妻少妇| 3wmmmm亚洲av在线观看| 亚洲av美国av| 亚洲狠狠婷婷综合久久图片| 亚洲aⅴ乱码一区二区在线播放| 韩国av一区二区三区四区| 成年版毛片免费区| 色哟哟哟哟哟哟| 最好的美女福利视频网| www.色视频.com| 成人国产综合亚洲| 在线观看66精品国产| 搡老熟女国产l中国老女人| 99久久久亚洲精品蜜臀av| 嫩草影视91久久| 男女下面进入的视频免费午夜| 深夜精品福利| 深爱激情五月婷婷| 午夜老司机福利剧场| 日韩有码中文字幕| 色综合欧美亚洲国产小说| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲av一区综合| 久久久久九九精品影院| 免费在线观看成人毛片| 欧美国产日韩亚洲一区| 一级毛片高清免费大全| 免费无遮挡裸体视频| 国内少妇人妻偷人精品xxx网站| 日韩成人在线观看一区二区三区| www.色视频.com| 国产精品爽爽va在线观看网站| 欧美黑人巨大hd| 日韩有码中文字幕| 欧美另类亚洲清纯唯美| 久久久精品大字幕| 亚洲无线观看免费| 国产精品久久久久久久电影 | 亚洲人与动物交配视频| 精品99又大又爽又粗少妇毛片 | 国产成人福利小说| 久久久久久久精品吃奶| 欧美色视频一区免费| 老司机福利观看| 日日摸夜夜添夜夜添小说| www日本在线高清视频| 婷婷精品国产亚洲av| 老熟妇仑乱视频hdxx| 热99在线观看视频| а√天堂www在线а√下载| 男人的好看免费观看在线视频| 美女高潮的动态| 亚洲天堂国产精品一区在线| 好看av亚洲va欧美ⅴa在| 在线天堂最新版资源| av中文乱码字幕在线| 18禁国产床啪视频网站| 久久香蕉国产精品| 精品久久久久久久末码| 国产精品美女特级片免费视频播放器| 蜜桃久久精品国产亚洲av| 色综合婷婷激情| 日本成人三级电影网站| 村上凉子中文字幕在线| 久久久久久久久久黄片| 男女之事视频高清在线观看| 青草久久国产| 精品久久久久久成人av| 天堂av国产一区二区熟女人妻| 日韩精品中文字幕看吧| 成人av一区二区三区在线看| av欧美777| 久久婷婷人人爽人人干人人爱| 亚洲熟妇熟女久久| 亚洲av电影在线进入| 禁无遮挡网站| 非洲黑人性xxxx精品又粗又长| 久久人妻av系列| 久久久久国产精品人妻aⅴ院| 欧美+亚洲+日韩+国产| 国产精品,欧美在线| 亚洲欧美日韩卡通动漫| 国产美女午夜福利| 男女那种视频在线观看| 天天一区二区日本电影三级| 嫩草影视91久久| 亚洲18禁久久av| 久久天躁狠狠躁夜夜2o2o| 成年女人毛片免费观看观看9| 日日摸夜夜添夜夜添小说| 我的老师免费观看完整版| www.色视频.com| 成人午夜高清在线视频| 中文字幕精品亚洲无线码一区| 制服人妻中文乱码| 757午夜福利合集在线观看| 久久欧美精品欧美久久欧美| 亚洲熟妇中文字幕五十中出| 国产精品久久视频播放| 色综合欧美亚洲国产小说| 在线国产一区二区在线| 国产99白浆流出| 精品国产美女av久久久久小说| 观看美女的网站| 日本撒尿小便嘘嘘汇集6| 窝窝影院91人妻| 亚洲精品在线美女| 黄色女人牲交| 亚洲欧美日韩卡通动漫| 国产精品亚洲av一区麻豆| a级一级毛片免费在线观看| 亚洲中文字幕一区二区三区有码在线看| 久久亚洲精品不卡| 免费大片18禁| 很黄的视频免费| 老司机福利观看| 精品久久久久久久久久免费视频| 在线天堂最新版资源| 色综合亚洲欧美另类图片| 亚洲五月天丁香| 亚洲人成网站在线播| 搞女人的毛片| 日韩精品青青久久久久久| 天堂动漫精品| 欧美成狂野欧美在线观看| 中文字幕人妻熟人妻熟丝袜美 | 中文字幕人成人乱码亚洲影| 国产私拍福利视频在线观看| 精品久久久久久久毛片微露脸| 热99re8久久精品国产| 99久久精品一区二区三区| 一进一出抽搐gif免费好疼| 变态另类丝袜制服| 在线观看日韩欧美| 伊人久久精品亚洲午夜| 色噜噜av男人的天堂激情| 亚洲国产精品成人综合色| 亚洲av电影不卡..在线观看| 欧美不卡视频在线免费观看| 国产不卡一卡二| 天天躁日日操中文字幕| 久久久久精品国产欧美久久久| 日本成人三级电影网站| 亚洲,欧美精品.| 日韩亚洲欧美综合| 国产欧美日韩精品一区二区| 久久精品人妻少妇| 黄色视频,在线免费观看| 国产精品亚洲av一区麻豆| 亚洲欧美日韩东京热| 日韩 欧美 亚洲 中文字幕| 久久亚洲真实| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品日产1卡2卡| 国产精品98久久久久久宅男小说| 丰满的人妻完整版| 久久久久九九精品影院| 高潮久久久久久久久久久不卡| 国产高清视频在线播放一区| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 9191精品国产免费久久| 美女大奶头视频| 欧美在线一区亚洲| 国产在视频线在精品| 国内精品美女久久久久久| 国产精品三级大全| 日本撒尿小便嘘嘘汇集6| 熟妇人妻久久中文字幕3abv| 最新中文字幕久久久久| 一本久久中文字幕| 美女被艹到高潮喷水动态| 亚洲真实伦在线观看| 男女之事视频高清在线观看| 国产精品 欧美亚洲| 日本成人三级电影网站| 精品乱码久久久久久99久播| 国产三级中文精品| 日韩欧美三级三区| 99久久精品国产亚洲精品| av在线蜜桃| 天堂√8在线中文| 女生性感内裤真人,穿戴方法视频| 嫁个100分男人电影在线观看| 色尼玛亚洲综合影院| 无限看片的www在线观看| 国产一区二区在线观看日韩 | 久久久久久久午夜电影| 亚洲狠狠婷婷综合久久图片| 亚洲最大成人手机在线| 国产精品永久免费网站| 国产精品久久久久久久电影 | 亚洲av熟女| 成人18禁在线播放| 欧美乱色亚洲激情| 欧美日本亚洲视频在线播放| а√天堂www在线а√下载| 在线免费观看不下载黄p国产 | 国产精品98久久久久久宅男小说| 日本与韩国留学比较| 久久久精品欧美日韩精品| 狂野欧美激情性xxxx| 一本综合久久免费| 又粗又爽又猛毛片免费看| 俄罗斯特黄特色一大片| 在线观看日韩欧美| 午夜精品一区二区三区免费看| 女同久久另类99精品国产91| 啦啦啦免费观看视频1| 国产伦人伦偷精品视频| 国产精品久久电影中文字幕| 久久精品国产亚洲av涩爱 | 热99在线观看视频| 国产伦人伦偷精品视频| 老汉色av国产亚洲站长工具| 精品免费久久久久久久清纯| 亚洲人与动物交配视频| 国产亚洲精品一区二区www| 日韩欧美精品v在线| 观看美女的网站| 精品人妻一区二区三区麻豆 | 中文在线观看免费www的网站| 日本五十路高清| 男人的好看免费观看在线视频| 在线国产一区二区在线| 欧美黑人欧美精品刺激| 老司机福利观看| 精华霜和精华液先用哪个| 国产精品美女特级片免费视频播放器| 无限看片的www在线观看| 亚洲精品一区av在线观看| 欧美zozozo另类| 欧美日韩亚洲国产一区二区在线观看| 熟妇人妻久久中文字幕3abv| 黄色女人牲交| 在线观看美女被高潮喷水网站 | 日韩精品中文字幕看吧| 午夜福利在线在线| 女人十人毛片免费观看3o分钟| 久久久久性生活片| 国产亚洲精品久久久久久毛片| 亚洲av成人av| 在线看三级毛片| 亚洲国产精品合色在线| 中文字幕人成人乱码亚洲影| 俺也久久电影网| 国产97色在线日韩免费| 18+在线观看网站| 丁香欧美五月| 琪琪午夜伦伦电影理论片6080| 成年人黄色毛片网站| 成人永久免费在线观看视频| 人人妻,人人澡人人爽秒播| 国产一区二区亚洲精品在线观看| 中出人妻视频一区二区| 婷婷亚洲欧美| 五月玫瑰六月丁香| 亚洲成av人片免费观看| 日日摸夜夜添夜夜添小说| www日本黄色视频网| 一个人免费在线观看电影| 成年免费大片在线观看| 欧美日韩精品网址| 免费电影在线观看免费观看| 国产成人aa在线观看| 午夜视频国产福利| 亚洲专区中文字幕在线| 69人妻影院| 精品人妻一区二区三区麻豆 | 亚洲熟妇熟女久久| 中出人妻视频一区二区| 日韩精品中文字幕看吧| 好男人电影高清在线观看| 丝袜美腿在线中文| 国产色爽女视频免费观看| 国产三级在线视频| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 悠悠久久av| 亚洲av电影在线进入| 亚洲第一欧美日韩一区二区三区| 国产一区二区三区视频了| 国产av不卡久久| 欧美日韩亚洲国产一区二区在线观看| 激情在线观看视频在线高清| 国产私拍福利视频在线观看| 国产真实乱freesex| 9191精品国产免费久久| 午夜影院日韩av| 亚洲最大成人手机在线| 午夜日韩欧美国产| 国产成人系列免费观看| tocl精华| 午夜免费男女啪啪视频观看 | 亚洲成a人片在线一区二区| 91九色精品人成在线观看| 色吧在线观看| 亚洲欧美日韩高清专用| 九色国产91popny在线| 国产精品综合久久久久久久免费| 国产一区二区激情短视频| 国产av一区在线观看免费| 国产极品精品免费视频能看的| 69人妻影院| svipshipincom国产片| 午夜两性在线视频| 超碰av人人做人人爽久久 | 日本撒尿小便嘘嘘汇集6| 熟妇人妻久久中文字幕3abv| 久久精品影院6| e午夜精品久久久久久久| 999久久久精品免费观看国产| 亚洲激情在线av| 一级黄片播放器| 不卡一级毛片| 亚洲av日韩精品久久久久久密| 免费看日本二区| 国产高潮美女av| 少妇的逼好多水| 又爽又黄无遮挡网站| av天堂中文字幕网| 国产精品爽爽va在线观看网站| 最近最新中文字幕大全电影3| 我的老师免费观看完整版| 少妇裸体淫交视频免费看高清| 淫妇啪啪啪对白视频| 国产男靠女视频免费网站| 麻豆国产av国片精品| 中文字幕av在线有码专区| 美女免费视频网站| 可以在线观看的亚洲视频| 国产精品乱码一区二三区的特点| 高清日韩中文字幕在线| 男人和女人高潮做爰伦理| 亚洲午夜理论影院| АⅤ资源中文在线天堂| 美女高潮的动态| 国产精品1区2区在线观看.| 午夜亚洲福利在线播放| 99久久成人亚洲精品观看| 中文字幕人成人乱码亚洲影| 给我免费播放毛片高清在线观看| 国产亚洲精品综合一区在线观看| 免费电影在线观看免费观看| 国内精品久久久久久久电影| 日韩精品中文字幕看吧| 有码 亚洲区| 婷婷丁香在线五月| 国产野战对白在线观看| 十八禁人妻一区二区| 白带黄色成豆腐渣| 国产一区二区激情短视频| 制服人妻中文乱码| 日本与韩国留学比较| 免费观看人在逋| 国产午夜精品论理片| 露出奶头的视频| 国产综合懂色| 高清日韩中文字幕在线| 国内久久婷婷六月综合欲色啪| 草草在线视频免费看| 成人av在线播放网站| 欧美+亚洲+日韩+国产| xxxwww97欧美| 亚洲成人精品中文字幕电影| 国产单亲对白刺激| 大型黄色视频在线免费观看| 老司机午夜十八禁免费视频| 国产三级中文精品| 麻豆国产av国片精品| 99久久久亚洲精品蜜臀av| 成人无遮挡网站| 在线天堂最新版资源| 给我免费播放毛片高清在线观看| 老司机深夜福利视频在线观看| www.999成人在线观看| 亚洲av电影不卡..在线观看| av天堂在线播放| 午夜精品久久久久久毛片777| 在线a可以看的网站| 一夜夜www| 国产精品亚洲美女久久久| 欧美zozozo另类| 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 亚洲色图av天堂|