Inferring the Material Properties of Granular Media for Robotic Tasks

by   Carolyn Matl, et al.

Granular media (e.g., cereal grains, plastic resin pellets, and pills) are ubiquitous in robotics-integrated industries, such as agriculture, manufacturing, and pharmaceutical development. This prevalence mandates the accurate and efficient simulation of these materials. This work presents a software and hardware framework that automatically calibrates a fast physics simulator to accurately simulate granular materials by inferring material properties from real-world depth images of granular formations (i.e., piles and rings). Specifically, coefficients of sliding friction, rolling friction, and restitution of grains are estimated from summary statistics of grain formations using likelihood-free Bayesian inference. The calibrated simulator accurately predicts unseen granular formations in both simulation and experiment; furthermore, simulator predictions are shown to generalize to more complex tasks, including using a robot to pour grains into a bowl, as well as to create a desired pattern of piles and rings. Visualizations of the framework and experiments can be viewed at t=11s



There are no comments yet.


page 1

page 4

page 5

page 6


Mobile 3D Printing Robot Simulation with Viscoelastic Fluids

The system design and algorithm development of mobile 3D printing robots...

STReSSD: Sim-To-Real from Sound for Stochastic Dynamics

Sound is an information-rich medium that captures dynamic physical event...

DiSECt: A Differentiable Simulator for Parameter Inference and Control in Robotic Cutting

Robotic cutting of soft materials is critical for applications such as f...

Learning physics-informed simulation models for soft robotic manipulation: A case study with dielectric elastomer actuators

Soft actuators offer a safe and adaptable approach to robotic tasks like...

DiSECt: A Differentiable Simulation Engine for Autonomous Robotic Cutting

Robotic cutting of soft materials is critical for applications such as f...

Modelling Arbitrary Complex Dielectric Properties – an automated implementation for gprMax

There is a need to accurately simulate materials with complex electromag...

Deep Probabilistic Surrogate Networks for Universal Simulator Approximation

We present a framework for automatically structuring and training fast, ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

A granular material is a collection of discrete macroscopic particles that primarily experience inelastic collisions and are usually unaffected by temperature [17]. Naturally, this definition spans a wide variety of materials, from sand and stones to powders and pills to seeds and cereals, making granular material one of the most manipulated substances in the construction, pharmaceutical, agriculture, and food industries [38]. As robots continue to permeate these industries, they will inevitably interact with granular materials in a variety of settings, such as pouring flour in the kitchen, scattering seeds on a farm, or traversing loose terrain [23]. Despite the ubiquity of granular materials, much of the physics behind their complex behaviors remains a mystery. Granular materials are often considered their own category of matter (separate from solids, liquids, and gases), and consequently, new constitutive laws must be defined. While recent work has delivered promising advances in the formulation of continuum models for granular materials [18, 3, 15, 33, 20, 19], researchers frequently rely on numerical techniques and simulations to corroborate and extend their predictions. The state-of-the-art numerical technique for simulating granular materials is the discrete element method (DEM), where each particle is represented independently. Although DEM has primarily been used in physics, it has also been leveraged by the robotics community for applications such as locomotion on sand or gravel [29, 22, 32] and motion planning for pouring [13]. However, calibrating the material parameters of the simulated grains is still an open problem. Most DEM calibration methods involve a tedious cycle of performing small-scale materials tests, manually tuning parameters, and iterating until the desired experimental behavior is captured. To reduce this burden, the physics community has posed the following challenge: how can particle-scale parameters (e.g., frictional properties) be determined from observations of macroscale behavior (e.g., granular pile shape) [1]? Answering this could lead to more efficient methods for refining simulations to match real macroscale behavior.

Fig. 1: Granular parameter inference framework. (Bottom branch): Multiple granular parameter sets are randomly sampled from a prior distribution , and for each, summary statistics are extracted from a simulation . The set of are used to learn a conditional density function for the parameters using BayesSim. (Top branch): The summary statistics of a depth image of a real granular formation are extracted. With the conditional density function, a posterior is estimated to infer the likeliest simulation granular parameters to match the macroscopic granular behavior.

