Material Behavior of Spherical Elastic-Plastic Granules under Diametrical Compression

Similar documents
Energy absorption during impact of wet particles in the fluidized bed: experiment and DEM simulation Sergiy Antonyuk

Fine adhesive particles A contact model including viscous damping

Shear dynamics simulations of high-disperse cohesive powder

Acoustic Analysis of the Impact of Moist Spherical Particles

C. Wassgren, Purdue University 1. DEM Modeling: Lecture 07 Normal Contact Force Models. Part II

EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY

Energy absorption during compression and impact of dry elastic-plastic spherical granules

URL: < / >

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT

Constitutive model development for granular and porous materials and modelling of particulate processes

Development and validation of calibration methods for discrete element modelling

This is the accepted version of a paper presented at 2014 IEEE Electrical Insulation Conference (EIC).

Computational Analysis for Composites

3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials

Stabilization of continuous fluidized bed spray granulation - a Lyapunov approach

Inter-particle force and stress models for wet and dry particulate flow at the intermediate flow regime

MULTI-SCALE SIMULATION APPROACH FOR THE FLUIDISED BED SPRAY GRANULATION PROCESS

Towards hydrodynamic simulations of wet particle systems

A multi-cell extension to the Barcelona Basic Model

DEM Study of Fluidized Bed Dynamics During Particle Coating in a Spouted Bed Apparatus

SUPER-INSULATED LONG-TERM HOT WATER STORAGE

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

Estimation of agglomerate properties from experiments for microscale simulations

IMPACT PROPERTY AND POST IMPACT VIBRATION BETWEEN TWO IDENTICAL SPHERES

Tangential load deflection behaviour at the contacts of soil particles

Numerical Investigation of the Effect of Recent Load History on the Behaviour of Steel Piles under Horizontal Loading

International Journal of Advancements in Research & Technology, Volume 5, Issue 3, March ISSN

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

Characterization of Thermo-mechanical Performances related to Breeder Pebble Beds by Discrete Element Method

Final Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13

Discrete Element Modelling of a Reinforced Concrete Structure

Module-4. Mechanical Properties of Metals

Introductory guide to measuring the mechanical properties of nanoobjects/particles

Experimental Investigation of Fully Plastic Contact of a Sphere Against a Hard Flat

Measuring the flow properties of powders. FT4 Powder Rheometer. freemantechnology

Effect of Applied Vibration on Silo Hopper Design

Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study

Impact Simulation of Extreme Wind Generated Missiles on Radioactive Waste Storage Facilities

Indentation Energy in Bending of Sandwich Beams with Composite Laminated Faces and Foam Core

Mechanical Properties of Materials

Computational Modelling of the Surface Roughness Effects on the Thermal-elastohydrodynamic Lubrication Problem

A DISCRETE ELEMENT METHODOLOGY FOR THE ANALYSIS OF COHESIVE GRANULAR BULK SOLID MATERIALS. Liz Ivoneth Del Cid

The University of Melbourne Engineering Mechanics

Contact, Adhesion and Rupture of Elastic Solids

The Effect of Coefficient of Restitution for Two Collision Models

Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 pebbles

The plastic behaviour of silicon subjected to micro-indentation

Review of Fluid Mechanics

Indentation tests of aluminium honeycombs

Elastic-plastic Contact of a Deformable Sphere Against a Rigid Flat for Varying Material Properties Under Full Stick Contact Condition

Cubzac-les-Ponts Experimental Embankments on Soft Clay

MHA042 - Material mechanics: Duggafrågor

Discrete Element Study of Aggregate Damage during Asphalt Compaction

Theory at a Glance (for IES, GATE, PSU)

Determination of single-particle drying kinetics in an acoustic levitator

Available online at ScienceDirect. Procedia Engineering 102 (2015 )

Table of Contents. Foreword... xiii Introduction... xv

The Influence of Contact Friction on the Breakage Behavior of Brittle Granular Materials using DEM

