The increasing sinusoidal simulations varied the box length L x from its initial length L x 0 according to the equation The results are shown in Fig. , which shows the fraction of broken bonds as a function of (a) amplitude and (b) frequency of the applied oscillation. Fig. shows the effect of varying frequency. At low amplitude (A = 1.1) the amount of damage is very low and approximately independent of ω. Put simply, if no damage occurs in an oscillation then the system can oscillate in perpertuity. The larger amplitude lines show that more damage has occurred by the end of the simulation; the damage is the least for the lowest frequency oscillations and the most for the highest frequency oscillations. The link between oscillation frequency and damage is roughly linear, although there are no samples between ω = 2 and ω = 5. In this range there is not an observable "tailing off" of the damage as ω increases, but instead damage increases monotonically with ω. The increased damage at higher frequencies is due to the increase in mechanical work being done to stretch the network, which can be used to break bonds within the network. We can also see the effect of increased mechanical work in Fig. , where a larger stretching amplitude also leads to a larger amount of damage at each stretching frequency. The linear relation could be explained as the stretching forces are given by dU dt ∝ Aω cos(ωt), which are maximised at ∝ Aω at peaks and Fig. Recording the fraction of broken bonds at the end of a simulation as a function of stretching frequency ω and amplitude A. The fraction of broken bonds at the end of the simulation increases approximately linearly with both frequency and amplitude, as the most damage is introduced at the maxima and minima of a stretching cycle. The energy inserted into the network at the maxima and minima of stretching is linearly related to the frequency and amplitude.