LS-DYNA软件中4种不同爆炸荷载施加方法的比较

周 清，张云海，齐 麟

(1.滨州市规划设计研究院， 山东 滨州 256600; 2.中国民航大学 机场学院， 天津 300300)

结构抗爆分析通常采用以数值分析为主、试验为辅的研究方法。通用动力学软件LS-DYNA被广泛地应用于抗爆数值分析中，LS-DYNA为用户提供了4种不同的爆炸荷载施加方法。为研究何种爆炸荷载施加方法最能反映真实的情况，首先选用具有代表性的钢筋混凝土板(以下简称RC板)抗爆试验文献作为研究对象，然后利用LS-DYNA软件建立与试验文献完全相同的数值分析模型并采用4种不同的方法施加爆炸荷载，最后通过与文献的试验结果进行比较选择出最合理的爆炸荷载施加方法。

1 LS-DYNA软件中4种不同的爆炸荷载施加方法

1.1 方法1：将爆炸荷载等效为三角形均布荷载直接施加于受爆面

图1 爆炸波冲击压力时程曲线

图2 简化后的三角形压力时程曲线

Pmax=1.407 17/Z+0.553 97/Z2-0.035 72/Z3+

0.000 625/Z4(0.05≤Z≤0.30)

(1)

Pmax=0.619 38/Z-0.032 62/Z2+0.213 24/Z3

(0.30≤Z≤1)

(2)

Pmax=0.0662/Z+0.405/Z2+0.3288/Z3

(1≤Z≤10)

其三，以今香日德为中心的道路：从这里通往境外的道路大体有五条④陈良伟：《丝绸之路河南道》，中国社会科学出版社，2002年，第200-205页。：

(3)

0.033 5Z4)×10-3×W1/3

(4)

1.2 方法2：直接调用LS-DYNA软件中的爆炸超压数据

LS-DYNA中存有由试验得到的TNT空中爆炸超压数据，可以利用关键字按以下方法直接调用。首先采用*LOAD_BLAST_ENHANCED关键字定义TNT当量、爆炸点及爆炸类型，考虑地面反射作用时，还应借助关键字*DEFINE_VETOR定义地面与爆炸点之间的法向向量；然后采用关键字*LOAD_SEGMENT与*LOAD_BLAST_SEGMENT_SET定义爆炸波的接触面。文献[4-6]采用该方法施加爆炸荷载。方法2同样不需要建立空气与炸药的数值分析模型，施加方法简单，计算时间比方法1略长。

1.3 方法3：建立炸药与空气数值分析模型并借助其状态方程施加爆炸荷载

方法3在抗爆数值分析中应用最广。采用方法3的建模方法时，需要对炸药和空气分别单独划分网格并赋予其各自的状态方程，且炸药网格和空气网格应在交界面上共节点。文献[7-10]采用该方法施加爆炸荷载。炸药需要借助关键字*MAT_HIGH_EXPLOSIVE_BURN模型配合Jones-Wilkens-Lee(JWL)状态方程定义[10]，以炸药爆轰产物的压力表示的JWL方程为：

(5)

式(5)中：P1为爆炸压力；V为相对体积；E0为初始内能密度；A、B、R1、R2、ω为状态方程参数。JWL方程参数与TNT材料参数取值如表1所示，表中ρ表示炸药密度，DT表示爆轰速度，Pcj表示爆轰压力。

空气采用*MAT_NULL配合*EOS_LINEAR_POLYNOMIAL状态方程定义[10]，其线性多项式状态方程为：

P2=C0+C1μ+C2μ2+C3μ3+

(6)

式(6)中：μ=ρ/ρ0-1，ρ表示空气的当前密度，ρ0是初始时刻的空气密度；P2为空气压力；C0～C6为多项式方程系数；E1为内能密度；V0为初始相对体积。各相关参数取值见表2。

表1 TNT炸药参数[10]

表2 空气参数[10]

方法3操作比以上方法1和方法2复杂，并且需要消耗较长的计算时间。

1.4 方法4：初始体积分数法施加爆炸荷载

方法4是在方法3基础上改进而得到的一种爆炸荷载施加方法，其施加方式详见文献[11-13]。该方法首先需要定义炸药、空气的材料参数与状态方程；然后借助关键字*ALE_MUTI-MATERIAL_GROUP定义一个包含炸药与空气的材料组；最后借助关键字*INITIAL_VOLUME_FRACTION_GEOMTRY进行填充。以图3为例说明方法3与方法4的区别，图3所建模型为1 m×1 m×1 m的正方体空气与半径0.2 m的球形炸药，当采用方法3建模时需分别单独建立空气与炸药的网格。而采用方法4建模时仅需要建立空气网格作为背景网格，然后利用关键字*INITIAL_VOLUME_FRACTION_GEOMTRY在指定区域将炸药材料填充至建好的空气网格中，并且不要求空气网格和炸药网格共节点。操作完成后，在炸药几何模型所占据的空间中，原网格中的空气被炸药材料所替代。关键字*INITIAL_VOLUME_FRACTION_GEOMTRY的定义如表3所示，共3行。第1行的4个关键字定义初始的空气材料，其中FMSID表示空气的PART号；FMIDITYP表示PART类型；BAMMG表示空气材料在材料组中的编号；NTRACE表示填充采样点数量，可采用默认值。第2行的3个关键字定义炸药材料，其中CNTTYP表示炸药的形状，6为球形；FILLOPT表示填充方式，0表示内部填充；FAMMG表示炸药材料在材料组中的编号。第3行定义炸药的位置及体积，其中X0、Y0、Z0表示炸药球心坐标，R0表示球半径。

图3 方法3与方法4示意图

