Next Article in Journal
Decontamination of Titanium Surface Using Different Methods: An In Vitro Study
Next Article in Special Issue
Neutron Shielding Performance of 3D-Printed Boron Carbide PEEK Composites
Previous Article in Journal
Role of Natural Stone Wastes and Minerals in the Alkali Activation Process: A Review
Previous Article in Special Issue
Preliminary Characterization of Glass/Alumina Composite Using Laser Powder Bed Fusion (L-PBF) Additive Manufacturing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete Element Simulation of the Effect of Roller-Spreading Parameters on Powder-Bed Density in Additive Manufacturing

1
Institute of Manufacturing Engineering, Huaqiao University, Xiamen 361021, China
2
National & Local Joint Engineering Research Center for Intelligent Manufacturing Technology of Brittle Material Products, Huaqiao University, Xiamen 361021, China
3
School of Mechanical Engineering, Xiangtan University, Xiangtan 411105, China
*
Author to whom correspondence should be addressed.
Submission received: 21 April 2020 / Revised: 7 May 2020 / Accepted: 11 May 2020 / Published: 15 May 2020
(This article belongs to the Special Issue Additive Manufacturing Materials and Their Applications)

Abstract

:
The powder-bed with uniform and high density that determined by the spreading process parameters is the key factor for fabricating high performance parts in Additive Manufacturing (AM) process. In this work, Discrete Element Method (DEM) was deployed in order to simulate Al2O3 ceramic powder roller-spreading. The effects of roller-spreading parameters include translational velocity Vs, roller’s rotational speed ω, roller’s diameter D, and powder layer thickness H on powder-bed density were analyzed. The results show that the increased translational velocity of roller leads to poor powder-bed density. However, the larger roller’s diameter will improve powder-bed density. Moreover, the roller’s rotational speed has little effect on powder-bed density. Layer thickness is the most significant influencing factor on powder-bed density. When layer thickness is 50 μm, most of particles are pushed out of the build platform forming a lot of voids. However, when the layer thickness is greater than 150 μm, the powder-bed becomes more uniform and denser. This work can provide a reliable basis for roller-spreading parameters optimization.

Graphical Abstract

1. Introduction

The powder-bed based additive manufacturing (AM), such as Inkjet Three-Dimensional (3D) Printing, selective laser sintering (SLS), and selective laser melting (SLM) are the most widely used additive manufacturing technologies [1,2]. The powder-bed density directly affects the effective powder thermal conductivity; consequently, the temperature distribution in melt pool and performance of AM parts [3]. The powder-bed with uniform and high density is the key factor to fabricate high-performance parts in AM processes [3,4]. Otherwise, the loose packing of powder-bed leads to defects, such as inner holes and poor performance of AM parts.
Currently, the measurement of powder-bed density is still facing challenges. Recently, the method of measuring powder-bed density by weighing the powder inside containers was proposed [5,6]. However, the method is rather uncertain due to heavily depending on the powder spreading process. It is worth mentioning that Usman et al. [4] presented a new method for the measurement of powder-bed density. The powder samples were bonded through an ultraviolet curable polymer at different locations. The compactness density of the powder-bed will be obtained by a nano-computing tomography (CT) scanning results. Although the method is high accuracy, the cost of the experiment is too expensive to be applied in industry.
With AM technology growing rapidly in industry and academia, there has been increasing interest in the application of Discrete Element Method (DEM) in AM processes during the past five years. Zohdi [7,8] and John C. [9] proposed DEM to simulate the process of AM, which establishes a more realistic particle contact model and consider the effects of multiple physical factors on the particle system. Subsequently, much work applied DEM to simulate powder spreading process [10,11,12,13,14,15,16,17]. These investigations indicate that DEM is a helpful tool for simulating powder mechanical behavior to overcome the experimental challenges in AM processes.
There are mainly two main methods for powder spreading: blade-spreading [10,16] and roller-spreading [13,18]. In the blade-spreading process, the straight movement of the blade is simple and easy to control; however, the powder bed cannot be compacted. For the roller-spreading process, the roller movement, including translation and counter-clockwise rotation, is relatively more complex than the blade. The advantage of roller-spreading is that it can get higher powder-bed density than that of the blade-spreading [16]. However, previous works focus on metal powder spreading by blade [6,12].
Eric J.R. [13] studied the effect of spreading conditions on the powder-bed roughness in roller-spreading process and showed that at the increasing the roller’ translational speed leads to an increase in the surface roughness of the powder bed. Nan [19] studied the effects of gap height and roller’ rotational speed on the evolving particle trajectory and spread layer uniformity by DEM. The results show that the segregation extent decreases as the gap height increases or the roller’ rotational speed decreases. However, the roller-spreading parameters of the above studies are simple and roller-spreading parameters are relatively complex, which has an important influence on powder flow and powder-bed densification mechanism. Thus, further research is needed in order to elucidate on the roller-spreading parameters for improving powder-bed density and processing efficiency.
In this work, DEM is used to simulate Al2O3 ceramic powder spreading by roller in AM processes. The Hertz–Mindlin with JKR model of Al2O3 ceramic powder is established. Simultaneously, the DEM parameters are calibrated. Subsequently, powders in the powder spreading process are divided into three zones: I. Powder free-flow zone, II. Powder compression zone, and III. Powder layer zone. A mathematical-statistical method is employed to quantitatively measure the powder-bed compactness density and uniformity of powder-bed density. Finally, the effect of roller- spreading parameters on Al2O3 ceramic powder-bed density is analyzed.

