The microstructural origin of strain hardening during plastic deformation in stage II deformation of face-centered cubic (fcc) metals can be attributed to the increase in dislocation density resulting in a formation of dislocation networks. Although this is a well known relation, the complexity of dislocation multiplication processes and details about the formation of dislocation networks have recently been revealed by discrete dislocation dynamics (DDD) simulations. It has been observed that dislocations, after being generated by multiplication mechanisms, show a limited expansion within their slip plane before they get trapped in the network by dislocation reactions. This mechanism involves multiple slip systems and results in a heterogeneous dislocation network, which is not reflected in most dislocation-based continuum models. We approach the continuum modeling of dislocation networks by using data science methods to provide a link between discrete dislocations and the continuum level. For this purpose, we identify relevant correlations that feed into a model for dislocation networks in a dislocation-based continuum theory of plasticity. As a key feature, the model combines the dislocation multiplication with the limitation of the travel distance of dislocations by formation of stable dislocation junctions. The effective mobility of the network is determined by a range of dislocation spacings which reproduces the scattering travel distances of generated dislocation as observed in DDD. The model is applied to a high-symmetry fcc loading case and compared to DDD simulations. The results show a physically meaningful microstructural evolution, where the generation of new dislocations by multiplication mechanisms is counteracted by a formation of a stable dislocation network. In conjunction with DDD, we observe a steady state interplay of the different mechanisms.