scholarly journals DEM simulation of anisotropic granular materials: elastic and inelastic behavior

2020 ◽  
Vol 22 (4) ◽  
Author(s):  
Giuseppina Recchia ◽  
Vanessa Magnanimo ◽  
Hongyang Cheng ◽  
Luigi La Ragione
Keyword(s):  
Small Strain ◽  
Contact Forces ◽  
Inelastic Behavior ◽  
Discrete Elements ◽  
Normal Contact ◽  
Dem Simulation ◽  
Linear Elastic ◽  
Loading And Unloading ◽  
Incremental Strain ◽  
Incremental Loading

AbstractIn this work, Discrete Elements Method simulations are carried out to investigate the effective stiffness of an assembly of frictional, elastic spheres under anisotropic loading. Strain probes, following both forward and backward paths, are performed at several anisotropic levels and the corresponding stress is measured. For very small strain perturbations, we retrieve the linear elastic regime where the same response is measured when incremental loading and unloading are applied. Differently, for a greater magnitude of the incremental strain a different stress is measured, depending on the direction of the perturbation. In the case of unloading probes, the behavior stays elastic until non-linearity is reached.Under forward perturbations, the aggregate shows an intermediate inelastic stiffness, in which the main contribution comes from the normal contact forces. That is, when forward incremental probes are applied the behavior of anisotropic aggregates is an incremental frictionless behavior. In this regime we show that contacts roll or slide so the incremental tangential contact forces are zero. Graphical Abstract

A Nonlinear Rail Vehicle Dynamics Computer Program SAMS/Rail: Part 1—Theory and Formulations

2009 ◽  
Author(s):  
Khaled E. Zaazaa ◽  
Brian Whitten ◽  
Brian Marquis ◽  
Erik Curtis ◽  
Magdy El-Sibaie ◽  
...  
Keyword(s):  
High Speed ◽  
Dynamic Performance ◽  
Three Dimensional ◽  
Contact Element ◽  
Contact Forces ◽  
Simulation Code ◽  
Vehicle Performance ◽  
Normal Contact ◽  
Advantages And Disadvantages ◽  
High Speed Trains

Accurate prediction of railroad vehicle performance requires detailed formulations of wheel-rail contact models. In the past, most dynamic simulation tools used an offline wheel-rail contact element based on look-up tables that are used by the main simulation solver. Nowadays, the use of an online nonlinear three-dimensional wheel-rail contact element is necessary in order to accurately predict the dynamic performance of high speed trains. Recently, the Federal Railroad Administration, Office of Research and Development has sponsored a project to develop a general multibody simulation code that uses an online nonlinear three-dimensional wheel-rail contact element to predict the contact forces between wheel and rail. In this paper, several nonlinear wheel-rail contact formulations are presented, each using the online three-dimensional approach. The methods presented are divided into two contact approaches. In the first Constraint Approach, the wheel is assumed to remain in contact with the rail. In this approach, the normal contact forces are determined by using the technique of Lagrange multipliers. In the second Elastic Approach, wheel/rail separation and penetration are allowed, and the normal contact forces are determined by using Hertz’s Theory. The advantages and disadvantages of each method are presented in this paper. In addition, this paper discusses future developments and improvements for the multibody system code. Some of these improvements are currently being implemented by the University of Illinois at Chicago (UIC). In the accompanying “Part 2” and “Part 3” to this paper, numerical examples are presented in order to demonstrate the results obtained from this research.

Estimation of Graphite Dust Production in a Pebble-Bed Type Nuclear Reflector

2014 ◽  
Author(s):  
Ji-Ho Kang ◽  
Eung Seon Kim ◽  
Seungyon Cho
Keyword(s):  
Adhesive Wear ◽  
Cell Model ◽  
Estimation Method ◽  
Unit Cell ◽  
Contact Forces ◽  
Dust Explosion ◽  
Dust Production ◽  
Pebble Bed ◽  
Normal Contact ◽  
Contact Points

