Granular Matter最新文献

The laser-beam powder bed fusion process for metals, commonly abbreviated as PBF-LB/M, is a widely used process for the additive manufacturing of parts. Numerical simulations are useful to identify optimal process parameters for different materials and to obtain detailed insights into process dynamics. The present work uses a single-phase incompressible Smoothed Particle Hydrodynamics (SPH) scheme to model PBF-LB/M which was found to reduce the required computational time and significantly stabilize the partially violent flow in the melt pool in comparison to a weakly compressible SPH approach. The laser-material interaction is realistically modelled by means of a ray tracing method. An approach to model the effective thermal coductivity of the powder bed is proposed. Excellent agreement between the simulation results and experimental X-ray analyses of the transition from conduction melting mode to keyhole mode including geometric properties of the vapor depression zone was found. These results prove the usability of SPH as a high precision simulation tool for PBF-LB/M.

Graphic abstract

The investigation of rockfill materials poses challenges due to their large particle size, associated high cost, and long laboratory testing duration. As a result, empirical correlations based on historical experimental studies are commonly used to design and analyse rockfill structures. However, the extensive use of rockfill in a wide range of applications and limited understanding of its mechanical behaviour emphasize the need for further research. These make it necessary to develop a robust technique capable of capturing key parameters such as particle shape and breakage, allowing for the simulation and study of large-scale assemblies with realistic boundary conditions. Given that the behaviour of rockfill is highly scale-dependent, primarily due to particle breakage, the simplified laboratory tests on the scaled-down assemblies can be misleading. Particle breakage is a fundamental phenomenon in the mechanical behaviour of rockfill and significantly affects shear strength, deformability, and porosity under different stress levels. The particle breakage is influenced by factors such as the rockfill’s maximum particle size, mineralogy, particle shape, gradation, and confining stresses. This study adopts a computationally efficient breakage method called the Modified Particle Replacement Method (MPRM) based on the Discrete Element Method. A Tile-Based Flexible Membrane (TBFM) for triaxial test modelling has been developed by employing segmental rectangular walls to create a deformable membrane. The effects of critical parameters, including particle shape, confining stress, membrane resolution, degree of flexibility, and the characteristic strength of the particles, are examined. The findings of the combined MPRM-TBFM approach demonstrate the significant influence of membrane flexibility on volumetric-related behaviour.

Graphical Abstract

Laser Powder Bed Fusion (L-PBF) has immense potential for the production of complex, lightweight, and high-performance components. The traditional optimization of process parameters is costly and time-intensive, due to reliance on experimental approaches. Current numerical analyses often model single-line scans, while it is necessary to model multiple fully scanned layers to optimize for bulk material quality. Here, we introduce a novel approach utilizing discrete element simulations with a ray tracing-modeled laser heat source. Our approach significantly reduces the cost and time consumption compared to conventional optimization methods. GPU acceleration enables efficient simulation of multiple layers, resulting in parameters optimized for bulk material. In a case study, parameters were optimized for AlSi10Mg in just 5 days, a process that would have taken over 8 months without GPU acceleration. Experimental validation affirms the quality of the optimized process parameters, achieving an optical density of 99.91%.

Graphical Abstract

Optimization using the accelerated simulation yielded an optimized parameter set within 5 days. This resulted in apart with an optical density of 99.91%.

We consider a dense aggregate of elastic, frictional particles isotropically compressed and next uniaxial strained at constant pressure. We show how failure can be predicted if fluctuations in the kinematics of contacting particles are introduced. We focus on the second order work and the possibility that at some stressed states it becomes negative under proper perturbations. Our analysis involves both a theoretical model and numerical simulations based upon the distinct element method (DEM). The theoretical model deals with contacting particles with incremental relative displacements that deviate from the average deformation in order to ensure their equilibrium. Because of this, the macroscopic stiffness tensor of the aggregate, that relates increments in stress with increments in strain, does not have the major symmetry. Consequently, in the hardening regime, we predict stressed states in which the second order work vanishes. The model seems transparent, and it makes clear and illustrative the role played by the fluctuations introduced in the kinematics of contacting particles in relation to the vanishing of second order work in an aggregate of compressed particles. The comparison with numerical simulations data supports the model.

Graphical Abstract

Statistical representation of the aggregate: conditional average.

We introduce a novel numerical method for the simulation of soft granular materials, in which the particles can undergo large strains under load without rupture. The proposed approach combines an explicit total Lagrangian formulation of Material Point Method (TLMPM) with the Contact Dynamics (CD) method. The TLMPM resolves particle bulk deformations whereas the CD treats contact interactions between soft particles. The efficiency and accuracy of this approach are illustrated by analyzing diametral compression of a soft circular particle and the compaction of an assembly of soft particles up to very high levels of packing fraction. We show that although the assembly undergoes a jamming transition, the particles continue to rearrange as they get increasingly distorted under load. Interestingly, as the packing fraction increases, a transition occurs to a regime fully governed by particle shape change. The evolution of the global stress as well as the connectivity of the particles as a function of the packing fraction are discussed and a predictive model relating stress to packing fraction beyond jamming transition is proposed.

