Anika> У такого астероида масса покоя будет 1800 петатонн... 
4/3*pi*40e3**3 = 2,68E14м³ = 2,11E15кг (2,12 тератонн).
Энергия на 40км/с будет E=mv²/2 = 1,69E24Дж = 4E11кт ТНТ = 403тератонны ТНТ.
Где-то ошибся

То есть по одному параметру ошибку сразу вижу, я плотность 5500 взял, а не 7874

Но это не даёт ошибку в 40 раз.
Сейчас ещё раз обратный расчёт проведу:
E = 10Тт ТНТ = 10E12 * 4,18E9 = 4,18E22 Дж.
m = 2*E/v² = 2*4,18E22/(40E3**2) = 5,23E13кг
V = m/ρ = 5,23E13/7874 = 6,64E9м³
R = (3*V/(4*pi))**(1/3) = 1166 метров.
Получается чуть более двух км. диаметром.
Upd: было сожрано парсером, подправил