No Arabic abstract
A discrete-element method (DEM) assembly of virtual particles is calibrated to approximate the behavior of a natural sand in undrained loading. The particles are octahedral, bumpy clusters of spheres that are compacted into assemblies of different densities. The contact model is a Jager generalization of the Hertz contact, which yields a small-strain shear modulus that is proportional to the square root of confining stress. Simulations made of triaxial extension and compression loading conditions and of simple shear produce behaviors that are similar to sand. Undrained cyclic shearing simulations are performed with nonuniform amplitudes of shearing pulses and with 24 irregular seismic shearing sequences. A methodology is proposed for quantifying the severities of such irregular shearing records, allowing the 24 sequences to be ranked in severity. The relative severities of the 24 seismic sequences show an anomalous dependence on sampling density. Four scalar measures are proposed for predicting the severity of a particular loading sequence. A stress-based scalar measure shows superior efficiency in predicting initial liquefaction and pore pressure rise.
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.
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.
Due to neutron irradiation, solid breeder blankets are subjected to complex thermo-mechanical conditions. Within one breeder unit, the ceramic breeder bed is composed of spherical-shaped lithium orthosilicate pebbles, and as a type of granular material, it exhibits strong coupling between temperature and stress fields. In this paper, we study these thermo-mechanical problems by developing a thermal discrete element method (Thermal-DEM). This proposed simulation tool models each individual ceramic pebble as one element and considers grain-scale thermo-mechanical interactions between elements. A small section of solid breeder pebble bed in HCPB is modelled using thousands of individual pebbles and subjected to volumetric heating profiles calculated from neutronics under ITER-relevant conditions. We consider heat transfer at the grain-scale between pebbles through both solid-to-solid contacts and the interstitial gas phase, and we calculate stresses arising from thermal expansion of pebbles. The overall effective conductivity of the bed depends on the resulting compressive stress state during the neutronic heating. The thermal-DEM method proposed in this study provides the access to the grain-scale information, which is beneficial for HCPB design and breeder material optimization, and a better understanding of overall thermo-mechanical responses of the breeder units under fusion-relevant conditions.
The New Horizons mission has returned stunning images of the bilobate Kuiper belt object (486958) Arrokoth. It is a contact binary, formed from two intact and relatively undisturbed predecessor objects joined by a narrow contact region. We use a version of pkdgrav, an N-body code that allows for soft-sphere collisions between particles, to model a variety of possible merger scenarios with the aim of constraining how Arrokoth may have evolved from two Kuiper belt objects into its current contact binary configuration. We find that the impact must have been quite slow (less than 5 m/s) and grazing (impact angles greater than 75 degrees) in order to leave intact lobes after the merger, in the case that both progenitor objects were rubble piles. A gentle contact between two bodies in a close synchronous orbit seems most plausible.
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.