Algorithms and capabilities of {it Solidity} to simulate interactions and packing of complex shapes. A number of numerical algorithms for simulation of particle packing have been proposed and used in a wide range of industries: mining, chemical engineering, pharmaceuticals, agriculture and food handling, etc. However, most of them can only deal with simple and regular shapes due to the complex and expensive numerical algorithms needed to simulate complex shapes. In this paper, a FEMDEM code, {it Solidity}, is used to more accurately capture the influence of complex shape. It combines deformable fracturing arbitrary-shaped particle interactions modelled by FEM with discrete particulate motion modelled by DEM. This paper will cover recent code optimisation for the contact force calculation with arbitrary body shape, parallelisation performance and discussion of results showing both deformable and rigid body versions of the code in different application scenarios. {it Solidity} also provides post-processing tools to analyse the particle packing structure in terms of local porosity and orientation distributions, contact forces, and coordination number, etc. Some examples of Platonic and Archimedean body packs are presented.