- •Содержание
- •Моделирование физико-механических компонентов наносистем
- •1. Введение в теорию моделирования физико-механических компонентов микро- и наносистем
- •1.1. Моделирование процессов в конструкциях микро- и наносистем
- •1.1.1. Классификация уравнений математической физики
- •1.1.2. Анализ численных методов решения
- •Метод граничных элементов
- •Классификация вариантов метода граничных элементов
- •Сравнение методов конечных и граничных элементов
- •2. Разработка дискретных моделей ито с использованием метода конечных разностей
- •2.1. Основные положения метода конечных разностей
- •2.2. Процедура построения разностной схемы
- •2.3. Оценка погрешности дискретной модели непрерывного процесса
- •2.4. Постановка задач расчета теплового процесса на дискретной модели
- •2.4.1 Уравнение теплопередачи тепла через элемент дискретной модели
- •2.4.2. Уравнение теплопроводности Фурье для дискретной модели блока
- •2.4.3. Моделирование на эвм тепловых процессов контактных соединений
- •3. Разработка дискретных моделей ито с использованием метода конечных элементов
- •3.1. Основные положения метода конечных элементов
- •3.2. Интерполяционные полиномы для дискретизированной области
- •3.3. Решение краевых задач методом конечных элементов.
- •Литература
- •105005, Москва, 2-я Бауманская, 5
Сравнение методов конечных и граничных элементов
Размерность задачи
Метод граничных элементов по сравнению с методом конечных элементов уменьшает размерность задачи на единицу. Это означает, что для двумерных задач нужно проводить интегрирование по контуру границы, а для трехмерных только по поверхности.
Нелинейная среда, анизотропия.
Если метод граничных элементов необходимо применить для решения нелинейный задач, то можно использовать подход, который применяется в методе конечных элементов. Вся область разделяется на множество подобластей и методом итераций находится решение. Однако при слишком большом количестве подобластей метод граничных элементов вырождается в метод конечных элементов, и становится менее эффективным с вычислительной точки зрения. Более сложный подход использует интегральное преобразование нелинейного ядра. Это означает, что система преобразуется в эквивалентную линейную систему путем интегрального преобразования. При этом для каждой области, обладающей нелинейными свойствами, записывается своё интегральное уравнение, и после, при формировании системы уравнений, проводиться согласование граничных условий. Аналогично, но уже с помощью линейного преобразования, решается проблема анизотропии некоторых областей.
Непрерывное решение внутри области
Метод граничных элементов моделирует только границу. После того как найдено решение на границе, решение внутри области будет непрерывным и дифференцируемым, когда как для метода конечных элементов непрерывность решения внутри области так же, как и производных решения, ограничивается порядком аппроксимирующей модели.
Распределение ошибки
Интегральные уравнения позволяют найти точное решение задачи. При применении метода граничных элементов ошибки возникают не только при дискретизации границы, но и при численном интегрировании и использовании одной из схем аппроксимации. Если использовать сложную процедуру численного интегрирования, аппроксимацию высших порядков и наиболее подробную дискретизацию границы, то можно значительно уменьшить ошибку. В любом случае численное интегрирование представляет собой более устойчивый и точный процесс, чем численное дифференцирование, что даёт дополнительные преимущества методу граничных элементов по сравнению с методом конечных элементов.
Ранее исследователи применяли только метод конечных элементов. Затем многие стали использовать метод граничных элементов. Практика показала, что значительный выигрыш даёт применение метода граничных элементов в трехмерном случае по сравнению с методом конечных элементов. В последнее время широко применяется гибридный вариант методов конечных и граничных элементов.
Рис.1.1.4. Сравнение методов конечных и граничных элементов.
Рис.1.1.5. Пример алгоритма численного решения прямой акустической задачи.
Гибридный вариант методов конечных и граничных элементов
Из сравнительного анализа методов конечных и граничных элементов очевидны недостатки каждого метода, однако, оптимальное сочетание двух методов позволяет создать мощный и эффективный инструмент для решения целого ряда практических задач.
Гибридный метод аналогичен методу граничных элементов для многосвязанных областей. Только в этом случае для одной или нескольких областей записывается формулировка метода конечных элементов. Формулировка метода конечных и граничных элементов записывается в матричном виде, используя вначале локальную нумерацию узлов и переходя затем к общей нумерации. Подстановка начальных и граничных условий, а также условий совместимости для узлов, принадлежащих одновременно нескольким областям. Условия совместимости определяются из физических принципов решаемой задачи, например равенство потенциала и потока на границе раздела.
Рис.1.1.6. Комплексный метод «среда-структура» (МКЭ + МГЭ).
Рис.1.1.7. Алгоритм решения задачи
В результате получается сильно разряженная матрица, имеющая явно выраженную блочно-линейную структуру и зонами совместимости. Оптимальная общая нумерация позволяет использовать эффективные алгоритмы решения систем линейных алгебраических уравнений. Области, обладающие линейными изотропными свойствами, а также тонкостенные структуры лучше моделировать методом граничных элементов. Нелинейность физических свойств моделируется методом конечных элементов. При этом выбирается некоторое начальное приближение для свойств области, и при определении величины физического поля эти характеристики корректируются.
В настоящее время разработан ряд программных комплексов на основе сеточных методов, позволяющих успешно их использовать инженерам-конструкторам электронных средств для анализа динамического поведения сложных конструкций.
Программные реализации метода конечных элементов. Системы конечно элементного анализа (FEA) в зависимости от возможностей подразделяются на системы среднего уровня (NASTRAN, PATRAN, COSMOS и др.) и так называемые "тяжелые " САПР (MARC, ANSYS). Однако вне зависимости от назначения построение конечно элементных систем выполняется по традиционной схеме, показанной на рис.3.2.2.2.3.
Препроцессор (один или несколько) выполняет построение геометрии модели, задание свойств материалов и разбиение модели на конечные элементы, т.е. производит подготовку задачи к решению. При этом препроцессоры современных систем позволяют импортировать геометрию модели из конструкторских САПР (AutoCAD, SolidWorks и т.д.). Подготовленная конечно элементная модель обрабатывается вычислителем (решателем) в соответствии с необходимым видом анализа. На заключительном этапе производится визуализация расчетных данных с помощью постпроцессора.
Конечно элементные вычислители (NASTRAN, MARC, COSMOS, ANSYS) позволяют проводить линейный и нелинейный динамический анализ при гармонических и случайных воздействиях на объект исследований, представленный конечно элементной моделью. В расчетах учитываются основные характеристики материалов, включая различные типы демпфирования:
- конструкционное демпфирование;
- внутреннее демпфирование в материалах;
-вязкоупругое демпфирование дискретных элементов;
в различных представлениях:
- эквивалентное вязкоупругое;
- часть от критического;
качественный показатель.
Рис.1.1.8. Функционально-компонентный состав САПР
Дадим краткую характеристику основных пакетов:
MSC/MARC представляет собой универсальную конечноэлементную программу для проведения углубленного анализа высоконелинейного поведения конструкций и решения задач теплопередачи. Дополняя нелинейные функции продуктов MSC/NASTRAN и MSC/DYTRAN, MSC/MARC проводит комплексный анализ задач, по условиям которых конструкции подвергаются большим перемещениям и поворотам, имеют высоконелинейные свойства материалов или свойства, зависящие от истории нагружения, а также анализ контактного взаимодействия конструкций. Применение современных формулировок и вычислительных методов обеспечивает надежность проектов и сокращает объем физического макетирования. MSC/MARC обладает широкими возможностями в области решения сложных нелинейных задач. MSC/MARC может работать на различных компьютерных платформах от рабочих станций до суперкомпьютеров. С помощью MARC, в частности, можно проводить динамический анализ, анализ разрушений, акустический анализ.
MSC/NASTRAN обеспечивает полный набор расчетов, включая расчеты собственных частот и форм колебаний, установившихся и неустановившихся процессов, акустики, частотных характеристик при воздействии случайных нагрузок, спектрального анализа.
Предусмотрена возможность моделирования практически всех видов материалов, в том числе композитных и гиперупругих. Расширенные функции включают технологию суперэлементов (подконструкций), модальный синтез и макроязык DMAP для создания пользовательских приложений. Наряду с расчетом конструкций MSC/NASTRAN может использоваться и для оптимизации проектов. Оптимизацию можно проводить для решения задач статики, устойчивости, установившихся и неустановившихся переходных процессов, собственных частот и форм колебаний, акустики. Вес, напряжения, перемещения, собственные частоты и многие другие характеристики могут рассматриваться либо в качестве целевых функций проекта (в этом случае их можно минимизировать или максимизировать), либо в качестве ограничений. Алгоритмы анализа чувствительности позволяют исследовать влияние различных параметров на поведение целевой функции и управлять процессом поиска оптимального решения.
(NASTRAN выпускается и другими производителями программного обеспечения; в частности, система UAI/ NASTRAN разработана американской компанией Universal Analytics, Inc (UAI)).
COSMOS/M DESIGNER II предназначена для расчетов прочности, теплоотдачи, устойчивости на основе метода конечных элементов. COSMOS/M DESIGNER II полностью интегрирована с MicroStation и MicroStation Modeler. Исходная модель строится непосредственно в MicroStation Modeler или импортируется из файлов форматов DXF, DWG, SAT, STEP. Расчетная модель (свойства материалов, нагрузки, расчетная сетка и т.д. ) создается полуавтоматически прямо в среде MicroStation Modeler и ассоциирована с геометрической моделью. Результаты анализа представляются в виде графиков, диаграмм; с возможностью наложения карты распределения параметра прямо на модель, создание мультипликации деформации детали и т.д. Особенность этой системы заключается в наличии блока быстрого расчета (Fast Finite Element, EFE), позволяющего резко сократить в 100 раз время расчета по сравнению с базовыми продуктами конечноэлементного анализа и совместить этапы проектирования и оптимизации модели без значительного увеличения рабочего времени. Возможно, также использовать расчетные блоки систем ANSYS и NASTRAN.
COSMOS/Works представляет собой систему инженерных расчётов, созданную специально для совместного использования с системой твёрдотельного параметрического моделирования SolidWorks. Система COSMOS/Works разработана американской компанией Structural Research & Analysis Corporation и поставляется в трёх вариантах конфигурации: базовой (Basic), промежуточной (Intermediate) и расширенной (Advanced); выполняющих следующие виды расчётов:
• статический расчёт напряжений и деформаций;
• расчёт продольного изгиба для определения критической нагрузки;
• частотный анализ для определения собственных частот и форм колебаний;
• линейный динамический анализ с использованием результатов частотного анализа;
• статический и динамический нелинейный структурный анализ.
COSMOS/Works позволяет передавать результаты вычислений в системы ANSYS и MSC/NASTRAN. Кроме того, пользователь может сохранять результаты расчёта в формате VRML, а также создавать AVI -файлы на основе анимационного представления нагружения детали.
ANSYS включает в себя следующие программы анализа поведения механических конструкций:
• ANSYS/Multiphysics — самый мощный, многоцелевой продукт компании, представляет собой программное средство анализа для широкого круга инженерных дисциплин, которое позволяет проводить расчетные исследования не только в таких отдельных областях знания, как прочность, распространение тепла, механика жидкостей и газов или электромагнетизм, но и решать связанные задачи. Эта программа обеспечивает оптимизацию проектных разработок на уровне, позволяющем моделировать инженерные проблемы в наиболее полной постановке.
• ANSYS/Mechanical предназначен для выполнения проектных разработок, анализа и оптимизации: решения сложных задач прочности конструкций, теплопередачи и акустики. Эта программа, обеспечивающая проверку правильности проектных работ, является мощным инструментом для определения перемещений, напряжений, усилий, температур и давлений, а также других важных параметров.
• ANSYS/Structural выполняет сложный прочностной анализ конструкций с учетом разнообразных нелинейностей, среди которых геометрическая и физическая нелинейности, нелинейное поведение конечных элементов и потеря устойчивости. Используется для точного моделирования поведения больших и сложных расчетных моделей.
• ANSYS/ED представляет собой программу, обладающую возможностью расчетного кода ANSYS/Multiphysics, но имеющую ограничения по размерам расчетной модели. Этот доступный пакет является идеальным для обучения и подготовки персонала.
FEMAP. В настоящий момент широкое распространение приобретает система FEMAP, рекомендованная производителями международной космической станции "Альфа" в качестве основного пре-постпроцессора для подготовки и визуализации расчетных данных в системах конечноэлементного анализа.
IDEAS - комплексная система автоматизированного проектирования моделирования ИТО. Подсистема моделирования базируется на сеточных методах (МКЭ, МГЭ, МБЭ и т.п.) и позволят проводить расчеты тепловых, акустических, вибрационных и электромагнитных полей, их распространения и влияния на конструкцию. Основная платформа Sun.
ИСКРА - Интегрированная система конечно-разностного анализа, система использующая МКР для решения задач оптимизации конструкций ИТО под воздействие волновых полей.
ИСПА - Интегрированная система прочностного анализа - базируется на МКЭ и используется для анализа вибрационных и тепловых полей в ИТО. Имеет развитый интерфейс с AutoCAD. Базовая версия реализована под Windows NT.