Abstract: A computer-implemented method for simulating a contact state of a rubber material with a contact surface is disclosed, wherein a rubber model of the rubber material is defined; the rubber material is modeled by a finite number of particles on which a motion equation is defined; a contact zone model of a contact zone including the contact surface is defined; the contact zone is modeled by a finite number of particles; the rubber model is contacted with the contact zone model and a deformation calculation of the rubber model is performed; and the motion equations defined on the particles of the rubber material have a term expressive of resilience of the rubber material.