Квант

33

Момент инерции в геометрии

Кандидат физико-математических наук
В. Н. ДУБРОВСКИЙ

Все хорошо знают, что математика постояно используется в физике. Режке встречается и потому менее известно обратное — использование физических понятий в математике. Один яркий пример такого рода — центр масс в геометрии — обсуждался в недавней стате М. Б. Балка и В. Г. Болтянского («Квант» № 4, с. 18). В нашей стате главную роль также играет понятие, занимаванное из механики, — момент инерции.

Но когда с понятием, взятым из физики, начинает работать математик, он быстро забывает его физический смысл и видит в нем лишь полезную математическую формулу. Так, понятие матерIALНОЙ ТОЧКИ А массы m, которое физик воспринимает как «тело массы m пренебрежимо малых размеров для данной задачи», математик определяет формально, как «пару (A, m)», состоящую из точки A (геометрической) и действительного числа m∈R». Физик дважды подумает, прежде чем воспользовется «отрицательной массой», а математик будет рассматривать как ни в чем не бывало «материальный точку»/(A, m), для которой m<0! И притом с пользователем для дела (см. цитированную выше статую, например задачу 5 на с. 19). Подобным же образом мы здесь определяем момент инерции чисто математическим выражением. А выведенная

34

затем формула Лагранжа для вычисления момента инерции станет ключом в своевобразную кунсткамеру геометрии, с множеством диковинных теorem и формул. Об этом наш расссказ.

Предваря возможный вопрос любозательного читателя, посямим все же физический смысл момента инерции. Это легко сделать, отправляясь от второго закона Ньютона. Привращаемым движением второй закон Ньютона записьается в виде M=IαM=I\alpha, где MM — момент прилагаемой силы, α\alpha — угловое ускорение, а II и есть как раз момент инерции (см. «Квант», 1971, № 1, с. 6). Сравнивая эту формулу с F=maF=ma, мы видим, что момент инерции представляется собой меру инерции врашающего тела, подобно тому, как масса представляется собой меру инерции при поступательном движении.

Что нужно знать о центре масс

Коротко напомним все необходимые нам сведения о центре масс (продобности можно найти в упомянутой статье М. Б. Балка и В. Г. Болтянского).

Рассмотрим систему из nn материалов точек (или точных масс), то есть nn пар (A1,m1),...,(An,mn)(A_1, m_1), ..., (A_n, m_n), где AiA_i — произвольные точки пространства, mim_i — действительные числа («массы», сосредоточенные в этих точках); более короткое и удобное обозначение для этой системы — m1A1,...,mnAnm_1A_1, ..., m_nA_n. Если общая масса m=m1+...+mn≠0m=m_1+...+m_n\neq 0, то

  1. существует и единственная такая точка ZZ, что

m1ZA→1+...+mnZA→n=0→;m_1\overrightarrow{ZA}_1+...+m_n\overrightarrow{ZA}_n=\overrightarrow{0};

эта точка называется центром масс системы m1A1,...,mnAnm_1A_1, ..., m_nA_n;

  1. для любой точки PP имеет место равенство

PZ→=1m(m1PA→1+...+mnPA→n);\overrightarrow{PZ}=\frac{1}{m}(m_1\overrightarrow{PA}_1+...+m_n\overrightarrow{PA}_n);

(1)

Рис. 1.

  1. (теорема о группировке) при замене части материалов точек системы m1A1,...,mnAnm_1A_1, ..., m_nA_n, например точек m1A1,...,mkAkm_1A_1, ..., m_kA_k, на их центр масс CC, в которой помещена их общая масса m′=m1+...+mkm'=m_1+...+m_k, положение центра масс всей системы не меняется, то есть центры масс системы m′C,mk+1Ak+1,...,mnAnm'C, m_{k+1}A_{k+1}, ..., m_nA_n и m1A1,...,mnAnm_1A_1, ..., m_nA_n совпадают.

Неопределенно из определения следует, что центр ZZ двух точных масс m1A1m_1A_1 и m2A2m_2A_2 находится по известному правилу рычаза: точка ZZ лежит на пряной части — коэффициент пропорциональности двух коллинеарных векторов).

