Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

DEM simulation of particle flow and separation in a vibrating flip-flow screen

Vibrating flip-flow screens (VFFS) with stretchable polyurethane sieve mats have been widely used in screening fine-grained materials in recent years. In this work, the discrete element method (DEM) is used to study the screening process in VFFS to explain particle flow and separation behavior at the particle scale. Unlike traditional vibrating screens, for VFFS, the amplitude response of each point on the elastic sieve mat is different everywhere. This study measures the kinematics of the elastic sieve mat under different conditions such as different stretched lengths and material loads. To establish the elastic sieve mat model in a DEM simulation, the continuous elastic sieve mat is discretized into multiple units, and the displacement signal of each unit tested is analyzed by Fourier series. The Fourier series analysis results of each unit are used as the setting parameters for motion. In this way, the movement of the elastic sieve mat is approximately simulated, and a DEM model of VFFS is produced. Through the simulation, the flow and separation of different-sized particles in VFFS are studied, and the reasonability of the simulation is verified by a pilot-scale screening experiment. The present study demonstrates the potential of the DEM method for the analysis of screening processes in VFFS.     

Soil flow analysis for grouser wheels based on a particle image velocimetry method

This paper presents an analysis, based on a particle image velocimetry method, of soil flow field beneath a grouser wheel traveling over loose soil. Although the grouser wheel is expected to have better traction and mobility over fine, loose soil, its interaction mechanisms with the soil remain to be elucidated. Thus, a particle image velocimetry-based soil flow analysis is conducted to directly observe soil behavior around the grouser wheel. In the experimental analysis, key parameters of the soil flow field, such as general shape, thickness, streamlines of the flow field, soil velocity on the streamlines, and soil failure angle are examined quantitatively. From the results, the soil flow shape periodically changes with wheel rotation, and this change appears, depending on wheel slip varying over time. Furthermore, the experimental result of the soil failure angle differs drastically from its typical theory. These results will contribute to modeling the mechanical interaction between the grouser wheel and soil.     

被动桩中土拱效应特征与影响参数研究

采用平面有限元方法,对被动桩在粘性、无粘性土体条件下的土拱形成机理进行了分析。从土拱效应的4个主要方面:土拱形状、桩周土体塑性(拉裂)区的分布、竖向位移等值线及桩后土体残余荷载分担比,探讨了改变桩、土参数,如粘聚力、内摩擦角、弹性模量、泊松比、剪胀角、桩土接触面特性时,土拱效应的变化规律。研究表明,土性不同时,土拱形成过程中土体的受力和破坏形式不同;泊松比、剪胀角、桩土接触特征对土拱效应的影响最明显;强度高、剪胀角大、桩土接触面粗糙或泊松比小的土体,易形成土拱。根据研究结果,指出Tom io Ito排桩塑性绕流土压力理论计算公式尚存不足。     

The crack tip zone shielding effect for ductile particle reinforced brittle materials

The crack tip zone shielding effect for the ductile particle reinforced brittle materials is analyzed by using a micromechanics constitutive theory. The theory is developed here to determine the elastoplastic constitutive behavior of the composite. The elastoplastic particles, with isotropic or kinematical hardening, are uniformly dispersed in the brittle elastic matrix. The method proposed is based on the Mori-Tanaka's concept of average stress in the composite. The macroscopic yielding condition and the incremental stress strain relation of the composite during plastic deformation are explicity given in terms of the macroscopioc applied stress and the microstructural parameters of the composite such as the volume fraction and yield stress of ductile particles, elastic constants of the two phases, etc. Finally, the contribution of the plastic deformation in the particles near a crack tip to the toughening of the composite is evaluated. The project supported by National Natural Science Foundation of China     

Micromechanics of breakage in sharp-edge particles using combined DEM and FEM

By combining DEM （Discrete Element Method） and FEM （Finite Element Method）, a model is established to simulate the breakage of twodimensional sharp-edge particles, in which the simulated particles are assumed to have no cracks. Particles can, however, crush during different stages of the numerical analysis, if stress-based breakage criteria are fulfilled inside the particles. With this model, it is possible to study the influence of particle breakage on macro- and micro-mechanical behavior of simulated angular materials. Two series of tests, with and without breakable particles, are simulated under different confining pressures based on conditions of biaxial tests. The results, presented in terms of micromechanical behavior for different confining pressures, are compared with macroparameters. The influence of particle breakage on microstructure of sharp-edge materials is discussed and the related confining pressure effects are investigated. Breakage of particles in rockfill materials are shown to reduce the anisotropy coefficients of the samples and therefore their strength and dilation behaviors.     

