- •1. Определение математической модели и математического моделирования
- •2. Основные этапы математического моделирования
- •3. Свойства математических моделей
- •4 Требования к математическим моделям
- •5 Классификация моделей
- •6. Иерархия мм и формы представления
- •7. Краевые задачи проектирования
- •9. Мм на микроуровне
- •13. Методика получения функциональных моделей
- •14. Метод получения топологических уравнений
- •15. Метод конечных элементов
- •16. Метод конечных разностей
- •17. Метод граничных элементов
- •18. Аналогии компонентных уравнений
- •19. Аналогии топологических уравнений
- •20. Получение эквивалентных схем технических объектов.
- •21. Аппроксимация табличных данных. Метод наименьших квадратов.
- •23. Метод Ритца-Галеркина
- •25. Табличный метод получения математических моделей систем
- •26. Узловой метод получения математических моделей систем.
- •28. Метод вращения Якоби
- •29. Методы решения систем линейных и нелинейных алгебраических уравнений.
- •30. Анализ в частотной области.
- •31. Сравнение методов конечныx элементов и конечных разностей
- •33. Математические модели дискретных устройств.
- •34. Многовариантный анализ.
- •35. Основные сведения из теории массового обслуживания
- •36. Имитационное моделирование смо
- •38. Геометрические модели
- •39. Методы и алгоритмы машинной графики.
7. Краевые задачи проектирования
Проектирование многих технических объектов связано с необходимостью анализа непрерывных физических процессов, математическим описанием которых являются дифференциальные уравнения в частных производных. Примером тому служат современные летательные аппараты, при проектировании и расчете которых широко используется анализ подобных моделей.
Примеры уравнений, составляющих основу моделей объектов на микроуровне. Первая важная задача проектирования летательного аппарата — определение прочности узлов и элементов конструкции при различных видах нагружения. Поэтому исследование напряженного состояния деталей конструкции и связанные с ним расчеты на прочность относятся к наиболее ответственным в самолетостроении.
Напряженное состояние деталей конструкции в зависимости от геометрии исследуемого узла, вида приложенной нагрузки и свойств материала описывается дифференциальными уравнениями различного вида. Любое из этих уравнений может быть получено из общего квазигармонического уравнения
(2.1.) где х, у, z — пространственные координаты; φ — искомая непрерывная функция; Кх, Ку,Kz — коэффициенты; Q— внешнее воздействие. В двухмерном случае при Кх = Ку=1 уравнение (2.1) сводится к уравнению, описывающему напряженное состояние, возникающее в поперечном сечении упругого однородного стержня под воздействием крутящего момента М:
д2φ/дх2 + д2φ/ду2 + 2Еθ = 0, (2.2)
где Е — модуль сдвига материала стержня; θ — угол закручивания на единицу длины, а функция φ — функция, связанная с напряжениями сдвига τх и τу уравнениями
τх = дφ/ду; τу = дφ/дх. (2.3)
В (1.2) в явном виде не входит крутящий момент, связанный с искомой функцией напряжения φ уравнением
М = 2 ,
где S — площадь рассматриваемого сечения.
Вторая важная задача проектирования летательного аппарата — изучение его аэродинамических свойств. Решение этой задачи связано с исследованием процессов обтекания газом поверхностей произвольной формы. Наиболее общими уравнениями, описывающими этот процесс, являются уравнения Навье — Стокса, которые в декартовой системе координат имеют вид
(2.4)
где и, v, ω— проекции вектора скорости; Fх , Fy ,Fz — проекции вектора силы на оси координат; ρ — плотность; р — давление; v = μ/ ρ(коэффициент вязкости).
Третья задача проектирования летательных аппаратов — расчет тепловых режимов работы деталей и узлов конструкции. Одним из основных аспектов задачи является определение температурных полей, имеющих место в конструкциях.
Примечание. Знание температурных полей необходимо для вычисления количества теплоты, подводимой к телу или отводимой от него. Кроме того, температурные поля влияют на распределение напряжений в конструкциях. Это обстоятельство особенно важно учитывать при проектировании вращающихся элементов летательных аппаратов.
Температурное поле в сплошной среде описывается уравнением теплопроводности. Последнее может быть получено из уравнения (1.1), если под функцией φ понимать температуру Т, а под коэффициентом К — коэффициент теплопроводности λ. Так, в двухмерном случае при условии, что коэффициенты теплопроводности λх и λу по соответствующим направлениям не зависят от координат, стационарное уравнение теплопроводности имеет вид
(2.5)
где Q — источник теплоты внутри тела, который считается положительным, если теплота подводится к телу. Сформулированные выше задачи — типичные для многих областей техники. Так, задачу исследования механических напряжений, возникающих в конструкциях, необходимо решать при проектировании мостов, арок, опор электропередачи и т. д. Рост быстроходности и удельной мощности тепловых двигателей вызывает необходимость более тщательного, чем ранее, исследования проблем механической прочности и тепловых режимов работы их деталей. Аналогичные проблемы возникают в автомобиле и турбиностроении. Проектирование дамб, плотин, дренажных и оросительных каналов невозможно без тщательного анализа течения грунтовых вод. Последняя задача является частным случаем сформулированной выше задачи о течении жидкостей и газов. В градостроительстве при проектировании системы водоснабжения городов необходим анализ течения грунтовых вод.
Анализ течения жидкого или газообразного теплоносителя на основе уравнений Навье — Стокса проводится при проектировании ядерных реакторов. Кроме того, особо важная роль при проектировании ядерных установок отводится расчету тепловыделяющей системы, математической моделью (ММ) которой является нестационарное уравнение теплопроводности. В этом случае в уравнении (1.6) дополнительно появляется член, описывающий изменение искомого температурного поля во времени. При анализе тепловых процессов в тепловыделяющих элементах (ТВЭЛах), например в высокотемпературных газоохлаждаемых реакторах, уравнение теплопроводности удобнее записывать в сферических координатах в виде
(2.6)
где Ср — удельная теплоемкость материала ТВЭЛов; R — удельная ядерная мощность.
8. Типы связей между подсистемами различной физической природы. Ранее были рассмотрены эквивалентные схемы однородных физических подсистем. Но реальный объект представляет собой совокупность разнородных физических подсистем. Согласно основным этапам получения ММС на макроуровне, после составления эквивалентных схем однородных подсистем следует установить связи между ними, т. е. определить их воздействие друг на друга. Можно выделить три типа связей: 1) трансформаторная При трансформаторном типе связи в одной подсистеме включается зависимый источник разности потенциалов. Этот источник зависит от разности потенциалов на зависимом источнике потока, установленном в другой подсистеме, который, в свою очередь, зависит от потока через первый источник. На рис. 10, а показана эквивалентная схема трансформаторной связи, где р и q - подсистемы различной природы. Если же источники в подсистемах поменять местами, то схема будет иметь такой вид, как изображено на рис. 10, б. Подобный вид связи уже встречался в ранее рассмотренных примерах эквивалентных схем однородных физических подсистем (рычаг, зубчатое зацепление колес). Для разнородных физических подсистем рассмотренная выше связь характерна при электромеханическом взаимодействии. 2) гираторная При гираторном типе связи в обеих физических подсистемах включаются зависимые источники одного вида: либо типа разности потенциалов, либо типа тока. Источники разности потенциалов зависят от потока через источник в другой подсистеме (рис. 12, а), источники потока зависят от разности потенциалов на источнике в другой подсистеме (рис. 12, б). Такой вид связи характерен при взаимодействии механической и гидравлической или пневматической подсистем. 3) через зависимые параметры элементов. ^ Связь через зависимые параметры элементов характерна при взаимодействии с тепловой подсистемой. В зависимости от значений температуры меняются параметры элементов взаимодействующей подсистемы, например вязкость жидкости, значение сопротивления резистора. Обратное влияние на тепловую подсистему учитывается введением в эквивалентную схему тепловой подсистемы зависимых источников теплового потока. При отсутствии массопереноса значение источника определяется мощностью, выделяемой элементом; при наличии массопереноса - переносимой массой. В случае моделирования влияния температуры на работу технических объектов без учета обратного воздействия предполагают, как будет изменяться температура во времени, и переносят эту зависимость на параметры элементов объекта, т. е. делают параметры зависимыми от времени. Аналогично моделируют старение элементов объекта при расчетах надежности.