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

    內爆炸載荷下泡沫鋁夾芯圓管塑性動力響應及能量耗散機理1)

    2023-11-16 06:42:30張?zhí)燧x劉志芳雷建銀王志華李世強
    力學學報 2023年10期
    關鍵詞:外管芯層內管

    張?zhí)燧x 劉志芳 雷建銀 王志華 李世強

    (太原理工大學機械與運載工程學院,太原 030024)

    (太原理工大學應用力學研究所,太原 030024)

    引言

    金屬與多孔材料結合構成的結構-功能一體化材料具有吸能性好、抗沖擊性強和輕質等優(yōu)點.在民用和軍事領域得到廣泛應用[1-9].

    長期以來,關于沖擊載荷作用下的結構變形機制和吸能特性的研究和取得的成果,已形成了基本的理論分析框架,也提出不同的實驗測試方法和數(shù)值計算模型[5,10-13].泡沫鋁方形夾芯板在爆炸載荷下動態(tài)響應的實驗測試和數(shù)值模擬研究表明: 泡沫芯層的能量吸收占比較大,且隨芯層密度和厚度的增大,能量吸收占比也隨之增大[5].為考慮爆炸載荷作用下夾芯結構的動力響應及分層梯度泡沫芯層分布對屈服軌跡和抗爆性的影響,可分別應用邊界屈服和內切屈服準則建立變形響應的近似解.實驗測試和數(shù)值模擬結果證明: 前者得出的結果與實驗測試和數(shù)值模擬結果吻合較好,后者較差[8].爆炸載荷作用下的3 階段理論模型考慮夾芯梁的一維流固耦合問題及結構的塑性彎曲和拉伸變形[14].剛塑性模型的假設多用于夾芯梁、板和殼的動力響應的分析,得到的理論預測值略高于實驗測試和有限元模擬結果[15-18].

    如何考慮梯度夾芯結構層間相互作用及結構的壓縮變形機制是理論建模的關鍵[19-24].材料迎爆面上受均布載荷作用的理論假設適用于起爆中心與結構距離較大的情形.在近距離爆炸載荷作用下須考慮結構動態(tài)響應的非線性特征.內爆載荷作用下泡沫夾芯圓柱殼結構動力響應分析屬于具有自由端邊界的近距離爆炸問題.涉及到強時空非線性和不均勻性[25-29].梯度泡沫鋁芯層圓柱殼的動力響應的有限元模擬表明: 在外管發(fā)生變形之前,梯度芯層壓潰變形表現(xiàn)出由低相對密度到高密度的逐層發(fā)展模式.而且梯度芯層的壓縮量大于均勻芯層的壓縮量[30].內爆載荷下梯度泡沫夾芯管的一維動力學模型將夾芯管的變形過程分為3 個階段: 流體與結構相互作用階段、芯層壓縮階段和外管的變形階段.分析表明: 炸藥量、芯層梯度與面板厚度對夾芯管動態(tài)響應產(chǎn)生影響;梯度泡沫夾芯管的吸能性能受內管壁厚的變化影響較大,受外管壁厚的變化影響較小;芯層的排布方式可改變結構的吸能性及外管變形特征[31].內爆載荷下夾芯圓管動態(tài)變形的4 階段理論模型考慮了外管未變形階段芯層在加速、減速運動下的壓縮、外管變形階段芯層二次壓縮和芯層密實化后結構整體變形的4 個變形過程.該理論解析模型能夠分析和解釋芯層壓縮和外管的變形階段應力的傳播過程.此外,基于相應的參數(shù)分析可對夾芯圓管進行進一步的優(yōu)化設計[32].

    已有內爆載荷作用下泡沫夾芯圓管動態(tài)響應的研究主要基于一維分析模型.但實際中人們更關注整體結構的響應.這就需要考慮環(huán)向膜力和軸向彎矩對夾芯圓管的變形行為和吸能機制的影響.

    作者所在課題組近期開展的研究工作涉及到泡沫鋁芯層軸向梯度分布對夾芯圓管的動態(tài)變形和吸能機制的影響[29].在實驗測試和數(shù)值模擬研究基礎上,假設內爆載荷具有空間呈三角形分布、時間服從指數(shù)衰減規(guī)律,建立考慮環(huán)向膜力和軸向彎矩對夾芯圓管動態(tài)響應的縱截面二維分析模型.通過選取不同幾何參數(shù),以無量綱形式確定了外管中心線的撓度,研究了夾芯圓管的內爆壓縮變形行為.如何在確定結構變形行為的基礎上進一步構造描述結構能量吸收機制的分析模型,并通過合理結構參數(shù)的選取獲得比吸能較大、圓管中心線撓度較小的最優(yōu)結構是一個亟待解決的問題.本文基于不同幾何參數(shù)和藥量的泡沫鋁夾芯圓管的內爆實驗觀察,確定圓管試件的變形模式.在考慮環(huán)向膜力和軸向彎矩的作用下,提出描述內爆載荷作用下內外管變形及吸能量機制的理論模型.基于三維Voronoi 算法構建了夾芯圓管有限元模型,研究幾何參數(shù)對結構的變形模式和抗爆性能的影響.

    1 實驗研究

    實驗試件的材料參數(shù)如表1 所示.不銹鋼的材料參數(shù)根據(jù)ISO 6892-1:2009 標準,使用拉伸測試方法獲得.利用線切割技術從100 mm 厚,相對密度為15%的閉孔泡沫鋁板上切割出泡沫鋁芯層,泡沫鋁基體材料為鋁合金.采用圖1 所示的實驗裝置進行泡沫鋁夾芯圓管內爆實驗.實驗裝置由球形乳化炸藥、泡沫填充夾層圓管、基座和支架組成,底座和支架與地面平行放置.為了避免爆炸后地面反射沖擊波的影響,泡沫填充的夾層圓管位于高度為1 m 的基座上.利用支架來保證炸藥在試件中心位置起爆,球形乳化炸藥的爆炸由雷管觸發(fā).將內外管與泡沫鋁芯層組裝到一起得到泡沫鋁夾芯圓管試樣,泡沫鋁與內外管之間無黏結.試件總高度L=100 mm,芯層壁厚為hC,內外管壁厚分別為hI和hE.實驗參數(shù)如表2 所示,其中,RI和RE分別為內外管半徑,W為乳化炸藥質量.芯層壁厚hC=RE-RI-hI由內外管半徑和內管壁厚確定.

    圖1 實驗裝置示意圖Fig.1 Schematic of the experimental setup

    表1 材料參數(shù)Table 1 Material parameters

    表2 實驗參數(shù)Table 2 Experiment parameters

    將表2 中7 組試樣進行內爆炸實驗,所得內外管跨中撓度如表3 所示,UI為內管撓度,UE為外管撓度.當炸藥質量相同時(Test 1-1 與Test 2-1、Test 1-2 與Test 2-2、Test 1-3 與Test 2-3),增大外管半徑使內管跨中撓度增大,芯層壓縮量增加,且外管跨中撓度減小.

    表3 實驗結果Table 3 Experiment results

    圖2 給出了試件1-2 的最終變形模式.從圖中可以看出,內外管的最終變形分為3 個區(qū)域: 塑性變形較大區(qū)域;圍繞塑性鉸剛性旋轉區(qū)域和邊界未變形區(qū)域[29].可見: 內外管的變形模式與文獻[33]報道的側向爆炸載荷作用下薄壁圓柱殼縱截面的變形模式相似;芯層的壓縮使內外管具有不同的未變形邊界.這個觀察結果為分析夾芯圓管在內爆載荷作用下的變形機制提供了科學依據(jù).

    圖2 試件1-2 實驗結果Fig.2 Experiment result of Test 1-2

    2 理論分析

    泡沫鋁夾芯圓管在內爆載荷下的動態(tài)響應可分為3 個階段: 爆炸載荷與內管相互作用,具有初始速度的內管壓縮芯層,芯層達到密實化后外管變形.取模型的軸向單位面積進行分析,其幾何尺寸及各個階段變形的受力分析如圖3 所示.其中,l為模型總高度的1/2,NI,NC和NE分別為內管、芯層及外管的環(huán)向膜力,Mx為軸向彎矩.

    圖3 結構變形過程: (a)軸向截面選取;(b)爆炸載荷作用于內管;(c)具有初始速度的內管壓縮芯層;(d)外管變形Fig.3 Deformation process of the structure: (a)Axial section selection;(b)Blast loading interacting with the internal tube;(c)Internal tube with the initial velocity impacting the foam core;(d)External tube deformation

    假設爆炸載荷在軸向空間呈三角形分布,在環(huán)向空間為均勻分布,在時間上為指數(shù)衰減.根據(jù)Henrych[34]對臨近爆炸載荷的分析,假設內爆炸下產(chǎn)生的脈沖載荷為

    其中,τ 為爆炸載荷持續(xù)時間,P0為爆炸載荷超壓峰值.關于爆炸載荷超壓峰值及超壓持續(xù)時間的計算方法,采用下式進行計算[35]

    其中,R為炸藥爆心到內管表面的距離,WT為TNT當量.根據(jù)山東圣世達化工有限責任公司提供的乳化炸藥當量實驗測試數(shù)據(jù),實驗用乳化炸藥與TNT炸藥質量關系式為WT=0.6W.

    基于泡沫夾芯圓管的變形特性,建立分析模型時采用了以下假設:

    (1)內外管假設為剛塑性模型,當材料達到屈服強度時發(fā)生塑性變形,卸載后結構不發(fā)生形變恢復;

    (2)忽略應變率和應變強化的影響,簡化理論推導過程,但會造成一定誤差;

    (3)假設結構在響應過程中以塑性變形為主,忽略內外管和芯層的破壞失效;

    (4)在芯層壓縮過程中,假設芯層在壓縮過程中以塑性變形為主并以恒定的平臺應力被壓縮.

    2.1 第1 階段: 爆炸載荷作用于內管

    根據(jù)圖3(a)中脈沖載荷函數(shù)分布及公式(1)中的函數(shù)假設,爆炸載荷作用于內管時內管獲得的總沖量為

    假設速度沿軸向分布為v0(1-x/l),內管的動量為

    其中,mI為內管的單位面積質量.結合式(7)和式(8)得出內管初速度

    2.2 第2 階段: 芯層壓縮

    如圖3(b)所示,根據(jù)爆炸載荷的大小,內管分為變形和未變形區(qū)域.當P0(1-x/l)≤NI時,內管不變形,當P0(1-x/l)>NI時,內管發(fā)生變形并壓縮芯層.內管膜力NI=σYIhI,其中,σYI為內管屈服強度.定義軸向方向的單位環(huán)向截面超壓=P0(1-x/l),在內管未變形區(qū)域,和NI存下以下關系

    假設內管受爆炸載荷產(chǎn)生的沖擊力與式(1)中的超壓在軸向的分布相同,內管的速度場近似等效為vI(x)=(1-x/l),uI為內管在初始位置隨時間變化的位移.在芯層壓縮階段,內管運動的微分方程為[36]

    其中,Qx為結構的剪力,σpl為芯層平臺應力,Nθ=NI+NC=N1,NC=σYChC為芯層環(huán)向膜力,σYC為芯層屈服強度.當x=0 時,剪力Qx=?Mx/?x=0 以滿足對稱性要求,且Mx=-M0以形成環(huán)向塑性鉸,M0為內管在純彎曲作用下的塑性彎矩.在x=l(1-NI/P0)=x1處,Mx=M0.因此,式(11)可變?yōu)?/p>

    初始條件:uI(0)=0,=P0τ/mI.因此可得

    式(14)可用于確定芯層壓縮階段某一時刻內管的位移和速度.在芯層壓縮過程中,若芯層壓縮至密實化,結構進入外管變形階段.因此,需要對內管最終位移與芯層量對比來確定結構在變形過程中是否進入外管變形階段.根據(jù)Li 等[37]提供的相對密度為15%泡沫鋁在準靜態(tài)壓縮載荷下的應力-應變曲線和Fan 等[38]對于多孔結構在壓縮過程中密實化的研究,確定15%泡沫鋁的密實化應變 εd為0.5.取公式(14)中=0,則內管在芯層壓縮階段的最終跨中撓度

    若uI1≤hCεd,外管不發(fā)生變形,泡沫夾芯管在芯層壓縮階段運動結束.若uI1>hCεd,泡沫芯層在第2 階段達到密實化,外管開始變形.此時式(15)中所得的內管撓度沒有考慮到外管變形的影響,對內管的最終撓度需進行進一步分析.

    2.3 第3 階段: 外管變形

    若外管發(fā)生變形,需確定芯層達到密實化時對應的時間,取uI=hCεd代入式(14)中,則

    在外管變形階段,芯層已達到密實化,內管繼續(xù)做減速運動.在外管變形階段,內管和泡沫芯層的運動使外管開始做加速運動.當外管速度與內管和芯層速度相同時,結構做整體減速運動直至變形結束.以內管的運動為研究對象,內管的加速度為

    其中,N2=NI+NC+NE,外管膜力NE=σYEhE,σYE為外管屈服強度.初始條件:uI2(t2)=hCεd,(t2)=-k1t2+P0τ/mI.結合式(18)和外管變形階段的初始條件得出內管在第3 階段的速度和位移

    當內管速度為0 時,結構整體運動結束,取=0代入式(19)中可得出結構變形結束所需時間

    在理論計算中忽略了內外管和芯層的耦合關系,因此,外管的位移uE=uI-hCεd.內管對應的變形輪廓

    取uf=hCεd代入式(22)中可得外管軸向變形邊界

    根據(jù)動量守恒定律,在外管變形階段,內/外管和泡沫芯層的共同速度為

    其中,mC和mE分別為泡沫芯層和外管的單位面積質量.結合式(17)和式(19),得到內外管和泡沫芯達到此速度所需的時間.

    2.4 能量吸收機制

    基于上述對于泡沫鋁夾芯圓管在內爆炸載荷下的3 階段變形分析,結合結構在響應過程中動能的變化及環(huán)向膜力和軸向彎矩的作用對其在響應過程中的能量吸收機理進行理論分析.如圖3(a)所示,以泡沫鋁夾芯圓管軸向單位面積為研究對象,內管的初動能為

    當內外管和芯層達到共同速度之后,結構開始做整體的減速運動.由動量守恒得出內外管和芯層的共同速度v1

    根據(jù)內外管在軸向方向上的變形位置x1和x2可以得出內外管達到共同速度時的動能

    內外管和芯層達到共同速度時泡沫鋁夾芯圓管軸向單位面積的能量吸收E1=-.

    在結構整體運動階段,由環(huán)向膜力和芯層動平臺應力所消耗在結構中的軸向單位面積的總能量E2為

    綜合考慮泡沫鋁夾芯圓管受到的環(huán)向膜力和軸向彎矩的影響可以得出,結構在整個運動過程中軸向單位面積吸收的總能量表達式

    根據(jù)表2 中實驗所得內外管最終撓度可知,泡沫鋁夾芯圓管在內爆炸載荷作用下變形較小.在內外管和芯層達到共同速度到整個泡沫夾芯管達到最終變形的過程中,由內外管和芯層的彎矩引起的能量吸收很小,可忽略不計.因此,式(32)可改寫為

    3 數(shù)值模擬及算法驗證

    3.1 有限元模型

    采用商業(yè)軟件ABAQUS/Explicit 有限元程序進行數(shù)值模擬.基于結構和載荷的對稱性,建立了泡沫鋁夾芯圓管的1/8 模型,如圖4 所示.1/8 模型的邊界條件定義為: 結構頂部為自由邊界,在其他邊緣施加對稱邊界條件(X-Z面約束Y方向位移和X,Z方向轉角,Y-Z面約束X方向位移和Y,Z方向轉角).內外管采用8 節(jié)點線性實體單元(C3D8R),利用三維Voronoi 算法[39]構建相對密度為15%泡沫鋁芯層,采用傳統(tǒng)的3 節(jié)點殼單元(S3R)進行網(wǎng)格劃分.為保持結構變形的穩(wěn)定性,在ABAQUS/explicit 中采用罰函數(shù)方法來加強約束,內外管與泡沫芯層之間定義硬接觸(hard contact),摩擦系數(shù)設置為0.02[39].采用表1 中的內外管和芯層的材料參數(shù)進行數(shù)值模擬.泡沫鋁芯層采用雙線性強化模型,內外管均為不銹鋼材料,在模擬中使用各向同性隨動強化塑性模型,應變率 ε 采用Cowper-Symonds 模型,屈服應力為 1 +(ε/C)1/q,其中C=100 s-1,q=10[33].球形爆炸載荷采用ConWep 算法來模擬結構變形.芯層和內外管的網(wǎng)狀尺寸定義為0.5 mm,進一步的細化并沒有明顯提高計算結果的精度.

    圖4 有限元模型Fig.4 Finite element model

    3.2 模型驗證

    圖5(a)為試件1-2 和試件2-3 實驗與數(shù)值模擬最終變形模態(tài)對比,數(shù)值模擬結果具有明顯的3 個變形區(qū)域,且內外管的未變形邊界不同,與圖1 實驗模態(tài)中觀察到的結果一致,二者吻合較好.以內外管中點撓度為對象來驗證數(shù)值模擬與理論分析模型的有效性,所得結果如圖5(b)所示.

    圖5 模型驗證Fig.5 Validation of the model

    在數(shù)值模擬中忽略了泡沫芯層在變形過程中的破壞以及內外管與芯層之間的耦合作用,低估結構的強度,導致大多數(shù)模擬結果略高于實驗值.在理論分析過程中,假設泡沫芯層達到密實化后,外管開始變形,已經(jīng)密實化的芯層不再進一步發(fā)生斷裂和失效破壞.此外,在結構3 階段變形過程中,采用解耦的辦法進行逐步分析,忽略了內外管與芯層之間的流固耦合作用.且理論分析中忽略了應變率和應變強化的影響,大多數(shù)理論預測的中點撓度略高于實驗結果.

    4 結果與討論

    4.1 參數(shù)無量綱化

    根據(jù)結構的幾何模型和材料強度進行無量綱參數(shù)的定義

    將無量綱參數(shù)代入式(13)和式(18)可得

    將內管和外管的跨中無量綱撓度定義為

    結合式(16)、式(20)、式(21)、式(35)和式(36)可得內外管中點無量綱撓度的最終表達式

    將無量綱參數(shù)代入式(25)、式(27)和式(28)并消除結構質量的影響,可得結構在各響應階段的比吸能

    4.2 參數(shù)分析

    圖6 列出了給定內外管和芯層材料參數(shù)的泡沫鋁夾芯圓管在TNT 藥量為30 g 下外管中點無量綱撓度隨結構壁厚變化的響應體.可以看出,利用泡沫芯層和內外管的壁厚可以有效地評估外管的無量綱撓度.內管壁厚對外管中點撓度影響較大,芯層和外管壁厚對外管中點撓度影響較小.當芯層和外管壁厚一定時,外管中點撓度隨內管壁厚的增加而減小.當 β1=0.053,β2=0.105,β3=0.053 時,外管無量綱中點撓度出現(xiàn)最大值,為0.3.由于外管撓度隨幾何參數(shù)的變化呈單調趨勢,因此,需結合總能量吸收來進一步得出結構的最優(yōu)構型.

    圖6 外管無量綱中點撓度響應體Fig.6 Response volume of nondimensional midpoint deflection of the external tube

    取內管半徑RI=38~44.5 mm,外管半徑RE=44.5~54 mm,內管壁厚hI=2~4 mm,外管壁厚hE=2 mm,在TNT 藥量為18 g 來研究結構比吸能和外觀中點撓度的最優(yōu)解,所得結果如圖7 所示.當內管半徑RI=38 mm,外管半徑RE=44.5 mm,內管壁厚hI=2 mm 時比吸能最大,為2.39 kJ/kg,此時對應的外管無量綱中點撓度為0.13.內管壁厚及內外管半徑的變化使結構在既定的炸藥質量下出現(xiàn)了許多外管無量綱撓度為0 的幾何構型.在外管無量綱撓度為0 的條件下,內管半徑RI=38 mm,外管半徑RE=53.5 mm,內管壁厚hI=2 mm 的比吸能最大,為1.15 kJ/kg.根據(jù)圖7 所列出的數(shù)據(jù),可以在實際工程中給定的外管撓度條件下找出比吸能最大的幾何構型.此外,也可以在給定比吸能條件下找出外管撓度最小的結構.

    圖7 比吸能與外管無量綱撓度最優(yōu)解Fig.7 The optimal solution of specific energy absorption and nondimensional midpoint deflection of the external tube

    4.3 梯度芯層

    如圖1 所示,芯層的相對密度沿軸向每20 mm變化一次.根據(jù)1/8 有限元模型的對稱性,在數(shù)值模擬中相對密度 ρs1的高度為10 mm,如圖8 所示.芯層的軸向梯度分布模型參數(shù)如表4 所示.其中,G0為非梯度模型對照組.

    圖8 不同泡沫芯層分布的比吸能和外管無量綱撓度Fig.8 Specific energy absorption and nondimensional deflection of external tube for various foam core distributions

    表4 梯度模型參數(shù)Table 4 Graded model parameters

    以內管半徑RI=38 mm,外管半徑RE=44.5 mm,內管壁厚hI=2 mm,外管壁厚hE=3 mm,TNT 藥量為30 g 為基礎來研究芯層軸向梯度分布對結構抗爆性能的影響,所得結果如圖8 所示.G5 梯度模型具有最高的比吸能,且外管無量綱中點撓度與其他混合梯度結構相同,但芯層吸能占比僅為33%.對于G1 梯度模型,外管無量綱中點撓度低于其他軸向梯度分布結構,且芯層能量吸收占比為38%.

    5 結論與展望

    采用實驗、理論和數(shù)值模擬相結合的方法,研究了泡沫鋁夾芯圓管在內爆炸載荷作用下的動力響應.建立了泡沫鋁夾芯圓管內爆炸響應的解析表達式.采用理論和有限元分析方法對泡沫鋁夾芯圓管的動力學行為、變形機理和能量吸收進行了詳細分析.詳細討論了外管半徑、內外管壁厚、炸藥量和芯層梯度分布對結構力學響應的影響.本文得出主要結論如下.

    (1)理論分析模型與泡沫鋁夾芯圓管在內爆炸載荷作用下的實驗和數(shù)值模擬結果吻合較好.所提出的動態(tài)顯示算法具有較好的可靠性.

    (2)通過內爆炸實驗得到了泡沫鋁夾芯圓管在內爆炸載荷作用下3 個不同變形區(qū)域: 塑性變形較大區(qū)域、圍繞塑性鉸剛性旋轉區(qū)域和邊界未變形區(qū)域.

    (3)內管壁厚對外管中點撓度影響較大,芯層和外管壁厚對外管中點撓度影響較小.通過改變結構的幾何參數(shù)得出了結構比吸能和外管中點撓度的最優(yōu)解集.

    (4)在所有芯層軸向梯度分布的泡沫鋁夾芯圓管中,芯層軸向分布由中心到邊緣遞減的梯度夾芯圓管模型的外管無量綱撓度最小,芯層能量吸收占比最大.

    本文所建立的理論解析模型中,內外管在受載變形時假設為理想剛塑性模型.而基體材料本身存在應力強化過程,且泡沫鋁夾芯圓管在受內爆炸載荷作用下的彈性變形效應[40]對結構的整體響應機制存在一定的影響,這也是所提出的理論預測值與實驗存在誤差的原因.關于如何結合基體材料的應力強化和結構的彈性效應來構建更為精確的動態(tài)顯示理論計算方法可展開更為深入的研究.

    文中未涉及關于徑向梯度泡沫鋁夾芯圓管在內爆炸載荷下的理論分析,徑向梯度芯層壓縮階段需考慮每一層泡沫芯的相互作用,可以假設梯度芯層為解耦的逐層壓潰模式[31].但如何確定徑向梯度的每一層芯層的變形邊界及應力波在梯度芯層中的傳播規(guī)律是需要解決的重點和難點.

    猜你喜歡
    外管芯層內管
    套管換熱器穩(wěn)定性設計計算的探討
    化肥設計(2022年2期)2022-04-28 09:23:52
    一種耐高溫能力強的橡膠輸送帶
    外管新規(guī)助力跨境支付業(yè)務升級提速
    中國外匯(2019年10期)2019-08-27 01:58:02
    空中爆炸載荷下梯度波紋夾層板抗爆性能仿真研究
    兵工學報(2017年6期)2017-07-10 10:28:52
    車用新型結構的夾層構件
    汽車文摘(2017年4期)2017-04-27 02:05:07
    下期要目
    化工設計(2017年2期)2017-03-05 10:04:30
    一種用于氣化系統(tǒng)的進料噴射器
    石油化工(2015年9期)2015-08-15 00:43:05
    煙草專賣內管委派制對縣級局內部監(jiān)管的思考
    新形勢下如何發(fā)揮專賣內管委派制的作用和效能
    芯層微孔發(fā)泡管材的制備
    久久精品国产亚洲av香蕉五月| 国产午夜精品论理片| av黄色大香蕉| 很黄的视频免费| 国产精品一区www在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 人妻久久中文字幕网| 桃色一区二区三区在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩高清专用| 免费看美女性在线毛片视频| 丰满乱子伦码专区| 免费电影在线观看免费观看| 97人妻精品一区二区三区麻豆| 国产高清视频在线播放一区| 国产精品1区2区在线观看.| 伊人久久精品亚洲午夜| 天堂av国产一区二区熟女人妻| 久久久国产成人精品二区| 久久久久九九精品影院| 舔av片在线| 国产精品人妻久久久影院| 日韩欧美一区二区三区在线观看| 亚洲精品日韩av片在线观看| 亚洲精品一区av在线观看| 国产三级中文精品| 日本免费一区二区三区高清不卡| 亚洲电影在线观看av| 国产精品1区2区在线观看.| 成人精品一区二区免费| 精品久久久久久成人av| 啦啦啦啦在线视频资源| 国产美女午夜福利| 久久久色成人| 国产亚洲av嫩草精品影院| 桃红色精品国产亚洲av| 狠狠狠狠99中文字幕| 男女边吃奶边做爰视频| 亚州av有码| www.色视频.com| 最好的美女福利视频网| 欧美又色又爽又黄视频| 99久久久亚洲精品蜜臀av| 中文亚洲av片在线观看爽| 麻豆成人av在线观看| 国产白丝娇喘喷水9色精品| 中文字幕av在线有码专区| 真实男女啪啪啪动态图| 国内久久婷婷六月综合欲色啪| 国产精品永久免费网站| 久久精品国产99精品国产亚洲性色| 麻豆av噜噜一区二区三区| 国产欧美日韩一区二区精品| 最近在线观看免费完整版| 一本一本综合久久| 国产欧美日韩一区二区精品| 国产淫片久久久久久久久| 可以在线观看的亚洲视频| 国产成人aa在线观看| 在线免费观看不下载黄p国产 | 午夜福利在线观看免费完整高清在 | 精品一区二区三区人妻视频| 欧美三级亚洲精品| 久久午夜福利片| 中国美女看黄片| 午夜免费激情av| 熟女人妻精品中文字幕| 日本撒尿小便嘘嘘汇集6| 国产精品伦人一区二区| 亚洲精品乱码久久久v下载方式| 日本三级黄在线观看| 精品人妻一区二区三区麻豆 | 一边摸一边抽搐一进一小说| 毛片一级片免费看久久久久 | 国内久久婷婷六月综合欲色啪| 日本一本二区三区精品| 精品国内亚洲2022精品成人| 日韩欧美在线乱码| 天天一区二区日本电影三级| 国产白丝娇喘喷水9色精品| 国产视频一区二区在线看| 亚洲中文字幕一区二区三区有码在线看| 日韩人妻高清精品专区| 亚洲精品亚洲一区二区| av在线观看视频网站免费| 3wmmmm亚洲av在线观看| 精品一区二区免费观看| 日本精品一区二区三区蜜桃| 国产黄a三级三级三级人| 欧美性猛交黑人性爽| 天堂动漫精品| 在线免费观看的www视频| 麻豆国产97在线/欧美| 18禁在线播放成人免费| 久久久久久久精品吃奶| 日韩欧美精品免费久久| 亚洲四区av| 精品一区二区三区视频在线| 亚洲成人久久爱视频| 99热这里只有是精品在线观看| 国产精品伦人一区二区| 久久精品国产鲁丝片午夜精品 | 人人妻,人人澡人人爽秒播| 亚洲精品在线观看二区| av天堂在线播放| 热99在线观看视频| 少妇被粗大猛烈的视频| 国内揄拍国产精品人妻在线| 日韩中字成人| 18禁黄网站禁片免费观看直播| 日韩中字成人| 国产一区二区三区av在线 | 1024手机看黄色片| 久久久久性生活片| 桃色一区二区三区在线观看| 精品无人区乱码1区二区| 一级av片app| 精品99又大又爽又粗少妇毛片 | 亚洲精品成人久久久久久| 村上凉子中文字幕在线| 搡老岳熟女国产| 亚洲成人精品中文字幕电影| 成人三级黄色视频| 亚洲精品日韩av片在线观看| 2021天堂中文幕一二区在线观| 1000部很黄的大片| 日本五十路高清| 欧美日韩精品成人综合77777| 欧美日本亚洲视频在线播放| 亚洲精品一区av在线观看| 久久国产乱子免费精品| 午夜免费成人在线视频| ponron亚洲| 日本免费一区二区三区高清不卡| 国产精品无大码| 日韩精品有码人妻一区| 丰满的人妻完整版| 免费av毛片视频| 成人国产麻豆网| 亚州av有码| 亚洲最大成人av| 亚洲av第一区精品v没综合| 日日干狠狠操夜夜爽| 99国产精品一区二区蜜桃av| 久久国内精品自在自线图片| 亚洲国产色片| 成人鲁丝片一二三区免费| 欧美成人免费av一区二区三区| 日韩欧美在线二视频| 久久久精品大字幕| 美女被艹到高潮喷水动态| netflix在线观看网站| 黄色一级大片看看| 国产老妇女一区| 一区二区三区高清视频在线| 久久欧美精品欧美久久欧美| 人人妻人人澡欧美一区二区| 男人舔奶头视频| 久久人人精品亚洲av| 在线观看一区二区三区| 免费av毛片视频| 久久中文看片网| 欧美一区二区国产精品久久精品| 琪琪午夜伦伦电影理论片6080| 免费大片18禁| 俄罗斯特黄特色一大片| 久久久久久久精品吃奶| 毛片一级片免费看久久久久 | 午夜福利成人在线免费观看| 精品99又大又爽又粗少妇毛片 | 精品一区二区免费观看| 麻豆国产97在线/欧美| 免费看美女性在线毛片视频| 午夜老司机福利剧场| 嫩草影院入口| 禁无遮挡网站| 简卡轻食公司| 内地一区二区视频在线| 午夜免费男女啪啪视频观看 | 一个人看视频在线观看www免费| 一个人看的www免费观看视频| 在线天堂最新版资源| 日本免费一区二区三区高清不卡| 精品人妻1区二区| 搡女人真爽免费视频火全软件 | 国产一级毛片七仙女欲春2| 亚洲美女黄片视频| 精品国产三级普通话版| 男女边吃奶边做爰视频| 一夜夜www| 真人一进一出gif抽搐免费| 韩国av在线不卡| 亚洲欧美精品综合久久99| 国产久久久一区二区三区| 中文字幕av在线有码专区| 动漫黄色视频在线观看| 国产不卡一卡二| 九九热线精品视视频播放| 午夜福利视频1000在线观看| 亚洲成a人片在线一区二区| 成人美女网站在线观看视频| 国内精品一区二区在线观看| 久久久午夜欧美精品| 丰满的人妻完整版| 亚洲狠狠婷婷综合久久图片| 亚洲av熟女| 日本撒尿小便嘘嘘汇集6| 国产大屁股一区二区在线视频| 国产伦在线观看视频一区| 国内少妇人妻偷人精品xxx网站| 99久久精品热视频| 嫩草影院精品99| 99久久精品国产国产毛片| 啦啦啦啦在线视频资源| 久久久精品大字幕| 国语自产精品视频在线第100页| 我的老师免费观看完整版| 午夜免费成人在线视频| 国产精品av视频在线免费观看| 美女高潮喷水抽搐中文字幕| 欧美性感艳星| 日韩av在线大香蕉| 日韩亚洲欧美综合| 久久99热6这里只有精品| 亚洲国产精品久久男人天堂| 熟妇人妻久久中文字幕3abv| 99视频精品全部免费 在线| av在线观看视频网站免费| 九九久久精品国产亚洲av麻豆| 国产一区二区亚洲精品在线观看| 欧美在线一区亚洲| 99精品久久久久人妻精品| 男人狂女人下面高潮的视频| 免费黄网站久久成人精品| 久久人妻av系列| 亚洲国产日韩欧美精品在线观看| 亚洲电影在线观看av| 人妻夜夜爽99麻豆av| 亚洲精品在线观看二区| 国产伦精品一区二区三区视频9| 国产亚洲91精品色在线| 国内毛片毛片毛片毛片毛片| 日韩欧美精品免费久久| 午夜激情福利司机影院| 日本三级黄在线观看| 国产激情偷乱视频一区二区| 黄色配什么色好看| 国产一级毛片七仙女欲春2| 亚洲国产精品sss在线观看| 亚洲综合色惰| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av第一区精品v没综合| www.www免费av| 久久久久性生活片| 老女人水多毛片| 国产伦在线观看视频一区| 99视频精品全部免费 在线| 国产精品综合久久久久久久免费| 精华霜和精华液先用哪个| 欧美日韩精品成人综合77777| 哪里可以看免费的av片| 69人妻影院| 国产淫片久久久久久久久| 精品99又大又爽又粗少妇毛片 | 搡女人真爽免费视频火全软件 | av.在线天堂| 无遮挡黄片免费观看| 色吧在线观看| 91av网一区二区| 亚洲欧美日韩无卡精品| 亚洲精品色激情综合| 听说在线观看完整版免费高清| 麻豆国产av国片精品| 99热这里只有是精品在线观看| а√天堂www在线а√下载| 久久天躁狠狠躁夜夜2o2o| 日韩大尺度精品在线看网址| 人人妻人人澡欧美一区二区| 亚洲av成人av| 中文字幕精品亚洲无线码一区| 国产伦一二天堂av在线观看| 永久网站在线| 亚洲18禁久久av| 可以在线观看的亚洲视频| 久久久久久久久久久丰满 | 亚洲精品亚洲一区二区| 亚洲性夜色夜夜综合| 亚洲av电影不卡..在线观看| 一级黄色大片毛片| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 尾随美女入室| 久久久成人免费电影| 久久久色成人| 免费av观看视频| 精品人妻视频免费看| 免费不卡的大黄色大毛片视频在线观看 | 国产在视频线在精品| 久久热精品热| 18禁在线播放成人免费| 亚洲国产精品久久男人天堂| 国产毛片a区久久久久| 97超视频在线观看视频| 免费在线观看成人毛片| 亚洲av二区三区四区| 九色成人免费人妻av| 国产精品伦人一区二区| 成人av在线播放网站| 少妇人妻精品综合一区二区 | 黄片wwwwww| 国内毛片毛片毛片毛片毛片| 内地一区二区视频在线| 色哟哟哟哟哟哟| 美女xxoo啪啪120秒动态图| av黄色大香蕉| 国内少妇人妻偷人精品xxx网站| 国国产精品蜜臀av免费| 极品教师在线视频| 精品一区二区三区人妻视频| 99九九线精品视频在线观看视频| 99热这里只有精品一区| 三级男女做爰猛烈吃奶摸视频| 久久久成人免费电影| 最近在线观看免费完整版| 亚洲国产日韩欧美精品在线观看| 少妇被粗大猛烈的视频| 嫩草影院精品99| 如何舔出高潮| 伦理电影大哥的女人| 婷婷六月久久综合丁香| 人妻少妇偷人精品九色| 欧美不卡视频在线免费观看| 久久久久久国产a免费观看| 又爽又黄a免费视频| 午夜福利18| 国产淫片久久久久久久久| 少妇高潮的动态图| 狂野欧美激情性xxxx在线观看| 日韩欧美国产在线观看| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线| 精品福利观看| 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看| 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| 国产一级毛片七仙女欲春2| 久久久久久久精品吃奶| 国产亚洲精品综合一区在线观看| 波多野结衣高清作品| 精品99又大又爽又粗少妇毛片 | 非洲黑人性xxxx精品又粗又长| 看免费成人av毛片| a级毛片a级免费在线| 国产在线精品亚洲第一网站| av国产免费在线观看| 97超级碰碰碰精品色视频在线观看| 永久网站在线| 午夜福利在线观看吧| 欧美激情久久久久久爽电影| 久久久久久久午夜电影| av在线老鸭窝| 国产精品av视频在线免费观看| 啦啦啦啦在线视频资源| 麻豆久久精品国产亚洲av| 69人妻影院| 国产探花在线观看一区二区| 日本 av在线| 国产美女午夜福利| h日本视频在线播放| 噜噜噜噜噜久久久久久91| 久久精品国产鲁丝片午夜精品 | 中文字幕精品亚洲无线码一区| 在线观看66精品国产| or卡值多少钱| 十八禁国产超污无遮挡网站| 欧美日韩精品成人综合77777| 国产成人aa在线观看| 国产乱人伦免费视频| 丰满乱子伦码专区| 又黄又爽又刺激的免费视频.| 日韩欧美精品v在线| 久99久视频精品免费| 又紧又爽又黄一区二区| 岛国在线免费视频观看| 天堂网av新在线| 日本五十路高清| 午夜免费激情av| 男女那种视频在线观看| 午夜福利在线在线| 在现免费观看毛片| 国产精品一区二区三区四区免费观看 | 999久久久精品免费观看国产| а√天堂www在线а√下载| 午夜视频国产福利| 日韩精品有码人妻一区| 日韩亚洲欧美综合| 久久久久久久久中文| 国产真实伦视频高清在线观看 | 中国美白少妇内射xxxbb| 久久精品久久久久久噜噜老黄 | 在线免费观看不下载黄p国产 | 久久天躁狠狠躁夜夜2o2o| 人人妻人人看人人澡| 色噜噜av男人的天堂激情| 成人av在线播放网站| 欧美最新免费一区二区三区| 亚洲无线在线观看| 在线播放国产精品三级| 亚洲真实伦在线观看| 久久午夜亚洲精品久久| 天天一区二区日本电影三级| 有码 亚洲区| 日本三级黄在线观看| 97碰自拍视频| 97热精品久久久久久| 亚洲午夜理论影院| 五月玫瑰六月丁香| 欧美高清性xxxxhd video| 欧美+亚洲+日韩+国产| 午夜福利视频1000在线观看| 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影| 精品久久久久久久久久免费视频| 久久精品国产自在天天线| 国产美女午夜福利| 日本-黄色视频高清免费观看| 欧美激情在线99| 日韩在线高清观看一区二区三区 | 国产视频内射| 成年版毛片免费区| 身体一侧抽搐| videossex国产| 国产av麻豆久久久久久久| 日韩欧美一区二区三区在线观看| 国产亚洲精品久久久com| 黄色日韩在线| av视频在线观看入口| 看黄色毛片网站| 日韩欧美精品v在线| 欧美xxxx性猛交bbbb| 九九在线视频观看精品| 男女之事视频高清在线观看| 国产在线男女| 99久久中文字幕三级久久日本| 女的被弄到高潮叫床怎么办 | 欧美激情在线99| 久久久久免费精品人妻一区二区| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 观看美女的网站| 精品福利观看| 1024手机看黄色片| 亚洲专区中文字幕在线| 极品教师在线视频| 亚洲综合色惰| 精品国内亚洲2022精品成人| 亚洲性夜色夜夜综合| 黄色一级大片看看| 嫁个100分男人电影在线观看| 免费在线观看影片大全网站| 啦啦啦啦在线视频资源| 丰满的人妻完整版| 联通29元200g的流量卡| 99九九线精品视频在线观看视频| 一个人看的www免费观看视频| 久久久久九九精品影院| 日本免费一区二区三区高清不卡| 国产69精品久久久久777片| 3wmmmm亚洲av在线观看| 最近最新中文字幕大全电影3| 99在线人妻在线中文字幕| 亚洲无线观看免费| 一个人观看的视频www高清免费观看| 久久香蕉精品热| 别揉我奶头~嗯~啊~动态视频| 亚洲乱码一区二区免费版| 色av中文字幕| 美女高潮喷水抽搐中文字幕| 国产亚洲精品av在线| 少妇猛男粗大的猛烈进出视频 | 天堂av国产一区二区熟女人妻| 我要搜黄色片| 在线观看66精品国产| 午夜爱爱视频在线播放| www日本黄色视频网| av黄色大香蕉| 91久久精品国产一区二区三区| 欧美精品啪啪一区二区三区| 毛片女人毛片| 国产av在哪里看| 别揉我奶头 嗯啊视频| 亚洲国产欧洲综合997久久,| 国产精品99久久久久久久久| 国产精品久久视频播放| 22中文网久久字幕| 久久草成人影院| 丰满人妻一区二区三区视频av| 天堂av国产一区二区熟女人妻| 国产高清有码在线观看视频| 亚洲一区二区三区色噜噜| 国产日本99.免费观看| 一a级毛片在线观看| 亚洲av电影不卡..在线观看| 久久精品国产自在天天线| 国产精品久久电影中文字幕| 欧美日韩精品成人综合77777| 香蕉av资源在线| 色视频www国产| 免费观看精品视频网站| 男女边吃奶边做爰视频| 人人妻人人看人人澡| 搡老妇女老女人老熟妇| 波野结衣二区三区在线| 亚洲精品久久国产高清桃花| 国语自产精品视频在线第100页| 一区二区三区免费毛片| 国产老妇女一区| 日韩高清综合在线| 999久久久精品免费观看国产| 日韩亚洲欧美综合| 亚洲专区中文字幕在线| 欧美色欧美亚洲另类二区| 午夜福利在线在线| 免费看日本二区| av天堂在线播放| 国产麻豆成人av免费视频| 午夜免费成人在线视频| 免费看av在线观看网站| avwww免费| 欧美成人免费av一区二区三区| 久久这里只有精品中国| 别揉我奶头 嗯啊视频| 九九久久精品国产亚洲av麻豆| 亚洲国产精品久久男人天堂| 亚州av有码| av专区在线播放| 成人特级av手机在线观看| 亚洲中文日韩欧美视频| 听说在线观看完整版免费高清| 亚洲va日本ⅴa欧美va伊人久久| 欧美性猛交黑人性爽| 人人妻人人看人人澡| 亚洲av第一区精品v没综合| 国产白丝娇喘喷水9色精品| 嫩草影院入口| 亚洲成人免费电影在线观看| 日本欧美国产在线视频| 人妻夜夜爽99麻豆av| 99riav亚洲国产免费| 在线天堂最新版资源| 俺也久久电影网| 国产精品国产高清国产av| 91在线精品国自产拍蜜月| 91狼人影院| а√天堂www在线а√下载| 亚洲色图av天堂| 一区二区三区免费毛片| 欧美三级亚洲精品| 1024手机看黄色片| 久久欧美精品欧美久久欧美| 久久精品国产亚洲网站| 深夜a级毛片| 亚洲男人的天堂狠狠| 搡女人真爽免费视频火全软件 | 制服丝袜大香蕉在线| 黄色一级大片看看| 国产精品久久视频播放| 亚洲第一区二区三区不卡| 熟女人妻精品中文字幕| 日本与韩国留学比较| 老司机午夜福利在线观看视频| 国产主播在线观看一区二区| 级片在线观看| 波多野结衣高清作品| 亚洲中文日韩欧美视频| 日韩欧美精品免费久久| 少妇人妻一区二区三区视频| 欧美潮喷喷水| 国模一区二区三区四区视频| 午夜影院日韩av| 亚洲人成网站高清观看| 亚洲av第一区精品v没综合| 午夜福利高清视频| 级片在线观看| 男女之事视频高清在线观看| 国内少妇人妻偷人精品xxx网站| 九九爱精品视频在线观看| 久久久久久大精品| 国产精品亚洲一级av第二区| 蜜桃亚洲精品一区二区三区| 国产一区二区亚洲精品在线观看| 欧美区成人在线视频| 蜜桃亚洲精品一区二区三区| 国产一区二区亚洲精品在线观看| 小说图片视频综合网站| 欧美中文日本在线观看视频| 久久精品影院6| av福利片在线观看| 又爽又黄无遮挡网站| 精品久久国产蜜桃| 国产精品亚洲一级av第二区| 如何舔出高潮| 淫秽高清视频在线观看| 熟女电影av网| 91精品国产九色| 精品久久国产蜜桃| 中文亚洲av片在线观看爽| 女的被弄到高潮叫床怎么办 | 伦理电影大哥的女人| 国产日本99.免费观看| 免费高清视频大片| 看十八女毛片水多多多| 桃色一区二区三区在线观看| 少妇的逼好多水|