In this study, an estimation method of graphite dust production in the pebble-bed type reflector region of Korean HCSB (Helium-Cooled Solid Breeder) TBM (Test Blanket Module) in the ITER (International Thermonuclear Experimental Reactor) project using FEM (Finite Element Method) was proposed and the amount of dust production was calculated. A unit-cell model of uniformly arranged pebbles was defined with appropriate thermal and mechanical loadings. A commercial FEM program, Abaqus V6.10 was used to model and solve the stress field under multiple contact constraints between pebbles in the unit-cell. Resulting normal contact forces and slip distances on contact points were applied into the Archard adhesive wear equation to calculate the amount of graphite dust. The friction effect on contact points was investigated. The calculation result showed that the amount of graphite dust production was estimated to 2.22∼3.67e−4 g/m3 which was almost linearly proportional to the friction coefficient. The analysis results will be used as the basis data for the consecutive study of dust explosion.

scholarly journals NUMERICAL MODELLING OF THE SOIL BEHAVIOUR BY USING NEWLY DEVELOPED ADVANCED MATERIAL MODEL

2017 ◽  
Vol 57 (1) ◽  
pp. 58-70 ◽  
Author(s):  
Jan Veselý
Keyword(s):  
Numerical Modelling ◽  
Plastic Material ◽  
Material Model ◽  
Small Strain ◽  
Theoretical Background ◽  
Linear Elastic ◽  
In Situ Data ◽  
Modified Cam Clay ◽  
Soil Behaviour

This paper describes a theoretical background, implementation and validation of the newly developed Jardine plastic hardening-softening model (JPHS model), which can be used for numerical modelling of the soils behaviour. Although the JPHS model is based on the elasto-plastic theory, like the Mohr-Coulomb model that is widely used in geotechnics, it contains some improvements, which removes the main disadvantages of the MC model. The presented model is coupled with an isotopically hardening and softening law, non-linear elastic stress-strain law, non-associated elasto-plastic material description and a cap yield surface. The validation of the model is done by comparing the numerical results with real measured data from the laboratory tests and by testing of the model on the real project of the tunnel excavation. The 3D numerical analysis is performed and the comparison between the JPHS, Mohr-Coulomb, Modified Cam-Clay, Hardening small strain model and monitoring in-situ data is done.

Mesodamage Evolution in Polycrystals

2005 ◽  
Vol 127 (2) ◽  
pp. 214-221 ◽  
Author(s):  
M. Chadli ◽  
A. Abdul-Latif
Keyword(s):  
Fatigue Life ◽  
Micromechanical Model ◽  
Small Strain ◽  
Cyclic Stress ◽  
Damage Variable ◽  
Thermodynamic Force ◽  
Inelastic Behavior ◽  
Loading Paths ◽  
Localized Plastic Deformation ◽  
Phase Angles

A micromechanical model of damaged elasto-inelastic behavior is proposed to predict the plastic fatigue life for fcc metallic polycrystals under multiaxial loading paths. This model is expressed in the time-dependent plasticity for a small strain assumption. In order to generalize and then to increase the model applicability (with respect to other works of the author) in describing the cyclic stress-strain evolution during plastic fatigue, it is therefore assumed that a damage variable initiates and then evolves at the grain level where the phenomenon of the localized plastic deformation occurs. The associated thermodynamic force of the damage variable is determined as a total granular energy (elastic and inelastic). The transition of the elastic strain from the single to the polycrystal, which is classically performed by averaging procedures in this type of modeling, is modified due to the coupling of such a strain with damage. The developed model is tested under different multiaxial cyclic loading situations (tension-compression and tension-torsion with different out-of-phase angles). The effects the loading paths and the grains aggregate type on the fatigue life are appropriately investigated. It is demonstrated that the model can correctly describe the overall and local damaged behavior of polycrystals.

Effects of Detwinning on the Inelasticity of AZ31B Sheets During Cyclic Loading and Unloading

2018 ◽  
Vol 10 (09) ◽  
pp. 1850095 ◽  
Author(s):  
H. Wang ◽  
D. Tang ◽  
D. Y. Li ◽  
Y. H. Peng ◽  
P. D. Wu
Keyword(s):  
Cyclic Loading ◽  
Magnesium Alloys ◽  
Plasticity Model ◽  
Inelastic Behavior ◽  
Deformation History ◽  
Cyclic Loading And Unloading ◽  
Polycrystal Plasticity ◽  
Self Consistent ◽  
Loading And Unloading ◽  
Previous Loading

