In this paper, we will present a technique of modeling contact among virtual rheological objects. A collision of virtual rheological objects, which are modeled as a lattice structure model, is expressed using artificial forces between the edges and the nodes. We will develop the time-integrated penalty method, where an artificial force involves a time-integral term.
