内部爆炸载荷作用下混凝土动力响应的数值模拟

 目前,采用实验方法对爆炸载荷破坏效应进行全面、准确描述仍是一个难题。运用ANSYS/LS-DYNA非线性显式动力学有限元程序,采用流固耦合算法,对混凝土在内部爆炸载荷作用下的空腔形成和发展规律等动力响应问题进行了数值模拟。通过实际算例,比较研究了混凝土敷设钢板与否对成坑大小、损伤区、压力云图以及应力时程曲线等的影响,进一步研究了有钢板加固时混凝土在爆炸荷载作用下的破坏机理。研究表明：采用流固耦合算法可以避免由于物质的大变形引起的单元畸变和网格消蚀等问题,提高了计算效率;很好地模拟了混凝土内部爆炸时,中心爆炸空腔压缩面和自由面的拉伸层裂破坏过程,且与理论计算值吻合较好;通过表面敷设钢板,可以减小混凝土的拉应力峰值,抑制混凝土中裂纹的发展,从而提高其防震塌能力以及结构的完整性。     

椭圆截面战斗部爆炸驱动破片作用过程的数值模拟

为研究爆轰驱动下椭圆截面自然破片杀伤战斗部壳体的膨胀破裂过程以及壳体破片径向速度分布,建立了椭圆截面战斗部三维模型。通过AUTODYN-3D软件,采用Lagrange算法模拟爆轰驱动下椭圆截面自然破片战斗部壳体的膨胀断裂过程,研究了端面单点中心起爆方式下短长轴断裂时间差与短长轴比的关系,以及不同起爆点、不同短长轴比和不同装填比(即装药与壳体质量之比)对椭圆截面战斗部径向破片速度分布的影响。结果表明:与端面中心单点起爆、端面长轴双点偏心起爆和端面短长轴四点偏心起爆相比,端面短轴双点偏心起爆方式对椭圆截面战斗部壳体破片径向速度的增益效果最好。装填比一定时,短、长轴断裂时间以及短、长轴断裂时间差与短长轴比呈线性关系,战斗部壳体膨胀过程中截面形状的实时短长轴比与加载时间呈线性关系;随着短长轴比的增大,战斗部壳体破片径向速度增益逐渐减小。短长轴比一定,装填比小于1时,破片速度随方位角增大呈正弦趋势上升,且短、长轴方向破片速度差与装填比呈线性关系。     

柱形容器内部爆炸流场的数值模拟

采用NND差分格式,同时利用点爆炸解作为计算初值,对圆柱形容器内部的二维轴对称非定常流场进行了数值模拟,清楚地反映了内部流场的变化过程,得到了容器壁面上各点的压力-时间关系。对计算结果进行了讨论,与文献中结果作了比较。     

水射流对煤体冲击的有限元与光滑粒子耦合法数值模拟

采用光滑粒子与有限元耦合算法对高压水射流冲击煤体进行了数值模拟,建立了长为5mm、半径为2mm的圆柱形水体以不同速度冲击煤体的计算模型。根据计算结果,分析了煤体在高压水射流冲击下的损伤模式、煤体中的应力波传播形式、煤体在高压水射流作用下的临界破煤压力。     

内部爆炸载荷作用下砌体墙结构的失效规律

为了探究钢筋混凝土框架-砌体墙建筑物的内爆毁伤效应,开展了内爆载荷下砌体墙结构动力响应和失效规律研究。采用数值模拟方法,并结合理论分析,研究了不同当量内爆载荷下砌体墙的失效机理以及失效形式转化过程。结果表明,小药量内部爆炸时,砌体墙的主要失效形式为弯曲导致的墙面开裂。随着药量增加,墙体边界在首道冲击波反射超压作用下发生以剪切变形为主导的失效。而在大药量条件下,首道冲击波超压使砌体墙材料达到极限抗压强度而发生压溃失效。研究结果可为砌体墙结构毁伤评估与防护设计提供技术参考。     

