Analysis of the flowability of cohesive powders using Distinct Element Method

Computer simulations using Distinct Element Method (DEM) have been carried out to investigate the effect of cohesion on the flowability of polydisperse particulate systems. For this purpose, two assemblies with different values of surface energy and made of 3000 spheres with the mechanical properties of glass beads were considered. The analysis of the flowability of the powders is presented in terms of the unconfined yield stress as a function of strain rate for different pre-consolidation loads. For values of the surface energy of 1.0 J/m2 and strain rates lower than 6 s− 1, the unconfined yield stress does not change significantly indicating a quasi-static behaviour of the particulate assemblies during the compression process. For larger strain rates, the unconfined yield stress varies with the power index of 1.2 of the strain rate. The influence of the pre-consolidating stress on the powder behaviour has also been investigated and a flow factor was obtained from the linear relationship between the unconfined yield stress and pre-consolidation stress. The computer simulations show qualitatively a good agreement with the experimental trends on highly cohesive powder flow behaviour

Numerical evaluation of three non-coaxial kinematic models using the distinct element method for elliptical granular materials

This is the accepted version of the following article: [Jiang, M. J., Liu, J. D., and Arroyo, M. (2016) Numerical evaluation of three non-coaxial kinematic models using the distinct element method for elliptical granular materials. Int. J. Numer. Anal. Meth. Geomech., 40: 2468–2488. doi: 10.1002/nag.2540.], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nag.2540/fullThis paper presents a numerical evaluation of three non-coaxial kinematic models by performing Distinct Element Method (DEM) simple shear tests on specimens composed of elliptical particles with different aspect ratios of 1.4 and 1.7. The models evaluated are the double-shearing model, the double-sliding free-rotating model and the double slip and rotation rate model (DSR2 model). Two modes of monotonic and cyclic simple shear tests were simulated to evaluate the role played by the inherent anisotropy of the specimens. The main findings are supported by all the DEM simple shear tests, irrespective of particle shape, specimen density or shear mode. The evaluation demonstrates that the assumption in the double-shearing model is inconsistent with the DEM results and that the energy dissipation requirements in the double-sliding free-rotating model appear to be too restrictive to describe the kinematic flow of elliptical particle systems. In contrast, the predictions made by the DSR2 model agree reasonably well with the DEM data, which demonstrates that the DSR2 model can effectively predict the non-coaxial kinematic behavior of elliptical particle systems.Peer ReviewedPostprint (author's final draft

Three-dimensional inelastic analysis for hot section components, BEST 3D code

The goal is the development of an alternative stress analysis tool, distinct from the finite element method, applicable to the engineering analysis of gas turbine engine structures. The boundary element method was selected for this development effort on the basis of its already demonstrated applicability to a variety of geometries and problem types characteristic of gas turbine engine components. Major features of the BEST3D computer program are described, and some of the significant developments carried out as part of the Inelastic Methods Contract are outlined

A new procedure for calculating contact stresses in gear teeth

A numerical procedure for evaluating and monitoring contact stresses in meshing gear teeth is discussed. The procedure is intended to extend the range of applicability and to improve the accuracy of gear contact stress analysis. The procedure is based upon fundamental solution from the theory of elasticity. It is an iterative numerical procedure. The method is believed to have distinct advantages over the classical Hertz method, the finite-element method, and over existing approaches with the boundary element method. Unlike many classical contact stress analyses, friction effects and sliding are included. Slipping and sticking in the contact region are studied. Several examples are discussed. The results are in agreement with classical results. Applications are presented for spur gears

Simulation of the effect of bond strength on the breakage pattern of agglomerates by Distinct Element Method

Chemical, pharmaceutical and food industries amongst many others use agglomerates either as intermediate or manufactured products. The mechanical strength of agglomerates under impact or shear deformation during handling and processing is of great interest to these industries for optimising product specification and functionality. The effect of surface energy on agglomerate behaviour under impact has been investigated using Distinct Element Method (DEM). Four different agglomerates were formed and impacted against a target along the direction of gravity for three different values of the surface energy (0.35, 3.5 and 35.0 J/m²). The agglomerate breakage pattern was influenced by the surface energy and a transition in the mode of failure of agglomerates was observed when the surface energy was varied. Based on the previous work, the surface energy is expressed in terms of Weber Number, We=(V-V0²)ρD/γ. Agglomerates showed extensive deformation under impact at the lowest value of surface energy (0.35 J/m²) and no evidence of fragmentation was found for any value of impact velocity. In this case the agglomerates behaved macroscopically in a ductile mode. At values of surface energy larger than 3.5 J/m² the agglomerates fragmented at the same time as local damage around the impact site occurred. This type of behaviour is typical of semi-brittle material failure. Therefore, the breakage pattern of agglomerates is influenced by the surface energy

Influence of measurement volume on predicted attrition by the distinct element method

During agitated drying and mixing processes, particle beds are exposed to shear deformation. This leads to particle attrition, the extent of which is dependent on the prevailing stresses and strains in the bed. The distributions of shear stresses and strain rates within the bed are highly non-uniform, requiring attention to localised conditions. Therefore a narrow angular sector of the bed is divided radially and vertically into a number of measurement cells, within which the stresses and strain rates are calculated throughout one rotation by the Distinct Element Method. These are then used in an empirical relationship of material breakage to predict the extent of attrition due to agitation. Here we investigate the influence of the measurement cell size on the estimated stresses and strain rates, and the subsequent effect on the predicted attrition. The measurement cell size is altered by varying the measurement sector size and the number of radial and vertical divisions within it. The median particle size is also varied to establish its influence on the predicted attrition. An increase in the average number of particles in a given cell, by varying the particle size or measurement cell dimensions, leads to a reduction in the estimated stresses and strain rates, and therefore a reduction in the predicted attrition. Comparison of the predicted attrition with the experimental breakage in the agitated vessel shows that the prediction method is accurate when the cell dimensions are comparable to the width of a naturally occurring shear band

ACCESS 3. Approximation concepts code for efficient structural synthesis: User's guide

A user's guide is presented for ACCESS-3, a research oriented program which combines dual methods and a collection of approximation concepts to achieve excellent efficiency in structural synthesis. The finite element method is used for structural analysis and dual algorithms of mathematical programming are applied in the design optimization procedure. This program retains all of the ACCESS-2 capabilities and the data preparation formats are fully compatible. Four distinct optimizer options were added: interior point penalty function method (NEWSUMT); second order primal projection method (PRIMAL2); second order Newton-type dual method (DUAL2); and first order gradient projection-type dual method (DUAL1). A pure discrete and mixed continuous-discrete design variable capability, and zero order approximation of the stress constraints are also included