Many key properties and applications of magnetic materials are strongly intertwined with the spatial distribution of magnetic moments inside the corresponding specimens. In addition to classical magnetic structures, magnetic skyrmions have raised interest in spintronics as carriers of information for future storage devices. In this talk, we present an advance in the mathematical modeling of magnetic skyrmions by analyzing the interplay of stochastic microstructures and chirality. Under the assumptions of stationarity and ergodicity, we characterize the Gamma-limit of a micromagnetic energy functional, including the Dzyaloshinkskii-Moriya contribution. Eventually, we present an explicit characterization of minimizers of the effective model in the case of magnetic multilayers. This talk is based on a joint work with E. Davoli and J. Ingmanns.