Shape optimization of impactor penetrating into concrete or limestone targets

The shape of the normally striking impactor that attains the maximum depth of penetration into a concrete or a limestone semi-infinite target for a given impact velocity is found. It is shown that the optimum shape is close to a blunt (in general case) cone and it is independent on the properties of the material of the target in the framework of the employed penetration model. The performance of some other typical shapes of the nose of the impactors (spherical-conic impactors, sharp-conic impactors, truncated-ogive impactors) are analyzed and compared with the optimal impactor.     

Optimal 3D impactors penetrating into layered targets

Optimization of 3D sharp high speed impactors with given form of a longitudinal contour, length, and volume, penetrating into layered ductile targets, both for conical and thin non-conical strikers using approximate models is studied. It is found that the impactor with the minimum drag moving in a homogenous target with a constant velocity penetrates to the maximal depth into a semi-infinite target and has the minimal ballistic limit when it penetrates into a finite thickness target, regardless of the distribution of the material properties of the target along its depth, the number of the layers, etc. Using the analogy with the hypersonic flow over the flying projectiles it is predicted that the optimal impactor should have a star-shaped form of the cross section. If an impactor has a polygonal cross sections allowing the inscribed circles, the ballistic limit and maximum depth of penetration are independent not only of the properties of the target but also of the form of the polygon in the cross section and equal to the corresponding values for the inscribed body of revolution.     

Penetration of rigid bodies into loose and layered media

Penetration and motion of rigid bodies in ground media attracts the researchers’ attention because of various problems arising as the technology evolves. In fact, there are two independent directions of studies in this field: (1) the problem of earth excavation when a rigid body of a definite shape slowly moves along a given trajectory in the ground; (2) an impact of a rapidly flying free rigid or deformable body against the ground. In the latter case, to which the proposed studies pertain, it is sometimes of interest to study the medium behavior and the motion of the free body, which moves in the medium after the impact owing to the kinetic energy of itself. In this field, a majority of studies deal with collision and penetration of bodies of various shapes into clay media. An extensive survey of these studies is given in [1]. After this survey appeared, numerous paper dealing with complicated collision conditions have been published [2]. Penetration in loose media has been studied much more rarely. The direct collision with fractured rock was studied in connection with the expected landing of spacecraft on other planets [3, 4]. In this case, the influence of grain dimensions and the density of the filling and vacuum on the penetration was studied for the initial velocities in the range of 1.7–10 m/s. On the other hand, in [5], the results of investigating the penetration of conic bodies in sand at entry velocities of 700–900 m/s are given; these velocities significantly exceed the speed of sound in this medium, which lies in the range of 100–200 m/s for dry sand. Analyzing the experimental results, the author came to the conclusion that it is necessary to use different representations of the drag force in the supersonic and subsonic modes. In the present paper, we do not consider the influence of the grain distribution, sand density, and filling methods on the penetration. But, as follows from the experiments whose results are described in [6] and [7], to represent the results of penetration of rigid bodies at velocities up to several hundreds of meters per second, in addition to the characteristics listed above, it is also required to describe the technology of the experiment preparation, because such media have the property of shape “memory.”     

Localized interaction models with non-constant friction for rigid penetrating impactors

We suggest approximate penetration models for rigid body penetration that take into account sliding velocity (SV) and pressure dependence of the friction coefficient (FC). It is showed that introducing variable FC in a localized interaction model (LIM) yields a model that belongs to the class of LIM. We developed a general method for determining the depth of penetration (DOP) using the piecewise linear approximation of the impactor’s generatrix. For some classes of SV dependent friction models we obtained analytical formulas for calculating the DOP. Using the experimental data available in the literature, we determined the dependencies of FC vs. pressure and SV. We conducted numerical modeling of penetration of a metal striker into metal and concrete shields employing models with variable and constant FC. Numerical simulations showed that taking into account variable FC strongly effects the DOP when FC changes appreciably for large velocities that are characteristic for the high-speed penetration.     

On the penetration of a rotating impactor into an elastic-plastic medium