This paper addresses this question by integrating a framework for likelihood-free Bayesian inference (BayesSim [37]) with a fast and efficient physics simulation platform capable of DEM simulations (NVIDIA’s Isaac Simulator [2]), as well as experimental observations using off-the-shelf depth cameras. Using this software and hardware framework, depicted in Figure 1, particle-scale material properties (i.e., coefficients of friction and restitution) of granular materials are inferred from simple, inexpensive observations of macroscale behavior in the real world. The calibrated simulator can then be leveraged to predict unseen granular formations in both simulation and experiment. This paper concludes by demonstrating how this research may be applied to robotics tasks by using the developed framework on an industrial robot to manipulate agricultural grains (i.e., couscous and barley). The main contributions of this paper are

  • A software and hardware framework that integrates likelihood-free Bayesian inference, an efficient physics simulator, and depth images of granular formations to estimate particle-scale material properties

  • A demonstration of the use of this framework to calibrate material properties of agricultural grains (couscous and barley) for a simulator, which is then evaluated against new scenarios in simulation and experiment

  • A robotic demonstration that uses the proposed framework to reason about agricultural grains in order to complete a simple manipulation task

Ii Related Work

Ii-a Granular Materials in Physics

The mechanics of granular materials is an active area of research within the physics community. Recent advancements present new constitutive relations of granular media that relate flow behavior to particle-scale properties like grain size [19], geometry [30], and surface friction [21]. Throughout this body of work, DEM simulators are consistently used to evaluate and extend theoretical explanations [12, 42, 21, 18, 30]. This emphasis on simulation encourages the development of frameworks that can infer difficult-to-measure material parameters (e.g., frictional coefficients of a grain) from macroscale behavior of collections of grains (e.g., pile shape). This work proposes one such framework, designed primarily for robotic applications due to its speed, efficiency, and minimal hardware requirements.

Ii-B Granular Materials in Robotics

Within robotics, the mechanics of granular materials has been investigated to design jamming grippers and manipulators [5, 41, 25, 9], as well as mobile robots that can traverse granular terrain [44, 23, 29, 22, 43, 32]. Recent work has explored the perception of granular materials [10, 7, 40] and their manipulation in scooping [10] and pouring [39] tasks. Such work has largely focused on the classification of granular materials rather than estimation of material parameters; furthermore, it has often been highly task-specific. However, granular parameter identification is important to enable physical reasoning of how different grains may behave in unseen scenarios while dynamically interacting with robots. Due to the complexity of granular physics, researchers in robotics have leveraged simulation to model robot-granular interactions [44, 43, 29, 32].

Ii-C Simulation of Granular Materials

The state-of-the-art numerical technique to simulate granular materials is the Discrete Element Method (DEM), where each grain is generally represented as an independent particle that interacts with others only through direct contact. Other granular simulation techniques include the Material Point Method (MPM) [16, 15], which constructs a background mesh to compute forces, and Position Based Dynamics (PBD) [31, 27, 28], which directly manipulates particle positions to enforce attachments and resolve collisions, typically trading off accuracy for speed. Although accurate, DEM often uses explicit time-integration, which can require impractically small time steps for stiff systems. In this work we model contact as a nonlinear complementarity problem (NCP) which is solved implicitly using a GPU-optimized non-smooth Newton solver in the NVIDIA Isaac Simulator platform [2, 26]. This enables both high accuracy and rapid simulation of

2000 spherical grains. Furthermore, the platform allows direct input of particle-scale parameters such as geometry and material properties, as well as a programmatic interface in which new testing scenarios can be rapidly created. The platform has recently been applied to reinforcement learning applications for non-DEM simulations

[6, 24].

Ii-D Automatic Calibration of Granular Material Simulation

Most techniques used to calibrate the material parameters for DEM simulation require tedious experiments and manual tuning until the desired experimental behavior is captured [11]. Thus, recent work aims to automate this process. For instance, in [14]

, a genetic algorithm and direct optimization are used to determine the frictional coefficients of sand grains by measuring discharge time from an hourglass. In


, an artificial neural network is used to infer multiple material parameters of iron ore sinter; the network greatly increases the speed of inference at the expense of training time. Recently, likelihood-free inference techniques, such as Approximate Bayesian Computation (ABC) methods, have been explored to increase the interpretability of results, as well as to enhance robustness to non-unique mappings between material properties and macroscale behavior. Specifically, in

[8], an iterative Bayesian filtering framework was used to infer the coefficients of friction of glass beads based on 100 measurements of porosity and stress per iteration. In the present work, we leverage a likelihood-free Bayesian inference framework called BayesSim [37]. The framework allows the complexities of the simulator to be fully abstracted from the inference process and requires only a single training phase followed by a prediction step. This is in contrast to multiple prediction steps necessary in typical Bayesian filtering. Because BayesSim is an -free method [34] (i.e., it does not rely on sampling with acceptance region defined by ), it is generally more sample-efficient and produces more accurate parametric approximations of the posteriors than ABC methods [34]. Furthermore, the preceding studies measured macroscopic behavior using highly-specific experimental setups and advanced laboratory equipment, which may be prohibitively expensive for many robotics applications. In contrast, the present work measures macroscopic behavior (e.g., granular formations from pouring grains through a funnel) using single images taken with inexpensive, off-the-shelf depth cameras.