Задачи

  1. Пусть a,b,ca, b, c — длины сторон BC,CA,ABBC, CA, AB треугольника ABCABC. Доказите, что центр его вписанной окружности является центром масс системы aA,bB,cCaA, bB, cC. (Указание: пользователь тем, что биссертиса треугольника делит его сторону на отрезки, пропорциональные прилежащим к ним сторонам ∗^*), доказите, что, например, барисентр материалов точек в bBbB и cCcC — основание биссертисы угла AA; примените теорему о группировке.)

  2. Рассмотрим три внешисанные окружности треугольника ABCABC — окружности, касающиеся (извне) одной из сторон треугольника и продолженной двух других сторон (рис. 1). Пусть A1,B1,C1A_1, B_1, C_1 — точки их касания со стороны BC,CA,ABBC, CA, AB. Доказите, что а) отрезки AA1,BB1AA_1, BB_1 и CC1CC_1 пересекаются в одной точке NN (точке Наделя), б) точка пересечения медиан MM треугольника лежит на отрезке ZNZN, где ZZ — центр вписанной окружности, и ∣ZM∣:∣MN∣=1:2|ZM|:|MN|=1:2.

Мoment инерции и формула Лагранжа

Мomentом инерции системы материалов точек

m1A1,m2A2,...,mnAnm_1A_1, m_2A_2, ..., m_nA_n

относительно фиксированной точки PP называется число

IP=m1∣PA1∣2+m2∣PA2∣2+...+mn∣PAn∣2.I_P=m_1|PA_1|^2+m_2|PA_2|^2+...+m_n|PA_n|^2.

Впервые это понятие было введено вликом Л. Эйлером (1707—1783).

Для любителей физики отметим, что в механике обычно рассматривается, собственно говоря, момент инерции системы точек относительно

*) См., например, статью Сефибекова С. Р. «Четыре доказательства теоремы о биссертисе» (Квант, 1983, № 8, с. 37).

**) См. также книгу Н. Н. Воробьева «Числа фибоначчи» (М.: Наука, 1978, с. 15).

35

Основные от точки AiA_i до этой оси (его произведение на угловую скорость вращения системы вокруг оси равно момент имплуса системы).

Он совпадает с IPI_P, когда точка PP и все точки AiA_i лежат в одной плоскости, а ось перепределения её и проходит через PP.

Формулу (2) можно записать короче, обозначив расстояние от AiA_i до PP через rir_i и воспользовавшись знаком суммы Σ∗\Sigma^*):

IP=∑imiri2.I_P = \sum_i m_i r_i^2.

(3)

Формула (3) и (2) можно придать другую форму. Пусть PP и QQ — две произвольные точки, тогда

IP=∑imi∣PAi∣2=∑imi(PQ→+QAi→)2I_P = \sum_i m_i |PA_i|^2 = \sum_i m_i (\overrightarrow{PQ} + \overrightarrow{QA_i})^2

=∑imi∣PQ∣2+2∑imi⋅PQ→⋅QAi→+∑imi∣QAi∣2,= \sum_i m_i |PQ|^2 + 2\sum_i m_i \cdot \overrightarrow{PQ} \cdot \overrightarrow{QA_i} + \sum_i m_i |QA_i|^2,

то есть

IP=IQ+m∣PQ∣2+2PQ→⋅(∑imiQAi→).I_P = I_Q + m |PQ|^2 + 2\overrightarrow{PQ} \cdot (\sum_i m_i \overrightarrow{QA_i}).

(4)

Если теперь в качестве QQ взять центр масс ZZ нашей системы, то, поскольку ∑imiZAi→=0→\sum_i m_i \overrightarrow{ZA_i} = \overrightarrow{0} (по определению центра масс), мы получим

IP=IZ+m∣PZ∣2.I_P = I_Z + m |PZ|^2.

(5)

В конце статьи мы показем, что момент инерции относительно центра масс выражается, если m≠0m \neq 0, следующей красивой формулой:

IZ=1m∑imimj∣AiAj∣2I_Z = \frac{1}{m} \sum_i m_i m_j |A_i A_j|^2

(6)

и, следовательно, для любой точки PP

IP=(1m∑imimj∣AiAj∣2)+m∣PZ∣2.I_P = \left( \frac{1}{m} \sum_i m_i m_j |A_i A_j|^2 \right) + m |PZ|^2.

(7)

Формулы (5)—(7) называются формулами Лаеранжа, на них основываются все дальнейшее.