2. Materials and Methods

2.1. Discrete Element Method

The DEM method with “soft ball model” was developed by Cundal and Strack [20]. It has been widely used to simulate powder behavior in AM processes [6,12]. In DEM equations of the linear motion and rotation motion of particle are governed by Newton’s second law:
m i d v i d t = j ( F i j n + F i j t ) + F i g
I i d ω i d t = j ( T i j s + T i j r )
where m i I i v i , and ω i are the mass, moment of inertia, translational velocity, and rotational velocity of particle i; F i j n and F i j t are the normal contact force and tangential contact force imposed on particle i by particle j or wall; F i g is the gravitational force of particle i; and, T i j s and T i j r are the torques on particle i caused by tangential force and rolling friction.
The Hertz–Mindlin JKR model [21] is taken into account simulations to account for adhesion force, which adds an additional normal cohesive force within the contact zone:
F i j , c n = 4 E 3 R α 3 4 π γ E α 3 2
The normal damping force is described as:
F i j , d n = 2 5 6 β S n m * v n r e l ¯
where E i , γ i , R i and E j , γ j , R j are the Young’s Modulus, Poisson ratio, and radius of contact element i and j, v n r e l ¯ is the normal relative velocity between contact element i and j, α is radius contact area. Table 1 defines the above parameters used in the contact model.
The tangential force between element i and j consists of two components: the tangential contact force and tangential damping force are defined, as following [22,23]:
F i j , c t = S t δ t
F i j , d t = 2 5 6 β S t m v t r e l ¯
The tangential force is limited by Coulomb friction according to F i j , c t < μ s F i j n .
The torque on particle i due to the tangential force is:
T i j , c s = R i F i j t
The torque acted on particle i due to the rolling friction is:
T i j r = μ r F i j , c n R i ω i
where μs and μr are the coefficient of static and rolling friction, respectively.
The model calculates attractive cohesion forces, even when the particles are not in real contact. The critical gap between particles δn,c and critical contact spot radius αn,c are given by [24]:
δ n , c = α c 2 R 4 π γ α c E
α n , c = 9 π γ R 2 2 E ( 3 4 2 2 )
When the particles are not in actual contact and the gap between particles is smaller than δn,c, the adhesion reaches a maximum value, called “pull-out force”, and the maximum cohesive force Fpullout is given by [24,25]:
F p u l l o u t = 3 2 π γ R

2.2. Calibration of DEM Parameters

For the DEM simulation of the spreading process, DEM parameters, including powder material property parameters and contact model parameters, play a crucial role in getting the correct results [26,27].
The Al2O3 ceramic powders are widely used in AM processes [28,29,30]. In this work, Al2O3 ceramic powders (Dongchao Technologies, Dongguan, China) are the research object. Particle size, density, shear modulus, and Poisson’s ratio are the powder material property parameters. The particles size and density are measured in order to calibrate the DEM parameters. A dynamic optical particle analyzer measures the particle size (Retsch Camsizer X2, Retsch technology GmbH, Haan, Germany). The particle size D (10), D (50), and D (90) are 24, 48, and 72 μm. Particles density is measured by a Gas Pycnometer (Bei Shi De Instrument-S & T, Beijing, China) with a value of 3820 kg/m3. The micromorphology of powders is observed by secondary electron imaging in the scanning electron microscope (SEM, Phenom ProX, Phenom-World, Eindhoven, Holland). Figure 1 shows the SEM micrograph for Al2O3 ceramic powders around ball shape. The shear modulus and Poisson’s ratio are drawn from reference literature [31]. According to the literature [32], Young’s modulus has little effect on the simulation results. The value of Young’s modulus is set as 3 GPa, which is decreased two orders of magnitude than the actual value in order to reduce computer time consumption.
The contact model parameters include the coefficient of restitution, static friction, and rolling friction. The static friction coefficients of particle-particle and particle-wall are measured by the Jenike shear test defined in ASTM D6773 [33], which has been widely used as a measure for powder flowability. Figure 2 shows the shear test results. The static friction coefficients of powder-powder and powder-wall are set as 0.34 and 0.20, respectively. It is very difficult to measure the rolling friction coefficient of fine powder. The Al2O3 ceramic is a typically hard and brittle material. The rolling resistance of Al2O3 ceramic powder is very weak. Therefore, the rolling friction coefficient is set as a small value of 0.05. The value of the restitution coefficient can reference literature [31]. According to the literature [15], the size of the DEM system has little effect on the angle of repose (AOR). The AOR defined in ISO-4490 [34] is a common technique, where powder flows through the funnel onto the plate and the slope angle of the developed cone to the base plate is the AOR and considered as a measure for powder flowability. The surface energy density γ of particle-particle and particle-wall is calibrated by comparing the simulations and experiments of the AOR, as shown in Figure 3a,b. Finally, simulations and experimental results of the avalanche angle fit well as shown in Figure 3c,d. Table 2 provides the DEM parameters of particles.