Iii Problem Statement

Iii-a Simplifying Assumptions

This section defines the key simplifications made in the simulator. First, we consider grains with a length scale of at least 1 mm, thus excluding powders. Given that the average size and mass of a grain can be measured, they are assumed as known inputs to the simulator. While natural grains will represent a distribution of sizes and masses, we make the simplification that all grains are of the same size and mass. We further assume that the grains are perfectly rigid, cohesionless, and experience negligible drag, so that the only sources of energy dissipation are friction and collisions with each other and the ground. The forces experienced during these interactions, depicted in Figure 2, are functions of the coefficients of restitution , sliding friction , and rolling friction for the grain-to-grain or grain-to-ground interactions. To simplify parameter search, we assume that the grain-to-grain and grain-to-ground coefficients are equivalent. Together, these three parameters enable the simulator to generate a diversity of realistic granular behavior. Using a subset of restricts the simulator’s expressiveness and may likely cause the real and simulated distributions of granular formations to be completely distinct. Finally, we represent all grains as spheres. While it has been shown that aggregate behavior is dependent on grain geometry and size, this coarse approximation makes collisions easier and faster to compute, and thus practical for robotic applications. Further, although the parameters are bounded by real-world values, it is not the goal of this paper to estimate physically-exact friction and restitution coefficients. Rather, the calibration framework aims to accurately simulate macroscopic behavior of the granular material. For non-spherical grains, larger or smaller valued material parameters may sufficiently compensate for a spherical approximation.

Iii-B Definitions


represent a vector of the grain simulation parameters, composed of the coefficients of sliding friction

, rolling friction , and restitution . Let represent the birds-eye-view depth image of the grains once they have reached steady-state. Variables and are used to distinguish real from simulated depth images, respectively, and , where is the simulation. Because of the high-dimensional space of depth images and

, we define the random variable

as a vector composed of statistics (such as average distance of grains from the centroid) to summarize the observed state . Let the likelihood function be defined as and the posterior be defined as . The likelihood function, defined implicitly by , is often intractable to evaluate or unavailable with such complex simulations, thus motivating the use of a likelihood-free inference approach.

Iii-C Objective

Given samples of and the corresponding parameters used to generate the simulations, we would like to approximate the posterior function such that, from a real observation , we can infer granular simulation parameters from the macroscopic behavior of the material, summarized by features . Because DEM simulations leverage various numerical solvers and are associated with complex systems of differential equations, we choose a simulator-agnostic inference framework that approximates the intractable posterior by sampling from the simulator. The likelihood-free Bayesian inference framework BayesSim [37], discussed in the following section, is used to approximate this posterior.

Iv Methods

The goal of this paper’s proposed framework is to estimate granular parameters from summary statistics using BayesSim, which approximates the posterior by sampling from a simulator . In this section, we discuss the BayesSim likelihood-free framework, simulator , physical experimental setup to observe , and extraction of summary statistics .

Fig. 2: The simulations capture sliding friction forces (), rolling friction forces (), and energy loss due to inelastic collisions for grain-grain and grain-ground interactions. These forces are parameterized by the coefficients of sliding friction (), rolling friction (), and restitution (), respectively.

Iv-a Granular Parameter Inference Using BayesSim

BayesSim directly learns a conditional density function, (where is the set of learnable parameters of ), by mapping summary statistics to a mixture of Gaussian components via random Fourier features [36]. The Mixture Density Random Fourier Feature (MDRFF) network learns from pairs of parameters and summary statistics, , generated by independently sampling from a predefined prior distribution, , and forward-simulating each sample, . For twice differentiability of , the Matérn kernel is used with , where increases with smoothness. The posterior is then parametrically approximated using the learned conditional density and a real observation . The main advantages of using BayesSim over classical optimization methods are that the complexities and inaccuracies of the DEM simulator and PDE solver are incorporated in the estimated posterior. The multimodality and uncertainty of the approximated posterior can be used to analyze the accuracy with respect to the parameters and robustness with respect to measurement noise. In addition, because BayesSim is an -free method, it improves on ABC methods by being more sample efficient and producing more accurate posterior approximations [34] with a limited number of samples. After the initial cost of running simulations and training the model, BayesSim can quickly produce posterior approximations, making it a compelling method for granular material inference problems in robotic applications.

Iv-B Simulation

