For the simulation of particles with large elastic-plastic deformation and breakage, two techniques have been traditionally used. The first technique is Finite Element Method (FEM) [1] while the second one is the Discrete Element Method (DEM) [2] with composite particles created by connecting simple DEM particles via elastic joints. Each approach has its advantages and disadvantages. FEM handles large elastic-plastic deformations easily but struggles with breakage simulations and requires large computational effort to assemble and solve stiffness matrix. DEM can handle particle breakage quickly and efficiently but struggles with large elastic-plastic deformations for shapes other than fibers and shells. A new approach was programmed in ANSYS Rocky that combines FEM and DEM. The approach is based on the algorithm described in [3]. The simulated particles are created from the simple particle shapes connected by joints like in the DEM approach. However, these simple particles are also allowed to deform subject to forces applied by the joints and particle-particle and particle-boundary contacts. A stiffness matrix is solved only on the element level making it ideal for the multi-GPU code implementation. Several examples of the new code usage will be presented. REFERENCES [1] ANSYS FEM Reference Manual. [2] ANSYS Rocky User Manual. [3] A. Potapov and C.S. Campbell, “A hybrid finite-element simulation of solid fracture”, Int. J. of Mod. Phys. C, 7 (02), 155-180 (1996).