毕程程;王志亮;郝士云;石恒;石高扬;李鸿儒;何爱林;
BI Cheng-cheng;WANG Zhi-liang;HAO Shi-yun;SHI Heng;SHI Gao-yang;LI Hong-ru;HE Ai-lin;School of Civil and Hydraulic Engineering,Hefei University of Technology;

• 1:合肥工业大学土木与水利工程学院

摘要(Abstract)：

为提高LS-DYNA软件模拟工程爆破时的精度,从理论上改进了确定极限时间步长Δte的方法,结合已有的数据对本方法的合理性进行了验证,并通过实例探讨了时间步长Δt对爆破模拟结果的影响。研究表明:Δte仅与炸药的爆速、炸药最小单元的最小边长以及计算时间步长比例因子有关,在网格划分时,炸药单元不能取过小或过于狭长;Δt对土体中爆炸空腔半径形成的影响不大,而对压力峰值的影响显著,且具有一定的规律性;随着比例距离Z的增加,压力峰值误差呈现由小到大,再由大到小的趋势。当计算区域范围在Z<0.5m·kg^(-1/3)时,Δt可取为1(-1/3)时,Δt可取为1³倍Δte;当Z≥0.5m·kg3倍Δte;当Z≥0.5m·kg^(-1/3)时,时间步长Δt可适当放宽,取为3(-1/3)时,时间步长Δt可适当放宽,取为3⁵倍Δte。对于一般的爆破模拟,考虑到经济性,Δt取为15倍Δte。对于一般的爆破模拟,考虑到经济性,Δt取为1³倍Δte即可满足大多数工程计算精度的要求。
In order to improve the precision of simulation blasting with LS-DYNA software,the method of determining limit time stepΔteis improved forward theoretically,the rationality of the method is verified using the existing data,and the influence of time stepΔt on blasting simulation results is discussed through engineering examples.Results shows thatΔteis only associated with the detonation velocity of explosive,the minimum side length of explosive unit and the scale factor of calculation time step.The explosive element cannot be too small or too narrow when the meshing grid is made.Δt has little effect on the formation of the cavity radius in the soil,but it exerts a significant effect on the peak pressure and has a certain regularity.With increasing the proportional distance Z,the pressure peak error increases firstly and then decreases.When the calculating scope is in the range of 03倍Δte即可满足大多数工程计算精度的要求。
In order to improve the precision of simulation blasting with LS-DYNA software,the method of determining limit time stepΔteis improved forward theoretically,the rationality of the method is verified using the existing data,and the influence of time stepΔt on blasting simulation results is discussed through engineering examples.Results shows thatΔteis only associated with the detonation velocity of explosive,the minimum side length of explosive unit and the scale factor of calculation time step.The explosive element cannot be too small or too narrow when the meshing grid is made.Δt has little effect on the formation of the cavity radius in the soil,but it exerts a significant effect on the peak pressure and has a certain regularity.With increasing the proportional distance Z,the pressure peak error increases firstly and then decreases.When the calculating scope is in the range of 0⁰.5 m·kg0.5 m·kg^(-1/3),Δt can be taken as 1(-1/3),Δt can be taken as 1³ timesΔte.When it exceeds 0.5 m·kg3 timesΔte.When it exceeds 0.5 m·kg^(-1/3),Δt can be appropriately increased,taking 3(-1/3),Δt can be appropriately increased,taking 3⁵ timesΔte.For the general blasting simulation,by taking into account the economy,Δt is taken as 15 timesΔte.For the general blasting simulation,by taking into account the economy,Δt is taken as 1³ timesΔte.And it may meet the requirements of calculation accuracy for most engineering practices.

关键词(KeyWords)： LS-DYNA;爆破;时间步长;爆速;最小单元
LS-DYNA;blasting;time step;detonation velocity;minimum element

Abstract:

Keywords:

基金项目(Foundation): 国家自然科学基金资助项目(51379147;51579062)

作者(Authors): 毕程程;王志亮;郝士云;石恒;石高扬;李鸿儒;何爱林;
BI Cheng-cheng;WANG Zhi-liang;HAO Shi-yun;SHI Heng;SHI Gao-yang;LI Hong-ru;HE Ai-lin;School of Civil and Hydraulic Engineering,Hefei University of Technology;

参考文献(References)：

