Esta es una pregunta difícil de resolver en general, pero he visto el siguiente tipo de enfoque adoptado con unos explícita timestepping "contacto algoritmos de" entre dos cuerpos en el de elementos finitos tipo de literatura si no hay necesidad de ser una distinción entre la estática y la dinámica de los coeficientes de fricción.
Puede ser visto como una forma de hacer lo que ha sido mencionado por DilithiumMatrix y David Hammen.
En primer lugar, hemos proyecto de impulso a $p$ de avance para cada cuerpo como si no hubiera fricción, por lo que
$p^{*} = p^{n} + \Delta t f^n$.
Si no hay deslizamiento de contacto, que significa al final de el paso, el cuerpo se mueve con la velocidad de la combinación de centro de masa (en este caso, supongo que el suelo no se va a mover por el cuerpo, por lo que podemos usar 0).
Entonces sabemos
$p^{n+1} = p^{*} + \Delta t f^n_c$.
Por lo tanto, la fuerza en el cuerpo para no-slip de contacto es $f_c = \frac{-p^*}{\Delta t}$ desde $p^{n+1} = 0$.
Ahora, supongo que usted sabe la fuerza normal (por ejemplo, debido a que el peso del cuerpo), en cuyo caso se puede ver si $|\mu f^n_{normal}| > |f_c|$, en cuyo caso, el objeto debe venir a una parada. De lo contrario, se debe calcular la fuerza de contacto como
$f_{contact} = -\mathrm{sign}(p^n) \mu f_{normal}$
Y el final con el impulso dado por $p^{n+1} = p^* + \Delta t f_{contact}$. O, con un total de
$p^{n+1} = \begin{cases} 0 & \text{ if } |\mu f_{normal}| \ge |\frac{-p^{n}}{\Delta t} - f^n| \\ p^{n} + \Delta t (f^n -\mathrm{sign}(p^n) \mu f_{normal}) & \text{otherwise}\end{cases}$.
Por favor, tenga en cuenta que he dejado en 1D y esto sólo representa el movimiento tangencial.
Usted podría ser capaz de establecer los umbrales para $\mu$ donde se muestra como dos valores diferentes, para obtener la estática y dinámica de los coeficientes de fricción pero no lo he probado. E. g. dado que el coeficiente de fricción estática es mayor que el dinámico, creo
$p^{n+1} = \begin{cases} 0 & \text{ if } |\mu_{static} f_{normal}| \ge |\frac{- p^{n}}{\Delta t} - f^n| \\ p^{n} + \Delta t (f^n -\mathrm{sign}(p^n) \mu_{dynamic} f_{normal}) & \text{otherwise}\end{cases}$
funcionará sin pasarse, pero todavía va a detener "demasiado pronto" si el objeto está en movimiento, pero el no-deslizamiento de la fuerza fue entre la estática y la dinámica de los umbrales.