Because of the complex connectivity of cross-linked polymers, generating structures for molecular simulations is a nontrivial task. In this work, a general methodology is presented for constructing post-cross-linked polymers by a new two-stage implementation of the Polymatic simulated polymerization algorithm, where linear polymers are first polymerized and then cross-linked. It is illustrated here for an example system of thermally cross-linked octene-styrene-divinylbenzene (OS-DVB) copolymers. In the molecular models, the degree of cross-linking is ranged from 0 to 100%, and the resulting structural and thermal properties are examined. The simulations reveal an increase in the free volume with higher cross-linking degrees. Shifts in the peaks of the structure factors, which are assigned to contributions from the backbone and side-chain atoms, correspond to the formation of larger free volume elements. Furthermore, the glass transition temperatures increase with higher degrees of cross-linking, while the thermal expansivity decreases. Comparisons with experimental results for similar systems are made when available. As demonstrated here, the presented methodology will provide an effective route to simulating post-cross-linked polymers for a variety of applications, which will enable an improved understanding of their structure-property relationships.