*) Этот знак позволяет кратко записывать суммы из многих слагаемых: вместе A1+A2+A3+⋯+AnA_1 + A_2 + A_3 + \dots + A_n мы пишем ∑iAi\sum_i A_i или просто ∑iAi\sum_i A_i (когда число слагаемых nn ясно из контекста). Аналогично, ∑imiri2\sum_i m_i r_i^2 — сокращение для m1r12+m2r22+⋯+mnrn2m_1 r_1^2 + m_2 r_2^2 + \dots + m_n r_n^2 и т. п.

Рис. 2.

Математика — это исчисство называется веши одинаково, а одинаковые — по-разному

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

Перепишем ее в «нейтральных» обозначениях, отвлекаясь от смысла входящих в нее величины. Пусть x1,…,xnx_1, \dots, x_n — какой-то набор векторов, c1,…,cnc_1, \dots, c_n — набор чисел, c=c1+⋯+cn≠0c = c_1 + \dots + c_n \neq 0, x0=(c1x1+⋯+cnxn)/cx_0 = (c_1 x_1 + \dots + c_n x_n) / c, yi=xi−x0y_i = x_i - x_0 при i=1,…,ni = 1, \dots, n. Тогда

c1x12+⋯+cnxn2=(c1y12+⋯+cnyn2)+cx02c_1 x_1^2 + \dots + c_n x_n^2 = (c_1 y_1^2 + \dots + c_n y_n^2) + cx_0^2

(8)

(xi2,yi2x_i^2, y_i^2 нужно понимать, как складные квадраты векторов xi,yix_i, y_i). Доказывать эту формулу нам уже не нужно, потому что всегда можно считать, что xi=PAi→x_i = \overrightarrow{PA_i}, тогда x0=PZ→x_0 = \overrightarrow{PZ}, где ZZ — магнитный системе материалов точек c1A1,…,cnAnc_1 A_1, \dots, c_n A_n, yi=ZAi→y_i = \overrightarrow{ZA_i} и формула (8) превращается в формулу Лаеранжа (5).

Теперь приведем еще три интерпретации этой формулы.

  1. Пусть x1,…,xnx_1, \dots, x_n — векторы перепределикляров, опущенных из точек A1,…,AnA_1, \dots, A_n на ось ll (рис. 2), c1,…,cnc_1, \dots, c_n — массы, помещенные в эти точки. Тогда x0x_0 — это вектор перепределикляра из их барицентра ZZ на ось ll (почему?), а yiy_i — векторы перепределикляров из точек AiA_i на ось l0l_0, проведенный через

36

точку ZZ паралельно ll. В этом случае (8) принимает вид теоремы Штейнера:

Il=Il0+md2I_l = I_{l_0} + md^2

(здесь IlI_l и Il0I_{l_0} — моменты инерции системы c1A1c_1A_1, …, cnAnc_nA_n относительно осей ll и l0l_0, m=cm=c — общая масса, dd — расстояние между осями).

  1. Пусть nn частиц A1A_1, …, AnA_n с массами m1m_1, …, mnm_n движутся равномерно и пря-молинейно со скоростью v⃗1\vec{v}_1, …, v⃗n\vec{v}_n относительно некоторой системы отчета SS. Тогда их центр масс движется в системе SS со скоростью v⃗0=(m1v⃗1+⋯+mnv⃗n)/m\vec{v}_0 = (m_1\vec{v}_1 + \dots + m_n\vec{v}_n) / m (m=m1+⋯+mnm = m_1 + \dots + m_n). Скорость частицы AiA_i в системе покоя центра масс по правлу сложения скоростей равна u⃗i=v⃗i−v⃗0\vec{u}_i = \vec{v}_i - \vec{v}_0. Пусть EE и E0E_0 — полные кинетические энергии всех части в системе SS и в системе центра масс, то есть E=m1v122+⋯+mnvn22E = \frac{m_1v_1^2}{2} + \dots + \frac{m_nv_n^2}{2}, E0=m1u122+⋯+mnun22E_0 = \frac{m_1u_1^2}{2} + \dots + \frac{m_nu_n^2}{2}. Полягая в (8) xi=v⃗ix_i = \vec{v}_i, ci=mi/2c_i = m_i/2, получим теорему Кёнца о кинетической энергии:

E=E0+mv022E = E_0 + \frac{mv_0^2}{2}

