Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userContact model for elastically anisotropic bodies and efficient implementation into the discrete-element method
68
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce a contact law for the normal force generated between two contacting, elastically anisotropic bodies of arbitrary geometry. The only requirement is that their surfaces be smooth and frictionless. This anisotropic contact law is obtained from a simplification of the exact solution to the continuum elasticity problem and takes the familiar form of Hertz contact law, with the only difference being the orientation-dependence of the material-specific contact modulus. The contact law is remarkably accurate when compared with the exact solution, for a wide range of materials and surface geometries. We describe a computationally efficient implementation of the contact law into a discrete element method code, taking advantage of the precomputation of the contact modulus over all possible orientations. Finally, we showcase two application examples based on real materials where elastic anisotropy of the particles induces noticeable effects on macroscopic behavior. Notably, the second example demonstrates the ability to engineer tunable vibrational band gaps in a one-dimensional granular crystal by mere rotation of the constituent spherical particles.
rate research
Read More
A flexible fiber model based on the discrete element method (DEM) is presented and validated for the simulation of uniaxial compression of flexible fibers in a cylindrical container. It is found that the contact force models in the DEM simulations have a significant impact on compressive forces exerted on the fiber bed. Only when the geometry-dependent normal contact force model and the static friction model are employed, the simulation results are in good agreement with experimental results. Systematic simulation studies show that the compressive force initially increases and eventually saturates with an increase in the fiber-fiber friction coefficient, and the fiber-fiber contact forces follow a similar trend. The compressive force and lateral shear-to-normal stress ratio increase linearly with increasing fiber-wall friction coefficient. In uniaxial compression of frictional fibers, more static friction contacts occur than dynamic friction contacts with static friction becoming more predominant as the fiber-fiber friction coefficient increases.
A fully parallel version of the contact dynamics (CD) method is presented in this paper. For large enough systems, 100% efficiency has been demonstrated for up to 256 processors using a hierarchical domain decomposition with dynamic load balancing. The iterative scheme to calculate the contact forces is left domain-wise sequential, with data exchange after each iteration step, which ensures its stability. The number of additional iterations required for convergence by the partially parallel updates at the domain boundaries becomes negligible with increasing number of particles, which allows for an effective parallelization. Compared to the sequential implementation, we found no influence of the parallelization on simulation results.
We present a statistical model which is able to capture some interesting features exhibited in the Brazilian test. The model is based on breakable elements which break when the force experienced by the elements exceed their own load capacity. In this model when an element breaks, the capacity of the neighboring elements are decreased by a certain amount assuming weakening effect around the defected zone. We numerically investigate the stress-strain behavior, the strength of the system, how it scales with the system size and also its fluctuation for both uniformly and weibull distributed breaking threshold of the elements in the system. We find that the strength of the system approaches its asymptotic value $sigma_c=1/6$ and $sigma_c=5/18$ for uniformly and Weibull distributed breaking threshold of the elements respectively. We have also shown the damage profile right at the point when the stress-strain curve reaches at its maximum and then it is compared with our experimental observations.
We investigate theoretically and numerically the use of the Least-Squares Finite-element method (LSFEM) to approach data-assimilation problems for the steady-state, incompressible Navier-Stokes equations. Our LSFEM discretization is based on a stress-velocity-pressure (S-V-P) first-order formulation, using discrete counterparts of the Sobolev spaces $H({rm div}) times H^1 times L^2$ respectively. Resolution of the system is via minimization of a least-squares functional representing the magnitude of the residual of the equations. A simple and immediate approach to extend this solver to data-assimilation is to add a data-discrepancy term to the functional. Whereas most data-assimilation techniques require a large number of evaluations of the forward-simulations and are therefore very expensive, the approach proposed in this work uniquely has the same cost as a single forward run. However, the question arises: what is the statistical model implied by this choice? We answer this within the Bayesian framework, establishing the latent background covariance model and the likelihood. Further we demonstrate that - in the linear case - the method is equivalent to application of the Kalman filter, and derive the posterior covariance. We practically demonstrate the capabilities of our method on a backward-facing step case. Our LSFEM formulation (without data) is shown to have good approximation quality, even on relatively coarse meshes - in particular with respect to mass-conservation and reattachment location. Adding limited velocity measurements from experiment, we show that the method is able to correct for discretization error on very coarse meshes, as well as correct for the influence of unknown and uncertain boundary-conditions.
We develop a method to efficiently construct phase diagrams using machine learning. Uncertainty sampling (US) in active learning is utilized to intensively sample around phase boundaries. Here, we demonstrate constructions of three known experimental phase diagrams by the US approach. Compared with random sampling, the US approach decreases the number of sampling points to about 20%. In particular, the reduction rate is pronounced in more complicated phase diagrams. Furthermore, we show that using the US approach, undetected new phase can be rapidly found, and smaller number of initial sampling points are sufficient. Thus, we conclude that the US approach is useful to construct complicated phase diagrams from scratch and will be an essential tool in materials science.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
