三点起爆控制参数对尾翼爆炸成型弹丸成型的影响

针对三点起爆控制参数（起爆直径、起爆同步误差）对尾翼爆炸成型弹丸（explosively formed projectile,EFP）成型的影响问题,理论分析了三点起爆条件下爆轰波的马赫碰撞,计算获得了不同起爆直径下马赫波的相关参数,利用LS-DYNA有限元软件通过数值模拟研究了不同起爆直径下三点起爆同步误差对EFP尾翼成型及飞行速度的影响规律。结果表明,起爆直径越大,马赫波在药型罩上的作用面积越小,马赫超压越大,形成的EFP长径比和速度越大;起爆直径为30、40、50 mm三点起爆成型装药形成较佳尾翼EFP应满足的最大起爆同步误差分别为50、100、150 ns;此外,尽量使中间起爆点起爆同步误差约为最大同步误差的一半,有利于降低尾翼EFP的侧向分速度,提高飞行稳定性。     

靶板在爆炸成型弹丸垂直侵彻下的层裂

为了准确掌握靶板层裂过程和规律，基于波动力学和基本假设，建立了爆炸成型弹丸（explosively formed projectile，EFP）垂直侵彻有限厚靶板时层裂的力学模型，得到了层裂点的表达式。研究结果表明:EFP速度为1 800 m/s、靶板厚度从35 mm增大到60 mm时，靶板背面弯月形层裂区厚度不断增大，弯月形层裂区长度不断减小；靶板厚度保持40 mm不变、EFP速度从1 600 m/s增大到1 900 m/s时，靶板背面层裂区厚度不断减小，弯月形层裂区长度不断增大。开展了EFP侵彻40 mm厚装甲钢靶板的实验，将实验结果和理论计算结果进行对比分析，两者吻合较好。     

三点起爆形成尾翼EFP的数值模拟和实验研究

利用LS-DYNA软件对三点起爆形成尾翼EFP的过程进行了数值模拟,深入研究了爆轰波传播过程中的波形结构和强度的变化规律以及药型罩材料在复合爆轰波作用下驱动变形的特性和规律,加深了对三点起爆条件下药型罩形成带尾翼EFP机理的认识。在此基础上,设计了三点同步起爆装置和EFP装药进行实验。研究结果表明:设计的三点起爆装置作用可靠,满足三点起爆EFP装药的设计要求;形成的EFP弹形稳定,与计算结果吻合较好;尾翼EFP飞行过程中的速度降减小,稳定性提高。     

含偏心起爆对EFP战斗部飞行特性的影响

为研究起爆不对称性对EFP战斗部飞行特性的影响, 对不同偏心量下?60 mm弧锥结合罩EFP战斗部进行飞行弹道实验。实验结果表明:偏心起爆条件下, 当相对偏心量小于3.3%时, EFP在网靶穿孔接近圆形, 弹丸飞行稳定; 起爆相对偏心量达到6.7%时, 弹丸飞行过程中摆动幅值增大, 降低了对目标的打击精度和毁伤效果。利用LS-DYNA及CFX非线性动力学有限元程序对不同起爆偏心量下成型EFP的空气动力学特性进行数值模拟, 描述了偏心起爆影响EFP成型对称性, 改变弹丸在飞行过程中流场的分布特征, 从而导致弹丸飞行过程中无规则运动的全过程。     

变截面爆炸成型弹丸垂直侵彻装甲钢板靶后破片质量模型

考虑爆炸成型弹丸(explosively-formed projectile，EFP)变截面的特性，基于流体力学Bernoulli方程和绝热剪切理论，改进了EFP垂直侵彻装甲钢板靶后破片质量模型，结合已有的试验数据和数值仿真方法检验了改进后模型的准确性。在此基础上，分析了靶板厚度和EFP着靶速度对靶板和EFP产生的靶后破片质量的影响规律。结果表明:相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随靶板厚度和EFP着靶速度的变化规律；当EFP着靶速度为1 650 m/s时，随着靶板厚度从30 mm增大到70 mm，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断增强；当靶板厚度为40 mm时，随着EFP着靶速度从1 650 m/s升高到1 860 m/s，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断减弱。     

