We present a general class of solutions to Einstein's field equations with two spacelike commuting Killing vectors by assuming the separation of variables of the metric components. The solutions can be interpreted as inhomogeneous cosmological models. We show that the singularity structure of the solutions varies depending on the different particular choices of the parameters and metric functions. There exist solutions with a universal big-bang singularity, solutions with timelike singularities in the Weyl tensor only, solutions with singularities in both the Ricci and the Weyl tensors, and also singularity-free solutions. We prove that the singularity-free solutions have a well-defined cylindrical symmetry and that they are generalizations of other singularity-free solutions obtained recently.