2.3. Simulation Conditions

The roller-spreading simulations include a feeding cylinder and a build platform. The width W and length L of the feeding cylinder and the build platform are 2 and 10 mm, as shown in Figure 4. The roller-spreading model is established as follows: initially, 10,000 particles are generated to form the powder-bed on the top of the feeding cylinder. Subsequently, the roller moves from left to right at speed Vs while it rotates counter-clockwise direction at speed ω. The powders are transported into the build platform forming powder layer by the roller. Table 3 shows the value of principal roller-spreading parameters. All of the simulation results presented in this work were carried out with DEM software, EDEM©.

2.4. Characterization of Powder-Bed Density

It is necessary to first characterize powder-bed density in order to study how roller-spreading parameters affect powder-bed density. The build platform is divided into 80 (4 × 20) grids as shown in Figure 3. The size of each grid is Lx (500 μm) × Ly (500 μm) × Lz, where Lz represents the thickness of powder layer formed in individual spreading process. We define average fraction φ ¯ and fraction standard deviation S of grids to characterize the powder-bed compactness density and the uniformity of powder-bed density, respectively. The calculation fraction of powder layer in the grid is given:
φ j = i n 4 3 π R i 3 L x L y L z
where R i is the radius of particle i, Lx (500 μm) × Ly (500 μm) are the lateral dimensions of the grid; n is the total number of particles in the j-th grid.
The average fraction φ ¯ of grids is:
φ ¯ = j = 1 N φ j N
where N is total number of grids.
The fraction standard deviation S of grids can be expressed as:
S = 1 N 1 i = j N ( φ j φ ¯ ) 2

3. Results and Discussion

Spreading powders are divided into three zones according to powder velocity or its normal force state, as shown in Figure 5.
(1)
Powder free-flow zone: the powder flow under the action of weight force and external force exerted by roller. The normal contact forces between particles are weak, but particles velocities are relatively higher in this zone. The particles in the zone that will fall into the build platform to form powder layer.
(2)
Powder compression zone: the powder-bed is compressed and deformed forming a considerably dense powder layer by roller. The particle normal contact force is strong and particles below the roller are squeezed in this zone with low velocity.
(3)
Powder layer zone: powder layer underneath the surface of the roller is formed and flattened by the translation motion of the roller. The particles in this zone keep stable.

3.1. The Effect of Roller’s Translational Velocity on Powder-Bed Density

Figure 6a shows the effect of translational velocity Vs on average fraction φ ¯ . The average fraction φ ¯ decreases with a linear trend from 0.51 to 0.43 with the Vs increasing. The fraction standard deviation S increases from 0.048 to 0.084 as the translational velocity increases, and the increasing trend becomes more significant when the Vs is greater than 120 mm/s, as shown in Figure 6b. The results show that the increased Vs lead to a poor powder-bed compactness density and uniformity of powder-bed density.
When the roller’s translational velocity is 160 mm/s, the powder will fly up leading to particle splash, which will lead to a poor powder-bed density, as shown in Figure 7. In addition, powder splash can contaminate the key parts of the AM machine. We have performed statistical analysis on the number of particles in the powder layer zone in order to further analyze the effect of Vs on the powder layer zone. Figure 8 shows that the number of particles decreases by 20% from 17,908 to 14,162 in powder layer zone when the Vs increases from 40 to 160mm/s. When the Vs increases, particles velocity also increases in the powder deformation zone, as shown in Figure 9. The high velocity of particles in the powder deformation zone crossing the gap between the roller and substrate will keep moving for a longer distance before residing on the substrate, leading to the number of particles decreasing in the powder layer zone [35]. Therefore, the increased Vs lead to a poor powder-bed density. Comparing the results of work [13], increasing the translational velocity of blade also leads to poor powder-bed quality.
According to the above results, the translation velocity should be as low as possible for improving the powder-bed density. However, this is in conflict with the production efficiency in engineering. In AM processes, the part is built up by layer-by-layer, and the translation velocity mainly determines the powder spreading time. It takes much time if the translation velocity is too low. Therefore, the minimum translational velocity should be limited to properly balance the powder-bed density and processing efficiency.

3.2. The Effect of Roller’s Rotational Speed on Powder-Bed Density