On the influence of the grain size distribution curve on P-wave velocity, constrained elastic modulus M max and Poisson s ratio of quartz sands

CONSOLIDATION OF SOIL

Analysis of contact deformation between a coated flat plate and a sphere and its practical application

Highly nonlinear solitary wave velocity measurement with a modified Michelson interferometer

Simulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran

Summary PHY101 ( 2 ) T / Hanadi Al Harbi

SHEAR-ASSISTED FLUIDIZED BED POWDER-COATING

Experimental research on electrical resistance of microcontacts

Modeling the dependence of the coefficient of restitution on the impact velocity in elasto-plastic collisions

Chapter 1 - Soil Mechanics Review Part A

Particle flow simulation of sand under biaxial test

Identification of model parameters from elastic/elasto-plastic spherical indentation

TC211 Workshop CALIBRATION OF RIGID INCLUSION PARAMETERS BASED ON. Jérôme Racinais. September 15, 2015 PRESSUMETER TEST RESULTS

Discrete element modeling of self-healing processes in damaged particulate materials

4.MECHANICAL PROPERTIES OF MATERIALS

Multiscale modeling of failure in ABS materials

CONTACT USING FRICTION HYSTERESIS UDC Adrian Kraft, Roman Pohrt

Characterization of particle impacts using acoustic emission

Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation

A contact model for sticking of adhesive mesoscopic particles

Micro-macro Modeling of Particle Crushing Based on Branch Lengths. Esmaeel Bakhtiary 1, and Chloé Arson 2

MOLECULAR SIEVES UOP MOLECULAR SIEVES*

Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions

OPTICAL DETECTION OF SHORT COLLISION TIMES

Particle Containment And Immobilization In Roll Good Materials

*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials

Chapter 7 Mixing and Granulation

Micromechanics of granular materials: slow flows

Loading σ Stress. Strain

Sticking of adhesive particles in elasto-plastic collisions

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

Geology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)

Numerical simulations for isostatic and die compaction of powder by the discrete element method

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

Micro-meso draping modelling of non-crimp fabrics

Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling

8. Contact Mechanics DE2-EA 2.1: M4DE. Dr Connor Myant 2017/2018

EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS

Chapter 6: Mechanical Properties of Metals. Dr. Feras Fraige

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Size Effects In the Crushing of Honeycomb Structures

Rock Mechanics Laboratory Tests for Petroleum Applications. Rob Marsden Reservoir Geomechanics Advisor Gatwick

Transcription:

Material Behavior of Spherical Elastic-Plastic Granules under Diametrical Compression Alexander Russell 1, Peter Müller 1 and Jürgen Tomas 1 1 Mechanical Process Engineering, Otto von Guericke University, D-39106 Magdeburg, Germany ABSTRACT The material behavior of dry, moist and wet spherical elastic-plastic granules subjected to multiple stressing conditions have been experimentally studied using diametrical compression tests. The force-displacement curves were approximated using the Hertz model for non-linear non-adhesive elastic contact behavior, the Stronge- Antonyuk model for non-linear non-adhesive elastic contact behavior with viscous effects and the Tomas-Antonyuk model for non-linear non-adhesive elastic-plastic contact behavior. The energy characteristics for single particle stressing by diametrical compression were analyzed and a corresponding energetic restitution coefficient was determined. Furthermore, the influences of particle size and moisture content on the material behavior and the energetic restitution coefficient have been discussed. 1. INTRODUCTION In general, particulate solids such as granules and agglomerates are stressed repeatedly by mechanical forces at random locations and intensities during almost every industrial mechanical operation such as processing, transportation, storage and handling. Several quality deteriorative processes advance during prolonged exposure to such forces. The nature, characteristics, intensities, velocities and frequencies of the stressing forces determine the rate and the magnitude of the deformation and the breakage behavior of the stressed material. Extensive theoretical, experimental and simulative efforts are necessary to realize a good understanding of the mechanisms and characteristics of the material behavior of particles under mechanical stressing conditions. Since several decades, quasi-static compression and dynamic impact have been the most popular single particle stressing techniques that have been used worldwide in the experimental approaches. In literature [1 5] one can find a lot of interesting research works striving to achieve an in-depth understanding of the material behavior of granular and particulate materials under different mechanical stressing conditions. In spite of such notable developments, there still exists the inability to characterize the material, comprehensively understand the interactions within the material in the micro level and above all manufacture designer materials with optimum properties for appropriate applications. Thus, the material behavior of granular and particulate materials continues to remain as an area requiring intense investigation. This work aims to accomplish a deeper understanding of the material behavior of spherical elastic-plastic granules under multiple stressing conditions as observed during fluidized bed granulation, drying, pneumatic conveying, mixing, etc., and strives to arrive at conclusive dependencies of granular compression behavior on particle size and moisture content. Further evaluations concerning contact deformation and energy characteristics are also presented in this paper. 2. EXPERIMENTAL 2.1. Material Sample Industrially manufactured synthetic hygroscopic spherical elastic-plastic zeolite 4A granules were chosen as test samples (Figure 1). They are produced by Chemiewerk Bad Köstritz, Germany and are marketed under the commercial name Köstrolith 4A. Their chemical formula is Na 2 O Al 2 O 3 2SiO 2 nh 2 O. They exhibit a uniform pore size of diameter 0.38 nm as well as possess a well-defined pore and channel system [6]. This well-defined pore and channel system give them the remarkable ability to selectively retain molecules of unique sizes and chemical properties within their porous network. For this reason, they are widely known as molecular sieves. They find their prime application in the drying of air, noble gases, aliphatic and aromatic liquids. Furthermore, they are used in the cleaning of technical gases and petrochemicals.

Figure 1 Digital image of the examined granules [7] In order to investigate the influence of particle size on the material behavior under quasi-static compression, three different size ranged granule samples were chosen. Due to the observed high value of sphericity (see Table 1), the granules have been assumed as perfectly homogeneous isotropic spheres and therefore the principles of continuum mechanics and contact mechanics for spheres have been comprehensively applied to arrive at the results described in this study. Similarly, to investigate the influence of moisture content, the granules were moistened using wet air in an air-conditioned climatic chamber at a temperature of 45 C and a relative humidity of 95 % to achieve a certain moist condition (where the pore saturation degree S ranges between zero and unity) and by temporarily immersing the granules in a water bath to achieve a totally wet condition (where the pore saturation degree S equals unity). Mean diameter d 50 in mm Table 1 Physical properties of the examined granules Moisture content -1 X W in kg H20 kg DS Sphericity ψ Porosity ε in % Pore saturation degree S Wet air Water bath Wet air Water bath 1.745 0.980 50.8 0.30 0.45 0.67 1 2.033 0.984 50.1 0.30 0.46 0.65 1 3.267 0.984 53.6 0.30 0.49 0.62 1 2.2. Measurement Technique A strength tester produced by Etewe Testing and Automation Systems, Germany was used to carry out the experiments described in this study. It is widely used by the departments of research and development in several chemical, pharmaceutical, food, ceramic and pigment industries to characterize granular products [7]. Particles of sizes ranging from 0.05 mm to 5 mm can be diametrically compressed with forces up to 2 kn at a wide range of loading velocities. Figure 2 depicts a schematic representation of the principle of quasi-static diametrical compression tests. Force levels with increasing intensities were chosen on the basis of their particle sizes and moisture contents such that the final force level was approximately equal to 90% of the average breakage force. The granules were subjected to diametrical loading and subsequent unloading at a unique contact position with a constant loading velocity of v L = 0.02 mm s -1. In order to achieve representative results, each of the granule samples were tested thirty times for each of the selected compressive force levels. Figure 2 Principle of quasi-static diametrical compression tests [8]

