72776

Выбор эргономически обоснованных параметров мобильного транспортного средства на основе оптимизированной модели его колебательной системы

Курсовая

Экономическая теория и математическое моделирование

Параметры пружин демпферов; длины моменты инерции функции дороги Для оптимизации было выбрано линейное ускорение массы М2 при варьировании параметров. Этот подход усложнит формулу кинетической энергии однако упростит формулы для потенциальной энергии и диссипативной функции.

Русский

2014-11-28

623.37 KB

5 чел.

Министерство образования Республики Беларусь

Белорусский Национальный Технический Университет

Факультет Информационных Технологий и Робототехники

Кафедра «Системы автоматизированного проектирования»

К У Р С О В А Я   Р А Б О Т А

По дисциплине "Компьютерные методы математического моделирования"

«Выбор эргономически обоснованных параметров мобильного транспортного средства на основе оптимизированной модели его колебательной системы»

         

                                                                       Выполнил:            студент группы 107522

                                                                                             Мартинович В. В.

                                                                          Руководитель:      Напрасников В. В.

Минск 2014

СОДЕРЖАНИЕ

Введение……………………………................…………………..………….………..………3

1. Постановка  задачи…………………………...............………....…….…………..…..……4

2. Вывод системы формул для расчёта оптимизационной модели автомобиля.........….…5

2.1. Расчётная схема многоопорной машины с указанием варьируемых параметров …...5

2.2. Вывод обобщенных перемещений………………………………………........................6

2.3. Составление выражений для удлинения и скорости удлинения упругих элементов..7

2.4. Полная кинетическая энергия системы. Ее выражение и частные производные…....8

2.5. Полная потенциальная энергия системы. Ее выражение и частные производные…..9

2.6. Диссипативная функция. Ее выражение и частные производные…...........................10

2.7. Составление системы уравнений Лагранжа 2-го рода ……………….........................11

2.8. Сведение системы ОДУ 2-го порядка к системе ОДУ 1-го порядка в канонической форме Коши..............................................................................................................................12

3. Тестирование полученной модели …………………..…………...…………...................13

3.1. Тест 1……………………………………………………………….…….........................13

3.2. Тест 2……………………………………………..……………………............................14

3.3. Тест 3……………………………………………….…….……………............................14

3.4. Тест 4………………………………………………..……….…………...........................15

3.5. Тест 5………………………………………………..……….…………...........................16

4. Оптимизация на основе параметров k3, l6……………………….......…………..............18

4.1 Разработка программы для решения системы ОДУ в канонической форме Коши средствами MathCAD, на основе предложенного алгоритма…………….……………….18

4.2 Изучение встроенной процедуры оптимизации в MathCAD …….........................…..19

4.3. Подготовка модели в виде пригодном к использованию функцией Minimize….......20

4.4. Выполнение оптимизационных вычислений………………………….........................21

4.5. Построение на одном графике ускорений верхней массы для исходных и найденных оптимизированных параметров..……....................................................................................22

5. Реализация в пакете MATLAB…………………………………………….…..................23

Заключение……………………………………………………………………...…................24

Приложение (MathCAD&MATLABFile)……………………………….….…….................25

Введение

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

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

Трудно назвать такую область техники, в которой не была бы актуальной проблема изучения упругих колебаний. Большое внимание исследователей привлечено к вопросам колебаний конструкций самых различных назначений: роторов турбин, валов двигателей внутреннего сгорания, турбинных лопаток, воздушных и гребных винтов, автомобилей и железнодорожных вагонов, кораблей, инженерных сооружений, перекрытий промышленных зданий, деталей, обрабатываемых на металлорежущих станках, вибротранспортеров и т.п. В ряде случаев колебания мешают нормальной эксплуатации или даже непосредственно угрожают прочности, постепенно подготавливая усталостное разрушение; в таких случаях теория может указать пути для уменьшения вредных колебаний. Наряду с этим она позволяет обосновать и оптимизировать технологические процессы, в которых колебания используются целенаправленно (например, в вибротранспортной технике).

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

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

3

1. ПОСТАНОВКА  ЗАДАЧИ

В курсовой работе необходимо:

  1.  разработать модель многоопорной колебательной системы (транспортного механизма);
  2.  на основе полученной модели произвести тестирование данного транспортного средства в зависимости от различных внутренних и внешних параметров, действующих на систему;
  3.  эргономически обоснованно реализовать выбор определенных параметров, влияющих на комфортабельность места водителя.