Critical dimensions are made equivalent in the simulated and physical scenes. The simulated funnel, which was modeled in 3D based on the physical funnel, is positioned at an equivalent distance above the ground surface. The spherical grains are scaled to match the diameter of a single grain of Israeli couscous along its longest axis (4 mm). For simplicity, at the start of the simulation, the 2000 grains are each spaced at 8 mm intervals in a 10 10 20 grid, where the bottom grains are just below the top of the funnel. The coefficient of restitution is the ratio of a grain’s speed before and after a collision; this coefficient ranges from 0 to 1, corresponding to perfectly inelastic and elastic collisions, respectively. The ranges of the coefficients of sliding friction and rolling friction that create a noticeable diversity in granular formations are 0 to 1 and to , respectively. For this reason, rolling friction is sampled from a logarithmic distribution. A simulated depth camera is positioned at a birds-eye viewpoint, 29 cm from the surface ground, giving it a 58 cm 58 cm field of view. The solver chosen is NVIDIA Isaac’s Preconditioned Conjugate Residual (PCR) solver, and 20 inner iterations with 4 outer iterations are used in experiments [26]. The simulator runs with time step sec, 10 substeps, and a relaxation factor of 0.75 for stability and convergence.

Iv-C Physical Setup

For the physical experiments, the funnel is positioned with its bottom tip 12 cm above the ground surface. The surface is covered in thin, moderate-friction black velvet to limit the spread of the grains upon impact, as well as to simplify segmentation. The primary grain used in the experiments is Israeli couscous, which has the approximate geometry of a rounded cylinder; barley, which has a similar density but is slightly larger and more ellipsoidal, is also used in the latter half of the experiments. An Intel RealSense D435 depth camera is positioned at 40 cm above the ground surface. This position was chosen to balance the trade-off between field-of-view and clarity of the depth image. Although the raw depth images are still noisy at 40 cm, eliminating the background using color segmentation results in well-defined point clouds of the granular formations. An ABB YuMi dual-arm robot is used for the robot experiments.

Iv-D Summary Statistic Extraction

A majority of DEM calibration work relies on extracting specific features (e.g., angle of repose) of a conical pile in order to infer a single parameter (e.g., ). However, initial experiments showed that pouring grains can result in piles or rings, which may be more informative of other material parameters. Moreover, for shallow formations, it may be difficult to measure angle of repose from a depth image. We take a more generalizable approach and propose using a vector of 16 summary statistics to serve as . To extract the summary statistics features from a granular formation, the depth image is first downsampled by factor . For , and for , to match resolutions. Because the are fairly noisy, a color threshold mask is used to segment pixels corresponding to the granular formations. These pixels are then reprojected into 3D using the camera intrinsics matrix. For the real depth image, a plane is fit to reprojected pixels corresponding to the ground surface. The orientation of the fitted ground plane informs the centering and re-orientation of the granular point cloud to account for camera tilt. Let represent the vector of radial distances of grain points in the horizontal plane from the centroid of the granular formation and be the vector of grain point heights. The summary statistics of a single depth image are chosen as , where , , , , and

are the mean, standard deviation, interquartile range, kurtosis, and distance correlation. A chi-distribution, chosen for its likeness in shape to a cross-section of a granular formation, is fit to

. Statistics , , and

denote degrees of freedom, shift, and scale, respectively, from a standardized chi probability density. These statistics were chosen to summarize the dispersion and statistical dependence of the grain locations. A preliminary ablation study suggested that subsets of

were sufficiently descriptive in capturing the granular distributions, unless the subsets contained less than three features or the features exclusively pertained to either or .

V Sensitivity Tests for System Design

Fig. 3: Sensitivity to observation error. An example is shown perturbed by different levels of noise. 50 different test images were generated by evenly sampling . The resulting approximated posteriors are plotted per true value of , where the highest peaks correspond to the brightest points. A perfect inference would result in a white line overlapping the red. While the model is fairly robust to the level of per-pixel Gaussian depth noise applied, accuracy is more affected by Gaussian blur. Posterior peaks become wider, and distributions can become multimodal with .

V-a Preliminary Tests for System Noise

When running the simulator to train BayesSim, there are two main sources of uncertainty. The first is purely algorithmic, due to the non-determinism of the simulator and numerical error. To measure algorithmic uncertainty, the simulator was run 1000 times with identical conditions except for sliding friction, where , evenly split. The measured summary statistics for each

had variances of at least two orders of magnitude smaller than their mean values. Second, an insufficient number of samples for training and interpolation uncertainty may contribute to poor performance. To test these, a training set was created by sampling

uniformly between 0.01 and 1. A 50-simulation testing set was generated by sampling between 0.01 and 1, independently of the training set. Supporting intuition, the accuracy of the friction coefficients inferred by BayesSim improved with more samples. For a training set of 500 samples, the error in the friction coefficient estimates is , while for 1000 samples, the error is reduced to .