3. RESULTS AND DISCUSSION 3.1. Force-Displacement Behavior On diametrical loading, the granule initially undergoes elastic deformation until the yield point as described by the Hertzian theory. Hertz described the non-adhesive quasi-static collision between two perfectly elastic isotropic solids mutually exerting a certain finite pressure on each other by means of an elliptical distribution of the contact pressure within the contact area [9]. In case of quasi-static diametrical compression of a spherical granule between two rigid flat walls, the normal elastic non-adhesive contact force in analogy with the Hertzian theory is given by F N,el = 0.33 E (1- ν 2 ) -1 d 0.5 s N 1.5 (1) where F N,el is the normal elastic non-adhesive contact force, s N the normal displacement, E the modulus of elasticity of the spherical granule, d its diameter and ν its Poisson s ratio. Figure 3 shows a good fit of the Hertz model for the observed force-displacement behavior within the yield limit. Figure 3 Observed force-displacement behavior and the suitable approximations applied to it On achieving the yield criterion at the yield point, plastic deformation adds on to the existing elastic deformation constituting the elastic-plastic deformation. Tomas defined the dimensionless contact area coefficient κ A for an adhesive contact of fine particle which analytically represents the ratio of the plastic deformation area A pl to the total deformation area A total including a certain constant normal elastic displacement according to [10] κ A = 0.67 + 0.33 A pl A total -1 = 1-0.33 s N,y 0.33 s N -0.33 (2) where s N,y is the normal contact displacement at the yield point (Figure 3). In case of quasi-static diametrical compression of a spherical granule between two rigid flat walls, the normal elastic-plastic non-adhesive contact force in analogy with the Tomas-Antonyuk model is given by [8, 11, 12] F N,el-pl = 0.5 λ el-pl E (1-ν 2 ) -1 d 0.5 s N,y 0.5 (1 0.33 s N,y 0.33 s N -0.33 ) s N (3) where F N,el-pl is the normal elastic-plastic non-adhesive contact force and λ el-pl the elastic-plastic fit parameter. However, Eqn. (3) represents the resultant normal elastic-plastic non-adhesive contact force comprising of a certain linear and a slightly non-linear normal elastic-plastic contact. The respective linear and non-linear normal elasticplastic non-adhesive contact force contributions namely, F N,lin,el-pl and F N,non-lin,el-pl to the overall normal elasticplastic non-adhesive contact force is given by F N,lin,el-pl = 0.5 λ el-pl E (1-ν 2 ) -1 d 0.5 s N,y 0.5 s N (4)

