Microstructure evolution of compressible granular systems under large deformations

We report three-dimensional particle mechanics static calculations that predict the microstructure evolution during die-compaction of elastic spherical particles up to relative densities close to one. We employ a nonlocal contact formulation that remains predictive at high levels of confinement by removing the classical assumption that contacts between particles are formulated locally as independent pair-interactions. The approach demonstrates that the coordination number depends on the level of compressibility, i.e., on Poisson's ratio, of the particles. Results also reveal that distributions of contact forces between particles and between particles and walls, although similar at jamming onset, are very different at full compaction. Particle–wall forces are in remarkable agreement with experimental measurements reported in the literature, providing a unifying framework for bridging experimental boundary observations with bulk behavior.     

Implicit frictional-contact model for soft particle systems

We introduce a novel numerical approach for the simulation of soft particles interacting via frictional contacts. This approach is based on an implicit formulation of the Material Point Method, allowing for large particle deformations, combined with the Contact Dynamics method for the treatment of unilateral frictional contacts between particles. This approach is both precise due to the treatment of contacts with no regularization and artificial damping parameters, and robust due to implicit time integration of both bulk degrees of freedom and relative contact velocities at the nodes representing the contact points. By construction, our algorithm is capable of handling arbitrary particle shapes and deformations. We illustrate this approach by two simple 2D examples: a Hertz contact and a rolling particle on an inclined plane. We also investigate the compaction of a packing of circular particles up to a solid fraction well above the jamming limit of hard particles. We find that, for the same level of deformation, the solid fraction in a packing of frictional particles is above that of a packing of frictionless particles as a result of larger particle shape change.     

粗糙表面接触研究进展

工程中的任何接触表面都不会是绝对光滑的,对于某些特殊接触问题必须考虑表面形貌的接触效应.粗糙表面接触研究是一门年轻的交叉学科.它主要研究粗糙表面接触力学,接触导电和导热,接触密封,接触振动和噪声,接触刚度,骨关节接触等问题,涉及固体力学,流体力学,随机过程及概率统计,计量学,电学和传热学,振动理论以及生物学等学科知识.本文综述了粗糙表面接触的研究进展.      

A microstructural elastoplastic model for unsaturated granular materials

The homogenization technique is used to obtain an elastoplastic stress–strain relationship for dry, saturated and unsaturated granular materials. Deformation of a representative volume of material is generated by mobilizing particle contacts in all orientations. In this way, the stress–strain relationship can be derived as an average of the mobilization behavior of these local contact planes. The local behavior is assumed to follow a Hertz–Mindlin’s elastic law and a Mohr–Coulomb’s plastic law. For the non-saturated state, capillary forces at the grain contacts are added to the contact forces created by an external load. They are calculated as a function of the degree of saturation, depending on the grain size distribution and on the void ratio of the granular assembly. Numerical simulations show that the model is capable of reproducing the major trends of a partially saturated granular assembly under various stress and water content conditions. The model predictions are compared to experimental results on saturated and unsaturated samples of silty sands under undrained triaxial loading condition. This comparison shows that the model is able to account for the influence of capillary forces on the stress–strain response of the granular materials and therefore, to reproduce the overall mechanical behavior of unsaturated granular materials.     

Study of the effect of contact force model on the dynamic response of mechanical systems with dry clearance joints: computational and experimental approaches

The main objective of this work is to present a computational and experimental study on the contact forces developed in revolute clearance joints. For this purpose, a well-known slider-crank mechanism with a revolute clearance joint between the connecting rod and slider is utilized. The intra-joint contact forces that are generated at these clearance joints are computed by considering several different elastic and dissipative approaches, namely those based on the Hertz contact theory and the ESDU tribology-based cylindrical contacts, along with a hysteresis-type dissipative damping. The normal contact force is augmented with the dry Coulomb’s friction force. In addition, an experimental apparatus is used to obtained some experimental data in order to verify and validate the computational models. From the outcomes reported in this paper, it is concluded that the selection of the appropriate contact force model with proper dissipative damping plays a significant role in the dynamic response of mechanical systems involving contact events at low or moderate impact velocities.     

Continuous contact force models for impact analysis in multibody systems