The ω of roller was set from 40 to 320 rpm in order to investigate the effect of roller’s rotational speed ω on powder-bed density, and Figure 10 shows the effect of ω on powder-bed density. With the increase of ω, the value of φ ¯ fluctuates around 0.49, which shows no significant changes. It can be seen from Figure 10b that the fraction standard deviation S also fluctuates around 0.058. The results show that the ω has little effect on the powder-bed density.
Figure 11 shows powder flow velocity under different roller’s rotational speed. When the particles near the roller reach the highest point under its action, particles begin to slide down entering the free-flow zone due to gravity. Particles velocity in the free-flow zone at 40 rpm is lower, while particles velocity shows no significant difference in the free-flow zone when ω increases from 80 to 320 rpm. In the powder deformation zone, particles velocity shows no significant difference when the ω increases from 40 to 320 rpm. Therefore, the ω has little effect on the powder-bed density. Usually, the higher particles velocity can overcome the adhesion energy of particles due to the greater kinetic energy. The powder flowability in the free-flow zone will be improved, which is easy to fall into the platform. It is known that adhesive and cohesive forces of fine powders lead to mostly unwanted agglomeration and poor powder flowability. Therefore, the higher ω is recommended in order to overcome poor powder flowability.

3.3. The Effect of Roller’s Diameter on Powder-Bed Density

This paper carried out DEM simulation on the roller-spreading process under different roller’s diameters in order to study the effect of roller structure parameters on the powder-bed density. It can be seen from Figure 12a that the φ ¯ increases from 0.47 to 0.50, when the diameter of the roller increases from 3 to 8 mm. It can be seen from Figure 12b that the fraction standard deviation S decreases. The results show that a larger roller’s diameter can get better powder-bed density.
Powders are divided into three zones according to powder velocity or normal contact force state, as shown in Figure 5. Figure 5b shows that the force in the particle system is disorganized. We can see that the particles are compacted by load of the roller. The force is attenuating from powder compression zone to powder free-flow zone. The particles are compressed and deformed forming a relatively dense powder layer in the powder compression zone. Therefore, the force in the powder compression zone is stronger than other zones. The number of normal contacts is quantitatively investigated in the powder free-flow zone and the powder compression zone, as shown in Figure 13. The high velocity powders with weak normal contact force in the range of 1.0 × 10−13 N to 1.0 × 10−7 N in the powder free-flow zone, as shown in Figure 5b. The strong normal contact force bigger than 1.0 × 10−7 N in the powder compression zone is the main driving force for the densification of the powder-bed. Figure 14 shows the cumulative number of particle normal contacts when normal contact force bigger than 1.0 × 10−7 N in the powder compression zone under different roller’s diameters. When D is 3 mm, the maximum magnitude of particle normal contact force is 6.0 × 10−6 N and number of particle normal contacts is 16,000. When D increases, the particle normal contact force is strengthening and the number of particle normal contacts is 18,000 when D = 8 mm. The contact zone between the roller and powder in the powder compression zone will increase. This is why increasing the roller’s diameter can lead to the densification of powder-bed.
The force of roller is discussed in the following in order to further analyze the force between the powder and roller. When the roller begins to enter the build platform, the total force of the roller in the vertical directions Fz is collected as shown in Figure 15. We can see that Fz show great changes when the roller moves over build platform. Obviously, the simulation results show the Fz increase with the diameter of roller, which means greater pressure on powder-bed. Therefore, the diameter of roller should be increased within a certain range in order to increase the powder-bed density.

3.4. The Effect of Powder Layer Thickness on Powder-Bed Density

Usually, the layer thickness of roller-spreading is from 50 to 175 μm. Figure 16a shows that the φ ¯ varies from 0.15 to 0.52 when the layer thickness increases from 50 to 175 μm. It can be seen from Figure 16b that the fraction standard deviation S is reduced as increasing the powder layer thickness. This finding suggests that the layer thickness is the most important factor that affects the powder-bed density.
When the layer thickness is 50 and 75 μm, the powder layer will have a lot of voids, as shown in Figure 17a,b. According to the literature [12], particle jamming becomes significantly stronger in the region before the blade when the layer thickness is too small. When the H is 50 and 75 μm, the particles are difficult to fall into the build platform due to particle jamming. Most of the particles are pushed out of the build platform by the roller. The powder-bed will have a lot of voids and this will affect the mechanical properties of the AM parts. When the powder layer thickness became smaller, the layer thickness has little effect on average fraction, as shown in the Figure 16a. However, when the layer thickness is greater than 150 μm, the powder-bed becomes more uniform and denser, as shown in the Figure 17c. This finding indicates that increasing layer thickness can obtain a good powder-bed compactness density and uniformity of powder-bed density.
From the perspective of powder-bed density, the layer thickness should be large. However, large layer thickness means that more particles need be bonded in the Inkjet 3D Printing process or melt in SLM processes, which leads to rough surface density and poor mechanical properties [36]. Roller-spreading is a key process in AM processes, and it needs to be related to the subsequent processing technology. Therefore, it is necessary to properly balance powder-bed density and the quality of AM parts.