F N,non-lin,el-pl = - 0.17 λ el-pl E (1-ν 2 ) -1 d 0.5 s N,y 0.83 s N 0.67 (5) Following a slight pause on reaching the maximum normal displacement s N,max (Figure 3) at the maximum normal force, diametrical unloading at the same constant velocity of v L = 0.02 mm s -1 was carried out. A sudden deviation from the characteristic elastic recovery of the unloading curve was observed just before the end of its course resulting in a faulty value of s N,U which is the normal displacement at unloading corresponding to zero normal force (Figure 3). This observation is due to the adhesion occurring at the contacts of the granule and the piston plates of the strength tester. This limitation was overcome by approximating the unloading curve based on the Stronge- Antonyuk model according to [12, 13] F N,U,non-lin,el = 0.33 E (1-ν 2 ) -1 µ D 0.5 d 0.5 (2 s N,max s N,y -1-1) 0.25 (s N s N,U ) 1.5 (6) where F N,U,non-lin,el is the normal elastic non-adhesive unloading force and µ D the coefficient that considers additional displacements due to viscous effects. 3.1. Energetic Restitution Coefficient The restitution coefficient e is a characteristic material parameter which is governed by several interdependent dissipative mechanisms [7]. It gives a vivid picture of the material behavior of a body at the collision event of a certain dynamic stressing condition. The restitution coefficient can be mathematically expressed as e = v rebound v impact -1 (7) where v rebound and v impact are the respective velocities of the body at the rebound and impact phases of a dynamic collision event. However, due to the dependency of the restitution coefficient on the energy characteristics occurring during the collision event, an energetic restitution coefficient e E is used in most cases to describe the material behavior during quasi-static collision events such as diametrical compression. The energetic restitution coefficient is given by [13] e E = W non-lin,el,r 0.5 W total,c -0.5 (8) where W non-lin,el,r is the elastic energy released at restitution and W total,c the total energy absorbed at compression which can be determined according to W total,c = W non-lin,el + W lin,el-pl + W non-lin,el-pl (9) W non-lin,el = 0.4 F N,el s N,y (10) W lin,el-pl = 0.5 F N,lin,el-pl (s N,max 2 - s N,y 2 ) s N,max -1 (11) W non-lin,el-pl = 0.6 F N,non-lin,el-pl (s N,max - s N,y 1.67 s N,max -0.67 ) (12) W non-lin,el,r = 0.4 F N,U,non-lin,el (s N,max - s N,U ) (13) where W non-lin,el is the energy absorbed at the non-linear elastic contact, W lin,el-pl the energy absorbed at the linear elastic-plastic contact and W non-lin,el-pl the energy absorbed at the non-linear elastic-plastic contact. No trending dependency of the energetic restitution coefficient is observed on the investigated particle size ranges (Figure 4). However, from the results furnished in this study it is seen that the energetic restitution coefficient remains approximately constant within the investigated range of particle sizes. With increasing moisture content the energetic restitution coefficient decreases significantly as shown in Figure 4. The deformation of pore geometry, swelling up of solid bridges, weakening of bridge bonds and many more such undesirable effects caused by the moisture present within the pores are responsible for this observation.

Figure 4 Energetic restitution coefficient expressed as a function of particle size and moisture content 4. CONCLUSIONS The material behavior of spherical elastic-plastic zeolite 4A granules of different particle sizes and moisture contents under compressive loading had been experimentally studied using single particle diametrical compression tests. Valuable information regarding their processing, transportation, storage and handling have been achieved. The experimental force-displacement behavior was approximated using the Hertz model for non-linear non-adhesive elastic contact behavior, the Stronge-Antonyuk model for non-linear non-adhesive elastic contact behavior with viscous effects and the Tomas-Antonyuk model for non-linear non-adhesive elastic-plastic contact behavior. The rate independent energy characteristics at compression and subsequent restitution was studied. An energetic restitution coefficient characterizing the material behavior under diametrical compression was determined. With increasing particle size, the energetic restitution coefficient approximately remains to be constant within the investigated range of the particle sizes used in this study. With increasing moisture content, the energetic restitution coefficient decreases significantly. This observed behavior is due to the absorption and release of inner stresses during wetting and subsequent drying processes where several localized dislocations, propagation of micro-cracks, deformation of pore network, swelling up of solid bridges, weakening of bridge bonds, etc., occur. As a consequence, the inelastic strain energy contributions to the total strain energy absorption during compression increases resulting in significant irreversible deformations and energy loss. ACKNOWLEDGEMENTS The authors express their sincere gratitude to the Ministry of Science and Economy, Saxony Anhalt, Germany, for the financial support of this study through the Graduate Scholarship of the Land. 5. NOMENCLATURE A pl : plastic deformation area [m 2 ] A total : total deformation area [m 2 ] d : particle diameter [m] d 50 : mean particle diameter [mm] E : modulus of elasticity [N m -2 ] e e E : restitution coefficient : energetic restitution coefficient F N,el : normal non-adhesive elastic contact force [N] : normal non-adhesive elastic-plastic contact force [N] F N,el-pl

