In the context of the de Broglie-Bohm pilot wave theory, numerical simulations for simple systems have shown that states that are initially out of quantum equilibrium - thus violating the Born rule - usually relax over time to the expected $|psi|^2$ distribution on a coarse-grained level. We analyze the relaxation of nonequilibrium initial distributions for a system of coupled one-dimensional harmonic oscillators in which the coupling depends explicitly on time through numerical simulations, focusing in the influence of different parameters such as the number of modes, the coarse-graining length and the coupling constant. We show that in general the system studied here tends to equilibrium, but the relaxation can be retarded depending on the values of the parameters, particularly to the one related to the strength of the interaction. Possible implications on the detection of relic nonequilibrium systems are discussed.