Расчет максимального значения восстанавливающей силы

Сибирский государственный университет путей сообщения

Домашнее задание по дисциплине «Математическое моделирование»

Задачи №1, №2

Разработал: студент гр. М-511

Ревнивцев

2008

Задача№1

В тупике железнодорожного пути установлен буфер (рисунок 1), имеющий упругий элемент с нелинейной жесткостью, восстанавливающая сила которого изменяется по закону

>>.

1 – вагон; 2 – буфер; 3 – демпфер

Рисунок 1 – Схема к решения задачи.

В направлении тупика движется вагон массой m и скоростью X. При столкновении вагона с упругим элементом, последний смещается на величину X1. В задаче также приняты следующие допущения: 1) масса буфера мала по сравнению с массой вагона; 2) после удара контакт между этими объектами сохраняется.

Восстанавливающая сила (закон изменения):

>> (1)

В задаче требуется определить:

- максимальное перемещение буфера;

- максимальное значение восстанавливающей силы;

- время, за которое восстанавливающая сила достигнет максимального значения.

На данную систему (рисунок 1) действуют силы: сила инерции движущегося вагона; сила демпфирования (или сила вязкого трения), пропорциональная скорости движения вагона; а также восстанавливающая сила упругого элемента-демпфера.

Сила инерции:

>> , (2)

Сила демпфирования:

>>, (3)

Сила упругости:

>> (4)

Для решения поставленной задачи следует решить обыкновенное дифференциальное уравнение второго порядка вида:

>> (5)

Заменим уравнение (5) системой уравнений первого порядка, для этого введем новую неизвестную функцию > > и перепишем исходное уравнение, представив его в виде системы из двух уравнений:

>> (6)

Решение проводим в системе MathCad с построением графических зависимостей: 1) скорости движения вагона от времени; 2) перемещения буфера от времени; 3) восстанавливающей силы от времени. Неизвестные выше исходные данные записываются непосредственно в программе.

Исходное уравнение имеет вид:

>> , (7)

где > > - коэффициент демпфирования, > >;

>> - масса вагона, кг;

>>- жесткость упругого элемента, Н/м;

>> - численный коэффициент, > >;

>> - скорость вагона при подходе к тупику, м/с;

Начальные условия:

> (8)>

Уравнение (8) решается в системе MathCad посредствам встроенной функцией – rkfixed:

Z:= rkfixed (y, 0, t, n. D) (9)

где Z – вектор неизвестных;

y – вектор начальных условий;

0 и t – интервал, на котором ищется решение;

n – количество точек на интервале.

>>


В ходе вычислений получена следующая система ответов:

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100fcaa0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004b1200005923000020454d4600000100180300001200000002000000000000000000000000000000d0120000101a0000cc0000001a010000000000000000000000000000651c0300504f0400160000000c000000180000000a0000001000000000000000000000000900000010000000520400005a080000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001100c09a110010000000249e1100a49b1100e4506032249e11001c9b1100100000008c9c1100089e11008a4f6032249e11001c9b11002000000080e02d311c9b1100249e110020000000ffffffff9c2fdf00fbe02d3127000300ffff0180ffff01800f020180ffffffff00000000000800000008000094e7120001000000000000005802000025000000372e9001cc00020f0502020204030204ef0200a07b20004000000000000000009f00000000000000430061006c006900620072000000000000000000e49b1100334e2f31282e0c32449f1100509b1100d532273106000000010000008c9b11008c9b11000076253106000000b49b11009c2fdf006476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000558c8740764087400000000057000000010000004c000000040000000000000000000000520400005a08000050000000200000003600000046000000280000001c0000004744494302000000ffffffffffffffff530400005a080000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c025901b200040000002e0118001c000000fb020300010000000000bc02000000cc0102022253797374656d0000000000000000000000000000000000000000000000000000040000002d010000040000002d01000004000000020101001c000000fb02f1ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0e0000000100040000000000b200580120840900040000002d010000040000002d010000030000000000

Колонка «0» - промежутки времени; колонка «1» - перемещение в каждый момент времени; колонка «2» - скорость вагона в каждый момент времени.

Определение восстанавливающей силы.

Расчет ведется при разбиении - > >

Функциональная зависимость в программе:

>>

0100090000032a0200000200a20100000000a201000026060f003a03574d46430100000000000100fcaa0000000001000000180300000000000018030000010000006c0000000000000000000000350000006f00000000000000000000004b1200005923000020454d4600000100180300001200000002000000000000000000000000000000d0120000101a0000cc0000001a010000000000000000000000000000651c0300504f0400160000000c000000180000000a0000001000000000000000000000000900000010000000520400005a080000250000000c0000000e000080250000000c0000000e000080120000000c00000001000000520000007001000001000000a4ffffff00000000000000000000000090010000000000cc04400022430061006c006900620072006900000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001100c09a110010000000249e1100a49b1100e4506032249e11001c9b1100100000008c9c1100089e11008a4f6032249e11001c9b11002000000080e02d311c9b1100249e110020000000ffffffff9c2fdf00fbe02d3127000300ffff0180ffff01800f020180ffffffff00000000000800000008000094e7120001000000000000005802000025000000372e9001cc00020f0502020204030204ef0200a07b20004000000000000000009f00000000000000430061006c006900620072000000000000000000e49b1100334e2f31282e0c32449f1100509b1100d532273106000000010000008c9b11008c9b11000076253106000000b49b11009c2fdf006476000800000000250000000c00000001000000250000000c00000001000000250000000c00000001000000180000000c0000000000000254000000540000000000000000000000350000006f00000001000000558c8740764087400000000057000000010000004c000000040000000000000000000000520400005a08000050000000200000003600000046000000280000001c0000004744494302000000ffffffffffffffff530400005a080000000000004600000014000000080000004744494303000000250000000c0000000e000080250000000c0000000e0000800e000000140000000000000010000000140000000400000003010800050000000b0200000000050000000c025901b200040000002e0118001c000000fb020300010000000000bc02000000cc0102022253797374656d0000000000000000000000000000000000000000000000000000040000002d010000040000002d01000004000000020101001c000000fb02f1ff0000000000009001000000cc0440002243616c6962726900000000000000000000000000000000000000000000000000040000002d010100040000002d010100040000002d010100050000000902000000020d000000320a0e0000000100040000000000b200580120840900040000002d010000040000002d010000030000000000

