С шаром, следственно, можно, в принципе, поступать по аналогии.
Нет, не совсем по аналогии, а ровно так, как Вы сказали в самом начале: считать распирание каждой следующей чутьсферки внутренним шаром и потом интегрировать по радиусу. Правда, у меня при этом получился несколько дикий (на мой взгляд) результат, я его даже и привести стесняюсь.
Я это и имел в виду: приравнять работу потере энергии поля

Второе слагаемое - это изменение (наверное, уменьшение) энергии поля внутри шара. Тут, похоже, свои тонкости. Логично вроде бы считать, что при возрастании

, распределение заряда внутри раздувшегося шара сохранит подобие; то есть каждая внутренняя частица шара пропорционально увеличит своё расстояние от центра. Полный заряд при этом, ес-нно, остаётся постоянным. Вообще-то энергия поля внутри равна

Я предлагаю ввести параметр

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

Мы можем рассматривать это выражение как функцию

. При этом
Замена

приводит интеграл к виду

Отсюда производная по

при

равна

Ну и наконец

Всё, устаааааааааал.