Graphical Abstract

The softening effect has been widely accepted as the fundamental mechanical property of the granular materials, which underlies some specific phenomena such as fluidization during vibration. In this paper, a series of resonance column experiments are performed to observe the modulus softening of granular materials. A statistical softening model is subsequently proposed and its applicability is verified through a quantitative analysis of the variation of the normalized modulus by changing the external confining pressure. The average potential energy in grain contact has been found to be a power-law scaling with grain size. An evolution model is further implemented to account for the experimental findings on the evolution of modulus of the granular system subjected to different confining pressures. The modulus evolution, including softening and recovery, can be captured by the unified evolution model.

Graphical Abstract

Shear modulus evolution

The focus of this study is the experimental characterization of cemented granular materials, with the aim of identifying the microscopic properties of the solid bonds and describing the extension to macroscopic mechanical strength of cemented samples. We chose to use artificially bonded granular materials, made of glass beads connected by solid paraffin bridges. The results of several sets of laboratory tests at different scales are presented and discussed. Micromechanical tests investigate the yield strength of single solid bonds between particles under traction, shearing, bending and torsion loading, as a function of variations in particle size, surface texture and binder content. Macro-scale tensile tests on cemented samples explore then the scale transition, including influence of confining walls through homothetic variations of the sample size. Despite the large statistical dispersion of the results, it was possible to derive and validate experimentally an analytical expression for micro tensile yield force as a function of the binder content, coordination number and grain diameter. In view of the data, an adhesive bond strength at the contact between bead and solid bond is deduced with very good accuracy and it is even reasonable to assume that the other threshold values (shear force, bending and torsion moments) are simply proportional to the tensile yield, thus providing a comprehensive 3D model of cemented bond. However, the considerable dispersion of the data at the sample scale prevents validation of the extended model for macroscopic yield stress. A final discussion examines the various factors that may explain intrinsic variability. By comparison with other more realistic systems studied in the literature in the context of bio-cementation, our artificial material nevertheless appears suitable for representing a cemented granular material. Being easy to implement, it could thus enable the calibration of discrete cohesion models for simulation of practical applications.

In this paper, a series of true triaxial tests with different intermediate principal stress ratios are conducted on both the lunar soil simulant and the sandy soils on earth using the discrete element method. An advanced discrete element servomechanism based on polyhedral specimen configuration is implemented such that true triaxial loading paths can be implemented under low confining pressure without introducing severe stress concentration. The high frictional angle and apparent cohesion of the lunar simulant are captured by employing a highly efficient contact model that fuses rolling resistance and van der Waals forces. The employed micro-scale parameters are calibrated based on the triaxial test results of the CSU-LRS-1 lunar soil simulant. The simulation results show that the lunar soil simulant exhibits lower shear strength with an increasing intermediate principal stress ratio. Generally, although the lunar soil simulant has a greater void ratio than that of sandy soils, the former exhibits significantly stronger shear-induced dilatancy and higher shear strength. The evolution of the load-bearing structure is quantified through a contact-normal-based fabric tensor. The interplay between internal structure evolution and external loadings can well explain the difference in mechanical behavior between lunar soil simulant and sandy soils on earth.

The collapse behaviour of granular materials is influenced by many factors, such as aspect ratio and inter-particle friction. However, the specific impact of basal to grain friction on column collapse remains poorly understood. In this study, we systematically analyse the effect of basal friction on gravity-driven granular column collapse using a validated smoothed particle hydrodynamics (SPH) model. The results show that such the basal friction coefficient does influence deposit geometry, deposit morphology, and energy conversion. To predict the run-out distance, we propose a modified formula that incorporates the basal friction coefficient, considering two extreme cases, i.e., μ = 0 and + ∞. The basal friction also exerts an influence on the final height, with higher friction coefficients resulting in greater final heights. As the friction coefficient increases, the aspect ratio corresponding to the maximum final height also increase. However, we observe a convergence of the effect of basal friction on the final height when μ > 0.5. Furthermore, the competition mechanism between the initial column aspect ratio and basal friction coefficient reveals two transition zones between the three main deposit regimes (regime I, regime II, and regime III). This suggests that the deposit regime can be influenced by basal friction. Additionally, an analysis of energy conversion supports many of the conclusions provided in the text and exhibits the interplay between pressure gradient and base friction. Our findings show the clear influence of basal friction on the collapse behaviour of granular materials and therefore should be carefully considered in future studies.

Graphical Abstract