ДРУГОЙ ВАРИАНТ ВЫКЛАДОК (лекции 2011).
Сумма двух последних слагаемых дает
n
- S l(t0j rj) (Fвj (0) + Fиj(0)) .
j=i
В самом деле:
n n n-1 n n
- S (j) S (k) + S (j) S (k) = - S (j) (k=j).
j=i k=j j=i k=j+1 j=i
j = i k = i,i+1,…,n j = i k = i+1,i+2,…,n
j = i +1 k = i+1,…,n j = i+1 k = i+2,…,n
…
j = n -1 k = n-1,n j = n-1 k = n
j = n k = n
После вычитания
j = i k = i= j
j = i +1 k = i+1 = j
…
j = n -1 k = n-1 = j
j = n k = n = j.
С учетом этого выражение для mi (0) примет вид:
n n n
mi (0) =- S (Mвj (0) + Mиj(0) ) - S l(t0j (vj (1-sj)qj + lj)) S (Fвk (0) + Fиk(0)) -
j=i j=i k=j
n
S l(t0j rj) (Fвj (0) + Fиj(0)).
j=i
i = 1,2,…,n.
Выражение для Mиj(0) имеет вид:
Mиj(0) = - (Ij w`j)(0) - l(wj(0)) (Ij wj)(0).
В последнем соотношении обозначено (Ijwj)(0) – вектор момента количества движения звена j в СК0. Его можно представить еще и так:
(Ij wj)(0) = t0j (Ij tj0 wj(0) ) = (t0j Ij tj0)wj(0) = Ij(0) wj(0),
где Ij(0) = t0j Ij tj0 - тензор инерции звена j в системе коодринат, оси которой параллельны осям СК0, а начало совмещено с центром масс звена j. Аналогично
(Ij w`j)(0) = Ij(0) w`j(0).
Справедливость выражения
Ij(0) = t0j Ij tj0
Следует из того, что
I(0) j = Sn i l(r (0)n i) m ni lT(r (0)n i) = Sn i l(t0j r n i) m ni lT(t0j r n i) =
Sn i t0j l(r n i) tj0 m ni t0j lT(r n i) tj0 = t0j {Sn i l(r n i) m ni lT(r n i)} tj0 = t0j Ij tj0 .
|
|
Диагональные элементы матрицы Ii - моменты инерции звена i относительно осей СКiс; недиагональные элементы – смешанные моменты инерции.
Подставив соотношение для Mиj(0) в формулу для расчета mi (0), получим:
n
mi (0) = S (Ij(0) w`j(0) + l(wj(0)) Ij(0) wj(0) - Mвj (0) ) +
j=i
n n
S l(t0j (vj (1-sj)qj + lj) S (mk V`k (0) - Fвk (0)) +
j =i k=j
n
S l(t0j rj) (mj V`j (0) - Fвj (0) ),
j=i
i = 1,2,…,n.
Последнее выражение можно преобразовать. Рассмотрим отдельно второе и третье слагаемые в этом выражении. Во втором слагаемом переменим обозначения (j<->k) для индексов суммирования, а затем изменим порядок суммирования:
n n n
S l(t0k (vk (1-sk)qk + lk) S (mj V`j (0) - Fвj (0)) + S l(t0j rj) (mj V`j (0) - Fвj (0) ) =
k=i j=k j =i
k = i j = i, i+1,…,n
k = i+1 j = i+1,…,n -> j = i,…,n; k = i,…,j.
…
k = n j = n
n j n
S S l(t0k (vk (1-sk)qk + lk) (mj V`j (0) - Fвj (0)) + S l(t0j rj) (mj V`j (0) - Fвj (0) ) =
j=i k=i j =i
n j j
S l (St0k (vk (1-sk)qk + lk) + t0j rj) (mj V`j (0) - Fвj (0)) = S l(Rji(0))(mj V`j (0) - Fвj (0)).
j =i k=i j =i
Здесь мы воспользовались тем, что
j
Rji(0) = t0jtj + S t0k (vk (1-sk)qk + lk)
k=i
и тем, что tj = rj.
С учетом выполненных преобразований выражение для расчета моментов реакции связей примет вид:
n
mi (0) = S (Ij(0) w`j(0) + l(wj(0)) Ij(0) wj(0) - Mвj (0) + l(Rji(0))(mj V`j (0) - Fвj (0))).
j=i
i = 1,2,…,n.
/*
Рассмотрим участок кинематической цепи, начиная от звена i. Согласно принципу Даламбера этот участок кинемсатической цепт находтся в равновесии под действием всех внешних сил и всех сил инерции, действующих на звенья от i до n и сил реакции связей действующих на звено i, что эквивалентно записи
n n
mi (0) + S (Mиj (0) + Mвj (0) ) + Sl(Rji(0))(Fиj + Fвj (0)) = 0,
j=i j=i
где
j
Rji(0) = t0jtj + S t0k (vk (1-sk)qk + lk)
k=i
Учитывая¸что
Mиj(0) = - Ij (0) w`j (0) - l(wj(0)) Ij(0) wj(0) ,
Fиj(0) = - mj V` (0)j
|
|
получим
n
mi (0) = S (Ij(0) w`j(0) + l(wj(0)) Ij(0) wj(0) - Mвj (0) + l(Rji(0))(mj V`j (0) - Fвj (0))).
j=i
*/
Геометрический смысл векторов Rji(0) был рассмотрен выше.
Последние выражения совместно с выражениями для расчета wj(0), w`j(0) позволяют определить моменты, действующие на звенья ИМ при его движении. Можно видеть, что к звену i приложены все внешние моменты, действующие на звенья с номерами от i до n, а также моменты сил инерции, обусловленные движением всех звеньев ИМ. На первое звено ИМ действуют все внешние моменты и моменты сил инерции всех звеньев. Момент реакции со стороны основания ИМ (действует на первое звено) равен m1(0). ИМ воздействует на основание с моментом (-m1(0)).
Ось шарнира с номером i совпадает с осью Zi-1. Поэтому для определения сил, действующих вдоль осей шарниров следует, во-первых, векторы сил реакций связей перевести в системы из СК0 в СК i-1 , а во-вторых, спроецировать получившиеся векторы на оси Zi-1:
n n
f Zi(i-1) = eT ti-1i ti0 S (mj V` (0)j - Fвj (0) ) = viT ti0 S (mj V` (0)j - Fвj (0) ),
j =i j =i
поскольку eT ti-1i = viT, (i = 1,2,…,n).