Computational & Technology Resources
an online resource for computational,
engineering & technology publications |
|
Civil-Comp Proceedings
ISSN 1759-3433 CCP: 89
PROCEEDINGS OF THE SIXTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: M. Papadrakakis and B.H.V. Topping
Paper 129
Development of a Continuum-Consistent Discrete Element Model for Solids Using the Triangle Finite Element V. Vadluga and R. Kacianauskas
Laboratory of Numerical Modelling, Vilnius Gediminas Technical University, Lithuania V. Vadluga, R. Kacianauskas, "Development of a Continuum-Consistent Discrete Element Model for Solids Using the Triangle Finite Element", in M. Papadrakakis, B.H.V. Topping, (Editors), "Proceedings of the Sixth International Conference on Engineering Computational Technology", Civil-Comp Press, Stirlingshire, UK, Paper 129, 2008. doi:10.4203/ccp.89.129
Keywords: discrete element, lattice, finite element, heterogeneous solids, dynamics, fracture.
Summary
The concept of the discrete element method (DEM) basically refers to the original
work of Cundall and Strack [1] which aims to describe the macroscopic behaviour
of granular assemblies.
Recently, the DEM has been extensively applied to the simulation of heterogeneous solids to study their dynamic deformation behaviour and fracture problems. It could be stated, however, that a unified approach consistent with classical (non-rotational) continuum concept is still at the stage of development. The paper addresses the development of the lattice-type DE model for planar classical continuum by combining it with the finite element technique. In the framework of the developed combined approach, their computation is based on the analogy between structural and continuum variables in the triangle, employing the natural concept of finite element method (FEM) suggested by Argyris at al. [2]. The above concept considerably modifies the lattice-type DEM model. As opposed to the extensive developments exploring a lattice with central interaction as suggested by Griffiths and Mustoe [3], we arrive at the model with angular interaction [4], expressed, however, in the alternative form. Unlike earlier developments based on a separate shear stiffness concept [3], here, shear stiffness works only in combination with axial stiffness. Consequently, the model demonstrates higher diversity of the Poisson's ratio. Secondly, it is not related to particular triangle shapes and may be extended to various meshes composed of triangles and to three-dimensional. From the physical point of view, mutual interaction is extended to the neighbouring particles, while the new approach requires modification of the algorithm and the DEM code. The developed DEM model was implemented into the original DEM code DEMMAT [5] and validated by simulating the dynamic behaviour of the plane stress problem against the 'exact' FE results. Brittle cracking with randomly distributed tensile strength properties of the material is also considered as an application. References
purchase the full-text of this paper (price £20)
go to the previous paper |
|