• [1]李利莎,谢清粮,郑全平,等.基于Lagrange、ALE和SPH算法的接触爆炸模拟计算[J].爆破,2011,28(1):18-22,27.LI L S,XIE Q L,ZHENG Q P,et al.Numerical simulation of contact Explosion based on lagrange ALE and SPH[J].Blasting,2011,28(1):18-22,27.
• [2]KONESHWARAN S,THAMBIRATNAM D P,GALLAGE C.Blast response of segmented bored tunnel using coupled SPH-FE method[J].Structures,2015(2):58-71.
• [3]石磊,杜修力,樊鑫.爆炸冲击波数值计算网格划分方法研究[J].北京工业大学学报,2010,36(11):1 465-1 470.SHI L,DU X L,FAN X.A study on the mesh generation method for numerical simulation of blast wave[J].Journal of Beijing University of Technology,2010,36(11):1 465-1 470.
• [4]赵蓓蕾,崔村燕,陈景鹏,等.网格密度对基于LSDYNA的冲击波数值模拟影响[J].自动化与仪器仪表,2015(9):227-228,231.ZHAO B L,CUI C Y,CHEN J P,et al.The effect of the mesh density on numerical simulation of the shock wave based on LS-DYNA[J].Automation and Instrumentation,2015(9):227-228,231.
• [5]高轩能,吴彦捷.TNT爆炸的数值计算及其影响因素[J].火炸药学报,2015,38(3):32-39.GAO X N,WU Y J.Numerical calculation and influence parameters for TNT explosion[J].Chinese Journal of Explosives and Propellants,2015,38(3):32-39.
• [6]ARTERO GUERRERO J,PERNAS SANCHEZ J,TEIXEIRA DIAS F.Blast wave dynamics:The influence of the shape of the explosive[J].Journal of Hazardous Materials,2017,331(5):189-199.
• [7]王建灵,郭炜,冯晓军.TNT、PBX和Hexel空中爆炸冲击波参数的实验研究[J].火炸药学报,2008,31(6):42-44.WANG J L,GUO W,FENG X J.Experimental research on the air explosion shock wave parameters of TNT、PBX and hexel[J].Chinese Journal of Explosives and Propellants,2008,31(6):42-44.
• [8]宋浦,杨凯,梁安定,等.国内外TNT炸药的JWL状态方程及其能量释放差异分析[J].火炸药学报,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 and Propellants,2013,36(2):42-45.
• [9]尚晓江,苏建宇,王化锋,等.ANSYS/LS-DYNA动力分析方法与工程实例[M].北京:中国水利水电出版社,2008.SHANG X J,SU J Y,WANG H F,et al.ANSYS/LS-DYNA dynamic analysis method and engineering example[M].Beijng:China Water and Power Press,2008.
• [10]张奇,张若京.ALE方法在爆炸数值模拟中的应用[J].力学季刊,2005,26(4):639-642.ZHANG Q,ZHANG R J.Numerical simulation of explosion using ALE method[J].Chinese Quarterly of Mechanics,2005,26(4):639-642.
• [11]郝士云,毕程程,王志亮.LS-DYNA模拟爆破中时间步长对计算结果影响探讨[J].爆破,2016,33(4):39-45.HAO S Y,BI C C,WANG Z L.Study on influence of time step on calculation result in blasting simulation with LS-DYNA[J].Blasting,2016,33(4):39-45.
• [12]WANG Z L,LI Y C,WANG J G.Numerical analysis of blast-induced wave propagation and spalling damage in a rock plate[J].International Journal of Rock Mechanics and Mining Sciences,2008,45(4):600-608.
• [13]LSTC.LS-Dyna keyword user's manual version 971[M].Livermore:Livermore Software Technology Corporation,2007.
• [14]时党勇,李裕春,张胜民.基于ANSYS/LS-DYNA 8.1进行显示动力分析[M].北京:清华大学出版社,2005.SHI D Y,LI Y C,ZHANG S M.The dynamic analysis based on ANSYS/LS-DYNA 8.1[M].Beijng:Tsinghua University Press,2005.
• [15]李夕兵,古德生,赖海辉,等.岩石与炸药波阻抗匹配的能量研究[J].中南矿冶学院学报,1992,23(1):18-23.LI X B,GU D S,LAI H H,et al.The energy analysis on matching of acoustic impedance between rock and explosive[J].Journal of Zhongnan Institute of Mining and Metallurgy,1992,23(1):18-23.
• [16]MAHMADI K,AQUELET N.Euler-Lagrange simulation of high pressure shock waves[J].Wave Motion,2015(54):28-42.
• [17]潘强,张继春,郭学彬.土体爆炸压密的数值模拟研究[J].矿业研究与开发,2011,31(2):101-104.PANG Q,ZHANG J C,GUO X B.Numerical simulation research on explosion compaction of soil[J].Mining Research and Development,2011,31(2):101-104.
• [18]王清洁,顾文彬,茼茂辉,等.半无限土介质中集团装药爆炸空腔的数值模拟[J].爆破,2002,19(3):17-19.WANG Q J,GU W B,MAI M H,et al.Numerical simulation of cylindrical charges explosion in semi-infinite soil medium[J].Blasting,2002,19(3):17-19.
• [19]顾文彬,叶序双,詹发民,等.球形装药半无限土介质中爆炸动力学分析[J].工程爆破,1999,5(1):5-10.GU W B,YE X S,ZHAN F M,et al.Dynamic analysis on spherical charges exploding in semi-infinite soil medium[J].Engineering Blasting,1999,5(1):5-10.
• [20]高金石,张奇.爆破理论与爆破优化[M].西安:西安地图出版社,1993.GAO J S,ZHANG Q Z.Blasting theory and blasting optimization[M].Xi'an:Xi'an Cartographic Publishing House,1993.