蜂窝遮弹层抗弹丸侵彻实验研究

为研究蜂窝遮弹层的抗弹丸侵彻性能,采用15 mm弹体对六边形单元蜂窝遮弹层结构进行了侵彻实验研究。实验结果表明,同钢筋混凝土遮弹层相比,蜂窝遮弹层的破坏仅发生在弹靶接触的六边形单元内,且破坏面积较小,同时弹体在侵彻蜂窝遮弹层过程中易发生偏航现象。应用应力波传播理论分析了实验结果,主要是六边形单元对其内的混凝土约束作用及其自身阻隔作用,使混凝土抗压强度和弹体在侵彻过程中受到的阻力增大,从而减小了弹体的破坏效应。     

粘弹层合板的稳态振动和层间应力

应用混合分层理论和Ressiner混合变分原理，在板厚方向取二次位移插值函数和三次、四次横向应力插值函数推导出粘弹层合板的动力学方程，得出简支粘弹层合板稳态振动的解。不仅得出与三层弹性板精确的自振频率吻合良好的解，而且对于粘弹层合板，所计算的自振频率和结构损耗因子也与三维结果吻合较好。计算了自由阻尼层合板对应的低阶法向位移响应幅值和层问横向应力的幅值。结果表明，较高的层间横向正应力是低频稳态振动中引起粘弹层合板分层破坏的主要因素，采用适当模量和厚度的粘弹性材料将有效地降低粘弹层合板的层间横向正应力的幅值。     

基于格子Boltzmann方法的挤出胀大的数值模拟

