Within the Discrete Element Method (DEM) framework, the multi-sphere model capable of creating irregular-shaped particles by several sub-spheres attracts increasing attention. Traditionally, the creation of irregular-shaped particles using the multi-sphere model requires manual determination of the position and diameter of each sub-sphere, which extremely increases the difficulty and workload of DEM modeling of irregular-shaped particles. To tackle this issue, a new automatic generation algorithm that helps to create arbitrary-shaped multi-sphere particles automatically was developed in this study. The quality of created multi-sphere particles can be easily controlled by adjusting the number of sub-spheres and size factor referred to as the size ratio of the smallest sub-sphere to the largest sub-sphere. To verify the effectiveness of the newly developed automatic generation algorithm, two verification tests were conducted in this study: (i) The packing behavior of irregular-shaped particles in a rectangular container and (ii) the mechanical behavior of irregular-shaped particles under the uniaxial compression. The fill heights of the simulations were in line with those of the experiments with an error of less than 8%. The mechanical behavior of irregular-shaped particles under the uniaxial compression can also be well reproduced. Besides, the influence of important parameters on simulation accuracy was also investigated. The results demonstrated that the newly developed automatic generation algorithm can accurately create complex-shaped particles with a small number of sub-spheres. Increasing the number of sub-spheres will improve the simulation accuracy but at the expense of computational efficiency. A larger size factor leads to higher computational efficiency in DEM simulation but is accompanied by a significant loss of accuracy.