4. Conclusions

In this work, DEM is used to simulate Al2O3 ceramic powder spreading by roller. The Hertz–Mindlin with JKR model of Al2O3 ceramic powder is established and calibrated. The powders are divided into three zones: I. Powder free-flow zone, II. Powder compression zone, and III. Powder layer zone. The effect of powder spreading parameters, including translational velocity Vs, roller’s rotational speed ω, roller’ diameter D, and thickness of powder layer H are on Al2O3 ceramic powder-bed density are analyzed.
(1)
The increased roller’s translational velocity lead to a reduced number of particles in powder layer zone. This will lead to poor powder-bed density. The roller’s rotational speed has little effect on the powder-bed density.
(2)
When the roller’s diameter increases, more particles are in the powder compression zone. The normal contact force becomes strong, which enhances densification of powder-bed.
(3)
Layer thickness is the most significant influencing factor on powder-bed density. When the layer thickness is 50 μm, most of the particles are pushed out of the build platform, forming a lot of voids. However, when the layer thickness is greater than 150 μm, the powder-bed becomes more uniform and denser.
The effect of roller-spreading parameters on the powder-bed density is studied in this paper. However, this paper does not optimize powder spreading parameters. In future work, an optimization model will be developed in order to improve the powder-bed density. The optimization model will be verified by experiments, which will provide a reliable basis for roller-spreading parameters optimization selection in AM.

Author Contributions

J.Z. and Y.T. performed and analyzed the simulations and wrote the manuscript with support from T.B. and Y.X. conducted the experiments. X.X. and S.J. supervised the whole work and discussed all results. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Fund of China (NSFC) (Grant No 11772135), High-level Personnel Fund of Huaqiao University (Z16X0110), Project for Postgraduates’ Innovative Fund in Scientific Research of Huaqiao University (18011080013), the Scientific Research Funds of Huaqiao University (605-50Y19023) and the Scientific Research Fund of Hunan Provincial Education Department (18B074).

Acknowledgments