(полная кинетическая энергия системы части равна энергии их движения относительно центра масс плоскости энергия движения центра масс).

Последний пример — из совершенного другой области — теории вероятности.

  1. Пусть в некотором эксперименте со случайным исходом измеряется некоторая величина XX, например, бросают игольную кость и определяют выпавшее число очков. Предположим, что, в зависимости от результата эксперимента, величина XX может принимать nn значений x1x_1, …, xnx_n (в примере с костью — 6 значений 1, 2, …, 6), причем значение xix_i принимается с вероятностью pip_i (числа pip_i неотрицательны и дают в сумме 1; в нашем примере, оченьindo, все значения равноверятны, то есть p1=p2=⋯=p6=1/6p_1 = p_2 = \dots = p_6 = 1/6). Число p1x1+⋯+pnxnp_1x_1 + \dots + p_nx_n называется средним значением (или математическим ожиданием) случайных XX и обозначается MXMX (в нашем примере MX=(1+⋯+6)/6=3,5MX = (1 + \dots + 6)/6 = 3,5). Таким образом, с точки зрения теории вероятностей выражения, входящие в (8) (при ci=pic_i = p_i, c=Σpi=1c = \Sigma p_i = 1), можно интерпретировать так*):

p1x12+⋯+pnxn2=MX2p_1x_1^2 + \dots + p_nx_n^2 = MX^2

— среднее значение случайной величины X2X^2 или, как говоря, второй момент (снова «момент»!) случайной величины XX; yi=xi−x0=xi−MXy_i = x_i - x_0 = x_i - MX — отклонения значений случайных XX от среднего значения;

p1y12+⋯+pnyn2=MX(X−MX)2p_1y_1^2 + \dots + p_ny_n^2 = MX(X - MX)^2

— среднее значение квадрата отклонения XX от MXMX; это очень важная характеристика случайной величины, показывающая степень ее «разброса» вокруг среднего значения, она называется диисперсей и обозначается DXDX; x02=(MX)2x_0^2 = (MX)^2.

В результате (8) превращается в формулу для вычисления дисперсии, имеющейся во всех учебниках теории вероятности (правда, не носшуюничих звучных имен):

DX=MX2−(MX)2DX = MX^2 - (MX)^2

Легко ли здесь узнать формулу Лагранжа?

Задча 3. Найдите дисперсию числа окнов при бросании кости.

Геометрические следствия

Формула Стоярта. Найдем расстояние dd от вершины AA треугольника ABCABC до точки DD на стороне BCBC, знаний длины aa, bb, cc сторон треугольника и длины отрезков BDBD и DCDC (рис. 3). Пусть p=∣BD∣:∣BC∣p = |BD| : |BC|, q=∣DC∣:∣BC∣q = |DC| : |BC|, тогда DD — это центр масс системы материалов точек qBqB, pCpC. Вычисляя момент инерции этой системы относительно точки AA, получим (см. (7))

  • Формула (8) верна и в случае, когда xix_i — числа: можно просто считать их векторами с координатами (xi;0;0)(x_i; 0; 0). С другой стороны, ее можно применять к случайным величинам, принимающим векторы значения.

38

Для доказательства заметим, что ZZ — это центр масс системы материальных точек aA,bB,cCaA, bB, cC (см. задачу 1). Найдем моменты инерции этой системы относительно ZZ (см. (6)):

IZ=abc2+bca2+cab2a+b+c=abc=2Rsin⁡A^⋅bc=4RS=4Rrp,I_Z = \frac{abc^2 + bca^2 + cab^2}{a+b+c} = abc = 2R \sin \hat{A} \cdot bc = 4RS = 4Rrp,

и относительно 0:

Π формуле (5)∣OZ∣2=(IO−IZ)/2p=R2−2Rr.\Pi \text{ формуле (5)} |OZ|^2 = (I_O - I_Z)/2p = R^2 - 2Rr.

Одна из красивейших теорем геометрии треугольника утверждает, что середины сторона треугольника, основания его высот и середины отрезков, соединяющих точку пересечения высот с вершинами, лежат на одной окружности — так называемой «окружности девятей точек» (см. «Квант» 1981, № 8, с. 34). Но еще большее впечатление произвольит.

Теорема Фейербаха. Окружность девятей точек касается четырех «замечательных» окружностей — вписанной и трех внешисанных.