Magnesium alloys exhibit significant inelastic behavior during unloading, especially when twinning and detwinning are involved. It is commonly accepted that noteworthy inelastic behavior will be observed during unloading if twinning occurs during previous loading. However, this phenomenon is not always observed for Mg sheets with strong rolled texture. Therefore, the inelasticity of AZ31B rolled sheets with different rolled textures during cyclic loading-unloading are investigated by elastic viscoplastic self-consistent polycrystal plasticity model. The incorporation of the twinning and detwinning model enables the treatment of detwinning, which plays an important role for inelastic behavior during unloading. The effects of texture, deformation history, and especially twinning and detwinning on the inelastic behaviors are carefully investigated and found to be remarkable. The simulated results are in agreement with the available experimental observations, which reveals that the inelastic behavior for strongly rolled sheets is very different than the extruded bars.

Analysis of Vibration Characteristics for Last Stage Blade With Friction Contact Surfaces of Steam Turbine

2011 ◽  
Author(s):  
Yutaka Yamashita ◽  
Koki Shiohata ◽  
Takeshi Kudo
Keyword(s):  
Dynamic Response ◽  
Steam Turbine ◽  
Contact Surface ◽  
Contact Forces ◽  
Analysis Method ◽  
Friction Damping ◽  
Normal Contact ◽  
Damping Characteristics ◽  
Last Stage ◽  
Contact Surfaces

Friction damping devices such as under platform dampers are installed for modern turbine blades to suppress dynamic vibrations of the blades. In order to secure the reliability of the blades, it is important to predict the dynamic response and friction damping characteristics accurately. In this present paper, the dynamic response and friction damping characteristics of a last stage blade (LSB) of a steam turbine with contact surfaces at the cover, tie-boss and blade root was investigated. Especially, it is focused on the effect of the non-uniform normal contact forces at the contact surface. To investigate the effect of non-uniform normal contact forces, an analysis method was developed. Analysis model of the LSB with contact surfaces was discretized by finite elements. Tangential forces at the contact surfaces were modeled by multi-DOF macro-slip modeling. The non-linear frequency responses of the LSB were obtained by using the harmonic balance method. Using this analysis method, the relationship between the contact surface behavior and the dynamic response was studied.

Foundation for Virtual Prototyping of Mechanical Power Management Functions in Actuators

2018 ◽  
Author(s):  
Jean-Charles Mare ◽  
Silvio Akitani
Keyword(s):  
Power Management ◽  
Virtual Prototyping ◽  
Contact Forces ◽  
Mechanical Power ◽  
System Level ◽  
Special Focus ◽  
Force Model ◽  
Electromechanical Actuators ◽  
Normal Contact ◽  
Management Functions

Beside the main functions related to the control and transformation of power, safety-critical electromechanical actuators require many additional functions for power routing, protection and limitation. In practice, these functions are implemented mechanically because their realization at motor drive level is not acceptable for performance and reliability reasons. Contact forces play a major role in these mechanical devices (e.g. endstop, lock, brake, torque limiter, etc.), being either functional to serve the need, or parasitic due to their alteration of performance. The virtual prototyping of such mechanical power management functions therefore requires normal and tangent forces to be modelled with the right level of realism and reduced complexity. This communication provides some proposals to be used as foundation for the system-level modelling and simulation of these types of mechanical power elements that can be found in electromechanical actuators. Special focus is given to the model architecting, decomposition and block-diagram implementation, using the example of normal contact forces. The illustrative example concerns an integrated landing gear extension/retraction electromechanical actuator which embeds free-fall and autolock features. It shows how a well implemented single model (e.g. generic normal contact force model) combined with a right model decomposition can meet various modelling needs (e.g. droppable end-stop, lock and shearable axial stop). The proposed models are made compatible for integration in a 2x1D mechanical model architecture (axial and rotational motion) developed by the authors in previous reported work.

Shakedown in Frictional Contact Problems

