The copolymerization reaction of styrene with divinylbenzene creates a cross-linked polymer, which one of the main applications is its use as a filling in chromatography columns (gel permeation or ion exchange, for example) and it can also be used as a biocatalyst support, precisely because it has cross-linked chains and supports extreme reaction conditions. This polymerization occurs by forming cross-links between the polymeric chains, resulting in high thermal and mechanical stability. A characteristic of this reaction is the formation of a gel fraction which, in fact, is an insoluble copolymer with a high degree of cross-linking in relation to the polymeric product that is in solution (sol fraction) formed at the end of the process. Therefore, knowing the chain size distribution of the gel fraction is a challenge in this system, considering that its high molecular weight makes it impossible to solubilize a sample for gel permeation chromatography analysis (GPC). The development of kinetic models capable of predicting the contribution of the soluble and gel fractions to the molar mass distribution in solution copolymerization of this copolymer is essential for the mitigation of this gel formation, allowing the stable operation of continuous industrial reactors. Based on what has been presented, the main objective of this study is to investigate the kinetics of cross-linked copolymer formation via free radical solution polymerization. Additionally, it was also intended to develop a model capable of predicting the molecular weight distribution of the sol and gel fractions. To achieve these goals, the Monte Carlo method was used, based on the division of the copolymer molecules in linear chains and in crosslinking groups. The developed model was validated with experimental data collected from other studies and predicted the conversion data, average and weight molar masses, molar mass distribution with distinction of the soluble fraction and gel, showing as a result that the soluble phase of this system represents only chains that they underwent only one crosslinking. This proposed model can be used for monitoring and optimization purposes, such as finding experimental conditions in which the gel fraction is generated in smaller amounts.