The authors would like to thank the editors, reviewers, and Lai for their comments and suggestions, which are very valuable in improving the quality of our paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brandão, A.; Gerard, R.; Gumpinger, J.; Beretta, S.; Makaya, A.; Pambaguian, L.; Ghidini, T. Challenges in Additive Manufacturing of Space Parts: Powder Feedstock Cross-Contamination and its Impact on End Products. Materials 2017, 10, 522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Rausch, A.; Küng, V.; Pobel, C.; Markl, M.; Körner, C. Predictive Simulation of Process Windows for Powder Bed Fusion Additive Manufacturing: Influence of the Powder Bulk Density. Materials 2017, 10, 1117. [Google Scholar] [CrossRef] [Green Version]
  3. Ali, U.; Mahmoodkhani, Y.; Imani Shahabad, S.; Esmaeilizadeh, R.; Liravi, F.; Sheydaeian, E.; Huang, K.Y.; Marzbanrad, E.; Vlasea, M.; Toyserkani, E. On the Measurement of Relative Powder-Bed Compaction Density in Powder-Bed Additive Manufacturing Processes. Mater. Des. 2018, 155, 495–501. [Google Scholar] [CrossRef]
  4. Zocca, A.; Colombo, P.; Gomes, C.M.; Günster, J. Additive Manufacturing of Ceramics: Issues, Potentialities, and Opportunities. J. Am. Ceram. Soc. 2015, 98, 1983–2001. [Google Scholar] [CrossRef]
  5. Wischeropp, T.M.; Emmelmann, C.; Brandt, M.; Pateras, A. Measurement of Actual Powder Layer Height and Packing Density in a Single Layer in Selective Laser Melting. Addit. Manuf. 2019, 28, 176–183. [Google Scholar] [CrossRef]
  6. Chen, H.; Wei, Q.; Wen, S.; Li, Z.; Shi, Y. Flow Behavior of Powder Particles in Layering Process of Selective Laser Melting: Numerical Modeling and Experimental Verification Based on Discrete Element Method. Int. J. Mach. Tools Manuf. 2017, 123, 146–159. [Google Scholar] [CrossRef]
  7. Zohdi, T.I. A Direct Particle-Based Computational Framework for Electrically Enhanced Thermo-Mechanical Sintering of Powdered Materials. Math. Mech. Solids 2013, 19, 93–113. [Google Scholar] [CrossRef]
  8. Zohdi, T.I. Additive Particle Deposition and Selective Laser Processing-A Computational Manufacturing Framework. Comput. Mech. 2014, 54, 171–191. [Google Scholar] [CrossRef]
  9. Steuben, J.C.; Iliopoulos, A.P.; Michopoulos, J.G. Discrete Element Modeling of Particle-Based Additive Manufacturing Processes. Comput. Method. Appl. M. 2016, 305, 537–561. [Google Scholar] [CrossRef] [Green Version]
  10. Haeri, S. Optimisation of Blade Type Spreaders for Powder Bed Preparation in Additive Manufacturing Using DEM Simulations. Powder Technol. 2017, 321, 94–104. [Google Scholar] [CrossRef] [Green Version]
  11. Lee, Y.S.; Nandwana, P.; Zhang, W. Dynamic Simulation of Powder Packing Structure for Powder Bed Additive Manufacturing. Int. J. Adv. Manuf. Technol. 2018, 96, 1–14. [Google Scholar] [CrossRef]
  12. Nan, W.; Ghadiri, M. Jamming during particle spreading in additive manufacturing. Powder Technol. 2018, 338, 253–262. [Google Scholar] [CrossRef]
  13. Parteli, E.J.R.; Pöschel, T. Particle-Based Simulation of Powder Application in Additive Manufacturing. Powder Technol. 2016, 288, 96–102. [Google Scholar] [CrossRef]
  14. Chen, H.; Wei, Q.; Zhang, Y.; Chen, F.; Shi, Y.; Yan, W. Powder-spreading mechanisms in powder-bed-based additive manufacturing: Experiments and computational modeling. Acta Mater. 2019, 179, 158–171. [Google Scholar] [CrossRef]
  15. Han, Q.; Gu, H.; Setchi, R. Discrete Element Simulation of Powder Layer Thickness in Laser Additive Manufacturing. Powder Technol. 2019, 352, 91–102. [Google Scholar] [CrossRef]
  16. Haeri, S.; Wang, Y.; Ghita, O.; Sun, J. Discrete Element Simulation and Experimental Study of Powder Spreading Process in Additive Manufacturing. Powder Technol. 2017, 306, 45–54. [Google Scholar] [CrossRef] [Green Version]
  17. Meier, C.; Weissbach, R.; Weinberg, J.; Wall, W.A.; John Hart, A. Modeling and Characterization of Cohesion in Fine Metal Powders with a Focus on Additive Manufacturing Process Simulations. Powder Technol. 2019, 343, 855–866. [Google Scholar] [CrossRef] [Green Version]
  18. Tan, Y.; Xiao, X.; Zhang, J.; Jiang, S. Determination of discrete element model contact parameters of nylon powder at SLS preheating temperature and its flow characteristics. Chin. J. Theor. Appl. Mech. 2019, 51, 56–63. [Google Scholar]
  19. Nan, W.; Pasha, M.; Ghadiri, M. Numerical simulation of particle flow and segregation during roller spreading process in additive manufacturing. Powder Technol. 2020, 364, 811–821. [Google Scholar] [CrossRef]
  20. Cundall, P.A.; Strack, O.D.L. Discussion: A discrete numerical model for granular assemblies. Géotechnique 1980, 30, 331–336. [Google Scholar] [CrossRef] [Green Version]
  21. Johnson, K.; Kendall, K. Surface Energy and the Contact of Elastic Solids. Proceedings of the Royal Society A: Mathematical. Phys. Eng. Sci. 1971, 324, 301–313. [Google Scholar]
  22. Di Renzo, A.; Di Maio, F.P. Comparison of Contact-Force Models for the Simulation of Collisions in DEM-based Granular Flow Codes. Chem. Eng. Sci. 2004, 59, 525–541. [Google Scholar] [CrossRef]
  23. Cleary, P.W. Predicting Charge Motion, Power Draw, Segregation and Wear in Ball Mills Using Discrete Element Methods. Miner. Eng. 1998, 11, 1061–1080. [Google Scholar] [CrossRef]
  24. Gilabert, F.A.; Roux, J.; Castellanos, A. Computer Simulation of Model Cohesive Powders: Influence of Assembling Procedure and Contact Laws on Low Consolidation States. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2007, 75, 011303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Mindlin, R.D. Vibrations of doubly-rotated-cut quartz plates with monoclinic symmetry. Int. J. Solids Struct. 1985, 21, 597–607. [Google Scholar] [CrossRef]
  26. Xiao, X.; Tan, Y.; Zhang, H.; Deng, R.; Jiang, S. Experimental and DEM studies on the Particle Mixing Performance in Rotating Drums: Effect of Area Ratio. Powder Technol. 2017, 314, 182–194. [Google Scholar] [CrossRef]
  27. Coetzee, C.J. Review: Calibration of the discrete element method. Powder Technol. 2017, 310, 104–142. [Google Scholar] [CrossRef]
  28. Peng, Z.; Luo, X.; Xie, Z.; An, D.; Yang, M. Effect of print path process on sintering behavior and thermal shock resistance of Al2O3 ceramics fabricated by 3D inkjet-printing. Ceram. Int. 2018, 44, 16766–16772. [Google Scholar] [CrossRef]
  29. Wu, D.; Liu, H.; Lu, F.; Ma, G.; Yan, S.; Niu, F.; Guo, D. Al2O3-YAG eutectic ceramic prepared by laser additive manufacturing with water-cooled substrate. Ceram. Int. 2019, 45, 4119–4122. [Google Scholar] [CrossRef]
  30. Hu, Y.; Ning, F.; Cong, W.; Li, Y.; Wang, X.; Wang, H. Ultrasonic Vibration-Assisted Laser Engineering Net Shaping of ZrO2-Al2O3 Bulk Parts: Effects on Crack Suppression, Microstructure, and Mechanical Properties. Ceram. Int. 2018, 44, 2752–2760. [Google Scholar] [CrossRef]
  31. Tahir, A.; Rasche, S.; Könke, C. Discrete Element Model Development of ZTA Ceramic Granular Powder Using Micro Computed Tomography. Adv. Powder Technol. 2018, 29, 3471–3482. [Google Scholar] [CrossRef]
  32. Hanley, K.J.; O’Sullivan, C. Analytical Study of the Accuracy of Discrete Element Simulations. Int. J. Numer. Meth. Eng. 2017, 109, 29–51. [Google Scholar] [CrossRef] [Green Version]
  33. ASTM_International. ASTM D6773-08, Standard Shear Test Method for Bulk Solids Using the Schulze Ring Shear Tester; ASTM_International: West Conshohocken, PA, USA, 2008. [Google Scholar]
  34. International Standardisation Organisation. ISO 4490:2014—Metallic Powders—Determination of Flow Rate by Means of a Calibrated Funnel (Hall Flowmeter); ISO: Geneva, Switzerland, 2014. [Google Scholar]
  35. Chen, H.; Chen, Y.; Liu, Y.; Wei, Q.; Shi, Y.; Yan, W. Packing quality of powder layer during counter-rolling-type powder spreading process in additive manufacturing. Int. J. Mach. Tools Manuf. 2020, 153, 103553. [Google Scholar] [CrossRef]
  36. Zheng, Y.; Zhang, K.; Liu, T.T.; Liao, W.H.; Zhang, C.D.; Shao, H. Cracks of Alumina Ceramics by Selective Laser Melting. Ceram. Int. 2019, 45, 175–184. [Google Scholar] [CrossRef]