Колонка «0» - значения восстанавливающей силы.

Далее в программе ведется построение необходимых графиков на интервале: > >.

График зависимости перемещения от времени

>>

График зависимости скорости от времени

>>

График зависимости силы от перемещения

>>

После проведения решений выведены лишь шестнадцать рассчитанных значений. В общем же получено 300 значений, что соответствую числу интервалов.

Проанализировав результаты получаем:

- максимальное перемещение буфера – 0,0782 метра;

- максимальное значение восстанавливающей силы - > >;

- время, за которое восстанавливающая сила достигнет максимального значения – 0,725 сек.

Задача №2

Данный планетарный редуктор (рисунок 2) представляет собой механическую вращательную систему, которая состоит из четырех подсистем. Связь между подсистемами осуществляется через зубчатое зацепление.

Опишем каждую из подсистем. Первая подсистема включает зубчатое колесо 1 (рисунок 2), которое находиться во внешнем зацеплении с сателлитом 2, вал «а» вращающийся в опоре 7 от привода 5, имеющий крутящий момент М>вх>. Вторая подсистема включает: сателлит 2, который имеет внешнее зацепление с подвижным зубчатым колесом 1 и внутреннее зацепление с неподвижным зубчатым колесом 3. Третья подсистема это неподвижное зубчатое колесо 3, которое находиться во внутреннем зацеплении с сателлитом 2. Четвертая подсистема состоит из сателлита 2, вращающегося вокруг оси вала «в» вместе с водилом 4, вала «б», вала «в», опор 8,9 и нагрузки 6.

Рисунок 2 – Схема планетарног редуктор 1- подвижное центральное колесо; 2 - сателлит; 3 - неподвижное центральное колесо; 4 - водило; 5 - привод; 6 - нагрузка; 7, 8, 9 - опоры валов; а, б, в – валы.

Исходные данные:

>> - радиусы делительных окружностей зубчатых колес 1-3;

>>- моменты инерции зубчатых колес 1-3 и водила относительно осей вращения;

>> - моменты инерции привода и нагрузки;

>> - моменты вращения сателлита 2 при вращении вместе с водилом;

>> - коэффициенты вязкого трения в опорах 7-9.

Потерями в опоре 9, массой вала «а» и податливостью вала «в» пренебречь.

Сателлит 2 совершает сложное движение: его абсолютное движение складывается из относительного – вокруг собственной оси и переносного – вращение вместе с водилом относительно оси вала «в». Пользуясь методом Виллиса, определим направление переносного движения, т.е направление вращения водило. Допустим, что все звенья передачи 1,2,3 и 4 жестко соединены друг с другом. Сообщая этой жесткой системе переносное вращательное движение вокруг оси вала «в» с угловой скоростью > > равной скорости вращения водила, но обратно по знаку. При таком движении водило окажется остановленным > >, в результате относительные угловые скорости зубчатых колес 1 и 3 будут равны:

>>(1)

>>(2)

где > >абсолютные угловые скорости колес 1 и 3; > >относительные угловые скорости колес 1 и 3; > > скорость вращения водила.

При > > планетарная передача превращается в простую зубчатую передачу в которой оси всех зубчатых колес неподвижны.

Тогда передаточное отношение передачи будет иметь вид:

>> (3)

Преобразуем данное уравнение:

>>

>>

А поскольку колесо 3 является неподвижным > >, то

>>(4)

Из (4) следует, что направление вращения водила совпадает с направлением вращения колеса 1, то есть скорости относительного и переносного вращения сателлита 2 противоположны по знаку.

Используя метод аналогий представим механическую систему в виде эквивалентных схем каждой из подсистемы (рисунок 3).

Запишем уравнения для источника > >:

>>,

где > >;

Запишем уравнение равновесия моментов:

>>;

>>,

Аналогично запишем уравнение для > >, > >, > >. Так,

>>,

где > >;

Запишем уравнение равновесия моментов:

>>;

>> (5)

Уравнение для > >:

>>, где > >;

>>;

>>. (6)

Уравнение для > >:

>> , где > >;

>>;

>>.

Запишем уравнения для зависимых источников угловых скоростей:

>>,

где > >;

>>=>> (7)

>>,

где > >;

>>=>>

Составим уравнение равновесия в узлах эквивалентных схем.

1 подсистема:

узел 8: > > (8)

узел 1: > > (9)

2 подсистема:

узел 7: > > (10)

узел 2: > > (11)

узел 8: > > (12)

3 подсистема:

узел 3: > > (13)

4 подсистема:

узел 9: > > (14)

узел 10: > > (15)

узел 4: > > (16)

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

>> (17)

>> или > > (18)

>> (19)

Получим следующие уравнения:

1 подсистема:

узел 8: > > (20)

узел 1: > > (21)

2 подсистема:

узел 7: > > (22)

узел 2: > > (23)

узел 8: > > (24)

3 подсистема:

Узел > > (25)

4 подсистема:

узел 9: > > (26)

узел 10: > > (27)

узел 4: > > (28)