One method for predicting the impact response of a multibody system is based on the assumption that the impacting bodies undergo local deformations and the contact forces are continuous. In a continuous analysis, the integration of the system equations of motion is carried out during the period of contact; therefore, a model for evaluating the contact forces is required. In this paper, two such contact force models are presented, both Hertzian in nature and based upon the direct-central impact of two solid particles.At low impact velocities, the energy dissipation during impact can be represented by material damping. A model is constructed based on the general trend of the Hertz contact law in conjuction with a hysteresis damping function. The unknown parameters are determined in terms of a given coefficient of restitution and the impact velocity. When local plasticity effects are the dominant factor accounting for the dissipation of energy at high impact velocities, a Hertzian contact force model with permanent indentation is constructed. Utilizing energy and momentum considerations, the unknown parameters in the model are again evaluated. The two particle models are generalized to an impact analysis between two bodies of a multibody system.     

Nonlocal implicit gradient-enhanced elasto-plasticity for the modelling of softening behaviour

An improved gradient-enhanced approach for softening elasto-plasticity is proposed, which in essence is fully nonlocal, i.e. an equivalent integral nonlocal format exists. The method utilises a nonlocal field variable in its constitutive framework, but in contrast to the integral models computes this nonlocal field with a gradient formulation. This formulation is considered ‘implicit’ in the sense that it strictly incorporates the higher-order gradients of the local field variable indirectly, unlike the common (explicit) gradient approaches. Furthermore, this implicit gradient formulation constitutes an additional partial differential equation (PDE) of the Helmholtz type, which is solved in a coupled fashion with the standard equilibrium condition. Such an approach is particularly advantageous since it combines the long-range interactions of an integral (nonlocal) model with the computational efficiency of a gradient formulation. Although these implicit gradient approaches have been successfully applied within damage mechanics, e.g. for quasi-brittle materials, the first attempts were deficient for plasticity. On the basis of a thorough comparison of the gradient-enhancements for plasticity and damage this paper rephrases the problem, which leads to a formulation that overcomes most reported problems. The two-dimensional finite element implementation for geometrically linear plain strain problems is presented. One- and two-dimensional numerical examples demonstrate the ability of this method to numerically model irreversible deformations, accompanied by the intense localisation of deformation and softening up to complete failure.     

Tangential interaction of elastic spherical particles in contact

The elastic interaction of spherical particles is studied. The distribution of the stress, normal to the contact plane, is determined by the rod model suggested recently, which is applicable in the more wide range of deformations as compared with the classical Hertz law. In the rod model context an inner part of compressed particles is regarded as an elastic cylindrical rod, which radius is equal to the contact radius. The rod reaction is added to the normal particle interaction corresponding with the Hertz solution. The resulting normal force passes into the Hertz solution for infinitesimal deformations and gives stronger particle repulsion for finite deformations. Here we solve the Mindlin problem for the rod model, i.e., derive the tangential interaction of initially compressed particles when a relative displacement takes place. The analytical expressions, which determine the total displacement of the sphere’s centers and the corresponding tangential force, are derived. So, the generalization of the classical Mindlin law is obtained for the rod model.     

微结构非局部效应的多重性模型及在微梁弯曲中的应用

基于Eringen非局部弹性理论,直接利用逐次逼近法推导了非局部应力场的精确表达,该精确的非局部应力可具体表示为一个无穷级数的形式. 然后以微梁的横向弯曲和纯弯曲变形为例,建立平衡方程并求解及分析了挠度受非局部效应的影响. 结果表明：根据所取非局部小尺度参数大小的不同,非局部微梁的弯曲挠度可低于也可以高于经典力学下的挠度,非局部效应的增大可提高亦可降低结构的抗弯刚度. 本文结果证明了Wang以及Lim等人分别提出的两种相反的非局部模型的各自正确性. 同时首次发现,弯曲挠度随着非局部效应的增大而上下波动且存在若干跳跃点,挠度是非局部小尺度参数的非单调函数,研究同时给出了一种确定材料非局部常数的建议途径.     

复合材料层合板静压痕接触的等效弹性模量

在复合材料层合板静压痕接触问题中,由于层合板接触区域的纤维走向发生了变化,导致沿厚度方向的弹性模量产生变化并随着凹坑深度的增加而增大。本文建立压头压入后的力学分析模型,得到了与凹坑深度相关的沿厚度方向的等效弹性模量;将等效弹性模量代入目前被广泛接受的修正Hertz接触理论,建立了新的接触力与凹坑深度的关系式。通过与静压痕试验中接触力和凹坑深度的关系进行比较,结果表明该接触关系式能较好地描述接触力和凹坑深度的关系。与文献修正理论相比,在压痕较小时,本文的修正结果较好。     

