We introduce a new class of inhomogeneous cosmological models as solutions to the Einstein-Maxwell equations in electrovacuum. The new models can be considered to be nonlinear perturbations, through an electromagnetic field, of the previously studied 'smooth Gowdy-symmetric generalised Taub-NUT solutions' in vacuum. Utilising methods from soliton theory, we analyse the effects of the Maxwell field on global properties of the solutions. In particular, we show existence of regular Cauchy horizons, and we investigate special singular cases in which curvature singularities form.