This research is devoted to the modeling of high-speed rectilinear penetration of a rigid axisymmetric body (impactor with a flat bluntness) into an elastic–plastic media with account for its rotation about the axis of symmetry. The body has an arbitrary shape of the meridian. The resistance to the motion is represented as the sum of the body drag and the contribution of friction. The dynamic system governing the body motion is derived and the qualitative and numerical analysis of the projectile movement and perforation of a slab are performed. The problem of shape optimization of impactor with a flat bluntness is studied using evolutionary algorithm.     

卵形头部刚性弹侵彻厚靶的半解析模型

基于确定靶体中速度势和速度场的方法分析刚性卵形头部弹体对有限厚靶的侵彻问题。推导了靶体中速度场与应力场的计算方法,利用据此编制的计算程序,计算了卵形头部钢弹体对铝靶的侵彻与穿透问题,给出了侵彻深度与剩余速度同初始碰撞速度的关系。结果表明,在对实验参数不经过任何调整的情况下,得到了同试验曲线相吻合的结果。可以看出该方法的有效性。     

弹体高速侵彻钢筋混凝土靶试验研究

为研究结构弹体对钢筋混凝土靶的高速侵彻破坏效应,利用口径35 mm弹道炮开展了1 030～1 520 m/s速度范围内的高速侵彻试验,获得了弹体的撞击速度、破坏形态、剩余长度、剩余质量和靶体中的侵彻深度及成坑尺寸等试验数据,分析了侵彻深度和侵彻机理随速度的变化关系。结果表明：在1 030～1 390 m/s的速度范围内,弹体头部磨蚀,磨蚀程度随侵彻速度增加而加剧,侵彻深度随撞击速度近似线性增大;撞击速度在1 390～1 480 m/s范围内,弹体头部严重磨蚀,侵彻深度随撞击速度增加而减小;撞击速度大于1 480 m/s后,弹体严重破碎,侵彻深度急剧下降。针对结构弹体高速侵彻过程中的破坏特点,将侵彻速度划分为刚体侵彻区、准刚体侵彻区、侵蚀体侵彻区和破碎体侵彻区,可为钻地弹结构设计提供参考。

     

Slender bodies of revolution with minimum wave drag in nonequilibrium supersonic flow

We examine the problem of finding the generatrix shape of a body of revolution which travels at supersonic speed and has minimum wave drag. We assume that any number of nonequilibrium processes can take place in the flow. The pressure distribution over the body surface is taken in the linear approximation [1, 2]. A survey of studies using linear theory to find bodies of revolution of optimal form in supersonic perfect gas flow can be found in [3]. The solution of the problem of finding the form of two-dimensional slender bodies of minimum wave drag in nonequilibrium supersonic flow was obtained in [4]. In the following we examine the optimization of only those bodies of revolution for which the leading point lies on the axis of symmetry.The author wishes to thank A. N. Kraiko for his helpful comments.     

Application of the angular superposition method to the contact problem on the compression of an elastic cylinder

The angular superposition method is used to construct an approximate solution of the contact problem on the compression of an elastic cylinder by two rigid plates. The solution thus obtained has a closed-form analytic expression and can be used in the entire domain of the cylinder cross-section. We analyze the absolute error, which takes the largest value near the points of contact between the plates and the cylinder, where the boundary conditions are discontinuous. According to the von Mises criterion, when moving into the depth of the cylinder from the contact site along the symmetry axis, the second invariant J ₂ of the stress deviator tensor first decreases and then, after attaining a minimum, increases and attains the largest value at a small depth, which agrees with Johnson’s photoelastic experiments and Dinnik’s computations. We present the graphs of the displacement and normal stress distributions over the contact site, the dependence of the compressing force on the displacements of rigid plates, and the dependence of the invariant J ₂ on the coordinate along the symmetry axis. If 640 computation points are chosen on the cylinder boundary and the Hertz law for the normal pressure on the contact site is used, then the error in the approximate solution near the endpoint of the contact site is approximately 55%, and if the proposed two-parameter normal law is used, then the error is of the order of 4%. On the free lateral surface of the cylinder boundary, we find the critical pointM*, which separates the cylinder contraction and extension parts.The contact problems are the most difficult problems, and their solution is complicated by the discontinuous boundary conditions [1–5]. In [6], the contact problem is solved by the Fourier method, which can be used only for bodies of classical shapes. In such cases, the problem can be reduced to solving coupled integral equations [7]. The interaction between the bandage and a cylindrical body is considered in [2, 6, 7]. In [8], the possibility of using the finite element method is investigated in the case of contact problems for a differential wheel with roughness of the contacting surfaces taken into account. In [9, 10], the method of homogeneous solutions is used to consider contact problems for a finite-dimensional elastic cylinder loaded on its end surfaces. Note that only error estimates are given in the literature cited above; the absolute error over the entire domain of the elastic body is not studied, although this is one of the important characteristics of the obtained approximate solution. A sufficiently complete survey of the literature in the field of contact interactions of elastic bodies is given in [3–5].In what follows, we propose to solve contact problems by the angular superposition method [11]. This method can be used for bodies of nonclassical shapes, which can be multiply connected, and the friction on the contact site can be taken into account. In the present paper, as a first example of applied character, we show how this method can be used in the simplest case. The multiple connectedness and the curvilinearity of the shape of the body, as well as taking into account the friction on the boundary, do not create new essential difficulties in this method.     