Мы доказем эту теорему только для вписанной окружности (остальные расстриваются аналогично — см. задачу 8), причем нам будет достаточно знать, что окружность девятей точек проходит через середины A1,B1A_1, B_1 и C1C_1 сторона треугольника ABCABC.

Треугольник A1B1C1A_1B_1C_1 подобен треугольнику ABCABC с коэффициентом 1/2, поэтому радиус его описанной окружности равен R/2R/2; ее центр обозначим через FF. Условие касания этой окружности с вписанной удобно записать и доказывать в виде:

∣FZ∣2=(R2−r)2.|FZ|^2 = \left(\frac{R}{2} - r\right)^2.

Поместим в вершины ABCABC треугольника массы a,b,ca, b, c (их центром масс будет точка ZZ), разобем каждую массу на две части: a=(p−b)+(p−c),b=(p−c)+(p−a),c=(p−a)+(p−b)a = (p-b) + (p-c), b = (p-c) + (p-a), c = (p-a) + (p-b) и попарно сгруппируем эти части, как показано на рисунке 6; центр масс при этом остается на месте. Центр масс системы (p−a)B,(p−a)C(p-a)B, (p-a)C находится в точке A1A_1, системы (p−b)A,(p−b)C(p-b)A, (p-b)C — в точке B1B_1, а системы (p−c)A,(p−c)B(p-c)A, (p-c)B — в точке C1C_1, следовательно, точка ZZ — это центр масс системы (p−a)A,(p−b)B1,(p−c)C1(p-a)A, (p-b)B_1, (p-c)C_1 (массы, полученные после группировки, мы одновременно уменшили вдвое). Найдем моменты инерции этой системы относительно точек FF и ZZ:

IF=((p−a)+(p−b))+(p−c))⋅(R2)2=R2a⋅p,I_F = ((p-a) + (p-b)) + (p-c)) \cdot \left(\frac{R}{2}\right)^2 = \frac{R^2}{a} \cdot p,

IZ=1p((p−a)(p−b)(c2)2)+(p−b)(p−c)(a2)2)+(p−a)(p−b)(b2)2).I_Z = \frac{1}{p} ((p-a)(p-b)\left(\frac{c}{2}\right)^2) + (p-b)(p-c)\left(\frac{a}{2}\right)^2) + (p-a)(p-b)\left(\frac{b}{2}\right)^2).

Задачка 7. Докажите, что для любого треугольника справедливо точекство (p−a)(p−b)c2+(p−b)(p−c)a2+(p−c)(p−a)b2=abcp2−4S2(p-a)(p-b)c^2 + (p-b)(p-c)a^2 + (p-c)(p-a)b^2 = abcp^2 - 4S^2.

Пользовать токдейством этой задачи (её решение приводится в конце журнала), завершить доказательство совсем легко: по формуле (5)

∣FZ∣2=IF−IZP=R24−abc4p+S2p2=R24−abc4RpR+(Sp)2=R24−rR+r2=(R2−r)2,|FZ|^2 = \frac{I_F - I_Z}{P} = \frac{R^2}{4} - \frac{abc}{4p} + \frac{S^2}{p^2} = \frac{R^2}{4} - \frac{abc}{4Rp} R + \left(\frac{S}{p}\right)^2 = \frac{R^2}{4} - rR + r^2 = \left(\frac{R}{2} - r\right)^2,

что и требовалось.

Задачка 8. Докажите утверждение теоремы Фейербаха для вписанных окружностей. (Указание: центр вписанной окружности, касающейся (извне) стороны aa, служит центром масс системы (−a)A,bB,cC(-a)A, bB, cC. Соберите эти массы в середина сторона, пользоваться равенствами a=(p−b)+(p−c),b=p−(p−b),c=p−(p−c).a = (p-b) + (p-c), b = p - (p-b), c = p - (p-c).)

Расстояние до центра масс и вывод формулы Лагранжа

Чтобы доказать формулы (7) и (6), мы выведем формулу, выражающего расстояние от произвольной точки PP до центра масс ZZ системы материальных точек m1A1,…,mnAnm_1A_1, \dots, m_nA_n через расстояния от этой точки до точек системы: ri=∣PAi∣r_i = |PA_i|, и попарные расстояния между

(Окончание см. на с. 48)