WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 4 | 5 ||

«Руководство по основным методам проведения анализа в программе ANSYS. ...»

-- [ Страница 6 ] --

12.3. Стратегии решения сложных задач Хотя методы SUPG и COLG обеспечивают более точное решение, чем метод MSU, они могут привести к паразитным колебаниям в решении. Это может привести к нереальному решению или затруднить сходимость. Вы можете столкнуться с отрицательными температурами при решении уравнения энергии. В этом случае следует улучшить сетку в проблематичной области. Вы можете также столкнуться с отрицательными значениями турбулентности. FLOTRAN попытается приблизить эти значения к реалистичным величинам, но это может привести к ситуации, при которой произойдет останов решения. В этом случае попробуйте использовать метод MSU на величины турбулентности, оставив SUPG или COLG на параметрах импульса. Решение уравнения многокомпонентного переноса должно быть ограничено диапазоном от 0 до 1.0, и колебания в SUPG или COLG решении могут затруднить решение задачи даже при включенной опции ограничения. Если сходимость не была достигнута, рекомендуется улучшить (детализировать) конечно – элементную сетку.

Руководство по решению сопряженных задач в ANSYS Глава 1. Сопряженный анализ 1.1. Определение сопряженного анализа.

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

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

1.2. Типы сопряженного анализа.

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

1.2.1. Последовательный метод.

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

последовательно-сопряженная физика и последовательно-слабая связь.

1.2.1.1. Последовательный сопряженный анализ физики.

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

1.2.1.2. Последовательный слабо-сопряженный анализ.

Последовательно слабо-сопряженный анализ в ANSYS позволяет реализовать взаимодействие жидкость - твердое тело между элементами ANSYS FLOTRAN и структурными, тепловыми элементами или элементами сопряженного поля.

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

1.2.2. Прямой метод.

Прямой метод обычно вовлекает один анализ, который использует сопряженный тип элемента, содержащий все необходимые степени свободы. Сопряжение управляется путем вычисления матриц элемента или векторов нагрузки элемента, которые содержат все необходимые термины. Примером этого анализа является пьезоэлектрический анализ, использующий элементы SOLID5, PLANE13, или SOLID98. Другим примером является MEMS анализ с элементом TRANS126.

1.2.3. Выбор между последовательным и прямым типом анализа.

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

Для сопряженных ситуаций, которые не выражают высокую степень нелинейности взаимодействия, последовательный метод более эффективен и гибок поскольку вы можете выполнить два анализа независимо друг от друга. Сопряжение может быть рекурсивным, при выполнении итераций между различными дисциплинами до тех пор, пока не будет достигнут желаемый уровень сходимости. Например, в анализе тепло – напряжение вы можете выполнить нелинейный нестационарный тепловой анализ после линейного статического анализа напряжений. Затем можете использовать узловые температуры любого шага нагружения теплового анализа в качестве нагрузок для анализа напряжений. В последовательном слабо-сопряженном анализе вы можете выполнить нелинейный нестационарный анализ взаимодействия жидкости с твердым телом используя жидкостные элементы FLOTRAN и структурные, тепловые или сопряженные элементы ANSYS.

Указанные типы анализа описаны в следующих главах:

• Последовательный сопряженный анализ.

• Последовательный слабо-сопряженный анализ взаимодействия жидкость – твердое тело.

• Прямой сопряженный анализ • Симуляция сопряженной циркуляционной физики.

• Акустика 1.3. Система единиц измерения.

В ANSYS вы должны быть уверены в том, что используете совместимую систему единиц для всех вводимых данных. Можете использовать любую непротиворечивую систему единиц. В случае анализа электромагнитного поля, просмотрите описание команды EMUNIT в ANSYS Commands Reference для получения дополнительной информации относительно соответствующих установок для проницаемости пространства и диэлектрической проницаемости.

Для микро – электромеханических систем (MEMS), лучше установить соответствующие единицы, поскольку размер компонент может составлять всего несколько микрон. Для удобства в таблицах с 1-1 по 1-8 приводятся переводные коэффициенты из стандартных единиц MKS в единицы µMKSV и µMSVfA.

