Finite strain, finite rotation quadratic tetrahedral element for the combined finite–discrete element method

Abstract In the past, the combined finite–discrete element was mostly based on linear tetrahedral finite elements. Locking problems associated with this element can seriously degrade the accuracy of their simulations. In this work an efficient ten‐noded quadratic element is developed in a format suitable for the combined finite–discrete element method (FEMDEM). The so‐called F‐bar approach is used to relax volumetric locking and an explicit finite element analysis is employed. A thorough validation of the numerical method is presented including five static and four dynamic examples with different loading, boundary conditions, and materials. The advantages of the new higher‐order tetrahedral element are illustrated when brought together with contact detection and contact interaction capability within a new fully 3D FEMDEM formulation. An application comparing stresses generated within two drop experiments involving different unit specimens called Vcross and VRcross is shown. The Vcross and VRcross units of ∼3.5 × 10 4 kg show very different stress generation implying different survivability upon collision with a deformable floor. The test case shows the FEMDEM method has the capability to tackle the dynamics of complex‐shaped geometries and massive multi‐body granular systems typical of concrete armour and rock armour layers. Copyright © 2009 John Wiley & Sons, Ltd.