Analysis of the dynamic radial compaction of granular media

The previously developed continual approximation is used to analyze the radial axisymmetric compaction of a granular medium in the presence of a rigid undeformable rod on the symmetry axis. It is shown that, during pulsed loading, high densities close to those corresponding to the nonporous state can be attained due to inertia effects. The influence of the initial radial dimensions of the rod-powder-medium system on the compaction process is analyzed. The problem is found to be scale invariant under various constraints imposed on the ratio of the characteristic dimensions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 6, pp. 181–189, November–December, 2008     

Effect of Wave Formation during Shock-Wave Compaction of Powders

The problem of shock-wave compaction of a metal powder enclosed in a metal container with a transverse partition is solved. A model of wave formation on the partition and in the compact adjacent to the partition is proposed; the model is based on the loss of strength in the powder due to collapsing of pores and to development of instability of the partition being compressed in the shock wave. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 119–130, January–February, 2006.     

颗粒材料中致密波结构研究

采用一维两相流模型与相应颗粒构形应力函数,研究了致密波的形成及其结构．用简化两相流模型系统地讨论致密波对有关因素的依赖关系．分析指出：小于基体材料音速的致密波仅能在非理想颗粒材料中存在,从波前到波后,所有状态物理量光滑过渡．大于基体材料音速的致密波,波头可能存在间断．应力函数与致密粘性确定后,致密波速度决定致密波结构、宽度、终态压实度．采用一维两相流模型模拟了活塞驱动颗粒床形成致密波这一动态过程．用线方法（MOL）对该方程组求数值解．计算表明,经过短暂的非稳态过程,颗粒床中形成一稳态致密波．分析了活塞速度与初始孔隙率对致密波结构的影响,并对简化两相流模型与两相流模型的计算结果进行了对比．     

Computation of compaction in compressible granular material

Equations modeling compaction in a mixture of granular high explosive and interstitial gas are solved numerically. Both phases are modeled as compressible, viscous fluids. This overcomes well known difficulties associated with computing shock jumps in the inviscid version of the equations, which cannot be posed in a fully conservative form. One-dimensional shock tube and piston-driven compaction solutions compare favorably with experiment and known analytic solutions. A simple two-dimensional extension is presented.     

Computational modeling of the cold compaction of ceramic powders

A finite-element model of the cold compaction of ceramic powders by uniaxial pressing is developed and validated by comparison with experimental data. The mechanical behavior of processing powders is assumed according to the modified Drucker-Prager cap model. The frictional effects and the mechanical behavior of tools involved in the process are taken into account. The proposed model allows evaluation of the density distribution into the processed part, as well as stress and strain fields. Variations of the density distribution due to the unloading and the ejection of the part are evaluated Published in Prikladnaya Mekhanika, Vol. 42, No. 10, pp. 135–143, October 2006.     

Radial expansion of a granular medium in spherical and cylindrical layers

An exact solution that describes the fields of displacements and stresses in an expanding spherical layer is constructed within the framework of the theory of small strains of a granular medium with rigid particles. For finite strains, the problem reduces to a nonlinear system of ordinary differential equations, which is solved by numerical methods. Similar solutions are found in the problem for a cylindrical layer. Based on these solutions, the effect of the dilatancy of the granular medium on the stress-strain state near expanding cavities is found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 190–196, May–June, 2009.     

A finite strain constitutive model for metal powder compaction using a unique and convex single surface yield function

Constitutive modelling of metal powder compaction processes is a challenge in view of realistic simulations. To this end, the article under consideration has two objectives: the first goal is to present a new unique and convex single surface yield function for pressure dependent materials, which is also applicable to other areas of granular materials such as soils or concrete. The flexibility is shown at various materials. The yield function is based on a log-interpolation of two known simple yield functions. A convexity proof of the new yield function is provided. The second objective is to propose a new rate-independent finite strain plasticity model for metal powder compaction, which is based on the multiplicative decomposition of the deformation gradient into an elastic and a plastic part with evolution equations for internal variables representing the basic behaviour of powder materials under compaction conditions. These variables are used for the evolution of the yield function in order to represent the compressible hardening behaviour of powder materials. On the basis of the constitutive model, the material parameters are identified at experimental data of copper powder.     

Fluidization of a granular medium in a viscous fluid under vertical vibration

The dynamics of a granular medium in a cavity filled with incompressible viscous fluid under harmonic vertical vibration are studied experimentally. The sand is fluidized in a relatively thin sublayer of the granular layer near the interface between the media. The fluidization is of the threshold type and is accompanied by intense parametric oscillations of the interface. For viscous fluids, the transition of the sand from a quasi-solid to a fluidized state and the reverse transition associated with a decrease in the oscillation rate occur with hysteresis. The nondimensional governing parameters determining the sand dynamics are established. The analysis is focused on the case of low nondimensional frequencies. Perm’, Paris. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 113–122, May–June, 2000.     

