定向式动能杆战斗部飞散实验研究

动能杆战斗部是一种新概念防空反导战斗部。采用软管X光闪光照相技术测试了定向式动能杆战斗部在不同装药厚度、不同装药宽度以及不同杆间距情况下动能杆的飞散情况。由动能杆的散布图像计算了动能杆的飞散速度及飞散角。研究结果给出了装药宽度、装药厚度、杆间距等对动能杆飞散形状、飞散速度和飞散角的影响的规律性认识。     

带隔板装药爆轰波马赫反射理论研究和数值模拟

基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。     

装药比对中心炸药抛撒钨颗粒群影响的实验研究

在炸药抛撒颗粒的过程中,颗粒与炸药质量之比对抛撒颗粒的动能特性有重要影响。对中心为8710药柱,周围为钨颗粒的圆柱形装药进行了实验。实验中装药中心药柱的质量为200g,钨颗粒的质量与药柱的质量比(简称装药比)在2～7之间。实验结果表明,在有效作用距离内(2m～6m)颗粒数密度随装药比的增加呈正比增加;在近距离上(2m,3m),随装药比增大颗粒对纸靶的最大穿深明显下降,但在较远的距离上(4m,5m,6m),各装药颗粒对纸靶的最大穿深比较接近,随装药比增大而下降的规律不明显。颗粒对纸靶的平均穿深随装药比变化的规律与最大穿深相似。纸靶单位面积上的总穿深对装药比不敏感,随装药比的增大略成线性增大的趋势。     

浅层水中沉底的两个装药爆炸的数值模拟研究

根据试验模型和试验结果,进行了浅层水中沉底的两个装药爆炸的数值模拟;通过与水中单个装药爆炸,以及无限水中两个装药同时爆炸的数值模拟结果的对比分析,研究了水底水面对沉底的两个装药同时爆炸产生的冲击波传播与相互作用的影响。结果表明:水底对冲击波压力峰值有削弱作用,水面使冲击作用冲量明显减小,冲击波相互作用压力叠加或多次冲击作用可提高爆炸威力。     

十字形内置破片定向战斗部破片的飞散特性

针对低附带弹药毁伤需求，设计了一种十字形内置破片定向战斗部，根据目标方位可选择不同起爆模式，进而控制破片的径向飞散特性，在目标区域内形成杀伤破片实现定向毁伤，在非目标区域内实现低附带毁伤。采用数值模拟研究了相邻2点起爆、相邻3点起爆两种模式下定向战斗部起爆时破片的驱动过程，给出了各个位置处破片的飞散速度、径向飞散角度等特征参数；制备了2发单元样弹并开展了地面静爆实验，通过高速摄影及靶板上破片的穿孔分布特征实测出破片的速度及径向飞散角，与数值模拟结果对比，验证了模拟的准确性。在此基础上，通过引入能量分配角建立了破片速度的修正公式，并根据模拟结果对公式参数进行了拟合分析。结果表明：相邻2点起爆、相邻3点起爆模式下，战斗部定向杀伤区破片径向飞散角分别控制在145°、65°以内，且该区域内的破片占破片总数的比例分别达到了50.4%、43%；同时，破片速度呈现梯次分布，介于535～770 m/s之间，对1.5 mm厚的Q235A钢板的穿甲率分别达到了94.4%、84.6%，可实现对轻型车辆类目标的毁伤，其余区域则为低附带安全区；基于能量分配模型求得的破片速度理论计算值与模拟值基本吻合。研究结果可...     

波形控制器对杀伤战斗部破片飞散特性影响研究