动能深侵彻弹的力学设计(Ⅱ):弹靶的相关力学分析与实例

介绍了作者关于动能深侵彻弹研究的部分工作。单纯动能弹侵彻存在侵彻深度上限，在技战术指标一定的前提条件下，弹体结构的优化可令动能深侵彻弹尽可能实现最大侵深。其中弹体结构的力学设计尤为重要，如弹形、弹材、质量比、长径比等。进一步开展弹体结构设计的力学研究，包括头部设计、战斗部后盖、装填比、弹形和弹材、侵彻效应中的尺度律、混凝土靶的设计，最后给出若干实例。     

On the stress state of shells penetrating into a deformable solid

The motion of an axisymmetric shell in a deformable solid medium is considered. It is assumed that the medium resistance is described by a two-term expression containing a constant term (the rigidity characteristic) and an inertial term quadratic with respect to the penetration velocity. A model of the impactor penetration with the normal interactions with the resisting medium taken into account is proposed. The membrane forces and the arising stresses are determined for decelerated motions of the impactor.     

Penetration of a copper rod into a sandy target

This paper presents the results of experimental and theoretical studies of high-velocity penetration of cylindrical copper rods into sand. The hydrodynamic Alekseevskii-Tate theory is modified to determine the penetration depth and wear velocity of the material of the rod penetrating into soil target in the plastic and hydrodynamic stages of penetration. The case where the target material is significantly less strong than the rod (impactor) material is considered.     

An engineering analysis of penetration of metal ball into fibre-reinforced composite targets

An engineering analysis of computing the penetration problem of a steel ball penetrating into fibre-reinforced composite targets is presented. Assume the metal ball is a rigid body, and the composite target is a transversely isotropic elasto-plastic material. In the analysis, a spherical cavity dilatation model is incorporated in the cylindrical cavity penetration method. Simulation results based on the modified model are in good agreement with the results for 3-D Kevlar woven （3DKW） composite anti-penetration experiments. Effects of the target material parameters and impact parameters on the penetration problem are also studied.     

Possibility of generating strong magnetic fields in conducting materials by the action of high-velocity penetrators

The possibility of increasing the intensity of a magnetic field produced previously in a conducting medium moving under the action of a high-velocity penetrating body is analyzed. A simplified model of the interaction of an impactor and a conducting target with a transverse magnetic field is constructed within the framework of a one-dimensional scheme. It is shown that the degree of increase in the field intensity is determined by the relation between the magnetic-field compressibility and diffusion factors, and the corresponding dimensionless determining parameters is determined. Magnetic-field compression is estimated for a perfectly conducting medium and media with real conductivity. The significance of the thermal and mechanical effects accompanying the penetration of an impactor into a target with a transverse magnetic field is assessed. Bauman Moscow State Technical University, Moscow 107005. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 3, pp. 13–18, May–June, 2000.     

Optimal shapes of axisymmetric bodies penetrating into soil

This paper presents the results of a study of the shapes of axisymmetric bodies with minimum drag and maximum depth of penetration into the plastic soils. Optimal shapes of bodies of revolution of given length and cross-sectional radius with generatrices represented by line segments are obtained by a modified method of local variations. The problem is solved using a binomial quadratic model of local interaction, including inertial and strength terms containing constant and Coulomb frictions. The drag forces and the penetration depth of cones and the obtained bodies of optimal shape are determined at different penetration velocities.     