The Effect of Limb Kinematics on the Speed of a Legged Robot on Granular Media

Achieving effective locomotion on diverse terrestrial substrates can require subtle changes of limb kinematics. Biologically inspired legged robots (physical models of organisms) have shown impressive mobility on hard ground but suffer performance loss on unconsolidated granular materials like sand. Because comprehensive limb–ground interaction models are lacking, optimal gaits on complex yielding terrain have been determined empirically. To develop predictive models for legged devices and to provide hypotheses for biological locomotors, we systematically study the performance of SandBot, a small legged robot, on granular media as a function of gait parameters. High performance occurs only in a small region of parameter space. A previously introduced kinematic model of the robot combined with a new anisotropic granular penetration force law predicts the speed. Performance on granular media is maximized when gait parameters utilize solidification features of the granular medium and minimize limb interference.     

Mechanism of compressive stress formation during weak shock waves impact with granular materials

 Certain aspects of wave propagation and the dynamic reaction of a granular material when subjected to a long-duration impulse load are studied. In the majority of studies published on this subject the unsteady pressure behavior at the end-wall covered by a layer of granular material was observed and documented. However, up to now little attention was given to explaining the physical mechanism of this process. Experimental results, obtained in the course of this study, regarding the pressure fields inside granular layers of different materials, clearly show that the compaction effect strongly depends on the characteristics of the medium. This phenomenon manifests itself by changing the gas-particle interaction in the course of the gas filtration, and by variation in the contribution of the different forces and effective stress, σ, to the energy exchange between the gas, the particles and the shock-tube wall. The material permeability,  f, the relative density, ν, and the particle response time, τ _p , are the most important parameters affecting the stress formation at the end-wall covered by the granular layer. In addition to the effect of the material parameters, the effective stress, σ, was found to strongly depend on the granular layer height, h. Based on detailed pressure measurements a qualitative analysis regarding the role of the particle rearrangement in the formation of the unsteady peak at the end-wall was performed. The phenomenology of the particle–particle interaction includes rotation and consolidation of the granules and movement or sliding of particle planes within the layer over each other. Most of these processes are frictional in their nature. They are related to the energy losses and affect the profile and magnitude of the compressive stress as measured at the shock-tube end-wall covered by the granular layer. Received: 10 June 1996/Accepted: 15 October 1996     

公路黄土地基处治前后沉降变形分析

在湿陷性黄土地区,黄土作为高速公路地基的主要原料其强度和变形特性不能满足要求,主要有振动碾压、冲击碾压和强夯3种有效且实用的处治方法。本文根据原位试验和室内试验的数据,进行有限元模拟,在不同高度路堤和不同加载方式情况下,地基沉降变形随着处治能级的增加而减小,路堤本身的沉降保持不变。对于地基处理的深度和沉降变形减小的效果,强夯最好,冲碾次之,振碾最弱。处治前后地表中心沉降变形值与加载路堤高度之间的关系均可用二阶多项式回归拟合。2～4m高的路堤采用振碾的方式处理地基,5～10m高的路堤采用冲碾方法,大于10m的路堤多用强夯或者直接修桥来代替,对于黄土地区高速公路建设与维护具有重要意义。     

振动摩擦机理及其非线性动力学特性

本文揭示了振动工况下散体物料的振动摩擦机理,给出了振动参数与土壤内摩擦力的关系曲线,并用数值试验方法模拟了复杂的振动摩擦压实过程,分别给出了振频和振幅对土壤孔隙比和轴向应力均值的关系曲线,验证了振动压实减摩机理的正确性.结果表明:松散物料在振动工况下,其颗粒间的摩擦力由静摩擦力转变为动摩擦力,随着振频和振幅的增加,摩擦阻力减小,从而土的内摩擦力减小,并且存在最优的振动强度使土壤在振动压实过程中的内摩擦力最小.建立了基于振动压路机的振动摩擦系统非线性动力学模型,利用渐进法得到了振动压实过程的共振响应和非共振响应.     

Cold compaction of ceramic powder: Computational analysis of the effect of pressing method and die shape

This paper deals with a computational analysis of the influence of the pressing method and part geometry on the final density distribution in the cold compaction process of ceramic alumina powders. The analysis is based on the explicit finite-element model proposed and validated in a previous study. The mechanical behavior of the processing material is described using a multisurface elastoplastic model, the modified Drucker-Prager/Cap model Published in Prikladnaya Mekhanika, Vol. 43, No. 10, pp. 129–134, October 2007.     

Sediment-discharge vector in a turbulent flow above an eroded bottom

