Уважаемые коллеги! Предлагаю обсудить...
Обозначим искомый радиус шара через r, диаметр заданного цилиндра – через D. Тогда центры шаров будут располагаться внутри равностороннего цилиндра с диаметром, равным D – 2r, и такой же высотой. Назовём такой цилиндр цилиндром центров. Отметим, что размеры цилиндра центров не фиксированы, а зависят от радиуса шаров. Рассмотрим несколько сечений этого цилиндра плоскостями.
1. Сечением цилиндра центров плоскостью, перпендикулярной к его оси, является окружность с диаметром, равным D – 2r. Интуиция подсказывает, что при плотной упаковке шаров их центры будут находиться в вершинах равностороннего треугольника . В итоге задача сводится к нахождению стороны равностороннего треугольника, вписанного в окружность с радиусом, равным (D – 2r)/2 = D/2 – r. Искомый радиус равен половине этой стороны, и, согласно условию, D = 1.
Известно, что для равностороннего треугольника со стороной, равной a, радиус R описанной окружности определяется по формуле R = a/√3. В нашем случае R = 1/2 – r, a = 2r. Следовательно, 1/2 – r = 2r/√3,
2r/√3 + r = 1/2, r(1 + 2/√3) = 1/2, r = 1/2 : (1 + 2/√3) = 1/2 • (√3 + 2)/√3 = √3/(2(√3 + 2)) = √3(2 – √3)/2 ≈ 0,23205. Полученный результат можно рассматривать как частный случай решения задачи о наиболее экономном заполнении круга k одинаковыми кругами при k = 3; об этом результате указано, в частности, в [1, с.108].
На рис. 1 ниже приведено сечение заданного цилиндра (не цилиндра центров!) плоскостью, перпендикулярной к его оси (нормальное сечение). Радиусы заштрихованных кругов равны r = √3(2 – √3)/2, диаметр большого круга равен D. Центры заштрихованных кругов лежат на окружности с диаметром, равным D – 2r = 1 – √3(2 – √3).
Рис. 1
Значение r
1 ≈ 0,23205 является первой оценкой искомого радиуса шара. Очевидно, что эта оценка не даёт ответа на вопрос задачи.
2. Осевое сечение цилиндра центров (сечение, проходящее через его ось) представляет собой квадрат со сторонами, равными D – 2r = 1 – 2r. Расположив центры шаров в трёх вершинах квадрата, получаем
1 – 2r = 2r (расстояние между центрами двух одинаковых шаров равно их диаметру), 4r = 1, r = 1/4 = 0,25000. Принимая
r
2 = 0,25000, видим, что r
2 > r
1. Но и эта оценка не даёт ответа на поставленный вопрос.
Установим, при каком положении равностороннего треугольника внутри квадрата площадь треугольника будет наибольшей. Для этого привлечём на помощь интуицию. Она подсказывает, что одну из вершин треугольника необходимо поместить в вершину квадрата так, чтобы высота треугольника, проведённая из этой вершины, совпала с диагональю квадрата. Выполним рис. 2.
Рис. 2
Из рисунка видно, что D – 2r = 2r • cos φ, φ = 45º – 60º/2 = 15º, cos 15º = (√3 + 1)/(2√2) [2, с. 107];
1 – 2r = 2r(√3 + 1)/(2√2), 2r(√3 + 1)/(2√2) + 2r = 1, 2r((√3 + 1)/(2√2) + 1) = 1, r = 1/(2((√3 + 1)/(2√2) + 1)) =
= 1/((√3 + 1)/√2 + 2) = √2/(2√2 + √3 + 1) ≈ 0,25433. Получили оценку r
3 ≈ 0,25433 > r
2.
3. Рассмотрим сечение цилиндра центров плоскостью P
α (рис. 3).
Рис. 3
Это сечение представляет собой эллипс, малая ось которого равна D – 2r, а большая равна √2(D – 2r). Соответственно малая полуось равна (D – 2r)/2, а большая (D – 2r)/√2. Поместим центр одного шара в начале малой оси, а центры двух других шаров по концам большой оси. Тогда квадрат минимального расстояния между шарами (первым и вторым, первым и третьим) равен 4r
2 = (D – 2r)
2/4 + (D – 2r)
2/2, откуда находим 4r
2 = (1 – 2r)
2/4 + (1 – 2r)
2/2 = 3(1 – 2r)
2/4, 2r = √3(1 – 2r)/2 = √3/2 – r√3, r(2 + √3) = √3/2,
r = √3/(2(2 + √3)) = √3(2 – √3)/2 = (2√3 – 3)/2 ≈ 0,23205. Вновь получили оценку r
4 = r
1 ≈ 0,23205.
Если поместить центры шаров вдоль большой оси эллипса, то получаем следующий результат:
4r = √2(D – 2r), 4r = √2(1 – 2r) = √2 – 2√2r, r(4 + 2√2) = √2, r = √2/(4 + 2√2) = 1/(2(√2 + 1)) = (√2 – 1)/2 ≈ 0,20711, то есть самую малую оценку из всех полученных, r
5 ≈ 0,20711 < r
1.
Есть основания предположить, что на этом можно остановиться, и оценка r
3 ≈ 0,25433 даёт ответ на поставленный вопрос. Обоснованием служит сравнение площадей сечений заданного цилиндра, соответствующих рассмотренным сечениям цилиндра центров. Имеем следующее:
1) в сечении круг с диаметром D = 1 и с площадью S
1 = π/4 ≈ 0,785; для него r = r
1 ≈ 0,23205;
2) в сечении квадрат со стороной D = 1 и с площадью S
2 = 1; для него r = r
3 ≈ 0,25433;
3) в сечении эллипс с полуосями a = 1/2, b = 1/√2 и с площадью S
3 = π/(2√2) ≈ 1,111; для него r = r
1 ≈ 0,23205.
Как видим, наибольший простор для размещения шаров предоставляет квадратное сечение заданного цилиндра. Несмотря на то, что эллипс имеет в данном случае бОльшую площадь (S
3 > S
2), радиус шаров для эллипса получается меньшим, чем для квадрата.
Если рассмотреть ещё, например, сечение цилиндра центров плоскостью P
β (рис. 4), то ему будет соответствовать сечение заданного цилиндра, представляющее собой половину эллипса, причём малая ось его равна D (соответственно, малая полуось равна D/2 = 1/2), а большая полуось равна √(D2/4 + D2) = D√3/2 = √3/2. Следовательно, площадь такого сечения равна S
4 = 1/2 • π • 1/2 • √3/2 = π√3/8 ≈ 0,680 < S
3 = 1,111 и S
4 < S
2.
Поэтому дальнейшего рассмотрения сечений заданного эллипса можно не проводить.
Рис. 4
В итоге получаем r = 0,25433.
Литература
1. Шклярский Д. О., Ченцов Н. Н., Яглом И. М. Геометрические оценки и задачи из комбинаторной геометрии. М.: Наука, 1974.
2. Цыпкин А. Г., Цыпкин Г. Г. Математические формулы. М.: Наука, 1985.
Никто не имеет замечаний?
С уважением.