Gasoline vapor/air premixed flame's unstretched laminar burning velocity
-
摘要: 为了研究汽油蒸气/空气的爆炸特性，介绍了汽油蒸气/空气预混火焰的无拉伸层流燃烧速率。通过实验研究了无拉伸层流燃烧速率，与汽油的主要组分异辛烷和正庚烷与空气的预混火焰层流燃烧速率做了对比，发现汽油蒸气/空气的无拉伸层流燃烧速率小于异辛烷和正庚烷与空气的预混火焰无拉伸层流燃烧速率，但无拉伸层流燃烧速率随当量比的变化规律相同，随着当量比增大，无拉伸层流燃烧速率呈先增大、再减小的变化趋势，最大值在当量比为1处取得。Abstract: In order to characterize the explosion of the gasoline vapor/air mixture, we introduced the unstretched laminar burning velocity of premixed gasoline vapor/air flame to the research of the mixture. The unstretched laminar burning velocity was experimentally studied and compared with those of two mixture gases consisting of such major compositions of gasoline as isooctane and n-heptane, respectively. The results show that the unstretched laminar burning velocity of the prepared gasoline vapor/air mixture is lower than those of the isooctane/air mixture gas and the n-heptane/air misture gas, but it takes on a similar tendency to change with the equivalence ratio, i.e. the unstretched laminar burning velocity firstly increases and then decreases with the increase of the equivalence ratio, and the maximum value can be obtained at the equivalence ratio of 1.
-
Key words:
- gasoline vapor /
- deflgration /
- premixed flame /
- overpressure /
- unstretched laminar /
- burning velocity
-
大型商用飞机撞击核电厂是一种潜在的威胁, 9·11事件以后尤其得到重视。对核电厂而言, 一旦发生破坏, 造成的危害不堪设想, 目前在核岛结构的设计中考虑大型商用飞机撞击的影响已经纳入规范。美国核管会在2009年将评价核设施结构防御飞机撞击能力作为联邦法规的新条款颁布[1], 同年发布NEI 07-13[2], 提供评价大型商用飞机撞击下核设施结构的响应的方法, 并在2011年正式发布了RG 1.217[3], 要求核设施申请者应提交大型商用飞机撞击下核设施结构响应的评价报告。在我国, 《核电厂厂址选择的外部人为事件》[4]中尚未强制要求评价核电厂结构防御大型商用飞机撞击能力, 但福岛核事故后, 对于一些超设计基准事件, 包括大型商用飞机撞击对核电厂的影响, 也越来越受到重视[5]。新项目审批重新开启时要求新项目须符合三代标准, 提高准入门槛, 按照国际最高安全要求新建核电项目。这说明, 新建的项目除应满足我国规范体系外, 还应按较高的国际标准进行校核验算, 这就要求新建的三代核电厂须评价核岛结构防御飞机撞击的能力。
近年来, 人们开展了在飞机撞击下混凝土结构响应的研究[6-7, 11-13], 这些研究多借鉴了其他冲击荷载(如弹体、弹丸及爆炸等)下混凝土结构响应的研究成果[6-10], 而飞机撞击混凝土结构的数值模拟及动力学分析仍需进一步深入。由于大型商用飞机撞击核电厂结构会产生两种效应, 局部破坏和整体振动, 即由飞机发动机引起的局部破坏(穿入、散裂、碎甲和穿透)和由整个飞机冲击引起的整体结构的反应分析, 因此, 需要进行两种典型的结构破坏模式的评估。由于飞机撞击混凝土结构具有短时强载的特征, 其破坏机制非常复杂[12-13]。整个撞击过程是一个瞬态、高应变率、复杂的破坏过程, 包含着几何、材料非线性方面的诸多问题, 也涉及接触面上伴随接触-撞击过程而产生的撞击能量的传递与交换、接触面的滑移、分离等非线性接触问题。在进行飞机撞击下混凝土结构响应的评价时, 通常先基于局部破坏的数值模拟, 验证选用的结构材料参数及失效准则的合理性, 再进行结构的整体响应分析。对于混凝土结构在飞机撞击下的局部破坏, 由于飞机属于硬碰撞物, 局部侵彻需要考虑[2, 10-11]。局部破坏由连续的3个阶段组成[2, 8], 飞射物穿入标靶, 标靶的散裂和碎甲, 飞射物穿透标靶。各状态定义如下:穿入(Penetration)为飞射物进入标靶; 散裂(Spalling)为标靶正前面材料的散裂(飞射物冲击的正面); 碎甲(Scabbing)为标靶后面材料的碎片(冲击面的背面); 穿透(Perforation)为飞射物完全穿透标靶。穿透速度是飞射物恰好完全穿透标靶并且没有穿出时的初始速度。剩余速度是当飞射物的初始速度大于穿透速度时, 飞射物从标靶中穿出的速度。在飞机撞击下混凝土结构响应数值模拟的精确性, 取决于混凝土材料本构方程及其参数选取的正确性, 因此数值模拟中如何选用合适的本构模型以及模型中参数, 已经成为数值模拟的关键。
本文中, 基于ANSYS/LS-DYNA的几种常用的材料参数模型, 进行在飞机撞击下混凝土结构局部破坏的数值分析, 并将数值分析结果与实验结果比较, 给出各个材料参数在模拟飞机撞击混凝土结构的适用性和准确性, 为整体评价飞机撞击下核电厂结构的防御能力提供可信的模拟方法和基础参数。
1. 基于ANSYS/LS-DYNA的数值模拟分析
1.1 实验模型
对于混凝土结构在飞机撞击下的局部破坏评价, 多采用侵彻动力分析方法。由于侵彻动力分析的准确性受多个参数的影响, 有很大的不确定性, 在进行侵彻动力分析中, 常辅以真实的实验过程及实验数据校正数值模拟结果的可靠性。比较典型的飞机撞击混凝土结构实验有, H.Tsubota等[16-18]描述的飞机撞击钢筋混凝土及钢板混凝土结构的实验。为评价飞机撞击混凝土结构的数值模拟结果, 本文中选取比较典型的飞机撞击钢筋混凝土响应的实验。飞机撞击钢筋混凝土实验中, 飞机垂直撞击钢筋混凝土靶板; 钢筋混凝土板的尺寸为1 500 mm×1 500 mm×600 mm, 钢筋间距为D3 mm@25 mm; 取1/7.5全尺寸飞机模型, 飞机模型由4种材料组成, 总重量为245.6 N。表 1为实验中飞机模型中4种材料的主要材料参数和重量; 图 1为实验中飞机撞击钢板混凝土结构的模型剖面示意图。
表 1 飞机模型的材料参数和重量Table 1. The parameters and the weight of aircraft model部件 材料 E/GPa Y/MPa εf W/N 机身外壳 玻璃 6.37 82.3 0.013 81.3 机身填充物 高密度泡沫 2.28 0.6 0.10 129.4 引擎 钢 206 797.4 0.20 28.1 赫氏支撑轴 碳纤维板 0.168 12.9 ∞ 6.8 
图 1 飞机撞击钢筋混凝土的实验模型剖面示意图Figure 1. Schematic section of experimental model of aircraft crashing on the concrete structures1.2 有限元模型
有限元模型建立过程中, 钢筋混凝土的位移协调分离式模型, 即钢筋采用空间梁单元beam 161划分网格, 混凝土采用三维实体solid 164单元, 钢筋与混凝土单元通过Constrained_Lagrange_In_Solid连接, 飞机模型中机身外壳和支撑轴采用壳shell 163单元, 机身填充物和引擎用实体solid 164单元。为保证求解过程的高效及计算结果的收敛性, 混凝土结构模型采用映射法划分网格, 飞机有限元模型均采用与混凝土结构模型一致的网格, 整个有限元模型如图 2所示。对于材料单元中无失效准则的, 通过Mat_Add_Erosion定义失效, 飞机模型中材料的失效准则均与实验模型中提供的参数一致, 采用混凝土材料的失效准则。飞机撞击混凝土结构的接触算法通过Contact_Eroding_Surface_To_Surface和Contact_Surfact_To_Surfact定义。计算过程中不同混凝土材料模型时, 飞机模型与混凝土结构间的接触算法不变, 同时为防止出现沙漏现象, 还加入了Control_Energy和Control_Hourglass。
图 2 飞机撞击混凝土结构的有限元网格划分模型Figure 2. Finite element mesh model of aircraft crashing on reinforced concrete1.3 混凝土本构模型
ANSYS/LS-DYNA广泛应用于爆炸和冲击荷载作用下结构响应的计算, 提供了丰富的混凝土材料模型[14-15], 如:Mat_Pseudo_Tensor(MAT016);Mat_Concrete_Damage_R3(MAT072R3);Mat_Winfrith_Concrete(MAT084);Mat_Johson_Holmquist_Concrete(MAT111);Mat_Soil_Concrete(MAT078);Mat_Brittle_Damage(MAT096);Mat_CSCM_Concrete(MAT159)。每个模型都有自己的优缺点, 需要进行分析才能选到合适的混凝土材料模型。本文中选取4种模型进行飞机撞击混凝土结构的模拟, 用以分析各个混凝土材料参数模拟飞机撞击混凝土结构的反应规律。
1.3.1损伤混凝土材料模型
MAT072R3为损伤混凝土材料模型, 是原72号混凝土材料模型的改进版, 采用3个剪切失效面, 同时考虑损伤和应变率的影响。由于材料模型的参数众多且复杂, 最后又对材料模型添加了参数自动生成功能, 该模型一般只需要混凝土抗压强度等较少的信息, 所以该材料模型适合于模拟混凝土爆炸以及碰撞载荷下的响应。
1.3.2连续帽盖混凝土模型
MAT159为连续帽盖混凝土模型, 该模型针对用于安全设施的混凝土模型, 如高速路旁的混凝土护栏、桥梁护栏等在受到碰撞时, 其变形和破坏能够很好地用该材料模型模拟, 它是地帽模型的扩展, 包括各向同性连续方程、屈服及硬化强度面、破坏方程, 同时也考虑了应变率效应。该材料模型主要用于模拟埋有钢筋的加强混凝土, 可用钢筋、混凝土耦合算法共同模拟钢筋混凝土。
1.3.3混凝土累积损伤模型
MAT111是混凝土累积损伤模型, Johnson-Holmquist-Concrete本构模型及损伤模型是针对混凝土在大应变、高应变率及高围压条件下提出的一种混凝土材料本构关系, 该模型考虑了材料损伤、应变率效应以及静水压力对于屈服应力的影响。它能较好地描述混凝土在高速撞击与侵彻下所产生的损伤、破碎及断裂(或层裂), 在数值模拟中得到了广泛的使用。
1.3.4混凝土Ottosen塑性模型
MAT084是基于剪切破坏提出的混凝土的塑性模型, 主要用于解决冲击荷载作用下钢筋混凝土的响应, 该模型允许每个单元的三向裂纹开展, 而且可在Ls-Prepost查看裂纹。该模型主要应用在模拟冲击荷载作用下混凝土结构的响应。
这些材料模型的理论公式及详细介绍可参见ANSYS/LSDYNA使用手册, 参见文献[14-15]。在第三代核电机组设计中评价飞机撞击混凝土安全壳结构时, 这些模型有很好的应用, 如AP1000核电机组设计中评价屏蔽厂房在大型商用飞机撞击下的响应时, 采用的是混凝土Ottosen塑性模型; EPR核电机组设计中评价混凝土安全壳结构在大型商用分析撞击下的响应时, 采用连续帽盖混凝土模型。
2. 数值分析结果
2.1 计算结果
飞机撞击混凝土钢板的数值计算过程为15 ms。5 ms左右前, 飞机引擎未到达混凝土靶板, 混凝土几乎没有破坏; 5 ms左右时, 飞机引擎开始撞击混凝土靶板; 9 ms左右, 整个飞机模型即将穿透混凝土靶板; 11 ms后, 飞机引擎穿透混凝土结构后飞离混凝土靶板。图 3为9 ms时靶板的正面破坏状态; 图 4为9 ms时靶板的背面的破坏状态; 图 5为9 ms时飞机穿透靶板的侧视图。由图可知, MAT084、MAT159、MAT111和MAT72R3均能模拟飞机撞击混凝土结构的局部效应, 飞机引擎撞击混凝土结构时均出现混凝土的破坏失效和钢筋的断裂的现象; 混凝土结构的背面破坏面积大于正面破坏面积; 相同撞击条件下, 根据混凝土结构破坏面积得到MAT072R3和MAT084两种模型计算结果较接近, MAT111和MAT159两种模型计算结果较接近, 且MAT072R3和MAT084模型比MAT111和MAT159模型破坏面积大。
2.2 与实验结果的比较
为更精确地说明4种混凝土材料模型在模拟飞机撞击混凝土结构的使用性, 图 6为飞机引擎和穿透机身速度的变化过程; 表 2列出了飞机撞击混凝土的数值分析结果与实验结果的对比。表中, dp为正面破坏面积直径, Ab为背面破坏面积。从图 6中可以看到, MAT084、MAT111和MAT72R3变化趋势一致, 穿透机身外壳的剩余速度在110 m/s左右, 引擎的剩余速度在60 m/s左右。表 2列出的实验结果机身外壳的剩余速度在110 m/s, 引擎的剩余速度在80 m/s。MAT084、MAT111和MAT72R3材料模型计算的机身外壳的剩余速度均与实验结果较接近; 对于引擎剩余速度, MAT72R3计算的68m/s较接近于实验结果。
表 2 数值分析结果与实验结果的对比Table 2. Comparison of numerical and the experimental results模型 v0/(m·s-1) vr/(m·s-1) dp/cm Ab/(cm·cm) 机身 引擎 实验 142 110 80 33 50×55 MAT084 142 108 52 30 39×42 MAT159 142 10 5 24 18×27 MAT111 142 113 58 18 30×33 MAT072R3 142 113 68 30 42×42 3. 结论
利用ANSYS/LSDYNA提供的可用于模拟冲击作用下混凝土结构性能的4种材料参数, 对日本Kobori综合研究公司的飞机撞击混凝土结构的实验进行了数值模拟分析, 比较了不同混凝土材料模型在飞机撞击下混凝土结构动力分析时, 混凝土结构靶板的正、背面破坏面积及飞机剩余速度等计算结果。得到:
(1) ANSYS/LSDYNA中混凝土材料模型较多, 各自有着其适用的领域, 在相同分析计算领域中, 不同的混凝土材料模型计算结果也不尽一致。选择合适的混凝土材料模型本构, 有利于提高数值分析计算结果的准确度。
(2) 飞机撞击混凝土结构数值模拟中, 采用的4种混凝土模型参数均能模拟飞机引擎穿入、散裂、碎甲等混凝土结构的局部破坏效应, 但计算的飞机撞击混凝土的正、背面破坏面积及剩余速度不同, MAT072R3和MAT084材料模型的计算结果与实验结果较接近, MAT111材料模型次之, 而MAT159材料模型有较大的差异, 正、背面破坏面积偏小, 剩余速度几乎降低为零。
(3) 由于数值分析计算结果存在一定的不确定性, 所以采用的MAT072R3、MAT084和MAT111混凝土材料模型参数均可为后续评估混凝土结构安全壳抵抗飞机撞击能力时提供基础参数, 也可应用于飞机撞击混凝土结构最不利撞击角度、飞机撞击混凝土结构的穿透公式及飞机撞击混凝土结构的安全壳的整体响应分析等方面的进一步研究中。
-
图 3 汽油蒸气体积分数为2.00%的汽油蒸气/空气预混火焰传播高速摄影图像
Figure 3. High-speed images of gasoline vapor/air premixed flame propagation at the gasoline vapor volume fraction of 2.00%
图 4 汽油蒸气体积分数不同的条件下汽油蒸气/空气预混火焰传播高速摄影图像
Figure 4. High-speed images of gasoline vapor/air premixed flame propagation at different gasoline vapor volume fractions
图 6 拉伸火焰速率与火焰拉伸率之间的关系
Figure 6. Relationship between stretched flame speed and flame stretch rate
-
[1] 刘文辉, 蒋新生, 周建忠, 等.不同环境条件油气爆炸极限测试模拟实验系统[J].后勤工程学院学报, 2013, 29(6):24-25. doi: 10.3969/j.issn.1672-7843.2013.06.005Liu Wenhui, Jiang Xinsheng, Zhou Jianzhong, et al. Simulation experimental system for explosion limit measurement of gasolineair mixture under different environmental conditions[J]. Journal of Logistical Engineering University, 2013, 29(6):24-25. doi: 10.3969/j.issn.1672-7843.2013.06.005 [2] Mannaa O, Mansour M S, Roberts W L, et al. Laminar burning velocities at elevated pressures for gasoline and gasoline surrogates associated with RON[J]. Combustion and Flame, 2015, 162(6):2311-2321. doi: 10.1016/j.combustflame.2015.01.004 [3] Lipatnikov A N, Shy S S, Li W Y. Experimental assessment of various methods of determination of laminar flame speed in experiments with expanding spherical flames with positive Markstein lengths[J]. Combustion and Flame, 2015, 162(7):2840-2854. doi: 10.1016/j.combustflame.2015.04.003 [4] Varea E, Modica V, Renou B, et al. Pressure effects on laminar burning velocities and Markstein lengths for Isooctane-Ethanol-Air mixtures[J]. Proceedings of the Combustion Institute, 2013, 34(1):735-744. doi: 10.1016/j.proci.2012.06.072 [5] Varea E, Modica V, Vandel A, et al. Measurement of laminar burning velocity and Markstein length relative to fresh gases using a new postprocessing procedure: Application to laminar spherical flames for methane, ethanol and isooctane/air mixtures[J]. Combustion and Flame, 2012, 159(2):577-590. doi: 10.1016/j.combustflame.2011.09.002 [6] 赵衡阳.气体与粉尘爆炸原理[M].北京:北京理工大学出版社, 1995:43. [7] Ai Yuhua, Zhou Zhen, Chen Zheng, et al. Laminar flame speed and Markstein length of syngas at normal and elevated pressures and temperatures[J]. Fuel, 2014, 137:339-345. doi: 10.1016/j.fuel.2014.08.004 [8] Bonhomme A, Selle L, Poinsot T. Curvature and confinement effects for flame speed measurements in laminar spherical and cylindrical flames[J]. Combustion and Flame, 2013, 160(7):1208-1214. doi: 10.1016/j.combustflame.2013.02.003 [9] Chao J, Lee J H S, Bauwens C R, et al. Laminar burning velocities of various silanes[J]. Journal of Loss Prevention in the Process Industries, 2015, 36:471-477. doi: 10.1016/j.jlp.2014.11.019 [10] Bosschaart K J, de Goey L P H. The laminar burning velocity of flames propagating in mixtures of hydrocarbons and air measured with the heat flux method[J]. Combustion and Flame, 2004, 136(3):261-269. doi: 10.1016/j.combustflame.2003.10.005 [11] Hayakawa A, Goto T, Mimoto R, et al. Laminar burning velocity and Markstein length of ammonia/air premixed flames at various pressures[J]. Fuel, 2015, 159:98-106. doi: 10.1016/j.fuel.2015.06.070 [12] 齐圣.受限空间油气爆炸与细水雾抑制实验与数值仿真研究[D].重庆: 后勤工程学院, 2014: 107. [13] Clavin P. Dynamic behavior of premixed flame fronts in laminar and turbulent flows[J]. Progress in Energy and Combustion Science, 1985, 11(1):1-59. doi: 10.1016/0360-1285(85)90012-7 [14] Kelley A P, Law C K. Nonlinear effects in the extraction of laminar flame speeds from expanding spherical flames[J]. Combustion and Flame, 2009, 156(9):1844-1851. doi: 10.1016/j.combustflame.2009.04.004 [15] Huang Zuohua, Zhang Yong, Zeng Ke, et al. Measurements of laminar burning velocities for natural gas-hydrogen-air mixtures[J]. Combustion and Flame, 2006, 146(1/2):302-311. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=2bcdfee224dec3ac28a2c4d5bab5510f [16] Liao S Y, Jiang D M, Gao J, et al. Measurements of Markstein numbers and laminar burning velocities for liquefied petroleum gas-air mixtures[J]. Fuel, 2004, 83(10):1281-1288. doi: 10.1016/j.fuel.2003.12.013 [17] Burke M P, Chen Z, Ju Y, et al. Effect of cylindrical confinement on the determination of laminar flame speeds using outwardly propagating flames[J]. Combustion and Flame, 2009, 156(4):771-779. doi: 10.1016/j.combustflame.2009.01.013 [18] Turns S R. An introduction to combustion: Concepts and applications[M]. 3ed. Boston: McGraw-Hill Higher Education, 2012: 218. [19] Zhang Peili, Du Yang, Wu Songlin, et al. Flame regime estimations of gasoline explosion in a tube[J]. Journal of Loss Prevention in the Process Industries, 2015, 33:304-310. doi: 10.1016/j.jlp.2015.01.010 [20] Smallbone A J, Liu W, Law C K, et al. Experimental and modeling study of laminar flame speed and non-premixed counterflow ignition of n-heptane[J]. Proceedings of the Combustion Institute, 2009, 32(1):1245-1252. doi: 10.1016/j.proci.2008.06.213 期刊类型引用(2)
1. 王亚磊, 郑立刚, 于水军, 朱小超, 李刚, 杜德朋, 窦增果. 约束端面对管内甲烷爆炸特性的影响. 爆炸与冲击. 2019(09): 139-148 . 
本站查看2. 郑立刚, 王亚磊, 于水军, 朱小超, 李刚, 杜德朋, 窦增果. NaHCO_3抑制瓦斯爆炸火焰与压力的耦合分析. 化工学报. 2018(09): 4129-4136 . 
百度学术其他类型引用(0)
-
下载:
下载:
百度学术
计量
- 文章访问数: 5097
- HTML全文浏览量: 1475
- PDF下载量: 203
- 被引次数: 2