本文通过采用三维动力有限元法和激光全息动光弹试验对半无限介质中柱状药包的应力场、台阶爆破应力场和若干爆破设计参数进行了初步数值分析研究。 首先阐述了三维动力有限元DYNA3D系统的数值计算方法，对有限元计算中的几个关键处理：沙漏控制、应力波与人工体积粘性、接触碰撞仅滑移（Sliding-only）边界、无反射边界等做了说明。为了应用三维动力有限元DYNA3D系统进行模拟计算，本文依据相似原理，采用1:10的几何相似比进行计算模拟。由于单值条件与物理相似条件在模拟上存在具体困难，因此完全满足现象相似的充要条件不太可能。在本文计算中主要应用了几何相似、边界条件相似和部分动力相似条件，在本论文中取C_μ=1，C_ι=10，C_ρ=1，C_Ε=1，C_ε=1，C_σ=1，C_ι=C_ι=10。由于岩体在动载作用下的破坏特征以及塑性特征的研究仍处在十分初浅的百家争鸣的阶段，尚未上升到数学描述阶段，故岩石在爆炸作用下的强动载荷作用特性还不是很清楚，为简化计算，本章将岩石材料假设为线弹性体，采用Mises屈服条件来描述岩石材料的破坏特性。 三维动力有限元数值分析和全息动光弹试验结果表明：半无限介质中柱状药包爆炸应力场在底部点火起爆后，应力波呈水滴状向外扩展，其波阵面与炮孔轴线所形成的角度基本不变，直到爆轰结束。此后有效应力场逐渐呈近似椭球体向四周传播并开始衰减。半无限介质中柱状药包的端部应力场与水平条形药包的应力场都存在着有效应力为零的“应力空洞”区，但是，条形药包端部应力场两端是对称分布的，而柱状药包的上下端部应力场是不对称的，且柱状药包底部零畸变范围要比顶部大，这个现象同样被动光弹试验所证实。双自由面的单孔台阶爆破数值计算表明，在应力波传播的初期，初始应力场的发展规律是和半无限体 铁道部科学研究院博士论文2001 中柱状药包爆炸应力场的发展规律是相似的，在有效应力场的传播过程 中，两端部的无畸变区由炮孔中心线处向最小抵抗线方向偏移。无超深 爆破时，由于底部受到柱状药包端部效应的影响，在炮孔附近距轴线相 等距离上，中部有效应力要比底部大1一2倍。 对深孔台阶爆破若干设计参数（炸药单耗、超深、填塞长度）的数 值计算表明，炸药单耗的变化对应力场有显著的影响，单耗增加42％时， 在抵抗线方向的有效应力相应增加1 ．57一3．1倍。超深改变对台阶爆破时 底部应力场的变化影响很大，台阶爆破最大超深约为0．3倍最小抵抗线 长度。这个结论与深孔台阶爆破的工程实践相吻合。数值计算验证了填 塞长度对深孔台阶爆破破碎效果的影响，虽然无填塞爆破对孔底抵抗线 位置的破坏应力基本没有影响，但对于台阶上部岩体中的应力场影响较 大，在相同条件下无填塞爆破填塞区中有效应力减小约60％，说明无填 塞不利于孔口区的破碎，从理论上证实了无填塞爆破对爆破能量是相当 大的浪费，并且得到最佳填塞长度Ls＝（0一1．0）万，这个结论与工程 实际经验也是相符的。