В качестве инструментариев к реализации поставленной задачи необходимо использовать следующие пакеты:

  1.  MathCAD;
  2.  MATLAB.

4

2. ВЫВОД СИСТЕМЫ ФОРМУЛ ДЛЯ РАСЧЁТА ОПТИМИЗАЦИОННОЙ МОДЕЛИ АВТОМОБИЛЯ

2.1. Расчётная схема многоопорной машины

с указанием варьируемых параметров

Кратко охарактеризуем её. Машина представляет собой корпус М1, который получает кинематическое воздействие от дороги , посредством двух опор. Эти опоры состоят из параллельно соединённой пружины и демпфера () и (). Сверху на корпусе М1 расположены массы М2, М3 и М4. Эти массы соединяется с корпусом М1 при помощи трех опор: соединённых параллельно пружины и демпфера (), () и (). Расстояние oт центра корпуса до опор () и () равны l1 и l2 соответственно. Расстояние от центра корпуса до опор () и () равно l6+l3+l4 и l4 соответственно. Расстояние от центра корпуса до опор () равно l5. Запишем исходные данные модели:

Параметры пружин, демпферов; длины, моменты инерции, функции дороги

                                                       

Для оптимизации было выбрано линейное ускорение массы М2 при варьировании параметров ,  .

5

2.2. Выбор обобщенных перемещений.

Целесообразно рассматривать относительные перемещения , , ,, , , которые отсчитываются относительно земли:

q1 — вертикальное перемещение массы корпуса относительно грунта;

q2 — угол поворота корпуса относительно горизонта;

q3 — вертикальное перемещение массы М2 относительно грунта;

q4 — угол поворота массы М2 относительно горизонта;

q5 — вертикальное перемещение массы М3 относительно грунта;

q6 — вертикальное перемещение массы М4 относительно грунта;

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

6

2.3. Составление выражений для удлинения и скорости удлинения упругих элементов

Перед началом записи выражений удлинения упругих элементов, определим следующие правило знаков: знак “-” будет браться, если упругий элемент сжимается, знак  “+“  будет браться, если упругий элемент растягивается.

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

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

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

7

2.4. Полная кинетическая энергия системы.

Ее выражение и частные производные

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

Вычислим частные производные вида:

Они понадобятся нам для составления уравнения Лагранжа 2-го рода.

8

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

Ее выражение и частные производные

Выражение для потенциальной энергии системы:

Определим частные производные вида:

которые понадобятся при составлении уравнений Лагранжа 2-го рода.

9

2.6. Диссипативная функция. Ее выражение и частные производные

Выражение для диссипативной функции:

Найдем частные производные вида:

10

2.7. Составление системы уравнений Лагранжа 2-го рода

В систему уравнений Лагранжа второго рода (количество уравнений равно количеству неизвестных, т.е. количество уравнений равно шести) будут входить уравнения типа:

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

Система Лагранжа 2-го рода:

11

2.8. Сведение системы ОДУ 2-го порядка к системе ОДУ 1-го порядка в канонической форме Коши

Для того чтобы свести систему ОДУ Лагранжа 2-го рода к системе ОДУ 1-го рода в канонической форме Коши, введём следующие новые переменные.

Заменив, получим:

Система значений вида  является системой, пригодной для решения пакетами MathCAD и MATLAB.

12

3. ТЕСТИРОВАНИЕ ПОЛУЧЕННЫХ УРАВНЕНИЙ

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

Тест № 1: принимаем функции дороги , равными нулю.

13

Результат – в отсутствии кинематического воздействия дороги система не колеблется, т.е. все перемещения (углы поворота) и скорость их изменения равны нулю.

Тест № 2: задаем очень большую жесткость для пружин, соединяющих массы М1 с М2 и М1 с М3 ) и задаем линейное перемещение массы М1  q1=2.

  

Результат: массы оказались «связанными», их колебания почти идентичны.

Тест № 3: разрываем связь между корпусом и массой М2 (k3=c3=0).

14

Результат – т.к. масса М2 получает кинематическое воздействие дороги через корпус, то разорвав связь между корпусом и массой М2, мы получаем отсутствие колебаний (поворотов) массы М2.

Тест № 4: добавляем демпферы.

15

Результат – наблюдаем угасание колебаний.

Тест № 5: разрываем связь между корпусом и массой М2 и М4 (k3=c3=k5=c5=0), принимаем функции дороги , .

16