Таблица 1.1. Механические переводные коэффициенты с MKS в uMKSV Механический Единица Размерность Помножить Для Размерность параметр MKS на это число получения единиц µMKSv Длина m m 106 µm µm Сила N (kg)(m)/(s)2 106 µN (kg)(µm)/(s) Время s s 1 s s Масса kg kg 1 kg kg Давление Pa (kg)/(m)(s)2 10-6 MPa (kg)/(µm)(s) Скорость m/s m/s 106 µm/s µm/s Ускорение m/(s)2 m/(s)2 106 µm/(s)2 µm/(s) Плотность kg/(m)3 kg/(m)3 10-18 kg/(µm)3 kg/(µm) Напряжение Pa kg/(m)(s)2 10-6 MPa kg/(µm)(s) Модуль Юнга Pa kg/(m)(s)2 10-6 MPa kg/(µm)(s) Мощность W (kg)(m)2/(s)3 1012 pW (kg)(µm)2/(s) Таблица 1.2. Тепловые коэффициенты перевода с MKS в uMKSV Тепловой Единица Размерность Помножить Для Размерность параметр MKS на число получения единиц µMKSv Теплопроводность W/m°K (kg)(m)/(°K)(s)3 106 pW/(µm)°K (kg)(µm)/(°K)(s) Плотность W/(m)2 kg/(s)3 1 pW/(µm)2 kg/(s) теплового потока Теплопроводность J/(kg)°K (m)2/(°K)(s)2 1012 pJ/(kg)°K (µm)2/(°K)(s) Тепловой поток W (kg)(m)2/(s)3 1012 pW (kg)(µm)2/(s) Объемное W/m3 (kg)/(m)(s)3 10-6 pW/(µm)3 kg/(µm)(s) энерговыделение Коэффициент W/(m)2°K kg/(s)2 1 pW/(µm)2°K kg/(s) конвекции Динамическая kg/(m)(s) kg/(m)(s) 10-6 kg/(µm)(s) kg/(µm)(s) вязкость Кинематическая (m)2/s (m)2/s 1012 (µm)2/s (µm)2/s вязкость Таблица 1.3. Электрические коэффициенты перевода для MKS в uMKSV Электрический Единица Размерность Помножить Для Размерность параметр MKS на число получения единиц µMKSv Ток A A 1012 pA pA Напряжение V (kg)(m)2/(A)(s)3 1 V (kg)(µm)2/(pA)(s) Электрический Единица Размерность Помножить Для Размерность параметр MKS на число получения единиц µMKSv Разряд C (A)(s) 1012 pC (pA)(s) Теплопроводность S/m (A)2(s)3/(kg)(m)3 106 pS/µm (pA)2(s)3/(kg)(µm) Сопротивление m (Kg)(m3/(A)2(s)3 10-6 T µm (kg)(µm)3/(pA)2(s) Диэлектрическая F/m (A)2(s)4/(kg)(m)3 106 pF/µm (pA)2(s)2/(kg)(µm) постоянная [1] Энергия J (kg)(m)2/(s)2 1012 pJ (kg)(µm)2/(s) Емкость F (A)2(s)4/(kg)(m)2 1012 pF (pA)2(s)4/(kg)(µm) Электрическое V/m (kg)(m)/(s)3(A) 10-6 V/µm (kg)(µm)/(s)3(pA) поле Плотность потока C/(m)2 (A)(s)/(m)2 1 pC/(µm)2 (pA)(s)/(µm) электрической индукции 1. Проницаемость пространства равна 8.854 x 10-6 pF/µm.

Таблица 1.4. Магнитные коэффициенты перевода с MKS в uMKSV Магнитный Единицы Размерность Помножить Для Размерность параметр MKS на число получения единиц µMKSv Поток Weber (kg)(m)2/(A)(s)2 1 Weber (kg)(µm)2/(pA)(s) Плотность Tesla kg/(A)(s)2 10-12 Tesla kg/(pA)(s) потока Напряженность A/m A/m 106 pA/µm pA/µm поля Ток A A 1012 pA pA Плотность тока A/(m)2 A/(m)2 1 pA/(µm)2 pA/(µm) Проницаемость H/m (kg)(m)/(A)2(s)2 10-18 TH/µm (kg)(µm)/(pA)2(s) [1] Индуктивность H (kg)(m)2/(A)2(s)2 10-12 TH (kg)(µm)2/(pA)2(s) 1. Проницаемость пространства равна 4 x 10-25 TH/µm.

Примечание В этих единицах может использоваться только постоянная проницаемость.

Таблица 1.5. Механические коэффициенты перевода с MKS в uMSVfA Механический MKS Размерность Помножить на Для Размерность параметр это число получения µMsvfa Механический MKS Размерность Помножить на Для Размерность параметр это число получения µMsvfa Длина m m 106 µm µm Сила N (kg)(m)/(s)2 109 nN (g)(µm)/(s) Время s s 1 s s Mass kg kg 103 g g Давление Pa (kg)/(m)(s)2 10-3 kPa g/(µm)(s) Скорость m/s m/s 106 µm/s µm/s Ускорение m/(s)2 m/(s)2 106 m/(s)2 µm/(s) Плотность kg/(m)3 kg/(m)3 10-15 g/(µm)3 g/(µm) Давление Pa kg/(m)(s)2 10-3 kPa g/(µm)(s) Модуль Юнга Pa kg/(m)(s)2 10-3 kPa g/(µm)(s) Мощность W (kg)(m)2/(s)3 1015 fW (g)(µm)2/(s) Таблица 1.6. Температурные коэффициенты перевода с MKS в uMSVfA Тепловой MKS Размерность Помножить Для Размерность параметр на это получения число µMsvfa Теплопроводность W/m°K (kg)(m)/(°K)(s)3 109 fW/(µm)°K (g)(µm)/(°K)(s) Плотность W/(m)2 kg/(s)3 103 fW/(µm)2 g/(s) теплового потока Теплоемкость J/(kg)°K (m)2/(°K)(s)2 1012 fJ/(g)°K (µm)2/(°K)(s) Тепловой поток W (kg)(m)2/(s)3 1015 fW (g)(µm)2/(s) Объемное W/m3 (kg)/(m)(s)3 10-3 fW/(µm)3 g/(µm)(s) энерговыделение Коэффициент W/(m)2°K kg/(s)2 103 fW/(µm)2°K g/(s) конвекции Динамическая kg/(m)(s) kg/(m)(s) 10-3 g/(µm)(s) g/(µm)(s) вязкость Кинематическая (m)2/s (m)2/s 1012 (µm)2/s (µm)2/s вязкость Таблица 1.7. Электрические коэффициенты перевода с MKS в uMSVfA Электрический MKS Размерность Помножить Для Размерность параметр на это получения число µMsvfa Ток A A 1015 fA fA Напряжение V (kg)(m)2/(A)(s)3 1 V (g)(µm)2/(fA)(s) Разряд C (A)(s) 1015 fC (fA)(s) Электропроводность S/m (A)2(s)3/(kg)(m)3 109 nS/µm (fA)2(s)3/(g)(µm) Электрический MKS Размерность Помножить Для Размерность параметр на это получения число µMsvfa Сопротивление m (kg)(m3/(A)2(s)3 10-9 - (g)(µm)3/(fA)2(s) Диэлектрическая F/m (A)2(s)4/(kg)(m)3 109 fF/µm (fA)2(s)2/(g)(µm) постоянная [1] Энергия J (kg)(m)2/(s)2 1015 fJ (g)(µm)2/(s) Емкость F (A)2(s)4/(kg)(m)2 1015 fF (fA)2(s)4/(g)(µm) Электрическое поле V/m (kg)(m)/(s)3(A) 10-6 V/µm (g)(µm)/(s)3(fA) Плотность потока C/(m)2 (A)(s)/(m)2 103 fC/(µm)2 (fA)(s)/(µm) электрической индукции 1. Диэлектрическая постоянная пространства равна 8.854 x 10-3 fF/µm.

Таблица 1.8. Магнитные коэффициенты перевода с MKS на uMKSVfA Магнитный MKS Размерность Помножить Для Размерность параметр на число получения µMKSv Поток Weber (kg)(m)2/(A)(s)2 1 Weber (g)(µm)2/(fA)(s) Плотность потока Tesla kg/(A)(s)2 10-12 - g/(fA)(s) Напряженность A/m A/m 109 fA/µm fA/µm поля Ток A A 1015 fA fA Плотность тока A/(m)2 A/(m)2 103 fA/(µm)2 fA/(µm) Проницаемость[1] H/m (kg)(m)/(A)2(s)2 10-21 - (g)(µm)/(fA)2(s) Индуктивность H (kg)(m)2/(A)2(s)2 10-15 - (g)(µm)2/(fA)2(s) 1. Проницаемость пространства равна 4 x 10-28 (g)(µm)/(fA)2(s)2.

Примечание В этих единицах может быть использована только постоянная проницаемость.

1.4. О маршрутах GUI и синтаксисе команд.

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

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

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

Глава 2. Анализ последовательно-сопряженной физики.

2.1. Определение анализа последовательно-сопряженной физики.

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

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

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

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

• Основные физические нагрузки не являются функцией другого физического анализа.

Такие нагрузки также называются номинальными граничными условиями.

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

Ниже приведены типичные задачи, которые можно решить при помощи программы ANSYS:

• Тепловые напряжения.

• Индукционный нагрев.

• Стационарное взаимодействие жидкости со структурой.

• Магнитно-структурное взаимодействие.

• Электростатическое – структурное взаимодействие.

• Проводимость – электростатика.

Программа ANSYS может выполнить мультифизический анализ с одной базой данных. Вся модель, при выполнении такого анализа, будет содержать один набор узлов и элементов. Что будут представлять эти элементы зависит от реализуемой физики.

2.2. Определение физической среды.

Программа ANSYS выполняет анализ последовательно-сопряженной физики, используя концепцию физической среды. Термин физической среды применим как к создаваемому вами файлу, который содержит все рабочие параметры и характеристики определенного физического анализа, так и к содержимому файла. Файл физической среды является ASCII файлом, который вы создаете, используя следующие методы:

Команда:

PHYSICS,WRITE,Title,Filename,Ext,-- GUI:

Main Menu> Preprocessor> Physics> Environment Main Menu> Solution> Physics> Environment Вы можете определить вплоть до девяти физических сред. Для каждой среды определяется уникальный заголовок командой PHYSICS. ANSYS присваивает каждой физической среде уникальный номер как часть расширения файла, определяющего эту среду. Рекомендуется использование заголовка, описывающего физику анализа. Заголовок также должен отличаться от заголовка анализа, определяемого командой /TITLE (Utility Menu> File> Change Title).

Команда PHYSICS,WRITE создает файл физической среды (например, Jobname.PH1) путем выборки следующей информации из базы данных ANSYS:

• Типы элементов и KEYOPT установки.

• Вещественные константы.

• Свойства материала.

• Система координат элемента.

• Опции решения анализа.

• Опции шага нагружения.

• Уравнения ограничения.

• Наборы сопряженных узлов.

• Заданные граничные условия и нагрузки.

• Установки предпочтения GUI.

• Заголовок анализа (/TITLE card) Можете также использовать команду PHYSICS,READ (Main Menu> Preprocessor> Physics> Environment> Read) для чтения файла физической среды, используя имя файла или заголовок, использованный при написании файла (этот заголовок включается в качестве комментария в начало файла физической среды). Перед чтением физического файла, программа ANSYS удаляет все граничные условия, нагрузки, узловые сопряжения, свойства материалов, опции анализа и уравнения ограничения хранящиеся в базе данных.

2.3. Основные процедуры анализа.

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

В косвенном методе используются разные базы данных и файлы результатов.

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

Рисунок 2.1. Поток данных в последовательном сопряженном анализе.

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

Рисунок 2.2. Поток данных в Последовательном сопряженном физическом анализе (используются физические среды).

Создавайте базу данных ANSYS помня о требованиях, предъявляемым к каждой физической среде. Перед созданием любой физической среды, назначьте номер типа элемента, номер материала, номер набора вещественных констант и номер набора элементных координат для каждого региона, площади или объема (посмотрите описание команд AATT и VATT).

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

Команда:

EMODIF GUI:

Main Menu> Preprocessor> Modeling> Move/Modify> Elements> Modify Attrib Косвенный метод идеален для одностороннего последовательного сопряжения, такого как типичный анализ тепло – напряжения. Приближение физической среды позволяет быстро переключаться между физическими средами, что идеально подходит для полностью сопряженных сценариев, требующих многократных прохождений между физическими решениями. Стационарное взаимодействие типа жидкость-структура или индукционный нагрев являются типичными примерами случаев, требующих использования физических сред.

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

• Выполните команду SAVE после создания физических сред и команду RESUME после каждого физического решения.

• Не записывайте результатов в базу данных (записывайте только в файл результатов).

Затем необходимо выполнить команду SET каждый раз перед считыванием данных из фала результатов в базу данных для постобработки. Для задания этой опции выполните команду /CONFIG,NOELAB,1 или вставьте линию "NO_ELDBW = 1" в файл config70.ans.

2.4. Передача нагрузок между дисциплинами.

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

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

Результаты анализа... Становятся нагрузками для следующего типа анализа...

Температуры теплового анализа или Массовая сила для структурного анализа FLOTRANа [TEMP, TBOT, TE2,... TTOP] или узловые нагрузки (температуры) для теплового анализа.

Силы статического, гармонического или Силовые нагрузки для структурного или нестационарного магнитного анализа флотрановского анализа.

Результаты анализа... Становятся нагрузками для следующего типа анализа...

[FORC] Силы, полученные в результате решения Силовые нагрузки для структурного электростатического анализа [FORC] анализа.

Энерговыделение полученное в результате Энерговыделение для теплового или решения магнитного анализа. [HGEN] флотрановского анализа.

Плотность электрического тока из анализа Плотность тока для анализа магнитного проводимости [JS] поля Давления из флотрановского анализа Поверхностные нагрузки (давления) для [PRES] структурного анализа (твердые или оболочечные элементы).

Реакционные нагрузки из любого анализа Силовые нагрузки для любого анализа [REAC] Плотности теплового потока из Поверхностные нагрузки (плотность флотрановского анализа [HFLU] потока) на элементы в тепловом анализе.

Плотности теплового потока из анализа Поверхностные нагрузки (плотность электромагнитного поля высокой частоты теплового потока) на элементы в тепловом [EHFLU] анализе.

Вычисленный во FLOTRANе коэффициент Поверхностные нагрузки (коэффициенты теплоотдачи и соответствующая теплоотдачи и температура среды) на температура среды [HFLM] элементы в тепловом анализе.

2.4.1. Совместимые типы элементов.

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

Базовая геометрия Базовая геометрия элемента имеет принятую по умолчанию конфигурацию, описываемую в “Руководстве по элементам ANSYS”. Для твердых элементов базовая геометрия включает четырехугольную, треугольную, шестигранную (блочную) и четырехгранную формы.

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

Порядок элемента Различают два порядка элементов ANSYS (включая p-элементы): элементы низкого (первого) порядка и элементы высокого (второго) порядка. Элементы высокого порядка имеют промежуточные узлы;

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

В мультифизических средах типы элемента должны поддерживать совместимую базовую геометрию. Если тип элемента допускает вырожденную геометрию, соответствующий тип элемента в других средах также должен поддерживать вырожденную геометрию. Например, SOLID92 (10-узловой четырехгранный структурный твердый элемент) совместим с SOLID (10-узловой четырехгранный тепловой твердый элемент). Однако, SOLID92 не совместим с вырожденной четырехгранной формой SOLID90 (20-узловой твердый тепловой элемент).

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

Ниже приведены нагрузки, которые можно считать с элементов первого или второго порядка и задать на элементы первого или второго порядка для решения следующей физической среды:

• Температуры [TEMP, TBOT, TE2,... TTOP] • Энерговыделение [HGEN] • Плотность тока [JS] • Давление [PRES] • Плотность теплового потока [HFLU] • Коэффициент теплоотдачи и температура граничной среды [HFLM] Нагрузки, требующие совместимости порядка элементов:

• Силовые нагрузки* [FORC] • Реакционные нагрузки [REAC] * Следующие электромагнитные элементы поддерживают структурные элементы первого и второго порядка: PLANE53, PLANE121, SOLID122, и SOLID123.

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

Структурный Тепловой Магнитный Электростатический Жидкостный Электрическая проводимость SOLID45 SOLID70 SOLID97, SOLID122[2] FLUID142 SOLID5, SOLID117[1] SOLID SOLID92 SOLID87 SOLID98, SOLID123[2] - SOLID HF119[1] SOLID95 SOLID90 SOLID117, SOLID122 - SOLID5, HF120 SOLID PLANE42 PLANE55 PLANE13, PLANE121[2] FLUID141 PLANE PLANE53[2] PLANE2 PLANE35 - - - PLANE82 PLANE77 PLANE53 PLANE121 - PLANE Структурный Тепловой Магнитный Электростатический Жидкостный Электрическая проводимость SHELL63, SHELL57, - - - SHELL SHELL181 SHELL SHELL91, SHELL132 - - SHELL LINK8 LINK33 - - FLUID116 LINK Примечание Если конечно элементная сетка использует вырожденную форму элемента, соответствующий тип элемента должен поддерживать ту же вырожденную форму. Например, если при построении к/э сетки используются пирамидальные элементы FLUID142, то элементы SOLID70 не совместимы, поскольку элементы SOLID70 не могут иметь пирамидальную форму. Для совместимости, элементы со степенью свободы VOLT должны иметь такую же реакционную силу (смотрите Совместимость элементов в Руководстве по электромагнитному анализу в ANSYS).

1. Поддерживает только элементы первого порядка.

2. Требуется KEYOPT опция для поддержания элементов первого порядка.

2.4.2. Типы файлов результатов, которые вы можете использовать.

Обычно в косвенном совместном анализе вы работаете с несколькими разными типами файлов результатов. Все файлы результатов для вашего анализа будут иметь одно и то же имя (рабочее имя, задаваемое командой /FILNAME (Utility Menu> File> Change Jobname)).

Однако, вы можете заметить различие этих файлов по их расширению.

Jobname.RFL Файл результатов FLOTRAN Jobname.RMG Файл результатов электромагнитного анализа.

Jobname.RTH Файл результатов теплового анализа.

Jobname.RST Все другие типы файлов результатов (структурного и анализа мультифизики).

2.4.3. Нестационарный жидкостно-структурный анализ.

В нестационарном жидкостно-структурном анализе вы можете выполнить структурный анализ в промежуточные моменты времени, соответствующие изменениям граничных условий жидкости. Например, необходимо выполнить структурный анализ в 2 секунды и входная скорость изменяется от 1.0 in/sec в 0.0 секунд до 5.0 in/sec в 4.0 сек. Вначале вы, как обычно выполняете структурный анализ в момент времени 2 с. Когда выполняется команда PHYSICS,READ,FLUID (Main Menu> Solution> Physics> Environment> Read) для возобновления жидкостного анализа, вы повторно задаете линейно изменяющуюся нагрузку.

Вы задаете входную граничную скорость 3.0 in/sec в 2.0 секунды и затем указываете на то, что это “старое” условие следующим образом:

Команда:

FLOCHECK, GUI:

Main Menu> Preprocessor> FLOTRAN Set Up> Flocheck Это означает, что входное граничное условие 3.0 in/sec в момент времени 2 секунды является начальной точкой для линейного изменения. Затем вы вводите финальную точку для линейного изменения, 5.0 in/sec в момент времени 4 секунды, и задаете линейно изменяющееся граничное условие следующим образом:

Команда:

FLDATA4,TIME,BC, GUI:

Main Menu> Preprocessor> FLOTRAN Set Up> Execution Ctrl Нестационарный анализ, как обычно, запускается командой SOLVE.

Более подробную информацию по заданию нестационарных граничных условий во FLOTRAN смотрите в Нестационарном анализе во FLOTRAN.

2.5. Выполнение Последовательного Совместного Анализа при использовании Физических Сред.

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

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

• Каждая площадь или объем твердой модели ANSYS имеет свои собственные специфичные требования к типу элемента, свойствам материала, и вещественным константам. Все объекты твердой модели должны иметь заданные номера типа элемента, набора вещественных констант, материала и элементной системы координат (из значения изменяются в соответствии с физической средой).

• Определенные группы площадей или объемов будут использованы в двух или более разных физических средах. Используемая конечно элементная сетка должна подходить для всех сред.

2. Создайте физическую среду. Выполняйте этот шаг для каждой физической дисциплины, которая является частью последовательного совместного анализа.

• При необходимости обращайтесь к различным разделам справочного описания ANSYS для того, чтобы определить, что вам следует задавать для определенного физического анализа.

• Задайте необходимые типы элементов, которые будут использоваться в физической симуляции (например, ET,1,141 или ET,2,142 для симуляции во FLOTRANе, ET,1,13 или ET,2,117 для магнитного решения и так далее).

Установите нулевой тип элемента (тип = 0, например ET,3,0) для использования в областях не связанных с данной физикой. Во время решения элементы нулевого типа игнорируются.

• Назначьте свойства материалов, набор вещественных констант и элементные системы координат в соответствии с установленными номерами атрибутов заданными ранее в модели.

• Назначьте номера атрибутов для типа элементов, материалов, вещественных констант и элементной системы координат площадей или объемов твердой модели (для этого воспользуйтесь командой AATT (Main Menu> Preprocessor> Meshing> Mesh Attributes> All Areas or Picked Areas) или командой VATT (Main Menu> Preprocessor> Meshing> Mesh Attributes> All Volumes or Picked Volumes)).

• Задайте номинальные нагрузки и граничные условия.

• Установите все опции решения.

• Выберете название физической среды и выполните команду PHYSICS,WRITE с этим названием. Например, в жидкостно-магнитном анализе вы могли бы использовать следующую команду для записи жидкостной физической среды:

Команда:

PHYSICS,WRITE,Fluids GUI:

Main Menu> Preprocessor> Physics> Environment> Write • Для создания следующей физической среды очистите базу данных текущей физической среды. Это можно выполнить следующим образом:

Команда:

PHYSICS,Clear GUI:

Main Menu> Preprocessor> Physics> Environment> Clear • Используя описанную выше процедуру, подготовьте следующую физическую среду.

• Для сохранения базы данных и указателей на файлы физики выполните команду SAVE.

Предположим, что имя мультифизического анализа "Induct", тогда первые два файла физической среды будут названы соответственно Induct.PH1 и Induct.PH2. Более подробную информацию по команде PHYSICS смотрите в Руководстве по командам ANSYS.

3. Проведем последовательный совместный анализ, выполнив по очереди каждый физический анализ.

4. /SOLU ! Вход в раздел решения 5. PHYSICS,READ,Magnetics ! Содержит магнитную среду 6. SOLVE 7. FINISH 8. /SOLU 9. PHYSICS,READ,Fluids 10. LDREAD,FORCE,,,,2,,rmg ! Магнитные силы Лоренца SOLVE Расширение, используемое в команде LDREAD связано со считываемым файлом результатов. Результаты теплового анализа считываются из файла Jobname.RTH. Все другие результаты кроме результатов магнитного и жидкостного анализа находились бы в файле Jobname.RST.

2.5.1. Обновление конечно элементной сетки.

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

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

Запуск симуляции этой задачи требует подстройки сетки поля таким образом, чтобы она совпадала с деформированной структурной сеткой. В ANSYS эта подстройка называется трансформацией сетки.

Для трансформации сетки используются следующие команды: DAMORPH (трансформирует смежные с площадями элементы), DVMORPH (трансформирует смежные с объемами элементы), DEMORPH (трансформирует выбранные элементы). Используйте опцию RMSHKY для задания одного из следующих способов трансформации сетки:

• Морфинг – программа передвигает узлы и элементы сетки “поля” до совпадения с деформированной структурной сеткой. В этом случае, программа не создает каких либо новых узлов или элементов.

• Повторное наложение сетки – удаляет сетку поля и заменяет ее новой к/э сеткой совпадающей с деформированной структурной сеткой. Повторное наложение сетки не меняет структурную сетку. Данный вид трансформации соединяет новую сетку поля с существующими узлами и элементами деформированной структурной сетки.

• Морфинг или повторное наложение сетки – программа вначале пытается трансформировать сетку поля. В случае неудачи, программа переключается в режим повторного наложения сетки в выбранной области поля. Эта последовательность принята по умолчанию.

Морфинг сетки оказывает влияние только на узлы и элементы. Морфинг не влияет на составляющие твердой модели (ключевые точки, линии, площади, объемы). Морфинг оставляет связь узлов и элементов с составляющими твердой модели. Узлы и элементы, прикрепленные к ключевым точкам, линиям и площадям, находящимся внутри выбранной для морфинга области могут, на самом деле, сдвинуть эти составляющие, однако, связь останется прежней.

Будьте внимательны при задании граничных условий и нагрузок на выбранную для морфинга область. Граничные условия и нагрузки, заданные на узлы и элементы подходят только для опции морфинга. Если граничные условия и нагрузки задаются непосредственно на узлы и элементы, то команды DAMORPH, DVMORPH, и DEMORPH требуют, чтобы эти нагрузки были удалены перед повторным наложением сетки. Граничные условия и нагрузки, заданные на элементы твердой модели будут корректно переданы на новую сетку.

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

Также будьте внимательны при задании начальных условий командой IC. Перед выполнением структурного анализа команды DAMORPH, DVMORPH, и DEMORPH требуют, чтобы начальные условия были удалены со всех узлов элементов нулевого типа в неструктурных областях. Для удаления начальных условий используйте команду ICDELE.

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

Смещения, полученные в результате проведения структурного анализа должны быть в базе данных перед выполнением команды трансформации. Результаты будут находиться в базе после структурного решения или после чтения файла результатов (команда SET в POST1).

Структурные узлы модели переместятся в деформированную позицию, определенную по вычисленным смещениям. Если вы проводите последующий структурный анализ, то в данном случае следует восстанавливать первоначальное положение структурных узлов. Это можно сделать, выбрав структурные узлы и после этого выполнив команду UPCOORD с FACTOR = -1. Команда:

UPCOORD,Factor GUI:

Main Menu> Solution> Load Step Opts> Other> Updt Node Coord Морфинг сетки поддерживает все двумерные модели с сеткой, состоящей из четырехугольных, треугольных элементов высокого или низкого порядка. Для двумерных моделей, все узлы и элементы должны быть в одной и той же плоскости. Не поддерживаются искривленные поверхности. В 3-D поддерживаются модели со следующими опциями морфинга и конфигурацией формы:

• Все четырехгранные элементы - (поддерживается морфинг и повторное неложение сетки) • Все блочные элементы - (поддерживается морфинг) • Все клинообразные элементы - (поддерживается морфинг) • Комбинация элементов пирамида - четырехгранник - (морфинг) • Комбинация элементов блок - клин - (морфинг) Морфинг наиболее вероятно подойдет для сеток с единообразными элементами (элементы созданные при помощи SMRTSIZE). Морфинг не применим к сильно деформированным элементам.

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

Модель бруска и воздуха Во многих случаях только часть модели требует проведения морфинга (то есть область, находящаяся в непосредственной близости от структурной области). В этом случае вам следует выбрать площади или объемы, расположенные в непосредственной близости от структурной модели. На приведенном ниже рисунке показан пример бруска расположенного в нескольких электростатических областях. Только площадь 3 требует трансформации сетки.

Для поддержания совместимости сетки с нетрансформированным регионом, алгоритм морфинга не влияет на узлы и элементы, которые граничат с трансформируемыми площадями или объемами. В этом примере, алгоритм не повлияет на узлы поверхности раздела 2 и 3 площадей.

Модель бруска с несколькими воздушными областями.

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

DAMORPH, DVMORPH, DEMORPH GUI:

Main Menu> Preprocessor> Meshing> Modify Mesh> Refine At> Areas Main Menu> Preprocessor> Meshing> Modify Mesh> Refine At> Volumes Main Menu> Preprocessor> Meshing> Modify Mesh> Refine At> Elements 2.5.2. Повторный запуск анализа при использовании метода физической среды.

Во многих последовательных совместных задачах существует необходимость повторного запуска одного из физических решений. Например, в индукционном нагреве вам необходимо повторно запустить нестационарный тепловой анализ в процессе последовательных сопряженных циклов. Для статического нелинейного структурного анализа сопряженного поля предпочтительно возобновить структурное решение вместо того, что бы запустить его сначала. Вы можете легко реализовать процедуру возобновления анализа, используя последовательный анализ сопряженного поля. Повторный запуск требует наличия EMAT, ESAV, и DB файлов для определенных физик. Вы можете изолировать EMAT и ESAV файлы для определенной физики при помощи команды /ASSIGN. Файл базы данных будет совместим с физикой, при использовании метода физической среды. Краткое описание процедуры возобновления решения:

1. Используйте команду /ASSIGN для переадресации EMAT и ESAV файлов перед решением области физики требующей возобновления решения.

2. Выполните повторный запуск решения.

3. Используйте команду /ASSIGN для переадресации EMAT и ESAV файлов к их принятым по умолчанию значениям для использования другими областями физики.

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

2.6. Настраиваемые параметры мультифизического решателя.

Выполнение Последовательного Сопряженного Анализа с Физическими Средами описывает метод последовательного сопряженного анализа, использующего приближение физических сред. Для многих задач этот метод может быть реализован при помощи макроса. Этот раздел подробно описывает настраиваемые макросы решения созданные для определенных задач.

2.6.1. Электростатический структурный сопряженный решатель.

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

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

2.6.1.1. Требования Макрос ESSOLV использует физические среды, как описано в Выполнении Последовательного Сопряженного анализа с Физическими средами. Таким образом, макрос требует, чтобы файл электростатической физики и файл структурной физики были созданы перед решением.

Решатель применим для двумерных моделей, включающих четырехгранные или треугольные элементы и трехмерных моделей использующих опцию формы элемента, подробно описанную в Обновлении к/э сетки. Решатель будет работать наилучшим образом в случае, если отношение сторон элементов конечно элементной модели (отношение длины к ширине граней элемента) не будет значительно отличаться от единицы. Элементы сгенерированные ANSYS при использовании SMRTSIZE в основном имеют приемлемую форму для этого типа симуляции.

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

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

При использовании команды ESSOLV со структурными оболочечными элементами (например, SHELL63, SHELL93), вы должны вручную задать флаг поверхности Максвелла на все элементы воздуха, окружающие оболочки перед записью файла электростатической физики. Используйте команду SFA для задания флага поверхности Максвелла на площади, представляющие оболочечные элементы. Это будет гарантировать то, что элементы воздуха с обоих сторон оболочек получат флаг поверхности Максвелла.

Примечание Если используются оболочки или структурные элементы низкого порядка, установите KEYOPT(7) = 1 для электростатических типов элементов что гарантированно обеспечит передачу сил.

2.6.1.2. Процедура Существует следующая процедура подготовки задачи, которая впоследствии будет решена при помощи макроса ESSOLV.

• Постройте твердую модель, включающую в себя электростатическую и структурную области. Наложите сетку на обе области.

• Создайте электростатическую физическую среду, назначив соответствующие типы элементов мешированной области, определив свойства материалов, задав граничные условия твердой модели и так далее. Для структурной области установите нулевой тип элемента (ET,,0). Запишите электростатическую физику в файл физики (PHYSICS,WRITE).

• Очистите электростатическую физику (PHYSICS,CLEAR) и установите физику для структурного анализа, выбрав соответствующий тип элементов, задав свойства материалов, определив граничные условия и нагрузки твердой модели, выбрав решатель, опции и так далее. Запишите структурную физику в файл физики PHYSICS,WRITE).

• Подготовьте твердую модель для морфинга (этот шаг необязателен и может быть выполнен интерактивно при помощи диалогового окна макроса ESSOLV). В двумерном анализе вы можете сгруппировать площади электростатической области, которая впоследствии будет трансформирована, в компонент площади. В трехмерном анализе ради удобства можно сгруппировать объемы в объемный компонент.

• Подготовьте компонент для исключения из морфа (этот шаг необязателен и может быть выполнен интерактивно в GUI после запуска ESSOLV). Можете выбрать линии (в двумерной модели) или площади (в трехмерной модели) электростатической области которую вы не хотите подвергать морфингу в процессе решения. Это зачастую необходимо выполнить в том случае, если морфинг будет расположен в небольшой по сравнению со всей моделью области.

• Если необходимо задать предварительные напряжения на структуру, подготовьте файл начальных напряжений. В дополнении создайте компонент элементов, которые содержаться в файле начальных напряжений. Формат файла начальных напряжений описан в Задание начальных напряжений в Руководстве по основному анализу в ANSYS.

• Выполните командный макрос ESSOLV.

2.6.1.3. Решение.

Командный макрос будет итерировать между электростатическим и структурным файлом физик до тех пор, пока изменение в структурном отклонении и накопленная электростатическая энергия в системе не достигнут определенного заранее значения. По умолчанию решение сходиться, когда максимальное структурное отклонение и накопленная электростатическая энергия будут меняться менее, чем на 0.5% между текущей и предыдущей итерациями. Вы можете изменить или деактивировать критерий сходимости, используя параметры ELECTOL и STRUTOL команды ESSOLV.

Для того чтобы убедиться в том, что ваша модель как следует подготовлена, вы можете выполнить анализ, состоящий из одного цикла, установив параметр MXLOOP = 1 в команде ESSOLV.

При успешном решении задачи макрос ESSOLV обновляет узловые координаты структуры и электростатической области. Для возврата к исходной геометрии, возобновите файл базы данных _ESSOLV.DB командой RESUME.

Ниже приведен алгоритм макроса ESSOLV.

2.6.1.4. Постобработка.

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

однако все вычисленные параметры будут точны, и вы можете просмотреть все остальные результаты. Если вам необходимо точное представление структурных смещений, возвратите структуру в исходное состояние, используя команду UPCOORD,-1.

Для постобработки структурных результатов выполните следующее:

1. Считайте файл структурной физической среды следующим образом:

Команда:

PHYSICS,READ GUI:

Main Menu> Preprocessor> Physics> Environment> Read 2. Считайте результаты из файла структурных результатов:

Команда:

SET GUI:

Main Menu> General Postproc> Read Results> Last Set 3. Если требуется точное представление смещений, выберете структурные узлы и возвратите структурную геометрию к исходному состоянию:

Команда:

NSEL GUI:

Utility Menu> Select> Entities Команда:

UPCOORD,-1.

GUI:

Main Menu> Preprocessor> Loads> Load Step Opts> Other> Updt Node Coord Если вы желаете вернуть деформированную геометрию, выполните команду UPCOORD,1. Для просмотра электростатических результатов выполните следующую процедуру.

Для постобработки электростатических результатов выполните следующее:

1. Считайте файл электростатической среды:

Команда:

PHYSICS,READ GUI:

Main Menu> Preprocessor> Physics> Environment> Read 2. Считайте результаты из файла электростатических результатов следующим образом:

Команда:

SET GUI:

Main Menu> General Postproc> Read Results> Last Set 2.6.1.5. ESSOLV мультирешения.

Зачастую возникает необходимость вычислить структурную деформацию в диапазоне заданных напряжений. Можете использовать ESSOLV для запуска последовательных решений. В этом случае используйте команды *DO и *ENDDO для задания цикла. Вы должны определить новые нагрузки и соответственно заново записать файл физики. Для второго и последующих вызовов ESSOLV, установите значение RUSEKY больше единицы. Не используйте команду SAVE для сохранения базы данных между вызовами ESSOLV.

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

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

2.6.2. Сопряженный решатель задач типа жидкость-структура.

Командный макрос FSSOLV является решателем, с помощью которого можно решить стационарную CFD (расчетная динамика жидкости) и структурную задачи. Этот решатель не применим к решению динамических задач взаимодействия жидкости со структурой (FSI).

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

2.6.2.1. Требования.

Макрос FSSOLV использует физические среды, как описано в Выполнении Последовательного Сопряженного анализа с Физическими средами. Таким образом, макрос требует, чтобы файл электростатической физики и файл структурной физики были созданы перед решением.

Решатель применим для двумерных моделей, включающих четырехгранные или треугольные элементы и трехмерных моделей использующих опцию формы элемента, подробно описанную в Обновлении к/э сетки. Решатель будет работать наилучшим образом в случае, если отношение сторон элементов конечно элементной модели (отношение длины к ширине граней элемента) не будет значительно отличаться от единицы. Элементы сгенерированные ANSYS при использовании SMRTSIZE в основном имеют приемлемую форму для этого типа симуляции.

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

2.6.2.2. Процедура Существует следующая процедура подготовки задачи, которая впоследствии будет решена при помощи макроса FSSOLV.

• Постройте твердую модель, включающую в себя жидкую и структурную области.

Наложите сетку на обе области.

• Создайте жидкостную физическую среду, назначив соответствующие типы элементов мешированной области, определив свойства материалов, задав граничные условия твердой модели и так далее. Для структурной области установите нулевой тип элемента (ET,,0). Запишите жидкостную физику в файл физики (PHYSICS,WRITE).

• Очистите жикостную физику (PHYSICS,CLEAR) и установите физику для структурного анализа, выбрав соответствующий тип элементов, задав свойства материалов, определив граничные условия и нагрузки твердой модели, выбрав решатель, опции и так далее. Запишите структурную физику в файл физики PHYSICS,WRITE).

