Rubber elasticity: A contact-probability model with harmonic entanglement constraints