Micromechanics of breakage in sharp-edge particles using combined DEM and FEM

By combining DEM (Discrete Element Method) and FEM (Finite Element Method),a model is established to simulate the breakage of two-dimensional sharp-edge particles,in which the simulated particles are assumed to have no cracks.Particles can,however,crush during different stages of the numerical analysis,if stress-based breakage criteria are fulfilled inside the particles.With this model,it is possible to study the influence of particle breakage on macro- and micro-mechanical behavior of simulated angular materials.Two series of tests,with and without breakable particles,are simulated under different confining pressures based on conditions of biaxial tests.The results,presented in terms of micromechanical behavior for different confining pressures,are compared with macroparameters.The influence of particle breakage on microstructure of sharp-edge materials is discussed and the related confining pressure effects are investigated.Breakage of particles in rockfill materials are shown to reduce the anisotropy coefficients of the samples and therefore their strength and dilation behaviors.     

多孔介质的流变模型研究

多孔介质在应力作用下具有弹性变形和黏性变形两种完全不同的变形机制，多孔介质的弹性变形是由介质的本体有效力所致，而黏性变表则是由介质的结构有效应力所致。多孔介质的总变形为弹性变形和黏性变形的叠加，计算多孔介质总应变量的流变模型必须同时采用本体有效应力和结构有效应力（双重有效应力），而传统的流变模型仅采用Terzaghi有效应力是不妥当的，它无法正确描述多孔的应变行为，采用了双重有效应力之后的流变模型，通过调节介质特性参数，可以拟合介质的实际应变行为，并且把多孔介质与普通固体联系了起来。     

基于持续同调的边坡土颗粒应力链失稳分析

采用颗粒离散元方法和持续同调理论,研究了内排土场堆叠至不同高度时的边坡稳定性。为便于研究,现采用一水平金属板向下施加压力,代替不同厚度土层的重力荷载,对边坡在竖向荷载作用下的失稳破坏过程进行了颗粒离散元模拟。研究了二维边坡土颗粒速度总矢量、边坡失稳破坏时滑移开裂面的角度以及边坡坡顶y方向的平均速度等宏观响应过程,并构建了自然堆积下边坡堆积体颗粒的法向力链无向网络模型。最后,用持续同调方法对边坡坡顶颗粒接触力链网络的拓扑特征进行分析,获得条码图,建立了岩体结构持续同调特征与失稳演化的关系。本文为研究边坡失稳拓扑识别提供了一种新方法,从而可以有效预测边坡的失稳破坏。     

Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression

Based on experimental results of brittle, intact sandstone under uniaxial compression, the micro-parameters were firstly confirmed by adopting particle flow code(PFC2D). Then, the validation of the simulated models were cross checked with the experimental results of brittle sandstone containing three parallel fissures under uniaxial compression. The simulated results agreed very well with the experimental results, including the peak strength, peak axial strain, and ultimate failure mode. Using the same microparameters, the numerical models containing a new geometry of three fissures are constructed to investigate the fissure angle on the fracture mechanical behavior of brittle sandstone under uniaxial compression. The strength and deformation parameters of brittle sandstone containing new three fissures are dependent to the fissure angle. With the increase of the fissure angle, the elastic modulus, the crack damage threshold,and the peak strength of brittle sandstone containing three fissures firstly increase and secondly decrease. But the peak axial strain is nonlinearly related to the fissure angle. In the entire process of deformation, the crack initiation and propagation behavior of brittle sandstone containing three fissures under uniaxial compression are investigated with respect tothe fissure angle. Six different crack coalescence modes are identified for brittle sandstone containing three fissures under uniaxial compression. The influence of the fissure angle on the length of crack propagation and crack coalescence stress is evaluated. These investigated conclusions are very important for ensuring the stability and safety of rock engineering with intermittent structures.     

Comparative study of different discrete element models and evaluation of equivalent micromechanical parameters