Результат – т.к. масса М3 находится посередине рамы и получает кинематическое воздействие дороги через корпус, то образовав связь между корпусом и массой М3, мы получаем только наличие поворотов рамы.

17

4. ОПТИМИЗАЦИЯ НА ОСНОВЕ ПАРАМЕТРОВ K3, L6

4.1. Разработка программы для решения системы ОДУ в канонической форме Коши средствами MathCAD, на основе предложенного алгоритма

Для решения системы ОДУ в каноническом виде Коши используем метод Рунге-Кутта 4-го порядка. Приращение переменной можно выразить по следующей формуле:

,

где:

Тогда программная реализация данного метода будет иметь вид:

В результате вычислений и всех тестов, получается полное совпадение результатов с результатами, полученными при помощи встроенной функции Rkadapt (Приложение 1).

18

4.2.Изучение встроенной процедуры оптимизации в MathCAD

Являясь мощным пакетом для различных математических вычислений, MathCAD включает в себя множество функций. Одной из таких функций является функция Minimize.

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

Перед началом блока оптимизации пишется ключевое слово Given, которое позволяет определить пакету MathCAD, где начинается блок минимизации. После ключевого слова Given перечисляются переменные с указанием граничных условий. По этим переменным MathCAD и будет искать минимальное значение функции. Например, следующая запись () означает, что MathCAD будет искать минимальное значение функции, изменяя параметр . После перечисления переменных следует вызов функции Minimize. Синтаксис у функции следующий:

Minimize(<имя_функции>,<параметр1>,<…>,<параметрn>),

где - <имя_функции> - имя той функции, минимальное значение которой необходимо найти, <параметр1>…<параметрn> - список параметров, по которым функция Minimize будет искать минимум.

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

Так же необходимо отметить, что чем больше входных параметров у функции Minimize, тем большее количество времени занимает поиск минимального значения. Отсюда следует вывод – количество входных параметров у функции Minimize (как было замечено выше, и у минимизированной функции тоже) должно быть минимально.

Функция Minimize возвращает матрицу-столбец с перечислением оптимальных параметров.

19

4.3. Подготовка модели в виде пригодном к использованию функцией Minimize

Однако, мы не можем использовать для минимизации непосредственно функцию вычисления ускорения, т.к. значение ускорения многократно уходит в отрицательную область, а значения функции должны быть неотрицательными. Поэтому для минимизации будем использовать функцию нахождения среднего квадрата ускорения. Эта функция имеет два входных параметра k3 и l6. Возвращает эта функция одно число – средний квадрат ускорений. (Приложение 2).

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

20

4.4. Выполнение оптимизационных вычислений.

Опишем блок минимизации:

где, Given – ключевое слово, указывающее пакету MathCAD начало блока минимизации, 0<=k3<=5, 0<=l6<=3 - ограничение на изменение входных параметров, Opt:=Minimize(USKOR,k3,l6) – вызов функции минимизации, USKOR – имя функции, используемой для вычисления среднего квадрата ускорения (разброса ускорений).

Просчитав возможные варианты и выбрав тот, при котором разброс ускорений имеет минимальное значение, функция Minimize вернёт нам следующую матрицу-столбец:

21

4.5. Построение на одном графике вертикального перемещения верхней массы для исходных и найденных оптимизированных параметров.

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

Как видно из рисунка, хотя оптимизированное перемещение меньше исходного,-  оптимизация была произведена успешно.

22

5. Реализация в пакете MATLAB

Все приведенные в данной работе тесты были также реализованы в среде MATLAB и имели успех. Итоговое значение (графики) реализации в двух различных средах оказалось полностью идентичным. Некоторые результаты тестов приведены ниже (в сравнении с пакетом MathCAD). Саму реализацию можно увидеть в Приложении 3 данной работы.

Тест № 1:                                                       Тест № 2:

                   

Тест № 3:                                                       Тест № 4:

            

Тест № 5:

23

ЗАКЛЮЧЕНИЕ

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

  1.  Получены уравнения, описывающие колебания многоопорной машины при кинематических воздействиях со стороны дороги;
  2.  Выполнено тестирование полученных уравнений на многочисленных примерах;
  3.  Построена оптимизационная модель, ориентированная на среду MATHCAD;
  4.  Выполнены оптимизационные вычисления;
  5.  Получили структурную схему для решения системы ОДУ в среде  MATLAB.

24

ПРИЛОЖЕНИЕ 1

25

ПРИЛОЖЕНИЕ 2