• Подготовьте твердую модель для морфинга (этот шаг необязателен и может быть выполнен интерактивно при помощи диалогового окна макроса FSSOLV). В двумерном анализе вы можете сгруппировать площади электростатической области, которая впоследствии будет трансформирована, в компонент площади. В трехмерном анализе ради удобства можно сгруппировать объемы в объемный компонент.

• Подготовьте компонент для исключения из морфа (этот шаг необязателен и может быть выполнен интерактивно в GUI после запуска FSSOLV). Можете выбрать линии (в двумерной модели) или площади (в трехмерной модели) электростатической области которую вы не хотите подвергать морфингу в процессе решения. Это зачастую необходимо выполнить в том случае, если морфинг будет расположен в небольшой по сравнению со всей моделью области.

• Выполните командный макрос FSSOLV.

2.6.2.3. Solution Макрос будет итерировать между жидкостным и структурным файлом физики до тех пор, пока изменение максимального структурного отклонения и общая равнодействующая сила жидкости, а также момент жидкости не достигнут предопределенного значения. По умолчанию решение сходится когда указанные параметры (максимальное структурное смещение, сила жидкости и величина момента) изменяются менее чем на 0.5% между текущим и предыдущим итеративным циклом. Вы можете изменить или деактивировать критерий сходимости, используя следующие параметры команды FSSOLV: FORCTOL, MOMETOL, и STRUTOL.