切割式多爆炸成形弹丸成形及对钢靶的穿甲效应

基于60mm 弧锥结合罩EFP装药,设计了一种在药型罩前适当位置安装可抛掷的十字形网栅的 切割式多爆炸成形弹丸战斗部结构,并进行了靶场静爆实验。由实验结果可知,该战斗部经网栅切割后能形 成5枚具有一定质量和方向性、可贯穿6mm 厚45钢靶的弹丸,有效提高了毁伤元的数量和毁伤面积。利用 LS-DYNA程序对弹丸侵彻45钢靶过程进行了数值模拟,分析了弹丸侵彻钢靶过程。通过对中心弹丸穿靶 模拟数据的处理,得到了其余尺寸弹丸侵彻不同厚度45钢靶的极限穿透速度计算公式,该公式可对相关切割 式多爆炸成形弹丸战斗部威力优化设计和评估提供依据。     

网络爆轰驱动飞片的设计技术

提出使用炸药网络爆轰技术,设计了一种新型的低速飞片击靶加载装置,结合数值模拟和实验测 量,对驱动飞片的速度、平面性及装置结构参数进行了系统分析。实验结果表明,无氧铜飞片可以使直径 100mm、厚度3mm 的飞片在飞行距离2mm 内达到稳定的飞行状态,击靶速度在200~350m/s范围内连 续可调,飞片击靶平面范围大、速度均匀、平面性好,为实现低冲击应力的一维应变加载提供了参考。     

平面分布式火力毁伤元的实验研究

单个EFP在毁伤面积上存在一定的局限性.为扩大EFP毁伤威力和毁伤面积,本文研制了一种高初速的EFP战斗部,设计出一种含有46个EFP的平面分布式火力毁伤元.首先采用X光实验技术对弹丸进行X光拍摄,通过对X光照片的判读获取弹丸的外形尺寸和飞行速度,完成单个EFP战斗部的优化设计;然后在靶场对优化后的EFP战斗部进行侵彻实验,通过靶场回收得到弹丸的外形尺寸,采用测速靶得到弹丸的飞行速度,进一步验证所设计变壁厚药型罩形成高速弹丸的稳定性;最后将EFP战斗部按一倍战斗部直径的间距集成在一个由聚氨酯材料形成的平板上,采用高速摄影技术对起爆后弹丸的飞行轨迹进行拍摄,通过对高速摄影录像的判读,获取弹丸的飞行轨迹,由此判断弹丸在飞行过程中是否发生发散.实验结果表明,所设计的EFP弹丸飞行速度约为2200m/s并具有较好的密实性,火力毁伤元起爆后所形成的弹幕没有发生发散,且具有较好的飞行稳定性.本文设计的平面分布式火力毁伤元可根据对毁伤面积的要求调整布置,扩大了EFP的应用范围.     

一种形成尾翼型聚能侵彻体新方法的实验研究

尾翼可以提高聚能侵彻体的飞行稳定性,如何形成有一定规律性的尾翼在聚能侵彻体的研究中具有重要意义。在药型罩上粘附隔板,利用隔板改变爆轰波波阵面的结构形状,使药型罩上的爆轰压力发生规律变化,从而发生规律性的变形,最终形成带有尾翼的聚能侵彻体,这是一种新的形成尾翼的方法。从药型罩微元压垮速度变化和药型罩表面爆轰压力变化的角度出发,对新方法形成尾翼的机理进行了初步探讨。通过合理设计隔板的几何尺寸对该方法形成尾翼型聚能侵彻体进行试验,试验回收到了带有尾翼的聚能侵彻体,说明新方法形成尾翼具有一定的可行性。     

带尾翼的翻转型爆炸成形弹丸的三维数值模拟

利用LS-DYNA三维动力有限元程序对三点起爆的翻转型EFP形成过程及三个尾翼进行了数值模拟。计算结果与实验结果的比较表明:数值计算结果和实验结果吻合较好,可为弹丸优化设计提供重要参考依据。对药型罩材料分别选择Johnson-Cook和Steinberg两种本构模型进行了数值模拟。计算结果表明:药型罩材料的本构模型选取对形成的尾翼效果有一定影响,Steinberg本构模型与实验结果符合更好。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).