2007 ◽  
Author(s):  
J. R. Barber ◽  
A. Klarbring ◽  
M. Ciavarella
Keyword(s):  
Frictional Contact ◽  
Elastic System ◽  
Contact Problems ◽  
Sufficient Condition ◽  
Normal Contact ◽  
Linear Elastic ◽  
Periodic Loading ◽  
Continuum Formulation ◽  
Complete Contact ◽  
The Continuum

If a linear elastic system with frictional interfaces is subjected to periodic loading, any slip which occurs generally reduces the tendency to slip during subsequent cycles and in some circumstances the system ‘shakes down’ to a state without slip. It has often been conjectured that a frictional Melan’s theorem should apply to this problem — i.e. that the existence of a state of residual stress sufficient to prevent further slip is a sufficient condition for the system to shake down. Here we discuss recent proofs that this is indeed the case for ‘complete’ contact problems if there is no coupling between relative tangential displacements at the interface and the corresponding normal contact tractions. By contrast, when coupling is present, the theorem applies only for a few special two-dimensional discrete cases. Counter-examples can be generated for all other cases. These results apply both in the discrete and the continuum formulation.

Unilateral Impact and Contact of Elastic Structures Using Lagrangian Multipliers

2011 ◽  
Author(s):  
Sebastian Tatzko
Keyword(s):  
Equations Of Motion ◽  
Structural Equations ◽  
Contact Forces ◽  
Integration Scheme ◽  
Lagrangian Multiplier ◽  
Elastic Structures ◽  
Lagrangian Multipliers ◽  
Linear Elastic ◽  
Time Stepping ◽  
Numerical Effort

This paper deals with linear elastic structures exposed to impact and contact phenomena. Within a time stepping integration scheme contact forces are computed with a Lagrangian multiplier approach. The main focus is turned on a simplified solving method of the linear complementarity problem for the frictionless contact. Numerical effort is reduced by applying a Craig-Bampton transformation to the structural equations of motion.

Effect of Housing Support on Bearing Dynamics

2015 ◽  
Vol 138 (1) ◽  
Author(s):  
Lijun Cao ◽  
Matthew D. Brouwer ◽  
Farshid Sadeghi ◽  
Lars-Erik Stacke
Keyword(s):  
Finite Element ◽  
Elastic Material ◽  
Contact Forces ◽  
Support Structure ◽  
Linear Elastic Material ◽  
Outer Race ◽  
Linear Elastic ◽  
Elastomeric Material ◽  
Bearing Dynamics ◽  
Element Method

The objective of this investigation was to determine the effect of housing support on bearing performance and dynamics. In order to achieve the objective, an existing dynamic bearing model (DBM) was coupled with flexible housing model to include the effect of support structure on bearing dynamics and performance. The DBM is based on the discrete element method, in which the bearing components are assumed to be rigid. To achieve the coupling, a novel algorithm was developed to detect contact conditions between the housing support and bearing outer race and then calculate contact forces based on the penalty method. It should be noted that although commercial finite element (FE) software such as abaqus is available to model flexible housings, combining these codes with a bearing model is quite difficult since the data transfer between the two model packages is time-consuming. So, a three-dimensional (3D) explicit finite element method (EFEM) was developed to model the bearing support structure for both linear elastic and nonlinear inelastic elastomeric materials. The constitutive relationship for elastomeric material is based on an eight chain model, which captures hyperelastic behavior of rubber for large strains. The viscoelastic property is modeled by using the generalized Maxwell-element rheological model to exhibit rate-dependent behaviors, such as creep and hysteresis on cyclic loading. The results of this investigation illustrate that elastomeric material as expected has large damping to reduce vibration and absorb energy, which leads to a reduction in ball–race contact forces and friction. A parametric study confirmed that the viscoelastic stress (VS) contributes significantly to the performance of the material, and without proper amount of viscoelasticity it loses its advantage in vibration reduction and exhibits linear elastic material characteristics. As expected, it is also demonstrated that housing supports made of linear elastic material provide minimal damping and rely on the bearing friction to dissipate energy. A study of housing support geometry demonstrates that bearing support plays a large role on the dynamic performance of the bearing. Motion of bearing outer race is closely related to the geometry and symmetry of the housing.

Sign in / Sign up

Export Citation Format

Share Document