The problem of determination of sediment discharge by a turbulent flow of a fluid above an eroded surface of an arbitrary relief with a finite slope of the bottom is considered. The surface of the bottom separates a stationary granular medium (sand) from a moving two-phase mixture of a fluid and solid particles. The medium is set into motion under the action of shear stress of the fluid. The medium obeys Coulomb's friction law for a granular medium and Prandtl's law of turbulent friction of the fluid. As a result of solving the boundary-value problem for the motion of a two-phase mixture of a fluid and solid particles, a generic formula for sediment discharges is derived. The sediment-discharge vector is expressed through the vector of shear stress on the bottom, the vector of the slope of the bottom, and the distribution function of the solid particles in the bottom layer for an arbitrary relief of the bottom with a finite slope. It is shown that the sediment discharge depends weakly on the detailed distribution of particles in the bottom layer. Conditions of failure of the bottom surface are obtained. The sediment-discharge formula allows one to derive a closed system of equations that determines the process of bottom erosion in the river or channel bed. Institute Problems of Mechanics, Russian Academy of Sciences, Moscow 117526. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 2, pp. 102–112, March–April, 2000.     

Flux in Porous Media with Memory: Models and Experiments

The classic constitutive equation relating fluid flux to a gradient in potential (pressure head plus gravitational energy) through a porous medium was discovered by Darcy in the mid 1800s. This law states that the flux is proportional to the pressure gradient. However, the passage of the fluid through the porous matrix may cause a local variation of the permeability. For example, the flow may perturb the porous formation by causing particle migration resulting in pore clogging or chemically reacting with the medium to enlarge the pores or diminish the size of the pores. In order to adequately represent these phenomena, we modify the constitutive equations by introducing a memory formalism operating on both the pressure gradient–flux and the pressure–density variations. The memory formalism is then represented with fractional order derivatives. We perform a number of laboratory experiments in uniformly packed columns where a constant pressure is applied on the lower boundary. Both homogeneous and heterogeneous media of different characteristic particle size dimension were employed. The low value assumed by the memory parameters, and in particular by the fractional order, demonstrates that memory is largely influencing the experiments. The data and theory show how mechanical compaction can decrease permeability, and consequently flux.     

Filtration in a Porous Granular Medium: 1. Simulation of Pore-Scale Particle Deposition and Clogging

This paper presents a numerical model for simulating the pore-scale transport and infiltration of dilute suspensions of particles in a granular porous medium under the action of hydrodynamic and gravitational forces. The formulation solves the Stokes’ flow equations for an incompressible fluid using a fixed grid, multigrid finite difference method and an embedded boundary technique for modeling particle–fluid coupling. The analyses simulate a constant flux of the fluid suspension through a cylindrical model pore. Randomly generated particles are collected within the model pore, initially through contact and attachment at the grain surface (pore wall) and later through mounding close to the pore inlet. Simple correlations have been derived from extensive numerical simulations in order to estimate the volume of filtered particles that accumulate in the pore and the differential pressure needed to maintain a constant flux through the pore. The results show that particle collection efficiency is correlated with the Stokes’ settling velocity and indirectly through the attachment probability with the particle–grain surface roughness. The differential pressure is correlated directly with the maximum mound height and indirectly with particle size and settling velocity that affect mound packing density. Simple modification factors are introduced to account for pore length and dip angle. These parameters are used to characterize pore-scale infiltration processes within larger scale network models of particle transport in granular porous media in a companion paper. Articlenote: Currently at GZA GeoEnvironmental Inc., 1 Edgewater Drive, Norwood, MA 02062, U.S.A.     

颗粒物质对风积沙压实特性影响研究

本文从颗粒物质的角度探讨了风积沙的压实特性和机理。研究表明,风积沙的重型击实试验曲线呈"双驼峰"形式,表明风积沙具有干压实的特性,水在风积沙的压实过程中扮演着不同的作用|室内振动试验表明振动加速度及振动时间影响着风积沙"巴西坚果效应"与"反巴西坚果效应"之间的转化,从而影响风积沙的振实密度,二者均存在一个产生最大振实干密度的最优值。现场试验通过控制机械参数(振幅A=0.4mm,频率f=48Hz)在压实5遍的时候,可使压实度达到96.8%,CBR值达到25.3,均能满足《公路路基设计规范》要求。     

Analysis of stability of a thin layer of granular material moving on an inclined plane

The stability of a layer of a granular medium on an inclined plane has been studied within the framework of the model of a non-Newtonian fluid with an index of 2, which ensures the experimentally found quadratic dependence of the shear stress on the shear rate. It is shown analytically and numerically that these flows are stable or unstable depending on the value of the generalized Reynolds number relative to the critical value equal to5 cot α. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 6, pp. 113–117, November–December, 1998.     

Theoretical model of acoustic emission under mechanical loading of rocks in the maximum compaction region

A theoretical model of changes in acoustic emission activity in a geomaterial under continuous or stepwise mechanical loading is justified. Based on this model, the experimentally found laws of emission in the region of the maximum compaction of rock samples with different rates of mechanical loading of these samples are analyzed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 4, pp. 145–152, July–August, 2006.