F N,lin,el-pl F N,non-lin,el F N,non-lin,el-pl F N,U,non-lin,el S s N s N,max s N,U s N,y : normal loading force at non-adhesive linear elastic-plastic contact [N] : normal loading force at non-adhesive non-linear elastic contact [N] : normal loading force at non-adhesive non-linear elastic-plastic contact [N] : normal unloading force at non-adhesive non-linear elastic contact [N] : pore saturation degree : normal contact displacement [m] : maximum normal contact displacement [m] : normal contact displacement at unloading corresponding to zero unloading force [m] : normal contact displacement at yield point [m] v impact : impact velocity [m s -1 ] v rebound : rebound velocity [m s -1 ] v L : loading velocity [mm s -1 ] W lin,el-pl : energy absorbed due to linear elastic-plastic contact [N m] W non-lin,el : energy absorbed due to non-linear elastic contact [N m] W non-lin,el-pl W non-lin,el,r : energy absorbed due to non-linear elastic-plastic contact [N m] : elastic energy restored at restitution [N m] W total,c X W : total energy absorbed at compression [N m] : moisture content [kg H2O kg -1 DS ] ε : porosity [%] κ A λ el-pl μ D ν ψ : dimensionless contact area coefficient : elastic-plastic fit parameter : coefficient for viscous effects : Poisson s ratio : sphericity 6. REFERENCES [1] Antonyuk, S., Tomas, J., Heinrich, S. and Mörl, L., Breakage behaviour of spherical granulates by compression, Chemical Engineering Science, Vol. 60, No. 14, 2005, pp. 4031-4044. [2] Aman, S., Tomas, J. and Kalman, H., Breakage probability of irregularly shaped particles, Chemical Engineering Science, Vol. 65, No. 5, 2010, pp. 1503-1512. [3] Rozenblat, Y., Grant, E., Levy, A., Kalman, H. and Tomas, J., Selection and breakage functions of particles under impact loads, Chemical Engineering Science, Vol. 71, 2012, pp. 56-66. [4] Müller, P., Seeger, M. and Tomas, J., Compression and breakage behavior of γ-al 2 O 3 granules, Powder Technology, Vol. 237, 2013, pp. 125-133. [5] Khanal, M., Schubert, W. and Tomas, J., DEM simulation of diametrical compression test on particle compounds, Granular Matter, Vol. 7, No. 2-3, 2005, pp. 83-90. [6] Schuman, K., Unger, B. and Brandt, A., Zeolithe als Sorptionsmittel (Zeolite as sorbent), Chemie Ingenieur Technik, Vol. 82, No. 6, 2010, pp. 929-940. [7] Russell, A., Müller, P. and Tomas, J., Multiple compression of moist spherical elastic-plastic zeolite 4A granules, Chemical Engineering and Technology, Vol. 36, No. 7, 2013, pp. 1240-1248. [8] Müller, P., Druck-, Stoß- und Bruchverhalten feuchter kugelförmiger Granulate (Compression, impact and breakage behavior of moist spherical granules), Dissertation (Ph.D. Thesis), Otto von Guericke University of Magdeburg, Docupoint Verlag, 2011, pp. 1-235. [9] Hertz, H., Über die Berührung fester elastischer Körper (About the contact of solid elastic bodies), Journal für die reine und angewandte Mathematik, Vol. 92, 1881, pp. 156-171. [10] Tomas, J., Assessment of mechanical properties of cohesive solids part 1: particle contact constitutive model, Particulate Science and Technology, Vol. 19, No. 2, 2001, pp. 95-110. [11] Tomas, J., Adhesion of ultrafine particles Energy absorption at contact, Chemical Engineering Science, Vol. 62, No. 21, 2007, pp. 5925-5939. [12] Antonyuk, S., Deformations- und Bruchverhalten von kugelförmigen Granulaten bei Druck- und Stoßbeanspruchung (Deformation and breakage behavior of spherical granules by compression and impact stressing), Dissertation (Ph.D. Thesis), Otto von Guericke University of Magdeburg, Docupoint Verlag, 2006, pp. 1-276. [13] Stronge, W.J., Impact Mechanics, Cambridge University Press, United Kingdom, 2000, pp. 1-280.