爆炸冲击波对建筑玻璃的影响.docxVIP

爆炸冲击波对建筑玻璃的影响 建筑玻璃的耐爆性已成为人们关注的中心主题。纯理论法难以量化和分析，实验研究成本高、风险。近年用数值模拟对玻璃抗爆性能的研究取得了一定进展。M.Larcher等建立了分层壳单元的夹层玻璃模型,定义了PVB胶的弹塑性失效准则,并通过实验验证了模型的正确性。T.Krauthammer等评估了爆炸冲击波负压阶段对玻璃板的影响。邓荣兵等通过LS-DYNA软件,采用多物质的ALE有限元法,研究了玻璃幕墙在爆炸冲击波作用下的动态响应过程,分析了玻璃幕墙的破坏情况。上述研究中,玻璃模型只是通过线弹性材料模型来模拟,并没有建立不同种类玻璃在爆炸冲击波作用下的计算模型,对玻璃的失效准则也只是定义了拉应力失效。 本文中,建立3种常用建筑玻璃在近场爆轰作用下的有限元分析模型,以冲击波超压对玻璃面板作用产生的拉应力和剪切应变为失效判据,并通过实验对模型进行对比验证。空气冲击波超压有自由场超压和壁面超压2种。壁面超压是指初始冲击波到达壁面(主要是地面)产生反射后壁面所承受的压力,该压力比自由场压力高。由于爆炸的破坏作用等原因,要测定自由场超压很困难,故本文中的超压值为壁面超压值。 1 空气波浪传播规律模拟 选用ANSYS/LS-DYNA作为计算工具。采用LS-DYNA对空气冲击波传播规律的模拟包括TNT炸药的起爆、爆轰波及空气冲击波的形成及传播、冲击波的相互作用等复杂的物理化学过程,几何模型如图1所示。 1.1 数值模拟方法 建立的几何模型主要包括TNT炸药、空气和玻璃。考虑到问题的对称性,为减少计算量,取1/2实体进行建模。其中,炸药采用柱形结构;空气场为长1.2 m、宽为0.25m与玻璃同宽的矩形结构。模型中对称面采用对称边界,非对称面采用透射边界,使用刚性壁面来模拟地面。模拟整个爆炸过程的数值方法主要包括Lagrange法、Euler法、ALE(arbitrary Lagrange-Euler)法、SPH(smoothed particle hydrodynamics)法。本文中采用ALE法实现炸药与空气场的耦合作用。首先,在结构边界运动的处理上,引进了Lagrange法的特点,能有效跟踪物质结构边界的运动;其次,在内部网格的划分上,吸收了Euler法的长处,使内部网格单元独立于物质实体,可以根据定义的参数,在求解过程中适当调整网格位置,使网格不出现严重畸变。 1.2 空气场、范围的耦合 模型中全部采用SOLID164八节点实体单元。为了模拟玻璃的裂纹扩展情况,用Lagrange单元描述玻璃,用多物质ALE单元描述空气和TNT炸药,采用Van Leer(二阶精度)方法进行计算。空气场与玻璃相连接的一层网格为ALE网格到Lagrange网格的过度层网格,实为Lagrange网格。当载荷过大时,会发生网格畸变,导致计算无法正常进行。本模拟中采用当过渡层网格变形在可控范围内并且冲击波正压段已作用到玻璃时,通过重启动删除炸药和空气场网格的方法保证计算正常进行,并且能提高计算效率,计算结果也在可接受范围。模拟中,通过命令*ALE_MUITI-MATERIAL_GROUP和*SECTION_SOLID_ALE定义空气和炸药的耦合及其算法;通过命令*CONSTRAINED_LAGRANGE_IN_SOLID_TITLE定义空气与玻璃面板之间的耦合作用。 1.3 材料模型和形状方程 1.3.1 爆爆过程状态方程 炸药为60g TNT,采用*HIGH_EXPLOSIVE_BURN模型作为炸药的材料模型,采用JWL状态方程描述TNT炸药,通过命令*INITIAL_DETONATION来控制炸药起爆的时间和爆心。JWL状态方程的表达式为: 式中:p为爆轰压力;V为爆轰产物体积相对于初始体积的比值;e0为初始体积内能;ω、A、B、R1和R2均为方程系数。该状态方程通过反应爆炸气体压力-体积关系来描述炸药的爆轰过程。 具体参数值为:炸药密度ρ=1.631g/cm3,爆速D=6.717km/s,爆轰波阵面压力pCJ=18.5GPa,e0=7GPa,V0=1.0,A=540.9GPa,B=9.373GPa,R1=4.5,R2=1.1,ω=0.35。 1.3.2 线性多项式状态方程 对空气采用*MAT_NULL材料模型以及线性多项式*EOS_LINEAR_POLYNOMIAL状态方程加以描述。线性多项式状态方程为: 式中:C0~C6为常量,e为初始体积内能,μ=1/V0-1,V0为相对体积。具体参数值为:空气密度ρ=0.001 3g/cm3,C0=-0.1 MPa,C1=C2=C3=0,C4=C5=0.4,C6=0,E=0.25 MPa,μ=1。 1.3.3 材料模型和强度模型 对玻璃面板采用框支约束,对玻璃材料四周采用固端约束。通过命令*MAT_A