По умолчанию максимальное количество глобальных итераций жидкости (команда FLDATA,ITER,EXEC) выполняемых в ходе CFD решения определяется установленным значением при записи файла физики. После начального жидкостно-структурного решения, максимальное количество итераций, необходимых для получения нормального решения может быть меньше. Для изменения этого значения используйте опцию FLUITER.

Для того, чтобы убедиться что ваша модель готова, выполните анализ, состоящий из одной итерации установив опцию MXLOOP = 1 в команде FSSOLV.

В случае успешного решения задачи, макрос FSSOLV обновляет узловые координаты структуры и жидкой области для выравнивания со структурной деформацией. Для отката к первоначальной геометрии, возобновите файл _FSSOLV.DB посредством команды RESUME.

Ниже приведен алгоритм, реализуемый макросом FSSOLV.

Пример стационарного жидкостно-структурного анализа, использующего макрос FSSOLV расположен в конце этой главы. В задаче вычисляются смещения тонкого алюминиевого бруска, расположенного в канале с текущей водой. Решение сходится после шести итераций, использующих принятые по умолчанию значения сходимости. Окончательный результат отклонения конца бруска показывает более чем 20% разницу по сравнению с начальным структурным решением, иллюстрируя важность решения данного типа задач рекурсивным методом.

2.6.2.4. Постобработка.

Результаты расчета структурной и CFD задачи можно посмотреть в постпроцессоре.

Помните о том, что FSSOLV обновляет геометрию структуры и жидкостной области.

Соответственно вы увидите результаты обновленной геометрии. Это может привести к неверному отображению смещений. Если вам необходимо получить точное представление о структурных смещениях, возвратите исходную структуру, выполнив команду UPCOORD,-1.

Для постобработки структурных результатов выполните следующее:

1. Считайте файл структурной среды:

Команда:

PHYSICS,READ GUI:

Main Menu> Preprocessor> Physics> Environment> Read 2. Считайте результаты структурного анализа:

Команда:

SET GUI:

Main Menu> General Postproc> Read Results> Last Set 3. Если требуется точное представление смещений, выберете структурные узлы и возвратите структурную геометрию к исходному состоянию используя следующие команды или маршруты GUI:

Команда:

NSEL GUI:

Utility Menu> Select> Entities Команда:

UPCOORD,-1.

GUI:

Main Menu> Preprocessor> Loads> Load Step Opts> Other> Updt Node Coord Если необходимо вернуться к деформированной геометрии, выполните команду UPCOORD,1. Для просмотра CFD результатов выполните следующую процедуру.

Для постобработки CFD результатов выполните следующее:

1. Считайте файл жидкостной среды:

Команда:

PHYSICS,READ GUI:

Main Menu> Preprocessor> Physics> Environment> Read 2. Считайте файл результатов CFD расчета:

Команда:

SET GUI:

Main Menu> General Postproc> Read Results> Last Set 2.6.2.5. FSSOLV мультирешения.

Зачастую необходимо вычислить структурные смещения в диапазоне заданных сил, действующих со стороны жидкости. Можете использовать FSSOLV для запуска последовательных решений. В этом случае, используйте команды *DO и *ENDDO для задания цикла. Вы должны определить новые нагрузки и переписать физический файл. Для второго и последующих вызовов FSSOLV установите значение RUSEKY больше 1. Не используйте команду SAVE для сохранения базы данных между вызовами FSSOLV.

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

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

2.7. Пример анализа теплоты-напряжения, использующего косвенный метод.

Описываемый в этом разделе пример демонстрирует простой анализ тепла – напряжения, использующего косвенный метод.

2.7.1. Описание задачи.

В данном примере, два длинных тонкостенных цилиндра, расположенные вокруг одной оси, поддерживаются при температуре (Ti) на внутренней поверхности и при температуре (To) на внешней поверхности. Необходимо определить распределение температур, осевое напряжение, и окружное напряжение в цилиндрах.

Material Properties Геометрические Внутренний цилиндр Внешний цилиндр параметры Нагрузка (сталь) (аллюминий) a =.1875 in. Ti = E = 30 x 106 psi E = 10.6 x 106 psi 200°F b =.40 in. To = 70°F =.65 x 10-5 in/in°F = 1.35 x 10-5 in/in°F c =.60 in. = 0.3 = 0. K = 2.2 btu/hr-in-°F K = 10.8 btu/hr-in-°F Основная процедура косвенного анализа:

1. Определение и решение тепловой задачи.

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

3. Считывание температур из файла результатов теплового расчета.

4. Решение структурной задачи.

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

/batch,list /show /title,тепловые напряжения в концентрических цилиндрах – косвенный метод /prep et,1,plane77,,,1 ! PLANE77 осесимметричная опция mp,kxx,1,2.2 ! теплопроводность стали mp,kxx,2,10.8 ! теплопроводность аллюминия rectng,.1875,.4,0,.05! модель rectng,.4,.6,0,. aglue,all numcmp,area asel,s,area,,1 атрибуты твердой модели aatt,1,1, asel,s,area,, aatt,2,1, asel,all esize,. amesh,all ! Наложить к/э сетку nsel,s,loc,x,. d,all,temp,200 ! Задание тепловых нагружений nsel,s,loc,x,. d,all,temp, nsel.all finish /solu solve finish /post path,radial,2 ! Определение имени пути и количества точек пути ppath,1,,.1875 ! Определить путь по расположению ppath,2,,. pdef,temp,temp ! ассоциация пути с температурой pasave,radial,filea ! сохранить путь во внешний файл plpath,temp ! Отобразить температурное решение finish /prep et,1,82,,,1 ! Переключиться на структурный элемент, SOLID mp,ex,1,30e6 ! Определить структурные свойства стали mp,alpx,1,.65e- mp,nuxy,1,. mp,ex,2,10.6e6 ! Определить структурные свойства аллюминия mp,alpx,2,1.35e- mp,nuxy,2,. nsel,s,loc,y,.05 ! Задать структурные граничные условия cp,1,uy,all nsel,s,loc,x,. cp,2,ux,all nsel,s,loc,y, d,all,uy, nsel,all finish /solu tref, ldread,temp,,,,,,rth ! чтение температур из результатов теплового решения solve finish /post paresu,radial,filea pmap,,mat pdef,sx,s,x ! радиальные напряжения pdef,sz,s,z ! окружное напряжение plpath,sx,sz ! отрисовка напряжений plpagm,sx,,node finish 2.8. Пример анализа тепло – напряжение использующего физические среды Данный раздел описывает решение той же задачи, которая была рассмотрена в предыдущей главе. На этот раз используется метод физической среды. В данном случае может быть нерациональным использование метода физической среды, поскольку задача достаточно проста. Однако, это описание позволит быстро переключаться между физическими средами при последующем моделировании.

Основные процедуры метода физической среды:

1. Определить тепловую задачу.

2. Записать тепловой файл физики.

3. Удалить граничные условия и опции.

4. Определить структурную задачу.

5. Записать структурный файл физики.

6. Считать тепловой файл физики.

7. Решить и выполнить постобработку тепловой задачи.

8. Считать структурный файл физики.

9. Считать температуры из файла тепловых результатов.

10. Решить и выполнить постобработку структурной задачи.

