Cálculo do tempo de colisão com aceleração

De Física Computacional
Ir para navegação Ir para pesquisar

Tendo em mãos a condição para que ocorra uma colisão: |ri(t+dt)rj(t+dt)|=σ, e que ri(t+dt)=ri+vidt+gdt22 e rj(t+dt)=rj+vjdt+gdt22, podemos elevar os dois lados da equação da condição ao quadrado, retirar o módulo uma vez que o resultado será sempre positivo, obtendo:
(ri(t+dt)rj(t+dt))2=σ2

ri2+vi2dt2+g2dt44+2rividt+2rigdt22+2vidt322(ri+vidt+gdt22)(rj+vjdt+gdt22)+rj2+vj2dt2+g2dt44+2rjvjdt+2rjgdt22+2vjdt32=σ2

ri2+vi2dt2+g2dt42+2rividt+2rigdt22+2vidt322rirj2rivjdt2rigdt222rjvidt2vivjdt22vigdt322rjgdt222vjgdt32g2dt42+rj2+vj2dt2+2rjvjdt+2rjgdt22+2vjgdt32=σ2

Reajeitando os termos antes da simplificação final:
ri22rirj+rj2+vi2dt22vivjdt2+vj2dt2+g2dt42g2dt42+2rigdt222rigdt22+2vidt322vigdt322rjgdt22+2rjgdt222vjgdt32+2vjgdt32+2ri(vivj)dt2rj(vivj)dt=σ2

Ao simplificar obtemos a mesma equação quadrática para se calcular o dtmin, onde ficou definido Δr=rirj e Δv=vivj na seção anterior:
(rirj)2+(vivj)2dt2+2(rirj)(vivj)dtσ2=0