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

    溫壓彈爆炸沖擊波作用下防護門破壞模式研究

    2016-09-18 02:45:46李秀地耿振剛苗朝陽
    振動與沖擊 2016年16期
    關(guān)鍵詞:防護門坑道沖擊波

    李秀地, 耿振剛, 苗朝陽, 楊 森

    (1.后勤工程學院 土木工程系 巖土力學與地質(zhì)環(huán)境保護重慶市重點實驗室,重慶 401311;2. 92303部隊,山東 青島 266404)

    ?

    溫壓彈爆炸沖擊波作用下防護門破壞模式研究

    李秀地1, 耿振剛1, 苗朝陽1, 楊森2

    (1.后勤工程學院 土木工程系 巖土力學與地質(zhì)環(huán)境保護重慶市重點實驗室,重慶401311;2. 92303部隊,山東青島266404)

    為研究溫壓彈爆炸沖擊波作用下防護門的破壞模式,利用ANSYS/LS-DYNA軟件建立了考慮材料應(yīng)變率效應(yīng)的分離式鋼筋混凝土防護門數(shù)值計算模型。對防護門加載溫壓彈爆炸沖擊波波形,分析了防護門的撓度、轉(zhuǎn)角和剪應(yīng)力,并與TNT的作用結(jié)果進行了對比。結(jié)果表明:溫壓彈爆炸沖擊波作用下,隨著距爆炸源距離的增加,防護門的破壞模式由直剪破壞、彎剪耦合破壞逐漸變?yōu)閺澢茐?;防護門的響應(yīng)與溫壓彈爆炸沖擊波的高沖量有關(guān),溫壓彈爆炸沖擊波作用下防護門的破壞程度大于TNT的作用結(jié)果;防護門跨中速度值可以用來快速評估防護門的破壞模式。

    爆炸力學;溫壓彈;鋼筋混凝土;防護門;破壞模式

    近年來以美國為首的西方國家對發(fā)展溫壓武器的興趣越來越濃厚,開發(fā)了包括BLU-118B溫壓炸彈、AGM-114M溫壓導(dǎo)彈等多種溫壓彈。與傳統(tǒng)高爆彈藥相比,溫壓彈獨特的爆炸毀傷效應(yīng)主要體現(xiàn)在長持續(xù)時間高壓沖擊波和熱殺傷,以及因爆炸耗氧造成的人員窒息等傷害,特別適于打擊坑道等封閉空間內(nèi)的人員和設(shè)備。研究溫壓彈爆炸效應(yīng)及其防護技術(shù)對提高坑道工程的戰(zhàn)時生存能力具有重要意義[1]??拥纼?nèi)的防護門距離口部一般較遠,熱殺傷不是主要威脅,而爆炸沖擊波將是防護門破壞的決定性因素。因此,本文對溫壓彈爆炸沖擊波作用下防護門的破壞模式進行研究。

    溫壓彈爆炸沖擊波研究方面,茍兵旺等[2-3]進行了溫壓炸藥與TNT坑道內(nèi)爆炸試驗,實測了坑道內(nèi)爆炸沖擊波波形,并進行了對比。防護門破壞模式研究方面,閻石等[4]研究了等沖量的TNT爆炸沖擊波不同加載速率下鋼筋混凝土板的破壞模式,結(jié)果表明隨著沖擊波超壓峰值的增大,鋼筋混凝土板破壞模式由彎曲破壞逐漸變?yōu)榧羟衅茐?;汪維等[5]研究了不同裝藥量TNT近爆作用下鋼筋混凝土板的破壞,結(jié)果表明隨著裝藥量的增大,鋼筋混凝土板的破壞模式由整體彎曲破壞逐漸變?yōu)榫植繘_切破壞。溫壓彈爆炸沖擊波作用下防護門的破壞模式研究幾乎處于空白,本文采用ANSYS/LS-DYNA[6]有限元軟件建立鋼筋混凝土防護門的分離式共節(jié)點計算模型,在試驗實測數(shù)據(jù)驗證的基礎(chǔ)上,分析溫壓彈爆炸沖擊波作用下防護門的破壞模式,并與TNT的作用結(jié)果進行對比。

    1 溫壓彈坑道內(nèi)爆炸沖擊波

    基于坑道堵口爆炸試驗實測數(shù)據(jù)[7]驗證,利用ANSYS/LS-DYNA軟件建立了坑道內(nèi)爆炸沖擊波三維數(shù)值計算模型。分別計算了一種RDX基含鋁溫壓炸藥(RDX/AL/黏結(jié)劑=68/28/4,文中記為TBE)和TNT各75 kg堵口爆炸時,坑道內(nèi)爆炸沖擊波的傳播規(guī)律,得到不同距離處爆炸沖擊波波形,如圖1所示。

    圖1 沖擊波波形Fig.1 Blast wave in tunnel

    由圖1可知,坑道內(nèi)溫壓炸藥爆炸沖擊波與TNT爆炸沖擊波波形均為鋸齒狀,且隨著傳播距離的增大,超壓逐漸衰減,正相作用時間逐漸增加。但溫壓炸藥爆炸沖擊波的傳播速度更快、超壓更高、正相作用時間更長。這是由于溫壓炸藥中添加的鋁粉等產(chǎn)生后燃反應(yīng),增大了溫壓炸藥的爆炸沖量,提高了溫壓炸藥的爆炸威力。本文計算條件下的溫壓炸藥沖擊波入射超壓峰值約為TNT的1.23倍~1.50倍,沖量約為TNT的1.42倍~1.62倍,這與文獻[2-3]坑道內(nèi)實測比值基本相同。

    2 防護門數(shù)值分析模型的建立與驗證

    2.1防護門計算模型

    防護門選取鋼筋混凝土單扇門,長2 m、寬1.5 m、厚0.27 m,門框厚0.06 m,門框與門體接觸寬度為0.05 m。雙層雙向布置φ18HRB335橫向受拉鋼筋和φ10HRB335縱向分布鋼筋,布置φ8HRB335單肢箍筋,受拉及分布鋼筋間距0.09 m,箍筋間距0.18 m。采用分離式共節(jié)點方式建模,由于對稱性,建立四分之一模型。網(wǎng)格尺寸取0.015 m,防護門的有限元模型如圖2所示。

    圖2 防護門有限元模型Fig.2 Finite model of blast door

    2.2材料模型

    混凝土選用MAT_CONCRETE_DAMAGE_ REL3材料。該模型考慮了混凝土材料的應(yīng)變率效應(yīng),僅需要輸入密度、泊松比和單軸抗壓強度三個參數(shù),便能較好的模擬爆炸沖擊作用下混凝土的動態(tài)響應(yīng)問題[8]。選用LS-DYNA中單軸抗壓強度為40 MPa的混凝土,混凝土密度取2.4×103kg/m3,泊松比取0.2。鋼筋使用塑性隨動模型MAT_PLASTIC_KINEMATIC,門框選擇剛體模型MAT_RIGID。混凝土的抗拉及抗壓強度動載增大系數(shù)DIF均采用歐洲混凝土規(guī)范CEB[9]推薦的公式計算,鋼筋的動載增大系數(shù)DIF采用Malvar[10]提出的公式計算。鋼筋材料參數(shù)如表1所示。

    表1 鋼筋材料參數(shù)

    2.3邊界條件

    防護門邊界條件通??紤]簡支或設(shè)置接觸。楊心宇等[11]通過對防護門在不同約束條件下的動力響應(yīng)進行分析,表明防護門配筋按四邊簡支計算偏于保守;陸新征等[12]通過設(shè)置門扇與門框、門軸與軸瓦兩種接觸研究了抗爆門的動力響應(yīng),并得到了抗爆門的反彈力。本文在考慮防護門的實際支承情況下設(shè)置了簡化接觸,即門扇與門框間設(shè)置Automatic Contact、Surface to Surface 接觸類型。

    2.4爆炸荷載的施加

    LS-DYNA可以通過流固耦合算法、CONWEP算法和直接加載法3種方法對防護門施加爆炸沖擊波荷載。CONWEP算法可以方便地計算自由場TNT爆炸沖擊波荷載,但無法計算坑道內(nèi)溫壓炸藥的爆炸沖擊波。流固耦合算法通常用于計算近爆作用下結(jié)構(gòu)動力響應(yīng)[5,13],防護門距坑道口部通常達到10 m以上,運用流固耦合算法模型單元數(shù)量十分驚人,計算甚至難以進行。直接加載法將沖擊波曲線作為面荷載直接施加于防護門上,模型簡單、概念清晰、計算效率高,已被廣泛應(yīng)用于研究鋼筋混凝土構(gòu)件在爆炸沖擊波作用下的動力響應(yīng)問題[12],本文運用直接加載法研究溫壓彈沖擊波作用下防護門的破壞模式。將圖1所示溫壓彈和TNT坑道內(nèi)入射沖擊波通過下式轉(zhuǎn)換為反射沖擊波[14]:

    (1)

    式中:ΔPm、ΔPrm分別為入射沖擊波和正反射沖擊波超壓峰值,MPa;n為正反射系數(shù),按表2取值;當ΔPm為表中范圍內(nèi)其它值時,用內(nèi)插法確定n值。

    表2 空氣沖擊波正反射系數(shù)[14]

    計算中將通過式(1)計算所得沖擊波正反射波形直接加載于防護門上,設(shè)置工況如表3所示。

    表3 計算工況

    2.5計算模型的驗證

    文獻[15]對鋼筋混凝土板在沖擊波作用下的動力響應(yīng)進行了研究,為驗證本文上述數(shù)值模型的可靠性,對文獻[15]中實驗NSC-NR進行模擬。采用直接加載法,按照實驗實際邊界條件,跨中位移時程曲線實驗及數(shù)值模擬計算結(jié)果,如圖3所示。

    圖3 計算模型的驗證Fig.3 Model validation

    由圖3可知,數(shù)值計算所得曲線與試驗實測曲線取得了較好的一致性;數(shù)值計算跨中撓度最大值為212.7 mm,與試驗實測值221 mm誤差為3.76%,表明本文提出的數(shù)值計算模型可以較好地模擬鋼筋混凝土構(gòu)件在爆炸沖擊波作用下的動力響應(yīng)。

    3 防護門破壞模式分析

    鋼筋混凝土防護門在爆炸荷載作用下的破壞模式通常有直剪破壞、彎曲破壞、彎剪耦合破壞等整體破壞和混凝土壓碎、開裂、破片等局部破壞[15],本文主要從整體破壞的角度分析防護門的破壞模式。以工況1、5和9作為典型荷載,研究溫壓彈堵口爆炸作用下,防護門在近、中和遠3個位置的破壞模式。

    3.1彎曲破壞分析

    以防護門支座處轉(zhuǎn)角作為其彎曲破壞特征進行研究,根據(jù)文獻[17],配單肢箍筋時,當支座轉(zhuǎn)角達到4°,鋼筋混凝土受彎構(gòu)件喪失抗彎承載能力;當支座轉(zhuǎn)角到達12°,構(gòu)件將因喪失整體性而完全破壞。本文計算以支座處轉(zhuǎn)角4°和12°為界,當支座處轉(zhuǎn)角小于4°時,防護門發(fā)生輕度彎曲變形,視為未破壞;當支座處轉(zhuǎn)角大于4°時,防護門發(fā)生嚴重彎曲破壞;當支座處轉(zhuǎn)角小于12°時,防護門雖遭受破壞,但抗住了一次打擊作用。3種工況下防護門背爆面跨中位移時程曲線如圖4所示。

    圖4 跨中位移時程曲線Fig.4 History curve of mid span displacement

    由圖4可知,防護門跨中撓度先增大后減小,最后殘余變形保持穩(wěn)定。溫壓彈堵口爆炸時防護門距口部5 m、15 m和25 m情況下,跨中最大撓度分別為72.3 cm、18.8 cm、8.7 cm。防護門支座處轉(zhuǎn)角θ可由下式計算:

    (2)

    式中:X為防護門跨中撓度;B為防護門寬度。

    將工況1、5和9計算得到的防護門跨中最大撓度分別代入式(2),可得防護門轉(zhuǎn)角分別為43.9°、14.8°和6.6°。因此,工況1和5條件下,防護門均因喪失整體性而完全破壞;工況9條件下,防護門雖未喪失整體性,但已發(fā)生嚴重彎曲破壞。

    3.2剪切破壞分析

    抗壓強度為40 MPa的混凝土,其抗剪強度為[16]:

    利用LS-PREPOST后處理軟件,以此抗剪強度為標準,可以清晰地觀察防護門受剪失效單元的位置。選取爆炸沖擊波作用初期(t=1 ms)和防護門達到最大撓度前典型剪切破壞的時刻,不同工況下防護門的剪切破壞形態(tài)如圖5所示。圖中淺色區(qū)域為受剪失效單元。

    圖5 防護門剪切破壞形態(tài)Fig.5 Shear failure mode of blast door

    由圖5(a)可知,對于工況1,在沖擊波作用初期(t=1 ms),防護門尚未達到最大撓度,支座處已產(chǎn)生較寬的貫通剪切裂縫,表明防護門發(fā)生了直接剪切破壞。達到最大撓度前,防護門支座處剪切裂縫進一步發(fā)展(t=11 ms),防護門門體沿支座處剪切裂縫發(fā)展方向,被沖擊波“推出”門框外。

    由圖5(b)可知,對于工況5,在沖擊波作用初期(t=1 ms),防護門支座處有較寬的剪切裂縫,但沒有貫通。隨著沖擊波作用時間的增加,防護門跨中剪跨區(qū)出現(xiàn)多條斜剪裂縫(t=13 ms),部分已經(jīng)貫通,有明顯的彎剪耦合破壞特征。

    由圖5(c)可知,對于工況9,在沖擊波作用初期(t=1 ms),防護門僅有少量支座附近單元發(fā)生剪切失效。隨著沖擊波作用時間的增加,防護門多個位置有剪切裂縫(t=15 ms),但裂縫寬度較窄,尚未完全貫通。

    由以上分析可知,工況1防護門主要破壞模式為直剪破壞,防護門被完全摧毀。工況5防護門主要破壞模式為彎剪耦合破壞,防護門破壞嚴重,工況9防護門主要破壞模式為彎曲破壞。

    4 溫壓彈與TNT作用下防護門破壞模式對比

    按照上述破壞模式分析方法,計算表3所示各工況下,防護門跨中撓度、轉(zhuǎn)角、達到最大位移時間Tmax、跨中最大速度Vmax和主要破壞模式分別如表4所示。

    表4 防護門計算結(jié)果

    由表4可知,隨著距爆炸源距離的增加,本文計算溫壓彈和TNT爆炸沖擊波作用下防護門的撓度、轉(zhuǎn)角和跨中最大速度均不斷地減小。相同位置處,溫壓彈爆炸沖擊波作用下防護門的撓度、轉(zhuǎn)角和跨中最大速度大于TNT作用的情況,防護門的破壞模式不同于TNT作用的情況。距坑道口部大于15 m時,防護門可以有效防護TNT爆炸沖擊波作用;但在溫壓彈爆炸沖擊波作用下,防護門仍然嚴重破壞。表明溫壓彈具有更高的爆炸威力。

    以工況4與工況7為例,分析溫壓彈爆炸沖擊波超壓峰值低于TNT情況下防護門的毀傷。工況4 TNT爆炸沖擊波超壓為工況7溫壓彈爆炸沖擊波超壓的1.43倍,但溫壓彈爆炸沖擊波作用下防護門最大位移卻大于TNT的作用。這是由于溫壓彈沖擊波沖量達TNT沖擊波沖量的1.62倍的結(jié)果,表明溫壓彈的高毀傷效能與其高沖量有關(guān)。

    由表4并結(jié)合圖4可知,不同破壞模式時,防護門跨中位移時程曲線初始段斜率(即跨中速度值)不同。大量計算表明,跨中速度值可以用來快速判斷防護門的破壞模式。本文計算條件下,跨中速度值大于18.4 m/s時,防護門發(fā)生直剪破壞;跨中速度值為10.0~18.4 m/s時,防護門發(fā)生彎剪耦合破壞;跨中速度值為3.8~10.0 m/s時,防護門發(fā)生彎曲破壞;跨中速度值小于3.8 m/s時,防護門不發(fā)生破壞。

    5 結(jié) 論

    基于試驗數(shù)據(jù)驗證,建立了防護門的數(shù)值計算模型,研究了溫壓彈爆炸沖擊波作用下防護門的破壞模式,并與TNT的作用結(jié)果進行了對比,得出以下結(jié)論:

    (1) 溫壓彈爆炸沖擊波作用下,隨著防護門距爆炸源距離的增大,防護門的破壞模式由近爆處的直剪破壞變?yōu)橹斜幍膹澕赳詈掀茐暮瓦h爆處的彎曲破壞。

    (2) 防護門的響應(yīng)與溫壓彈爆炸沖擊波的高沖量有關(guān),溫壓彈爆炸沖擊波作用下坑道內(nèi)相同位置處防護門的破壞程度大于TNT作用的情況,防護門的破壞模式與TNT作用時不盡相同。

    (3) 防護門的破壞模式可以通過跨中速度值快速評估。本文計算條件下,跨中速度值大于18.4 m/s時,防護門發(fā)生直剪破壞;跨中速度值為10.0~18.4 m/s時,防護門發(fā)生彎剪耦合破壞;跨中速度值為3.8~10.0 m/s時,防護門發(fā)生彎曲破壞;跨中速度小于3.8 m/s時,防護門不發(fā)生破壞。

    [1] 周豐峻. 信息化戰(zhàn)爭條件下作戰(zhàn)工程保障和工程兵多樣化軍事行動任務(wù)探索[C]//中國土木工程學會防護工程分會第十二次學術(shù)年會論文集.張家界,2010:1-6.

    [2] 茍兵旺,李芝榮,閆瀟敏,等. 復(fù)雜坑道內(nèi)溫壓炸藥沖擊波實驗研究[J]. 火工品,2014(2):41-45.

    GOU Bingwang, LI Zhirong, YAN Xiaomin, et al. Experimental study on shock wave effects of thermo-baric explosive in complex tunnel[J]. Initiators and Pyrotechnics, 2014(2):41-45.

    [3] 孔霖,蘇健軍,李芝榮,等. 不同裝藥坑道內(nèi)爆炸沖擊波傳播規(guī)律的試驗研究[J]. 火工品,2012(3):21-24.

    KONG Lin, SU Jianjun, LI Zhirong, et al. Test study on explosion shock wave propagation of different explosives inside tunnels[J]. Initiators and Pyrotechnics, 2012(3):21-24.

    [4] 閻石,張亮,王丹. 鋼筋混凝土板在爆炸荷載作用下的破壞模式分析[J]. 沈陽建筑大學學報:自然科學版,2005, 21(3):177-180.

    YAN Shi, ZHANG Liang, WANG Dan. Failure mode analysis for RC slab under explosive loads[J]. Journal of Shenyang Architecture and Civil Engineering University:Nature Science Edition, 2005, 21(3):177-180.

    [5] 汪維,張舵,盧芳云,等. 方形鋼筋混凝土板的近場抗爆性能[J]. 爆炸與沖擊,2012,32(3): 251-258.

    WANG Wei, ZHANG Duo, LU Fangyun, et al. Anti-explosion performances of square reinforced concrete slabs under close-in explosions[J]. Explosion and Shock Waves,2012,32(3):251-258.

    [6] LS-DYNA keyword uer’s manual[M]. Livemore, California: Livemore Software Technology Corporation, 2007, 5, Version 971.

    [7] 吳云泉, 許金余, 陳海天,等. 坑道口部爆炸沖擊波壓力傳播規(guī)律[J]. 空軍工程大學學報:自然科學版,2009,10(5):90-94.

    WU Yunquan, XU Jinyu, CHEN Haitian, et al. Study of shock wave pressure propagation laws of the explosion at the tunnel entrance[J]. Journal of Air Force Engineering University:Natural Science Edition, 2009, 10(5):90-94.

    [8] TU Z G, LU Y. Evaluation of typical concrete material models used in hydrocodes for high dynamic response simulations[J]. International Journal of Impact Engineering, 2009,36:132-146.

    [9] DU B C E. CEB-FIP model code 1990: Design code[M]. London: Thomas Telford Ltd., 1990.

    [10] MALVAR L J.Review of static and dynamic properties of steel reinforcing bars[J].ACI Materials Journal, 1998, 95(5):609-614.

    [11] 楊心宇,趙躍堂,周澤鑫. 平板防護門不同約束形式下的動力響應(yīng)[J]. 爆破, 2014, 2(31):47-56.

    YANG Xinyu, ZHAO Yuetang, ZHOU Zexin. Dynamic response of blast door under different forms of constrains[J]. Blasting, 2014, 2(31):47-56.

    [12] 陸新征,江見鯨. 抗爆門結(jié)構(gòu)考慮接觸影響的動力有限元分析[J]. 力學與實踐, 2003,25(2):24-26.

    LU Xinzheng, JIANG Jianjing. Safety assessment of blast-resistant doors using dynamic finite element method and contact analysis[J]. Mechanics in Engineering, 2003,25(2):24-26.

    [13] 龔順風,鄧歡,朱升波,等. 近爆作用下鋼筋混凝土板動態(tài)破壞的數(shù)值模擬研究[J]. 振動與沖擊,2012,31(2):20-24.

    GONG Shunfeng, DENG Huan, ZHU Shengbo, et al. Numerical simulation for dynamic failure of a reinforced concrete slab under close-in explosion[J]. Journal of Vibration and Shock, 2012,31(2):20-24.

    [14] 方秦,柳錦春. 地下防護結(jié)構(gòu)[M]. 北京:中國水利水電出版社,2010.

    [15] THIAGARAJAN G, KADAMBI A V, ROBERT S, et al. Experimental and finite element analysis of doubly reinforced concrete slabs subjected to blast loads[J]. International Journal of Impact Engineering, 2015,75:162-173.

    [16] TRIVEDI N, SINGH R K. Prediction of impact induced failure modes in reinforced concrete slabs through nonlinear transient dynamic finite element simulation[J]. Annals of Nuclear Energy, 2013,56:109-121.

    [17] TM5-1300. The design of structures to resist the effects of accidental explosions[M].Washington DC:US Department of the Army, Navy, and Air Force, 1990.

    Failure mode of a blast door subjected to the explosion wave of a thermobaric bomb

    LI Xiudi1, GENG Zhengang1, MIAO Chaoyang1, YANG Sen2

    (1. Chongqing Key Laboratory of Geomechanics & Geoenvironment Protection in Department of Civil Engineering,Logistical Engineering University, Chongqing 401311, China;2. Unit 92303, Qingdao 266404, China)

    To reveal the failure mode of a blast door subjected to the explosion wave of a thermobaric bomb, a separate reinforced concrete numerical model considering strain rate effects of a blast door was established with the dynamic finite element software ANSYS/LS-DYNA. The explosion wave of a thermobaric bomb was loaded on the blast door and the deflection and shear stress of the blast door were analyzed and compared with the shock wave of TNT. The results show that with the distance from explosion increasing, the failure mode of blast door under the thermobaric bomb gradually changes from punching shear failure and mixed flexural-shear failure to flexural failure. The dynamic response of the blast door is related to the high impulse of the thermobaric bomb, and the damage of the blast door under the thermobaric bomb is more serious than that under TNT. The velocity of the mid span can be used to estimate the failure mode of the blast door.

    mechanics of explosion; thermobaric bomb; reinforced concrete; blast door; failure mode

    國家部委基金(CY213J009)

    2015-09-29修改稿收到日期:2015-12-31

    李秀地 男,博士,教授,1970年生

    O383

    A

    10.13465/j.cnki.jvs.2016.16.032

    猜你喜歡
    防護門坑道沖擊波
    點陣夾芯結(jié)構(gòu)拱形防護門抗爆性能研究
    水下明挖隧道防護門設(shè)置對疏散救援影響研究
    北方交通(2022年8期)2022-08-15 09:48:12
    一種智能數(shù)控車削單元防護門控制系統(tǒng)的設(shè)計
    武漢沖擊波
    中國公路(2019年10期)2019-06-28 03:05:08
    戰(zhàn)時指揮所
    能源物聯(lián)網(wǎng)沖擊波
    能源(2018年10期)2018-12-08 08:02:34
    一種火箭起飛后快速關(guān)閉的液氫防護塔
    醫(yī)生集團沖擊波
    超聲雙探頭聯(lián)合定位法在體外沖擊波碎石術(shù)中的應(yīng)用
    半干旱區(qū)石質(zhì)坡面坑道式造林整地小氣候效應(yīng)觀測初報
    防護林科技(2013年1期)2013-08-08 11:20:18
    欧美日韩中文字幕国产精品一区二区三区| 婷婷丁香在线五月| 久久久国产成人免费| 他把我摸到了高潮在线观看| 身体一侧抽搐| 亚洲欧美激情综合另类| 香蕉国产在线看| 国产不卡一卡二| 91字幕亚洲| 九九在线视频观看精品| 97超级碰碰碰精品色视频在线观看| 最近最新免费中文字幕在线| 99久久精品热视频| 亚洲 国产 在线| 亚洲熟女毛片儿| 色老头精品视频在线观看| 人妻丰满熟妇av一区二区三区| 一边摸一边抽搐一进一小说| 欧美激情在线99| 极品教师在线免费播放| 国产精品一区二区三区四区久久| 久久国产乱子伦精品免费另类| 国产精品一区二区精品视频观看| 日韩欧美在线乱码| 日韩欧美国产在线观看| 99久国产av精品| 在线免费观看不下载黄p国产 | 亚洲激情在线av| 村上凉子中文字幕在线| 老熟妇仑乱视频hdxx| 白带黄色成豆腐渣| 欧美日韩国产亚洲二区| 高潮久久久久久久久久久不卡| av片东京热男人的天堂| 非洲黑人性xxxx精品又粗又长| 在线观看午夜福利视频| 亚洲aⅴ乱码一区二区在线播放| 18禁国产床啪视频网站| 亚洲av成人不卡在线观看播放网| 久久香蕉精品热| 国产亚洲av高清不卡| 小蜜桃在线观看免费完整版高清| 中文在线观看免费www的网站| 精品午夜福利视频在线观看一区| 男女视频在线观看网站免费| 精品久久久久久成人av| 天天躁日日操中文字幕| 国产三级黄色录像| 国产精品爽爽va在线观看网站| 九色成人免费人妻av| 国产精品免费一区二区三区在线| 村上凉子中文字幕在线| 男人舔女人的私密视频| 一进一出抽搐gif免费好疼| 亚洲中文字幕日韩| 曰老女人黄片| 欧美日韩瑟瑟在线播放| 欧美绝顶高潮抽搐喷水| 天堂av国产一区二区熟女人妻| 成人一区二区视频在线观看| www.精华液| 亚洲男人的天堂狠狠| 精品99又大又爽又粗少妇毛片 | 午夜免费激情av| 久久久国产成人免费| 好男人在线观看高清免费视频| 国产成人影院久久av| 一区二区三区国产精品乱码| 午夜久久久久精精品| 99热这里只有精品一区 | 国产伦精品一区二区三区四那| 高清在线国产一区| 熟女少妇亚洲综合色aaa.| 日韩有码中文字幕| 国产精品美女特级片免费视频播放器 | 脱女人内裤的视频| 国产又黄又爽又无遮挡在线| 日韩欧美国产在线观看| 国产 一区 欧美 日韩| 国产又黄又爽又无遮挡在线| 中文在线观看免费www的网站| 黑人欧美特级aaaaaa片| 久久久久精品国产欧美久久久| 少妇的丰满在线观看| 欧美日韩亚洲国产一区二区在线观看| 黄色片一级片一级黄色片| 久久九九热精品免费| 午夜精品久久久久久毛片777| 国产真实乱freesex| 日韩有码中文字幕| 精品熟女少妇八av免费久了| 别揉我奶头~嗯~啊~动态视频| 精品欧美国产一区二区三| 久9热在线精品视频| x7x7x7水蜜桃| 午夜福利成人在线免费观看| 国产久久久一区二区三区| 国产精品永久免费网站| 日韩大尺度精品在线看网址| 亚洲 欧美一区二区三区| 国产在线精品亚洲第一网站| 国产精品野战在线观看| 一区二区三区国产精品乱码| 亚洲一区二区三区不卡视频| 特级一级黄色大片| 久久精品亚洲精品国产色婷小说| 免费看光身美女| www日本在线高清视频| 无限看片的www在线观看| 成人欧美大片| 九色国产91popny在线| 一级毛片高清免费大全| 国产成+人综合+亚洲专区| 国产亚洲精品一区二区www| 国产精品亚洲美女久久久| 日韩国内少妇激情av| 日韩精品中文字幕看吧| 国产精品99久久久久久久久| 国产乱人伦免费视频| 国产午夜福利久久久久久| 国产精品久久久久久久电影 | 波多野结衣巨乳人妻| 99久久精品一区二区三区| 国语自产精品视频在线第100页| 国产精品99久久久久久久久| 超碰成人久久| 级片在线观看| 毛片女人毛片| 老汉色av国产亚洲站长工具| 午夜福利高清视频| 一个人看视频在线观看www免费 | 我要搜黄色片| 好男人电影高清在线观看| 一级毛片女人18水好多| bbb黄色大片| 俄罗斯特黄特色一大片| 999久久久国产精品视频| 国产伦在线观看视频一区| 成人亚洲精品av一区二区| a级毛片在线看网站| 在线免费观看不下载黄p国产 | 精品久久蜜臀av无| 亚洲天堂国产精品一区在线| 人妻丰满熟妇av一区二区三区| 国产欧美日韩精品一区二区| 国产1区2区3区精品| 一区二区三区国产精品乱码| 岛国在线观看网站| 久久久久亚洲av毛片大全| 午夜福利18| 亚洲av成人不卡在线观看播放网| 国产野战对白在线观看| 91av网一区二区| 亚洲av免费在线观看| 久久精品影院6| 欧美日韩乱码在线| 性色avwww在线观看| 三级毛片av免费| 欧美极品一区二区三区四区| 亚洲avbb在线观看| 精品国内亚洲2022精品成人| 在线国产一区二区在线| 免费在线观看视频国产中文字幕亚洲| 亚洲精品一卡2卡三卡4卡5卡| 老汉色∧v一级毛片| 在线视频色国产色| 久久99热这里只有精品18| 国产成年人精品一区二区| 午夜福利免费观看在线| 国产成人精品久久二区二区91| 免费看十八禁软件| 久久性视频一级片| 欧美日韩一级在线毛片| 欧美成人免费av一区二区三区| 国内精品久久久久精免费| 欧美乱色亚洲激情| av福利片在线观看| 午夜福利高清视频| 色噜噜av男人的天堂激情| 国产伦精品一区二区三区视频9 | 19禁男女啪啪无遮挡网站| 最新中文字幕久久久久 | 日韩中文字幕欧美一区二区| 无限看片的www在线观看| 啦啦啦韩国在线观看视频| 亚洲男人的天堂狠狠| 美女 人体艺术 gogo| 我要搜黄色片| 九色成人免费人妻av| 在线视频色国产色| 12—13女人毛片做爰片一| 国产毛片a区久久久久| 久久人妻av系列| 熟女少妇亚洲综合色aaa.| 免费在线观看亚洲国产| 毛片女人毛片| 这个男人来自地球电影免费观看| 19禁男女啪啪无遮挡网站| 午夜免费观看网址| 亚洲片人在线观看| 嫩草影院精品99| 母亲3免费完整高清在线观看| 后天国语完整版免费观看| 熟女少妇亚洲综合色aaa.| 9191精品国产免费久久| 亚洲熟女毛片儿| 午夜免费成人在线视频| 久久国产乱子伦精品免费另类| 亚洲18禁久久av| 久久久久久久精品吃奶| 久久久色成人| www.精华液| 搡老熟女国产l中国老女人| 国内精品美女久久久久久| 啦啦啦免费观看视频1| 欧美性猛交╳xxx乱大交人| www.www免费av| 国产单亲对白刺激| 精品久久久久久久久久久久久| 亚洲国产日韩欧美精品在线观看 | 欧美不卡视频在线免费观看| 久久久久久九九精品二区国产| 久久精品国产清高在天天线| 91av网站免费观看| 女警被强在线播放| 国产综合懂色| 国产高清videossex| 亚洲成av人片在线播放无| 久久久久久九九精品二区国产| 欧美黑人欧美精品刺激| 精品欧美国产一区二区三| 欧美最黄视频在线播放免费| 国产极品精品免费视频能看的| 男女床上黄色一级片免费看| 后天国语完整版免费观看| 精品国产乱子伦一区二区三区| 悠悠久久av| 国产成人福利小说| 国产美女午夜福利| 欧美高清成人免费视频www| 婷婷精品国产亚洲av在线| av在线天堂中文字幕| 美女高潮的动态| 岛国在线观看网站| 久久婷婷人人爽人人干人人爱| 欧美黑人巨大hd| 国产伦精品一区二区三区四那| 丁香欧美五月| 久久久久性生活片| 国内少妇人妻偷人精品xxx网站 | 日韩欧美国产一区二区入口| 美女cb高潮喷水在线观看 | 黄片大片在线免费观看| 免费无遮挡裸体视频| 九九热线精品视视频播放| 亚洲天堂国产精品一区在线| 中文字幕人成人乱码亚洲影| 黄片大片在线免费观看| 日韩国内少妇激情av| 国产午夜精品久久久久久| 亚洲第一电影网av| 在线观看一区二区三区| www国产在线视频色| 一区二区三区高清视频在线| 久久九九热精品免费| 国产精品一区二区三区四区免费观看 | 亚洲精品美女久久av网站| 蜜桃久久精品国产亚洲av| 精品电影一区二区在线| 99riav亚洲国产免费| 看片在线看免费视频| 999久久久国产精品视频| 色综合站精品国产| 色综合欧美亚洲国产小说| 又黄又爽又免费观看的视频| 97超视频在线观看视频| 免费在线观看亚洲国产| 老司机在亚洲福利影院| 91麻豆精品激情在线观看国产| 十八禁人妻一区二区| 99热这里只有是精品50| 国产精品影院久久| 最近最新中文字幕大全免费视频| 最近最新中文字幕大全电影3| 亚洲电影在线观看av| 国产精品一区二区三区四区久久| e午夜精品久久久久久久| 不卡av一区二区三区| 草草在线视频免费看| 国产亚洲精品久久久com| 一级毛片精品| 午夜福利18| 69av精品久久久久久| 国产 一区 欧美 日韩| 国产精品香港三级国产av潘金莲| 精品午夜福利视频在线观看一区| 一级a爱片免费观看的视频| 美女 人体艺术 gogo| 麻豆久久精品国产亚洲av| 国产精品久久久久久人妻精品电影| 久久久久久久久中文| 久久精品影院6| 一本精品99久久精品77| 久久亚洲精品不卡| 午夜成年电影在线免费观看| 国内揄拍国产精品人妻在线| 亚洲欧美精品综合久久99| 男女午夜视频在线观看| 夜夜夜夜夜久久久久| 亚洲 国产 在线| 久久久国产成人免费| 夜夜看夜夜爽夜夜摸| 亚洲国产欧美一区二区综合| 成人国产综合亚洲| 欧美在线一区亚洲| 成人国产综合亚洲| 99久久无色码亚洲精品果冻| 国语自产精品视频在线第100页| 亚洲精品中文字幕一二三四区| 亚洲国产日韩欧美精品在线观看 | 免费在线观看成人毛片| 婷婷精品国产亚洲av在线| 可以在线观看的亚洲视频| 亚洲一区二区三区色噜噜| 在线观看午夜福利视频| 国产精品99久久99久久久不卡| 欧美成人性av电影在线观看| 色老头精品视频在线观看| 黄频高清免费视频| 欧美日韩瑟瑟在线播放| 欧美在线一区亚洲| 国产精品1区2区在线观看.| 搡老妇女老女人老熟妇| 亚洲一区二区三区不卡视频| 日本a在线网址| 色综合婷婷激情| 亚洲国产欧美人成| 久久久久久久久中文| 女生性感内裤真人,穿戴方法视频| 国产成人av激情在线播放| 在线免费观看的www视频| 在线播放国产精品三级| 国产精品,欧美在线| 国产精品永久免费网站| 国产精品一及| 国产在线精品亚洲第一网站| 免费高清视频大片| 99精品在免费线老司机午夜| 九色成人免费人妻av| 高潮久久久久久久久久久不卡| 黄片小视频在线播放| 级片在线观看| 国产精品免费一区二区三区在线| 制服人妻中文乱码| 丰满人妻一区二区三区视频av | 看黄色毛片网站| 亚洲国产欧美人成| 亚洲午夜理论影院| av在线天堂中文字幕| netflix在线观看网站| 国产精品九九99| 精品一区二区三区av网在线观看| 长腿黑丝高跟| 小蜜桃在线观看免费完整版高清| 亚洲性夜色夜夜综合| netflix在线观看网站| 露出奶头的视频| 成人欧美大片| 精品久久久久久久毛片微露脸| 久久久久久久久中文| 国产精品爽爽va在线观看网站| 中文字幕av在线有码专区| 琪琪午夜伦伦电影理论片6080| 国产1区2区3区精品| 国产午夜精品久久久久久| 丁香六月欧美| 久久久久久大精品| 香蕉久久夜色| 日日干狠狠操夜夜爽| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 亚洲专区国产一区二区| 男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 免费观看的影片在线观看| 999久久久精品免费观看国产| 男人和女人高潮做爰伦理| 亚洲 欧美 日韩 在线 免费| 成人av在线播放网站| 男人舔女人的私密视频| 国产97色在线日韩免费| 国产一区在线观看成人免费| 国产精品日韩av在线免费观看| 国产1区2区3区精品| a级毛片a级免费在线| 国产亚洲精品av在线| 精品国产乱码久久久久久男人| 搞女人的毛片| 亚洲国产欧美一区二区综合| 精品国产亚洲在线| 精品熟女少妇八av免费久了| 波多野结衣高清无吗| 成年免费大片在线观看| 最新中文字幕久久久久 | 91av网一区二区| 欧美极品一区二区三区四区| 午夜免费观看网址| 免费观看精品视频网站| 中文字幕熟女人妻在线| 欧美日韩瑟瑟在线播放| 丝袜人妻中文字幕| 999精品在线视频| 亚洲人成电影免费在线| 村上凉子中文字幕在线| 国产精品久久电影中文字幕| 久久久久久久久久黄片| 看片在线看免费视频| 天天添夜夜摸| 欧美不卡视频在线免费观看| 桃红色精品国产亚洲av| ponron亚洲| 欧美色欧美亚洲另类二区| 丁香欧美五月| 国产人伦9x9x在线观看| www.999成人在线观看| 国产美女午夜福利| 最好的美女福利视频网| 国产激情久久老熟女| 欧美午夜高清在线| 国产欧美日韩精品一区二区| 精品国产乱码久久久久久男人| 日韩精品中文字幕看吧| 美女cb高潮喷水在线观看 | 婷婷丁香在线五月| 欧美午夜高清在线| 免费一级毛片在线播放高清视频| 91av网站免费观看| av福利片在线观看| 午夜精品久久久久久毛片777| 国产精品影院久久| 99热这里只有精品一区 | www.精华液| 久久久久久久久久黄片| 在线国产一区二区在线| 欧美在线一区亚洲| 不卡一级毛片| 国内精品久久久久久久电影| 一级a爱片免费观看的视频| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 一区二区三区国产精品乱码| 成人午夜高清在线视频| 精品国产超薄肉色丝袜足j| 国产成人av教育| 国产精品久久久av美女十八| 国产综合懂色| 亚洲av电影在线进入| 欧美绝顶高潮抽搐喷水| 国产成人影院久久av| 日本免费一区二区三区高清不卡| 免费看日本二区| 亚洲美女视频黄频| 精品国产乱子伦一区二区三区| 成年女人永久免费观看视频| 老司机深夜福利视频在线观看| 精品久久久久久久末码| 国内少妇人妻偷人精品xxx网站 | 亚洲18禁久久av| 9191精品国产免费久久| 日本 欧美在线| 我的老师免费观看完整版| xxxwww97欧美| 国产精品久久久久久精品电影| 看免费av毛片| 久久久精品欧美日韩精品| 国产精品一区二区三区四区免费观看 | 亚洲国产色片| 国产野战对白在线观看| 级片在线观看| 一个人免费在线观看的高清视频| 757午夜福利合集在线观看| 欧美成人一区二区免费高清观看 | 国产v大片淫在线免费观看| 亚洲激情在线av| 亚洲av中文字字幕乱码综合| 亚洲av熟女| 女同久久另类99精品国产91| 在线免费观看不下载黄p国产 | 啦啦啦韩国在线观看视频| 日韩欧美免费精品| 熟女电影av网| 两人在一起打扑克的视频| 黄色女人牲交| 日日摸夜夜添夜夜添小说| 国产精品精品国产色婷婷| 亚洲欧美日韩东京热| 久久性视频一级片| 久久人人精品亚洲av| 天堂影院成人在线观看| 亚洲狠狠婷婷综合久久图片| 三级男女做爰猛烈吃奶摸视频| 久99久视频精品免费| 国产一区在线观看成人免费| 日韩大尺度精品在线看网址| 久久中文字幕一级| 成人特级黄色片久久久久久久| 免费在线观看成人毛片| 国产高清激情床上av| 一级毛片精品| 免费人成视频x8x8入口观看| 久久久久久人人人人人| 熟女少妇亚洲综合色aaa.| 一二三四社区在线视频社区8| 又黄又粗又硬又大视频| 亚洲成av人片免费观看| 国产精品久久久久久亚洲av鲁大| 亚洲国产中文字幕在线视频| 19禁男女啪啪无遮挡网站| 老司机深夜福利视频在线观看| 亚洲精华国产精华精| 身体一侧抽搐| 老司机福利观看| 哪里可以看免费的av片| 日本熟妇午夜| 一级作爱视频免费观看| 男女下面进入的视频免费午夜| 亚洲 欧美一区二区三区| 亚洲av熟女| 国产高潮美女av| 成熟少妇高潮喷水视频| 国产一区二区在线av高清观看| 国产精品99久久久久久久久| 成年人黄色毛片网站| 午夜福利18| 亚洲人成网站高清观看| 午夜亚洲福利在线播放| 999久久久精品免费观看国产| 夜夜夜夜夜久久久久| 国产伦精品一区二区三区视频9 | 欧美日韩瑟瑟在线播放| 成人av在线播放网站| avwww免费| 99精品在免费线老司机午夜| 日日干狠狠操夜夜爽| 99久久久亚洲精品蜜臀av| 欧美日韩亚洲国产一区二区在线观看| 国产精品九九99| 国语自产精品视频在线第100页| 狂野欧美白嫩少妇大欣赏| 国产蜜桃级精品一区二区三区| 日韩欧美在线乱码| 母亲3免费完整高清在线观看| 伊人久久大香线蕉亚洲五| 在线播放国产精品三级| 亚洲狠狠婷婷综合久久图片| 日韩中文字幕欧美一区二区| 免费看日本二区| 亚洲精品美女久久av网站| 国内精品一区二区在线观看| 国产精品美女特级片免费视频播放器 | 狂野欧美白嫩少妇大欣赏| 国产主播在线观看一区二区| 一区二区三区国产精品乱码| 性色avwww在线观看| 亚洲成人精品中文字幕电影| 99热这里只有是精品50| 亚洲欧美日韩卡通动漫| 又粗又爽又猛毛片免费看| 色尼玛亚洲综合影院| 亚洲狠狠婷婷综合久久图片| 精品国产亚洲在线| 成年女人毛片免费观看观看9| 宅男免费午夜| 一个人看视频在线观看www免费 | 一级毛片女人18水好多| 91在线精品国自产拍蜜月 | 亚洲国产欧美人成| 国产精品99久久99久久久不卡| 久久久久性生活片| 男女床上黄色一级片免费看| 成人高潮视频无遮挡免费网站| 国产精品美女特级片免费视频播放器 | e午夜精品久久久久久久| 美女黄网站色视频| 国产精品乱码一区二三区的特点| 熟女人妻精品中文字幕| 黄色女人牲交| 国产免费男女视频| 操出白浆在线播放| 黄色女人牲交| 亚洲专区字幕在线| 1024手机看黄色片| 欧美在线黄色| 欧美高清成人免费视频www| 国产亚洲av嫩草精品影院| 国产不卡一卡二| 亚洲午夜理论影院| 免费观看的影片在线观看| 国产一区二区在线观看日韩 | 国内精品一区二区在线观看| 色吧在线观看| 不卡av一区二区三区| 亚洲人与动物交配视频| 免费av毛片视频| 免费看光身美女| 午夜激情欧美在线| 99精品欧美一区二区三区四区| 成人国产综合亚洲| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩精品亚洲av| 精品福利观看| 夜夜夜夜夜久久久久| 国产午夜精品久久久久久| 久久国产精品影院| 首页视频小说图片口味搜索|