We develop and solve a new molecular model for nonlinear elasticity of entangled polymer networks. This model combines and generalizes several succeseful ideas introduced over the years in the field of the rubber elasticity. The topological constraints imposed by the neighboring network chains on a given network are represented by the confining potential that changes upon network deformation. This topological potential restricts fluctuations of the network chain to the nonaffinely deformed confining tube. Network chains are allowed to fluctuate and redistribute their length along the contour of their confining tubes. The dependence of the stress sigma on the elongation coefficient A for the uniaxially deformed network is usially represented in the form of the Mooney stress, f*(1/lambda) = sigma/(lambda - 1/lambda(2)). We find a simple expression for the Mooney stress, f*(1/lambda) = G(c) + G(e)/(0.74lambda + 0.61lambda(-1/2) - 0.35), where G(c) and G(e) are phantom and entangled network moduli. This allows one to analyze the experimental data in the form of the universal plot and to obtain the two moduli G(c) and G(e) related to the densities of cross-links and entanglements of the individual networks. The predictions of our new model are in good agreement with experimental data for uniaxially deformed polybutadiene, poly(dimethylsiloxane), and natural rubber networks, as well as with recent computer simulations.