圆柱形容器中心内爆的流固耦合计算.docVIP

圆柱形容器中心内爆的流固耦合计算 张亚军　　 徐胜利 中国科学技术大学力学和机械工程系， 合肥 230026 摘 要 针对中心装药的冲击波加载，本文研究椭球封头圆柱壳的流固耦合问题。基于法向位移、速度和力建立流固界面边界条件。采用FVM求解积分ALE型无粘流体力学方程组，采用FEM求解虚功原理，得到内爆流场冲击波和壳体应力波的传播过程。结果表明：本文方法和计算软件可研究内爆引起的壳体流固耦合问题。尽管壳体结构和装药位置相对简单，但冲击波在壳内传播和内壁面多次反射产生的波系结构相当复杂，压力云图刻画了冲击波演化图像。壳体所受加载特征是：前期表现为冲击波超压，后期表现为冲量，壳体变形依赖冲击波加载历史。应力波在壳体内外壁面反射和透射、并沿壳体传播。封头顶点、封头和侧壁交线以及距中心最近的中环面出现较大等效应力，也是结构易破坏的部位。 关键词　流固耦合; 圆柱壳; FVM; FEM; 三维数值计算 引言 单层和多层爆炸容器在爆炸与冲击实验中有着广泛应用。爆炸冲击加载的圆柱壳动力学响应研究实质上是爆炸冲击载荷作用下的流固耦合问题。典型结构是带封头的单层圆柱壳，特别是小药量炸药爆轰情况。耦合计算需采用Euler方法描述流体、Lagrange方法描述固体，这比高压下固体采用流体弹塑性模型计算困难。Redekop在螺旋坐标系下计算了钢螺旋板受冲击载荷作用的动态响应[1]。采用非线性有限元，Ewing等[2]计算了薄壁墙受冲击载荷作用的响应问题。方法是：用FVM方法和流体动力学方程，对爆炸场进行数值模拟，用FEM方法对结构响应进行数值分析。霍宏发对单、双层椭球封头厚壁圆柱壳进行了较为系统的实验研究[3]，采用商业软件ANSYS-LSDYNA计算了爆炸流场，并解耦分析了封头和圆柱壳的振动模式。本文目的是：发展固体和流体的三维解算器和计算软件，耦合求解流体守恒方程和固体控制方程，解决爆炸载荷的精确确定和动边界处理等问题，致力于开展真正意义上的流、固耦合计算研究，认识冲击载荷作用下的流固耦合现象和机理。 控制方程及数值方法 本文研究对象是带椭球封头的圆柱壳。与文献[4]类似，其过程包括冲击波产生、传播和壳体动力响应三个相互关联的阶段。本文求解思路和文献[4]类似。由于本文重点关注冲击载荷作用下的流、固耦合问题，爆源采用简化模型：由点爆炸自相似解给出爆炸流场计算的初场，描述炸药爆炸冲击波产生。采用三维ALE型Euler方程描述爆炸中场，具体地 　　　　　　　 (1) 其中，，，，和为控制体积和表面积。。分别代表密度、压力和x、y、z方向速度分量和单位体积总能。状态方程取 　　　　　　　　 (2) 其中，为比热比。采用FVM方法和PPM格式[4,5]离散方程(1)，控制体单元为六面体。要说明的是：在轴线处作了挖空的近似处理。 壳体计算采用Lagrange坐标系下的虚功原理。具体为 (3) 其中，为质量密度，为加速度，为Cauchy应力，是当前位置的体力密度，是作用在的表面力，为虚位移。方程（3）求解采用位移模式下、线性基函数的FEM方法。 当冲击波在壳体内壁面反射，流固边界开始运动。在流固内界面，流体和固体网格粘结，流体网格速度等于壳体内壁面运动速度，壳体受力等于流体网格的压力。对中心装药，壳体封头顶点仅考虑z向自由度。沿各坐标轴，流场和壳体取对称边界条件。壳体外边界采用自由边界。耦合计算中，以壳体稳定性条件要求的时间步长为基准，先计算爆炸流场，再计算壳体。壳体隐式时间积分采用Newmark方法，使壳体计算时间和流场同步。 计算结果分析和讨论 为与实验结果对比，本文选择文献[3]的实验装置进行数值模拟。圆柱壳材料为普通碳钢。取物理域的1/8作为计算域（网格见图1），中心装药为1kg球形TNT。壳体壁厚53mm、内半径1060mm。封头长、短轴比为2：1、封头顶点到中心距离为2120mm。采用率形式的线弹性本构，杨氏模量，泊松比。流场和壳体计算网格分别取31×31×31、7×31×61。自相似解结束压比取为15。 图2给出了本文和文献[4]得到的壳体中环面（离中心最近）无量纲超压时间曲线。时间和压力参考量分别为160微秒、1.013×105pa。图2表明：计算得到的压力时间曲线与实测结果超压峰值差别很小。 图1　计算模型及网格分布（无量纲尺寸） 计算结果 b) 实验结果 图2 壳体中环面内壁超压时间曲线 图3给出了压力和密度云图表示的爆炸场演化过程。图3(a)给出了自相似解映射到三维流场的图像，由于三维和一维计算网格的差异，因此，自相似解映射后的冲击波阵面不太光滑。图3(b)表明：第一道球形冲击波