Regularity for Shearable Nonlinearly Elastic Rods in Obstacle Problems

Based on the Cosserat theory describing planar deformations of shearable nonlinearly elastic rods we study the regularity of equilibrium states for problems where the deformations are restricted by rigid obstacles. We start with the discussion of general conditions modeling frictionless contact. In particular we motivate a contact condition that, roughly speaking, requires the contact forces to be directed normally, in a generalized sense, both to the obstacle and to the deformed shape of the rod. We show that there is a jump in the strains in the case of a concentrated contact force, i.e., the deformed shape of the rod has a corner. Then we assume some smoothness for the boundary of the obstacle and derive corresponding regularity for the contact forces. Finally we compare the results with the case of unshearable rods and obtain interesting qualitative differences. (Accepted January 21, 1998)     

An elasto-plastic model for granular materials with microstructural consideration

In this paper, we have extended the granular mechanics approach to derive an elasto-plastic stress–strain relationship. The deformation of a representative volume of the material is generated by mobilizing particle contacts in all orientations. Thus, the stress–strain relationship can be derived as an average of the mobilization behavior of these local contact planes. The local behavior is assumed to follow a Hertz–Mindlin’s elastic law and a Mohr–Coulomb’s plastic law. Essential features such as continuous displacement field, inter-particle stiffness, and fabric tensor are discussed. The predictions of the derived stress–strain model are compared to experimental results for sand under both drained and undrained triaxial loading conditions. The comparisons demonstrate the ability of this model to reproduce accurately the overall mechanical behavior of granular media and to account for the influence of key parameters such as void ratio and mean stress. A part of this paper is devoted to the study of anisotropic specimens loaded in different directions, which shows the model capability of considering the influence of inherent anisotropy on the stress–strain response under a drained triaxial loading condition.     

Dynamics of flexible mechanical systems with contact-impact and plastic deformations

A computer based formulation for the analysis of mechanical systems is investigated as a feasible method to predict the impact response of complex structural systems. A general methodology for the dynamic analysis of rigid-flexible multibody systems using a number of redundant Cartesian coordinates and the method of the Lagrange multipliers is presented. The component mode synthesis is then used to reduce the number of flexible degrees of freedom. In many impact situations, the individual structural members are overloaded giving rise to plastic deformations in highly localized regions, called plastic hinges. This concept is used by associating revolute nonlinear actuators with constitutive relations corresponding to the collapse behavior of the structural components. The contact of the system components is described using a continuous force model based on the Hertz contact law with hysteresis damping. The effect and importance of structural damping schemes in flexible bodies are also addressed here. Finally, the validity of this methodology is assessed by comparing the results of the proposed models with those obtained in different experimental tests where: a beam collides transversally with a rigid block; a torque box impacts a rigid barrier.     

Fault diagnosis of a rotor bearing system using response surface method

Dynamic analysis of a high-speed rotor bearing systems is challenged by their highly nonlinear and complex properties. Hence, an approximate response surface method (RSM) is utilized to analyze the effects of design and operating parameters on the vibration signature of a rotor-bearing system. This paper focuses on accurate performance prediction, which is essential to the design of high performance rotor bearing system. It considers distributed defects such as internal radial clearance and surface waviness of the bearing components. In the mathematical formulation the contacts between the rolling elements and the races are considered as nonlinear springs, whose stiffnesses are obtained by using Hertzian elastic contact deformation theory. The governing differential equations of motion are obtained by using Lagrange's equations. In terms of the feature that the nonlinear bearing forces act on the system, a reduction method and corresponding integration technique is used to increase the numerical stability and decrease computer time for system analysis. Parameters effects are analyzed together and its influence considered with DOE and Surface Response Methodology are used to predict dynamic response of a rotor-bearing system.     

Modeling intermittent contact for flexible multibody systems

This paper consists of two parts. The first part presents a complementarity based recursive scheme to model intermittent contact for flexible multibody systems. A recursive divide-and-conquer framework is used to explicitly impose the bilateral constraints in the entire system. The presented approach is an extension of the hybrid scheme for rigid multibody systems to allow for small deformations in form of local mode shapes. The normal contact and frictional complementarity conditions are formulated at position and velocity level, respectively, for each body in the system. The recursive scheme preserves the essential characteristics of the contact model and formulates a minimal size linear complementarity problem at logarithmic cost for parallel implementation. For a certain class of contact problems in flexible multibody systems, the complementarity based time-stepping scheme requires prohibitively small time-steps to retain accuracy. Modeling intermittent contact for this class of contact problems motivated the development of an iterative scheme. The second part of the paper describes this iterative scheme to model unilateral constraints for a multibody system with relatively fewer contacts. The iterative scheme does not require a traditional complementarity formulation and allows the use of any higher order integration methods. A comparison is then made between the traditional complementarity formulation and the presented iterative scheme via numerical examples.     

