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

    大型激波管產(chǎn)生的長持時沖擊波的特性研究

    2024-01-01 00:00:00程建華孟兵兵王輝頓志林張新生郭帥房
    爆破器材 2024年5期

    [摘 要] 為得到長持時沖擊波,在大型激波管高壓段2個位置預(yù)設(shè)TNT和預(yù)充高壓氣體,進(jìn)行了壓力測試,并得到了測點處的壓力時程曲線;又采用ANSYS/LS-DYNA軟件對上述工況進(jìn)行了數(shù)值仿真和參數(shù)校正;然后,對不同工況產(chǎn)生的沖擊波特性進(jìn)行仿真研究;最后,明確了不同裝藥位置、裝藥量、預(yù)充高壓氣體壓力和泄壓口長度對沖擊波特性的影響。結(jié)果表明:當(dāng)采用2處TNT裝藥,裝藥量和泄壓口長度不變時,峰值超壓和超壓時長隨著預(yù)充高壓氣體壓力的增大而增大,峰值到達(dá)時間不一定降低;當(dāng)預(yù)充高壓氣體壓力和泄壓口長度一定,2處TNT裝藥量增幅不大時,峰值超壓和峰值到達(dá)時間會降低,但超壓時間增加;同時增加TNT裝藥量和預(yù)充高壓氣體壓力,會提高峰值超壓,但對超壓時長的影響不明顯;當(dāng)TNT裝藥位置為1個,裝藥量和泄壓口長度一定時,預(yù)充高壓氣體壓力增大會增加超壓時長和峰值壓力,與2處裝藥位置時規(guī)律一致,但沖擊波的峰值到達(dá)時間降低;當(dāng)其余條件不變時,減小泄壓口長度會增加峰值超壓壓力,但是峰值到達(dá)時間卻呈先減小、后增大的趨勢。

    [關(guān)鍵詞] 激波管;長持時沖擊波;TNT;預(yù)充高壓氣體

    [分類號] O383

    Characteristics of Long Duration Shock Waves Generated by Large Shock Tubes

    CHENG Jianhua①, MENG Bingbing②, WANG Hui①, DUN Zhilin①, ZHANG Xinsheng①, GUO shuaifang②

    ① School of Civil Engineering, Henan Polytechnic University (Henan Jiaozuo, 454000)

    ② School of Resources and Environment, Henan Polytechnic University (Henan Jiaozuo, 454000)

    [ABSTRACT] In order to obtain long-duration" shock waves, TNT and pre-filled high-pressure gas were pre-set at two positions in the high-pressure section of a large shock tube. Experiments were conducted and pressure time history curves were obtained at monitoring points. ANSYS/LS-DYNA software was used to calibrate the numerical simulation parameters of the above working conditions. Then, the characteristics of shock waves generated under different working conditions were simulated and studied. Finally, the effects of different charge positions, charge weights, pre-filled gas pressures, and length of the pressure relief port on the shock wave characteristics were clarified. The results indicate that when two TNT charges are used and the charge weight and the length of the pressure relief port remain unchanged, the peak overpressure and overpressure duration increase with the increase of pre-filled high-pressure gas pressure, and the peak overpressure arrival time does not necessarily decrease. When the pre-filled high-pressure gas pressure and the length of the pressure relief port are constant, and the increase in TNT charge at two locations is not significant, the peak overpressure and peak overpressure

    arrival time will decrease, but the overpressure time will increase. Simultaneously increasing TNT charge and pre-filled high-pressure gas will increase the peak overpressure, but the effect on the duration of overpressure is not significant. When TNT is placed at one location with constant TNT charge and length of the pressure relief port, an increase in pre-filled high-pressure gas pressure will increase the duration of overpressure and peak pressure, consistent with the pattern of two loading positions, but the peak overpressure arrival time will decrease. When all the other conditions remain constant, reducing the length of the pressure relief port will increase peak overpressure, but arrival time of peak overpressure exhibits a trend of first decreasing and then increasing.

    [KEYWORDS] shock tube; long duration shock wave; TNT; pre-filled high-pressure gas

    doi:10.3969/j.issn.1001-8352.2024.05.004

    0 引言

    核爆往往伴隨著強烈的沖擊波、光輻射以及地震等破壞效應(yīng),由此會產(chǎn)生危害較大的裂隙區(qū)、振動區(qū)等[1-2]。核爆沖擊波相較于其他爆炸載荷具有長持時的特性。若想建立針對核爆沖擊波的防護工程,則有必要對長持時沖擊波的特性進(jìn)行研究。

    激波管常作為實現(xiàn)超壓沖擊波的關(guān)鍵工具。通常由高壓段、試驗段和膜片等組成,膜片安裝在高壓段和試驗段之間。驅(qū)動方式主要有高壓氣體和炸藥爆破。當(dāng)高壓段產(chǎn)生的沖擊波壓力超過膜片設(shè)定壓力時,膜片將破裂,沖擊波接著向試驗段傳播,形成激波。通過對激波管內(nèi)高壓段的驅(qū)動源和管型設(shè)計,有望模擬出長持時沖擊波。

    學(xué)者們采用現(xiàn)場試驗和數(shù)值仿真預(yù)測相結(jié)合的方法對激波管產(chǎn)生的沖擊波特性進(jìn)行了大量研究。發(fā)現(xiàn)沖擊波特性與高壓段內(nèi)的炸藥質(zhì)量[3-4]、預(yù)充氣體壓力[5-6]、膜片數(shù)量和破壞性態(tài)[7-8]、激波管截面變化[9-11]等參數(shù)息息相關(guān)。研究表明,沖擊波經(jīng)過管壁的反射和傳播后,才會逐漸形成平面沖擊波,并且峰值超壓隨著傳播距離的增加而減??;激波管的形狀(試驗段長度、直徑、角度)也會對沖擊波型產(chǎn)生影響;增大炸藥質(zhì)量也可以增加峰值超壓和波速。有的學(xué)者根據(jù)試驗數(shù)據(jù)對激波管產(chǎn)生的沖擊波超壓進(jìn)行了預(yù)測,以期為設(shè)計激波管時提供數(shù)據(jù)支撐[5,12]。另一部分學(xué)者借助激波管形成的沖擊波對面部裝備[13-14]、混凝土試樣[15-16]和其他樣品[17]等進(jìn)行了模擬爆炸沖擊試驗,為防護工程的建設(shè)提供參考。但是大多數(shù)研究的沖擊波是處于短持時階段,對于長持時沖擊波的形成以及特性的研究較少。

    本文中,通過在大型激波管內(nèi)采用TNT和預(yù)充高壓氣體驅(qū)動的方式模擬出長持時沖擊波。首先,對2處裝藥位置和預(yù)充高壓氣體壓力產(chǎn)生的沖擊波壓力時程曲線進(jìn)行監(jiān)測;然后,用ANSYS/LS-DYNA軟件對實際激波管建立有限元模型,并根據(jù)試驗結(jié)果進(jìn)行參數(shù)驗證,對沖擊波壓力云圖的演化進(jìn)行了研究;最后,研究了不同裝藥量、裝藥位置、泄壓口長度和預(yù)充氣體壓力對沖擊波特征的影響。研究成果可為實現(xiàn)長持時沖擊波提供參考。

    1 激波管內(nèi)脈沖壓力測試

    1.1 試驗設(shè)備

    試驗系統(tǒng)包括兩部分:一部分是激波管;另一部分是對沖擊波的測試系統(tǒng)。試驗激波管由高壓段、試驗段、膜片和封堵門組成,真實尺寸見圖1。封堵門用來模擬響應(yīng)結(jié)構(gòu)。其中,高壓段長度為33.00 m,最大截面直徑為3.00 m;高壓段和試驗段之間設(shè)有直徑1.18 m的膜片;試驗段內(nèi)部為空氣,長度為17.10 m,最大截面直徑為2.30 m;封堵門直徑為18.30 m,并且與試驗段之間設(shè)有長度為0.10 m的泄壓口。高壓段和試驗段的變截面尺寸見圖1。

    為了對激波管內(nèi)形成的沖擊波壓力特征進(jìn)行研究,在距泄壓口0.20 m和4.60 m處設(shè)置測點P1和P2。采用PCB-109C11傳感器對測點進(jìn)行壓力時程曲線的監(jiān)測,量程為0~690 MPa,采樣頻率為500 kHz,靈敏度為0.01 mV/kPa。采用成都泰測公司的Blast-Pro-1601爆破測試儀進(jìn)行數(shù)據(jù)采集。該設(shè)備內(nèi)置計算機和電池,不需要外部供電,采樣頻率可達(dá)1 MHz,觸發(fā)延時0~200 ms,可以滿足本次測試的需求。

    1.2 試驗方案

    為了得到預(yù)想的長持時沖擊波特征,激波管的高壓段設(shè)置1#和2#共2處TNT裝藥,并且預(yù)充高壓氣體。通過調(diào)整TNT質(zhì)量和預(yù)充高壓氣體壓力來改變沖擊波的特性,激波管高壓段、試驗段的尺寸保持不變,并使用傳感器對沖擊波特征進(jìn)行監(jiān)測。試驗方案見表1。由于TNT炸藥應(yīng)用及現(xiàn)場試驗的復(fù)雜性,共進(jìn)行2個工況的試驗,后續(xù)不同工況的沖擊波特性通過數(shù)值仿真手段研究。

    1.3 試驗結(jié)果

    分析2種工況的壓力時程曲線(圖2)。不同工況下,相同測點的壓力時程曲線形態(tài)基本一致。經(jīng)過高壓段和試驗段激波管管壁的整形,由高壓段內(nèi)TNT和高壓氣體產(chǎn)生的沖擊波壓力較為緩慢地增至峰值;隨后,又出現(xiàn)了短暫的壓力平臺期;然后,衰減至大氣壓。衰減過程中,由于試驗段管壁擴張段的影響,沖擊波波形出現(xiàn)波動,但沒有出現(xiàn)負(fù)壓區(qū)。因此可以得出,沖擊波在試驗段傳播后形成了平面波,同一橫截面上不同測點的壓力隨時間變化基本一致。

    對比分析發(fā)現(xiàn),上述2種工況均可以產(chǎn)生1 600 ms的長持時沖擊波,當(dāng)激波管的結(jié)構(gòu)不變時,增大TNT裝藥量和預(yù)充氣體壓力可以增加沖擊波的峰值壓力,但沖擊波的超壓時長基本不變。

    2 長持時沖擊波數(shù)值仿真

    2.1 數(shù)值模型的建立

    由于激波管模型龐大且具有軸對稱特性,為節(jié)省計算空間,建立1/4數(shù)值模型,同時采用cm-g-ms 單位制進(jìn)行計算。在數(shù)值模擬中,忽略膜片對流體的影響,通過合并高壓段和試驗段連接處的節(jié)點對膜片簡化,使得炸藥的爆轟產(chǎn)物和預(yù)充高壓氣體理想地從高壓段流出。為了提高計算效率和保證計算結(jié)果的可靠性,網(wǎng)格劃分過程中在模型中心軸附近加密處理。劃分單元160 568個。劃分的有限元模型見圖3。

    模型施加的邊界條件為:對所建1/4對稱模型的對稱位置施加對稱約束。即在YOZ平面上施加X方向位移約束UX=0;在XOZ平面上施加Y方向位移約束UY=0;在出口處的管壁上施加全約束,在封堵門端面上施加UZ=0的位移約束。

    2.2 炸藥模型

    由圖1得出,2處炸藥裝藥均位于激波管中心線上,分別距膜片2.70 m(2#)和5.20 m(1#)。選擇在炸藥中心起爆,空氣和炸藥均采用歐拉網(wǎng)格。TNT炸藥選擇*Mat_High_Explosive_Burn高能炸藥模型。炸藥產(chǎn)生的爆轟產(chǎn)物采用JWL狀態(tài)方程描述,即*Eos_JWL。JWL狀態(tài)方程定義了壓力和比容之間的變化關(guān)系。

    p=A1-ωρ1R1ρ0e-R1ρ0β1+B1-ωρ1R2ρ0e-R2β0A+ωEρ1。(1)

    式中:p為壓力;E為單位體積炸藥的初始內(nèi)能;A、B、R1、R2、ω、β1、β0" 均為炸藥方程參數(shù)。

    TNT的性能參數(shù)具體見表2[18]。

    2.3 預(yù)充高壓氣體和空氣狀態(tài)方程

    模型中的激波管高壓段采用炸藥和預(yù)充高壓氣體同時驅(qū)動,試驗段內(nèi)為空氣。因此,需要對預(yù)充高壓氣體和空氣分別定義。預(yù)充高壓氣體和試驗段內(nèi)的空氣均選擇*Mat_Null材料模型,二者均作為理想氣體處理,線性多項式狀態(tài)方程為*Eos_Linear_Polynomial。

    p=c0+c1μ+c2μ2+c3μ3+(c4+c5μ+c6μ2+)e0。(2)

    式中:μ=ρ/ρ0-1,ρ和ρ0分別為理想氣體的密度和初始密度;e0為初始單位體積內(nèi)能;c0~c6為與氣體性質(zhì)相關(guān)的常數(shù),對于理想氣體,c0=c1=c2=c3

    =c6=0,c4=c5=0.4 。

    對于空氣和高壓氣體,ρ0和e0的取值見表3。表3中,n為高壓氣體壓強和空氣壓強的比值。

    3 仿真結(jié)果

    3.1 數(shù)值仿真驗證

    為驗證建立的激波管有限元模型預(yù)測沖擊波特性的可行性,對工況1-1和1-2分別在測點P1和P2進(jìn)行數(shù)值模擬結(jié)果和試驗數(shù)據(jù)的對比,測得的壓力時程曲線以及峰值壓力見圖4和表4。由圖4和表4可知,2種工況在測點的試驗曲線和模擬曲線基本一致,峰值超壓的試驗值和模擬值的相對誤差小于15%,表明該有限元模型、高壓氣體、TNT等參數(shù)的選取和網(wǎng)格的劃分均具有科學(xué)性,可以用來預(yù)測不同工況下激波管產(chǎn)生的沖擊波超壓曲線。" 從圖4還可得出,工況1-1和工況1-2下產(chǎn)生的沖擊波的超壓時間達(dá)到了1 600 ms,實現(xiàn)了長持時沖擊波。在工況1-1下,P2處的峰值壓力達(dá)到了208.21 kPa,通過增加1#和2#位置的TNT裝藥量和預(yù)充高壓氣體的壓力(工況1-2),可以有效提高超壓的峰值壓力(278.31 kPa)和整體壓力。同時發(fā)現(xiàn),仿真得到的沖擊波波形更加復(fù)雜,二次峰值略高于試驗數(shù)據(jù)。這是由于仿真過程為理想條件,高能流體不會在試驗段泄漏,會不斷在激波管壁面反射疊加,這對沖擊波壓力的疊加是有利的;但在試驗過程中,無法達(dá)到仿真的試驗環(huán)境。同時,發(fā)現(xiàn)峰值到達(dá)時間的試驗值提前于模擬值,證明試驗段內(nèi)的激波速度的試驗值高于模擬值??赡苁菍げü軐嶋H剖面圖進(jìn)行了簡化處理導(dǎo)致的。由于測量過程中出現(xiàn)了傳感器誤差,導(dǎo)致圖4(b)工況1-2的模擬值和試驗值均出現(xiàn)了較小的誤差與雜波??赡苁怯捎谠摲N工況下裝藥量和預(yù)充高壓氣體壓力較大,對傳感器的固定起到了一些干擾。

    3.2 沖擊波壓力和等壓面演化

    對激波管高壓段中的高能流體產(chǎn)生的壓力場演化進(jìn)行研究,工況1-1的壓力云圖和等壓面云圖演化見圖5。

    由圖5(a)看出:當(dāng)t=0 ms時,高壓段內(nèi)壓力為0.42 MPa,證明預(yù)充高壓氣體成功;而隨著時間推進(jìn),TNT開始起爆,爆轟氣體和預(yù)充高壓氣體相互作用,產(chǎn)生在高壓段傳播的稀疏波和向試驗段演化的激波,沖破膜片后,在試驗段形成較低壓和低速激波;并且發(fā)現(xiàn)稀疏波在向激波轉(zhuǎn)化的過程中,兩者之間出現(xiàn)了接觸面,根據(jù)相容關(guān)系,接觸面的兩側(cè)壓力相等;通過觀察高壓段內(nèi)的稀疏波特征得出,稀疏波經(jīng)過壁面的不斷反射,逐漸追趕激波管軸線傳播的球形波陣面;隨著時間推進(jìn),稀疏波前沿的波形更加復(fù)雜,但也會慢慢被激波管管壁整形成較為平整的平面波。由此可見,激波管產(chǎn)生的沖擊波不會立即生成平面波,而是通過一定距離演化后,才會形成等壓平面波。

    平面波的演化過程可由圖5(b)更清楚地表現(xiàn)。當(dāng)t=1.0 ms時,2團由炸藥爆炸產(chǎn)生的高能流體在該處形成了等壓面;當(dāng)t=6.0 ms時,高壓段內(nèi)炸藥和隔膜片之間形成了等壓沖擊波,隨著時間延長,等壓面逐漸在高壓段內(nèi)形成,并向試驗段傳遞;當(dāng)t=926.0 ms時,高壓段內(nèi)已經(jīng)沒有沖擊波等壓面,這可能是由于沖擊波在激波管內(nèi)的反射作用造成的,而沖擊波通過試驗段內(nèi)的整形,慢慢形成等壓平面波,依次向泄壓口推進(jìn)。

    3.3 預(yù)充高壓氣體壓力和TNT質(zhì)量對沖擊波特性的影響

    為了探究不同預(yù)充高壓氣體壓力和TNT質(zhì)量對超壓時程的影響,對工況2-1、2-2和2-3分別進(jìn)行了仿真。仿真方案見表5。測點P2處得到的壓力時程曲線如圖6所示。對比2個裝藥位置下,5種工況產(chǎn)生的沖擊波特性,見圖7。為提高仿真效率,將仿真時間縮短至600 ms,并根據(jù)600 ms的超壓判斷超壓時間。在600 ms時,超壓越大,則相應(yīng)的超壓時間越長。

    對比圖2中的工況1-1和圖6中的工況2-2的仿真壓力時程曲線發(fā)現(xiàn),相同TNT質(zhì)量下(1#和2#處均是400 g TNT),降低高壓段內(nèi)預(yù)充氣體的壓力,峰值壓力降低,峰值到達(dá)時間延長,證明激波波速有所降低。對比工況2-1和2-3發(fā)現(xiàn),在低TNT裝藥量前提下(1#和2#處裝藥300 g),預(yù)充高壓氣體壓力越低,試驗段內(nèi)激波傳播速度和峰值壓力越低,超壓時間減少。證明在泄壓口和試驗段長度等其他參數(shù)不變的情況下,高壓段裝藥量和預(yù)充高壓氣體壓力均會影響超壓時間。

    3.4 預(yù)充高壓氣體壓力和TNT裝藥位置對沖擊波特性的影響

    為了研究單一裝藥位置條件下,預(yù)充高壓氣體壓力對沖擊波特性的影響,控制在1#處裝藥250 g,對預(yù)充氣體壓力進(jìn)行調(diào)節(jié)。 具體方案見表6。

    單一裝藥位置下,不同預(yù)充壓力的沖擊波特性見圖8。并且提取了沖擊波的峰值壓力和峰值到達(dá)時間,見圖9??梢钥闯?,在相同裝藥量的情況下,沖擊波峰值壓力隨著預(yù)充氣體壓力的增加而增大,工況3-2的最大峰值壓力達(dá)到189.0 kPa。3種工況下的沖擊波波形基本一致,證明在其他條件不變的情況下,改變預(yù)充高壓氣體壓力不會影響沖擊波的形態(tài),只會改變沖擊波的峰值壓力。當(dāng)600 ms時,工況3-2的沖擊波壓力最大,相應(yīng)地,該工況的超壓時長也會最大。同時發(fā)現(xiàn),峰值到達(dá)時間與預(yù)充高壓氣體壓力負(fù)相關(guān),預(yù)充高壓氣體壓力越大,峰值到達(dá)時間越長,沖擊波傳播速度越慢??赡苁穷A(yù)充高壓氣體壓力越高,容易在激波管內(nèi)形成壅流,反而影響了沖擊波傳播速度,但超壓持續(xù)時間增加。

    3.5 泄壓口長度對沖擊波特性的影響

    為探究不同泄壓口長度對沖擊波特性的影響,制定了如表7的試驗方案。提取的壓力時程曲線、峰值壓力與峰值到達(dá)時間見圖10和圖11??梢钥闯?,在裝藥量和預(yù)充高壓氣體一致的情況下,3種工況下的沖擊波壓力時程曲線也基本相同,峰值到達(dá)時間和峰值壓力差別較大。隨著泄壓口長度的增大,沖擊波的峰值壓力逐漸降低,但峰值到達(dá)時間呈現(xiàn)先減小、后增大的趨勢,證明沖擊波的波速演化也呈相同規(guī)律。同時還發(fā)現(xiàn),仿真時間為600 ms時的沖擊波壓力與峰值壓力的順序一致,證明超壓時間與泄壓口長度呈負(fù)相關(guān),泄壓口長度越大,超壓時長越短。

    4 結(jié) 論

    采用大型激波管來實現(xiàn)長持時沖擊波,運用試驗和數(shù)值仿真相結(jié)合的方式,對不同裝藥位置、裝藥量、泄壓口長度和預(yù)充高壓氣體壓力產(chǎn)生的沖擊波特性進(jìn)行了研究,得出以下結(jié)論:

    1)當(dāng)裝藥位置為2個,裝藥量和泄壓口長度不變時,峰值超壓和超壓時長隨著預(yù)充高壓氣體壓力增大而增大,峰值到達(dá)時間可能會增加,也可能縮短;當(dāng)預(yù)充高壓氣體壓力和泄壓口長度一定,2處TNT裝藥量增幅不大時,峰值超壓和峰值到達(dá)時間會降低,但超壓時間增加。

    2)當(dāng)裝藥位置為1個且裝藥量一定時,超壓時長和超壓峰值壓力隨著預(yù)充高壓氣體壓力的增加而增大,但沖擊波的峰值到達(dá)時間卻呈現(xiàn)相反規(guī)律,證明沖擊波的傳播速度與預(yù)充高壓氣體壓力負(fù)相關(guān)。

    3)當(dāng)裝藥位置、裝藥量和預(yù)充高壓氣體壓力一定時,沖擊波的峰值壓力隨著試驗段內(nèi)的泄壓口變小而增大,但峰值到達(dá)時間卻呈現(xiàn)先減小、后增大的趨勢,證明沖擊波的傳播速度與泄壓口的長度存在最佳匹配參數(shù)。

    研究了影響長持時沖擊波特性的TNT裝藥量和預(yù)充高壓氣體壓力,取得了基本認(rèn)識,其他影響因素(激波管截面尺寸、多個裝藥位置)對沖擊波特性的影響尚不明確,以及該激波管模型對防護工程的響應(yīng)規(guī)律尚未研究,可作為下一步研究方向。

    參考文獻(xiàn)

    [1] 郭乙木, 楊慶大. 核爆炸條件下地沖擊波傳播的巖體阻尼研究[J]. 巖土力學(xué), 1995(3): 55-61.

    GUO Y M, YANG Q D. Study on rock damping in case of nuclear explosive wave propagation [J]. Rock and Soil Mechanics, 1995(3): 55-61.

    [2] 任輝啟, 黃魁, 吳祥云, 等. 地面目標(biāo)空氣沖擊波動壓毀傷研究進(jìn)展[J]. 防護工程, 2021, 43(1): 1-9.

    REN H Q, HUANG K, WU X Y, et al. Research progress of dynamic pressure damage to ground targets by air shock wave [J]. Protective Engineering, 2021, 43(1): 1-9.

    [3] 任輝啟, 王世合, 周松柏, 等. 大型爆炸波模擬裝置研制及其應(yīng)用[C]//第十六屆全國激波與激波管學(xué)術(shù)會議論文集. 洛陽, 2014: 10-22.

    REN H Q, WANG S H, ZHOU S B, et al. The development and application of large blast wave simulator [C]//The 16th National Conference on Shock Waves and Shock Tubes. Luoyang, 2014: 10-22.

    [4] 張軍, 黃含軍, 王軍評, 等. 炸藥驅(qū)動式爆炸管的載荷計算[J]. 裝備環(huán)境工程, 2021, 18(5): 21-27.

    ZHANG J, HUANG H J, WANG J P, et al. Simulation on the blast load inside the explosively drived shock tube[J]. Equipment Environmental Engineering, 2021, 18(5): 21-27.

    [5] 周岳蘭, 裴魯, 龍仁榮, 等. 激波管內(nèi)壓力脈沖演化特性及模擬空爆沖擊波的方法研究[J]. 兵工學(xué)報, 2023, 44(12): 3815-3825.

    ZHOU Y L, PEI L, LONG R R, et al. Study on the evolution characteristics of pressure pulse in shock tube and method of simulating air explosion shock wave [J]. Acta Armamentarii, 2023, 44(12): 3815-3825.

    [6] 張坤玉, 陳德, 吳昊. 高壓氣體驅(qū)動激波管的數(shù)值模擬與參數(shù)影響分析[J]. 高壓物理學(xué)報, 2023, 37(3): 7-17.

    ZHANG K Y, CHEN D, WU H. Numerical simulation and parametric analysis of high-pressure gas-driven shock tube [J]. Chinese Journal of High Pressure Physics, 2023, 37(3): 7-17.

    [7] 汪維, 任輝啟, 康冬麗, 等. 大型爆炸波模擬裝置大氣室多膜片驅(qū)動器驅(qū)動數(shù)值模擬研究[C]//第十六屆全國激波與激波管學(xué)術(shù)會議." 洛陽, 2014: 605-612.

    WANG W, REN H Q, KANG D L, et,al. Numerical simulation research of large blast wave simulator loaded by multiple diaphragm atmospheric chamber driver [C]//The 16th National Conference on Shock Waves and Shock Tubes. Luoyang, 2014: 605-612.

    [8] 何起光, 張偉, 陳小偉, 等. 激波管聚酯膜片變形過程分析[J]. 爆炸與沖擊, 2019, 39(3): 33201.

    HE Q G, ZHANG W, CHEN X W, et al. Analysis on the deformation process of PET shock tube diaphragm[J]. Explosion and Shock Waves, 2019, 39(3): 33201.

    [9] 鄭監(jiān), 盧芳云, 陳榮. 柱形裝藥條件下錐形水中爆炸激波管內(nèi)的沖擊波特性[J]. 爆炸與沖擊, 2021, 41(10): 75-86.

    ZHENG J, LU F Y, CHEN R. Shock wave characteristics in a conical water explosion shock tube under cylindrical charge condition [J]. Explosion and Shock Waves, 2021, 41(10): 75-86.

    [10] 楊軍, 薛斌. 激波管管長對階躍壓力波形的影響分析[J]. 振動與沖擊, 2019, 38(3): 252-257.

    YANG J, XUE B. Effects of shock tube length on step pressure waveform [J]. Journal of Vibration and Shock, 2019, 38(3): 252-257.

    [11] STEWART J B. Approximating a free-field blast environment in the test section of an explosively driven conical shock tube [J]. Shock Waves, 2019, 29(2): 355-360.

    [12] 白世杰, 梁興雨, 王昆. 用于激波增強的變截面激波管設(shè)計[J]. 燃燒科學(xué)與技術(shù), 2023, 29(4): 406-413.

    BAI S J, LIANG X Y, WANG K. Design of variable cross section for shock enhancement in shock tubes [J]. Journal of Combustion Science and Technology, 2023, 29(4): 406-413.

    [13] 康越, 張仕忠, 張遠(yuǎn)平, 等. 基于激波管評價的單兵頭面部裝備沖擊波防護性能研究[J]. 爆炸與沖擊, 2021, 41(8): 176-188.

    KANG Y, ZHANG S Z, ZHANG Y P, et al. Research on anti-shockwave performance of the protective equipment for the head of a soldier based on shock tube evaluation [J]. Explosion and Shock Waves, 2021, 41(8): 176-188.

    [14] 蔡志華, 賀葳, 汪劍輝, 等. 爆炸波致顱腦損傷力學(xué)機制與防護綜述[J]. 兵工學(xué)報, 2022, 43(2): 467-480.

    CAI Z H, HE W, WANG J H, et al. Review on mechanical mechanism of blast-induced traumatic brain injury and protection technology [J]. Acta Armamentarii, 2022, 43(2): 467-480.

    [15] LI Y, AOUDE H. Influence of steel fibers on the static and blast response of beams built with high-strength concrete and high-strength reinforcement[J]. Engineering Structures, 2020, 221: 111031.

    [16] 陳德, 吳昊, 徐世林, 等. 單向砌體填充墻激波管試驗和動力行為分析[J]. 爆炸與沖擊, 2023, 43(8): 134-152.

    CHEN D, WU H, XU S L, et al. Shock tube tests and dynamic behavior analyses on one-way masonry-infilled walls [J]. Explosion and Shock Waves, 2023, 43(8): 134-152.

    [17] 田銳, 魏剛, 張涵哲, 等. 基于空氣激波管的GFRP層合板抗爆性能研究[J]. 復(fù)合材料科學(xué)與工程, 2023(3): 68-75.

    TIAN R, WEI G, ZHANG H Z, et al. Study of the blast resistance of GFRP laminates based on air shock tube [J]. Composites Science and Engineering, 2023(3): 68-75.

    [18] 宋浦, 楊凱, 梁安定, 等. 國內(nèi)外TNT炸藥的JWL狀態(tài)方程及其能量釋放差異分析[J]. 火炸藥學(xué)報, 2013, 36(2): 42-45.

    SONG P, YANG K, LIANG A D, et al. Difference analysis on JWL-EOS and energy release of different TNT charge[J]. Chinese Journal of Explosives amp; Propellants, 2013, 36(2): 42-45.

    美女被艹到高潮喷水动态| 久久久精品欧美日韩精品| 大型黄色视频在线免费观看| 欧美bdsm另类| 美女大奶头视频| 美女黄网站色视频| 美女xxoo啪啪120秒动态图| 一级毛片电影观看 | 三级经典国产精品| 日韩高清综合在线| 欧美成人免费av一区二区三区| 国产精品免费一区二区三区在线| 成年女人看的毛片在线观看| 国产成人精品久久久久久| 欧美日本亚洲视频在线播放| 国产精品人妻久久久影院| 中文字幕av成人在线电影| 久久精品国产自在天天线| 男人舔女人下体高潮全视频| 91久久精品电影网| 亚洲一区二区三区色噜噜| 亚洲精品日韩av片在线观看| 校园春色视频在线观看| 国产精品久久视频播放| 在线免费十八禁| 精品午夜福利视频在线观看一区| 床上黄色一级片| 亚洲欧美精品自产自拍| 直男gayav资源| 成人国产麻豆网| 久久热精品热| 成人特级黄色片久久久久久久| 伦理电影大哥的女人| 熟女电影av网| 变态另类成人亚洲欧美熟女| 国产精品电影一区二区三区| 久久久久久久久久成人| 国产精品电影一区二区三区| .国产精品久久| 国产午夜精品久久久久久一区二区三区 | 亚洲中文日韩欧美视频| 婷婷精品国产亚洲av| 大型黄色视频在线免费观看| 日韩制服骚丝袜av| 看非洲黑人一级黄片| 亚洲精品色激情综合| 观看美女的网站| 国产精华一区二区三区| 成人无遮挡网站| 日韩欧美免费精品| а√天堂www在线а√下载| 激情 狠狠 欧美| 国产欧美日韩精品亚洲av| 美女黄网站色视频| 欧美精品国产亚洲| 欧美日本视频| 搞女人的毛片| 看免费成人av毛片| 又黄又爽又免费观看的视频| 精品日产1卡2卡| 99久久九九国产精品国产免费| 国内少妇人妻偷人精品xxx网站| 热99在线观看视频| 99国产精品一区二区蜜桃av| 日本撒尿小便嘘嘘汇集6| 啦啦啦观看免费观看视频高清| h日本视频在线播放| 国产三级在线视频| 搡女人真爽免费视频火全软件 | 一本一本综合久久| 国产三级在线视频| 国产精品久久久久久精品电影| 国产成人福利小说| 美女内射精品一级片tv| 国产视频内射| 国产aⅴ精品一区二区三区波| 国产又黄又爽又无遮挡在线| 午夜精品一区二区三区免费看| 熟妇人妻久久中文字幕3abv| 亚洲激情五月婷婷啪啪| 国内揄拍国产精品人妻在线| 嫩草影院精品99| 男人的好看免费观看在线视频| 内地一区二区视频在线| av专区在线播放| 午夜视频国产福利| 精品少妇黑人巨大在线播放 | 99在线人妻在线中文字幕| 国产一区亚洲一区在线观看| 六月丁香七月| 麻豆av噜噜一区二区三区| 国产探花在线观看一区二区| 尾随美女入室| 日本-黄色视频高清免费观看| 日韩制服骚丝袜av| 国产黄a三级三级三级人| 久久午夜福利片| 乱系列少妇在线播放| 精品人妻一区二区三区麻豆 | 久久久午夜欧美精品| 国产精品av视频在线免费观看| 99九九线精品视频在线观看视频| 亚洲国产精品成人久久小说 | 精品午夜福利视频在线观看一区| 亚洲国产精品久久男人天堂| 99热网站在线观看| 免费观看人在逋| 成年女人毛片免费观看观看9| av专区在线播放| 黄色一级大片看看| 99久久久亚洲精品蜜臀av| 午夜视频国产福利| 激情 狠狠 欧美| 欧美在线一区亚洲| 欧美性感艳星| 日本撒尿小便嘘嘘汇集6| 哪里可以看免费的av片| 午夜福利在线观看吧| 亚洲aⅴ乱码一区二区在线播放| 免费观看人在逋| 国内精品宾馆在线| 热99re8久久精品国产| 久久久精品大字幕| 亚洲欧美精品自产自拍| 国产蜜桃级精品一区二区三区| 成人午夜高清在线视频| 日本欧美国产在线视频| 精品久久久久久久久久久久久| 国内精品一区二区在线观看| 亚洲精品一区av在线观看| 久久精品国产亚洲av涩爱 | 99视频精品全部免费 在线| 美女被艹到高潮喷水动态| 变态另类丝袜制服| 亚洲精品国产成人久久av| 久久久久九九精品影院| 国产人妻一区二区三区在| 日韩欧美国产在线观看| 国产高清不卡午夜福利| 午夜a级毛片| 超碰av人人做人人爽久久| 免费在线观看成人毛片| av专区在线播放| ponron亚洲| 亚洲熟妇中文字幕五十中出| 久久天躁狠狠躁夜夜2o2o| 午夜爱爱视频在线播放| 午夜精品国产一区二区电影 | 五月玫瑰六月丁香| 在线观看美女被高潮喷水网站| 国产视频内射| 精品人妻一区二区三区麻豆 | 两性午夜刺激爽爽歪歪视频在线观看| 成人精品一区二区免费| 国产精品野战在线观看| 免费大片18禁| 久久综合国产亚洲精品| 丰满的人妻完整版| 久久精品91蜜桃| 午夜福利在线在线| 亚洲自偷自拍三级| 国产欧美日韩一区二区精品| 免费观看在线日韩| 免费人成在线观看视频色| 免费高清视频大片| 在线国产一区二区在线| 一进一出抽搐gif免费好疼| 人人妻人人澡欧美一区二区| 国产精品日韩av在线免费观看| 精品一区二区三区视频在线| 十八禁国产超污无遮挡网站| 色5月婷婷丁香| 亚洲美女搞黄在线观看 | 国产免费男女视频| 国产高清视频在线播放一区| 99在线人妻在线中文字幕| 成人无遮挡网站| 成年女人毛片免费观看观看9| 日韩,欧美,国产一区二区三区 | 啦啦啦韩国在线观看视频| 日本-黄色视频高清免费观看| 色综合亚洲欧美另类图片| 插逼视频在线观看| 97超碰精品成人国产| 久99久视频精品免费| 熟女人妻精品中文字幕| 亚洲中文字幕日韩| 桃色一区二区三区在线观看| 人妻少妇偷人精品九色| 欧美zozozo另类| avwww免费| 亚洲精品影视一区二区三区av| 麻豆乱淫一区二区| 日日摸夜夜添夜夜添av毛片| 在线观看av片永久免费下载| 色5月婷婷丁香| 国产三级中文精品| 久久久久国内视频| 亚洲精品久久国产高清桃花| 国产极品精品免费视频能看的| 精品99又大又爽又粗少妇毛片| 欧美性猛交╳xxx乱大交人| 中文字幕av成人在线电影| 男人舔女人下体高潮全视频| 精品久久久久久久人妻蜜臀av| 国产极品精品免费视频能看的| 亚洲一区二区三区色噜噜| av在线亚洲专区| 午夜福利视频1000在线观看| 亚洲自偷自拍三级| 草草在线视频免费看| 亚洲性久久影院| 亚洲人成网站在线观看播放| 国模一区二区三区四区视频| av卡一久久| 国产高潮美女av| 日本撒尿小便嘘嘘汇集6| 国产精品久久视频播放| 中文字幕熟女人妻在线| 免费观看在线日韩| avwww免费| 淫秽高清视频在线观看| 亚洲成a人片在线一区二区| 国产又黄又爽又无遮挡在线| 波野结衣二区三区在线| 色综合色国产| 国产不卡一卡二| 国产精品一区二区三区四区久久| 国产高清视频在线播放一区| 亚洲自偷自拍三级| 男人的好看免费观看在线视频| 老女人水多毛片| 一个人看的www免费观看视频| 亚洲av免费在线观看| 日本一二三区视频观看| 精品欧美国产一区二区三| av国产免费在线观看| 人妻制服诱惑在线中文字幕| 女生性感内裤真人,穿戴方法视频| 床上黄色一级片| 麻豆精品久久久久久蜜桃| 中文字幕av在线有码专区| 日日摸夜夜添夜夜添小说| 中文字幕av成人在线电影| 欧美性猛交黑人性爽| 亚洲国产欧美人成| 九九久久精品国产亚洲av麻豆| 成人二区视频| 一级黄片播放器| 国产成人福利小说| 韩国av在线不卡| 性欧美人与动物交配| 午夜爱爱视频在线播放| 国产亚洲91精品色在线| 97在线视频观看| 真实男女啪啪啪动态图| 深爱激情五月婷婷| 97人妻精品一区二区三区麻豆| 亚洲成人中文字幕在线播放| 亚洲一区二区三区色噜噜| 亚洲最大成人av| 国产成人影院久久av| 日韩大尺度精品在线看网址| 国产精品亚洲一级av第二区| 我的女老师完整版在线观看| 高清毛片免费观看视频网站| 亚洲国产欧洲综合997久久,| 看非洲黑人一级黄片| 超碰av人人做人人爽久久| 一级av片app| 在线观看免费视频日本深夜| 亚洲电影在线观看av| 国产高清视频在线观看网站| 一进一出抽搐gif免费好疼| 色吧在线观看| 免费人成视频x8x8入口观看| 男人狂女人下面高潮的视频| 国产美女午夜福利| 中文资源天堂在线| 你懂的网址亚洲精品在线观看 | 免费av不卡在线播放| 日韩av在线大香蕉| 久久久久久久久久黄片| 国产精品久久视频播放| 国产aⅴ精品一区二区三区波| 成人综合一区亚洲| 嫩草影视91久久| 深爱激情五月婷婷| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 久久久久久国产a免费观看| 国产日本99.免费观看| 成年版毛片免费区| 人人妻人人澡欧美一区二区| 女生性感内裤真人,穿戴方法视频| 欧美日韩在线观看h| 国产精品爽爽va在线观看网站| 久久国内精品自在自线图片| 淫妇啪啪啪对白视频| 精品一区二区三区人妻视频| 综合色丁香网| 精品午夜福利视频在线观看一区| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 国产久久久一区二区三区| 亚洲成人久久性| 免费大片18禁| 久久久a久久爽久久v久久| 人妻夜夜爽99麻豆av| 床上黄色一级片| 久久久色成人| 成人综合一区亚洲| 男人舔奶头视频| 成年女人永久免费观看视频| 欧美丝袜亚洲另类| 18+在线观看网站| 男人舔女人下体高潮全视频| 一级毛片我不卡| 亚洲性久久影院| 亚洲av免费在线观看| 狂野欧美白嫩少妇大欣赏| 成人一区二区视频在线观看| 国产一区二区在线观看日韩| 1024手机看黄色片| 美女被艹到高潮喷水动态| 国产爱豆传媒在线观看| 亚州av有码| 美女 人体艺术 gogo| 国产激情偷乱视频一区二区| 啦啦啦韩国在线观看视频| www日本黄色视频网| 男人的好看免费观看在线视频| 18+在线观看网站| 一个人看的www免费观看视频| 国产黄a三级三级三级人| 色视频www国产| 国产亚洲精品av在线| 亚洲精品在线观看二区| 99久久成人亚洲精品观看| 欧美成人免费av一区二区三区| 亚洲欧美精品自产自拍| 午夜亚洲福利在线播放| 精品久久久久久久久久免费视频| 91久久精品电影网| 综合色av麻豆| 晚上一个人看的免费电影| 欧美一区二区亚洲| 最近在线观看免费完整版| 久久精品国产亚洲av香蕉五月| 狂野欧美白嫩少妇大欣赏| 狂野欧美激情性xxxx在线观看| 天堂影院成人在线观看| 无遮挡黄片免费观看| 99riav亚洲国产免费| 亚洲在线自拍视频| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 免费观看在线日韩| 亚洲专区国产一区二区| 亚洲无线观看免费| 欧美高清性xxxxhd video| 日韩人妻高清精品专区| 亚洲人成网站在线观看播放| 国产黄色视频一区二区在线观看 | 久久久色成人| 99热网站在线观看| 亚洲av免费高清在线观看| 午夜福利视频1000在线观看| 亚洲婷婷狠狠爱综合网| 国产精品久久久久久久久免| 搡老熟女国产l中国老女人| 在线免费观看的www视频| 精品欧美国产一区二区三| 成人高潮视频无遮挡免费网站| 日韩中字成人| 国产麻豆成人av免费视频| 天天一区二区日本电影三级| 国国产精品蜜臀av免费| 久久人人爽人人片av| 三级国产精品欧美在线观看| 哪里可以看免费的av片| 一级a爱片免费观看的视频| av女优亚洲男人天堂| 又黄又爽又免费观看的视频| 中文字幕av在线有码专区| 亚洲国产精品sss在线观看| 日韩 亚洲 欧美在线| 亚洲专区国产一区二区| 中国国产av一级| 国产美女午夜福利| 久久6这里有精品| 亚洲av免费高清在线观看| 91精品国产九色| av在线蜜桃| 欧美在线一区亚洲| av专区在线播放| 精品久久久久久久久亚洲| 欧美xxxx性猛交bbbb| 永久网站在线| 色综合站精品国产| 亚洲第一电影网av| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看 | 观看美女的网站| 亚洲在线自拍视频| 日韩精品有码人妻一区| 久久精品国产亚洲网站| 国产免费一级a男人的天堂| av卡一久久| 日韩欧美精品免费久久| 亚洲中文日韩欧美视频| 国产精品免费一区二区三区在线| 欧美中文日本在线观看视频| 成人特级黄色片久久久久久久| 一级毛片久久久久久久久女| 久久久a久久爽久久v久久| 亚洲av中文字字幕乱码综合| 99国产极品粉嫩在线观看| 国产精品一区www在线观看| 九九在线视频观看精品| 中国国产av一级| 18+在线观看网站| 十八禁网站免费在线| 国产大屁股一区二区在线视频| 热99在线观看视频| 高清毛片免费观看视频网站| 大香蕉久久网| 美女xxoo啪啪120秒动态图| 最近最新中文字幕大全电影3| av中文乱码字幕在线| 国产麻豆成人av免费视频| 91久久精品电影网| 麻豆国产av国片精品| 亚洲欧美精品自产自拍| 俺也久久电影网| 一级黄色大片毛片| 成人无遮挡网站| .国产精品久久| 成人无遮挡网站| 国产精品亚洲一级av第二区| av视频在线观看入口| 日本熟妇午夜| 一夜夜www| 91狼人影院| 99九九线精品视频在线观看视频| 国产在线男女| 精品久久国产蜜桃| 亚洲五月天丁香| 在线观看66精品国产| 久久韩国三级中文字幕| 色综合色国产| 丰满乱子伦码专区| 成人精品一区二区免费| 色综合站精品国产| 国产在线男女| 久久久精品94久久精品| 午夜激情欧美在线| 国内精品美女久久久久久| 亚洲va在线va天堂va国产| 看免费成人av毛片| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 少妇的逼好多水| 天堂网av新在线| 亚洲欧美成人综合另类久久久 | 国产精品久久久久久久电影| 六月丁香七月| 精品不卡国产一区二区三区| 亚洲婷婷狠狠爱综合网| av福利片在线观看| 一级毛片我不卡| 欧美最新免费一区二区三区| 亚洲18禁久久av| 精品人妻偷拍中文字幕| 国产成人一区二区在线| or卡值多少钱| 精品久久久久久久人妻蜜臀av| 欧美成人免费av一区二区三区| 午夜福利在线观看免费完整高清在 | 六月丁香七月| 一个人看的www免费观看视频| a级毛片免费高清观看在线播放| 午夜福利成人在线免费观看| 啦啦啦观看免费观看视频高清| 久久精品久久久久久噜噜老黄 | 黄色日韩在线| 夜夜爽天天搞| 日韩国内少妇激情av| 色综合站精品国产| 亚洲成av人片在线播放无| 俄罗斯特黄特色一大片| 国产精品99久久久久久久久| 乱人视频在线观看| 国产精品不卡视频一区二区| 偷拍熟女少妇极品色| 国产亚洲精品综合一区在线观看| 亚洲无线在线观看| 蜜桃亚洲精品一区二区三区| 日韩成人伦理影院| 国产一区二区在线av高清观看| 一级毛片aaaaaa免费看小| 日日撸夜夜添| 晚上一个人看的免费电影| 国内揄拍国产精品人妻在线| 久久精品综合一区二区三区| 久久精品夜色国产| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| eeuss影院久久| 99热这里只有精品一区| 久久人人精品亚洲av| 人妻少妇偷人精品九色| 在线观看av片永久免费下载| 亚洲最大成人手机在线| 在线观看66精品国产| 毛片女人毛片| 成年女人看的毛片在线观看| 99在线视频只有这里精品首页| 身体一侧抽搐| 亚洲无线在线观看| 18+在线观看网站| 午夜亚洲福利在线播放| 亚洲欧美成人精品一区二区| 国产高清有码在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 一本一本综合久久| 国产精品久久久久久亚洲av鲁大| 日本黄色片子视频| 无遮挡黄片免费观看| 国产精品永久免费网站| 亚洲真实伦在线观看| 成人二区视频| 国产蜜桃级精品一区二区三区| 99热这里只有是精品50| 成年女人永久免费观看视频| 国产午夜精品论理片| 嫩草影院新地址| 婷婷亚洲欧美| 一级毛片我不卡| 91久久精品电影网| 麻豆av噜噜一区二区三区| 精品久久久久久久久亚洲| 国产又黄又爽又无遮挡在线| 欧美精品国产亚洲| 国产精品人妻久久久久久| 嫩草影院新地址| 欧美日韩一区二区视频在线观看视频在线 | 日本熟妇午夜| 欧美xxxx黑人xx丫x性爽| 亚洲中文字幕一区二区三区有码在线看| 别揉我奶头~嗯~啊~动态视频| 一区二区三区高清视频在线| 日韩 亚洲 欧美在线| 成年版毛片免费区| 男人和女人高潮做爰伦理| 人人妻人人澡欧美一区二区| 99热这里只有精品一区| 久久久久久久久久黄片| 中文字幕av在线有码专区| 69av精品久久久久久| 你懂的网址亚洲精品在线观看 | 亚洲18禁久久av| 97超碰精品成人国产| 中国国产av一级| 黄色日韩在线| 亚洲欧美中文字幕日韩二区| 麻豆久久精品国产亚洲av| 国产日本99.免费观看| 一a级毛片在线观看| 国产在视频线在精品| 深夜a级毛片| 精品一区二区三区视频在线| av天堂在线播放| 欧洲精品卡2卡3卡4卡5卡区| 天堂√8在线中文| 亚洲欧美成人精品一区二区| 丝袜美腿在线中文| 熟女电影av网| 少妇人妻精品综合一区二区 | 日本黄色视频三级网站网址| 久久久国产成人免费| 神马国产精品三级电影在线观看| 老司机福利观看| 少妇人妻一区二区三区视频| 国内久久婷婷六月综合欲色啪| 男人舔女人下体高潮全视频| 国产aⅴ精品一区二区三区波| 男女边吃奶边做爰视频| 97超碰精品成人国产| 亚洲高清免费不卡视频| 国产高清视频在线播放一区| 久久九九热精品免费| 99国产极品粉嫩在线观看| 欧美成人精品欧美一级黄| 国产精品一及| 久久久久久久亚洲中文字幕| 女生性感内裤真人,穿戴方法视频| 国内精品一区二区在线观看| 美女内射精品一级片tv| 欧美三级亚洲精品| 国产亚洲欧美98| 国产精品久久久久久精品电影| 精品一区二区三区av网在线观看| a级毛片a级免费在线| 免费看光身美女| 精品久久久久久久末码| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| 亚洲国产欧洲综合997久久,| 欧美人与善性xxx| 春色校园在线视频观看| 成人av一区二区三区在线看| 极品教师在线视频| 波多野结衣高清作品| 亚洲中文日韩欧美视频|