Figure 1. Scanning electron microscope (SEM) micrograph of Al2O3 ceramic powders.
Figure 1. Scanning electron microscope (SEM) micrograph of Al2O3 ceramic powders.
Materials 13 02285 g001
Figure 2. Shear test results of powder-powder and powder-wall.
Figure 2. Shear test results of powder-powder and powder-wall.
Materials 13 02285 g002
Figure 3. Comparing simulations and experiments: (a) Simulations results of angle of repose (AOR); (b) Experimental results of AOR; (c) Simulations results of avalanche angle; and, (d) Experimental results of avalanche angle.
Figure 3. Comparing simulations and experiments: (a) Simulations results of angle of repose (AOR); (b) Experimental results of AOR; (c) Simulations results of avalanche angle; and, (d) Experimental results of avalanche angle.
Materials 13 02285 g003
Figure 4. DEM simulation of roller-spreading process.
Figure 4. DEM simulation of roller-spreading process.
Materials 13 02285 g004
Figure 5. Powders are divided into three zones when D = 5 mm, Vs = 80 mm/s, ω = 200 rpm, and H = 150 μm at t = 0.28 s.
Figure 5. Powders are divided into three zones when D = 5 mm, Vs = 80 mm/s, ω = 200 rpm, and H = 150 μm at t = 0.28 s.
Materials 13 02285 g005
Figure 6. The effect of roller’s translational velocity on powder-bed density when D = 5 mm, ω = 200 rpm, H =150 μm.
Figure 6. The effect of roller’s translational velocity on powder-bed density when D = 5 mm, ω = 200 rpm, H =150 μm.
Materials 13 02285 g006
Figure 7. The powder flow pattern when Vs is 160 mm/s at t = 0.14 s.
Figure 7. The powder flow pattern when Vs is 160 mm/s at t = 0.14 s.
Materials 13 02285 g007
Figure 8. Number of particles in powder layer zone.
Figure 8. Number of particles in powder layer zone.
Materials 13 02285 g008
Figure 9. Powder flow velocity under different roller’s translational velocity.
Figure 9. Powder flow velocity under different roller’s translational velocity.
Materials 13 02285 g009
Figure 10. The effect of roller’ rotational speed on powder-bed density when D = 5 mm, Vs = 80 mm/s, H = 150 μm.
Figure 10. The effect of roller’ rotational speed on powder-bed density when D = 5 mm, Vs = 80 mm/s, H = 150 μm.
Materials 13 02285 g010
Figure 11. Powder flow velocity under different roller’s rotational speed at t= 0.28 s.
Figure 11. Powder flow velocity under different roller’s rotational speed at t= 0.28 s.
Materials 13 02285 g011
Figure 12. The effect of roller’s diameters on powder-bed density when Vs = 80 mm/s, ω = 200 rpm, H = 150 μm.
Figure 12. The effect of roller’s diameters on powder-bed density when Vs = 80 mm/s, ω = 200 rpm, H = 150 μm.
Materials 13 02285 g012
Figure 13. Number distribution of the normal particle contact force when D = 5 mm, Vs = 80 mm/s, ω = 200 rpm, H = 150 μm at t = 0.28 s.
Figure 13. Number distribution of the normal particle contact force when D = 5 mm, Vs = 80 mm/s, ω = 200 rpm, H = 150 μm at t = 0.28 s.
Materials 13 02285 g013
Figure 14. Cumulative numerical distribution of the normal particle contact force bigger than 1.0 × 10−7 N in the powder compression zone under different roller’s diameters when Vs = 80 mm/s, ω = 200 rpm, H = 150 μm at t = 0.28 s.
Figure 14. Cumulative numerical distribution of the normal particle contact force bigger than 1.0 × 10−7 N in the powder compression zone under different roller’s diameters when Vs = 80 mm/s, ω = 200 rpm, H = 150 μm at t = 0.28 s.
Materials 13 02285 g014
Figure 15. The total force of the roller in the vertical directions under different roller’s diameters.
Figure 15. The total force of the roller in the vertical directions under different roller’s diameters.
Materials 13 02285 g015
Figure 16. The effect of powder layer thickness on powder-bed density when D = 5 mm, Vs = 80 mm/s, ω = 200 rpm.
Figure 16. The effect of powder layer thickness on powder-bed density when D = 5 mm, Vs = 80 mm/s, ω = 200 rpm.
Materials 13 02285 g016
Figure 17. The powder-bed characteristics under different layer thickness values: (a) H = 50 μm; (b) H = 75 μm; and, (c) H = 150 μm.
Figure 17. The powder-bed characteristics under different layer thickness values: (a) H = 50 μm; (b) H = 75 μm; and, (c) H = 150 μm.
Materials 13 02285 g017
Table 1. List of equations used in the contact model.
Table 1. List of equations used in the contact model.
ParameterEquations
Normal stiffness constant S n = 2 E * R * δ n
Damping ratio β = I n e I n 2 e + π 2
Tangential stiffness S t = 8 G R δ n
Equivalent shear modulus G = [ ( 1 γ i 2 ) G i + ( 1 γ j 2 ) G j ] 1
Equivalent Young’s Modulus E * = [ ( 1 γ i 2 ) E i + ( 1 γ j 2 ) E j ] 1
Equivalent mass m * = [ 1 m i + 1 m j ] 1
Equivalent radius R * = [ 1 R i + 1 R j ] 1
Table 2. Discrete Element Method (DEM) parameters of Al2O3 ceramic particles.
Table 2. Discrete Element Method (DEM) parameters of Al2O3 ceramic particles.
ParameterValue
Particle density (kg/m3)3820
Particle size D (50) (μm)48
Particle shear modulus (GPa)3
Particle Poisson’s ratio0.3
Roller density (kg/m3)7800
Roller shear modulus (GPa)80
Roller Poisson’s ratio0.3
Static friction coefficient of particle-wall0.20
Rolling friction coefficient of particle-wall0.05
Restitution coefficient of particle-wall0.52
Surface energy density of particle-wall (mJ/m2)0.17
Static friction coefficient of particle-particle0.34
Rolling friction coefficient of particle-particle0.05
Restitution coefficient of particle-particle0.50
Surface energy density of particle-particle (mJ/m2)0.15
Table 3. The value of roller-spreading parameters.
Table 3. The value of roller-spreading parameters.
ParametersValue
Roller’s translational velocity Vs (mm/s)40/60/80/100/120/140/160
Roller’s rotational velocity ω (rpm)40/80/120/160/200/240/280/320
Roller’s diameter D (mm)3/4/5/6/7/8
Powder layer thickness H (μm)50/75/100/125/150/175

Share and Cite

MDPI and ACS Style

Zhang, J.; Tan, Y.; Bao, T.; Xu, Y.; Xiao, X.; Jiang, S. Discrete Element Simulation of the Effect of Roller-Spreading Parameters on Powder-Bed Density in Additive Manufacturing. Materials 2020, 13, 2285. https://0-doi-org.brum.beds.ac.uk/10.3390/ma13102285

AMA Style

Zhang J, Tan Y, Bao T, Xu Y, Xiao X, Jiang S. Discrete Element Simulation of the Effect of Roller-Spreading Parameters on Powder-Bed Density in Additive Manufacturing. Materials. 2020; 13(10):2285. https://0-doi-org.brum.beds.ac.uk/10.3390/ma13102285

Chicago/Turabian Style

Zhang, Jiangtao, Yuanqiang Tan, Tao Bao, Yangli Xu, Xiangwu Xiao, and Shengqiang Jiang. 2020. "Discrete Element Simulation of the Effect of Roller-Spreading Parameters on Powder-Bed Density in Additive Manufacturing" Materials 13, no. 10: 2285. https://0-doi-org.brum.beds.ac.uk/10.3390/ma13102285

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop