• <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
    在线十欧美十亚洲十日本专区| 日韩欧美一区视频在线观看| 少妇精品久久久久久久| 午夜福利,免费看| 丰满人妻熟妇乱又伦精品不卡| 成年人午夜在线观看视频| 91大片在线观看| 免费久久久久久久精品成人欧美视频| 亚洲五月婷婷丁香| 国产av精品麻豆| 美女高潮喷水抽搐中文字幕| 午夜免费成人在线视频| 精品人妻1区二区| 国产在线视频一区二区| 亚洲人成电影观看| 国产精品久久电影中文字幕 | 国产在线精品亚洲第一网站| 午夜福利影视在线免费观看| 日韩欧美国产一区二区入口| 精品人妻在线不人妻| 久久人妻福利社区极品人妻图片| 人人妻人人澡人人看| 久久人人爽av亚洲精品天堂| 又紧又爽又黄一区二区| h视频一区二区三区| 91av网站免费观看| 精品久久久久久电影网| 啦啦啦中文免费视频观看日本| 麻豆乱淫一区二区| 在线永久观看黄色视频| 欧美精品人与动牲交sv欧美| 欧美成人免费av一区二区三区 | 成年版毛片免费区| 国产日韩欧美在线精品| 免费观看av网站的网址| 精品久久久久久电影网| 欧美日韩国产mv在线观看视频| 亚洲全国av大片| 一二三四社区在线视频社区8| 久久精品91无色码中文字幕| 亚洲欧洲日产国产| 国产精品国产高清国产av | 岛国毛片在线播放| 黄片大片在线免费观看| svipshipincom国产片| 黑人操中国人逼视频| 999久久久精品免费观看国产| 国产又爽黄色视频| 99国产精品免费福利视频| 亚洲av国产av综合av卡| 国产区一区二久久| 国产国语露脸激情在线看| 亚洲综合色网址| 老司机午夜福利在线观看视频 | 色综合婷婷激情| 国产精品久久久久成人av| 水蜜桃什么品种好| 老熟妇仑乱视频hdxx| 午夜福利,免费看| 久热爱精品视频在线9| 女人被躁到高潮嗷嗷叫费观| 啦啦啦免费观看视频1| 亚洲,欧美精品.| 成人亚洲精品一区在线观看| 日韩熟女老妇一区二区性免费视频| 纯流量卡能插随身wifi吗| 成人av一区二区三区在线看| 国产又色又爽无遮挡免费看| 国产欧美日韩精品亚洲av| 欧美日韩亚洲综合一区二区三区_| 中文字幕人妻丝袜制服| 亚洲熟女精品中文字幕| 婷婷丁香在线五月| 大码成人一级视频| 久久国产精品人妻蜜桃| 成人特级黄色片久久久久久久 | 日韩一区二区三区影片| 国产成人一区二区三区免费视频网站| 老汉色∧v一级毛片| 啦啦啦在线免费观看视频4| 日本一区二区免费在线视频| 亚洲人成77777在线视频| 欧美变态另类bdsm刘玥| 80岁老熟妇乱子伦牲交| 久久婷婷成人综合色麻豆| 欧美久久黑人一区二区| 欧美国产精品一级二级三级| 精品午夜福利视频在线观看一区 | 国产无遮挡羞羞视频在线观看| 水蜜桃什么品种好| 国产一区二区在线观看av| a级毛片在线看网站| 久久久国产精品麻豆| tube8黄色片| 亚洲va日本ⅴa欧美va伊人久久| 天堂8中文在线网| 久久国产亚洲av麻豆专区| 9热在线视频观看99| 三上悠亚av全集在线观看| 国产成人系列免费观看| 国产成人精品久久二区二区91| 国产精品免费大片| 一二三四在线观看免费中文在| 国产主播在线观看一区二区| 视频在线观看一区二区三区| 欧美久久黑人一区二区| 中文字幕制服av| 亚洲精品一二三| 青草久久国产| 亚洲第一欧美日韩一区二区三区 | 国产av国产精品国产| a级毛片黄视频| 纯流量卡能插随身wifi吗| 亚洲视频免费观看视频| 精品福利永久在线观看| 操美女的视频在线观看| 久久人人爽av亚洲精品天堂| 精品少妇一区二区三区视频日本电影| 美女国产高潮福利片在线看| 搡老乐熟女国产| 国精品久久久久久国模美| 欧美国产精品一级二级三级| www日本在线高清视频| 飞空精品影院首页| kizo精华| 久久久久久久精品吃奶| av福利片在线| 久久人人爽av亚洲精品天堂| 变态另类成人亚洲欧美熟女 | 久久久久精品人妻al黑| 国产精品一区二区在线观看99| 一边摸一边抽搐一进一小说 | 国产精品免费大片| 在线观看舔阴道视频| 免费在线观看影片大全网站| 国产精品国产av在线观看| 男女下面插进去视频免费观看| 视频区图区小说| www.熟女人妻精品国产| 人人澡人人妻人| 咕卡用的链子| 激情在线观看视频在线高清 | 国产精品久久久久久精品电影小说| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 99在线人妻在线中文字幕 | 免费在线观看影片大全网站| 精品少妇一区二区三区视频日本电影| 美女视频免费永久观看网站| 在线永久观看黄色视频| 香蕉丝袜av| 欧美精品一区二区免费开放| 午夜激情久久久久久久| 精品国产国语对白av| 精品一区二区三区视频在线观看免费 | 精品一区二区三区视频在线观看免费 | 51午夜福利影视在线观看| 久久久精品免费免费高清| 深夜精品福利| 捣出白浆h1v1| 人人妻人人澡人人爽人人夜夜| 国产精品九九99| 欧美午夜高清在线| 成在线人永久免费视频| 女人高潮潮喷娇喘18禁视频| 国产1区2区3区精品| 亚洲伊人久久精品综合| 黄网站色视频无遮挡免费观看| 精品久久久久久久毛片微露脸| 久久精品国产a三级三级三级| 色播在线永久视频| 黑人欧美特级aaaaaa片| 国产男女内射视频| 91大片在线观看| 在线十欧美十亚洲十日本专区| 成人av一区二区三区在线看| 国产日韩欧美视频二区| 在线天堂中文资源库| 色精品久久人妻99蜜桃| 狂野欧美激情性xxxx| 国产成人一区二区三区免费视频网站| 国产成+人综合+亚洲专区| a级毛片在线看网站| 天天操日日干夜夜撸| 亚洲黑人精品在线| 日韩欧美三级三区| 男女床上黄色一级片免费看| 露出奶头的视频| 久久国产精品大桥未久av| 亚洲一码二码三码区别大吗| 亚洲天堂av无毛| 国产精品欧美亚洲77777| 成人精品一区二区免费| 国产单亲对白刺激| 国产精品秋霞免费鲁丝片| 午夜久久久在线观看| 极品人妻少妇av视频| 18禁黄网站禁片午夜丰满| 国产高清激情床上av| 久热这里只有精品99| 一进一出好大好爽视频| 黄色丝袜av网址大全| 91国产中文字幕| 亚洲精品自拍成人| av超薄肉色丝袜交足视频| 亚洲成人免费电影在线观看| 亚洲全国av大片| avwww免费| 成年人午夜在线观看视频| 亚洲色图 男人天堂 中文字幕| 大陆偷拍与自拍| 欧美精品人与动牲交sv欧美| 欧美精品一区二区免费开放| 久久久久视频综合| 热99re8久久精品国产| 制服人妻中文乱码| 亚洲综合色网址| 欧美人与性动交α欧美精品济南到| 97人妻天天添夜夜摸| 人人妻人人澡人人看| 国产区一区二久久| 香蕉丝袜av| 69精品国产乱码久久久| 久久久欧美国产精品| 51午夜福利影视在线观看| 亚洲国产av影院在线观看| 亚洲色图综合在线观看| 日日摸夜夜添夜夜添小说| 成年人午夜在线观看视频| 黄色成人免费大全| 少妇 在线观看| 伊人久久大香线蕉亚洲五| 伊人久久大香线蕉亚洲五| 日韩视频在线欧美| 国产一区有黄有色的免费视频| 日韩成人在线观看一区二区三区| 巨乳人妻的诱惑在线观看| 欧美+亚洲+日韩+国产| 一个人免费看片子| 亚洲少妇的诱惑av| 国产亚洲av高清不卡| 免费黄频网站在线观看国产| 日韩精品免费视频一区二区三区| 国产一卡二卡三卡精品| 大陆偷拍与自拍| 欧美国产精品一级二级三级| 亚洲一区中文字幕在线| 亚洲情色 制服丝袜| 脱女人内裤的视频| 国产精品一区二区精品视频观看| 久久人人97超碰香蕉20202| 一本久久精品| 日韩免费高清中文字幕av| 亚洲国产成人一精品久久久| 久久精品91无色码中文字幕| 少妇裸体淫交视频免费看高清 | 最近最新免费中文字幕在线| 老司机午夜十八禁免费视频| 亚洲天堂av无毛| 久久亚洲精品不卡| 极品教师在线免费播放| 精品国产一区二区三区久久久樱花| 亚洲九九香蕉| 一级,二级,三级黄色视频| 黄片大片在线免费观看| 亚洲精品中文字幕一二三四区 | 久9热在线精品视频| 亚洲精华国产精华精| 如日韩欧美国产精品一区二区三区| 首页视频小说图片口味搜索| 国产精品亚洲av一区麻豆| 久久久精品94久久精品| 亚洲欧洲日产国产| 亚洲欧美日韩高清在线视频 | 国产成人啪精品午夜网站| 脱女人内裤的视频| 国产精品1区2区在线观看. | 99在线人妻在线中文字幕 | 国产精品久久久久成人av| 国产精品99久久99久久久不卡| 纯流量卡能插随身wifi吗| 成年动漫av网址| 亚洲精品国产精品久久久不卡| 麻豆av在线久日| 狂野欧美激情性xxxx| 国产99久久九九免费精品| 久久影院123| 国精品久久久久久国模美| netflix在线观看网站| 曰老女人黄片| 亚洲精品一卡2卡三卡4卡5卡| 999久久久精品免费观看国产| 丁香六月欧美| 搡老熟女国产l中国老女人| 热99国产精品久久久久久7| 国产亚洲av高清不卡| 国产男靠女视频免费网站| 国产不卡av网站在线观看| 啦啦啦 在线观看视频| 国产成人精品无人区| 婷婷丁香在线五月| 高清欧美精品videossex| 中文字幕制服av| 久久亚洲精品不卡| av有码第一页| 99国产精品一区二区蜜桃av | 国产高清国产精品国产三级| a级毛片在线看网站| 亚洲熟女精品中文字幕| 国产精品秋霞免费鲁丝片| 国产激情久久老熟女| 午夜日韩欧美国产| www.自偷自拍.com| 亚洲久久久国产精品| 1024视频免费在线观看| 久久国产精品大桥未久av| a在线观看视频网站| 首页视频小说图片口味搜索| 久久 成人 亚洲| 欧美成人午夜精品| 无遮挡黄片免费观看| 考比视频在线观看| 美女高潮到喷水免费观看| 肉色欧美久久久久久久蜜桃| 一夜夜www| av片东京热男人的天堂| av国产精品久久久久影院| 国产精品一区二区精品视频观看| 国产亚洲一区二区精品| 久久精品国产a三级三级三级| 天堂动漫精品| 国产精品久久久久久人妻精品电影 | 亚洲va日本ⅴa欧美va伊人久久| 亚洲色图av天堂| 九色亚洲精品在线播放| 日韩中文字幕欧美一区二区| 巨乳人妻的诱惑在线观看| 五月开心婷婷网| 国产精品久久久久久人妻精品电影 | a级毛片黄视频| 国产在线一区二区三区精| 亚洲精品粉嫩美女一区| 深夜精品福利| 国产精品久久久久久精品古装| 亚洲视频免费观看视频| 极品教师在线免费播放| 99香蕉大伊视频| av福利片在线| 少妇的丰满在线观看| 麻豆国产av国片精品| 亚洲av电影在线进入| 日韩精品免费视频一区二区三区| 国产一卡二卡三卡精品| 久久精品成人免费网站| 老司机午夜十八禁免费视频| 国产精品亚洲av一区麻豆| 午夜成年电影在线免费观看| 日本av手机在线免费观看| 黑人操中国人逼视频| 日韩大码丰满熟妇| 精品国产一区二区久久| 国产在线精品亚洲第一网站| 在线亚洲精品国产二区图片欧美| 在线观看免费视频日本深夜| 老熟妇乱子伦视频在线观看| 在线观看人妻少妇| 极品教师在线免费播放| 亚洲国产av影院在线观看| 久久国产亚洲av麻豆专区| 超碰97精品在线观看| 无遮挡黄片免费观看| 一个人免费在线观看的高清视频| 制服诱惑二区| 777米奇影视久久| 国产不卡av网站在线观看| 午夜视频精品福利| 国产日韩欧美在线精品| 国产欧美日韩精品亚洲av| 久久中文字幕一级| 人人妻人人澡人人看| 久久99热这里只频精品6学生| 色婷婷久久久亚洲欧美| 久久久久久久大尺度免费视频| 高清视频免费观看一区二区| 中亚洲国语对白在线视频| 色婷婷av一区二区三区视频| 精品国产乱码久久久久久男人| 亚洲一区二区三区欧美精品| 黄频高清免费视频| 女性生殖器流出的白浆| 亚洲成av片中文字幕在线观看| 男女免费视频国产| 18禁黄网站禁片午夜丰满| 久久热在线av| 国产精品国产高清国产av | 巨乳人妻的诱惑在线观看| 99国产综合亚洲精品| 色婷婷av一区二区三区视频| 高潮久久久久久久久久久不卡| 激情在线观看视频在线高清 | 国产成人精品无人区| 日韩中文字幕视频在线看片| 久久午夜综合久久蜜桃| 欧美成人午夜精品| 午夜老司机福利片| 午夜久久久在线观看| 黄色视频在线播放观看不卡| 动漫黄色视频在线观看| svipshipincom国产片| 99精品在免费线老司机午夜| 欧美 亚洲 国产 日韩一| 757午夜福利合集在线观看| 国产在线一区二区三区精| 亚洲男人天堂网一区| 亚洲五月婷婷丁香| 一级毛片女人18水好多| 国产精品久久久久成人av| 欧美激情 高清一区二区三区| 在线av久久热| 淫妇啪啪啪对白视频| 人人妻,人人澡人人爽秒播| 日韩大片免费观看网站| 一区二区日韩欧美中文字幕| 天堂8中文在线网| 精品人妻熟女毛片av久久网站| 国产成人免费无遮挡视频| 日韩视频一区二区在线观看| 亚洲自偷自拍图片 自拍| cao死你这个sao货| 久久久精品区二区三区| 视频区图区小说| 国产亚洲精品第一综合不卡| 午夜激情久久久久久久| 黄片小视频在线播放| 黄色a级毛片大全视频| 五月开心婷婷网| 一区二区三区乱码不卡18| 高清黄色对白视频在线免费看| 变态另类成人亚洲欧美熟女 | 男女边摸边吃奶| 青青草视频在线视频观看| 考比视频在线观看| 成人国语在线视频| 12—13女人毛片做爰片一| 一级片免费观看大全| 精品久久久精品久久久| 99精品欧美一区二区三区四区| 亚洲国产av影院在线观看| 黑人操中国人逼视频| 国产亚洲精品一区二区www | 999久久久国产精品视频| 99re6热这里在线精品视频| 人人妻人人澡人人看| 精品卡一卡二卡四卡免费| 久久久国产欧美日韩av| 免费观看av网站的网址| 国产人伦9x9x在线观看| 亚洲国产欧美日韩在线播放| 午夜福利免费观看在线| 99热网站在线观看| a级毛片黄视频| 亚洲av美国av| 欧美精品一区二区免费开放| 一区二区三区国产精品乱码| 涩涩av久久男人的天堂| 啪啪无遮挡十八禁网站| 久久青草综合色| 欧美日韩亚洲高清精品| videosex国产| 亚洲成国产人片在线观看| 又紧又爽又黄一区二区| 国产伦理片在线播放av一区| av福利片在线| 日本av手机在线免费观看| 久久久久网色| 乱人伦中国视频| 午夜福利视频在线观看免费| 欧美黑人精品巨大| 日韩人妻精品一区2区三区| av福利片在线| 午夜免费鲁丝| 欧美av亚洲av综合av国产av| 国产成人一区二区三区免费视频网站| 不卡一级毛片| 十八禁人妻一区二区| 女性被躁到高潮视频| 久久国产精品男人的天堂亚洲| 麻豆国产av国片精品| 天堂8中文在线网| 久久精品aⅴ一区二区三区四区| 桃花免费在线播放| 激情视频va一区二区三区| 亚洲精华国产精华精| 成年人午夜在线观看视频| 日韩欧美免费精品| 热re99久久国产66热| 亚洲成av片中文字幕在线观看| 亚洲精品国产色婷婷电影| 午夜免费成人在线视频| 波多野结衣一区麻豆| 亚洲五月色婷婷综合| 黄色视频,在线免费观看| 国产1区2区3区精品| 日韩三级视频一区二区三区| av天堂久久9| 巨乳人妻的诱惑在线观看| a在线观看视频网站| 亚洲熟妇熟女久久| 另类亚洲欧美激情| 日韩视频一区二区在线观看| 免费不卡黄色视频| 一边摸一边抽搐一进一出视频| 久热这里只有精品99| 国产黄频视频在线观看| 两人在一起打扑克的视频| 亚洲午夜理论影院| 国产淫语在线视频| 在线观看免费日韩欧美大片| 国产成人精品久久二区二区免费| 久久狼人影院| a级片在线免费高清观看视频| 一边摸一边抽搐一进一小说 | 丝袜喷水一区| 十分钟在线观看高清视频www| 成人av一区二区三区在线看| 日韩三级视频一区二区三区| av国产精品久久久久影院| 亚洲精品一二三| 黄色丝袜av网址大全| 人妻 亚洲 视频| 国产精品久久电影中文字幕 | 高清毛片免费观看视频网站 | 午夜91福利影院| 精品一区二区三卡| 一本色道久久久久久精品综合| 精品亚洲乱码少妇综合久久| 男女下面插进去视频免费观看| 无限看片的www在线观看| 欧美日韩亚洲综合一区二区三区_| 成人三级做爰电影| 在线观看免费视频网站a站| 精品一区二区三区四区五区乱码| 久久久久久免费高清国产稀缺| 久久久久国内视频| 国产精品免费视频内射| 久久久久久久久久久久大奶| 久久精品国产综合久久久| 亚洲中文av在线| 久久国产精品人妻蜜桃| 人妻一区二区av| 国产在线精品亚洲第一网站| 宅男免费午夜| 黑人欧美特级aaaaaa片| 欧美激情 高清一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲,欧美精品.| aaaaa片日本免费| 黄色视频在线播放观看不卡| 99精品欧美一区二区三区四区| bbb黄色大片| 另类精品久久| 国产高清视频在线播放一区| 久久中文看片网| 人人妻人人澡人人看| 蜜桃在线观看..| 18禁黄网站禁片午夜丰满| 飞空精品影院首页| 亚洲人成伊人成综合网2020| 欧美精品高潮呻吟av久久| 成年人黄色毛片网站| 性少妇av在线| 亚洲精品在线美女| 免费看十八禁软件| 不卡一级毛片| 日本vs欧美在线观看视频| 午夜福利视频在线观看免费| 法律面前人人平等表现在哪些方面| 男女床上黄色一级片免费看| 亚洲中文日韩欧美视频| 纵有疾风起免费观看全集完整版| 在线观看www视频免费| 久久这里只有精品19| av在线播放免费不卡| 国产aⅴ精品一区二区三区波| 亚洲精品美女久久久久99蜜臀| 亚洲精品久久午夜乱码| av有码第一页| 国产成人精品久久二区二区免费| 丝袜美腿诱惑在线| 女人爽到高潮嗷嗷叫在线视频| 国产男女超爽视频在线观看| 真人做人爱边吃奶动态| 午夜福利在线免费观看网站| 久久 成人 亚洲| 精品国产一区二区三区久久久樱花| 精品国产国语对白av| 亚洲精品中文字幕一二三四区 | 国产欧美日韩一区二区三| 精品一区二区三区av网在线观看 | 国产一区二区在线观看av| 搡老岳熟女国产| 国产人伦9x9x在线观看| 91精品三级在线观看| 在线十欧美十亚洲十日本专区| 国产成人av教育| 免费在线观看完整版高清| 亚洲午夜精品一区,二区,三区| 亚洲人成伊人成综合网2020| 日本精品一区二区三区蜜桃| 色婷婷久久久亚洲欧美| 亚洲国产av新网站| www.999成人在线观看| 777久久人妻少妇嫩草av网站| 国产精品秋霞免费鲁丝片| 中文字幕人妻丝袜一区二区|