Premium
Penalty function method for combined finite–discrete element systems comprising large number of separate bodies
Author(s) -
Munjiza A.,
Andrews K. R. F.
Publication year - 2000
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/1097-0207(20001220)49:11<1377::aid-nme6>3.0.co;2-b
Subject(s) - finite element method , discrete element method , discretization , mixed finite element method , computer science , extended finite element method , finite element limit analysis , penalty method , context (archaeology) , algorithm , mathematical optimization , mathematics , computational science , topology (electrical circuits) , mathematical analysis , engineering , structural engineering , mechanics , physics , paleontology , biology , combinatorics
Large‐scale discrete element simulations, the combined finite–discrete element method, DDA as well as a whole range of related methods, involve contact of a large number of separate bodies. In the context of the combined finite–discrete element method, each of these bodies is represented by a single discrete element which is then discretized into finite elements. The combined finite–discrete element method thus also involves algorithms dealing with fracture and fragmentation of individual discrete elements which result in ever changing topology and size of the problem. All these require complex algorithmic procedures and significant computational resources, especially in terms of CPU time. In this context, it is also necessary to have an efficient and robust algorithm for handling mechanical contact. In this work, a contact algorithm based on the penalty function method and incorporating contact kinematics preserving energy balance, is proposed and implemented into the combined finite element code. Copyright © 2000 John Wiley & Sons, Ltd.