Comparative studies of different discrete element models of a rock-type material are presented. The discrete element formulation employs spherical particles with the cohesive interaction model combining linear elastic behaviour with brittle failure. Numerical studies consisted in simulation of the uniaxial compression test. Two cylindrical specimens with particle size distributions yielding different degree of heterogeneity have been used. Macroscopic response produced by different discrete element models has been compared. The main difference between the compared models consists in the evaluation of micromechanical constitutive parameters. Two approaches are compared. In the first approach, the contact stiffness and strength parameters depend on the local particle size, while in the second approach, global uniform contact parameters are assumed for all the contacting pairs in function of average geometric measures characterizing the particle assembly. The size dependent contact parameters are calculated as functions of geometric parameters characterizing each contacting particle pair. As geometric scaling parameters, the arithmetic and harmonic means, as well as the minimum of the radii of two contacting particles are considered. Two different models with size dependent contact parameters are formulated. The performance of these models is compared with that of the discrete element model with global uniform contact parameters. Equivalence between the models with size dependent and uniform contact parameters has been checked. In search of this equivalence, different methods of evaluation of global uniform parameters have been studied. The contact stiffness has been evaluated in terms of the average radius of the particle assembly or in terms of the averages of the arithmetic and harmonic means of the contact pair radii, the geometric parameters used in the evaluation of the contact stiffness in the size-dependent models. The uniform contact strengths have been determined as functions of the averages of radii squares, squares of arithmetic radii means or squares of minimum radii of the contacting pairs.For the more homogenous specimen, the models with local size dependent parameters and models with global uniform parameters give similar response. The models with uniform parameters evaluated according to the averages of the geometric parameters used in the evaluation of local parameters ensure better agreement with the respective models with size-dependent parameters than the models with uniform parameters evaluated according to the particle radii. Simulations using the more heterogenous specimen reveal differences between the considered models. There are significant differences in stress–strain curves as well as in the failure pattern. The models with local size-dependent parameters are more sensitive to the change of heterogeneity than the model with global uniform parameters.     

Multiscale approaches to describe mechanical responses induced by particle removal in granular materials

A multi-scale approach is proposed as an attempt to assess internal erosion induced effects on the mechanical properties of a granular medium. When submitted to internal flow, some particles tend to be removed from the initial granular material, leading to drastic changes in the microstructure. The medium is changing during its lifetime and this cannot be ignored in attempts to model its mechanical behavior. In this first analysis, the degradation of a granular assembly is simulated through the progressive removal of its finest particles for both isotropic and anisotropic stress states. A discrete element model as well as an analytical micromechanical model are used to compare induced deformations and properties changes along the erosion process. The shear strength and flow properties appear to be strongly modified by this extraction phase: both models highlight a change from a dilatant to a contractant behavior with the degradation as well as strong dependency of the overall stability on the mobilized friction level. In particular, results show how failure is triggered when particle removal occurs for mobilized frictions greater than the one reached at the critical state.     

Analysis of soil flow and traction mechanics for lunar rovers over different types of soils using particle image velocimetry

Grouser wheels have been used in planetary rovers to improve mobility performance on sandy terrains. The biggest difference between a wheel with and without grousers is the soil behavior beneath the wheel as the grousers shovel the soil. By analyzing the soil flow, we gain insight into the mechanics dominating the interaction between the wheel and the soil, directly impacting performance. As the soil flow varies depending on the soil properties, the effects of soil type on soil behavior and wheel-traveling performance should be studied. This paper reveals the difference in soil flow and wheel performance on cohesive and non-cohesive soils. We conducted a series of single wheel tests over different types of soils under several wheel-traveling conditions. Soil flow was visualized by using particle image velocimetry (PIV). The experimental results indicate that soil flow characteristics highly depend on the shear strength of the soil. The cohesive soil exhibited lower fluidity due to its higher shear strength. At the same time, the wheel displayed a higher traveling performance over the cohesive soil, that is, a lower slip ratio.     

一种土的非线性弹性本构模型参数的反演方法

旨在提出一种土的非线性弹性本构模型参数反演的方法。以现今普遍实行的地基载荷试验为基础,依据遗传算法的组合优化理论,采用正演计算和遗传算法优化相结合的方式,建立了土层非线性弹性本构模型参数反演的方法;并依据某黄土场地地基载荷试验数据,实施了黄土土层非线性弹性本构模型参数反演的全过程。计算结果表明,所建立的方法可以实现土层非线性弹性本构模型中相互关联的多个参数的组合优化,并在对初始值要求较低的情况下,可以获得良好的参数反演精度。从而为土的变形特性分析和土与其中及相邻结构的共同作用分析,提供了较好的土体本构模型参数的确定方法。     

