黄土和砂土岩中填实爆炸辐射弹性波的对比研究

研究地下爆炸弹性区的震动特性,关键是获得场地介质与爆炸能量耦合作用下辐射弹性波的实验参数。对于不易加工成大尺寸模型的砂土岩,为研究其填实爆炸下辐射弹性波的特征,采用0.125 g TNT微型炸药球作为爆炸源,以塑性区可置换的?1 370 mm×1 200 mm黄土样品作为提供应力波传播路径的载体,用波阻抗近似相等的重塑黄土和砂土岩样品分别作为源区介质,对比分析了两种介质中微药量填实爆炸辐射的弹性波传播特征。实验结果表明:在测试范围内,两种介质中填实爆炸激发的弹性应力波粒子速度(位移)峰值的衰减规律、波形的主频变化规律一致;砂土岩中爆炸辐射的弹性波粒子速度(位移)峰值整体高于黄土、粒子速度波形的半高宽和主频低于黄土;砂土岩中爆炸耦合的向外传播的弹性波能量比黄土大。实测结果反映,黄土和砂土岩中填实爆炸弹性波能量耦合强度的差别。     

The constitutive equation of a dilute suspension of spherical microcapsules

This paper studies the rheological behavior of a dilute suspension of spherical microcapsules, i.e. spherical thin elastic membranes filled with an incompressible liquid. Previous results obtained for the motion of such capsules freely suspended in a simple shear flow are first extended to general linear shear flows. Then the stresslet term in the external flow is computed to order ?², where ?, assumed to be small with respect to unity, is the ratio of viscous to elastic forces acting on the particle. The resulting constitutive equation is of the viscoelastic type and is similar to the one obtained for liquid droplets. It predicts that the microcapsule suspension exhibits a shear dependent viscosity and normal stress effects. The exact dependency of these phenomena on the microscopic parameters of the suspension is explicitly provided by the model.     

Multiple scattering of elastic waves in metal-matrix composite materials with high volume concentration of particles

A new approach is proposed to investigate the propagation of compressional (P) and shear (SV) waves in metal-matrix composite materials with high volume concentration of particles. The theory of quasicrystalline approximation and Waterman's T matrix formalism are employed to treat the multiple scattering resulting from the particles in composites. The addition theorem for spherical Bessel functions is used to accomplish the translation between different coordinate systems. The analytical expression of the Percus–Yevick correlation function is also given. Closed form solutions for the effective propagation constants and the dynamic effective elastic modulus of materials are obtained in the low frequency limit. At higher frequencies, only numerical results of them are presented. Numerical examples show that the phase velocities of P and SV waves in the composite materials with low volume concentration in the low frequency are in good agreement with the results in previous literatures. The effects of the incident wave number, the volume fraction of particles and the material properties of the particles and matrix on the dynamic effective elastic modulus are also examined.     

Constitutive model for predicting dynamic interactions between soil ejecta and structural panels

A constitutive model is developed for the high-rate deformation of an aggregate comprising of mono-sized spherical particles with a view to developing an understanding of dynamic soil-structure interactions in landmine explosions. The constitutive model accounts for two regimes of behaviour. When the particle assembly is widely dispersed (regime I), the contacts between particles are treated as collisions, analogous to those between molecules in a gas or liquid. At high packing densities (regime II) the contacts are semi-permanent and consolidation is dominated by particle deformation and inter-particle friction. Regime I is modelled by extending an approach proposed by Bagnold (1954. Experiments on a gravity-free dispersion of large solid particles in a Newtonian fluid under shear. Proceedings of the Royal Society of London A 225, 49-63) to a general strain history comprising volumetric and deviatoric deformation. For regime II, classical soil mechanics models (such as Drucker-Prager) are employed. The overall model is employed to investigate the one-dimensional impact of sand against a rigid stationary target. The calculations illustrate that, unlike single-particle impact, the momentum transmitted to a rigid target is insensitive to the particle co-efficient of restitution, but strongly dependent on initial density. The constitutive model is also used to examine the spherical expansion of a shell of sand (both dry and water saturated). In line with initial experimental observations, the wet sand is predicted to form clumps while the dry sand fully disperses. The model shows that this clumping of explosively loaded wet sand exerts higher pressures on nearby targets compared to equivalent dry sand explosions.     

Mechanism of shear attenuation near the interface under combined compression and shear impact loading

We investigate the compressional/shear coupling plastic wave propagation characteristics analytically for ideal elastic–plastic materials in both stress and particle velocity spaces, focusing on the shear wave attenuation near the interface occurring in pressure–shear plate impact experiments. The results show that the shear attenuation is strongly associated with the wave propagation characteristics of the coupling waves. In the stress space, as the shear stress increases, an adjustment of the stress components is observed and the final stress state along the wave path is a combined pure shear- and hydrostatic pressure-state. In the particle velocity space, the wave structures with different loading and maximal transverse particle velocity are obtained. The maximal transverse particle velocity varies with the longitudinal velocity and forms a boundary line. Once the loading transverse velocity exceeds this line, a transverse particle velocity discontinuity occurs at the impact interface. If the bonding strength is sufficiently high, there will be a shear band in the target in the extreme vicinity of the interface.     

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

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

     