/batch,list /show /title,тепловое напряжение в концентрических цилиндрах – метод физической среды /prep et,1,plane77,,,1 ! PLANE77 осе симметричная опция mp,kxx,1,2.2 ! Теплопроводность стали mp,kxx,2,10.8 ! Теплопроводность аллюминия rectng,.1875,.4,0,.05 ! Модель rectng,.4,.6,0,. aglue, all numcmp,area asel,s,area,,11 !Атрибуты твердой модели aatt,1,1, asel,s,area,, aatt,2,1, asel,all esize,. amesh,all ! Наложить сетку nsel,s,loc,x,. d,all,temp,200 ! Температурные нагрузки nsel,s,loc,x,. d,all,temp, nsel,all physics,write,thermal ! Записать файл тепловой физики physics,clear ! Удалить все граничные условия и опции et,1,82,,,1 ! Переключиться на структурный элемент, SOLID mp,ex,1,30e6 ! Задать структурные свойства стали mp,aplx,1,.65e- mp,nuxy,1,. mp,ex,2,10.6e6 ! Задать структурные свойства алюминия mp,alpx,2,1.35e- mp,nuxy,2,. nsel,s,loc,y,.05 ! Задать структурные граничные условия cp,1,uy,all nsel,s,loc,x,. cp,2,ux,all nsel,s,loc,y, d,all,uy, nsel,all tref, physics,write,struct ! Записать файл структурной физики save ! Решить базу данных finish /solu physics,read,thermal ! Считать файл тепловой физики solve ! Решить тепловую задачу finish /post path,radial,2 ! Задать имя пути и количество точек ppath,1,,.1875 ! Задать путь по положению ppath,2,,. pdef,temp,temp ! Накладываем температуру на путь pasave,radial,filea ! Сохранить путь во внешний файл plpath.temp ! Построить график температур finish /solu physics,read,struct ! Считать файл структурной физики ldread,temp,,,,,,rth ! Считать температуры solve ! Решить структурную задачу finish /post paresu,raidal,filea ! Восстановить путь pmap,,mat pdef,sx,s,x pdef,sz,s,z plpath,sx,sz ! Построить график напряжений plpagm,sx,,node finish Профиль напряжений Радиальные напряжения.

2.9. Пример жидкостно – структурного анализа использующего физические среды.

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

2.9.1. Постановка задачи По каналу, в котором размещено резиновое кольцо, течет вода. Скорость воды на входе в канал составляет 0.35 m/sec. Цель задачи заключается в определении потерь давления и отклонения кольца. В конце раздела приведен листинг APDL программы, позволяющей решить данную задачу посредством ANSYS.

2.9.2. Процедура.

Построить для дальнейшего анализа жидкостно-структурную модель. Для данного примера можно смоделировать три области: (а) кольцо, (b) небольшая область жидкости вокруг кольца, требующая морфинга, (с) остальная жидкая область.

Диаграмма анализируемой преграды в канале.

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

Жидкостный анализ использует нулевой тип элементов для кольца, и структурный анализ использует нулевой тип элементов для жидкости.

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

2.9.2.1. Построение модели.

Построим модель, включающую жидкую область и кольцевую область.

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

Атрибуты физической среды Набор вещественных Область Тип Материал констант Кольца 2 2 Жидкость 1 1 2.9.2.2. Создание жидкостной физической среды.

Для этого назначьте типы элементов и определите свойства материала для жидкой области.

Свойства воды задаются посредством команд FLDATA. Задаются параметры управления решением, такие как количество итераций в начальном FLOTRAN анализе. Включаем опцию турбулентности. Дальнейшее смотрите в листинге программы.

Жидкая физическая среда Область Тип Материал Набор вещественных констант Кольцо Нулевой тип (0) Отсутствует Отсутствует Жидкость FLUID141 Вязкость, плотность Отсутствует • Задайте соответствующие номинальные граничные условия и нагрузки на жидкую область.

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

(кольцо) Команда CM,GASKET,NODES GUI:

Utility Menu> Select> Comp/Ass'y> Create Component • Запишите физическую среду жидкости в файл.

Команда:

PHYSICS,WRITE,FLUID,FLUID GUI:

Main Menu> Preprocessor> Physics> Environment> Write 2.9.2.3. Создание структурной физической среды • Удалите всю информацию определенную для жидкой среды с целью подготовки модели для структурной среды.

Команда:

PHYSICS,CLEAR GUI:

Main Menu> Preprocessor> Physics> Environment> Clear • Измените тип элементов областей с жидкого на структурный. FLUID141 должен стать PLANE42. Установите нулевой тип элементов для жидкой области, поскольку эта область не требуется для структурной физической среды.

Структурная физическая среда Вещественная Область Тип Материал константа Кольцо HYPER56 Mooney-Rivlin нет Жидкость Нулевой тип нет нет • Установите структурные свойства для каждой области, требуемой для структурного анализа.

• Задайте структурные граничные условия.

Номинальные граничные условия структурной физики • Определите шаг нагружения и опции решения.

• Запишите в файл структурную физическую среду (например., PHYSICS,WRITE,STRUC,STRUC).

2.9.2.4. Цикл решения жидкость/структура В интерактивном или в пакетном режиме (смотрите листинг программы), вводится цикл решения жидкость-структура. В этом случае, максимальное смещение кольца (maximum gasket displacement (MGD)) становится общим монитором сходимости. Задача считается решенной в том случае, если изменение MGD между двумя последовательными структурными решениями меньше величины сходимости.

Начальный анализ во FLOTRAN должен решается до тех пор пока не сойдется.

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

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

Этапы, составляющие цикл решения:

1. Чтение жидкой физической среды.

Команда:

PHYSICS,READ,fluid GUI:

Main Menu> Solution> Physics> Environment> Read 2. По необходимости установите требуемый параметр FLOTRAN (например количество глобальных итераций) Команда:

FLDATA2,ITER,EXEC, GUI:

Main Menu> Solution> FLOTRAN Setup> Execution Control 3. Решение во FLOTRAN.

Команда:

SOLVE GUI:

Main Menu> Solution> Run FLOTRAN 4. Чтение структурной среды.

Команда:

PHYSICS,READ,struc GUI:

Main Menu> Solution> Physics> Environment> Read 5. При необходимости используйте команду /ASSIGN для рестарта структурного решения.

Команда:

/ASSIGN,esave,struc,esav /ASSIGN,emat,struc,emat GUI:

Utility Menu> File> ANSYS File Options 6. Возвратите узлы в исходное положение для последующего структурного анализа и будущего морфинга.

Не выполняйте этот шаг в первом цикле решения типа жидкость-структура.

Команда:

PARSAV,ALL GUI:

Utility Menu> Parameters> Save Parameters Команда:

RESUME GUI:

Utility Menu> File> Resume Jobname.db Команда:

PARRES GUI:

Utility Menu> Parameters> Restore Parameters 7. Рестарт анализа.

Не выполняйте этот шаг в первом цикле решения типа жидкость-структура.

Команда:

ANTYPE,STATIC,REST GUI:

Main Menu> Solution> Restart 8. Выберете узлы/элементы на которые будут заданы давления полученные в результате решения задачи во FLOTRAN.

9. Выполните команду LDREAD.

Команда:

LDREAD,PRES,LAST,,,,,rfl GUI:

Main Menu> Solution> Define Loads> Apply> Structural> Pressure> From Fluid Analy 10. Выберете опцию не использовать мультифреймовые файлы рестарта.

Команда:

RESCONTROL,,NONE GUI:

Main Menu> Solution> Nonlinear> Restart Control 11. Решите структурную задачу и сохраните базу данных.

Команда:

SOLVE GUI:

Main Menu> Solution> Solve Команда:

SAVE GUI:

Utility Menu> File> Save as Jobname.db 12. Выполните морфинг сетки в небольшой области жидкости вокруг кольца (имя компонента AREA2).

Команда:

DAMORPH,AREA2,, GUI:

Main Menu> Preprocessor> Meshing> Modify Mesh> Refine At> Areas 13. Сравните величину сдвига сетки с последним значением. Выберете именованный компонент GASKET и выведете список узловых координат.

14. Проверьте сходимость сравнив максимальные значения последовательных смещений кольца.

15. Просмотрите графики в файле gasket.grph.

2.9.3. Результаты Цикл решения жидкость-структура выполнялся до тех пор, пока не была достигнута сходимость. Допуск сходимости составлял 0.5%. Для первого анализа во FLOTRAN, оказалось достаточно 400 глобальных итераций для достижения сходимости. В цикле жидкость структура, количество итераций для оставшегося решения во FLOTRAN установлено равным 100.

На нижеприведенных рисунках показаны линии потока и контуры давления. Качественно результаты выглядят идентично для недеформированного (первого анализа) и деформированного (окончательного анализа) случаев.

Линии потока возле кольца Контуры давления Также на нижеприведенном рисунке показаны нагрузки в кольце, полученные в результате последнего анализа. Максимальная нагрузка в последнем анализа приблизительно на 25% меньше, чем максимальная нагрузка первого анализа. Это указывает на то, что рассмотрение эффекта смещенной геометрии на поле потока показало существенное отличие результатов.

Перед построением графиков структурных результатов структурные узлы должны быть возвращены в их исходное положение.

Команда:

UPCOORD,- GUI:

Main Menu> Solution> Load Step Opts> Other> Updt Node Coord Напряжения в кольце /Batch,list /prep /sho,gasket,grph shpp,off ET,1,141 ! Жидкость – статичная сетка ET,2,56, ! Гиперэластичный элемент !!!!!!! Мультифизика – взаимодействие жидкости со структурой !!!!!!! Деформация кольца в поле потока.

!

!!!!!!! Графики записываются в файл gasket.grph.

!

! - Вода течет в вертикальной трубе через резиновое кольцо ! - Определить равновесное положение кольца и результирующее поле потока !

! | | ! | | ! |----------| Граница "области морфинга жидкости" ! | | ! | | Кольцо ! | | ! |----------| Граница "области морфинга жидкости" (sf) ! | | !

!! 1. Полностью построить модель:

!! Жидкая область – статичная сетка !!

!! Кольцо оставляет отверстие в центре трубы !! Деформируемая жидкая область – определяемая пользователем область вокруг !! кольца. Сетка жидкости здесь будет деформирована и обновлена !! по мере деформации кольца.

!!

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

yent = 0.0 ! Y координата входа в трубу dyen = 1.0 ! Длина недеформируемой геометрии ysf1 = yent+dyen ! Y координата области входа в деформируемую жидкость dsf1 = 0. ygas = ysf1+dsf1 ! Y координата основания кольца dg = 0.02 ! толщина кольца dg2=dg/2.

ytg = ygas+dg ! Y координата исходного положения верхней части кольца dsf2 = 0.5 ! толщина области нисходящего потока ysf2 = ytg + dsf2! Y верха морфируемой жидкой области нисходящего потока dyex = 6.0 ! длина выхода жидкости x = 0. ! ось симметрии dgasr =.20 ! начальный промежуток кольца piper = 0.3 ! радиус трубы xrgap = piper-dgasr !!

!!! Создание геометрии !!

rect,xrgap,piper,ygas,ytg ! A1:Кольцо (ключевые точки 1-4) rect,x,piper,ysf1,ysf2 ! A2: Морфируемая жидкая область rect,x,piper,yent,ysf1 ! A3: Жидкая область со статичной сеткой rect,x,piper,ysf2,ysf2+dyex ! A4: Жидкая сетка со статичной областью aovlap,all k,22,xrgap+dg2,ygas+dg rarc = dg2*1. larc,1,4,22,rarc al,6, adelete, al,6,3,22,7,8,5,21, !! Данные разбиения сетки ngap = 10 ! Число элементов вдоль промежутка ngas = 10 ! Число элементов вдоль кольца rgas = -2 ! Коэффициент разбивка вдоль кольца nflu = ngap+ngas ! Количество элементов вдоль жидкой области raflu = -3 ! Элементы жидкости возле стен и центра nenty =8 ! Элементы вдоль потока - вход raent =5 ! Коэффициент разбивки вдоль входной области nfl1 = 20 ! Элементы вдоль потока – первая морфируемая жидкость.

nthgas = 4 ! Элементы кольца nfl2 = 3 ! Элементы вдоль потока – вторая морфируемая жидкость.

next = 30 ! Элементы вдоль потока – выходная область rext = rafls = 12 ! Начальный коэффициент разбивки элементов - морф.жидкость lesize,1,,,ngas,rgas lesize,3,,,ngas,rgas nfl11= nfl1*2+ lsel,s,,,2,4,2 ! (Модифицирует lesize линии 8 при изменении сетки кольца) lesize,all,,,nthgas alls lesize,5,,,nflu,raflu lesize,7,,,nflu,raflu lesize,9,,,nflu,raflu lesize,15,,,nflu,raflu lesize,18,,,nenty,1./raent lesize,17,,,nenty,1./raent lesize,21,,,nfl1,rafls lesize,8,,,nfl11,-1./(rafls+3) lesize,22,,,nfl1,rafls lesize,19,,,next,rext lesize,20,,,next,rext !!! AATT,MAT,REAL,TYPE - Установка атрибутов площадей asel,s,,,1, aatt,2,2,2 ! Кольцо (материал 2) asel,s,,, cm,area2,area alist ! Список площадей, выбранных для дальнейшего морфинга asel,a,,,5, aatt,1,1,1 ! Жидкая площадь (материал 1) alls eshape, asel,u,,,2, amesh,all eshape, asel,s,,,2, amesh,all !----------------- !!!!! Создание элементного графика и запись в gasket.grph asel,s,,,1, esla,s /Title, Начальная сетка для кольца и соседних с ним областей eplot /ZOOM,1,RECT,0.3,-0.6,0.4,-0. alls !----------------- !!