为提升杀伤战斗部破片轴向飞散的集中度，提高战斗部的轴向杀伤威力，提出使用波形控制器控制破片的飞散方向。基于爆轰波在波形控制器界面发生反射的规律以及Shapiro公式，设计了波形控制器的形状，使用LS-DYNA有限元软件和ALE(arbitrary Lagrange-Euler)算法对破片的飞散过程进行数值计算，结合战斗部原理样机静爆试验，验证了使用波形控制器改善破片飞散特性方法的合理性。对比了有无波形控制器时破片飞散过程的差异，对无波形控制器以及波形控制器材料分别为尼龙、聚氨酯和聚四氟乙烯（polytetrafluoroethylene,PTFE）时杀伤战斗部的破片飞散速度和破片飞散角规律进行了分析。结果表明：波形控制器可以减小战斗部中心和两端位置的破片飞散速度大小差异，使中心到两端位置的破片飞散方向角变化均匀，破片在轴向的分布更加均匀；不同材料的波形控制器对破片飞散特性影响不同，波形控制器的使用减小了破片飞散角，增大了破片分布密度，提升了破片飞散的集中度。破片飞散角数值计算值与试验计算值误差在6.53%之内，与无波形控制器的杀伤战斗部原理样机相比，含波形控制器且材料为尼龙、聚氨酯和P...     

钨合金杆材料属性与侵彻性能的关系

通过对钨合金杆侵彻半无限厚铝合金靶问题的数值模拟,研究了高速侵彻中动能杆材料属性与侵彻性能的影响关系.计算采用的弹塑性流体本构模型中,用来描述材料属性的主要参数有:失效应变、屈服强度、剪切模量等.数值结果表明,失效应变是影响侵深的主要参数,而屈服强度和剪切模量对侵深的影响很小.另外,在此基础上,对着速在785 m/s～1924 m/s范围内的动能杆侵彻问题进行了数值模拟,计算结果与试验数据符合较好.     

球形和等长径圆柱装药在爆炸近区内载荷差异的实验研究

为了研究爆炸近区内装药形状对爆炸载荷的影响,运用压杆法测量了爆炸容器内壁的压力。通过 实验发现,对于圆柱装药径向方向上的正反射爆炸载荷,等长径圆柱装药的压力峰值、冲量与相同当量下球形 装药对应量的比最大超过2。还通过实验验证了数值模拟,并运用数值模拟获得2种装药形状下爆炸载荷的 静水压、密度、粒子速度,得到爆炸近区内2种装药形状载荷的差异以及该差异随爆心距离的变化规律。     

多点起爆下鼓形战斗部的威力特性

为调控对地弹药破片杀伤威力场,研究了一种鼓形战斗部在静爆和动爆下的威力特性。采用数值模拟研究了端面中心单点、中心单点两种起爆方式下,鼓形战斗部相较于同口径圆柱形战斗部在静爆下的破片威力特性及动爆下对地面装甲车辆目标的毁伤面积。在此基础上,调整鼓形战斗部起爆方式为偏心两线同时起爆、偏心两线序贯起爆及偏心两线同时-序贯起爆,计算了不同偏心起爆下鼓形战斗部静爆时的破片速度、飞散角和动爆时对车辆目标的毁伤面积及有效破片落地动能分布。研究表明,相比于同口径的圆柱形战斗部结构,鼓形战斗部的破片飞散角增大了55.98%,对地面军用车辆的毁伤面积最大增大了59.3%;相对于偏心两线同时起爆,偏心两线同时-序贯起爆的鼓形战斗部破片飞散角增大了18.0%,破片飞散的离散程度提高了11.48%;相对于装药中心单点起爆,偏心两线序贯起爆下鼓形战斗部的毁伤面积受炸高影响较小,在落角50°、落速200 m/s、炸高为9 m时的毁伤面积达47.15 m²。通过调整战斗部的结构和起爆方式,可有效增大破片的飞散角,增大破片对目标的覆盖面积,提高战斗部的毁伤效能。     

偏心对称起爆战斗部破片初速的增益

在极坐标下基于弹塑性基本方程建立了壳体膨胀的运动方程,且考虑了两点对称偏心起爆的碰撞效应,利用Whitham方法对两对称爆轰波的碰撞叠加进行了计算,得到了三波点的迹线、马赫波超压及马赫杆的高度等参数。通过联合求解,导出了马赫波区破片的初速计算公式,并利用AUTODYN软件进行了数值模拟,理论计算与数值模拟结果符合较好,验证了理论模型的可靠性。计算结果表明,两点对称偏心起爆时定向区破片初速增益超过30%,起爆点的夹角变化对破片的初速大小及飞散偏转影响并不明显。     

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

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

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


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.     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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).