Com δn>0 confirmado, o algoritmo precisa direcionar a força ao longo da linha de contato e devolver essa informação à integração do movimento.
5
Vetor unitário normal (n̂ij)
Para direcionar a força ao longo da linha de afastamento, divide-se o vetor do Passo 2 pelo seu módulo: n̂ij = (Xi−Xj)/dij
6
Força de contato
A sobreposição do Passo 3 entra na equação não-linear de Hertz: F⃗n,ij = [(4/3·E*·√R*)·δn1,5]·n̂ij (componente elástica; o mesmo δn é derivado no tempo para o amortecedor).
7
Nova aceleração
A força resultante (somada à gravidade e às forças de outros vizinhos) é dividida pela massa: a⃗i = ΣF⃗ij/mi + g⃗ A partir da aceleração, atualiza-se a velocidade, avança-se o tempo e retorna-se ao Passo 1.
🔁
Resumo: o que na realidade é a deformação do material é traduzido no computador como a subtração geométrica entre a distância dos centros e a soma dos raios (δn = Ri+Rj−dij). A cada incremento de tempo, a distância dij encurta, o δn aumenta e gera uma força de repulsão de Hertz cada vez maior, até frear as partículas.
Prof. Cláudio Roberto Duarte — UFU · DEM & LIGGGHTSFonte: dem_detalhamento.php — Tópico 26