!!!!! 2. Создание физической среды для жидкости et,1, et,2,0 ! Кольцо становится нулевым элементом vin=3.5e-1 ! Входная скорость воды (meters/second) !! Опции управлением CFD решением flda,solu,flow, flda,solu,turb, flda,iter,exec, flda,outp,sumf, !! Данные CFD свойств flda,prot,dens,constant flda,prot,visc,constant flda,nomi,dens,1000. ! плотность воды 1000 kg/m flda,nomi,visc,4.6E-4 ! 4.6E-4 kg-s/m (вязкость воды) flda,conv,pres,1.E-8 ! Величина сходимости уравнения давления !! CFD Граничные условия (Задаются на твердую модель) lsel,s,,,8,17, lsel,a,,, dl,all,,vx,0.,1 ! симметрия средней линии lsel,s,,, dl,all,,vx,0., dl,all,,vy,vin,1 ! входное условие lsel,s,,, lsel,a,,,18, lsel,a,,,21, dl,all,,vx,0.,1 ! внешняя стенка dl,all,,vy,0., lsel,s,,,1,3, lsel,a,,, dl,all,,vx,0.,1 ! кольцо dl,all,,vy,0., lsel,s,,, dl,15,,pres,0.,1 ! давление на выходе !!! создадим именованный компонент узлов нижней части кольца lsel,s,,, nsll,, cm,gasket,node nlist ! список начальный узловых координат основания кольца /com, +++++++++ исходные координаты кольца -------- alls /title, жидкостный анализ physics,write,fluid,fluid !!

!!!!! 3. создадим структурную физическую среду !!

physics,clear et,1,0 ! нулевой элемент для жидкой области et,2,56 ! элемент кольца - материал mp,ex,2,2.82E+6 ! модуль Юнга для резины mp,nuxy,2,0.49967 ! коэффициент Пуассона для резины tb,mooney, tbdata,1,0.293E+6 ! константы Mooney-Rivlin tbdata,2,0.177E+6 ! " " " lsel,s,,, nsll,, d,all,ux,0.

d,all,uy,0. ! зафиксируем конец кольца alls /title,структурный анализ finish /solu antype,static nlgeom,on cnvtol,f,,,,- physics,write,struc,struc physics,clear save !!

!!!!! 4. цикл взаимодействия жидкость-структура !!

loop=25 ! максимально допустимое количество циклов toler=0.005 ! допуск сходимости для максимального смещения *dim,dismax,array,loop ! определим массив величин максимального смещения *dim,strcri,array,loop ! определим массив значений сходимости *dim,index,array,loop *do,i,1,loop ! запуск решений жидкость -> структура /solu physics,read,fluid ! чтение жидкой среды *if,i,ne,1,then flda,iter,exec,100 ! выполнить 100 глобальных итераций для *endif ! каждой новой геометрии solve ! решить во FLOTRAN fini ! end of fluid portion physics,read,struc ! чтение структурной среды /assign,esave,struc,esav ! файлы для нелинейного структурного рестарта /assign,emat,struc,emat *if,i,gt,1,then ! цикл структурного рестарта parsave,all ! сохранение параметров для проверки сходимости resume ! восстановить базу данных для того чтобы восстановить ! исходное положение узлов parresume ! восстановить параметры необходимые для проверки сходимости /prep antype,stat,rest fini *endif /solu solc,off lsel,s,,,1,3,2 ! выбор соответствующих линий для задания давлений lsel,a,,,6 ! на поверхность кольца nsll,, esel,s,type,, ldread,pres,last,,,,,rfl ! задание поверхностных нагрузок (давлений) из Flotran alls rescontrol,,none ! не использовать мультифреймовый рестарт solve *if,i,eq,1,then save ! сохранить исходные координаты узлов при первом запуске *endif fini /post cmsel,s,gasket nsort,u,sum,1, *get,dismax(i),sort,0,max ! получить максимальную величину смещения strcri(i)=toler*dismax(i) alls fini /prep mkey=2 ! выбор уровня морфинга сетки жидкости damorph,area2,,mkey ! морфинг жидкой области !---------------- !!!!! Построение графика и запись его в файл gasket.grph fini /prep et,1, asel,s,,,1, esla,s /Title, EPLOT после DAMORPH,area2,,%mkey% шаг номер %i% eplot alls !----------------- cmsel,s,gasket nlist ! список обновленных координат основания кольца для проверки /com, +++++++++ обновленные координаты кольца -------- alls fini /assign,esav /assign,emat !!!! проверка критерия сходимости imax= i index(i)=i *if,i,gt,1,then strcri(i)=abs(dismax(i)-dismax(i-1))-toler*dismax(i-1) *if,strcri(i),le,0,then strcri(i)= *exit ! остановить цикл при достижении сходимости *endif *endif *enddo !!!!! конец расчетного цикла save ! сохранение узловых координат деформированной геометрии !!!!! печать сходиомости *vwrite (/'номер цикла. Макс-ое сиещение структурная сходимость') /nopr *vlen,imax *vwrite,index(1),dismax(1),strcri(1) (f7.0,2e17.4) finish !!!!! постобработка результатов !!! 1. результаты Flotran.

physics,read,fluid /post set,last /Title, Flotran: линии течения возле кольца plnsol,strm /Title, Flotran: контуры давления plnsol,pres fini !!! 2. структурные результаты.

physics,read,struc /post set,last upcoord,-1 ! Возврат к исходному положению узлов /Title, структурные результаты: напряжения von Mises plnsol,s,eqv,1, fini !/exit,nosave 2.10. Пример анализа индукционного нагрева использующего физические среды.

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

2.10.1. Описание задачи.

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

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

Осесимметричная одномерная полоска области индукционного нагрева.

2.10.2. Процедура.

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

Блок схема решения Рассмотрим процедуру решения задачи индукционного нагрева.

2.10.2.1. Шаг 1: Задание свойств Свойства физической среды Область Тип Материал Набор вещественных констант Брус 1 2 Катушка 2 3 Область Тип Материал Набор вещественных констант Воздух 2 1 Поверхность бруса 3 2 2.10.2.2. Шаг 2: Построение модели Постройте модель. Задайте свойства материалов. Поверхность бруса будет использована для задания поверхностного элемента симулирующего лучистый теплообмен.

2.10.2.3. Шаг 3: Создание электромагнитной среды.

Создайте электромагнитную физическую среду определив тип элемента и свойства материалов, как показано в нижеприведенной таблице.

Электромагнитная физическая среда.

Область Тип Материал Вещественная константа Брус PLANE13 MURX(T), RSVX(T) Нет Катушка PLANE13 MURX Нет Воздух PLANE13 MURX Нет Поверхность бруса Нулевой Type (0) Нет Нет • Задайте номинальные граничные условия и нагрузки.

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

• Запишите файл электромагнитной физической среды.

Команда:

PHYSICS,WRITE GUI:

Main Menu> Preprocessor> Physics> Environment 2.10.2.4. Шаг 4: Создайте тепловую физическую среду.

Создайте тепловую физическую среду:

• Удалите номинальные граничные условия и повторно выставьте опции.

Команда:

PHYSICS,CLEAR GUI:

Main Menu> Preprocessor> Physics> Environment • Смените тип элементов с электромагнитного на тепловой так же как и KEYOPT опции. Задайте нулевой тип элемента для воздуха и катушки (в данном случае в тепловом анализе рассматривается брус).

Тепловая физическая среда.

Область Тип Материал Вещественная конст KXX(T), Брус PLANE55 Нет ENTH(T) Нулевой тип Type Катушка Нет Нет (0) Нулевой тип Type Воздух Нет e Нет (0) Поверхность Постоянная Стефана SURF151 EMIS бруса Больцмана • Задайте тепловые свойства и вещественные константы.

• Установите соответствующие номинальные граничные условия и нагрузки, как показано ниже.

Номинальные тепловые граничные условия.

• Задайте соответствующие опции шага нагружения и опции решения.

Команда:

PHYSICS,WRITE GUI:

Main Menu> Preprocessor> Physics> Environment> Write Main Menu> Solution> Physics> Environment> Write 2.10.2.5. Шаг 5: Запись тепловой физической среды.

Сохраните файл тепловой физической среды.

Команда:

PHYSICS,READ GUI:

Main Menu> Solution> Physics> EnvironmentMain Menu> Solution> Physics> Environment> Read 2.10.2.6. Шаг 6: Подготовка DO цикла Подготовьте *DO цикл, который будет итерировать между гармоническим электромагнитным и нестационарным тепловым анализом.

Считайте электромагнитную физическую среду и решите гармонический анализ.

Команда:

SOLVE GUI:

Main Menu> Solution> Solve> Current LS • Считайте температуры для вычисления свойств материалов.

Команда:

LDREAD GUI:

Main Menu> Preprocessor> Define Loads> Apply> Structural> Temperature> From Therm Analy • Для того, чтобы впоследствии использовать файлы ESAV и EMAT в нестационарном тепловом анализе используйте команду /ASSIGN.

Команда:

/ASSIGN GUI:

Utility Menu> File> ANSYS File Options • Считайте тепловую физическую среду.

Команда:

PHYSICS,READ GUI:

Main Menu> Preprocessor> Physics> Environment> Read • Считайте энерговыделения из электромагнитного анализа.

Команда:

LDREAD GUI:

Main Menu> Preprocessor> Define Loads> Apply> Thermal> Heat Generat> From Mag Analy • Решите нестационарный тепловой анализ для заданного приращения времени.

• Возврат файла к исходному состоянию.

Команда:

/ASSIGN GUI:

Utility Menu> File> ANSYS File Options 2.10.2.7. Шаг 7: Повтор предыдущего шага.

Повторит предыдущий шаг для следующего временного интервала t.

2.10.2.8. Шаг 8: Постобработка результатов.

Обработка результатов выполненного расчета.

2.10.2.9. Листинг программы.

Текст после восклицательного знака является комментарием.

/batch,list /filenam,induc /prep shpp,off /title, индукционный нагрев цилиндрического стержня /com, et,1,13,,,1 ! PLANE13, осе симметричный, степень свободы AZ et,2,13,,, et,3,151,,,1,1,1 ! SURF151, тепловой, лучистый r,3,0 ! Набор вещественных констант для SURF row=.015 ! внешний радиус ric=.0175 ! внутренний радиус катушки roc=.0200 ! внешний радиус катушки ro=.05 ! внешний радиус модели t=.001 ! толщина модели freq=150000 ! частота (Гц.) pi=4*atan(1) ! число пи cond=.392e7 ! максимальная проводимость muzero=4e-7*pi ! проницаемость среды mur=200 ! максимальная относительная проницаемость skind=sqrt(1/(pi*freq*cond*muzero*mur)) ! глубина проникновения поля ftime=3 ! финальное время tinc=.05 ! шаг временного приращения для гармонического анализа time=0 ! начальное время delt=.01 !

emunit,mks ! единицы измерения в магнитном анализе mp,murx,1,1 ! относительная проницаемость воздуха mp,murx,3,1 ! относительная проницаемость катушки mptemp,1,25.5,160,291.5,477.6,635,698 ! температуры для относительной прониц ти mptemp,7,709,720.3,742,761, mpdata,murx,2,1,200,190,182,161,135,104 ! относительная проницаемость стали mpdata,murx,2,7,84,35,17,1, mptemp mptemp,1,0,125,250,375,500,625 ! температуры для сопротивления mptemp,7,750,875, mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e- mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6 ! сопротивления стали rectng,0,row,0,t ! брус rectng,row,ric,0,t ! воздушный промежуток rectng,ric,roc,0,t ! катушка rectng,roc,ro,0,t ! окружающая среда (воздух) aglue,all numcmp,area ksel,s,loc,x,row ! выбираем ключевые точки по внешнему радиусу kesize,all,skind/2 ! задаем размер к/э сетки в половину глубины проникновения ! поля.

ksel,s,loc,x,0 ! выбираем ключевые точки в центре kesize,all,40*skind ! задаем размер сетки lsel,s,loc,y,t/2 ! выбираем вертикальные линии lesize,all,,,1 ! задаем 1 разбиение в глубину lsel,all asel,s,area,, aatt,2,1,1 ! устанавливаем атрибуты бруска asel,s,area,, aatt,3,1,2 ! задаем атрибуты катушки asel,s,area,,2,4, aatt,1,1,2 ! атрибуты воздуха asel,all mshape,0,2d mshk, amesh,1 ! накладываем сетку на брус lsel,s,loc,y, lsel,a,loc,y,t lsel,u,loc,x,row/ lesize,all,. lsel,all amesh,all ! накладываем сетку на оставшиеся площади n ! создаем пространственный узел для SURF *get,nmax,node,,num,max lsel,s,loc,x,row type, real, mat, lmesh,all *get,emax,elem,,num,max emodif,emax,3,nmax ! модифицируем элемент для того, чтобы добавить ! пространственный узел для учета лучистого теплообмена et,3,0 ! установка 3 типа в нулевой элемент nsel,s,loc,x d,all,az,0 ! задание граничных условий.

