Fujiyabu et al. [1] experimentally reported that mixing of 3-arm star prepolymers into 4-arm analogs improves the toughness of the resultant polymer networks compared to that of the base network composed of only 4-arm star polymers. The present study investigated the mechanism of this phenomenon based on phantom chain simulations for polymer networks consisting of mixtures of star branch prepolymers with equal arm length and different arm numbers, (f(1), f(2))= (3, 4), (3, 6) and )3, 8), for various f(1) = 3 prepolymer fractions phi(3). These networks were created via end-linking reactions between prepolymers traced by a Brownian dynamics scheme, and the network structure was stored at different conversion ratios phi(c), ranging from 0.6 to 0.9. The cycle rank of the gelated networks xi is entirely consistent with that given by mean-field theory, demonstrating that the examined network structure is statistically valid. The networks were stretched with energy minimization until the break, and fracture characteristics including strain at break epsilon(b), stress at break phi(b), work for fracture W-b, and the ratio of broken strands phi(bb), were obtained. Plots of epsilon(b), sigma(b) = phi(bb), and W-b = phi(bb) versus xi roughly follow the master curves reported for the base networks without mixing, implying that the change of fracture properties by the mixing of f(1) = 3 mainly corresponds to a decrease in.. The mixing slightly suppresses sigma(b) = phi(bb) and W-b = phi(bb) for large f(2) cases compared to the base networks because of a biased breakage at the network strands without extenders, which are prepolymers with only two reacted arms. The analysis also revealed a new master curve for the xi-dependence of the molecular weight of broken strands.