Motion of a cylinder in a viscous electroconductive medium with a magnetic field

The method of force sources is used to consider the planar problem of the motion of a circular cylinder in a viscous electroconductive medium with a magnetic field. The conventional and magnetic Reynolds numbers are assumed to be small. Expressions are obtained for the hydrodynamic reaction forces of the medium, acting on the moving cylinder. It is shown that as a result of the flow anisotropy in the medium, caused by the magnetic field, in addition to the resistance forces on bodies moving at an angle to the field, there are deflecting forces perpendicular to the velocity vector. The velocity field disturbances at great distances from the moving cylinder are determined.The problems of viscous electroconductive flow about solid bodies in the presence of a magnetic field constitute one of the divisions of magnetohydrodynamics. Motion of an electroconductive medium in a magnetic field gives rise to inductive electromagnetic fields and currents which interact with the velocity and pressure hydrodynamic fields in the medium [1, 2]. Under conditions of sufficiently strong interaction, the number of independent flow similarity parameters in MHD is considerably greater than in conventional hydrodynamics. This circumstance complicates the theoretical analysis of MHD flow about bodies, and therefore we must limit ourselves to consideration of individual particular flow cases.Here we consider the linear problem of the motion of an infinite circular cylinder in a viscous incompressible medium with finite electroconductivity located in a uniform magnetic field.There are many studies devoted to the flow of a viscous electroconductive medium with a magnetic field about solid bodies (see, for example, [3–5]). Because of this, some of the results obtained here include previously known results, which will be indicated below. In contrast to the cited studies, the examination is made by the method of force sources, suggested in [6]. This method permits obtaining integral equations for the distribution of the forces acting on the surface of the moving body. Their solution is obtained for small Reynolds and Hartmann numbers. Then the nature of the velocity disturbances at great distances from the body are determined. These results are compared with conventional viscous flow about a cylinder in the Oseen approximation.     

Complete integrability of the equations of motion of a spatial pendulum in a medium flow with rotational derivatives of the torque produced by the medium taken into account

We construct a nonlinearmodel of the mediumaction on a rigid body taking into account the dependence of the force arm on the reduced angular velocity of the body. In this case, the moment of the action force itself is also a function of the angle of attack. As experimental data processing for the motion of homogeneous circular cylinders in water has shown, it is necessary to take these facts into account in modeling.Studying the model of interaction between the spatial pendulum and the medium, we found a new case of complete integrability in elementary functions. This allowed us to find several qualitative analogies between the motions of bodies that are free in the resisting environment and the oscillations of bodies partially fixed in the homogeneous flow of incoming medium.     

刚体单元及其在多体系统动力学中的应用

多体系统动力学分析软件要求人工输入形状复杂物体的质量、质心位置和转动惯量,而实际上这些参量并不容易获得。本文探索了一种以组成物体的刚体单元为基本要素的新方法,并结合实际需要具体构造了刚性四面体和刚性梁单元。以刚体单元为基础并内嵌网格剖分模块的分析软件能够自动获得这些参数,从而具备处理任何复杂系统的能力。仿真结果的对比分...     

不同重力场下颗粒冲击过程的离散元分析

开展不同重力场下颗粒材料冲击动力学研究有助于加深对颗粒运动机制的理解和深空探测任务的实施。本文采用离散元模拟对颗粒介质受球形冲击物的冲击过程进行了数值模拟,并与地球重力下冲击的试验结果进行对比验证。在此基础上,进一步研究了重力加速度对冲击物动力学的影响规律。计算结果表明,在所有重力加速度下,冲击物的穿透深度d与冲击速度v₀的关系可以用Poncelet模型表达;d与冲击物下落的总高度H表现为d～Hⁿ的幂律关系,当H<10 m时,d与H的幂率标度为0.322,而H>10 m时,d与H的幂率标度下降到0.211。此外,穿透深度小于冲击物半径时,重力加速度对冲击物减速过程无影响。在所有的重力加速度下,当冲击速度大于5 m/s时,冲击物的持续碰撞时间t_c是恒定的,且与重力的-1/2次方呈正比。