Microstructures and corresponding properties of compacted powders ultimately depend on the mechanical response of individual particles. In principle, computational simulations can predict the results of powder compaction processes, but the selection of appropriate models for both particle–particle interactions and particle deformations across all relevant length scales remain nontrivial tasks, especially in material systems lacking detailed mechanical property information. The work presented here addresses these issues by conducting uniaxial compressions in situ inside of a scanning electron microscope to characterize the mechanical response of individual micron-sized particles of a molecular crystal, hexanitrohexaazaisowurtzitane (CL-20). This experimental approach enabled the collection of quantitative force and displacement data alongside simultaneous imaging to capture morphology changes. The results reveal information about elastic deformation, yield, plastic deformation, creep, and fracture phenomena. Accordingly, this work demonstrates a generalizable approach for assessing the mechanical response of individual micron-sized molecular crystal particles and utilizing those responses in particle-level models. Graphic abstract: [Figure not available: see fulltext.].
By generalizing a geometric argument for frictionless spheres, a model is proposed for the jamming density φJ of mechanically stable packings of bidisperse, frictional spheres. The monodisperse, μs-dependent jamming density φJmono(μs) is the only input required in the model, where μs is the coefficient of friction. The predictions of the model are validated by robust estimates of φJ obtained from computer simulations of up to 107 particles for a wide range of μs, and size ratios up to 40:1. Although φJ varies nonmonotonically with the volume fraction of small spheres fs for all μs, its maximum value φJ,max at an optimal fmaxs are both μs dependent. The optimal fmaxs is characterized by a sharp transition in the fraction of small rattler particles.
Two key mechanical processes exist in the formation of powder compacts. These include the complex kinematics of particle rearrangement as the powder is densified and particle deformation leading to mechanical failure and fragmentation. Experiments measuring the time varying forces across a densifying powder bed have been performed in powders of microcrystalline cellulose with mean particle sizes between 0.4 and 1.2 mm. In these experiments, diagnostics measured the applied and transmitted loads and the bulk powder density. Any insight into the particle behavior must be inferred from deviations in the smoothly increasing stress-density compaction relationship. By incorporating a window in the compaction die body, simultaneous images of particle rearrangement and fracture at the confining window are captured. The images are post-processed in MATLAB® to track individual particle motion during compression. Complimentary discrete element method (DEM) simulations are presented and compared to experiment. The comparison provides insight into applying DEM methods for simulating large or permanent particle deformation and suggests areas for future study.
Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.
Particle characteristics can drastically influence the process-structure-property-performance aspects of granular materials in compression. We aim to computationally simulate the mechanical processes of stress redistribution in compacts including the kinematics of particle rearrangement during densification and particle deformation leading to fragmentation. Confined compression experiments are conducted with three sets of commercial microcrystalline cellulose particles nearly spherical in shape with different mean particle size. Experimentally measured compression curves from tall powder columns are fitted with the Kenkre et al. (J. of American Chemical Society, Vol. 79, No. 12) model. This model provides a basis to derive several common two-parameter literature models and as a framework to incorporate statistical representations of critical particle behaviors. We focus on the low-stress compression data and the model comparisons typically not discussed in the literature. Additional single particle compressions report fracture strength with particle size for comparison to the apparent particle strength extracted from bulk compression data.
Using random walk analyses we explore diffusive transport on networks obtained from contacts between isotropically compressed, monodisperse, frictionless sphere packings generated over a range of pressures in the vicinity of the jamming transition p→0. For conductive particles in an insulating medium, conduction is determined by the particle contact network with nodes representing particle centers and edges contacts between particles. The transition rate is not homogeneous, but is distributed inhomogeneously due to the randomness of packing and concomitant disorder of the contact network, e.g., the distribution of the coordination number. A narrow escape time scale is used to write a Markov process for random walks on the particle contact network. This stochastic process is analyzed in terms of spectral density of the random, sparse, Euclidean and real, symmetric, positive, semidefinite transition rate matrix. Results show network structures derived from jammed particles have properties similar to ordered, euclidean lattices but also some unique properties that distinguish them from other structures that are in some sense more homogeneous. In particular, the distribution of eigenvalues of the transition rate matrix follow a power law with spectral dimension 3. However, quantitative details of the statistics of the eigenvectors show subtle differences with homogeneous lattices and allow us to distinguish between topological and geometric sources of disorder in the network.