Micromechanical study of elastic moduli of loose granular materials

In micromechanics of the elastic behaviour of granular materials, the macro-scale continuum elastic moduli are expressed in terms of micro-scale parameters, such as coordination number (the average number of contacts per particle) and interparticle contact stiffnesses in normal and tangential directions. It is well-known that mean-field theory gives inaccurate micromechanical predictions of the elastic moduli, especially for loose systems with low coordination number. Improved predictions of the moduli are obtained here for loose two-dimensional, isotropic assemblies. This is achieved by determining approximate displacement and rotation fields from the force and moment equilibrium conditions for small sub-assemblies of various sizes. It is assumed that the outer particles of these sub-assemblies move according to the mean field. From the particle displacement and rotation fields thus obtained, approximate elastic moduli are determined. The resulting predictions are compared with the true moduli, as determined from the discrete element method simulations for low coordination numbers and for various values of the tangential stiffness (at fixed value of the normal stiffness). Using this approach, accurate predictions of the moduli are obtained, especially when larger sub-assemblies are considered. As a step towards an analytical formulation of the present approach, it is investigated whether it is possible to replace the local contact stiffness matrices by a suitable average stiffness matrix. It is found that this generally leads to a deterioration of the accuracy of the predictions. Many micromechanical studies predict that the macroscopic bulk modulus is hardly influenced by the value of the tangential stiffness. It is shown here from the discrete element method simulations of hydrostatic compression that for loose systems, the bulk modulus strongly depends on the stiffness ratio for small stiffness ratios.     

基于有限元法的轮轨蠕滑特性研究

轮轨滚动接触蠕滑率／力理论是轮轨相互作用一系列问题研究的基础．现有的几种计算蠕滑力的理论模型均建立在Hertz接触条件和半空间假设的基础上，已经无法完成对复杂的接触问题的进一步研究．使用有限元参数二次规划法来求解轮轨的三维弹性／弹塑性接触问题，得出了在不同的轴重、牵引力矩、摩擦系数、踏面形状、横向力条件下的轮轨接触力．提出用轮轨接触的轮周位移计算蠕滑率的新方法，并对在各种参数下所得到的蠕滑率进行了分析比较．     

Frictional contact of two rotating elastic disks

We study the problem of constrained uniform rotation of two precompressed elastic disks made of different materials with friction forces in the contact region taken into account. The exact solution of the problem is obtained by the Wiener-Hopf method.An important stage in the study of rolling of elastic bodies is the Hertz theory [1] of contact interaction of elastic bodies with smoothly varying curvature in the contact region under normal compression. Friction in the contact region is assumed to be negligible. If there are tangential forces and the friction in the contact region is taken into account, then the picture of contact interaction of elastic bodies changes significantly. Although the normal contact stress distribution strictly follows the Hertz theory for bodies with identical elastic properties and apparently slightly differs from the Hertz diagram for bodies made of different materials, the presence of tangential stresses results in the splitting of the contact region into the adhesion region and the slip region. This phenomenon was first established by Reynolds [2], who experimentally discovered slip regions near points of material entry in and exit from the contact region under constrained rolling of an aluminum cylinder on a rubber base. The theoretical justification of the partial slip phenomenon in the contact region, discovered by Reynolds [2], can be found in Carter [3] and Fromm [4]. Moreover, Fromm presents a complete solution of the problem of constrained uniform rotation of two identical disks. Apparently, Fromm was the first to consider the so-called “clamped” strain and postulated that slip is absent at the point at which the disk materials enter the contact region.Ishlinskii [5, 6] gave an engineering solution of the problem on slip in the contact region under rolling friction. Considering the problem on a rigid disk rolling on an elastic half-plane, we model this problem by an infinite set of elastic vertical rods using Winkler-Zimmermann type hypotheses. Numerous papers of other authors are surveyed in Johnson’s monograph [7].The exact solution of the problem on the constrained uniform rotation of precompressed rigid and elastic disks under the assumptions of Fromm’s theory is contained in the papers [8, 9]. In the present paper, we generalize the solution obtained in [8, 9] to the case of two elastic disks made of different materials.