Как написано в википедии в статье "Общая теория относительности": «Если гравитационная масса точно равна инерционной, то в выражении для ускорения тела, на которое действуют лишь гравитационные силы, обе массы сокращаются. Поэтому ускорение тела, а следовательно, и его траектория не зависит от массы и внутреннего строения тела. Если же все тела в одной и той же точке пространства получают одинаковое ускорение, то это ускорение можно связать не со свойствами тел, а со свойствами самого́ пространства в этой точке.»
Но также можно предположить, что не только ускорения но и само движение тела связано с состоянием пространства в той области пространства, в которой находится тело.
Рассмотрим неподвижное сферическое тело. Вокруг него существует гравитационное поле. Его эквипотенциальная поверхность представляет собой сферу, в центре которой находится тело. Если тело поместить в центр координат, то уравнение такой поверхности

Вокруг движущегося тела также существует гравитационное поле. Но поскольку существует максимальная скорость передачи взаимодействий, то эквипотенциальная поверхность уже не будет симметрична относительно тела. В первом приближении ее уравнение (тело в центре координат и движется вдоль оси х) будет

, где

– скорость тела,

– максимальная скорость передачи взаимодействий.
Именно эта асимметрия гравитационного поля и вызывает движение тела. Тело деформирует пространство, пространство воздействует на тело.
Поскольку масса сейчас вроде не играет роли гравитационного заряда, то во избежание путаницы введем физическую величину «гравитационный заряд»

.
Если предположить, что закон всемирного тяготения через гравитационный заряд выражается следующим образом

, то единичным гравитационным зарядом будет обладать тело, сила притяжения к которому на единичном расстоянии будет равна гравитационной постоянной

В качестве гравитационного потенциала будем использовать величину

, тогда радиус поверхности с единичным потенциалом вокруг тела с единичным гравитационным зарядом будет равен единице. А уравнение такой поверхности вокруг любого сферического тела

Импульс тела с точностью до множителя будет равен площади сферического пояса единичной эквипотенциальной поверхности, заключенного между плоскостями, перпендикулярными направлению движения и проходящими через тело и через центр поверхности.

. Кинетическая энергия, также с точностью до множителя, будет равна площади сферы с радиусом, равным расстоянию между телом и центром единичной эквипотенциальной поверхности.

Рассмотрим неупругое соударение двух тел с единичным гравитационным зарядом, одно из которых покоится. В момент соударения единичную эквипотенциальную поверхность можно построить исходя из следующих уравнений


где

– расстояние точки единичной эквипотенциальной поверхности до неподвижного тела

– расстояние точки единичной эквипотенциальной поверхности до центра эквипотенциальной поверхности движущегося тела

– угол между линией соединяющей точку поверхности с неподвижным телом и направлением движения второго тела

– угол между

и


– гравитационный потенциал в точке.
Систему уравнений я решал численным методом при помощи итерационной процедуры. Затем находил плоскость перпендикулярную направление движения второго тела, которая делит полученную поверхность на две части равной площади. А потом находил площадь сферического пояса между этой плоскостью и параллельной ей плоскостью, проходящей через неподвижное тело. Скорость двух тел как единого целого после столкновения как доля от скорости света

будет равна площади данного сферического пояса деленной на половину площади единичной эквипотенциальной поверхности.
В результате получилось следующее:
При скорости движущегося тела

0,1 скорость после столкновения 0,05000283. При 0,3 - 0,150762; 0,5 - 0,252908; 0,7 - 0,351875; 0,9 - 0,427935.
В первом приближении вполне соответствует закону сохранения импульса.