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

    燃料拋撒的初始速率與結(jié)構(gòu)動(dòng)態(tài)響應(yīng)數(shù)值模擬

    2015-05-10 00:33:11陳明生白春華李建平
    含能材料 2015年4期
    關(guān)鍵詞:殼體間距燃料

    陳明生, 白春華, 李建平

    (北京理工大學(xué)爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100081)

    1 引 言

    中心拋撒藥驅(qū)動(dòng)固體、液體或固液混合物的拋撒技術(shù),目前廣泛應(yīng)用于消防[1]、爆破除塵[2]和云霧爆轟[3]等領(lǐng)域。從起爆到殼體膨脹、破裂為破片的過程[4],以及燃料拋撒、碎解、霧化的過程[5],是爆轟驅(qū)動(dòng)燃料拋撒[6]的兩個(gè)外在表現(xiàn)形式。

    試驗(yàn)研究[7-8]發(fā)現(xiàn)殼體結(jié)構(gòu)在沖擊下發(fā)生變形和斷裂,其破片斷面主要是由韌性剪切斷裂造成的。裂紋由殼體外側(cè)向內(nèi)側(cè)徑向發(fā)展,當(dāng)裝藥結(jié)構(gòu)內(nèi)部沖擊壓力大于殼體內(nèi)表面張緊環(huán)所能承受的極限應(yīng)力時(shí),材料由塑性變形到完全斷裂的情況就會(huì)發(fā)生[9]。也就是說在燃料從裂隙中泄壓之前,裝置結(jié)構(gòu)對(duì)燃料的拋撒運(yùn)動(dòng)狀態(tài)具有約束作用。因此適當(dāng)減小殼體材料的屈服強(qiáng)度和厚度可提高燃料近場(chǎng)的拋撒速度[10]。除了削弱殼體的約束,增加中心拋撒藥質(zhì)量也可以達(dá)到相同的提高效果,但在一定范圍內(nèi)燃料拋散的最終云霧半徑與拋撒藥質(zhì)量無關(guān)[11]。對(duì)于圓柱體裝藥結(jié)構(gòu),燃料在拋撒近場(chǎng)階段的加速過程可認(rèn)為是一維柱對(duì)稱運(yùn)動(dòng)[12]?;谥鶎?duì)稱關(guān)系,張奇[13]推導(dǎo)了燃料近場(chǎng)拋撒最大速度的計(jì)算方法; Gardner[14-15]和Glass[16]構(gòu)造了燃料拋撒的近場(chǎng)、遠(yuǎn)場(chǎng)以及近遠(yuǎn)場(chǎng)耦合模型,隨后Singh[17]對(duì)近場(chǎng)模型進(jìn)行了深化,研究了液滴的初始粒徑分布; 丁玨[18-19]建立了拋撒近場(chǎng)階段的初期一維單相流動(dòng)模型和后期兩相流數(shù)學(xué)模型,并給出了變質(zhì)量運(yùn)動(dòng)邊界的處理方法; 薛社生[20]利用有限差分程序模擬了遠(yuǎn)場(chǎng)的燃料拋散過程。

    以往研究的裝藥結(jié)構(gòu)多為圓柱體[21-22],而針對(duì)非圓截面殼體結(jié)構(gòu)的燃料拋撒研究較少。與圓柱體結(jié)構(gòu)不同,非圓截面殼體在外形上的不規(guī)則性,使其無法簡(jiǎn)化成單純的柱對(duì)稱模型。非圓截面裝置在各個(gè)方向上的燃料拋撒初始速率和殼體應(yīng)力分布是不同的。本研究針對(duì)一種扇形截面的裝藥結(jié)構(gòu),在中心拋撒藥驅(qū)動(dòng)下拋撒燃料的初始速率以及殼體的動(dòng)態(tài)響應(yīng)進(jìn)行了數(shù)值模擬,并在試驗(yàn)結(jié)果對(duì)比的基礎(chǔ)上,分析中心拋撒藥類型及其位置、不同殼體刻槽形式對(duì)燃料拋撒的影響規(guī)律。

    2 模擬參數(shù)

    2.1 結(jié)構(gòu)模型

    殼體的材料為鋁,厚度為3mm,其結(jié)構(gòu)橫截面為扇形,模型與尺寸如圖1所示。拋撒藥與弧線對(duì)應(yīng)的圓心距離為185 mm。拋撒藥的起爆位置為下端。模型中空氣區(qū)域尺寸為1 m×1 m×2 m??諝庥虿捎脽o反射邊界條件。計(jì)算中采用流-固耦合算法[23]。殼體、中心管與端蓋采用拉氏網(wǎng)格,拋撒藥、燃料和空氣采用歐拉網(wǎng)格。模型中的體單元形狀為六面體,殼單元形狀為四邊形。單元總數(shù)量約為250萬(wàn)個(gè),其中單元最大邊長(zhǎng)為15 mm,最小邊長(zhǎng)為1 mm。為描述方便,以幾何形心為中心對(duì)殼體橫截面進(jìn)行方向標(biāo)定,如圖2所示。外殼刻槽總數(shù)為18道,分別為0°方向2道,45°方向4道,180°方向8道,315°方向4道??滩坶g距為40 mm,刻槽深度為2 mm。

    圖1 結(jié)構(gòu)模型

    Fig.1 Structure model

    圖2 方向示意

    Fig.2 Direction view

    2.2 材料參數(shù)

    計(jì)算過程涉及的材料包括拋撒藥、空氣、燃料和殼體等物質(zhì)。材料模型與數(shù)值取自文獻(xiàn)[23-24]。拋撒藥采用JWL狀態(tài)方程,參數(shù)見表1。

    空氣材料選用LINEAR_POLYNOMIAL狀態(tài)方程。質(zhì)量密度ρ為1.29×10-3g·cm-3,單位體積內(nèi)能E0為2.5×105J·m-3,材料常數(shù)C4為0.4,C5為0.4。

    燃料用水代替,采用GRUNEISEN狀態(tài)方程。質(zhì)量密度ρ為1.0 g·cm-3,材料常數(shù)C為1.65,S1為1.92,S2為-0.096,γ為0.35。

    端蓋材料為45#鋼,選用PLASTIC_KINEMATIC模型。質(zhì)量密度ρ為7.83 g·cm-3,彈性模量E為210 GPa,泊松比ν為0.3,屈服應(yīng)力σs為0.355 GPa,屈服段模量Et為10 GPa,失效應(yīng)變量δ為1.2。

    中心管材料為25#鋼,選用PLASTIC_KINEMATIC模型。質(zhì)量密度ρ為7.83 g·cm-3,彈性模量E為210 GPa,泊松比ν為0.3,屈服應(yīng)力σs為0.275 GPa,屈服段模量Et為10 GPa,失效應(yīng)變量δ為0.8。

    殼體材料為鋁,選用PLASTIC_KINEMATIC模型。質(zhì)量密度ρ為2.7 g·cm-3,彈性模量E為6 GPa,泊松比ν為0.3,屈服應(yīng)力σs為0.262 GPa,屈服段模量Et為1 GPa,失效應(yīng)變量δ為0.6。

    為便于分析,針對(duì)分析對(duì)象做三點(diǎn)說明。第一,在燃料單元的選擇方面,由于形成云霧的燃料拋撒方向主要為徑向,且整體軸向中部橫截面的燃料拋撒效果最為顯著,因此在統(tǒng)計(jì)中僅對(duì)燃料中部橫截面周邊單元的水平方向拋撒速率進(jìn)行分析。同時(shí)為降低隨機(jī)誤差,取某一單元及其相鄰共5個(gè)單元的平均速率。第二,在拋撒速率的選取方面,根據(jù)燃料拋撒的近場(chǎng)特性并作為遠(yuǎn)場(chǎng)初始條件[26],選擇最大速率作為分析數(shù)據(jù)。其主要原因是由于在殼體破壞后燃料開始出現(xiàn)剝離破碎,計(jì)算模型中使用的水狀態(tài)方程顯然已經(jīng)不再適用。因此,需要將總計(jì)算時(shí)間控制在殼體破壞后的一定時(shí)間內(nèi)。通過前期試算發(fā)現(xiàn)取總計(jì)算時(shí)長(zhǎng)為2 ms時(shí),模擬結(jié)果與試驗(yàn)情況較為相符。第三,在殼體破壞視角方面,考慮既展現(xiàn)全面觀察殼體的破裂情況又節(jié)省篇幅,觀測(cè)視角采用俯視。

    2.3 模擬與試驗(yàn)驗(yàn)證

    為了與試驗(yàn)驗(yàn)證中的高速攝像系統(tǒng)拍攝視角保持一致,針對(duì)模擬結(jié)果選取0°方向0~2 ms時(shí)刻對(duì)應(yīng)的殼體響應(yīng)過程,如圖3所示。其中拋撒藥為梯恩梯(TNT)。

    為了對(duì)數(shù)值模擬的可靠性進(jìn)行驗(yàn)證,下面進(jìn)行試驗(yàn)對(duì)比。試驗(yàn)的主要條件為: 殼體下端板距離地面1.0 m。試驗(yàn)采用HG-100K高速攝像系統(tǒng),以2000 fps的拍攝速度在0°方向進(jìn)行圖像采集,拍攝點(diǎn)與拋撒中心距離為100 m。試驗(yàn)中采集的燃料拋撒過程如圖4所示。

    表1 拋撒藥的計(jì)算參數(shù)[24-25]

    Table 1 Parameters for explosives

    explosiveρ/g·cm-3D/m·s-1pCJ/GPaA/GPaB/GPaR1R2w0E0/MJ·kg-1TNT1.63693021373.8 2.74.20.90.357RDX1.69831030.15850184.61.30.3810HMX1.89991042778.3 7.14.21.00.3010.5

    Note:ρis density,Dis detonation velocity,pCJis detonation pressure at C-J state,A,B,R1,R2andw0are material consants,E0is kinetic energy per unit mass.

    a. 0 msb. 0.5 msc. 1.0 msd. 1.5 mse. 2.0 ms

    圖3 殼體應(yīng)力模擬結(jié)果(×1011Pa)

    Fig.3 Simulation results for shell stress(×1011Pa)

    a. 0 msb. 0.5 msc. 1.0 msd. 1.5 mse. 2.0 ms

    圖4 燃料拋撒與殼體破裂的試驗(yàn)結(jié)果

    Fig.4 Experimental results for fuel dispersion and shell breakage

    由圖3和圖4可知,受到拋撒藥下端起爆的影響,殼體下部最先出現(xiàn)外凸,0°~45°方向刻槽開裂,燃料沿著裂縫向外拋撒。由于計(jì)算中采用水的狀態(tài)模型代替了試驗(yàn)中使用的燃料,因此導(dǎo)致二者在殼體破裂的均勻性方面略有差異。但總體上看,模擬結(jié)果與試驗(yàn)的殼體破裂趨勢(shì)較為一致。

    針對(duì)模擬與試驗(yàn)結(jié)果,將單側(cè)燃料邊界在水平方向的膨脹速率進(jìn)行對(duì)比,如圖5所示。

    由圖5可知,在中心拋撒藥的爆轟驅(qū)動(dòng)下燃料邊界迅速膨脹; 殼體材料在應(yīng)力波作用下失效破裂燃料從裂隙中泄壓,膨脹速率下降至150 m·s-1左右; 之后由于爆轟產(chǎn)物向外傳播,邊界膨脹速率又開始出現(xiàn)上升。模擬值與試驗(yàn)結(jié)果在總體趨勢(shì)上是一致的,從而驗(yàn)證了數(shù)值模擬的正確性。

    圖5 邊界速率的模擬結(jié)果與試驗(yàn)值對(duì)比

    Fig.5 Comparison of edge velocity between simulation and experiment

    3 中心拋撒藥的影響

    3.1 拋撒藥類型

    為研究中心拋撒藥類型對(duì)拋撒效果的影響,分別針對(duì)相同質(zhì)量的梯恩梯(TNT)、黑索今(RDX)和奧克托今(HMX)進(jìn)行計(jì)算,具體參數(shù)見表1。在中心拋撒藥為TNT、與圓心距離為185 mm、刻槽間距為40 mm且刻槽深度為2 mm的條件下,計(jì)算得到各個(gè)方向的燃料速率時(shí)程曲線,圖6所示。在不同拋撒藥條件下計(jì)算得到的燃料最大拋撒速率如圖7所示,殼體應(yīng)力與破壞情況如圖8所示。

    由圖6可知,在中心拋撒藥起爆后,燃料由靜止開始加速,拋撒速率快速上升; 之后殼體受到應(yīng)力波沖擊并破裂,燃料泄壓并受到空氣阻力作用,燃料速率進(jìn)入減速階段。隨著中心拋撒藥產(chǎn)生的爆轟產(chǎn)物向外擴(kuò)散推動(dòng)燃料,燃料速率又有所上升。燃料的拋撒速率變化主要為先快速上升,后下降、振蕩,再逐漸上升趨勢(shì)。

    在不同方向上的燃料最大速率不同,180°方向最大,45°方向其次,90°方向最小。最大速率出現(xiàn)時(shí)間上0°方向最早,45°和180°方向其次,90°方向最晚。這主要是由于沖擊波達(dá)到各個(gè)方向最外層燃料的時(shí)刻以及殼體破裂的位置不同引起的。

    圖6 不同方向TNT的燃料拋撒速率時(shí)程曲線

    Fig.6 Dispersal velocity of fuel vs time for TNT at different directions

    圖7 不同拋撒藥條件下的拋撒燃料最大速率

    Fig.7 Maximum dispersal velocity of fuel for different explosives

    a. TNTb. RDXc. HMX

    圖8 不同拋撒藥條件下的殼體應(yīng)力與破壞情況(2 ms,1×1011Pa)

    Fig.8 Shell stress and breakage for different explosives(2 ms, 1×1011Pa)

    由圖7可知,中心拋撒藥對(duì)燃料初始拋撒速率的增強(qiáng)作用由大到小排序?yàn)镽DX>TNT>HMX,而由圖8可知在殼體破裂的理想程度上排序?yàn)門NT>HMX>RDX。由此可見,低爆速、低爆壓中心拋撒藥驅(qū)動(dòng)燃料拋撒的速率較低,但可保障殼體沿著預(yù)制刻槽破裂。適當(dāng)提高拋撒藥的爆速、爆壓雖然給予了燃料較高的拋撒速率,但導(dǎo)致殼體出現(xiàn)橫向斷裂,加劇殼體的破壞(扭曲、碎片等)嚴(yán)重程度。而過高的能量密度使得殼體過早破裂,燃料的最大拋撒速率反而降低了,尤其是45°和180°方向。4個(gè)方向的拋撒速率平均下降18%。這可能是由于上述兩個(gè)方向的殼體最先與沖擊波相遇并破裂,而高爆速高爆壓使得燃料從裂隙中泄壓時(shí)間提前,燃料加速過程太短引起的。

    3.2 拋撒藥位置

    由于殼體結(jié)構(gòu)的特殊性,針對(duì)不同拋撒藥(中心管)位置進(jìn)行計(jì)算。在平面幾何層面上,若使得橫截面被參考坐標(biāo)軸分割的4個(gè)“子面積”均相等,那么拋撒藥中心位置應(yīng)為189.6 mm。因此取拋撒藥中心與圓弧對(duì)應(yīng)圓心的距離分別為175, 185, 195 mm共3種情況進(jìn)行計(jì)算,獲得4個(gè)方向燃料的最大速率如圖9所示,殼體應(yīng)力與破壞情況如圖10所示。其中拋撒藥為TNT,刻槽間距為40 mm,刻槽深度為2 mm。

    圖9 不同拋撒藥中心與圓弧對(duì)應(yīng)圓心的距離條件下的燃料最大速率

    Fig.9 Maximum velocity of fuel for different distances between central explosive and circular center

    a. 175 mmb. 185 mmc. 195 mm

    圖10 不同拋撒藥中心與圓弧對(duì)應(yīng)圓心的距離條件下的殼體應(yīng)力與破壞情況(2 ms,1×1011Pa)

    Fig.10 Shell stress and breakage for different distances between central explosive and circular center(2 ms,1×1011Pa)

    在0°方向3種計(jì)算條件下的速率為110~130 m·s-1。當(dāng)拋撒藥中心與圓弧對(duì)應(yīng)圓心的距離為175 mm時(shí),180°方向?qū)?yīng)的速率最大,為231 m·s-1。隨著距離的增加,拋撒藥爆轟產(chǎn)生的沖擊波到達(dá)0°和45°方向殼體的時(shí)刻也逐漸滯后,而到達(dá)180°和90°方向的時(shí)刻不斷提前。與其他兩種情況相比,距離為185 mm的計(jì)算值在4個(gè)方向的速率始終保持較高水平。當(dāng)拋撒藥中心與圓弧對(duì)應(yīng)圓心的距離由185 mm增加為195 mm時(shí),0°,45°和90°方向燃料最大速率減小,而180°方向燃料最大速率增大。

    由圖10可知,在殼體破壞方面,隨著拋撒藥中心與圓弧對(duì)應(yīng)圓心的距離的增加,45°和180°方向殼體橫向斷裂現(xiàn)象逐漸凸顯。尤其當(dāng)距離為195 mm時(shí),0°方向殼體未沿著預(yù)制刻槽位置開裂。由此可見,拋撒藥中心與殼體各個(gè)方向周邊的距離差異不應(yīng)過大。

    4 刻槽條件的影響

    4.1 刻槽間距

    刻槽的主要目的是為了調(diào)整燃料的能量輸出,便于燃料從預(yù)制刻槽產(chǎn)生的裂隙中泄壓,同時(shí)避免殼體橫向斷裂引起燃料拋撒形成畸形云霧。下面分別針對(duì)刻槽間距為0 mm(無刻槽)、40 mm和80 mm共3種情況進(jìn)行計(jì)算。其中拋撒藥為TNT,拋撒藥與圓弧對(duì)應(yīng)圓心的距離為185 mm,刻槽深度為2 mm。獲得4個(gè)方向燃料的最大速率如圖11所示,殼體應(yīng)力與破壞情況如圖12所示。

    由圖11和圖12可知,當(dāng)刻槽間距的不同時(shí)殼體破裂位置不同,無刻槽殼體不產(chǎn)生“條形”裂隙,不同方向上燃料的最大速率也相對(duì)較小,在180°方向的殼體出現(xiàn)了“尖端”,其附近的殼體受到的應(yīng)力較大。在40 mm刻槽間距條件下,燃料拋撒最大速率較高,180°方向的速率為283 m·s-1。當(dāng)刻槽間距為80 mm時(shí),0°方向的殼體向內(nèi)彎曲,180°方向的殼體產(chǎn)生2道失效裂隙。因此較小的刻槽間距能夠?yàn)樽罱K云霧區(qū)域內(nèi)的燃料分布均勻性提供有利條件[21]。

    圖11 不同刻槽間距條件下的燃料最大速率

    Fig.11 Maximum velocity of fuel with different groove spacing for TNT

    a. 0 mmb. 40 mmc. 80 mm

    圖12 不同刻槽間距條件下的殼體應(yīng)力與破壞情況(2 ms,1×1011Pa)

    Fig.12 Shell stress and breakage with different groove spacing for TNT (2 ms,1×1011Pa)

    4.2 刻槽深度

    增加刻槽深度相當(dāng)于削弱刻槽位置的殼體強(qiáng)度,從而使得殼體在沖擊波作用下薄弱殼體處的破裂時(shí)刻提前??滩凵疃冗^小引起殼體在中心拋撒藥驅(qū)動(dòng)下發(fā)生不規(guī)則斷裂; 刻槽深度過大則可能引起殼體失穩(wěn),并在燃料裝填與保存過程中發(fā)生泄漏??紤]到殼體厚度為3 mm,針對(duì)刻槽深度分別為0 mm(無刻槽)、1 mm和2 mm共3種計(jì)算結(jié)果進(jìn)行對(duì)比。其中拋撒藥為TNT,拋撒藥與圓弧對(duì)應(yīng)圓心的距離為185 mm,刻槽間距為40 mm。4個(gè)方向燃料的最大速率如圖13所示,殼體應(yīng)力與破壞情況如圖14所示。

    圖13 不同刻槽深度條件下的燃料最大速率

    Fig.13 The maximum velocity of fuel with different groove depth for TNT

    由圖13可知,隨著刻槽深度的減小,同一方向上的燃料拋撒速率逐漸降低(0°方向除外)。其中在90°方向速率降低至102 m·s-1。而由圖14可知當(dāng)刻槽深度由2 mm降低至1 mm時(shí),殼體形變量減小、0°方向殼體出現(xiàn)應(yīng)力集中并局部橫向斷裂; 當(dāng)刻槽深度為0 mm(無刻槽)時(shí),殼體形變量進(jìn)一步減小,且在2 ms計(jì)算范圍內(nèi)殼體未出現(xiàn)破裂。所以在不穿透殼體的前提下,較大的刻槽深度一方面對(duì)燃料拋撒最大速率起到增強(qiáng)作用,另一方面可減少殼體周邊的應(yīng)力集中。

    a. 0 mmb. 1 mmc. 2 mm

    圖14 不同刻槽深度條件下的殼體應(yīng)力與破壞情況(2 ms,1×1011Pa)

    Fig.14 Shell stress and breakage with different groove depth for TNT(2 ms,1×1011Pa)

    5 結(jié) 論

    本文基于有限元分析法,對(duì)扇形殼體結(jié)構(gòu)的拋撒燃料過程進(jìn)行了數(shù)值研究,主要結(jié)論如下:

    (1) 在中心拋撒藥的爆轟驅(qū)動(dòng)下,燃料邊界膨脹速率先快速上升,后小幅下降再重新上升。與試驗(yàn)結(jié)果對(duì)比發(fā)現(xiàn),數(shù)值模擬結(jié)果符合實(shí)際情況。

    (2) 不同中心拋撒藥作用下,雖然殼體整體破壞趨勢(shì)較為一致,但在能量密度過高的拋撒藥驅(qū)動(dòng)下燃料最大速率反而降低、殼體局部橫向斷裂以及和殼體碎片增多。另外,拋撒藥的中心位置不應(yīng)過度偏離橫截面的幾何形心。

    (3) 無刻槽殼體在2 ms內(nèi)未破裂且燃料速率較低,不利于燃料拋撒; 刻槽間距40 mm的計(jì)算結(jié)果優(yōu)于80 mm和無刻槽情況。當(dāng)刻槽間距40 mm且刻槽深度為2 mm時(shí),殼體應(yīng)力集中現(xiàn)象較少,且燃料在180°方向的最大速率為283 m·s-1。

    參考文獻(xiàn):

    [1] Stefanski K, Lewandowski D, Dygdala R S, et al. Explosive formation and spreading of water-spray cloud-experimental development and model analyses[J].CentralEuropeanJournalofEnergeticMaterials, 2009, 6(3-4): 291-302.

    [2] 劉鋒. 爆炸水霧降除爆破拆除粉塵的研究[D]. 淮南: 安徽理工大學(xué), 2011.

    LIU Feng. Research on blasting water mist reducing dust resulting from blasting demolition[D]. Huainan: Anhui University of Science and Technology, 2011.

    [3] 白春華, 梁慧敏, 李建平, 等. 云霧爆轟[M]. 北京: 科學(xué)出版社, 2012: 1-6.

    [4] 張奇, 郭彥懿, 白春華, 等. 中心藥量對(duì)燃料的拋散作用[J]. 火炸藥學(xué)報(bào), 2001 (1): 17-19.

    ZHANG Qi, GUO Yan-yi, BAI Chun-hua, et al. Dispersal affection of center high explosive charge to fuel[J].ChineseJournalofExplosives&Propellants. 2001 (1): 17-19.

    [5] 薛社生, 劉家驄, 彭金華. 液體燃料爆炸拋撒過程分析[J]. 南京理工大學(xué)學(xué)報(bào), 1998, 22(1): 34-38.

    XUE She-sheng, LIU Jia-cong, PENG Jin-hua. Analysis of explosive dispersion for liquid fuel[J].JournalofNanjingUniversityofScienceandTechnology, 1998, 22(1): 34-38.

    [6] 薛社生, 劉家驄, 秦承森, 等. 燃料爆炸拋撒成霧的實(shí)驗(yàn)與數(shù)值研究[J]. 爆炸與沖擊, 2001, 21(4): 272-276.

    XUE She-sheng, LIU Jia-cong, QIN Cheng-sen, et al. Experimental and numerical investigation on explosive dispersal and cloud forming[J].ExplosionandShockWaves, 2001, 21(4): 272-276.

    [7] 蔣治海, 龍新平, 韓勇, 等. 炸藥爆炸驅(qū)動(dòng)殼體破裂及液體噴射過程試驗(yàn)研究[J]. 含能材料, 2011, 19(3): 321-324.

    JIANG Zhi-hai,LONG Xin-ping,HAN Yong, et al. Experimental studies on shell fracturing and jet forming and developing process driven by detonation[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2011, 19(3): 321-324.

    [8] 高重陽(yáng). 內(nèi)部爆炸載荷下柱殼結(jié)構(gòu)破裂問題的研究[D]. 北京: 清華大學(xué), 2000.

    GAO Chong-yang. Research on the rupture of cylindrical shell structure under internal explosive loading[D]. Beijing: Tsinghua University , 2000.

    [9] Chester R H, Rodney F R. Fracture behavior of tubular bombs[J].JournalofAppliedPhysics, 1968, 39(3): 1856-1862.

    [10] 周保順, 王少龍, 羅永鋒, 等. 殼體對(duì)燃料近區(qū)拋撒速度影響的數(shù)值模擬[J]. 火炸藥學(xué)報(bào), 2010, 33(2): 53-56.

    ZHOU Bao-shun, WANG Shao-long, LUO Yong-feng, et al. Numerical simulation on influence of shell on fuel near-field dispersal velocity[J].ChineseJournalofExplosives&Propellants, 2010, 33(2): 53-56.

    [11] 張奇, 覃彬, 白春華, 等. 中心裝藥對(duì)FAE燃料成霧特性影響的試驗(yàn)分析[J]. 含能材料, 2007, 15(5): 447-450.

    ZHANG Qi, QIN Bin, BAI Chun-hua, et al. Effection of total energy of center explosive charge on fuel dispersal characteristics feature[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2007, 15(5): 447-450.

    [12] 薛社生, 劉家驄, 彭金華. 液體燃料爆炸拋撒的近場(chǎng)階段研究[J]. 南京理工大學(xué)學(xué)報(bào), 1997, 21(4): 333-336.

    XUE She-sheng, LIU Jia-cong, PENG Jin-hua. The Solution of near-field dispersion for liquid fuel FAE[J].JournalofNanjingUniversityofScienceandTechnology, 1997, 21(4): 333-336.

    [13] 張奇, 白春華, 劉慶明, 等. 不耦合裝藥爆炸作用下燃料空氣炸藥的近區(qū)拋撒[J] . 彈道學(xué)報(bào), 2000, 12( 2): 22-25.

    ZHANG Qi, BAI Chun-hua, LIU Qing-ming, et al. Near field dispersal of FAE fuel under load of uncoupling charge in center pipe[J].JournalofBallistics, 2000, 12(2): 22-25.

    [14] Gardner D R. Near-field Dispersal Modeling for Liquid Fuel-air-explosive[R]. Sandia National Laboratories Report, SAND90-0686, 1990.

    [15] Gardner D R, Glass M W. A Coupled near-field dispersal model for fuel-air-explosives[R]. Sandia National Laboratories Report, SAND90-0687, 1990.

    [16] Glass W M. Far-field Dispersal modeling for fuel-air-explosives[R]. Sandia National Laboratories Report, SAND90-0528, 1990.

    [17] Singh S K, Singh V. Extended near-field modeling and droplet size distribution for fuel-air explosive warhead[J].DefenceScienceJournal, 2001, 51(3): 303-314.

    [18] 丁玨, 劉家驄, 彭金華. 液體燃料爆炸拋撒初期的一維數(shù)值模擬[J]. 彈道學(xué)報(bào), 1999, 11(3): 29-33.

    DING Jue, LIU Jia-cong, PENG Jin-hua. One-dimensional numerical simulation for the initial regime of liquid fuel explosive dispersal[J].JournalofBallistics, 1999, 11(3): 29-33.

    [19] 丁玨, 劉家驄. 液體燃料爆炸拋撒和FAE形成過程的數(shù)值模擬[J]. 南京理工大學(xué)學(xué)報(bào), 2000, 24(2): 168-171.

    DING Jue, LIU Jia-cong. Numerical simulation on the process of explosive dispersal for forming FAE cloud[J].JournalofNanjingUniversityofScienceandTechnology, 2000, 24(2): 168-171.

    [20] 薛社生. 液體燃料的爆炸拋撒研究[D]. 南京: 南京理工大學(xué), 1997.

    XUE She-sheng. Investigation of explosion dispersion for liquid fuels[D]. Nanjing: Nanjing University of Science and Technology, 1997.

    [21] 惠君明, 張?zhí)? 郭學(xué)永. FAE裝置參數(shù)對(duì)燃料拋撒與爆炸威力影響的實(shí)驗(yàn)研究[J]. 高壓物理學(xué)報(bào), 2004, 18(2): 103-108.

    HUI Jun-ming, ZHANG Tao, GUO Xue-yong. Experimental study on the influence of device parameters of FAE to fuel dispersion and explosive power[J].ChineseJournalofHighPressurePhysics, 2004, 18(2): 103-108.

    [22] 郭學(xué)永. 云爆戰(zhàn)斗部基礎(chǔ)技術(shù)研究[D]. 南京: 南京理工大學(xué), 2006.

    GUO Xue-yong. Basic technical research on fuel air explosive warhead[D]. Nanjing: Nanjing University of Science and Technology, 2006.

    [23] 李裕春, 時(shí)黨勇, 趙遠(yuǎn). ANSYS 11.0/LS-DYNA基礎(chǔ)理論與工程實(shí)踐[M]. 北京: 中國(guó)水利水電出版社, 2008: 280-289.

    [24] 白金澤. LS-DYNA3D基礎(chǔ)理論與實(shí)例分析[M]. 北京: 科學(xué)出版社, 2005: 89-100.

    [26] ZHANG Qi, BAI Chun-hua, LIU Qing-ming, et al. Study on near field dispersal of fuel air explosive[J].JournalofBeijingInstituteofTechnology, 1999, 8(2): 113-118.

    猜你喜歡
    殼體間距燃料
    宋代君臣殿上間距考論
    減速頂殼體潤(rùn)滑裝置的研制與構(gòu)想
    來自沙特的新燃料
    生物燃料
    導(dǎo)彈燃料知多少
    軍事文摘(2020年14期)2020-12-17 06:27:16
    高速公路指揮中心小間距LED應(yīng)用探討
    汽車變速箱殼體零件自動(dòng)化生產(chǎn)線
    油泵殼體的沖壓工藝及模具設(shè)計(jì)
    算距離
    TH型干燥筒殼體的焊接
    焊接(2016年2期)2016-02-27 13:01:22
    久久人人爽av亚洲精品天堂| 嫁个100分男人电影在线观看| 69av精品久久久久久| 夜夜爽天天搞| bbb黄色大片| 国产激情久久老熟女| 可以在线观看毛片的网站| 日韩精品中文字幕看吧| 纯流量卡能插随身wifi吗| 日本欧美视频一区| 欧美激情极品国产一区二区三区| 欧美在线一区亚洲| 操出白浆在线播放| 亚洲免费av在线视频| 午夜成年电影在线免费观看| 欧美日韩国产mv在线观看视频| 757午夜福利合集在线观看| 一边摸一边做爽爽视频免费| 亚洲成人精品中文字幕电影 | 国产精品一区二区三区四区久久 | 无遮挡黄片免费观看| 老熟妇乱子伦视频在线观看| 亚洲自偷自拍图片 自拍| netflix在线观看网站| 91精品三级在线观看| 久热这里只有精品99| 男女做爰动态图高潮gif福利片 | 男女床上黄色一级片免费看| 免费女性裸体啪啪无遮挡网站| 国产精华一区二区三区| 老司机福利观看| 国产成人免费无遮挡视频| 欧美性长视频在线观看| 精品国产亚洲在线| 久久亚洲精品不卡| 中文字幕另类日韩欧美亚洲嫩草| 国产成人av教育| 高清在线国产一区| 欧美乱色亚洲激情| 伊人久久大香线蕉亚洲五| 一a级毛片在线观看| 久久精品人人爽人人爽视色| 别揉我奶头~嗯~啊~动态视频| 成人三级做爰电影| 91在线观看av| 欧美日韩亚洲高清精品| 欧美激情 高清一区二区三区| 免费女性裸体啪啪无遮挡网站| 久99久视频精品免费| 免费av中文字幕在线| 亚洲av片天天在线观看| 国产视频一区二区在线看| 日本五十路高清| 亚洲视频免费观看视频| 亚洲五月天丁香| 亚洲成a人片在线一区二区| 欧美最黄视频在线播放免费 | 日韩欧美三级三区| 久久亚洲真实| 伦理电影免费视频| 国产亚洲av高清不卡| 美女大奶头视频| 五月开心婷婷网| 国产成人影院久久av| 婷婷丁香在线五月| 亚洲专区字幕在线| 动漫黄色视频在线观看| 一进一出好大好爽视频| 亚洲人成电影观看| 国产精品永久免费网站| 老司机在亚洲福利影院| 亚洲国产看品久久| 嫁个100分男人电影在线观看| 亚洲avbb在线观看| 国产av一区在线观看免费| 成人永久免费在线观看视频| 午夜精品国产一区二区电影| 18禁国产床啪视频网站| 亚洲一码二码三码区别大吗| 久久香蕉激情| 91精品国产国语对白视频| 久久久久久免费高清国产稀缺| 精品国产一区二区三区四区第35| а√天堂www在线а√下载| 极品人妻少妇av视频| 欧美最黄视频在线播放免费 | av在线播放免费不卡| 亚洲狠狠婷婷综合久久图片| 制服诱惑二区| 国产在线精品亚洲第一网站| 老司机深夜福利视频在线观看| 999久久久精品免费观看国产| 婷婷六月久久综合丁香| 日韩人妻精品一区2区三区| 99久久国产精品久久久| 亚洲成国产人片在线观看| 午夜免费成人在线视频| 亚洲 欧美 日韩 在线 免费| 老汉色∧v一级毛片| 男女做爰动态图高潮gif福利片 | 亚洲,欧美精品.| 大陆偷拍与自拍| 1024视频免费在线观看| 欧美乱妇无乱码| 久99久视频精品免费| 亚洲欧美激情综合另类| 国产成人欧美在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产人伦9x9x在线观看| 国产乱人伦免费视频| 亚洲欧美日韩无卡精品| 国产免费现黄频在线看| 中文字幕色久视频| 女人爽到高潮嗷嗷叫在线视频| 国产在线观看jvid| 黑丝袜美女国产一区| 久久这里只有精品19| 久久国产亚洲av麻豆专区| 黄网站色视频无遮挡免费观看| 国产在线观看jvid| 国产有黄有色有爽视频| 色综合欧美亚洲国产小说| 一级a爱视频在线免费观看| 久久国产精品男人的天堂亚洲| 黄色女人牲交| 国产一区二区三区在线臀色熟女 | 久久人人精品亚洲av| 亚洲狠狠婷婷综合久久图片| 99久久综合精品五月天人人| a级毛片黄视频| 国产一区在线观看成人免费| 久久精品91无色码中文字幕| e午夜精品久久久久久久| 中文字幕高清在线视频| 国产av在哪里看| 在线观看舔阴道视频| 在线观看免费视频日本深夜| 99热只有精品国产| 99在线视频只有这里精品首页| 怎么达到女性高潮| 久久久久亚洲av毛片大全| 在线观看免费视频网站a站| 大型av网站在线播放| 精品日产1卡2卡| 国产精品乱码一区二三区的特点 | 久热这里只有精品99| 国产精品久久久久久人妻精品电影| 91麻豆av在线| 一级黄色大片毛片| 国产精品国产av在线观看| 国产精品 国内视频| 亚洲男人的天堂狠狠| 超碰成人久久| 国产精品国产av在线观看| www.www免费av| 最近最新中文字幕大全电影3 | 国产主播在线观看一区二区| 国产精华一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看影片大全网站| av在线天堂中文字幕 | 国产av一区在线观看免费| 精品少妇一区二区三区视频日本电影| 十八禁网站免费在线| 欧美乱妇无乱码| 精品国内亚洲2022精品成人| 悠悠久久av| 国产精品二区激情视频| 日本wwww免费看| 国产成人欧美在线观看| 757午夜福利合集在线观看| 啦啦啦免费观看视频1| 国产亚洲精品久久久久5区| 欧美一级毛片孕妇| 在线观看66精品国产| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 麻豆一二三区av精品| www.自偷自拍.com| 51午夜福利影视在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 水蜜桃什么品种好| 亚洲精品粉嫩美女一区| 女人被狂操c到高潮| 老司机亚洲免费影院| 午夜福利欧美成人| 日韩有码中文字幕| 97人妻天天添夜夜摸| 国产一区二区在线av高清观看| 免费在线观看亚洲国产| 黄色成人免费大全| 美女高潮喷水抽搐中文字幕| 国产麻豆69| 久久精品亚洲精品国产色婷小说| av中文乱码字幕在线| 久久人妻熟女aⅴ| 婷婷丁香在线五月| 黄色a级毛片大全视频| 国产精品99久久99久久久不卡| 国产高清激情床上av| 国产精品久久视频播放| 国产激情欧美一区二区| 午夜视频精品福利| 免费观看人在逋| 成人三级做爰电影| 国产精品一区二区精品视频观看| 深夜精品福利| 欧美成人性av电影在线观看| 夫妻午夜视频| 亚洲成人免费电影在线观看| 美女高潮到喷水免费观看| 一个人观看的视频www高清免费观看 | 亚洲一区二区三区不卡视频| 国产精品香港三级国产av潘金莲| 久久久久九九精品影院| 亚洲在线自拍视频| 成人国产一区最新在线观看| 久久99一区二区三区| 国产片内射在线| 日本黄色视频三级网站网址| 十八禁网站免费在线| 一进一出抽搐动态| 午夜福利影视在线免费观看| 丰满的人妻完整版| 久久香蕉激情| 精品久久久精品久久久| 1024香蕉在线观看| 国产精品秋霞免费鲁丝片| 日日爽夜夜爽网站| 变态另类成人亚洲欧美熟女 | 夜夜躁狠狠躁天天躁| 丰满饥渴人妻一区二区三| 亚洲熟妇熟女久久| 日本wwww免费看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品中文字幕在线视频| 国产免费现黄频在线看| 大码成人一级视频| 亚洲人成网站在线播放欧美日韩| 亚洲狠狠婷婷综合久久图片| 精品午夜福利视频在线观看一区| 亚洲人成网站在线播放欧美日韩| 国产免费现黄频在线看| 后天国语完整版免费观看| 99久久99久久久精品蜜桃| 日本欧美视频一区| 制服诱惑二区| 又黄又粗又硬又大视频| 亚洲久久久国产精品| 久久久国产成人免费| 欧美不卡视频在线免费观看 | 色老头精品视频在线观看| 五月开心婷婷网| av天堂久久9| 搡老岳熟女国产| 在线观看日韩欧美| 国产高清videossex| 亚洲熟女毛片儿| 久久久水蜜桃国产精品网| 国产成人av教育| 91精品国产国语对白视频| 国产成人系列免费观看| 国产成人啪精品午夜网站| 亚洲伊人色综图| 欧美日韩视频精品一区| 久热爱精品视频在线9| 男女床上黄色一级片免费看| 超色免费av| 我的亚洲天堂| 国产又爽黄色视频| 波多野结衣一区麻豆| 咕卡用的链子| 99国产精品免费福利视频| 亚洲免费av在线视频| 97人妻天天添夜夜摸| 亚洲精品久久成人aⅴ小说| 51午夜福利影视在线观看| 国产欧美日韩综合在线一区二区| 亚洲欧美日韩高清在线视频| 黄网站色视频无遮挡免费观看| 久久精品aⅴ一区二区三区四区| 欧美日韩乱码在线| 女人精品久久久久毛片| 亚洲国产欧美网| ponron亚洲| 亚洲精品美女久久av网站| 国产精品偷伦视频观看了| 免费av中文字幕在线| 老司机在亚洲福利影院| 国产在线观看jvid| 99国产精品一区二区蜜桃av| 老司机午夜福利在线观看视频| 好看av亚洲va欧美ⅴa在| 国产不卡一卡二| 一级作爱视频免费观看| 中文字幕高清在线视频| 久久精品aⅴ一区二区三区四区| 久99久视频精品免费| 9热在线视频观看99| 久久久水蜜桃国产精品网| 国产亚洲欧美精品永久| 身体一侧抽搐| 日日干狠狠操夜夜爽| 国产成人av激情在线播放| 日韩成人在线观看一区二区三区| 午夜免费激情av| 99久久99久久久精品蜜桃| 一进一出好大好爽视频| 97超级碰碰碰精品色视频在线观看| 九色亚洲精品在线播放| 一个人免费在线观看的高清视频| 欧美中文综合在线视频| 男女午夜视频在线观看| 免费看十八禁软件| 一区二区三区国产精品乱码| av视频免费观看在线观看| av福利片在线| 岛国在线观看网站| 19禁男女啪啪无遮挡网站| 熟女少妇亚洲综合色aaa.| 亚洲国产欧美日韩在线播放| 中亚洲国语对白在线视频| 免费观看精品视频网站| av在线播放免费不卡| 国产成人av教育| 一级a爱视频在线免费观看| 亚洲av日韩精品久久久久久密| 成人18禁在线播放| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| av国产精品久久久久影院| 国产又色又爽无遮挡免费看| 在线观看免费午夜福利视频| 亚洲国产欧美网| 999精品在线视频| 午夜免费观看网址| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 精品久久久精品久久久| 欧美一区二区精品小视频在线| av有码第一页| 夜夜夜夜夜久久久久| 免费在线观看视频国产中文字幕亚洲| 国产99白浆流出| 深夜精品福利| 大香蕉久久成人网| 性欧美人与动物交配| 国产亚洲精品一区二区www| 搡老岳熟女国产| 亚洲一码二码三码区别大吗| 精品欧美一区二区三区在线| 亚洲国产精品合色在线| 在线观看一区二区三区| 一级,二级,三级黄色视频| 丰满人妻熟妇乱又伦精品不卡| 一二三四在线观看免费中文在| 午夜福利一区二区在线看| 亚洲七黄色美女视频| 一级片'在线观看视频| 丰满饥渴人妻一区二区三| 欧美日韩乱码在线| 国产亚洲精品综合一区在线观看 | 香蕉国产在线看| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区激情短视频| www.精华液| 亚洲精品成人av观看孕妇| 在线观看一区二区三区| 成年女人毛片免费观看观看9| 12—13女人毛片做爰片一| 99国产精品一区二区蜜桃av| 女人爽到高潮嗷嗷叫在线视频| 亚洲成av片中文字幕在线观看| 中出人妻视频一区二区| 大型黄色视频在线免费观看| 国产熟女午夜一区二区三区| 麻豆成人av在线观看| 久久人妻福利社区极品人妻图片| 日日摸夜夜添夜夜添小说| 国产精品久久久人人做人人爽| 免费看十八禁软件| 超碰成人久久| 一进一出好大好爽视频| 午夜久久久在线观看| 男女午夜视频在线观看| 一级,二级,三级黄色视频| 欧美性长视频在线观看| 91麻豆精品激情在线观看国产 | 久久狼人影院| 两人在一起打扑克的视频| 欧美亚洲日本最大视频资源| 交换朋友夫妻互换小说| 欧美人与性动交α欧美软件| 成在线人永久免费视频| 成人精品一区二区免费| 88av欧美| 老司机午夜福利在线观看视频| 亚洲五月色婷婷综合| 亚洲成人国产一区在线观看| 这个男人来自地球电影免费观看| 高清欧美精品videossex| 在线永久观看黄色视频| 在线观看午夜福利视频| av天堂久久9| 男人舔女人的私密视频| 亚洲五月天丁香| 老司机亚洲免费影院| 视频区欧美日本亚洲| 精品国产乱子伦一区二区三区| 亚洲熟妇中文字幕五十中出 | 国产精品爽爽va在线观看网站 | 日韩欧美一区二区三区在线观看| 欧美日韩一级在线毛片| 欧美日韩黄片免| 中亚洲国语对白在线视频| 国产免费男女视频| 亚洲av电影在线进入| 久久精品亚洲精品国产色婷小说| 一个人观看的视频www高清免费观看 | 国产成+人综合+亚洲专区| 久久人妻熟女aⅴ| 欧美日韩av久久| 又大又爽又粗| 欧美日韩亚洲高清精品| 欧美丝袜亚洲另类 | 久久精品影院6| 91精品国产国语对白视频| 99精国产麻豆久久婷婷| 欧美日韩瑟瑟在线播放| 视频在线观看一区二区三区| 亚洲精品国产区一区二| 亚洲欧美激情综合另类| 久久狼人影院| 久久久久久久久免费视频了| 淫秽高清视频在线观看| 欧美精品啪啪一区二区三区| 午夜免费激情av| 一进一出抽搐gif免费好疼 | 在线永久观看黄色视频| 免费在线观看日本一区| av在线天堂中文字幕 | 叶爱在线成人免费视频播放| 又黄又爽又免费观看的视频| 日韩高清综合在线| 亚洲 国产 在线| 性欧美人与动物交配| 国产1区2区3区精品| 亚洲国产看品久久| 色综合欧美亚洲国产小说| 天天躁狠狠躁夜夜躁狠狠躁| 老司机午夜十八禁免费视频| 老熟妇仑乱视频hdxx| 青草久久国产| 亚洲人成网站在线播放欧美日韩| av超薄肉色丝袜交足视频| 免费在线观看视频国产中文字幕亚洲| 精品久久久久久久毛片微露脸| 国产精品电影一区二区三区| 村上凉子中文字幕在线| 无人区码免费观看不卡| 每晚都被弄得嗷嗷叫到高潮| 国产色视频综合| 亚洲午夜精品一区,二区,三区| 99国产综合亚洲精品| 国产99久久九九免费精品| 在线免费观看的www视频| 国内久久婷婷六月综合欲色啪| 国产又色又爽无遮挡免费看| 女人被狂操c到高潮| 最好的美女福利视频网| 99热国产这里只有精品6| 亚洲精品国产区一区二| 嫁个100分男人电影在线观看| 在线观看66精品国产| 国产精品久久久久久人妻精品电影| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av成人av| 精品久久久久久久久久免费视频 | 中文字幕人妻熟女乱码| 久久久久久免费高清国产稀缺| 国产色视频综合| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 精品一区二区三区视频在线观看免费 | 国产精品偷伦视频观看了| 女同久久另类99精品国产91| 亚洲aⅴ乱码一区二区在线播放 | 啦啦啦免费观看视频1| 成人亚洲精品av一区二区 | av片东京热男人的天堂| 日韩欧美国产一区二区入口| 无人区码免费观看不卡| 午夜精品国产一区二区电影| 日韩有码中文字幕| 大型黄色视频在线免费观看| 18禁黄网站禁片午夜丰满| 久久久久久人人人人人| 无人区码免费观看不卡| 人人妻,人人澡人人爽秒播| 老鸭窝网址在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看www视频免费| 别揉我奶头~嗯~啊~动态视频| 日韩免费av在线播放| 国产成人av教育| 一边摸一边做爽爽视频免费| 一级片免费观看大全| 手机成人av网站| 19禁男女啪啪无遮挡网站| 91精品国产国语对白视频| 国产精品1区2区在线观看.| 夜夜爽天天搞| 国产亚洲精品久久久久久毛片| 亚洲国产欧美一区二区综合| 欧美人与性动交α欧美软件| 男女下面插进去视频免费观看| 亚洲国产精品合色在线| 日韩av在线大香蕉| 国产1区2区3区精品| 一个人免费在线观看的高清视频| 亚洲第一欧美日韩一区二区三区| 美女大奶头视频| 最新在线观看一区二区三区| 日本免费一区二区三区高清不卡 | 午夜精品久久久久久毛片777| 精品少妇一区二区三区视频日本电影| 午夜免费激情av| 国产99白浆流出| 精品熟女少妇八av免费久了| www.自偷自拍.com| 极品人妻少妇av视频| 黄网站色视频无遮挡免费观看| 午夜精品在线福利| 亚洲国产欧美日韩在线播放| 99久久国产精品久久久| 亚洲 国产 在线| 首页视频小说图片口味搜索| 日韩 欧美 亚洲 中文字幕| 高清在线国产一区| 日韩欧美一区二区三区在线观看| 正在播放国产对白刺激| 99热只有精品国产| 两个人看的免费小视频| 中文字幕人妻丝袜一区二区| 免费一级毛片在线播放高清视频 | 性色av乱码一区二区三区2| 亚洲自偷自拍图片 自拍| cao死你这个sao货| 国产精品98久久久久久宅男小说| 国产精品 国内视频| 亚洲中文字幕日韩| 在线观看免费日韩欧美大片| 热re99久久精品国产66热6| 久9热在线精品视频| 黄色 视频免费看| 日韩有码中文字幕| 99国产极品粉嫩在线观看| 夜夜夜夜夜久久久久| 男人舔女人下体高潮全视频| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 欧美大码av| www.熟女人妻精品国产| 亚洲片人在线观看| 最近最新中文字幕大全电影3 | 国产亚洲欧美98| 国产野战对白在线观看| 成人国语在线视频| 午夜福利欧美成人| 久久热在线av| 国产欧美日韩一区二区三| 久久午夜综合久久蜜桃| 黄片小视频在线播放| 搡老乐熟女国产| 亚洲精品美女久久av网站| 精品国产国语对白av| а√天堂www在线а√下载| www日本在线高清视频| 一夜夜www| 欧美乱码精品一区二区三区| 男女下面插进去视频免费观看| 香蕉久久夜色| 99国产精品一区二区蜜桃av| 亚洲成人免费电影在线观看| 国产精品美女特级片免费视频播放器 | 国产成人欧美| 亚洲欧美日韩无卡精品| 国产精品永久免费网站| 日韩三级视频一区二区三区| 国产亚洲精品综合一区在线观看 | 亚洲男人的天堂狠狠| 深夜精品福利| 9191精品国产免费久久| 国产成人精品无人区| 亚洲成人免费电影在线观看| 在线观看日韩欧美| 日本欧美视频一区| 怎么达到女性高潮| 香蕉丝袜av| 在线观看免费视频日本深夜| 精品一区二区三区四区五区乱码| 日本 av在线| 日日夜夜操网爽| 久久人妻福利社区极品人妻图片| 欧美国产精品va在线观看不卡| 国产成人精品无人区| 九色亚洲精品在线播放| 精品免费久久久久久久清纯| 久久精品亚洲熟妇少妇任你| 色精品久久人妻99蜜桃| 老司机亚洲免费影院| 久久精品亚洲精品国产色婷小说| av在线天堂中文字幕 | 天天躁夜夜躁狠狠躁躁| 很黄的视频免费|