V-B Sensitivity to Observation Error

A key assumption of our system is that the summary statistics of the real and simulated depth images will lie within the same distribution. However, the two might not align or overlap given drastically different observations (e.g., due to a low-quality experimental depth image). To test the effect of observation noise on inference, we simulate two types of camera noise on a 50-sample testing set, where is evenly sampled from 0.01 to 1 with all other parameters equivalent. An XY Gaussian blur with is convolved with the depth image, and pixel-wise Gaussian noise with is added to each depth value, with ranges chosen to mimic realistic noise magnitudes. These two perturbations are applied independently, and the resulting posteriors per true friction value are illustrated in Figure 3. As shown, inference of is fairly robust to pixel-wise noise, but is sensitive to high values of Gaussian blur.

V-C Propagation of Error

Given an observation , the approximated posterior encodes the uncertainty of the inferred parameters . Sensitivity tests were performed to investigate how this uncertainty propagates to errors in the summary statistics . Fixing all other parameters, , , and required errors of , respectively, to perturb the height and average radial distance to more than the width of a single grain. For real-world context, manufacturing applications may require that the diameter of a granular pile be predicted to within 5% of the true value. (For example, for a pile with base diameter equal to 50 times a grain diameter of 4 mm, an error of 1 cm may be allowable). With the proposed system, this specification is fulfilled by the optimal values, with allowable variation of that is within the performance achieved in the following experiments.

Vi Results

Vi-a Sim-to-sim Parameter Estimation

Fig. 4: Examples comparing simulated depth images of testing set and simulated depth images from running a forward simulation with inferred parameters. Table shows corresponding errors in inferred parameters, as well as L2 errors between summary statistics.

The first objective was to test the accuracy of the trained BayesSim model on simulated depth images. Given a simulated input image , we aimed to quantify how close the inferred parameters were to . Furthermore, to evaluate the propagation of error in , a forward simulation was run to obtain a new pile . The difference between the piles was evaluated via the L2-norm error between and , both standardized over the summary statistics of the training dataset. Note that, for some applications, evaluating particular subsets of the summary statistics might be more useful (e.g., height of the granular pile or the angle of repose). To enable generality, all 16 summary statistics defined in Section IV-D are used to calculate the L2-norm error. We avoid evaluating accuracy via direct comparisons of the downsampled depth images, as the computed error would be highly sensitive to misalignments of the scattered grains rather than reflect meaningful differences in spatial distributions. The 50 test grain formations were simulated by randomly sampling and uniformly and logarithmically from their given ranges. The average error of the inferred was , while that of and were and , respectively. As in Section V-C, these error bounds correspond to negligible changes in the resulting statistics. As shown in Figure 4, forward simulations of the inferred parameters for the 50 test formations generate visually equivalent shapes. To consider as a baseline for the following experiments, the average L2-norm error between the 50 pairs of was . Sim-to-sim comparisons also enable a principled approach to testing sensitivity to model error. For instance, the grains are spherical, when in reality most grains are not. We tested sensitivity to mismatched geometry by creating a testing set of cylindrical and cube-shaped grains, where the diameter and length of each grain were 4 mm, the diameter of the training set spheres. Each testing set consisted of 50 different formations which varied in , while and were fixed at 0.5 and , respectively. Parameters , , and were inferred from the test formations and forward simulated with spherical grains. Comparing (with cylindrical and cubic grains) with (using spherical grains) the L2-norm errors were and , respectively. Spherical grains can approximate formations made with cubes better than cylinders, possibly due to the asymmetry of cylinder geometry. Depending on the application, a spherical representation may be sufficient.

Vi-B Real-to-sim Parameter Estimation

The goal of this experiment was to test whether material parameters of real-world granular formations could be inferred to produce accurate simulations of macroscopic behavior. The initial test set was composed of real depth images of 2000 grains of Israeli couscous, poured through a funnel fixed at 12 cm above the ground surface. Ten different pours were executed, and for each, ten consecutive depth images were collected and averaged to mitigate temporal noise from the camera. For each of the averaged images , summary statistics were computed, and , , and were inferred. Forward simulation was executed for each parameter set, and the L2-error was computed between the summary statistics of the real and forward-simulated granular formations. From the ten parameter sets inferred, the five corresponding to the lowest errors were averaged. The inferred parameters for couscous were , , and , and the parameters for barley were , , and . Forward simulation with the estimated parameters captures similar macroscopic behavior as the five real-world pours at 12 cm, with an average L2-norm error of and for the couscous and barley formations, respectively. However, it is critical to test how well these parameters generalize to pours at different heights. Five more pours of each grain were performed, with each pour at a different height: 2, 4, 6, 8, 10 cm off the ground. For each of these pours, a forward simulation was run with the same grain parameters while shifting the funnel height. Figure 5 illustrates that the inferred parameters generalize reasonably well, with relatively low L2 errors. This suggests that it may be possible under this framework to reason about pour heights to create desired pile and ring shapes given a specific type of grain. We test this hypothesis in the following experiment.