将单相格子Boltzmann方法(lattice Boltzmann method, LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型．分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合．对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致．表明给出的LBM可以捕捉挤出胀大的瞬态效应．     

钨合金弹体对混凝土靶的超高速侵彻机理

为研究钨合金弹体超高速侵彻混凝土靶的相关机理,构建了适用于超高速撞击的金属强度模型、失效模型和混凝土的本构模型,对93钨合金弹体超高速撞击混凝土靶问题进行了数值模拟。开展了钨合金弹体超高速撞击混凝土靶实验,分析了靶板成坑特性,研究了侵彻总深度和残余弹体长度随撞击速度的变化规律,理论分析了长杆钨弹超高速撞击混凝土的侵彻模型和混凝土靶内的应力波传播。得到以下主要结论:（1）利用金属及混凝土的新本构模型获得的超高速撞击混凝土靶的破坏形貌数值模拟结果与实验结果一致;（2）超高速撞击条件下混凝土靶成坑为“弹坑＋弹洞”形,成坑体积与弹体动能近似成正比;（3）超高速撞击条件下,侵彻深度随弹速提高呈现先增大后减小的现象,高速段侵深降低是弹体经历销蚀侵彻后“刚体侵彻阶段”减少造成的;（4）建立的钨合金超高速撞击混凝土侵彻分析模型,可用来预估侵彻深度、残余弹长、蘑菇头直径等参数;（5）采用建立的超高速撞击混凝土靶内应力波传播理论模型得到的计算结果与实验结果吻合较好。     

特征线法在三维弹／粘塑性波数值解析中的应用

根据广义特征理论建立了求解三维非线性应力波传播问题的特征关系式，并采用特征线法对纵向冲击载荷下弹／粘塑性矩形截面棒的三维应力波传播过程进行了数值模拟。     

不均匀磁电弹板中的厚度扭曲波的传播

压电-压磁复合材料或结构的许多应用是与弹性波的传播密切相关的,这要求人们首先从理论的角度弄清楚弹性波的传播规律。本文研究由多个不均匀磁电弹介质组成的薄板厚度扭曲波的传播性质,从磁电弹全耦合场三维方程出发得到了其精确解,根据所得到的解分析了波的传播特征,这些结果对于理解和设计谐振器、滤波器以及声波元件提供了有价值的理论基础。     

材料的弹粘塑性本构关系和应力波的演化

试验表明，大多数工程材料在冲击载荷作用之下的变形一般都同时包含有可恢复的瞬态性弹性变形和不可恢复的粘滞性塑性变形，即其本构关系可以用弹粘塑性模型来描述。本文从内变量理论出发，探讨了时率相关材料的弹粘塑性本构关系的一般特性，建立了增量型的弹粘塑性本构关系的一般理论框架和普适的表达式，并且对两种最常用的本构模型——Bodner-Partom模型和Johnson-Cook模型给出了在一维应变条件下的具体形式。通过计算和讨论一维应变粘塑性靶板中冲击波的衰减机制和应力波的演化规律，特别是考察各种粘塑性本构模型中的材料参数对冲击波的衰减和应力波的演化的影响，得出了一些可以直接应用或具有一定借鉴价值的结果，为研究应力波的其他衰减机制以及在人防工程中智能防护层设计时新材料的选取奠定了基础。     

用动光弹研究受冲击载荷的阶梯轴中的应力波的传播

本文用动光弹方法研究了受冲击载荷的阶梯轴中的应力波的传播。得到了反映应力波传播全过程的动态条纹图。本文还将截面变化处的入射、反射和透射应力脉冲之比的实验结果与初等理论进行了比较,对动光弹的实验装置进行了较详细的描述,并对试验中的有关问题进行了讨论。     

无限粘弹介质中圆孔时变问题的解析分析

本文推导粘弹介质中圆孔孔径时变时的应力和位移.由粘弹解与弹性解的对应关系得到粘弹时变应力解.用直接解方程法求径向位移,最终归结为求解关于待定函数的l阶非齐次微分方程.将半径时变函数泰勒展开,用幂级数解法得到一般情况下的解.在寻找定解条件时,采用了对待定函数的光滑化处理,认为在t=0的微小邻域内函数仍满足微分方程,通过积分得到与待定系数数目相同的定解条件,从而获得本问题径向位移解析解.对Maxwell粘弹模型的求解证明了该法的可靠性.文中解适用于任意粘弹模型和孔径任意时变的情况.     

45#钢的损伤演化方程和层裂准则研究

通过唯象分析和细观物理统计相结合的方法建立了一种韧性材料的损伤演化方程。在试验结果和内变量理论的基础上得到了45#钢的含损伤热—粘塑性本构关系。用有限差分方法计算了45#钢平板撞击所产生的应力波传播规律、损伤演化规律及层裂进程。通过自由面速度历史的数值模拟,并基于计算结果与试验结果间的最佳一致性,得到了损伤演化方程中的材料参数和极限损伤数值,并以此为依据建立了材料的应力率层裂准则。     

PREDICTION OF THE VISCOELASTIC PROPERTIES OF THE EQUIVALENT PARTICLE FOR THE INTERCALATED MULTI-LAYER STACK OF NANOPLASTICS

The aim of this paper is to apply the asymptotic homogenization method to deter- mining analytically and numerically the transversely isotropic viscoelastic relaxation moduli of the equivalent particle for the intercalated multi-layer stack of intercalated type nanoplastics.A two-phase multilayered material containing n layers is considered.The matrix is assumed to be an isotropic viscoelastic standard linear body and the reinforcement is assumed to be an isotropic elastic body.Final explicit analytical formulae for the effective elastic moduli of the multilay- ered material are derived first;and then the correspondence principle is employed to obtain the homogenized relaxation moduli of the equivalent intercalated particle.A numerical example is given.Final explicit analytical formulae in the time domain derived here make it convenient to estimate the influence of all the particle parameters of micro-structural details on the effective properties of the equivalent intercalated particle.The results of this paper can also be applied to multi-layer composites.     

A combined viscoelastic–viscoplastic behavior of particle reinforced composites

This study formulates a micromechanical model for predicting effective viscoelastic–viscoplastic responses of composites. The studied composites consist of solid spherical particle reinforcements dispersed in a homogeneous matrix. The particle constituent is assumed linear elastic, while the matrix exhibits combined viscoelastic–viscoplastic responses. The Schapery integral model is used for the 3D isotropic non-linear viscoelastic responses. Two viscoplastic models are considered: the Perzyna model, having a rate-independent yield surface and an overstress function, and the Valanis endochronic model based on an irreversible thermodynamics without a yield surface. The Valanis model is suitable for materials when viscoplastic responses occur at early loadings (small stress levels). A unit-cell model with four particle and polymer sub-cells is generated to obtain homogenized responses of the particle-reinforced composites. Available micromechanical models and experimental data in the literature are used to verify the proposed micromechanical model in predicting effective time-dependent and inelastic responses of composites. Field variables in the homogenized composites are compared to the ones in heterogeneous composites. The heterogeneous composites, having detailed particle geometries, are modeled using finite element (FE) method.     

线黏弹性球面发散应力波的频率响应特性

基于线黏弹性球面波Laplace域的理论解, 得到了不同传播距离处粒子速度、粒子位移、应力、应变等力学量的传递函数。以标准线性固体模型为例, 重点讨论了粒子速度频率响应函数的传播特征, 指出随着传播距离的增加, 粒子速度幅频响应函数的高频响应会低于低频响应, 而在理想弹性条件下, 粒子速度幅频响应函数的高频响应一直高于低频响应。以弹性半径为0.025 m的空腔爆炸为例, 采用Laplace数值逆变换方法对粒子速度波形的演化进行了分析, 给出了粒子速度强间断幅值及粒子速度峰值随传播距离变化的衰减规律曲线, 指出黏弹性介质中粒子速度幅值的衰减曲线介于理想弹性介质中粒子速度幅值衰减曲线和黏弹性介质中粒子速度强间断幅值衰减曲线之间。     

The propagation of initially plane waves in nonhomogeneous viscoelastic media

In this study, the propagation of an initially plane wave in a linear isotropic nonhomogeneous viscoelastic medium, where the nonhomogeneity varies transversely to the direction of propagation, is investigated. For this purpose, first the propagation of waves in a linear isotropic viscoelastic medium of arbitrary inhomogeneity is studied by employing the notion of singular surfaces. The characteristic equation governing wave velocities, and the growth and decay equations describing the change of the strength of the discontinuity as the wave front moves are obtained.In the second part of this work, the propagation of initially plane waves is studied for three types of inhomogeneities by employing the findings established in the first part. The first kind of inhomogeneity considered is of axisymmetrical type where the wave propagation velocity depends on the radial coordinate only, increasing linearly up to a certain radial distance and remaining constant thereafter. The second kind is also axisymmetrical with a wave velocity distribution decreasing linearly till a given value of the radial coordinate. In the third one, the wave velocity is assumed to vary linearly over a given interval along a certain coordinate axis only, which is perpendicular to the direction of propagation, and remain constant outside. The ray and wave front analyses are carried out and the decay or growth of stress and velocity discontinuities are studied for each of the three cases.     

Load motion along a beam on a viscoelastic half-space

An expression is derived for equivalent foundation of a viscoelastic half-space interacting with an Euler–Bernoulli beam. It is shown that this equivalent viscoelastic foundation depends on frequencies and wave numbers of the waves in the beam. The real and imaginary part of it substantially varies for phase velocities in between the Rayleigh and shear waves velocities. Radiation of elastic waves occurs for velocities larger than some velocity in that interval. The steady-state beam displacements due to a uniformly moving constant load are calculated for different velocities. The maximum displacement under the load takes place for a velocity of order of the Rayleigh waves velocity.     

Perturbation of a dynamic planar crack moving in a model viscoelastic solid

Weight functions, which give stress intensity factors in terms of applied loading, are constructed, for three-dimensional time-dependent loading of a semi-infinite crack, propagating at uniform speed. Both a model problem, governed by a scalar wave equation, and the full vectorial problem for Mode I loading, are considered. The medium through which the crack propagates is viscoelastic; the approach is general but explicit formulae are given when the medium is a Maxwell fluid. The weight functions are exploited to develop formulae for the first-order perturbations of stress intensity factors when the crack edge is no longer straight but becomes slightly wavy. Implications for stability, and for “crack front waves” in the case of the Mode I problem, are discussed.     

Gravity waves and Rayleigh Taylor instability on a Jeffrey-fluid

A theory for linear surface gravity waves on a semi-infinite layer of viscoelastic fluid described by a Jeffrey model is presented. Results are given for the decay rate and the phase velocity as a function of the parameters of the fluid: a nondimensional time constant, and a ratio of the retardation time to the relaxation time. At small wave numbers the behavior is Newtonian. In other cases depending on the nondimensional parameters, a number of possible other behaviors exist. The influence of the non-dimensional parameters on the growth rate of Rayleigh-Taylor instability is also discussed.     

Propagation of a stress pulse in a viscoelastic rod

Experimental work is reported on the propagation of a stress pulse in a viscoelastic waveguide. The data obtained are compared with results of analysis using one-dimensional wave-propagation theory. The waveguide used in this work is a low-density polyethylene rod 1/2 in. in diameter and 30-in. long. Stress input to the waveguide and the resulting particle velocity at three stations are measured using a crystal stress transducer, two Faraday-principle velocity transducers and a capacitor transducer. The experiment is described mathematically as a boundary-value problem formulated in terms of the one-dimensional equation of motion, the strain-displacement relationship, a hereditary constitutive equation and the stress-boundary condition. Fourier transform and inversion yield an integral expression for velocity which is evaluated numerically at three stations using measured values for the stress-boundary condition, material attenuation and phase velocity. The analytical results compare favorably with the experimental data. The one-dimensional theory appears adequate to describe pulse propagation of this type. The attenuation and phase velocity used here are found to be a linear function and a logarithmic increasing function of frequency respectively.     

孤立波与半无限复合材料体耦合作用研究

基于一维颗粒链中产生的高度非线性孤立波,研究孤立波与半无限复合材料体的耦合作用。根据赫兹定律推导了一维颗粒链中颗粒间相互作用的运动微分方程,建立了颗粒链与半无限复合材料体的接触模型。对于颗粒与复合材料的接触,采用已有文献中修正后的赫兹定律,研究了高度非线性孤立波与半无限复合材料体的耦合力学作用机理,推导了颗粒链与半无限复合材料体的相互耦合运动微分方程组,通过数值计算,得到了各颗粒的内力、速度、位移曲线。分析了材料属性对回弹孤立波出现的时间、幅值的影响。结果表明:随着纤维方向弹性模量的增大,次级回弹波出现的时间和波幅都逐渐增大,随着垂直纤维方向弹性模量的增大,次级回弹波出现的时间先减小后增大,次级回弹波的幅值逐渐减小直至消失。     

Wave Propagation Parallel to the Layers in Elastic or Viscoelastic Layered Composites

Transient waves propagating parallel to the layers in a linear elastic or viscoelastic layered composite are studied. A step load in time is applied at the boundary x = 0 and the head-of-the-pulse asymptotic solution is obtained for large x and large time t. For viscoelastic composites the interaction between the dissipation and the dispersion is controlled by a parameter γ that contains the material mismatch of the layers and the distance: propagated by the waves. As the distance increases, so does γ, and the oscillatory response diminishes. For elastic composites, we show how the oscillatory response depends on the mismatch of the material properties and the thicknesses of the layers. We show that there are composites other than the one with zero mismatch for which the oscillatory response is almost nonexistent.