粒子入射角度对冷喷涂涂层形成的影响

采用有限元数值计算方法，研究冷喷涂材料改性过程中铜粒子与铜基体非垂直碰撞的变形行为。针对单个粒子相同入射速度不同入射角度的碰撞条件，探讨粒子与基体的结合强度、侵彻深度以及绝热剪切失稳的发生条件。结果表明:随着粒子入射角度的增大，侵彻深度逐渐减少，粒子与基体的结合强度逐渐减弱。发生绝热剪切失稳的条件是入射速度的法向分量大于碰撞过程的临界速度。     

Role of particle stiffness and inter-particle sliding friction in milling of particles

Discrete element method （DEM） has been used to investigate the effects of particle elastic modulus and coefficient of inter-particle sliding friction on milling of mineral particles. An autogeneous mill of 600 mm diameter and 320 mm length with 14,500 particles has been selected for the simulation. Various mill performance parameters, for example, particle trajectories, collision frequency, collision energy and mill power have been evaluated to understand the effects of particle elastic modulus and inter-particle sliding friction during milling of particles. For the given model, it has been concluded that at high energy range, as the elastic modulus and particle sliding friction increase the energy dissipated among the particles increases. The collision frequency increases with the increase in elastic modulus, however, this trend is not clearly observed with increasing inter-particle sliding friction. The power draw of the mill increases with the increase in fraction of mill critical speed.     

非饱和土力学理论的研究进展

回顾了非饱和土有效应力的发展,目前普遍认同采用两个应力变量来建立本构模型,且对基质吸力中毛细和粘吸两部分作用进行了阐述。分析了非饱和土强度问题,包括抗剪强度和抗拉强度。讨论了非饱和土的本构模型问题,包括基于净应力和基质吸力的弹塑性模型,基于Bishop有效应力和基质吸力的水力力学耦合弹塑性模型,以及双孔隙结构的模型。最后探讨了热力学方法和多孔介质理论在非饱和土中的应用,基于多孔介质理论在多场耦合条件下土体复杂的行为是当前值得研究的问题。     

非比例加载下两相介质弹塑性特性的有限元分析

在Ｉｌｙｕｓｈｉｎ五维应变空间下,利用弹塑性有限变形的有限元法。研究了两相介质在非比例加载下的弹塑性特性,提出了应力模、延迟角随应变路径、转折角变化的近似计算公式,计算结果还表明,工程材料的许多复杂本构特性都是由第二相介质引起的．     

Axisymmetrical analytical solution for vertical vibration of end-bearing pile in saturated viscoelastic soil layer

Based on elasticity and the theory of saturated porous media, and regarding the pile and the soil as a single phase elastic and a saturated viscoelastic media, respectively, the dynamical behavior of vertical vibration of an end-bearing pile in a saturated viscoelastic soil layer is investigated in the frequency domain using the Helmholtz decomposition and variable separation method. The axisymmetrical analytical solutions for vertical vibrations of the pile in a saturated viscoelastic soil layer are obtained, and the analytical expression of the dynamical complex stiffness of the pile top is presented. Responses of dynamic stiffness factor and equivalent damping of pile top with respect to the frequency are shown in figures using a numerical method. Effects of the saturated soil parameters, modulus ratio of the pile to soil, slenderness ratio of pile and pile＇s Poisson ratio, etc. on the stiffness factor and damping are examined. It is shown that, due to the effect of the transversal deformation of the pile and the radial force of the saturated viscoelastic soil acting on the pile, the dynamic stiffness factor and the damping derived from the axisymmetrical solution are greatly different from those derived from the classical Euler-Bernoulli rod model, especially at some specific excitation frequencies. Therefore, there are limitations on applicability of the Euler-Bernoulli rod model in analyzing verticai vibration of the pile. More accurate analysis should be based on a three-dimensional model.     

On dissipative effects of paper web adhesion strength

This work is concerned with the adhesion strength between a paper web and a metal roll surface, which is a common situation in paper machines world-wide. It is shown that the classic expression relating the work of adhesion to the peeling angle and web tension is, in general, insufficient. An improved model is suggested to take into account the energy dissipation due to elastic–plastic deformation behavior of wet paper materials. To judge the model, an industrially relevant example of wet newsprint and a mild steel surface is studied. It is found that the agreement between theory and experimental observations is excellent. A key result is that elastic–plastic material behavior must always be included for wet paper materials in peeling processes.