相对于方法3，方法4的优越性体现在以下2个方面：1) 不需要单独建立炸药单元，仅需在K文件中将炸药进行替换，可以非常方便地定义不同药量的爆炸荷载。2) 可以定义各种形状如圆形、方形、柱形等炸药数值分析模型，避免了不规则炸药体网格划分的困难。

表3 关键字*INITIAL_VOLUME_FRACTION_GEOMTRY的应用

2 试验文献的选择与4种爆炸荷载施加方法的比较

孙文彬[14]对RC板进行了抗爆试验研究，笔者选取该文献作为参考对象。试验采用的2块RC板尺寸均为1 300 mm×1 000 mm×120 mm。板设计为单向板，长边简支，短边自由。采用双层双向配筋，受力的长边配筋为10@100，非受力的短边配筋为10@200。板保护层厚度取值25mm。试验钢筋采用HRB500，其屈服强度为560 MPa，极限强度为605 MPa，弹性模量为230 GPa。采用C40混凝土，抗压强度为48 MPa，弹性模量为48 GPa。试验在爱丁堡国防科技组织的爆炸仓内进行，爆炸点位于RC板中心正上部0.6 m处，爆炸Ⅰ、爆炸Ⅱ的TNT当量为0.079 kg、2.09 kg。爆炸Ⅰ比例距离Z=1.4 m/kg1/3，为了观测RC板在弹性阶段的动力响应。爆炸Ⅱ比例距离Z=0.47 m/kg1/3，为了观测RC板在塑性阶段的爆炸响应和破坏特征。试验布置及采用方法3建立的数值分析模型如图4所示。

图4 试验布置及数值分析模型示意图

结合试验提供的材料数据，采用LS-DYNA的72#材料模型模拟混凝土[15]。当采用国际单位时，72#模型定义C40强度混凝土材料的主要参数如表4所示。其中，MID为材料标识号；ρ为材料密度；A0为初始抗压强度，负值；RSIZE为单位转换数值：m/inches；UCF为单位转换数值：Pa/PSI；LCRATE为所要调用的曲线号。

表4 C40混凝土参数(72#模型)

采用LS-DYNA的3#模型(*MAT_PLASTIC_KINEMATIC)模拟HRB500钢筋[15]，钢筋的主要参数采用表5定义。其中，ρ表示密度；fy为屈服应力；HP为硬化参数；ν为泊松比；E和Etan分别为弹性模型和切线模量；C、P为Cowper-Symonds应变率参数；FS为失效应变。

添加关键字*MAT_ADD_EROSION定义混凝土单元的失效，定义最大抗压强度fc max≤40 MPa与最大抗拉强度ft max≤4 MPa为失效准则，当材料强度大于以上任意数值时，混凝土单元将被删除。钢筋采用3#模型自带的失效准则定义，当钢筋的失效应变大于0.2时，钢筋单元将被删除。炸药与空气定义为多物质欧拉几何实体(EULER)，结构体定义为拉格朗日几何实体(LAGRANGE)，借助于关键字*CONSTRAINED_LAGRANGE_IN_SOILD实现爆炸波与结构体的接触，接触类型可采用无侵蚀的罚函数耦合。

表5 HRB500钢筋参数

2.1 爆炸Ⅰ，4种爆炸荷载施加方法的比较

表6 爆炸I数据汇总

图5 爆炸I，RC板跨中位移时程曲线

通过对表6与图5分析可知，爆炸Ⅰ时：1) 方法1得到的数值分析结果与试验吻合最好，其最大加速度amax、最大位移fmax、最大位移出现时间tmax略大于试验结果。2) 方法2得到的最大加速度amax、最大位移fmax较试验结果偏小，而最大位移出现时间tmax与试验结合相吻合。3) 方法3、方法4得到的各数据均大于试验结果。由于试验结果的RC板在爆炸Ⅰ时最大位移fmax仅为1.12 mm，方法3、方法4比试验结果分别增大0.93 mm、1.53 mm。考虑到试验产生的误差，其计算结果仍然是合理的。(4)除方法4产生少量塑性变形以外，其余3种方法的RC板均处于弹性阶段，与试验结果相一致。以上说明，爆炸荷载较小时4种爆炸荷载施加方法均能较为准确地反映RC板在弹性阶段的动力响应。

2.2 爆炸Ⅱ，4种爆炸荷载施加方法的比较

爆炸Ⅱ的TNT当量为2.09 kg，爆炸点位于RC板中心正上部0.6 m处，比例距离Z=0.47 m/kg1/3。爆炸Ⅱ的RC板发生明显的破坏，图6为RC板顶与底的破坏图。由相关文献可知，板跨中表现为明显的弯曲破坏，支座未出现剪切破坏特征。板顶跨中出现120～150 mm的混凝土压碎区，板底出现2条宽度6～8 mm的主裂缝。板底形成近似正方形的坑，坑边长约为550 mm，周长约为2 200 mm，最大深度约为65 mm。板跨中产生不可恢复的永久变形，最大挠度为50 mm。

图6 爆炸Ⅱ，试验的RC板破坏图

图7 爆炸Ⅱ，4种爆炸荷载施加方法的RC板破坏图

表7 爆炸Ⅱ深坑数据

图8 爆炸Ⅱ，RC板跨中位移时程曲线

3 结论

1) 当爆炸荷载较小时，RC板处于弹性阶段，4种方法均能较为准确地模拟结构构件的动力响应。

2)当爆炸荷载较大时，RC板处于塑性与破坏阶段。方法3可以较为精确地模拟RC板的破坏特征，方法2与方法4可以近似地模拟板的破坏特征，但与试验结果存在一定的误差，不应采用方法1施加爆炸荷载。