深水爆炸冲击波作用下圆柱壳动态响应影响因素的数值模拟研究

 运用通用有限元程序ABAQUS,对圆柱壳在深水爆炸冲击波作用下的动态响应进行了数值模拟,研究了圆柱壳所处深度、爆心方位和预应力对于圆柱壳动态响应的影响。研究结果表明：在同样的水下爆炸冲击波作用下,随着深度的增大,圆柱壳的毁伤逐渐加重,而且响应过程中的变形逐渐增大,低频响应阶段的速度曲线趋于平坦;爆距不变时,不同的爆心方位在不同深度上的毁伤效果存在一定差异;预应力的存在使圆柱壳的毁伤加重,并且预应力对圆柱壳动态响应的影响与爆心方位有关。研究结果对于潜体结构的生命力评估具有一定参考价值。     

水下爆炸特性的二维数值模拟研究

 对带有3 mm左右厚铝壳或钢壳的柱形TNT、RS-211、T/Γ（35/65）以及ROT901装药在水下采用一端点起爆后,水中冲击波超压分布进行了二维数值模拟研究。对该四种炸药的每一种,分别考虑了药量为1 kg、3 kg、5 kg和8 kg四种TNT当量的柱形装药情况。柱形装药的长径比为L/D=1.5左右。计算使用的程序为二维欧拉多流体网格法流体力学计算程序MFIC。给出了水中冲击波超压的空间和时间分布;对铝壳和钢壳两种情况的结果进行了分析、比较;同时给出了缩比模型与1∶1模型之间相应物理量之间的缩比关系。     

物理质团法:一个普适的数值方法体系

无网格方法根据分布于近邻空间各个方向的微元体物理信息构造离散方程,显著降低了空间导数计算对微元体本身及微元体之间拓扑结构的条件限制,极大提高了拉氏方法的大变形计算能力.由于不能利用微元体的完备几何信息,不容易构造符合物理的无网格算法,对那些物理参数存在间断的模型对象,难以获得稳定和准确的计算结果.本文基于对物理规律及数值模拟发展趋势的分析,提出符合物理且具有强普适性的无网格方法体系.基于该方法的一维算例表明,即使物理参数存在强烈间断,数值结果也能很好地逼近问题的真解.     

水中障碍物爆炸毁伤效应的二维数值模拟

 建立了水下障碍物的爆炸毁伤效应的物理模型和计算模型。建立了守恒方程组、混凝土的动态弹塑性本构关系和破坏准则,并确定了其在大应变、高应变率及高压条件下的的屈服准则。基于多物质流体的Euler算法,用面向对象的C++语言自行编制了NM-MMIC通用多物质二维流体弹塑性程序,对不同炸点位置、不同炸药量对轨条砦爆炸破坏过程进行数值模拟,并采用自行开发的后处理软件VISC2D实现了可视化,定量地反应了障碍物迎爆面所承受的压力变化过程,确定了轨条砦毁伤的最佳药量,并对毁伤程度进行了评价,得到了爆炸破障的毁伤规律。研究结果对设计和评价反水中障碍物的弹药具有一定的参考价值。     

串联弹前级装药对后级随进弹爆炸作用的研究

 为获得串联战斗部随进弹在受到前级装药爆炸作用时,质量对其所受最大过载及动能衰减的影响规律,采用非线性动力学软件AUTODYN,基于Euler-Lagrange耦合算法,对前级装药的爆炸过程中,不同弹径、不同质量的随进弹穿越爆轰场进行了数值模拟。得出了随进弹在前级爆炸作用下所受最大过载和动能衰减的变化规律。基于数值模拟结果建立了随进弹剩余速度与质量的无量纲关系式,计算结果得到了实验验证。     

