Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hybrid material-point spheropolygon-element method for solid and granular material interaction
94
0
0.0
(
0
)
Added by
Amos Jiang
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Capturing the interaction between objects that have an extreme difference in Young s modulus or geometrical scale is a highly challenging topic for numerical simulation. One of the fundamental questions is how to build an accurate multi-scale method with optimal computational efficiency. In this work, we develop a material-point-spheropolygon discrete element method (MPM-SDEM). Our approach fully couples the material point method (MPM) and the spheropolygon discrete element method (SDEM) through the exchange of contact force information. It combines the advantage of MPM for accurately simulating elastoplastic continuum materials and the high efficiency of DEM for calculating the Newtonian dynamics of discrete near-rigid objects. The MPM-SDEM framework is demonstrated with an explicit time integration scheme. Its accuracy and efficiency are further analysed against the analytical and experimental data. Results demonstrate this method could accurately capture the contact force and momentum exchange between materials while maintaining favourable computational stability and efficiency. Our framework exhibits great potential in the analysis of multi-scale, multi-physics phenomena
rate research
Read More
A monolithic coupling between the material point method (MPM) and the finite element method (FEM) is presented. The MPM formulation described is implicit, and the exchange of information between particles and background grid is minimized. The reduced information transfer from the particles to the grid improves the stability of the method. Once the residual is assembled, the system matrix is obtained by means of automatic differentiation. In such a way, no explicit computation is required and the implementation is considerably simplified. When MPM is coupled with FEM, the MPM background grid is attached to the FEM body and the coupling is monolithic. With this strategy, no MPM particle can penetrate a FEM element, and the need for computationally expensive contact search algorithms used by existing coupling procedures is eliminated. The coupled system can be assembled with a single assembly procedure carried out element by element in a FEM fashion. Numerical results are reported to display the performances and advantages of the methods here discussed.
Mixtures of fluids and granular sediments play an important role in many industrial, geotechnical, and aerospace engineering problems, from waste management and transportation (liquid--sediment mixtures) to dust kick-up below helicopter rotors (gas--sediment mixtures). These mixed flows often involve bulk motion of hundreds of billions of individual sediment particles and can contain both highly turbulent regions and static, non-flowing regions. This breadth of phenomena necessitates the use of continuum simulation methods, such as the material point method (MPM), which can accurately capture these large deformations while also tracking the Lagrangian features of the flow (e.g. the granular surface, elastic stress, etc.). Recent works using two-phase MPM frameworks to simulate these mixtures have shown substantial promise; however, these approaches are hindered by the numerical limitations of MPM when simulating pure fluids. In addition to the well-known particle ringing instability and difficulty defining inflow/outflow boundary conditions, MPM has a tendency to accumulate quadrature errors as materials deform, increasing the rate of overall error growth as simulations progress. In this work, we present an improved, two-phase continuum simulation framework that uses the finite volume method (FVM) to solve the fluid phase equations of motion and MPM to solve the solid phase equations of motion, substantially reducing the effect of these errors and providing better accuracy and stability for long-duration simulations of these mixtures.
In order to accelerate implementation of hyperelastic materials for finite element analysis, we developed an automatic numerical algorithm that only requires the strain energy function. This saves the effort on analytical derivation and coding of stress and tangent modulus, which is time-consuming and prone to human errors. Using the one-sided Newton difference quotients, the proposed algorithm first perturbs deformation gradients and calculate the difference on strain energy to approximate stress. Then, we perturb again to get difference in stress to approximate tangent modulus. Accuracy of the approximations were evaluated across the perturbation parameter space, where we find the optimal amount of perturbation being $10^{-6}$ to obtain stress and $10^{-4}$ to obtain tangent modulus. Single element verification in ABAQUS with Neo-Hookean material resulted in a small stress error of only $7times10^{-5}$ on average across uniaxial compression and tension, biaxial tension and simple shear situations. A full 3D model with Holzapfel anisotropic material for artery inflation generated a small relative error of $4times10^{-6}$ for inflated radius at $25 kPa$ pressure. Results of the verification tests suggest that the proposed numerical method has good accuracy and convergence performance, therefore a good material implementation algorithm in small scale models and a useful debugging tool for large scale models.
We investigate the bulldozing motion of a granular sandpile driven forwards by a vertical plate. The problem is set up in the laboratory by emplacing the pile on a table rotating underneath a stationary plate; the continual circulation of the bulldozed material allows the dynamics to be explored over relatively long times, and the variation of the velocity with radius permits one to explore the dependence on bulldozing speed within a single experiment. We measure the time-dependent surface shape of the dune for a range of rotation rates, initial volumes and radial positions, for four granular materials, ranging from glass spheres to irregularly shaped sand. The evolution of the dune can be separated into two phases: a rapid initial adjustment to a state of quasi-steady avalanching perpendicular to the blade, followed by a much slower phase of lateral spreading and radial migration. The quasi-steady avalanching sets up a well-defined perpendicular profile with a nearly constant slope. This profile can be scaled by the depth against the bulldozer to collapse data from different times, radial positions and experiments onto common master curves that are characteristic of the granular material and depend on the local Froude number. The lateral profile of the dune along the face of the bulldozer varies more gradually with radial position, and evolves by slow lateral spreading. The spreading is asymmetrical, with the inward progress of the dune eventually arrested and its bulk migrating to larger radii. A one-dimensional depth-averaged model recovers the nearly linear perpendicular profile of the dune, but does not capture the finer nonlinear details of the master curves. A two-dimensional version of the model leads to an advection-diffusion equation that reproduces the lateral spreading and radial migration.
This paper studies a model of two-phase flow with an immersed material viscous interface and a finite element method for numerical solution of the resulting system of PDEs. The interaction between the bulk and surface media is characterized by no-penetration and slip with friction interface conditions. The system is shown to be dissipative and a model stationary problem is proved to be well-posed. The finite element method applied in this paper belongs to a family of unfitted discretizations. The performance of the method when model and discretization parameters vary is assessed. Moreover, an iterative procedure based on the splitting of the system into bulk and surface problems is introduced and studied numerically.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