Mass and heat transfer during solid sphere dissolution in a non-uniform fluid flow

We studied a nonisothermal dissolution of a solvable solid spherical particle in an axisymmetric non-uniform fluid flow when the concentration level of the solute in the solvent is finite (finite dilution of solute approximation). It is shown that simultaneous heat and mass transfer during solid sphere dissolution in a uniform fluid flow, axisymmetric shear flow, shear-translational flow and flow with a parabolic velocity profile can be described by a system of generalized equations of convective diffusion and energy. Solutions of diffusion and energy equations are obtained in an exact analytical form. Using a general solution the asymptotic solutions for heat and mass transfer problem during spherical solid particle dissolution in a uniform fluid flow, axisymmetric shear flow, shear-translational flow and flow with parabolic velocity profile are derived. Theoretical results are in compliance with the available experimental data on falling urea particles dissolution in water and for solid sphere dissolution in a shear flow.     

基于ZWT方程的线黏弹性球面波分析

针对时间尺度为1~100 s的爆炸冲击加载,忽略低频Maxwell体的松弛效应,同时为简化分析,忽略非线性弹簧效应项,推导了三维应力状态下的线黏弹性ZWT本构关系。基于球面波的基本动力学方程,结合ZWT线黏弹性本构关系,得到了以位移u表征的三阶波动方程。利用该方程分析固体介质中线黏弹性球面波传播过程中的吸收和弥散现象,得知:高频球面波的衰减因子趋于常数,相速趋于高频下的纵波速度;低频球面波的衰减因子和圆频率的平方成正比,其相速趋于低频下的纵波速度;低频球面波的纵波波速低于高频球面波的纵波波速,两者的比值和介质的泊松比、弹性模量及Maxwell体弹性模量相关。     

水击驻波场中分散相颗粒的运动分析

根据水击在管段内形成的驻波场现象,分析了流体内分散相颗粒受到的驻波作用力;运用李雅普诺夫稳定判据研究了颗粒积聚与分离的机理;考虑到颗粒运动方程的严重刚性而很难进行数值求解,采用相空间和非对称分析方法获得了分散相颗粒的运动轨迹近似解,并进行了实验验证。结果表明:水击驻波场中分散相颗粒的受力方程中惯性项对颗粒初始运动速率的影响不可忽略;在水击驻波波节的±λ/4范围内,分散相颗粒经过一定的时间会发生积聚,其运动速度呈对称分布,最大速度出现在3λ/8位置处;随着分散相颗粒粒径和密度等物性参数以及水击驻波的频率和连续相初始速度的增大,颗粒达到平衡位置的时间呈减小趋势,且连续相的初始速度对颗粒到达波节时间的影响显著。     

Rotation of dispersed particles in a vortex field

An expression is obtained for the angular velocity of a spherical dispersed particle in a viscous fluid in an external vortex field with an harmonic time dependence. This expression is then used for investigating a system of two rotating dispersed particles whose rotation is the result of the interaction of the particles in the field of an incident sound wave. It is found that such a system possesses a rather interesting nontrivial property: under certain conditions it has a resonant frequency at which the rotation of the particles relative to the fluid is most intense.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 186–188, July–August, 1992.     

Microscale Numerical Simulation of the Permeability Reduction due to Trapping of Suspended Fine Particles Within Sand Sediments

The reduction in permeability of sediments due to blockages caused by the trapping of suspended particles is a common concern for the extraction processes of oil or natural gas. In this study, the effect of trapped fine particles in sand sediments is studied numerically using a three-dimensional lattice Boltzmann method. The geometrical properties of larger, immobile, sand grains are digitally extracted by the spherical harmonics series expansions of CT scans of real sand grains. The migrating fine particles are assumed to be spherical in shape with their volumes following a log-normal distribution. These fine particles, together with larger frame sands, are positioned, without overlapping, within a microscopic, cubic, domain with periodic boundaries. The remaining empty volume is filled with water and imposing a pressure gradient simulates the flow of fluid through the sediment. As a result of fine particles becoming trapped by the frame sand, the initial porosity of which is 0.589, the absolute permeability of the system is reduced by approximately 60?C90?%, corresponding to fine particle saturations of 0.15?C0.29, respectively. The permeability change due to the trapping of fine particles is also modelled theoretically using not only volume saturations but also specific surface areas of both the frame sands and the fine particles with a coefficient of proportionality.     

