Колебания пусковой установки


                              Схема установки:


                                                      Рис.1

Задание на проект:

Пусковая установка находится на корабле, совершающем колебания (угол [pic]
- стационарная функция известного вида.)
В момент времени t = tк производится пуск ракеты.

Требуется:
1. Получить уравнение малых колебаний ракеты с направляющей с учетом
   воздействия со стороны корабля.
2. Определить закон изменения момента управляющего двигателя Мупр(t),
   обеспечивающего минимум среднего значения угловой скорости пусковой
   установки к заданному моменту времени t = tк. Мощность двигателя
   ограничена ( | Мупр.|[pic])



                              Расчетная схема:


                                                           Рис.2

Где точка А считается центром масс платформы с ракетой.
[pic] и [pic] - кинематическое возбуждение точек основания
[pic]- угол подъема платформы в стационарном состоянии
[pic] - приращение угла (считается малым)


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



[pic]                   [pic]      [pic]
                                 [pic]



                                                     [pic]
       [pic]
                                                                    [pic]

                                    Рис.3



              Где [pic], [pic] или с учетом малости воздействия
                                                   [pic], [pic]

                 Тогда возмущающие функции будут иметь вид:
                     [pic]                          (1)
                      [pic]                         (2)


Кинетическая энергия системы:

                          [pic]                (3)
[pic]- абсолютная скорость центра масс платформы,
[pic]- момент инерции платформы с ракетой, относительно центра масс.

 По теореме косинусов:       [pic]      (4),                            где
                                    [pic]

Таким образом, кинетическая энергия системы запишется в виде:

                      [pic]                         (5)

Потенциальная энергия системы:

Поскольку перемещения системы считаются малыми, а пружина обладает
достаточной жесткостью, потенциальной энергией силы тяжести пренебрегаем.
То есть потенциальная энергия системы будет потенциальной энергией,
накопленной в пружине.

                          [pic]                (6)

С учетом (1) и (2) получаем:

                            [pic]             (7)

Для записи уравнения движения воспользуемся уравнением    Лагранжа:

                            [pic]            (8)

                     [pic]                           (9)

                         [pic]                  (10)

                          Учитывая, что [pic] получим:

            [pic][pic]                                       (11)
                              [pic]        (12)


Подставляя (11) и (12) в уравнение Лагранжа, получим следующее:

                       [pic]                     (13)

Уравнение движения будет иметь вид:

[pic]                                      (14)

Или, с учетом управляющего момента:

      [pic]                    (15)

Считаем, что на систему действуют функция: [pic]
где  А –амплитуда, а [pic]-частота вынуждающих функций.

Уравнение движения можно переписать в виде:

            [pic]                       (16)

                                 где   [pic]
Решение этого дифференциального уравнения состоит из двух частей:
1. Решение однородного дифференциального равнения
2. Частное решение неоднородного уравнения

Решение однородного уравнения имеет вид:

                            [pic][pic]      (17)

Частное решение неоднородного уравнения при произвольном воздействии будет
выглядеть так:

    [pic]         (18)


          Тогда общее решение дифференциального уравнения:


                                 [pic]  (19)

Выражение для скорости:


[pic](20)



Компенсирующий двигатель включается в момент времени [pic] .
Он работает до момента времени  [pic] . Мощность двигателя – ограничена.
Интегрирование начинаем в момент времени [pic], но т.к. [pic] функция
известного вида, а начальный момент времени - произвольный, то не важно, с
какого момента начинать интегрирование, поэтому, начальный момент времени
принимаем
нулевым. Исходя из подобных соображений, начальные условия так же считаем
нулевыми, т.е.          [pic]

              Таким образом, приходим к выражению для скорости:

                               [pic]     (21)

В момент пуска ракеты угловая скорость вращения платформы должна быть
минимальной, в идеале – нулевой, поэтому:

               [pic]     (22)

Если добиться нулевого значения угловой скорости не представляется
возможным, то потребуем нахождения угловой скорости в заданных пределах
[pic]

Идеология решения такой задачи такова: Разобьем подинтегральное выражение
на два интеграла. Тогда выражение для скорости запишется в следующем виде:

[pic]
                                     (23)

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


[pic]
               Функция управляющего момента будет иметь такой вид:


                                           [pic]       (23)

                                                      где   [pic]



         Область, ограничивающая управляющий момент:


                      [pic]



                   [pic]                                    [pic]
      [pic]


                             Рис 4.


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



Задаемся следующими параметрами установки:

[pic]


Тогда остальные параметры будут вычисляться по формулам:

              [pic]      [pic]

                                     [pic]

Амплитудное значение возмущающей функции: [pic]

Рассмотрим наиболее «неприятный» случай – когда частота возмущающей функции
совпадает с собственной частотой системы, т.е. [pic]

График возмущающей функции   [pic]:


                             Рис. 5



Момент времени пуска ракеты [pic]
График управляющего момента:



                                  Рис.6



Обозначим функцию [pic]
График этой функции:



                             Рис. 7



                График функции [pic]



Рис.8



    Представим функции [pic][pic] на одном графике:



                             Рис 8.



Тогда для амплитудного значения управляющего момента [pic]        [pic],
времени пуска ракеты равным 15 с, для заданных параметров установки,
значение угловой скорости равно:[pic]
Амплитудное значение управляющего момента было получено путем подбора, при
этом выбиралось минимальное значение угловой скорости.



Для сравнения, выберем управляющий момент двигателя в виде произведения
амплитудного значения момента на косинус с обратным знаком, т.е. :
[pic]

Вид функции управляющего момента представлен на рисунке:



Рис.9


С учетом наложения [pic]:



Рис.10


Тогда для получения аналогичного значения угловой скорости, а именно для
[pic] потребуется амплитудное значение управляющего момента, большее в два
раза, чем при предыдущем случае, а именно  [pic].

Таким образом, обеспечив вид функции управляющего момента двигателя, как
         [pic]
можно добиться минимального значения угловой скорости к моменту времени
пуска ракеты, используя двигатель мощностью меньшей, чем в случае чисто
косинусоидального вида управляющего момента.

По найденному амплитудному значению управляющего момента можно подобрать
двигатель требуемой мощности, закон движения которого будет изменяться
согласно предъявляемым выше требованиям.
-----------------------
[pic]



[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]

[pic]