我要投搞

标签云

收藏小站

爱尚经典语录、名言、句子、散文、日志、唯美图片

当前位置:斗牛棋牌下载 > 爆速 >

爆炸中应力波理论分析及数值模拟

归档日期:07-25       文本归类:爆速      文章编辑:爱尚语录

  爆炸中应力波理论分析及数值模拟_建筑/土木_工程科技_专业资料。爆炸中应力波理论分析及数值模拟 摘要:利用质量守恒定理以及动量守恒原理,对爆炸过程进行分 析,推导出应力波在爆炸过程中的传播规律:应力波的幅值,波形和 传播速度都会随着介质到重要中心的距离的变化而改变

  爆炸中应力波理论分析及数值模拟 摘要:利用质量守恒定理以及动量守恒原理,对爆炸过程进行分 析,推导出应力波在爆炸过程中的传播规律:应力波的幅值,波形和 传播速度都会随着介质到重要中心的距离的变化而改变, 并且呈现衰 减趋势。并用 ANSYS 模拟球形装药的应力波传播,对上述传播规律进 行说明。 关键词:爆炸 应力波 数值模拟 The Theoretical Analysis and Numerical Modeling of Explosive Stress Wave Abstract:Analysis explosion process with the law of the law of conservation of energy and the law of conservation of mass.Propagation rule of stress wave in the explosion process is deduced.The rule suggests amplitude,waveform and wave velocity all change along with the change of media’s distance to the center of the explosion,and show a trend of attenuation. Simulate stress wave of spherical charge by ANSYS and prove the rule mentioned above. Keywords: Explosion , Stress Wave,Numerical modeling 爆炸时炸药会突然在物理和化学性质上发生巨大变化, 同时伴随 着巨大能量的释放,在爆炸冲击波向外传播是对周围介质进行作用, 所以能够认为是应力波在介质中传播的过程。 随着介质中质点距离爆 炸中心的距离的不同,应力波呈现出不同的特性,在炸药中传播的是 爆轰波,附近介质中为冲击波,随距离增大变为塑性波和弹性波。如 果想要更加深入的认识和理解在工程实际中出现的现象, 以及如何应 对和处理子相关的问题,首先要理解爆炸过程中的冲击波的特性才 行。因为自古以来就有爆炸现象,并且在国防和工农业等各个领域被 广泛应用,所以深受重视。仅仅依靠爆破和抗震等积累经验进行爆破 实践的指导是不充分的, 所以需要对爆破中各种应力波的传播特性进 行分析和模拟。 首先运用质量和动量守恒原理对爆炸过程进行分析, 通过对爆炸 应力波传播规律的推到可以看出, 在介质与装药中心之间的距离不断 增大的过程中,应力波的幅值衰减,波形在变化,传播速度碱小。利 用 LS-DYNA 对球形装药爆炸后应力波传播过程进行模拟分析, 可以较 好的验证这种规律, 所以可以通过模拟对规律进行更加深刻的理解和 认识。 1.应力波特性与距离关系 考察一维应力波(质点运动方向与波的传播方向相同的应力波) 在介质的传播规律,如图 1 所示。 A ? ? d? , ? ? d? , du B C, ? , ? , u A1 B1 图1 一维应力波传播模型 设波速为 C,在各向同性的均匀固体中传播,设某时刻 t 应力波 刚到达 AA1 位置处,假设应力波到达之前介质质点的相关参数为:密 度? 、 应力 ? 、 速度为零。 在经过了 dt 的时间之后, 应力波作用过后, 介质质点的参数变化为:密度 ? ? d? ,应力 ? ? d? ,速度为 du,应力 波的波阵面运动到 BB1。对控制体 AA1BB1 应用质量守恒和动量守恒定 律,因为控制面 AB 和 A1B1 上没有质量交换和动量交换的过程,又因 为在整个过程中压力可以自动达到平衡, 所以能够得到下面的关系表 达式: ?du ? Cd? ?Cdu ? d? (1) (2) 而且应力与体积应变 ? 的关系: d? ? f (?) (3) 由(1)(2)式可得 C2 ? d? / d? (4) 比容 v 与密度的关系式为: v ? 1 ? d? dv ? ? ?2 (5) 将(5)式带入(4)式得 C2 ? d? d? ? ?v 2 2 ? dv dv (6) 又因为 d? ? ? dv v (7) 所以有 C 2 ? v(? v d? 1 d? )d? ? v ? dv d? ? d? (8) (9) 即[1] : C? 1 d? ? d? 2.讨论 图 2 [ 2 ] 给出了相对体积变形和应力波之间的关系曲线。当固体中 的应力 ? 变化的时候,在介质中传播的应力波的幅值、形状和速度都 会随距离的不同而变化。 图2 压应力和相对体积关系 d? 为常数,说明介质处于线) ? ? ? a 时,由图 2 可见 态,应力波为线弹性波,波速 C ? 和泊松比。 E (1 ? ? ) 式中 E, ? 杨氏模量 ? (1 ? ? )(1 ? 2? ) (2) ? a ? ? ? ? b 时从图 2 中的曲线可以看出,在压应力增量很小 的情况下,介质中的质点就会出现相对较大的变形和位移增量,这个 区域就是塑形变形区域,此时材料已经进入屈服阶段,此时介质中传 播的仅仅是弹性纵波,不存在横波。在这个阶段曲线斜率降低很多, 表明了塑性波的传播速度与弹性波的速度相比较有一个较大程度的 降低。 (3) ? b ? ? ? ? c 时,由图 2 曲线可得,此时 d? 增大,但是仍然 d? 小于第一阶段,说明它的传播速度小于弹性波速(声速),但是大于 塑性波的波速。然而压应力值却大于弹性压应力极限,所以此时传播 的就是亚声速冲击波。 (4) ? ? ? c 时,图 2 中的曲线变得更加陡峭,在这个阶段 d? 比 d? 阶段(3)大很多,说明它的波速比弹性波更大,压应力的最大值已 经超过了介质的塑性极限,此时就是超声速冲击波。 通过以上分析可得,应力波的传播速度和波形等是由固体中的 压应力决定的。炸药介质内部 压应力与爆炸压力相接近,此 (R ? Rc) 时传播的应力波为冲击波,波形十分陡峭,其传播速度大于声速;当 距离增加到 R c ? R ? Rb 时,冲击波的传播速度以及压应力的最大值都 有所减小,此时为亚声速传播的冲击波区;当 R b ? R ? Ra 时,应力波 进一步衰减成弹塑性应力波的双波结构, 切塑性波的传播速度小于弹 性波的传播速度;当 R a ? R 时,应力波进一步衰减,塑性波消失,仅 剩下弹性波在介质中传播。应力波形图如图 3 [ 3 ] 所示。 图 3 典型应力波形图 传播的数值模拟 3. 无限介质中爆炸应力波 用 LS-DYNA 模拟球形装药在无限介质中的爆炸, 以此来验证上述 应力波的传播规律。 因为球体具有对称性, 所以可以建立实体的 18 模 型,在进行单元网格划分时采用轴对称形式 [ 4 ] ,模型如图 4、图 5 所 示。 图 4 装药中心示意图 图 5 网格划分模型图 LS-DYNA 可以对爆炸过程进行计算,而且可以对爆炸与介质之间 的相互作用进行耦合计算,对于计算爆炸地震波的生成传播,以及它 对结构的作用是一种较为可靠的方法。 将水泥作为传播介质,其密度为 2.13g cm3 ,屈服强度 39.8Mpa,弹 m ? s ?1 。 性波传播速度 3233 建立半径为 100cm 的模型进行计算, 由理论 计算可得弹性波到达边界的时间约为 309?s ,取计算时间为 300?s ,, 忽略波的反射效应。 图6 单元应力— 时间曲线 节点应力— 时间曲线 所示。计算结果显示:距离中心最近的单元 和节点(图)的应力和压力是最大的,并且远大于材料的屈服极限, 随着时间的变化演变成冲击波形式;在离开爆源更远一些的位置 (图),尽管介质的应力和压力已经变得较为平缓,但它们的最大值 依然大于介质的弹性极限(图);在远离爆炸中心的地方(图),应 力波变为弹性波。这些模拟结果与理论分析的结果吻合的较好,说明 模拟手段可以有效的预见爆炸中与应力波传播有关的现象。 参考文献: [1]王礼立,朱兆祥.应力波基础(第二版).国防工业出版社,1991 [2]熊建国.爆炸动力学及其应用.科学出版社,1987 [3]郭伟国.应力波基础简明教程.西北工业大学出版社,2007.4 [4]时党勇,李裕春,张胜民。基于 ANSYS/LS-DYNA8.1 进行显式动力 分析.清华大学出版社,2005.

本文链接:http://azizburney.com/baosu/358.html