Extinction of elastic wave energy due to scattering in a viscoelastic medium

Forward scattering theorem for elastic longitudinal and shear wave scatterings by an arbitrary-shaped three-dimensional object embedded in a viscoelastic medium is derived. It is shown that the formulae for extinction cross-sections of an object in an energy-absorbing medium are formally the same with those of the object in the lossless elastic medium. Numerical calculations are executed for the longitudinal wave scattering in an epoxy matrix by a spherical inclusion with different material properties. The condition of negative extinction is examined with the causality constraint on the viscoelastic medium taken into account. It is found that the negative extinction occurs in the Rayleigh limit when the attenuation of the medium is sufficiently high and, more restrictedly, the wave speed in the object is larger than that in the medium, while it occurs less likely in the high frequency range considered in this paper (0<ka<100).     

Experimental study of planar shock wave interactions with dense packed sand wall

A dense packed sand wall is impacted by a planar shock wave in a horizontal shock tube to study the shock-sand wall interaction. The incident shock Mach number ranges from 2.18 to 2.38. A novel device for actively rupturing diaphragm is designed for the driver section of the shock tube. An apparatus for loading particles is machined by the electrical discharge cutting technique to create a dense packed particle wall. High-speed schlieren imaging system and synchronized pressure measurement system are used together to capture the wave structures and particle cloud velocity. The dynamic evolution model from dense packed particles to dense gas–solid cloud at the initial driving stage is established. The blockage and permeation effects of the sand wall work together and influence each other. The high pressure gas behind the incident shock wave blocked by the sand wall pushes the upstream front of the wall forward like a piston. Meanwhile, the high speed gas permeating through the sand wall drags the sands of the most downstream layer forward. The incident shock strength, initial sand wall thickness and particle diameter are varied respectively to investigate the shock attenuation and the wall acceleration. Increasing the sands diameter or mixing in small diameter sands can significantly attenuate the incident shock. The smaller particles or the particles in thinner wall can be dispersed into a larger range in the process of transform from dense packed particles to dense gas–solid cloud. Moreover, the stronger incident shock can disperse the particles into a larger region.     

基于IS–FEM的颗粒曳力系数计算与试验验证

基于浸没光滑有限元模型（immersed smooth finite element model,IS–FEM）,计算球形与非球形颗粒曳力系数。设计颗粒曳力图像测量试验,验证IS–FEM模拟精度。颗粒相的运动行为基于连续介质理论的光滑有限元法求解;流体控制方程通过特征分解的半隐式有限元法求解;颗粒与流体相通过非贴体网格交换数据。结果表明,球形颗粒流场特征对称分布,非球形颗粒稳定沉降时长轴与重力方向垂直。不同雷诺数下球形颗粒的曳力系数计算值与Stokes曳力系数一致,非球形颗粒曳力系数高于等效球形颗粒,IS–FEM计算值与沉降试验吻合良好。

     

Non-linear modulation of SH waves in a two-layered plate and formation of surface SH waves

Non-linear modulation of shear horizontal (SH) waves in a two-layered elastic plate of uniform thickness is considered. Both layers are assumed to be homogeneous, isotropic and incompressible elastic and having different mechanical properties. The problem is investigated by a perturbation method and in the analysis it is assumed that between the linear shear velocities of the top layer, c₁, and the bottom layer, c₂, the inequality c₁<c₂ is valid. In the layered structure then an SH wave exists if the wave velocity c of the wave satisfies either the condition c₁<c?c₂ or the one c₁<c₂?c. Here the problem is examined under the former condition and it is shown that the non-linear modulation of SH waves is governed by a non-linear Schrödinger equation. In this case the formation of surface SH (Love) waves is also revealed if the top layer is thinner when compared with the bottom layer. Then the stability condition is discussed and the existence of bright (envelope) and dark solitons are manifested.     

Discrete element modeling of acoustic emission in rock fracture

The acoustic emission (AE) features in rock fracture are simulated numerically with discrete element model (DEM). The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.     

近壁泊肃叶流中活性粒子迁移规律研究

研究活性粒子在剪切流中的迁移规律对实现颗粒分离和过程强化均具有重要意义.基于耗散粒子动力学理论,建立了描述微通道内近壁泊肃叶流中活性粒子迁移运动的数学模型,考察了活性粒子圆周运动角速度、手性诱导角速度、直行运动速度和转向扩散系数对大肠杆菌和常规活性粒子横向迁移速度和受迫转向频率的影响规律,并确定近壁剪切流中活性粒子横向迁移的形成机制.结果表明,近壁剪切流场中大肠杆菌的横向迁移速度随剪切速率增大先快速增加继而趋于稳定;大肠杆菌横向迁移速度随圆周运动角速度增大而减小,随手性诱导角速度、直行运动速度和转向扩散系数的增大而增大;大肠杆菌的受迫转向频率受圆周运动角速度、直行运动速度和转向扩散系数的影响小,而随手性诱导角速度的增大而加快;相比大肠杆菌,常规活性粒子横向迁移速度显著减小、受迫转向频率明显变慢,二者受直行运动速度和转向扩散系数的影响规律与大肠杆菌类似.直行运动是活性粒子形成横向迁移运动的前提,其他运动参数和结构参数均可一定程度促进或抑制活性粒子在近壁剪切流场中的横向迁移.     