26

ПРИЛОЖЕНИЕ 3

27


 

А также другие работы, которые могут Вас заинтересовать

37103. Коллективизация сельского хозяйства и ее последствия 19.56 KB
  Коллективиза́ция это процесс объединения единоличных крестьянских хозяйств в коллективные хозяйства колхозы в СССР. объединению подлежало 25 процентов крестьянских хозяйств. Основу аграрного сектора экономики страны составляли мелкие крестьянские хозяйства которые имели полунатуральный характер.
37104. Начальный период Великой Отечественной войны. Причины поражения Красной Армии на первом этапе войны 15.43 KB
  Причины поражения Красной Армии на первом этапе войны. 22 июня 1941 года нацистская Германия вероломно нарушив договор о ненападении внезапно без объявления войны нанесла по Советскому Союзу мощный удар. Этот день вошел в историю нашей страны трагической датой стал днем начала неимоверно тяжелой войны советского народа против фашизма справедливо названной Великой Отечественной войной.
37105. Причины, цена и историческое значение победы советского народа в Великой Отечественной войне 17.16 KB
  Решения партии лежали в основе постановлений Государственного Комитета Обороны Президиума Верховного Совета СССР и Совета Народных Комиссаров СССР приказов Народного комиссара обороны и Верховного Главнокомандующего а также призывов и лозунгов обращенных к советскому народу и воинам Вооруженных Сил. В руках Государственного Комитета Обороны сосредоточилась вся полнота власти в государстве. Все граждане и все партийные советские комсомольские и военные органы обязаны были беспрекословно выполнять решения и распоряжения Государственного...
37106. Послевоенное устройство мира. Формирование двух мировых систем и их противостояние 18.75 KB
  Особенностью послевоенного периода стали народнодемократические социалистические революции в странах Восточной Европы и ряде стран Азии которые при поддержке СССР приступили к строительству социализма. Образовалась мировая система социализма во главе с СССР. В годы мировой войны сложилась единая антифашистская коалиция – союз великих держав – СССР США Великобритания и Франция. Основой конфронтации стали отношения между двумя сверхдержавами – СССР и США.
37107. СССР в период восстановления разрушенного войной хозяйства 23.22 KB
  В 1943 году для организации восстановительных работ в Совете Народных Комиссаров СССР был образован специальный комитет а в Госплане СССР управление обеспечивающее планирование деятельности всех ведомств причастных к этой работе. Постановление ЦК ВКП б и СНК СССР от 21 августа 1943 г. В Госплане СССР и научных учреждениях были разработаны краткосрочные и долгосрочные планы восстановления народного хозяйства как в целом так и по отдельным отраслям. На сессии Верховного Совета СССР был утвержден пятилетний план восстановления и развития...
37108. Разоблачение культа личности Сталина. XX съезд КПСС 15.81 KB
  Именно Хрущеву было поручено прочитать доклад и лично встретить непредсказуемую реакцию участников съезда. Созванный за восемь месяцев до уставного срока в связи с насущной необходимостью подвести итоги произошедших после смерти Сталина изменений и дискуссии о выборе курса съезд завершился знаменитым секретным докладом Хрущева. в последний день работы XX съезда на закрытом заседании с докладом О культе личности и его последствиях выступил первый секретарь ЦК КПСС Н. В докладе были обнародованы и осуждены факты массовых репрессий...
37109. Политика либерализации советского общества. «Оттепель». Н.С. Хрущев – политик и человек 18.45 KB
  Хрущев – политик и человек. Хрущев. После недолгой борьбы к власти пришел Никита Сергеевич Хрущев. Резкий решительный неосторожный в словах и поступках Хрущев прошел все ступени партийной работы возглавлял крупные парторганизации.
37110. Нарастание застойных явлений в жизни СССР в 70-е – начале 80-х гг. и их сущность 27.5 KB
  Сохранилась жёстко централизованная система власти. Костяк этой власти составили аппаратчики. К власти пришло третье поколение советских правителей типичных аппаратчиков исполнителей. В моменты ослабления центрального руководства незаметные исполнители почувствовали мощь власти и те блага которые можно приобрести с её помощью.
37111. Перестройка политической системы СССР 24.95 KB
  В истории СССР начался этап получивший название перестройка . По решению XIX партконференции учреждается новый высший орган законодательной власти Съезд народных депутатов СССР и соответствующие республиканские съезды. Был введен пост Президента СССР первым и последним Президентом СССР стал М.