点起爆炸药驱动平板飞片运动的数值模拟研究

 为研究点起爆发散爆轰波驱动飞片的运动特性,自行研制了二维流体动力学程序TDY2D,采用由实验所标定炸药的状态方程参数,对3种长径比相差较大的点起爆炸药驱动平板飞片运动的模型进行了数值模拟计算,并对结果进行了分析。计算所得的飞片运动历史和击靶波形均与实验符合较好,具有较高的计算精度,从而验证了该程序的准确性和有效性。     

缺陷对圆管膨胀断裂影响的实验研究

采用高速分幅摄影技术,通过对比实验,获得了清晰的无缺陷和预制缺陷金属圆管膨胀断裂的演化图像.实验结果表明,缺陷大大加速了圆管的断裂过程;有缺陷圆管的产物开始泄漏的时间比无缺陷圆管提前5.5μs,采用应力/应变集中法很好地解释了这种现象.另外,针对预制缺陷圆管膨胀断裂时出现的裂块较小而数量较多、无缺陷圆管膨胀断裂时裂块相...     

HMX炸药燃烧转爆轰数值模拟

以两相流模型为基础,气相产物状态方程采用基于统计物理的类CHEQ的计算结果,建立了HMX炸药的燃烧转爆轰数学模型。采用CE/SE方法模拟了颗粒度为125μm的HMX炸药的燃烧转爆轰过程,得到了爆轰参数及流场变化规律。模拟了装填密度对HMX炸药燃烧转爆轰的影响,并与实验进行了对比。数值模拟结果表明,在相同的点火条件下,爆轰成长距离在一定范围内随装填密度呈"U"形变化。     

Numerical investigations of the effect of oblique impact on particle deformation in cold spraying by the SPH method

In this study, a systematic examination of the oblique impacting of copper particles in cold spraying was conducted by using the smoothed particle hydrodynamics (SPH) method compared to the Lagrangian method. 3D models were employed owing to the asymmetric characteristic of the oblique impacting. It is found that in the oblique impact, the additional tangential component of particle velocity along the substrate surface could create a tensile force and decrease the total contact area and bonding strength between the particle and the substrate. The simulation results compare fairly well to the experiment results. Meanwhile, the asymmetric deformation can result in the focus of the shear friction on a small contact zone at one side, which may rise the interfacial temperature and thus facilitate the occurrence of the possible shear instability. Therefore, there probably exists an angle range, where the deposition efficiency may be promoted rather than the normal angle. Moreover, the particle deformation behavior simulated by the SPH method is well comparable to that simulated by the Lagrangian method and the experimental results, which indicates the applicability of the SPH method for simulating the impact process in cold spraying besides the previously used Arbitrary Lagrangian Eulerian (ALE) method.     

隔板深度对激波聚焦起爆影响的数值模拟研究

考虑几何结构参数对激波聚焦触发爆轰波的复杂影响,对H2/Air预混气的环形射流激波聚焦起爆现象开展了数值模拟研究,详细分析了不同隔板深度条件下的激波聚焦过程、流场演化特征以及爆轰波参数变化规律.研究结果表明,凹腔内激波聚焦诱导的局部爆炸以及隔板前缘处射流形成"卷吸涡"是引起爆轰波触发的两个重要机制,而隔板深度是影响环形...     

钢筋混凝土靶板在弹丸冲击及爆炸载荷下响应的数值模拟

 采用HJC混凝土损伤本构模型及LS-DYNA的流固耦合算法,分别对钢筋混凝土靶板在弹丸冲击和爆炸载荷作用下的响应进行了有限元数值模拟,其中模拟参数由实验数据拟合重新获取。将模拟结果与实验结果和经验公式进行对比分析,结果表明：数值模拟再现了弹体贯穿靶板过程中的开坑、隧道及漏斗碎裂区,计算得到的弹体弹道极限及残余速度与实验数据吻合较好;此外,数值模拟也很好地再现了炸药爆炸后冲击波的传播过程以及爆炸载荷作用下混凝土的破坏情况,模拟结果与实验现象具有良好的一致性。