Close tidal encounters among large planetesimals and moons should have been more common than grazing or normal impacts. Using a mass spring model within an N-body simulation, we simulate the deformation of the surface of an elastic spherical body caused by a close parabolic tidal encounter with a body that has similar mass as that of the primary body. Such an encounter can induce sufficient stress on the surface to cause brittle failure of an icy crust and simulated fractures can extend a large fraction of the radius of body. Strong tidal encounters may be responsible for the formation of long graben complexes and chasmata in ancient terrain of icy moons such as Dione, Tethys, Ariel and Charon.

We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number $\,k_2\,$ divided by the tidal quality factor $\,Q\,$) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of $\,k_2/Q\,$, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky (2012a; CeMDA 112$\,:\,$283) for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalised to inhomogeneous and/or non-spherical bodies.