A numerical algorithm is presented for the filling process of a cylindrical column with equilateral cylinders. The process is based on simplified mechanics - the elements are added one by one until the mechanical equilibrium is reached. The final structure is examined with respect to the global and local porosity distribution. Oscillating radial porosity profile is obtained in accordance with experimental data.