Fig. 5: Comparisons of real depth images of piles poured at different heights and simulated depth images generated from running a forward simulation at the corresponding heights, using parameters inferred from a height of 12 cm. L2 errors are listed below the depth images.

Vi-C Robotic Demonstrations

The first demonstration evaluated whether the presented framework could enable an industrial robot to pour granular material into a desired shape, as may be necessary in a kitchen or factory. A new BayesSim model was trained on 1000 simulated piles. The simulator parameters , , and were set to their inferred values for couscous (, , and , respectively), and the funnel height was uniformly sampled between 1 and 13 cm. Two summary statistics were computed from the resulting granular formation (i.e., the 5th and 50th percentile of granular radial distance, which approximated the inner and outer radii of a ring). BayesSim was then evaluated on a testing set to infer the funnel height from these two summary statistics; the height was estimated within cm. For the demonstration, summary statistics were chosen to correspond to a desired pattern of two concentric rings with an inner pile. BayesSim inferred corresponding funnel heights of 27.1, 10.1, 1.5 cm. The robot was commanded to pour couscous from these heights, resulting in the desired pattern (Figure 6). The second demonstration evaluated how closely the presented framework could predict undesirable spilling during a granular pouring task. The calibrated simulator was used to simulate the pouring of couscous and barley into a cereal bowl, with the simulation relaxation constant increased from to to ensure accurate simulation of high-speed collisions. The number of grains that spilled out of the bowl was counted. The industrial robot was then used to precisely repeat the experiment in the real-world, pouring grains into a velvet-lined cereal bowl to match simulator conditions. As tabulated in Figure 6, in most cases, the calibrated simulator predicted the number of spilled grains across different pouring heights with surprising accuracy. (Please see the supplementary video.)

Fig. 6: Two granular manipulation tasks were tested. (Left): The height of the funnel is inferred to create a desired ring shape. Performance is demonstrated by creating a pattern of three concentric rings. (Right): The calibrated simulator is used to recreate a real-life scene of pouring grains into a bowl. At various heights, both the simulations of couscous and barley reasonably estimate the number of grains that leave the cereal bowl.

Vii Conclusions

In this paper, the material parameters of granular materials were inferred using a new framework combining likelihood-free Bayesian inference, efficient simulation, and simple experiments. The use of GPU-based simulation and off-the-shelf depth cameras may be particularly appealing for robotics applications. Simulation-to-simulation inference was highly accurate, and simulation-to-experiment inference trailed closely in performance. Robotics demos showed that the inferred parameters generalized well to different pouring heights, and furthermore, that a robot can effectively reason about granular material to pour desired granular formations and predict its behavior in dynamic scenarios. There are numerous exciting opportunities for future work in perceiving and reasoning about granular materials. To improve inference, observations can be amended to include temporal sequences of depth images, with additional summary statistics to capture transient dynamics. To enhance simulation, NVIDIA’s Isaac Simulator can be applied to model cohesion, allowing liquid-solid interactions (e.g., capillary bridging); in addition, simulated particle geometries could be represented as superquadrics [35], enabling more accurate approximation of complex real-world grain geometries. Finally, more elaborate robotic manipulation tasks can be explored, such as moving the end effector along a trajectory to create an asymmetric grain trail, or scooping and pouring granular material into and out of assorted containers.


