Vibration intensity of multi milliseconds blasting in underground mine based on LS-DYNA
-
摘要: 某地下矿山上方地表分布有土坯房,一层、二层砖混结构民房以及水稻田.该矿山生产过程中产生的爆破震动易影响周边民房的稳定.为研究该矿山井下生产爆破时地表的震动强度,确保民房震动强度安全,利用LS-DYNA有限元方法对某次具有代表性的生产爆破下的地表震动强度进行了分析研究.通过对该次具有代表性的生产爆破数值计算与爆破震动现场测试结果对比分析表明,采用ANSYS/LS-DYNA的数值计算结果与实际测试值偏差不大,数值模拟分析方法对爆破震动强度的预测具有一定的指导价值.
-
关键词:
- ANSYS/LS-DYNA /
- 多段毫秒延时爆破 /
- 爆破震动
Abstract: There are adobe houses, one or two layers of brick housing and paddy field over the underground mine. The stability of the houses is influenced by blasting vibration generated during the production process of the mines. In order to predict vibration intensity and ensure the safety of surface buildings in an underground mine production multi millisecond blasting, a numerical simulation analysis of the vibration of a certain representative production blasting is conducted by LS-DYNA software. The comparative analysis of the numerical results and the actual measurements shows that the result of ANSYS/LS-DYNA software element software has little difference with the actual test values. ANSYS/LS-DYNA software has a certain guiding value to forecast blasting vibration intensity.-
Keywords:
- ANSYS/LS-DYNA /
- multi milliseconds blasting /
- blasting vibration
-
0 引言
多头螺旋管作为近几年来兴起的新型高效换热元件,已然成为传热行业研究的热点,相比对称型良好的传统换热管材,如光管、翅片管等,螺旋管内的流体在进行管内外对流换热的过程中,涉及一个更为复杂的旋转紊流的三维流动[1].考虑到多头螺旋管复杂、特殊的结构及换热过程流体复杂的流动状态,文中选用数值模拟的方法对其换热过程进行研究,通过模拟计算,得到螺旋管内流体紊流流动过程中的温度场与流场分布,为其强化换热的研究提供理论依据.
1 换热过程的有限元模型
实体模型是建立有限元模型的基础,考虑到多头螺旋管非轴对称的结构特性,文中选用了三维实体建模软件Pro/E 来建立多头螺旋管的模型,然后以特定格式导入到ANSYS 中进行数值模拟分析.表 1示出了实验管材的基本结构参数.
表 1 多头螺旋管结构参数根据表 1 给出的管材结构参数,建立多头螺旋管的实体模型,选用当量直径为内径20 mm, 管长为200 mm 的光管,由于文中研究的是管内紊流流体的流场及温度场变化,所以在建立模型时只需建立流体的实体模型,这使得模型得到了科学的简化.图 1及图 2 分别为螺旋管换热过程的几何模型及有限元模型[2].
2 边界条件及载荷施加
固体壁面处确定为无滑移条件,为了简化计算模型,假定换热过程,液体的物性参数为常值,管外的高温空气以恒定的壁面热流率形式施加[3-5],具体的载荷及流体物性参数如表 2 所示.
表 2 载荷及流体物性参数3 模拟结果及分析
3.1 温度场结果分析
图 3及图 4 分别示出了光管及六头螺旋管换热过程中的温度分布,从图 3、图 4 不难看出,低温的水在流动过程中,温度逐渐升高,越靠近管壁,温度越高,但由于流体存在轴向的速度,各个部分的温度并不相同,沿着流动方向温度分布呈现峰状,低温流体区域逐步递减,并且相对光管,螺旋管的峰形更尖锐,这也从侧面说明螺旋管内流体的换热更剧烈,效率更高.临近出口端,流体的温度分布基本趋向均匀.
图 5和图 6 将温度分布以数据的形式反映在数学曲线上,其中横坐标代表从入口端到出口端的管的轴向长度,纵坐标代表对应管长截面上的温度值.设温升的相对提高为ε,Tg为光管中流体从入口端到出口端的温升,Tl为螺旋管内流体入口端与出口端的温差[6].则根据曲线中的数据不难得到:
相对光管,螺旋管内流体的温升相对提高了近2.6 倍.这从客观数据上证明了螺旋管的强化传热作用.
3.2 流场结果分析
图 7及图 8 表征了光管和螺旋管中流体运动过程的速度分布,从图 7、图 8 可以看出,沿管材的径向呈现速度由近壁处向中心逐渐增大,且近壁处速度升幅剧烈,中间部位速度变化较为平缓,速度几近峰值.相对光管,螺旋管内紊流的流体,紊流的无序特征变得明显,特别是在近壁处,流体已出现明显的流体回旋运动,即所谓的有涡流动[7],正是由于这种三维有涡流动的存在,使得螺旋管内壁不易结垢,从而起到自除垢的作用.
图 9和图 10 以数值曲线的方式反映了出口流体沿径向上合速度的变化,不难发现无论光管还是螺旋管,其管内流体的合速度分布均沿管的轴向成中心对称,相比之下,螺旋管内流体近壁处速度变化更剧烈.根据曲线上的速度值可得出光管出口流体的平均速度约为0.6 m/s, 螺旋管出口流体的平均速度约为0.3 m/s.这主要是因为螺旋管旋转的凸筋结构增加了流体流动的阻力,但是在流体流动过程中,螺旋管正是由于其具有这种特殊的螺旋凸筋结构,破坏了流体传热的边界层,才使得其强化传热强度得到了提高[7-14],同时,根据过增元院士的场协同理论可知,速度矢量与热流矢量的夹角越小,速度场与热流场的协同性越好,则其强化传热性能越好.相对光管,螺旋管内流体靠近管壁处其速度表现为轴向、径向和周向上三维流动,这种有涡流动大大改善了速度场和热流场的协同效应,从而起到明显的强化传热作用[15-19].
4 结论
(1)在等热流的情况下,多头螺旋管的换热性能优于光管,以文中研究的六头螺旋管为例,其流体出口温升较之基准管相对提高近2.6 倍.
(2)流体流动过程中,流体流速的变化与到管壁的距离成反比,近壁处速度变化剧烈.相比光管内流体,多头螺旋管内流体流动阻力更大,其出口速度相对降低了50 %.
(3)多头螺旋管特有凸筋结构的存在,一方面破坏了流体传热的边界层,并且改善了速度场与热流场的协同性,使得强化传热效率得到提高,另一方面,有涡流动对管壁的冲刷作用,使得多头螺旋管管壁不易结垢.
-
表 1 各测点实测振速峰值
测点 斜距/m 最大合成速度/(cm·s-1) 三维方向及其对应速度/(cm·s-1) 主震频率/Hz 测点1 84.91 1.136 T 1.11 35.5 V 1.04 42.0 L 9.93 15.8 测点2 171.93 0.291 T 0.190 8.63 V 0.245 10.9 L 0.241 9.25 测点3 258.61 0.116 T 0.063 5 10.4 V 0.010 2 30.7 L 0.010 2 27.3 测点4 285.31 0.282 T 0.216 27.0 V 0.216 20.0 L 0.279 37.0 测点5 298.40 振速小于0.051 cm/s, 没有触发测振仪 表 2 花岗岩主要物理力学参数表
密度/(kg·m-3) 单轴抗压强度/MPa 泊松比 弹性模量/MPa 纵波速/(m·s-1) 2 566 142 0.3 4.66×104 4 500 表 3 炸药主要参数表
密度/(kg·m-3) 爆速/(m·s-1) A B R1 R2 1 000 3 600 2.140×1011 1.82×109 4.15 0.95 表 4 土体主要参数表
密度/(kg·m-3) G K a0 a1 a2 Pc 1.8×103 1.601×103 1.328×1010 3.3×109 1.31×104 0.123 2 0.0 表 5 各测点模拟计算振速峰值
测点 斜距/m 最大合成速度/(m·s-1) 三维分速度方向 三维方向速度/(m·s-1) 测点1 84.91 2.70 T 2.54 V 0.861 L 0.308 测点2 171.93 0.611 T 0.492 V 0.264 L 0.351 测点3 258.61 0.492 T 0.394 V 0.346 L 0.462 测点4 285.31 0.329 T 0.283 V 0.301 L 0.317 测点5 298.40 0.091 T 0.012 8 V 0.010 7 L 0.069 -
[1] 石少卿, 康建功, 江敏, 等.L S-DYNA在爆炸与冲击领域内的工程应用[M].北京:中国建筑工业出版社, 2011. [2] 张智超, 刘汉龙, 陈育民, 等.爆破地震的数值模拟及爆破振动规律分析[J].郑州大学学报(工学版), 2012, 33(5):10-15. http://www.cnki.com.cn/Article/CJFDTOTAL-ZZGY201205004.htm [3] 夏祥, 李海波, 李俊如, 等.岩体爆生裂纹的数值模拟[J].岩土力学, 2006, 27 (11): 1987-1991. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200611026.htm [4] 梁斌, 陈忠富, 卢永刚.柱形装药在混凝土中爆炸波传播数值模拟[J].弹箭与制导学报, 2008, 28(4):93-95. http://www.cnki.com.cn/Article/CJFDTOTAL-DJZD200804027.htm [5] 肖明, 张雨霆, 陈俊涛, 等.地下洞室开挖爆破围岩松动圈的数值分析计算[J].岩土力学, 2010, 31(8): 2613-2618. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201008044.htm [6] 杨海书, 林从谋, 林丽群, 等.复杂结构体系下隧道爆破振动对房屋影响的试验研究[J].山东科技大学学报(自然科学版), 2011(2):65-79. http://www.cnki.com.cn/Article/CJFDTOTAL-SDKY201102010.htm [7] 张智超, 陈育民, 刘汉龙.微差爆破模拟天然地震的数值分析与效果评价[J].岩土力学, 2013(1): 265-274. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201301039.htm [8] 李宁, 张承客, 周钟.边坡爆破开挖对邻近已有洞室影响研究[J].岩石力学与工程学报, 2012(增刊2): 3471-3476. http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2012S2004.htm [9] 肖文芳, 肖文涛, 房泽法.地铁隧道掘进爆破地表振动效应研究[J].工程爆破, 2012, 18(2):46-49. http://www.cnki.com.cn/Article/CJFDTOTAL-GCBP201202014.htm [10] 史秀志, 曾志林, 田建军, 等.深井开采爆破对巷道影响的数值模拟分析[J].工程爆破, 2010, 16(2):10-14. http://www.cnki.com.cn/Article/CJFDTOTAL-GCBP201002004.htm [11] 时党勇, 李裕春, 张胜民.基于ANSYS/LS-DYNA8.1进行显式动力分析[M].北京:清华大学出版社, 2005. [12] 白金泽. LS-DYNA3D理论基础与实例分析[M].北京:科学出版社, 2005. [13] 何涛. ANSYS10.0/LS-DYNA非线性有限元分析实例指导教程[M].北京:机械工业出版社, 2007. [14] 尚晓江, 苏建宇, 魏久安, 等.ANSYS/LS-DYNA动力分析方法与工程实例[M].北京:北京水利水电出版社, 2006. [15] Schneidersr. A-grid algorithm for the generation of hexahedral element meshes[J].Engineering with computers, 1996(12):168-177. doi: 10.1007/BF01198732