Love wave propagation in heterogeneous micropolar media

The present paper framed to study the impact of heterogeneity on propagation of Love wave in a heterogeneous micropolar layer over an elastic inhomogeneous stratum, when both rigidity and density are assumed to vary linearly with depth. The equations of motion have been formulated separately for layer and half-space under suitable boundary conditions. Analytical solution for the dispersion equation has been obtained using method of separation of variables by means of the Airy function and Whittaker function. Some particular cases have also been investigated. Further, as a special case the velocity equation for isotropic layer over a homogeneous half-space coincides with the standard result of Love wave. Numerical calculations of frequency relation have been performed and depicted by means of graphs to exhibit the substantial impact of heterogeneity, micropolar parameters and wave number on the phase velocity of Love wave. The wave velocity is strongly influenced by these parameters.     

Hydromagnetic Stokes flow in a rotating fluid with suspended small particles

An initial value investigation is made of the motion of an incompressible, viscous conducting fluid with embedded small spherical particles bounded by an infinite rigid non-conducting plate. Both the plate and the fluid are in a state of solid body rotation with constant angular velocity about an axis normal to the plate. The flow is generated in the fluid-particle system due to non-torsional oscillations of a given frequency superimposed on the plate in the presence of a transverse magnetic field. The operational method is used to derive exact solutions for the fluid and the particle velocities, and the wall shear stress. The small and the large time behaviour of the solutions is discussed in some detail. The ultimate steady-state solutions and the structure of the associated boundary layers are determined with physical implications. It is shown that rotation and magnetic field affect the motion of the fluid relatively earlier than that of the particles when the time is small. The motion for large times is set up through inertial oscillations of frequency equal to twice the angular velocity of rotation. The ultimate boundary layers are established through inertial oscillations. The shear stress at the plate is calculated for all values of the frequency parameter. The small and large-time behaviour of the shear stress is discussed. The exact solutions for the velocity of fluid and the wall shear stress are evaluated numerically for the case of an impulsively moved plate. It is found that the drag and the lateral stress on the plate fluctuate during the non-equilibrium process of relaxation if the rotation is large. The present analysis is very general in the sense that many known results in various configurations are found to follow as special cases.     

Impact comminution of solids due to local kinetic energy of high shear strain rate: I. Continuum theory and turbulence analogy

The modeling of high velocity impact into brittle or quasibrittle solids is hampered by the unavailability of a constitutive model capturing the effects of material comminution into very fine particles. The present objective is to develop such a model, usable in finite element programs. The comminution at very high strain rates can dissipate a large portion of the kinetic energy of an impacting missile. The spatial derivative of the energy dissipated by comminution gives a force resisting the penetration, which is superposed on the nodal forces obtained from the static constitutive model in a finite element program. The present theory is inspired partly by Grady's model for expansive comminution due to explosion inside a hollow sphere, and partly by analogy with turbulence. In high velocity turbulent flow, the energy dissipation rate gets enhanced by the formation of micro-vortices (eddies) which dissipate energy by viscous shear stress. Similarly, here it is assumed that the energy dissipation at fast deformation of a confined solid gets enhanced by the release of kinetic energy of the motion associated with a high-rate shear strain of forming particles. For simplicity, the shape of these particles in the plane of maximum shear rate is considered to be regular hexagons. The particle sizes are assumed to be distributed according to the Schuhmann power law. The condition that the rate of release of the local kinetic energy must be equal to the interface fracture energy yields a relation between the particle size, the shear strain rate, the fracture energy and the mass density. As one experimental justification, the present theory agrees with Grady's empirical observation that, in impact events, the average particle size is proportional to the (−2/3) power of the shear strain rate. The main characteristic of the comminution process is a dimensionless number B_a (Eq. (37)) representing the ratio of the local kinetic energy of shear strain rate to the maximum possible strain energy that can be stored in the same volume of material. It is shown that the kinetic energy release is proportional to the (2/3)-power of the shear strain rate, and that the dynamic comminution creates an apparent material viscosity inversely proportional to the (1/3)-power of that rate. After comminution, the interface fracture energy takes the role of interface friction, and it is pointed out that if the friction depends on the slip rate the aforementioned exponents would change. The effect of dynamic comminution can simply be taken into account by introducing the apparent viscosity into the material constitutive model, which is what is implemented in the paper that follows.