In this contribution, we have constructed new analytical solutions for initial data of the Einstein equations. Such solutions are valuable for gaining a better understanding of problems involving strong gravitational and/or electromagnetic interactions in general relativity. In this process we have examined an inhomogeneous cosmological model consisting a lattice of regularly arranged, charged black holes with initial data corresponding to the maximum expansion of a cosmological solution. We have also refined the method in such a way that the values of the mass and charge of the sources can be prescribed beforehand subject to certain constraints dictated by the field equations.
