Contact force models for non-spherical particles with different surface properties : a review

This paper reviews the state-of-the-art contact force models for non-spherical particles, which describe the relationship between the contact force and the geometrical, material, and mechanical properties of the contacting particles. The review aims to select better contact force models to improve the current simplified contact force models used in discrete element method (DEM) simulations. First, the contact force models based on the extension of the classical Hertz theory are reviewed, in which a recent accurate geometrical contact force model is highlighted. Secondly, the research on the effects of different variables such as elastoplasticity, viscoelasticity, adhesion and surface roughness on contact force are reviewed respectively and then incorporated into the accurate geometrical contact force model. Thirdly, tangential force models for non-spherical particles in contact under various loading regimes are reviewed as well. Based on the review, a full set of improved contact force models for DEM is recommended. These contact force models can more accurately predict the contact force and contact area for non-spherical particles, which can be beneficial to the DEM simulation in emerging areas, such as nanoparticles and additive manufacturing

Discrete element modelling of the dynamic behaviour of non-spherical particulate materials

PhD ThesisA numerical model based on the discrete element (DE) method, for modelling the flow of irregularly shaped, smooth-surfaced particles in a 3-D system is presented. An existing DE program for modelling the contact between spherical particles in periodic space (without real walls or boundaries) was modified to model non-spherical particles in a system with containing walls. The new model was validated against analytical calculations of single particle movements and also experimentally against data from physical experiments using synthetic non-spherical particles at both a particle and bulk scale. It was then used to study the effect of particle shape on the flow behaviour of assemblies of particles with various aspect ratios discharging from a flat-bottomed hopper. The particles were modelled using the Multi-Sphere Method (MSM) which is based on the CSG (Constructive Solid Geometry) technique for construction of complex solids by combining primitive shapes. In this method particle geometry is approximated using overlapping spheres of arbitrary diameter which are fixed in position relative to each other. The contact mechanics and contact detection method are the same as those used for spheres, except that translation and rotation of element spheres are calculated with respect to the motion of the whole particle....Numerical simulations of packing and flow of particles from a flat-bottomed hopper with a range of aspect ratios were performed to investigate the effect of particle shape on packing and flow behaviour of a particulate assembly. It was found that the particle shape influenced both bed structure and flow characteristics such as flow pattern, shear band strength and the occurrence of bridging. The flow of the bed of spherical particles was smoother than the flow of beds of elongated particles in which flow was fluctuating and there was more resistance to shear.Ministry of Culture and Higher Education of IRAN: University of Mashhad

Compliant contact force models in multibody dynamics : evolution of the Hertz contact theory

Over the last decades, several compliant contact force models have been proposed. However, no complete and systematic comparison has been done on these models, which provides information on their range of application and accuracy for use in different contact scenarios. Thus, the selection of an appropriate model for a given contact problem is still an important and challenging issue to be addressed. The Hertzian contact theory remains the foundation for almost all of the available force models, but by itself, it is not appropriate for most impacts in practice, due to the amount of energy dissipated during the impact. A good number of contact force models have been offered that augment the Hertzian law with a damping term to accommodate the energy loss during the impact process for small or moderate impact velocities. In this work, the main issues associated with the most common compliant contact force models of this type are analyzed. Results in terms of the dynamic simulations of multibody systems are presented, which allow for the comparison of the similarities and differences among the models considered.Fundação para a Ciência e a Tecnologia (FCT) - DACHOR - Multibody Dynamics and Control of Hybrid Active Orthoses (MIT-Pt/BSHHMS/0042/2008), BIOJOINTS - Development of advanced biological joint models for human locomotion biomechanics (PTDC/EMEPME/ 099764/2008), SFRH/BD/40164/2007, SFRH/BD/64477/200

Fast Monte Carlo Simulations for Quality Assurance in Radiation Therapy

Monte Carlo (MC) simulation is generally considered to be the most accurate method for dose calculation in radiation therapy. However, it suffers from the low simulation efficiency (hours to days) and complex configuration, which impede its applications in clinical studies. The recent rise of MRI-guided radiation platform (e.g. ViewRay’s MRIdian system) brings urgent need of fast MC algorithms because the introduced strong magnetic field may cause big errors to other algorithms. My dissertation focuses on resolving the conflict between accuracy and efficiency of MC simulations through 4 different approaches: (1) GPU parallel computation, (2) Transport mechanism simplification, (3) Variance reduction, (4) DVH constraint. Accordingly, we took several steps to thoroughly study the performance and accuracy influence of these methods. As a result, three Monte Carlo simulation packages named gPENELOPE, gDPMvr and gDVH were developed for subtle balance between performance and accuracy in different application scenarios. For example, the most accurate gPENELOPE is usually used as golden standard for radiation meter model, while the fastest gDVH is usually used for quick in-patient dose calculation, which significantly reduces the calculation time from 5 hours to 1.2 minutes (250 times faster) with only 1% error introduced. In addition, a cross-platform GUI integrating simulation kernels and 3D visualization was developed to make the toolkit more user-friendly. After the fast MC infrastructure was established, we successfully applied it to four radiotherapy scenarios: (1) Validate the vender provided Co60 radiation head model by comparing the dose calculated by gPENELOPE to experiment data; (2) Quantitatively study the effect of magnetic field to dose distribution and proposed a strategy to improve treatment planning efficiency; (3) Evaluate the accuracy of the build-in MC algorithm of MRIdian’s treatment planning system. (4) Perform quick quality assurance (QA) for the “online adaptive radiation therapy” that doesn’t permit enough time to perform experiment QA. Many other time-sensitive applications (e.g. motional dose accumulation) will also benefit a lot from our fast MC infrastructure