nsel,all esel,s,mat,, bfe,all,js,,,,15e6 ! плотность тока в катушке esel,all finish /solu antyp,harm harfrq, physics,write,emag ! запись файла электромагнитной среды finish /prep lsclear,all ! очистка всех граничных условий и опций et,1,55,,,1 ! осесимметричный тепловой элемент PLANE et,2,0 ! нулевой тип элемента для катушки и воздуха et,3,151,,,1,1,1 ! радиационный элемент SURF keyopt,3,9, r,3,1,5.67e-8 ! угловой коэффициент, постоянная Стефана-Больцмана mptemp mptemp,1,0,730,930,1000 ! температуры для теплопроводности mpdata,kxx,2,1,60.64,29.5,28, mptemp ! температуры для энтальпии mptemp,1,0,27,127,327,527, mptemp,7,765,765.001, mpdata,enth,2,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e mpdata,enth,2,7,3.548547e9,3.548556e9,4.3520e mp,emis,2,.68 ! степень черноты finish /solu antype,trans toffst, tunif,100 ! начальная температура модели d,nmax,temp,25 ! температура внешней среды cnvtol,heat,1 ! допуск сходимости kbc,1 ! пошаговые нагрузки trnopt,full autos,on ! автоматический выбор временного шага deltim,1e-5,1e-6,delt,on outres,basic,all ! сохранить всю информацию по шагу нагружения physics,write,thermal ! записать файл тепловой среды finish *do,i,1,ftime/tinc ! цикл решения time=time+tinc ! приращение времени physics,read,emag ! считывание файла э/м среды /solu *if,i,eq,1,then tunif,100 ! начальная температура *else ldread,temp,last,,,,,rth ! считать температуры из теплового анализа *endif solve ! решить гармонический анализ finish physics,read,thermal ! считать файл тепловой среды /assign,esav,therm,esav ! переключить файлы для использования в тепловом ! рестарте /assign,emat,therm,emat /solu *if,i,gt,1,then antype,trans,rest ! тепловой рестарт *endif time,time ! время в конце теплового анализа esel,s,mat,,2 ! выбрать брусок ldread,hgen,,,,2,,rmg ! задать энерговыделение, полученное из э/м анализа esel,all solve finish /assign,esav ! привести файлы к исходному состоянию /assign,emat *enddo ! конец цикла решения finish save ! сохранить базу данных /post26 ! постобработчик /show nsol,2,1,temp,,tempcl ! сохранить температуры на центральной линии бруса nsol,3,2,temp,,tempsurf ! сохранить температуры по радиусу бруса plvar,2,3 ! построить график зависимости тем-ры от времени prvar,2,3 ! распечатать температуры finish 2.10.2.10. Результаты Зависимость температуры бруска от времени 2.11. Пример структурно – электростатического анализа.

Пример, описанный в этом разделе демонстрирует структурно-электростатический анализ, выполненный при помощи макроса ESSOLV.

2.11.1. Постановка задачи.

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

Кремниевый брусок над заземленным затвором.

2.11.2. Построение и решение модели (командный метод) Текст после восклицательного знака является комментарием.

/batch,list /prep7, Silicon beam deflection from an applied voltage et,1,121 ! Элемент PLANE121 для воздуха et,2,121 ! Временный элемент для бруса emunit,epzro,8.854e-6 ! Проницаемость среды, единицы µMKSV mp,perx,2,1 ! Относительная диэлектрическая проницаемость для ! воздуха bl=150 ! Длина бруса (µm) bh=2.0 ! Высота бруса glc=bl/2 ! Центральное расположение заземления gl=90 ! Длина проводящей зоны заземления gh=1.5 ! Толщина затвора gap=4.5 ! Воздушный промежуток vltg=120 ! Напряжение rectng,0,bl,gap,bh+gap ! Создание модели rectng,glc-gl/2,glc+gl/2,-gh, rectng,-10,170,-20, aovlap,all asel,s,area,,1 ! Площадь бруса aatt,1,, asel,s,area,,4 ! Площадь для воздуха cm,air,area ! Сгруппировать площади воздуха в элемент aatt,2,, allsel,all smrtsiz, amesh,1 ! Мешировать брус mshape, amesh,4 ! Мешировать воздух треугольниками asel,s,area,, lsla,s dl,all,,volt,vltg ! Задать напряжение на брус asel,s,area,, lsla,s dl,all,,volt,0 ! Заземленный проводник (не мешированый) allsel,all et,1,0 ! Задать структуре нулевой тип элемента physics,write,ELECTROS ! Записать файл электростатической физики physics,clear ! Очистить физику et,1,82,,,2 ! Задание типа элементов для бруска et,2,0 ! Для воздуха нулевой тип элементов mp,ex,1,170e3 ! Установка коэффициента µN/(µm)** mp,nuxy,1,0. dl,4,,ux,0 ! Ограничивающие условия бруса dl,4,,uy, allsel,all finish physics,write,STRUCTURE ! Запись структурного физического файла ESSOLV,'ELECTROS','STRUCTURE',2,0,'air',,,,10 ! Решение сопряженной задачи finish physics,read,ELECTROS ! Чтение файла электростатической физики /post set,last ! Извлечение резульатов esel,s,mat,,2 ! Выбор элементов воздуха etable,fx,fmag,x ! Извлечение электростатических сил etable,fy,fmag,y ssum ! Сумма сил finish 2.12. Пример жидкостно-структурного анализа, использующего FSSOLV Пример, описанный в этом разделе демонстрирует жидкостно-структурный анализ, выполненный при помощи макроса FSSOLV.

2.12.1. Постановка задачи.

В этом примере тонкий брус расположен внутри канала по которому течет вода. Разница давления между верхом и низом бруса вызывают колебания бруска.

Тонкий алюминиевый брус, обтекаемый водой.

2.12.2. Создание и решение модели (командный метод).

Текст за восклицательным знаком является комментарием.

/BATCH,list /Title, Тонкий алюминиевый брусок, обтекаемый водой.

!!

!!!!! 0. Определение значений параметров.

!!

fact=1.0e-3 ! переводной коэффициент length=.75*fact ! длина бруска d0 =.005*fact ! толщина бруска d1 =.04*fact ! поперечный размер нижней части канала d2 =.02*fact ! поперечный размер верхней части канала dcen=d1+0.5*d0 ! центральное расположение стержня widt=d0+d1+d2 ! общая ширина модели dl =widt elsi=0.5*d0 ! размер элемента для меширования modx=70.0e9 ! модуль Юнга nust=0.3 ! коэффицинт Пуассона denw=1.0e3 ! плотность воды visw=1.0e-3 ! вязкость воды vin =8.85e-3 ! скорость воды на входе (метры/секунда) /PREP !!

!!!!! 1. создание геометрии и наложение конечно элементной сетки !!

et,1,141 ! двумерный Flotran элемент для симуляции жидкости et,2,141 ! временный элемент для бруса rect,,length+dl,,widt rect,,length,d1,d1+d aovlap,all numcmp,area ASEL,S,loc,y,dcen ASEL,INVERT AATT,1,, type, ASEL,INVERT aatt,2,, type, alls esize,elsi amesh,all alls fini /prep !!

!!!!! 2. создание физической среды для жидкости !!

et,1, et,2,0 ! назначить нулевой тип элемента для бруса !! средства управления CFD решением flda,solu,flow, flda,solu,turb, flda,advm,mome,supg flda,advm,turb,supg flda,iter,exec, flda,outp,sumf, !! информация по CFD свойствам flda,prot,dens,constant flda,prot,visc,constant flda,nomi,dens,denw ! 1000 kg/m3 плотность воды flda,nomi,visc,visw ! 4.6E-4 kg-s/m (вязкость воды) flda,conv,pres,1.E-9 ! величина сходимости для уравнений давления !! CFD граничные условия lsel,s,loc,x, lsel,u,loc,y,dcen dl,all,,vy,0., dl,all,,vx,vin,1 ! входное условие lsel,s,loc,y, lsel,a,loc,y,widt dl,all,,vx,0.,1 ! внешние стенки dl,all,,vy,0., asel,s,type,, lsla,s dl,all,,vx,0.,1 ! структурные граничные условия dl,all,,vy,0., lsel,s,loc,x,length+dl dl,all,,pres,0.,1 ! давление на выходе asel,s,type,, cm,area2,area ! компонент жидкости для морфинга alls fini /title,Анализ жидкости physics,write,fluid,fluid physics,clear /prep !!

!!!!! 3. создание физической среды для структуры !!

et,1, et,2,42 ! структурный тип элемента для бруска mp,ex,2,modx mp,nuxy,2,nust LSEL,S,LOC,Y,dcen LSEL,R,LOC,X, dl,all,,UX,0 ! ограничивающее условие для бруска при X= dl,all,,UY, alls fini /title,Structural Analysis physics,write,struc,struc physics,clear !!

!!!!! 4. Запуск макроса FSSOLV с рестартом и принятым по умолчанию допуском !!!!! сходимости !!

FSSOLV,'fluid','struc',2,'area2',,,,,25,40,1, /exit,nosa Глава 3. Последовательный слабосопрояженный анализ.

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

Последовательный слабо сопряженный метод использует элементы ANSYS/FLOTRAN для жидкости и структурные, тепловые или сопряженные элементы для твердых частей модели.

Обычно задачи в которых рассматривается взаимодействие жидкость – твердое тело включают следующее:

• Задачи, которые включают в себя взаимодействие тепловых, структурных или электрических полей с полем течения.

• Задачи, которые включают в себя течение, вызванное структурными изменениями.

• Задачи, описывающие течение, вызванное деформацией твердого тела.

Ниже приведен алгоритм слабосопряженного метода решения задач взаимодействия жидкости с твердым телом.

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

3.2. Граница твердо жидкого раздела.

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

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

Примечание.

На границе раздела твердой части модели с жидкостью необходимо иметь два набора составляющих твердой модели и узлов: один набор для жидкой области и второй набор для твердой части модели.

Параметры передаваемые через границу раздела.

Передача нагрузки со стороны Передача нагрузки со стороны Интерполяционный жидкости к твердой части твердой части модели к метод модели жидкости • Сила по трем направлениям • Смещения координат • Скорость к/э сетки Консервативный • Тепловая мощность • Температуры • Смещения • Плотность силы по трем • Скорость к/э сетки направлениям координат Неконсервативный • Плотность теплового потока • Температуры 3.3. Элементы жидкости и твердой части модели.

Анализ взаимодействия жидкости с твердой частью модели использует элементы жидкости ANSYS/FLOTRAN: FLUID141 или FLUID142. Дополнительную информацию по элементам FLOTRAN и их свойствам смотрите в ANSYS Elements Reference и ANSYS CFD Flotran Analysis Guide.

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

Структурные элементы двумерные трехмерные балочные оболочечные HYPER PLANE2 SOLID45 BEAM3 SHELL63 HYPER PLANE42 SOLID92 BEAM23 SHELL93 HYPER PLANE82 SOLID95 BEAM188 SHELL181 HYPER PLANE182 SOLID185 BEAM PLANE183 SOLID SOLID Тепловые элементы двумерные трехмерные оболочечные PLANE35 SOLID70 SHELL PLANE55 SOLID SOLID Сопряженные элементы двумерные трехмерные оболочечные PLANE13 SOLID5 SHELL PLANE67 SOLID SOLID SOLID 3.4. Анализ взаимодействия жидкости с твердой частью модели.

Процедура выполнения анализа взаимодействия жидкости с твердой частью модели состоит из пяти основных этапов:

1. Настройка жидкостного анализа и анализа твердой части модели.

2. Отметка поверхности твердо жидкого раздела.

3. Определение опции решения твердо жидкого анализа.

4. Получение решения.

5. Постобработка результатов.

3.4.1. Настройка жидкостного анализа и анализа твердой части модели.

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

Информацию по настройке анализа жидкости смотрите в ANSYS CFD Flotran Analysis Guide.

Информацию по настройке анализа твердой модели смотрите в ANSYS Structural Analysis Guide и ANSYS Thermal Analysis Guide.

3.4.2. Отметка поверхности твердо жидкого раздела.

Следующий шаг заключается в отметке определенным номером поверхностей твердо жидкого раздела где передается нагрузка. Для этого используется семейство SF команд (SF, SFA, SFE, или SFL) и метка поверхностной нагрузки FSIN. Установите флаг раздела твердо жидких поверхностей дважды: для стороны жидкости и для стороны твердой части модели.