The authors were supported in part by the National Science Foundation Graduate Research Fellowship. We thank Miles Macklin, Renato Gasoto, and Hammad Mazhar for their tireless assistance with the NVIDIA Isaac Simulator, as well as Matthew Matl, Karl Van Wyk, and Qian Wan for their insightful feedback and suggestions.


  • [1] Granular material by design. Accessed:2019-09-15.
  • [2] NVIDIA Isaac: Virtual simulator for robots. Accessed: 2019-09-11.
  • [3] L Anand and C Gu. Granular materials: constitutive equations and strain localization. Journal of the Mechanics and Physics of Solids, 48(8):1701–1733, 2000.
  • [4] L Benvenuti, C Kloss, and S Pirker. Identification of DEM simulation parameters by artificial neural networks and bulk experiments. Powder technology, 291:456–465, 2016.
  • [5] Eric Brown, Nicholas Rodenberg, John Amend, Annan Mozeika, Erik Steltz, Mitchell R Zakin, Hod Lipson, and Heinrich M Jaeger. Universal robotic gripper based on the jamming of granular material. Proceedings of the National Academy of Sciences, 107(44):18809–18814, 2010.
  • [6] Yevgen Chebotar, Ankur Handa, Viktor Makoviychuk, Miles Macklin, Jan Issac, Nathan Ratliff, and Dieter Fox. Closing the sim-to-real loop: Adapting simulation randomization with real world experience. In 2019 International Conference on Robotics and Automation (ICRA), pages 8973–8979. IEEE, 2019.
  • [7] Carolyn L Chen, Jeffrey O Snyder, and Peter J Ramadge. Learning to identify container contents through tactile vibration signatures. In 2016 IEEE International Conference on Simulation, Modeling, and Programming for Autonomous Robots (SIMPAR), pages 43–48. IEEE, 2016.
  • [8] Hongyang Cheng, Takayuki Shuku, Klaus Thoeni, Pamela Tempone, Stefan Luding, and Vanessa Magnanimo. An iterative Bayesian filtering framework for fast and automated calibration of DEM models. Computer Methods in Applied Mechanics and Engineering, 350:268–294, 2019.
  • [9] Matteo Cianchetti, Tommaso Ranzani, Giada Gerboni, Thrishantha Nanayakkara, Kaspar Althoefer, Prokar Dasgupta, and Arianna Menciassi. Soft robotics technologies to address shortcomings in today’s minimally invasive surgery: the STIFF-FLOP approach. Soft Robotics, 1(2):122–131, 2014.
  • [10] Samuel Clarke, Travers Rhodes, Christopher G Atkeson, and Oliver Kroemer. Learning audio feedback for estimating amount and flow of granular material. In Conference on Robot Learning, pages 529–550, 2018.
  • [11] CJ Coetzee. Calibration of the discrete element method. Powder Technology, 310:104–142, 2017.
  • [12] Peter A Cundall and Otto DL Strack. A discrete numerical model for granular assemblies. geotechnique, 29(1):47–65, 1979.
  • [13] Karola Dierichs, Tobias Schwinn, and Achim Menges. Robotic pouring of aggregate structures. In Rob| Arch 2012, pages 196–205. Springer, 2013.
  • [14] Huy Q Do, Alejandro M Aragón, and Dingena L Schott. Automated discrete element method calibration using genetic and optimization algorithms. In EPJ Web of Conferences, volume 140, page 15011. EDP Sciences, 2017.
  • [15] Sachith Dunatunga and Ken Kamrin. Continuum modelling and simulation of granular flows through their many phases. Journal of Fluid Mechanics, 779:483–513, 2015.
  • [16] Yuanming Hu, Yu Fang, Ziheng Ge, Ziyin Qu, Yixin Zhu, Andre Pradhana, and Chenfanfu Jiang. A moving least squares material point method with displacement discontinuity and two-way rigid body coupling. ACM Transactions on Graphics (TOG), 37(4):150, 2018.
  • [17] Heinrich M Jaeger and Sidney R Nagel. Physics of the granular state. Science, 255(5051):1523–1531, 1992.
  • [18] Pierre Jop, Yoël Forterre, and Olivier Pouliquen. A constitutive law for dense granular flows. Nature, 441(7094):727, 2006.
  • [19] Ken Kamrin. Quantitative rheological model for granular materials: The importance of particle size. Handbook of Materials Modeling: Applications: Current and Emerging Materials, pages 1–24, 2018.
  • [20] Ken Kamrin. Non-locality in granular flow: Phenomenology and modeling approaches. Frontiers in Physics, 7:116, 2019.
  • [21] Ken Kamrin and Georg Koval. Effect of particle surface friction on nonlocal constitutive behavior of flowing granular media. Computational Particle Mechanics, 1(2):169–176, 2014.
  • [22] Margaret A Knuth, JB Johnson, MA Hopkins, RJ Sullivan, and JM Moore. Discrete element modeling of a Mars Exploration Rover wheel in granular material. Journal of Terramechanics, 49(1):27–36, 2012.
  • [23] Chen Li, Tingnan Zhang, and Daniel I Goldman. A terradynamics of legged locomotion on granular media. science, 339(6126):1408–1412, 2013.
  • [24] Jacky Liang, Viktor Makoviychuk, Ankur Handa, Nuttapong Chentanez, Miles Macklin, and Dieter Fox. GPU-accelerated robotic simulation for distributed reinforcement learning. arXiv preprint arXiv:1810.05762, 2018.
  • [25] Andrea J Liu and Sidney R Nagel. Granular and jammed materials. Soft Matter, 6(13):2869–2870, 2010.
  • [26] Miles Macklin, Kenny Erleben, Matthias Müller, Nuttapong Chentanez, Stefan Jeschke, and Viktor Makoviychuk. Non-smooth Newton methods for deformable multi-body dynamics. ACM Transactions on Graphics (TOG), 38(5):1–20, 2019.
  • [27] Miles Macklin and Matthias Müller. Position based fluids. ACM Transactions on Graphics (TOG), 32(4):104, 2013.
  • [28] Miles Macklin, Matthias Müller, Nuttapong Chentanez, and Tae-Yong Kim. Unified particle physics for real-time applications. ACM Transactions on Graphics (TOG), 33(4):153, 2014.
  • [29] Ryan D Maladen, Yang Ding, Paul B Umbanhowar, and Daniel I Goldman. Undulatory swimming in sand: experimental and simulation studies of a robotic sandfish. The International Journal of Robotics Research, 30(7):793–805, 2011.
  • [30] Hammad Mazhar, Jonas Bollmann, Endrina Forti, Andreas Praeger, Tim Osswald, and Dan Negrut. Studying the effect of powder geometry on the selective laser sintering process. Society of Plastics Engineers (SPE) ANTEC, 2014.
  • [31] Matthias Müller, Bruno Heidelberger, Marcus Hennix, and John Ratcliff. Position based dynamics. Journal of Visual Communication and Image Representation, 18(2):109–118, 2007.
  • [32] Dan Negrut, Daniel Melanz, Hammad Mazhar, David Lamb, Paramsothy Jayakumar, and Michael Letherwood. Investigating through simulation the mobility of light tracked vehicles operating on discrete granular terrain. SAE International Journal of Passenger Cars-Mechanical Systems, 6(2013-01-1191):369–381, 2013.
  • [33] Sia Nemat-Nasser and Juhua Zhang. Constitutive relations for cohesionless frictional granular materials. International Journal of Plasticity, 18(4):531–547, 2002.
  • [34] George Papamakarios and Iain Murray. Fast -free inference of simulation models with Bayesian conditional density estimation. In Advances in Neural Information Processing Systems, pages 1028–1036, 2016.
  • [35] Alexander Podlozhnyuk, Stefan Pirker, and Christoph Kloss.

    Efficient implementation of superquadric particles in discrete element method within an open-source framework.

    Computational Particle Mechanics, 4(1):101–118, 2017.
  • [36] Ali Rahimi and Benjamin Recht. Random features for large-scale kernel machines. In Advances in Neural Information Processing Systems, pages 1177–1184, 2008.
  • [37] Fabio Ramos, Rafael Carvalhaes Possas, and Dieter Fox. BayesSim: adaptive domain randomization via probabilistic inference for robotics simulators. arXiv preprint arXiv:1906.01728, 2019.
  • [38] Patrick Richard, Mario Nicodemi, Renaud Delannay, Philippe Ribiere, and Daniel Bideau. Slow relaxation and compaction of granular systems. Nature materials, 4(2):121, 2005.
  • [39] Connor Schenck, Jonathan Tompson, Dieter Fox, and Sergey Levine. Learning robotic manipulation of granular media. arXiv preprint arXiv:1709.02833, 2017.
  • [40] Jivko Sinapov, Connor Schenck, and Alexander Stoytchev. Learning relational object categories using behavioral exploration and multimodal perception. In 2014 IEEE International Conference on Robotics and Automation (ICRA), pages 5691–5698. IEEE, 2014.
  • [41] Elliot Thompson-Bean, Oliver Steiner, and Andrew McDaid. A soft robotic exoskeleton utilizing granular jamming. In 2015 IEEE International Conference On Advanced Intelligent Mechatronics (AIM), pages 165–170. IEEE, 2015.
  • [42] Qiong Zhang and Ken Kamrin. Microscopic description of the granular fluidity field in nonlocal flow modeling. Physical review letters, 118(5):058001, 2017.
  • [43] Tingnan Zhang, Feifei Qian, Chen Li, Pierangelo Masarati, Aaron M Hoover, Paul Birkmeyer, Andrew Pullin, Ronald S Fearing, and Daniel I Goldman. Ground fluidization promotes rapid running of a lightweight robot. The International Journal of Robotics Research, 32(7):859–869, 2013.
  • [44] Yifan Zhu, Laith Abdulmajeid, and Kris Hauser. A data-driven approach for fast simulation of robot locomotion on granular media. In 2019 International Conference on Robotics and Automation (ICRA), pages 7653–7659. IEEE, 2019.