При выполнении команд SF, SFA, SFE, или SFL используйте один и тот же номер раздела для поверхностей жидкости и твердого тела, где передается нагрузка. Нагрузка передается между поверхностями жидкости и твердой части модели, имеющих один и тот же номер раздела.

3.4.3. Задание опций решения твердо жидкого анализа.

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

3.4.3.1. Установка основных опций анализа.

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

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

Команда FSAN,OFF отключает эту опцию (по умолчанию опция анализа твердо жидкого взаимодействия выключена).

Используйте команду FSOR для установки последовательности анализа жидкости и твердой части модели. Для того, чтобы выполнить первым анализ жидкости выполните команду со следующими аргументами FSOR,FLUID,SOLID. По умолчанию анализ жидкости следует за анализом твердой части модели. Физика задачи определяет порядок решения. Наиболее существенно влияющий на решение анализ должен быть выполнен в первую очередь.

Используйте команду FSTR для задания типа анализа (стационарный, нестационарный). По умолчанию задан нестационарный анализ. Команда FSTR имеет более высокий приоритет, чем команды ANTYPE или FLDATA1.

Команда FSIN используется для определения метода интерполяции для передачи нагрузки через границу твердо жидкого раздела.

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

Консервативная интерполяция:

• Элементы низкого порядка, как на стороне жидкости, так и на стороне твердой части модели.

• Количество узлов на стороне жидкой части модели превышает количество узлов на стороне твердой части модели.

Неконсервативная интерполяция:

• Элементы высокого порядка на стороне твердой части модели.

• Количество узлов на стороне твердой части модели превышает количество узлов на стороне жидкости.

GUI эквиваленты приведенных команд:

GUI:

Main Menu>Preprocessor>FSI Set Up>Options Main Menu>Solution>FSI Set Up>Options 3.4.3.2. Задание времени окончания анализа и временного шага.

Необходимо определить время окончания анализа и значение временного шага для анализа твердо жидкого взаимодействия. Для задания времени окончания анализа используйте команду FSTI. Значение времени окончания анализа должно быть кратно величине временного шага. Последовательный слабосопряженный метод поддерживает только постоянный шаг приращения времени. По умолчанию приращение временного шага и время окончания анализа равны единице.

Для задания времени окончания анализа используйте следующее:

Команда:

FSTI GUI:

Main Menu>Preprocessor>FSI Set Up>Time Main Menu>Solution>FSI Set Up>Time Для задания шага временного приращения:

Команда:

FSDT GUI:

Main Menu>Preprocessor>FSI Set Up>Time Main Menu>Solution>FSI Set Up>Time Вы должны также задать шаг временного приращения для анализа твердой и жидкой частей модели используя команды DELTIM and FLDATA4,TIME,STEP,Value соответственно. Шаг временного приращения для анализа твердой и жидкой частей модели должен быть меньше или равен шагу временного приращения для твердо жидкостного анализа. Анализ твердо жидкого взаимодействия итерирует между анализом жидкости и твердого тела так, что нагрузка передается через границу твердо жидкого раздела через приращения времени, заданное командой FSDT.

3.4.3.3. Задание максимального количества итераций.

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

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

Команда:

FSIT GUI:

Main Menu>Preprocessor>FSI Set Up>Iterations Main Menu>Solution>FSI Set Up>Iterations По умолчанию выставлено 10 итераций.

3.4.3.4. Определение величин сходимости.

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

Для задания нормы сходимости используйте следующее:

Команда:

FSCO GUI:

Main Menu>Preprocessor>FSI Set Up>Convergence Main Menu>Solution>FSI Set Up>Convergence По умолчанию величина сходимости равна 0.0001.

Показатель сходимости.

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

Conv = log (x/min)/(log max/min) где:

x = L2 норма каждого параметра, передаваемого через границу твердо жидкого раздела.

min = критерий сходимости, задаваемый командой FSCO max = пороговая величина расхождения (для вычислений принимается равной 100) Показатель сходимости используется для нормирования критерия сходимости для каждого параметра, передаваемого через границу раздела. Сходимость для параметра, передаваемого через границу твердо жидкого раздела определяется по следующим условиям:

Когда x min, Conv 0 и удовлетворен критерий сходимости.

Когда x max, Conv 1 и решение расходится.

3.4.3.5. Задания частоты сохранения результатов расчета.

Вы имеете возможность сохранять выходные данные каждый временной шаг или каждый N ый временной шаг. Выходные данные с заданной частотой записываются в файл результатов анализа жидкости (Jobname.RFL) и в файл результатов анализа твердой части модели (Jobname.RST или Jobname.RTH).

Для установки частоты сохранения выходных данных используйте следующее:

Команда:

FSOU GUI:

Main Menu>Preprocessor>FSI Set Up>Convergence Main Menu>Solution>FSI Set Up>Convergence Вы можете использовать команды FLDATA2, FLDATA4, FLDATA4A, FLDATA5, или OUTRES для установки частоты записи выходных данных как для анализа жидкости, так и для анализа твердой части модели. Команда задания частоты сохранения выходных данных FSOU перезаписывает частоту сохранения, заданную другими командами.

3.4.3.6. Установка величин релаксации.

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

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

Для задания значений релаксации используйте следующее:

Команда:

FSRE GUI:

Main Menu>Preprocessor>FSI Set Up>Relaxation Main Menu>Solution>FSI Set Up>Relaxation 3.4.3.7. Составление списка опций решения.

Для проверки состояния опций анализа твердо жидкого взаимодействия, используйте следующее:

Команда:

FSST GUI:

Main Menu>Preprocessor>FSI Set Up>Status Main Menu>Solution>FSI Set Up>Status 3.4.4. Получение решения.

Для запуска решения используйте следующую команду:

Команда:

SOLVE GUI:

Main Menu>Solution>Current LS Можно повторно запустить анализ твердо жидкого взаимодействия с последнего временного шага:

Команда:

FSRS GUI:

Main Menu>Solution>FSI Set Up>Restart Примечание Анализ твердо жидкого взаимодействия не поддерживает мультифреймовый рестарт доступный в структурном анализе.

3.4.5. Постобработка результатов.

Для постобработки результатов анализа твердо жидкого взаимодействия нужно возобновить базу данных. Можно просмотреть результаты, используя стандартные команды ANSYS POST1 и POST26. Не поддерживается одновременная постобработка результатов анализа жидкости и анализа твердой части модели. Результаты анализа жидкости сохраняются в файле Jobname.RFL и результаты анализа твердой части модели сохраняются в файле Jobname.RST или Jobname.RTH. Информацию по постобработке смотрите в An Overview of Postprocessing в ANSYS Basic Analysis Guide.

3.5. На заметку.

Существует несколько моментов о которых не следует забывать:

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

• Оболочечные и балочные элементы не могут иметь элементы жидкости с обеих сторон.

• Решатель DDS не поддерживает анализ твердо жидкого взаимодействия.

• Для задач, требующих небольшую сжимаемость в жидкой области, вы можете задать модуль объемного сжатия жидкости, используя команду FLDATA16,BULK,BETA.

• Для ослабления колебаний структурной области задайте конструкционное демпфирование для твердой области модели, используя команды ALPHAD and BETAD.

3.6. Простой двумерный анализ твердо жидкого взаимодействия (пакетный или командный метод).

3.6.1. Постановка задачи.

В данном примере рассматривается твердо жидкое взаимодействие в простом поршне с неоднородной конечно элементной сеткой на границе раздела. Движение поршня описывается синусоидальным законом. Рисунок демонстрирует к/э сетку и граничные условия. Необходимо определить распределение скорости и смещений на стороне жидкости.

Простой поршень.

3.6.2. Исходные данные.

Размеры поршня и цилиндра:

Длина цилиндра = 1. Диаметр цилиндра = 0. Высота поршня = 0. Свойства твердого материала:

Модуль Юнга = 1.0e Плотность = Коэффициент Пуассона = 0. Свойства жидкости:

Модуль объемного сжатия = Плотность = 1. Вязкость = 0. Теплопроводность = 1. Теплоемкость = 1. Движение поршня синусоидальное с максимальным смещением равным 0.85.

3.6.3. Результаты.

Смещение жидкости Скорость жидкости 3.6.4. Листинг программы.

/BATCH MAXDISP = 0. *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C *DEL,_FNC_C *SET,_FNCNAME,'uxdsp' *DIM,_FNC_C1,, *DIM,_FNC_C2,, *SET,_FNC_C1(1),MAXDISP *SET,_FNC_C2(1),2. ! /INPUT,time.func *DIM,%_FNCNAME%,TABLE,6,4, ! Начало уравнения: a*sin(omega*{TIME}) %_FNCNAME%(0,0,1)= 0.0, - %_FNCNAME%(2,0,1)= 0. %_FNCNAME%(3,0,1)= %_FNC_C1(1)% %_FNCNAME%(4,0,1)= %_FNC_C2(1)% %_FNCNAME%(5,0,1)= 0. %_FNCNAME%(6,0,1)= 0. %_FNCNAME%(0,1,1)= 1.0, -1, 0, 1, 18, 3, %_FNCNAME%(0,2,1)= 0.0, -1, 9, 1, -1, 0, %_FNCNAME%(0,3,1)= 0, -2, 0, 1, 17, 3, - %_FNCNAME%(0,4,1)= 0.0, 99, 0, 1, -2, 0, ! Конец уравнения: a*sin(omega*{TIME}) !--> /PREP ET,1, KEYOPT,1,4, ET,2, RECT,0.1,1.0,,0. TYPE, LESIZE,4,,, LESIZE,2,,, LESIZE,1,,, LESIZE,3,,, ALLSEL AMESH,ALL ALLSEL,BELOW,AREA NSEL,R,LOC,X,0. SF,ALL,FSIN, ALLSEL,BELOW,AREA NSEL,R,LOC,X, D,ALL,UX, D, ALL, UY, D, ALL, PRES, ALLSEL, BELOW, AREA NSEL,R,LOC,Y, D, ALL, UY, D, ALL, VY, ALLSEL, BELOW, AREA NSEL,R,LOC,Y,0. D,ALL,UY, D,ALL,VY, ALLSEL RECT,,0.1,,0. LESIZE,5,,, LESIZE,7,,, LESIZE,8,,, LESIZE,6,,, TYPE, MAT, ALLSEL ASEL,S,,, ALLSEL,BELOW,AREA AMESH,ALL ALLSEL,BELOW,AREA NSEL,R,LOC,X,0. SF,ALL,FSIN, ALLSEL,BELOW,AREA NSEL,R,EXT D,ALL,UX,%uxdsp% D,ALL,UY, ALLSEL ! Свойства твердого материала MP,EX,1,2.0E MP,DENS,1, MP,PRXY,1,0. FLDATA1,SOLU,FLOW, FLDATA1,SOLU,TEMP, FLDATA1,SOLU,TRAN, FLDATA1,SOLU,ALE, FLDATA2,TIME,GLOB,10, ! Свойства жидкости FLDATA7,PROT,DENS,AIR-SI FLDATA7,PROT,VISC,AIR-SI FINISH /SOLU FSAN,ON ! включение FSI анализа FSCO,ALL,1.0e-2 ! сходимость для всех параметров границы раздела FSRE,ALL,1.0 ! релаксация для всех параметров границы раздела FSOU,1 ! частота сохранения выходных данных DELTIM,0.0125 ! приращение времени для SOLID анализа FLDA, TIME, STEP, 0.0125 ! приращение времени для FLUID анализа FSIT, 10 ! число итераций для каждого временного шага FSTI, 1.0 ! время окончания сопряженного анализа FSDT,0.0125 ! приращение времени для сопряженного анализа FSOR,SOLID ! порядок решения сопряж-го анализа (SOLID первый) FSTR,SOLID, ! нестационарный анализ для SOLID FSTR,FLUID, ! нестационарный анализ для FLUID SAVE SOLVE SAVE FINISH /SOLU FSRS,TIME,1.0 ! время рестарта для сопряженного анализа FSTI,1.25 ! время окончания сопряженного анализа FSDT,0.0125 ! приращение времени для сопряженного анализа SOLVE SAVE FINISH /POST FILE,,rfl NSOL,2,23,UX NSOL,3,23,VX PRVAR, PLVAR, PRVAR, PLVAR, FINISH 3.7. Простой трехмерный анализ твердо жидкого взаимодействия (пакетный или командный метод) 3.7.1. Постановка задачи В этом примере рассматривается анализ взаимодействия течения, омывающего твердый объект. Известна скорость на входе для первого решения. При рестарте входная скорость задается равной нулю. Необходимо определить распределение скорости и смещение жидкости.

3.7.2. Исходные данные.

Свойства твердого материала:

Модуль Юнга = 2.0e9 N/m Плотность = 7800 Kg/m Коэффициент Пуассона = 0. Свойства жидкости:

Pages:     | 1 |   ...   | 4 | 5 ||



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.